Box.shape is not implemented

[Kylotan]

I am trying to perform a query and limit it to an area around a coordinate. Since I could not find a way to do it with a radius (around doesn't do intersections with Ways, apparently), I've tried to create a bounding box and pass that as a shape into an intersecting call.

Unfortunately, the docs here (https://docs.geodesk.com/python/Box#properties) seem to be wrong and Box.shape doesn't seem to exist.

import geodesk
my_box = geodesk.Box(1, 2, 3, 4)
its_shape = my_box.shape

This gives AttributeError: Attribute not found for shape.

I don't see any mention of 'shape' in PyBox.h so I'm unsure why it's in the docs. (Same goes for buffered - that doesn't appear anywhere either.)

---

(I would appreciate any sort of workaround that can find ways within N metres of a given longitude and latitude, e.g. by generating some other shape I can pass to intersecting)

[clarisma]

Thanks for reporting this!

The missing shape attribute is an old issue (#6), high time for me to implement it. Will schedule it this week.

around doesn't do intersections with Ways, apparently

An around query should do exactly what you have in mind. It should pick up features of any type. Do you have a specific case where it is failing?

[Kylotan]

Thanks for the rapid response! Regarding the around/ways question, here's my test case:

import geodesk
features = geodesk.Features(geodesk_filename) # replace with path to europe.gol or united_kingdom.gol
query = """nwr[amenity=pub,bar,cafe,nightclub,"night club","Arts Centre","Beer garden"]"""
# https://www.openstreetmap.org/#map=19/51.551002/-0.140316 - "The Assembly House"
lat, lng = 51.551002, -0.140316  
all_results = features(query).around(meters=250, lat=lat, lon=lng)

# The query finds things
assert all_results.count > 0

# It does find 'Lilac' next door, so this assert passes
# https://www.openstreetmap.org/node/12251483338
assert any("Lilac" in feature.str("name") for feature in all_results)

# It doesn't find 'The Assembly House', so this assert fails
# https://www.openstreetmap.org/way/973132221
assert any("The Assembly House" in feature.str("name") for feature in all_results)

I've tried this with several venues and it always returns Nodes but never Ways.

[clarisma]

You will need to add an a type to the query, otherwise areas won't be returned:

na[amenity=pub,bar,cafe]

• w only returns lineal ways.
• a returns ways and relations that represent areas
• r returns relations that aren't areas

That's a common confusion for users who are accustomed to the Overpass query syntax. See GOQL Reference.

GOQL interprets n, w, a, r as geometric types, rather than primitive types. To find out which to use, you can view the wiki page for a given OSM tag, such as amenity=pub, and look for these symbols on the right:

n / w / a / r

So, for pubs you would use na, since they can be nodes or areas.

[Kylotan]

Thank you for the explanation. With the changes you suggested, the original query works fine.

I think that OSM Wiki page is more confusing than useful however! They need to fix their terminology around Ways and Areas because it is completely contradictory at the moment.