假設爲超過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。 –