generate_noise_pic.m

an implementation of Notch_Filter in matlab for image processing

function [noise_pic]=generate_Periodic_Noise(pic)

pic=double(pic);

A=10;
[M,N]=size(pic);
u0=M/2;
v0=0;

noise=sinusoidal_noise(A,u0,v0,M,N);
noise=double(noise);

noise_pic=pic+noise;

noise_pic=255*(noise_pic-min(min(noise_pic)))/(max(max(noise_pic))-min(min(noise_pic)));
noise_pic=uint8(noise_pic);
imwrite(noise_pic,'noise_pic.bmp');




% Calculate noise and Pic_noise spectrum
%for x=1:size(noise,1)
%    for y=1:size(noise,2)
%        noise(x,y)=(-1)^(x+y)*noise(x,y);
%        noise_pic_temp(x,y)=(-1)^(x+y)*double(noise_pic(x,y));
        %    end
    %end

%noise_spectrum=log10(abs(fft2(noise))+1);
%noise_spectrum=255*(noise_spectrum-min(min(noise_spectrum)))/(max(max(noise_spectrum))-min(min(noise_spectrum)));

%noise_pic_spectrum=log10(abs(fft2(noise_pic_temp))+1);
%noise_pic_spectrum=255*(noise_pic_spectrum-min(min(noise_pic_spectrum)))/(max(max(noise_pic_spectrum))-min(min(noise_pic_spectrum)));

%noise=255*(noise-min(min(noise)))/(max(max(noise))-min(min(noise)));

%figure(1);
%subplot(2,2,1);imshow(uint8(noise));
%subplot(2,2,2);imshow(uint8(noise_pic));
%subplot(2,2,3);imshow(uint8(noise_spectrum));
%subplot(2,2,4);imshow(uint8(noise_pic_spectrum));