[原创]MATLAB斜温图、高空测风绘制函数

沙颖凯

这里有Skew-T的最新版本，相比于之前的内容有不少优化： [http://bbs.06climate.com/forum.p ... ad&tid=27252&extra=] 本帖已不再更新。
______________________________________________________________________________________________
刚刚完成，新鲜出炉的斜温图（T-log P）和单站测风图的MATLAB程序，按列输入探空数据即可使用，tskew_plot绘制斜温图，wind_polar绘制单站测风图。

tskew_plot.m (6.92 KB, 下载次数: 255)

2014-7-26 21:42 上传

点击文件名下载附件

wind_polar.m (1.84 KB, 下载次数: 199)

2014-7-26 21:42 上传

点击文件名下载附件

windbarb_profile.m (6.4 KB, 下载次数: 190)

2014-7-26 21:42 上传

点击文件名下载附件



tskew_plot和wind_polar参考MIT Open Course Ware使用的程序，以及GrADS的斜温图函数plotskew.gs编写的。目的就是在MATLAB环境下实现plotskew.gs的效果。

考虑到实际应用的灵活性，这里把单站测风和斜温图的脚本分开了。wind_polar可以设定测风数据的层次，默认输出1000, 850, 700, 500, 300hPa的数据，若提供的数据不包含这些层次（一般都包含）会默认线性插值到该层次。

tskew_plot画风廓线的时候会用到windbarb_profile.m，我不记得这是谁写的了，大概是在Research Gate的某个地方找到的。总之，脚本里面列出了主要参考的东西，也感谢所有这些参考程序的贡献者。

以下作一些简单的介绍，最后会附上测试代码，数据和效果图。

tskew_plot(pz, tz, rhz, wdir, wsp, Filter)
% ====================================================================== %
%   pz: pressure, hPa%   tz: Temperature, degC
%   rhz: Relative Humidity,  (0-1);%   wdir: Wind Direction, deg (0-360);
%   wsp: Wind Speed, m/s
%   Filter: Apply 1-D digital filter to smooth the wind obs., 1 = do;
%       0 = not do, 0 is default. See MATLAB function "filter.m".
% ======================================================================= %

wind_polar(pz, wdir, wsp, level)
% Polar plot for Upper-Air wind monitoring/sounding
% ======================================================================= %
%   pz: pressure, hPa
%   wdir: Wind Direction, deg (0-360);
%   wsp: Wind Speed, m/s
% ======================================================================= %
测试数据
soundings_RF6.mat (19.94 KB, 下载次数: 130)

2014-7-26 21:43 上传

点击文件名下载附件

Sounding Data via U. Wyoming[http://weather.uwyo.edu/upperair/sounding.html]

1. clear all
   
2. close all
   
3. %% Import data
   
4. load('soundings_RF6.mat')
   
5. sounding=sound_oak.rf61700Z;
   
6. pz=sounding(:, 1); % Pressure (hPa)
   
7. tz=sounding(:, 3); % Temperature (degC)
   
8. tdz=sounding(:, 4);% Dew-Point Temperature (degC)
   
9. rhz=sounding(:, 5)/100; % Relative Humidity
   
10. wdir=sounding(:, 7); % Wind Direction (deg)
    
11. wsp=sounding(:, 8)*0.514; % Wind Speed (m/s)
    
12. %% Plot
    
13. fig1=figure('Color', 'w');
    
14. title('\fontsize{8}Skew-T Plot', 'FontWeight', 'bold')
    
15. tskew_plot(pz, tz, rhz, wdir, wsp, 0);
    
16. %fig2=figure('Color', 'w');
    
17. %title('\fontsize{8}Upper-wind Polar Plot', 'FontWeight', 'bold')
    
18. %wind_polar(pz, wdir, wsp)

复制代码

结果
斜温图：

单站测风：

_____________________________________________________________________________________
2014/7/28
想把斜温图里风向标改在上面的，可自行修改"windbarb_profile"中变量y2在不同风向角的取值，也可参考Mapping Toolbox中风向标的函数[http://www.mathworks.com/matlabc ... 1-wind-barb-plotter]自己重写一个。

freekiller

不错，谢谢分享。 Matlab绘制风向杆的时候是向下的，儿一般习惯是短线在上。这个问题我一直没解决

沙颖凯

看在我不收贡献+不设权限+原创+代码公开+数据齐全的份上，求加精，求奖励贡献，我还有好多书想下呢，贡献快到碗里来 (ΦωΦ)。

顺便推销自己的帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=25219 同样求贡献

沙颖凯

freekiller 发表于 2014-7-26 23:11
不错，谢谢分享。 Matlab绘制风向杆的时候是向下的，儿一般习惯是短线在上。这个问题我一直没解决

没有这种说法，这里用的是逆时针，东风向上，GrADS用的好象是顺时针，就是西风向上。风羽一直向上的用不了。

freekiller

沙颖凯 发表于 2014-7-27 11:00
没有这种说法，这里用的是逆时针，东风向上，GrADS用的好象是顺时针，就是西风向上。风羽一直向上的用不 ...

我不需要说法，你自己看呗

沙颖凯

freekiller 发表于 2014-7-27 11:33
我不需要说法，你自己看呗

这不就是西风向上东风向下的么，我习惯叫这种是顺时针的。短线一直在上你怎么画气旋。
不过我觉得无所谓，顺时针和逆时针不影响图的效果。

freekiller

沙颖凯 发表于 2014-7-27 11:56
这不就是西风向上东风向下的么，我习惯叫这种是顺时针的。短线一直在上你怎么画气旋。
不过我觉得无所谓 ...

你没明白我的意思，我的意思是你西风bar应向上
           

沙颖凯

freekiller 发表于 2014-7-27 12:05
你没明白我的意思，我的意思是你西风bar应向上

把windbarb_profile.m里面"line([x2 x3],[y2 y3],'linewidth',width,'color',color)"这种语句的x和y的终点修改下，应该是把计算距离的符号变成和原来相反的。

大山与大海

厉害厉害　呵呵　

不一样的生活

好东西，共分享。。。。。。