超声探头几何参数自动校正方法和装置
技术领域
本发明涉及超声诊断设备技术领域,特别是涉及一种超声探头几何参数自动校正方法和装置。
背景技术
医用超声成像设备因其无创、实时、便捷等优点,在现代医学成像中占有重要位置。优质的超声图像能为临床提供更准确的诊断信息。超声图像的质量与超声探头、主机***以及显示器等密切相联。其中,超声探头作为超声的发射和接收组件,是超声成像设备中重要组成部分。描述超声探头几何特征的参数,在超声成像设备中被用于扫描控制、发射和接收聚焦、DSC变换等过程。超声探头几何参数的准确性,是影响超声图像的质量的关键因素之一。如果超声探头的实际几何参数与超声成像设备中所使用的参数有较大偏离,会导致发射和接收聚焦不能达到最优,DSC变换存在较大误差,从而降低图像的横向分辨率和几何位置精度,即降低了图像质量。为了获得最优的图像质量,在超声成像设备中需要使用能够提供最优图像质量的超声探头几何参数。
不同类型的探头具有不同的几何参数。即使是同一种探头同一批次的两只探头,由于探头在生产过程中不可避免地引入误差,也会存在几何参数上微小差别。探头生产商提供给超声设备生产厂商的探头规格书会包含探头的几何参数,但该参数只是统计均值,并非所有探头严格符合该数值。在超声设备生产过程中,需要对每只探头在均值参数的基础上进行校正,达到最优参数,以提供最优质量的图像。一种解决方案是,人为输入各参数的经验值,根据某种法则,判定最优的一组几何参数,作为该探头的几何参数。但这种解决方案需要耗费大量人力和时间,在实际生产过程中是不可取的。另一种解决方案是提供一种装置,能够自动搜索判定出最优的一组几何参数,高效完成该任务。
专利CN102871685B提供了一种超声探头几何参数的校正方法,该方法采集设定的几何参数对应的感兴趣区域内I、Q数据,并对I、Q数据求幅度后进行二维傅里叶变换,求取某一范围内的频率能量值,以能量值作为判断依据找到最优几何参数。该方法确实能够解决上述问题。实际情况中,中间过程的I、Q数据对超声诊断设备用户是不可见的,而最终显示在显示设备上的超声图像对用户来说更为直观。此外,超声图像的分辨率和几何位置精度也是通过最终显示的图像而非中间过程的数据来进行测量的。根据最终显示的超声图像的质量作为判断依据更适合用于自动寻找最优的探头几何参数。
发明内容
本发明的目的在于克服上述不足,提供一种超声探头几何参数的自动校正方法和装置,使超声诊断设备能使用最优超声探头几何参数得到最优质量的图像。
本发明的目的是通过以下技术方案实现的:
一种超声探头几何参数的自动校正方法,其特殊之处在于:它包括如下步骤:
步骤S1:确定超声探头几何参数列表PL,记PL={Pk|k=1,2,L,N};
步骤S2:超声成像设备使用默认的超声探头几何参数Pdef,连接超声探头与超声成像设备,将超声探头固定于超声体模上方,使超声成像设备处于正常工作状态,调节设备参数,使显示设备能清晰显示体模上靶点的B模式图像,设置感兴趣区域,使得所述感兴趣区域包含一个或多个靶点;
步骤S3:依次读取几何参数列表PL中每一项Pk,将该项参数用于超声成像设备对超声体模进行成像,其他参数设置与步骤S2中相同,然后采集显示图像中感兴趣区域图像数据,根据感兴趣区域图像计算图像优劣度量Mk;
步骤S4:从所有参数项求得的感兴趣区域图像优劣度量集合ML={Mk|k=1,2,L,N},找到最优的图像优劣度量Mbest,该最优图像优劣度量Mbest对应的参数Pbest即为该探头校正得到的最优几何参数,然后将该最优几何参数作为所述探头的几何参数应用于所述超声成像设备的成像;
在上述步骤中,所述超声探头几何参数为曲率半径,或者为阵元间距,或者为扫描宽度,或者为曲率半径和阵元间距,或者为扫描宽度和阵元间距。
优选地,步骤S2中所述默认的超声探头几何参数Pdef采用探头供应商提供的规格值。
优选地,步骤S2中所述超声体模为标准超声仿组织体模或自制体模,所述超声体模包含若干靶线。
进一步地,步骤S3中所述图像优劣度量Mk采用所述感兴趣区域中所有靶点的面积之和来表示,所述面积为分割结果中该靶点的像素个数,所述根据感兴趣区域图像计算图像优劣度量Mk,具体为:在感兴趣区域图像中寻找到所有靶点,并计算所有靶点的面积之和;步骤S4中所述最优的图像优劣度量等同于该所有靶点所有面积之和中的最小值。
进一步地,步骤S3中所述图像优劣度量Mk采用感兴趣区域中所有靶点的平均面积来表示,所述面积为分割结果中该靶点的像素个数,所述根据感兴趣区域图像计算图像优劣度量Mk,具体为:在感兴趣区域图像中寻找到所有靶点,并计算所有靶点的平均面积;步骤S4中所述最优的图像优劣度量Mbest为该所有靶点平均面积中的最小值。
更进一步地,所述寻找所有靶点过程具体为:以像素值大于某个阈值为准则对图像进行分割,分割的结果视为靶点;所述阈值人为设定或自适应确定。
本发明还提供一种超声探头几何参数自动校正装置,其连接到超声成像设备,包括几何参数列表计算及存储模块、采集模块、图像优劣度计算模块和参数选取模块,其特殊之处在于:
所述超声探头几何参数为曲率半径,或者为阵元间距,或者为扫描宽度,或者为曲率半径和阵元间距,或者为扫描宽度和阵元间距;
所述几何参数列表计算及存储模块,用于计算待校正几何参数列表并进行存储;
所述采集模块,用于依次读取几何参数列表并将各几何参数设置到超声成像设备,采集到图像;
所述图像优劣度计算模块,用于在所述采集到的图像中寻找到靶点,并计算图像优劣度;
所述参数选取模块,用于在所有参数项求得的图像优劣度量集合中选择最优的图像优劣度量及其对应的参数,并将该最优几何参数作为所述探头的几何参数应用于所述超声成像设备的成像。
本发明提供的超声探头几何参数的自动校正方法和装置适用于探头,包括并不限于凸阵、微凸阵、线阵、相控阵。
本发明的有益效果是:
本发明提供的超声探头几何参数的自动校正方法和装置,在探头规格中的几何参数基础上,对每只探头都能自动自适应寻找到提供最优图像质量的几何参数,将该最优几何参数提供给超声成像设备,可以得到最佳分辨率和几何位置精度的超声图像,从而提高图像质量。
附图说明
图1为本发明提供的一种超声探头几何参数的自动校正方法的流程图;
图2为本发明提供的一种超声探头几何参数的自动校正装置的结构示意图;
图3为实施例一的超声探头几何参数列表;
图4为实施例一采用默认的超声探头几何参数对体模的成像及设置的感兴趣区域;
图5为实施例一采用探头几何参数列表P1项对体模成像的感兴趣区域图像;
图6为对图5的分割结果;
图7为实施例二的超声探头几何参数列表;
图8为实施例二采用默认的超声探头几何参数对体模的成像及设置的感兴趣区域;
图9为实施例二采用探头几何参数列表P1项对体模成像的感兴趣区域图像;
图10为对图9的分割结果;
图11为实施例三的超声探头几何参数列表;
图12为实施例三采用默认的超声探头几何参数对体模的成像及设置的感兴趣区域;
图13为实施例三采用探头几何参数列表P1,1项对体模成像的感兴趣区域图像;
图14为对图13的分割结果。
图中,几何参数列表计算及存储模块110、采集模块120、图像优劣度计算模块130、参数选取模块140。
具体实施方式
下面详细描述本发明的实施例。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
实施例一
下面以几何参数为超声探头曲率半径,其出厂规格为60mm的凸阵对本发明实施例的超声探头几何参数的自动校正方法进行详细说明。
本发明提供的一种超声探头几何参数的自动校正方法,包括如下步骤:
步骤S1:确定超声探头几何参数列表PL={Pk|k=1,2,L,N},本实施例中Pk仅包含曲率半径取值Ri,即Pk={Ri};
参数列表PL中的各项Ri值由探头的工程参数决定。可以由如下两种方式产生:
方式1:Ri直接由经验参数给出;
方式2:给出Ri取值范围[Rmin,Rmax]及步进值Rstep,生成Ri的值:
Ri=Rmin+(i-1)*Rstep (1)
Rmin≤Ri≤Rmax,Rmin按经验可取Rmin=0.9*Rdef,Rmax按经验可取Rmax=1.1*Rdef,Rdef为供应商提供的探头曲率半径出厂规格值。步进值Rstep按经验可按式(2)取值:
曲率半径出厂规格为60mm的凸阵探头Rdef=60mm,使用方式2确定Ri过程如下:确定曲率半径取值范围为[Rmin,Rmax]=[54,66],步进值Rstep=0.1,则由式(1)得
Ri=Rmin+(i-1)*Rstep=54+(i-1)*0.1,i=1,2,L,N,N=121 (3)
此时参数列表PL各项值如图3所示。参数列表PL取值并不限于此。
步骤S2:超声成像设备使用默认的超声探头几何参数Pdef,连接超声探头与超声成像设备,将超声探头固定于超声体模上方,使超声成像设备处于正常工作状态,调节设备参数,使显示设备能清晰显示体模上靶点的B模式图像,设置感兴趣区域;
默认的超声探头几何参数Pdef采用探头供应商提供的规格值,即取Pdef={Rdef=60}。此处使用标准仿组织体模进行成像。调整超声成像设备参数如深度、焦点位置、增益、TGC等,使其对体模上的靶点清晰成像,成像结果如图4所示,靶点处具有高亮度。设置感兴趣区域如图4中的虚线方框所示。
步骤S3:依次读取几何参数列表PL中每一项Pk(k=1,2,L,N),将该项参数用于超声成像设备对超声体模进行成像,其他参数设置同步骤S2中相同。采集显示图像中感兴趣区域图像数据,根据感兴趣区域图像计算图像优劣度量Mk(k=1,2,L,N);
读取P1={R1=54},将此值用于对体模成像,得到感兴趣区域图像如图5所示。像素值阈值Th按经验取值128,使用像素值阈值Th对感兴趣区域图像进行分割,分割结果如图6所示,白色区域为靶点部分。一般地,体模中含有若干靶点,此时分割结果中会包含所有靶点。
计算图像优劣度量可采用如下两种方式:
方式1:求感兴趣区域中所有靶点面积之和;
方式2:求感兴趣区域中所有靶点的平均面积。
某一靶点的面积定义为分割结果中该靶点的像素个数。
选取上述任一方式求得P1对应的图像优劣度量M1。
依次读取P2,P3,…,P121,按照读取P1后相同的操作计算得到对应的图像优劣度量M2,M3,…,M121。
步骤S4:对所有参数项求得的图像优劣度量集合ML={M1,M2,L,MN},找到最优的图像优劣度量Mbest,该最优图像优劣度量Mbest对应的几何参数Pbest即为该探头校正得到的最优几何参数。并将该最优几何参数作为所述探头的几何参数应用于所述超声成像设备的成像。
如果计算图像优劣度量是按照上述步骤S3中五种方式中任一种进行的,那么此步骤最优图像优劣度量Mbest即为ML中最小的值,其对应的几何参数Pbest即为该探头校正得到的最优曲率半径参数。
几何参数为扫描宽度时的超声探头几何参数的自动校正方法的实施与实施例一相同,不再一一详细描述。
实施例二
下面以几何参数为超声探头阵元间距,其出厂规格为0.78mm的凸阵对本发明实施例的超声探头几何参数的自动校正方法进行详细说明。
本发明提供的一种超声探头几何参数的自动校正方法,包括如下步骤:
步骤S1:确定超声探头几何参数列表PL={Pk|k=1,2,L,N},本实施例中Pk仅包含阵元间距取值Ci,即Pk={Ci};
参数列表PL中的各项Ci值由探头的工程参数决定。可以由如下两种方式产生:
方式1:Ci直接由经验参数给出;
方式2:给出Ci取值范围[Cmin,Cmax]及步进值Cstep,生成Ci的值:
Ci=Cmin+(i-1)*Cstep (4)
Cmin≤Ci≤Cmax,Cmin按经验可取Cmin=0.9*Cdef,Cmax按经验可取Cmax=1.1*Cdef,Cdef为供应商提供的探头曲率半径出厂规格值。步进值Cstep按经验可按式(5)取值:
阵元间距出厂规格为0.78mm的凸阵探头Cdef=0.78mm,使用方式2确定Ci过程如下:确定曲率半径取值范围为[Cmin,Cmax]=[0.702,0.858],步进值Cstep=0.001,则由式(4)得
Ci=Cmin+(i-1)*Cstep=0.702+(i-1)*0.001,i=1,2,L,N,N=157 (6)
此时参数列表PL各项值如图7所示。参数列表PL取值并不限于此。
步骤S2:超声成像设备使用默认的超声探头几何参数Pdef,连接超声探头与超声成像设备,将超声探头固定于超声体模上方,使超声成像设备处于正常工作状态,调节设备参数,使显示设备能清晰显示体模上靶点的B模式图像,设置感兴趣区域;
默认的超声探头几何参数Pdef采用探头供应商提供的规格值,即取Pdef={Cdef=0.78}。此处使用标准仿组织体模进行成像。调整超声成像设备参数如深度、焦点位置、增益、TGC等,使其对体模上的靶点清晰成像,成像结果如图8所示,靶点处具有高亮度。设置感兴趣区域如图8中的虚线方框所示。
步骤S3:依次读取几何参数列表PL中每一项Pk(k=1,2,L,N),将该项参数用于超声成像设备对超声体模进行成像,其他参数设置同步骤S2中相同。采集显示图像中感兴趣区域图像数据,根据感兴趣区域图像计算图像优劣度量Mk(k=1,2,L,N);
读取P1={C1=0.702},将此值用于成像,得到感兴趣区域图像如图9所示。像素值阈值Th按经验取值128,使用像素值阈值Th对感兴趣区域图像进行分割,分割结果如图10所示,白色区域为靶点部分。一般地,体模中含有若干靶点,此时分割结果中会包含所有靶点。
计算图像优劣度量可采用如下两种方式:
方式1:求感兴趣区域中所有靶点面积之和;
方式2:求感兴趣区域中所有靶点的平均面积。
某一靶点的面积定义为分割结果中该靶点的像素个数。
选取上述任一方式求得P1对应的图像优劣度量M1。
依次读取P2,P3,…,P157,按照读取P1后相同的操作计算得到对应的图像优劣度量M2,M3,…,M157。
步骤S4:对所有参数项求得的图像优劣度量集合ML={M1,M2,L,MN},找到最优的图像优劣度量Mbest,该最优图像优劣度量Mbest对应的几何参数Pbest即为该探头校正得到的最优几何参数。并将该最优几何参数作为所述探头的几何参数应用于所述超声成像设备的成像。
如果计算图像优劣度量是按照上述步骤S3中五种方式中任一种进行的,那么此步骤最优图像优劣度量Mbest即为ML中最小的值,其对应的几何参数Pbest即为该探头校正得到的最优阵元间距参数。
实施例三
下面以几何参数为超声探头曲率半径和阵元间距,曲率半径出厂规格为60mm,阵元间距出厂规格为0.78mm的凸阵对本发明实施例的超声探头几何参数的自动校正方法进行详细说明。
本发明提供的一种超声探头几何参数的自动校正方法,包括如下步骤:
步骤S1:确定超声探头几何参数列表PL={Pi,j|i=1,2,L,N1,j=1,2,L,N2},本实施例中Pi,j包含曲率半径取值Ri和阵元间距取值Cj,即Pi,j={Ri,Cj};
参数列表PL中的各项Ri值和Cj值由探头的工程参数决定。可以由如下两种方式产生:
方式1:Ri、Cj直接由经验参数给出;
方式2:给出Ri取值范围[Rmin,Rmax]及步进值Rstep,生成Ri的值;给出Cj取值范围[Cmin,Cmax]及步进值Cstep,生成Cj的值:
Ri=Rmin+(i-1)*Rstep (7)
Cj=Cmin+(j-1)*Cstep (8)
Rmin≤Ri≤Rmax,Rmin按经验可取Rmin=0.9*Rdef,Rmax按经验可取Rmax=1.1*Rdef,Rdef为供应商提供的探头曲率半径出厂规格值。步进值Rstep按经验可按式(9)取值:
Cmin≤Ci≤Cmax,Cmin按经验可取Cmin=0.9*Cdef,Cmax按经验可取Cmax=1.1*Cdef,Cdef为供应商提供的探头曲率半径出厂规格值。步进值Cstep按经验可按式(10)取值:
曲率半径出厂规格为60mm的凸阵探头Rdef=60mm,使用方式2确定Ri过程如下:确定曲率半径取值范围为[Rmin,Rmax]=[54,66],步进值Rstep=0.1,则由式(7)得
Ri=Rmin+(i-1)*Rstep=54+(i-1)*0.1,i=1,2,L,N1,N1=121 (11)
阵元间距出厂规格为0.78mm的凸阵探头Cdef=0.78mm,使用方式2确定Cj过程如下:确定曲率半径取值范围为[Cmin,Cmax]=[0.702,0.858],步进值Cstep=0.001,则由式(8)得
Ci=Cmin+(i-1)*Cstep=0.702+(i-1)*0.001,i=1,2,L,N2,N2=157 (12)
此时参数列表PL各项值如图11所示。参数列表PL取值并不限于此。
步骤S2:超声成像设备使用默认的超声探头几何参数Pdef,连接超声探头与超声成像设备,将超声探头固定于超声体模上方,使超声成像设备处于正常工作状态,调节设备参数,使显示设备能清晰显示体模上靶点的B模式图像,设置感兴趣区域;
默认的超声探头几何参数Pdef采用探头供应商提供的规格值,即取Pdef={Rdef=60,Cdef=0.78}。此处使用标准仿组织体模进行成像。调整超声成像设备参数如深度、焦点位置、增益、TGC等,使其对体模上的靶点清晰成像,成像结果如图12所示,靶点处具有高亮度。设置感兴趣区域如图12中的虚线方框所示。
步骤S3:依次读取几何参数列表PL中每一项Pi,j(i=1,2,L,N1,j=1,2,L,N2),将该项参数用于超声成像设备对超声体模进行成像,其他参数设置同步骤S2中相同。采集显示图像中感兴趣区域图像数据,根据感兴趣区域图像计算图像优劣度量Mi,j(i=1,2,L,N1,j=1,2,L,N2);
读取P1,1={R1=54,C1=0.702},将此值用于成像,得到感兴趣区域图像如图13所示。像素值阈值Th按经验取值128,使用像素值阈值Th对感兴趣区域图像进行分割,分割结果如图14所示,白色区域为靶点部分。一般地,体模中含有若干靶点,此时分割结果T中会包含所有靶点,此时分割结果中会包含所有靶点。
计算图像优劣度量可采用如下两种方式:
方式1:求感兴趣区域中所有靶点面积之和;
方式2:求感兴趣区域中所有靶点的平均面积。
某一靶点的面积定义为分割结果中该靶点的像素个数。
选取上述任一方式求得P1,1对应的图像优劣度量M1,1。
依次读取P1,2,P1,3,…,P121,157,按照读取P1后相同的操作计算得到对应的图像优劣度量M1,2,M1,3,…,M121,157。
步骤S4:对所有参数项求得的图像优劣度量集合
ML={Mi,j|i=1,2,L,N1,j=1,2,L,N2},找到最优的图像优劣度量Mbest,该最优图像优劣度量Mbest对应的几何参数Pbest即为该探头校正得到的最优几何参数。并将该最优几何参数作为所述探头的几何参数应用于所述超声成像设备的成像。如果计算图像优劣度量是按照上述步骤S3中五种方式中任一种进行的,那么此步骤最优图像优劣度量Mbest即为ML中最小的值,其对应的几何参数Pbest即为该探头校正得到的最优曲率半径和阵元间距参数。
几何参数为扫描宽度和阵元间距时的超声探头几何参数的自动校正方法的实施与实施例三相同,不再一一详细描述。
实施例四
相应地,基于同一发明构思,本发明还提供一种超声探头几何参数的自动校正装置,其连接到超声成像设备,如图2所示,包括几何参数列表计算及存储模块110、采集模块120、图像优劣度计算模块130和参数选取模块140,其中:
所述超声探头几何参数或者为曲率半径,或者为阵元间距,或者为扫描宽度,或者为曲率半径和阵元间距,或者为扫描宽度和阵元间距;
所述几何参数列表计算及存储模块110,用于计算待校正几何参数列表并进行存储;
所述采集模块120,用于依次读取几何参数列表并将各几何参数设置到超声成像设备,采集到图像;
所述图像优劣度计算模块130,用于在所述采集到的图像中寻找到靶点,并计算图像优劣度;
所述参数选取模块140,用于在所有参数项求得的图像优劣度量集合中选择最优的图像优劣度量及其对应的参数,并将该最优几何参数作为所述探头的几何参数应用于所述超声成像设备的成像。
本发明实施例四中的超声探头几何参数的自动校正装置,其校正工作过程与实施例一、二、三相同,因此,在本发明实施例中不再一一详细描述。