[发明专利]一种基于小波域正则化的相位重建方法

摘要

本发明涉及一种基于小波域正则化的相位重建方法,包括在待成像的物平面位置放置刀口器具，采集获得刀口图像，得到二维点扩散函数。放置成像物体,对物体成像，获得相衬成像结果。对相衬成像结果图进行2层小波变换。对各小波通道的信号分别进行傅里叶变换得到，并分别进行正则化处理，获得吸收项在小波通道的分量。将吸收项在各小波通道的分量分别进行傅里叶逆变换，得到空间域的各通道吸收项。进行小波逆变换，即可获得小波域正则化后的相衬图像吸收项。获得相位信息重建结果。

主权项

一种基于小波域正则化的相位重建方法,包括下列步骤：①X射线类同轴相衬成像参数设置：包括光源到物体的距离以及此对应的物体到探测器的距离；②设置数字放射成像系统的曝光参数，在待成像的物平面位置放置刀口器具，采集获得刀口图像，从图像中获取不同位置的刀口截面曲线，而后将刀口截面曲线进行平均，再对平均曲线求导数，可获得对应的表征系统特性的线扩散函数，将该线扩散函数在直角坐标系中旋转一周，计算得到二维点扩散函数是空间坐标，而后对其求傅里叶变换得到相衬成像系统的恶化传递函数为空间频率；③放置成像物体,对物体成像，获得相衬成像结果④对相衬成像结果图进行2层小波变换，得到2层共8个通道的小波变换结果a,b∈{L,H}，J＝1,2为小波分解层数；⑤对各小波通道的信号分别进行傅里叶变换得到而后对分别进行正则化处理，获得吸收项的在小波通道的分量这里的正则化参数αJab设定为各通道信噪比的倒数，表示对取共轭，为空间频率,其中为相衬成像系统的恶化传递函数，λ为X射线波长，re＝2.82×10‑15m为经典电子半径，σKN为康普顿散射截面，M为成像放大倍数，R2为物体到探测器的距离；⑥将吸收项在各小波通道的分量分别进行傅里叶逆变换，得到空间域的各通道吸收项综合所有通道的吸收项进行小波逆变换，即可获得小波域正则化后的相衬图像吸收项⑦利用公式即可获得相位信息重建结果