调用np.divide计算NDVI出现uint16溢出、值超理论范围问题

问题原因 #

异常由uint16无符号整数的算术溢出导致：
uint16的取值范围为0~65535，不支持负数。当img8像素值小于img4对应像素值时，img8 - img4的负结果会触发下溢出，被自动转换为极大正整数，最终分子大于分母，得到超出[-1,1]范围的异常NDVI值，和你观测到的运行警告、异常结果完全对应。

解决方法 #

核心逻辑为计算NDVI前先将数组转换为支持负数的数值类型，推荐以下两种实现：

• 方案1：直接转换为浮点型计算，实现最简洁

# 所有运算自动提升为float32类型，避免溢出
ndvi = np.divide(img8.astype(np.float32) - img4, img8.astype(np.float32) + img4)

• 方案2：先转有符号整型做加减，再转浮点计算，适合对内存占用敏感的场景

# int16支持负数，内存占用和原uint16一致
img8_int = img8.astype(np.int16)
img4_int = img4.astype(np.int16)
ndvi = (img8_int - img4_int) / (img8_int + img4_int).astype(np.float32)

效果验证 #

以你给出的异常像素值为例，转换类型后计算结果为(1462 - 1652)/(1462 + 1652) ≈ -0.061，符合NDVI的理论取值范围。

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