终显示的影像编写脚本计算不同类型的指数，如归一化水体指数（Normalized Difference Water Index, NDWI）、归一化建筑指数（Normalized Difference Built-up Index, NDBI）等，并在地图上显示结果。

编写脚本来计算和显示不同类型的指数，比如归一化水体指数（NDWI）和归一化建筑指数（NDBI），通常涉及到遥感图像处理和地理信息系统（GIS）技术。以下是一个简化的步骤：

1. **导入所需库**：首先，你需要导入像GDAL（用于读取和操作地理空间数据）、NumPy（数值计算库）和matplotlib（可视化工具）这样的库。

import gdal
import numpy as np
import matplotlib.pyplot as plt

2. **读取遥感图像**：使用GDAL库加载需要分析的卫星或航空影像数据。

image = gdal.Open('path_to_image.tif')
band1 = image.GetRasterBand(1)
band2 = image.GetRasterBand(2)

3. **计算NDWI**：NDWI通常是基于近红外波段（例如B8）和绿波段（例如B4）的差值除以它们的和。公式大致如下：

ndwi = (band1_data - band2_data) / (band1_data + band2_data)

4. **计算NDBI**：NDBI则可能涉及短波红外（SWIR，例如B5）和近红外的比值。类似地：

ndbi = (band5_data - band1_data) / (band5_data + band1_data)

5. **归一化**：将计算出的数据进行归一化，以便于后续的地图显示。

6. **创建RGB图或伪彩色图**：结合可见光波段（例如红、绿、蓝）和处理后的指数数据，创建一个可以显示在地图上的图像。

rgb_array = ... # 使用可见光波段构建RGB数组
colorized_ndwi = ... # 将NDWI映射到颜色范围
colorized_ndbi = ...

fig, ax = plt.subplots()
ax.imshow(rgb_array)
ax.imshow(colorized_ndwi, cmap='RdYlGn', alpha=0.5)
ax.imshow(colorized_ndbi, cmap='coolwarm', alpha=0.5)

7. **保存和显示**：最后，保存图像并显示结果。

plt.show()