Geo Pandas Canada

import pandas as pd
import geopandas
import matplotlib.pyplot as plt
from geopy.geocoders import Nominatim

geolocator = Nominatim(user_agent="rajarcsp")

def get_lat_lang(city):
    location = geolocator.geocode(city)
    return location.latitude, location.longitude

to_lat, to_long = get_lat_lang('Toronto')
mt_lat, mt_long = get_lat_lang('Montreal')
wl_lat, wl_long = get_lat_lang('Waterloo')

df = pd.DataFrame(
    {'City': ['Toronto', 'Montreal', 'Waterloo', ],
     'Country': ['Canada', 'Canada', 'Canada'],
     'Latitude': [to_lat, mt_lat, wl_lat],
     'Longitude': [to_long, mt_long, wl_long]})

gdf = geopandas.GeoDataFrame(
    df, geometry=geopandas.points_from_xy(df.Longitude, df.Latitude))

print(gdf.head())

world = geopandas.read_file(geopandas.datasets.get_path('naturalearth_lowres'))

# We restrict to South America.
ax = world[world.continent == 'North America'].plot(
    color='white', edgecolor='black')

# We can now plot our GeoDataFrame.
gdf.plot(ax=ax, color='red')

plt.show()