Application: gvSIG desktop » gvSIG geoprocess

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

 *

 * Copyright (C) 2010 Generalitat Valenciana.

 *

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

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

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

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

 *

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

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

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 * GNU General Public License for more details.

 *

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

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

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

 */

package org.gvsig.geoprocess.algorithm.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.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.exception.DataException;

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

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

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

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

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

import org.gvsig.geoprocess.sextante.AbstractSextanteGeoProcess;

/**

 * 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        SELECTED_GEOM     = "SELECTED_GEOM";

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

        /*

         * (non-Javadoc)

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

         */

        public void defineCharacteristics(){

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

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

        // setGeneratesUserDefinedRasterOutput(false);

                try {

                        m_Parameters.addInputVectorLayer(LAYER1, 

                                                                                                Sextante.getText("Input_layer"), 

                                                                                                IVectorLayer.SHAPE_TYPE_WRONG, 

                                                                                                true);

                        m_Parameters.addInputVectorLayer(LAYER2, 

                                                                                                Sextante.getText( "Input_layer"), 

                                                                                                IVectorLayer.SHAPE_TYPE_WRONG, 

                                                                                                true);

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

                        m_Parameters.addBoolean(NEAREST, Sextante.getText("use_the_nearest"), false);

                        m_Parameters.addString(FUNCTION_LIST, Sextante.getText("Function_list"));

                } catch (RepeatedParameterNameException e) {

                        Sextante.addErrorToLog(e);

                }

                addOutputVectorLayer(RESULT,

                                                                Sextante.getText("Spatialjoin"),

                                                                OutputVectorLayer.SHAPE_TYPE_UNDEFINED);

        }

        /*

         * (non-Javadoc)

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

         */

        public boolean processAlgorithm() throws GeoAlgorithmExecutionException {

                IVectorLayer layer1 = m_Parameters.getParameterValueAsVectorLayer(LAYER1);

                IVectorLayer layer2 = m_Parameters.getParameterValueAsVectorLayer(LAYER2);

                boolean selectedGeom = m_Parameters.getParameterValueAsBoolean(SELECTED_GEOM);

                boolean nearest = m_Parameters.getParameterValueAsBoolean(NEAREST);

                String functionList = m_Parameters.getParameterValueAsString(FUNCTION_LIST);

                loadSummary(functionList);

                gvVectorLayer lyr1 = null;

                gvVectorLayer lyr2 = null;

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

                        lyr2 = ((gvVectorLayer)layer2);

                        lyr1 = ((gvVectorLayer)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(), "JSIRTree");

                        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 {

                        FeatureSet features = null;

                        features = lyr1.getFeatureStore().getFeatureSet();

                        FeatureType featureType1 = features.getDefaultFeatureType();

                        features = lyr2.getFeatureStore().getFeatureSet();

                        FeatureType featureType2 = features.getDefaultFeatureType();

                        GeometryOperation operation = null;

                        FeatureStore outFeatStore = null;

                        if(nearest) {

                                outFeatStore = buildOutPutStoreFromUnion(featureType1, 

                                                featureType2, 

                                                lyr1.getShapeType(), 

                                                Sextante.getText("SpatialJoin"), 

                                                RESULT, 

                                                "DIST", 

                                                Double.class);

                                operation = new SpatiallyIndexedSpatialJoinOperation(lyr2, indexName + "_idx");

                        } else {

                                outFeatStore = buildSpatialJoinOutPutStore(featureType1,

                                                lyr1.getShapeType(), 

                                                Sextante.getText("SpatialJoin"), 

                                                RESULT);

                                IDissolveRule rule = new DissolveRule(0, funcMap);

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

                                operation = new IntersectsSpatialJoinOperation(lyr2, indexName + "_idx", summary);

                        }

                        operation.setProgressModel(this);

                        operation.computesGeometryOperation(lyr1.getFeatureStore(), 

                                        outFeatStore, 

                                        attrNames, 

                                        selectedGeom, 

                                        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

         */

        @SuppressWarnings("unchecked")

        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 ((gvVectorLayer)output).getFeatureStore();

                } catch (UnsupportedOutputChannelException e) {

                        Sextante.addErrorToLog(e);

        } catch (GeoAlgorithmExecutionException e) {

            Sextante.addErrorToLog(e);

        }

                return null;

        }

}