R GIS – How to Crop Simple Features Object

cliprsf

Is there a function to crop sf map object, similar to maptools::pruneMap(lines, xlim= c(4, 10), ylim= c(10, 15)) used for SpatialPolygon or SpatialLine?

I am considering st_intersection() but there can be proper way.

Best Answer

Since today, there is a st_crop function in the github version of sf (devtools::install_github("r-spatial/sf"), probably on CRAN in the near future too).

Just issue:

st_crop(nc, c(xmin=-82, xmax=-80, ymin=35, ymax=36))

The vector must be named with xmin xmax ymin ymax (in whichever order).

You can also use any object that can be read by st_bbox as cropping limits, which is very handy.

Related Solutions

[GIS] R – Crop a pixel image produced by spatstat with a boundary shape

The best approach is to coerce your spatstat "image" object to an sp object and use the over function or clip in the raster package. Here is a function to perform the coercion to an sp SpatialGridDataFrame.

as.SpatialGridDataFrame.im <- function(from){
    offset=c(from$xrange[1] + 0.5 * from$xstep, from$yrange[1] + 0.5 * from$ystep)
      cellsize = c(diff(from$xrange)/from$dim[2], diff(from$yrange)/from$dim[1])
dim = c(from$dim[2], from$dim[1])
  gt = GridTopology(offset, cellsize, dim)
    m = t(from$v[nrow(from$v):1, ])
  data = data.frame(v = as.vector(m))
SpatialGridDataFrame(gt, data)
}

[GIS] Clipping raster brick object in R

In order to write a raster stack or brick, the Interleave band can be used in the option argument.

    writeRaster(ndvi_stack, filename="<filename>.tif", options="INTERLEAVE=BAND", overwrite=TRUE)

Best Answer

Related Solutions

[GIS] R – Crop a pixel image produced by spatstat with a boundary shape

[GIS] Clipping raster brick object in R

Related Question