From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 125/CH6/EX6.5/Fig6_5.sce                           |  37 +++++++++++++++++++++
 .../Fig6_7_Degraded_NoiseImage_Restoration.JPG     | Bin 0 -> 47700 bytes
 2 files changed, 37 insertions(+)
 create mode 100755 125/CH6/EX6.5/Fig6_5.sce
 create mode 100755 125/CH6/EX6.5/Fig6_7_Degraded_NoiseImage_Restoration.JPG

(limited to '125/CH6/EX6.5')

diff --git a/125/CH6/EX6.5/Fig6_5.sce b/125/CH6/EX6.5/Fig6_5.sce
new file mode 100755
index 000000000..c82622ead
--- /dev/null
+++ b/125/CH6/EX6.5/Fig6_5.sce
@@ -0,0 +1,37 @@
+//Caption:Scilab code performs inverse filtering 
+//Degrade the image by means of a known blur
+//Apply inverse filter to the blurred image and see the restored image
+//Fig6.5
+//page 330
+clc;
+close;
+x =imread('E:\DIP_JAYARAMAN\Chapter6\flower2.jpg');
+x=double(rgb2gray(x));
+[M N]=size(x);
+h = zeros(M,N);
+for i = 1:11
+    for j = 1:11
+        h(i,j) = 1/121;
+    end
+end
+sigma = sqrt(4*10^(-7));
+freqx = fft2d(x); //Fourier transform of input image
+freqh = fft2d(h);//Fourier transform of degradation
+y = real(ifft2d(freqh.*freqx));
+freqy = fft2d(y);
+powfreqx = freqx.^2/(M*N);
+alpha = 0.5; //Indicates inverse filter 
+freqg = ((freqh.')').*abs(powfreqx)./(abs(freqh.^2).*abs(powfreqx)+alpha*sigma^2);
+Resfreqx = freqg.*freqy;
+Resa = real(ifft2d(Resfreqx));
+x = uint8(x);
+y = uint8(y);
+Resa = uint8(Resa)
+ShowImage(x,'Original Image')
+title('Original Image')
+figure
+ShowImage(y,'Degraded Image')
+title('Degraded Image')
+figure
+ShowImage(Resa,'Restored Image')
+title('Restored Image')
\ No newline at end of file
diff --git a/125/CH6/EX6.5/Fig6_7_Degraded_NoiseImage_Restoration.JPG b/125/CH6/EX6.5/Fig6_7_Degraded_NoiseImage_Restoration.JPG
new file mode 100755
index 000000000..8b8ff8a19
Binary files /dev/null and b/125/CH6/EX6.5/Fig6_7_Degraded_NoiseImage_Restoration.JPG differ
-- 
cgit