v.proj

Hello GRASS Users,

Thanks again to J. Evenden and C. Weiss for their help with my SPCS questions.

I have succeeded in setting up a location/mapset in the NYSPCS, and now
I want to bring in a vector map, done in the UTM system in another location/
mapset, to this new one using v.proj. I answered all the prompts from
v.proj, i.e. name, database, location,...etc. pertaining to the map I want to
bring in, but I get an error message:

  major axis, radius = 0, or not spec'ed
  cannot initialize pj

The PROJ_INFO and PROJ_UNITS files in my new location were created by GRASS when
I set it up. Those in the location where I want to grab the UTM map seem okay,
they list the ellipsoid, etc. I'm not clear what is missing or where it should be.
I have not been able to locate much documentation on what the PROJ_INFO files should
look like.

All help is greatly appreciated.

Gary

GaryO@hydroqual.com

Ouch, we solved this one before and I'm trying to dredge it from
memory. I *believe* the problem is that the PROJ_INFO was created
with an older version of GRASS than is needed by v.proj. We have
been going in the other direction (SP to UTM) and we needed to
create a brand new UTM location under GRASS4.1. This made the
later version of the PROJ_INFO, which I then copied back over the
old one in our main location. I think it is the "a" and "es" that
were missing.

Our current PROJ_INFO looks like:
  name: UTM
  proj: utm
  ellps: clark66
  a: 6378206.4000000004
  es: 0.0067686580
  zone: 10

and our State Plane (California) looks like:
  name: State Plane
  proj: lcc
  a: 0.63782064e+07
  es: 0.6768657997291094e-02
  x_0: 0.6096012192024384e+06
  y_0: 0
  lon_0: 122dw
  lat_0: 37d40'n
  lat_1: 39d50'n
  lat_2: 38d20'n

I believe that these work for us.

Sue Huse
REGIS
UC Berkeley