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;

import com.vividsolutions.jts.algorithm.CGAlgorithms;

import com.vividsolutions.jts.geom.Coordinate;

import com.vividsolutions.jts.io.geojson.GeoJsonWriter;

import com.vividsolutions.jts.operation.buffer.BufferParameters;

import com.vividsolutions.jts.operation.distance.DistanceOp;

import com.vividsolutions.jts.operation.overlay.snap.GeometrySnapper;

import com.vividsolutions.jts.operation.valid.IsValidOp;

import com.vividsolutions.jts.operation.valid.TopologyValidationError;

import java.awt.Rectangle;

import java.awt.Shape;

import java.awt.geom.AffineTransform;

import java.awt.geom.Rectangle2D;

import java.util.logging.Level;

import org.apache.commons.codec.binary.Hex;

import org.apache.commons.lang3.StringUtils;

import org.cresques.cts.IProjection;

import org.gvsig.fmap.crs.CRSFactory;

import org.gvsig.fmap.geom.Geometry;

import org.gvsig.fmap.geom.GeometryLocator;

import org.gvsig.fmap.geom.GeometryManager;

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

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

import org.gvsig.fmap.geom.complex.Complex;

import org.gvsig.fmap.geom.exception.CreateGeometryException;

import org.gvsig.fmap.geom.jts.operation.towkb.OGCWKBEncoder;

import org.gvsig.fmap.geom.jts.operation.towkb.PostGISEWKBEncoder;

import org.gvsig.fmap.geom.jts.operation.towkt.EWKTWriter;

import org.gvsig.fmap.geom.jts.primitive.Envelope2D;

import org.gvsig.fmap.geom.jts.primitive.Envelope3D;

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

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

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

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

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

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

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

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

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

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

import org.gvsig.fmap.geom.type.GeometryType;

import org.slf4j.Logger;

import org.slf4j.LoggerFactory;

/**

 * @author gvSIG Team

 *

 */

@SuppressWarnings("UseSpecificCatch")

public abstract class AbstractGeometry implements GeometryJTS {

    protected static final Logger LOGGER = LoggerFactory.getLogger(AbstractGeometry.class);

    /**

     *

     */

    private static final long serialVersionUID = 4999326772576222293L;

    private GeometryType geometryType;

    private IProjection projection;

    /**

     *

     * @param type

     * @param subtype

     */

    public AbstractGeometry(int type, int subtype) {

        try {

            this.geometryType = GeometryLocator.getGeometryManager().getGeometryType(type, subtype);

        } catch (Exception e) {

            // TODO: Ver de crear una excepcion para esto o si hay alguna en el API.

            throw new RuntimeException(e);

        }

    }

    public void setGeometryType(GeometryType geometryType) {

        this.geometryType = geometryType;

    }

    @Override

    public GeometryType getGeometryType() {

        return geometryType;

    }

    protected GeometryManager getManager() {

        return GeometryLocator.getGeometryManager();

    }

    /*

     * (non-Javadoc)

     *

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

     */

    @Override

    public boolean contains(Geometry geometry) throws GeometryOperationNotSupportedException,

            GeometryOperationException {

        if (!(geometry instanceof GeometryJTS)) {

            return false;

        }

        com.vividsolutions.jts.geom.Geometry theJTS = this.getJTS();

        if(this instanceof Aggregate){

            theJTS = theJTS.union();

        }

        com.vividsolutions.jts.geom.Geometry otherJTS = ((GeometryJTS)geometry).getJTS();

        if(geometry instanceof Aggregate){

            otherJTS = otherJTS.union();

        }

        return theJTS.contains(otherJTS);

    }

    /*

     * (non-Javadoc)

     *

     * @see java.lang.Comparable#compareTo(java.lang.Object)

     */

    @Override

    public int compareTo(Object o) {

        return getJTS().compareTo(((GeometryJTS) o).getJTS());

    }

    /*

     * (non-Javadoc)

     *

     * @see java.awt.Shape#contains(double, double)

     */

    @Override

    public boolean contains(double x, double y) {

        notifyDeprecated("Calling deprecated method of geometry contains from shape interface");

        Shape shp = getShape();

        return shp.contains(x, y);

    }

    /*

     * (non-Javadoc)

     *

     * @see java.awt.Shape#intersects(double, double, double, double)

     */

    @Override

    public boolean intersects(double x, double y, double w, double h) {

        notifyDeprecated("Calling deprecated method of geometry intersects from shape interface");

        Shape shp = this.getShape();

        return shp.intersects(x, y, w, h);

    }

    /*

     * (non-Javadoc)

     *

     * @see java.awt.Shape#contains(double, double, double, double)

     */

    @Override

    public boolean contains(double x, double y, double w, double h) {

        notifyDeprecated("Calling deprecated method of geometry contains from shape interface");

        Shape shp = this.getShape();

        return shp.contains(x, y, w, h);

    }

    /*

     * (non-Javadoc)

     *

     * @see java.awt.Shape#contains(java.awt.geom.Point2D)

     */

    @Override

    public boolean contains(java.awt.geom.Point2D p) {

        notifyDeprecated("Calling deprecated method of geometry contains from shape interface");

        Shape shp = this.getShape();

        return shp.contains(p);

    }

    /*

     * (non-Javadoc)

     *

     * @see java.awt.Shape#contains(java.awt.geom.Rectangle2D)

     */

    @Override

    public boolean contains(Rectangle2D r) {

        notifyDeprecated("Calling deprecated method of geometry contains from shape interface");

        Shape shp = this.getShape();

        return shp.contains(r);

    }

    /*

     * (non-Javadoc)

     *

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

     */

    @Override

    public double distance(org.gvsig.fmap.geom.Geometry other) throws GeometryOperationNotSupportedException,

            GeometryOperationException {

        return getJTS().distance(((GeometryJTS) other).getJTS());

    }

    /*

     * (non-Javadoc)

     *

     * @see

     * org.gvsig.fmap.geom.Geometry#isWithinDistance(org.gvsig.fmap.geom.Geometry

     * , double)

     */

    @Override

    public boolean isWithinDistance(org.gvsig.fmap.geom.Geometry other, double distance)

            throws GeometryOperationNotSupportedException, GeometryOperationException {

        return com.vividsolutions.jts.operation.distance.DistanceOp.isWithinDistance(this.getJTS(),

                ((GeometryJTS) other).getJTS(), distance);

    }

    /*

     * (non-Javadoc)

     *

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

     */

    @Override

    public boolean overlaps(org.gvsig.fmap.geom.Geometry geometry) throws GeometryOperationNotSupportedException,

            GeometryOperationException {

        // TODO: this method can be implemented throw invokeOperation

        return getJTS().overlaps(((GeometryJTS) geometry).getJTS());

    }

    @Override

    public boolean coveredBy(Geometry geometry) throws GeometryOperationNotSupportedException,

            GeometryOperationException {

        return getJTS().coveredBy(((GeometryJTS) geometry).getJTS());

    }

    @Override

    public boolean covers(Geometry geometry) throws GeometryOperationNotSupportedException, GeometryOperationException {

        // TODO: this method can be implemented throw invokeOperation

        return getJTS().covers(((GeometryJTS) geometry).getJTS());

    }

    @Override

    public boolean crosses(Geometry geometry) throws GeometryOperationNotSupportedException, GeometryOperationException {

        // TODO: this method can be implemented throw invokeOperation

        return getJTS().crosses(((GeometryJTS) geometry).getJTS());

    }

    @Override

    public boolean intersects(Geometry geometry) throws GeometryOperationNotSupportedException,

            GeometryOperationException {

        return getJTS().intersects(((GeometryJTS) geometry).getJTS());

    }

    @Override

    public boolean touches(Geometry geometry) throws GeometryOperationNotSupportedException, GeometryOperationException {

        // TODO: this method can be implemented throw invokeOperation

        return getJTS().touches(((GeometryJTS) geometry).getJTS());

    }

    @Override

    public boolean disjoint(Geometry geometry) throws GeometryOperationNotSupportedException,

            GeometryOperationException {

        // TODO: this method can be implemented throw invokeOperation

        return getJTS().disjoint(((GeometryJTS) geometry).getJTS());

    }

    @Override

    public boolean within(Geometry geometry) throws GeometryOperationNotSupportedException, GeometryOperationException {

        // TODO: this method can be implemented throw invokeOperation

        return getJTS().within(((GeometryJTS) geometry).getJTS());

    }

    @Override

    public double area() throws GeometryOperationNotSupportedException, GeometryOperationException {

        return getJTS().getArea();

    }

    @Override

    public double perimeter() throws GeometryOperationNotSupportedException, GeometryOperationException {

        return getJTS().getLength();

    }

    /*

     * (non-Javadoc)

     *

     * @see org.gvsig.fmap.geom.Geometry#intersects(java.awt.geom.Rectangle2D)

     */

    @Override

    public boolean intersects(Rectangle2D r) {

        double x = r.getMinX();

        double y = r.getMinY();

        double w = r.getWidth();

        double h = r.getHeight();

        return fastIntersects(x, y, w, h);

    }

    /*

     * (non-Javadoc)

     *

     * @see org.gvsig.fmap.geom.Geometry#fastIntersects(double, double, double,

     * double)

     */

    @Override

    public boolean fastIntersects(double x, double y, double w, double h) {

        com.vividsolutions.jts.geom.Geometry rect; 

        com.vividsolutions.jts.geom.Coordinate[] coord = new com.vividsolutions.jts.geom.Coordinate[5];

        coord[0] = new com.vividsolutions.jts.geom.Coordinate(x, y);

        coord[1] = new com.vividsolutions.jts.geom.Coordinate(x + w, y);

        coord[2] = new com.vividsolutions.jts.geom.Coordinate(x + w, y + w);

        coord[3] = new com.vividsolutions.jts.geom.Coordinate(x, y + w);

        coord[4] = new com.vividsolutions.jts.geom.Coordinate(x, y);

        rect = JTSUtils.createJTSPolygon(coord);

        return getJTS().intersects(rect);

    }

    /*

     * (non-Javadoc)

     *

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

     */

    @Override

    public Envelope getEnvelope() {

        if (is3D()) {

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

            double minx = Double.POSITIVE_INFINITY;

            double miny = Double.POSITIVE_INFINITY;

            double minz = Double.POSITIVE_INFINITY;

            double maxx = Double.NEGATIVE_INFINITY;

            double maxy = Double.NEGATIVE_INFINITY;

            double maxz = Double.NEGATIVE_INFINITY;

            double x;

            double y;

            double z;

            for (Coordinate coordinate : coordinates) {

                x = coordinate.x;

                y = coordinate.y;

                z = coordinate.z;

                minx = Math.min(x, minx);

                miny = Math.min(y, miny);

                minz = Math.min(z, minz);

                maxx = Math.max(x, maxx);

                maxy = Math.max(y, maxy);

                maxz = Math.max(z, maxz);

            }

            if (minx <= maxx && miny <= maxy && minz <= maxz) {

                Point min = new Point3D(minx, miny, minz);

                Point max = new Point3D(maxx, maxy, maxz);

                return new Envelope3D(min, max, this.getProjection());

            }

            return new Envelope3D(this.getProjection());

        } else {

            com.vividsolutions.jts.geom.Envelope envelope = getJTS().getEnvelopeInternal();

            return new Envelope2D(envelope, this.getProjection());

        }

    }

    /*

     * (non-Javadoc)

     *

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

     */

    @Override

    public boolean isSimple() {

        return this.getJTS().isSimple();

    }

    /*

     * (non-Javadoc)

     *

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

     */

    @Override

    public boolean isCCW() throws GeometryOperationNotSupportedException, GeometryOperationException {

        return CGAlgorithms.isCCW(this.getJTS().getCoordinates());

    }

    /*

     * (non-Javadoc)

     *

     * @see org.gvsig.fmap.geom.Geometry#invokeOperation(int,

     * org.gvsig.fmap.geom.operation.GeometryOperationContext)

     */

    @Override

    public Object invokeOperation(int index, GeometryOperationContext ctx)

            throws GeometryOperationNotSupportedException, GeometryOperationException {

        return getManager().invokeOperation(index, this, ctx);

    }

    /*

     * (non-Javadoc)

     *

     * @see org.gvsig.fmap.geom.Geometry#invokeOperation(java.lang.String,

     * org.gvsig.fmap.geom.operation.GeometryOperationContext)

     */

    @Override

    public Object invokeOperation(String opName, GeometryOperationContext ctx)

            throws GeometryOperationNotSupportedException, GeometryOperationException {

        return getManager().invokeOperation(opName, this, ctx);

    }

    /*

     * (non-Javadoc)

     *

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

     */

    @Override

    public int getType() {

        return this.getGeometryType().getType();

    }

    @Override

    public Object convertTo(String format) throws GeometryOperationNotSupportedException, GeometryOperationException {

        if( StringUtils.isEmpty(format) ) {

            throw new IllegalArgumentException("Can't accept null as format name.");

        }

        format = format.trim().toLowerCase();

        switch(format) {

            case "jts":

                return this.getJTS();

            case "wkb":

                return this.convertToWKB();

            case "hexwkb":

                return this.convertToHexWKB();

            case "ewkb":

                return this.convertToEWKB();

            case "wkt":

                return this.convertToWKT();

            case "gml":

                return GMLUtils.geometry2GML(this);

            case "json":

            case "geojson":

                return this.convertToGeoJson();

            default:

                throw new IllegalArgumentException("Format '"+format+"' not supported");

        }

    }

    @Override

    public byte[] convertToWKB() throws GeometryOperationNotSupportedException, GeometryOperationException {

        try {

            return new OGCWKBEncoder().encode(this);

        } catch (Exception e) {

            throw new GeometryOperationException(e);

        }

    }

    @Override

    public String convertToHexWKB() throws GeometryOperationNotSupportedException, GeometryOperationException {

        try {

            byte[] bytes = new OGCWKBEncoder().encode(this);

            return Hex.encodeHexString(bytes);

        } catch (Exception e) {

            throw new GeometryOperationException(e);

        }

    }

    public String convertToGeoJson() throws GeometryOperationNotSupportedException, GeometryOperationException {

        try {

            GeoJsonWriter writer = new GeoJsonWriter();

            writer.setEncodeCRS(true);

            return writer.write(this.getJTS());

        } catch (Exception e) {

            throw new GeometryOperationException(e);

        }

    }

    @Override

    public String convertToHexEWKB() throws GeometryOperationNotSupportedException, GeometryOperationException {

        try {

            byte[] bytes = new PostGISEWKBEncoder().encode(this);

            return Hex.encodeHexString(bytes);

        } catch (Exception e) {

            throw new GeometryOperationException(e);

        }

    }

    @Override

    public byte[] convertToWKBQuietly() {

        try {

            return new OGCWKBEncoder().encode(this);

        } catch (Exception e) {

            return null;

        }

    }

    @Override

    public String convertToHexWKBQuietly() {

        try {

            byte[] bytes = new OGCWKBEncoder().encode(this);

            return Hex.encodeHexString(bytes);

        } catch (Exception e) {

            return null;

        }

    }

    @Override

    public byte[] convertToWKB(int srs) throws GeometryOperationNotSupportedException, GeometryOperationException {

        /*

             * No se sabe si la especificaci?n de OGC soporta SRS. OGCWKBEncoder no.

         */

        try {

            return new OGCWKBEncoder().encode(this);

        } catch (Exception e) {

            throw new GeometryOperationException(e);

        }

    }

    @Override

    public byte[] convertToWKBForcingType(int srs, int type) throws GeometryOperationNotSupportedException,

            GeometryOperationException {

        /*

             * No se sabe si la especificaci?n de OGC soporta SRS. OGCWKBEncoder no.

         */

        Geometry geom = this;

        if (this.getType() != type) {

            com.vividsolutions.jts.geom.Geometry jts = getJTS();

            jts = JTSUtils.convertTypes(jts, this.getType(), type);

            geom = JTSUtils.createGeometry(this.getProjection(), jts);

        }

        try {

            return new OGCWKBEncoder().encode(geom);

        } catch (Exception e) {

            throw new GeometryOperationException(e);

        }

    }

    @Override

    public byte[] convertToEWKB() throws GeometryOperationNotSupportedException, GeometryOperationException {

        try {

            return new PostGISEWKBEncoder().encode(this);

        } catch (Exception e) {

            throw new GeometryOperationException(e);

        }

    }

    @Override

    public byte[] convertToEWKB(int srs) throws GeometryOperationNotSupportedException, GeometryOperationException {

        /*

             * No se sabe si la especificaci?n de OGC soporta SRS. OGCWKBEncoder no.

         */

        try {

            return new PostGISEWKBEncoder().encode(this);

        } catch (Exception e) {

            throw new GeometryOperationException(e);

        }

    }

    @Override

    public byte[] convertToEWKBForcingType(int srs, int type) throws GeometryOperationNotSupportedException,

            GeometryOperationException {

        /*

             * No se sabe si la especificaci?n de OGC soporta SRS. OGCWKBEncoder no.

         */

        Geometry geom = this;

        if (this.getType() != type) {

            com.vividsolutions.jts.geom.Geometry jts = getJTS();

            jts = JTSUtils.convertTypes(jts, this.getType(), type);

            geom = JTSUtils.createGeometry(this.getProjection(), jts);

        }

        try {

            return new PostGISEWKBEncoder().encode(geom);

        } catch (Exception e) {

            throw new GeometryOperationException(e);

        }

    }

    @Override

    public String convertToWKT() throws GeometryOperationNotSupportedException, GeometryOperationException {

        int subType = getGeometryType().getSubType();

        EWKTWriter writer; // = null;

        switch (subType) {

            case Geometry.SUBTYPES.GEOM3D:

                writer = new EWKTWriter(3, false);

                break;

            case Geometry.SUBTYPES.GEOM2DM:

                writer = new EWKTWriter(3, true);

                break;

            case Geometry.SUBTYPES.GEOM3DM:

                writer = new EWKTWriter(4, true);

                break;

            default:

                writer = new EWKTWriter(2, false);

                break;

        }

        com.vividsolutions.jts.geom.Geometry jts = getJTS();

        return writer.write(jts);

    }

    @Override

    public String convertToWKTQuietly() {

        try {

            return this.convertToWKT();

        } catch(Exception ex) {

            return null;

        }

    }

    @Override

    public org.gvsig.fmap.geom.Geometry buffer(double distance) throws GeometryOperationNotSupportedException,

        GeometryOperationException {

        if( distance==0 ) {

          return this;

        }

        return JTSUtils.createGeometry(this.getProjection(), getJTS().buffer(distance));

    }

    @Override

    public org.gvsig.fmap.geom.Geometry buffer(double distance, int joinStyle, boolean capButt) throws GeometryOperationNotSupportedException,

        GeometryOperationException {

        if( distance==0 ) {

          return this;

        }

        int quadrantSegments = JTSUtils.calculateQuadrantSegments(joinStyle);

        return JTSUtils.createGeometry(this.getProjection(), getJTS().buffer(

                distance, 

                quadrantSegments, 

                capButt ? BufferParameters.CAP_FLAT : BufferParameters.CAP_ROUND

        ));

    }

    /*

     * (non-Javadoc)

     *

     * @see org.gvsig.fmap.geom.Geometry#snapTo(org.gvsig.fmap.geom.Geometry,

     * double)

     */

    @Override

    public org.gvsig.fmap.geom.Geometry snapTo(org.gvsig.fmap.geom.Geometry other, double snapTolerance)

            throws GeometryOperationNotSupportedException, GeometryOperationException {

        Geometry result; // = null;

        GeometrySnapper snapper = new GeometrySnapper(getJTS());

        com.vividsolutions.jts.geom.Geometry jts_result = snapper.snapTo(((GeometryJTS) other).getJTS(), snapTolerance);

        result = JTSUtils.createGeometry(this.getProjection(), jts_result);

        return result;

    }

    /*

     * (non-Javadoc)

     *

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

     */

    @Override

    public Point getInteriorPoint() throws GeometryOperationNotSupportedException, GeometryOperationException {

        com.vividsolutions.jts.geom.Geometry geometry = getJTS();

        com.vividsolutions.jts.geom.Point point = geometry.getInteriorPoint();

        Geometry result = JTSUtils.createGeometry(this.getProjection(), point);

        return (Point) result;

    }

    /*

     * (non-Javadoc)

     *

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

     */

    @Override

    public boolean isValid() {

        return getJTS().isValid();

    }

    @Override

    public ValidationStatus getValidationStatus() {

        DefaultValidationStatus status = new DefaultValidationStatus(ValidationStatus.VALID, null);

        com.vividsolutions.jts.geom.Geometry jtsgeom; // = null;

        try {

            jtsgeom = this.getJTS();

            IsValidOp validOp = new IsValidOp(jtsgeom);

            if (!validOp.isValid() ) {

                status.setValidationError(validOp.getValidationError());

            }

        } catch (Throwable ex) {

            status.setStatusCode(ValidationStatus.CURRUPTED);

            status.setMesage("The geometry is corrupted.");

            if (this instanceof OrientableSurface) {

                int vertices = ((OrientableSurface) this).getNumVertices();

                if (vertices < 3) {

                    status.setStatusCode(ValidationStatus.TOO_FEW_POINTS);

                    status.setMesage(TopologyValidationError.errMsg[TopologyValidationError.TOO_FEW_POINTS]);

                }

            } else if (this instanceof OrientableCurve) {

                int vertices = ((OrientableCurve) this).getNumVertices();

                if (vertices < 2) {

                    status.setStatusCode(ValidationStatus.TOO_FEW_POINTS);

                    status.setMesage(TopologyValidationError.errMsg[TopologyValidationError.TOO_FEW_POINTS]);

                }

            }

        }

        return status;

    }

    /*

     * (non-Javadoc)

     *

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

     */

    @Override

    public org.gvsig.fmap.geom.Geometry makeValid() {

        try {

            ValidationStatus vs = this.getValidationStatus();

            if (vs.isValid()) {

                return this;

            }

            Geometry g; // = null;

            switch (vs.getStatusCode()) {

                case Geometry.ValidationStatus.RING_SELF_INTERSECTION:

                case Geometry.ValidationStatus.SELF_INTERSECTION:

                    g = this.buffer(0);

                    if (g.isValid()) {

                        return g;

                    }

                    break;

                case Geometry.ValidationStatus.TOO_FEW_POINTS:

                    if (this instanceof OrientableCurve) {

                        int vertices = ((OrientableCurve) this).getNumVertices();

                        if (vertices < 2) {

                            return null; // new

                            // DefaultNullGeometry(this.getGeometryType());

                        }

                    }

                    if (this instanceof OrientableSurface) {

                        int vertices = ((OrientableSurface) this).getNumVertices();

                        if (vertices < 3) {

                            return null; // new

                            // DefaultNullGeometry(this.getGeometryType());

                        }

                    }

            }

        } catch (Exception ex) {

            return null;

        }

        return null;

    }

    /*

     * (non-Javadoc)

     *

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

     */

    @Override

    public Rectangle2D getBounds2D() {

        com.vividsolutions.jts.geom.Envelope envInternal = getJTS().getEnvelopeInternal();

        return new Rectangle2D.Double(envInternal.getMinX(), envInternal.getMinY(), envInternal.getWidth(),

                envInternal.getHeight());

    }

    /*

     * (non-Javadoc)

     *

     * @see java.awt.Shape#getBounds()

     */

    @Override

    public Rectangle getBounds() {

        return this.getShape().getBounds();

    }

    /*

     * (non-Javadoc)

     *

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

     */

    @Override

    public Shape getInternalShape() {

        return getShape();

    }

    @Override

    public void rotate(double radAngle, double basex, double basey) {

        AffineTransform at = new AffineTransform();

        at.rotate(radAngle, basex, basey);

        this.transform(at);

    }

    @Override

    public void move(double dx, double dy) {

        AffineTransform at = new AffineTransform();

        at.translate(dx, dy);

        this.transform(at);

    }

    @Override

    public void scale(Point basePoint, double sx, double sy) {

        AffineTransform at = new AffineTransform();

        at.setToTranslation(basePoint.getX(), basePoint.getY());

        at.scale(sx, sy);

        at.translate(-basePoint.getX(), -basePoint.getY());

        this.transform(at);

    }

    @Override

    public Geometry[] closestPoints(Geometry other) throws GeometryOperationNotSupportedException,

            GeometryOperationException {

        Point[] points;

        Coordinate[] jts_points = DistanceOp.nearestPoints(getJTS(), ((GeometryJTS) other).getJTS());

        points = new Point[jts_points.length];

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

            try {

                points[i] = JTSUtils.createPoint(this.getGeometryType(), this.getProjection(), jts_points[i]);

            } catch (CreateGeometryException e) {

                throw new GeometryOperationException(e);

            }

        }

        return (Geometry[]) points;

    }

    @Override

    public Geometry convexHull() throws GeometryOperationNotSupportedException, GeometryOperationException {

        return JTSUtils.createGeometry(this.getProjection(), getJTS().convexHull(), null);

    }

    @Override

    public Geometry difference(Geometry other) throws GeometryOperationNotSupportedException,

            GeometryOperationException {

        com.vividsolutions.jts.geom.Geometry otherJTS = ((GeometryJTS) other).getJTS();

        if(other instanceof Aggregate){

            otherJTS = otherJTS.union();

        }

        return JTSUtils.createGeometry(this.getProjection(), getJTS().difference(otherJTS), null);

    }

    @Override

    public Geometry intersection(Geometry other) throws GeometryOperationNotSupportedException,

            GeometryOperationException {

        com.vividsolutions.jts.geom.Geometry otherJTS = ((GeometryJTS) other).getJTS();

        if(other instanceof Aggregate){

            otherJTS = otherJTS.union();

        }

        return JTSUtils.createGeometry(this.getProjection(), getJTS().intersection(otherJTS), null);

    }

    @Override

    public Geometry union(Geometry other) throws GeometryOperationNotSupportedException, GeometryOperationException {

        try {

            com.vividsolutions.jts.geom.Geometry jts = getJTS();

            com.vividsolutions.jts.geom.Geometry otherJts = ((GeometryJTS)other).getJTS();

            if(jts.isValid() && otherJts.isValid()){

                return JTSUtils.createGeometry(this.getProjection(), jts.union(otherJts), this.getGeometryType());

            }

            MultiPrimitive geom = this.getManager().createMultiPrimitive(geometryType);

            geom.addPrimitives(this);

            geom.addPrimitives(other);

            return geom;

        } catch (CreateGeometryException ex) {

            throw new GeometryOperationException(

                    this.getGeometryType().getType(), 

                    this.getManager().getGeometryOperationCode(Geometry.OPERATIONS.UNION), 

                    ex);

        }

    }

    @Override

    public org.gvsig.fmap.geom.primitive.Point centroid() throws GeometryOperationNotSupportedException,

            GeometryOperationException {

        try {

            return JTSUtils.createPoint(this.getGeometryType(), this.getProjection(), getJTS().getCentroid().getCoordinate());

        } catch (CreateGeometryException e) {

            throw new GeometryOperationException(e);

        }

    }

    protected void notifyDeprecated(String message) {

        LOGGER.info(message);

    }

    @Override

    public boolean ensureOrientation(boolean ccw) throws GeometryOperationNotSupportedException, GeometryOperationException {

        if (ccw != isCCW()) {

            flip();

            return true;

        }

        return false;

    }

    /* (non-Javadoc)

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

     */

    @Override

    public boolean out(Geometry geometry) throws GeometryOperationNotSupportedException, GeometryOperationException {

        GeometryJTS otherJtsGeom = (GeometryJTS) geometry;

        return (!contains(otherJtsGeom) && !intersects(otherJtsGeom));

    }

    @Override

    public boolean equals(Object obj) {

        if (obj instanceof GeometryJTS) {

            return this.getJTS().equals(((GeometryJTS) obj).getJTS());

        }

        return false;

    }

    @Override

    public String toString() {

        return this.getGeometryType().getFullName();

    }

    @Override

    public IProjection getProjection() {

        return this.projection;

    }

    @Override

    public void setProjectionIffNull(IProjection projection) {

        if (this.projection == null) {

            this.projection = projection;

        }

    }

    @Override

    public void setProjection(IProjection projection) {

        this.projection = projection;

    }

    @Override

    public void setProjection(String projection) {

        IProjection proj = CRSFactory.getCRS("EPSG:4326");

        this.setProjection(proj);

    }

    @Override

    @SuppressWarnings("CloneDoesntCallSuperClone")

    public Geometry clone() throws CloneNotSupportedException {

        return this.cloneGeometry();

    }

    @Override

    public Geometry boundary() {

        return JTSUtils.createGeometry(this.getProjection(), getJTS().getBoundary(), null);

    }

}