不同时间和高度风速的变化图画图出错

半颗麦芽糖

王老师，现在我有一个txt格式的数据，第一列是时次，第二列是高度，第三列是风速，想要画一张不同时间和高度风速的变化图，自己编了一个脚本，但是画图的时候内容一直都是纯蓝色的，我考虑到风速之间间隔太小的问题也相应调小了色标，但是还是一样，附上脚本、数据和图，请王老师指导，多谢！ 11.txt (25 KB, 下载次数: 30)

2016-12-24 18:52 上传

点击文件名下载附件
fn = 'F:/11.txt'
ncol = numasciicol(fn)
nrow = numasciirow(fn)
a = asciiread(fn,shape=(nrow,ncol))
time = a[:,0]
height = a[:,1]
aod = a[:,2]
axesm()
#levs = arange(0.200,0.700,0.010)
layer = contourf(time, height, aod)
colorbar(layer)

Y.L.Zhang

来顶顶顶……

MeteoInfo

等值线分析需要二维数组。

1. fn = 'C:/Temp/11.txt'
   
2. ncol = numasciicol(fn)
   
3. nrow = numasciirow(fn)
   
4. a = asciiread(fn,shape=(nrow,ncol))
   
5. time = a[:,0]
   
6. height = a[:,1]
   
7. aod = a[:,2]
   
8. tn = 144
   
9. hn = len(aod) / tn
   
10. tt = time[:144]
    
11. h = height[::144]
    
12. aod = aod.reshape((hn, tn))
    
13. layer = contourf(tt, h, aod)
    
14. colorbar(layer)

fn = 'C:/Temp/11.txt'<br /> ncol = numasciicol(fn)<br /> nrow = numasciirow(fn)<br /> a = asciiread(fn,shape=(nrow,ncol))<br /> time = a[:,0]<br /> height = a[:,1]<br /> aod = a[:,2]<br /> tn = 144<br /> hn = len(aod) / tn<br /> tt = time[:144]<br /> h = height[::144]<br /> aod = aod.reshape((hn, tn))<br /> layer = contourf(tt, h, aod)<br /> colorbar(layer)

readtable也可以：

1. fn = 'C:/Temp/11.txt'
   
2. table = readtable(fn, format='%i%i%f', headerlines=-1, \
   
3.     colnames=['Time','Height','WindSpeed'])
   
4. tt = table['Time']
   
5. h = table['Height']
   
6. ws = table['WindSpeed']
   
7. tn = 144
   
8. hn = len(ws) / tn
   
9. tt = tt[:144]
   
10. h = h[::144]
    
11. ws = ws.reshape((hn, tn))
    
12. layer = contourf(tt, h, ws)
    
13. colorbar(layer)

fn = 'C:/Temp/11.txt'<br /> table = readtable(fn, format='%i%i%f', headerlines=-1, \<br />     colnames=['Time','Height','WindSpeed'])<br /> tt = table['Time']<br /> h = table['Height']<br /> ws = table['WindSpeed']<br /> tn = 144<br /> hn = len(ws) / tn<br /> tt = tt[:144]<br /> h = h[::144]<br /> ws = ws.reshape((hn, tn))<br /> layer = contourf(tt, h, ws)<br /> colorbar(layer)

半颗麦芽糖

MeteoInfo 发表于 2016-12-24 23:12
等值线分析需要二维数组。

王老师太给力了，32个赞

bbs

学习了~~~

半颗麦芽糖

MeteoInfo 发表于 2016-12-24 23:12
等值线分析需要二维数组。

请教一下王老师这句h = h[::144]是什么意思？

MeteoInfo

半颗麦芽糖 发表于 2016-12-25 15:35
请教一下王老师这句h = h[::144]是什么意思？

间隔144取数组要素，其实就是将数组抽稀。目的是取出没有重复的高度数组。

半颗麦芽糖

MeteoInfo 发表于 2016-12-25 18:26
间隔144取数组要素，其实就是将数组抽稀。目的是取出没有重复的高度数组。

多谢王老师指导，学习了

zhksoft

crainkuang

帖主，能分享一下你的txt吗？没钱下不了啊。如果能够附上你做的图更好。因为我想利用这个软件做一个大尺度格局的风场时空分布图。邮箱是 kuangfl@126.com。翘首以盼，谢谢了！