ENVI中打开MODIS数据及简单处理
- 格式:doc
- 大小:12.00 KB
- 文档页数:1
Envi调用MODISReprojectionTool对MODIS产品进行批处理拼接1熟悉MRTMODIS产品的类型不同,一景HDF格式的影像包含的波段也各不相同。
MRT处理时需要选择处理波段,0表示不作处理,1表示处理,第一要确定影像的波段数。
1、拼接Mrtmosaic.exe程序用来拼接影像。
调用方式为:mrtmosaic -i "g:\n%1.txt" -s " 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 " -o "g:\temp.hdf"命令行关心如下Usage: mrtmosaic -i input_filenames_file -t -h -o output_filename-s spectral_subset "b1 b2 ... bN"-g filename for the log file需要指定三个参数:-I 指定输入文件,能够采取两种方式1、1.hdf 2.hdf 3.hdf2、把影像的完整路径储存到txt文件中,作为input参数。
如-I MOD092018001.txt。
-s 指定需要处理的波段,同样能够采纳两种方式1、直截了当给出,如–s “0 0 0 0 0 0 0 0 00 0 0”,注意英文引号;2、指定一个txt路径,让程序读取;-o 指定输出路径,一样直截了当给出1、直截了当给出,如–o g:\tmp.hdf。
注意直截了当存为HDF格式,便于后续处理。
2、重投影、裁切Resample.exe用来重投影、裁切是MRT程序的核心。
调用方式为:命令行关心如下:Usage: resample -p parameter_file [options]Options that override parameter file specifications:-i input_file_name-o output_file_name-r resampling_type [NN BI CC NONE]-t projection_type [AEA ER GEO HAM IGH ISIN LA LCC MERCAT MOL PS SIN TM UTM]-j projection_parameter_list "p1 p2 (15)-s spectral_subset "b1 b2 ... bN"If using the -s switch, the SDSs should be represented as anarray of 0s and 1s. A '1' specifies to process that SDS;'0' specifies to skip that SDS. Unspecified SDSs will not be processed.If the -s switch is not specified, then all SDSs will be processed.-a spatial_subset_type [INPUT_LAT_LONG INPUT_LINE_SAMPLE OUTPUT_PROJ_COORDS]-l spatial_subset "ULlat ULlong LRlat LRlong"-or- "ULline ULsample LRline LRsample (0-based)"-or- "ULprojx ULprojy LRprojx LRprojy"NOTE: line/sample must be specified for the highest resolution of all SDSs specified to be processed in the product.-u UTM_zone-x pixel_size-g filename for the log file能够只指定1个参数:用–p 读入prm参数文件,进行处理。
美国RSI公司(Research Systems Inc.)的产品ENVI能很好地支持HDF数据格式。
ENVI (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影像。
利用ENVI处理Modis数据的方法
1、对于hdf格式modis数据,首先打开hkm(500m)数据,这个数据有两个数据组(反射率与辐射率,要用到辐射率数据),都有1-7这样7个波段。
再打开qkm(250m)数据,这个数据有两个数据组(反射率与辐射率,要用到辐射率数据),都有1-2这样2个波段。
2、数据制作时真彩色数据处理波段组合:
qkm band1 250m
hkm band4 500m
hkm band3 500m
3、需要首先将modis两个辐射率数据组(hkm500m辐射率数据和qkm250m辐射率数据),进行自动纠正,map/geoference modis ,同时系统默认会做双眼皮除去处理,
hkm500m辐射率数据做纠正时需要可以选择143波谱子集,并重采样成250米。
生成hkm-143jz文件
qkm500m辐射率数据做纠正时可以选择1波谱子集,生成qkmjz文件
4、再对纠正后的数据裁成两幅一样大小的数据
5、波段组合:
qkmjz band1 250m
hkm-143jz band4 250m
hkm-143jz band3 250m
效果图如右图所示:。
一般处理步骤:
1几何校正
Map----georeference MODIS--------选择波段--------选择投影方式(经纬度)------OK(选择辐射率数据,波段1-7)
2.拼接
3.大气校正(flaash)
用envi自带的flaash模块做大气校正,注意,输入文件的要求:1)数据是经过定标后的辐射率数据,2)数据要有中心波长值,可通过edit header输入,3)支持四种数据类型:浮点型等等,4)数据储存类型:标准栅格文件,且必须是BIP或BIL,5)适用的数据光谱范围是0.4-2500μm。
1.1打开经过几何校正的数据,输入中心波段,直接导入输好的txt文件;
1.2修改数据类型,点开basic tools---convert data把数据类型改成BIL,选择
in place直接覆盖原文件;
1.3点开Spectral-----Flaash做大气校正,输入辐射率文件,
2创建ROI区域
3导出所需数据。
遥感图像处理系统ENVI及其在MODIS数据处理中的应用摘要:介绍了美国遥感图像处理系统ENVI的主要功能,包括数据接口、交互式分析功能、波谱工具、矢量处理功能、及遥感图像处理功能等。
并将MODIS 影像(HDF数据格式)在ENVI中进行了处理分析,包括数据读入、影像几何校正、监督分类和制图。
应用证明,ENVI可以处理最先进的卫星数据格式,在高级遥感图像处理和高光谱影像处理方面具有优势。
关键词:软件;ENVI;MODIS;应用引言ENVI软件是一种遥感图像处理软件,是处理、分析并显示多光谱数据、高光谱数据和雷达数据的高级工具,直观且功能强大,获2001年美国权威机构NIMA遥感软件测评第一。
可广泛应用于地质、林业、农业、模式识别、军事、自然资源勘探、水/海岸资源管理、环境和土地利用管理等部门。
EOS-MODIS影像具有36个波段和250~1000m的地面分辨率,利用ENVI软件可对其进行处理和分析。
1遥感图像处理系统ENVI主要功能1.1数据输入输出格式在数据接口方面,ENVI支持众多数据格式,可方便的读入信息源和输出成果,在ENVI中处理分析的专题图可在多种常用RS/GIS软件中显示。
ENVI 支持输入格式:TIFF、JPEG、BMP等通用图像格式;ARC/Info Images(.bil)、ArcViewShape(.shp)、AutoCAD DXF、MapInfo(.mid)等矢量格式;LandsatTM、SPOT、IKONOS、A VHRR、MODIS、Radar、Thermal、Military等遥感数据格式;PCI(.pix)、ERMapper、ERDAS IMAGINE等遥感软件格式;及ASCII、DOQ等数据格式。
ENVI支持的输出格式为ARC/Info Images(.bil)、ArcView Shape(.shp)、ASCII、BMP、ERDAS7.5(.lan)、ERMAPPER、GIF、PCI、RGB 等。
ENVI处理MODIS的步骤全过程〔图解〕制作人:杨鹏〔中国科学院某某光机所〕一:打开MODIS影像数据需要对MODIS数据〔或者产品数据〕进展几何校正,同时需要.shp格式的行政区矢量文件。
然后才能对自己的感兴趣区进展剪裁。
具体步骤如下:1:打开ENVI2:打开所用数据使用ENVI软件打开所用数据〔包括:HDF、TXT、SHP等格式〕,在本例中打开HDF 格式的栅型数据格式。
File—Open External File—Generic Formats—HDF打开数据后如下列图:每个数据集都有不同的数据,可参照MODIS-HDF格式〔可用HDF Explorer查看,里面有详细的说明〕查看自己需要的数据在哪个数据集,然后打开。
依次进展下去可以看到如下对话框:然后选中自己需要的波段,点击Load Band。
就打开了所需的数据图像。
可以看出Cursor Location中显示的是波段里数据的DN值,需要查找偏移量与缩放因子然后进展适当的运算才能求得所需的量〔反射率、发射率等〕。
3:打开图像的另一种方法针对此问题可以直接用另一种形式打开直接可以看到所需量的数据,无需转换。
其过程如下:File--Open Image File以打开31波段的辐射量为例:从图中可以看出,其显示的值直接就是辐射量,而不是DN值,可以直接进展用了,无需转换。
二:对MODIS数据进展几何校正1:打开MODIS影像〔步骤见上〕2:对MODIS数据进展几何校正,步骤如下Basic Tools--Preprocessing--Data-Specific Utilities--MODIS--Georeferance Data选中点击后,进入下一个对话框:在此可以对整个数据集直接批量处理,也可以对单个波段进展处理,随便自己的需要选择,点OK继续,可以看到如下:不需要任何选择,点OK继续。
在此将经纬度格式选好,然后选Memory这样的目的是直接预览无需保存。
MODIS产品下载与数据处理小结遥感基础与处理2010-02-15 23:32:16 阅读2098 评论34 字号:大中小订阅以下是我个人总结的关于MODIS产品下载与数据处理方面的一些内容,供想要入门的朋友参考。
从USGS网站网站可以浏览并下载需要的MODIS数据,不过另一些提供下载的网站还有::22000/WebAccess/drill?attrib=home&next=group (最佳下载渠道)/data/modis/order_data.html (拥有几乎全部方式的下载链接)https:///api/ (WIST 搜索和下载MODIS产品的一个网站)ftp:// (FTP下载方式的MODIS产品网站)其中MODIS产品的介绍可以从网站https:///lpdaac/products/modis_products_table 了解。
无论何种方式下载,都需要等待很长时间,与网速有很大关系,因为需要一景一景的下载,文件大小根据所选数据产品的类型有关。
一、MODIS产品数据的打开数据下载后,在ENVI、ERDAS IMAGINE等遥感图像处理软件中都可以打开。
如果下载的是HDF格式,在ENVI中直接使用[File] ----> [Open Image File]命令打开,选择HDF格式就行。
打开后,产品自动加载入“Available Band List”中,自动识别地理参考信息。
在ERDAS IMAGINE中可以使用[Import]工具选择相应类型打开(这里不作介绍)。
二、MODIS产品数据的感兴趣提取在ENVI中打开后,需要对影像进行感兴趣区域的提取。
这里只介绍一种方法,即使用已有的矢量多边形数据提取影像中的感兴趣区域。
一般这种情况是根据某国家(省、市、地区等)行政边界提取某地区的影像数据。
首先在选中的产品数据的打开窗口(Image Display)中添加矢量边界层([Overlay]--->[Vectors]),然后选择[Basic Tools]中的[Masking]菜单下的[Build Mask]工具,打开后,首先选择你打开影像数据的窗口(Select Input Display),如果你影像在Display #1中打开,就选择选项中的Display #1即可。
E N V I处理M O D I S步骤(图解)ENVI处理MODIS的步骤全过程(图解)制作人:杨鹏(中国科学院安徽光机所)fga12345678@一:打开MODIS影像数据需要对MODIS数据(或者产品数据)进行几何校正,同时需要.shp格式的行政区矢量文件。
然后才能对自己的感兴趣区进行剪裁。
具体步骤如下:1:打开ENVI2:打开所用数据使用ENVI软件打开所用数据(包括:HDF、TXT、SHP等格式),在本例中打开HDF格式的栅型数据格式。
File—Open External File—Generic Formats—HDF打开数据后如图所示:每个数据集都有不同的数据,可参照MODIS-HDF格式(可用HDF Explorer 查看,里面有详细的说明)查看自己需要的数据在哪个数据集,然后打开。
依次进行下去可以看到如下对话框:然后选中自己需要的波段,点击Load Band。
就打开了所需的数据图像。
可以看出Cursor Location中显示的是波段里数据的DN值,需要查找偏移量与缩放因子然后进行适当的运算才能求得所需的量(反射率、发射率等)。
3:打开图像的另一种方法针对此问题可以直接用另一种形式打开直接可以看到所需量的数据,无需转换。
其过程如下: File--Open Image File以打开31波段的辐射量为例:从图中可以看出,其显示的值直接就是辐射量,而不是DN值,可以直接进行用了,无需转换。
二:对MODIS数据进行几何校正1:打开MODIS影像(步骤见上)2:对MODIS数据进行几何校正,步骤如下Basic Tools--Preprocessing--Data-Specific Utilities--MODIS--Georeferance Data选中点击后,进入下一个对话框:在此可以对整个数据集直接批量处理,也可以对单个波段进行处理,随便自己的需要选择,点OK继续,可以看到如下:不需要任何选择,点OK继续。
一般说来,用ENVI打开MODIS HDF数据有以下几种方式:
第一种是直接用File->Open Image File打开,主要是针对Level1B数据和Level2数据的部分波段。
以MOD021KM数据为例,采用这种方式打开得到的图像是定标后的反射率、辐射亮度以及发射率数据,即图像灰度具有明确的物理含义,不需要再进行波段运算进行定标。
这种方法打开数据速度快,但是适用的数据有限,打开后得到的图像波段也有限。
比如MOD02数据中也有经纬度、太阳/传感器天顶角、方位角波段,用这种方式就无法打开。
第二种是是用File->Open External File->Generic Formats->HDF打开,可打开各种产品。
该方法实际上是打开HDF文件,特别是像MODIS的很多陆地产品,如地表反射率、LAI、LST、BRDF/Albedo等(就是文件名中带有h??v??的),都需要用这种方式打开。
打开之后用户还需要选择HDF文件中的数据集(dataset),如果是多波段还需要指定数据格式(BSQ\BIP\BIL)。
采用这种方式打开HDF文件可以获取文件中所有数据集的信息,打开得到的波段也是未做过定标的,需要从HDF文件中查找定标系数通过波段运算手工定标。
查看HDF数据集属性可以通过Basic Tools->Preprocessing->Data-Specific Utilities->View HDF Dataset Attributes实现。
另外通过File->Open External File->EOS->MODIS也可以打开部分MODIS数据,它与第一种打开方式一样,这里不再重复。
关于MODIS数据的几何校正,对于Level1B和Level2级产品,由于其HDF文件中一般都含有经纬度波段,可采用GLT的方法对其进行校正。
相应的菜单是Map->Georeference from Input Geometry->Build GLT和Map->Georeference from Input Geometry->Georeference from GLT。
用GLT方法校正需要注意输入的经纬度图波段要与待校正的数据波段行列大小一致。
在Map菜单下还有一个Georeference MODIS功能,可以对采用Open Image File方式打开得到的MODIS数据波段进行校正。
通常对Level1B数据采用这种方法进行,因为速度快,而且不需要生成GLT临时文件。
但这种方法存在一个问题,就是对很多无法通过Open Image File方式打开的数据波段失效。