数据简介
转载请在文章开头注明出处。
环境一号卫星HJ-1A星上搭载的HSI超光谱成像仪能完成对地刈宽为50公里、地面像元分辨率为100米、光谱范围（0.45-0.95µm）和115个光谱谱段的推扫成像。
操作软件ENVI5.3.1：https://www.ixxin.cn/2016/12/14/envi-5-3-sp1pojie/
预处理步骤
读取影像
ENVI5.3.1直接open HSI数据的.H5文件即可。
点击ImageData，双击Banddata即可加载到右边。
同时获取以下属性信息：
以下信息用于FLAASH：
Scalefactor：100（定标值系数）
MapInformation→SceneCenterLat（纬度） SceneCenterLong（经度）
MapInformation→PixelSpacing（像元大小）
ImageAttributes→SceneStartTime（数据获取时间）
以下信息用于几何配准：
MapInformation:
Projection(投影)
Utm-zone(带号)
ProductUpperLeftX(左上角X坐标)
ProductUpperLeftY(左上角Y坐标)
单击Open Rasters按钮，打开图像文件。
点击File——save as 保存为ENVI标准格式。
几何配准
需要在Classic中进行。
加载上一步保存的影像。
右键数据，选择Edit Header。
点击Edit Attributs。
选择Map info
依次选择投影信息，椭球体：WGS84，填写带号，填写左上角坐标，填写像元大小100.（所有参数上一步均已获得）
条带去除
HSI前面波段数据质量很差，有很多条带，如图：
第一步，将条带旋转到垂直状态：
用到工具：rotate/flip data
角度为：-14.39
第二步去除条带：
用到的工具：spear vertical removal
参数设置为100
去条带后对比：
第三步，旋转到原始状态
用到工具：rotate/flip data
角度为：14.39
FLAASH大气校正
第一步，数据转换
FLAASH需要BIL存储方式
在Toolbox工具箱中，双击Raster Management→Convert Interleave工具，选择刚打开的图像文件，在Convert File Parameters中，设置Output Interleave：BIL，选择输出文件名“HSI_BIL.dat”和路径，单击OK按钮。
第二步：FLAASH大气校正
（1）在Toolbox工具箱中，双击Radiometric Correction/Atmospheric Correction Module /FLAASH Atmospheric Correction工具，单击Input Radiance Image按钮，选择上一步保存的“HSI_BIL.dat”图像。
（2）在打开的Enter ASCII Filename containing wavelengths and fwhm对话框中，选择含有波段信息的txt文件HJ1A-HSI-BandInfo.txt，并在Input ASCII Filename对话框中手动输入Wavelength和FWHM的行数，如图所示。
提示：HJ1A-HIS数据的文件夹中包含一个文本文件，如本示例中的“HJ1A-HSI-2-65-A2-20110820-L20000597997_BandInfo.txt”，其中包括了每个波段的波谱范围和波段宽度。“HJ1A-HSI-BandInfo.txt”就是基于这个文件生成的，即第一列和第二例的平均值作为中心波长（wavelenth），第三列作为波段宽度（FWHM）。“HJ1A-HSI-BandInfo.txt”下载地址：http://pan.baidu.com/s/1o64sVzG
（3）在出现的Radiance Scale Factors对话中选中Use single scale factor for all bands，缩放系数填写：1000。
提示：下载的HJ1A-HIS数据已经经过了定标，为了减少储存空间而扩大了100倍。单位是100* W/(m2 *µm * sr)。
（4）单击Output Reflectance File按钮，选择输出文件名和路径。
（5）单击Output Directory for FLAASH Files按钮，选择文件输出目录。
（6）传感器与图像信息
lLat：40.4616° Lon：116.195°
l传感器类型（Sensor Type）：UNKNOWN-HIS
l传感器飞行高度Sensor Altitude：650
l图像区域平均海拔（Ground Elevation）：0.5
l成像时间（Flight Date）：2011-08-20 Flight Time：03:03:47.96。
（7）大气参数
l模型（Atmospheric Model）：Mid-Latitude Summer
提示：根据图像中心纬度40°和8月份成像从帮助表中查询获取。
l执行水汽反演（Water Retrieval）：Yes。
lWater Absorption Feature：820nm
l气溶胶模型（Aerosol Model）：Urban
l气溶胶反演（Aerosol Retrieval）：None
提示：HJ1A-HIS数据缺少短波红外波段，采用能见度值参与气溶胶计算。
（8） 其他按照默认：
• 初始能见度（Initial Visibility）：30。
• Spectral Polishing：yes
• Width（number of bands）：5
• Wavelength Recalibration：No
（9） 单击高光谱设置（Hyperspectral Settings）按钮，选择默认设置。
（10）单击Advanced Settings按钮，打开高级设置面板。修改以下几个参数，其他按照默认：Tile Size：100
（11）参数如图所示。单击Apply按钮，执行FLAASH处理。
结果：
后语
参考文章：http://blog.sina.com.cn/s/blog_764b1e9d0102uxde.html
打赏直接扫描二维码，打赏名单请看：https://www.ixxin.cn/about.html
