《第四次实验》word版.doc

数字图像处理实验报告（四）一 实验目的：熟悉Visual c+的使用, 学会用哈夫曼编码程序对图像进行处理二 实验过程：（1） 程序编程 void CImgcxView:OnHuffmanCoding() / 循环变量if (!checkInit()return ;LONG i;LONG j;LONG k;/ 灰度计数int nNs256;/ 灰度概率分布float fPs256;/ 映射关系int iMap256;/ Huffman编码CString m_strCode 256;/ 变量初始化memset(nNs, 0, sizeof(nNs);/ 指向DIB的指针/LPBYTE lpDIB;/ 指向DIB象素指针/LPBYTE lpDIBBits;/ 锁定DIB/lpDIB = (LPBYTE) :GlobalLock(HGLOBAL) hDIB);/ 对各像素进行灰度分布统计for (i = 0; i m_nYSize; i +)for (j = 0; j m_nXSize; j +)/ 对各像素进行灰度统计unsigned char V = *(unsigned char *)m_pImage + m_nXSize * i + j);nNsV+;/ 计算灰度分布密度for(i = 0; i 256; i+)fPsi = nNsi / (m_nYSize * m_nXSize * 1.0f);/ 初始化for (i = 0; i 256; i +)iMapi = i;/ 用冒泡法对fPs进行排序for (j = 0; j 256 - 1; j +)for (i = 0; i fPsi + 1)/ 互换float fTemp = fPsi;fPsi = fPsi + 1;fPsi + 1 = fTemp;/ 更新映射关系for (k = 0; k 256; k +)/ 判断是否是fPsi的子节点if (iMapk = i)/ 映射到节点i+1iMapk = i + 1;else if (iMapk = i + 1)/ 映射到节点iiMapk = i;/ 计算哈夫曼编码表for (i = 0; i 0)break;/ 开始编码for (i = i; i 256 - 1; i +)/ 更新m_strCodefor (k = 0; k 256; k +)/ 判断是否是fPsi的子节点并编码字符串if (iMapk = i)m_strCodek = 1 + m_strCodek;else if (iMapk = i + 1)m_strCodek = 0 + m_strCodek;/ 概率最小的两个概率相加，结果保存到fPsi + 1fPsi + 1 += fPsi;/ 改变映射关系for (k = 0; k 256; k +)/ 判断是否是fPsi的子节点if (iMapk = i)/ 映射到节点i+1iMapk = i + 1;/ 重新排序for (j = i + 1; j fPsj + 1)/ 互换float fTemp = fPsj;fPsj = fPsj + 1;fPsj + 1 = fTemp;/ 更新映射关系for (k = 0; k 256; k +)/ 判断是否是fPsi的子节点if (iMapk = j)/ 映射到节点j+1iMapk = j + 1;else if (iMapk = j + 1)/ 映射到节点jiMapk = j;else/ 退出循环break;/ 对DIB进行编码压缩/ 位掩码unsigned char Mask8 = 128, 64, 32, 16, 8, 4, 2, 1;CString strTemp = ;/ 打开文件CFile file;file.Open(XXX, CFile:modeCreate | CFile:modeReadWrite);/ 对各像素进行编码for (i = 0; i m_nYSize; i +)for (j = 0; j 256)dounsigned char T = 0;for (k = 0; k min(8, strTemp.GetLength(); k+)if (strTemp.Mid(k, 1) = 1)T |= Maskk;file.Write(&T, 1);loop+;len -= 8;if (strTemp.GetLength() =0);strTemp=;/ 关闭文件file.Close();/ 恢复光标EndWaitCursor();/ 返回TRUEreturn ;/ TODO: Add your command handler code here/bool CImgcxView:checkInit(void) if(!m_bFileOpen)/没有读入图像 AfxMessageBox(Please load a grey image);return 0;if(:GetImageType() != 8)AfxMessageBox(Please