java gdal波段合成

你可以使用Java的GDAL库来实现波段合成。下面是一个简单的示例代码：

import org.gdal.gdal.*;
import org.gdal.gdalconst.*;

public class BandMerge {

  public static void main(String[] args) {
    // Register all GDAL drivers
    gdal.AllRegister();

    // Open the input files
    Dataset input1 = gdal.Open("input1.tif", gdalconst.GA_ReadOnly);
    Dataset input2 = gdal.Open("input2.tif", gdalconst.GA_ReadOnly);

    // Get the number of bands in each input file
    int bands1 = input1.getRasterCount();
    int bands2 = input2.getRasterCount();

    // Create the output file
    String[] options = new String[]{"-of", "GTiff"};
    Dataset output = gdal.GetDriverByName("GTiff").Create(
      "output.tif",
      input1.getRasterXSize(),
      input1.getRasterYSize(),
      bands1 + bands2,
      gdalconst.GDT_Float32,
      options
    );

    // Copy the geotransform and projection from the first input file
    output.SetGeoTransform(input1.GetGeoTransform());
    output.SetProjection(input1.GetProjection());

    // Copy the data from the input files
    for (int i = 1; i <= bands1; i++) {
      Band band = input1.GetRasterBand(i);
      output.WriteRaster(0, 0, input1.getRasterXSize(), input1.getRasterYSize(), band.ReadRaster(0, 0, input1.getRasterXSize(), input1.getRasterYSize()));
    }
    for (int i = 1; i <= bands2; i++) {
      Band band = input2.GetRasterBand(i);
      output.WriteRaster(0, 0, input2.getRasterXSize(), input2.getRasterYSize(), band.ReadRaster(0, 0, input2.getRasterXSize(), input2.getRasterYSize()));
    }

    // Close the datasets
    input1.delete();
    input2.delete();
    output.delete();
  }

}

这个代码示例会将两个输入文件中的所有波段合成为一个输出文件。你可以根据自己的需求进行修改和扩展。