1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39
| function [Gx, Gy, norme, normediv] = gradient( Im )
%GRADIENT Summary of this function goes here
% Detailed explanation goes here
Im = double(Im);
s = size(Im);
Gx = zeros(size(Im));
Gy = zeros(size(Im));
norme= zeros(size(Im));
normediv= zeros(size(Im));
divIm = zeros(size(Im));
for i=1 : s(1) -1
for j=1 : s(2)
Gx(i,j)= Im(i+1, j) - Im(i,j);
end
end
for i=1 : s(1)
for j=1 : s(2) -1
Gy(i,j)= Im(i, j+1) - Im(i,j);
end
end
for i=1 : s(1)
for j=1 : s(2)
norme(i,j)= sqrt(Gx(i,j)^2 + Gy(i,j)^2);
end
end
for i=1 : s(1)
for j=1 : s(2)
divIm(i,j) = Gx(i,j) + Gy(i,j);
normediv(i,j) = abs(divIm(i,j));
end
end
end |
Partager