发明名称 多相水平集的正电子断层扫描重建方法
摘要 本发明公开了一种多相水平集的正电子断层扫描重建方法,先获取投影数据,再确定水平集函数的个数,用K均值聚类的方法得到图像初始浓度值,经计算得到一个用于重建图像的目标函数,然后将该目标函数对每一个像素求偏导数并将初始图像中像素点分别对水平集函数求偏导数,将以上得到的两个偏导数相乘得到当前水平集进化的修正值,将水平集函数进化后作为下次迭代时的初始水平集,更新重建图像,得到此次迭代的估计图像,并作为下次迭代时的初始图像,对每个区域求均值,得到浓度估计值,并作为下次迭代时的浓度估计值,再循环以上步骤,直到重建后的图像收敛,本发明具有消除噪声,保持边缘,消除边缘伪影等优点。
申请公布号 CN1640361A 申请公布日期 2005.07.20
申请号 CN200510037622.2 申请日期 2005.01.06
申请人 东南大学 发明人 朱宏擎;周键;舒华忠;罗立民;李松毅
分类号 A61B6/03;A61B8/13 主分类号 A61B6/03
代理机构 南京经纬专利商标代理有限公司 代理人 陆志斌
主权项 1.一种多相水平集的正电子断层扫描重建方法,其特征在于采用下列步骤:1)获取投影数据,对获得的断层面投影数据用现有的断层扫描(CT)用的滤波反投影(FBP)方法进行重建,得到一个初始的图像,选取水平集进化的步长Δt,2)根据初始图像的浓度组成,确定水平集函数的个数,3)用K均值聚类的方法得到图像的初始浓度值,4)对现有的加权最小二乘方法的目标函数加上经离散后的初始图像的全变分的β倍,且该β为0~1之间中的任意一个数,得到一个用于重建图像的目标函数,再将该目标函数对每一个像素求偏导数,5)将初始图像中每一个像素点分别对每一个水平集函数求偏导数,6)将步骤4)和步骤5)的两个偏导数相乘,得到当前水平集进化的修正值7)对每个水平集函数,将初始水平集函数减去步长Δt倍的水平集进化修正值得到进化后的水平集函数,并将这些进化后的水平集函数作为下次迭代时的初始水平集,8)利用海信函数,更新重建图像,得到此次迭代的估计图像,并将此次迭代的估计图像作为下次迭代时的初始图像,并对此次迭代的估计图像中的每个区域求均值,得到每个区域的浓度估计值,并将这些浓度估计值作为下次迭代时的浓度估计值,再返回到第4步,直到重建后的图像收敛。
地址 210096江苏省南京市四牌楼2号