165
1. 打开数据文件：datafile1=fopen('Terrain.hsi')。

2. 读入文件头的信息：cur=fread(datafile1,4,'int32')。

执行结果为（cur 的四个值：图像宽
度，图像长度，波段数，像素所占字节数）：307，500，210，2。

3. cur=fread(datafile1,inf,'int16')。

可通过”size(cur)”获得cur 的大小，这里为32235000
（=307*500*210）。

4. 将一维数组变换成210x （307*500）大小的矩阵：array=reshape(cur,210,307*500)。

Array
矩阵的每一列代表一个像素点在不同波段的值的大小，每一行为一个波段对应的全部数据。

5. 对每个波段求其标准差。

先将包含图像数据的矩阵转换成500*307)x210的矩阵，
使每一列的数据为一波段的全部数据。

1.变换矩阵
stdv=reshape(array',500*307,210);2.求方差stdv=std （stdv ）；3.画出标准差形成的曲线（图1）：
050100150200250
050
100
150
200
250
300
350
400
图1
由图像可大致分析出1-22、102-109、137-151、194-210可能为无用数据。

6. 到这里就可以用array 中的数据画出任意波段的图像。

例如，选取有用数据20个波
段的图像：a.提取该波段的全部数据并将它转成307x500的矩阵：
pic=reshape(array(175,1:500*307),307,500);b.将矩阵内的数据显示出来：imshow （pic,[]）（图2）：
图2
7. 求出它对应的直方图：hist(pic)。

（画出矩阵的直方图）如图2。

050
100
150
200
250
图3。