How to get the lat long polygon for a given airport

What I need to do is to get the bounding box polygon for all airports worldwide. The only way I'm aware of do this, is by using esy.osm.pbf. I've written the following code:

osm = esy.osm.pbf.File(file_path)
lat_lons = OrderedDict()
index = 0
print("There are " + str(len(airport.refs)) + " nodes to go through.")
for ref in airport.refs:
    print(str(index))
    index = index + 1
    nodes = [entry for entry in osm if entry.id == ref]
    lonlat = nodes[0].lonlat
    print(lonlat)
    print(type(lonlat[1]))
    print(type(lonlat[0]))
    lat_lons[ref] = ([lonlat[1], lonlat[0]])
return lat_lons

But it takes much too long to do this, that's just for a single airport, and I need to do it for all airports worldwide. Is there a better way of fetching the latitudes and longitudes for the polygons for all airports worldwide?