CN110840484B - 自适应匹配最优声速的超声成像方法、装置及超声设备 - Google Patents
自适应匹配最优声速的超声成像方法、装置及超声设备 Download PDFInfo
- Publication number
- CN110840484B CN110840484B CN201911184446.3A CN201911184446A CN110840484B CN 110840484 B CN110840484 B CN 110840484B CN 201911184446 A CN201911184446 A CN 201911184446A CN 110840484 B CN110840484 B CN 110840484B
- Authority
- CN
- China
- Prior art keywords
- sound velocity
- roi
- depth
- optimal
- phase coherence
- 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 86
- 238000005070 sampling Methods 0.000 claims abstract description 122
- 238000000034 method Methods 0.000 claims abstract description 68
- 238000004364 calculation method Methods 0.000 claims abstract description 34
- 230000008569 process Effects 0.000 claims abstract description 32
- 238000002604 ultrasonography Methods 0.000 claims description 22
- 230000015572 biosynthetic process Effects 0.000 claims description 12
- 238000003786 synthesis reaction Methods 0.000 claims description 12
- 238000012285 ultrasound imaging Methods 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 210000004872 soft tissue Anatomy 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 101100421134 Schizosaccharomyces pombe (strain 972 / ATCC 24843) sle1 gene Proteins 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 239000000523 sample Substances 0.000 description 7
- 210000001519 tissue Anatomy 0.000 description 6
- 238000013441 quality evaluation Methods 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 5
- 230000009471 action Effects 0.000 description 3
- 230000003111 delayed effect Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000001303 quality assessment method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/58—Testing, adjusting or calibrating the diagnostic device
- A61B8/585—Automatic set-up of the device
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种自适应匹配最优声速的超声成像方法、装置及超声设备,通过获取待成像图像的通道数据,并根据预设边界信息从成像区域中划分出ROI区域,然后在每个声速下根据对应声速计算出ROI区域中每个采样点到各个阵元的延时时间,从而得到与每个采样点对应的多个延时时间,并根据各个延时时间从通道数据中得到相应的接收延时信号,并根据每个声速下与每个采样点对应的多个接收延时信号得到各声速下与每个采样点对应的相位相干因子,根据各个相位相干因子得到ROI区域的最优声速;本申请不需要进行横向傅里叶变换,且相位相干因子的计算过程简单、计算量小,因此本发明在使用过程中计算简单高效,有利于提高自适应匹配最优声速的效率。
Description
技术领域
本发明实施例涉及超声成像技术领域,特别是涉及一种自适应匹配最优声速的超声成像方法、装置及超声设备。
背景技术
超声成像技术在医学中的应用十分广泛,超声图像为医生对病情的诊断和治疗提供重要参考依据。在医疗超声成像过程中,声速是重要的成像参数之一,由于超声波在人体不同组织中传播时,声速会差异较大,若用于波束合成的声速为偏离真实声速的超声波在经过局部组织时,就会出现相位偏差,导致组织的偏移,影响图像的成像质量,从而影响医疗诊断的准确性。因此,在超声成像过程中匹配出最优声速是提高超声图像成像质量的必要手段。
目前,如专利US 20150073276A1—《aberration correction using channeldata in ultrasound imaging system》,中描述了一种超声通道数据的处理及分析来评估当前成像质量最优声速的方法,该专利中主要是将对应一帧B图像的通道ADC数据,经过一系列不同迭代声速的波束合成处理后,得到各个迭代声速下的波束合成后的RF数据;根据系统预设的ROI深度方向的上下边界线以及横向的左右边界线,截取出ROI框范围内的RF数据,对这部分数据进行横向傅里叶变换,得到各个声速下不同深度采样点的横向空间频谱;设定低频阈值频率,以归一化处理后的零到低频阈值频率的积分和作为该深度采样点的图像质量评估因子,该质量评估因子越大表明图像的横向分辨率越好。用ROI框内各个深度采样点的频谱加权求和后的频谱,作为该ROI区域的平均频谱,该频谱的归一化后的低频部分积分和就是ROI区域图像质量评估因子。绘制各个迭代声速的ROI区域图像质量评估因子的曲线,对该曲线拟合平滑处理之,就可以求得该曲线的极大值,以该极大值所对应的迭代声速作为当前成像的最优声速。但是,专利US 20150073276A1中提出的方法需要在ROI区域各深度采样点处做横向傅里叶变换,对于实时系统成像来说,运算量较大,影响整体工作效率。
鉴于此,如何提供一种解决上述技术问题的自适应匹配最优声速的超声成像方法、装置及超声设备成为本领域技术人员目前需要解决的问题。
发明内容
本发明实施例的目的是提供一种自适应匹配最优声速的超声成像方法、装置及超声设备,在使用过程中计算简单高效,有利于提高自适应匹配最优声速的效率。
为解决上述技术问题,本发明实施例提供了一种自适应匹配最优声速的超声成像方法,包括:
获取待成像图像的通道数据;
依据预设边界信息从成像区域中划分出ROI区域;
在每个声速下,根据相应的声速得到所述ROI区域内每个采样点到各个阵元的延时时间,并依据每个所述延时时间从所述通道数据中得到每个所述采样点对应各个所述阵元的接收延时信号;
根据各个所述接收延时信号得到每个所述声速下各个所述采样点的相位相干因子;
根据各个所述相位相干因子,得到所述ROI区域的最优声速。
可选的,各个所述采样点为每个所述声速下每条接收线上的不同深度的采样点;
所述根据各个所述相位相干因子,得到所述ROI区域的最优声速的过程包括:
对各条所述接收线上位于同一深度的采样点的相位相干因子求加权平均,得到每个所述声速下不同深度处的相位相干因子;
将每个所述声速下不同深度处的相位相干因子求加权平均,得到每个所述声速的相位相干因子;
将每个所述声速的相位相干因子进行拟合得到拟合曲线,并将所述拟合曲线的极大值对应的声速作为所述ROI区域的最优声速。
可选的,所述根据各个所述接收延时信号得到每个所述声速下各个所述采样点的相位相干因子的过程包括:
依据与每个所述声速下、各条所述接收线上每个所述采样点对应的各个接收延时信号及相位相干因子计算关系式,得到每个所述采样点的相位相干因子,其中,所述相位相干因子计算关系式为:
可选的,所述ROI区域为一个;所述方法还包括:
将所述ROI区域的最优声速作为所述成像区域的全局最优声速;
依据所述成像区域的全局最优声速进行波束合成,得到与所述成像区域对应的全局最优超声图像。
可选的,所述ROI区域为多个,各个所述ROI区域沿所述成像区域的深度方向排列,各个所述ROI区域依次相连构成所述成像区域;所述方法还包括:
分别根据每个所述ROI区域各自对应的最优声速进行波束合成,得到与每个所述ROI区域各自对应的子超声图像;
将各个所述子超声图像进行拼接,得到与所述成像区域对应的全局最优超声图像。
可选的,所述分别根据每个所述ROI区域各自对应的最优声速进行波束合成,得到与每个所述ROI区域各自对应的子超声图像的过程为:
依据与第k个ROI区域对应的最优声速,按照预设规则对所述第k个ROI区域在深度方向的起始位置、结束位置和相邻采样点的间距进行调节;其中,k=1,2,3,…K,K为ROI区域总数量;
依据与所述第k个ROI区域对应的最优声速以及调节后的起始位置、结束位置和相邻采样点的间距进行波束合成,得到与所述第k个ROI区域对应的子超声图像;其中,所述预设规则为:
start_segk=depth_startk*speedoptk/1540;
end_segk=depth_endk*speedoptk/1540;
samplestep_segk=samplestep*speedoptk/1540;
其中,1540为理论上超声在软组织内传播的平均速度,depth_startk和depth_endk分别表示第k个ROI区域的起始位置和结束位置;start_segk和end_segk分别表示第k个ROI区域调节后的起始位置和结束位置;samplestep表示声速1540下相邻采样点的间距;speedoptk表示第k个ROI区域的最优声速;samplestep_segk表示与第k个ROI区域对应的调节后的相邻采样点的间距。
可选的,所述ROI区域为整个所述成像区域,所述方法还包括:
根据每个所述声速下不同深度处的相位相干因子,得到与每个所述深度处的最大相位相干因子对应的目标声速,以得到与每个所述深度各自对应的最优声速;
依据与每个所述深度各自对应的最优声速建立深度与最优声速的关系曲线;
将所述深度与最优声速的关系曲线沿深度方向划分为M段,并计算出每一段的平均最优声速,所述M为多个;
依据各段的平均最优声速,得到每一段的局部声速;
依据与每段各自对应的局部声速,建立深度与局部声速的关系曲线。
可选的,所述依据各段的平均最优声速,得到每一段的局部声速的过程包括:
依据每段的平均最优声速以及局部声速计算关系式,得到各段的局部声速;其中,所述局部声速计算关系式为:
本发明实施例相应的提供了一种自适应匹配最优声速的超声成像装置,包括:
获取模块,用于获取待成像图像的通道数据;
划分模块,用于依据预设边界信息从成像区域中划分出ROI区域;
延时模块,用于在每个声速下,根据相应的声速得到所述ROI区域内每个采样点到各个阵元的延时时间,并依据每个所述延时时间从所述通道数据中得到每个所述采样点对应各个所述阵元的接收延时信号;
计算模块,用于根据各个所述接收延时信号得到每个所述声速下各个所述采样点的相位相干因子;
分析模块,用于根据各个所述相位相干因子,得到所述ROI区域的最优声速。
本发明实施例还提供了一种超声设备,包括:
存储器,用于存储计算机程序;
处理器,用于执行所述计算机程序时实现如上述所述自适应匹配最优声速的超声成像方法的步骤。
本发明实施例提供了一种自适应匹配最优声速的超声成像方法、装置及超声设备,本申请中在进行超声成像自适应匹配最优声速时,通过获取待成像图像的通道数据,并根据预设边界信息从成像区域中划分出ROI区域,然后在每个声速下,根据对应声速计算出ROI区域中每个采样点到各个阵元的延时时间,从而得到与每个采样点对应的多个延时时间,并根据各个延时时间从通道数据中得到相应的接收延时信号,并根据每个声速下与每个采样点对应的多个接收延时信号得到各声速下与每个采样点对应的相位相干因子,并进一步根据各个相位相干因子得到ROI区域的最优声速;该过程中不需要进行横向傅里叶变换,且相位相干因子的计算过程简单、计算量小,因此本发明在使用过程中计算简单高效,有利于提高自适应匹配最优声速的效率。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对现有技术和实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种自适应匹配最优声速的超声成像方法的流程示意图;
图2为本发明实施例提供的一种线阵下ROI区域划分方式示意图;
图3为本发明实施例提供的一种线阵的接收延时计算示意图;
图4为本发明实施例提供的一种根据相位相干因子计算最优声速的流程示意图;
图5为与图4对应的最优声速计算示意图;
图6为本发明实施例提供的一种硬件实时计算SCF的LUT表示意图;
图7为本发明实施例提供的一种线阵下多ROI区域的划分示意图;
图8为本发明实施例提供的一种最优声速与局部声速的示意图;
图9为与图8中对应的线阵探头的成像区域中各部分组织的声速示意图;
图10为本发明实施例提供的一种自适应匹配最优声速的超声成像装置的结构框图。
具体实施方式
本发明实施例提供了一种自适应匹配最优声速的超声成像方法、装置及超声设备,在使用过程中计算简单高效,有利于提高自适应匹配最优声速的效率。
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参照图1,图1为本发明实施例提供的一种自适应匹配最优声速的超声成像方法的流程示意图。该方法包括:
S11:获取待成像图像的通道数据;
具体的,例如待成像图像为B图像,则获取一帧B图像的通道数据,其中,一般通过弱聚焦发射波束方式,16次发射的数据就可以通过多波束合成的方式组成一帧B图像。另外,具体如何获取一帧B图像的通道数据为现有技术,本申请对此不做详述。
S12:依据预设边界信息从成像区域中划分出ROI区域;
需要说明的是,成像区域确定了待成像图像成像的位置,在进行ROI区域划分之前,可以预先设置好用来圈定ROI区域的边界信息,其中,预设边界信息包括横向起始百分比参数fraction_start、横向结束百分比参数fraction_end、深度方向起始边界线depth_start和深度方向结束边界线depth_end,以线阵探头为例,本实施例中提供的一种ROI区域划分方式可参照图2,其中,图2中上端的一排为线阵探头的各个阵元。
另外,在不同声速下进行划分ROI区域时,由于不同声速下相同的采样对象之间的间距可能会发生拉伸或者压缩,为了使在每个声速下所划分出的ROI区域中所包含的采样对象一一对应相同,因此,可以分别根据相应的声速对所划分的ROI区域在深度方向的起始边界线和结束边界线的位置进行调节,例如对于声速cl来说,划分出的ROI区域在深度方向的起始边界线位置为depth_startl=(cl/c0)*depth_start,结束边界线位置为depth_endl=(cl/c0)*depth_end,其中,c0为系统声速,其数值为1540m/s。
S13:在每个声速下,根据相应的声速得到ROI区域内每个采样点到各个阵元的延时时间,并依据每个延时时间从通道数据中得到每个采样点对应各个阵元的接收延时信号;
需要说明的是,可以预先设置好迭代声速起始值和结束值,以及迭代声速步距,例如迭代声速的起始值和结束值分别为1440m/s和1580m/s,也即迭代声速范围为1440~1580m/s,迭代声速步距为10m/s,因此,在对迭代声速过程中会得到各个声速,并且在每个声速下均计算出ROI区域内每个采样点到各个阵元的延时时间,然后根据相应的延时时间从通道数据中获取对应的接收延时信号,然后再根据每个采样点对应的各个接收延时信号计算出相应采样点的相位相干因子。也即,在声速迭代过程中,每个声速均计算出在该声速下各个采样点的相位相干因子后,在进行下一个声速下各个采样点的相位相干因子计算,直至最后一个声速下各个采样点的相位相干因子计算完成。
可以理解的是,在进行延时时间计算时,根据ROI区域内的采样点到接收孔径中的某个阵元之间的声程除以对应的声速,即可得到采样点到该阵元的延时时间,然后根据该延时时间从通道数据的信号数据中获取与该延时时间对应的接收延时信号,也即,获取出在相应声速下对应采样点到阵元的接收延时信号。
例如,共有3个声速分别为c1、c2、c3,ROI区域内的采样点有N个,接收孔径中阵元数量共为128个,则在声速c1下计算出N个采样点各自到各个阵元的延时时间,则每个采样点对应128个延时时间,相应的根据每个延时时间能够从通道数据中得到对应的接收延时信号,则每个采样点对应128个接收延时信号;同样,在声速c2下每个采样点对应得到与128个阵元各自对应的延时时间,每个采样点对应128个接收延时信号,在声速c3下,同样一个采样点对应128个接收延时信号。
还需要说明的是,为了进一步保证在不同声速下,所划分出的ROI区域中包含的采样点一一对应相同,在进行延时时间计算时,在每个声速下可以根据当前声速对深度方向的相邻采样点的间距进行调节,例如,在声速cl下深度方向的相邻采样点的间距可以调节为(cl/2)/Fs,其中,Fs为系统采样率,在系统默认声速1540m/s下,深度方向采样间距为(c0/2)/Fs。
S14:根据各个接收延时信号得到每个声速下各个采样点的相位相干因子;
S15:根据各个相位相干因子,得到ROI区域的最优声速。
由于图像聚焦越好,相位相干因子的值越接近于1,反之越接近于0,因此,本申请中可以将相位相干因子用于图像质量评估。具体的,在上述得到每个声速下各个采样点的各个接收延时信号后,根据相应声速下每个采样点的各个接收延时信号,能够计算出在该声速下每个采样点的相位相干因子,然后根据每个声速下每个采样点的相位相干因子就能够得出ROI区域的最优声速。由于根据接收延时信号计算出采样点的相位相干因子的计算过程要比现有技术中对ROI区域范围内的RF数据进行横向傅里叶变换简单的多,因此,本申请中的计算过程会更加简单,确定出最优声速的效率会更高。
在上述实施例的基础上,本实施例对技术方案作了进一步的说明和优化。
具体如下:
进一步的,各个采样点为每个声速下每条接收线上的不同深度的采样点;
具体的,在实际应用中接收孔径中的每个阵元会对应一条接收线,并且每条接收线上沿深度方向分布有不同深度的采样点,也即,本实施例中的各个采样点是由每个声速下位于每条接收线上不同深度处的采样点构成的,并且在一个声速下,各条接收线上的采样点数量相同,其中一条接收线上的各个采样点的深度不同,每条接收线上的各个采样点在深度方向的分布情况相同。
则,在进行延时时间计算时,本实施例以线阵为例进行说明,可参照图3,图3中以一条接收线上的第n个采样点到第i个阵元为例进行说明,第n个采样点到第i个阵元的延时时间可以由第n个采样点到第i个阵元之间的声程除以对应的声速得到,其中,第n个采样点到第i个阵元之间的声程为l1+l2,接收线上的每个采样点到各个阵元之间的延时时间均按照上述方法得到,从而得到各条接收线上的各采样点到各阵元的延时时间,进一步得到各个接收延时信号及对应的相位相干因子。
则,在上述实施例中的S15中根据各个相位相干因子,得到ROI区域的最优声速的过程可以包括以下步骤,具体请参照图4:
S150:对各条接收线上位于同一深度的采样点的相位相干因子求加权平均,得到每个声速下不同深度处的相位相干因子;
具体的,请参照图5,其中,列方向是采样点深度方向,横向是接收线排列方向,页方向是迭代声速排布方向,其中,图5中的Speed为声速。在接收线排列方向分布有各条接收线,在声速迭代方向分布有各个声速,在每个声速对应的各条接收线上分布有多个沿深度方向的各个采样点(如图5中的a部分),也即,在S14中得到的是三维数组,本实施例中通过沿接收线排列方向将每个声速下、各条接收线上处于同一深度的各个采样点的相位相干因子进行加权求平均,即可的到每个声速下不同深度处的相位相干因子,也即得到二维数组(如图5中的b部分)。
S151:将每个声速下不同深度处的相位相干因子求加权平均,得到每个声速的相位相干因子;
具体的,在得到每个声速下不同深度处的相位相干因子后,将每个声速下的各个深度处的相位相干因子进行加权平均,即可得到与相应声速对应的相位相干因子;也即,将图5中的b部分中每条线上的各个深度对应的相位相干因子进行加权平均,得到与对应声速的相位相干因子,例如共有L个声速,则得到与每个声速各自对应的相位相干因子,也即得到一个一维数组(如图5中的c部分)。
S152:将每个声速的相位相干因子进行拟合得到拟合曲线,并将拟合曲线的极大值对应的声速作为ROI区域的最优声速。
具体的,为了得到更加优选的声速,本实施例中可以将每个声速各自对应的相位相干因子进行曲线拟合,得到对应的拟合曲线,从该拟合曲线中选择出极大值对应的声速,并将该声速作为ROI区域的最优声速。
当然,在实际应用中还可以直接从每个声速各自对应的相位相干因子中选择出最大相位相干因子对应的声速,并将该声速作为ROI区域的最优声速。
更进一步的,上述实施例中的S14中的根据各个接收延时信号得到每个声速下各个采样点的相位相干因子的过程,具体可以包括以下步骤:
依据与每个声速下、各条接收线上每个采样点对应的各个接收延时信号及相位相干因子计算关系式,得到每个采样点的相位相干因子,其中,相位相干因子计算关系式为:
需要说明的是,在计算各个采样点的相位相干因子时,对于每个声速下的每条接收线,接收线上的每个采样点对应多个接收延时信号,则采样点的相位相干因子可以根据与该采样点对应的各个接收延时信号及上述相位相干因子计算关系式得到,其中,当时,该接收延时信号的符号位bi(n)为-1,当该接收延时信号的符号位bi(n)为+1。由上述相位相干因子计算关系式可知,本申请中的相位相干因子计算方式简单,并且在计算过程中不需要进行阈值设置,能够进一步降低人为因素的影响,使分析结果更加准确。
另外,还需要说明的是计算相位相干因子SCF的优势不仅在于计算量小,而且其可以兼顾算法在硬件环境下实施,例如,对于在硬件上实时计算接收延时信号的系统而言,SCF可以制作成查询表,从而避免平方运算,其中,正数的二进制格式符号位为0,负数的二进制格式符号位为1,将其符号位直接相加可以得到二进制格式符号位之和,可以将该和记为SB,将SB作为该表的查询地址,根据SB与SCF表达式的关系,可以知道该查询表一共有EleNum+1个元素,且数值关于EleNum/2这条直线对称,也就是说只需要存储EleNum/2个元素组成的数组即可,详细的LUT表(也即查询表)格式如图6所示。
例如,接收线上深度方向第n个采样点的各个阵元对应的接收延时信号均为正数,则采用二进制格式表示各个接收延时信号的符号位,则每个接收延时信号对应的符号位均为0,各个符号位之和SB=0,根据预设的对应关系在SCF表达式中每个接收延时信号的符号位均为1,因此,SB=0对应的SCF表达式中bi(n)均为1,此时SCF(n)等于1。另外,当接收线上深度方向第n个采样点的各个阵元对应的接收延时信号均为负数时,则采用二进制格式表示各个接收延时信号的符号位,则每个接收延时信号对应的符号位均为1,各个符号位之和SB=EleNum,根据预设的对应关系在SCF表达式中每个接收延时信号的符号位均为-1,因此,SB=EleNum对应的SCF表达式中bi(n)均为-1,此时SCF(n)也等于1,SCF表达式中bi(n)均为1对应的SCF值与SCF表达式中bi(n)均为-1对应的SCF值相等,也即,SB=0对应的SCF值与SB=EleNum对应的SCF值相同,按照上述方式,将SB的各个具体数值与SCF的值进行对应,形成查询表,并且查询表中各个元素的数值关于EleNum/2这条直线对称。
当然,若在上位机上实现相位相干因子的计算方法,则可以使用相干因子CF(Coherence Factor)及相位相干因子PCF等CF的变形版本来作为图像质量评估因子,其中,相干因子CF及相位相干因子PCF的计算方法可以参照现有技术中已有的计算方法,在实际应用中具体采用哪种本申请不作特殊限定,能够实现本申请的目的即可。
在上述实施例的基础上,本实施例对技术方案作了进一步优化和说明,具体如下:
进一步的,从成像区域中划分出的ROI区域为一个;则,该方法还可以包括以下步骤:
将ROI区域的最优声速作为成像区域的全局最优声速;
依据成像区域的全局最优声速进行波束合成,得到与成像区域对应的全局最优超声图像。
可以理解的是,在划分一个ROI区域时,需要以该ROI区域的最优声速表征整个成像区域的最优声速,因此,在得到ROI区域的最优声速后,直接将该最优声速作为成像区域的全局最优声速即可,并根据该全局最优声速进行波束合成,得到与该成像区域对应的全局最优超声图像,也即得到待成像图像的全局最优超声图像。
进一步的,从成像区域中划分出的ROI区域为多个,并且各个ROI区域沿成像区域的深度方向排列,各个ROI区域依次相连构成成像区域;则,该方法还可以包括以下步骤:
分别根据每个ROI区域各自对应的最优声速进行波束合成,得到与每个ROI区域各自对应的子超声图像;
将各个子超声图像进行拼接,得到与成像区域对应的全局最优超声图像。
可以理解的是,在对成像区域进行ROI区域划分时,可以沿深度方向将整个成像区域划分为多个ROI区域,并且各个ROI区域依次相连,如图7所示,将成像区域划分成3个ROI区域,分别为seg0、seg1和seg2,其中,start_seg1为seg1区域在深度方向的起始线位置,end_seg1为seg1区域在深度方向的结束线位置。当然,ROI区域的具体划分数量可以根据实际需要进行确定,本实施例不作特殊限定。
具体的,每个ROI区域均按照上述实施例中所提供的方法得到对应的最优声速,为了得到整个成像区域的全局最优图像,可以分别根据各个ROI区域的最优声速进行波束合成,得到与每个ROI区域对应的子超声图像,也即先得到与每个ROI区域各自对应的局部最优超声图像,然后将各个子超声图像(局部最优超声图像)按照对应的顺序进行拼接,即可得到与整个成像区域对应的全局最优超声图像。
由于相邻的两个ROI区域的声速不连续,若直接将相邻两个ROI区域对应的子超声图像相连接,则会造成两个子超声图像相连的位置呈现出不连贯的伪影,也就是不同声速下,图像会出现一定程度的拉伸或压缩,声速较小时图像会被压缩,声速较大时图像呈现出拉伸现象,例如图7中的seg0最优声速为1480m/s,seg1最优声速为1540m/s,seg0的底部会有seg1的图像出现,从而造成图像重复出现的现象。本实施例中为了解决图像的拉伸和压缩效应对图像拼接产生的影响,在根据每个ROI区域的最优声速进行波束合成时,可以对相应ROI区域的在深度方向的起始位置、结束位置和相邻采样点的间距进行调节,具体如下:
上述分别根据每个ROI区域各自对应的最优声速进行波束合成,得到与每个ROI区域各自对应的子超声图像的过程,具体可以为:
依据与第k个ROI区域对应的最优声速,按照预设规则对第k个ROI区域在深度方向的起始位置、结束位置和相邻采样点的间距进行调节;其中,k=1,2,3,…K,K为ROI区域总数量;
依据与第k个ROI区域对应的最优声速以及调节后的起始位置、结束位置和相邻采样点的间距进行波束合成,得到与第k个ROI区域对应的子超声图像;其中,预设规则为:
start_segk=depth_startk*speedoptk/1540;
end_segk=depth_endk*speedoptk/1540;
samplestep_segk=samplestep*speedoptk/1540;
其中,1540为理论上超声在软组织内传播的平均速度,depth_startk和depth_endk分别表示第k个ROI区域的起始位置和结束位置;start_segk和end_segk分别表示第k个ROI区域调节后的起始位置和结束位置;samplestep表示声速1540下相邻采样点的间距,具体可以为(1540/2)/Fs;speedoptk表示第k个ROI区域的最优声速;samplestep_segk表示与第k个ROI区域对应的调节后的相邻采样点的间距。
需要说明的是,按照上述方法根据每个ROI区域的最优声速进行波束合成时,对每个ROI区域的在深度方向的起始位置、结束位置和相邻采样点的间距进行调节,然后再将得到的与每个ROI区域对应的子超声图像拼接起来,从而得到更加精确的全局最优超声图像。另外,在拼接完成后将图像深度标尺(也即图像的深度方向点坐标)按照1540m/s的声速进行设置,以便很好的兼容后续的图像测量等操作,也即,在最后显示图像时均使用1540m/s的声速。
由于声速在不同组织介质中的传播速度不同,因此为了进一步为医生的诊断提供更加完善的参考数据,辅助医生识别不同深度部位的组织介质,本申请中还可以将整个成像区域作为一个ROI区域进行分析,并得到深度索引-局部声速图,以便辅助医生观察扫查面的解剖信息。具体的:
在对成像区域进行ROI区域划分时,将整个成像区域划分为一个ROI区域;
则,该方法还可以包括以下步骤:
根据每个声速下不同深度处的相位相干因子,得到与每个深度处的最大相位相干因子对应的目标声速,以得到与每个深度各自对应的最优声速;
依据与每个深度各自对应的最优声速建立深度与最优声速的关系曲线;
将深度与最优声速的关系曲线沿深度方向划分为M段,并计算出每一段的平均最优声速,M为多个;
依据各段的平均最优声速,得到每一段的局部声速;
依据与每段各自对应的局部声速,建立深度与局部声速的关系曲线。
需要说明的是,具体可以参照图5的b部分,在得到每个声速下不同深度处的相位相干因子后,可以将同一个深度处的各个相位相干因子中的最大相位相干因子作为与该深度对应的目标声速(也即与该深度对应的最优声速),从而得到每个深度的最优深度。然后根据每个深度各自的最优声速能够得到深度-最优声速曲线,再将该深度-最优声速曲线沿深度方向划分为M段,为了提高结果准确性,本实施例中在进行迭代声速时,迭代步距应小于10m/s,例如可以设置为2m/s,并且各段等间距划分,M可以采用较大的数值,例如在深度方向有多少个采样点就划分多少段,当然,也可以每一段包括多个采样点,具体可以根据实际需要进行确定,本申请不作特殊限定。具体的,根据每一段中各个采样点对应的声速计算出每一段的平均最优声速(各个声速的平均值),进一步得到每一段的局部声速,从而得到深度-局部声速曲线。
具体的,在将深度-最优声速曲线沿深度方向划分为M段之前,还可以先对深度-最优声速曲线进行平滑降噪处理,从而使得到的结果更加准确。另外,当每个采样间隔为第一段时,则该段的平均最优声速即为对应深度处对应的最优声速,若每一段包括多个采样间隔,则可以根据各个深度处的最优声速进行加权平均,得到对应段的平均最优声速,具体采用哪种方式得到每一段的平均最优声速,本申请不做特殊限定。
其中,依据各段的平均最优声速,得到每一段的局部声速的过程包括:
依据每段的平均最优声速以及局部声速计算关系式,得到各段的局部声速;其中,局部声速计算关系式为:
需要说明的是,由于在对第一段扫查过程中,超声波直接进入第一段,此时第一段的平均最优声速就是第一段的局部声速,而对于第二段或者之后的各段而言,超声波会经过当前段之前的各段后才能进入当前段,因此,当前段的平均最优声速与其成像介质的声速(也即局部声速)不相等。由和可知,每一段的平均最优声速均与位于其前面各段的局部声速有关,因此,可以根据和得到第m段的局部声速也即可以进一步根据得到的深度-最优声速曲线反推出每一段的局部声速,具体可以根据第m段的平均最优声速和第m-1段的平均最优声速计算得到第m段的局部声速,其中,每一段的平均最优声速均可以根据深度-最优声速曲线得到。在得到各段的局部声速后,建立深度与局部声速的关系曲线,例如将深度-最优声速曲线沿深度方向分成3段,具体请参照图8,其中,图中的曲线为深度-最优声速曲线,折线为深度-局部声速曲线,可得到图9中成像区域由三种声速分别为1480m/s、1540m/s以及1580m/s的物质组成。
在上述实施例的基础上,本发明实施例相应的提供了一种自适应匹配最优声速的超声成像装置,具体请参照图10。该装置包括:
获取模块21,用于获取待成像图像的通道数据;
划分模块22,用于依据预设边界信息从成像区域中划分出ROI区域;
延时模块23,用于在每个声速下,根据相应的声速得到ROI区域内每个采样点到各个阵元的延时时间,并依据每个延时时间从通道数据中得到每个采样点对应各个阵元的接收延时信号;
计算模块24,用于根据各个接收延时信号得到每个声速下各个采样点的相位相干因子;
分析模块25,用于根据各个相位相干因子,得到ROI区域的最优声速。
需要说明的是,本实施例中提供的自适应匹配最优声速的超声成像装置在进行超声成像自适应匹配最优声速时,通过获取模块21获取待成像图像的通道数据,通过划分模块22根据预设边界信息从成像区域中划分出ROI区域,然后延时模块23在每个声速下,根据对应声速计算出ROI区域中每个采样点到各个阵元的延时时间,从而得到与每个采样点对应的多个延时时间,并根据各个延时时间从通道数据中得到相应的接收延时信号,计算模块24再根据每个声速下与每个采样点对应的多个接收延时信号得到各声速下与每个采样点对应的相位相干因子,分析模块25进一步根据各个相位相干因子得到ROI区域的最优声速;本申请不需要进行横向傅里叶变换,且相位相干因子的计算过程简单、计算量小,因此本发明在使用过程中计算简单高效,有利于提高自适应匹配最优声速的效率。
另外,对于本实施例中所涉及到的自适应匹配最优声速的超声成像方法的具体介绍请参照上述方法实施例,本申请在此不再赘述。
在上述实施例的基础上,本发明实施例还提供了一种超声设备,该设备包括:
存储器,用于存储计算机程序;
处理器,用于执行计算机程序时实现如上述自适应匹配最优声速的超声成像方法的步骤。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
还需要说明的是,在本说明书中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其他实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (10)
1.一种自适应匹配最优声速的超声成像方法,其特征在于,包括:
获取待成像图像的通道数据;
依据预设边界信息从成像区域中划分出ROI区域;其中,所述预设边界信息包括横向起始百分比参数fraction_start、横向结束百分比参数fraction_end、深度方向起始边界线depth_start和深度方向结束边界线depth_end;
在每个声速下,根据相应的声速得到所述ROI区域内每个采样点到各个阵元的延时时间,并依据每个所述延时时间从所述通道数据中得到每个所述采样点对应各个所述阵元的接收延时信号;
根据各个所述接收延时信号得到每个所述声速下各个所述采样点的相位相干因子;
根据各个所述相位相干因子,得到所述ROI区域的最优声速;其中:
在不同声速下进行划分ROI区域时,依据对应声速cl对所划分的ROI区域在深度方向的起始边界线和结束边界线的位置进行调节,其中,调节后的ROI区域在深度方向的起始边界线位置为depth_startl=(cl/c0)*depth_start,结束边界线位置为depth_endl=(cl/c0)*depth_end,其中,c0为系统声速。
2.根据权利要求1所述的自适应匹配最优声速的超声成像方法,其特征在于,各个所述采样点为每个所述声速下每条接收线上的不同深度的采样点;
所述根据各个所述相位相干因子,得到所述ROI区域的最优声速的过程包括:
对各条所述接收线上位于同一深度的采样点的相位相干因子求加权平均,得到每个所述声速下不同深度处的相位相干因子;
将每个所述声速下不同深度处的相位相干因子求加权平均,得到每个所述声速的相位相干因子;
将每个所述声速的相位相干因子进行拟合得到拟合曲线,并将所述拟合曲线的极大值对应的声速作为所述ROI区域的最优声速。
4.根据权利要求1至3任意一项所述的自适应匹配最优声速的超声成像方法,其特征在于,所述ROI区域为一个;所述方法还包括:
将所述ROI区域的最优声速作为所述成像区域的全局最优声速;
依据所述成像区域的全局最优声速进行波束合成,得到与所述成像区域对应的全局最优超声图像。
5.根据权利要求1至3任意一项所述的自适应匹配最优声速的超声成像方法,其特征在于,所述ROI区域为多个,各个所述ROI区域沿所述成像区域的深度方向排列,各个所述ROI区域依次相连构成所述成像区域;所述方法还包括:
分别根据每个所述ROI区域各自对应的最优声速进行波束合成,得到与每个所述ROI区域各自对应的子超声图像;
将各个所述子超声图像进行拼接,得到与所述成像区域对应的全局最优超声图像。
6.根据权利要求5所述的自适应匹配最优声速的超声成像方法,其特征在于,所述分别根据每个所述ROI区域各自对应的最优声速进行波束合成,得到与每个所述ROI区域各自对应的子超声图像的过程为:
依据与第k个ROI区域对应的最优声速,按照预设规则对所述第k个ROI区域在深度方向的起始位置、结束位置和相邻采样点的间距进行调节;其中,k=1,2,3,…K,K为ROI区域总数量;
依据与所述第k个ROI区域对应的最优声速以及调节后的起始位置、结束位置和相邻采样点的间距进行波束合成,得到与所述第k个ROI区域对应的子超声图像;其中,所述预设规则为:
start_segk=depth_startk*speedoptk/1540;
end_segk=depth_endk*speedoptk/1540;
samplestep_segk=samplestep*speedoptk/1540;
其中,1540为理论上超声在软组织内传播的平均速度,depth_startk和depth_endk分别表示第k个ROI区域的起始位置和结束位置;start_segk和end_segk分别表示第k个ROI区域调节后的起始位置和结束位置;samplestep表示声速1540下相邻采样点的间距;speedoptk表示第k个ROI区域的最优声速;samplestep_segk表示与第k个ROI区域对应的调节后的相邻采样点的间距。
7.根据权利要求2所述的自适应匹配最优声速的超声成像方法,其特征在于,所述ROI区域为整个所述成像区域,所述方法还包括:
根据每个所述声速下不同深度处的相位相干因子,得到与每个所述深度处的最大相位相干因子对应的目标声速,以得到与每个所述深度各自对应的最优声速;
依据与每个所述深度各自对应的最优声速建立深度与最优声速的关系曲线;
将所述深度与最优声速的关系曲线沿深度方向划分为M段,并计算出每一段的平均最优声速,所述M为多个;
依据各段的平均最优声速,得到每一段的局部声速;
依据与每段各自对应的局部声速,建立深度与局部声速的关系曲线。
9.一种自适应匹配最优声速的超声成像装置,其特征在于,包括:
获取模块,用于获取待成像图像的通道数据;
划分模块,用于依据预设边界信息从成像区域中划分出ROI区域;其中,所述预设边界信息包括横向起始百分比参数fraction_start、横向结束百分比参数fraction_end、深度方向起始边界线depth_start和深度方向结束边界线depth_end;
延时模块,用于在每个声速下,根据相应的声速得到所述ROI区域内每个采样点到各个阵元的延时时间,并依据每个所述延时时间从所述通道数据中得到每个所述采样点对应各个所述阵元的接收延时信号;
计算模块,用于根据各个所述接收延时信号得到每个所述声速下各个所述采样点的相位相干因子;
分析模块,用于根据各个所述相位相干因子,得到所述ROI区域的最优声速;其中:
在不同声速下进行划分ROI区域时,依据对应声速cl对所划分的ROI区域在深度方向的起始边界线和结束边界线的位置进行调节,其中,调节后的ROI区域在深度方向的起始边界线位置为depth_startl=(cl/c0)*depth_start,结束边界线位置为depth_endl=(cl/c0)*depth_end,其中,c0为系统声速。
10.一种超声设备,其特征在于,包括:
存储器,用于存储计算机程序;
处理器,用于执行所述计算机程序时实现如权利要求1至8任意一项所述自适应匹配最优声速的超声成像方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911184446.3A CN110840484B (zh) | 2019-11-27 | 2019-11-27 | 自适应匹配最优声速的超声成像方法、装置及超声设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911184446.3A CN110840484B (zh) | 2019-11-27 | 2019-11-27 | 自适应匹配最优声速的超声成像方法、装置及超声设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110840484A CN110840484A (zh) | 2020-02-28 |
CN110840484B true CN110840484B (zh) | 2022-11-11 |
Family
ID=69605349
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911184446.3A Active CN110840484B (zh) | 2019-11-27 | 2019-11-27 | 自适应匹配最优声速的超声成像方法、装置及超声设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110840484B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112370073A (zh) * | 2020-10-12 | 2021-02-19 | 深圳华声医疗技术股份有限公司 | 超声声速校正方法、超声信号处理设备及存储介质 |
CN112641468B (zh) * | 2020-12-15 | 2022-11-04 | 青岛海信医疗设备股份有限公司 | 超声成像方法、超声图像处理方法、超声设备及存储介质 |
CN115712120A (zh) * | 2022-11-30 | 2023-02-24 | 武汉联影医疗科技有限公司 | 超声成像方法、装置、计算机设备和存储介质 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1974672B9 (en) * | 2007-03-28 | 2014-04-16 | Kabushiki Kaisha Toshiba | Ultrasonic imaging apparatus and ultrasonic velocity optimization method |
JP2010099452A (ja) * | 2008-09-25 | 2010-05-06 | Fujifilm Corp | 超音波診断装置及び超音波診断方法 |
JP5528083B2 (ja) * | 2009-12-11 | 2014-06-25 | キヤノン株式会社 | 画像生成装置、画像生成方法、及び、プログラム |
WO2012043569A1 (ja) * | 2010-09-27 | 2012-04-05 | 富士フイルム株式会社 | 超音波画像生成装置、超音波画像生成方法およびプログラム |
JP2012192133A (ja) * | 2011-03-18 | 2012-10-11 | Fujifilm Corp | 超音波診断装置および超音波画像生成方法 |
US9052268B2 (en) * | 2011-03-22 | 2015-06-09 | Fujifilm Corporation | Ultrasound diagnostic apparatus and method of producing ultrasound image |
JP5281107B2 (ja) * | 2011-03-25 | 2013-09-04 | 富士フイルム株式会社 | 超音波診断装置および超音波画像生成方法 |
JP2016539707A (ja) * | 2013-12-09 | 2016-12-22 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 画像情報に基づく画像合成 |
CN109754407B (zh) * | 2019-01-10 | 2021-06-01 | 青岛海信医疗设备股份有限公司 | 超声图像处理方法、装置及设备 |
-
2019
- 2019-11-27 CN CN201911184446.3A patent/CN110840484B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110840484A (zh) | 2020-02-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110840484B (zh) | 自适应匹配最优声速的超声成像方法、装置及超声设备 | |
JP4787358B2 (ja) | 超音波診断装置 | |
EP3581961A1 (en) | Method and apparatus for ultrasound imaging with improved beamforming | |
JP7358457B2 (ja) | 超音波画像による脂肪層の識別 | |
EP2633818B1 (en) | Ultrasonic diagnostic apparatus | |
US9081097B2 (en) | Component frame enhancement for spatial compounding in ultrasound imaging | |
EP3278734B1 (en) | Beamforming device, ultrasonic imaging device, and beamforming method allowing simple spatial smoothing operation | |
KR20150112057A (ko) | 공간 일관성 기초 초음파 신호 처리 모듈 및 그에 의한 초음파 신호 처리 방법 | |
Long et al. | Incoherent clutter suppression using lag-one coherence | |
JP2022518345A (ja) | 低コスト、低サイズおよび低電力で超音波イメージングにおける画質を維持する方法 | |
Jansen et al. | Enhanced radon domain beamforming using deep-learning-based plane wave compounding | |
JP2020523143A (ja) | 超音波画像を処理するための方法及びシステム | |
CN111388012B (zh) | 用于检测组织硬度的方法、设备及系统 | |
Du et al. | A novel transcranial ultrasound imaging method with diverging wave transmission and deep learning approach | |
JP7171625B2 (ja) | 超音波画像を処理するための方法及びシステム | |
US20230061869A1 (en) | System and methods for beamforming sound speed selection | |
Kim et al. | Learning based approach for speed-of-sound adaptive rx beamforming | |
Sharifzadeh et al. | Phase Aberration Correction: A Deep Learning-Based Aberration to Aberration Approach | |
Kwon et al. | Estimation and suppression of side lobes in medical ultrasound imaging systems | |
CN112370073A (zh) | 超声声速校正方法、超声信号处理设备及存储介质 | |
EP3603526B1 (en) | Ultrasonic diagnostic device and method for controlling ultrasonic diagnostic device | |
EP4390841A1 (en) | Image acquisition method | |
EP4258016A1 (en) | A computer-implemented method for beamforming of ultrasound channel data | |
US20230360225A1 (en) | Systems and methods for medical imaging | |
JP2018157982A (ja) | 超音波診断装置及びプログラム |
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 |