Application: gvSIG desktop » gvSIG geoprocess

/**

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2012 gvSIG Association.

 *

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

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

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

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

 *

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

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

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 * GNU General Public License for more details.

 *

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

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

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

 * MA  02110-1301, USA.

 *

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

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

 */

package org.gvsig.geoprocess.algorithm.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 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.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.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.dispose.DisposableIterator;

/**

 * 

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

 */

public class IntersectsSpatialJoinOperation extends GeometryOperation {

	/**

	 * Specialized instance in nearest neighbor search.

	 */

	private FeatureIndex              index                 = null;

	private Summary                   summary               = null;

	private FeatureStore              storeOverlay          = null;

	private FeatureSelection          featureSelection      = null;

	public IntersectsSpatialJoinOperation(FlyrVectIVectorLayer targetLayer, String indexName, Summary summary, AbstractSextanteGeoProcess p) {

		super(p);

		FeatureIndexes indexes = targetLayer.getFeatureStore().getIndexes();

		index = indexes.getFeatureIndex(indexName);

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

			Envelope rect = g.getEnvelope();

			FeatureSet featSet = index.getMatchFeatureSet(rect);

			boolean first = true;

			int numReg = 0;

			DisposableIterator it = featSet.iterator();

			while(it.hasNext()) {

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

				if( featureSelection != null && 

					!featureSelection.isEmpty() && 

					!featureSelection.isSelected(feat))

					continue;

				if(first) {

					summary.loadDefaultSummarizes(feat);

					first = false;

				}

				summary.updateValues(feat);

				numReg ++;

			}

			buildNewFeature(featureInput, summary, numReg);

		} 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.ArrayList;

import java.util.Iterator;

import java.util.List;

import com.vividsolutions.jts.geom.Geometry;

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

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

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

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

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

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

import org.gvsig.tools.dispose.DisposableIterator;

/**

 *

 * @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 FeatureIndex              index                 = null;

	public SpatiallyIndexedSpatialJoinOperation(FlyrVectIVectorLayer targetLayer, String indexName, AbstractSextanteGeoProcess p) {

		super(p);

		FeatureIndexes indexes = targetLayer.getFeatureStore().getIndexes();

		index = indexes.getFeatureIndex(indexName);

		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(org.gvsig.fmap.geom.Geometry g, Feature featureInput) {

	    boolean addedFeature = false;

		if(g == null)

			return lastEditFeature;

		if(featureSelection == null) {

			try {

				featureSelection = storeOverlay.getFeatureSelection();

			} catch (DataException e1) {

				//Sin selecci?n tiramos

			}

		}

		try {

			Geometry gJts = GeometryUtil.geomToJTS(g);

            Envelope rect = null;

		    rect = g.getEnvelope();

			FeatureSet featSetMatch = index.getMatchFeatureSet(rect);

			FeatureSet featSetNearest = index.getNearestFeatureSet(numOfNeighbours, rect);

			double nearestDistance = Double.MAX_VALUE;

            DisposableIterator it = featSetMatch.iterator();

            //Recorremos el featSetMatch para calcular la nearestDistance

            while(it.hasNext()) {

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

                if( featureSelection != null &&

                    !featureSelection.isEmpty() &&

                    !featureSelection.isSelected(feat))

                    continue;

                List geomList = feat.getGeometries();

                if(geomList == null) {

                    org.gvsig.fmap.geom.Geometry g2 = feat.getDefaultGeometry();

                    Geometry g2Jts = GeometryUtil.geomToJTS(g2);

                    double dist = gJts.distance(g2Jts);

                    if (dist <= nearestDistance) {

                        nearestDistance = dist;

                    }

                    continue;

                }

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

                while(itGeom.hasNext()) {

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

                    Geometry g2Jts = GeometryUtil.geomToJTS(g2);

                    double dist = gJts.distance(g2Jts);

                    if (dist <= nearestDistance) {

                        nearestDistance = dist;

                    }

                }

            }

            it.dispose();

            //Recorremos el featSetNearest para calcular la nearestDistance

            it = featSetNearest.iterator();

            while(it.hasNext()) {

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

                if( featureSelection != null &&

                    !featureSelection.isEmpty() &&

                    !featureSelection.isSelected(feat))

                    continue;

                List geomList = feat.getGeometries();

                if(geomList == null) {

                    org.gvsig.fmap.geom.Geometry g2 = feat.getDefaultGeometry();

                    Geometry g2Jts = GeometryUtil.geomToJTS(g2);

                    double dist = gJts.distance(g2Jts);

                    if (dist <= nearestDistance) {

                        nearestDistance = dist;

                    }

                    continue;

                }

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

                while(itGeom.hasNext()) {

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

                    Geometry g2Jts = GeometryUtil.geomToJTS(g2);

                    double dist = gJts.distance(g2Jts);

                    if (dist <= nearestDistance) {

                        nearestDistance = dist;

                    }

                }

            }

            it.dispose();

            //Recorremos el featSetMatch para obtener las geometr?as que est?n a una distancia igual o menor a la nearestDistance

            it = featSetMatch.iterator();

            while (it.hasNext()) {

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

                if (featureSelection != null && !featureSelection.isEmpty()

                    && !featureSelection.isSelected(feat))

                    continue;

                List geomList = feat.getGeometries();

                if (geomList == null) {

                    org.gvsig.fmap.geom.Geometry g2 = feat.getDefaultGeometry();

                    Geometry g2Jts = GeometryUtil.geomToJTS(g2);

                    double dist = gJts.distance(g2Jts);

                    if (dist <= nearestDistance) {

                        buildFeature(featureInput, feat, new Double(

                            nearestDistance), g);

                        addedFeature= true;

                    }

                     continue;

                }

                Iterator<org.gvsig.fmap.geom.Geometry> itGeom =

                    geomList.iterator();

                while (itGeom.hasNext()) {

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

                    Geometry g2Jts = GeometryUtil.geomToJTS(g2);

                    double dist = gJts.distance(g2Jts);

                    if (dist <= nearestDistance) {

                        buildFeature(featureInput, feat, new Double(

                            nearestDistance), g);

                        addedFeature= true;

                    }

                }

            }

            it.dispose();

            //Recorremos el featSetNearest para obtener las geometr?as que est?n a una distancia igual o menor a la nearestDistance

            //si no se ha a?adido todav?a ninguna feature

            if (!addedFeature) {

                it = featSetNearest.iterator();

                while (it.hasNext()) {

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

                    if (featureSelection != null && !featureSelection.isEmpty()

                        && !featureSelection.isSelected(feat))

                        continue;

                    List geomList = feat.getGeometries();

                    if (geomList == null) {

                        org.gvsig.fmap.geom.Geometry g2 =

                            feat.getDefaultGeometry();

                        Geometry g2Jts = GeometryUtil.geomToJTS(g2);

                        double dist = gJts.distance(g2Jts);

                        if (dist <= nearestDistance) {

                            buildFeature(featureInput, feat, new Double(

                                nearestDistance), g);

                            addedFeature = true;

                        }

                        continue;

                    }

                    Iterator<org.gvsig.fmap.geom.Geometry> itGeom =

                        geomList.iterator();

                    while (itGeom.hasNext()) {

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

                        Geometry g2Jts = GeometryUtil.geomToJTS(g2);

                        double dist = gJts.distance(g2Jts);

                        if (dist <= nearestDistance && !addedFeature) {

                            buildFeature(featureInput, feat, new Double(

                                nearestDistance), g);

                            addedFeature = true;

                        }

                    }

                }

                it.dispose();

            }

		} catch(FeatureIndexException e) {

			Sextante.addErrorToLog(e);

		} catch (DataException 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) {

	}

}

/**

 * 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 org.gvsig.fmap.dal.DALLocator;

import org.gvsig.fmap.dal.DataManager;

import org.gvsig.fmap.dal.DataTypes;

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

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

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

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

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

/**

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

		IVectorLayer layer2 = m_Parameters.getParameterValueAsVectorLayer(LAYER2);

		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 lyr1 = null;

		FlyrVectIVectorLayer lyr2 = null;

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

			lyr2 = ((FlyrVectIVectorLayer)layer2);

			lyr1 = ((FlyrVectIVectorLayer)layer1);

		} else 

			return false;

		DataManager dataManager = DALLocator.getDataManager();

		//Builds a new JSIRTree index. To do that we have to replace the default index and 

		//when the operation has finished then it restores the original default index 

		String indexName = "GEOMETRY";

		try {

			indexName = lyr2.getFeatureStore().getDefaultFeatureType().getDefaultGeometryAttributeName();

			FeatureAttributeDescriptor fat = lyr2.getFeatureStore().getDefaultFeatureType().getAttributeDescriptor(indexName);

			String defaultIndex = dataManager.getDefaultFeatureIndexProviderName(

					fat.getDataType().getType());

			dataManager.setDefaultFeatureIndexProviderName(

					fat.getDataType().getType(), "JSIRSpatialIndexProvider");

			lyr2.getFeatureStore().createIndex(

					lyr2.getFeatureStore().getDefaultFeatureType(), indexName, indexName + "_idx");

			dataManager.setDefaultFeatureIndexProviderName(

					fat.getDataType().getType(), defaultIndex);

		} catch (DataException e) {

			Sextante.addErrorToLog(e);

		}

		//Builds the output and computes the operation

		try {

			FeatureType featureType1 = lyr1.getFeatureStore().getDefaultFeatureType();

			FeatureType featureType2 = lyr2.getFeatureStore().getDefaultFeatureType();

			GeometryOperation operation = null;

			FeatureStore outFeatStore = null;

            if (nearest) {

                outFeatStore =

                    buildOutPutStoreFromUnion(featureType1, featureType2,

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

                        RESULT, "DIST", Double.class);

                operation =

                    new SpatiallyIndexedSpatialJoinOperation(lyr2, indexName

                        + "_idx", this);

            } else {

                outFeatStore =

                    buildSpatialJoinOutPutStore(featureType1,

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

                        RESULT);

                IDissolveRule rule = new DissolveRule(0, funcMap);

                Summary summary =

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

                operation =

                    new IntersectsSpatialJoinOperation(lyr2,

                        indexName + "_idx", summary, this);

            }

            operation.setTaskStatus(getStatus());

			operation.computesGeometryOperation(lyr1.getFeatureStore(), 

					outFeatStore, 

					attrNames, 

					selectedGeomInput, 

					selectedGeomOverlay, 

					true);		

		} catch (DataException e) {

			Sextante.addErrorToLog(e);

		}

		return true;

	}

	/**

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

    }

}

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.

#

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' ?>

<?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.33</version>

	</parent>

	<dependencies>

		<dependency>

		    <groupId>org.gvsig</groupId>

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