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

[Talk-cz] Import of farmland from LPIS

Vlákno 22.7. - 30.7.2014, počet zpráv: 31


22.7.2014 03:28:38 (#1)
gravatar

Pavel Machek

<pavel at ucw.cz>
1056 1226
Hi! I'd like to start import of LPIS farmland database, as we have very good coverage of houses, forests and water, but farmland is very good at places and completely missing at different places. Import page is at https://wiki.openstreetmap.org/wiki/LPIS , import script is attached to this email. Best regards, Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html ------------- další část --------------- #!/usr/bin/python # http://www.lpis.cz/index.html # Mapa: # http://eagri.cz/public/app/lpisext/lpis/verejny/ # Vyzadat data na: # http://eagri.cz/public/app/lpisext/lpis/verejny/exportDat.html # Vyzadat: Pudni bloky # Souradny systrem: wgs-84 # Aha, a ty cisla katastralnich uzemi jde dokonce vykoukat z osm, jako ref u boundary=administrative, admin_level=10 # A... nefunguje jim dobre captcha, takze jde stahnout vic uzemi jednim formularem. # doubek: 631035 # babice: 600601 # for A in *.zip; do unzip $A; done import sys import shapefile from pyproj import * print "<?xml version='1.0' encoding='UTF-8'?>" print "<osm version='0.6' generator='pyshp'>" #p1 = Proj(init='esri:102067') p1 = Proj(init='EPSG:4326') p2 = Proj(init='EPSG:4326') sf = shapefile.Reader(sys.argv[1]) global field field = {} def a_field(num, name): global field i1, i2, i3, i4 = sf.fields[num] if name != i1: print 'Unexpected field name <', name field[name] = num a_field(1, 'ID_FB') a_field(3, 'NKODFB') a_field(6, 'PLATNYOD') a_field(7, 'VYMERAM') a_field(8, 'KULTURA') a_field(9, 'KULTURA_KL') a_field(10, 'KULTURAOD') a_field(11, 'EKO') a_field(21, 'VYSKA') a_field(22, 'SVAZITOST') a_field(26, 'KU_KOD') global nodeid nodeid = 0 def convert(point): lon, lat = point lon, lat = transform(p1, p2, lon, lat) #lat += 50.013082018919185-50.013834 #lon += 14.748617781670555-14.749757 return lon, lat def write_point(point): global nodeid lon, lat = convert(point) tags = '<tag k="created_by" v="shpupload"/>' tags += '<tag k="source" v="lpis"/>' nodeid -= 1 print '<node id="%d" lon="%f" lat="%f\">%s</node>' % ( nodeid, lon, lat, tags ) return nodeid def attr(shrec, name): return shrec.record[field[name]] for shrec in sf.shapeRecords(): shape = shrec.shape pts = [] for point in shape.points: pts += [ write_point(point) ] nodeid -= 1 print '<way id="%d">' % nodeid print ' <tag k="created_by" v="pyshp"/>' print ' <tag k="id_fb" v="%s"/>' % attr(shrec, 'ID_FB') print ' <tag k="source" v="lpis"/>' print ' <tag k="lpis:kultura" v="%d"/>' % attr(shrec, 'KULTURA') kul = attr(shrec, 'KULTURA') if kul == 2: print ' <tag k="landuse" v="farmland"/>' elif kul == 3: print ' <tag k="landuse" v="farmland"/><tag k="crop" v="hop"/>' elif kul == 30: print ' <tag k="landuse" v="farmland"/><tag k="crop" v="hop"/>' elif kul == 31: print ' <tag k="landuse" v="farmland"/><tag k="crop" v="hop"/>' elif kul == 41: print ' <tag k="landuse" v="vineyard"/><tag k="barrier" v="fence"/>' elif kul == 61: print ' <tag k="landuse" v="orchard"/><tag k="barrier" v="fence"/>' elif kul == 62: print ' <tag k="landuse" v="orchard"/>' elif kul == 7: print ' <tag k="landuse" v="meadow"/><tag k="meadow" v="agricultural"/>' elif kul == 71: print ' <tag k="landuse" v="meadow"/><tag k="meadow" v="agricultural"/><tag k="barrier" v="fence"/>' elif kul == 72: print ' <tag k="landuse" v="meadow"/><tag k="meadow" v="agricultural"/>' elif kul == 91: print ' <tag k="landuse" v="forest"/><tag k="barrier" v="fence"/>' elif kul == 92: print ' <tag k="landuse" v="farm"/><tag k="crop" v="vegetables"/>' elif kul == 97: print ' <tag k="landuse" v="reservoir"/>' elif kul == 98: print ' <tag k="landuse" v="scrub"/><tag k="note" v="rychle rostouci dreviny"/>' elif kul == 99: print ' <tag k="landuse" v="forest"/>' else: print ' <tag k="landuse" v="unknown_farmland_%d"/>' % kul print ' <tag k="ele" v="%d"/>' % attr(shrec, 'VYSKA') print ' <tag k="ref" v="%s"/>' % attr(shrec, 'NKODFB') for pt in pts: print ' <nd ref="%d"/>' % pt print '</way>' print '</osm>'

22.7.2014 03:38:39 (#2)
gravatar

Pavel Machek

<pavel at ucw.cz>
1056 1226
Hi! zobrazit citaci
> I'd like to start import of LPIS farmland database, as we have very > good coverage of houses, forests and water, but farmland is very good > at places and completely missing at different places. > > Import page is at https://wiki.openstreetmap.org/wiki/LPIS , import > script is attached to this email.
Is someone subscribed to the imports@ list who could forward? I already created another openstreetmap account for the import, but I don't think I want to subscribe to the mailing list. Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

22.7.2014 03:40:32 (#3)
gravatar

hanoj

<ehanoj at gmail.com>
718
zobrazit citaci
> The data are freely available to anyone by law as they were created by government agency using tax payers' money. There are no copyrights > on the data. The basic registers are created according to the Lawe No. 111/2009 Sb.. The possibility of extraction and use of data for OSM is > based on § 62 of Law 111/2009 Sb. Sorry, we are not aware of any English translations of Czech laws.
*** Ehm, a na to jsi přišel jak? ha hanoj 2014-07-22 15:28 GMT+02:00 Pavel Machek <pavel na ucw.cz>: zobrazit citaci
> Hi! > > I'd like to start import of LPIS farmland database, as we have very > good coverage of houses, forests and water, but farmland is very good > at places and completely missing at different places. > > Import page is at https://wiki.openstreetmap.org/wiki/LPIS , import > script is attached to this email. > > Best regards, > Pavel > -- > (english) http://www.livejournal.com/~pavelmachek > (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html > > _______________________________________________ > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz >

22.7.2014 05:45:14 (#4)
gravatar

Petr Vejsada

<osm at propsychology.cz>
516
Ahoj, tak jsem si "stáhl" jedno vybrané KÚ, nahrál do Postgre a zkoumám Dokonce se to někde i plusmínus kryje s landuse, co je v OSM :-). Díky za objevení tohoto zdroje dat a ukázku, jak se k nim dostat. Ty .shp jsou dobře zpracovatelné No a teď mám asi milion myšlenek, tak jen některé: - o licenci ať koumá někdo jiný, ale proč je tam ta kapča? Proč to nechtějí dát celé najednou? - 13 tisíc KÚ je strašně moc. Import, dělaný tak, jak navrhuješ, nebude hotový nikdy. Po zkušenostech z importu adres - nedělej si iluze, že se toho budou účastnit stovky lidí. Paretovo pravidlo platí i zde. - dělat to manuálně bude neuvěřitelně pracné, tedy natáhnout do JOSM tak, jak to je - merging data - myslím, že je opravdu naprosto nutné to propojit s existujícími polygony landuse, co jsou v OSM Hraju si a myšlenky jsem mezitím zapomněl ;-) http://pedro.poloha.net/osm/lpis.png http://pedro.poloha.net/osm/lpis-osmonly.png Zelené tečkované jsou pastviny, ostře žluté tečkované je školka. Ty červené čáry (nejsou to čáry, ale polygony), jsou z krajinných prvků. V přiloženém DOCu jsem se nedočetl, co ty krajinné prvky znamenají. Bylo by škoda to nevyužít. Nějaké meze či co? Chce to promyslet to spojování s existujícími daty v OSM a promyslet aktualizace. Taky vymyslet, jak to dostat celé a jak to aktualizovat. Dále - co všechno importovat - je tam flag Ekologické zemědělství, zajímá nás? V Postgisu by si s tím šlo docela pěkně hrát, naskriptovat aktualizace atd. a určitě i nějaké testování a propojování s existujícími landuse. Každopádně je to velmi zajímavé. Nějaká domluva o přebírání dat bez kapči všeho najednou? Šlo by to zpracovávat z těch mailů. Představuji i na to separátní adresu, jakmile přijde mail, hned se obsah nahraje do Postgre :-) -- Petr Dne Út 22. července 2014 15:28:38, Pavel Machek napsal(a): zobrazit citaci
> Import page is at https://wiki.openstreetmap.org/wiki/LPIS , import > script is attached to this email.

22.7.2014 10:55:23 (#5)
gravatar

Pavel Machek

<pavel at ucw.cz>
1056 1226
On Tue 2014-07-22 15:40:32, hanoj wrote: zobrazit citaci
> > The data are freely available to anyone by law as they were created by government agency using tax payers' money. There are no copyrights > > on the data. The basic registers are created according to the Lawe No. 111/2009 Sb.. The possibility of extraction and use of data for OSM is > > based on § 62 of Law 111/2009 Sb. Sorry, we are not aware of any English translations of Czech laws. > > *** Ehm, a na to jsi přišel jak?
Ty data poskytuje Ministerstvo Zemedelstvi, AFAICT, takze jak by to melo byt jinak? Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

22.7.2014 11:07:38 (#6)
gravatar

Pavel Machek

<pavel at ucw.cz>
1056 1226
Ahoj! zobrazit citaci
> tak jsem si "stáhl" jedno vybrané KÚ, nahrál do Postgre a zkoumám Dokonce se > to někde i plusmínus kryje s landuse, co je v OSM :-). Díky za objevení tohoto > zdroje dat a ukázku, jak se k nim dostat. Ty .shp jsou dobře zpracovatelné > > No a teď mám asi milion myšlenek, tak jen některé: > > - o licenci ať koumá někdo jiný, ale proč je tam ta kapča? Proč to nechtějí > dát celé najednou?
Tvrdi, ze jim to zatezuje procesor... zobrazit citaci
> - 13 tisíc KÚ je strašně moc. Import, dělaný tak, jak navrhuješ, nebude hotový > nikdy. Po zkušenostech z importu adres - nedělej si iluze, že se toho budou > účastnit stovky lidí. Paretovo pravidlo platí i zde.
Co je Paretovo pravidlo? A vic pracny nez budovy to nebude, ne? (Ty se delaj uplne rucne). Ja jsem ted mapoval landuse v okoli rican rucne, a tohle je urcite lepsi. zobrazit citaci
> - dělat to manuálně bude neuvěřitelně pracné, tedy natáhnout do JOSM tak, jak > to je
No, konverzni skript mam, nekde je pribalenej. Ziskat ty zip fily je obtizny, zbytek uz je docela dobrej. zobrazit citaci
> - merging data - myslím, že je opravdu naprosto nutné to propojit s > existujícími polygony landuse, co jsou v OSM
No, prvni otazka je proc (zas tak moc tech dat neni, a kdyz to propojime, bude se to mizerne editovat) a hned druha je jak... Mezi lesem a polem byva uzky pruh ne-pole. Podobne mezi plotem a polem. Kdyz to spojime, nebude to odpovidat realite, a pokud nekde ten pruh bude, bude dost tezky ho beznym editorem vyrobit... zobrazit citaci
> Hraju si a myšlenky jsem mezitím zapomněl ;-) > > http://pedro.poloha.net/osm/lpis.png > http://pedro.poloha.net/osm/lpis-osmonly.png > > Zelené tečkované jsou pastviny, ostře žluté tečkované je školka. > > Ty červené čáry (nejsou to čáry, ale polygony), jsou z krajinných prvků. V > přiloženém DOCu jsem se nedočetl, co ty krajinné prvky znamenají. Bylo by > škoda to nevyužít. Nějaké meze či co?
Krajinne prvky jsem zatim nezkoumal, me zajimaly hlavne louky. To bych nechal jako separatni problem. zobrazit citaci
> Chce to promyslet to spojování s existujícími daty v OSM a promyslet > aktualizace. Taky vymyslet, jak to dostat celé a jak to aktualizovat. Dále - > co všechno importovat - je tam flag Ekologické zemědělství, zajímá > nás?
No, dostat to cele by bylo cool. Jinak ta captcha se nemusi zadavat, staci jen menit cislo, takze to neni tak strasny. Aktualizace.. no vzhledem k captche bych ocekaval ze se udelaj rucne kdyz se dana oblast zmeni. Ekologicke zemedelstvi -- snazil jsem se googlit, ale zda se ze na to nemame tagy. Co jsem zkusil a nemusi byt dobry napad -- dost uzitecny by bylo kdyby z toho slo vykoukat automaticky barrier=fence. zobrazit citaci
> V Postgisu by si s tím šlo docela pěkně hrát, naskriptovat aktualizace atd. a > určitě i nějaké testování a propojování s existujícími landuse. > > Každopádně je to velmi zajímavé. Nějaká domluva o přebírání dat bez kapči > všeho najednou? Šlo by to zpracovávat z těch mailů. Představuji i na to > separátní adresu, jakmile přijde mail, hned se obsah nahraje do > Postgre :-)
Jestli bys byl schopen vytvorit adresu ktera to dostane a rozbali, to by bylo fajn... Jinak na nezanedbatelny casti uzemi uz nekdo ty pole udelal, a jsem si celkem jistej ze s timhle se to bude delat lip nez uplne rucne. Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

23.7.2014 12:05:13 (#7)
gravatar

Petr Vejsada

<osm at propsychology.cz>
516
Ahoj, Dne Út 22. července 2014 23:07:38, Pavel Machek napsal(a): zobrazit citaci
> > - 13 tisíc KÚ je strašně moc. Import, dělaný tak, jak navrhuješ, nebude > > hotový nikdy. Po zkušenostech z importu adres - nedělej si iluze, že se > > toho u účastnit stovky lidí. Paretovo pravidlo platí i zde. > > Co je Paretovo pravidlo? A vic pracny nez budovy to nebude, ne? (Ty se > delaj uplne rucne).
to pravidlo říká, že v rámci tohoto Matrixu 20 % lidí udělá 80 % práce a naopak, 80 % lidí udělá jen 20 % práce toho celku (populace). Na importu adres pracuje/pracovalo dohromady i se mnou 14 lidí. Z toho 8 zpracovalo více než 10 tis. adres. Víc pracné než budovy to nebude, jistěže ne. U budov taky IMO nemáme šanci ve stávající sestavě ten "import" dokončit, kdežto u té zemědělské půdy to reálné je. Jen si nemyslím, že je to reálné ručně. zobrazit citaci
> Ja jsem ted mapoval landuse v okoli rican rucne, a tohle je urcite lepsi. > > > - dělat to manuálně bude neuvěřitelně pracné, tedy natáhnout do JOSM tak, > > jak to je > > No, konverzni skript mam, nekde je pribalenej. Ziskat ty zip fily je > obtizny, zbytek uz je docela dobrej.
nj, rozumím, docela v poho, když děláš okolí Říčan. Kolikrát se Říčany vejdou do ČR? ;-) zobrazit citaci
> > - merging data - myslím, že je opravdu naprosto nutné to propojit s > > existujícími polygony landuse, co jsou v OSM > > No, prvni otazka je proc (zas tak moc tech dat neni, a kdyz to > propojime, bude se to mizerne editovat) a hned druha je jak... > > Mezi lesem a polem byva uzky pruh ne-pole. Podobne mezi plotem a > polem. Kdyz to spojime, nebude to odpovidat realite, a pokud nekde ten > pruh bude, bude dost tezky ho beznym editorem vyrobit...
Já jsem tím nemyslel lepit polygony k sobě tam, kde v realitě nejsou nalepené. Myslel jsem tím třeba to, že tam, kde je les, rybník, řeka, továrna, dálnice, obytná zóna a kdoví co ještě, takže tam nemůže být zároveň zemědělská půda. Toto ocheckovat strojově a v případě kolize řešit ručně, nějak podobně, jako se to dělá u adres. Myslím, že se to celkem osvědčilo. Pak je otázkou, co se vzájemně sousedícími polygony téhož. I na tom obrázku, co jsem dával odkaz, je vidět, že pastvina sousedí s pastvinou. Měly by to tedy být dvě pastviny vedle sebe nebo to spojit do jedné pastviny? To záleží na tom, jak by se naplánovala budoucí aktualizace. Nestudoval jsem to tak podrobně, abych věděl, zda tam je nějaké datum, jak je to se všelijakými ID, zda máme možnost se dovědět, zda byl nějaký polygon smazán a jak se to dovíme atd. Prostě aby to bylo dlouhodobě funkční a aktuální. Co třeba s tímhle? http://pedro.poloha.net/osm/pastvina-les.png Možná je tam mezi lesem a pastvinou opravdu území ničeho, možná ne, ale je tam průnik pastviny a lesa, což je blbost. To je třeba řešit. A tohle už je komplikovanější: http://pedro.poloha.net/osm/pastvina-les-2.png Je zřejmé, že něco z toho není přesné. Jak rozhodnout co to je? Les nebo pastvina? zobrazit citaci
> Krajinne prvky jsem zatim nezkoumal, me zajimaly hlavne louky. To bych > nechal jako separatni problem.
Jo, to můžeme. zobrazit citaci
> > Chce to promyslet to spojování s existujícími daty v OSM a promyslet > > aktualizace. Taky vymyslet, jak to dostat celé a jak to aktualizovat. Dále > > - co všechno importovat - je tam flag Ekologické zemědělství, zajímá nás? > > No, dostat to cele by bylo cool. Jinak ta captcha se nemusi zadavat, > staci jen menit cislo, takze to neni tak strasny.
Takže je triviální naskriptovat něco, co si zažádá o všechna KÚ, ale to bych nedělal; spíš bych se domluvil. Nejspíš by rychle opravili kapču a zabanovali nás. zobrazit citaci
> Aktualizace.. no vzhledem k captche bych ocekaval ze se udelaj rucne > kdyz se dana oblast zmeni.
? nebo ten strojový čas/zátěž, ne? zobrazit citaci
> Jestli bys byl schopen vytvorit adresu ktera to dostane a rozbali, to > by bylo fajn...
To není problém a můžu to rovnou převést tím tvým skriptem. Když jsme u toho - kolik lidí má na Windows nainstalovaný Python (z toho kolik lidí má pyproj a shapefile), respektive kolik lidí je schopno si toto nainstalovat (na Windows?) V této konferenci jistě více než obvykle ;-), ale asi ne úplně každý, kdo by chtěl mapovat. No a když takový mail udělám, tak co mám pak s těmi rozbalenými (zkonvertovanými) soubory dělat? zobrazit citaci
> Jinak na nezanedbatelny casti uzemi uz nekdo ty pole udelal, a jsem si > celkem jistej ze s timhle se to bude delat lip nez uplne rucne.
To určitě ano, jen pořád pochybuji, že to bude stačit k dokončení celého importu v nějakém konečném čase. -- -p-

23.7.2014 01:17:02 (#8)
gravatar

Pavel Machek

<pavel at ucw.cz>
1056 1226
Ahoj! zobrazit citaci
> to pravidlo říká, že v rámci tohoto Matrixu 20 % lidí udělá 80 % práce a > naopak, 80 % lidí udělá jen 20 % práce toho celku (populace). > > Na importu adres pracuje/pracovalo dohromady i se mnou 14 lidí. Z toho 8 > zpracovalo více než 10 tis. adres.
Dobra prace :-). zobrazit citaci
> Víc pracné než budovy to nebude, jistěže ne. U budov taky IMO nemáme šanci ve > stávající sestavě ten "import" dokončit, kdežto u té zemědělské půdy to reálné > je. Jen si nemyslím, že je to reálné ručně.
Ja nevim co je "rucne". Kdyz mam ty zip fily, tak to muzu skriptem rozbalit, konvertovat do osm, kouknout na to v josm, a upload. Rozumna davka je imo tak 10 oblasti najednou.. zobrazit citaci
> > Ja jsem ted mapoval landuse v okoli rican rucne, a tohle je urcite lepsi. > > > > > - dělat to manuálně bude neuvěřitelně pracné, tedy natáhnout do JOSM tak, > > > jak to je > > > > No, konverzni skript mam, nekde je pribalenej. Ziskat ty zip fily je > > obtizny, zbytek uz je docela dobrej. > > nj, rozumím, docela v poho, když děláš okolí Říčan. Kolikrát se Říčany vejdou > do ČR? ;-)
20x. Tak to jsem cekal horsi, to bude za chvili hotovy :-). (30km okoli Rican me zajima natolik ze jsem zacal ty pole mapovat rucne klikanim podle bingu... a to je fakt dost otrava.) zobrazit citaci
> Já jsem tím nemyslel lepit polygony k sobě tam, kde v realitě nejsou nalepené. > Myslel jsem tím třeba to, že tam, kde je les, rybník, řeka, továrna, > dálnice,
Aha, ok. zobrazit citaci
> obytná zóna a kdoví co ještě, takže tam nemůže být zároveň zemědělská půda. > Toto ocheckovat strojově a v případě kolize řešit ručně, nějak podobně, jako > se to dělá u adres. Myslím, že se to celkem osvědčilo.
Predstavoval jsem si, ze by na to mel rvat josm-i validator, bohuzel nerve. (A teda zrovna prekryv lesa a pastviny se dela -- zvirata se pak maj kde schovat). zobrazit citaci
> Pak je otázkou, co se vzájemně sousedícími polygony téhož. I na tom obrázku, > co jsem dával odkaz, je vidět, že pastvina sousedí s pastvinou. Měly by to > tedy být dvě pastviny vedle sebe nebo to spojit do jedné pastviny? > To záleží
Urcite nespojovat. Kdyz jsou 2 pastviny, je mezi nima plot / elektrickej ohradnik. Kdyz jsou 2 pole, na kazdym muze rust neco jinyho, a i kdyz tam nahodou roste to samy, je mensi prasarna jit po tom uzkym pruhu mezi nima. zobrazit citaci
> na tom, jak by se naplánovala budoucí aktualizace. Nestudoval jsem to tak > podrobně, abych věděl, zda tam je nějaké datum, jak je to se všelijakými ID, > zda máme možnost se dovědět, zda byl nějaký polygon smazán a jak se to dovíme > atd. Prostě aby to bylo dlouhodobě funkční a aktuální.
ID tam jsou, delam z nich ref=. Je tam datum aktualizace. zobrazit citaci To muze byt i v poradku.. i kdyz nejspis neni, to se nam proste vymstila generalizace lesu. Posunout les podle ortofota.. zobrazit citaci
> Možná je tam mezi lesem a pastvinou opravdu území ničeho, možná ne, ale je tam > průnik pastviny a lesa, což je blbost. To je třeba řešit. > > A tohle už je komplikovanější: > http://pedro.poloha.net/osm/pastvina-les-2.png > > Je zřejmé, že něco z toho není přesné. Jak rozhodnout co to je? Les nebo > pastvina?
Les je blbe, delala se tam moc silna generalizace. zobrazit citaci
> > > Chce to promyslet to spojování s existujícími daty v OSM a promyslet > > > aktualizace. Taky vymyslet, jak to dostat celé a jak to aktualizovat. Dále > > > - co všechno importovat - je tam flag Ekologické zemědělství, zajímá nás?
(Jinak k tomu ekologickymu zemedelstvi jsem nenasel odpovidajici tag na osm, takze to prozatim ignoruju). zobrazit citaci
> > No, dostat to cele by bylo cool. Jinak ta captcha se nemusi zadavat, > > staci jen menit cislo, takze to neni tak strasny. > > Takže je triviální naskriptovat něco, co si zažádá o všechna KÚ, ale to bych > nedělal; spíš bych se domluvil. Nejspíš by rychle opravili kapču a zabanovali > nás.
Umis to naskriptovat? Me by to chvili trvalo. Tech uzemi je jen par tisic, kdyz to udelame o weekendu, myslim ze se nebudou s banovanim obtezovat, zvlast kdyz tam dame par sekund delay aby to ten server uplne neutavilo. I kdyz, dohoda by asi byla lepsi :-). zobrazit citaci
> > Jestli bys byl schopen vytvorit adresu ktera to dostane a rozbali, to > > by bylo fajn... > > To není problém a můžu to rovnou převést tím tvým skriptem. Když jsme u toho - > kolik lidí má na Windows nainstalovaný Python (z toho kolik lidí má pyproj a > shapefile), respektive kolik lidí je schopno si toto nainstalovat (na Windows?) > V této konferenci jistě více než obvykle ;-), ale asi ne úplně každý, kdo by > chtěl mapovat.
Ten pyproj tam nakonec neni potreba. shapefile se instaluje pomoci pip-u. A sorry, Windows neznam a nedoporucuju. zobrazit citaci
> No a když takový mail udělám, tak co mám pak s těmi rozbalenými > (zkonvertovanými) soubory dělat?
Dat na ftp? Pro zacatek by se hodilo nezkonvertovany, protoze skript bude jeste potrebovat par uprav. zobrazit citaci
> > Jinak na nezanedbatelny casti uzemi uz nekdo ty pole udelal, a jsem si > > celkem jistej ze s timhle se to bude delat lip nez uplne rucne. > > To určitě ano, jen pořád pochybuji, že to bude stačit k dokončení celého > importu v nějakém konečném čase.
No, kdyz jsem s osm zacinal, byla tu velka bila plocha. Kreslit pole rucne je fakt otrava, a lidi to stejne delaj. Asi bych se toho nebal :-). [Samozrejme by byla moznost to tam proste nahrat vsechno, s tim ze kazdy opravi prekryvy (etc) na uzemi ktery ho zajima. Rychlejc bychom dosahli kompletniho pokryti za cenu kvality dat... povetsinou kvality okolnich lesu.] Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

23.7.2014 02:07:16 (#9)
gravatar

Martin Švec - OSM

<osm at maatts.cz>
109
Ahoj, cvičně jsem si taky zkusil tři KÚ, první dojmy: * Skript by měl slučovat duplicitní uzly (triviální). * Skript by měl eliminovat uzel v cestě, pokud je totožný jako předcházející uzel (triviální). * Musí se ručně řešit multipolygony. * Spousta sousedících polygonů má překryvy, i v rámci jednoho KU. * LPIS obsahuje jen zemědělskou půdu, takže v mapě vznikají "ostrůvky" polygonů místo souvislé plochy. * Podezřele hodně luk je označeno jako oplocených (kul. 71) -- skript předpokládá, že každá pastvina má ohradník/plot? * Jak už zaznělo, budou potíže na rozhraních s existujícími lesy. * LPIS na řadě míst nekopíruje hranice parcel, takže budou nesoulady s daty RUIANu, plochami kreslenými podle KM atd. * Od pohledu mi přijde, že parcely RUIANu by byly kvalitnější zdroj pro plošné mapování. Na druhou stranu, už teď skript ušetří mraky ručního obkreslování, stačí jen v JOSM doopravit problémy. :-) (V příloze je skript s přidaným primitivním odstraňováním duplicitních uzlů. Nicméně, Python vůbec neznám a z jeho syntaxe dostávám osypky, takže zbytek nechám pythonistům :-)) Martin On 22.7.2014 15:28, Pavel Machek wrote: zobrazit citaci
> Hi! > > I'd like to start import of LPIS farmland database, as we have very > good coverage of houses, forests and water, but farmland is very good > at places and completely missing at different places. > > Import page is at https://wiki.openstreetmap.org/wiki/LPIS , import > script is attached to this email. > > Best regards, > Pavel > > > _______________________________________________ > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz
------------- další část --------------- #!/usr/bin/python # http://www.lpis.cz/index.html # Mapa: # http://eagri.cz/public/app/lpisext/lpis/verejny/ # Vyzadat data na: # http://eagri.cz/public/app/lpisext/lpis/verejny/exportDat.html # Vyzadat: Pudni bloky # Souradny systrem: wgs-84 # Aha, a ty cisla katastralnich uzemi jde dokonce vykoukat z osm, jako ref u boundary=administrative, admin_level=10 # A... nefunguje jim dobre captcha, takze jde stahnout vic uzemi jednim formularem. # doubek: 631035 # babice: 600601 # for A in *.zip; do unzip $A; done import sys import shapefile from pyproj import * print "<?xml version='1.0' encoding='UTF-8'?>" print "<osm version='0.6' generator='pyshp'>" #p1 = Proj(init='esri:102067') p1 = Proj(init='EPSG:4326') p2 = Proj(init='EPSG:4326') sf = shapefile.Reader(sys.argv[1]) global field field = {} def a_field(num, name): global field i1, i2, i3, i4 = sf.fields[num] if name != i1: print 'Unexpected field name <', name field[name] = num a_field(1, 'ID_FB') a_field(3, 'NKODFB') a_field(6, 'PLATNYOD') a_field(7, 'VYMERAM') a_field(8, 'KULTURA') a_field(9, 'KULTURA_KL') a_field(10, 'KULTURAOD') a_field(11, 'EKO') a_field(21, 'VYSKA') a_field(22, 'SVAZITOST') a_field(26, 'KU_KOD') global nodeid nodeid = 0 global nodemap nodemap = {} def convert(point): lon, lat = point lon, lat = transform(p1, p2, lon, lat) #lat += 50.013082018919185-50.013834 #lon += 14.748617781670555-14.749757 return lon, lat def write_point(point): global nodeid global nodemap lon, lat = convert(point) nodekey = '%f/%f' % (lon, lat) if nodekey in nodemap: return nodemap[nodekey] tags = '<tag k="created_by" v="shpupload"/>' tags += '<tag k="source" v="lpis"/>' nodeid -= 1 print '<node id="%d" lon="%f" lat="%f\">%s</node>' % ( nodeid, lon, lat, tags ) nodemap[nodekey] = nodeid; return nodeid def attr(shrec, name): return shrec.record[field[name]] for shrec in sf.shapeRecords(): shape = shrec.shape prevpt = 6666 pts = [] for point in shape.points: curpt = write_point(point) if prevpt != curpt: pts += [ curpt ] prevpt = curpt nodeid -= 1 print '<way id="%d">' % nodeid print ' <tag k="created_by" v="pyshp"/>' print ' <tag k="id_fb" v="%s"/>' % attr(shrec, 'ID_FB') print ' <tag k="source" v="lpis"/>' print ' <tag k="lpis:kultura" v="%d"/>' % attr(shrec, 'KULTURA') kul = attr(shrec, 'KULTURA') if kul == 2: print ' <tag k="landuse" v="farmland"/>' elif kul == 3: print ' <tag k="landuse" v="farmland"/><tag k="crop" v="hop"/>' elif kul == 30: print ' <tag k="landuse" v="farmland"/><tag k="crop" v="hop"/>' elif kul == 31: print ' <tag k="landuse" v="farmland"/><tag k="crop" v="hop"/>' elif kul == 41: print ' <tag k="landuse" v="vineyard"/><tag k="barrier" v="fence"/>' elif kul == 61: print ' <tag k="landuse" v="orchard"/><tag k="barrier" v="fence"/>' elif kul == 62: print ' <tag k="landuse" v="orchard"/>' elif kul == 7: print ' <tag k="landuse" v="meadow"/><tag k="meadow" v="agricultural"/>' elif kul == 71: print ' <tag k="landuse" v="meadow"/><tag k="meadow" v="agricultural"/><tag k="barrier" v="fence"/>' elif kul == 72: print ' <tag k="landuse" v="meadow"/><tag k="meadow" v="agricultural"/>' elif kul == 91: print ' <tag k="landuse" v="forest"/><tag k="barrier" v="fence"/>' elif kul == 92: print ' <tag k="landuse" v="farm"/><tag k="crop" v="vegetables"/>' elif kul == 97: print ' <tag k="landuse" v="reservoir"/>' elif kul == 98: print ' <tag k="landuse" v="scrub"/><tag k="note" v="rychle rostouci dreviny"/>' elif kul == 99: print ' <tag k="landuse" v="forest"/>' else: print ' <tag k="landuse" v="unknown_farmland_%d"/>' % kul print ' <tag k="ele" v="%d"/>' % attr(shrec, 'VYSKA') print ' <tag k="ref" v="%s"/>' % attr(shrec, 'NKODFB') for pt in pts: print ' <nd ref="%d"/>' % pt print '</way>' print '</osm>'

23.7.2014 02:43:20 (#10)
gravatar

Petr Vejsada

<osm at propsychology.cz>
516
Ahoj, Dne St 23. července 2014 01:17:02, Pavel Machek napsal(a): zobrazit citaci
> Ja nevim co je "rucne". Kdyz mam ty zip fily, tak to muzu skriptem > rozbalit, konvertovat do osm, kouknout na to v josm, a upload. Rozumna > davka je imo tak 10 oblasti najednou..
"Ručně" říkám přesně tomu, co tady popisuješ ;-). "Automaticky" říkám tomu, že pokud se na tom shodneme, upraví se polygon lesa sám bez hledání a klikání. Pořád mám v qgisu natažené to jedno KÚ a vidím tam velké překryvy v landuse=residential, vidím tam díry ve školce (skript vůbec neřeší), školka zasahuje do obytné oblasti, takže v těch dírách je co? Obytná oblast? Nevím, žádný dům tam nestojí. Oproti tomu jsem si všiml, že v louce je díra pro sloup elektrického vedení a dokonce to sedí na sloup, co je v OSM, super :-). Je tam i mezera pro silnici, ve které je silnice, taky super, ale na jednom kousku silnice leze do louky, takže někde bude nějaká nepřesnost. Pak tam vidím oblasti, které se částečně shodují s daty v OSM - v OSM je zemědělská půda. Nekryje se to přesně (tím nemyslím centimetry), ale prostě se to liší. Co s tím? zobrazit citaci
> > nj, rozumím, docela v poho, když děláš okolí Říčan. Kolikrát se Říčany > > vejdou do ČR? ;-) > > 20x. Tak to jsem cekal horsi, to bude za chvili hotovy :-). > > (30km okoli Rican me zajima natolik ze jsem zacal ty pole mapovat > rucne klikanim podle bingu... a to je fakt dost otrava.)
Aha, kruh 60km v průměru, no odsejpá ti to. zobrazit citaci
> > Co třeba s tímhle? > > http://pedro.poloha.net/osm/pastvina-les.png > > To muze byt i v poradku.. i kdyz nejspis neni, to se nam proste > vymstila generalizace lesu. Posunout les podle ortofota..
No, myslím, že buďto je les nebo je pastvina. Obojímu dohromady se říká step ;-). zobrazit citaci
> (Jinak k tomu ekologickymu zemedelstvi jsem nenasel odpovidajici tag > na osm, takze to prozatim ignoruju).
No OSM je obecná mapa, tak tam nemusí být úplně všechno zobrazit citaci
> > Takže je triviální naskriptovat něco, co si zažádá o všechna KÚ, ale to > > bych nedělal; spíš bych se domluvil. Nejspíš by rychle opravili kapču a > > zabanovali nás. > > Umis to naskriptovat? Me by to chvili trvalo. Tech uzemi je jen par > tisic, kdyz to udelame o weekendu, myslim ze se nebudou s banovanim > obtezovat, zvlast kdyz tam dame par sekund delay aby to ten server > uplne neutavilo.
Nezkoumal jsem, podle toho, jak to popisuješ, tak je to jen načíst z DB čísla katastrálních území a pro každé poslat jeden POST s příslušnými parametry na URL, kde se ty požadavky zadávají, ne? zobrazit citaci
> I kdyz, dohoda by asi byla lepsi :-).
Možná jsem příliš ovlivněný prací na těch adresách, jen mi přijde fajn možnost zjistit, že je někde něco aktualizované a hned to promítnout do OSM, takže bych nechtěl riskovat opravu captcha. zobrazit citaci
> > No a když takový mail udělám, tak co mám pak s těmi rozbalenými > > (zkonvertovanými) soubory dělat? > > Dat na ftp? Pro zacatek by se hodilo nezkonvertovany, protoze skript > bude jeste potrebovat par uprav.
No to by šlo, jen pořád nevím, k čemu by to bylo? Aby si to lidi mohli stáhnout z FTP namísto vyplnění toho formuláře? Nevíme, kdy si kdo co stáhne a je zbytečné, aby soubor ležel půl roku na FTP, když po půl roce bude jiný, aktuální.? zobrazit citaci
> [Samozrejme by byla moznost to tam proste nahrat vsechno, s tim ze > kazdy opravi prekryvy (etc) na uzemi ktery ho zajima. Rychlejc bychom > dosahli kompletniho pokryti za cenu kvality dat... povetsinou kvality > okolnich lesu.]
... a bude landuse=farmland over landuse=farmland :-( -- -p-

23.7.2014 02:52:01 (#11)
gravatar

Petr Vejsada

<osm at propsychology.cz>
516
Ahoj, Dne St 23. července 2014 02:07:16, Martin Švec - OSM napsal(a): zobrazit citaci
> * LPIS na řadě míst nekopíruje hranice parcel, takže budou nesoulady s daty > RUIANu, plochami kreslenými podle KM atd. * Od pohledu mi přijde, že > parcely RUIANu by byly kvalitnější zdroj pro plošné mapování.
Tohle chce prozkoumat - v RUIAN u parcel je druh_pozemku a zpusob_vyuziti_pozemku , každopádně je tam nejen zemědělská půda, ale i např. hřbitov, sportovní plochy, vodní nádrže, dopravní komunikace, ... Hranice tvoří hranice parcel, jasně. Tak to chce zjistit, jak je to v těch případech, kdy LPIS nekopíruje hranice parcel. V realitě je přeci možné, abych neměl pole na celé parcele, ale jen na části. Případně oba zdroje dat spojit dohromady.. -- -p-

23.7.2014 09:29:54 (#12)
gravatar

hanoj

<ehanoj at gmail.com>
718
zobrazit citaci
> On Tue 2014-07-22 15:40:32, hanoj wrote: >> > The data are freely available to anyone by law as they were created by government agency using tax payers' money. There are no copyrights >> > on the data. The basic registers are created according to the Lawe No. 111/2009 Sb.. The possibility of extraction and use of data for OSM is >> > based on § 62 of Law 111/2009 Sb. Sorry, we are not aware of any English translations of Czech laws. >> >> *** Ehm, a na to jsi přišel jak? > > Ty data poskytuje Ministerstvo Zemedelstvi, AFAICT, takze jak by to > melo byt jinak?
*** Zacal bych tim, ze si ten odkazovany zakon prectes. Ve skutecnosti je to uplne jinak: http://www.zakonyprolidi.cz/cs/1997-252#p3ab-2 http://www.zakonyprolidi.cz/cs/2000-121#p3-1-a hanoj

23.7.2014 12:06:21 (#13)
gravatar

Martin Švec - OSM

<osm at maatts.cz>
109
Dne 23.7.2014 2:43, Petr Vejsada napsal(a): zobrazit citaci
> Pořád mám v qgisu natažené to jedno KÚ a vidím tam velké překryvy v > landuse=residential, vidím tam díry ve školce (skript vůbec neřeší), školka > zasahuje do obytné oblasti, takže v těch dírách je co? Obytná oblast? Nevím, > žádný dům tam nestojí. Oproti tomu jsem si všiml, že v louce je díra pro sloup > elektrického vedení a dokonce to sedí na sloup, co je v OSM, super :-). Je tam > i mezera pro silnici, ve které je silnice, taky super, ale na jednom kousku > silnice leze do louky, takže někde bude nějaká nepřesnost. Pak tam vidím > oblasti, které se částečně shodují s daty v OSM - v OSM je zemědělská půda. > Nekryje se to přesně (tím nemyslím centimetry), ale prostě se to liší. Co s > tím?
Z toho co jsem viděl mi přijde, že LPIS data vznikly stejně jako OSM -- ručně se obkreslovaly KM a ortofoto položené přes sebe. Na volných plochách polygony LPISu víceméně odpovídají KM (s odchylkou centimetrů až metrů). Podél hranic lesa, kolem vesnic, kolem objektů uprostřed pole apod. je polygon nakreslený "od oka" podle ortofoto snímků. Což odpovídá podstatě OSM, ale bude se to špatně přilepovat k jiným datům. Martin

24.7.2014 12:01:04 (#14)
gravatar

Pavel Machek

<pavel at ucw.cz>
1056 1226
Ahoj! zobrazit citaci
> cvičně jsem si taky zkusil tři KÚ, první dojmy: > > * Skript by měl slučovat duplicitní uzly (triviální). > * Skript by měl eliminovat uzel v cestě, pokud je totožný jako > předcházející uzel (triviální).
Diky za skript! zobrazit citaci
> * Musí se ručně řešit multipolygony.
Co jsem videl na zernovce, tak maji zajimavy pristup. Uprostred pastviny je pristresek, a oni ho rucne z toho polygonu vyndali. Tohle bych radeji smazal / udelal z toho pristresku building=yes. zobrazit citaci
> * Spousta sousedících polygonů má překryvy, i v rámci jednoho KU.
Jako ze LPIS se prekryva sam se sebou? zobrazit citaci
> * LPIS obsahuje jen zemědělskou půdu, takže v mapě vznikají > "ostrůvky" polygonů místo souvislé plochy.
To se asi ocekava, ne? Lesy a vodu uz mame, tak ted pole a louky. zobrazit citaci
> * Podezřele hodně luk je označeno jako oplocených (kul. 71) -- > skript předpokládá, že každá pastvina má ohradník/plot?
Tak, ale ted jsem to vypnul... ono to moc nesedi. zobrazit citaci
> * LPIS na řadě míst nekopíruje hranice parcel, takže budou nesoulady > s daty RUIANu, plochami kreslenými podle KM atd.
Mozna ze to pole neni "do kraje" parcely? zobrazit citaci
> * Od pohledu mi přijde, že parcely RUIANu by byly kvalitnější zdroj >pro plošné mapování.
Jsou v RUIANu udaje co kde je? Je nekde konverzni skript/example vystup? zobrazit citaci
> Na druhou stranu, už teď skript ušetří mraky ručního obkreslování, > stačí jen v JOSM doopravit problémy. :-)
Presne. zobrazit citaci
> (V příloze je skript s přidaným primitivním odstraňováním > duplicitních uzlů. Nicméně, Python vůbec neznám a z jeho syntaxe > dostávám osypky, takže zbytek nechám pythonistům :-))
Tak doufam ze osypky rychle pominou :-). Dik, Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

24.7.2014 12:14:51 (#15)
gravatar

Pavel Machek

<pavel at ucw.cz>
1056 1226
On Wed 2014-07-23 09:29:54, hanoj wrote: zobrazit citaci
> > On Tue 2014-07-22 15:40:32, hanoj wrote: > >> > The data are freely available to anyone by law as they were created by government agency using tax payers' money. There are no copyrights > >> > on the data. The basic registers are created according to the Lawe No. 111/2009 Sb.. The possibility of extraction and use of data for OSM is > >> > based on § 62 of Law 111/2009 Sb. Sorry, we are not aware of any English translations of Czech laws. > >> > >> *** Ehm, a na to jsi přišel jak? > > > > Ty data poskytuje Ministerstvo Zemedelstvi, AFAICT, takze jak by to > > melo byt jinak? > *** Zacal bych tim, ze si ten odkazovany zakon prectes. Ve skutecnosti > je to uplne jinak: > > http://www.zakonyprolidi.cz/cs/1997-252#p3ab-2 > http://www.zakonyprolidi.cz/cs/2000-121#p3-1-a
Kazdopadne, public domain to je, ze? Bude tohle jasnejsi? The data are freely available, as they are excluded from copyright law protection under "publicly available list" exception. http://www.zakonyprolidi.cz/cs/2000-121#p3-1-a http://www.zakonyprolidi.cz/cs/1997-252#p3ab-2 . U.S. equivalent of this would be "public domain". Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

24.7.2014 01:03:24 (#16)
gravatar

Pavel Machek

<pavel at ucw.cz>
1056 1226
Ahoj! zobrazit citaci
> Dne St 23. července 2014 01:17:02, Pavel Machek napsal(a): > > > Ja nevim co je "rucne". Kdyz mam ty zip fily, tak to muzu skriptem > > rozbalit, konvertovat do osm, kouknout na to v josm, a upload. Rozumna > > davka je imo tak 10 oblasti najednou.. > > "Ručně" říkám přesně tomu, co tady popisuješ ;-). "Automaticky" říkám tomu, že > pokud se na tom shodneme, upraví se polygon lesa sám bez hledání a klikání.
No, ja bych rad aby import nejak probehl, aspon pro okoli Rican (:-), protoze me nebavi rucne obklikavat louky. Myslim, ze "import rucne" je zvladnutelny, a to "automaticky" bych asi programoval dost nerad -- neumim postgis a neznam ty spravny knihovny. Na druhou stranu budu rad kdyz to nekdo v tom postgisu udela za me. zobrazit citaci
> Pořád mám v qgisu natažené to jedno KÚ a vidím tam velké překryvy v > landuse=residential, vidím tam díry ve školce (skript vůbec neřeší), školka > zasahuje do obytné oblasti, takže v těch dírách je co? Obytná oblast? Nevím, > žádný dům tam nestojí. Oproti tomu jsem si všiml, že v louce je díra pro sloup > elektrického vedení a dokonce to sedí na sloup, co je v OSM, super :-). Je tam > i mezera pro silnici, ve které je silnice, taky super, ale na jednom kousku > silnice leze do louky, takže někde bude nějaká nepřesnost. Pak tam vidím > oblasti, které se částečně shodují s daty v OSM - v OSM je zemědělská půda. > Nekryje se to přesně (tím nemyslím centimetry), ale prostě se to liší. Co s > tím?
No, nahrazovani existujicich dat by se asi melo delat rucne -- z ucty k tem co ty data vyrobili. zobrazit citaci
> > > nj, rozumím, docela v poho, když děláš okolí Říčan. Kolikrát se Říčany > > > vejdou do ČR? ;-) > > > > 20x. Tak to jsem cekal horsi, to bude za chvili hotovy :-). > > > > (30km okoli Rican me zajima natolik ze jsem zacal ty pole mapovat > > rucne klikanim podle bingu... a to je fakt dost otrava.) > > Aha, kruh 60km v průměru, no odsejpá ti to.
Ne, nerikam ze to tam je hotovy -- klikat podle Bingu by to byla otrava -- ale tam bych to urcite chtel, a i to tam trochu znam. zobrazit citaci
> > > Takže je triviální naskriptovat něco, co si zažádá o všechna KÚ, ale to > > > bych nedělal; spíš bych se domluvil. Nejspíš by rychle opravili kapču a > > > zabanovali nás. > > > > Umis to naskriptovat? Me by to chvili trvalo. Tech uzemi je jen par > > tisic, kdyz to udelame o weekendu, myslim ze se nebudou s banovanim > > obtezovat, zvlast kdyz tam dame par sekund delay aby to ten server > > uplne neutavilo. > > Nezkoumal jsem, podle toho, jak to popisuješ, tak je to jen načíst z DB čísla > katastrálních území a pro každé poslat jeden POST s příslušnými parametry na > URL, kde se ty požadavky zadávají, ne?
Ne, tak jednoduchy to neni, tohle bych umel. Je tam nejakej javaskript; kdyz se stranka reloaduje, da to novou capchu. Ale bez reloadu (rucne) tam jde dat jiny cislo, kliknout "provest export", zmenit cislo KU, "provest export"... Ale lynxem to nenaskriptuju. zobrazit citaci
> > I kdyz, dohoda by asi byla lepsi :-). > > Možná jsem příliš ovlivněný prací na těch adresách, jen mi přijde fajn možnost > zjistit, že je někde něco aktualizované a hned to promítnout do OSM, takže > bych nechtěl riskovat opravu captcha.
No, uvadej tam helpdesk na mze.cz ... zobrazit citaci
> > > No a když takový mail udělám, tak co mám pak s těmi rozbalenými > > > (zkonvertovanými) soubory dělat? > > > > Dat na ftp? Pro zacatek by se hodilo nezkonvertovany, protoze skript > > bude jeste potrebovat par uprav. > > No to by šlo, jen pořád nevím, k čemu by to bylo? Aby si to lidi mohli > stáhnout z FTP namísto vyplnění toho formuláře? Nevíme, kdy si kdo co stáhne > a je zbytečné, aby soubor ležel půl roku na FTP, když po půl roce bude jiný, > aktuální.?
Aby lidi co vypnej formular nemuseli jeste zarizovat emailovou adresu s rozbalovanim. ... ok, asi by bylo jednodussi rozbalovat to z mbox formatu. zobrazit citaci
> > [Samozrejme by byla moznost to tam proste nahrat vsechno, s tim ze > > kazdy opravi prekryvy (etc) na uzemi ktery ho zajima. Rychlejc bychom > > dosahli kompletniho pokryti za cenu kvality dat... povetsinou kvality > > okolnich lesu.] > > ... a bude landuse=farmland over landuse=farmland :-(
No, to by se doladilo rucne :-). Na druhou stranu importovat automaticky vse, co se nekryje vubec s nicim, by mohl byt dobry zacatek. Jinak zda se ze zmeny se daj sehnat ze: http://eagri.cz/public/web/mze/farmar/elektronicka-vymena-dat/prehled-vystavenych-sluzeb/ A zda-se ze je i nejaky wfs server, takze to treba nebude nutny tahat pres web+mail... http://eagri.cz/public/web/mze/farmar/LPIS/novinky/nove-wms-a-wfs-sluzby.html Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

24.7.2014 02:11:55 (#17)
gravatar

Petr Vejsada

<osm at propsychology.cz>
516
Ahoj, ani nic nequotuji, jen co se mi zdá, tak podle toho zákona, co sem dával odkaz hanoj, jsou v LPIS jen ty pozemky, u kterých to jejich uživatelé chtějí. Kdo tam nechce být, tak tam není. Chápu správně? Udělal jsem narychlo vizualizaci RUIAN, pouze druh pozemku. V tabulkách je ještě způsob využití pozemku, což je údaj upřesňující. Tím jsem se teď nezabýval. Obrázky jsou na http://pedro.poloha.net/osm/landuse vždy OSM odklad a pak RUIAN obarvené parcely podle druhu pozemku. Zdá se, že to nevypadá zle, na to, že je to jen to hrubé členění. Je tam i screenshot z qgisu, kde je vysvětlení, co která barva znamená. S qgisem moc neumím a tak nevím, jak vyexportovat legendu; proto ten screenshot. Jsou tam vidět zahrady i trávníky u domů, lesy i orná půda, koryto potoka atd. S tím upřesněním by to bylo samozřejmě ještě lepší, ale to je 30 kategorií a nevím, jaké barvičky použít :-) Kombinací druh_pozemku a zpusob_vyuziti_pozemku je v RUIAN 60. Dávám to jen ke zvážení, ale teď se na to opravdu z časových důvodů nebudu vrhat, i když mě to mrzí a láká mě to. -- Petr

26.7.2014 12:26:34 (#18)
gravatar

Pavel Machek

<pavel at ucw.cz>
1056 1226
Ahoj! zobrazit citaci
> ani nic nequotuji, jen co se mi zdá, tak podle toho zákona, co sem dával odkaz > hanoj, jsou v LPIS jen ty pozemky, u kterých to jejich uživatelé chtějí. Kdo > tam nechce být, tak tam není. Chápu správně?
Jistej si nejsem, ale taky mam ten pocit. Na druhou stranu -- je to kvuli dotacim, tak tam vetsina pozemku bude. zobrazit citaci
> Udělal jsem narychlo vizualizaci RUIAN, pouze druh pozemku. V tabulkách je > ještě způsob využití pozemku, což je údaj upřesňující. Tím jsem se teď > nezabýval. > > Obrázky jsou na http://pedro.poloha.net/osm/landuse > > vždy OSM odklad a pak RUIAN obarvené parcely podle druhu pozemku. > > Zdá se, že to nevypadá zle, na to, že je to jen to hrubé členění.
Pekne, moc pekne. Kdyz to budu chtit zopakovat vlastnim qgisem -- jakou adresu musim zadat? zobrazit citaci
> Je tam i screenshot z qgisu, kde je vysvětlení, co která barva znamená. S > qgisem moc neumím a tak nevím, jak vyexportovat legendu; proto ten screenshot. > > Jsou tam vidět zahrady i trávníky u domů, lesy i orná půda, koryto potoka atd. > S tím upřesněním by to bylo samozřejmě ještě lepší, ale to je 30 kategorií a > nevím, jaké barvičky použít :-)
Kdyby se z toho podarilo dostat udaje o zahradkach a plotech, hned by mapa vypadala zas o kousek lip :-) Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

26.7.2014 03:26:49 (#19)
gravatar

Petr Vejsada

<osm at propsychology.cz>
516
Ahoj, Dne So 26. července 2014 00:26:34, Pavel Machek napsal(a): zobrazit citaci [...] zobrazit citaci
> Pekne, moc pekne. Kdyz to budu chtit zopakovat vlastnim qgisem -- > jakou adresu musim zadat?
na to potřebuješ zdroj těch RUIAN dat, třeba Postgis, .shp, ... Příslušný .qgs je na http://pedro.poloha.net/osm/ruian.qgs , musíš si ho upravit. zobrazit citaci
> Kdyby se z toho podarilo dostat udaje o zahradkach a plotech, hned by > mapa vypadala zas o kousek lip :-)
to nejsou ploty na tom obrázky, to jsou rozhraní jednotlivých polygonů. Což nevylučuje plot, ale ani neoznačuje. -- Petr

26.7.2014 09:05:20 (#20)
gravatar

Martin Švec - OSM

<osm at maatts.cz>
109
Ahoj, pro zajímavost, tvůj skript jsem při zkoumání LPISu nakonec zavrhl, jelikož jsem našel tohle: http://wiki.openstreetmap.org/wiki/Ogr2osm Stačí doplnit translation funkci (vypůjčeno z pův. skriptu, v příloze) a veškerou srandu udělá jeden příkaz: ogr2osm.py PLPIS_233406_KU_KOD_742376.shp -t trn-lpis.py -o ujezd.lpis.osm Mj. je tím vyřešeno slučování uzlů, multipolygony a nejspíš i další problémy. Martin ------------- další část --------------- A non-text attachment was scrubbed... Name: trn-lpis.py Type: text/x-python Size: 1477 bytes Desc: [žádný popis není k dispozici] URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20140726/a7f09f22/attachment.py>

28.7.2014 02:53:55 (#21)
gravatar

Pavel Machek

<pavel at ucw.cz>
1056 1226
Hi! I'd like to start import of LPIS farmland database, as we have very good coverage of houses, forests and water, but farmland is very good at places and completely missing at different places. Import page is at https://wiki.openstreetmap.org/wiki/LPIS , import script will probably be ogr2osm + script below. Best regards, Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html ------------- další část --------------- A non-text attachment was scrubbed... Name: trn-lpis.py Type: text/x-python Size: 1477 bytes Desc: [žádný popis není k dispozici] URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20140728/ea96a956/attachment.py>

28.7.2014 09:39:04 (#22)
gravatar

Martin Švec - OSM

<osm at maatts.cz>
109
Ahoj, pár věcí: * Opravil jsem skript, aby se díval i do pole KULTURA_KL, teď už taguje zelinářskou zahradu atd. * Když narazí na neznámý kód kultury, vyhodí výjmku místo tagu. * Vyházel jsem zbývající ploty -- byly tam z nějakého důvodu? * crop=vegetable nebo crop=vegetables ? taginfo zná víc těch prvních * kul. 98 (rychle rostoucí dřeviny): landuse=scrub ??? * Proč se taguje id_fb? Je to k něčemu dobrý? Zkoušel jsem nanečisto v JOSM napojování většího kusu dat na existující landuse polygony, je to pěkná drbačka, překryvů a děr je víc než dost. Pokud Marián dodělá tracer, který by přidával políčka jedno po druhým, byl by asi praktičtější. (A pohrávám si s myšlenkou dodělat do JOSM funkci, která přilepí úsek jedné cesty ke druhé mezi určenými body. Fakt nic takovýho neexistuje nebo jen špatně hledám?) Btw, rozjel jsem PostGIS a zkoušel totéž s parcelami RUIANu. Když se vhodně sloučí parcely podle typů jak psal Petr Vejsada, tak je to taky použitelné. Proti LPISu RUIAN líp "pokrývá celou plochu KÚ", ale LPIS zase obecně víc odpovídá Bingu. Problémů s napojováním na okolí je u obou zhruba stejně. (Docela pěkně vychází parcely RUIANu ve městech a vesnicích. Zkusím z toho někdy zmapovat vnitřnosti u pár vesnic, kolik klikací práce to ušetří.) Martin On 28.7.2014 14:53, Pavel Machek wrote: zobrazit citaci
> Hi! > > I'd like to start import of LPIS farmland database, as we have very > good coverage of houses, forests and water, but farmland is very good > at places and completely missing at different places. > > Import page is at https://wiki.openstreetmap.org/wiki/LPIS , import > script will probably be ogr2osm + script below. > > Best regards, > Pavel > > > > > _______________________________________________ > Talk-cz mailing list > Talk-cz na openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz
------------- další část --------------- A non-text attachment was scrubbed... Name: trn-lpis.py Type: text/x-python Size: 1584 bytes Desc: [žádný popis není k dispozici] URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20140728/9d5b09a7/attachment.py>

28.7.2014 10:35:11 (#23)
gravatar

Pavel Machek

<pavel at ucw.cz>
1056 1226
Ahoj! zobrazit citaci
> pro zajímavost, tvůj skript jsem při zkoumání LPISu nakonec zavrhl, jelikož jsem našel tohle: > > http://wiki.openstreetmap.org/wiki/Ogr2osm > > Stačí doplnit translation funkci (vypůjčeno z pův. skriptu, v příloze) a veškerou srandu udělá jeden příkaz: > > ogr2osm.py PLPIS_233406_KU_KOD_742376.shp -t trn-lpis.py -o ujezd.lpis.osm
No, stahnul jsem si ogr2osm.py, ale nechce fungovat: git clone git://github.com/pnorman/ogr2osm.git pavel na amd:~/g/ogr2osm$ ./ogr2osm.py File "./ogr2osm.py", line 538 featuresmap = {feature.geometry : feature for feature in Feature.features} ^ SyntaxError: invalid syntax pavel na amd:~/g/ogr2osm$ python3 ./ogr2osm.py File "./ogr2osm.py", line 320 reproject = lambda(geometry): None ^ SyntaxError: invalid syntax pavel na amd:~/g/ogr2osm$ To druhy jsem zvladnul nejak osalit, jenze nemam potrebny balicky pro python3: from osgeo import ogr from osgeo import osr . Jestli ma nekdo napad... Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

28.7.2014 10:42:46 (#24)
gravatar

Pavel Machek

<pavel at ucw.cz>
1056 1226
On Mon 2014-07-28 22:35:11, Pavel Machek wrote: zobrazit citaci
> Ahoj! > > > pro zajímavost, tvůj skript jsem při zkoumání LPISu nakonec zavrhl, jelikož jsem našel tohle: > > > > http://wiki.openstreetmap.org/wiki/Ogr2osm > > > > Stačí doplnit translation funkci (vypůjčeno z pův. skriptu, v příloze) a veškerou srandu udělá jeden příkaz: > > > > ogr2osm.py PLPIS_233406_KU_KOD_742376.shp -t trn-lpis.py -o ujezd.lpis.osm > > No, stahnul jsem si ogr2osm.py, ale nechce fungovat: > > git clone git://github.com/pnorman/ogr2osm.git > > pavel na amd:~/g/ogr2osm$ ./ogr2osm.py > File "./ogr2osm.py", line 538 > featuresmap = {feature.geometry : feature for feature in > Feature.features} > ^ > SyntaxError: invalid syntax > pavel na amd:~/g/ogr2osm$ python3 ./ogr2osm.py > File "./ogr2osm.py", line 320 > reproject = lambda(geometry): None > ^ > SyntaxError: invalid syntax > pavel na amd:~/g/ogr2osm$ > > To druhy jsem zvladnul nejak osalit, jenze nemam potrebny balicky pro python3:
Tak vyreseno, ted by to melo bezet pod python 2.6 i 3+ Pavel diff --git a/ogr2osm.py b/ogr2osm.py index 587a5ae..0184262 100755 --- a/ogr2osm.py +++ b/ogr2osm.py @@ -317,13 +317,13 @@ def getTransform(layer): # Some python magic: skip reprojection altogether by using a dummy # lamdba funcion. Otherwise, the lambda will be a call to the OGR # reprojection stuff. - reproject = lambda(geometry): None + reproject = lambda geometry: None else: destSpatialRef = osr.SpatialReference() # Destionation projection will *always* be EPSG:4326, WGS84 lat-lon destSpatialRef.ImportFromEPSG(4326) coordTrans = osr.CoordinateTransformation(spatialRef, destSpatialRef) - reproject = lambda(geometry): geometry.Transform(coordTrans) + reproject = lambda geometry: geometry.Transform(coordTrans) return reproject @@ -535,7 +535,7 @@ def output(): nodes = [geom for geom in Geometry.geometries if type(geom) == Point] ways = [geom for geom in Geometry.geometries if type(geom) == Way] relations = [geom for geom in Geometry.geometries if type(geom) == Relation] - featuresmap = {feature.geometry : feature for feature in Feature.features} + featuresmap = dict([(feature.geometry, feature) for feature in Feature.features]) # Open up the output file with the system default buffering with open(options.outputFile, 'w', -1) as f: -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

28.7.2014 10:53:53 (#25)
gravatar

Martin Švec - OSM

<osm at maatts.cz>
109
On 28.7.2014 22:35, Pavel Machek wrote: zobrazit citaci
> To druhy jsem zvladnul nejak osalit, jenze nemam potrebny balicky pro python3: > > from osgeo import ogr > from osgeo import osr > > . Jestli ma nekdo napad... > Pavel
Mně to fungovalo samo od sebe... Koukám do systému, ogr.py/osr.py mi do Gentoo zavlekla GDAL knihovna s python USE flagem. A jako aktivní python mám 2.7. Martin

28.7.2014 10:59:51 (#26)
gravatar

Marián Kyral

<mkyral at email.cz>
2389 2735
Dne 28.7.2014 21:39, Martin S(vec - OSM napsal(a): zobrazit citaci
> Ahoj, pár ve(cí: > > * Opravil jsem skript, aby se díval i do pole KULTURA_KL, ted( uz( > taguje zelinár(skou zahradu atd. > * Kdyz( narazí na neznámý kód kultury, vyhodí výjmku místo tagu. > * Vyházel jsem zbývající ploty -- byly tam z ne(jakého du*vodu? > * crop=vegetable nebo crop=vegetables ? taginfo zná víc te(ch prvních > * kul. 98 (rychle rostoucí dr(eviny): landuse=scrub ??? > * Proc( se taguje id_fb? Je to k ne(c(emu dobrý? > > Zkous(el jsem nanec(isto v JOSM napojování ve(ts(ího kusu dat na > existující landuse polygony, je to pe(kná drbac(ka, pr(ekryvu* a de(r > je víc nez( dost. Pokud Marián dode(lá tracer, který by pr(idával > políc(ka jedno po druhým, byl by asi praktic(te(js(í.
Dode(lám, dode(lám. Sám to uz( nutne( potr(ebuji ;-) Jen nevím kdy to bude. Momentálne( r(es(ím konflikty v poslední verzi Traceru a zatím se mi ne(jak nedar(í. Je tam ne(co s(patne(, ale ne a ne pr(ijít na to co :-( zobrazit citaci
> (A pohrávám si s mys(lenkou dode(lat do JOSM funkci, která pr(ilepí > úsek jedné cesty ke druhé mezi urc(enými body. Fakt nic takovýho > neexistuje nebo jen s(patne( hledám?) >
Moz(ná je trochu problém, jak poznat co se k c(emu má pr(ilepit, pr(ípadne( jak to jednodus(e zadat. Marián zobrazit citaci
> Btw, rozjel jsem PostGIS a zkous(el totéz( s parcelami RUIANu. Kdyz( > se vhodne( slouc(í parcely podle typu* jak psal Petr Vejsada, tak je > to taky pouz(itelné. Proti LPISu RUIAN líp "pokrývá celou plochu KÚ", > ale LPIS zase obecne( víc odpovídá Bingu. Problému* s napojováním na > okolí je u obou zhruba stejne(. (Docela pe(kne( vychází parcely RUIANu > ve me(stech a vesnicích. Zkusím z toho ne(kdy zmapovat vnitr(nosti u > pár vesnic, kolik klikací práce to us(etr(í.) > > Martin > > On 28.7.2014 14:53, Pavel Machek wrote: >> Hi! >> >> I'd like to start import of LPIS farmland database, as we have very >> good coverage of houses, forests and water, but farmland is very good >> at places and completely missing at different places. >> >> Import page is at https://wiki.openstreetmap.org/wiki/LPIS , import >> script will probably be ogr2osm + script below. >> >> Best regards, >> Pavel >> >> >> >> >> _______________________________________________ >> 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
------------- dal?í ?ást --------------- HTML p?íloha byla odstran?na... URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20140728/32669bc0/attachment.html>

28.7.2014 11:03:40 (#27)
gravatar

Pavel Machek

<pavel at ucw.cz>
1056 1226
Ahoj! zobrazit citaci
> * Opravil jsem skript, aby se díval i do pole KULTURA_KL, teď už taguje zelinářskou zahradu atd. > * Když narazí na neznámý kód kultury, vyhodí výjmku místo tagu. > * Vyházel jsem zbývající ploty -- byly tam z nějakého důvodu?
No, pro me je informace o plotech docela dulezita.. mam pocit ze vinice a chmelnice plot proste maji. Ale asi je to trochu drze... zobrazit citaci
> * crop=vegetable nebo crop=vegetables ? taginfo zná víc těch prvních
Tak vemem crop=vegerable... zobrazit citaci
> * kul. 98 (rychle rostoucí dřeviny): landuse=scrub ???
http://wiki.openstreetmap.org/wiki/Tag:natural%3Dscrub Aha, ma tam byt natural, ale jinak to myslim sedi. Opravim. zobrazit citaci
> * Proč se taguje id_fb? Je to k něčemu dobrý?
Asi neni. Vyhodim ho. zobrazit citaci
> Zkoušel jsem nanečisto v JOSM napojování většího kusu dat na > existující landuse polygony, je to pěkná drbačka, překryvů a děr je > víc než dost. Pokud Marián dodělá tracer, který by přidával políčka > jedno po druhým, byl by asi praktičtější.
No, ja bych uplne landuse=farmland napriklad s lesama nespojoval.. ono vedle toho lesa byva ruzne siroky pruh ne-lesa ale jeste taky ne-pole. zobrazit citaci
> (A pohrávám si s myšlenkou > dodělat do JOSM funkci, která přilepí úsek jedné cesty ke druhé mezi > určenými body. Fakt nic takovýho neexistuje nebo jen špatně hledám?)
Nevim o ni. zobrazit citaci
> Btw, rozjel jsem PostGIS a zkoušel totéž s parcelami RUIANu. Když se > vhodně sloučí parcely podle typů jak psal Petr Vejsada, tak je to > taky použitelné. Proti LPISu RUIAN líp "pokrývá celou plochu KÚ", > ale LPIS zase obecně víc odpovídá Bingu. Problémů s napojováním na > okolí je u obou zhruba stejně. (Docela pěkně vychází parcely RUIANu > ve městech a vesnicích. Zkusím z toho někdy zmapovat vnitřnosti u > pár vesnic, kolik klikací práce to ušetří.)
No, kdyby z toho byly nejaky soubory, treba pro okoli rican, tak se na rad podivam. Kdyby se v osm objevily ploty, tak me to hodne potesi... a ploty by mely jit vykoukat z mistni znalosti + hranic pozemku. Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

29.7.2014 12:28:46 (#28)
gravatar

Martin Švec - OSM

<osm at maatts.cz>
109
On 28.7.2014 23:03, Pavel Machek wrote: zobrazit citaci
> No, pro me je informace o plotech docela dulezita.. mam pocit ze > vinice a chmelnice plot proste maji. Ale asi je to trochu drze...
Mně se nelíbily ploty kolem a napříč loukami, kde vím že tam určitě nejsou. U vinic+chmelnic... mno, tam je pravděpodobnost plotu asi docela vysoká. Těžko říct... zobrazit citaci
> http://wiki.openstreetmap.org/wiki/Tag:natural%3Dscrub > > Aha, ma tam byt natural, ale jinak to myslim sedi. Opravim.
Znova jsem se nad tím zamyslel... Záleží, co se myslí "porostem rychle rostoucích dřevin". Natural=scrub je neobhospodařovaná příroda, landuse=* je obhospodařovaná krajina. Jestli jsou to plantáže na biomasu, landuse by tomu slušel víc. (Taginfo landuse=scrub zná, ale většinou na uzlech (wtf, jednotlivé keře?). Co jsem namátkou našel na polygonech, tak je to chybně otagovaný natural=scrub.) Update: jo, je to biomasa: viz par. 3i písm. j) zákona č. 252/1997 Sb. zobrazit citaci
>> Zkoušel jsem nanečisto v JOSM napojování většího kusu dat na >> existující landuse polygony, je to pěkná drbačka, překryvů a děr je >> víc než dost. Pokud Marián dodělá tracer, který by přidával políčka >> jedno po druhým, byl by asi praktičtější. > No, ja bych uplne landuse=farmland napriklad s lesama nespojoval.. ono > vedle toho lesa byva ruzne siroky pruh ne-lesa ale jeste taky ne-pole.
To je asi individuální... Když jsou mezi lesem a polem tři metry křovin a šouší, dávám pruh natural=scrub. Když jsou na okraji jen koleje od traktoru a rovnou les, spojuju přímo. Nemám rád v mapě bílé škvíry :-) Ty pruhy křoví jsou občas přímo v KM. zobrazit citaci
>> Btw, rozjel jsem PostGIS a zkoušel totéž s parcelami RUIANu. Když se >> vhodně sloučí parcely podle typů jak psal Petr Vejsada, tak je to >> taky použitelné. Proti LPISu RUIAN líp "pokrývá celou plochu KÚ", >> ale LPIS zase obecně víc odpovídá Bingu. Problémů s napojováním na >> okolí je u obou zhruba stejně. (Docela pěkně vychází parcely RUIANu >> ve městech a vesnicích. Zkusím z toho někdy zmapovat vnitřnosti u >> pár vesnic, kolik klikací práce to ušetří.) > No, kdyby z toho byly nejaky soubory, treba pro okoli rican, tak se na > rad podivam. Kdyby se v osm objevily ploty, tak me to hodne > potesi... a ploty by mely jit vykoukat z mistni znalosti + hranic pozemku.
Pokusím se vyrobit ukázku, až to trochu učešu. Ale teď se k tomu pár dnů nedostanu. Zatím hlavně googlím, co vůbec PostGIS nad daty dokáže. Btw, ty ploty... To je nějaká specifická úchylka? :-)) Martin

29.7.2014 12:58:49 (#29)
gravatar

Pavel Machek

<pavel at ucw.cz>
1056 1226
Ahoj! zobrazit citaci
> >No, pro me je informace o plotech docela dulezita.. mam pocit ze > >vinice a chmelnice plot proste maji. Ale asi je to trochu drze... > > Mně se nelíbily ploty kolem a napříč loukami, kde vím že tam určitě nejsou. U vinic+chmelnic... mno, tam je pravděpodobnost plotu asi docela vysoká. Těžko říct... >
Jo, tak ty uz jsem zrusil. Mel jsem pocit, ze kolem "stale pastviny" by mel byt plot, ale realita se mnou nesouhlasila. zobrazit citaci
> >http://wiki.openstreetmap.org/wiki/Tag:natural%3Dscrub > > > >Aha, ma tam byt natural, ale jinak to myslim sedi. Opravim. > > Znova jsem se nad tím zamyslel... Záleží, co se myslí "porostem rychle rostoucích dřevin". Natural=scrub je neobhospodařovaná příroda, landuse=* je obhospodařovaná krajina. Jestli jsou to plantáže na biomasu, landuse by tomu slušel víc. (Taginfo landuse=scrub zná, ale většinou na uzlech (wtf, jednotlivé keře?). Co jsem namátkou našel na polygonech, tak je to chybně otagovaný natural=scrub.) >
landuse=scrub neni problem. zobrazit citaci
> >>Btw, rozjel jsem PostGIS a zkoušel totéž s parcelami RUIANu. Když se > >>vhodně sloučí parcely podle typů jak psal Petr Vejsada, tak je to > >>taky použitelné. Proti LPISu RUIAN líp "pokrývá celou plochu KÚ", > >>ale LPIS zase obecně víc odpovídá Bingu. Problémů s napojováním na > >>okolí je u obou zhruba stejně. (Docela pěkně vychází parcely RUIANu > >>ve městech a vesnicích. Zkusím z toho někdy zmapovat vnitřnosti u > >>pár vesnic, kolik klikací práce to ušetří.) > >No, kdyby z toho byly nejaky soubory, treba pro okoli rican, tak se na > >rad podivam. Kdyby se v osm objevily ploty, tak me to hodne > >potesi... a ploty by mely jit vykoukat z mistni znalosti + hranic pozemku. > > Pokusím se vyrobit ukázku, až to trochu učešu. Ale teď se k tomu pár dnů nedostanu. Zatím hlavně googlím, co vůbec PostGIS nad daty dokáže. Btw, ty ploty... To je nějaká specifická úchylka? :-)) >
No... z mapy me v zasade zajima jedina vec. "Da se tudy projit?". No, a je jenom nekolik malo veci, co se projit nedaji: Budova, plot a repkovy pole. Budovy uz v mape mape, s polema se snazim pohnout, no tak zbyvaj ty ploty :-). Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

29.7.2014 01:06:19 (#30)
gravatar

Dalibor Jelínek

<dalibor at dalibor.cz>
415 1552
Ahoj, co myslis tim, ze landuse=scrub neni problem? Ja bych pouziti tohoto tagu jako problem videl, protoze nam to nic nevykresli a navic si lidi budou myslet, ze je to preklep z natural=scrub. Skutecne se jedna o nejake plantaze biomasy? Pokud jo, pak by tomu asi slusel jiny tag a to asi landuse=farmland + nejaky crop= Ale spis bych si tipoval, ze to bude proste jen krovi, ne? Zdravi, Dalibor zobrazit citaci
> -----Original Message----- > From: Pavel Machek [mailto:pavel na ucw.cz]
zobrazit citaci
> > > >http://wiki.openstreetmap.org/wiki/Tag:natural%3Dscrub > > > > > >Aha, ma tam byt natural, ale jinak to myslim sedi. Opravim. > > > > Znova jsem se nad tím zamyslel... Záleží, co se myslí "porostem rychle > > rostoucích dřevin". Natural=scrub je neobhospodařovaná příroda, > > landuse=* je obhospodařovaná krajina. Jestli jsou to plantáže na > > biomasu, landuse by tomu slušel víc. (Taginfo landuse=scrub zná, ale > > většinou na uzlech (wtf, jednotlivé keře?). Co jsem namátkou našel na > > polygonech, tak je to chybně otagovaný natural=scrub.) > > > > landuse=scrub neni problem. >

30.7.2014 12:28:26 (#31)
gravatar

Pavel Machek

<pavel at ucw.cz>
1056 1226
Ahoj! zobrazit citaci
> co myslis tim, ze landuse=scrub neni problem? > Ja bych pouziti tohoto tagu jako problem videl, > protoze nam to nic nevykresli a navic si lidi budou > myslet, ze je to preklep z natural=scrub.
Myslel jsem tim, ze je pro me jednoduchy to zmenit. zobrazit citaci
> Skutecne se jedna o nejake plantaze biomasy? > Pokud jo, pak by tomu asi slusel jiny tag a to > asi landuse=farmland + nejaky crop=
Jsou to "rychle rostouci dreviny". zobrazit citaci
> Ale spis bych si tipoval, ze to bude proste jen krovi, ne?
Ja bych to spis videl na ty plantaze... Asi by slo "landuse=forest, forest=scrub", ale mozna bude nejjednodussi nechat tam "natural=scrub" i kdyz to neni uplne presny... Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

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