CT平行束和扇形束算法的转换培训课件.pptVIP

;平行束图像重建;中心切片定理: 二维图像的一维投影（线积分）的傅立叶变换，恰好等于图像本身的二维傅立叶变换的一个特定截面。; 通过转动投影方向，可以得到各个方向上傅立叶变换的特定截面，从而获得整个二维平面的傅立叶变换，最后由傅立叶逆变换得到重建的图像;平行束图像重建;将直角坐标（u，v）转换成极坐标（w, θ）.;具体的推导;交换积分次序：;扇形束图像重建;对与扇形束成像，我们并没有相应的中心切片定理。我们只好想个别的办法来推导扇形束的图像重建算法。这个办法就是把扇形束的成像问题转化成平行光束的成像问题，把平行光束图像重建的算法修正一下然后应用于解决扇形束的成像问题中。;11;12;等角度扇形重建算法; 已知平行束重建算法为：  转化为极坐标后得： ;利用 将平行束的变量dsdθ换成扇形束的变量dγdβ， 其中雅克比因子为 ，这样可得;16;斜坡滤波器卷积核的一个特殊性质：;18; 假如我们现在想重建某点 的值，我们先确定一个β，即确定了源的位置，由于重建点的位置也是确定的，故D?和γ?均为确定的值。然后用卷积核 对不同γ角的信号进行卷积。 当这个 角度的卷积完成后，我们对0°-180°所有角度的β做 一个积分，即背投影过程。 这就是扇形束的滤波背投影算法。;短扫描;根据同样的道理，当扇形束探测器旋转2π, 每条投影射线也都被测到了两次。冗余的数据可由下面这个表达式给出    由于数据冗余，在扇形束数据采集中没有必要让探测器做2π 全扫描。扫描角度(β) 可以小于2π, 这种扫描方式叫做短扫描。; 角度β 的最小取值范围取决于数据采集时物体与探测器之间的几何关系。角度β的最小取值区间可能小于π (下图左)，可能等于π (中)，也可能大于π ( 右)。确定扫描区间的原则是，我们感兴趣的物体中的每一点都要有180°的角度覆盖。  要注意的是，在扇形束短扫描成像中，并不是所有的线积分都被刚好测到一次。有些线积分被测到一次，而另一些线积分会被测到两次。即使在???描角度β 的范围小于π 的情形, 还是有一些线积分会被测到两次的。其实，任何直线，只要它与扇形的焦点轨迹有两个交点，那么沿这条线的线积分就被测到了两次(图3.10)。;  其实，任何直线，只要它与扇形的焦点轨迹有两个交点，那么沿这条线的线积分就被测到了两次。为了获得足够的数据来做断层成像，我们要求过物体的每一条线的线积分都要至少被测到一次。对于冗余的数据，在图像重建时需要施加适当的权函数来处理。举例来说，若一个线积分被测到了两次，对它们要进行加权，而且它们的权因子之和一定要是1。;数据冗余条件;主要参考资料： [1] CT重建算法；闫镔，李磊 编著；2014，科学出版社 [2] CT重建算法；曾根生；2009年，高等教育出版社;平行束图像重建;中心切片定理: 二维图像的一维投影（线积分）的傅立叶变换，恰好等于图像本身的二维傅立叶变换的一个特定截面。; 通过转动投影方向，可以得到各个方向上傅立叶变换的特定截面，从而获得整个二维平面的傅立叶变换，最后由傅立叶逆变换得到重建的图像;平行束图像重建;将直角坐标（u，v）转换成极坐标（w, θ）.;具体的推导;交换积分次序：;扇形束图像重建;对与扇形束成像，我们并没有相应的中心切片定理。我们只好想个别的办法来推导扇形束的图像重建算法。这个办法就是把扇形束的成像问题转化成平行光束的成像问题，把平行光束图像重建的算法修正一下然后应用于解决扇形束的成像问题中。;11;12;等角度扇形重建算法; 已知平行束重建算法为：  转化为极坐标后得： ;利用 将平行束的变量dsdθ换成扇形束的变量dγdβ， 其中雅克比因子为 ，这样可得;16;斜坡滤波器卷积核的一个特殊性质：;18; 假如我们现在想重建某点 的值，我们先确定一个β，即确定了源的位置，由于重建点的位置也是确定的，故D?和γ?均为确定的值。然后用卷积核 对不同γ角的信号进行卷积。 当这个 角度的卷积完成后，我们对0°-180°所有角度的β做 一个积分，即背投影过程。 这就是扇形束的滤波背投影算法。;短扫描;根据同样的道理，当扇形束探测器旋转2π, 每条投影射线也都被测到了两次。冗余的数据可由下面这个表达式给出    由于数据冗余，在扇形束数据采集中没有必要让探