> > Co myslíte?

>>> Co myslíte? *** moc tomu nerozumim, ale WMS CUZK KM má standadní >>> GetFeatureInfo např. >>> http://openlayers.org/dev/examples/getfeatureinfo-control.html [2] > > hanoj >

> Dobrý den, > http://services.cuzk.cz/wms.asp?REQUEST=GetFeatureInfo&SERVICE=WMS&VERSION=1.1.1&LAYERS=DEF_PARCELY&SRS=EPSG:102067&BBOX=-823894,-1066830,-823041,-1066262&WIDTH=1279&HEIGHT=853&QUERY_LAYERS=DEF_PARCELY&X=577&Y=350 > http://services.cuzk.cz/wms.asp?REQUEST=GetFeatureInfo&SERVICE=WMS&VERSION=1.1.1&LAYERS=DEF_BUDOVY&SRS=EPSG:102067&BBOX=-823662,-1066899,-822809,-1066330&WIDTH=1279&HEIGHT=853&QUERY_LAYERS=DEF_BUDOVY&X=722&Y=381 > > J. Veselý > > ---------- Původní zpráva ---------- > Od: Marián Kyral <mkyral na email.cz> > Datum: 28. 1. 2014 > Předmět: Re: [Talk-cz] RUIAN info > >> Dne 28.1.2014 10:28, hanoj napsal: >> >>>>> Co myslíte? *** moc tomu nerozumim, ale WMS CUZK KM má standadní >>>>> GetFeatureInfo např. >>>>> http://openlayers.org/dev/examples/getfeatureinfo-control.html [2] >>> >>> hanoj >>> >> >> Díky za tip. To vypadá slibně. Bohužel o WMS toho moc nevím. takže si to >> budu muset nastudovat. Nějaký tip na přehlednou a snadno pochopitelnou >> dokumentaci? A nejlépe i nějaký funkční příklad url pro CUZK, abych si s >> tím mohl hrát. >> >> Díky, >> Marián >> >> _______________________________________________ >> Talk-cz mailing list >> Talk-cz na openstreetmap.org >> https://lists.openstreetmap.org/listinfo/talk-cz > > _______________________________________________ > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz [1]

> Díky, > > večer na to kouknu. Zatím to ale vypadá, že tam těch informací moc není. U parcely v podstatě pouze odkaz do VDP, u budovy je ještě navíc číslo popisné. To nejdůležitější pro mapování, způsob využití, je až ve VDP. Takže jeden dotaz pro zjištění ruian id a druhý do VDP pro zjištění způsobu využití. Navíc už to je html a to bude potřeba parsovat.No čekal jsem víc :-( > > Marián > > Dne 28.1.2014 12:35, JV napsal: Dobrý den, > http://services.cuzk.cz/wms.asp?REQUEST=GetFeatureInfo&SERVICE=WMS&VERSION=1.1.1&LAYERS=DEF_PARCELY&SRS=EPSG:102067&BBOX=-823894,-1066830,-823041,-1066262&WIDTH=1279&HEIGHT=853&QUERY_LAYERS=DEF_PARCELY&X=577&Y=350 > http://services.cuzk.cz/wms.asp?REQUEST=GetFeatureInfo&SERVICE=WMS&VERSION=1.1.1&LAYERS=DEF_BUDOVY&SRS=EPSG:102067&BBOX=-823662,-1066899,-822809,-1066330&WIDTH=1279&HEIGHT=853&QUERY_LAYERS=DEF_BUDOVY&X=722&Y=381 > > J. Veselý > > ---------- Původní zpráva ---------- > Od: Marián Kyral <mkyral na email.cz> > Datum: 28. 1. 2014 > Předmět: Re: [Talk-cz] RUIAN info > Dne 28.1.2014 10:28, hanoj napsal: > >>>> Co myslíte? *** moc tomu nerozumim, ale WMS CUZK KM má standadní >>>> GetFeatureInfo např. >>>> http://openlayers.org/dev/examples/getfeatureinfo-control.html [2] >> >> hanoj >> > > Díky za tip. To vypadá slibně. Bohužel o WMS toho moc nevím. takže si to > budu muset nastudovat. Nějaký tip na přehlednou a snadno pochopitelnou > dokumentaci? A nejlépe i nějaký funkční příklad url pro CUZK, abych si s > tím mohl hrát. > > Díky, > Marián > > _______________________________________________ > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz > > _______________________________________________ > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz [1]

> Petr Vejsada mi laskavě zřídil přístup na jeho server, takže jsem se mohl porozhlídnout v databázi RUIAN. Trochu jsem si pohrál a momentálně dokáži vygenerovat následující údaje ve formátu json: > > { > "coordinates":{ > "lat":"49.66918", > "lon":"18.3651988" > }, > "source":"cuzk:ruian", > "stavebni_objekt":{ > "ruian_id":"20381174", > "cislo_domovni":"{46}", > "pocet_podlazi":"1", > "zpusob_vyuziti":"stavba pro administrativu", > "plati_od":"22.12.2013" > }, > "parcela":{ > "ruian_id":"3410187802", > "druh_pozemku":"zastavu011bnu00e1 plocha a nu00e1dvou0159u00ed", > "zpusob_vyuziti":null, > "plati_od":"01.06.2007" > } > } > > Předpokládám, že tam ještě přidám minimálně adresu. Taktéž nebude problém přidat odkaz do VDP, případně i jinam. Jako parametr jsou pouze souřadnice LatLon. To je pro mne mnohem jednodušší, než u WMS zjišťovat pixely v bboxu. I když, možná to JOSM má někde jednoduše přístupné. V každém případě, získané údaje z WMS neodpovídají očekávaní. > > Zkusím spíchnout první nástřel pluginu pro JOSM, který zobrazí ten JSON. Pak hodím nějaký screenshot a zkušební jar. > > Marián > > Dne 28.1.2014 14:08, Marián Kyral napsal: > > Díky, > > večer na to kouknu. Zatím to ale vypadá, že tam těch informací moc není. U parcely v podstatě pouze odkaz do VDP, u budovy je ještě navíc číslo popisné. To nejdůležitější pro mapování, způsob využití, je až ve VDP. Takže jeden dotaz pro zjištění ruian id a druhý do VDP pro zjištění způsobu využití. Navíc už to je html a to bude potřeba parsovat.No čekal jsem víc :-( > > Marián > > Dne 28.1.2014 12:35, JV napsal: Dobrý den, > http://services.cuzk.cz/wms.asp?REQUEST=GetFeatureInfo&SERVICE=WMS&VERSION=1.1.1&LAYERS=DEF_PARCELY&SRS=EPSG:102067&BBOX=-823894,-1066830,-823041,-1066262&WIDTH=1279&HEIGHT=853&QUERY_LAYERS=DEF_PARCELY&X=577&Y=350 > http://services.cuzk.cz/wms.asp?REQUEST=GetFeatureInfo&SERVICE=WMS&VERSION=1.1.1&LAYERS=DEF_BUDOVY&SRS=EPSG:102067&BBOX=-823662,-1066899,-822809,-1066330&WIDTH=1279&HEIGHT=853&QUERY_LAYERS=DEF_BUDOVY&X=722&Y=381 > > J. Veselý > > ---------- Původní zpráva ---------- > Od: Marián Kyral <mkyral na email.cz> > Datum: 28. 1. 2014 > Předmět: Re: [Talk-cz] RUIAN info > Dne 28.1.2014 10:28, hanoj napsal: > >>>> Co myslíte? *** moc tomu nerozumim, ale WMS CUZK KM má standadní >>>> GetFeatureInfo např. >>>> http://openlayers.org/dev/examples/getfeatureinfo-control.html [2] >> >> hanoj >> > > Díky za tip. To vypadá slibně. Bohužel o WMS toho moc nevím. takže si to > budu muset nastudovat. Nějaký tip na přehlednou a snadno pochopitelnou > dokumentaci? A nejlépe i nějaký funkční příklad url pro CUZK, abych si s > tím mohl hrát. > > Díky, > Marián > > _______________________________________________ > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz > > _______________________________________________ > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz [1]

> Cau, > > to je hezky, to uz bych rad mel v JOSM. :-) > > Bude se mi to hodit, kdyz mapuju plochy, pak je informace > > z katastru o typu pozemku celkem uzitecna. > > Zdravi, > > Dalibor > > FROM: Marián Kyral [mailto:mkyral na email.cz] > SENT: Friday, January 31, 2014 12:00 AM > TO: OpenStreetMap Czech Republic > SUBJECT: Re: [Talk-cz] RUIAN info > > Ahoj, > > Tak nakonec bych to viděl takhle: > > http://www.kyralovi.cz/tmp/josm/pointinfo_building.png [2] > > http://www.kyralovi.cz/tmp/josm/pointinfo_garden.png [3] > > http://www.kyralovi.cz/tmp/josm/pointinfo_swamp.png [4] > > Ještě adresa a myslím, že to bude. Mělo by to být lehce rozšířitelné i na jiné zdroje - stačí udělat nový modul a nějakou konfigurační stránku. Jen mne teď jiný zdroj nenapadá. Maximálně tak wikipedie - třeba zobrazení informací o obci - jestli by to k něčemu bylo. > > Marián > > Dne 29.1.2014 00:00, Marián Kyral napsal: > > Petr Vejsada mi laskavě zřídil přístup na jeho server, takže jsem se mohl porozhlídnout v databázi RUIAN. Trochu jsem si pohrál a momentálně dokáži vygenerovat následující údaje ve formátu json: > > { > "coordinates":{ > "lat":"49.66918", > "lon":"18.3651988" > }, > "source":"cuzk:ruian", > "stavebni_objekt":{ > "ruian_id":"20381174", > "cislo_domovni":"{46}", > "pocet_podlazi":"1", > "zpusob_vyuziti":"stavba pro administrativu", > "plati_od":"22.12.2013" > }, > "parcela":{ > "ruian_id":"3410187802", > "druh_pozemku":"zastavu011bnu00e1 plocha a nu00e1dvou0159u00ed", > "zpusob_vyuziti":null, > "plati_od":"01.06.2007" > } > } > > Předpokládám, že tam ještě přidám minimálně adresu. Taktéž nebude problém přidat odkaz do VDP, případně i jinam. Jako parametr jsou pouze souřadnice LatLon. To je pro mne mnohem jednodušší, než u WMS zjišťovat pixely v bboxu. I když, možná to JOSM má někde jednoduše přístupné. V každém případě, získané údaje z WMS neodpovídají očekávaní. > > Zkusím spíchnout první nástřel pluginu pro JOSM, který zobrazí ten JSON. Pak hodím nějaký screenshot a zkušební jar. > > Marián > > Dne 28.1.2014 14:08, Marián Kyral napsal: > > Díky, > > večer na to kouknu. Zatím to ale vypadá, že tam těch informací moc není. U parcely v podstatě pouze odkaz do VDP, u budovy je ještě navíc číslo popisné. To nejdůležitější pro mapování, způsob využití, je až ve VDP. Takže jeden dotaz pro zjištění ruian id a druhý do VDP pro zjištění způsobu využití. Navíc už to je html a to bude potřeba parsovat.No čekal jsem víc :-( > > Marián > > Dne 28.1.2014 12:35, JV napsal: > > Dobrý den, > http://services.cuzk.cz/wms.asp?REQUEST=GetFeatureInfo&SERVICE=WMS&VERSION=1.1.1&LAYERS=DEF_PARCELY&SRS=EPSG:102067&BBOX=-823894,-1066830,-823041,-1066262&WIDTH=1279&HEIGHT=853&QUERY_LAYERS=DEF_PARCELY&X=577&Y=350 [5] > http://services.cuzk.cz/wms.asp?REQUEST=GetFeatureInfo&SERVICE=WMS&VERSION=1.1.1&LAYERS=DEF_BUDOVY&SRS=EPSG:102067&BBOX=-823662,-1066899,-822809,-1066330&WIDTH=1279&HEIGHT=853&QUERY_LAYERS=DEF_BUDOVY&X=722&Y=381 [6] > > J. Veselý > > ---------- Původní zpráva ---------- > Od: Marián Kyral <mkyral na email.cz> > Datum: 28. 1. 2014 > Předmět: Re: [Talk-cz] RUIAN info > > Dne 28.1.2014 10:28, hanoj napsal: > >>>> Co myslíte? *** moc tomu nerozumim, ale WMS CUZK KM má standadní >>>> GetFeatureInfo např. >>>> http://openlayers.org/dev/examples/getfeatureinfo-control.html [7] [2] >> >> hanoj >> > > Díky za tip. To vypadá slibně. Bohužel o WMS toho moc nevím. takže si to > budu muset nastudovat. Nějaký tip na přehlednou a snadno pochopitelnou > dokumentaci? A nejlépe i nějaký funkční příklad url pro CUZK, abych si s > tím mohl hrát. > > Díky, > Marián > > _______________________________________________ > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz [1] > > _______________________________________________ > > Talk-cz mailing list > > Talk-cz na openstreetmap.org > > https://lists.openstreetmap.org/listinfo/talk-cz [1]

> Ještě upozorním, že RUIAN data nejsou dostupná všude. Byl jsem nemile > překvapen, když jsem zjistil, že chybí kus Ostravy. To jsem nečekal :-D

> Ještě by se mi líbila jedna věc: Pokud kliknu na plochu, kde leží nějaká > komunikace, tak by se mohl zobrazit název této komunikace (pokud má - > třeba ulice). > > Teoreticky tuším, jak by se to dalo řešit: > 1) Pokud je parcela typu, který mne zajímá - silnice, ostatní > komunikace... tak: > 2) Najdu tečny s hranicí parcely na kterých daný bod leží. > 3) Vezmu kratší tečnu a hledám ulici, která tuto tečnu protíná. > 4) Pokud takováto ulice existuje, vrátím její název.

> Může se stát, že těch ulic bude více - křižovatky "Y". Pak bych asi > nevrátil nic.

> Ahoj, > > Dne Ne 2. února 2014 14:32:19, Marián Kyral napsal(a): > >> Ještě by se mi líbila jedna věc: Pokud kliknu na plochu, kde leží >> nějaká >> komunikace, tak by se mohl zobrazit název této komunikace (pokud má - >> třeba ulice). >> >> Teoreticky tuším, jak by se to dalo řešit: >> 1) Pokud je parcela typu, který mne zajímá - silnice, ostatní >> komunikace... tak: >> 2) Najdu tečny s hranicí parcely na kterých daný bod leží. >> 3) Vezmu kratší tečnu a hledám ulici, která tuto tečnu protíná. >> 4) Pokud takováto ulice existuje, vrátím její název. > > Odkud chceš brát název ulice? Z OSM nebo RUIAN? No to je jedno, není > třeba > dělat tečny (tečna je přímka a tedy neexistuje kratší/delší tečna). >

> Není problém najít nejbližší ulici (komunikaci) k bodu a třeba i > takovou, > která leží uvnitř plochy. > > Takže klik na bod, najdeme parcelu a nejbližší komunikaci, která > zároveň leží > (st_intersects) alespoň částí uvnitř parcely. Počet takových komunikací > se > bude limitně blížit jedné ;). > > Jestli ti jde jen o název ulice, tak stačí najít nejbližší ulici. > > Nejrychlejší hledání je asi > > select nazev from rn_ulice order by definicni_cara

> <-> > <muj_bod_kam_jsem:kliknul typu geometry point a stejného SRID, jako má > deiniční čára> limit 1;

>> Může se stát, že těch ulic bude více - křižovatky "Y". Pak bych asi >> nevrátil nic. > > Skoro vždy je něco nejbližší, s výjimkou, kdy je stejně daleko více než > jedna > ulice. Spíš nevracet nic, když to "nejbližší" bude třeba kilometr > daleko.

> > -- > Petr > > > _______________________________________________ > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz

> > <-> > > <muj_bod_kam_jsem:kliknul typu geometry point a stejného SRID, jako má > > deiniční čára> limit 1; > > Huh, cože?

> Dne Ne 2. února 2014 18:30:04, Marián Kyral napsal(a): > >> > <-> >> > <muj_bod_kam_jsem:kliknul typu geometry point a stejného SRID, jako má >> > deiniční čára> limit 1; >> >> Huh, cože? > > To je finta Postgisu. Normálně nemůže existovat index, ve kterém bys > hledal > podle vzdálenosti od určitého bodu či jiné geometrie, protože ten bod > není v > době vytváření indexu znám. Tohle to umí i tak > > http://boundlessgeo.com/2011/09/indexed-nearest-neighbour-search-in-postgis/ >

> -- > Petr > > _______________________________________________ > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz

> Ahoj, > > Dne Ne 2. února 2014 19:28:12, Marián Kyral napsal(a): > > definiční čára, mně blbne na klávesnici F, fakt, nedělám si srandu... > > Geometry typ prostě obsahuje údaje o tvaru nějakého objektu. Zobrazení: > > třeba > select st_astext(definicni_cara) > > Geometrie může být bod, čára, polygon, multipolygon, oblouky a též směs všeho > uvedeného. > > SRID = identifikace projekce, ve kreré je geometrie vyjádřena. Používáme v OSM > dvě - 900913 pro zobrazování a 4326, což jsou ty stupně, jak je známe z > běžného života. > > RUIAN mám uložen v 900913, aby s tím Mapnik neměl tolik práce. > > select st_srid(definicni_cara) from ruian.rn_ulice limit 1; > > st_srid > --------- > 900913 > (1 řádka) > > Když chceš najít nejbližší ulici k bodu, vyjádřenému v "občanských" > souřadnicích, tedy 4326, je potřeba to převést do 900913, protože takto jsou > uložena data (platí jen pro tu databázi, co mám na serveru; někdo jiný to mlže > mít jinak). > > To se dělá transformací - select > st_transform(<geometrie_v_lidských_souřadnicích>,<do jakého systému>) > > Bod jakožto data typu geometrie vytvoříš třeba funkcí > st_makepoint(<lon>,<lat>). Tato funkce vrátí data typu geometry a je to point. > > Není však jasné, v jakém souřadnicovém systému to vlastně je. Proto použijeme > st_setsrid(<geometry>,<srid>) > > Takže do PG zadáme bod třeba > > select st_setsrid(st_makepoint(14,50),4326) > > > Lidsky čitelné totéž: > > select st_astext(st_setsrid(st_makepoint(14,50),4326)); > st_astext > -------------- > POINT(14 50) > (1 řádka) > > > V jakém souřadnicovém systému to vlastně máme? > > select st_srid(st_setsrid(st_makepoint(14,50),4326)); > st_srid > --------- > 4326 > (1 řádka) > > > Ten samý bod si zobrazíme v lidsky čitelné formě v souřadnicovém systému > 900913: > > select st_astext(st_transform(st_setsrid(st_makepoint(14,50),4326),900913)); > st_astext > ------------------------------------------ > POINT(1558472.87110583 6446275.84101716) > (1 řádka) > > A konečně se dostáváme k cíli: > > select nazev from ruian.rn_ulice order by definicni_cara <-> > st_transform(st_setsrid(st_makepoint(14,50),4326),900913) limit 1; > nazev > ---------- > Na Zámku > (1 řádka) > > > Jak je ta ulice daleko? > > Funkce st_distance(geometry,geometry) > > Ta ráda vrací vzdálenost v radiánech na zemském povrchu, osobně dávám přednost > metrům ;-). V metrech nám to řekne, když geometrie budou geografie. Geography > je podobný datový typ, vyjadřuje se v souřadnicovém systému 4326. > > Toho se docílí přetypováním na typ geography. > > Celý select pak vypadá: > > select nazev,st_distance( (st_transform(definicni_cara,4326))::geography, > (st_setsrid(st_makepoint(14,50),4326))::geography ) from ruian.rn_ulice order > by definicni_cara <-> st_transform(st_setsrid(st_makepoint(14,50),4326),900913) > limit 1; > > nazev | st_distance > ----------+-------------- > Na Zámku | 26.405619556 > (1 řádka) > > Nejbližší bod ulice Na Zámku je od nás vzdálen 26 metrů a 40 centimetrů. > > Cvičení: > > Jak si zobrazíme lidsky čitelné souřadnice adresního bodu z RUIAN? > > A: > select st_astext(st_transform(definicni_bod,4326)) from ruian.rn_adresni_misto > where kod=21411409; > > st_astext > ------------------------------------------ > POINT(13.6752996130858 49.2886006596294) > (1 řádka) > > > Učebnice Postgisu: > > http://workshops.boundlessgeo.com/postgis-intro/ > > je IMO super, ale nedá se to za 10 minut > > -- > Petr > > > _______________________________________________ > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz >

> Ahoj, > > Dne Ne 2. února 2014 19:28:12, Marián Kyral napsal(a): > > definiční čára, mně blbne na klávesnici F, fakt, nedělám si srandu... > > Geometry typ prostě obsahuje údaje o tvaru nějakého objektu. Zobrazení: > > třeba > select st_astext(definicni_cara) > > Geometrie může být bod, čára, polygon, multipolygon, oblouky a též směs > všeho > uvedeného. > > SRID = identifikace projekce, ve kreré je geometrie vyjádřena. > Používáme v OSM > dvě - 900913 pro zobrazování a 4326, což jsou ty stupně, jak je známe z > běžného života. > > RUIAN mám uložen v 900913, aby s tím Mapnik neměl tolik práce. > > select st_srid(definicni_cara) from ruian.rn_ulice limit 1; > > st_srid > --------- > 900913 > (1 řádka) > > Když chceš najít nejbližší ulici k bodu, vyjádřenému v "občanských" > souřadnicích, tedy 4326, je potřeba to převést do 900913, protože takto > jsou > uložena data (platí jen pro tu databázi, co mám na serveru; někdo jiný > to mlže > mít jinak). > > To se dělá transformací - select > st_transform(<geometrie_v_lidských_souřadnicích>,<do jakého systému>) > > Bod jakožto data typu geometrie vytvoříš třeba funkcí > st_makepoint(<lon>,<lat>). Tato funkce vrátí data typu geometry a je to > point. > > Není však jasné, v jakém souřadnicovém systému to vlastně je. Proto > použijeme > st_setsrid(<geometry>,<srid>) > > Takže do PG zadáme bod třeba > > select st_setsrid(st_makepoint(14,50),4326) > > > Lidsky čitelné totéž: > > select st_astext(st_setsrid(st_makepoint(14,50),4326)); > st_astext > -------------- > POINT(14 50) > (1 řádka) > > > V jakém souřadnicovém systému to vlastně máme? > > select st_srid(st_setsrid(st_makepoint(14,50),4326)); > st_srid > --------- > 4326 > (1 řádka) > > > Ten samý bod si zobrazíme v lidsky čitelné formě v souřadnicovém > systému > 900913: > > select > st_astext(st_transform(st_setsrid(st_makepoint(14,50),4326),900913)); > st_astext > ------------------------------------------ > POINT(1558472.87110583 6446275.84101716) > (1 řádka) > > A konečně se dostáváme k cíli: > > select nazev from ruian.rn_ulice order by definicni_cara <-> > st_transform(st_setsrid(st_makepoint(14,50),4326),900913) limit 1; > nazev > ---------- > Na Zámku > (1 řádka) > > > Jak je ta ulice daleko? > > Funkce st_distance(geometry,geometry) > > Ta ráda vrací vzdálenost v radiánech na zemském povrchu, osobně dávám > přednost > metrům ;-). V metrech nám to řekne, když geometrie budou geografie. > Geography > je podobný datový typ, vyjadřuje se v souřadnicovém systému 4326. > > Toho se docílí přetypováním na typ geography. > > Celý select pak vypadá: > > select nazev,st_distance( > (st_transform(definicni_cara,4326))::geography, > (st_setsrid(st_makepoint(14,50),4326))::geography ) from ruian.rn_ulice > order > by definicni_cara <-> > st_transform(st_setsrid(st_makepoint(14,50),4326),900913) > limit 1; > > nazev | st_distance > ----------+-------------- > Na Zámku | 26.405619556 > (1 řádka) > > Nejbližší bod ulice Na Zámku je od nás vzdálen 26 metrů a 40 > centimetrů. > > Cvičení: > > Jak si zobrazíme lidsky čitelné souřadnice adresního bodu z RUIAN? > > A: > select st_astext(st_transform(definicni_bod,4326)) from > ruian.rn_adresni_misto > where kod=21411409; > > st_astext > ------------------------------------------ > POINT(13.6752996130858 49.2886006596294) > (1 řádka) > > > Učebnice Postgisu: > > http://workshops.boundlessgeo.com/postgis-intro/ > > je IMO super, ale nedá se to za 10 minut > > -- > Petr > > > _______________________________________________ > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz

> Ahoj, > tak jsem se dopracoval k tomuto selectu: > > select u.kod, u.nazev, ST_asText(st_transform(u.definicni_cara,4326)) > from ( select kod, nazev, definicni_cara > from ruian.rn_ulice > order by definicni_cara <-> > > st_transform(st_setsrid(st_makepoint(18.36564928953012,49.670527512403766),4 > 326),900913) LIMIT 100) as u > where st_distance( (st_transform(u.definicni_cara,4326))::geography, > (st_setsrid(st_makepoint(18.3657227215035,49.66980665513853),4326))::geograp > hy ) < 10 > ; > > Jediná věc mne zarazila. Musel jsem u subselectu nastavit "LIMIT 100". S > nastavením "LIMIT 1" mi to vrátilo výsledek jen někdy. Klikal jsme podél > ulice, na jednom místě mi to ulici vypsalo a o kousek dál už zase ne.

> Teď to vypadá, že to funguje slušně. Jen na křižovatkách to někdy vrátí > vedlejší ulici. Ale to bych neviděl jako moc velký problém. Inteligentní > uživatel klikne kousek dál od křižovatky.

> Ahoj, > > Dne Út 4. února 2014 19:27:37, Marián Kyral napsal(a): > >> Ahoj, >> tak jsem se dopracoval k tomuto selectu: >> >> select u.kod, u.nazev, ST_asText(st_transform(u.definicni_cara,4326)) >> from ( select kod, nazev, definicni_cara >> from ruian.rn_ulice >> order by definicni_cara <-> >> >> st_transform(st_setsrid(st_makepoint(18.36564928953012,49.670527512403766),4 >> 326),900913) LIMIT 100) as u >> where st_distance( (st_transform(u.definicni_cara,4326))::geography, >> (st_setsrid(st_makepoint(18.3657227215035,49.66980665513853),4326))::geograp >> hy ) < 10 >> ; >> >> Jediná věc mne zarazila. Musel jsem u subselectu nastavit "LIMIT 100". >> S >> nastavením "LIMIT 1" mi to vrátilo výsledek jen někdy. Klikal jsme >> podél >> ulice, na jednom místě mi to ulici vypsalo a o kousek dál už zase ne. > > Já tě tím nechtěl pro začátek zatěžovat, ale teď vidím, že jsem měl, > no, > alespoň sis na to přišel sám. > > V tom odkazu, co jsem posílal, je psáno, že třídění pomocí <-> je > přibližné, > protože geometrie (čára) se převede na bbox. Hledáš-li vzdálenost bodu > od > bodu, pak bbox bodu je totéž jako bod samotný, ale bbox čáry je prostě > obdelník, takový, aby se do něho ta čára vešla. > > Tak, jak jsi to vyřešil, to řeším i já, tedy podobně. Na konec toho > druhého > selectu dávám order by st_distance(...) limit 1. st_distance už třídí > exaktně, > sice pomalu, ale u 100 položek to zase nevadí. U 20M položek (parcely) > je to > sakra znát. > > Takže ten tvůj select by v mém podání vypadal: > > ........ 326),900913) LIMIT 100) as u > where st_distance( (st_transform(u.definicni_cara,4326))::geography, > (st_setsrid(st_makepoint(18.3657227215035,49.66980665513853),4326))::geography > ) < 10 order by st_distance( > (st_transform(u.definicni_cara,4326))::geography, > (st_setsrid(st_makepoint(18.3657227215035,49.66980665513853),4326))::geography > limit 1 > > to vrátí s vysokou pravděpodobností tu správnou ulici, ledaže by ani > těch 100 > řádku, nalezených aproximovaným výběrem z bboxů, neobsahovalo ten řádek > se > skutečně nejbližší ulicí. > >> Teď to vypadá, že to funguje slušně. Jen na křižovatkách to někdy >> vrátí >> vedlejší ulici. Ale to bych neviděl jako moc velký problém. >> Inteligentní >> uživatel klikne kousek dál od křižovatky. > > protože tam nemáš ten order_by st_distance. Vrátí ti to prostě první > řádek, > který vyhovuje podmínce, že st_distance < 10 a to může být klidně i ta > druhá > ulice. > > -- > Petr > > > _______________________________________________ > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz