二值图像腐蚀（vc实现）

为了处理上的方便，将彩色图像首先进行了灰度化、二值化处理。

void CISLSView::OnErosionBin()
{
//程序编制：李立宗  lilizong@gmail.com
//2012-8-17
//图像腐蚀，二值图像腐蚀。
	//int structure[3][3]={0,1,0,1,1,1,0,1,0};
	int structure[3][3]={1,1,1,1,1,1,1,1,1};
	erosionBin(structure);
}

void CISLSView::erosionBin(int structure[3][3])
{
	if(myImage1.IsNull())
		OnOpenResourceFile();
	if(!myImage2.IsNull())
		myImage2.Destroy();
	if(myImage2.IsNull()){
		myImage2.Create(myImage1.GetWidth(),myImage1.GetHeight(),24,0);
	}
	if(myImage3.IsNull()){
		myImage3.Create(myImage1.GetWidth(),myImage1.GetHeight(),24,0);
	}
	//COLORREF pixel; 
	int maxY = myImage1.GetHeight();
	int maxX=myImage1.GetWidth();
	byte* pRealData;
	byte* pRealData2;
	byte* pRealData3;
	pRealData=(byte*)myImage1.GetBits();
	pRealData2=(byte*)myImage2.GetBits();
	pRealData3=(byte*)myImage3.GetBits();
	int pit=myImage1.GetPitch();
	int pit2=myImage2.GetPitch();
	int pit3=myImage3.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 bitCount3=myImage3.GetBPP()/8;
	int tempR,tempG,tempB;
	//float temp,tempX,tempY;
	int temp;
	float u0,u1;   //均值
	float w0,w1;   //概率
	float sum0,sum1;  //像素和
	int optIndex,optT;   //最优阈值，及其所在像素的值
	float fVaria,fMaxVaria=0;   //临时方差,最大方差

	//int pixelR[256],pixelG[256],pixelB[256];
	int pixel[256]={0};   //不要忘记初始化
	//灰度化
	for (int y=0; y<maxY; y++) {
		for (int x=0; x<maxX; x++) {
			temp=*(pRealData+pit*(y)+(x)*bitCount);
			if(bitCount==3)
			{
				tempR=*(pRealData+pit*(y)+(x)*bitCount);
				tempG=*(pRealData+pit*(y)+(x)*bitCount+1);
				tempB=*(pRealData+pit*(y)+(x)*bitCount+2);
				temp=(int)(tempR*0.49+tempG*0.31+tempB*0.2);
				//temp=(int)((tempR+tempG+tempB)/3);
			}
			*(pRealData3+pit3*(y)+(x)*bitCount3)=temp;
			*(pRealData3+pit3*(y)+(x)*bitCount3+1)=temp;
			*(pRealData3+pit3*(y)+(x)*bitCount3+2)=temp;

		}
	}
	//二值化
	for (int y=0; y<maxY; y++) {
		for (int x=0; x<maxX; x++) {
			temp=*(pRealData3+pit3*(y)+(x)*bitCount3);
			if(temp>128)
				temp=255;
			else
				temp=0;
			*(pRealData3+pit3*(y)+(x)*bitCount3)=temp;
			*(pRealData3+pit3*(y)+(x)*bitCount3+1)=temp;
			*(pRealData3+pit3*(y)+(x)*bitCount3+2)=temp;

		}
	}	
	//CString str;
	//str.Format(TEXT("%d"),Td);
	//AfxMessageBox(str);
	//
	for (int y=0; y<maxY; y++) {
		for (int x=0; x<maxX; x++) {

			*(pRealData2+pit2*(y)+(x)*bitCount2)=255;
			*(pRealData2+pit2*(y)+(x)*bitCount2+1)=255;
			*(pRealData2+pit2*(y)+(x)*bitCount2+2)=255;

		}
	}
	int m,n;
	for (int y=1; y<maxY-1; y++) {
		for (int x=1; x<maxX-1; x++) {
			temp=*(pRealData3+pit3*(y)+(x)*bitCount3);

			for(m=0;m<3;m++)
				for(n=0;n<3;n++)
				{	if(structure[m][n]==0)
				continue;
			temp=*(pRealData3+pit3*(y-m+1)+(x+n-1)*bitCount3);
			if(temp==0)
			{
				*(pRealData2+pit2*(y)+(x)*bitCount2)=0;
				*(pRealData2+pit2*(y)+(x)*bitCount2+1)=0;
				*(pRealData2+pit2*(y)+(x)*bitCount2+2)=0;
				break;
			}

				}
		}
	}
	Invalidate();
}