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.slf4j.Logger;

import org.slf4j.LoggerFactory;

import com.vividsolutions.jts.geom.Geometry;

import java.util.Iterator;

/**

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

        Iterator 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.getCopy();

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

        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 {

                summary.updateValues(elem.feat);

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

            }

        }

    }

}