图像平滑-加权平均

void CSDIELSView::OnSmoothingWeight()
{
//程序编制： 李立宗   lilizong@gmail.com
//2012-8-8
	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;
		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++) {
				tempR=(int)((float)(
					*(pRealData+pit*(y-1)+(x-1)*bitCount)*M[0][0]+*(pRealData+pit*(y-1)+(x)*bitCount)*M[0][1]+*(pRealData+pit*(y-1)+(x+1)*bitCount)*M[0][2]
				+	*(pRealData+pit*(y)+(x-1)*bitCount)*M[1][0]+*(pRealData+pit*(y)+(x)*bitCount)*M[1][1]+*(pRealData+pit*(y)+(x+1)*bitCount)*M[1][2]
				+*(pRealData+pit*(y+1)+(x-1)*bitCount)*M[2][0]+*(pRealData+pit*(y+1)+(x)*bitCount)*M[2][1]+*(pRealData+pit*(y+1)+(x+1)*bitCount)*M[2][2]	
				)/sum+0.5);
				if(bitCount==1)
				{
					tempG=tempR;
					tempB=tempR;
				}
				else
				{
					tempG=(int)((float)(
					*(pRealData+pit*(y-1)+(x-1)*bitCount+1)*M[0][0]+*(pRealData+pit*(y-1)+(x)*bitCount+1)*M[0][1]+*(pRealData+pit*(y-1)+(x+1)*bitCount+1)*M[0][2]
				+	*(pRealData+pit*(y)+(x-1)*bitCount+1)*M[1][0]+*(pRealData+pit*(y)+(x)*bitCount+1)*M[1][1]+*(pRealData+pit*(y)+(x+1)*bitCount+1)*M[1][2]
				+*(pRealData+pit*(y+1)+(x-1)*bitCount+1)*M[2][0]+*(pRealData+pit*(y+1)+(x)*bitCount+1)*M[2][1]+*(pRealData+pit*(y+1)+(x+1)*bitCount+1)*M[2][2]	
				)/sum+0.5);
					tempB=(int)((float)(
					*(pRealData+pit*(y-1)+(x-1)*bitCount+2)*M[0][0]+*(pRealData+pit*(y-1)+(x)*bitCount+2)*M[0][1]+*(pRealData+pit*(y-1)+(x+1)*bitCount+2)*M[0][2]
				+	*(pRealData+pit*(y)+(x-1)*bitCount+2)*M[1][0]+*(pRealData+pit*(y)+(x)*bitCount+2)*M[1][1]+*(pRealData+pit*(y)+(x+1)*bitCount+2)*M[1][2]
				+*(pRealData+pit*(y+1)+(x-1)*bitCount+2)*M[2][0]+*(pRealData+pit*(y+1)+(x)*bitCount+2)*M[2][1]+*(pRealData+pit*(y+1)+(x+1)*bitCount+2)*M[2][2]	
				)/sum+0.5);

				}

				*(pRealData2+pit2*y+x*bitCount2)=tempR;
				*(pRealData2+pit2*y+x*bitCount2+1)=tempG;
				*(pRealData2+pit2*y+x*bitCount2+2)=tempB;
			}
		}
		Invalidate();
}