modis terra卫星过境时间_好文速递：IDL程序去除MODIS遥感影像镶嵌图中的黑色条带...

采用基于地理坐标的镶嵌方法和剪切线的方法对两幅MODIS遥感影像做镶嵌时, 在影像接缝处出现了黑色条带, 条带值显示为-NaN或者0。主要针对这一现象介绍在IDL7.0环境下编写的Delete NaN程序和Delete 0程序的应用。该程序实现了消除黑色条带, 并找回了丢失的数据, 该程序运行速度快, 操作简便, 特别适用于MODIS遥感数据的影像镶嵌图。

1 IDL编程语言

IDL( InteractiveDataLanguage)交互式数据语言是面向矩阵的且用于数据可视化研究与应用开发的第四代计算机语言 是进行二维及多维数据可视化分析及应用开发的理想软件工具。IDL致力于科学数据的可视化和分析, 是跨平台应用开发的最佳选择。它集可视、交互分析、大型商业开发为一体, 为用户提供了完善、灵活、有效的开发环境。

IDL具有快速分析超大规模数据的能力, 具有高级图像处理能力、交互式二维和三维图形技术、面向对象的编程方式、OpenGL图形加速功能、集成数学分析与统计软件包、完善的信号处理和图像处理功能、灵活的数据输入输出方式、跨平台图形用户界面工具包、连接ODBC兼容数据库及多种外部程序连接工具。

IDL还为用户提供了可视化数据分析的解决方案, 既可以让科学研究人员交互式浏览和分析数据, 又为程序员提供了快速程序原型开发并跨平台发布的高级编程工具。IDL使科学家无需编写大量的传统程序就可直接研究数据。

随着MODIS遥感数据的广泛应用, 为了快速、精确、简便地获取目标区域的MODIS遥感数据, 经过预处理后的MODIS遥感图像需要进行镶嵌。在基于地理坐标的镶嵌中, 使用剪切线方法, 两幅图像接缝处出现了黑色条带, 且黑色条带处的数据丢失并显示值为-NaN或者0 。为了消除黑色条带, 并找回丢失的数据, 我们编写了IDL程序( Delete NaN.pro和Delete 0.pro) , 该程序在IDL7.0环境下运行, 效果良好。本文介绍该程序及该程序在利用MODIS数据反演的比值植被指数RVI、亮温温度的镶嵌图中的应用。

2 程序Delete NaN的代码实现与结果保存

该程序使用了替代方法, 把一幅图中的-NaN值用另一幅图中对应位置的数据替换, 替换后得到的图像数据输出到计算机D盘目录下, 并以文件名Data1.dat保存。具体代码如下:index=WHERE( ～ FINITE( Data1) )Data1[ index] =Data2[ index] //用遥感图像2中的数据替换遥感图像1中的-NaN值My OutFileName=`D: Data1.dat' //替换后的遥感图像以路径D: Data1.dat保存。 3 程序Delete 0的代码实现与结果保存该程序也使用了替代方法, 主要功能是消除镶嵌图中黑色条带处的0值, 输出结果以文件名Data2.dat。保存在D盘。具体代码如下:index=WHERE( Data1 eq0)Data1[ index] =Data2[ index] //用遥感图像2中的数据替换遥感图像1中的0值。My OutFileName=`D: Data2.dat' //替换后的遥感图像以路径D: Data2.dat保存。 4 程序Delete NaN的应用

以TERRA卫星2006年3月2日凌晨2点45分和2点50分的台湾地区的MODIS数据为例来说明程序Delete NaN的应用。本文通过MODIS第1和2波段来计算比值植被指数RVI, 对于MODIS图像中的任何一个像元。

把经过RVI运算得到的两幅RVI图像RVI 02 45和RVI 02 50进行两次镶嵌, 镶嵌时不使用剪切线, 得到两幅镶嵌图 :一次是RVI 02 45在上面, 镶嵌后得到图像RVI 1;另一次是RVI 02 50在上面, 镶嵌后得到图像RVI 2。然后在IDL7.0环境下运行程序Delete NaN.pro, 显现打开第一幅图像的对话框, 选择镶嵌后的图像RVI 1;然后又出现了打开第二幅图像的对话框, 选择图像RVI 2, 运行结束后,在计算机的D盘根目录下生成了Data1.dat和Data1.hdr两个文件。用ENVI4.5遥感图像处理软件打开Data1.dat文件, 可以看到镶嵌图中原是黑色条带的地方, 黑色条带被消除, 黑色条带区域原是-NaN值的地方也找回了丢失值。

5 程序Delete 0的应用

在使用TERRA卫星2008年8月30日凌晨4点25分和4点30分的四川攀枝花地区的MODIS遥感数据反演的第31波段的亮温温度的镶嵌图中, 程序Delete 0也得到了很好的应用, 过程同上, 结果见图4中的T31。第31波段的亮温温度用下式计算

式中:I31是MODIS第31波段的热辐射强度( W· m-2 · μm-1 · sr-1 ), T31的单位是开尔文( K), I31用下式计算：

 I31 =radiance scalesB ( SIB.T.FS-radiance offsetsB )

式中:SIB.T.FS是记录值, 定标参数radiance scalesB和radiance offsetsB与第31波段有关, 可以使用HDFexplorer软件或者HDFview软件查取。