CN101606851A - 超声波诊断装置 - Google Patents
超声波诊断装置 Download PDFInfo
- Publication number
- CN101606851A CN101606851A CNA200910147985XA CN200910147985A CN101606851A CN 101606851 A CN101606851 A CN 101606851A CN A200910147985X A CNA200910147985X A CN A200910147985XA CN 200910147985 A CN200910147985 A CN 200910147985A CN 101606851 A CN101606851 A CN 101606851A
- Authority
- CN
- China
- Prior art keywords
- signal
- diagnostic ultrasound
- target location
- ultrasound equipment
- received signal
- 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
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
- 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/52023—Details of receivers
- G01S7/5203—Details of receivers for non-pulse systems, e.g. CW systems
- G01S7/52031—Extracting wanted echo signals
-
- 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/8959—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using coded signals for correlation purposes
-
- 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/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- 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
-
- 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/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/50—Systems of measurement, based on relative movement of the target
- G01S15/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
- G01S15/586—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种超声波诊断装置。FSK调制器20通过基于模式发生器24所提供的周期信号序列的频移键控(FSK)的方式,对RF波振荡器22所提供的RF波进行FSK调制处理,从而产生FSK连续波。从FSK调制器20输出的连续波在延迟电路26I和延迟电路26Q中被延迟,然后作为参考信号提供给接收混频器30的每个混频器。延迟电路26I和延迟电路26Q中的每一个以根据目标位置的深度而定的延迟量对连续波进行延迟,并输出延迟的参考信号。因此,在增加来自目标位置的接收信号与参考信号之间相关性的情况下执行解调处理,并通过多普勒信息分析单元44选择性地提取目标位置的多普勒信息。
Description
技术领域
本发明涉及一种超声波诊断装置,更具体地,涉及一种采用调制连续波的超声波诊断装置。
背景技术
连续波多普勒是已知的一种使用连续波的超声波诊断装置技术。在连续波多普勒技术中,数兆赫兹的正弦波形的发射波连续地辐射到活体中,然后连续接收活体中的反射波。反射波包括活体内的运动元素(例如血流)产生的多普勒频移信息。因此,通过提取多普勒频移信息并对其进行频率分析,能够形成例如反映运动元素的速度信息的多普勒波形。
为了快速获得速度测量值,使用连续波的连续波多普勒技术通常优于使用脉冲波的脉冲多普勒。在这种情况下,本申请的发明人进行了连续波多普勒技术的相关研究。在其中一个研究成果中,本申请的发明人提出了公开在JP2005-253949A中的调频连续波多普勒(FMCW多普勒)的相关技术。
由此典型的连续波多普勒设备(即不使用FMCW多普勒的设备)不能执行位置测量。在这点上,本申请的发明人在JP2006-14916A中提出了一种技术,通过使用FMCW多普勒,这种技术不但能测量活体内组织的速度,还能测量活体内的组织的位置。
上述公布所述的FMCW多普勒技术是使新形式的超声波诊断成为可能的革命性技术。本申请的发明人持续研究并改进这项里程碑式的技术。
发明内容
考虑到上述情形,提出了本发明,其有利地提供了一种通过使用连续波提取目标位置的活体内(in-vivo)信息的改进技术。
为了获得上述优点,根据本发明的一个方面,提供了一种超声波诊断装置,其包括:发射信号处理单元,其输出发射信号,所述发射信号为基于周期信号序列进行了数字调制的连续波;发射/接收单元,其向活体发射与发射信号对应的发射波,并从活体接收与发射波关联的接收波,从而获得接收信号;接收信号处理单元,其通过使用基于发射信号形成的参考信号,对接收信号进行解调处理,从而获得解调信号;以及活体内信息提取单元,其从解调信号中提取活体内信息,其中选择性地提取活体内目标位置的活体内信息。
在以上方面中,通过调整从活体内目标位置获得的接收信号的周期信号序列与参考信号的周期信号序列之间的相关性,可以提取例如作为与参考信号的相关度较高的信号分量的来自目标位置的接收信号。另外,通过使用带通滤波器或低通滤波器从接收信号中提取多普勒信息,可以实现例如对目标位置的活体内信息的选择性提取。此处,在以上方面中,频移键控、相移键控、幅移键控等等可用作数字调制处理。此外,理想情况是,参考信号的波形和发射信号的波形完全相同。但是,参考信号和发射信号也可以是认为其波形基本相同的对应关系。
根据本发明,可以通过使用数字调制的连续波实现目标位置的多普勒信息的提取。
附图说明
参照下图,将详细描述本发明的优选实施例,其中:
图1是图解说明根据本发明的超声波诊断装置的整体结构的功能框图;
图2是说明发射信号的图,其中所述发射信号是通过FSK形成的连续波;
图3是说明FSK调制器输出的连续波的相位变化的图;
图4是说明FSK调制器输出的连续波的频谱的图;
图5是显示参考信号的相位与接收信号的相位匹配的状态的图;
图6是显示将参考信号的相位延迟T的状态的图;
图7是显示将参考信号的相位延迟2T的状态的图;
图8是图解说明参考波的相位和乘法器的输出之间的关系的曲线图;
图9是图解说明通过混频处理得到的解调信号的频谱的图;
图10是图解说明模式发生器的具体示例的图;
图11是说明延迟电路和接收混频器的变化示例的图。
具体实施方式
下面将参照附图详细描述本发明的优选实施例。
图1是图解说明根据本发明优选实施例的超声波诊断装置的整体结构的功能框图。发射换能器10向活体内连续发射发射波,接收换能器12从活体内连续接收反射波。因此,通过不同的换能器进行发射和接收,从而实现了采用称为连续波多普勒方法技术的发射/接收。本实施例中采用的连续波是通过FSK调制器20形成的数字调制的连续波。
通过基于模式发生器24提供的周期信号序列的频移键控(Frequency Shift Keying,FSK)的方式,FSK调制器20对RF(射频)波振荡器22提供的RF波进行调制处理,从而产生FSK连续波。下文将描述通过频移键控(FSK)形成的连续波的波形。FSK调制器20向功率放大器14输出经过数字(FSK)调制的连续波。
功率放大器14对FSK调制的连续波进行功率放大,然后将经过功率放大的FSK调制的连续波提供给发射换能器10。发射换能器10发射与FSK调制的连续波对应的发射波。之后,接收换能器12连续接收来自活体内的反射波。
前置放大器16对接收换能器12所提供的接收波信号进行接收处理,例如低噪声放大,从而产生接收RF信号,该接收RF信号被输出到接收混频器30。接收混频器30是用于对接收RF信号进行正交检波以产生复基带信号的电路,该接收混频器30由两个混频器32和34组成。每个混频器均为将接收RF信号与预设的参考信号混频的电路。
基于数字调制(FSK调制)的连续波(即发射信号)产生提供给接收混频器30的每个混频器的参考信号。具体而言,从FSK调制器20输出的连续波分别经延迟电路26I和26Q延迟。将经延迟电路26I延迟的连续波提供给混频器32,将经延迟电路26Q延迟的连续波提供给混频器34。
延迟电路26I和延迟电路26Q的每一个以根据目标位置的深度而定的延迟量对连续波进行延迟,从而输出延迟的参考信号。延迟电路26I和延迟电路26Q中的每一个均可由例如n级移位寄存器构成。在这种情况下,从移位寄存器的n级抽头(tap)中选择用于与目标位置的深度对应的延迟量的抽头,从而输出与目标位置的深度对应的参考信号(即延迟的连续波)。
这里,延迟电路26I和延迟电路26Q通过将各连续波的相位相对于彼此移动π/2来进行延迟处理。因此,混频器32输出同相(in phase)信号分量(I信号分量),混频器34输出正交(quadrature)信号分量(Q信号分量)。之后,通过设置在接收混频器30下游的LPF(低通滤波器)36和38分别去除同相信号分量和正交信号分量中的高频分量,从而能够提取检波后仅具有必要带宽的解调信号。
如将在下面的本发明的技术原理说明中所详述的,接收混频器输出信号(即解调信号)包括大量来自目标位置的接收信号分量,该接收混频器输出信号是接收RF信号与参考信号在每个混频器中进行混频的结果。LPF 36和38提取来自目标位置的接收信号分量中所包含的直流(低频)信号分量(如果目标位置出现移动元素,则包含多普勒信号分量)。
FFT电路(快速傅立叶变换电路)40和42对每个解调信号(同相信号分量和正交信号分量)执行FFT操作。因此,解调信号在FFT电路40和42的每一个中变换为频谱。这里,FFT电路40和42输出的频谱是以频谱数据的形式提供的,其频率分辨率δf取决于电路设置条件等等。
多普勒信息分析单元44从变换成频谱的解调信号中提取包含极性的多普勒信息。这时,由于延迟电路26I和延迟电路26Q已根据活体内目标位置的深度调整了参考信号和接收信号之间的延迟关系,所以可以选择性地提取目标位置的多普勒信息。将参照本发明的原理说明来详细描述延迟调整与目标位置的多普勒信息的提取之间的关系。多普勒信息分析单元44针对活体内每个深度(每个位置)提取多普勒信息,从而对沿超声波束(声射线)的每个深度计算活体内组织的速度,并实时输出结果。这里,通过扫描超声波束,可以以二维或三维的方式计算在活体内每个位置的组织的速度。
显示处理单元46基于每个深度(位置)的活体内组织的速度建立多普勒波形或曲线图,该多普勒波形或曲线图包括例如与深度速度有关的信息,并令显示单元48实时显示以这种方式建立的多普勒波形和曲线图。这里,系统控制单元50控制图1所示的超声波诊断装置中的每个单元。具体而言,系统控制单元50执行发射控制、接收控制、显示控制等等。
如上文所概述,根据本实施例,发射并接收与经过FSK调制的连续波对应的超声波,以获得接收信号,并根据活体内目标位置的深度调整参考信号与接收信号之间的延迟关系,之后执行解调处理,来自目标位置的接收信号与参考信号之间的相关度因此增加,从而可以选择性地提取目标位置的活体内信息。这里,将详细描述选择性地提取目标位置的多普勒信息的原理。
图2说明的发射信号是通过频移键控(FSK)方式形成的连续波。图2(I)图解说明了模式发生器(在图1中以附图标记24表示)输出的周期信号序列的示例。该模式发生器产生例如值随机变化的二进制码(伪随机信号),如图2(I)所示。可以使用以脉冲压缩等等方式实现的码序列,例如PN(伪噪声)序列、M序列或格雷(Gorey)序列,作为是周期信号序列的伪随机信号的码序列。
图2(II)图解说明了在FSK调制器(图1中以附图标记20表示)中形成的调制连续波(即发射信号)。FSK调制器通过基于如图2(I)所示的二进制码的频移键控的方式对RF波(载波)进行FSK调制处理,从而产生是连续波的发射信号。FSK调制器通过例如在二进制码为“1”的位(bit)周期设置频率f1,并在二进制码为“0”的位周期设置频率f2,从而产生图2(II)所示的发射信号。
如此,例如,发射换能器(在图1中以附图标记10表示)发射超声波,该超声波是与图2(II)所示的发射信号对应的连续波,并且通过接收换能器(图1中以附图标记12表示)得到来自活体的接收信号。
图3是说明FSK调制器输出的连续波的相位变化的图。在图3中,载波的相位始终用固定值0表示,并表示了多个对应于频率f1或f2的相位旋转。为了对FSK调制信号进行解调,只需要区分信号的频率。在图2(II)所示的连续波的情况下,可通过区分频率f1和频率f2对图2(II)所示的信号进行解调以获得图2(I)所示的二进制码。如果只区分频率就足够了,则不需要特别地确定对应于符号的每个频率的相位。例如,如果只需要区分频率,那么符号过渡时的非连续的相位变化不会造成任何问题。然而,符号过渡时的非连续相位会因此造成频带扩展。因此,通过确定每个频率和1位的时间长度T,从而使相位在符号过渡时连续变化,则可在较窄的带宽中实现FSK调制和解调。
图3图解说明了满足上述条件时的连续波的相位变化。在图3所示的连续波的相位变化中,当二进制码的符号变化时,相位连续变化;在二进制码为“1”的位周期中(在频率f1的周期中),连续波的相位变化为(即提前),在二进制码为“0”的位周期中(在频率f2的周期中),连续波的相位变化为(即延迟)。
图3所示的示例性连续波(FSK调制波)呈现这样的特性:不论符号如何变化,连续波的相位都连续地变化,从而使相位总在图3所示的曲线图的上升虚线或下降虚线上移动。进一步地,通过将对应于符号的频率空间设置为某较小值,可以得到窄带的FSK连续波。频率空间越窄,1位中的相位旋转越小。通过设置频率空间使得该相位旋转量为π/2,即通过将连续波的相位在1位周期内改变+π/2或-π/2,可以在保持窄带宽的前提下,实现数字错误位保持减少的数字信号的发射和接收。这种情况下的FSK方法特别称为MSK(最小频移键控,Minimum Shift Keying)方法。
图4是说明FSK调制器输出的连续波的频谱的图。在图4中,水平轴表示频率,垂直轴表示每个频率分量的电功率,图4所示的曲线图显示了FSK信号(FSK连续波)频谱的变化,其中通过频移键控的方法,根据调制度,该FSK信号的频谱相对于频率为f0载波(RF波)移动f1或f2。这里,FSK的调制度m定义如下:
通过上述公式表示的调制度m是频率差(f1和f2的差)和比特率fb(bit/s)的比,并与模拟FM调制的调制度β的意义相似。增大调制度m对应于设置较大的频率差,使得频谱扩展。因此希望在考虑例如超声波探头的频带等条件下设置调制度m。
根据本实施例,对FSK调制器(图1中以附图标记20表示)形成的发射信号进行延迟处理,从而形成参考信号,之后在接收混频器(图1中以附图标记30表示)中使用该参考信号对接收信号进行混频处理(即,参考信号与接收信号相乘)。在该混频处理中,来自对应于经过延迟处理的参考信号相位的深度(即目标位置的深度)的接收信号与参考信号的相关性增加到最大,而来自除目标位置的深度以外的深度的接收信号和参考信号之间的相关性大大减小。由此,根据本实施例的超声波诊断装置具有位置选择性功能。
图5到图7是说明根据本实施例的位置选择性的图,每幅图图解说明了每个来自目标位置的接收信号(接收波)和参考信号(参考波)的相位的时间变化,以及来自目标位置的接收信号和参考信号之间的相位差的时间变化。在图5到图7的每一幅图中,相位是在伪随机信号的1位周期T期间的连续波相位的变化。进一步地,每幅图也图解说明了在的情况下,相位差的时间变化。
乘法器(即图1中的接收混频器30)根据接收信号和参考信号之间的相位差输出结果(例如电压)。例如,当上述信号之间的相位差是π/2,3π/2,-π/2或-3π/2时,乘法器的输出值为“0”,当上述信号之间的相位差为0或2π时,乘法器的输出值为“+1”,当上述信号之间的相位差为π或-π时,乘法器的输出值为“-1”。
这里,尽管图5到图7所示示例中使用了在7位周期中变化的伪随机模式,在该装置应用时,希望使用在100到200位周期中变化的伪随机模式。
图5图解说明了参考信号的相位与来自目标位置的接收信号的相位一致的情况。当参考信号的相位与来自目标位置的接收信号的相位一致时,上述信号之间的相位差总为0。假设当参考信号和接收信号之间的相位差为0时,乘法器(即图1中的接收混频器30)的输出为“+1”,当上述信号的相位彼此一致时,7位周期中可得到的乘法器输出的相加值为“+7”。
同时,图6和图7图解说明了参考信号的相位与来自目标位置的接收信号的相位不一致的情况。
具体地,图6图解说明了参考信号(粗体虚线)的相位相对于来自目标位置的接收信号(实线)的相位延迟T的情况。这种情况下上述信号之间的相位差如图6下半部分所示,其中在的情况下乘法器(即图1中的接收混频器30)的输出总是“0”,且7位周期中可得到的乘法器输出的相加值也为“0”。
图7图解说明了参考信号(粗体虚线)的相位相对于来自目标位置的接收信号(实线)的相位延迟2T的情况。这种情况下上述信号之间的相位差如图7下半部分所示,其中在的情况下乘法器(即图1中的接收混频器30)的输出是随机重复的“+1”或“-1”,且7位周期中可得到的乘法器输出的相加值为“+1”。
如上文所述,当参考信号的相位与来自目标位置的接收信号的相位一致时,乘法器输出的相加值为“+7”,该值为极大值,而当参考信号的相位和接收信号的相位彼此不一致时,乘法器输出的相加值为“0”或“+1”,该值大大减小了。
图8是图解说明参考信号的相位和乘法器输出之间的关系的图。图8所示的曲线图对应于图5到7所示的示例。在图8中,曲线图的垂直轴表示7位周期中可得到的乘法器输出的相加值,水平轴表示相对于接收信号的参考信号(参考波)的相位,即接收信号和参考信号的相位差。这里,图8的水平轴给出的是当作为参考时,相位的相对值(即,当时,假设为“1”)。
当参考信号和接收信号之间的相位差为0时,如上文参考图5所述,乘法器输出的相加值为“+7”。因此,在图8中,参考波相位为“0”时得到的乘法器输出为“+7”。进一步地,如上文参考图6和图7所述,当相位差为时,乘法器输出的相加值为“0”,当相位差为时,乘法器输出的相加值为“+1”。因此,在图8中,当参考波相位为1时,所得到的乘法器输出值为“0”,当参考波相位为2时,所得到的乘法器输出值为“+1”。
图8也图解说明了通过进一步将参考信号的相位变为3、4、5、6、7等等,所得到的乘法器输出值。如图所示,当相位差为0时,乘法器输出值(7位周期内的相加值)为“+7”,该值与其它相位差的情况相比明显较大。进一步地,当参考波的相位差的相对值移动7后,由于将具有7位周期的伪随机信号移动了对应于一个周期的量,因此此时乘法器的输出与相位差为0时相同,也为明显的大值“+7”。
由上文可知,根据本实施例,通过调整参考信号的相位使其与到目标的延迟时间相对应,可以选择性地检测反射波功率以及只与目标深度有关的多普勒信息。同时,这种选择性随着伪随机模式长度的增加而更加敏锐。进一步地,随着模式长度的增加,该选择性逼近与脉冲波多普勒(PW多普勒)的距离选通相似的特性。换言之,可以在保持连续波的优点的同时,获得脉冲波的优点。
图9是说明根据本实施例的解调信号的图,并图解说明了通过混频器处理获得的解调信号的频谱。图9所示的解调信号对应于当参考信号与接收信号的相关度最大时,参考信号与接收信号的相乘结果。换言之,图9所示的解调信号是来自目标位置的接收信号,和相位与目标位置的深度匹配的参考信号的相乘结果。
图9所示的解调信号包含直流(低频)信号分量,和频率为RF信号频率f0两倍的谐波分量。这里,如图2(II)所示,经过FSK调制的连续波的频率是例如频率f1和频率f2的随机重复。当频率f1和频率f2之间的差比频率f0小时,在直流(DC)附近,以重叠方式检测到对应于频率f1的多普勒频谱和对应于频率f2的多普勒频谱,从而得到如图9所示的频谱。
出现的多普勒信号依附于直流信号分量以及频率为RF信号频率f0两倍的谐波分量。这里,由于LPF(图1中以附图标记36和38表示)除去了谐波分量,因此只提取直流(低频)信号分量。只有图9中所示的直流信号分量的频谱被发送到FFT电路(图1中以附图标记40和42表示)。之后,在多普勒信息分析单元(图1中以附图标记44表示)中,从图9所示的直流信号分量的频谱中提取多普勒信号,以及基于多普勒频移量等,计算例如目标位置存在的血液的流速。由于接收混频器(图1中以附图标记30表示)中进行正交检波,因此也可以确定流速的极性。同时,通过从直流信号分量的频谱中提取杂波信号,也可以得到例如存在于目标位置的血管壁的位置。
这里,超声波在活体内传播时的衰减特性称为频率相关衰减(FDA),已知其基本上与频率以及在活体内的软组织中的传播距离成正比。FDA影响解调信号频谱的形状。更具体而言,与从身体表面附近获得的解调信号的频谱相比,从较深的位置获得的解调信号的频谱经受较大的衰减。因此,可以对FFT电路(图1中用附图标记40和42表示)输出的频谱进行针对活体内的频率相关衰减的补偿处理。例如,从较深的位置获得的信号被较大程度地放大,从而抵消衰减量随位置(深度)的变化。还可以对接收混频器(图1中用附图标记30表示)的输出进行针对频率相关衰减的补偿处理。
图10图解说明了模式发生器(图1中用附图标记24表示)的具体示例。在图10所示的示例电路结构中,将晶体振荡器(OSC)输出的载波(RF信号)变换为矩形信号并且n分频,将所得的信号输入到移位寄存器。此外,将从移位寄存器的中间级选取的输出n1至n3和移位寄存器的最后一级的输出进行异或,从而可以得到PN模式。当移位寄存器具有n级时,PN模式的长度为(2n-1)位。例如,当n=10时,可以简单地配置1023位的PN模式。在这种情况下,接收信号和参考信号的相位彼此完全相同时得到的乘法器输出是1位时的1023倍,而其他位置的输出最多是1位时的数倍。因此可以显著增强选择能力。图10中的电路可以被配置为数字电路,在易于获得集成电路(IC)方面也是有利的。
基于模式发生器输出的周期信号序列,FSK调制器(图1中以附图标记20表示)通过进行频移键控(FSK)来形成连续波。FSK调制器可采用预先准备好的频率为f1的载波和频率为f2的载波,并根据例如图10所示的电路得到的PN模式适当地选择频率为f1的载波或频率为f2的载波,从而在时间轴上将所选择的载波彼此相连,由此形成FSK连续波。
图11是说明延迟电路(图1中用附图标记26I和26Q表示)和接收混频器(图1中用附图标记30表示)的变化示例的图。在图11所示的电路中,通过n级移位寄存器对参考信号进行延迟,从移位寄存器的各抽头n1至nn同时输出延迟时间互不相等的多个延迟参考信号。之后,在与移位寄存器的各抽头n1至nn对应的n个混频器中的每一个中,将接收信号与每个延迟参考信号相乘,从而同时(并行)输出与多个延迟参考信号对应的多个相乘结果。通过预先将在活体内沿深度方向排列的多个目标位置与移位寄存器的各抽头n1至nn关联,可以同时获得多个目标位置的多普勒信息。例如,可以配置成可同时检测超声波束上所有位置的组织信息。
虽然已对本发明的优选实施例及一些变化示例进行了描述,但它们仅仅是各方面的例证性示例,因而并不限制本发明的范围。
例如,虽然在上述实施例中,使用频移键控(FSK)对连续波进行数字调制,对于本领域技术人员显而易见的是,也可以使用相移键控(PSK)、幅移键控(ASK)等等来代替FSK作为数字调制方法。这里,也可以将数字调制的连续波的数据存储在存储器等等中,并基于从存储器读取的数据来产生连续波。
Claims (15)
1、一种超声波诊断装置,包括:
发射信号处理单元,其输出发射信号,所述发射信号为基于周期信号序列进行了数字调制的连续波;
发射/接收单元,其向活体发射与发射信号对应的发射波,并从活体接收与发射波关联的接收波,从而获得接收信号;
接收信号处理单元,其通过使用基于发射信号形成的参考信号,对接收信号进行解调处理,从而获得解调信号;以及
活体内信息提取单元,其从解调信号中提取活体内信息,
其中
选择性地提取活体内目标位置的活体内信息。
2、如权利要求1所述的超声波诊断装置,其中
发射信号处理单元输出的发射信号是通过频移键控形成的连续波。
3、如权利要求2所述的超声波诊断装置,其中
发射信号处理单元输出的发射信号是通过基于周期信号序列的频移键控的方式改变频率而形成的连续波。
4、如权利要求3所述的超声波诊断装置,其中
发射信号处理单元输出的发射信号是通过在连续变化相位的同时改变频率而形成的连续波。
5、如权利要求4所述的超声波诊断装置,其中
发射信号处理单元输出的发射信号是通过在周期信号序列的每个信号的周期T内,对于每个信号以+π/2或-π/2来改变相位而形成的连续波。
6、如权利要求1所述的超声波诊断装置,其中
接收信号处理单元通过使用与发射信号的波形一致的参考信号对接收信号进行解调处理。
7、如权利要求1所述的超声波诊断装置,其中
通过对从活体内目标位置获得的接收信号的周期信号序列与被调整的参考信号的周期信号序列之间的相关性进行解调处理,选择性地提取目标位置的活体内信息。
8、如权利要求7所述的超声波诊断装置,其中
根据活体内目标位置的深度进行延迟处理以调整接收信号与参考信号之间的延迟关系,从而调整相关性。
9、如权利要求8所述的超声波诊断装置,其中
通过延迟处理,增强接收信号的周期信号序列与参考信号的周期信号序列之间的相关性。
10、如权利要求9所述的超声波诊断装置,其中
通过将参考信号以根据目标位置的深度而定的延迟量进行延迟,使从目标位置获得的接收信号的信号序列模式与参考信号的信号序列模式相互匹配。
11、如权利要求1所述的超声波诊断装置,其中
发射信号处理单元输出的发射信号是通过基于周期信号序列的频移键控的方式改变频率而形成的连续波,以及
接收信号处理单元通过使用与发射信号波形一致的参考信号对接收信号进行解调处理。
12、如权利要求11所述的超声波诊断装置,其中
通过对从活体内目标位置获得的接收信号的周期信号序列与被调整的参考信号的周期信号序列之间的相关性进行的解调处理,选择性地提取目标位置的多普勒信息。
13、如权利要求12所述的超声波诊断装置,其中
通过将参考信号以根据目标位置的深度而定的延迟量进行延迟,使从目标位置获得的接收信号的信号序列模式与参考信号的信号序列模式相互匹配。
14、如权利要求1所述的超声波诊断装置,其中
活体内信息提取单元提取与解调信号中所包含的直流信号分量对应的多普勒信号分量,作为活体内信息。
15、如权利要求1所述的超声波诊断装置,进一步包括
延迟处理单元基于多个不同的延迟量对参考信号进行延迟,从而形成与多个目标位置对应的多个延迟参考信号;
其中
接收信号处理单元通过使用多个参考信号对接收信号进行解调处理,从而形成与多个目标位置对应的多个解调信号;以及
活体内信息提取单元基于多个解调信号提取在活体内沿深度方向排列的多个目标位置的活体内信息。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008-156770 | 2008-06-16 | ||
JP2008156770A JP5654198B2 (ja) | 2008-06-16 | 2008-06-16 | 超音波診断装置 |
JP2008156770 | 2008-06-16 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101606851A true CN101606851A (zh) | 2009-12-23 |
CN101606851B CN101606851B (zh) | 2012-10-03 |
Family
ID=41059705
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910147985XA Expired - Fee Related CN101606851B (zh) | 2008-06-16 | 2009-06-12 | 超声波诊断装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US8343055B2 (zh) |
EP (1) | EP2135556B1 (zh) |
JP (1) | JP5654198B2 (zh) |
CN (1) | CN101606851B (zh) |
AT (1) | ATE536806T1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103549975A (zh) * | 2013-10-14 | 2014-02-05 | 华南理工大学 | 基于跳频超声信号分离的结石检测装置及其方法 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8100832B2 (en) * | 2007-04-27 | 2012-01-24 | Hitachi Aloka Medical, Ltd. | Ultrasound diagnostic apparatus |
JP5629052B2 (ja) * | 2008-06-03 | 2014-11-19 | 日立アロカメディカル株式会社 | 超音波診断装置 |
JP5654198B2 (ja) | 2008-06-16 | 2015-01-14 | 日立アロカメディカル株式会社 | 超音波診断装置 |
JP6116853B2 (ja) * | 2011-11-30 | 2017-04-19 | 東芝メディカルシステムズ株式会社 | 超音波診断装置及び画像処理方法 |
US8737555B2 (en) * | 2011-12-22 | 2014-05-27 | Landis+Gyr Technologies, Llc | Digital signal processing for PLC communications having communication frequencies |
US10270491B2 (en) | 2017-08-31 | 2019-04-23 | Landis+Gyr Technologies, Llc | Power-line communication systems AMD methods having location-extendable collector for end-point data |
US10340980B1 (en) | 2018-05-07 | 2019-07-02 | Landis+Gyr Technologies, Llc | Time synchronization apparatuses and methods for power-distribution systems and the like |
US11255948B1 (en) * | 2019-06-08 | 2022-02-22 | Robert Dickerman | Determining frequency spectra and polarity of frequency for imbalanced quadrature signals including for Doppler system target speed and direction |
FR3101955B1 (fr) * | 2019-10-09 | 2021-10-29 | Bh Tech | Dispositif et procede de mesure de distance par ultrasons |
Family Cites Families (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4167879A (en) * | 1978-06-06 | 1979-09-18 | Panametrics, Inc. | Method and apparatus for examining a solid |
US4176351A (en) | 1978-08-18 | 1979-11-27 | Raytheon Company | Method of operating a continuous wave radar |
FR2447041A1 (fr) * | 1979-01-19 | 1980-08-14 | Inst Nat Sante Rech Med | Perfectionnements aux velocimetres doppler a bruit pseudo-aleatoires |
US4578677A (en) * | 1983-09-23 | 1986-03-25 | The United States Of America As Represented By The Secretary Of The Navy | Range doppler coupling magnifier |
CN86203861U (zh) | 1986-05-31 | 1987-11-11 | 西安交通大学 | 一种扩展连续波多普勒血流仪功能的装置 |
US5022400A (en) * | 1988-02-08 | 1991-06-11 | Walters Glenn A | Large time bandwidth echographic signal processor |
US4918706A (en) | 1988-12-28 | 1990-04-17 | Sperry Marine Inc. | Spread spectrum long loop receiver |
US5224482A (en) * | 1991-04-08 | 1993-07-06 | Kabushiki Kaisha Toshiba | Ultrasound high velocity flow correlation measurement using coded pulses |
JPH04351093A (ja) * | 1991-05-28 | 1992-12-04 | Mitsubishi Heavy Ind Ltd | 超音波テレメーター |
JPH08204766A (ja) * | 1995-01-25 | 1996-08-09 | Katsuyoshi Azeyanagi | 時間制限波形の周波数分析検出方式 |
US6179781B1 (en) * | 1999-03-31 | 2001-01-30 | Acuson Corporation | Medical diagnostic ultrasound method and apparatus for improving doppler processing |
US6491631B2 (en) * | 2001-01-11 | 2002-12-10 | General Electric Company | Harmonic golay-coded excitation with differential pulsing for diagnostic ultrasound imaging |
JP2002136522A (ja) * | 2000-11-02 | 2002-05-14 | Japan Science & Technology Corp | 超音波測定装置 |
US7094204B2 (en) * | 2002-08-23 | 2006-08-22 | Siemens Medical Solutions Usa, Inc. | Coded excitation imaging for use with bipolar, unipolar and other waveforms |
US6953434B2 (en) * | 2002-09-24 | 2005-10-11 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus to enhance ultrasound contrast imaging using stepped-chirp waveforms |
US6960169B2 (en) * | 2003-05-19 | 2005-11-01 | Siemens Medical Solutions Usa, Inc. | Spread spectrum coding for ultrasound contrast agent imaging |
JP4688466B2 (ja) * | 2003-10-10 | 2011-05-25 | アロカ株式会社 | 超音波診断装置 |
EP1695665A4 (en) | 2003-12-02 | 2011-05-11 | Hitachi Medical Corp | ULTRASOUND DEVICE |
JP4490128B2 (ja) * | 2004-01-29 | 2010-06-23 | アロカ株式会社 | 超音波診断装置 |
JP4652731B2 (ja) | 2004-07-01 | 2011-03-16 | アロカ株式会社 | 超音波診断装置 |
US7338450B2 (en) * | 2004-08-27 | 2008-03-04 | General Electric Company | Method and apparatus for performing CW doppler ultrasound utilizing a 2D matrix array |
US8372010B2 (en) * | 2004-10-20 | 2013-02-12 | Kabushiki Kaisha Toshiba | Ultrasonic doppler diagnosis device |
JP2006288974A (ja) * | 2005-04-14 | 2006-10-26 | Kanazawa Inst Of Technology | 超音波診断装置 |
US7887487B2 (en) * | 2005-07-11 | 2011-02-15 | Siemens Medical Solutions Usa, Inc. | Ultrasound diagnostic flow imaging with coded excitation |
JP2007330541A (ja) | 2006-06-15 | 2007-12-27 | Kanazawa Inst Of Technology | 超音波診断装置 |
US8100832B2 (en) | 2007-04-27 | 2012-01-24 | Hitachi Aloka Medical, Ltd. | Ultrasound diagnostic apparatus |
JP5654198B2 (ja) | 2008-06-16 | 2015-01-14 | 日立アロカメディカル株式会社 | 超音波診断装置 |
-
2008
- 2008-06-16 JP JP2008156770A patent/JP5654198B2/ja not_active Expired - Fee Related
-
2009
- 2009-06-03 AT AT09007364T patent/ATE536806T1/de active
- 2009-06-03 EP EP09007364A patent/EP2135556B1/en not_active Not-in-force
- 2009-06-09 US US12/480,874 patent/US8343055B2/en not_active Expired - Fee Related
- 2009-06-12 CN CN200910147985XA patent/CN101606851B/zh not_active Expired - Fee Related
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103549975A (zh) * | 2013-10-14 | 2014-02-05 | 华南理工大学 | 基于跳频超声信号分离的结石检测装置及其方法 |
CN103549975B (zh) * | 2013-10-14 | 2016-08-10 | 华南理工大学 | 基于跳频超声信号分离的结石检测装置及其方法 |
Also Published As
Publication number | Publication date |
---|---|
ATE536806T1 (de) | 2011-12-15 |
EP2135556A1 (en) | 2009-12-23 |
US20090312636A1 (en) | 2009-12-17 |
EP2135556B1 (en) | 2011-12-14 |
CN101606851B (zh) | 2012-10-03 |
JP5654198B2 (ja) | 2015-01-14 |
US8343055B2 (en) | 2013-01-01 |
JP2009297351A (ja) | 2009-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101606851B (zh) | 超声波诊断装置 | |
CN101596115B (zh) | 超声波诊断装置 | |
US8100832B2 (en) | Ultrasound diagnostic apparatus | |
EP2383583A2 (en) | Acoustic Distance Measurement System Having Cross Talk Immunity | |
JP5459963B2 (ja) | 超音波診断装置 | |
JP2010279531A (ja) | 超音波診断装置 | |
CN111413416B (zh) | 基于本振信号的正交解调方法、装置、系统及超声系统 | |
JP2010110503A (ja) | 超音波診断装置 | |
JP5325503B2 (ja) | 超音波診断装置 | |
JP2010057812A (ja) | 超音波診断装置 | |
JP2010269004A (ja) | 超音波診断装置 | |
EP4040182A1 (en) | Radar signal processing device, radar device, radar signal processing method and non-transitory computer-readable medium | |
JP2009261749A (ja) | 超音波診断装置 | |
JP2011240006A (ja) | 超音波診断装置 | |
JP2009297444A (ja) | 超音波診断装置 | |
JP2011217898A (ja) | 超音波診断装置 | |
JP2011104241A (ja) | 超音波診断装置 | |
JP2011036599A (ja) | 超音波診断装置 | |
JP2010075289A (ja) | 超音波診断装置 | |
JP2011217897A (ja) | 超音波診断装置 | |
JP2010188115A (ja) | 超音波診断装置 | |
JP2009247365A (ja) | 超音波診断装置 | |
JP2010125246A (ja) | 超音波診断装置 | |
JP2010162082A (ja) | 超音波診断装置 | |
JP2011011041A (ja) | 超音波診断装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C53 | Correction of patent for invention or patent application | ||
CB02 | Change of applicant information |
Address after: Tokyo, Japan, Japan Applicant after: Hitachi Aloka Medical Ltd. Address before: Tokyo, Japan, Japan Applicant before: Arockar K. K. |
|
COR | Change of bibliographic data |
Free format text: CORRECT: APPLICANT; FROM: AROCKAR K.K. TO: HITACHI ALOKA MEDICAL LTD. |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121003 Termination date: 20150612 |
|
EXPY | Termination of patent right or utility model |