Re: [Talk-cz] LPIS import
Ahoj, Dne Po 21. července 2014 01:00:19, Pavel Machek napsal(a): > Maji tam klic podle kteryho to pujde aktualizovat, a normalne > zemedelska puda relativne trvanliva je. Spojit s jinym landusem -- > spis radej ne, ne? Maji to docela podrobny, tj pole zacina na kraji cesty. no to asi jo, kde je pole, tam bude pole, dokud nějaký developer neuplatí úředníky a nepostaví tam mrakodrapy. Myslel jsem spíš proměny pěstovaných plodin,ale to asi do OSM dávat nebudeme. Prostě pole je pole a hotovo. Spojovat s jiným landuse - no vlastně to, co navrhuješ, t.j. jen zemědělskou půdu, to není špatné řešení.. Nač si komplikovat život. Asi už je opravdu čas, abych udělal TMS vrstvu landuse z dat z RUIAN. Jen tak si to vizualizovat a teprve uvidíme podle obrázků, zda by to k něčemu mohlo být dobré. > > Určitě najdeme spoustu míst, kde podle RUIAN je sídliště a podle EAGRI > > louka či pole (po zkušenostech s importem adres ;-)). > > No, tohle bude krasne videt na bingu :-). Zná vůbec někdo zpoždění Bingu? T.j. zda lze zjistit datum snímků? Možná je to triviální, ale teď nevím. > > > Nakolik budeme která data brát jako referenční či jak se rozhodneme, co > > převezmeme. Možná existují i další zdroje pro landuse.? > > > > Co budeme dělat se stávajícími polygony landuse. Dělat díry v polygonech z > > EAGRI? > > No, asi zalezi co bude lepsi. V oblastech, co jsem mapoval ja, je LPIS > lepsi nez rucni trasovani z ortofota... Takze pekne vymazat a nahrat > znova. No tady narazíme na to, že ne vše je zemědělská půda. Myslíš vymazat landuse=residential či industrial? Asi ne. -- Petr ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] LPIS import
Ahoj, Dne Po 21. července 2014 01:10:53, Pavel Machek napsal(a): > > A ted.. umel by nekdo automaticky z mbox souboru s mailama vytahnout > > zip-ovy prilohy? ano, ripmime, celé roky jsem si takhle posílal mailem data z meteostanice na zahradě. Na server přišel od meteostanice mail s přílohou, server předhodil mail ripmime, ten to rozbalil a nahrál do Postgresu. Není třeba to tahat z mboxu, ale mailer (v mém případě sendmail/procmail, nesmějte se, no) to rovnou předhodí ripmime a rovnou dál, tedy zpracování okamžitě. > Dalsi vec co by se hodila.. vytahnout nejakym xapi seznam vsech ref od > > boundary=administrative, admin_level=10 . To se potom zadava do jejich > webu... Originál je v RUIAN, seznam katastrálních území, je jich něco přes 13 tisíc. Je to totéž, co admin-level 10. Není problém. > Dobrou noc, Brou! -- Petr ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] LPIS import
Ahoj! > A ted.. umel by nekdo automaticky z mbox souboru s mailama vytahnout > zip-ovy prilohy? Dalsi vec co by se hodila.. vytahnout nejakym xapi seznam vsech ref od boundary=administrative, admin_level=10 . To se potom zadava do jejich webu... Ja bych na to asi casem prisel, ale pro dnesek bylo dost programovani... Dobrou noc, 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] LPIS import
Ahoj!
> > A ted ten skript.
Tak dobra zprava, ten web to umi dat we wgs-84.
Spatna je, ze to posila mailem, dle cisla katastralniho uzemi. Aha, a
ty cisla katastralnich uzemi jde dokonce vykoukat z osm, jako ref u
boundary=administrative, admin_level=10. Bohuzel jsou to maly uzemi.
> :-), normálka ;). Landuse je opravdu hodně potřeba. Vždycky tak trochu
> "závidím" karviňákům tu jejich barevnou mapu. Přemýšlel jsem o možnosti
> využít
> parcely z RUIAN. Tam je nejen zemědělská půda, ale i leccos ostatního -
> zastavěná plocha, silnice, dálnice a co já vím, co ještě. Zemědělská půda tam
> určitě není tak podrobně jako v eagri/lpis.
> Máš představu, nakolik jsou údaje trvanlivé, t.j. zda příští rok ještě budou
> platit? Jak aktualizovat? Jak to spojit s jiným landuse, neb nejen
> zemědělskou
> půdou povrch naší matičky pokryt jest.
Maji tam klic podle kteryho to pujde aktualizovat, a normalne
zemedelska puda relativne trvanliva je. Spojit s jinym landusem --
spis radej ne, ne? Maji to docela podrobny, tj pole zacina na kraji cesty.
> Dle právě provedeného výpočtu má v RUIAN digitální mapu 64 procent území
> republiky, takže sice dobré, ale ne dost.
>
> Určitě najdeme spoustu míst, kde podle RUIAN je sídliště a podle EAGRI louka
> či pole (po zkušenostech s importem adres ;-)).
No, tohle bude krasne videt na bingu :-).
> Nakolik budeme která data brát jako referenční či jak se rozhodneme, co
> převezmeme. Možná existují i další zdroje pro landuse.?
>
> Co budeme dělat se stávajícími polygony landuse. Dělat díry v polygonech z
> EAGRI?
No, asi zalezi co bude lepsi. V oblastech, co jsem mapoval ja, je LPIS
lepsi nez rucni trasovani z ortofota... Takze pekne vymazat a nahrat
znova.
A ted.. umel by nekdo automaticky z mbox souboru s mailama vytahnout
zip-ovy prilohy?
Pavel
(Aktualni verse importeru v priloze).
--
(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 ""
print ""
#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 = ''
tags += ''
nodeid -= 1
print '%s' % ( 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 '' % nodeid
print ' '
print ' ' % attr(shrec, 'ID_FB')
print ' '
print ' ' % attr(shrec, 'KULTURA')
kul = attr(shrec, 'KULTURA')
if kul == 2:print ' '
elif kul == 3: print ' '
elif kul == 30: print ' '
elif kul == 31: print ' '
elif kul == 41: print ' '
elif kul == 61: print ' '
elif kul == 62: print ' '
elif kul == 7: print ' '
elif kul == 71: print ' '
elif kul == 72: print ' '
elif kul == 91: print ' '
elif kul == 92: print ' '
elif kul == 97: print ' '
elif kul == 98: print ' '
elif kul == 99: print ' '
else: print ' ' % kul
print ' ' % attr(shrec, 'VYSKA')
print ' ' % attr(shrec, 'NKODFB')
for pt in pts:
print ' ' % pt
print ''
print ''
___
Talk-cz mailing list
Talk-cz@openstreetmap.org
https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] LPIS import
Ahoj, Dne Ne 20. července 2014 23:53:53, Pavel Machek napsal(a): > A ted ten skript. :-), normálka ;). Landuse je opravdu hodně potřeba. Vždycky tak trochu "závidím" karviňákům tu jejich barevnou mapu. Přemýšlel jsem o možnosti využít parcely z RUIAN. Tam je nejen zemědělská půda, ale i leccos ostatního - zastavěná plocha, silnice, dálnice a co já vím, co ještě. Zemědělská půda tam určitě není tak podrobně jako v eagri/lpis. Máš představu, nakolik jsou údaje trvanlivé, t.j. zda příští rok ještě budou platit? Jak aktualizovat? Jak to spojit s jiným landuse, neb nejen zemědělskou půdou povrch naší matičky pokryt jest. Dle právě provedeného výpočtu má v RUIAN digitální mapu 64 procent území republiky, takže sice dobré, ale ne dost. Určitě najdeme spoustu míst, kde podle RUIAN je sídliště a podle EAGRI louka či pole (po zkušenostech s importem adres ;-)). Nakolik budeme která data brát jako referenční či jak se rozhodneme, co převezmeme. Možná existují i další zdroje pro landuse.? Co budeme dělat se stávajícími polygony landuse. Dělat díry v polygonech z EAGRI? Otázek bude jistě mnoho a mnoho, děkuji, že jsi toto téma rozstřelil. Opravdu se těším na barevnou mapu! -- Petr ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] LPIS import
On Sun 2014-07-20 23:53:26, Pavel Machek wrote:
> Ahoj!
>
> V priloze je skript na import zemedelske pudy z LPIS.
A ted ten skript.
--
(english) http://www.livejournal.com/~pavelmachek
(cesky, pictures)
http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html
#!/usr/bin/python
import shapefile
from pyproj import *
print ""
print ""
#p1 = Proj(init='my:2065')
p1 = Proj(init='esri:102067')
#p1 = Proj(init='EPSG:4326')
p2 = Proj(init='EPSG:4326')
sf = shapefile.Reader("PLPIS_223941_KU_KOD_631035")
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)
return lon, lat
def write_point(point):
global nodeid
lon, lat = convert(point)
tags = ''
tags += ''
nodeid -= 1
print '%s' % ( 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 '' % nodeid
print ' '
print ' ' % attr(shrec, 'ID_FB')
print ' '
kul = attr(shrec, 'KULTURA')
if kul == 2:print ' '
elif kul == 3: print ' '
elif kul == 30: print ' '
elif kul == 31: print ' '
elif kul == 41: print ' '
elif kul == 61: print ' '
elif kul == 62: print ' '
elif kul == 7: print ' '
elif kul == 71: print ' '
elif kul == 72: print ' '
elif kul == 91: print ' '
elif kul == 92: print ' '
elif kul == 97: print ' '
elif kul == 98: print ' '
elif kul == 99: print ' '
else: print ' ' % kul
for pt in pts:
print ' ' % pt
print ''
print ''
___
Talk-cz mailing list
Talk-cz@openstreetmap.org
https://lists.openstreetmap.org/listinfo/talk-cz
[Talk-cz] LPIS import
Ahoj! V priloze je skript na import zemedelske pudy z LPIS. Bohuzel, tak jak to je, tak jsou data posunuta cca o 100 metru. Netusite nekdo, cim by to mohlo byt? (Na druhou stranu, asi nebude problem posunout to zpet pokud je posun tak nejak konstantni...) Mejte se, 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
