Shape/cornerstat.m at master · nbul/Shape · 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
%% Corner distribution
X2 = 3; % binning interval for corners

image7 = figure;
corner_rot2=corner_rot2(~cellfun(@isempty, corner_rot2));
for ll=1:length(corner_rot2)
    if corner_rot2{ll}
        scatter(corner_rot2{ll}(:,1),corner_rot2{ll}(:,2));
        hold on
    end
end
hold off

print(image7, '-dtiff', '-r300', 'all_corners.tif');

%% Distance and angle for each corner
corner_angle = zeros(1,1);
distance_corner = zeros(1,1);
pp = 0;

for l=1:length(corner_rot2)
    for m=1:length(corner_rot2{l})
        pp=pp+1;
        distance_corner(pp) = sqrt(corner_rot2{l}(m,2)*corner_rot2{l}(m,2) +...
            corner_rot2{l}(m,1)*corner_rot2{l}(m,1));

        a = corner_rot2{l}(m,2)/distance_corner(pp);
        b = corner_rot2{l}(m,1)/distance_corner(pp);
        corner_angle(pp)= atan2d(a,b);
        if corner_angle(pp)<0
            corner_angle(pp)= corner_angle(pp)+360;
        end
    end
end

%% Number of corners by angle
number_corners = zeros(1,360/X2);
for j=X2:X2:360
    for q=1:length(corner_angle)
        if corner_angle(q)<=j && corner_angle(q)>j-X2
            number_corners(j/X2)=number_corners(j/X2)+1;
        end
    end
end
image8=figure;
plot(X2:X2:360,number_corners, 'LineWidth', 2);
set(gca,'FontSize',18);
ylim([0 max(number_corners)*1.1]);
xlim([0 400]);
print(image8, '-dtiff', '-r300', 'number_corners.tif');

%% Distance of corners by angleN= zeros(360, 100);
curve_corners = struct([]);
gof_corners = struct([]);
distance_temp = struct([]);

radius_corners = zeros(1, 360/X2);
sigma1_corners = zeros(1, 360/X2);
pvalue_corners = zeros(1, 360/X2);
N_corners= zeros(X2, 100);

for i=X2:X2:360
    distance_temp{i/X2}= zeros(1,1);
    for l=1:length(distance_corner)
        if corner_angle(l)<=i && corner_angle(l)>i-X2
            distance_temp{i/X2} = [distance_temp{i/X2}; distance_corner(l)];
        end
    end
    distance_temp{i/X2} =  distance_temp{i/X2}(distance_temp{i/X2}~=0);
    [N_corners(i/X2,:), bins] = histc(distance_temp{i/X2},binrange);
    [curve_corners{i/X2},gof_corners{i/X2}] = fit(bincenter',N_corners(i/X2,:)','gauss1');
    radius_corners(i/X2) = curve_corners{i/X2}.b1;
    sigma1_corners(i/X2) = sqrt(curve_corners{i/X2}.c1*curve_corners{i/X2}.c1/2);
    pvalue_corners(i/X2) = gof_corners{i/X2}.rsquare;
end
%% Radius of corners circle - radius_corners
image9=figure;
cd(sum_dir);
plot(X2:X2:360, radius_corners, 'LineWidth', 2);
set(gca,'FontSize',18);
ylim([0 max(radius_corners)*1.1]);
xlim([0 400]);
print(image9, '-dtiff', '-r300', 'corners_radius.tif');
%% Number of corners as a function of radius
image10=figure;
radius_vs_corner = [radius', number_corners'];
radius_vs_corner = sortrows(radius_vs_corner,1);
[f_rc,gof_rc] = fit(radius_vs_corner(:,1),radius_vs_corner(:,2),'exp1');

number_corners_fit = f_rc.a*exp(f_rc.b*(min(radius_vs_corner(:,1)-1):1:(max(radius_vs_corner(:,1))+1)));
plot((min(radius_vs_corner(:,1)-1):1:(max(radius_vs_corner(:,1))+1)),number_corners_fit, 'r', 'LineWidth',2);
hold on;
scatter(radius_vs_corner(:,1),radius_vs_corner(:,2), 'b');
set(gca,'FontSize',18);
ylim([0 max(radius_vs_corner(:,2))*1.1]);
xlim([0 max(radius_vs_corner(:,1))*1.1]);
hold on
annotation('textbox',[.2 .8 .0 .0],'String',{['R_s_q_u_a_r_e = ', num2str(gof_rc.rsquare)],...
    ['f = ', num2str(f_rc.a), ' * exp(', num2str(f_rc.b), ' * x)'] },'FitBoxToText','on','FontSize', 14);
hold off;
print(image10, '-dtiff', '-r300', 'radius_vs_corners.tif');