MCD19A2大气遥感数据处理

Question

有没有人会处理MCD19A2数据的啊？批量处理，将时间、aod、aodqa、经纬度读取出来。本人刚研一，最近也是刚刚了解数据结构，python基础几乎为零，从网上找代码学习也只是稀里糊涂的。时间是读取Oribt_time_stamp里面的。

Answer 1

题主你好,我新写了篇文章看不知道能否帮到你，是用PYTHON 处理 MCD19A2(C6),文章链接：http://t.csdn.cn/ZX35B
目前已经出了C6.1产品精度应该会有较大的提升（https://lpdaac.usgs.gov/products/mcd19a2v061/）

不知道您是否使用GEE，GEE可以不用下载繁多的数据直接在服务器云端上处理速度比自己批量处理要快的多，以下为GEE代码示例：

//所选区域范围 
var Bounds = 
    ee.Geometry.Polygon(
         [[[53.39306671552197, 54.598371442447],
          [53.39306671552197, 0.7485052185944432],
          [163.08056671552197, 0.7485052185944432],
          [163.08056671552197, 54.598371442447]]], null, false);
var MAIAC_ImageCollection = ee.ImageCollection('MODIS/006/MCD19A2_GRANULES')//选择数据集
      .filterDate("2022-01-01", "2022-01-02")//日期筛选
      .filter(ee.Filter.eq('SATELLITE','A'))   //T代表Terra A 代表 Aqua
      .filterBounds(Bounds);//区域范围筛选

var QA_MAIAC_ImageCollection = MAIAC_ImageCollection.map(function(img) {
  //质量筛选
  //bitmask质量筛选函数
  function bitwiseExtract(value, fromBit, toBit) {
    if (toBit === undefined) toBit = fromBit
    var maskSize = ee.Number(1).add(toBit).subtract(fromBit)
    var mask = ee.Number(1).leftShift(maskSize).subtract(1)
    return value.rightShift(fromBit).bitwiseAnd(mask)
  }
  //QA筛选
  var qa = img.select("AOD_QA");//选择QA波段
  var aod550=img.select('Optical_Depth_055').multiply(0.001)//乘以比例因子
  var qamask=bitwiseExtract(qa,8,11).eq(0)//8-11位等于0的也就是 Bits 8-11: QA for AOD  0: Best quality 最好质量的 GEE可以在数据集里查找
  var Best_Aod550=aod550.updateMask(qamask).rename("Best_Aod550")//掩膜提取
  //addband 增加波段
  img=img.addBands(Best_Aod550)
  return img
})
print(QA_MAIAC_ImageCollection)
Map.addLayer(QA_MAIAC_ImageCollection.select('Best_Aod550').mean(), {min: 0,max: 1,palette: ['black', 'blue', 'purple', 'cyan', 'green', 'yellow', 'red']}, 'Best_Aod550');//添加地图显示