数字图象处理Chaper5 图像重建.pptVIP

滤波反投影重建图像 Original Image A Radon Transform B 0o 179o MATLAB Example of Computed Tomography using Filtered Back-Projection (FBP) Recovered image C B=radon(A,0:179); C=iradon(B,0:179,nearest, Hann); load trees I=ind2gray(X,map); A=imcrop(I,[1,1,128,128]); figure; imshow (A); B=radon(A,0:179); BI=ind2gray(B,map); figure; imshow (BI,128) BS=B(1:128,1); figure; plot(BS) C=iradon(B,0:179,nearest, Hann); figure; imshow (C) MATLAB program for Projection Back Projection Load an existing color image Transform it into gray image Crop it Display it Radon transform it Transform into gray image Display projections as an image Extract one projection Plot one projection Inverse Radon Display recovered image 投影滤波重建 使用扇形射线束滤波反投影重建 基本的扇形射线几何原理 使用扇形射线束滤波反投影重建 扇形射线上任意一点的极坐标表示 扇形束滤波反投影 Algebraic Reconstruction Techniques (ART) Iterative Reconstruction with ART 图像和Radon变换 正弦图 f(x,y) 0o 179o MATLAB Example of the Radon Transform R g(s,q) x y q s g(s,0) s B=radon(A,0:179); Full circle in RT Shifted full circle in RT Empty circle in RT Shifted empty circle in RT Thin stick in RT Shifted thin stick in RT Full triangle in RT Shifted full triangle in RT Full square in RT Shifted full square in RT Empty square in RT Shifted empty square in RT | x |2/3 + | y |2/3 ≤ 1 in RT | x | + | y | ≤ 1 in RT | x |3/2 + | y |3/2 ≤ 1 in RT | x |2 + | y |2 ≤ 1 in RT | x |6 + | y |6 ≤ 1 in RT | x |n + | y |n ≤ 1 for n ? ? in RT 2D Gaussian in RT Shifted 2D Gaussian in RT 从正弦图到图像：反投影法 原始图像 正弦图（投影数据） 反投影重建 傅里叶切片定理（Projection-Slice Theorem） If g(s,q) is the Radon transform of f(x,y), then G(ns ,q) = F(ns cosq , ns sinq) G(ns,q) = 1D Fourier transform of g(s,q) w.r.t. s F(nx,ny) = 2D Fourier transform of f(x,y) x s q s g(s,q) q nx ny The 1D FT of the projection of an image = Slice of the 2D FT of the image F(nx,ny) f(x,y) y Theory of pure RT and IRT Radon transform Inverse Radon transform 斜坡滤波器。不可积？ 从正弦图到图像：滤波反投影法 滤波反投影法：滤波器设计 加窗限制带宽的斜坡滤波器（降低模糊和振铃效应） 滤波反投影重建图像 使用斜坡滤波器 使用汉明窗斜坡滤波器 第5章 图像复原与重建Image Resto