fixed missing reprojection of pixel coordinates in mapper module

EOA-team · Oct 18, 2022 · 2c896a5 · 2c896a5
1 parent 279b7bb
commit 2c896a5
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 0 deletions.
diff --git a/eodal/operational/mapping/mapper.py b/eodal/operational/mapping/mapper.py
@@ -556,6 +556,11 @@ def _get_observation(
                 res["sensing_date"] = scenes_date["sensing_date"].values[0]
                 res["scene_id"] = scenes_date["scene_id"].values[0]
                 res['sensing_time'] = scenes_date['sensing_time'].values[0]
+                # make sure the result is projected into the target EPSG code, otherwise
+                # the resulting GeoDataFrame contains wrong coordinates, i.e., the
+                # coordinates were not projected into the target CRS
+                if res.crs.to_epsg() != scenes_date.target_crs.unique()[0]:
+                    res.to_crs(epsg=scenes_date.target_crs.unique()[0], inplace=True)
             # or the feature
             else:
                 if sensor.lower() == 'sentinel1':

diff --git a/eodal/operational/mapping/sentinel2.py b/eodal/operational/mapping/sentinel2.py
@@ -205,6 +205,10 @@ def _read_multiple_scenes(
                 if feature_gdf.empty:
                     continue
                 res = feature_gdf
+                # make sure the coordinates are reprojected if necessary
+                if res.crs.to_epsg() != scenes_date.target_crs.unique()[0]:
+                    res.to_crs(epsg=scenes_date.target_crs.unique()[0], inplace=True)
+
                 if isinstance(candidate_scene, (pd.Series, gpd.GeoSeries)):
                     res["sensing_date"] = candidate_scene["sensing_date"]
                     res['sensing_time'] = candidate_scene['sensing_time']