2015年7月6日

MeteoInfoLab脚本示例:AIRS Swath HDF数据

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

posted @ 2015-07-06 21:22 yaqiang 阅读(642) 评论(0) 推荐(0)

2015年7月5日

MeteoInfoLab脚本示例:TRMM 2A12 HDF数据

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

posted @ 2015-07-05 19:12 yaqiang 阅读(538) 评论(0) 推荐(0)

2015年7月4日

MeteoInfoLab脚本示例:FY-3A AOD HDF数据

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

posted @ 2015-07-04 19:37 yaqiang 阅读(711) 评论(3) 推荐(0)

2015年7月3日

MeteoInfoLab脚本示例:读取文本文件绘制散度图

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

posted @ 2015-07-03 15:38 yaqiang 阅读(1341) 评论(0) 推荐(0)

2015年6月30日

MeteoInfoLab脚本示例:站点数据绘制等值线

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

posted @ 2015-06-30 11:12 yaqiang 阅读(3851) 评论(0) 推荐(0)

MeteoInfoLab脚本示例:添加南海脚图

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

posted @ 2015-06-30 11:10 yaqiang 阅读(774) 评论(0) 推荐(0)

MeteoInfoLab脚本示例:地图投影

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

posted @ 2015-06-30 11:08 yaqiang 阅读(1255) 评论(0) 推荐(0)

MeteoInfoLab脚本示例:加载地图图层

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

posted @ 2015-06-30 11:06 yaqiang 阅读(2967) 评论(0) 推荐(0)

2015年6月26日

MeteoInfoLab脚本示例:站点数据散点图

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

posted @ 2015-06-26 14:05 yaqiang 阅读(1358) 评论(2) 推荐(0)

MeteoInfoLab脚本示例:inpolygon

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

posted @ 2015-06-26 11:52 yaqiang 阅读(1235) 评论(0) 推荐(0)

导航