How to replace NA's in a raster object

r spatial raster
ils picture ils · Aug 15, 2012 · Viewed 16.7k times · Source

I need to replace the NA's in the raster object (r) from the example below.

library(raster)
filename <- system.file("external/test.grd", package="raster")
r <- raster(filename)

I also tried to remove these these (and place the result in a data.frame), but to no avail.

dfr <- as.data.frame(r, na.rm=T)
summary(dfr)
# test       
# Min.   : 128.4  
# 1st Qu.: 293.2  
# Median : 371.4  
# Mean   : 423.2  
# 3rd Qu.: 499.8  
# Max.   :1805.8  
# NA's   :6097

Answer

Robert Hijmans picture Robert Hijmans · Feb 16, 2013

A more memory safe approach (for large files) would be to use reclassify:

library(raster)
filename <- system.file("external/test.grd", package="raster")
r <- raster(filename)
rna <- reclassify(r, cbind(NA, 250))

Related questions

convert matrix to raster in R

I have a matrix data with spatial coordinates and one variable. The spatial resolution is 1000 meters. > str(dat1) > List of 3 > $ x: num [1:710] 302340 303340 304340 305340 306340 ... > $ y: num [1:1241] 5431470 5432470 5433470 5434470 5435470 ... > $ z: num [1:710, 1:1241] 225 225 225 225 225 ... I want to convert it into raster format. > …

Read More
How to remove the box frame in "plot.raster" in R package "raster"

I need to remove the box frame around the figure in R package "raster", but I cannot figure out which argument I should change. The example is as follows: library(raster) r <- raster(nrows=10, ncols=10) r <- …

Read More
Bring rasters to same extent by filling with NA - in R

I have several cropped rasters with different geometry/outlines. Specifically spatial yield maps from several years of the same field, but the extent varies - the measurements were not always overall the whole field, but in some years only part …

Read More