查看“︁投影切片定理”︁的源代码

在数学上，二维情况下的'''投影切片定理'''（{{Lang-en|projection-slice theorem}}）（或称'''中心切片定理'''，{{Lang-en|central slice theorem}}、'''傅里叶切片定理'''，{{Lang-en|Fourier slice theorem}}）表明以下两个运算的结果相等：

* 将二维函数<math>f\left(\mathbf{r}\right)</math>先投影到一维的线上（即进行[[拉东变换]]），并对投影结果进行[[傅里叶变换]]。
* 对相同的函数先进行二维傅里叶变换，然后通过平行于投影线的原点对其进行切片。

以算子形式表示，令：

* <math>F_1</math>和<math>F_2</math>分别为一维和二维傅里叶变换算子
* <math>P_1</math>为投影算子（即将二维函数投影为一维的线）
* <math>S_1</math>为切片算子（从函数中通过原点提取一维切片）

则有：

<math>F_1 P_1 = S_1 F_2</math>

以上结论可以推广到高维情况。

该定理可以应用于医学[[CT扫描]]中，此时，“投影”是对体内器官的X光成像。对成像结果的傅里叶变换可以看作是体内器官三维密度的傅里叶变换的切片，而这些切片通过插值可以构造密度的完整傅里叶变换。对得到的完整结果应用傅里叶逆变换可以得到目标体的密度。这一技术由{{Link-en|罗纳德·布雷斯韦尔|Ronald_N._Bracewell}}于1956年为射电天文学问题开发。<ref>{{cite journal |last=Bracewell |first=Ronald N. |year=1956 |title=Strip integration in radio astronomy |url=https://www.publish.csiro.au/ph/pdf/ph560198 |journal=Australian Journal of Physics |volume=9 |issue=2 |pages=198–217 |bibcode=1956AuJPh...9..198B |doi=10.1071/PH560198 |doi-access=free |access-date=2024-04-22 |archive-date=2023-08-14 |archive-url=https://web.archive.org/web/20230814220211/https://www.publish.csiro.au/ph/pdf/ph560198 |dead-url=no }}</ref>

== ''N''维的投影切片定理 ==
在''N''维情况下，投影切片定理表明''N''维函数<math>f\left( \mathbf{r} \right)</math>投影于''m''维[[欧几里得空间#与流形的关系|线性子流形]]的傅里叶变换等同于该函数''N''维傅里叶变换的''m''维切片，切片由''m''维的线性子流形组成，该流形穿过傅里叶空间中的原点，并平行于投影子流形。以算子形式表示该定理，则有：

: <math>F_mP_m=S_mF_N.\,</math>

== 傅里叶切片定理的推广 ==
除了推广到''N''维空间外，投影切片定理还可以通过改变基函数得到进一步推广。<ref name="NgFourierSlicePhotography">{{cite journal |last=Ng |first=Ren |year=2005 |title=Fourier Slice Photography |url=https://graphics.stanford.edu/papers/fourierphoto/fourierphoto-600dpi.pdf |journal=ACM Transactions on Graphics |volume=24 |issue=3 |pages=735–744 |doi=10.1145/1073204.1073256 |access-date=2024-04-22 |archive-date=2024-05-18 |archive-url=https://web.archive.org/web/20240518035112/http://www.graphics.stanford.edu/papers/fourierphoto/fourierphoto-600dpi.pdf |dead-url=no }}</ref>为了方便表示，我们将基的变化表示为矩阵<math>B</math>，该矩阵为大小为<math>N \times N</math>的可逆矩阵。广义傅里叶切片定理则可以表示为：

<math>F_m P_m B = S_m \frac{B^{-T}}{|B^{-T}|} F_N</math>

其中<math>B^{-T}=(B^{-1})^T</math>是对基变换矩阵<math>B</math>求逆的转置。

== 二维情况下的证明 ==
[[File:ProjectionSlice.png|center|frame|二维投影切片定理的图示。<math>f\left(\mathbf{r}\right) 
</math>和<math>F \left( \mathbf{k} \right) 
</math> 是二维傅里叶变换对。<math>f\left(\mathbf{r}\right) 
</math>在 x 轴上的投影是<math>f\left(\mathbf{r}\right) 
</math>沿平行于 y 轴的方向的线积分，记为<math>p\left(x \right) 
</math>。通过<math>F \left( \mathbf{k} \right) 
</math>的切片位于<math>k_x 
</math> 轴上，该轴平行于<math>x 
</math> 轴并，记为 <math>s\left(k_x\right) 
</math>。投影切片定理指出<math>p\left(x \right) 
</math> 和<math>s\left(k_x\right) 
</math> 是一维傅里叶变换对。]]
二维的情况下的投影切片定理容易证明。此处取投影线作为 x 轴，由于可以通过平移和旋转变换投影线到x轴上，所该选取方式具有普遍性。


设<math>f\left(x,y\right) 
</math>为二维函数，其在x轴方向方向上的投影<math>p\left(x \right) 
</math>为：

: <math>p(x)=\int_{-\infty}^\infty f(x,y)\,dy.</math>

<math>f\left(x,y\right) 
</math>的二维傅里叶变换为：

: <math>
F(k_x,k_y)=\int_{-\infty}^\infty \int_{-\infty}^\infty
f(x,y)\,e^{-2\pi i(xk_x+yk_y)}\,dxdy.
</math>

在傅里叶域的切片<math>s(k_x)</math>为：

: <math>s(k_x)=F(k_x,0)
=\int_{-\infty}^\infty \int_{-\infty}^\infty f(x,y)\,e^{-2\pi ixk_x}\,dxdy
</math>
::: <math>=\int_{-\infty}^\infty
\left[ \int_{-\infty}^\infty f(x,y)\,dy\right]\,e^{-2\pi ixk_x} dx 
</math>
::: <math>=\int_{-\infty}^\infty p(x)\,e^{-2\pi ixk_x} dx
</math>

其正是<math>p\left(x \right) 
</math>的傅里叶变换。更高维的证明可以从以上例子中推广得到。

== FHA（Fourier-Hankel-Abel）循环 ==
如果二维函数 <math>f\left(\mathbf{r}\right) 
</math>是圆对称的，它可以表示为<math>f\left(r\right) 
</math>，其中<math>r=\left| \mathbf{r} \right| 
</math>。此时，投影到任何投影线上都是<math>f\left(r\right) 
</math>的{{Link-en|阿贝尔变换|Abel_transform}}。<math>f\left(\mathbf{r}\right) 
</math> 的二维傅里叶变换等价于<math>f\left(r\right) 
</math>零阶[[漢克爾變換|汉克尔变换]]给出的圆对称函数，因此它也将表示通过原点的任何切片。投影切片定理指出，投影的傅里叶变换等同于傅里叶变换的切片，以算子形式表示为：

: <math>F_1 A_1 = H,</math>

其中<math>A_1 
</math>表示阿贝尔变换算子，将二维圆对称函数投影到一维的线上，<math>F_1 
</math>表示一维傅里叶变换算子，<math>H 
</math>表示零阶汉克尔变换算子。

== 扇形束和锥形束CT的推广 ==
投影切片定理适用于具有平行束投影的CT图像重建，但其并不适用于扇形束或锥形束的CT。该定理于1995年被Shuang-ren Zhao扩展到了扇形束和锥形束CT图像重建。<ref name="ZhaoFSliceThoerem2">{{cite book|author=Zhao S.R. and H.Halling|title=1995 IEEE Nuclear Science Symposium and Medical Imaging Conference Record|chapter=A new Fourier method for fan beam reconstruction|volume=2|year=1995|pages=1287–91|doi=10.1109/NSSMIC.1995.510494|isbn=978-0-7803-3180-8|s2cid=60933220}}</ref>

== 参见 ==

* [[拉東變換#與傅立葉變換的關係|拉东变换#与傅里叶变换的关系]]

== 参考文献 ==
{{Reflist}}
[[Category:图像处理]]
[[Category:积分变换]]