Windows下MRT批處理MODIS數(shù)據(jù)

1、Windows下MRT批處理MODIS數(shù)據(jù)1.MRT安裝   LP DAAC最新發(fā)布的MRT支持windows平臺下處理MODIS影像，能夠完全不需要通過cygwin來模擬LINUX系統(tǒng)。MRT下載地址如下：在成功注冊之后，就可以點擊下載Windows NT+32-bit     解壓之后得到MRT_Win.zip、mrt_install.bat、reg_set.exe、unzip.exe四個文件，運行mrt_install.bat文件，然后按照提示輸入MRT和JAVA安裝路徑就可以正確安裝并運行MRT。ps:安裝前，需要確定電腦是否安裝了ja

2、va，若沒有需要去下載最新版的java，然后記住安裝路徑，以便之后的輸入。  2.數(shù)據(jù)準(zhǔn)備   本人的MRT安裝在 c:modistools目錄下，運行c:modistoolsbinModistool.bat 進入MRT GUI界面，選擇一副影像制作批量處理需要的*.prm文件。     影像中參數(shù)設(shè)置如圖，需要注意的是，在設(shè)置輸出影像時需要確定輸出影像的格式如E:mA2000049.tif ， 最重要的是要點saveparameter file 保存A200004

3、9.prm文件，保存后不需要run，直接退出MRT GUI即可。     將保存的A2000049.prm文件放到需要處理的MODIS *.hdf格式的影像數(shù)據(jù)的文件目錄中，如E:m中。  3.cmd命令實現(xiàn)MRT批處理    運行cmd 命令，將工作目錄設(shè)置到 c:modistoolsbin 中，即MRT安裝目錄中的 bin 文件夾中    輸入 java -jar MRTBatch.jar -d e:m -p  

4、e:mA2000049.prm -o e:m 其中，-d 表示的是影像數(shù)據(jù)存儲的目錄，-p 表示經(jīng)過MRT GUI處理的prm文件路徑，-o 表示輸出路徑。這串命令表示的是對所有的影像數(shù)據(jù)批處理得到每個影像的拼接和重采樣的 prm 文件。    運行成功并得到所有影像的 prm 文件后，繼續(xù)在輸入 MRTBatch.bat 命令，執(zhí)行這個bat文件，即進行影像的批處理。 上圖中是顯示的是敲擊的命令，以及運行得到的結(jié)果，兩個MOD13Q1.的prm文件就是第

5、一個語句命令生成，下面的tif影像即為重新投影得到的結(jié)果。因為本人所用的數(shù)據(jù)是MOD13Q1 包括12個波段的數(shù)據(jù)，在prm文件中的 -s '1 1  1 1 1 1 1 1 1 1 1 1' 語句表示要輸出的波段，這個在MRT GUI界面中左邊的中間部分可以選擇自己想要的波段。    在輸出影像的其他波段的時候，應(yīng)該將利用MRT GUI生成的prm文件中SPECTRAL_SUBSET = ( 1 )#ORIG_SPECTRAL_SUBSET = ( 0 0 0 0 0 0 0 0 0 0 0 1 )

6、0;修改為SPECTRAL_SUBSET = ( 0 0 0 0 0 0 0 0 0 0 0 1 )#ORIG_SPECTRAL_SUBSET = ( 0 0 0 0 0 0 0 0 0 0 0 1 )。其中，1表示想輸出的波段。     ps：如果在輸入MRTBatch.bat時出現(xiàn) error，說沒找到環(huán)境變量、MRTDATADIR沒有被定義的錯誤?？梢栽谧员碇蠬KEY_CURRENT_USEREnvironment中設(shè)置環(huán)境變量，設(shè)置如下  MRT_HOME c:modistools  MRT_DATA_D

7、IR c:modistoolsdata   Pathc:modistoolsbin   設(shè)置完成后，建議重啟電腦之后再進行操作。   如果出現(xiàn) java不是內(nèi)部或外部命令的時候，建議關(guān)掉cmd，然后重新打開即可解決問題。 若還是不能解決問題，建議重新安裝java。  4.基于MRT的Matlab實現(xiàn)MODIS數(shù)據(jù)重采樣   matlab純屬個人練手的結(jié)果，其中不足的地方希望得到高手的指教。本代碼能夠借助MRT對影像進行重采樣，不過首先還是需要利用MRT GUI生成一個prm文

8、件，作為批量操作的樣本。影像數(shù)據(jù)存儲在e:mdata , prm文件生成的路徑為 e:mprm 。這個程序的*.m文件以及MATlLAB的工作路徑需要放在MRT安裝目錄的bin文件夾中才能夠順利的運行。關(guān)于這一點，本人希望能夠得到改進，但是目前還不知道怎樣解決才比較合適，求高手指導(dǎo)。 代碼如下:   clear all；clc;workingpath='e:m'x=workingpath,'data' mylist=dir(x);filenum=length(mylist);batrunfile=fop

9、en('c:modistoolsbinrunbfi.bat','wt'); for i=3:filenum    inputfilename=mylist(i).name(11:23);     sourcename=workingpath,'A2000049.prm'%   copyname=workingpath,'prm',inputfilename,'_resample','.prm'   

10、copyfile(sourcename,copyname);        copyname1=fopen(copyname,'wt');        myinputxt1='INPUT_FILENAME = ',x,'',mylist(i).name;    myinputxt2='SPECTRAL_SUBSET = ( 1 0 0 0 0 0 0 0 0 00 0 )'    myinp

11、utxt3='SPATIAL_SUBSET_TYPE = INPUT_LAT_LONG'    myinputxt4='SPATIAL_SUBSET_UL_CORNER = (39.999999996 95.7759995131 )'    myinputxt5='SPATIAL_SUBSET_LR_CORNER = (19.999999998 130.540728914 )'    myinputxt6='OUTPUT_FILENAME =',x,'',

12、inputfilename,'.tif'    myinputxt7='RESAMPLING_TYPE = NEAREST_NEIGHBOR'    myinputxt8='OUTPUT_PROJECTION_TYPE = GEO'    myinputxt9='OUTPUT_PROJECTION_PARAMETERS = ('    myinputxt10=' 0.0 0.0 0.0'    myinputxt11

13、=' 0.0 0.0 0.0'    myinputxt12=' 0.0 0.0 0.0'    myinputxt13=' 0.0 0.0 0.0'    myinputxt14=' 0.0 0.0 0.0 )'    myinputxt15='DATUM = WGS84'        fprintf(copyname1,'n%sn',myinputxt1);

14、0;   fprintf(copyname1,'n%sn',myinputxt2);    fprintf(copyname1,'n%sn',myinputxt3);    fprintf(copyname1,'n%sn',myinputxt4);    fprintf(copyname1,'%sn',myinputxt5);    fprintf(copyname1,'n%sn',myinputxt6); 

15、60; fprintf(copyname1,'n%sn',myinputxt7);    fprintf(copyname1,'n%sn',myinputxt8);    fprintf(copyname1,'n%sn',myinputxt9);    fprintf(copyname1,'%sn',myinputxt10);    fprintf(copyname1,'%sn',myinputxt11);    fprintf(copyname1,'%sn',myinputxt12);    fprintf(copyname1,'%sn',myinputxt13);    fprintf(copyname1,'%snn',myinputxt14);    fprintf(copyname1,'%sn',myinputxt15);        y= 'resamp