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
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
|
//Ex10_22
// Document Thresholding Using Moving Averages
// Version : Scilab 5.4.1
// Operating System : Window-xp, Window-7
//Toolbox: Image Processing Design 8.3.1-1
//Toolbox: SIVP 0.5.3.1-2
//Reference book name : Digital Image Processing
//book author: Rafael C. Gonzalez and Richard E. Woods
clc;
close;
clear;
xdel(winsid())//to close all currently open figure(s).
a=imread("Ex10_22.tif");
a1=im2double(a);
figure,ShowImage(a1,'Gray Image');
title('Original Image','color','blue','fontsize',4);
[M,N]=size(a);
Threshold = CalculateOtsuThreshold(a1);
Thresh_Image=im2bw(a1,Threshold);
figure,ShowImage(Thresh_Image,'Binary Image');
title('Thresholded Image with Otsu Method','color','blue','fontsize',4);
mask=zeros(1,20);
array=[];
for i=1:M
if(pmodulo(i,2)==0)
array=[array mtlb_fliplr(a1(i,:))];
else
array=[array a1(i,:)];
end
end
disp('first');
for i=1:length(array)
for j=1:length(mask)
if(j<length(mask)) then
mask(j)=mask(j+1);
else
mask(j)=array(i);
end
end
avg(1,i)=sum(mask)/length(mask);
end
disp('Second');
len=1;
for i=1:M
if(pmodulo(i,2)==0)
b(i,:)=avg(len:len+N-1);
len=len+N;
else
b(i,:)=avg(len:len+N-1);
len=len+N;
end
end
disp('Last');
b=0.5*b;
for i=1:M
for j=1:N
if(b(i,j)>a1(i,j)) then
c(i,j)=0;
else
c(i,j)=1;
end
end
end
figure,ShowImage(c,'Binary Image');
title('Local Thresholding Using Moving Average','color','blue','fontsize',4);
|
