发明名称 |
基于全广义变分正则化的低剂量X射线CT图像重建方法 |
摘要 |
一种基于全广义变分正则化的低剂量X射线CT图像重建方法,包括如下步骤:(1)获取CT设备的系统参数和低剂量扫描协议下对数变换后的投影数据;(2)对投影数据进行Anscombe变换,将服从Poission分布的投影数据转化为近似服从方差为1的Gaussian分布数据u;(3)建立基于全广义变分最小化的理想数据恢复模型,使用Chambolle‑Pock算法求解得到恢复后的投影数据f;(4)对步骤(3)得到的恢复后的投影数据f进行Anscombe逆变换,再通过滤波反投影算法得到CT重建图像。该方法可以在投影数据不满足分段常数假设的前提下去除图像的噪声和条形伪影,同时较好地保持图像的空间分辨率。 |
申请公布号 |
CN105844678A |
申请公布日期 |
2016.08.10 |
申请号 |
CN201610427491.7 |
申请日期 |
2016.06.15 |
申请人 |
赣南师范学院 |
发明人 |
牛善洲;李楠;吴恒;马建华;喻高航 |
分类号 |
G06T11/00(2006.01)I |
主分类号 |
G06T11/00(2006.01)I |
代理机构 |
北京科亿知识产权代理事务所(普通合伙) 11350 |
代理人 |
赵蕊红 |
主权项 |
一种基于全广义变分正则化的低剂量X射线CT图像重建方法,其特征在于,包括如下步骤:(1)获取CT设备的系统参数和低剂量扫描协议下对数变换后的投影数据q;(2)对步骤(1)获取的投影数据q进行Anscombe变换,将服从复合Poission分布的投影数据q转化为近似服从方差为1的Gaussian分布数据u;(3)对步骤(2)获取的数据u建立基于全广义变分最小化的理想数据恢复模型,使用Chambolle‑Pock算法求解得到恢复后的投影数据f;(4)对步骤(3)得到的恢复后的投影数据f进行Anscombe逆变换,再通过滤波反投影算法得到CT重建图像。 |
地址 |
341000 江西省赣州市经济技术开发区赣南师范学院 |