遥感水体反演

基本步骤：数据读取-辐射 定标-几何校正-大气校正（波普响应函数）

1、搭建环境，将拓展工具（ENVI_HJ1A1B_Tools.sav）复制到安装目录extension下；

2、数据下载：中国资源卫星应用中心

详细步骤：

1、双击打开“ENVI_HJ1A1B_Tools.sav”→选择“ccd(多光谱数据)”→点击选择“input path”选择多光谱数据文件→点击search→选择该文件→output→选择calibration→apply；

2、裁剪相关区域：打开envi classic→windows→available fist list→file→open new file→（选择hj1b-ccd1-20091006-cal-sub.img)→选择图像右键→选择load true color→选择classic界面下的file→save as→ envi standard→import→spatial subset→image→选择研究区域→OK→OK→choose命名”----subset“

3、几何校正：打开envi →打开”——subset“和 “TM-baseimage.img”→”geometric  correction"→registration→image to image→选择 TM-baseimage band4→OK→——subset band4→否→OK→options→clear→选择控制点（20个左右）→在options→file warp→选择”----subset“→choose“----subset-jz“

4、打开”681_HJ1BCCD1.txt"在excel中粘贴为数值，并出入折线图，然后插入相关公式

5、打开classic→window→start new plot window→file →input data→ASCII→波普响应函数txt→column2-5→wavelength:nanometer→OK→file→save plot →spectral library（波谱库）→select all items→OK→choose→hj1b-ccd1-lib.

6、数据格式转换：打开“----subset-jz“→raster management→convert interleave→选择“BIP\BIL"→choose“----subset-jz-BIP\BIL“

7、大气校正：打开下载数据目录下的XML元数据查看里面的拍摄时间及中央经纬线计算出度分秒

在工具箱radiometric correction →FLAASH atmospheric correction→ input”----subset-jz-BIP\BIL“→use single →factor“10”→output“”----subset-jz-BIP\BIL-flash“→分别输入lat\lan,altitude ,ground elevation ,flight time等数据→点击底部的multispectral setting→file function file（过滤文件）→hj1b-ccd1-lib. →OK→default→”460：2100nm→wrap band1→OK

8.链接视图

在envi的工具栏中选择“view”→create→打开data manager→“----subset-jz-BIP\BIL“打开数据格式转换过的文件→右击选择”true color“→点击视图的另一个视图框→导入”----subset-jz-BIP\BIL-flash“大气校正后的数据→右击”true color“→”view“→”选择“link view"→ link all→OK→display→profiles→spectral

9.叶绿素反演

在工具箱中打开”band math “→输入公式”float（b4)/b3→add to list→OK→选择大气校正的b3和b4波段→choose“b3b4”

10.采样点的加载

打开classic经典界面→window→available band list→打开"b3b4"数据→在load band界面→选择“overlay”→region of interest →ROI-TYPE→input points from ASCII→打开实测数据中的反演点文件→在x points “3”→comprise→选择“individual”→选择“geographic lat/lon”投影→OK→file→output rois to ascii →b3b4→选择edit output ascii→选择(1,4.5打勾）→保存“反演点”用excel处理出反演点的公式

11、建立mask

打开数据“b3b4”和taihu-vector.shp数据→在工具箱中打开”build mask“→选择”b3b4"→choose→保存“taihu-mask"→在工具箱中 选择：”apply mask“→”b3b4”→select band mask→保存“taihu-mask-caijian→忽略零值

12.打开"band math"→输入用excel处理出反演点的公式→add to list→OK→b3b4→choose→taihu-result→OK→忽略零值