Application: gvSIG desktop » gvSIG geoprocess

/* gvSIG. Desktop Geographic Information System.

 *

 * Copyright ? 2007-2017 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.

 */

package org.gvsig.geoprocess.algorithm.dispersepoints;

import java.io.File;

import java.util.ArrayList;

import java.util.Collection;

import java.util.HashSet;

import java.util.Iterator;

import java.util.List;

import java.util.ListIterator;

import java.util.NavigableSet;

import java.util.Set;

import org.mapdb.DB;

import org.mapdb.DBMaker;

import es.unex.sextante.core.Sextante;

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

import org.gvsig.fmap.geom.Geometry;

import org.gvsig.fmap.geom.Geometry.SUBTYPES;

import org.gvsig.fmap.geom.GeometryLocator;

import org.gvsig.fmap.geom.GeometryManager;

import org.gvsig.fmap.geom.SpatialIndex;

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

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

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

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

import org.gvsig.tools.ToolsLocator;

import org.gvsig.tools.folders.FoldersManager;

/**

 * Disperse points operation

 *

 * @author fdiaz</a>

 */

public class DispersePointsOperation extends GeometryOperation {

    private SpatialIndex index;

    private double scatterRadius;

    private double matchDistance;

    NavigableSet<Integer> processed;

    /**

     * Builds an instance of this operation.

     *

     * @param index

     * @param scatterRadius

     * @param matchDistance

     * @param p

     *

     */

    public DispersePointsOperation(SpatialIndex index, double scatterRadius, double matchDistance,

        AbstractSextanteGeoProcess p) {

        super(p);

        this.index = index;

        this.scatterRadius = scatterRadius;

        this.matchDistance = matchDistance;

        FoldersManager folderManager = ToolsLocator.getFoldersManager();

        File dbfile = folderManager.getUniqueTemporaryFile("dispersepoints.mapdb");

        // Creariamos la BBDD

        DB db = DBMaker.fileDB(dbfile).fileMmapEnable().make();

        // Y dentro de ella una tabla "processed" para nuestro set.

        processed = (NavigableSet<Integer>) db.treeSet("processed").create();

    }

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

        boolean addedFeature = false;

        if (g == null)

            return lastEditFeature;

        try {

            FeatureReference reference = feature.getReference();

            if (!processed.contains(reference.hashCode())) {

                final List<FeatureReference> group = new ArrayList<FeatureReference>();

                Iterator<?> iterator = index.query(g.buffer(matchDistance));

                while (iterator.hasNext()) {

                    FeatureReference ref = (FeatureReference) iterator.next();

                    if (!processed.contains(ref.hashCode())) {

                        Feature feat = ref.getFeature().getCopy();

                        Geometry g2 = feat.getDefaultGeometry();

                        double dist = g.distance(g2);

                        if (dist <= matchDistance) {

                            group.add(ref);

                            processed.add(ref.hashCode());

                        }

                    }

                }

                if (group.size() == 1) {

                    buildFeature(feature, g);

                } else {

                    buildFeatures(group, (Point) g);

                }

            }

        } catch (Exception e) {

            Sextante.addErrorToLog(e);

        }

        return lastEditFeature;

    }

    /**

     * @param group

     * @throws DataException

     * @throws CreateGeometryException

     */

    private void buildFeatures(List<FeatureReference> group, Point center) throws CreateGeometryException, DataException {

        if (group.isEmpty()) {

            return;

        }

        GeometryManager geomManager = GeometryLocator.getGeometryManager();

        int points = group.size();

        double step = 2 * Math.PI / points;

        int i = 0;

        for (FeatureReference ref : group) {

            Feature feature = ref.getFeature().getCopy();

            buildFeature(feature,

                geomManager.createPoint(

                    center.getX() + (Math.cos(i * step) * this.scatterRadius),

                    center.getY() + (Math.sin(i * step) * this.scatterRadius),

                    SUBTYPES.GEOM2D));

            i++;

        }

    }

    /*

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

        invoke(g, (Feature) feature);

    }

    /**

     * Builds a feature and adds it to the output file

     *

     * @param feat

     * @param value

     * @param g

     * @throws DataException

     */

    private void buildFeature(Feature feat, org.gvsig.fmap.geom.Geometry g) throws DataException {

        EditableFeature outFeat = persister.getOutputFeatureStore().createNewFeature(feat);

        try {

            persister.addFeature(outFeat, g);

        } catch (CreateGeometryException e) {

            Sextante.addErrorToLog(e);

        }

    }

}