Пример #1
0
def admin_UGANDA_OSM(filename, shapeu, admins):
    logo.starting("Attributes read", parseosm.getNbRelation())
    for relationid in parseosm.getIterRelation():
        logo.progress()
        relation = parseosm.getRelation(relationid)

        admins[relationid] = { "name" : relation['tags']['name'],
                               "level" : int(relation['tags']['admin_level']),
                               "inner" : set(),
                               "outer" : set(),
                             }
        if 'old_name' in relation['tags']:
            admins[relationid]["old_name"] = relation['tags']['old_name']

        wayids = [ data[1] for data in relation['members']
                   if data[0] == 'way' and data[2] in ('outer', 'inner') ]
        lineset = set()
        for wayid in wayids:
            pntinring = []
            linegeometry = parseosm.getGeometryWay(wayid)
            for lon, lat in linegeometry:
                pointid = shapeu.getPoint(lon, lat)
                if pointid is not None:
                    pntinring.append(pointid)

            for pnt in xrange(1, len(pntinring)):
                if pntinring[pnt-1] ==  pntinring[pnt]:
                    # If 2 coordinates after rounding give the same point id
                    continue
                segment = shapeu.getSegment(pntinring[pnt-1], pntinring[pnt])
                lineset.add(shapeu.getLine(segment))

        # Update each administrative level
        admins[relationid]["outer"].update(lineset)
    logo.ending()
Пример #2
0
def read_UGANDA_OSM(filename, shapeu):
    parseosm.parse_xml(open(filename).read())
    shapeutil.precision = 14  # don't do aggressive rounding
    shapeutil.testnearest = []  # nor neighbour hack

    # Read each polygon and build the connection arrays (point, segment, line)
    logo.starting("Geometry read", parseosm.getNbRelation())
    for relationid in parseosm.getIterRelation():
        logo.progress()
        relation = parseosm.getRelation(relationid)
        wayids = [ data[1] for data in relation['members']
                   if data[0] == 'way' and data[2] in ('outer', 'inner') ]

        for wayid in wayids:
            linegeometry = parseosm.getGeometryWay(wayid)
            lon1, lat1 = linegeometry[0]
            for lon2, lat2 in linegeometry[1:]:
                shapeu.makeSegment(lon1, lat1, lon2, lat2)
                lon1, lat1 = lon2, lat2
    logo.ending()