eo-tides 0.6.0__tar.gz → 0.6.1__tar.gz

{eo_tides-0.6.0 → eo_tides-0.6.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: eo-tides
-Version: 0.6.0
+Version: 0.6.1
 Summary: Tide modelling tools for large-scale satellite earth observation analysis
 Author: Robbi Bishop-Taylor, Stephen Sagar, Claire Phillips, Vanessa Newey
 Author-email: Robbi.BishopTaylor@ga.gov.au
@@ -32,10 +32,11 @@ Requires-Dist: pandas>=2.2.0
 Requires-Dist: psutil>=5.8.0
 Requires-Dist: pyogrio>=0.10.0
 Requires-Dist: pyproj>=3.7.0
-Requires-Dist: pyTMD==2.2.1
+Requires-Dist: pyTMD>=2.2.2
 Requires-Dist: scikit-learn>=1.4.0
 Requires-Dist: scipy>=1.14.1
 Requires-Dist: shapely>=2.0.6
+Requires-Dist: timescale>=0.0.3
 Requires-Dist: tqdm>=4.55.0
 Requires-Dist: xarray>=2022.3.0
 Provides-Extra: notebooks

{eo_tides-0.6.0 → eo_tides-0.6.1}/eo_tides/model.py

@@ -19,6 +19,7 @@ import numpy as np
 import pandas as pd
 import pyproj
 import pyTMD
+import timescale.time
 from tqdm import tqdm
 
 from .utils import DatetimeLike, _set_directory, _standardise_models, _standardise_time, idw
@@ -95,7 +96,7 @@ def _model_tides(
     lon, lat = transformer.transform(x.flatten(), y.flatten())
 
     # Convert datetime
-    timescale = pyTMD.time.timescale().from_datetime(time.flatten())
+    ts = timescale.time.Timescale().from_datetime(time.flatten())
 
     try:
         # Read tidal constants and interpolate to grid points
@@ -136,10 +137,10 @@ def _model_tides(
     # Compute delta times based on model
     if pytmd_model.corrections in ("OTIS", "ATLAS", "TMD3", "netcdf"):
         # Use delta time at 2000.0 to match TMD outputs
-        deltat = np.zeros_like(timescale.tt_ut1)
+        deltat = np.zeros_like(ts.tt_ut1)
     else:
         # Use interpolated delta times
-        deltat = timescale.tt_ut1
+        deltat = ts.tt_ut1
 
     # In "one-to-many" mode, extracted tidal constituents and timesteps
     # are repeated/multiplied out to match the number of input points and
@@ -149,7 +150,7 @@ def _model_tides(
     points_repeat = len(x) if mode == "one-to-many" else 1
     time_repeat = len(time) if mode == "one-to-many" else 1
     t, hc, deltat = (
-        np.tile(timescale.tide, points_repeat),
+        np.tile(ts.tide, points_repeat),
         hc.repeat(time_repeat, axis=0),
         np.tile(deltat, points_repeat),
     )

{eo_tides-0.6.0 → eo_tides-0.6.1}/eo_tides/utils.py

@@ -181,9 +181,9 @@ def _clip_model_file(
     """
     Clips tide model netCDF datasets to a bounding box.
 
-    If the bounding box crosses 0 degrees longitude (e.g. Greenwich),
-    the function will clip the dataset into two parts and concatenate
-    them along the x-dimension to create a continuous result.
+    If the bounding box crosses 0 degrees longitude (e.g. Greenwich prime
+    meridian), the dataset will be clipped into two parts and concatenated
+    along the x-dimension to create a continuous result.
 
     Parameters
     ----------
@@ -226,20 +226,22 @@ def _clip_model_file(
     xcoords = nc[xcoord].compute()
     ycoords = nc[ycoord].compute()
 
-    # If data falls within 0-360 degree bounds, then clip directly
-    if (bbox.left >= 0) & (bbox.right <= 360):
+    # Convert longitudes to 0-360 convention
+    left = bbox.left % 360
+    right = bbox.right % 360
+
+    # If left coordinate is smaller than right, bbox does not cross
+    # zero longitude and can be clipped directly
+    if left <= right:  # bbox does not cross 0
         nc_clipped = nc.sel({
             ydim: (ycoords >= bbox.bottom) & (ycoords <= bbox.top),
-            xdim: (xcoords >= bbox.left) & (xcoords <= bbox.right),
+            xdim: (xcoords >= left) & (xcoords <= right),
         })
 
-    # If bbox crosses zero longitude, extract left and right
-    # separately and then combine into one concatenated dataset
-    elif (bbox.left < 0) & (bbox.right > 0):
-        # Convert longitudes to 0-360 range
-        left = bbox.left % 360
-        right = bbox.right % 360
-
+    # If left coordinate is larger than right, bbox crosses zero longitude.
+    # If so, extract left and right separately and then combine into one
+    # concatenated dataset
+    elif left > right:  # bbox crosses 0
         # Extract data from left of 0 longitude, and convert lon
         # coords to -180 to 0 range to enable continuous interpolation
         # across 0 boundary
@@ -357,7 +359,7 @@ def clip_models(
         model_files = model_database[m].get("model_file", [])
         grid_file = model_database[m].get("grid_file", [])
 
-        # Convert to list if strings and combine
+        # Convert to list of strings and combine
         model_files = model_files if isinstance(model_files, list) else [model_files]
         grid_file = grid_file if isinstance(grid_file, list) else [grid_file]
         all_files = model_files + grid_file

{eo_tides-0.6.0 → eo_tides-0.6.1}/eo_tides.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: eo-tides
-Version: 0.6.0
+Version: 0.6.1
 Summary: Tide modelling tools for large-scale satellite earth observation analysis
 Author: Robbi Bishop-Taylor, Stephen Sagar, Claire Phillips, Vanessa Newey
 Author-email: Robbi.BishopTaylor@ga.gov.au
@@ -32,10 +32,11 @@ Requires-Dist: pandas>=2.2.0
 Requires-Dist: psutil>=5.8.0
 Requires-Dist: pyogrio>=0.10.0
 Requires-Dist: pyproj>=3.7.0
-Requires-Dist: pyTMD==2.2.1
+Requires-Dist: pyTMD>=2.2.2
 Requires-Dist: scikit-learn>=1.4.0
 Requires-Dist: scipy>=1.14.1
 Requires-Dist: shapely>=2.0.6
+Requires-Dist: timescale>=0.0.3
 Requires-Dist: tqdm>=4.55.0
 Requires-Dist: xarray>=2022.3.0
 Provides-Extra: notebooks

{eo_tides-0.6.0 → eo_tides-0.6.1}/eo_tides.egg-info/requires.txt

@@ -7,10 +7,11 @@ pandas>=2.2.0
 psutil>=5.8.0
 pyogrio>=0.10.0
 pyproj>=3.7.0
-pyTMD==2.2.1
+pyTMD>=2.2.2
 scikit-learn>=1.4.0
 scipy>=1.14.1
 shapely>=2.0.6
+timescale>=0.0.3
 tqdm>=4.55.0
 xarray>=2022.3.0
 

{eo_tides-0.6.0 → eo_tides-0.6.1}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "eo-tides"
-version = "0.6.0"
+version = "0.6.1"
 description = "Tide modelling tools for large-scale satellite earth observation analysis"
 authors = [
     { name = "Robbi Bishop-Taylor" },
@@ -44,10 +44,11 @@ dependencies = [
     "psutil>=5.8.0",
     "pyogrio>=0.10.0",
     "pyproj>=3.7.0",
-    "pyTMD==2.2.1",
+    "pyTMD>=2.2.2",
     "scikit-learn>=1.4.0",
     "scipy>=1.14.1",
     "shapely>=2.0.6",
+    "timescale>=0.0.3",
     "tqdm>=4.55.0",
     "xarray>=2022.3.0",
 ]

{eo_tides-0.6.0 → eo_tides-0.6.1}/tests/test_utils.py

@@ -3,9 +3,11 @@ import tempfile
 from datetime import datetime
 
 import numpy as np
+import odc.geo.geom
 import pandas as pd
 import pytest
 
+from eo_tides.model import model_tides
 from eo_tides.utils import _standardise_models, _standardise_time, clip_models, idw, list_models
 
 
@@ -53,17 +55,96 @@ def test_standardise_models(model, ensemble_models, exp_process, exp_request, ex
 
 
 def test_clip_models():
+    # Set input and output paths
+    in_dir = "tests/data/tide_models"
     out_dir = pathlib.Path("tests/data/tide_models_clipped")
 
+    # Clip models to bbox
     clip_models(
-        input_directory="tests/data/tide_models",
+        input_directory=in_dir,
         output_directory=out_dir,
         bbox=(122.27, -18.07, 122.29, -18.05),
     )
 
+    # Assert that files were exported for all available models
     output_files = set([i.stem for i in out_dir.iterdir()])
     assert output_files == set(["GOT5", "EOT20", "hamtide"])
 
+    # Set modelling location
+    x, y = 122.28, -18.06
+    time = pd.date_range(start="2000-01", end="2001-03", freq="5h")
+
+    # Model using unclipped vs clipped files
+    df_unclipped = model_tides(
+        x=x,
+        y=y,
+        time=time,
+        model="HAMTIDE11",
+        directory=in_dir,
+        crop=False,
+    )
+    df_clipped = model_tides(
+        x=x,
+        y=y,
+        time=time,
+        model="HAMTIDE11",
+        directory=out_dir,
+        crop=False,
+    )
+
+    # Verify both produce the same results
+    assert np.allclose(df_unclipped.tide_height, df_clipped.tide_height)
+
+
+# Test clipping across multiple global locations using synthetic HAMTIDE11 data
+@pytest.mark.parametrize(
+    "bbox, name",
+    [
+        ((-166, 14, -151, 29), "hawaii"),  # entirely W of prime meridian
+        ((-13, 49, 6, 60), "uk"),  # crossing prime meridian
+        ((105.292969, -47.872144, 160.312500, -5.266008), "aus"),  # entirely E of prime meridian
+        ((-256.640625, 7.013668, -119.794922, 63.391522), "pacific"),  # crossing antimeridian
+    ],
+)
+def test_clip_models_bbox(bbox, name):
+    # Set input and output paths
+    in_dir = "tests/data/tide_models_synthetic/"
+    out_dir = f"tests/data/tide_models_synthetic_{name}/"
+
+    # Clip models to input bbox
+    clip_models(
+        input_directory=in_dir,
+        output_directory=out_dir,
+        bbox=bbox,
+        model="HAMTIDE11",
+        overwrite=True,
+    )
+
+    # Set modelling location based on bbox centroid
+    x, y = odc.geo.geom.BoundingBox(*bbox, crs="EPSG:4326").polygon.centroid.xy
+    time = pd.date_range(start="2000-01", end="2001-03", freq="5h")
+
+    # Model using unclipped vs clipped files
+    df_unclipped = model_tides(
+        x=x,
+        y=y,
+        time=time,
+        model="HAMTIDE11",
+        directory=in_dir,
+        crop=False,
+    )
+    df_clipped = model_tides(
+        x=x,
+        y=y,
+        time=time,
+        model="HAMTIDE11",
+        directory=out_dir,
+        crop=False,
+    )
+
+    # Verify both produce the same results
+    assert np.allclose(df_unclipped.tide_height, df_clipped.tide_height)
+
 
 @pytest.mark.parametrize(
     "input_value, expected_output",