ENVI对国产卫星支持补丁.doc

《ENVI对国产卫星支持补丁.doc》由会员分享，可在线阅读，更多相关《ENVI对国产卫星支持补丁.doc（10页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、【精品文档】如有侵权，请联系网站删除，仅供学习与交流ENVI对国产卫星支持补丁.精品文档.ENVI5.1-IDL8.3官方更新对国产卫星的原生支持NVI5.1-IDL8.3 r6的官方更新补丁， 安装好ENVI5.1之后，解压覆盖文件即可。下载地址：下载后将后缀.exe1改成.exe注：R6补丁其他工具存在一些BUG，恢复原版文件：该补丁主要新增对国产卫星高分一号、资源一号02C、资源三号传感器的原生支持，包括数据读取、正射校正、图像融合、flaash大气校正等。另外修正了一些BUG。下面我们以高分一号影像数据为例，来说明这次更新对国产卫星数据的支持。一、数据读取在ENVI中，选择file-O

2、pen As-CRESDA-GF-1，直接选择.xml文件打开影像数据。补丁更新了ENVI可以识别的传感器列表，加入了国产卫星。如下图所示：图1 新增对国产卫星传感器的支持在Data manager中可以看到ENVI自动识别了相应的RPC文件。如下图所示。 图2 Data Manager中显示的RPC信息在Layer Manager中选中影像数据，单击右键，选择View Metadata，可以查看元数据信息，包括传感器类型、RPC信息、坐标系统、光谱信息、获取时间等，都可以自动读取。如下图所示，是多光谱图像各个波段的详细信息。 图3 Metadata查看这里读入了多光谱和全色影像。如下图所示。

3、图4 多光谱影像（上）和全色影像（下）二、正射校正ENVI自动识别了RPC文件，直接选择ENVI中的正射校正的workflow进行正射校正。由于缺少控制点信息这里直接使用无控制点正射校正功能。1. 在Toolbox中，选择Geometric Correction/Orthorectification/RPC Orthorectification Workflow，打开正射校正工作流。选择待校正影像文件。点击Next。2. 进入Orthorectification参数设置。全色图像像元大小为2，多光谱设为8；重采样选择三次卷积，其他默认。如下图所示。图5 正射校正参数（全色图像）3. 单击Fin

4、ish执行正射校正。结果如下图所示。图6 正射校正结果（左-多光谱图像，右-全色图像）三、图像融合1. 在Toolbox中，选择Image Sharpening/Gram-Schmidt Pan Sharpening，分别选择全色和多光谱数据进行融合处理，单击Next。2. 进入Pan Sharpening Parameters对话框，如下图，Sensor（传感器）选项下拉菜单中，已经添加了GF-1、ZY-3、ZY-1-02C选项。这里选择GF-1，三次卷积，选择好输出路径和文件名，单击OK。完成图像融合。图7 图像融合支持国产卫星传感器3. 融合效果如下图所示。图8 融合效果示例（上-多光谱

5、图像，下-融合结果） 四、辐射定标1. 在Toolbox中，选择Radiometric Correction/Radiometric Calibration，打开辐射定标工作流，选择多光谱影像作为Input File，点击Ok。2. 进入Radiometric Calibration对话框。Calibration Type（定标类型）有两个下拉选项：Radiance（辐射亮度）和Reflectance（反射率），这里选择Radiance；单击Apply FLAASH Setting，因为下面我们将进行Flaash大气校正，所以这里应用为Flaash大气校正对图像文件的要求：Output Int

6、erleave（输出数据存储类型）为BIL格式，Output Data Type（输出数据类型）为Float浮点型，Scale Factor（缩放尺度）为0.10；最后选择输出路径和文件名，点击OK。完成定标。图9 辐射定标参数设置3. 如下图所示，得到辐射亮度数数据。图10 辐射亮度结果五、Flaash大气校正1. 在Toolbox中，选择Radiometric Correction/Atmospheric Correction Module/FLAASH Atmospheric Correction，打开FLAASH大气校正的工作流。2. 在FLAASH大气校正参数面板，有三大部分参数，首

7、先设置第一部分：选择输入文件为经过辐射定标后的多光谱文件，选择输出路径和文件名。3. FLAASH大气校正参数面板第二部分，传感器基本信息设置： Scene Center Location（成像中心点经纬度）自动从影像中获取； Sensor Altitude（传感器高度）：在Sensor Type（传感器类型）中选择Multispectral/GF-1，可以自动获取，注：传感器类型支持国产卫星，如下图所示； Pixel Size（像元大小）：8m； Ground Elevation（成像区域平均高度）：127.95m，约等于0.128km； Flight Date（成像时间）：可以在数据上点击

8、右键，选择View Metadata，从Time中读取为2013-06-13T10:35:07Z，实际上是2013年6月13日10点35分07秒，减去8换算为GMT时间。图11 FLAASH大气校正支持国产卫星传感器4. FLAASH大气校正参数面板第三部分：大气模型和气溶胶类型的选择等。整个FLAASH大气校正参数设置如下图所示。图12 FLAASH大气校正参数 5. 单击Multispectral Setting按钮，可以看到Filter Function File已经自动读取入光谱响应曲线“gf1.sli”。6. 单击Advanced Settings，在高级设置中： Tile Size

9、 默认的是Cash size 的大小，手动改为50-100Mb（根据内存大小设定），单击OK；讨论：一般中低分辨率的影像近似天顶角（Zenith Angle）为180，方位角（Azimuth Angle）为0，即垂直观测。ENVI FLAASH的天顶角：90-180度，方位角：-180 and 180，没有查到高分一号的角度说明，这里选择默认。图13 大气校正高级参数设置7. 设置好后，在大气校正模块面板中，单击Apply。8. 大气校正完成后，检查大气校正的结果，分别加载校正前后的图像（选择CIR假彩色方式加载，可以更好的识别植被），查看典型地物的大气校正前后的光谱曲线，大致可以看出大气校正

10、之后消除了大气散射的影响。如下图所示。 图14 高分一号FAASH大气校正结果图15 大气校正前后的植被波谱曲线（上-校正前，下-校正后）ENVI下高分一号PMS相机多光谱数据大气校正(2013-11-15 11:32:04) 高分一号 PMS相机的多光谱数据拥有很好的空间分辨率（8米）和幅宽（60km），包括蓝、绿、红、近红外4个波段。本文介绍利用ENVI FLAASH工具完成高分一号PMS多光谱图像的大气校正，主要分为以下三步：注：本文是在ENVI5.0下和ENVI classic完成，其他版本基本类似。1 数据定标利用以下公式可将GF-1卫星各载荷的通道观测值计数值DN转换为卫星载荷入瞳

11、处等效表观辐亮度数据。式中：Gain为定标斜率；DN 为卫星载荷观测值； Bias为定标截距。定标参数参照表1，来源于中国资源卫星应用中心：表1 高分一号星PMS相机的定标系数卫星载荷波段号GainBiasPMS1PAN0.1886-13.127Band10.20824.6186Band20.16724.8768Band30.17484.8924Band40.1883-9.4771PMS2PAN0.1878-7.9731Band10.20727.5348Band20.17763.9395Band30.177-1.7445Band40.1909-7.2053WFV1Band10.1709-0.0

12、039Band20.1398-0.0047Band30.1195-0.0030Band40.1338-0.0274WFV2Band10.15885.5303Band20.1515-13.642Band30.1251-15.382Band40.1209-7.985WFV3Band10.155612.28Band20.1700-7.9336Band30.1392-7.031Band40.1354-4.3578WFV4Band10.18193.6469Band20.1762-13.54Band30.1463-10.998Band40.1522-12.142进行数据定标计算时候，我们可以有两种方法选择

13、使用：一、Bandmath以下是利用Bandmath工具进行定标，这种方法灵活，不过需要一个一个波段定标。1) 有了绝对定标参数和定标公式，选择Toolbox/Band Ratio/Band Math工具很容易进行传感器定标。2) 由于是单个波段文件的定标，选择Toolbox/Raster Management/Layer Stacking将定标后的单波段文件组合成一个多波段文件。二、Apply Gain and Offset工具我们可以使用Apply Gain and Offset工具。1) 启动Toolbox/Radiometric Correction/Apply Gain and Of

14、fset。2) 在Gain and Offset Values中填入 增益Gain值和偏移Offset。3) 选择输出路径和文件名执行定标。图1 Gain and Offset Values面板计算好之后，继续执行以下步骤（原数据的储存顺序是BIP）：(1) 选择Toolbox/Raster Management/ Edit ENVI Header -Edit Attributes-Wavelength，将每个波段的中心波长输入b1(502nm)，b2(576nm)，b3(680nm) ，b4(810nm)讨论：由于没有找到官方公布的中心波长信息，这里取与波谱响应中反射率为1的波长。目前为止，

15、影像数据准备工作完成了。2 波谱响应函数波段响应函数，英文名为spectral response function（SRF）或者叫Relative Spectral Response (RSR) ，与宽波段传感器出现。我们知道每一个波段都有一个波段范围，比如HJ-b1(475nm)波段为520-430纳米，实际上传感器的感光元件在这个波段范围内的每一点所感应的强度都是不一样。在成像中，原则上讲应该根据波段响应函数来进行加权平均，但由于处理起来比较麻烦，而且一般的精度要求不太高，所以大多数图像都是直接取了波段范围内的中点值来运算。波段响应函数是描述一定波长范围内（超出波段范围）的量子效应，当需要

16、精确计算像元响应时候，比如大气校正反演真实地表反射率，就需要使用波谱响应函数。任何传感器在设计时都会给出严格的波段响应函数，高分一号也不例外。在中国资源卫星应用中心可下载：将数据拷贝到文本文件中。如下图中为波谱响应函数数字表达的一部分，第一列表示波长，后面五列分别表示一个全色和4个多光谱对应波长的波谱响应值。图2 高分一号波谱响应值在ENVI中，使用波谱曲线来描述波谱响应函数，也就是以波长作为x轴，波谱响应值作为y轴，存储格式为ENVI波谱库文件(.sli)。下面介绍ENVI classic中的波谱响应函数的制作。(1) 启动ENVI Classic(2) 选择Window-Start New

17、 Plot Window，ENVI Plot Window窗口中，选择File-Input Data-ASCII,如图3所示，自动将第一列作为X轴，后面3-6列作为Y轴。单击OK。(3) 如图所示，生成了4条曲线。选择Edit-Data Parameters，更改每一条曲线的名称：b1,b2,b3,b4，便于区分。(4) 选择File- Save Plot As-Spectral Library，将波谱曲线保存为波谱库文件.sli。图3 导入ASCII文件图4 高分一号波谱响应函数的曲线表达注：如果出现以下情况，原因是没有将多光谱波段1和波段2空白部分填补0。图5 错误情况3 FLAASH大气

18、校正（1） 工具箱/Radiometric Correction/Atmospheric Correction Module/FLAASH Atmospheric Correction打开FLAASH大气校正模块；（2） 点击Input Radiance Image，前面处理好的数据，在Radiance Scale Factors面板中选择Use single scale factor for all bands，由于定标的辐射量数据与FLAASH的辐射亮度的单位相差10倍，所以在此Single scale factor选择默认：10，单击OK；注：定标后的单位是Wm2sr1m1，与FLAAS

19、H要求的单位（W）/（cm2*nm*sr）相差10倍关系，因此在Radiance Scale Factors中输入10缩放系数。（3） 设置输出文件及路径设置；（4） 传感器基本信息设置：l 成像中心点经纬度FLAASH自动从影像中获取。l 传感器高度（Sensor Altitude）：505km注：没有得到官方正式数据l 像元大小（pixel Size）：8ml 成像区域平均高度可以通过统计DEM数据获取l 成像时间：从数据头文件中读取（.XML），位置为：2013-09-08 13:14:07；减去8换算GMT时间（5） 大气模型和气溶胶模型，根据经纬度和影像区域选择。（6） 气溶胶反演方

20、法选择None（缺少短波红外），能见度设置为40km。图6 FLAASH大气校正参数设置（7） 单击Multispectral Setting按钮，在Filter Function File 导入光谱响应曲线“gf-1pms.sli”，单击OK；D:Program FilesExelisENVI51classicfilt_funcgf1.sli（8） 单击Advanced Settings，在高级设置中l Tile Size 默认的是Cash size 的大小，手动改为50-100Mb（根据内存大小设定），单击OK；l 从xml文件中得到： 97.086779.1983讨论：一般中低分辨率的影像近似天顶角：180和方位角：0，即垂直观测。ENVI FLAASH的天顶角：90-180度，方位角：-180 and 180，没有查到高分一号的角度说明，这里选择默认。图7 大气校正高级参数设置（9） 设置好后，在大气校正模块面板中，单击Apply。（10） 大气校正完成后，检查大气校正的结果，分别加载校正前后的图像（选择CIR假彩色方式加载，可以更好的识别植被），查看典型地物的大气校正前后的光谱曲线。4 总结从上面可以看到，对于未知的多光谱传感器的FLAASH大气校正，关键是获取波谱响应函数，而波谱响应函数一般都是公开的。类似的方法可用于高分一号宽幅传感器的FLAASH大气校正。