Changeset 22758 for issm/trunk/src/m/coordsystems/xy2ll.m

Timestamp:

05/10/18 10:24:27 (7 years ago)

Author:

Mathieu Morlighem

Message:

merged trunk-jpl and trunk for revision 22757

Location:

issm/trunk

Files:

• . (modified) (1 prop)
• src (modified) (1 prop)
• src/m/coordsystems/xy2ll.m (modified) (5 diffs)

issm/trunk/src/m/coordsystems/xy2ll.m

@@ -1 @@
-function [lat,lon] = xy2ll(x,y,sgn,central_meridian,standard_parallel)
+function [lat,lon,scale_factor] = xy2ll(x,y,sgn,central_meridian,standard_parallel)
 %XY2LL - converts xy to lat long
 %
@@ -5 +5 @@
 %   latitude and longitude Stereographic (X,Y) coordinates for the polar
 %   regions.
+%   Optional scale factor provides the scaling factor needed to correct projection error
+%   in areas and volumes
 %   Author: Michael P. Schodlok, December 2003 (map2xy.m)
 %
 %   Usage:
 %      [lat,lon] = xy2ll(x,y,sgn);
+%      [lat,lon,scale_factor] = xy2ll(x,y,sgn);
 %      [lat,lon] = xy2ll(x,y,sgn,central_meridian,standard_parallel);
 %
@@ -29 +32 @@
         end
 else
+        help xy2ll
+        error('bad usage');
+end
+
+if nargout~=3 & nargout~=2,
         help xy2ll
         error('bad usage');
@@ -64 +72 @@
 [res1,res2] = find(rho(:) <= 0.1);
 if res1
-        lat(res1) = 90. * sgn;
+        lat(res1) = pi/2. * sgn;
         lon(res1) = 0.0;
 end
@@ -71 +79 @@
 lat = lat * 180. / pi;
 lon = lon - delta; 
+if nargout==3,
+        m=((1+sin(abs(slat)*pi/180))*ones(length(lat),1)./(1+sin(abs(lat)*pi/180)));
+        scale_factor=(1./m).^2;
+end