Re: [FreeGeoCZ] transformace mezi souradnymi ref. sys.

[Michal <michal AT lightcomp.cz>]

hanoj wrote:
> Nicmene jsem narazil na problemy s GDAL, resp. ogr2ogr, gdalwarp, kde
> nejsem schopen [4],[6] pouzit (mimo jeden pripad[5]) parametr s
> transf. klicem +towgs84, tak aby transformace probehla s danou
> presnosti (<1m).
>
> Muzou mi zde zkusenejsi napomoci k uspesnemu pouzivani GDAL v bash?

Ja jsem nedavno resil stejny problem pri prevodu S-JTSK -> google. U me 
pomohlo prevadet pres longlat. T.j.

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=longlat +datum=WGS84" \
  s-jtsk.tif tmp_longlat.tif

gdalwarp \
  -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" \
  tmp_longlat.tif google.tif


Michal