D=-dsco NAME="軍福路 interpolated contours" \
-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-02-16T04:33:14+0000
## Last-Updated: 2025-02-27T02:26:07+0000
##     Update #: 133
Project=junfu_contours
Viewer=viking
#$(Project).view:
f=../junfu.kmz
o=-q $< -dialect SQLite -sql
DelaunayTriangulation.vik:
DelaunayTriangulation.kml:
input_contours.vik:
%.vik:%.kml; $(Viewer) $<
%.kml:%.csv Makefile
	ogr2ogr $@ $< -f LIBKML -sql 'SELECT Name FROM $(basename $<)' $(DS)
	cat $@
input_contours.kml:input_contours.csv Makefile
	ogr2ogr -f LIBKML $@ $o \
	"SELECT geometry FROM $(basename $<)"
	cat $@
DelaunayTriangulation.kml:input_contours.csv Makefile
	ogr2ogr -f LIBKML $@ $o \
	"SELECT ST_DelaunayTriangulation(st_collect(geometry)) AS GEOMETRY FROM $(basename $<)"
	cat $@
DelaunayTriangulation.csv:input_contours.csv
	ogr2ogr $@ -lco GEOMETRY=AS_WKT $o \
	"SELECT ST_DelaunayTriangulation(geometry), 1 AS Name FROM $(basename $<)"
DelaunayTriangulation.info:input_contours.csv
	ogrinfo $o \
	"SELECT ST_DelaunayTriangulation(geometry) FROM $(basename $<)"
# https://gis.stackexchange.com/questions/490443/why-does-st-dissolvepoints-give-me-nothing/490621
input_contours.csv:input_contours_raw.csv
	ogr2ogr $@ $< -lco GEOMETRY=AS_WKT -dialect SQLite -sql \
	'SELECT geometry AS WKT FROM "$(basename $<)"'
input_contours_raw.csv:$f #OK!
	ogr2ogr -f CSV $@ $< -lco GEOMETRY=AS_WKT -dim XYZ -zfield name -dialect SQLite \
	-sql 'SELECT name,geometry FROM "SELECT" WHERE name+0 != 17'
.SECONDARY:
clean:; rm *.csv *.kml