Как нарисовать круг с долготой, широтой и радиусом (км) с помощью Cartopy в Python ⇐ Python

[Anonymous]

Код: Выделить всё

#!/usr/bin/env python

import os, sys
import pandas as pd
import cartopy
import cartopy.crs as ccrs
import matplotlib.pyplot as plt
import shapely.geometry as sgeom
import numpy as np
from cartopy.geodesic import Geodesic

if __name__ == '__main__':

stn = pd.read_csv('obs_station.csv')
gd = Geodesic()

lcc = ccrs.LambertConformal(central_longitude=126., central_latitude=38.)
fig = plt.figure(figsize=(7,7))
ax = fig.add_subplot(111, projection=lcc)
ax.coastlines(resolution='50m')
geoms = []
for lon, lat in zip(stn['longitude'], stn['latitude']):
cp = gd.circle(lon=lon, lat=lat, radius=250000.)
geoms.append(sgeom.Polygon(cp))
ax.add_geometries(geoms, crs=lcc, edgecolor='r')
ax.set_extent([120., 133., 30., 43.])
plt.show()

Файл obs_station.csv содержит несколько координат долготы и широты.
Используя приведенный выше код, я пытаюсь нарисовать круги определенного радиуса (250 км).
Но на карте ничего нет, как показано ниже. На карте отображается только береговая линия.
В чем проблема?
Результат: