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

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

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

                                                           to.toGeo().getCS());

                mtDatum = tr.getMathTransform();

            }

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

                CoordSys geogcs = from.toGeo();

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

                                                           geogcs.getCS());

                mt = tr.getMathTransform();

                if (mtDatum != null) {

                    mt2 = mtDatum;

                }

                geogcs = to.toGeo();

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

                                                           to.getCS());

                if (mt2 == null) {

                    mt2 = tr.getMathTransform();

                } else {

                    mt3 = tr.getMathTransform();

                }

            } else {

                if (from.projCS == null) {

                    mt = mtDatum;

                }

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

                                                           to.getCS());

                if (mt == null) {

                    mt = tr.getMathTransform();

                    if (mtDatum != null) {

                        mt2 = mtDatum;

                    }

                } else {

                    mt2 = tr.getMathTransform();

                }

            }

        } catch (CannotCreateTransformException e) {

            // TODO Bloque catch generado autom?ticamente

            e.printStackTrace();

        }

    }

    public IProjection getPOrig() {

        return from;

    }

    public IProjection getPDest() {

        return to;

    }

    public ICoordTrans getInverted() {

        if (invertedCT == null)

            invertedCT = new CoordTrans(to, from);

        return invertedCT;

    }

    public Point2D convert(Point2D ptOrig, Point2D ptDest) {

        CoordinatePoint pt1 = new CoordinatePoint(ptOrig);

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

        ptDest = null;

        try {

            mt.transform(pt1, pt2);

            ptDest = pt2.toPoint2D();

            if (mt2 != null) {

                mt2.transform(pt2, pt1);

                ptDest = pt1.toPoint2D();

                if (mt3 != null) {

                    mt3.transform(pt1, pt2);

                    ptDest = pt2.toPoint2D();

                }

            }

            /*} catch (MismatchedDimensionException e) {

                    // TODO Bloque catch generado autom?ticamente

                    e.printStackTrace();

            */

        } catch (TransformException e) {

            // TODO Bloque catch generado autom?ticamente

            e.printStackTrace();

        }

        return ptDest;

    }

    public String toString() {

        return tr.toString();

    }

    /* (non-Javadoc)

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

     */

    public Rectangle2D convert(Rectangle2D rect) {

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

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

        Point2D pt3 = new Point2D.Double(rect.getMinX(), rect.getMaxY());

        Point2D pt4 = new Point2D.Double(rect.getMaxX(), rect.getMinY());

        convert(pt1, pt1);

        convert(pt2, pt2);

        convert(pt3, pt3);

        convert(pt4, pt4);

        double min_x = Math.min(

            Math.min(pt1.getX(), pt2.getX()),

            Math.min(pt3.getX(), pt4.getX()));

        double min_y = Math.min(

            Math.min(pt1.getY(), pt2.getY()),

            Math.min(pt3.getY(), pt4.getY()));

        double max_x = Math.max(

            Math.max(pt1.getX(), pt2.getX()),

            Math.max(pt3.getX(), pt4.getX()));

        double max_y = Math.max(

            Math.max(pt1.getY(), pt2.getY()),

            Math.max(pt3.getY(), pt4.getY()));

        return new Rectangle2D.Double(

            min_x, min_y, max_x - min_x, max_y - min_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.impl.cts.gt2;

import java.util.Locale;

import org.cresques.cts.IDatum;

import org.geotools.cs.CoordinateSystemFactory;

import org.geotools.cs.GeographicCoordinateSystem;

import org.geotools.cs.HorizontalDatum;

import org.opengis.referencing.FactoryException;

/**

 * Datum (y Ellipsoid) de GeoTools2.

 *

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

 */

public class CSDatum implements IDatum {

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

    /*

     * INSERT INTO epsg_coordinatereferencesystem VALUES (

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

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

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

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

     *

     *  DX (m) = -270.933

            DY (m) =  115.599

            DZ (m) = -360.226

            EX (") = -5.266

            EY (") = -1.238

              EZ (")  =  2.381

            FE (ppm) = -5.109

     */

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

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

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

    /*

    # La Canoa

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

    no_defs <>

    # PSAD56

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

     */

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

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

                                  "TOWGS84[-270.933,115.599,-360.226,-5.266,-1.238,2.381,-5.109],AUTHORITY[\"EPSG\",\"6288\"]]";

    private static String line9 = "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\",\"8903\"]]"+

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

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

        "AXIS[\"Lat\",NORTH],AXIS[\"Long\",EAST]," +

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

    private static String line10 = "GEOGCS[\"RGF93\"," +

		"DATUM[\"Reseau Geodesique Francais 1993\"," +

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

	    "TOWGS84[0,0,0,0,0,0,0],AUTHORITY[\"EPSG\",\"6171\"]],"+

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

        "UNIT[\"DMSH\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9108\"]]," +

        "AXIS[\"Lat\",NORTH],AXIS[\"Long\",EAST]," +

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

    private static String line11 = "GEOGCS[\"ETRS89\","+

        "DATUM[\"European_Terrestrial_Reference_System_1989\","+

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

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

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

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

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

        "UNIT[\"degree\",0.01745329251994328,"+

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

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

    //," +

    //"TOWGS84[0,0,0,0,0,0,0]

    private static String line12 =

    	"GEOGCS[\"Mars 2000\","+

	    	"DATUM[\"D_Mars_2000\","+

	    		"SPHEROID[\"Mars_2000_IAU_IAG\",3396190.0, 169.89444722361179],"+

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

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

	    	"UNIT[\"Decimal_Degree\",0.0174532925199433]]";

    public final static CSDatum wgs84 = new CSDatum(line1);

    public final static CSDatum ed50 = new CSDatum(line2);

    public final static CSDatum ntf = new CSDatum(line3);

    public final static CSDatum d73 = new CSDatum(line4);

    public final static CSDatum nad27 = new CSDatum(line5);

    public final static CSDatum nad83 = new CSDatum(line6);

    public final static CSDatum lomaQuintana = new CSDatum(line7);

    public final static CSDatum laCanoa = new CSDatum(line8);

    public static CSDatum etrs89 = null;

    public static CSDatum ntfParis = null;

    public static CSDatum posgar = null;

    public static CSDatum rgf93 = null;

    public static CSDatum mars = null;

    public static CSDatum moon = null;

    private String sGeo1 = "GEOGCS[\"WGS 84\",";

    private String sGeo2 =

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

        "UNIT[\"DMSH\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9108\"]]," +

        "AXIS[\"Lat\",NORTH],AXIS[\"Long\",EAST]," +

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

    private HorizontalDatum datum = null;

    static {

    	try {

			ntfParis =  new CSDatum().fromWKT(line9);

			rgf93 =  new CSDatum().fromWKT(line10);

			etrs89 = new CSDatum().fromWKT(line11);

			posgar = new CSDatum().fromWKT(

				"GEOGCS[\"POSGAR\","+

				"DATUM[\"POSGAR\"," +

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

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

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

                "UNIT[\"DMSH\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9108\"]]," +

                "AXIS[\"Lat\",NORTH],AXIS[\"Long\",EAST]," +

                "AUTHORITY[\"EPSG\",\"4172\"]]");

			mars = new CSDatum().fromWKT(line12);

			moon = new CSDatum().fromWKT(

				"GEOGCS[\"Moon 2000\"," +

				"DATUM[\"D_Moon_2000\"," +

				"SPHEROID[\"Moon_2000_IAU_IAG\",1737400.0, 0.0]," +

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

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

				"UNIT[\"Decimal_Degree\",0.0174532925199433]]");

		} catch (FactoryException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

    }

    public CSDatum() {

    }

    public CSDatum(HorizontalDatum datum) {

    	this.datum = datum;

    }

    public CSDatum(String sDatum) {

        try {

            fromWKT(sGeo1 + sDatum + sGeo2);

        } catch (FactoryException e) {

            // TODO Bloque catch generado autom?ticamente

            e.printStackTrace();

        }

    }

    public CSDatum fromWKT(String s) throws FactoryException {

        datum = ((GeographicCoordinateSystem) CoordinateSystemFactory.getDefault()

                .createFromWKT(s)).getHorizontalDatum();

        return this;

    }

    public String getName(Locale loc) {

        return datum.getName().toString();

    }

    HorizontalDatum getDatum() {

        return datum;

    }

    public double getESemiMajorAxis() {

        return datum.getEllipsoid().getSemiMajorAxis();

    }

    public double getEIFlattening() {

        return datum.getEllipsoid().getInverseFlattening();

    }

    public String toString() {

        return datum.toString();

    }

}