我刚刚从Basemap转换到cartopy,用于绘制星空图,我更喜欢它。
(主要原因是在某些电脑上Basemap出现故障,我无法修复。)
我唯一遇到的问题是如何得到一个Tissot圆(用于显示望远镜的视野锥形)。
这是一些绘制随机星星的示例代码(我在真实场景中使用目录):
import matplotlib.pyplot as plt
from cartopy import crs
import numpy as np
# create some random stars:
n_stars = 100
azimuth = np.random.uniform(0, 360, n_stars)
altitude = np.random.uniform(75, 90, n_stars)
brightness = np.random.normal(8, 2, n_stars)
fig = plt.figure()
ax = fig.add_subplot(1,1,1, projection=crs.NorthPolarStereo())
ax.background_patch.set_facecolor('black')
ax.set_extent([-180, 180, 75, 90], crs.PlateCarree())
plot = ax.scatter(
azimuth,
altitude,
c=brightness,
s=0.5*(-brightness + brightness.max())**2,
transform=crs.PlateCarree(),
cmap='gray_r',
)
plt.show()
我该如何在该图像中添加一个特定半径的Tissot圆(以度为单位)? https://en.wikipedia.org/wiki/Tissot%27s_indicatrix
radius
的n_points
个点吗?因此,fwd
会计算给定方向上给定坐标加上半径的端点? - MaxNoe