Application: gvSIG desktop

/*

 * Cresques Mapping Suite. Graphic Library for constructing mapping applications.

 *

 * Copyright (C) 2004-5.

 *

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

 * 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., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 *

 * For more information, contact:

 *

 * cresques@gmail.com

 */

package org.cresques.cts.gt2;

import org.cresques.cts.ICoordTrans;

import org.cresques.cts.IProjection;

import org.geotools.ct.*;

//Geotools dependencies

import org.geotools.pt.*;

//OpenGIS dependencies

import org.opengis.referencing.operation.TransformException;

import java.awt.geom.Point2D;

import java.awt.geom.Rectangle2D;

//import org.geotools.pt.MismatchedDimensionException;

/**

 * Transforma coordenadas entre dos sistemas

 * @see org.creques.cts.CoordSys

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

 */

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

        pt1 = convert(pt1, null);

        pt2 = convert(pt2, null);

        rect = new Rectangle2D.Double();

        rect.setFrameFromDiagonal(pt1, pt2);

        return rect;

    }

}