RE: [FreeGeoCZ] Konverze z S-JTSK pomoci gdalwarp

[Štěpán Kafka <kafka AT email.cz>]

NO, obávám se, že něco posunutého mají na CENII, taky jsme s tím měli potíže

RNDr. Štěpán Kafka
Help Service Remote Sensing
http://www.bnhelp.cz
e-mail: kafka AT email.cz
tel.: 327 514 118
mob.: 724 254 731

> -----Original Message-----
> From: freegeocz-bounces AT fsv.cvut.cz [mailto:freegeocz-
> bounces AT fsv.cvut.cz] On Behalf Of Michal
> Sent: Thursday, April 02, 2009 6:13 PM
> To: freegeocz AT fsv.cvut.cz
> Subject: [FreeGeoCZ] Konverze z S-JTSK pomoci gdalwarp
> 
> 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