MRO projekt 3:46:45am dzień 1

img = imread('test7.png');
bw = im2bw(img, 0.9);
bw = ~bw;

% matlabpool open
rots = 1:16; % todo: zwiększyć do 1:360. Niech się komp męczy, a co!
rot2area = zeros(size(rots));
rotang=360/size(rots,2);
parfor i=rots
    imm = imrotate(bw, i*rotang);

    % szumy won
    L = bwlabel(imm);
    s = regionprops(L, 'Area');
    area_values = [s.Area];
    idx = find(5 < area_values);
    imm = ismember(L, idx);

    % pole otoczki
    [r,c] = find(imm);
    height = max(r)-min(r);
    width = max(c)-min(c);
    rot2area(i) = height * width;
end

% dobieramy najlepsze obroty

[~,rots] = sort(rot2area);
rots = rots(1:4);

imm = imrotate(bw, rots(4)*rotang);

% odszumiamy raz jeszcze
L = bwlabel(imm);
s = regionprops(L, 'Area');
area_values = [s.Area];
idx = find(5 < area_values);
imm = ismember(L, idx);

imshow(imm)

% % szumy won
% L = bwlabel(bw);
% s = regionprops(L, 'Area');
% area_values = [s.Area];
% idx = find(5 < area_values);
% bw = ismember(L, idx);
% imshow(bw)
%
% % oblicz otoczkę wszystkiego
% [r,c] = find(bw);
% height = max(r)-min(r);
% width = max(c)-min(c);
% area = height * width;

%rots = 1:2;
%t_im1rs = zeros(size(rots,2),1);
%for angle = rots
%    t_im1rs(angle) = 0; %imrotate(t_im1, angle);
%    tmp = regionprops(imrotate(im, angle), 'BoundingBox');
%    [ul_corner, width] = tmp.BoundingBox
%end