摘要：MeteoInfoLab的georead函数提供了读取shape文件、image文件（JPG、PNG等，需要有相应的地理定位文件）文件生成图层的功能（事实上shaperead也是同样的功能，不过函数名太过具体）。对于Polygon图层可以设置其颜色透明度，可以用makecolors函数生成颜色列表，... 阅读全文

摘要：示例数据：ftp://ftp.bom.gov.au/anon/sample/catalogue/Satellite/IDE00220.201507140300.nc该数据的分辨率很高（22000*22000），如果全部读入一个数组需要很大的内存，可以在读入时设置step（比如4，以间隔4读取数据），... 阅读全文

摘要：MODIS卫星很多陆面数据都是Sinusoidal投影，数据被分为一个个10*10度（赤道地区）的瓦片（http://modis-land.gsfc.nasa.gov/MODLAND_grid.html），数据格式通常是HDF-EOS格式，作为EOS格式MeteoInfo可以自动读出数据的坐标信息（... 阅读全文

摘要：利用比湿和温度计算相对湿度的函数是qair2rh(qair, temp, press=1013.25)，三个参数分别是比湿、温度和气压，气压有一个缺省值1013.25，因此计算地面相对湿度的时候也可以不给气压参数。下面的例子打开CFS模式的地面结果，并从比湿和温度数组计算出相对湿度数组并绘图。脚本程... 阅读全文

摘要：在MeteoInfo中添加了Google网络地图功能。打开MeteoInfo软件，选中图层管理区的一个Map Frame（New Map Frame），点击鼠标右键，在弹出菜单中点击Add Web Layer菜单添加一个网络图层，缺省是OpenStreetMap。双击网络图层名，在图层属性对话框中选... 阅读全文

摘要：需要温度和风场U/V分量格点数据，计算中主要用到cdiff函数，结果用GrADS验证一致。脚本程序：print 'Open data files...'f_air = addfile('D:/Temp/nc/air.2011.nc')f_uwnd = addfile('D:/Temp/nc/uwnd... 阅读全文

摘要：这里用美国做例子，有一个美国区域的格点温度场数据（usgrid.data），需要计算出每个州（state）的平均温度。当然需要有一个包含各州行政区域的shape文件了（相关文件可以在此帖中下载：http://bbs.06climate.com/forum.p ... d=11070&fromuid=... 阅读全文

摘要：在脚本中使用中文需要指明是unicode编码，即在含有中文的字符串前加u，比如：u'中文'。还需要将字体指定为一种中文字体。详见下面的例子。脚本程序：x = [1,2,3,4]y = [1,4,9,16]plot(x, y, 'r--*')ylabel(u'Y轴', fontname=u'宋体', ... 阅读全文

摘要：OMI卫星格点数据的例子，全球臭氧柱总量分布。脚本程序：#Add data filefolder = 'D:/Temp/hdf/'fns = 'OMI-Aura_L3-OMTO3e_2005m1214_v002-2006m0929t143855.he5'fn = folder + fnsf = ad... 阅读全文

摘要：这个例子读取OMI卫星Swath数据中的CloudFaction变量并绘图。脚本程序：#Add data filefolder = 'D:/Temp/hdf/'fns = 'OMI-Aura_L2-OMNO2_2008m0720t2016-o21357_v003-2008m0721t101450.h... 阅读全文

摘要：TOMS (Total Ozone Mapping Spectrometer)数据是全球臭氧观测。脚本程序：#Add data filefolder = 'D:/Temp/hdf/'fns = 'TOMS-EP_L3-TOMSEPL3_2000m0101_v8.HDF'fn = folder + f... 阅读全文

摘要：利用Unidata netCDF Java库对远程文件的读取能力（OpenDAP, ADDE, THREDDS等），可以读取远程文件并绘图。脚本程序：fn = 'http://monsoondata.org:9090/dods/model'f = addfile(fn)v = f['ps']ps =... 阅读全文

摘要：MERRA是NOAA的一种再分析资料，HDF数据遵循COARDS协议，读取比较简单。脚本程序：#Add data filefolder = 'D:/Temp/hdf/'fns = 'MERRA300.prod.assim.tavg1_2d_slv_Nx.20150101.SUB.hdf'fn = f... 阅读全文

摘要：TRMM 3B43是卫星观测月平均降水量产品，是HDF的格点数据。需要注意的是数据中降水变量维的顺序里经度维在前纬度维在后，这与通常的设置（纬度维在前经度维在后）相反，需要对获取的二维数组进行转置，使用transpose函数。数据的经纬度数组也需要定义（文件中没有）。数据中的0值为缺测值（fill_... 阅读全文

摘要：用U/V分量数据计算涡度和散度，计算涡度的函数是hcurl，计算散度的函数是hdivg，参数都是U, V。脚本程序：f = addfile('D:/Temp/GrADS/model.ctl')u = f['U'][0,0,:,:]v = f['V'][0,0,:,:]vort = hcurl(u, ... 阅读全文

摘要：AIRS Grid HDF数据是HDF4 EOS格式，数据地理坐标信息可以被MeteoInfo自动识别，脚本程序更为简单。需要注意的是读取数据时Y轴是反向的（卫星数据通常如此）。脚本程序：#Add data filefolder = 'D:/Temp/hdf/'fns = 'AIRS.2002.08... 阅读全文

摘要：脚本程序：#Add data filefn = 'D:/Temp/hdf/FY2C_CLC_MLT_NOM_20070730_1800.hdf'f = addfile(fn)#Get data variablev = f['FY2C_Hourly_Cloud_Classification']#Set... 阅读全文

摘要：全球长时间序列中尺度气旋数据（http://cioss.coas.oregonstate.edu/eddies/）有netCDF格式，散点数据类型，只有一个很大的维Nobs = 2590938。尝试读取了部分经纬度和track数据，并绘制散度图。脚本程序：fn = 'D:/Temp/nc/track... 阅读全文

摘要：例子中的AIRS Swath HDF数据在Polar Stereographic（南极）投影中接近矩形，需要先从数据中读出经纬度及相关数据数组，利用surfacem函数绘制Swath数据（散点），在surfacem函数中经纬度数据会被投影到目的投影坐标系（在axesm函数中定义），并从2维散点数据生... 阅读全文

摘要：TRMM 2A12 HDF数据是卫星观测的SWATH数据（轨道数据），比格点数据处理起来要麻烦一些。数据的经纬度保存在geolocation变量中，需要先将经纬度数据读出来（均为2维数组），然后读取云水含量数据（cldWater）。虽然都是2维数组，但并不是格点数据，相当于2维的散点数据，点数相当多... 阅读全文

摘要：FY3A卫星有HDF格式的AOD产品数据，全球范围，分辨率为0.05度。读取数据文件变量后要重新设定X/Y维，数据是Y轴反向的，且需要除以10000得到AOD值。脚本程序：#Add data filefn = 'D:/Temp/hdf/FY3A_MERSI_GBAL_L3_ASL_MLT_GLL_2... 阅读全文

摘要：MeteoInfoLab中读取文本文件数据的函数是asciiread，获取文本文件行、列数的函数是numasciirow和numasciicol，和NCL中函数名一致，但都是小写字母。本例中的示例数据文件来自此帖中的数据flood.dat：http://bbs.06climate.com/forum... 阅读全文