Application: gvSIG desktop

package org.gvsig.fmap;

import java.io.File;

import java.net.URL;

import junit.framework.Test;

import junit.framework.TestCase;

import junit.framework.TestSuite;

import org.cresques.cts.IProjection;

import org.gvsig.fmap.crs.CRSFactory;

import org.gvsig.fmap.drivers.formats.dbf.reading.TestDbf;

import org.gvsig.fmap.drivers.formats.dgn.reading.TestDgn;

import org.gvsig.fmap.drivers.iterators.DBFeatureIteratorTest;

import org.gvsig.fmap.drivers.iterators.FeatureIteratorTest;

import org.gvsig.fmap.drivers.iterators.PerformanceFeatureIteratorTest;

import org.gvsig.fmap.drivers.iterators.spatialindex.SpatialIndexTest;

import org.gvsig.fmap.mapcontext.MapContext;

import org.gvsig.fmap.mapcontext.ViewPort;

import org.gvsig.fmap.mapcontext.exceptions.LoadLayerException;

import org.gvsig.fmap.mapcontext.layers.FLayer;

import org.gvsig.fmap.mapcontext.layers.LayerFactory;

import org.gvsig.fmap.mapcontext.layers.LayersIteratorTest;

import org.gvsig.fmap.mapcontext.rendering.legend.TestAbstractIntervalLegend;

import org.gvsig.fmap.mapcontext.rendering.legend.TestIClassifiedLegend;

import org.gvsig.fmap.mapcontext.rendering.legend.TestILegend;

import org.gvsig.fmap.mapcontext.rendering.symbols.TestCartographicSupportForSymbol;

import org.gvsig.fmap.mapcontext.rendering.symbols.TestCartographySupport;

import org.gvsig.fmap.mapcontext.rendering.symbols.TestDrawFills;

import org.gvsig.fmap.mapcontext.rendering.symbols.TestDrawLines;

import org.gvsig.fmap.mapcontext.rendering.symbols.TestDrawMarkers;

import org.gvsig.fmap.mapcontext.rendering.symbols.TestISymbol;

import org.gvsig.fmap.mapcontext.rendering.symbols.TestMultiLayerSymbol;

import org.gvsig.fmap.mapcontrol.tools.AreaListenerTest;

import com.vividsolutions.jts.operation.overlay.SnappingOverlayOperationTest;

public class AllTests extends TestCase{

	/**

	 * The EPSG:4326 projection

	 */

	public static IProjection TEST_DEFAULT_PROJECTION =

		CRSFactory.getCRS("EPSG:4326");

	/**

	 * The EPSG:23030 projection

	 */

	public static IProjection TEST_DEFAULT_MERCATOR_PROJECTION =

		CRSFactory.getCRS("EPSG:23030");

	/**

	 * The EPSG:23029 projection

	 */

	public static IProjection test_newProjection =

		CRSFactory.getCRS("EPSG:23029");

	public static Test suite() {

		TestSuite suite = new TestSuite("Test for com.iver.cit.gvsig.fmap");

		//$JUnit-BEGIN$

		/* LayersIterator */

			suite.addTestSuite(LayersIteratorTest.class);

		/* Symbols (jaume) */

			// integration tests

				// ISymbol

				suite.addTestSuite(TestISymbol.class);

				suite.addTestSuite(TestMultiLayerSymbol.class);

				// CartographicSupport

				suite.addTestSuite(TestCartographicSupportForSymbol.class);

		//

		/*

		 * Feature iterators

		 * */

		suite.addTestSuite(FeatureIteratorTest.class);

		/*

		 * Other Tests present in FMap (cesar)

		 * Remove them from here and the src-test dir if they are not

		 * useful anymore.

		 */

		suite.addTestSuite(AreaListenerTest.class);

		suite.addTestSuite(DBFeatureIteratorTest.class);

		suite.addTestSuite(PerformanceFeatureIteratorTest.class);

		suite.addTestSuite(SnappingOverlayOperationTest.class);

		suite.addTestSuite(SpatialIndexTest.class);

		suite.addTestSuite(TestAbstractIntervalLegend.class);

		suite.addTestSuite(TestCartographicSupportForSymbol.class);

		suite.addTestSuite(TestDbf.class);

		suite.addTestSuite(TestDgn.class);

		suite.addTestSuite(TestDrawFills.class);

		suite.addTestSuite(TestDrawLines.class);

		suite.addTestSuite(TestDrawMarkers.class);

		suite.addTestSuite(TestIClassifiedLegend.class);

		suite.addTestSuite(TestILegend.class);

		suite.addTestSuite(TestMultiLayerSymbol.class);

		//$JUnit-END$

		return suite;

	}

//// jaume

//// PASTED FROM FeatureIteratorTest.java to be globally accessible

	static final String fwAndamiDriverPath = "../_fwAndami/gvSIG/extensiones/com.iver.cit.gvsig/drivers";

	private static File baseDataPath;

	private static File baseDriversPath;

	public static void setUpDrivers() {

		try {

			URL url = AllTests.class.getResource("testdata");

			if (url == null)

				throw new Exception("No se encuentra el directorio con datos de prueba");

			baseDataPath = new File(url.getFile());

			if (!baseDataPath.exists())

				throw new Exception("No se encuentra el directorio con datos de prueba");

			baseDriversPath = new File(fwAndamiDriverPath);

			if (!baseDriversPath.exists())

				throw new Exception("Can't find drivers path: " + fwAndamiDriverPath);

			LayerFactory.setDriversPath(baseDriversPath.getAbsolutePath());

			if (LayerFactory.getDM().getDriverNames().length < 1)

				throw new Exception("Can't find drivers in path: " + fwAndamiDriverPath);

		} catch (Exception e) {

			throw new RuntimeException(e);

		}

	}

	public static FLayer newLayer(String fileName,

			String driverName)

	throws LoadLayerException {

		File file = new File(baseDataPath, fileName);

		return LayerFactory.createLayer(fileName,

				driverName,

				file, TEST_DEFAULT_MERCATOR_PROJECTION);

	}

	public static MapContext newMapContext(IProjection projection) {

		ViewPort vp = new ViewPort(projection);

		return new MapContext(vp);

	}

//// end past

}

/* gvSIG. Sistema de Informaci?n Geogr?fica de la Generalitat Valenciana

 *

 * Copyright (C) 2004 IVER T.I. and Generalitat Valenciana.

 *

 * This program is free software; you can redistribute it and/or

 * modify it under the terms of the GNU General Public License

 * as published by the Free Software Foundation; either version 2

 * of the License, or (at your option) any later version.

 *

 * This program 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 General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 *

 * For more information, contact:

 *

 *  Generalitat Valenciana

 *   Conselleria d'Infraestructures i Transport

 *   Av. Blasco Ib??ez, 50

 *   46010 VALENCIA

 *   SPAIN

 *

 *      +34 963862235

 *   gvsig@gva.es

 *      www.gvsig.gva.es

 *

 *    or

 *

 *   IVER T.I. S.A

 *   Salamanca 50

 *   46005 Valencia

 *   Spain

 *

 *   +34 963163400

 *   dac@iver.es

 */

package org.gvsig.fmap.drivers.formats.dgn.reading;

import java.io.File;

import org.gvsig.fmap.drivers.formats.dgn.reading.DGNElemCore;

import org.gvsig.fmap.drivers.formats.dgn.reading.DGNReader;

import org.gvsig.fmap.drivers.formats.dgn.reading.DgnMemoryDriver;

import junit.framework.TestCase;

import com.hardcode.gdbms.driver.exceptions.InitializeDriverException;

import com.hardcode.gdbms.driver.exceptions.OpenDriverException;

import com.hardcode.gdbms.driver.exceptions.ReadDriverException;

public class TestDgn extends TestCase {

	DgnMemoryDriver driver = new DgnMemoryDriver();

	File file = new File("c:/LP01C.dgn");

	public static void main(String[] args) {

		junit.textui.TestRunner.run(TestDgn.class);

	}

	protected void setUp() throws Exception {

		super.setUp();

	}

	/*

	 * Test method for 'com.iver.cit.gvsig.fmap.drivers.dgn.DgnMemoryDriver.initialize()'

	 */

	public void testInitialize() {

		try {

			driver.open(file);

			driver.initialize();

			DGNReader reader = new DGNReader(file.getAbsolutePath());

				reader.DGNGotoElement(100);

				DGNElemCore elemento = reader.DGNReadElement();

				reader.DGNDumpElement(reader.getInfo(), elemento, null);

		}  catch (OpenDriverException e) {

			e.printStackTrace();

		} catch (InitializeDriverException e) {

			e.printStackTrace();

		} catch (ReadDriverException e) {

			e.printStackTrace();

		}

	}

}

/* gvSIG. Sistema de Informaci?n Geogr?fica de la Generalitat Valenciana

 *

 * Copyright (C) 2004 IVER T.I. and Generalitat Valenciana.

 *

 * This program is free software; you can redistribute it and/or

 * modify it under the terms of the GNU General Public License

 * as published by the Free Software Foundation; either version 2

 * of the License, or (at your option) any later version.

 *

 * This program 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 General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 *

 * For more information, contact:

 *

 *  Generalitat Valenciana

 *   Conselleria d'Infraestructures i Transport

 *   Av. Blasco Ib??ez, 50

 *   46010 VALENCIA

 *   SPAIN

 *

 *      +34 963862235

 *   gvsig@gva.es

 *      www.gvsig.gva.es

 *

 *    or

 *

 *   IVER T.I. S.A

 *   Salamanca 50

 *   46005 Valencia

 *   Spain

 *

 *   +34 963163400

 *   dac@iver.es

 */

package org.gvsig.fmap.drivers.formats.dbf.reading;

import java.io.File;

import org.gvsig.fmap.drivers.formats.dbf.reading.DBFDriver;

import junit.framework.TestCase;

import com.hardcode.gdbms.driver.exceptions.OpenDriverException;

public class TestDbf extends TestCase {

	DBFDriver driver = new DBFDriver();

	File file = new File("c:/prueba.dbf");

	public static void main(String[] args) {

		junit.textui.TestRunner.run(TestDbf.class);

	}

	protected void setUp() throws Exception {

		super.setUp();

	}

	/*

	 * Test method for 'com.iver.cit.gvsig.fmap.drivers.dbf.DBFDriver.initialize()'

	 */

	public void testInitialize() {

		try {

			driver.open(file);

			System.out.println(driver.getFields());//FieldValue(1,0));

		} catch (OpenDriverException e) {

			e.printStackTrace();

		}

	}

}

/*

 * Created on 17-abr-2007

 *

 * gvSIG. Sistema de Informaci?n Geogr?fica de la Generalitat Valenciana

 *

 * Copyright (C) 2004 IVER T.I. and Generalitat Valenciana.

 *

 * This program is free software; you can redistribute it and/or

 * modify it under the terms of the GNU General Public License

 * as published by the Free Software Foundation; either version 2

 * of the License, or (at your option) any later version.

 *

 * This program 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 General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 *

 * For more information, contact:

 *

 *  Generalitat Valenciana

 *   Conselleria d'Infraestructures i Transport

 *   Av. Blasco Ib??ez, 50

 *   46010 VALENCIA

 *   SPAIN

 *

 *      +34 963862235

 *   gvsig@gva.es

 *      www.gvsig.gva.es

 *

 *    or

 *

 *   IVER T.I. S.A

 *   Salamanca 50

 *   46005 Valencia

 *   Spain

 *

 *   +34 963163400

 *   dac@iver.es

 */

/* CVS MESSAGES:

*

* $Id$

* $Log$

* Revision 1.11  2007-09-20 08:08:29  jaume

* ReadExpansionFileException removed from this context

*

* Revision 1.10  2007/06/29 13:07:01  jaume

* +PictureLineSymbol

*

* Revision 1.9  2007/06/27 20:18:35  azabala

* *** empty log message ***

*

* Revision 1.8  2007/06/07 10:20:38  azabala

* includes closeIterator

*

* Revision 1.7  2007/06/07 09:31:42  azabala

* *** empty log message ***

*

* Revision 1.5  2007/05/29 19:11:03  azabala

* *** empty log message ***

*

* Revision 1.4  2007/05/23 16:54:29  azabala

* new test for bug 2510 phpcollab task (infinite loop for poly-valencia.shp)

*

* Revision 1.3  2007/05/08 08:48:27  jaume

* GRAPHIC TESTING FOR ISYMBOLS!!!

*

* Revision 1.2  2007/04/25 15:16:10  azabala

* tests for featureiterators with mysql database

*

* Revision 1.1  2007/04/19 18:12:56  azabala

* *** empty log message ***

*

*

*/

package org.gvsig.fmap.drivers.iterators;

import java.awt.Dimension;

import java.awt.Toolkit;

import java.awt.geom.Rectangle2D;

import java.awt.image.BufferedImage;

import java.io.File;

import java.net.URL;

import java.util.prefs.Preferences;

import junit.framework.TestCase;

import org.cresques.cts.IProjection;

import org.geotools.resources.geometry.XRectangle2D;

import org.gvsig.fmap.core.features.IFeature;

import org.gvsig.fmap.core.geometries.IGeometry;

import org.gvsig.fmap.core.geometries.utils.FConverter;

import org.gvsig.fmap.crs.CRSFactory;

import org.gvsig.fmap.drivers.reading.ConnectionJDBC;

import org.gvsig.fmap.drivers.reading.DBLayerDefinition;

import org.gvsig.fmap.drivers.reading.DefaultJDBCDriver;

import org.gvsig.fmap.drivers.reading.ICanReproject;

import org.gvsig.fmap.drivers.reading.IFeatureIterator;

import org.gvsig.fmap.drivers.reading.IVectorialJDBCDriver;

import org.gvsig.fmap.mapcontext.MapContext;

import org.gvsig.fmap.mapcontext.ViewPort;

import org.gvsig.fmap.mapcontext.exceptions.LoadLayerException;

import org.gvsig.fmap.mapcontext.layers.FLayer;

import org.gvsig.fmap.mapcontext.layers.LayerFactory;

import org.gvsig.fmap.mapcontext.layers.vectorial.FLyrVect;

import com.hardcode.gdbms.driver.exceptions.ReadDriverException;

import com.iver.utiles.swing.threads.Cancellable;

import com.vividsolutions.jts.geom.Geometry;

import com.vividsolutions.jts.geom.GeometryFactory;

import com.vividsolutions.jts.io.ParseException;

import com.vividsolutions.jts.io.WKTReader;

public class FeatureIteratorTest extends TestCase {

	static final String fwAndamiDriverPath = "../_fwAndami/gvSIG/extensiones/com.iver.cit.gvsig/drivers";

	private static File baseDataPath;

	private static File baseDriversPath;

	public static String SHP_DRIVER_NAME = "gvSIG shp driver";

	public  static String DXF_DRIVER_NAME = "gvSIG DXF Memory Driver";

	//TODO MOVER TODO LO ESTATICO A UNA CLASE AUXILIAR QUE NO SEA JUNIT

	static IProjection PROJECTION_DEFAULT =

		CRSFactory.getCRS("EPSG:23030");

	static IProjection newProjection =

		CRSFactory.getCRS("EPSG:23029");

	static{

		try {

			doSetup();

		} catch (Exception e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

	}

	private static void doSetup() throws Exception{

		URL url = FeatureIteratorTest.class.getResource("testdata");

		if (url == null)

			throw new Exception("No se encuentra el directorio con datos de prueba");

		baseDataPath = new File(url.getFile());

		if (!baseDataPath.exists())

			throw new Exception("No se encuentra el directorio con datos de prueba");

		baseDriversPath = new File(fwAndamiDriverPath);

		if (!baseDriversPath.exists())

			throw new Exception("Can't find drivers path: " + fwAndamiDriverPath);

		LayerFactory.setDriversPath(baseDriversPath.getAbsolutePath());

		if (LayerFactory.getDM().getDriverNames().length < 1)

			throw new Exception("Can't find drivers in path: " + fwAndamiDriverPath);

	}

	protected void setUp() throws Exception {

		super.setUp();

		doSetup();

	}

	public static  FLayer newLayer(String fileName,

									   String driverName)

								throws LoadLayerException {

		File file = new File(baseDataPath, fileName);

		return LayerFactory.createLayer(fileName,

										driverName,

										file, PROJECTION_DEFAULT);

	}

	public static FLyrVect newJdbcLayer(String layerName) throws LoadLayerException{

	        String dbURL = "jdbc:mysql://localhost:3306/datos";

	        String user = "root";

	        String pwd = "root";

	        String tableName = layerName;

	        //Connection conn;

	        try {

				Class.forName("com.mysql.jdbc.Driver");

			} catch (ClassNotFoundException e1) {

				throw new LoadLayerException(layerName, e1);

			}

			try {

				ConnectionJDBC conn = new ConnectionJDBC();

				conn.setDataConnection(dbURL, user, pwd);

		        String fidField = "gid";

		        String geomField = "geom";

		        String[] fields = new String[1];

		        fields[0] = "gid";

		        String whereClause = "";

		        IVectorialJDBCDriver driver = (IVectorialJDBCDriver)

		        	LayerFactory.getDM().getDriver("mySQL JDBC Driver");

		        String strEPSG = "23030";

		        DBLayerDefinition lyrDef = new DBLayerDefinition();

		        lyrDef.setName(layerName);

		        lyrDef.setTableName(tableName);

		        lyrDef.setWhereClause(whereClause);

		        lyrDef.setFieldGeometry(geomField);

		        lyrDef.setFieldNames(fields);

		        lyrDef.setFieldID(fidField);

		        lyrDef.setSRID_EPSG(strEPSG);

		        lyrDef.setConnection(conn);

		        if (driver instanceof ICanReproject)

		        {

		            ((ICanReproject)driver).setDestProjection(strEPSG);

		        }

		        ((DefaultJDBCDriver)driver).setData(conn, lyrDef);

		        IProjection proj = newProjection;//ahora, si no reproyecta la bbdd reproyecta el iterador

//		        if (driver instanceof ICanReproject)

//		        {

//		            proj = CRSFactory.getCRS("EPSG:" + ((ICanReproject)driver).getSourceProjection());

//		        }

		        return (FLyrVect) LayerFactory.createDBLayer(driver, layerName, proj);

			}catch(Exception e){

				e.printStackTrace();

			}

			return null;

	}

	//test para chequear la rapidez de los metodos de verificacion de candidatos en consultas

	//a partir de rectangulo

	public void test0(){

		String pol1 = "MULTIPOLYGON (((413434.195990259 4790622.699703139, 413441.7000601477 4790624.699360616, 413442.6978540559 4790620.700045662, 413435.2020303979 4790618.700388186, 413434.195990259 4790622.699703139, 413434.195990259 4790622.699703139)))";

		String pol2 = "MULTIPOLYGON (((411585.59723499016 4791781.20126231, 411575.39664767997 4791791.999412685, 411586.29816459515 4791802.297648691, 411596.4987519054 4791791.499498316, 411585.59723499016 4791781.20126231, 411585.59723499016 4791781.20126231)))";

		WKTReader reader = new WKTReader(new GeometryFactory());

		try {

			Geometry geo1 = reader.read(pol1);

			Geometry geo2 = reader.read(pol2);

			IGeometry igeo1 = FConverter.jts_to_igeometry(geo1);

			IGeometry igeo2 = FConverter.jts_to_igeometry(geo2);

			double xmin = 410000;

			double xmax = 415000;

			double ymin = 4790000;

			double ymax = 4793000;

			Rectangle2D rect = new Rectangle2D.Double(xmin, ymin, (xmax-xmin), (ymax-ymin));

			//PrecciseSpatialCheck

			 long t0 = System.currentTimeMillis();

			 int NUM_ITERATIONS = 50000;

			 for(int i = 0; i < NUM_ITERATIONS;i++){

				 igeo1.fastIntersects(rect.getMinX(),

						rect.getMinY(), rect.getWidth(), rect.getHeight());

				 igeo2.fastIntersects(rect.getMinX(),

							rect.getMinY(), rect.getWidth(), rect.getHeight());

			 }

			 long t1 = System.currentTimeMillis();

			 System.out.println((t1-t0)+" con precissespatialcheck");

			//FastSpatialCheck

			 Rectangle2D b1 = igeo1.getBounds2D();

			 Rectangle2D b2 = igeo2.getBounds2D();

			 long t2 = System.currentTimeMillis();

			 for(int i = 0; i < NUM_ITERATIONS;i++){

				 XRectangle2D.intersectInclusive(rect, b1);

				 XRectangle2D.intersectInclusive(rect, b2);

			 }

			 long t3 = System.currentTimeMillis();

			 System.out.println((t3-t2)+" con fastspatialcheck");

			 assertTrue((t3-t2) < (t1 - t0));

			//Spatial index check

		} catch (ParseException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

	}

	//pruebas de iteracion para shp (vectorialfileadapter) y dxf (vectorialadapter)

	public void test1() {

		try {

			//pruebas de reproyeccion y seleccion de numero de campos

			FLyrVect lyr = (FLyrVect) newLayer("Cantabria.shp", SHP_DRIVER_NAME);

			lyr.setAvailable(true);

			//iteration selecting two fields

			IFeatureIterator iterator = lyr.getSource().getFeatureIterator(new String[]{"ID", "AREA"}, null);

			IFeature feature = iterator.next();

			assertTrue(feature != null);

			//iteration with reprojection

			iterator = lyr.getSource().getFeatureIterator((String[])null, newProjection);

			IFeature feature2 = iterator.next();

			//test of field restriction

			assertTrue(feature.getAttributes().length == 2);

			//si pasamos null, no devuelve ning?n campo

			assertTrue(feature2.getAttributes().length == 0);

			assertTrue(!feature.getGeometry().toJTSGeometry().equals(feature2.getGeometry().toJTSGeometry()));

			iterator.closeIterator();

			//pruebas de iteracion espacial

			FLyrVect lyr2 = (FLyrVect) newLayer("Edificaciones.shp", SHP_DRIVER_NAME);

			lyr2.setAvailable(true);

			double xmin = 410000;

			double xmax = 415000;

			double ymin = 4790000;

			double ymax = 4793000;

			Rectangle2D rect = new Rectangle2D.Double(xmin, ymin, (xmax-xmin), (ymax-ymin));

			//fast iteration

			long t0 = System.currentTimeMillis();

			iterator = lyr2.getSource().getFeatureIterator(rect, null, newProjection, true);

			int numFeatures = 0;

			while(iterator.hasNext()){

				feature = iterator.next();

				numFeatures++;

			}

			long t1 = System.currentTimeMillis();

			//vemos si hay diferencia entre el fast = true y el fast = false

			iterator = lyr2.getSource().getFeatureIterator(rect, null, newProjection, false);

			int numFeatures3 = 0;

			while(iterator.hasNext()){

				feature = iterator.next();

				numFeatures3++;

			}

			long t2 = System.currentTimeMillis();

			assertTrue(numFeatures3 <= numFeatures);

			iterator.closeIterator();

			System.out.println((t1-t0)+" en la iteracion rapida");

			System.out.println("Recuperados "+numFeatures);

			System.out.println((t2-t1)+" en la iteracion lenta");

			System.out.println("Recuperados "+numFeatures3);

			int numFeatures2 = 0;

			iterator = lyr2.getSource().getFeatureIterator();

			while(iterator.hasNext()){

				feature = iterator.next();

				numFeatures2++;

			}

			assertTrue(numFeatures2 > numFeatures);

			iterator.closeIterator();

			//pruebas de iteracion en base a criterios alfanumericos

			final String expr = "select * from " + lyr.getRecordset().getName() + " where area < 700;";

			iterator = lyr.getSource().getFeatureIterator(expr, null);

			numFeatures2 = 0;

			while(iterator.hasNext()){

				feature = iterator.next();

				numFeatures2++;

			}

			iterator.closeIterator();

			FLayer lyr53245 = LayerFactory.createLayer("/home/jaume/callestodo.shp",

					SHP_DRIVER_NAME,

					new File("/home/jaume/callestodo.shp"), PROJECTION_DEFAULT);

			iterator = ((FLyrVect) lyr53245).getSource().getFeatureIterator("select * from " + ((FLyrVect) lyr53245).getRecordset().getName() +" where (TIPCALLE = 'C' or  TIPCALLE='AV' or TIPCALLE='PL');", null);

			while(iterator.hasNext()){

				feature = iterator.next();

				numFeatures++;

			}

			iterator.closeIterator();

			iterator = lyr.getSource().getFeatureIterator();

			while(iterator.hasNext()){

				feature = iterator.next();

				numFeatures++;

			}

			assertTrue(numFeatures2 > 0);

			assertTrue(numFeatures > numFeatures2);

			iterator.closeIterator();

			//pruebas con el driver dxf en vez del shp

			FLyrVect lyr3 = (FLyrVect) newLayer("es1120003_2_02.dxf", DXF_DRIVER_NAME);

			lyr3.setAvailable(true);

			iterator = lyr3.getSource().getFeatureIterator("select * from "+lyr3.getRecordset().getName()+" where ID > 17;", null);

			numFeatures2 = 0;

			while(iterator.hasNext()){

				feature = iterator.next();

				numFeatures2++;

			}

			iterator.closeIterator();

			assertTrue(numFeatures2 > 0);

			//queda por probar los iteradores contra base de datos

			//en base a criterios alfanumericos y espaciales

		} catch (LoadLayerException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		} catch (ReadDriverException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		} 

	}

//	public void test3(){

//		try {

//			FLyrVect layer = (FLyrVect) newLayer("poly-valencia.shp", SHP_DRIVER_NAME);

//			IFeatureIterator iterator = layer.getSource().getFeatureIterator();

//			int numFeatures = 0;

//			while(iterator.hasNext()){

//				iterator.next();

//				numFeatures++;

//			}

//			assert(layer.getSource().getShapeCount() == numFeatures);

//

//		} catch (LoadLayerException e) {

//			// TODO Auto-generated catch block

//			e.printStackTrace();

//		} catch (ReadDriverException e) {

//			// TODO Auto-generated catch block

//			e.printStackTrace();

//		} catch (ExpansionFileReadException e) {

//			// TODO Auto-generated catch block

//			e.printStackTrace();

//		}

//

//	}

	//test of featureIterator with layer poly-valencia using a spatial filter.

	//jaume found a bug with an infinite loop iteration

	public void test4(){

		try {

			FLyrVect layer = (FLyrVect) newLayer("poly-valencia.shp", SHP_DRIVER_NAME);

			Rectangle2D extent = layer.getFullExtent();

			String[] fields = layer.getRecordset().getFieldNames();

			IFeatureIterator iterator = layer.

							getSource().

							getFeatureIterator(extent,fields, null, true);

			int numFeatures = 0;

			while(iterator.hasNext()){

				iterator.next();

				numFeatures++;

			}

			iterator.closeIterator();

			assert(layer.getSource().getShapeCount() == numFeatures);

		} catch (LoadLayerException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		} catch (ReadDriverException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

	}

	public void test5(){

		try {

			int WIDTH = 300;

			int HEIGHT = 200;

			FLyrVect layer = (FLyrVect) newLayer("municipios.shp", SHP_DRIVER_NAME);

			ViewPort vp = new ViewPort(PROJECTION_DEFAULT);

			vp.setImageSize(new Dimension(WIDTH, HEIGHT));

			vp.setExtent(layer.getFullExtent());

			MapContext mapa = new MapContext(vp);

			mapa.getLayers().addLayer(layer);

			BufferedImage img = new BufferedImage(WIDTH, HEIGHT,

					BufferedImage.TYPE_INT_ARGB);

			long t0 = System.currentTimeMillis();

			mapa.draw(img, img.createGraphics(),mapa.getScaleView());

			long t1 = System.currentTimeMillis();

			System.out.println((t1-t0)+" en dibujar con MapContext.draw");

			long t2 = System.currentTimeMillis();

			layer.draw(img, img.createGraphics(),vp, new Cancellable(){

				public boolean isCanceled() {

					return false;

				}

				public void setCanceled(boolean canceled) {

				}}, getScaleView(vp));

			long t3 = System.currentTimeMillis();

			System.out.println((t3-t2)+" en dibujar con layer._draw");

		} catch (LoadLayerException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		} catch (ReadDriverException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		} 

	}

	//TODO Mover spatialindextest a esta clase

	public void test7(){

	}

	private long getScaleView(ViewPort viewPort) {

		double[] CHANGE = { 100000, 100, 1, 0.1, 160934.4,

				91.44, 30.48, 2.54, 1/8.983152841195214E-4 };

		Preferences prefsResolution = Preferences.userRoot().node( "gvsig.configuration.screen" );

		Toolkit kit = Toolkit.getDefaultToolkit();

		double dpi = prefsResolution.getInt("dpi",kit.getScreenResolution());

		IProjection proj = viewPort.getProjection();

		if (viewPort.getImageSize() == null)

			return -1;

		if (viewPort.getAdjustedExtent() == null) {

			return 0;

		}

		if (proj == null) {

			double w = ((viewPort.getImageSize().getWidth() / dpi) * 2.54);

			return (long) (viewPort.getAdjustedExtent().getWidth() / w * CHANGE[viewPort.getMapUnits()]);

		}

		return Math.round(proj.getScale(viewPort.getAdjustedExtent().getMinX(),

				viewPort.getAdjustedExtent().getMaxX(), viewPort.getImageSize()

						.getWidth(), dpi));

	}

}

package org.gvsig.fmap.drivers.iterators.spatialindex;

import java.awt.geom.Point2D;

import java.awt.geom.Rectangle2D;

import java.io.File;

import java.net.URL;

import java.util.List;

import junit.framework.TestCase;

import org.cresques.cts.IProjection;

import org.gvsig.fmap.core.geometries.IGeometry;

import org.gvsig.fmap.crs.CRSFactory;

import org.gvsig.fmap.drivers.exceptions.ExpansionFileReadException;

import org.gvsig.fmap.drivers.exceptions.VisitorException;

import org.gvsig.fmap.drivers.iterators.spatialindex.PersistentRTreeJsi;

import org.gvsig.fmap.drivers.iterators.spatialindex.QuadtreeGt2;

import org.gvsig.fmap.drivers.iterators.spatialindex.SpatialIndexException;

import org.gvsig.fmap.mapcontext.exceptions.LoadLayerException;

import org.gvsig.fmap.mapcontext.layers.FBitSet;

import org.gvsig.fmap.mapcontext.layers.FLayer;

import org.gvsig.fmap.mapcontext.layers.LayerFactory;

import org.gvsig.fmap.mapcontext.layers.vectorial.FLyrVect;

import com.hardcode.gdbms.driver.exceptions.ReadDriverException;

public class SpatialIndexTest extends TestCase {

	// TODO MOVER TODO LO ESTATICO A UNA CLASE AUXILIAR QUE NO SEA JUNIT

	String fwAndamiDriverPath = "../_fwAndami/gvSIG/extensiones/com.iver.cit.gvsig/drivers";

	File baseDataPath;

	File baseDriversPath;

	String SHP_DRIVER_NAME = "gvSIG shp driver";

	IProjection PROJECTION_DEFAULT = CRSFactory.getCRS("EPSG:23030");

	FLyrVect cantabria;

	FLyrVect edificaciones;

	public SpatialIndexTest(String name) {

		super(name);

	}

	protected void setUp() throws Exception {

		super.setUp();

		URL url = SpatialIndexTest.class.getResource("testdata");

		if (url == null)

			throw new Exception(

					"No se encuentra el directorio con datos de prueba");

		baseDataPath = new File(url.getFile());

		if (!baseDataPath.exists())

			throw new Exception(

					"No se encuentra el directorio con datos de prueba");

		baseDriversPath = new File(fwAndamiDriverPath);

		if (!baseDriversPath.exists())

			throw new Exception("Can't find drivers path: "

					+ fwAndamiDriverPath);

		LayerFactory.setDriversPath(baseDriversPath.getAbsolutePath());

		if (LayerFactory.getDM().getDriverNames().length < 1)

			throw new Exception("Can't find drivers in path: "

					+ fwAndamiDriverPath);

		cantabria = (FLyrVect) newLayer("Cantabria.shp", SHP_DRIVER_NAME);

		edificaciones = (FLyrVect) newLayer("Edificaciones.shp", SHP_DRIVER_NAME);

	}

	public FLayer newLayer(String fileName, String driverName)

			throws LoadLayerException {

		File file = new File(baseDataPath, fileName);

		return LayerFactory.createLayer(fileName, driverName, file,

				PROJECTION_DEFAULT);

	}

	protected void tearDown() throws Exception {

		super.tearDown();

		PROJECTION_DEFAULT = null;

	}

	public void testSpatialIndexFullExtent() throws ReadDriverException,

			ExpansionFileReadException, VisitorException, LoadLayerException {

		Rectangle2D rect2D = cantabria.getFullExtent();

		FBitSet bitset = cantabria.queryByRect(rect2D);

		assertTrue(bitset.cardinality() != 0);

		double x = rect2D.getCenterX();

		double y = rect2D.getCenterY();

		rect2D = new Rectangle2D.Double(x, y, 100d, 100d);

		bitset = cantabria.queryByRect(rect2D);

		assertTrue(bitset.cardinality() != 0);

	}

	public void testPersistentRtreeJsi() throws LoadLayerException,

			ReadDriverException, ExpansionFileReadException,

			SpatialIndexException {

		try {

			PersistentRTreeJsi rtree = new PersistentRTreeJsi(baseDataPath

					.getAbsolutePath()

					+ "/Cantabria", true);

			int numShapes = cantabria.getSource().getShapeCount();

			for (int i = 0; i < numShapes; i++) {

				IGeometry g = cantabria.getSource().getShape(i);

				if (g == null)

					System.out.println("geometria a null");

				rtree.insert(g.getBounds2D(), i);

			}

			rtree.flush();

			rtree = null;

			rtree = new PersistentRTreeJsi(baseDataPath.getAbsolutePath()

					+ "/Cantabria", false);

			Rectangle2D rect2D = cantabria.getFullExtent();

			List results = rtree.query(rect2D);

			assertTrue(results.size() != 0);

		} catch (Exception e) {

			e.printStackTrace();

		}

	}

	public void testPerformanceRtreeJsiQuadtreeGt2() throws LoadLayerException, ExpansionFileReadException, SpatialIndexException, ReadDriverException {

		String fileName = baseDataPath.getAbsolutePath()

							+ "/Edificaciones.shp";

		long t0 = System.currentTimeMillis();

		QuadtreeGt2	spatialIndex = new QuadtreeGt2(fileName,

														"NM", 

										edificaciones.getSource().getFullExtent(),

										edificaciones.getSource().getShapeCount(), 

										true);

		int numShapes = edificaciones.getSource().getShapeCount();

		for (int i = 0; i < numShapes; i++) {

			IGeometry g = edificaciones.getSource().getShape(i);

			spatialIndex.insert(g.getBounds2D(), i);

		}

		spatialIndex.flush();

		long t1 = System.currentTimeMillis();

		System.out.println((t1-t0)+ " en generar Edificaciones.qix");

		t0 = System.currentTimeMillis();

		PersistentRTreeJsi rtree = new PersistentRTreeJsi(baseDataPath

				.getAbsolutePath()

				+ "/Edificaciones", true);

		numShapes = edificaciones.getSource().getShapeCount();

		for (int i = 0; i < numShapes; i++) {

			IGeometry g = edificaciones.getSource().getShape(i);

			rtree.insert(g.getBounds2D(), i);

		}

		t1 = System.currentTimeMillis();

		System.out.println((t1-t0)+" en generar Edificaciones.rix");

		rtree.flush();

		t0 = System.currentTimeMillis();

		spatialIndex = new QuadtreeGt2(fileName,

				"NM", 

				edificaciones.getSource().getFullExtent(),

				edificaciones.getSource().getShapeCount(), 

				false);

		t1 = System.currentTimeMillis();

		System.out.println((t1-t0)+" en cargar Edificaciones.qix");

		t0 = System.currentTimeMillis();

		rtree = new PersistentRTreeJsi(baseDataPath

				.getAbsolutePath()

				+ "/Edificaciones", false);

		t1 = System.currentTimeMillis();

		System.out.println((t1-t0)+" en cargar Edificaciones.rix");

		Rectangle2D rect2D = edificaciones.getFullExtent();

		t0 = System.currentTimeMillis();

		List results = rtree.query(rect2D);

		int nb1 = results.size();

		t1 = System.currentTimeMillis();

		System.out.println((t1-t0)+" en resolver query con RIX");

		System.out.println("encontrados "+nb1);

		t0 = System.currentTimeMillis();

		results = spatialIndex.query(rect2D);

		int nb2 = results.size();

		t1 = System.currentTimeMillis();

		System.out.println((t1-t0)+" en resolver query con QIX");

		System.out.println("encontrados "+nb2);

	}

	public void testNearestRectanglePoint() throws SpatialIndexException, ExpansionFileReadException, ReadDriverException{

		long t0 = System.currentTimeMillis();

		PersistentRTreeJsi rtree = new PersistentRTreeJsi(baseDataPath

				.getAbsolutePath()

				+ "/Edificaciones", false);

		long t1 = System.currentTimeMillis();

		System.out.println((t1-t0)+" en cargar Edificaciones.rix");

		Rectangle2D edifFullExtent = edificaciones.getFullExtent();

		double xc = edifFullExtent.getCenterX();

		double yc = edifFullExtent.getCenterY();

		double width = 200;

		int numberOfNearest = 20;

		Point2D p = new Point2D.Double(xc, yc);

		Rectangle2D query = new Rectangle2D.Double(xc, yc, width, width);

		t0 = System.currentTimeMillis();

		List solution = rtree.findNNearest(numberOfNearest, p);

		t1 = System.currentTimeMillis();

		System.out.println((t1-t0)+" en obtener los 20 mas cercanos al punto");

		assert(solution.size() > 0);

		t0 = System.currentTimeMillis();

		solution = rtree.findNNearest(numberOfNearest, query);

		t1 = System.currentTimeMillis();

		System.out.println((t1-t0)+" en obtener los 20 mas cercanos al rectangulo");

		assert(solution.size() > 0);

	}

}

PROJCS["European 1950 Based UTM, Zone 30 North, Meter",GEOGCS["ERP50-W",DATUM["D_ERP50-W",SPHEROID["International_1924",6378388,297.0000000000601]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-3],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["Meter",1]]

/*

 * Created on 28-may-2007

 *

 * gvSIG. Sistema de Informaci?n Geogr?fica de la Generalitat Valenciana

 *

 * Copyright (C) 2004 IVER T.I. and Generalitat Valenciana.

 *

 * This program is free software; you can redistribute it and/or

 * modify it under the terms of the GNU General Public License

 * as published by the Free Software Foundation; either version 2

 * of the License, or (at your option) any later version.

 *

 * This program 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 General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 *

 * For more information, contact:

 *

 *  Generalitat Valenciana

 *   Conselleria d'Infraestructures i Transport

 *   Av. Blasco Ib??ez, 50

 *   46010 VALENCIA

 *   SPAIN

 *

 *      +34 963862235

 *   gvsig@gva.es

 *      www.gvsig.gva.es

 *

 *    or

 *

 *   IVER T.I. S.A

 *   Salamanca 50

 *   46005 Valencia

 *   Spain

 *

 *   +34 963163400

 *   dac@iver.es

 */

/* CVS MESSAGES:

*

* $Id$

* $Log$

* Revision 1.6  2007-09-20 08:08:29  jaume

* ReadExpansionFileException removed from this context

*

* Revision 1.5  2007/06/29 13:07:01  jaume

* +PictureLineSymbol

*

* Revision 1.4  2007/06/07 10:20:38  azabala

* includes closeIterator

*

* Revision 1.3  2007/06/07 09:31:42  azabala

* *** empty log message ***

*

* Revision 1.2  2007/05/30 20:12:41  azabala

* fastIteration = true optimized.

*

* Revision 1.1  2007/05/29 19:11:03  azabala

* *** empty log message ***

*

*

*/

package org.gvsig.fmap.drivers.iterators;

import java.awt.geom.Rectangle2D;

import junit.framework.TestCase;

import org.cresques.cts.ICoordTrans;

import org.gvsig.fmap.core.features.IFeature;

import org.gvsig.fmap.drivers.exceptions.ExpansionFileReadException;

import org.gvsig.fmap.drivers.iterators.SpatialQueryFeatureIterator;

import org.gvsig.fmap.drivers.iterators.spatialindex.ISpatialIndex;

import org.gvsig.fmap.drivers.reading.IFeatureIterator;

import org.gvsig.fmap.mapcontext.exceptions.LoadLayerException;

import org.gvsig.fmap.mapcontext.layers.vectorial.FLyrVect;

import com.hardcode.gdbms.driver.exceptions.ReadDriverException;

/**

 * Tests to probe feature iteration methods.

 *

 * These test are not functional-test (performance).

 * @author azabala

 *

 */

public class PerformanceFeatureIteratorTest extends TestCase {

   static FLyrVect lyr;

	static{

			try {

				lyr = (FLyrVect) FeatureIteratorTest.newLayer("ejesc.shp",

										FeatureIteratorTest.SHP_DRIVER_NAME);

				lyr.setAvailable(true);

			} catch (LoadLayerException e) {

				// TODO Auto-generated catch block

				e.printStackTrace();

			}

	}

	public void test1() {

		try {

			//pruebas de iteracion espacial

			System.out.println("TEST 1: ESPACIAL CON FULL EXTENT Y REPROYECCI?N");

			Rectangle2D rect = lyr.getFullExtent();

			IFeature feature = null;

			//fast iteration

			long t0 = System.currentTimeMillis();

			ISpatialIndex spatialIndex = lyr.getSource().getSpatialIndex();

			lyr.getSource().setSpatialIndex(null);

			//Sin indice espacial, rapida

			//si pedimos reproyeccion, el rectangulo de consulta debe estar en la proyeccion

			//de destino

			ICoordTrans trans = FeatureIteratorTest.PROJECTION_DEFAULT.getCT(FeatureIteratorTest.newProjection);

			rect = trans.convert(rect);

			IFeatureIterator iterator = lyr.getSource().getFeatureIterator(rect,

															null,

															FeatureIteratorTest.newProjection,

															true);

			int numFeatures = 0;

			while(iterator.hasNext()){

				feature = iterator.next();

				numFeatures++;

			}

			long t1 = System.currentTimeMillis();

			iterator.closeIterator();

			//sin indice espacial, lenta

			iterator = lyr.getSource().getFeatureIterator(rect,

															null,

								FeatureIteratorTest.newProjection,

												false);

			int numFeatures2 = 0;

			while(iterator.hasNext()){

				feature = iterator.next();

				numFeatures2++;

			}

			long t2 = System.currentTimeMillis();

			iterator.closeIterator();

			lyr.getSource().setSpatialIndex(spatialIndex);

			long t3 = System.currentTimeMillis();

			//con indice espacial rapida

			iterator = lyr.getSource().getFeatureIterator(rect,

															null,

								FeatureIteratorTest.newProjection,

												true);

			int numFeatures3 = 0;

			while(iterator.hasNext()){

				feature = iterator.next();

				numFeatures3++;

			}

			long t4 = System.currentTimeMillis();

			iterator.closeIterator();

			//con indice espacial lenta

			iterator = lyr.getSource().getFeatureIterator(rect,

					null,

					FeatureIteratorTest.newProjection,

					false);

			int numFeatures4 = 0;

			while(iterator.hasNext()){

			feature = iterator.next();

			numFeatures4++;

			}

			long t5 = System.currentTimeMillis();

			iterator.closeIterator();

			System.out.println((t1-t0)+" en la iteracion rapida sin indice espacial");

			System.out.println("Recuperados "+numFeatures);

			System.out.println((t4-t3)+" en la iteracion rapida con indice espacial");

			System.out.println("Recuperados "+numFeatures3);

			System.out.println((t2-t1)+" en la iteracion lenta sin indice espacial");

			System.out.println("Recuperados "+numFeatures2);

			System.out.println((t5-t4)+" en la iteracion lenta con indice espacial");

			System.out.println("Recuperados "+numFeatures4);

		} catch (ReadDriverException e) {

			// TODO Auto-generated catch block