E
Learn EoR

電波干渉計 Visibility とフーリエ変換の直感的理解

MWA(基線長 55m)のデータを例に、visibilityの物理的意味とCNNの周波数方向への畳み込みがなぜ不適切かを視覚的に解説します。

1. フーリエ変換とは何か

フーリエ変換の本質は「ある信号を、さまざまな周期(スケール)の波の重ね合わせとして分解する」操作です。下のスライダーで各波の振幅を変えてみてください。

元の信号(3つの波の重ね合わせ)
低周波 (f=2)0.6
中周波 (f=5)0.4
高周波 (f=11)0.3
フーリエ変換の結果

棒の高さ=その周波数成分の振幅。スライダーと連動します。

2. 1次元 → 2次元フーリエ変換

1次元では「時間 → 周波数」でしたが、2次元では「空間の画像 → 空間周波数」になります。2次元の「波」とは空間的な縞模様であり、任意の画像はさまざまな方向・細かさの縞模様の重ね合わせとして分解できます。

2次元の「波」=空間の縞模様

低周波・横縞
(大スケール)

+

高周波・横縞
(小スケール)

+

中周波・斜め縞
(中スケール)

=

重ね合わせ
(これが「画像」)

フーリエ空間 (u-v平面)

各点が一つの縞模様に対応。中心からの距離 = 空間周波数(縞の細かさ)

u-v空間の点をドラッグして縞模様を確認

点を動かすと右の縞模様がリアルタイムに変わります。赤線=縞の方向、緑矢印=k-ベクトル(波の伝播方向)

u-v 空間(点をドラッグ)
対応する縞模様
kx(横方向)
3.0
ky(縦方向)
2.0
空間周波数 |k|
3.6
縞の傾き(水平から)
56°

中心に近いほど太い縞(大スケール)・遠いほど細かい縞(小スケール)。点の位置の角度が縞の法線方向を決める。kx=0なら横縞、ky=0なら縦縞。

3. 干渉計と Visibility の関係

干渉計の2本のアンテナの相関出力(visibility)は、天球上の輝度分布の2次元フーリエ成分をサンプリングしたものです(van Cittert-Zernike定理)。

基線ベクトル b を持つアンテナ対が測定するu-v座標は:

u = b / λ  (基線長 ÷ 観測波長)
実空間(アンテナ配置)基線 b = 55m入射電波(平面波)フーリエ空間 (u-v平面)uvb/λu = b/λ同じ基線(55m)でも観測周波数が変わると λ が変わり、u-v上の位置が変わる

4. 周波数チャンネルとスケールの対応

「スケール」とは天球上の角度的な大きさのことです。角度スケール θ ≈ λ/b であり、u-v平面の中心から遠いほど小さな(細かい)角度スケールに対応します。

観測周波数150 MHz
波長 λ
2.00 m
u = b/λ
27.5 λ
角度スケール θ ≈ λ/b
2.1°

スライダーを動かすと、同じ55m基線でも周波数によって感度を持つスケールが変わることがわかります

5. なぜCNNの畳み込みは周波数方向に不適切か

CNNの畳み込み層は並進不変性を仮定しています:同じカーネルをデータ上でスライドさせても同じパターンが同じ意味を持つ、という前提です。

問題:Visibilityデータの周波数軸上では、チャンネルごとに対応するu値(=角度スケール)が異なります。80MHzでは u≈14.7λ(θ≈3.9°)、300MHzでは u≈55.0λ(θ≈1.0°)。周波数方向にカーネルをスライドさせると、カーネルが「見ている」物理スケールが連続的に変わり、並進不変性が成り立ちません。

赤い枠がCNNカーネル。周波数方向にスライドすると対応する角度スケールが変わってしまう

解決の方向性:各チャンネルのu値に応じてカーネルを適応的に変化させるか、あらかじめu-v空間(物理的に等しいスケール空間)にグリッディングしてから処理する、などのアプローチが考えられます。

まとめ

フーリエ変換=信号を波の重ね合わせに分解する操作

Visibility=天球輝度分布の2次元フーリエ成分のサンプル

u-v座標 = b/λ → 同じ基線でも周波数が変わればu-v上の位置が変わる

スケール(θ ≈ λ/b)=その基線・波長で感度を持つ天球上の角度サイズ

CNNが不適切な理由=周波数方向の並進不変性が物理的に成立しない