demtools 0.1.0__tar.gz

demtools-0.1.0/MANIFEST.in

@@ -0,0 +1,3 @@
+include LICENSE
+include README.md
+include requirements.txt

demtools-0.1.0/PKG-INFO

@@ -0,0 +1,86 @@
+Metadata-Version: 2.4
+Name: demtools
+Version: 0.1.0
+Summary: A collection of tools for working with DEM (Digital Elevation Model) raster files
+Home-page: https://github.com/AaronOET/demtools
+Author: aaronchh
+Author-email: aaronchh <aaronhsu219@gmail.com>
+License: MIT
+Project-URL: Homepage, https://github.com/AaronOET/demtools
+Classifier: Programming Language :: Python :: 3
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Operating System :: OS Independent
+Requires-Python: >=3.8
+Description-Content-Type: text/markdown
+Requires-Dist: numpy>=1.20.0
+Dynamic: author
+Dynamic: home-page
+Dynamic: requires-python
+
+# DEMTOOLS
+
+A collection of Python tools for working with DEM (Digital Elevation Model) raster files.
+
+> **GDAL Installation**: GDAL is required for this package. For conda environments, use `conda install gdal` to install GDAL. For non-conda environments, download the appropriate wheel file from [https://github.com/cgohlke/geospatial-wheels/releases](https://github.com/cgohlke/geospatial-wheels/releases) to install GDAL.
+
+## Installation
+
+```bash
+pip install -e .
+```
+
+## Features
+
+- **chgnodata**: Convert (replace) the nodata value of GeoTIFF files, updating pixel data accordingly
+- **setnodata**: Assign a nodata value to GeoTIFF files without modifying pixel data (metadata only)
+- **csv2tif**: Convert a plain-numeric CSV raster grid to a GeoTIFF file
+
+## Usage
+
+### Command-line
+
+```bash
+# List all available commands
+demtools-info
+
+# Convert nodata values for all TIF files in the current directory
+chgnodata -a
+
+# Convert nodata of a single file to -9999
+chgnodata -i dem.tif -v -9999
+
+# Set nodata metadata for all TIF files (no pixel data change)
+setnodata -a -v -9999
+
+# Convert a CSV grid to GeoTIFF
+csv2tif -i grid.csv -o dem.tif --xll 250000 --yll 2500000 --cellsize 5 --epsg 32648
+```
+
+### Python API
+
+```python
+from demtools import chgnodata, setnodata, csv2tif
+
+# Convert nodata value in a raster file
+chgnodata.convert_nodata_value("dem.tif", output_nodata=-9999)
+
+# Assign nodata metadata without touching pixel values
+setnodata.define_nodata_value("dem.tif", nodata_value=-9999)
+
+# Convert CSV grid to GeoTIFF
+csv2tif.csv_to_tif(
+    input_csv="grid.csv",
+    output_tif="dem.tif",
+    xll=250000,
+    yll=2500000,
+    cellsize=5,
+    nodata=-999,
+    epsg=32648,
+)
+```
+
+## Notes
+
+- `chgnodata` and `setnodata` automatically create a `RAS_BAK/` directory with backup copies before modifying files.
+- `csv2tif` expects a plain numeric CSV (no headers), with rows ordered from north to south.
+- GDAL must be installed separately via conda or a pre-built wheel; it is not listed in `requirements.txt` as it cannot be reliably installed via pip on all platforms.

demtools-0.1.0/README.md

@@ -0,0 +1,67 @@
+# DEMTOOLS
+
+A collection of Python tools for working with DEM (Digital Elevation Model) raster files.
+
+> **GDAL Installation**: GDAL is required for this package. For conda environments, use `conda install gdal` to install GDAL. For non-conda environments, download the appropriate wheel file from [https://github.com/cgohlke/geospatial-wheels/releases](https://github.com/cgohlke/geospatial-wheels/releases) to install GDAL.
+
+## Installation
+
+```bash
+pip install -e .
+```
+
+## Features
+
+- **chgnodata**: Convert (replace) the nodata value of GeoTIFF files, updating pixel data accordingly
+- **setnodata**: Assign a nodata value to GeoTIFF files without modifying pixel data (metadata only)
+- **csv2tif**: Convert a plain-numeric CSV raster grid to a GeoTIFF file
+
+## Usage
+
+### Command-line
+
+```bash
+# List all available commands
+demtools-info
+
+# Convert nodata values for all TIF files in the current directory
+chgnodata -a
+
+# Convert nodata of a single file to -9999
+chgnodata -i dem.tif -v -9999
+
+# Set nodata metadata for all TIF files (no pixel data change)
+setnodata -a -v -9999
+
+# Convert a CSV grid to GeoTIFF
+csv2tif -i grid.csv -o dem.tif --xll 250000 --yll 2500000 --cellsize 5 --epsg 32648
+```
+
+### Python API
+
+```python
+from demtools import chgnodata, setnodata, csv2tif
+
+# Convert nodata value in a raster file
+chgnodata.convert_nodata_value("dem.tif", output_nodata=-9999)
+
+# Assign nodata metadata without touching pixel values
+setnodata.define_nodata_value("dem.tif", nodata_value=-9999)
+
+# Convert CSV grid to GeoTIFF
+csv2tif.csv_to_tif(
+    input_csv="grid.csv",
+    output_tif="dem.tif",
+    xll=250000,
+    yll=2500000,
+    cellsize=5,
+    nodata=-999,
+    epsg=32648,
+)
+```
+
+## Notes
+
+- `chgnodata` and `setnodata` automatically create a `RAS_BAK/` directory with backup copies before modifying files.
+- `csv2tif` expects a plain numeric CSV (no headers), with rows ordered from north to south.
+- GDAL must be installed separately via conda or a pre-built wheel; it is not listed in `requirements.txt` as it cannot be reliably installed via pip on all platforms.

demtools-0.1.0/demtools/__init__.py

@@ -0,0 +1,17 @@
+"""
+DEMTOOLS - A collection of tools for working with DEM (Digital Elevation Model) raster files.
+"""
+
+__version__ = '0.1.0'
+
+__all__ = [
+    'chgnodata',
+    'csv2tif',
+    'setnodata',
+    'describe',
+]
+
+from . import chgnodata
+from . import csv2tif
+from . import setnodata
+from . import describe

demtools-0.1.0/demtools/chgnodata.py

@@ -0,0 +1,184 @@
+#!/usr/bin/env python3
+"""
+chgnodata - Convert nodata values for TIF raster files
+Uses GDAL/OGR Python bindings
+"""
+
+import glob
+import argparse
+import os
+import sys
+from osgeo import gdal, gdalconst
+
+def convert_nodata_value(input_file, output_nodata=-999):
+    """
+    Convert nodata value of a raster file to specified value
+
+    Args:
+        input_file (str): Path to input TIF file
+        output_nodata (float): New nodata value (default: -999)
+    """
+    print(f"Processing: {input_file}")
+
+    # Open the input dataset
+    src_ds = gdal.Open(input_file, gdalconst.GA_ReadOnly)
+    if src_ds is None:
+        print(f"Error: Could not open {input_file}")
+        return False
+
+    # Get raster band (assuming single band, modify if multi-band)
+    src_band = src_ds.GetRasterBand(1)
+
+    # Get current nodata value
+    current_nodata = src_band.GetNoDataValue()
+    print(f"  Current nodata value: {current_nodata}")
+
+    if current_nodata == output_nodata:
+        print(f"  Nodata value already set to {output_nodata}, skipping...")
+        src_ds = None
+        return True
+
+    # Read the data
+    data = src_band.ReadAsArray()
+
+    # Get raster properties
+    rows, cols = data.shape
+    geotransform = src_ds.GetGeoTransform()
+    projection = src_ds.GetProjection()
+    data_type = src_band.DataType
+
+    # Create backup directory
+    backup_dir = "RAS_BAK"
+    if not os.path.exists(backup_dir):
+        os.makedirs(backup_dir)
+        print(f"  Created backup directory: {backup_dir}")
+
+    backup_filename = os.path.basename(input_file)
+    backup_file = os.path.join(backup_dir, backup_filename)
+
+    # Create backup
+    print(f"  Creating backup: {backup_file}")
+    gdal.GetDriverByName('GTiff').CreateCopy(backup_file, src_ds)
+
+    # Close source dataset
+    src_ds = None
+
+    # Create new dataset (overwrite original)
+    driver = gdal.GetDriverByName('GTiff')
+    out_ds = driver.Create(input_file, cols, rows, 1, data_type)
+
+    # Set geotransform and projection
+    out_ds.SetGeoTransform(geotransform)
+    out_ds.SetProjection(projection)
+
+    # Get output band
+    out_band = out_ds.GetRasterBand(1)
+
+    # Replace current nodata values with new nodata value
+    if current_nodata is not None:
+        data[data == current_nodata] = output_nodata
+
+    # Write the data
+    out_band.WriteArray(data)
+
+    # Set the new nodata value
+    out_band.SetNoDataValue(output_nodata)
+
+    # Flush data to disk
+    out_band.FlushCache()
+    out_ds = None
+
+    print(f"  Successfully converted nodata value to {output_nodata}")
+    return True
+
+
+def main():
+    parser = argparse.ArgumentParser(
+        prog='chgnodata',
+        description='Convert nodata values for TIF raster files',
+        formatter_class=argparse.RawDescriptionHelpFormatter,
+        epilog="""
+Examples:
+  chgnodata -a                       # Process all *.tif in current directory (nodata=-999)
+  chgnodata -a -v -9999              # Process all *.tif, set nodata to -9999
+  chgnodata -i dem.tif               # Process a single file (nodata=-999)
+  chgnodata -i dem.tif -v -9999      # Process a single file, set nodata to -9999
+        """
+    )
+
+    # Mutually exclusive group: -a or -i
+    group = parser.add_mutually_exclusive_group(required=True)
+    group.add_argument(
+        '-a', '--all',
+        action='store_true',
+        help='Process all *.tif files in the current directory'
+    )
+    group.add_argument(
+        '-i', '--input',
+        metavar='FILE',
+        help='Process a single TIF file'
+    )
+
+    parser.add_argument(
+        '-v', '--value',
+        dest='nodata_value',
+        type=float,
+        default=-999,
+        metavar='NODATA',
+        help='New nodata value to set (default: -999)'
+    )
+
+    args = parser.parse_args()
+    new_nodata = args.nodata_value
+
+    print(f"chgnodata: converting nodata values to {new_nodata}")
+    print("=" * 60)
+
+    if args.all:
+        tif_files = glob.glob("*.tif")
+        if not tif_files:
+            print("No TIF files found in current directory.")
+            sys.exit(0)
+
+        print(f"Found {len(tif_files)} TIF file(s):")
+        for f in tif_files:
+            print(f"  - {f}")
+        print(f"\nProcessing files...")
+        print("-" * 40)
+
+        success_count = 0
+        for tif_file in tif_files:
+            if convert_nodata_value(tif_file, new_nodata):
+                success_count += 1
+            print()
+
+        print("Conversion complete!")
+        print(f"Successfully processed: {success_count}/{len(tif_files)} files")
+        if success_count < len(tif_files):
+            print("Some files failed to process. Check error messages above.")
+
+    elif args.input:
+        input_file = args.input
+        if not os.path.isfile(input_file):
+            print(f"Error: File not found: {input_file}")
+            sys.exit(1)
+
+        print("-" * 40)
+        success = convert_nodata_value(input_file, new_nodata)
+        print()
+        if success:
+            print("Conversion complete!")
+        else:
+            print("Conversion failed.")
+            sys.exit(1)
+
+
+if __name__ == "__main__":
+    gdal.UseExceptions()
+
+    try:
+        main()
+    except Exception as e:
+        print(f"Error: {e}")
+        import traceback
+        traceback.print_exc()

demtools-0.1.0/demtools/csv2tif.py

@@ -0,0 +1,157 @@
+#!/usr/bin/env python3
+"""
+csv2tif - Convert a CSV raster grid to a GeoTIFF file
+Uses GDAL/OGR Python bindings
+
+Usage:
+    csv2tif -i input.csv
+    csv2tif -i input.csv -o output.tif
+    csv2tif -i input.csv -o output.tif --xll 0.0 --yll 0.0 --cellsize 1.0
+    csv2tif -i input.csv --epsg 32648
+"""
+
+import argparse
+import os
+import sys
+
+import numpy as np
+from osgeo import gdal, osr
+
+
+def csv_to_tif(input_csv, output_tif, xll=0.0, yll=0.0, cellsize=1.0,
+               nodata=-999, epsg=None):
+    """
+    Convert a CSV raster grid to a GeoTIFF.
+
+    The CSV is expected to be a plain numeric grid (no headers) where rows
+    correspond to raster rows from top (north) to bottom (south).
+
+    Args:
+        input_csv  (str):   Path to input CSV file.
+        output_tif (str):   Path to output GeoTIFF file.
+        xll        (float): X coordinate of the lower-left corner (default 0.0).
+        yll        (float): Y coordinate of the lower-left corner (default 0.0).
+        cellsize   (float): Cell/pixel size in map units (default 1.0).
+        nodata     (float): NoData value to assign (default -999).
+        epsg       (int):   EPSG code for the spatial reference (optional).
+    """
+    print(f"Input CSV  : {input_csv}")
+    print(f"Output TIF : {output_tif}")
+
+    # --- Read CSV -------------------------------------------------------
+    try:
+        data = np.genfromtxt(input_csv, delimiter=",", dtype=np.float32)
+    except Exception as exc:
+        print(f"Error reading CSV: {exc}")
+        sys.exit(1)
+
+    if data.ndim == 1:
+        # Single-row CSV — treat as one-row raster
+        data = data.reshape(1, -1)
+
+    rows, cols = data.shape
+    print(f"Grid size  : {cols} cols x {rows} rows")
+
+    # --- Build GeoTransform ---------------------------------------------
+    # GDAL geotransform: (x_top_left, pixel_width, 0, y_top_left, 0, -pixel_height)
+    x_top_left = xll
+    y_top_left = yll + rows * cellsize
+    geotransform = (x_top_left, cellsize, 0.0, y_top_left, 0.0, -cellsize)
+
+    # --- Create output GeoTIFF -----------------------------------------
+    driver = gdal.GetDriverByName("GTiff")
+    out_ds = driver.Create(output_tif, cols, rows, 1, gdal.GDT_Float32)
+    if out_ds is None:
+        print(f"Error: Could not create {output_tif}")
+        sys.exit(1)
+
+    out_ds.SetGeoTransform(geotransform)
+
+    # Spatial reference
+    srs = osr.SpatialReference()
+    if epsg is not None:
+        srs.ImportFromEPSG(epsg)
+        out_ds.SetProjection(srs.ExportToWkt())
+        print(f"Projection : EPSG:{epsg}")
+    else:
+        print("Projection : none (no EPSG provided)")
+
+    # Write band
+    out_band = out_ds.GetRasterBand(1)
+    out_band.SetNoDataValue(nodata)
+    out_band.WriteArray(data)
+    out_band.FlushCache()
+
+    out_ds.FlushCache()
+    out_ds = None
+
+    print(f"NoData     : {nodata}")
+    print(f"Done -> {output_tif}")
+
+
+def main():
+    parser = argparse.ArgumentParser(
+        prog='csv2tif',
+        description="Convert a CSV raster grid to a GeoTIFF",
+        formatter_class=argparse.RawDescriptionHelpFormatter,
+        epilog="""
+Examples:
+  csv2tif -i input.csv
+  csv2tif -i input.csv -o output.tif
+  csv2tif -i input.csv -o output.tif --xll 250000 --yll 2500000 --cellsize 5
+  csv2tif -i input.csv --epsg 32648
+        """
+    )
+    parser.add_argument(
+        "-i", "--input", required=True,
+        help="Input CSV file path"
+    )
+    parser.add_argument(
+        "-o", "--output", default=None,
+        help="Output TIF file path (default: same name as input with .tif extension)"
+    )
+    parser.add_argument(
+        "--xll", type=float, default=0.0,
+        help="X coordinate of lower-left corner (default: 0.0)"
+    )
+    parser.add_argument(
+        "--yll", type=float, default=0.0,
+        help="Y coordinate of lower-left corner (default: 0.0)"
+    )
+    parser.add_argument(
+        "--cellsize", type=float, default=1.0,
+        help="Cell size in map units (default: 1.0)"
+    )
+    parser.add_argument(
+        "-n", "--nodata", type=float, default=-999,
+        help="NoData value to assign (default: -999)"
+    )
+    parser.add_argument(
+        "--epsg", type=int, default=None,
+        help="EPSG code for spatial reference (optional)"
+    )
+
+    args = parser.parse_args()
+
+    if not os.path.isfile(args.input):
+        print(f"Error: Input file not found: {args.input}")
+        sys.exit(1)
+
+    output_tif = args.output
+    if output_tif is None:
+        base = os.path.splitext(args.input)[0]
+        output_tif = base + ".tif"
+
+    csv_to_tif(
+        input_csv=args.input,
+        output_tif=output_tif,
+        xll=args.xll,
+        yll=args.yll,
+        cellsize=args.cellsize,
+        nodata=args.nodata,
+        epsg=args.epsg,
+    )
+
+
+if __name__ == "__main__":
+    main()

demtools-0.1.0/demtools/describe.py

@@ -0,0 +1,78 @@
+#!/usr/bin/env python
+"""
+Command-line utility to display descriptions of demtools functionality.
+"""
+
+import argparse
+import sys
+from textwrap import dedent
+
+TOOL_DESCRIPTIONS = {
+    'chgnodata': """
+        Convert nodata values for TIF raster files.
+
+        This tool reads one or more GeoTIFF files, replaces the existing nodata
+        value with a new one (re-writing pixel data), and saves a backup of the
+        original file in the RAS_BAK/ directory.
+
+        Examples:
+            chgnodata -a                  # Process all *.tif (nodata=-999)
+            chgnodata -a -v -9999         # Process all *.tif, set nodata to -9999
+            chgnodata -i dem.tif          # Process a single file (nodata=-999)
+            chgnodata -i dem.tif -v -9999 # Process a single file, set nodata to -9999
+    """,
+    'setnodata': """
+        Define a nodata value for TIF raster files without modifying pixel data.
+
+        This tool updates only the nodata *metadata* of one or more GeoTIFF files
+        without touching the actual pixel values. A backup is saved in RAS_BAK/.
+
+        Examples:
+            setnodata -a                  # Process all *.tif (nodata=-999)
+            setnodata -a -v -9999         # Process all *.tif, define nodata as -9999
+            setnodata -i dem.tif          # Process a single file (nodata=-999)
+            setnodata -i dem.tif -v -9999 # Process a single file, define nodata as -9999
+    """,
+    'csv2tif': """
+        Convert a plain-numeric CSV raster grid to a GeoTIFF file.
+
+        The CSV must contain only numeric values arranged in rows (top→bottom) and
+        columns (left→right), with no header row. Georeferencing parameters
+        (lower-left corner coordinates, cell size, EPSG) can be provided via options.
+
+        Examples:
+            csv2tif -i grid.csv
+            csv2tif -i grid.csv -o dem.tif
+            csv2tif -i grid.csv --xll 250000 --yll 2500000 --cellsize 5 --epsg 32648
+            csv2tif -i grid.csv -n -9999
+    """,
+}
+
+
+def main():
+    parser = argparse.ArgumentParser(
+        prog='demtools-info',
+        description='Display descriptions of demtools commands',
+        formatter_class=argparse.RawDescriptionHelpFormatter,
+    )
+    parser.add_argument(
+        'tool',
+        nargs='?',
+        choices=list(TOOL_DESCRIPTIONS.keys()),
+        help='Tool name to describe (omit to list all tools)',
+    )
+    args = parser.parse_args()
+
+    if args.tool:
+        print(f"\n--- {args.tool} ---")
+        print(dedent(TOOL_DESCRIPTIONS[args.tool]))
+    else:
+        print("\nDEMTOOLS — DEM raster utility commands\n")
+        for tool, desc in TOOL_DESCRIPTIONS.items():
+            first_line = dedent(desc).strip().splitlines()[0]
+            print(f"  {tool:<14} {first_line}")
+        print("\nRun 'demtools-info <tool>' for details on a specific command.")
+
+
+if __name__ == "__main__":
+    main()

demtools-0.1.0/demtools/setnodata.py

@@ -0,0 +1,164 @@
+#!/usr/bin/env python3
+"""
+setnodata - Define (assign) a nodata value for TIF raster files without modifying pixel data
+Uses GDAL/OGR Python bindings
+"""
+
+import glob
+import argparse
+import os
+import sys
+from osgeo import gdal, gdalconst
+
+def define_nodata_value(input_file, nodata_value=-999):
+    """
+    Define a nodata value for a raster file without changing any pixel data.
+    Only the nodata metadata is updated in the file.
+
+    Args:
+        input_file (str): Path to input TIF file
+        nodata_value (float): Nodata value to assign (default: -999)
+    """
+    print(f"Processing: {input_file}")
+
+    # Open the input dataset in read-only to inspect current state
+    src_ds = gdal.Open(input_file, gdalconst.GA_ReadOnly)
+    if src_ds is None:
+        print(f"Error: Could not open {input_file}")
+        return False
+
+    # Get raster band (assuming single band, modify if multi-band)
+    src_band = src_ds.GetRasterBand(1)
+
+    # Get current nodata value
+    current_nodata = src_band.GetNoDataValue()
+    print(f"  Current nodata value: {current_nodata}")
+
+    if current_nodata == nodata_value:
+        print(f"  Nodata value already set to {nodata_value}, skipping...")
+        src_ds = None
+        return True
+
+    # Close read-only handle before backup/update
+    src_ds = None
+
+    # Create backup directory
+    backup_dir = "RAS_BAK"
+    if not os.path.exists(backup_dir):
+        os.makedirs(backup_dir)
+        print(f"  Created backup directory: {backup_dir}")
+
+    backup_filename = os.path.basename(input_file)
+    backup_file = os.path.join(backup_dir, backup_filename)
+
+    # Create backup before modifying
+    print(f"  Creating backup: {backup_file}")
+    src_ds = gdal.Open(input_file, gdalconst.GA_ReadOnly)
+    gdal.GetDriverByName('GTiff').CreateCopy(backup_file, src_ds)
+    src_ds = None
+
+    # Re-open in update mode to set nodata metadata only (no pixel data changed)
+    ds = gdal.Open(input_file, gdalconst.GA_Update)
+    if ds is None:
+        print(f"Error: Could not open {input_file} for update")
+        return False
+
+    band = ds.GetRasterBand(1)
+    band.SetNoDataValue(nodata_value)
+    band.FlushCache()
+    ds = None
+
+    print(f"  Successfully defined nodata value as {nodata_value}")
+    return True
+
+
+def main():
+    parser = argparse.ArgumentParser(
+        prog='setnodata',
+        description='Define a nodata value for TIF raster files without modifying pixel data',
+        formatter_class=argparse.RawDescriptionHelpFormatter,
+        epilog="""
+Examples:
+  setnodata -a                       # Process all *.tif in current directory (nodata=-999)
+  setnodata -a -v -9999              # Process all *.tif, define nodata as -9999
+  setnodata -i dem.tif               # Process a single file (nodata=-999)
+  setnodata -i dem.tif -v -9999      # Process a single file, define nodata as -9999
+        """
+    )
+
+    # Mutually exclusive group: -a or -i
+    group = parser.add_mutually_exclusive_group(required=True)
+    group.add_argument(
+        '-a', '--all',
+        action='store_true',
+        help='Process all *.tif files in the current directory'
+    )
+    group.add_argument(
+        '-i', '--input',
+        metavar='FILE',
+        help='Process a single TIF file'
+    )
+
+    parser.add_argument(
+        '-v', '--value',
+        dest='nodata_value',
+        type=float,
+        default=-999,
+        metavar='NODATA',
+        help='Nodata value to assign (default: -999)'
+    )
+
+    args = parser.parse_args()
+    new_nodata = args.nodata_value
+
+    print(f"setnodata: defining nodata value as {new_nodata}")
+    print("=" * 60)
+
+    if args.all:
+        tif_files = glob.glob("*.tif")
+        if not tif_files:
+            print("No TIF files found in current directory.")
+            sys.exit(0)
+
+        print(f"Found {len(tif_files)} TIF file(s):")
+        for f in tif_files:
+            print(f"  - {f}")
+        print(f"\nProcessing files...")
+        print("-" * 40)
+
+        success_count = 0
+        for tif_file in tif_files:
+            if define_nodata_value(tif_file, new_nodata):
+                success_count += 1
+            print()
+
+        print("Done!")
+        print(f"Successfully processed: {success_count}/{len(tif_files)} files")
+        if success_count < len(tif_files):
+            print("Some files failed to process. Check error messages above.")
+
+    elif args.input:
+        input_file = args.input
+        if not os.path.isfile(input_file):
+            print(f"Error: File not found: {input_file}")
+            sys.exit(1)
+
+        print("-" * 40)
+        success = define_nodata_value(input_file, new_nodata)
+        print()
+        if success:
+            print("Done!")
+        else:
+            print("Failed.")
+            sys.exit(1)
+
+
+if __name__ == "__main__":
+    gdal.UseExceptions()
+
+    try:
+        main()
+    except Exception as e:
+        print(f"Error: {e}")
+        import traceback
+        traceback.print_exc()

demtools-0.1.0/demtools.egg-info/PKG-INFO

@@ -0,0 +1,86 @@
+Metadata-Version: 2.4
+Name: demtools
+Version: 0.1.0
+Summary: A collection of tools for working with DEM (Digital Elevation Model) raster files
+Home-page: https://github.com/AaronOET/demtools
+Author: aaronchh
+Author-email: aaronchh <aaronhsu219@gmail.com>
+License: MIT
+Project-URL: Homepage, https://github.com/AaronOET/demtools
+Classifier: Programming Language :: Python :: 3
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Operating System :: OS Independent
+Requires-Python: >=3.8
+Description-Content-Type: text/markdown
+Requires-Dist: numpy>=1.20.0
+Dynamic: author
+Dynamic: home-page
+Dynamic: requires-python
+
+# DEMTOOLS
+
+A collection of Python tools for working with DEM (Digital Elevation Model) raster files.
+
+> **GDAL Installation**: GDAL is required for this package. For conda environments, use `conda install gdal` to install GDAL. For non-conda environments, download the appropriate wheel file from [https://github.com/cgohlke/geospatial-wheels/releases](https://github.com/cgohlke/geospatial-wheels/releases) to install GDAL.
+
+## Installation
+
+```bash
+pip install -e .
+```
+
+## Features
+
+- **chgnodata**: Convert (replace) the nodata value of GeoTIFF files, updating pixel data accordingly
+- **setnodata**: Assign a nodata value to GeoTIFF files without modifying pixel data (metadata only)
+- **csv2tif**: Convert a plain-numeric CSV raster grid to a GeoTIFF file
+
+## Usage
+
+### Command-line
+
+```bash
+# List all available commands
+demtools-info
+
+# Convert nodata values for all TIF files in the current directory
+chgnodata -a
+
+# Convert nodata of a single file to -9999
+chgnodata -i dem.tif -v -9999
+
+# Set nodata metadata for all TIF files (no pixel data change)
+setnodata -a -v -9999
+
+# Convert a CSV grid to GeoTIFF
+csv2tif -i grid.csv -o dem.tif --xll 250000 --yll 2500000 --cellsize 5 --epsg 32648
+```
+
+### Python API
+
+```python
+from demtools import chgnodata, setnodata, csv2tif
+
+# Convert nodata value in a raster file
+chgnodata.convert_nodata_value("dem.tif", output_nodata=-9999)
+
+# Assign nodata metadata without touching pixel values
+setnodata.define_nodata_value("dem.tif", nodata_value=-9999)
+
+# Convert CSV grid to GeoTIFF
+csv2tif.csv_to_tif(
+    input_csv="grid.csv",
+    output_tif="dem.tif",
+    xll=250000,
+    yll=2500000,
+    cellsize=5,
+    nodata=-999,
+    epsg=32648,
+)
+```
+
+## Notes
+
+- `chgnodata` and `setnodata` automatically create a `RAS_BAK/` directory with backup copies before modifying files.
+- `csv2tif` expects a plain numeric CSV (no headers), with rows ordered from north to south.
+- GDAL must be installed separately via conda or a pre-built wheel; it is not listed in `requirements.txt` as it cannot be reliably installed via pip on all platforms.

demtools-0.1.0/demtools.egg-info/SOURCES.txt

@@ -0,0 +1,17 @@
+MANIFEST.in
+README.md
+pyproject.toml
+requirements.txt
+setup.py
+demtools/__init__.py
+demtools/chgnodata.py
+demtools/csv2tif.py
+demtools/describe.py
+demtools/setnodata.py
+demtools.egg-info/PKG-INFO
+demtools.egg-info/SOURCES.txt
+demtools.egg-info/dependency_links.txt
+demtools.egg-info/entry_points.txt
+demtools.egg-info/requires.txt
+demtools.egg-info/top_level.txt
+demtools/tests/__init__.py

demtools-0.1.0/demtools.egg-info/dependency_links.txt

@@ -0,0 +1 @@
+

demtools-0.1.0/demtools.egg-info/entry_points.txt

@@ -0,0 +1,5 @@
+[console_scripts]
+chgnodata = demtools.chgnodata:main
+csv2tif = demtools.csv2tif:main
+demtools-info = demtools.describe:main
+setnodata = demtools.setnodata:main

demtools-0.1.0/demtools.egg-info/requires.txt

@@ -0,0 +1 @@
+numpy>=1.20.0

demtools-0.1.0/demtools.egg-info/top_level.txt

@@ -0,0 +1 @@
+demtools

demtools-0.1.0/pyproject.toml

@@ -0,0 +1,31 @@
+[build-system]
+requires = ["setuptools>=61.0", "wheel"]
+build-backend = "setuptools.build_meta"
+
+[project]
+name = "demtools"
+version = "0.1.0"
+authors = [
+    { name = "aaronchh", email = "aaronhsu219@gmail.com" },
+]
+description = "A collection of tools for working with DEM (Digital Elevation Model) raster files"
+readme = "README.md"
+license = { text = "MIT" }
+requires-python = ">=3.8"
+dependencies = [
+    "numpy>=1.20.0",
+]
+classifiers = [
+    "Programming Language :: Python :: 3",
+    "License :: OSI Approved :: MIT License",
+    "Operating System :: OS Independent",
+]
+
+[project.urls]
+Homepage = "https://github.com/AaronOET/demtools"
+
+[project.scripts]
+demtools-info = "demtools.describe:main"
+chgnodata = "demtools.chgnodata:main"
+setnodata = "demtools.setnodata:main"
+csv2tif = "demtools.csv2tif:main"

demtools-0.1.0/requirements.txt

@@ -0,0 +1 @@
+numpy>=1.20.0

demtools-0.1.0/setup.cfg

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

demtools-0.1.0/setup.py

@@ -0,0 +1,34 @@
+from setuptools import setup, find_packages
+
+with open("README.md", "r", encoding="utf-8") as fh:
+    long_description = fh.read()
+
+with open("requirements.txt", "r", encoding="utf-8") as f:
+    requirements = f.read().splitlines()
+
+setup(
+    name="demtools",
+    version="0.1.0",
+    author="aaronchh",
+    author_email="aaronhsu219@gmail.com",
+    description="A collection of tools for working with DEM (Digital Elevation Model) raster files",
+    long_description=long_description,
+    long_description_content_type="text/markdown",
+    url="https://github.com/AaronOET/demtools",
+    packages=find_packages(),
+    classifiers=[
+        "Programming Language :: Python :: 3",
+        "License :: OSI Approved :: MIT License",
+        "Operating System :: OS Independent",
+    ],
+    python_requires=">=3.8",
+    install_requires=requirements,
+    entry_points={
+        "console_scripts": [
+            "demtools-info=demtools.describe:main",
+            "chgnodata=demtools.chgnodata:main",
+            "setnodata=demtools.setnodata:main",
+            "csv2tif=demtools.csv2tif:main",
+        ],
+    },
+)