Application: gvSIG desktop

/* gvSIG. Desktop Geographic Information System.

 *

 * Copyright ? 2007-2015 gvSIG Association

 *

 * 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., 51 Franklin Street, Fifth Floor, Boston,

 * MA  02110-1301, USA.

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

package org.gvsig.fmap.geom.jts.primitive.curve.arc;

import java.awt.geom.AffineTransform;

import java.awt.geom.PathIterator;

import com.vividsolutions.jts.geom.Coordinate;

import org.gvsig.fmap.geom.Geometry;

import org.gvsig.fmap.geom.GeometryException;

import org.gvsig.fmap.geom.aggregate.MultiLine;

import org.gvsig.fmap.geom.aggregate.MultiPoint;

import org.gvsig.fmap.geom.aggregate.MultiPolygon;

import org.gvsig.fmap.geom.jts.aggregate.MultiLine3D;

import org.gvsig.fmap.geom.jts.aggregate.MultiPoint3D;

import org.gvsig.fmap.geom.jts.aggregate.MultiPolygon3D;

import org.gvsig.fmap.geom.jts.primitive.curve.line.Line3D;

import org.gvsig.fmap.geom.jts.primitive.point.Point3D;

import org.gvsig.fmap.geom.jts.primitive.point.PointJTS;

import org.gvsig.fmap.geom.jts.primitive.surface.polygon.Polygon3D;

import org.gvsig.fmap.geom.jts.util.ArrayListCoordinateSequence;

import org.gvsig.fmap.geom.jts.util.JTSUtils;

import org.gvsig.fmap.geom.operation.GeometryOperationException;

import org.gvsig.fmap.geom.operation.GeometryOperationNotSupportedException;

import org.gvsig.fmap.geom.primitive.Line;

import org.gvsig.fmap.geom.primitive.Point;

import org.gvsig.fmap.geom.primitive.Polygon;

/**

 * @author fdiaz

 *

 */

public class Arc2DZ extends AbstractArc {

    /**

     *

     */

    private static final long serialVersionUID = -5691008010954470053L;

    private double zValue = Double.NaN;

    /**

     * @param subtype

     */

    public Arc2DZ() {

        super(Geometry.SUBTYPES.GEOM3D);

    }

    /*

     * (non-Javadoc)

     *

     * @see

     * org.gvsig.fmap.geom.primitive.Arc#setPoints(org.gvsig.fmap.geom.primitive

     * .Point, double, double, double)

     */

    public void setPoints(Point center, double radius, double startAngle, double angleExt) {

        center = fixPoint(center);

        double diameter = radius * 2;

        double x;

        double y;

        double start;

        double extent;

        double angleOffset;

        final double pi2 = Math.PI * 2;

        angleOffset = 0;

        if (angleExt < 0) {

            angleOffset = pi2 + angleExt;

            angleExt = Math.abs(angleExt);

        }

        x = center.getX() - radius;

        y = center.getY() - radius;

        if (angleExt > 0 && (angleExt % pi2) == 0) {

            start = 0;

            extent = 360;

        } else {

            angleExt = angleExt % pi2; // Asumimos que aqui angleExt es siempre

                                       // positivo

            startAngle = Math.abs(startAngle) % pi2;

            start = Math.toDegrees(pi2 - startAngle + angleOffset);

            extent = -Math.toDegrees(pi2 - angleExt);

        }

        double angleStart = Math.toRadians(-start);

        double initX = x + (Math.cos(angleStart) * 0.5 + 0.5) * diameter;

        double initY = y + (Math.sin(angleStart) * 0.5 + 0.5) * diameter;

        init = new Point3D(initX, initY, ((Point3D) center).getZ());

        double angleEnd = Math.toRadians(-start - extent);

        double endX = x + (Math.cos(angleEnd) * 0.5 + 0.5) * diameter;

        double endY = y + (Math.sin(angleEnd) * 0.5 + 0.5) * diameter;

        end = new Point3D(endX, endY, ((Point3D) center).getZ());

        double angleMiddle = Math.toRadians(-start - extent) / 2;

        double middleX = x + (Math.cos(angleMiddle) * 0.5 + 0.5) * diameter;

        double middleY = y + (Math.sin(angleMiddle) * 0.5 + 0.5) * diameter;

        middle = new Point3D(middleX, middleY, ((Point3D) center).getZ());

    }

    /*

     * (non-Javadoc)

     *

     * @see org.gvsig.fmap.geom.primitive.Arc#getCenterPoint()

     */

    public Point getCenterPoint() {

        ((PointJTS) init).getJTS();

        Point3D center = new Point3D(JTSUtils.getCircumcentre(init, middle, end));

        return center;

    }

    /*

     * (non-Javadoc)

     *

     * @see org.gvsig.fmap.geom.Geometry#cloneGeometry()

     */

    public Geometry cloneGeometry() {

        Arc2DZ arc2D = new Arc2DZ();

        arc2D.setPoints((Point)init.cloneGeometry(), (Point)middle.cloneGeometry(), (Point)end.cloneGeometry());

        return arc2D;

    }

    /*

     * (non-Javadoc)

     *

     * @see

     * org.gvsig.fmap.geom.jts.primitive.curve.line.AbstractLine#fixPoint(org

     * .gvsig.fmap.geom.primitive.Point)

     */

    @Override

    protected Point fixPoint(Point point) {

        if (point instanceof Point3D) {

            Point3D point3D = (Point3D) point;

            if (Double.isNaN(zValue) && !(Double.isNaN(point3D.getZ()))) {

                zValue = point3D.getZ();

                if (init!=null) {

                    ((Point3D) init).setZ(zValue);

                }

                if (middle!=null) {

                    ((Point3D) middle).setZ(zValue);

                }

                if (end!=null) {

                    ((Point3D) end).setZ(zValue);

                }

                return point3D;

            }

        }

        return new Point3D(point.getX(), point.getY(), Double.isNaN(zValue)?0:zValue);

    }

    /*

     * (non-Javadoc)

     *

     * @see org.gvsig.fmap.geom.jts.GeometryJTS#getJTS()

     */

    public com.vividsolutions.jts.geom.Geometry getJTS() {

        PathIterator pi = getPathIterator(null);

        ArrayListCoordinateSequence coordinates = new ArrayListCoordinateSequence();

        double coords[] = new double[6];

        while (!pi.isDone()) {

            switch (pi.currentSegment(coords)) {

            case PathIterator.SEG_MOVETO:

                coordinates.add(new Coordinate(coords[0], coords[1], zValue));

                break;

            case PathIterator.SEG_LINETO:

                coordinates.add(new Coordinate(coords[0], coords[1], zValue));

                break;

            case PathIterator.SEG_QUADTO:

                coordinates.add(new Coordinate(coords[0], coords[1], zValue));

                coordinates.add(new Coordinate(coords[2], coords[3], zValue));

                break;

            case PathIterator.SEG_CUBICTO:

                coordinates.add(new Coordinate(coords[0], coords[1], zValue));

                coordinates.add(new Coordinate(coords[2], coords[3], zValue));

                coordinates.add(new Coordinate(coords[4], coords[5], zValue));

                break;

            case PathIterator.SEG_CLOSE:

                coordinates.add(coordinates.get(0));

                break;

            }

            pi.next();

        }

        return JTSUtils.createJTSLineString(coordinates);

    }

    /* (non-Javadoc)

     * @see org.gvsig.fmap.geom.primitive.Line#toPoints()

     */

    public MultiPoint toPoints() throws GeometryException {

        MultiPoint multiPoint = new MultiPoint3D();

        Coordinate[] coordinates = getJTS().getCoordinates();

        multiPoint.ensureCapacity(coordinates.length);

        for (int i = 0; i < coordinates.length; i++) {

            multiPoint.addPoint(new Point3D(coordinates[i]));

        }

        return multiPoint;

    }

    /* (non-Javadoc)

     * @see org.gvsig.fmap.geom.primitive.Line#toLines()

     */

    public MultiLine toLines() throws GeometryException {

        MultiLine multiLine = new MultiLine3D();

        Line line = new Line3D(getJTS().getCoordinates());

        multiLine.addPrimitive(line);

        return multiLine;

    }

    /* (non-Javadoc)

     * @see org.gvsig.fmap.geom.primitive.Line#toPolygons()

     */

    public MultiPolygon toPolygons() throws GeometryException {

        MultiPolygon multiPolygon = new MultiPolygon3D();

        Polygon polygon = new Polygon3D(getJTS().getCoordinates());

        multiPolygon.addPrimitive(polygon);

        return multiPolygon;

    }

}