[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