摘要：在眼科疾病檢測中，為了對被檢測者進行快速、準確、自動化的瞳孔定位，提出一種改進徑向對稱變換的瞳孔中心點定位算法。首先利用灰度積分投影法結合最大類間方差法，完成對人眼圖像的粗分割，并根據多團塊篩選條件提取出只包含瞳孔的感興趣區域(RegionOfInterest，ROI)。然后對ROI采用最小外接矩形結合灰度級形態學線性濾波方法，完成搜索半徑范圍的設置。最后，利用改進的徑向對稱變換算法進行瞳孔中心點定位。實驗結果表明：本文算法的定位誤差在8pixel以內，平均定位時間為0.366s，能夠適應人眼圖像中噪聲干擾、采集不完整等大量非理性狀態，滿足多種紅外眼科疾病檢測設備對瞳孔定位算法的要求。

　　關鍵詞：眼科檢查；瞳孔定位；圖像處理；徑向對稱變換

　　1引言

　　目前，全球范圍內視覺健康問題嚴峻[1]：視力損傷或失明人數高達22億人以上，其中至少10億人的視力損傷問題本可預防或尚待解決。有效的眼科疾病檢測對眼部疾病的預防和治療具有重要作用。臨床常用的眼科疾病檢測儀器，如：光學相干層析成像儀(OpticalCoherenceTomo-graphy，OCT)[2-4]、眼底相機[5-6]等在檢測之前，首先需要將儀器成像光軸對準瞳孔中心點，以保證光源發出的光經瞳孔入射到被檢測者眼內。在實際操作中，盡管醫護人員可通過手動調整實現對準，但是手動對準對操作者的熟練度和被檢測者的配合度均有一定的要求，較為繁瑣耗時。因此，自動化人眼對準顯得尤為重要[7]。自動化人眼對準裝置利用瞳孔定位算法獲取瞳孔中心點坐標，進而借助步進電機實現成像光軸與瞳孔中心的對準，具有實時、高效、操作簡單等特點。其中，瞳孔定位算法的性能對定位結影響最大，因此研究快速、準確、不受噪聲和測量條件（如：局部反射光斑、眼瞼眼睫毛遮擋、離焦模糊）影響的瞳孔定位算法具有重要實用價值。

　　2傳統徑向對稱變換算法

　　徑向對稱變換算法是以廣義對稱變換為基礎發展而來的一種基于梯度信息的目標檢測算子。該算法根據被檢測區域的梯度幅值信息和梯度方向信息構造出評價區域對稱性的新指標—徑向對稱度，再采用投票的方式計算被檢測區域中每一點在給定的搜索半徑n處的徑向對稱度。隨著搜索半徑n的遞增，具有徑向對稱特點的目標區域在對稱中心處可以累加到較大的徑向對稱度，從而完成目標區域的檢測。

　　3改進的瞳孔定位算法

　　本文首先利用灰度積分投影融合最大類間方差的方法對人眼圖像進行粗分割，并結合多團塊篩選條件提取出只包含瞳孔的ROI，有效避免了算法由于大范圍運算導致的計算成本高、易受噪聲干擾等問題。然后根據ROI的最小外接矩形和瞳孔區域間的幾何關系，結合灰度級形態學線性濾波，完成了半徑范圍的設置，實現了參數自適應。最后，在設置的搜索半徑內對ROI進行最大徑向對稱度計算，完成瞳孔中心點定位。算法流程如圖1所示。

　　4定位實驗與結果分析

　　為了驗證本文算法的有效性，利用計算機硬件配置為IntelCorei5-6300HQCPU，主頻為2.3GHz，4GB內存，實驗平臺為MATLABR2016b進行實驗。利用實驗室自研的瞳孔自動化對準裝置采集人眼圖像，實驗包括了CCD相機反射像明顯、離焦模糊、眼瞼睫毛遮擋嚴重等情況，圖像大小為752pixel×480pixel，結果如圖6（彩圖見期刊電子版）所示。圖6(a)~6(d)中紅色框標注的區域為ROI區域。瞳孔中心定位結果如圖6(e)~6(h)所示，其中紅色十字（圖中和綠色十字重合）表示人工標定的瞳孔中心點（由專家手動標注，取多次標注的平均值），綠色十字表示算法標定的瞳孔中心點，白色圓圈為算法標定的瞳孔邊界。從提取的ROI圖像來看，灰度積分投影融合最大類間方差方法有效提取出了瞳孔區域，ROI的灰度值范圍為[0,50]，邊界灰度值范圍為[30,50]，滿足感興趣的瞳孔-虹膜邊界灰度值范圍，為后續定位打下了基礎。從定位結果圖像來看，本文算法對存在各類噪聲的人眼圖像的定位結果與人工標定的中心點基本重合，說明本文算法具有較高的定位精度，對噪聲有良好的魯棒性，適用于多種情況下的瞳孔中心點定位。

　　5結論

　　本文根據紅外眼科疾病檢測設備對于眼部對準成像的需求，提出了一種基于徑向對稱變換的瞳孔中心點定位算法。首先利用灰度積分投影結合最大類間方差對人眼圖像進行粗分割，并通過多團塊篩選條件，完成了對人眼圖像的ROI分割。再根據ROI的最小外接矩形和瞳孔區域間的幾何關系，并結合灰度級形態學線性濾波，減小了睫毛遮擋對參數設置的影響，確定了搜索半徑范圍。最后，利用改進的徑向對稱變換算法進行瞳孔中心點定位。

　　參考文獻:

　　[1]FLAXMANSR,BOURNERRA,RESNIKOFFS,etal..Globalcausesofblindnessanddistancevisionimpairment1990-2020:asystematicreviewandmeta-analysis[J].TheLancetGlobalHealth,2017,5(12):e1221-e1234.

　　[2]蔡懷宇，張瑋茜，韓曉袍，等眼軸及眼前節SS-OCT--體化成像系統[J].紅外與激光工程,2020,49(4):0413006.

　　CAIHY,ZHANGWQ,HANXY,etal..IntegratedimagingsystemofeyeaxisandanteriorsegmentbasedonSS-OCT[J].InfraredcandLaserEngineering,2020,49(4):0413006.(inChinese)

　　[3]劉珊珊,汪毅,張瑋茜,等.基于SS_OCT的大范圍眼軸長度測量系統[J].光子學報,2019,48(5):512002.

　　LIUSHSH,WANGY,ZHANGWQ,etal..Large-scaleaxiallengthmeasuringsystembasedonSS-OCT[J].ActaPhotonicaSinica,2019,48(5):512002.(inChinese)

　　蔡懷宇1，史玉1*，婁世良1，汪毅1，陳文光2，陳曉冬