zoukankan      html  css  js  c++  java
  • MeteoInfoLab脚本示例:SeaWiFS HDF Grid数据

    SeaWiFS HDF Grid数据读取,特别是涉及到了文件的众多属性数据的读取,数据取对数后绘图。

    脚本程序:

    #Add data file
    f = addfile('D:/Temp/hdf/S1999001.L3m_DAY_CHL_chlor_a_9km.hdf')
    #Get data variable
    vname = 'l3m_data'
    v = f[vname]
    #Set x/y
    ny = f.attrvalue('Number_of_Lines')[0]
    nx = f.attrvalue('Number_of_Columns')[0]
    sx = f.attrvalue('Westernmost_Longitude')[0]
    ex = f.attrvalue('Easternmost_Longitude')[0]
    sy = f.attrvalue('Southernmost_Latitude')[0]
    ey = f.attrvalue('Northernmost_Latitude')[0]
    x = linspace(sx, ex, nx)
    y = linspace(sy, ey, ny)
    #Set x/y dimensions
    v.setdim('Y', y, 0)
    v.setdim('X', x, 1)
    #Get data array
    fillv = v.attrvalue('Fill')[0]
    scale = v.attrvalue('Slope')[0]
    offset = v.attrvalue('Intercept')[0]
    data = v[::-1,:] * scale + offset
    data.fill_value = fillv
    data = log(data)
    #Plot
    axesm()
    world = shaperead('D:/Temp/map/country1.shp')
    geoshow(world)
    levs = arange(-6, 4, 0.2)
    layer = imshowm(data, levs)
    colorbar(layer)
    title('SeaWiFS - ' + vname + ' (log scale)')

  • 相关阅读:
    小工具
    git
    git如何做个人构建
    菜鸟教程
    Xftp和Xshell
    IDEA
    webStorm
    HBuilder
    chrome浏览器
    Vue-Server-Renderer
  • 原文地址:https://www.cnblogs.com/yaqiang/p/4865470.html
Copyright © 2011-2022 走看看