我想创建一个带有四个波段的图像子集。因此,我在R中使用了裁剪功能。
A<-raster("L8_stacked.tif")
subset<-extent(c(639451, 660104, 5469254, 5489566))
B<-crop(A,subset)
结果,我在.tif文件中得到了只有一个波段的栅格。我是否必须定义其他选项才能获得具有4个波段的子集图像?
正如其他人在注释中已经指出的那样,该raster()
函数返回一个(单个)RasterLayer对象。如果要使用多层栅格对象,则需要使用stack()
或brick()
函数将图像加载到RIe中:
A <- stack("L8_stacked.tif")
如果随后将crop()
扩展区与函数配合使用,则结果应为栅格堆栈,其中包含原始图像中的所有波段。
要了解有关光栅包的更多信息,请阅读本文档。
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句