envi遥感图像处理之大气纠正.docVIP

大气校正说明文档 步骤一：辐射定标 本实验采用的是绝对辐射定标，直接建立遥感影像DN值与接收到的能量的之间的关系。 建立关系所采用的公式是： 其中， ， Lmaxλ和Lminλ通过参看遥感影像的头文件进行确定。fullDNrange取的是255。 具体操作如下： 打开遥感影像文件及其头文件 根据头文件信息计算gain和offset的值 在envi的Basic Tools中打开 Band Math像，将本步骤采用的公式写入band math中，计算出L。 至此，就完成了遥感影像的辐射定标过程。 步骤二：未进行大气校正所得到的反射率 本步骤讲述如何从经过辐射定标的遥感影像直接生成地物的反射率的影像，制作该影像的目的是为了与后面经过大气校正后的影像进行对比。 本步骤所采用的公式是： 其中，L是由上步所算出来的，d指的是实际的日地距离，单位是天文距离，ESUN指的是太阳平均辐射强度，θ为太阳天顶角。 d值可以由观测时间查阅相关资料获得。ESUN值也可以由相关资料获取。可以从头文件中获得。 具体操作如下： 查阅相关资料，确定参数θ、d、ESUN 在envi的Basic Tools中打开 Band Math像，将本步骤采用的公式写入band math中 确定变量b2为上步所算的L，并由此计算出未进行大气校正的反射率。 由此，我们就得到了未经大气校正的反射率。 步骤三：进行大气校正，得到地物反射率 由于大气的影响，会使得遥感影像的反射率发生较大的变化，为了得到地表的真实的反射情况，我们就需要对所得的遥感影像进行大气校正。 本步骤采用的公式是： 和以及是由6S模型经过计算得到的，是在第一步中计算得到的。 具体操作步骤如下： 确定6S模型的参数。 Table 5 6s Parameters 输入值 含义 注释 卫星参数 7 Landsat卫星 卫星参数 11 22 2.5 116.62 23.11 month, day, hour, long, lat 成像月份、日期、时间（格林威治时间）、图像中心点精度、纬度 季节特征 3 Midlatitude winter 中纬度冬季 气溶胶类型 1 continental 陆地 可视距离 23 Visibility, km 能见度（晴朗天气） 目标物高度 -0.1 Target altitude in km 表示位于海平面 卫星参数 -1000 Satellite case 传感器搭载于卫星上 卫星波段 参数 28（band4） Spectral band definition TM的波段取值： 25: band1 26: band2 27: band3 28: band4 29: band 5 30: band7 不包括热红外波段 地表性质 0 Homogeneous case 均一地表 BRDF影响 0 No BRDF effect 没有BRDF效应 目标特性 0 Uniform target 目标物均一 图像反射率 0.5 Reflectance = 0.5 图像辐亮度 100 根据此表，确定6S模型在各个波段上的参数。 在DOS窗口中运行6S模型程序，并依次输入上步所确定的参数 由6S模型得到参数和以及 在envi的Basic Tools中打开 Band Math像，将本步骤采用的公式写入band math中，由此计算出该遥感影像经过大气校正以后的反射率。 这样，我们就得到了经过大气校正的地表的反射率了。 步骤四：将大气校正前后的遥感影像进行对比 本步骤是让我们将遥感影像的前后进行对照，来看看进行大气校正前后遥感影像之间的差异。 本步骤所采用的方法是用未校正的反射率减去校正后的反射率，得到一幅新的影像，从新的影像上去看二者之间的差异。 具体步骤是：在envi的Basic Tools中打开 Band Math，在band math中输入：abs(b1-b2)。然后得到了新的图像，从新的图像上来看二者之间的差异。 至此，我们完成了遥感影像的辐射定标和计算未进行大气校正的地物的反射率和进行大气校正了的地物的反射率以及二者之间的差异。 本实验有关的数据： 使用的数据为广东省汕头市的ETM+影像，成像时间为2001年11月22日，2:28:18.000（格林威治时间）。 波段列表： (um) 蓝波段 30m TM2：0.52~0.60 (um) 绿波段 30m TM3：0.63~0.69 (um) 红波段 30m TM4：0.76~0.90 (um) 近红外波段 30m TM5：1.55~1.75 (um) 短波红外波段 30m TM7：2.08~2.35 (um) 短波红外波段 30m ETM+