### Makefile ---
N=Evanston Illinois address grid
DESC1=Using PLSS data.
## Author: Dan Jacobson https://www.jidanni.org/
## Copyright: https://www.gnu.org/licenses/gpl.html
## Created: 2024-02-18T00:16:28+0000
## Last-Updated: 2024-03-11T03:27:25+0000
##     Update #: 357

U = ../../../../utilities
export PERLLIB = $U
export PATH := $U:$(PATH)
include $U/m1.makefile
# For each target point use gdaltransform with the four -gcp of its governing points.
B=~/Downloads/PLSSFirstDivision.geojsonl
G=egrep
T=evanston_il
C= > $@.t && test -s $@.t && mv $@.t $@
# NO ragged Lakeshore sections please!
P=IL030410N0130E0SN(11|12|13|14|24|25)|IL030410N0140E0SN(19|30)|IL030420N0130E0SN(33|34)
#There are even "sections" for the Ouilmette Reservation!
$T.vik:
$T.csv: corners.csv corners_plus.csv pid2addr
	./pid2addr corners.csv corners_plus.csv $C
corners_plus.csv:corners.csv
	./pid2addr --false-corners $< $C
corners.csv: $B $U/PLSS_corners_4.pm # Makefile
	$G '$P' $< |\
	perl -MPLSS_corners_4 -we PLSS_corners_4::parse_seq > /tmp/$@
	test -s /tmp/$@
	sort -u -o $@.t /tmp/$@
	perl -F, -anwle '$$h{$$F[2]}++; END{for(keys %h){$(\
	)warn "Better do something about the double values for $$_" $(\
	)unless $$h{$$_}==1;}}' $@.t
	mv $@.t $@
clean:; rm *.csv *.t  #not *.kmz

# OLDER 2023 abandoned approaches:
v:
	gmt grdinfo tmp.grd
#	false
	gmt math -T1/3/1 T = -o1,1,1 |\
	gmt xyz2grd -R-126.000/-100.000/30.500/50.000 -I0.5 -Gtmp.grd
export CampusDr_SheridanRd=-87.67440 42.05014
export CentralSt_AsburyAv=-87.6893661 42.0642361
export CentralSt_GreeleyAv=-87.73243 42.06466 #NW corner of area that needs separate chunk with oblique axes!
export ChicagoAv_DavisSt=-87.6796573 42.0461939
export ChicagoAv_SheridanRd=-87.6772906 42.0511062
export Howard_StAsburyAv=-87.6901786 42.0195323
export LakeSt_SheridanRd=-87.67314 42.04398
export LakeSt_ShermanAv=-87.68192 42.04409
export UniversityPl_ShermanAv=-87.68177 42.05126
rid=123419 #OSM relation id
# Tilted downtown area.
evanston.view:
%.view:downtown.gmt east.gmt southwest.gmt
	./$* > /tmp/$*.kml
	viking --map 13 --external /tmp/$*.kml
%.gmt:%.cvs
%.gmt:%.kml
	ogr2ogr -nlt POLYGON $@ $?
downtown.csv: #counter clockwise
	/bin/echo -e "\
	$(LakeSt_ShermanAv)\n\
	$(LakeSt_SheridanRd)\n\
	$(CampusDr_SheridanRd)\n\
	$(ChicagoAv_SheridanRd)\n\
	$(UniversityPl_ShermanAv)\n\
	$(LakeSt_ShermanAv)" > $@
$(rid).kml:; ../fetch_boundary $(rid)
#Whole city
boundary_file=$rid.gmt