D=-dsco NAME="鄰別圖 Neighborhood map" \
-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-04-26T10:22:38+0000
## Last-Updated: 2025-05-06T02:12:47+0000
##     Update #: 64
Q=qingfuli
H=hepingli
d=~/Downloads/d/data
h=$d/$H.csv
q=$d/$Q.csv
N="addr:neighbourhood"
F='SELECT ST_ConcaveHull(ST_Collect(GEOMETRY)), $N AS name \
FROM $(notdir $(basename $<)) GROUP BY $N ORDER BY $N+0'
%.kmz:$d/%.csv
	ogr2ogr -q $@ $< -dialect SQLite -sql $F -of LIBKML $D \
	-oo X_POSSIBLE_NAMES=@lon -oo Y_POSSIBLE_NAMES=@lat
$H.vik:
$Q.vik:
%.vik:%.kmz; viking $<
## Obtaining the data from OpenStreetMap:
A=Heping Village, Beitun
#A=Qingfu Village
o=csv(::lon,::lat,"addr:neighbourhood")
Query=https://overpass-turbo.eu/?Q=[out:$o][timeout:25];\
{{geocodeArea:$A}}->.searchArea;\
node["addr:neighbourhood"~"鄰"](area.searchArea);out geom;
manually_download_data:; $${BROWSER?} '$(Query)'
#then press Run and Export and "raw data directly from Overpass API"
.PRECIOUS: %.kmz