Application: gvSIG desktop

/**

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

	}

}

<?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.impl;

import org.cresques.ProjectionLibrary;

import org.cresques.impl.cts.ProjectionPool;

import org.gvsig.fmap.crs.CRSFactory;

import org.gvsig.tools.library.AbstractLibrary;

import org.gvsig.tools.library.LibraryException;

public class CresquesCtsLibrary extends AbstractLibrary {

    public void doRegistration() {

        registerAsImplementationOf(ProjectionLibrary.class,1);

    }

    protected void doInitialize() throws LibraryException {

        CRSFactory.registerCRSFactory( new ProjectionPool());

    }

    protected void doPostInitialize() throws LibraryException {

        // Nothing to do

    }

}

/**

 * 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.impl.cts.wkt;

import java.util.ArrayList;

/**

 * Generates a WKT from CRS parameters.

 * as, for example in

 * "GEOGCS[\"ETRS89\","+

        "DATUM[\"European_Terrestrial_Reference_System_1989\","+

            "SPHEROID[\"GRS 1980\",6378137,298.257222101,"+

                "AUTHORITY[\"EPSG\",\"7019\"]],"+

            "AUTHORITY[\"EPSG\",\"6258\"]," +

            "TOWGS84[0,0,0,0,0,0,0]],"+

        "PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],"+

        "UNIT[\"degree\",0.01745329251994328, AUTHORITY[\"EPSG\",\"9122\"]],"+

        "AUTHORITY[\"EPSG\",\"4258\"]]";

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

 *

 */

public abstract class WKT {

	static class Authority {

		String name;

		String code;

		public Authority(String n, String c) {

			name = n;

			code = c;

		}

		public String toWKT() {

			return "AUTHORITY[\""+name+"\",\""+code+"\"]";

		}

	}

	static class Parameter {

		protected String name;

		private String value;

		protected Authority authority; // optional

		public Parameter(String n) {

			name = n;

		}

		public Parameter(String n, String v) {

			name = n;

			value = v;

		}

		public Parameter(String n, Authority a) {

			name = n;

			authority = a;

		}

		public String toWKT() {

			String txt = "PARAMETER[\""+name+"\","+value;

			if (authority != null)

				txt += ","+authority.toWKT();

			return txt + "]";

		}

	}

	static class Spheroid {

		String name;

		double semiMajor;

		double inverseFlattening;

		Authority authority = null; // optional

		public Spheroid(String n, double s, double f) {

			name = n;

			semiMajor = s;

			inverseFlattening = f;

		}

		public Spheroid(String n, double s, double f, Authority a) {

			name = n;

			semiMajor = s;

			inverseFlattening = f;

			authority = a;

		}

		public String toWKT() {

			String txt = "SPHEROID[\""+name+"\","+

				semiMajor+","+inverseFlattening;

			if (authority != null)

				txt += ","+authority.toWKT();

			return txt + "]";

		}

	}

	static class ToWGS84 {

		double dx = 0;

		double dy = 0;

		double dz = 0;

		double ex = 0;

		double ey = 0;

		double ez = 0;

		double ppm = 0;

		public ToWGS84() {

		}

		public ToWGS84(double dx, double dy, double dz, double ex, double ey, double ez, double ppm) {

			this.dx = dx;

			this.dy = dy;

			this.dz = dz;

			this.ex = ex;

			this.ey = ey;

			this.ez = ez;

			this.ppm = ppm;

		}

		public String toWKT() {

			return "TOWGS84["+dx+","+dy+","+dz+","+ex+","+ey+","+ez+","+ppm+"]";

		}

	}

	static class Datum {

		String name;

		Spheroid spheroid;

		ToWGS84 towgs84 = new ToWGS84(); // optional

		Authority authority; // optional

		public Datum(String n, Spheroid s) {

			init (n, s, new ToWGS84(), null);

		}

		public Datum(String n, Spheroid s, ToWGS84 t) {

			init (n, s, t, null);

		}

		public Datum(String n, Spheroid s, Authority a) {

			init (n, s, new ToWGS84(), a);

		}

		public Datum(String n, Spheroid s, ToWGS84 t, Authority a) {

			init (n, s, t, a);

		}

		private void init(String n, Spheroid s,

			ToWGS84 t, Authority a) {

			name = n;

			spheroid = s;

			towgs84 = t;

			authority = a;

		}

		public String toWKT() {

			String txt = "DATUM[\""+name+"\","+spheroid.toWKT();

			if (towgs84 != null)

				txt += ","+towgs84.toWKT();

			if (authority != null)

				txt += ","+authority.toWKT();

			return txt+"]";

		}

	}

	static class PrimeM {

		String name;

		double longitude;

		Authority authority = null; // optional

		public PrimeM(String n, double l) {

			name = n;

			longitude = l;

		}

		public PrimeM(String n, double l, Authority a) {

			name = n;

			longitude = l;

			authority = a;

		}

		public String toWKT() {

			String txt = "PRIMEM[\""+name+"\","+longitude;

			if (authority != null)

				txt += ","+authority.toWKT();

			return txt + "]";

		}

	}

	static class Unit {

		String name;

		double conversionFactor;

		Authority authority; // optional

		public Unit(String n, double l) {

			name = n;

			conversionFactor = l;

		}

		public Unit(String n, double l, Authority a) {

			name = n;

			conversionFactor = l;

			authority = a;

		}

		public String toWKT() {

			String txt = "UNIT[\""+name+"\","+conversionFactor;

			if (authority != null)

				txt += ","+authority.toWKT();

			return txt + "]";

		}

	}

	static class Orientation {

		public final static int NORTH = 0;

		public final static int SOUTH = 1;

		public final static int EAST = 2;

		public final static int WEST = 3;

		public final static int UP = 4;

		public final static int DOWN = 5;

		public final static int OTHER = 6;

		protected final String [] txt = {

			"NORTH","SOUTH","EAST","WEST","UP","DOWN","OTHER"

		};

		protected int to = -1;

		public Orientation(int to) {

			this.to = to;

		}

	}

	static class Axis extends Orientation {

		String name;

		public Axis(String name, int to) {

			super(to);

			this.name = name;

		}

		public String toWKT() {

			return "AXIS[\""+name+"\","+txt[to]+ "]";

		}

	}

	static class TwinAxes {

		Axis axe1;

		Axis axe2;

		public TwinAxes(Axis a1, Axis a2) {

			axe1 = a1;

			axe2 = a2;

		}

		public String toWKT() {

			return axe1.toWKT()+","+axe2.toWKT();

		}

	}

	static class GeogCS extends WKT {

		String name;

		Datum datum;

		PrimeM primeMeridian;

		Unit unit;

		TwinAxes axes = null; // optional

		Authority authority = null; // optional

		public GeogCS(String n, Datum d, PrimeM p, Unit u) {

			init(n, d, p, u, null, null);

		}

		public GeogCS(String n, Datum d, PrimeM p, Unit u,

				Authority a) {

			init(n, d, p, u, null, a);

		}

		public GeogCS(String n, Datum d, PrimeM p, Unit u,

				TwinAxes t) {

			init(n, d, p, u, t, null);

		}

		public GeogCS(String n, Datum d, PrimeM p, Unit u,

				TwinAxes t, Authority a) {

			init(n, d, p, u, t, a);

		}

		private void init(String n, Datum d, PrimeM p, Unit u,

			TwinAxes t, Authority a) {

			name = n;

			datum = d;

			primeMeridian = p;

			unit = u;

			axes = t;

			authority = a;

		}

		public String toWKT() {

			String txt = "GEOGCS[\""+name+"\","+datum.toWKT()+

				","+primeMeridian.toWKT()+","+unit.toWKT();

			if (axes != null)

				txt += ","+axes.toWKT();

			if (authority != null)

				txt += ","+authority.toWKT();

			return txt+"]";

		}

	}

	static class Projection extends Parameter {

		public Projection(String n) {

			super(n);

		}

		public Projection(String n, Authority a) {

			super(n,a);

		}

		public String toWKT() {

			String txt = "PROJECTION[\""+name+"\"";

			if (authority != null)

				txt += ","+authority.toWKT();

			return txt + "]";

		}

	}

	static class ProjCS extends Parameter {

		GeogCS geogcs;

		Projection proj;

		ArrayList params = new ArrayList();

		Unit unit;

		TwinAxes axes = null; // optional

		public ProjCS(String n, GeogCS cs, Projection prj, Unit u ) {

			super(n);

			Parameter [] p = {};

			init(n, cs, prj, p, u, null);

		}

		public ProjCS(String n, GeogCS cs, Projection prj,

				Parameter [] p, Unit u ) {

			super(n);

			init(n, cs, prj, p, u, null);

		}

		public ProjCS(String n, GeogCS cs, Projection prj,

				Unit u, TwinAxes t, Authority a ) {

			super(n, a);

			Parameter [] p = {};

			init(n, cs, prj, p, u, t);

		}

		public ProjCS(String n, GeogCS cs, Projection prj,

				Parameter [] p, Unit u, TwinAxes t, Authority a ) {

			super(n, a);

			init(n, cs, prj, p, u, t);

		}

		private void init(String n, GeogCS cs, Projection prj,

			Parameter [] p, Unit u, TwinAxes t) {

			geogcs = cs;

			proj = prj;

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

				params.add(p[i]);

			unit = u;

			axes = t;

		}

		public void add(Parameter p) {

			params.add(p);

		}

		public String toWKT() {

			String txt = "PROJCS[\""+name+"\","+geogcs.toWKT()+

				","+proj.toWKT();

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

				txt += ","+((Parameter)params.get(i)).toWKT();

			txt += ","+unit.toWKT();

			if (axes != null)

				txt += ","+axes.toWKT();

			if (authority != null)

				txt += ","+authority.toWKT();

			return txt+"]";

		}

	}

	abstract public String toWKT();

}

/**

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

import java.util.ArrayList;

import java.util.Iterator;

import java.util.Map;

import java.util.TreeMap;

import org.cresques.cts.ICRSFactory;

import org.cresques.cts.IProjection;

import org.cresques.impl.cts.gt2.CSDatum;

import org.cresques.impl.cts.gt2.CSGaussPt;

import org.cresques.impl.cts.gt2.CSLambertCC;

import org.cresques.impl.cts.gt2.CSMercator;

import org.cresques.impl.cts.gt2.CSUTM;

import org.cresques.impl.cts.gt2.CoordSys;

/**

 * Pool de proyeccions (cs+datum) conocidas.

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

 */

public class ProjectionPool implements ICRSFactory {

    static TreeMap data = null;

    static {

        CoordSys cs = null;

        data = new TreeMap();

        cs = (new CSUTM(CSDatum.wgs84, 30)).toGeo();

        cs.setAbrev("EPSG:4326"); // WGS84 (World Geodesic Datum)

        data.put(cs.getAbrev(), cs);

        data.put("CRS:84", cs); // CRS:84 = EPSG:4326

        cs = (new CSUTM(CSDatum.ed50, 30)).toGeo();

        cs.setAbrev("EPSG:4230"); // Datum Europeu Internacional ED50

        data.put(cs.getAbrev(), cs);

        cs = (new CSUTM(CSDatum.d73, 30)).toGeo();

        cs.setAbrev("EPSG:4274"); // Datum 73 de Lisboa

        data.put(cs.getAbrev(), cs);

        cs = (new CSUTM(CSDatum.nad27, 30)).toGeo();

        cs.setAbrev("EPSG:4267"); // NAD 27

        data.put(cs.getAbrev(), cs);

        cs = (new CSUTM(CSDatum.nad83, 30)).toGeo();

        cs.setAbrev("EPSG:4269"); // NAD 83

        data.put(cs.getAbrev(), cs);

        cs = (new CSUTM(CSDatum.lomaQuintana, 30)).toGeo();

        cs.setAbrev("EPSG:4288"); // PSAD 56 'Loma Quintana'

        data.put(cs.getAbrev(), cs);

        cs = (new CSUTM(CSDatum.laCanoa, 30)).toGeo();

        cs.setAbrev("EPSG:4247"); // PSAD 56 'Loma Quintana'

        data.put(cs.getAbrev(), cs);

        cs = (new CSUTM(CSDatum.ntfParis, 30)).toGeo();

        cs.setAbrev("EPSG:4807"); // NTF Paris

        data.put(cs.getAbrev(), cs);

        cs = (new CSUTM(CSDatum.etrs89, 30)).toGeo();

        cs.setAbrev("EPSG:4258"); // ETRS 89

        data.put(cs.getAbrev(), cs);

        for (int i = 1; i <= 60; i++) {

            String huso = Integer.toString(i);

            if (i < 10) {

                huso = "0" + huso;

            }

            cs = new CSUTM(CSDatum.wgs84, i);

            cs.setAbrev("EPSG:326" + huso);

            data.put(cs.getAbrev(), cs);

            cs = new CSUTM(CSDatum.ed50, i);

            cs.setAbrev("EPSG:230" + huso);

            data.put(cs.getAbrev(), cs);

            if (i>2 && i<=23) {

                cs = new CSUTM(CSDatum.nad27, i);

                cs.setAbrev("EPSG:267" + huso);

                data.put(cs.getAbrev(), cs);

                cs = new CSUTM(CSDatum.nad83, i);

                cs.setAbrev("EPSG:269" + huso);

                data.put(cs.getAbrev(), cs);

            }

            if (i>27 && i<39) {

                cs = new CSUTM(CSDatum.etrs89, i);

                cs.setAbrev("EPSG:258" + huso);

                data.put(cs.getAbrev(), cs);

            }

        }

        cs = CSGaussPt.hgd73;

        cs.setAbrev("EPSG:27492"); // Projec??o Gauss do Datum 73 de Lisboa (no EPSG found)

        data.put(cs.getAbrev(), cs);

        cs = new CSUTM(CSDatum.d73, 29);

        cs.setAbrev("EPSG:27429"); // Projec??o Gauss do Datum 73 de Lisboa (no EPSG found)

        data.put(cs.getAbrev(), cs);

        for (int i = 18; i <= 21; i++) {

            String huso = Integer.toString(i);

            if (i < 10) {

                huso = "0" + huso;

            }

            cs = new CSUTM(CSDatum.laCanoa, i);

            cs.setAbrev("EPSG:247" + huso);

            data.put(cs.getAbrev(), cs);

        }

        for (int i = 16; i <= 22; i++) {

            String huso = Integer.toString(i);

            // Psad56 Ecuador

            cs = new CSUTM(CSDatum.lomaQuintana, i);

            cs.setAbrev("EPSG:288" + huso);

            data.put(cs.getAbrev(), cs);

            cs = new CSUTM(CSDatum.lomaQuintana, i, "S");

            cs.setAbrev("EPSG:288" + (i+60) );

            data.put(cs.getAbrev(), cs);

        }

        //		cs = new CSLambertCC(CSDatum.nad27, -105D, 49D, 49D, 77D, 0, 0);

        //		cs.setAbrev("LCCCan");

        //		data.put(cs.getAbrev(), cs);

        /* Para el server WMS de canad?:

         * EPSG:42101

         * EPSG:42304

         * EPSG:4269

         */

        /*

         * 42101,PROJCS["WGS 84 / LCC Canada",

         * GEOGCS["WGS 84",DATUM["WGS_1984",

         * SPHEROID["WGS_1984",6378137,298.257223563]],

         * PRIMEM["Greenwich",0],UNIT["Decimal_Degree",0.0174532925199433]],

         *

         * PROJECTION["Lambert_Conformal_Conic_2SP"],

         * PARAMETER["central_meridian",-95.0],

         * PARAMETER["latitude_of_origin",0],

         * PARAMETER["standard_parallel_1",49.0],

         * PARAMETER["standard_parallel_2",77.0],

         * PARAMETER["false_easting",0.0],

         * PARAMETER["false_northing",-8000000.0],

         * UNIT["Meter",1],AUTHORITY["EPSG","42101"]]

         */

        cs = new CSLambertCC(CSDatum.wgs84, -95, 0, 49, 77, 0, -8000000.0);

        cs.setAbrev("EPSG:42101");

        data.put(cs.getAbrev(), cs);

        /* 42304,PROJCS["NAD83 / NRCan LCC Canada",

         * GEOGCS["NAD83",DATUM["North_American_Datum_1983",

         * SPHEROID["GRS_1980",6378137,298.257222101]],

         * PRIMEM["Greenwich",0],

         * UNIT["Decimal_Degree",0.0174532925199433]],

         *

         * PROJECTION["Lambert_Conformal_Conic_2SP"],

         * PARAMETER["central_meridian",-95.0],

         * PARAMETER["latitude_of_origin",49.0],

         * PARAMETER["standard_parallel_1",49.0],

         * PARAMETER["standard_parallel_2",77.0],

         * PARAMETER["false_easting",0.0],

         * PARAMETER["false_northing",0.0],

         * UNIT["Meter",1],AUTHORITY["EPSG","42304"]]

         */

        cs = new CSLambertCC(CSDatum.nad83, -95, 49, 49, 77, 0, 0);

        cs.setAbrev("EPSG:42304");

        data.put(cs.getAbrev(), cs);

        /*

         * EPSG:26915 - NAD83 / UTM zone 15N

         * EPSG:31466 - Gau?-Kr?ger band 2

         * EPSG:31467 - Gau?-Kr?ger band 3

         * EPSG:4314  - DHDN

         */

        /*

         * 27572=PROJCS["NTF (Paris) / Lambert zone II",

             GEOGCS["NTF (Paris)",

                    DATUM["Nouvelle_Triangulation_Francaise_Paris",

                          SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936269,

                                   AUTHORITY["EPSG","7011"]],

                          TOWGS84[-168,-60,320,0,0,0,0],

                          AUTHORITY["EPSG","6807"]],

                    PRIMEM["Paris",2.33722917,AUTHORITY["EPSG","8903"]],

                    UNIT["grad",0.01570796326794897,AUTHORITY["EPSG","9105"]],

                    AUTHORITY["EPSG","4807"]],

             PROJECTION["Lambert_Conformal_Conic_1SP"],

                        PARAMETER["latitude_of_origin",52],

                        PARAMETER["central_meridian",0],

                        PARAMETER["scale_factor",0.99987742],

                        PARAMETER["false_easting",600000],

                        PARAMETER["false_northing",2200000],

                        UNIT["metre",1,AUTHORITY["EPSG","9001"]],

                        AUTHORITY["EPSG","27572"]]

         */

        cs = new CSLambertCC(CSDatum.ntfParis, 0, 46.79999999999995, 0.99987742, 600000, 2200000);

        cs.setAbrev("EPSG:27572");

        data.put(cs.getAbrev(), cs);

        cs = new CSLambertCC(CSDatum.ntfParis, 0, 52, 0.99987742, 600000, 2200000);

        cs.setAbrev("EPSG:27582");

        data.put(cs.getAbrev(), cs);

        /*

         * # RGF93

         * <4171> +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs  <>

         * # RGF93 / Lambert-93

         */

        cs = (new CSUTM(CSDatum.rgf93, 30)).toGeo();

        cs.setAbrev("EPSG:4171"); // NTF Paris

        data.put(cs.getAbrev(), cs);

        /*

         *  <2154> +proj=lcc +lat_1=49 +lat_2=44 +lat_0=46.5 +lon_0=3

         *  +x_0=700000 +y_0=6600000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0

         *  +units=m +no_defs  <>

         */

        cs = new CSLambertCC(CSDatum.rgf93, 3.0, 46.5, 49.0, 44.0, 700000.0, 6600000.0);

        cs.setAbrev("EPSG:2154");

        data.put(cs.getAbrev(), cs);

        cs = new CSMercator(CSDatum.wgs84);

        cs.setAbrev("EPSG:54004");

        data.put(cs.getAbrev(), cs);

        cs.setAbrev("EPSG:9804");

        data.put(cs.getAbrev(), cs);

        // Lo que faltaba: ?planetas!

        cs = (new CSUTM(CSDatum.moon, 30)).toGeo();

        cs.setAbrev("IAU2000:30100"); // Moon

        data.put(cs.getAbrev(), cs);

        cs = (new CSUTM(CSDatum.mars, 30)).toGeo();

        cs.setAbrev("IAU2000:49900"); // Mars

        data.put(cs.getAbrev(), cs);

        /*

         * CRSs argentinos.

         * coordenadas geograficas

			PosGAr 4172

			PosGAr98 4190

			coordenadas proyectadas

			POSGAR 94/Argentina 1 22191

			POSGAR 94/Argentina 2 22192

			POSGAR 94/Argentina 3 22193

			POSGAR 94/Argentina 4 22194

			POSGAR 94/Argentina 5 22195

			POSGAR 94/Argentina 6 22196

			POSGAR 94/Argentina 7 22197

			POSGAR 98/Argentina 1 22181

			POSGAR 98/Argentina 2 22182

			POSGAR 98/Argentina 3 22183

			POSGAR 98/Argentina 4 22184

			POSGAR 98/Argentina 5 22185

			POSGAR 98/Argentina 6 22186

			POSGAR 98/Argentina 7 22187

			4221 GEOGCS["GCS_Campo_Inchauspe",DATUM["D_Campo_Inchauspe",SPHEROID["International_1924",6378388,297]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]

			22191	EPSG	22191	PROJCS["Argentina_Zone_1",GEOGCS["GCS_Campo_Inchauspe",DATUM["D_Campo_Inchauspe",SPHEROID["International_1924",6378388,297]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",1500000],PARAMETER["False_Northing",0],PARAMETER["Central_Meridian",-72],PARAMETER["Scale_Factor",1],PARAMETER["Latitude_Of_Origin",-90],UNIT["Meter",1]]	+proj=tmerc +lat_0=-90 +lon_0=-72 +k=1.000000 +x_0=1500000 +y_0=0 +ellps=intl +units=m

			22192	EPSG	22192	PROJCS["Argentina_Zone_2",GEOGCS["GCS_Campo_Inchauspe",DATUM["D_Campo_Inchauspe",SPHEROID["International_1924",6378388,297]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",2500000],PARAMETER["False_Northing",0],PARAMETER["Central_Meridian",-69],PARAMETER["Scale_Factor",1],PARAMETER["Latitude_Of_Origin",-90],UNIT["Meter",1]]	+proj=tmerc +lat_0=-90 +lon_0=-69 +k=1.000000 +x_0=2500000 +y_0=0 +ellps=intl +units=m

			22193	EPSG	22193	PROJCS["Argentina_Zone_3",GEOGCS["GCS_Campo_Inchauspe",DATUM["D_Campo_Inchauspe",SPHEROID["International_1924",6378388,297]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",3500000],PARAMETER["False_Northing",0],PARAMETER["Central_Meridian",-66],PARAMETER["Scale_Factor",1],PARAMETER["Latitude_Of_Origin",-90],UNIT["Meter",1]]	+proj=tmerc +lat_0=-90 +lon_0=-66 +k=1.000000 +x_0=3500000 +y_0=0 +ellps=intl +units=m

			22194	EPSG	22194	PROJCS["Argentina_Zone_4",GEOGCS["GCS_Campo_Inchauspe",DATUM["D_Campo_Inchauspe",SPHEROID["International_1924",6378388,297]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",4500000],PARAMETER["False_Northing",0],PARAMETER["Central_Meridian",-63],PARAMETER["Scale_Factor",1],PARAMETER["Latitude_Of_Origin",-90],UNIT["Meter",1]]	+proj=tmerc +lat_0=-90 +lon_0=-63 +k=1.000000 +x_0=4500000 +y_0=0 +ellps=intl +units=m

			22195	EPSG	22195	PROJCS["Argentina_Zone_5",GEOGCS["GCS_Campo_Inchauspe",DATUM["D_Campo_Inchauspe",SPHEROID["International_1924",6378388,297]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",5500000],PARAMETER["False_Northing",0],PARAMETER["Central_Meridian",-60],PARAMETER["Scale_Factor",1],PARAMETER["Latitude_Of_Origin",-90],UNIT["Meter",1]]	+proj=tmerc +lat_0=-90 +lon_0=-60 +k=1.000000 +x_0=5500000 +y_0=0 +ellps=intl +units=m

			22196	EPSG	22196	PROJCS["Argentina_Zone_6",GEOGCS["GCS_Campo_Inchauspe",DATUM["D_Campo_Inchauspe",SPHEROID["International_1924",6378388,297]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",6500000],PARAMETER["False_Northing",0],PARAMETER["Central_Meridian",-57],PARAMETER["Scale_Factor",1],PARAMETER["Latitude_Of_Origin",-90],UNIT["Meter",1]]	+proj=tmerc +lat_0=-90 +lon_0=-57 +k=1.000000 +x_0=6500000 +y_0=0 +ellps=intl +units=m

			22197	EPSG	22197	PROJCS["Argentina_Zone_7",GEOGCS["GCS_Campo_Inchauspe",DATUM["D_Campo_Inchauspe",SPHEROID["International_1924",6378388,297]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",7500000],PARAMETER["False_Northing",0],PARAMETER["Central_Meridian",-54],PARAMETER["Scale_Factor",1],PARAMETER["Latitude_Of_Origin",-90],UNIT["Meter",1]]	+proj=tmerc +lat_0=-90 +lon_0=-54 +k=1.000000 +x_0=7500000 +y_0=0 +ellps=intl +units=m

         */

        cs = new CoordSys(CSDatum.posgar);

        cs.setAbrev("EPSG:4172"); // Posgar

        data.put(cs.getAbrev(), cs);

        cs = new CoordSys(

			"GEOGCS[\"GCS_Campo_Inchauspe\"," +

				"DATUM[\"D_Campo_Inchauspe\"," +

					"SPHEROID[\"International_1924\",6378388,297],"+

					"TOWGS84[0,0,0,0,0,0,0]]," +

				"PRIMEM[\"Greenwich\",0]," +

				"UNIT[\"Degree\",0.017453292519943295]]");

        cs.setAbrev("EPSG:4221"); // Campo Inchauspe

        data.put(cs.getAbrev(), cs);

       for (int i=1; i<=7; i++) {

	        cs = new CoordSys(

	        	"PROJCS[\"Argentina_Zone_"+i+"\"," +

	        		"GEOGCS[\"GCS_Campo_Inchauspe\"," +

	        			"DATUM[\"D_Campo_Inchauspe\"," +

		        			"SPHEROID[\"International_1924\",6378388,297],"+

		        			"TOWGS84[0,0,0,0,0,0,0]]," +

		        		"PRIMEM[\"Greenwich\",0]," +

		        		"UNIT[\"Degree\",0.017453292519943295]]," +

		        	"PROJECTION[\"Transverse_Mercator\"]," +

		        	"PARAMETER[\"False_Easting\",1500000]," +

		        	"PARAMETER[\"False_Northing\",0]," +

		        	"PARAMETER[\"Central_Meridian\","+(-75+3*i)+"]," +

		        	"PARAMETER[\"Scale_Factor\",1]," +

		        	"PARAMETER[\"Latitude_Of_Origin\",-90]," +

		        	"UNIT[\"Meter\",1]]");

	        cs.setAbrev("EPSG:2219"+i); // Posgar

	        data.put(cs.getAbrev(), cs);

        }

       /*

        * pendiente de a?adir:

        *

        * EPSG 3003: Montemario / Italy Zone 1 - Pendiente de a?adir

      	* EPSG 4149: CH1903 - Pendiente de a?adir

      	* Ecuador:

      	*   PSAD56, Geo, UTM 16/17S,16/17N

      	*/

    }

    /**

     * Mete una nueva proyeccion en la Pool.

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

     * @param proj Proyeccion

     */

    public static void add(String name, IProjection proj) {

        data.put(name, proj);

    }

    /**

     * 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) {

        IProjection proj = null;

        if (ProjectionPool.data.containsKey(name)) {

            proj = (IProjection) ProjectionPool.data.get(name);

        } else {

        	// Consultation to remote EPSG database

        	// if (right)

        	//    buil new IProjection from GML

        	// else

            System.err.println("ProjectionPool: Key '" + name + "' not set.");

        }

        return proj;

    }

    public static Iterator iterator() {

        ArrayList projs = new ArrayList();

        Iterator iter = data.entrySet().iterator();

        while (iter.hasNext()) {

            projs.add(((Map.Entry) iter.next()).getValue());

        }

        return projs.iterator();

    }

	public boolean doesRigurousTransformations() {

		return false;

	}

    /* (non-Javadoc)

     * @see org.cresques.cts.ICRSFactory#get(java.lang.String, java.lang.String)

     */

    @Override

    public IProjection get(String format, String value) {

        // TODO Auto-generated method stub

        return null;

    }

}

/**

 * 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.impl.cts.gt2;

import java.awt.Color;

import java.awt.Graphics2D;

import java.awt.geom.Point2D;

import java.awt.geom.Rectangle2D;

import org.cresques.cts.ICoordTrans;

import org.cresques.cts.IDatum;

import org.cresques.cts.IProjection;

import org.cresques.geo.ViewPortData;

import org.geotools.cs.CoordinateSystem;

import org.geotools.cs.CoordinateSystemFactory;

import org.geotools.cs.GeographicCoordinateSystem;

import org.geotools.cs.ProjectedCoordinateSystem;

import org.opengis.referencing.FactoryException;

import org.gvsig.fmap.crs.CRSFactory;

/**

 * Sistema de Coordenadas (Proyecci?n).

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

 */

public class CoordSys implements IProjection {

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

    protected CoordinateSystemFactory csFactory = CoordinateSystemFactory.getDefault();

    protected CSDatum datum = null;

    protected GeographicCoordinateSystem geogCS = null;

    protected ProjectedCoordinateSystem projCS = null;

    protected String abrev = "";

    Color gridColor = basicGridColor;

    public CoordSys(CSDatum datum) {

        this.datum = datum;

        this.geogCS = new GeographicCoordinateSystem(datum.getName(null), datum.getDatum());

	}

    public CoordSys(String wkt) {

    	try {

        	//((GeographicCoordinateSystem)

    		CoordinateSystem cs = CoordinateSystemFactory.getDefault().createFromWKT(wkt);

//			 ).getHorizontalDatum();

    		if (cs instanceof GeographicCoordinateSystem)

    			geogCS = (GeographicCoordinateSystem) cs;

    		if (cs instanceof ProjectedCoordinateSystem) {

    			projCS = (ProjectedCoordinateSystem) cs;

    			geogCS = projCS.getGeographicCoordinateSystem();

    		}

    		datum = new CSDatum(geogCS.getHorizontalDatum());

    	} catch (FactoryException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

    }

    /**

     * Crea un nuevo CoordSys geogr?fico a partir de uno proyectado.

     * Si el actual es geogr?fico retorna el mismo.

     * @return

     */

    public CoordSys toGeo() {

        CoordSys coordSys = new CoordSys((CSDatum) getDatum());

        if (geogCS != null) {

            coordSys.geogCS = this.geogCS;

        } else {

            coordSys.geogCS = projCS.getGeographicCoordinateSystem();

        }

        return coordSys;

    }

    public IDatum getDatum() {

        return datum;

    }

    public CoordinateSystem getCS() {

        if (projCS != null) {

            return projCS;

        }

        return geogCS;

    }

    public CoordinateSystem getGeogCS() {

        if (geogCS != null) {

            return geogCS;

        }

        return projCS.getGeographicCoordinateSystem();

    }

    /* (no Javadoc)

     * @see org.cresques.cts.IProjection#createPoint(double, double)

     */

    public Point2D createPoint(double x, double y) {

        return new Point2D.Double(x, y);

    }

    public String toString() {

        if (projCS != null) {

            return projCS.toString();

        }

        return geogCS.toString();

    }

    public void setAbrev(String abrev) {

        this.abrev = abrev;

    }

    public String getAbrev() {

        return abrev;

    }

    /* (no Javadoc)

     * @see org.cresques.cts.IProjection#drawGrid(java.awt.Graphics2D, org.cresques.geo.ViewPortData)

     */

    public void drawGrid(Graphics2D g, ViewPortData vp) {

        // TODO Ap?ndice de m?todo generado autom?ticamente

    }

    public void setGridColor(Color c) {

        gridColor = c;

    }

    public Color getGridColor() {

        return gridColor;

    }

    /* (no Javadoc)

     * @see org.cresques.cts.IProjection#toGeo(java.awt.geom.Point2D)

     */

    public Point2D toGeo(Point2D pt) {

    //TODO VCN Esto si no lo comento no me trasforma el punto en coordenadas geogr?ficas.

//        if (getGeogCS() == geogCS) {

//            return pt;

//        } else {

            CoordTrans ct = new CoordTrans(this, toGeo());

            return ct.convert(pt, null);

//        }

    }

    /* (no Javadoc)

     * @see org.cresques.cts.IProjection#fromGeo(java.awt.geom.Point2D, java.awt.geom.Point2D)

     */

    public Point2D fromGeo(Point2D gPt, Point2D mPt) {

        // TODO Ap?ndice de m?todo generado autom?ticamente

        return null;

    }

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

        double scale = 0D;

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

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

            // 1852.0 metros x minuto de meridiano

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

                    w; // pixels

        }

        return scale;

    }

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

    	double w =0;

		double h =0;

		double wExtent =0;

		double hExtent =0;

    	if (projCS!=null) {

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

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

			wExtent =w * scale*distanceUnits/ mapUnits;

			hExtent =h * scale*distanceUnits/ mapUnits;

		}else {

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

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

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

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

		}

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

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

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

    	return  rec;

    }

    public boolean isProjected() {

    	return projCS != null;

    }

	public ICoordTrans getCT(IProjection dest) {

		return new CoordTrans(this, (CoordSys) dest);

	}

	public String getFullCode() {

		return getAbrev();

	}

	public Object clone() throws CloneNotSupportedException {

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

	    }

    /* (non-Javadoc)

     * @see org.cresques.cts.IProjection#export(java.lang.String)

     */

    @Override

    public String export(String format) {

        // TODO Auto-generated method stub

        return null;

    }

}

/**

 * 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.impl.cts.gt2;

import org.cresques.cts.UTM;

import org.geotools.cs.AxisInfo;

import org.geotools.cs.Projection;

import org.geotools.units.Unit;

import org.gvsig.fmap.crs.CRSFactory;

import org.opengis.referencing.FactoryException;

/**

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

 */

public class CSUTM extends CoordSys implements UTM {

    public CSUTM(CSDatum datum, int zone) {

        super(datum);

        init(datum, zone, "N");

    }

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

        super(datum);

        init(datum, zone, ns);

    }

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

        Unit linearUnit = Unit.METRE;

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

        params.setParameter("semi_major",

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

        params.setParameter("semi_minor",

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

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

        params.setParameter("latitude_of_origin", 0.0);

        params.setParameter("scale_factor", 0.9996);

        params.setParameter("false_easting", 500000.0);

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

        	params.setParameter("false_northing", 10000000.0);

        else

        	params.setParameter("false_northing", 0.0);

        try {

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

                                                               "Transverse_Mercator",

                                                               params);

            projCS = csFactory.createProjectedCoordinateSystem(projection.getName()

                                                                         .toString(),

                                                               geogCS,

                                                               projection,

                                                               linearUnit,

                                                               AxisInfo.X,

                                                               AxisInfo.Y);

        } catch (FactoryException e) {

            // TODO Bloque catch generado autom?ticamente

            e.printStackTrace();

        }

    }

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

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

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

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

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

                    w; // pixels

        }

        return scale;

    }

    public String toString() {

        return projCS.toString();

    }

    public Object clone() throws CloneNotSupportedException {

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

     }

}

/**

 * 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.impl.cts.gt2;

import org.geotools.cs.AxisInfo;

import org.geotools.cs.GeographicCoordinateSystem;

import org.geotools.cs.Projection;

import org.geotools.measure.Angle;

import org.geotools.units.Unit;

import org.opengis.referencing.FactoryException;

/**

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

 *

 */

public class CSGaussPt extends CoordSys {

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

    public CSGaussPt(CSDatum datum) {

        super(datum);

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

                                                datum.getDatum());

        Unit linearUnit = Unit.METRE;

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

        params.setParameter("semi_major",

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

        params.setParameter("semi_minor",

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

        params.setParameter("central_meridian",

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

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

        params.setParameter("latitude_of_origin",

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

        params.setParameter("scale_factor", 1.0);

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

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

        try {

            Projection projection = csFactory.createProjection("GAUSS",

                                                               "Oblique_Stereographic",

                                                               params);

            projCS = csFactory.createProjectedCoordinateSystem(projection.getName()

                                                                         .toString(),

                                                               geogCS,

                                                               projection,

                                                               linearUnit,

                                                               AxisInfo.X,

                                                               AxisInfo.Y);

        } catch (FactoryException e) {

            // TODO Bloque catch generado autom?ticamente

            e.printStackTrace();

        }

    }

    public String toString() {

        return projCS.toString();

    }

}

/**

 * 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.impl.cts.gt2;

import java.awt.geom.Point2D;

import java.awt.geom.Rectangle2D;

import org.cresques.cts.ICoordTrans;

import org.cresques.cts.IProjection;

import org.geotools.ct.CannotCreateTransformException;

import org.geotools.ct.CoordinateTransformation;

import org.geotools.ct.CoordinateTransformationFactory;

import org.geotools.ct.MathTransform;

import org.geotools.pt.CoordinatePoint;

import org.opengis.referencing.operation.TransformException;

//import org.geotools.pt.MismatchedDimensionException;

/**

 * Transforma coordenadas entre dos sistemas

 * @see org.creques.cts.CoordSys

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

 */

class CoordTrans implements ICoordTrans {

    private CoordinateTransformationFactory trFactory = CoordinateTransformationFactory.getDefault();

    private CoordinateTransformation tr = null;

    private MathTransform mt = null;

    private MathTransform mt2 = null;

    private MathTransform mt3 = null;

    private MathTransform mtDatum = null;

    private CoordSys from = null;

    private CoordSys to = null;

    private ICoordTrans invertedCT = null;

    public CoordTrans(CoordSys from, CoordSys to) {

        this.from = from;

        this.to = to;

        // Si los dos CoordSys son proyectados, entonces hay

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

        // Si hay cambio de datum son 3 las transformaciones

        try {