主权项 |
1、一种变尺度的组织位移估计方法,包括以下步骤:1)从组织压缩前、后的二维射频信号中分别取出第一条扫描线的数据,设为s1(n)和s2(n),n表示该两条扫描线上的数据序号,1≤n≤nmax;2)从该扫描描线数据s1(n)中取一小段尺度为T(1)的数据d1,其数据个数为U,U=round(T(1)×U1),其中,T的单位为mm,U1代表1mm的组织对应的数据个数,round(·)代表四舍五入的取整操作,T(1)为尺度随扫描线序号k变化的函数T(k)在k=1时的取值,该数据d1的序号n1在1≤n1≤U的范围内选择;在τ1到τ2确定的搜索范围内计算该小段数据与扫描线数据s2(n)的互相关函数R(τ)如下<math> <mrow> <mi>R</mi> <mrow> <mo>(</mo> <mi>τ</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <munderover> <mi>Σ</mi> <mrow> <mi>m</mi> <mo>=</mo> <msub> <mi>n</mi> <mn>1</mn> </msub> </mrow> <mrow> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>+</mo> <mi>U</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>s</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <msub> <mi>s</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mi>τ</mi> <mo>)</mo> </mrow> </mrow> <msqrt> <munderover> <mi>Σ</mi> <mrow> <mi>m</mi> <mo>=</mo> <msub> <mi>n</mi> <mn>1</mn> </msub> </mrow> <mrow> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>+</mo> <mi>U</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msubsup> <mi>s</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mo>·</mo> <munderover> <mi>Σ</mi> <mrow> <mi>m</mi> <mo>=</mo> <msub> <mi>n</mi> <mn>1</mn> </msub> </mrow> <mrow> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>+</mo> <mi>U</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msubsup> <mi>s</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mi>τ</mi> <mo>)</mo> </mrow> </msqrt> </mfrac> <mrow> <mo>(</mo> <msub> <mi>τ</mi> <mn>1</mn> </msub> <mo>≤</mo> <mi>τ</mi> <mo>≤</mo> <msub> <mi>τ</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </math> 其中m为计算过程中表示数据序号的循环变量,τ1为0,τ2为对组织施加的压缩量,以采样数据的个数表示;3)确定该互相关函数R(τ)的最大值对应的位置t1,t1就是数据d1在组织压缩后的位移;4)依次从该扫描线数据s1(n)中取一小段尺度为T(1)即数据个数为U的数据d2、d3、…、dN,每段数据的序号依次错开V个采样数据,直到再错开V个采样数据将超出s1(n)的范围,按步骤2、3相同的方法依次得到各段数据对应的位移t2、t3、…、tN,其中N为小段数据的总数;则位移序列t1、t2、…、tN为第一条扫描线数据s1(n)对应的组织的位移估计;5)利用与步骤1-4相同的方法,依次得到第2、3、…、M条扫描线数据对应的组织的位移估计,其中M为表示探头的扫描线总数;其特征在于,所述的第k条扫描线数据采用的尺度为T(k),且如果M为奇数,T(k)满足T(1)=T(M)=Tmax,<math> <mrow> <mi>T</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>M</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <mo>=</mo> <mi>T</mi> <mi>min</mi> <mo>,</mo> </mrow> </math> 当1≤k≤(M+1)/2时,T(k)随着k的增大而减小,当(M+1)/2≤k≤M时,T(k)随着k的增大而增大;如果M为偶数,T(k)满足T(1)=T(M)=Tmax,T(M/2)=T(M/2+1)=Tmin,当1≤k≤M/2时,T(k)随着k的增大而减小,M/2+1≤k≤M时,T(k)随着k的增大而增大;其中,Tmax和Tmin分别为变尺度方法的最大尺度和最小尺度。 |