private void setSceneRasterStartAndStopTime(Product product) {
final Product sourceProduct = getSourceProduct();
final ProductData.UTC startTime = sourceProduct.getStartTime();
final ProductData.UTC stopTime = sourceProduct.getEndTime();
final ProductSubsetDef subsetDef = getSubsetDef();
if (startTime != null
&& stopTime != null
&& subsetDef != null
&& subsetDef.getRegion() != null) {
final double height = sourceProduct.getSceneRasterHeight();
final Rectangle region = subsetDef.getRegion();
final double regionY = region.getY();
final double regionHeight = region.getHeight();
final double dStart = startTime.getMJD();
final double dStop = stopTime.getMJD();
final double vPerLine = (dStop - dStart) / (height - 1);
final double newStart = vPerLine * regionY + dStart;
final double newStop = vPerLine * (regionHeight - 1) + newStart;
product.setStartTime(new ProductData.UTC(newStart));
product.setEndTime(new ProductData.UTC(newStop));
} else {
product.setStartTime(startTime);
product.setEndTime(stopTime);
}
}
@Override
public void decode(Modis35ProfileReadContext ctx, Product p) throws IOException {
NetcdfFile ncFile = ctx.getNetcdfFile();
p.setStartTime(
TimeUtils.getSceneRasterTime(
ncFile, Modis35Constants.START_DATE_ATT_NAME, Modis35Constants.START_TIME_ATT_NAME));
p.setEndTime(
TimeUtils.getSceneRasterTime(
ncFile, Modis35Constants.STOP_DATE_ATT_NAME, Modis35Constants.STOP_TIME_ATT_NAME));
}
static boolean setStartTime(final Product product) {
if (product.getStartTime() == null) {
final ProductData.UTC utc = getTimeAttrValue(product, ProductAttributes.START_TIME);
if (utc != null) {
product.setStartTime(utc);
return true;
}
}
return false;
}
protected Product createProduct(
GCTileFile refGcFile, String prodName, String prodType, int width, int height)
throws IOException {
final Product product = new Product(prodName, prodType, width, height);
product.setFileLocation(new File(refGcFile.getFilePath()));
product.setStartTime(refGcFile.getStartDate());
product.setEndTime(refGcFile.getEndDate());
addGeoCoding(product);
addBands(product, refGcFile);
addIndexCodingAndBitmasks(product.getBand("SM"));
refGcFile.readMetadata(product.getMetadataRoot());
return product;
}
private static Product createProduct() throws Exception {
Product product = new Product("p", "t", 32, 256);
final TiePointGrid lat =
new TiePointGrid("lat", 2, 2, 0f, 0f, 32f, 256f, new float[] {+40f, +40f, -40f, -40f});
final TiePointGrid lon =
new TiePointGrid("lon", 2, 2, 0f, 0f, 32f, 256f, new float[] {-80f, +80f, -80f, +80f});
product.addTiePointGrid(lat);
product.addTiePointGrid(lon);
product.setGeoCoding(new TiePointGeoCoding(lat, lon));
product.setPreferredTileSize(32, 16);
product.setStartTime(ProductData.UTC.parse("2003-01-01", "yyyy-MM-dd"));
product.setEndTime(ProductData.UTC.parse("2003-01-02", "yyyy-MM-dd"));
return product;
}
@Override
protected Product readProductNodesImpl() throws IOException {
final String s = getInput().toString();
final File file0 = new File(s);
final File dir = file0.getParentFile();
final S2FilenameInfo fni0 = S2FilenameInfo.create(file0.getName());
if (fni0 == null) {
throw new IOException();
}
Header metadataHeader = null;
final Map<Integer, BandInfo> fileMap = new HashMap<Integer, BandInfo>();
if (dir != null) {
File[] files =
dir.listFiles(
new FilenameFilter() {
@Override
public boolean accept(File dir, String name) {
return name.endsWith(Sentinel2ProductReaderPlugIn.JP2_EXT);
}
});
if (files != null) {
for (File file : files) {
int bandIndex = fni0.getBand(file.getName());
if (bandIndex >= 0 && bandIndex < WAVEBAND_INFOS.length) {
final S2WavebandInfo wavebandInfo = WAVEBAND_INFOS[bandIndex];
BandInfo bandInfo =
new BandInfo(
file, bandIndex, wavebandInfo, imageLayouts[wavebandInfo.resolution.id]);
fileMap.put(bandIndex, bandInfo);
}
}
}
File[] metadataFiles =
dir.listFiles(
new FilenameFilter() {
@Override
public boolean accept(File dir, String name) {
return name.startsWith("MTD_") && name.endsWith(".xml");
}
});
if (metadataFiles != null && metadataFiles.length > 0) {
File metadataFile = metadataFiles[0];
try {
metadataHeader = Header.parseHeader(metadataFile);
} catch (JDOMException e) {
BeamLogManager.getSystemLogger()
.warning("Failed to parse metadata file: " + metadataFile);
}
} else {
BeamLogManager.getSystemLogger().warning("No metadata file found");
}
}
final ArrayList<Integer> bandIndexes = new ArrayList<Integer>(fileMap.keySet());
Collections.sort(bandIndexes);
if (bandIndexes.isEmpty()) {
throw new IOException("No valid bands found.");
}
String prodType = "S2_MSI_" + fni0.procLevel;
final Product product =
new Product(
String.format("%s_%s_%s", prodType, fni0.orbitNo, fni0.tileId),
prodType,
imageLayouts[S2Resolution.R10M.id].width,
imageLayouts[S2Resolution.R10M.id].height);
try {
product.setStartTime(ProductData.UTC.parse(fni0.start, "yyyyMMddHHmmss"));
} catch (ParseException e) {
// warn
}
try {
product.setEndTime(ProductData.UTC.parse(fni0.stop, "yyyyMMddHHmmss"));
} catch (ParseException e) {
// warn
}
if (metadataHeader != null) {
SceneDescription sceneDescription = SceneDescription.create(metadataHeader);
int tileIndex = sceneDescription.getTileIndex(fni0.tileId);
Envelope2D tileEnvelope = sceneDescription.getTileEnvelope(tileIndex);
Header.Tile tile = metadataHeader.getTileList().get(tileIndex);
try {
product.setGeoCoding(
new CrsGeoCoding(
tileEnvelope.getCoordinateReferenceSystem(),
imageLayouts[S2Resolution.R10M.id].width,
imageLayouts[S2Resolution.R10M.id].height,
tile.tileGeometry10M.upperLeftX,
tile.tileGeometry10M.upperLeftY,
tile.tileGeometry10M.xDim,
-tile.tileGeometry10M.yDim,
0.0,
0.0));
} catch (FactoryException e) {
// todo - handle e
} catch (TransformException e) {
// todo - handle e
}
}
for (Integer bandIndex : bandIndexes) {
final BandInfo bandInfo = fileMap.get(bandIndex);
final Band band = product.addBand(bandInfo.wavebandInfo.bandName, ProductData.TYPE_UINT16);
band.setSpectralWavelength((float) bandInfo.wavebandInfo.centralWavelength);
band.setSpectralBandwidth((float) bandInfo.wavebandInfo.bandWidth);
band.setSpectralBandIndex(bandIndex);
band.setSourceImage(new DefaultMultiLevelImage(new Jp2MultiLevelSource(bandInfo)));
}
product.setNumResolutionLevels(imageLayouts[0].numResolutions);
return product;
}
