簡介這份資源聚焦虹膜分割與虹膜提取的工程實(shí)現(xiàn)面向具備一定圖像處理基礎(chǔ)、希望用C與OpenCV落地生物特征識(shí)別項(xiàng)目的開發(fā)者與學(xué)習(xí)者。包內(nèi)提供可直接編譯運(yùn)行的Visual Studio工程核心代碼以cpp源文件形式給出并配有測(cè)試圖片與一份方法講解PPT便于對(duì)照理解分割流程與算法思路。資源共39個(gè)文件包含10張bmp與10張jpg測(cè)試圖像、2個(gè)cpp實(shí)現(xiàn)文件以及sln、vcproj、exe、pdb等工程與調(diào)試文件另有pdf說明文檔壓縮包約3.26MB結(jié)構(gòu)完整、開箱即用。目前已有1124人學(xué)習(xí)下載適合作為課程設(shè)計(jì)、畢業(yè)設(shè)計(jì)或虹膜識(shí)別入門階段的參考實(shí)現(xiàn)讀者可借此快速搭建實(shí)驗(yàn)環(huán)境觀察分割效果并在此基礎(chǔ)上做算法改進(jìn)與二次開發(fā)。1. 虹膜分割到底在分什么從一張人眼圖到可用的掩膜很多人第一次拿到虹膜分割任務(wù)以為就是把瞳孔那一圈摳出來。真上手才發(fā)現(xiàn)一張人眼圖里要分的東西至少有三層瞳孔、虹膜環(huán)帶、以及把這兩者包住的眼瞼和睫毛干擾區(qū)。虹膜分割代碼要輸出的不是一張好看的圖而是一張二值掩膜告訴后續(xù)特征提取模塊「哪些像素屬于虹膜、哪些是噪聲」。這份 C 和 OpenCV 實(shí)現(xiàn)的資源核心價(jià)值就在于把這條鏈路用原生代碼跑通不依賴深度學(xué)習(xí)框架編譯完就能在本地對(duì)單張圖或批量圖做虹膜提取。它適合兩類人一類是做生物特征識(shí)別、想理解傳統(tǒng)圖像處理流程的從業(yè)者另一類是需要把虹膜定位嵌進(jìn) C 工程、又不想引入龐大推理依賴的開發(fā)者。整份代碼圍繞灰度化、邊緣檢測(cè)、霍夫圓變換、掩膜生成幾個(gè)環(huán)節(jié)展開參數(shù)都暴露在源碼里改起來直觀。下面按「先跑通、再調(diào)參、最后避坑」的順序拆。2. 環(huán)境搭建與第一張掩膜把 C 工程跑起來2.1 依賴選型為什么是 OpenCV 而不是別的虹膜分割的經(jīng)典路線是「邊緣 圓擬合」這條路線上 OpenCV 提供的HoughCircles、Canny、GaussianBlur幾乎覆蓋了全部需求。選 OpenCV 而不是自己寫卷積理由很實(shí)際霍夫圓變換內(nèi)部涉及累加器空間和半徑搜索手寫一版能跑但調(diào)參成本極高而 OpenCV 的實(shí)現(xiàn)已經(jīng)做了非極大值抑制和梯度篩選參數(shù)語義清晰。C 版本相比 Python 的優(yōu)勢(shì)在于批量處理上千張圖時(shí)沒有解釋器開銷且能直接編進(jìn)已有的桌面或嵌入式工程。常見做法是用 CMake 管理工程把 OpenCV 作為外部依賴鏈接進(jìn)來。我一般會(huì)先確認(rèn)本機(jī) OpenCV 版本因?yàn)?4.x 和 3.x 在HoughCircles的參數(shù)默認(rèn)值上有細(xì)微差別尤其是param2的敏感度。# 查看本機(jī) OpenCV 版本決定后續(xù)鏈接的庫名 pkg-config --modversion opencv4 # 若輸出 4.x則 CMake 里用 find_package(OpenCV 4 REQUIRED)這段命令的作用是確認(rèn)版本號(hào)。如果輸出為空說明pkg-config路徑?jīng)]配好需要手動(dòng)指定OpenCV_DIR。參數(shù)上沒什么可調(diào)的純粹是環(huán)境探測(cè)。2.2 編譯與運(yùn)行從源碼到可執(zhí)行文件工程結(jié)構(gòu)通常是src/放源碼、include/放頭文件、CMakeLists.txt在根目錄。下面是一個(gè)最小可用的 CMake 配置鏈接 OpenCV 的核心和圖像處理模塊。cmake_minimum_required(VERSION 3.10) project(IrisSegmentation) # 只鏈接需要的模塊減少編譯體積 find_package(OpenCV REQUIRED COMPONENTS core imgproc highgui) add_executable(iris_seg src/main.cpp src/iris_segment.cpp) target_include_directories(iris_seg PRIVATE ${OpenCV_INCLUDE_DIRS} include) target_link_libraries(iris_seg PRIVATE ${OpenCV_LIBS})邏輯說明COMPONENTS里只寫core、imgproc、highgui是因?yàn)楹缒し指钣貌坏絛nn、videoio這些模塊少鏈接能加快編譯。target_include_directories把自定義頭文件目錄加進(jìn)去避免源碼里寫一長串相對(duì)路徑。參數(shù)上REQUIRED保證找不到 OpenCV 時(shí)直接報(bào)錯(cuò)而不是靜默跳過。編譯命令就三步mkdir build cd build cmake .. make -j4-j4是按四核并行編譯機(jī)器核多可以改大。編譯完成后build/下會(huì)生成iris_seg可執(zhí)行文件運(yùn)行時(shí)傳入圖片路徑即可。2.3 第一張掩膜的生成流程主流程分四步讀圖轉(zhuǎn)灰度、高斯模糊降噪、霍夫圓檢測(cè)定位瞳孔和虹膜外邊界、根據(jù)兩個(gè)圓生成環(huán)帶掩膜。下面這段是核心函數(shù)的骨架。// 輸入原始 BGR 圖輸出虹膜環(huán)帶掩膜 cv::Mat segmentIris(const cv::Mat src) { cv::Mat gray, blurred; cv::cvtColor(src, gray, cv::COLOR_BGR2GRAY); // 高斯核 5x5sigma 由 OpenCV 自動(dòng)推算先壓掉睫毛高頻噪聲 cv::GaussianBlur(gray, blurred, cv::Size(5, 5), 0); std::vectorcv::Vec3f circles; // param1 是 Canny 高閾值param2 是累加器閾值越小越容易檢出假圓 cv::HoughCircles(blurred, circles, cv::HOUGH_GRADIENT, 1, gray.rows / 8, 200, 50, 20, 200); cv::Mat mask cv::Mat::zeros(gray.size(), CV_8UC1); if (circles.size() 2) { // 半徑大的當(dāng)虹膜外邊界小的當(dāng)瞳孔邊界 cv::Vec3f pupil circles[0], iris circles[1]; if (pupil[2] iris[2]) std::swap(pupil, iris); cv::circle(mask, cv::Point(pupil[0], pupil[1]), pupil[2], 0, -1); // 挖掉瞳孔 cv::circle(mask, cv::Point(iris[0], iris[1]), iris[2], 255, -1); // 填充虹膜 } return mask; }邏輯說明先轉(zhuǎn)灰度是因?yàn)榛舴驁A只吃單通道。高斯模糊的核大小 5x5 是經(jīng)驗(yàn)值太大虹膜邊緣會(huì)糊太小睫毛噪聲壓不住。HoughCircles里dp1表示累加器分辨率與原圖一致minDist設(shè)為行高的八分之一防止同一個(gè)圓被檢出多次。param1200是 Canny 高閾值param250是累加器閾值這兩個(gè)是調(diào)參重點(diǎn)。最后用「先填白再挖黑」的方式生成環(huán)帶比逐像素判斷半徑更省事。參數(shù)怎么改如果檢不出圓先把param2降到 30 左右如果檢出太多假圓把param2提到 70 以上。minDist太小會(huì)導(dǎo)致瞳孔和虹膜被當(dāng)成同一個(gè)圓太大則可能漏掉半徑接近的雙圓。3. 參數(shù)調(diào)優(yōu)與魯棒性讓掩膜在臟圖上也能用3.1 霍夫圓參數(shù)的實(shí)際影響HoughCircles的五個(gè)關(guān)鍵參數(shù)里dp、minDist、param1、param2、minRadius/maxRadius每一個(gè)都會(huì)顯著改變結(jié)果。下面這張表是我在幾十張不同質(zhì)量人眼圖上試出來的經(jīng)驗(yàn)區(qū)間不是理論值是踩坑踩出來的。參數(shù)作用保守取值激進(jìn)取值翻車表現(xiàn)dp累加器分辨率12取 2 時(shí)小瞳孔直接漏檢minDist圓心最小間距rows/8rows/4太小雙圓合并太大漏檢param1Canny 高閾值200100太低邊緣碎圓擬合發(fā)散param2累加器閾值5030太低滿屏假圓minRadius最小半徑2010太小把噪點(diǎn)當(dāng)瞳孔調(diào)參順序建議先固定dp1、minDistrows/8然后調(diào)param2讓圓的數(shù)量穩(wěn)定在 2 到 4 個(gè)最后用minRadius/maxRadius把明顯不合理的圓篩掉。param1一般不動(dòng)除非圖像對(duì)比度特別低。3.2 預(yù)處理對(duì)分割結(jié)果的影響原始人眼圖里睫毛、反光、眼瞼陰影都是干擾源。直接上霍夫圓十張里能有三張檢出正確雙圓就不錯(cuò)了。常見做法是在灰度化之后加一步形態(tài)學(xué)開運(yùn)算把細(xì)長的睫毛去掉同時(shí)保留虹膜邊緣。cv::Mat denoise(const cv::Mat gray) { cv::Mat blurred, opened; cv::GaussianBlur(gray, blurred, cv::Size(7, 7), 0); // 3x3 矩形核做開運(yùn)算先腐蝕掉睫毛再膨脹回來 cv::Mat kernel cv::getStructuringElement(cv::MORPH_RECT, cv::Size(3, 3)); cv::morphologyEx(blurred, opened, cv::MORPH_OPEN, kernel); return opened; }邏輯說明核從 5x5 提到 7x7 是為了讓睫毛這種細(xì)結(jié)構(gòu)在模糊階段就衰減掉。開運(yùn)算的核不能大3x3 足夠再大虹膜外邊緣會(huì)被啃掉一圈導(dǎo)致半徑偏小。這一步的代價(jià)是瞳孔邊緣也會(huì)略微變鈍所以后續(xù)param2要相應(yīng)降低 5 到 10。另一個(gè)常見干擾是紅外圖里的高光斑。高光在灰度圖上表現(xiàn)為一塊高亮區(qū)域霍夫圓會(huì)把它的邊緣當(dāng)成圓弧。處理辦法是在開運(yùn)算之后加一步閾值截?cái)喟殉^ 240 的像素壓到 240削弱高光梯度。3.3 雙圓匹配邏輯的邊界處理HoughCircles返回的圓是無序的不能默認(rèn)第一個(gè)是瞳孔、第二個(gè)是虹膜。除了按半徑排序還要加兩個(gè)約束圓心距離不能太遠(yuǎn)、半徑比要在合理范圍。正常正視人眼里瞳孔半徑約為虹膜半徑的 0.3 到 0.6圓心偏移不超過虹膜半徑的 0.2。bool isValidPair(const cv::Vec3f a, const cv::Vec3f b) { float rMin std::min(a[2], b[2]); float rMax std::max(a[2], b[2]); float ratio rMin / rMax; float dx a[0] - b[0], dy a[1] - b[1]; float dist std::sqrt(dx * dx dy * dy); // 半徑比 0.3~0.6圓心偏移不超過大圓半徑的 0.2 return ratio 0.3f ratio 0.6f dist rMax * 0.2f; }邏輯說明這個(gè)函數(shù)在候選圓列表里兩兩配對(duì)只保留滿足幾何約束的組合。如果沒有任何一對(duì)通過說明這張圖質(zhì)量太差應(yīng)該跳過而不是硬分割。參數(shù) 0.3 和 0.6 是統(tǒng)計(jì)值遇到側(cè)視或斜視人眼可以放寬到 0.25 和 0.7但放寬后假陽性會(huì)上升需要配合后續(xù)掩膜面積過濾。4. 避坑與排查五個(gè)真實(shí)翻車記錄4.1 現(xiàn)象編譯通過但運(yùn)行時(shí)報(bào)缺少 opencv_world 庫原因CMake 鏈接的是模塊化庫但運(yùn)行時(shí)動(dòng)態(tài)鏈接器找不到.so路徑。解決在~/.bashrc里把 OpenCV 的lib目錄加進(jìn)LD_LIBRARY_PATH或者編譯時(shí)加-Wl,-rpath把路徑寫死進(jìn)可執(zhí)行文件。我一般用后者省得換終端就失效。4.2 現(xiàn)象同一張圖兩次運(yùn)行結(jié)果不一樣原因HoughCircles內(nèi)部有并行累加多線程下浮點(diǎn)累加順序不同會(huì)導(dǎo)致結(jié)果微抖。解決在調(diào)用前設(shè)cv::setNumThreads(1)犧牲一點(diǎn)速度換確定性。批量處理時(shí)這個(gè)設(shè)置尤其重要否則日志對(duì)不上。4.3 現(xiàn)象掩膜邊緣有一圈鋸齒原因cv::circle畫圓時(shí)默認(rèn)用 8 鄰域連接半徑大時(shí)邊緣呈階梯狀。解決把cv::circle的最后一個(gè)參數(shù)lineType設(shè)為cv::LINE_AA開啟抗鋸齒。代價(jià)是掩膜邊緣會(huì)有灰度過渡值后續(xù)如果要做嚴(yán)格二值判斷需要再閾值化一次。4.4 現(xiàn)象瞳孔被挖掉后虹膜環(huán)帶斷裂原因瞳孔圓和虹膜圓圓心不重合時(shí)先填白再挖黑會(huì)在偏移方向留下一塊空白。解決不要用兩次cv::circle改成逐像素判斷「到瞳孔圓心距離大于瞳孔半徑且到虹膜圓心距離小于虹膜半徑」這樣即使圓心偏移也能得到連續(xù)環(huán)帶。4.5 現(xiàn)象批量處理到某張圖直接崩潰原因circles為空時(shí)訪問circles[0]越界。解決在取圓之前加if (circles.size() 2) return emptyMask;返回全黑掩膜并記錄文件名。血淚經(jīng)驗(yàn)是批量腳本一定要把失敗文件名寫進(jìn)日志否則幾千張圖里崩一張回頭找都找不到。5. 進(jìn)階技巧用掩膜面積和圓度做自動(dòng)質(zhì)檢跑通分割只是第一步真正上線要解決的是「怎么知道這張掩膜是對(duì)的」。我的習(xí)慣是在生成掩膜后加一道自動(dòng)質(zhì)檢用兩個(gè)指標(biāo)篩掉明顯錯(cuò)誤的結(jié)果掩膜面積占比和環(huán)帶圓度。掩膜面積占比指虹膜環(huán)帶像素?cái)?shù)除以整圖像素?cái)?shù)。正視人眼里這個(gè)值通常在 0.02 到 0.15 之間太小說明圓檢小了太大說明把眼瞼也框進(jìn)去了。圓度用環(huán)帶外接矩形的寬高比衡量正常虹膜接近正圓寬高比應(yīng)在 0.8 到 1.25 之間超出說明檢出的圓被拉扁了。// 返回 true 表示掩膜通過質(zhì)檢 bool qualityCheck(const cv::Mat mask) { int area cv::countNonZero(mask); double ratio static_castdouble(area) / (mask.rows * mask.cols); if (ratio 0.02 || ratio 0.15) return false; std::vectorstd::vectorcv::Point contours; cv::findContours(mask, contours, cv::RETR_EXTERNAL, cv::CHAIN_APPROX_SIMPLE); if (contours.empty()) return false; cv::Rect box cv::boundingRect(contours[0]); double aspect static_castdouble(box.width) / box.height; return aspect 0.8 aspect 1.25; }邏輯說明countNonZero統(tǒng)計(jì)掩膜里白色像素?cái)?shù)除以總像素得到占比。findContours取最外層輪廓boundingRect拿到外接矩形寬高比就是圓度的粗略度量。參數(shù) 0.02、0.15、0.8、1.25 都是經(jīng)驗(yàn)閾值換數(shù)據(jù)集要重新統(tǒng)計(jì)。這套質(zhì)檢不能保證 100% 準(zhǔn)確但能把明顯翻車的圖擋在后續(xù)流程之外省下大量人工復(fù)核時(shí)間。從那以后我每次做虹膜分割都強(qiáng)制先跑一遍質(zhì)檢再進(jìn)特征提取寧可多花幾毫秒也不想讓一張爛掩膜污染整個(gè)批次的統(tǒng)計(jì)結(jié)果。希望幫到你。本文還有配套的精品資源點(diǎn)擊獲取