D=-dsco NAME="軍福路分布 Junfu roads distribution" \ -dsco DESCRIPTION="See $(subst $(HOME)/,https://www.,$(PWD))" ## Author: Dan Jacobson https://www.jidanni.org/ ## Copyright: https://www.gnu.org/licenses/gpl.html ## Created: 2025-01-27T05:25:49+0000 ## Last-Updated: 2025-02-02T06:38:20+0000 ## Update #: 436 Query=https://overpass-turbo.eu/?Q=[out:json][timeout:25];{{geocodeArea:Taichung}}->.searchArea;\ way[highway][name~"軍福"](area.searchArea);out geom;&C=24.180237;120.713608;14 Project=junfu Viewer=viking $(Project).view: %.view:%.kmz if pidof -q $(Viewer); then killall --wait --signal SIGHUP $(Viewer); fi; \ >> /tmp/nohup2.out nohup $(Viewer) $< & sleep 11 %.kmz:%.geojson ogr2ogr $@ -q -f LIBKML $< -dialect SQLite -sql \ 'SELECT Name, ST_LineMerge(ST_Collect(geometry)) FROM $(basename $<) $(\ )GROUP BY Name ORDER BY Name+0' $D %.geojson: ~/Downloads/d/data/%.geojson # A hack, as I don't have SQLite REGEXP() installed. perl -pwle 's/"name:en\D+(\d+)[^"]+"/"name":"$$1"/||s/.*"name.*//' $< > $@ manually_download_data:; $${BROWSER?} '$(Query)' .PRECIOUS: %.kmz # https://gis.stackexchange.com/questions/489971/best-way-to-put-downloaded-openstreetmap-roads-back-together-with-ogr/ # Yes we are aware that the one-way loop at the train station at the # west end of Junfu 19th Road is difficult to automatically join into # one road. # Yes I could use ST_Length to print out a report of all roads length.