1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
|
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 2004-2006 - INRIA - Farid Belahcene
// This file must be used under the terms of the CeCILL.
// This source file is licensed as described in the file COPYING, which
// you should have received as part of this distribution. The terms
// are also available at
// http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt
function pie(varargin)
// Copyright INRIA
// This function draws a pie, if size of x is N then pie function draws a pie with N parts, the area of the ith part is equal to (x(i)/sum(x))*( surface of the unit cercle).
//
//
// syntax : pie(x[,sp[,txt]])
// Input : The input arguments must have the same size
// x : a scalar or a vector of positive reals.
// sp : a scalar or a vector of reals. The sp vector allows to cut one or several parts of the pie
// txt : a cell or a vector of strings. The txt vector allows to write a text for each part of the pie
varlist = varargin;
// Input arguments must have the same length
for i=2:size(varlist)
if size(varlist(i-1),"*") <> size(varlist(i),"*") then
error(msprintf(gettext("%s: Wrong size for input arguments: Matrices of same size expected.\n"), "pie"));
end
end
txt = [];
esp = [];
// Detect and set the RHS arguments
// case : pie(x)
if size(varlist) == 1 then
if type(varlist(1))==1 & and(varlist(1)>0) then
x = varlist(1)
else
error(msprintf(gettext("%s: Wrong type for input argument #%d: A real scalar or vector expected.\n"), "pie", 1));
end
// case : pie(x,sp) or pie(x,txt)
elseif size(varlist) == 2 then
if type(varlist(1)) == 1 & and(varlist(1)>0) then
x = varlist(1);
else
error(msprintf(gettext("%s: Wrong type for input argument #%d: A real scalar or vector expected.\n"), "pie", 1));
end
if type(varlist(2)) == 1 | type(varlist(2)) == 4
esp = varlist(2);
elseif type(varlist(2)) == 10 | iscellstr(varlist(2)) then
txt = varlist(2);
end
// case : pie(x,sp,txt)
elseif size(varlist) == 3 then
if type(varlist(1)) == 1 & and(varlist(1)>0) then
x = varlist(1);
else
error(msprintf(gettext("%s: Wrong type for input argument #%d: A real scalar or vector expected.\n"), "pie", 1));
end
if type(varlist(2)) == 1 | type(varlist(2)) == 4
esp = varlist(2);
else
error(msprintf(gettext("%s: Wrong type for input argument #%d: A real scalar or vector expected.\n"), "pie", 2));
end
if type(varlist(3)) == 10
txt = varlist(3);
elseif iscellstr(varlist(3)) then
for j=1:size(varlist(3),"*")
txt(j) = varlist(3).entries(j);
end
else
error(msprintf(gettext("%s: Wrong size for input argument #%d: Vector of strings expected.\n"),"pie", 3));
end
else
error(msprintf(gettext("%s: Wrong number of input argument(s): %d to %d expected.\n"), "pie", 1, 3));
end
if isempty(x) then
error(msprintf(gettext("%s: Wrong type for input argument #%d: A real scalar or vector expected.\n"), "pie", 1));
end
// xi and yi represents the coordinates of each polyline
// iesp is the index of the part of i which are separated of the pie
// teta_1 is the start angle of the arc polyline, teta_2 is the end angle of the arc polyline
iesp = find(esp>0 | esp<0);
teta_2 = 0 + %pi/2;
yi = [];
xi = [];
CurColor = 0;
drawlater();
// Create a close polyline for every parts of pie, the polyline inside color is determinated by the plot colormap
for i=1:size(x,"*")
xi = [];
yi = [];
teta_1 = teta_2;
teta_i = (x(i)/sum(x))*2*%pi;
teta_2 = teta_1 + teta_i
if size(x,"*") <> 1 then
xi(1) = 0;
yi(1) = 0;
xi(2) = cos(teta_1);
yi(2) = sin(teta_1);
else
xi(1) = cos(teta_1);
yi(1) = sin(teta_1);
end
inter = 1/(100*x(i)/sum(x));
k = inter;
while k<1
xi($+1) = cos((1-k)*teta_1 +k*teta_2);
yi($+1) = sin((1-k)*teta_1 +k*teta_2);
k= k + inter;
end
xi($+1) = cos(teta_2);
yi($+1) = sin(teta_2);
xfpolys(xi,yi);
e = gce();
a = gca();
ei = e.children;
if or(i == iesp) then
ei.x_shift = ones(1,size(xi,"*")) * (1/10) * cos((teta_2+teta_1)/2);
ei.y_shift = ones(1,size(yi,"*")) * (1/10) * sin((teta_2+teta_1)/2);
if txt <> []
xstring(cos((teta_2+teta_1)/2)*1.2+ei.x_shift(1)-0.1*(cos((teta_2+teta_1)/2)<0),sin((teta_2+teta_1)/2)*1.2+ei.y_shift(1),txt(i));
else
xstring(cos((teta_2+teta_1)/2)*1.2+ei.x_shift(1)-0.1*(cos((teta_2+teta_1)/2)<0),sin((teta_2+teta_1)/2)*1.2+ei.y_shift(1),string(round((x(i)/sum(x))*100)) + "%");
end
else
if txt <> []
xstring(cos((teta_2+teta_1)/2)*1.2-0.1*(cos((teta_2+teta_1)/2)<0),sin((teta_2+teta_1)/2)*1.2,txt(i));
else
xstring(cos((teta_2+teta_1)/2)*1.2-0.1*(cos((teta_2+teta_1)/2)<0),sin((teta_2+teta_1)/2)*1.2,string(round((x(i)/sum(x))*100)) + "%");
end
end
[Color,CurColor] = setDefaultColor(CurColor);
ei.background = Color;
ei.fill_mode = "on";
a.isoview = "on";
a.box = "off";
a.axes_visible = "off";
// Update data_bounds
a.data_bounds = [min(-1.3,a.data_bounds(1,1)) min(-1.3,a.data_bounds(1,2));max(1.3,a.data_bounds(2,1)) max(1.3,a.data_bounds(2,2))];
end
drawnow();
endfunction
|
