transformez 0.1.0__py3-none-any.whl

transformez/__init__.py

@@ -0,0 +1,125 @@
+# -*- coding: utf-8 -*-
+
+__version__ = "0.1.0"
+__author__ = "Matthew Love"
+__credits__ = "CIRES"
+
+import os
+
+# --- fix the PROJ_LIB path to work with rasterio/pyproj
+def _find_proj_lib():
+    """Locate the best available PROJ_LIB path."""
+    
+    try:
+        import rasterio
+        # Common path in wheels: site-packages/rasterio/proj_data
+        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
+            
+        # Linux wheels often put it in .libs adjacent to the package
+        # e.g. site-packages/rasterio.libs/proj.db
+        parent = os.path.dirname(os.path.dirname(rasterio.__file__))
+        libs = glob.glob(os.path.join(parent, 'rasterio.libs*'))
+        if libs:
+            # Look inside the libs folder
+            for root, _, files in os.walk(libs[0]):
+                if 'proj.db' in files:
+                    return root
+    except ImportError:
+        pass
+
+    # Try PyProj's bundled data (Fallback)
+    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()
+
+# Unset conflicting system/conda vars if we found a better one
+if 'PROJ_LIB' in os.environ:
+    del os.environ['PROJ_LIB']
+
+if target_proj_lib:
+    os.environ['PROJ_LIB'] = target_proj_lib
+    # print(f"DEBUG: PROJ_LIB set to {target_proj_lib}")
+
+# # --- PROJ_LIB Fix (Run before any geo-imports) ---
+# # This prevents "PROJ: proj_create_from_database: Cannot find proj.db" errors
+# # when conflicting Conda/System PROJ installations exist.
+# try:
+#     # 1. We must import pyproj first to let it find its own bundled data
+#     import os
+#     import pyproj
+    
+#     # 2. Get the valid data directory from the wheel
+#     proj_lib = pyproj.datadir.get_data_dir()
+#     #print(proj_lib)
+#     # 3. Force the environment variable to use this valid path
+#     #    (Overriding any bad global/Conda defaults)
+#     os.environ['PROJ_LIB'] = proj_lib
+    
+# except ImportError:
+#     # If pyproj isn't installed yet (e.g. during pip install), skip this.
+#     pass
+
+# --- End PROJ_LIB Fix ---
+
+from .hooks import TransformezHook
+#from .modules import DatumGridFetcher
+from fetchez.hooks.registry import HookRegistry
+# form fetchez.registry import FetchezRegistry # Implicitly handled via setup_fetchez passing cls
+
+def setup_fetchez(registry_cls):
+    """Called by fetchez when loading plugins.
+    Registers modules, hooks, and presets.
+    """
+
+    # module should gather necessary grids to do the transformation
+    # Register the Module with Fetchez
+    # registry_cls.register_module(
+    #     mod_key='datum_grids',
+    #     mod_cls=DatumGridFetcher,
+    #     metadata={
+    #         'desc': 'Fetch NOAA VDatum and PROJ-CDN grids for a region',
+    #         'tags': ['vdatum', 'geoid', 'transformation']
+    #     }
+    # )
+
+    # Register Fetchez Hooks
+    HookRegistry.register_hook(TransformezHook)
+    
+    # Register Global Presets
+    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 GTX shift grid.",
+        hooks=[
+            {"name": "transformez", "args": {}} 
+        ]
+    )
+
+
+    # "transform-pipeline": {
+    #     "help_text": "Generate shift grid based on region, then apply it to files.",
+    #     "hooks": [
+    #         {
+    #             "name": "transformez", 
+    #             "args": {"stage": "pre", "datum_in": "5703", "output_grid": "/tmp/shift.gtx"}
+    #         },
+    #         {
+    #             "name": "transformez", 
+    #             "args": {"stage": "file", "apply": "True", "output_grid": "/tmp/shift.gtx"}
+    #         },
+    #         {
+    #              "name": "audit" 
+    #         }
+    #     ]
+    # }

transformez/cli.py

@@ -0,0 +1,228 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+
+"""
+transformez.cli
+~~~~~~~~~~~~~
+
+The transformez CLI.
+Generates vertical transformation grids using the Fetchez-DLIM ecosystem.
+
+:copyright: (c) 2010-2026 Regents of the University of Colorado
+:license: MIT, see LICENSE for more details.
+"""
+
+import os
+import sys
+import argparse
+import logging
+
+import pyproj
+import rasterio
+
+from . import __version__
+from .transform import VerticalTransform
+from .definitions import Datums
+from .grid_engine import plot_grid, GridWriter, GridEngine
+
+import fetchez
+from fetchez import spatial
+from fetchez import utils
+from fetchez.spatial import parse_region, Region
+
+logging.basicConfig(level=logging.INFO, format='[ %(levelname)s ] %(name)s: %(message)s', stream=sys.stderr)
+logger = logging.getLogger(__name__)
+
+logging.getLogger('fetchez').setLevel(logging.WARNING)
+
+def parse_compound_datum(datum_arg):
+    """Parse a datum string that might contain a geoid override.
+    Format: "EPSG" or "EPSG:GEOID" or "NAME:GEOID"
+    """
+    if ':' in str(datum_arg):
+        parts = str(datum_arg).split(':')
+        datum = Datums.get_vdatum_by_name(parts[0])
+        geoid = parts[1]
+        return datum, geoid
+    else:
+        return Datums.get_vdatum_by_name(datum_arg), None
+
+    
+def get_grid_info(filename):
+    """Extract region, resolution, and SRS from a raster using Rasterio."""
+    
+    with rasterio.open(filename) as ds:
+        bounds = ds.bounds # left, bottom, right, top
+        width = ds.width
+        height = ds.height
+        gt = ds.transform.to_gdal() # (c, a, b, f, d, e)
+        srs_wkt = ds.crs.to_wkt() if ds.crs else None
+        
+        return {
+            'te': (bounds.left, bounds.bottom, bounds.right, bounds.top),
+            'region': (bounds.left, bounds.right, bounds.bottom, bounds.top),
+            'nx': width,
+            'ny': height,
+            'gt': gt,
+            'srs_wkt': srs_wkt
+        }
+
+
+def transformez_cli():
+    parser = argparse.ArgumentParser(
+        description=f'%(prog)s ({__version__}): Generate a vertical transformation grid',
+        formatter_class=argparse.RawTextHelpFormatter,
+        epilog="CUDEM home page: <http://cudem.colorado.edu>"
+    )
+
+    input_grp = parser.add_mutually_exclusive_group(required=True)
+    input_grp.add_argument('-R', '--region', help=spatial.region_help_msg())
+    input_grp.add_argument('--dem', help='Input DEM to transform. Automatically sets Region and Resolution.')
+    
+    sel_grp = parser.add_argument_group('Geospatial Selection')
+    sel_grp.add_argument('-E', '--increment', help='Grid resolution (e.g. 0.0001 or 1s) (Required if not using --dem).')
+
+    datum_group = parser.add_argument_group('Datum Configuration')    
+    datum_group.add_argument('-I', '--vdatum_in', default='5703', 
+                        help='Input vertical datum. Format: "EPSG" or "EPSG:GEOID" (e.g. "5703:g2012a")')
+    datum_group.add_argument('-O', '--vdatum_out', default='7662', 
+                        help='Output vertical datum. Format: "EPSG" or "EPSG:GEOID"')
+    datum_group.add_argument('--epoch-in', type=float, default=1997.0, 
+                        help='Input coordinate epoch (decimal year).')
+    datum_group.add_argument('--epoch-out', type=float, default=1997.0, 
+                        help='Output coordinate epoch (decimal year).')
+
+    proc_group = parser.add_argument_group('Processing Options')    
+    proc_group.add_argument('--preview', action='store_true', help='Plot the transformation grid (matplotlib) before processing.')
+    proc_group.add_argument('--output', help='Output filename (default auto-generated).')
+    
+    sys_group = parser.add_argument_group('System & Logging')    
+    sys_group.add_argument('-D', '--cache-dir', help='Directory for storing temporary grids.')
+    sys_group.add_argument('-k', '--keep-cache', action='store_true', help='Do not delete temporary files after run.')
+    sys_group.add_argument('-l', '--list-epsg', action='store_true', help='List supported EPSG codes/names and exit.')
+    sys_group.add_argument('-q', '--quiet', action='store_true', help='Suppress log output.')
+    sys_group.add_argument('--version', action='version', version=f'%(prog)s {__version__}')
+
+    # Handle negative coordinates in arguments
+    fixed_argv = spatial.fix_argparse_region(sys.argv[1:])
+    args = parser.parse_args(fixed_argv)
+
+    if args.quiet:
+        logger.setLevel(logging.WARNING)
+
+    if args.list_epsg:
+        def _print_epsg(title, data):
+            print(f'{title}:')
+            for key, val in data.items():
+                print(f'  {key}\t{val["name"]}')
+        
+        _print_epsg('HTDP EPSG', Datums.HTDP)
+        _print_epsg('CDN EPSG', Datums.CDN)
+        _print_epsg('Tidal EPSG', Datums.TIDAL)
+        sys.exit(0)
+
+    cache_dir = args.cache_dir or os.path.join(os.path.expanduser('~'), '.transformez')
+    if not os.path.exists(cache_dir):
+        os.makedirs(cache_dir)
+
+    epsg_in, geoid_in = parse_compound_datum(args.vdatum_in)
+    epsg_out, geoid_out = parse_compound_datum(args.vdatum_out)
+
+    if args.dem:
+        if not os.path.exists(args.dem):
+            logger.error(f'Input DEM not found: {args.dem}')
+            sys.exit(1)
+            
+        logger.info(f'Extracting grid info from DEM: {args.dem}')
+        info = get_grid_info(args.dem)
+        logger.info(f'Grid info is: {info}')
+        
+        # Override region and dimensions from DEM
+        # info['te'] is (w, s, e, n)
+        region_obj = spatial.Region.from_list(info['region'])
+        nx = info['nx']
+        ny = info['ny']
+
+        if info.get('srs_wkt'):
+            pass
+        
+        # Default Output Name for DEM
+        if not args.output:
+            base, ext = os.path.splitext(args.dem)
+            dst_grid = f'{base}_trans_{epsg_out}{ext}'
+        else:
+            dst_grid = args.output
+            
+    elif args.region:
+        # parse_region returns a list of Region objects
+        # we're only using one for now.
+        these_regions = parse_region(args.region)
+        region_obj = these_regions[0]
+
+        try:
+            if '/' in args.increment:
+                inc_x, inc_y = [float(x) for x in args.increment.split('/')]
+            else:
+                inc_val = utils.str2inc(args.increment)
+                inc_x, inc_y = inc_val, inc_val
+                
+            width = region_obj.width
+            height = region_obj.height
+            nx = int(width / inc_x)
+            ny = int(height / inc_y)
+        except Exception as e:
+            logger.error(f'Invalid increment: {args.increment}. {e}')
+            sys.exit(1)
+        
+        if not args.output:
+            base, ext = 'transformez_trans', 'tif' # Default to TIF now
+            v_out_str = str(epsg_out)
+            if geoid_out:
+                v_out_str += f"_{geoid_out}"
+            dst_grid = f"{base}_{epsg_in}_{v_out_str}.{ext}"
+        else:
+            dst_grid = args.output                
+    else:
+        logger.error("Region or DEM is required.")
+        sys.exit(1)
+            
+    # Initialize Vertical Transform
+    vt = VerticalTransform(
+        region=region_obj, 
+        nx=nx,
+        ny=ny,
+        epsg_in=epsg_in, 
+        epsg_out=epsg_out,
+        geoid_in=geoid_in,
+        geoid_out=geoid_out,
+        epoch_in=args.epoch_in,
+        epoch_out=args.epoch_out,
+        cache_dir=cache_dir,
+    )
+    
+    logger.info(f"Generating shift grid: {epsg_in} -> {epsg_out}")
+    shift_array, _ = vt._vertical_transform(vt.epsg_in, vt.epsg_out)
+    
+    if shift_array is not None:
+        if args.preview:
+            plot_grid(
+                shift_array, 
+                region=region_obj, 
+                title=f"Shift: {epsg_in} -> {epsg_out}"
+            )
+
+        if args.dem:
+            # Apply to DEM
+            logger.info(f"Applying transformation to DEM...")
+            GridEngine.apply_vertical_shift(args.dem, shift_array, dst_grid)
+        else:
+            # Just write the grid
+            logger.info(f"Saving transformation grid to: {dst_grid}")            
+            GridWriter.write(dst_grid, shift_array, region_obj)            
+        return            
+    else:
+        logger.error("Failed to generate transformation grid.")
+
+        
+if __name__ == '__main__':
+    transformez_cli()

transformez/definitions.py

@@ -0,0 +1,187 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+
+"""
+transformez.definitions
+~~~~~~~~~~~~~
+
+This file contains the various vertical datum transformation references 
+and definitions.
+
+:copyright: (c) 2010-2026 Regents of the University of Colorado
+:license: MIT, see LICENSE for more details.
+"""
+
+import logging
+
+logger = logging.getLogger(__name__)
+
+class Datums:
+    """Class to manage vertical datum definitions and lookups."""
+
+    # =========================================================================
+    # Vertical Datum References
+    # =========================================================================
+    SURFACES = {
+        # --- Tidal Datums ---
+        1089: {'name': 'mllw', 'description': 'Mean Lower Low Water', 'uncertainty': 0, 'epsg': 5866},
+        5866: {'name': 'mllw', 'description': 'Mean Lower Low Water', 'uncertainty': 0, 'epsg': 5866},
+        1091: {'name': 'mlw',  'description': 'Mean Low Water', 'uncertainty': 0, 'epsg': 1091},
+        5869: {'name': 'mhhw', 'description': 'Mean Higher High Water', 'uncertainty': 0, 'epsg': 5869},
+        5868: {'name': 'mhw', 'description': 'Mean High Water', 'uncertainty': 0, 'epsg': 5868},
+        5714: {'name': 'msl', 'description': 'Mean Sea Level', 'uncertainty': 0, 'epsg': 5714},
+        5713: {'name': 'mtl', 'description': 'Mean Tide Level', 'uncertainty': 0, 'epsg': 5713},
+        
+        # --- Hydraulic / River Datums ---
+        # Columbia River Datum (No standard EPSG, using 0 placeholder or custom)
+        0:    {'name': 'crd', 'description': 'Columbia River Datum', 'uncertainty': 0, 'epsg': 0},
+        
+        # IGLD 1985 (Dynamic Height)
+        5609: {'name': 'IGLD85', 'description': 'International Great Lakes Datum 1985', 'uncertainty': 0, 'epsg': 5609},
+        
+        # IGLD Low Water Datum (Chart Datum for Lakes)
+        # VDatum uses 'LWD_IGLD85' string
+        9000: {'name': 'LWD_IGLD85', 'description': 'IGLD85 Low Water Datum', 'uncertainty': 0, 'epsg': 5609},
+        
+        # --- Legacy Vertical ---
+        # NGVD29 is often best handled via VDatum (VERTCON) if PROJ isn't configured
+        5702: {'name': 'NGVD29', 'description': 'National Geodetic Vertical Datum 1929', 'uncertainty': 0.05, 'epsg': 5702},
+    }
+
+    HTDP = {
+        4269: {'name': 'NAD_83(2011/CORS96/2007)', 'description': '(North American plate fixed)', 'htdp_id': 1, 'uncertainty': .02, 'epoch': 1997.0},
+        6781: {'name': 'NAD_83(2011/CORS96/2007)', 'description': '(North American plate fixed)', 'htdp_id': 1, 'uncertainty': .02, 'epoch': 1997.0},
+        6319: {'name': 'NAD_83(2011/CORS96/2007)', 'description': '(North American plate fixed)', 'htdp_id': 1, 'uncertainty': .02, 'epoch': 1997.0},
+        6321: {'name': 'NAD_83(PA11/PACP00)', 'description': '(Pacific plate fixed)', 'htdp_id': 2, 'uncertainty': .02, 'epoch': 1997.0},
+        6324: {'name': 'NAD_83(MA11/MARP00)', 'description': '(Mariana plate fixed)', 'htdp_id': 3, 'uncertainty': .02, 'epoch': 1997.0},
+        4979: {'name': 'WGS_84(original)', 'description': '(NAD_83(2011) used)', 'htdp_id': 4, 'uncertainty': 0, 'epoch': 1997.0},
+        7815: {'name': 'WGS_84(original)', 'description': '(NAD_83(2011) used)', 'htdp_id': 4, 'uncertainty': 0, 'epoch': 1997.0},
+        7816: {'name': 'WGS_84(original)', 'description': '(NAD_83(2011) used)', 'htdp_id': 4, 'uncertainty': 0, 'epoch': 1997.0},
+        7656: {'name': 'WGS_84(G730)', 'description': '(ITRF91 used)', 'htdp_id': 5, 'uncertainty': 0, 'epoch': 1997.0},
+        7657: {'name': 'WGS_84(G730)', 'description': '(ITRF91 used)', 'htdp_id': 5, 'uncertainty': 0, 'epoch': 1997.0},
+        7658: {'name': 'WGS_84(G873)', 'description': '(ITRF94 used)', 'htdp_id': 6, 'uncertainty': 0, 'epoch': 1997.0},
+        7659: {'name': 'WGS_84(G873)', 'description': '(ITRF94 used)', 'htdp_id': 6, 'uncertainty': 0, 'epoch': 1997.0},
+        7660: {'name': 'WGS_84(G1150)', 'description': '(ITRF2000 used)', 'htdp_id': 7, 'uncertainty': 0, 'epoch': 1997.0},
+        7661: {'name': 'WGS_84(G1150)', 'description': '(ITRF2000 used)', 'htdp_id': 7, 'uncertainty': 0, 'epoch': 1997.0},
+        7662: {'name': 'WGS_84(G1674)', 'description': '(ITRF2008 used)', 'htdp_id': 8, 'uncertainty': 0, 'epoch': 2000.0},
+        7663: {'name': 'WGS_84(G1674)', 'description': '(ITRF2008 used)', 'htdp_id': 8, 'uncertainty': 0, 'epoch': 2000.0},
+        7664: {'name': 'WGS_84(G1762)', 'description': '(IGb08 used)', 'htdp_id': 9, 'uncertainty': 0, 'epoch': 2000.0},
+        7665: {'name': 'WGS_84(G1762)', 'description': '(IGb08 used)', 'htdp_id': 9, 'uncertainty': 0, 'epoch': 2000.0},
+        7666: {'name': 'WGS_84(G2139)', 'description': '(ITRF2014=IGS14=IGb14 used)', 'htdp_id': 10, 'uncertainty': 0, 'epoch': 1997.0},
+        7667: {'name': 'WGS_84(G2139)', 'description': '(ITRF2014=IGS14=IGb14 used)', 'htdp_id': 10, 'uncertainty': 0, 'epoch': 1997.0},
+        4910: {'name': 'ITRF88', 'description': '', 'htdp_id': 11, 'uncertainty': 0, 'epoch': 1988.0},
+        4911: {'name': 'ITRF89', 'description': '', 'htdp_id': 12, 'uncertainty': 0, 'epoch': 1988.0},
+        7901: {'name': 'ITRF89', 'description': '', 'htdp_id': 12, 'uncertainty': 0, 'epoch': 1988.0},
+        7902: {'name': 'ITRF90', 'description': '(PNEOS90/NEOS90)', 'htdp_id': 13, 'uncertainty': 0, 'epoch': 1988.0},
+        7903: {'name': 'ITRF91', 'description': '', 'htdp_id': 14, 'uncertainty': 0, 'epoch': 1988.0},
+        7904: {'name': 'ITRF92', 'description': '', 'htdp_id': 15, 'uncertainty': 0, 'epoch': 1988.0},
+        7905: {'name': 'ITRF93', 'description': '', 'htdp_id': 16, 'uncertainty': 0, 'epoch': 1988.0},
+        7906: {'name': 'ITRF94', 'description': '', 'htdp_id': 17, 'uncertainty': 0, 'epoch': 1988.0},
+        7907: {'name': 'ITRF96', 'description': '', 'htdp_id': 18, 'uncertainty': 0, 'epoch': 1996.0},
+        7908: {'name': 'ITRF97', 'description': 'IGS97', 'htdp_id': 19, 'uncertainty': 0, 'epoch': 1997.0},
+        7909: {'name': 'ITRF2000', 'description': 'IGS00/IGb00', 'htdp_id': 20, 'uncertainty': 0, 'epoch': 2000.0},
+        7910: {'name': 'ITRF2005', 'description': 'IGS05', 'htdp_id': 21, 'uncertainty': 0, 'epoch': 2000.0},
+        7911: {'name': 'ITRF2008', 'description': 'IGS08/IGb08', 'htdp_id': 22, 'uncertainty': 0, 'epoch': 2000.0},
+        7912: {'name': 'ELLIPSOID', 'description': 'IGS14/IGb14/WGS84/ITRF2014 Ellipsoid', 'htdp_id': 23, 'uncertainty': 0, 'epoch': 2000.0},
+        1322: {'name': 'ITRF2020', 'description': 'IGS20', 'htdp_id': 24, 'uncertainty': 0, 'epoch': 2000.0},
+        
+    }
+
+    CDN = {
+        # CONUS / Alaska / Hawaii / PR / VI
+        5703: {'name': 'NAVD88 height', 'vdatum_id': 'navd88:m:height', 'default_geoid': 'g2018', 'ellipsoid': 6319},
+        6360: {'name': 'NAVD88 height (usFt)', 'default_geoid': 'g2018'},
+        8228: {'name': 'NAVD88 height (Ft)', 'default_geoid': 'g2018'},
+
+        # Puerto Rico
+        6641: {'name': 'PRVD02 height', 'vdatum_id': 'prvd02:m:height', 'default_geoid': 'g2018', 'ellipsoid': 6319},
+        
+        # Virgin Islands
+        6642: {'name': 'VIVD09 height', 'vdatum_id': 'vivd09:m:height', 'default_geoid': 'g2018', 'ellipsoid': 6319},
+
+        # Canada (CGVD2013 uses CGG2013 geoid)
+        # Note: You need to ensure 'CGG2013' is fetchable via your fetcher or map it to a filename
+        6647: {'name': 'CGVD2013(CGG2013)', 'vdatum_id': 'cgvd2013:m:height', 'default_geoid': 'CGG2013'},
+
+        # Global EGM
+        3855: {'name': 'EGM2008 height', 'vdatum_id': 'egm2008:m:height', 'default_geoid': 'egm2008'},
+        5773: {'name': 'EGM96 height', 'vdatum_id': 'egm96:m:height', 'default_geoid': 'egm96'},
+        
+        # # Ellipsoidal (Hubs) - No Geoid needed
+        # 6319: {'name': 'NAD83(2011)', 'vdatum_id': 'nad83_2011:m:height'},
+        # 4979: {'name': 'WGS84', 'vdatum_id': 'wgs84:m:height'},
+    }
+    
+    GEOIDS = {
+        # Standard PROJ-CDN Geoids (Default provider is 'proj')
+        'g2018':   {'name': 'geoid 2018', 'uncertainty': .0127, 'provider': 'proj'},
+        'g2012b':  {'name': 'geoid 2012b', 'uncertainty': .017,  'provider': 'proj'},
+        'geoid09': {'name': 'geoid 2009', 'uncertainty': .05,   'provider': 'proj'},
+        
+        # New XGEOIDs via VDatum (Provider is 'vdatum')
+        'xgeoid20b': {'name': 'xgeoid20b', 'uncertainty': .02, 'provider': 'vdatum'},
+        'xgeoid19b': {'name': 'xgeoid19b', 'uncertainty': .02, 'provider': 'vdatum'},
+
+        'egm2008': {'name': 'EGM2008', 'uncertainty': 0, 'provider': 'proj'},
+        'egm96':   {'name': 'EGM96',   'uncertainty': 0, 'provider': 'proj'},
+
+        'CGG2013': {'name': 'CGG2013', 'uncertainty': 0.01, 'provider': 'proj'}, 
+    }
+
+    @classmethod
+    def get_default_geoid(cls, epsg):
+        """Return default geoid for a generic CDN EPSG, or None."""
+        
+        try:
+            e_int = int(epsg)
+        except:
+            return None
+            
+        if e_int in cls.CDN:
+            return cls.CDN[e_int].get('default_geoid')
+        return None
+
+    
+    @classmethod
+    def get_vdatum_by_name(cls, datum_name):
+        """Return the vertical datum EPSG based on the vertical datum name."""
+        
+        if datum_name is None:
+            return None
+        
+        try:
+            datum_int = int(datum_name)
+        except (ValueError, TypeError):
+            datum_int = None
+        
+        for frame_set in [cls.SURFACES, cls.HTDP, cls.CDN]:
+            if datum_int in frame_set:
+                return datum_int
+            
+            for epsg, info in frame_set.items():
+                if str(datum_name).lower() in info['name'].lower():
+                    return epsg
+        
+        return None
+
+
+    @classmethod
+    def get_vdatum_id(cls, epsg):
+        """Retrieve the NOAA VDatum CLI string for an EPSG."""
+        if epsg in cls.SURFACES: return cls.SURFACES[epsg].get('vdatum_id')
+        if epsg in cls.CDN: return cls.CDN[epsg].get('vdatum_id')
+        if epsg == 6319: return 'nad83_2011:m:height'
+        return None
+
+    
+    @classmethod
+    def get_frame_type(cls, epsg):
+        """Identify which frame set an EPSG belongs to."""
+        
+        if epsg in cls.SURFACES:
+            return 'surface'
+        if epsg in cls.HTDP:
+            return 'htdp'
+        if epsg in cls.CDN:
+            return 'cdn'
+        return None