Gao Ruiyang, Tsui Po-Hsiang, Li Sinan, Bin Guangyu, Tai Dar-In, Wu Shuicai, Zhou Zhuhuang
Department of Biomedical Engineering, College of Chemistry and Life Science, Beijing University of Technology, Beijing, China.
Department of Medical Imaging and Radiological Sciences, College of Medicine, Chang Gung University, Taoyuan, Taiwan; Division of Pediatric Gastroenterology, Department of Pediatrics, Chang Gung Memorial Hospital at Linkou, Taoyuan, Taiwan; Liver Research Center, Chang Gung Memorial Hospital, Linkou, Taoyuan, Taiwan; Research Center for Radiation Medicine, Chang Gung University, Taoyuan, Taiwan.
Comput Methods Programs Biomed. 2024 Nov;256:108374. doi: 10.1016/j.cmpb.2024.108374. Epub 2024 Aug 13.
Ultrasound information entropy imaging is an emerging quantitative ultrasound technique for characterizing local tissue scatterer concentrations and arrangements. However, the commonly used ultrasound Shannon entropy imaging based on histogram-derived discrete probability estimation suffers from the drawbacks of histogram settings dependence and unknown estimator performance. In this paper, we introduced the information-theoretic cumulative residual entropy (CRE) defined in a continuous distribution of cumulative distribution functions as a new entropy measure of ultrasound backscatter envelope uncertainty or complexity, and proposed ultrasound CRE imaging for tissue characterization.
We theoretically analyzed the CRE for Rayleigh and Nakagami distributions and proposed a normalized CRE for characterizing scatterer distribution patterns. We proposed a method based on an empirical cumulative distribution function estimator and a trapezoidal numerical integration for estimating the normalized CRE from ultrasound backscatter envelope signals. We presented an ultrasound normalized CRE imaging scheme based on the normalized CRE estimator and the parallel computation technique. We also conducted theoretical analysis of the differential entropy which is an extension of the Shannon entropy to a continuous distribution, and introduced a method for ultrasound differential entropy estimation and imaging. Monte-Carlo simulation experiments were performed to evaluate the estimation accuracy of the normalized CRE and differential entropy estimators. Phantom simulation and clinical experiments were conducted to evaluate the performance of the proposed normalized CRE imaging in characterizing scatterer concentrations and hepatic steatosis (n = 204), respectively.
The theoretical normalized CRE for the Rayleigh distribution was π/4, corresponding to the case where there were ≥10 randomly distributed scatterers within the resolution cell of an ultrasound transducer. The theoretical normalized CRE for the Nakagami distribution decreased as the Nakagami parameter m increased, corresponding to that the ultrasound backscattered statistics varied from pre-Rayleigh to Rayleigh and to post-Rayleigh distributions. Monte-Carlo simulation experiments showed that the proposed normalized CRE and differential entropy estimators can produce a satisfying estimation accuracy even when the size of the test samples is small. Phantom simulation experiments showed that the proposed normalized CRE and differential entropy imaging can characterize scatterer concentrations. Clinical experiments showed that the proposed ultrasound normalized CRE imaging is capable to quantitatively characterize hepatic steatosis, outperforming ultrasound differential entropy imaging and being comparable to ultrasound Shannon entropy and Nakagami imaging.
This study sheds light on the theory and methodology of ultrasound normalized CRE. The proposed ultrasound normalized CRE can serve as a new, flexible quantitative ultrasound envelope statistics parameter. The proposed ultrasound normalized CRE imaging may find applications in quantified characterization of biological tissues. Our code will be made available publicly at https://github.com/zhouzhuhuang.
超声信息熵成像作为一种新兴的定量超声技术,用于表征局部组织散射体的浓度和排列。然而,基于直方图离散概率估计的常用超声香农熵成像存在直方图设置依赖性和估计器性能未知的缺点。在本文中,我们引入了在累积分布函数连续分布中定义的信息论累积剩余熵(CRE),作为超声背向散射包络不确定性或复杂性的新熵度量,并提出了用于组织表征的超声CRE成像。
我们从理论上分析了瑞利分布和 Nakagami 分布的 CRE,并提出了一种归一化 CRE 来表征散射体分布模式。我们提出了一种基于经验累积分布函数估计器和梯形数值积分的方法,用于从超声背向散射包络信号中估计归一化 CRE。我们提出了一种基于归一化 CRE 估计器和平行计算技术的超声归一化 CRE 成像方案。我们还对作为香农熵向连续分布扩展的微分熵进行了理论分析,并介绍了一种超声微分熵估计和成像方法。进行了蒙特卡罗模拟实验,以评估归一化 CRE 和微分熵估计器的估计精度。分别进行了仿体模拟和临床实验,以评估所提出的归一化 CRE 成像在表征散射体浓度和肝脂肪变性方面的性能(n = 204)。
瑞利分布的理论归一化 CRE 为π/4,对应于超声换能器分辨率单元内有≥10 个随机分布散射体的情况。Nakagami 分布的理论归一化 CRE 随着 Nakagami 参数 m 的增加而降低,这对应于超声背向散射统计从瑞利前分布变为瑞利分布再变为瑞利后分布。蒙特卡罗模拟实验表明,即使测试样本量较小,所提出的归一化 CRE 和微分熵估计器也能产生令人满意的估计精度。仿体模拟实验表明,所提出的归一化 CRE 和微分熵成像能够表征散射体浓度。临床实验表明,所提出的超声归一化 CRE 成像能够定量表征肝脂肪变性,优于超声微分熵成像,与超声香农熵成像和 Nakagami 成像相当。
本研究阐明了超声归一化 CRE 的理论和方法。所提出的超声归一化 CRE 可作为一种新的、灵活的定量超声包络统计参数。所提出的超声归一化 CRE 成像可能在生物组织的定量表征中找到应用。我们的代码将在 https://github.com/zhouzhuhuang 上公开提供。