[Octave-patch-tracker] [patch #10022] [octave forge] (mapping) Auxillary Latitudes

[Philip Nienhuis]

Follow-up Comment #19, patch #10022 (project octave):

Thanks very much.
Now I get:

>> iso2geo (geo2iso ([-90 : 10 : 90]))
ans =
 Columns 1 through 17:
  -90  -80  -70  -60  -50  -40  -30  -20  -10    0   10   20   30   40   50  
60   70
 Columns 18 and 19:
   80   90


which looks good!
But:
I'm astounded because of this in geo2iso:

:
  elseif (strncmpi (angleUnit, "degrees", min (length (angleUnit), 7)))
:   phi = deg2rad (phi);

:
  chi = geo2con (phi);
:


where in geo2iso.m:
* I'd expect the call 2 geo2con to also contain input args E and "radians", as
geo2con also needs the ellipsoid and higher up, phi was morphed into radians;
* geo2con.m expects "degrees" by default.

If I change the call 2 geocon() into

  chi = geo2con (phi, E, "radians");


I get:

>> geo2iso ([-90 : 10 : 90], "wgs84", "degrees")
ans =
 Columns 1 through 6:
         -Inf  -1.3921e+02  -9.9071e+01  -7.5123e+01  -5.7614e+01 
-4.3465e+01
 Columns 7 through 12:
  -3.1281e+01  -2.0288e+01  -9.9846e+00  -6.3611e-15   9.9846e+00  
2.0288e+01
 Columns 13 through 18:
   3.1281e+01   4.3465e+01   5.7614e+01   7.5123e+01   9.9071e+01  
1.3921e+02
 Column 19:
          Inf


which is the exact same output as Matlab.

So the new iso2geo.m still yields erroneous values. But that is easily solved
as in iso2geo you likewise need to feed args 'E, "radians"' to the input args
of con2geo() ...

OK I'll polish style and Texinfo and push geo2iso.m & iso2geo.m later on. With
that done it's time for a new mapping pkg release.


    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/patch/?10022>

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/