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.

 */

/*

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

import java.util.ArrayList;

import java.util.List;

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

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

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

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

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

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

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

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

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

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

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

import org.gvsig.tools.dispose.DisposableIterator;

import org.slf4j.Logger;

import org.slf4j.LoggerFactory;

import com.vividsolutions.jts.geom.Geometry;

/**

 * Dissolve operation

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

 */

public class DissolveOperationFast extends GeometryOperation {

        private static Logger                                  logger                      = LoggerFactory.getLogger(DissolveOperationFast.class.getName());

        private EditableFeature                  lastEditFeature  = null;

        private ArrayList<Element>               featureList      = null;       

        private IDissolveRule                    rule             = null;

        private Summary                          summary          = null; 

        class Element {

                public int                 id              = -1;

                public Feature             feat            = null;

                public List<Element>       overlapList     = new ArrayList<Element>();

                public boolean             insertedToJoin  = false;

                public Geometry            jtsGeom         = null;

        }

        class NodeTree {

                public Element       element     = null;

                public int           pos         = 0;

                public NodeTree       parent      = null;

                public NodeTree(Element node, NodeTree parent) {

                        this.element = node;

                        this.parent = parent;

                }

                public Element getNext() {

                        if(pos < element.overlapList.size())

                                return element.overlapList.get(pos++);

                        return null;

                }

        }

        public DissolveOperationFast(IDissolveRule rule, AbstractSextanteGeoProcess p) {

                super(p);

                this.rule = rule;

                featureList = new ArrayList<Element>();

        }

        /*

         * (non-Javadoc)

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

         */

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

                return null;

        }

        /*

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

                if(g == null)

                        return;

        }

        /*

         * (non-Javadoc)

         * @see org.gvsig.geoprocess.algorithm.base.core.IOperation#getResult()

         */

        public Object getResult() {

                return lastEditFeature;

        }

        /**

         * Computes a complete operation over the input FeatureStore. The result of this operation

         * is stored in the output FeatureStore. 

         * @param inFeatStore

         *        Input FeatureStore

         * @param outFeatStore

         *        Output FeatureStore

         * @param attrNames

         *        List of attributes to build the output feature store

         * @param selectedGeom

         *        If it is true only the selected geometries will be processed

         * @throws DataException

         */

        @SuppressWarnings({ "deprecation" })

        public void computesGeometryOperation(FeatureStore inFeatStore,

                                                                        FeatureStore outFeatStore,

                                                                        String[] attrNames,

                                                                        boolean selectedGeomInput,

                                                                        boolean selectedGeomOutput,

                                                                        boolean closeOutStore) throws DataException {

                this.inFeatureStore = inFeatStore;

                FeatureSet featuresSet = null;

                featuresSet = inFeatStore.getFeatureSet();

                setFeatureStore(outFeatStore, attrNames);

                DisposableIterator it = null;

                if(selectedGeomInput) {

            FeatureSelection ds = inFeatStore.getFeatureSelection();

            it = ds.iterator();

            numberOfFeatures = (int) ds.getSelectedCount();

                } else {

                        it = featuresSet.iterator();

                        numberOfFeatures = (int)featuresSet.getSize();

                }

        if (status != null && process != null) {

            status.setRangeOfValues(0, numberOfFeatures);

            process.setProgress(0, numberOfFeatures);

        }

        //Crear lista de elementos

                int iCount = 0;

                while( it.hasNext() && !process.getTaskMonitor().isCanceled()) {

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

                        Element el = new Element();

                        el.feat = feat;

                        el.id = iCount;

                        featureList.add(el);

            if (status != null && process != null) {

                status.setCurValue(iCount);

                process.setProgress(iCount, numberOfFeatures);

            }

                        iCount ++;

                }

                it.dispose();

                //Crear listas de solapes para cada feature

                iCount = 0;

                while (iCount < featureList.size() && !process.getTaskMonitor().isCanceled()) {

                        Element elem1 = featureList.get(iCount);

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

                        elem1.jtsGeom = GeometryUtil.geomToJTS(geom1);

                        if (status != null) 

                status.setCurValue(iCount);

            if(process != null) 

                process.setProgress(iCount, numberOfFeatures);

                        for (int i = iCount + 1; i < featureList.size(); i++) {

                                Element elem2 = featureList.get(i);

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

                                elem2.jtsGeom = GeometryUtil.geomToJTS(geom2);

                                if(rule.verifyIfDissolve(elem1.jtsGeom, elem2.jtsGeom, elem1.feat, elem2.feat)) {

                                        elem1.overlapList.add(elem2);

                                        elem2.overlapList.add(elem1);

                                }

                        }

                        iCount ++;

                }

                //Se calculan las listas de geometrias a unir

                //Para cada feature se obtiene su lista de elementos que solapan y para 

                //cada elemento que solapa se obtiene su lista. Finalmente todas se unen y 

                //y se hace un insert de una feature nueva

                List<Geometry> listResult = new ArrayList<Geometry>();

                iCount = 0;

                int iFeat = 0;

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

                while (iCount < featureList.size() && !process.getTaskMonitor().isCanceled()) {

                        Element elem1 = featureList.get(iCount);

                        summary.loadDefaultSummarizes(elem1.feat);

                        if (status != null) 

                status.setCurValue(iCount);

            if(process != null) 

                process.setProgress(iCount, numberOfFeatures);

            if(!elem1.insertedToJoin) {

                                buildListToJoin(elem1, iFeat);

                                iFeat ++;

                        }

                        /*if(!elem1.insertedToJoin) {

                                elem1.insertedToJoin = true;

                                listResult.clear();

                                //org.gvsig.fmap.geom.Geometry dalGeom = elem1.feat.getDefaultGeometry();

                                //Geometry jtsGeom = GeometryUtil.geomToJTS(dalGeom);

                                listResult.add(elem1.jtsGeom);

                                buildListToJoin(listResult, elem1.overlapList);

                                Geometry newGeom = computesUnion(listResult, elem1);//GeometryUtil.geometryUnion(listResult, elem1.feat.getDefaultGeometry().getGeometryType().getType());

                                try {

                                        addFeatureToOutput(newGeom, elem1.feat, iFeat);

                                } catch (CreateGeometryException e) {

                                        logger.info("Error a?adiendo geometr?a", e);

                                } catch (DataException e) {

                                        logger.info("Error a?adiendo geometr?a", e);

                                }

                                iFeat ++;

                        }*/

                        iCount ++;

                }

                if(closeOutStore && persister != null)

                        persister.end();

        }

        /**

         * Computes the union of the list of geometries

         * @param listResult

         * @param elem1

         * @return

         */

        private Geometry computesUnion(List<Geometry> listResult, Element elem1) {

                int splitValue = 500;

                Geometry newGeom = null;

                if(listResult.size() > splitValue) {

                        List<List<Geometry>> list = splitList(listResult, splitValue);

                        List<Geometry> result = new ArrayList<Geometry>();

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

                                Geometry aux = GeometryUtil.geometryUnion(list.get(i), elem1.feat.getDefaultGeometry().getGeometryType().getType());

                                result.add(aux);

                        }

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

                                newGeom = GeometryUtil.geometryUnion(result, elem1.feat.getDefaultGeometry().getGeometryType().getType());

                        }

                } else {

                        newGeom = GeometryUtil.geometryUnion(listResult, elem1.feat.getDefaultGeometry().getGeometryType().getType());        

                }

                return newGeom;

        }

        private Geometry computesUnion3(List<Geometry> listResult) {

                Geometry newGeom = null;

                int iCount = 0;

                int max = listResult.size();

                if(process != null)

                        process.setName("Generating union");

                while(listResult.size() > 1) {

                        List<Geometry> list = new ArrayList<Geometry>();

                        if (status != null)  

                status.setCurValue(iCount);

            if(process != null)

                process.setProgress(iCount, max);

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

                                if(i == (listResult.size() - 1))

                                        list.add(listResult.get(i));

                                else {

                                        newGeom = JTSFacade.union(listResult.get(i), listResult.get(i + 1));

                                        list.add(newGeom);

                                }

                        }

                        listResult = list;

                }

                return newGeom;

        }

        /**

         * Splits the array of geometries to compute its union because JTS cannot support

         * a lot of geometries

         * @param list

         * @param n

         * @return

         */

        private List<List<Geometry>> splitList(List<Geometry> list, int n) {

                int elements = (int)(list.size() / n);

                List<List<Geometry>> l = new ArrayList<List<Geometry>>();

                for (int i = 0; i < elements; i++) {

                        l.add(list.subList((i * n), (i * n) + n));

                }

                if(elements * n < list.size()) {

                        l.add(list.subList((elements * n), list.size()));

                }

                return l;

        }

        /**

         * Adds a feature to the output

         * @param newGeom

         * @param feat

         * @param newFeatID

         * @throws DataException

         * @throws CreateGeometryException

         */

        private void addFeatureToOutput(Geometry newGeom, Feature feat, int newFeatID) throws DataException, CreateGeometryException {

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

                result.setDouble(0, newFeatID);

                result.set(1, feat.get(rule.getIndexField()));

                result.setGeometry("GEOMETRY", GeometryUtil.jtsToGeom(newGeom));

                summary.loadEditableFeature(result);

                lastEditFeature = persister.addFeature(result, result.getDefaultGeometry());

        }

        /**

         * Builds the union of all lists 

         * @param listResult

         * @param listToJoin

         */

        private void buildListToJoin(List<Geometry> listResult, List<Element> listToJoin) {

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

                        Element elem = listToJoin.get(i);

                        if(!elem.insertedToJoin) {

                                elem.insertedToJoin = true;

                                buildListToJoin(listResult, elem.overlapList);

                                summary.updateValues(elem.feat);

                                //org.gvsig.fmap.geom.Geometry dalGeom = elem.feat.getDefaultGeometry();

                                //Geometry jtsGeom = GeometryUtil.geomToJTS(dalGeom);

                                listResult.add(elem.jtsGeom);

                        }

                }

        }

        /**

         * Builds the union of all lists 

         * @param listResult

         * @param listToJoin

         */

        private void buildListToJoin(Element elem, int iFeat) {

                Geometry newGeom = null;

                if(elem.overlapList.size() == 0) {

                        if(!elem.insertedToJoin)

                                elem.insertedToJoin = true;

                        try {

                                addFeatureToOutput(elem.jtsGeom, elem.feat, iFeat);

                        } catch (CreateGeometryException e) {

                                logger.info("Error a?adiendo geometr?a", e);

                        } catch (DataException e) {

                                logger.info("Error a?adiendo geometr?a", e);

                        }

                } else {

                        List<Geometry> listResult = new ArrayList<Geometry>();

                        NodeTree subtree = new NodeTree(elem, null);

                        //Hacemos un recorrido en profundidad del ?rbol para a?adir

                        //todos los elementos a la lista de geometrias a unir sin

                        //repetir ninguna.

                        while (subtree != null) {

                                if(!subtree.element.insertedToJoin) {

                                        listResult.add(subtree.element.jtsGeom);

                                        summary.updateValues(subtree.element.feat);

                                        subtree.element.insertedToJoin = true;

                                }

                                boolean back = false;

                                Element l = subtree.getNext();

                                if(l == null) 

                                        back = true;

                                while(!back && l.insertedToJoin) {

                                        l = subtree.getNext();

                                        if(l == null) 

                                                back = true;

                                }

                                if(back) {

                                        subtree = subtree.parent;

                                        continue;

                                }

                                subtree = new NodeTree(l, subtree);

                        }

                        newGeom = computesUnion3(listResult);

                        try {

                                addFeatureToOutput(newGeom, elem.feat, iFeat);

                        } catch (DataException e) {

                                logger.info("Imposible insertar en la tabla", e);

                        } catch (CreateGeometryException e) {

                                logger.info("Error a?adiendo geometr?a", e);

                        }

                }

        }

}