Application: gvSIG desktop

/* gvSIG. Geographic Information System of the Valencian Government

 *

 * Copyright (C) 2007-2008 Infrastructures and Transports Department

 * of the Valencian Government (CIT)

 *

 * 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.

 *

 */

/*

 * AUTHORS (In addition to CIT):

 * 2008 Prodevelop S.L  main development

 */

package org.gvsig.geocoding.geommatches;

import java.awt.geom.PathIterator;

import java.util.ArrayList;

import java.util.HashMap;

import java.util.Iterator;

import java.util.List;

import java.util.Map;

import org.gvsig.fmap.dal.feature.Feature;

import org.gvsig.fmap.geom.Geometry;

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

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

import org.gvsig.fmap.geom.util.UtilFunctions;

import com.vividsolutions.jts.geom.Coordinate;

import com.vividsolutions.jts.linearref.LengthIndexedLine;

/**

 *

 * This class searchs positions over lines (ways)

 *

 * @author <a href="mailto:jsanz@prodevelop.es"> Jorge Gaspar Sanz Salinas</a>

 * @author <a href="mailto:vsanjaime@prodevelop.es"> Vicente Sanjaime Calvet</a>

 *

 */

public class SimpleLineMatcher extends GeomMatcher {

        /**

         * This method gets the position over lines (Streets) or next to the lines

         *

         * @return

         */

        @SuppressWarnings("unchecked")

        public Geometry[] match() {

                // Group the geometries by one field

                HashMap geomstreets1 = groupFeaturesByAttribute(this

                                .getGroupFieldPosition(), (List) this

                                .getFeatures().get(0));

                Iterator it = geomstreets1.entrySet().iterator();

                Point2D[] points = new Point2D[geomstreets1.size()];

                int position = 0;

                while (it.hasNext()) {

                        Map.Entry e = (Map.Entry) it.next();

                        ArrayList geomss = (ArrayList) e.getValue();

                        // Insert geometries of one street into array

                        Geometry[] arrGeoms = new Geometry[geomss.size()];

                        for (int i = 0; i < geomss.size(); i++) {

                                arrGeoms[i] = (Geometry) geomss.get(i);

                        }

                        // Parse geometries to JTS

                        com.vividsolutions.jts.geom.Geometry[] geomsjts = parseGeomsToJTS(

                                        geomss.get(0).getClass(), arrGeoms);

                        // Make the UNION of the all geometries with the same attribute

                        com.vividsolutions.jts.geom.Geometry geomstreetjts = unionLinesJTS(geomsjts);

                        // Get the point of the geometry (LINE)

                        Point2D point = getLinesPosition(geomstreetjts, getParams());

                        points[position] = point;

                        position++;

                }

                return points;

        }

        /**

         * Calculate the position inside sigle line

         *

         * @param geometryJts

         *            line of street

         * @param params

         *            Parameters to position the point over the lines

         * @return

         */

        private Point2D getLinesPosition(

                        com.vividsolutions.jts.geom.Geometry geometryJts,

                        MatchesParams params) {

                double totaldistance = geometryJts.getLength();

                LengthIndexedLine lenline = new LengthIndexedLine(geometryJts);

                Long distance = params.getDistance();

                Integer relativedist = params.getRelativeDistance();

                Double offset = params.getOffset();

                Coordinate coors = null;

                // first calculate with the distance

                if (distance != null) {

                        if (lenline.isValidIndex(distance)) {

                                coors = lenline.extractPoint(distance);

                                // with offset

                                if (offset != null && offset.doubleValue() > 0) {

                                        return getOffsetPosition(new Point2D(coors.x, coors.y),

                                                        offset);

                                }

                                // Without offset

                                else {

                                        return new Point2D(coors.x, coors.y);

                                }

                        }

                }

                // Calculate percentage of distance

                else if (relativedist != null) {

                        if (totaldistance > 0) {

                                double dist = (relativedist * totaldistance) / 100.0;

                                coors = lenline.extractPoint(dist);

                                // with offset

                                if (offset != null && offset.doubleValue() > 0) {

                                        return getOffsetPosition(new Point2D(coors.x, coors.y),

                                                        offset);

                                }

                                // Without offset

                                else {

                                        return new Point2D(coors.x, coors.y);

                                }

                        }

                }

                // Any other case

                else {

                        if (totaldistance > 0) {

                                double dist = totaldistance / 2.0;

                                coors = lenline.extractPoint(dist);

                                // with offset

                                if (offset != null && offset.doubleValue() > 0) {

                                        return getOffsetPosition(new Point2D(coors.x, coors.y),

                                                        offset);

                                }

                                // Without offset

                                else {

                                        return new Point2D(coors.x, coors.y);

                                }

                        }

                }

                return null;

        }

        /**

         * Calculate the position in the line with offset

         *

         * @param linePoint

         *            middle point over the liner

         * @param dist

         *            perpendicular distance since the line

         * @return

         */

        @SuppressWarnings("unchecked")

        private Point2D getOffsetPosition(Point2D linePoint, double dist) {

                List feats = (List) getFeatures().get(0);

                Geometry[] geoms = new Geometry[feats.size()];

                Iterator it = feats.iterator();

                int cont = 0;

                while (it.hasNext()) {

                        geoms[cont] = (Geometry) ((Feature) it.next())

                                        .getGeometry("GEOMETRY");

                        cont++;

                }

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

                        double[] coords = new double[6];

                        Point2D inicio = null;

                        Point2D fin = null;

                        boolean inter = geoms[j].intersects(linePoint.getBounds2D());

                        if (inter) {

                                PathIterator iter = ((Curve2D) geoms[j]).getGeneralPathX()

                                                .getPathIterator(null);

                                iter.currentSegment(coords);

                                inicio = new Point2D(coords[0], coords[1]);

                                iter.next();

                                iter.currentSegment(coords);

                                fin = new Point2D(coords[0], coords[1]);

                                log.debug("Return point with offset");

                                return getPerpendicularPoint(inicio, fin, linePoint, dist);

                        }

                }

                log.debug("Return original point without offset");

                return linePoint;

        }

        /**

         *

         * This method calculates a point perpendicular to the line at a distance

         *

         * @param inicio

         *            initial point of the line

         * @param fin

         *            end point of the line

         * @param linePoint

         *            middle point over line

         * @param dist

         *            perpendicular distance to the line

         * @return

         */

        private Point2D getPerpendicularPoint(Point2D inicio, Point2D fin,

                        Point2D linePoint, double dist) {

                java.awt.geom.Point2D p1 = new java.awt.geom.Point2D.Double(inicio

                                .getX(), inicio.getY());

                java.awt.geom.Point2D p2 = new java.awt.geom.Point2D.Double(fin.getX(),

                                fin.getY());

                java.awt.geom.Point2D perpPoint = new java.awt.geom.Point2D.Double(

                                linePoint.getX(), linePoint.getY());

                java.awt.geom.Point2D[] p = UtilFunctions.getPerpendicular(p1, p2,

                                perpPoint);

                java.awt.geom.Point2D unit = UtilFunctions.getUnitVector(p[0], p[1]);

                java.awt.geom.Point2D res = new java.awt.geom.Point2D.Double(perpPoint

                                .getX()

                                + (unit.getX() * dist), perpPoint.getY() + (unit.getY() * dist));

                return new Point2D(res.getX(), res.getY());

        }

}