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 java.awt.geom.Point2D;

import java.awt.geom.Rectangle2D;

import org.cresques.cts.ICoordTrans;

import org.cresques.cts.IProjection;

//Geotools dependencies

import org.geotools.pt.*;

import org.geotools.ct.*;

//OpenGIS dependencies

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>

 */

public class CoordTrans implements ICoordTrans {

        private CoordinateTransformationFactory trFactory = CoordinateTransformationFactory.getDefault();

        private CoordinateTransformation tr = null;

        private MathTransform mt = null, mt2 = null, mt3 = null, mtDatum = null;

        private CoordSys from = null, to = 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() {

                ICoordTrans cti = null;

                return cti;

        }

        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;

        }

}