Read ESRI Grid File

Read an ESRI ASCII grid file

import os

import pooch

from PVGeo.grids import EsriGridReader

Download a sample ESRI grid file

url = "https://raw.githubusercontent.com/pyvista/vtk-data/master/Data/esri_grid.dem.zip"
file_paths = pooch.retrieve(url=url, known_hash=None, processor=pooch.Unzip())
file_path = [f for f in file_paths if os.path.basename(f) == "esri_grid.dem"][0]
file_path

'/home/runner/.cache/pooch/fb1e1863d6713886b69896153e1b360e-esri_grid.dem.zip.unzip/esri_grid.dem'

dem = EsriGridReader().apply(file_path)
dem

Header	Data Arrays
ImageData Information N Cells 4875263 N Points 4879680 X Bounds 3.666e+05, 3.732e+05 Y Bounds 3.733e+06, 3.740e+06 Z Bounds 0.000e+00, 0.000e+00 Dimensions 2208, 2210, 1 Spacing 3.000e+00, 3.000e+00, 3.000e+00 N Arrays 1	Name Field Type N Comp Min Max Data Points float64 1 -3.936e+02 3.988e+02

Apply a filter to the DEM to have realistic topography

warped = dem.warp_by_scalar()
warped.plot(cmap="terrain", clim=[-100, 400])

Static Scene

Interactive Scene