改进的中值滤波函数

void CSDIELSView::OnMedianFilterImprove()
	{
//程序编制：李立宗
//lilizong@gmail.com
//2012-8-10
		if(myImage1.IsNull())
			OnOpenResourceFile();
		if(!myImage2.IsNull())
			myImage2.Destroy();
		if(myImage2.IsNull()){
			myImage2.Create(myImage1.GetWidth(),myImage1.GetHeight(),24,0);
		}
		//COLORREF pixel; 
		int maxY = myImage1.GetHeight();
		int maxX=myImage1.GetWidth();
		byte* pRealData;
		byte* pRealData2;
		pRealData=(byte*)myImage1.GetBits();
		pRealData2=(byte*)myImage2.GetBits();
		int pit=myImage1.GetPitch();
		int pit2=myImage2.GetPitch();
		//需要注意,pit和pit2的值并不一样,所以如果使用一个值,会导致不同的结果出现
		//CString str;
		//str.Format(TEXT("%d"),pit);
		//MessageBox(str);
		//str.Format(TEXT("%d"),pit2);
		//MessageBox(str);
		int bitCount=myImage1.GetBPP()/8;
		int bitCount2=myImage2.GetBPP()/8;
		int tempR,tempG,tempB;
		int temp,tempX,tempY;
		//int M[3][3]={{1,2,1},{2,4,2},{1,2,1}};
		int M[3][3]={{1,1,1},{1,2,1},{1,1,1}};
		int sum=0;
		int m;
		int pixel[9];
		for(int i=0;i<3;i++)
			for(int j=0;j<3;j++)
				sum=sum+M[i][j];
		//	tempR=tempG=tempG=0;
		//说明：将生产的图像作为24位图处理。
		for (int y=1; y<maxY-1; y++) {
			for (int x=1; x<maxX-1; x++) {

				m=0;
				for(int i=y-1;i<=y+1;i++)
					for(int j=x-1;j<=x+1;j++)
					{
						pixel[m++]= *(pRealData+pit*(i)+(j)*bitCount)  ;
					}

					tempR=getValue(pixel);

					if(bitCount==1)
					{
						tempG=tempR;
						tempB=tempR;
					}
					else
					{				m=0;
					for(int i=y-1;i<=y+1;i++)
						for(int j=x-1;j<=x+1;j++)
						{
							pixel[m++]= *(pRealData+pit*(i)+(j)*bitCount+1)  ;
						}

						tempG=getValue(pixel);
						tempG=pixel[8];
						m=0;
						for(int i=y-1;i<=y+1;i++)
							for(int j=x-1;j<=x+1;j++)
							{
								pixel[m++]= *(pRealData+pit*(i)+(j)*bitCount+2)  ;
							}
							tempB=getValue(pixel);
					}
					*(pRealData2+pit2*y+x*bitCount2)=tempR;
					*(pRealData2+pit2*y+x*bitCount2+1)=tempG;
					*(pRealData2+pit2*y+x*bitCount2+2)=tempB;
			}
		}
		Invalidate();
	}

定义的函数getValue如下：

int getValue(int a[9])
	{
		int i,j,temp;
		int current;
		current=a[4];
		for(i=0;i<8;i++)
			for(j=i+1;j<9;j++) /*注意循环的上下限*/
				if(a[i]>a[j]) {
					temp=a[i];
					a[i]=a[j];
					a[j]=temp;
				}
				if(current==a[0]||current==a[8])
					return a[4];
				else
					return current;
	}