CN103565473A - 图像处理模块和图像处理方法 - Google Patents
图像处理模块和图像处理方法 Download PDFInfo
- Publication number
- CN103565473A CN103565473A CN201310319539.9A CN201310319539A CN103565473A CN 103565473 A CN103565473 A CN 103565473A CN 201310319539 A CN201310319539 A CN 201310319539A CN 103565473 A CN103565473 A CN 103565473A
- Authority
- CN
- China
- Prior art keywords
- spread function
- point spread
- image
- signal
- beam shaping
- 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.)
- Pending
Links
- 238000012545 processing Methods 0.000 title claims abstract description 40
- 238000003672 processing method Methods 0.000 title abstract description 3
- 239000000284 extract Substances 0.000 claims abstract description 8
- 238000007493 shaping process Methods 0.000 claims description 106
- 238000002604 ultrasonography Methods 0.000 claims description 48
- 230000007613 environmental effect Effects 0.000 claims description 32
- 238000000034 method Methods 0.000 claims description 30
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 claims description 5
- 238000003384 imaging method Methods 0.000 abstract description 47
- 239000000523 sample Substances 0.000 abstract description 28
- 230000006870 function Effects 0.000 description 409
- 238000010586 diagram Methods 0.000 description 31
- 238000011084 recovery Methods 0.000 description 13
- 238000013480 data collection Methods 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 5
- 231100000915 pathological change Toxicity 0.000 description 5
- 230000036285 pathological change Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 4
- 239000012141 concentrate Substances 0.000 description 4
- 238000012805 post-processing Methods 0.000 description 4
- 210000001519 tissue Anatomy 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- NBIIXXVUZAFLBC-UHFFFAOYSA-N Phosphoric acid Chemical compound OP(O)(O)=O NBIIXXVUZAFLBC-UHFFFAOYSA-N 0.000 description 2
- 238000000205 computational method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000010409 thin film Substances 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8977—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52046—Techniques for image enhancement involving transmitter or receiver
- G01S7/52049—Techniques for image enhancement involving transmitter or receiver using correction of medium-induced phase aberration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Acoustics & Sound (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Image Processing (AREA)
Abstract
本发明公开一种图像处理模块和图像处理方法。所述超声成像设备包括:超声探头单元,将超声波发射到对象的目标部分,接收被对象的目标部分反射的回波超声波,并将接收的回波超声波转换为超声信号;波束成形单元,对超声信号执行波束成形,并输出波束成形的超声信号;图像生成单元,基于波束成形的超声信号估计至少一个点扩散函数,从而生成图像。图像生成单元可从点扩散函数数据库提取至少一个点扩散函数以决定至少一个点扩散函数,从而生成图像。可选择地,图像生成单元可决定第一点扩散函数,使用第一点扩散函数决定第二点扩散函数,并使用第二点扩散函数生成图像。
Description
本申请要求于2012年7月27日提交到美国专利商标局的第61/676,548号美国专利申请以及于2013年6月3日提交到韩国知识产权局的第10-2013-0063737号韩国专利申请的权益,所述申请的公开通过引用合并于此。
技术领域
实施例涉及一种图像处理模块、图像生成方法和使用所述图像生成方法的超声成像设备。
背景技术
超声成像设备将超声波发射到对象(诸如,人体)中的特定目标部分,收集目标部分反射的超声波,并使用关于收集的超声波的信息获取对象中的各种组织或结构的断层扫描图像,例如,各种内部器官和软组织的断层扫描图像。超声成像设备便宜、小型,不发出辐射(诸如,X射线),并实时再现对象内部的图像。为此,超声成像设备已经在医疗领域中被广泛地使用。
超声成像设备如下获取图像。使用超声成像设备的超声探头的超声波生成单元(诸如,换能器)将电信号转换为超声波。转换的超声波被发射到目标部分。诸如换能器的超声传感器接收目标部分反射的超声波,并将接收的超声波转换为电信号,以获取多通道超声信号。对超声信号执行波束成形。即,由超声传感器收集的超声信号之间的时间差被修正,并且预定权重被添加到各个超声信号以增强来自特定位置的信号并衰减来自另一位置的信号,然后超声信号被集中。随后,超声成像设备基于波束成形的超声信号获取超声图像。
发明内容
本发明的一方面在于提供一种图像处理模块、图像生成方法以及使用所述图像生成方法的超声成像设备,其中,在波束成形之后生成超声图像期间,容易并简单地计算点扩散函数,从而使用少量资源快速地决定点扩散函数。
本发明的另一方面在于提供一种图像处理模块、图像生成方法以及使用所述图像生成方法的超声成像设备,其中,根据点扩散函数的快速决定,更快速地执行由超声成像设备执行的超声图像的生成,此外,在快速地生成图像的同时适当地决定点扩散函数,从而生成高分辨率超声图像。
本发明的另外的方面将在以下描述中被部分地阐明,并且部分地通过描述将是明显的,或者可通过实施本发明被得知。
根据本发明的一方面,一种图像处理模块可包括:波束成形单元,对输入信号执行波束成形,并输出波束成形的信号;点扩散函数数据库,包括基于用于波束成形的信号的至少一个环境变量预先获取的至少一个二维点扩散函数;图像生成单元,从点扩散函数数据库选择至少一个二维点扩散函数,并使用波束成形的信号和选择的点扩散函数执行解卷积,从而生成图像。所述至少一个环境变量可包括超声波的声速或到对象中的目标部分的距离。图像生成单元可基于所述至少一个环境变量对波束成形的信号分类,并基于分类结果从点扩散函数数据库提取至少一个点扩散函数。
根据本发明的另一方面,一种图像处理模块可包括:波束成形单元,对输入信号执行波束成形,并输出波束成形的信号;图像生成单元,基于波束成形的信号估计至少一个第一点扩散函数,基于估计的至少一个第一点扩散函数估计第二点扩散函数,并使用波束成形的信号和估计的第二点扩散函数生成图像。
图像生成单元可使用波束成形的信号和估计的至少一个第一点扩散函数估计第二点扩散函数。图像生成单元可使用第一点扩散函数和至少一个第二点扩散函数估计变量来估计第二点扩散函数。所述至少一个第二点扩散函数估计变量包括从超声波的声速、声速的改变、到对象中的目标部分的距离、输入信号生成模块的排列形式和每个通道的输入信号的衰减率中选择的至少一个。
根据本发明的另一方面,一种超声成像设备可包括:超声探头单元,将超声波发射到对象的目标部分,接收被对象的目标部分反射的回波超声波,并将接收的回波超声波转换为超声信号;波束成形单元,对超声信号执行波束成形,并输出波束成形的超声信号;图像生成单元,基于波束成形的超声信号估计至少一个点扩散函数,从而生成图像。所述超声成像设备还可包括:点扩散函数数据库,包括用于波束成形的超声信号的至少一个点扩散函数。所述超声成像设备可从点扩散函数数据库调用至少一个点扩散函数,以估计第一点扩散函数。所述超声成像设备可使用估计的至少一个第一点扩散函数来估计第二点扩散函数,并使用波束成形的超声信号和估计的第二点扩散函数来生成图像。
根据本发明的另一方面,一种图像生成方法可包括:将超声波发射到对象的目标部分,接收被对象的目标部分反射的回波超声波,并将接收的回波超声波转换为超声信号;输出基于超声信号执行波束成形的超声信号;基于波束成形的超声信号从点扩散函数数据库提取至少一个点扩散函数;使用所述至少一个点扩散函数从波束成形的超声信号生成图像。
根据本发明的另一方面,一种图像生成方法可包括:将超声波发射到对象的目标部分,接收被对象的目标部分反射的回波超声波,并将接收的回波超声波转换为超声信号;输出基于超声信号执行波束成形的超声信号;基于波束成形的超声信号估计至少一个第一点扩散函数;使用估计的至少一个第一点扩散函数估计第二点扩散函数;使用估计的第二点扩散函数从波束成形的信号生成图像。
附图说明
通过以下结合附图对实施例的描述,本发明的这些和/或其他方面将变得清楚和更容易理解,在附图中:
图1是示出图像生成单元的构造的示图;
图2是示出点扩散函数的示图;
图3是示出理想图像和射频(RF)图像与解卷积之间的关系的示图;
图4A是示出理想图像和RF图像之间的关系的示图;
图4B是示出根据深度的基于RF信号的目标部分ob1的超声图像的示例的示图;
图4C是示出用于说明对象的深度的超声图像的示图;
图5是示出图像生成单元的一个实施例的构造的示图;
图6是示出点扩散函数数据库的一个实施例的示图;
图7是示出图像生成单元的另一实施例的构造的示图;
图8和图9是示出根据声速的点扩散函数的示图;
图10是示出图像生成单元的另一实施例的构造的示图;
图11是示出超声成像设备的一个实施例的透视图;
图12是示出超声成像设备的实施例的构造的示图;
图13是示出超声探头单元的一个实施例的平面图;
图14是示出波束成形单元的一个实施例的构造的示图;
图15是示出超声成像设备的图像生成单元的一个实施例的构造的示图;
图16是示出第二点扩散函数的示图;
图17至图19是示出基于输入图像的理想图像恢复方法的各种实施例的流程图;
图20至图22是示出超声成像设备的控制方法的各种实施例的流程图。
具体实施方式
现在将详细参考本发明的实施例,附图中示出本发明的实施例的示例,其中,相同的标号始终表示相同的元件。
图1是示出图像生成单元的构造的示图。
如图1所示,图像生成单元10基于输入信号d生成图像信号m,并输出图像信号m。输入信号d可以是通过诸如声波或电磁波的波生成的信号。例如,在输入信号d是声波信号的情况下,输入信号d可以是从可听频率的声波或大于可听频率(即,20kHz)的频率的超声波获取的信号。另一方面,在输入信号d是电磁波信号的情况下,输入信号d可以是从在雷达中使用的微波(具有10cm至100cm的波长)等获取的信号。
为了从输入信号d生成图像信号m,图像生成单元10估计至少一个点扩散函数(PSF)并使用估计的结果执行解卷积,以生成和获取与将被拍摄的原始对象相同或近似的图像。图像生成单元10以图像信号m的形式输出生成和获取的图像。
点扩散函数是与通过图像拍摄设备的拍摄而获取的用于生成最终图像数据的图像数据相关的函数。点扩散函数主要用于恢复理想图像数据。
图2是示出点扩散函数的示图。如图2所示,在图像拍摄设备的对象图像获取期间,由于图像拍摄设备的技术或物理特性或者噪声η,图像拍摄设备可输出与理想图像o不同的信号,例如,RF信号d(诸如来自超声成像设备的超声信号)。
换句话说,图像拍摄设备获取的RF信号d是具有添加的噪声η并且由于图像拍摄设备的技术或物理特性而变形的理想图像。
图3是示出在超声成像设备中理想图像和RF图像与解卷积之间的关系的示图。图3的最左部分示出人体内组织的理想形状。如果如图3所示,理想图像被给定为fR,则在图3的中间部分由超声成像设备的超声探头收集的且经过波束成形的超声图像被表示为gR。即,理想图像与RF图像不同,将参照图4A至图4C对此进行更加详细的描述。
图4A是示出理想图像和RF图像之间的关系的示图。图4A所示的基于输入信号的图像是使用超声成像设备获取的超声图像的图解示例。如果理想状态下的目标部分ob1的理想图像o被显示为如图4A的左侧部分所示,则基于输入信号d(诸如RF信号)的目标部分ob1的图像被显示为如图4A的右侧部分所示。具体地讲,基于输入信号的图像中的目标部分ob1被显示为好像理想图像中的目标部分ob1垂直和水平地变宽一样。即,基于输入信号的图像与理想图像明显不同。因此,在基于输入信号d(即,RF信号)的图像没有改变地被恢复的情况下,对象看起来与其实际形状不同。
理想图像和基于输入信号的图像可以根据深度而彼此进一步不同。图4B是示出根据深度的基于RF信号的目标部分ob1的超声图像的示例的示图,图4C是示出用于说明对象的深度的超声图像的示图。
如果如图4B所示目标部分ob1和图像数据收集单元(诸如超声探头单元)之间的距离较短,例如,人体内的病变位于如图4C所示的第一深度,则基于输入信号的目标部分ob1的图像与相同对象的理想图像相同或非常相似。然而,如果目标部分ob1和图像数据收集单元之间的距离增加,例如,人体内的病变位于如图4C所示的第四深度或第五深度,则基于输入信号的目标部分ob1的图像沿横向扩展,因此与同样的目标部分ob1的理想图像明显不同。即,理想图像o中的目标部分ob1与基于输入信号的图像中的目标部分ob1根据数据收集单元与目标部分ob1之间的距离在形状上彼此进一步不同。
因此,在使用RF信号d恢复理想图像o的情况下,可修正理想图像o和基于RF信号d的图像之间的不同,以获得目标部分ob1的正确图像。在这种情况下,假设在理想图像o和获取的RF信号d之间实现了预定关系,则使用预定函数修正RF信号d以恢复图像。这时使用的函数是点扩散函数H。
理想图像o、点扩散函数H、噪声η和输入信号d(即,RF信号)之间的关系可通过等式1来表示。
[等式1]
d=Hf+η
其中,d指示输出的RF信号,H指示点扩散函数,f指示理想图像的信号,η指示噪声。
假设不存在噪声,则RF信号d可通过理想图像o和点扩散函数H的乘积来表示。因此,如果适用于测量的RF信号d的点扩散函数H已知,则可获得与测量的RF信号d相应的理想图像o。换句话说,如果点扩散函数H和RF信号d已知,则可恢复与对象相同或近似的图像。
图像生成单元10使用RF信号(即,输入信号d)和适合于输入信号d的点扩散函数H来生成与理想图像o相同或非常相似的图像,并沿与图2所示的方向相反的方向输出与生成的图像相应的图像信号m。换句话说,图像生成单元10通过如图3的中间部分所示的解卷积添加并合成适用于RF信号gR的点扩散函数hR,以生成与理想图像o(fR)相同或近似的恢复的图像。
为了基于输入信号d恢复与对象相同或近似的图像,图像生成单元10使用包括一维或二维点扩散函数的点扩散函数(PSF)数据库12或使用基于一维点扩散函数估计二维点扩散函数的方法,在沿各个方向不降低分辨率的情况下估计点扩散函数。
图像生成单元10使用估计的点扩散函数H转换通过拍摄对象所获取的输入信号d(诸如RF信号),从而使基于输入信号d的图像具有与原始对象相同或相似的形式或形状。即,图像生成单元10使用点扩散函数恢复获取的输入信号d的图像,并输出恢复的图像的图像信号m。
图5是示出图像生成单元的一个实施例的构造的示图。
如图5所示,图像生成单元10可包括点扩散函数选择单元11和解卷积单元14。
点扩散函数选择单元11检索点扩散函数数据库12,以从点扩散函数数据库12选择和调用至少一个点扩散函数。图6是示出点扩散函数数据库12的一个实施例的示图。如图6所示,点扩散函数数据库12包括至少一组点扩散函数。
在一个实施例中,点扩散函数数据库12可仅包括一维点扩散函数或仅包括二维点扩散函数。可选择地,点扩散函数数据库12可包括一维点扩散函数和二维点扩散函数。在其他实施例中,点扩散函数数据库12可包括更高维数的点扩散函数,诸如三维点扩散函数或四维点扩散函数。
存储在点扩散函数数据库12中的点扩散函数可以是基于至少一个环境变量预先计算的点扩散函数。在一个实施例中,环境变量可以是波的速度(图6的水平方向)。波的速度可以是声波(诸如超声波)的速度(声速)。在另一实施例中,环境变量可以是数据收集单元(诸如超声探头单元)和目标部分之间的距离(图6的垂直方向)。例如,数据收集单元和目标部分之间的距离可以是人体内的特定目标部分(诸如病变)的深度。
当然,存储在点扩散函数数据库12中的点扩散函数可以是基于两个环境变量预先计算的点扩散函数。在这种情况下,两个环境变量可以是如图6所示的声速和目标部分的深度。换句话说,如图6所示,点扩散函数数据库12可包括根据环境变量(诸如声速和深度)分别获取的多组点扩散函数。
假设轴向上的一维点扩散函数彼此相似,则为了快速估计二维点扩散函数,点扩散函数数据库12可包括通过根据各种环境变量(诸如超声波的速度和目标部分ob1的深度)的实时测量而获取的至少一个二维点扩散函数,下文将对此进行描述。
点扩散函数选择单元11从如上所述构造的点扩散函数数据库12选择至少一个点扩散函数。在一个实施例中,点扩散函数选择单元11可从点扩散函数数据库12选择二维点扩散函数。可选择地,点扩散函数选择单元11可基于各种环境变量对波束成形的信号分类,并基于分类结果从点扩散函数数据库12选择二维点扩散函数。
在这种情况下,点扩散函数选择单元11可根据输入信号d选择至少一个点扩散函数。如果目标部分ob1位于人体内不是很深的位置,例如图4C的第一深度,则点扩散函数选择单元11从点扩散函数数据库12选择与第一深度相应的点扩散函数。在拍摄位于图4C的第四深度或第五深度的目标部分ob1的情况下,点扩散函数选择单元11从点扩散函数数据库12选择与第四深度或第五深度相应的点扩散函数。
点扩散函数选择单元11将选择的点扩散函数或关于选择的点扩散函数的信息发送到解卷积单元14。
解卷积单元14使用点扩散函数选择单元11选择的点扩散函数对输入信号d进行解卷积。因此,从输入信号d(即,图2所示的RF信号d)获取与理想图像o相同或近似的图像信号m。
在点扩散函数选择单元11从点扩散函数数据库12选择二维点扩散函数的情况下,解卷积单元14基于选择的二维点扩散函数对输入信号d进行二维解卷积。在点扩散函数选择单元11从点扩散函数数据库12选择一维点扩散函数的情况下,解卷积单元14可对输入信号d进行一维解卷积。
因此,图像生成单元10可从输入信号d生成与理想图像o相同或相似的图像。生成的图像作为图像信号m被输出。
如上所述,图像生成单元10可根据输入信号d从点扩散函数数据库12选择至少一个二维点扩散函数,并使用选择的二维点扩散函数对输入信号d进行解卷积,以基于输入信号d生成与理想图像o相同或相似的图像。
如图2以及等式1所示,执行从RF信号d获取理想图像o的处理以计算点扩散函数的反函数。在点扩散函数是二维的情况下,当直接计算二维点扩散函数的反函数时,计算复杂度增加。
在使用一维点扩散函数的情况下,计算量减小。然而,由于一维点扩散函数的限制,特定区域的分辨率降低。在一个方向上,例如,在垂直方向上,分辨率没有降低。在另一个方向上,例如,在横向上,分辨率会降低。
然而,图像生成单元10从点扩散函数数据库12选择适用于输入信号d的二维点扩散函数,并使用选择的点扩散函数执行二维解卷积。结果,可避免获取反函数的复杂的计算处理,从而简化恢复图像的计算并提高分辨率。
图7是示出图像生成单元的另一实施例的构造的示图。
如图7所示,图像生成单元10可包括第一点扩散函数选择单元11a、第二点扩散函数估计单元13和解卷积单元14。
以与如上所述的点扩散函数选择单元11相似的方式,第一点扩散函数选择单元11a可从附加的点扩散函数数据库12选择至少一个点扩散函数,并将选择的点扩散函数发送到第二点扩散函数估计单元13。
在一个实施例中,点扩散函数数据库12可包括基于至少一个环境变量预先计算的一维点扩散函数。例如,如图4C所示,点扩散函数数据库12可包括根据环境变量(诸如声速和深度)分别获取的点扩散函数。
在一个实施例中,第一点扩散函数选择单元11a可从点扩散函数数据库12选择至少一个一维点扩散函数。
在这种情况下,假设轴向上的一维点扩散函数彼此相似,则为了快速估计二维点扩散函数,第一点扩散函数选择单元11a可根据各种环境变量(诸如超声波的速度和目标部分ob1的深度)选择横向上的一维点扩散函数。第一点扩散函数选择单元11a选择的一维点扩散函数被发送到第二点扩散函数估计单元13。
第二点扩散函数估计单元13使用接收的一维点扩散函数估计二维点扩散函数。此时,发送到第二点扩散函数估计单元13的点扩散函数可以是横向上的一维点扩散函数。
一维点扩散函数可以按轴向或横向被划分。轴向上的一维点扩散函数基于深度或声速彼此近似。
图8和图9是示出根据声速和深度的点扩散函数的示图。图8和图9的左侧部分以曲线图的形式示出在超声成像设备中使用的一维点扩散函数(1DPSF)的示例。图8是示出在超声波具有1540m/s的声速的情况下通过将超声波发射到具有不同深度的三个目标部分所获取的超声波的一维点扩散函数的曲线图。图9是示出在超声波具有1450m/s的声速的情况下通过将超声波发射到具有不同深度的三个目标部分所获取的超声波的一维点扩散函数的曲线图。
如图8和图9的左侧部分所示,虽然超声波的声速彼此不同,但是一维点扩散函数在形状上彼此非常相似。此外,即使在目标部分的深度彼此不同(深度1至深度3)的情况下,一维点扩散函数也具有非常相似的模式。
因此,虽然超声波的声速彼此不同并且目标部分的深度彼此不同,但是一维点扩散函数在形状上彼此非常相似。
图8和图9的右侧部分示出基于深度的二维点扩散函数。与图8和图9的左侧部分所示的一维点扩散函数不同,二维点扩散函数根据深度彼此非常不同。即,与轴向上的点扩散函数不同,横向上的点扩散函数根据深度或声速具有不同的形状。
假设轴向上的一维点扩散函数是一致的,则因为轴向上的一维点扩散函数根据超声波的声速或深度没有明显的改变,所以第二点扩散函数估计单元13可基于第一点扩散函数选择单元11a选择的横向上的一维点扩散函数,使用计算方法来估计二维点扩散函数。
在一个实施例中,第二点扩散函数估计单元13可考虑输入信号d的特性、预先输入的设置或各种变量(诸如第二点扩散函数估计变量)、以及第一点扩散函数选择单元11a选择的横向上的一维点扩散函数,来估计二维点扩散函数。
例如,在第二点扩散函数估计单元13中使用的第二点扩散函数估计变量可以是用于对输入信号d执行波束成形的波束成形系数。
在第二点扩散函数估计变量被应用于超声成像设备的情况下,第二点扩散函数估计变量可以是诸如超声波的声速、声速的改变、到对象中的目标部分的距离、输入信号生成模块的排列形式或每个通道的输入信号的衰减率的变量。第二点扩散函数估计单元13可使用这种变量来估计更正确的二维点扩散函数。
解卷积单元14使用估计的二维点扩散函数执行解卷积,以生成和输出图像信号m。
因此,图像生成单元10可根据输入信号d从点扩散函数数据库12选择至少一个一维点扩散函数,基于选择的一维点扩散函数估计二维点扩散函数,并根据估计的二维点扩散函数执行解卷积,以基于输入信号d生成与理想图像o相同或相似的图像。
图10是示出图像生成单元的另一实施例的构造的示图。
如图10所示,图像生成单元10可包括第一点扩散函数估计单元11b、第二点扩散函数估计单元13和解卷积单元14。
第一点扩散函数估计单元11b可基于输入信号d估计适当的一维点扩散函数。在估计一维点扩散函数的情况下,与估计二维点扩散函数的情况相比,考虑更少的输入信号d的变量。因此,与直接估计二维点扩散函数相比,决定点扩散函数所需的计算量减小。
在一个实施例中,第一点扩散函数估计单元11b可估计一个方向或多个方向上的一个点扩散函数。可选择地,第一点扩散函数估计单元11b可估计一个方向或多个方向上的多个点扩散函数。
在一个实施例中,第一点扩散函数估计单元11b可使用至少一个滤波器来估计第一点扩散函数。在这种情况下,滤波器可以是最小二乘滤波器(LSF),诸如最小的最小二乘滤波器或倒频谱滤波器。
在一个实施例中,第一点扩散函数估计单元11b可估计横向上的点扩散函数。
第二点扩散函数估计单元13使用第一点扩散函数估计单元11b估计的第一点扩散函数来估计第二点扩散函数。
在一个实施例中,假设如上参照图8和图9所述,轴向上的点扩散函数根据声速和深度几乎一致,则第二点扩散函数估计单元13可基于第一点扩散函数估计单元11b估计的横向上的点扩散函数,使用计算方法来估计二维点扩散函数。
如图10所示,第二点扩散函数估计单元13可考虑输入信号d的特性、预先输入的设置或各种变量、以及第一点扩散函数估计单元11b估计的横向上的一维点扩散函数,来估计第二点扩散函数。
以与上述相同的方式,第二点扩散函数估计单元13可使用波束成形系数来估计第二点扩散函数。此外,在第二点扩散函数估计单元13用于超声成像设备中的情况下,第二点扩散函数估计单元13可使用另一变量(诸如超声波的声速或到对象中的目标部分的距离)来估计第二点扩散函数。
解卷积单元14使用估计的二维点扩散函数执行解卷积,以生成和输出图像信号m。
如上所述,图像生成单元10可在不检索附加的数据库的情况下,根据输入信号d来估计至少一个一维点扩散函数,基于估计的一维点扩散函数估计二维点扩散函数,并使用估计的二维点扩散函数来对输入信号d进行解卷积,以基于输入信号d恢复与理想图像o相同或相似的图像。
以下,将参照图11至图18描述超声成像设备的实施例。
图11是示出超声成像设备的一个实施例的透视图,图12是示出超声成像设备的实施例的构造的示图。如图11和图12所示,超声成像设备包括超声探头单元P和主体M,超声探头单元P收集关于对象ob的内部的超声信号,主体M从超声探头单元P收集的超声信号生成超声图像。
超声探头单元P可包括超声生成器P11和超声接收器P12,超声生成器P11基于施加的功率生成超声波,并将生成的超声波发射到对象ob中的目标部分ob1,超声接收器P12接收对象ob的目标部分ob1反射的回波超声波,以将接收的回波超声波转换为电信号。
超声生成器P11在安装在主体M中的超声生成控制器310的控制下,根据施加到超声生成器P11的脉冲信号或AC电流而振动,以生成超声波。
超声生成器P11生成的超声波被对象中的目标部分ob1反射。超声接收器P12接收回波超声波,并根据接收的回波超声波的频率而振动,以将接收的回波超声波转换为预定的电信号。结果,超声接收器P12可输出超声信号x。在超声成像设备连接到光声成像设备的情况下,超声接收器P12可接收由于发射激光而从目标部分ob1生成的诸如超声波的声波。
超声生成器P11和超声接收器P12的功能可由设置在超声探头单元P的一端的超声换能器P10执行。图13是示出超声探头单元的一个实施例的平面图。如图13所示,超声换能器P10安装在超声探头单元P的一端。
换能器是用于将一种形式的能量(诸如电能)转换为另一种形式的能量(诸如波能或光能)的装置。超声换能器P10执行波能和电能之间的转换。具体地讲,超声换能器P10根据输入的预定脉冲电流而振动,以生成超声波。此外,超声换能器P10根据从外部接收的超声波(诸如回波超声波)而振动,以生成预定脉冲的电信号(以下,称为超声信号)。因此,超声换能器P10可执行超声生成器P11和超声接收器P12的功能。
更具体地讲,超声换能器P10从电源单元311(诸如外部电源装置或内部电存储装置,例如,电池)接收AC电流,并且超声换能器P10的压电振动器或薄膜根据施加的AC电流而振动,以生成超声波。另一方面,超声换能器P10的压电振动器或薄膜根据接收的超声波而振动,以生成振动频率的频率的AC电流,从而将超声波转换为超声信号。如图12和图13所示,转换的超声信号通过多个通道c1至c10被发送到主体M的波束成形单元100。
如图13所示,多个超声换能器P10可安装在超声探头单元P的一端。例如,64或128个超声换能器P10可安装在超声探头单元P的一端。在多个超声换能器P10安装在超声探头单元P的一端的情况下,超声信号通过与超声换能器P10的数量(例如,64或128个超声换能器P10)相应的多个通道被发送到波束成形单元100。
超声换能器P10的示例可包括使用磁体的磁致伸缩效应的磁致伸缩超声换能器、使用压电材料的压电效应的压电超声换能器、以及使用成百上千个微机械薄膜的振动来发送和接收超声波的电容式微机械超声换能器。此外,超声换能器P10的示例可包括用于根据电信号生成超声波或根据超声波生成电信号的其他种类的换能器。
在一个实施例中,主体M还可包括波束成形单元100、图像生成单元200、点扩散函数数据库220、系统控制器300和图像处理单元320。
主体M的波束成形单元100从超声接收器P12接收超声信号,并对接收的超声信号执行波束成形。具体地讲,波束成形单元100对多个通道数据执行波束成形,以估计在特定空间内反射的波的大小。
图14是示出波束成形单元的一个实施例的构造的示图。如图14所示,波束成形单元100可包括时间差修正部分110和集中部分120。
参照图14,超声换能器P10(具体地讲,超声接收器P12)接收目标部分ob1反射的回波超声波。然而,安装在超声探头单元P中的超声换能器T1至T6与目标部分ob1之间的距离不同,而声速几乎一致。因此,即使在回波超声波被同一目标部分ob1反射的情况下,各个超声换能器T1至T6也在不同时间从同一目标部分ob1接收回波超声波。换句话说,各个超声换能器T1至T6输出的超声信号基于由同一目标部分ob1反射的回波超声波而具有预定时间差。为此,可首先修正各个超声换能器T1至T6生成的超声信号之间的时间差。
波束成形单元100的时间差修正部分110修正超声信号之间的时间差。例如,如图13和图14所示,时间差修正部分110将输入到特定通道的超声信号的发送延迟预定程度,以使输入到各个通道的超声信号同时被发送到集中部分120。
波束成形单元100的集中部分120集中时间差已经被修正的超声信号。
集中部分120将预定权重(即,波束成形系数)添加到每个输入的超声信号,以增强来自特定位置的信号并衰减来自另一位置的信号,然后集中超声信号。结果,可生成基于用户需求或用户便利性的超声图像。
在这种情况下,集中部分120可使用不考虑超声接收器P12输出的超声信号的波束成形系数来集中超声信号(数据独立型波束成形方法)。可选择地,可基于输入的超声信号计算最优的波束成形系数,并可使用计算的波束成形系数来集中超声信号(自适应波束成形方法)。
可通过等式2来表示这种波束成形处理。
[等式2]
其中,n是指示目标部分ob1的位置的值(即,目标部分ob1的特定深度),m指示在目标部分ob1的位置n的第m通道,wm指示添加到第m通道的超声信号的波束成形系数w。同时,Δm是用于时间差修正部分110延迟超声信号的发送时间的时间延迟值。根据等式2,集中部分120集中各个通道的时间差已经被修正的超声信号,以输出波束成形的超声信号d。
如图12和图14所示,通过波束成形单元100执行波束成形的超声信号d被发送到图像生成单元200。
图15是示出超声成像设备的图像生成单元的一个实施例的构造的示图。
如图15所示,图像生成单元200可包括点扩散函数选择单元210和解卷积单元240。
点扩散函数选择单元210从点扩散函数数据库220选择至少一个点扩散函数。在这种情况下,点扩散函数数据库220可包括一维点扩散函数和/或二维点扩散函数。
如上所述,当在超声成像设备中恢复超声图像时对点扩散函数影响明显的因素包括超声波的速度(即,声速)以及产生或反射超声波的目标部分与收集单元之间的距离(例如,人体内的病变的深度)。因此,点扩散函数数据库220可包括使用对点扩散函数影响明显的环境变量(诸如超声波的速度和目标位置的深度)预先计算和测量的一维点扩散函数或二维点扩散函数,以恢复超声波。
在一个实施例中,在点扩散函数选择单元210选择一维点扩散函数的情况下,图像生成单元200还可包括第二点扩散函数估计单元230,第二点扩散函数估计单元230基于点扩散函数选择单元210选择的一维点扩散函数估计二维点扩散函数。
在这种情况下,点扩散函数选择单元210基于波束成形的信号d,根据对象ob内的目标部分ob1的深度或超声波的速度,从点扩散函数数据库220选择适当的一维点扩散函数。在这种情况下,一维点扩散函数可以是目标部分ob1的横向上的点扩散函数。
第二点扩散函数估计单元230基于点扩散函数选择单元210选择的一维点扩散函数估计二维点扩散函数。假设轴向上的一维点扩散函数是一致的,则因为如上所述,轴向上的点扩散函数根据声速或目标部分ob1的深度没有明显的改变,所以第二点扩散函数估计单元230可使用点扩散函数选择单元210选择的一维点扩散函数(例如,横向上的点扩散函数)估计并获取二维点扩散函数。
图16是示出第二点扩散函数的示图。第二点扩散函数估计单元230可如图16的左侧部分所示组合选择的一维点扩散函数和附加的系数(诸如第二点扩散函数估计系数),以估计如图16的右侧部分所示的第二点扩散函数。
如图16所示,第二点扩散函数估计系数可以是波束成形系数。可使用从超声换能器P10的排列结构、声速、深度和信号衰减率中选择的至少一个来计算波束成形系数。另一方面,第二点扩散函数估计系数可以是从超声换能器P10的排列结构、声速、深度和信号衰减率中选择的至少一个。为此,第二点扩散函数估计单元230可从波束成形单元100接收波束成形的超声信号d或关于波束成形的超声信号d的信息。
第二点扩散函数估计单元230可将如上所述估计的二维点扩散函数发送到解卷积单元240,解卷积单元240可对波束成形的信号d进行二维解卷积。
在图像生成单元200的另一实施例中,点扩散函数选择单元210可从点扩散函数数据库220选择与输入的波束成形的超声信号d相应的适当的二维点扩散函数。
在这种情况下,假设轴向上的一维点扩散函数相似,则点扩散函数数据库220的至少一个二维点扩散函数可以是根据各种环境变量实际测量并获取的点扩散函数。
在选择和决定二维点扩散函数之后,点扩散函数选择单元210可将决定的点扩散函数直接发送到解卷积单元240,解卷积单元240可对波束成形的信号d进行二维解卷积。
在另一实施例中,图像生成单元200可包括第一点扩散函数估计单元211和第二点扩散函数估计单元230。第一点扩散函数估计单元211在不检索点扩散函数数据库220的情况下基于波束成形的超声信号d估计一维点扩散函数。
在这种情况下,第一点扩散函数估计单元211可使用至少一个滤波器(诸如最小的最小二乘滤波器或倒频谱滤波器)来估计一维点扩散函数。
估计的一维点扩散函数被发送到第二点扩散函数估计单元230。第二点扩散函数估计单元230基于估计的一维点扩散函数估计二维点扩散函数。以如上所述的相同的方式,假设轴向上的一维点扩散函数是一致的,则第二点扩散函数估计单元230可使用第一点扩散函数估计单元211估计的一维点扩散函数来估计二维点扩散函数。
第二点扩散函数估计单元230将估计的二维点扩散函数发送到解卷积单元240。
解卷积单元240使用从点扩散函数选择单元210或第二点扩散函数估计单元230接收的二维点扩散函数对波束成形的信号d进行二维解卷积,以恢复与理想图像o相同或相似的图像。
恢复的图像的图像信号可被发送到图像处理单元320、存储单元321或显示单元dp。
图像处理单元320可基于图像生成单元200生成的图像来生成可由用户识别的超声图像,或者对图像生成单元200生成的图像执行预定的图像处理。
在一个实施例中,图像处理单元320可基于从图像生成单元200输出的图像信号m以各种模式生成超声图像。A模式和B模式可用作超声图像模式。A模式是使用幅值显示超声图像的模式。具体地讲,A模式是这样的模式:将目标部分t显示为到超声探头单元10的距离或时间,并将反射强度显示为幅值。B模式是在屏幕上将回波超声波的大小显示为亮度的模式。在以B模式生成超声图像的情况下,用户可仅使用超声图像容易直观地检查对象ob内的组织或结构。为此,通常使用B模式。图4C中示出了B模式下的超声图像。
此外,图像处理单元320可根据用户意图或为用户方便而修正图像生成单元200生成的超声图像。例如,图像处理单元320可修正亮度、照度、对比度或颜色,从而用户可清楚地观看超声图像中的组织。在从图像生成单元200输出多个超声图像的情况下,可使用输出的超声图像生成三维立体超声图像。
存储单元321临时或持续存储图像生成单元200生成的图像信号以及图像处理单元320修正的图像的图像信号。
显示单元dp将超声图像呈现给用户。根据实施例,显示单元dp可将图像生成单元200解卷积的超声图像直接显示给用户,或者将通过图像处理单元320对图像生成单元200生成的图像执行预定图像处理而获取的超声图像显示给用户。在这种情况下,显示在显示单元dp上的超声图像可以是B模式下的超声图像或三维立体超声图像。
同时,超声成像设备的主体M还可包括超声生成控制器310。超声生成控制器310从系统控制器300接收控制命令,根据接收的控制命令生成预定的脉冲信号,并将生成的脉冲信号发送到超声生成器P11。超声生成器P11根据接收的脉冲信号而振动,以生成超声波。此外,超声生成控制器310可生成用于电连接到超声生成器P11的电源单元311的附加控制信号,并将生成的控制信号发送到电源单元311。电源单元311根据接收的控制信号将预定的AC电流施加到超声生成器P11。
系统控制器300控制超声探头单元P、波束成形单元100、图像生成单元200、图像处理单元320、存储单元321和显示单元dp的操作。
根据实施例,系统控制器300可根据预定设置或通过附加的输入单元i输入的用户指令或命令生成用于超声探头单元P等的预定控制命令。
输入单元i接收用于控制超声成像设备的用户指令或命令。输入单元i可包括各种用户接口,诸如键盘、鼠标、轨迹球、触摸屏和开关。可选择地,输入单元i可以是连接到主体M的工作站。
在如上所述的超声成像设备的实施例中,超声换能器P10安装在超声探头单元P中,波束成形单元100、图像生成单元200、超声生成控制器310和图像处理单元320安装在主体M中。然而,根据实施例,安装在主体M中一些组件可设置在超声探头单元P中。
以下,将参照图17至图19描述基于输入图像的理想图像恢复方法的各种实施例。
图17是示出使用点扩散函数的图像生成方法的一个实施例的流程图。如图17所示,首先输入信号d。在一个实施例中,输入信号d可以是超声成像设备的波束成形的信号,诸如波束成形的超声信号(s410)。
随后,分析输入信号d(s411)。在这种情况下,另外添加到输入信号d的元信息可用于分析输入信号d。可基于诸如用于获取输入信号d的超声波的声速或数据收集单元和目标部分ob1之间的距离(例如,目标部分ob1的深度)的环境变量来分析输入信号d。
根据输入信号d的分析结果选择适当的二维点扩散函数。选择的二维点扩散函数可以是从存储在点扩散函数数据库12中的多个点扩散函数中选择的至少一个。
从点扩散函数数据库12调用选择的二维点扩散函数(s412)。存储在点扩散函数数据库12中的点扩散函数可以是根据环境变量(诸如如图6所示的超声波的声速或数据收集单元和目标部分ob1之间的距离)预先测量的点扩散函数。
使用选择的二维点扩散函数对输入信号d执行二维解卷积(s413)。结果,获取并输出与理想图像相同或非常相似的图像信号m(s414)。
图18是示出使用点扩散函数的图像生成方法的另一实施例的流程图。参照图18,首先输入信号d(s420)。以与如上所述的相同方式分析输入信号d(s421)。输入信号d可以是波束成形的信号。此外,可基于如上所述的环境变量分析输入信号d。例如,环境变量可以是声速或深度。
根据输入信号d的分析结果选择适当的一维点扩散函数。选择的一维点扩散函数可以是存储在点扩散函数数据库12中的至少一个点扩散函数。在一个实施例中,存储在点扩散函数数据库12中的点扩散函数可以是根据环境变量(诸如超声波的声速或数据收集单元和目标部分ob1之间的距离)预先测量的一维点扩散函数。例如,一维点扩散函数可以是横向上的点扩散函数。从点扩散函数数据库12调用选择的一维点扩散函数(s422)。
使用调用的一维点扩散函数来估计二维点扩散函数。因为如上所述,轴向上的点扩散函数根据声速或距离没有明显的改变,所以可基于横向上的一维点扩散函数估计二维点扩散函数(s423)。
同时,根据实施例,可使用预先定义的用于估计二维点扩散函数的第二点扩散函数估计变量或者用户选择的第二点扩散函数估计变量。例如,第二点扩散函数估计变量可以是波束成形系数。更具体地讲,第二点扩散函数估计变量可以是诸如超声波的声速、声速的改变、到对象中的目标部分的距离、输入信号生成模块的排列形式或每个通道的输入信号的衰减率的变量。
使用估计的二维点扩散函数对输入信号d执行解卷积(s424),以恢复图像。其后,获取恢复的图像的图像信号(s425)。
图19是示出使用点扩散函数的图像生成方法的另一实施例的流程图。如图19所示,首先输入信号d(s430)。如上所述,输入信号d可以是波束成形的信号。
根据输入信号d的分析结果估计一维点扩散函数(s431)。此时,可使用至少一个滤波器针对输入信号d估计一维点扩散函数。可使用最小的最小二乘滤波器或倒频谱滤波器。
当如上所述估计了一维点扩散函数时,基于估计的一维点扩散函数估计二维点扩散函数(s432)。在这种情况下,可以以与参照图18描述的实施例中相似的方式估计二维点扩散函数。
当获取了二维点扩散函数时,使用获取的二维点扩散函数对输入信号d执行二维解卷积(s433)。结果,恢复图像并获取恢复的图像的图像信号m(s434)。
以下,将参照图20至图22描述超声成像设备的控制方法的各种实施例。
图20是示出超声成像设备的控制方法的一个实施例的流程图。如图20所示,超声成像设备的超声探头单元P将超声波发射到对象的目标部分ob1(s440)。目标部分ob1完全或部分地反射发射的超声波。超声探头单元P接收反射的超声波(即,回波超声波),并将接收的超声波转换为电信号,以输出超声信号(s441)。
修正输出的超声信号之间的时间差,然后集中超声信号以对超声信号执行波束成形(s442)。分析波束成形的超声信号(s443),并选择适合于波束成形的超声信号d的点扩散函数。
在这种情况下,诸如超声波的声速或目标部分ob1的深度的环境变量可用于分析波束成形的超声信号d。为了分析波束成形的超声信号d,可使用添加到超声信号d的各种信息以及波束成形的超声信号d。作为分析的结果,可决定用于获取超声信号的超声波的声速和目标部分ob1的深度。
根据分析结果,选择适当的点扩散函数。选择的点扩散函数可存储在点扩散函数数据库220中。例如,如图6所示,当根据分析决定了用于获取超声信号的超声波的声速和目标部分ob1的深度时,从点扩散函数数据库220选择与决定的声速和深度相应的点扩散函数。从点扩散函数数据库220调用点扩散函数(s444)。
使用调用的点扩散函数对波束成形的超声信号d进行二维解卷积(s445)。结果,恢复并获取超声图像(s446)。对恢复和获取的超声图像执行图像后处理,以按照需要调整图像的亮度或对比度(s447),然后将经图像处理后的超声图像输出到诸如监视器的显示装置dp(s448)。
图21是示出超声成像设备的控制方法的另一实施例的流程图。如图21所示,超声探头单元P将超声波发射到对象的目标部分ob1(s450),接收被目标部分ob1反射的回波超声波,并将接收的超声波转换为超声信号(s451)。随后,对超声信号执行诸如时间差修正和集中的波束成形(s452)。
分析波束成形的超声信号(s453),并根据分析结果选择适当的一维点扩散函数。一维点扩散函数可以是存储在点扩散函数数据库220中的点扩散函数。此外,诸如对象ob中的目标部分ob1的深度(例如,人体内病变的深度)的关于目标部分ob1的位置的信息可用于选择适当的一维点扩散函数。此外,还可使用关于超声探头单元生成的超声波的信息或关于接收的回波超声波的声速的信息。同时,选择的点扩散函数可以是横向上的点扩散函数(s454)。
基于一维点扩散函数估计二维点扩散函数(s455)。假设用于实际恢复图像的轴向上的点扩散函数与预定的轴向上的点扩散函数相同,则因为如上所述,轴向上的点扩散函数在每个声速彼此近似,所以使用选择的横向上的点扩散函数来计算和估计二维点扩散函数。
在一个实施例中,可使用由用户预先定义或选择的二维点扩散函数估计变量(例如,波束成形系数)来估计二维点扩散函数。超声波的声速、声速的改变、到对象中的目标部分的距离、输入信号生成模块的排列形式或每个通道的输入信号的衰减率可用作二维点扩散函数估计变量。
使用二维点扩散函数对波束成形的超声信号d进行二维解卷积(s456),以恢复超声图像(s457)。恢复的超声图像通过预定的图像后处理被修正(s458),并被输出到显示装置(s459)。
图22是示出超声成像设备的控制方法的另一实施例的流程图。如图22所示,超声探头单元P将超声波发射到对象的目标部分ob1(s460),接收被目标部分ob1反射的回波超声波,并将接收的超声波转换为超声信号(s461)。随后,以与前述实施例中相同的方式对超声信号执行诸如时间差修正和集中的波束成形(s462)。
分析波束成形的超声信号d(s463),以估计用于波束成形的超声信号d的一维点扩散函数(s464)。可使用至少一个滤波器(诸如最小二乘滤波器或倒频谱滤波器)来估计一维点扩散函数。可选择地,各种传统方法可用于估计一维点扩散函数。此时,可以不估计轴向上的一维点扩散函数,而可估计横向上的一维点扩散函数。
基于估计的一维点扩散函数估计二维点扩散函数(s465)。以与前述实施例中相同的方式,假设用于实际恢复图像的轴向上的点扩散函数与预定的轴向上的点扩散函数相同,则使用估计的横向上的点扩散函数来计算和估计二维点扩散函数。在这种情况下,如上所述的二维点扩散函数估计变量可用于估计二维点扩散函数。
使用二维点扩散函数对波束成形的超声信号d进行二维解卷积(s466)。结果,恢复与理想图像o相同或非常相似的超声图像(s467)。可通过预定的图像后处理来修正恢复的超声图像(s468)。将恢复的或后处理的超声图像输出到显示装置(s469)。
从以上描述清楚的是,在图像处理模块、使用该图像处理模块的超声成像设备以及使用该超声成像设备的超声图像生成方法中,当恢复图像(例如,基于超声信号生成图像)时,可决定适当的点扩散函数。
在由超声成像设备执行的超声图像的生成期间可快速决定点扩散函数,从而提高超声图像的恢复速度并因此快速生成超声图像。
除了快速生成超声图像外,还可决定适当的点扩散函数,从而生成高分辨率的超声图像。
可基于一维点扩散函数高速获取二维点扩散函数,从而容易地获取高分辨率和高质量的超声图像。
因此,诸如医生或超声检查者的用户可使用超声成像设备快速准确地获取的超声图像来更准确地诊断病人。
此外,决定点扩散函数的方法可应用于雷达或声学信号处理,从而根据实施例改善雷达或声学信号处理设备的操作。
虽然已经示出和描述了本发明的一些实施例,但是本领域技术人员将理解,在不脱离本发明的原理和精神的情况下,可对这些实施例进行改变,本发明的范围由权利要求及其等同物限定。
Claims (20)
1.一种图像处理模块,包括:
波束成形单元,对输入信号执行波束成形,并输出波束成形的信号;
点扩散函数(PSF)数据库,包括基于用于波束成形的信号的至少一个环境变量预先获取的至少一个二维点扩散函数;
图像生成单元,从点扩散函数数据库选择至少一个二维点扩散函数,并使用波束成形的信号和选择的点扩散函数执行解卷积,从而生成图像。
2.根据权利要求1所述的图像处理模块,其中,所述至少一个环境变量包括超声波的声速或到对象中的目标部分的距离。
3.根据权利要求1所述的图像处理模块,其中,图像生成单元基于所述至少一个环境变量对波束成形的信号分类,并基于分类结果从点扩散函数数据库选择至少一个二维点扩散函数。
4.一种图像处理模块,包括:
波束成形单元,对输入信号执行波束成形,并输出波束成形的信号;
图像生成单元,基于波束成形的信号估计至少一个第一点扩散函数,基于估计的至少一个第一点扩散函数估计第二点扩散函数,并使用波束成形的信号和估计的第二点扩散函数执行解卷积,从而生成图像。
5.根据权利要求4所述的图像处理模块,还包括:
点扩散函数数据库,包括用于波束成形的信号的至少一个点扩散函数,
其中,基于至少一个环境变量预先获取所述至少一个点扩散函数。
6.根据权利要求5所述的图像处理模块,其中,所述至少一个环境变量包括超声波的声速或到对象中的目标部分的距离。
7.根据权利要求5所述的图像处理模块,其中,图像生成单元从点扩散函数数据库提取至少一个点扩散函数以估计至少一个第一点扩散函数。
8.根据权利要求7所述的图像处理模块,其中,图像生成单元基于所述至少一个环境变量对波束成形的信号分类,并基于分类结果从点扩散函数数据库提取至少一个点扩散函数以估计至少一个第一点扩散函数。
9.根据权利要求4所述的图像处理模块,其中,图像生成单元使用波束成形的信号和估计的至少一个第一点扩散函数估计第二点扩散函数。
10.根据权利要求4所述的图像处理模块,其中,图像生成单元使用至少一个滤波器估计第一点扩散函数。
11.根据权利要求10所述的图像处理模块,其中,所述至少一个滤波器包括最小二乘滤波器(LSF)或倒频谱滤波器。
12.根据权利要求10所述的图像处理模块,其中,图像生成单元使用第一点扩散函数和至少一个第二点扩散函数估计变量来估计第二点扩散函数。
13.根据权利要求12所述的图像处理模块,其中,所述至少一个第二点扩散函数估计变量包括波束成形系数。
14.根据权利要求12所述的图像处理模块,其中,所述至少一个第二点扩散函数估计变量包括从超声波的声速、声速的改变、到对象中的目标部分的距离、输入信号生成模块的排列形式和每个通道的输入信号的衰减率中选择的至少一个。
15.根据权利要求4所述的图像处理模块,其中,第一点扩散函数是一维点扩散函数,第二点扩散函数是二维点扩散函数。
16.一种图像生成方法,包括:
将超声波发射到对象的目标部分,接收被对象的目标部分反射的回波超声波,并将接收的回波超声波转换为超声信号;
输出基于超声信号执行波束成形的超声信号;
基于波束成形的超声信号从点扩散函数数据库提取至少一个点扩散函数;
使用所述至少一个点扩散函数执行解卷积以从波束成形的超声信号生成图像,
其中,点扩散函数数据库包括基于用于波束成形的超声信号的至少一个环境变量预先获取的至少一个点扩散函数。
17.根据权利要求16所述的图像生成方法,其中,所述至少一个环境变量包括超声波的声速或到对象中的目标部分的距离。
18.根据权利要求16所述的图像生成方法,其中,基于波束成形的超声信号从点扩散函数数据库提取至少一个点扩散函数的步骤包括:基于所述至少一个环境变量对波束成形的超声信号分类,并基于分类结果从点扩散函数数据库提取至少一个点扩散函数。
19.根据权利要求16所述的图像生成方法,其中,使用所述至少一个点扩散函数执行解卷积以从波束成形的超声信号生成图像的步骤包括:
使用提取的至少一个点扩散函数估计另一个点扩散函数;
使用估计的另一个点扩散函数执行解卷积以从波束成形的超声信号生成图像。
20.根据权利要求19所述的图像生成方法,其中,所述至少一个环境变量包括超声波的声速或到对象中的目标部分的距离。
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201261676548P | 2012-07-27 | 2012-07-27 | |
US61/676,548 | 2012-07-27 | ||
KR1020130063737A KR102091536B1 (ko) | 2012-07-27 | 2013-06-03 | 영상 처리 모듈, 상기 영상 처리 모듈을 이용한 초음파 영상 장치 및 상기 초음파 영상 장치를 이용한 초음파 영상 생성 방법 |
KR10-2013-0063737 | 2013-06-03 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103565473A true CN103565473A (zh) | 2014-02-12 |
Family
ID=48906111
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310319539.9A Pending CN103565473A (zh) | 2012-07-27 | 2013-07-26 | 图像处理模块和图像处理方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US10660605B2 (zh) |
EP (1) | EP2690462B1 (zh) |
JP (1) | JP2014023934A (zh) |
CN (1) | CN103565473A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106132308A (zh) * | 2014-03-31 | 2016-11-16 | 皇家飞利浦有限公司 | 用于利用使用肋间空间的相干复合的声成像的系统和方法 |
CN111551943A (zh) * | 2020-05-19 | 2020-08-18 | 中国科学院声学研究所 | 一种基于damas2的稀疏布阵高分辨三维声成像方法及系统 |
CN111885964A (zh) * | 2018-03-08 | 2020-11-03 | 皇家飞利浦有限公司 | 识别超声图像的特性的系统和方法 |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014233602A (ja) * | 2013-06-05 | 2014-12-15 | キヤノン株式会社 | 被検体情報取得装置、被検体情報取得装置の制御方法 |
KR101832835B1 (ko) * | 2013-07-11 | 2018-02-28 | 삼성전자주식회사 | 영상 처리 모듈, 초음파 영상 장치, 영상 처리 방법 및 초음파 영상 장치의 제어 방법 |
KR20150010860A (ko) * | 2013-07-19 | 2015-01-29 | 삼성전자주식회사 | 초음파 영상 장치 및 그 제어 방법 |
KR20150050728A (ko) * | 2013-10-31 | 2015-05-11 | 삼성전자주식회사 | 초음파 영상 장치 및 그 제어 방법 |
KR20150118495A (ko) * | 2014-04-14 | 2015-10-22 | 삼성전자주식회사 | 초음파 프로브, 초음파 영상 장치 및 초음파 영상 장치를 제어하는 방법 |
KR20150118732A (ko) * | 2014-04-15 | 2015-10-23 | 삼성전자주식회사 | 초음파 장치 및 그 제어 방법 |
US9786057B2 (en) * | 2014-09-19 | 2017-10-10 | Lasertec Coporation | Inspection apparatus, coordinate detection apparatus, coordinate detection method, and wavefront aberration correction method |
US20160180143A1 (en) * | 2014-12-17 | 2016-06-23 | Intel Corporation | Eye tracking with 3d eye position estimations and psf models |
JP6912561B2 (ja) * | 2016-09-16 | 2021-08-04 | アンスティチュ ナショナル ドゥ ラ サンテ エ ドゥ ラ ルシェルシュ メディカル | 血液を伴う試料をイメージングするための方法および関連装置 |
CN112114045B (zh) * | 2020-08-03 | 2023-09-19 | 中国科学院大学 | 一种面向超声显微镜的分辨率增强方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1600274A (zh) * | 2003-09-24 | 2005-03-30 | 株式会社东芝 | 超析像处理装置与医用图像诊断装置 |
US20100054576A1 (en) * | 2008-08-26 | 2010-03-04 | Kazuhiro Tsujita | Image processing apparatus, image processing method, and image processing program |
EP2159753B1 (en) * | 2008-08-25 | 2010-10-27 | Sony Corporation | Image processing apparatus, imaging apparatus, image processing method, and program |
CN102388402A (zh) * | 2010-02-10 | 2012-03-21 | 松下电器产业株式会社 | 图像处理装置及图像处理方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6928182B1 (en) | 1998-10-15 | 2005-08-09 | Kui Ming Chui | Imaging |
US7720268B2 (en) | 2005-07-15 | 2010-05-18 | Siemens Corporation | System and method for ultrasound specific segmentation using speckle distributions |
US20070083114A1 (en) * | 2005-08-26 | 2007-04-12 | The University Of Connecticut | Systems and methods for image resolution enhancement |
KR100764414B1 (ko) | 2006-03-09 | 2007-10-05 | 삼성전기주식회사 | Psf 선택 모듈, 디지털 자동 초점 조절 장치 및 psf선택 방법 |
US7783433B2 (en) * | 2007-05-22 | 2010-08-24 | Honeywell International Inc. | Automated defect detection of corrosion or cracks using SAFT processed Lamb wave images |
KR101412752B1 (ko) * | 2007-11-26 | 2014-07-01 | 삼성전기주식회사 | 디지털 자동 초점 영상 생성 장치 및 방법 |
US8289438B2 (en) * | 2008-09-24 | 2012-10-16 | Apple Inc. | Using distance/proximity information when applying a point spread function in a portable media device |
JP2011035514A (ja) | 2009-07-30 | 2011-02-17 | Sony Corp | 画像処理装置および方法、並びにプログラム |
WO2012162413A2 (en) * | 2011-05-23 | 2012-11-29 | University Of Pittsburgh-Of The Commonwealth System Of Higher Education | Methods and apparatuses for measuring tissue stiffness changes using ultrasound elasticity imaging |
JP2013005258A (ja) * | 2011-06-17 | 2013-01-07 | Panasonic Corp | ブレ補正装置、ブレ補正方法及び帳票 |
KR101894391B1 (ko) * | 2011-10-05 | 2018-09-04 | 삼성전자주식회사 | 진단영상 생성장치, 의료영상시스템 및 빔포밍 수행방법 |
-
2013
- 2013-07-26 CN CN201310319539.9A patent/CN103565473A/zh active Pending
- 2013-07-26 JP JP2013155690A patent/JP2014023934A/ja active Pending
- 2013-07-29 EP EP13178349.0A patent/EP2690462B1/en active Active
- 2013-07-29 US US13/952,927 patent/US10660605B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1600274A (zh) * | 2003-09-24 | 2005-03-30 | 株式会社东芝 | 超析像处理装置与医用图像诊断装置 |
EP2159753B1 (en) * | 2008-08-25 | 2010-10-27 | Sony Corporation | Image processing apparatus, imaging apparatus, image processing method, and program |
US20100054576A1 (en) * | 2008-08-26 | 2010-03-04 | Kazuhiro Tsujita | Image processing apparatus, image processing method, and image processing program |
CN102388402A (zh) * | 2010-02-10 | 2012-03-21 | 松下电器产业株式会社 | 图像处理装置及图像处理方法 |
Non-Patent Citations (3)
Title |
---|
HO-CHUL SHIN ET AL.: "ESTIMATION OF AVERAGE SPEED OF SOUND USING DECONVOLUTION OF MEDICAL ULTRASOUND DATA", 《ULTRASOUND IN MED. & BIOL.》 * |
HO-CHUL SHIN ET AL.: "Estimation of speed of sound in dual-layered media using medical ultrasound image deconvolution", 《ULTRASONICS》 * |
JIRIK R ET AL.: "Two-dimensional blind Bayesian deconvolution of medical ultrasound images", 《IEEE TRANSACTIONS ON ULTRASONICS》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106132308A (zh) * | 2014-03-31 | 2016-11-16 | 皇家飞利浦有限公司 | 用于利用使用肋间空间的相干复合的声成像的系统和方法 |
US10746859B2 (en) | 2014-03-31 | 2020-08-18 | Koninklijke Philips N.V. | System and method for acoustic imaging with coherent compounding using intercostal spaces |
CN111885964A (zh) * | 2018-03-08 | 2020-11-03 | 皇家飞利浦有限公司 | 识别超声图像的特性的系统和方法 |
CN111551943A (zh) * | 2020-05-19 | 2020-08-18 | 中国科学院声学研究所 | 一种基于damas2的稀疏布阵高分辨三维声成像方法及系统 |
CN111551943B (zh) * | 2020-05-19 | 2022-07-12 | 中国科学院声学研究所 | 一种基于damas2的稀疏布阵高分辨三维声成像方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
JP2014023934A (ja) | 2014-02-06 |
US10660605B2 (en) | 2020-05-26 |
EP2690462B1 (en) | 2020-12-23 |
EP2690462A3 (en) | 2014-09-03 |
EP2690462A2 (en) | 2014-01-29 |
US20140031689A1 (en) | 2014-01-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103565473A (zh) | 图像处理模块和图像处理方法 | |
CN106875372B (zh) | 用于在医学图像中将结构分割的方法和系统 | |
US11488298B2 (en) | System and methods for ultrasound image quality determination | |
US20120065510A1 (en) | Ultrasound system and method for calculating quality-of-fit | |
US10856849B2 (en) | Ultrasound diagnostic device and ultrasound diagnostic device control method | |
US20130172749A1 (en) | Providing doppler spectrum images corresponding to at least two sample volumes in ultrasound system | |
EP2609866B1 (en) | Providing motion mode image in ultrasound system | |
KR101832835B1 (ko) | 영상 처리 모듈, 초음파 영상 장치, 영상 처리 방법 및 초음파 영상 장치의 제어 방법 | |
US11272906B2 (en) | Ultrasonic imaging device and method for controlling same | |
US20210369236A1 (en) | Methods and systems for performing transvalvular pressure quantification | |
CN110392553A (zh) | 用于定位声学传感器的定位设备和系统 | |
US20160310109A1 (en) | Ultrasonic imaging apparatus and control method therefor | |
US9375199B2 (en) | Ultrasonic imaging apparatus and image display method thereof | |
CN116709994A (zh) | 用于监测胎儿的心跳的装置 | |
CN112334072B (zh) | 用于执行彩色多普勒超声成像的方法和系统 | |
CN113316420A (zh) | 用于监测心脏的功能的方法和系统 | |
EP3914161B1 (en) | Methods and systems for investigating blood vessel characteristics | |
CN110420036B (zh) | 多类型探头兼容的数据处理方法、装置及可读存储介质 | |
KR102091536B1 (ko) | 영상 처리 모듈, 상기 영상 처리 모듈을 이용한 초음파 영상 장치 및 상기 초음파 영상 장치를 이용한 초음파 영상 생성 방법 | |
US20150015720A1 (en) | Image processor, ultrasonic imaging device, and image processing method | |
CN112930145A (zh) | 超声控制单元 | |
CN113382685B (zh) | 用于研究血管特性的方法和系统 | |
CN113543721B (zh) | 用于采集复合3d超声图像的方法和系统 | |
US20210128108A1 (en) | Loosely coupled probe position and view in ultrasound imaging | |
EP4132364A1 (en) | Methods and systems for obtaining a 3d vector flow field |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140212 |