clc;
clear;
close all;
img = imread('crop/1.png');
img = im2double(img);
[~,~,imxy] = my_gradient(img);
s1 = mean(imxy(:));
disp(s1);
function [imx,imy,imnorm] = my_gradient(im)
imx = [im(:,end,:) - im(:, 1,:),-diff(im,1,2)];
imy = [im(end,:,:) - im(1, :,:);-diff(im,1,1)];
imnorm = sqrt(imx.^2 + imy.^2);
end