如何使用 python 创建力矩 1 贴图（也称为速度场或速度贴图）？

Question

我需要使用 Python 创建时刻 1 地图，但我无法使用 SpectralCube 包来生成它（.moment(order=1) 方法）。问题是我卡住了，基本上我写的所有东西都是垃圾。有没有人有矩图的经验并且可以提供帮助？ 我正在使用立方体并使用 astropy 处理拟合。 为了让您确定我的意思，我将从 CASA 文档中进行解释： Moments explained

Answer 1

Whosebug 不是问这个问题的最佳场所；您可能会从天文学专用论坛上提问。

也就是说，您的 moment-1 地图看起来像垃圾的最可能原因是您包含大量噪音。由于第一时刻是强度加权坐标，如果你有一个包含大量噪声像素的大立方体，你可能只会得到噪声。一个好的方法是mask the data before computing the moments。如果您仅使用在噪声之上具有显着检测的像素，则矩图通常会干净。