本章节主要参考《python地理空间分析指南》第六章

---

一、图像波段变换

波段变换最常用的地方就是进行图像显示，例如使用假彩色图像以凸显植被信息。图像波段变换即将图像波段的组合顺序重新排列并显示的方法，主要使用GDAL库。下面以一个例子进行实验，首先给出数据的下载地址：

http://git.io/vqs41

打开原图可以看到，植被明显呈现红色（植被在绿光波段反射率高，而图像将绿光波段放到了红色通道），说明该图为假彩色影像：

在本示例中，将使用gdal_array模块加载图片到numpy数组中，然后将其另存为新的geotiff文件。需要注意的是，numpy引用数组的方式是（y，x）（行，列）型的，常用的电子表格及其他软件则是（x，y）（行，列）型。

下面给出波段变换的代码：

from osgeo import gdal_array

# 源图片的名称
src = "FalseColor.tif"
# 将源图片载入到数组中
arr = gdal_array.LoadFile(src)
# 交换波段1和波段2的位置，使用“高级分片”功能直接对波段进行重新排列
output = gdal_array.SaveArray(arr[[1, 0, 2], :], "swap.tif", format="GTiff",prototype=src)
# 取消输出，避免在某些平台上损坏文件
output = None

波段变换后的影像为：

二、创建直方图

 

三、图像裁剪

四、图像分类

五、图像特征提取

六、变化监测