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

---
 3176/CH5/EX5.2/Ex5_2.sce | 214 +++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 214 insertions(+)
 create mode 100644 3176/CH5/EX5.2/Ex5_2.sce

(limited to '3176/CH5/EX5.2/Ex5_2.sce')

diff --git a/3176/CH5/EX5.2/Ex5_2.sce b/3176/CH5/EX5.2/Ex5_2.sce
new file mode 100644
index 000000000..607616bb8
--- /dev/null
+++ b/3176/CH5/EX5.2/Ex5_2.sce
@@ -0,0 +1,214 @@
+//Ex5_2
+// Illustration of Mean Filters
+//To impliment the Following Mean Restoration filter 
+//                  (I)Arithmetic  (II)Geometric  (III)Harmonic (IV)Contra Harmonic
+
+// 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).
+
+
+function [f]=arithmetic_mean(v,m,n)
+        w=fspecial('average',m); 
+        f=imfilter(v,w); 
+endfunction
+
+function [f]=geometric_mean1(g,m,n);//gmean1() is used to filter an image using Geometric mean filter 
+    size1=m;
+    q=m*n;
+    g=double(g);
+    [nr,nc]=size(g);
+    temp=zeros(nr+2*floor(size1/2),nc+2*floor(size1/2));
+    temp(ceil(size1/2):nr+ceil(size1/2)-1,ceil(size1/2):nc+ceil(size1/2)-1)=g(1:$,1:$)
+    temp=temp+1;
+    for i=ceil(size1/2):nr+ceil(size1/2)-1
+        for j=ceil(size1/2):nc+ceil(size1/2)-1
+            t=temp(i-floor(size1/2):1:i+floor(size1/2),j-floor(size1/2):1:j+floor(size1/2)) ;
+            temp2(i,j)=prod(t);                
+        end
+    end
+    temp3=temp2.^(1/q);
+    nn=temp3(ceil(size1/2):nr+ceil(size1/2)-1,ceil(size1/2):nc+ceil(size1/2)-1)
+    f1=nn-1;
+    f=mat2gray(f1)
+endfunction
+
+function [f]=geometric_mean2(g,m,n);//gmean2() is used to filter an image using Geometric mean filter 
+    size1=m;
+    q=m*n;
+    [nr,nc]=size(g);
+    temp=zeros(nr+2*floor(size1/2),nc+2*floor(size1/2));
+    temp(ceil(size1/2):nr+ceil(size1/2)-1,ceil(size1/2):nc+ceil(size1/2)-1)=g(1:$,1:$)
+    for i=ceil(size1/2):nr+ceil(size1/2)-1
+        for j=ceil(size1/2):nc+ceil(size1/2)-1
+            t=temp(i-floor(size1/2):1:i+floor(size1/2),j-floor(size1/2):1:j+floor(size1/2)) ;
+            temp2(i,j)=geomean(t);                   
+        end
+    end
+    nn=temp2(ceil(size1/2):nr+ceil(size1/2)-1,ceil(size1/2):nc+ceil(size1/2)-1)
+    f=mat2gray(nn)
+endfunction
+
+function [f]=Harmonic_mean(g,m,n) //harmean1() is used to filter an image using Harmonic mean filter.
+    size1=m;
+    d=m*n;
+    g=double(g);
+    [nr,nc]=size(g);
+    temp=zeros(nr+2*floor(size1/2),nc+2*floor(size1/2));
+    temp(ceil(size1/2):nr+ceil(size1/2)-1,ceil(size1/2):nc+ceil(size1/2)-1)=g(1:$,1:$);
+
+    for i=ceil(size1/2):nr+ceil(size1/2)-1
+        for j=ceil(size1/2):nc+ceil(size1/2)-1
+            t=temp(i-floor(size1/2):1:i+floor(size1/2),j-floor(size1/2):1:j+floor(size1/2)) ;
+            t1=ones(m,n)./(t+%eps);
+            t2=sum(t1);
+            temp2(i,j)=d/t2;            
+        end
+    end
+    nn=temp2(ceil(size1/2):nr+ceil(size1/2)-1,ceil(size1/2):nc+ceil(size1/2)-1);
+    f=mat2gray(nn);
+endfunction
+
+function [f]=Contra_Harmonic_mean(g,m,n,Q) //charmean1() is use to filter an image using Contra Harmonic mean filter
+    size1=m;
+    d=m*n;
+    g=double(g);
+    [nr,nc]=size(g);
+    temp=zeros(nr+2*floor(size1/2),nc+2*floor(size1/2));
+    temp(ceil(size1/2):nr+ceil(size1/2)-1,ceil(size1/2):nc+ceil(size1/2)-1)=g(1:$,1:$)
+    disp(Q)
+    for i=ceil(size1/2):nr+ceil(size1/2)-1
+        for j=ceil(size1/2):nc+ceil(size1/2)-1
+            t=temp(i-floor(size1/2):1:i+floor(size1/2),j-floor(size1/2):1:j+floor(size1/2)) ;
+            d1=(t+%eps).^Q;
+            n1=(t+%eps).^(Q+1);
+            d2=sum(d1);
+            n2=sum(n1);
+            temp2(i,j)=n2/(d2);
+        end
+    end
+    nn=temp2(ceil(size1/2):nr+ceil(size1/2)-1,ceil(size1/2):nc+ceil(size1/2)-1)
+    f=nn;
+endfunction
+
+/////////////////////////////////////    Main Programm    ////////////////////
+
+gray=imread("Ex5_2.tif");
+//gray=rgb2gray(a);
+//gray=im2double(gray);
+figure,ShowImage(gray,'Gray Image');
+title('Original Image');
+[M,N]=size(gray);
+
+/////////////////////////////////////    Arithmetical Mean Filter    ////////////////////
+v=imnoise(gray,'gaussian',0,0.02);
+figure,ShowImage(v,'Noisy Image');
+title('Image with Gaussian Noise');
+m=3;n=3;
+[f]=arithmetic_mean(v,m,n);
+figure,ShowImage(f,'Recovered Image');
+title('Recovered Image with Arithmetical Mean Filter');
+
+/////////////////////////////////////    Geometric Mean Filter    ////////////////////
+v=imnoise(gray,'gaussian',0,0.02);
+figure,ShowImage(v,'Noisy Image');
+title('Image with Gaussian Noise');
+m=3;n=3;
+[f]=geometric_mean1(v,m,n);
+figure,ShowImage(f,'Recovered Image');
+title('Recovered Image with Geometric Mean Filter');
+
+
+///////////////////////////////////////    Geometric Mean Filter    ////////////////////
+//v=imnoise(gray,'gaussian',0,0.02);
+//figure,ShowImage(v,'Noisy Image');
+//title('Image with Gaussian Noise');
+//m=3;n=3;
+//[f]=geometric_mean2(v,m,n);
+//figure,ShowImage(f,'Recovered Image');
+//title('Recovered Image with Geometric Mean Filter');
+
+
+///////////////////////////////////////    Harmonic Mean Filter    ////////////////////
+//temp(1:M,1:N)=0.5;
+//r3=imnoise(temp,'salt & pepper',0.1);   // Generate salt & pepper Noise 
+//gray_noise_salt=gray;                  // Add salt Noise Only
+//[r c]=find(r3==1);
+//        for i=1:length(r)
+//                gray_noise_salt(r(i),c(i)) = 255; 
+//        end
+//figure,ShowImage(gray_noise_salt,'Noisy Image');
+//title('Image with Salt Noise');
+//m=3;n=3;
+//[f]=Harmonic_mean(gray_noise_salt,m,n);
+//figure,ShowImage(f,'Recovered Image');
+//title('Recovered Image with Harmonic Mean Filter');
+//
+////////////////////////////////  Contra_Harmonic Mean Filter (Pepper)   ////////////////////
+temp(1:M,1:N)=0.5;
+r3=imnoise(temp,'salt & pepper',0.05);    //Generate salt & pepper Noise 
+gray_noise_pepper=gray;                   //Add pepper Noise Only
+[r c]=find(r3==0);                        //Find pepper Noise Only
+        for i=1:length(r)
+                gray_noise_pepper(r(i),c(i)) = 0; 
+        end
+figure,ShowImage(gray_noise_pepper,'Noisy Image');
+title('Image with pepper Noise');
+m=3;n=3;Q=1.5;
+[f]=Contra_Harmonic_mean(gray_noise_pepper,m,n,Q);
+figure,ShowImage(f,'Recovered Image');
+title('Recovered Image with Contra Harmonic Mean Filter[ Q=1.5 ]');
+
+///////////////////////////////////    Contra_Harmonic Mean Filter (Salt)   ////////////////////
+temp(1:M,1:N)=0.5;
+r3=imnoise(temp,'salt & pepper',0.1);    //Generate salt & pepper Noise 
+gray_noise_salt=gray;                   //Add salt Noise Only
+[r c]=find(r3==1);
+        for i=1:length(r)
+                gray_noise_salt(r(i),c(i)) = 255; 
+        end
+figure,ShowImage(gray_noise_salt,'Noisy Image');
+title('Image with Salt Noise');
+m=3;n=3;Q=-1.5;
+[f]=Contra_Harmonic_mean(gray_noise_salt,m,n,Q);
+figure,ShowImage(f,'Recovered Image');
+title('Recovered Image with Contra Harmonic Mean Filter[ Q=-1.5 ]');
+
+
+//////////////////////////////  Contra_Harmonic Mean Filter (Pepper)   ////////////////////
+temp(1:M,1:N)=0.5;
+r3=imnoise(temp,'salt & pepper',0.05);    //Generate salt & pepper Noise 
+gray_noise_pepper=gray;                  // Add pepper Noise Only
+[r c]=find(r3==0);                        //Find pepper Noise Only
+        for i=1:length(r)
+                gray_noise_pepper(r(i),c(i)) = 0; 
+        end
+figure,ShowImage(gray_noise_pepper,'Noisy Image');
+title('Image with pepper Noise');
+m=3;n=3;Q=-1.5;
+[f]=Contra_Harmonic_mean(gray_noise_pepper,m,n,Q);
+figure,ShowImage(f,'Recovered Image');
+title('Recovered Image with Contra Harmonic Mean Filter[ Q=-1.5 ]');
+
+///////////////////////////////////    Contra_Harmonic Mean Filter (Salt)   ////////////////////
+temp(1:M,1:N)=0.5;
+r3=imnoise(temp,'salt & pepper',0.1);    //Generate salt & pepper Noise 
+gray_noise_salt=gray;                   //Add salt Noise Only
+[r c]=find(r3==1);
+        for i=1:length(r)
+                gray_noise_salt(r(i),c(i)) = 255; 
+        end
+figure,ShowImage(gray_noise_salt,'Noisy Image');
+title('Image with Salt Noise');
+m=3;n=3;Q=1.5;
+[f]=Contra_Harmonic_mean(gray_noise_salt,m,n,Q);
+figure,ShowImage(f,'Recovered Image');
+title('Recovered Image with Contra Harmonic Mean Filter[ Q=1.5 ]');
-- 
cgit