transformez 0.4.2__tar.gz → 0.5.0__tar.gz

{transformez-0.4.2 → transformez-0.5.0}/CHANGELOG.md

@@ -5,6 +5,21 @@ All notable changes to this project will be documented in this file.
 The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/),
 and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
 
+## [0.4.4] - 2026-06-25
+
+### Added
+- support for projected input rasters in the raster command
+- add 'save_shift' to transform_raster api/cli
+
+### CHANGED
+- rety failed downloads, such as FES
+
+## [0.4.3] - 2026-05-04
+
+### Added
+- new logo
+- force htdp 3.5.0
+
 ## [0.3.5] - 2026-04-08
 
 ### Added

{transformez-0.4.2 → transformez-0.5.0}/CITATION.cff

@@ -14,6 +14,6 @@ authors:
 
 website: https://transformez.readthedocs.io
 title: "Transformez"
-version: 0.3.6
-date-released: 2026-04-08
+version: 0.5.0
+date-released: 2026-06-06
 url: "https://github.com/continuous-dems/transformez"

{transformez-0.4.2 → transformez-0.5.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: transformez
-Version: 0.4.2
+Version: 0.5.0
 Summary: A standalone utility for vertical elevation datum transformations.
 Project-URL: Homepage, https://github.com/continuous-dems/transformez
 Project-URL: Issues, https://github.com/continuous-dems/transformez/issues
@@ -47,16 +47,16 @@ Provides-Extra: preview
 Requires-Dist: matplotlib; extra == 'preview'
 Description-Content-Type: text/markdown
 
-<!-- <p align="center"> -->
-<!-- 	<a href="https://github.com/continuous-dems"> -->
-<!-- 		<img src="https://github.com/continuous-dems/fetchez/blob/modules/docs/source/_static/continuous_dems_logo.svg" height="80" alt="Continuous DEMs Logo"> -->
-<!-- 	</a> -->
-<!-- </p> -->
+<p align="center">
+	<a href="https://github.com/continuous-dems">
+		<img src="https://raw.githubusercontent.com/continuous-dems/transformez/refs/heads/main/docs/source/_static/continuous_dems_logo_mini_inf.svg" height="80" alt="Continuous DEMs Logo">
+	</a>
+</p>
 <h1 align="center">Transformez</h1>
 <p align="center"><strong>Global vertical datum transformations, simplified.</strong></p>
 
 <p align="center">
-  <a href="https://github.com/continuous-dems/transformez"><img src="https://img.shields.io/badge/version-0.3.6-blue.svg" alt="Version"></a>
+  <a href="https://github.com/continuous-dems/transformez"><img src="https://img.shields.io/badge/version-0.5.0-blue.svg" alt="Version"></a>
   <a href="LICENSE"><img src="https://img.shields.io/badge/license-MIT-green.svg" alt="License"></a>
   <a href="https://www.python.org/"><img src="https://img.shields.io/badge/python-3.12+-yellow.svg" alt="Python"></a>
   <a href="https://badge.fury.io/py/transformez"><img src="https://badge.fury.io/py/transformez.svg" alt="PyPI version"></a>
@@ -127,4 +127,5 @@ Would you like to know more? Check out our [Official Documentation](https://tran
 ## ⚖ License
 
 This project is licensed under the MIT License - see the [LICENSE](https://github.com/ciresdem/transformez/blob/main/LICENSE) file for details.
+
 Copyright (c) 2010-2026 Regents of the University of Colorado

{transformez-0.4.2 → transformez-0.5.0}/README.md

@@ -1,13 +1,13 @@
-<!-- <p align="center"> -->
-<!-- 	<a href="https://github.com/continuous-dems"> -->
-<!-- 		<img src="https://github.com/continuous-dems/fetchez/blob/modules/docs/source/_static/continuous_dems_logo.svg" height="80" alt="Continuous DEMs Logo"> -->
-<!-- 	</a> -->
-<!-- </p> -->
+<p align="center">
+	<a href="https://github.com/continuous-dems">
+		<img src="https://raw.githubusercontent.com/continuous-dems/transformez/refs/heads/main/docs/source/_static/continuous_dems_logo_mini_inf.svg" height="80" alt="Continuous DEMs Logo">
+	</a>
+</p>
 <h1 align="center">Transformez</h1>
 <p align="center"><strong>Global vertical datum transformations, simplified.</strong></p>
 
 <p align="center">
-  <a href="https://github.com/continuous-dems/transformez"><img src="https://img.shields.io/badge/version-0.3.6-blue.svg" alt="Version"></a>
+  <a href="https://github.com/continuous-dems/transformez"><img src="https://img.shields.io/badge/version-0.5.0-blue.svg" alt="Version"></a>
   <a href="LICENSE"><img src="https://img.shields.io/badge/license-MIT-green.svg" alt="License"></a>
   <a href="https://www.python.org/"><img src="https://img.shields.io/badge/python-3.12+-yellow.svg" alt="Python"></a>
   <a href="https://badge.fury.io/py/transformez"><img src="https://badge.fury.io/py/transformez.svg" alt="PyPI version"></a>
@@ -78,4 +78,5 @@ Would you like to know more? Check out our [Official Documentation](https://tran
 ## ⚖ License
 
 This project is licensed under the MIT License - see the [LICENSE](https://github.com/ciresdem/transformez/blob/main/LICENSE) file for details.
+
 Copyright (c) 2010-2026 Regents of the University of Colorado

{transformez-0.4.2 → transformez-0.5.0}/src/transformez/_version.py

@@ -1,4 +1,4 @@
 # This file is auto-generated by Hatchling. As such, do not:
 #   - modify
 #   - track in version control e.g. be sure to add to .gitignore
-__version__ = VERSION = '0.4.2'
+__version__ = VERSION = '0.5.0'

{transformez-0.4.2 → transformez-0.5.0}/src/transformez/api.py

@@ -113,6 +113,8 @@ def generate_grid(
     increment: Union[str, float],
     datum_in: str,
     datum_out: str,
+    epoch_in: str = "2010.0",
+    epoch_out: str = "2010.0",
     decay_pixels: int = 100,
     out_fn: Optional[str] = None,
     cache_dir: Optional[str] = None,
@@ -126,6 +128,8 @@ def generate_grid(
         increment: Resolution (e.g., '3s' or 0.0008333).
         datum_in: Source datum (e.g., 'mllw', '5703').
         datum_out: Target datum (e.g., '4979', '6319').
+        epoch_in: Source epoch (e.g., '2010.0')
+        epoch_out: Target epoch (e.g., '2010.0')
         decay_pixels: Set the pixel decay in case extrapolation is required.
         out_fn: If provided, saves the grid to this file (.tif or .gtx).
         cache_dir: Path to store downloaded grids.
@@ -167,6 +171,8 @@ def generate_grid(
         epsg_out=epsg_out,
         geoid_in=geoid_in,
         geoid_out=geoid_out,
+        epoch_in=epoch_in,
+        epoch_out=epoch_out,
         decay_pixels=decay_pixels,
         cache_dir=cache_dir,
         use_stations=use_stations,
@@ -196,6 +202,7 @@ def transform_raster(
     z_unit_in: Optional[str] = "auto",
     z_unit_out: Optional[str] = "auto",
     use_stations: bool = False,
+    save_shift: bool = False,
     verbose: bool = False,
 ) -> Optional[str]:
     """Apply a vertical datum transformation directly to an existing raster file.
@@ -209,6 +216,8 @@ def transform_raster(
         cache_dir: Path to store downloaded grids.
         z_unit_in: Input DEM z units.
         z_unit_out: Output DEM z units.
+        use_stations: Generate the shift grid from available tide stations,
+        safe_shift: Save the generated shift raster to disk.
         verbose: Enable debug logging.
 
     Returns:
@@ -216,16 +225,48 @@ def transform_raster(
     """
 
     import rasterio
+    from rasterio.warp import transform_bounds, reproject, Resampling
+    from rasterio.transform import from_bounds
+    from fetchez.spatial import Region
 
     if not os.path.exists(input_raster):
         logger.error(f"Input raster not found: {input_raster}")
         return None
 
     with rasterio.open(input_raster) as src:
-        bounds = src.bounds
-        region_obj = Region(bounds.left, bounds.right, bounds.bottom, bounds.top)
+        native_crs = src.crs
+        native_bounds = src.bounds
+        native_transform = src.transform
         nx, ny = src.width, src.height
 
+    is_projected = native_crs.is_projected if native_crs else False
+    if is_projected:
+        logger.info(
+            f"Projected CRS detected ({native_crs}). Extracting WGS84 envelope..."
+        )
+        w, s, e, n = transform_bounds(native_crs, "EPSG:4326", *native_bounds)
+
+        buffer = 0.05
+        region_obj = Region(w - buffer, e + buffer, s - buffer, n + buffer)
+        logger.info(f"Using WGS84 region: {region_obj}")
+
+        inc_deg = 3.0 / 3600.0
+        vt_nx = int((region_obj.xmax - region_obj.xmin) / inc_deg)
+        vt_ny = int((region_obj.ymax - region_obj.ymin) / inc_deg)
+    else:
+        region_obj = Region(
+            native_bounds.left,
+            native_bounds.right,
+            native_bounds.bottom,
+            native_bounds.top,
+        )
+        vt_nx, vt_ny = nx, ny
+
+    # with rasterio.open(input_raster) as src:
+    #     bounds = src.bounds
+    #     region_obj = Region(bounds.left, bounds.right, bounds.bottom, bounds.top)
+    #     nx, ny = src.width, src.height
+
     epsg_in, geoid_in = _parse_datum(datum_in)
     epsg_out, geoid_out = _parse_datum(datum_out)
 
@@ -266,6 +307,47 @@ def transform_raster(
         logger.error("Failed to generate shift array for the raster bounds.")
         return None
 
+    if is_projected:
+        logger.info("Warping shift grid back to native raster projection...")
+        wgs_transform = from_bounds(
+            region_obj.xmin,
+            region_obj.ymin,
+            region_obj.xmax,
+            region_obj.ymax,
+            vt_nx,
+            vt_ny,
+        )
+        native_shift_array = np.zeros((ny, nx), dtype=np.float32)
+
+        reproject(
+            source=shift_array,
+            destination=native_shift_array,
+            src_transform=wgs_transform,
+            src_crs="EPSG:4326",
+            dst_transform=native_transform,
+            dst_crs=native_crs,
+            resampling=Resampling.bilinear,
+        )
+
+        shift_array = native_shift_array
+
+    if save_shift:
+        shift_fn = f"{os.path.splitext(output_raster)[0]}_shiftgrid.tif"
+        logger.info(f"Saving aligned shift grid to {shift_fn}...")
+        with rasterio.open(
+            shift_fn,
+            "w",
+            driver="GTiff",
+            height=ny,
+            width=nx,
+            count=1,
+            dtype=shift_array.dtype,
+            crs=native_crs,
+            transform=native_transform,
+            nodata=-9999.0,
+        ) as dst:
+            dst.write(shift_array, 1)
+
     success = GridEngine.apply_vertical_shift(
         input_raster,
         shift_array,

{transformez-0.4.2 → transformez-0.5.0}/src/transformez/cli.py

@@ -19,6 +19,12 @@ from fetchez.cli import setup_logging
 
 from transformez import api
 
+TRANSFORMEZ_COMMANDS = {
+    "Execution": ["run", "grid", "raster"],
+    "Discovery": ["list"],
+    "External": ["htdp", "vdatum"],
+}
+
 
 class TransformezMainGroup(FetchezMainGroup):
     """A custom Click Group that handles deprecated aliases."""
@@ -68,7 +74,7 @@ class TransformezMainGroup(FetchezMainGroup):
 @click.group(
     name="transform",
     cls=TransformezMainGroup,
-    fetchez_commands=["run", "grid", "raster", "list", "htdp", "vdatum"],
+    fetchez_commands=TRANSFORMEZ_COMMANDS,
 )
 @click.version_option(package_name="transformez")
 @click.option("--verbose", is_flag=True, help="Enable verbose debug logging.")
@@ -296,6 +302,11 @@ def transform_grid(
     is_flag=True,
     help="Force RBF interpolation using live tide stations instead of global satellite models.",
 )
+@click.option(
+    "--save-shift",
+    is_flag=True,
+    help="Save the aligned vertical shift grid to disk alongside the output DEM.",
+)
 def transform_raster(
     input_file,
     input_datum,
@@ -305,6 +316,7 @@ def transform_raster(
     out,
     decay_pixels,
     use_stations,
+    save_shift,
 ):
     """Apply a vertical datum shift to an existing DEM."""
 
@@ -320,6 +332,7 @@ def transform_raster(
         z_unit_in=in_units,
         z_unit_out=out_units,
         use_stations=use_stations,
+        save_shift=save_shift,
         verbose=True,
     )
 
@@ -395,12 +408,17 @@ def htdp_group():
 
 
 @htdp_group.command("install", cls=FetchezMainCommand)
-def install_htdp():
-    """Downloads and compiles the HTDP executable."""
+@click.option(
+    "--version",
+    default="3.5.0",
+    help="HTDP version to install (e.g., 3.3.0, 3.5.0, 3.6.0)",
+)
+def install_htdp(version):
+    """Downloads and compiles the HTDP executable from github."""
 
     from transformez.htdp import install_htdp_binary
 
-    install_htdp_binary()
+    install_htdp_binary(version=version)
 
 
 # --- VDATUM CLI GROUP ---

{transformez-0.4.2 → transformez-0.5.0}/src/transformez/definitions.py

@@ -151,11 +151,11 @@ class Datums:
             "epoch": 1997.0,
         },
         4979: {
-            "name": "WGS_84(original)",
+            "name": "WGS_84(G2139)",
             "description": "(NAD_83(2011) used)",
-            "htdp_id": 4,
+            "htdp_id": 10,
             "uncertainty": 0,
-            "epoch": 1997.0,
+            "epoch": 2020.0,
         },
         7815: {
             "name": "WGS_84(original)",

{transformez-0.4.2 → transformez-0.5.0}/src/transformez/grid_engine.py

@@ -22,9 +22,17 @@ from rasterio.transform import from_bounds
 from scipy import ndimage
 from scipy.interpolate import Rbf
 
+os.environ["HDF5_USE_FILE_LOCKING"] = "FALSE"
+
 logger = logging.getLogger(__name__)
 
 
+class GridCorruptionError(Exception):
+    """Raised when a fetched grid is corrupted and needs re-downloading."""
+
+    pass
+
+
 def plot_grid(grid_array, region, title="Vertical Shift Preview"):
     """Plot the transformation grid using Matplotlib."""
 
@@ -139,6 +147,7 @@ class GridEngine:
 
             try:
                 with rasterio.open(fn) as src:
+                    # logger.info(f"{fn}: {src}")
                     src_data = src.read(1).astype(np.float32)
                     src_nodata = src.nodata
 
@@ -166,8 +175,57 @@ class GridEngine:
                     mosaic[valid_mask] = temp_buffer[valid_mask]
 
             except Exception as e:
-                logger.warning(f"Failed to reproject {fn}: {e}")
+                error_msg = str(e)
+
+                if any(
+                    err in error_msg
+                    for err in [
+                        "-101",
+                        "HDF error",
+                        "not recognized as a supported file format",
+                        "RasterioIOError",
+                    ]
+                ):
+                    logger.error(f" CRITICAL: Corrupted grid chunk detected in {fn}!")
+
+                    real_path = fn.split(":")[1] if fn.startswith("netcdf:") else fn
+                    if os.path.exists(real_path):
+                        logger.warning(
+                            f"Auto-deleting corrupted cache file to force re-fetch: {real_path}"
+                        )
+                        try:
+                            os.remove(real_path)
+                        except OSError:
+                            pass
+
+                    raise GridCorruptionError(f"Corrupted file deleted: {real_path}")
+
+                logger.exception(f"Failed to reproject {fn}: {e}")
                 continue
+            # except Exception as e:
+            #     error_msg = str(e)
+
+            #     if "-101" in error_msg or "HDF error" in error_msg:
+            #         logger.error(f" CRITICAL: Corrupted NetCDF chunk detected in {fn}!")
+
+            #         # Extract the real file path from the GDAL netcdf string
+            #         real_path = fn.split(":")[1] if fn.startswith("netcdf:") else fn
+
+            #         if os.path.exists(real_path):
+            #             logger.warning(f"Auto-deleting corrupted cache file: {real_path}")
+            #             os.remove(real_path)
+
+            #         raise RuntimeError(
+            #             f"Transformation aborted to prevent math corruption. "
+            #             f"The corrupted file has been deleted. Please re-run your command to fetch a fresh copy!"
+            #         )
+
+            #     # For all other normal errors, log and continue as usual
+            #     logger.exception(f"Failed to reproject {fn}: {e}")
+            #     continue
+            # except Exception as e:
+            #     logger.exception(f"Failed to reproject {fn}: {e}")
+            #     continue
 
         # Fill inland areas (decaying to 0) before we clear the remaining NaNs
         # mosaic = GridEngine.fill_nans(mosaic, decay_pixels=decay_pixels)

{transformez-0.4.2 → transformez-0.5.0}/src/transformez/htdp.py

@@ -16,6 +16,7 @@ import sys
 import os
 import subprocess
 import tempfile
+import shutil
 import logging
 import urllib.request
 import zipfile
@@ -40,11 +41,16 @@ HAS_HTDP = utils.cmd_check(f"htdp{ae}", htdp_cmd)
 class HTDP:
     """Wrapper for the NGS HTDP software."""
 
-    def __init__(self, htdp_bin: str = "htdp", verbose: bool = True):
+    def __init__(
+        self, htdp_bin: str = "htdp", version: str = "3.5.0", verbose: bool = True
+    ):
         self.htdp_bin = htdp_bin
+        self.version = version
         self.verbose = verbose
         if not HAS_HTDP:
-            logger.error("HTDP is not installed or not in PATH.")
+            logger.error(
+                f"HTDP {self.version} is not installed or not in PATH. Run 'transformez htdp install --version {self.version}"
+            )
 
     def run_grid(self, region, nx, ny, epsg_in, epsg_out, epoch_in, epoch_out):
         """Main entry point called by transform.py.
@@ -90,8 +96,17 @@ class HTDP:
                 for y_idx, lat in enumerate(lats):
                     for x_idx, lon in enumerate(lons):
                         # "Lat Lon Height TextID"
-                        # We use PNT_x_y tags to robustly map output back to grid
-                        f.write(f'{lat:.9f} {lon:.9f} 0.000 "PNT_{x_idx}_{y_idx}"\n')
+                        # We use PNT_x_y tags to map output back to grid
+                        # htdp_lon = lon * -1.0
+                        if lon < 0:
+                            # Western Hemisphere (e.g., -124.5 -> 124.5 W)
+                            htdp_lon = abs(lon)
+                        else:
+                            # Eastern Hemisphere (e.g., 139.0 E -> 221.0 W)
+                            htdp_lon = 360.0 - lon
+                        f.write(
+                            f'{lat:.9f} {htdp_lon:.9f} 0.000 "PNT_{x_idx}_{y_idx}"\n'
+                        )
 
             # Write Control File
             self._write_control(
@@ -107,7 +122,6 @@ class HTDP:
                 return np.zeros((ny, nx))
 
             grid = self._read_grid(out_fn, (ny, nx))
-
             return grid
 
     def _read_grid(self, filename: str, shape: Tuple[int, int]) -> np.ndarray:
@@ -244,7 +258,61 @@ def download_htdp(target_dir=None):
             logger.error(f"Failed to download or extract HTDP source: {e}")
 
 
-def install_htdp_binary():
+def install_htdp_binary(version="3.5.0"):
+    """Downloads and automatically compiles a specific HTDP version from GitHub."""
+
+    cache_dir = os.path.join(os.getcwd(), "transformez_cache", "bin")
+    os.makedirs(cache_dir, exist_ok=True)
+
+    # Format the GitHub Tag (e.g., "v3.5.0")
+    v_tag = f"v{version}" if not version.startswith("v") else version
+    clean_version = v_tag.replace("v", "")
+
+    if sys.platform == "win32":
+        logger.info(f"Downloading Windows HTDP executable ({v_tag})...")
+        url = f"https://github.com/noaa-ngs/HTDP/releases/download/{v_tag}/htdp.exe"
+        exe_path = os.path.join(cache_dir, f"htdp_{clean_version}.exe")
+
+        try:
+            urllib.request.urlretrieve(url, exe_path)
+            logger.info(f"HTDP {clean_version} installed successfully to: {exe_path}")
+        except Exception as e:
+            logger.error(f"Failed to download {v_tag} from GitHub: {e}")
+
+    else:
+        logger.info(f"Downloading Unix HTDP source code ({v_tag})...")
+        url = f"https://github.com/noaa-ngs/HTDP/archive/refs/tags/{v_tag}.zip"
+        zip_path = os.path.join(cache_dir, f"htdp_{clean_version}.zip")
+
+        try:
+            urllib.request.urlretrieve(url, zip_path)
+
+            extract_dir = os.path.join(cache_dir, f"htdp_src_{clean_version}")
+            with zipfile.ZipFile(zip_path, "r") as z:
+                z.extractall(extract_dir)
+            os.remove(zip_path)
+
+            # GitHub extracts to a folder named Repository-Tag (e.g., HTDP-3.5.0)
+            source_folder = os.path.join(extract_dir, f"HTDP-{clean_version}")
+            fortran_file = os.path.join(source_folder, "htdp.f")
+            out_bin = os.path.join(cache_dir, f"htdp_{clean_version}")
+
+            logger.info(f"Compiling HTDP {clean_version} with gfortran...")
+            subprocess.run(
+                ["gfortran", "-o", out_bin, fortran_file],
+                check=True,
+                capture_output=True,
+            )
+            logger.info(f"HTDP compiled successfully to: {out_bin}")
+
+            # Clean up the raw source files
+            shutil.rmtree(extract_dir)
+
+        except Exception as e:
+            logger.error(f"Failed to install HTDP {clean_version}: {e}")
+
+
+def _install_htdp_binary():
     """Downloads and automatically compiles HTDP."""
 
     cache_dir = os.path.join(os.getcwd(), "transformez_cache", "bin")

{transformez-0.4.2 → transformez-0.5.0}/src/transformez/transform.py

@@ -110,6 +110,7 @@ class VerticalTransform:
             region=self.region.to_list(),
             outdir=self.cache_dir,
             threads=2,
+            check_size=True,
             **kwargs,
         )
 
@@ -206,7 +207,9 @@ class VerticalTransform:
 
         return valid
 
-    def _get_grid(self, provider, name):
+    def _get_grid(self, provider, name, max_retries=3):
+
+        from .grid_engine import GridCorruptionError
 
         if not name:
             return np.zeros((self.ny, self.nx))
@@ -217,98 +220,119 @@ class VerticalTransform:
         if "geoid=" in name:
             name = name.split("=")[1]
 
-        files = self.fetch_grid(provider, datatype=name, query=name)
-        if provider == "vdatum":
-            import rasterio
-            from datetime import datetime
-
-            def get_vdatum_date(gtx_path):
-                """Finds and parses the release date from VDatum metadata files."""
+        for attempt in range(max_retries):
+            files = self.fetch_grid(provider, datatype=name, query=name)
+            if provider == "vdatum":
+                import rasterio
+                from datetime import datetime
 
-                dir_name = os.path.dirname(gtx_path)
-                meta_files = [
-                    f for f in os.listdir(dir_name) if f.endswith((".met", ".inf"))
-                ]
+                def get_vdatum_date(gtx_path):
+                    """Finds and parses the release date from VDatum metadata files."""
 
-                if not meta_files:
-                    return datetime(1970, 1, 1)
-
-                meta_path = os.path.join(dir_name, meta_files[0])
-                try:
-                    with open(meta_path, "r") as f:
-                        content = f.read().splitlines()
+                    dir_name = os.path.dirname(gtx_path)
+                    meta_files = [
+                        f for f in os.listdir(dir_name) if f.endswith((".met", ".inf"))
+                    ]
 
-                    if not content:
+                    if not meta_files:
                         return datetime(1970, 1, 1)
 
-                    first_line = content[0]
-
-                    # Parse the first line "#Mon Jul 08 10:27:07 EDT 2019"
-                    if first_line.startswith("#"):
-                        parts = first_line.replace("#", "").split()
-                        if len(parts) >= 6:
-                            year = int(parts[-1])
-                            day = int(parts[2])
-                            month_map = {
-                                "Jan": 1,
-                                "Feb": 2,
-                                "Mar": 3,
-                                "Apr": 4,
-                                "May": 5,
-                                "Jun": 6,
-                                "Jul": 7,
-                                "Aug": 8,
-                                "Sep": 9,
-                                "Oct": 10,
-                                "Nov": 11,
-                                "Dec": 12,
-                            }
-                            month = month_map.get(parts[1][:3].title(), 1)
-                            return datetime(year, month, day)
-
-                    # Fallback to scanning for "released_date="
-                    for line in content:
-                        if "released_date=" in line:
-                            date_str = line.split("=")[1].strip()
-                            m, d, y = map(int, date_str.split("/"))
-                            return datetime(y, m, d)
+                    meta_path = os.path.join(dir_name, meta_files[0])
+                    try:
+                        with open(meta_path, "r") as f:
+                            content = f.read().splitlines()
+
+                        if not content:
+                            return datetime(1970, 1, 1)
+
+                        first_line = content[0]
+
+                        # Parse the first line "#Mon Jul 08 10:27:07 EDT 2019"
+                        if first_line.startswith("#"):
+                            parts = first_line.replace("#", "").split()
+                            if len(parts) >= 6:
+                                year = int(parts[-1])
+                                day = int(parts[2])
+                                month_map = {
+                                    "Jan": 1,
+                                    "Feb": 2,
+                                    "Mar": 3,
+                                    "Apr": 4,
+                                    "May": 5,
+                                    "Jun": 6,
+                                    "Jul": 7,
+                                    "Aug": 8,
+                                    "Sep": 9,
+                                    "Oct": 10,
+                                    "Nov": 11,
+                                    "Dec": 12,
+                                }
+                                month = month_map.get(parts[1][:3].title(), 1)
+                                return datetime(year, month, day)
+
+                        # Fallback to scanning for "released_date="
+                        for line in content:
+                            if "released_date=" in line:
+                                date_str = line.split("=")[1].strip()
+                                m, d, y = map(int, date_str.split("/"))
+                                return datetime(y, m, d)
+
+                    except Exception as e:
+                        logger.debug(f"Failed to parse date from {meta_path}: {e}")
 
-                except Exception as e:
-                    logger.debug(f"Failed to parse date from {meta_path}: {e}")
+                    return datetime(1970, 1, 1)
 
-                return datetime(1970, 1, 1)
+                def sort_key(filepath):
+                    # Time Sorting (Oldest -> Newest)
+                    date_val = get_vdatum_date(filepath)
 
-            def sort_key(filepath):
-                # Time Sorting (Oldest -> Newest)
-                date_val = get_vdatum_date(filepath)
+                    try:
+                        with rasterio.open(filepath) as src:
+                            b = src.bounds
+                            area = (b.right - b.left) * (b.top - b.bottom)
+                    except Exception:
+                        area = float("inf")
 
-                try:
-                    with rasterio.open(filepath) as src:
-                        b = src.bounds
-                        area = (b.right - b.left) * (b.top - b.bottom)
-                except Exception:
-                    area = float("inf")
+                    return (date_val.timestamp(), -area)
 
-                return (date_val.timestamp(), -area)
+                files.sort(key=sort_key, reverse=True)
 
-            files.sort(key=sort_key, reverse=True)
+            if not files:
+                return np.zeros((self.ny, self.nx))
 
-        if not files:
-            return np.zeros((self.ny, self.nx))
+            try:
+                if provider == "seanoe" or provider == "fes":
+                    var_name = (
+                        "lat_elevation" if "lat" in name.lower() else "msl_elevation"
+                    )
+                    nc_path = f"netcdf:{files[0]}:{var_name}"
+                    return GridEngine.load_and_interpolate(
+                        [nc_path],
+                        self.region,
+                        self.nx,
+                        self.ny,
+                        decay_pixels=self.decay_pixels,
+                    )
 
-        if provider == "seanoe" or provider == "fes":
-            var_name = "lat_elevation" if "lat" in name.lower() else "msl_elevation"
-            nc_path = f"netcdf:{files[0]}:{var_name}"
-            return GridEngine.load_and_interpolate(
-                [nc_path], self.region, self.nx, self.ny, decay_pixels=self.decay_pixels
-            )
+                return GridEngine.load_and_interpolate(
+                    files, self.region, self.nx, self.ny, decay_pixels=self.decay_pixels
+                )
+            except GridCorruptionError:
+                if attempt < max_retries - 1:
+                    logger.warning(
+                        f"Download corruption detected. Retrying fetch (Attempt {attempt + 2}/{max_retries})..."
+                    )
+                    continue
+                else:
+                    logger.error(
+                        "Max retries reached. Could not secure an uncorrupted grid."
+                    )
+                    return np.zeros((self.ny, self.nx))
 
-        return GridEngine.load_and_interpolate(
-            files, self.region, self.nx, self.ny, decay_pixels=self.decay_pixels
-        )
+        return np.zeros((self.ny, self.nx))
 
     def _get_htdp_shift(self, epsg_from, epsg_to, epoch_from, epoch_to):
-        """Calculate Frame Shift via HTDP."""
+        """Calculate Frame Shift via HTDP with Fallback."""
 
         if epsg_from == epsg_to and epoch_from == epoch_to:
             return np.zeros((self.ny, self.nx))
@@ -317,12 +341,41 @@ class VerticalTransform:
 
         try:
             logger.info(f"    [HTDP] Frame Shift: EPSG:{epsg_from} -> EPSG:{epsg_to}")
-            tool = htdp.HTDP(verbose=False)
+            tool = htdp.HTDP(version="3.5.0", verbose=False)
+
+            # Attempt 1: Full Cross-Epoch Shift (Datum Shift + Crustal Velocity)
+            # This seems to fail more often than not!
             grid = tool.run_grid(
                 self.region, self.nx, self.ny, epsg_from, epsg_to, epoch_from, epoch_to
             )
+
+            # --- FALLBACK ---
+            if not np.any(grid):
+                logger.warning(
+                    f"    [HTDP WARNING] Cross-epoch shift failed (likely outside modeled velocity region for {epoch_from} -> {epoch_to})."
+                )
+                logger.warning(
+                    f"    [HTDP WARNING] Falling back to static datum shift at Output Epoch {epoch_to}."
+                )
+
+                # Attempt 2: Static Shift (Datum Shift Only)
+                grid = tool.run_grid(
+                    self.region,
+                    self.nx,
+                    self.ny,
+                    epsg_from,
+                    epsg_to,
+                    epoch_to,
+                    epoch_to,
+                )
+
             if np.any(grid):
                 logger.info(f"    [HTDP] Component Shift (Mean: {np.mean(grid):.3f}m)")
+            else:
+                logger.error(
+                    "    [HTDP FATAL] Both dynamic and static shifts failed returning zeros."
+                )
+
             return grid
 
         except Exception as e:

{transformez-0.4.2 → transformez-0.5.0}/src/transformez/utils.py

@@ -16,6 +16,7 @@ import subprocess
 import logging
 import numpy as np
 import rasterio
+import shutil
 
 logger = logging.getLogger(__name__)
 
@@ -95,3 +96,41 @@ class RasterQuery:
             results[valid] = self.data[rows[valid], cols[valid]]
 
         return results
+
+
+def export_cache(cache_dir=None, output_name="transformez_offline_cache"):
+    """Packs the local transformez cache into a ZIP file for offline use."""
+
+    if cache_dir is None:
+        cache_dir = os.path.join(os.getcwd(), "transformez_cache")
+
+    if not os.path.exists(cache_dir):
+        logger.error(f"[EXPORT FATAL] Cache directory not found at: {cache_dir}")
+        logger.error("Run a transformation to populate the cache before exporting.")
+        return None
+
+    # Determine output path
+    out_path = os.path.abspath(output_name)
+
+    logger.info("-" * 60)
+    logger.info(f"Packing offline cache bundle from: {cache_dir}")
+    logger.info(
+        "This may take a minute depending on the size of your downloaded grids..."
+    )
+
+    try:
+        # shutil.make_archive(base_name, format, root_dir)
+        zip_path = shutil.make_archive(out_path, "zip", cache_dir)
+
+        # Get human-readable file size
+        size_mb = os.path.getsize(zip_path) / (1024 * 1024)
+
+        logger.info(
+            f"Successfully exported offline cache bundle: {zip_path} ({size_mb:.1f} MB)"
+        )
+        logger.info("-" * 60)
+        return zip_path
+
+    except Exception as e:
+        logger.error(f"[EXPORT FATAL] Failed to export cache: {e}")
+        return None