Read Surfer Grid File

Read an Surfer ASCII grid file

import pooch

from PVGeo.grids import SurferGridReader

Download a sample Surfer grid file

fname = "surfer-grid.grd"
url = "https://github.com/OpenGeoVis/PVGeo/raw/main/tests/data/{}".format(fname)
file_path = pooch.retrieve(url=url, known_hash=None)

dem = SurferGridReader().apply(file_path)
dem

Header	Data Arrays
ImageData Information N Cells 40001 N Points 40404 X Bounds 0.000e+00, 1.105e+06 Y Bounds 0.000e+00, 9.050e+05 Z Bounds 0.000e+00, 0.000e+00 Dimensions 222, 182, 1 Spacing 5.000e+03, 5.000e+03, 1.000e+00 N Arrays 1	Name Field Type N Comp Min Max Data Points float64 1 -2.519e+02 1.754e+02

Apply a filter to the DEM to have realistic topography

warped = dem.warp_by_scalar(scale_factor=300.0)
warped.plot(cmap="terrain")

Static Scene

Interactive Scene