我用我自己的数据计算了这些指数,确实发现了些问题。
数据为大气校正后的sentinel 2地表反射率影像,假彩色显示为

如果由SAVI、NDBI和MNDWI计算IBI,IBI的数值变得很不正常,图像线性拉伸呈现灰白色调,灰度直方图出现异常值


后来用Landsat 8数据也出现了这种情况,看来并不是偶然。翻阅徐涵秋老师的论文(Xu, Hanqiu. 2008. A new index for delineating built-up land features in satellite imagery. International Journal of Remote Sensing, 29, 4269–4276.),尝试用整合后的IBI计算公式:

重新计算的IBI就正常多了,直方图也没问题


稍微调整一下阈值分割显示建筑物区域

虽然有很多裸露土壤的影响,但是建筑物显示是正常的。
为什么公式整合前后出现这样的差异我也不大清楚,可能是整合前公式的分母存在太多偏0值,可能是MNDWI的值太大,不是搞这个研究领域的没去深究。楼主用后面那个公式试试看结果怎样。