“4-专题：基于环境小卫星的湖泊水质遥感监测\3-大气校正\太 湖矢量数据”
3、叶绿素反演
叶绿素反演
本环节最重要的是地面实测点与星上 点的反演模型建立，涉及到定量遥感 的知识。 在前人研究了大量的算法和模型的基 础上，本专题选择经验模型之一波段 值模型（BNIR/BRED）来进行叶绿素 的反演，做为示例。 波段比值模型：
度。数据获ຫໍສະໝຸດ 说明 环境小卫星CCD-1B数据，在环保部卫星环境应用中心免费下 载获取，环境小卫星CCD数据下载地址： 环境小卫星的读取补丁下载地址： /s/blog_764b1e9d0101a9dh.html TM数据下载地址为： 大气校正用到的环境小卫星数据波谱响应函数下载地址： /n16/n1115/n1522/n2118/index. html
制图输出
启动32位envi，选择菜单File→Print，调用ArcGIS的制图工具。 如果没有安装ArcGIS，可以使用ENVI Classic的QuickMap。
专题总结
该专题从原始的 HJ1B-CCD 数据开始，在 ENVI 下经过一系列 数据预处理，得到具有正确地理坐标和真实光谱信息的太湖水 面数据，结合数学模型和定量遥感技术，反演太湖叶绿素浓度 。得到太湖叶绿素浓度成果图。 掌握知识：
反演出的参数a和b，带入反演模型：
Chla= 0.0417* BNIR/BRED+0.0167 叶绿素反演 Basic Tool->Bandmath， 表达式：0.0417*b1 + 0.0167 设置b1为b4/b3图像
得到叶绿素反演图。
4、成果后处理及应用
4.1 反演结果验证
2.1 数据读取和定标
安装环境小卫星数据读取和定标补丁
Home\Program Files\Exelis\ENVI51\Extensions目录下
数据读取和定标
主菜单->File->Open As->HJ-1A1B
2.2 几何校正（一、目的）
这一步是对环境数据进行几何校正，使其具有精确的地理信息 ，基准影像可以是地形图，也可以是已经过校正的其他中高分 辨率影像，本专题中采用已经过校正的TM数据作为基准影像。 目的是让环境小卫星影像拥有精确的坐标，可以与其他数据向 匹配，如实测GPS数据。
波谱范围：400－2500nm
2.3 大气校正（四、FLAASH基本参数设置）
基本参数设置
2.3 大气校正（五、FLAASH多光谱参数设置）
多光谱参数设置
2.3 大气校正（六、FLAASH高级参数设置）
高级参数设置
2.3 大气校正（七、结果检验）
结果检验 大气校正前后典型地物（植被）光谱曲线对比
将验证点对应的叶绿素反演值导出来，与验证点的实测值在 Excel表中一一对应。 同样可以进行反演模型建立，以R2对进行结果验证。
4.2 反演结果制图输出
彩色显示
在Layer Manager中，右键选择Raster Color Slices，在弹出的 对话框中设置颜色表，单击OK,将结果以彩色方式显示。
专题：基于环境小卫星的湖泊 水质遥感监测
技术支持邮箱：ENVI-IDL@ 技术支持热线：400-819-2881-7 官方技术博客：/enviidl 官方技术qq群：148564800
1、专题背景
专题背景
太湖是中国第三大淡水湖， 也是富营养化比较严重的湖 泊之一。悬浮物、叶绿素 a 、有色可溶性有机物等是水 质监测的重要参数。
Chla=a*（BNIR/BRED）+b
模型应用于图像
模型修正 星上点与地面实测点建模
3.1 采样点实测数据整理
这一步需要整理实地调查数据，将水面调查点的空间位置及实 测chla值一一对应。 实测点归类——反演点、验证点。 按照标准格式：点号、水面调查点的经纬度、叶绿素含量，以 四列形式保存为.txt和Excel文件。
ENVI主模块 FLAASH大气校正扩展模块 IDL开发的环境小卫星数据读取补丁 Excel相关功能
处理流程
数 据 预 处 理
HJ1B-CCD1原始数据
数据读取 辐射定标 几何校正 大气校正 太湖区裁剪
环境小卫星数据 处理补丁
基准影像 波谱响应函数 太湖矢量数据
星上点与地面实测点建模 叶 绿 素 浓 度 反 演
国产卫星数据的处理 预处理流程的安排 数据处理各个环节在ENVI下的操作及参数选择 定量遥感技术实现
2.2 几何校正 (二、工程区裁剪）
由于整景数据范围非常大，而我们的工程区只是其中一小部分 ，所以在进行几何配准之前，将太湖以及周边区域裁剪出来。
2.2 几何校正（三、图像配准）
以TM作为基准影像对环境小卫星图像进行图像配准
2.3 大气校正（一、概述）
获得地物真实的光谱信息是遥感反演的前提，大气校正可消除 大气影响，还原地物的真实信息，是定量遥感数据预处理中必 不可少的环节。
2.3 大气校正（二、制作波谱响应曲线）
环境小卫星提供了波谱响应函数，以文本形式提供，第一列表 示波长（ nm ），后面四列分别表示 4 个波段对应波长的波谱 响应值。需要制作波谱曲线来描述波谱响应函数，用于大气校 正
2.3 大气校正——（三、FLAASH大气校正）
数据准备 FLAASH对图像文件有以下几个要求：

这些水质参数浓度的变化，会引起水体生物光学特性和水面反射率的改 变，利用遥感技术，根据水体光谱特性与水质参数浓度间的关系，反演 水质参数，可以实现湖泊水质的高频率、大范围、准实时监测
专题概述
本专题详细介绍了利用环境小卫星 CCD-1B图像反演太湖叶绿 素a浓度的完整流程，涉及环境小卫星的数据读取、辐射定标 、几何校正、大气校正、反演模型的建立、遥感反演过程、反 演结果验证等内容。 使用模块和工具：
2.4 裁剪太湖区
利 用 ENVI FX 提 供 的 面 向 对 象 图 像 分 割 工 具 ， 首 先 获 得 Shapefile，然后利用ENVI的矢量编辑功能获取太湖矢量文件。 保存在 “4-专题：基于环境小卫星的湖泊水质遥感监测\3-大气 校正\太湖矢量数据”
2.4 裁剪太湖区
利用获得的Shapefile进行图像裁剪 太湖矢量数据：
数据是经过定标后的辐射亮度（辐射率）数据，单位是： （μ W）/（cm2*nm*sr）。 数据带有中心波长（wavelenth）值，如果是高光谱还必须有波段宽度（
FWHM）,这两个参数都可以通过编辑头文件信息输入（Edit Header）。
数据类型
• 支持四种数据类型：浮点型（ floating ）、长整型 (long integer ) 、整型（ integer）和无符号整型 (unsigned int)。数据存储类型： ENVI标准栅格格式 文件，且是BIP或者BIL。
模型修正
模型应用于图像
叶绿素浓度分布图
结果分析
流程说明
根据环境小卫星CCD数据特点及太湖水质反演技术要求，采用 的技术路线为：先对环境小卫星CCD数据进行数据预处理：数 据读取、辐射定标、几何校正、大气校正、太湖区裁剪， 利 用波段比值法对实测的叶绿素 a浓度数据建立反演模型，将模 型应用于太湖水面区域影像，反演出整个太湖区的叶绿素 a浓
3.2 采样点星上数据获取
生成波段比值影像
Basic Tool->Bandmath，表达式：float(b4)/b3
在波段比值影像上导入采样点数据，将影像上对应位置的像素 值导出
3.3 模型参数反演
在Excel中选中BNIR/BRED值与叶绿素a实测值，绘制散点图
3.4 叶绿素反演
进入专题
数据：
“4-专题：基于环境小卫星的湖泊水质遥感监测”
成果后期 叶绿素反演 处理及应 用
数据预处理
2、数据预处理
图像预处理流程
HJ1B-CCD1原始数据
数据读取 环境小卫星数据处理补丁 辐射定标 几何校正 大气校正 太湖区裁剪 基准影像 波谱响应函数
太湖矢量数据
太湖水面HJ1B-CCD1数据