Changeset 28013 for issm/trunk/src/m/coordsystems/ll2xy.m

Timestamp:

11/15/23 12:14:04 (16 months ago)

Author:

Mathieu Morlighem

Message:

merged trunk-jpl and trunk for revision 28011

Location:

issm/trunk

Files:

• . (modified) (1 prop)
• src (modified) (1 prop)
• src/m/coordsystems/ll2xy.m (modified) (1 diff)

issm/trunk/src/m/coordsystems/ll2xy.m

@@ -40 +40 @@
 end
 
-%Move to CoordTransform now...
-if exist('CoordTransform_matlab')==3
-        disp('Calling CoordTransform instead, make sure to change your MATLAB script');
-        if sgn==+1
-                if delta==45 && slat==70      %BedMachine
-                        [x y]=CoordTransform(lat,lon,'EPSG:4326','EPSG:3413'); return;
-                elseif delta==0 && slat==75   %IBCAO
-                        [x y]=CoordTransform(lat,lon,'EPSG:4326','EPSG:3996'); return;
-                elseif delta==39 && slat==71  %Bamber
-                        Bamber_proj = '+proj=stere +lat_0=90 +lat_ts=71 +lon_0=-39 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs';
-                        [x y]=CoordTransform(lat,lon,'EPSG:4326',Bamber_proj); return;
-                else
-                        error('not supported yet');
-                end
-        else
-                assert(delta==0); assert(slat ==71);
-                [x y]=CoordTransform(lat,lon,'EPSG:4326','EPSG:3031');
-                return;
-        end
+%Choose ellipsoid
+if 0
+        %Hughes ellipsoid
+        re   = 6378.273*10^3; % Radius of the earth in meters
+        ex2 = .006693883;     % Eccentricity of the Hughes ellipsoid squared
+else
+        %WGS84 ellipsoid
+        re = 6378137;         % Radius of the earth in meters
+        f  = 1./298.257223563;% Earth flattening
+        ex2 = 2*f-f^2;        % Eccentricity squared
 end
 
-% Conversion constant from degrees to radians
-cde  = 57.29577951;
-% Radius of the earth in meters
-re   = 6378.273*10^3;
-% Eccentricity of the Hughes ellipsoid squared
-ex2   = .006693883;
-% Eccentricity of the Hughes ellipsoid
-ex    =  sqrt(ex2);
+% Eccentricity 
+ex = sqrt(ex2);
 
 latitude  = abs(lat) * pi/180.;