【正文】 sel detection include segmentation technique, tracking techniques and fitting techniques. This thesis arms to study the minimalpath based tracking algorithm and the corresponding implementation techniques. The main contributions are the follows: 1) presented a brief review of various published detection methods of retinal vessels。盡我所知，除文中特別加以標(biāo)注和致謝的地方外，不包含其他人或組織已經(jīng)發(fā)表或公布過的研究成果，也不包含我為獲得 及其它教育機構(gòu)的學(xué)位或?qū)W歷而使用過的材料。臨床上,視網(wǎng)膜血管網(wǎng)絡(luò)對高血壓、糖尿病、動脈硬化、腎炎等疾病的診斷、治療、愈后評價都具有重要意義 [23]。目前視網(wǎng)膜血管的檢測方法已有不少，但是基于 Canny 算子與數(shù)學(xué)形態(tài)學(xué)濾波相結(jié)合的方法算法尚未見發(fā)表。該方法既可有效地抑制圖像中的噪聲，有可保持邊緣的連通性，是一種較理想的邊緣檢測算子。圖象的邊緣部分集中了圖象的大部分信息，圖象邊緣的確定與提取對于整個圖象場景的識別與理解是非常重要的，同時也是圖象浙江省“新苗人才計劃”項目 視網(wǎng)膜血管分割系統(tǒng)設(shè)計與實現(xiàn) 項目研究報告6分割所依賴的重要特征，邊緣檢測主要是圖象的灰度變化的度量、檢測和定位。 基于 Laplacian算子的視網(wǎng)膜圖像的分割方法Laplacian 算子是常用的邊緣檢測算子。圖11（a）表示離散拉普拉斯算子的模板，圖11（b）表示其擴展模板，從模板形式容易看出，如果在圖像中一個較暗的區(qū)域中出現(xiàn)了一個亮點，那么用拉普拉斯運算就會使這個亮點變得更亮。 11(a) 11(b) 基于 Sobel算子的視網(wǎng)膜圖像的分割方法Sobel 算子在圖像上求出二維的空間梯度，和梯度高的區(qū)域與邊緣的區(qū)域相浙江省“新苗人才計劃”項目 視網(wǎng)膜血管分割系統(tǒng)設(shè)計與實現(xiàn) 項目研究報告7一致。而 LOG 算子比 Laplacian 算子多做的就是它是先對圖像進行了 Gauss 平滑。圖像采集主要完成視網(wǎng)膜原始圖像的獲取及保存任務(wù)，然后經(jīng)過預(yù)處理階段進行圖像的去噪、平滑等操作，形成標(biāo)準(zhǔn)視網(wǎng)膜數(shù)據(jù)庫，再對這些數(shù)據(jù)庫中的圖像進行算法處理，然后保存處理以后的圖片以進行不同的算法之間的比較。以一個人采集為例：第一步，將 Nommyd7 免散瞳眼底相機跟計算機連接好。具體命名方法如下：每幅圖像用二個字段表示 ，即表示為 AB. bmp 的形式。這樣命名，對于較小的環(huán)境而言，是有好處的，可以直接通過命名，得知這是誰的視網(wǎng)膜。這也是我們這個系統(tǒng)第一步要做的事情。算法：1) 取得圖像文件的拷貝文件；2) 取得圖像的數(shù)據(jù)區(qū)指針，得到 R，G，B 三個分值；3) 對 R，G，B 三個分量做比較，gray=max（R，G，B） ；4) 將 gray 返回給 R，G，B 三個分量；5) 顯示圖像。并在此基礎(chǔ)上提出了一個實用算法。(3)對同一邊緣僅有唯一的響應(yīng)，即單個邊緣產(chǎn)生多個響應(yīng)的概率要低，并且虛假邊緣應(yīng)受到極大的抑制。一般通過計算梯度幅值完成。通俗的講，高斯濾波就是對整幅圖像進行加權(quán)平均的過程，每一個像素點的值，都由其本身和鄰域內(nèi)的其他像素值經(jīng)過加權(quán)平均后得到。 計算梯度與梯度幅值這一部分是以 F(x,y)為基礎(chǔ)的，要求出 F(x,y)的梯度則必須要求出 F(x,y)的 x 偏導(dǎo)和 y 偏導(dǎo)。X偏導(dǎo)與Y偏導(dǎo)分別為一下兩個。M[i,j]是梯度幅值。其實可以說是4個方向，因為這個鄰域是對稱。在每一個點上, 鄰域的中心像素 m [ i, j ]與沿梯度方向的 2 個梯度幅值的插值結(jié)果進行比較。然后就是拿中間的那個值與8和4這兩個位置的像素進行比較，如果不是最好的則將其設(shè)置為0，不然則這個像素就是邊緣的極值點。 Canny算法的缺點(1) Canny算法在求梯度的時候，由于模板的不同會有很大的不同。膨脹和腐蝕是數(shù)學(xué)形態(tài)學(xué)的2種最基本運算，由這2種基本運算可以組成多種復(fù)合運算，產(chǎn)生各種形態(tài)學(xué)實用算法。(x,y) ∈D } ○ + f b其中，f(x,y)和 b(x,y)表示定義在二維離散空間 D 和 D 上的兩個離散函fb數(shù)，f(x,y)表示輸入圖像，b(x,y)表示結(jié)構(gòu)元素； 腐蝕：是膨脹的對偶運算，表示用結(jié)構(gòu)元素 b 對圖像 f 進行檢測，尋找圖像內(nèi)部可以容納該元素的區(qū)域，若原點在結(jié)構(gòu)元素內(nèi)部，則腐蝕具有收縮圖像作用；若原點不在結(jié)構(gòu)元素內(nèi)部，則可用于填充圖像內(nèi)部的孔洞；同樣，其離散形式的變換如下：(f ?b) (s,t)=min{f(s+x,t+y)b(x,y)︱(s+x,t+y) ∈D 。這個數(shù)學(xué)形態(tài)學(xué)主要做的是前期的工作，而 Canny算法本身應(yīng)該說是不變的。本文中數(shù)學(xué)形態(tài)學(xué)主要用于前期的圖像去噪。因為數(shù)學(xué)形態(tài)學(xué)是基于信息的幾何特征，預(yù)先定義了結(jié)構(gòu)元素，所以能達到既能消除噪聲又能很好的保留圖像的信息。 浙江省“新苗人才計劃”項目 視網(wǎng)膜血管分割系統(tǒng)設(shè)計與實現(xiàn) 項目研究報告16第 3章 實驗結(jié)果 實驗結(jié)果實驗的結(jié)果是通過 5 個人的視網(wǎng)膜進行邊緣識別。雖然從這張照片中看不去數(shù)學(xué)形態(tài)學(xué)的很大的優(yōu)勢，但是如果當(dāng)照片原圖沒有那么好的情況下，則優(yōu)勢就很明顯了。本文主要的研究工作和成果體現(xiàn)在如下幾個方面:1).對 Sobel，Reberts，Laplacian，LOG，Canny 幾種常見算法進行了比較，然后通過對不同質(zhì)量照片的出來，來發(fā)現(xiàn)它們的優(yōu)點和缺點。并對數(shù)學(xué)形態(tài)學(xué)也有了一定的掌握的，數(shù)學(xué)形態(tài)學(xué)是一門新興的、以形態(tài)為基礎(chǔ)對圖像進行分析的學(xué)科,它用具有一定形態(tài)的結(jié)構(gòu)元素去度量和提取圖像中的對應(yīng)形狀以實現(xiàn)對圖像的分析和識別、邊緣檢測等圖像預(yù)處理問題中有著明顯的優(yōu)勢。尤其是在一些核心算法的研究上，如何針對視網(wǎng)膜圖像的特點，提出更高效的邊緣提取算法，是后續(xù)研究工作的重點。作為新興的生物特征身份識別技術(shù)，視網(wǎng)膜識別憑借其高度的準(zhǔn)確性，已經(jīng)得到學(xué)術(shù)界和企業(yè)界越來越廣泛的關(guān)注。private:int m_imgWidthOut。ImgSegment(CSize size,int nBitCount,LPRGBQUAD lpColorTable,unsigned char *pImgData)。void Canny()。//求梯度的函數(shù) //閾值估計void EstimateThreshold(unsigned char * imgIn, int nWidth, int nHeight, int *pnThdHigh,int *pnThdLow, unsigned char * imgOut, double dRatioHigh, double dRationLow)。private://模板卷積函數(shù)void TemplateEdge(unsigned char * imgIn,int width,int height,int nBitCount,int *mask,int maskW,int maskH,unsigned char *amp。G_mask[0] = 2。浙江省“新苗人才計劃”項目 視網(wǎng)膜血管分割系統(tǒng)設(shè)計與實現(xiàn) 項目研究報告23G_mask[4] = 2。G_mask[8] = 9。G_mask[12] = 15。G_mask[16] = 9。G_mask[20] = 2。G_mask[24] = 2。Grad_maskX[1] = 0。Grad_maskX[5] = 1。Grad_maskY[0] = 0。Grad_maskY[4] = 0。浙江省“新苗人才計劃”項目 視網(wǎng)膜血管分割系統(tǒng)設(shè)計與實現(xiàn) 項目研究報告24Grad_maskY[8] = 0。m_lpColorTableOut = NULL。memcpy(thism_lpColorTableOut,m_lpColorTable,sizeof(RGBQUAD) * thism_nColorTableLengthOut)。//臨時變量 unsigned char * T_ImgDataOut_G。T_ImgDataOut_G = new unsigned char[LineByte * thism_imgWidthOut]。浙江省“新苗人才計劃”項目 視網(wǎng)膜血管分割系統(tǒng)設(shè)計與實現(xiàn) 項目研究報告25//thism_lpImgDataOut = new unsigned char[LineByte * thism_imgWidthOut]。//梯度幅值for(int i =0 。//閾值化thisHysteresis(thism_lpImgDataOut,thism_imgWidth,thism_imgHeight,T_ImgDataOut_S)。delete [] T_ImgDataOut_Y。//輸出圖像的顏色表LPRGBQUAD m_lpColorTableOut。int *m_maskBuf。Morphology(CSize size, int nBitCount, LPRGBQUAD lpColorT