Re: [Talk-cz] Hranice obcí
2010/1/18 Mike Crash : > To je šílené schéma, vůbec si nedovedu představit, jak tohle dělat v OSM. *** to schéma jenom říká, že ruzna varianty kardinality 1:1 až 1:N funguji pouze u: NUTS0->NUTS1->NUTS2->NUTS3->NUTS4->LAU1->LAU2->katastralni uzemi(k.u.) * okres by slo nahradit POU a ORP (coz je mozna zbytecne) * take nektere NUTS by slo vynechat, nebot maji pro CR kardinalitu vzdy 1:1 to je na zamysleni > Pokusně jsem označil dvě vesnice Borkovice a Mažice na [1], je to tak > správně? Jde o hranice katastrálních území stejného jména. *** chybne to neni, ale bylo by velmi dobre mit vsechny hranice jako admin_level = 10 tj. k.u. Teprve z nich tvorit nadrazene hranice, tj. obce, okresy ... a nejaky zdroj a znacku o stavu source=cuzk:prehledky date:import_state=2009-06-21 tagy jsem dopsal na wiki http://wiki.openstreetmap.org/wiki/User:Hanoj/grass hanoj ___ Talk-cz mailing list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Hranice obcí
To je šílené schéma, vůbec si nedovedu představit, jak tohle dělat v OSM. Pokusně jsem označil dvě vesnice Borkovice a Mažice na [1], je to tak správně? Jde o hranice katastrálních území stejného jména. [1] http://www.openstreetmap.org/?lat=49.2156&lon=14.6364&zoom=13&layers=B000FTF On 17.1.2010 14:01, hanoj wrote: > 2010/1/17 Martin Kupec : >> Bohuzel u nas se vyskytuje jak, ze obec ma vic katastralnich uzemi, tak, >> ze obec je fyzicky mensi nez cele jeji katastralni uzemi. > *** tak to Martine myslim neni, ani nebylo to takto popsano u Broumovska: > * vazba 1:N obec : katastralni uzemi > * vzaba N:1 casti_obce : katastralni uzemi > > >>> jakým způsobem se má definovat hranice obce tagy boundary a admin_level? >>> Nevím, jestli má být hranice obce podle katastrální mapy, takže bude >>> mnohem větší než samotná obec, nebo se má obtáhnout jen obec podle UHUL? > *** rozlisujte "zastavene uzemi obce" a "spravni uzemi obce". > "zastavene uzemi obce" je jen urbanisticky termin, z pohledu > administrativniho cleneni nema zadny vyznam. > > a nastudujte si: > http://www.czso.cz/csu/rso.nsf/i/schema_soustavy > > zdravi > hanoj > > ___ > Talk-cz mailing list > Talk-cz@openstreetmap.org > http://lists.openstreetmap.org/listinfo/talk-cz > ___ Talk-cz mailing list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import adres z katastralni mapy
Tak jsem na wiki [1] vytvoril stranku s informacemi o importu. Zaroven jsem jsem updatoval program [2], tak aby generoval tagy FIXME pro ty adresni body, kde prirazeni mezi databazi a mapou neni jednoznacne. >Martin Kupec >Dalsi faze bude sehnat si od hanoje vektorizovane obrysy >katastralnich uzemi(nejak bojuju s GRASSem, takze si to nejak >nejsem schopen udelat sam) a pospojovat je na polygony a pridat >jim podle polohy nazvy. S mapou od hanoje jsem si dneska hral a vysledek vypada celkem nadejne, takze mapu s polygony snad budu schopny dodat. Chci to jen jeste trochu doladit ... [1] http://wiki.openstreetmap.org/wiki/Import_Adres_ČR [2] http://lkabrt.aspone.cz/osm/cuzk.zip ___ Talk-cz mailing list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import dibavodu
> tak jsem na to koukal a nebude to tak jednoduche jak si asi vsichni myslime. > OSM nejak strasne moc pokrocilo, protoze co jsem tak nahodne koukal tak > vetsi nadrze uz tam skoro vsechny jsou. Bude to myslim velka prace to > procistit. Chceme tedy import delat drive nez nejakou tu detekci hotoveho? *** Pavel rad skace po hlave do metrove vody. > Asi by bylo dobre udelat 2 datasety. Jeden kde se nic neprotina se muze > naimportovat ihned a druhy se muze po castech rucne zkontrolovat a az pote > importovat. Jinak kvalitativne mi prisel lepsi Dibavod, kdyz jsem videl to > same v OSM, ale tezko rict. 2 sety dobra cesta zdravi hanoj ___ Talk-cz mailing list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import dibavodu
Ahoj, tak jsem na to koukal a nebude to tak jednoduche jak si asi vsichni myslime. OSM nejak strasne moc pokrocilo, protoze co jsem tak nahodne koukal tak vetsi nadrze uz tam skoro vsechny jsou. Bude to myslim velka prace to procistit. Chceme tedy import delat drive nez nejakou tu detekci hotoveho? Asi by bylo dobre udelat 2 datasety. Jeden kde se nic neprotina se muze naimportovat ihned a druhy se muze po castech rucne zkontrolovat a az pote importovat. Jinak kvalitativne mi prisel lepsi Dibavod, kdyz jsem videl to same v OSM, ale tezko rict. Dale tagy davam k way takoveto (pro nadrze): Prijde mi dobre dat tam i to A05, ale neni to asi nutnost kdyz se bude hledat i podle reservoir. Co myslite? Jinak reky budou velky problem, ktery zatim netusim jak udelat. Waterway je ok, ale riverbank oni nemaji. Maji skutecne jen linie brehu a jeste jsem nevykoukal jak je spojit dohromady do polygonu. Nevidim tam relace. Takze bych urcite zacal nadrzemi, tam je polygonu 71922 a 1483946 nodu. Reknete co s tou duplicitou. Tomas Pavel Machek napsal(a): Ahoj! ten sqlite neni zadnej standard, to jsem si udelal vlastni format, aby se s tim dalo lepe pracovat... To co je v priloze jsem jeste nekde na disku nasel. Nemam se toho ujmout? Nerad bych, aby to dopadlo spatne. Ted se mi to nechce cist jestli uz to delas, ale jestli budes chtit muzu se na to podivat. Jestli to udelas budu velmi rad. Nakonec bych to asi zvladnul, ale Ty s tim mas podstatne vic zkusenosti. Take jsem premyslel o tom oznaceni vod, ktere jsou v pruniku. Asi by to nemusel byt problem. Slo by to udelat tak, ze nactu vsechny vody v cechach a pres nejaky ctvereckovy spatialindex se udela prunik. Pokud bude prunik s dibavod, muze se entita oznacit nejakym atributem (obsolete=true), ktery si pak kazdy muze vyfiltrovat a pripadne smazat nebo naopak atribut odstranit. To bych take asi mohl udelat. Ja myslim ze to oznaceni duplicit je dalsi krok ktery pujde udelat pres xapi kdykoliv pozdejc... a mozna bych s tim import nezdrzoval... Pavel ___ Talk-cz mailing list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Hranice obcí
> Dne neděle 17 Leden 2010 14:01:57 hanoj napsal(a): >> *** rozlisujte "zastavene uzemi obce" a "spravni uzemi obce". >> "zastavene uzemi obce" je jen urbanisticky termin, z pohledu >> administrativniho cleneni nema zadny vyznam. > > Dobře, a podle čeho máme teda "obtáhnout" a jakým způsobem máme otagovat > hranici města nebo vesnice, tak aby se vykreslila při nízkých zoomech > mapnikem, tak jak je to mapříklad u Olomouci? *** bud musis specifikovat jakou hranici z vyse popsanych mas na mysli, nebo poslat href. To v tve otazce neni... jediny zdroj hranic, ktery mame, v tomto pripade katastralnich (vrstva prehledky): http://wms.cuzk.cz/wms.asp?SERVICE=WMS&request=GetCapabilities to co z katastralnich hranic muzeme sjednocenim udelat je popsano zde: http://www.czso.cz/csu/rso.nsf/i/schema_soustavy vysledek sjednoceni je popsan tady(zjednodusene, Praha si veci rada komplikuje): http://wiki.openstreetmap.org/wiki/Key:admin_level#admin_level hanoj ___ Talk-cz mailing list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Pár dotazů, názorů
OK, vypadá to, že to tam už někdo přidal. Je možné, že se zatím nebudou vykreslovat, ale každopádně máme standard. Já to tedy ve svých skriptech změním. On Sun, 17 Jan 2010 18:26:19 +0100, Vojta wrote: > Zdravím, na http://topo.geofabrik.de/symbols_en.html se objevily i české > symboly. Vzhledem k tomu že je na to odkaz přímo v textu na Wiki, tak > počítám že se to může v klidu využívat..? Jdu to zkusit na pár svých > trasách. -- Petr Dlouhý ___ Talk-cz mailing list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Pár dotazů, názorů
Zdravím, na http://topo.geofabrik.de/symbols_en.html se objevily i české symboly. Vzhledem k tomu že je na to odkaz přímo v textu na Wiki, tak počítám že se to může v klidu využívat..? Jdu to zkusit na pár svých trasách. Vojta Pro některé české symboly tam zatím nic není, psal jsem to do diskuze na té stránce, ale zatím to nikdo nepřidal. Můžeme to na tu stránku přidat, ale napřed bych se shodl na názvech těch symbolů. Ve svém skriptu jsem dočasně použil názvy "L", "bottom_halfcircle", "right_up_triangle" a "upside_T", ale je to první co mě napadlo, tak jestli někdo vymyslí lepší názvy, tak bych byl rád. Jinak tvar těch symbolů je na [1], takže ten překlad není až tak nutný. [1] http://topo.geofabrik.de/symbols_en.html ___ Talk-cz mailing list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Hranice obcí
Dne neděle 17 Leden 2010 14:01:57 hanoj napsal(a): > *** rozlisujte "zastavene uzemi obce" a "spravni uzemi obce". > "zastavene uzemi obce" je jen urbanisticky termin, z pohledu > administrativniho cleneni nema zadny vyznam. Dobře, a podle čeho máme teda "obtáhnout" a jakým způsobem máme otagovat hranici města nebo vesnice, tak aby se vykreslila při nízkých zoomech mapnikem, tak jak je to mapříklad u Olomouci? -- Ing. Radek Bartoň Faculty of Information Technology Department of Computer Graphics and Multimedia Brno University of Technology E-mail: black...@post.cz Web: http://blackhex.no-ip.org Jabber: black...@jabber.cz ___ Talk-cz mailing list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import dibavodu
Ahoj! >ten sqlite neni zadnej standard, to jsem si udelal vlastni format, aby se >s tim dalo lepe pracovat... To co je v priloze jsem jeste nekde na disku >nasel. Nemam se toho ujmout? Nerad bych, aby to dopadlo spatne. Ted se mi >to nechce cist jestli uz to delas, ale jestli budes chtit muzu se na to >podivat. Jestli to udelas budu velmi rad. Nakonec bych to asi zvladnul, ale Ty s tim mas podstatne vic zkusenosti. >Take jsem premyslel o tom oznaceni vod, ktere jsou v pruniku. Asi by to >nemusel byt problem. Slo by to udelat tak, ze nactu vsechny vody v cechach >a pres nejaky ctvereckovy spatialindex se udela prunik. Pokud bude prunik >s dibavod, muze se entita oznacit nejakym atributem (obsolete=true), ktery >si pak kazdy muze vyfiltrovat a pripadne smazat nebo naopak atribut >odstranit. To bych take asi mohl udelat. Ja myslim ze to oznaceni duplicit je dalsi krok ktery pujde udelat pres xapi kdykoliv pozdejc... a mozna bych s tim import nezdrzoval... 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 http://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Hranice obcí
2010/1/17 Martin Kupec : > Bohuzel u nas se vyskytuje jak, ze obec ma vic katastralnich uzemi, tak, > ze obec je fyzicky mensi nez cele jeji katastralni uzemi. *** tak to Martine myslim neni, ani nebylo to takto popsano u Broumovska: * vazba 1:N obec : katastralni uzemi * vzaba N:1 casti_obce : katastralni uzemi >> jakým způsobem se má definovat hranice obce tagy boundary a admin_level? >> Nevím, jestli má být hranice obce podle katastrální mapy, takže bude >> mnohem větší než samotná obec, nebo se má obtáhnout jen obec podle UHUL? *** rozlisujte "zastavene uzemi obce" a "spravni uzemi obce". "zastavene uzemi obce" je jen urbanisticky termin, z pohledu administrativniho cleneni nema zadny vyznam. a nastudujte si: http://www.czso.cz/csu/rso.nsf/i/schema_soustavy zdravi hanoj ___ Talk-cz mailing list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import dibavodu
Ahoj,
ten sqlite neni zadnej standard, to jsem si udelal vlastni format, aby
se s tim dalo lepe pracovat... To co je v priloze jsem jeste nekde na
disku nasel. Nemam se toho ujmout? Nerad bych, aby to dopadlo spatne.
Ted se mi to nechce cist jestli uz to delas, ale jestli budes chtit muzu
se na to podivat.
Jinak ano byl to dibavod.
Take jsem premyslel o tom oznaceni vod, ktere jsou v pruniku. Asi by to
nemusel byt problem. Slo by to udelat tak, ze nactu vsechny vody v
cechach a pres nejaky ctvereckovy spatialindex se udela prunik. Pokud
bude prunik s dibavod, muze se entita oznacit nejakym atributem
(obsolete=true), ktery si pak kazdy muze vyfiltrovat a pripadne smazat
nebo naopak atribut odstranit. To bych take asi mohl udelat.
Tomas
Pavel Machek napsal(a):
Ahoj!
Kdysi davno jsi psal:
- Importovany linie jednotlivych toku 5.212.525 bodu 250.612 linii
- Importovany nadrze 1.486.406 bodu 72.026 ploch 822 multipolygonu z
toho 1.522 der (ostrovu)
- Importovany nazvy vcetne identifikatoru pro pripadny merge v
- budoucnu
- Sjednoceny sdilene nody (stejna lokace), ktere neni mozne
simplifikovat 187.921 bodu
- Export relaci pro polygony. Smer polygonu urcuje outer/inner.
- Simplifikace dle zadane hodnoty
Subject byl "import dat povodi Labe" -- jsou to ta sama data jako
dibavod (http://www.dibavod.cz/index.php?id=27) nebo ne?
Python skripty jsem rozchodil, ale ony na vstupu potrebuji .sqlite, a
ja to neumim ze shapefilu do .sqlite dostat :-(. Podle nejakyho navodu
jsem zkusil:
spatialite> .loadshp A04zvm_Melioracni_kanaly A04 UTF-8 2065 the_geom
load shapefile error: cannot open shapefile 'A04zvm_Melioracni_kanaly'
cause: 'A04zvm_Melioracni_kanaly.dbf' contains unsupported
data types
spatialite>
Jeste me napadlo oficialne kontaktovat urad s dotazem o tom
zakonu. Oni
musi oficialne odpovedet a pak bysme 100% vedeli jak to tedy je. Mohl
by
se toho nekdo ujmout?
Zda se ze tohle zvladnul hanoj :).
Pavel
import struct, datetime, itertools, string
def dbfreader(f):
"""Returns an iterator over records in a Xbase DBF file.
The first row returned contains the field names.
The second row contains field specs: (type, size, decimal places).
Subsequent rows contain the data records.
If a record is marked as deleted, it is skipped.
File should be opened for binary reads.
"""
# See DBF format spec at:
# http://www.pgts.com.au/download/public/xbase.htm#DBF_STRUCT
numrec, lenheader, lenrecord = struct.unpack('import struct, dbf, cPickle, time
import sqlite3, os.path, math
NULL_SHAPE = 0
POINT_SHAPE = 1
POLYLINE_SHAPE = 3
POLYGON_SHAPE = 5
def pnInPoly(pts, pt):
c = False
j = len(pts) - 1
for i in xrange(len(pts)):
if ((pts[i][1] <= pt[1]) and (pt[1] < pts[j][1])) or ((pts[j][1] <=
pt[1]) and (pt[1] < pts[i][1])):
if pt[0] < (float(pts[j][0] - pts[i][0]) * (pt[1] - pts[i][1]) /
(pts[j][1] - pts[i][1]) + pts[i][0]):
c = not c
j = i
return c
def reader(filename, records = -1):
f = open(filename, 'rb')
f.seek(100)
while 1:
try:
(number, length) = struct.unpack('>ii', f.read(8))
except:
print "end of file..."
break
record = f.read(length * 2)
if ord(record[0]) == NULL_SHAPE:
# Null shape
assert (len(record) == 4)
yield (number, 0, None)
elif ord(record[0]) == POINT_SHAPE:
# Point shape
assert (len(record) == 20)
(typ, x, y) = struct.unpack('= len(parts):
stop = -1
else:
stop = parts[i + 1]
current = parts[i]
polygonpart = []
while current != stop:
if current >= len(points):
break
polygonpart.append(points[current])
current += 1
polygon.append(polygonpart)
yield (number, typ, polygon)
else:
raise Exception('Unknown shape')
records -= 1
if records == 0:
break
f.close()
def isIn(index, pt):
x = pt[0] / 100
y = pt[1] / 100
for poly in index.get((x, y), []):
for p in poly[1]:
if pnInPoly(p, pt):
return poly[0]
return None
def jtsk2wgs84(X, Y):
# Prepocet vstupnich udaju
H = 245
# Vypocet zemepisnych souradnic z rovinnych souradnic
a = 6377397.15508
e = 0.081696831215303
n = 0.97992470462083
konst_u_ro = 12310230.12797036
sinUQ = 0.863499969506341
cosUQ = 0.504348889819882
sinVQ = 0.420215144586493
cosVQ = 0.907424504992097
alfa = 1.000597498371542
k = 1.003419163966575
ro = math.sqrt(X * X + Y * Y)
epsilon = 2 * math.atan(Y / (ro + X))
D = ep
Re: [Talk-cz] Hranice obcí
Bohuzel u nas se vyskytuje jak, ze obec ma vic katastralnich uzemi, tak, ze obec je fyzicky mensi nez cele jeji katastralni uzemi. Nedavno to tady nekdo rozebiral pro Broumovsko. Martin Kupec On Sun, Jan 17, 2010 at 12:25:19PM +0100, Mike Crash wrote: > Ahoj, > > jakým způsobem se má definovat hranice obce tagy boundary a admin_level? > Nevím, jestli má být hranice obce podle katastrální mapy, takže bude > mnohem větší než samotná obec, nebo se má obtáhnout jen obec podle UHUL? > Pak to ale nebude sedět v admin_level (viz [1]), kde by to mělo být od > největšího k nejmenšímu. Tam je katastrální území admin_level=10, ale > obec je admin_level=8, ale bylo by menší než katastrální území. Takže > jedno z toho je špatně. > > Např. v Praze je to tak, že celá Praha je admin_level=8 a části (např. > Běchovice) je 10, s čímž souhlasím. U malých měst a vesnic tedy bude > hranice obce (=8) podle katastrální mapy a katastrální území (=10) jako > takové nebude definováno? > > [1] http://wiki.openstreetmap.org/wiki/Key:admin_level#admin_level > > ___ > Talk-cz mailing list > Talk-cz@openstreetmap.org > http://lists.openstreetmap.org/listinfo/talk-cz ___ Talk-cz mailing list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz
[Talk-cz] Hranice obcí
Ahoj, jakým způsobem se má definovat hranice obce tagy boundary a admin_level? Nevím, jestli má být hranice obce podle katastrální mapy, takže bude mnohem větší než samotná obec, nebo se má obtáhnout jen obec podle UHUL? Pak to ale nebude sedět v admin_level (viz [1]), kde by to mělo být od největšího k nejmenšímu. Tam je katastrální území admin_level=10, ale obec je admin_level=8, ale bylo by menší než katastrální území. Takže jedno z toho je špatně. Např. v Praze je to tak, že celá Praha je admin_level=8 a části (např. Běchovice) je 10, s čímž souhlasím. U malých měst a vesnic tedy bude hranice obce (=8) podle katastrální mapy a katastrální území (=10) jako takové nebude definováno? [1] http://wiki.openstreetmap.org/wiki/Key:admin_level#admin_level ___ Talk-cz mailing list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz
