Application: gvSIG desktop

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

 *

 * Copyright (C) 2006 Instituto de Desarrollo Regional 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

 *

 *   Instituto de Desarrollo Regional (Universidad de Castilla La-Mancha)

 *   Campus Universitario s/n

 *   02071 Alabacete

 *   Spain

 *

 *   +34 967 599 200

 */

package es.idr.teledeteccion.connection;

import java.io.BufferedReader;

import java.io.File;

import java.io.IOException;

import java.io.InputStreamReader;

import java.sql.Connection;

import java.sql.ResultSet;

import java.sql.SQLException;

import java.sql.Statement;

import java.util.logging.Logger;

import org.geotools.referencing.factory.epsg.HSQLDataSource;

import org.geotools.referencing.factory.iau2000.FactoryUsingHSQL;

import org.hsqldb.jdbc.jdbcDataSource;

/**

 * Clase para la conexi?n con la base de datos de hsqldb.

 * Establece el driver necesario, as? como la cadena de

 * conexi?n a la base de datos de la EPSG y la IAU2000

 * 

 * @author Jos? Luis G?mez Mart?nez (jolugomar@gmail.com)

 *

 */

public class EpsgConnection extends jdbcDataSource {

	Connection connect;

	public EpsgConnection() {	

	/*	try {

			Class.forName("org.hsqldb.jdbcDriver");

		} catch (ClassNotFoundException e) {

			e.printStackTrace();

		}*/

	}

	public void setConnectionEPSG() {

		HSQLDataSource ds = new HSQLDataSource();

		try {

			connect = ds.getConnection();

		} catch (SQLException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

	}

	/**

	 * Establece la conexi?n con la base de datos de la IAU2000

	 *

	 */

	public void setConnectionIAU2000() {

		setDatabase("jdbc:hsqldb:file:gvSIG/extensiones/org.gvsig.crs.extension/iau2000");

		setUser("sa");

		try {

			connect = super.getConnection();

		} catch (SQLException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

	}

	/**

	 * Establece la conexi?n con la base de datos de ESRI

	 *

	 */

	public void setConnectionEsri() {

		setDatabase("jdbc:hsqldb:file:gvSIG/extensiones/org.gvsig.crs.extension/esri");

		setUser("sa");

		try {

			connect = super.getConnection();

		} catch (SQLException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

	}

	/**

	 * Establece la conexi?n con la base de datos de USR

	 *

	 */

	public void setConnectionUsr() {

		setDatabase("jdbc:hsqldb:file:gvSIG/extensiones/org.gvsig.crs.extension/usr");

		setUser("sa");

		try {

			connect = super.getConnection();

		} catch (SQLException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

	}

	public Connection getConnection(){

		return connect;

	}

	public void shutdown() throws SQLException {

        Statement st = connect.createStatement();

        // db writes out to files and performs clean shuts down

        // otherwise there will be an unclean shutdown

        // when program ends

        st.execute("SHUTDOWN");

        connect.close();    // if there are no other open connection

    }

	public synchronized void update(String expression) throws SQLException {

        Statement st = null;

        st = connect.createStatement();    // statements

        int i = st.executeUpdate(expression);    // run the query

        if (i == -1) {

            System.out.println("db error : " + expression);

        }

        st.close();

    }

	/*

	*//**

	 * Establece la conexi?n con la base de datos de la EPSG

	 *

	 *//*

	public void setConnectionEPSG() {

		try {			

			connect =  DriverManager.getConnection("jdbc:hsqldb:file:gvSIG/extensiones/org.gvsig.crs.extension/db_epsg", "sa", "");

		} catch (SQLException e1) {

			e1.printStackTrace();

		}

	}

	*//**

	 * Establece la conexi?n con la base de datos de la IAU2000

	 *

	 *//*

	public void setConnectionIAU2000() {

		try {			

			connect =  DriverManager.getConnection("jdbc:hsqldb:file:gvSIG/extensiones/org.gvsig.crs.extension/db_iau2000", "sa", "");			

		} catch (SQLException e1) {

			e1.printStackTrace();

		}

	}

	*//**

	 * Establece la conexi?n con la base de datos de ESRI

	 *

	 *//*

	public void setConnectionEsri() {

		try {			

			connect =  DriverManager.getConnection("jdbc:hsqldb:file:/home/jlgomez/gvSIGF2/_fwAndami/gvSIG/extensiones/org.gvsig.crs.extension/db_esri", "sa", "");			

		} catch (SQLException e1) {

			e1.printStackTrace();

		}

	}

	*//**

	 * Establece la conexi?n con la base de datos de USR

	 *

	 *//*

	public void setConnectionUsr() {

		try {			

			connect =  DriverManager.getConnection("jdbc:hsqldb:file:gvSIG/extensiones/org.gvsig.crs.extension/db_usr", "sa", "");			

		} catch (SQLException e1) {

			e1.printStackTrace();

		}

	}

	public Connection getConnection(){

		return connect;

	}

	public void shutdown() throws SQLException {

        Statement st = connect.createStatement();

        // db writes out to files and performs clean shuts down

        // otherwise there will be an unclean shutdown

        // when program ends

        st.execute("SHUTDOWN");

        connect.close();    // if there are no other open connection

    }*/

	 /**

     * Returns {@code true} if the database contains data. This method returns {@code false}

     * if an empty EPSG database has been automatically created by HSQL and not yet populated.

     */

    private static boolean dataExists(final Connection connection) throws SQLException {

        final ResultSet tables = connection.getMetaData().getTables(

                null, null, "IAU2000_%", new String[] {"TABLE"});

        final boolean exists = tables.next();

        tables.close();

        return exists;

    }

}

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

 *

 * Copyright (C) 2006 Instituto de Desarrollo Regional 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

 *

 *   Instituto de Desarrollo Regional (Universidad de Castilla La-Mancha)

 *   Campus Universitario s/n

 *   02071 Alabacete

 *   Spain

 *

 *   +34 967 599 200

 */

package es.idr.teledeteccion.connection;

import java.io.Serializable;

import java.sql.Connection;

import java.sql.ResultSet;

import java.sql.ResultSetMetaData;

import java.sql.SQLException;

import java.sql.Statement;

/**

 * Clase para realizar las consultas a las bases de datos utilizadas

 * 

 * @author Jos? Luis G?mez Mart?nez (jolugomar@gmail.com)

 *

 */

public class Query implements Serializable{

	/**

	 * 

	 */

	private static final long serialVersionUID = 1L;

	/**

	 * Acepta la sentencia y la conexi?n a la base de datos actual,

	 * realiza una selecci?n y devuelve los datos en un ResultSet

	 * @param sentence

	 * @param conn

	 * @return

	 */

	public static synchronized ResultSet select(String sentence, Connection conn){

		Statement st = null;

		ResultSet rs = null;

		try {

			st = conn.createStatement();

			rs = st.executeQuery(sentence);

			st.close();

		} catch (SQLException e) {

			e.printStackTrace();

		}		

		return rs;

	}	

}

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

 *

 * Copyright (C) 2006 Instituto de Desarrollo Regional 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

 *

 *   Instituto de Desarrollo Regional (Universidad de Castilla La-Mancha)

 *   Campus Universitario s/n

 *   02071 Alabacete

 *   Spain

 *

 *   +34 967 599 200

 */

package org.gvsig.crs;

/**

 * Clase que captura las excepciones cuando se definen los CRS

 * 

 * @author Miguel Garc?a Jim?nez (garciajimenez.miguel@gmail.com)

 *

 */

public class CrsException extends Exception {

	/**

	 * 

	 */

	private static final long serialVersionUID = 1L;

	public CrsException(Exception e) {

		super(e);

	}

}

/* gvSIG. Geographic Information System of the Valencian Government

 *

 * Copyright (C) 2007-2008 Infrastructures and Transports Department

 * of the Valencian Government (CIT)

 *

 * 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., 51 Franklin Street, Fifth Floor, Boston,

 * MA  02110-1301, USA.

 *

 */

package org.gvsig.crs;

import java.io.File;

import java.util.Arrays;

import java.util.Iterator;

import org.cresques.ProjectionLibrary;

import org.geotools.factory.FactoryRegistry;

import org.geotools.referencing.operation.MathTransformProvider;

import org.gvsig.andami.messages.NotificationManager;

import org.gvsig.fmap.crs.CRSFactory;

import org.gvsig.tools.library.AbstractLibrary;

import org.gvsig.tools.library.LibraryException;

/**

 * @author gvSIG Team

 * @version $Id$

 *

 */

public class JCRSLibrary extends AbstractLibrary {

    @Override

    public void doRegistration() {

        registerAsImplementationOf(ProjectionLibrary.class, 10);

    }

    @Override

    protected void doInitialize() throws LibraryException {

//    	cleanGeotoolsDatabase();

    	initializeGeotoolsProviders();

        CRSFactory.cp = new CrsFactory();

    }

    @Override

    protected void doPostInitialize() throws LibraryException {

        // nothing to do 

    }

    //Hemos cambiado el jar de epsg-hsql y tal vez esto ya no sea necesario

//	private void cleanGeotoolsDatabase(){

//        

//		// Parche para limpiar los archivos temporales de la base de datos cacheada de EPSG cuando ha podido quedar corrupta 

//		String tempFolder =  System.getProperty("java.io.tmpdir") ;

//        String pathSeparator = System.getProperty("file.separator") ;

//        String cachedDatabasesDir = tempFolder+pathSeparator+"Geotools"+pathSeparator+"Cached databases";

//        

//        File fileProperties = new File(cachedDatabasesDir+pathSeparator+"EPSG.properties");

//        File fileScript = new File(cachedDatabasesDir+pathSeparator+"EPSG.script");

//        File fileData = new File(cachedDatabasesDir+pathSeparator+"EPSG.data");

//        

//        boolean canDelete = true;

//        if ((fileProperties.exists() || fileScript.exists()) && !fileData.exists()){

//        	File cachedDatabasesDirFile = new File(cachedDatabasesDir);

//        	String[] fileNames = cachedDatabasesDirFile.list();

//        	for (int i = 0; i < fileNames.length; i++) {

//				File file = new File(cachedDatabasesDir+pathSeparator+fileNames[i]);

//				canDelete &= file.delete();

//			}

//        	if (!canDelete){

//        		NotificationManager.showMessageWarning("La base de datos de EPSG ha quedado corrupta, por favor salga de la aplicaci?n y borre manualmente el directorio "+cachedDatabasesDir, null);

//        	}

//        }

//        // Fin del parche

//        

//        /* Parche para prevenir una posible corrupci?n de la base de datos de EPSG que se

//         * produce en windows al tener m?s de una instancia de gvSIG abierta. 

//         */

//        if (fileProperties.exists()){

//        	if (!fileProperties.delete()){

//        		NotificationManager.showMessageWarning("La base de datos de EPSG ha podido quedar corrupta, por favor salga de la aplicaci?n y borre manualmente el directorio "+cachedDatabasesDir, null);

//        	}

//        }

//        // Fin del parche

//            

//

//	}

	private void initializeGeotoolsProviders(){

		/*

		 * Se eliminan del registro de geotools las proyecciones que aporta libJCrs para asegurar que

		 * son estas ?ltimas las que se utilizan.

		 */

		final Class[] categories = {org.geotools.referencing.operation.MathTransformProvider.class};

		FactoryRegistry registry = new FactoryRegistry(Arrays.asList(categories));

		Iterator providers = registry.getServiceProviders(MathTransformProvider.class);

		Iterator providers2 = null;

		MathTransformProvider provider = null;

		MathTransformProvider provider2 = null;

		while (providers.hasNext()){

			provider = (MathTransformProvider) providers.next();

			if(provider.nameMatches("IDR")){

				providers2 = registry.getServiceProviders(MathTransformProvider.class);

				while (providers2.hasNext()){

					provider2 = (MathTransformProvider) providers2.next();

					if(provider2.nameMatches(provider.getName().toString()) && !provider2.nameMatches("IDR")) {

						registry.deregisterServiceProvider(provider2, categories[0]);

					}

				}

			}

		}

		//TODO: Solo para depurar: BORRAR

		providers = registry.getServiceProviders(MathTransformProvider.class);

		while (providers.hasNext()) {

            provider = (MathTransformProvider) providers.next();

            System.out.println(provider.toString());

		}

		/*

		 *

		 */

	}

}

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

 *

 * Copyright (C) 2006 Instituto de Desarrollo Regional 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

 *

 *   Instituto de Desarrollo Regional (Universidad de Castilla La-Mancha)

 *   Campus Universitario s/n

 *   02071 Alabacete

 *   Spain

 *

 *   +34 967 599 200

 */

package org.gvsig.crs;

import java.awt.Color;

import java.awt.Graphics2D;

import java.awt.geom.Point2D;

import java.awt.geom.Rectangle2D;

import org.cresques.cts.ICoordTrans;

import org.cresques.cts.IDatum;

import org.cresques.cts.IProjection;

import org.cresques.geo.ViewPortData;

import org.gvsig.crs.proj.CrsProj;

import org.gvsig.crs.proj.JNIBaseCrs;

import org.gvsig.crs.proj.OperationCrsException;

import org.gvsig.fmap.crs.CRSFactory;

/**

 * Clase que construye el CRS a partir de la cadena WKT

 * @author Jos? Luis G?mez Mart?nez (jolugomar@gmail.com)

 * @author Diego Guerrero Sevilla (diego.guerrero@uclm.es)

 * @author Miguel Garc?a Jim?nez (garciajimenez.miguel@gmail.com)

 *

 */

public class Crs implements ICrs {

	private static final Color basicGridColor = new Color(64, 64, 64, 128);

	private Proj4 proj4;

	private String proj4String;

	private String trans;

	//private String transOrigin = "";

	private String abrev;

	private String name = "";

	private CrsProj crsProj;

	private CrsProj crsBase = null;

	private CrsWkt crsWkt;

	private int epsg_code = 23030;

	String sourceTrParams = null;

	String targetTrParams = null;

	String wkt = null;

	Color gridColor = basicGridColor;

	CRSDatum datum = null;	

	public Crs(int epsgCode, int aut) throws CrsException

	{

		String strEpsgCode = "";

		if (aut == 1) {

			strEpsgCode="EPSG:"+epsgCode;			

		} else if (aut == 2){

			strEpsgCode="ESRI:"+epsgCode;

		} else if (aut == 3){

			strEpsgCode="IAU2000:"+epsgCode;

		} else if (aut == 4) {

			strEpsgCode="USR:"+epsgCode;

		} else System.out.println("Error, autorithy err?neo");

		crsWkt=new CrsWkt(strEpsgCode);

		setWKT(crsWkt.getWkt());

	}

	/**

	 * Construye un CRS a partir del c?digo EPSG o de la cadena WKT.

	 * 

	 * @param code

	 * @throws CrsException

	 */

	public Crs(String code) throws CrsException {

		String fullCode;

		setWKT(code);

		if(code.charAt(0) == 'E' || code.charAt(0) == 'G' || code.charAt(0) == 'P')

			fullCode = code;

		else fullCode = "EPSG:"+code;

		String cod = "";

		if(code.length() < 15 ) {

			code = code.substring(code.indexOf(":")+1);

			try {

				//Creamos el objeto que tendra los diferentes parametros de la cadena Wkt

				crsWkt = new CrsWkt(fullCode);				

				setName(fullCode);

				setAbrev(fullCode);

			} catch (NumberFormatException e) {

				// TODO Auto-generated catch block

				e.printStackTrace();

			}

		}else {

			String aux;

			String code2 = "";

			for(int i = 0; i < code.length(); i++) {

				aux = ""+code.charAt(i);

				if(!aux.equals(" ")) {

					code2+=aux;

				}else {

					code2 += "";

				}

			}

			crsWkt = new CrsWkt(code2);

    		setName(fullCode);

    		//setAbrev(crsWkt.getName());

    		setAbrev(crsWkt.getAuthority()[0]+":"+crsWkt.getAuthority()[1]);	    	

    	}

		//	Asignar el datum y el crs base (en el caso de ser projectado):

		if (!(crsWkt.getSpheroid()[1].equals("")||crsWkt.getSpheroid()[2].equals(""))){

			double eSemiMajorAxis = Double.valueOf(crsWkt.getSpheroid()[1]).doubleValue();

			double eIFlattening = Double.valueOf(crsWkt.getSpheroid()[2]).doubleValue();

			datum = new CRSDatum(eSemiMajorAxis,eIFlattening);

			//crs base (en el caso de ser projectado):

			/*if(this.isProjected()){

				String proj4Base = "+proj=latlong +a=" + crsWkt.getSpheroid()[1] + " +rf=" + crsWkt.getSpheroid()[2] ;

				crsBase = new CrsProj(proj4Base);

			}*/

		}

	}	

	/**

	 *Construye un CRS a partir del c?digo EPSG o de la cadena WKT.

	 * 

	 * @param epsg_cod

	 * @param code

	 * @throws CrsException

	 */

	public Crs(int epsg_cod, String code) throws CrsException {

		String fullCode;

		setWKT(code);

		setCode(epsg_cod);

		if(code != null || code.charAt(0) == 'E' || code.charAt(0) == 'G' || code.charAt(0) == 'P')

			fullCode = code;

		else fullCode = "EPSG:"+code;

		String cod = "";

		if(code.length() < 15 ) {

			code = code.substring(code.indexOf(":")+1);

			try {

				//Creamos el objeto que tendra los diferentes parametros de la cadena Wkt

				crsWkt = new CrsWkt(fullCode);				

				setName(fullCode);

				setAbrev(fullCode);

			} catch (NumberFormatException e) {

				// TODO Auto-generated catch block

				e.printStackTrace();

			}

		}else {

			String aux;

			String code2 = "";

			for(int i = 0; i < code.length(); i++) {

				aux = ""+code.charAt(i);

				if(!aux.equals(" ")) {

					code2+=aux;

				}else {

					code2 += "";

				}

			}

			crsWkt = new CrsWkt(code2);

			/*

			 * Arreglo temporal para ver si funcionan de la iau2000 las proyecciones

			 * cilindrica equidistante y oblicua cilindrica equidistante

			 * que no estan en gdal, por lo que asignaremos directamente su cadena

			 * en proj4 para trabajar con ellas

			 */

			if (!crsWkt.getProjection().equals("Equidistant_Cylindrical") && !crsWkt.getProjection().equals("Oblique_Cylindrical_Equal_Area")){		    		    		

	    		setName(fullCode);

	    		//setAbrev(crsWkt.getName());

	    		setAbrev(crsWkt.getAuthority()[0]+":"+crsWkt.getAuthority()[1]);		    	 

			}

			else if (crsWkt.getProjection().equals("Equidistant_Cylindrical")){

				String spheroid1 = crsWkt.getSpheroid()[1];

				String spheroid2 = crsWkt.getSpheroid()[2];

				String centralMeridian = "";

				String falseNorthing = "";

				String falseEasting = "";

				String standardParallel1 = "0.0";

				for (int i=0; i<crsWkt.getParam_name().length;i++){

					if (crsWkt.getParam_name()[i].equals("Central_Meridian"))

						centralMeridian = crsWkt.getParam_value()[i];	

					if (crsWkt.getParam_name()[i].equals("False_Easting"))

						falseEasting = crsWkt.getParam_value()[i];					

					if(crsWkt.getParam_name()[i].equals("False_Northing"))

						falseNorthing = crsWkt.getParam_value()[i];					

					if (crsWkt.getParam_name()[i].equals("Standard_Parallel_1"))

						standardParallel1 = crsWkt.getParam_value()[i];					

				}

				if (spheroid2.equals("0.0")){

					proj4String = "+proj=eqc +a=" + spheroid1 +

					" +lon_0="+ centralMeridian + " +x_0="+falseEasting+" +y_0="+falseNorthing +

					" +lat_ts="+standardParallel1;

				} else {

					proj4String = "+proj=eqc +a="+ spheroid1 +" +rf=" + spheroid2 +

					" +lon_0="+ centralMeridian + " +x_0="+falseEasting+" +y_0="+falseNorthing +

					" +lat_ts="+standardParallel1;

				}				

				setName(fullCode);

	    		//setAbrev(crsWkt.getName());

	    		setAbrev(crsWkt.getAuthority()[0]+":"+crsWkt.getAuthority()[1]);

			}

			else if (crsWkt.getProjection().equals("Oblique_Cylindrical_Equal_Area")){

				String spheroid1 = crsWkt.getSpheroid()[1];

				String spheroid2 = crsWkt.getSpheroid()[2];

				String centralMeridian = "";

				String falseNorthing = "";

				String falseEasting = "";

				String standardParallel1 = "";

				String standardParallel2 = "0.0";

				for (int i=0; i<crsWkt.getParam_name().length;i++){					

					if (crsWkt.getParam_name()[i].equals("Central_Meridian"))

						centralMeridian = crsWkt.getParam_value()[i];	

					if (crsWkt.getParam_name()[i].equals("False_Easting"))

						falseEasting = crsWkt.getParam_value()[i];					

					if(crsWkt.getParam_name()[i].equals("False_Northing"))

						falseNorthing = crsWkt.getParam_value()[i];					

					if (crsWkt.getParam_name()[i].equals("Standard_Parallel_1"))

						standardParallel1 = crsWkt.getParam_value()[i];

					if (crsWkt.getParam_name()[i].equals("Standard_Parallel_2"))

						standardParallel2 = crsWkt.getParam_value()[i];		

				}

				if (spheroid2.equals("0.0")){

					proj4String = "+proj=ocea +a="+ spheroid1  +

					" +lon_0="+ centralMeridian + " +x_0="+falseEasting+" +y_0="+falseNorthing +

					" +lat_1="+standardParallel1+" +lat_2="+standardParallel2+" +lon_1=long_1" +

					" +lon_2=long_2 +no_defs";

				} else {

					proj4String = "+proj=ocea +a="+ spheroid1 + " +rf=" + spheroid2 +

					" +lon_0="+ centralMeridian + " +x_0="+falseEasting+" +y_0="+falseNorthing +

					" +lat_1="+standardParallel1+" +lat_2="+standardParallel2+" +lon_1=long_1" +

					" +lon_2=long_2 +no_defs";

				}				

				setName(fullCode);

	    		//setAbrev(crsWkt.getName());

	    		setAbrev(crsWkt.getAuthority()[0]+":"+crsWkt.getAuthority()[1]);

			}

		}

		//	Asignar el datum:

		if (!(crsWkt.getSpheroid()[1].equals("")||crsWkt.getSpheroid()[2].equals(""))){

			double eSemiMajorAxis = Double.valueOf(crsWkt.getSpheroid()[1]).doubleValue();

			double eIFlattening = Double.valueOf(crsWkt.getSpheroid()[2]).doubleValue();

			datum = new CRSDatum(eSemiMajorAxis,eIFlattening);

			// Crs base (en el caso de ser projectado):

			/*if(this.isProjected()){

				String proj4Base = "+proj=latlong +a=" + crsWkt.getSpheroid()[1] + " +rf=" + crsWkt.getSpheroid()[2] ;

				crsBase = new CrsProj(proj4Base);

			}*/

		}

	}

	/**

	 * Construye un CRS a partir del c?digo EPSG o de la cadena WKT.

	 * En el caso de WKT le a?ade a la cadena proj4 el string params.

	 * 

	 * @param epsg_cod

	 * @param code

	 * @param params

	 * @throws CrsException

	 */

	public Crs(int epsg_cod, String code,String params) throws CrsException {

		String fullCode;

		setCode(epsg_cod);

		setWKT(code);

		if(code.charAt(0) == 'E' || code.charAt(0) == 'G' || code.charAt(0) == 'P')

			fullCode = code;

		else fullCode = "EPSG:"+code;

		String cod = "";

		if(code.length() < 15 ) {

			code = code.substring(code.indexOf(":")+1);

			try {

				//Creamos el objeto que tendra los diferentes parametros de la cadena Wkt

				crsWkt = new CrsWkt(fullCode);

				setName(fullCode);

				setAbrev(fullCode);

			} catch (NumberFormatException e) {

				// TODO Auto-generated catch block

				e.printStackTrace();

			}

		}else {

			String aux;

			String code2 = "";

			for(int i = 0; i < code.length(); i++) {

				aux = ""+code.charAt(i);

				if(!aux.equals(" ")) {

					code2+=aux;

				}else {

					code2 += "";

				}

			}

			crsWkt = new CrsWkt(code2);

    		setName(fullCode);

    		//setAbrev(crsWkt.getName());

    		setAbrev(crsWkt.getAuthority()[0]+":"+crsWkt.getAuthority()[1]);

		}

		//	Asignar el datum:

		if (!(crsWkt.getSpheroid()[1].equals("")||crsWkt.getSpheroid()[2].equals(""))){

			double eSemiMajorAxis = Double.valueOf(crsWkt.getSpheroid()[1]).doubleValue();

			double eIFlattening = Double.valueOf(crsWkt.getSpheroid()[2]).doubleValue();

			datum = new CRSDatum(eSemiMajorAxis,eIFlattening);

			// Crs base (en el caso de ser projectado):

			/*if(this.isProjected()){

				String proj4Base = "+proj=latlong +a=" + crsWkt.getSpheroid()[1] + " +rf=" + crsWkt.getSpheroid()[2] ;

				crsBase = new CrsProj(proj4Base);

			}*/

		}

	}

/*	public Crs(CrsEpsg source) throws CrsException {

		crsProj4 = source;

	}*/

	/**

	 * @param code 

	 */

	public void setTrans(String code) {

		trans = code;

		changeTrans(trans);

	}

	/**

	 * 

	 * @param code

	 */

	public void changeTrans(String code) {

		getCrsProj().changeStrCrs(code);

	}

	/**

	 * 

	 */

	public String getAbrev() {

		return abrev;

	}

	/**

	 * 

	 * @param code

	 */

	protected void setAbrev(String code) {

		abrev = code;

	}

	/**

	 * 

	 * @param nom

	 */

	public void setName(String nom) {

		name = nom;

	}

	/**

	 * @return

	 */

	public IDatum getDatum() {

		return datum;

	}

	/**

	 * @return

	 */

	public CrsWkt getCrsWkt() {

		return crsWkt;

	}

	/**

	 * 

	 * @param wkt

	 */

	public void setWKT(String wkt){

		this.wkt = wkt;

	}

	/**

	 * @return

	 */

	public String getWKT(){

		return this.wkt;

	}

	/**

	 * 

	 * @param params

	 */

	public void setTransformationParams(String sourceParams, String targetParams){

		this.sourceTrParams = sourceParams;

		this.targetTrParams = targetParams;

	}

	/**

	 * Devuelve los parametros de la transformacion del crs fuente

	 * @return

	 */

	public String getSourceTransformationParams() {

		return this.sourceTrParams;

	}

	/**

	 * Devuelve los parametros de la transformacion del crs destino

	 * @return

	 */

	public String getTargetTransformationParams() {

		return this.targetTrParams;

	}

	/**

	 * @return

	 */

	public String getTransformationParams(){

		return this.sourceTrParams;

	}

	/**

	 * 

	 * @return

	 */

	public CrsProj getCrsProj() {

		if (crsProj == null)

			try {

				crsProj = new CrsProj(getProj4String());

			} catch (CrsException e) {

				// TODO Auto-generated catch block

				e.printStackTrace();

			} 

		return crsProj;

	}

	/**

	 * @param x

	 * @param y

	 * @return

	 */

	public Point2D createPoint(double x, double y) {

		return new Point2D.Double(x,y);

	}

	/**

	 * @param g

	 * @param vp

	 */

	public void drawGrid(Graphics2D g, ViewPortData vp) {

		// TODO Auto-generated method stub

	}

	/**

	 * @param c

	 */

	public void setGridColor(Color c) {

		gridColor = c;

	}

	/**

	 * @return

	 */

	public Color getGridColor() {

		return gridColor;

	}

	/**

	 * @param dest

	 * @return

	 */

	public ICoordTrans getCT(IProjection dest) {

		try {

			if (dest == this)

				return null;

			COperation operation = null;

			if(((ICrs)dest).getSourceTransformationParams() != null || ((ICrs)dest).getTargetTransformationParams() != null) 

				operation = new COperation(this, (ICrs)dest,((ICrs)dest).getTargetTransformationParams(),((ICrs)dest).getSourceTransformationParams());

			else

				operation = new COperation(this, (ICrs)dest,sourceTrParams,targetTrParams);

			return operation;

		}catch (CrsException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

			return null;

		}

		/*

		try {

			operation = new COperation(this, (ICrs)dest);

		} catch (CrsException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

		if (getTransformationParams()!=null){

			if (isTargetParams())

				operation.setParamsCrsProj(new CrsProj(crsDest.getProj4String()+getTransformationParams()), true);

			else

				operation.setParamsCrsProj(new CrsProj(getProj4String()+getTransformationParams()), false);

			return operation;

		}

		return operation;	*/	

	}

	/**

	 * @param pt

	 * @return

	 */

	public Point2D toGeo(Point2D pt) {

		if (isProjected()){

			double x[] = {pt.getX()};

			double y[] = {pt.getY()};

			double z[] = {0D};

			try {

				JNIBaseCrs.operate( x , y, z,

						getCrsProj(),crsBase);

			} catch (OperationCrsException e) {

				// TODO Auto-generated catch block

				e.printStackTrace();

			}

			return new Point2D.Double(x[0],y[0]);

		}

		else

			return pt;

	}

	/**

	 * @param gPt

	 * @param mPt

	 * @return

	 */

	public Point2D fromGeo(Point2D gPt, Point2D mPt) {

		// TODO Auto-generated method stub

		return null;

	}

	/**

	 * @return

	 */

	public boolean isProjected() {

		return !getCrsProj().isLatlong();		

	}

	/**

	 * @param minX

	 * @param maxX

	 * @param width

	 * @param dpi

	 * @return

	 */

	public double getScale(double minX, double maxX, double width, double dpi) {

		double scale = 0D;

		//Esto ahora se hace fuera

//        if (!isProjected()) { // Es geogr?fico; calcula la escala.

//            scale = ((maxX - minX) * // grados

//

//            // 1852.0 metros x minuto de meridiano

//            (dpi / 2.54 * 100.0 * 1852.0 * 60.0)) / // px / metro

//                    width; // pixels

//        }

//        else{

        	 scale = ((maxX - minX) * // metros

                    (dpi / 2.54 * 100.0)) / // px / metro

                    width; // pixels

//        }

        return scale;

	}

	public double getScale(double minX, double maxX, double minY, double maxY, double width, double dpi) {

		double scale = 0D;

		double incX = (maxX-minX);

		if (!isProjected()) {

			double a = getDatum().getESemiMajorAxis();

			double invF = getDatum().getEIFlattening();

			double meanY = (minY+maxY)/2.0;

			double radius = 0.0;

			if (invF == Double.POSITIVE_INFINITY){

				radius = a;

			}

			else{

				double e2 = 2.0/invF-Math.pow(1.0/invF,2.0);

				radius = a/Math.sqrt(1.0-e2*Math.pow(Math.sin(meanY*Math.PI/180.0),2.0))*Math.cos(meanY*Math.PI/180.0);

			}

			incX *= Math.PI/180.0*radius;

		}

		scale = (incX * // metros

				(dpi / 2.54 * 100.0)) / // px / metro

					width; // pixels

        return scale;

	}

	/**

	 * 

	 * @param epsg_cod

	 */

	public void setCode(int epsg_cod){

		epsg_code = epsg_cod;

	}

	/**

	 * @return

	 */

	public int getCode() {

		// TODO Auto-generated method stub

		return epsg_code;

	}

	public String getFullCode() {

		if (this.sourceTrParams == null && this.targetTrParams == null)

			return getAbrev();

		String sourceParams = "";

		String targetParams = "";

		if (sourceTrParams != null)

			sourceParams = sourceTrParams;

		if (targetTrParams != null)

			targetParams = targetTrParams;

		return getAbrev()+":proj@"+sourceParams+"@"+targetParams;

	}

    public Rectangle2D getExtent(Rectangle2D extent,double scale,double wImage,double hImage,double mapUnits,double distanceUnits,double dpi) {

    	double w =0;

		double h =0;

		double wExtent =0;

		double hExtent =0;

//    	if (isProjected()) {

			w = ((wImage / dpi) * 2.54);

			h = ((hImage / dpi) * 2.54);

			wExtent =w * scale*distanceUnits/ mapUnits;

			hExtent =h * scale*distanceUnits/ mapUnits;

			//En las no proyectadas, ahora viene el calculo hecho

			//al haber incluido los grados en el array DISTANCETRANS2METER del MapContext

//		}else {

//			w = ((wImage / dpi) * 2.54);

//			h = ((hImage / dpi) * 2.54);

//			wExtent =(w*scale*distanceUnits)/ (mapUnits*1852.0*60.0);

//			hExtent =(h*scale*distanceUnits)/ (mapUnits*1852.0*60.0);

//		}

    	double xExtent = extent.getCenterX() - wExtent/2;

		double yExtent = extent.getCenterY() - hExtent/2;

		Rectangle2D rec=new Rectangle2D.Double(xExtent,yExtent,wExtent,hExtent);

    	return  rec;

    }

	private Proj4 getProj4() {

		if (proj4 == null)

			try {

				proj4 = new Proj4();

			} catch (CrsException e) {

				// TODO Auto-generated catch block

				e.printStackTrace();

			}

		return proj4;

	}

	/** 

	 * 

	 * @return Cadena proj4 Correspondiente al CRS.

	 * @throws CrsException 

	 */

	public String getProj4String() throws CrsException {

		if (proj4String == null)

			proj4String = getProj4().exportToProj4(this);			

		return proj4String;

	}

	public Object clone() throws CloneNotSupportedException {

	       return CRSFactory.getCRS( this.getFullCode() );

	    }

}

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

 *

 * Copyright (C) 2006 Instituto de Desarrollo Regional 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

 *

 *   Instituto de Desarrollo Regional (Universidad de Castilla La-Mancha)

 *   Campus Universitario s/n

 *   02071 Alabacete

 *   Spain

 *

 *   +34 967 599 200

 */

package org.gvsig.crs.ogr;

import java.sql.ResultSet;

import java.sql.SQLException;

import es.idr.teledeteccion.connection.EpsgConnection;

import es.idr.teledeteccion.connection.Query;

/**

 * Clase que consigue los par?metros necesarios de la base de

 * datos de la EPSG una vez escogido un CRS, y que ser?n

 * los que se utilicen para la generaci?n de la cadena wkt

 * correcta para generar el CRS.

 * 

 * @author Jos? Luis G?mez Mart?nez (jolugomar@gmail.com)

 *

 */

public class CrsEPSG {

	int epsg_code;

	boolean crs_source;

	int source_code;

	int projection_code;	

	int coord_op_code;	

	int datum_code;

	int ellipsoid_code, prime_meridian_code;

	String PROJCS = null;

	String GEOGCS = null;

	String DATUM = null;

	String[] SPHEROID = null;

	String[] PRIMEM = null;

	String UNIT_A = "degree";

	String UNIT_B = "Meter";

	double unit_BValue = 1.0;

	String PROJECTION = null;

	//ArrayList PARAMETER = null;

	String[] param_name = null;

	String[] param_value = null;

	String[] AUTHORITY = null;

	EpsgConnection connect;

	private int coordSysCode;

	public CrsEPSG(){		

	}	

	public CrsEPSG(int code, boolean source_yn, int source_cod, int coord_op_cod, EpsgConnection conn){

		epsg_code = code;

		crs_source = source_yn;

		source_code = source_cod;		

		coord_op_code = coord_op_cod;

		connect = conn;

		SPHEROID = new String[3];

		PRIMEM = new String[2];

		AUTHORITY = new String[2];

		initialize();

	}

	public void initialize(){

		ResultSet result;

		ResultSet result2;

		/*

		 * Si es base solamente realizamos una consulta, y ya tenemos el Datum

		 * (ver en caso de ser fuente que PROJCS no tiene que estar definido)

		 */

		if (crs_source){

			String sentence = "SELECT coord_ref_sys_name, datum_code, coord_sys_code " +

							  "FROM epsg_coordinatereferencesystem " +	                              

							  "WHERE coord_ref_sys_code = " + epsg_code;

			result = Query.select(sentence,connect.getConnection());

			try {

				while (result.next()){

					GEOGCS = result.getString("coord_ref_sys_name");

					datum_code = Integer.parseInt(result.getString("datum_code"));

					coordSysCode = result.getInt("coord_sys_code");

				}

			} catch (SQLException e) {

				e.printStackTrace();

			}			

		}

		else{

			String sentence = "SELECT coord_ref_sys_name, coord_sys_code " +

							  "FROM epsg_coordinatereferencesystem " +	                              

							  "WHERE coord_ref_sys_code = " + epsg_code;

			result = Query.select(sentence,connect.getConnection());

			String sentence2 = "SELECT coord_ref_sys_name, datum_code " +

							  "FROM epsg_coordinatereferencesystem " +	                              

							  "WHERE coord_ref_sys_code = " + source_code;

			result2 = Query.select(sentence2,connect.getConnection());

			try {

				while (result.next()){

					PROJCS = result.getString("coord_ref_sys_name");

					coordSysCode = result.getInt("coord_sys_code");

				}

				while (result2.next()){

					datum_code = Integer.parseInt(result2.getString("datum_code"));

					GEOGCS = result2.getString("coord_ref_sys_name");

				}

			} catch (SQLException e) {

				e.printStackTrace();

			}

		}

		/*

		 * hasta aqui hemos rellenado projcs, geocs y tenemos el codigo de datum

		 * vamos a obtener los c?digos del elipsoide y el prime_meridian

		 */

		String sentence = "SELECT datum_name, ellipsoid_code, prime_meridian_code " +

						  "FROM epsg_datum " +	                              

						  "WHERE datum_code = " + datum_code;

		result = Query.select(sentence,connect.getConnection());

		try {

			while (result.next()){

				DATUM = result.getString("datum_name");

				ellipsoid_code = Integer.parseInt(result.getString("ellipsoid_code"));	

				prime_meridian_code = Integer.parseInt(result.getString("prime_meridian_code"));

			}			

		} catch (SQLException e) {

			e.printStackTrace();

		}

		/*

		 * Ahora vamos a coger varios campos del elipsoide y operaremos con ellos, si la

		 * informaci?n de los mismos no esta en las unidades de medida que utilizaremos.

		 */

		sentence = "SELECT ellipsoid_name, semi_major_axis, inv_flattening, uom_code, " +

					"semi_minor_axis, ellipsoid_shape " +

					  "FROM epsg_ellipsoid " +	                              

					  "WHERE ellipsoid_code = " + ellipsoid_code;

		result = Query.select(sentence,connect.getConnection());

		SPHEROID = getEllipsoid(result);

		sentence = "SELECT prime_meridian_name, greenwich_longitude, uom_code " +

					  "FROM epsg_primemeridian " +	                              

					  "WHERE prime_meridian_code = " + prime_meridian_code;

		result = Query.select(sentence,connect.getConnection());

		PRIMEM = getPrimeMeridian(result);

		AUTHORITY = getAuthority(epsg_code);

		/*

		 * si tiene proyecci?n porque es proyectado, la buscamos

		 */

		if (!crs_source){			

			sentence = "SELECT coord_op_method_code " +

						"FROM epsg_coordoperation " +

						"WHERE coord_op_code = " + coord_op_code;

			result = Query.select(sentence,connect.getConnection());

			try {

				while (result.next()){

					projection_code = result.getInt("coord_op_method_code");