raster

My input raster consists of multiple layers, each with an image area surrounded by no data values. These layers don't overlap completely, and I am trying to output a file which consists of only the intersection of all bands (the zone which has no NoData values on any layer). The following works for a few layers, but not for the 50+ that I have in my real files (at least 3000x3000 pixels): library(raster) fin = "D:\\temp\\all_modes.pix" fout = "D:\\temp\\test.pix" inbands = stack(fin, bands = c(3:20)) NAvalue(inbands) = 0 # Not great: #out = all(is.na(inbands) == FALSE) * inbands #writeRaster

I need to use the function extract() to do a weighted average extraction from a raster using a grid cell of equal sized squares. My polygon grid is in UTM21n and the raster is in GCS WGS84 datum D. Do I have to reproject the raster before using it into extract()? Or will the function handle it properly? You can find the source code of function extract for SpatialPolygons here . The code starts with the following snippet: setMethod('extract', signature(x='Raster', y='SpatialPolygons'), function(x, y, fun=NULL, na.rm=FALSE, weights=FALSE, cellnumbers=FALSE, small=FALSE, df=FALSE, layer, nl,

I have a 2D footprint result for an eddy-covariance system on a 1100*1100 square meter, with the domain of domain = c(-100,1000,-100, 1000) the cell size of the domain (each grid) is 2 meters with the origin (0,0) accordingly on this location: Lon = -97.191391 #longtitude Lat = 36.055935 #latitude An example data (9.01 Mb) is attached here. FFP.rds then I can plot a Bing map on my site as this: library(OpenStreetMap) library(rgdal) map <- openmap(c(36.05778,-97.19250), c(36.05444,-97.18861),type='bing') plot(map) and I use image.plot() from package "fields" to plot the footprint with these

I need to get the cell numbers within a distance (e.g., 10 km) of one point using R, but I didn't figure out how to handle it for the raster data. library(raster) r <- raster(ncols=10, nrows=10) cellFromXY(r, c(2,2)) # get the cell number of one point How to get the cell numbers within a distance of one point? Use extract with the buffer option and cellnumbers=TRUE r <- raster(ncols=100, nrows=100) r[]<-runif(ncell(r)) xy <- cbind(-50, seq(-80, 80, by=20)) extract(r, xy[1:3,], buffer=1000000,cellnumbers=TRUE) 来源： https://stackoverflow.com/questions/13996256/get-cell-number-within-a-distance-of

I am trying to plot some global maps (raster files) and I have some problems in setting up a good color scale for my data. What I would like to do is to plot my data using a divergent palette (e.g. cm.colors ), and I would like to center the color "white" of such scale with the value zero, but without having to set symmetric values in the scale (i.e. the same value both negative and positive, i.e. limits=c(-1,1) ). Additionally, I would like to plot all values above and/or below a certain value all with the same color. In other words, if we suppose that my map has a range of -100 to 150, I