Application: gvSIG desktop

/*

 * Created on 09-may-2006

 *

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

 *

 * Copyright (C) 2004 IVER T.I. and 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.

 *

 * For more information, contact:

 *

 *  Generalitat Valenciana

 *   Conselleria d'Infraestructures i Transport

 *   Av. Blasco Ib??ez, 50

 *   46010 VALENCIA

 *   SPAIN

 *

 *      +34 963862235

 *   gvsig@gva.es

 *      www.gvsig.gva.es

 *

 *    or

 *

 *   IVER T.I. S.A

 *   Salamanca 50

 *   46005 Valencia

 *   Spain

 *

 *   +34 963163400

 *   dac@iver.es

 */

/* CVS MESSAGES:

 *

 * $Id: AdjacencyDissolveVisitor.java 10626 2007-03-06 16:55:54Z caballero $

 * $Log$

 * Revision 1.3  2007-03-06 16:47:58  caballero

 * Exceptions

 *

 * Revision 1.2  2006/07/21 11:06:06  azabala

 * fixed bug 604: empty dist field in buffered dissolved features

 *

 * Revision 1.1  2006/06/20 18:20:45  azabala

 * first version in cvs

 *

 * Revision 1.1  2006/05/24 21:11:14  azabala

 * primera version en cvs despues de refactoring orientado a crear un framework extensible de geoprocessing

 *

 *

 */

package com.iver.cit.gvsig.geoprocess.impl.dissolve.fmap;

import java.awt.geom.Rectangle2D;

import java.util.ArrayList;

import java.util.List;

import java.util.Stack;

import com.hardcode.gdbms.driver.exceptions.ReadDriverException;

import com.hardcode.gdbms.engine.data.driver.DriverException;

import com.hardcode.gdbms.engine.values.DoubleValue;

import com.hardcode.gdbms.engine.values.Value;

import com.hardcode.gdbms.engine.values.ValueFactory;

import com.iver.cit.gvsig.exceptions.expansionfile.ExpansionFileReadException;

import com.iver.cit.gvsig.exceptions.visitors.ProcessVisitorException;

import com.iver.cit.gvsig.exceptions.visitors.VisitorException;

import com.iver.cit.gvsig.fmap.core.IGeometry;

import com.iver.cit.gvsig.geoprocess.core.fmap.FeatureProcessor;

import com.iver.cit.gvsig.geoprocess.core.fmap.XTypes;

import com.vividsolutions.jts.geom.Geometry;

/**

 * <p>

 * This Visitor generates a dissolve of input layer based in adjacency criteria:

 * it dissolves two polygons, if and only if they are adjacent.

 * </p>

 * TODO: To dissolve buffers from radial rings create a DissolveVisitor

 * to compare FROM-TO Fields.

 *

 * @author azabala

 *

 */

public class AdjacencyDissolveVisitor extends DissolveVisitor {

        /**

         * Reference to the buffer distance of the visited buffered feature.

         */

        private DoubleValue currentBufferDistance = null;

        public AdjacencyDissolveVisitor(String dissolveField,

                        FeatureProcessor processor) {

                super(dissolveField, processor);

        }

        protected boolean verifyIfDissolve(DissolvedFeature f1, DissolvedFeature f2) {

                Geometry geo1 = f1.getJtsGeometry();

                Geometry geo2 = f2.getJtsGeometry();

                return geo1.intersects(geo2);

        }

        /**

         * Creates a new IFeature with util info for dissolve geoprocess

         * (it ignore non numerical values, etc)

         *

         * @param g

         * @param index

         * @return

         * @throws DriverException

         */

        protected DissolvedFeature createFeature(IGeometry g, int index) {

                DissolvedFeature solution = new DissolvedFeature(g, null, index);

                return solution;

        }

        /**

         * We overwrite this method because we are not interested in save sumarization

         * function values with dissolved features.

         * Instead, we want to add buffer distance like an attribute of them.

         */

        public void visit(IGeometry g, int index) throws VisitorException, ProcessVisitorException {

                if(g == null)

                        return;

                if(g.getGeometryType() != XTypes.POLYGON &&

                                g.getGeometryType() != XTypes.MULTI)

                        return;

                if (!dissolvedGeometries.get(index)) {

                        try {

                                int fieldIndex = recordset.getFieldIndexByName("DIST");

                                currentBufferDistance = (DoubleValue) recordset.getFieldValue(index, fieldIndex);

                                // if we havent dissolved this feature

                                Stack toDissol = new Stack();// stack for adjacent features

                                DissolvedFeature feature = createFeature(g, index);

                                toDissol.push(feature);

                                ArrayList geometries = new ArrayList();

                                Value[] values = dissolveGeometries(toDissol, geometries);

                                Geometry geometry = union(geometries);

                                Value[] valuesWithFID = new Value[values.length + 1];

                                System.arraycopy(values, 0, valuesWithFID, 1, values.length);

                                valuesWithFID[0] = ValueFactory.createValue(fid);

                                DissolvedFeature dissolved = new DissolvedFeature(null,valuesWithFID, fid/*index*/);

                                dissolved.setJtsGeometry(geometry);

                                this.featureProcessor.processFeature(dissolved);

                                fid++;

                                resetFunctions();

                        } catch (ReadDriverException e) {

                                throw new ProcessVisitorException(recordset.getName(),e,

                                        "Error al procesar las geometrias a fusionar durante dissolve");

                        } catch (ExpansionFileReadException e) {

                                throw new ProcessVisitorException(recordset.getName(),e,

                                "Error al procesar las geometrias a fusionar durante dissolve");

                        } catch (VisitorException e) {

                                throw new ProcessVisitorException(recordset.getName(),e,

                                "Error al procesar las geometrias a fusionar durante dissolve");

                        }

                }// if

        }

        /**

         * We overwrite this method to ignore sumarization values and to

         * add buffer distance to the attributes of the result features.

         * @throws VisitorException

         * @throws ExpansionFileReadException

         * @throws ReadDriverException

         */

        protected Value[] dissolveGeometries(Stack toDissol,

                        List geometries) throws

                        ReadDriverException, ExpansionFileReadException, VisitorException {

                IndividualGeometryDissolveVisitor visitor = null;

                DissolvedFeature feature = null;

                while (toDissol.size() != 0) {

                        feature = (DissolvedFeature) toDissol.pop();

                        // flags this idx (to not to process in future)

                        dissolvedGeometries.set(feature.getIndex());

                        if (visitor == null) {

                                visitor = new IndividualGeometryDissolveVisitor(feature,

                                                dissolvedGeometries, toDissol,

                                                numericField_sumarizeFunction);

                                visitor.setDissolveField(this.dissolveField);

                        } else {

                                visitor.setProcessedFeature(feature);

                        }

                        Rectangle2D bounds = feature.getGeometry().getBounds2D();

                        double xmin = bounds.getMinX();

                        double ymin = bounds.getMinY();

                        double xmax = bounds.getMaxX();

                        double ymax = bounds.getMaxY();

                        double magnify = 15d;

                        Rectangle2D query = new Rectangle2D.Double(xmin - magnify, ymin

                                        - magnify, (xmax - xmin) + magnify, (ymax - ymin) + magnify);

                        strategy.process(visitor, query);

                        //al final de toda la pila de llamadas recursivas,

                        //geometries tendr? todas las geometrias que debemos dissolver

                        geometries.add(feature.getJtsGeometry());

                }// while

                Value[] values = new Value[1];

                values[0] = currentBufferDistance;

                return values;

        }

}