Application: gvSIG desktop » gvSIG projections

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

<!--

    gvSIG. Desktop Geographic Information System.

    Copyright (C) 2007-2013 gvSIG Association.

    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 3

    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.

    For any additional information, do not hesitate to contact us

    at info AT gvsig.com, or visit our website www.gvsig.com.

-->

<!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. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2013 gvSIG Association.

 *

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

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

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

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;

	}

}

/**

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2013 gvSIG Association.

 *

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

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

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

package org.gvsig.fmap.crs;

import org.cresques.cts.ICRSFactory;

import org.cresques.cts.IProjection;

import org.slf4j.Logger;

import org.slf4j.LoggerFactory;

/**

 * Fabrica de CRS.

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

 * @author luisw

 *

 */

public class CRSFactory {

	private static Logger logger = LoggerFactory.getLogger(CRSFactory.class);

	/**

	 * @deprecated use {@link #registerCRSFactory(ICRSFactory)}

	 */

	public static ICRSFactory cp = null;

	private static ICRSFactory currentcp = null;

	public static void registerCRSFactory(ICRSFactory cp) {

		logger.info("Register CRSFactory " + cp.getClass().getName() + "/" + cp.toString() + ".");

		CRSFactory.currentcp = cp;

		CRSFactory.cp = cp;

	}

	public static ICRSFactory getCRSFactory() {

		return currentcp;

	}

	public static IProjection getCRS(String code) {

		if( cp != currentcp ) {

			logger.warn("Changed CRSFactory without calling registerCRSFactory.");

			registerCRSFactory(cp);

		}

		return currentcp.get(code);

	}

}

/**

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2013 gvSIG Association.

 *

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

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

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

/*

 * 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. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2013 gvSIG Association.

 *

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

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

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

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.addDynFieldObject(FIELD_SOURCE_PROJECTION).setDescription(

            DESCRIPTION_SOURCE_PROJECTION).setMandatory(true).setClassOfValue(IProjection.class);

        definition.addDynFieldObject(FIELD_DESTINATION_PROJECTION).setDescription(

            DESCRIPTION_DESTINATION_PROJECTION).setMandatory(true).setClassOfValue(IProjection.class);

    }

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

    }

}

/**

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2013 gvSIG Association.

 *

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

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

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

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. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2013 gvSIG Association.

 *

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

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

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

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. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2013 gvSIG Association.

 *

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

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

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

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

import org.slf4j.Logger;

import org.slf4j.LoggerFactory;

public class ProjectionLibrary extends AbstractLibrary {

	private static final Logger logger = LoggerFactory.getLogger(ProjectionLibrary.class);

    public void doRegistration() {

        registerAsAPI(ProjectionLibrary.class);

        require(ToolsLibrary.class);

    }

	protected void doInitialize() throws LibraryException {

	}

	protected void doPostInitialize() throws LibraryException {

		if (CRSFactory.getCRSFactory()==null) {

			logger.warn("has not registered an implementation of " + this.getClass().getName()+".");

		}

		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. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2013 gvSIG Association.

 *

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

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

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

package org.cresques;

import java.util.Locale;

import org.gvsig.tools.ToolsLocator;

import org.gvsig.tools.i18n.I18nManager;

/**

* Bridge class to provide internationalization services to the library.

* It uses the gvsig-i18n library as a backend, and includes its

* necessary initialization.

* 

*/

public class Messages {

	private static I18nManager manager = null;

	/**

	 * Loads the translations in the dictionary. It initializes the backend

	 * gvsig-i18n library

	 *

	 */

	private static void init() {

		manager = ToolsLocator.getI18nManager(); 

		manager.addResourceFamily("org.cresques.resources.i18n.text", Messages.class.getClassLoader(), Messages.class.getClass().getName());

	}

	/**

	 * Gets the translation associated with the provided translation key.

	 * 

	 * @param key The translation key which identifies the target text

	 * @return The translation associated with the provided translation key.

	 */

	public static String getText(String key) {

		if (manager == null ) {

			init();

		}

		return manager.getTranslation(key);

	}

}

/**

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2013 gvSIG Association.

 *

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

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

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

package org.cresques.px;

import java.awt.geom.Point2D;

import java.awt.geom.Rectangle2D;

import java.text.DecimalFormat;

/**

 *	Clase que getiona el extent de una imagen

 *	

 *  @author Luis W.Sevilla (sevilla_lui@gva.es)

 */

public class Extent {

    Point2D min = null;

    Point2D max = null;

    /**

     * Constructor sin par?metros

     */

    public Extent() {

        min = new Point2D.Double(999999999.0, 999999999.0);

        max = new Point2D.Double(-999999999.0, -999999999.0);

    }

    /**

     * Constructor 

     * @param pt1	punto que representa la esquina superior izquierda

     * @param pt2	punto que representa la esquina inferior derecha

     */

    public Extent(Point2D pt1, Point2D pt2) {

        newExtent(pt1.getX(), pt1.getY(), pt2.getX(), pt2.getY());

    }

    /**

     * Contructor

     * @param x1 punto que representa la coordenada X de la esquina superior izquierda

     * @param y1 punto que representa la coordenada Y de la esquina superior izquierda

     * @param x2 punto que representa la coordenada X de la esquina inferior derecha

     * @param y2 punto que representa la coordenada Y de la esquina inferior derecha

     */

    public Extent(double x1, double y1, double x2, double y2) {

        newExtent(x1, y1, x2, y2);

    }

    /**

     * Constructor

     * @param r	Rectangulo 2D

     */

    public Extent(Rectangle2D r) {

        newExtent(r.getX(), r.getY(), r.getX() + r.getWidth(),

                  r.getY() + r.getHeight());

    }

    /**

     * Constructor de copia

     * @param ext	Objeto Extent

     */

    public Extent(Extent ext) {

        newExtent(ext.minX(), ext.minY(), ext.maxX(), ext.maxY());

    }

    /**

     * Crea un objeto extent identico y lo retorna

     * @return Objeto extent

     */

    public Object clone() {

        Extent e = (Extent) clone();

        e.min = (Point2D) min.clone();

        e.max = (Point2D) max.clone();

        return e;

    }

    private void newExtent(double x1, double y1, double x2, double y2) {

        min = new Point2D.Double(Math.min(x1, x2), Math.min(y1, y2));

        max = new Point2D.Double(Math.max(x1, x2), Math.max(y1, y2));

    }

    /**

     * Obtiene la coordenada X m?nima

     * @return valor de la coordenada X m?nima

     */

    public double minX() {

        return min.getX();

    }

    /**

     * Obtiene la coordenada Y m?nima

     * @return valor de la coordenada X m?nima

     */

    public double minY() {

        return min.getY();

    }

    /**

     * Obtiene la coordenada X m?xima

     * @return valor de la coordenada X m?xima

     */

    public double maxX() {

        return max.getX();

    }

    /**

     * Obtiene la coordenada Y m?xima

     * @return valor de la coordenada Y m?xima

     */

    public double maxY() {

        return max.getY();

    }

    /**

     * Obtiene el punto m?nimo

     * @return m?nimo

     */

    public Point2D getMin() {

        return min;

    }

    /**

     * Obtiene el punto m?ximo

     * @return m?ximo

     */

    public Point2D getMax() {

        return max;

    }

    public boolean isAt(Point2D pt) {

        if (pt.getX() < minX()) {

            return false;

        }

        if (pt.getX() > maxX()) {

            return false;

        }

        if (pt.getY() < minY()) {

            return false;

        }

        if (pt.getY() > maxY()) {

            return false;

        }

        return true;

    }

    public double width() {

        return Math.abs(maxX() - minX());

    }

    public double height() {

        return Math.abs(maxY() - minY());

    }

    /**

     * Verifica un punto, y modifica el extent si no est? incluido

     */

    public void add(Point2D pt) {

        if (pt == null) {

            return;

        }

        min.setLocation(Math.min(pt.getX(), minX()), Math.min(pt.getY(), minY()));

        max.setLocation(Math.max(pt.getX(), maxX()), Math.max(pt.getY(), maxY()));

    }

    public void add(Extent ext) {

        if (ext == null) {

            return;

        }

        min.setLocation(Math.min(ext.minX(), minX()),

                        Math.min(ext.minY(), minY()));

        max.setLocation(Math.max(ext.maxX(), maxX()),

                        Math.max(ext.maxY(), maxY()));

    }

    /**

     * Obtiene la escala

     * @param width	Ancho

     * @param height	Alto

     * @return

     */

    public double[] getScale(int width, int height) {

        return getScale((double) width, (double) height);

    }

    public double[] getScale(double width, double height) {

        double[] scale = new double[2];

        scale[0] = ((float) width) / width();

        scale[1] = ((float) height) / height();

        return scale;

    }

    public Rectangle2D toRectangle2D() {

        return new Rectangle2D.Double(minX(), minY(), width(), height());

    }

    public String toString() {

        DecimalFormat format = new DecimalFormat("####.000");

        return "Extent: (" + format.format(minX()) + "," +

               format.format(minY()) + "), (" + format.format(maxX()) + "," +

               format.format(maxY()) + ")";

    }

    public interface Has {

        public Extent getExtent();

    }

}

*.dfPackage

*.wmf

<!--

    gvSIG. Desktop Geographic Information System.

    Copyright (C) 2007-2013 gvSIG Association.

    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 3

    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.

    For any additional information, do not hesitate to contact us

    at info AT gvsig.com, or visit our website www.gvsig.com.

-->

<html>

	<body>Pixel: Clases para dibujar.

</body>

</html>

/**

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2013 gvSIG Association.

 *

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

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

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

package org.cresques.cts;

/**

 * @author "Luis W. Sevilla" (sevilla_lui@gva.es)

 */

public interface IDatum {

    /**

     * Semieje Mayor del elipsoide.

     * @return

     */

    public double getESemiMajorAxis();

    /**

     * Aplanamiento del elipsoide.

     * @return

     */

    public double getEIFlattening();

}

/**

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2013 gvSIG Association.

 *

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

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

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

package org.cresques.cts;

import org.cresques.geo.ViewPortData;

import java.awt.Color;

import java.awt.Graphics2D;

import java.awt.geom.Point2D;

import java.awt.geom.Rectangle2D;

import org.gvsig.tools.lang.Cloneable;

/**

 *

 * @author "Luis W. Sevilla" (sevilla_lui@gva.es)

 */

public interface IProjection extends Cloneable{

    public IDatum getDatum();

    public Point2D createPoint(double x, double y);

    // TODO Quitar si no son necesarias.

    public String getAbrev();

    /**

     * Devuelve getAbrev() mas los parametros de transformacion si los hay

     * ej.: (EPSG:23030:proj@+proj...@...)

     * 

     * @return getAbrev() o getAbrev()+parametros

     */

    public String getFullCode();

    public void drawGrid(Graphics2D g, ViewPortData vp);

    public void setGridColor(Color c);

    public Color getGridColor();

    /**

     * Crea un ICoordTrans para transformar coordenadas

     * desde el IProjection actual al dest.

     * @param dest

     * @return

     */

    public ICoordTrans getCT(IProjection dest);

    public Point2D toGeo(Point2D pt);

    public Point2D fromGeo(Point2D gPt, Point2D mPt);

    public boolean isProjected();

    /**

     * First two parameters must be in meters.

     * This should be changed (map units should be used) and then

     * change the places where this method is used.

     * 

     * @param minX in meters

     * @param maxX in meters

     * @param width in pixels (dots)

     * @param dpi dots per inch

     * @return Scale denominator ( the "X" in "1 : X" )

     */

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

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

}

/**

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2013 gvSIG Association.

 *

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

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

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

package org.cresques.cts;

public interface UTM extends IProjection {

}

/**

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2013 gvSIG Association.

 *

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

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

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

package org.cresques.cts;

import java.awt.geom.Point2D;

import java.awt.geom.Rectangle2D;

/**

 * @author "Luis W. Sevilla" (sevilla_lui@gva.es)

 */

public interface ICoordTrans {

    public IProjection getPOrig();

    public IProjection getPDest();

    public Point2D convert(Point2D ptOrig, Point2D ptDest);

    /**

     * Reprojects the input rectangle and returns the minimum rectangle

     * containing the result of that reprojection (which was not

     * necessarily a rectangle).

     * 

     * Reprojecting the diagonal an assuming that

     * the result is the diagonal of the rectangle to be returned is wrong

     * 

     * Reprojecting the four corners and searching for min and max

     * to create the resulting rectangle is not completely accurate

     * but provides an acceptable approximation 

     * 

     * @param rect

     * @return

     */

    public Rectangle2D convert(Rectangle2D rect);

    public ICoordTrans getInverted();

}

/**

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2013 gvSIG Association.

 *

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

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

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

package org.cresques.cts;

public interface ICRSFactory {

	public boolean doesRigurousTransformations();

    /**

     * Devuelve una proyeccion a partir de una cadena.

     * @param name abreviatura de la proyecccion (i.e. EPSG:23030)

     * @return Proyeccion si existe

     */

    public IProjection get(String name);

}

/**

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2013 gvSIG Association.

 *

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

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

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

package org.cresques.cts;

import java.awt.geom.Point2D;

/**

 * Operaciones relacionadas con las proyecciones y sistemas

 * de coordenadas.

 *

 * cmartinez: Esta clase no deber?a formar parte de una API, pero

 * se deja hasta que se aborde el refactoring de libProjection.

 *

 * @author Luis W. Sevilla (sevilla_lui@gva.es)

 */

public class GeoCalc {

    IProjection proj;

    /**

     *

     * @param proj

     */

    public GeoCalc(IProjection proj) {

        this.proj = proj;

    }

    /* (non-Javadoc)

	 * @see org.cresques.impl.cts.GeoCalc#distanceGeo(java.awt.geom.Point2D, java.awt.geom.Point2D)

	 */

    public double distanceGeo(Point2D pt1, Point2D pt2) {

        double R2 = Math.pow(proj.getDatum().getESemiMajorAxis(), 2);

        double dLat = Math.toRadians(pt2.getY() - pt1.getY());

        double dLong = Math.toRadians(pt2.getX() - pt1.getX());

        double alfa = Math.toRadians(pt1.getY());

        double alfa2 = Math.toRadians(pt2.getY());

        if (Math.abs(alfa2) < Math.abs(alfa)) {

            alfa = alfa2;

        }

        double ds2 = (R2 * dLat * dLat) +

                     (R2 * Math.cos(alfa) * Math.cos(alfa) * dLong * dLong);

        return Math.sqrt(ds2);

    }

    /* (non-Javadoc)

	 * @see org.cresques.impl.cts.GeoCalc#distanceEli(java.awt.geom.Point2D, java.awt.geom.Point2D)

	 */

    public double distanceEli(Point2D pt1, Point2D pt2) {

        double lat1 = Math.toRadians(pt1.getY());

        double lon1 = -Math.toRadians(pt1.getX());

        double lat2 = Math.toRadians(pt2.getY());

        double lon2 = -Math.toRadians(pt2.getX());

        double F = (lat1 + lat2) / 2D;

        double G = (lat1 - lat2) / 2D;

        double L = (lon1 - lon2) / 2D;

        double sing = Math.sin(G);

        double cosl = Math.cos(L);

        double cosf = Math.cos(F);

        double sinl = Math.sin(L);

        double sinf = Math.sin(F);

        double cosg = Math.cos(G);

        double flat = 1D / proj.getDatum().getEIFlattening();

        double S = (sing * sing * cosl * cosl) + (cosf * cosf * sinl * sinl);

        double C = (cosg * cosg * cosl * cosl) + (sinf * sinf * sinl * sinl);

        double W = Math.atan2(Math.sqrt(S), Math.sqrt(C));

        double R = Math.sqrt((S * C)) / W;

        double H1 = ((3D * R) - 1D) / (2D * C);

        double H2 = ((3D * R) + 1D) / (2D * S);

        double D = 2D * W * proj.getDatum().getESemiMajorAxis();

        return (D * ((1D + (flat * H1 * sinf * sinf * cosg * cosg)) -

               (flat * H2 * cosf * cosf * sing * sing)));

    }

    /**

     * Algrothims from Geocentric Datum of Australia Technical Manual

     *

     * http://www.anzlic.org.au/icsm/gdatum/chapter4.html

     *

     * This page last updated 11 May 1999

     *

     * Computations on the Ellipsoid

     *

     * There are a number of formulae that are available

     * to calculate accurate geodetic positions,

     * azimuths and distances on the ellipsoid.

     *

     * Vincenty's formulae (Vincenty, 1975) may be used

     * for lines ranging from a few cm to nearly 20,000 km,

     * with millimetre accuracy.

     * The formulae have been extensively tested

     * for the Australian region, by comparison with results

     * from other formulae (Rainsford, 1955 & Sodano, 1965).

     *

     * * Inverse problem: azimuth and distance from known

     *                 latitudes and longitudes

     * * Direct problem: Latitude and longitude from known

     *                 position, azimuth and distance.

     * * Sample data

     * * Excel spreadsheet

     *

     * Vincenty's Inverse formulae

     * Given: latitude and longitude of two points

     *                 (phi1, lembda1 and phi2, lembda2),

     * Calculate: the ellipsoidal distance (s) and

     * forward and reverse azimuths between the points (alpha12, alpha21).

     */

    /* (non-Javadoc)

	 * @see org.cresques.impl.cts.GeoCalc#distanceVincenty(java.awt.geom.Point2D, java.awt.geom.Point2D)

	 */

    public double distanceVincenty(Point2D pt1, Point2D pt2) {

        return distanceAzimutVincenty(pt1, pt2).dist;

    }

    /**

     * Returns the distance between two geographic points on the ellipsoid

     *        and the forward and reverse azimuths between these points.

     *       lats, longs and azimuths are in decimal degrees, distance in metres

     *        Returns ( s, alpha12,  alpha21 ) as a tuple

     * @param pt1

     * @param pt2

     * @return

     */

    protected GeoData distanceAzimutVincenty(Point2D pt1, Point2D pt2) {

        GeoData gd = new GeoData(0, 0);

        double f = 1D / proj.getDatum().getEIFlattening();

        double a = proj.getDatum().getESemiMajorAxis();

        double phi1 = pt1.getY();

        double lembda1 = pt1.getX();

        double phi2 = pt2.getY();

        double lembda2 = pt2.getX();

        if ((Math.abs(phi2 - phi1) < 1e-8) &&

                (Math.abs(lembda2 - lembda1) < 1e-8)) {

            return gd;

        }

        double piD4 = Math.atan(1.0);

        double two_pi = piD4 * 8.0;

        phi1 = (phi1 * piD4) / 45.0;

        lembda1 = (lembda1 * piD4) / 45.0; // unfortunately lambda is a key word!

        phi2 = (phi2 * piD4) / 45.0;

        lembda2 = (lembda2 * piD4) / 45.0;

        double b = a * (1.0 - f);

        double TanU1 = (1 - f) * Math.tan(phi1);

        double TanU2 = (1 - f) * Math.tan(phi2);

        double U1 = Math.atan(TanU1);

        double U2 = Math.atan(TanU2);

        double lembda = lembda2 - lembda1;

        double last_lembda = -4000000.0; // an impossibe value

        double omega = lembda;

        // Iterate the following equations,