@Override protected void addAbstractedMetadataHeader(final MetadataElement root) throws IOException { final MetadataElement absRoot = AbstractMetadata.addAbstractedMetadataHeader(root); final MetadataElement origProdRoot = AbstractMetadata.addOriginalProductMetadata(root); addManifestMetadata(getProductName(), absRoot, origProdRoot, false); acqMode = absRoot.getAttributeString(AbstractMetadata.ACQUISITION_MODE); setSLC(absRoot.getAttributeString(AbstractMetadata.SAMPLE_TYPE).equals("COMPLEX")); // get metadata for each band addBandAbstractedMetadata(absRoot, origProdRoot); addCalibrationAbstractedMetadata(origProdRoot); addNoiseAbstractedMetadata(origProdRoot); }
public static ProductData.UTC getTime(final MetadataElement elem, final String tag) { String start = elem.getAttributeString(tag, AbstractMetadata.NO_METADATA_STRING); start = start.replace("T", "_"); return AbstractMetadata.parseUTC(start, Sentinel1Constants.sentinelDateFormat); }
private static MetadataElement findElementContaining( final MetadataElement parent, final String elemName, final String attribName, final String attValue) { final MetadataElement[] elems = parent.getElements(); for (MetadataElement elem : elems) { if (elem.getName().equalsIgnoreCase(elemName) && elem.containsAttribute(attribName)) { String value = elem.getAttributeString(attribName); if (value != null && value.equalsIgnoreCase(attValue)) return elem; } } return null; }
@Override protected void addGeoCoding(final Product product) { TiePointGrid latGrid = product.getTiePointGrid(OperatorUtils.TPG_LATITUDE); TiePointGrid lonGrid = product.getTiePointGrid(OperatorUtils.TPG_LONGITUDE); if (latGrid != null && lonGrid != null) { setLatLongMetadata(product, latGrid, lonGrid); final TiePointGeoCoding tpGeoCoding = new TiePointGeoCoding(latGrid, lonGrid, Datum.WGS_84); product.setGeoCoding(tpGeoCoding); return; } final MetadataElement absRoot = AbstractMetadata.getAbstractedMetadata(product); final String acquisitionMode = absRoot.getAttributeString(AbstractMetadata.ACQUISITION_MODE); int numOfSubSwath; switch (acquisitionMode) { case "IW": numOfSubSwath = 3; break; case "EW": numOfSubSwath = 5; break; default: return; } String[] bandNames = product.getBandNames(); Band firstSWBand = null, lastSWBand = null; boolean firstSWBandFound = false, lastSWBandFound = false; for (String bandName : bandNames) { if (!firstSWBandFound && bandName.contains(acquisitionMode + 1)) { firstSWBand = product.getBand(bandName); firstSWBandFound = true; } if (!lastSWBandFound && bandName.contains(acquisitionMode + numOfSubSwath)) { lastSWBand = product.getBand(bandName); lastSWBandFound = true; } } if (firstSWBand == null && lastSWBand == null) return; final GeoCoding firstSWBandGeoCoding = bandGeocodingMap.get(firstSWBand); final int firstSWBandHeight = firstSWBand.getRasterHeight(); final GeoCoding lastSWBandGeoCoding = bandGeocodingMap.get(lastSWBand); final int lastSWBandWidth = lastSWBand.getRasterWidth(); final int lastSWBandHeight = lastSWBand.getRasterHeight(); final PixelPos ulPix = new PixelPos(0, 0); final PixelPos llPix = new PixelPos(0, firstSWBandHeight - 1); final GeoPos ulGeo = new GeoPos(); final GeoPos llGeo = new GeoPos(); firstSWBandGeoCoding.getGeoPos(ulPix, ulGeo); firstSWBandGeoCoding.getGeoPos(llPix, llGeo); final PixelPos urPix = new PixelPos(lastSWBandWidth - 1, 0); final PixelPos lrPix = new PixelPos(lastSWBandWidth - 1, lastSWBandHeight - 1); final GeoPos urGeo = new GeoPos(); final GeoPos lrGeo = new GeoPos(); lastSWBandGeoCoding.getGeoPos(urPix, urGeo); lastSWBandGeoCoding.getGeoPos(lrPix, lrGeo); final float[] latCorners = { (float) ulGeo.getLat(), (float) urGeo.getLat(), (float) llGeo.getLat(), (float) lrGeo.getLat() }; final float[] lonCorners = { (float) ulGeo.getLon(), (float) urGeo.getLon(), (float) llGeo.getLon(), (float) lrGeo.getLon() }; ReaderUtils.addGeoCoding(product, latCorners, lonCorners); AbstractMetadata.setAttribute(absRoot, AbstractMetadata.first_near_lat, ulGeo.getLat()); AbstractMetadata.setAttribute(absRoot, AbstractMetadata.first_near_long, ulGeo.getLon()); AbstractMetadata.setAttribute(absRoot, AbstractMetadata.first_far_lat, urGeo.getLat()); AbstractMetadata.setAttribute(absRoot, AbstractMetadata.first_far_long, urGeo.getLon()); AbstractMetadata.setAttribute(absRoot, AbstractMetadata.last_near_lat, llGeo.getLat()); AbstractMetadata.setAttribute(absRoot, AbstractMetadata.last_near_long, llGeo.getLon()); AbstractMetadata.setAttribute(absRoot, AbstractMetadata.last_far_lat, lrGeo.getLat()); AbstractMetadata.setAttribute(absRoot, AbstractMetadata.last_far_long, lrGeo.getLon()); // add band geocoding final Band[] bands = product.getBands(); for (Band band : bands) { band.setGeoCoding(bandGeocodingMap.get(band)); } }
private void addBandAbstractedMetadata( final MetadataElement absRoot, final MetadataElement origProdRoot) throws IOException { MetadataElement annotationElement = origProdRoot.getElement("annotation"); if (annotationElement == null) { annotationElement = new MetadataElement("annotation"); origProdRoot.addElement(annotationElement); } // collect range and azimuth spacing double rangeSpacingTotal = 0; double azimuthSpacingTotal = 0; boolean commonMetadataRetrieved = false; double heightSum = 0.0; int numBands = 0; final String annotFolder = getRootFolder() + "annotation"; final String[] filenames = listFiles(annotFolder); if (filenames != null) { for (String metadataFile : filenames) { final Document xmlDoc = XMLSupport.LoadXML(getInputStream(annotFolder + '/' + metadataFile)); final Element rootElement = xmlDoc.getRootElement(); final MetadataElement nameElem = new MetadataElement(metadataFile); annotationElement.addElement(nameElem); AbstractMetadataIO.AddXMLMetadata(rootElement, nameElem); final MetadataElement prodElem = nameElem.getElement("product"); final MetadataElement adsHeader = prodElem.getElement("adsHeader"); final String swath = adsHeader.getAttributeString("swath"); final String pol = adsHeader.getAttributeString("polarisation"); final ProductData.UTC startTime = getTime(adsHeader, "startTime"); final ProductData.UTC stopTime = getTime(adsHeader, "stopTime"); final String bandRootName = AbstractMetadata.BAND_PREFIX + swath + '_' + pol; final MetadataElement bandAbsRoot = AbstractMetadata.addBandAbstractedMetadata(absRoot, bandRootName); final String imgName = FileUtils.exchangeExtension(metadataFile, ".tiff"); imgBandMetadataMap.put(imgName, bandRootName); AbstractMetadata.setAttribute(bandAbsRoot, AbstractMetadata.SWATH, swath); AbstractMetadata.setAttribute(bandAbsRoot, AbstractMetadata.polarization, pol); AbstractMetadata.setAttribute(bandAbsRoot, AbstractMetadata.annotation, metadataFile); AbstractMetadata.setAttribute(bandAbsRoot, AbstractMetadata.first_line_time, startTime); AbstractMetadata.setAttribute(bandAbsRoot, AbstractMetadata.last_line_time, stopTime); if (AbstractMetadata.isNoData(absRoot, AbstractMetadata.mds1_tx_rx_polar)) { AbstractMetadata.setAttribute(absRoot, AbstractMetadata.mds1_tx_rx_polar, pol); } else { AbstractMetadata.setAttribute(absRoot, AbstractMetadata.mds2_tx_rx_polar, pol); } final MetadataElement imageAnnotation = prodElem.getElement("imageAnnotation"); final MetadataElement imageInformation = imageAnnotation.getElement("imageInformation"); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.data_take_id, Integer.parseInt(adsHeader.getAttributeString("missionDataTakeId"))); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.slice_num, Integer.parseInt(imageInformation.getAttributeString("sliceNumber"))); rangeSpacingTotal += imageInformation.getAttributeDouble("rangePixelSpacing"); azimuthSpacingTotal += imageInformation.getAttributeDouble("azimuthPixelSpacing"); AbstractMetadata.setAttribute( bandAbsRoot, AbstractMetadata.line_time_interval, imageInformation.getAttributeDouble("azimuthTimeInterval")); AbstractMetadata.setAttribute( bandAbsRoot, AbstractMetadata.num_samples_per_line, imageInformation.getAttributeInt("numberOfSamples")); AbstractMetadata.setAttribute( bandAbsRoot, AbstractMetadata.num_output_lines, imageInformation.getAttributeInt("numberOfLines")); AbstractMetadata.setAttribute( bandAbsRoot, AbstractMetadata.sample_type, imageInformation.getAttributeString("pixelValue").toUpperCase()); heightSum += getBandTerrainHeight(prodElem); if (!commonMetadataRetrieved) { // these should be the same for all swaths // set to absRoot final MetadataElement generalAnnotation = prodElem.getElement("generalAnnotation"); final MetadataElement productInformation = generalAnnotation.getElement("productInformation"); final MetadataElement processingInformation = imageAnnotation.getElement("processingInformation"); final MetadataElement swathProcParamsList = processingInformation.getElement("swathProcParamsList"); final MetadataElement swathProcParams = swathProcParamsList.getElement("swathProcParams"); final MetadataElement rangeProcessing = swathProcParams.getElement("rangeProcessing"); final MetadataElement azimuthProcessing = swathProcParams.getElement("azimuthProcessing"); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.range_sampling_rate, productInformation.getAttributeDouble("rangeSamplingRate") / Constants.oneMillion); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.radar_frequency, productInformation.getAttributeDouble("radarFrequency") / Constants.oneMillion); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.line_time_interval, imageInformation.getAttributeDouble("azimuthTimeInterval")); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.slant_range_to_first_pixel, imageInformation.getAttributeDouble("slantRangeTime") * Constants.halfLightSpeed); final MetadataElement downlinkInformationList = generalAnnotation.getElement("downlinkInformationList"); final MetadataElement downlinkInformation = downlinkInformationList.getElement("downlinkInformation"); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.pulse_repetition_frequency, downlinkInformation.getAttributeDouble("prf")); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.range_bandwidth, rangeProcessing.getAttributeDouble("processingBandwidth") / Constants.oneMillion); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.azimuth_bandwidth, azimuthProcessing.getAttributeDouble("processingBandwidth")); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.range_looks, rangeProcessing.getAttributeDouble("numberOfLooks")); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.azimuth_looks, azimuthProcessing.getAttributeDouble("numberOfLooks")); if (!isTOPSAR() || !isSLC()) { AbstractMetadata.setAttribute( absRoot, AbstractMetadata.num_output_lines, imageInformation.getAttributeInt("numberOfLines")); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.num_samples_per_line, imageInformation.getAttributeInt("numberOfSamples")); } addOrbitStateVectors(absRoot, generalAnnotation.getElement("orbitList")); addSRGRCoefficients(absRoot, prodElem.getElement("coordinateConversion")); addDopplerCentroidCoefficients(absRoot, prodElem.getElement("dopplerCentroid")); commonMetadataRetrieved = true; } ++numBands; } } // set average to absRoot AbstractMetadata.setAttribute( absRoot, AbstractMetadata.range_spacing, rangeSpacingTotal / (double) numBands); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.azimuth_spacing, azimuthSpacingTotal / (double) numBands); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.avg_scene_height, heightSum / filenames.length); }
static void addManifestMetadata( final String productName, final MetadataElement absRoot, final MetadataElement origProdRoot, boolean isOCN) { final String defStr = AbstractMetadata.NO_METADATA_STRING; final int defInt = AbstractMetadata.NO_METADATA; final MetadataElement XFDU = origProdRoot.getElement("XFDU"); final MetadataElement informationPackageMap = XFDU.getElement("informationPackageMap"); final MetadataElement contentUnit = informationPackageMap.getElement("contentUnit"); AbstractMetadata.setAttribute(absRoot, AbstractMetadata.PRODUCT, productName); final String descriptor = contentUnit.getAttributeString("textInfo", defStr); AbstractMetadata.setAttribute(absRoot, AbstractMetadata.SPH_DESCRIPTOR, descriptor); AbstractMetadata.setAttribute(absRoot, AbstractMetadata.antenna_pointing, "right"); final MetadataElement metadataSection = XFDU.getElement("metadataSection"); final MetadataElement[] metadataObjectList = metadataSection.getElements(); for (MetadataElement metadataObject : metadataObjectList) { final String id = metadataObject.getAttributeString("ID", defStr); if (id.endsWith("Annotation") || id.endsWith("Schema")) { // continue; } else if (id.equals("processing")) { final MetadataElement processing = findElement(metadataObject, "processing"); final MetadataElement facility = processing.getElement("facility"); final MetadataElement software = facility.getElement("software"); final String org = facility.getAttributeString("organisation"); final String name = software.getAttributeString("name"); final String version = software.getAttributeString("version"); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.ProcessingSystemIdentifier, org + ' ' + name + ' ' + version); final ProductData.UTC start = getTime(processing, "start"); AbstractMetadata.setAttribute(absRoot, AbstractMetadata.PROC_TIME, start); } else if (id.equals("acquisitionPeriod")) { final MetadataElement acquisitionPeriod = findElement(metadataObject, "acquisitionPeriod"); final ProductData.UTC startTime = getTime(acquisitionPeriod, "startTime"); final ProductData.UTC stopTime = getTime(acquisitionPeriod, "stopTime"); AbstractMetadata.setAttribute(absRoot, AbstractMetadata.first_line_time, startTime); AbstractMetadata.setAttribute(absRoot, AbstractMetadata.last_line_time, stopTime); } else if (id.equals("platform")) { final MetadataElement platform = findElement(metadataObject, "platform"); String missionName = platform.getAttributeString("familyName", "Sentinel-1"); final String number = platform.getAttributeString("number", defStr); if (!missionName.equals("ENVISAT")) missionName += number; AbstractMetadata.setAttribute(absRoot, AbstractMetadata.MISSION, missionName); final MetadataElement instrument = platform.getElement("instrument"); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.SWATH, instrument.getAttributeString("swath", defStr)); String acqMode = instrument.getAttributeString("mode", defStr); if (acqMode == null || acqMode.equals(defStr)) { final MetadataElement extensionElem = instrument.getElement("extension"); if (extensionElem != null) { final MetadataElement instrumentModeElem = extensionElem.getElement("instrumentMode"); if (instrumentModeElem != null) acqMode = instrumentModeElem.getAttributeString("mode", defStr); } } AbstractMetadata.setAttribute(absRoot, AbstractMetadata.ACQUISITION_MODE, acqMode); } else if (id.equals("measurementOrbitReference")) { final MetadataElement orbitReference = findElement(metadataObject, "orbitReference"); final MetadataElement orbitNumber = findElementContaining(orbitReference, "OrbitNumber", "type", "start"); final MetadataElement relativeOrbitNumber = findElementContaining(orbitReference, "relativeOrbitNumber", "type", "start"); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.ABS_ORBIT, orbitNumber.getAttributeInt("orbitNumber", defInt)); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.REL_ORBIT, relativeOrbitNumber.getAttributeInt("relativeOrbitNumber", defInt)); AbstractMetadata.setAttribute( absRoot, AbstractMetadata.CYCLE, orbitReference.getAttributeInt("cycleNumber", defInt)); String pass = orbitReference.getAttributeString("pass", defStr); if (pass.equals(defStr)) { final MetadataElement extension = orbitReference.getElement("extension"); final MetadataElement orbitProperties = extension.getElement("orbitProperties"); pass = orbitProperties.getAttributeString("pass", defStr); } AbstractMetadata.setAttribute(absRoot, AbstractMetadata.PASS, pass); } else if (id.equals("measurementFrameSet")) { } else if (id.equals("generalProductInformation")) { MetadataElement generalProductInformation = findElement(metadataObject, "generalProductInformation"); if (generalProductInformation == null) generalProductInformation = findElement(metadataObject, "standAloneProductInformation"); String productType = "unknown"; if (isOCN) { productType = "OCN"; } else { if (generalProductInformation != null) productType = generalProductInformation.getAttributeString("productType", defStr); } AbstractMetadata.setAttribute(absRoot, AbstractMetadata.PRODUCT_TYPE, productType); if (productType.contains("SLC")) { AbstractMetadata.setAttribute(absRoot, AbstractMetadata.SAMPLE_TYPE, "COMPLEX"); } else { AbstractMetadata.setAttribute(absRoot, AbstractMetadata.SAMPLE_TYPE, "DETECTED"); AbstractMetadata.setAttribute(absRoot, AbstractMetadata.srgr_flag, 1); } } } }
@Override protected void addBands(final Product product) { String bandName; boolean real = true; Band lastRealBand = null; String unit; final MetadataElement absRoot = AbstractMetadata.getAbstractedMetadata(product); for (Map.Entry<String, ImageIOFile> stringImageIOFileEntry : bandImageFileMap.entrySet()) { final ImageIOFile img = stringImageIOFileEntry.getValue(); final String imgName = img.getName().toLowerCase(); final MetadataElement bandMetadata = absRoot.getElement(imgBandMetadataMap.get(imgName)); final String swath = bandMetadata.getAttributeString(AbstractMetadata.swath); final String pol = bandMetadata.getAttributeString(AbstractMetadata.polarization); final int width = bandMetadata.getAttributeInt(AbstractMetadata.num_samples_per_line); final int height = bandMetadata.getAttributeInt(AbstractMetadata.num_output_lines); String tpgPrefix = ""; String suffix = pol; if (isSLC() && isTOPSAR()) { suffix = swath + '_' + pol; tpgPrefix = swath; } int numImages = img.getNumImages(); if (isSLC()) { numImages *= 2; // real + imaginary } for (int i = 0; i < numImages; ++i) { if (isSLC()) { for (int b = 0; b < img.getNumBands(); ++b) { if (real) { bandName = "i" + '_' + suffix; unit = Unit.REAL; } else { bandName = "q" + '_' + suffix; unit = Unit.IMAGINARY; } final Band band = new Band(bandName, ProductData.TYPE_INT16, width, height); band.setUnit(unit); product.addBand(band); bandMap.put(band, new ImageIOFile.BandInfo(band, img, i, b)); AbstractMetadata.addBandToBandMap(bandMetadata, bandName); if (real) { lastRealBand = band; } else { ReaderUtils.createVirtualIntensityBand(product, lastRealBand, band, '_' + suffix); } real = !real; // add tiepointgrids and geocoding for band addTiePointGrids(band, imgName, tpgPrefix); // reset to null so it doesn't adopt a geocoding from the bands product.setGeoCoding(null); } } else { for (int b = 0; b < img.getNumBands(); ++b) { bandName = "Amplitude" + '_' + suffix; final Band band = new Band(bandName, ProductData.TYPE_INT32, width, height); band.setUnit(Unit.AMPLITUDE); product.addBand(band); bandMap.put(band, new ImageIOFile.BandInfo(band, img, i, b)); AbstractMetadata.addBandToBandMap(bandMetadata, bandName); SARReader.createVirtualIntensityBand(product, band, '_' + suffix); // add tiepointgrids and geocoding for band addTiePointGrids(band, imgName, tpgPrefix); } } } } }