叶绿素a反演研究
【实验目的及要求】
通过实验,综合利用遥感图像处理方法进行叶绿素a反演。要求掌握辐射定标、大气校正、波谱提取、统计分析和波段运算等相关操作。
【实验环境】
ENVI5.0及以上版本或ENVI Classic
【实验过程]
(实验步骤、记录、数据、分析)
1、加载数据
加载影像数据TM5的MTL文件以及水库范围shp文件。
2、步骤
(1)辐射定标
在toolbox中点击radiometric calibration,选择MultiSpectral的数据。并进行输出设置。
(2)大气校正
在工具栏找到FLAASH大气校正,输入辐射定标文件,出现界面选择第二项并点击ok。
–选择输出路径及文件命名。根据月份及经纬度选择模型为MLS(根据影像月份及经纬度选取),以及中心高程为0.416km(根据实际影像获取)。
–在多光谱设置中,点击defaults,选择第一项,ok。
(3)点位监测及对应波谱提取
打开波谱曲线,将大伙房chla监测数据的经纬度复制在go to,并用英文逗号隔开,回车。打开新的波谱曲线,将第一个波谱的曲线拉到新的曲线图中并进行命名。依次将75个经纬度输入。最后点击Export,输出为ASCII。
(4)统计分析
将保存下来的数据进行整理,并将数据对应复制到大伙房75个chla监测数据,进行不同波段及波段组合的值和监测数据的相关性,找到一个最好的波段或组合,建立方程式。相关性最好为B4/B3,方程式为线性方程y=18.7x-16.908。
(5)范围裁剪
工具栏选择Subset Data from ROIs。选取shp数据,ROI选择yes,输入值为255,并选择输出路径及名称。
(6)波段计算
在工具栏中选择band math,将统计分析中选出的方程输入到band math中,并add to list,最后ok。出来的界面对B3、B4进行对应选择。出来的结果进行颜色显示。
【叶绿素a反演结果】
【结论】
通过利用遥感图像处理方法进行叶绿素a反演,操作有辐射定标、大气校正、波谱提取、统计分析和波段运算,通过辐射定标和大气校正将辐射亮度值或者表面反射率转换成实际的地表反射率。通过经纬度对影像波谱进行提取,并与chla数据建立模型,应用选好的模型进行波段运算进行叶绿素浓度反演。通过颜色显示可以看出叶绿素浓度高低程度,紫色显示叶绿素浓度最低,红色显示叶绿素浓度高。