| gvSIG desktop 1

/* gvSIG. Sistema de Informaci?n Geogr?fica de la Generalitat Valenciana

 *

 * Copyright (C) 2004 IVER T.I. and Generalitat Valenciana.

 *

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

 *

 *  Generalitat Valenciana

 *   Conselleria d'Infraestructures i Transport

 *   Av. Blasco Ib??ez, 50

 *   46010 VALENCIA

 *   SPAIN

 *

 *      +34 963862235

 *   gvsig@gva.es

 *      www.gvsig.gva.es

 *

 *    or

 *

 *   IVER T.I. S.A

 *   Salamanca 50

 *   46005 Valencia

 *   Spain

 *

 *   +34 963163400

 *   dac@iver.es

 */

package com.iver.cit.gvsig.fmap.tools;

import java.awt.Cursor;

import java.awt.Image;

import java.awt.Point;

import java.awt.Toolkit;

import java.awt.geom.Point2D;

import javax.swing.ImageIcon;

import org.cresques.cts.IProjection;

import com.iver.cit.gvsig.fmap.MapContext;

import com.iver.cit.gvsig.fmap.MapControl;

import com.iver.cit.gvsig.fmap.ViewPort;

import com.iver.cit.gvsig.fmap.crs.CRSFactory;

import com.iver.cit.gvsig.fmap.tools.Events.MeasureEvent;

import com.iver.cit.gvsig.fmap.tools.Listeners.PolylineListener;

import com.iver.cit.gvsig.fmap.tools.geo.Geo;

/**

 * Implementaci?n de la interfaz MeasureListener como herramienta para medir el ?rea.

 * @author   Vicente Caballero Navarro

 */

public class AreaListenerImpl implements PolylineListener {

        private final Image iarea = new ImageIcon(MapControl.class.getResource(

                                "images/AreaCursor.gif")).getImage();

        private Cursor cur = Toolkit.getDefaultToolkit().createCustomCursor(iarea,

                        new Point(16, 16), "");

        protected MapControl mapCtrl;

        protected MeasureEvent event;

        /**

         * Crea un nuevo AreaListenerImpl.

         *

         * @param mc MapControl.

         */

        public AreaListenerImpl(MapControl mc) {

                this.mapCtrl = mc;

        }

        /**

         * @see com.iver.cit.gvsig.fmap.tools.Listeners.PolylineListener#points(com.iver.cit.gvsig.fmap.tools.Events.MeasureEvent)

         */

        public void points(MeasureEvent event) {

                this.event = event;

                double dist = 0;

                double distAll = 0;

                ViewPort vp = mapCtrl.getMapContext().getViewPort();

                for (int i = 0; i < (event.getXs().length - 1); i++) {

                        dist = 0;

                        Point p = new Point(event.getXs()[i].intValue(),

                                        event.getXs()[i].intValue());

                        Point p2 = new Point(event.getXs()[i + 1].intValue(),

                                        event.getXs()[i + 1].intValue());

                        ///dist = vp.toMapDistance((int) p.distance(p2));

                        dist = vp.distanceWorld(p, p2);

                        distAll += dist;

                }

                System.out.println("Per?metro = " + distAll + " ?rea = " +

                        (returnArea(vp.toMapPoint(

                                        new Point2D.Double(

                                                event.getXs()[event.getXs().length - 2].doubleValue(),

                                                event.getYs()[event.getYs().length - 2].doubleValue())))));

        }

        protected double returnArea(Point2D point) {

                Double[] xs=event.getXs();

                Double[] ys=event.getYs();

                if (mapCtrl.getProjection().isProjected()) {

                        return returnCoordsArea(xs,ys,point);

                }

                return returnGeoCArea(xs,ys,point);

        }

        /**

         * Calcula el ?rea.

         *

         * @param aux ?ltimo punto.

         *

         * @return ?rea.

         */

        public double returnCoordsArea(Double[] xs,Double[] ys, Point2D point) {

                Point2D aux=point;

                double elArea = 0.0;

                Point2D pPixel;

                Point2D p = new Point2D.Double();

                Point2D.Double pAnt = new Point2D.Double();

                ViewPort vp = mapCtrl.getMapContext().getViewPort();

                for (int pos = 0; pos < xs.length-1; pos++) {

                        pPixel = new Point2D.Double(xs[pos].doubleValue(),

                                        ys[pos].doubleValue());

                        p = pPixel;//vp.toMapPoint(pPixel);

                        if (pos == 0) {

                                pAnt.x = aux.getX();

                                pAnt.y = aux.getY();

                        }

                        elArea = elArea + ((pAnt.x - p.getX()) * (pAnt.y + p.getY()));

                        pAnt.setLocation(p);

                }

                elArea = elArea + ((pAnt.x - aux.getX()) * (pAnt.y + aux.getY()));

                elArea = Math.abs(elArea / 2.0);

                return (elArea*(Math.pow(MapContext.CHANGEM[vp.getMapUnits()],2)));

        }

public static void main(String[] args) {

        IProjection projectionUTM = CRSFactory.getCRS("EPSG:23030");

        ViewPort vpUTM = new ViewPort(projectionUTM);

        MapControl mcUTM=new MapControl();

        mcUTM.setMapContext(new MapContext(vpUTM));

        AreaListenerImpl areaListenerUTM=new AreaListenerImpl(mcUTM);

        IProjection projectionGeo = CRSFactory.getCRS("EPSG:4230");

        ViewPort vpGeo = new ViewPort(projectionGeo);

        MapControl mcGeo=new MapControl();

        mcGeo.setMapContext(new MapContext(vpGeo));

        AreaListenerImpl areaListenerGeo=new AreaListenerImpl(mcGeo);

        Double[] xsUTMCaseta=new Double[] {new Double(547508.77),new Double(547517.73),new Double(547512.65)};

        Double[] ysUTMCaseta=new Double[] {new Double(4704333.97),new Double(4704331.3),new Double(4704315.2)};

        double areaUTMCaseta=areaListenerUTM.returnCoordsArea(xsUTMCaseta,ysUTMCaseta,new Point2D.Double(547512.65,4704315.2));

        Double[] xsGeoCaseta=new Double[] {new Double(-2.42192383),new Double(-2.42181545),new Double(-2.42187771)};

        Double[] ysGeoCaseta=new Double[] {new Double(42.48914909),new Double(42.48912295),new Double(42.48897922)};

        double areaGeoCCaseta=areaListenerGeo.returnGeoCArea(xsGeoCaseta,ysGeoCaseta,new Point2D.Double(-2.42187771,42.48897922));

        System.out.println("AreaUTMCaseta = "+ areaUTMCaseta);

        System.out.println("AreaGeoCCaseta = "+ areaGeoCCaseta);

        Double[] xsUTM=new Double[] {new Double(731292),new Double(731901),new Double(730138)};

        Double[] ysUTM=new Double[] {new Double(4351223),new Double(4350768),new Double(4349232)};

        double areaUTM=areaListenerUTM.returnCoordsArea(xsUTM,ysUTM,new Point2D.Double(730138,4349232));

        Double[] xsGeo=new Double[] {new Double(-0.31888183),new Double(-0.31173131),new Double(-0.33268401)};

        Double[] ysGeo=new Double[] {new Double(39.27871741),new Double(39.27464327),new Double(39.26117368)};

        double areaGeoC=areaListenerGeo.returnGeoCArea(xsGeo,ysGeo,new Point2D.Double(-0.33268401,39.26117368));

        System.out.println("AreaUTM = "+ areaUTM);

        System.out.println("AreaGeoC = "+ areaGeoC);

        Double[] xsUTMspain=new Double[] {new Double(-12806),new Double(1025790),new Double(-31353.14)};

        Double[] ysUTMspain=new Double[] {new Double(4793276.43),new Double(4719090.94),new Double(4125607.02)};

        double areaUTMspain=areaListenerUTM.returnCoordsArea(xsUTMspain,ysUTMspain,new Point2D.Double(730138,4349232));

        Double[] xsGeospain=new Double[] {new Double(-9.22743872),new Double(3.33087936),new Double(-9.01458587),new Double(-9.22743872)};

        Double[] ysGeospain=new Double[] {new Double(43.02384666),new Double(42.38528811),new Double(37.06396689),new Double(43.02384666)};

        double areaGeospainC=areaListenerGeo.returnGeoCArea(xsGeospain,ysGeospain,new Point2D.Double(-9.01458587,37.06396689));

        System.out.println("AreaUTMSpain = "+ areaUTMspain);

        System.out.println("AreaGeoSpainC = "+ areaGeospainC);

}

        public double returnGeoCArea(Double[] xs,Double[] ys,Point2D point) {

                double[] lat=new double[xs.length];

                double[] lon=new double[xs.length];

                for (int K= 0; K < xs.length; K++){

                        lon[K]= xs[K].doubleValue()/Geo.Degree;

                        lat[K]= ys[K].doubleValue()/Geo.Degree;

                }

                return (Geo.sphericalPolyArea(lat,lon,xs.length-1)*Geo.SqM);///1.29132441;//Esto es una constante para ajustar el resultado, que por alg?n motivo no es correcto.

        }

        /**

         * @see com.iver.cit.gvsig.fmap.tools.Listeners.ToolListener#getCursor()

         */

        public Cursor getCursor() {

                return cur;

        }

        /**

         * @see com.iver.cit.gvsig.fmap.tools.Listeners.PolylineListener#pointFixed(com.iver.cit.gvsig.fmap.tools.Events.MeasureEvent)

         */

        public void pointFixed(MeasureEvent event) {

        }

        /**

         * @see com.iver.cit.gvsig.fmap.tools.Listeners.PolylineListener#polylineFinished(com.iver.cit.gvsig.fmap.tools.Events.MeasureEvent)

         */

        public void polylineFinished(MeasureEvent event) {

        }

        /**

         * @see com.iver.cit.gvsig.fmap.tools.Listeners.ToolListener#cancelDrawing()

         */

        public boolean cancelDrawing() {

                return false;

        }

}