CN103181780A - 信号分析方法、超声波图像分析方法以及超声波成像系统 - Google Patents
信号分析方法、超声波图像分析方法以及超声波成像系统 Download PDFInfo
- Publication number
- CN103181780A CN103181780A CN2012100732034A CN201210073203A CN103181780A CN 103181780 A CN103181780 A CN 103181780A CN 2012100732034 A CN2012100732034 A CN 2012100732034A CN 201210073203 A CN201210073203 A CN 201210073203A CN 103181780 A CN103181780 A CN 103181780A
- Authority
- CN
- China
- Prior art keywords
- signal
- mode function
- analyzed
- essential mode
- measured
- 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
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 64
- 238000003384 imaging method Methods 0.000 title claims abstract description 34
- 238000003703 image analysis method Methods 0.000 title abstract 2
- 238000012545 processing Methods 0.000 claims abstract description 44
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 5
- 238000012216 screening Methods 0.000 claims abstract 3
- 238000002604 ultrasonography Methods 0.000 claims description 43
- 238000006243 chemical reaction Methods 0.000 claims description 23
- 239000000463 material Substances 0.000 claims description 17
- 238000001914 filtration Methods 0.000 claims description 7
- 239000000203 mixture Substances 0.000 claims description 5
- 238000010009 beating Methods 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 4
- 210000004204 blood vessel Anatomy 0.000 description 22
- 210000004369 blood Anatomy 0.000 description 11
- 239000008280 blood Substances 0.000 description 11
- 238000000034 method Methods 0.000 description 7
- 238000007405 data analysis Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 239000000523 sample Substances 0.000 description 4
- 238000007493 shaping process Methods 0.000 description 4
- 210000001519 tissue Anatomy 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 241001465754 Metazoa Species 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000017531 blood circulation Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 241000272525 Anas platyrhynchos Species 0.000 description 1
- 208000012661 Dyskinesia Diseases 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 210000003743 erythrocyte Anatomy 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 230000017311 musculoskeletal movement, spinal reflex action Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000001766 physiological effect Effects 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 235000014347 soups Nutrition 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- 210000003462 vein Anatomy 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/06—Measuring blood flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/488—Diagnostic techniques involving Doppler signals
-
- 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/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
Abstract
一种信号分析方法、超声波图像分析方法以及超声波成像系统。所述信号分析方法包括如下步骤。接收待分析信号。利用经验模态分解及迭代的方式对于待分析信号进行筛选,以分解出至少一本质模态函数。对本质模态函数进行希尔伯特变换。经变换后的本质模态函数包括相位信息。利用相位处理手段来处理经变换后的本质模态函数,以获得包括角频率信息的本质模态函数。上述信号分析方法可应用于超声波成像系统,以辨识超声波图像信息。
Description
技术领域
本发明涉及一种信号分析方法,且特别涉及一种超声波图像信号的分析方法以及其应用。
背景技术
在已知技术中,超声波成像已广泛应用于各种领域中,例如工业、军事、医学等各种不同领域范围。其中,超声波成像系统于医学应用时,可用来测量人体组织的生理特性,诸如血管或血液的流速等。其利用多普勒(Doppler)原理将超声波能量传送到待测的人体组织区域,再接收由待测区域反射的超声波能量值。根据所反射的超声波能量,超声波成像系统可显示此待测区域的二维超声波图像。然而,对于欲检测血管中血液的低流速或较小血管管径而言,其信息必须搭配特定的信号分析方法来分析待测区域所反射的超声波能量值始可获得。
已知的信号分析方法是在频率域中对反射的超声波能量进行分析。然而,此种频率域(frequency domain)方法存在有部分缺失,不易克服。详细而言,根据多普勒原理,当超声波成像系统将超声波发射至待测区域内的待检测物质,例如血管中的红血球,而由该待检测物质反射的回波信号(echo),其频率会偏移,此频率偏移量与待检测物质的速度在发射方向上的分量成正比。在脉冲波(pulsed wave)多普勒模式中,超声波成像系统的探头向待测区域发射一系列的短脉冲,超声波成像系统所接收的反射信号可表示为二维的数据集,其中一维代表脉冲发射索引(即慢时间轴(slow time axis)),另一维则为飞行时间(即快时间轴(fast time axis))。慢轴上的信号即携有多普勒频移的讯息。因此,分析慢轴信号的相位即可获得待检测物质的图像信息。
已知的信号分析方法通常是计算两相邻慢轴的自相关(selfcorrelation)函数,以估计相位偏移(phase shift),并据此计算待检测物质的速度等信息。但此已知的信号分析方法是将慢轴上的信号经傅立叶变换(Fourier transform)后始进行计算,也就是在频率域中进行信号的分析。然而,目前已知的信号分析方法存在如底下的缺失。一是多普勒信号容易受到低频信号的干扰,例如血管壁的脉动、心脏脉搏、呼吸或不自主运动。而此方法必须利用wall filter高通滤波器来滤除噪声,始可进行后续的分析。设计理想的高通滤波器并非易事,从而提高了频率域分析方法的困难度。二是频率域分析方法的核心算法-傅立叶变换,其本质上为积分变换,时间域信号一旦经过变换,所有的时变讯息将完全遗失。因此已知的信号分析方法无法提供实际测量情况下,特定瞬时的频率信息。
发明内容
本发明的一实施例提出一种时间域(Time domain)的信号分析方法。所述信号分析方法包括如下步骤:(a)接收待分析信号。(b)利用经验模态分解(Empirical Mode Decomposition,EMD),以迭代方式(iterative method)对待分析信号进行筛选(iteratively sift),以分解出至少一本质模态函数(intrinsicfunction,IMF)。(c)对本质模态函数进行希尔伯特变换(normalized Hilberttransform)。经变换后的本质模态函数包含相位信息(phase information)。(d)利用相位处理手段来处理经变换后的本质模态函数,以获得包含角频率(angular frequency)信息的本质模态函数。
本发明的另一实施例提出一种超声波图像分析方法,其适用于超声波成像系统。所述超声波图像分析方法包括如下步骤:(a)接收待分析信号。(b)利用经验模态分解,以迭代的方式对于待分析信号进行筛选,以分解出至少一本质模态函数。(c)对本质模态函数进行希尔伯特变换。经变换后的本质模态函数包含相位信息。(d)利用相位处理手段来处理经变换后的本质模态函数,以获得包含角频率信息的本质模态函数。(e)比对经处理后的本质模态函数与相关待分析信号的超声波图像,以辨识超声波图像信息。
本发明的另一实施例提出一种超声波成像系统,其包括信号收发模块、信号处理模块以及图像显示模块。信号收发模块接收待分析信号。信号处理模块利用经验模态分解,以迭代方式对待分析信号进行筛选,以分解出至少一本质模态函数。信号处理模块对本质模态函数进行希尔伯特变换。经变换后的本质模态函数包含相位信息。信号处理模块利用相位处理手段来处理经变换后的本质模态函数,以获得包含角频率信息的本质模态函数。信号处理模块比对经处理后的本质模态函数与相关待分析信号的超声波图像,以辨识超声波图像的图像信息。图像显示模块根据信号处理模块对于待分析信号的处理结果来显示对应一待测区域的该超声波图像。
为让本发明的上述特征和优点能更明显易懂,下文特举实施例,并配合附图作详细说明如下。
附图说明
图1绘示本发明实施例的信号分析方法步骤流程图。
图2绘示本发明实施例的超声波成像系统方块示意图。
图3A绘示本发明实施例的仿体的超声波图像。
图3B绘示图3A的仿体的待分析信号的本质模态函数。
图4A绘示本发明另一实施例的仿体的超声波图像。
图4B绘示图4A的仿体待分析信号的本质模态函数。
【主要元件符号说明】
100:超声波成像系统
110:信号收发模块
112:换能器阵列单元
114:驱动单元
120:信号处理模块
122:数据处理单元
124:数据分析单元
130:图像显示模块
132:显示屏幕
200:待测区域
S110、S120、S130、S140:信号分析方法的步骤
W1、W2:血管的大小
D1、D2:血管壁与仿体表面的距离
具体实施方式
图1绘示本发明实施例的信号分析方法的步骤流程图。请参考图1,本实施例的信号分析方法在时间域处理信号,以获得时变(time-varying)的频率信息。首先,在步骤S110中,接收待分析信号。此待分析信号包含例如由待测区域内的待检测物质所反射的慢轴(Slow time axis)信号,其包括多普勒频移信息(doppler shift)。接着,在步骤S120中,利用经验模态分解,以迭代方式对待分析信号进行筛选,以分解出多个不同频段的本质模态函数。经验模态分解方法可以参考US20100092028A1、US6901353B1等文献,并不在此限。此步骤分解所得的结果包含至少一本质模态函数。在实施例中,如果根据此至少一本质模态函数即可分析获得待检测物质的图像信息时,则本发明并不局限于分解出多个不同频段的本质模态函数。在步骤S130中,对不同频段的本质模态函数进行希尔伯特变换。而经变换后的各本质模态函数包含相位信息。继之,在步骤S140中,利用相位处理手段来处理经变换后的各本质模态函数,以获得包含角频率信息的各本质模态函数。
在本实施例中,步骤S140的相位处理手段包括搭配曲线拟合(curve fitting)来滤除角频率信息中的突波(surge)等噪声,亦即角频率跳动信息的成分,以得到稳定的角频率信息。详细而言,本实施例的相位处理手段包括对各本质模态函数的相位信息进行微分,以获得对应的角频率信息。亦即,如果各本质模态函数的相位信息为ψ,角频率信息为ω,此相位处理手段首先对各本质模态函数的相位信息ψ进行微分,以获得角频率信息ω=dψ/dt。接着,此相位处理手段利用曲线拟合(curve fitting)来滤除角频率信息ω中的部分成分,例如突波等噪声。
简单来说,本实施例的信号分析方法直接分析待测区域所反射的慢轴信号,先对其进行经验模态分解,再利用希尔伯特变换修整所得的本质模态函数的包络线,进而直接求出各本质模态函数的相位信息,以计算其角频率信息。
此外,在一实施例中,所述信号分析方法还包括比对分析所得的本质模态函数与相关待分析信号的二维图像,以辨识二维图像上的图像信息,此图像信息包含待检测物质的位置、大小及移动速度。由于本实施例的信号分析方法采用时间域方法,因此其分析结果包含时变的频率信息,可即时掌握待测区域的实际测量情况。
本实施例的信号分析方法可广泛地应用在分析经希尔伯特变换修整包络线后的本质模态函数,并利用相位处理手段来获得稳定的角频率信息,此信号分析方法并不限于任何形式的信号检测暨处理系统。为更清楚地了解本发明,以下将配合附图,以本实施例的信号分析方法应用在超声波成像系统的范例实施例来作详细说明。
图2绘示本发明实施例的超声波成像系统的方块示意图。请参考图2,本实施例的超声波成像系统100包括信号收发模块110、信号处理模块120以及图像显示模块130。信号收发模块110包括换能器阵列单元112及驱动单元114。在此,信号收发模块110例如包括超声波成像系统100前端的超声波探头。在本实施例中,驱动单元114用以提供驱动信号至换能器阵列单元112内的各换能器,以触发各该换能器将超声波信号传递至待测区域200。换能器阵列单元112被触发后产生超声波信号,并发射该超声波信号至待测区域200。待测区域200通常可为生物体组织区域或其传输媒介(例如:人体组织区域/动物组织区域或是其传输媒介),在本实施例中例如是血管及在其中流动的血液。接着,由待测区域200反射出的超声波反射波再被换能器阵列单元112所接收,以使其后的信号处理模块120取得待测区域的原始数据,即待分析信号。
在本实施例中,换能器阵列单元112在接收到待测区域200所反射的反射波之后,可将该等模拟的反射波信号变换为数字信号。换句话说,本实施例的换能器阵列单元112可还包括模拟数字变换器,用以将该等模拟的反射波信号变换为数字信号,但本发明并不限于此。在其他实施例中,模拟数字变换功能也可由驱动单元114内部的电路来执行,或者由信号收发模块110与信号处理模块120之间的介面电路来执行,本发明并不加以限制。
接着,在本实施例中,信号处理模块120包括数据处理单元122以及数据分析单元124。在取得待分析信号之后,数据处理单元122对待分析信号的数据进行波束成形与聚焦操作。如前所述,换能器阵列单元112中的各换能器用以接收待测区域200反射出的反射波。对各该换能器而言,其延迟时间各不相同,因此需要一个波束成形电路来对各该反射波的图像信号进行不同的延迟。是以,数据处理单元122包括波束成形电路来适当地将微分延迟引入到所接收的各图像信号中,以动态地聚焦信号,进而产生出精确的待测区域二维超声波图像。
之后,数据处理单元122将波束成形后的信号自极坐标变换为直角坐标,以获得用直角坐标来表示信号强度及相位的I/Q数据(I/Q data)。进而,数据处理单元122将I/Q数据进行格式变换,以使I/Q数据可供图像显示模块130在其显示屏幕132上显示。在本实施例中,二维的超声波图像为超声波成像系统100在B扫描模式下检测待测区域200所得。
另一方面,为了进一步辨识超声波图像的图像信息,本实施例的数据分析单元124利用图1的信号分析方法对待测区域200反射的待分析信号进行经验模态分解。接着,数据分析单元124利用希尔伯特变换修整待分析信号的本质模态函数的包络线,进而在时间域利用前述的相位处理手段直接求出各本质模态函数的相位信息,以计算其角频率信息。之后,数据分析单元124比对经处理后的各本质模态函数与待分析信号的超声波图像,以辨识超声波图像的图像信息。在此,待分析信号为超声波成像系统100在A扫描模式下检测待测区域200所得。超声波图像的图像信息例如是生物体组织(例如:人、动物等)中血管的位置、大小及其中血液的流速。
另外,本公开的超声波图像的分析方法的各步骤的相关细节均可以由图1至图2实施例的叙述中获致足够的教示、建议与实施说明,因此不再赘述。另一方面,本公开例的超声波图像分析方法可用软件或可编程逻辑阵列(FPGA,Field Programmable Gate Array)芯片或其他软硬件等方式来实现,本发明并不加以限制。
图3A绘示本发明实施例的仿体(phantom)的超声波图像。图3B绘示图3A的仿体待分析信号的本质模态函数。其中图3A的横轴表示仿体的深度,纵轴表示仿体的横向位置;图3B的横轴表示仿体的横向位置,纵轴表示仿体各频段经正规化(normalized)的IMF强度。请参考图3A及图3B,图3A的仿体的超声波图像例如是由图2的超声波成像系统100在B扫描模式下检测待测区域200所得,其中仿体被配置于待测区域200中可模拟活体的软组织,因此适于体外实验的超声波成像。在图3A中,超声波成像系统100的超声波探头对待测区域的仿体发射超声波信号,以检测其二维超声波图像。此二维超声波图像包括血管位置、均匀仿体介质以及检测时所使用的标定点,其细节如图3A所绘示。
图3B所例示的不同频段的本质模态函数,为例如通过图1的信号分析方法处理所得的本质模态函数。如前所述,利用此方法可直接在时间域中处理各本质模态函数,以滤除信号突波等噪声。在本实施例中,比对图3B的多个本质模态函数与图3A的超声波图像,可辨识超声波图像上的血管的位置、大小及其中血液的流速。在实务操作上,超声波成像系统100的探头由例如图3A中的左侧向仿体发射超声波信号来获得反馈信号。因此,在图3B中,各本质模态函数在不同介质与界面中会产生信号变化值,如图中虚线所圈之处以及两直虚线包围住的区域。该等虚线所圈之处即是对应图3A的血管壁的图像,而两直虚线包围住的区域即是血管中血液流动的区域,其中宽度W1代表血管的大小;D1代表血管壁与仿体表面的距离,即血管的位置。此外,根据图3B中两直虚线包围住的区域的各本质模态函数的强度也可计算出血管内血液的流速。在本实施例中,血液的流速例如是每秒2.0毫升(2.0ml/s)。因此,利用此方法可即时的判断浅层及深层血管的位置与大小,并且可用于判断不同血管的流速值,以区别其种类,例如动脉、静脉或微血管。
图4A绘示本发明另一实施例的仿体超声波图像。图4B绘示图4A的仿体待分析信号的本质模态函数。其中图4A的横轴表示仿体的深度,纵轴表示仿体的位置;图4B的横轴表示仿体的位置,纵轴表示仿体各频段经正规化(normalized)的IMF强度。请参考图4A及图4B,图4A的二维超声波图像包括血管位置、均匀仿体介质以及检测时所使用的标定点,其细节如图4A所绘示。同样地,图4B所例示的不同频段的本质模态函数是通过图1的信号分析方法处理所得的本质模态函数。在本实施例中,比对图4B所例示的多个本质模态函数与图4A的超声波图像,同样可辨识出超声波图像上的血管的位置、大小及其中血液的流速。在图4B中的虚线所圈之处的信号变化值即是对应图4A的血管壁的图像,而两直虚线包围住的区域即是血管中血液流动的区域,其中宽度W2代表血管的大小;D2代表血管壁与仿体表面的距离,即血管的位置。此外,比较图3B与图4B中各对应的本质模态函数的强度可知,两者的强度并不相同,代表图3A与图4A的血管内血液的流速并不相同。在本实施例中,根据图4B中两直虚线包围住的区域的各本质模态函数的强度来计算,可得其血液的流速为例如每秒4.0毫升(4.0ml/s)。
综上所述,在本发明的范例实施例中,信号分析方法利用信号在不同介质与界面的不同特性,搭配频段拆解,分析信号传递至不同介质与界面所产生的频率变化值,以准确描绘出待测物的位置及大小。此外,所述信号分析方法利用相位处理手段,搭配曲线拟合在时间域中滤除角频率信息中的突波等噪声,以得到稳定的角频率信息,并可即时掌握待测区域的实际测量情况。
虽然本发明已以实施例公开如上,然其并非用以限定本发明,本领域技术人员,在不脱离本发明的精神和范围内,当可作些许的更动与润饰,故本发明的保护范围当视所附权利要求书所界定者为准。
Claims (24)
1.一种时间域的信号分析方法,包括:
接收待分析信号;
利用经验模态分解,以迭代方式对该待分析信号进行筛选,以分解出至少一本质模态函数;
对该至少一本质模态函数进行希尔伯特变换,其中经变换后的该至少一本质模态函数包含相位信息;以及
利用相位处理手段来处理经变换后的该至少一本质模态函数,以获得包含角频率信息的该至少一本质模态函数。
2.如权利要求1所述的信号分析方法,还包括:
比对经处理后的该至少一本质模态函数与相关该待分析信号的二维图像,以辨识该二维图像上的图像信息。
3.如权利要求1所述的信号分析方法,其中利用该相位处理手段来处理经变换后的该至少一本质模态函数的步骤在时间域内进行。
4.如权利要求1所述的信号分析方法,其中利用该相位处理手段处理经变换后的该至少一本质模态函数的步骤包括:
对该至少一本质模态函数的该相位信息进行微分,以获得该角频率信息;以及
利用曲线拟合手段来滤除该角频率跳动信息的成分。
5.如权利要求1所述的信号分析方法,其中在利用经验模态分解筛选该待分析信号的步骤中分解所得的该至少一本质模态函数包括多个不同频段的本质模态函数。
6.如权利要求2所述的信号分析方法,还包括:
发射检测信号至待测区域;以及
接收该待测区域所反射的该待分析信号,
其中相关该待分析信号的该二维图像包括该待测区域的图像。
7.如权利要求6所述的信号分析方法,其中该待测区域内含有待检测物质,该待测区域所反射的该待分析信号包含该待检测物质的多普勒频移信息。
8.如权利要求7所述的信号分析方法,其中该二维图像上的图像信息包含该待检测物质的移动速度及位置。
9.一种超声波图像的分析方法,适用于超声波成像系统,该分析方法包括:
接收待分析信号;
利用经验模态分解,以迭代方式对该待分析信号进行筛选,以分解出至少一本质模态函数;
对该至少一本质模态函数进行希尔伯特变换,其中经变换后的该至少一本质模态函数包含相位信息;
利用相位处理手段来处理经变换后的该至少一本质模态函数,以获得包含角频率信息的该至少一本质模态函数;以及
比对经处理后的该至少一本质模态函数与相关该待分析信号的超声波图像,以辨识该超声波图像信息。
10.如权利要求9所述的超声波图像分析方法,其中利用该相位处理手段来处理经变换后的该至少一本质模态函数的步骤在时间域内进行。
11.如权利要求9所述的超声波图像分析方法,其中利用该相位处理手段处理经变换后的该至少一本质模态函数的步骤包括:
对该至少一本质模态函数的该相位信息进行微分,以获得该角频率信息;以及
利用曲线拟合手段来滤除该角频率跳动信息的成分。
12.如权利要求9所述的超声波图像分析方法,其中在利用经验模态分解筛选该待分析信号的步骤中分解所得的该至少一本质模态函数包括多个不同频段的本质模态函数。
13.如权利要求9所述的超声波图像分析方法,还包括:
发射检测信号至待测区域;以及
接收该待测区域所反射的该待分析信号,
其中相关该待分析信号的该超声波图像包括该待测区域的图像。
14.如权利要求13所述的超声波图像分析方法,其中该待测区域含有待检测物质,该待测区域所反射的该待分析信号包含该待检测物质的多普勒频移信息。
15.如权利要求14所述的超声波图像分析方法,其中该超声波图像的图像信息包含该待检测物质的位置、大小及移动速度。
16.如权利要求9所述的超声波图像分析方法,其中该待分析信号为在A扫描模式或B扫描模式下检测待测区域所得。
17.一种超声波成像系统,包括:
信号收发模块,接收待分析信号;
信号处理模块,利用经验模态分解,以迭代方式对该待分析信号进行筛选,以分解出至少一本质模态函数;对该至少一本质模态函数进行希尔伯特变换,其中经变换后的该至少一本质模态函数包含相位信息;利用相位处理手段来处理经变换后的该至少一本质模态函数,以获得包含角频率信息的该至少一本质模态函数;以及比对经处理后的该至少一本质模态函数与相关该待分析信号的超声波图像,以辨识该超声波图像信息;以及
图像显示模块,根据该信号处理模块对该待分析信号的处理结果来显示对应待测区域的该超声波图像。
18.如权利要求17所述的超声波成像系统,其中该信号分析模块在时间域内利用该相位处理手段来处理经变换后的该至少一本质模态函数。
19.如权利要求17所述的超声波成像系统,其中该信号分析模块利用该相位处理手段处理经变换后的该至少一本质模态函数时,对该至少一本质模态函数的该相位信息进行微分,以获得该角频率信息;以及利用曲线拟合手段来滤除该角频率跳动信息的成分。
20.如权利要求17所述的超声波成像系统,其中该信号分析模块分解所得的该至少一本质模态函数包含多个不同频段的本质模态函数。
21.如权利要求17所述的超声波成像系统,其中信号收发模块更发射检测信号至该待测区域;以及接收该待测区域所反射的该待分析信号,其中相关该待分析信号的该超声波图像包含该待测区域的图像。
22.如权利要求21所述的超声波成像系统,其中该待测区域含有待检测物质,该待测区域所反射的该待分析信号包括该待检测物质的多普勒频移信息。
23.如权利要求22所述的超声波成像系统,其中该超声波图像的图像信息包括该待检测物质的位置、大小及移动速度。
24.如权利要求17所述的超声波成像系统,其中该待分析信号为超声波成像系统在A扫描模式或B扫描模式下检测该待测区域所得。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW100148921A TWI482613B (zh) | 2011-12-27 | 2011-12-27 | 訊號分析方法、超音波影像分析方法以及超音波成像系統 |
TW100148921 | 2011-12-27 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103181780A true CN103181780A (zh) | 2013-07-03 |
CN103181780B CN103181780B (zh) | 2015-08-05 |
Family
ID=48654600
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210073203.4A Active CN103181780B (zh) | 2011-12-27 | 2012-03-19 | 信号分析方法、超声波图像分析方法以及超声波成像系统 |
Country Status (3)
Country | Link |
---|---|
US (1) | US8837798B2 (zh) |
CN (1) | CN103181780B (zh) |
TW (1) | TWI482613B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108697354A (zh) * | 2016-02-23 | 2018-10-23 | 梅约医学教育与研究基金会 | 超声血流成像 |
Families Citing this family (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103576131B (zh) * | 2013-11-02 | 2016-02-10 | 中国人民解放军海军航空工程学院 | 距离扩展目标内禀模态特征能量智能融合检测方法 |
CN103605120B (zh) * | 2013-11-02 | 2016-03-02 | 中国人民解放军海军航空工程学院 | 距离扩展目标内禀模态局部奇异值有序融合检测方法 |
CN103558022B (zh) * | 2013-11-08 | 2017-02-08 | 国网福建晋江市供电有限公司 | 基于振动信号聚类的配电开关机械状态诊断方法 |
TWI506583B (zh) * | 2013-12-10 | 2015-11-01 | 國立中央大學 | 分析系統及其方法 |
KR102243022B1 (ko) * | 2014-03-05 | 2021-04-21 | 삼성메디슨 주식회사 | 선택 정보에 기초하여 관심 영역에 포함된 혈류에 대한 정보를 출력하는 방법, 장치 및 시스템. |
CN104138277B (zh) * | 2014-07-29 | 2016-03-16 | 哈尔滨工业大学 | 定量超声系统中基于经验模态分解的组织微观结构检测方法 |
US10271821B2 (en) | 2014-12-23 | 2019-04-30 | Industrial Technology Research Institute | Method of ultrasound imaging and ultrasound scanner |
CN105069440A (zh) * | 2015-08-19 | 2015-11-18 | 杭州电子科技大学 | 基于na-memd和互信息的脑电imf分量选择方法 |
TWI567732B (zh) * | 2015-10-05 | 2017-01-21 | Jin-Feng Lin | A Time-Frequency Energy Distribution Characteristic Analysis Method for Underwater Audio Signal Based on Empirical Mode Disassembly |
CN105676202B (zh) * | 2016-01-22 | 2018-04-17 | 西安电子科技大学 | 基于局部均值分解的目标识别特征提取方法 |
CN107480412A (zh) * | 2016-06-06 | 2017-12-15 | 杭州慧牧科技有限公司 | 一种畜牧秤称重方法 |
US11286765B2 (en) * | 2016-12-07 | 2022-03-29 | Halliburton Energy Services, Inc. | Measuring invisible lost time in drilling operations |
EP3665475B1 (en) | 2017-08-10 | 2023-11-29 | Mayo Foundation for Medical Education and Research | Shear wave elastography with ultrasound probe oscillation |
CN107703464A (zh) * | 2017-09-13 | 2018-02-16 | 武汉科技大学 | 一种动态磁通信号处理方法 |
CN108182415B (zh) * | 2017-12-29 | 2021-05-04 | 哈尔滨工业大学 | 基于自滤波变频经验模态分解获得时频分布的方法 |
TWI682169B (zh) * | 2018-03-29 | 2020-01-11 | 佳世達科技股份有限公司 | 超音波成像方法 |
CN110175508B (zh) * | 2019-04-09 | 2021-05-07 | 杭州电子科技大学 | 一种应用于超声波局放检测的特征值提取方法 |
CN110796047B (zh) * | 2019-10-18 | 2022-08-09 | 哈尔滨工业大学 | 一种基于机器学习的自适应稀疏时频分析方法 |
CN111079710B (zh) * | 2019-12-31 | 2023-04-18 | 江苏理工学院 | 一种基于改进ceemd滚动轴承信号的多层降噪方法 |
CN111275022B (zh) * | 2020-03-19 | 2023-05-16 | 山东宜佳成新材料有限责任公司 | 基于遗忘因子型经验模态分解的污渍检测分析方法及应用 |
CN111985342B (zh) * | 2020-07-23 | 2022-08-26 | 河海大学 | 一种基于经验模态分解的海杂波时间相关性分析方法 |
CN111899724A (zh) * | 2020-08-06 | 2020-11-06 | 中国人民解放军空军预警学院 | 基于希尔伯特黄变换的语音特征系数提取方法及相关设备 |
CN113040741A (zh) * | 2021-03-01 | 2021-06-29 | 中国农业大学 | 作物根区多频电阻抗成像方法及系统 |
CN115255405B (zh) * | 2022-09-23 | 2022-12-23 | 相国新材料科技江苏有限公司 | 一种增材制造设备的智能控制方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6381559B1 (en) * | 1996-08-12 | 2002-04-30 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Empirical mode decomposition apparatus, method and article of manufacture for analyzing biological signals and performing curve fitting |
TW594032B (en) * | 2003-07-03 | 2004-06-21 | Chung Shan Inst Of Science | ROI segmentation processing method for digital ultrasound image |
US6901353B1 (en) * | 2003-07-08 | 2005-05-31 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Computing Instantaneous Frequency by normalizing Hilbert Transform |
TWI252929B (en) * | 2004-08-31 | 2006-04-11 | Bai-Chi Li | Ultrasonic imaging method and device thereof |
CN101822548A (zh) * | 2010-03-19 | 2010-09-08 | 哈尔滨工业大学(威海) | 基于相关性分析与经验模态分解的超声信号去噪方法 |
Family Cites Families (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5582171A (en) | 1994-07-08 | 1996-12-10 | Insight Medical Systems, Inc. | Apparatus for doppler interferometric imaging and imaging guidewire |
JP3406106B2 (ja) | 1995-02-06 | 2003-05-12 | ジーイー横河メディカルシステム株式会社 | 超音波画像表示方法および超音波診断装置 |
US6311130B1 (en) * | 1996-08-12 | 2001-10-30 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Computer implemented empirical mode decomposition method, apparatus, and article of manufacture for two-dimensional signals |
US5983162A (en) * | 1996-08-12 | 1999-11-09 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Computer implemented empirical mode decomposition method, apparatus and article of manufacture |
US6738734B1 (en) * | 1996-08-12 | 2004-05-18 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Empirical mode decomposition apparatus, method and article of manufacture for analyzing biological signals and performing curve fitting |
US6068599A (en) | 1997-07-14 | 2000-05-30 | Matsushita Electric Industrial Co., Ltd. | Blood vessel puncturing device using ultrasound |
US6464641B1 (en) | 1998-12-01 | 2002-10-15 | Ge Medical Systems Global Technology Company Llc | Method and apparatus for automatic vessel tracking in ultrasound imaging |
US6423006B1 (en) | 2000-01-21 | 2002-07-23 | Siemens Medical Solutions Usa, Inc. | Method and apparatus for automatic vessel tracking in ultrasound systems |
US6482161B1 (en) | 2000-06-29 | 2002-11-19 | Acuson Corporation | Medical diagnostic ultrasound system and method for vessel structure analysis |
US6862558B2 (en) * | 2001-02-14 | 2005-03-01 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Empirical mode decomposition for analyzing acoustical signals |
US7137963B2 (en) | 2002-08-26 | 2006-11-21 | Flowcardia, Inc. | Ultrasound catheter for disrupting blood vessel obstructions |
US7464006B1 (en) * | 2003-10-09 | 2008-12-09 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Analyzing nonstationary financial time series via hilbert-huang transform (HHT) |
TWI268149B (en) | 2004-12-29 | 2006-12-11 | Ind Tech Res Inst | Intravascular ultrasound probing device |
EP2712553A3 (en) | 2005-01-11 | 2014-09-17 | Volcano Corporation | Vascular image co-registration |
US20060184029A1 (en) | 2005-01-13 | 2006-08-17 | Ronen Haim | Ultrasound guiding system and method for vascular access and operation mode |
KR100752333B1 (ko) | 2005-01-24 | 2007-08-28 | 주식회사 메디슨 | 3차원 초음파 도플러 이미지의 화질 개선 방법 |
US7941298B2 (en) * | 2006-09-07 | 2011-05-10 | DynaDx Corporation | Noise-assisted data analysis method, system and program product therefor |
JP5014051B2 (ja) | 2007-10-09 | 2012-08-29 | 株式会社ユネクス | 血管超音波画像測定方法 |
US9826959B2 (en) | 2008-11-04 | 2017-11-28 | Fujifilm Corporation | Ultrasonic diagnostic device |
TWI399195B (zh) | 2008-12-16 | 2013-06-21 | Ind Tech Res Inst | 動態立體血管超音波影像系統及其操作方法 |
US20100207942A1 (en) | 2009-01-28 | 2010-08-19 | Eigen, Inc. | Apparatus for 3-d free hand reconstruction |
DE102009012821A1 (de) | 2009-03-13 | 2010-09-16 | Wittenstein Ag | Verfahren zur simultanen Bestimmung von Position und Geschwindigkeit von Objekten |
JP5606143B2 (ja) | 2009-06-08 | 2014-10-15 | 株式会社東芝 | 超音波診断装置、画像処理装置、画像処理方法および画像表示方法 |
US8861811B2 (en) | 2009-08-19 | 2014-10-14 | Merck Sharp & Dohme Corp. | System and method for segmenting M-mode ultrasound images showing blood vessel wall motion over time |
DE102009043069A1 (de) | 2009-09-25 | 2011-04-07 | Siemens Aktiengesellschaft | Visualisierungsverfahren und Bildgebungssystem |
US9445780B2 (en) | 2009-12-04 | 2016-09-20 | University Of Virginia Patent Foundation | Tracked ultrasound vessel imaging |
-
2011
- 2011-12-27 TW TW100148921A patent/TWI482613B/zh active
-
2012
- 2012-03-19 CN CN201210073203.4A patent/CN103181780B/zh active Active
- 2012-07-19 US US13/553,793 patent/US8837798B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6381559B1 (en) * | 1996-08-12 | 2002-04-30 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Empirical mode decomposition apparatus, method and article of manufacture for analyzing biological signals and performing curve fitting |
TW594032B (en) * | 2003-07-03 | 2004-06-21 | Chung Shan Inst Of Science | ROI segmentation processing method for digital ultrasound image |
US6901353B1 (en) * | 2003-07-08 | 2005-05-31 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Computing Instantaneous Frequency by normalizing Hilbert Transform |
TWI252929B (en) * | 2004-08-31 | 2006-04-11 | Bai-Chi Li | Ultrasonic imaging method and device thereof |
CN101822548A (zh) * | 2010-03-19 | 2010-09-08 | 哈尔滨工业大学(威海) | 基于相关性分析与经验模态分解的超声信号去噪方法 |
Non-Patent Citations (1)
Title |
---|
YUFENG ZHANG,ECT: "The Removal of Wall Components in Doppler Ultrasound Signals by Using Empirical Mode Decomposition Algorithm", 《IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING》, vol. 54, no. 9, 31 December 2007 (2007-12-31) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108697354A (zh) * | 2016-02-23 | 2018-10-23 | 梅约医学教育与研究基金会 | 超声血流成像 |
CN108697354B (zh) * | 2016-02-23 | 2022-12-06 | 梅约医学教育与研究基金会 | 超声血流成像 |
Also Published As
Publication number | Publication date |
---|---|
TWI482613B (zh) | 2015-05-01 |
US8837798B2 (en) | 2014-09-16 |
CN103181780B (zh) | 2015-08-05 |
TW201325555A (zh) | 2013-07-01 |
US20130163839A1 (en) | 2013-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103181780B (zh) | 信号分析方法、超声波图像分析方法以及超声波成像系统 | |
CN102458260B (zh) | 测量心脏收缩功能的方法和装置 | |
US7535797B2 (en) | High-resolution ultrasound displacement measurement apparatus and method | |
US8684934B2 (en) | Adaptively performing clutter filtering in an ultrasound system | |
KR102051293B1 (ko) | 의료 초음파 전단파 이미징에서의 분류 프리프로세싱 | |
US20150148673A1 (en) | Ultrasonic diagnostic apparatus and elasticity evaluation method | |
CN104510499B (zh) | 医学超声成像中的剪切波检测 | |
KR102025328B1 (ko) | 평면파 합성을 이용한 초음파 벡터 도플러 영상의 생성 장치 및 방법 | |
CN103908289B (zh) | 消除剪切波中的背景噪声的方法和相应的超声成像系统 | |
US9261485B2 (en) | Providing color doppler image based on qualification curve information in ultrasound system | |
CN108324262B (zh) | 用于动脉脉搏测量的方法和装置 | |
CN107049361B (zh) | 使用剪切波的声音速度成像 | |
US10448925B2 (en) | Ultrasonic diagnostic apparatus and method for reducing clutter | |
US9746322B2 (en) | Thickness measuring device and thickness measuring method | |
JP2015062016A (ja) | 多数の超音波パルス形状を使用した検出装置 | |
WO2020042020A1 (zh) | 一种超声弹性检测设备及剪切波弹性成像方法、装置 | |
CN114502063B (zh) | 用于确定感兴趣对象的材料类型的方法和系统 | |
CN107427281B (zh) | 超声波诊断装置 | |
CN110392553B (zh) | 用于定位声学传感器的定位设备和系统 | |
EP3613350A1 (en) | Systems and method for performing pulse wave velocity measurements | |
EP3482692A1 (en) | Spectrum analysis device, spectrum analysis method, and ultrasonic imaging device | |
WO2016210246A2 (en) | Ultrasound image-based concentration measurement | |
WO2015008852A1 (en) | Doppler ultrasound diagnosis device, doppler ultrasound diagnosis method, and program | |
KR102132072B1 (ko) | 초음파 시스템 및 그 동작 방법 | |
Rodriguez-Molares et al. | Reconstruction of specular reflectors by iterative image source localization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |