[发明专利]迭代间曲率流滤波加权最小二乘正电子发射断层成像方法

摘要

本发明公开了一种用于构建成像图像的迭代间曲率流滤波加权最小二乘正电子发射断层成像方法，首先获取观测数据，确定初始图像的规模，将观测数据中的各个分量进行逐个与1比较，与误差方差中的各个分量逐个对应相除等后将中间图像进行曲率流滤波，可得到作为初始图像的滤波图像，重复这个过程直到重建后的图像收敛，可得到最终成像的图像；本发明结合当前PET成像数据的特点，在现有最小二乘方法迭代过程中加入曲率流滤波机制，以此来抑制图像噪化，提高成像图像的质量，实现较为简单，处理速度较快，对抑制重建图像的背景干扰也有着良好的效果。

主权项

1.一种用于构建成像图像的迭代间曲率流滤波加权最小二乘正电子发射断层成像方法，其特征在于采用下列步骤：1)获取观测数据，记录其规模M，并将它保存在向量y＝[y₁，Λ，y_M]^T中；根据待重建图像的尺寸要求，确定初始图像的规模N，给定它的初始灰度值大于1，并将其保存在向量f＝[f₁，Λ，f_N]^T中，2)根据观测数据y和初始图像f的大小，计算M×N的系统概率矩阵A，3)将观测数据y中的各个分量进行逐个与1比较：如果比1大则保留，比1小则设置为1，并将比较结果保存在向量d中，d＝[d₁，Λ，d_M]^T，作为观测数据的误差方差，4)将观测数据y中的各个分量与误差方差d中的各个分量逐个对应相除得到修正的观测投影数据y~=[y~1,Λ,y~M]T,]]>其中y~i=yi/di,]]>i＝1，Λ，M，5)将系统概率矩阵A的转置与修正的观测数据相乘得到向量g，g=ATy~,]]>g＝[g₁，Λ，g_N]^T，6)将系统概率矩阵A和初始图像f相乘，得到前向投影p，p＝[p₁，Λ，p_M]^T，再将该前向投影p中的各个分量与误差方差d中的各个分量逐个对应相除得到修正的前向投影p~=[p~1,Λ,p~M]T,]]>其中p~i=pi/di,]]>i＝1，Λ，M，然后利用系统概率矩阵A的转置与该前向投影相乘得到反投影h：h=ATp~,]]>h＝[h₁，Λ，h_N]^T，7)将向量g的各个分量与反投影h中的各个分量逐个对应相除得到图像修正值c：c＝[c₁，Λ，c_N]^T，其中c_j＝g_j/h_j，j＝1，Λ，N，8)将该图像修正值c中的各个分量与初始图像f中的各个分量逐个对应相乘得到中间图像f^=[f^1,Λ,f^N]T,]]>其中f^j=cj×fj,]]>j＝1，Λ，N，9)将中间图像进行至少一次曲率流滤波，可得到滤波图像10)将这个滤波图像作为初始图像，返回到第6步，重复这个过程直到重建后的图像收敛，可得到最终成像的图像。