def find(self, lat, lon):
    'Find every NUTS region for this lat, lon'
    p = geometry.Point(lon, lat)
    nuts = []
    for region in self.shapes['features']:
        s = geometry.shape(region['geometry'])
        if s.contains(p):
            nuts.append(region['properties'])
    return sorted(nuts, key=(lambda row: row['LEVL_CODE']))
Find every NUTS region for this lat, lon
def find(self, lat, lon):
    
    p = geometry.Point(lon, lat)
    nuts = []
    for region in self.shapes['features']:
        s = geometry.shape(region['geometry'])
        if s.contains(p):
            nuts.append(region['properties'])
    return sorted(nuts, key=(lambda row: row['LEVL_CODE']))