Cartopy投影变换 - 从整数值转换为地图坐标

时间:2017-10-03 11:51:52

标签: python cartopy

我有一张2D图像及其投影信息。

我无法弄清楚如何将像素数转换为正确的投影。

例如

In [23]: image.shape
Out[23]: (1351, 1143)

In [24]: projection
Out[24]: _EPSGProjection(32611)

我现在想要在半径10

的指数(500,700)处绘制一个圆圈

我想我可以做这样的事情来得到正确的预测:

projection.transform_point(500, 700, ???) 

但我不知道基于像素的坐标参考系统应该是什么。

我是否完全错过了有关投影的明显内容?

这是一个基本的例子:

image=np.random.rand(500, 500)

crs_latlon = ccrs.PlateCarree()

subplot_kw = dict(projection=crs_latlon)

fig, ax= plt.subplots(1, 1, figsize=(10, 10), subplot_kw=subplot_kw)

ax.imshow(image, extent=(-20, 20, -15, 15))
ax.set_extent((-30, 30, -25, 25))

ax.coastlines()
plt.show()

这绘制了非洲海岸随机数据的图像。

我现在如何在原始图像的索引处引用某些内容?

e.g。画一个以图像为中心的圆圈(200,400)?

0 个答案:

没有答案