function [ S_geo_x,S_geo_y ] = S_geo_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 quadratischen 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);
                S_Nenner_x = S_Nenner_x + A;
            else
                delta_x = 1/2 * laenge_A + (j-1) * laenge_A;
                S_Zaehler_x = S_Zaehler_x + (delta_x * A);
                S_Nenner_x = S_Nenner_x + A;
            end
            %Zeile
            if i == 1
                delta_y = 1/2 * laenge_A;
                S_Zaehler_y = S_Zaehler_y + (delta_y * A);
                S_Nenner_y = S_Nenner_y + A;
            else
                delta_y = 1/2 * laenge_A + (i-1) * laenge_A;
                S_Zaehler_y = S_Zaehler_y + (delta_y * A);
                S_Nenner_y = S_Nenner_y + A;
            end
        end
        
        
    end
end

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

end