/*
* GeoTools - The Open Source Java GIS Toolkit
* http://geotools.org
*
* (C) 2002-2008, 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.data.store;
import java.util.Iterator;
import org.geotools.data.crs.ReprojectFeatureResults;
import org.geotools.geometry.jts.GeometryCoordinateSequenceTransformer;
import org.geotools.referencing.CRS;
import org.geotools.referencing.ReferencingFactoryFinder;
import org.opengis.feature.simple.SimpleFeature;
import org.opengis.referencing.crs.CoordinateReferenceSystem;
import org.opengis.referencing.operation.MathTransform2D;
import com.vividsolutions.jts.geom.LineString;
import com.vividsolutions.jts.geom.Point;
public class ReprojectFeatureResultsTest extends
FeatureCollectionWrapperTestSupport {
CoordinateReferenceSystem target;
GeometryCoordinateSequenceTransformer transformer;
protected void setUp() throws Exception {
super.setUp();
target = CRS.parseWKT(
"PROJCS[\"BC_Albers\",GEOGCS[\"GCS_North_American_1983\",DATUM[\"North_American_Datum_1983\",SPHEROID[\"GRS_1980\",6378137,298.257222101],TOWGS84[0,0,0]],PRIMEM[\"Greenwich\",0],UNIT[\"Degree\",0.017453292519943295]],PROJECTION[\"Albers_Conic_Equal_Area\"],PARAMETER[\"False_Easting\",1000000],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",-126],PARAMETER[\"Standard_Parallel_1\",50],PARAMETER[\"Standard_Parallel_2\",58.5],PARAMETER[\"Latitude_Of_Origin\",45],UNIT[\"Meter\",1],AUTHORITY[\"EPSG\",\"42102\"]]"
);
MathTransform2D tx = (MathTransform2D) ReferencingFactoryFinder.getCoordinateOperationFactory(null)
.createOperation(crs,target).getMathTransform();
transformer = new GeometryCoordinateSequenceTransformer();
transformer.setMathTransform( tx );
transformer.setCoordinateReferenceSystem( target );
}
public void testNormal() throws Exception {
Iterator reproject = new ReprojectFeatureResults( delegate, target ).iterator();
Iterator reader = delegate.iterator();
while( reader.hasNext() ) {
SimpleFeature normal = (SimpleFeature) reader.next();
SimpleFeature reprojected = (SimpleFeature) reproject.next();
Point p1 = (Point) normal.getAttribute( "defaultGeom" );
Point p2 = (Point) reprojected.getAttribute( "defaultGeom" );
if(p1 != null) {
assertEquals( crs, p1.getUserData() );
assertEquals( target, p2.getUserData() );
p1 = (Point) transformer.transform( p1 );
assertTrue( p1.equals( p2 ) );
} else {
assertNull(p2);
}
LineString l1 = (LineString) normal.getAttribute( "otherGeom" );
LineString l2 = (LineString) reprojected.getAttribute( "otherGeom" );
if(l1 != null) {
l1 = (LineString) transformer.transform( l1 );
assertTrue( l1.equals( l2 ) );
} else {
assertNull(l2);
}
}
}
public void testBounds() throws Exception {
ReprojectFeatureResults rfr = new ReprojectFeatureResults( delegate, target );
rfr.getBounds();
}
}