[Proj] Convert a planetographic latitude to a planetocentric latitude

[Even Rouault]

On samedi 18 août 2018 20:11:24 CEST Kristian Evers wrote:
> Jean-Christophe,
> 
> You can do that transformation with PROJ, yes. I can see why you get
> confused though. It is not really apparent but to do this type of
> transformation cs2cs needs some hint that it is doing a change in reference
> frame (which this essentially is). You can do that by adding +towgs84=0,0,0
> to both CRS1 and CRS2. When I do that I get the following:
 
> $ echo 0 80 0 | cs2cs +proj=latlong +R=3396190 +towgs84=0,0,0 +to
> +proj=latlong +a=3396190 +b=3376200 +towgs84=0,0,0
 0dE     80d6'51.824"N
> 19392.342
> 
> It is not unreasonable to expect cs2cs to figure this out on it’s own. I’ll
> give it some thoughts the next couple of days and see if I can come up with
> a simple solution that handles this particular case.

Kristian,

cs2cs *used* to do this implicit +towgs84=0,0,0 in versions earlier than 4.6, 
and in 4.6 this was changed:

4.6.0 Release Notes
-------------------

 o MAJOR: Rework pj_transform() to avoid applying ellipsoid to ellipsoid
   transformations as a datum shift when no datum info is available.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com