首先插入R包raster
install.packages("raster")
library(raster)
raster为旧版本,使用新的 terra 版本体验可能更好
这里由于我使用的空间操作不多,所以这里运用raster。
链接: raster包系列帮助
链接:terra版本帮助
创建简单栅格数据
r <- raster(ncol=10, nrow=10, xmx=-80, xmn=-150, ymn=20, ymx=60)
#创建名字为r的栅格,行列号为10*10,X、Y轴的范围为80-150,20-60
r
values(r)<-1:ncell(r)
plot(r)
结果
class : RasterLayer
dimensions : 10, 10, 100 (nrow, ncol, ncell)
resolution : 7, 4 (x, y)
extent : -150, -80, 20, 60 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +no_defs
读取单张影像
filename <- system.file("external/lux.shp", package="raster")
filename
## [1] "C:/soft/R/R-devel/library/raster/external/lux.shp"
这是帮助文件的代码,不知道为啥我读取不了,所以用了下面这个:
f<-r"(external/lux.shp)"
r1<-raster(f)
读取成功!
再来读取一连串,非常好用!
tif_file_name <- r"(E:\02_Project\01_Chlorophyll\ClimateZone\Split\A_LCC0.TIF)"
tif_file <- raster(tif_file_name)
引用来自https://blog.csdn.net/zhebushibiaoshifu/article/details/128485386
完成!