brightness_temperature

astropy.units.equivalencies.brightness_temperature(frequency, beam_area=None)

定义Jy/sr和“亮度温度”之间的转换， \(T_B\) ，在凯尔文。亮度温度是射电天文学中非常常用的单位。例如，见“射电天文学工具”（Wilson 2009）eqn 8.16和eqn 8.19（这些页面可在 google books ）

\(T_B \equiv S_\nu / \left(2 k \nu^2 / c^2 \right)\)

如果输入为Jy/beam或Jy（假设它来自单个光束），则光束面积对于该计算至关重要：亮度温度与光束面积成反比。

参数:

频率 ： Quantity数量

被观察者 spectral 相等的 Unit （例如，频率或波长）。这个变量被命名为“频率”，因为它在射电天文学中更常用。向后兼容性注意：以前版本的亮度-温度等效性使用了关键字 disp ，不再受支持。

beam_area : Quantity [:ref: 'solid angle' ]数量 [:ref: 'solid angle']

以角度单位表示的光束面积，即甾体当量

实例

Arecibo C波段光束：

>>> import numpy as np
>>> from astropy import units as u
>>> beam_sigma = 50*u.arcsec
>>> beam_area = 2*np.pi*(beam_sigma)**2
>>> freq = 5*u.GHz
>>> equiv = u.brightness_temperature(freq)
>>> (1*u.Jy/beam_area).to(u.K, equivalencies=equiv)  
<Quantity 3.526295144567176 K>

VLA合成梁：

>>> bmaj = 15*u.arcsec
>>> bmin = 15*u.arcsec
>>> fwhm_to_sigma = 1./(8*np.log(2))**0.5
>>> beam_area = 2.*np.pi*(bmaj*bmin*fwhm_to_sigma**2)
>>> freq = 5*u.GHz
>>> equiv = u.brightness_temperature(freq)
>>> (u.Jy/beam_area).to(u.K, equivalencies=equiv)  
<Quantity 217.2658703625732 K>

任何通用表面亮度：

>>> surf_brightness = 1e6*u.MJy/u.sr
>>> surf_brightness.to(u.K, equivalencies=u.brightness_temperature(500*u.GHz)) 
<Quantity 130.1931904778803 K>