private static void setLatLongMetadata(
final Product product,
final MetadataElement absRoot,
final String tagLat,
final String tagLon,
final float x,
final float y) {
final PixelPos pixelPos = new PixelPos(x, y);
final GeoPos geoPos = new GeoPos();
if (product.getGeoCoding() == null) return;
product.getGeoCoding().getGeoPos(pixelPos, geoPos);
final MetadataAttribute lat = absRoot.getAttribute(tagLat);
if (lat != null) lat.getData().setElemDouble(geoPos.getLat());
final MetadataAttribute lon = absRoot.getAttribute(tagLon);
if (lon != null) lon.getData().setElemDouble(geoPos.getLon());
}
public void addGeocoding(Product product) throws IOException {
double pixelX = 0.5;
double pixelY = 0.5;
double easting;
double northing;
double pixelSizeX;
double pixelSizeY;
boolean pixelRegistered = true;
if (productReader.getProductType() == SeadasProductReader.ProductType.ANCNRT) {
pixelRegistered = false;
}
if (productReader.getProductType() == SeadasProductReader.ProductType.OISST) {
Variable lon = ncFile.findVariable("lon");
Variable lat = ncFile.findVariable("lat");
Array lonData = lon.read();
// TODO: handle the 180 degree shift with the NOAA products - need to modify
// SeaDasFileReader:readBandData
// Below is a snippet from elsewhere in BEAM that deals with this issue...
// SPECIAL CASE: check if we have a global geographic lat/lon with lon from 0..360 instead of
// -180..180
// if (isShifted180(lonData)) {
// // if this is true, subtract 180 from all longitudes and
// // add a global attribute which will be analyzed when setting up the
// image(s)
// final List<Variable> variables = ncFile.getVariables();
// for (Variable next : variables) {
// next.getAttributes().add(new Attribute("LONGITUDE_SHIFTED_180", 1));
// }
// for (int i = 0; i < lonData.getSize(); i++) {
// final Index ii = lonData.getIndex().set(i);
// final double theLon = lonData.getDouble(ii) - 180.0;
// lonData.setDouble(ii, theLon);
// }
// }
final Array latData = lat.read();
final int lonSize = lon.getShape(0);
final Index i0 = lonData.getIndex().set(0);
final Index i1 = lonData.getIndex().set(lonSize - 1);
pixelSizeX =
(lonData.getDouble(i1) - lonData.getDouble(i0)) / (product.getSceneRasterWidth() - 1);
easting = lonData.getDouble(i0);
final int latSize = lat.getShape(0);
final Index j0 = latData.getIndex().set(0);
final Index j1 = latData.getIndex().set(latSize - 1);
pixelSizeY =
(latData.getDouble(j1) - latData.getDouble(j0)) / (product.getSceneRasterHeight() - 1);
// this should be the 'normal' case
if (pixelSizeY < 0) {
pixelSizeY = -pixelSizeY;
northing = latData.getDouble(latData.getIndex().set(0));
} else {
northing = latData.getDouble(latData.getIndex().set(latSize - 1));
}
northing -= pixelSizeX / 2.0;
easting += pixelSizeY / 2.0;
try {
product.setGeoCoding(
new CrsGeoCoding(
DefaultGeographicCRS.WGS84,
product.getSceneRasterWidth(),
product.getSceneRasterHeight(),
easting,
northing,
pixelSizeX,
pixelSizeY,
pixelX,
pixelY));
} catch (FactoryException e) {
throw new IllegalStateException(e);
} catch (TransformException e) {
throw new IllegalStateException(e);
}
} else {
String east = "Easternmost_Longitude";
String west = "Westernmost_Longitude";
String north = "Northernmost_Latitude";
String south = "Southernmost_Latitude";
Attribute latmax = ncFile.findGlobalAttributeIgnoreCase("geospatial_lat_max");
if (latmax != null) {
east = "geospatial_lon_max";
west = "geospatial_lon_min";
north = "geospatial_lat_max";
south = "geospatial_lat_min";
} else {
latmax = ncFile.findGlobalAttributeIgnoreCase("upper_lat");
if (latmax != null) {
east = "right_lon";
west = "left_lon";
north = "upper_lat";
south = "lower_lat";
}
}
final MetadataElement globalAttributes =
product.getMetadataRoot().getElement("Global_Attributes");
easting = (float) globalAttributes.getAttribute(east).getData().getElemDouble();
float westing = (float) globalAttributes.getAttribute(west).getData().getElemDouble();
pixelSizeX = Math.abs(easting - westing) / product.getSceneRasterWidth();
northing = (float) globalAttributes.getAttribute(north).getData().getElemDouble();
float southing = (float) globalAttributes.getAttribute(south).getData().getElemDouble();
if (northing < southing) {
mustFlipY = true;
northing = (float) globalAttributes.getAttribute(south).getData().getElemDouble();
southing = (float) globalAttributes.getAttribute(north).getData().getElemDouble();
}
pixelSizeY = Math.abs(northing - southing) / product.getSceneRasterHeight();
if (pixelRegistered) {
northing -= pixelSizeY / 2.0;
westing += pixelSizeX / 2.0;
} else {
pixelX = 0.0;
pixelY = 0.0;
}
try {
product.setGeoCoding(
new CrsGeoCoding(
DefaultGeographicCRS.WGS84,
product.getSceneRasterWidth(),
product.getSceneRasterHeight(),
westing,
northing,
pixelSizeX,
pixelSizeY,
pixelX,
pixelY));
} catch (FactoryException e) {
throw new IllegalStateException(e);
} catch (TransformException e) {
throw new IllegalStateException(e);
}
}
}
private static void updateMetadata(
final Product sourceProduct, final Product targetProduct, ProductSubsetDef subsetDef)
throws IOException {
try {
final MetadataElement root = targetProduct.getMetadataRoot();
if (root == null) return;
final MetadataElement absRoot = root.getElement("Abstracted_Metadata");
if (absRoot == null) return;
boolean nearRangeOnLeft = isNearRangeOnLeft(targetProduct);
final MetadataAttribute firstLineTime = absRoot.getAttribute("first_line_time");
if (firstLineTime != null) {
final ProductData.UTC startTime = targetProduct.getStartTime();
if (startTime != null) firstLineTime.getData().setElems(startTime.getArray());
}
final MetadataAttribute lastLineTime = absRoot.getAttribute("last_line_time");
if (lastLineTime != null) {
final ProductData.UTC endTime = targetProduct.getEndTime();
if (endTime != null) lastLineTime.getData().setElems(endTime.getArray());
}
final MetadataAttribute totalSize = absRoot.getAttribute("total_size");
if (totalSize != null) totalSize.getData().setElemUInt(targetProduct.getRawStorageSize());
if (nearRangeOnLeft) {
setLatLongMetadata(targetProduct, absRoot, "first_near_lat", "first_near_long", 0.5f, 0.5f);
setLatLongMetadata(
targetProduct,
absRoot,
"first_far_lat",
"first_far_long",
targetProduct.getSceneRasterWidth() - 1 + 0.5f,
0.5f);
setLatLongMetadata(
targetProduct,
absRoot,
"last_near_lat",
"last_near_long",
0.5f,
targetProduct.getSceneRasterHeight() - 1 + 0.5f);
setLatLongMetadata(
targetProduct,
absRoot,
"last_far_lat",
"last_far_long",
targetProduct.getSceneRasterWidth() - 1 + 0.5f,
targetProduct.getSceneRasterHeight() - 1 + 0.5f);
} else {
setLatLongMetadata(
targetProduct,
absRoot,
"first_near_lat",
"first_near_long",
targetProduct.getSceneRasterWidth() - 1 + 0.5f,
0.5f);
setLatLongMetadata(targetProduct, absRoot, "first_far_lat", "first_far_long", 0.5f, 0.5f);
setLatLongMetadata(
targetProduct,
absRoot,
"last_near_lat",
"last_near_long",
targetProduct.getSceneRasterWidth() - 1 + 0.5f,
targetProduct.getSceneRasterHeight() - 1 + 0.5f);
setLatLongMetadata(
targetProduct,
absRoot,
"last_far_lat",
"last_far_long",
0.5f,
targetProduct.getSceneRasterHeight() - 1 + 0.5f);
}
final MetadataAttribute height = absRoot.getAttribute("num_output_lines");
if (height != null) height.getData().setElemUInt(targetProduct.getSceneRasterHeight());
final MetadataAttribute width = absRoot.getAttribute("num_samples_per_line");
if (width != null) width.getData().setElemUInt(targetProduct.getSceneRasterWidth());
final MetadataAttribute offsetX = absRoot.getAttribute("subset_offset_x");
if (offsetX != null && subsetDef.getRegion() != null)
offsetX.getData().setElemUInt(subsetDef.getRegion().x);
final MetadataAttribute offsetY = absRoot.getAttribute("subset_offset_y");
if (offsetY != null && subsetDef.getRegion() != null)
offsetY.getData().setElemUInt(subsetDef.getRegion().y);
final MetadataAttribute slantRange = absRoot.getAttribute("slant_range_to_first_pixel");
if (slantRange != null) {
final TiePointGrid srTPG = targetProduct.getTiePointGrid("slant_range_time");
if (srTPG != null) {
final double slantRangeTime;
if (nearRangeOnLeft) {
slantRangeTime = srTPG.getPixelDouble(0, 0) / 1000000000.0; // ns to s
} else {
slantRangeTime =
srTPG.getPixelDouble(targetProduct.getSceneRasterWidth() - 1, 0)
/ 1000000000.0; // ns to s
}
final double halfLightSpeed = 299792458.0 / 2.0;
final double slantRangeDist = slantRangeTime * halfLightSpeed;
slantRange.getData().setElemDouble(slantRangeDist);
}
}
setSubsetSRGRCoefficients(sourceProduct, targetProduct, subsetDef, absRoot, nearRangeOnLeft);
} catch (Exception e) {
throw new IOException(e);
}
}