Application: gvSIG desktop

bin

bin-test

package org.cresques.cts;

import junit.framework.Test;

import junit.framework.TestSuite;

public class AllTests {

	public static Test suite() {

		TestSuite suite = new TestSuite("Test for org.cresques.cts");

		//$JUnit-BEGIN$

		//$JUnit-END$

		return suite;

	}

}

#Translations for language [fr]

#Mon Oct 30 09:38:21 CET 2006

Aceptar=Accepter

Aplicar=Appliquer

Block_Size_=

Cancelar=Annuler

Compression_=

datum=

Generate_Tfw_=

Interleave_=

Photometric_=

Progressive_=

projection=

reference_system=

zone=

#Translations for language [de]

#Mon Oct 30 09:38:21 CET 2006

Aceptar=OK

Aplicar=Anwenden

Block_Size_=

Cancelar=Abbrechen

Compression_=

datum=

Generate_Tfw_=

Interleave_=

Photometric_=

Progressive_=

projection=

reference_system=

zone=

#Translations for language [eu]

#Mon Oct 30 09:38:21 CET 2006

Aceptar=Ados

Aplicar=Aplikatu

Block_Size_=

Cancelar=Utzi

Compression_=

datum=

Generate_Tfw_=

Interleave_=

Photometric_=

Progressive_=

projection=

reference_system=

zone=

#Translations for language [it]

#Tue Nov 07 12:30:01 CET 2006

Aceptar=Accetta

Aplicar=Applica

Block_Size_=Dimensioni blocco\:

Cancelar=Cancella

Compression_=Compressione

datum=Datum

Generate_Tfw_=Genera Tfw\:

Interleave_=Interlacciato\:

Photometric_=Fotometrica\:

Progressive_=Progressivo\:

projection=Proiezione

reference_system=Sistema di riferimento

zone=Fuso

#Translations for language [es]

#Mon Oct 30 09:38:21 CET 2006

Aceptar=Aceptar

Aplicar=Aplicar

Block_Size_=Tama\u00f1o de bloque\:

Cancelar=Cancelar

Compression_=Compresi\u00f3n\:

datum=Datum

Generate_Tfw_=Generar Tfw\:

Interleave_=Entrelazado\:

Photometric_=Fotom\u00e9trica\:

Progressive_=Progresivo\:

projection=Proyecci\u00f3n

reference_system=Sistema de Referencia

zone=Huso

#Translations for language [en]

#Mon Oct 30 09:38:21 CET 2006

Aceptar=Accept

Aplicar=Apply

Block_Size_=Block Size\:

Cancelar=Cancel

Compression_=Compression\:

datum=Datum

Generate_Tfw_=Generate Tfw\:

Interleave_=Interleave\:

Photometric_=Photometric\:

Progressive_=Progressive\:

projection=Projection

reference_system=Reference System

zone=Zone

#Translations for language [gl]

#Mon Oct 30 09:38:21 CET 2006

Aceptar=Aceptar

Aplicar=Aplicar

Block_Size_=

Cancelar=Cancelar

Compression_=

datum=

Generate_Tfw_=

Interleave_=

Photometric_=

Progressive_=

projection=

reference_system=

zone=

#Translations for language [ca]

#Mon Oct 30 09:38:21 CET 2006

Aceptar=Acceptar

Aplicar=Aplicar

Block_Size_=

Cancelar=Cancel\u00b7lar

Compression_=

datum=

Generate_Tfw_=

Interleave_=

Photometric_=

Progressive_=

projection=Projecci\u00f3

reference_system=Sistema de refer\u00e8ncia

zone=

#Translations for language [pt]

#Mon Oct 30 09:38:21 CET 2006

Aceptar=Aceitar

Aplicar=Aplicar

Block_Size_=

Cancelar=Cancelar

Compression_=

datum=

Generate_Tfw_=

Interleave_=

Photometric_=

Progressive_=

projection=

reference_system=

zone=

#Translations for language [cs]

#Mon Nov 06 09:04:31 CET 2006

Aceptar=Budi\u017e

Aplicar=Pou\u017e\u00edt

Block_Size_=

Cancelar=Zru\u0161it

Compression_=

datum=Datum

Generate_Tfw_=

Interleave_=

Photometric_=

Progressive_=

projection=Projekce

reference_system=Prostorov\u00fd referen\u010dn\u00ed syst\u00e9m

zone=Z\u00f3na

<?xml version="1.0" encoding="UTF-8" ?>

<!DOCTYPE log4j:configuration SYSTEM "log4j.dtd">

<log4j:configuration xmlns:log4j="http://jakarta.apache.org/log4j/">

	<appender name="CONSOLE" class="org.apache.log4j.ConsoleAppender">

		<layout class="org.apache.log4j.PatternLayout">

			<param name="ConversionPattern" value="%d{HH:mm:ss,SSS} %-5p [%c{2}.%M()]\n  %m%n" />

		</layout>

	</appender>

	<category name="org.gvsig.tools">

		<priority value="DEBUG" />

	</category>

	<category name="org.cresques">

		<priority value="DEBUG" /> 

	</category>

	<category name="org.gvsig.fmap.crs">

		<priority value="DEBUG" /> 

	</category>

	<root>

		<priority value="INFO" />

		<appender-ref ref="CONSOLE" />

	</root>

</log4j:configuration>

/* 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.

 * 

 */

/*

 * AUTHORS (In addition to CIT):

 * 2009 {DiSiD Technologies}  {{Task}}

 */

package org.gvsig.fmap.crs.persistence;

import org.cresques.cts.IProjection;

import org.gvsig.fmap.crs.CRSFactory;

import org.gvsig.tools.dataTypes.DataTypes;

import org.gvsig.tools.dynobject.DynStruct;

import org.gvsig.tools.persistence.AbstractSinglePersistenceFactory;

import org.gvsig.tools.persistence.PersistenceManager;

import org.gvsig.tools.persistence.PersistentState;

import org.gvsig.tools.persistence.exception.PersistenceException;

/**

 * Factory to persist {@link IProjection} objects. The information about the

 * {@link IProjection} which will be persisted is its full code, containing all

 * required information to be able to reconstruct the {@link IProjection} object

 * through the {@link CRSFactory#getCRS(String)} method.

 * <p>

 * <strong>NOTE:</strong>To activate this factory, it must be instanced and

 * registered in the {@link PersistenceManager}. This will be usually performed

 * from the project Library.

 * </p>

 * 

 * @author 2009- <a href="cordinyana@gvsig.org">C?sar Ordi?ana</a> - gvSIG team

 * @author 2010- <a href="jjdelcerro@gvsig.org">Joaquin Jose del Cerro</a> -

 *         gvSIG team

 */

public class ProjectionPersistenceFactory extends

		AbstractSinglePersistenceFactory {

	private static final String FIELD_FULLCODE = "fullcode";

	private static final String DESCRIPTION_FULLCODE = "Projection abbreviature";

	private static final String DYNCLASS_NAME = "Projection";

	private static final String DYNCLASS_DESCRIPTION = "Projection";

	/**

	 * Creates a new {@link ProjectionPersistenceFactory}.

	 */

	public ProjectionPersistenceFactory() {

		super(IProjection.class, DYNCLASS_NAME, DYNCLASS_DESCRIPTION, null,

				null);

		DynStruct definition = this.getDefinition();

		definition.addDynField(FIELD_FULLCODE).setDescription(

				DESCRIPTION_FULLCODE).setMandatory(true).setType(

				DataTypes.STRING);

	}

	public Object createFromState(PersistentState state)

			throws PersistenceException {

		return CRSFactory.getCRS(state.getString(FIELD_FULLCODE));

	}

	public void saveToState(PersistentState state, Object obj)

			throws PersistenceException {

		state.set(FIELD_FULLCODE, ((IProjection) obj).getFullCode());

	}

}

/* 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.fmap.crs.persistence;

import org.cresques.DataTypes;

import org.cresques.cts.ICoordTrans;

import org.cresques.cts.IProjection;

import org.gvsig.tools.dynobject.DynStruct;

import org.gvsig.tools.persistence.AbstractSinglePersistenceFactory;

import org.gvsig.tools.persistence.PersistentState;

import org.gvsig.tools.persistence.exception.PersistenceException;

/**

 * @author gvSIG Team

 * @version $Id$

 *

 */

public class CoordTransPersistenceFactory extends AbstractSinglePersistenceFactory {

    private static final String FIELD_SOURCE_PROJECTION = "source";

    private static final String DESCRIPTION_SOURCE_PROJECTION  = "Source projection";

    private static final String FIELD_DESTINATION_PROJECTION = "destination";

    private static final String DESCRIPTION_DESTINATION_PROJECTION = "Destination Projection";

    private static final String DYNCLASS_NAME = "CoordTrans";

    private static final String DYNCLASS_DESCRIPTION = "Coordinates Transformation";

    /**

     * Creates a new {@link CoordTransPersistenceFactory}.

     */

    public CoordTransPersistenceFactory() {

        super(ICoordTrans.class, DYNCLASS_NAME, DYNCLASS_DESCRIPTION, null, null);

        DynStruct definition = this.getDefinition();

        definition.addDynField(FIELD_SOURCE_PROJECTION).setDescription(

            DESCRIPTION_SOURCE_PROJECTION).setMandatory(true).setType(

                DataTypes.CRS);

        definition.addDynField(FIELD_DESTINATION_PROJECTION).setDescription(

            DESCRIPTION_DESTINATION_PROJECTION).setMandatory(true).setType(

                DataTypes.CRS);

    }

    public Object createFromState(PersistentState state)

            throws PersistenceException {

        IProjection sourceProjection = (IProjection)state.get(FIELD_SOURCE_PROJECTION);

        IProjection destinationProjection = (IProjection)state.get(FIELD_DESTINATION_PROJECTION);      

        return sourceProjection.getCT(destinationProjection);

    }

    public void saveToState(PersistentState state, Object obj)

            throws PersistenceException {

        state.set(FIELD_SOURCE_PROJECTION, ((ICoordTrans) obj).getPOrig());

        state.set(FIELD_DESTINATION_PROJECTION, ((ICoordTrans) obj).getPDest());

    }

}

package org.gvsig.fmap.crs;

import org.cresques.cts.ICRSFactory;

import org.cresques.cts.IProjection;

/**

 * Fabrica de CRS.

 * Centraliza las peticiones de creaci?n de objetos CRS de todo fmap.

 * @author luisw

 *

 */

public class CRSFactory {

	public static ICRSFactory cp = null;

	public static IProjection getCRS(String code) {

		return cp.get(code);

	}

}

/* 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.cresques;

/**

 * @author gvSIG Team

 * @version $Id$

 * 

 */

public interface DataTypes extends org.gvsig.tools.dataTypes.DataTypes {

    /*

	 * 

	 */

    public static final int CRS = OBJECT + 1;

}

/* 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.cresques.coerce;

import org.cresques.cts.IProjection;

import org.gvsig.tools.dataTypes.CoercionException;

import org.gvsig.tools.dataTypes.DataTypesManager.Coercion;

/**

 * Convert a Projection to String.

 * 

 * Support convert:

 * - Projection to String (do nothing)

 * 

 * @author gvSIG Team

 * @version $Id$

 * 

 */

public class CoerceToString implements Coercion {

	Coercion previous = null;

	public CoerceToString() {

		// Do nothing

	}

	public CoerceToString(Coercion previous) {

		this.previous = previous;

	}

	public Object coerce(Object value) throws CoercionException {

		try {

			if( value == null || value instanceof String ) {

				return value;

			}

			if( value instanceof IProjection ) {

				return ((IProjection)value).getFullCode();

			}

			if( previous != null ) {

				return previous.coerce(value);

			}

		} catch (Exception e) {

			throw new CoercionException(e);

		}

		throw new CoercionException();

	}

}

/* 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.cresques.coerce;

import org.cresques.cts.IProjection;

import org.gvsig.fmap.crs.CRSFactory;

import org.gvsig.tools.dataTypes.CoercionException;

import org.gvsig.tools.dataTypes.DataTypesManager.Coercion;

/**

 * Convert a string value of projection code to IProjection

 * 

 * @author gvSIG Team

 * @version $Id$

 * 

 */

public class CoerceToCRS implements Coercion {

	public Object coerce(Object value) throws CoercionException {

		try {

			if( value == null || value instanceof IProjection) {

				return value;

			}

			Object proj = CRSFactory.getCRS(value.toString());

			if( proj == null ) {

				throw new CoercionException("Can't convert value '"+value.toString()+"' to Projection.");

			}

			return proj;

		} catch (Exception e) {

			throw new CoercionException(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.cresques;

import org.cresques.coerce.CoerceToCRS;

import org.cresques.coerce.CoerceToString;

import org.cresques.cts.IProjection;

import org.gvsig.fmap.crs.CRSFactory;

import org.gvsig.fmap.crs.persistence.CoordTransPersistenceFactory;

import org.gvsig.fmap.crs.persistence.ProjectionPersistenceFactory;

import org.gvsig.tools.ToolsLibrary;

import org.gvsig.tools.ToolsLocator;

import org.gvsig.tools.dataTypes.DataTypesManager;

import org.gvsig.tools.library.AbstractLibrary;

import org.gvsig.tools.library.LibraryException;

public class ProjectionLibrary extends AbstractLibrary {

    @Override

    public void doRegistration() {

        registerAsAPI(ProjectionLibrary.class);

        require(ToolsLibrary.class);

    }

	protected void doInitialize() throws LibraryException {

	}

	protected void doPostInitialize() throws LibraryException {

		if (CRSFactory.cp==null) {

			throw new CRSFactoryNotRegisteredException();

		}

		DataTypesManager dataTypesManager = ToolsLocator.getDataTypesManager();

        dataTypesManager.addtype(DataTypes.CRS, "CRS", "CRS",

            IProjection.class, new CoerceToCRS());

        dataTypesManager.setCoercion(DataTypes.STRING, new CoerceToString(

            dataTypesManager.getCoercion(DataTypes.STRING)));

		// Register the PersistenceFactory to be able to persist 

		// IProjection objects

		ToolsLocator.getPersistenceManager().registerFactory(

				new ProjectionPersistenceFactory());

		ToolsLocator.getPersistenceManager().registerFactory(

            new CoordTransPersistenceFactory());

	}

}

/* 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.cresques.impl;

import org.cresques.ProjectionLibrary;

import org.cresques.impl.cts.ProjectionPool;

import org.gvsig.fmap.crs.CRSFactory;

import org.gvsig.tools.library.AbstractLibrary;

import org.gvsig.tools.library.LibraryException;

public class CresquesCtsLibrary extends AbstractLibrary {

    @Override

    public void doRegistration() {

        registerAsImplementationOf(ProjectionLibrary.class);

    }

    protected void doInitialize() throws LibraryException {

        CRSFactory.cp = new ProjectionPool();

    }

    protected void doPostInitialize() throws LibraryException {

        // Nothing to do

    }

}

/*

 * Cresques Mapping Suite. Graphic Library for constructing mapping applications.

 *

 * Copyright (C) 2004-7.

 *

 * 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:

 *

 * cresques@gmail.com

 */

package org.cresques.impl.cts.gt2;

import org.cresques.cts.UTM;

import org.geotools.cs.AxisInfo;

import org.geotools.cs.Projection;

import org.geotools.units.Unit;

import org.gvsig.fmap.crs.CRSFactory;

import org.opengis.referencing.FactoryException;

/**

 * @author "Luis W. Sevilla" <sevilla_lui@gva.es>

 */

public class CSUTM extends CoordSys implements UTM {

    public CSUTM(CSDatum datum, int zone) {

        super(datum);

        init(datum, zone, "N");

    }

    public CSUTM(CSDatum datum, int zone, String ns) {

        super(datum);

        init(datum, zone, ns);

    }

    public void init(CSDatum datum, int zone, String ns) {

        Unit linearUnit = Unit.METRE;

        javax.media.jai.ParameterList params = csFactory.createProjectionParameterList("Transverse_Mercator");

        params.setParameter("semi_major",

                            datum.getDatum().getEllipsoid().getSemiMajorAxis());

        params.setParameter("semi_minor",

                            datum.getDatum().getEllipsoid().getSemiMinorAxis());

        params.setParameter("central_meridian", (double) ((zone * 6) - 183));

        params.setParameter("latitude_of_origin", 0.0);

        params.setParameter("scale_factor", 0.9996);

        params.setParameter("false_easting", 500000.0);

        if (ns.toUpperCase().compareTo("S") == 0)

        	params.setParameter("false_northing", 10000000.0);

        else

        	params.setParameter("false_northing", 0.0);

        try {

            Projection projection = csFactory.createProjection("UTM" + zone,

                                                               "Transverse_Mercator",

                                                               params);

            projCS = csFactory.createProjectedCoordinateSystem(projection.getName()

                                                                         .toString(),

                                                               geogCS,

                                                               projection,

                                                               linearUnit,

                                                               AxisInfo.X,

                                                               AxisInfo.Y);

        } catch (FactoryException e) {

            // TODO Bloque catch generado autom?ticamente

            e.printStackTrace();

        }

    }

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

        double scale = super.getScale(minX, maxX, w, dpi);

        if (projCS != null) { // Es geogr?fico; calcula la escala.

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

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

                    w; // pixels

        }

        return scale;

    }

    public String toString() {

        return projCS.toString();

    }

    public Object clone() throws CloneNotSupportedException {

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

     }

}

/*

 * Cresques Mapping Suite. Graphic Library for constructing mapping applications.

 *

 * Copyright (C) 2004-5.

 *

 * 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:

 *

 * cresques@gmail.com

 */

package org.cresques.impl.cts.gt2;

import org.geotools.cs.AxisInfo;

import org.geotools.cs.GeographicCoordinateSystem;

import org.geotools.cs.Projection;

import org.geotools.measure.Angle;

import org.geotools.units.Unit;

import org.opengis.referencing.FactoryException;

/**

 * @author "Luis W. Sevilla" <sevilla_lui@gva.es>

 *

 */

public class CSGaussPt extends CoordSys {

    public static CSGaussPt hgd73 = new CSGaussPt(CSDatum.d73);

    public CSGaussPt(CSDatum datum) {

        super(datum);

        geogCS = new GeographicCoordinateSystem(datum.getName(null),

                                                datum.getDatum());

        Unit linearUnit = Unit.METRE;

        javax.media.jai.ParameterList params = csFactory.createProjectionParameterList("Oblique_Stereographic");

        params.setParameter("semi_major",

                            datum.getDatum().getEllipsoid().getSemiMajorAxis());

        params.setParameter("semi_minor",

                            datum.getDatum().getEllipsoid().getSemiMinorAxis());

        params.setParameter("central_meridian",

                            (new Angle("8?07'54.862\"W")).degrees());

        //params.setParameter("longitude_of_origin", (new Angle("8?07'54.862\"")).degrees());

        params.setParameter("latitude_of_origin",

                            (new Angle("39?40'N")).degrees());

        params.setParameter("scale_factor", 1.0);

        params.setParameter("false_easting", 0D); //180598D);

        params.setParameter("false_northing", 0D); //-86990D);

        try {

            Projection projection = csFactory.createProjection("GAUSS",

                                                               "Oblique_Stereographic",

                                                               params);

            projCS = csFactory.createProjectedCoordinateSystem(projection.getName()

                                                                         .toString(),

                                                               geogCS,

                                                               projection,

                                                               linearUnit,

                                                               AxisInfo.X,

                                                               AxisInfo.Y);

        } catch (FactoryException e) {

            // TODO Bloque catch generado autom?ticamente

            e.printStackTrace();

        }

    }

    public String toString() {

        return projCS.toString();

    }

}

/*

 * Cresques Mapping Suite. Graphic Library for constructing mapping applications.

 *

 * Copyright (C) 2004-5.

 *

 * 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:

 *

 * cresques@gmail.com

 */

package org.cresques.impl.cts.gt2;

import java.awt.geom.Point2D;

import java.awt.geom.Rectangle2D;

import org.cresques.cts.ICoordTrans;

import org.cresques.cts.IProjection;

import org.geotools.ct.CannotCreateTransformException;

import org.geotools.ct.CoordinateTransformation;

import org.geotools.ct.CoordinateTransformationFactory;

import org.geotools.ct.MathTransform;

import org.geotools.pt.CoordinatePoint;

import org.opengis.referencing.operation.TransformException;

//import org.geotools.pt.MismatchedDimensionException;

/**

 * Transforma coordenadas entre dos sistemas

 * @see org.creques.cts.CoordSys

 * @author "Luis W. Sevilla" <sevilla_lui@gva.es>

 */

class CoordTrans implements ICoordTrans {

    private CoordinateTransformationFactory trFactory = CoordinateTransformationFactory.getDefault();

    private CoordinateTransformation tr = null;

    private MathTransform mt = null;

    private MathTransform mt2 = null;

    private MathTransform mt3 = null;

    private MathTransform mtDatum = null;

    private CoordSys from = null;

    private CoordSys to = null;

    private ICoordTrans invertedCT = null;

    public CoordTrans(CoordSys from, CoordSys to) {

        this.from = from;

        this.to = to;

        // Si los dos CoordSys son proyectados, entonces hay

        // que hacer dos transformaciones (pasar por geogr?ficas)

        // Si hay cambio de datum son 3 las transformaciones

        try {

            if (from.getDatum() != to.getDatum()) {

                tr = trFactory.createFromCoordinateSystems(from.toGeo().getCS(),

                                                           to.toGeo().getCS());

                mtDatum = tr.getMathTransform();

            }

            if ((from.projCS != null) && (to.projCS != null)) {

                CoordSys geogcs = from.toGeo();

                tr = trFactory.createFromCoordinateSystems(from.getCS(),

                                                           geogcs.getCS());

                mt = tr.getMathTransform();

                if (mtDatum != null) {

                    mt2 = mtDatum;

                }

                geogcs = to.toGeo();

                tr = trFactory.createFromCoordinateSystems(geogcs.getCS(),

                                                           to.getCS());

                if (mt2 == null) {

                    mt2 = tr.getMathTransform();

                } else {

                    mt3 = tr.getMathTransform();

                }

            } else {

                if (from.projCS == null) {

                    mt = mtDatum;

                }

                tr = trFactory.createFromCoordinateSystems(from.getCS(),

                                                           to.getCS());

                if (mt == null) {

                    mt = tr.getMathTransform();

                    if (mtDatum != null) {

                        mt2 = mtDatum;

                    }

                } else {

                    mt2 = tr.getMathTransform();

                }

            }

        } catch (CannotCreateTransformException e) {

            // TODO Bloque catch generado autom?ticamente

            e.printStackTrace();

        }

    }

    public IProjection getPOrig() {

        return from;

    }

    public IProjection getPDest() {

        return to;

    }

    public ICoordTrans getInverted() {

        if (invertedCT == null)

            invertedCT = new CoordTrans(to, from);

        return invertedCT;

    }

    public Point2D convert(Point2D ptOrig, Point2D ptDest) {

        CoordinatePoint pt1 = new CoordinatePoint(ptOrig);

        CoordinatePoint pt2 = new CoordinatePoint(0D, 0D);

        ptDest = null;

        try {

            mt.transform(pt1, pt2);

            ptDest = pt2.toPoint2D();

            if (mt2 != null) {

                mt2.transform(pt2, pt1);

                ptDest = pt1.toPoint2D();

                if (mt3 != null) {

                    mt3.transform(pt1, pt2);

                    ptDest = pt2.toPoint2D();

                }

            }

            /*} catch (MismatchedDimensionException e) {

                    // TODO Bloque catch generado autom?ticamente

                    e.printStackTrace();

            */

        } catch (TransformException e) {

            // TODO Bloque catch generado autom?ticamente

            e.printStackTrace();

        }

        return ptDest;

    }

    public String toString() {

        return tr.toString();

    }

    /* (non-Javadoc)

     * @see org.cresques.cts.ICoordTrans#convert(java.awt.geom.Rectangle2D)

     */

    public Rectangle2D convert(Rectangle2D rect) {

        Point2D pt1 = new Point2D.Double(rect.getMinX(), rect.getMinY());

        Point2D pt2 = new Point2D.Double(rect.getMaxX(), rect.getMaxY());

        pt1 = convert(pt1, null);

        pt2 = convert(pt2, null);

        rect = new Rectangle2D.Double();

        rect.setFrameFromDiagonal(pt1, pt2);

        return rect;

    }

}

/*

 * Cresques Mapping Suite. Graphic Library for constructing mapping applications.

 *

 * Copyright (C) 2004-6.

 *

 * 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:

 *

 * cresques@gmail.com

 */

package org.cresques.impl.cts.gt2;

import java.util.Locale;

import org.cresques.cts.IDatum;

import org.geotools.cs.CoordinateSystemFactory;

import org.geotools.cs.GeographicCoordinateSystem;

import org.geotools.cs.HorizontalDatum;

import org.opengis.referencing.FactoryException;

/**

 * Datum (y Ellipsoid) de GeoTools2.

 *

 * @author "Luis W. Sevilla" <sevilla_lui@gva.es>

 */

public class CSDatum implements IDatum {

    private static String line1 = "DATUM[\"WGS_1984\"," +

                                  "SPHEROID[\"WGS 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]]," +

                                  "TOWGS84[0,0,0,0,0,0,0],AUTHORITY[\"EPSG\",\"6326\"]]";

    private static String line2 = "DATUM[\"European_Datum_1950\"," +

                                  "SPHEROID[\"International 1924\",6378388,297,AUTHORITY[\"EPSG\",\"7022\"]]," +

                                  "TOWGS84[-84,-107,-120,0,0,0,0],AUTHORITY[\"EPSG\",\"6230\"]]";

    private static String line3 = "DATUM[\"Nouvelle_Triangulation_Francaise\"," +

                                  "SPHEROID[\"Clarke 1880 (IGN)\",6378249.2,293.466021293627, AUTHORITY[\"EPSG\",\"7011\"]]," +

                                  "TOWGS84[-168,-60,320,0,0,0,0],AUTHORITY[\"EPSG\",\"6275\"]]";

    private static String line4 = "DATUM[\"Datum 73\"," +

                                  "SPHEROID[\"International 1924\",6378388,297,AUTHORITY[\"EPSG\",\"7022\"]]," +

                                  "TOWGS84[-87,-98,-121,0,0,0,0],AUTHORITY[\"EPSG\",\"4274\"]]";

    private static String line5 = "DATUM[\"North_American_Datum_1927\"," +

                                  "SPHEROID[\"Clarke 1866\",6378206.4,294.978698213901,AUTHORITY[\"EPSG\",\"7008\"]]," +

                                  "TOWGS84[-3,142,183,0,0,0,0],AUTHORITY[\"EPSG\",\"6267\"]]";

    private static String line6 = "DATUM[\"North_American_Datum_1983\"," +

                                  "SPHEROID[\"GRS 1980\",6378137,298.257222101,AUTHORITY[\"EPSG\",\"7019\"]]," +

                                  "TOWGS84[0,0,0,0,0,0,0],AUTHORITY[\"EPSG\",\"6269\"]]";

    /*

     * INSERT INTO epsg_coordinatereferencesystem VALUES (

     * 4288, 'Loma Quintana', 1313, 'geographic 2D', 6422, 6288,

     * Null, Null, Null, Null, 'Geodetic survey.',

     * 'Superseded by La Canoa (code 4247).', '',

     * 'EPSG', '2004/01/06', '2003.37', 1, 0 );

     *

     *  DX (m) = -270.933

            DY (m) =  115.599

            DZ (m) = -360.226

            EX (") = -5.266

            EY (") = -1.238

              EZ (")  =  2.381

            FE (ppm) = -5.109

     */

    private static String line7 = "DATUM[\"Loma Quintana\"," +

                                  "SPHEROID[\"International 1924\",6378388,297,AUTHORITY[\"EPSG\",\"7022\"]]," +

                                  "TOWGS84[0,0,0,0,0,0,0],AUTHORITY[\"EPSG\",\"6288\"]]";

    /*

    # La Canoa

    <4247> +proj=longlat +ellps=intl +towgs84=-273.5,110.6,-357.9,0,0,0,0

    no_defs <>

    # PSAD56

    <4248> +proj=longlat +ellps=intl +towgs84=-288,175,-376,0,0,0,0 no_defs <>

     */

    private static String line8 = "DATUM[\"La Canoa\"," +

                                  "SPHEROID[\"International 1924\",6378388,297,AUTHORITY[\"EPSG\",\"7022\"]]," +