水环境遥感3教学课件.ppt

导出点的反演值 这样就将水面调查点与BNIR/BRED对应的值导出来了，如下： 我们会发现，导出来的经纬度与输入实测点的经纬度不完全一致，这是因为一般影像中像元的坐标是取中心点的经纬度，而我们实测的点位不一定恰好对应着影像像元的中心点，所以当输入的经纬度与影像上单个像元的经纬度不一致时，就会采用就近原则，与最临近的像素点匹配，输出该点的经纬度坐标。环境小卫星的空间分辨率是30米，换算成经纬度在几秒的范围。 将相应的BNIR/BRED值导入Excel表中，与实测值一一对应。 第三步：模型参数反演 （1）有了实测的叶绿素a含量值与BNIR/BRED值的对应表，在 Excel 中选BNIR/BRED值与叶绿素 a 实测值绘制散点图。 在 Excel 中选BNIR/BRED值与叶绿素a实测值绘制散点图 在散点图上选中散点，单击右键-添加趋势线，打开设置趋势线格式面板，勾选线性，显示公式，显示R平方值。点击“关闭”按钮，线性回归方程和R方值在散点图上显示。 得到反演模型为： y=0.0422x + 0.0171，R2=0.7929。 第四步：叶绿素反演 有了参数a和b，反演模型可表达为： Chla= 0.0422* BNIR/BRED +0.0171 使用Basic Tool-Band math，在 Enter an expression下面输入表达式： 0.0422*b1 + 0.0171，单击 Add to List，单击 OK，在Variables to Bands Pairings面板中，选择b1为 b4/b3的图像，设置输出路径和文件名，单击 OK，计算得到叶绿素反演图。 成果后期处理与应用 反演结果输出 在display 中，选择Tools-Color Mapping-ENVI Color Tables，将结果以彩色方式显示。 打开Overlay-Annotation 注记面板，打开工具 Object-Color Ramp，鼠标在图中空白处点击，添加一个色带图例。 叶绿素反演结果 可以看到，利用遥感手段可以快速、大范围的获取水质参数信息，除了能反演出水质参数含量外，还可以获取空间分布信息。同时，环境小卫星2天的重返周期以及可以免费下载数据，可以满足遥感在水质参数监测中的业务化应用要求。 第三步，几何校正（图像配准） 下面以TM作为基准影像对环境小卫星图像进行图像配准。 （1）打开基准影像TM_baseimage.img。 （2）选择主菜单- Map- Registration -Automatic Registration：Image to Image，选择基准影像（TM_baseimage.img）。 （3）选择被配准影像 HJ1B-CCD1Cal-sub.img。 （4）选择控制点，并调整控制点，直到总的RMS Error小于1个像素时，完成控制点的选择。点击Ground Control Points Selection上的File-Save GCPs to ASCII，保存控制点。 （5）在Ground Control Points Selection面板上，选择Options-Warp File （As Image Map），选择校正文 件。在校正参数面板中，投影参数默认， （6）在X和Y的像元大小输入30米，按回车，图像输出大小自动更改重采样方法选择Nearest Neighbor（保持图像辐射值不变） （7）选择输出路径和文件名HJ1B-CCD1Cal-sub-jz.img，单击OK 按钮 几何校正输出参数设置 （8）打开校正后的影像和基准影像，在显示校正后影像的窗口中，右键选择Geographic Link 命令，选择显示图像的两个窗口，打开十字光标查看校正结果 。 第四步：大气校正 环境小卫星提供了波谱响应函数，以文本形式提供，第一列表示波长（nm），后面四列分别表示 4 个波段对应波长的波谱响应值。 为了做大气校正，需要制作波谱曲线来描述波谱响应函数，用于大气校正。 注：波谱响应函数与CCD型号相对应，这里的 CCD类型是环境 1B 星CCD1（可以从文件名中看出） 一、 制作波谱曲线 （1）主菜单 Window-Start New Plot Window，打开 ENVI Plot Window 面板，在波谱绘制窗口中，选择 File-Input Data-ASCII，导入“环境1B星 CCD1 光谱响应曲线.txt”文本文件，如图，自动将第一列作为x轴，后面4列作为y轴，波长单位选择Nanometers，单击OK。 制作波谱响应函数文件 （2）如图，在绘制窗口生成了4条曲线，选择 Edit-Data Para