我试图为两个保存为单带的多波段images.Output运行马赛克。

Question

这是我得到的最后输出，我应该把最后的输出作为一个包含两个波段的文件：

​

​

以下是我使用的代码：

A11 <-brick("E:/Official/PROJECTS/R_Progrm/1.tif") // to read multiband image
B11<-brick("E:/Official/PROJECTS/R_Progrm/3.tif")  // To read multiband image
mos1 <- mosaic(A11,B11,fun=max,tolerance=0.5, 
filename="Mosaic_new",overwrite=TRUE)
plot(mos1,main="Mosaic_new1")
writeRaster(x=mos1,file="E:/Official/PROJECTS/R_Progrm/M11.tif",options="INTERLEAVE=BAND",format="GTiff",datatype="FLT8S",overwrite=TRUE)

Answer 1

在你的问题中所显示的情节，是显示你的输出图像的两个波段。因此，您的代码及其输出不应该有任何问题。如果问题与将所有波段可视化为RGB图像有关，则必须修改绘图函数的参数，这意味着必须提供波段组合。例如：

plotRGB(a, r = 4, g = 3, b = 2, axes=TRUE, main="3 Band Color Composite Image")
box(col="white")

另外，您也可以尝试下面给出的代码，这些代码对我来说很好，我希望它能解决您的问题。

a <- stack("Path to first raster")
b <- stack("Path to second raster")
rast.list <- list(a,b)     
rast.list$fun <- mean
rast.mosaic <- do.call(mosaic,rast.list)
plot(rast.mosaic)
writeRaster(rast.mosaic,"Output_Raster_Name",format="GTiff",overwrite=TRUE)

Answer 2

rm(list = ls())
gc() 
memory.limit(size= 2000)
library(rgdal)
library(raster)
install.packages("gdalUtils")
library(gdalUtils)
library(sp)
setwd("E:/Official/PROJECTS/R_Progrm/MOs/")
list.files()
file1=file.path(getwd(), "", "1.tif")
gdal_setInstallation()
valid_install <- !is.null(getOption("gdalUtils_gdalPath"))
if(require(raster) && require(rgdal) && valid_install)
{
 layer1 <- file.path(getwd(), "", "1.tif")
 layer2 <- file.path(getwd(), "", "3.tif")
 file_list=c(layer1,layer2)
 mosaic_rasters(gdalfile=file_list,dst_dataset="E:/Official/PROJECTS/R_Progrm/MOs//test_mosaic.GTiff",separate=TRUE,of="GTiff",verbose=TRUE)
gdalinfo("test_mosaic.GTiff")
}