org.matsim.contrib.socnetgen.sna.graph.spatial.analysis.EdgeLengthMedian.java Source code

Introduction

Here is the source code for org.matsim.contrib.socnetgen.sna.graph.spatial.analysis.EdgeLengthMedian.java

Source

/* *********************************************************************** *
 * project: org.matsim.*
 * EdgeLengthMedian.java
 *                                                                         *
 * *********************************************************************** *
 *                                                                         *
 * copyright       : (C) 2011 by the members listed in the COPYING,        *
 *                   LICENSE and WARRANTY file.                            *
 * email           : info at matsim dot org                                *
 *                                                                         *
 * *********************************************************************** *
 *                                                                         *
 *   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.                                   *
 *   See also COPYING, LICENSE and WARRANTY file                           *
 *                                                                         *
 * *********************************************************************** */
package org.matsim.contrib.socnetgen.sna.graph.spatial.analysis;

import gnu.trove.map.hash.TObjectDoubleHashMap;

import java.util.Set;

import org.apache.commons.math.stat.descriptive.DescriptiveStatistics;
import org.matsim.contrib.common.gis.CRSUtils;
import org.matsim.contrib.socnetgen.sna.gis.DistanceCalculatorFactory;
import org.matsim.contrib.socnetgen.sna.graph.Vertex;
import org.matsim.contrib.socnetgen.sna.graph.spatial.SpatialVertex;

/**
 * @author illenberger
 *
 */
public class EdgeLengthMedian extends AbstractSpatialProperty {

    @Override
    public TObjectDoubleHashMap<Vertex> values(Set<? extends Vertex> vertices) {
        @SuppressWarnings("unchecked")
        Set<? extends SpatialVertex> spatialVertices = (Set<? extends SpatialVertex>) vertices;

        TObjectDoubleHashMap<Vertex> values = new TObjectDoubleHashMap<Vertex>();

        for (SpatialVertex vertex : spatialVertices) {
            if (vertex.getPoint() != null) {
                double val = edgeLengthMedian(vertex);
                if (!Double.isNaN(val))
                    values.put(vertex, val);
            }
        }

        return values;
    }

    private double edgeLengthMedian(SpatialVertex vertex) {
        if (calculator == null) {
            calculator = DistanceCalculatorFactory
                    .createDistanceCalculator(CRSUtils.getCRS(vertex.getPoint().getSRID()));
        }

        DescriptiveStatistics stats = new DescriptiveStatistics();

        for (SpatialVertex neighbor : vertex.getNeighbours()) {
            if (neighbor.getPoint() != null) {
                if (vertex.getPoint().getSRID() == neighbor.getPoint().getSRID()) {
                    stats.addValue(calculator.distance(vertex.getPoint(), neighbor.getPoint()));
                } else {
                    throw new RuntimeException("Points do not share the same coordinate reference system.");
                }
            }
        }

        if (stats.getN() > 1)
            return stats.getPercentile(50);
        else
            return Double.NaN;
    }
}