diff options
| author | Thomas Knudsen <thokn@sdfe.dk> | 2017-11-16 20:01:39 +0100 |
|---|---|---|
| committer | Thomas Knudsen <thokn@sdfe.dk> | 2017-11-16 20:01:39 +0100 |
| commit | e6b218cf4fce2917cdf529988dc9ef9d296c17c1 (patch) | |
| tree | ce2bca5473bd1174d998c28886c754225c11f4dc /src | |
| parent | 93685def8f4b23d220d04bc658a0ff27e93149e6 (diff) | |
| download | PROJ-e6b218cf4fce2917cdf529988dc9ef9d296c17c1.tar.gz PROJ-e6b218cf4fce2917cdf529988dc9ef9d296c17c1.zip | |
Switch pj_factors to use proj_geoc_lat
Diffstat (limited to 'src')
| -rw-r--r-- | src/pj_factors.c | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/src/pj_factors.c b/src/pj_factors.c index 69d59e49..e48e8961 100644 --- a/src/pj_factors.c +++ b/src/pj_factors.c @@ -14,6 +14,7 @@ int pj_factors(LP lp, PJ *P, double h, struct FACTORS *fac) { double cosphi, t, n, r; int err; + PJ_COORD coo = {{lp.lam, lp.phi, 0, 0}}; err = proj_errno_reset (P); @@ -38,7 +39,7 @@ int pj_factors(LP lp, PJ *P, double h, struct FACTORS *fac) { /* If input latitudes are geocentric, convert to geographic */ if (P->geoc) - lp.phi = atan(P->rone_es * tan(lp.phi)); + lp = proj_geoc_lat (P, PJ_INV, coo).lp; /* If latitude + one step overshoots the pole, move it slightly inside, */ /* so the numerical derivative still exists */ |