Application: gvSIG desktop » gvSIG geoprocess

#/bin/bash

#USER=

#PASS=

REPO="https://devel.gvsig.org/svn/gvsig-geoprocess"

INCREMENT=3 #1:major number 2:minor number 3:bug number

PARAMETERS_LIN="-Dnative-platform=linux -Dnative-compiler=gcc4 -Dnative-arch=i386 -Dnative-distribution=all -Dlibrary-type=dynamic -Dnative-classifier=linux-all-gcc4-i386-dynamic -Dnative_classifier=linux-all-gcc4-i386-dynamic -Dmaven.test.skip"

VERSION_TO_FIX=`cat pom.xml | grep "<version>" | head -n 1 | sed 's/<version>//g' | sed 's/-SNAPSHOT<\/version>//g'`;

MAJOR_NUMBER=`echo $VERSION_TO_FIX | cut -d"." -f 1`;

MINOR_NUMBER=`echo $VERSION_TO_FIX | cut -d"." -f 2`;

BUG_NUMBER=`echo $VERSION_TO_FIX | cut -d"." -f 3`;

VERSION_TO_FIX=$MAJOR_NUMBER"."$MINOR_NUMBER"."$BUG_NUMBER

if [ $INCREMENT -eq 1 ]; then

  MAJOR_NUMBER=$((MAJOR_NUMBER+1));

fi

if [ $INCREMENT -eq 2 ]; then

  MINOR_NUMBER=$((MINOR_NUMBER+1));

fi

if [ $INCREMENT -eq 3 ]; then

  BUG_NUMBER=$((BUG_NUMBER+1));

fi

NEXT_DEV_VERSION=$MAJOR_NUMBER"."$MINOR_NUMBER"."$BUG_NUMBER"-SNAPSHOT"

echo Fixed version:$VERSION_TO_FIX

echo New development version:$NEXT_DEV_VERSION

#VERSION_TO_FIX=2.2.2

#NEXT_DEV_VERSION=2.2.3-SNAPSHOT

#mvn release:prepare -B -DdryRun=false -Dtag=v$VERSION_TO_FIX -DreleaseVersion=$VERSION_TO_FIX -DdevelopmentVersion=$NEXT_DEV_VERSION -DpreparationGoals="clean install" -Dusername=$USER -Dpassword=$PASS -DignoreSnapshots=true -Darguments=$PARAMETERS_LIN

read -p "¿Prepare de $i correcto? ¿Continuar? s/n: " CONTINUE

if [ "$CONTINUE" == "n" ]; then

  mvn release:rollback

  exit;

fi

mvn release:perform -Dusername=$USER -Dpassword=$PASS -Darguments=$PARAMETERS_LIN

read -p "¿Perform de $i correcto? ¿Continuar? s/n: " CONTINUE

if [ "$CONTINUE" == "n" ]; then

  mvn release:rollback

  exit;

fi

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 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.gvsig.tools.dispose.DisposableIterator;

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;

		DisposableIterator 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 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) {

		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) {

			/*try {

				Shape shape1 = g1.getShape();

				Area area1 = new Area(shape1);

				Shape shape2 = g2.getShape();

				Area area2 = new Area(shape2);

				area1.intersect(area2);

				PathIterator it = area1.getPathIterator(null);

				Surface s = null;

				MultiSurface ms = (MultiSurface)geomManager.create(TYPES.MULTISURFACE, SUBTYPES.GEOM2D);

				float coords[] = new float[6];

				while(!it.isDone()) {

					switch (it.currentSegment(coords)) {

					case PathIterator.SEG_MOVETO:

						s = (Surface)geomManager.create(TYPES.SURFACE, SUBTYPES.GEOM2D);

					case PathIterator.SEG_LINETO:

						org.gvsig.fmap.geom.primitive.Point p = (org.gvsig.fmap.geom.primitive.Point)geomManager.create(TYPES.POINT, SUBTYPES.GEOM2D);

						p.setX(coords[0]);

						p.setY(coords[1]);

						s.addVertex(p);

						break;

					case PathIterator.SEG_CLOSE:

						p = (org.gvsig.fmap.geom.primitive.Point)geomManager.create(TYPES.POINT, SUBTYPES.GEOM2D);

						p.setX(coords[0]);

						p.setY(coords[1]);

						s.addVertex(p);

						ms.addSurface(s);

						break;

					}

					it.next();

				}

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

			} catch (CreateGeometryException e1) {

			} catch (DataException e1) {

			}*/

			errorInfo = true;

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

		}

		return lastEditFeature;

	}

	private boolean typesMatch(int dalType, Geometry newGeom) {

		return 

		   ((dalType == TYPES.MULTICURVE && 

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

		   (dalType == TYPES.MULTIPOINT && 

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

		   (dalType == TYPES.MULTISURFACE && 

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

	       (dalType == TYPES.CURVE && newGeom instanceof LineString) ||

	       (dalType == TYPES.SURFACE && newGeom instanceof Polygon) ||

	       (dalType == TYPES.POINT && newGeom instanceof Point));

	}

//	private MultiCurve convertPolygonJTS2LineJTS(Polygon p) {

//		

//	}

	/**

	 * 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) {

		/*if(g == null)

			return;

		Geometry jtsGeom = Converter.geometryToJts(g);

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

			return;

		if(jtsGeom.intersects(overlaysGeom)) {

			try {

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

				persister.addFeature(feature, newGeom);

			} catch (CreateGeometryException e) {

				Sextante.addErrorToLog(e);

			} catch (DataException e) {

				Sextante.addErrorToLog(e);

			}

		} */

	}

}

/**

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

/**

 * 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());

		//La de puntos se genera siempre

		outFeatStorePoint =

			buildOutPutStoreFromUnion(featureType1, featureType2,

					IVectorLayer.SHAPE_TYPE_POINT,

					getTranslation("Intersection_point"), RESULT_POINT);

		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_POLYGON, getTranslation("Intersection_polygon"), RESULT_POL);

			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_LINE, getTranslation("Intersection_line"), RESULT_LINE);

			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();

	}

}

<?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.3</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 int attributePosition = 0;

    public FieldDistance(int attributePosition) {

        this.attributePosition = attributePosition;

    }

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

        double attr = 0D;

        Object obj = feature.get(attributePosition);

        if (obj instanceof Double)

            attr = ((Double) obj).doubleValue();

        else

            if (obj instanceof Integer)

                attr = ((Integer) obj).doubleValue();

            else

                if (obj instanceof Float)

                    attr = ((Float) obj).doubleValue();

                else

                    return 0;

        return getInInternalUnits(attr, 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();

        int attributePosition = m_Parameters.getParameterValueAsInt(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)

                distance = new FieldDistance(attributePosition);

            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 ++;

        		}

			}

        }