r.in.nasadem

Creates a DEM from 1 arcsec NASADEM tiles.

r.in.nasadem [-z] [version=string] [layer=string] output=name [username=string] [password=string] [url=string] [local=name] [region=xmin,ymin,xmax,ymax] [method=string] [resolution=float] [memory=integer] [--overwrite] [--verbose] [--quiet] [--qq] [--ui]

Example:

r.in.nasadem output=name

grass.script.run_command("r.in.nasadem", version="NASADEM_HGT.001", layer="hgt", output, username=None, password=None, url="https://e4ftl01.cr.usgs.gov/MEASURES", local=None, region=None, method="bilinear_f", resolution=None, memory=300, flags=None, overwrite=False, verbose=False, quiet=False, superquiet=False)

Example:

gs.run_command("r.in.nasadem", output="name")

Parameters

Command linePython (grass.script)

version=string
    NASADEM version
    Allowed values: NASADEM_HGT.001
    Default: NASADEM_HGT.001
layer=string
    NASADEM layer to import
    Allowed values: hgt, num, swb
    Default: hgt
    hgt: Void-filled DEM
    num: Number of scenes
    swb: Updated SRTM water body data
output=name [required]
    Name for output raster map
username=string
    Username for authentication
password=string
    Password for authentication
url=string
    Base URL to fetch NASADEM tiles
    Default: https://e4ftl01.cr.usgs.gov/MEASURES
local=name
    Local folder with NASADEM tiles
    Use local folder instead of URL to retrieve NASADEM tiles
region=xmin,ymin,xmax,ymax
    Import subregion only (default is current region)
    Format: xmin,ymin,xmax,ymax - usually W,S,E,N
method=string
    Resampling method to use for reprojection (required if location projection not longlat)
    Allowed values: nearest, bilinear, bicubic, lanczos, bilinear_f, bicubic_f, lanczos_f
    Default: bilinear_f
    nearest: nearest neighbor
    bilinear: bilinear interpolation
    bicubic: bicubic interpolation
    lanczos: lanczos filter
    bilinear_f: bilinear interpolation with fallback
    bicubic_f: bicubic interpolation with fallback
    lanczos_f: lanczos filter with fallback
resolution=float
    Resolution of output raster map (required if location projection not longlat)
memory=integer
    Maximum memory to be used (in MB)
    Cache size for raster rows
    Default: 300
-z
    Create zero elevation for missing tiles
--overwrite
    Allow output files to overwrite existing files
--help
    Print usage summary
--verbose
    Verbose module output
--quiet
    Quiet module output
--qq
    Very quiet module output
--ui
    Force launching GUI dialog

version : str, optional
    NASADEM version
    Allowed values: NASADEM_HGT.001
    Default: NASADEM_HGT.001
layer : str, optional
    NASADEM layer to import
    Allowed values: hgt, num, swb
    hgt: Void-filled DEM
    num: Number of scenes
    swb: Updated SRTM water body data
    Default: hgt
output : str, required
    Name for output raster map
    Used as: output, raster, name
username : str, optional
    Username for authentication
password : str, optional
    Password for authentication
url : str, optional
    Base URL to fetch NASADEM tiles
    Default: https://e4ftl01.cr.usgs.gov/MEASURES
local : str, optional
    Local folder with NASADEM tiles
    Use local folder instead of URL to retrieve NASADEM tiles
    Used as: input, dir, name
region : tuple[float, float, float, float] | list[float] | str, optional
    Import subregion only (default is current region)
    Format: xmin,ymin,xmax,ymax - usually W,S,E,N
    Used as: xmin,ymin,xmax,ymax
method : str, optional
    Resampling method to use for reprojection (required if location projection not longlat)
    Allowed values: nearest, bilinear, bicubic, lanczos, bilinear_f, bicubic_f, lanczos_f
    nearest: nearest neighbor
    bilinear: bilinear interpolation
    bicubic: bicubic interpolation
    lanczos: lanczos filter
    bilinear_f: bilinear interpolation with fallback
    bicubic_f: bicubic interpolation with fallback
    lanczos_f: lanczos filter with fallback
    Default: bilinear_f
resolution : float, optional
    Resolution of output raster map (required if location projection not longlat)
memory : int, optional
    Maximum memory to be used (in MB)
    Cache size for raster rows
    Default: 300
flags : str, optional
    Allowed values: z
    z
        Create zero elevation for missing tiles
overwrite: bool, optional
    Allow output files to overwrite existing files
    Default: False
verbose: bool, optional
    Verbose module output
    Default: False
quiet: bool, optional
    Quiet module output
    Default: False
superquiet: bool, optional
    Very quiet module output
    Default: False

DESCRIPTION

r.in.nasadem imports all NASADEM tiles covering the current region or region extents given with the region option into GRASS and patches the tiles together.

Cited from https://lpdaac.usgs.gov/news/release-nasadem-data-products/:

"NASADEM extends the legacy of the Shuttle Radar Topography Mission (SRTM) by improving the digital elevation model (DEM) height accuracy and data coverage as well as providing additional SRTM radar-related data products. The improvements were achieved by reprocessing the original SRTM radar signal data and telemetry data with updated algorithms and auxiliary data not available at the time of the original SRTM processing."

In total, 14,520 NASADEM tiles are available.

r.in.nasadem imports data from NASADEM_HGT.001, or from a local copy of these data provided with the local option.

There are three different layers available in the NASADEM_HGT.001 product: hgt, num, and swb, each of these layers can be imported with r.in.nasadem with the layer option. Importantly, for the NASADEM tiles download a user registration is needed at https://urs.earthdata.nasa.gov/users/new

In the earthdata user profile, two specific applications must be approved in the "My application" tab:

"LP DAAC Data Pool" application, and
"Earthdata Search" application.

EXAMPLE

Import of NASADEM_HGT.001 covering the current computational region:

# run in LatLong location - Sicily East, Italy
g.region n=39 s=37 w=14 e=16 res=0:00:01 -p

# use own credentials here
r.in.nasadem user="my_nasa_user" password="my_nasa_pw" output=nasadem_sicily_1arc memory=2000
r.univar nasadem_sicily_1arc

Figure: Eta volcano (Sicily, Italy) shown in NVIZ

REFERENCES

M. Neteler, 2005. SRTM and VMAP0 data in OGR and GRASS. GRASS Newsletter, Vol.3, pp. 2-6, June 2005. ISSN 1614-8746.

AUTHORS

Markus Metz
Reprojection support: Anika Bettge, mundialis

SOURCE CODE

Available at: r.in.nasadem source code (history)
Latest change: Friday Feb 21 10:10:05 2025 in commit 7d78fe3