speckimg/mergeimages.m at main · comediaLKB/speckimg · 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
function [Reached_Pattern,Global_image] = mergeimages(Maximum_intensity,O,EstimatedNumber,xx,yy,xpixel,ypixel,num)
%% contact: leizhu201806@hotmail.com // version 07/2021
%%
Maximum_intensity = Maximum_intensity;% the maximum value of each image
O = O;% the partial images
EstimatedNumber = EstimatedNumber;% the estimated rank
xx = xx;% the position of maximum value of each image in X coordinate
yy = yy;% the position of maximum value of each image in Y coordinate
xpixel = xpixel;% the size of image in X coordinate
ypixel = ypixel;% the size of image in Y coordinate
% close all
Normlized_Maximum_intensity = zeros(size(Maximum_intensity));
Size = 500;
Global_image = zeros(Size,Size); %  Create the empty matrix
Central_Xposition = zeros(EstimatedNumber,1); %  Create the empty matrix to store the position of each emitter in X coordinate
Central_Yposition = zeros(EstimatedNumber,1); %  Create the empty matrix to store the position of each emitter in Y coordinate
ff = 0;
for k=1:EstimatedNumber
    if k == 1
        num = num;
        SNormlized_Maximum_intensity = Maximum_intensity(:,num)/max(Maximum_intensity(:,num));
        [aa,indices] = sort(SNormlized_Maximum_intensity,'descend');
        Known_Finger = indices;
        Reached_Pattern(k) = Known_Finger(1);
        for zz = 1:1
            % calculate the relative postion
            global_image = ones(size(Global_image));
            Central_Xposition(Reached_Pattern(k)) = floor(Size/2)+xx(Reached_Pattern(k),Reached_Pattern(k))...
                -xx(Reached_Pattern(k),Reached_Pattern(k))-1;
            Central_Yposition(Reached_Pattern(k)) = floor(Size/2)+yy(Reached_Pattern(k),Reached_Pattern(k))...
                -yy(Reached_Pattern(k),Reached_Pattern(k))-1;
            hh = histogram(O{Reached_Pattern(k)},100);
            zz = hh.Values;
            background = hh.BinEdges(find(zz == max(zz))+1);
            global_image = global_image*background;
            global_image(Central_Xposition(Reached_Pattern(k))+(-floor(xpixel/2):floor(xpixel/2)-1),...
                Central_Yposition(Reached_Pattern(k))+(-floor(ypixel/2):floor(ypixel/2)-1))...
                = O{Reached_Pattern(k)};
            Global_image = Global_image + global_image; % merge them into the blobal image.
            ff = ff+1
        end
%         figure
%         imshow(Global_image,[]);colormap hot;
    else
        for zz = 2:(EstimatedNumber)
            NextPSF = indices(zz,1);
            check = ismember(Reached_Pattern,NextPSF);
            if check == 0
                Reached_Pattern(k) = NextPSF;
                break
            elseif check == 1
            end
        end
        if SNormlized_Maximum_intensity(NextPSF)<= 0.01
            Reached_Pattern(k) = [];
            break
        end
        global_image = ones(size(Global_image));
        temp = Reached_Pattern(k-1);
        Central_Xposition(Reached_Pattern(k)) = Central_Xposition(temp)-...
            0*floor(Size/2)+xx(Reached_Pattern(k),temp)-xx(temp,temp);
        Central_Yposition(Reached_Pattern(k)) = Central_Yposition(temp)-...
            0*floor(Size/2)+yy(Reached_Pattern(k),temp)-yy(temp,temp);
        hh = histogram(O{Reached_Pattern(k)},100);
        zz = hh.Values;
        background = hh.BinEdges(find(zz == max(zz))+1);
        global_image = global_image*background;
        global_image((Central_Xposition(Reached_Pattern(k))-floor(xpixel/2)):(Central_Xposition(Reached_Pattern(k))+floor(xpixel/2)-1),...
            (Central_Yposition(Reached_Pattern(k))-floor(ypixel/2)):(Central_Yposition(Reached_Pattern(k))+floor(ypixel/2))-1)...
            = O{Reached_Pattern(k)};
        Global_image = Global_image + global_image; % merge them into the blobal image.
        ff = ff+1
        imshow(Global_image,[]);colormap hot;
        SNormlized_Maximum_intensity = Maximum_intensity(:,Reached_Pattern(k))/max(Maximum_intensity(:,Reached_Pattern(k)));
        [aa,indices] = sort(SNormlized_Maximum_intensity,'descend');
    end
end