爽 宋
北京工业大学 环境与生命学部(北京 100124),
Faculty of Environment and Life, Beijing University of Technology, Beijing 100124, P.R.China
智能化生理测量与临床转化北京市国际科研合作基地(北京 100124),
Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing 100124, P.R.China
英华 张
北京工业大学 环境与生命学部(北京 100124),
Faculty of Environment and Life, Beijing University of Technology, Beijing 100124, P.R.China
著黄 周
北京工业大学 环境与生命学部(北京 100124),
Faculty of Environment and Life, Beijing University of Technology, Beijing 100124, P.R.China
智能化生理测量与临床转化北京市国际科研合作基地(北京 100124),
Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing 100124, P.R.China
水才 吴
北京工业大学 环境与生命学部(北京 100124),
Faculty of Environment and Life, Beijing University of Technology, Beijing 100124, P.R.China
智能化生理测量与临床转化北京市国际科研合作基地(北京 100124),
Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing 100124, P.R.China
北京工业大学 环境与生命学部(北京 100124),
Faculty of Environment and Life, Beijing University of Technology, Beijing 100124, P.R.China
智能化生理测量与临床转化北京市国际科研合作基地(北京 100124),
Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing 100124, P.R.China
首都医科大学附属北京佑安医院(北京 100069),
Beijing You'an Hospital, Capital Medical University, Beijing 100069, P.R.China
Corresponding author.
1.2.2. 超声背散射零差 K 成像
超声参数成像可直观显示组织定征结果
[
12
]
。为了将超声背散射零差 K 成像用于猪肝组织微波消融消融监测,我们通过 Matlab 2016a 软件实现超声背散射信号零差 K 分布参数值估算及参数成像。
超声背散射零差 K 成像流程如
所示:(1)对采集到的原始超声背散射信号进行希尔伯特变换,取模后得到所监测组织的包络信号。包络信号经对数压缩和数字扫描变换后可重建出 B 超成像。(2)将采集到的原始超声背散射信号进行噪声辅助互相关算法(noise-assisted correlation algorithm,NCA)
[
19
,
22
]
预处理后,再进行包络检测。然后利用滑动窗口法
[
12
,
19
]
估算超声背散射零差 K 模型散射子数目参数(
μ
和
α
)。最后,对数字扫描变换后得到的估算值矩阵使用颜色映射,得到
μ
或
α
参数图像,即为超声零差 K 成像。
由于在无回声区域存在的噪声信号会导致零差 K 成像伪影,为避免噪声的影响,我们引入 NCA
[
19
,
22
]
。NCA 是将随机人工白噪声分 2 次分别加入每帧超声背散射信号,然后计算这 2 帧包含额外噪声的超声背散射信号的相关系数,最后将相关系数小于阈值所对应的背散射信号的振幅设为零,即用零值代替无回声区域的噪声信号。本文通过 Matlab 软件中的 awgn 函数向原始超声射频信号中添加白噪声,awgn 函数的输入参数
snr
设置为 10。将阈值设置为 80%,将幅值为零的超声背散射信号估算出的零差 K 模型参数也设置为零,就能得到经 NCA 处理后的超声零差 K 成像。
滑动窗口是用一个方形窗口在包络信号中滑动,通过窗口内的超声背散射数据估计
μ
或
α
值,这些估值被赋给窗口中心的像素。窗口以一定的重叠率滑过整个包络信号后,得到
μ
或
α
的估算值矩阵。由于使用滑动窗口算法得到的估算值矩阵的尺寸会小于原始包络信号矩阵的尺寸,因此还需要对
μ
和
α
的估算值矩阵进行数字扫描变换。滑动窗口的大小会影响超声零差 K 成像的分辨率和参数估值的稳定性
[
19
]
,较大的窗口可获得更稳定的参数估算,但会降低成像分辨率;较小的窗口可获得更高的成像分辨率,但会降低参数估算的稳定性。综合考虑成像的分辨率和参数估算的稳定性,我们参阅文献
[
12
,
19
]
并反复实验,最终选择窗口边长等于入射超声脉冲长度的 3 倍,窗口的重叠率为 50%。