Application: gvSIG desktop » gvSIG projections

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

<project xmlns="http://maven.apache.org/POM/4.0.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://maven.apache.org/POM/4.0.0 http://maven.apache.org/xsd/maven-4.0.0.xsd">

    <modelVersion>4.0.0</modelVersion>

    <groupId>org.gvsig</groupId>

    <artifactId>org.gvsig.projection.api</artifactId>

    <version>2.0.47</version>

    <packaging>jar</packaging>

    <name>${project.artifactId}</name>

    <licenses>

        <license>

            <name>GPLv3</name>

            <url>http://www.gnu.org/licenses/gpl.html</url>

            <distribution>repo</distribution>

            <comments>

                "GPL" stands for "General Public License". The most widespread such

                license is the GNU General Public License, or GNU GPL for short.

                This can be further shortened to "GPL", when it is understood that

                the GNU GPL is the one intended.

            </comments>

        </license>

    </licenses>

    <organization>

        <name>gvSIG Association</name>

        <url>http://www.gvsig.com/</url>

    </organization>

    <scm>

        <connection>scm:svn:https://devel.gvsig.org/svn/gvsig-jcrs/org.gvsig.projection/tags/org.gvsig.projection.api-2.0.47</connection>

        <developerConnection>scm:svn:https://devel.gvsig.org/svn/gvsig-jcrs/org.gvsig.projection/tags/org.gvsig.projection.api-2.0.47</developerConnection>

        <url>https://devel.gvsig.org/redmine/projects/gvsig-jcrs/repository/show/org.gvsig.projection/tags/org.gvsig.projection.api-2.0.47</url>

    </scm>

    <repositories>

        <repository>

            <id>gvsig-public-http-repository</id>

            <name>gvSIG maven public HTTP repository</name>

            <url>http://devel.gvsig.org/m2repo/j2se</url>

            <releases>

                <enabled>true</enabled>

                <updatePolicy>daily</updatePolicy>

                <checksumPolicy>warn</checksumPolicy>

            </releases>

            <snapshots>

                <enabled>true</enabled>

                <updatePolicy>daily</updatePolicy>

                <checksumPolicy>warn</checksumPolicy>

            </snapshots>

        </repository>

    </repositories>

    <distributionManagement>

        <site>

            <id>gvsig-repository</id>

            <url>dav:https://devel.gvsig.org/sites/${project.artifactId}/${project.version}</url>

        </site>

        <repository>

            <id>gvsig-repository</id>

            <name>gvSIG maven repository</name>

            <url>dav:https://devel.gvsig.org/m2repo/j2se</url>

            <uniqueVersion>false</uniqueVersion>

        </repository>

        <snapshotRepository>

            <id>gvsig-repository</id>

            <name>gvSIG maven repository</name>

            <url>dav:https://devel.gvsig.org/m2repo/j2se</url>

            <uniqueVersion>false</uniqueVersion>

        </snapshotRepository>

    </distributionManagement>

    <issueManagement>

        <system>redmine</system>

        <url>https://devel.gvsig.org/redmine/projects/gvsig-jcrs</url>

    </issueManagement>

    <mailingLists>

        <mailingList>

            <name>User List</name>

            <subscribe>http://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_usuarios</subscribe>

            <unsubscribe>http://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_usuarios</unsubscribe>

            <post>gvsig_usuarios@listserv.gva.es</post>

            <archive>http://www.gvsig.org/web/community/mailing-lists/gvsig-nabble</archive>

            <otherArchives>

                <otherArchive>http://www.mail-archive.com/gvsig_usuarios@runas.cap.gva.es</otherArchive>

                <otherArchive>http://dir.gmane.org/index.php?prefix=gmane.comp.gis.gvsig</otherArchive>

            </otherArchives>

        </mailingList>

        <mailingList>

            <name>Developer List</name>

            <subscribe>http://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_desarrolladores</subscribe>

            <unsubscribe>http://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_desarrolladores</unsubscribe>

            <post>gvsig_desarrolladores@listserv.gva.es</post>

            <archive>http://www.gvsig.org/web/community/mailing-lists/gvsig-nabble</archive>

            <otherArchives>

                <otherArchive>http://www.mail-archive.com/gvsig_desarrolladores@runas.cap.gva.es</otherArchive>

                <otherArchive>http://dir.gmane.org/index.php?prefix=gmane.comp.gis.gvsig</otherArchive>

            </otherArchives>

        </mailingList>

        <mailingList>

            <name>International List</name>

            <subscribe>http://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_internacional</subscribe>

            <unsubscribe>http://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_internacional</unsubscribe>

            <post>gvsig_internacional@listserv.gva.es</post>

            <archive>http://www.gvsig.org/web/community/mailing-lists/gvsig-nabble</archive>

            <otherArchives>

                <otherArchive>http://www.mail-archive.com/gvsig_internacional@runas.cap.gva.es</otherArchive>

                <otherArchive>http://dir.gmane.org/index.php?prefix=gmane.comp.gis.gvsig</otherArchive>

            </otherArchives>

        </mailingList>

    </mailingLists>

    <developers>

        <developer>

            <id>jjdelcerro</id>

            <name>Joaquín José del Cerro</name>

            <email>jjdelcerro@gvsig.com</email>

            <roles>

                <role>Architect</role>

                <role>Developer</role>

            </roles>

        </developer>

        <developer>

            <id>cordinayana</id>

            <name>Cèsar Ordiñana</name>

            <email>cordinyana@gvsig.com</email>

            <roles>

                <role>Architect</role>

                <role>Developer</role>

            </roles>

        </developer>

    </developers>

    <dependencyManagement>

        <dependencies>

            <dependency>

                <groupId>junit</groupId>

                <artifactId>junit</artifactId>

                <version>3.8.1</version>

            </dependency>

            <dependency>

                <groupId>org.gvsig</groupId>

                <artifactId>org.gvsig.tools.lib</artifactId>

                <version>${gvsig.tools.version}</version>

            </dependency>

        </dependencies>

    </dependencyManagement>

    <dependencies>

        <dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.tools.lib</artifactId>

        </dependency>

        <dependency>

            <groupId>junit</groupId>

            <artifactId>junit</artifactId>

            <scope>test</scope>

        </dependency>            

    </dependencies>

    <build>

        <extensions>

            <extension>

                <groupId>org.apache.maven.wagon</groupId>

                <artifactId>wagon-webdav-jackrabbit</artifactId>

                <version>2.4</version>

            </extension>

        </extensions>

        <pluginManagement>

            <plugins>

                <plugin>

                    <groupId>org.apache.maven.plugins</groupId>

                    <artifactId>maven-surefire-plugin</artifactId>

                    <version>2.7.1</version>

                    <configuration>

                        <forkMode>once</forkMode>

                    </configuration>

                </plugin>

                <plugin>

                    <groupId>org.apache.maven.plugins</groupId>

                    <artifactId>maven-compiler-plugin</artifactId>

                    <version>2.3.2</version>

                    <configuration>

                        <source>1.5</source>

                        <target>1.5</target>

                        <encoding>ISO-8859-1</encoding>

                    </configuration>

                </plugin>

                <plugin>

                    <groupId>org.apache.maven.plugins</groupId>

                    <artifactId>maven-resources-plugin</artifactId>

                    <version>2.4.3</version>

                    <configuration>

                        <encoding>ISO-8859-1</encoding>

                    </configuration>

                </plugin>

                <plugin>

                    <groupId>org.apache.maven.plugins</groupId>

                    <artifactId>maven-javadoc-plugin</artifactId>

                    <version>2.9</version>

                    <configuration>

                        <docTitle>${project.name} ${project.version}</docTitle>

                        <encoding>ISO-8859-1</encoding>

                        <footer />

                        <header />

                        <keywords>true</keywords>

                        <detectLinks>true</detectLinks>

                        <javaApiLinks>

                            <property>

                                <name>api_1.4</name>

                                <value>http://download.oracle.com/javase/1.4.2/docs/api/</value>

                            </property>

                            <property>

                                <name>api_1.5</name>

                                <value>http://download.oracle.com/javase/1.5.0/docs/api/</value>

                            </property>

                            <property>

                                <name>api_1.6</name>

                                <value>http://download.oracle.com/javase/6/docs/api/</value>

                            </property>

                        </javaApiLinks>

                    </configuration>

                </plugin>

                <plugin>

                    <groupId>org.apache.maven.plugins</groupId>

                    <artifactId>maven-site-plugin</artifactId>

                    <configuration>

                        <chmod>true</chmod>

                        <chmodMode>a+rwX</chmodMode>

                    </configuration>

                </plugin>

                <plugin>

                    <groupId>org.apache.maven.plugins</groupId>

                    <artifactId>maven-jar-plugin</artifactId>

                    <version>2.3.1</version>

                    <configuration>

                        <archive>

                            <manifestEntries>

                                <Built-By>gvSIG Team</Built-By>

                            </manifestEntries>

                        </archive>

                    </configuration>

                </plugin>

                <plugin>

                    <groupId>org.codehaus.mojo</groupId>

                    <artifactId>properties-maven-plugin</artifactId>

                    <version>1.0-alpha-2</version>

                </plugin>

                <plugin>

                    <groupId>org.codehaus.gmaven</groupId>

                    <artifactId>gmaven-plugin</artifactId>

                    <version>1.4</version>

                </plugin>

                <plugin>

                    <groupId>org.apache.maven.plugins</groupId>

                    <artifactId>maven-antrun-plugin</artifactId>

                    <version>1.7</version>

                </plugin>

                <plugin>

                    <groupId>org.codehaus.mojo</groupId>

                    <artifactId>animal-sniffer-maven-plugin</artifactId>

                    <version>1.9</version>

                </plugin>

                <plugin>

                    <groupId>org.apache.maven.plugins</groupId>

                    <artifactId>maven-release-plugin</artifactId>

                    <version>2.1</version>

                </plugin>

                <plugin>

                    <groupId>org.apache.maven.plugins</groupId>

                    <artifactId>maven-assembly-plugin</artifactId>

                    <version>2.2-beta-5</version>

                </plugin>

                <plugin>

                    <groupId>org.apache.maven.plugins</groupId>

                    <artifactId>maven-clean-plugin</artifactId>

                    <version>2.4.1</version>

                </plugin>

                <plugin>

                    <groupId>org.codehaus.mojo</groupId>

                    <artifactId>exec-maven-plugin</artifactId>

                    <version>1.2</version>

                </plugin>

                <plugin>

                    <!-- invoke as "mvn license:format" to fix the header -->

                    <!-- invoke as "mvn license:check" to check the header -->

                    <!-- See http://code.google.com/p/maven-license-plugin/wiki/Configuration -->

                    <groupId>com.mycila.maven-license-plugin</groupId>

                    <artifactId>maven-license-plugin</artifactId>

                    <version>1.10.b1</version>

                    <configuration>

                        <header>/src/main/config/header.txt</header>

                        <failIfMissing>false</failIfMissing>

                        <aggregate>false</aggregate>

                        <strictCheck>true</strictCheck>

                        <encoding>ISO-8859-1</encoding>

                        <quiet>false</quiet>

                        <excludes>

                            <exclude>**/*.properties</exclude>

                            <exclude>**/*.txt</exclude>

                            <exclude>**/*.jpg</exclude>

                            <exclude>**/*.jpeg</exclude>

                            <exclude>**/*.svg</exclude>

                            <exclude>**/*.png</exclude>

                            <exclude>**/*.gvssym</exclude>

                            <exclude>**/*.def</exclude>

                            <exclude>**/*.msg</exclude>

                        </excludes>

                    </configuration>

                </plugin>

                <plugin>

                    <groupId>org.apache.maven.plugins</groupId>

                    <artifactId>maven-source-plugin</artifactId>

                    <version>2.2.1</version>

                </plugin>

            </plugins>

        </pluginManagement>

        <plugins>

            <plugin>

                <groupId>org.apache.maven.plugins</groupId>

                <artifactId>maven-jar-plugin</artifactId>

                <configuration>

                </configuration>

                <executions>

                    <!-- Generates a jar file only with the test classes -->

                    <execution>

                        <goals>

                            <goal>test-jar</goal>

                        </goals>

                    </execution>

                </executions>

            </plugin>

            <plugin>

                <groupId>org.apache.maven.plugins</groupId>

                <artifactId>maven-release-plugin</artifactId>

                <configuration>

                    <tagBase>https://devel.gvsig.org/svn/gvsig-jcrs/${gvsig.projection.svnroot}/tags/</tagBase>

                    <!-- This is a 2.0 beta 6 bug. By default the goals are 'clean verify', so the modules

                    are not installed, so the dependencies between modules are not met. With this configuration

                    we force the module installation. -->

                    <preparationGoals>clean install</preparationGoals>

                    <autoVersionSubmodules>true</autoVersionSubmodules>

                </configuration>

            </plugin>

            <plugin>

                <!-- to continue building  when test fail -->

                <groupId>org.apache.maven.plugins</groupId>

                <artifactId>maven-surefire-plugin</artifactId>

                <configuration>

                    <testFailureIgnore>true</testFailureIgnore>

                </configuration>

            </plugin>

        </plugins>

    </build>

    <reporting>

        <excludeDefaults>true</excludeDefaults>

        <plugins>

            <plugin>

                <groupId>org.apache.maven.plugins</groupId>

                <artifactId>maven-project-info-reports-plugin</artifactId>

                <version>2.7</version>

            </plugin>

            <plugin>

                <groupId>org.apache.maven.plugins</groupId>

                <artifactId>maven-javadoc-plugin</artifactId>

                <version>2.9</version>

            </plugin>

            <plugin>

                <groupId>org.apache.maven.plugins</groupId>

                <artifactId>maven-jxr-plugin</artifactId>

                <version>2.3</version>

                <configuration>

                    <linkJavadoc>true</linkJavadoc>

                </configuration>

            </plugin>

        </plugins>

    </reporting>

    <properties>

        <project.build.sourceEncoding>ISO-8859-1</project.build.sourceEncoding>

        <localRepository>${user.home}/.m2/repository</localRepository>

        <site-repository>dav:https://devel.gvsig.org/sites/${project.artifactId}/${project.version}</site-repository>

        <gvsig.projection.svnroot>org.gvsig.projection</gvsig.projection.svnroot>

        <gvsig.tools.version>3.0.252</gvsig.tools.version>

    </properties>

</project>

<?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.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,

        //  until there is no significant change in lembda

        double Sin_sigma = 0;

        // Iterate the following equations,

        //  until there is no significant change in lembda

        double Cos_sigma = 0;

        double Cos2sigma_m = 0;

        double alpha = 0;

        double sigma = 0;

        double sqr_sin_sigma = 0;

        while ((last_lembda < -3000000.0) ||

                   ((lembda != 0) &&

                   (Math.abs((last_lembda - lembda) / lembda) > 1.0e-9))) {

            sqr_sin_sigma = Math.pow(Math.cos(U2) * Math.sin(lembda), 2) +

                            Math.pow(((Math.cos(U1) * Math.sin(U2)) -

                                     (Math.sin(U1) * Math.cos(U2) * Math.cos(lembda))),

                                     2);

            Sin_sigma = Math.sqrt(sqr_sin_sigma);

            Cos_sigma = (Math.sin(U1) * Math.sin(U2)) +

                        (Math.cos(U1) * Math.cos(U2) * Math.cos(lembda));

            sigma = Math.atan2(Sin_sigma, Cos_sigma);

            double Sin_alpha = (Math.cos(U1) * Math.cos(U2) * Math.sin(lembda)) / Math.sin(sigma);

            alpha = Math.asin(Sin_alpha);

            Cos2sigma_m = Math.cos(sigma) -

                          ((2 * Math.sin(U1) * Math.sin(U2)) / Math.pow(Math.cos(alpha),

                                                                        2));

            double C = (f / 16) * Math.pow(Math.cos(alpha), 2) * (4 +

                       (f * (4 - (3 * Math.pow(Math.cos(alpha), 2)))));

            last_lembda = lembda;

            lembda = omega +

                     ((1 - C) * f * Math.sin(alpha) * (sigma +

                     (C * Math.sin(sigma) * (Cos2sigma_m +

                     (C * Math.cos(sigma) * (-1 +

                     (2 * Math.pow(Cos2sigma_m, 2))))))));

        }

        double u2 = (Math.pow(Math.cos(alpha), 2) * ((a * a) - (b * b))) / (b * b);

        double A = 1 +

                   ((u2 / 16384) * (4096 +

                   (u2 * (-768 + (u2 * (320 - (175 * u2)))))));

        double B = (u2 / 1024) * (256 +

                   (u2 * (-128 + (u2 * (74 - (47 * u2))))));

        double delta_sigma = B * Sin_sigma * (Cos2sigma_m +

                             ((B / 4) * ((Cos_sigma * (-1 +

                             (2 * Math.pow(Cos2sigma_m, 2)))) -

                             ((B / 6) * Cos2sigma_m * (-3 +

                             (4 * sqr_sin_sigma)) * (-3 +

                             (4 * Math.pow(Cos2sigma_m, 2)))))));

        double s = b * A * (sigma - delta_sigma);

        double alpha12 = Math.atan2((Math.cos(U2) * Math.sin(lembda)),

                                    ((Math.cos(U1) * Math.sin(U2)) -

                                    (Math.sin(U1) * Math.cos(U2) * Math.cos(lembda))));

        double alpha21 = Math.atan2((Math.cos(U1) * Math.sin(lembda)),

                                    ((-Math.sin(U1) * Math.cos(U2)) +

                                    (Math.cos(U1) * Math.sin(U2) * Math.cos(lembda))));

        if (alpha12 < 0.0) {

            alpha12 = alpha12 + two_pi;

        }

        if (alpha12 > two_pi) {

            alpha12 = alpha12 - two_pi;

        }

        alpha21 = alpha21 + (two_pi / 2.0);

        if (alpha21 < 0.0) {

            alpha21 = alpha21 + two_pi;

        }

        if (alpha21 > two_pi) {

            alpha21 = alpha21 - two_pi;

        }

        alpha12 = (alpha12 * 45.0) / piD4;

        alpha21 = (alpha21 * 45.0) / piD4;

        return new GeoData(0, 0, s, alpha12, alpha21);

    }

    /**

     * Vincenty's Direct formulae

     * Given: latitude and longitude of a point (phi1, lembda1) and

     * the geodetic azimuth (alpha12)

     * and ellipsoidal distance in metres (s) to a second point,

     *

     * Calculate: the latitude and longitude of the second point (phi2, lembda2)

     * and the reverse azimuth (alpha21).

     */

	/**

	 * Returns the lat and long of projected point and reverse azimuth

	 *        given a reference point and a distance and azimuth to project.

	 *  lats, longs and azimuths are passed in decimal degrees.

	 * Returns ( phi2,  lambda2,  alpha21 ) as a tuple

	 * @param pt

	 * @param azimut

	 * @param dist

	 * @return

	 */

    protected GeoData getPointVincenty(Point2D pt, double azimut, double dist) {

        GeoData ret = new GeoData(0, 0);

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

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

        double phi1 = pt.getY();

        double lembda1 = pt.getX();

        double alpha12 = azimut;

        double s = dist;

        double piD4 = Math.atan(1.0);

        double two_pi = piD4 * 8.0;

        phi1 = (phi1 * piD4) / 45.0;

        lembda1 = (lembda1 * piD4) / 45.0;

        alpha12 = (alpha12 * piD4) / 45.0;

        if (alpha12 < 0.0) {

            alpha12 = alpha12 + two_pi;

        }

        if (alpha12 > two_pi) {

            alpha12 = alpha12 - two_pi;

        }

        double b = a * (1.0 - f);

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

        double U1 = Math.atan(TanU1);

        double sigma1 = Math.atan2(TanU1, Math.cos(alpha12));

        double Sinalpha = Math.cos(U1) * Math.sin(alpha12);

        double cosalpha_sq = 1.0 - (Sinalpha * Sinalpha);

        double u2 = (cosalpha_sq * ((a * a) - (b * b))) / (b * b);

        double A = 1.0 +

                   ((u2 / 16384) * (4096 +

                   (u2 * (-768 + (u2 * (320 - (175 * u2)))))));

        double B = (u2 / 1024) * (256 +

                   (u2 * (-128 + (u2 * (74 - (47 * u2))))));

        // Starting with the approximation

        double sigma = (s / (b * A));

        double last_sigma = (2.0 * sigma) + 2.0; // something impossible

        // Iterate the following three equations

        //  until there is no significant change in sigma

        // two_sigma_m , delta_sigma

        double two_sigma_m = 0;

        while (Math.abs((last_sigma - sigma) / sigma) > 1.0e-9) {

            two_sigma_m = (2 * sigma1) + sigma;

            double delta_sigma = B * Math.sin(sigma) * (Math.cos(two_sigma_m) +

                                 ((B / 4) * (Math.cos(sigma) * ((-1 +

                                 (2 * Math.pow(Math.cos(two_sigma_m), 2))) -

                                 ((B / 6) * Math.cos(two_sigma_m) * (-3 +

                                 (4 * Math.pow(Math.sin(sigma), 2))) * (-3 +

                                 (4 * Math.pow(Math.cos(two_sigma_m), 2))))))));

            last_sigma = sigma;

            sigma = (s / (b * A)) + delta_sigma;

        }

        double phi2 = Math.atan2(((Math.sin(U1) * Math.cos(sigma)) +

                                 (Math.cos(U1) * Math.sin(sigma) * Math.cos(alpha12))),

                                 ((1 - f) * Math.sqrt(Math.pow(Sinalpha, 2) +

                                                      Math.pow((Math.sin(U1) * Math.sin(sigma)) -

                                                               (Math.cos(U1) * Math.cos(sigma) * Math.cos(alpha12)),

                                                               2))));

        double lembda = Math.atan2((Math.sin(sigma) * Math.sin(alpha12)),

                                   ((Math.cos(U1) * Math.cos(sigma)) -

                                   (Math.sin(U1) * Math.sin(sigma) * Math.cos(alpha12))));

        double C = (f / 16) * cosalpha_sq * (4 + (f * (4 - (3 * cosalpha_sq))));

        double omega = lembda -

                       ((1 - C) * f * Sinalpha * (sigma +

                       (C * Math.sin(sigma) * (Math.cos(two_sigma_m) +

                       (C * Math.cos(sigma) * (-1 +

                       (2 * Math.pow(Math.cos(two_sigma_m), 2))))))));

        double lembda2 = lembda1 + omega;

        double alpha21 = Math.atan2(Sinalpha,

                                    ((-Math.sin(U1) * Math.sin(sigma)) +

                                    (Math.cos(U1) * Math.cos(sigma) * Math.cos(alpha12))));

        alpha21 = alpha21 + (two_pi / 2.0);

        if (alpha21 < 0.0) {

            alpha21 = alpha21 + two_pi;

        }

        if (alpha21 > two_pi) {

            alpha21 = alpha21 - two_pi;

        }

        phi2 = (phi2 * 45.0) / piD4;

        lembda2 = (lembda2 * 45.0) / piD4;

        alpha21 = (alpha21 * 45.0) / piD4;

        ret.pt = new Point2D.Double(lembda2, phi2);

        ret.azimut = alpha21;

        return ret;

    }

    /* (non-Javadoc)

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

	 */

    public double surfaceSphere(Point2D pt1, Point2D pt2, Point2D pt3) {

        double sup = -1;

        double A = distanceGeo(pt1, pt2);

        double B = distanceGeo(pt2, pt3);

        double C = distanceGeo(pt3, pt1);

        sup = (((A + B + C) - Math.toRadians(180D)) * Math.PI * proj.getDatum()

                                                                    .getESemiMajorAxis()) / Math.toRadians(180D);

        return sup;

    }

    /*

     * F?rmulas de Vincenty's.

     * (pasadas de http://wegener.mechanik.tu-darmstadt.de/GMT-Help/Archiv/att-8710/Geodetic_py

     * http://www.icsm.gov.au/icsm/gda/gdatm/index.html

     */

    protected class GeoData {

        Point2D pt;

        double azimut;

        double revAzimut;

        double dist;

        public GeoData(double x, double y) {

            pt = new Point2D.Double(x, y);

            azimut = revAzimut = dist = 0;

        }

        public GeoData(double x, double y, double dist, double azi, double rAzi) {

            pt = new Point2D.Double(x, y);

            azimut = azi;

            revAzimut = rAzi;

            this.dist = dist;

        }

    }

}

<!--

    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>Clases relacionadas con el manejo de proyecciones.

</body>

</html>

/* gvSIG. Desktop Geographic Information System.

 *

 * Copyright ? 2007-2015 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 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.

 *

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

 *

 */

public class CoordTransRuntimeException extends RuntimeException {

    /**

     *

     */

    private static final long serialVersionUID = 1348981767734366776L;

    IProjection source;

    IProjection target;

    double x;

    double y;

    /**

     *

     */

    public CoordTransRuntimeException(IProjection source, IProjection target, double x, double y, Throwable cause) {

        super("Error reprojecting point (" + x + "," + y + ") from " + source.getAbrev() + " to " + target.getAbrev()

            + ".", cause);

        this.source = source;

        this.target = target;

        this.x = x;

        this.y = y;

    }

    /**

    *

    */

    public CoordTransRuntimeException(IProjection source, IProjection target, double x, double y) {

        this(source, target, x, y, null);

    }

    /**

     * @return the source

     */

    public IProjection getSource() {

        return source;

    }

    /**

     * @return the target

     */

    public IProjection getTarget() {

        return target;

    }

    /**

     * @return the x

     */

    public double getX() {

        return x;

    }

    /**

     * @return the y

     */

    public double getY() {

        return y;

    }

}

/**

 * 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 java.awt.Color;

import java.awt.Graphics2D;

import java.awt.geom.Point2D;

import java.awt.geom.Rectangle2D;

import org.cresques.geo.ViewPortData;

import org.gvsig.tools.lang.Cloneable;

/**

 *

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

 */

public interface IProjection extends Cloneable{

    /**

     * Return the string representation of projection in the format.

     *

     * Return null if can't export to format

     *

     * @param format

     * @return

     */

    public String export(String format);

    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) throws CoordTransRuntimeException;

    /**

     * 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) throws CoordTransRuntimeException;

    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 static final String FORMAT_WKT = "wkt";

    public static final String FORMAT_PROJ4 = "proj4";

    public static final String FORMAT_WKT_ESRI = "wkt_esri";

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

    /**

     * Devuelve una proyecci?n a partir de una cadena en uno de los formatos admitidos.

     * Devuelve null si no es capaz de crear la proyecci?n.

     *

     * @param format

     * @param value

     * @return

     */

    public IProjection get(String format, String value);

}

*.dfPackage

*.wmf

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