[FreeGeoCZ] Konverze z S-JTSK pomoci gdalwarp

[Michal <michal AT lightcomp.cz>]

Zdravim.

Postahoval jsem si mapy CR 1:10000 (puvodem z geoportal.cenia.cz) z 
rapidshare:
http://rapidshare.com/users/PXQC3K
http://rapidshare.com/users/WQUS78
http://rapidshare.com/users/NOEMBM

Chtel jsem si je prekonvertovat na TMS dlazdice, abych je mohl pouzit 
treba v maemo-mapper nebo tangogps, a nejak jsem narazil. Soubory maji v 
nazvu S-JTSK souradnice, takze map soubor vubec nepotrebuji. Na 
geoportal.cenia.cz jsem si overil, ze jsou souradnice v poradku. 
Napriklad soubor 13-22_1043000-1034000_602500-593500.png jsem 
konvertoval nasledovne:

# 256 png -> truecolor tiff
pngtopnm 13-22_1043000-1034000_602500-593500.png | pnmtotiff -truecolor 
> tmp.tif

gdal_translate \
  -a_srs '+proj=krovak +lat_0=49.5 +lon_0=24.83333333333333 +alpha=0 
+k=0.9999 +x_0=0 +y_0=0 +ellps=bessel 
+towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56 +units=m +no_defs' \
  -gcp 0 0 -602500 -1034000 \
  -gcp 0 10799 -602500 -1043000 \
  -gcp 10799 0 -593500 -1034000 \
  -gcp 10799 10799 -593500 -1043000 \
  tmp.tif tmp_ref.tif

gdalwarp \
  -s_srs '+proj=krovak +lat_0=49.5 +lon_0=24.83333333333333 +alpha=0 
+k=0.9999 +x_0=0 +y_0=0 +ellps=bessel 
+towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56 +units=m +no_defs' \
  -t_srs '+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 
+x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs' \
  -r cubic \
  -tps \
  tmp_ref.tif 13-22_1043000-1034000_602500-593500.tif

gdal2tiles.py 13-22_1043000-1034000_602500-593500.tif




Vysledek je mozne videt na:

http://michal.etc.cz/ZM10/13-22_1043000-1034000_602500-593500
http://michal.etc.cz/ZM10/05-44_975000-966000_789500-780500
http://michal.etc.cz/ZM10/23-51_1128000-1119000_849000-840000
http://michal.etc.cz/ZM10/34-20_1221500-1212500_585500-576500

Vse je trochu posunute a zda se mi, ze kazda mapa o trochu jinou 
vzdalenost. Nevim jestli delam neco spatne, nebo je chyba v gdalwarp, 
ale v qlandkartegt ma vysledny geotiff divne souradnice. Levy horni roh 
souboru 13-22_1043000-1034000_602500-593500.png je:

S-JTSK        -602500 -1034000
qlandkartegt  50.324466667 16.357533333
cs2cs         50.323727056 16.356153468

Ten rozdil +- odpovida tomu posunu. Tusi nekdo, kde by mohl byt problem? 
Pokud by si s tim chtel nekdo pohrat, tak puvodni png soubory je mozne 
stahnout z:

http://michal.etc.cz/ZM10/13-22_1043000-1034000_602500-593500.png
http://michal.etc.cz/ZM10/05-44_975000-966000_789500-780500.png
http://michal.etc.cz/ZM10/23-51_1128000-1119000_849000-840000.png
http://michal.etc.cz/ZM10/34-20_1221500-1212500_585500-576500.png

Michal