[Proj] How does scaling function?

[Michael Speth]

Greetings,
   We have recently discovered a new scaling problem within the HEALPix projection separate from the previous scaling error.  We have identified the problem which can be identified in the code below.

FORWARD(e_healpix_forward); /* ellipsoidal */
    double beta = auth_lat(lp.phi, P->e, 0);
    lp.phi = beta;
    return healpix_sphere(lp, P);

Here the HEALPix projection is being called on an ellipsoidal model of the earth.  The following is necessary:

Convert from Ellipsoidal to Sphere
Convert from Sphere to Plane

We need proj.4 to scale based on the sphere model not the ellipsoidal model.  How does scaling work (I have seen in proj.c that there are pre and post scaling functions).  Is it possible to set a configuration option in P so that proj.4 scales according to a Sphere and not an Ellipsoidal model?

Thanks

________________________________

Please consider the environment before printing this email
Warning: This electronic message together with any attachments is confidential. If you receive it in error: (i) you must not read, use, disclose, copy or retain it; (ii) please contact the sender immediately by reply email and then delete the emails.
The views expressed in this email may not be those of Landcare Research New Zealand Limited. http://www.landcareresearch.co.nz
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20121123/cf3b8b14/attachment.html>