%%%%%%%%%%%%%%%
%%请不要改动
warning('off','all');
graphics_toolkit('gnuplot')
pkg load image
%%%%%%%%%%%%%%%%%%
figure(1)
%%%%%%%%% Begin %%%%%%%%%
F=imread('./img/moon_bw.tif');
N=ndims(F);
if N>2
F=rgb2gray(F);
end;
F=im2double(F);
Fn=imnoise(F,'salt & pepper',0.02);
figure(1)
h=fspecial('average',3);
G=imfilter(Fn,h);
subplot(1,2,1);
imshow(Fn);
subplot(1,2,2);
imshow(G);
%%%%%%%%% End %%%%%%%%%
sa=pwd;
print(1,'-djpeg','./pictures/step1/picture1.jpg');
run('./task1/test1.m');
system('python3 ./task1/test1py.py');
第2关:中值滤波
%%%%%%%%%%%%%%%
%%请不要改动
warning('off','all');
graphics_toolkit('gnuplot')
pkg load image
%%%%%%%%%%%%%%%%%%
figure(1)
%%%%%%%%% Begin %%%%%%%%%
I=imread('./img/pcbxray_gray.tif');
N=ndims(I);
if N>2
I=rgb2gray(I);
end;
F=im2double(I);
Fn=imnoise(F,'salt & pepper',0.1);
figure(1)
subplot(1,2,1);
imshow(Fn);
subplot(1,2,2);
G=medfilt2(Fn,[3,3]);
imshow(G);
%%%%%%%%% End %%%%%%%%%
sa=pwd;
print(1,'-djpeg','./pictures/step2/picture1.jpg');
run('./task2/test1.m');
system('python3 ./task2/test1py.py');
第3关:锐化滤波——roberts算子
%%%%%%%%%%%%%%%
%%请不要改动
warning('off','all');
graphics_toolkit('gnuplot')
pkg load image
%%%%%%%%%%%%%%%%%%
figure(1)
%%%%%%%%% Begin %%%%%%%%%
I=imread('img/moon_bw.tif');
N=ndims(I);
if N>2
I=rgb2gray(I);
end;
I=im2double(I);
h1=[-1 0;0 1];
h2=[0 -1;1 0];
I1=imfilter(I,h1);
I2=imfilter(I,h2);
If=abs(I1)+abs(I2);
G=I+If;
figure(1)
subplot(1,3,1);
imshow(I);
subplot(1,3,2);
imshow(If);
subplot(1,3,3);
imshow(G);
%%%%%%%%% End %%%%%%%%%
sa=pwd;
print(1,'-djpeg','./pictures/step3/picture1.jpg');
run('./task3/test1.m');
system('python3 ./task3/test1py.py');
第4关:锐化滤波——laplacian算子
%%%%%%%%%%%%%%%
%%请不要改动
warning('off','all');
graphics_toolkit('gnuplot')
pkg load image
%%%%%%%%%%%%%%%%%%
figure(1)
%%%%%%%%% Begin %%%%%%%%%
F=imread('./img/moon_bw.tif');
N=ndims(F);
if N>2
F=rgb2gray(F);
end;
F=im2double(F);
h=fspecial('laplacian',0);
G=imfilter(F,h);
G1=F-G;
figure(1)
subplot(1,3,1);
imshow(F,[]);
subplot(1,3,2);
imshow(G,[]);
subplot(1,3,3);
imshow(G1,[]);
%%%%%%%%% End %%%%%%%%%
sa=pwd;
print(1,'-djpeg','./pictures/step4/picture1.jpg');
run('./task4/test1.m');
system('python3 ./task4/test1py.py');