find_cross_xy.m

function [x,y] = find_cross_xy(I,ROI)
% FIND_CROSS_XY returns the estimated coordinate of the hand drawn
% cross
%
% Syntax:
%    [x,y] = find_cross_xy(I,ROI)
%
% Inputs:
%    I - Image I
%    ROI - Rectange Region of Interest (where the cross is likely to be)
%
% Outputs:
%    x,y coords of the cross


I=createMask_black(I);
I=imcrop(I,ROI);


BW=~I;
%  figure,imshow(BW);
BW=imclearborder(BW); % entfernen der Randobjekte
%   imshow(BW);

%% Linien detektieren
[H,T,R] = hough(BW);




P  = houghpeaks(H,10,'threshold',ceil(0.3*max(H(:))));
x = T(P(:,2)); y = R(P(:,1));
% plot(x,y,'s','color','white');

lines = houghlines(BW,T,R,P,'FillGap',1,'MinLength',1);
%   figure, imshow(I), hold on
max_len = 0;

    function []= plotlines()
for k = 1:length(lines)
   xy = [lines(k).point1; lines(k).point2];
   plot(xy(:,1),xy(:,2),'LineWidth',2,'Color','green');

   % Plot beginnings and ends of lines
   plot(xy(1,1),xy(1,2),'x','LineWidth',2,'Color','yellow');
   plot(xy(2,1),xy(2,2),'x','LineWidth',2,'Color','red');

   % Determine the endpoints of the longest line segment
   len = norm(lines(k).point1 - lines(k).point2);
   if ( len > max_len)
      max_len = len;
      
   end
end
    end


%calculate the cross coords
upper_limit=100000;
lower_limit=0;
x_min=upper_limit;% large pixelcount 
y_min=upper_limit;
x_max=lower_limit;% small 
y_max=lower_limit;

x_sum=0;
y_sum=0;
%% Bestimmen des average Punkts um ROI mithilfe eines Radius zu bestimmen
for k= 1:length(lines)
    x_sum=x_sum+lines(k).point1(1)+lines(k).point2(1);
    y_sum=y_sum+lines(k).point1(2)+lines(k).point2(2);
    
    
end    

x_avg=x_sum/(2*length(lines));
y_avg=y_sum/(2*length(lines));

radius=300;
errorNo=0;
%% Bestimme den Mittelpunkt des Kreuzes. hier anhand der Mitte der Extremwerte
%  avg der Werte auch denkbar.(war teilweise weniger robust)
function[]= calculateIntersection()
    errorNo=errorNo+1;
for k= 1:length(lines)
  
    %% eingrenzen des ROI fuer Linien
    if (abs(lines(k).point1(1)-x_avg)<radius && abs(lines(k).point1(2)-y_avg)<radius)
        %% betrachte Horizontale Linien
        if (abs(lines(k).theta)>80)
            if(lines(k).point1(2)<y_min)
                    y_min=lines(k).point1(2);
            end
                
                if(lines(k).point1(2)>y_max)
                    y_max=lines(k).point1(2);
                end
        else
        
        %% Betrachte Vertikale
            if(abs(lines(k).theta)<10)
                if(lines(k).point1(1)<x_min)
                    x_min=lines(k).point1(1);
                end
                
                if(lines(k).point1(1)>x_max)
                    x_max=lines(k).point1(1);
                end
                

            end
        end
    end


end
if(x_max==lower_limit || x_min==upper_limit)
%     error('ERROR:  Kreuz nicht im ROI oder zu klein / ungenau')   


    figure; imshow(I);
    
    if(errorNo==2)
    title('ErrorLvl 2: Zoome in das Kreuz! und drücke Enter');

        zoom on;
        waitfor(gcf, 'CurrentCharacter', char(13))
        zoom reset
        zoom off
            title('wähle das Kreuz aus mit linker Maustaste!');

        [x,y] = ginput(1)       
    else
        title('ERROR Lvl1 :X Koord Kreuz nicht im ROI wähle Kreuz wähle ROI (Grob) manuell');
        [x_avg,y_avg]=ginput(1);
        calculateIntersection();
    end
    
    
    
end

if(y_max==lower_limit || y_min==upper_limit)
%     error('ERROR:  Kreuz nicht im ROI oder zu klein / ungenau')
   figure; imshow(I);
    
    if(errorNo==2)
    title('ErrorLvl 2: Zoome in das Kreuz! und drücke Enter');

        zoom on;
        waitfor(gcf, 'CurrentCharacter', char(13))
        zoom reset
        zoom off
            title('wähle das Kreuz aus mit linker Maustaste!');

        [x,y] = ginput(1)       
    else
        title('ERROR Lvl1 :Y Koord Kreuz nicht im ROI wähle Kreuz wähle ROI (Grob) manuell');
        [x_avg,y_avg]=ginput(1);
        calculateIntersection();
    end
end

%% Wenn Fehler auftritt soll manuell ein ROI ausgewählt werden können.
if(errorNo<2)
   y=(y_max+y_min)/2;
   x=(x_max+x_min)/2;
end
end

 calculateIntersection();
 plot(x,y,'x','LineWidth',4,'Color','yellow');

close gcf
%% Berechne die Koordinaten im gesamten (urspruenglichen) Bild
 x=ROI(1)+x;
y=ROI(2)+y;
% Swap x und y. Grund nicht ganz klar. funktioniert so.
% c=x;
% x=y;
% y=c;


end