%   BHPF
%   2/09
%   HE

close all; clear;
im=imread('C:\Hector\EE_552\Images\scene5\VIS_0109.bmp');
[M,N,c]=size(im);

[x,y]=meshgrid(1:M,1:N);
K1=2^8-1; K2=0;
in=(K1-K2)*exp(-((x-M/2).^2+(y-N/2).^2)./(2*50^2))+K2; in=in';
imR=double(im(:,:,1)).*in;
figure; imagesc(imR); colormap(gray);
save 'Test1_data1' imR
figure; imagesc(imR./in); colormap(gray);
figure; plot(in(M/2,:));