⚠️HamishBowman: /* Creating a Voronoi diagram */

2010-10-08T08:05:55Z

Creating a Voronoi diagram

← Older revision		Revision as of 08:05, 8 October 2010
Line 62:		Line 62:

	* Overlay the Voronoi diagram with the buffer:		* Overlay the Voronoi diagram with the buffer:
	v.overlay ain=local_circuits_nzgd2k_voronoi bin=local_circuits_nzgd2k_buf150km out=local_circuits_nzgd2k_voronoi_buf150km op=and		v.overlay ain=local_circuits_nzgd2k_voronoi bin=local_circuits_nzgd2k_buf150km \
			out=local_circuits_nzgd2k_voronoi_buf150km op=and

	== Rendering ==		== Rendering ==

⚠️HamishBowman: /* Prep work */

2010-10-08T08:05:05Z

Prep work

← Older revision		Revision as of 08:05, 8 October 2010
Line 43:		Line 43:
	v.in.ascii in=nzgd2k_epsg_xy.dat \		v.in.ascii in=nzgd2k_epsg_xy.dat \
	out=local_circuits_nzgd2k x=6 y=5 \		out=local_circuits_nzgd2k x=6 y=5 \
	columns='datum varchar(8), projection varchar(40), epsg_code integer, proj4_str varchar(255), lat_0 double, lon_0 double' ~~--o~~		columns='datum varchar(8), projection varchar(40), epsg_code integer, proj4_str varchar(255), lat_0 double, lon_0 double'

	== Creating a Voronoi diagram ==		== Creating a Voronoi diagram ==

⚠️HamishBowman: /* Rendering */

2010-10-08T08:04:07Z

Rendering

← Older revision		Revision as of 08:04, 8 October 2010
Line 71:		Line 71:

	* ...		* ...
	g.manual htmlmapdriver		g.manual {{cmd\|htmlmapdriver}}

⚠️HamishBowman: init

2010-10-08T08:02:36Z

init

New page

This page will describe an example of using GRASS 6's HTMLMAP driver to create a clickable web page based on vector data.

Specifically, we'll create a page with a map showing all of the local Transverse Mercator EPSG codes, and clicking on the map will show you the code and its parameters.

There's a little prep work needed to extract the EPSG codes into a GRASS vector map. The details of doing this are not important to the overall example.

== Prep work ==

* Parse the epsg file for NZGD2000 projections and make into a flat file, with lat_0 and lon_0 extracted:

grep -A1 NZGD2000 /usr/share/proj/epsg | grep -v '^--$' | \
grep -B1 '+lat_0' | grep -v '^--$' | \
sed -e 's| / $[^ ]*$ $[^ ]*$| / \1^\2|' \
-e 's| / $[^ ]*$ $[^ ]*$| / \1^\2|' \
-e 's| / $[^ ]*$ $[^ ]*$| / \1^\2|' \
-e 's| / $[^ ]*$ $[^ ]*$| / \1^\2|' | \
tr '\n' ' ' | tr '#' '\n' | \
sed -e 's/^ //' -e 's/ <> $//' -e 's/ \/ / /' \
-e 's/ /|/' -e 's/ /|/' -e 's/ /|/' \
-e 's/|</|/' -e 's/>|/|/' -e 's/\^/ /g' \
> nzgd2k_epsg.dat

cat nzgd2k_epsg.dat | sed -e 's/.*+lat_0=//' -e 's/ .*//' \
> lat0s.dat

cat nzgd2k_epsg.dat | sed -e 's/.*+lon_0=//' -e 's/ .*//' \
> lon0s.dat

paste -d'|' nzgd2k_epsg.dat lat0s.dat lon0s.dat > nzgd2k_epsg_xy.dat

* Some of these projections have an undefined (or 0 value) for +lat_0. To make things easy we can specify a nearby value by hand:
cut -d'|' -f2,5 nzgd2k_epsg_xy.dat | grep '|0$' | cut -f1 -d'|'

Auckland Islands TM 2000 (50.75 degrees S)
Campbell Island TM 2000 (52.5 degrees S)
Antipodes Islands TM 2000 (49.75 degrees S)
Raoul Island TM 2000 (29.25 degrees S)
Chatham Island Circuit 2000 (44 degrees S)

* Then we can import it into a GRASS lat/lon WGS84 location:
v.in.ascii in=nzgd2k_epsg_xy.dat \
out=local_circuits_nzgd2k x=6 y=5 \
columns='datum varchar(8), projection varchar(40), epsg_code integer, proj4_str varchar(255), lat_0 double, lon_0 double' --o

== Creating a Voronoi diagram ==

* Create a LCC location to work in (NZCS2000, epsg code <3851>)
* Use the {{cmd|v.proj}} module to bring in the projection center points:
v.proj in=local_circuits_nzgd2k loc=ll_wgs84 mapset=nz

* Set the bounding region for the diagram with the {{cmd|g.region}} module:
g.region north=9380000 south=4600000 east=5960000 west=470000 res=10000

* Run {{cmd|v.voronoi}}:
v.voronoi in=local_circuits_nzgd2k out=local_circuits_nzgd2k_voronoi

* Create a 150km buffer around the points:
v.buffer in=local_circuits_nzgd2k out=local_circuits_nzgd2k_buf150km dist=150000
: ''As these are all tmerc, maybe we should use the v.buffer angle= option to allow further use north/south?''

* Overlay the Voronoi diagram with the buffer:
v.overlay ain=local_circuits_nzgd2k_voronoi bin=local_circuits_nzgd2k_buf150km out=local_circuits_nzgd2k_voronoi_buf150km op=and

== Rendering ==

* Grab a 1:10 million scale vector coastline of the world from http://www.naturalearthdata.com
** Import the shapefile into a suitable location with {{cmd|v.in.ogr}}.
** Then bring it over into the LCC location with {{cmd|v.proj}}.

* ...
g.manual htmlmapdriver

d.save -o > render_commands.sh
d.mon start=HTMLMAP
. render_commands.sh
d.mon stop=HTMLMAP

* TBC ...

HTMLMAP driver - Revision history

⚠️HamishBowman: /* Creating a Voronoi diagram */

⚠️HamishBowman: /* Prep work */

⚠️HamishBowman: /* Rendering */

⚠️HamishBowman: init