t1_5

function [mImage, oImage, eImage] = sobelFunc(image, threshold)
    sxMat = [-1 0 1; -2 0 2; -1 0 1];
    syMat = [1 2 1; 0 0 0; -1 -2 -1];

    convX = conv2(image, sxMat, 'same');
    convY = conv2(image, syMat, 'same');

    mImage = abs(convX) + abs(convY);
    oImage = atan2(convY, convX);

    eImage = zeros(size(mImage, 1), size(mImage, 2));
    for i = 1 : size(mImage, 1)
        for j = 1 : size(mImage, 2)
            if mImage(i, j) > threshold
                eImage(i, j) = 255;
            end
        end
    end
end