VC++图像处理算法(Sobel浮雕效果马赛克高通滤波).doc

/Sobel/具体操作的时候稍微做些修改即可void CDibView:OnMENUSobel() /灰度图像数据的获得 HANDLE data1handle; LPBITMAPINFOHEADER lpBi; CDibDoc *pDoc=GetDocument(); HDIB hdib; unsigned char *hData; unsigned char *data; hdib=pDoc-m_hDIB; BeginWaitCursor(); lpBi=(LPBITMAPINFOHEADER)GlobalLock(HGLOBAL)hdib); hData= lpbi +* (LPDWORD)lpbi + 256*sizeof(RGBQUAD); /得到指向位图像素值的指针 pDoc-SetModifiedFlag(TRUE);/设修改标志为TRUE data1handle=GlobalAlloc(GMEM_SHARE,WIDTHBYTES(lpBi-biWidth*8)*lpBi-biHeight); /申请存放处理后的像素值的缓冲区 data=(unsigned char*)GlobalLock(HGLOBAL)data1handle); AfxGetApp()-BeginWaitCursor(); int i,j,buf,buf1,buf2; for( j=0; jbiHeight; j+)/以下循环求(x,y)位置的灰度值 for( i=0; ibiWidth; i+) if(i-1)=0)&(i+1)biWidth)&(j-1)=0)&(j+1)biHeight) /对于图像四周边界处的向素点不处理 buf1=(int)*(hData+(i+1)*WIDTHBYTES(lpBi-biWidth*8)+(j-1) +2*(int)*(hData+(i+1)*WIDTHBYTES(lpBi-biWidth*8)+(j) +(int)(int)*(hData+(i+1)*WIDTHBYTES(lpBi-biWidth*8)+(j+1); buf1=buf1-(int)(int)*(hData+(i-1)*WIDTHBYTES(lpBi-biWidth*8)+(j-1) -2*(int)(int)*(hData+(i-1)*WIDTHBYTES(lpBi-biWidth*8)+(j) -(int)(int)*(hData+(i-1)*WIDTHBYTES(lpBi-biWidth*8)+(j+1); /x方向加权微分 buf2=(int)(int)*(hData+(i-1)*WIDTHBYTES(lpBi-biWidth*8)+(j+1) +2*(int)(int)*(hData+(i)*WIDTHBYTES(lpBi-biWidth*8)+(j+1) +(int)(int)*(hData+(i+1)*WIDTHBYTES(lpBi-biWidth*8)+(j+1); buf2=buf2-(int)(int)*(hData+(i-1)*WIDTHBYTES(lpBi-biWidth*8)+(j-1) -2*(int)(int)*(hData+(i)*WIDTHBYTES(lpBi-biWidth*8)+(j-1) -(int)(int)*(hData+(i+1)*WIDTHBYTES(lpBi-biWidth*8)+(j-1); /y方向加权微分 buf=abs(buf1)+abs(buf2);/求梯度 if(buf255) buf=255; if(bufbiWidth+j)=(BYTE)0; 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); GlobalFree(date1handle); EndWaitCursor(); Invalidate(TRUE); /Sobel/数字图象的马赛克效果函数/void ImageMasaic(CImage *pImgm,int n) /马赛克效果函数， pImgm 待处理位图指针，n马赛克矩形的大小。 structIMAGEPARAMENT P; /自定义的结构体int Dx,Dy,i,j,s,t;BYTE*buf,*buf1;GetImageParament(pImgm,&P); /读取图像参数（自定义函数） Dx=P.nWidth/n*n; /确定最后点的位置Dy=P.nHeight/n*n; /确定最下行的位置for (j=0;jGetPixelAddress(0,j); /获取行首地址 for (i=0;iDx;i+=n)/循环2 for (s=1;sn;s+)/循环3 memcpy(&buf(i+s)*P.nBytesPerPixel,&bufi*P.nBytesPerPixel,P.nBytesPerPixel); /复制像素值给后面的n-1个点 for (t=1;tGetPixelAddress(0,j+t);/获取行首地址 memcpy(buf1,buf,P.nBytesPerLine);/整行复制像素值 /数字图象的马赛克效果函数/数字图象的马浮雕效果函数/voidCDibView:OnFDImage()/产生浮雕效果图函数HANDLEdata1handle;LPBITMAPINFOHEADERlpBi;CDibDoc*pDoc=GetDocument();HDIBhdib;unsignedchar*hData;unsignedchar*data;hdib=pDoc-GetHDIB(); BeginWaitCursor();lpBi=(LPBITMAPINFOHEADER)GlobalLock(HGLOBAL)hdib);hData=(unsignedchar*)FindDIBBits(LPSTR)lpBi);pDoc-SetModifiedFlag(TRUE);data1handle=GlobalAlloc(GMEM_SHARE,WIDTHBYTES(lpBi-biWidth*8)*lpBi-biHeight);/声明一个缓冲区用来暂存处理后的图像数据data=(unsignedchar*)GlobalLock(HGLOBAL)data1handle);AfxGetApp()-BeginWaitCursor();inti,j,buf;for(i=lpBi-biHeight;i=2;i-)for(j=lpBi-biWidth;j=2;j-)/浮雕处理buf=*(hData+(lpBi-biHeight-i)*WIDTHBYTES(lpBi-biWidth*8)+j)-*(hData+(lpBi-biHeight-i+1)*WIDTHBYTES(lpBi-biWidth*8)+j-1)+128;if(buf255)buf=255;if(bufbiHeight-i)*WIDTHBYTES(lpBi-biWidth*8)+j)=(BYTE)buf;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);/显示图像/数字图象的马浮雕效果函数/基于模板的高通滤波/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 mask33=-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=0) & (j+t)=0) & (i+s)biHeight) & (j+t)biWidth)sumw += mask1+s1+t;sum+=*(hData+(i+s)*WIDTHBYTES(lpBi-biWidth*8)+(j+t)*mask1+s1+t;if(sumw=0) sumw=1; sum/=sumw;if(sum255)sum=255;if(sumbiWidth*8)+j)=sum;for( j=0; jbiHeight; j+)for(