[GRASS5] V.proj error

This is an error that I haven’t seen come across the list. I’m trying to reproject vectors from a latlon nad27 location into a wgs84 utm location. I’m using a 9 July 2005 MAC OSX binary. Any thoughts?

v.proj input=chev_heb_roads location=mrss_latlon mapset=Apache-Sitgreaves output=chevheber_roads

Input Projection Parameters: +proj=latlong +a=6378206.4 +rf=294.9786982 +no_defs +nadgrids=/Applications/Grass/grass61cvs.app/Contents/Resources/grass-6.1.cvs/etc/nad/conus
Input Unit Factor: 1

Output Projection Parameters: +proj=utm +zone=12 +a=6378137 +rf=298.257223563 +no_defs +towgs84=0.000,0.000,0.000
Output Unit Factor: 1

Creating vector file…
pj_transform() failed
cause: failed to load NAD27-83 correction file
Error in pj_do_transform


C. Michael Barton, Professor of Anthropology
School of Human Evolution and Social Change
PO Box 872402
Arizona State University
Tempe, AZ 85287-2402
USA

Phone: 480-965-6262
Fax: 480-965-7671
www: <www.public.asu.edu/~cmbarton>

This is an error that I haven¹t seen come across the list. I¹m trying
to reproject vectors from a latlon nad27 location into a wgs84 utm
location. I¹m using a 9 July 2005 MAC OSX binary. Any thoughts?

v.proj input=chev_heb_roads location=mrss_latlon
mapset=Apache-Sitgreaves output=chevheber_roads

Input Projection Parameters: +proj=latlong +a=6378206.4
+rf=294.9786982 +no_defs
+nadgrids=/Applications/Grass/grass61cvs.app/Contents/Resources/grass
-6.1.cv s/etc/nad/conus
Input Unit Factor: 1

Output Projection Parameters: +proj=utm +zone=12 +a=6378137
+rf=298.257223563 +no_defs +towgs84=0.000,0.000,0.000
Output Unit Factor: 1

Creating vector file...
pj_transform() failed
cause: failed to load NAD27-83 correction file
Error in pj_do_transform

this has to do with Paul's recent fix to make `g.proj -j` display the
full path name to the nadgrid file. Does the
/Applications/Grass/grass61cvs.app/Contents/Resources/grass-6.1.cvs/etc/nad/conus

file exist? I saw the same problem with r(?).proj when I changed the
string to $GISBASE locally, but it has been working fine for me using
the full path (in linux).

try using cs2cs from the command line with some test points:

IN_PROJ="`g.proj -jf`"
OUT_PROJ="+proj=longlat +towgs84=0.000,0.000,0.000"

echo "lon lat" | cs2cs -f %.7f $IN_PROJ +to $OUT_PROJ

Hamish

Hamish,

The conus file does exist in the following location

/Applications/Grass/grass61cvs.app/Contents/Resources/grass-6.1.cvs/etc/nad/
conus

Michael

On 7/21/05 11:52 PM, "Hamish" <hamish_nospam@yahoo.com> wrote:

This is an error that I haven¹t seen come across the list. I¹m trying
to reproject vectors from a latlon nad27 location into a wgs84 utm
location. I¹m using a 9 July 2005 MAC OSX binary. Any thoughts?

v.proj input=chev_heb_roads location=mrss_latlon
mapset=Apache-Sitgreaves output=chevheber_roads

Input Projection Parameters: +proj=latlong +a=6378206.4
+rf=294.9786982 +no_defs
+nadgrids=/Applications/Grass/grass61cvs.app/Contents/Resources/grass
-6.1.cv s/etc/nad/conus
Input Unit Factor: 1

Output Projection Parameters: +proj=utm +zone=12 +a=6378137
+rf=298.257223563 +no_defs +towgs84=0.000,0.000,0.000
Output Unit Factor: 1

Creating vector file...
pj_transform() failed
cause: failed to load NAD27-83 correction file
Error in pj_do_transform

this has to do with Paul's recent fix to make `g.proj -j` display the
full path name to the nadgrid file. Does the
/Applications/Grass/grass61cvs.app/Contents/Resources/grass-6.1.cvs/etc/nad/co
nus

file exist? I saw the same problem with r(?).proj when I changed the
string to $GISBASE locally, but it has been working fine for me using
the full path (in linux).

try using cs2cs from the command line with some test points:

IN_PROJ="`g.proj -jf`"
OUT_PROJ="+proj=longlat +towgs84=0.000,0.000,0.000"

echo "lon lat" | cs2cs -f %.7f $IN_PROJ +to $OUT_PROJ

Hamish

____________________
C. Michael Barton, Professor of Anthropology
School of Human Evolution and Social Change
PO Box 872402
Arizona State University
Tempe, AZ 85287-2402
USA

Phone: 480-965-6262
Fax: 480-965-7671
www: <www.public.asu.edu/~cmbarton>

On 7/21/05 11:52 PM, "Hamish" <hamish_nospam@yahoo.com> wrote:

Input Projection Parameters: +proj=latlong +a=6378206.4
+rf=294.9786982 +no_defs
+nadgrids=/Applications/Grass/grass61cvs.app/Contents/Resources/grass
-6.1.cv s/etc/nad/conus
Input Unit Factor: 1

Output Projection Parameters: +proj=utm +zone=12 +a=6378137
+rf=298.257223563 +no_defs +towgs84=0.000,0.000,0.000
Output Unit Factor: 1

Creating vector file...
pj_transform() failed
cause: failed to load NAD27-83 correction file
Error in pj_do_transform

this has to do with Paul's recent fix to make `g.proj -j` display the
full path name to the nadgrid file. Does the
/Applications/Grass/grass61cvs.app/Contents/Resources/grass-6.1.cvs/etc/nad/co
nus

file exist? I saw the same problem with r(?).proj when I changed the
string to $GISBASE locally, but it has been working fine for me using
the full path (in linux).

try using cs2cs from the command line with some test points:

IN_PROJ="`g.proj -jf`"
OUT_PROJ="+proj=longlat +towgs84=0.000,0.000,0.000"

echo "lon lat" | cs2cs -f %.7f $IN_PROJ +to $OUT_PROJ

Hello Michael, did you do this? I can't reproduce the problem using spearfish. If you can reproduce it using cs2cs it might make things clearer. E.g. does the following work:

echo "-103.7 44.4" | \
cs2cs +proj=longlat +a=6378206.4 +rf=294.9786982 \
+nadgrids=/Applications/Grass/grass61cvs.app/Cont\
ents/Resources/grass-6.1.cvs/etc/nad/conus \
+to +proj=utm +zone=12 +a=6378137 +rf=298.257223563 \
+towgs84=0.000,0.000,0.000

(should give 1081318.90 4942279.92 0.00)

Paul