Application: gvSIG desktop » gvSIG projections

package es.idr.projection.test;

import org.gvsig.crs.CrsException;

import org.gvsig.crs.CrsFactory;

import org.gvsig.crs.ICrs;

import org.gvsig.crs.Proj4;

import org.opengis.referencing.FactoryException;

import org.opengis.referencing.NoSuchAuthorityCodeException;

public class TestGeotools21 {

	/**

	 * @param args

	 */

	public static void main(String[] args) {

		try {	        	

    		String strEpsgCode;

			// Notas:

			// - transverse mercator south orientated no esta soportada por Proj4

    		// - Buscar ejemplo de la CylindricalEqualArea

    		// - Buscar ejemplo de la AzimuthalEquidistant

    		// - Buscar ejemplo de la Gnomonic

    		// - Buscar ejemplo de la Bonne

    		// - Buscar ejemplo de la Equirectangular

    		// - Buscar ejemplo de la EquidistantConic

    		// - Buscar ejemplo de la MillerCylindrical

    		// - Buscar ejemplo de la Mollweide

    		// - Buscar ejemplo de la EckertIV

    		// - Buscar ejemplo de la EckertVI

    		// - Buscar ejemplo de la Robinson

    		// - Buscar ejemplo de la VanDerGrinten

    		// - Buscar ejemplo de la Sinusoidal

    		// - Buscar ejemplo de la GallStereographic

    		// - Buscar ejemplo de la Goode

			//strEpsgCode="EPSG:2065"; // Krovak

    		//strEpsgCode="EPSG:2052";  // tm south orientated No en mi Proj4

    		//strEpsgCode="EPSG:2084";  // tm south orientated

    		//strEpsgCode="EPSG:32662";  // Plate Carree  -> Funciona, problema elipsoide

    											// Tomara como radio el semieje mayor del elipsoide

    		//strEpsgCode="EPSG:27200";  // NewZealandMapGrid -> Mirar lo de +datum

    		// He tenido que incluir un mecanismo para valores concretos de parametros

    		strEpsgCode="EPSG:32761";  // CassiniSoldner -> Falla +to_meter=..

    		//strEpsgCode="EPSG:2099";

    		//strEpsgCode="EPSG:2200";  // en lugar de stere pone sterea,igual 2171

    									// 2036 - falla, 32661

    		//strEpsgCode="EPSG:2163";    // Laea esferica, funciona

    		//strEpsgCode="EPSG:3035";    // ETRS-LAEA, funciona

    		//strEpsgCode="EPSG:29100";	// Polyconic, falla en geotools

    		//strEpsgCode="EPSG:2964";	// Albers-conic-equal-area

    		                            // debería fallar por cambio de unidad

    		//strEpsgCode="EPSG:24571";    // Oblicua de mercator, rf mas decimales

    		                               // quitar segundo alpha

    		//strEpsgCode="EPSG:29700";    // Oblicua de mercator, rf mas decimales

    		//strEpsgCode="EPSG:21780";		//somer, sin alpha

    		ICrs crs=null;

    		//CoordinateReferenceSystem source;

    		//source = CRS.decode(strEpsgCode);

        	crs = new CrsFactory().getCRS(strEpsgCode); 

        	System.out.println(crs.getWKT());        	

			System.out.println(crs.getProj4String()); //proj4.exportToProj4(crs));

		} catch (CrsException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

	}

}

package es.idr.test;

import org.geotools.referencing.CRS;

import org.opengis.referencing.FactoryException;

import org.opengis.referencing.NoSuchAuthorityCodeException;

import org.opengis.referencing.crs.CoordinateReferenceSystem;

public class TestGeotools21 {

	/**

	 * @param args

	 */

	public static void main(String[] args) {

		try {	        	

    		String strEpsgCode="EPSG:27572";

    		CoordinateReferenceSystem source;

    		source = CRS.decode(strEpsgCode);

        	System.out.println(source.toWKT());

		} catch (NoSuchAuthorityCodeException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		} catch (FactoryException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

	}

}

package es.idr.test.connection;

import java.sql.ResultSet;

import java.sql.SQLException;

import es.idr.teledeteccion.connection.EpsgConnection;

import es.idr.teledeteccion.connection.Query;

public class TestConnection {

	public TestConnection () {

		testEpsgDataBase();

		testEsriDataBase();

		testIau2000DataBase();

	}

	/**

	 * @param args

	 */

	public static void main(String[] args) {

		// TODO Auto-generated method stub

		TestConnection tc = new TestConnection();

	}

	private void testEpsgDataBase() {

		EpsgConnection connect = new EpsgConnection();

		connect.setConnectionEPSG();

		String sentence = "SELECT coord_ref_sys_name " +

							"FROM epsg_coordinatereferencesystem " +	                              

							"WHERE coord_ref_sys_code = 23030";

		ResultSet result;

        try {			

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

			result.next();

			System.out.println(result.getString("coord_ref_sys_name"));

		} catch (SQLException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

	}

	private void testEsriDataBase (){

		EpsgConnection connect = new EpsgConnection();

		connect.setConnectionEsri();

		String sentence = "SELECT esri_proj " +							 

						  "FROM ESRI " +	                              

						  "WHERE esri_code = 23030";

		ResultSet result = null;

    	try {

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

			result.next();

			System.out.println(result.getString("esri_proj"));

		} catch (SQLException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

	}

	private void testIau2000DataBase() {

		EpsgConnection connect = new EpsgConnection();

		connect.setConnectionIAU2000();

		String sentence = "SELECT iau_geog " +							 

						  "FROM IAU2000 " +	                              

						  "WHERE iau_code = 19900";

		ResultSet result = null;    	

		try {

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

			result.next();

			System.out.println(result.getString("iau_geog"));

		} catch (SQLException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

	}

}

package es.idr.test.connection;

import java.sql.Connection;

import java.sql.DriverManager;

import java.sql.SQLException;

import java.sql.Statement;

public class UsrConnection {

	Connection connect;

	public UsrConnection() {

		try {

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

		} catch (ClassNotFoundException e) {

			e.printStackTrace();

		}

		try {			

			connect =  DriverManager.getConnection("jdbc:hsqldb:file:/home/jlgomez/gvSIGF2/gt2-usr-hsql/src/org/geotools/referencing/factory/usr/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

    }

	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();

    }    

}

package es.idr.test;

import java.sql.Connection;

import java.sql.PreparedStatement;

import java.sql.ResultSet;

import java.sql.SQLException;

import java.util.ArrayList;

import javax.units.Unit;

import org.geotools.referencing.CRS;

import org.geotools.referencing.crs.DefaultGeographicCRS;

import org.geotools.referencing.crs.DefaultProjectedCRS;

import org.geotools.referencing.datum.DefaultGeodeticDatum;

import org.gvsig.crs.Crs;

import org.gvsig.crs.CrsException;

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

import org.opengis.referencing.FactoryException;

import org.opengis.referencing.NoSuchAuthorityCodeException;

import org.opengis.referencing.crs.CoordinateReferenceSystem;

import org.opengis.referencing.crs.GeographicCRS;

public class ValidateEpsgCRSs {

	public ValidateEpsgCRSs() {

	}

	/**

	 * @param args

	 */

	public static void main(String[] args) {

		HSQLDataSource ds = new HSQLDataSource();

		Connection conn;

		try {

			conn = ds.getConnection();

			final PreparedStatement stmt;

	        stmt = conn.prepareStatement("SELECT COORD_REF_SYS_CODE"

	                                      + " FROM epsg_coordinatereferencesystem");

	        ResultSet result;

			result = stmt.executeQuery();

	        ArrayList res = new ArrayList();

	        ArrayList units = new ArrayList();

	        ArrayList validos = new ArrayList();

	        while (result.next()) {

	        	int cod = result.getInt("COORD_REF_SYS_CODE");

        		String strEpsgCode="EPSG:"+cod;		

        		CoordinateReferenceSystem source;

        			try {

						source = CRS.decode(strEpsgCode);

						validos.add(String.valueOf(cod));

						if (source instanceof DefaultProjectedCRS) {

							String str;

							boolean esta = false;

							DefaultProjectedCRS sour = (DefaultProjectedCRS) source;

							Unit u = ((DefaultGeodeticDatum) sour.getDatum()).getPrimeMeridian().getAngularUnit();

							units.add(u.toString());

							/*for (int j = 0; j< units.size(); j++) {

					        	if (units.get(j).equals(u.toString())){

					        		esta = true;

					        		break;

					        	}						        	

					        }

							if (!esta){

								units.add(u.toString());

								esta = false;

							}*/

							for (int i=0; i< sour.getConversionFromBase().getParameterValues().values().size();i++) {

								str = sour.getConversionFromBase().getParameterValues().values().get(i).toString();

								u = sour.getConversionFromBase().getParameterValues().parameter(str.split("=")[0]).getUnit();

								units.add(u.toString());

							}

						}

						if (source instanceof DefaultGeographicCRS) {

							DefaultGeographicCRS sour = (DefaultGeographicCRS) source;

							boolean esta = false;

							Unit u = ((DefaultGeodeticDatum) sour.getDatum()).getPrimeMeridian().getAngularUnit();

							units.add(u.toString());							

						}

					} catch (NoSuchAuthorityCodeException e) {						

						res.add(String.valueOf(cod)+" --> "+e);

						// TODO Auto-generated catch block

						//e.printStackTrace();

					} catch (FactoryException e) {

						res.add(String.valueOf(cod)+" --> "+e);

						// TODO Auto-generated catch block

						//e.printStackTrace();

					}

	        }

	        for (int i = 0; i< res.size(); i++) {

	        	System.out.println(res.get(i));	        	

	        }	       

	        for (int i = 0; i< units.size(); i++) {

	        	System.out.println(units.get(i));

	        }

	        System.out.println("Numero CRSs validos: " + validos.size());

	        System.out.println("Numero CRSs fallidos: " + res.size());

	        System.out.println("Numero medidas: " + units.size());

		} catch (SQLException e1) {

			// TODO Auto-generated catch block

			e1.printStackTrace();

		}	

	}

}

			/*CoordinateReferenceSystem sourceCRS = CRS.decode("EPSG:4326");

			CoordinateReferenceSystem targetCRS = CRS.decode("EPSG:23032");

			MathTransform transform = CRS.findMathTransform(sourceCRS, targetCRS);

			System.out.println("Dimensiones: "+transform.getSourceDimensions());

			double points[] = {2.5,3.0};

			System.out.println("Puntos fuente: "+ points[0]+"  "+points[1]);

			transform.transform(points, 0, points, 0, 1);

			System.out.println("Puntos destino: "+ points[0]+"  "+points[1]);

			int epsgCount = 0;

			System.out.println("CRSs Disponibles en EPSG: \n");

			Set codes = CRS.getSupportedCodes("EPSG");

			for (Iterator iter = codes.iterator(); iter.hasNext();) {

			String code = (String) iter.next();

			System.out.println("EPSG:"+code);

			epsgCount++;

			}

			System.out.println("N�mero de CRSs en EPSG: "+epsgCount);

			try {

			System.out.println("Descripci�n de EPSG:23030:  \n"+CRS.getAuthorityFactory(false).getDescriptionText("EPSG:23030"));

			} catch (FactoryRegistryException e) {

//			 TODO Auto-generated catch block

			e.printStackTrace();

			} catch (FactoryException e) {

//			 TODO Auto-generated catch block

			e.printStackTrace();

			}

		} catch (NoSuchAuthorityCodeException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		} catch (FactoryException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}*/  

package es.idr.test;

import java.sql.SQLException;

import es.idr.teledeteccion.connection.EpsgConnection;

public class InsertUserDataCrs {

	/**

	 * @param args

	 */

	public static void main(String[] args) {

		// TODO Auto-generated method stub

		EpsgConnection conn = new EpsgConnection();

		conn.setConnectionUsr();

		String sentence = "INSERT INTO USR VALUES(2000,'PROJCS[\"Anguilla_1957_British_West_Indies_Grid\",GEOGCS[\"GCS_Anguilla_1957\",DATUM[\"D_Anguilla_1957\",SPHEROID[\"Clarke_1880_RGS\",6378249.145,293.465]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"False_Easting\",400000.0],PARAMETER[\"False_Northing\",0.0],PARAMETER[\"Central_Meridian\",-62.0],PARAMETER[\"Scale_Factor\",0.9995],PARAMETER[\"Latitude_Of_Origin\",0.0],UNIT[\"Meter\",1.0]]','Anguilla_1957_British_West_Indies_Grid','GCS_Anguilla_1957','D_Anguilla_1957')";

		try {

			conn.update(sentence);

			conn.shutdown();

		} catch (SQLException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

	}

}

package es.idr.test;

import java.io.BufferedWriter;

import java.io.FileWriter;

import java.io.IOException;

import java.sql.Connection;

import java.sql.PreparedStatement;

import java.sql.ResultSet;

import java.sql.SQLException;

import java.util.ArrayList;

import javax.units.ConversionException;

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

import org.geotools.referencing.wkt.UnformattableObjectException;

import org.gvsig.crs.Crs;

import org.gvsig.crs.CrsException;

import org.gvsig.crs.CrsFactory;

import org.gvsig.crs.CrsGT;

import org.gvsig.crs.ICrs;

import org.gvsig.crs.Proj4;

import org.opengis.referencing.NoSuchAuthorityCodeException;

import org.opengis.referencing.crs.CoordinateReferenceSystem;

public class TestWKTandProj4Epsg {

	public static void main(String[] args) {

		HSQLDataSource ds = new HSQLDataSource();

		Connection conn;

		try {

			conn = ds.getConnection();

			int indice = 0;

			int fallo = 0;

			final PreparedStatement stmt;

	        stmt = conn.prepareStatement("SELECT COORD_REF_SYS_CODE, COORD_REF_SYS_NAME"

	        		          				+" FROM epsg_coordinatereferencesystem ");

	        		          				//+" WHERE COORD_REF_SYS_CODE>29372");

	        ResultSet result;

			result = stmt.executeQuery();

	        String filename = "salida.txt";

	        BufferedWriter bw =  new BufferedWriter(new FileWriter(filename));

	        while (result.next()) {

	        	ICrs crs = null;

	        	int cod = result.getInt("COORD_REF_SYS_CODE");

	        	try {

					crs = new CrsFactory().getCRS("EPSG:"+cod);

					if (crs instanceof CrsGT) {

//						//continue;

						if (crs.getProj4String() == null){

							bw.write("C?digo: "+cod+"\tNombre: "+result.getString("COORD_REF_SYS_NAME")+"\n");

				        	bw.write("Cadena WKT: "+crs.getWKT()+"\n");

				        	bw.write("Cadena Proj4: No cadena proj4 \n");

				        	bw.write("Fallo: Fallo generando en cadena proj4 \n");		        	

				        	bw.write("\n\n");

						}

					}						

					else {

						if (crs != null) {

							bw.write("C?digo: "+cod+"\tNombre: "+result.getString("COORD_REF_SYS_NAME")+"\n");

				        	bw.write("Cadena WKT: "+crs.getWKT()+"\n");

				        	bw.write("Cadena Proj4: "+crs.getProj4String()+"\n");

				        	bw.write("Fallo: Creado con crs antiguo \n");		        	

				        	bw.write("\n\n");

						}

						else {

							bw.write("C?digo: "+cod+"\tNombre: "+result.getString("COORD_REF_SYS_NAME")+"\n");

							bw.write("Fallo: Falla tanto geotools como crs antiguo \n");		        	

				        	bw.write("\n\n");

						}

					}

				} catch (CrsException e) {

					bw.write("C?digo: "+cod+"\tNombre: "+result.getString("COORD_REF_SYS_NAME")+"\n");

		        	bw.write("Cadena WKT: No cadena wkt \n");

		        	bw.write("Cadena Proj4: No cadena proj4\n");

		        	bw.write("Fallo: "+e+" \n");		        	

		        	bw.write("\n\n");

				} catch (NumberFormatException e) {

					bw.write("C?digo: "+cod+"\tNombre: "+result.getString("COORD_REF_SYS_NAME")+"\n");

		        	bw.write("Cadena WKT: No cadena wkt \n");

		        	bw.write("Cadena Proj4: No cadena proj4\n");

		        	bw.write("Fallo: "+e+" \n");		        	

		        	bw.write("\n\n");

				} catch (ConversionException e) {

					bw.write("C?digo: "+cod+"\tNombre: "+result.getString("COORD_REF_SYS_NAME")+"\n");

		        	bw.write("Cadena WKT: No cadena wkt \n");

		        	bw.write("Cadena Proj4: No cadena proj4\n");

		        	bw.write("Fallo: "+e+" \n");		        	

		        	bw.write("\n\n");

				}catch (UnformattableObjectException e) {					

				}				

	        }

	        bw.flush();

	        bw.close();

		} catch (SQLException e1) {

			// TODO Auto-generated catch block

			e1.printStackTrace();

		} catch (IOException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}	

	}

}

package es.idr.test;

import java.sql.Connection;

import java.sql.PreparedStatement;

import java.sql.ResultSet;

import java.sql.SQLException;

import java.util.ArrayList;

import javax.units.Unit;

import org.geotools.referencing.CRS;

import org.geotools.referencing.crs.DefaultGeographicCRS;

import org.geotools.referencing.crs.DefaultProjectedCRS;

import org.geotools.referencing.datum.DefaultGeodeticDatum;

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

import org.opengis.referencing.FactoryException;

import org.opengis.referencing.NoSuchAuthorityCodeException;

import org.opengis.referencing.crs.CoordinateReferenceSystem;

import es.idr.teledeteccion.connection.EpsgConnection;

import es.idr.teledeteccion.connection.Query;

public class ValidateEsriCRSs {

	public ValidateEsriCRSs() {

	}

	/**

	 * @param args

	 */

	public static void main(String[] args) {

		//HSQLDataSource ds = new HSQLDataSource();

		EpsgConnection conn = new EpsgConnection();

		//Connection conn;

		conn.setConnectionEsri();

		try {

			//conn = ds.getConnection();

			//final PreparedStatement stmt;

	        /*stmt = conn.prepareStatement("SELECT esri_code"

	                                      + " FROM esri");

	        */

	        String sentence = "SELECT esri_code,esri_wkt " +							 

								  "FROM ESRI";

	        ResultSet result;

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

			//result = stmt.executeQuery();

	        ArrayList res = new ArrayList();

	        ArrayList units = new ArrayList();

	        ArrayList validos = new ArrayList();

	        while (result.next()) {

	        	int cod = result.getInt("esri_code");

	        	String wkt = result.getString("esri_wkt");

        		//String strEpsgCode="ESRI:"+cod;		

        		CoordinateReferenceSystem source;

        			/*try {

						source = CRS.parseWKT(wkt);

						validos.add(String.valueOf(cod));

						if (source instanceof DefaultProjectedCRS) {

							String str;

							boolean esta = false;

							DefaultProjectedCRS sour = (DefaultProjectedCRS) source;

							Unit u = ((DefaultGeodeticDatum) sour.getDatum()).getPrimeMeridian().getAngularUnit();

							units.add(u.toString());

							for (int j = 0; j< units.size(); j++) {

					        	if (units.get(j).equals(u.toString())){

					        		esta = true;

					        		break;

					        	}						        	

					        }

							if (!esta){

								units.add(u.toString());

								esta = false;

							}

							for (int i=0; i< sour.getConversionFromBase().getParameterValues().values().size();i++) {

								str = sour.getConversionFromBase().getParameterValues().values().get(i).toString();

								u = sour.getConversionFromBase().getParameterValues().parameter(str.split("=")[0]).getUnit();

								units.add(u.toString());

							}

						}

						if (source instanceof DefaultGeographicCRS) {

							DefaultGeographicCRS sour = (DefaultGeographicCRS) source;

							boolean esta = false;

							Unit u = ((DefaultGeodeticDatum) sour.getDatum()).getPrimeMeridian().getAngularUnit();

							units.add(u.toString());							

						}

					} catch (NoSuchAuthorityCodeException e) {

						res.add(String.valueOf(cod)+" --> "+e);

						// TODO Auto-generated catch block

						//e.printStackTrace();

					} catch (FactoryException e) {

						res.add(String.valueOf(cod)+" --> "+e);

						// TODO Auto-generated catch block

						//e.printStackTrace();

					}*/

	        }

	        /*for (int i = 0; i< res.size(); i++) {

	        	System.out.println(res.get(i));

	        }*/	

	        for (int i = 0; i< units.size(); i++) {

	        	System.out.println(units.get(i));

	        }

	        System.out.println("Numero CRSs validos: " + validos.size());

	        System.out.println("Numero medidas: " + units.size());

	        System.out.println("Numero CRSs fallidos: " + res.size());

		} catch (SQLException e1) {

			// TODO Auto-generated catch block

			e1.printStackTrace();

		}	

	}

}

package es.idr.crs.crsGT;

import java.util.Arrays;

import java.util.Iterator;

import org.geotools.factory.FactoryRegistry;

import org.geotools.referencing.operation.MathTransformProvider;

public class ProvidersTest {

	public static void main(String[] args) {

		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]);

				}

			}

		}

		providers = registry.getServiceProviders(MathTransformProvider.class);

		while (providers.hasNext()) {

            provider = (MathTransformProvider) providers.next();

            if (provider.nameMatches("IDR"))

            	System.out.println("-IDR-"+provider.toString());

            else

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

		}

	}

}

package es.idr.crs.crsGT;

import org.cresques.cts.IProjection;

import org.gvsig.crs.CrsException;

import org.gvsig.crs.CrsGT;

import org.gvsig.crs.ICrs;

import org.gvsig.crs.Proj4;

import org.gvsig.crs.repository.EpsgRepositoryGT;

import org.gvsig.crs.repository.EsriRepositoryGT;

import org.gvsig.crs.repository.Iau2000RepositoryGT;

import org.gvsig.crs.repository.UsrRepositoryGT;

public class CrsGTTest {

	public CrsGTTest() {

		epsgTest();

		esriTest();

		iau2000Test();

		usrTest();

	}

	/**

	 * @param args

	 */

	public static void main(String[] args) {

		CrsGTTest g = new CrsGTTest();

	}

	public void epsgTest() {

		EpsgRepositoryGT repo = new EpsgRepositoryGT();

		IProjection crs=null;

		crs = repo.getCrs("23030");

		if (crs!=null){

			System.out.println(crs.getAbrev());

			System.out.println(((ICrs)crs).getCode());

			System.out.println(((ICrs)crs).getWKT());

			System.out.println(((ICrs)crs).getDatum().getEIFlattening());

			System.out.println(((ICrs)crs).getDatum().getESemiMajorAxis());

			System.out.println(((ICrs)crs).isProjected()+"\n\n");

			System.out.println("----------- CrsWkt:--------------");

			System.out.println(((ICrs)crs).getCrsWkt().getAuthority()[0]+":"+((ICrs)crs).getCrsWkt().getAuthority()[1]);

			System.out.println(((ICrs)crs).getCrsWkt().getName());

			System.out.println(((ICrs)crs).getCrsWkt().getDatumName());

			System.out.println(((ICrs)crs).getCrsWkt().getGeogcs());

			System.out.println(((ICrs)crs).getCrsWkt().getProjcs());

			System.out.println(((ICrs)crs).getCrsWkt().getProjection());

			System.out.println(((ICrs)crs).getCrsWkt().getPrimen()[0]+" -- "+ ((ICrs)crs).getCrsWkt().getPrimen()[1]+"\n\n");

			try {

				Proj4 proj4 = new Proj4();

				System.out.println(proj4.exportToProj4(((CrsGT)crs).getCrsGT()));

			} catch (CrsException e) {

				// TODO Auto-generated catch block

				e.printStackTrace();

			}

		}

	}

	public void esriTest() {

		EsriRepositoryGT repo = new EsriRepositoryGT();

		IProjection crs=null;

		crs = repo.getCrs("23030");

		if (crs!=null){

			System.out.println(crs.getAbrev());

			System.out.println(((ICrs)crs).getCode());

			System.out.println(((ICrs)crs).getWKT());

			System.out.println(((ICrs)crs).getDatum().getEIFlattening());

			System.out.println(((ICrs)crs).getDatum().getESemiMajorAxis());

			System.out.println(((ICrs)crs).isProjected()+"\n\n");

			System.out.println("----------- CrsWkt:--------------");

			System.out.println(((ICrs)crs).getCrsWkt().getAuthority()[0]+":"+((ICrs)crs).getCrsWkt().getAuthority()[1]);

			System.out.println(((ICrs)crs).getCrsWkt().getName());

			System.out.println(((ICrs)crs).getCrsWkt().getDatumName());

			System.out.println(((ICrs)crs).getCrsWkt().getGeogcs());

			System.out.println(((ICrs)crs).getCrsWkt().getProjcs());

			System.out.println(((ICrs)crs).getCrsWkt().getProjection());

			System.out.println(((ICrs)crs).getCrsWkt().getPrimen()[0]+" -- "+ ((ICrs)crs).getCrsWkt().getPrimen()[1]+"\n\n");

			try {

				Proj4 proj4 = new Proj4();

				System.out.println(proj4.exportToProj4(((CrsGT)crs).getCrsGT()));

			} catch (CrsException e) {

				// TODO Auto-generated catch block

				e.printStackTrace();

			}

		}

	}

	public void iau2000Test() {

		Iau2000RepositoryGT repo = new Iau2000RepositoryGT();

		IProjection crs=null;

		crs = repo.getCrs("19900");

		if (crs!=null){

			System.out.println(crs.getAbrev());

			System.out.println(((ICrs)crs).getCode());

			System.out.println(((ICrs)crs).getWKT());

			System.out.println(((ICrs)crs).getDatum().getEIFlattening());

			System.out.println(((ICrs)crs).getDatum().getESemiMajorAxis());

			System.out.println(((ICrs)crs).isProjected()+"\n\n");

			System.out.println("----------- CrsWkt:--------------");

			System.out.println(((ICrs)crs).getCrsWkt().getAuthority()[0]+":"+((ICrs)crs).getCrsWkt().getAuthority()[1]);

			System.out.println(((ICrs)crs).getCrsWkt().getName());

			System.out.println(((ICrs)crs).getCrsWkt().getDatumName());

			System.out.println(((ICrs)crs).getCrsWkt().getGeogcs());

			System.out.println(((ICrs)crs).getCrsWkt().getProjcs());

			System.out.println(((ICrs)crs).getCrsWkt().getProjection());

			System.out.println(((ICrs)crs).getCrsWkt().getPrimen()[0]+" -- "+ ((ICrs)crs).getCrsWkt().getPrimen()[1]+"\n\n");

			try {

				Proj4 proj4 = new Proj4();

				System.out.println(proj4.exportToProj4(((CrsGT)crs).getCrsGT()));

			} catch (CrsException e) {

				// TODO Auto-generated catch block

				e.printStackTrace();

			}

		}

	}

	public void usrTest() {

		UsrRepositoryGT repo = new UsrRepositoryGT();

		IProjection crs=null;

		crs = repo.getCrs("2000");

		if (crs!=null){

			System.out.println(crs.getAbrev());

			System.out.println(((ICrs)crs).getCode());

			System.out.println(((ICrs)crs).getWKT());

			System.out.println(((ICrs)crs).getDatum().getEIFlattening());

			System.out.println(((ICrs)crs).getDatum().getESemiMajorAxis());

			System.out.println(((ICrs)crs).isProjected()+"\n\n");

			System.out.println("----------- CrsWkt:--------------");

			System.out.println(((ICrs)crs).getCrsWkt().getAuthority()[0]+":"+((ICrs)crs).getCrsWkt().getAuthority()[1]);

			System.out.println(((ICrs)crs).getCrsWkt().getName());

			System.out.println(((ICrs)crs).getCrsWkt().getDatumName());

			System.out.println(((ICrs)crs).getCrsWkt().getGeogcs());

			System.out.println(((ICrs)crs).getCrsWkt().getProjcs());

			System.out.println(((ICrs)crs).getCrsWkt().getProjection());

			System.out.println(((ICrs)crs).getCrsWkt().getPrimen()[0]+" -- "+ ((ICrs)crs).getCrsWkt().getPrimen()[1]+"\n\n");

			try {

				Proj4 proj4 = new Proj4();

				System.out.println(proj4.exportToProj4(((CrsGT)crs).getCrsGT()));

			} catch (CrsException e) {

				// TODO Auto-generated catch block

				e.printStackTrace();

			}

		}

	}

}

package org.gvsig.crs.jscrollpanel;

import java.awt.Dimension;

import javax.swing.BorderFactory;

import javax.swing.JFrame;

import javax.swing.JScrollPane;

import javax.swing.JTable;

import javax.swing.table.DefaultTableModel;

import javax.swing.table.TableColumn;

public class ScrollPanelTest {

	JTable jTable = null;

	JScrollPane jScrollPane = null;

	public ScrollPanelTest() {

		super();

		JFrame jFrame = new JFrame();

		jFrame.getContentPane().add(getJScrollPane());	

		jFrame.setSize(500, 150);

		jFrame.setVisible(true);

		// TODO Auto-generated constructor stub

	}

	public JTable getJTable() {

		if (jTable == null) {

			String[] columnNames= {"UNA COLUMNA", "DOS COLUMNAS"};

			Object[][]data = {};			

			DefaultTableModel dtm = new DefaultTableModel(data, columnNames);

			jTable = new JTable(dtm);

			jTable.setAutoResizeMode(JTable.AUTO_RESIZE_OFF);

			TableColumn column = null;

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

			       column = jTable.getColumnModel().getColumn(i);

			       if (i == 0) {

			           column.setPreferredWidth(600);

			       }

			   }

			//jTable.setPreferredSize(new Dimension(800, 100));			

		}

		return jTable;

	}

	private JScrollPane getJScrollPane() {

		if (jScrollPane == null) {

			jScrollPane = new JScrollPane(getJTable(),JScrollPane.VERTICAL_SCROLLBAR_AS_NEEDED, JScrollPane.HORIZONTAL_SCROLLBAR_ALWAYS);

			jScrollPane.setPreferredSize(new Dimension(500,150));

			jScrollPane.setBorder(

				    BorderFactory.createCompoundBorder(

					BorderFactory.createCompoundBorder(

							BorderFactory.createTitledBorder("TABLA GEN?RICA DE 2 CAMPOS"),

							BorderFactory.createEmptyBorder(5,5,5,5)),

							jScrollPane.getBorder()));

			jScrollPane.setViewportView(getJTable());

		}

		return jScrollPane;

	}

	public static void main(String[] args){

		ScrollPanelTest al = new ScrollPanelTest();

	}

}

package org.gvsig.crs.repository;

import org.gvsig.crs.Crs;

public class RepositoryTest {

	public static void main(String[] args){

		EpsgRepository epsgRep= new EpsgRepository();

		Crs crs = (Crs)epsgRep.getCrs("23030");

		if (crs != null){

			String abrev = crs.getAbrev();

			System.out.println(abrev);

		}

		else

			System.out.println("el c?digo no existe en el repositorio.");

	}

}

- COpetation.getInverted() no est? teniendo en cuenta si hay par?metros de transformaci?n -HECHO-

- Resolver el tema del cacheado de CRSs en los repositorios.

- Por qu? no est?n implementados  getPOrig() y getPDest() de COperation? -HECHO-

- Plantearse la necesidad de ICOperation:

	- COperation ya no la implementa. La conservo preventivamente. -HECHO-

- Crs y CrsGT: No construir los objetos CrsProj en el constructor, hacerlo la primera vez que se pidan.

- Habr?a que reducir el tiempo de arranque de la extensi?n.

- Depuraciones para detectar ineficiencias.

	- COperation.operate se llama para cada punto: optimizar su contenido (creaci?n de objetos). -HECHO-

	- Hacer que CrsGT.getCT no construya nuevos CrsProj cada vez que se llama.

- libJCrs y windows?????????????'

- ICrs: mirar qu? m?todos no tienen por qu? estar en la interfaz.

- Dependencias de gvSIG en libJCrs:

	- COperation.convert(...) utiliza: NotificationManager.addError(e) de andami;

	- Persistencia: depende de PluginServices de andami. -> subir la persistencia a la extension. -HECHO-

	- dependencia de CMS: IProjection, ICoordTrans, IDatum, ICRSFactory, org.cresques.geo.ViewPortData.

- Crs.getScale(): double invF = getDatum().getEIFlattening(); if (invF == Double.POSITIVE_INFINITY){}

	?Qu? valor da getDatum().getEIFlattening() en el caso de una esfera? geoltools da Double.POSITIVE_INFINITY.

- Busqueda IAU por zona: en mayusculas no encuentra. -HECHO-

- Persistencia de transformaciones en F2. -HECHO-

- Persistencia en el proyecto: Implementar IProyectio.getFullCode() y lo que esto implica. -HECHO-

- OGR y WINDOWS....!!!!???? 

- Revisar classpath de extJCRS (?log4j-1.2.8.jar?,etc...)

- Cuando CRS fuente y destino son iguales -> no operar, devolver el mismo punto.

BUGS:

- Cargas un shp 23030 diciendole que es 4258 con transf. EPSG (7) y salta una excepcci?n de proj4: controlar o subir estos

	errores (es un error como el que silenciamos de nadgrids).

	En este caso el error es debido a que las coordenadas se salen del rango del CRS elegido. Esta situaci?n ArcMap la maneja desavilitando

	la capa y dejando de visualizarla (algo parecedo a cuando gvSIG le pone la x de erroneo a una capa).

- Cuando hay un fallo al elegir el CRS de la vista y despu?s guardas el proyecto, parece ser que se guarda como null (o algo as?)

	y falla al abrirlo. Posibles soluciones: en esos casos guardar el CRS por defecto.

- Nomenclator: Cuando pulsas localizar, debe dibuajar el punto en la vista y hacer zoom a el. Lo dibuja pero aparece una

	excepci?n y no hace el zoom.

org.geotools.referencing.operation.projection.IdrAitoff$Provider

org.geotools.referencing.operation.projection.IdrAzimuthalEquidistant$Provider

org.geotools.referencing.operation.projection.IdrAzimuthalEquidistant$Provider_Modified

org.geotools.referencing.operation.projection.IdrAzimuthalEquidistant$Provider_Guam

org.geotools.referencing.operation.projection.IdrBonne$Provider

org.geotools.referencing.operation.projection.IdrBonne$Provider_SouthOrientated

org.geotools.referencing.operation.projection.IdrCassiniSoldner$Provider

org.geotools.referencing.operation.projection.IdrCassiniSoldner$Provider_Hyperbolic

org.geotools.referencing.operation.projection.IdrCrasterParabolic$Provider

org.geotools.referencing.operation.projection.IdrCylindricalEqualArea$Provider

org.geotools.referencing.operation.projection.IdrEckertI$Provider

org.geotools.referencing.operation.projection.IdrEckertII$Provider

org.geotools.referencing.operation.projection.IdrEckertIII$Provider

org.geotools.referencing.operation.projection.IdrEckertIV$Provider

org.geotools.referencing.operation.projection.IdrEckertV$Provider

org.geotools.referencing.operation.projection.IdrEckertVI$Provider

org.geotools.referencing.operation.projection.IdrEquidistantConic$Provider

org.geotools.referencing.operation.projection.IdrEquidistantCylindrical$Provider

org.geotools.referencing.operation.projection.IdrEquirectangular$Provider

org.geotools.referencing.operation.projection.IdrFlatPolarQuartic$Provider

org.geotools.referencing.operation.projection.IdrGallStereographic$Provider

org.geotools.referencing.operation.projection.IdrGnomonic$Provider

org.geotools.referencing.operation.projection.IdrGoode$Provider

org.geotools.referencing.operation.projection.IdrHammerAitoff$Provider

org.geotools.referencing.operation.projection.IdrKrovak$Provider

org.geotools.referencing.operation.projection.IdrLabordeMadagascar$Provider

org.geotools.referencing.operation.projection.IdrLambertAzimuthalEqualArea$Provider

org.geotools.referencing.operation.projection.IdrLambertConformalWestOrientated$Provider

org.geotools.referencing.operation.projection.IdrLambertConicNearConformal$Provider

org.geotools.referencing.operation.projection.IdrLoximuthal$Provider

org.geotools.referencing.operation.projection.IdrMercator$Provider1SP

org.geotools.referencing.operation.projection.IdrMillerCylindrical$Provider

org.geotools.referencing.operation.projection.IdrMollweide$Provider

org.geotools.referencing.operation.projection.IdrNearSidedPerspective$Provider

org.geotools.referencing.operation.projection.IdrNewZealandMapGrid$Provider

org.geotools.referencing.operation.projection.IdrObliqueMercator$Provider

org.geotools.referencing.operation.projection.IdrObliqueMercator$Provider_Hotine

org.geotools.referencing.operation.projection.IdrObliqueMercator$Provider_TwoPoint

org.geotools.referencing.operation.projection.IdrObliqueMercator$Provider_Hotine_TwoPoint

org.geotools.referencing.operation.projection.IdrOrthographic$Provider

org.geotools.referencing.operation.projection.IdrPlateCarree$Provider

org.geotools.referencing.operation.projection.IdrPolyconic$Provider

org.geotools.referencing.operation.projection.IdrQuarticAuthalic$Provider

org.geotools.referencing.operation.projection.IdrRobinson$Provider

org.geotools.referencing.operation.projection.IdrSinusoidal$Provider

org.geotools.referencing.operation.projection.IdrTransverseMercator$Provider

org.geotools.referencing.operation.projection.IdrTransverseMercator$Provider_SouthOrientated

org.geotools.referencing.operation.projection.IdrTransverseMercator$Provider_TransverseMercatorZonedGridSystem

org.geotools.referencing.operation.projection.IdrTunisiaMiningGrid$Provider

org.geotools.referencing.operation.projection.IdrVanDerGrinten$Provider

org.geotools.referencing.operation.projection.IdrWinkelI$Provider

org.geotools.referencing.operation.projection.IdrWinkelII$Provider

org.geotools.referencing.operation.projection.IdrWinkelTripel$Provider

org.gvsig.crs.JCRSLibrary

es.idr.teledeteccion.connection.esri.FactoryUsingHSQL

es.idr.teledeteccion.connection.usr.FactoryUsingHSQL

es.idr.teledeteccion.connection.iau2000.FactoryUsingHSQL

/*

 *    GeoTools - OpenSource mapping toolkit

 *    http://geotools.org

 *    (C) 2005-2006, GeoTools Project Managment Committee (PMC)

 * 

 *    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 es.idr.teledeteccion.connection.iau2000;

// J2SE dependencies

import java.io.IOException;

import java.sql.Connection;

import java.sql.PreparedStatement;

import java.sql.ResultSet;

import java.sql.SQLException;

import java.sql.Statement;

import java.text.ParseException;

import java.util.Hashtable;

import java.util.Set;

import org.geotools.metadata.iso.citation.CitationImpl;

import org.geotools.referencing.CRS;

import org.geotools.referencing.wkt.Parser;

import org.gvsig.crs.CrsFactory;

import org.hsqldb.jdbc.jdbcDataSource;

import org.opengis.metadata.citation.Citation;

import org.opengis.referencing.FactoryException;

import org.opengis.referencing.IdentifiedObject;

import org.opengis.referencing.crs.CRSAuthorityFactory;

import org.opengis.referencing.crs.CRSFactory;

import org.opengis.referencing.crs.CompoundCRS;

import org.opengis.referencing.crs.CoordinateReferenceSystem;

import org.opengis.referencing.crs.DerivedCRS;

import org.opengis.referencing.crs.EngineeringCRS;

import org.opengis.referencing.crs.GeocentricCRS;

import org.opengis.referencing.crs.GeographicCRS;

import org.opengis.referencing.crs.ImageCRS;

import org.opengis.referencing.crs.ProjectedCRS;

import org.opengis.referencing.crs.TemporalCRS;

import org.opengis.referencing.crs.VerticalCRS;

import org.opengis.util.InternationalString;

public class FactoryUsingHSQL extends jdbcDataSource implements CRSAuthorityFactory {

    /**

     * The connection to the EPSG database.

     */

    protected Connection connection;

    //  object factory

    protected CRSFactory crsFactory;

    /** Cache of parsed CoordinateReferenceSystem WKT by IAU2000_NUMBER */

    private Hashtable cache = new Hashtable();

    /**

     * The database name.

     *

     * @since 2.3

     */

    public static final String DATABASE_NAME = "IAU2000";

    public FactoryUsingHSQL() {

    	 try {

             loadDefault();

         }

         catch( IOException oops ){

             System.err.println("Could not load "+ oops );

         }

    	//this(FactoryFinder.getCRSFactory(null));

    }   

	protected FactoryUsingHSQL(CRSFactory factory) {

		this.crsFactory = factory;

        try {

            loadDefault();

        }

        catch( IOException oops ){

            System.err.println("Could not load "+ oops );

        }

	}

	protected void loadDefault() throws IOException {

		String db = "jdbc:hsqldb:file:" + CrsFactory.getDataBaseFolder().getAbsolutePath() + "/iau2000";

		setDatabase(db);

		setUser("sa");

		try {

			connection = super.getConnection();

		} catch (SQLException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

    }

	public CoordinateReferenceSystem createCoordinateReferenceSystem(String code) throws FactoryException {

		// TODO Auto-generated method stub

		CoordinateReferenceSystem returnValue = null;

		String wkt = null;

		String IAU2000_NUMBER = code.substring(code.indexOf(":")+1, code.length());

		int cod = Integer.parseInt(IAU2000_NUMBER);

		if( cache.contains( IAU2000_NUMBER ) ){

            Object value = cache.get( IAU2000_NUMBER );

            if( value instanceof Throwable ){

                throw new FactoryException( "WKT for "+code+" could not be parsed", (Throwable) value );

            }

            if( value instanceof CoordinateReferenceSystem){

                return (CoordinateReferenceSystem) value;

            }            

        }

		try {        	

        	final PreparedStatement stmt;

            /*stmt = connection.prepareStatement("SELECT IAU_WKT"

                                          + " FROM IAU2000"

                                          + " WHERE IAU_CODE = " + code.substring(code.indexOf(":")+1, code.length()));

            ResultSet result = stmt.executeQuery();

            */

            String sentence = "SELECT IAU_WKT"

            + " FROM IAU2000"

            + " WHERE IAU_CODE = " + code.substring(code.indexOf(":")+1, code.length());

        	Statement st = null;

    		ResultSet result = null;

        	st = connection.createStatement();

			result = st.executeQuery(sentence);

			st.close();        	

        	while (result.next()) {

                wkt    = result.getString("iau_wkt"); 	

            }

            result.close();

            if( wkt.indexOf( cod ) == -1){

            	wkt = wkt.trim();

            	wkt = wkt.substring(0, wkt.length()-1 );

            	wkt += ",AUTHORITY[\"IAU2000\",\""+cod+"\"]]";

            }

        } catch (SQLException exception) {

            System.out.println(exception);

        }

        Parser parser = new Parser();

		try {

			returnValue = parser.parseCoordinateReferenceSystem(wkt);

		} catch (ParseException e) {

			// TODO Auto-generated catch block

			System.out.println("Cadena WKT no ha podido ser parseada");

			e.printStackTrace();

		}

        return returnValue;		

	}

	public CompoundCRS createCompoundCRS(String arg0) throws FactoryException {

		// TODO Auto-generated method stub

		return null;

	}

	public DerivedCRS createDerivedCRS(String arg0) throws FactoryException {

		// TODO Auto-generated method stub

		return null;

	}

	public EngineeringCRS createEngineeringCRS(String arg0) throws FactoryException {

		// TODO Auto-generated method stub

		return null;

	}

	public GeographicCRS createGeographicCRS(String arg0) throws FactoryException {

		// TODO Auto-generated method stub

		return null;

	}

	public GeocentricCRS createGeocentricCRS(String arg0) throws FactoryException {

		// TODO Auto-generated method stub

		return null;

	}

	public ImageCRS createImageCRS(String arg0) throws FactoryException {

		// TODO Auto-generated method stub

		return null;

	}

	public ProjectedCRS createProjectedCRS(String arg0) throws FactoryException {

		// TODO Auto-generated method stub

		return null;

	}

	public TemporalCRS createTemporalCRS(String arg0) throws FactoryException {

		// TODO Auto-generated method stub

		return null;

	}

	public VerticalCRS createVerticalCRS(String arg0) throws FactoryException {

		// TODO Auto-generated method stub

		return null;

	}

	public Citation getAuthority() {		

		Citation citation = CitationImpl.createCitation("IAU2000");

		citation.getIdentifiers().add("IAU2000");

		return citation;

	}

	public Set getAuthorityCodes(Class arg0) throws FactoryException {

		// TODO Auto-generated method stub

		return null;

	}

	public InternationalString getDescriptionText(String arg0) throws FactoryException {

		// TODO Auto-generated method stub

		return null;

	}

	public IdentifiedObject createObject(String code) throws FactoryException {

		return createCoordinateReferenceSystem(code);

	}

	public Citation getVendor() {

		return CitationImpl.createCitation("IAU2000");

	}

}

/*

 *    GeoTools - OpenSource mapping toolkit

 *    http://geotools.org

 *    (C) 2005-2006, GeoTools Project Managment Committee (PMC)