void GrayTrans::Binary()
{
long i,j;
unsigned char pixel;
long lHistogram[256];//直方图数组
//定义阈值，最大最小灰度值，平均阈值
unsigned char iThreshold,iNewThreshold,iMaxGrayValue,iMinGrayValue,
iMean1GrayValue,iMean2GrayValue;
long lP1,lP2,lS1,lS2;
int iIterationTimes;//迭代次数
if (m_nBitCount==8)
{
if (m_pImgDataOut!=NULL)
{
delete []m_pImgDataOut;
m_pImgDataOut=NULL;
}
if (m_lpColorTableOut!=NULL)
{
delete []m_lpColorTableOut;
m_lpColorTableOut=NULL;
}
m_nBitCountOut=m_nBitCount;
m_nColorTableLengthOut=ComputeColorTableLength(m_nBitCountOut);
m_lpColorTableOut=new RGBQUAD[m_nColorTableLengthOut];
memcpy(m_lpColorTableOut, m_lpColorTable,
sizeof(RGBQUAD)*m_nColorTableLengthOut);
m_imgWidthOut=m_imgWidth;
m_imgHeightOut=m_imgHeight;
int lineByte=(m_imgWidth*m_nBitCount/8+3)/4*4;//每行像素占字节数
m_pImgDataOut=new unsigned char[lineByte*m_imgHeight]; //为输出图像数据开辟空间
memset(m_pImgDataOut, (BYTE)255, m_imgWidth* m_imgHeight);
for (i = 0; i < 256;i++)
{
lHistogram[i]=0;
}
iMaxGrayValue = 0;
iMinGrayValue = 255;
//求直方图数组，求出最大，最小灰度值
for (i = 0;i < m_imgWidth ;i++)
{
for(j = 0;j pixel)
{
iMinGrayValue = pixel;
}
if(iMaxGrayValue < pixel)
{
iMaxGrayValue = pixel;
}
}
}
iNewThreshold = (iMinGrayValue + iMaxGrayValue)/2;
iThreshold = 0;
//迭代求出阈值
for(iIterationTimes = 0; iThreshold != iNewThreshold &&
iIterationTimes < 100;iIterationTimes ++)
{
iThreshold = iNewThreshold;
lP1 =0;
lP2 =0;
lS1 = 0;
lS2 = 0;
for (i = iMinGrayValue;i < iThreshold;i++)
{
lP1 += lHistogram[i]*i;
lS1 += lHistogram[i];
}
iMean1GrayValue = (unsigned char)(lP1 / lS1);
for (i = iThreshold+1;i < iMaxGrayValue;i++)
{
lP2 += lHistogram[i]*i;
lS2 += lHistogram[i];
}
iMean2GrayValue = (unsigned char)(lP2 / lS2);
iNewThreshold =  (iMean1GrayValue + iMean2GrayValue)/2;
}
//二值化
for (int i=0; i