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

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

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

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

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

import com.vividsolutions.jts.geom.Geometry;

import es.unex.sextante.core.Sextante;

/**

 * Atomic dissolve operation

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

 */

public class AtomicDissolveOperation extends GeometryOperation {

        private IDissolveRule                    rule             = null;

        private ArrayList<Geometry>              geometryList     = null;

        private ArrayList<Feature>               featList         = null;

        private Feature                          feature          = null;

        private EditableFeature                  result           = null;

        private FeatureStore                     outFeatureStore  = null;

        private int                              newFeatID        = 0;

        private Summary                          summary          = null;                       

        public AtomicDissolveOperation(AbstractSextanteGeoProcess p) {

                super(p);

        }

        /**

         * Sets IDissolveCriteria. This object will choose what geometry

         * will be dissolved   

         * @param criteria

         */

        public void setCriteria(IDissolveRule rule) {

                this.rule = rule;

        }

        /**

         * Sets the list of geometries

         * @param geoList

         */

        public void setGeometryList(ArrayList<Geometry> geoList) {

                this.geometryList = geoList;

        }

        /**

         * Sets the list of geometries

         * @param geoList

         */

        public void setFeatureList(ArrayList<Feature> featList) {

                this.featList = featList;

        }

        /**

         * Sets the current feature

         * @param feat

         */

        public void setFeature(Feature feat) {

                this.feature = feat;

        }

        /**

         * Sets the output FeatureStore

         * @param out

         * @throws DataException 

         */

        public void setFeatureStore(FeatureStore out, String[] attrNames) {

                this.outFeatureStore = out;

        }

        /*

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

                Geometry g1 = geometryList.get(0);

                if(g1 == null || g instanceof NullGeometry) 

                        return null;

                com.vividsolutions.jts.geom.Geometry g2 = GeometryUtil.geomToJTS(g);

                //Si son iguales las geometrias no se opera

                if(g1 == null || g1.compareTo(g2) == 0)

                        return null;

                if(rule.verifyIfDissolve(g1, g2, this.feature, feature)) {

                        geometryList.add(g2);

                        Geometry newGeom = GeometryUtil.geometryUnion(geometryList, g.getGeometryType().getType());

                        try {

                                result = outFeatureStore.createNewFeature();

                                result.setDouble(0, newFeatID);

                                result.set(1, feature.get(rule.getFieldName()));

                                summary.updateValues(feature);

                                summary.loadEditableFeature(result);

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

                                geometryList.clear();

                                geometryList.add(newGeom);

                                return result;

                        } catch (DataException e) {

                                Sextante.addErrorToLog(e);

                        }

                }

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

        }

        /*

         * (non-Javadoc)

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

         */

        public Object getResult() {

                newFeatID ++;

                return result;

        }

        /**

         * Computes this operation without build an output layer. When the operation

         * finish the result is got with getResult method.

         * @param selectedGeom

         * @throws DataException

         */

        public void computesOperation(boolean selectedGeom) throws DataException {

                result = null;

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

                summary.loadDefaultSummarizes(this.feature);

                int i = featList.size() - 1;

                while(i >= 0) {

                        if(process.getTaskMonitor().isCanceled()) 

                                return;

                        Feature feature = featList.get(i);

                        org.gvsig.fmap.geom.Geometry geom = feature.getDefaultGeometry();

                        EditableFeature ef = invoke(geom, feature);

                        if(ef != null) {

                                featList.remove(i);

                                if (process != null) 

                                        process.setProgress(numberOfFeatures - featList.size(), numberOfFeatures);

                                if(rule instanceof AdjacencyDissolveRule)

                                        i = featList.size();

                        }

                        i--;

                }

                //Si el resultado es nulo es que no se ha hecho un dissolve porque no hab?a geometrias que unir

                //por lo que el resultado ser?n las geometrias de entrada

                if(result == null) {

                        try {

                                Geometry newGeom = GeometryUtil.geometryUnion(geometryList, feature.getDefaultGeometry().getGeometryType().getType());

                                result = outFeatureStore.createNewFeature();

                                result.setDouble(0, newFeatID);

                                result.set(1, feature.get(rule.getFieldName()));

                                summary.updateValues(feature);

                                summary.loadEditableFeature(result);

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

                        } catch (DataException e) {

                                Sextante.addErrorToLog(e);

                        }

                }

        }

        /**

         * Computes this operation without build an output layer. When the operation

         * finish the result is got with getResult method.

         * @param selectedGeom

         * @throws DataException

         */

//        public void computesOperation(boolean selectedGeom) throws DataException {

//                FeatureSet featuresSet = null;

//                featuresSet = inFeatureStore.getFeatureSet();

//                result = null;

//                DisposableIterator it = null;

//                

//                summary = new Summary(rule, outFeatureStore.getDefaultFeatureType());

//                summary.loadDefaultSummarizes(this.feature);

//

//                if(selectedGeom) {

//            FeatureSelection ds = inFeatureStore.getFeatureSelection();

//            it = ds.iterator();

//        } else {

//                        it = featuresSet.iterator();

//        }

//                

//                featAnalized = 0;

//                while( it.hasNext() ) {

//                        Feature feature = (Feature)it.next();

//                        if(analizedFeats[featAnalized]) {

//                                featAnalized ++;

//                                continue;

//                        }

//                        

//                        List geomList = feature.getGeometries();

//                        

//                        if(geomList == null) {

//                                org.gvsig.fmap.geom.Geometry geom = feature.getDefaultGeometry();

//                                invoke(geom, feature);

//                                featAnalized ++;

//                                continue;

//                        }

//

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

//                        while(itGeom.hasNext()) {

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

//                                invoke(g, feature);

//                        }

//                        featAnalized ++;

//                }

//                

//                //Si el resultado es nulo es que no se ha hecho un dissolve porque no hab?a geometrias que unir

//                //por lo que el resultado ser?n las geometrias de entrada

//                if(result == null) {

//                        try {

//                                Geometry newGeom = GeometryUtil.geometryUnion(geometryList, feature.getDefaultGeometry().getGeometryType().getType());

//                                result = outFeatureStore.createNewFeature();

//                                result.setDouble(0, newFeatID);

//                                result.set(1, feature.get(rule.getIndexField()));

//                                summary.updateValues(feature);

//                                summary.loadEditableFeature(result);

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

//                        } catch (DataException e) {

//                                Sextante.addErrorToLog(e);

//                        }

//                }

//        }

}