利用Envi Idl和MRT 对MODIS影像批量转投影

119 阅读 0 评论 79 点赞

我是靠谱客的博主矮小花瓣，这篇文章主要介绍利用Envi Idl和MRT 对MODIS影像批量转投影，现在分享给大家，希望可以做个参考。

1、Envi IDL

复制代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
pro enviprojection
  COMPILE_OPT IDL2
  envi, /restore_base_save_files
  envi_batch_init, log_file='batch.txt'
  cd,'F:MCD12Q2'
  HDFFiles = FILE_Search("*.HDF")
  FileCount = N_ELEMENTS(HDFFiles)
  IF FileCount EQ 0 THEN RETURN
  o_proj = ENVI_PROJ_CREATE(/geographic)
  FOR NX =0,FileCount -1 DO BEGIN
    FileName = HDFFiles[NX]
    envi_open_file, FileName , r_fid=fid
    if (fid eq -1) then begin
      envi_batch_exit
      return
    endif
    envi_file_query, fid[0], dims=dims, nb=nb
    pos  = lindgen(nb)
    indexstr = strpos( STRUPCASE(FileName),".HDF")
    out_name = strmid(FileName,0,indexstr) + "_Geo.raw"
    o_pixel_size = [1000,1000]  ;
    envi_convert_file_map_projection, fid=fid, $
      pos=pos, dims=dims, o_proj=o_proj, $
      o_pixel_size=o_pixel_size, grid=[10,10], $
      out_name=