/*
* GeoTools - The Open Source Java GIS Toolkit
* http://geotools.org
*
* (C) 2007 - 2016, Open Source Geospatial Foundation (OSGeo)
*
* This library is free software; you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public
* License as published by the Free Software Foundation;
* version 2.1 of the License.
*
* This library is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* Lesser General Public License for more details.
*/
package org.geotools.gce.imagemosaic;
import it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReaderSpi;
import java.awt.Rectangle;
import java.awt.RenderingHints;
import java.awt.geom.AffineTransform;
import java.awt.image.RenderedImage;
import java.io.File;
import java.io.FileNotFoundException;
import java.io.IOException;
import java.net.URL;
import javax.imageio.ImageReadParam;
import javax.imageio.spi.ImageReaderSpi;
import javax.media.jai.ImageLayout;
import javax.media.jai.JAI;
import org.geotools.coverage.grid.GridGeometry2D;
import org.geotools.coverage.grid.io.AbstractGridFormat;
import org.geotools.coverage.grid.io.footprint.MultiLevelROI;
import org.geotools.data.DataUtilities;
import org.geotools.factory.Hints;
import org.geotools.gce.imagemosaic.GranuleDescriptor.GranuleOverviewLevelDescriptor;
import org.geotools.geometry.GeneralEnvelope;
import org.geotools.geometry.jts.ReferencedEnvelope;
import org.geotools.referencing.CRS;
import org.geotools.referencing.crs.DefaultGeographicCRS;
import org.geotools.referencing.crs.DefaultProjectedCRS;
import org.geotools.referencing.operation.transform.AffineTransform2D;
import org.geotools.referencing.wkt.Parser;
import org.geotools.test.TestData;
import it.geosolutions.jaiext.utilities.ImageLayout2;
import org.junit.Assert;
import org.junit.Test;
import org.opengis.geometry.BoundingBox;
import org.opengis.parameter.GeneralParameterValue;
import org.opengis.parameter.ParameterValue;
import org.opengis.referencing.FactoryException;
import org.opengis.referencing.crs.CoordinateReferenceSystem;
import org.opengis.referencing.datum.PixelInCell;
import org.opengis.referencing.operation.NoninvertibleTransformException;
/**
* Testing {@link GranuleDescriptor} class.
*
* @author Daniele Romagnoli, GeoSolutions SAS
* @author Stefan Alfons Krueger (alfonx), Wikisquare.de : Support for jar:file:foo.jar/bar.properties URLs
*
* @source $URL$
*/
public class GranuleTest extends Assert {
private static final double DELTA = 10E-6;
private static final double DELTASCALE = 10E-2;
private static CoordinateReferenceSystem WGS84;
static{
try{
WGS84 = CRS.decode("EPSG:4326",true);
} catch (FactoryException fe){
WGS84 = DefaultGeographicCRS.WGS84;
}
}
private static final ReferencedEnvelope TEST_BBOX = new ReferencedEnvelope(12.139578206197234,15.036279855058655,40.5313698832181,42.5511689138571,WGS84);
private static final ImageReaderSpi spi = new TIFFImageReaderSpi();
public GranuleTest() {
}
@Test
public void testGranuleLevels() throws FileNotFoundException, IOException {
//get some test data
final URL testUrl= TestData.url(this, "/overview/0/D220161A.tif");
testUrl.openStream().close();
//Create a GranuleDescriptor
final GranuleDescriptor granuleDescriptor = new GranuleDescriptor(DataUtilities.urlToFile(testUrl).getAbsolutePath()
, TEST_BBOX, spi, (MultiLevelROI) null);
assertNotNull(granuleDescriptor.toString());
//Get a GranuleOverviewLevelDescriptor
final GranuleOverviewLevelDescriptor granuleOverviewLevelDescriptor = granuleDescriptor.getLevel(2);
assertNotNull(granuleOverviewLevelDescriptor);
final int h = granuleOverviewLevelDescriptor.getHeight();
final int w = granuleOverviewLevelDescriptor.getWidth();
assertEquals(47, h);
assertEquals(35, w);
final double scaleX = granuleOverviewLevelDescriptor.getScaleX();
final double scaleY = granuleOverviewLevelDescriptor.getScaleY();
assertEquals("ScaleX not equal", scaleX, 4.0d, DELTASCALE);
assertEquals("ScaleY not equal", scaleY, 3.9788d, DELTASCALE);
final Rectangle rect = granuleOverviewLevelDescriptor.getBounds();
assertEquals(rect.x, 0);
assertEquals(rect.y, 0);
assertEquals(rect.width, 35);
assertEquals(rect.height, 47);
final AffineTransform btlTransform = granuleOverviewLevelDescriptor.getBaseToLevelTransform();
final double[] baseMatrix = new double[6];
btlTransform.getMatrix(baseMatrix);
assertEquals("m00 not equal", baseMatrix[0], 4.0d, DELTASCALE);
assertEquals("m10 not equal", baseMatrix[1], 0.0d, DELTA);
assertEquals("m01 not equal", baseMatrix[2], 0.0d, DELTA);
assertEquals("m11 not equal", baseMatrix[3], 3.9788d, DELTASCALE);
assertEquals("m02 not equal", baseMatrix[4], 0.0d, DELTA);
assertEquals("m12 not equal", baseMatrix[5], 0.0d, DELTA);
final AffineTransform2D g2wtTransform = granuleOverviewLevelDescriptor.getGridToWorldTransform();
final double[] g2wMatrix = new double[6];
g2wtTransform.getMatrix(g2wMatrix);
assertEquals("m00 not equal", g2wMatrix[0], 0.08276290425318347d, DELTASCALE);
assertEquals("m10 not equal", g2wMatrix[1], 0.0d, DELTA);
assertEquals("m01 not equal", g2wMatrix[2], 0.0d, DELTA);
assertEquals("m11 not equal", g2wMatrix[3], -0.04297444746040424d, DELTASCALE);
assertEquals("m02 not equal", g2wMatrix[4], 12.139578206197234d, DELTA);
assertEquals("m12 not equal", g2wMatrix[5], 42.5511689138571d, DELTA);
}
@Test
public void testLoadRaster() throws FileNotFoundException, IOException, NoninvertibleTransformException {
//get some test data
final File testMosaic = TestData.file(this, "/rgb");
assertTrue(testMosaic.exists());
final URL testUrl= TestData.url(this, "/rgb/global_mosaic_12.png");
testUrl.openStream().close();
final GranuleDescriptor granuleDescriptor = new GranuleDescriptor(DataUtilities.urlToFile(testUrl).getAbsolutePath()
, TEST_BBOX, spi, (MultiLevelROI) null);
final GranuleOverviewLevelDescriptor granuleOverviewLevelDescriptor = granuleDescriptor.getLevel(0);
assertNotNull(granuleOverviewLevelDescriptor);
final Hints crsHints = new Hints(Hints.DEFAULT_COORDINATE_REFERENCE_SYSTEM, DefaultGeographicCRS.WGS84);
final ImageMosaicReader reader = (ImageMosaicReader) new ImageMosaicFormat().getReader(testMosaic,crsHints);
assertNotNull(reader);
final RasterManager manager = reader.getRasterManager(reader.getGridCoverageNames()[0]);
// use imageio with defined tiles
final ParameterValue<Boolean> useJai = AbstractGridFormat.USE_JAI_IMAGEREAD.createValue();
useJai.setValue(false);
final ParameterValue<String> tileSize = AbstractGridFormat.SUGGESTED_TILE_SIZE.createValue();
tileSize.setValue("10,10");
// Creating a request
final RasterLayerRequest request = new RasterLayerRequest(new GeneralParameterValue[] {useJai ,tileSize},manager);
final ImageReadParam readParameters = new ImageReadParam();
readParameters.setSourceRegion(new Rectangle(0,0,50,50));
final AffineTransform2D gridToWorldTransform = granuleOverviewLevelDescriptor.getGridToWorldTransform();
ImageLayout layout = new ImageLayout2().setTileGridXOffset(0).setTileGridYOffset(0).setTileHeight(10).setTileWidth(10);
RenderingHints rHints = new RenderingHints(JAI.KEY_IMAGE_LAYOUT, layout);
Hints hints = new Hints(rHints);
final RenderedImage raster = granuleDescriptor.loadRaster(readParameters, 0, TEST_BBOX, gridToWorldTransform.inverse(),
request, hints).getRaster();
assertEquals(raster.getWidth(), 50);
assertEquals(raster.getHeight(), 50);
AffineTransform translate = new AffineTransform(gridToWorldTransform);
translate.preConcatenate(AffineTransform.getTranslateInstance(2, 2));
final RenderedImage translatedRaster = granuleDescriptor.loadRaster(readParameters, 0, TEST_BBOX, new AffineTransform2D(translate).inverse(), request, hints).getRaster();
assertEquals(translatedRaster.getWidth(), 50);
assertEquals(translatedRaster.getHeight(), 50);
}
static final String NZTM_WKT_NE = "PROJCS[\"NZGD2000 / New Zealand Transverse Mercator 2000\", \n"
+ " GEOGCS[\"NZGD2000\", \n"
+ " DATUM[\"New Zealand Geodetic Datum 2000\", \n"
+ " SPHEROID[\"GRS 1980\", 6378137.0, 298.257222101, AUTHORITY[\"EPSG\",\"7019\"]], \n"
+ " TOWGS84[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], \n"
+ " AUTHORITY[\"EPSG\",\"6167\"]], \n"
+ " PRIMEM[\"Greenwich\", 0.0, AUTHORITY[\"EPSG\",\"8901\"]], \n"
+ " UNIT[\"degree\", 0.017453292519943295], \n"
+ " AXIS[\"Geodetic latitude\", NORTH], \n"
+ " AXIS[\"Geodetic longitude\", EAST], \n"
+ " AUTHORITY[\"EPSG\",\"4167\"]], \n"
+ " PROJECTION[\"Transverse_Mercator\", AUTHORITY[\"EPSG\",\"9807\"]], \n"
+ " PARAMETER[\"central_meridian\", 173.0], \n"
+ " PARAMETER[\"latitude_of_origin\", 0.0], \n"
+ " PARAMETER[\"scale_factor\", 0.9996], \n"
+ " PARAMETER[\"false_easting\", 1600000.0], \n"
+ " PARAMETER[\"false_northing\", 10000000.0], \n"
+ " UNIT[\"m\", 1.0], \n"
+ " AXIS[\"Northing\", NORTH], \n"
+ " AXIS[\"Easting\", EAST], \n"
+ " AUTHORITY[\"EPSG\",\"2193\"]]";
static final String NZTM_WKT_EN = "PROJCS[\"NZGD2000 / New Zealand Transverse Mercator 2000\", \n"
+ " GEOGCS[\"NZGD2000\", \n"
+ " DATUM[\"New Zealand Geodetic Datum 2000\", \n"
+ " SPHEROID[\"GRS 1980\", 6378137.0, 298.257222101, AUTHORITY[\"EPSG\",\"7019\"]], \n"
+ " TOWGS84[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], \n"
+ " AUTHORITY[\"EPSG\",\"6167\"]], \n"
+ " PRIMEM[\"Greenwich\", 0.0, AUTHORITY[\"EPSG\",\"8901\"]], \n"
+ " UNIT[\"degree\", 0.017453292519943295], \n"
+ " AXIS[\"Geodetic longitude\", EAST], \n"
+ " AXIS[\"Geodetic latitude\", NORTH], \n"
+ " AUTHORITY[\"EPSG\",\"4167\"]], \n"
+ " PROJECTION[\"Transverse_Mercator\"], \n"
+ " PARAMETER[\"central_meridian\", 173.0], \n"
+ " PARAMETER[\"latitude_of_origin\", 0.0], \n"
+ " PARAMETER[\"scale_factor\", 0.9996], \n"
+ " PARAMETER[\"false_easting\", 1600000.0], \n"
+ " PARAMETER[\"false_northing\", 10000000.0], \n"
+ " UNIT[\"m\", 1.0], \n"
+ " AXIS[\"Easting\", EAST], \n"
+ " AXIS[\"Northing\", NORTH], \n"
+ " AUTHORITY[\"EPSG\",\"2193\"]]";
@Test
public void testCRS_NorthingEasting() throws Exception {
// get some test data
final File testMosaic = TestData.file(this, "/crs_nztm");
assertTrue(testMosaic.exists());
final Parser parser = new Parser();
final DefaultProjectedCRS crs_EN = (DefaultProjectedCRS) parser.parseObject(NZTM_WKT_EN);
final DefaultProjectedCRS crs_NE = (DefaultProjectedCRS) parser.parseObject(NZTM_WKT_NE);
final ImageMosaicReader reader = (ImageMosaicReader) new ImageMosaicFormat()
.getReader(testMosaic);
assertNotNull(reader);
final RasterManager manager = reader.getRasterManager(reader.getGridCoverageNames()[0]);
// FIXME: somehow when run under JUnit the bounds end up as (y,x) rather than (x,y). Works
// fine in GeoServer. Hack it :(
final ReferencedEnvelope mosaicBounds = new ReferencedEnvelope(1587997.8835, 1612003.2265,
6162000.4515, 6198002.1165, crs_EN);
manager.spatialDomainManager.coverageBBox = mosaicBounds;
manager.spatialDomainManager.coverageEnvelope = new GeneralEnvelope(mosaicBounds);
// set up the request (north-east version)
final ReferencedEnvelope requestBBoxNE = new ReferencedEnvelope(6154440.101350001,
6204842.43235, 1583436.86902, 1617044.34782, crs_NE);
final ParameterValue<GridGeometry2D> requestedBBox = AbstractGridFormat.READ_GRIDGEOMETRY2D
.createValue();
requestedBBox.setValue(new GridGeometry2D(PixelInCell.CELL_CENTER, reader
.getOriginalGridToWorld(PixelInCell.CELL_CENTER), requestBBoxNE, null));
final RasterLayerRequest requestNE = new RasterLayerRequest(
new GeneralParameterValue[] { requestedBBox }, manager);
BoundingBox checkCropBBox = requestNE.spatialRequestHelper.getComputedBBox();
assertNotNull(checkCropBBox);
assertEquals(
"ReferencedEnvelope[1587997.8835 : 1612003.2265, 6162000.4515 : 6198002.1165]",
checkCropBBox.toString());
// set up the request (east-north version)
final ReferencedEnvelope requestBBoxEN = new ReferencedEnvelope(1583436.86902,
1617044.34782, 6154440.101350001, 6204842.43235, crs_EN);
requestedBBox.setValue(new GridGeometry2D(PixelInCell.CELL_CENTER, reader
.getOriginalGridToWorld(PixelInCell.CELL_CENTER), requestBBoxEN, null));
final RasterLayerRequest requestEN = new RasterLayerRequest(
new GeneralParameterValue[] { requestedBBox }, manager);
checkCropBBox = requestEN.spatialRequestHelper.getComputedBBox();
assertNotNull(checkCropBBox);
assertEquals(
"ReferencedEnvelope[1587997.8835 : 1612003.2265, 6162000.4515 : 6198002.1165]",
checkCropBBox.toString());
}
}