[Talk-cz] silnice 1 a 2 tridy

Pavel Machek pavel na ucw.cz
Pondělí Říjen 8 19:44:35 UTC 2007


Ahoj!

> tak silnice mám, ale taky mám dost jiné práce. Je tady někdo schopný a
> ochotný, kdo by uměl převézt silnice do formátu OSM? Momentálně je to v
> S-JTSK, Esri Shape file, ale jsem schopný dodat WGS84 v libovolném
> formátu podporovaném knihovou OGR [1]

Super! No, mam tu shp2osm skript, takze nejjednodussi pro me by asi
byl wgs84 shape file. Teda dokopat ten skript do _pouzitelny_ podoby
zrejme bude dost prace...
								Pavel
-- 
(english) http://www.livejournal.com/~pavelmachek
(cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html
------------- další část ---------------
#!/bin/bash
ID_NODE=0
> nodes.tmp
ID_SEGMENT=0
> segments.tmp

shpdump /data/gis/READ-ONLY/*/cr-shp-wgs84/linie/silnice_useky.shp | grep -v Shapefile | grep -v Bounds | grep -v nVertices | grep -v to | sed 's/(//' | sed 's/)//' | sed 's/,//g' | grep 0 | while true; do
    read LON LAT C D RING || exit

    ID_NODE=$[$ID_NODE + 1]
    echo "<node lat='$LAT' lon='$LON' id='$ID_NODE' >" >> nodes.tmp
    echo "<tag k='c' v='$C' />" >> nodes.tmp
    echo "<tag k='d' v='$D' />" >> nodes.tmp
    echo "</node>" >> nodes.tmp

    if [ .$RING = .Ring ]; then
	echo -n '|'
    else
	ID_SEGMENT=$[$ID_SEGMENT + 1]
	echo "<segment from='$[$ID_NODE - 1]' to='$ID_NODE'  id='$ID_SEGMENT'>" >> segments.tmp
	echo "</segment>" >> segments.tmp
	fi
#    echo $LON '|' $LAT '|' $C '|' $D '|' $RING
    echo -n "."
done
echo

echo "<osm version='0.3'>" > result.osm
cat nodes.tmp >> result.osm
cat segments.tmp >> result.osm
echo "</osm>" >> result.osm


Další informace o konferenci talk-cz