雷登變換

數學上,雷登變換是一種積分變換,這個變換將二維平面函數變換成一個定義在二維空間上的一個線性函數(的意思是對做雷登變換),而的值為函數對該條線做積分的值。以右圖為例,黃色區域即是線則是代表

雷登變換將函數 映射到
本圖是將下圖做雷登變換後得到的影像,越亮的區域代表值越大,黑色的區域為0。
原始函數是白色區域為1,黑色區域為0。

雷登變換是Johann Radon在西元1917年提出[1],他也同時提出雷登變換的反變換公式,以及三次空間的雷登變換公式。 三次空間雷登變換,是對一個平面積分(對線積分則是 X-ray transform)。而在不久之後,更高維度的歐幾里得空間的雷登變換被提出,更詳盡的廣義雷登變換要查Integral geometry。 在複數上有和雷登變換相似的Penrose transform,雷登變換被廣泛的應用在斷層掃描,雷登反變換可以從斷層掃描的剖面圖重建出投影前的函數。

簡介

若函數表示一個未知的密度,對做雷登變換,相當於得到投影後的訊號,舉例來說:相當於人體組織,斷層掃描的輸出訊號相當於經過雷登變換的。 因此,可以用雷登反變換從投影後的密度函數,重建原始的密度函數,它也是重建斷層掃描的數學理論基礎,另一個被廣為人知名詞的是三維重建

雷登變換後的訊號稱作"""正弦圖"""(sinogram),因為一個偏離中心的點的雷登變換是一個正弦曲線(sinusoid)。所以對一些小點的雷登變換,會看起來像很多不同振福、相位的正弦函數重疊在一起。

雷登變換可以應用在:X射線電腦斷層掃描條碼掃描器、macromolecular assemblies電子顯微鏡例如:病毒Reflection seismology蛋白質複合體,而且也是雙曲線 偏微分方程(hyperbolic partial differential equations)的解。

定義

令密度函數是一個的定義域為 的緊緻台(compact support)。令為雷登變換的運算子(operator),則是一個定義在 空間中的直線,它的定義如下

可以把直線 改寫成一個弧長的參數式

是直線和原點的距離,而是垂直於的法線和軸的夾角, 接下來,我們可以令當作平面上的新座標系統,把這個座標變換帶入到雷登變換得到

更進一步,我們可以把推廣到歐幾里得空間,對一個緊緻台(compact support)的連續函數做雷登變換後的函數是定義在 超平面上,

積分的對象是自然超平面測度(natural hypersurface measure),而是原本的的高維推廣。可以觀察到對裡的任意元素, 都是某個軌跡方程式的解

是一個單位向量且屬於,n維的雷登變換可以改寫成定義在 上的函數

也可以藉由其他方式將雷登變換推廣,也就是對的k維仿射子空間作(k-dimensional affine subspaces)積分。 而這種推廣雷登變換的特殊情況被廣泛應用在X射線電腦斷層掃描,他的做法是對一條直線積分。

與傅立葉變換的關係

主條目:Projection-slice theorem

雷登變換和傅立葉變換之間有很強的關聯性。單變數的傅立葉變換的定義是

而雙變數的傅立葉變換是

把雷登變換的運算子的表記從 改成 。根據Projection-slice theorem學說,

因此一個初始函數沿著一條線傾角的二維的傅立葉變換,相當於對雷登變換做一維的傅立葉變換。這個結果可以推廣到n維

對偶變換

對偶雷登變換是雷登變換的埃爾米特伴隨。令在空間上的函數,而對偶雷登變換的運算子定義為。作用在

積分的範圍是所有和相交的超平面集合,而測度(measure)是集合特殊的機率測度(Probability measure), 當對著旋轉時,的值不會改變

對於一個二維的雷登變換,其對偶變換是

在影像處理的文章中,對偶變換經常被稱作反向傳播算法(back propagation) [2], 因為

交結性質

根據拉普拉斯算子的定義是

這是一個旋轉不變性的二階微分算子,在空間,半徑的二階倒數(second derivative)

也是旋轉不變性。 而雷登變換與其對偶變換屬於交結運算子(intertwining operator),是因為

重建方法

重建處理是指從投影影像重建一個影像,或是一個函數。重建處理是一種逆問題(inverse problem)。

雷登反變換公式

對於二維雷登變換,最常被使用的解析公式(analytical formula),是Filtered Backprojection Formula或雷登反變換公式,反變換公式為

[3]

函數滿足[4],卷積核 (convolution kernel) 在一些文章中稱作Ramp filter。

不適定問題 (ill-posedness)

直覺上,反變換公式應該和微分類似,。我們可以看的出來反變換公式 的行為類似微分。大致上來說,這個反變換公式把目標奇異化(singular);要如何量化雷登反轉化的不適定問題 (ill-posedness)呢?首先可以寫出

即是前面定義的反變換運算子,且伴隨著(adjoint to)雷登變換,因此,上式變成

複數指數函數,是固有函數 (eigenfunction) , 而特徵值 (eigenvalue)為的奇異值 (singular values) 是, 因為這些奇異值 (singular values)會趨近於0,所以是無界的(unbounded) [5]

反變換公式

外顯(explicit)且計算效率好的雷登反變換公式,以及他的對偶是存在的。n維的反雷登變換可以由[6]

其中

拉普拉斯算子(Laplacian),偽微分算子(pseudodifferential operator)

傅立葉變換的運算子(operator)。

參見

注釋

  1. . [2017-06-29]. (原始内容存档于2017-07-19).
  2. . [2017-06-29]. (原始内容存档于2017-07-19).
  3. (PDF). [2017-06-29]. (原始内容 (PDF)存档于2018-11-25).
  4. (PDF). [2017-06-29]. (原始内容 (PDF)存档于2018-11-25).
  5. (PDF). [2017-06-29]. (原始内容 (PDF)存档于2018-11-25).
  6. Helgason 1984,Theorem I.2.13

參考

  • Deans, Stanley R., , New York: John Wiley & Sons, 1983.
  • Helgason, Sigurdur, , Mathematical Surveys and Monographs 39 2nd, Providence, R.I.: American Mathematical Society, 2008, ISBN 978-0-8218-4530-1, MR 2463854.
  • Helgason, Sigurdur, , Academic Press, 1984, ISBN 0-12-338301-3.
  • Herman, Gabor T., 2nd, Springer, 2009, ISBN 978-1-85233-617-2.
  • Minlos, R.A., , Hazewinkel, Michiel (编), , Springer, 2001, ISBN 978-1-55608-010-4.
  • Natterer, Frank, , Classics in Applied Mathematics 32, Society for Industrial and Applied Mathematics, ISBN 0-89871-493-1
  • Natterer, Frank; Wübbeling, Frank, , Society for Industrial and Applied Mathematics, ISBN 0-89871-472-9.
  • Radon, Johann, , Berichte über die Verhandlungen der Königlich-Sächsischen Akademie der Wissenschaften zu Leipzig, Mathematisch-Physische Klasse [Reports on the proceedings of the Royal Saxonian Academy of Sciences at Leipzig, mathematical and physical section] (Leipzig: Teubner), 1917, (69): 262–277; Translation: Radon, J.; Parks, P.C. (translator), , IEEE Transactions on Medical Imaging, 1986, 5 (4): 170–176, PMID 18244009, doi:10.1109/TMI.1986.4307775.
  • Roerdink, J.B.T.M., , Hazewinkel, Michiel (编), , Springer, 2001, ISBN 978-1-55608-010-4.
  • 埃里克·韦斯坦因. . MathWorld..


This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.