Statistics
| Revision:

svn-gvsig-desktop / trunk / org.gvsig.desktop / org.gvsig.desktop.compat.cdc / org.gvsig.fmap.geometry / org.gvsig.fmap.geometry.operation / src / main / java / org / gvsig / fmap / geom / operation / utils / PointGetAngle.java @ 40559

History | View | Annotate | Download (3.24 KB)

1 40559 jjdelcerro
/**
2
 * gvSIG. Desktop Geographic Information System.
3 40435 jjdelcerro
 *
4 40559 jjdelcerro
 * Copyright (C) 2007-2013 gvSIG Association.
5
 *
6 40435 jjdelcerro
 * This program is free software; you can redistribute it and/or
7
 * modify it under the terms of the GNU General Public License
8 40559 jjdelcerro
 * as published by the Free Software Foundation; either version 3
9 40435 jjdelcerro
 * of the License, or (at your option) any later version.
10 40559 jjdelcerro
 *
11 40435 jjdelcerro
 * This program is distributed in the hope that it will be useful,
12
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14
 * GNU General Public License for more details.
15 40559 jjdelcerro
 *
16 40435 jjdelcerro
 * You should have received a copy of the GNU General Public License
17
 * along with this program; if not, write to the Free Software
18 40559 jjdelcerro
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
19 40435 jjdelcerro
 * MA  02110-1301, USA.
20 40559 jjdelcerro
 *
21
 * For any additional information, do not hesitate to contact us
22
 * at info AT gvsig.com, or visit our website www.gvsig.com.
23 40435 jjdelcerro
 */
24
/*
25
 * AUTHORS (In addition to CIT):
26
 * 2009 {Iver T.I.}   {Task}
27
 */
28
29
package org.gvsig.fmap.geom.operation.utils;
30
31
import org.gvsig.fmap.geom.Geometry;
32
import org.gvsig.fmap.geom.GeometryLocator;
33
import org.gvsig.fmap.geom.GeometryManager;
34
import org.gvsig.fmap.geom.Geometry.TYPES;
35
import org.gvsig.fmap.geom.operation.GeometryOperation;
36
import org.gvsig.fmap.geom.operation.GeometryOperationContext;
37
import org.gvsig.fmap.geom.operation.GeometryOperationException;
38
import org.gvsig.fmap.geom.operation.GeometryOperationNotSupportedException;
39
import org.gvsig.fmap.geom.operation.distance.PointDistance;
40
import org.gvsig.fmap.geom.primitive.Point;
41
import org.gvsig.fmap.geom.type.GeometryType;
42
import org.gvsig.tools.locator.LocatorException;
43
44
/**
45
 * @author <a href="mailto:jpiera@gvsig.org">Jorge Piera</a>
46
 */
47
public class PointGetAngle extends GeometryOperation{
48
    public static final String NAME = "getAngle";
49
        public static final int CODE = GeometryLocator.getGeometryManager().
50
            getGeometryOperationCode(NAME);
51
52
        /* (non-Javadoc)
53
         * @see org.gvsig.fmap.geom.operation.GeometryOperation#getOperationIndex()
54
         */
55
        public int getOperationIndex() {
56
                return CODE;
57
        }
58
59
        /* (non-Javadoc)
60
         * @see org.gvsig.fmap.geom.operation.GeometryOperation#invoke(org.gvsig.fmap.geom.Geometry, org.gvsig.fmap.geom.operation.GeometryOperationContext)
61
         */
62
        public Object invoke(Geometry geom, GeometryOperationContext ctx)
63
        throws GeometryOperationException {
64
65
                Geometry geom2 = (Geometry)ctx.getAttribute("geom");
66
                GeometryType geomType2 = geom2.getGeometryType();
67
68
                if ((TYPES.POINT != geom.getType()) && TYPES.POINT != geomType2.getType()){
69
                        throw new UnsupportedOperationException("The distance only can be execudet between two points");
70
                }
71
72
                Point start = (Point)geom;
73
                Point end = (Point)geom;
74
75
                double distance;
76
                try {
77
                        distance = ((Double)GeometryLocator.getGeometryManager().invokeOperation(PointDistance.CODE, geom, ctx)).doubleValue();
78
                } catch (LocatorException e) {
79
                        throw new GeometryOperationException(e);
80
                } catch (GeometryOperationNotSupportedException e) {
81
                        throw new GeometryOperationException(e);
82
                }
83
84
                double angle = Math.acos((end.getX() - start.getX()) / distance);
85
86
                if (start.getY() > end.getY()) {
87
                        angle = -angle;
88
                }
89
90
                if (angle < 0) {
91
                        angle += (2 * Math.PI);
92
                }
93
94
                return new Double(angle);
95
        }
96
97
}