Как нарисовать круг с долготой, широтой и радиусом (км) в Cartopy of 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 км).
Но на карте ничего нет, как показано ниже. Покажите карту только с береговыми линиями.
Я не знаю, в чем проблема.
Результат: