function [ S_pho_x,S_pho_y ] = S_pho_xy( M_Ascii,von_zeile,bis_zeile,von_spalte,bis_spalte )
%Bestimmung des geometrischen Schwerpunkts in x- und y-Richtung einer Messregion
%   Detailed explanation goes here

%Flächeninhalt eines Pixels
A = 1;
laenge_A = sqrt(A);  %Seitenlänge des Pixels

%Setzen der Variablen wegen selbst-inkrementieren
S_Zaehler_x = 0;
S_Nenner_x = 0;
S_Zaehler_y = 0;
S_Nenner_y = 0;

%3.Schleife die jeden Punkt der input-Matrix anwählt
for i=von_zeile:1:bis_zeile      %Schleife für Anfang bis Ende der einzelnen Zeilen 
    for j=von_spalte:1:bis_spalte     %Durchläuft alle Spalten der jeweiligen Zeile
        
        %3.1 Abgleich ob der Wert der einzelnen Elemente >= q_005
        %entspricht BZW größer Schwellenwert z.B. 1 cd/m2
        if isnan(M_Ascii(i,j))
            continue
        else
        %Berechnet den Zaehler und Nenner für den Schwepunkt in x- und
        %y-Richtung
            %Anpassug der Berechnung an die Matrixhandhabung in Matlab (1.
            %Zeile und 1. Spalte gesondert berechnen)
            %Spalte
            if j == 1
                delta_x = 1/2 * laenge_A;
                S_Zaehler_x = S_Zaehler_x + (delta_x * A * M_Ascii(i,j));
                S_Nenner_x = S_Nenner_x + (A * M_Ascii(i,j));
            else
                delta_x = 1/2 * laenge_A + (j-1) * laenge_A;
                S_Zaehler_x = S_Zaehler_x + (delta_x * A * M_Ascii(i,j));
                S_Nenner_x = S_Nenner_x + (A * M_Ascii(i,j));
            end
            %Zeile
            if i == 1
                delta_y = 1/2 * laenge_A;
                S_Zaehler_y = S_Zaehler_y + (delta_y * A * M_Ascii(i,j));
                S_Nenner_y = S_Nenner_y + (A * M_Ascii(i,j));
            else
                delta_y = 1/2 * laenge_A + (i-1) * laenge_A;
                S_Zaehler_y = S_Zaehler_y + (delta_y * A * M_Ascii(i,j));
                S_Nenner_y = S_Nenner_y + (A * M_Ascii(i,j));
            end
        end
            
        
    end
end

%Berechnung des Schwerpunktes in x- und y-Richtung
S_pho_x = S_Zaehler_x / S_Nenner_x;
S_pho_y = S_Zaehler_y / S_Nenner_y;

end