CN104997533A - 超声探头几何参数自动校正方法和装置 - Google Patents

超声探头几何参数自动校正方法和装置 Download PDF

Info

Publication number
CN104997533A
CN104997533A CN201510348976.2A CN201510348976A CN104997533A CN 104997533 A CN104997533 A CN 104997533A CN 201510348976 A CN201510348976 A CN 201510348976A CN 104997533 A CN104997533 A CN 104997533A
Authority
CN
China
Prior art keywords
image
geometric parameter
ultrasonic probe
parameter
area
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.)
Granted
Application number
CN201510348976.2A
Other languages
English (en)
Other versions
CN104997533B (zh
Inventor
杜春宁
胡立钢
彭江
段西尧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong Qiyou Medical Co.,Ltd.
Original Assignee
WUHAN SOUNSHINE ELECTRONIC ENGINEERING CO LTD
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by WUHAN SOUNSHINE ELECTRONIC ENGINEERING CO LTD filed Critical WUHAN SOUNSHINE ELECTRONIC ENGINEERING CO LTD
Priority to CN201510348976.2A priority Critical patent/CN104997533B/zh
Publication of CN104997533A publication Critical patent/CN104997533A/zh
Application granted granted Critical
Publication of CN104997533B publication Critical patent/CN104997533B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4444Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4444Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
    • A61B8/4461Features of the scanning mechanism, e.g. for moving the transducer within the housing of the probe
    • A61B8/4466Features of the scanning mechanism, e.g. for moving the transducer within the housing of the probe involving deflection of the probe
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4483Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
    • A61B8/4488Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer the transducer being a phased array
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/58Testing, adjusting or calibrating the diagnostic device

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Gynecology & Obstetrics (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明提供一种超声探头几何参数的自动校正方法,它包括如下步骤:确定超声探头几何参数列表;使用默认的超声探头几何参数,调节设备参数,设置感兴趣区域;读取几何参数列表中每一项,将该项参数用于超声成像设备对超声体模进行成像,采集显示图像中感兴趣区域图像数据;根据感兴趣区域图像计算图像优劣度量;另外本发明还提供与其相关的超声探头几何参数的自动校正装置;本发明在探头规格中的几何参数基础上,对每只探头都能自动自适应寻找到提供最优图像质量的几何参数,将该最优几何参数提供给超声成像设备,可以得到最佳分辨率和几何位置精度的超声图像,从而提高图像质量。

Description

超声探头几何参数自动校正方法和装置
技术领域
本发明涉及超声诊断设备技术领域,特别是涉及一种超声探头几何参数自动校正方法和装置。
背景技术
医用超声成像设备因其无创、实时、便捷等优点,在现代医学成像中占有重要位置。优质的超声图像能为临床提供更准确的诊断信息。超声图像的质量与超声探头、主机系统以及显示器等密切相联。其中,超声探头作为超声的发射和接收组件,是超声成像设备中重要组成部分。描述超声探头几何特征的参数,在超声成像设备中被用于扫描控制、发射和接收聚焦、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)取值:
R s t e p = { 0.05 , R d e f ≤ 20 0.1 , R d e f > 20 - - - ( 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)取值:
R s t e p = { 0.0005 , C d e f ≤ 0.5 0.001 , C d e f > 0.5 - - - ( 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)取值:
R s t e p = { 0.05 , R d e f ≤ 20 0.1 , R d e f > 20 - - - ( 9 )
Cmin≤Ci≤Cmax,Cmin按经验可取Cmin=0.9*Cdef,Cmax按经验可取Cmax=1.1*Cdef,Cdef为供应商提供的探头曲率半径出厂规格值。步进值Cstep按经验可按式(10)取值:
R s t e p = { 0.0005 , C d e f ≤ 0.5 0.001 , C d e f > 0.5 - - - ( 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,用于在所有参数项求得的图像优劣度量集合中选择最优的图像优劣度量及其对应的参数,并将该最优几何参数作为所述探头的几何参数应用于所述超声成像设备的成像。
本发明实施例四中的超声探头几何参数的自动校正装置,其校正工作过程与实施例一、二、三相同,因此,在本发明实施例中不再一一详细描述。

Claims (7)

1.一种超声探头几何参数的自动校正方法,其特征在于:它包括如下步骤:
步骤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即为该探头校正得到的最优几何参数,然后将该最优几何参数作为所述探头的几何参数应用于所述超声成像设备的成像;
在上述步骤中,所述超声探头几何参数为曲率半径,或者为阵元间距,或者为扫描宽度,或者为曲率半径和阵元间距,或者为扫描宽度和阵元间距。
2.根据权利要求1所述的超声探头几何参数的自动校正方法,其特征在于:步骤S2中所述默认的超声探头几何参数Pdef采用探头供应商提供的规格值。
3.根据权利要求1所述的超声探头几何参数的自动校正方法,其特征在于:步骤S2中所述超声体模为标准超声仿组织体模或自制体模,所述超声体模包含若干靶线。
4.根据权利要求1所述的超声探头几何参数的自动校正方法,其特征在于:步骤S3中所述图像优劣度量Mk采用所述感兴趣区域中所有靶点的面积之和来表示,所述面积为分割结果中该靶点的像素个数,所述根据感兴趣区域图像计算图像优劣度量Mk,具体为:在感兴趣区域图像中寻找到所有靶点,并计算所有靶点的面积之和;步骤S4中所述最优的图像优劣度量等同于该所有靶点所有面积之和中的最小值。
5.根据权利要求1所述的超声探头几何参数的自动校正方法,其特征在于:步骤S3中所述图像优劣度量Mk采用感兴趣区域中所有靶点的平均面积来表示,所述面积为分割结果中该靶点的像素个数,所述根据感兴趣区域图像计算图像优劣度量Mk,具体为:在感兴趣区域图像中寻找到所有靶点,并计算所有靶点的平均面积;步骤S4中所述最优的图像优劣度量Mbest为该所有靶点平均面积中的最小值。
6.根据权利要求4或5所述的超声探头几何参数的自动校正方法,其特征在于:所述寻找所有靶点过程具体为:以像素值大于某个阈值为准则对图像进行分割,分割的结果视为靶点;所述阈值人为设定或自适应确定。
7.一种超声探头几何参数自动校正装置,其连接到超声成像设备,包括几何参数列表计算及存储模块(110)、采集模块(120)、图像优劣度计算模块(130)和参数选取模块(140),其特征在于:
所述超声探头几何参数为曲率半径,或者为阵元间距,或者为扫描宽度,或者为曲率半径和阵元间距,或者为扫描宽度和阵元间距;
所述几何参数列表计算及存储模块(110),用于计算待校正几何参数列表并进行存储;
所述采集模块(120),用于依次读取几何参数列表并将各几何参数设置到超声成像设备,采集到图像;
所述图像优劣度计算模块(130),用于在所述采集到的图像中寻找到靶点,并计算图像优劣度;
所述参数选取模块(140),用于在所有参数项求得的图像优劣度量集合中选择最优的图像优劣度量及其对应的参数,并将该最优几何参数作为所述探头的几何参数应用于所述超声成像设备的成像。
CN201510348976.2A 2015-06-23 2015-06-23 超声探头几何参数自动校正方法和装置 Active CN104997533B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510348976.2A CN104997533B (zh) 2015-06-23 2015-06-23 超声探头几何参数自动校正方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510348976.2A CN104997533B (zh) 2015-06-23 2015-06-23 超声探头几何参数自动校正方法和装置

Publications (2)

Publication Number Publication Date
CN104997533A true CN104997533A (zh) 2015-10-28
CN104997533B CN104997533B (zh) 2017-11-14

Family

ID=54370598

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510348976.2A Active CN104997533B (zh) 2015-06-23 2015-06-23 超声探头几何参数自动校正方法和装置

Country Status (1)

Country Link
CN (1) CN104997533B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108464845A (zh) * 2018-01-22 2018-08-31 苏州佳世达电通有限公司 一种超声波探头的异常侦测方法和超声波诊断系统
CN111212223A (zh) * 2020-01-10 2020-05-29 深圳奥比中光科技有限公司 一种设定成像参数的方法、系统及计算机可读存储介质
WO2020192403A1 (zh) * 2019-03-28 2020-10-01 无锡声美达医学技术有限公司 一种仪器测试方法、装置及设备
CN113545798A (zh) * 2020-04-23 2021-10-26 深圳市理邦精密仪器股份有限公司 一种颅内压测量方法及超声设备
CN114098809A (zh) * 2021-11-29 2022-03-01 无锡海鹰电子医疗系统有限公司 一种超声诊断仪发射系统及其快速优化方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008144452A1 (en) * 2007-05-16 2008-11-27 Verathon Inc. System and method for ultrasonic harmonic imaging
CN102028499A (zh) * 2010-12-24 2011-04-27 飞依诺科技(苏州)有限公司 超声探头自适应成像系统
CN202069602U (zh) * 2011-04-14 2011-12-14 武汉超信电子工程有限公司 具有机器参数校正功能的超声诊断仪
CN102871685A (zh) * 2012-09-14 2013-01-16 飞依诺科技(苏州)有限公司 超声探头几何参数的校正方法和装置及系统
CN103356241A (zh) * 2013-08-02 2013-10-23 中国十九冶集团有限公司职工医院 二维超声设备成像质量评估系统及评估方法
CN103635148A (zh) * 2011-07-07 2014-03-12 浜松光子学株式会社 生物体测量装置以及图像制作方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008144452A1 (en) * 2007-05-16 2008-11-27 Verathon Inc. System and method for ultrasonic harmonic imaging
CN102028499A (zh) * 2010-12-24 2011-04-27 飞依诺科技(苏州)有限公司 超声探头自适应成像系统
CN202069602U (zh) * 2011-04-14 2011-12-14 武汉超信电子工程有限公司 具有机器参数校正功能的超声诊断仪
CN103635148A (zh) * 2011-07-07 2014-03-12 浜松光子学株式会社 生物体测量装置以及图像制作方法
CN102871685A (zh) * 2012-09-14 2013-01-16 飞依诺科技(苏州)有限公司 超声探头几何参数的校正方法和装置及系统
CN103356241A (zh) * 2013-08-02 2013-10-23 中国十九冶集团有限公司职工医院 二维超声设备成像质量评估系统及评估方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108464845A (zh) * 2018-01-22 2018-08-31 苏州佳世达电通有限公司 一种超声波探头的异常侦测方法和超声波诊断系统
WO2020192403A1 (zh) * 2019-03-28 2020-10-01 无锡声美达医学技术有限公司 一种仪器测试方法、装置及设备
CN111212223A (zh) * 2020-01-10 2020-05-29 深圳奥比中光科技有限公司 一种设定成像参数的方法、系统及计算机可读存储介质
CN111212223B (zh) * 2020-01-10 2021-01-22 奥比中光科技集团股份有限公司 一种设定成像参数的方法、系统及计算机可读存储介质
WO2021138988A1 (zh) * 2020-01-10 2021-07-15 深圳奥比中光科技有限公司 一种设定成像参数的方法、系统及计算机可读存储介质
CN113545798A (zh) * 2020-04-23 2021-10-26 深圳市理邦精密仪器股份有限公司 一种颅内压测量方法及超声设备
CN114098809A (zh) * 2021-11-29 2022-03-01 无锡海鹰电子医疗系统有限公司 一种超声诊断仪发射系统及其快速优化方法

Also Published As

Publication number Publication date
CN104997533B (zh) 2017-11-14

Similar Documents

Publication Publication Date Title
CN104997533A (zh) 超声探头几何参数自动校正方法和装置
US20180161015A1 (en) Variable speed of sound beamforming based on automatic detection of tissue type in ultrasound imaging
EP2325672A1 (en) Spatial compound imaging in an ultrasound system
US20120289831A9 (en) Ultrasound diagnostic apparatus and method for tracing movement of tissue
US11666310B2 (en) Ultrasound diagnostic apparatus and method for controlling ultrasound diagnostic apparatus using predetermined imaging conditions for B-mode image generation
US9310472B2 (en) Focal point information determination method and apparatus, and ambient sound velocity obtaining method and apparatus
US20220079551A1 (en) Methods and systems for valve regurgitation assessment
KR20110064852A (ko) 부 관심영역에 기초하여 3차원 초음파 영상을 제공하는 초음파 시스템 및 방법
CN106659480A (zh) 超声波诊断图像生成装置以及方法
EP1911402B1 (en) Imaging diagnosis device, measurement point setting method, and program
JP7292370B2 (ja) 胎児体重推定を実施するための方法およびシステム
EP3338643A1 (en) Ultrasonic diagnostic device and method for controlling ultrasonic diagnostic device
CN113573645B (zh) 用于调整超声探头的视场的方法和系统
CN115103634A (zh) 超声波诊断装置、超声波诊断装置的控制方法及超声波诊断装置用处理器
US10792014B2 (en) Ultrasound inspection apparatus, signal processing method for ultrasound inspection apparatus, and recording medium
US9223012B2 (en) Ultrasound image enhancement based on entropy information
CN115330681A (zh) 一种超声图像测量方法、系统、智能终端及存储介质
KR101097539B1 (ko) 볼륨 데이터를 보정하는 초음파 장치 및 방법
CN104684486A (zh) 超声波诊断装置、超声波图像生成方法及程序
US20230360225A1 (en) Systems and methods for medical imaging
US11974884B2 (en) Method and system for dynamically adjusting imaging parameters during an ultrasound scan
EP4159139A1 (en) System and method for segmenting an anatomical structure
US20140121523A1 (en) Measuring apparatus and method thereof
US20210106313A1 (en) Methods and systems for ultrasound image display
CN117392050A (zh) 学习模型、超声波诊断装置及其系统、图像诊断装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20201019

Address after: 430074 East Lake New Technology Development Zone, Wuhan City, Hubei Province

Patentee after: WUHAN YOUKEY BIO-MEDICAL ELECTRONICS Co.,Ltd.

Address before: 430079 Hubei city of Wuhan province East Lake Development Zone Industrial Park Road, Dingxin building two floor

Patentee before: WUHAN QSONO ELECTRONICS Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20231025

Address after: 528000, Unit 601, 6th Floor, Block B, Building 1, Block B1, Hantian Science and Technology City, Dongping Road, Pingxi Shanghai Village, Guicheng Street, Nanhai District, Foshan City, Guangdong Province

Patentee after: Guangdong Qiyou Medical Co.,Ltd.

Address before: 430074 No.1, 2nd floor, building 7, zone B, 818 Gaoxin Avenue, Donghu New Technology Development Zone, Wuhan City, Hubei Province

Patentee before: WUHAN YOUKEY BIO-MEDICAL ELECTRONICS CO.,LTD.