MeteoInfoLab脚本示例：读取文本文件绘制散点图

MeteoInfo

MeteoInfoLab中读取文本文件数据的函数是asciiread，获取文本文件行、列数的函数是numasciirow和numasciicol，和NCL中函数名一致，但都是小写字母。本例中的示例数据文件来自此帖中的数据flood.dat：http://bbs.06climate.com/forum.php?mod=viewthread&tid=29257&extra=page%3D1，该数据文件没有文件头（直接是数据），共三列，分别是经度、纬度和值。先获取文件的行数和列数，然后用asciiread读取文件数据到一个多维数组，数据的维设定需要给定：shape=(nrow,ncol)。然后从数组中取出经度、纬度和值的一维数组（lon, lat, v），用scatterm函数绘制散点图。

脚本程序：

1. fn = 'D:/Temp/ascii/flood.dat'
   
2. ncol = numasciicol(fn)
   
3. nrow = numasciirow(fn)
   
4. a = asciiread(fn,shape=(nrow,ncol))
   
5. lon = a[:,0]
   
6. lat = a[:,1]
   
7. v = a[:,2]
   
8. axesm()
   
9. mlayer = shaperead('D:/Temp/map/bou2_4p.shp')
   
10. geoshow(mlayer, facecolor=(220,220,220), edgecolor=None)
    
11. layer = scatterm(lon, lat, v, edge=False, size=2)
    
12. xlim(109.2, 111)
    
13. ylim(19.8, 21.8)
    
14. title('Flood plain area')
    
15. ylabel('Latitude')
    
16. xlabel('Longitude')
    
17. colorbar(layer)

fn = 'D:/Temp/ascii/flood.dat'<br /> ncol = numasciicol(fn)<br /> nrow = numasciirow(fn)<br /> a = asciiread(fn,shape=(nrow,ncol))<br /> lon = a[:,0]<br /> lat = a[:,1]<br /> v = a[:,2]<br /> axesm()<br /> mlayer = shaperead('D:/Temp/map/bou2_4p.shp')<br /> geoshow(mlayer, facecolor=(220,220,220), edgecolor=None)<br /> layer = scatterm(lon, lat, v, edge=False, size=2)<br /> xlim(109.2, 111)<br /> ylim(19.8, 21.8) <br /> title('Flood plain area')<br /> ylabel('Latitude')<br /> xlabel('Longitude')<br /> colorbar(layer)

MeteoInfo

半颗麦芽糖 发表于 2016-12-5 23:50
王老师，如何横纵坐标是别的量该如何处理呢？我照这个画出来横纵坐标都默认表示为经纬度（加了°E），如果 ...

把脚本和截图贴出来看看

oubiwang

Mark,感觉很快就会用到了

chary

王老师，
您给的示例脚本中shape是不是要事先定义一下？

chary

不好意思 没问题 是我搞错了

jzd19851102

王老师，如果数据是多维数组 能通过shape设定吗
比如是nrow*ncol*m*n    m是时次  n是变量数

MeteoInfo

jzd19851102 发表于 2016-6-8 09:26
王老师，如果数据是多维数组 能通过shape设定吗
比如是nrow*ncol*m*n    m是时次  n是变量数

可以

chennizi

MeteoInfo 发表于 2016-6-8 12:06
可以

王老师，我发现这些脚本绘图的图形，都是一张一张，要点击画图区域右键选择保存的，1、请问有没有一条语句可以在画图结束后直接将图形保存为某种格式存在某个地方？2、在meteoinfolab里面输入脚本后，我都是点击工具栏里面的黑色三角形来运行程序画图的，能不能直接在console里面输入某个命令来运行，我试过输入run script，不可以运行？3、感觉这些脚本都是一次只能出一张图啊，能不能在脚本里面加入循环语句，让一个脚本读许多类似文件出几张不同的图，同时把图一张一张存在一个地方？

MeteoInfo

chennizi 发表于 2016-6-13 10:12
王老师，我发现这些脚本绘图的图形，都是一张一张，要点击画图区域右键选择保存的，1、请问有没有一条语 ...

参考这里：http://www.meteothinker.com/docs ... s/plot/savefig.html

半颗麦芽糖

王老师，如何横纵坐标是别的量该如何处理呢？我照这个画出来横纵坐标都默认表示为经纬度（加了°E），如果横坐标是时间，纵坐标是高度该用什么语句？