用R语言加载单幅或多幅dat栅格数据

首先插入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
在这里插入图片描述
完成!

  • 7
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值