[GRASS5] Datum transformation with r.proj

Joel Peter William Pitt pittj2 at lincoln.ac.nz
Tue Oct 5 21:18:07 EDT 2004


Hi,

I was curious about whether r.proj and v.proj automatically do datum 
transformation
if it seems necessary. I followed the output and it seems to be doing it:

Input Projection Parameters: +proj=nzmg +lat_0=-41.0000000000 
+lon_0=173.0000000000 +x_0=2510000.0000000000 +y_0=6023150.0000000000 
+a=6378388 +rf=297 +no_defs +nadgrids=nzgd2kgrid0005.gsb
Input Unit Factor: 1

Output Projection Parameters: +proj=tmerc +lat_0=0.0000000000 
+lon_0=173.0000000000 +k_0=0.9996000000 +x_0=1600000.0000000000 
+y_0=10000000.0000000000 +a=6378137 +rf=298.257222101 +no_defs 
+towgs84=0.000,0.000,0.000
Output Unit Factor: 1
Input:
Cols:   45 (10200)
Rows:   50 (14900)
North: 5444600.000000 (6790000.000000)
South: 5439600.000000 (5300000.000000)
West:  2006400.000000 (1990000.000000)
East:  2010900.000000 (3010000.000000)
ew-res: 100.000000
ns-res: 100.000000

Output:
Cols:   41 (41)
Rows:   46 (46)
North: 4882159.892804 (4882159.892804)
South: 4877561.531759 (4877561.531759)
West:  1096524.283564 (1096524.283564)
East:  1100623.252896 (1100623.252896)
ew-res: 99.974862
ns-res: 99.964371
Allocating memory and reading input map...  100%
Projecting...  100%

-----------
Both locations have their projections and datums set correctly using
entries already in the etc/datumtransform.table file.

However, when I check the my maps, they are still off centre by the
amount expected by no datum transformation....

(Note, perhaps the man page for r.proj should be changed to
mention that it automatically does datum transformation,
rather than "supporting" it. Since I wasn't sure if there were
meant to be some missing parameters or what....)

Any help appreciated, thanks
Joel




More information about the grass-dev mailing list