%matplotlib inline
import requests
import geopandas
import matplotlib.pyplot as plt
source_url = ("https://www.naturalearthdata.com/"\
"http//www.naturalearthdata.com/"\
"download/10m/cultural/ne_10m_airports.zip"
)
source_url
air = geopandas.read_file(source_url)
air.head()
- Project airports to Web Mercator
air = air.to_crs(epsg=3857)
- Extract XY coordinates from a point
GeoDataFrame
air["x"] = air.geometry.x
air["y"] = air.geometry.y
plt.scatter(air["x"], air["y"])
air.to_file('airports_clean.geojson', driver='GeoJSON')