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

[Talk-cz] Skript na parcely RUIANu

Vlákno 25.7. - 26.7.2014, počet zpráv: 5


25.7.2014 09:22:23 (#1)
gravatar

Martin Švec - OSM

<osm at maatts.cz>
109
Ahoj, obrázky parcel od Petra Vejsady mi nějak nedaly spát, tak jsem včera v noci ubastlil skript na konverzi parcel do OSM. Skriptík nemá jiné ambice než vyzkoušet si natažení parcel do JOSM a vůbec se seznámit s VFR formátem RUIANu, takže nečekejte zázraky. Mj. nefungují multipolygony, beru jen polygony typu gml:LinearRing apod. Tagování jsem narychlo opsal z http://wiki.openstreetmap.org/wiki/Cs:RUIAN, knihovny na zpracování GML apod. jsem zatím nehledal. Vyžaduje perl, proj, Geo::Proj4 (http://search.cpan.org/dist/Geo-Proj4/), a grid jezek_czech08.llb přejmenovaný na "czech" (vyšťouráno tady v archívu, http://freegis.fsv.cvut.cz/gwiki/S-JTSK_/_Grid). Pár obrázků: http://www.maatts.cz/obrazek/3/lomnice-ruian-parcely-png/ http://www.maatts.cz/obrazek/3/blansko-ruian-parcely-png/ http://www.maatts.cz/obrazek/3/brno-ruian-parcely-png/ Vyzkoušeno na uvedených třech obcích = na čtvrté klidně může zhavarovat. ;-) BACHA - ve velkých městech skript žere mraky paměti, optimalizacema jsem se netrápil. :-)) Třeba Brno (kód obce 582786) má ve výsledku 1.122.529 uzlů a 251.328 parcel, a skriptík na to potřebuje aspoň 12GB RAM a/nebo hodně swapu na SSD. Martin ------------- další část --------------- A non-text attachment was scrubbed... Name: par2osm.pl Type: text/x-perl Size: 10580 bytes Desc: [žádný popis není k dispozici] URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20140725/c664e1c3/attachment.pl>

25.7.2014 09:54:00 (#2)
gravatar

Marián Kyral

<mkyral at email.cz>
2501 2837
Ahoj, vypadá to hezky, ale fakt chceme mít v OSM přesný obraz KM? Sem tam si zemědělci něco přiorají, občas zase kus nechají a ten brzy zaroste nějakým křovím nebo lesem. Někde jsem zase viděl velkou zahradu, o kterou se už majitel nechtěl starat, tak jí kus pronajal sousedovi a ten si o to rozšířil své pole. Jak už tu padlo, v KM je zanesen právní stav, ale v OSM potřebujeme skutečnost. Z tohohle pohledu mi přijde pLPIS lepší. Plus to, že je dostupný pro celou republiku. Marián On 25. července 2014 21:22:23 CEST, "Martin Švec - OSM" <osm na maatts.cz> wrote: zobrazit citaci
>Ahoj, > >obrázky parcel od Petra Vejsady mi nějak nedaly spát, tak jsem včera v >noci ubastlil skript na konverzi parcel do OSM. > >Skriptík nemá jiné ambice než vyzkoušet si natažení parcel do JOSM a >vůbec se seznámit s VFR formátem RUIANu, takže nečekejte zázraky. Mj. >nefungují multipolygony, beru jen polygony typu gml:LinearRing apod. >Tagování jsem narychlo opsal z >http://wiki.openstreetmap.org/wiki/Cs:RUIAN, knihovny na zpracování >GML apod. jsem zatím nehledal. > >Vyžaduje perl, proj, Geo::Proj4 >(http://search.cpan.org/dist/Geo-Proj4/), a grid jezek_czech08.llb >přejmenovaný na "czech" (vyšťouráno tady v archívu, >http://freegis.fsv.cvut.cz/gwiki/S-JTSK_/_Grid). > >Pár obrázků: > >http://www.maatts.cz/obrazek/3/lomnice-ruian-parcely-png/ >http://www.maatts.cz/obrazek/3/blansko-ruian-parcely-png/ >http://www.maatts.cz/obrazek/3/brno-ruian-parcely-png/ > >Vyzkoušeno na uvedených třech obcích = na čtvrté klidně může >zhavarovat. ;-) > >BACHA - ve velkých městech skript žere mraky paměti, optimalizacema >jsem se netrápil. :-)) Třeba Brno (kód obce 582786) má ve výsledku >1.122.529 uzlů a 251.328 parcel, a skriptík na to potřebuje aspoň 12GB >RAM a/nebo hodně swapu na SSD. > >Martin > > > > >------------------------------------------------------------------------ > >_______________________________________________ >Talk-cz mailing list >Talk-cz na openstreetmap.org >https://lists.openstreetmap.org/listinfo/talk-cz
-- Odesláno z mého telefonu s Androidem pomocí pošty K-9 Mail. Omluvte prosím moji stručnost. ------------- další část --------------- HTML příloha byla odstraněna... URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20140725/8c10f562/attachment.html>

25.7.2014 10:04:50 (#3)
gravatar

Martin Švec - OSM

<osm at maatts.cz>
109
No, poté co jsem si to vyzkoušel jsem si prakticky jistý, že 1:1 hromadný převod všech parcel do OSM není průchozí :-) Ty miliony rozdrobených parcelek by nadělaly solidní binec v databázi i na mapě. A jak říkáš, zrovna LPIS u zemědělské půdy líp kopíruje realitu. Spíš je to o tom hledat způsoby, jak ty RUIAN data rozumně využít, když už jsou k dispozici. Martin Dne 25.7.2014 21:54, Marián Kyral napsal(a): zobrazit citaci
> Ahoj, vypadá to hezky, ale fakt chceme mít v OSM přesný obraz KM? > Sem tam si zemědělci něco přiorají, občas zase kus nechají a ten > brzy zaroste nějakým křovím nebo lesem. Někde jsem zase viděl velkou > zahradu, o kterou se už majitel nechtěl starat, tak jí kus pronajal > sousedovi a ten si o to rozšířil své pole. > > Jak už tu padlo, v KM je zanesen právní stav, ale v OSM potřebujeme > skutečnost. Z tohohle pohledu mi přijde pLPIS lepší. Plus to, že je > dostupný pro celou republiku. > > Marián > > On 25. července 2014 21:22:23 CEST, "Martin Švec - OSM" > <osm na maatts.cz> wrote: > > Ahoj, > > obrázky parcel od Petra Vejsady mi nějak nedaly spát, tak jsem včera v > noci ubastlil skript na konverzi parcel do OSM. > > Skriptík nemá jiné ambice než vyzkoušet si natažení parcel do JOSM a > vůbec se seznámit s VFR formátem RUIANu, takže nečekejte zázraky. Mj. > nefungují multipolygony, beru jen polygony typu gml:LinearRing apod. > Tagování jsem narychlo opsal z > http://wiki.openstreetmap.org/wiki/Cs:RUIAN, knihovny na zpracování > GML apod. jsem zatím nehledal. > > Vyžaduje perl, proj, Geo::Proj4 > (http://search.cpan.org/dist/Geo-Proj4/), a grid jezek_czech08.llb > přejmenovaný na "czech" (vyšťouráno tady v archívu, > http://freegis.fsv.cvut.cz/gwiki/S-JTSK_/_Grid). > > Pár > obrázků: > > http://www.maatts.cz/obrazek/3/lomnice-ruian-parcely-png/ > http://www.maatts.cz/obrazek/3/blansko-ruian-parcely-png/ > http://www.maatts.cz/obrazek/3/brno-ruian-parcely-png/ > > Vyzkoušeno na uvedených třech obcích = na čtvrté klidně může > zhavarovat. ;-) > > BACHA - ve velkých městech skript žere mraky paměti, optimalizacema > jsem se netrápil. :-)) Třeba Brno (kód obce 582786) má ve výsledku > 1.122.529 uzlů a 251.328 parcel, a skriptík na to potřebuje aspoň 12GB > RAM a/nebo hodně swapu na SSD. > > Martin > > > ---------------------------------------------------------------------- > > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz > > > -- > Odesláno z mého telefonu s Androidem pomocí pošty K-9 Mail. Omluvte > prosím moji stručnost. > > > _______________________________________________ > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz

25.7.2014 10:25:24 (#4)
gravatar

Petr Vejsada

<osm at propsychology.cz>
516
Ahoj, ze 40 milionů drobných parcelek se dá jedním příkazem udělat 60 multipolygonů select druh_pozemku,zpusob_vyuziti_pozemku,st_union(hranice) from ruian.rn_parcela group by druh_pozemku,zpusob_vyuziti_pozemku a z nich dalším jedním či dvěma příkazy se dá udělat N polygonů, které budou sdružovat sousedící parcely se stejným druh_pozemku a zpusob_vyuziti_pozemku st_dump(hranice) ... atd. a toto zkombinovat s LPIS. LPIS by měl mát přednost, ano, a tam, kde LPIS není, protože se majitel půdy nezaregistroval či jde o landuse mimo zemědělství, tam použit RUIAN. ? Dne Pá 25. července 2014 22:04:50, Martin Švec - OSM napsal(a): zobrazit citaci
> No, poté co jsem si to vyzkoušel jsem si prakticky jistý, že 1:1 > hromadný převod všech parcel do OSM není průchozí :-) Ty miliony > rozdrobených parcelek by nadělaly solidní binec v databázi i na mapě. > A jak říkáš, zrovna LPIS u zemědělské půdy líp kopíruje realitu. > > Spíš je to o tom hledat způsoby, jak ty RUIAN data rozumně využít, > když už jsou k dispozici. >
-- Petr

26.7.2014 12:00:29 (#5)
gravatar

Martin Švec - OSM

<osm at maatts.cz>
109
Ahoj, ano, tak bych si představoval to rozumné využití RUIANu. Akorát mě nenapadlo, že je to až tak snadné. Asi je čas naučit se PostGIS ;-) Martin On 25.7.2014 22:25, Petr Vejsada wrote: zobrazit citaci
> Ahoj, > > ze 40 milionů drobných parcelek se dá jedním příkazem udělat 60 multipolygonů > > select druh_pozemku,zpusob_vyuziti_pozemku,st_union(hranice) from > ruian.rn_parcela group by druh_pozemku,zpusob_vyuziti_pozemku > > a z nich dalším jedním či dvěma příkazy se dá udělat N polygonů, které budou > sdružovat sousedící parcely se stejným druh_pozemku a zpusob_vyuziti_pozemku > st_dump(hranice) ... atd. > > a toto zkombinovat s LPIS. LPIS by měl mát přednost, ano, a tam, kde LPIS > není, protože se majitel půdy nezaregistroval či jde o landuse mimo > zemědělství, tam použit RUIAN. ? > > Dne Pá 25. července 2014 22:04:50, Martin Švec - OSM napsal(a): > >> No, poté co jsem si to vyzkoušel jsem si prakticky jistý, že 1:1 >> hromadný převod všech parcel do OSM není průchozí :-) Ty miliony >> rozdrobených parcelek by nadělaly solidní binec v databázi i na mapě. >> A jak říkáš, zrovna LPIS u zemědělské půdy líp kopíruje realitu. >> >> Spíš je to o tom hledat způsoby, jak ty RUIAN data rozumně využít, >> když už jsou k dispozici. >> > -- > 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