/*
 *    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.equalsExact( 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.equalsExact( l2 ) );
			} else {
			    assertNull(l2);
			}
		}
		
	}
	
	public void testBounds() throws Exception {
	    ReprojectFeatureResults rfr = new ReprojectFeatureResults( delegate, target );
	    rfr.getBounds();
	}
}