/*
* GeoTools - The Open Source Java GIS Toolkit
* http://geotools.org
*
* (C) 2010, 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.grid;
import java.util.Map;
import com.vividsolutions.jts.geom.Polygon;
import org.geotools.data.simple.SimpleFeatureCollection;
import org.geotools.data.simple.SimpleFeatureIterator;
import org.geotools.data.simple.SimpleFeatureSource;
import org.geotools.geometry.jts.ReferencedEnvelope;
import org.geotools.referencing.CRS;
import org.geotools.referencing.crs.DefaultGeographicCRS;
import org.opengis.feature.simple.SimpleFeature;
import org.opengis.referencing.FactoryException;
import org.opengis.referencing.crs.CoordinateReferenceSystem;
import org.junit.Test;
import static org.junit.Assert.*;
/**
* Unit tests for the Grids class.
*
* @author mbedward
* @since 2.7
*
* @source $URL$
* @version $Id$
*/
public class GridsSquareTest extends TestBase {
private final ReferencedEnvelope bounds = new ReferencedEnvelope(0, 90, 0, 100, null);
private final double sideLen = 9.0;
private final int expectedRows = (int) (bounds.getHeight() / sideLen);
private final int expectedCols = (int) (bounds.getWidth() / sideLen);
private final int expectedNumElements = expectedRows * expectedCols;
private final ReferencedEnvelope expectedBounds = new ReferencedEnvelope(0, sideLen * expectedCols, 0, sideLen * expectedRows, null);
@Test
public void createGrid() throws Exception {
SimpleFeatureSource gridSource = Grids.createSquareGrid(bounds, sideLen);
assertGridSizeAndIds(gridSource);
SimpleFeatureIterator iter = gridSource.getFeatures().features();
try {
Polygon poly = (Polygon) iter.next().getAttribute("element");
assertEquals(5, poly.getCoordinates().length);
} finally {
iter.close();
}
}
@Test
public void createDensifiedGrid() throws Exception {
final int vertexDensity = 10;
SimpleFeatureSource gridSource = Grids.createSquareGrid(bounds, sideLen, sideLen / vertexDensity);
assertGridSizeAndIds(gridSource);
SimpleFeatureIterator iter = gridSource.getFeatures().features();
try {
Polygon poly = (Polygon) iter.next().getAttribute("element");
assertTrue(poly.getCoordinates().length - 1 >= 4 * vertexDensity);
} finally {
iter.close();
}
}
@Test(expected=IllegalArgumentException.class)
public void createGrid_InvalidBounds() {
Grids.createSquareGrid(ReferencedEnvelope.EVERYTHING, sideLen);
}
@Test(expected=IllegalArgumentException.class)
public void createGrid_NullBounds() {
Grids.createSquareGrid(null, sideLen);
}
@Test(expected=IllegalArgumentException.class)
public void createGrid_BadSideLength() {
Grids.createSquareGrid(bounds, 0);
}
@Test(expected=IllegalArgumentException.class)
public void incompatibleCRS() {
CoordinateReferenceSystem boundsCRS = DefaultGeographicCRS.WGS84;
CoordinateReferenceSystem builderCRS;
try {
builderCRS = CRS.parseWKT(getSydneyWKT());
} catch (FactoryException ex) {
throw new IllegalStateException("Error in test code");
}
GridFeatureBuilder builder = new GridFeatureBuilder(createFeatureType(builderCRS)) {
@Override
public void setAttributes(GridElement el, Map<String, Object> attributes) {
throw new UnsupportedOperationException("Not supported yet.");
}
};
Grids.createSquareGrid(new ReferencedEnvelope(150, 151, -33, -34, boundsCRS), sideLen, sideLen, builder);
}
private void assertGridSizeAndIds(SimpleFeatureSource gridSource) throws Exception {
SimpleFeatureCollection grid = gridSource.getFeatures();
assertEquals(expectedNumElements, grid.size());
assertTrue(expectedBounds.boundsEquals2D(grid.getBounds(), TOL));
boolean[] flag = new boolean[expectedNumElements + 1];
int count = 0;
SimpleFeatureIterator iter = grid.features();
try {
while (iter.hasNext()) {
SimpleFeature f = iter.next();
int id = (Integer) f.getAttribute("id");
assertFalse(id == 0);
assertFalse(flag[id]);
flag[id] = true;
count++ ;
}
} finally {
iter.close();
}
assertEquals(expectedNumElements, count);
}
}