diff options
Diffstat (limited to '3176/CH5/EX5.10/Ex5_10.sce')
-rw-r--r-- | 3176/CH5/EX5.10/Ex5_10.sce | 44 |
1 files changed, 44 insertions, 0 deletions
diff --git a/3176/CH5/EX5.10/Ex5_10.sce b/3176/CH5/EX5.10/Ex5_10.sce new file mode 100644 index 000000000..dddcbe872 --- /dev/null +++ b/3176/CH5/EX5.10/Ex5_10.sce @@ -0,0 +1,44 @@ +//Ex5_10
+//Image Bluring Due to Motion
+// 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).
+
+gray=imread("Ex5_10.png");
+gray=im2double(rgb2gray(gray));
+//gray=im2double(imresize(a,0.5));
+figure,ShowImage(gray,'Gray Image');
+title('Original Image','color','blue','fontsize',4);
+[M,N]=size(gray);
+
+h=fft2(gray);//fft2() is used to find 2-Dimensional Fast Fourier Transform of an matrix
+i=log(1+abs(h));
+in=fftshift(i);//fftshift() is used to rearrange the fft output, moving the zero frequency to the center of the spectrum.
+inm=mat2gray(in)
+
+a=0.1;b=0.1;T=1; // Motion and Exposure Value
+for u=1:M
+ for v=1:N
+ H(u,v)=(T/(%pi*(u*a+v*b)))*(sin(%pi*(u*a+v*b)))*exp(-%i*%pi*(u*a+v*b)); //Motion Blure Function
+ end
+end
+
+n=h.*H;//Multiply the Original Spectrum with the Degradation Function.
+Image_filter=abs(ifft(n));
+Image_filter=mat2gray(Image_filter)
+figure,ShowImage(Image_filter,'Filtered Image');
+title('Motion Blure Image','color','blue','fontsize',4);
+
+
+
+
+
+
|