[GIS] Python: Geopy’s Nominatim reverse geocoding: Access Type (e.g. amenity)

nominatimopenstreetmappython

I have some corrdinates (lat/long) and I want to get further information on them. When using the nominatim homepage, I get quite a good overview of the available information, which exist. But when using my python code, I only get the basic address information. In the end, I want to add especially the amenity information, which I will later on use to filter my data.

The picture shows, what I want (in yellow): https://i.stack.imgur.com/Gn2C6.png

And his is what I do in python and what I get:

p = geopy.point.Point(48.705293, 8.984022)

r = Nominatim(user_agent="my-application").reverse(p)

print(r)

>>Logistikfläche, Calwer Straße, Sindelfingen (Stadt), Sindelfingen, Landkreis Böblingen, Baden-Württemberg, 71063, Deutschland

Best Answer

As far as I know, you can show bunch of data from Nominatim reverse geocode. This is an example:

from geopy.geocoders import Nominatim

locator = Nominatim(user_agent=”myGeocoder”)
coordinates = “53.480837, -2.244914”
location = locator.reverse(coordinates)
location.raw

If you just wanna look for the address, simply type location.address. Whereas, raw element of location is used to show a lot of data like this:

{‘address’: {‘building’: ‘Eagle Insurance Buildings’, 
‘city’: ‘Manchester’, 
‘country’: ‘United Kingdom’, 
‘country_code’: ‘gb’, 
‘county’: ‘Greater Manchester’, 
‘house_number’: ‘68’, 
‘postcode’: ‘M2 4JG’, 
‘road’: ‘Cross Street’, 
‘state’: ‘England’, 
‘state_district’: ‘North West England’, 
‘suburb’: ‘City Centre’}, 
‘boundingbox’: [‘53.480856’, ‘53.4810634’, ‘-2.2451761’, ‘-2.2449576’], 
‘display_name’: ‘Eagle Insurance Buildings, 68, Cross Street, City Centre, Manchester, Greater Manchester, North West England, England, M2 4JG, United Kingdom’,
‘lat’: ‘53.4809597’, 
‘licence’: ‘Data © OpenStreetMap contributors, ODbL 1.0. https://osm.org/copyright', 
‘lon’: ‘-2.24506682746292’, 
‘osm_id’: 37139875, 
‘osm_type’: ‘way’, 
‘place_id’: 86197794}

Hope it helps, I refer to this.

Related Solutions

OpenStreetMap – How to Perform Reverse Geocoding with Nominatim

Up to now Nominatim has no support for finding street intersections. However, the developers will welcome your help implementing it.

Perhaps there are other solutions for you, as explained here:

If you know the two streets for which you want to find the intersection, then you can retrieve them via two geocoding-queries and compare their nodes for equality (a road intersection in OSM is defined by a node being part of two or more ways).
If you just know the rough location, then you can grab all surrounding roads (e.g. using the Overpass API) and again start comparing each road with each other one for equal nodes.

Python Pandas – Reverse Geocoding of Pandas DataFrame with Lat/Long Columns

Here the Nominatim geocoder (free to choose) from the GeoPy geocoding Python library was used, for more details, please check the documentation.

Coordinates of point features should be passed as a pair (see this thread for more details), otherwise, you may get this error:

ValueError: Must be a coordinate pair or Point.

Therefore from geopy.point import Point was additionally imported.

The code below also uses the vectorize of NumPy.

import numpy as np
import pandas as pd
from geopy.geocoders import Nominatim
from geopy.point import Point

geolocator = Nominatim(user_agent="test")

def reverse_geocoding(lat, lon):
    try:
        location = geolocator.reverse(Point(lat, lon))
        return location.raw['display_name']
    except:
        return None

df = pd.DataFrame({'id': [1,2,3,4,5],
                   'Latitude': [30.197535, 34.895699, 33.636700, 33.636700, 32.733601],
                   'Longitude': [-97.662015, -82.218903, -84.428101, -84.428101, -117.190000]
                    })

df['address'] = np.vectorize(reverse_geocoding)(df['Latitude'], df['Longitude'])

print(df)

will result in:

   id   Latitude   Longitude                                            address
0   1  30.197535  -97.662015  Austin-Bergstrom International Airport, 3600, ...
1   2  34.895699  -82.218903  Greenville-Spartanburg International Airport, ...
2   3  33.636700  -84.428101  Hartsfield–Jackson Atlanta International Airpo...
3   4  33.636700  -84.428101  Hartsfield–Jackson Atlanta International Airpo...
4   5  32.733601 -117.190000  San Diego International Airport, North Harbor ...