## 軍福路 interpolated contour(s). See
## https://gis.stackexchange.com/questions/490024/how-to-add-one-missing-contour-line-with-gdal-ogr
## 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-04-07T00:07:15+0000
##     Update #: 318
v:h.kml; viking $^
e.csv:
	echo 'WKT\
"LINESTRING Z (120.71991 24.17197 12,120.72250 24.17178 12,120.72420 24.17151 12)"\
"LINESTRING Z (120.72420 24.17225 13,120.72251 24.17254 13,120.71991 24.17265 13)"\
"LINESTRING Z (120.71991 24.17409 15,120.72420 24.17399 15)"\
"LINESTRING Z (120.72420 24.17504 16,120.71991 24.17497 16)"'|tr -d \\\\ > $@
f.csv:e.csv Makefile
	ogr2ogr $@ $< -lco GEOMETRY=AS_WKT \
	-sql 'SELECT GEOMETRY AS WKT FROM $(basename $<)' \
	-s_srs EPSG:4326 -t_srs EPSG:3826 -dialect SQLite
g.tif:f.csv; gdal_grid -q -a linear $< $@
h.csv:g.tif; gdal_contour -q -a Name $< $@ -lco GEOMETRY=AS_WKT -i 1
h.kml:h.csv
	ogr2ogr $@ $< -f LIBKML -sql \
	'SELECT Name FROM $(basename $<)' \
	-s_srs EPSG:3826 -t_srs EPSG:4326
clean:; rm *.csv *.kml *.tif
.SECONDARY: