image = imread("lenna.png");

image = cast(image, 'double');

% Gray
gray_image = cvt_rgb2gray(image);
gray_image = cast(gray_image,'double');

% Filters
sobel_e = [
    [-1,  0,  1]
    [-2,  0,  2]
    [-1,  0,  1]
    ];

sobel_w = [
    [ 1,  0, -1]
    [ 2,  0, -2]
    [ 1,  0, -1]
    ];

% Filtering
filtered_e = apply_kernel(gray_image, sobel_e);
filtered_w = apply_kernel(gray_image, sobel_w);

filtered_eabs = abs(filtered_e);
filtered_e(:,:,:) = filtered_e(:,:,:) - filtered_eabs(:,:,:);

filtered_wabs = abs(filtered_w);
filtered_w(:,:,:) = filtered_w(:,:,:) - filtered_wabs(:,:,:);

normalized_e = normalize_image(filtered_e);
normalized_w = normalize_image(filtered_w);

combined = combine_images(filtered_e, filtered_w, 1);
combined = normalize_image(combined);










#####



function [outputArg1] = combine_images(image, image2, value)
%BRIGHTNESS Summary of this function goes here
%   Detailed explanation goes here

[row, col, ~] = size(image);

target = image;

for i=1:row
    for j=1:col

        target(i,j) = (image(i,j) + value*image2(i,j));

    end
end
% target = normalizeimage(target);

outputArg1 = target;
end