Changeset 27912

Timestamp:

09/20/23 09:56:10 (18 months ago)

Author:

Mathieu Morlighem

Message:

CHG: added fill holes option

File:

• issm/trunk-jpl/src/m/modeldata/interpFromGeotiff.m (modified) (2 diffs)

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

@@ -1 @@
-function dataout = interpFromGeotiff(geotiffname,X,Y,nanValue),
+function dataout = interpFromGeotiff(geotiffname,X,Y,nanValue,fillholes)
+%INTERPFROMGEOTIFF - interpolate field in geotiff onto list of points
+%
+%   Usage:
+%      dataout = interpFromGeotiff(geotiffname,X,Y,nanValue,fillholes)
+%      dataout = interpFromGeotiff(geotiffname,X,Y);
+
 
 if nargin < 4
         nanValue = 10^30;
+        fillholes = false;
+end
+if nargin < 5
+        fillholes = false;
 end
 
@@ -64 +74 @@
                 data(find(data<=nanValue))=NaN;
         end
-
+        if fillholes
+                disp('Filling holes');
+                data = inpaint_nans(data);
+                disp('done');
+        end
 end