请各位达人帮下忙，小弟问下灰度图象处理的实现 - VC++/MFC - C++论坛

得分:0

　　一、 BMP位图操作
　　BMP位图包括位图文件头结构BITMAPFILEHEADER、位图信息头结构BITMAPINFOHEADER、位图颜色表RGBQUAD和位图像素数据四部分。处理位图时要根据文件的这些结构得到位图文件大小、位图的宽、高、实现调色板、得到位图像素值等等。对于256级灰度图像每个像素用8bit表示颜色的索引值，这里要注意的一点是在BMP位图中，位图的每行像素值要填充到一个四字节边界，即位图每行所占的存储长度为四字节的倍数，不足时将多余位用0填充。
　　在处理图像应用程序的文档类（CdibDoc.h）中声明如下宏及公有变量：
　　#define WIDTHBYTES(bits) (((bits) + 31) / 32 * 4)//计算图像每行象素所占的字节数目
　　HANDLE m_hDIB;//存放位图数据的句柄
　　CPalette* m_palDIB;//指向调色板Cpalette类的指针
　　CSize m_sizeDoc; //初始化视图的尺寸
　　1、读取灰度BMP位图
　　根据BMP位图文件的结构，操作BMP位图文件读入数据，重载了文挡类的OnOpenDocument函数如下：

BOOL CDibDoc::OnOpenDocument(LPCTSTR lpszPathName)
{
　CFile file;
　CFileException fe;
　if (!file.Open(lpszPathName, CFile::modeRead | CFile::shareDenyWrite, &fe))
　{
　　AfxMessageBox("文件打不开");
　　return FALSE;
　　}//打开文件
　DeleteContents();//删除文挡
　BeginWaitCursor();
　BITMAPFILEHEADER bmfHeader;//定义位图文件头结构
　DWORD dwBitsSize;
　HANDLE hDIB;
　LPSTR pDIB;
　BITMAPINFOHEADER *bmhdr;//指向位图信息头结构的指针
　dwBitsSize = file.GetLength();//得到文件长度
　if (file.Read((LPSTR)&bmfHeader, sizeof(bmfHeader)) !=
　　　　　　　　　　　　　sizeof(bmfHeader))
　　return FALSE;
　　if (bmfHeader.bfType != 0x4d42) //检查是否为BMP文件
　　　return FALSE;
　　　hDIB=(HANDLE) ::GlobalAlloc(GMEM_MOVEABLE |
　　　　　　　　　　　　　　　　 GMEM_ZEROINIT, dwBitsSize);
　　　//申请缓冲区
　　　if (hDIB == 0)
　　　{
　　　　return FALSE;
　　　}
　　　pDIB = (LPSTR) ::GlobalLock((HGLOBAL)hDIB);
　　　//得到申请的缓冲区的指针
　　　if (file.ReadHuge(pDIB, dwBitsSize - sizeof(BITMAPFILEHEADER)) !=
　　　　　dwBitsSize - sizeof(BITMAPFILEHEADER) )
　　　{
　　　　::GlobalUnlock((HGLOBAL)hDIB);
　　　　hDIB=NULL;
　　　　return FALSE;
　　　　}//读数据，包括位图信息、位图颜色表、图像像素的灰度值
　　　bmhdr=(BITMAPINFOHEADER*)pDIB;//为指向位图信息头结构的指针付值
　　　::GlobalUnlock((HGLOBAL)hDIB);
　　　if ((*bmhdr).biBitCount!=8) //验证是否为8bit位图
　　　　return FALSE;
　　　　m_hDIB=hDIB;
　　　　InitDIBData();
　　　　//自定义函数，根据读入的数据得到位图的宽、高、颜色表
　　　　// 来得到初始化视的尺寸、生成调色板
　　　　EndWaitCursor();
　　　　SetPathName(lpszPathName);//设置存储路径
　　　　SetModifiedFlag(FALSE); // 设置文件修改标志为FALSE
　　　　return TRUE;
　　　}

2、灰度位图数据的存储
　　为了将图像处理后所得到的像素值保存起来，重载了文档类的OnSaveDocument函数，其具体实现如下：

BOOL CDibDoc::OnSaveDocument(LPCTSTR lpszPathName)
{
　CFile file;
　CFileException fe;
　BITMAPFILEHEADER bmfHdr; // 位图文件头结构
　LPBITMAPINFOHEADER lpBI; //指向位图信息结构的指针
　DWORD dwDIBSize;
　if (!file.Open(lpszPathName, CFile::modeCreate |
　　　CFile::modeReadWrite | CFile::shareExclusive, &fe))
　{
　　AfxMessageBox("文件打不开");
　}//打开文件
　BOOL bSuccess = FALSE;
　BeginWaitCursor();
　lpBI = (LPBITMAPINFOHEADER) ::GlobalLock((HGLOBAL) m_hDIB);
　if (lpBI == NULL)
　　return FALSE;
　　dwDIBSize = *(LPDWORD)lpBI + 256*sizeof(RGBQUAD);
　　 // Partial Calculation
　　DWORD dwBmBitsSize;//BMP文件信息结构所占的字节数
　　dwBmBitsSize=WIDTHBYTES((lpBI->biWidth)*((DWORD)lpBI->biBitCount)) 　*lpBI->biHeight;// 存储时位图所有像素所占的总字节数
　　dwDIBSize += dwBmBitsSize;
　　lpBI->biSizeImage = dwBmBitsSize; // 位图所有像素所占的总字节数
　　//以下五句为文件头结构填充值
　　bmfHdr.bfType =0x4d42; // 文件为"BMP"类型
　　bmfHdr.bfSize = dwDIBSize + sizeof(BITMAPFILEHEADER);//文件总长度
　　bmfHdr.bfReserved1 = 0;
　　bmfHdr.bfReserved2 = 0;
　　bmfHdr.bfOffBits = (DWORD)sizeof(BITMAPFILEHEADER) + lpBI->biSize
　　　　　　　　　　　　+ 256*sizeof(RGBQUAD);
　　//位图数据距问件头的偏移量
　　file.Write((LPSTR)&bmfHdr, sizeof(BITMAPFILEHEADER));//写文件头
　　file.WriteHuge(lpBI, dwDIBSize);
　　//将位图信息（信息头结构、颜色表、像素数据）写入文件
　　::GlobalUnlock((HGLOBAL) m_hDIB);
　　EndWaitCursor();
　　SetModifiedFlag(FALSE); // back to unmodified
　　return TRUE;
　}

二、调色板的操作
　　灰度图像要正确显示，必须实现逻辑调色板和系统调色板，通过在主框架类中处理Windows定义的消息WM_QUERYNEWPALETTE 、WM_PALETTECHANGED及视图类中处理自定义消息WM_DOREALIZE（该消息在主框架窗口定义如下：#define WM_REALIZEPAL (WM_USER+100)）来实现调色板的操作。

void CMainFrame::OnPaletteChanged(CWnd* pFocusWnd)
{ //总实现活动视的调色板
　CMDIFrameWnd::OnPaletteChanged(pFocusWnd);
　CMDIChildWnd* pMDIChildWnd = MDIGetActive();
　if (pMDIChildWnd == NULL)
　　return
　　CView* pView = pMDIChildWnd->GetActiveView();
　　ASSERT(pView != NULL);
　　SendMessageToDescendants(WM_DOREALIZE, (WPARAM)pView->m_hWnd);
　　//通知所有子窗口系统调色板已改变
　}
BOOL CMainFrame::OnQueryNewPalette()//提供实现系统调色板的机会
　{
　　// 实现活动视的调色板
　　CMDIChildWnd* pMDIChildWnd = MDIGetActive();
　　if (pMDIChildWnd == NULL)
　　　return FALSE; // no active MDI child frame (no new palette)
　　　CView* pView = pMDIChildWnd->GetActiveView();
　　　ASSERT(pView != NULL);
　　　//通知活动视图实现系统调色板
　　　pView->SendMessage(WM_DOREALIZE, (WPARAM)pView->m_hWnd);
　　　return TRUE;
　}
　LRESULT CDibView::OnDoRealize(WPARAM wParam, LPARAM)//实现系统调色板
　{
　　ASSERT(wParam != NULL);
　　CDibDoc* pDoc = GetDocument();
　　if (pDoc->m_hDIB == NULL)
　　　return 0L; // must be a new document
　　　CPalette* pPal = pDoc->m_palDIB;
　　　//调色板的颜色表数据在InitDIBData()函数中实现
　　　if (pPal != NULL)
　　　{
　　　　CMainFrame* pAppFrame = (CMainFrame*) AfxGetApp()->m_pMainWnd;
　　　　ASSERT_KINDOF(CMainFrame, pAppFrame);
　　　　CClientDC appDC(pAppFrame);
　　　　CPalette* oldPalette = appDC.SelectPalette(pPal, ((HWND)wParam) != m_hWnd);
　　　　//只有活动视才可以设为"FALSE"，
　　　　// 即根据活动视的调色板设为"前景"调色板
　　　if (oldPalette != NULL)
　　　　{
　　　　　UINT nColorsChanged = appDC.RealizePalette();//实现系统调色板
　　　　　if (nColorsChanged > 0)
　　　　　　pDoc->UpdateAllViews(NULL);//更新视图
　　　　　　appDC.SelectPalette(oldPalette, TRUE);
　　　　　　//将原系统调色板置为逻辑调色板
　　　　　}
　　　　else
　　　　{
　　　　　TRACE0("\tSelectPalette failed in
　　　　　　　　　 CDibView::OnPaletteChanged\n");
　　　　　}
　　　　}

　　注：在调用API函数显示位图时，不要忘记设置逻辑调色板，即"背景"调色板，否则位图将无法正确显示。

三、图像的数字化处理
　　通过以上读文件的操作，已经得到图像数据，由于得到的数据包括多余信息，所以在进行数字图像处理时要进一步删除多余信息，只对位图的像素进行操作，以基于模板的高通滤波为例来讲述数字图像处理的实现：

void CDibView::OnMENUHighPass()
{ HANDLE data1handle;
　LPBITMAPINFOHEADER lpBi;
　CDibDoc *pDoc=GetDocument();
　HDIB hdib; unsigned char *hData; unsigned char *data;
　hdib=pDoc->GetHDIB();
　BeginWaitCursor();
　lpBi=(LPBITMAPINFOHEADER)GlobalLock((HGLOBAL)hdib);
　hData=(unsigned char*)FindDIBBits((LPSTR)lpBi);
　pDoc->SetModifiedFlag(TRUE);
　data1handle=GlobalAlloc(GMEM_SHARE,WIDTHBYTES(lpBi->biWidth*8)*lpBi->biHeight);
　data=(unsigned char*)GlobalLock((HGLOBAL)data1handle);
　AfxGetApp()->BeginWaitCursor();
　int i,j,s,t,ms=1;
　int sum=0,sumw=0;
　int mask[3][3]={{-1,-1,-1},{-1,9,-1},{-1,-1,-1}};
　for(i=0; ibiHeight; i++)
　　for(j=0; jbiWidth; j++)
　　　{
　　　　sumw=0; sum=0;
　　　　for(s=(-ms); s<=ms; s++)
　　　　　for(t=(-ms); t<=ms; t++)
　　　　　　if(((i+s)>=0) && ((j+t)>=0) && ((i+s)biHeight) && ((j+t)biWidth))
　　　　　　{
　　　　　　　sumw += mask[1+s][1+t];
　　　　　　sum+=*(hData+(i+s)*WIDTHBYTES(lpBi->biWidth*8)+(j+t))*mask[1+s][1+t];
}
　　　　　if(sumw==0) sumw=1; sum/=sumw;
　　　　　if(sum>255)sum=255;
　　　　　if(sum<0)sum=0;
　　　　　　　*(data+i*WIDTHBYTES(lpBi->biWidth*8)+j)=sum;
　　　　　　}
　　　　　for( j=0; jbiHeight; j++)
　　　　　　for( i=0; ibiWidth; i++)
*(hData+i*WIDTHBYTES(lpBi->biWidth*8)+j)=*(data+i*WIDTHBYTES(lpBi->biWidth*8)+j);
　　　　　　AfxGetApp()->EndWaitCursor();
　　　　　　GlobalUnlock((HGLOBAL)hdib);
　　　　　　GlobalUnlock(data1handle);
　　　　　　EndWaitCursor();
　　　　　Invalidate(TRUE);
　}

四、图像的基本操作处理
　　1、图像平移
　　图像平移只是改变图像在屏幕上的位置，图像本身并不发生变化。假设原图像区域左上角坐标为(x0, y0)，右下角坐标为(x1, y1)，将图像分别沿x和y轴平移dx和dy，则新图像的左上角坐标为(x0+dx, y0+dy)，右下角坐标为(x1+dx, y1+dy)。坐标平移变换公式为：
　　x1 = x + dx
　　y1 = y + dy
　　在屏幕上实现图像的移动分为四个步骤：
　　⑴ 读原图像到缓冲区；
　　⑵ 擦除视图上原图像；
　　⑶ 计算平移后的新坐标。
　　⑷ 利用API函数：：StretchDIBits()在新的左上角坐标位置处重新显示原图像。
　　其中，擦除原图像的方法与图形变换中擦除原图形的方法一致，在实现中仍采用XOR异或方式画图擦除原图像。对于新坐标值的计算还需要考虑边界情况，不要在图像平移后超出允许的屏幕范围。
　　2、图像颠倒
　　图像颠倒是指把定义好的图像区域上下翻转地显示在屏幕上。分析图像颠倒的过程，可发现每行的图像信息都保持不变，而只是改变了行的顺序，将第一行与最后的第n行相互交换，第二行与第n － 1行交换……，依此类推，从而实现了图像的颠倒。只需采用按行交换的方式，即可方便地修改缓冲区内容，实现图像的颠倒。基本步骤如下
　　(1)将原图像读入缓冲区，并擦除原图像；
　　(2) 计算图像的高度，即行数height；计算图像宽度width；根据宽度、高度生成新缓冲区；

　　(3)把第一行与最末行交换，第2行与第n－1行交换……，依此类推，直至全部交换完毕。既原图中的(x、y)点，在新生成的图象中对应为x1=x，y1=height-1-y。把原图中的象素值读入新缓冲区的(x1，y1)点处。
　　(4)把交换后的图像缓冲区内容重新显示在屏幕上。
　　3、图像镜像变换
　　镜像变换是指将指定区域的图像左右翻转地显示在屏幕。分析镜像变换过程可以发现：每行图像信息的处理方式是相同的，而且行顺序不发生变化，只是每一行的像素信息按从左到右的顺序进行了左右颠倒，从而实现了镜像变换。因此，采用按行逐点变换的方式实现图像的镜像。
　　给出原图中的任意点（x, y）镜像变换后的新坐标（x1, y1）的坐标变换公式：
　　x1 = width-x-1
　　y1 = y
　　根据以上公式，对各个像素点计算新坐标后，把原图中的象素值读入新缓冲区的(x1，y1)点处。
　　4、图像任意角度的旋转
　　图像旋转是指把定义的图像绕某一点以逆时针或顺时针方向旋转一定的角度，通常是指绕图像的中心以逆时针方向旋转。
　　首先根据旋转的角度、图象对角线的长度计算旋转后的图像的最大宽度、高度，根据旋转后图象最大的宽度、高度生成新的缓冲区，假设图像的左上角为（left, top）,右下角为（right, bottom），则图像上任意点（x, y）绕其中心（xcenter, ycenter）逆时针旋转angle角度后，新的坐标位置（x1, y1）的计算公式为：
　　xcenter = (width+1)/2+left;
　　ycenter = (height+1)/2+top;
　　x1 = (x-xcenter) cosθ- (y - ycenter) sinθ+xcenter;
　　y1 = (x-xcenter) sinθ+ (y- ycenter) cosθ+ ycenter;
　　与图像的镜像变换相类似，把原图中的象素值读入新缓冲区的(x1，y1)点处。注意在新缓冲区中与原图没有对应的象素点的值用白色代替。