遥感数据处理之MCD19A2数据的处理

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

      这一段时间一直在做modis气溶胶产品——MCD19A2这个数据的处理工作，总结了一些东西，在这里分享一下。

先解众惑：为什么使用MCD19A2数据处理出来的结果会有3个或4个或5个波段。

这里使用HDF Exploer这个软件，下载链接：

HDFExploer~~~_MCD19A2-互联网文档类资源-CSDN下载

从下图可以看到打开轨道时间数据集可以看到有4个时间段，末尾为T的表示TERRA星在UTC时间上午2点20和4点拍摄的，末尾为A的表示AQUA星拍摄的，它是将4个时间段拍摄的产品数据相当于做了一个波段合成，一个时间对应一个波段，所以如果拍摄了4次，就会有一个4波段的数据，3次就有一个3波段的数据。

 modis数据一般可以用MCTK和MRT这两个工具进行处理，下面介绍一下使用MCTK处理MCD19A2的流程：

使用MCTK处理。

ps：本次使用的环境为ENVI5.3。

（1）.mctk插件的下载与安装可以参考Modis Conversion Toolkit (MCTK) 工具及用户手册_u014248147的博客-CSDN博客_mctk，也可以在ENVI中的app store进行安装。

（2）.打开mctk软件，选择一个MCD19A2数据，在弹出的第一个对话框内，选择grid1KM，然后根据以下步骤进行处理。

这里附上MCTK用户指导的链接，本来我想直接上传到这里，但是CSDN一直提示我该资源已存在，无法上传，所以这里附上github上的原文链接：

https://github.com/dawhite/MCTK/blob/master/MCTK%20User's%20Guide.pdf

 结果如下，就如上文所说，一个波段对应一个时间拍摄的数据。

                                                                                                                                                                                                        在下方附上IDL里面使用mctk对MCD19A2数据进行批处理代码

pro mcd19_process
  compile_opt idl2
  e=envi(/headless)
  file='E:\MCD19A2\'  ;file为存储MCD19A2数据的目录
  ds=file_search(file,'*.hdf',count=n)

  for i=0,n-1 do begin
    modis_file =ds[i]
    output_location= 'F:\MCD19A2_1\'  ;输出路径
    base_name=file_basename(ds[i])
    output_rootname1= base_name
    name_split=STRSPLIT(base_name, '.', /EXTRACT)
    mod_version=name_split[0]

    ;  grid_name1 = 'MOD_Grid_MOD17A2H'
    sd_names1 = ['Optical_Depth_055']   ;我这里是对550nm处的AOD进行处理
    grid_name1 = 'Grid1km'    ;打开的是Grid1KM的数据集
    out_method= 1            ;对应插件中选择投影输出类型
    output_projection = envi_proj_create(/geographic)
    out_ps1= [0.01015055,0.00835245]  ;定义输出像元，可以从ENVI中View metadata中查看
    interpolation_method = 2   ;对应插值的方法
    ;Set reprojection background and any native fill values to NaN
    nan_fill= float('NaN')    ;无效值填充为NaN
    bridges = mctk_create_bridges()
    convert_modis_data,in_file=modis_file,out_path=output_location, out_root=output_rootname1, $
      gd_name=grid_name1, sd_names=sd_names1,out_method=out_method, out_proj=output_projection, $
      out_ps=out_ps1, interp_method=interpolation_method,background=nan_fill, fill_replace_value=nan_fill, $
      r_fid_array=r_fid_array, r_fname_array=r_fname_array,bridges=bridges,msg=msg
    mctk_destroy_bridges, bridges

  endfor

end

 如果要对MCD19A2数据进行后续处理（多个波段进行有效值均值合成为单波段数据），可以跳转到下篇这篇博客。

遥感数据处理之MCD19A2数据的后续处理_三十二号星期八的博客-CSDN博客

如有错误，欢迎指正，万分感谢！