Change nodata in rtc dem

CHANGELOG.md

@@ -6,6 +6,12 @@ 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).
 
+## [6.2.4]
+
+### Changed
+- Use the envelope of the map-overlay.kml to determine the input dem.tif
+- Remove the NoData setting in the output dem.tif
+
 ## [6.2.3]
 
 ### Added

hyp3_gamma/dem.py

@@ -80,6 +80,26 @@ def shift_for_antimeridian(dem_file_paths: List[str], directory: Path) -> List[s
     return shifted_file_paths
 
 
+def get_envelope_geometry(geometry):
+    # get the envelope of the geometry
+    geometry_out = ogr.Geometry(ogr.wkbMultiPolygon)
+    geoms = [g for g in geometry]
+    boxes = []
+    for g in geoms:
+        boxes.append(g.GetEnvelope())
+
+    minlat = min([x[2] for x in boxes])
+    maxlat = max([x[3] for x in boxes])
+
+    for item in boxes:
+        wkt = f'POLYGON (({item[0]}  {minlat}, {item[0]} {maxlat}, {item[1]} {maxlat}, {item[1]} {minlat},' \
+               f' {item[0]} {minlat}))'
+        poly = ogr.CreateGeometryFromWkt(wkt)
+        geometry_out.AddGeometry(poly)
+
+    return geometry_out
+
+
 def prepare_dem_geotiff(output_name: str, geometry: ogr.Geometry, pixel_size: float = 30.0):
     """Create a DEM mosaic GeoTIFF covering a given geometry.
 
@@ -99,11 +119,14 @@ def prepare_dem_geotiff(output_name: str, geometry: ogr.Geometry, pixel_size: fl
         with TemporaryDirectory() as temp_dir:
             temp_path = Path(temp_dir)
 
+            geometry = get_envelope_geometry(geometry)
+
             centroid = geometry.Centroid()
             dem_file_paths = get_dem_file_paths(geometry.Buffer(0.15))
 
             if geometry.GetGeometryName() == 'MULTIPOLYGON':
                 centroid = get_centroid_crossing_antimeridian(geometry)
+
                 dem_file_paths = shift_for_antimeridian(dem_file_paths, temp_path)
 
             dem_vrt = temp_path / 'dem.vrt'

hyp3_gamma/rtc/rtc_sentinel.py

@@ -165,6 +165,7 @@ def prepare_dem(safe_dir: str, dem_name: str, bbox: List[float] = None, dem: str
             geometry = ogr.CreateGeometryFromWkt(wkt)
         else:
             geometry = get_geometry_from_kml(f'{safe_dir}/preview/map-overlay.kml')
+
         prepare_dem_geotiff(dem_tif, geometry, pixel_size)
         run(f'dem_import {dem_tif} {dem_image} {dem_par} - - $DIFF_HOME/scripts/egm2008-5.dem '
             f'$DIFF_HOME/scripts/egm2008-5.dem_par - - - 1')
@@ -390,7 +391,8 @@ def rtc_sentinel_gamma(safe_dir: str, resolution: float = 30.0, radiometry: str
     if include_dem:
         with NamedTemporaryFile() as temp_file:
             run(f'data2geotiff dem_seg.par dem_seg 2 {temp_file.name}')
-            gdal.Translate(f'{product_name}/{product_name}_dem.tif', temp_file.name, outputType=gdalconst.GDT_Int16)
+            gdal.Translate(f'{product_name}/{product_name}_dem.tif', temp_file.name, noData="none",
+                           outputType=gdalconst.GDT_Int16)
     if include_inc_map:
         run(f'data2geotiff dem_seg.par corrected.inc_map 2 {product_name}/{product_name}_inc_map.tif')
     if include_scattering_area: