Jp2ExeOpImage(BandInfo bandInfo, MultiLevelModel imageModel, int level) throws IOException {
super(
DataBuffer.TYPE_USHORT,
bandInfo.imageLayout.width,
bandInfo.imageLayout.height,
getTileDim(bandInfo, level),
null,
ResolutionLevel.create(imageModel, level));
final File resolvedFile = bandInfo.imageFile.getCanonicalFile();
if (!resolvedFile.exists()) {
throw new FileNotFoundException("File not found: " + bandInfo.imageFile);
}
if (resolvedFile.getParentFile() == null) {
throw new IOException("Can't determine package directory");
}
final File cacheDir =
new File(
new File(SystemUtils.getApplicationDataDir(), "jopenjpeg/cache"),
resolvedFile.getParentFile().getName());
cacheDir.mkdirs();
if (!cacheDir.exists() || !cacheDir.isDirectory() || !cacheDir.canWrite()) {
throw new IOException("Can't access package cache directory");
}
this.bandInfo = bandInfo;
this.imageFile = resolvedFile;
this.cacheDir = cacheDir;
this.openFiles = new HashMap<File, Jp2File>();
this.locks = new HashMap<File, Object>();
}
private static void initGeoTools() {
// Must store EPSG database in BEAM home, otherwise it will be deleted from default temp
// location (Unix!, Windows?)
File epsgDir = new File(SystemUtils.getApplicationDataDir(true), EPSG_DATABASE_DIR_NAME);
System.setProperty(HsqlEpsgDatabase.DIRECTORY_KEY, epsgDir.getAbsolutePath());
}
/** The product reader for Radarsat2 products. */
public class Radarsat2ProductReader extends AbstractProductReader {
protected Radarsat2ProductDirectory dataDir = null;
private static final String lutsigma = "lutSigma";
private static final String lutgamma = "lutGamma";
private static final String lutbeta = "lutBeta";
private boolean isAscending;
private boolean isAntennaPointingRight;
private static final boolean flipToSARGeometry =
System.getProperty(SystemUtils.getApplicationContextId() + ".flip.to.sar.geometry", "false")
.equals("true");
/**
* Constructs a new abstract product reader.
*
* @param readerPlugIn the reader plug-in which created this reader, can be <code>null</code> for
* internal reader implementations
*/
public Radarsat2ProductReader(final ProductReaderPlugIn readerPlugIn) {
super(readerPlugIn);
}
/**
* Closes the access to all currently opened resources such as file input streams and all
* resources of this children directly owned by this reader. Its primary use is to allow the
* garbage collector to perform a vanilla job.
*
* <p>
*
* <p>This method should be called only if it is for sure that this object instance will never be
* used again. The results of referencing an instance of this class after a call to <code>close()
* </code> are undefined.
*
* <p>
*
* <p>Overrides of this method should always call <code>super.close();</code> after disposing this
* instance.
*
* @throws java.io.IOException if an I/O error occurs
*/
@Override
public void close() throws IOException {
if (dataDir != null) {
dataDir.close();
dataDir = null;
}
super.close();
}
/**
* Provides an implementation of the <code>readProductNodes</code> interface method. Clients
* implementing this method can be sure that the input object and eventually the subset
* information has already been set.
*
* <p>
*
* <p>This method is called as a last step in the <code>readProductNodes(input, subsetInfo)</code>
* method.
*
* @throws java.io.IOException if an I/O error occurs
*/
@Override
protected Product readProductNodesImpl() throws IOException {
Product product;
try {
final File fileFromInput = ReaderUtils.getFileFromInput(getInput());
dataDir = createDirectory(fileFromInput);
dataDir.readProductDirectory();
product = dataDir.createProduct();
addCalibrationLUT(product, fileFromInput.getParentFile());
product.getGcpGroup();
product.setFileLocation(fileFromInput);
product.setProductReader(this);
product.setModified(false);
final MetadataElement absRoot = AbstractMetadata.getAbstractedMetadata(product);
isAscending = absRoot.getAttributeString(AbstractMetadata.PASS).equals("ASCENDING");
isAntennaPointingRight =
absRoot.getAttributeString(AbstractMetadata.antenna_pointing).equals("right");
} catch (Exception e) {
Debug.trace(e.toString());
final IOException ioException = new IOException(e.getMessage());
ioException.initCause(e);
throw ioException;
}
return product;
}
protected Radarsat2ProductDirectory createDirectory(final File fileFromInput) {
return new Radarsat2ProductDirectory(fileFromInput, fileFromInput.getParentFile());
}
/**
* Read the LUT for use in calibration
*
* @param product the target product
* @param folder the folder containing the input
* @throws IOException if can't read lut
*/
private static void addCalibrationLUT(final Product product, final File folder)
throws IOException {
final MetadataElement absRoot = AbstractMetadata.getAbstractedMetadata(product);
final boolean isAscending =
absRoot.getAttributeString(AbstractMetadata.PASS).equals("ASCENDING");
final boolean isAntennaPointingRight =
absRoot.getAttributeString(AbstractMetadata.antenna_pointing).equals("right");
final boolean flipLUT =
flipToSARGeometry
&& ((isAscending && !isAntennaPointingRight)
|| (!isAscending && isAntennaPointingRight));
final File sigmaLUT = new File(folder, lutsigma + ".xml");
final File gammaLUT = new File(folder, lutgamma + ".xml");
final File betaLUT = new File(folder, lutbeta + ".xml");
final MetadataElement origProdRoot = AbstractMetadata.getOriginalProductMetadata(product);
readCalibrationLUT(sigmaLUT, lutsigma, origProdRoot, flipLUT);
readCalibrationLUT(gammaLUT, lutgamma, origProdRoot, flipLUT);
readCalibrationLUT(betaLUT, lutbeta, origProdRoot, flipLUT);
}
private static void readCalibrationLUT(
final File file, final String lutName, final MetadataElement root, final boolean flipLUT)
throws IOException {
if (!file.exists()) return;
final org.jdom.Document xmlDoc = XMLSupport.LoadXML(file.getAbsolutePath());
final Element rootElement = xmlDoc.getRootElement();
final Element offsetElem = rootElement.getChild("offset");
final double offset = Double.parseDouble(offsetElem.getValue());
final Element gainsElem = rootElement.getChild("gains");
double[] gainsArray = StringUtils.toDoubleArray(gainsElem.getValue().trim(), " ");
if (flipLUT) {
double tmp;
for (int i = 0; i < gainsArray.length / 2; i++) {
tmp = gainsArray[i];
gainsArray[i] = gainsArray[gainsArray.length - i - 1];
gainsArray[gainsArray.length - i - 1] = tmp;
}
}
final MetadataElement lut = new MetadataElement(lutName);
root.addElement(lut);
final MetadataAttribute offsetAttrib =
new MetadataAttribute("offset", ProductData.TYPE_FLOAT64);
offsetAttrib.getData().setElemDouble(offset);
lut.addAttribute(offsetAttrib);
final MetadataAttribute gainsAttrib =
new MetadataAttribute("gains", ProductData.TYPE_FLOAT64, gainsArray.length);
gainsAttrib.getData().setElems(gainsArray);
lut.addAttribute(gainsAttrib);
}
/** {@inheritDoc} */
@Override
protected void readBandRasterDataImpl(
int sourceOffsetX,
int sourceOffsetY,
int sourceWidth,
int sourceHeight,
int sourceStepX,
int sourceStepY,
Band destBand,
int destOffsetX,
int destOffsetY,
int destWidth,
int destHeight,
ProductData destBuffer,
ProgressMonitor pm)
throws IOException {
final ImageIOFile.BandInfo bandInfo = dataDir.getBandInfo(destBand);
if (bandInfo != null && bandInfo.img != null) {
if (isAscending) {
readAscendingRasterBand(
sourceOffsetX,
sourceOffsetY,
sourceStepX,
sourceStepY,
destBuffer,
destOffsetX,
destOffsetY,
destWidth,
destHeight,
0,
bandInfo.img,
bandInfo.bandSampleOffset,
isAntennaPointingRight);
} else {
readDescendingRasterBand(
sourceOffsetX,
sourceOffsetY,
sourceStepX,
sourceStepY,
destBuffer,
destOffsetX,
destOffsetY,
destWidth,
destHeight,
0,
bandInfo.img,
bandInfo.bandSampleOffset,
isAntennaPointingRight);
}
}
}
public void readAscendingRasterBand(
final int sourceOffsetX,
final int sourceOffsetY,
final int sourceStepX,
final int sourceStepY,
final ProductData destBuffer,
final int destOffsetX,
final int destOffsetY,
final int destWidth,
final int destHeight,
final int imageID,
final ImageIOFile img,
final int bandSampleOffset,
final boolean isAntennaPointingRight)
throws IOException {
final Raster data;
synchronized (dataDir) {
final ImageReader reader = img.getReader();
final ImageReadParam param = reader.getDefaultReadParam();
param.setSourceSubsampling(
sourceStepX, sourceStepY, sourceOffsetX % sourceStepX, sourceOffsetY % sourceStepY);
final RenderedImage image = reader.readAsRenderedImage(0, param);
if (flipToSARGeometry) {
if (isAntennaPointingRight) { // flip the image up side down
data =
image.getData(
new Rectangle(
destOffsetX,
img.getSceneHeight() - destOffsetY - destHeight,
destWidth,
destHeight));
} else { // flip the image upside down, then flip it left to right
data =
image.getData(
new Rectangle(
img.getSceneWidth() - destOffsetX - destWidth,
img.getSceneHeight() - destOffsetY - destHeight,
destWidth,
destHeight));
}
} else {
data = image.getData(new Rectangle(destOffsetX, destOffsetY, destWidth, destHeight));
}
}
final DataBuffer dataBuffer = data.getDataBuffer();
final SampleModel sampleModel = data.getSampleModel();
final int sampleOffset = imageID + bandSampleOffset;
if (flipToSARGeometry) {
final int[] dArray = new int[destWidth * destHeight];
sampleModel.getSamples(
0, 0, destWidth, destHeight, imageID + bandSampleOffset, dArray, dataBuffer);
int srcStride, destStride;
if (isAntennaPointingRight) { // flip the image upside down
for (int r = 0; r < destHeight; r++) {
srcStride = r * destWidth;
destStride = (destHeight - r - 1) * destWidth;
for (int c = 0; c < destWidth; c++) {
destBuffer.setElemIntAt(destStride + c, dArray[srcStride + c]);
}
}
} else { // flip the image upside down, then flip it left to right
for (int r = 0; r < destHeight; r++) {
srcStride = r * destWidth;
destStride = (destHeight - r) * destWidth;
for (int c = 0; c < destWidth; c++) {
destBuffer.setElemIntAt(destStride - c - 1, dArray[srcStride + c]);
}
}
}
} else { // no flipping is needed
sampleModel.getSamples(
0, 0, destWidth, destHeight, sampleOffset, (int[]) destBuffer.getElems(), dataBuffer);
}
}
public void readDescendingRasterBand(
final int sourceOffsetX,
final int sourceOffsetY,
final int sourceStepX,
final int sourceStepY,
final ProductData destBuffer,
final int destOffsetX,
final int destOffsetY,
final int destWidth,
final int destHeight,
final int imageID,
final ImageIOFile img,
final int bandSampleOffset,
final boolean isAntennaPointingRight)
throws IOException {
final Raster data;
synchronized (dataDir) {
final ImageReader reader = img.getReader();
final ImageReadParam param = reader.getDefaultReadParam();
param.setSourceSubsampling(
sourceStepX, sourceStepY, sourceOffsetX % sourceStepX, sourceOffsetY % sourceStepY);
final RenderedImage image = reader.readAsRenderedImage(0, param);
if (flipToSARGeometry && isAntennaPointingRight) { // flip the image left to right
data =
image.getData(
new Rectangle(
img.getSceneWidth() - destOffsetX - destWidth,
destOffsetY,
destWidth,
destHeight));
} else {
data = image.getData(new Rectangle(destOffsetX, destOffsetY, destWidth, destHeight));
}
}
final DataBuffer dataBuffer = data.getDataBuffer();
final SampleModel sampleModel = data.getSampleModel();
final int sampleOffset = imageID + bandSampleOffset;
if (flipToSARGeometry && isAntennaPointingRight) { // flip the image left to right
final int[] dArray = new int[destWidth * destHeight];
sampleModel.getSamples(0, 0, destWidth, destHeight, sampleOffset, dArray, dataBuffer);
int srcStride, destStride;
for (int r = 0; r < destHeight; r++) {
srcStride = r * destWidth;
destStride = r * destWidth + destWidth;
for (int c = 0; c < destWidth; c++) {
destBuffer.setElemIntAt(destStride - c - 1, dArray[srcStride + c]);
}
}
} else { // no flipping is needed
sampleModel.getSamples(
0, 0, destWidth, destHeight, sampleOffset, (int[]) destBuffer.getElems(), dataBuffer);
}
}
}