findFootprintArea.m

function [xmin,ymin,xmax,ymax] = findFootprintArea(I)
% FINDFOOTPRINTAREA finds the approximated Rectangle around the footprint
%                   
% Syntax:
%    [xmin,ymin,xmax,ymax] = findFootprintArea(I)


%% Settings
% Puffer um auch gesamten Bereich zu erhalten. (da Rauschfilter den verkleinert)
puffer=90;
%% load Image
img_a=I;
%  figure,imshow(img_a);

BW=createMask_pink(img_a);
%   imshow(BW);

%% Filtere Rauschen 

seD = strel('diamond',2);
BWfinal = imerode(BW,seD);
BWfinal = imerode(BWfinal,seD);
%    imshow(BWfinal), title('segmented image');
 
seD = strel('diamond',1);
BWfinal = imerode(BWfinal,seD);
BWfinal = imerode(BWfinal,seD);
BWfinal = imerode(BWfinal,seD);
BWfinal = imerode(BWfinal,seD);
BWfinal = imerode(BWfinal,seD);
BWfinal = imerode(BWfinal,seD);
BWfinal = imerode(BWfinal,seD);
BWfinal = imerode(BWfinal,seD);
BWfinal = imerode(BWfinal,seD);
BWfinal = imerode(BWfinal,seD);
BWfinal = imerode(BWfinal,seD);
BWfinal = imerode(BWfinal,seD);
  imshow(BWfinal), title('3x segmented image');


%% Bestimme Umriss
[x,y]=find(BWfinal==true);
xmax=max(y)+puffer;
xmin=min(y)-puffer-60; % xmin manuell angepasst
ymax=max(x)+puffer;
ymin=min(x)-puffer;

 rectangle('Position',[xmin,ymin,xmax-xmin+puffer,ymax-ymin+puffer],'EdgeColor','b','LineWidth',3);
end