Due to issues with Ghostscript, I updated GMT this morning to dev 6.5.0, after switching, grdclip() no longer clipped values correctly.
grid = pygmt.datasets.load_earth_relief(resolution=‘03s’, region = [min_max_lon[0], min_max_lon[1], min_max_lat[0], min_max_lat[1]])
print(f’Before ** Grid min {grid.min()}, Grid Max {grid.max()}’)
grid = pygmt.grdclip(grid, above=[0.0,10000.0])
print(f’After ** Grid min {grid.min()}, Grid Max {grid.max()}’)
Before ** Grid min <xarray.DataArray ‘z’ ()>
array(-3475.2976, dtype=float32), Grid Max <xarray.DataArray ‘z’ ()>
array(1173., dtype=float32)
After ** Grid min <xarray.DataArray ‘z’ ()>
array(-3475.2976, dtype=float32), Grid Max <xarray.DataArray ‘z’ ()>
array(1173., dtype=float32)
If I revert to GMT 6.4, this works correctly; however, transparency does not work due to incompatibility with Ghostscript.