平面反向诺莫投影

Question

我有一个关于旋转的二级场，主要是研究星系的形状扭曲，公式如下：

γ=γ1+iγ2=|γ|exp(−2iϕ)

这里的 ϕ 是方向角。

如果这个量是在天空坐标上测量的，天空坐标是球面坐标，假设测量的中心点是 (α0 ,δ0)，而我们要测量的另一个点是相对于这个中心点的 (α ,δ) 天空坐标。那么，我该如何使用正射投影将球面坐标转换为笛卡尔坐标，来测量 tangential（切向）和 cross（交叉）分量的 γ，并用 healpy 或 matplotlib 来展示这个过程呢？同时，mathematica 中也有类似的功能？

Answer 1

看看这个叫做 basemap 的 matplotlib 工具包。

还有它在 gnomonic 投影 和 坐标转换 上的例子。

一旦你创建了一个叫 m 的 Basemap 对象，你就可以把 纬度, 经度 的坐标转换成地图上的坐标，方法是：

x, y = m(latitude, longitude)

而你也可以进行反向转换，方法是：

latitude, longitude = m(x, y, inverse=True)