assertTrue(driver.canProcess(DriverCapabilities.CONNECT, source, null));
LOGGER.info("ACCEPTED: " + source.toString());
// getting access to the file
CoverageAccess access = null;
try {
access = driver.process(DriverCapabilities.CONNECT, source, null, null, null);
if (access == null) {
throw new IOException("Unable to connect");
}
// get the names
final List<Name> names = access.getNames(null);
for (Name name : names) {
// get a source
final CoverageSource gridSource = access.access(name, null, AccessType.READ_ONLY, null, null);
if (gridSource == null) {
throw new IOException("Unable to access");
}
LOGGER.info("Connected to coverage: " + name.toString());
// TEMPORAL DOMAIN
final TemporalDomain temporalDomain = gridSource.getTemporalDomain();
if (temporalDomain == null) {
LOGGER.info("Temporal domain is null");
} else {
// temporal crs
LOGGER.info("TemporalCRS: " + temporalDomain.getCoordinateReferenceSystem());
// print the temporal domain elements
for (DateRange tg : temporalDomain.getTemporalElements(true, null)) {
LOGGER.info("Global Temporal Domain: " + tg.toString());
}
// print the temporal domain elements with overall = true
StringBuilder overallTemporal = new StringBuilder("Temporal domain element (overall = true):\n");
for (DateRange tg : temporalDomain.getTemporalElements(false, null)) {
overallTemporal.append(tg.toString()).append("\n");
}
LOGGER.info(overallTemporal.toString());
}
// VERTICAL DOMAIN
final VerticalDomain verticalDomain = gridSource.getVerticalDomain();
if (verticalDomain == null) {
LOGGER.info("Vertical domain is null");
} else {
// vertical crs
LOGGER.info("VerticalCRS: " + verticalDomain.getCoordinateReferenceSystem());
// print the Vertical domain elements
for (NumberRange<Double> vg : verticalDomain.getVerticalElements(true, null)) {
LOGGER.info("Vertical domain element: " + vg.toString());
}
// print the Vertical domain elements with overall = true
StringBuilder overallVertical = new StringBuilder("Vertical domain element (overall = true):\n");
for (NumberRange<Double> vg : verticalDomain.getVerticalElements(false, null)) {
overallVertical.append(vg.toString()).append("\n");
}
LOGGER.info(overallVertical.toString());
}
// HORIZONTAL DOMAIN
final SpatialDomain spatialDomain = gridSource.getSpatialDomain();
if (spatialDomain == null) {
LOGGER.info("Horizontal domain is null");
} else {
// print the horizontal domain elements
final CoordinateReferenceSystem crs2D = spatialDomain.getCoordinateReferenceSystem2D();
assert crs2D != null;
final MathTransform2D g2w = spatialDomain.getGridToWorldTransform(null);
assert g2w != null;
final Set<? extends BoundingBox> spatialElements = spatialDomain.getSpatialElements(true, null);
assert spatialElements != null && !spatialElements.isEmpty();
final StringBuilder buf = new StringBuilder();
buf.append("Horizontal domain is as follows:\n");
buf.append("G2W:").append("\t").append(g2w).append("\n");
buf.append("CRS2D:").append("\t").append(crs2D).append("\n");
for (BoundingBox bbox : spatialElements) {
buf.append("BBOX:").append("\t").append(bbox).append("\n");
}
LOGGER.info(buf.toString());
}
CoverageReadRequest readRequest = new CoverageReadRequest();
// //
//
// Setting up a limited range for the request.
//
// //
LinkedHashSet<NumberRange<Double>> requestedVerticalSubset = new LinkedHashSet<NumberRange<Double>>();
SortedSet<? extends NumberRange<Double>> verticalElements = verticalDomain
.getVerticalElements(false, null);
final int numLevels = verticalElements.size();
final Iterator<? extends NumberRange<Double>> iterator = verticalElements.iterator();
for (int i = 0; i < numLevels; i++) {
NumberRange<Double> level = iterator.next();
if (i % (numLevels / 5) == 1) {
requestedVerticalSubset.add(level);
}
}
readRequest.setVerticalSubset(requestedVerticalSubset);
SortedSet<DateRange> requestedTemporalSubset = new DateRangeTreeSet();
SortedSet<? extends DateRange> temporalElements = temporalDomain.getTemporalElements(false, null);
final int numTimes = temporalElements.size();
Iterator<? extends DateRange> iteratorT = temporalElements.iterator();
for (int i = 0; i < numTimes; i++) {
DateRange time = iteratorT.next();
if (i % (numTimes / 5) == 1) {
requestedTemporalSubset.add(time);
}
}
readRequest.setTemporalSubset(requestedTemporalSubset);
CoverageResponse response = gridSource.read(readRequest, null);
if (response == null || response.getStatus() != Status.SUCCESS || !response.getExceptions().isEmpty()) {
throw new IOException("Unable to read");
}
final Collection<? extends Coverage> results = response.getResults(null);
int index = 0;
for (Coverage c : results) {
GridCoverageResponse resp = (GridCoverageResponse) c;
GridCoverage2D coverage = resp.getGridCoverage2D();
String title = coverage.getSampleDimension(0).getDescription().toString();
// Crs and envelope
if (isInteractiveTest) {
// ImageIOUtilities.visualize(coverage.getRenderedImage(), "tt",true);
coverage.show(title + " " + index++);
} else {
PlanarImage.wrapRenderedImage(coverage.getRenderedImage()).getTiles();
}
final StringBuilder buffer = new StringBuilder();
buffer.append("GridCoverage CRS: ")
.append(coverage.getCoordinateReferenceSystem2D().toWKT())
.append("\n");
buffer.append("GridCoverage GG: ")
.append(coverage.getGridGeometry().toString()).append("\n");
LOGGER.info(buffer.toString());
}
gridSource.dispose();
}
} catch (Throwable t) {
if (LOGGER.isLoggable(Level.WARNING)) {
LOGGER.log(Level.WARNING, t.getLocalizedMessage(), t);
}
} finally {
if (access != null) {
try {
access.dispose();
} catch (Throwable t) {
// Does nothing
}
}
}