细菌分割与统计程序_程序讲解20111011.pptVIP

double * count=new double[256]; memset(count,0,256*sizeof(double)); int i,j; for(i = 0; i cyDIB; i++) for(j = 0; j cxDIB; j++) { lpSrc = (unsigned char*)lpDIBBits + lLineBytes * (cyDIB - 1 - i) + j; count[*lpSrc]=count[*lpSrc]+1; } double max=-1; int indexmax; for(i=0;i256;i++) if(count[i]max) { max=count[i]; indexmax=i; } for(i=0;i256;i++) count[i]=count[i]/max; // 创建蓝色画笔对象 CPen* pPenRed = new CPen; pPenRed-CreatePen(PS_SOLID,1,RGB(0,0,155)); dc-SelectObject(pPenRed); for(i=0;i256;i++) { dc-MoveTo (cxDIB+20+i,355); dc-LineTo (cxDIB+20+i,int(355-330*count[i])); } dc-TextOut (cxDIB+20,360,0); dc-TextOut (cxDIB+20+50-1,360,50); dc-TextOut (cxDIB+20+100-2,360,100); dc-TextOut (cxDIB+20+150-2,360,150); dc-TextOut (cxDIB+20+200-2,360,200); dc-TextOut (cxDIB+20+255-2,360,255); ::DeleteObject (font); } 8.二值化函数的编写 建立“细胞计数程序”菜单，下拉菜单建控件 ID_IMG_Threshold2Binary 1图像二值化 在Doc类中对控件“1图像二值化”添加消息处理函数，编写代码。 函数主要完成：统计图像灰度直方图 计算阈值（阈值的选取方法要灵活） 用计算的阈值对图像二值化 void CImagetestDoc::OnIMGThreshold2Binary() { unsigned char *lpSrc; LPSTR lpDIB = (LPSTR) ::GlobalLock((HGLOBAL)m_hDIB); int cxDIB = (int) ::DIBWidth(lpDIB); int cyDIB = (int) ::DIBHeight(lpDIB); LPSTR lpDIBBits=::FindDIBBits (lpDIB); // 计算图像每行的字节数 long lLineBytes = WIDTHBYTES(cxDIB * 8); //首先利用直方图信息，搜索到可以将白色背景与黑色的细菌分开的灰度阈值 int * count1=new int[256]; int * count2=new int[64]; int * count3=new int[64]; memset(count1,0,256*sizeof(int)); memset(count2,0,64*sizeof(int)); memset(count3,0,64*sizeof(int)); int i,j; int threshold; for(i = 0; i cyDIB; i++) { // 每列 for(j = 0; j cxDIB; j++) { // 指向DIB第i行，第j个象素的指针（这里的行为从上到下的） lpSrc = (unsigned char*)lpDIBBits + lLineBytes * (cyDIB - 1 - i) + j; count1[*lpSrc]=count1[*lpSrc]+1; } } //将中间结果写入到文本文件以便后续分析 FILE *fp; fp=fopen(count1.txt,w+); for(i = 0; i 256; i++) fprintf(fp,%d ,count1[i]); fclose(fp); //对count1序列进行分段累积，形成count2 f