如何从OpenStreetMap数据中找到街道交叉口列表？

Question

我正在寻找一种从OpenStreetMap（OSM）数据中准确检索街道交叉口的方法。我知道有人问过和回答过类似的问题，但我可以从建议的方法中检索到的数据并不十分准确。

首先，我知道以下问题：

• 如何从OpenStreetMap查找十字路口？

上述问题的答案建议：

“查询给定边界框中的所有方式，并查找由两个或多个方式共享的节点，如另一个答案中所述。”

我遵循这个建议，编写了一个python脚本，从我从OpenStreetMap下载的xml文件（osm文件）中提取节点元素。代码如下：

try:
    from xml.etree import cElementTree as ET
except ImportError, e:
    from xml.etree import ElementTree as ET

def extract_intersections(osm, verbose=True):
    # This function takes an osm file as an input. It then goes through each xml 
    # element and searches for nodes that are shared by two or more ways.
    # Parameter:
    # - osm: An xml file that contains OpenStreetMap's map information
    # - verbose: If true, print some outputs to terminal.
    # 
    # Ex) extract_intersections('WashingtonDC.osm')
    #
    tree = ET.parse(osm)
    root = tree.getroot()
    counter = {}
    for child in root:
        if child.tag == 'way':
            for item in child:
                if item.tag == 'nd':
                    nd_ref = item.attrib['ref']
                    if not nd_ref in counter:
                        counter[nd_ref] = 0
                    counter[nd_ref] += 1

    # Find nodes that are shared with more than one way, which
    # might correspond to intersections
    intersections = filter(lambda x: counter[x] > 1,  counter)

    # Extract intersection coordinates
    # You can plot the result using this url.
    # http://www.darrinward.com/lat-long/
    intersection_coordinates = []
    for child in root:
        if child.tag == 'node' and child.attrib['id'] in intersections:
            coordinate = child.attrib['lat'] + ',' + child.attrib['lon']
            if verbose:
                print coordinate
            intersection_coordinates.append(coordinate)

    return intersection_coordinates

如果我使用从OSM导出的数据运行此代码（例如，我使用从出口区域导出的数据：Min Lat： 38.89239，Max Lat： 38.89981，Min Lon：-77.03212，Max Lon：-77.02119.)，它会打印出看起来喜欢：

38.8966440,-77.0259810
38.8973430,-77.0280900
38.9010391,-77.0270309
38.8961050,-77.0319620
...

（我用过http://www.darrinward.com/lat-long/ 显然，数据包含一些不是交点的节点（它们可能是朝向两个STEET的存储）

我做错了什么，还是这是我能从OSM得到的最好的“交叉”数据？我感谢你的帮助和意见。

最好的

Answer 1

第一个Tipp：

不要只与谷歌地图进行比较，主要是与OpenStreetMap可视化进行比较。特别是复杂的街道交叉口，尽管它们代表同一条道路，但可以建立不同的模型。

2） 当前位置如果你真的使用了正确的方式：那是步行道，混合了街道吗？有各种不同的类型，具有不同的网络属性：可供车辆访问，等等。在谷歌地图中，白色道路是指可供车辆访问的道路

3） 进一步看，如果你没有得到房子的多边形混合。