我需要使用Cartopy绘制地图。任务要求如下:
使用正射投影,以汉堡为中心,绘制汉堡到东京的大圆路径。
import cartopy
import cartopy.crs as ccrs
ax = plt.axes(projection=ccrs.Orthographic(central_longitude = long_ham, central_latitude = lat_ham))
ax.coastlines()
long_ham, lat_ham = 9.993682, 53.551086
long_tok, lat_tok = 139.839478, 35.652832
plt.plot([long_ham, long_tok], [lat_ham, lat_tok],
color='blue', linewidth=2, marker='o',
transform=ccrs.Geodetic(),
)
plt.text(long_ham , lat_ham , 'Hamburg',
horizontalalignment='right',
transform=ccrs.Geodetic())
plt.text(long_tok , lat_tok , 'Tokyo',
horizontalalignment='left',
transform=ccrs.Geodetic())
# Save the plot by calling plt.savefig() BEFORE plt.show()
plt.savefig('great circle distance.pdf')
plt.savefig('great circle distance.png')
plt.show()
有没有办法显示整个地球而不是当前获取的局部地图?