Cut hole in raster with shapefile and sample at closest edges

Geographic Information Systems Asked on January 3, 2022

I have a raster in which I’d like to cut a “hole” with a shapefile. I have points close to the boundary, but inside the hole, and I would like to find the closest point on the edge. To give a concrete example, I have a global dataset (1′ resolution) and a shapefile of the coast of Africa. My data points are locations of cities near the coast, but some of these points are not exactly on the coast, so I’d like to find the closest point outside the shapefile, i.e. closest point over water.

Normally, I’d use: gdalwarp -of -GTiff -cutline Coast.shp -crop_to_cutline -dstnodata 0.0 Data.tif Clipped_Data.tif but this obviously only gives me what’s inside the shapefile; I can’t find an inverse clipping for this command.

I did some digging and found this, which allowed me to set everything inside the shapefile to 0, but when I sample using gmt grdtrack it just sets the values to 0 (which makes sense). Also, when I tried crop=True, invert=True I got: ValueError: crop and invert cannot both be True. So I reverted to using crop=False, invert=True.

So, how can I sample this raster at the closest point outside the shapefile? I think if I can get the hole to be set to NaN, I can use the -T flag in gmt grdtrack.

One Answer

Take a look at mapproject -L option. I think it does what you want and no need to use a grid.

Answered by J Luis on January 3, 2022

Add your own answers!

Related Questions

Zonal statistics (mode, variety) in PostGIS

0  Asked on September 1, 2020 by albine-pro


Converting GRIB2 to GeoTIFF

1  Asked on September 1, 2020 by arnpry


Trouble installing the ODC-Core in docker on Ubuntu 18.04?

0  Asked on September 1, 2020 by sargam-panwar


Lizmap URL highligt and zoom to feature on layer

0  Asked on August 31, 2020 by geom


Split and export quadrant of a big map

1  Asked on August 30, 2020 by pileoni


Bulk Snapping Points to Line Vertex in QGIS?

2  Asked on August 30, 2020 by npt_rich


Kilometer to degree and back

1  Asked on August 29, 2020 by salrad


Writing a 3-band image with rasterio

2  Asked on August 28, 2020 by sam-murphy


Batch export rasters as rendered images

1  Asked on August 28, 2020


Ask a Question

Get help from others!

© 2023 All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP