Import adres a vzdálenost AM od SO na Talk-cz

Petr Vejsada osm na propsychology.cz #m589dd6 24. 5. 2014 16:26

Ahoj, přestalo mě bavit dohledávat domy, ke kterým patří adresa, když jsou kilometry daleko. Tak jsem upravil algoritmus, aby v případě, že je AM od SO dále než 50 metrů, použil definiční bod SO namísto souřadnic AM. Algoritmus moc nefungoval a tak jsem se pustil do pátrání, jak je to možné. Inu - není žádného nesmyslu, který by mohl existovat jen teoreticky. Moje fantazie není tak bohatá, abych dokázal vymyslet nesmysl, který by přede mnou ČÚZK již neuvedl do praxe. Takže tu máme definiční bod stavebního objektu, který leží 221 km daleko od jeho hranic (slovy dvě stě dvacet jeden kilometr). pedro=> select kod,distance_meters(definicni_bod,hranice) as distance from ruian.rn_stavebni_objekt where hranice is not NULL and definicni_bod is not NULL order by distance desc limit 10; kod | distance ----------+------------------ 51833077 | 221482.323142855 13549308 | 8630.303164072 41460022 | 6718.926763008 14596784 | 6509.639499045 22863940 | 6323.752735814 27885631 | 6323.517911398 22863966 | 6316.406526901 1606387 | 6313.425573178 22863796 | 6304.353506391 27885712 | 6292.054487993 (10 řádek) pedro=> select count(*) from ruian.rn_stavebni_objekt where distance_meters(hranice,definicni_bod) > 1000 and hranice is not NULL and definicni_bod is not NULL; count ------- 459 (1 řádka) Praha nebo Brno, komu na tom záleží, co. Už mě to fakt s***.

hanoj ehanoj na gmail.com #md067aa 24. 5. 2014 21:07

kod | distance ----------+------------------ 51833077 | 221482.323142855

*** jedna se o tyto dva nize? Ja tam totiz problem nevidim... http://vdp.cuzk.cz/vdp/ruian/stavebniobjekty/51833077 http://vdp.cuzk.cz/vdp/ruian/adresnimista/41487770 ha hanoj

Petr Vejsada osm na propsychology.cz #m42fa08 24. 5. 2014 21:20

Ahoj, ano, to je ono. Adresní místo v tomto případě má geometrii a ta je shodná s definičním bodem. Až na tu "maličkost", že hranice, tedy vlastní polygon budovy, se nenachází v Klatovech, ale 221km odsud, v Poličce. pedro=# select local_textgeom(definicni_bod) from ruian.rn_stavebni_objekt where kod=51833077; local_textgeom ----------------------------------- 49.3928060531079 13.2536505761335 (1 řádka) pedro=# select local_textgeom(st_centroid(hranice)) from ruian.rn_stavebni_objekt where kod=5183307> local_textgeom ----------------------------------- 49.7136924259116 16.2748977303676 (1 řádka) zeměpisná délka se liší o "zanedbatelné" 3 stupně...

hanoj ehanoj na gmail.com #mada175 24. 5. 2014 21:56

ano, to je ono. Adresní místo v tomto případě má geometrii a ta je shodná s definičním bodem. Až na tu "maličkost", že hranice, tedy vlastní polygon budovy, se nenachází v Klatovech, ale 221km odsud, v Poličce.

*** IMHO so 51833077 polygon budovy nema: http://pastebin.com/AvUQxjvf http://vdp.cuzk.cz/vymenny_format/soucasna/20140430_OB_555771_UZSZ.xml.gz k.ú. Klatovy, kam "Kal č.p.25" spadá nemá ještě DKM, tudíž ani polygony v RUIAN. http://vdp.cuzk.cz/vdp/ruian/katastralniuzemi/665797 ha hanoj

Petr Vejsada osm na propsychology.cz #m22a549 24. 5. 2014 22:25

Ahoj, UZSZ - vybral sis základní datovou sadu. Ta polygony neobsahuje, je třeba použít sadu UKSH. Bez ohledu na to, zda DKM v Klatovech je či není, tak soubor 20140430_OB_555771_UKSH.xml.gz tam ten polygon prostě má a leží v Poličce. http://pastebin.com/6SeTEnUb

hanoj ehanoj na gmail.com #m77c699 24. 5. 2014 22:52

Uz to vidim (UZSZ neni UKSH), mas recht. Ta geometrie so 51833077 je ukradena tomuto domu so 51474468: http://pastebin.com/G0RXydht http://vdp.cuzk.cz/vdp/ruian/stavebniobjekty/51474468 hezky ;) ha hanoj

Petr Souček soucekp na atlas.cz #mc91381 26. 5. 2014 19:33

Dobrý den, děkuji za objevení tohoto problému. Již to máme zapsané v našem interním helpdesku a problém řešíme. Ozvu se hned, jak zjistíme příčinu a opravíme ji. Vypadá to na nějakou chybu v párování budov z ISKN a stavebních objektů v ISÚI, které se provádí až v RÚIAN. Petr Souček