clip at -1000, not 0

scottstanie · Dec 16, 2018 · 0d0e55b · 0d0e55b
1 parent e5c250d
commit 0d0e55b
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 3 deletions.
diff --git a/sardem/loading.py b/sardem/loading.py
@@ -59,8 +59,10 @@ def load_elevation(filename):
             dem_img = data.reshape((1201, 1201))
         else:
             raise ValueError("Invalid .hgt data size: must be square size 1201 or 3601")
-        # TODO: do we always want this? Not sure
-        dem_img = np.clip(dem_img, 0, None)
+        # TODO: Verify that the min real value will be above -1000
+        min_valid = -1000
+        # Set NaN values to 0
+        dem_img[dem_img < min_valid] = 0
 
     return dem_img
 

diff --git a/setup.py b/setup.py
@@ -5,7 +5,7 @@
 
 setuptools.setup(
     name="sardem",
-    version="0.2.4",
+    version="0.2.5",
     author="Scott Staniewicz",
     author_email="scott.stanie@utexas.com",
     description="Create upsampled DEMs for InSAR processing",