CN117503203B - 一种用于超声环阵成像的相位畸变校正方法和系统 - Google Patents
一种用于超声环阵成像的相位畸变校正方法和系统 Download PDFInfo
- Publication number
- CN117503203B CN117503203B CN202410005094.5A CN202410005094A CN117503203B CN 117503203 B CN117503203 B CN 117503203B CN 202410005094 A CN202410005094 A CN 202410005094A CN 117503203 B CN117503203 B CN 117503203B
- Authority
- CN
- China
- Prior art keywords
- sound velocity
- array element
- imaging
- array
- imaging point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 122
- 238000000034 method Methods 0.000 title claims abstract description 77
- 238000012937 correction Methods 0.000 title claims abstract description 38
- 230000004075 alteration Effects 0.000 title claims abstract description 28
- 238000009826 distribution Methods 0.000 claims abstract description 59
- 230000005284 excitation Effects 0.000 claims abstract description 40
- 238000007781 pre-processing Methods 0.000 claims abstract description 17
- 238000009499 grossing Methods 0.000 claims abstract description 14
- 238000002604 ultrasonography Methods 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 13
- 238000012800 visualization Methods 0.000 claims description 6
- 125000004122 cyclic group Chemical group 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 230000002238 attenuated effect Effects 0.000 claims description 3
- 238000012935 Averaging Methods 0.000 claims description 2
- 238000012285 ultrasound imaging Methods 0.000 description 10
- 210000001519 tissue Anatomy 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 230000006870 function Effects 0.000 description 4
- 230000009897 systematic effect Effects 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 210000000481 breast Anatomy 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 230000001934 delay Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009213 extracorporeal shockwave therapy Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000005865 ionizing radiation Effects 0.000 description 1
- 210000005075 mammary gland Anatomy 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/13—Tomography
- A61B8/14—Echo-tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
- A61B8/4488—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer the transducer being a phased array
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4053—Scaling of whole images or parts thereof, e.g. expanding or contracting based on super-resolution, i.e. the output image resolution being higher than the sensor resolution
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Gynecology & Obstetrics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
一种用于超声环阵成像的相位畸变校正方法和系统,其方法包括:设置超声环阵的阵元激发和接收模式,采集超声回波信号;对各通道接收数据进行数据预处理,将拾取的直达波信号区间赋值为零;设置先验的声速范围来计算不同声速下的相干因子,得到每个成像点在每个先验声速下对应的相干因子值;将得到的相干因子值分别对水平和垂直方向进行插值,再根据最大化相干因子分别得到水平和垂直分布的声速;将水平和垂直分布的声速相加和平滑处理,得到最终的声速分布;根据成像点位置、激发阵元位置以及接收阵元位置的几何距离和得到的声速分布来计算延迟时间,通过延迟时间将各通道信号进行对齐;将对齐的各通道信号进行聚焦,得到高分辨率超声成像结果。
Description
技术领域
本发明属于超声成像技术,具体涉及一种用于超声环阵成像的相位畸变校正方法和系统。
背景技术
超声断层扫描成像技术是一种无电离辐射、高帧率的医学影像模态,它通过采集生物体软组织的透射和散射超声波形信号,将其反演为组织结构图像。超声成像以实时、安全、低成本等优点被广泛应用于临床。现有的临床超声成像探头是线阵或凸阵换能器,它们具有便携性以及更具成本效益的优势。但由于线阵或凸阵视野的受限,检测生物体组织的分辨率是各向异性的,横向分辨率低于轴向分辨率。其次,图像质量在很大程度上取决于生物体内部组织的深度,因为来自深度区域的回波强度取决于沿传播路径衰减的变化。此外,扫描图像的质量高度依赖于超声科医生的经验和培训水平。环阵换能器拥有360°的视野角度,超声波束可以聚焦在每个点上。此外,环阵换能器内的任何区域都有靠近该区域的检测元素,因此可以抑制图像质量的深度依赖性。环阵换能器不仅可以同时接收不同方向的回波,还可以同时接收散射波和透射波,利用超声环阵成像技术可以获得生物体组织的结构信息。
现有的超声环阵成像技术通过标准的延时叠加(delay-and-sum,DAS)波束形成方法得到超声图像,但由于延时叠加方法无法区分主瓣和副瓣的信号,且副瓣的信号含有干扰波和杂波,所成图像质量不佳尚未达到临床使用的标准。通过假设传播介质是各向同性的,根据恒定声速和已知激发阵元、成像点和接收阵元间的几何距离计算各阵元时间延迟,因此,恒定声速的选取误差影响最终的图像质量。在实际超声环阵成像时,通常将成像物体置于水中,假设水温保持不变且水温不会受到检测物体的影响,将水的声速近似为一个恒定声速,这与实际物理场景有差异而导致建模的系统性误差。一般情况下使用恒定声速可以得到较好的成像分辨率,但是当重建生物体不同组织间声速差异大、非均质性高时,文献(Pinton G F, Trahey G E, Dahl J J. Sources of image degradation infundamental and harmonic ultrasound imaging using nonlinear, full-wavesimulations[J]. IEEE transactions on ultrasonics, ferroelectrics, andfrequency control, 2011, 58(4): 754-765.)指出利用恒定声速进行波束形成方法时图像会出现局部结构变形、重影等模糊情况,造成超声图像分辨率低。文献(Kim H, Song I,Kang J, et al. Phase aberration correction for ultrasound imaging guidedextracorporeal shock wave therapy (ESWT): Feasibility study[J]. Ultrasonics,2023, 132: 107011.)中对比了不同恒定声速下的波束形成结果,通过视觉评价发现,若选取的恒定声速与真实声速的误差大时,重建图像的整体质量严重下降,成像目标的物理位置以及界面出现了明显扭曲。反之,当恒定声速误差小时,成像目标的空间分辨率得到提高,但随着成像深度的增加,重建图像的相位像差影响逐渐增大,图像质量和目标位置也随着深度的增加而逐渐失真。
现有技术的缺点表现为:
1)在实际情况下,通过水的声速近似得到恒定声速,这与实际物理场景有差异而导致建模的系统性误差,将导致通道信号的失准,造成超声成像的分辨率降低;
2)在遇到复杂介质或声速差异较大时,尤其当作用于存在显著的声速不均匀性介质时,则在利用恒定声速计算几何延迟后,高阶相位畸变会导致信号失准,随着成像深度的增加,重建图像的相位像差影响逐渐增大,超声图像出现局部结构变形、重影等分辨率降低情况。
发明内容
本发明要解决现有超声环阵成像技术的问题,提出一种用于超声环阵成像的相位畸变校正方法和系统。
本发明提出的一种用于超声环阵成像的相位畸变校正方法,包括以下步骤:
S1,设置超声环阵的阵元激发和接收模式,实现超声回波信号采集;
S2,对各通道接收数据进行数据预处理,利用切除法将拾取的直达波信号区间赋值为零;
S3,设置先验的声速范围来计算不同声速下的相干因子,得到每个成像点在每个先验声速下对应的相干因子值;
S4,将得到的相干因子值分别对水平和垂直方向进行插值,再根据最大化相干因子来分别得到水平和垂直分布的声速;
S5,将水平和垂直分布的声速相加和平滑处理,得到最终的声速分布;
S6,根据成像点、激发阵元以及接收阵元的几何距离和得到的声速分布来计算延迟时间,通过延迟时间将各通道信号进行对齐;
S7,使用延迟叠加(DAS)波束形成方法将对齐的各通道信号进行聚焦,从而得到高分辨率超声成像结果;
进一步,在步骤S1中,设置环阵换能器每次仅有一个阵元激发,所有阵元接收信号,该过程将重复进行直至所有阵元都依次激发完毕,在这种模式下,激发次数等于环阵换能器的阵元数目。
进一步,在步骤S1中,利用超声设备进行数据采集时,在每个阵元激发之后,在相邻阵元激发时需要有一定的时间间隔,以确保环阵内的超声波衰减到可忽略的水平。
进一步,在步骤S2中,对环阵换能器各通道接收到的回波信号进行预处理,主要去除激发阵元传播到接收阵元位置的直达波信号。数据预处理根据直达波信号能量强,在回波数据曲线上振幅大的特征,来确定直达波截取区间范围。通过找到接收回波信号上最大值位置,再将直达波区间范围信号归零。
进一步,在步骤S3中,根据成像介质设定先验的声速范围为1460-1620m/s,声速增量为1m/s。根据已知成像点位置、激发阵元位置以及接收阵元位置的几何距离和设定的先验声速计算延迟时间,从各通道信号中提取延迟信号,此时每个成像点处的声速均为恒定声速。延迟时间/>由公式(1)得到:
式中,表示/>序号激发阵元位置,/>表示/>序号接收阵元位置,/>表示坐标/>成像点位置,/>表示坐标/>成像点处的声速。
进一步,在步骤S3中,一个激发阵元对应多个接收阵元,故能计算出多组延迟时间,从而能提取到多个延迟信号。将从每个接收阵元提取的延迟信号代入下述公式(2),计算该声速下的相干因子:
式中,表示坐标/>成像点位置的相干因子, N表示接收阵元的总数目,表示具体激发阵元序号,/>表示具体接收阵元序号,表示坐标/>成像点位置在序号为m的激发阵元在序号为n接收阵元上提取的延迟信号,/>为序号m的激发阵元经过坐标/>成像点位置到达序号为n接收阵元上的时间。
进一步,在步骤S3中,对设置所有先验声速进行循环求解,在每个成像点位置均可以得到每个声速对应的相干因子值。
进一步,在步骤S4中,为了进一步提高求解精度,在相干因子和声速分布上做一维插值处理。考虑到直接求解全局声速会得到离散的声速分布,于是对每个声速下的二维相干因子矩阵分别在水平方向和垂直方向上求均值。
进一步,在步骤S4中,根据最大化相干因子原则,每个成像点位置选择相干因子最大值时对应的声速,从而得到水平和垂直分布的声速。
进一步,在步骤S5中,将得到水平和垂直分布的声速进行相加求平均,再将得到的声速做进一步的平滑处理。进一步平滑策略主要是进行9点平均加权滑动窗口平滑方法,见下述公式(3)
式中,表示坐标/>成像点位置的声速。
进一步,在步骤S6中,利用步骤S5中得到的声速分布,在每个成像点位置,利用声速分布代替恒定声速代入公式(1)计算延迟时间。
本发明的第二个方面涉及一种用于超声环阵成像的相位畸变校正系统,包括数据采集模块、数据预处理模块、估计声速分布模块、相位校正模块、超声成像模块以及图像可视化模块,其中,
数据采集模块,用于设置超声环阵的阵元激发和接收模式,利用接收阵元记录接收到的回波信号;
数据预处理模块,用于对各通道接收数据进行预处理,主要是利用切除法将拾取的直达波信号区间赋值为零;
估计声速分布模块,用于将先验声速范围和计算相干因子联系起来,从而得到每个成像点位置的声速分布。
相位校正模块,用于已知成像点位置、激发阵元位置以及接收阵元位置的几何距离和每个成像点位置的声速分布计算校正后的延迟时间,通过计算得到的延迟时间将各通道信号进行对齐;
超声成像模块,用于将对齐的各通道信号使用延迟叠加(DAS)波束形成方法进行聚焦,从而得到高分辨率超声成像结果;
图像可视化模块,用于显示相位校正后的高分辨率超声成像结果。
本发明的第三个方面涉及一种计算机可读存储介质,其上存储有程序,该程序被处理器执行时,实现本发明的一种用于超声环阵成像的相位畸变校正方法。
本发明的第四个方面涉及一种用于超声环阵成像的相位畸变校正装置,包括存储器和一个或多个处理器,所述存储器中存储有可执行代码,所述一个或多个处理器执行所述可执行代码时,用于实现本发明的一种用于超声环阵成像的相位畸变校正方法。
本发明的创新点是:
1)相比于现有的恒定声速估计方法,本发明利用回波数据计算在不同声速下延时信号的相干因子,使用插值、相干因子最大化原则以及平滑方法求解每个成像点位置的声速分布。
2)相比于现有超声环阵波束形成方法技术,本发明利用声速分布来代替恒定声速来计算延迟时间,可减少由实际物理建模中产生的系统性误差。
本发明基于所有合成孔径通道数据和相干因子来估计声速分布,用声速分布代替恒定声速计算各阵元时间延迟,通过标准的延时叠加波束形成方法,解决现有超声环阵成像技术在遇到生物体组织差异大或非均质性高时产生的相位畸变问题,从而提高超声图像的分辨率。
本发明的有益效果是:与人为测量恒定声速方法相比,本发明得到声速分布基于所有合成孔径通道数据,得到的声速分布可以减少与实际物理场景有差异而导致建模的系统性误差。进一步的,利用声速分布来代替恒定声速,再结合已知几何距离来计算延迟时间,可以解决在遇到复杂结构和声速不均匀介质时各通道信号不能很好的聚焦问题,可以有效避免信号失准,从而提高超声图像的分辨率。另外,本发明是估计每个成像点位置上的声速分布,与采样时间间隔无关,可适用于实际不同采样率的采集数据。
附图说明
图1是本发明方法的流程图;
图2是本发明的数据预处理模块的流程图;
图3是本发明估计声速分布模块的流程图;
图4(a)—图4(b)是本发明测试例中的一组结果展示:图4 (a)为基于恒定声速的超声成像结果;图4 (b)为基于相位校正的超声成像结果;
图5是本发明的系统架构图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
需要说明的是,在不冲突的情况下,下述的实施例及实施方式中的特征可以相互组合。
实施例1
本实施例提出一种用于超声环阵成像的相位畸变校正方法,实施例选用公开数据集OA-breast乳腺数据的二维切片,基于matlab平台和k-wave工具箱进行数值仿真。图1是该方法的整体流程图,具体包括以下步骤:
步骤S1,设置超声环阵的阵元激发和接收模式,实现超声回波信号采集。
在本实施例中,具体的仿真参数如下:超声环阵阵元数目为512,超声环阵半径为40mm,超声波激发的中心频率为3MHz,矩形计算域边长,计算网格点数为。设置超声环阵换能器每次一个阵元激发所有阵元接收信号,该过程将重复进行直至所有阵元都依次激发完毕。在这种模式下,激发次数等于超声环阵换能器的阵元数目512。
步骤S2,对各通道接收数据进行数据预处理,利用切除法将拾取的直达波信号区间赋值为零,具体流程如图2所示。
在本实施例中,对环阵换能器各通道接收到的回波信号进行预处理,包括去除激发阵元传播到接收阵元位置的直达波信号。由于直达波信号的能量强,在回波数据曲线上表现为振幅大。首先借助matlab的plot函数绘制某一接收阵元的回波信号,再通过图像上的位置索引来确定直达波的区间宽度。对于每个接收阵元上的回波信号,通过最大值max函数找到振幅最强位置,再应用切除法将直达波区间范围信号置零。
步骤S3,设置先验的声速范围来计算不同声速下的相干因子,得到每个成像点在每个预设声速下对应的相干因子值,具体流程见图3。
在本实施例中,根据已知OA-breast二维切片成像介质,设置先验声速范围为1460-1620m/s,增量为1m/s。其中成像区域大小为,成像区域网格为。对每个先验声速,利用已知成像点位置、激发阵元位置以及接收阵元位置的几何距离计算延迟时间/>,再由延迟时间从各接收阵元中抽取信号,此时每个成像点处的声速均为恒定声速。延迟时间/>由下述公式(1)计算得到:
式中,表示/>序号激发阵元位置,/>表示/>序号接收阵元位置,,/>表示坐标/>成像点位置,/>表示坐标/>成像点处的声速。
在本实施例中,由于一个激发阵元对应阵元数目512个接收元,故能计算出512组延迟时间,从而能提取到512个延迟信号,再将提取的延迟信号代入下述公式(2),计算该声速下的相干因子:
式中,表示坐标/>成像点位置的相干因子,N表示接收阵元的总数目,表示具体激发阵元序号,/>表示具体接收阵元序号,表示坐标/>成像点位置在序号为m的激发阵元在序号为n接收阵元上提取的延迟信号,/>为序号m的激发阵元经过坐标/>成像点位置到达序号为n接收阵元上的时间。
在本实施例中,对所有声速进行循环求解,在每个成像点位置均可以得到每个声速对应的相干因子值,因此最终得到矩阵。
求解相干因子部分的伪代码如下所示:
步骤S4,将得到的相干因子值分别对水平和垂直方向进行平滑处理,再根据最大化相干因子来分别得到水平和垂直分布的声速;
在本实施例中,在每个声速下,先对得到矩阵抽取的二维切片,对二维切片分别在水平方向和垂直方向上对相干因子求平均。为了进一步提高求解声速精度,使用interp1函数分别对相干因子和先验声速在一维上做插值。再根据最大化相干因子原则,每个成像点位置选择相干因子最大值时对应的声速,从而得到水平和垂直分布的声速。
步骤S5,将水平和垂直分布的声速相加和平滑处理,得到最终的声速分布。
在本实施例中,将得到水平分布和垂直分布声速进行相加求平均,此时得到的声速分布存在离散特征,需要将得到的声速在全局上做进一步的平滑处理。进一步平滑策略主要是进行9点平均加权滑动窗口平滑方法,见下述公式(3)
式中,表示坐标/>成像点位置的声速。
步骤S6,根据成像点、激发阵元以及接收阵元的几何距离和得到的声速分布来计算延迟时间,通过延迟时间将各通道信号进行对齐;
在本实施例中,进一步,在步骤S6中,利用步骤S5中得到的声速分布,在已知激发阵元、成像点和接收阵元的几何距离计算各阵元时间延迟,具体的是针对不同成像点,利用每个成像点的声速分布代替恒定声速代入公式(1)计算时间延迟,通过延迟时间将各通道信号进行对齐。
步骤S7,使用延迟叠加(DAS)波束形成方法将对齐的各通道信号进行聚焦,从而得到高分辨率超声成像结果。
在本实施例中,为了清楚的说明基于相位校正的超声图像分辨率的提升,与如图4(a)基于恒定声速的超声成像结果进行对比,可以看到利用恒定声速得到延时信号有一定的失准,尤其是乳腺内部两个高速散射点位置,可以看到两个位置的点不能聚焦,以及部分乳腺结构呈现的对比度差,导致超声成像结果分辨率低。图4(b)为基于相位校正的超声成像结果,可以看出结构轮廓更加清晰,尤其在两个高速散射点位置聚焦,使得成像分辨率得到提升。
实施例2
本实施例中,一种用于超声环阵成像的相位畸变校正系统,用于实现实施例1的一种用于超声环阵成像的相位畸变校正方法,可以分成多个模块,每个模块具有不同的功能。如图5系统架构图所示,该系统包括数据采集模块、数据预处理模块、估计声速分布模块、相位校正模块、超声成像模块以及图像可视化模块。
数据采集模块,用于设置超声环阵的阵元激发和接收模式,记录接收阵元接收到的回波信号。此处的数据采集可以是利用仿真数据,也可以是在超声CT设备上采集到的原始数据。设置每次仅有一个阵元激发,所有阵元接收信号,该过程将重复进行直至所有阵元都依次激发完毕。若利用超声设备进行数据采集时,在每个阵元激发之后,为防止多次反射波的干扰,在相邻阵元激发间隔一定的时间,确保环阵内的超声波衰减到可忽略的水平。
数据预处理模块,用于对各通道接收数据进行预处理,主要是利用切除法将拾取的直达波信号区间赋值为零。根据直达波信号能量强、在回波数据曲线上振幅大,来确定直达波截取区间范围,然后将该区间范围信号置零。
估计声速分布模块,用于设置先验声速范围来计算不同声速下的相干因子,再根据最大化相干因子原则来得到每个成像点位置的声速分布。根据成像介质设定先验声速范围,然后计算不同声速和已知成像点位置、激发阵元位置以及接收阵元位置的几何距离下的延迟时间,从各通道信号中提取延迟信号。根据从每个接收阵元提取的延迟信号计算相干因子,再进行插值、最大化相干因子原则以及平滑处理得到最终声速分布。
相位校正模块,用于已知成像点位置、激发阵元位置以及接收阵元位置的几何距离和每个成像点位置的声速分布计算校正后的延迟时间,通过计算得到的延迟时间将各通道信号进行对齐。像差波前由相位或时间延迟组成,这些延迟与由于组织非均质性引起的介质声速的局部变化有关,通过相位畸变校正可减少像差对成像结果的干扰。
超声成像模块,用于将对齐的各通道信号使用延迟叠加(DAS)波束形成方法进行聚焦,从而得到高分辨率超声成像结果。
图像可视化模块,用于显示相位校正后的高分辨率超声成像结果。在视觉上直接对比通过测量水温得到的恒定声速下的超声图像与基于相位校正的超声图像,体现一种用于超声环阵成像的相位畸变校正方法可以提高超声成像分辨率。
实施例3
本实施例涉及一种计算机可读存储介质,其上存储有程序,该程序被处理器执行时,实现实施例1的一种用于超声环阵成像的相位畸变校正方法。
实施例4
本实施例涉及一种用于超声环阵成像的相位畸变校正装置,包括存储器和一个或多个处理器,所述存储器中存储有可执行代码,所述一个或多个处理器执行所述可执行代码时,用于实现实施例1的一种用于超声环阵成像的相位畸变校正方法。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白,本发明的实施例可提供为方法、系统或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本发明,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本发明中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
以上所述仅为本发明的实施例而已,并不用于限制本发明。对于本领域技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (10)
1.一种用于超声环阵成像的相位畸变校正方法,其特征在于,包括如下步骤:
S1、设置超声环阵的阵元激发和接收模式,实现超声回波信号采集;
S2、对各通道接收数据进行数据预处理,利用切除法将拾取的直达波信号区间赋值为零;
S3、设置先验的声速范围来计算不同声速下的相干因子,得到每个成像点在每个先验声速下对应的相干因子值;
根据成像介质设定先验的声速范围为1460-1620m/s,声速增量为1m/s;根据已知成像点位置、激发阵元位置以及接收阵元位置的几何距离和设定的先验声速计算延迟时间,从各通道信号中提取延迟信号,此时每个成像点处的声速均为恒定声速;延迟时间/>由公式(1)得到:/>式中,/>表示/>序号激发阵元位置,/>表示/>序号接收阵元位置,/>表示坐标/>成像点位置,/>表示坐标/>成像点处的声速;
一个激发阵元对应多个接收阵元,故能计算出多组延迟时间,从而能提取到多个延迟信号;将从每个阵元提取的延迟信号代入下述公式(2),计算该声速下的相干因子:式中,/>表示坐标/>成像点位置的相干因子,/>代表接收阵元的数目,/>代表具体激发阵元序号,/>代表具体接收阵元序号,/>表示坐标/>成像点位置在序号为m的激发阵元在序号为n接收阵元上提取的延迟信号,/>为序号m的激发阵元经过坐标/>成像点位置到达序号为n接收阵元上的时间;
S4、将得到的相干因子值分别对水平和垂直方向进行插值,再根据最大化相干因子来分别得到水平和垂直分布的声速;
S5、将水平和垂直分布的声速相加和平滑处理,得到最终的声速分布;
S6、根据成像点、激发阵元以及接收阵元的几何距离和得到的声速分布来计算延迟时间,通过延迟时间将各通道信号进行对齐;
S7、使用延迟叠加(DAS)波束形成方法将对齐的各通道信号进行聚焦,从而得到高分辨率超声成像结果。
2.根据权利要求1所述的方法,其特征在于,在步骤S1中,设置环阵换能器每次仅有一个阵元激发,所有阵元接收信号,该过程将重复进行直至所有阵元都依次激发完毕,在这种模式下,激发次数等于环阵换能器的阵元数目。
3.根据权利要求2所述的方法,其特征在于,在步骤S1中,利用超声设备进行数据采集时,在每个阵元激发之后,在相邻阵元激发时需要有一定的时间间隔,以确保环阵内的超声波衰减到可忽略的水平。
4.根据权利要求1所述的方法,其特征在于,在步骤S2中,对环阵换能器各通道接收到的回波信号进行预处理,主要去除激发阵元直接传播到接收阵元位置的直达波信号;数据预处理根据直达波信号能量强,在回波数据曲线上振幅大的特征,来确定直达波截取区间范围;通过找到接收回波信号上最大值位置,再将直达波区间范围信号归零。
5.根据权利要求1所述的方法,其特征在于,步骤S4中,为了进一步提高求解声速精度,在相干因子和声速分布上做一维插值处理;考虑到直接求解全局声速会得到离散的声速分布,于是对每个声速下的二维相干因子矩阵分别在水平方向和垂直方向上求均值;再根据最大化相干因子原则,每个成像点位置选择相干因子最大值时对应的声速,从而得到水平和垂直分布的声速。
6.根据权利要求1所述的方法,其特征在于,步骤S5中,将得到水平和垂直分布的声速进行相加求平均,再将得到的声速做进一步的平滑处理;进一步平滑策略主要是进行9点平均加权滑动窗口平滑方法,见下述公式(3):式中,/>表示坐标/>成像点位置的声速。
7.根据权利要求1所述的方法,其特征在于,步骤S6中,利用步骤S5中得到的声速分布代替恒定声速计算延迟时间,通过延迟时间将各通道信号进行对齐。
8.一种用于超声环阵成像的相位畸变校正系统,其特征在于,包括:数据采集模块、数据预处理模块、估计声速分布模块、相位校正模块、超声成像模块以及图像可视化模块,其中,
数据采集模块,用于设置超声环阵的阵元激发和接收模式,利用接收阵元记录接收到的回波信号;
数据预处理模块,用于对各通道接收数据进行预处理,主要是利用切除法将拾取的直达波信号区间赋值为零;
估计声速分布模块,用于将先验声速范围和计算相干因子联系起来,从而得到每个成像点位置的声速分布;根据成像介质设定先验的声速范围为1460-1620m/s,声速增量为1m/s;根据已知成像点位置、激发阵元位置以及接收阵元位置的几何距离和设定的先验声速计算延迟时间,从各通道信号中提取延迟信号,此时每个成像点处的声速均为恒定声速;延迟时间/>由公式(1)得到:/>式中,/>表示/>序号激发阵元位置,/>表示/>序号接收阵元位置,/>表示坐标/>成像点位置,/>表示坐标/>成像点处的声速;
一个激发阵元对应多个接收阵元,故能计算出多组延迟时间,从而能提取到多个延迟信号;将从每个阵元提取的延迟信号代入下述公式(2),计算该声速下的相干因子:式中,/>表示坐标/>成像点位置的相干因子, />代表接收阵元的数目,/>代表具体激发阵元序号,/>代表具体接收阵元序号,/>表示坐标/>成像点位置在序号为m的激发阵元在序号为n接收阵元上提取的延迟信号,/>为序号m的激发阵元经过坐标/>成像点位置到达序号为n接收阵元上的时间;
相位校正模块,用于已知成像点位置、激发阵元位置以及接收阵元位置的几何距离和每个成像点位置的声速分布计算校正后的延迟时间,通过计算得到的延迟时间将各通道信号进行对齐;
超声成像模块,用于将对齐的各通道信号使用延迟叠加(DAS)波束形成方法进行聚焦,从而得到高分辨率超声成像结果;
图像可视化模块,用于显示相位校正后的高分辨率超声成像结果。
9.一种计算机可读存储介质,其特征在于,其上存储有程序,该程序被处理器执行时,实现权利要求1-7中任一项所述的一种用于超声环阵成像的相位畸变校正方法。
10.一种用于超声环阵成像的相位畸变校正装置,其特征在于,包括存储器和一个或多个处理器,所述存储器中存储有可执行代码,所述一个或多个处理器执行所述可执行代码时,用于实现权利要求1-7中任一项所述的一种用于超声环阵成像的相位畸变校正方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410005094.5A CN117503203B (zh) | 2024-01-03 | 2024-01-03 | 一种用于超声环阵成像的相位畸变校正方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410005094.5A CN117503203B (zh) | 2024-01-03 | 2024-01-03 | 一种用于超声环阵成像的相位畸变校正方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117503203A CN117503203A (zh) | 2024-02-06 |
CN117503203B true CN117503203B (zh) | 2024-03-22 |
Family
ID=89749703
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410005094.5A Active CN117503203B (zh) | 2024-01-03 | 2024-01-03 | 一种用于超声环阵成像的相位畸变校正方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117503203B (zh) |
Citations (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5331964A (en) * | 1993-05-14 | 1994-07-26 | Duke University | Ultrasonic phased array imaging system with high speed adaptive processing using selected elements |
WO1996023387A1 (en) * | 1995-01-23 | 1996-08-01 | Commonwealth Scientific And Industrial Research Organisation | Phase and/or amplitude aberration correction for imaging |
US5566675A (en) * | 1995-06-30 | 1996-10-22 | Siemens Medical Systems, Inc. | Beamformer for phase aberration correction |
US5793701A (en) * | 1995-04-07 | 1998-08-11 | Acuson Corporation | Method and apparatus for coherent image formation |
US5865753A (en) * | 1997-05-27 | 1999-02-02 | Friauf; Walter S. | Coherent detection ultrasound system |
US5910115A (en) * | 1997-09-22 | 1999-06-08 | General Electric Company | Method and apparatus for coherence filtering of ultrasound images |
US6508764B1 (en) * | 2000-10-31 | 2003-01-21 | Koninklijke Philips Electronics N.V. | Aberration correction apparatus and methods |
CN101077306A (zh) * | 2007-07-05 | 2007-11-28 | 中国计量科学研究院 | 一种相控阵超声二维阵列三维成像方法及其应用 |
WO2012008217A1 (ja) * | 2010-07-14 | 2012-01-19 | 株式会社日立メディコ | 超音波画像の画像復元方法およびその装置並びに超音波診断装置 |
WO2013082455A1 (en) * | 2011-12-01 | 2013-06-06 | Maui Imaging, Inc. | Motion detection using ping-based and multiple aperture doppler ultrasound |
JP5504357B1 (ja) * | 2013-01-09 | 2014-05-28 | 日立アロカメディカル株式会社 | 超音波診断装置 |
WO2016073976A1 (en) * | 2014-11-07 | 2016-05-12 | Tessonics Corporation | An ultrasonic adaptive beamforming method and its application for transcranial imaging |
JP2017104476A (ja) * | 2015-03-09 | 2017-06-15 | 炭 親良 | ビームフォーミング方法、計測イメージング装置、及び、通信装置 |
CN106901693A (zh) * | 2017-01-19 | 2017-06-30 | 南京大学 | 基于环形超声换能器阵列的大景深光声‑声速双模成像法 |
JP2019030623A (ja) * | 2016-11-15 | 2019-02-28 | 炭 親良 | ビームフォーミング方法、計測イメージング装置、及び、通信装置 |
CN109754407A (zh) * | 2019-01-10 | 2019-05-14 | 青岛海信医疗设备股份有限公司 | 超声图像处理方法、装置及设备 |
WO2019219549A1 (en) * | 2018-05-15 | 2019-11-21 | Koninklijke Philips N.V. | Synthetic transmit focusing ultrasound system with speed of sound aberration correction |
JP2020039841A (ja) * | 2017-11-24 | 2020-03-19 | 炭 親良 | 計測イメージング装置 |
WO2020081962A1 (en) * | 2018-10-19 | 2020-04-23 | Duke University | Methods, systems and computer program products for ultrasound imaging using coherence contribution |
AU2020104259A4 (en) * | 2020-12-07 | 2021-03-11 | Shanxi University | Ultrafast dynamics imaging method and system based on single-molecule quantum coherence |
CN114755311A (zh) * | 2022-06-14 | 2022-07-15 | 之江实验室 | 一种超声探头阵元一致性检测的方法及测量装置 |
CN114964581A (zh) * | 2022-06-01 | 2022-08-30 | 合肥工业大学智能制造技术研究院 | 一种基于超声相控阵聚焦原理的应力检测方法 |
WO2023091751A1 (en) * | 2021-11-22 | 2023-05-25 | The Trustees Of Columbia University In The City Of New York | Systems and methods for cavitation mapping with spatial-temporal parallel processing |
CN116482231A (zh) * | 2023-04-12 | 2023-07-25 | 中国工程物理研究院化工材料研究所 | 一种材料内部缺陷高信噪比超声成像方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8088067B2 (en) * | 2002-12-23 | 2012-01-03 | Insightec Ltd. | Tissue aberration corrections in ultrasound therapy |
US7207942B2 (en) * | 2003-07-25 | 2007-04-24 | Siemens Medical Solutions Usa, Inc. | Adaptive grating lobe suppression in ultrasound imaging |
US7744532B2 (en) * | 2004-03-31 | 2010-06-29 | Siemens Medical Solutions Usa, Inc. | Coherence factor adaptive ultrasound imaging methods and systems |
US8254654B2 (en) * | 2007-10-31 | 2012-08-28 | University Of Southern California | Sidelobe suppression in ultrasound imaging using dual apodization with cross-correlation |
CA2789129C (en) * | 2010-02-08 | 2017-08-22 | Dalhousie University | Ultrasound imaging system using beamforming techniques for phase coherence grating lobe suppression |
KR102103137B1 (ko) * | 2012-03-26 | 2020-04-22 | 마우이 이미징, 인코포레이티드 | 가중 인자들을 적용함으로써 초음파 이미지 품질을 향상시키는 시스템들 및 방법들 |
US11397167B2 (en) * | 2016-11-22 | 2022-07-26 | The Board Of Trustees Of The Leland Stanford Junior University | Local speed of sound estimation method for medical ultrasound |
EP4034849A1 (en) * | 2019-09-25 | 2022-08-03 | Institut National De La Sante Et De La Recherche Medicale - Inserm | A method for measuring the speed of sound in liver with a moving probe and associated methods and devices |
JP2023114623A (ja) * | 2022-02-07 | 2023-08-18 | キヤノンメディカルシステムズ株式会社 | 超音波診断装置 |
-
2024
- 2024-01-03 CN CN202410005094.5A patent/CN117503203B/zh active Active
Patent Citations (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5331964A (en) * | 1993-05-14 | 1994-07-26 | Duke University | Ultrasonic phased array imaging system with high speed adaptive processing using selected elements |
WO1996023387A1 (en) * | 1995-01-23 | 1996-08-01 | Commonwealth Scientific And Industrial Research Organisation | Phase and/or amplitude aberration correction for imaging |
US5793701A (en) * | 1995-04-07 | 1998-08-11 | Acuson Corporation | Method and apparatus for coherent image formation |
US5566675A (en) * | 1995-06-30 | 1996-10-22 | Siemens Medical Systems, Inc. | Beamformer for phase aberration correction |
US5865753A (en) * | 1997-05-27 | 1999-02-02 | Friauf; Walter S. | Coherent detection ultrasound system |
US5910115A (en) * | 1997-09-22 | 1999-06-08 | General Electric Company | Method and apparatus for coherence filtering of ultrasound images |
US6508764B1 (en) * | 2000-10-31 | 2003-01-21 | Koninklijke Philips Electronics N.V. | Aberration correction apparatus and methods |
CN101077306A (zh) * | 2007-07-05 | 2007-11-28 | 中国计量科学研究院 | 一种相控阵超声二维阵列三维成像方法及其应用 |
WO2012008217A1 (ja) * | 2010-07-14 | 2012-01-19 | 株式会社日立メディコ | 超音波画像の画像復元方法およびその装置並びに超音波診断装置 |
WO2013082455A1 (en) * | 2011-12-01 | 2013-06-06 | Maui Imaging, Inc. | Motion detection using ping-based and multiple aperture doppler ultrasound |
JP5504357B1 (ja) * | 2013-01-09 | 2014-05-28 | 日立アロカメディカル株式会社 | 超音波診断装置 |
WO2016073976A1 (en) * | 2014-11-07 | 2016-05-12 | Tessonics Corporation | An ultrasonic adaptive beamforming method and its application for transcranial imaging |
JP2017104476A (ja) * | 2015-03-09 | 2017-06-15 | 炭 親良 | ビームフォーミング方法、計測イメージング装置、及び、通信装置 |
JP2019030623A (ja) * | 2016-11-15 | 2019-02-28 | 炭 親良 | ビームフォーミング方法、計測イメージング装置、及び、通信装置 |
CN106901693A (zh) * | 2017-01-19 | 2017-06-30 | 南京大学 | 基于环形超声换能器阵列的大景深光声‑声速双模成像法 |
JP2020039841A (ja) * | 2017-11-24 | 2020-03-19 | 炭 親良 | 計測イメージング装置 |
WO2019219549A1 (en) * | 2018-05-15 | 2019-11-21 | Koninklijke Philips N.V. | Synthetic transmit focusing ultrasound system with speed of sound aberration correction |
WO2020081962A1 (en) * | 2018-10-19 | 2020-04-23 | Duke University | Methods, systems and computer program products for ultrasound imaging using coherence contribution |
CN109754407A (zh) * | 2019-01-10 | 2019-05-14 | 青岛海信医疗设备股份有限公司 | 超声图像处理方法、装置及设备 |
AU2020104259A4 (en) * | 2020-12-07 | 2021-03-11 | Shanxi University | Ultrafast dynamics imaging method and system based on single-molecule quantum coherence |
WO2023091751A1 (en) * | 2021-11-22 | 2023-05-25 | The Trustees Of Columbia University In The City Of New York | Systems and methods for cavitation mapping with spatial-temporal parallel processing |
CN114964581A (zh) * | 2022-06-01 | 2022-08-30 | 合肥工业大学智能制造技术研究院 | 一种基于超声相控阵聚焦原理的应力检测方法 |
CN114755311A (zh) * | 2022-06-14 | 2022-07-15 | 之江实验室 | 一种超声探头阵元一致性检测的方法及测量装置 |
CN116482231A (zh) * | 2023-04-12 | 2023-07-25 | 中国工程物理研究院化工材料研究所 | 一种材料内部缺陷高信噪比超声成像方法 |
Non-Patent Citations (6)
Title |
---|
Joint generalized coherence factor and minimum variance beamformer for synthetic aperture ultrasound imaging;Zhengfeng Lan;《IEEE transactions on ultrasonics, ferroelectrics, and frequency control》;20211231;第1167-1183页 * |
冒秋琴 ; 陈尧 ; 石文泽 ; 卢超 ; 李秋锋 ; .频域相位相干合成孔径聚焦超声成像研究.仪器仪表学报.(02),第137-147页. * |
太赫兹稀疏阵列相干因子栅旁瓣抑制方法;宋叶江;崔振茂;安健飞;喻洋;成彬彬;;红外与毫米波学报;20200815(04);第40-46页 * |
最小方差波束形成与广义相干系数融合的医学超声成像方法;吴文焘;蒲杰;吕D;;声学学报;20110115(01);第68-74页 * |
结合子空间投影的幅度相位估计波束形成算法;佘黎煌;闫鑫;张石;;东北大学学报(自然科学版);20180715(07);第18-21页 * |
非均匀组织医学超声发射声场仿真;解卓丽;周浩;郑音飞;;声学学报;20131115(06);第21-26页 * |
Also Published As
Publication number | Publication date |
---|---|
CN117503203A (zh) | 2024-02-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11446006B2 (en) | Adjusting measurements of the effects of acoustic radiation force for background motion effects | |
Deng et al. | Ultrasonic shear wave elasticity imaging sequencing and data processing using a verasonics research scanner | |
KR101868381B1 (ko) | 의료용 초음파 이미징에서의 전단파 정보의 해석 | |
KR102332464B1 (ko) | 분석 데이터로부터 전단파 추정 | |
JP5771758B2 (ja) | 超音波診断装置 | |
Gammelmark et al. | 2-D tissue motion compensation of synthetic transmit aperture images | |
US20140005548A1 (en) | Ultrasonic acoustic radiation force excitation for ultrasonic material property measurement and imaging | |
KR20170019327A (ko) | 음향 방사력 이미징에서의 적응식 모션 추정 | |
WO2015128656A1 (en) | Methods and systems for measuring properties with ultrasound | |
CN112754529B (zh) | 基于频域迁移的超声平面波成像方法、系统及存储介质 | |
Gauss et al. | Wavefront estimation in the human breast | |
KR20150112057A (ko) | 공간 일관성 기초 초음파 신호 처리 모듈 및 그에 의한 초음파 신호 처리 방법 | |
CN110840488A (zh) | 一种基于剪切波的成像方法、系统及装置 | |
Noda et al. | Ultrasound imaging with a flexible probe based on element array geometry estimation using deep neural network | |
JP7253540B2 (ja) | 患者適応型のせん断波生成による超音波せん断波イメージング | |
US20220280138A1 (en) | Ultrasonic imaging of acoustic attenuation coefficients with elevation compounding | |
Lou et al. | Spatial smoothing coherence factor for ultrasound computed tomography | |
CN117503203B (zh) | 一种用于超声环阵成像的相位畸变校正方法和系统 | |
CN111388012B (zh) | 用于检测组织硬度的方法、设备及系统 | |
CN113424073A (zh) | 材料非线性体积弹性的超声估算 | |
CN115869010A (zh) | 用于波束形成声速选择的系统和方法 | |
EP2853918B1 (en) | Shear wave detection in medical ultrasound imaging | |
Wang et al. | Multiline acquisition beamforming for ultrasound computed tomography | |
KR101610877B1 (ko) | 공간 일관성 기초 초음파 신호 처리 모듈 및 그에 의한 초음파 신호 처리 방법 | |
Feigin et al. | Computing Speed-of-Sound From Ultrasound: User-Agnostic Recovery and a New Benchmark |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |