« zpět na výpis měsíce |

[Talk-cz] RUIAN info

Vlákno 27.1. - 4.2.2014, počet zpráv: 21


27.1.2014 08:40:12 (#1)
gravatar

Marián Kyral

<mkyral at email.cz>
2462 2837
Ahoj, v souvislosti s RUIAN mne napadla ještě jedna věc. Líbil by se mi nějaký nástroj na principu traceru, který by ale místo tracování tvaru objektu vrátil informace o daném bodě z RUIAN. Tedy, kliknu v JOSM do mapy a otevře se mi okno, ve kterém budou infomace z RUIAN. U pozemku způsob využití, případně vlastník, u budov, typ budovy, vlastník a info o pozemku na kterém ta budova stojí. Nevím přesně, jaké informace jsou v RUIAN dostupné. Mohl by to být samostatný zásuvný modul, případně bych to mohl přidat jako další režim do Tracer zásuvného modulu. Jediné, co potřebuji, je nějaký server, odkud bych tahal data. Formát předpokládám json, ale klidně i jiný rozumný. Co myslíte? Marián

27.1.2014 09:34:02 (#2)
gravatar

Pavel Kwiecien

<pavel.kwiecien at seznam.cz>
48 5000
Ahoj, U stavebních objektů to jsou například následující položky (možná nejsou všechny): CastObce, CisloDomov, GlobalniId, Id, IdTransakc, IsknBudova, Kod, PlatiOd, PocetBytu, PripojeniK, PripojeniP, PripojeniV, RizeniId, TypOchrany, TypStavebn, VybaveniVy, ZpusobVyta, ZpusobVyuz, ZpusobyO_1, ZpusobyOch. U parcel to jsou následující položky (možná nejsou všechny): KmenoveCis, VymeraParc, ZpusobyVyu, DruhCislov, DruhPozemk, Kod, PlatiOd, IdTransakc, RizeniId, Pododdelen. Kromě způsobu využití u budov a pozemků to jsou z pohledu OSM asi dost nezajímavé položky. Zdraví Pavel Kwiecien
---------- Původní zpráva ---------- Od: Marián Kyral <mkyral na email.cz> Datum: 27. 1. 2014 Předmět: [Talk-cz] RUIAN info "Ahoj, v souvislosti s RUIAN mne napadla ještě jedna věc. Líbil by se mi nějaký nástroj na principu traceru, který by ale místo tracování tvaru objektu vrátil informace o daném bodě z RUIAN. Tedy, kliknu v JOSM do mapy a otevře se mi okno, ve kterém budou infomace z RUIAN. U pozemku způsob využití, případně vlastník, u budov, typ budovy, vlastník a info o pozemku na kterém ta budova stojí. Nevím přesně, jaké informace jsou v RUIAN dostupné. Mohl by to být samostatný zásuvný modul, případně bych to mohl přidat jako další režim do Tracer zásuvného modulu. Jediné, co potřebuji, je nějaký server, odkud bych tahal data. Formát předpokládám json, ale klidně i jiný rozumný. Co myslíte? Marián _______________________________________________ Talk-cz mailing list Talk-cz na openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz" ------------- další část --------------- HTML příloha byla odstraněna... URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20140127/2129a78e/attachment.html>

28.1.2014 07:22:46 (#3)
gravatar

JV

<j.v.2 at seznam.cz>
57
Dobrý den, v nestrukturovaném tvaru to lze zobrazit přímo z veřejného dálkového přístupu: http://vdp.cuzk.cz/vdp/ruian/stavebniobjekty/id_stavebniho_objektu http://vdp.cuzk.cz/vdp/ruian/parcely/id_parcely a podobně. Nezobrazuje se tam přímo vlastník (RUIAN neobsahuje údaje o vlastnících), jsou tam ale odkazy do Nahlížení do KN. J. Veselý
---------- Původní zpráva ---------- Od: Marián Kyral <mkyral na email.cz> Datum: 27. 1. 2014 Předmět: [Talk-cz] RUIAN info "Ahoj, v souvislosti s RUIAN mne napadla ještě jedna věc. Líbil by se mi nějaký nástroj na principu traceru, který by ale místo tracování tvaru objektu vrátil informace o daném bodě z RUIAN. Tedy, kliknu v JOSM do mapy a otevře se mi okno, ve kterém budou infomace z RUIAN. U pozemku způsob využití, případně vlastník, u budov, typ budovy, vlastník a info o pozemku na kterém ta budova stojí. Nevím přesně, jaké informace jsou v RUIAN dostupné. Mohl by to být samostatný zásuvný modul, případně bych to mohl přidat jako další režim do Tracer zásuvného modulu. Jediné, co potřebuji, je nějaký server, odkud bych tahal data. Formát předpokládám json, ale klidně i jiný rozumný. Co myslíte? Marián _______________________________________________ Talk-cz mailing list Talk-cz na openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz" ------------- další část --------------- HTML příloha byla odstraněna... URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20140128/8e505cbc/attachment.html>

28.1.2014 10:28:53 (#4)
gravatar

hanoj

<ehanoj at gmail.com>
718
zobrazit citaci
> > Co myslíte?
*** moc tomu nerozumim, ale WMS CUZK KM má standadní GetFeatureInfo např. http://openlayers.org/dev/examples/getfeatureinfo-control.html hanoj ------------- další část --------------- HTML příloha byla odstraněna... URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20140128/4cac76b4/attachment.html>

28.1.2014 12:23:39 (#5)
gravatar

Marián Kyral

<mkyral at email.cz>
2462 2837
Dne 28.1.2014 10:28, hanoj napsal: zobrazit citaci
>>> 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

28.1.2014 12:35:59 (#6)
gravatar

JV

<j.v.2 at seznam.cz>
57
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: zobrazit citaci
>>> 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" ------------- další část --------------- HTML příloha byla odstraněna... URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20140128/443474f5/attachment.html>

28.1.2014 02:08:31 (#7)
gravatar

Marián Kyral

<mkyral at email.cz>
2462 2837
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: zobrazit citaci
> 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]
Links: ------ [1] https://lists.openstreetmap.org/listinfo/talk-cz ------------- další část --------------- HTML příloha byla odstraněna... URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20140128/023a62fc/attachment.html>

29.1.2014 12:00:52 (#8)
gravatar

Marián Kyral

<mkyral at email.cz>
2462 2837
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: zobrazit citaci
> 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]
_______________________________________________ Talk-cz mailing list Talk-cz na openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz [1] Links: ------ [1] https://lists.openstreetmap.org/listinfo/talk-cz ------------- další část --------------- HTML příloha byla odstraněna... URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20140129/4a85ec1c/attachment.html>

31.1.2014 12:00:00 (#9)
gravatar

Marián Kyral

<mkyral at email.cz>
2462 2837
Ahoj, Tak nakonec bych to viděl takhle: http://www.kyralovi.cz/tmp/josm/pointinfo_building.png http://www.kyralovi.cz/tmp/josm/pointinfo_garden.png http://www.kyralovi.cz/tmp/josm/pointinfo_swamp.png 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: zobrazit citaci
> 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]
_______________________________________________ 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] Links: ------ [1] https://lists.openstreetmap.org/listinfo/talk-cz ------------- další část --------------- HTML příloha byla odstraněna... URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20140131/f7c2b5fc/attachment.html>

31.1.2014 07:37:33 (#10)
gravatar

Dalibor Jelínek

<dalibor at dalibor.cz>
415 1552
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 http://www.kyralovi.cz/tmp/josm/pointinfo_garden.png http://www.kyralovi.cz/tmp/josm/pointinfo_swamp.png 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":"zastav\u011bn\u00e1 plocha a n\u00e1dvo\u0159\u00ed", "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 <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> &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 <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> &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 <mailto:mkyral na email.cz> > Datum: 28. 1. 2014 Předmět: Re: [Talk-cz] RUIAN info Dne 28.1.2014 10:28, hanoj napsal: zobrazit citaci
>>> 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 <mailto:Talk-cz na openstreetmap.org> https://lists.openstreetmap.org/listinfo/talk-cz _______________________________________________ Talk-cz mailing list Talk-cz na openstreetmap.org <mailto:Talk-cz na openstreetmap.org> https://lists.openstreetmap.org/listinfo/talk-cz _______________________________________________ Talk-cz mailing list Talk-cz na openstreetmap.org <mailto:Talk-cz na openstreetmap.org> https://lists.openstreetmap.org/listinfo/talk-cz _______________________________________________ Talk-cz mailing list Talk-cz na openstreetmap.org <mailto:Talk-cz na openstreetmap.org> https://lists.openstreetmap.org/listinfo/talk-cz ------------- další část --------------- HTML příloha byla odstraněna... URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20140131/e01bd921/attachment.html>

31.1.2014 11:45:50 (#11)
gravatar

Marián Kyral

<mkyral at email.cz>
2462 2837
Jo jo. Přesně k tomu jsem to dělal. On je sice typ pozemku značen přímo na Katastrální mapě, ale já si ty značky nepamatuji :-( Během víkendu bych to chtěl dotáhnout do konce, hodit na github a zeptat se na možnost začlenění v josm-dev konferenci. Ale uvidím jak mi to vyjde, rodinka mi už rozplánovala celý víkend a možná se k tomu nedostanu. Marián Dne 31.1.2014 07:37, Dalibor Jelínek napsal: zobrazit citaci
> 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]
_______________________________________________ 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] _______________________________________________ Talk-cz mailing list Talk-cz na openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz [1] Links: ------ [1] https://lists.openstreetmap.org/listinfo/talk-cz [2] http://www.kyralovi.cz/tmp/josm/pointinfo_building.png [3] http://www.kyralovi.cz/tmp/josm/pointinfo_garden.png [4] http://www.kyralovi.cz/tmp/josm/pointinfo_swamp.png [5] http://services.cuzk.cz/wms.asp?REQUEST=GetFeatureInfo&amp;SERVICE=WMS&amp;VERSION=1.1.1&amp;LAYERS=DEF_PARCELY&amp;SRS=EPSG:102067&amp;BBOX=-823894,-1066830,-823041,-1066262&amp;WIDTH=1279&amp;HEIGHT=853&amp;QUERY_LAYERS=DEF_PARCELY&amp;X=577&amp;Y=350 [6] http://services.cuzk.cz/wms.asp?REQUEST=GetFeatureInfo&amp;SERVICE=WMS&amp;VERSION=1.1.1&amp;LAYERS=DEF_BUDOVY&amp;SRS=EPSG:102067&amp;BBOX=-823662,-1066899,-822809,-1066330&amp;WIDTH=1279&amp;HEIGHT=853&amp;QUERY_LAYERS=DEF_BUDOVY&amp;X=722&amp;Y=381 [7] http://openlayers.org/dev/examples/getfeatureinfo-control.html ------------- další část --------------- HTML příloha byla odstraněna... URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20140131/081c663a/attachment.html>

2.2.2014 02:32:19 (#12)
gravatar

Marián Kyral

<mkyral at email.cz>
2462 2837
Ahoj, dokončil jsem první, testovací verzi. http://www.kyralovi.cz/tmp/josm/pointinfo_beta.png Ke stažení zde: http://www.kyralovi.cz/tmp/josm/pointinfo.jar Zdrojové kódy: https://github.com/mkyral/josm-pointInfo Umí to v podstatě vše, co jsem chtěl. Zatím se mi nedaří přinutit okno "Connecting server" aby po skončení zmizelo. Takže pořád vespod straší. Taky jsem zatím neřešil kompletní překlad do češtiny - stejný problém jako u Tracer Pluginu. Je možné, že bude problém s čísly orientačními v případě, že jich je na jedné budově více. Prosím o otestování a nahlášení problémových míst. 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 Plugin by měl být lehce rozšiřitelný o další moduly. Jako vstup jsou LatLon souřadnice, jako výstup se očekává html text. Co dělá modul uvnitř, odkud stahuje a jak je zpracuje je jeho věc. Jediné co bude potřeba, je dopsat nějakou konfigurační stránku, kde si budu uživatel vybrat, o které moduly má zájem. 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. Omlouvám se, pokud není geometrické názvosloví správně, už je to pár let, co jsem ze školy. Snad je to pochopitelné. Bohužel vůbec netuším, jak tohle udělat s daty v RUIAN databázi. Uvítám jakýkoli tip jak na to. Díky, Marián

2.2.2014 05:31:45 (#13)
gravatar

Petr Vejsada

<osm at propsychology.cz>
516
Ahoj, Dne Ne 2. února 2014 14:32:19, Marián Kyral napsal(a): zobrazit citaci
> 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
:-), i v Praze jsou díry zobrazit citaci
> 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; zobrazit citaci
> 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

2.2.2014 06:30:04 (#14)
gravatar

Marián Kyral

<mkyral at email.cz>
2462 2837
Dne 2.2.2014 17:31, Petr Vejsada napsal: zobrazit citaci
> 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). >
Z RUIAN. A za ty tečny se omlouvám, myslel jsem úsečky, které spojují protilehlé strany parcely a které procházejí bodem X. V podstatě by výsledkem téměř vždy měla být jen jedna ulice a to ta, která protíná úsečku nejblíže bodu X. zobrazit citaci
> 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
Tomu ještě rozumím ;-) zobrazit citaci
> <-> > <muj_bod_kam_jsem:kliknul typu geometry point a stejného SRID, jako má > deiniční čára> limit 1;
Huh, cože? zobrazit citaci
>> 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.
To by šlo. zobrazit citaci
> > -- > Petr > > > _______________________________________________ > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz

2.2.2014 06:45:21 (#15)
gravatar

Petr Vejsada

<osm at propsychology.cz>
516
Dne Ne 2. února 2014 18:30:04, Marián Kyral napsal(a): zobrazit citaci
> > <-> > > <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

2.2.2014 07:28:12 (#16)
gravatar

Marián Kyral

<mkyral at email.cz>
2462 2837
Dne 2.2.2014 18:45, Petr Vejsada napsal: zobrazit citaci
> 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/ >
OK. Tak to by bylo k tomu záhadnému operátoru "<->". A teď ještě "typu geometry point a stejného SRID, jako má deiniční čára" Typ geometry point je 'POINT( lon lat)' ? A jak zjistím, jaké SRID má ta čára? Koukal jsem na definici na SRID na wiki, ale moc jsem to nepobral. Ještě upozorňuji, že PostGis vůbec, ale vůbec a opravdu ani trochu neznám. Nějaký tip na jednoduchý dokument vysvětlující základy? Co jsem zatím našel moc užitečné nebylo :-( Marián zobrazit citaci
> -- > Petr > > _______________________________________________ > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz

2.2.2014 08:23:29 (#17)
gravatar

Petr Vejsada

<osm at propsychology.cz>
516
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

2.2.2014 10:49:23 (#18)
gravatar

Marián Kyral

<mkyral at email.cz>
2462 2837
Moc díky. Myslím, že pro začátek mi to stačí. Marián Dne 2.2.2014 20:23, Petr Vejsada napsal(a): zobrazit citaci
> 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 >

4.2.2014 07:27:37 (#19)
gravatar

Marián Kyral

<mkyral at email.cz>
2462 2837
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),4326),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 ; 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. Díky, Marián Dne 2.2.2014 20:23, Petr Vejsada napsal: zobrazit citaci
> 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

4.2.2014 09:26:45 (#20)
gravatar

Petr Vejsada

<osm at propsychology.cz>
516
Ahoj, Dne Út 4. února 2014 19:27:37, Marián Kyral napsal(a): zobrazit citaci
> 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í. zobrazit citaci
> 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

4.2.2014 09:56:22 (#21)
gravatar

Marián Kyral

<mkyral at email.cz>
2462 2837
Super, teď už to funguje správně. Taky mě ten order by mohl napadnout :-D Marián Dne 4.2.2014 21:26, Petr Vejsada napsal: zobrazit citaci
> 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

« zpět na výpis měsíce