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

import java.util.ArrayList;

import java.util.List;

import java.util.Stack;

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

import org.gvsig.fmap.geom.GeometryLocator;

import org.gvsig.fmap.geom.GeometryManager;

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

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

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

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

import com.vividsolutions.jts.geom.Geometry;

/**

 * Fuse spatially operation

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

 */

public class FuseSpatiallyOperation extends GeometryOperation {

        protected GeometryManager                geomManager           = GeometryLocator.getGeometryManager();

        private int                              id                    = 0;

        private List<Feature>                    featList              = null;

        private FeatureStore                     outFeatureStoreTable  = null;

        private String                           idField               = "FID";

        public FuseSpatiallyOperation(List<Feature> f, FeatureStore outFeatStoreTable, String[] fieldNames, String idField) throws DataException {

                this.featList = f;

                this.outFeatureStoreTable = outFeatStoreTable;

                this.idField = idField;

        }

        public FuseSpatiallyOperation(Stack<Feature> f, FeatureStore outFeatStoreTable, String[] fieldNames) throws DataException {

                this.featList = f;

                this.outFeatureStoreTable = outFeatStoreTable;

        }

        /*

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

                try {

                        Geometry jtsInputGeom = GeometryUtil.geomToJTS(g);

                        insertInTable(feature);

                        for (int i = featList.size() - 1; i >= 0; i--) {

                                Feature f = featList.get(i);

                                g = f.getDefaultGeometry();

                                Geometry jtsGeom = GeometryUtil.geomToJTS(g); //Multigeometry support

                                if(jtsGeom.intersects(jtsInputGeom)) {

                                        jtsInputGeom = jtsInputGeom.union(jtsGeom);

                                        featList.remove(i);

                                        insertInTable(f);

                                }

                        }

                        /*int i = featList.size() - 1;

                        while(i >= 0) {

                                Feature f = featList.get(i);

                                geom = f.getDefaultGeometry();

                                Geometry jtsGeom = GeometryUtil.geomToJTS(geom); //Multigeometry support

                                if(jtsInputGeom.intersects(jtsGeom)) {

                                        jtsInputGeom = jtsInputGeom.union(jtsGeom);

                                        featList.remove(i);

                                        insertInTable(f);

                                        i = featList.size();

                                }

                                i --;

                        }*/

                        //Cuando no genera dos capas mete los campos de la feature de entrada en la de salida perdiendose el resto

                        if(outFeatureStoreTable == null) { 

                                String[] fieldNames = persister.getFieldNamesWithoutGeom();

                                ArrayList<String> fields = new ArrayList<String>();

                                ArrayList<Object> values = new ArrayList<Object>();

                                fields.add(idField);

                                values.add(id);

                                for (int j = 0; j < fieldNames.length; j++) {

                                        Object obj = feature.get(fieldNames[j]);

                                        if(obj != null && fieldNames[j].compareTo(idField) != 0) {

                                                fields.add(fieldNames[j]);

                                                values.add(obj);

                                        }

                                }

                                lastEditFeature = persister.addFeature(jtsInputGeom, fields, values);

                        } else

                                lastEditFeature = persister.addFeature(jtsInputGeom, idField, id);

                } catch (DataException e) {

                        return null;

                } catch (CreateGeometryException e) {

                        return null;

                }

                id++;

                return lastEditFeature;

        }

        private void insertInTable(Feature f) throws DataException {

                if(outFeatureStoreTable == null)

                        return;

                EditableFeature edFeat = outFeatureStoreTable.createNewFeature();

                edFeat.set(idField, id);

                FeatureAttributeDescriptor[] attrList = f.getType().getAttributeDescriptors();

                for (int j = 0; j < attrList.length; j++) {

                        if(attrList[j].getName().compareTo("GEOMETRY") != 0) {

                                edFeat.set(attrList[j].getName(), f.get(attrList[j].getName()));

                        }

                }

                outFeatureStoreTable.insert(edFeat);

        }

        /**

         * Removes duplicate fields

         * @param names

         * @return

         */

        public static String[] checkFields(String[] names) {

                if(names.length <= 1)

                        return names;

                int cont = 0;

                int i = 1;

                while(i < names.length) {

                        if(names[0].compareTo(names[i]) == 0) {

                                names[0] = "FID_" + cont;

                                i = 0;

                                cont ++;

                        }

                        i ++;

                }

                return names;

        }

        /*

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

        }

        public DALFeaturePersister getWriter() {

                return persister;

        }

}