There are some functions in the 'raster' package. Perhaps something like: library(raster) x <- raster("vegetation file") extract(x, houses, buffer=100) -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/Buffering-in-RS-Image-tp5887532p5890394.html Sent from the R-sig-geo mailing list archive at Nabble.com.