Dne 26.12.2016 v 23:55 Petr Morávek [Xificurk] napsal(a):
> 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.
>
Snažím se to napasovat na KM. Konkrétně na tuhle definici:
<map>
<tag key='name' value='český CUZK: barevný'/>
<tag key='type' value='wms'/>
<tag key='url'
value='http://services.cuzk.cz/wms/wms.asp?FORMAT=image/png&VERSION=1.1.1&SERVICE=WMS&REQUEST=GetMap&LAYERS=DEF_PARCELY,DEF_BUDOVY,RST_PK_I,RST_KMD_I,dalsi_p_mapy_i,obrazy_parcel_i,parcelni_cisla_i,hranice_parcel_barevne,omp_i&STYLES=&SRS={proj}&WIDTH={width}&HEIGHT={height}&BBOX={bbox}&TRANSPARENT=true'/>
<tag key='pixel-per-eastnorth' value='9.045115508227324'/>
<tag key='max-zoom' value='100'/>
<tag key='projections' value=''/>
</map>
Žádnou projekci tam nemám, předpokládám, že se použije nějaká výchozí.
Marián
> S pozdravem,
> Petr Morávek aka Xificurk
>
> _______________________________________________
> Talk-cz mailing list
> Talk-cz@openstreetmap.org
> https://lists.openstreetmap.org/listinfo/talk-cz
_______________________________________________
Talk-cz mailing list
Talk-cz@openstreetmap.org
https://lists.openstreetmap.org/listinfo/talk-cz
