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>'

Hi! zobrazit citaci 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

zobrazit citaci *** 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

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

On Tue 2014-07-22 15:40:32, hanoj wrote: zobrazit citaci 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

Ahoj! zobrazit citaci Tvrdi, ze jim to zatezuje procesor... zobrazit citaci 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 No, konverzni skript mam, nekde je pribalenej. Ziskat ty zip fily je obtizny, zbytek uz je docela dobrej. zobrazit citaci 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 Krajinne prvky jsem zatim nezkoumal, me zajimaly hlavne louky. To bych nechal jako separatni problem. zobrazit citaci 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 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

Ahoj, Dne Út 22. července 2014 23:07:38, Pavel Machek napsal(a): 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. 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 nj, rozumím, docela v poho, když děláš okolí Říčan. Kolikrát se Říčany vejdou do ČR? ;-) 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, 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 Jo, to můžeme. 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. zobrazit citaci ? nebo ten strojový čas/zátěž, ne? zobrazit citaci 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 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-

Ahoj! zobrazit citaci Dobra prace :-). 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.. zobrazit citaci 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 Aha, ok. zobrazit citaci 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 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 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 Les je blbe, delala se tam moc silna generalizace. zobrazit citaci (Jinak k tomu ekologickymu zemedelstvi jsem nenasel odpovidajici tag na osm, takze to prozatim ignoruju). zobrazit citaci 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 Ten pyproj tam nakonec neni potreba. shapefile se instaluje pomoci pip-u. A sorry, Windows neznam a nedoporucuju. zobrazit citaci Dat na ftp? Pro zacatek by se hodilo nezkonvertovany, protoze skript bude jeste potrebovat par uprav. zobrazit citaci 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

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 ------------- 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>'

Ahoj, Dne St 23. července 2014 01:17:02, Pavel Machek napsal(a): zobrazit citaci "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 Aha, kruh 60km v průměru, no odsejpá ti to. zobrazit citaci No, myslím, že buďto je les nebo je pastvina. Obojímu dohromady se říká step ;-). zobrazit citaci No OSM je obecná mapa, tak tam nemusí být úplně všechno zobrazit citaci 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 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 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 ... a bude landuse=farmland over landuse=farmland :-( -- -p-

Ahoj, Dne St 23. července 2014 02:07:16, Martin Švec - OSM napsal(a): zobrazit citaci 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-

zobrazit citaci *** 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

Dne 23.7.2014 2:43, Petr Vejsada napsal(a): zobrazit citaci 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

Ahoj! zobrazit citaci Diky za skript! zobrazit citaci 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 Jako ze LPIS se prekryva sam se sebou? zobrazit citaci To se asi ocekava, ne? Lesy a vodu uz mame, tak ted pole a louky. zobrazit citaci Tak, ale ted jsem to vypnul... ono to moc nesedi. zobrazit citaci Mozna ze to pole neni "do kraje" parcely? zobrazit citaci Jsou v RUIANu udaje co kde je? Je nekde konverzni skript/example vystup? zobrazit citaci Presne. zobrazit citaci 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

On Wed 2014-07-23 09:29:54, hanoj wrote: zobrazit citaci 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

Ahoj! zobrazit citaci 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 No, nahrazovani existujicich dat by se asi melo delat rucne -- z ucty k tem co ty data vyrobili. zobrazit citaci 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 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 No, uvadej tam helpdesk na mze.cz ... zobrazit citaci Aby lidi co vypnej formular nemuseli jeste zarizovat emailovou adresu s rozbalovanim. ... ok, asi by bylo jednodussi rozbalovat to z mbox formatu. zobrazit citaci 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

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

Ahoj! zobrazit citaci Jistej si nejsem, ale taky mam ten pocit. Na druhou stranu -- je to kvuli dotacim, tak tam vetsina pozemku bude. zobrazit citaci Pekne, moc pekne. Kdyz to budu chtit zopakovat vlastnim qgisem -- jakou adresu musim zadat? zobrazit citaci 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

Ahoj, Dne So 26. července 2014 00:26:34, Pavel Machek napsal(a): zobrazit citaci [...] zobrazit citaci 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 to nejsou ploty na tom obrázky, to jsou rozhraní jednotlivých polygonů. Což nevylučuje plot, ale ani neoznačuje. -- Petr

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>

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>

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 ------------- 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>

Ahoj! zobrazit citaci 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

On Mon 2014-07-28 22:35:11, Pavel Machek wrote: zobrazit citaci 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

On 28.7.2014 22:35, Pavel Machek wrote: zobrazit citaci 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

Dne 28.7.2014 21:39, Martin S(vec - OSM napsal(a): zobrazit citaci 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 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 ------------- dal?í ?ást --------------- HTML p?íloha byla odstran?na... URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20140728/32669bc0/attachment.html>

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... zobrazit citaci Tak vemem crop=vegerable... zobrazit citaci http://wiki.openstreetmap.org/wiki/Tag:natural%3Dscrub Aha, ma tam byt natural, ale jinak to myslim sedi. Opravim. zobrazit citaci Asi neni. Vyhodim ho. zobrazit citaci 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 Nevim o ni. zobrazit citaci 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

On 28.7.2014 23:03, Pavel Machek wrote: zobrazit citaci 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 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 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 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

Ahoj! zobrazit citaci Jo, tak ty uz jsem zrusil. Mel jsem pocit, ze kolem "stale pastviny" by mel byt plot, ale realita se mnou nesouhlasila. zobrazit citaci landuse=scrub neni problem. zobrazit citaci 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

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 zobrazit citaci

Ahoj! zobrazit citaci Myslel jsem tim, ze je pro me jednoduchy to zmenit. zobrazit citaci Jsou to "rychle rostouci dreviny". zobrazit citaci 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