Application: gvSIG desktop » gvSIG geoprocess

#

# gvSIG. Desktop Geographic Information System.

#

# Copyright (C) 2007-2012 gvSIG Association.

#

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

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

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

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

#

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

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

# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

# GNU General Public License for more details.

#

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

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

# Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,

# MA  02110-1301, USA.

#

# For any additional information, do not hesitate to contact us

# at info AT gvsig.com, or visit our website www.gvsig.com.

#

basic_vect_algorithms=Capas vectoriales

Spatialjoin=Enlace espacial

Input_layer=Capa de entrada

selected_geometries=Geometrias seleccionadas (Capa de entrada 1)

use_the_nearest=Usar el m?s cercano

input_layers=Capas de entrada

input_layers_join=Capa de revestimiento

summary_function=Funci?n resumen

use_the_nearest=Usar el m?s cercano

Selected_geometries_input_layer=Geometr?as seleccionadas (capa de entrada)

Selected_geometries_overlay_layer=Geometr?as seleccionadas (capa de union)

#

# 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

Spatialjoin=Spatial Join

Input_layer=Input cover

selected_geometries=Selected geometries (input layer 1)

use_the_nearest=Use the nearest

input_layers=Input coverages

input_layers_join=Overlay layer

summary_function=Summary function

Selected_geometries_input_layer=Selected features (input layer)

Selected_geometries_overlay_layer=Selected features (overlay layer)

<?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, al igual que un join entre tablas, permite transferir los atributos de una capa a otra en base a una caracter&#237;stica com&#250;n. A diferencia del join de las bases de datos relacionales, en este caso la caracter&#237;stica com&#250;n no es un que un campo de las dos tablas tome el mismo valor (la clave del join), sino que los elementos relacionados de las dos capas cumplan unos criterios espaciales.&#10;&#10;El geoproceso Enlace Espacial implementado por la extensi&#243;n de geoprocesamiento de gvSIG permite seguir dos tipos de criterios espaciales para establecer el enlace espacial:&#10;&#10;&lt;UL&gt;&#10;&lt;LI&gt; Vecino m&#225;s pr&#243;ximo (relaci&#243;n 1-&gt;1). Asigna a un elemento de la capa origen los atributos del elemento m&#225;s pr&#243;ximo de la capa enlazada. En el caso de que el elemento m&#225;s pr&#243;ximo intersecte (o est&#233; contenido para el caso de pol&#237;gonos) al elemento original, habiendo por tanto varias intersecciones, el algoritmo tomar&#225; el primer elemento analizado de las posibles intersecciones.&lt;/LI&gt;&#10;&lt;LI&gt;Contenido en (relaci&#243;n 1-&gt;M). Relaciona un elemento de la capa origen con varios elementos de la capa destino (en concreto, con aquellos que son intersectados). En este caso la capa origen no heredar&#225; los atributos de la capa relacionada, sino que la operativa ser&#225; muy parecida a la del geoproceso Disolver. Para los M elementos relacionados con un elemento de la capa origen, se dar&#225; al usuario la posibilidad de escoger una o varias funciones resumen (media, m&#237;nimo, m&#225;ximo, sumatorio) que se aplicar&#225;n sobre los atributos num&#233;ricos de la capa enlazada.&lt;/LI&gt;&#10;&lt;/UL&gt;&#10;&#10;Este Geoproceso puede ser de gran utilidad para la realizaci&#243;n de an&#225;lisis de corredor. Por ejemplo:&#10;&lt;UL&gt;&#10;&lt;LI&gt;Qu&#233; zonas urbanos no tienen una parada de autobuses en un radio de 500 m.&lt;/LI&gt;&#10;&lt;LI&gt;Qu&#233; zonas urbanas carecen de colegios en un radio de 1000 m.&lt;/LI&gt;&#10;&lt;LI&gt;Qu&#233; pozos incumplen la normativa al no respetar la distancia m&#237;nima entre dos consecutivos.&lt;/LI&gt;&#10;&lt;/UL&gt;&#10;Etc." description="Descripci&#243;n" type="0">

			<image description="" file="spatialjoindesc.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="LAYER1" text="" description="Capa de entrada" type="3">

		</element>

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

		</element>

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

		</element>

		<element name="NEAREST" text="" description="Usar el m&#225;s cercano" type="3">

		</element>

		<element name="FUNCTION_LIST" text="" description="Function_list" type="3">

		</element>

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

		</element>

		<element name="RESULT" text="" description="Spatial Join" 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 transfer attributes from a source layer to a target layer (very similar to an alphanumerical join). The difference with a traditional join is that the transfer criteria is not the value of one (or many) given attributes. The criteria to transfer attribute values is spatial: intersection or proximity.&#10;&#10;If users select intersection like transfer criteria, each feature of the target layer could intersect many features of the source layer (a 1-N spatial relationship). If users select proximity (nearest neighbour) we have a 1-N spatial relationship." description="Descripci&#243;n" type="0">

			<image description="" file="spatialjoindesc.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="LAYER1" text="" description="Capa de entrada" type="3">

		</element>

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

		</element>

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

		</element>

		<element name="NEAREST" text="" description="Usar el m&#225;s cercano" type="3">

		</element>

		<element name="FUNCTION_LIST" text="" description="Function_list" type="3">

		</element>

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

		</element>

		<element name="RESULT" text="" description="Spatial Join" type="2">

		</element>

	</help>

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

org.gvsig.geoprocess.algorithm.spatialjoin.SpatialJoinLibrary

/**

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

import java.util.ArrayList;

import java.util.HashMap;

import java.util.Iterator;

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.exceptions.UnsupportedOutputChannelException;

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

import es.unex.sextante.outputs.OutputVectorLayer;

import org.gvsig.fmap.dal.DALLocator;

import org.gvsig.fmap.dal.DataManager;

import org.gvsig.fmap.dal.DataTypes;

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

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

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

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

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

import org.gvsig.fmap.geom.Geometry;

import org.gvsig.fmap.geom.GeometryLocator;

import org.gvsig.fmap.geom.GeometryManager;

import org.gvsig.fmap.geom.SpatialIndex;

import org.gvsig.fmap.geom.SpatialIndexFactory;

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

import org.gvsig.geoprocess.algorithm.dissolve.DissolveRule;

import org.gvsig.geoprocess.algorithm.dissolve.IDissolveRule;

import org.gvsig.geoprocess.algorithm.dissolve.Summary;

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

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

import org.gvsig.tools.dynobject.DynObject;

import org.gvsig.tools.exception.BaseException;

import org.gvsig.tools.visitor.VisitCanceledException;

import org.gvsig.tools.visitor.Visitor;

/**

 * Spatial join algorithm

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

 */

public class SpatialJoinAlgorithm extends AbstractSextanteGeoProcess {

	public static final String        NEW_FIELD             = "NUM_RELA";

	public static final String        RESULT                = "RESULT";

	public static final String        LAYER1                = "LAYER1";

	public static final String        LAYER2                = "LAYER2";

	public static final String        SELECTGEOM_INPUT      = "SELECTGEOM_INPUT";

	public static final String        SELECTGEOM_OVERLAY    = "SELECTGEOM_OVERLAY";

	public static final String        NEAREST               = "NEAREST";

	public static final String        FUNCTION_LIST         = "FUNCTION_LIST";

	public static final String        Summary[]             = {"Min", "Max", "Sum", "Avg"};

	private boolean                   funcList[]            = new boolean[Summary.length];

	private HashMap<String, String>   funcMap               = new HashMap<String, String>();

	public void defineCharacteristics(){

        setName(getTranslation("Spatialjoin"));

        setGroup(getTranslation("basic_vect_algorithms"));

        // setGeneratesUserDefinedRasterOutput(false);

		try {

			m_Parameters.addInputVectorLayer(LAYER1,

                getTranslation("Input_layer"),

												IVectorLayer.SHAPE_TYPE_WRONG,

												true);

			m_Parameters.addInputVectorLayer(LAYER2,

                getTranslation("Input_layer"),

												IVectorLayer.SHAPE_TYPE_WRONG,

												true);

            m_Parameters.addBoolean(SELECTGEOM_INPUT,

            		getTranslation("Selected_geometries_input_layer"), false);

            m_Parameters.addBoolean(SELECTGEOM_OVERLAY,

            		getTranslation("Selected_geometries_overlay_layer"), false);

            m_Parameters.addBoolean(NEAREST, getTranslation("use_the_nearest"), false);

            m_Parameters.addString(FUNCTION_LIST,

                getTranslation("Function_list"));

		} catch (RepeatedParameterNameException e) {

			Sextante.addErrorToLog(e);

		}

		addOutputVectorLayer(RESULT,getTranslation("Spatialjoin"),

								OutputVectorLayer.SHAPE_TYPE_UNDEFINED);

	}

	public boolean processAlgorithm() throws GeoAlgorithmExecutionException {

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

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

    	}

		IVectorLayer layer1 = m_Parameters.getParameterValueAsVectorLayer(LAYER1); //Capa de entrada

		IVectorLayer layer2 = m_Parameters.getParameterValueAsVectorLayer(LAYER2); //Capa de revestimiento

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

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

		boolean nearest = m_Parameters.getParameterValueAsBoolean(NEAREST);

		String functionList = m_Parameters.getParameterValueAsString(FUNCTION_LIST);

		loadSummary(functionList);

		FlyrVectIVectorLayer inputLayer = null;

		FlyrVectIVectorLayer overlayLayer = null;

		if(layer2 instanceof FlyrVectIVectorLayer && layer1 instanceof FlyrVectIVectorLayer) {

			overlayLayer = ((FlyrVectIVectorLayer)layer2); //Capa de revestimiento

			inputLayer = ((FlyrVectIVectorLayer)layer1); //Capa de entrada

		} else

			return false;

		DataManager dataManager = DALLocator.getDataManager();

		//Builds the output and computes the operation

		try {

			FeatureType featureTypeInputLayer = inputLayer.getFeatureStore().getDefaultFeatureType(); //Capa de entrada

			FeatureType featureTypeOverlayLayer = overlayLayer.getFeatureStore().getDefaultFeatureType(); //Capa de revestimiento

			GeometryOperation operation = null;

			FeatureStore outFeatStore = null;

            if (nearest) {

                outFeatStore =

                    buildOutPutStoreFromUnion(featureTypeInputLayer, featureTypeOverlayLayer,

                        inputLayer.getShapeType(), getTranslation("SpatialJoin"),

                        RESULT, "DIST", Double.class);

                operation =

                    new SpatiallyIndexedSpatialJoinOperation(overlayLayer, createIndex(overlayLayer.getFeatureStore(), selectedGeomOverlay), this);

            } else {

                outFeatStore =

                    buildSpatialJoinOutPutStore(featureTypeInputLayer,

                        inputLayer.getShapeType(), getTranslation("SpatialJoin"),

                        RESULT);

                IDissolveRule rule = new DissolveRule(0, funcMap);

                Summary summary =

                    new Summary(rule, outFeatStore.getDefaultFeatureType());

                operation =

                    new IntersectsSpatialJoinOperation(overlayLayer,

                        createIndex(overlayLayer.getFeatureStore(), selectedGeomOverlay), summary, this);

            }

            operation.setTaskStatus(getStatus());

			operation.computesGeometryOperation(inputLayer.getFeatureStore(),

					outFeatStore,

					attrNames,

					selectedGeomInput,

					selectedGeomOverlay,

					true);

		} catch (Exception e) {

			Sextante.addErrorToLog(e);

        }

		return true;

	}

	private SpatialIndex createIndex(FeatureStore store, boolean selectionOnly) throws BaseException{

	    GeometryManager geomManager = GeometryLocator.getGeometryManager();

	    SpatialIndexFactory factory = geomManager.getSpatialIndexFactory(geomManager.SPATIALINDEX_DEFAULT_RTREE);

	    DynObject parameters = factory.createParameters();

	    SpatialIndex index = (SpatialIndex) factory.create(parameters, geomManager);

	    final SpatialIndex wrappedIndex = store.wrapSpatialIndex(index);

	    Visitor visitor = new Visitor() {

            @Override

            public void visit(Object obj) throws VisitCanceledException, BaseException {

                Feature f=(Feature) obj;

                Geometry g = f.getDefaultGeometry();

                FeatureReference ref = f.getReference();

                wrappedIndex.insert(g, ref);

            }

        };

        if(selectionOnly){

            store.getFeatureSelection().accept(visitor);

        } else {

            store.accept(visitor);

        }

	    return wrappedIndex;

	}

	/**

	 * Checks if the parameter is in Summary list

	 * @param it

	 * @return the position in the list

	 */

	private int isInList(String it) {

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

			if(Summary[i].compareTo(it) == 0)

				return i;

		}

		return -1;

	}

	/**

	 * Loads the list of functions to use

	 * @param functionList

	 */

	private void loadSummary(String functionList) {

		String[] attrList = functionList.split(";");

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

			String[] func = attrList[i].split(",");

			for (int j = 1; j < func.length; j++) {

				int pos = isInList(func[j]);

				if(pos != -1) {

					funcList[pos] = true;

					funcMap.put(Summary[pos], func[0]);

				}

			}

		}

	}

	/**

	 * Builds the output FeatureStore

	 * @param featureType

	 * @return FeatureStore

	 */

	protected FeatureStore buildSpatialJoinOutPutStore(FeatureType featureType1,

											int shapeType,

											String sextanteLayerName,

											String sextanteLayerLabel) {

		ArrayList<Class> typesList = new ArrayList<Class>();

		ArrayList<String> attr = new ArrayList<String>();

		Iterator it = featureType1.iterator();

		while( it.hasNext() ) {

			FeatureAttributeDescriptor attribute = (FeatureAttributeDescriptor)it.next();

			if (attribute.getDataType().getType() != DataTypes.GEOMETRY) {

				attr.add(attribute.getName());

				typesList.add(attribute.getObjectClass());

			}

		}

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

			if(funcList[i]) {

				String fieldName = funcMap.get(Summary[i]);

				if(fieldName.length() >= 6)

					fieldName = fieldName.substring(0, 7);

				attr.add(fieldName + "_" + Summary[i]);

				typesList.add(Double.class);

			}

		}

		attr.add(NEW_FIELD);

		typesList.add(Integer.class);

		attrNames = new String[attr.size()];

		attr.toArray(attrNames);

		Class[] types = new Class[typesList.size()];

		typesList.toArray(types);

		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;

	}

    @Override

    public Class<? extends GeoAlgorithmParametersPanel> getCustomParametersPanelClass() {

        return SpatialJoinParametersPanel.class;

    }

}

/**

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2012 gvSIG Association.

 *

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

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

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

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

 *

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

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

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 * GNU General Public License for more details.

 *

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

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

 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,

 * MA  02110-1301, USA.

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

package org.gvsig.geoprocess.algorithm.spatialjoin;

import java.awt.BorderLayout;

import java.awt.Color;

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

import javax.swing.ComboBoxModel;

import javax.swing.DefaultComboBoxModel;

import javax.swing.JCheckBox;

import javax.swing.JComboBox;

import javax.swing.JLabel;

import javax.swing.JPanel;

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

import org.gvsig.geoprocess.sextante.gui.algorithm.AlgorithmOutputPanel;

import org.gvsig.gui.beans.table.TableContainer;

import org.gvsig.gui.beans.table.exceptions.NotInitializeException;

import es.unex.sextante.core.GeoAlgorithm;

import es.unex.sextante.core.ObjectAndDescription;

import es.unex.sextante.core.OutputObjectsSet;

import es.unex.sextante.core.ParametersSet;

import es.unex.sextante.core.Sextante;

import es.unex.sextante.dataObjects.IVectorLayer;

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

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

import es.unex.sextante.gui.core.SextanteGUI;

import es.unex.sextante.outputs.Output;

/**

 * Panel for dissolve algorithm

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

 */

public class SpatialJoinParametersPanel extends GeoAlgorithmParametersPanel implements ActionListener {

	private static final long                serialVersionUID   = 1L;

	private GeoAlgorithm                     m_Algorithm        = null;

	private JComboBox                        layersCombo        = null;

	private JComboBox                        layersJoinCombo    = null;

	private JCheckBox                        selectionInput     = null;

	private JCheckBox                        selectionOverlay   = null;

	private JCheckBox                        nearest       = null;

	private AlgorithmOutputPanel             output             = null;

	private final String[]                   columnNames        = { "Min", "Max", "Sum", "Avg", "Field ID" };

	private final int[]                      columnWidths       = { 35, 35, 35, 35, 334 };

	private TableContainer                   table              = null;

	private OutputChannelSelectionPanel      outputChannelSelectionPanel;

	private JPanel                           outputPanel;

	public SpatialJoinParametersPanel() {

		super();

	}

    public void init(GeoAlgorithm algorithm) {

    	m_Algorithm = algorithm;

    	initGUI();

    }

	private void initGUI() {

		GridBagLayout gbl = new GridBagLayout();

		this.setLayout(gbl);

		this.setBorder(BorderFactory.createLineBorder(Color.gray));

		GridBagConstraints gbc = new GridBagConstraints();

		gbc.fill = GridBagConstraints.HORIZONTAL;

		gbc.weightx = 1.0;

		gbc.gridx = 0;

		gbc.gridy = 0;

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

		this.add(getComboPanel(GeoProcessLocator.getGeoProcessManager().getTranslation("input_layers"), getLayers1Combo()), gbc);

		gbc.gridy = 1;

		this.add(getComboPanel(GeoProcessLocator.getGeoProcessManager().getTranslation("input_layers_join"), getLayers2Combo()), gbc);

		gbc.gridy = 2;

		this.add(getSelectionInputCheck(), gbc);

		gbc.gridy = 3;

		this.add(getSelectionOverlayCheck(), gbc);

		gbc.gridy = 4;

		this.add(getNearestCheck(), gbc);

		gbc.gridy = 5;

		this.add(new JLabel(GeoProcessLocator.getGeoProcessManager().getTranslation("summary_function")), gbc);

		gbc.gridy = 6;

		gbc.fill = GridBagConstraints.BOTH;

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

		gbc.weighty = 1.0;

		this.add(getRadioButtonTable(), gbc);

		gbc.gridy = 7;

		gbc.fill = GridBagConstraints.HORIZONTAL;

		gbc.weighty = 0.0;

		this.add(getOutputChannelSelectionPanel(), gbc);

		initTable();

		getRadioButtonTable().setVisible(false);

	}

	/**

	 * Gets the output panel (SEXTANTE)

	 * @return

	 */

	private JPanel getOutputChannelSelectionPanel() {

		if(outputPanel == null) {

			try {

				outputPanel = new JPanel();

				outputPanel.setLayout(new BorderLayout());

				final OutputObjectsSet ooSet = m_Algorithm.getOutputObjects();

				final Output out = ooSet.getOutput(SpatialJoinAlgorithm.RESULT);

				outputChannelSelectionPanel = new OutputChannelSelectionPanel(out, m_Algorithm.getParameters());

				outputPanel.add(new JLabel(GeoProcessLocator.getGeoProcessManager().getTranslation("Spatialjoin") +

						" [" + GeoProcessLocator.getGeoProcessManager().getTranslation("Vectorial") +"]                        "), BorderLayout.WEST);

				outputPanel.add(outputChannelSelectionPanel, BorderLayout.CENTER);

			} catch (final Exception e) {

				Sextante.addErrorToLog(e);

			}

		}

		return outputPanel;

	}

	/**

	 * Gets the output panel (DAL)

	 * @return

	 */

	@SuppressWarnings("unused")

	private AlgorithmOutputPanel getAlgorithmOutputPanel() {

		if(output == null) {

			output = new AlgorithmOutputPanel();

		}

		return output;

	}

	/**

	 * Gets a new JPanel with the text and JComboBox

	 * @param text

	 * @param combo

	 * @return

	 */

	public JPanel getComboPanel(String text, JComboBox combo) {

		JPanel panel = new JPanel();

		GridBagLayout gbl = new GridBagLayout();

		panel.setLayout(gbl);

		GridBagConstraints gbc = new GridBagConstraints();

		gbc.fill = GridBagConstraints.NONE;

		gbc.weightx = 0;

		gbc.gridx = 0;

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

		JLabel label = new JLabel(text);

		label.setPreferredSize(new Dimension(180, 18));

		panel.add(label, gbc);

		gbc.fill = GridBagConstraints.HORIZONTAL;

		gbc.weightx = 1.0;

		gbc.gridx = 1;

		gbc.anchor = GridBagConstraints.EAST;

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

		panel.add(combo, gbc);

		return panel;

	}

	/**

	 * Gets a ComboBox

	 * @return

	 */

	public JComboBox getLayers1Combo() {

		if(layersCombo == null) {

			layersCombo = new JComboBox();

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

			ComboBoxModel comboModel = new DefaultComboBoxModel(getLayerList());

			layersCombo.setModel(comboModel);

			layersCombo.addActionListener(this);

		}

		return layersCombo;

	}

	/**

	 * Gets a ComboBox

	 * @return

	 */

	public JComboBox getLayers2Combo() {

		if(layersJoinCombo == null) {

			layersJoinCombo = new JComboBox();

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

			ComboBoxModel comboModel = new DefaultComboBoxModel(getLayerList());

			layersJoinCombo.setModel(comboModel);

			layersJoinCombo.addActionListener(this);

		}

		return layersJoinCombo;

	}

	/**

	 * Gets a CheckBox

	 * @return

	 */

	public JCheckBox getSelectionInputCheck() {

		if(selectionInput == null) {

			selectionInput = new JCheckBox(GeoProcessLocator.getGeoProcessManager().getTranslation("Selected_geometries_input_layer"));

		}

		return selectionInput;

	}

	/**

	 * Gets a CheckBox

	 * @return

	 */

	public JCheckBox getSelectionOverlayCheck() {

		if(selectionOverlay == null) {

			selectionOverlay = new JCheckBox(GeoProcessLocator.getGeoProcessManager().getTranslation("Selected_geometries_overlay_layer"));

		}

		return selectionOverlay;

	}

	/**

	 * Gets a CheckBox

	 * @return

	 */

	public JCheckBox getNearestCheck() {

		if(nearest == null) {

			nearest = new JCheckBox(GeoProcessLocator.getGeoProcessManager().getTranslation("use_the_nearest"));

			nearest.setSelected(true);

			nearest.addActionListener(this);

		}

		return nearest;

	}

	/**

	 * Gets the summary table

	 * @return TableContainer

	 */

	public TableContainer getRadioButtonTable() {

		if (table == null) {

			table = new TableContainer(columnNames, columnWidths, null);

			table.setModel("ARGBBandSelectorModel");

			table.setControlVisible(false);

			table.initialize();

		}

		return table;

	}

	//------------------------------------------------------------

	public void actionPerformed(ActionEvent e) {

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

			if(!getNearestCheck().isSelected())

				getRadioButtonTable().setVisible(true);

			else

				getRadioButtonTable().setVisible(false);

		}

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

			initTable();

			//If the second layer has not numerical fields only the nearest method can be applied

			IVectorLayer lyr = getSelectedVectorLayer2();

			Class[] types = lyr.getFieldTypes();

			boolean hasNumericField = false;

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

				if(types[i] == Integer.class

						|| types[i] == Double.class

						|| types[i] == Float.class

						|| types[i] == Short.class) {

					hasNumericField = true;

					break;

				}

			}

			if(!hasNumericField){

				getNearestCheck().setSelected(true);

			}

		}

	}

	/**

	 * Adds to the table one entry for each field

	 */

	private void initTable() {

		try {

			getRadioButtonTable().removeAllRows();

			if(getSelectedVectorLayer2() == null)

				return;

			Class<?>[] types = getSelectedVectorLayer2().getFieldTypes();

			for (int i = 0; i < getSelectedVectorLayer2().getFieldCount(); i++) {

				if(Number.class.isAssignableFrom(types[i])) {

					addTableRow(getSelectedVectorLayer2().getFieldName(i));

				}

			}

		} catch (NotInitializeException e) {

			Sextante.addErrorToLog(e);

		}

	}

	/**

	 * A?ade una banda a la tabla bandas de la imagen asignandole un nombre y

	 * valor a los checkbox

	 * @param bandName Nombre de la banda

	 * @throws NotInitializeException

	 */

	private void addTableRow(String fieldName) throws NotInitializeException {

		Object[] row = {	new Boolean(false),

							new Boolean(false),

							new Boolean(false),

							new Boolean(false),

							fieldName };

		getRadioButtonTable().addRow(row);

	}

	@Override

	public void assignParameters() {

		try {

			ParametersSet params = m_Algorithm.getParameters();

			params.getParameter(SpatialJoinAlgorithm.LAYER1).setParameterValue(getSelectedVectorLayer());

			params.getParameter(SpatialJoinAlgorithm.LAYER2).setParameterValue(getSelectedVectorLayer2());

			params.getParameter(SpatialJoinAlgorithm.SELECTGEOM_INPUT).setParameterValue(getSelectionInputCheck().isSelected());

			params.getParameter(SpatialJoinAlgorithm.SELECTGEOM_OVERLAY).setParameterValue(getSelectionOverlayCheck().isSelected());

			params.getParameter(SpatialJoinAlgorithm.NEAREST).setParameterValue(getNearestCheck().isSelected());

			params.getParameter(SpatialJoinAlgorithm.FUNCTION_LIST).setParameterValue(getValues());

			OutputObjectsSet ooSet = m_Algorithm.getOutputObjects();

			Output out = ooSet.getOutput(SpatialJoinAlgorithm.RESULT);

			//Reponer estas l?neas para cambiar el panel de salida y comentar la siguiente

			//AlgorithmOutputPanel fsp = getAlgorithmOutputPanel();

			//out.setOutputChannel(new CompositeSourceOutputChannel(fsp.getOutputParameters()));

	         out.setOutputChannel(outputChannelSelectionPanel.getOutputChannel());

		} catch (Exception e) {

			Sextante.addErrorToLog(e);

        }

	}

	@Override

	public void setOutputValue(String arg0, String arg1) {

	}

	@Override

	public void setParameterValue(String arg0, String arg1) {

	}

	/**

	 * Gets the input layer list

	 * @return

	 */

	private ObjectAndDescription[] getLayerList() {

		IVectorLayer[] layers = SextanteGUI.getInputFactory()

					.getVectorLayers(IVectorLayer.SHAPE_TYPE_WRONG);

		ObjectAndDescription[] oad = new ObjectAndDescription[layers.length];

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

			oad[i] = new ObjectAndDescription(layers[i].getName(), layers[i]);

		return oad;

	}

	/**

	 * Gets the selected vector layer in the JComboBox

	 * @return

	 */

	private IVectorLayer getSelectedVectorLayer() {

		if(layersCombo.getSelectedItem() != null)

			return (IVectorLayer)((ObjectAndDescription)layersCombo.getSelectedItem()).getObject();

		return null;

	}

	/**

	 * Gets the selected vector layer in the JComboBox

	 * @return

	 */

	private IVectorLayer getSelectedVectorLayer2() {

		if(layersJoinCombo.getSelectedItem() != null)

			return (IVectorLayer)((ObjectAndDescription)layersJoinCombo.getSelectedItem()).getObject();

		return null;

	}

	/**

	 * Gets the field list of the selected layer

	 * @return

	 */

	public String[] getFieldList() {

		IVectorLayer layer = getSelectedVectorLayer2();

		String[] data = new String[layer.getFieldCount()];

		for (int i = 0; i < layer.getFieldCount(); i++)

			data[i] = layer.getFieldName(i);

		return data;

	}

	/**

	 * Formats the content of the table

	 * @return

	 */

	private String getValues() {

    	String str = "";

    	try {

			for (int i = 0; i < getRadioButtonTable().getRowCount(); i++) {

				str = str + (String)getRadioButtonTable().getModel().getValueAt(i, 4) + ",";

				for (int j = 0; j < getRadioButtonTable().getModel().getColumnCount() - 1; j++)

					if(((Boolean)getRadioButtonTable().getModel().getValueAt(i, j)).booleanValue())

						str = str + SpatialJoinAlgorithm.Summary[j] + ",";

				str = str.substring(0, str.length() - 1) + ";";

			}

		} catch (NotInitializeException e) {

			Sextante.addErrorToLog(e);

		}

    	if (str.isEmpty()) {

    		return str;

    	}

    	return str.substring(0, str.length() - 1);

    }

}

/**

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

import java.util.List;

import es.unex.sextante.core.Sextante;

import java.util.Iterator;

import org.apache.commons.lang3.mutable.MutableBoolean;

import org.apache.commons.lang3.mutable.MutableInt;

import org.jfree.base.config.ModifiableConfiguration;

import org.gvsig.fmap.dal.DataTypes;

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.FeatureAttributeDescriptor;

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

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

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

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.dal.feature.exception.FeatureIndexException;

import org.gvsig.fmap.geom.Geometry;

import org.gvsig.fmap.geom.SpatialIndex;

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

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

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

import org.gvsig.geoprocess.algorithm.dissolve.Summary;

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

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

import org.gvsig.tools.exception.BaseException;

import org.gvsig.tools.visitor.VisitCanceledException;

import org.gvsig.tools.visitor.Visitor;

/**

 *

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

 */

public class IntersectsSpatialJoinOperation extends GeometryOperation {

	/**

	 * Specialized instance in nearest neighbor search.

	 */

	private SpatialIndex              index                 = null;

	private Summary                   summary               = null;

	private FeatureStore              storeOverlay          = null;

	private FeatureSelection          featureSelection      = null;

	public IntersectsSpatialJoinOperation(FlyrVectIVectorLayer targetLayer, SpatialIndex index, Summary summary, AbstractSextanteGeoProcess p) {

		super(p);

	    this.index = index;

		this.summary = summary;

		storeOverlay = targetLayer.getFeatureStore();

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.geoprocess.algorithm.base.core.GeometryOperation#invoke(org.gvsig.fmap.geom.Geometry, org.gvsig.fmap.dal.feature.Feature)

	 */

	@SuppressWarnings("deprecation")

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

		if(g == null)

			return lastEditFeature;

		if(featureSelection == null) {

			try {

				featureSelection = storeOverlay.getFeatureSelection();

			} catch (DataException e1) {

				//Sin selecci?n tiramos

			}

		}

		try {

	          final MutableBoolean first = new MutableBoolean(true);

	          final MutableInt numReg = new MutableInt(0);

            index.query(g, new Visitor() {

                @Override

                public void visit(Object obj) throws VisitCanceledException, BaseException {

                    FeatureReference ref = (FeatureReference) obj;

                    Feature feat = ref.getFeature();

                    Geometry g2 = feat.getDefaultGeometry();

                    if(g.intersects(g2)){

                        if(first.getValue()) {

                            summary.loadDefaultSummarizes(feat);

                            first.setFalse();

                        }

                        summary.updateValues(feat);

                        numReg.increment();

                    }

                }

            });

			buildNewFeature(featureInput, summary, numReg.getValue());

		} catch(FeatureIndexException e) {

			Sextante.addErrorToLog(e);

		} catch (DataException e) {

			Sextante.addErrorToLog(e);

		}

		return lastEditFeature;

	}

	/**

	 * Builds a new output feature and adds it to the output file

	 * @param feat1

	 * @param feat2

	 * @param value

	 * @param g

	 * @throws DataException

	 */

	@SuppressWarnings("unchecked")

	private void buildNewFeature(Feature feat, Summary summary, int numReg) throws DataException {

		EditableFeature outFeat = persister.getOutputFeatureStore().createNewFeature();

		FeatureAttributeDescriptor[] attrs = feat.getType().getAttributeDescriptors();

		//Loads the old feature

		int pos = 0;

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

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

				outFeat.set(pos, feat.get(i));

				pos ++;

			}

		}

		//Loads the statistics

		summary.loadEditableFeature(outFeat);

		//Loads the new field

		outFeat.set(SpatialJoinAlgorithm.NEW_FIELD, new Integer(numReg));

		//Saves the geometry

		try {

			List l = feat.getGeometries();

			if(l == null) {

				persister.addFeature(outFeat, feat.getDefaultGeometry());

			} else {

				EditableFeature editFeat = null;

				for (int i = 0; i < l.size(); i++) {

					if(editFeat == null) {

						editFeat = persister.addFeature(outFeat, (org.gvsig.fmap.geom.Geometry)l.get(i));

					} else {

						persister.addGeometryToExistingFeature(editFeat, (org.gvsig.fmap.geom.Geometry)l.get(i));

					}

				}

			}

		} catch (CreateGeometryException e) {

			Sextante.addErrorToLog(e);

		}

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.geoprocess.algorithm.base.core.GeometryOperation#invoke(org.gvsig.fmap.geom.Geometry, org.gvsig.fmap.dal.feature.EditableFeature)

	 */

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

	}

}

/**

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

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

import org.gvsig.i18n.Messages;

import org.gvsig.tools.library.LibraryException;

/**

 * Initialization of SpatialJoinLibrary library.

 * 

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

 */

public class SpatialJoinLibrary extends AlgorithmAbstractLibrary {

    @Override

    protected void doInitialize() throws LibraryException {

    }

    @Override

    protected void doPostInitialize() throws LibraryException {

        Messages.addResourceFamily(

            "org.gvsig.geoprocess.algorithm.spatialjoin.spatialjoin",

            SpatialJoinLibrary.class.getClassLoader(), SpatialJoinLibrary.class

                .getClass().getName());

        registerGeoProcess(new SpatialJoinAlgorithm());

    }

}

/**

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

import java.util.Iterator;

import org.apache.commons.lang3.mutable.MutableDouble;

import org.apache.commons.lang3.mutable.MutableObject;

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.FeatureIndex;

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

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

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.fmap.geom.GeometryLocator;

import org.gvsig.fmap.geom.GeometryManager;

import org.gvsig.fmap.geom.SpatialIndex;

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

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

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

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

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

import org.gvsig.tools.exception.BaseException;

import org.gvsig.tools.visitor.VisitCanceledException;

import org.gvsig.tools.visitor.Visitor;

/**

 *

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

 */

public class SpatiallyIndexedSpatialJoinOperation extends GeometryOperation {

	private FeatureStore               storeOverlay        = null;

	private FeatureSelection           featureSelection    = null;

	/**

	 * Geometry.distance() is a costly operation. Thats the reason for we are

	 * only using a nearest neighbor. TODO SpatialIndex works with Rectangle2D,

	 * and this may be a simplification that drives to errors. Make additional

	 * probes to use a number of default neighbors

	 */

	static final int                  DEFAULT_NUM_NEIGBOURS = 1;

	/**

	 * Number of neighbors that nearestFinder must found.

	 */

	int                               numOfNeighbours       = DEFAULT_NUM_NEIGBOURS;

	/**

	 * Specialized instance in nearest neighbor search.

	 */

	private SpatialIndex              index                 = null;

	public SpatiallyIndexedSpatialJoinOperation(FlyrVectIVectorLayer targetLayer, SpatialIndex index, AbstractSextanteGeoProcess p) {

		super(p);

		this.index = index;

		storeOverlay = targetLayer.getFeatureStore();

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.geoprocess.algorithm.base.core.GeometryOperation#invoke(org.gvsig.fmap.geom.Geometry, org.gvsig.fmap.dal.feature.Feature)

	 */

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

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

	    final MutableObject<Feature> foundFeature = new MutableObject<Feature>(null);

        boolean addedFeature = false;

        if (g == null)

            return lastEditFeature;

        try {

            final MutableDouble nearestDistance = new MutableDouble(Double.POSITIVE_INFINITY);

            index.query(g, new Visitor() {

                @Override

                public void visit(Object obj) throws VisitCanceledException, BaseException {

                    FeatureReference ref = (FeatureReference) obj;

                    Feature feat = ref.getFeature();

                    Geometry g2 = feat.getDefaultGeometry();

                    double dist = g.distance(g2);

                    if (dist <= nearestDistance.getValue()) {

                        nearestDistance.setValue(dist);

                        foundFeature.setValue(feat.getCopy());

                    }

                }

            });

            Iterator<?> iterator = index.queryNearest(g);

            while (iterator.hasNext()) {

                FeatureReference ref = (FeatureReference) iterator.next();

                Feature feat = ref.getFeature();

                Geometry g2 = feat.getDefaultGeometry();

                double dist = g.distance(g2);

                if (dist <= nearestDistance.getValue()) {

                    nearestDistance.setValue(dist);

                    foundFeature.setValue(feat.getCopy());

                }

            }

            if(foundFeature!=null){

                buildFeature(featureInput, foundFeature.getValue(), new Double(nearestDistance.getValue()), g);

                addedFeature = true;

            }

        } catch (Exception e) {

            Sextante.addErrorToLog(e);

        }

        return lastEditFeature;

    }

	/**

	 * Builds a feature and adds it to the output file

	 * @param feat1

	 * @param feat2

	 * @param value

	 * @param g

	 * @throws DataException

	 */

	private void buildFeature(Feature feat1, Feature feat2, Object value, org.gvsig.fmap.geom.Geometry g) throws DataException {

		EditableFeature outFeat = persister.getOutputFeatureStore().createNewFeature();

		int sizeFeat1 = feat1.getType().size() - 1;

		for (int i = 0; i < sizeFeat1; i++)

			outFeat.set(i, feat1.get(i));

		for (int i = sizeFeat1; i < sizeFeat1 + feat2.getType().size() - 1; i++)

			outFeat.set(i, feat2.get(i - sizeFeat1));

		outFeat.set(outFeat.getType().size() - 2, value);

		try {

			persister.addFeature(outFeat, g);

		} catch (CreateGeometryException e) {

			Sextante.addErrorToLog(e);

		}

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.geoprocess.algorithm.base.core.GeometryOperation#invoke(org.gvsig.fmap.geom.Geometry, org.gvsig.fmap.dal.feature.EditableFeature)

	 */

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

	}

}

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

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

  <modelVersion>4.0.0</modelVersion>

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

  <packaging>jar</packaging>

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

	<parent>

		<groupId>org.gvsig</groupId>

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

		<version>2.2.107</version>

	</parent>

	<dependencies>

		<dependency>

		    <groupId>org.gvsig</groupId>

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

            <scope>compile</scope>

   		</dependency>

   		<dependency>

		    <groupId>org.gvsig</groupId>

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

            <scope>compile</scope>

   		</dependency>

	</dependencies>

</project>

org.gvsig.geoprocess.algorithm.lateralbuffer.LateralBufferLibrary

#

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

#

lateral_buffer=Buffer lateral

left=Izquierda

right=Derecha

select_lateral=Seleccionar lateral

SplitLines._Introduccion_de_datos=Cortar l\u00edneas. Introducci\u00f3n de datos

split_lines=Cortar l\u00edneas

distance_section=Distancia de la secci\u00f3n

Mensaje_split=Geoproceso split....

Mensaje_procesando_split=Partiendo geometr\u00edas...

new_geoprocess=Geoprocesos nuevos

concentrical_rings=N\u00famero de anillos concentricos

ignore_invalid_lines=Ignorar l\u00edneas inv\u00e1lidas

Use_simplified_geometries=Usar geometr\u00edas simplificadas

#

# 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