matlab实现了图像通信中的全搜索算法与三步搜索算法

用matlab实现了图像通信中的全搜索算法与三步搜索算法这两种运动估值算法，给出了运动矢量图，对两种方法的计算复杂度和搜索性能进行了客观的比较。

clear all;

I1=imread（‘claire1.bmp‘）;  %read the first frame
I2=imread（‘claire2.bmp‘）;  %read the second frame


dm=7;
%给图像扩边，每个边都扩dm大小
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%         |-----|-------------------------------|------|
%         |  B  |              C                |  D   |
%         |-----|-------------------------------|------|
%         |     |                               |      |
%         |     |                               |      |
%         |  E  |               A               |  F   |   ----> 图像
%         |     |                               |      |
%         |     |                               |      |
%         |     |                               |      |
%         |-----|-------------------------------|------|
%         |  G  |              H                |  I   |
%         |-----|-------------------------------|------|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

I1=double（I1）;
I2=double（I2）;
[rownum colnum] = size（I1）;


%为扩的边赋值  请自己写程序代码
% v=0;
% II=[v*ones（rownumdm） I1 v*ones（rownumdm）];
% II=[v*ones（dmcolnum+2*dm）;II;v*ones（dmcolnum+2*dm）];

II=zeros（rownum+2*dmcolnum+2*dm）;
II（dm+1:dm+rownumdm+1:dm+colnum）=I1;
for i=1:dm
    II（idm+1:dm+colnum）=II（dm+1dm+1:dm+colnum）;
    II（rownum+dm+idm+1:dm+colnum）=II（dm+rownumdm+1:dm+colnum）;
end
for j=1:dm
    II（1:rownum+2*dmj）=II（1:rownum+2*dmdm+1）;
    II（1:rownum+2*dmcolnum+dm+j）=II（1:rownum+2*dmdm+colnum）;
end

%下面进行全搜索算法
 blocksize=16;
 rowblocks =rownum/blocksize;
 colblocks =colnum/blocksize;
 A=99999999999999999999;         %为了找到最小的均方误差，A用于设定一个很大的初值
 Eij=0;
 xrecord=ones（1616）;            %xrecord，yrecord用于存放匹配块的块号，即运动矢量
 yrecord=ones（1616）;
 diff=zeros（256256）;           %这幅图像的大小为256*256，diff用于存放像素差值
 tic
                for x=0:（rowblocks-1）         %x表示行中第几个子块
                    row=x*blocksize;
                   for y=0:（colblocks-1）         %y表示列中第几个子块
                       col=y*blocksize;
%                        tempx=x*blocksize+1:（x+1）*blocksize;
%                        tempy=y*blocksize+1:（y+1）*blocksize;

                                 for p=-dm:dm
                                    for q=-dm:dm      %（p，q）表示在x，y对应子块在前一帧中所的搜索位置
                                        Eij=0;
                                        Eij=sum（sum（（I2（row+1:row+blocksizecol+1:col+blocksize）-II（row+dm+p+1:row+dm+p+blocksizecol+dm+q+1:col+dm+q+blocksize））.^2））/（blocksize^2）; 
                                        if Eij                                                    A=Eij;
                                                    xrecord（x+1y+1）=p;
                                                    yrecord（x+1y+1）=q;
                                            end   
                                    end
                                end
                                A=999999

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      66614  2006-05-19 07:56  图像通信中的全搜索与三步搜索算法的matlab实现\claire1.BMP

     文件      66614  2006-05-19 07:56  图像通信中的全搜索与三步搜索算法的matlab实现\claire2.BMP

     文件       6353  2009-05-23 12:58  图像通信中的全搜索与三步搜索算法的matlab实现\FS.m

     文件      14847  2009-05-23 11:33  图像通信中的全搜索与三步搜索算法的matlab实现\hs_err_pid2092.log

     文件      14973  2009-04-26 08:48  图像通信中的全搜索与三步搜索算法的matlab实现\hs_err_pid2240.log

     文件      14805  2009-04-26 06:00  图像通信中的全搜索与三步搜索算法的matlab实现\hs_err_pid2520.log

     文件      14929  2009-04-26 07:14  图像通信中的全搜索与三步搜索算法的matlab实现\hs_err_pid608.log

     文件      14714  2009-04-26 06:44  图像通信中的全搜索与三步搜索算法的matlab实现\hs_err_pid888.log

     文件       6764  2009-04-26 09:29  图像通信中的全搜索与三步搜索算法的matlab实现\SEE.m

     文件       7571  2009-05-23 12:57  图像通信中的全搜索与三步搜索算法的matlab实现\S_3SS.m

     文件     389632  2009-05-23 12:03  图像通信中的全搜索与三步搜索算法的matlab实现\测试结果文档.doc

     文件        138  2010-11-24 19:54  图像通信中的全搜索与三步搜索算法的matlab实现\说明.txt

     目录          0  2010-11-24 20:00  图像通信中的全搜索与三步搜索算法的matlab实现

----------- ---------  ---------- -----  ----

               617954                    13