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

import org.geotools.cs.AxisInfo;

import org.geotools.cs.GeographicCoordinateSystem;

import org.geotools.cs.Projection;

import org.geotools.units.Unit;

import org.opengis.referencing.FactoryException;

/**

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

 *

 */

public class CSLambertCC extends CoordSys {

    public CSLambertCC(CSDatum datum, double meridian, double origin,

                       double sp1, double sp2, double est, double nort) {

        super(datum);

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

                                                datum.getDatum());

        Unit linearUnit = Unit.METRE;

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

        params.setParameter("semi_major",

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

        params.setParameter("semi_minor",

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

        params.setParameter("central_meridian", meridian);

        params.setParameter("latitude_of_origin", origin);

        params.setParameter("standard_parallel_1", sp1);

        params.setParameter("standard_parallel_2", sp2);

        params.setParameter("false_easting", est);

        params.setParameter("false_northing", nort);

        try {

            Projection projection = csFactory.createProjection("Lambert",

                                                               "Lambert_Conformal_Conic_2SP",

                                                               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 CSLambertCC(CSDatum datum, double meridian, double origin,

                    double sf, double est, double nort) {

            super(datum);

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

                            datum.getDatum());

            Unit linearUnit = Unit.METRE;

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

            params.setParameter("semi_major",

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

            params.setParameter("semi_minor",

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

            params.setParameter("central_meridian", meridian);

            params.setParameter("latitude_of_origin", origin);

            params.setParameter("scale_factor", sf);

            params.setParameter("false_easting", est);

            params.setParameter("false_northing", nort);

            try {

                    Projection projection = csFactory.createProjection("Lambert",

                                    "Lambert_Conformal_Conic_1SP",

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

    }

}