[GRASS5] m.proj module added to GRASS 6.1-cvs

Hi,

I have (finally) written for GRASS 6 the replacement m.proj2 from GRASS
5.
Same name and functionality, totally different code.

I've called it m.proj. It's a shell script front end to PROJ.4's cs2cs
utility, can read data from stdin or a file, and output to stdout or a
file. There are easy switches (-i and -o) for the common task of
converting to/from WGS84 LL and the current GIS map projection.

e.g., to turn a lon/lat coord into a projected vector map:
G61> echo "170.54321 -45.98765" | m.proj -i | v.in.ascii output=tmp_point

For me it is quite useful. There is probably a bug or 3 in a corner
somewhere, so please test and for now use with care.

It's in CVS now:
  http://freegis.org/cgi-bin/viewcvs.cgi/grass6/scripts/m.proj/m.proj?rev=HEAD

Here's the help page:
  http://grass.itc.it/grass61/manuals/html61_user/m.proj.html

[dev note: help page doesn't show up on the "misc modules" help list?]

I would be grateful if someone wiser in the ways of PROJ4 could have a
look at that help page and check for errors, especially in the EXAMPLES
section where it talks about using custom projection paramters. I wasn't
sure what was GRASS 5's modified proj quirks and what were PROJ4's.

+towgs84=a,b,c counts as a de facto +datum= declaration, and so datum
transform always happens if it's used, correct?

enjoy, and thanks in advance for any fact checking,
Hamish

----------------------------------------
m.proj

Description:
Convert coordinates from one projection to another (cs2cs frontend).

Usage:
m.proj [-iodv] [input=filename] [output=filename] [fs=character]
   [proj_in=string] [proj_out=string]

Flags:
  -i Use LL WGS84 as input and current location as output projection
  -o Use current location as input and LL WGS84 as output projection
  -d Output long/lat in decimal degrees or other projections with many decimal places
  -v Verbose mode (print projection parameters and filenames to stderr)

Parameters:
     input Input coordinate file (omit to read from stdin)
    output Output coordinate file (omit to send to stdout)
        fs Field separator
             default: |
   proj_in Input projection parameters (PROJ.4 style)
  proj_out Output projection parameters (PROJ.4 style)

(attachments)

m.proj.png

On Sat, Apr 22, 2006 at 09:55:32PM +1200, Hamish wrote:
...

Here's the help page:
  http://grass.itc.it/grass61/manuals/html61_user/m.proj.html

[dev note: help page doesn't show up on the "misc modules" help list?]

New modules appear on the Web each Saturday Italian morning.

At time of writing your mail it must have been updated on
the web for 1-2 hours (I don't remember exactly at which time
I have issued the cronjob):

http://www.timeanddate.com/worldclock/fixedtime.html?month=4&day=22&year=2006&hour=21&min=55&sec=0&p1=264

Markus