Application: gvSIG desktop » gvSIG geoprocess

org.gvsig.geoprocess.algorithm.intersection.IntersectionLibrary

#

# 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.

#

basic_vect_algorithms=Capas vectoriales

Input_layer=Capa de entrada

Overlays_layer=Capa de intersecci?n

Selected_geometries=Geometrias seleccionadas

Intersection=Intersecci?n

Selected_geometries_input_layer_inters=Geom. seleccionadas (Capa entrada)

Selected_geometries_overlay_layer_inters=Geom. seleccionadas (Capa intersecci?n)

problems_with_some_geometries=El proceso ha terminado pero ocurri? un error en alguna geometr?a. Es posible que el resultado no sea del todo exacto. Consulte el log.

#

# 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.

#

basic_vect_algorithms=Vector layers tools

Input_layer=Input cover

Overlays_layer=Intersection cover

Selected_geometries=Selected geometries

Intersection=Intersection

Selected_geometries_input_layer_inters=Selected features (Input cover)

Selected_geometries_overlay_layer_inters=Selected features (Intersection cover)

problems_with_some_geometries=The process has ended but an error happened in some geometry. Maybe the result not be totally accuracy. See the log.

<?xml version='1.0' encoding='ISO-8859-1' standalone='yes' ?>

<!--

    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.

-->

	<help>

		<element name="DESCRIPTION" text='Este geoproceso opera sobre dos capas, la "capa de entrada" y la "capa de solape", las geometr&#237;as de estas capas pueden ser de pol&#237;gonos, l&#237;neas o puntos.&#10;&#10;Para cada geometr&#237;a de la capa de entrada, calcula la interseccion con las diferentes geometr&#237;as de la capa de solape, originando un nuevo elemento por cada intersecci&#243;n. Este elemento tomar&#225; todos los atributos alfanum&#233;ricos de las geometr&#237;as que lo originaron (de entrada y solape). Por este motivo (modela zonas del espacio que cumplen la condici&#243;n de pertenecer a los dos pol&#237;gonos que lo han originado) a este geoproceso se le conoce como operador AND espacial.&#10;&#10;Un ejemplo de aplicaci&#243;n de este geoproceso ser&#237;a para, dada una capa de usos del suelo Corine 2000, y una capa del mapa geol&#243;gico nacional, obtener una capa de pol&#237;gonos con informaci&#243;n homog&#233;nea de uso del suelo y material geol&#243;gico.' description="Descripci&#243;n" type="0">

			<image description="" file="intersectdesc.png">

			</image>

		</element>

		<element name="ADDITIONAL_INFO" text="" description="Informaci&#243;n adicional" type="0">

		</element>

		<element name="EXTENSION_AUTHOR" text="Nacho Brodin" description="Algoritmo creado por" type="0">

		</element>

		<element name="HELP_AUTHOR" text="" description="Ayuda creada por" type="0">

		</element>

		<element name="USER_NOTES" text="" description="Notas de usuario" type="0">

		</element>

		<element name="LAYER" text="" description="Capa de entrada" type="3">

		</element>

		<element name="INTER" text="" description="Capa de revestimiento" type="3">

		</element>

		<element name="CHECK" text="" description="Geometrias seleccionadas" type="3">

		</element>

		<element name="OUTPUT_DESCRIPTION" text="" description="Descripci&#243;n" type="2">

		</element>

		<element name="RESULT" text="" description="Intersecci&#243;n" type="2">

		</element>

	</help>

<?xml version='1.0' encoding='ISO-8859-1' standalone='yes' ?>

<?xml version='1.0' encoding='ISO-8859-1' standalone='yes' ?>

<!--

    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.

-->

	<help>

		<element name="DESCRIPTION" text='This geoprocess works with two layers: the input layer and the overlay layer.&#10;&#10;For each geometry in the input layer, it computes its intersections with the geometries of the overlay layer. For each intersection, it adds a new feature with the intersection as geometry and all attributes of the features which originates it. It is called "Spatial AND", because features of the result layer models a space common to both layers.' description="Descripci&#243;n" type="0">

			<image description="" file="intersectdesc.png">

			</image>

		</element>

		<element name="ADDITIONAL_INFO" text="" description="Informaci&#243;n adicional" type="0">

		</element>

		<element name="EXTENSION_AUTHOR" text="Nacho Brodin" description="Algoritmo creado por" type="0">

		</element>

		<element name="HELP_AUTHOR" text="" description="Ayuda creada por" type="0">

		</element>

		<element name="USER_NOTES" text="" description="Notas de usuario" type="0">

		</element>

		<element name="LAYER" text="" description="Capa de entrada" type="3">

		</element>

		<element name="INTER" text="" description="Capa de revestimiento" type="3">

		</element>

		<element name="CHECK" text="" description="Geometrias seleccionadas" type="3">

		</element>

		<element name="OUTPUT_DESCRIPTION" text="" description="Descripci&#243;n" type="2">

		</element>

		<element name="RESULT" text="" description="Intersecci&#243;n" type="2">

		</element>

	</help>

<?xml version='1.0' encoding='ISO-8859-1' standalone='yes' ?>

/**

 * 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 org.gvsig.geoprocess.algorithm.base.core.AlgorithmAbstractLibrary;

import org.gvsig.i18n.Messages;

import org.gvsig.tools.library.LibraryException;

/**

 * Initialization of IntersectionLibrary library.

 * 

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

 */

public class IntersectionLibrary extends AlgorithmAbstractLibrary {

    @Override

    protected void doInitialize() throws LibraryException {

    }

    @Override

    protected void doPostInitialize() throws LibraryException {

        Messages.addResourceFamily(

            "org.gvsig.geoprocess.algorithm.intersection.intersection",

            IntersectionLibrary.class.getClassLoader(),

            IntersectionLibrary.class.getClass().getName());

        registerGeoProcess(new IntersectionAlgorithm());

    }

}

/**

 * 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 javax.swing.JOptionPane;

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

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;

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.tools.namestranslator.NamesTranslator;

/**

 * Intersection algorithm

 * <UL>

 * <LI>Pol-Pol: 3 layers (polygon, point, line)</LI>

 * <LI>Pol-Line: 2 layers (point, line)</LI>

 * <LI>Pol-Point: 1 layer (point)</LI>

 * <LI>Line-Point: 1 layer (point)</LI>

 * <LI>Line-Line: 2 layers (point, line)</LI>

 * <LI>Point-Point: 1 layer (point)</LI>

 * </UL>

 *

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

 */

public class IntersectionAlgorithm extends AbstractSextanteGeoProcess {

    public static final String RESULT_POL = "RESULT_POL";

    public static final String RESULT_POINT = "RESULT_POINT";

    public static final String RESULT_LINE = "RESULT_LINE";

    public static final String LAYER = "LAYER";

    public static final String INTER = "INTER";

    public static final String SELECTGEOM_INPUT = "SELECTGEOM_INPUT";

    public static final String SELECTGEOM_OVERLAY = "SELECTGEOM_OVERLAY";

    /*

	 * (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(SELECTGEOM_INPUT,

                    getTranslation("Selected_geometries_input_layer_inters"), false);

            m_Parameters.addBoolean(SELECTGEOM_OVERLAY,

                    getTranslation("Selected_geometries_overlay_layer_inters"), false);

        } catch (RepeatedParameterNameException e) {

            Sextante.addErrorToLog(e);

        }

        addOutputVectorLayer(RESULT_POL, getTranslation("Intersection_polygon"),

                OutputVectorLayer.SHAPE_TYPE_POLYGON);

        addOutputVectorLayer(RESULT_LINE, getTranslation("Intersection_line"),

                OutputVectorLayer.SHAPE_TYPE_LINE);

        addOutputVectorLayer(RESULT_POINT, getTranslation("Intersection_point"),

                OutputVectorLayer.SHAPE_TYPE_POINT);

    }

    /*

	 * (non-Javadoc)

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

     */

    public boolean processAlgorithm() throws GeoAlgorithmExecutionException {

        if (existsOutPutFile(IntersectionAlgorithm.RESULT_LINE, 0)) {

            throw new GeoAlgorithmExecutionException(getTranslation("file_exists"));

        }

        if (existsOutPutFile(IntersectionAlgorithm.RESULT_POINT, 0)) {

            throw new GeoAlgorithmExecutionException(getTranslation("file_exists"));

        }

        if (existsOutPutFile(IntersectionAlgorithm.RESULT_LINE, 0)) {

            throw new GeoAlgorithmExecutionException(getTranslation("file_exists"));

        }

        IVectorLayer inter = m_Parameters.getParameterValueAsVectorLayer(INTER);

        IVectorLayer layer = m_Parameters.getParameterValueAsVectorLayer(LAYER);

        boolean selectedGeomInput = m_Parameters.getParameter(SELECTGEOM_INPUT).getParameterValueAsBoolean();

        boolean selectedGeomOverlay = m_Parameters.getParameter(SELECTGEOM_OVERLAY).getParameterValueAsBoolean();

        boolean error = false;

        try {

            error = computesIntersection(layer, inter, layer.getShapeType(), selectedGeomInput, selectedGeomOverlay);

        } catch (DataException e) {

            Sextante.addErrorToLog(e);

            return false;

        }

        if (getTaskMonitor().isCanceled()) {

            return false;

        }

        if (error) {

            JOptionPane.showMessageDialog(null,

                    getTranslation("problems_with_some_geometries"), "Error",

                    JOptionPane.WARNING_MESSAGE);

        }

        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

     * @throws GeoAlgorithmExecutionException

     */

    private boolean computesIntersection(IVectorLayer layer,

            IVectorLayer overlay,

            int shapeType,

            boolean selectedGeomInput,

            boolean selectedGeomOverlay) throws DataException, GeoAlgorithmExecutionException {

        FeatureStore storeLayer = null;

        FeatureStore storeOverlay = null;

        if (layer instanceof FlyrVectIVectorLayer

                && overlay instanceof FlyrVectIVectorLayer) {

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

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

        } else {

            return false;

        }

        FeatureType featureType1 = storeLayer.getDefaultFeatureType();

        FeatureType featureType2 = storeOverlay.getDefaultFeatureType();

        FeatureStore outFeatStorePol = null;

        FeatureStore outFeatStoreLine = null;

        FeatureStore outFeatStorePoint = null;

        IntersectionOperation operation = new IntersectionOperation(storeOverlay, this);

        operation.setTaskStatus(getStatus());

        this.namesTranslator = NamesTranslator.createTrimTranslator(11);

        //La de puntos se genera siempre

        outFeatStorePoint

                = buildOutPutStoreFromUnion(featureType1, featureType2,

                        IVectorLayer.SHAPE_TYPE_POINT,

                        getTranslation("Intersection_point"), RESULT_POINT);

        if (outFeatStorePoint != null) {

            getStatus().setTitle("Point");

            operation.computesGeometryOperation(storeLayer, outFeatStorePoint,

                    attrNames, selectedGeomInput, selectedGeomOverlay, true);

        }

        //La de pol?gonos solo si es intersecci?n entre pol?gonos

        if (isPolygon(storeLayer) && isPolygon(storeOverlay)) {

            outFeatStorePol

                    = buildOutPutStoreFromUnion(featureType1, featureType2,

                            IVectorLayer.SHAPE_TYPE_MULTIPOLYGON, getTranslation("Intersection_polygon"), RESULT_POL);

            if (outFeatStorePol != null) {

                getStatus().setTitle("Polygon");

                operation.computesGeometryOperation(storeLayer,

                        outFeatStorePol, attrNames, selectedGeomInput,

                        selectedGeomOverlay, true);

            }

        }

        //La capa de l?neas se genera cuando ning?na de las dos es de puntos

        if (!isPoint(storeLayer) && !isPoint(storeOverlay)) {

            outFeatStoreLine

                    = buildOutPutStoreFromUnion(featureType1, featureType2,

                            IVectorLayer.SHAPE_TYPE_MULTILINE, getTranslation("Intersection_line"), RESULT_LINE);

            if (outFeatStoreLine != null) {

                getStatus().setTitle("Line");

                operation.computesGeometryOperation(storeLayer,

                        outFeatStoreLine, attrNames, selectedGeomInput,

                        selectedGeomOverlay, true);

            }

        }

        if (outFeatStorePol == null) {

            getNewVectorLayer(RESULT_POL, getTranslation("Null_polygon"),

                    OutputVectorLayer.SHAPE_TYPE_POLYGON, new Class[]{Integer.class}, new String[]{""});

        }

        if (outFeatStoreLine == null) {

            getNewVectorLayer(RESULT_LINE, getTranslation("Null_line"),

                    OutputVectorLayer.SHAPE_TYPE_LINE, new Class[]{Integer.class}, new String[]{""});

        }

        if (outFeatStorePoint == null) {

            getNewVectorLayer(RESULT_POINT, getTranslation("Null_point"),

                    OutputVectorLayer.SHAPE_TYPE_POINT, new Class[]{Integer.class}, new String[]{""});

        }

        return operation.getErrorInfo();

    }

}

/**

 * 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 java.util.Iterator;

import java.util.List;

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

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

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

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

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

import org.gvsig.fmap.geom.GeometryLocator;

import org.gvsig.fmap.geom.GeometryManager;

import org.gvsig.fmap.geom.Geometry.TYPES;

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

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.operation.fromjts.FromJTS;

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

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

import org.gvsig.geoprocess.algorithm.base.core.GeometryOperation;

import org.gvsig.geoprocess.algorithm.base.util.GeometryUtil;

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

import org.slf4j.Logger;

import org.slf4j.LoggerFactory;

import com.vividsolutions.jts.geom.Geometry;

import com.vividsolutions.jts.geom.GeometryCollection;

import com.vividsolutions.jts.geom.LineString;

import com.vividsolutions.jts.geom.MultiLineString;

import com.vividsolutions.jts.geom.MultiPoint;

import com.vividsolutions.jts.geom.MultiPolygon;

import com.vividsolutions.jts.geom.Point;

import com.vividsolutions.jts.geom.Polygon;

import com.vividsolutions.jts.precision.EnhancedPrecisionOp;

import es.unex.sextante.core.Sextante;

/**

 * Builds a geometry with the intersection between two layers

 *

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

 */

public class IntersectionOperation extends GeometryOperation {

    private FeatureStore storeOverlay = null;

    protected GeometryManager geomManager = GeometryLocator.getGeometryManager();

    private Logger log = LoggerFactory.getLogger(IntersectionOperation.class);

    private boolean errorInfo = false;

    public IntersectionOperation(FeatureStore overlayLayer, AbstractSextanteGeoProcess p) {

        super(p);

        this.storeOverlay = overlayLayer;

    }

    public boolean getErrorInfo() {

        return errorInfo;

    }

    /**

     * Computes intersection between the geometry and the overlay layer. The

     * fields of the intersected features will be added.

     *

     * @param g

     * @param featureInput

     * @return

     */

    @SuppressWarnings({"unchecked", "deprecation"})

    public EditableFeature invoke(org.gvsig.fmap.geom.Geometry g, Feature featureInput) {

        if (g == null) {

            return lastEditFeature;

        }

        FeatureSet features = null;

        Iterator it = null;

        try {

            if (selectedGeomOverlay) {

                features = (FeatureSet) storeOverlay.getSelection();

            } else {

                features = storeOverlay.getFeatureSet();

            }

            it = features.iterator();

        } catch (DataException e) {

            Sextante.addErrorToLog(e);

            return lastEditFeature;

        }

        while (it.hasNext()) {

            Feature featureOverlay = (Feature) it.next();

            List<org.gvsig.fmap.geom.Geometry> geomList = featureOverlay.getGeometries();

            if (geomList == null) {

                org.gvsig.fmap.geom.Geometry geom = featureOverlay.getDefaultGeometry();

                lastEditFeature = intersection(g, geom, featureInput, featureOverlay);

                continue;

            }

            Iterator<org.gvsig.fmap.geom.Geometry> itGeom = geomList.iterator();

            while (itGeom.hasNext()) {

                org.gvsig.fmap.geom.Geometry geom = itGeom.next();

                lastEditFeature = intersection(g, geom, featureInput, featureOverlay);

            }

        }

//		it.dispose();

        return lastEditFeature;

    }

    private EditableFeature intersection(org.gvsig.fmap.geom.Geometry g1,

            org.gvsig.fmap.geom.Geometry g2,

            Feature featureInput,

            Feature featureOverlay) {

        if (g1 == null || g2 == null) {

            return null;

        }

        Geometry overlaysGeom = null;

        Geometry jtsGeom = null;

        int outPutType = TYPES.SURFACE;

        try {

            outPutType = persister.getOutputFeatureStore().getDefaultFeatureType()

                    .getDefaultGeometryAttribute().getGeomType().getType();

        } catch (DataException e2) {

            Sextante.addErrorToLog(e2);

        }

        try {

            if (((g1 instanceof Surface && g2 instanceof Curve)

                    || (g2 instanceof Surface && g1 instanceof Curve))

                    && outPutType == TYPES.MULTIPOINT) {

                org.gvsig.fmap.geom.Geometry overGeom = (g1 instanceof Surface) ? g1 : g2;

                jtsGeom = (g1 instanceof Surface) ? GeometryUtil.geomToJTS(g2) : GeometryUtil.geomToJTS(g1);

                overlaysGeom = (Geometry) overGeom.invokeOperation("toJTSLineString", null);

            } else {

                if (g1 instanceof Surface

                        && g2 instanceof Surface

                        && (outPutType == TYPES.MULTIPOINT)) {

                    jtsGeom = (Geometry) g1.invokeOperation("toJTSLineString", null);

                    overlaysGeom = (Geometry) g2.invokeOperation("toJTSLineString", null);

                } else {

                    jtsGeom = GeometryUtil.geomToJTS(g1);

                    overlaysGeom = GeometryUtil.geomToJTS(g2);

                }

            }

            if (!jtsGeom.getEnvelope().intersects(overlaysGeom.getEnvelope())) {

                return lastEditFeature;

            }

            if (jtsGeom.intersects(overlaysGeom)) {

                Geometry newGeom = EnhancedPrecisionOp.intersection(jtsGeom, overlaysGeom);

                if (!newGeom.isEmpty()) {

                    if (typesMatch(outPutType, newGeom) || newGeom instanceof GeometryCollection) {

                        lastEditFeature = persister.addFeature(featureInput, featureOverlay, newGeom);

                    } else {

                        //Para intersecciones entre pol?gonos cuando la salida es de tipo l?nea

                        //la generamos a partir del pol?gono resultante de la intersecci?n

                        if (g1 instanceof Surface

                                && g2 instanceof Surface

                                && outPutType == TYPES.MULTICURVE

                                && (newGeom instanceof Polygon || newGeom instanceof MultiPolygon)) {

                            GeometryOperationContext ctx = new GeometryOperationContext();

                            ctx.setAttribute(FromJTS.PARAM, newGeom);

                            org.gvsig.fmap.geom.Geometry newDalGeom = (org.gvsig.fmap.geom.Geometry) geomManager.invokeOperation(FromJTS.NAME, ctx);

                            newGeom = (Geometry) newDalGeom.invokeOperation("toJTSLineString", null);

                        }

                        lastEditFeature = persister.addFeature(featureInput, featureOverlay, newGeom);

                    }

                }

            }

        } catch (CreateGeometryException e) {

            Sextante.addErrorToLog(e);

        } catch (DataException e) {

            Sextante.addErrorToLog(e);

        } catch (GeometryOperationNotSupportedException e1) {

            Sextante.addErrorToLog(e1);

        } catch (GeometryOperationException e1) {

            Sextante.addErrorToLog(e1);

        } catch (com.vividsolutions.jts.geom.TopologyException e) {

            errorInfo = true;

            log.info("Problems operating intersection: ", e);

        }

        return lastEditFeature;

    }

    private boolean typesMatch(int dalType, Geometry newGeom) {

        return ((geomManager.isSubtype(TYPES.MULTICURVE, dalType)

                && (newGeom instanceof MultiLineString || newGeom instanceof LineString))

                || (geomManager.isSubtype(TYPES.MULTIPOINT, dalType)

                && (newGeom instanceof MultiPoint || newGeom instanceof Point))

                || (geomManager.isSubtype(TYPES.MULTISURFACE, dalType)

                && (newGeom instanceof Polygon || newGeom instanceof MultiPolygon))

                || (geomManager.isSubtype(TYPES.CURVE, dalType) && newGeom instanceof LineString)

                || (geomManager.isSubtype(TYPES.SURFACE, dalType) && newGeom instanceof Polygon)

                || (geomManager.isSubtype(TYPES.POINT, dalType) && newGeom instanceof Point));

    }

    /**

     * clips feature's geometry with the clipping geometry, preserving feature's

     * original attributes. If feature's geometry doesn't touch clipping

     * geometry, it will be ignored.

     */

    public void invoke(org.gvsig.fmap.geom.Geometry g, EditableFeature featureInput) {

    }

}

<?xml version="1.0" encoding="UTF-8"?>

<project xmlns="http://maven.apache.org/POM/4.0.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://maven.apache.org/POM/4.0.0 http://maven.apache.org/maven-v4_0_0.xsd">

  <modelVersion>4.0.0</modelVersion>

  <artifactId>org.gvsig.geoprocess.algorithm.intersection</artifactId>

  <packaging>jar</packaging>

  <name>org.gvsig.geoprocess.algorithm.intersection</name>

	<parent>

		<groupId>org.gvsig</groupId>

		<artifactId>org.gvsig.geoprocess.algorithm</artifactId>

		<version>2.2.182</version>

	</parent>

	<dependencies>

		<dependency>

		    <groupId>org.gvsig</groupId>

   			<artifactId>org.gvsig.geoprocess.algorithm.base</artifactId>

            <scope>compile</scope>

   		</dependency>

	</dependencies>

</project>

/**

 * 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.buffer;

import org.cresques.cts.IProjection;

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

import org.gvsig.fmap.mapcontext.MapContext;

/**

 * Abstract distance implementation.

 * 

 * @author gvSIG Team

 * @version $Id$

 */

public abstract class AbstractDistance implements IDistance {

    public static final double EARTH_RADIUS = 6378137d;

    protected int numberOfRings = 1;

    protected Feature feature = null;

    public void setFeature(Feature feature) {

        this.feature = feature;

    }

    /**

     * For computing with geodetic coordinates:

     * returns the angular measure (in radians)

     * for a distance over the earth along a

     * meridiam.

     * Because this consideration is an approximation,

     * we consideer the eart like an sphere (not an

     * ellipsoid)

     * 

     * @param dist

     *            distance in meters

     * @return arc of meridian whose length is the specified

     *         distance

     */

    private double toSexaAngularMeasure(double dist) {

        /*

         * dist = 6378km(terrestrial radius) -> 2PI

         * passed distance -> incognite

         */

        return Math.toDegrees((2 * Math.PI * dist) / EARTH_RADIUS);

    }

    /**

     * Converts a distance entered by user in the GUI in the same distance

     * in internal units (measure units)

     */

    protected double getInInternalUnits(double userEntryDistance,

        IProjection proj, int distanceUnits, int mapUnits) {

       /* double[] trans2Meter = MapContext.getDistanceTrans2Meter();

        double distInInternalUnits =

            (userEntryDistance / trans2Meter[mapUnits])

                * trans2Meter[distanceUnits];

        if ((proj != null) && !(proj.isProjected()))

            distInInternalUnits = toSexaAngularMeasure(distInInternalUnits);

        return distInInternalUnits;*/

    	/*

    	 * En caso de que el sistema de referencia sea proyectado se usa la distancia

    	 * que el usuario ha definido, ya que esta se supone en metros. En caso de que

    	 * las coordenadas sean en geogr?ficas se toma la distancia que el usuario ha

    	 * definido en metros en el ecuador. Esto supone una desviaci?n a medida que 

    	 * nos alejamos del ecuador. Una posible soluci?n ser?a obtener esa distancia 

    	 * para la posici?n en el planeta en la que se encuentra la cartograf?a en la

    	 * proyecci?n en la que est? definida. 

    	 */

    	if ((proj != null) && !(proj.isProjected()))

    		return userEntryDistance / 111352D;

        return userEntryDistance;

    }

    public void setNumberOfRings(int n) {

    	numberOfRings = n;

    }

}

/**

 * 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.buffer;

import org.gvsig.geoprocess.algorithm.base.core.AlgorithmAbstractLibrary;

import org.gvsig.i18n.Messages;

import org.gvsig.tools.library.LibraryException;

/**

 * Initialization of BufferLibrary library.

 * 

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

 * @author gvSIG Team

 */

public class BufferLibrary extends AlgorithmAbstractLibrary {

    @Override

    protected void doInitialize() throws LibraryException {

        // Nothing to do

    }

    @Override

    protected void doPostInitialize() throws LibraryException {

        Messages.addResourceFamily(

            "org.gvsig.geoprocess.algorithm.buffer.buffer", BufferLibrary.class

                .getClassLoader(), BufferLibrary.class.getClass().getName());

        registerGeoProcess(new BufferAlgorithm());

    }

}

/**

 * 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.buffer;

import org.cresques.cts.IProjection;

/**

 * Computes a constant size of each geometry built

 * 

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

 */

public class FieldDistance extends AbstractDistance {

    private String attributeName = null;

    public FieldDistance(String attributeName) {

        this.attributeName = attributeName;

    }

    public double getBufferDistance(IProjection projection, int distanceUnits, int mapUnits) {

        double value = 0D;

        Object obj;

        obj = feature.get(attributeName);

        if (obj instanceof Number) {

            value = ((Number) obj).doubleValue();

        } else {

            return 0;

        }

        return getInInternalUnits(value, projection, distanceUnits, mapUnits) * numberOfRings;

    }

}

/**

 * 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.buffer;

import java.io.File;

import javax.swing.JOptionPane;

import org.cresques.cts.IProjection;

import org.gvsig.fmap.dal.DALLocator;

import org.gvsig.fmap.dal.DataManager;

import org.gvsig.fmap.dal.DataStoreParameters;

import org.gvsig.fmap.dal.DataTypes;

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

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

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

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

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

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

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

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

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

import org.gvsig.fmap.dal.serverexplorer.filesystem.FilesystemServerExplorer;

import org.gvsig.fmap.dal.serverexplorer.filesystem.FilesystemServerExplorerParameters;

import org.gvsig.fmap.geom.Geometry;

import org.gvsig.geoprocess.lib.api.GeoProcessLocator;

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

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

import es.unex.sextante.additionalInfo.AdditionalInfoNumericalValue;

import es.unex.sextante.core.Sextante;

import es.unex.sextante.dataObjects.IVectorLayer;

import es.unex.sextante.exceptions.GeoAlgorithmExecutionException;

import es.unex.sextante.exceptions.NullParameterValueException;

import es.unex.sextante.exceptions.OptionalParentParameterException;

import es.unex.sextante.exceptions.RepeatedParameterNameException;

import es.unex.sextante.exceptions.UndefinedParentParameterNameException;

import es.unex.sextante.exceptions.UnsupportedOutputChannelException;

import es.unex.sextante.exceptions.WrongParameterIDException;

import es.unex.sextante.exceptions.WrongParameterTypeException;

import es.unex.sextante.gui.algorithm.GeoAlgorithmParametersPanel;

import es.unex.sextante.outputs.OutputVectorLayer;

/**

 * Geoprocess that computes a buffer area around each feature's geometry of the

 * input layer. <br>

 * All the points interior to this buffer area has to be at a distance inferior

 * to "buffer distance" to the original geometry. This buffer distance could be

 * constant, or it could be a function of the value of a feature attribute.<br>

 * 

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

 * @author gvSIG Team

 */

public class BufferAlgorithm extends AbstractSextanteGeoProcess {

    public static final String RESULT                   = "RESULT";

    public static final String LAYER                    = "LAYER";

    public static final String SELECTED_GEOM            = "SELECTED_GEOM";

    public static final String DISTANCE                 = "DISTANCE";

    public static final String FIELD                    = "FIELD";

    public static final String DISSOLVE                 = "DISSOLVE";

    public static final String ROUND_BORDER             = "ROUND_BORDER";

    public static final String AREA                     = "AREA";

    public static final String RING_NUMBER              = "RING_NUMBER";

	/**

	 * For polygonal buffers, only compute exterior buffers

	 * (is the default operation, it applies to any geometry type)

	 */

	public static final byte BUFFER_OUTSIDE_POLY        = 0;

	/**

	 * For polygonal buffers, only compute interior buffers

	 */

	public static final byte BUFFER_INSIDE_POLY         = 1;

	/**

	 * For polygonal buffers, compute interior and exterior buffers

	 */

	public static final byte BUFFER_INSIDE_OUTSIDE_POLY = 2;

    public static final byte SOURCE_FIELDS              = 0;

    public static final byte FID_DIST_FIELDS            = 1;

    public static final byte FID_FROM_TO_FIELDS         = 2;

    protected byte           tableFields                = -1; 

    private int              inflArea                   = BUFFER_OUTSIDE_POLY;

    protected boolean        dissolve                   = false;

    protected int            rings                      = 0;

    protected boolean        selectedGeom               = false;

    protected FeatureStore   inputStore                 = null;

    protected FeatureStore   outputStore                = null;

    protected IVectorLayer   sextanteInputLayer         = null;

    protected boolean        round_border               = false; 

    public static String[] sOptions = {

    	GeoProcessLocator.getGeoProcessManager().getTranslation("poly_out"), 

    	GeoProcessLocator.getGeoProcessManager().getTranslation("poly_in"),

    	GeoProcessLocator.getGeoProcessManager().getTranslation("poly_inandout") 

    };

    public void defineCharacteristics() {

        setName(getTranslation("Buffer"));

        setGroup(getTranslation("basic_vect_algorithms"));

        // setGeneratesUserDefinedRasterOutput(false);

        try {

            m_Parameters.addInputVectorLayer(LAYER, getTranslation("Input_layer"), 

            	IVectorLayer.SHAPE_TYPE_WRONG, true);

            m_Parameters.addBoolean(SELECTED_GEOM, getTranslation("Selected_geometries"), false);

            m_Parameters.addNumericalValue(DISTANCE, getTranslation("area_distance"), 0,

                AdditionalInfoNumericalValue.NUMERICAL_VALUE_DOUBLE);

            m_Parameters.addTableField(FIELD, getTranslation("area_field"), "LAYER");

            m_Parameters.addBoolean(DISSOLVE, getTranslation("Dissolve_entities"), false);

            m_Parameters.addBoolean(ROUND_BORDER, getTranslation("Round_border"), true);

            m_Parameters.addSelection(AREA, getTranslation("Builds_influence_area"), sOptions);

            m_Parameters.addSelection(RING_NUMBER, getTranslation("Number_of_rings"),

                new String[] { "1", "2", "3" });

        } catch (RepeatedParameterNameException e) {

            Sextante.addErrorToLog(e);

        } catch (UndefinedParentParameterNameException e) {

            Sextante.addErrorToLog(e);

        } catch (OptionalParentParameterException e) {

            Sextante.addErrorToLog(e);

        }

        addOutputVectorLayer(RESULT, getTranslation("Buffer"),

            OutputVectorLayer.SHAPE_TYPE_POLYGON);

    }

    protected void readParameters() throws WrongParameterTypeException, NullParameterValueException, WrongParameterIDException {

    	sextanteInputLayer = m_Parameters.getParameterValueAsVectorLayer(LAYER);

    	selectedGeom = m_Parameters.getParameter(SELECTED_GEOM).getParameterValueAsBoolean();

    	rings = m_Parameters.getParameterValueAsInt(RING_NUMBER);

        dissolve = m_Parameters.getParameter(DISSOLVE).getParameterValueAsBoolean();

        round_border = m_Parameters.getParameter(ROUND_BORDER).getParameterValueAsBoolean();

    }

    public boolean processAlgorithm() throws GeoAlgorithmExecutionException {

    	if(existsOutPutFile(BufferAlgorithm.RESULT, 0)) {

    		throw new GeoAlgorithmExecutionException(getTranslation("file_exists"));

    	}

        double distanceValue = m_Parameters.getParameter(DISTANCE).getParameterValueAsDouble();

        String attributeName = m_Parameters.getParameterValueAsString(FIELD);

        inflArea = m_Parameters.getParameterValueAsInt(AREA);

        readParameters();

        if (sextanteInputLayer.getShapeType() != IVectorLayer.SHAPE_TYPE_POLYGON

            && inflArea != BUFFER_OUTSIDE_POLY) {

            JOptionPane.showMessageDialog(null,

                getTranslation("Wrong_type_for_this_shapetype"), "Error",

                JOptionPane.WARNING_MESSAGE);

            inflArea = BUFFER_OUTSIDE_POLY;

        }

        if (sextanteInputLayer instanceof FlyrVectIVectorLayer)

            inputStore = ((FlyrVectIVectorLayer) sextanteInputLayer).getFeatureStore();

        else

            return false;

        try {

            // Object to compute the distance

            IDistance distance = null;

            if (distanceValue == 0) {

                //int featureAttributePosition = ((FlyrVectIVectorLayer) sextanteInputLayer).getFeatureIndexByFieldIndex(attributePosition);

                distance = new FieldDistance(attributeName);//featureAttributePosition);

            } else {

                distance = new ConstantDistance(distanceValue);

            }

            // Object to compute the buffer operation

            BufferOperation operation = null;

            switch (inflArea) {

            case BUFFER_OUTSIDE_POLY:

                operation = new OutBufferOperation(distance, inputStore, this, getTableFieldsStructure());

                break;

            case BUFFER_INSIDE_POLY:

                operation = new InBufferOperation(distance, inputStore, this, getTableFieldsStructure());

                break;

            case BUFFER_INSIDE_OUTSIDE_POLY:

                operation = new InOutBufferOperation(distance, inputStore, this, getTableFieldsStructure());

                break;

            }

            operation.setTypeOfCap(round_border ? BufferOperation.CAP_ROUND : BufferOperation.CAP_SQUARE);

            operation.setGeoProcess(this, 100);

        	// Builds the output FeatureStore

            outputStore = buildOutPutStore(IVectorLayer.SHAPE_TYPE_POLYGON,

        					getTranslation("Buffer"), RESULT);

        	if(!dissolve) {

        		computesBufferAlgWithoutDissolve(operation);

        	} else {

        		computesBufferAlgWithDissolve(operation);

        	}

        } catch (DataException e) {

            Sextante.addErrorToLog(e);

            return false;

        }

		if(getTaskMonitor().isCanceled())

			return false;

        return true;

    }

    private void computesOneRingBufferWithDissolve(

    		BufferOperation operation, 

    		FuseOperation fuseOp, 

    		String idFile,

    		int initialPosition,

    		int ringPosition) throws ValidateDataParametersException, DataException {

    	String file = System.getProperty("java.io.tmpdir") + File.separator + (idFile) + ".shp";

		FeatureStore outAuxFeatStore = buildTemporalStore(org.gvsig.fmap.geom.Geometry.TYPES.SURFACE, file, 

				inputStore.getDefaultFeatureType().getDefaultSRS());

		operation.getDistance().setNumberOfRings(ringPosition + 1);

		operation.setNumberOfRadialBuffers(1);

		operation.setTaskStatus(getStatus());

		m_Task.setProgressText(getTranslation("calc_buffer") + " [Ring:" + (rings - ringPosition) + "]");

		operation.computesGeometryOperation(inputStore, outAuxFeatStore,

				attrNames, selectedGeom, false, true);

		outAuxFeatStore = open(file, inputStore.getDefaultFeatureType().getDefaultSRS());

		m_Task.setProgressText(getTranslation("fuse_spatially") + " [Ring:" + (rings - ringPosition) + "]");

		fuseOp.computesGeometryOperation(outAuxFeatStore, outputStore, attrNames, ringPosition == initialPosition ? true : false);

		outAuxFeatStore.dispose();

    }

    protected void computesBufferAlgWithDissolve(BufferOperation operation) throws DataException {

    	long fileTmpId = System.currentTimeMillis();

    	FuseOperation fuseOp = new FuseOperation(this);

    	if(inflArea == BUFFER_OUTSIDE_POLY || inflArea == BUFFER_INSIDE_OUTSIDE_POLY) {

    		try {

    			for (int i = rings; i >= 0 ; i--) {

    				computesOneRingBufferWithDissolve(

    						operation, 

    						fuseOp, 

    						fileTmpId + i + "",

    						rings,

    						i);

    			}

    		} catch (ValidateDataParametersException e) {

    			Sextante.addErrorToLog(e);

    		}

    	}

    	if(inflArea == BUFFER_INSIDE_POLY) {

    		try {

    			for (int i = 0; i < rings + 1 ; i++) {

    				computesOneRingBufferWithDissolve(

    						operation, 

    						fuseOp, 

    						fileTmpId + i + "",

    						0,

    						i);

    			}

    		} catch (ValidateDataParametersException e) {

    			Sextante.addErrorToLog(e);

    		}

    	}

    	fuseOp.end();

    }

    protected void computesBufferAlgWithoutDissolve(BufferOperation operation) throws DataException {

    	operation.setNumberOfRadialBuffers(rings + 1);

    	operation.setTaskStatus(getStatus());

    	operation.computesGeometryOperation(inputStore, 

    			outputStore,

    			attrNames, 

    			selectedGeom, 

    			false, 

    			true);

    }

    /**

     * Gets the constant with the table structure 

     * @return

     */

    protected byte getTableFieldsStructure() {

    	if(tableFields == -1) {

            if(!dissolve) {

            	tableFields = SOURCE_FIELDS;

        	} else {

                if (inflArea == BUFFER_INSIDE_OUTSIDE_POLY) 

                	tableFields = FID_FROM_TO_FIELDS;

                else

                	tableFields = FID_DIST_FIELDS;

        	}

    	}

    	return tableFields;

    }

    /**

     * Open a <code>DataStore</code>

     * @param file

     * @param proj

     * @return

     * @throws InitializeException

     * @throws ProviderNotRegisteredException

     * @throws ValidateDataParametersException

     */

    private FeatureStore open(String file, IProjection proj) throws InitializeException, ProviderNotRegisteredException, ValidateDataParametersException {

    	DataManager manager = DALLocator.getDataManager();

    	DataStoreParameters params = manager.createStoreParameters("Shape");

    	params.setDynValue("shpfile", file);

    	params.setDynValue("crs", proj);

    	FeatureStore featureStore = (FeatureStore) manager.openStore(params.getDataStoreName(), params);

    	return featureStore;

    }

    /**

     * Builds the output FeatureStore

     * 

     * @param featureType

     * @return FeatureStore

     * @throws DataException 

     */

    protected FeatureStore buildOutPutStore(int shapeType, 

    		String sextanteLayerName, 

    		String sextanteLayerLabel) throws DataException {

    	FeatureType featureType = inputStore.getDefaultFeatureType();

    	if(getTableFieldsStructure() == SOURCE_FIELDS) {

    		return super.buildOutPutStore(featureType, shapeType,sextanteLayerName, sextanteLayerLabel);

    	} 

    	Class<?>[] types = null;

    	if (getTableFieldsStructure() == FID_FROM_TO_FIELDS) {

    		types = new Class[] { Integer.class, Double.class, Double.class };

    		attrNames = new String[] { "FID", "FROM", "TO" };

    	} 

    	if (getTableFieldsStructure() == FID_DIST_FIELDS) {

    		types = new Class[] { Integer.class, Double.class };

    		attrNames = new String[] { "FID", "DIST" };

    	}

    	try {

    		IVectorLayer output =

    				getNewVectorLayer(sextanteLayerLabel, sextanteLayerName,

    						shapeType, types, attrNames);

    		return ((FlyrVectIVectorLayer) output).getFeatureStore();

    	} catch (UnsupportedOutputChannelException e) {

    		Sextante.addErrorToLog(e);

    	} catch (GeoAlgorithmExecutionException e) {

    		Sextante.addErrorToLog(e);

    	}

    	return null;

    }

    /**

     * Builds the output FeatureStore

     * 

     * @param featureType

     * @return FeatureStore

     * @throws DataException 

     * @throws ValidateDataParametersException 

     */

    protected FeatureStore buildTemporalStore(int shapeType, 

    		String file, 

    		IProjection crs) throws DataException, ValidateDataParametersException {

    	FeatureType featureType = inputStore.getDefaultFeatureType();

        int[] types = null;

        if(getTableFieldsStructure() == SOURCE_FIELDS) {

        	FeatureAttributeDescriptor[] desc = featureType.getAttributeDescriptors();

        	types = new int[desc.length - 1];

        	attrNames = new String[desc.length - 1];

        	int attrCount = 0;

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

        		if(desc[i].getType() != DataTypes.GEOMETRY) {

        			types[attrCount] = desc[i].getType();

        			attrNames[attrCount] = desc[i].getName();

        			attrCount ++;

        		}

			}

        }

        if (getTableFieldsStructure() == FID_FROM_TO_FIELDS) {

            types = new int[] { java.sql.Types.INTEGER, java.sql.Types.DOUBLE, java.sql.Types.DOUBLE };

            attrNames = new String[] { "FID", "FROM", "TO" };

        }

        if (getTableFieldsStructure() == FID_DIST_FIELDS) {

            types = new int[] { java.sql.Types.INTEGER, java.sql.Types.DOUBLE };

            attrNames = new String[] { "FID", "DIST" };

        }

        return create(file, types, attrNames, crs);

    }

	@SuppressWarnings("deprecation")

	public FeatureStore create(String sFilename,

			int[] types, String[] attrNames, IProjection crs) {

		try {

			DataManager manager = DALLocator.getDataManager();

			FilesystemServerExplorerParameters explorerParams =

			    (FilesystemServerExplorerParameters) manager.createServerExplorerParameters("FilesystemExplorer");

			explorerParams.setRoot(new File(sFilename).getParent());

			FilesystemServerExplorer explorer = (FilesystemServerExplorer) manager.createServerExplorer(explorerParams);

			NewFeatureStoreParameters newParams = (NewFeatureStoreParameters) explorer.getAddParameters(new File(sFilename));

			EditableFeatureType ft = newParams.getDefaultFeatureType();

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

				ft.add(attrNames[i], types[i], 7);

			}

			ft.add("GEOMETRY", DataTypes.GEOMETRY).setGeometryType(

					Geometry.TYPES.SURFACE).setGeometrySubType(Geometry.SUBTYPES.GEOM2D);

			newParams.setDefaultFeatureType(ft);

			newParams.setDynValue("CRS", crs);

			newParams.setDynValue("geometryType", Geometry.TYPES.SURFACE);

			manager.newStore("FilesystemExplorer",

                newParams.getDataStoreName(), newParams, true);

			FeatureStore featureStore = (FeatureStore) manager.openStore(newParams.getDataStoreName(), newParams);			

			featureStore.edit(FeatureStore.MODE_APPEND);

			return featureStore;

		} catch (Exception e) {

			Sextante.addErrorToLog(e);

		}

		return null;

	}

    @Override

    public Class<? extends GeoAlgorithmParametersPanel> getCustomParametersPanelClass() {

        return BufferParametersPanel.class;

    }

}

/**

 * 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.buffer;

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

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

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

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

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

import org.gvsig.geoprocess.algorithm.base.util.GeometryUtil;

import org.gvsig.geoprocess.algorithm.base.util.JTSFacade;

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

import com.vividsolutions.jts.geom.Geometry;

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

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

import com.vividsolutions.jts.simplify.TopologyPreservingSimplifier;

import es.unex.sextante.core.Sextante;

/**

 * Buffer operation

 * 

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

 * @author gvSIG Team

 */

public class InBufferOperation extends BufferOperation {