我正在使用 Cartopy(TransverseMercator
投影)绘制 map 。我在轴坐标 (p_a = (0.1, 0.9)
) 中有一个点,我需要它的纬度/经度坐标(即 PlateCarree
投影中的坐标)。我怎样才能做到这一点?
最佳答案
为此,您需要将点从轴坐标转换为显示坐标,然后转换为数据坐标,最后转换为纬度/经度坐标。因此,您需要从 matplotlib 进行转换和 cartopy .
点 p_a = (0.1, 0.9)
似乎在有效的纬度/经度坐标之外(对于默认的 ccrs.TransverseMercator()
)。因此我使用 p_a = (0.6, 0.6)
import cartopy.crs as ccrs
import matplotlib.pyplot as plt
proj = ccrs.TransverseMercator()
proj_cart = ccrs.PlateCarree()
f, ax = plt.subplots(subplot_kw=dict(projection=proj))
ax.coastlines()
# define point
p_a = (0.6, 0.6)
# plot point in Axes coordinates
ax.plot(*p_a, transform=ax.transAxes, marker='o', ms=10)
# convert from Axes coordinates to display coordinates
p_a_disp = ax.transAxes.transform(p_a)
# convert from display coordinates to data coordinates
p_a_data = ax.transData.inverted().transform(p_a_disp)
# convert from data to cartesian coordinates
p_a_cart = proj_cart.transform_point(*p_a_data, src_crs=proj)
# make sure we are correct
ax.plot(*p_a_cart, transform=proj_cart, marker='x', ms=10)
这会产生下图:
关于python - Cartopy:将点从轴坐标转换为纬度/经度坐标,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/56662941/