維納濾波C++代碼

1、精選優(yōu)質(zhì)文檔-傾情為你奉上Wiener2濾波C+代碼(2011-03-18 10:17:00)標簽：/code written by Kay() 2011/03/17BOOL WienerFilter(WORD *pData, int iWid, int iHgh, int ihoodW, int ihoodH)if (pData = NULL)return FALSE;int iSize = iWid*iHgh;float *pKernel=new floatihoodW*ihoodH;for ( int i = 0 ; i ihoodW*ihoodH ; i + ) pKerneli =

2、(float)1/(ihoodW*ihoodH);float *pMean=new floatiSize;float *pVal =new floatiSize;for ( i = 0 ; i iSize ; i + ) pMeani = (float)(pDatai*pDatai);Filter2(pVal,pMean,iWid,iHgh,pKernel,ihoodW,ihoodH);Filter2(pMean,pData,iWid,iHgh,pKernel,ihoodW,ihoodH);for ( i = 0 ; i iSize ; i + ) pVali = pVali-(pMeani*

3、pMeani);/temp for meandouble dMean=0;for ( i = 0 ; i iSize ; i + )dMean+=(double)pVali/iHgh;dMean=dMean/iWid;for ( i = 0 ; i 0?(pVali-(float)dMean):0;float ftemp = (pDatai-pMeani)*pVali/(pVali+(float)dMean)+pMeani;if (ftemp (GRAYLEVEL-1)ftemp = GRAYLEVEL-1;pDatai = (WORD)ftemp;delete pKernel;delete

4、pMean;delete pVal;return TRUE;void Filter2(float *pRes,WORD *pData, int iWid, int iHgh, float *pKernel, int iKWid, int iKHgh)float ftemp=0;for ( int j = 0 ; j iHgh ; j + )for ( int i = 0 ; i iWid ; i + )ftemp=0;for ( int x = 0 ; x iKWid ; x + )for ( int y = 0 ; y iKHgh ; y + )int ix=(x-iKWid/2+i)0?0

5、:(x-iKWid/2+i);int iy=(y-iKHgh/2+j)iWid) ix=iWid-1;if (iy+1)iHgh) iy=iHgh-1;ftemp+=pDataix+iy*iWid*pKernelx+y*iKWid;pResj*iWid+i = ftemp;void Filter2(float *pRes,float *pData, int iWid, int iHgh, float *pKernel, int iKWid, int iKHgh)float ftemp=0;for ( int j = 0 ; j iHgh ; j + )for ( int i = 0 ; i iWid ; i + )ftemp=0;for ( int x = 0 ; x iKWid ; x + )for ( int y = 0 ; y iKHgh ; y + )int ix=(x-iKWid/2+i)0?0:(x-iKWid/2+i);int iy=(y-iKHgh/2+j)iWid) ix=iWid-1;if (iy+