-
Notifications
You must be signed in to change notification settings - Fork 2
/
prepare
21 lines (17 loc) · 861 Bytes
/
prepare
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
#!/bin/bash
rm -rvf japan
mkdir -p build japan
if [ ! -f build/polbnda_jpn.shp ]; then
curl -o build/gm-jpn-all_u_2_2.zip 'http://www1.gsi.go.jp/geowww/globalmap-gsi/download/data/gm-japan/gm-jpn-all_u_2_2'
unzip -od build -j build/gm-jpn-all_u_2_2.zip gm-jpn-all_u_2_2/polbnda_jpn.shp gm-jpn-all_u_2_2/polbnda_jpn.dbf
chmod a-x build/polbnda_jpn.*
fi
geo2topo -q 1e5 -n municipalities=<( \
shp2json -n build/polbnda_jpn.shp \
| ndjson-map '(d.id = d.properties.adm_code, delete d.properties, d)' \
| ndjson-filter '(d.id != "UNK")' \
| geoproject -n 'd3.geoAzimuthalEqualArea().center([138.5, 36 - 18 / 60]).rotate([-138.5, -36 + 18 / 60]).scale([1700]).translate([2500, -1930])') \
| toposimplify -f -p 0.001 \
| topomerge prefectures=municipalities -k 'd.id.slice(0, 2)' \
| topomerge country=prefectures \
> japan/japan.json