DS =\
-dsco NAME="Dr. A. Santos Ave., Parañaque, Philippines addresses" \
-dsco DESCRIPTION="See $(subst $(HOME)/,https://www.,$(PWD))"
### Makefile ---
## Author: Dan Jacobson https://www.jidanni.org/
## Copyright: https://www.gnu.org/licenses/gpl.html
## Created: 2025-02-24T02:00:52+0000
## Last-Updated: 2025-04-10T19:17:20+0000
##     Update #: 175
##
# https://overpass-turbo.eu/?Q=%7B%7BgeocodeArea%3APara%C3%B1aque%7D%7D-%3E.q%3B%0Anw%5B%22addr%3Ahousenumber%22%5D%5B%22addr%3Astreet%22%7E%22Santos%20Ave%22%5D%28area.q%29%3Bout%20geom%3B
d=~/Downloads
f=$d/santos.geojson
P=dr_a_santos_ave
q=slight_angle
v:$q.vik #P.vik
%.vik:%.kmz; viking $<
#We got house outlines(w:way), and address points (n:node), from our Overpass Turbo OSM query.
#Let's make them all into points for easier mapping:
raw.csv:$f
	ogr2ogr $@ $< -lco GEOMETRY=AS_WKT -lco STRING_QUOTING=IF_NEEDED -dialect SQLite -q -sql \
	'SELECT "addr:housenumber" AS name, ST_Centroid(GEOMETRY) AS "WKT" FROM $(notdir $(basename $<))'
filtered.csv:raw.csv
	perl -F, -anwle 'if(/WKT/){print;next}; if($$F[-1] =~ /8\d{3}/){print "$$F[0],$$&"}' $< > $@
$P.kmz:filtered.csv
	ogr2ogr $(DS) -f LIBKML $@ $< -q -sql \
	'SELECT name FROM $(basename $<)'
# Test to see if A. Santos Ave. uses the same grid as Malate:
$q.optfile:raw.csv
	ogr2ogr /vsistdout/ -f CSV $< -lco STRING_QUOTING=IF_NEEDED \
	-sql 'SELECT X(GEOMETRY) AS P, Name AS Q, X(GEOMETRY), Y(GEOMETRY) FROM $(basename $<)' \
	-dialect SQLite | egrep 8165\|842 | sed 'y/,/ /;s/^/-gcp /' > $@
	set 1928 14.570815 120.986704; echo $$3 $$1 $$3 $$2|sed 's/^/-gcp /' >> $@ #Bocobo
# Drape contours over the wedge. Our own version of "gdal_create -outsize 2 2".
# (Makes four points per line, we actually only want the middle two, but will get clipped anyway):
K=gdal_contour -q -a Name $< $@ -lco GEOMETRY=AS_WKT
%.csv:%.xyz; $K -i 100
%.single_axis.csv:%.xyz; $K -fl 5000
clean:; rm *.xyz *.kml *.csv *.optfile
W=120.95
E=121.05
S=8500
N=1500
$q.%.xyz: Makefile; perl -wle "for([$W,$N],[$E,$N],[$W,$S],[$E,$S]){print qq(@\$$_ @\$$_[$*])}" > $@
$q.%.kml:$q.%.csv $q.optfile
	ogr2ogr -f LIBKML $@ $< --optfile $q.optfile \
	-sql 'SELECT Name FROM "$(basename $<)"' #-clipdst $B.csv
$q.kmz: $(guts) $q.1.kml; ogrmerge -f LIBKML -overwrite_ds -o $@ $^
.SECONDARY:
.PRECIOUS: %.kmz