具有多个波段的R剪辑光栅

具有多个波段的R剪辑光栅,r,crop,subset,raster,r-raster,R,Crop,Subset,Raster,R Raster,我想创建一个具有四个波段的图像子集。因此,我在R中使用crop函数 A<-raster("L8_stacked.tif") subset<-extent(c(639451, 660104, 5469254, 5489566)) B<-crop(A,subset) A正如其他注释中已经指出的那样,raster()函数返回一个(单个)RasterLayer对象。如果需要多层光栅对象,则需要使用stack()或brick()函数将图像加载到R中。例如: A <- stack(

我想创建一个具有四个波段的图像子集。因此,我在R中使用crop函数

A<-raster("L8_stacked.tif")
subset<-extent(c(639451, 660104, 5469254, 5489566))
B<-crop(A,subset)

A正如其他注释中已经指出的那样,
raster()
函数返回一个(单个)RasterLayer对象。如果需要多层光栅对象,则需要使用
stack()
brick()
函数将图像加载到R中。例如:

A <- stack("L8_stacked.tif")

A你在你的工作区里用
B
能得到什么?(就此而言,您是否验证了
A
包含您期望的四个频带/层?)A本身只有一个频带。试试nlayers(A)-对吗?从砖(“L8_stacked.tif”)开始,然后裁剪。