[發(fā)明專利]一種基于圖像識別的錐細胞密度的計算方法

有效

• 申請?zhí)枺?/td>: 201310714382.X
• 申請日: 2013-12-23
• 公開（公告）號: CN103745257A
• 公開（公告）日: 2014-04-23
• 發(fā)明（設(shè)計）人: 廖娜;王勤美;陳浩;厲以宇;李超宏
• 申請（專利權(quán)）人: 溫州醫(yī)科大學眼視光器械有限公司;溫州醫(yī)科大學眼視光研究院
• 主分類號: G06M11/00
• 分類號: G06M11/00;G06F19/00;A61B3/12;A61B3/14
• 代理公司: 溫州金甌專利事務(wù)所(普通合伙) 33237
• 代理人: 黃肇平

1.一種基于圖像識別的錐細胞密度的計算方法，其特征在于：其包括以下步驟：

通過掃描設(shè)備獲得一段眼底錐細胞的視頻圖像,截取視頻圖像,選擇待操作的視網(wǎng)膜區(qū)域，并將其保存為圖像文件，選擇圖像文件中待計數(shù)錐細胞密度的區(qū)域作為候選區(qū)域；

（2）將候選區(qū)域內(nèi)的像素大小（pixel²）通過I=????????????????????????????????????????????????換算為視網(wǎng)膜實際大小（mm²），其中I為視網(wǎng)膜的實際大小，x為眼軸長度，為角膜頂點到眼第二節(jié)點的距離，為視角；

（3）識別待計數(shù)的候選區(qū)域內(nèi)的錐細胞，并對錐細胞進行數(shù)量統(tǒng)計和計算相鄰錐細胞間距離；

（4）將步驟（3）中得到的總錐細胞數(shù)除以步驟（2）中的到的視網(wǎng)膜實際面積獲得步驟（1）中候選區(qū)域內(nèi)的錐細胞密度。

2.根據(jù)權(quán)利要求1所述的一種基于圖像識別的錐細胞密度的計算方法，其特征在于：所述掃描設(shè)備為眼底自適應共焦激光掃描檢眼鏡。

3.根據(jù)權(quán)利要求1所述的一種基于圖像識別的錐細胞密度的計算方法，其特征在于：在所述步驟（2）中，將眼球看作是“四個共軸球面折射系統(tǒng)”組成的光學系統(tǒng)，包括角膜前后表面，和晶狀體前后表面；利用近軸光學單球面折射和共軸球面折射系統(tǒng)的成像原理計算全眼系統(tǒng)的主點、焦點及節(jié)點，計算實際視網(wǎng)膜大小。

4.根據(jù)權(quán)利要求1所述的一種基于圖像識別的錐細胞密度的計算方法，其特征在于：在所述步驟（2）中，應用Gullstrand-Emsley?No.1模型眼計算視網(wǎng)膜的放大率，?

M=??????????????，并對掃描設(shè)備對視網(wǎng)膜放大率的影響進行校正，其中M為放大率，?P為試戴鏡的度數(shù)，d為鏡眼距，y為角膜頂點到入瞳的距離。

5.根據(jù)權(quán)利要求1所述的一種基于圖像識別的錐細胞密度的計算方法，其特征在于：所述步驟（3）中，包括以下步驟：步驟一、將掃描設(shè)備獲取的錐細胞圖由空間域轉(zhuǎn)換為頻域上的功率譜圖，通過獲得二維序列I＇，其中T1，T2是水平和垂直方向上的采樣周期，是廣義的噪聲；

步驟二、讀取圖像文件后手動選取5個以上相鄰的錐細胞，作為樣本錐細胞利用周期圖法估計錐細胞響應頻率；

步驟三、將步驟二中獲得的響應頻率作為低通濾波的截止頻率對圖像文件進行降噪處理；

步驟四、錐細胞識別，查找臨近區(qū)域最大值的位置，并將圖像轉(zhuǎn)為二進制圖像，最大值的位置即對應錐細胞中心；

步驟五、對二進制圖像膨脹腐蝕操作后標記連通域內(nèi)的質(zhì)心，得到錐細胞中心的橫縱坐標值，將其存入文本；

步驟六、人工修正，手動刪除或添加自動識別中錯誤的錐細胞，在文本中自動去除或增加這些點的坐標；

步驟七、最終得到圖像中總的錐細胞數(shù)。

6.根據(jù)權(quán)利要求5所述的一種基于圖像識別的錐細胞密度的計算方法，其特征在于：所述步驟四采用Matlab二維窗口設(shè)計函數(shù)創(chuàng)建一個近似圓對稱的帶通濾波器，通頻帶由步驟二中取得的錐細胞響應頻率決定，輸出濾波器的頻率響應曲線，，其中f為濾波器的頻率，I＇為二維序列，h是濾波器的脈沖信號。

7.根據(jù)權(quán)利要求6所述的一種基于圖像識別的錐細胞密度的計算方法，其特征在于：通過應用一個窗口刪除損壞的像素，，其中w為窗口。