Posuny - pokus o zpřesnění na Talk-cz

Petr Vejsada osm na propsychology.cz #m73bf2d 4. 2. 2014 18:31

Ahoj, udělal jsem experimentální vrstvu s budovami, která, pokud jsem něco nezvoral, by měla být podle Xificurka s použitím gridu, ale možná jsem fakt něco nedomyslel, pže to dopadlo nic moc. http://pedro.poloha.net/mapa , url vrstvy je http://tile.poloha.net/temp_budovy/z/x/y.png Jak jsem postupoval: - zavedl jsem fiktivní SRID 999 do spatial_ref_sys, kde jsem změnil proj4text tak, jak to má Xificurk - stáhnul jsem si grid podle Petrova odkazu (btw: Petře, pokud to čteš, ta zdrojová forma by nebyla?) - zkopíroval tabulku s geometriemi budov - updatnul geometrii takto: update temp_budovy set hranice= (st_transform(st_setsrid(st_transform(hranice,5514),999),900913)) což by mělo dělat, že se z 900913 geometrie přepočítá zpět na 5514 a podle gridu se zase přepočítá zpátky do 900913. Je to podle mých očí jen horší. Proč?

Marián Kyral mkyral na email.cz #m9681e5 4. 2. 2014 20:21

Ahoj, udělal jsem experimentální vrstvu s budovami, která, pokud jsem něco nezvoral, by měla být podle Xificurka s použitím gridu, ale možná jsem fakt něco nedomyslel, pže to dopadlo nic moc. http://pedro.poloha.net/mapa , url vrstvy je http://tile.poloha.net/temp_budovy/z/x/y.png Jak jsem postupoval: - zavedl jsem fiktivní SRID 999 do spatial_ref_sys, kde jsem změnil proj4text tak, jak to má Xificurk - stáhnul jsem si grid podle Petrova odkazu (btw: Petře, pokud to čteš, ta zdrojová forma by nebyla?) - zkopíroval tabulku s geometriemi budov - updatnul geometrii takto: update temp_budovy set hranice= (st_transform(st_setsrid(st_transform(hranice,5514),999),900913)) což by mělo dělat, že se z 900913 geometrie přepočítá zpět na 5514 a podle gridu se zase přepočítá zpátky do 900913. Je to podle mých očí jen horší. Proč?

Jak na to tak koukám, tak se to posunulo na opačnou stranu. Místo dolů a vlevo je nová vrstva posunuta nahoru a vpravo. Někde tam změň znaménko :D Marián

"Petr Morávek [Xificurk]" petr na pada.cz #m628862 4. 2. 2014 21:16

Ahoj,

Ahoj, udělal jsem experimentální vrstvu s budovami, která, pokud jsem něco nezvoral, by měla být podle Xificurka s použitím gridu, ale možná jsem fakt něco nedomyslel, pže to dopadlo nic moc. http://pedro.poloha.net/mapa , url vrstvy je http://tile.poloha.net/temp_budovy/z/x/y.png Jak jsem postupoval: - zavedl jsem fiktivní SRID 999 do spatial_ref_sys, kde jsem změnil proj4text tak, jak to má Xificurk - stáhnul jsem si grid podle Petrova odkazu (btw: Petře, pokud to čteš, ta zdrojová forma by nebyla?)

bohužel git.zcu.cz se zdá opravdu mrtvý... a já bohužel ten zdrojový soubor nikde na disku nemám, ale podařilo se mi ho dohledat v archivu: http://web.archive.org/web/20091003020944/http://git.zcu.cz/grid/czech.lla Nicméně je to z roku 2009, tak nevím jestli to je stejná verze.

- zkopíroval tabulku s geometriemi budov - updatnul geometrii takto: update temp_budovy set hranice= (st_transform(st_setsrid(st_transform(hranice,5514),999),900913)) což by mělo dělat, že se z 900913 geometrie přepočítá zpět na 5514 a podle gridu se zase přepočítá zpátky do 900913. Je to podle mých očí jen horší. Proč?

Já teda do těch transformací moc nevidím, ale není možné, že se tam kumuluje nějaká chyba tím převodem tam a zpátky? Já si teď u sebe pustím nový čerstvý import RUIANu a pak můžem porovnat výstup na nějakém konkrétním stavebním objektu, což? Zdraví, Petr Morávek aka Xificurk

Petr Vejsada osm na propsychology.cz #m1ab4b1 5. 2. 2014 19:33

Ahoj,

- zavedl jsem fiktivní SRID 999 do spatial_ref_sys, kde jsem změnil proj4text tak, jak to má Xificurk - stáhnul jsem si grid podle Petrova odkazu (btw: Petře, pokud to čteš, ta zdrojová forma by nebyla?)
bohužel git.zcu.cz se zdá opravdu mrtvý... a já bohužel ten zdrojový soubor nikde na disku nemám, ale podařilo se mi ho dohledat v archivu: http://web.archive.org/web/20091003020944/http://git.zcu.cz/grid/czech.lla Nicméně je to z roku 2009, tak nevím jestli to je stejná verze.

Zkompiloval jsem, binárky jsou rozdílné (diff], což nic neznamená. Asi nemá smysl dál zkoumat.

(st_transform(st_setsrid(st_transform(hranice,5514),999),900913)) Je to podle mých očí jen horší. Proč?
Já teda do těch transformací moc nevidím, ale není možné, že se tam kumuluje nějaká chyba tím převodem tam a zpátky?

Asi ne. Vytvořil jsem nové schema a naimportovat Aš, Šluknov, Vyšší Brod, Třinec a Jablunkov a podle oka mi to přijde úplně stejné a úplně stejně špatně, jako v předchozím případě. To je opravdu zvláštní. Definici projekce 5514 mám od tebe, grid mám od tebe a s výsledkem nejsem spokojen, kdežto ty se svými daty ano. Je to horší než bez gridu. A mám stejný pocit, jako Marián, že grid to koriguje v podstatě správně, správným směrem, jen s opačnou orientací.

Já si teď u sebe pustím nový čerstvý import RUIANu a pak můžem porovnat výstup na nějakém konkrétním stavebním objektu, což?

Tož to bychom měli. V jakém souřadnicovém systému máš RUIAN? Já v 900913 kvůli vykreslování dlaždic. Mohu ti někde vystavit nějakou část, asi stačí kod a hranice a nechat si třeba v qgisu ukázat st_diff

"Petr Morávek [Xificurk]" petr na pada.cz #m813507 5. 2. 2014 19:48

Ahoj,

(st_transform(st_setsrid(st_transform(hranice,5514),999),900913)) Je to podle mých očí jen horší. Proč?
Já teda do těch transformací moc nevidím, ale není možné, že se tam kumuluje nějaká chyba tím převodem tam a zpátky?
Asi ne. Vytvořil jsem nové schema a naimportovat Aš, Šluknov, Vyšší Brod, Třinec a Jablunkov a podle oka mi to přijde úplně stejné a úplně stejně špatně, jako v předchozím případě. To je opravdu zvláštní. Definici projekce 5514 mám od tebe, grid mám od tebe a s výsledkem nejsem spokojen, kdežto ty se svými daty ano. Je to horší než bez gridu. A mám stejný pocit, jako Marián, že grid to koriguje v podstatě správně, správným směrem, jen s opačnou orientací.

Ahaaa! Teď mi to teprve došlo - on byl RUIAN zpočátku (a v té době jsem zkoumal přesnost) v jiné projekci. Původně to bylo EPSG:2065 a nyní EPSG:5514, která se pokud vím liší právě v prohození souřadnic a změně znamének. Já blbec ten grid slepě bez kontroly převzal... takže ten grid je teď úplně mimo... chtělo by to najít zdrojová data a přegenerovat pro 5514.

Já si teď u sebe pustím nový čerstvý import RUIANu a pak můžem porovnat výstup na nějakém konkrétním stavebním objektu, což?
Tož to bychom měli. V jakém souřadnicovém systému máš RUIAN? Já v 900913 kvůli vykreslování dlaždic. Mohu ti někde vystavit nějakou část, asi stačí kod a hranice a nechat si třeba v qgisu ukázat st_diff

Já to nechávám v originále a transformuji až on-demand, podle toho na co to používám. (Pokud chci počítat vzdálenosti dvou objektů, tak je lepší to dělat v 5514.) Zdraví, Petr Morávek aka Xificurk

Marián Kyral mkyral na email.cz #mdb4564 5. 2. 2014 19:54

Ahoj,
(st_transform(st_setsrid(st_transform(hranice,5514),999),900913)) Je to podle mých očí jen horší. Proč?
Já teda do těch transformací moc nevidím, ale není možné, že se tam kumuluje nějaká chyba tím převodem tam a zpátky?
Asi ne. Vytvořil jsem nové schema a naimportovat Aš, Šluknov, Vyšší Brod, Třinec a Jablunkov a podle oka mi to přijde úplně stejné a úplně stejně špatně, jako v předchozím případě. To je opravdu zvláštní. Definici projekce 5514 mám od tebe, grid mám od tebe a s výsledkem nejsem spokojen, kdežto ty se svými daty ano. Je to horší než bez gridu. A mám stejný pocit, jako Marián, že grid to koriguje v podstatě správně, správným směrem, jen s opačnou orientací.
Ahaaa! Teď mi to teprve došlo - on byl RUIAN zpočátku (a v té době jsem zkoumal přesnost) v jiné projekci. Původně to bylo EPSG:2065 a nyní EPSG:5514, která se pokud vím liší právě v prohození souřadnic a změně znamének. Já blbec ten grid slepě bez kontroly převzal... takže ten grid je teď úplně mimo... chtělo by to najít zdrojová data a přegenerovat pro 5514.

Pokud se nepodaří to přegenerování, tak by možná šlo zjistit odchylku nové pozice od staré a o tu pak korigovat na opačnou stranu. Záleží jak moc náročné by to bylo. Marián

"Petr Morávek [Xificurk]" petr na pada.cz #mf0d750 5. 2. 2014 20:46

Ahoj,
(st_transform(st_setsrid(st_transform(hranice,5514),999),900913)) Je to podle mých očí jen horší. Proč?
Já teda do těch transformací moc nevidím, ale není možné, že se tam kumuluje nějaká chyba tím převodem tam a zpátky?
Asi ne. Vytvořil jsem nové schema a naimportovat Aš, Šluknov, Vyšší Brod, Třinec a Jablunkov a podle oka mi to přijde úplně stejné a úplně stejně špatně, jako v předchozím případě. To je opravdu zvláštní. Definici projekce 5514 mám od tebe, grid mám od tebe a s výsledkem nejsem spokojen, kdežto ty se svými daty ano. Je to horší než bez gridu. A mám stejný pocit, jako Marián, že grid to koriguje v podstatě správně, správným směrem, jen s opačnou orientací.
Ahaaa! Teď mi to teprve došlo - on byl RUIAN zpočátku (a v té době jsem zkoumal přesnost) v jiné projekci. Původně to bylo EPSG:2065 a nyní EPSG:5514, která se pokud vím liší právě v prohození souřadnic a změně znamének. Já blbec ten grid slepě bez kontroly převzal... takže ten grid je teď úplně mimo... chtělo by to najít zdrojová data a přegenerovat pro 5514.

Hm, tak tím to nebude... obě transformace mi dávaj v absolutních číslech stejné výsledky. A až na tu magii se znaménky stále vychází správně test uvedený na http://grass.fsv.cvut.cz/gwiki/S-JTSK-Grid Takže asi jediným zbývajícím zdrojem nepřesnosti je opravdu: "Tato verze je vytvořena pro variantu ETRS89(ETRF1989) tj. platné do 2.1.2011.[1] Vzhledem k posunu souřadnic od té doby platné ETRS89(ETRF2000) neposkytuje níže popsané řešení deklarovanou přesnost." Petr

Petr Vejsada osm na propsychology.cz #mf03445 5. 2. 2014 21:40

Ahoj,

Ahaaa! Teď mi to teprve došlo - on byl RUIAN zpočátku (a v té době jsem zkoumal přesnost) v jiné projekci. Původně to bylo EPSG:2065 a nyní EPSG:5514, která se pokud vím liší právě v prohození souřadnic a změně znamének. Já blbec ten grid slepě bez kontroly převzal... takže ten grid je teď úplně mimo... chtělo by to najít zdrojová data a přegenerovat pro 5514.

Jo, tak to bude ono. Můžeme zkusit tu tabulku z archivu z roku 2009, dohledat (bývalého) majitele gitu, proj-4.8.0/src/projects.h, přepsat do pl/pgsql to, na co tu odkazoval pan Veselý z ČÚZK - mimochodem - ten program ve Fortranu by asi hodně pomohl; ten není k dispozici?)

Já si teď u sebe pustím nový čerstvý import RUIANu a pak můžem porovnat výstup na nějakém konkrétním stavebním objektu, což?

Tak to už asi nemusíme.

Petr Vejsada osm na propsychology.cz #m1fef6d 5. 2. 2014 21:45

Ahoj,

Pokud se nepodaří to přegenerování, tak by možná šlo zjistit odchylku nové pozice od staré a o tu pak korigovat na opačnou stranu. Záleží jak moc náročné by to bylo.

To právě nejde. Od toho je ta tabulka. Na každém místě je ta odchylka jiná. Nejlepším řešením by bylo zprovoznit to http://www.cuzk.cz/Zememerictvi/Geodeticke-zaklady-na-uzemi-CR/GNSS/Nova-realizace-systemu-ETRS89-v-CR.aspx kde, bohužel, není ten zdroják ve Fortranu. Ten se asi prodává? (pane Veselý)?

Jiří Veselý j.v.2 na seznam.cz #md6002b 5. 2. 2014 22:26

Dobrý večer, to nevím, ale poslal jsem dotaz. Když nic jiného, vzorečky jsou i v tom PDF ;-) J. Veselý

Petr Vejsada osm na propsychology.cz #m6b9ad2 5. 2. 2014 23:00

Dobrý večer,

Dobrý večer, to nevím, ale poslal jsem dotaz. Když nic jiného, vzorečky jsou i v tom PDF ;-)

no, přeci jenom přepis z Fortranu do PL/PGSQL by šel asi o mnoho rychleji. Děkujeme moc za vstřícnost!