Changeset 26958

Timestamp:

04/15/22 12:30:25 (3 years ago)

Author:

Mathieu Morlighem

Message:

CHG: switching to CoordTransform to avoid EPSG problems

Location:

issm/trunk-jpl/src/m/coordsystems

Files:

• ll2xy.m (modified) (1 diff)
• xy2ll.m (modified) (1 diff)

issm/trunk-jpl/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
+                assert(delta==45); assert(slat ==70);
+                [x y]=CoordTransform(lat,lon,'EPSG:4326','EPSG:3413');
+                return;
+        else
+                assert(delta==0); assert(slat ==71);
+                [x y]=CoordTransform(lat,lon,'EPSG:4326','EPSG:3031');
+                return;
+        end
+end
+
 % Conversion constant from degrees to radians
 cde  = 57.29577951;

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

@@ -34 +34 @@
         help xy2ll
         error('bad usage');
+end
+
+%Move to CoordTransform now...
+if exist('CoordTransform_matlab')==3
+        disp('Calling CoordTransform instead, make sure to change your MATLAB script');
+        if sgn==+1
+                assert(delta==45); assert(slat ==70);
+                [lat lon]=CoordTransform(x, y,'EPSG:3413','EPSG:4326');
+                return;
+        else
+                assert(delta==0); assert(slat ==71);
+                [lat lon]=CoordTransform(x, y,'EPSG:3031','EPSG:4326');
+                return;
+        end
 end