[Proj] WGS84 to Sphere Inconsistency Between Proj Versions
Frank Warmerdam
warmerdam at pobox.com
Thu Aug 21 23:38:50 EDT 2008
Faron Anslow wrote:
> I just ran this:
>
> cs2cs +proj=latlong +datum=WGS84 +to +proj=lcc +lat_1=50 +lat_2=50
> +lat_0=50 +lon_0=-107 +a=6371200.0000000000 +es=0.0 +f=0.0
> +towgs84=0,0,0 -r
>
> on:
> 70.933 -8.667
>
> and got:
> 2873633.37 4593659.18 -12148.43
>
> which is my original matlab answer and that from the old proj4.4. Now,
> the question is if forcing the datum shift with +towgs=0,0,0 is
> appropriate? Doesn't +towgs48 conflict with/override the spherical
> projection definition?
Faron,
I can't imagine any situation other than an effort to match old answers
where it makes sense to apply a plain lat/long on a sphere to lat/long on an
ellipsoid conversion the way this is doing.
If you want the lat/long values computed from the lcc projection based on
a sphere to treated as WGS84, then use +nadgrids=@null (or in 4.6.0 just
omit a datum specifier for the lcc projection). This is *likely* want
you want.
Best regards,
--
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up | Frank Warmerdam, warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush | Geospatial Programmer for Rent
More information about the Proj
mailing list