# dongmao.make -- analysis of mountains to my east
# Copyright       : http://www.fsf.org/copyleft/gpl.html
# Author          : Dan Jacobson -- http://jidanni.org/location/
# Created On      : Thu Nov 18 10:00:00 2004
# Last Modified On: Thu Nov 18 12:24:10 2004
# Update Count    : 23

dongmao: #2004.10
	d.mon start=PNG select=PNG
	d.frame -e
	a=taizhong y=2675560 z=846;\
	g.region rast=$$a s=2660000 n=2680000 w=236000 e=300000; \
	d.3d -l map=$$a elevation=$$a.bend color=color \
	to_coordinate=250000,$$y,$$z \
	from_coordinate=235410,$$y,$$z exaggeration=1 field=44 lines=1
	d.mon stop=PNG
#We bend to get the earth's curvature. E.g.,
$Q.bend:
	g.region rast=$Q
	xc=235410 yc=2675560;\
	r.mapcalc "$@=round($Q-((x()-$${xc?})^2+(y()-$${yc?})^2)*.0000000675)"
# Comments: the GIS image has lines at due east. This must be due to
# GRASS sweeping around a circle from an origin there...
baigu:#白姑大山
	d.frame -e
	a=taizhong y=2675560 z=846;\
	g.region rast=$$a s=2675000 n=2678000 w=258000 e=262000;\
	d.3d -l map=$$a elevation=$$a.bend color=color \
	to_coordinate=250000,$$y,$$z\
	from_coordinate=235410,$$y,$$z exaggeration=1 field=44 lines=1
#東卯山 is obvious, skipped
jiuxian:#九仙山、基隆山
	d.frame -e
	a=taizhong y=2675560 z=846;\
	g.region rast=$$a s=2673000 n=2674500 w=255000 e=257000;\
	d.3d -l map=$$a elevation=$$a.bend color=color \
	to_coordinate=250000,$$y,$$z\
	from_coordinate=235410,$$y,$$z exaggeration=1 field=44 lines=1
jilong:#基隆山 proof
	d.frame -e
	a=taizhong y=2675560 z=846;\
	g.region rast=$$a s=2673900 n=2674500 w=255000 e=256000;\
	d.3d -l map=$$a elevation=$$a.bend color=color \
	to_coordinate=255624,2674247,2938 \
	from_coordinate=235410,$$y,$$z exaggeration=1 field=6 lines=1
#Photo by Sander Borghuis