Re: [Talk-cz] Import of farmland from LPIS
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 -Original Message- From: Pavel Machek [mailto:pa...@ucw.cz] 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. ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
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. Myslel jsem tim, ze je pro me jednoduchy to zmenit. 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. 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 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
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: 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@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz def getKultura(attrs): if 'KULTURA_KL' in attrs and attrs['KULTURA_KL'] != '': return int(attrs['KULTURA_KL']) if 'KULTURA' in attrs and attrs['KULTURA'] != '': kul = int(attrs['KULTURA']) return -1; def filterTags(attrs): if not attrs: return tags = {} tags['source'] = 'lpis' if 'ID_FB' in attrs: tags['id_fb'] = attrs['ID_FB']; if 'VYSKA' in attrs: tags[ele] = %d % int(float(attrs['VYSKA'])) if 'NKODFB' in attrs: tags[ref] = attrs['NKODFB'] kul = getKultura(attrs) if kul == 2:tags[landuse] = farmland elif kul == 3: tags[landuse] = farmland; tags[crop] = hop elif kul == 30: tags[landuse] = farmland; tags[crop] = hop elif kul == 31: tags[landuse] = farmland; tags[crop] = hop elif kul == 41: tags[landuse] = vineyard elif kul == 61: tags[landuse] = orchard elif kul == 62: tags[landuse] = orchard elif kul == 7: tags[landuse] = meadow; tags[meadow] = agricultural elif kul == 71: tags[landuse] = meadow; tags[meadow] = agricultural elif kul == 72: tags[landuse] = meadow; tags[meadow] = agricultural elif kul == 91: tags[landuse] = forest elif kul == 92: tags[landuse] = farm; tags[crop] = vegetables elif kul == 97: tags[landuse] = reservoir elif kul == 98: tags[landuse] = scrub; tags[note]=rychle rostouci dreviny elif kul == 99: tags[landuse] =forest else: raise Exception (unknown farmland: %s % kul) return tags ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
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@amd:~/g/ogr2osm$ ./ogr2osm.py File ./ogr2osm.py, line 538 featuresmap = {feature.geometry : feature for feature in Feature.features} ^ SyntaxError: invalid syntax pavel@amd:~/g/ogr2osm$ python3 ./ogr2osm.py File ./ogr2osm.py, line 320 reproject = lambda(geometry): None ^ SyntaxError: invalid syntax pavel@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 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
On Mon 2014-07-28 22:35:11, Pavel Machek wrote: 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@amd:~/g/ogr2osm$ ./ogr2osm.py File ./ogr2osm.py, line 538 featuresmap = {feature.geometry : feature for feature in Feature.features} ^ SyntaxError: invalid syntax pavel@amd:~/g/ogr2osm$ python3 ./ogr2osm.py File ./ogr2osm.py, line 320 reproject = lambda(geometry): None ^ SyntaxError: invalid syntax pavel@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 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
On 28.7.2014 22:35, Pavel Machek wrote: 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 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
Dne 28.7.2014 21:39, Martin S(vec - OSM napsal(a): 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 :-( (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 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@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
Ahoj! * 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... * crop=vegetable nebo crop=vegetables ? taginfo zná víc těch prvních Tak vemem crop=vegerable... * 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. * Proč se taguje id_fb? Je to k něčemu dobrý? Asi neni. Vyhodim ho. 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. (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. 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 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
On 28.7.2014 23:03, Pavel Machek wrote: 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... 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. 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. 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 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
Ahoj, Dne So 26. července 2014 00:26:34, Pavel Machek napsal(a): Obrázky jsou na http://pedro.poloha.net/osm/landuse [...] 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. 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 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
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 def filterTags(attrs): if not attrs: return tags = {} tags['source'] = 'lpis' if 'ID_FB' in attrs: tags['id_fb'] = attrs['ID_FB']; if 'VYSKA' in attrs: tags[ele] = %d % int(float(attrs['VYSKA'])) if 'NKODFB' in attrs: tags[ref] = attrs['NKODFB'] kul = -1 if 'KULTURA' in attrs and attrs['KULTURA'] != '': kul = int(attrs['KULTURA']) if kul == 2:tags[landuse] = farmland elif kul == 3: tags[landuse] = farmland; tags[crop] = hop elif kul == 30: tags[landuse] = farmland; tags[crop] = hop elif kul == 31: tags[landuse] = farmland; tags[crop] = hop elif kul == 41: tags[landuse] = vineyard; tags[barrier] = fence elif kul == 61: tags[landuse] = orchard; tags[barrier] = fence elif kul == 62: tags[landuse] = orchard elif kul == 7: tags[landuse] = meadow; tags[meadow] = agricultural elif kul == 71: tags[landuse] = meadow; tags[meadow] = agricultural elif kul == 72: tags[landuse] = meadow; tags[meadow] = agricultural elif kul == 91: tags[landuse] = forest; elif kul == 92: tags[landuse] = farm; tags[crop] = vegetables elif kul == 97: tags[landuse] = reservoir elif kul == 98: tags[landuse] = scrub; tags[note]=rychle rostouci dreviny elif kul == 99: tags[landuse] =forest else: tags[landuse] =unknown_farmland_%s % kul return tags ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
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ě? Jistej si nejsem, ale taky mam ten pocit. Na druhou stranu -- je to kvuli dotacim, tak tam vetsina pozemku bude. 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? 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 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
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 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
Dne 23.7.2014 2:43, Petr Vejsada napsal(a): 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 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
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í). Diky za skript! * 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. * Spousta sousedících polygonů má překryvy, i v rámci jednoho KU. Jako ze LPIS se prekryva sam se sebou? * 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. * 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. * 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? * 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? Na druhou stranu, už teď skript ušetří mraky ručního obkreslování, stačí jen v JOSM doopravit problémy. :-) Presne. (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 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
Ahoj! 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. 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. 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. 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. 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 helpd...@mze.cz ... 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. [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 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
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 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
[Talk-cz] Import of farmland from LPIS
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 #!/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' ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
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. 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 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
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 pa...@ucw.cz: 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@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
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): Import page is at https://wiki.openstreetmap.org/wiki/LPIS , import script is attached to this email. ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
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? 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@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
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? Tvrdi, ze jim to zatezuje procesor... - 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. - 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. - 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... 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. 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. 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 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
Ahoj, Dne Út 22. července 2014 23:07:38, Pavel Machek napsal(a): - 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ě. 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? ;-) - 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? Krajinne prvky jsem zatim nezkoumal, me zajimaly hlavne louky. To bych nechal jako separatni problem. Jo, to můžeme. 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. Aktualizace.. no vzhledem k captche bych ocekaval ze se udelaj rucne kdyz se dana oblast zmeni. ? nebo ten strojový čas/zátěž, ne? 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? 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- ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
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: 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@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz #!/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 = 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
Re: [Talk-cz] Import of farmland from LPIS
Ahoj, Dne St 23. července 2014 02:07:16, Martin Švec - OSM napsal(a): * 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- ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz