在R中镶嵌多波段光栅时的错误

8yoxcaq7 于 2023-05-20 发布在其他

关注(0)|答案(3)|浏览(138)

我试图将多个多波段数据集镶嵌到R中的一个大型多波段栅格中。我试图使用此代码进行此操作，但它返回一个单波段图像。

library(sp)  
library(raster)  
library(rgdal)  
setwd("C:\\Projects\\Rice-fallow_4_states\\Bihar\\S1")  
x <- list.files(".\\", pattern='tif$',
                    full.names = TRUE) # list all the rasters  
X1<- as.list(x)    

allrasters1 <- lapply(X1, raster)  
x.mosaic <- do.call(merge,allrasters1)   
names(x.mosaic)  
x.mosaic  
plot(x.mosaic)

来源：https://stackoverflow.com/questions/76259972/error-in-mosaicing-multi-band-raster-in-r

3条答案

按热度按时间

qltillow1#

我建议你应该使用terra代替raster包。terra比raster快得多。可以使用以下代码

library(terra) 

setwd("C:\\Projects\\Rice-fallow_4_states\\Bihar\\S1")  
x <- list.files(".\\", pattern='tif$',
                full.names = TRUE) # list all the rasters  

allrasters1 <- lapply(x, rast)  
x.mosaic <- do.call(mosaic, allrasters1)  
names(x.mosaic)  
x.mosaic  
plot(x.mosaic)

赞(0）回复(0）举报 2023-05-20

jjjwad0x2#

这就是如何使用Terra（“栅格”的替换）来实现这一点

library(terra)
ff <- list.files(".\\", pattern='tif$',
                    full.names = TRUE) 
x <- sprc(ff) 
m <- mosaic(x)

赞(0）回复(0）举报 2023-05-20

dfuffjeb3#

在R中使用raster()加载光栅时，它只加载图像的一个层。如果要加载整个多波段图像，请使用stack()或brick()。我相信brick()是最适合多波段图像的，但我不经常处理这些，所以我还没有找到stack()之间有意义的区别。

赞(0）回复(0）举报 2023-05-20

我来回答

在R中镶嵌多波段光栅时的错误

3条答案

相关问题

热门标签

最新问答