区域增长法（vc实现）

通过在图像上单击，获取当前点像素，然后根据当前像素进行像素区间的选择。

void CISLSView::OnLButtonDown(UINT nFlags, CPoint point){
//程序编制：李立宗  [e-mail]
//2012-8-14
if(flag)
{
flag = FALSE;
CISLSDoc *pDoc = GetDocument();
regionGrow(point);
}
//CView::OnLButtonDown(nFlags, point);
}


// 区域增长
void CISLSView::regionGrow(CPoint point)
{
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;

//}
//}
//计算开始
//取得点取点的值
int Tpoint=*(pRealData+pit*(point.y)+(point.x)*bitCount);
int Td;
if(TFlag)
{
Td=myT;
TFlag=FALSE;
}
else
Td=10;
//CString str;
//str.Format(TEXT("%d"),Td);
//AfxMessageBox(str);
for (int y=0; y<maxY; y++) {
for (int x=0; x<maxX; x++) {
temp=*(pRealData3+pit3*(y)+(x)*bitCount3);
if(temp-Tpoint>Td)
temp=255;
else
temp=0;
*(pRealData2+pit2*(y)+(x)*bitCount2)=temp;
*(pRealData2+pit2*(y)+(x)*bitCount2+1)=temp;
*(pRealData2+pit2*(y)+(x)*bitCount2+2)=temp;

}
}
Invalidate();

}

void CISLSView::OnRegionGrowDemo()
{
if(myImage1.IsNull())
OnOpenResourceFile();
AfxMessageBox(TEXT("请在源图像内单击初始点"));
flag=TRUE;
}


void CISLSView::OnRegionGrowCustome()
{
if(myImage1.IsNull())
OnOpenResourceFile();
CMyThreshlod myDlg;
if(myDlg.DoModal () == IDOK)
myT=myDlg.myT;
else
myT=0;
//CString str;
//str.Format(TEXT("%d"),myT);
//AfxMessageBox(str);
TFlag=TRUE;
AfxMessageBox(TEXT("请在源图像内单击初始点"));
flag=TRUE;
}