Ahoj! Chtel bych importovat melioracni kanaly... Predstavuju si to zhruba takhle: <way id="-1"> <nd ref="-1" /> <nd ref="-2" /> <nd ref="-3" /> <nd ref="-4" /> <nd ref="-5" /> <tag k="source" v="dibavod A04" /> <tag k="waterway" v="ditch" /> <tag k="dibavod:id" v="135" /> </way> (Doufam ze melioracni kanal odpovida ditch...) Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

Ahoj, On Tue 19-10-10 22:05:45, Pavel Machek wrote: [...] zobrazit citaci Cesky preklad Map_Features napovida "waterway=drain", nicmene po precteni anglicke verze[1], srovnani s "drain"[2] a vysvetleni slova "ditch"[3] se priklanim k "waterway=ditch". Libor [1] http://wiki.openstreetmap.org/wiki/Tag:waterway%3Dditch [2] http://wiki.openstreetmap.org/wiki/Tag:waterway%3Ddrain [3] http://dictionary.cambridge.org/dictionary/british/ditch_1

Ahoj! zobrazit citaci Tak melioracni kanaly jsou uploadovane. Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

Ahoj! Chtel bych importovat baziny, a predstavuju si to nejak takhle... <way id="-1000001"> <nd ref="-1" /> <nd ref="-2" /> <nd ref="-3" /> <nd ref="-4" /> <nd ref="-5" /> <nd ref="-6" /> <nd ref="-7" /> <nd ref="-8" /> <nd ref="-9" /> <nd ref="-10" /> <nd ref="-11" /> <nd ref="-12" /> <nd ref="-14" /> <nd ref="-15" /> <nd ref="-16" /> <nd ref="-1" /> <tag k="source" v="dibavod A06" /> <tag k="natural" v="wetland" /> <tag k="dibavod:id" v="42" /> </way> Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

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. -- masox 2010/10/21 Pavel Machek <pavel na ucw.cz> zobrazit citaci ------------- další část --------------- HTML příloha byla odstraněna... URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20101022/a1eb2d04/attachment.html>

2010/10/22 Jan Masopust <masopust.jan na gmail.com>: zobrazit citaci 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

Právě řeším povodí Malše od Budějovic proti proudu, konroluju podle cenia ortofota, a zatím jediný meliorační kanál nesedí. Buď vede souběžně lépe pasující potok, nebo vede napříč polem či dokonce zástavbou a na fotce není ani známka. Tedy všechny zatím mažu. BTW, co s potokem, který už neexistuje (zastavěná oblast)? J&D Dne 22. října 2010 17:37 Michal Grézl <michal.grezl na openstreetmap.cz> napsal(a): zobrazit citaci -- -- Ing. Jan Dudík

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

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

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ádra<9E>í) Tok naz : Plačický potok Voda typ : P Hlgp vuv : 103010170 Tok id : 104550000100 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('<iddddii', 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] / 1000000 y = pt[1] / 1000000 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))) pom = (pom - 1) / (pom + 1) if abs(pom - sinB) < 1e-15: break Bjtsk = math.atan(pom / math.sqrt(1 - pom * pom)) # Pravouhle souradnice ve S-JTSK a = 6377397.15508 f_1 = 299.152812853 e2 = 1 - (1 - 1 / f_1) * (1 - 1/f_1) ro = a / math.sqrt(1 - e2 * math.sin(Bjtsk) * math.sin(Bjtsk)) x = (ro + H) * math.cos(Bjtsk) * math.cos(Ljtsk) y = (ro + H) * math.cos(Bjtsk) * math.sin(Ljtsk) z = ((1 - e2) * ro + H) * math.sin(Bjtsk) # Pravouhle souradnice v WGS-84 dx = 570.69 dy = 85.69 dz = 462.84 wz = -5.2611 / 3600 * math.pi / 180 wy = -1.58676 / 3600 * math.pi / 180 wx = -4.99821 / 3600 * math.pi / 180 m = 3.543e-6 xn = dx + (1 + m) * (x + wz * y - wy * z) yn = dy + (1 + m) * (-wz * x + y + wx * z) zn = dz + (1 + m) * (wy * x - wx * y + z) # Geodeticke souradnice v systemu WGS-84 a = 6378137.0 f_1 = 298.257223563 a_b = f_1 / (f_1 - 1) p = math.sqrt(xn * xn + yn * yn) e2 = 1 - (1 - 1 / f_1) * (1 - 1/f_1) theta = math.atan(zn * a_b / p) st = math.sin(theta) ct = math.cos(theta) t = (zn + e2 * a_b * a * st * st * st) / (p - e2 * a * ct * ct * ct) B = math.atan(t) L = 2 * math.atan(yn / (p + xn)) H = math.sqrt(1 + t * t) * (p - a / math.sqrt(1 + (1 - e2) * t * t)) # Format vystupnich hodnot return (180 * L / math.pi , 180 * B / math.pi) def createTables(cur): cur.execute("""CREATE TABLE IF NOT EXISTS Nodes(ID_NODE INTEGER PRIMARY KEY AUTOINCREMENT, X INTEGER NOT NULL, Y INTEGER NOT NULL)""") cur.execute("""CREATE TABLE IF NOT EXISTS Lines(ID_LINE INTEGER PRIMARY KEY AUTOINCREMENT, ID_DIBAVOD INTEGER)""") cur.execute("""CREATE TABLE IF NOT EXISTS Areas(ID_AREA INTEGER PRIMARY KEY AUTOINCREMENT, ID_DIBAVOD INTEGER)""") cur.execute("""CREATE TABLE IF NOT EXISTS Relations(ID_RELATION INTEGER PRIMARY KEY AUTOINCREMENT)""") cur.execute("""CREATE TABLE IF NOT EXISTS Dibavod(ID_DIBAVOD INTEGER PRIMARY KEY AUTOINCREMENT, NAME TEXT, OBJ_ID INTEGER NOT NULL)""") cur.execute("""CREATE TABLE IF NOT EXISTS LineNodes(ID_LINE INTEGER, ID_NODE INTEGER, SEQ INTEGER)""") cur.execute("""CREATE TABLE IF NOT EXISTS LineRelations(ID_RELATION INTEGER, ID_LINE INTEGER, SEQ INTEGER)""") cur.execute("""CREATE TABLE IF NOT EXISTS AreaNodes(ID_AREA INTEGER, ID_NODE INTEGER, SEQ INTEGER)""") cur.execute("""CREATE TABLE IF NOT EXISTS AreaRelations(ID_RELATION INTEGER, ID_AREA INTEGER, SEQ INTEGER)""") def importVody(prefix, nameColumn, idColumn): global allNodes con = sqlite3.connect("vody.sqlite") con.isolation_level = None cur = con.cursor() createTables(cur) cur.execute("BEGIN") total = 0 allLines = [] lastTyp = None for (number, typ, obj) in reader(os.path.join('data', prefix + '.shp')): total += 1 if typ in (POLYLINE_SHAPE, POLYGON_SHAPE): if typ == POLYLINE_SHAPE: tablePrefix = "Line" elif typ == POLYGON_SHAPE: tablePrefix = "Area" if lastTyp is None: lastTyp = typ assert lastTyp == typ lineIDs = [] for line in obj: nodeIDs = [] if typ == POLYGON_SHAPE: line = line[:-1] for pt in line: pt2 = jtsk2wgs84(-pt[1], -pt[0]) pt2 = (int(pt2[0] * 10000000 + 0.5), int(pt2[1] * 10000000 + 0.5)) if pt2 in allNodes: nodeID = allNodes[pt2] else: cur.execute("insert into Nodes(X, Y) values (?, ?)", pt2) nodeID = cur.lastrowid allNodes[pt2] = nodeID if len(nodeIDs) > 0 and nodeID == nodeIDs[-1]: pass else: nodeIDs.append(nodeID) cur.execute("insert into %ss(ID_DIBAVOD) values (NULL)" % tablePrefix) lineID = cur.lastrowid lineIDs.append(lineID) for i, nodeID in enumerate(nodeIDs): cur.execute("insert into %sNodes(ID_%s, ID_NODE, SEQ) values (?, ?, ?)" % (tablePrefix, tablePrefix), (lineID, nodeID, i + 1)) allLines.append(lineIDs) if len(obj) > 1: print "Relations: ", len(obj) cur.execute("insert into Relations VALUES(NULL)") relationID = cur.lastrowid for i, lineID in enumerate(lineIDs): cur.execute("insert into %sRelations(ID_RELATION, ID_%s, SEQ) values (?, ?, ?)" % (tablePrefix, tablePrefix), (relationID, lineID, i + 1)) # break elif typ == POINT_SHAPE: ( x, y ) = obj ( lon, lat ) = jtsk2wgs84(-y, -x) print "<node id='%d' lon='%f' lat='%f' />" % ( -total, lon, lat ) else: raise Exception('Unknown shape') #if total >= 10: # break #break if total % 200 == 0: print total print total f = open(os.path.join('data', prefix + '.dbf'), "rb") db = dbf.dbfreader(f) fieldnames = db.next() fieldspecs = db.next() if nameColumn != '': inazev = fieldnames.index(nameColumn) else: inazev = None iid = fieldnames.index(idColumn) i = 0 for row in db: n = None if inazev: n = unicode(row[inazev], 'cp1250').strip() if len(n) == 0: n = None did = long(row[iid]) cur.execute("INSERT INTO Dibavod(NAME, OBJ_ID) VALUES (?, ?)", (n, did)) idDibavod = cur.lastrowid for line in allLines[i]: cur.execute("UPDATE %ss SET ID_DIBAVOD=? WHERE ID_%s=?" % (tablePrefix, tablePrefix), (idDibavod, line)) i += 1 if i % 200 == 0: print i if i == len(allLines): break f.close() cur.execute("COMMIT") con.close() if __name__ == '__main__': allNodes = {} if 1: print "<?xml version='1.0' encoding='UTF-8'?>" print "<osm version='0.5' generator='shapefile'>" #importVody('A02_Vodni_tok_JU', 'NAZ_TOK', 'UTOKJ_ID') #importVody('A05_Vodni_nadrze', 'NAZ_NA', 'NADR_GID') #importVody('A01_Vodni_tok_CEVT', 'NAZ_TOK', 'TOK_ID') #importVody('A04zvm_Melioracni_kanaly', '', 'LCRM_ID') #importVody('A06_Bazina_mocal', '', 'ET_ID') #importVody('I01zvm_Jezy', 'TOK_ID', 'OT06_ID') importVody('C02_Koupaci_oblasti', '', '') -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

Pavel Machek napsal(a): zobrazit citaci 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 ------------- další část --------------- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 261 bytes Desc: OpenPGP digital signature URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20101024/a3f84749/attachment.sig>

Ahoj! zobrazit citaci No, ale nic lepsiho zatim neni. Vyrobit novy tag my neprijde vhodny. zobrazit citaci Ten uz tam je - typicky rybnik nebo reka. zobrazit citaci Zatim mam: <?xml version='1.0' encoding='UTF-8'?> <osm version='0.5' generator='shapefile'> <node id='-1' lon='15.793556' lat='50.165569'> <tag k='name' v='písník Březhrad (u nádraží)' /> <tag k='source' v='dibavod C02 koupaci oblasti' /> <tag k='dibavod:kraj' v='PARDUBICKÝ' /> <tag k='dibavod:tok_naz' v='Plačický potok' /> <tag k='dibavod:voda_typ' v='P' /> <tag k='dibavod:tok_id' v='104550000100' /> <tag k='dibavod:pozn' v='není přímé spojení s vodním tokem' /> <tag k='dibavod:orp' v='Pardubice' /> <tag k='dibavod:obec' v='Opatovice nad Labem' /> <tag k='sport' v='swimming' /> </node> ...co takhle pridat swimming=natural? Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

Alespoň to, případně něco jiného. Každopádně by mělo jít poznat, že to je přírodní koupaliště, abychom byli schopní to převést na nějaký konkrétnější tag, až vznikne (i kdybychom tam měli přidat nějaký vymyšlený tag). Tag sport=swimming tam samozřejmě může zůstat. On Sun, 24 Oct 2010 19:09:44 +0200, Pavel Machek <pavel na ucw.cz> wrote: zobrazit citaci -- Petr Dlouhý

On Sun 2010-10-24 20:04:47, Petr Dlouhý wrote: zobrazit citaci Je to na serveru :-). Bohuzel do uploadu bazin nekdo zasahnul, takze to pujde hur :-(. Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

Dne 24.10.2010 23:54, Pavel Machek napsal(a): zobrazit citaci To chce ty data rozdrbavat na pidi kousky, aby import bloku netrval dyl nez par minut ... Jeste importujes nebo uz plati "krles" ??? Jinak sem narazil namatkou na useky potoka, kde nesedi nazev importovany a nazev vygrabnuty z km. A ted babo rad, protoze az takovy znalec mistnich nazvu nejsem.

2010/10/25 jzvc <jzvc na tpfree.fdns.net>: zobrazit citaci to uz sem zazil take s temi nazvy potoku, ja bych tam nechal oba v nejake rozeznatelne forme, treba name:dibavod a name:katastr nebo taknejak. A do name nektery z nich vybrat. -- Michal Grézl http://openstreetmap.cz

Ahoj! zobrazit citaci Uz jsem to nejak vybojoval. 2 baziny se neuploadly, cert je vem. zobrazit citaci To neni uplne trivialni, protoze na 'svech' pak vznikaji nespojeny cesty co je potreba rucne upravit. Hezky by bylo kdyby josm uploadnul cestu vzdy ve chvili kdy ma potrebne body. Bohuzel to nedela :-(. (Namet na zapoctak: udelatko co to delat bude :-). zobrazit citaci Mezi lon 14 a 15 jeste budu par potoku importovat, ale uz je mozny opravovat. Ty co uz tam jsou jsou finalni. Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

Ahoj! zobrazit citaci No, ta data fakt jsou nejaka divna. Ty soubehy s potokama jsou divny. Na druhou stranu melioracni ryhu si predstavuju... jako ryhu ve ktery obcas je voda. Necekal bych ze to bude videt na ortofotu. Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

no já jsem zatím všechny importované meliorační kanály v okolí nechanic smazal, nebot na uvedených místech ve skutečnosti nic není, ani ta tebou zmiňovaná rýha. Možná se bude jednat o meliorační trubky položené do země někdy v 50 letech Pražák zobrazit citaci

Dne 25. října 2010 22:36 Pavel Machek <pavel na ucw.cz> napsal(a): zobrazit citaci No, ona by i ta rýha měla být vidět na ortofoto, protože přes rýhu traktory neorají ani nesečou. Už jsem ale narazil na několik, které možná k něčemu budou (v parcích) J&D -- -- Ing. Jan Dudík

Ahoj, jak to vypadá? U nás je spousta chyb a stále spousta nezapojených node. Už můzeme opravovat bez omezení nebo se na importu ještě pracuje? Mike On 25.10.2010 15:25, Pavel Machek wrote: zobrazit citaci