| gvSIG desktop 1

/*

 * Created on 10-abr-2006

 *

 * gvSIG. Sistema de Informaci?n Geogr?fica de la Generalitat Valenciana

 *

 * Copyright (C) 2004 IVER T.I. and Generalitat Valenciana.

 *

 * 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., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 *

 * For more information, contact:

 *

 *  Generalitat Valenciana

 *   Conselleria d'Infraestructures i Transport

 *   Av. Blasco Ib??ez, 50

 *   46010 VALENCIA

 *   SPAIN

 *

 *      +34 963862235

 *   gvsig@gva.es

 *      www.gvsig.gva.es

 *

 *    or

 *

 *   IVER T.I. S.A

 *   Salamanca 50

 *   46005 Valencia

 *   Spain

 *

 *   +34 963163400

 *   dac@iver.es

 */

/* CVS MESSAGES:

 *

 * $Id:

 * $Log:

 */

package org.gvsig.topology.topologyrules.jtsisvalidrules;

import org.gvsig.fmap.core.FGeometryUtil;

import org.gvsig.jts.JtsUtil;

import org.gvsig.jts.SnapLineStringSelfIntersectionChecker;

import org.gvsig.topology.AbstractTopologyRule;

import org.gvsig.topology.Messages;

import org.gvsig.topology.Topology;

import org.gvsig.topology.TopologyError;

import org.gvsig.topology.TopologyRuleDefinitionException;

import org.gvsig.topology.WrongLyrForTopologyException;

import com.hardcode.gdbms.driver.exceptions.ReadDriverException;

import com.iver.cit.gvsig.fmap.core.FMultiPoint2D;

import com.iver.cit.gvsig.fmap.core.IFeature;

import com.iver.cit.gvsig.fmap.layers.FLyrVect;

import com.iver.utiles.XMLEntity;

import com.vividsolutions.jts.geom.Coordinate;

import com.vividsolutions.jts.geom.Geometry;

import com.vividsolutions.jts.geom.GeometryCollection;

import com.vividsolutions.jts.geom.LineString;

import com.vividsolutions.jts.geom.MultiPolygon;

import com.vividsolutions.jts.geom.Polygon;

public class PolygonMustNotHaveSelfIntersectedRings extends AbstractTopologyRule{

        private double snapTolerance;

        public PolygonMustNotHaveSelfIntersectedRings(Topology topology, FLyrVect lyr, double snapTolerance){

                super(topology, lyr);

                this.snapTolerance = snapTolerance;

        }

        public PolygonMustNotHaveSelfIntersectedRings(){}

        public PolygonMustNotHaveSelfIntersectedRings(FLyrVect lyr, double snapTolerance){

                this(null, lyr, snapTolerance);

        }

        public String getName() {

                return Messages.getText("POLYGON_MUST_NOT_HAVE_INTERSECTED_RINGS");

        }

        public void checkPreconditions() throws TopologyRuleDefinitionException {

                int shapeType;

                try {

                        shapeType = this.originLyr.getShapeType();

                        int numDimensions = FGeometryUtil.getDimensions(shapeType);

                        if(numDimensions != 2)

                                throw new WrongLyrForTopologyException("MustNotHaveSelfIntersectedRings solo aplica sobre capas de dimension 2");

                } catch (ReadDriverException e) {

                        throw new TopologyRuleDefinitionException(

                                        "Error al tratar de verificar el tipo de geometria");

                }

        }

        public void validateFeature(IFeature feature) {

                Geometry jtsGeo = feature.getGeometry().toJTSGeometry();

                if (jtsGeo instanceof Polygon) {

                        Polygon polygon = (Polygon) jtsGeo;

                        checkPolygon(polygon, feature);

                } else if (jtsGeo instanceof MultiPolygon) {

                        MultiPolygon multiPoly = (MultiPolygon) jtsGeo;

                        for (int i = 0; i < multiPoly.getNumGeometries(); i++) {

                                Polygon polygon = (Polygon) multiPoly.getGeometryN(i);

                                checkPolygon(polygon, feature);

                        }

                } else if (jtsGeo instanceof GeometryCollection) {

                        MultiPolygon multiPoly = JtsUtil

                                        .convertIfPossible((GeometryCollection) jtsGeo);

                        for (int i = 0; i < multiPoly.getNumGeometries(); i++) {

                                Polygon polygon = (Polygon) multiPoly.getGeometryN(i);

                                checkPolygon(polygon, feature);

                        }

                }

        }

        private  void checkPolygon(Polygon polygon, IFeature feature){

                LineString shell = polygon.getExteriorRing();

                SnapLineStringSelfIntersectionChecker checker =

                        new SnapLineStringSelfIntersectionChecker(shell, snapTolerance);

                if(checker.hasSelfIntersections()){

                        Coordinate[] selfIntersections = checker.getSelfIntersections();

                        addError(selfIntersections, feature);

                }

                int numHoles = polygon.getNumInteriorRing();

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

                        LineString hole = polygon.getInteriorRingN(i);

                        checker = new SnapLineStringSelfIntersectionChecker(hole, snapTolerance);

                        if(checker.hasSelfIntersections()){

                                Coordinate[] selfIntersections = checker.getSelfIntersections();

                                addError(selfIntersections, feature);

                        }//if

                }//for

        }

        private void addError(Coordinate[] selfIntersections, IFeature feature){

                double[] x = new double[selfIntersections.length];

                double[] y = new double[selfIntersections.length];

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

                        x[i] = selfIntersections[i].x;

                        y[i] = selfIntersections[i].y;

                }

                FMultiPoint2D errorGeo = new FMultiPoint2D(x, y);

                TopologyError error = new TopologyError(errorGeo, this, feature, topology );

                addTopologyError(error);

        }

        public XMLEntity getXMLEntity(){

                XMLEntity xml = super.getXMLEntity();

                xml.putProperty("snapTolerance", snapTolerance);

                return xml;

        }

        public void setXMLEntity(XMLEntity xml){

                super.setXMLEntity(xml);

                if(xml.contains("snapTolerance")){

                        snapTolerance = xml.getDoubleProperty("snapTolerance");

                }

        }

}