视音频处理技术实验3.doc

浙江理工大学视音频信号处理 实验报告 20132014学年第一学期 学 院信息学院班 级数字媒体与技术11(1)姓 名张志辉学 号2011329700132任课教师张海翔（上课时间：周五1-2节） 数字媒体技术专业 2013年 11 月视音频信号处理实 验 报 告学号：2011329700132 姓名：张志辉 班级：数字媒体与技术11(1)一、实验名称视频数据处理二、实验内容1、请编写图像矩阵产生程序，具体要求为：a) 共产生两个图像矩阵A256256,和B256256b) 图像矩阵中的点A(I,j) = i*256+j,B(I,j) = i*256+j,；c) 特别地，当25=i25+16, 20=j20+16, A(i,j) = i - 20;d) 特别地，当32=i32+16, 16=j16+16, B(i,j) = i - 27;2、 以题1中得到的图像数据基础上，以图像A为参考帧，请用全搜索方法计算图像B中16*16图像块MBlock(32,16)的运动矢量。3、 以题1中得到的图像数据基础上，以图像A为参考帧，请用三步搜索方法计算图像B中16*16图像块MBlock(32,16)的运动矢量。4、 给定一个8*8的DCT系数矩阵C88,12-8-10000000-100000-1-10000000-100000000000000000000000000000000000000请编写程序将DCT系数按Zigzag扫描顺序转换为系数序列，并用游程编码对系数序列进行编码，以的格式输出的中间符号。三、实验步骤14实验1：/ Ani.cpp : 定义应用程序的入口点。/#include stdafx.h#include Ani.h#define MAX_LOADSTRING 100/ 全局变量:HINSTANCE hInst;/ 当前实例TCHAR szTitleMAX_LOADSTRING;/ 标题栏文本TCHAR szWindowClassMAX_LOADSTRING;/ 主窗口类名/ 此代码模块中包含的函数的前向声明:ATOMMyRegisterClass(HINSTANCE hInstance);BOOLInitInstance(HINSTANCE, int);LRESULT CALLBACKWndProc(HWND, UINT, WPARAM, LPARAM);INT_PTR CALLBACKAbout(HWND, UINT, WPARAM, LPARAM);int APIENTRY _tWinMain(_In_ HINSTANCE hInstance, _In_opt_ HINSTANCE hPrevInstance, _In_ LPTSTR lpCmdLine, _In_ int nCmdShow)UNREFERENCED_PARAMETER(hPrevInstance);UNREFERENCED_PARAMETER(lpCmdLine); / TODO: 在此放置代码。MSG msg;HACCEL hAccelTable;/ 初始化全局字符串LoadString(hInstance, IDS_APP_TITLE, szTitle, MAX_LOADSTRING);LoadString(hInstance, IDC_ANI, szWindowClass, MAX_LOADSTRING);MyRegisterClass(hInstance);/ 执行应用程序初始化:if (!InitInstance (hInstance, nCmdShow)return FALSE;hAccelTable = LoadAccelerators(hInstance, MAKEINTRESOURCE(IDC_ANI);/ 主消息循环:while (GetMessage(&msg, NULL, 0, 0)if (!TranslateAccelerator(msg.hwnd, hAccelTable, &msg)TranslateMessage(&msg);DispatchMessage(&msg);return (int) msg.wParam;/ 函数: MyRegisterClass()/ 目的: 注册窗口类。/ATOM MyRegisterClass(HINSTANCE hInstance)WNDCLASSEX wcex;wcex.cbSize = sizeof(WNDCLASSEX);wcex.style= CS_HREDRAW | CS_VREDRAW;wcex.lpfnWndProc= WndProc;wcex.cbClsExtra= 0;wcex.cbWndExtra= 0;wcex.hInstance= hInstance;wcex.hIcon= LoadIcon(hInstance, MAKEINTRESOURCE(IDI_ANI);wcex.hCursor= LoadCursor(NULL, IDC_ARROW);wcex.hbrBackground= (HBRUSH)(COLOR_WINDOW+1);wcex.lpszMenuName= MAKEINTRESOURCE(IDC_ANI);wcex.lpszClassName= szWindowClass;wcex.hIconSm= LoadIcon(wcex.hInstance, MAKEINTRESOURCE(IDI_SMALL);return RegisterClassEx(&wcex);/ 函数: InitInstance(HINSTANCE, int)/ 目的: 保存实例句柄并创建主窗口/ 注释:/ 在此函数中，我们在全局变量中保存实例句柄并/ 创建和显示主程序窗口。/BOOL InitInstance(HINSTANCE hInstance, int nCmdShow) HWND hWnd; hInst = hInstance; / 将实例句柄存储在全局变量中 hWnd = CreateWindow(szWindowClass, szTitle, WS_OVERLAPPEDWINDOW, CW_USEDEFAULT, 0, CW_USEDEFAULT, 0, NULL, NULL, hInstance, NULL); if (!hWnd) return FALSE; ShowWindow(hWnd, nCmdShow); UpdateWindow(hWnd); return TRUE;/ 函数: WndProc(HWND, UINT, WPARAM, LPARAM)/ 目的: 处理主窗口的消息。/ WM_COMMAND- 处理应用程序菜单/ WM_PAINT- 绘制主窗口/ WM_DESTROY- 发送退出消息并返回/LRESULT CALLBACK WndProc(HWND hWnd, UINT message, WPARAM wParam, LPARAM lParam)int wmId, wmEvent;PAINTSTRUCT ps;HDC hdc;switch (message)case WM_COMMAND:wmId = LOWORD(wParam);wmEvent = HIWORD(wParam);/ 分析菜单选择:switch (wmId)case IDM_ABOUT:DialogBox(hInst, MAKEINTRESOURCE(IDD_ABOUTBOX), hWnd, About);break;case IDM_EXIT:DestroyWindow(hWnd);break;default:return DefWindowProc(hWnd, message, wParam, lParam);break;case WM_PAINT:hdc = BeginPaint(hWnd, &ps);char A256256,B256256;int j; FILE *fa, *fb;for(int i=0;i256;i+)for(j=0; j256;j+)if(i=25)&(j=20)Aij=i-20;SetPixel(hdc,i,j,i-20);elseAij=i*256+j;SetPixel(hdc,i,j,i*256+j);if(i=32)&(j=16)Bij=i-27;SetPixel(hdc,i+300,j,i-27);elseBij=i*256+j;SetPixel(hdc,i+300,j,i*256+j);fa=fopen(A.raw,wb);fb=fopen(B.raw,wb);fwrite(A,sizeof(A),1,fa);fwrite(B,sizeof(B),1,fb);fclose(fa);fclose(fb);EndPaint(hWnd, &ps);break;case WM_DESTROY:PostQuitMessage(0);break;default:return DefWindowProc(hWnd, message, wParam, lParam);return 0;/ “关于”框的消息处理程序。INT_PTR CALLBACK About(HWND hDlg, UINT message, WPARAM wParam, LPARAM lParam)UNREFERENCED_PARAMETER(lParam);switch (message)case WM_INITDIALOG:return (INT_PTR)TRUE;case WM_COMMAND:if (LOWORD(wParam) = IDOK | LOWORD(wParam) = IDCANCEL)EndDialog(hDlg, LOWORD(wParam);return (INT_PTR)TRUE;break;return (INT_PTR)FALSE;实验2：#include #include #include int main()int A256256,B256256;int i,j,a1,b1,a,b,c,d,p=0;for(i=0;i256;i+)for(j=0;j=20 & j=25 & i=16 & j=32 & i48) Bij=i-27;else Bij=i*256+j;a=32;b=16;a1=32;b1=16;for(i=0;i240;i+)for(j=0;j240;j+) if(abs(B3216-Aij)1)&(abs(B4731-Ai+15j+15)1)&(abs(B4716-Ai+15j)1)&(abs(B3231-Aij+15)1) a1=i; b1=j; p=1; if(p=1) break;printf(运动矢量为:(%d,%d)n,32-a1,16-b1);return 0;实验3：#include #include #include int main()int A256256,B256256;int i,j,a1,b1,a,b,c,d,count=1;for(i=0;i256;i+)for(j=0;j=20 & j=25 & i=16 & j=32 & i0;k-)d=100;printf(第%d步搜索:n,count);for(i=a1-k;i=(a1-k)&i=(b1-k)&j=(b1+k);j=j+k) if(i!=a1|j!=b1)printf(A%d%d=%d ,i,j,Aij);c=B3216-Aij;if(abs(c)d | d=100)d=abs(c);a=i;b=j;else printf( );printf(n);a1=a;b1=b;count+;printf(n);printf(运动矢量为:(%d,%d)n,32-a1,16-b1);return 0;实验4：#include#include#define N 8using namespace std;void Zigzag(int aNN) int squa=64; for (int i=0;iN;i+) for (int j=0;jN;j+) int s=i+j; if(sN) aij=s*(s+1)/2+(s%2!=0)?i:j); else int sn = (N-1-i) + (N-1-j); aij=squa-sn*(sn+1)/2-(N-(sn%2!=0)? i:j); void Rle(int bN*N) int i,x,m=0; for(i=0;i64;i+) if(bi!=0) x=bi; i+; while(bi=0) m+; i+; coutx,m ; i-; m=0; int main() int ANN=0; int BN*N=0; int CNN=12,-8,-1, 0, 0, 0, 0, 0, 0, 0,-1, 0, 0, 0, 0, 0, -1,-1, 0, 0, 0, 0, 0, 0, 0,-1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,