一种三维超声断层成像系统原始信号质量判断方法技术方案

本发明专利技术属于超声断层成像领域，具体涉及一种三维超声断层成像系统原始信号质量判断方法，包括：分别获得渡越时间法质量矩阵和能量分布法质量矩阵；对渡越时间法质量矩阵和能量分布法质量矩阵按照规则处理，得到新的质量矩阵。本发明专利技术可以高效地对大数据量的原始信号质量进行判断，同时该算法支持并行计算，可以很好地进行加速，做到近实时判断。做到近实时判断。做到近实时判断。

【技术实现步骤摘要】
一种三维超声断层成像系统原始信号质量判断方法


[0001]本专利技术属于超声断层成像领域，具体涉及一种三维超声断层成像系统原始信号质量判断方法。

技术介绍

[0002]超声断层成像系统一次同时重建反射和透射图像，原始信号的质量对于图像质量非常重要。因为超声断层成像系统并行采集通道数多，采集链路硬件和控制机制复杂，一次采集数据量庞大(82GB)，现没有一种定量的高效快速原始信号质量判断方法。
[0003]基于此，申请人提出一种基于渡越时间和能量分布的三维超声断层成像系统原始信号质量判断方法。

技术实现思路

[0004]为了弥补现有技术的不足，本专利技术提供一种三维超声断层成像系统原始信号质量判断方法的技术方案。
[0005]一种三维超声断层成像系统原始信号质量判断方法，包括：
[0006]步骤1，分别获得渡越时间法质量矩阵和能量分布法质量矩阵；
[0007]步骤2，对渡越时间法质量矩阵和能量分布法质量矩阵按照规则处理，得到新的质量矩阵。
[0008]进一步地，所述步骤1中获得渡越时间法质量矩阵的操作包括：
[0009]S1，取一段原始AScan信号，求得该原始AScan信号的信号包络，并计算该原始AScan信号的理论到达时间；
[0010]S2，基于理论到达时间截取一段包络信号，计算该截取段包络信息的离散系数；
[0011]S3，判断离散系数是否大于0；
[0012]若大于0，则对截取的包络信号进行处理，得到实际渡越时间；
[0013]若小于等于0，则判断实际渡越时间为0；
[0014]S4，将实际渡越时间与理论渡越时间进行对比，标定该AScan信号质量好坏；
[0015]S5，遍历所有AScan信号，得到渡越时间法质量矩阵。
[0016]进一步地，所述S1中信号包络的求法包括：使用一个模板信号进行匹配滤波，模板信号为理论的发射信号或者实际采集到的信号，然后去掉直流分量，做希尔伯特变换得到包络。
[0017]进一步地，所述S1中理论到达时间的计算包括：发射阵元的空间坐标为(x
tx
,y
tx
,z
tx
)，接收阵元的空间坐标为(x
rx
，y
rx
，z
rx
)，理论到达时间其中c为声速。
[0018]进一步地，所述S2中截取一段包络信号的操作包括：基于t
direct
和窗长度winLen，取(t
direct
‑
winLen/2,t
direct
+winLen)段包络信号。
[0019]进一步地，所述S2中离散系数的计算公式为：离散系数＝信号方差/信号均值。
[0020]进一步地，所述S3中对截取的包络信号进行处理包括：以2*包络信号离散系数为阈值，保留大于该值部分，其余置为0；对剩余的值排序，得到最大的两个值val1和val2；
[0021]若val2>0，且val2比val1更早到达，则val2的到达时间即为信号实际到达时间，记为t
real
；
[0022]如果val1>0，则将val1的到达时间记为t
real
。
[0023]进一步地，所述S4中将实际渡越时间与理论渡越时间进行对比包括：如果|t
real
‑
t
direct
|>winLen/2，则标定该AScan信号质量为坏；否则标定AScan信号质量为好。
[0024]进一步地，所述步骤1中获得能量分布法质量矩阵的操作包括：
[0025]H1，取一个阵元发射所有阵元接收的AScan信号，并计算所有AScan信号的理论能量分布；
[0026]H2，求每一个AScan信号的绝对值的最大值，做归一化得到一个实际能量分布向量；
[0027]H3，计算理论能量分布与实际能量分布的相对误差，相对误差绝对值小于等于20％的发射阵元记为好，反之记为坏；
[0028]H4，遍历所有发射阵元，得到能量分布法质量矩阵。
[0029]进一步地，所述H1中AScan信号的理论能量分布计算方法包括：
[0030]a，分别计算发射阵元和接收阵元的连线与发射阵元法向、发射阵元X方向、接收阵元法向以及接收阵元X方向的夹角，其中，X方向的夹角是指连线在XY平面的投影线与X方向的夹角；
[0031]b，发射阵元通过其两个夹角和三维指向性，获得其指向性系数，接收阵元通过其两个夹角和三维指向性，获得其指向性系数，发射阵元的指向性系数与接收阵元的指向性系数相乘；
[0032]c，不变发射阵元，遍历所有接收阵元，计算相乘后的指向性系数，得到一位系数向量；
[0033]d，对该系数向量进行归一化处理。
[0034]进一步地，所述步骤2具体包括：
[0035]若渡越时间法质量矩阵中的AScan信号与能力分布法质量矩阵中的发射阵列同时为坏，则记信号为坏；
[0036]否则记信号为好；
[0037]得到新的质量矩阵。
[0038]与现有技术相比，本专利技术的有益效果是：
[0039]本专利技术可以高效地对大数据量的原始信号质量进行判断，同时该算法支持并行计算，可以很好地进行加速，做到近实时判断。
附图说明
[0040]图1为本专利技术流程图；
[0041]图2为本专利技术中获得渡越时间法质量矩阵的操作流程图；
[0042]图3为本专利技术中获得能量分布法质量矩阵的操作流程图；
[0043]图4为本专利技术中AScan信号的理论能量分布计算方法流程图；
[0044]图5为本专利技术中模拟计算圆形超声换能器三维指向性的XZ平面的能量分布曲线图；
[0045]图6为本专利技术中发射接收阵元连接线与法线及坐标系X方向夹角θ,β定义示意图；
[0046]图7为本专利技术中半球形孔径中发射阵元接收阵元连线与各自法向夹角θ角三维定义示意图；
[0047]图8为本专利技术最终获得的质量判断矩阵示意图，行表示发射阵元序号，列表示接收阵元序号，矩阵值为1或0，1表示质量判断结果为差，白色表示；0表示质量判断结果为好，黑色表示。
具体实施方式
[0048]在本专利技术的描述中，需要理解的是，术语“一端”、“另一端”、“外侧”、“上”、“内侧”、“水平”、“同轴”、“中央”、“端部”、“长度”、“外端”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本专利技术和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本专利技术的限制。
[0049]下面结合附图对本专利技术作进一步说明。
[0050]请参阅图1
‑
8，一种三维超声断层成像系统原始信号质量判断方法，包本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种三维超声断层成像系统原始信号质量判断方法，其特征在于，包括：步骤1，分别获得渡越时间法质量矩阵和能量分布法质量矩阵；步骤2，对渡越时间法质量矩阵和能量分布法质量矩阵按照规则处理，得到新的质量矩阵。2.根据权利要求1所述的一种三维超声断层成像系统原始信号质量判断方法，其特征在于，所述步骤1中获得渡越时间法质量矩阵的操作包括：S1，取一段原始AScan信号，求得该原始AScan信号的信号包络，并计算该原始AScan信号的理论到达时间；S2，基于理论到达时间截取一段包络信号，计算该截取段包络信息的离散系数；S3，判断离散系数是否大于0；若大于0，则对截取的包络信号进行处理，得到实际渡越时间；若小于等于0，则判断实际渡越时间为0；S4，将实际渡越时间与理论渡越时间进行对比，标定该AScan信号质量好坏；S5，遍历所有AScan信号，得到渡越时间法质量矩阵。3.根据权利要求2所述的一种三维超声断层成像系统原始信号质量判断方法，其特征在于，所述S1中信号包络的求法包括：使用一个模板信号进行匹配滤波，模板信号为理论的发射信号或者实际采集到的信号，然后去掉直流分量，做希尔伯特变换得到包络。4.根据权利要求2所述的一种三维超声断层成像系统原始信号质量判断方法，其特征在于，所述S1中理论到达时间的计算包括：发射阵元的空间坐标为(x
tx
,y
tx
,z
tx
)，接收阵元的空间坐标为(x
rx
，y
rx
，z
rx
)，理论到达时间其中c为声速。5.根据权利要求4所述的一种三维超声断层成像系统原始信号质量判断方法，其特征在于，所述S2中截取一段包络信号的操作包括：基于t
direct
和窗长度winLen，取(t
direct
‑
winLen/2,t
direct
+winLen)段包络信号。6.根据权利要求4所述的一种三维超声断层成像系统原始信号质量判断方法，其特征在于，所述S2中离散系数的计算公式为：离散系数＝信号方差/信号均值。7.根据权利要求2所述的一种三维超声断层成像系统原始信号质量判断方法，其特征在于，所述S3中对截取的包络信...

【专利技术属性】
技术研发人员：雷晓旭，韩春林，黄跃龙，
申请(专利权)人：浙江衡玖医疗器械有限责任公司，
类型：发明
国别省市：