[发明专利]基于全波形反演的医学成像方法、计算机设备及存储介质

说明书

本申请涉及一种基于全波形反演的医学成像方法、计算机设备及存储介质，本申请的技术方案通过确定初始速度模型；获取高频超声波换能器获取人体部位的实际观测数据；并根据所述实际观测数据和所述初始速度模型进行反演得到目标速度模型；最后基于所述目标速度模型进行人体部位的医学成像，把石油勘探中近年来热门FWI（Full Waveform Inversion）全波形反演的技术引入超声医学影像领域，采集的超声波数据，进行反演，得到目标介质的物理属性（例如速度、密度），该属性可以进行相关的医学分析，为目标介质的准确成像奠定基础。

技术领域

本申请涉及医学影像技术领域，特别是涉及一种基于全波形反演的医学成像方法、计算机设备及存储介质。

背景技术

现代医学影像领域，成像技术主要有核磁共振(MRI)、CT、超声波等手段，其中核磁共振成本高、设备体积庞大(动辄千万的造假，数百吨钢铁屏蔽强磁场，磁场铁芯也动辄数十吨乃至百吨)，对检测目标也有特殊要求(比如不能含有金属)；CT虽然精度高，但是辐射剂量大，设备体积大，成本也相对较高；超声波体积小，成本低，但是只是利用了简单的反射波成像，存在精度不足，无法穿透特殊的部位，以及对特殊部位的杂乱波场无法准确成像等缺点(例如头颅，超声波穿透能力弱，波场在头颅内来回震荡，难以识别反射和多次杂波干扰，无法精确成像。

发明内容

基于此，有必要针对上述技术问题，提供一种能够解决现有技术中精度不足问题的技术方案。为此，本申请提供了一种基于全波形反演的医学成像方法、计算机设备及存储介质。

本申请提供的基于全波形反演的医学成像方法，包括：一种基于全波形反演的医学成像方法，其特征在于，包括：确定初始速度模型；获取高频超声波换能器获取人体部位的实际观测数据；根据所述实际观测数据和所述初始速度模型进行反演得到目标速度模型；基于所述目标速度模型进行人体部位的医学成像。

可选地，所述根据所述实际观测数据和所述初始速度模型进行反演得到目标速度模型包括：获取由所述初始速度模型模拟产生正演数据，记录所有时刻的模拟波场数据U；计算相应接收点位置上模拟波场数据与实际观测数据的波场残差：Δd＝d_obs-d_mod，其中d_obs表示实际观测数据，d_mod表示模拟波场数据，Δd表示波场残差；根据所述波场残差计算瞬时梯度场；根据瞬时梯度场确定全区域梯度；根据所述全区域梯度更新初始速度模型以得到目标速度模型。

可选地，所述根据所述波场残差计算瞬时梯度场包括：将所述波场残差进行延拓反传，得到每一时刻残差波场：|d_obs-S×d_mod|，其中，S表示：振幅系数，d_obs表示实际观测数据，d_mod表示模拟波场数据；计算正传波场二阶偏导数将两个量做内积得到瞬时梯度：

可选地，所述根据瞬时梯度场确定全区域梯度包括：将所有时刻的瞬时梯度场叠加，得到单炮梯度；将所有单炮梯度叠加，得到全区域梯度。

可选地，根据所述全区域梯度更新初始速度模型以得到目标速度模型包括：将所述全区域梯度进行规则化；用规则化后的所述全区域梯度更新初始速度模型。

可选地，所述用规则化后的所述全区域梯度更新初始速度模型以后，还包括：判断波场残差是否满足更新后的初始速度模型的目标函数；如果满足，则输出更新后初始速度模型为目标速度模型。

可选地，所述获取由所述初始速度模型模拟产生正演数据包括：获取实际观测数据；对所述实际观测数据在初始速度模型上进行逆时偏移成像得到成像结果；对所述成像结果进行反偏移得到反偏移数据；使用所述反偏移数据作为所述正演数据。

本申请还提供了一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现权利要求以上任一项所述方法的步骤。