abstreet/import.sh

99 lines
3.1 KiB
Bash
Raw Normal View History

#!/bin/bash
set -e
# First prepare input.
2018-10-03 17:47:59 +03:00
function get_if_needed {
if [ ! -f $2 ]; then
wget $1 -O $2;
fi
}
2018-10-20 05:20:42 +03:00
mkdir -p data/input data/raw_maps
# TODO refactor a variant for .zips?
if [ ! -f data/input/N47W122.hgt ]; then
get_if_needed \
https://dds.cr.usgs.gov/srtm/version2_1/SRTM1/Region_01/N47W122.hgt.zip \
data/input/N47W122.hgt.zip;
unzip -d data/input data/input/N47W122.hgt.zip;
rm -f data/input/N47W122.hgt.zip;
fi
if [ ! -d data/input/google_transit_2018_18_08/ ]; then
get_if_needed \
https://metro.kingcounty.gov/GTFS/google_transit_2018_18_08.zip \
data/input/google_transit_2018_18_08.zip;
unzip -d data/input/google_transit_2018_18_08 data/input/google_transit_2018_18_08.zip;
rm -f data/input/google_transit_2018_18_08.zip;
fi
if [ ! -f data/input/traffic_signals.kml ]; then
# From https://data.seattle.gov/Transportation/Traffic-Signals/nr6x-wnd5
get_if_needed \
http://data-seattlecitygis.opendata.arcgis.com/datasets/ff97a6eb8ac84356beea09138c6e1ec3_0.kml \
data/input/traffic_signals.kml;
fi
2018-10-15 05:06:57 +03:00
if [ ! -f data/input/neighborhoods.geojson ]; then
# https://data.seattle.gov/dataset/Neighborhoods/2mbt-aqqx in GeoJSON, not SHP
get_if_needed \
https://github.com/seattleio/seattle-boundaries-data/raw/master/data/neighborhoods.geojson \
data/input/neighborhoods.geojson;
2018-10-15 05:06:57 +03:00
fi
if [ ! -f data/input/Seattle.osm ]; then
get_if_needed \
http://download.bbbike.org/osm/bbbike/Seattle/Seattle.osm.gz \
data/input/Seattle.osm.gz;
gunzip data/input/Seattle.osm.gz;
fi
for poly in `ls data/polygons/`; do
name=`basename -s .poly $poly`;
if [ ! -f data/input/$name.osm ]; then
osmosis \
--read-xml enableDateParsing=no file=data/input/Seattle.osm \
--bounding-polygon file=data/polygons/$name.poly completeWays=true \
--write-xml data/input/$name.osm
fi
done
if [ ! -f data/shapes/parcels ]; then
# From https://gis-kingcounty.opendata.arcgis.com/datasets/king-county-parcels--parcel-area/geoservice
# TODO This isn't a direct link
#get_if_needed https://opendata.arcgis.com/datasets/8058a0c540434dadbe3ea0ade6565143_439.kml data/input/King_County_Parcels__parcel_area.kml;
cd kml
time cargo run --release -- \
--input=../data/input/King_County_Parcels__parcel_area.kml \
--output=../data/shapes/parcels
cd ..
fi
if [ ! -f data/shapes/blockface ]; then
# From http://data-seattlecitygis.opendata.arcgis.com/datasets/blockface
get_if_needed https://opendata.arcgis.com/datasets/a1458ad1abca41869b81f7c0db0cd777_0.kml data/input/blockface.kml;
cd kml
time cargo run --release -- \
--input=../data/input/blockface.kml \
--output=../data/shapes/blockface
cd ..
fi
cd convert_osm
2018-10-31 01:00:32 +03:00
for poly in `ls ../data/polygons/`; do
name=`basename -s .poly $poly`;
cargo run --release -- \
--osm=../data/input/$name.osm \
--elevation=../data/input/N47W122.hgt \
--traffic_signals=../data/input/traffic_signals.kml \
--parcels=../data/shapes/parcels \
--parking_shapes=../data/shapes/blockface \
--gtfs=../data/input/google_transit_2018_18_08 \
--neighborhoods=../data/input/neighborhoods.geojson \
--output=../data/raw_maps/$name.abst
done