python 遥感图像分类

图像分类的常用方法一般包括监督分类和非监督分类两大类，本文采用阈值切片法对遥感影像进行分类，使用数据为landsat8遥感影像，下载地址：

http://github.com/GeospatialPython/Learn/blob/master/thermal.zip?raw=true

该影像为单波段遥感影像

利用切片法对上述影像进行分类，代码如下：

from osgeo import gdal_array

#输入文件
src = "./thermal/thermal.tif"

# 输出文件名
tgt = "classified.jpg"

#使用gdal库加载图片到numpy
srcArr = gdal_array.LoadFile(src)

#根据类别数目将直方图分割成20个颜色区间
classes = gdal_array.numpy.histogram(srcArr, bins = 20)[1]

#颜色查找表的记录数4必须是len(classes)+1,声明RGN元组
lut =  [[255, 0, 0], [191, 48, 48], [166, 0, 0], [255, 64, 64], [255, 155, 155],
       [255, 116, 0], [191, 113, 48], [255, 178, 115], [0,