Application: gvSIG desktop » gvSIG geoprocess

<?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>1.0.0-SNAPSHOT</version>

	</parent>

	<dependencies>

		<dependency>

		    <groupId>org.gvsig</groupId>

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

   			<version>1.0.0-SNAPSHOT</version>

            <scope>compile</scope>

   		</dependency>

	</dependencies>

</project>

<assembly>

<!-- Empty distribution file -->

</assembly>

#maven.buildNumber.plugin properties file

#Thu Aug 18 10:21:29 CEST 2011

buildNumber=2042

gvSIG_Algorithms=Algoritmos gvSIG

Input_layer=Capa de entrada

Overlays_layer=Capa de revestimiento

Selected_geometries=Geometrias seleccionadas

Intersection=Intersecci?n

gvSIG_Algorithms=gvSIG algorithms

Input_layer=Input cover

Overlays_layer=Overlay layer

Selected_geometries=Selected geometries

Intersection=Intersection

org.gvsig.geoprocess.algorithm.intersection.IntersectionLibrary

/*

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

 *

 * Copyright (C) 2010 Generalitat Valenciana.

 *

 * This program is free software; you can redistribute it and/or

 * modify it under the terms of the GNU General Public License

 * as published by the Free Software Foundation; either version 2

 * of the License, or (at your option) any later version.

 *

 * This program is distributed in the hope that it will be useful,

 * but WITHOUT ANY WARRANTY; without even the implied warranty of

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 * GNU General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 *

 */

package org.gvsig.geoprocess.algorithm.intersection;

import es.unex.sextante.core.Sextante;

import es.unex.sextante.dataObjects.IVectorLayer;

import es.unex.sextante.exceptions.GeoAlgorithmExecutionException;

import es.unex.sextante.exceptions.RepeatedParameterNameException;

import es.unex.sextante.outputs.OutputVectorLayer;

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

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

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

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

import org.gvsig.geoprocess.core.gvVectorLayer;

import org.gvsig.geoprocess.sextante.AbstractSextanteGeoProcess;

/**

 * Intersection algorithm

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

 */

public class IntersectionAlgorithm extends AbstractSextanteGeoProcess {

	public static final String  RESULT    = "RESULT";

	public static final String  LAYER     = "LAYER";

	public static final String  INTER     = "INTER";

	public static final String  CHECK     = "CHECK";

	/*

	 * (non-Javadoc)

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

	 */

	public void defineCharacteristics() {

		setName(Sextante.getText("Intersection"));

		setGroup(Sextante.getText("gvSIG_Algorithms"));

        // setGeneratesUserDefinedRasterOutput(false);

		try {

			m_Parameters.addInputVectorLayer(LAYER, 

												Sextante.getText("Input_layer"), 

												IVectorLayer.SHAPE_TYPE_WRONG, 

												true);

			m_Parameters.addInputVectorLayer(INTER, 

												Sextante.getText("Overlays_layer"), 

												IVectorLayer.SHAPE_TYPE_WRONG, 

												true);

			m_Parameters.addBoolean(CHECK, Sextante.getText("Selected_geometries"), false);

		} catch (RepeatedParameterNameException e) {

			Sextante.addErrorToLog(e);

		}

		addOutputVectorLayer(RESULT,

								Sextante.getText("Intersection"),

								OutputVectorLayer.SHAPE_TYPE_UNDEFINED);

	}

	/*

	 * (non-Javadoc)

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

	 */

	public boolean processAlgorithm() throws GeoAlgorithmExecutionException {

		IVectorLayer inter = m_Parameters.getParameterValueAsVectorLayer(INTER);

		IVectorLayer layer = m_Parameters.getParameterValueAsVectorLayer(LAYER);

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

		try {

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

		} catch (DataException e) {

			Sextante.addErrorToLog(e);

			return false;

		}

		return true;

	}

	/**

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

	 * @param layer

	 *        Input layer

	 * @param templateGeometry

	 * @param shapeType

	 *        Output shape type

	 * @param selectedGeom

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

	 */

	private void computesIntersection(IVectorLayer layer,

								IVectorLayer overlay,

								int shapeType, 

								boolean selectedGeom) throws DataException {

		FeatureStore storeLayer = null;

		FeatureStore storeOverlay = null;

		if(layer instanceof gvVectorLayer && 

			overlay instanceof gvVectorLayer) {

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

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

		} else

			return;

		FeatureSet features1 = null;

		features1 = storeLayer.getFeatureSet();

		FeatureType featureType1 = features1.getDefaultFeatureType();

		FeatureSet features2 = null;

		features2 = storeOverlay.getFeatureSet();

		FeatureType featureType2 = features2.getDefaultFeatureType();

		FeatureStore outFeatStore = null;

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

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

			outFeatStore = buildOutPutStoreFromUnion(featureType1, featureType2, IVectorLayer.SHAPE_TYPE_POINT, 

													Sextante.getText("Intersection"), RESULT);

		else

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

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

				outFeatStore = buildOutPutStoreFromUnion(featureType1, featureType2, IVectorLayer.SHAPE_TYPE_POINT, 

														Sextante.getText("Intersection"), RESULT);

			else

				outFeatStore = buildOutPutStoreFromUnion(featureType1, featureType2, shapeType, 

														Sextante.getText("Intersection"), RESULT);

		IntersectionOperation operation = new IntersectionOperation(storeOverlay);

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

	}

}

/*

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

 *

 * Copyright (C) 2010 Generalitat Valenciana.

 *

 * This program is free software; you can redistribute it and/or

 * modify it under the terms of the GNU General Public License

 * as published by the Free Software Foundation; either version 2

 * of the License, or (at your option) any later version.

 *

 * This program is distributed in the hope that it will be useful,

 * but WITHOUT ANY WARRANTY; without even the implied warranty of

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 * GNU General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 */

package org.gvsig.geoprocess.algorithm.intersection;

import java.util.Iterator;

import java.util.List;

import com.vividsolutions.jts.geom.Geometry;

import com.vividsolutions.jts.precision.EnhancedPrecisionOp;

import es.unex.sextante.core.Sextante;

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.exception.CreateGeometryException;

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

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

import org.gvsig.tools.dispose.DisposableIterator;

/**

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

	public IntersectionOperation(FeatureStore overlayLayer) {

		this.storeOverlay = overlayLayer;

	}

	/**

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

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

		if(g == null)

			return lastEditFeature;

		Geometry jtsGeom = GeometryUtil.geomToJTS(g);

		FeatureSet features = null;

		DisposableIterator it = null;

		try {

			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(jtsGeom, 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(jtsGeom, geom, featureInput, featureOverlay);

			}

		}

		it.dispose();

		return lastEditFeature;

	}

	private EditableFeature intersection(Geometry jtsGeom, 

											org.gvsig.fmap.geom.Geometry overGeom, 

											Feature featureInput, 

											Feature featureOverlay) {

		Geometry overlaysGeom = GeometryUtil.geomToJTS(overGeom);

		try {

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

				return lastEditFeature;

			if(jtsGeom.intersects(overlaysGeom)) {

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

				if(!newGeom.isEmpty())

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

			}

		} catch (CreateGeometryException e) {

			Sextante.addErrorToLog(e);

		} catch (DataException e) {

			Sextante.addErrorToLog(e);

		}

		return lastEditFeature;

	}

	/**

	 * 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. Sistema de Informaci�n Geogr�fica de la Generalitat Valenciana

 *

 * Copyright (C) 2010 Generalitat Valenciana.

 *

 * This program is free software; you can redistribute it and/or

 * modify it under the terms of the GNU General Public License

 * as published by the Free Software Foundation; either version 2

 * of the License, or (at your option) any later version.

 *

 * This program is distributed in the hope that it will be useful,

 * but WITHOUT ANY WARRANTY; without even the implied warranty of

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 * GNU General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 *

 */

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(

            IntersectionAlgorithm.class,

            getLanguageStrings("org.gvsig.geoprocess.algorithm.intersection.intersection"));

    }

}

#maven.buildNumber.plugin properties file

#Thu Aug 18 10:21:15 CEST 2011

buildNumber=2042

/*

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

 *

 * Copyright (C) 2010 Generalitat Valenciana.

 *

 * This program is free software; you can redistribute it and/or

 * modify it under the terms of the GNU General Public License

 * as published by the Free Software Foundation; either version 2

 * of the License, or (at your option) any later version.

 *

 * This program is distributed in the hope that it will be useful,

 * but WITHOUT ANY WARRANTY; without even the implied warranty of

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 * GNU General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 *

 */

package org.gvsig.geoprocess.algorithm.clip;

import es.unex.sextante.core.Sextante;

import es.unex.sextante.dataObjects.IVectorLayer;

import es.unex.sextante.exceptions.GeoAlgorithmExecutionException;

import es.unex.sextante.exceptions.RepeatedParameterNameException;

import es.unex.sextante.outputs.OutputVectorLayer;

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

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

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

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

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

import org.gvsig.geoprocess.core.gvVectorLayer;

import org.gvsig.geoprocess.sextante.AbstractSextanteGeoProcess;

/**

 * Clip algorithm

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

 */

public class ClipAlgorithm extends AbstractSextanteGeoProcess {

	public static final String  RESULT    = "RESULT";

	public static final String  LAYER     = "LAYER";

	public static final String  CLIP      = "CLIP";

	public static final String  CHECK     = "CHECK";

	/*

	 * (non-Javadoc)

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

	 */

	public void defineCharacteristics(){

		setName(Sextante.getText("Clip"));

		setGroup(Sextante.getText("gvSIG_Algorithms"));

        // setGeneratesUserDefinedRasterOutput(false);

		try {

			m_Parameters.addInputVectorLayer(LAYER, 

												Sextante.getText("Input_layer"), 

												IVectorLayer.SHAPE_TYPE_WRONG, 

												true);

			m_Parameters.addInputVectorLayer(CLIP, 

												Sextante.getText("Clip_layer"), 

												IVectorLayer.SHAPE_TYPE_WRONG, 

												true);

			m_Parameters.addBoolean(CHECK, Sextante.getText("Selected_geometries"), false);

		} catch (RepeatedParameterNameException e) {

			Sextante.addErrorToLog(e);

		}

		addOutputVectorLayer(RESULT,

								Sextante.getText( "Clip"),

								OutputVectorLayer.SHAPE_TYPE_UNDEFINED);

	}

	/*

	 * (non-Javadoc)

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

	 */

	public boolean processAlgorithm() throws GeoAlgorithmExecutionException {

		org.gvsig.fmap.geom.Geometry clippingGeometry = null;

		IVectorLayer clip = m_Parameters.getParameterValueAsVectorLayer(CLIP);

		IVectorLayer layer = m_Parameters.getParameterValueAsVectorLayer(LAYER);

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

		try {

			clippingGeometry = ScalableUnionOperation.joinLayerGeometries(clip);

		} catch (Exception e) {

			Sextante.addErrorToLog(e);

			return false;

		}

		FeatureStore storeLayer = null;

		if(layer instanceof gvVectorLayer && clippingGeometry != null)

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

		else

			return false;

		try {

			FeatureSet features = null;

			features = storeLayer.getFeatureSet();

			FeatureType featureType = features.getDefaultFeatureType();

			FeatureStore outFeatStore = buildOutPutStore(featureType, layer.getShapeType(), Sextante.getText("Clip"), RESULT);

			ClipOperation operation = new ClipOperation(clippingGeometry);

			operation.setProgressModel(this);

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

		} catch (DataException e) {

			Sextante.addErrorToLog(e);

			return false;

		}

		return true;

	}

}

/*

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

 *

 * Copyright (C) 2010 Generalitat Valenciana.

 *

 * This program is free software; you can redistribute it and/or

 * modify it under the terms of the GNU General Public License

 * as published by the Free Software Foundation; either version 2

 * of the License, or (at your option) any later version.

 *

 * This program is distributed in the hope that it will be useful,

 * but WITHOUT ANY WARRANTY; without even the implied warranty of

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 * GNU General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 */

package org.gvsig.geoprocess.algorithm.clip;

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.geom.exception.CreateGeometryException;

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

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

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

import com.vividsolutions.jts.geom.Geometry;

import es.unex.sextante.core.Sextante;

/**

 * This class analyzes all features of a layer, clipping its geometries

 * with the convex hull of another layer.

 * If the geometry of the feature analyzed doesnt intersect with the convex

 * hull geometry, the clipvisitor will ignore it.

 *

 * It intersects, computes intersection and creates a new feature with

 * the same attributes and the new intersection geometry.

 * 

 * @author Nacho Brodin (nachobrodin@gmail.com)

 */

public class ClipOperation extends GeometryOperation {

	/**

	 * Clipping geometry: the convex hull of the

	 * clipping layer

	 */

	private Geometry                         clippingConvexHull   = null;

	private DALFeaturePersister              persister            = null;

	public ClipOperation(org.gvsig.fmap.geom.Geometry clip) {

		this.clippingConvexHull = GeometryUtil.geomToJTS(clip);

	}

	/**

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

	 * feature's original attributes.

	 * If feature's geometry doesnt touch clipping geometry, it will be

	 * ignored.

	 */

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

		if(g == null)

			return lastEditFeature;

		Geometry jtsGeom = GeometryUtil.geomToJTS(g);

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

			return lastEditFeature;

		if(jtsGeom.intersects(clippingConvexHull)) {

			try {

				Geometry newGeom = jtsGeom.intersection(clippingConvexHull);

				lastEditFeature = persister.addFeature(feature, newGeom);

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

				Sextante.addErrorToLog(e);

				if(! jtsGeom.isValid()) {

					System.out.println("La geometria de entrada no es valida");

					jtsGeom = GeometryUtil.removeDuplicatesFrom(jtsGeom);

				}

				if(! clippingConvexHull.isValid()) {

					System.out.println("La geometria de recorte no es valida");

					clippingConvexHull = GeometryUtil.removeDuplicatesFrom(clippingConvexHull);

				}

				try {

					Geometry newGeom = jtsGeom.intersection(clippingConvexHull);

					lastEditFeature = persister.addFeature(feature, newGeom);

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

					Sextante.addErrorToLog(ee);

				} catch (CreateGeometryException ee) {

					Sextante.addErrorToLog(ee);

				} catch (DataException ee) {

					Sextante.addErrorToLog(ee);

				}

			} catch (CreateGeometryException e) {

				Sextante.addErrorToLog(e);

			} catch (DataException e) {

				Sextante.addErrorToLog(e);

			}

		}

		return lastEditFeature;

	}

	/**

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

	 * feature's original attributes.

	 * If feature's geometry doesnt touch clipping geometry, it will be

	 * ignored.

	 */

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

		if(g == null)

			return;

		lastEditFeature = feature;

		Geometry jtsGeom = GeometryUtil.geomToJTS(g);

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

			return;

		if(jtsGeom.intersects(clippingConvexHull)) {

			try {

				Geometry newGeom = jtsGeom.intersection(clippingConvexHull);

				persister.addFeature(feature, newGeom);

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

				Sextante.addErrorToLog(e);

				if(! jtsGeom.isValid()) {

					System.out.println("La geometria de entrada no es valida");

					jtsGeom = GeometryUtil.removeDuplicatesFrom(jtsGeom);

				}

				if(! clippingConvexHull.isValid()) {

					System.out.println("La geometria de recorte no es valida");

					clippingConvexHull = GeometryUtil.removeDuplicatesFrom(clippingConvexHull);

				}

				try {

					Geometry newGeom = jtsGeom.intersection(clippingConvexHull);

					persister.addFeature(feature, newGeom);

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

					Sextante.addErrorToLog(ee);

				} catch (CreateGeometryException ee) {

					Sextante.addErrorToLog(ee);

				} catch (DataException ee) {

					Sextante.addErrorToLog(ee);

				}

			} catch (CreateGeometryException e) {

				Sextante.addErrorToLog(e);

			} catch (DataException e) {

				Sextante.addErrorToLog(e);

			}

		}

	}

	/**

	 * Ends the edition and closes the FeatureStore

	 */

	public void end() {

		persister.end();

	}

	public String getProcessDescription() {

		return "Clipping features agaisnt a clip geometry";

	}

}

/*

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

 *

 * Copyright (C) 2010 Generalitat Valenciana.

 *

 * This program is free software; you can redistribute it and/or

 * modify it under the terms of the GNU General Public License

 * as published by the Free Software Foundation; either version 2

 * of the License, or (at your option) any later version.

 *

 * This program is distributed in the hope that it will be useful,

 * but WITHOUT ANY WARRANTY; without even the implied warranty of

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 * GNU General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 *

 */

package org.gvsig.geoprocess.algorithm.clip;

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

import org.gvsig.i18n.Messages;

import org.gvsig.tools.library.LibraryException;

/**

 * Initialization of ClipLibrary library.

 * 

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

 */

public class ClipLibrary extends AlgorithmAbstractLibrary {

    @Override

    protected void doInitialize() throws LibraryException {

    }

    @Override

    protected void doPostInitialize() throws LibraryException {

        Messages.addResourceFamily("org.gvsig.geoprocess.algorithm.clip.clip",

            ClipLibrary.class.getClassLoader(), ClipLibrary.class.getClass()

                .getName());

        registerGeoProcess(ClipAlgorithm.class,

            getLanguageStrings("org.gvsig.geoprocess.algorithm.clip.clip"));

    }

}

org.gvsig.geoprocess.algorithm.clip.ClipLibrary

gvSIG_Algorithms=gvSIG algorithms

Clip_layer=Clip cover

Clip=Clip

Input_layer=Input cover

Selected_geometries=Selected features

gvSIG_Algorithms=Algoritmos gvSIG

Clip_layer=Capa de recorte

Clip=Recorte

Input_layer=Capa de entrada

Clip_layer=Capa de recorte

Selected_geometries=Geometrias seleccionadas

<?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.clip</artifactId>

  <packaging>jar</packaging>

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

	<parent>

		<groupId>org.gvsig</groupId>

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

		<version>1.0.0-SNAPSHOT</version>

	</parent>

	<dependencies>

		<dependency>

		    <groupId>org.gvsig</groupId>

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

   			<version>1.0.0-SNAPSHOT</version>

            <scope>compile</scope>

   		</dependency>

	</dependencies>

</project>

<assembly>

<!-- Empty distribution file -->

</assembly>

package org.gvsig.geoprocess.algorithm.convexhull;

import java.awt.Dimension;

import java.awt.GridBagConstraints;

import java.awt.GridBagLayout;

import java.awt.Insets;

import java.awt.event.ActionEvent;

import java.awt.event.ActionListener;

import java.util.ArrayList;

import java.util.Iterator;

import java.util.List;

import javax.swing.JButton;

import javax.swing.JComboBox;

import es.unex.sextante.core.GeoAlgorithm;

import es.unex.sextante.core.Sextante;

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

import org.slf4j.Logger;

import org.slf4j.LoggerFactory;

import org.gvsig.fmap.dal.DALLocator;

import org.gvsig.fmap.dal.DataManager;

import org.gvsig.fmap.dal.DataStoreParameters;

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

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

import org.gvsig.fmap.mapcontrol.swing.dynobject.DynObjectEditor;

import org.gvsig.tools.service.ServiceException;

/**

 * @deprecated

 * @author Nacho Brodin (nachobrodin@gmail.com)

 */

public class ConvexHullParametersPanel extends GeoAlgorithmParametersPanel implements ActionListener {

    private static final long serialVersionUID = 3260891769310434508L;

    private static final Logger LOG =

        LoggerFactory.getLogger(ConvexHullParametersPanel.class);

	private JComboBox                        combo               = null;

	private JButton                          button              = null;

	private ArrayList<DataStoreParameters>   paramList           = new ArrayList<DataStoreParameters>();

	public ConvexHullParametersPanel() {

		super();

	}

    public void init(GeoAlgorithm algorithm) {

    	initGUI();

    }

	private void initGUI() {

		GridBagLayout gbl = new GridBagLayout();

		this.setLayout(gbl);

		GridBagConstraints gbc = new GridBagConstraints();

		gbc.fill = GridBagConstraints.HORIZONTAL;

		gbc.weightx = 1.0;

		gbc.gridx = 0;

		gbc.gridy = 0;

		gbc.insets = new Insets(0, 2, 0, 0);

		this.add(getCombo(), gbc);

		gbc.fill = GridBagConstraints.NONE;

		gbc.weightx = 0;

		gbc.gridx = 1;

		gbc.insets = new Insets(0, 5, 0, 2);

		this.add(getButton(), gbc);

	}

	/**

	 * Gets a ComboBox

	 * @return

	 */

	public JComboBox getCombo() {

		if(combo == null) {

			combo = new JComboBox();

			loadProviderList(combo, paramList);

			combo.addActionListener(this);

			combo.setPreferredSize(new Dimension(0, 18));

		}

		return combo;

	}

	/**

	 * Gets a JButton

	 * @return

	 */

	public JButton getButton() {

		if(button == null) {

			button = new JButton("...");

			button.setPreferredSize(new Dimension(60, 18));

			button.addActionListener(this);

		}

		return button;

	}

	@Override

    public void assignParameters() {

        // Nothing to do

	}

	@Override

	public void setOutputValue(String arg0, String arg1) {

	}

	@Override

	public void setParameterValue(String arg0, String arg1) {

	}

	@SuppressWarnings("unchecked")

	private void loadProviderList(JComboBox c, ArrayList<DataStoreParameters> paramList) {

		try {

		DataManager manager = DALLocator.getDataManager();

		List list = manager.getStoreProviders();

		Iterator it = list.iterator();

		c.removeAllItems();

		paramList.clear();

		while(it.hasNext()) {

			try {

				String provider = it.next().toString();

				DataStoreParameters param = manager.createStoreParameters(provider);

				c.addItem(provider);

				paramList.add(param);

			} catch (InitializeException e1) {

				Sextante.addErrorToLog(e1);

			} catch (ProviderNotRegisteredException e1) {

				Sextante.addErrorToLog(e1);

			}

		}

		}catch(Exception e){}

	}

	public void actionPerformed(ActionEvent e) {

		if(e.getSource() == getButton()) {

			int index = getCombo().getSelectedIndex();

            try {

                DynObjectEditor editor = new DynObjectEditor(paramList.get(index));

                editor.editObject(true);

            } catch (ServiceException ex) {

                LOG.error(

                    "Error creating a Swing component for the DynObject: "

                        + paramList.get(index), ex);

                Sextante.addErrorToLog(ex);

            }

		}

		if(e.getSource() == getCombo()) {

		}

	}

}

/*

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

 *

 * Copyright (C) 2010 Generalitat Valenciana.

 *

 * This program is free software; you can redistribute it and/or

 * modify it under the terms of the GNU General Public License

 * as published by the Free Software Foundation; either version 2

 * of the License, or (at your option) any later version.

 *

 * This program is distributed in the hope that it will be useful,

 * but WITHOUT ANY WARRANTY; without even the implied warranty of

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 * GNU General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 *

 */

package org.gvsig.geoprocess.algorithm.convexhull;

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

import org.gvsig.i18n.Messages;

import org.gvsig.tools.library.LibraryException;

/**

 * Initialization of ConvexHullLibrary library.

 * 

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

 */

public class ConvexHullLibrary extends AlgorithmAbstractLibrary {

    @Override

    protected void doInitialize() throws LibraryException {

    }

    @Override

    protected void doPostInitialize() throws LibraryException {

        Messages.addResourceFamily(

            "org.gvsig.geoprocess.algorithm.convexhull.convexhull",

            ConvexHullLibrary.class.getClassLoader(), ConvexHullLibrary.class

                .getClass().getName());

        registerGeoProcess(

            ConvexHullAlgorithm.class,

            getLanguageStrings("org.gvsig.geoprocess.algorithm.convexhull.convexhull"));

    }

}

/*

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

 *

 * Copyright (C) 2010 Generalitat Valenciana.

 *

 * This program is free software; you can redistribute it and/or

 * modify it under the terms of the GNU General Public License

 * as published by the Free Software Foundation; either version 2

 * of the License, or (at your option) any later version.

 *

 * This program is distributed in the hope that it will be useful,

 * but WITHOUT ANY WARRANTY; without even the implied warranty of

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 * GNU General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 *

 */

package org.gvsig.geoprocess.algorithm.convexhull;

import java.util.Iterator;

import java.util.List;

import es.unex.sextante.core.Sextante;

import es.unex.sextante.dataObjects.IVectorLayer;

import es.unex.sextante.exceptions.GeoAlgorithmExecutionException;

import es.unex.sextante.exceptions.RepeatedParameterNameException;

import es.unex.sextante.outputs.OutputVectorLayer;

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

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

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

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

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

import org.gvsig.fmap.geom.Geometry;

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

import org.gvsig.geoprocess.core.gvVectorLayer;

import org.gvsig.geoprocess.sextante.AbstractSextanteGeoProcess;

import org.gvsig.tools.dispose.DisposableIterator;

/**

 * Convex Hull Algorithm

 * @author Nacho Brodin (nachobrodin@gmail.com)

 */

public class ConvexHullAlgorithm extends AbstractSextanteGeoProcess {

	public static final String  LAYER     = "LAYER";

	public static final String  RESULT    = "RESULT";

	public static final String  CHECK     = "CHECK";

	public void defineCharacteristics(){

		setName(Sextante.getText("Convex Hull"));

		setGroup(Sextante.getText("gvSIG_Algorithms"));

        // setGeneratesUserDefinedRasterOutput(false);

		try {

			m_Parameters.addInputVectorLayer(LAYER, 

												Sextante.getText("Input_layer"), 

												IVectorLayer.SHAPE_TYPE_WRONG, 

												true);

			addOutputVectorLayer(RESULT,

									Sextante.getText("ConvexHull"),

									OutputVectorLayer.SHAPE_TYPE_POLYGON);

			m_Parameters.addBoolean(CHECK, Sextante.getText("Selected_geometries"), false);

		} catch (RepeatedParameterNameException e) {

			Sextante.addErrorToLog(e);

		}

	}

	/*

	 * (non-Javadoc)

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

	 */

	@SuppressWarnings("unchecked")

	public boolean processAlgorithm() throws GeoAlgorithmExecutionException {

		IVectorLayer input = m_Parameters.getParameterValueAsVectorLayer(LAYER);

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

		FeatureStore store = null;

		if(input instanceof gvVectorLayer)

			store = ((gvVectorLayer)input).getFeatureStore();

		else

			return false;

		ConvexHullOperation convexHullOperation = new ConvexHullOperation();

		FeatureSet features = null;

		try {

			features = store.getFeatureSet();

			DisposableIterator it = null;

			if(selectedGeom) {

                FeatureSelection ds = store.getFeatureSelection();

                it = ds.iterator();

				//if(!it.hasNext())Mensaje de no hay geometrias seleccionadas

			} else

				it = features.iterator();

			int numberOfFeatures = (int)features.getSize();

			int iCount = 0;

			while( it.hasNext() ) {

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

				List geomList = feature.getGeometries();

				setProgress(iCount, numberOfFeatures);

				iCount ++;

				if(geomList == null) {

					Geometry geom = feature.getDefaultGeometry();

					if(geom != null)

						convexHullOperation.invoke(geom);

					continue;

				}

				Iterator<Geometry> itGeom = geomList.iterator();

				while(itGeom.hasNext()) {

					Geometry g = itGeom.next();

					convexHullOperation.invoke(g);

				}

			}

			Geometry g = convexHullOperation.getGeometry();

			if(g == null)

				return false;

			String[] sNames = {"ID"};

			Class [] types = {Integer.class};

			IVectorLayer output = getNewVectorLayer(RESULT,

													Sextante.getText("ConvexHull"),

													OutputVectorLayer.SHAPE_TYPE_POLYGON, types, sNames);

			com.vividsolutions.jts.geom.Geometry jtsGeom = GeometryUtil.geomToJTS(g);

			output.addFeature(jtsGeom, new Object[]{new Integer(0)});

			it.dispose();

		} catch (DataException e) {

			Sextante.addErrorToLog(e);

		}

		return !m_Task.isCanceled();

	}

}

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

 *

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

 *

 * This program is free software; you can redistribute it and/or

 * modify it under the terms of the GNU General Public License

 * as published by the Free Software Foundation; either version 2

 * of the License, or (at your option) any later version.

 *

 * This program is distributed in the hope that it will be useful,

 * but WITHOUT ANY WARRANTY; without even the implied warranty of

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 * GNU General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 */

package org.gvsig.geoprocess.algorithm.convexhull;

import javax.swing.JFrame;

import javax.swing.UIManager;

public class TestPanel {

	private JFrame 			            frame = new JFrame();

	@SuppressWarnings("deprecation")

	private ConvexHullParametersPanel	hp = null;

	public TestPanel() {

		super();

		initialize();

	}

	public static void main(String[] args){

		try {

			UIManager.setLookAndFeel("com.jgoodies.looks.plastic.PlasticXPLookAndFeel");

		} catch( Exception e ) {

			System.err.println( "No se puede cambiar al LookAndFeel");

		}

		new TestPanel();

	}

	@SuppressWarnings("deprecation")

	private void initialize() {

		frame.setSize(new java.awt.Dimension(370, 50));

		hp = new ConvexHullParametersPanel();

		hp.init(null);

		frame.setContentPane(hp);

		frame.setResizable(true);

		frame.setTitle("Panel");

		frame.setVisible(true);

		frame.addWindowListener(new java.awt.event.WindowAdapter() {

			public void windowClosing(java.awt.event.WindowEvent e) {

				System.exit(0);

			}

		});

	}

}

/*

 * Created on 16-feb-2006

 *

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

 *

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

 *

 * This program is free software; you can redistribute it and/or

 * modify it under the terms of the GNU General Public License

 * as published by the Free Software Foundation; either version 2

 * of the License, or (at your option) any later version.

 *

 * This program is distributed in the hope that it will be useful,

 * but WITHOUT ANY WARRANTY; without even the implied warranty of

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 * GNU General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 *

 * For more information, contact:

 *

 *  Generalitat Valenciana

 *   Conselleria d'Infraestructures i Transport

 *   Av. Blasco Ib��ez, 50

 *   46010 VALENCIA

 *   SPAIN

 *

 *      +34 963862235

 *   gvsig@gva.es

 *      www.gvsig.gva.es

 *

 *    or

 *

 *   IVER T.I. S.A

 *   Salamanca 50

 *   46005 Valencia

 *   Spain

 *

 *   +34 963163400

 *   dac@iver.es

 */

/* CVS MESSAGES:

*

* $Id: ScalableConvexHullVisitor.java 10626 2007-03-06 16:55:54Z caballero $

* $Log$

* Revision 1.2  2007-03-06 16:47:58  caballero

* Exceptions

*

* Revision 1.1  2006/06/20 18:20:45  azabala

* first version in cvs

*

* Revision 1.2  2006/06/02 18:21:28  azabala

* *** empty log message ***

*

* Revision 1.1  2006/05/24 21:13:31  azabala

* primera version en cvs despues de refactoring orientado a crear un framework extensible de geoprocessing

*

* Revision 1.3  2006/03/15 18:31:06  azabala

* *** empty log message ***

*

* Revision 1.2  2006/03/07 21:01:33  azabala

* *** empty log message ***

*

* Revision 1.1  2006/03/06 19:48:39  azabala

* *** empty log message ***

*

* Revision 1.2  2006/03/05 19:57:48  azabala

* *** empty log message ***

*

* Revision 1.1  2006/02/17 16:32:50  azabala

* *** empty log message ***

*

*

*/

package org.gvsig.geoprocess.algorithm.convexhull;

import com.vividsolutions.jts.geom.Geometry;

import com.vividsolutions.jts.geom.GeometryCollection;

import com.vividsolutions.jts.geom.GeometryFactory;

import es.unex.sextante.core.Sextante;

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

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

import org.gvsig.fmap.geom.GeometryLocator;

import org.gvsig.fmap.geom.GeometryManager;

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.geoprocess.algorithm.base.util.GeometryUtil;

/**

 * Convex Hull operation

 * @author Nacho Brodin (nachobrodin@gmail.com)

 */

public class ConvexHullOperation {

	public static final String  NAME             = "ConvexHull";

	public static final int     CODE             = GeometryLocator.getGeometryManager().getGeometryOperationCode(NAME);

	Geometry                    geometry         = null;

	GeometryFactory             geomFact         = null;

	private GeometryManager     geometryManager  = null;

	public ConvexHullOperation() {

		geometry = null;

		geometryManager = GeometryLocator.getGeometryManager();

	}

	public void setFeatureStore(FeatureStore out, String[] attrNames) throws DataException {

	}

	/**

	 * Returns FMap convex hull geometry.

	 * @return

	 */

	public org.gvsig.fmap.geom.Geometry getGeometry() {

		if(geometry == null)

			return null;

		GeometryOperationContext ctx = new GeometryOperationContext();

		ctx.setAttribute(FromJTS.PARAM, geometry);

		try {

			return (org.gvsig.fmap.geom.Geometry)geometryManager.invokeOperation(FromJTS.NAME, ctx);

		} catch (GeometryOperationNotSupportedException e) {

			Sextante.addErrorToLog(e);

			return null;

		} catch (GeometryOperationException e) {

			Sextante.addErrorToLog(e);

			return null;

		}

	}

	public org.gvsig.fmap.geom.Geometry invoke(org.gvsig.fmap.geom.Geometry g) {

		if(g == null)

			return null;

		com.vividsolutions.jts.geom.Geometry actualGeometry = GeometryUtil.geomToJTS(g);

		if(geometry == null)

			geometry = actualGeometry;

		else {

			com.vividsolutions.jts.geom.Geometry[] geoms = new com.vividsolutions.jts.geom.Geometry[2];

			geoms[0] = geometry;

			geoms[1] = actualGeometry;

			if(geomFact == null)

				geomFact = new GeometryFactory();

			GeometryCollection gc = geomFact.createGeometryCollection(geoms);

			geometry = gc.convexHull();	

		}

		return getGeometry();

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.fmap.geom.operation.GeometryOperation#getOperationIndex()

	 */

	public int getOperationIndex() {

		return CODE;

	}

}

org.gvsig.geoprocess.algorithm.convexhull.ConvexHullLibrary

gvSIG_Algorithms=gvSIG algorithms

Input_layer=Input cover

ConvexHull=Convex Hull

Selected_geometries=Selected features

gvSIG_Algorithms=Algoritmos gvSIG

Input_layer=Capa de entrada

ConvexHull=Convex Hull

Selected_geometries=Geometrias seleccionadas

<?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.convexhull</artifactId>