一种骨骼超声系统中编码增强的聚焦超声骨组织微观结构检测方法
【专利摘要】本发明公开了一种骨骼超声系统中编码增强的聚焦超声骨组织微观结构检测方法，其步骤如下：一、采用超声背散射检测方法，依据Faran圆柱模型，将骨组织微观结构的检测等效为平均散射子间距的测量；二、采用高频单阵元聚焦换能器，利用互补Golay序列对激励超声换能器发射宽时信号，以发射信号作为匹配滤波函数对回波信号进行解码压缩，Golay序列对长度可调，利用解码压缩背散射信号估计平均散射子间距。该方法在利用超声背散射技术测量骨微结构MSS中提高了超声分辨率，同时减少了声衰减对测量结果的影响。本发明应用于生物医学检测、组织定征领域。
【专利说明】一种骨骼超声系统中编码增强的聚焦超声骨组织微观结构 检测方法
【技术领域】
[0001] 本发明涉及一种骨骼超声系统中骨组织微观结构检测方法，尤其涉及一种编码增 强的聚焦超声骨组织微观结构检测方法。
【背景技术】
[0002] 生物组织主要由周期性排布的胶原纤维组成，可被等效为周期性分布的散射子。 现有的很多模型以此为基础用来对生物组织进行表征，如用来对肝脏疾病、血管厚度和骨 小梁间距等进行定征。通过测量平均散射子间距（MSS)，来对组织疾病做出诊断。
[0003] 超声背散射技术已被证明能够有效地对MSS进行测量。生物组织背散射信号包括 来自规则散射兀的相干成分和来自弥散散射兀的非相干成分。对MSS的测量主要利用了其 中的相干散射兀成分，而弥散散射兀成分则对MSS的测量起干扰作用。现有的MSS测量技术 主要有周期图估计、倒谱估计和自相关谱估计等方法。但是，在对MSS进行估计的过程中， 受背散射信号序列短和噪声大的干扰，上述估计方法的鲁棒性较低。为此，一些方法如自回 归倒谱估计方法、周期性小波变换等方法被用来提高MSS估计精度。
[0004] 另外两个影响MSS估计精度的主要因素是分辨率和声衰减。目前在骨微结构MSS 测量中常用的超声频率为〇. 5MHz左右，超声波长大于6mm，在对骨微结构（0. 15?1. 2mm) 进行测量时，频率低，分辨率不足。此外，松质骨组织为多孔固体互联基质，超声在传播过程 中会由于多重散射而产生大量衰减，特别是频率越高，声衰减越大，使得背散射回波信号信 噪比较低，而现有技术没有有效地解决分辨率和声衰减对背散射信号MSS估计的影响。
【发明内容】
[0005] 为了解决现有检测技术频率低、分辨率不足、无法有效应对组织中大量声衰减及 噪声的问题，本发明提供了一种骨骼超声系统中编码增强的聚焦超声骨组织微观结构检测 方法。该方法在利用超声背散射技术测量骨微结构MSS中提高了超声分辨率，同时减少了 声衰减对测量结果的影响。
[0006] 本发明的目的是通过以下技术方案实现的：
[0007] -种骨骼超声系统中编码增强的聚焦超声骨组织微观结构检测方法，包括如下步 骤：
[0008] -、采用超声背散射检测方法，依据Faran圆柱模型，将骨组织微观结构的检测等 效为平均散射子间距的测量；
[0009] 二、采用高频单阵元聚焦换能器，利用互补Golay序列对激励超声换能器发射宽 时信号，以发射信号作为匹配滤波函数对回波信号进行解码压缩，Golay序列对长度可调， 利用解码压缩背散射信号估计平均散射子间距。
[0010] 所述步骤一中，超声背散射检测的测量方式为单一探头的背散射测量方式，测量 环境为水浸式测量。
[0011] 所述步骤一中，测量参数为散射子间距。根据骨微结构显微CT三维重建结果，将 骨微结构等效为间距排列的Faran圆柱散射子，利用该等效散射子间距来表征骨微结构。
[0012] 所述步骤一中，，采用的超声激励频率为0. 7?2MHz可调，激励的超声能量水平低 于 100mW/cm2。
[0013] 所述步骤二中，在用互补Golay序列对激励超声换能器时，根据被测骨组织厚度， 通过编程选择构造 4 位[(+1，+1，+1，-1) / (+1，+1，-1，+1) ]、8 位[(+1，+1，+1，-1，+1，+1，-1， +1) / (+1，+1，+1，-1，-1，-1，+1，-1)]或 10 位[(+1，+1，-1，+1，-1，+1，-1，-1，+1，+1)/ (+1， +1，-1，+1，+1，+1，+1，+1，-1，-1)]互补序列对[Golayl/Golay2]来激励超声换能器。超声 换能器采用聚焦换能器，对同一生物组织聚焦点分别采用不同Golay序列两次照射，发射 波持续时间Λ tin。< Lb/vb，Lb骨组织厚度，vb为骨组织声速。
[0014] 所述步骤二中，在激励超声换能器时，采用正弦信号作为基础信号调制Golay互 补序列。
[0015] 所述步骤二中，在对背散射回波信号进行解码压缩过程中，采集Golayl和Golay2 输入信号F inc;(t)F' in。⑴作为匹配滤波函数，则获得背散射信号为：
[0016] Fcep(t) = cov[Finc(t)Fbsc(t)]+cov[F/inc(t)F/bsc(t)],
[0017] 其中，Fbs。⑴F' bs。⑴分别为Golayl和Golay2的背散射回波信号。
[0018] 本发明中，对解码压缩背散射信号采用基于AR模型的倒谱估计方法，提取背散射 信号(t)倒谱峰值，利用该峰值估计骨组织微观结构信息MSS。
[0019] 本发明提供的骨骼超声系统中编码增强的聚焦超声骨组织微观结构检测方法，应 用于生物医学检测、组织定征领域，与现有检测技术及方法相比，具有如下优势：
[0020] (1)激励信号时，激励信号持续时间长于传感器脉冲响应时间，在不增加峰值声功 率的情况下增加超声信号携带能量，提高主瓣峰值，抑制旁瓣，提高骨组织微观结构检测的 分辨率。
[0021] (2)减少声衰减和随机噪声对测量结果的影响，提高MSS估计结果鲁棒性。
【专利附图】
【附图说明】
[0022] 图1为本发明中基于骨微结构显微CT三维重建的等效散射子示意图；
[0023] 图2为本发明仿真试验中沿波传播方向上声速分布示意图；
[0024] 图3为本发明仿真试验中沿波传播方向上组织材料特性分布示意图；
[0025] 图4为本发明中正弦信号调制Golay相位编码序列示意图；
[0026] 图5为本发明中Golay编码增强的骨微结构检测方法框图；
[0027] 图6为本发明中编码增强的单次编码背散射回波信号；
[0028] 图7为本发明中经过解码压缩的背散射回波信号；
[0029] 图8为本发明中未经编码增强的MSS估计结果；
[0030] 图9为本发明中经编码增强的MSS估计结果；
[0031] 图10为本发明中存在噪声（SNR = 45dB)情况下，未经编码增强的MSS估计结果；
[0032] 图11为本发明中存在噪声（SNR = 45dB)情况下，经编码增强的MSS估计结果。
【具体实施方式】
[0033] 下面结合附图对本发明的技术方案作进一步的说明，但并不局限于此，凡是对本 发明技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围，均应涵盖 在本发明的保护范围中。
[0034] 编码增强的聚焦超声骨组织微观结构检测方法采用超声背散射检测方法，采用单 一超声换能器，测量环境为水浸式测量，检测对象为能够表征骨组织微观结构特征的等效 模型MSS。等效模型如图1所示，首先，利用显微CT松质骨微观结构进行扫描和三维重建， 并得到特定位置特定路径上的骨小梁分布信息，根据该信息利用Faran圆柱模型去等效骨 组织微观结构，测量对象为该Faran圆柱模型的等效散射子间距MSS。
[0035] 图2-3给出了本发明Faran圆柱模型中沿波传播方向上组织材料特性及声速分 布示意图。散射子沿超声波传播方向等间距排列，散射子密度和散射子间介质密度分别为 938kg/m 3和1300kg/m3,声速分布分别为1500m/s和3300m/s，超声激励频率为0. 7?2MHz 可调，激励的超声能量水平低于lOOmW/cm2。超声诱导声压通过下式进行计算：
[0036] p = Ζ ω ξ，
[0037] 式中z = cQ*P为尸阻抗，cQ为超尸传播速度，Ρ为介质患'度；ω =2Jif为角频 率，ξ为质点位移。
[0038] 在用互补Golay序列对激励换能器时，根据被测骨组织厚度，利用FPGA编程控制 选择构造4位、8位或10位互补序列对[Golayl/Golay2]来激励超声换能器。Golay互补序 列对的构造方式如下：已知η位互补序列对为A[n]、B[n]，以A[n]、B[n]为基构造2n位互 补序列对{A[n] |B[n]}、{A[n] |?B[n]}，其中?B[n]为B[n]的补余序列对。本发明中以 2 位序列对 A[2] = [+1，+1]、B[2] = [+1，-1]为基，构造的 4 位 A[4]B[4]、8 位 A[8]B[8]、 10位A [10] B [10]序列对为：
[0039] [ (+1，+1，+1，_1) / (+1，+1，_1，+1)];
[0040] [ (+1，+1，+1，-1，+1，+1，-1，+1) / (+1，+1，+1，-1，-1，-1，+1，-1)];
[0041 ] [(+1, +1, -1, +1, -1, +1, -1, -1, +1, +1) / (+1, +1, -1, +1, +1, +1, +1, +1, -1, -1)]。
[0042] 超声检测过程中，通过控制面板选择Golay编码类型，控制FPGA调用相应程序激 励换能器。
[0043] 在利用Golay序列激励超声换能器时，通过对相位编码序列进行调制。调制过程 采用正弦信号作为基础序列，Golay相位编码序列中+1采用相位为0正弦信号进行调制，-1 采用相位为η的正弦信号进行调制。基础序列频率为聚焦换能器中心频率f P。4位编码 调制结果如图4所示。
[0044] 利用聚焦换能器对同一生物组织聚焦点分别采用经过编码单元和调制单元的脉 冲信号激励换能器，发射波持续时间Δ tin。< Lb/vb，Lb为骨组织厚度，vb为骨组织声速。编 码回波信号F bsc;(t)、F' bsc;(t)采用数据采集卡采集。同时在解码压缩过程中，利用数据采 集卡分别采集经过调制的两次编码发射信号F inc;(t)、F' inc;(t)作为匹配滤波函数，该匹配 滤波函数与对应编码回波信号做相关后求和，得到用于骨组织MSS估计的背散射信号：
[0045] Fcep(t) = cov[Finc(t)Fbsc(t)]+cov[F/inc ⑴ bsc (t)]。
[0046] 该MSS测量过程如图5所示。其中，数据卡采集到的单次编码回波信号Fbsc;(t)、 F' bsc;(t)如图6所示；经过解码压缩的背散射回波信号F_(t)如图7所示。
[0047]利用该背散射信号，采用AR倒谱估计方法，背散射信号Frap (t)倒谱峰值即为骨组 织微观结构信息MSS。未经编码增强的MSS估计结果和编码增强的MSS估计结果如图8-9 所示，可以看出经编码增强后，AR倒谱估计能够得到准确的MSS值，正确倒谱峰值得到增 强，分辨率得到提高，同时对AR倒谱干扰峰值起到抑制作用；在SNR = 45dB的高斯噪声影 响下，未经编码增强的MSS估计结果和编码增强的MSS估计结果如图10-11所示，可以看出 经编码增强MSS估计方法倒谱峰值明显，对高斯噪声具有比较强的抗干扰能力。
【权利要求】
1. 一种骨骼超声系统中编码增强的聚焦超声骨组织微观结构检测方法，其特征在于所 述骨组织微观结构检测方法步骤如下： 一、 采用超声背散射检测方法，依据Faran圆柱模型，将骨组织微观结构的检测等效为 平均散射子间距的测量； 二、 采用高频单阵元聚焦换能器，利用互补Golay序列对激励超声换能器发射宽时信 号，以发射信号作为匹配滤波函数对回波信号进行解码压缩，利用解码压缩背散射信号估 计平均散射子间距。
2. 根据权利要求1所述的骨骼超声系统中编码增强的聚焦超声骨组织微观结构检测 方法，其特征在于所述步骤一中，超声背散射检测的测量方式为单一探头的背散射测量方 式，测量环境为水浸式测量。
3. 根据权利要求1所述的骨骼超声系统中编码增强的聚焦超声骨组织微观结构检测 方法，其特征在于所述步骤一中，根据骨微结构显微CT三维重建结果，将骨微结构等效为 间距排列的Faran圆柱散射子，利用该等效散射子间距来表征骨微结构。
4. 根据权利要求1所述的骨骼超声系统中编码增强的聚焦超声骨组织微观结构检测 方法，其特征在于所述步骤一中，采用的超声激励频率为0. 7?2MHz可调，激励的超声能量 水平低于100mW/cm2。
5. 根据权利要求1所述的骨骼超声系统中编码增强的聚焦超声骨组织微观结构检测 方法，其特征在于所述步骤二中，在用互补Golay序列对激励超声换能器时，根据被测骨组 织厚度，选择构造4位[(+1，+1，+1， -1)/(+1，+1，-1，+1)]、8位[(+1，+1，+1，-1，+1，+1， -1， +1) / (+1，+1，+1，-1，-1，-1，+1，-1)]或 10 位[(+1，+1，-1，+1，-1，+1，-1，-1，+1，+1)/ (+1， +1，-1，+1，+1，+1，+1，+1，-1，-1)]互补序列对[Golayl/Golay2]，对同一生物组织聚焦点 分别采用Golayl、Golay2两次激励超声换能器，发射波持续时间Λ tin。< Lb/vb，Lb骨组织 厚度，vb为骨组织声速。
6. 根据权利要求1或5所述的骨骼超声系统中编码增强的聚焦超声骨组织微观结构检 测方法，其特征在于所述步骤二中，在激励超声换能器时，采用正弦信号作为基础信号调制 Golay互补序列。
7. 根据权利要求1所述的骨骼超声系统中编码增强的聚焦超声骨组织微观结构检测 方法，其特征在于所述步骤二中，在对背散射回波信号进行解码压缩过程中，采集Golayl 和Golay2输入信号Finc;(t)F' inc;(t)作为匹配滤波函数，则获得背散射信号为： Fcep (t) = COV [Finc (t) Fbsc (t) ] +cov [F;inc (t) F;bsc (t)], 其中，Fbs。⑴F' bs。⑴分别为Golayl和Golay2的背散射回波信号。
8. 根据权利要求1所述的骨骼超声系统中编码增强的聚焦超声骨组织微观结构检测 方法，其特征在于所述步骤二中，对解码压缩背散射信号采用基于AR模型的倒谱估计方 法，提取背散射信号(t)倒谱峰值，利用该峰值估计骨组织微观结构信息MSS。
【文档编号】A61B8/00GK104146729SQ201410366115
【公开日】2014年11月19日 申请日期:2014年7月29日 优先权日:2014年7月29日
【发明者】沈毅, 潘文磊, 金晶, 王艳, 章欣 申请人:哈尔滨工业大学