Changeset 26701

Timestamp:

12/02/21 14:05:48 (3 years ago)

Author:

Cheng Gong

Message:

ADD:interpret MEaSUREs Greenland Ice velocity to mesh

Location:

issm/trunk-jpl/src/m/contrib

Files:

• chenggong (added)
• chenggong/interpFromMEaSUREsGeotiff.m (added)
• morlighem/modeldata/interpFromGeotiff.m (modified) (3 diffs)

issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpFromGeotiff.m

@@ -1 @@
-function dataout = interpFromGeotiff(geotiffname,X,Y),
+function dataout = interpFromGeotiff(geotiffname,X,Y,nanValue),
+
+if nargin < 4
+        nanValue = 10^30;
+end
 
 usemap = 0;
@@ -21 +25 @@
         %Get image info
         Tinfo = imfinfo(geotiffname);
-        N     = Tinfo.Width;
-        M     = Tinfo.Height;
-        dx    = Tinfo.ModelPixelScaleTag(1);
-        dy    = Tinfo.ModelPixelScaleTag(2);
-        minx  = Tinfo.ModelTiepointTag(4);
-        maxy  = Tinfo.ModelTiepointTag(5);
+        N     = Tinfo(1).Width;
+        M     = Tinfo(1).Height;
+        dx    = Tinfo(1).ModelPixelScaleTag(1);
+        dy    = Tinfo(1).ModelPixelScaleTag(2);
+        minx  = Tinfo(1).ModelTiepointTag(4);
+        maxy  = Tinfo(1).ModelTiepointTag(5);
 
         %Generate vectors
@@ -55 +59 @@
                 data=double(flipud(imread(geotiffname)));
         end
-        data(find(abs(data)>10^30))=NaN;
+        if nanValue > 0
+                data(find(abs(data)>=nanValue))=NaN;
+        else 
+                data(find(data<=nanValue))=NaN;
+        end
 
 end