现在的位置: 首页 > 综合 > 正文

【技术类】【ArcGIS对国产卫星的支持2：高分一号卫星】篇5、实时表观反射率产品生产

2018年02月18日 ⁄ 综合 ⁄ 共 1487字 ⁄ 字号小中大 ⁄ 评论关闭

表观反射率就是指大气层顶的反射率，辐射定标的结果之一，又称视反射率。英文表示为：apparent reflectance （=地表反射率+大气反射率）

分发商提供的数据基本都是灰度值（DN值），必须转换为反射率才能进行运算（比如NDVI运算），否则得到的影像产品是不严密的。因此，我们需要在进行专题产品生产前对影像进行表观反射率产品生产。

ArcGIS内置近40种实时处理函数，通过表观反射率函数（Apparent ReflectanceFunction）可以实时表观反射率产品。

并非所有影像都可以进行此类产品生产，需具备以下条件：

- 影像采集时间和影像的太阳高度角

- 每个波段增益，偏置，以及太阳辐照度

增益和偏置是影像计算发射率的两个基本参数，如果这两个参数存在偏差，辐射定标、辐射校正及大气校正的结果都会出现错误。我们可以根据USGS公布的结果以及影像头文件获取增益和偏置值。辐照度指单位面积接收到的辐射通量。

表观反射率产品生产

第一步 展开栅格函数处理链，选中Table of Contents窗口中的影像数据->右键properties->Function，展开栅格函数处理链。

第二步 添加表观反射率函数

右键GeometricFunction->Insert-> Apparent Reflectance Function

     在添加函数之前，可以利用“像元检查”工具进行影像灰度值查看，以对比函数处理后的影像数据。
     参数设置：
     - 勾选Albedo选项，增加反照率进行表观反射率计算。
     - 各波段增益值、偏置值以及太阳辐照度已内置，无须手动输入。

第三步 删除拉伸函数

右键Stretch Funciton->Remove

此处的影像拉伸基于原始影像的统计值实现，由于影像经过表观反射率函数计算后统计值会发生改变，因此此处如果继续保留拉伸函数，会影响表观反射率函数计算后得到的影像输出结果。

      第四步     点击应用，实现反射率产品实时生成
     - 在影像分析窗口（Image Analysis）选中DRA选项，对影像进行实时拉伸显示    （此处的影像拉伸基于经过表观反射率计算后得到的影像统计值）。
     - 再次利用“像元检查”工具进行像元值查看，可以看到影像灰度值已经转换为反射率值。