Re: [Talk-cz] dotaz na ceniaortofotomapu
Dne 22. října 2010 23:40 Jachym Cepicky jachym.cepi...@gmail.com napsal(a): Mám zprávy, že cenia před pár dny vyměnila data taky, pokud vím, je tam arcgis a ten má dost problémy při převodu z jtsk-wgs (tedy v nastavení na cenii určitě). existuje mapa s půlmetrovým pixelem - jakou používáte? *** Adresa s CENIA ortofoto 0,5m/px sice existuje, ale ted je na ni mapa s 1m/px ha hanoj ___ Talk-cz mailing list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] import dibavod A04zvm_Melioracni_kanaly
Pkud se týká melioračních kanálů v okolí Nechanic tak kanály na naimportovaných místech ve skutečnosti nejsou, mhlo by se snad jednat o meliorační potrubí položené někdy v 50 letech. Na povrchu v uvedených místech není nic viděrt aproto jsem tyto kanály smazal Pražák Původní zpráva Od: Michal Grézl michal.gr...@openstreetmap.cz Předmět: Re: [Talk-cz] import dibavod A04zvm_Melioracni_kanaly Datum: 22.10.2010 17:38:10 2010/10/22 Jan Masopust masopust@gmail.com: Ahoj, To je sice hezké, ale na nic. Tyto data jsou dosti nepřesné (asi +- 30m), ale hlavně se kryjí s importovanými potoky. Tzn., že tam, kde bylo zmapováno jsou data třikrát. Nevím, jestli je případ, kde je importován kanál a zároveň není importován potok. Podle mě ne (kdyžtak mě opravte). Takže by se s tím mělo něco udělat. Nevím, jestli stojí za to mít v mapě takto pochybná data. Já bych navrhoval u všech importovaných kanálů změnit tag třeba na waterway=FIXMEditch, aby se to nerendrovalo a pak to vše ručně projít. stejne se to musi projit kvuli potokum, ale fakt sou ty ditch uplne mimo, mam je nechat nebo smazat? bude se s tim neco delat? -- Michal Grézl http://openstreetmap.cz ___ 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] import dibavod I01 -- jezy
Ahoj! Za chvili hodlam zacit importovat jezy. Tenhle import je generovan ponekud jednodussim hackem, ale jsou to jenom body - melo by to byt ok. ?xml version='1.0' encoding='UTF-8'? osm version='0.5' generator='JOSM' node id='-13350' action='modify' timestamp='2010-10-23T20:17:15Z' visible='true' lat='48.98367' lon='14.789454' tag k='source' v='dibavod I01' / tag k='waterway' v='weir' / /node node id='-13348' action='modify' timestamp='2010-10-23T20:17:15Z' visible='true' lat='48.986945' lon='14.778233' tag k='source' v='dibavod I01' / tag k='waterway' v='weir' / /node -- (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] dotaz na ceniaortofotomapu
Mám zprávy, že cenia před pár dny vyměnila data taky, pokud vím, je tam arcgis a ten má dost problémy při převodu z jtsk-wgs ... j Zdeněk Pražák píše v Pá 22. 10. 2010 v 09:29 +0200: V minulosti jsem ke kontrole zakreslení rybníků používal ceniaortofotomapu. Při rozlišení okolo 20 m poskytovala ostré záběry ve kterých byly zřetelně vidět břehy. Včera jsem chtěl zkontrolovat zakreslení nově importovaného potoka a zjistil jsem, že při rozlišení 20 m je mapa neostrá a potoky nejsou při tomto rozlišení zřetelné. Změnilo se u této ortofotomapy něco v nastaveních, že nyní poskytuje tak neostré záběry? Pražák ___ Talk-cz mailing list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz -- Jachym Cepicky e-mail: jachym.cepicky gmail com URL: http://les-ejk.cz PGP Public key: http://les-ejk.cz/pgp/JachymCepicky.pgp ___ Talk-cz mailing list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz
[Talk-cz] import dibavod C02 -- koupaci oblasti
Ahoj!
Uvazuju jak importovat koupaci oblasti.
Tenhle import je generovan ponekud jednodussim hackem, ale jsou to
jenom body - melo by to byt ok.
Aktualni verse je:
?xml version='1.0' encoding='UTF-8'?
osm version='0.5' generator='JOSM'
node id='-264' action='modify' timestamp='2010-10-23T20:41:21Z'
visible='true' lat='49.650985' lon='16.604522'
tag k='note' v='koupaci oblast ve volne prirode dle DIBAVODu' /
tag k='source' v='dibavod C02' /
tag k='sport' v='swimming' /
/node
Je nejake vhodnejsi tagovani?
Hmm v databazi to vypada takhle...
Rekr id: KO530901
Kraj : PARDUBICKÝ
Koup naz : písník Březhrad (u nádra9Eí)
Tok naz: Plačický potok
Voda typ : P
Hlgp vuv : 103010170
Tok id : 10455100
Pozn : není přímé spojení s vodním tokem
Idvt : 10100651
Nadr gid : 103010170003
Orp: Pardubice
Orp id : 5309
Obec : Opatovice nad Labem
Obec id: 575429
...a uvazuju ze ten nazev koupaliste by se asi hodil... Zatim to
vyrabim *strasnym* hackem (prilozen), ktery byl puvodne urcen na
plneni databaze, ne vyrobu xmlka... Bohuzel se v nem moc nevyznam :-(.
Napady?
Pavel
#!/usr/bin/python
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 /osm
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('idd', record)
yield (number, typ, (x, y))
elif ord(record[0]) in (POLYLINE_SHAPE, POLYGON_SHAPE):
# Polygon or polyline shape
(typ, bbleft, bbtop, bbright, bbbottom, numParts, numPoints) =
struct.unpack('iii', record[:44])
record = record[44:]
parts = []
for i in xrange(numParts):
parts.append(struct.unpack('i', record[i*4:(i+1)*4])[0])
record = record[numParts * 4:]
points = []
for i in xrange(numPoints):
points.append(struct.unpack('dd', record[i * 16:(i + 1) * 16]))
polygon = []
for i in xrange(len(parts)):
if i + 1 = 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 = epsilon / n
S = 2 * math.atan(math.exp(1 / n * math.log(konst_u_ro / ro))) - math.pi / 2
sinS = math.sin(S)
cosS = math.cos(S)
sinU = sinUQ * sinS - cosUQ * cosS * math.cos(D)
cosU = math.sqrt(1 - sinU * sinU)
sinDV = math.sin(D) * cosS / cosU
cosDV = math.sqrt(1 - sinDV * sinDV)
sinV = sinVQ * cosDV - cosVQ * sinDV
cosV = cosVQ * cosDV + sinVQ * sinDV
Ljtsk = 2 * math.atan(sinV / (1 + cosV)) / alfa
t = math.exp(2 / alfa * math.log((1 + sinU) / cosU / k))
pom = (t - 1) / (t + 1)
while True:
sinB = pom
pom = t * math.exp(e * math.log((1 + e * sinB) / (1 - e * sinB)))
Re: [Talk-cz] import dibavod C02 -- koupaci oblasti
Pavel Machek napsal(a): Ahoj! Uvazuju jak importovat koupaci oblasti. Tenhle import je generovan ponekud jednodussim hackem, ale jsou to jenom body - melo by to byt ok. Aktualni verse je: ?xml version='1.0' encoding='UTF-8'? osm version='0.5' generator='JOSM' node id='-264' action='modify' timestamp='2010-10-23T20:41:21Z' visible='true' lat='49.650985' lon='16.604522' tag k='note' v='koupaci oblast ve volne prirode dle DIBAVODu' / tag k='source' v='dibavod C02' / tag k='sport' v='swimming' / /node Je nejake vhodnejsi tagovani? Před časem jsem hledal na wiki, jak koupací oblasti tagovat, ale celkem jsem narazil... zdá se, že je to zatím neprošlápnutá oblast. Osobně si myslím, že by se tag sport měl používat pro označení míst, kde se daný sport skutečně provozuje jako sport (aspoň na amaterské úrovni). Tj. sport=swimming je vhodný pro klasický plavečák (leisure=swimming_pool), méně vhodný pro akvapark s tobogány (leisure=water_park) a nevhodný pro písák/lom (???), kam se jezdí lidi v létě vykoupat. Na wiki je obecně ke sportu uvedeno Since this is a non-physical tag it should be combined with one of these (physical) tags, tj. skutečně by se tam měl přihodit i nějaký tag označující objekt, ve kterém se to plavání provozuje. Pro přírodní koupaliště jsem nic rozumného nenašel, možná jedině leisure=swimming_natural podle http://wiki.openstreetmap.org/wiki/Proposed_features/Swimming_pool Ale tagwatch říká, že se to momentálně nepoužívá. Petr signature.asc Description: OpenPGP digital signature ___ Talk-cz mailing list Talk-cz@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-cz
