D, I, II versus Planet071003 na Talk-cz

hanoj enemy na mail.muni.cz #m69e57c 14. 10. 2007 9:24

Ahoj, nize popsan postup mapove algebry na Planet a nove silnice: 0. Jachymova data obsahuji (zaokrouhleno na tisice): 14000 prvku coz je 22000 km. 1. Z planet071003 byly vytazeny highway [1] 2. pomoci obalove zony 50m na kazdou stranu vytvorena ochranna zona. 3. Jachymova data byla rozbita na usecky (184000 prvku) 4. A pote byla dotazovana vuci ochranne zone: [2] * z toho 142000 bylo shledano jako nic nekrizici (zelena) * a 42000 bylo shledano jako konfliktni (cervena) 5. data od Jachyma byla seskupena zpet podle puvodniho deleni na puvodnich 14000 prvku Mozne problemy: 1) cast silnice nebyla importovana ac pouze krizila existujici OSMsilnici, ale sama osobe neni v OSM. 2) casti silnice byly importovany, protoze byla mimo ochrannou zonu OSM, ac jiz tataz v OSM existuje. 3) usek silnice byl navrzen k importu, ale byl rozdelen na vice casti, ac jej reprezentuje jeden prvek v shapefilu. 1 prvek = "-----" po cisteni je 1 prvek = "- ---" 4) spatne serazene segmenty silnic OSM vytvari v shapefile uzavrene okruhy misto lomenych car, tudiz jsou vyrazeny casti, ktere v OSM nejsou. Patrno napr. na silnici I/6 a I/3, vidno i na [1]. [1] http://dhost.info/krchov/files/osm/planet071003-highway.jpg [2] http://dhost.info/krchov/files/osm/import.jpg hoj ha

Michal Grézl michal.grezl na gmail.com #m8ff174 14. 10. 2007 14:55

Ahoj, nize popsan postup mapove algebry na Planet a nove silnice: 0. Jachymova data obsahuji (zaokrouhleno na tisice): 14000 prvku coz je 22000 km. 1. Z planet071003 byly vytazeny highway [1] 2. pomoci obalove zony 50m na kazdou stranu vytvorena ochranna zona. 3. Jachymova data byla rozbita na usecky (184000 prvku) 4. A pote byla dotazovana vuci ochranne zone: [2] * z toho 142000 bylo shledano jako nic nekrizici (zelena) * a 42000 bylo shledano jako konfliktni (cervena) 5. data od Jachyma byla seskupena zpet podle puvodniho deleni na puvodnich 14000 prvku Mozne problemy: 1) cast silnice nebyla importovana ac pouze krizila existujici OSMsilnici, ale sama osobe neni v OSM. 2) casti silnice byly importovany, protoze byla mimo ochrannou zonu OSM, ac jiz tataz v OSM existuje. 3) usek silnice byl navrzen k importu, ale byl rozdelen na vice casti, ac jej reprezentuje jeden prvek v shapefilu. 1 prvek = "-----" po cisteni je 1 prvek = "- ---" 4) spatne serazene segmenty silnic OSM vytvari v shapefile uzavrene okruhy misto lomenych car, tudiz jsou vyrazeny casti, ktere v OSM nejsou. Patrno napr. na silnici I/6 a I/3, vidno i na [1]. [1] http://dhost.info/krchov/files/osm/planet071003-highway.jpg [2] http://dhost.info/krchov/files/osm/import.jpg hoj ha

Takze kdy zacnes s uploadem toho zeleneho? A kde se da stahnout to cervene, aby se mohli vytvorit ty krizici ale neexistujici, o kterych si mluvil v moznych problemech 1)? Vypada to znacne impozantne, asi zacnu setrit na novy pocitac, ponevadz pri tomhle pokryti bude muj stavajici srot nepouzitelny v josm;)

hanoj enemy na mail.muni.cz #mca6bdb 14. 10. 2007 17:45

Takze kdy zacnes s uploadem toho zeleneho? A kde se da stahnout to cervene, aby se mohli vytvorit ty krizici ale neexistujici, o kterych si mluvil v moznych problemech 1)?

*** spatne jsem se vyjadril, ja jsem provedl jen selekci ANO/NE v shapefilu a poskytl ho Pavlovi...

Vypada to znacne impozantne, asi zacnu setrit na novy pocitac, ponevadz pri tomhle pokryti bude muj stavajici srot nepouzitelny v josm;)

*** to asi nepomuze neb ta pomalost tkvi v jave. Mam moznost porovnani vykonu AMD Athalon XP 2GHz s Radeon 9200 (4roky zpet) a Intel Dual Core 2x2.5GHz s Nvidia FX 1500 (1 rok zpet) a zadna vyrazna zmena se nekona, limity (mnozstvi nactenych dat) stejne. ha hanoj

Michal Grézl michal.grezl na gmail.com #m345333 14. 10. 2007 18:43

Vypada to znacne impozantne, asi zacnu setrit na novy pocitac, ponevadz pri tomhle pokryti bude muj stavajici srot nepouzitelny v josm;)
*** to asi nepomuze neb ta pomalost tkvi v jave. Mam moznost porovnani vykonu AMD Athalon XP 2GHz s Radeon 9200 (4roky zpet) a Intel Dual Core 2x2.5GHz s Nvidia FX 1500 (1 rok zpet) a zadna vyrazna zmena se nekona, limity (mnozstvi nactenych dat) stejne.

Java neni pomala, pomalej je program, ja mam 1ghz stroj a napriklad ten omalovankovej plugin je naprosto nepouzitelnej, kdezto na 3ghz to jede uplne normalne. Dalsi rychlostni problem je pri velkem mnozstvi naloudovanych gpx dat, to zase muze byt jen tezko chyba javy, spis spatne napsana metoda prochazejici jednotlive body a vykreslujici je, na rychlejsi masine je mnozstvi gpx ktere to zvladne radove vetsi. Ale to je uz hodne OT.

Pavel Machek pavel na ucw.cz #mfaa79c 14. 10. 2007 21:57

AhoJ!

Ahoj, nize popsan postup mapove algebry na Planet a nove silnice: 0. Jachymova data obsahuji (zaokrouhleno na tisice): 14000 prvku coz je 22000 km. 1. Z planet071003 byly vytazeny highway [1] 2. pomoci obalove zony 50m na kazdou stranu vytvorena ochranna zona. 3. Jachymova data byla rozbita na usecky (184000 prvku) 4. A pote byla dotazovana vuci ochranne zone: [2] * z toho 142000 bylo shledano jako nic nekrizici (zelena) * a 42000 bylo shledano jako konfliktni (cervena) 5. data od Jachyma byla seskupena zpet podle puvodniho deleni na puvodnich 14000 prvku Mozne problemy: 1) cast silnice nebyla importovana ac pouze krizila existujici OSMsilnici, ale sama osobe neni v OSM. 2) casti silnice byly importovany, protoze byla mimo ochrannou zonu OSM, ac jiz tataz v OSM existuje. 3) usek silnice byl navrzen k importu, ale byl rozdelen na vice casti, ac jej reprezentuje jeden prvek v shapefilu. 1 prvek = "-----" po cisteni je 1 prvek = "- ---" 4) spatne serazene segmenty silnic OSM vytvari v shapefile uzavrene okruhy misto lomenych car, tudiz jsou vyrazeny casti, ktere v OSM nejsou. Patrno napr. na silnici I/6 a I/3, vidno i na [1].

Vysledek jsem dal na http://atrey.karlin.mff.cuni.cz/~pavel/outgoing/cervena.osm.bz2 http://atrey.karlin.mff.cuni.cz/~pavel/outgoing/zelena.osm.bz2 Bohuzel je tam nejaky problem -- ve way se nejak opakuji body -- takze to jeste neni na primy import... ...a stejne bude potreba udelat import po castech, dle meho odhadu bych to cele najednou uploadoval zhruba dva mesice. Vzal jsem z puvodnich dat (nemaji problemy s opakovanim bodu) ctverec o 0.5 stupne kolem 51 / 13.5 -- maly, pro testovani. Konflikty jsem vyresil rucne a ted se asi pustim do uploadu. Rozdeleni na ctverecky je spatny -- vzniknou tam chyby na svech -- ale bojim se ze to je jedina realna cesta jak to uploadovat. (Nebo ma nekdo napad?) Pavel

AhoJ!

Ahoj, nize popsan postup mapove algebry na Planet a nove silnice: 0. Jachymova data obsahuji (zaokrouhleno na tisice): 14000 prvku coz je 22000 km. 1. Z planet071003 byly vytazeny highway [1] 2. pomoci obalove zony 50m na kazdou stranu vytvorena ochranna zona. 3. Jachymova data byla rozbita na usecky (184000 prvku) 4. A pote byla dotazovana vuci ochranne zone: [2] * z toho 142000 bylo shledano jako nic nekrizici (zelena) * a 42000 bylo shledano jako konfliktni (cervena) 5. data od Jachyma byla seskupena zpet podle puvodniho deleni na puvodnich 14000 prvku Mozne problemy: 1) cast silnice nebyla importovana ac pouze krizila existujici OSMsilnici, ale sama osobe neni v OSM. 2) casti silnice byly importovany, protoze byla mimo ochrannou zonu OSM, ac jiz tataz v OSM existuje. 3) usek silnice byl navrzen k importu, ale byl rozdelen na vice casti, ac jej reprezentuje jeden prvek v shapefilu. 1 prvek = "-----" po cisteni je 1 prvek = "- ---" 4) spatne serazene segmenty silnic OSM vytvari v shapefile uzavrene okruhy misto lomenych car, tudiz jsou vyrazeny casti, ktere v OSM nejsou. Patrno napr. na silnici I/6 a I/3, vidno i na [1].

Vysledek jsem dal na http://atrey.karlin.mff.cuni.cz/~pavel/outgoing/cervena.osm.bz2 http://atrey.karlin.mff.cuni.cz/~pavel/outgoing/zelena.osm.bz2 Bohuzel je tam nejaky problem -- ve way se nejak opakuji body -- takze to jeste neni na primy import... ...a stejne bude potreba udelat import po castech, dle meho odhadu bych to cele najednou uploadoval zhruba dva mesice. Vzal jsem z puvodnich dat (nemaji problemy s opakovanim bodu) ctverec o 0.5 stupne kolem 51 / 13.5 -- maly, pro testovani. Konflikty jsem vyresil rucne a ted se asi pustim do uploadu. Rozdeleni na ctverecky je spatny -- vzniknou tam chyby na svech -- ale bojim se ze to je jedina realna cesta jak to uploadovat. (Nebo ma nekdo napad?) Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

Pavel Machek pavel na ucw.cz #m377d66 14. 10. 2007 22:29

Ahoj!

Takze kdy zacnes s uploadem toho zeleneho? A kde se da stahnout to

Vzorek je k dispozici v 0.5 stupne okoli 51 / 13.5. Pavel

Michal Grézl michal.grezl na gmail.com #mba670d 15. 10. 2007 15:08

AhoJ!
Ahoj, nize popsan postup mapove algebry na Planet a nove silnice: 0. Jachymova data obsahuji (zaokrouhleno na tisice): 14000 prvku coz je 22000 km. 1. Z planet071003 byly vytazeny highway [1] 2. pomoci obalove zony 50m na kazdou stranu vytvorena ochranna zona. 3. Jachymova data byla rozbita na usecky (184000 prvku) 4. A pote byla dotazovana vuci ochranne zone: [2] * z toho 142000 bylo shledano jako nic nekrizici (zelena) * a 42000 bylo shledano jako konfliktni (cervena) 5. data od Jachyma byla seskupena zpet podle puvodniho deleni na puvodnich 14000 prvku Mozne problemy: 1) cast silnice nebyla importovana ac pouze krizila existujici OSMsilnici, ale sama osobe neni v OSM. 2) casti silnice byly importovany, protoze byla mimo ochrannou zonu OSM, ac jiz tataz v OSM existuje. 3) usek silnice byl navrzen k importu, ale byl rozdelen na vice casti, ac jej reprezentuje jeden prvek v shapefilu. 1 prvek = "-----" po cisteni je 1 prvek = "- ---" 4) spatne serazene segmenty silnic OSM vytvari v shapefile uzavrene okruhy misto lomenych car, tudiz jsou vyrazeny casti, ktere v OSM nejsou. Patrno napr. na silnici I/6 a I/3, vidno i na [1].
Vysledek jsem dal na http://atrey.karlin.mff.cuni.cz/~pavel/outgoing/cervena.osm.bz2 http://atrey.karlin.mff.cuni.cz/~pavel/outgoing/zelena.osm.bz2

Muzu teda otevrit soubor cervena a rucne pridat to co v nem bude navic oproti stavajicimu OSM? Pochopil sem to cele dobre?:)

Bohuzel je tam nejaky problem -- ve way se nejak opakuji body -- takze to jeste neni na primy import... ...a stejne bude potreba udelat import po castech, dle meho odhadu bych to cele najednou uploadoval zhruba dva mesice. Vzal jsem z puvodnich dat (nemaji problemy s opakovanim bodu) ctverec o 0.5 stupne kolem 51 / 13.5 -- maly, pro testovani. Konflikty jsem vyresil rucne a ted se asi pustim do uploadu. Rozdeleni na ctverecky je spatny -- vzniknou tam chyby na svech -- ale bojim se ze to je jedina realna cesta jak to uploadovat. (Nebo ma nekdo napad?) Pavel

no ja bych to treba uploadoval po ctvercich a pak to rucne spojil, nahral bych si tam nejakou referencni caru at je presne videt kde to konci. Nebo ty data rozsekat primo po cestach s patricnymi nody, mam za to, ze prislusny script uz existuje, pokud ne nebyl by problem ho napsat. Tohle by bylo mnohem schudnejsi urcite. Tzn vysledny soubor by obsahoval 1 (ci vice) way + nody co k nemu patri.

On 10/14/07, Pavel Machek <pavel na ucw.cz> wrote:

AhoJ!
Ahoj, nize popsan postup mapove algebry na Planet a nove silnice: 0. Jachymova data obsahuji (zaokrouhleno na tisice): 14000 prvku coz je 22000 km. 1. Z planet071003 byly vytazeny highway [1] 2. pomoci obalove zony 50m na kazdou stranu vytvorena ochranna zona. 3. Jachymova data byla rozbita na usecky (184000 prvku) 4. A pote byla dotazovana vuci ochranne zone: [2] * z toho 142000 bylo shledano jako nic nekrizici (zelena) * a 42000 bylo shledano jako konfliktni (cervena) 5. data od Jachyma byla seskupena zpet podle puvodniho deleni na puvodnich 14000 prvku Mozne problemy: 1) cast silnice nebyla importovana ac pouze krizila existujici OSMsilnici, ale sama osobe neni v OSM. 2) casti silnice byly importovany, protoze byla mimo ochrannou zonu OSM, ac jiz tataz v OSM existuje. 3) usek silnice byl navrzen k importu, ale byl rozdelen na vice casti, ac jej reprezentuje jeden prvek v shapefilu. 1 prvek = "-----" po cisteni je 1 prvek = "- ---" 4) spatne serazene segmenty silnic OSM vytvari v shapefile uzavrene okruhy misto lomenych car, tudiz jsou vyrazeny casti, ktere v OSM nejsou. Patrno napr. na silnici I/6 a I/3, vidno i na [1].
Vysledek jsem dal na http://atrey.karlin.mff.cuni.cz/~pavel/outgoing/cervena.osm.bz2 http://atrey.karlin.mff.cuni.cz/~pavel/outgoing/zelena.osm.bz2

Muzu teda otevrit soubor cervena a rucne pridat to co v nem bude navic oproti stavajicimu OSM? Pochopil sem to cele dobre?:)

Bohuzel je tam nejaky problem -- ve way se nejak opakuji body -- takze to jeste neni na primy import... ...a stejne bude potreba udelat import po castech, dle meho odhadu bych to cele najednou uploadoval zhruba dva mesice. Vzal jsem z puvodnich dat (nemaji problemy s opakovanim bodu) ctverec o 0.5 stupne kolem 51 / 13.5 -- maly, pro testovani. Konflikty jsem vyresil rucne a ted se asi pustim do uploadu. Rozdeleni na ctverecky je spatny -- vzniknou tam chyby na svech -- ale bojim se ze to je jedina realna cesta jak to uploadovat. (Nebo ma nekdo napad?) Pavel

no ja bych to treba uploadoval po ctvercich a pak to rucne spojil, nahral bych si tam nejakou referencni caru at je presne videt kde to konci. Nebo ty data rozsekat primo po cestach s patricnymi nody, mam za to, ze prislusny script uz existuje, pokud ne nebyl by problem ho napsat. Tohle by bylo mnohem schudnejsi urcite. Tzn vysledny soubor by obsahoval 1 (ci vice) way + nody co k nemu patri.

Pavel Machek pavel na ucw.cz #m6a2283 15. 10. 2007 15:58

Ahoj!

0. Jachymova data obsahuji (zaokrouhleno na tisice): 14000 prvku coz je 22000 km. 1. Z planet071003 byly vytazeny highway [1] 2. pomoci obalove zony 50m na kazdou stranu vytvorena ochranna zona. 3. Jachymova data byla rozbita na usecky (184000 prvku) 4. A pote byla dotazovana vuci ochranne zone: [2] * z toho 142000 bylo shledano jako nic nekrizici (zelena) * a 42000 bylo shledano jako konfliktni (cervena) 5. data od Jachyma byla seskupena zpet podle puvodniho deleni na puvodnich 14000 prvku

Tak mam novou versi skriptu, naucil jsem ho vzit z jednoho .shp seznam toho co uploadovat, a z druheho .shp seznam co vynechat, a zda se, ze se chova slusne... Pavel #!/usr/bin/perl use Geo::ShapeFile; # ../shpupload; cat out.osm | cstocs il2 utf8 > result-full.osm %excluded_nodes = {}; createExcludes("sil_cz12x-NE.shx"); open(OUT, ">out.osm"); print OUT "<?xml version='1.0' encoding='UTF-8'?>\n"; print OUT "<osm version='0.5' generator='shpupload'>\n"; %nodes = {}; $nodeid = -1; $way = ""; createCoasts("sil_cz12x.shx", 50.0, 14.5, 0.5); print OUT "</osm>\n"; sub createExcludes() { my ($Filename) = @_; my $shapefile = new Geo::ShapeFile($Filename); printf "%d shapes excluded\n", $shapefile->shapes(); for(1 .. $shapefile->shapes()) { $shapenum = $_; my $shape = $shapefile->get_shp_record($shapenum); $LastValid = 0; print LOG "Shape\n"; foreach $Point($shape->points()){ $Long = $Point->X(); $Lat = $Point->Y(); $Node = sprintf("lon=%f lat=%f", $Long, $Lat); $excluded_nodes{$Node} = 1; } } } sub createCoasts(){ my ($Filename, $Lat, $Long, $Size) = @_; $X1 = $Long - $Size; $X2 = $Long + $Size; $Y1 = $Lat - $Size; $Y2 = $Lat + $Size; open(LOG, ">log.txt"); print LOG "Area $X1 to $X2, $Y1 to $Y2\n"; my $shapefile = new Geo::ShapeFile($Filename); printf "%d shapes\n", $shapefile->shapes(); for(1 .. $shapefile->shapes()) { $shapenum = $_; my $shape = $shapefile->get_shp_record($shapenum); $LastValid = 0; print LOG "Shape\n"; foreach $Point($shape->points()){ $Long = $Point->X(); $Lat = $Point->Y(); $Node = sprintf("lon=%f lat=%f", $Long, $Lat); $InArea = ($Lat > $Y1 && $Lat < $Y2 && $Long > $X1 && $Long < $X2) && !$excluded_nodes{$Node}; if($InArea){ if (!$LastValid) { my %record = $shapefile->get_dbf_record($shapenum); print %record; print "\n"; $way .= "<way id='$nodeid'>\n"; $way .= " <tag k=\"created_by\" v=\"shpupload\"/>\n"; $kod = $record{'HSKOD'}; if (($kod =~ "S1") || ($kod =~ "S1T")) { $way .= " <tag k=\"highway\" v=\"primary\"/>\n"; } if (($kod =~ "S2") || ($kod =~ "S2T")) { $way .= " <tag k=\"highway\" v=\"secondary\"/>\n"; } if ($kod =~ "D[DR].*") { $way .= " <tag k=\"highway\" v=\"motorway\"/>\n"; } if ($kod =~ ".*T") { $way .= " <tag k=\"tunnel\" v=\"true\"/>\n"; $way .= " <tag k=\"layer\" v=\"-1\"/>\n"; } $way .= " <tag k=\"name\" v=\"$record{'NAZEV'}_$record{'OZNACENI'}\"/>\n"; $way .= " <tag k=\"ref\" v=\"$record{'OZNACENI'}\"/>\n"; if ($record{'POZNAMKA'}) { $way .= " <tag k=\"note\" v=\"$record{'POZNAMKA'}\"/>\n"; } $way .= " <tag k=\"source\" v=\"HELP SERVICE - REMOTE SENSING spol. s r.o. http://www.bnhelp.cz\"/>\n"; $nodeid--; } $Node = uploadNode($Lat, $Long); printf LOG "Node #%d: %f, %f\n", $Node, $Lat, $Long; $LastNode = $Node; $CountB++; $LastValid = 1; } else { if ($LastValid) { print OUT "$way</way>\n"; $way = ""; } $LastValid = 0; } $CountA++; } if ($LastValid) { print OUT "$way</way>\n"; $way = ""; } } print "Uploading $CountB of $CountA\n"; print LOG "Complete\n"; close LOG; print "Upload complete\n"; } sub uploadNode(){ ($Lat, $Long) = @_; $Tags = "<tag k=\"created_by\" v=\"shpupload\"/>"; $Tags .= "<tag k=\"source\" v=\"HELP SERVICE - REMOTE SENSING spol. s r.o. http://www.bnhelp.cz\"/>"; $Node = sprintf("lon=%f lat=%f", $Long, $Lat); if ($nodes{$Node} < 0) { $way .= " <nd ref='$nodes{$Node}' />\n"; } else { $way .= " <nd ref='$nodeid' />\n"; $nodes{$Node} = $nodeid; $Node = sprintf("<node id=\"$nodeid\" lon=\"%f\" lat=\"%f\">$Tags</node>", $Long, $Lat); print OUT " $Node\n"; $nodeid--; } return($response); }

Ahoj!

0. Jachymova data obsahuji (zaokrouhleno na tisice): 14000 prvku coz je 22000 km. 1. Z planet071003 byly vytazeny highway [1] 2. pomoci obalove zony 50m na kazdou stranu vytvorena ochranna zona. 3. Jachymova data byla rozbita na usecky (184000 prvku) 4. A pote byla dotazovana vuci ochranne zone: [2] * z toho 142000 bylo shledano jako nic nekrizici (zelena) * a 42000 bylo shledano jako konfliktni (cervena) 5. data od Jachyma byla seskupena zpet podle puvodniho deleni na puvodnich 14000 prvku

Tak mam novou versi skriptu, naucil jsem ho vzit z jednoho .shp seznam toho co uploadovat, a z druheho .shp seznam co vynechat, a zda se, ze se chova slusne... Pavel #!/usr/bin/perl use Geo::ShapeFile; # ../shpupload; cat out.osm | cstocs il2 utf8 > result-full.osm %excluded_nodes = {}; createExcludes("sil_cz12x-NE.shx"); open(OUT, ">out.osm"); print OUT "<?xml version='1.0' encoding='UTF-8'?>\n"; print OUT "<osm version='0.5' generator='shpupload'>\n"; %nodes = {}; $nodeid = -1; $way = ""; createCoasts("sil_cz12x.shx", 50.0, 14.5, 0.5); print OUT "</osm>\n"; sub createExcludes() { my ($Filename) = @_; my $shapefile = new Geo::ShapeFile($Filename); printf "%d shapes excluded\n", $shapefile->shapes(); for(1 .. $shapefile->shapes()) { $shapenum = $_; my $shape = $shapefile->get_shp_record($shapenum); $LastValid = 0; print LOG "Shape\n"; foreach $Point($shape->points()){ $Long = $Point->X(); $Lat = $Point->Y(); $Node = sprintf("lon=%f lat=%f", $Long, $Lat); $excluded_nodes{$Node} = 1; } } } sub createCoasts(){ my ($Filename, $Lat, $Long, $Size) = @_; $X1 = $Long - $Size; $X2 = $Long + $Size; $Y1 = $Lat - $Size; $Y2 = $Lat + $Size; open(LOG, ">log.txt"); print LOG "Area $X1 to $X2, $Y1 to $Y2\n"; my $shapefile = new Geo::ShapeFile($Filename); printf "%d shapes\n", $shapefile->shapes(); for(1 .. $shapefile->shapes()) { $shapenum = $_; my $shape = $shapefile->get_shp_record($shapenum); $LastValid = 0; print LOG "Shape\n"; foreach $Point($shape->points()){ $Long = $Point->X(); $Lat = $Point->Y(); $Node = sprintf("lon=%f lat=%f", $Long, $Lat); $InArea = ($Lat > $Y1 && $Lat < $Y2 && $Long > $X1 && $Long < $X2) && !$excluded_nodes{$Node}; if($InArea){ if (!$LastValid) { my %record = $shapefile->get_dbf_record($shapenum); print %record; print "\n"; $way .= "<way id='$nodeid'>\n"; $way .= " <tag k=\"created_by\" v=\"shpupload\"/>\n"; $kod = $record{'HSKOD'}; if (($kod =~ "S1") || ($kod =~ "S1T")) { $way .= " <tag k=\"highway\" v=\"primary\"/>\n"; } if (($kod =~ "S2") || ($kod =~ "S2T")) { $way .= " <tag k=\"highway\" v=\"secondary\"/>\n"; } if ($kod =~ "D[DR].*") { $way .= " <tag k=\"highway\" v=\"motorway\"/>\n"; } if ($kod =~ ".*T") { $way .= " <tag k=\"tunnel\" v=\"true\"/>\n"; $way .= " <tag k=\"layer\" v=\"-1\"/>\n"; } $way .= " <tag k=\"name\" v=\"$record{'NAZEV'}_$record{'OZNACENI'}\"/>\n"; $way .= " <tag k=\"ref\" v=\"$record{'OZNACENI'}\"/>\n"; if ($record{'POZNAMKA'}) { $way .= " <tag k=\"note\" v=\"$record{'POZNAMKA'}\"/>\n"; } $way .= " <tag k=\"source\" v=\"HELP SERVICE - REMOTE SENSING spol. s r.o. http://www.bnhelp.cz\"/>\n"; $nodeid--; } $Node = uploadNode($Lat, $Long); printf LOG "Node #%d: %f, %f\n", $Node, $Lat, $Long; $LastNode = $Node; $CountB++; $LastValid = 1; } else { if ($LastValid) { print OUT "$way</way>\n"; $way = ""; } $LastValid = 0; } $CountA++; } if ($LastValid) { print OUT "$way</way>\n"; $way = ""; } } print "Uploading $CountB of $CountA\n"; print LOG "Complete\n"; close LOG; print "Upload complete\n"; } sub uploadNode(){ ($Lat, $Long) = @_; $Tags = "<tag k=\"created_by\" v=\"shpupload\"/>"; $Tags .= "<tag k=\"source\" v=\"HELP SERVICE - REMOTE SENSING spol. s r.o. http://www.bnhelp.cz\"/>"; $Node = sprintf("lon=%f lat=%f", $Long, $Lat); if ($nodes{$Node} < 0) { $way .= " <nd ref='$nodes{$Node}' />\n"; } else { $way .= " <nd ref='$nodeid' />\n"; $nodes{$Node} = $nodeid; $Node = sprintf("<node id=\"$nodeid\" lon=\"%f\" lat=\"%f\">$Tags</node>", $Long, $Lat); print OUT " $Node\n"; $nodeid--; } return($response); } -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html

Michal Grézl michal.grezl na gmail.com #m8d374d 15. 10. 2007 18:15

Ahoj!
Takze kdy zacnes s uploadem toho zeleneho? A kde se da stahnout to
Vzorek je k dispozici v 0.5 stupne okoli 51 / 13.5. Pavel

neni to spis 50deg a 13.5 deg?

Pavel Machek pavel na ucw.cz #m6b2a96 15. 10. 2007 21:42

Ahoj!

Takze kdy zacnes s uploadem toho zeleneho? A kde se da stahnout to
Vzorek je k dispozici v 0.5 stupne okoli 51 / 13.5.
neni to spis 50deg a 13.5 deg?

Bylo to 51 / 13.5. Ted uz by melo byt uploadovano: 49 / 12.5 51 / 13.5 50 / 14.5 51 / 15.5 51 / 16.5 (Chtel jsem nejdriv neco "maleho" na otestovani, a pak jsem pridal okoli Prahy). Pavel

BH singularita na gmail.com #m7953fe 15. 10. 2007 21:53

Takze kdy zacnes s uploadem toho zeleneho? A kde se da stahnout to
Vzorek je k dispozici v 0.5 stupne okoli 51 / 13.5.
neni to spis 50deg a 13.5 deg?
Bylo to 51 / 13.5. Ted uz by melo byt uploadovano:

Tam jsem se mrknul a bylo tam kus nemecka. Kdyz jsem se mrknul do mapy co je na 51 / 13.5 ("Loc: 51°N, 13°30'E dist:10km" na mapy.cz) tak na mne vykoukly Drazdany . Kde je chyba a co se vlastne kam naimportovalo? :) No v nejhorsim z toho bude mala invaze a par ceskych silnic bude frklych v nemecku :) Martin

Pavel Machek pavel na ucw.cz #mb2edce 15. 10. 2007 22:01

Takze kdy zacnes s uploadem toho zeleneho? A kde se da stahnout to
Vzorek je k dispozici v 0.5 stupne okoli 51 / 13.5.
neni to spis 50deg a 13.5 deg?
Bylo to 51 / 13.5. Ted uz by melo byt uploadovano:
Tam jsem se mrknul a bylo tam kus nemecka. Kdyz jsem se mrknul do mapy co je na 51 / 13.5 ("Loc: 51°N, 13°30'E dist:10km" na mapy.cz) tak na mne vykoukly Drazdany . Kde je chyba a co se vlastne kam naimportovalo? :)

Trochu vic na jih by neco melo byt. Pul stupne :-).

No v nejhorsim z toho bude mala invaze a par ceskych silnic bude frklych v nemecku :)

Skript se o invazi do Nemecka pokusil, ale vcas jsem mu to zatrhnul. (Je to zajimavy, v Nemecku je zrejme osm dost aktivni, takze nam Nemci mapuji pohranici... doufam ze nebude nejaka vetsi invaze ;-). Takovy vic stredocesky priklad je 50.354 14.702. Pavel

Jachym Cepicky jachym.cepicky na gmail.com #m085e93 15. 10. 2007 22:24

jenom poznámka: OSM je v Německu dost populární. Kromě jiného je podporováno GRASS Anvender Verein (http://grass-verein.de) jachym Pavel Machek píše v Po 15. 10. 2007 v 22:01 +0200:

signature.asc

BH singularita na gmail.com #mf46421 15. 10. 2007 22:38

Trochu vic na jih by neco melo byt. Pul stupne :-).

Aha, uz tam kus vidim, cip cech tam je :)

Takovy vic stredocesky priklad je 50.354 14.702.

Koukal jsem, BTW proc je v name obsah hodnoty ref s podtrzitkem? name=_274 ref=274 Martin

Pavel Machek pavel na ucw.cz #m4388d0 15. 10. 2007 22:42

Trochu vic na jih by neco melo byt. Pul stupne :-).
Aha, uz tam kus vidim, cip cech tam je :)
Takovy vic stredocesky priklad je 50.354 14.702.
Koukal jsem, BTW proc je v name obsah hodnoty ref s podtrzitkem? name=_274 ref=274

No, v puvodni databazi bylo pole nazev... takze jsem udelal import jako $way .= " <tag k=\"name\"v=\"$record{'NAZEV'}_$record{'OZNACENI'}\"/>\n"; ...aby tam neco bylo... Pavel

hanoj enemy na mail.muni.cz #m24231a 16. 10. 2007 8:56

*** to sice jo, ale z vyse uvedeneho postupu plyne, ze ma mnozina ANO a NE prunik, takze budou vznikat prekryvy OSM / BNhelp.
Naopak, ne? Vyhazim vsechno co je v NE -> vysledek bude ze obcas zahodi o bod vic nez by mel... coz neni uplne pekny, ale zas zadna velka katastrofa...

*** jo mas pravdu, spatne jsem to precet...

*** jeste jsem vyrobil jednu verzi, za pokus to stoji
Zkusil jsem to, vysledek je v priloze... obcas to spoji zacatek a konec cesty hranou, nebo tak neco.

*** to nastane tehdy, kdyz puvodni segment pretina jina silnice a rozdeli ho vejpul, ale on byl zpatky seskupen v jeden, viz puvodni mail bod 4. Provest tyto opravy je ale asi podstatne snazsi nez dokreslovat ty vice chybejici. Taktez nerozumim proc pouzivas tag "name=_274" s informaci ktera je v ref? Na cislo komunikace je ref, na nazev ulice name. ha hanoj

Pavel Machek pavel na ucw.cz #mf33064 16. 10. 2007 11:06

Taktez nerozumim proc pouzivas tag "name=_274" s informaci ktera je v ref? Na cislo komunikace je ref, na nazev ulice name.

Ono v tom importu byl i nazev, ktery jsem dal pred "_".... jenze v ty databaze je vyplnenej velmi ridce, jestli vubec nekdy :-(. Takze name proste vyhodim. Pavel