Application: gvSIG desktop

/**

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2013 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 3

 * 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.primitive.impl;

import java.awt.geom.Rectangle2D;

import java.util.ArrayList;

import java.util.List;

import com.vividsolutions.jts.geom.GeometryFactory;

import com.vividsolutions.jts.geom.LineString;

import com.vividsolutions.jts.geom.LinearRing;

import com.vividsolutions.jts.geom.MultiPolygon;

import com.vividsolutions.jts.geom.Polygon;

import com.vividsolutions.jts.geom.PrecisionModel;

import com.vividsolutions.jts.geom.TopologyException;

import com.vividsolutions.jts.precision.GeometryPrecisionReducer;

import org.cresques.cts.IProjection;

import org.gvsig.fmap.geom.Geometry;

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

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

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

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

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

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

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

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

import org.gvsig.fmap.geom.util.Converter;

/**

 * Pol�gono 2D.

 *

 * @author Vicente Caballero Navarro

 */

public class Surface2D extends OrientableSurface2D implements Surface {

        private static final long serialVersionUID = -8448256617197415743L;

        /**

         * The constructor with the GeometryType like and argument 

         * is used by the {@link GeometryType}{@link #create()}

         * to create the geometry

         * @param type

         * The geometry type

         */

        public Surface2D(GeometryType geomType) {

                super(geomType);

        }

        /**

         * Constructor used in the {@link Geometry#cloneGeometry()} method.

         * @param geomType

         * @param id

         * @param projection

         * @param gpx

         */

        Surface2D(GeometryType geomType, String id, IProjection projection, GeneralPathX gpx) {

                super(geomType, id, projection, gpx);

        }        

        /*

         * (non-Javadoc)

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

         */

        public int getShapeType() {

                return TYPES.SURFACE;

        }

        /*

         * (non-Javadoc)

         * @see org.gvsig.fmap.geom.primitive.FShape#cloneFShape()

         */

        public FShape cloneFShape() {

                return new Surface2D(getGeometryType(), id, projection, (GeneralPathX) gp.clone());

        }

        /*

         * (non-Javadoc)

         * @see org.gvsig.fmap.geom.primitive.impl.OrientablePrimitive2D#intersects(java.awt.geom.Rectangle2D)

         */

        public boolean intersects(Rectangle2D r) {

                return gp.intersects(r);

        }

        public SurfaceAppearance getSurfaceAppearance() {

                // TODO Auto-generated method stub

                return null;

        }

        public void setSurfaceAppearance(SurfaceAppearance app) {

                // TODO Auto-generated method stub

        }

    public Geometry union(Geometry other)

        throws GeometryOperationNotSupportedException,

        GeometryOperationException {

        com.vividsolutions.jts.geom.Geometry this_jts, other_jts, union_jts;

        Geometry resp = null;

        try {

            this_jts = getJTS();

            other_jts = Converter.geometryToJts(other);

            union_jts = this_jts.union(other_jts);

            if (union_jts instanceof Polygon) {

                double tole = union_jts.getArea();

                tole = tole * 1e-12;

                union_jts = removeInvalidHoles((Polygon) union_jts, tole);

            } else {

                if (union_jts instanceof MultiPolygon) {

                    double tole = union_jts.getArea();

                    tole = tole * 1e-12;

                    union_jts =

                        removeInvalidHoles((MultiPolygon) union_jts, tole);

                }

            }

            resp = Converter.jtsToGeometry(union_jts);

            return resp;

        } catch (CreateGeometryException e) {

            throw new GeometryOperationException(e);

        } catch (TopologyException e) {

            /*

             * If JTS throws this, we'll try to simplify them

             */

            try {

                this_jts = getJTS();

                other_jts = Converter.geometryToJts(other);

                PrecisionModel pm = new PrecisionModel(1000000);

                GeometryPrecisionReducer gpr = new GeometryPrecisionReducer(pm);

                this_jts = gpr.reduce(this_jts);

                other_jts = gpr.reduce(other_jts);

                union_jts = this_jts.union(other_jts);

                if (union_jts instanceof Polygon) {

                    double tole = union_jts.getArea();

                    tole = tole * 1e-12;

                    union_jts = removeInvalidHoles((Polygon) union_jts, tole);

                } else {

                    if (union_jts instanceof MultiPolygon) {

                        double tole = union_jts.getArea();

                        tole = tole * 1e-12;

                        union_jts =

                            removeInvalidHoles((MultiPolygon) union_jts, tole);

                    }

                }

                resp = Converter.jtsToGeometry(union_jts);

                return resp;

            } catch (Exception exc) {

                throw new GeometryOperationException(exc);

            }

        }

    }

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

        return Converter.geometryToJts(this);

    }

    private MultiPolygon removeInvalidHoles(MultiPolygon mpo, double tol) {

        GeometryFactory gf = new GeometryFactory(mpo.getPrecisionModel());

        int npo = mpo.getNumGeometries();

        Polygon[] pos = new Polygon[npo];

        for (int i=0; i<npo; i++) {

            pos[i] = removeInvalidHoles((Polygon) mpo.getGeometryN(i), tol);

        }

        return gf.createMultiPolygon(pos);

    }

    private Polygon removeInvalidHoles(Polygon po, double tol) {

        GeometryFactory gf = new GeometryFactory(po.getPrecisionModel());

        int nholes = po.getNumInteriorRing();

        List validholes = new ArrayList();

        for (int i=0; i<nholes; i++) {

            LineString ls = po.getInteriorRingN(i);

            LinearRing lr = toLinearRing(ls, gf);

            if (getLinearRingArea(lr, gf) > tol) {

                validholes.add(lr);

            }

        }

        if (validholes.size() < nholes) {

            LinearRing[] holes = (LinearRing[]) validholes.toArray(new LinearRing[0]);

            return gf.createPolygon(

                toLinearRing(po.getExteriorRing(), gf), holes);

        } else {

            return po;

        }

    }

    private double getLinearRingArea(LinearRing lr, GeometryFactory gf) {

        Polygon po = gf.createPolygon(lr);

        double resp = po.getArea(); 

        return resp;

    }

    private LinearRing toLinearRing(LineString ls, GeometryFactory gf) {

        return gf.createLinearRing(ls.getCoordinateSequence());

    }

}