transformez 0.2.0__tar.gz → 0.2.2__tar.gz

transformez-0.2.2/CHANGELOG.md

@@ -0,0 +1,14 @@
+# Changelog
+
+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.1.0] - 2026-02-10
+### Added
+
+### Changed
+- Now uses raserio
+- Renamed project to `transformez`.
+- Refactored and decoupled from old cudem.vdatums

transformez-0.2.2/CITATION.cff

@@ -0,0 +1,19 @@
+cff-version: 1.2.0
+message: "If you use this software, please cite it as below."
+authors:
+  - family-names: Love
+    given-names: Matthew
+  - family-names: Amante
+    given-names: Christopher
+  - family-names: MacFerrin
+    given-names: Michael
+  - family-names: Lim
+    given-names: Elliot
+  - family-names: Fisher
+    given-names: Matt
+  
+website: https://ciresdem.github.io/transformez/
+title: "Transformez"
+version: 0.1.0
+date-released: 2026-02-10
+url: "https://github.com/ciresdem/transformez"

{transformez-0.2.0/src/transformez.egg-info → transformez-0.2.2}/PKG-INFO

@@ -1,8 +1,11 @@
 Metadata-Version: 2.4
 Name: transformez
-Version: 0.2.0
+Version: 0.2.2
 Summary: A standalone utility for vertical elevation datum transformations.
-Author: Matthew Love
+Project-URL: Homepage, https://github.com/continuous-dems/transformez
+Project-URL: Issues, https://github.com/continuous-dems/transformez/issues
+Author-email: Matthew Love <matthew.love@colorado.edu>, Christopher Amante <christopher.amante@colorado.edu>, Elliot Lim <elliot.lim@colorado.edu>, Michael MacFerrin <michael.macferrin@colorado.edu>
+Maintainer-email: Matthew Love <matthew.love@colorado.edu>
 License: MIT License
         
         Copyright (c) 2010-2026 Regents of the University of Colorado
@@ -24,21 +27,19 @@ License: MIT License
         LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
         OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
         SOFTWARE.
-Project-URL: Homepage, https://github.com/ciresdem/transformez
-Project-URL: Issues, https://github.com/ciresdem/transformez/issues
+License-File: AUTHORS.md
+License-File: LICENSE
 Keywords: Geospatial
-Classifier: Programming Language :: Python :: 3
-Classifier: License :: OSI Approved :: MIT License
 Classifier: Operating System :: OS Independent
+Classifier: Programming Language :: Python :: 3
 Classifier: Topic :: Scientific/Engineering :: GIS
-Description-Content-Type: text/markdown
-License-File: LICENSE
-License-File: AUTHORS.md
+Requires-Python: >=3.8
+Requires-Dist: fetchez>0.3.3
 Requires-Dist: numpy<2.0.0
 Requires-Dist: pyproj
 Requires-Dist: rasterio
-Requires-Dist: fetchez>0.3.3
-Dynamic: license-file
+Requires-Dist: scipy
+Description-Content-Type: text/markdown
 
 # 🌍 Transformez ↕
 
@@ -124,4 +125,4 @@ shift, unc = vt._vertical_transform(vt.epsg_in, vt.epsg_out)
 ## 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
+Copyright (c) 2010-2026 Regents of the University of Colorado

transformez-0.2.2/pyproject.toml

@@ -0,0 +1,75 @@
+[build-system]
+requires = ["hatchling>=1.5.0", "hatch-vcs"]
+build-backend = "hatchling.build"
+
+[project]
+name = 'transformez'
+description = 'A standalone utility for vertical elevation datum transformations.'
+readme = "README.md"
+requires-python = ">=3.8"
+license = { file = "LICENSE" }
+dynamic = ["version"]
+authors = [
+    {name = "Matthew Love", email = "matthew.love@colorado.edu"},
+    {name = "Christopher Amante", email = "christopher.amante@colorado.edu"},
+    {name = "Elliot Lim", email = "elliot.lim@colorado.edu"},
+    {name = "Michael MacFerrin", email = "michael.macferrin@colorado.edu"},
+]
+maintainers = [
+    {name = "Matthew Love", email = "matthew.love@colorado.edu"},
+]
+classifiers = [
+    'Programming Language :: Python :: 3',
+    'Operating System :: OS Independent',
+    'Topic :: Scientific/Engineering :: GIS',
+]
+dependencies = [
+    'numpy<2.0.0',
+    'pyproj',
+    'rasterio',
+    'fetchez>0.3.3',
+    'scipy',
+]
+
+keywords = ["Geospatial"]
+
+[dependency-groups]
+test = [
+    "pytest",
+    "pytest-cov",
+]
+
+[project.entry-points."fetchez.plugins"]
+transformez = "transformez"
+
+[project.scripts]
+transformez = 'transformez.cli:transformez_cli'
+
+[project.urls]
+Homepage = "https://github.com/continuous-dems/transformez"
+Issues = "https://github.com/continuous-dems/transformez/issues"
+
+# --- Build configuration ---
+
+[tool.hatch.version]
+source = "vcs"
+
+[tool.hatch.build.hooks.version]
+path = "src/transformez/_version.py"
+
+[tool.hatch.build.targets.wheel]
+packages = ["src/transformez"]
+
+[tool.hatch.build.targets.sdist]
+include = [
+  "src/transformez",
+  "tests",
+  "CHANGELOG.md",
+  "CITATION.cff",
+  "CONTRIBUTING.md",
+  "AUTHORS.md",
+]
+exclude = [
+  ".github",
+]
+

transformez-0.2.2/src/transformez/__init__.py

@@ -0,0 +1,101 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+
+"""
+transformez
+~~~~~~~~~~~~~
+
+:copyright: (c) 2010-2026 Regents of the University of Colorado
+:license: MIT, see LICENSE for more details.
+"""
+
+__author__ = "Matthew Love"
+__credits__ = "CIRES"
+
+try:
+    from transformez._version import __version__
+except ImportError:
+    # Fallback when using the package from source without installing
+    # in editable mode with pip (nobody should do this):
+    # <https://pip.pypa.io/en/stable/topics/local-project-installs/#editable-installs>
+    import warnings
+
+    warnings.warn(
+        "Importing 'transformez' outside a proper installation."
+        " It's highly recommended to install the package from a stable release or"
+        " in editable mode.",
+        stacklevel=2,
+    )
+    __version__ = "dev"
+
+
+import os
+import glob
+# from .hooks import TransformezHook
+# from fetchez.hooks.registry import HookRegistry
+from fetchez.registry import FetchezRegistry
+
+from .modules import TransformezMod
+
+def _find_proj_lib():
+    """Locate the best available PROJ_LIB path."""
+
+    try:
+        import rasterio
+        r_path = os.path.join(os.path.dirname(rasterio.__file__), "proj_data")
+        if os.path.exists(os.path.join(r_path, "proj.db")):
+            return r_path
+
+        parent = os.path.dirname(os.path.dirname(rasterio.__file__))
+        libs = glob.glob(os.path.join(parent, "rasterio.libs*"))
+        if libs:
+            for root, _, files in os.walk(libs[0]):
+                if "proj.db" in files:
+                    return root
+    except ImportError:
+        pass
+
+    try:
+        import pyproj
+        p_path = pyproj.datadir.get_data_dir()
+        if os.path.exists(os.path.join(p_path, "proj.db")):
+            return p_path
+    except ImportError:
+        pass
+
+    return None
+
+target_proj_lib = _find_proj_lib()
+
+if "PROJ_LIB" in os.environ:
+    del os.environ["PROJ_LIB"]
+
+if target_proj_lib:
+    os.environ["PROJ_LIB"] = target_proj_lib
+
+def setup_fetchez(registry_cls):
+    """Called by fetchez when loading plugins.
+
+    Registers modules, hooks, and presets.
+    """
+
+    registry_cls.register_module(
+        'transformez',
+        TransformezMod,
+        metadata={
+            'desc': 'Generate vertical datum shift grids on-demand.',
+            "tags": ["vdatum", "transformation", "shift-grid"],
+            "category": "Tools"
+        }
+    )
+
+    # HookRegistry.register_hook(TransformezHook)
+    # from fetchez.presets import register_global_preset
+    # register_global_preset(
+    #     name="make-shift-grid",
+    #     help_text="Download datum grids and composite them into a single shift grid.",
+    #     hooks=[
+    #         {"name": "transformez", "args": {}}
+    #     ]
+    # )
+setup_fetchez(FetchezRegistry)

transformez-0.2.2/src/transformez/_version.py

@@ -0,0 +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.2.2'

{transformez-0.2.0 → transformez-0.2.2}/src/transformez/cli.py

@@ -73,7 +73,7 @@ def list_supported_datums():
 
 def transformez_cli():
     parser = argparse.ArgumentParser(
-        description="Global Vertical Datum Transformer",
+        description=f"{utils.CYAN}%(prog)s{utils.RESET} ({__version__}) :: Global Vertical Datum Transformer",
         epilog="Examples:\n"
                "  transformez -R -166/-164/63/64 -I mllw -O 4979\n"
                "  transformez input_dem.tif -I mllw -O 5703:geoid=g2012b",
@@ -100,6 +100,9 @@ def transformez_cli():
     grp_sys.add_argument('--list-datums', action='store_true', help='List supported datums.')
     grp_sys.add_argument('--cache-dir', help='Override cache directory.')
     grp_sys.add_argument('--verbose', action='store_true', help='Enable debug logging.')
+    grp_sys.add_argument(
+        "-v", "--version", action="version", version=f"%(prog)s {__version__}"
+    )
 
     args = parser.parse_args()
 

{transformez-0.2.0 → transformez-0.2.2}/src/transformez/grid_engine.py

@@ -165,6 +165,10 @@ class GridWriter:
     def write(filename, data, region):
         """Write a vertical shift grid using Rasterio."""
 
+        dirname = os.path.dirname(filename)
+        if dirname and not os.path.exists(dirname):
+            os.makedirs(dirname)
+
         if not filename.endswith('.tif'):
             filename = os.path.splitext(filename)[0] + '.tif'
 

transformez-0.2.2/src/transformez/modules.py

@@ -0,0 +1,108 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+
+"""
+transformez.modules
+~~~~~~~~~~~~~
+
+Some modules for `fetchez`
+
+:copyright: (c) 2010-2026 Regents of the University of Colorado
+:license: MIT, see LICENSE for more details.
+"""
+
+import os
+import logging
+from fetchez import core, cli
+from transformez.transform import VerticalTransform
+from transformez.grid_engine import GridWriter
+from transformez.definitions import Datums
+
+logger = logging.getLogger(__name__)
+
+
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+
+"""
+transformez.modules.transformez_mod
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+"""
+
+logger = logging.getLogger(__name__)
+
+@cli.cli_opts(
+    help_text="Generate a vertical shift grid (e.g. MLLW to NAVD88).",
+    src_datum="Source Datum (e.g. 'mllw', '5703', '4979').",
+    dst_datum="Destination Datum (e.g. '5703:geoid=g2012b').",
+    increment="Grid resolution (default: 3s).",
+    output_name="Optional output filename override."
+)
+class TransformezMod(core.FetchModule):
+    """A dynamic Fetchez module that generates vertical shift grids on demand.
+
+    Usage:
+      ... transformez --src-datum mllw --dst-datum 5703
+    """
+
+    def __init__(self, src_datum='5703', dst_datum='4979', increment='3s', output_name=None, **kwargs):
+        super().__init__(name="transformez", **kwargs)
+        self.src_datum = src_datum
+        self.dst_datum = dst_datum
+        self.increment = increment
+        self.output_name = output_name
+
+        s_name = str(self.src_datum).replace(':', '_')
+        d_name = str(self.dst_datum).replace(':', '_')
+        w, e, s, n = self.region
+        self.dst_fn = os.path.join(self._outdir, f"shift_{s_name}_to_{d_name}_{w}_{s}.tif")
+
+    def run(self):
+        from fetchez import utils
+        try:
+            inc_val = utils.str2inc(self.increment)
+            nx = int(self.region.width / inc_val)
+            ny = int(self.region.height / inc_val)
+        except Exception:
+            logger.warning(f"Invalid increment '{self.increment}', defaulting to 3s (~0.000833).")
+            # Default roughly 3 arc-seconds (approx 90m)
+            nx = int(self.region.width / 0.00083333333)
+            ny = int(self.region.height / 0.00083333333)
+
+        def parse_d(d_str):
+            if ':' in str(d_str):
+                parts = d_str.split(':')
+                geoid = parts[1].split('=')[1] if 'geoid=' in parts[1] else parts[1]
+                return parts[0], geoid
+            return d_str, None
+
+        epsg_in, geoid_in = parse_d(self.src_datum)
+        epsg_out, geoid_out = parse_d(self.dst_datum)
+
+        vt = VerticalTransform(
+            region=self.region,
+            nx=nx, ny=ny,
+            epsg_in=epsg_in, epsg_out=epsg_out,
+            geoid_in=geoid_in, geoid_out=geoid_out,
+        )
+
+        logger.info(f"Generating shift grid: {self.src_datum} -> {self.dst_datum}...")
+        shift_array, _ = vt._vertical_transform(vt.epsg_in, vt.epsg_out)
+
+        if shift_array is None:
+            logger.error("Transformation failed (No coverage or invalid datums).")
+            return
+
+        GridWriter.write(self.dst_fn, shift_array, self.region)
+
+        self.add_entry_to_results(
+            url=f"file://{self.dst_fn}",
+            dst_fn=self.dst_fn,
+            data_type="gtiff",
+            meta={
+                "src_datum": self.src_datum,
+                "dst_datum": self.dst_datum,
+                "generator": "transformez"
+            }
+        )

{transformez-0.2.0 → transformez-0.2.2}/src/transformez/transform.py

@@ -70,6 +70,7 @@ class VerticalTransform:
 
     def _get_native_ellipsoid(self, epsg, ref_type):
         """Helper to identify the native frame of a datum."""
+
         if ref_type == 'surface':
             # NOAA VDatum = NAD83, Global = WGS84
             region = Datums.SURFACES[epsg].get('region')
@@ -80,7 +81,7 @@ class VerticalTransform:
         elif ref_type == 'htdp':
             # If it's a Frame, it is its own native ellipsoid
             return epsg
-        return WGS84_EPSG # Default fallback
+        return WGS84_EPSG # Default
 
     def fetch_grid(self, module_name, **kwargs):
         """Generic fetcher wrapper."""
@@ -158,34 +159,58 @@ class VerticalTransform:
             logger.error(f"    [HTDP] Failed: {e}")
             return np.zeros((self.ny, self.nx))
 
+    def _fetch_geoid_with_fallback(self, target_geoid):
+        """Fetches a geoid grid. If the primary geoid lacks coverage (e.g., GEOID18 in AK),
+        it automatically falls back to older, compatible models.
+        """
+
+        # Ordered list of preferred US geoids (Newest to Oldest)
+        us_geoids = ['g2018', 'g2012b', 'geoid09']
+
+        if target_geoid in us_geoids:
+            start_idx = us_geoids.index(target_geoid)
+            geoids_to_try = us_geoids[start_idx:]
+        else:
+            geoids_to_try = [target_geoid]
+
+        for g in geoids_to_try:
+            geoid_def = Datums.GEOIDS.get(g, {})
+            provider = geoid_def.get('provider', 'proj')
+            grid = self._get_grid(provider, g)
+
+            if np.any(grid):
+                if g != target_geoid and self.verbose:
+                    logger.info(f"    [Geoid Fallback] '{target_geoid}' lacks coverage here. Falling back to '{g}'.")
+                return grid, g
+
+        return np.zeros((self.ny, self.nx)), target_geoid
+
     # =========================================================================
     # Chains
     # =========================================================================
     def _get_vdatum_chain(self, datum_name, geoid_name):
         """Builds shift: Tidal -> [NAD83 Native]."""
-
         total_shift = np.zeros((self.ny, self.nx))
         desc = []
 
         # Tidal -> LMSL
-        if datum_name not in ["msl", "5714", "lmsl"]:
-            grid = self._get_grid("vdatum", datum_name)
-            if not np.any(grid):
-                return None, f"Missing Tidal Grid: {datum_name}"
-
+        if datum_name not in ['msl', '5714', 'lmsl']:
+            grid = self._get_grid('vdatum', datum_name)
+            if not np.any(grid): return None, f"Missing Tidal Grid: {datum_name}"
             total_shift += grid
             desc.append(f"({datum_name}->LMSL)")
 
         # LMSL -> Ortho (TSS)
-        tss = self._get_grid("vdatum", "tss")
+        tss = self._get_grid('vdatum', 'tss')
         total_shift += tss
         desc.append("TSS(LMSL->NAVD88)")
 
-        # Ortho -> NAD83 (Geoid)
-        actual_geoid = geoid_name if geoid_name else "g2018"
-        geoid = self._get_grid("proj", actual_geoid)
-        total_shift += geoid
-        desc.append(f"Geoid({actual_geoid}->NAD83)")
+        # Ortho -> NAD83 (Smart Geoid Fallback)
+        actual_geoid = geoid_name if geoid_name else 'g2018'
+        geoid_grid, used_geoid = self._fetch_geoid_with_fallback(actual_geoid)
+
+        total_shift += geoid_grid
+        desc.append(f"Geoid({used_geoid}->NAD83)")
 
         return total_shift, " + ".join(desc)
 
@@ -219,8 +244,13 @@ class VerticalTransform:
         mss_name = model_def["grids"].get("mss")
         if mss_name:
             mss_grid = self._get_grid(provider, mss_name)
+            if provider == 'seanoe' or 'fes' in model.lower():
+                mss_grid -= 0.70
+                desc.append("MSS->WGS84(TP_Corr)")
+            else:
+                desc.append("MSS->Ellipsoid")
+
             total_shift += mss_grid
-            desc.append("MSS->Ellipsoid")
 
         if not desc:
             return total_shift, "Global Chain (Empty)"
@@ -232,50 +262,38 @@ class VerticalTransform:
     # =========================================================================
     def _step_to_hub(self, epsg, ref_type, geoid=None, epoch=None):
         shift = np.zeros((self.ny, self.nx))
-        if epsg == self.hub_epsg:
-            return shift, "Already at Hub"
+        if epsg == self.hub_epsg: return shift, "Already at Hub"
 
-        # Determine the Native Ellipsoid of this specific Input
         native_epsg = self._get_native_ellipsoid(epsg, ref_type)
-
-        # Calculate Input -> Native
         chain_shift = None
         chain_desc = ""
 
-        if ref_type == "surface":
-            datum_name = Datums.SURFACES[epsg]["name"]
-            region_tag = Datums.SURFACES[epsg].get("region")
+        if ref_type == 'surface':
+            datum_name = Datums.SURFACES[epsg]['name']
+            region_tag = Datums.SURFACES[epsg].get('region')
 
-            if region_tag == "usa":
+            if region_tag == 'usa':
                 s, d = self._get_vdatum_chain(datum_name, geoid)
-                # Fallback
                 if s is None:
-                    # Switch to Global Chain (Native=WGS84)
                     native_epsg = WGS84_EPSG
                     proxy_name = Datums.get_global_proxy(datum_name)
                     if proxy_name:
-                        s, d = self._get_global_chain(proxy_name, model="fes2014")
-                        d = f"Global({proxy_name}) [Proxy] -> WGS84"
-
+                        s, d = self._get_global_chain(proxy_name, model='fes2014')
+                        # d = f"Global({proxy_name}) [Proxy] -> WGS84"
                 chain_shift, chain_desc = s, d
 
-            elif region_tag == "global":
+            elif region_tag == 'global':
                  chain_shift, chain_desc = self._get_global_chain(datum_name)
 
-        elif ref_type == "cdn":
-            # Ortho -> Native
-            target_geoid = geoid if geoid else "g2018"
-            geoid_def = Datums.GEOIDS.get(target_geoid, {})
-            provider = geoid_def.get("provider", "proj")
-            chain_shift = self._get_grid(provider, target_geoid)
-            chain_desc = f"Ortho(via {target_geoid}) -> Frame({native_epsg})"
+        elif ref_type == 'cdn':
+            target_geoid = geoid if geoid else 'g2018'
+            chain_shift, used_geoid = self._fetch_geoid_with_fallback(target_geoid)
+            chain_desc = f"Ortho(via {used_geoid}) -> Frame({native_epsg})"
 
-        elif ref_type == "htdp":
-            # Frame is already Native
+        elif ref_type == 'htdp':
             chain_shift = np.zeros((self.ny, self.nx))
             chain_desc = f"Frame({epsg})"
 
-        # --- Native -> Hub ---
         if chain_shift is not None:
             if native_epsg != self.hub_epsg:
                 htdp_shift = self._get_htdp_shift(native_epsg, self.hub_epsg, epoch, self.epoch_out)
@@ -287,12 +305,9 @@ class VerticalTransform:
 
     def _step_from_hub(self, epsg, ref_type, geoid=None, epoch=None):
         shift = np.zeros((self.ny, self.nx))
-        if epsg == self.hub_epsg:
-            return shift, "Remain at Hub"
+        if epsg == self.hub_epsg: return shift, "Remain at Hub"
 
         native_epsg = self._get_native_ellipsoid(epsg, ref_type)
-
-        # --- Hub -> Native ---
         total_out = np.zeros((self.ny, self.nx))
         desc_parts = []
 
@@ -301,26 +316,25 @@ class VerticalTransform:
              total_out += htdp_shift
              desc_parts.append(f"Hub({self.hub_epsg}->{native_epsg})")
 
-        # --- Native -> Output ---
-        if ref_type == "surface":
-            datum_name = Datums.SURFACES[epsg]["name"]
-            chain_geoid = geoid if geoid else "g2018"
+        if ref_type == 'surface':
+            datum_name = Datums.SURFACES[epsg]['name']
+            chain_geoid = geoid if geoid else 'g2018'
             s, d = self._get_vdatum_chain(datum_name, chain_geoid)
             if s is None:
                 return np.zeros((self.ny, self.nx)), "FAILED Output Chain"
 
-            # Subtract because Chain is Input->Native
             total_out -= s
             desc_parts.append(f"Native -> VDatum({datum_name})")
 
-        elif ref_type == "cdn":
-            target_geoid = geoid if geoid else "g2018"
-            geoid_def = Datums.GEOIDS.get(target_geoid, {})
-            provider = geoid_def.get("provider", "proj")
-            geoid_grid = self._get_grid(provider, target_geoid)
+        elif ref_type == 'cdn':
+            target_geoid = geoid if geoid else 'g2018'
+            geoid_grid, used_geoid = self._fetch_geoid_with_fallback(target_geoid)
+
+            if not np.any(geoid_grid):
+                logger.warning(f"Geoid {target_geoid} (and fallbacks) not found/covered.")
 
             total_out -= geoid_grid
-            desc_parts.append(f"Native -> Ortho(via {target_geoid})")
+            desc_parts.append(f"Native -> Ortho(via {used_geoid})")
 
         return total_out, " + ".join(desc_parts)
 

transformez-0.2.2/src/transformez/utils.py

@@ -0,0 +1,96 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+
+"""
+transformez.utils
+~~~~~~~~~~~~~
+
+This holds various utility functions.
+
+:copyright: (c) 2010-2026 Regents of the University of Colorado
+:license: MIT, see LICENSE for more details.
+"""
+
+import os
+import subprocess
+import logging
+import numpy as np
+import rasterio
+
+logger = logging.getLogger(__name__)
+cmd_exists = lambda x: any(os.access(os.path.join(path, x), os.X_OK)
+                           for path in os.environ['PATH'].split(os.pathsep))
+
+
+def run_cmd(args):
+    """Standalone replacement for utils.run_cmd using subprocess."""
+
+    logger.info(f"Running: {' '.join(args) if isinstance(args, list) else args}")
+
+    result = subprocess.run(
+        args,
+        shell=False if isinstance(args, list) else True,
+        capture_output=True,
+        text=True
+    )
+    return result.stdout, result.returncode
+
+
+def cmd_check(cmd_str, cmd_vers_str):
+    """check system for availability of 'cmd_str'"""
+
+    if cmd_exists(cmd_str):
+        cmd_vers, status = run_cmd(f"{cmd_vers_str}")
+        return cmd_vers.rstrip()
+    return b"0"
+
+
+class RasterQuery:
+    """Raster query for point clouds.
+    Pre-loads raster data and inverse transform to rapidly query (X, Y) arrays.
+    """
+
+    def __init__(self, filename, default_nodata=0.0):
+        if not filename or not os.path.exists(filename):
+            raise FileNotFoundError(f"Raster not found: {filename}")
+
+        self.default_nodata = default_nodata
+
+        with rasterio.open(filename) as src:
+            self.data = src.read(1)
+            self.transform = src.transform
+            self.inv_transform = ~src.transform
+            self.bounds = src.bounds
+            self.width = src.width
+            self.height = src.height
+
+            if src.nodata is not None:
+                self.data[self.data == src.nodata] = self.default_nodata
+            self.data = np.nan_to_num(self.data, nan=self.default_nodata)
+
+    def query(self, x, y):
+        """query the raster at given X, Y numpy arrays using pure vectorized affine math."""
+
+        q_x = np.asarray(x).copy()
+        q_y = np.asarray(y)
+
+        if self.bounds.left < 0 and np.any(q_x > 180):
+            q_x = np.where(q_x > 180, q_x - 360, q_x)
+        elif self.bounds.left >= 0 and np.any(q_x < 0):
+            q_x = np.where(q_x < 0, q_x + 360, q_x)
+
+        cols_f, rows_f = self.inv_transform * (q_x, q_y)
+
+        cols = np.floor(cols_f).astype(int)
+        rows = np.floor(rows_f).astype(int)
+
+        valid = (
+            (rows >= 0) & (rows < self.height) &
+            (cols >= 0) & (cols < self.width)
+        )
+
+        results = np.full_like(q_x, self.default_nodata, dtype=self.data.dtype)
+        if np.any(valid):
+            results[valid] = self.data[rows[valid], cols[valid]]
+
+        return results

transformez-0.2.2/tests/test_cli.py

@@ -0,0 +1,36 @@
+import subprocess
+import sys
+
+# CMD will run Transformez
+CMD = [sys.executable, "-m", "transformez.cli"]
+
+
+def run_transformez(args):
+    """Run transformez and return result."""
+
+    return subprocess.run(CMD + args, capture_output=True, text=True)
+
+
+def test_help():
+    """Does the help menu work?"""
+
+    result = run_transformez(["--help"])
+    assert result.returncode == 0
+
+
+def test_version():
+    """Does version print?"""
+
+    result = run_transformez(["--version"])
+    assert result.returncode == 0
+
+
+def test_list_modules():
+    """Can we list datums without crashing?"""
+
+    result = run_transformez(["--list-datums"])
+    assert result.returncode == 0
+    assert "lat" in result.stdout
+    assert "mllw" in result.stdout
+    assert "NAVD88" in result.stdout
+    assert "g2018" in result.stdout

transformez-0.2.0/PKG-INFO

@@ -1,127 +0,0 @@
-Metadata-Version: 2.4
-Name: transformez
-Version: 0.2.0
-Summary: A standalone utility for vertical elevation datum transformations.
-Author: Matthew Love
-License: MIT License
-        
-        Copyright (c) 2010-2026 Regents of the University of Colorado
-        
-        Permission is hereby granted, free of charge, to any person obtaining a copy
-        of this software and associated documentation files (the "Software"), to deal
-        in the Software without restriction, including without limitation the rights
-        to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
-        copies of the Software, and to permit persons to whom the Software is
-        furnished to do so, subject to the following conditions:
-        
-        The above copyright notice and this permission notice shall be included in all
-        copies or substantial portions of the Software.
-        
-        THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
-        IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
-        FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
-        AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
-        LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
-        OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
-        SOFTWARE.
-Project-URL: Homepage, https://github.com/ciresdem/transformez
-Project-URL: Issues, https://github.com/ciresdem/transformez/issues
-Keywords: Geospatial
-Classifier: Programming Language :: Python :: 3
-Classifier: License :: OSI Approved :: MIT License
-Classifier: Operating System :: OS Independent
-Classifier: Topic :: Scientific/Engineering :: GIS
-Description-Content-Type: text/markdown
-License-File: LICENSE
-License-File: AUTHORS.md
-Requires-Dist: numpy<2.0.0
-Requires-Dist: pyproj
-Requires-Dist: rasterio
-Requires-Dist: fetchez>0.3.3
-Dynamic: license-file
-
-# 🌍 Transformez ↕
-
-**Global vertical datum transformations, simplified.**
-
-*Transformez Les Données*
-
-> 🚀 **v0.2.0:** Now supporting global tidal transformations via FES2014 & SEANOE.
-
-**Transformez** is a standalone Python engine for converting geospatial data between vertical datums (e.g., `MLLW` ↔ `NAVD88` ↔ `Ellipsoid`).
-
----
-
-## Installation
-
-```bash
-pip install transformez
-```
-
-*Requires [htdp](https://geodesy.noaa.gov/TOOLS/Htdp/Htdp.shtml) to be in your system PATH for frame transformations.*
-
-## Usage
-
-**Generate a vertical shift grid for anywhere on Earth.**
-
-```bash
-# Transform MLLW to WGS84 Ellipsoid in Norton Sound, AK
-# (Where NOAA has no coverage!)
-transformez -R -166/-164/63/64 -E 3s \
-    --input-datum mllw \
-    --output-datum 4979 \
-    --output shift_ak.tif
-```
-
-**Transform a raster directly.** Transformez reads the bounds/resolution from the file.
-
-```bash
-transformez --dem input_bathymetry.tif \
-    --input-datum "mllw" \
-    --output-datum "5703:geoid=geoid12b" \
-    --output output_navd88.tif
-```
-
-**Integrate directly into your download pipeline.**
-
-```bash
-# Download GEBCO and shift EGM96 to WGS84 on the fly
-fetchez gebco ... --hook transformez:datum_in=5773,datum_out=4979
-```
-
-## Python API
-
-```python
-from transformez.transform import VerticalTransform
-from fetchez.spatial import Region
-
-# Define a region in India (Bay of Bengal)
-region = Region(80, 85, 10, 15)
-
-# Initialize Transformer
-# Requesting "MLLW" in India triggers the Global Fallback automatically
-vt = VerticalTransform(
-    region=region,
-    nx=1000, ny=1000,
-    epsg_in="mllw",       # Will resolve to FES2014 LAT
-    epsg_out="epsg:4979"  # WGS84 Ellipsoid
-)
-
-# Generate Shift
-shift, unc = vt._vertical_transform(vt.epsg_in, vt.epsg_out)
-```
-
-## Supported Datums
-
-* **Tidal**: mllw, mhhw, msl, lat
-
-* **Ellipsoidal**: 4979 (WGS84), 6319 (NAD83 2011)
-
-* **Orthometric**: 5703 (NAVD88), egm2008, egm96
-
-* **Geoids**: g2018, g2012b, geoid09, xgeoid20b
-
-## 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.2.0/pyproject.toml

@@ -1,41 +0,0 @@
-[build-system]
-requires = ['setuptools>=61.0']
-build-backend = 'setuptools.build_meta'
-
-[project]
-name = 'transformez'
-version = '0.2.0'
-description = 'A standalone utility for vertical elevation datum transformations.'
-readme = "README.md"
-authors = [{ name = 'Matthew Love' }]
-license = { file = 'LICENSE' }
-classifiers = [
-    'Programming Language :: Python :: 3',
-    'License :: OSI Approved :: MIT License',
-    'Operating System :: OS Independent',
-    'Topic :: Scientific/Engineering :: GIS',
-]
-dependencies = [
-    'numpy<2.0.0',
-    'pyproj',
-    'rasterio',
-    'fetchez>0.3.3',
-]
-
-keywords = ["Geospatial"]
-
-[project.entry-points."fetchez.plugins"]
-transformez = "transformez"
-
-[project.scripts]
-transformez = 'transformez.cli:transformez_cli'
-
-[project.urls]
-Homepage = "https://github.com/ciresdem/transformez"
-Issues = "https://github.com/ciresdem/transformez/issues"
-
-[tool.setuptools.packages.find]
-where = ["src"]  # Look for packages inside the 'src' folder
-include = ["transformez*"]  # Include the package named 'transformez'
-exclude = ["tests*", "docs*", "scripts*", "examples*"]
-namespaces = false

transformez-0.2.0/setup.cfg

@@ -1,4 +0,0 @@
-[egg_info]
-tag_build = 
-tag_date = 0
-

transformez-0.2.0/setup.py

@@ -1,2 +0,0 @@
-from setuptools import setup
-setup()

transformez-0.2.0/src/transformez/__init__.py

@@ -1,75 +0,0 @@
-#!/usr/bin/env python
-# -*- coding: utf-8 -*-
-
-"""
-transformez
-~~~~~~~~~~~~~
-
-:copyright: (c) 2010-2026 Regents of the University of Colorado
-:license: MIT, see LICENSE for more details.
-"""
-
-__version__ = "0.2.0"
-__author__ = "Matthew Love"
-__credits__ = "CIRES"
-
-import os
-import glob
-from .hooks import TransformezHook
-from fetchez.hooks.registry import HookRegistry
-# from fetchez.registry import FetchezRegistry
-
-def _find_proj_lib():
-    """Locate the best available PROJ_LIB path."""
-
-    try:
-        import rasterio
-        r_path = os.path.join(os.path.dirname(rasterio.__file__), "proj_data")
-        if os.path.exists(os.path.join(r_path, "proj.db")):
-            return r_path
-
-        parent = os.path.dirname(os.path.dirname(rasterio.__file__))
-        libs = glob.glob(os.path.join(parent, "rasterio.libs*"))
-        if libs:
-            for root, _, files in os.walk(libs[0]):
-                if "proj.db" in files:
-                    return root
-    except ImportError:
-        pass
-
-    try:
-        import pyproj
-        p_path = pyproj.datadir.get_data_dir()
-        if os.path.exists(os.path.join(p_path, "proj.db")):
-            return p_path
-    except ImportError:
-        pass
-
-    return None
-
-target_proj_lib = _find_proj_lib()
-
-if "PROJ_LIB" in os.environ:
-    del os.environ["PROJ_LIB"]
-
-if target_proj_lib:
-    os.environ["PROJ_LIB"] = target_proj_lib
-
-def setup_fetchez(registry_cls):
-    """Called by fetchez when loading plugins.
-
-    Registers modules, hooks, and presets.
-    """
-
-    HookRegistry.register_hook(TransformezHook)
-
-    from fetchez.presets import register_global_preset
-    register_global_preset(
-        name="make-shift-grid",
-        help_text="Download datum grids and composite them into a single shift grid.",
-        hooks=[
-            {"name": "transformez", "args": {}}
-        ]
-    )
-
-# setup_fetchez(FetchezRegistry)

transformez-0.2.0/src/transformez/utils.py

@@ -1,43 +0,0 @@
-#!/usr/bin/env python
-# -*- coding: utf-8 -*-
-
-"""
-transformez.utils
-~~~~~~~~~~~~~
-
-This holds various utility functions.
-
-:copyright: (c) 2010-2026 Regents of the University of Colorado
-:license: MIT, see LICENSE for more details.
-"""
-
-import os
-import subprocess
-import logging
-
-logger = logging.getLogger(__name__)
-cmd_exists = lambda x: any(os.access(os.path.join(path, x), os.X_OK)
-                           for path in os.environ['PATH'].split(os.pathsep))
-
-
-def run_cmd(args):
-    """Standalone replacement for utils.run_cmd using subprocess."""
-
-    logger.info(f"Running: {' '.join(args) if isinstance(args, list) else args}")
-
-    result = subprocess.run(
-        args,
-        shell=False if isinstance(args, list) else True,
-        capture_output=True,
-        text=True
-    )
-    return result.stdout, result.returncode
-
-
-def cmd_check(cmd_str, cmd_vers_str):
-    """check system for availability of 'cmd_str'"""
-
-    if cmd_exists(cmd_str):
-        cmd_vers, status = run_cmd(f"{cmd_vers_str}")
-        return cmd_vers.rstrip()
-    return b"0"

transformez-0.2.0/src/transformez.egg-info/SOURCES.txt

@@ -1,22 +0,0 @@
-AUTHORS.md
-LICENSE
-README.md
-pyproject.toml
-setup.py
-src/transformez/__init__.py
-src/transformez/cli.py
-src/transformez/definitions.py
-src/transformez/grid_engine.py
-src/transformez/hooks.py
-src/transformez/htdp.py
-src/transformez/spatial.py
-src/transformez/srs.py
-src/transformez/transform.py
-src/transformez/utils.py
-src/transformez/vdatum.py
-src/transformez.egg-info/PKG-INFO
-src/transformez.egg-info/SOURCES.txt
-src/transformez.egg-info/dependency_links.txt
-src/transformez.egg-info/entry_points.txt
-src/transformez.egg-info/requires.txt
-src/transformez.egg-info/top_level.txt

transformez-0.2.0/src/transformez.egg-info/dependency_links.txt

@@ -1 +0,0 @@
-

transformez-0.2.0/src/transformez.egg-info/entry_points.txt

@@ -1,5 +0,0 @@
-[console_scripts]
-transformez = transformez.cli:transformez_cli
-
-[fetchez.plugins]
-transformez = transformez

transformez-0.2.0/src/transformez.egg-info/requires.txt

@@ -1,4 +0,0 @@
-numpy<2.0.0
-pyproj
-rasterio
-fetchez>0.3.3

transformez-0.2.0/src/transformez.egg-info/top_level.txt

@@ -1 +0,0 @@
-transformez