CT-算法重建-马建-Lesson01-CT-Recon.pptVIP

小结 (Summarize) CT图像重建的历史（ A brief history ） 基本术语（Some basic terms） 图像重建思想（Reconstruction ideas） 图像重建算法（Reconstruction algorithms) 作业(Assignment) 什么叫做投影（projection）？给出某一物体f(x,y)在θ角度投影的数学表达式？ 什么是Radon变换？ 什么是弦图（sinogram）？正常情况下弦图的形状？在第三代CT中如果某个探测器损坏，其弦图应该是什么样的？ 傅立叶切片定理是怎样表述的？并证明之。(水平或垂直方向) 5. 用迭代法重建下列矩阵 2 5 3 4 上机作业 一. 基本概念操作 Look for phantom, radon and iradon i Matlab Help在Matlab上学习Phantom/radon/iradon函数； Make a phantom size (128,128)设计一个体模，大小为128x128； Do radon transformation (CT projections) with theta1 = (0:10:180) and thera2 = (0:1:180)分别以间隔10和间隔1度做radon； Notice the sinograms.显示各自的弦图； 上机作业 二. 在matlab上运行下列命令，给出注释 ph = phantom(256); figure,imshow(ph,[]) ； theta = 0:0.5:179; sino = radon(ph,theta); figure,imshow(sino,[]); * * * * * * * * * * * * * * * * * * * * 由弦图重建图像 弦图 (已知) CT图像 (待求) 直接矩阵求解法（Direct matrix inversion） 迭代法（ iterative） 傅里叶重建法（Fourier reconstruction） 反投影法（Back-projection） 滤波反投影法（Filtered back projection） Algorithms 直接矩阵求解法 2x2 matrix A11 A12 Po (0) = 3 2 A21 A22 Po (1) = 7 4 P9o (0) P9o (1) =4 =6 A11 + A12 = 3 A21 + A22 = 7 A11 + A21 = 4 A12 + A22 = 6 1 1 0 0 A11 3 0 0 1 1 A12 7 1 0 1 0 A21 4 1 0 0 1 A22 6 F * A = P A = P F^(-1) = 设待重建图像大小为：28 *28 = 256 * 256， 则，矩阵F的大小为： 216 * 216 = 232 因此，直接求矩阵的反变换很难实现 直接矩阵求解法的缺点 迭代法 给出初始矩阵 用初始矩阵形成投影 把待重建物体投影与模拟投影进行比较 若误差满足要求，迭代停止 迭代法的基本步骤 迭代法的一个实列 比较精确，但速度慢，多用于核医学设备的图像重建和低剂量CT的图像重建 迭代法的缺点 傅立叶变换重建方法 傅立叶变换 傅立叶切片定律及证明 傅立叶重建的基本步骤 傅立叶变换 傅立叶切片定律 一个物体的1D投影的傅立叶变换精确地等于物体2D傅立叶变换在同一角度的直线。换言之：图像沿某一方向的投影，经过1D 傅立叶变换之后，对应2D傅立叶变换平面的一条线 傅立叶切片定律的示意图 傅立叶切片定律的示意图 傅立叶切片定律的示意图 傅立叶切片定律的示意图 由投影到图像 x y ? g(s,?) μ(x,y) wx wy ? s 1D-Fourier transform F[g(s)] 中心切片定理(Central Slice Theorem) 1D FT F[g(s)] 2D FT F[μ(x,y)] s Radon transform (i.e., measurement) at view angle 90? g(s, 90?? wx wy wx wy 2D spatial domain of μ(x,y) 2D frequency domain of μ(x,y) 中心切片定理(Central Slice Theorem) x y ? g(s,?) μ(x,y) wx wy ?