请问如何用MATLAB实现radon变换？

theta=0:180theta是radon变换的角度，可以自由设置，这里是从0°到180°一共181个角度，如果只求等分的60个角度的radon变换，可以改成theta=0:3:177；如果只求一个特殊角度的radon变换可以写成theta=x这里x就是希望的角度兆察。

R = radon(I, theta)

R是存储radon变换的值，它是一个矩阵，列数是theta的个数，表示每一个角度生拦猜纳成一列，行数是被处理的矩阵（I）对角线的长度。

R一般用于inverse radon transform，程序如下：

I1=iradon(R,theta)

imshow(I1)

iradon命令有简没很多参数，具体使用方法请参考http://www.mathworks.com/access/helpdesk/help/toolbox/images/index.html?/access/helpdesk/help/toolbox/images/iradon.html&http://www.mathworks.com/cgi-bin/texis/webinator/search?pr=Whole_site&db=MSS&prox=page&rorder=750&rprox=750&rdfreq=500&rwfreq=500&rlead=250&sufs=0&order=r&whole=Whole_site&entire_flag=1&is_summary_on=1&ResultCount=10&query=inverse+radon+transform

：）

拉东变换曾称拉当变换。拉东变换是奥地利数学家拉东（Radon，1917年）在数学研究中首先推导出建立图像的理论。

数学原理

拉东变换(或称经腊猛典Radon变换誉局清)是由奥地利数学家J.Radon于1917年提出来的. 作为积分几何学的基石 ,它为一大类图像重构( 层析成像)问题提供了一个统一的数学基础 ,已被广泛应用于物理、医学庆前、天文、分子生物、 材料科学、核磁共振、无损检测、地球物理等方面.

在地震资料处理中,积分路径为线性的Radon变换又称τ- p变换(τ记作垂直双程旅行时间-截断时间,p为射线参数-斜率),或称为倾斜叠丛纯加.

Ψ(p,τ)=∑Φ(x,τ+px)

τ=t－px(2.73)

p=dt/dx=1/v*=sinθ/v

Fourier投影定理(又称切片投影祥祥定理)证明Radon变换和Fourier变换有明确的对等谨郑搏关系,凡能用Fourier变换解决的问题都能用Radon变换解决

---