-
Notifications
You must be signed in to change notification settings - Fork 185
/
quiver.m
69 lines (51 loc) · 1.48 KB
/
quiver.m
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
function h = quiver(grains,dir,varargin)
% plot directions at grain centers
%
% Syntax
% quiver(grains,dir,'color','r')
%
% mtexdata forsterite
% grains = calcGrains(ebsd('indexed'))
% quiver(grains,grains.meanRotation.axis,'color','r')
%
% Input
% grains - @grain2d
% dir - @vector3d
%
% Options
% antipodal - plot directions or axes
% maxHeadSize - size of the arrow
%
pos = grains.centroid;
if ~check_option(varargin,'noScaling')
dir = 0.2*grains.diameter .* normalize(dir) * ...
get_option(varargin,'autoScaleFactor',1);
varargin = ['linewidth',2,'autoScale','off',varargin];
else
varargin = ['linewidth',2,'autoScaleFactor',0.25,varargin];
end
% get color
c = get_option(varargin,'color');
if check_option(varargin,'antipodal') || dir.antipodal
varargin = [{'MaxHeadSize',0},varargin];
% double color;
if isnumeric(c) && length(c) == length(dir), c = [c;c]; end
pos = [pos;pos];
dir = [dir(:);-dir(:)];
else
%varargin = [{'MaxHeadSize',5},varargin];
end
% if different color are given - separate them
[c,~,id] = unique(c,'rows');
if length(id) == length(dir)
varargin = delete_option(varargin,'color',1);
for i = 2:size(c,1)
hold on
h(i) = optiondraw(quiver3(pos.x(id == i),pos.y(id == i),pos.z(id == i),...
dir.x(id == i),dir.y(id == i),dir.z(id == i)),varargin{:},'color',c(i,:)); %#ok<AGROW>
end
else
h = optiondraw(quiver3(pos.x,pos.y,pos.z,dir.x,dir.y,dir.z),varargin{:});
end
if nargout == 0, clear h; end
end