Dne 24.12.2016 v 00:04 Ha Noj napsal(a): >> 1) Dělám něco špatně já? > *** GRID je spočítán mezi ETRS89 a JTSK, takže by tam mělo být ještě > něco mezi ETRS89 a WGS84. Ale nepotěším tě, neboť to dělá jen cca 30 cm. > > echo "18.845223 49.547998" | cs2cs +init=epsg:4326 +to +init=epsg:4258 > +towgs84="0.0473,0.0467,-0.0253,0.000891,0.00539,-0.008772,-0.00158" -f > "%.6f" | cs2cs +init=epsg:4258 +to +init=epsg:5514 > +nadgrids=./jezek08_jtskcz.llb -f "%.3f" > -432758.771 -1136759.330 -0.009
Uff, musím říct, že tvoje zmínka o ETRS89 (a rozdílnosti oproti WGS84) pro mne otevřela nová dvířka, která si nejsem jistý, že jsem chtěl otevřít :-)) Zajímalo by mne odkud pochází tvoje transformační parametry towgs84? Vůbec se mi je nepodařilo nikde vygooglit a obecně jsem ve spoustě zdrojů narážel na doporučení považovat oba systémy za "shodné". A když se podívám na definici v /usr/share/proj/epsg tak tam je: > # ETRS89 > <4258> +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs <> Snažil jsem se dál empiricky dopátrat, jaké transformace provádí ČÚZK. Jako dataset jsem použil adresní místa - pro každou obec (>6000) jsem náhodně vybral jedno adresní místo (to by snad mělo zajistit reprezentativní vzorek). Z WFS jsem pak stáhl souřadnice v EPSG:4258 a EPSG:4326, na ně vyzkoušel různé transformace a výsledek porovnával s původními souřadnicemi v EPSG:5514 z WFS. A mám pocit, že čím víc se v tom rýpu, tím větší zmatek v tom mám :/ 1) EPSG:4258 -> EPSG:5514 pomocí gridu, např. echo "16.914199 49.148716" | cs2cs +init=epsg:4258 +to +init=epsg:5514 +nadgrids=czech -f "%.3f" -576504.892 -1168238.275 -0.000 Výsledná odchylka na souboru adresních bodů: 4 +/- 2 cm (max. 15 cm) HURÁ! Zdá se, že grid skutečně funguje skvěle a souhlasí s tím, co posílá ČÚZK. 2) EPSG:4258 -> EPSG:5514 pomocí sedmiprvkové transformace, např. echo "16.914199 49.148716" | cs2cs +init=epsg:4258 +to +init=epsg:5514 +towgs84="570.8,85.7,462.8,4.998,1.587,5.261,3.56" -f "%.3f" -576504.701 -1168238.186 -44.382 Výsledná odchylka na souboru adresních bodů: 16 +/- 11 cm (max. 82 cm) Podle očekávání dává sedmiprvková transformace horší výsledky. 3) EPSG:4326 -> EPSG:5514 pomocí gridu, např. echo "16.914193 49.148715" | cs2cs +init=epsg:4326 +to +init=epsg:5514 +nadgrids=czech -f "%.3f" -576505.338 -1168238.341 0.000 Výsledná odchylka na souboru adresních bodů: 18 +/- 16 cm (max. 116 cm) Tohle je to, co používám teď na import administrativních hranic (AFAIK to samé je na poloha.net) a řeším, proč se ta čísla v některých případech rozcházejí o více jak metr. 4) EPSG:4326 -> EPSG:4258 -> EPSG:5514 pomocí gridu, např. echo "16.914193 49.148715" | cs2cs +init=epsg:4326 +to +init=epsg:4258 +towgs84="0.0473,0.0467,-0.0253,0.000891,0.00539,-0.008772,-0.00158" -f "%.8f" | cs2cs +init=epsg:4258 +to +init=epsg:5514 +nadgrids=czech -f "%.3f" -576505.114 -1168238.150 -0.009 Výsledná odchylka na souboru adresních bodů: 31 +/- 10 cm (max. 86 cm) Nápad s přidáním transformace WGS84<->ETRS89 sice stáhnul maximální chybu o těch 30 cm, na druhou stranu průměrná chyba šla nahoru, takže tohle asi taky nebude správná cesta (?) 5) EPSG:4326 -> EPSG:5514 pomocí sedmiprvkové transformace, např. echo "16.914193 49.148715" | cs2cs +init=epsg:4326 +to +init=epsg:5514 +towgs84="570.8,85.7,462.8,4.998,1.587,5.261,3.56" -f "%.3f" -576505.148 -1168238.251 -44.382 Výsledná odchylka na souboru adresních bodů: 14 +/- 7 cm (max. 36 cm) Tohle je další překvapení - tahle cesta nejlépe reprodukuje transformaci ČÚZK přestože víme, že rozhodně není nejpřesnější. 6) EPSG:4326 -> EPSG:4258 -> EPSG:5514 pomocí sedmiprvkové transformace, např. echo "16.914193 49.148715" | cs2cs +init=epsg:4326 +to +init=epsg:4258 +towgs84="0.0473,0.0467,-0.0253,0.000891,0.00539,-0.008772,-0.00158" -f "%.8f" | cs2cs +init=epsg:4258 +to +init=epsg:5514 +towgs84="570.8,85.7,462.8,4.998,1.587,5.261,3.56" -f "%.3f" -576504.923 -1168238.061 -44.391 Výsledná odchylka na souboru adresních bodů: 24 +/- 9 cm (max. 49 cm) Tohle už uvádím jen pro úplnost. Abych pravdu řekl, tak vůbec netuším, jaký z toho všeho udělat závěr. Vypadá to, že grid, co jsme používali doteď, stále "funguje" správně. A naopak jsou to data ČÚZK v EPSG:4326, kterým se nedá moc věřit. @Marián: Ty se snažíš vyřešit zarovnání dat, co posílá Petr z poloha.net (hádám, že vznikají transformací ze zdrojového EPSG:5514 do latlon pomocí gridu) s čím konkrétně? S WMS podklady? A ty teď taháš v jaké projekci? EPSG:4326? Jestli jo, tak to zkus vyměnit za EPSG:4258. S pozdravem, Petr Morávek aka Xificurk _______________________________________________ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz