跳转到内容

User:NtuEEsaber/沙盒

维基百科,自由的百科全书
雷登轉換將函數 映射到
本圖是將下圖做雷登轉換後得到的影像,越亮的區域代表值越大,黑色的區域為0。
原始函數是白色區域為1,黑色區域為0。

數學上,雷登變換是一種積分變換,這個轉換將函數 轉換成一個定義在二維空間上的函數,而在某線上的值等於對該條線做線積分的值,

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

簡介

[编辑]

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

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

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

定義

[编辑]

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

對於一個弧長 的線,可以把直線 變成一個參數式

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

更進一步,我們可以把推廣到,對一個緊緻台(compact support)的連續函數,轉換後的函數是定義在 超平面上,

對所有

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

也能把仿射子空間,而這種推廣雷登變換的特殊情況被廣泛應用在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]

參見

[编辑]

注釋

[编辑]
  1. ^ https://en.wikipedia.org/wiki/Radon_transform#CITEREFRadon1917
  2. ^ https://en.wikipedia.org/wiki/Radon_transform#CITEREFRoerdink2001
  3. ^ http://statweb.stanford.edu/~candes/math262/Lectures/Lecture09.pdf
  4. ^ http://statweb.stanford.edu/~candes/math262/Lectures/Lecture10.pdf
  5. ^ http://statweb.stanford.edu/~candes/math262/Lectures/Lecture10.pdf

參考

[编辑]