fig6_8_ebma.m

these are the Mtlab functions to plot the figures in the book Communication system (by Haykin)

%for performing EBMA between two images
%main file

%use the following if the program will be used as a function
%function[]=ebma(FileName1, FileName2,dx,dy, MVname, PImgName);
%frame1=fread(fopen(FileName1),[dx,dy]);
%frame2=fread(fopen(FileName2),[dx,dy]);

dx=352;dy=288;
frame1=fread(fopen('foreman69.Y'),[dx,dy]);
frame2=fread(fopen('foreman72.Y'),[dx,dy]);

img1=zeros(dy,dx);img2=zeros(dy,dx);
img1=frame1'; img2=frame2';
figure(1);
image(img1), colormap(gray(256)); axis image;
figure(2);
image(img2), colormap(gray(256)); axis image;

bsize=16; srange=16;
pimg=img1;
bdy=ceil(dy/bsize); bdx=ceil(dx/bsize);
mvy=zeros(bdy,bdx); mvx=zeros(bdy,bdx);

[mvy,mvx,pimg]=matching_half(img1, img2, dy, dx, bsize, srange, srange, mvy, mvx, pimg);
figure(3);
image(pimg), colormap(gray(256)); axis image; axis off;
print -deps2 foreman_EBMA_pimg.eps

figure(4);
imagesc(abs(img2-pimg)), colormap(gray(256)); axis image;
figure(5);
plot_MV_function(dy,dx,bsize,mvy,mvx);set(gca,'XTick',[],'YTick',[]);
print -deps2 foreman_EBMA_MV.eps

mse=sum(sum((img2-pimg).*(img2-pimg)))/dx/dy;
psnr=10*log10(255*255/mse);
%29.85