[Octave-patch-tracker] [patch #10021] [octave forge] (mapping) meridianfwd

[Philip Nienhuis]

Follow-up Comment #5, patch #10021 (project octave):

Also,

>> [lat, lon] = geodeticfwd (0, 0, [0 90], 2 * pi() / 2, "Unit Sphere",
"degrees")
lat =
   7.0165e-15   4.2962e-31

lon =
            0   7.0165e-15


where I'd expect lon to be [180 180] (going from equator to antipode on
equator, over the North Pole and along the equator).
In fact any azimuth value should lead to this result.
E.g.,

>> [latr, lonr] = reckon (0, 0, 180, [-90 : 10: 90])
latr =
   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0

lonr =
 Columns 1 through 14:
  -180  -180  -180  -180  -180  -180  -180  -180  -180  -180  -180  -180  -180
 -180
 Columns 15 through 19:
  -180  -180  -180  -180  -180


Moreover,

>> [lat, lon] = geodeticfwd (0, 0, [-90 : 30 : 90], pi()/2, "Unit Sphere",
"degrees")
lat =
 Columns 1 through 6:
   3.5082e-15   3.0000e+01   6.0000e+01   9.0000e+01   6.0000e+01  
3.0000e+01
 Column 7:
   3.5082e-15

lon =
 Columns 1 through 4:
  -90.0000 +       0i  -90.0000 +       0i  -90.0000 +  0.0000i         0 +   
   0i
 Columns 5 through 7:
   90.0000 -  0.0000i   90.0000 +       0i   90.0000 +       0i

>> [latr, lonr] = reckon (0, 0, 90, [-90 : 30: 90])
latr =
         0   30.0000   60.0000   90.0000   60.0000   30.0000         0

lonr =
  -90  -90  -90    0   90   90   90


I'd first try to let geodeticfwd.m match reckon.m's result for a Unit Sphere
(note: radius =1.000).
Latitude seem to match, longitudes don't.

    _______________________________________________________

Reply to this item at:

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

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