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

    }

}