专利名称：基于正交多频超声稳态回波信号的成像方法与装置的制作方法
技术领域：
本发明涉及超声波成像技术，具体涉及利用正交频分复用(OFDM)调制方式构造相互 正交的连续多频超声信号，并通过测量该多频超声信号的稳态回波信号进行成像的方法与 装置。
背景技术：
目前通用的超声回波成像方法主要有两种基于单脉冲和基于单频连续波。基于单脉 冲的成像方法，是利用超声波在传播路径上遇到介质的不均匀界面时会发生反射的物理特 性来实现的。发送端向人体组织发送一个脉冲，由于人体组织和脏器具有不同的声阻抗， 使超声波在不均匀界面上反射回波，检测出这些回波信号，对其放大和处理，最后在显示 器上显示，即可得到脉冲回波的图像。由发送脉冲和接收的回波脉冲的时间差来估算人体 组织间的分布。
基于单频连续波的成像方法，是在发送端发送某一频率的连续波，通过比较发射超声 波和其回波之间的相位差来获得回波相对于发射波的延时，从而估算不均匀界面的距离。 但是，当有多个反射界面时，回波是由多个不同相位的连续波叠加形成的，会互相干扰， 使得难以精确估计每个界面的回波与发射波的相位差。基于单频连续波的成像方法多用于 利用多普勒频移来检测回波波源的速度分布。
基于脉冲的成像方法，因为每个脉冲是瞬时波，每个脉冲的宽度(持续时间)相当短， 致使接收回波脉冲的检测难度大，容易漏检；若增加瞬时脉冲的宽度，又会给接收回波的
装置增加定时(定位)的难度，定位的偏差会带来延时估计比较大的误差。同时，超声探 头不能在同一时间发射和接收超声波，只有当探头发射完脉冲以后，处于静止期间，才能
进行接收。假设有两个界面A和B相距很近，当一束脉冲射入时， 一部分被A面反射， 另一部分传到B面被反射，如果脉冲前沿经B反射到A时，A的反射脉冲后沿还没有离开， 这样经A反射的脉冲后沿就和经B反射的脉冲前沿衔接起来，重叠成一个脉冲，不能分辨 出两个脉冲，那么A、 B两点将不能被分辨。所以，只有当脉冲宽度小于往返两点之间所 需的时间时，回波显示器上才能显示出两个独立的回波信号。脉冲宽度愈小，纵向分辨力 愈高。但是縮短脉冲宽度又会减小发射超声的能量，影响灵敏度。用单频连续波来获得超 声波成像的方法，适合于只有单个反射界面的情况，因为此时只有一个回波，相位估算的 精度高。对于有多个反射界面的情况，得到的是多个反射回波的叠加，此时很难估算出每 个界面的反射回波。
由此可见，基于脉冲成像的方法会带来纵向分辨力和灵敏度之间的矛盾，而用单频连 续波来获得超声波成像的方法又很难估算出多个界面的反射回波。并且这两种方法都易受 噪声的干扰，而本发明基于正交多频超声稳态回波信号的成像方法与装置，可以很好地解
决上述问题。因为该方法接收的反射回波是一个多频的稳态连续波，由于发射的是多频调 制波，携带信息多，通过设计专门的算法，能够分解出每个界面的反射回波，并且可消除 噪声的影响；由于是连续波，能量稳定，不需要定时，因而检测容易，不存在定时上的误 差，同时还能够兼顾纵向分辨力和灵敏度。
发明内容
本发明的目的在于克服现有技术存在的测量误差大、难以兼顾纵向分辨力和灵敏度等 不足，提供基于正交多频超声稳态回波信号的成像方法与装置，该方法利用正交频分复用 (OFDM)调制方式来构造相互正交的多频连续基带信号，通过检测稳态的回波信号，计 算回波信号相对于发射信号的延时来成像。本发明装置计算精度高，抗干扰性好，抗噪能 力强，安装简单，使用方便，且成本低廉。本发明通过如下技术方案实现
基于正交多频超声稳态回波信号的成像方法，包括如下步骤利用正交频分复用
(OFDM)调制方式构造相互正交的多频连续基带信号，并用超声载波(载波频率为,)
将该多频连续基带信号调制为多频超声连续发射信号；超声波探头将该多频超声连续发射 信号发送出去，相当于发送多路单频超声波连续信号，每一路信号经过待成像对象的每一 个不均匀界面都会有反射；超声探头检测反射回来的稳态回波信号，该回波信号相当于多 路反射波的叠加信号；该叠加信号经解调后得到接收的多频连续基带信号，对该多频连续 基带信号进行正交分解，利用每路信号互为正交的特性分解出每路接收信号，列出关于每 个不均匀界面反射波的幅度和相位的方程组；通过解方程组求出不均匀反射界面反射波的 幅度和相位，计算每个界面反射回波的幅度和相对于发射波的延时(位置)来成像。
用于上述方法的装置包括超声波探头、控制单元及分别与控制单元连接的操作单元、 收发单元、处理单元和显示单元，超声波探头与收发单元相连，处理单元还分别与收发单
元和显示单元连接；
操作单元为人机交互模块，用于操作者选择测量信号的超声波频率和对显示单元的显
示方式进行选择；
控制单元根据操作单元发出的操作指令来控制收发单元、处理单元和显示单元； 收发单元通过将驱动信号传给超声波探头来发送一个超声波，超声波探头也将收到的
回波信号传送给收发单元；
处理单元用于产生多频正交连续信号和调制解调用的超声载波，并将解调出的多频连
续基带信号经过A/D采样、量化成数字信号，然后做正交分解，列出关于每个界面反射波
的幅度和相位的方程组，再进行计算处理；
显示单元对处理单元处理后的结果根据操作者选择的显示方式来显示成像，所述显示
方式包括幅度调制和辉度调制方式；
上述装置中，所述处理单元包括数字信号处理器、调制解调器、振荡器、若干分频器 和A/D转换器，处理单元完成信号产生和处理的工作根据控制单元提供的信息产生对应
的多频连续基带信号，基带信号经过调制成为超声波信号并传给收发单元；对收发单元传 送过来的回波信号进行解调；对解调后的信号进行处理，计算每个界面反射回波的幅度和 相对于发射波的延时(位置)。
上述装置中，所述处理单元包括三个分频器分频器之一控制输出多声频信号频率， 分频器之二控制输出超声载频信号频率，分频器之三控制输出A/D采样时钟信号频率。多 频连续基带信号经超声频载波信号调制后通过有线或无线方式传送到收发单元，再由收发 单元发送给超声波探头。
本发明的工作过程为通过操作单元发出操作指令，然后控制单元发出相应的控制给处 理单元，处理单元立即给本次操作产生特定的多频基带信号和调制解调超声频信号，即配 置分频器1 (控制输出多声频信号频率)、分频器2 (控制输出超声载频信号频率)和分频 器3 (控制输出A/D采样时钟信号频率)的控制参数。随后多频连续基带信号经超声频载 波信号调制后通过有线或无线方式传送到收发单元，再由收发单元发送给探头。
然后收发单元等待控制单元发出接收回波的命令，收发单元接收到稳态的回波后给处理 单元，处理单元对稳态回波先解调，然后做正交分解，列出关于每个界面反射波的幅度和 相位的方程组。
最后通过解方程组求出非均匀反射界面反射波的幅度和相位。
若采用振幅调制的方式(A型超声诊断仪所用的方式)来显示成像，则横坐标为被测组 织的深度和不均匀界面之间的距离间隔(计算出来的相位转换成延时，再转换成距离)，纵 坐标表示回波的幅度，移动探头就可以获得不同位置处的回波信息。若采用辉度调制的方 式(B型超声诊断仪所用的方式)来显示成像则需要快速扫描逐次获得不同位置的深度方 向所有界面的反射回波。 一帧扫描完成后，便可以获得一幅扫描方向决定的垂直平面二维 超声断层图像。显示的成像根据回波的强弱(即是振幅的大小)来控制荧光屏上的光点辉 度，根据不均匀界面的位置(由相位求得)来控制光点的位置。
与现有技术相比，本发明具有以下优点-
1、 本发明测量精度高。根据超声协会所定标准，超声波在人体软组织中平均声速为 1540m/s,基于超声波脉冲成像方法中脉冲宽度一般是几ns，假设lps，则接收回波脉冲定 时最大误差为lns，带来的测量距离误差是1540m/s^ns^.54mm，所以其测量精度1.54mm。 而本发明的测量精度可以提高很多，假设控制单元的A/D采样率为lOMHz (lxl06Hz)， 测量误差为一个样点，那么测量精度为0.154 mm。而且还可以通过设计数字信号处理算法 进一步提高相位的估计精度，从而提高测量精度。
2、 本发明抗干扰性好，利用每路信号相互正交的特性来消除相互之间的干扰。而且因 为噪声和每路信号不相关，所以还能消除噪声的影响。
3、 由于本发明接收到的反射回波是一个稳态连续波，能量稳定，所以检测容易，不需 要定时，因而不存在定时上的误差，既可以提高纵向分辨力，又可以保证测量的灵敏度。
4、 由于本发明采用单探头，所以对换能器的要求降低，使得成本降低。现代社会人们 对健康、疾病越来越关注，成本的降低使得该装置很适合农村、社区和家庭医疗使用。
5、 由于本发明的成像精度不再依赖于测量的精度，而是与A/D采样率、计算精度相关， 所以一方面可以提高采样率来提高测量精度，另外一方面可以把测得的稳态连续回波存储 或者通过网络发送到大型的处理机来处理计算，通过增加计算量来提高幅度和相位的估计 精度，得到更为精确的成像结果，诊断出细微的病理变换，提高诊断水平。
6、 随着现代处理器(处理机)计算处理能力的大大提高，本发明的信号处理过程中求 解(公式IO)或者(公式ll)的幅度和相位不再是难题，所以该技术是可行的，且适用性强。
7、 本发明装置简单，安装使用方便，可以应用于探头、收发单元与处理单元有线连接 和无线连接的场合。
图1为本发明装置在实施方式中的组成框图2为图1所示装置的处理单元采用有线工作方式时的组成框图； 图3图1所示装置的处理单元采用无线工作方式时的组成框图 图4a为实施方式中多频正交连续信号调制前的频谱示意图； 图4b为实施方式中多频正交连续信号调制后的频谱示意图； 图5为本发明装置在实施方式中的工作流程示意图6为图5中处理单元求解每个非均匀界面反射波的幅度和相位的具体流程；
具体实施例方式
下面结合附图对本发明的具体实施方式
作进一步说明。
如图1所示，本发明装置包括超声波探头、操作单元、控制单元、收发单元、处理单 元和显示单元。如图2所示，处理单元包括数字信号处理器(DSP或者单片处理器ARM等)、 调制解调器、振荡器、若干分频器和模数转化器(A/D)。振荡器使用20M晶体振荡器，分 频器用锁相环电路实现，调制器和解调器用乘法电路实现，模数转化器和数字信号处理器 选用一块TI的DSP芯片TMS 320C2812实现。
超声波探头即为一个超声波转换器，由一种压电材料，即压电陶瓷如钛酸钡、钛酸铅、 锆钛酸铅(PZT)等组成。超声波探头用来与成像对象接触。超声波探头与收发单元相连， 收发单元通过将驱动信号传给超声波探头来发送超声波，超声波探头也将收到的回波信号 传送给收发单元。
操作单元也即是人机交互接口，由操作者根据实际需要选择操作，如选择检査人体的 不同部位或者组织。因为超声波测量不同组织所需要的超声波频率是不一样的，通常用于 眼、乳腺等浅表部位疾病探査的超声频率为十几至几十兆赫兹，而用于腹部疾病探査频率 为5兆赫兹左右。
控制单元根据操作单元发出的操作指令来控制收发单元、处理单元和显示单元。
处理单元结构如图2所示，包括数字信号处理器(DSP或者单片处理器ARM等)、调制 解调器、振荡器、若干分频器和模数转化器(A/D)。处理单元主要完成信号产生和处理的工 作1)根据控制单元提供的信息产生对应的多频连续基带信号，基带信号经过调制成为超 声波信号并传给收发单元；2)对收发单元传送过来的回波信号进行解调；3)对解调后的 信号进行处理，计算每个界面反射回波的幅度和相对于发射波的延时(位置)。
显示单元对处理单元处理后的结果根据操作者选择的显示方式来直观地显示，提供给
操作者观察，该实施例可以选择用幅度调制或辉度调制的方式来显示成像。显示设备可选 用14寸液晶显示屏。
收发单元和控制单元、处理单元有两种连接方式有线连接和无线连接方式。在如图 3所示无线连接方式中，控制单元、处理单元通过射频收发模块来传送相应信息，收发单 元也应包括射频收发模块。
本发明利用OFDM调制方式来构造相互正交的连续多频超声信号，通过检测稳态的回 波信号，计算回波信号相对于发射信号的延时来成像，其原理如下
设发射的单频超声波连续信号为"W = ^eQS(2^^)，该信号透过人体组织在第一非均匀介
质界面反射回来的回波信号经过测量为浙a)"cos(2Wi)。其中,。是发射信号频率，伊'
是相位差。根据数字信号处理理论，可以计算出发射和接收信号的相位差"为
A =maxarg{〈s(0，"G,A》} (公式1)
〈"〉为两个信号的互相关操作。那么超声波发射端即探头和人体组织第一非均匀介质界面 的距离为
/,=^^rv (公式2)
其中v为超声波在人体组织内的传播速度。
当有多个非均匀介质界面，根据超声波的透射和反射原理，设每个界面反射回来的回波
信号为^,W"cos(2Wi)，若有N个界面则接收到的回波信号为
々,anw) = |;4COS(2;r/。,-a) (公式3)
这一个方程无法解出N个界面对应的相位差A,为了解出N个界面对应的相位差
A，…^V和幅度4…^w，至少需要构造M22W个方程。我们利用OFDM调制方式的原理 来构造M = 个相互正交的多频超声连续发射信号，这样来构造M个方程。
指数信号y;w^'w和/^)^'w在[o，r]区间正交，令/。/"a/，则正交性需要满足 & = V2' ^ = V2",'-,2)'^ = s—yr) e,攀=。(公式4)
当A/r-"， ^ = ^就可以得到&=0。 OFDM调制中用到的是""。基于这个原理产生M个 相互正交的信号。设观察时间周期为r-lww， A/ = + = l^fe， M个正交信号的频率组 C/i,…/M)可以从集合F:("^/^里面选择(注意这里所选的频率要远远小于载波频率,)， n为正数。这样，这M个信号在[O，r]区间是相互正交的。M个正交信号的叠加即为需要发 送的多频连续基带信号,(^,..^) = |>。3(2;^,)，它可以看作是M路信号的发送，每一路
,=1
发送的是一个单频信号。
多频连续基带信号 ^ 产生后，选取合适的超声载波(载波频率为
,)对其进行调制。这相当于对M路的每一路信号进行调制，调制后的信号成为多频超声 连续发射信号，由收发单元发送出去。每一路信号经过每一个不均匀界面都会有反射。
假设第/路发送的基带信号是cos(2;r々)，其频率/=《=込/,然后假设超声波会遇到 N个不均匀的界面，这样就会收到N个反射波(解调后的)的叠加信号
《(C) = 4 cos(2;nyj (f — t, )) + cos(2;nyj (卜r2)) + + ^ cos(2;nyj (f — )) =4 cos(2;r々—) + ^2 cos(2;r 乂/ - /^ 2) +…+ Jw cos(2;r 乂f 一 )
(公式5)
其中^cos(2;ry;(卜r,》-4cos(2;ry;卜/^)是发射波发射出去后遇到的第一个不均匀界面的
反射波(解调后的)，4是幅度，r,是延迟时间，/^-2"/巧二2;ny^是相位差。以此类推，
^vCOS(2;ry;(f-^)) = 4^03(2^1/;/-/^^)是发射波发射出去后遇到的第N个不均匀界面的
反射波(解调后的)，4v是第N界面的幅度，^是第N界面反射波相对于发射波的延迟时
间，*w = 2;r= 2;ny;^是第N界面反射波相对于发射波的相位差。
收发单元接收到的回波是M路反射波的叠加，经解调之后得到多频连续基带信号 w,y;，…/M)。超声波在发送、传输、接收的时候会不可避免地混入噪声"(O，所以解调后
的多频连续基带信号为
<formula>formula see original document page 10</formula>先对此基带信号d(/， /,…/M)做理论上的分析。对d(/，,，…/M)在[o,r]上做正交分解，
利用每路信号互为正交的特性可以分解出每路的接收信号。(公式6)与第i路发送基带信号 cos(2;ry;/)在[o，r]区间上求互相关，则可以得到第i路的接收信号
<formula>formula see original document page 10</formula>(公式7)
因为每路信号是互为正交的，所以当"/时，f《",/》cos(2;r_/;0^ = 0。噪声和信号 是不相关的，因此f G)cos(2;ry;0^ = 0，所以(公式7)第三个等号成立。从这里可以看出，
正交信号抗干扰性好，并且还能够消除噪声。把(公式5)代入(公式7)，可以得到g,々,y;)的
展开表达式
<formula>formula see original document page 10</formula>(公式8)
由(公式7)可知，d(/，yj，…/M)cos(2;r々)通过积分得到一个确定的值，即
<formula>formula see original document page 10</formula>(公式9)
因为多频连续基带信号调制到高频载波处是一个窄带信号，即,》4/~， 虽
然每路信号的频率不同，但是都接近于载波频率，所以每路超声波的传播速度相同，每个
不均匀界面反射波的延迟时间r相同，幅度A也相同。由此M路可以得到M个方程
<formula>formula see original document page 10</formula>(公式10)
牟cos(M伊j) + 4 cos(M伊2) +…+ A cos(Mpw)=《 若不均匀界面的个数N是已知的，那就构造M-2N个相互正交的多频超声连续信号，
这样来构造2N个方程。由2N个方程就可以解出2N个参数。
(公式IO)方程组是一个非线性方程组，其求解方法有很多，现简介两种求解方法。 求解方法一采取穷举方法来求解，这种方法比较适用于反射界面不多的情形。由于超
声波的传播有衰减，每个界面的幅度满足-1<4，4，一Ar<l，且传播距离越远，相位差越 大，因此每个界面的反射回波相位差满足0< a <…<& <;r 。
举例把(0,;r)以;r/100的间隔等分，^,^，…，^在(0，;r)区间100个点值里面穷举， 且满足0<^<%<...<^^<;2"。每一轮计算过程1)炉p^,…，^在(0,;r)区间IOO个点值
中取满足要求的值；2)把^，^，…，^的值代入(公式10)，得到超定线性方程组；3)对超
定方程组求解，判断该超定方程组是否存在数学解。先去掉超定方程组中一^相关的方程， 然后求解4，4，…^^。若该超定方程组不存在数学解，则本轮计算结束；4)判断所求的幅
度是否满足物理要求，即是否满足-1<4，4，一^^<1，若不满足，则该轮的^,^,…,^v取
值不正确，重新选择新的值来进行下一轮的计算过程。
(公式IO)方程组的求解方法二把非线性方程转换成切比雪夫多项式线性方程组。0
次切比雪夫多项式rQ(x) = l ，设1次切比雪夫多项式
7^)-cos(a)^,7]02)二cos02)-;)C2…7](Xw):cos(^)-^ ， 贝'j
a = cos(;c,), % = wc cos(x2)…= "atcos(;^)。由切比雪夫多项式性质知，M次切比雪
夫多项式4(;0 = cos(M^),rM(jc2) = cos(M%)'..rM(;Cw) = cos(M^)，则可以把(公式10)转 换成为切比雪夫多项式线性方程组
4r2(Xl) + 4r2(X2) — .. + 4vr2(xw) = 52' (公式U)
.47^ ")+471/ (a)+…+4v^ Ow) =
可以采用牛顿法、赛德尔等成熟的数值计算方法解出4，4，…4^， x,，a，…^，进而解出
w，^2，…，^。
当不均匀反射界面的个数N是未知的，则根据实际情况估计人体组织的不均匀反射界 面个数来对N赋一个初始值L，该值应该大于实际的N值，即是Z〉iV。通过构造M-2L
个相互正交的多频超声连续信号，构造出M个方程。此时M(方程个数)〉2N(未知数个数)， 对该超定方程组求解，首先需要去掉一些冗余的方程，即是把相关的方程去掉，从而得到 不均匀反射界面的个数N和剩余的2N个方程。最后对剩余的2N个方程求解，求解的方法 和N已知时相同。
通过解方程组求出的不均匀界面反射波的幅度和相位。由(公式2)相位和位置的关系得 到不均匀界面间的距离分布。若采用振幅调制的方式(A型超声诊断仪所用的方式)来显 示成像，则横坐标为被测组织的深度和不均匀界面之间的距离间隔(计算出来的相位转换 成延时，再转换成距离)，纵坐标表示回波的幅度，移动探头就可以获得不同位置处的回波 信息。若采用辉度调制的方式(B型超声诊断仪所用的方式)来显示成像则需要快速扫描 逐次获得不同位置的深度方向所有界面的反射回波。 一帧扫描完成后，便可以获得一幅扫 描方向决定的垂直平面二维超声断层图像。显示的成像根据回波的强弱(即是振幅的大小) 来控制荧光屏上的光点辉度，根据不均匀界面的位置(由相位求得)来控制光点的位置。 举例探头所测的每个探测点的回波的观测时间设为r-l/7W，探测点与探测点之间的观测
时间保护间隙为2;=0.25; "假设每帧成像要显示100个探测点的信息，则需要(1+0.25)
"00-125ms的时间，所以单个探头的帧速为8帧/s,可以达到低速的成像水平。随着探头 个数(探头阵元个数)的增加，成像的帧速成线性增加，如只需要10个阵元，就能达到 80帧/s的高速成像水平。
在如图3所示无线连接方式中，处理单元和收发单元包含同有线连接方式相同的电路和 器件，并且还包括工作在2.40的射频收发模块。
本发明的工作步骤如图5所示，具体如下
步骤1检测操作单元是否有操作请求，有操作请求则立即传送到控制单元。 应用举例启动仪器，选择超声波要检査人体什么部位或者什么组织，是眼部浅表组
织，还是腹部的组织等。
步骤2控制单元接收到请求信号后，发出相应的控制信息给处理单元、收发单元和显
示单元。
控制单元控制收发单元什么时候发送超声波信号，什么时候接收；控制单元控制显示 单元的显示方式，可以选择以辉度调制或者幅度调制的方式来显示。根据操作请求附加显 示病人的相关信息，回波图像局部放大等。控制单元给处理单元的信息也主要是根据操作 请求来选择相应的信息，如操作单元选择检査浅表部分，则选择超声频率为十几、几十MHz， 若选择深层部位，则选择超声频率为5MHz左右。还要确定或者初步设置这种检査的非均 匀界面的最大个数N，选择2N个多频基带信号的频率和调制解调的超声载波频率。然后这 些相关信息一起给处理单元。
步骤3 处理单元配置各个分频器控制参数，产生多频正交连续基带信号和调制解调 用的超声载波。多频正交连续信号的频谱示意图如图4a所示。
本发明可以根据不同情况灵活选用2N个频率来产生多频正交连续基带信号，即从频率 集合F:^A/fen为正整数)里面任意选择2N个不同的频率，然后根据已经确定的载波频
率来产生载波，用来调制和解调。
步骤4 多频正交连续基带信号通过超声频载波调制后通过有线或无线方式发送到收 发单元，收发单元再通过与其相连的探头发射出去。调制后的多频超声连续发射信号的频 谱示意图如图4b所示。
步骤5 收发单元收到稳态的回波后，传送到处理单元，处理单元解调出多频基带信 号，经A/D采样、量化后成为数字信号。
解调出的多频连续基带信号经过A/D采样、量化成数字信号，送到数字信号处理单元 (DSP)计算处理，A/D的采样频率由DSP控制。'
步骤6 数字信号处理计算过程，主要求解(公式10)或者(公式11)，得到每个非 均匀界面反射波的幅度和相位，本实施例采用穷举的方法来求解方程组(公式10)，具体 流程参见图6。
数字信号处理单元(DSP)对A/D转化后的数字信号做正交分解，列出2N个方程。若 非均匀界面的个数N是已知的，则用穷举的方法求解方程个数为2N的方程组，得到N个 非均匀界面的幅度和相位。若非均匀界面的个数N是未知的，则需要先把线性相关的方程 组去掉，得到非均匀界面的个数M和剩余的2M个方程，然后令N=M，求解方程个数为 2N的方程组，得到N个非均匀界面的幅度和相位。穷举的时候根据每个非均匀界面反射波
的幅度和相位的限定条件-1<4，4，一^<1和0<^<%< <^<"来縮短穷举的计算时间。
步骤7 对处理单元得到的结果非均匀界面的个数、每个非均匀界面反射波的幅度、 相位信息(延时或者位置信息)以直观的形式在显示单元显示，方便和用户的交互。该实 施例可以对每个界面反射的回波以幅度调制或者辉度调制的方式来显示。
权利要求
1、基于正交多频超声稳态回波信号的成像方法，其特征在于包括如下步骤利用正交频分复用调制方式构造相互正交的多频连续基带信号，并用超声载波将该多频连续基带信号调制为多频超声连续发射信号；超声波探头将该多频超声连续发射信号发送出去，相当于发送多路单频超声波连续信号，每一路信号经过待成像对象的每一个不均匀界面都会有反射；超声探头检测反射回来的稳态回波信号，该回波信号相当于多路反射波的叠加信号；该叠加信号经解调后得到接收的多频连续基带信号，对该多频连续基带信号进行正交分解，利用每路信号互为正交的特性分解出每路接收信号，列出关于每个不均匀界面反射波的幅度和相位的方程组；通过解方程组求出不均匀反射界面反射波的幅度和相位，计算每个界面反射回波的幅度和相对于入射界面的位置来成像。
2、 根据权利要求1所述的基于正交多频超声稳态回波信号的成像方法，其特征在于采 用幅度调制或辉度调制方式来显示成像，若采用振幅调制的方式来显示成像，则横坐标为 被测组织的深度和不均匀界面之间的距离间隔，纵坐标表示回波的幅度，移动探头就可以 获得不同位置处的回波信息；若采用辉度调制的方式来显示成像，则需要快速扫描逐次获 得不同位置的深度方向所有界面的反射回波， 一帧扫描完成后，便可以获得一幅扫描方向 决定的垂直平面二维超声断层图像，显示的成像根据回波的振幅大小来控制荧光屏上的光 点辉度，根据不均匀界面的位置来控制光点的位置。
3、 根据权利要求1所述的基于正交多频超声稳态回波信号的成像方法，其特征在于当发射的单频超声波连续信号为s(O = ^COS(2;^(/)，该信号透过第一非均匀介质界面反射 回来的回波信号经过测量为d("^)=爿cos(2;r/(; —a)，其中/0是发射信号频率，^是相位差；当有N个非均匀介质界面，根据超声波的透射和反射原理，则接收到的回波信号为 ， A,…Pw) = X！ 4 cos(2W i) ，' =1为了解出N个界面对应的相位差^，…^和幅度4…4v，需要构造M22iV个方程；利用正交频分复用调制方式构造M个相互正交的多频超声连续基带信号，并构造M个方程。
4、 根据权利要求3所述的基于正交多频超声稳态回波信号的成像方法，其特征在于取 M-2N，利用正交频分复用调制方式构造M个正交信号，该M个正交信号叠加为需要发w送的多频连续基带信号W" 乂 ，…/a/ ) = S cos(2;r々)。,-=1
5、 一种用于实现权利要求1 4任一项所述方法的装置，其特征在于包括超声波探头、 控制单元及分别与控制单元连接的操作单元、收发单元、处理单元和显示单元，超声波探 头与收发单元相连，处理单元还分别与收发单元和显示单元连接。所述操作单元为人机交互模块，用于操作者选择测量信号的超声波频率和对显示单元 的显示方式进行选择；所述控制单元根据操作单元发出的操作指令来控制收发单元、处理单元和显示单元；所述收发单元通过将驱动信号传给超声波探头来发送一个超声波，超声波探头也将收 到的回波信号传送给收发单元；所述处理单元用于产生多频正交连续信号和调制解调用的超声载波，并将解调出的多 频连续基带信号经过A/D采样、量化成数字信号，然后做正交分解，列出关于每个界面反 射波的幅度和相位的方程组，再进行计算处理；所述显示单元对处理单元处理后的结果根据操作者选择的显示方式来显示成像，所述 显示方式包括幅度调制和辉度调制方式。
6、 根据权利要求5所述的装置，其特征在于所述处理单元包括数字信号处理器、调制 解调器、振荡器、若干分频器和A/D转换器，处理单元完成信号产生和处理的工作根据 控制单元提供的信息产生对应的多频基带信号，基带信号经过调制成为超声波信号并传给 收发单元；对收发单元传送过来的回波信号进行解调；对解调后的信号进行处理，计算每 个界面反射回波的幅度和相位，再把幅度和相位信息转换成为幅度和位置信息，在显示单 元按照所选的显示方式来显示成像。
7、 根据权利要求6所述的装置，其特征在于处理单元包括三个分频器，分频器之一控 制输出多声频信号频率，分频器之二控制输出超声载频信号频率，分频器之三控制输出A/D 采样时钟信号频率。
8、 根据权利要求7所述的装置，其特征在于多频连续基带信号经超声频载波信号调制 后通过有线或无线方式传送到收发单元，再由收发单元发送给超声波探头。
全文摘要
本发明公开了基于正交多频超声稳态回波信号的成像方法与装置，该方法利用正交频分复用调制方式构造相互正交的多频连续基带信号，经调制为多频超声连续发射信号后由超声波探头发送出去，超声探头检测反射回来的稳态回波信号，该信号相当于多路反射波的叠加信号；对经解调后得到接收的多频连续基带信号进行正交分解，利用每路信号互为正交的特性分解出每路接收信号，计算每个界面反射回波的幅度和相对于入射界面的位置来成像。该装置包括超声波探头、控制单元及分别与控制单元连接的操作单元、收发单元、处理单元和显示单元，超声波探头与收发单元相连，处理单元还分别与收发单元和显示单元连接；本发明成像精确、抗干扰性好；采用单探头，成本降低。
文档编号G01S15/89GK101357068SQ20081019851
公开日2009年2月4日 申请日期2008年9月12日 优先权日2008年9月12日
发明者宁更新, 燕 曹, 岗 韦 申请人:华南理工大学