dist-s1-enumerator 1.0.5__tar.gz → 1.0.7__tar.gz

{dist_s1_enumerator-1.0.5 → dist_s1_enumerator-1.0.7}/CHANGELOG.md

@@ -6,6 +6,19 @@ 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).
 
+## [1.0.7] - 2025-10-24
+
+### Added
+* Updated print statements for multi-window to show the max number of pre-images per year.
+
+
+## [1.0.6] - 2025-10-09
+
+### Added
+* Tests for MGRS tiles and burst geometries ensuring they overlap the antimeridian if and only if they are multipolygons.
+* Updated readme with information about these mgrs/burst geometries and the dateline.
+* Updated burst geometries using geopackage geometry from opera_adt/burst_db. See this issue: https://github.com/opera-adt/burst_db/issues/120 - some dateline geometries were cut off. The new geometries are larger and the table went from 24 MB to 37 MB. See this notebook: https://github.com/OPERA-Cal-Val/dist-s1-research/blob/dev/marshak/Zb_update_burst_table/Update%20Burst%20Geometries%20Using%20Geopackage.ipynb
+
 ## [1.0.5] - 2025-09-29
 
 ### Fixed

{dist_s1_enumerator-1.0.5/src/dist_s1_enumerator.egg-info → dist_s1_enumerator-1.0.7}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: dist-s1-enumerator
-Version: 1.0.5
+Version: 1.0.7
 Summary: Enumeration and ops library for the OPERA DIST-S1 project
 Author-email: "Richard West, Charlie Marshak, Talib Oliver-Cabrera, and Jungkyo Jung" <charlie.z.marshak@jpl.nasa.gov>
 License: Apache-2.0
@@ -125,6 +125,15 @@ There are two category of tests: unit tests and integration tests. The former ca
 The integration tests that are the most time consuming are represented by the notebooks and are run only upon a release PR.
 These notebook tests are tagged with `notebooks` and can be excluded from the other tests with `pytest tests -m 'not notebooks'`.
 
+# Remarks about the Dateline/Dateline and Geometry
+
+The antimeridian (or dateline) is the line at the -180 longitude mark that global CRS tiles are wrapped by standard global reference systems.
+The geometries of the bursts and the MGRS tiles in this package are all in `epsg:4326` (standard lon/lat). 
+The geometries are all between -180 and 180 so those geometries that cross the antimeridian/dateline are generally wrapped.
+For MGRS tiles, the statement that a geometry overlaps the antimeridian occurs if and only if the geometry is a Polygon.
+The same is true for burst geometries.
+See `test_antimeridian_crossing` in [`tests/test_mgrs_burst_data.py`](tests/test_mgrs_burst_data.py).
+
 # Contributing
 
 We welcome contributions to this open-source package. To do so:

{dist_s1_enumerator-1.0.5 → dist_s1_enumerator-1.0.7}/README.md

@@ -102,6 +102,15 @@ There are two category of tests: unit tests and integration tests. The former ca
 The integration tests that are the most time consuming are represented by the notebooks and are run only upon a release PR.
 These notebook tests are tagged with `notebooks` and can be excluded from the other tests with `pytest tests -m 'not notebooks'`.
 
+# Remarks about the Dateline/Dateline and Geometry
+
+The antimeridian (or dateline) is the line at the -180 longitude mark that global CRS tiles are wrapped by standard global reference systems.
+The geometries of the bursts and the MGRS tiles in this package are all in `epsg:4326` (standard lon/lat). 
+The geometries are all between -180 and 180 so those geometries that cross the antimeridian/dateline are generally wrapped.
+For MGRS tiles, the statement that a geometry overlaps the antimeridian occurs if and only if the geometry is a Polygon.
+The same is true for burst geometries.
+See `test_antimeridian_crossing` in [`tests/test_mgrs_burst_data.py`](tests/test_mgrs_burst_data.py).
+
 # Contributing
 
 We welcome contributions to this open-source package. To do so:

{dist_s1_enumerator-1.0.5 → dist_s1_enumerator-1.0.7}/src/dist_s1_enumerator/dist_enum.py

@@ -123,7 +123,10 @@ def enumerate_one_dist_s1_product(
     if lookback_strategy == 'immediate_lookback':
         # Add 5 minutes buffer to ensure we don't include post-images in pre-image set.
         print('Searching for pre-images for immediate_lookback products')
-        print(f'Lookback days {params.delta_lookback_days} and window days {params.delta_window_days}')
+        print(
+            f'Lookback days {params.delta_lookback_days} and window days {params.delta_window_days} '
+            f'with max pre-images per burst {params.max_pre_imgs_per_burst}'
+        )
         post_date_min = df_rtc_post.acq_dt.min() - pd.Timedelta(seconds=300)
         earliest_lookback = params.delta_window_days + params.delta_lookback_days
         latest_lookback = params.delta_lookback_days
@@ -146,7 +149,10 @@ def enumerate_one_dist_s1_product(
         df_rtc_pre_list = []
         zipped_data = list(zip(params.delta_lookback_days, params.max_pre_imgs_per_burst))
         print('Searching for pre-images for multi_window baseline')
-        print(f'Lookback days {params.delta_lookback_days} and window days {params.delta_window_days}')
+        print(
+            f'Lookback days {params.delta_lookback_days} and window days {params.delta_window_days} '
+            f'with max pre-images per burst {params.max_pre_imgs_per_burst}'
+        )
         for delta_lookback_day, max_pre_img_per_burst in tqdm(
             zipped_data,
             desc='Windows',

{dist_s1_enumerator-1.0.5 → dist_s1_enumerator-1.0.7/src/dist_s1_enumerator.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: dist-s1-enumerator
-Version: 1.0.5
+Version: 1.0.7
 Summary: Enumeration and ops library for the OPERA DIST-S1 project
 Author-email: "Richard West, Charlie Marshak, Talib Oliver-Cabrera, and Jungkyo Jung" <charlie.z.marshak@jpl.nasa.gov>
 License: Apache-2.0
@@ -125,6 +125,15 @@ There are two category of tests: unit tests and integration tests. The former ca
 The integration tests that are the most time consuming are represented by the notebooks and are run only upon a release PR.
 These notebook tests are tagged with `notebooks` and can be excluded from the other tests with `pytest tests -m 'not notebooks'`.
 
+# Remarks about the Dateline/Dateline and Geometry
+
+The antimeridian (or dateline) is the line at the -180 longitude mark that global CRS tiles are wrapped by standard global reference systems.
+The geometries of the bursts and the MGRS tiles in this package are all in `epsg:4326` (standard lon/lat). 
+The geometries are all between -180 and 180 so those geometries that cross the antimeridian/dateline are generally wrapped.
+For MGRS tiles, the statement that a geometry overlaps the antimeridian occurs if and only if the geometry is a Polygon.
+The same is true for burst geometries.
+See `test_antimeridian_crossing` in [`tests/test_mgrs_burst_data.py`](tests/test_mgrs_burst_data.py).
+
 # Contributing
 
 We welcome contributions to this open-source package. To do so:

{dist_s1_enumerator-1.0.5 → dist_s1_enumerator-1.0.7}/tests/test_mgrs_burst_data.py

@@ -1,5 +1,5 @@
 import pytest
-from shapely.geometry import Point
+from shapely.geometry import LineString, MultiPolygon, Point, Polygon
 
 from dist_s1_enumerator.exceptions import NoMGRSCoverage
 from dist_s1_enumerator.mgrs_burst_data import (
@@ -161,3 +161,20 @@ def test_all_bursts_in_lut() -> None:
     df_merged = df_bursts.merge(df_mgrs_lut, on='jpl_burst_id', indicator=True, how='left')
     burst_ids_not_in_lut = df_merged[df_merged['_merge'] == 'left_only'].jpl_burst_id.unique().tolist()
     assert len(burst_ids_not_in_lut) == 0
+
+
+def test_antimeridian_crossing() -> None:
+    df_mgrs = get_mgrs_table()
+    df_burst = get_burst_table()
+    antimeridian_0 = LineString(coordinates=((-180, 90), (-180, -90))).buffer(0.00000001)
+    antimeridian_1 = LineString(coordinates=((180, 90), (180, -90))).buffer(0.00000001)
+
+    for df in [df_mgrs, df_burst]:
+        for antimeridian in [antimeridian_0, antimeridian_1]:
+            ind_anti = df.geometry.intersects(antimeridian)
+            df_antimerid = df[ind_anti].reset_index(drop=True)
+            any_polygons = (df_antimerid.geometry.map(lambda geo: isinstance(geo, Polygon))).sum()
+            assert any_polygons == 0
+            df_antimerid_not = df[~ind_anti].reset_index(drop=True)
+            any_multis = (df_antimerid_not.geometry.map(lambda geo: isinstance(geo, MultiPolygon))).sum()
+            assert any_multis == 0