Here is my best shot:
library(raster) ## An example raster logo <- raster(system.file("external/rlogo.grd", package="raster")) ## Clip out the lower-left 50*100 pixel rectangle as a new raster 'r' cropWithRowCol <- function(r, rows, cols) { cc <- cellFromRowColCombine(r, rownr=rows, colnr=cols) crop(r, rasterFromCells(r, cc, values=FALSE)) } r <- cropWithRowCol(logo, nrow(logo) - 49:0, 1:100) ## Extract multipliers used to appropriately size the selected device w <- ncol(r)/max(dim(r)) h <- nrow(r)/max(dim(r)) ## Set up appropriately sized device with no borders and required coordinate system ## png("eg.png", width=480*w, height=480*h) dev.new(width=5*w, height=5*h) plot.new() par(mar=c(0,0,0,0), oma=c(0,0,0,0)) plot.window(xlim=extent(r)[1:2], ylim=extent(r)[3:4], xaxs="i",yaxs="i") ## Finally, plot the (sub-)raster to it plot(r, axes=FALSE, legend=FALSE, add=TRUE) ## dev.off()
(Please remember that in an interactively resizable device, a resizing device will spoil the constructed map format.)

source share