## Makefile --- Las Animas Co., Colorado addressing and road numbering grid.
## Author: Dan Jacobson https://www.jidanni.org/
## Copyright: https://www.gnu.org/licenses/gpl.html
## Created: 2023-11-22T14:00:05+0000
## Last-Updated: 2023-12-31T10:02:52+0000
##     Update #: 84
## Now using PLSS as basis.
U = ../../../../../utilities
export PERLLIB = $U
export PATH := $U:$(PATH)
N=_POSSIBLE_NAMES=field_
KO = 	-f LIBKML -dsco NAME="Las Animas Co. Colorado address and road name grid" \
	-dsco DESCRIPTION="Using BLM PLSS data. See $(subst $(HOME)/,https://www.,$(PWD))"
CO=	-oo HEADERS=NO -oo KEEP_GEOM_COLUMNS=NO -oo X$N1 -oo Y$N2 \
	-sql 'SELECT field_3 AS Name FROM "$(basename $<)"'
l.vik:
l.csv:pid2addr n.csv; ./$^ > $@
s=sqlite3
C= ~/Downloads/BLM_CO_PLSS_Points_Corners.csv
n.csv:c.sqlite # [04] get half miles to show nature of .. 95, 95.5, 97, 97,5..
	perl -wle 'my %h=(X=>"CO060310S0(5[^0]|6.)0W0_[1-6][04]0100", $(\
	)Y=>"CO060(2[6-9]|3.)0S0(53|61)0W0_100[1-6][04]0"); $(\
	)for(sort keys %h){print "SELECT XCOORD,YCOORD,POINTID FROM t WHERE POINTID REGEXP $(\
	)\"$$h{$$_}\" ORDER BY $${_}COORD + 0;"}'| $s -csv $< | grep -v CO0602[67]0S0530W0 > $@
c.sqlite:$C #Big file...
	echo .import --csv $< t | $s $@
%.kmz %.kml:%.csv; ogr2ogr $@ $? $(KO) $(CO)
%.vik:%.kmz
	cd /tmp && nohup viking -e $(PWD)/$? &
	sleep 19
.PRECIOUS: %.kmz
#Wow, two furlongs distance! (16 points per square mile), etc.:
ppc:$C
	perl -F, -anwle '$$k{$$F[5]}++; END{for(sort keys %k){next unless $$k{$$_}>=50; $(\
	)printf "%s %s\n", $$_, $$k{$$_};}}' $?