zoukankan      html  css  js  c++  java
  • MeteoInfoLab脚本示例:FY-3A AOD HDF数据

    FY3A卫星有HDF格式的AOD产品数据,全球范围,分辨率为0.05度。读取数据文件变量后要重新设定X/Y维,数据是Y轴反向的,且需要除以10000得到AOD值。

    脚本程序:

    #Add data file
    fn = 'D:/Temp/hdf/FY3A_MERSI_GBAL_L3_ASL_MLT_GLL_20140331_AOAM_5000M_MS.HDF'
    f = addfile(fn)
    #Get data variable
    v = f['Aerosol_Optical_Thickness_of_MERSI_550nm']
    #Set x/y
    x = arange1(-30.0,7200, 0.05)
    y = arange1(-90.0,3600, 0.05)
    #Set x/y dimensions
    v.setdim('Y', y, 0)
    v.setdim('X', x, 1)
    #Get data array
    aod = v[:,:]
    #aod = v[[15, 55],[70, 140]]
    aod = aod[::-1,:] / 10000.
    aod = aod[[15, 55],[70, 140]]
    #Plot
    axesm()
    world = shaperead('D:/Temp/map/country1.shp')
    china = shaperead('D:/Temp/map/bou2_4p.shp')
    geoshow(china, edgecolor='gray')
    geoshow(world)
    levs = arange1(0, 20, 0.05)
    cols = makecolors(21, cmap='matlab_jet')
    cols[0] = 'w'
    layer = imshowm(aod, levs, colors=cols)
    colorbar(layer)
    title('FY-3A AOD (2014-03-31)')

  • 相关阅读:
    docker镜像
    docker常用命令
    docker基础
    跨站脚本漏洞(XSS)基础
    Session、Cookie与Token
    linux之curl工具
    ssl证书与java keytool工具
    mysql主从复制
    linux之平均负载(学习笔记非原创)
    mysql8.0忘记密码如何操作?
  • 原文地址:https://www.cnblogs.com/yaqiang/p/4621153.html
Copyright © 2011-2022 走看看