图像处理：形状特征提取算法

void CMy2011302590043View::OnShapeFeatures() 
{
// TODO: Add your command handler code here
CMy2011302590043Doc* pDoc = GetDocument();
    ASSERT_VALID(pDoc);
DWORD m_dwSizeImage = pDoc->m_Bmpfile.m_Rows*pDoc->m_Bmpfile.m_Cols*pDoc->m_Bmpfile.m_PxlBytes;
LPBYTE m_lpImageCopy = new BYTE[m_dwSizeImage];
DWORD nWidth=pDoc->m_Bmpfile.m_Cols;
DWORD nHeight=pDoc->m_Bmpfile.m_Rows;
WORD wBitCount=pDoc->m_Bmpfile.m_PxlBytes*8;
DWORD lRowBytes=nWidth*((DWORD)(wBitCount/8));
LPBYTE lpData=pDoc->m_Bmpfile.m_pImgDat;
LPBYTE lpOldBits=m_lpImageCopy;


int stop=0;
int counter=0;//counter用于计算连通数
int i,j,t,present;
DWORD m,n;
BYTE *p_temp;
p_temp=new BYTE[lRowBytes*nHeight];
memset(p_temp,255,lRowBytes*nHeight);

const int T=50;
if (wBitCount == 8)
{

//让第一行和第一列像素值为255，
for(i=0;i<nWidth;i++) *(lpData+(nHeight-1)*lRowBytes+i)=255;
for(j=0;j<nHeight;j++) *(lpData+(nHeight-1-j)*lRowBytes)=255;

   
    for(j=1;j<nHeight-1;j++)
   {
if(stop==1)
  break;
for(i=1;i<nWidth-1;i++)
{
      if(counter>255)
       {
       AfxMessageBox("连通区域数目太多，请减少样本个数");//
       stop=1;
       return;
       }
      if(*(lpData+(nHeight-j-1)*lRowBytes+i)<T)//若第（ij）个像素=0！！！！位图扫描顺序从下到上，从左到右
       {
       if(*(lpData+(nHeight-j-1+1)*lRowBytes+i+1)<T)//若右上像素=0
       {
       *(p_temp+(nHeight-j-1)*lRowBytes+i)=*(p_temp+(nHeight-j-1+1)*lRowBytes+i+1);//将右上像素标号赋予当前像素
       present=*(p_temp+(nHeight-j-1+1)*lRowBytes+i+1);//且将右上像素标号赋给present
      
       if(*(lpData+(nHeight-j-1)*lRowBytes+i-1)<T&&*(p_temp+(nHeight-j-1)*lRowBytes+i-1)!=present)//若左像素=0且标号不等于右上
       {
       int temp=*(p_temp+(nHeight-j-1)*lRowBytes+i-1);//即有一个等价冲突，将左像素的连通标号赋值给temp
       if(temp<present)//如果左像素标号小于右上像素
       {
       present=temp;//将左的标号赋给present
       temp=*(p_temp+(nHeight-j-1+1)*lRowBytes+i+1);//将右上标号赋给temp
       }
         counter--;
         for(m=1;m<nHeight;m++)//这两个for循环的目的是将之前和右上一个连通区域的标号都改为左像素的标号，并将大于右上的标号都减一