【正文】 9] = 4。G_mask[20] = 2。G_mask[21] = 4。G_mask[22] = 5。G_mask[23] = 4。G_mask[24] = 2。int Grad_maskX[9]。int Grad_maskY[9]。Grad_maskX[0] = 0。Grad_maskX[1] = 0。Grad_maskX[2] = 0。Grad_maskX[3] = 1。Grad_maskX[4] = 0。Grad_maskX[5] = 1。Grad_maskX[6] = 0。Grad_maskX[7] = 0。Grad_maskX[8] = 0。Grad_maskY[0] = 0。Grad_maskY[1] = 1。Grad_maskY[2] = 0。Grad_maskY[3] = 0。Grad_maskY[4] = 0。Grad_maskY[5] = 0。Grad_maskY[6] = 0。Grad_maskY[7] = 1。浙江省“新苗人才計劃”項目 視網(wǎng)膜血管分割系統(tǒng)設(shè)計與實現(xiàn) 項目研究報告24Grad_maskY[8] = 0。if(thism_lpImgDataOut != NULL){delete []m_lpImgDataOut。m_lpImgDataOut = NULL。}if(thism_lpColorTableOut != NULL){delete []m_lpColorTableOut。m_lpColorTableOut = NULL。}thism_nBitCountOut = m_nBitcount。thism_nColorTableLengthOut = thisComputeColorTableLength(thism_nBitCountOut)。if(m_nColorTableLengthOut != 0){thism_lpColorTableOut = new RGBQUAD[thism_nColorTableLengthOut]。memcpy(thism_lpColorTableOut,m_lpColorTable,sizeof(RGBQUAD) * thism_nColorTableLengthOut)。}thism_imgHeightOut = thism_imgHeight。thism_imgWidthOut = thism_imgWidth。int LineByte = (thism_imgWidth * thism_nBitcount/8 + 3 ) / 4 * 4。//臨時變量 unsigned char * T_ImgDataOut_G。unsigned char * T_ImgDataOut_X。unsigned char * T_ImgDataOut_Y。unsigned char * T_ImgDataOut_S。T_ImgDataOut_G = new unsigned char[LineByte * thism_imgWidthOut]。T_ImgDataOut_X = new unsigned char[LineByte * thism_imgWidthOut]。T_ImgDataOut_Y = new unsigned char[LineByte * thism_imgWidthOut]。T_ImgDataOut_S = new unsigned char[LineByte * thism_imgWidthOut]。浙江省“新苗人才計劃”項目 視網(wǎng)膜血管分割系統(tǒng)設(shè)計與實現(xiàn) 項目研究報告25//thism_lpImgDataOut = new unsigned char[LineByte * thism_imgWidthOut]。//Gauss 平滑thisTemplateEdge(m_pImgData,thism_imgWidth,thism_imgHeight,thism_nBitcount,G_mask,5,5,T_ImgDataOut_G,false)。//求 X的梯度thisTemplateEdge(T_ImgDataOut_G,thism_imgWidth,thism_imgHeight,thism_nBitcount,Grad_maskX,3,3,T_ImgDataOut_X,true)。//求 Y的梯度thisTemplateEdge(T_ImgDataOut_G,thism_imgWidth,thism_imgHeight,thism_nBitcount,Grad_maskY,3,3,T_ImgDataOut_Y,true)。//梯度幅值for(int i =0 。 i LineByte * thism_imgWidthOut。 i++){thism_lpImgDataOut[i] = thisComputeGrad(T_ImgDataOut_X[i],T_ImgDataOut_Y[i])。}//非極大值限制thisNoMaxSuppress(thism_lpImgDataOut,T_ImgDataOut_X,T_ImgDataOut_Y,thism_imgWidth,thism_imgHeight,T_ImgDataOut_S)。//閾值化thisHysteresis(thism_lpImgDataOut,thism_imgWidth,thism_imgHeight,T_ImgDataOut_S)。memcpy(thism_lpImgDataOut,T_ImgDataOut_S ,LineByte * thism_imgWidthOut)。delete [] T_ImgDataOut_G。delete [] T_ImgDataOut_X。delete [] T_ImgDataOut_Y。浙江省“新苗人才計劃”項目 視網(wǎng)膜血管分割系統(tǒng)設(shè)計與實現(xiàn) 項目研究報告26delete [] T_ImgDataOut_S。}浙江省“新苗人才計劃”項目 視網(wǎng)膜血管分割系統(tǒng)設(shè)計與實現(xiàn) 項目研究報告27附錄 2 數(shù)學形態(tài)學程序源代碼：class Morphology : public ImgCenterDib {public://輸出圖像每像素的位數(shù)int m_nBitCountOut。//輸出圖像位圖數(shù)據(jù)指針unsigned char * m_lpImgDataOut。//輸出圖像的顏色表LPRGBQUAD m_lpColorTableOut。private:int m_imgWidthOut。int m_imgHeightOut。int m_nColorTableLengthOut。int *m_maskBuf。int m_maskW。int m_maskH。public:Morphology()。Morphology(CSize size, int nBitCount, LPRGBQUAD lpColorTable, unsigned char *pImgData)。CSize GetDimensions()。//灰值腐蝕浙江省“新苗人才計劃”項目 視網(wǎng)膜血管分割系統(tǒng)設(shè)計與實現(xiàn) 項目研究報告28void GrayErosion()。//灰值膨脹void GrayDilation()。//灰值開void GrayOpen()。//灰值閉void GrayClose()。//灰值形態(tài)學梯度void GrayGradient()。 ~Morphology()。//結(jié)構(gòu)元素（模板）輸入函數(shù)，在調(diào)用其它形態(tài)學處理函數(shù)前必須先調(diào)用此函數(shù)void InputMask(int *mask,int maskW, int maskH)。private://灰值腐蝕基本運算 ,灰值開、閉、形態(tài)學梯度等操作都要調(diào)用這個函數(shù)void BasicErosionForGray(unsigned char *imgBufIn, unsigned char *imgBufOut,int imgWidth,int imgHeight,int *maskBuf, int maskW, int maskH)。//灰值膨脹基本運算 ,灰值開、閉、形態(tài)學梯度操作等都要調(diào)用這個函數(shù)void BasicDilationForGray(unsigned char *imgBufIn, unsigned char *imgBufOut,int imgWidth,int imgHeight,int *maskBuf, int maskW, int maskH)。}。浙江省“新苗人才計劃”項目 視網(wǎng)膜血管分割系統(tǒng)設(shè)計與實現(xiàn) 項目研究報告29致 謝衷心的感謝張永平教授給我們耐心的指導和幫助，張教授嚴謹?shù)闹螌W態(tài)度、豐富的研究經(jīng)驗、開闊的思維方式都給我們以很大的啟迪和教益。這不僅豐富了我們的知識而且更開闊了我們的視野。從項目的選題、研究過程的逐步推進、問題的解決到論文的撰寫、修改等全過程都是在張教授的悉心指導下完成的。在此期間，老師給予了我們很多關(guān)心、寬容和支持，同時，感謝蘇老師和所有幫助過我們的教師，正是在老師們的幫助下，才使得課題順利的完成。在此向他們表示由衷的謝意。浙江省“新苗人才計劃”項目 視網(wǎng)膜血管分割系統(tǒng)設(shè)計與實現(xiàn) 項目研究報告30數(shù)據(jù)庫圖像