跳轉到

6.2 MRI 與遙望星空:k 空間裡的人體與宇宙#

1973 年,兩篇論文幾乎同時發表,都宣稱能用數學方法從間接測量中「重建」出影像:一篇關於人體(核磁共振成像,MRI),一篇關於宇宙(電波天文學的影像合成)。兩者背後是同一套數學——傅立葉轉換的逆問題。這一節講這兩個故事,以及它們共同的核心概念:k 空間。

從共振到影像:MRI 的誕生#

核磁共振的物理基礎在 1930–40 年代建立:拉比(Isidor Rabi,1944 年諾貝爾獎)測得原子核在磁場中的共振;1946 年布洛赫與柏塞爾分別在史丹佛與哈佛測得凝聚態物質的核磁共振訊號(1952 年諾貝爾物理獎)。核磁共振譜學(NMR)隨即成為化學分析的主力工具——問題是它只能給你「一管樣品的平均資訊」,不是影像。

1971 年,雷蒙德·達馬迪安(Raymond Damadian)發現老鼠的正常組織與腫瘤組織的核磁共振弛豫時間不同,提出用 NMR 偵測疾病的構想。1973 年,紐約州立大學石溪分校的化學家保羅·勞特伯(Paul Lauterbur)在《Nature》發表關鍵一步:在主磁場上疊加線性梯度磁場,使空間中每個位置的共振頻率不同——於是訊號的頻譜直接編碼了空間資訊。他稱之為「以梯度磁場做核磁共振的影像重建」(zeugmatography)。

同年,英國諾丁漢的物理學家彼得·曼斯菲爾德(Peter Mansfield)發展出利用梯度與回波做快速空間選擇的方法,後來演變成至今仍是 MRI 主力的「回波平面成像」(EPI)。兩人共享 2003 年諾貝爾生理學或醫學獎。

k 空間:頻率域的另一個名字#

勞特伯與曼斯菲爾德的方法,在 1970 年代中期被 Richard Ernst(蘇黎世 ETH,1991 年諾貝爾化學獎)與其同事系統化為今天通行的形式:二维傅立葉成像。在 MRI 中,掃描器並不直接測量人體斷層的影像,而是測量影像的二維傅立葉轉換的取樣——這個頻率資料平面被稱為 k 空間。掃描的過程,就是用梯度脈衝一點一點地描過 k 空間;成像就是一次逆傅立葉轉換。

\[\text{影像}(x,y) = \iint S(k_x,k_y)\, e^{i(k_x x + k_y y)}\, dk_x dk_y\]

k 空間的中心是低頻(影像的整體輪廓與對比),邊緣是高頻(細節與銳利度)。這也解釋了 MRI 的實務規則:k 空間取樣不足時,影像會出現環狀偽影(吉布斯現象的醫學版!);運動會在 k 空間留下痕跡,重建成幽靈般的殘影。傅立葉分析的每個性質,都在掃描室裡有了臨床對應。

電波望遠鏡:同一套數學,望向宇宙#

1973 年的另一篇里程碑出自電波天文學。自 1940 年代以來,天文學家想用分布在各地的電波望遠鏡做「干涉」:兩台相距甚遠的望遠鏡同時觀測同一目標,訊號相關後得到可見度函數——在數學上,這正是天空亮度分布的傅立葉轉換的取樣(范西特–澤尼克定理,Van Cittert–Zernike,1930 年代)。基線越長,取樣到的空間頻率越高,解析度越好。

英國劍橋的萊爾(Martin Ryle,1974 年諾貝爾物理獎)發展了地球自轉合成(earth-rotation aperture synthesis):地球轉動時,基線相對天空的方向緩慢掃過 k 空間,一整夜得到一整片取樣。1974 年的論文(如著名的 Very Long Baseline Interferometry 與合成成像方法的成熟)與後續幾十年,把這條路推到極致。今日事件視界望遠鏡(EHT)用橫跨地球的望遠鏡陣列拍出 2019 年 M87 黑洞影像——解析度所需基線相當於地球直徑——正是這條系譜的最新成員:資料相關得到 k 空間取樣,再用(含正則化的)傅立葉反演重構影像。

人體與黑洞,共享同一個數學座標系:k 空間。差別只在距離的尺度,不在原理。

醫療、科研與產業的共生#

MRI 的商業化是另一段制度史:1970 年代達馬迪安創辦 Fonar 公司爭取專利與市場,奇異(GE)、西門子等醫療器材巨頭投入研發;1980 年代初 FDA 批准,MRI 迅速成為醫院標配。全球每年數以億計的掃描,每張背後都是即時的 FFT。而天文干涉成像則由國家級科研機構(美國國家電波天文台、歐洲南天天文台)支撐,經費來自公部門對基礎科學的長期投資——同一套數學,兩種截然不同的資助與治理模式,都成功了。

傅立葉的方程,就這樣同時寫進了醫院的收費系統與諾貝爾獎的名單。但在讚嘆之餘,理論家們早已注意到傅立葉的盲區:它擅長回答「含哪些頻率」,卻拙於回答「何時含這些頻率」。下一步的突破——小波——正由此而來。