用于超声环阵成像的相位畸变校正方法和系统.pdf

《用于超声环阵成像的相位畸变校正方法和系统.pdf》由会员分享，可在线阅读，更多相关《用于超声环阵成像的相位畸变校正方法和系统.pdf（18页完成版）》请在专利查询网上搜索。

1、(19)国家知识产权局(12)发明专利申请(10)申请公布号 (43)申请公布日 (21)申请号 202410005094.5(22)申请日 2024.01.03(71)申请人 之江实验室地址 311121 浙江省杭州市余杭区中泰街道科创大道之江实验室(72)发明人 林宏翔荣超蓝郑枫吴光未李劲松(74)专利代理机构 杭州天正专利事务所有限公司 33201专利代理师 王兵杨东炜(51)Int.Cl.A61B 8/00(2006.01)A61B 8/14(2006.01)G06T 5/80(2024.01)G06T 3/4053(2024.01)(54)发明名称一种用于超声环阵成像的相位畸变校正方。

2、法和系统(57)摘要一种用于超声环阵成像的相位畸变校正方法和系统，其方法包括：设置超声环阵的阵元激发和接收模式，采集超声回波信号；对各通道接收数据进行数据预处理，将拾取的直达波信号区间赋值为零；设置先验的声速范围来计算不同声速下的相干因子，得到每个成像点在每个先验声速下对应的相干因子值；将得到的相干因子值分别对水平和垂直方向进行插值，再根据最大化相干因子分别得到水平和垂直分布的声速；将水平和垂直分布的声速相加和平滑处理，得到最终的声速分布；根据成像点位置、激发阵元位置以及接收阵元位置的几何距离和得到的声速分布来计算延迟时间，通过延迟时间将各通道信号进行对齐；将对齐的各通道信号进行聚焦，得到高分。

3、辨率超声成像结果。权利要求书3页 说明书10页 附图4页CN 117503203 A2024.02.06CN 117503203 A1.一种用于超声环阵成像的相位畸变校正方法，其特征在于，包括如下步骤：S1、设置超声环阵的阵元激发和接收模式，实现超声回波信号采集；S2、对各通道接收数据进行数据预处理，利用切除法将拾取的直达波信号区间赋值为零；S3、设置先验的声速范围来计算不同声速下的相干因子，得到每个成像点在每个先验声速下对应的相干因子值；S4、将得到的相干因子值分别对水平和垂直方向进行插值，再根据最大化相干因子来分别得到水平和垂直分布的声速；S5、将水平和垂直分布的声速相加和平滑处理，得到最。

4、终的声速分布；S6、根据成像点、激发阵元以及接收阵元的几何距离和得到的声速分布来计算延迟时间，通过延迟时间将各通道信号进行对齐；S7、使用延迟叠加(DAS)波束形成方法将对齐的各通道信号进行聚焦，从而得到高分辨率超声成像结果。2.根据权利要求1所述的方法，其特征在于，在步骤S1中，设置环阵换能器每次仅有一个阵元激发，所有阵元接收信号，该过程将重复进行直至所有阵元都依次激发完毕，在这种模式下，激发次数等于环阵换能器的阵元数目。3.根据权利要求2所述的方法，其特征在于，在步骤S1中，利用超声设备进行数据采集时，在每个阵元激发之后，在相邻阵元激发时需要有一定的时间间隔，以确保环阵内的超声波衰减到可忽。

5、略的水平。4.根据权利要求1所述的方法，其特征在于，在步骤S2中，对环阵换能器各通道接收到的回波信号进行预处理，主要去除激发阵元直接传播到接收阵元位置的直达波信号；数据预处理根据直达波信号能量强，在回波数据曲线上振幅大的特征，来确定直达波截取区间范围；通过找到接收回波信号上最大值位置，再将直达波区间范围信号归零。5.根据权利要求1所述的方法，其特征在于，步骤S3中，根据成像介质设定先验的声速范围为14601620m/s，声速增量为1m/s；根据已知成像点位置、激发阵元位置以及接收阵元位置的几何距离和设定的先验声速计算延迟时间，从各通道信号中提取延迟信号，此时每个成像点处的声速均为恒定声速；延迟。

6、时间 由公式(1)得到：式中，表示激发阵元位置，表示 序号接收阵元位置，表示坐标成像点位置，表示坐标成像点处的声速。6.根据权利要求5所述的方法，其特征在于，步骤S3中，一个激发阵元对应多个接收阵元，故能计算出多组延迟时间，从而能提取到多个延迟信号；将从每个阵元提取的延迟信号代入下述公式(2)，计算该声速下的相干因子：式中，表示坐标成像点位置的相干因子，N代表接收阵元的数目，权利要求书1/3 页2CN 117503203 A2代表具体激发阵元序号代表具体接收阵元序号，表示坐标成像点位置在序号为m的激发阵元在序号为n接收阵元上提取的延迟信号，为序号m的激发阵元经过坐标成像点位置到达序号为n接收阵。

7、元上的时间。7.根据权利要求1所述的方法，其特征在于，步骤S4中，为了进一步提高求解声速精度，在相干因子和声速分布上做一维插值处理；考虑到直接求解全局声速会得到离散的声速分布，于是对每个声速下的二维相干因子矩阵分别在水平方向和垂直方向上求均值；再根据最大化相干因子原则，每个成像点位置选择相干因子最大值时对应的声速，从而得到水平和垂直分布的声速。8.根据权利要求1所述的方法，其特征在于，步骤S5中，将得到水平和垂直分布的声速进行相加求平均，再将得到的声速做进一步的平滑处理；进一步平滑策略主要是进行9点平均加权滑动窗口平滑方法，见下述公式(3)：式中，表示坐标()成像点位置的声速。9.根据权利要求。

8、1所述的方法，其特征在于，步骤S6中，利用步骤S5中得到的声速分布代替恒定声速计算延迟时间，通过延迟时间将各通道信号进行对齐。10.一种用于超声环阵成像的相位畸变校正系统，其特征在于，包括：数据采集模块、数据预处理模块、估计声速分布模块、相位校正模块、超声成像模块以及图像可视化模块，其中，数据采集模块，用于设置超声环阵的阵元激发和接收模式，利用接收阵元记录接收到的回波信号；数据预处理模块，用于对各通道接收数据进行预处理，主要是利用切除法将拾取的直达波信号区间赋值为零；估计声速分布模块，用于将先验声速范围和计算相干因子联系起来，从而得到每个成像点位置的声速分布；相位校正模块，用于已知成像点位置、。

9、激发阵元位置以及接收阵元位置的几何距离和每个成像点位置的声速分布计算校正后的延迟时间，通过计算得到的延迟时间将各通道信号进行对齐；超声成像模块，用于将对齐的各通道信号使用延迟叠加(DAS)波束形成方法进行聚焦，从而得到高分辨率超声成像结果；图像可视化模块，用于显示相位校正后的高分辨率超声成像结果。11.一种计算机可读存储介质，其特征在于，其上存储有程序，该程序被处理器执行时，实现权利要求19中任一项所述的一种用于超声环阵成像的相位畸变校正方法。12.一种用于超声环阵成像的相位畸变校正装置，其特征在于，包括存储器和一个或多个处理器，所述存储器中存储有可执行代码，所述一个或多个处理器执行所述可执行。

10、代码权利要求书2/3 页3CN 117503203 A3时，用于实现权利要求19中任一项所述的一种用于超声环阵成像的相位畸变校正方法。权利要求书3/3 页4CN 117503203 A4一种用于超声环阵成像的相位畸变校正方法和系统技术领域0001本发明属于超声成像技术，具体涉及一种用于超声环阵成像的相位畸变校正方法和系统。背景技术0002超声断层扫描成像技术是一种无电离辐射、高帧率的医学影像模态，它通过采集生物体软组织的透射和散射超声波形信号，将其反演为组织结构图像。超声成像以实时、安全、低成本等优点被广泛应用于临床。现有的临床超声成像探头是线阵或凸阵换能器，它们具有便携性以及更具成本效益的优。

11、势。但由于线阵或凸阵视野的受限，检测生物体组织的分辨率是各向异性的，横向分辨率低于轴向分辨率。其次，图像质量在很大程度上取决于生物体内部组织的深度，因为来自深度区域的回波强度取决于沿传播路径衰减的变化。此外，扫描图像的质量高度依赖于超声科医生的经验和培训水平。环阵换能器拥有360 的视野角度，超声波束可以聚焦在每个点上。此外，环阵换能器内的任何区域都有靠近该区域的检测元素，因此可以抑制图像质量的深度依赖性。环阵换能器不仅可以同时接收不同方向的回波，还可以同时接收散射波和透射波，利用超声环阵成像技术可以获得生物体组织的结构信息。0003现有的超声环阵成像技术通过标准的延时叠加（delayands。

12、um，DAS）波束形成方法得到超声图像，但由于延时叠加方法无法区分主瓣和副瓣的信号，且副瓣的信号含有干扰波和杂波，所成图像质量不佳尚未达到临床使用的标准。通过假设传播介质是各向同性的，根据恒定声速和已知激发阵元、成像点和接收阵元间的几何距离计算各阵元时间延迟，因此，恒定声速的选取误差影响最终的图像质量。在实际超声环阵成像时，通常将成像物体置于水中，假设水温保持不变且水温不会受到检测物体的影响，将水的声速近似为一个恒定声速，这与实际物理场景有差异而导致建模的系统性误差。一般情况下使用恒定声速可以得到较好的成像分辨率，但是当重建生物体不同组织间声速差异大、非均质性高时，文献（Pinton G F,。

13、Trahey G E,Dahl J J.Sources of image degradation in fundamental and harmonic ultrasound imaging using nonlinear,fullwave simulationsJ.IEEE transactions on ultrasonics,ferroelectrics,and frequency control,2011,58(4):754765.）指出利用恒定声速进行波束形成方法时图像会出现局部结构变形、重影等模糊情况，造成超声图像分辨率低。文献（Kim H,Song I,Kang J,et al.。

14、Phase aberration correction for ultrasound imaging guided extracorporeal shock wave therapy (ESWT):Feasibility studyJ.Ultrasonics,2023,132:107011.）中对比了不同恒定声速下的波束形成结果，通过视觉评价发现，若选取的恒定声速与真实声速的误差大时，重建图像的整体质量严重下降，成像目标的物理位置以及界面出现了明显扭曲。反之，当恒定声速误差小时，成像目标的空间分辨率得到提高，但随着成像深度的增加，重建图像的相位像差影响逐渐增大，图像质量和目标位置也随着深度的增。

15、加而逐渐失真。说明书1/10 页5CN 117503203 A50004现有技术的缺点表现为：00051）在实际情况下，通过水的声速近似得到恒定声速，这与实际物理场景有差异而导致建模的系统性误差，将导致通道信号的失准，造成超声成像的分辨率降低；00062）在遇到复杂介质或声速差异较大时，尤其当作用于存在显著的声速不均匀性介质时，则在利用恒定声速计算几何延迟后，高阶相位畸变会导致信号失准，随着成像深度的增加，重建图像的相位像差影响逐渐增大，超声图像出现局部结构变形、重影等分辨率降低情况。发明内容0007本发明要解决现有超声环阵成像技术的问题，提出一种用于超声环阵成像的相位畸变校正方法和系统。00。

16、08本发明提出的一种用于超声环阵成像的相位畸变校正方法，包括以下步骤：0009S1，设置超声环阵的阵元激发和接收模式，实现超声回波信号采集；0010S2，对各通道接收数据进行数据预处理，利用切除法将拾取的直达波信号区间赋值为零；0011S3，设置先验的声速范围来计算不同声速下的相干因子，得到每个成像点在每个先验声速下对应的相干因子值；0012S4，将得到的相干因子值分别对水平和垂直方向进行插值，再根据最大化相干因子来分别得到水平和垂直分布的声速；0013S5，将水平和垂直分布的声速相加和平滑处理，得到最终的声速分布；0014S6，根据成像点、激发阵元以及接收阵元的几何距离和得到的声速分布来计算。

17、延迟时间，通过延迟时间将各通道信号进行对齐；0015S7，使用延迟叠加(DAS)波束形成方法将对齐的各通道信号进行聚焦，从而得到高分辨率超声成像结果；0016进一步，在步骤S1中，设置环阵换能器每次仅有一个阵元激发，所有阵元接收信号，该过程将重复进行直至所有阵元都依次激发完毕，在这种模式下，激发次数等于环阵换能器的阵元数目。0017进一步，在步骤S1中，利用超声设备进行数据采集时，在每个阵元激发之后，在相邻阵元激发时需要有一定的时间间隔，以确保环阵内的超声波衰减到可忽略的水平。0018进一步，在步骤S2中，对环阵换能器各通道接收到的回波信号进行预处理，主要去除激发阵元传播到接收阵元位置的直达波。

18、信号。数据预处理根据直达波信号能量强，在回波数据曲线上振幅大的特征，来确定直达波截取区间范围。通过找到接收回波信号上最大值位置，再将直达波区间范围信号归零。0019进一步，在步骤S3中，根据成像介质设定先验的声速范围为14601620m/s，声速增量为1m/s。根据已知成像点位置、激发阵元位置以及接收阵元位置的几何距离和设定的先验声速计算延迟时间，从各通道信号中提取延迟信号，此时每个成像点处的声速均为恒定声速。延迟时间由公式(1)得到：说明书2/10 页6CN 117503203 A600200021式中，表示序号激发阵元位置，表示序号接收阵元位置，表示坐标成像点位置，表示坐标成像点处的声速。。

19、0022进一步，在步骤S3中，一个激发阵元对应多个接收阵元，故能计算出多组延迟时间，从而能提取到多个延迟信号。将从每个接收阵元提取的延迟信号代入下述公式(2)，计算该声速下的相干因子：00230024式中，表示坐标成像点位置的相干因子，N表示接收阵元的总数目，表 示 具体 激 发阵 元 序 号，表 示 具体 接 收阵 元 序 号，表示坐标成像点位置在序号为m的激发阵元在序号为n接收阵元上提取的延迟信号，为序号m的激发阵元经过坐标成像点位置到达序号为n接收阵元上的时间。0025进一步，在步骤S3中，对设置所有先验声速进行循环求解，在每个成像点位置均可以得到每个声速对应的相干因子值。0026进一步。

20、，在步骤S4中，为了进一步提高求解精度，在相干因子和声速分布上做一维插值处理。考虑到直接求解全局声速会得到离散的声速分布，于是对每个声速下的二维相干因子矩阵分别在水平方向和垂直方向上求均值。0027进一步，在步骤S4中，根据最大化相干因子原则，每个成像点位置选择相干因子最大值时对应的声速，从而得到水平和垂直分布的声速。0028进一步，在步骤S5中，将得到水平和垂直分布的声速进行相加求平均，再将得到的声速做进一步的平滑处理。进一步平滑策略主要是进行9点平均加权滑动窗口平滑方法，见下述公式(3)00290030式中，表示坐标成像点位置的声速。0031进一步，在步骤S6中，利用步骤S5中得到的声速分。

21、布，在每个成像点位置，利用声速分布代替恒定声速代入公式(1)计算延迟时间。说明书3/10 页7CN 117503203 A70032本发明的第二个方面涉及一种用于超声环阵成像的相位畸变校正系统，包括数据采集模块、数据预处理模块、估计声速分布模块、相位校正模块、超声成像模块以及图像可视化模块，其中，0033数据采集模块，用于设置超声环阵的阵元激发和接收模式，利用接收阵元记录接收到的回波信号；0034数据预处理模块，用于对各通道接收数据进行预处理，主要是利用切除法将拾取的直达波信号区间赋值为零；0035估计声速分布模块，用于将先验声速范围和计算相干因子联系起来，从而得到每个成像点位置的声速分布。0。

22、036相位校正模块，用于已知成像点位置、激发阵元位置以及接收阵元位置的几何距离和每个成像点位置的声速分布计算校正后的延迟时间，通过计算得到的延迟时间将各通道信号进行对齐；0037超声成像模块，用于将对齐的各通道信号使用延迟叠加(DAS)波束形成方法进行聚焦，从而得到高分辨率超声成像结果；0038图像可视化模块，用于显示相位校正后的高分辨率超声成像结果。0039本发明的第三个方面涉及一种计算机可读存储介质，其上存储有程序，该程序被处理器执行时，实现本发明的一种用于超声环阵成像的相位畸变校正方法。0040本发明的第四个方面涉及一种用于超声环阵成像的相位畸变校正装置，包括存储器和一个或多个处理器，所。

23、述存储器中存储有可执行代码，所述一个或多个处理器执行所述可执行代码时，用于实现本发明的一种用于超声环阵成像的相位畸变校正方法。0041本发明的创新点是：00421)相比于现有的恒定声速估计方法，本发明利用回波数据计算在不同声速下延时信号的相干因子，使用插值、相干因子最大化原则以及平滑方法求解每个成像点位置的声速分布。00432)相比于现有超声环阵波束形成方法技术，本发明利用声速分布来代替恒定声速来计算延迟时间，可减少由实际物理建模中产生的系统性误差。0044本发明基于所有合成孔径通道数据和相干因子来估计声速分布，用声速分布代替恒定声速计算各阵元时间延迟，通过标准的延时叠加波束形成方法，解决现有。

24、超声环阵成像技术在遇到生物体组织差异大或非均质性高时产生的相位畸变问题，从而提高超声图像的分辨率。0045本发明的有益效果是：与人为测量恒定声速方法相比，本发明得到声速分布基于所有合成孔径通道数据，得到的声速分布可以减少与实际物理场景有差异而导致建模的系统性误差。进一步的，利用声速分布来代替恒定声速，再结合已知几何距离来计算延迟时间，可以解决在遇到复杂结构和声速不均匀介质时各通道信号不能很好的聚焦问题，可以有效避免信号失准，从而提高超声图像的分辨率。另外，本发明是估计每个成像点位置上的声速分布，与采样时间间隔无关，可适用于实际不同采样率的采集数据。附图说明0046图1是本发明方法的流程图；说明。

25、书4/10 页8CN 117503203 A80047图2是本发明的数据预处理模块的流程图；0048图3是本发明估计声速分布模块的流程图；0049图4(a)图4(b)是本发明测试例中的一组结果展示：图4 (a)为基于恒定声速的超声成像结果；图4 (b)为基于相位校正的超声成像结果；0050图5是本发明的系统架构图。具体实施方式0051下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。0052需。

26、要说明的是，在不冲突的情况下，下述的实施例及实施方式中的特征可以相互组合。0053实施例10054本实施例提出一种用于超声环阵成像的相位畸变校正方法，实施例选用公开数据集OAbreast乳腺数据的二维切片，基于matlab平台和kwave工具箱进行数值仿真。图1是该方法的整体流程图，具体包括以下步骤：0055步骤S1，设置超声环阵的阵元激发和接收模式，实现超声回波信号采集。0056在本实施例中，具体的仿真参数如下：超声环阵阵元数目为512，超声环阵半径为40mm，超声波激发的中心频率为3MHz，矩形计算域边长，计算网格点数为。设置超声环阵换能器每次一个阵元激发所有阵元接收信号，该过程将重复进行。

27、直至所有阵元都依次激发完毕。在这种模式下，激发次数等于超声环阵换能器的阵元数目512。0057步骤S2，对各通道接收数据进行数据预处理，利用切除法将拾取的直达波信号区间赋值为零，具体流程如图2所示。0058在本实施例中，对环阵换能器各通道接收到的回波信号进行预处理，包括去除激发阵元传播到接收阵元位置的直达波信号。由于直达波信号的能量强，在回波数据曲线上表现为振幅大。首先借助matlab的plot函数绘制某一接收阵元的回波信号，再通过图像上的位置索引来确定直达波的区间宽度。对于每个接收阵元上的回波信号，通过最大值max函数找到振幅最强位置，再应用切除法将直达波区间范围信号置零。说明书5/10 页。

28、9CN 117503203 A900590060步骤S3，设置先验的声速范围来计算不同声速下的相干因子，得到每个成像点在每个预设声速下对应的相干因子值，具体流程见图3。0061在本实施例中，根据已知OAbreast二维切片成像介质，设置先验声速范围为14601620m/s，增量为1m/s。其中成像区域大小为，成像区域网格为。对每个先验声速，利用已知成像点位置、激发阵元位置以及接收阵元位置的几何距离计算延迟时间，再由延迟时间从各接收阵元中抽取信号，此时每个成像点处的声速均为恒定声速。延迟时间由下述公式(1)计算得到：00620063式中，表示序号激发阵元位置，表示序号接收阵元位置，表示坐标成像点。

29、位置，表示坐标成像点处的声速。0064在本实施例中，由于一个激发阵元对应阵元数目512个接收元，故能计算出512组延迟时间，从而能提取到512个延迟信号，再将提取的延迟信号代入下述公式(2)，计算该声速下的相干因子：0065说明书6/10 页10CN 117503203 A100066式中，表示坐标成像点位置的相干因子，N表示接收阵元的总数目，表 示具体 激发阵 元序号，表示具体 接收阵 元序号，表示坐标成像点位置在序号为m的激发阵元在序号为n接收阵元上提取的延迟信号，为序号m的激发阵元经过坐标成像点位置到达序号为n接收阵元上的时间。0067在本实施例中，对所有声速进行循环求解，在每个成像点位。

30、置均可以得到每个声速对应的相干因子值，因此最终得到矩阵。0068求解相干因子部分的伪代码如下所示：00690070步骤S4，将得到的相干因子值分别对水平和垂直方向进行平滑处理，再根据最大化相干因子来分别得到水平和垂直分布的声速；0071在 本 实 施 例 中，在 每 个 声 速 下，先 对 得 到矩 阵 抽 取的二维切片，对二维切片分别在水平方向和垂直方向上对相干因子求平均。为了进一步提高求解声速精度，使用interp1函数分别对相干因子和先验声速在一维上做插值。再根据最大化相干因子原则，每个成像点位置选择相干因子最大值时对应的声速，从而说明书7/10 页11CN 117503203 A11得。

31、到水平和垂直分布的声速。0072步骤S5，将水平和垂直分布的声速相加和平滑处理，得到最终的声速分布。0073在本实施例中，将得到水平分布和垂直分布声速进行相加求平均，此时得到的声速分布存在离散特征，需要将得到的声速在全局上做进一步的平滑处理。进一步平滑策略主要是进行9点平均加权滑动窗口平滑方法，见下述公式(3)00740075式中，表示坐标成像点位置的声速。0076步骤S6，根据成像点、激发阵元以及接收阵元的几何距离和得到的声速分布来计算延迟时间，通过延迟时间将各通道信号进行对齐；0077在本实施例中，进一步，在步骤S6中，利用步骤S5中得到的声速分布，在已知激发阵元、成像点和接收阵元的几何距。

32、离计算各阵元时间延迟，具体的是针对不同成像点，利用每个成像点的声速分布代替恒定声速代入公式(1)计算时间延迟，通过延迟时间将各通道信号进行对齐。0078步骤S7，使用延迟叠加(DAS)波束形成方法将对齐的各通道信号进行聚焦，从而得到高分辨率超声成像结果。0079在本实施例中，为了清楚的说明基于相位校正的超声图像分辨率的提升，与如图4(a)基于恒定声速的超声成像结果进行对比，可以看到利用恒定声速得到延时信号有一定的失准，尤其是乳腺内部两个高速散射点位置，可以看到两个位置的点不能聚焦，以及部分乳腺结构呈现的对比度差，导致超声成像结果分辨率低。图4(b)为基于相位校正的超声成像结果，可以看出结构轮廓。

33、更加清晰，尤其在两个高速散射点位置聚焦，使得成像分辨率得到提升。0080实施例20081本实施例中，一种用于超声环阵成像的相位畸变校正系统，用于实现实施例1的一种用于超声环阵成像的相位畸变校正方法，可以分成多个模块，每个模块具有不同的功能。如图5系统架构图所示，该系统包括数据采集模块、数据预处理模块、估计声速分布模块、相位校正模块、超声成像模块以及图像可视化模块。0082数据采集模块，用于设置超声环阵的阵元激发和接收模式，记录接收阵元接收到的回波信号。此处的数据采集可以是利用仿真数据，也可以是在超声CT设备上采集到的原始数据。设置每次仅有一个阵元激发，所有阵元接收信号，该过程将重复进行直至所有。

34、阵元都依次激发完毕。若利用超声设备进行数据采集时，在每个阵元激发之后，为防止多次反射波的干扰，在相邻阵元激发间隔一定的时间，确保环阵内的超声波衰减到可忽略的水平。0083数据预处理模块，用于对各通道接收数据进行预处理，主要是利用切除法将拾取的直达波信号区间赋值为零。根据直达波信号能量强、在回波数据曲线上振幅大，来确定直达波截取区间范围，然后将该区间范围信号置零。说明书8/10 页12CN 117503203 A120084估计声速分布模块，用于设置先验声速范围来计算不同声速下的相干因子，再根据最大化相干因子原则来得到每个成像点位置的声速分布。根据成像介质设定先验声速范围，然后计算不同声速和已知。

35、成像点位置、激发阵元位置以及接收阵元位置的几何距离下的延迟时间，从各通道信号中提取延迟信号。根据从每个接收阵元提取的延迟信号计算相干因子，再进行插值、最大化相干因子原则以及平滑处理得到最终声速分布。0085相位校正模块，用于已知成像点位置、激发阵元位置以及接收阵元位置的几何距离和每个成像点位置的声速分布计算校正后的延迟时间，通过计算得到的延迟时间将各通道信号进行对齐。像差波前由相位或时间延迟组成，这些延迟与由于组织非均质性引起的介质声速的局部变化有关，通过相位畸变校正可减少像差对成像结果的干扰。0086超声成像模块，用于将对齐的各通道信号使用延迟叠加(DAS)波束形成方法进行聚焦，从而得到高分。

36、辨率超声成像结果。0087图像可视化模块，用于显示相位校正后的高分辨率超声成像结果。在视觉上直接对比通过测量水温得到的恒定声速下的超声图像与基于相位校正的超声图像，体现一种用于超声环阵成像的相位畸变校正方法可以提高超声成像分辨率。0088实施例30089本实施例涉及一种计算机可读存储介质，其上存储有程序，该程序被处理器执行时，实现实施例1的一种用于超声环阵成像的相位畸变校正方法。0090实施例40091本实施例涉及一种用于超声环阵成像的相位畸变校正装置，包括存储器和一个或多个处理器，所述存储器中存储有可执行代码，所述一个或多个处理器执行所述可执行代码时，用于实现实施例1的一种用于超声环阵成像的。

37、相位畸变校正方法。0092还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。0093本领域技术人员应明白，本发明的实施例可提供为方法、系统或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可。

38、用程序代码的计算机可用存储介质（包括但不限于磁盘存储器、CDROM、光学存储器等）上实施的计算机程序产品的形式。0094本发明可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本发明，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。0095本发明中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其。

39、他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的说明书9/10 页13CN 117503203 A13部分说明即可。0096以上所述仅为本发明的实施例而已，并不用于限制本发明。对于本领域技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。说明书10/10 页14CN 117503203 A14图1说明书附图1/4 页15CN 117503203 A15图2说明书附图2/4 页16CN 117503203 A16图3图4（a）说明书附图3/4 页17CN 117503203 A17图4（b）图5说明书附图4/4 页18CN 117503203 A18。