一种采用双尺度的生物组织位移估计方法技术

本发明专利技术为了解决选取大尺度和小尺度的矛盾，提出一种采用双尺度的组织位移估计方法，在组织压缩量比较大的情况下，选择大尺度和小尺度进行组织位移估计相结合的方法，对这两种组织位移估计的优缺点实现了互补，提高了组织位移估计的精度；在利用组织位移求组织对应的应变时，采用最小二乘拟合法，大大提高了组织应变估计的精度。

【技术实现步骤摘要】

本专利技术涉及超声弹性成像
，具体地涉及一种生物组织位移估计方法。
技术介绍
在超声弹性成像中，最关键的问题在于对组织的位移分布进行估计。选择压缩前 信号中的小段数据用来跟踪对应的小段组织的位移，称为跟踪波段，其长度称为尺度。选择 合适的尺度在超声弹性成像的位移估计中非常重要。当对组织施加的压缩量比较小的时 (如小于1 ，所选的尺度越大，包含的信息越多，位移估计的精度越高，对随机噪声的干 扰也越不敏感；当对组织施加的压缩量比较大的时(如1-5% )，由于组织被压缩，组织压缩 前、后形变较大，尺度越大位移估计的精度越低，同时，尺度越大又越不容易受到随机噪声 的干扰；而当尺度较小的时候，组织压缩前、后的形变较小，位移估计的精度越高，然而此时 的位移估计又易受到随机噪声的干扰。弹性成像技术最早由Jonathan Ophir于1991年提出，后在理论和技术上都获得 了飞速发展。由于其可获得独特的组织特征信息，使得其在某些应用场合比如乳腺，对于恶 性肿瘤和良性肿瘤的诊断具有很高的特异性和敏感性，从而在临床上获得了广泛应用。生物组织的弹性反映了其硬度，与病灶的生物学特征密切相关，但是其和组织的 散射并无直接关系，因此，B超图像不能很好的反应组织的这一生物学特征。临床病理学的 研究证明，脂肪、乳腺、纤维化组织、非浸润性癌变和浸润性癌变组织的弹性系数存在着明 显差异，这意味着弹性系数对于病变组织的定性分析具有重要价值。由于不同组织和同一种组织不同病变情况下弹性系数的差异，在外力作用下，其 发生形变的程度就会不同。超声弹性成像的基本原理为将超声探头嵌于一块挤压板中，沿 着探头的纵向压缩组织，分别采集组织压缩前、后的射频信号；组织被压缩时，组织内将会 产生一个沿压缩方向的应变，如果组织内部弹性模量分布不均勻，组织内的应变分布也会 有所差异；弹性模量较大的区域，引起的应变比较小；反之，弹性模量较小的区域，相应的 应变比较大。通过一些方法估计出组织内部不同位置的位移，从而计算出组织内部的应变 分布情况，用来间接描述组织内部的弹性模量分布，从而描述组织的生理、病理状态。对于二维超声弹性成像，一般采用线阵的B型超声探头，采集组织压缩前、后的每 条扫描线的射频信号，分别按照超声弹性成像原理计算出每一条扫描线对应组织的一维应 变分布。最后把所有扫描线对应的一维应变分布按扫描线顺序组成一个二维应变分布，以 灰度图或伪彩图显示，用来间接描述组织内部的弹性模量分布。现有的超声弹性成像步骤为1.使用全数字B型超声仪器(线阵探头)得到待测生物组织(或弹性仿体)压缩 前的一副数字化的二维射频信号；2.手持该B型超声仪器探头(或利用步进电动机)，沿着探头的纵向对该组织施 加一个微小的挤压(控制在之内)，得到组织压缩后的一幅数字化二维射频信号；3.确定需要计算的范围并设为R，从步骤1和2得到的组织压缩前后的二维射频4信号中选择一条扫描线L (L-R≥1且L+R≤Lmax，Lmax为扫描线的最大值)，将压缩前序号在 L-R:L+R的扫描线数据设为PreS(L)，同理压缩后的扫描线数据设为PostS(L)；4.从压缩前信号I^reS(L)中取一小段长度为T的数据(I1，数据个数为U，U = round (Τ · U1),其中T的单位为mm，U1代表Imm的组织对应的数据个数，round代表取整操作，该 数据Cl1的序号从Ii1到ni+U_l，ni可在1≤n1≤U的范围内选择；在确定的搜索范围τ 到 τ2内，求该小段数据与压缩后数据d2的绝对值差和，其中(12在^到τ 2内取，其规模和大 小与Cl1相同，计算公式如下ε (i) =Σ I Cl1-Cl2 (i) I τ 1≤i ≤ τ2其中τ工为0，τ2为对组织施加的压缩量，以采样数据的个数表示；5.确定步骤4得到的绝对值差和函数ε (i)的最小值对应的位置t1;、就是数据 Cl1在组织压缩后的位移；6.依次从压缩前数据I^reS(L)中取小段长度为T即数据个数为U的数据d2，d3， A dn，每段数据的序号依次错开V个采样数据1≤V≤U,直到在错开V个采样数据将超出 PreS(L)的范围，按步骤4、5相同的方法依次得到各段数据对应的位移t2，t3，Atn，其中η 为小段数据的总数；则位移序列t1;t2，t3，Atn为压缩前数据ft~eS(L)对应的组织的位移估 计；7.利用步骤3-6相同的方法，依次得到第L+l，L+2，ΛLmax-L扫描线数据对应的生 物组织的位移估计；8.对ft~eS(L)对应的组织的位移估计序列t1; t2，t3，Atn求差分，得到生物组织 PreS(L)对应组织的应变分布，计算公式如下s1= (t2-t1)/v,S2 = (t3-t2)/v ,∧s(n-1) = (tn-tn-1)/v其中，Sl，s2，S3, A Slri分别为屯，屯，d3，A Cllri对应的组织应变；9.利用步骤8相同的方法，依次得到第R+l，R+2，Λ Lmax-R条扫描线数据对应生物 组织的应变分布；10.将步骤9得到的Lmax_2 · R条扫描线数据对应的应变分布插值成Lmax条，按照 扫描线的顺序组合成一个二维数据，并以灰度图或伪彩图的形式表示出来，得到组织的二 维应变分布图。在超声弹性成像中，选择压缩前的信号的中的小段数据用来跟踪对应的小段组织 的位移，称为跟踪波段。其长度称为跟踪波段程度，或者尺度。选择合适的尺度在超声弹性 成像的位移估计中非常重要。当对组织施加的压缩量比较小(如小于1%)的时候，尺度越 大，包含的信息越多，位移估计的精度越高，对随机噪声的干扰也越不敏感。当对组织施加 的压缩量比较大(如1-5%)的时候，尺度越大，由于组织被压缩，组织压缩前后信号重合部 分越小，利用绝对值差和最小值位置进行位移估计的精度越低，同时，尺度越大又越不容易 受到随机噪声的干扰；而当尺度较小的时候，组织压缩前、后信号重合的部分较大，位移估 计的精度越高，然而此时的位移估计有对随机噪声的干扰比较敏感。
技术实现思路
本专利技术为了解决对组织施加的压缩量比较大的时候大尺度和小尺度的矛盾，提出 一种采用两种尺度的组织位移估计方法。在组织压缩量比较大时，将选择大尺度和小尺度 相结合的组织位移估计方法，可对这两种组织位移估计的优缺点实现扬长避短，提高组织 位移的估计精度。本专利技术实现专利技术目的采用的技术方案是1、，包括以下步骤a.使用数字化超声仪器测量并得到待测生物组织在压缩前的一副数字化的二维 射频信号；b.沿着数字化超声仪器的线阵探头的纵向对待测生物组织施加一个微小的挤压， 压缩量控制在之内，得到压缩后的一幅数字化二维射频信号；c.确定需要计算的范围并设为R，从步骤a和b得到的待测生物组织压缩前后的 二维射频信号中各选择一条扫描线L(L-R彡1且L+R彡Lmax，Lfflax为扫描线的最大值)，将压 缩前的序号在L-R: L+R的扫描线数据设为I^reS (L)，压缩后的序号在L-R: L+R的扫描线数据 设为 PostS (L)；d.从压缩前的扫描线数据I^reS(L)中取一小段长度为Ta的数据(I1，数据个数为 Ua, Ua = round (Ta .U1)，其中Ta的单位为HinuU1代表Imm的生本文档来自技高网...

【技术保护点】
一种采用双尺度的生物组织位移估计方法，包括以下步骤：ａ．使用数字化超声仪器测量并得到待测生物组织在压缩前的一副数字化的二维射频信号；ｂ．沿着数字化超声仪器的线阵探头的纵向对待测生物组织施加一个微小的挤压，压缩量控制在１％之内，得到压缩后的一幅数字化二维射频信号；ｃ．确定需要计算的范围并设为Ｒ，从步骤ａ和ｂ得到的待测生物组织压缩前后的二维射频信号中各选择一条扫描线Ｌ（Ｌ－Ｒ≥１且Ｌ＋Ｒ≤Ｌ↓［ｍａｘ］，Ｌ↓［ｍａｘ］为扫描线的最大值），将压缩前的序号在Ｌ－Ｒ：Ｌ＋Ｒ的扫描线数据设为ＰｒｅＳ（Ｌ），压缩后的序号在Ｌ－Ｒ：Ｌ＋Ｒ的扫描线数据设为ＰｏｓｔＳ（Ｌ）；ｄ．从压缩前的扫描线数据ＰｒｅＳ（Ｌ）中取一小段长度为Ｔ↓［ａ］的数据ｄ↓［１］，数据个数为Ｕ↓［ａ］，Ｕ↓［ａ］＝ｒｏｕｎｄ（Ｔ↓［ａ］.Ｕ↓［１］），其中Ｔ↓［ａ］的单位为ｍｍ，Ｕ↓［１］代表１ｍｍ的生物组织对应的数据个数，ｒｏｕｎｄ代表取整操作，该数据ｄ↓［１］的序号从ｎ↓［１］到ｎ↓［１］＋Ｕ↓［ａ］－１，ｎ↓［１］可在１≤ｎ↓［１］≤Ｕ↓［ａ］的范围内选择；在确定的搜索范围τ↓［１］到τ↓［２］内，求该小段数据ｄ↓［１］与压缩后数据ｄ↓［２］的绝对值差和，其中ｄ↓［２］在τ↓［１］到τ↓［２］内取，其规模和大小与ｄ↓［１］相同，计算公式如下：ε′（ｉ）＝∑｜ｄ↓［１］－ｄ↓［２］（ｉ）｜τ↓［１］≤ｉ≤τ↓［２］其中τ↓［１］为０，τ↓［２］为对组织施加的压缩量，以采样数据的个数表示；ｅ．确定步骤ｄ得到的绝对值差和函数ε′（ｉ）的最小值对应的位置ｔ↓［１］′，ｔ↓［１］′就是利用尺度Ｔ↓［ａ］进行位移估计时数据ｄ↓［１］在待测组织压缩后的位移；ｆ．从压缩前的扫描线数据ＰｒｅＳ（Ｌ）中的数据段ｄ↓［１］的中央位置再取尺度为Ｔ↓［ｂ］、序号从ｎ↓［１］＋ｒｏｕｎｄ（Ｕ↓［ａ］－Ｕ↓［ｂ］／２）到ｎ↓［１］＋ｒｏｕｎｄ（Ｕ↓［ａ］－Ｕ↓［ｂ］／２）＋Ｕ↓［ｂ］－１的一小段数据ｅ↓［１］，该数据的个数为Ｕ↓［ｂ］，Ｕ↓［ｂ］＝ｒｏｕｎｄ（Ｔ↓［ｂ］×Ｕ↓［１］），其中Ｔ↓［ｂ］的单位为ｍｍ，并且满足Ｔ↓［ｂ］≤Ｔ↓［ａ］，在τ↓［３］到τ↓［４］确定的搜索范围内求该小段数据与扫描线数据ｅ↓［２］的绝对值差和，其中ｅ↓［２］在τ↓［３］到τ↓［４］内取，其规模和大小与ｄ↓［２］相同，计算公式如下：ε″（ｉ）＝∑｜ｅ↓［１］－ｅ↓［２］（ｉ）｜τ↓［３］≤ｉ≤τ↓...

【技术特征摘要】

【专利技术属性】
技术研发人员：齐育彬，
申请(专利权)人：深圳市蓝韵实业有限公司，
类型：发明
国别省市：94[中国|深圳]