区域增长法（vc实现）

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

void CISLSView::OnLButtonDown(UINT nFlags, CPoint point)
  {
//程序编制：李立宗
//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;
  }