各种基本运动估计算法代码

各种运动估计算法基础，基本思想是将图像序列的每一帧分成许多互不重叠的宏块，并认为宏块内所有象素的位移量都相同，然后对每个宏块到参考帧某一给定特定搜索范围内根据一定的匹配准则找出与当前块最相似的块，即匹配块，匹配块与当前块的相对位移即为运动矢量。视频压缩的时候，只需保存运动矢量和残差数据就可以完全恢复出当前块。

#include “stdio.h“
#include “stdlib.h“
#include “malloc.h“
#include “string.h“
#include “sys/timeb.h“
#include “math.h“
#include “inf.h“
#define MIN（ab） （（a）<（b））?（a）:（b）;
#define MAX（ab） （（a）>（b））?（a）:（b）;
int SAD（const int oxconst int oyconst int dxconst int dyconst int heightconst int widthuint32 &best_sad）
{
	const int rx=ox+dxry=oy+dy;
	if（ abs（dx）>MAX_MOTION || abs（dy）>MAX_MOTION || flag_search[dx][dy]>0 ）
		return 0;
	if（ rx<0 || ry<0 || rx+height>XX || ry+width>YY ）
		return 0;
	uint32 sad=0;
	frame_info.frame_pot++;
	const uint8 *p1=¤t_frame[ox][oy]*p2=&ref_frame[rx][ry];
	for（int i=0;i	{
		for（int j=0;j		{
			sad+=abs（*（p1++）-*（p2++））;
		}
		p1+=（YY-width）;p2+=（YY-width）;
	}
	flag_search[dx][dy]=1+sad;
	if（sad>=best_sad）
		return 0;
	best_sad=sad;
	return 1;
}
void rebuilt（const int xconst int yconst int heightconst int width）
{
	const int dx=frame_info.mv[x][y].dxdy=frame_info.mv[x][y].dy;
	const int ox=x*BLOCK_HEIGTHoy=y*BLOCK_WIDTH;
	const int rx=ox+dxry=oy+dy;
	int tmp;
	const uint8 *p1=¤t_frame[ox][oy]*p2=&ref_frame[rx][ry];
	if（ abs（dx）>MAX_MOTION || abs（dy）>MAX_MOTION ）
		exit（0）;
	if（ rx<0 || ry<0 || rx+height>XX || ry+width>YY ）
		exit（0）;
	for（int i=0;i	{
		for（int j=0;j		{
			tmp=*（p1++）-*（p2++）;
			frame_info.frame_sse+=tmp*tmp;
		}
		p1+=（YY-width）;p2+=（YY-width）;	
	}
}
void search_FS（const int xconst int yconst int heigthconst int width）
{
	const int ox=x*BLOCK_HEIGTHoy=y*BLOCK_WIDTH;
	uint32 sad=0xffffff;
	MV mv;
	for（int i=-MAX_MOTION;i<=MAX_MOTION;i++）
	{
		for（int j=-MAX_MOTION;j<=MAX_MOTION;j++）
		{
			if（SAD（oxoyijheigthwidthsad）==1）
			{
				mv.dx=i;mv.dy=j;
			}
		}
	}
	frame_info.mv[x][y]=mv;
	frame_info.sad[x][y]=sad;
	frame_info.frame_sad+=sad;
}
#define PATTERN_SEARCH（patternnumflag） \
	do\
	{\
		mvx=mv.dx;mvy=mv.dy;\
		for（int i=0;i		{\
			if（SAD（oxoymvx+pattern[i][0]mvy+pattern[i][1]heigthwidthsad）==1）\
			{\
				mv.dx=mvx+pattern[i][0];mv.dy=mvy+pattern[i][1];\
			}\
		}\
	}while （ （mv.dx!=mvx || mv.dy!=mvy） && flag==1 ）;
void search_4SS（const int xconst int yconst int heigthconst int width）
{
	const int ox=x*BLOCK_HEIGTHoy=y*BLOCK_WIDTH;
	const int L4SS[9][2]={{00}{02}{-22}{-20}{-2-2}{0-2}{2-2}{20}{22}};
	const int S4SS[9][2]={{00}{01}{-11}{-10}{-1-1}{0-1}{1-1}{10}{11}};
	uint32 sad=0xffffff;
	MV mv={00};int mvxmvy;

	PATTERN_SEARCH（L4SS91）
	PATTERN_SEARCH（S4SS90）

	frame_info.mv[x][y]=mv;
	frame_info.sad[x][y]=sad;
	frame_info.frame_sad+=sad;
}
void search_BBGDS（const int xconst int yconst int heigthconst int width）
{
	const int ox=x*BLOCK_HEIGTHoy=y*BLOCK_WIDTH;
	const int BBGDS[9][2]={{00}{01}{-11}{-10}{-1-1}{0-1}{1-1}{10}{11}};
	uint32 sad=0xffffff;
	MV mv={00};int mvxmvy;

	PATTERN_SEARCH（BBGDS91）

	frame_info.mv[x][y]=mv

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

     文件      10312  2011-09-21 09:11  运动估计方法比较\Debug\BuildLog.htm

     文件      60828  2011-09-21 09:11  运动估计方法比较\Debug\main.obj

     文件          0  2010-06-25 12:41  运动估计方法比较\Debug\main.sbr

     文件         65  2011-09-21 09:11  运动估计方法比较\Debug\mt.dep

     文件      41984  2011-03-12 16:33  运动估计方法比较\Debug\vc60.idb

     文件      53248  2011-03-12 16:33  运动估计方法比较\Debug\vc60.pdb

     文件      60416  2011-09-17 16:53  运动估计方法比较\Debug\vc80.idb

     文件      61440  2011-09-17 16:53  运动估计方法比较\Debug\vc80.pdb

     文件      82944  2010-06-25 12:41  运动估计方法比较\Debug\zhang.bsc

     文件     572928  2011-09-21 09:11  运动估计方法比较\Debug\zhang.exe

     文件        406  2011-09-19 15:16  运动估计方法比较\Debug\zhang.exe.embed.manifest

     文件        472  2011-09-19 15:16  运动估计方法比较\Debug\zhang.exe.embed.manifest.res

     文件        381  2011-09-21 09:11  运动估计方法比较\Debug\zhang.exe.intermediate.manifest

     文件    1187036  2011-09-21 09:11  运动估计方法比较\Debug\zhang.ilk

     文件       1118  2011-09-19 15:16  运动估计方法比较\inf.h

     文件      18260  2011-09-21 09:11  运动估计方法比较\main.cpp

     文件       4329  2006-05-20 14:41  运动估计方法比较\zhang.dsp

     文件        535  2006-05-20 10:07  运动估计方法比较\zhang.dsw

     文件     723968  2011-09-21 15:03  运动估计方法比较\zhang.ncb

     文件      48640  2011-01-06 13:35  运动估计方法比较\zhang.opt

     文件       3377  2011-03-12 16:33  运动估计方法比较\zhang.plg

     文件        875  2011-09-19 15:15  运动估计方法比较\zhang.sln

    ..A..H.     11776  2011-09-21 15:03  运动估计方法比较\zhang.suo

     文件       5212  2011-09-19 15:15  运动估计方法比较\zhang.vcproj

     文件       1413  2011-09-19 10:49  运动估计方法比较\zhang.vcproj.admin-PC.admin.user

     文件       1241  2011-09-17 17:11  运动估计方法比较\说明.txt

     目录          0  2011-09-21 09:11  运动估计方法比较\Debug

     目录          0  2011-09-21 09:11  运动估计方法比较

     文件      60416  2011-09-21 09:11  运动估计方法比较\Debug\vc90.idb

     文件      69632  2011-09-21 09:11  运动估计方法比较\Debug\vc90.pdb

............此处省略14个文件信息