叶绿素a计算方案.docVIP

ENVI支持下利用环境小卫星数据进行水质监测方案 a、有色可溶性有机物等是水质监测的重要参数，这些水质参数浓度的变化，会引起水体生物光学特性和水面反射率的改变，利用遥感技术，根据水体光谱特性与水质参数浓度间的关系，反演水质参数，可以实现湖泊水质的高频、大范围、准实时监测。 技术流程 图1是基于环境小卫星CCD数据的叶绿素浓度反演流程图，基于其他多光谱数据（如TM等）可采用类似的流程。 采用的技术路线为：先对环境小卫星CCD数据进行数据预处理：数据读取、辐射定标、几何校正、大气校正、水域区裁剪；利用波段比值法对实测的叶绿素a浓度数据建立反演模型，将模型应用于水域区影像，反演出整个水域的叶绿素a浓度。 叶绿素a浓度反演模型有很多，如下表为前人研究的较成熟模型。 比值指数 Bnir/Bred Chla = a*X+b Chla = a*ln(X) +b Chla = a*eb*X Chla = a*Xb 归一化植被指数 X=( Bnir - Bred)/( Bnir + Bred) Chla = a*X + b Chla = a*ln(X) +b 这里使用比值指数中的Chla = a*X+b模型。 这个流程中主要使用的工具包括：ENVI环境小卫星数据处理工具包（可免费获取），波段运算，图像自动配准工具，ENVI中FLAASH大气校正扩展模块，Microsoft Excel 散点图工具等。 ? 图1 流程图 各个步骤关键技术点 一、???????? 本方案用到的环境小卫星CCD-1B数据，该数据是可以在环保部卫星环境应用中心免费下载获取（下载地址：）。根据水体中叶绿素a浓度的季节变化情况，夏季和秋季（6月-11月）水体中的叶绿素a浓度较高，应选择这个时间段内的数据。 二、???????? 网上免费获取的 HJ-1B卫星CCD的分发格式为Geotiff，每一个波段为一个Geotiff文件，并提供一个元数据说明（.XML），包含波段信息、定标参数、成像时间等。 一般的流程是依次打开每个波段，并用元数据说明中提供的定标参数（包括定标公式、参数、单位）用Band Math工具进行手动定标，再用Layer stcaking功能将4个波段的图像合成为一个多波段的数据文件。 这里使用“ENVI环境小卫星数据处理工具包”，直接读取CCD数据的.XML文件，输出结果为一个多波段的ENVI标准栅格文件，并带有中心波长等信息，可选择输出定标结果（辐射亮度数据）。工具包下载地址为：/ESRI/viewthread.php?tid=75575 三、???????? 使用TM数据为基准影像，对环境小卫星数据进行几何校正(图像配准)，目的是让野外实测数据能准确定位到影像上。 由于整景CCD数据范围非常大（400x400Km），而我们的工程区只是其中小部分，所以在进行几何配准之前，将水域区以及周边一定区域裁剪出来。同样也将TM影像裁剪相同区域。 使用ENVI下的自动图像配准工具（主菜单-Map-Registration-Automatic Registration：Image to Image），自动选择同名点完成精确配准。 四、???????? 对于环境小卫星数据，提供光谱响应函数即可FLAASH大气校正模块中进行大气校正，波谱响应函数一般是以文本格式提供，下载地址： /n16/n1115/n1522/n2118/index.html。 使用主菜单Window-Start New Plot Window，选择File-Input Data-ASCII，直接导入文本文件，得到波谱响应函数的波谱曲线格式，可用于FLAASH大气校正工具中。 五、???????? 这一步需要结合实地调查数据，将水面调查点与实测叶绿素浓度在空间上一一对应起来。实地调查数据包括：水面调查点的经纬度、叶绿素含量。以四列形式保存为.txt和Excel文件，以下格式： 序号 纬度 经度 chla实测含量（mg/L） 1 31.516334 120.182215 0.0531 2 31.457170 120.041978 0.0535 3 31.459692 120.177602 0.0496 …… 利用Bandmath工具，得到比值指数(b4/b3)。在ROI Tool中，选择ROI_Type-Input Points from ASCII，选择文本格式的实地调查数据反演点，将实测点经纬度位置以ROI方式显示在比值指数图像中，在ROI TOOL中，选择File-Output ROIs to ASCII，就很容易的将实测数据中的每一组