Sobel算子边缘检测（vc实现）

void CISLSView::OnSobel()
 {
//程序编制：李立宗  lilizong@gmail.com
//2012-8-12
  if(myImage1.IsNull())
   OnOpenResourceFile();
  if(!myImage2.IsNull())
   myImage2.Destroy();
  if(myImage2.IsNull()){
   myImage2.Create(myImage1.GetWidth(),myImage1.GetHeight(),24,0);
  }
  float tempArray[9]={-1,-2,-1,0,0,0,1,2,1};
  DoTemplate(&myImage1,&myImage2,3,3,1, 1, tempArray, 3);
  float tempArray1[9]={1,0,-1,2,0,-2,1,0,-1};
  DoTemplate(&myImage1,&myImage3,3,3,1, 1, tempArray1, 3);
  getMaxValue(&myImage2,&myImage3);
  Invalidate();
 }
void CISLSView::DoTemplate(CImage *myImage1,CImage *myImage2,int tempX, int tempY, int tempCenterX, int tempCenterY, float *tempArray, float tempCoef)
 {
  /*CString str;
  str.Format(TEXT("%f"),tempArray[2]);
  MessageBox(str);*/
  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=255,tempG=255,tempB=0;
  float temp;
  // tempR=tempG=tempG=0;
  //说明：将生产的图像作为24位图处理。
  for (int y=tempY; y<maxY-tempY; y++) {
   for (int x=tempX; x<maxX-tempX; x++) {
    temp=0;
    for(int i=0;i<tempX;i++)
     for(int j=0;j<tempY;j++)
      temp=temp+(*(pRealData+pit*(y-tempCenterY+i)+(x-tempCenterX+j)*bitCount))*tempArray[i*tempX+j];
    temp=tempCoef*temp;
    temp=fabs(temp);
    if(temp>255)
     tempR=255;
    else
     tempR=(int)(temp+0.5);
    if(bitCount==1)
    {
     tempG=tempR;
     tempB=tempR;
    }
    else
    {
     temp=0;
     for(int i=0;i<tempX;i++)
      for(int j=0;j<tempY;j++)
       temp=temp+(*(pRealData+pit*(y-tempCenterY+i)+(x-tempCenterX+j)*bitCount+1))*tempArray[i*tempX+j];
     temp=tempCoef*temp;
     temp=fabs(temp);
     if(temp>255)
      tempG=255;
     else
      tempG=(int)(temp+0.5);
     temp=0;
     for(int i=0;i<tempX;i++)
      for(int j=0;j<tempY;j++)
       temp=temp+(*(pRealData+pit*(y-tempCenterY+i)+(x-tempCenterX+j)*bitCount+2))*tempArray[i*tempX+j];
     temp=tempCoef*temp;
     temp=fabs(temp);
     if(temp>255)
      tempB=255;
     else
      tempB=(int)(temp+0.5);
    }
    *(pRealData2+pit2*y+x*bitCount2)=tempR;
    *(pRealData2+pit2*y+x*bitCount2+1)=tempG;
    *(pRealData2+pit2*y+x*bitCount2+2)=tempB;
   }
  }
  // Invalidate();
 }
 void CISLSView::getMaxValue(CImage *myImage1, CImage *myImage2)
 {
  //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=255,tempG=255,tempB=0;
  int tempR1,tempR2,tempG1,tempG2,tempB1,tempB2;
  float temp;
  // tempR=tempG=tempG=0;
  //说明：将生产的图像作为24位图处理。
  for (int y=0; y<maxY; y++) {
   for (int x=0; x<maxX; x++) {
    tempR1=*(pRealData+pit*(y)+(x)*bitCount);
    tempR2=*(pRealData2+pit2*y+x*bitCount2);
    if(tempR1<tempR2)
     tempR1=tempR2;
    if(bitCount==1)
    {
     tempG=tempR;
     tempB=tempR;
    }
    else
    {
     tempG1=*(pRealData+pit*(y)+(x)*bitCount+1);
     tempG2=*(pRealData2+pit2*y+x*bitCount2+1);
     if(tempG1<tempG2)
      tempG1=tempG2;
     tempB1=*(pRealData+pit*(y)+(x)*bitCount+2);
     tempB2=*(pRealData2+pit2*y+x*bitCount2+2);
     if(tempB1<tempB2)
      tempB1=tempB2;
    }
    *(pRealData+pit*(y)+(x)*bitCount)=tempR1;
    *(pRealData+pit*(y)+(x)*bitCount+1)=tempG1;
    *(pRealData+pit*(y)+(x)*bitCount+2)=tempB1;
   }
  }
 }