Statistics
| Revision:

svn-gvsig-desktop / trunk / org.gvsig.desktop / org.gvsig.desktop.compat.cdc / org.gvsig.fmap.geometry / org.gvsig.fmap.geometry.jts / src / main / java / org / gvsig / fmap / geom / jts / primitive / curve / arc / AbstractArc.java @ 42464

History | View | Annotate | Download (16.1 KB)

1
/* gvSIG. Desktop Geographic Information System.
2
 *
3
 * Copyright ? 2007-2015 gvSIG Association
4
 *
5
 * This program is free software; you can redistribute it and/or
6
 * modify it under the terms of the GNU General Public License
7
 * as published by the Free Software Foundation; either version 2
8
 * of the License, or (at your option) any later version.
9
 *
10
 * This program is distributed in the hope that it will be useful,
11
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
12
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
13
 * GNU General Public License for more details.
14
 *
15
 * You should have received a copy of the GNU General Public License
16
 * along with this program; if not, write to the Free Software
17
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
18
 * MA  02110-1301, USA.
19
 *
20
 * For any additional information, do not hesitate to contact us
21
 * at info AT gvsig.com, or visit our website www.gvsig.com.
22
 */
23
package org.gvsig.fmap.geom.jts.primitive.curve.arc;
24

    
25
import java.awt.Shape;
26
import java.awt.geom.AffineTransform;
27
import java.awt.geom.PathIterator;
28

    
29
import org.cresques.cts.CoordTransRuntimeException;
30
import org.cresques.cts.ICoordTrans;
31

    
32
import org.gvsig.fmap.geom.Geometry;
33
import org.gvsig.fmap.geom.jts.gputils.DefaultGeneralPathX;
34
import org.gvsig.fmap.geom.jts.primitive.curve.AbstractCurve;
35
import org.gvsig.fmap.geom.jts.primitive.point.PointJTS;
36
import org.gvsig.fmap.geom.jts.util.UtilFunctions;
37
import org.gvsig.fmap.geom.operation.GeometryOperationException;
38
import org.gvsig.fmap.geom.operation.GeometryOperationNotSupportedException;
39
import org.gvsig.fmap.geom.primitive.Arc;
40
import org.gvsig.fmap.geom.primitive.GeneralPathX;
41
import org.gvsig.fmap.geom.primitive.Point;
42

    
43

    
44
/**
45
 * @author fdiaz
46
 *
47
 */
48
public abstract class AbstractArc extends AbstractCurve implements Arc {
49

    
50
    /**
51
     *
52
     */
53
    private static final long serialVersionUID = 454301669807892457L;
54

    
55
    /**
56
     * @param type
57
     * @param subtype
58
     */
59
    protected AbstractArc(int subtype) {
60
        super(Geometry.TYPES.ARC, subtype);
61
    }
62

    
63
    protected Point init;
64

    
65
    /**
66
     * This is the middle point (belongs to the arc), not the center
67
     * of the circle/ellipse
68
     */
69
    protected Point middle;
70
    protected Point end;
71

    
72

    
73
    /* (non-Javadoc)
74
     * @see org.gvsig.fmap.geom.primitive.Curve#setPoints(org.gvsig.fmap.geom.primitive.Point, org.gvsig.fmap.geom.primitive.Point)
75
     */
76
    public void setPoints(Point initialPoint, Point endPoint) {
77
        String message = "Calling deprecated method setPoints of a arc";
78
        notifyDeprecated(message);
79
        throw new UnsupportedOperationException(message);
80
    }
81

    
82
    /* (non-Javadoc)
83
     * @see org.gvsig.fmap.geom.primitive.OrientablePrimitive#getCoordinateAt(int, int)
84
     */
85
    public double getCoordinateAt(int index, int dimension) {
86
        String message = "Calling deprecated method setPoints of a arc";
87
        notifyDeprecated(message);
88
        throw new UnsupportedOperationException(message);
89
    }
90

    
91
    /* (non-Javadoc)
92
     * @see org.gvsig.fmap.geom.primitive.OrientablePrimitive#setCoordinateAt(int, int, double)
93
     */
94
    public void setCoordinateAt(int index, int dimension, double value) {
95
        String message = "Calling deprecated method setPoints of a arc";
96
        notifyDeprecated(message);
97
        throw new UnsupportedOperationException(message);
98
    }
99

    
100
    /* (non-Javadoc)
101
     * @see org.gvsig.fmap.geom.primitive.OrientablePrimitive#addVertex(org.gvsig.fmap.geom.primitive.Point)
102
     */
103
    public void addVertex(Point point) {
104
        String message = "Calling deprecated method setPoints of a arc";
105
        notifyDeprecated(message);
106
        throw new UnsupportedOperationException(message);
107
    }
108

    
109
    /* (non-Javadoc)
110
     * @see org.gvsig.fmap.geom.primitive.OrientablePrimitive#addVertex(double, double)
111
     */
112
    public void addVertex(double x, double y) {
113
        String message = "Calling deprecated method setPoints of a arc";
114
        notifyDeprecated(message);
115
        throw new UnsupportedOperationException(message);
116
    }
117

    
118
    /* (non-Javadoc)
119
     * @see org.gvsig.fmap.geom.primitive.OrientablePrimitive#addVertex(double, double, double)
120
     */
121
    public void addVertex(double x, double y, double z) {
122
        String message = "Calling deprecated method setPoints of a arc";
123
        notifyDeprecated(message);
124
        throw new UnsupportedOperationException(message);
125
    }
126

    
127
    /* (non-Javadoc)
128
     * @see org.gvsig.fmap.geom.primitive.OrientablePrimitive#removeVertex(int)
129
     */
130
    public void removeVertex(int index) {
131
        String message = "Calling deprecated method setPoints of a arc";
132
        notifyDeprecated(message);
133
        throw new UnsupportedOperationException(message);
134
    }
135

    
136
    /* (non-Javadoc)
137
     * @see org.gvsig.fmap.geom.primitive.OrientablePrimitive#getVertex(int)
138
     */
139
    public Point getVertex(int index) {
140
        String message = "Calling deprecated method setPoints of a arc";
141
        notifyDeprecated(message);
142
        throw new UnsupportedOperationException(message);
143
    }
144

    
145
    /* (non-Javadoc)
146
     * @see org.gvsig.fmap.geom.primitive.OrientablePrimitive#getNumVertices()
147
     */
148
    public int getNumVertices() {
149
        String message = "Calling deprecated method setPoints of a arc";
150
        notifyDeprecated(message);
151
        throw new UnsupportedOperationException(message);
152
    }
153

    
154
    /* (non-Javadoc)
155
     * @see org.gvsig.fmap.geom.primitive.OrientablePrimitive#insertVertex(int, org.gvsig.fmap.geom.primitive.Point)
156
     */
157
    public void insertVertex(int index, Point p) {
158
        String message = "Calling deprecated method setPoints of a arc";
159
        notifyDeprecated(message);
160
        throw new UnsupportedOperationException(message);
161
    }
162

    
163
    /* (non-Javadoc)
164
     * @see org.gvsig.fmap.geom.primitive.OrientablePrimitive#setVertex(int, org.gvsig.fmap.geom.primitive.Point)
165
     */
166
    public void setVertex(int index, Point p) {
167
        String message = "Calling deprecated method setPoints of a arc";
168
        notifyDeprecated(message);
169
        throw new UnsupportedOperationException(message);
170
    }
171

    
172
    /* (non-Javadoc)
173
     * @see org.gvsig.fmap.geom.primitive.OrientablePrimitive#setGeneralPath(org.gvsig.fmap.geom.primitive.GeneralPathX)
174
     */
175
    public void setGeneralPath(GeneralPathX generalPathX) {
176
        String message = "Calling deprecated method setPoints of a arc";
177
        notifyDeprecated(message);
178
        throw new UnsupportedOperationException(message);
179
    }
180

    
181
    /* (non-Javadoc)
182
     * @see org.gvsig.fmap.geom.primitive.OrientablePrimitive#addMoveToVertex(org.gvsig.fmap.geom.primitive.Point)
183
     */
184
    public void addMoveToVertex(Point point) {
185
        String message = "Calling deprecated method setPoints of a arc";
186
        notifyDeprecated(message);
187
        throw new UnsupportedOperationException(message);
188
    }
189

    
190
    /* (non-Javadoc)
191
     * @see org.gvsig.fmap.geom.primitive.OrientablePrimitive#closePrimitive()
192
     */
193
    public void closePrimitive() {
194
        String message = "Calling deprecated method setPoints of a arc";
195
        notifyDeprecated(message);
196
        throw new UnsupportedOperationException(message);
197
    }
198

    
199
    /* (non-Javadoc)
200
     * @see org.gvsig.fmap.geom.primitive.OrientablePrimitive#ensureCapacity(int)
201
     */
202
    public void ensureCapacity(int capacity) {
203
        // TODO Auto-generated method stub
204

    
205
    }
206

    
207
    /* (non-Javadoc)
208
     * @see org.gvsig.fmap.geom.Geometry#reProject(org.cresques.cts.ICoordTrans)
209
     */
210
    public void reProject(ICoordTrans ct) {
211
        //FIXME: Esto solo ser?a correcto para transformaciones de traslaci?n, rotaci?n y escala
212
        // Ser?a incorrecto para las de deformaci?n en cizallamiento
213

    
214
        try {
215
        init.reProject(ct);
216
        middle.reProject(ct);
217
        end.reProject(ct);
218
        } catch (CoordTransRuntimeException e){
219
            //Si ha fallado la reproyecci?n de alguno de los puntos, ponemos todas las coordenadas a 0
220
            init.setX(0);
221
            init.setY(0);
222
            middle.setX(0);
223
            middle.setY(0);
224
            end.setX(0);
225
            end.setY(0);
226
        }
227
    }
228

    
229
    /* (non-Javadoc)
230
     * @see org.gvsig.fmap.geom.Geometry#transform(java.awt.geom.AffineTransform)
231
     */
232
    public void transform(AffineTransform at) {
233
        //FIXME: Esto solo ser?a correcto para transformaciones de traslaci?n, rotaci?n y escala
234
        // Ser?a incorrecto para las de deformaci?n en cizallamiento
235

    
236
        init.transform(at);
237
        middle.transform(at);
238
        end.transform(at);
239
    }
240

    
241
    /* (non-Javadoc)
242
     * @see org.gvsig.fmap.geom.Geometry#getDimension()
243
     */
244
    public int getDimension() {
245
        return init.getDimension();
246
    }
247

    
248
    /* (non-Javadoc)
249
     * @see org.gvsig.fmap.geom.Geometry#getShape(java.awt.geom.AffineTransform)
250
     */
251
    public Shape getShape(AffineTransform affineTransform) {
252
        return new DefaultGeneralPathX(getPathIterator(affineTransform),false,0);
253
    }
254

    
255
    /* (non-Javadoc)
256
     * @see org.gvsig.fmap.geom.Geometry#getShape()
257
     */
258
    public Shape getShape() {
259
        return getShape(null);
260
    }
261

    
262
    /* (non-Javadoc)
263
     * @see org.gvsig.fmap.geom.jts.GeometryJTS#is3D()
264
     */
265
    public boolean is3D() {
266
        return ((PointJTS)init).is3D();
267
    }
268

    
269
    /**
270
     * @param initialPoint
271
     * @return
272
     */
273
    protected abstract Point fixPoint(Point point);
274

    
275
    /* (non-Javadoc)
276
     * @see org.gvsig.fmap.geom.primitive.Arc#setPoints(org.gvsig.fmap.geom.primitive.Point, org.gvsig.fmap.geom.primitive.Point, org.gvsig.fmap.geom.primitive.Point)
277
     */
278
    public void setPoints(Point startPoint, Point midPoint, Point endPoint) {
279
        init = fixPoint(startPoint);
280
        middle = fixPoint(midPoint);
281
        end = fixPoint(endPoint);
282
    }
283

    
284
    /* (non-Javadoc)
285
     * @see org.gvsig.fmap.geom.primitive.Arc#getInitPoint()
286
     */
287
    public Point getInitPoint() {
288
        return init;
289
    }
290

    
291
    /* (non-Javadoc)
292
     * @see org.gvsig.fmap.geom.primitive.Arc#getEndPoint()
293
     */
294
    public Point getEndPoint() {
295
        return end;
296
    }
297

    
298
    /* (non-Javadoc)
299
     * @see org.gvsig.fmap.geom.primitive.Arc#getMiddlePoint()
300
     */
301
    public Point getMiddlePoint() {
302
        return middle;
303
    }
304

    
305
    /**
306
     * Leaves the angle between PI and -PI
307
     * @param angle (radians)
308
     * @return
309
     */
310
    protected double normalizeAngle(double angle) {
311
        if (angle > -Math.PI && angle <= Math.PI) {
312
            return angle;
313
        }
314

    
315
        if (angle == Double.NEGATIVE_INFINITY || angle == Double.POSITIVE_INFINITY) {
316
            return 0;
317
        }
318

    
319
        double abs_ang = Math.abs(angle);
320
        double remove = Math.floor(abs_ang / (2 * Math.PI));
321
        remove = remove * 2 * Math.PI;
322
        double resp = 0;
323

    
324
        if (angle > 0) {
325
            resp = angle - remove;
326
            if (resp > Math.PI) {
327
                // final adjustment
328
                resp = resp - 2 * Math.PI;
329
            }
330
        } else {
331
            resp = angle + remove;
332
            if (resp <= -Math.PI) {
333
                // final adjustment
334
                resp = resp + 2 * Math.PI;
335
            }
336
        }
337

    
338
        return resp;
339
    }
340

    
341

    
342
    /* (non-Javadoc)
343
     * @see org.gvsig.fmap.geom.primitive.Arc#setPointsStartExt(org.gvsig.fmap.geom.primitive.Point, double, double, double)
344
     */
345
    public void setPointsStartExt(Point center, double radius, double startAngle, double angleExt) {
346
        setPoints(center, radius, startAngle, angleExt);
347
    }
348

    
349
    /* (non-Javadoc)
350
     * @see org.gvsig.fmap.geom.primitive.Arc#setPointsStartEnd(org.gvsig.fmap.geom.primitive.Point, double, double, double)
351
     */
352
    public void setPointsStartEnd(Point center, double radius, double startAngle, double endAngle) {
353

    
354
        if (startAngle == endAngle) {
355
            setPointsStartExt(center, radius, startAngle, 0);
356
        } else {
357

    
358
            /*
359
             * Normalize then force clockwise:
360
             */
361
            double norm_start = normalizeAngle(startAngle);
362
            double norm_end = normalizeAngle(endAngle);
363
            double ang_ext = 0;
364

    
365
            // clockwise
366
            // ang_ext must be positive
367
            if (norm_start >= norm_end) {
368
                ang_ext = norm_start - norm_end;
369
            } else {
370
                ang_ext = 2 * Math.PI - (norm_end - norm_start);
371
            }
372
            setPointsStartExt(center, radius, startAngle, ang_ext);
373

    
374
            // finally call other method with ang_ext
375
        }
376
    }
377

    
378
    /* (non-Javadoc)
379
     * @see org.gvsig.fmap.geom.Geometry#getGeneralPath()
380
     */
381
    public GeneralPathX getGeneralPath() {
382

    
383
       GeneralPathX gp = new DefaultGeneralPathX(getPathIterator(null, getManager().getFlatness()), is3D(), 0.0);
384
        return gp;
385
    }
386

    
387

    
388
    /* (non-Javadoc)
389
     * @see org.gvsig.fmap.geom.Geometry#getPathIterator(java.awt.geom.AffineTransform)
390
     */
391
    public PathIterator getPathIterator(AffineTransform at) {
392
        return getPathIterator(at, getManager().getFlatness());
393
    }
394

    
395

    
396
    /* (non-Javadoc)
397
     * @see org.gvsig.fmap.geom.Geometry#getPathIterator(java.awt.geom.AffineTransform, double)
398
     */
399
    public PathIterator getPathIterator(AffineTransform at, double flatness) {
400

    
401
        java.awt.geom.Point2D.Double p1 = new java.awt.geom.Point2D.Double(init.getX(), init.getY());
402
        java.awt.geom.Point2D.Double p2 = new java.awt.geom.Point2D.Double(middle.getX(), middle.getY());
403
        java.awt.geom.Point2D.Double p3 = new java.awt.geom.Point2D.Double(end.getX(), end.getY());
404

    
405
        java.awt.geom.Arc2D arco = UtilFunctions.createArc(p1, p2, p3);
406
        if (arco == null) {
407
            logger.info("Did not set arc points (probably aligned points): " + p1.getX() + " " + p1.getY() + " :: "
408
                + p2.getX() + " " + p2.getY() + " :: " + p3.getX() + " " + p3.getY());
409
            throw new IllegalArgumentException("Did not set arc points (probably aligned points).");
410
        }
411

    
412
        return arco.getPathIterator(at, flatness);
413
    }
414

    
415
    /* (non-Javadoc)
416
     * @see org.gvsig.fmap.geom.jts.GeometryJTS#flip()
417
     */
418
    public void flip() throws GeometryOperationNotSupportedException, GeometryOperationException {
419
        Point aux = init;
420
        init = end;
421
        end = aux;
422
    }
423

    
424

    
425
    /* (non-Javadoc)
426
     * @see org.gvsig.fmap.geom.primitive.Arc#getStartAngle()
427
     */
428
    public double getStartAngle() throws GeometryOperationNotSupportedException, GeometryOperationException {
429
        return getAngle(getCenterPoint(), getInitPoint());
430
    }
431

    
432
    /* (non-Javadoc)
433
     * @see org.gvsig.fmap.geom.primitive.Arc#getEndAngle()
434
     */
435
    public double getEndAngle() throws GeometryOperationNotSupportedException, GeometryOperationException {
436
        return getAngle(getCenterPoint(), getEndPoint());
437
    }
438

    
439
    private double getAngle(Point start, Point end) throws GeometryOperationNotSupportedException, GeometryOperationException {
440
        double angle = Math.acos((end.getX() - start.getX()) / start.distance(end));
441

    
442
        if (start.getY() > end.getY()) {
443
            angle = -angle;
444
        }
445

    
446
        if (angle < 0) {
447
            angle += (2 * Math.PI);
448
        }
449

    
450
        return angle;
451
    }
452

    
453

    
454
    /* (non-Javadoc)
455
     * @see org.gvsig.fmap.geom.Geometry#offset(double)
456
     */
457
    public Geometry offset(double distance) throws GeometryOperationNotSupportedException, GeometryOperationException {
458
        // TODO Auto-generated method stub
459
        Point center = getCenterPoint();
460
        double radius = center.distance(init);
461
        double scale = (radius+distance)/radius;
462
        AffineTransform at = getScaleAffineTransform(center, scale);
463
        Geometry cloned = this.cloneGeometry();
464
        cloned.transform(at);
465
        return cloned;
466
    }
467

    
468
    protected AffineTransform getScaleAffineTransform(Point center, Double scale)
469
        throws GeometryOperationNotSupportedException,
470
        GeometryOperationException {
471

    
472
        AffineTransform translate =
473
            AffineTransform
474
                .getTranslateInstance(-center.getX(), -center.getY());
475

    
476
        AffineTransform scaleTransform = AffineTransform.getScaleInstance(scale,scale);
477

    
478
        AffineTransform inverseTranslate =
479
            AffineTransform.getTranslateInstance(center.getX(), center.getY());
480
        AffineTransform at = new AffineTransform(translate);
481

    
482
        at.preConcatenate(scaleTransform);
483
        at.preConcatenate(inverseTranslate);
484
        return at;
485
    }
486
}