如何使用python创建矩1图（也称为谱图最大值坐标）？

Question

如何使用python（astropy）创建矩1图（也称为谱图最大值坐标）？

astropy

可以用来绘制第9时刻吗？我使用

.moment(order=9)

获得的数据看起来有点奇怪，我在

astropy

文档中只看到了三个阶矩图。

astropy

只能绘制最多三个订单吗？

Answer 1

moment

s 是沿着给定强度和光谱坐标的光谱轴计算的。光谱维度是表征时空几何和拓扑的量，如速度、频率、波长或能量。

因此，对于第

0

，谱线上的积分强度是（立方）单位乘以谱轴单位（例如 K k/ms）：moment

import astropy.units as u
from astropy.utils import data
from spectral_cube import SpectralCube
fn = data.get_pkg_data_filename('tests/data/example_cube.fits', 'spectral_cube')
cube = SpectralCube.read(fn)
print(cube)
moment_0 = cube.moment(order=0)
print(moment_0)

st

1

是谱线的强度加权速度，因此单位与谱线单位相同（例如km/s）。

moment

nd

2

是速度色散或宽度或谱线，因此单位是谱轴平方（例如 km^2 / s^2）。根据定义，不存在第

moment

第

9

：

https://spectral-cube.readthedocs.io/en/latest/api/spectral_cube.SpectralCube.html#spectral_cube.SpectralCube.moment

，https： //spectral-cube.readthedocs.io/en/latest/moments.html，https://spectral-cube.readthedocs.io/en/latest/ 并不是说你不能有第

moment

th

9

，但它会像 km^9 / s^9。