Wolbachia_malaria_AWM/plotDash.m at main · zhuolinqu/Wolbachia_malaria_AWM · GitHub

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
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%       Before you use call this function,
%       I recommand you to set your figure's xlim and ylim (i.e., axis range)
%       and size of the figure
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [p] = plotDash(h, x, y, spec_, num_, varargin)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   h: Handle of the figure
%       ex) h = figure(1)
%   x: data x
%   y: data y
%   spec_: spec of line - (relative length of)(solid, blank, ...)
%          Default option: -1 --> dashed line
%       ex) spec_ = [1, 0.5] || spec_ =  [10, 5] || ...
%       ex) spec_ = [1, 0.1] || ...
%       ex) spec_ = [1, 0.1, 0.1, 0.1] || ...
%   num_: How many times that pattern will be repeated for xlim or ylim
%   vargargin: "Pair(s)" of options for plot
%       ex) varargin{1} = 'color', varargin{2} = 'g'

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   User input - Points per pattern
    n_node = 50;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    %   Size of x
    nx = length(x);
    ny = length(y);
    if (nx ~= ny)
        p = plot(0, 0);
        delete(p);
        return;
    else
        n = nx;
    end

    %   Size of plot
    tmp = get(h, 'InnerPosition');
    W = tmp(3);
    D = tmp(4);
    SF = W / 560;
    R  = D / W * 28.55 / 27.1;
    if (sum(size(h.CurrentAxes)) == 0)
        p = plot(0, 0);
        delete(p);
        DX = max(x) - min(x);
        DY = max(y) - min(y);
    else
        tmp = h.CurrentAxes.XLim;
        DX = tmp(2) - tmp(1);
        tmp = h.CurrentAxes.YLim;
        DY = tmp(2) - tmp(1);
        clear tmp
    end

    %   Normalization
    r = zeros(n, 2);
    for i = 1:n
        r(i, 1) = x(i) / DX * SF;
        r(i, 2) = y(i) / DY * SF * R;
    end

    %   Curve length
    s  = zeros(n    , 1);
    ds = zeros(n - 1, 1);
    for i = 1:(n - 1)
        ds(i)    = norm(r(i + 1, :) - r(i, :));
        s(i + 1) = s(i) + ds(i);
    end

    %   Linearly splined data - User input: 50
    ds_min = min(ds);
    ds_min = min(ds_min, 1/(n_node*num_));
    s_s = 0:ds_min:s(end);
    x_s = interp1(s, x, s_s);
    y_s = interp1(s, y, s_s);

    %   Default line spec = dashed
    if (spec_(1) < 0)
        spec_ = [10, 6];
        num_ = 26;
    end

    %   spec_: absolute --> relative
    spec_ = spec_ / sum(spec_) / num_;
    m = length(spec_);
    l_spec = sum(spec_);

    %   Set T/F flag for index
    k = length(s_s);
    flag_idx = true(k, 1);
    for i = 1:k
        curr = rem(s_s(i), l_spec);
        for j = 2:2:m
            flag = (curr > sum(spec_(1:(j-1))));
            flag = (curr < sum(spec_(1:(j-0)))) && flag;

            if (flag)
                flag_idx(i) = false;
            end
        end
        flag_idx(i);
    end

    %   Set NaN to x_s and y_s
    x_s(~flag_idx) = NaN;
    y_s(~flag_idx) = NaN;

    %   Plot graph
    p = plot(h.CurrentAxes, x_s, y_s,'Marker','none');
    %   Apply nargin
    for i = 1:2:(nargin - 5)
        set(p, varargin{i}, varargin{i + 1});
    end

end