1.摘要
阿尔茨海默病（AD）的早期诊断高度依赖于脑病理图像的精确分割，但传统多阈值图像分割方法在噪声抑制和空间结构信息利用方面存在不足，难以应对 AD 图像中复杂纹理与高信息密度的问题。为此，本文提出了一种融合灰度强度、非局部均值和局部熵的三维 Rényi 熵模型，通过联合直方图同时表征灰度、空间与纹理特征，从而更全面地刻画图像不确定性。针对高维阈值优化难题，本文设计了一种量子混合电鳗觅食优化算法（QHEEFO），引入量子隧穿策略、量子控制因子和对数增强扰动机制，以提升全局搜索能力并避免早熟收敛。
2.基于 NLM 和局部熵的三维联合直方图
传统基于二维联合直方图的图像分割方法在医学图像中对复杂结构和模糊边界的表征能力有限。非局部均值滤波（NLM）利用图像的非局部自相似性，在抑制噪声的同时有效保留细节。将待处理的当前图像块的中心像素与搜索窗口内的所有像素进行比较，根据图像块之间的相似度来确定权重，并对图像进行加权平均运算：
$$ NL(i) = \frac{\sum_{j\in\omega(i,j)}\omega(i,j)\times I(j)}{\sum_{j\in\omega(i,j)}\omega(i,j)} $$
其中，权重系数：
$$ \omega(i,j)=\exp\left(-\frac{|\mu(i)-\mu(j)|}{\sigma^2}\right) $$
尽管 NLM 能有效建模非局部相似性，但其对复杂纹理和结构不确定性的刻画仍然不足。为此，引入局部熵作为第三特征通道，用于度量像素邻域灰度分布的不确定性，其定义为：
$$ e(i,j)=-\sum_{g=0}^{G-1}p_g(i,j)\ln\left(p_g(i,j)+\varepsilon\right) $$
将原始灰度 $t(i,j)$、NLM 滤波结果 $s(i,j)$ 与局部熵 $e(i,j)$ 组成三元特征向量 $(t,s,e)$，并构建三维联合直方图 $P(t,s,e)$。
3.3D Rényi 熵多阈值分割
传统基于二维联合直方图的 Rényi 熵多阈值分割方法在医学图像中难以有效表征复杂纹理和微观结构，尤其在组织边界模糊、纹理交织的情况下判别能力明显不足。为此，本文提出一种基于局部熵增强的三维（3D）Rényi 熵多阈值分割模型，引入原始灰度、非局部均值（NLM）和局部熵三种特征，构建更高判别力的 3D 联合特征空间。
对输入图像分别计算灰度图、NLM 滤波图和局部熵图，并为每个像素构建三元特征向量 $(t,s,e)$，进而通过量化生成三维联合直方图 $P(t,s,e)$。多阈值分割被建模为最大化 3D Rényi 熵的优化问题，其目标函数定义为：
$$ \arg\max_{(t_k,s_k,e_k)}R_\alpha=\sum_{k=1}^L\frac{1}{1-\alpha}\ln\left(\sum_{t=0}^{t_k}\sum_{s=0}^{s_k}\sum_{e=0}^{e_k}\left(\frac{P_{t,s,e}}{P_{Nk}}\right)^\alpha\right) $$
4.改进算法
量子隧穿策略（QTS）
EEFO 算法仅依据个体适应度进行位置更新，个体间交互完全由适应度差异驱动，在多峰或高维问题中灵活性不足，易导致搜索停滞和早熟收敛。为此，本文引入基于量子隧穿效应的改进策略，通过量子因子 $\alpha$ 与对数变换，使搜索个体能够跨越低适应度区域而不依赖适应度梯度，从而显著增强全局搜索能力并有效缓解早熟收敛问题：
$$ \begin{cases} x_i(t+1) &= x_j(t) + \alpha \times | \bar{x}(t) - x_j(t) | \times \log\left(\frac{1}{u}\right), & p_1 > 0.5, \ \text{if } fit(x_j(t)) < fit(x_i(t)) \ x_i(t+1) &= x_j(t) - \alpha \times | x_r(t) - x_i(t) | \times \log\left(\frac{1}{u}\right), & p_1 \le 0.5, \ \text{if } fit(x_j(t)) < fit(x_i(t)) \ x_i(t+1) &= x_i(t) + \alpha \times | \bar{x}(t) - x_j(t) | \times \log\left(\frac{1}{u}\right), & p_1 > 0.5, \ \text{if } fit(x_j(t)) > fit(x_i(t)) \ x_i(t+1) &= x_i(t) - \alpha \times | x_r(t) - x_i(t) | \times \log\left(\frac{1}{u}\right), & p_1 \le 0.5, \ \text{if } fit(x_j(t)) > fit(x_i(t)) \end{cases} $$
