Čistě matematická transformace dosahuje značné nepřesnosti v závislosti na lokalitě až 70 metrů. http://grass.fsv.cvut.cz/gwiki/Chyba_při_transformaci_z_WGS84_do_S-JTSK
Převod GML se dělá pomocí ogr2ogr, případně lze GDAL knihovnu začlenit do Céčkového programu nebo do Python skriptu bez problémů pomocí API: http://gdal.org/gdal_tutorial.html MK ----- Original Message ----- From: Pavel Machek [mailto:pa...@ucw.cz] To: OpenStreetMap Czech Republic [mailto:talk-cz@openstreetmap.org] Sent: Sat, 30 Jun 2012 12:45:23 +0200 Subject: Re: [Talk-cz] Data RUIAN - výměnný formát > On Sat 2012-06-23 04:45:21, Jan Bilak wrote: > Díky. Nevíte o nějakých > knihovnách (.NET, Java, JavaScript, C, C++, > ...) pro transformaci pomocí > toho S-JTSK gridu? Nebo, pokud nejsou > přímo knihovny, ve kterých > opensource programech s vhodnou licencí by > tato transformace šla > najít? Ja tu mam: gdalwarp -s_srs "+proj=krovak +a=6377397.155 > +rf=299.1528128 +no_defs +towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56" > -t_srs "+proj=latlong +a=6378137 +rf=298.257223563 > +no_defs +towgs84=0.000,0.000,0.000" > /data/gis/READ-ONLY/cechy.tif /tmp/delme.tiff a pak pascalovej zdrojak: { > Copyright 2005 Zdenek Hrdina, distribute under GPLv2 } procedure > jtsk_wgs( X,Y,Hel:double; var B,L,H:double); {Vypocet zemepisnych souradnic > v systemu WGS-84 z rovinnych souradnic S-JTSK a elipsoidicke > vysky} procedure transformace_BLH(var B,L,H: double); {Transformace > zemepisnych souradnic z JTSK do WGS} var > lat,lon,alt,x1,y1,z1,x2,y2,z2:double; ... Poslu nebo by mel jit > vygooglit. -- (english) http://www.livejournal.com/~pavelmachek (cesky, > pictures) > http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html _______________________________________________ Talk-cz > mailing > list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz _______________________________________________ Talk-cz mailing list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz
