Image_processing/Task4_code.m at master · murtuzaal2/Image_processing · GitHub

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
40
41
42
43
44
45
46
47
48
49
50
51
52
clear;
u0 = imread('test.png');
u0 = rgb2gray(u0);
u0 = im2double(u0);
u0 = imnoise(u0,'gaussian',0.1);
u = u0;

u1 = u;

lambda = 0.01;

for t = 1:100;

for i = 2:299;
for j = 2:299;

ux(i,j) = (u(i+1,j) - u(i-1,j))^2/2;
uy(i,j) = (u(i,j+1) - u(i,j-1))^2/2;


g(i,j) = 1/(1+((ux(i,j)+uy(i,j))/(lambda ^2)));
end
end

for i = 2:298;
for j = 2:298;

Gx1(i,j) = (g(i,j)+g(i+1,j))/2;
Gx2(i,j) = (g(i,j)+g(i-1,j))/2;
Gy1(i,j) = (g(i,j)+g(i,j+1))/2;
Gy2(i,j) = (g(i,j)+g(i,j-1))/2;

end
end

for i = 2:298;
for j = 2:298;

u1(i,j) = u(i,j) + (Gx1(i,j)*(u(i+1,j)-u(i,j))) - (Gx2(i,j)*(u(i,j)-u(i-1,j))) + (Gy1(i,j)*(u(i,j+1)-u(i,j)))-(Gy2(i,j)*(u(i,j)-u(i,j-1)));

end
u = u1;

end

end

figure
subplot(1,2,1),
imshow(u0);
subplot(1,2,2),
imshow(u1);