2
假设为超过4频段我手动创建的像素值，我想将它们存储在一个TIFF文件。如何将超过4个波段的像素阵列保存到java中的Tiff/GeoTiff文件中？
那些波段可以用于R，G，B，温度（温度值不在0到255的范围内，因此我使用int代替pexils的字节）等等，即任何可以从卫星采取
现在我想将这些像素保存到一个tiff文件。在java中有一个BufferedImage类，它有很多类型，如：TYPE_4BYTE_ABGR，TYPE_BYTE_GRAY等。然而，它们都不适用于多频段多于4个频段。有TYPE_CUSTOM，但是当指定它并试图将数据保存到Tiff文件时，它会给你一个异常，因为它不支持写操作（仅用于读操作，即它可以读取文件并将其类型设置为TYPE_CUSTOM，如果它不明白该类型，但不能将该文件写入未理解的类型）。
以下代码适用于3个乐队，即使不正确（它不会显示彩色图像，它看起来像是一个失真的图像，但缺少线条）但是对于4个以上的乐队，我该怎么做？
ImageOutputStream ios = ImageIO.createImageOutputStream(os);
Iterator<ImageWriter> writers = ImageIO.getImageWritersByFormatName("tiff");
ImageWriter writer = writers.next();
writer.setOutput(ios);
int index = 0;
int[] pixels = new int[width*height*numberOfBands];
for (int i = 0; i < width; i++) {
for (int j = 0; j < height; j++) {
for (int k = 0; k < numberOfBands; k++) {
pixels[index++] = //any values;
}
}
}
DataBuffer dataBuffer = new DataBufferInt(pixels, pixels.length);
// Create Raster
WritableRaster writableRaster = Raster.createBandedRaster
(dataBuffer, width, height,
width, // scanlineStride
new int[numberOfBands], // bankIndices,
new int[numberOfBands], // bandOffsets,
null); // location
// Create the image
BufferedImage bufferedImage = new BufferedImage
(width, height, BufferedImage.TYPE_BYTE_RGB);
bufferedImage.setData(writableRaster);
IIOImage iioImage = new IIOImage(bufferedImage, null, null);
ImageWriteParam param = writer.getDefaultWriteParam();
writer.write(null, iioImage, param);
我使用GeoTools的方式
编辑： 根据@iant我改变了代码，但它是只给空白的透明背景，甚至我保留了相同的band数，即3个乐队。 @iant你可以检查下面的代码。
package examples;
import java.awt.image.WritableRaster;
import java.io.File;
import java.io.IOException;
import javax.media.jai.RasterFactory;
import org.geotools.coverage.CoverageFactoryFinder;
import org.geotools.coverage.grid.GridCoordinates2D;
import org.geotools.coverage.grid.GridCoverage2D;
import org.geotools.coverage.grid.GridCoverageFactory;
import org.geotools.coverage.grid.GridEnvelope2D;
import org.geotools.coverage.grid.GridGeometry2D;
import org.geotools.coverage.grid.io.AbstractGridFormat;
import org.geotools.coverage.grid.io.OverviewPolicy;
import org.geotools.gce.geotiff.GeoTiffFormat;
import org.geotools.gce.geotiff.GeoTiffReader;
import org.opengis.coverage.grid.GridCoverageWriter;
import org.opengis.parameter.GeneralParameterValue;
import org.opengis.parameter.ParameterValue;
public class CreateTiffImageTest2 {
public static void main(String[] args) throws IOException {
File file = new File("/home/mosab/Desktop/input/tif.tif");
ParameterValue<OverviewPolicy> policy = AbstractGridFormat.OVERVIEW_POLICY.createValue();
policy.setValue(OverviewPolicy.IGNORE);
ParameterValue<String> gridsize = AbstractGridFormat.SUGGESTED_TILE_SIZE.createValue();
ParameterValue<Boolean> useJaiRead = AbstractGridFormat.USE_JAI_IMAGEREAD.createValue();
useJaiRead.setValue(true);
GeoTiffReader geoTiffReader = new GeoTiffReader(file);
GridCoverage2D cov = geoTiffReader.read(new GeneralParameterValue[] { policy, gridsize, useJaiRead });
GridGeometry2D geometry = cov.getGridGeometry();
GridEnvelope2D gridEnvelope = geometry.getGridRange2D();
int w = (int) gridEnvelope.getWidth();
int h = (int) gridEnvelope.getHeight();
WritableRaster writableRaster = RasterFactory.createBandedRaster(java.awt.image.DataBuffer.TYPE_DOUBLE, w, h, 3,
null);
double[] data = new double[3];
double[] dest = new double[3];
for (int i = 0; i < w; i++) {
for (int j = 0; j < h; j++) {
GridCoordinates2D coord = new GridCoordinates2D(i, j);
cov.evaluate(coord, dest);
data[0] = dest[0];
data[1] = dest[1];
data[2] = dest[2];
writableRaster.setPixel(i, j, data);
}
float perc = 100.0f * i/w;
if (i % 100 == 0) {
System.out.println("done " + perc);
}
}
// Wrap the raster as a Coverage
GridCoverageFactory factory = CoverageFactoryFinder.getGridCoverageFactory(null);
GridCoverage2D gc = factory.create("name", writableRaster, cov.getEnvelope());
File out = new File("/home/mosab/Desktop/input/tifgen.tif");
GeoTiffFormat format = new GeoTiffFormat();
GridCoverageWriter writer = format.getWriter(out);
try {
writer.write(gc, null);
writer.dispose();
} catch (IllegalArgumentException | IOException e) {
// TODO Auto-generated catch block
e.printStackTrace();
}
}
}
更新2：
import java.awt.image.BufferedImage;
import java.awt.image.WritableRaster;
import java.io.File;
import java.io.IOException;
import javax.imageio.ImageIO;
import javax.media.jai.RasterFactory;
import org.geotools.coverage.CoverageFactoryFinder;
import org.geotools.coverage.grid.GridCoverage2D;
import org.geotools.coverage.grid.GridCoverageFactory;
import org.geotools.coverage.grid.io.AbstractGridFormat;
import org.geotools.coverage.grid.io.GridFormatFinder;
import org.geotools.factory.Hints;
import org.geotools.gce.geotiff.GeoTiffFormat;
import org.geotools.geometry.jts.ReferencedEnvelope;
import org.geotools.referencing.CRS;
import org.opengis.coverage.grid.GridCoverageWriter;
import org.opengis.geometry.MismatchedDimensionException;
import org.opengis.referencing.FactoryException;
import org.opengis.referencing.NoSuchAuthorityCodeException;
import org.opengis.referencing.crs.CoordinateReferenceSystem;
public class Test2 {
public static void print(Object o) {
System.out.println(o);
}
public static void main(String[] args)
throws MismatchedDimensionException, NoSuchAuthorityCodeException, FactoryException, IOException {
File out = new File("/home/mosab/Desktop/input/1.tif");
BufferedImage img = ImageIO.read(out);
// ColorModel colorModel = img.getColorModel(
WritableRaster raster = img.getRaster();
int w = img.getWidth();
int h = img.getHeight();
print("width = " + w);
print("heigh = " + h);
int numBands = raster.getNumBands();
WritableRaster writableRaster = RasterFactory.createBandedRaster(java.awt.image.DataBuffer.TYPE_INT, w, h, 3,
null);
//as I said pixels are created manually but I used here pixels from an image to check the approach
int[] data = new int[3];
for (int i = 0; i < w; i++) {
for (int j = 0; j < h; j++) {
for (int k = 0; k < numBands; k++) {
data[k] = raster.getSample(i, j, k);
}
writableRaster.setPixel(i, j, data);
}
}
GridCoverageFactory factory = CoverageFactoryFinder.getGridCoverageFactory(null);
CoordinateReferenceSystem crs = CRS.decode("EPSG:27700");
int llx = 500000;
int lly = 105000;
ReferencedEnvelope referencedEnvelope = new ReferencedEnvelope(llx, llx + (w * 10), lly, lly + (h * 10), crs);
GridCoverage2D gc = factory.create("name", writableRaster, referencedEnvelope);
AbstractGridFormat format = GridFormatFinder.findFormat(out);
Hints hints = null;
if (format instanceof GeoTiffFormat) {
hints = new Hints(Hints.FORCE_LONGITUDE_FIRST_AXIS_ORDER, Boolean.TRUE);
}
File out1 = new File("/home/mosab/Desktop/input/tifgen.tif");
GridCoverageWriter writer = format.getWriter(out1);
try {
writer.write(gc, null);
writer.dispose();
} catch (IllegalArgumentException | IOException e) {
// TODO Auto-generated catch block
e.printStackTrace();
}
}
}
说明：我用INT类型和尺寸3的阵列，因为原始图像具有RGB频带。
如果你有这样的特殊要求，你可能会考虑拥有自己的文件格式。 – Piglet
我会努力让它在使用3个乐队之前尝试4 –
@Piglet可以使用TIff/Geotiff格式完成，因为它是为此而设计的。顺便提一下，我之前提到了4个以上的乐队，因为在4的情况下，一些软件会认为它是ARGB，即+ aplha。 –