Application: gvSIG desktop

package org.cresques.geo;

public class Ellipsoid {

        //'WGS72':        (6378135.0, 298.26, 1972),

        public static Ellipsoid wgs72 = new Ellipsoid("WGS72", 6378135.0, 298.26);

        //'WGS84':        (6378137.0, 298.257223563, 1984),

        public static Ellipsoid wgs84 = new Ellipsoid("WGS84", 6378137.0, 298.257223563);

        //'GRS80':        (6378137.0, 298.257222101, 1980), // , 6356752.314

        public static Ellipsoid grs80 = new Ellipsoid("GRS80", 6378137.0, 298.257222101); 

        //'Hayford':        (6378388.0, 297.0, 1909),

        public static Ellipsoid hayford = new Ellipsoid("Hayford", 6378388.0, 297.0);

        //'Struve':        (6378298.0, 299.73, 1860),

        public static Ellipsoid struve = new Ellipsoid("Struve", 6378298.0, 299.73);

        // ed50 ... ?hayford?

        public static Ellipsoid ed50 = new Ellipsoid("ED50", 6378388.0, 297.0);

        // Clarcke 1866        a = 6378206.4;        f = 294.9786982;

        public static Ellipsoid clarke66 = new Ellipsoid("CLARKE 1866", 6378206.4, 294.9786982);

        private String pName = null;

        private double pMajor = 0.0;

        private double pFlat = 0.0;

        public double a, b, f, e, e2, ep, ep2;

        public Ellipsoid(String name, double major, double flat) {

                pName = name;

                pMajor = major;

                pFlat = flat;

                a = pMajor;

                f = pFlat;

                ainvto();

        }

        private void ainvto()

        {

                b = (a * (pFlat - 1.0D)) / pFlat;

                e2 = (2D * pFlat - 1.0D) / (pFlat * pFlat);

                e = Math.sqrt(e2);

                ep2 = e2 / (1.0D - e2);

                ep = Math.sqrt(ep2);

        }

        private void abto()

        {

                if(Math.abs(a - b) < 9.9999999999999995E-021D)

                {

                        pFlat = 0.0D;

                        e2 = 0.0D;

                } else

                {

                        pFlat = a / (a - b);

                        e2 = (2D * pFlat - 1.0D) / (pFlat * pFlat);

                }

                e = Math.sqrt(e2);

                ep2 = e2 / (1.0D - e2);

                ep = Math.sqrt(ep2);

        }

        public double [] getParam() {

                double [] elipar = new double[8];

                elipar[1] = pMajor;

                elipar[2] = pFlat;

                elipar[3] = b;

                elipar[4] = e;

                elipar[5] = e2;

                elipar[6] = ep;

                elipar[7] = ep2;

                return elipar;

        }

}