Merge pull request #13 from cirrusasf/jz_pfa_split

v0.7.3 fix a calculation of min and max height of DEM to speed up the RTC process

Author: cirrusasf

CHANGELOG.md

@@ -6,6 +6,10 @@ The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/),
 and this project adheres to [PEP 440](https://www.python.org/dev/peps/pep-0440/)
 and uses [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
 
+## [0.7.3]
+#### Fixed
+* correly calcualte the min and max values of the dem raster to produce geogrid, largely imporve the RTC process speed.
+
 ## [0.7.2]
 ### Fixed
 * clip the RTC output of UMBRA to the shape of what is defined in the json file. 

src/multirtc/define_geogrid.py

@@ -136,15 +136,9 @@ def generate_geogrids(
         min_height = isce3.core.MINIMUM_HEIGHT
         max_height = isce3.core.MAXIMUM_HEIGHT
     else:
-        dem_raster = isce3.io.Raster(str(dem_path))
-        dem = isce3.geometry.DEMInterpolator()
-        # FIXME: figure out why load with bounds doesn't work
-        # minx, miny, maxx, maxy = slc.footprint.bounds
-        # dem.load_dem(dem_raster, min_x=minx, max_x=maxx, min_y=miny, max_y=maxy)
-        dem.load_dem(dem_raster)
-        dem.compute_min_max_mean_height()
-        min_height = dem.min_height
-        max_height = dem.max_height
+        with rasterio.open(str(dem_path)) as src:
+            min_height = (src.stats()[0]).min
+            max_height = (src.stats()[0]).max
 
     geogrid = isce3.product.bbox_to_geogrid(
         slc.radar_grid,