R語言,raster包静浴,alignextent對齊函數堰氓,animate動畫函數,aproxNA缺失填補函數马绝。
alignExtent兩個對象對齊
# Tue Sep 08 11:19:11 2020 -
# 字符編碼:UTF-8
# R 版本:R x64 4.0.2
# cgh163email@163.com
#個人筆記不負責任
# —— 拎了個梨??
rm(list=ls());gc()
.rs.restartR()
require(raster)
# Tue Sep 08 11:26:41 2020 -alignExtent兩個對象對齊-----------------------------
# 將范圍對象與另一個對象對齊可確保它具有相同的原點和分辨率豆赏。 由于數據不精確,只能用于調整對象富稻。 alignExtent不應用于強制匹配確實不匹配的數據
r <- raster()
e <- extent(-10.1, 9.9, -20.1, 19.9)
e
ea <- alignExtent(e, r) # e,r兩個對象對齊
ea
r
extent(r)
animate對Raster *對象的圖層進行動畫處理
# Tue Sep 08 13:47:47 2020 --animate對Raster *對象的圖層進行動畫處理----------------------------
b <- brick(system.file("external/rlogo.grd", package="raster")) # 創(chuàng)建一個RasterBrick對象
plot(b)
animate(b, n=2) #基礎
# 多參數:
animate(b,
n = 2, # 循環(huán)2次
pause = 1 # 層之間的停頓應間隔多長時間掷邦?
,main = 'aac' # 每個圖層的標題
,zlim = 2 # 要繪制的值的范圍
,maxpixels = 2 # 整數>0。用于圖的最大像元數
)
# Tue Sep 08 13:57:29 2020 - diy
fp.grd <- brick('惠東水底山環(huán)線徒步.grd')
plot(fp.grd)
animate(fp.grd,n = 5,main = 'aaac')
aproxNA通過在圖層之間進行插值來估算值為NA的像元值
# Tue Sep 08 14:02:12 2020 --aproxNA通過在圖層之間進行插值來估算值為NA的像元值----------------------------
r <- raster(ncols=5, nrows=5)
r1 <- setValues(r, runif(ncell(r)))
plot(r1)
r2 <- setValues(r, runif(ncell(r)))
plot(r2)
r3 <- setValues(r, runif(ncell(r)))
plot(r3)
r4 <- setValues(r, runif(ncell(r)))
plot(r4)
r5 <- setValues(r, NA)
plot(r5)
r6 <- setValues(r, runif(ncell(r)))
plot(r6)
r1[6:10] <- NA
r2[5:15] <- NA
r3[8:25] <- NA
s <- stack(r1,r2,r3,r4,r5,r6)
plot(s)
s[1:5] <- NA
plot(s)
x1 <- approxNA(s)
plot(x1)
x2 <- approxNA(s, rule=2)
plot(x2)
x3 <- approxNA(s, rule=2, z=c(1,2,3,5,14,15))
plot(x3)
animate(s,n=2)
animate(x3,n=2)
# Tue Sep 08 14:12:20 2020 -End
缺失填補
缺失填補