brightness_temperature

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

Jy/Srと“明るい温度”の換算を定義しました \(T_B\) ケルビンにいます。亮温は電波天文学で非常によく使われる単位である.例えば、“Tools of Radio Astronomy”(Wilson 2009)EQN 8.16およびEQN 8.19(これらのページは google books ）。

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

入力がJy/ビームまたはJy(単一ビームからのものと仮定する)である場合、ビーム面積は、輝度温度がビーム面積に反比例するという計算に重要である。

パラメータ

周波数、周波数 ： Quantity量

観察されたのは spectral 等価物. Unit (例えば、周波数または波長)。この変数は電波天文学でよりよく使われるため“周波数”と命名された。後方互換性説明：従来バージョンの輝度温度当量使用キーワード disp これ以上支持されていません

beam_area量 [“立体角度”]

角度単位で表されるビーム面積、すなわち立体当量

実例.

アレシボ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>