Raster clipping issue

Question

I need to clip a raster using a shapefile (many polygons), but I want it to clip only the pixels that are completely inside the polygons. I mean, I do not want to include the pixels that are part outside the boundaries of the polygons, is there any way to do this?

I know I can convert the raster into a shp and clip it but I do not want the pixels to split either; I need only whole pixels covered by the polygons.

Ideas anyone?

Eva

curtvprice · Answer

If the raster cells are much bigger than your polygons, a straight-forward approach would be to convert your clipped raster cells (some of which are totally inside, some partially outside) to polygon using the Raster To Polygon tool, and use the Select By Location tool to select the polygons that are entirely inside (CONTAINED BY) your clip polygons. You could then leave your cells as polygons or convert back to raster using Polygon To Raster (using the extent and snap raster and cell size of your clipped raster so it will line up).

Dan's other suggestion is also a good approach -- take your clip polygons, buffer them with a negative value of half the diagonal of a cell (cellsize * 1.41 * 0.5), and clip by that polygon. If the cells are not as huge as I envisioned from your description, this would be a better approach, as zillions of cell polygons are kind of hard to work with.

EvaPintado · Answer

Thank you Dan. It is not a question of liking or not the jagged look, I need to have only the pixels completely inside the polygons (which are many and very irregular in shape) to do some calculations. Thank you anyway, will keep looking to see if can  filter them this in any way.

DanPatterson_Retired · Answer

You can't change the rule, if the centre is in the polygon, then you have to live with the fact that it will 'appear' that part of a 'cell' is split by the boundary.  Consider a perfectly horizontal or vertical line.... that feature will connect the centers of the cells.  Shifting will only move the problem in the other direction.  Or consider a polygon square.  In order for a raster cell to be completely within the square, then all 4 corners would have to be within the bounds of the polygon.  The only way this can happen is if the raster representation of the square polygon is shrunk by 1 cell.  If you don't like the jagged look, then you can use a smaller cell size along with limiting the zoom extent of a layer.  It is a difference between raster and vector data that is not surmountable even if using hexagonal raster cells or cells where the X and Y sizes can vary.