MeteoInfo

摘要： 例子中的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... 阅读全文

摘要： 站点数据绘制等值线需要首先将站点数据插值为格点数据，MeteoInfo中提供了反距离权法（IDW）和cressman两个方法，其中IDW方法可以有插值半径的选项。这里示例读取一个MICAPS第一类数据（地面全要素观测），获取6小时累积降水数据（Precipitation6h），然后用站点数据的gri... 阅读全文

摘要： 添加南海脚图实际上就是增加一个地图坐标系（Axesm），在用axesm函数创建地图坐标系的时候可以指定position参数（即位置和大小），第1和第2个参数是坐标系的左下角x, y位置，后两个参数是坐标系的宽度和高度，单位都是Figure的归一化坐标（0 - 1）。由于地图坐标系在Figure缩放时... 阅读全文

摘要： 在用axesm函数创建地图坐标系的时候可以指定地图投影（设置projinfo参数），地图投影可以通过projinfo函数来创建，里面的参数依据proj4投影字符串，可以参考此网页：http://remotesensing.org/geotiff/proj_list。投影后的地图坐标系显示范围的设置可... 阅读全文

摘要： 应用最广泛的的地图数据应该是shape格式，网络上有很多免费下载资源。MeteoInfoLab中读取shape文件的函数是shaperead，参数即文件名，返回数据包含图形和属性信息的图层对象。矢量图层按照图元类型通常分点、线、面三种，都可以通过geoshow函数加入到地图坐标系（Axesm）中。在... 阅读全文

摘要： 这里演示从micaps第一类数据（地面全要素观测）中读取一个变量（用DimDataFile类的stationdata方法），然后maskout掉中国区域之外的数据，利用scatterm函数绘制散点图。脚本程序：f = addfile_micaps('D:/Temp/micaps/10101414.0... 阅读全文

摘要： inpollygon函数是用来判断带坐标（x/y）的数据是否在某个或者一组多边形（Polygon）中，返回的结果中如果做多边形内则值为1，否则值为-1。下面一个例子演示了利用一个shape文件和inpolygon函数生成这种0、1数据。需要下载最新的MeteoInfo版本（1.2.6R1）：http... 阅读全文