close all; clear all; clc;

a = imread('coins.png');

bin = a>80;
SE = strel('disk',3);
bin = imopen(bin, SE);

aseg = bwlabel(bin);
N = max(aseg(:));
props = regionprops(bin,a,'all');

bez = zeros(size(a));
nowy = zeros(size(a));

for k=1:N
   temp = (aseg==k);

   if( props(k).MajorAxisLength > 51 && props(k).MinorAxisLength>51 )
      bez = bez | temp;
      mediana = median(props(k).PixelValues);
      odch = double((props(k).PixelValues));
      odchylenie = std(odch);

      for y=1:246
          for x=1:300
          end
      end


   end

end

subplot(121), imshow(a);
subplot(122), imshow(bez);