DlgHuffman.cpp
资源名称:图像识别.rar [点击查看]
上传用户:xjt2008yy
上传日期:2010-01-18
资源大小:272k
文件大小:5k
源码类别:
生物技术
开发平台:
Visual C++
- // DlgCodingHuffman.cpp : implementation file
- //
- #include "stdafx.h"
- #include "imageProcessing.h"
- #include "DlgCoding.h"
- #include <math.h>
- #ifdef _DEBUG
- #define new DEBUG_NEW
- #undef THIS_FILE
- static char THIS_FILE[] = __FILE__;
- #endif
- /////////////////////////////////////////////////////////////////////////////
- // CDlgHuffman dialog
- CDlgHuffman::CDlgHuffman(CWnd* pParent /*=NULL*/)
- : CDialog(CDlgHuffman::IDD, pParent)
- {
- //{{AFX_DATA_INIT(CDlgHuffman)
- m_dEntropy = 0.0;
- m_dCodLength = 0.0;
- m_dRatio = 0.0;
- //}}AFX_DATA_INIT
- }
- void CDlgHuffman::DoDataExchange(CDataExchange* pDX)
- {
- CDialog::DoDataExchange(pDX);
- //{{AFX_DATA_MAP(CDlgHuffman)
- DDX_Control(pDX, IDC_LST_Table, m_lstTable);
- DDX_Text(pDX, IDC_EDIT1, m_dEntropy);
- DDX_Text(pDX, IDC_EDIT2, m_dCodLength);
- DDX_Text(pDX, IDC_EDIT3, m_dRatio);
- //}}AFX_DATA_MAP
- }
- BEGIN_MESSAGE_MAP(CDlgHuffman, CDialog)
- //{{AFX_MSG_MAP(CDlgHuffman)
- //}}AFX_MSG_MAP
- END_MESSAGE_MAP()
- /////////////////////////////////////////////////////////////////////////////
- // CDlgHuffman message handlers
- BOOL CDlgHuffman::OnInitDialog()
- {
- // 调用默认得OnInitDialog()函数
- CDialog::OnInitDialog();
- // 循环变量
- LONG i;
- LONG j;
- LONG k;
- // 中间变量
- double dT;
- // 字符串变量
- CString str2View;
- // 控件ListCtrl的ITEM
- LV_ITEM lvItem;
- // 保存控件ListCtrl中添加的ITEM编号
- int nItem2View;
- // 保存计算中间结果的数组
- double * dTemp;
- // 数组用来存放灰度值和其位置之间的映射
- int * n4Turn;
- // 初始化变量
- m_dEntropy = 0.0;
- m_dCodLength = 0.0;
- // 分配内存
- m_strCode = new CString[nColorNum];
- n4Turn = new int[nColorNum];
- dTemp = new double[nColorNum];
- // 令dTemp赋值
- // 开始的灰度值按灰度值大小排列相同
- for (i = 0; i < nColorNum; i ++)
- {
- dTemp[i] = dProba[i];
- n4Turn[i] = i;
- }
- // 用冒泡法对进行灰度值出现的概率排序
- // 同时改变灰度值位置的映射关系
- for (j = 0; j < nColorNum - 1; j ++)
- {
- for (i = 0; i < nColorNum - j - 1; i ++)
- {
- if (dTemp[i] > dTemp[i + 1])
- {
- dT = dTemp[i];
- dTemp[i] = dTemp[i + 1];
- dTemp[i + 1] = dT;
- // 将i和i+1灰度的位置值互换
- for (k = 0; k < nColorNum; k ++)
- {
- if (n4Turn[k] == i)
- n4Turn[k] = i + 1;
- else if (n4Turn[k] == i + 1)
- n4Turn[k] = i;
- }
- }
- }
- }
- /*******************************************************
- 计算哈夫曼编码表
- *******************************************************/
- // 从概率大于0处开始编码
- for (i = 0; i < nColorNum - 1; i ++)
- {
- if (dTemp[i] > 0)
- break;
- }
- for (; i < nColorNum - 1; i ++)
- {
- // 更新m_strCode
- for (k = 0; k < nColorNum; k ++)
- {
- // 灰度值是否i
- if (n4Turn[k] == i)
- {
- // 灰度值较小的码字加1
- m_strCode[k] = "1" + m_strCode[k];
- }
- else if (n4Turn[k] == i + 1)
- {
- // 灰度值较小的码字加0
- m_strCode[k] = "0" + m_strCode[k];
- }
- }
- // 概率最小的两个概率相加,保存在dTemp[i + 1]中
- dTemp[i + 1] += dTemp[i];
- // 改变映射关系
- for (k = 0; k < nColorNum; k ++)
- {
- // 将位置为i的灰度值i改为灰度值i+1
- if (n4Turn[k] == i)
- n4Turn[k] = i + 1;
- }
- // 重新排序
- for (j = i + 1; j < nColorNum - 1; j ++)
- {
- if (dTemp[j] > dTemp[j + 1])
- {
- // 互换
- dT = dTemp[j];
- dTemp[j] = dTemp[j + 1];
- dTemp[j + 1] = dT;
- // // 将i和i+1灰度的位置值互换
- for (k = 0; k < nColorNum; k ++)
- {
- if (n4Turn[k] == j)
- n4Turn[k] = j + 1;
- else if (n4Turn[k] == j + 1)
- n4Turn[k] = j;
- }
- }
- else
- // 退出循环
- break;
- }
- }
- // 计算图像熵
- for (i = 0; i < nColorNum; i ++)
- {
- if (dProba[i] > 0)
- {
- // 计算图像熵
- m_dEntropy -= dProba[i] * log(dProba[i]) / log(2.0);
- }
- }
- // 计算平均码字长度
- for (i = 0; i < nColorNum; i ++)
- {
- // 累加
- m_dCodLength += dProba[i] * m_strCode[i].GetLength();
- }
- // 计算编码效率
- m_dRatio = m_dEntropy / m_dCodLength;
- // 保存变动
- UpdateData(FALSE);
- /*************************************************
- 输出计算结果
- *************************************************/
- // 设置CListCtrl控件样式
- m_lstTable.ModifyStyle(LVS_TYPEMASK, LVS_REPORT);
- // 给List控件添加Header
- m_lstTable.InsertColumn(0, "灰度值", LVCFMT_LEFT, 60, 0);
- m_lstTable.InsertColumn(1, "概率值", LVCFMT_LEFT, 78, 0);
- m_lstTable.InsertColumn(2, "哈夫曼编码", LVCFMT_LEFT, 110, 1);
- m_lstTable.InsertColumn(3, "码字长度", LVCFMT_LEFT, 78, 2);
- // 设置样式为文本
- lvItem.mask = LVIF_TEXT;
- // 添加显示
- for (i = 0; i < nColorNum; i ++)
- {
- // 第一列显示
- lvItem.iItem = m_lstTable.GetItemCount();
- str2View.Format("%u",i);
- lvItem.iSubItem = 0;
- lvItem.pszText= (LPTSTR)(LPCTSTR)str2View;
- nItem2View = m_lstTable.InsertItem(&lvItem);
- // 其它列显示
- lvItem.iItem = nItem2View;
- // 添加灰度值的概率值
- lvItem.iSubItem = 1;
- str2View.Format("%f",dProba[i]);
- lvItem.pszText = (LPTSTR)(LPCTSTR)str2View;
- m_lstTable.SetItem(&lvItem);
- // 添加哈夫曼编码
- lvItem.iSubItem = 2;
- lvItem.pszText = (LPTSTR)(LPCTSTR)m_strCode[i];
- m_lstTable.SetItem(&lvItem);
- // 添加码字长度
- lvItem.iSubItem = 3;
- str2View.Format("%u",m_strCode[i].GetLength());
- lvItem.pszText = (LPTSTR)(LPCTSTR)str2View;
- m_lstTable.SetItem(&lvItem);
- }
- // 内存释放
- delete n4Turn;
- delete dTemp;
- // 返回TRUE
- return TRUE;
- }