Update fillsinks.m

Removed bug when calling fillsinks with an additional sink raster.
csdms-contrib · Nov 8, 2015 · e1c2f82 · e1c2f82
1 parent 41ec27f
commit e1c2f82
Showing 1 changed file with 4 additions and 3 deletions.
diff --git a/@GRIDobj/fillsinks.m b/@GRIDobj/fillsinks.m
@@ -133,8 +133,6 @@
  I = (imdilate(Inan,ones(3)) & ~Inan) | I;
  end
 
- SINKS = SINKS | I;
-
  % Refine markers by identifying the lowest value in each sink
  STATS = regionprops(SINKS,dem,'MinIntensity','PixelIdxList');
  SINKS = false(size(SINKS));
@@ -145,10 +143,13 @@
 
  marker = -inf(size(dem),class(dem));
  marker(SINKS) = -dem(SINKS);
+ marker(I) = -dem(I);
  demfs = -imreconstruct(marker,-dem);
 
 end
 
 % nans in the dem are set to nan again
-demfs(Inan) = nan;
+try
+ demfs(Inan) = nan;
+end
 DEM.Z = demfs;