高光谱图像使用kmean算法进行无监督分类

from spectral import *
import matplotlib.pyplot as plt

#读取文件，未加载到内存中，使用索引后再读取索引数据到内存中
data=envi.open(r"20230518123937/20230518123937.hdr",r"20230518123937/20230518123937.os")

#将整个图片数据加载到内存中，默认使用bip格式，即R-行，C-列，B-波段
bip=data.open_memmap()

#截取部分区域，第3001行到3300行和第1列到300列的区域
bip_pixels=bip[3000:3300,0:300,:]

#k近邻聚类，（图片数据，类别数，最大迭代次数）
(m, c) = kmeans(bip_pixels,5, 300000)
plt.figure()
for i in range(c.shape[0]):
    plt.plot(c[i])
plt.grid()
plt.show()