窗口法估计瞳孔中心位置

1 题目

窗口法估计瞳孔中心位置。

 

2 算法原理

 

资料来源：田启川, 潘泉, 梁彦等. 瞳孔中心指导的虹膜边界快速定位算法. 系统仿真学报, 第18卷第17期, 2006年7月.

 

3 实现代码

 /* *Get postion of pupil center with window * *by LiChengyu, 2010-1-5 */ void main() { byte [,]f = LoadImg(); if (f==null) return; ShowImg("f",f); int w = f.GetLength(0); int h = f.GetLength(1); float sum = 0.0f; float temp = 100000.0f;//set "temp" as a large number, to make "temp" //get the value of "sum" in the first loop int pos_x = 0;//store the x, y axis position of pupil center int pos_y = 0; int size = 30;//the size of window, equal to the diameter of pupil //Get position of pupil center for(int x = 0; x < w-size; x += 10) for(int y = 0; y < h-size; y += 10) { sum = 0; //calculate sum of pixels in the window for(int i = x; i < x+size; i++) for(int j = y; j < y+size; j++) { sum += f[i,j]; } if(temp > sum) { temp = sum; pos_x = x; pos_y = y; } }//end of for pos_x += size / 2; pos_y += size / 2; f[pos_x,pos_y] = 255; f[pos_x+1,pos_y] = 255; f[pos_x-1,pos_y] = 255; f[pos_x,pos_y+1] = 255; f[pos_x,pos_y-1] = 255; ShowImg("f with center",f); } 

 

注：该代码可以在窗口求和部分优化，以步长为10进行窗口移动，每两个相邻窗口有20列或行数据相同，不必重复计算。

 

代码执行环境：MyImgPro 西安交大图像所 梁毅军开发

 

4 结果

图1 瞳孔原图

图2 瞳孔中心

5 评价

优点：不需二值化处理，简单易行。

缺点：窗口宽度要与瞳孔直径一直，否则无法获取瞳孔中心的准确位置；窗口大小不能过大，否则计算量增大很多。