Application: gvSIG desktop » gvSIG geoprocess

/**

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2012 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.geoprocess.algorithm.intersection;

import es.unex.sextante.core.Sextante;

import es.unex.sextante.dataObjects.IVectorLayer;

import es.unex.sextante.exceptions.GeoAlgorithmExecutionException;

import es.unex.sextante.exceptions.RepeatedParameterNameException;

import es.unex.sextante.outputs.OutputVectorLayer;

import org.gvsig.fmap.dal.exception.DataException;

import org.gvsig.fmap.dal.feature.FeatureSet;

import org.gvsig.fmap.dal.feature.FeatureStore;

import org.gvsig.fmap.dal.feature.FeatureType;

import org.gvsig.geoprocess.lib.sextante.AbstractSextanteGeoProcess;

import org.gvsig.geoprocess.lib.sextante.dataObjects.FlyrVectIVectorLayer;

/**

 * Intersection algorithm

 * @author <a href="mailto:nachobrodin@gmail.com">Nacho Brodin</a>

 */

public class IntersectionAlgorithm extends AbstractSextanteGeoProcess {

        public static final String  RESULT    = "RESULT";

        public static final String  LAYER     = "LAYER";

        public static final String  INTER     = "INTER";

        public static final String  CHECK     = "CHECK";

        /*

         * (non-Javadoc)

         * @see es.unex.sextante.core.GeoAlgorithm#defineCharacteristics()

         */

        public void defineCharacteristics() {

        setName(getTranslation("Intersection"));

        setGroup(getTranslation("basic_vect_algorithms"));

        // setGeneratesUserDefinedRasterOutput(false);

                try {

                        m_Parameters.addInputVectorLayer(LAYER, 

                getTranslation("Input_layer"),

                                                                                                IVectorLayer.SHAPE_TYPE_WRONG, 

                                                                                                true);

                        m_Parameters.addInputVectorLayer(INTER, 

                getTranslation("Overlays_layer"),

                                                                                                IVectorLayer.SHAPE_TYPE_WRONG, 

                                                                                                true);

            m_Parameters.addBoolean(CHECK,

                getTranslation("Selected_geometries"), false);

                } catch (RepeatedParameterNameException e) {

                        Sextante.addErrorToLog(e);

                }

                addOutputVectorLayer(RESULT,

 getTranslation("Intersection"),

                                                                OutputVectorLayer.SHAPE_TYPE_UNDEFINED);

        }

        /*

         * (non-Javadoc)

         * @see es.unex.sextante.core.GeoAlgorithm#processAlgorithm()

         */

        public boolean processAlgorithm() throws GeoAlgorithmExecutionException {

                IVectorLayer inter = m_Parameters.getParameterValueAsVectorLayer(INTER);

                IVectorLayer layer = m_Parameters.getParameterValueAsVectorLayer(LAYER);

                boolean selectedGeom = m_Parameters.getParameter(CHECK).getParameterValueAsBoolean();

                try {

                        computesIntersection(layer, inter, layer.getShapeType(), selectedGeom);

                } catch (DataException e) {

                        Sextante.addErrorToLog(e);

                        return false;

                }

                return true;

        }

        /**

         * Builds a layer with the intersection between the input layer and the templateGeometry

         * @param layer

         *        Input layer

         * @param templateGeometry

         * @param shapeType

         *        Output shape type

         * @param selectedGeom

         *        If it's true only selected geometries will be computed

         */

        private void computesIntersection(IVectorLayer layer,

                                                                IVectorLayer overlay,

                                                                int shapeType, 

                                                                boolean selectedGeom) throws DataException {

                FeatureStore storeLayer = null;

                FeatureStore storeOverlay = null;

                if(layer instanceof FlyrVectIVectorLayer && 

                        overlay instanceof FlyrVectIVectorLayer) {

                        storeLayer = ((FlyrVectIVectorLayer)layer).getFeatureStore();

                        storeOverlay = ((FlyrVectIVectorLayer)overlay).getFeatureStore();

                } else

                        return;

                FeatureSet features1 = null;

                features1 = storeLayer.getFeatureSet();

                FeatureType featureType1 = features1.getDefaultFeatureType();

                FeatureSet features2 = null;

                features2 = storeOverlay.getFeatureSet();

                FeatureType featureType2 = features2.getDefaultFeatureType();

                FeatureStore outFeatStore = null;

                //Si las dos capas son de lineas la resultante es de l?neas

        if (isLine(storeLayer) && isLine(storeOverlay))

            outFeatStore =

                buildOutPutStoreFromUnion(featureType1, featureType2,

                    IVectorLayer.SHAPE_TYPE_POINT,

                    getTranslation("Intersection"), RESULT);

        else

        // Si alguna de las dos capas es de puntos la resultante tambi?n lo es

        if (isPoint(storeLayer) || isPoint(storeOverlay))

            outFeatStore =

                buildOutPutStoreFromUnion(featureType1, featureType2,

                    IVectorLayer.SHAPE_TYPE_POINT,

                    getTranslation("Intersection"), RESULT);

        else

            outFeatStore =

                buildOutPutStoreFromUnion(featureType1, featureType2,

                    shapeType, getTranslation("Intersection"), RESULT);

                IntersectionOperation operation = new IntersectionOperation(storeOverlay);

                operation.computesGeometryOperation(storeLayer, outFeatStore, attrNames, selectedGeom, true);

        }

}