[go: up one dir, main page]
Menu ▾ ▴

[57f2a1]: / anal / pp_merge.m  Maximize  Restore  History

Download this file

152 lines (119 with data), 4.7 kB 

  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
148
149
150
151
152
% Merge power profiles with similar resolution r_w



function [r_pprange_merged,r_pperr_merged,r_ppw_merged,r_pp_merged,r_pprofile_id] = pp_merge(r_pprange,r_pperr,r_ppw,r_pp,name_expr)



if nargin<2, error('Not enough input parameters.'); end



if nargin<5, name_expr = []; end

if nargin<4, warning('r_pp not as input --> r_pp_merged = Nan vector'), r_pp = []; end

if nargin<3, r_ppw = []; end



if isempty(r_pp), r_pp = NaN(length(r_pprange),1); end



% figure(1)

% subplot(3,2,1)

% plot(r_pprange,'.')

% ylabel('r\_pprange')

% title('Original profiles')

% subplot(3,2,3)

% plot(r_ppw,'.')

% ylabel('r\_ppw')



if isempty(find(diff(r_pprange)<0,1)) && isempty(find(abs(diff(r_pprange))>10*mean(abs(diff(r_pprange))),1))  % if there are no negative or large steps in the range, no merging is needed

    r_pp_merged      = r_pp;

    r_pperr_merged   = r_pperr;

    r_pprange_merged = r_pprange;

    r_ppw_merged     = r_ppw; 

    r_pprofile_id    = ones(length(r_pprange),1); 

    return

end



if isempty(r_ppw) && isempty(r_pperr)

    r_pp_merged      = r_pp;

    r_pperr_merged   = r_pperr;

    r_pprange_merged = r_pprange;

    r_ppw_merged     = r_ppw; 

    r_pprofile_id    = ones(length(r_pprange),1); 

    return

end





if ~isempty(r_ppw)

    

    rres   = max(1,min(r_ppw)); %range resolution given by p_dtau or smallest volume

    r_ww   = unique(round(r_ppw/rres));



    if isempty(diff(r_ww))    % only one (1) value in r_ww

        rres = 1;

        r_dw = 1;

    else

        r_dw   = unique(diff(r_ww))+1;

    end



    rres_x = round(r_ppw/rres/r_dw(1));

    ranges = round(r_pprange/rres);%+1e-3*rres_x;

    r_pp_merged      = [];

    r_pperr_merged   = [];

    r_ppw_merged     = [];

    r_pprange_merged = [];

    r_pprofile_id    = [];



    n_profiles = length(unique(rres_x));                 % number of profiles after merging

    merge_pprofiles = sort(unique(rres_x),'descend');    % sort from high to low 'r_ppw widths'



    id = 0;



    for kk = 1:n_profiles



        ii = find(rres_x==merge_pprofiles(kk));

        r_pp_tmp      = r_pp(ii);

        r_pperr_tmp   = r_pperr(ii);

        r_ppw_tmp     = r_ppw(ii);

        r_pprange_tmp = r_pprange(ii);

        ranges_tmp    = ranges(ii);

    

        a3 = [];

    

        for range=unique(ranges_tmp)'    % must be a row vector    

            aa = find(ranges_tmp==range);

    

            if length(aa)>1

       

                w2=1./r_pperr_tmp(aa).^2;

       

                r_pp_tmp(aa(1))      = sum(r_pp_tmp(aa).*w2)/sum(w2);

                r_pperr_tmp(aa(1))   = 1/sqrt(sum(w2));

                r_ppw_tmp(aa(1))     = sum(r_ppw_tmp(aa).*w2)/sum(w2);

                r_pprange_tmp(aa(1)) = sum(r_pprange_tmp(aa).*w2)/sum(w2);

       

                a3 = [a3 aa(2:end)'];

            end

        end

    

        r_pp_tmp(a3)      = [];

        r_pperr_tmp(a3)   = [];

        r_pprange_tmp(a3) = [];

        r_ppw_tmp(a3)     = []; 

    

        [r_pprange_sorted,jj] = sort(r_pprange_tmp);

        r_pp_sorted           = r_pp_tmp(jj);  

        r_pperr_sorted        = r_pperr_tmp(jj);  

        r_ppw_sorted          = r_ppw_tmp(jj);  

    

        % Also separate if there is a big jump in the range:

    

        diffrange_tmp = diff(r_pprange_sorted);

        mean_diff     = mean(abs(diffrange_tmp));

        aa = find(abs(diffrange_tmp)>10*mean_diff);                 % find sudden big steps in the range 

    

        if isempty(aa)

            id = id + 1;

            r_pprofile_id = [r_pprofile_id id*ones(1,length(r_ppw_sorted))];

        else

            bb = [0 aa' length(r_ppw_sorted)];

            for gg = 1:(length(aa)+1)

                id = id + 1;

                r_pprofile_id = [r_pprofile_id id*ones(1,bb(gg+1)-bb(gg))];

            end

        end

   

        r_pprange_merged = [r_pprange_merged r_pprange_sorted'];

        r_pp_merged      = [r_pp_merged r_pp_sorted'];

        r_pperr_merged   = [r_pperr_merged r_pperr_sorted'];

        %r_ppw_merged     = [r_ppw_merged r_ppw_sorted(jj)']; 

        r_ppw_merged     = [r_ppw_merged r_ppw_sorted']; 

   

    end

    r_pp_merged      = r_pp_merged';

    r_pperr_merged   = r_pperr_merged';

    r_pprange_merged = r_pprange_merged';

    r_ppw_merged     = r_ppw_merged'; 

    r_pprofile_id    = r_pprofile_id'; 



elseif isempty(r_ppw) && ~isempty(name_expr)



    % special merging

    

else

    error('At least either r_ppw or name_expr is needed as input')

end



% figure(1)

% subplot(3,2,2)

% plot(r_pprange_merged,'.b')

% title('Merged profiles')

% ylabel('merged r\_pprange')

% subplot(3,2,4)

% plot(r_ppw_merged,'.b')

% ylabel('merged r\_ppw')

% subplot(3,2,6)

% plot(r_pprofile_id,'.b')

% ylabel('power profile id')

% ylim([0 id+1])