save_hdfeos5: UTM support for --subset mode

insarlab · Mar 18, 2024 · c64d497 · c64d497
1 parent 0e8fe23
commit c64d497
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/mintpy/save_hdfeos5.py b/src/mintpy/save_hdfeos5.py
@@ -243,6 +243,10 @@ def get_output_filename(metadata, suffix=None, update_mode=False, subset_mode=Fa
         lat0 = lat1 + float(metadata['Y_STEP']) * int(metadata['LENGTH'])
         lon1 = lon0 + float(metadata['X_STEP']) * int(metadata['WIDTH'])
 
+        # convert UTM to lat/lon
+        if not metadata.get('Y_UNIT', 'degrees').lower().startswith('deg'):
+            [lat0, lat1], [lon0, lon1] = ut.utm2latlon(metadata, easting=[lon0, lon1], northing=[lat0, lat1])
+
         lat0Str = f'N{round(lat0*1e3):05d}'
         lat1Str = f'N{round(lat1*1e3):05d}'
         lon0Str = f'E{round(lon0*1e3):06d}'