使用HoughCircles来检测和测量瞳孔和虹膜

Question

我正在尝试使用OpenCV，更具体地说是它的HoughCircles来检测和测量瞳孔和虹膜，目前我一直在玩这个函数中的一些变量，因为它要么返回0个圆圈，要么过量。下面是我正在使用的代码和测试图像。

测量光圈的代码：

 eye1 = [self increaseIn:eye1 Contrast:2 andBrightness:0];

cv::cvtColor(eye1, eye1, CV_RGBA2RGB);
cv::bilateralFilter(eye1, eye2, 75, 100, 100);

cv::vector<cv::Vec3f> circles;

cv::cvtColor(eye2, eye1, CV_RGBA2GRAY);

cv::morphologyEx(eye1, eye1, 4, cv::getStructuringElement(cv::MORPH_RECT,cv::Size(3, 3)));
cv::threshold(eye1, eye1, 0, 255, cv::THRESH_OTSU);

eye1 = [self circleCutOut:eye1 Size:50];

cv::GaussianBlur(eye1, eye1, cv::Size(7, 7), 0);

cv::HoughCircles(eye1, circles, CV_HOUGH_GRADIENT, 2, eye1.rows/4);

测量瞳孔的代码：

eye1 = [self increaseBlackPupil:eye1];
    cv::Mat eye2 = cv::Mat::zeros(eye1.rows, eye1.cols, CV_8UC3);

    eye1 = [self increaseIn:eye1 Contrast:2 andBrightness:0];

    cv::cvtColor(eye1, eye1, CV_RGBA2RGB);
    cv::bilateralFilter(eye1, eye2, 75, 100, 100);

    cv::threshold(eye2, eye1, 25, 255, CV_THRESH_BINARY);

    cv::SimpleBlobDetector::Params params;
    params.minDistBetweenBlobs = 75.0f;
    params.filterByInertia = false;
    params.filterByConvexity = false;
    params.filterByCircularity = false;
    params.filterByArea = true;
    params.minArea = 50;
    params.maxArea = 500;

    cv::Ptr<cv::FeatureDetector> blob_detector = new cv::SimpleBlobDetector(params);
    blob_detector->create("SimpleBlob");
    cv::vector<cv::KeyPoint> keypoints;
    blob_detector->detect(eye1, keypoints);

我知道图像很粗糙，我也一直试图找到一种方法来清理它并使边缘更清晰。

所以我的问题是明白：我可以做些什么来调整函数HoughCircles中的参数或更改图像以检测虹膜和瞳孔？

由于

Answer 1

好的，没有经过太多的实验，我的理解是你在使用霍夫圆检测器之前只对图像应用了双边滤镜。

在我看来，你需要在这个过程中加入一个阈值步骤。

我拍摄了您在帖子中提供的样本图片并使其经历了以下步骤：

1. 转换为灰度

2. 形态梯度

3. 阈值

4. hough circle detection。

在阈值处理步骤之后，我只得到左眼的以下图像：



这是代码：

灰阶：

   cvCvtColor(im_rgb,im_rgb,CV_RGB2GRAY);

形态：

      cv::morphologyEx(im_rgb,im_rgb,4,cv::getStructuringElement(cv::MORPH_RECT,cv::Size(size,size)));

阈值：

    cv::threshold(im_rgb, im_rgb, low, high, cv::THRESH_OTSU);

霍夫圈检测：

    cv::vector<cv::Vec3f> circles;
    cv::HoughCircles(im_rgb, circles, CV_HOUGH_GRADIENT, 2, im_rgb.rows/4); 

现在我打印时：

    NSLog(@"Found %ld cirlces", circles.size());

我明白了：

     "Found 1 cirlces"

希望这有帮助。