将矢量(或弧)绘制到玫瑰图上。 MATLAB

时间:2014-05-14 17:03:41

标签: matlab plot rose-diagram

我有两个数据集。一个详细列出angles(我正在绘制玫瑰图):

angles
-0.8481065519
0.0367932161
2.6273740453
...
n

另一个,从这组角度详述方向统计:

angle,error
-0.848106563,0.8452778824

angle基本上定义方向均值,error是圆形方差,angle

的任何一侧基本上都是错误栏

到目前为止,我已经使用这组角度绘制了一个玫瑰直方图,如下:

h = rose(angles,36)

我想创建一个方向统计量angle的图(它不需要长度/幅度 - 只需要圆圈图的边缘),周围有error。 举个例子:

An example of what I want to do.

我在Matlab手工添加了这些行。如果可能的话,也许在弧内也可能有阴影。或者,(并且可能是优选的)将在玫瑰图箱上方(因此它不覆盖数据)只有一条条带中心线(显示angle和{{1}的阴影周围}。

提前致谢。

1 个答案:

答案 0 :(得分:3)

这个怎么样?

%// Data
angles = 2*pi*.8*randn(1,1e4);
angle = -0.848106563;
error = 0.8452778824;

%// Plot rose
rose(angles, 36);
axis image %// make axis square
hold on

%// Plot mean
a = axis;
a = a(2); %// size of axis
plot([0 cos(angle)*a], [0 sin(angle)*a], 'r')

%// Plot error as many shaded triangles that compose a circular wedge
t = linspace(-error/2+angle,error/2+angle,100); %// increase "100" if needed
for k = 1:numel(t)-1
    h = patch([0 cos(t(k))*a cos(t(k+1))*a 0], ...
        [0 sin(t(k))*a sin(t(k+1))*a 0], [.5 0 0], 'edgecolor', 'none');
        %// change color [.5 0 0] to something else if desired. Note also alpha
    set(h,'Facealpha',.3) %// make transparent
end     

%// Place rose on top by rearranging order of axis children
ch = get(gca,'children');
set(gca,'children',[ch(2:end); ch(1)]);

enter image description here

为此,您需要使用具有透明度的图形渲染器。因此,您可能需要调整图形的renderer属性。