From 1d3508f66fddf3f8c6416e59ab9fbde405dbcf25 Mon Sep 17 00:00:00 2001
From: typebrook <typebrook@gmail.com>
Date: Sun, 17 Nov 2019 09:18:07 +0800
Subject: update

---
 scripts/osm/osm.village.nat_ref.makefile | 59 ++++++++++++++++++++++++++++++++
 1 file changed, 59 insertions(+)
 create mode 100644 scripts/osm/osm.village.nat_ref.makefile

(limited to 'scripts/osm/osm.village.nat_ref.makefile')

diff --git a/scripts/osm/osm.village.nat_ref.makefile b/scripts/osm/osm.village.nat_ref.makefile
new file mode 100644
index 0000000..3937144
--- /dev/null
+++ b/scripts/osm/osm.village.nat_ref.makefile
@@ -0,0 +1,59 @@
+data/taiwan-latest.osm.pbf:
+	mkdir -p data
+	curl -o $@ http://download.geofabrik.de/asia/taiwan-latest.osm.pbf
+
+data/village.zip:
+	mkdir -p data
+	curl -o $@ -L http://data.moi.gov.tw/MoiOD/System/DownloadFile.aspx\?DATA\=B8AF344F-B5C6-4642-AF46-1832054399CE
+
+data/VILLAGE_MOI_1081007.shp: data/village.zip
+	@if [ ! -e $@ ]; then unzip $< -d data && rm data/*xml data/*xlsx; fi
+
+village.csv: data/taiwan-latest.osm.pbf
+	ogr2ogr $@ $< \
+		-lco GEOMETRY=AS_XY \
+		-dialect sqlite \
+		-sql "SELECT osm_id, name, other_tags, ST_PointOnSurface(geometry) \
+              FROM multipolygons \
+              WHERE admin_level = '9'"
+
+village.no_nat_ref.csv: village.csv
+	grep -v nat_ref $< > $@
+
+village.with_nat_ref.csv: village.csv
+	(head -1 $<; grep nat_ref $<) |\
+	sed -r "s/\"\"\".*nat_ref\"\"=>\"\"([^\"]+).*\"\"\"/\1/g" |\
+	sed '1s/other_tags/nat_ref/'> $@
+
+matched.csv: data/VILLAGE_MOI_1081007.shp village.no_nat_ref.csv
+	ogr2ogr $@ $(word 2,$^) \
+		-oo X_POSSIBLE_NAMES=X -oo Y_POSSIBLE_NAMES=Y \
+		-dialect sqlite \
+		-sql "SELECT osm.osm_id, gov.* \
+		      FROM 'village.no_nat_ref' osm, '$<'.VILLAGE_MOI_1081007 gov \
+		      WHERE osm.name = gov.VILLNAME AND Intersects(gov.geometry, osm.geometry)"
+
+matched.by_ref.csv: data/VILLAGE_MOI_1081007.shp village.with_nat_ref.csv
+	ogr2ogr $@ $(word 2,$^) \
+		-oo X_POSSIBLE_NAMES=X -oo Y_POSSIBLE_NAMES=Y \
+		-dialect sqlite \
+		-sql "SELECT osm.osm_id, gov.* \
+		      FROM 'village.with_nat_ref' osm, '$<'.VILLAGE_MOI_1081007 gov \
+		      WHERE osm.nat_ref = gov.VILLCODE"
+
+confilct.list: matched.csv
+	cat $< | cut -d',' -f2 | sort | uniq -d | xargs -I {} grep {} $<
+
+origin.osm: matched.csv data/taiwan-latest.osm.pbf
+	cat $< |\
+	sed 1d |\
+	cut -d'"' -f2 |\
+	sed -nr 's/.*/r\0/p' |\
+	osmium getid $(word 2,$^) --id-file - --output-format=osm > $@
+
+change.list: matched.csv
+	cat $< |\
+	sed 1d |\
+	awk -F',' '{print "relation", $$1, "nat_ref", $$2}' > $@
+
+# sed -i -r 's/([0-9]+ +)(.+)$/\1"name:en" "\2"/' diff.eng
-- 
cgit v1.2.3-70-g09d2