ualsg / hdb3d-code Goto Github PK

query = """
[out:json][timeout:30];
// area-of-interest here
area[name='Cape Peninsula University of Technology (Bellville Campus)']->.a;

// --when areas have duplicate names given the world has a limited amount of uniquely named places
// --main area
//area[name="Western Cape"]->.b;
// -- target area ~ can be way or relation
//relation(area.b)[name="Woodstock"];
//map_to_area -> .a;

(
  (
    // I want all buildings
    way[building](area.a);

    // plus every building:part
    way["building:part"](area.a);
    // and relation type=multipolygon ~ to removed courtyards from buildings
    relation["building"]["type"="multipolygon"](area.a);
  );
-
  // excluding buildings with relation type=building role=outline
  // to remove outlines that surround building:part
  (
    // for every way in the input set select the relations of which it is an "outline" member
    rel(bw:"outline")["type"="building"];
    // back to the ways with role "outline"
    way(r:"outline");
  );
);
out body;
>;
out skel qt;
"""

url = "http://overpass-api.de/api/interpreter"
r = requests.get(url, params={'data': query})
#rr = r.read()
gj = osm2geojson.json2geojson(r.json())