Matlab中含零元素共生矩阵的熵计算方法问询

你提到的这个点其实是GLCM熵计算里很常见的细节，先明确一下GLCM熵的标准定义公式：

熵 = $-\sum_{i,j} c_{ij} \log_2(c_{ij})$（这里通常用2为底的对数，也有用自然对数的，取决于具体分析标准）

这里的核心问题就是当$c_{ij}=0$时，$0 \times \log(0)$这个项的处理——从数学极限的角度来说，$\lim_{x \to 0^+} x\log(x) = 0$，所以我们可以直接把所有$c_{ij}=0$对应的项视为0，根本不会对最终的熵值产生贡献，完全不用担心熵会无定义。

在Matlab里实现的时候，为了避免直接计算log(0)出现警告或者NaN值，你可以用这个实用的小技巧：

• 先确保你的GLCM是归一化后的概率矩阵（熵计算通常基于概率分布，而非原始计数）
• 通过逻辑索引把0元素对应的对数结果置0，再进行计算

示例代码如下：

% 假设glcm是你已生成的共生矩阵
glcm_normalized = glcm / sum(glcm(:)); % 归一化为概率矩阵
log_glcm = log2(glcm_normalized); % 计算对数
log_glcm(glcm_normalized == 0) = 0; % 将0元素对应的对数结果置0
entropy = -sum(glcm_normalized(:) .* log_glcm(:)); % 计算熵值

这种做法既完全符合熵的数学定义，又巧妙避开了log(0)的报错问题。你之前的误解主要是没注意到公式里是0乘以log(0)，而非单独计算log(0)——前者在数学上是有明确取值的，后者才会无意义~

内容的提问来源于stack exchange，提问作者Postermaestro