java gdal 拼接多波段影像

Java可以使用GDAL的Java绑定库来拼接多波段影像。以下是一个简单的示例代码：

import org.gdal.gdal.Dataset;
import org.gdal.gdal.gdal;
import org.gdal.gdalconst.gdalconstConstants;
import org.gdal.osr.SpatialReference;

public class ImageMosaic {

    public static void main(String[] args) {
        // 注册GDAL驱动
        gdal.AllRegister();

        // 打开第一幅影像
        Dataset src1 = gdal.Open("image1.tif", gdalconstConstants.GA_ReadOnly);
        double[] geoTransform1 = src1.GetGeoTransform();
        SpatialReference srs1 = new SpatialReference(src1.GetProjection());

        // 打开第二幅影像
        Dataset src2 = gdal.Open("image2.tif", gdalconstConstants.GA_ReadOnly);
        double[] geoTransform2 = src2.GetGeoTransform();
        SpatialReference srs2 = new SpatialReference(src2.GetProjection());

        // 计算新影像的范围和分辨率
        double minX = Math.min(geoTransform1[0], geoTransform2[0]);
        double maxX = Math.max(geoTransform1[0] + geoTransform1[1] * src1.getRasterXSize(),
                geoTransform2[0] + geoTransform2[1] * src2.getRasterXSize());
        double minY = Math.min(geoTransform1[3] + geoTransform1[5] * src1.getRasterYSize(),
                geoTransform2[3] + geoTransform2[5] * src2.getRasterYSize());
        double maxY = Math.max(geoTransform1[3], geoTransform2[3]);
        double pixelSize = Math.min(geoTransform1[1], geoTransform2[1]);

        // 创建新影像
        Dataset dst = gdal.GetDriverByName("GTiff").Create("mosaic.tif", 
                (int) ((maxX - minX) / pixelSize), (int) ((maxY - minY) / pixelSize), 
                src1.getRasterCount() + src2.getRasterCount(), gdalconstConstants.GDT_Float32);
        dst.SetGeoTransform(new double[] { minX, pixelSize, 0, maxY, 0, -pixelSize });
        dst.SetProjection(srs1.ExportToWkt());

        // 拼接第一幅影像
        for (int i = 1; i <= src1.getRasterCount(); i++) {
            dst.AddBand(gdalconstConstants.GDT_Float32);
            dst.GetRasterBand(i).WriteRaster(0, 0, src1.getRasterXSize(), src1.getRasterYSize(),
                    src1.GetRasterBand(i).ReadRaster(0, 0, src1.getRasterXSize(), src1.getRasterYSize(), null),
                    src1.getRasterXSize(), src1.getRasterYSize(), gdalconstConstants.GDT_Float32, 0, 0);
        }

        // 拼接第二幅影像
        for (int i = 1; i <= src2.getRasterCount(); i++) {
            dst.AddBand(gdalconstConstants.GDT_Float32);
            dst.GetRasterBand(i + src1.getRasterCount()).WriteRaster(0, 0, src2.getRasterXSize(), src2.getRasterYSize(),
                    src2.GetRasterBand(i).ReadRaster(0, 0, src2.getRasterXSize(), src2.getRasterYSize(), null),
                    src2.getRasterXSize(), src2.getRasterYSize(), gdalconstConstants.GDT_Float32, 0, 0);
        }

        // 关闭影像
        src1.delete();
        src2.delete();
        dst.delete();
    }

}

这个示例代码假设两幅影像的分辨率相同，如果分辨率不同，则需要对代码进行修改。