import earthpy.plot as ep
country_data_url = "https://www.naturalearthdata.com/http//www.naturalearthdata.com/download/50m/cultural/ne_50m_admin_0_countries.zip"
et.data.get_data(url=country_data_url)

# Create a path to the countries shapefile
countries_path = os.path.join("data", "earthpy-downloads",
                              "ne_10m_admin_0_countries", 
                              "ne_10m_admin_0_countries.shp")


gpd.read_file(countries_path)
# Create a path to the countries shapefile
countries_path = os.path.join("data", "earthpy-downloads",
                              "ne_10m_admin_0_countries", 
                              "ne_10m_admin_0_countries.shp")
f, ax1 = plt.subplots(figsize=(10,6))
im = countries.plot(column='POP_EST',
              ax=ax1)
ep.colorbar(im)
plt.show()