|
@@ -34,7 +34,7 @@ Projections are defined by precise mathematical relations, so the method
|
|
|
of projecting coordinates from a geographic reference frame
|
|
|
(latitude-longitude) into a projected cartesian reference frame (eg metres)
|
|
|
is governed by these equations. Inverse projections can also be achieved.
|
|
|
-The public-domain Unix software package <i>PROJ</i> [1] has been designed to
|
|
|
+The public-domain Unix software package <i>PROJ.4</i> [1] has been designed to
|
|
|
perform these transformations, and the user's manual contains a detailed
|
|
|
description of over 100 useful projections. This also includes a
|
|
|
programmers library of the projection methods to support other software
|
|
@@ -69,7 +69,9 @@ prior to import or after export, using software such as <i>PROJ.4</i>'s
|
|
|
<i><a href="http://proj.maptools.org/">cs2cs</a></i> [1]. This is an easy
|
|
|
method for converting an ASCII file containing a list of coordinate points,
|
|
|
since there is no topology to be preserved and <i>cs2cs</i> can be used to
|
|
|
-process simple lists using a one-line command.
|
|
|
+process simple lists using a one-line command. The <em>m.proj</em> module
|
|
|
+provides a handy front end to <tt>cs2cs</tt>.
|
|
|
+
|
|
|
<p>
|
|
|
The format of files containing vector maps with <b>lines</b> and <b>arcs</b> is
|
|
|
generally more complex, as parts of the data stored in the files will describe
|
|
@@ -180,7 +182,7 @@ then use that to set the region bounds before performing the reprojection:
|
|
|
|
|
|
<div class="code"><pre>
|
|
|
# calculate where output map will be
|
|
|
-GRASS> r.proj in=elevation location=ll_wgs84 mapset=user1 -p
|
|
|
+GRASS> r.proj input=elevation location=ll_wgs84 mapset=user1 -p
|
|
|
Source cols: 8162
|
|
|
Source rows: 12277
|
|
|
Local north: -4265502.30382993
|
|
@@ -189,7 +191,7 @@ Local west: 14271663.19157564
|
|
|
Local east: 14409956.2693866
|
|
|
|
|
|
# same calculation, but in a form which can be cut and pasted into a g.region call
|
|
|
-GRASS> r.proj in=elevation location=ll_wgs84 mapset=user1 -g
|
|
|
+GRASS> r.proj input=elevation location=ll_wgs84 mapset=user1 -g
|
|
|
n=-4265502.30382993 s=-4473453.15255565 w=14271663.19157564 e=14409956.2693866 rows=12277 cols=8162
|
|
|
|
|
|
GRASS> g.region n=-4265502.30382993 s=-4473453.15255565 \
|
|
@@ -225,7 +227,7 @@ cols: 8136
|
|
|
cells: 99535824
|
|
|
|
|
|
# finally, perform the reprojection
|
|
|
-GRASS> r.proj in=elevation location=ll_wgs84 mapset=user1 memory=800
|
|
|
+GRASS> r.proj input=elevation location=ll_wgs84 mapset=user1 memory=800
|
|
|
</pre></div>
|
|
|
|
|
|
|
|
@@ -284,6 +286,7 @@ Springer-Verlag, Berlin, 2nd edition.
|
|
|
<a href="g.proj.html">g.proj</a>,
|
|
|
<a href="g.setproj.html">g.setproj</a>,
|
|
|
<a href="i.rectify.html">i.rectify</a>,
|
|
|
+<a href="m.proj.html">m.proj</a>,
|
|
|
<a href="r.support.html">r.support</a>,
|
|
|
<a href="r.stats.html">r.stats</a>,
|
|
|
<a href="v.proj.html">v.proj</a>,
|