2020-07-18 数据预处理

数据转换 HDF --> TIFF

• MOD09GA、MCD43A1
  • 工具：MRT
  • 拼接：H25V4, H26V4, H26V5, H27V5四个区域
  • 投影转换：geographic
  • 分辨率：0.01 decimal degrees
• MOD11_L2
  • 工具：MRTSwath IDL
  • ditto
  • 对每个hdf文件，运行IDL转成dat，然后用ENVI打开dat转成tiff

研究区域的确定

• 计算出研究区域对应的经纬度BBOX
  • minLon = 105.78432759485935
  • maxLon = 111.88621521850175
  • minLat = 36.95825861049653
  • maxLat = 41.45847714860637
• 对每一组数据
  • L2 Swath数据转换后的图像：根据经纬度BBOX计算对应图像中的BBOX坐标
  • L3 网格化的数据转换后的图像：直接使用0812数据计算的BBOX坐标
• 图像中任意（x, y）点的经纬度计算

  • from osgeo import gdal
    ds = gdal.Open('path/to/file')
    width = ds.RasterXSize
    height = ds.RasterYSize
    gt = ds.GetGeoTransform()
    # bbox
    minLon = gt[0]
    minLat = gt[3] + width*gt[4] + height*gt[5] 
    maxLon = gt[0] + width*gt[1] + height*gt[2]
    maxLat = gt[3]
    # any point (xy以左上角为原点)
    lat = gt[3] + x*gt[4] + y*gt[5] 
    lon = gt[0] + x*gt[1] + y*gt[2]
    

• 由经纬度坐标计算图像坐标，即以上公式的反过程