[Proj] Newbie problem

[Gerald Evenden]

The reason for failure is no ellipsoid figure is given: needs either 
+ellps= or
+a= and +rf= (or equivalent).

Also, since you want UTM, suggest using +proj=utm this will insure all 
other
requirements are properly set up.

proj (or lproj) +proj=utm +ellps=WGS84 +zone=16 <<EOF
-88 45
EOF

yields

421184.70       4983436.77

+proj=utm will also accept +lon_0=-88 instead of +zone and will adjust 
to
the proper zone cm (-87).  For south latitudes use +south option.

On Apr 11, 2005, at 12:23 AM, Tiger Wings wrote:

> Hi
>
> I've a geotiff metdata as below:
>
> PCS = 26916 (NAD83 / UTM zone 16N)
> Projection = 16016 (UTM zone 16N)
> Projection Method: CT_TransverseMercator
>    ProjNatOriginLatGeoKey: 0.000000 (  0d 0' 0.00"N)
>    ProjNatOriginLongGeoKey: -87.000000 ( 87d 0' 0.00"W)
>    ProjScaleAtNatOriginGeoKey: 0.999600
>    ProjFalseEastingGeoKey: 500000.000000 m
>    ProjFalseNorthingGeoKey: 0.000000 m
> GCS: 4269/NAD83
> Datum: 6269/North American Datum 1983
> Ellipsoid: 7019/GRS 1980 (6378137.00,6356752.31)
> Prime Meridian: 8901/Greenwich (0.000000/  0d 0' 0.00"E)
> Projection Linear Units: 9001/metre (1.000000m)
>
> When I tried command like this:
> %>proj +proj=tmerc +lat_0=0 +lon_0=-87 k=0.9996 +x_0=500000 +y_0=0
>
> I was told:
> Rel. 4.4.9, 29 Oct 2004
> <proj>:
> projection initialization failure
> cause: major axis or radius = 0 or not given
> program abnormally terminated
>
> Did I miss something?
>
> Thanks a lot.
>
> T. W.
> _______________________________________________
> Proj mailing list
> Proj at xserve.flids.com
> http://xserve.flids.com/mailman/listinfo/proj
>
>
_____________________________________
Jerry and the low riders: Daisy Mae and Joshua.
"The whole religious complexion of the modern world is due to the
absence from Jerusalem of a lunatic asylum." Havelock Ellis, 1914