def test_grdclip_no_outgrid(grid):
"""
Test the below and above parameters for grdclip with no set outgrid.
"""
temp_grid = grdclip(grid=grid, below=[-1500, -1800], above=[-200, 40])
assert temp_grid.dims == ("lat", "lon")
assert temp_grid.gmt.gtype == 1 # Geographic grid
assert temp_grid.gmt.registration == 1 # Pixel registration
npt.assert_allclose(temp_grid.min(), -1800)
npt.assert_allclose(temp_grid.max(), 40)
def test_grdclip_no_outgrid(grid, expected_grid):
"""
Test the below and above parameters for grdclip with no set outgrid.
"""
temp_grid = grdclip(
grid=grid, below=[550, -1000], above=[700, 1000], region=[-53, -49, -19, -16]
)
assert temp_grid.dims == ("lat", "lon")
assert temp_grid.gmt.gtype == 1 # Geographic grid
assert temp_grid.gmt.registration == 1 # Pixel registration
xr.testing.assert_allclose(a=temp_grid, b=expected_grid)
def test_grdclip_outgrid(grid):
"""
Test the below and above parameters for grdclip and creates a test outgrid.
"""
with GMTTempFile(suffix=".nc") as tmpfile:
result = grdclip(grid=grid,
outgrid=tmpfile.name,
below=[-1500, -1800],
above=[-200, 40])
assert result is None # return value is None
assert os.path.exists(path=tmpfile.name) # check that outgrid exists
result = (grdinfo(grid=tmpfile.name, force_scan=0,
per_column="n").strip().split())
assert int(result[4]) == -1800 # minimum value
assert int(result[5]) == 40 # maximum value
def test_grdclip_outgrid(grid, expected_grid):
"""
Test the below and above parameters for grdclip and creates a test outgrid.
"""
with GMTTempFile(suffix=".nc") as tmpfile:
result = grdclip(
grid=grid,
outgrid=tmpfile.name,
below=[550, -1000],
above=[700, 1000],
region=[-53, -49, -19, -16],
)
assert result is None # return value is None
assert os.path.exists(path=tmpfile.name) # check that outgrid exists
temp_grid = load_dataarray(tmpfile.name)
assert temp_grid.dims == ("lat", "lon")
assert temp_grid.gmt.gtype == 1 # Geographic grid
assert temp_grid.gmt.registration == 1 # Pixel registration
xr.testing.assert_allclose(a=temp_grid, b=expected_grid)
-2000 m via the ``below`` parameter. """ import pygmt fig = pygmt.Figure() # Define region of interest around Iceland region = [-28, -10, 62, 68] # Load sample grid (3 arc minute global relief) in target area grid = pygmt.datasets.load_earth_relief(resolution="03m", region=region) # Plot original grid fig.basemap(region=region, projection="M12c", frame=["f", '+t"original grid"']) fig.grdimage(grid=grid, cmap="oleron") # Shift plot origin of the second map by "width of the first map + 0.5 cm" # in x direction fig.shift_origin(xshift="w+0.5c") # Set all grid points < 0 m to a value of -2000 m. grid = pygmt.grdclip(grid, below=[0, -2000]) # Plot clipped grid fig.basemap(region=region, projection="M12c", frame=["f", '+t"clipped grid"']) fig.grdimage(grid=grid) fig.colorbar(frame=["x+lElevation", "y+lm"], position="JMR+o0.5c/0c+w8c") fig.show()