Come sostituire NA in un oggetto raster

Question

Devo sostituire gli NA nell'oggetto raster (r) dall'esempio seguente.

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

Ho anche provato a rimuovere questi questi (e inserire il risultato in un data.frame), ma senza alcun risultato.

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 1

io non sono sicuro che abbia senso per rimuovere NA valori da un oggetto raster, ma si può facilmente sostituire.

Ad esempio:

oldpar <- par(mfrow=c(1, 2)) 
plot(r) 
r[is.na(r)] <- 250 
plot(r) 
par(oldpar) 

---

Se davvero si vuole, è possibile estrarre i valori raster in un vettore e quindi rimuovere i valori NA. (. Anche se, dal momento che si perdono le informazioni spaziali, non riesco a vedere come questo può essere utile)

r <- raster(filename) 

r <- values(r) 
head(r) 
[1] NA NA NA NA NA NA 

head(na.omit(r)) 
[1] 633.686 712.545 654.162 604.442 857.256 755.506 

Answer 2

Un approccio sicuro più memoria (per file di grandi dimensioni) sarebbe quella di utilizzare riclassificare:

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