微分几何,大哥们帮帮忙

\u5fae\u5206\u51e0\u4f55,\u5927\u54e5\u4eec\u5e2e\u5e2e\u5fd9

\u53bb\u90ae\u7bb1\u627e\u627e\uff0c\u6211\u4e0d\u8bb0\u5f97\u505a\u6ca1\u505a

\u4e3b\u66f2\u7387 \u4e3b\u65b9\u5411\u80af\u5b9a\u5bf9,\u53ef\u4ee5\u5199\u5171\u8f6d

(一)傅里叶反变换重建

傅里叶反变换重建法是一种重要的变换重建方法。它主要包括以下三个步骤:①建立数学模型,其中已和未知量都是连续实数的函数;②利用反变换公式解未知量;③调节反变换公式以适应离散、有噪声应用的需求。注意上面第二步中在解未知量时理论上可以有多个等价的公式来解。而第三步中,由于离散化时可采用不同的近似,所以理论上等价的方法对实际数据应用的结果会不同。在具体应用中,所测量到的数据对应于许多个离散点(s,θ)上的估计值g(s,θ),而所重建的图像也是一个离散的数组。

设在s和θ上都均匀采样的情况。考虑在N个相差△θ的角度上的投影,在每个角上用M个间距为△s的射线测量,定义整数M+和M-为:�

(6.181)

为了保证一系列射线{(m△s,n△θ):M�-≤m≤M+,1≤n≤N}覆盖单位圆,需要选△θ=π/N和△s=1/M+。此时g(m△s,n△θ)为平行投影的射线数据。设图像区被一个直角网格覆盖,其中K+和K-用类似于式(6.181)的方法定义(K为X方向上的点数),L+和L-也用类似于式(6.181)的方法定义(L为Y方向上的点数)。根据这些定义,一个重建算法就是要通过M×N个测量值g(m△s,n△θ)估计出在K×L个采样点的f(K△x,l△y)。

1.傅里叶变换投影定理变换方法的基础是傅里叶变换投影定理。设G(R,θ)是g(s,θ)对应第一个变量s的一维傅里叶变换,即

(6.182)

F(X,Y)是f(x,y)的二维傅里叶变换:

(6.183)

那么可以证明如下投影定理

G(R,θ)=F(Rcosθ,Rsinθ) (6.184)

即f(x,y)以θ角进行投影的傅里叶变换等于f(x,y)的傅里叶变换在傅里叶空间(R,θ)处的值。换句话说,f(x,y)在与X轴成θ角的直线上投影的傅里叶变换是f(x,y)的傅里叶变换在朝向角θ上的一个截面。

2.傅里叶反变换重建公式根据傅里叶变换投影定理很容易就可以得到傅里叶反变换重建的公式。对式(6.184)两边在直角坐标系中取傅里叶反变换

(6.185)

注意到G(·)是g(·)的傅里叶变换,所以上式是给出g(s,θ),计算f(x,y)的一个重建公式。这也是第一个在投影重建中得到应用的技术。实用中需要加一个窗,以把积分区限制在傅里叶XY平面上的一个有限区域W,这样得到f(x,y)的一个带限逼近fw(x,y)。

(6.186)

现在考虑计算G(·)。实际中G(·)只在一系列θ=θn角取值(θn代表nΔθ)。G(R,θn)可用在一系列采样点(mΔs,θn)对g(·)的求和得到。

如果令R=kΔR(k为整数,ΔR为采样),取ΔR=1/(MΔs),则

G∑(ΔR,θ△n)=Δs(mΔs,θn)exp〔-j2π(mΔs)〕(6.187)

根据G∑(kΔR,θn)可对任意的(X,Y)插值出G(X2+Y2)1/2,arctgY\X。然后由式(6.186)可知

(6.188)

于是fw(x,y)可由下式确定

(6.189)

进一步,如果令△x=1/(U△X),△y=1/(V△Y),则

(6.190)

图6-134右半部给出在(傅里叶)XY平面上的(k△R,θn)位置,用“·”表示。在这些采样位置上可用式(6.187)计算G。这些点在平面上形成一个极坐标模式,各点在极坐标系中是等距分布的。左半部中“+”点表示位置(u△X,υ△Y),需要知道这些点的傅里叶变换以根据式(6.190)计算fw,(u△X,υ△Y)。这些点在平面上形成一个直角坐标网格。可根据在极坐标系中的已知值进行插值以得到这些网格点的估计值。上式只给出有限个(U个k,V个l)估计值fw(k△x,l△y)。对这些参数的双重求和可用FFT快速计算。图6-134(a)右半部的极坐标模式对傅里叶平面的采样不是很有效,一种可能的改进办法是采用图6-134(b)的模式。�

图6-134傅里叶空间的直角和极坐标网络

基于傅里叶反变换的重建技术主要有三个步骤:①对以角θn(n=1,2,…,N)的投影进行一维傅里叶变换;②在傅里叶空间从极坐标向直角坐标插值;③进行二维傅里叶反变换以得到重建图像,由于第三步需要用到二维变换,所以不能根据所获得的部分投影数据重建图像,而必须在获得全部投影数据后才能再重建图像。

题目做好了,在你邮箱里

  • 璋佽兘缁欐垜涓浜涚瀛﹀鐨勪簨杩?
    绛旓細1953骞达綖1954骞村湪鍖椾含鍥涗腑浠绘暀锛屽洜鍙i娇涓嶆竻锛岃鎷掔粷涓婅鍙版巿璇撅紝鍙彲鎵规敼浣滀笟锛屽悗琚滃仠鑱屽洖涔″吇鐥呪濓紝璋冨洖鍘﹂棬澶у浠昏祫鏂欏憳锛屽悓鏃剁爺绌舵暟璁猴紝瀵圭粍鍚堟暟瀛︿笌鐜颁唬缁忔祹绠$悊銆佺瀛﹀疄楠屻佸皷绔妧鏈佷汉绫荤敓娲荤殑瀵嗗垏鍏崇郴绛夐棶棰樹篃浣滀簡鐮旂┒銆1956骞磋皟鍏ヤ腑鍥界瀛﹂櫌鏁板鐮旂┒鎵銆1980骞村綋閫変腑绉戦櫌鐗╃悊瀛︽暟瀛﹂儴濮斿憳(鐜板湪鐨...
  • 寰垎鍑犱綍,澶у摜浠府甯繖
    绛旓細A I=dr^2>0
  • 扩展阅读:吉利几何e萤火虫缺点 ... 几何a为什么不建议买 ... 良辰好景知几何沈宴清 ... 几何e萤火虫值得买吗 ... 几何e为什么没人买 ... 微分几何答案完整版陈维桓 ... 几何c ... 几何画板手机版 ... 微分几何第二版陈维恒 ...

    本站交流只代表网友个人观点,与本站立场无关
    欢迎反馈与建议,请联系电邮
    2024© 车视网