[GRASS5] V.proj error

Hamish hamish_nospam at yahoo.com
Fri Jul 22 02:52:08 EDT 2005


> 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




More information about the grass-dev mailing list