本文选取2001年5月20日中国北部及蒙古地区(经纬度范围:92.49°- 116.97°,33.88°- 41.23°)的一景MODIS数据进行分析,主要从读取数据、分析经纬度波段信息、第一、四、三波段融合显示、影像地理校正几方面对该景数据进行了分析。

美国ESRI公司的产品ENVI能很好地支持HDF数据格式。ENV(The Environment for Visualizing Images)遥感影像处理软件,是分析、处理并显示多光谱数据、高光谱数据和雷达数据的高级工具。ENVI能接受大量的传感器数据,是世界目前唯一能较好全面支持HDF科学数据格式的遥感影像软件。ENVI可以直接读取HDF格式(如图2所示),并能识别HDF格式中所包含的所有文件信息(如图3所示)。ENVI 打开HDF格式文件后,会自动将该数据文件所包含的所有图像信息、属性信息、文本信息作为波段列于一个波段列表中,用户可以清晰地浏览每一波段的详细信息,包括波段名称、图像波段波长、波段大小、数据类型及文件内插方式等多种信息。方便用户显示图像,并对各种属性及文本文件作各种分析。

本文选取2001年5月20日中国北部及蒙古地区(经纬度范围:92.49°- 116.97°,33.88°- 41.23°)的一景MODIS数据进行分析,主要从读取数据、分析经纬度波段信息、第一、四、三波段融合显示、影像地理校正几方面对该景数据进行了分析,具体步骤如下:

(1)、数据读取:打开ENVI,在主菜单中选择File\Open External File\Generic Formats\HDF,选择文件“MODO2QKM_03.hdf”,表示是该景MODIS数据的250米数据文件,从下图中可以看到,该文件中除两个影像波段外,还包含经度波段、纬度波段、热红外探测器的噪声信息、反射率变化参数等信息。

(2)、1、4、3波段影像融合:MODIS数据的第一、四、三波段的波段宽度分别为0.62μm ~ 0.67μm 、0.545μm ~ 0.565μm、0.459μm ~ 0.479μm,近似于可见光的红、绿、蓝波段,所以第一、四、三波段组合比较接近真彩色,故常选用这三个波段来表示MODIS影像。此处用同样方式打开500米数据文件,该文件共包含五个影像波段,将影像融合所需要的第3和第4波段进行重采样,即将其空间分辨率由500米重采样为250米,并与步骤(1)中第1波段组合,进行彩色方式显示。为提高成果影像的空间分辨率,笔者又将143波段组合影像进行对比度调整输出后,与真实空间分辨率为250米的第一波段进行影像融合(用HIS融合法),得到了几何清晰度更高的143波段融合影像(如图4所示)。图4中左侧为1、4、3波段彩色组合显示及局部放大,右侧为143波段组合输出后又与1波段进行融合的结果,可以很明显地看到,右侧的影像细节非常突出。体现了具有较高分辨率的第一波段的优势。
   (3)、影像地理校正,由于MODIS数据本身带有详细的经纬度波段信息,这种地理信息以波段的形式存放,如图5中的灰度波段所示,该灰度影像每一象素的灰度值记录的是空间分辨率为1公里的MODIS数据中对应象素点的经纬度信息,这种详细的地理信息可以使影像不需要选择大量地面控制点就可以作精纠正,而且精度会比选控制点的方法更高。ENVI软件提供了“Georeference from Input Geometry(用既定地理信息校正影像)”功能,即用现成的“IGM”( Input Geometry) 对影像进行地理校正。图5为校正前用影像关联方式在影像上查询每点的经纬度信息。
    影像地理校正可以分两步进行,首先建立GLT(Geometry Lookup File),从ENVI主菜单中选择Map\ Georeference from Input Geometry\Build GLT,然后分别按提示选择经度波段和纬度波段,在确认地图投影转换对话框后,会出现如图6所示的GLT文件参数设置:理论上“Output Pixel Size”应为250米,但由于影像边缘变形等原因,由经纬度波段计算出来的象素大小通常不是250米,而是接近250米。

“Output Rotation”是指校正后公里网格旋转的角度,目的是通过旋转公里来减少输出影像的尺寸,以减少数据量。如果希望出图时公里网仍将是正南正北的,可以将此参数改为0。在确认上述两个参数后就可以开始进行GLT计算。

Line波段和Sample波段中的每点的象素值分别表示在原始影像(校正前影像)上对应点的行列数。从图上看,影像呈网格状,用鼠标查询象素值会发现数值有正有负,如果GLT值为正,表示该点正好对应于原始影像上的某一象素点,如果为负值,则表示该点是由原始影像的象素点经过最近邻法插值得到。计算GLT(地理信息查询表)的过程比较慢,它是校正过程中最重要的一个环节。
    接着选择Map\ Georeference from Input Geometry\Georeference from GLT,选择影像文件和GLT文件,就会很快生成校正结果,如图8所示(投影类型是UTM46区)。并可在此基础上叠合经纬网、公里网和矢量信息等。
    此外,ENVI的波段运算功能可以按指定公式进行波段运算,该功能可以非常方便地用来计算植被指数等参数。但由于MODIS的植被指数计算比较复杂,至今笔者还没有从文献中找到合适的计算方法,故此处不再详述。综上所述,RSI公司的ENVI能全面支持EOS/MODIS数据,为HDF科学数据格式的管理和分析提供了最佳解决方案。

Envi处理MODIS流程相关推荐

  1. ENVI水体提取流程实践与问题

    ENVI水体提取流程实践与问题 卫片解译中,水体属于比较好判别的一类.有水体指数.阈值法.监督分类等一系列方法进行提取,今天介绍利用水体指数和ENVI提取水体的过程. 1.所用数据 Landsat8卫 ...

  2. 正弦曲线投影及ENVI下Modis数据坐标投影定义

    正弦曲线投影是一种等面积的伪圆柱投影.规定纬线投影为平行直线,经线投影为对称于中央经线的正弦曲线,同一纬线上经距相等,纬距向两极缩小.主要用于小比例尺世界地图. 在ENVI中对应的正弦曲线投影名称叫s ...

  3. 关于MODIS数据说明及简单处理

    MODIS数据产品与应用 http://www.doc88.com/p-5953424728565.html MODIS数据说明 - 我若不勇敢,谁替我坚强 - 博客频道 - CSDN.NET htt ...

  4. 使用MODIS Level 1B 1KM 数据反演AOD实验流程

    1.实验思路 实验的整体思路: 第一部分是MODIS数据预处理工作. 数据预处理工作主要包括对 MODIS LIB数据进行去除蝴蝶结效应.几何校正.波段合成.裁剪等操作 第二部分为城市地区AOD及PM ...

  5. 使用GDAL获取HDF等数据集中的图像

    使用GDAL读取HDF,NetCDF等数据集中的数据时,一般需要两个步骤,第一,获取数据集中的子数据集:第二,通过第一步获取的子数据集读取图像数据.一般的HDF图像中会有很多个子数据集,比如常用的MO ...

  6. 影像信息提取之—— 监督与非监督分类

     基于光谱的影像的分类可分为监督与非监督分类,这类分类方法适合于中低分辨率的数据,根据其原理有基于传统统计分析的.基于神经网络的.基于模式识别的等. 本专题以ENVI的监督与非监督分类的实际操作为 ...

  7. MODISL1B数据FLAASH大气校正

           MODIS数据是著名的大尺度低空间分辨率数据源,一般用于大范围监测.宏观信息的快速提取.在MODIS数据预处理时,最常做的是几何校正.对MODIS数据进行几何校正的工具很多,诸如MCTK ...

  8. IDL学习——处理自带经纬度文件的遥感影像,以哨兵5P数据为例

    IDL学习--处理自带经纬度文件的遥感影像,以哨兵5P数据为例 自带经纬度文件处理流程 1.寻找经纬度文件 2.构建GLT文件 3.使用GLT文件进行几何校正 4.批处理过程中遇到的问题 最近一直在做 ...

  9. envi反演水质参数_ENVI5.2中基于MODIS数据的海表温度反演

    劈窗算法最初是为反演海面温度开发的,具体地说是针对NOAA/AVHRR的4和5通道设计的,后来也被用来反演地表温度,这种算法较成熟,精度也高.劈窗算法以地表热辐射传导方程为基础,利用10~13μm 大 ...

最新文章

  1. oracle 行级死锁_ORACLE死锁的分类
  2. 虚拟化系列-Windows server 2012 Remote桌面与应用
  3. 深度学习核心技术精讲100篇(七十五)-集成学习
  4. 【机器视觉】 par_join算子
  5. HYSBZ - 2160 拉拉队排练(回文自动机)
  6. java单列_Java 单例模式
  7. html怎么给变量添加样式,通过CSS变量修改样式
  8. Java虚拟机(五)——程序计数器
  9. android如何不自动获取焦点,Android 如何让EditText不自动获取焦点
  10. MySQL 5.6 for Windows 解压缩版配置安装(转)
  11. HDU 1358 (所有前缀中的周期串) Period
  12. Linux下常用网络配置命令
  13. 联想计算机拆机,Lenovo(联想)Y470笔记本电脑拆机清灰图文教程
  14. 今日头条android+x86,GitHub - teajoyus/TouTiao: 今日头条 For Android
  15. graphene-django开发实践——登录为例
  16. 【转载:80个Python经典资料(教程+源码+工具)汇总】
  17. Linux网络编程 - 在服务器端运用进程间通信之管道(pipe)
  18. AD学习记录03-规则
  19. 卷积神经网络(CNN)和循环神经网络(RNN)有什么区别?
  20. 宿主软件中文版-Cakewalk Sonar Producer Edition 8.5.3 WiN

热门文章

  1. 第三周 AVI、MP4、WAV文件格式
  2. 建立基于IP地址访问的网站
  3. 手把手介绍Manjaro中添加Google输入法---亲身测试,成功安装
  4. iOS 四舍六入五成双算法
  5. 帧中继网络与NBMA/P2MA
  6. 【数据采集平台】教程-单页面采集
  7. 阿里云elasticsearch试用套路
  8. 录简写与缩略语 问题解决工具及要点
  9. 台式电脑接路由器步骤_台式电脑怎么连接无线路由器 台式电脑无线路由器
  10. 一个简易版的新闻应用(同时兼容手机和平板)