我有一张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)?