快速中值濾波算法(編輯修改稿)
2025-07-17 06:39 本頁面
【文章內(nèi)容簡介】 *PI*p/N。 for(k=j。kN1。k+=(1l)) { if (FFTorIFFT) { //若做傅立葉正變換 cosw =cos(w)。 sinw =sin(w)。 } else{ //傅立葉反變換 cosw =cos(w)。 sinw =sin(w)。 } X = Tr[k] + Tr[k+B]*cosw Ti[k+B] * sinw。 Y = Ti[k] + Tr[k+B]*sinw + Ti[k+B] * cosw。 XX = Tr[k] Tr[k+B]*cosw + Ti[k+B] * sinw。 YY = Ti[k] Tr[k+B]*sinw Ti[k+B] * cosw。 Tr[k] = X。 Ti[k] = Y。 Tr[k+B] = XX。 Ti[k+B] = YY。 } } }}void CImageProcessingDoc::OnFft(){// TODO: Add your mand handler code here int i,j。 int ii,jj。 float temp。 float *Tr。 float *Ti。 Tr = new float[nCols]。 Ti = new float[nCols]。 if ( FourierDataR) { delete FourierDataR。 FourierDataR = NULL。 } if ( FourierDataI) { delete FourierDataI。 FourierDataR = NULL。 } FourierDataR = new float[nRows*nCols]。 FourierDataI = new float[nRows*nCols]。 for(i=0。inRows。i++){ for(j=0。jnCols。j++){ //圖像數(shù)據(jù)先給傅立葉變換數(shù)組 FourierDataR[i*nCols+j] = (float) mSourceData[i*nCols+j]。 FourierDataI[i*nCols+j] = 。 } for (i=0。inRows。i++){ //每行進行傅立葉變換 for (j=0。jnCols。j++){ Tr[j] = FourierDataR[i*nCols + j]。 Ti[j] = FourierDataI[i*nCols + j]。 } FFTandIFFT(Tr,Ti,nCols,1)。 for (j=0。jnCols。j++){ FourierDataR[i*nCols + j] = Tr[j]。 FourierDataI[i*nCols + j] = Ti[j]。 } } delete Tr。 delete Ti。 Tr = new float[nRows]。 Ti = new float[nRows]。 for(j=0。jnCols。j++){ //每列進行傅立葉變換 for (i=0。inRows。i++){ Tr[i] = FourierDataR[i*nCols + j]。 Ti[i] = FourierDataI[i*nCols + j]。 } FFTandIFFT(Tr,Ti,nRows,1)。 for (i=0。inRows。i++){ FourierDataR[i*nCols + j] = Tr[i]。 FourierDataI[i*nCols + j] = Ti[i]。 } } for (i=0。inRows。i++){ for (j=0。jnCols。j++){ temp = sqrt(FourierDataR [i*nCols+j]*FourierDataR [i*nCols+j] +FourierDataI [i*nCols+j]*FourierDataI[i*nCols+j] )。 temp /= 100。 if(temp ) temp = 。 ii = nRows 1 (inRows/2?i+nRows/2:inRows/2)。 jj = (jnCols/2)?(j+nCols/2):(jnCols/2)。 //將變換后現(xiàn)實的原點調(diào)整在中心位置 pResultData[3*ii*nCols+3*jj] = (int) temp。 pResultData[3*ii*nCols+3*jj+1] = (int) temp。 pResultData[3*ii*nCols+3*jj+2] = (int) temp。 } //更新顯示 UpdateAllViews(NULL)。 delete FourierDataR。 delete FourierDataI。 FourierDataI = NULL。 FourierDataR = NULL。 return。 }void CImageProcessingDoc::OnSaltPepperNoice(){// TODO: Add your mand handler code here // TODO: Add your mand handler code here int x。 int y。 Salt_Pepper_Noise(mSourceData,nCols,nRows)。 for(y = 0。 y nRows。 y++){ for(x = 0。 x nCols。 x++){ pSourceData[3*y*nCols+3*x] = (unsigned char) mSourceData[y*nCols+x]。 pSourceData[3*y*nCols+3*x+1] = (unsigned char) mSourceData[y*nCols+x]。 pSourceData[3*y*nCols+3*x+2] = (unsigned char) mSourceData[y*nCols+x]。 } UpdateAllViews(NULL)。 }void CImageProcessingDoc::OnRandomNo
高考資料相關(guān)推薦
鄂ICP備17016276號-1