以下是我個(gè)人總結(jié)的關(guān)于MODIS產(chǎn)品下載與數(shù)據(jù)處理方面的一些內(nèi)容，供想要入門的朋友參考。

從USGS網(wǎng)站 http://glovis.網(wǎng)站可以瀏覽并下載需要的MODIS數(shù)據(jù)，不過(guò)另一些提供下載的網(wǎng)站還有：

http://e4eil01.cr.:22000/WebAccess/drill?attrib=home&next=group （最佳下載渠道）

http:///data/modis/order_data.html （擁有幾乎全部方式的下載鏈接）

https://wist.echo./api/ (WIST 搜索和下載MODIS產(chǎn)品的一個(gè)網(wǎng)站）

ftp://e4ftl01u.ecs. （FTP下載方式的MODIS產(chǎn)品網(wǎng)站）

其中MODIS產(chǎn)品的介紹可以從網(wǎng)站 https://lpdaac./lpdaac/products/modis_products_table 了解。

無(wú)論何種方式下載，都需要等待很長(zhǎng)時(shí)間，與網(wǎng)速有很大關(guān)系，因?yàn)樾枰痪耙痪暗南螺d，文件大小根據(jù)所選數(shù)據(jù)產(chǎn)品的類型有關(guān)。

一、MODIS產(chǎn)品數(shù)據(jù)的打開

數(shù)據(jù)下載后，在ENVI、ERDAS IMAGINE等遙感圖像處理軟件中都可以打開。如果下載的是HDF格式，在ENVI中直接使用 [File] ----> [Open Image File]命令打開，選擇HDF格式就行。打開后，產(chǎn)品自動(dòng)加載入“Available Band List”中，自動(dòng)識(shí)別地理參考信息。在ERDAS IMAGINE中可以使用[Import]工具選擇相應(yīng)類型打開（這里不作介紹）。

二、MODIS產(chǎn)品數(shù)據(jù)的感興趣提取

在ENVI中打開后，需要對(duì)影像進(jìn)行感興趣區(qū)域的提取。這里只介紹一種方法，即使用已有的矢量多邊形數(shù)據(jù)提取影像中的感興趣區(qū)域。一般這種情況是根據(jù)某國(guó)家（省、市、地區(qū)等）行政邊界提取某地區(qū)的影像數(shù)據(jù)。首先在選中的產(chǎn)品數(shù)據(jù)的打開窗口（Image Display）中添加矢量邊界層（[Overlay]--->[Vectors]），然后選擇[Basic Tools]中的[Masking]菜單下的[Build Mask]工具，打開后，首先選擇你打開影像數(shù)據(jù)的窗口（Select Input Display），如果你影像在Display #1中打開，就選擇選項(xiàng)中的Display #1即可。點(diǎn)擊[OK]后跳出的對(duì)話框中的[Options]菜單下選擇[Import EVFs]命令（當(dāng)然，這個(gè)根據(jù)你定義感興趣區(qū)域的數(shù)據(jù)選擇，如果是ROI，則選擇[Import ROIs]），選中矢量邊界文件（.evf格式），再確定輸出結(jié)果的方式[File]或[Memory]。一般我選后者。點(diǎn)擊[Apply]后生成MASK文件。利用生成的MASK文件使用[Apply Mask]工具進(jìn)行掩膜操作，即只提取感興趣區(qū)域。

在[Apply Mask Input File]窗口中，選擇原始影像文件，即打開的影像。選中文件后，進(jìn)入[Spatial Subset]中，點(diǎn)擊[ROI/EVF]按鈕，選擇感興趣矢量文件，和剛才定義MASK所用的矢量evf文件相同。[OK]后推出，再點(diǎn)擊[Select Mask Band]按鈕，選擇剛剛生成的MASK文件，完成后，[OK]即可完成感興趣影像的提取。忘記提到的是，MASK掉的部分默認(rèn)值為0.0。

三、MODIS產(chǎn)品數(shù)據(jù)的重新投影

MODIS產(chǎn)品數(shù)據(jù)的投影一般是Sinusoidal，不符合我們正常使用數(shù)據(jù)的查看方式，因此需要進(jìn)行重新投影。一般我選擇UTM，WGS84形式的。以河南省為例，我通過(guò)上一步的感興趣提取操作后，對(duì)感興趣區(qū)域進(jìn)行重新投影。在ENVI中的[Map]菜單下選擇[Convert Map Projection]命令，進(jìn)入對(duì)話框，點(diǎn)擊[Chang Proj...]按鈕，選擇UTM，[Datum]為WGS84，[Zone]為50。[OK]后，再[OK]，完成重新投影。如果打開的影像發(fā)現(xiàn)是全白色的，只有一些區(qū)域是黑色的，那么你將紅色框調(diào)整到某個(gè)位置，進(jìn)行圖像增強(qiáng)就可以了，如線性拉伸、直方圖均衡等。

如果輸出的重投影影像看上去還是感興趣小，圖像大的時(shí)候，你可以再針對(duì)這個(gè)數(shù)據(jù)做一次MASK操作。記得[Spatial Subset]一定要定義好，和上同。

 今天就說(shuō)到這里了，附上河南省的MODIS產(chǎn)品圖（EVI的）。

四、MRT批處理MODIS產(chǎn)品數(shù)據(jù)

 Modis Reprojection Tool 可以批處理Modis產(chǎn)品，如mosaic，resample，reprojection。在windows平臺(tái)下可以在Dos命令框中運(yùn)行自己創(chuàng)建的批處理*.bat文件實(shí)現(xiàn)批處理。值得注意的是：MRT工具讀取文件的路徑不可以包含中文字符，因此，請(qǐng)將文件保存在英文字母的路徑下。這里只以Global 0.05Deg MOD11C3產(chǎn)品為例，提取2001~2007年北京地區(qū)的圖像。使用的命令是resample。

Resample -p *.prm -I *.hdf -o *.tif

*.prm文件如下：

 INPUT_FILENAME = C:\TDDOWNLOAD\MOD11C3\MOD11C3.A2001001.005.2006360120227.hdf

 SPECTRAL_SUBSET = ( 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 )

 SPATIAL_SUBSET_TYPE = INPUT_LAT_LONG

 SPATIAL_SUBSET_UL_CORNER = ( 41.06 115.41 )

SPATIAL_SUBSET_LR_CORNER = ( 39.42 117.5 )

 OUTPUT_FILENAME = C:\TDDOWNLOAD\MOD11C3\output\MOD11C3.A2001001.005.tif

 RESAMPLING_TYPE = BILINEAR

 OUTPUT_PROJECTION_TYPE = GEO

 OUTPUT_PROJECTION_PARAMETERS = (

 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.0 0.0 )

 DATUM = WGS84

 OUTPUT_PIXEL_SIZE = 0.005

批處理命令可在EXCEL中制作完成，如圖

五、使用MCTK進(jìn)行MODIS數(shù)據(jù)的批處理

由于MODIS數(shù)據(jù)量比較大，有時(shí)候批處理使用第四部分的方法還是比較麻煩?？梢允褂肕CTK工具（MODIS CONVERSION TOOLKIT)在IDL編程環(huán)境下進(jìn)行批處理。網(wǎng)站上有MCTK工具的下載鏈接。下面提供本人最近處理一批MOD13Q1數(shù)據(jù)的批處理代碼（僅供參考）：
PRO BATCH_MOD13Q1_PROCESS
    COMPILE_OPT IDL2
    ENVI, /RESTORE_BASE_SAVE_FILES
    ENVI_BATCH_INIT, LOG_FILE='BATCH.LOG'
    PRINT, 'START : ',SYSTIME()
   
    ROOT_DIR = 'G:\MODIS\HENAN_RICE_SEASON\MOD13Q1\'
    FNS = FILE_SEARCH(ROOT_DIR,'*.HDF',COUNT = COUNT)
    PRINT, 'There ara totally', COUNT,' images.'
   
    OUTPUT_LOCATION = 'G:\MODIS\HENAN_RICE_SEASON\MOD13Q1\OUTPUT_LOCATION\'
    GRID_NAME = 'MODIS_Grid_16DAY_250m_500m_VI'
    SD_NAME_EVI = ['250m 16 days EVI']
    SD_NAME_QA = ['250m 16 days VI Quality']
    SD_NAME_NIR = ['250m 16 days NIR reflectance']
    SD_NAME_MIR = ['250m 16 days MIR reflectance']
    ;OUTPUT_METHOD = 1  ;REPROJECTED
    OUTPUT_METHOD = 0 ;Standard
   
    ;投影轉(zhuǎn)換設(shè)定
    ;UNITS = ENVI_TRANSLATE_PROJECTION_UNITS('Meters')
    ;OUTPUT_PROJECTION = ENVI_PROJ_CREATE(/UTM,ZONE=49,UNITS=UNITS)
    ;OUTPUT_PS_X = 250
    ;OUTPUT_PS_Y = 250
   
    ;INTERPOLATION_METHOD = 8 ;TRIANGULATION WITH NEAREST NEIGHBOR
   
    ;導(dǎo)出EVI, MIR, NIR, QA
    FOR i = 0, COUNT-1  DO BEGIN
        FILENAME = FNS[i]
        A = STRPOS(FILENAME,'.')
        OUTPUT_ROOT_NAME_NIR = 'NIR_'+ STRMID(FILENAME,A+1,8)
        CONVERT_MODIS_DATA, IN_FILE = FILENAME, $
            OUT_PATH = OUTPUT_LOCATION, OUT_ROOT=OUTPUT_ROOT_NAME_NIR, $
            /HIGHER_PRODUCT, /GRID, GD_NAME=GRID_NAME,SD_NAME = SD_NAME_NIR, $
            OUT_METHOD = OUTPUT_METHOD, $;OUT_PROJ = OUTPUT_PROJECTION, $
            OUT_PS_X = OUTPUT_PS_X, OUT_PS_Y = OUTPUT_PS_Y, $
            NUM_X_PTS = 50, $
            NUM_Y_PTS=50, $;INTERP_METHOD = INTERPOLATION_METHOD, $
            BACKGROUND='0', FILL_REPLACE_VALUE='0',$
            R_FID_ARRAY=R_FID_ARRAY, R_FNAME_ARRAY=R_FNAME_ARRAY, /NO_MSG
           
        OUTPUT_ROOT_NAME_MIR = 'MIR_' + STRMID(FILENAME,A+1,8)
        CONVERT_MODIS_DATA, IN_FILE = FILENAME, $
            OUT_PATH = OUTPUT_LOCATION, OUT_ROOT=OUTPUT_ROOT_NAME_MIR, $
            /HIGHER_PRODUCT, /GRID, GD_NAME=GRID_NAME,SD_NAME = SD_NAME_MIR, $
            OUT_METHOD = OUTPUT_METHOD, $;OUT_PROJ = OUTPUT_PROJECTION, $
            OUT_PS_X = OUTPUT_PS_X, OUT_PS_Y = OUTPUT_PS_Y, $
            NUM_X_PTS = 50, $
            NUM_Y_PTS=50, $;INTERP_METHOD = INTERPOLATION_METHOD, $
            BACKGROUND='0', FILL_REPLACE_VALUE='0',$
            R_FID_ARRAY=R_FID_ARRAY, R_FNAME_ARRAY=R_FNAME_ARRAY, /NO_MSG
           
        OUTPUT_ROOT_NAME_EVI = 'EVI_' + STRMID(FILENAME,A+1,8)
        CONVERT_MODIS_DATA, IN_FILE = FILENAME, $
            OUT_PATH = OUTPUT_LOCATION, OUT_ROOT=OUTPUT_ROOT_NAME_EVI, $
            /HIGHER_PRODUCT, /GRID, GD_NAME=GRID_NAME,SD_NAME = SD_NAME_EVI, $
            OUT_METHOD = OUTPUT_METHOD, $;OUT_PROJ = OUTPUT_PROJECTION, $
            OUT_PS_X = OUTPUT_PS_X, OUT_PS_Y = OUTPUT_PS_Y, $
            NUM_X_PTS = 50, $
            NUM_Y_PTS=50, $;INTERP_METHOD = INTERPOLATION_METHOD, $
            BACKGROUND='0', FILL_REPLACE_VALUE='0',$
            R_FID_ARRAY=R_FID_ARRAY, R_FNAME_ARRAY=R_FNAME_ARRAY, /NO_MSG
           
        OUTPUT_ROOT_NAME_QA = 'QA_' + STRMID(FILENAME,A+1,8)
        CONVERT_MODIS_DATA, IN_FILE = FILENAME, $
            OUT_PATH = OUTPUT_LOCATION, OUT_ROOT=OUTPUT_ROOT_NAME_QA, $
            /HIGHER_PRODUCT, /GRID, GD_NAME=GRID_NAME,SD_NAME = SD_NAME_QA, $
            OUT_METHOD = OUTPUT_METHOD, $;OUT_PROJ = OUTPUT_PROJECTION, $
            OUT_PS_X = OUTPUT_PS_X, OUT_PS_Y = OUTPUT_PS_Y, $
            NUM_X_PTS = 50, $
            NUM_Y_PTS=50, $;INTERP_METHOD = INTERPOLATION_METHOD, $
            BACKGROUND='0', FILL_REPLACE_VALUE='0',$
            R_FID_ARRAY=R_FID_ARRAY, R_FNAME_ARRAY=R_FNAME_ARRAY, /NO_MSG
    ENDFOR

    PRINT, 'END : ', SYSTIME()
    ENVI_BATCH_EXIT
   
END