使用OpenCV增強手掌靜脈

Question

我正試圖在OpenCV中實現一種算法，以提供手掌靜脈模式的細節。我基於我在互聯網上找到的名爲「A Contactless Biometric System Using Palm Print and Palm Vein Features」的論文。我感興趣的部分是章節3.2預處理。那裏展示了所涉及的步驟。

我想用OpenCV來實現這個實現，但是直到現在我仍然堅持下去。特別是他們對低通濾波器的響應使用拉普拉斯濾波器來隔離主脈，但無論我嘗試的參數如何，我的結果都會非常嘈雜！

任何幫助將不勝感激！

Answer 1

好吧，終於我已經明白瞭如何去做。這是我的代碼：

#include <opencv2/core/core.hpp> 
#include <opencv2/highgui/highgui.hpp> 
#include <opencv2/imgproc/imgproc.hpp> 

#define THRESHOLD 150 
#define BRIGHT 0.7 
#define DARK 0.2 

using namespace std; 
using namespace cv; 

int main() 
{ 

    // Read source image in grayscale mode 
    Mat img = imread("roi.png", CV_LOAD_IMAGE_GRAYSCALE); 

    // Apply ??? algorithm from https://stackoverflow.com/a/14874992/2501769 
    Mat enhanced, float_gray, blur, num, den; 
    img.convertTo(float_gray, CV_32F, 1.0/255.0); 
    cv::GaussianBlur(float_gray, blur, Size(0,0), 10); 
    num = float_gray - blur; 
    cv::GaussianBlur(num.mul(num), blur, Size(0,0), 20); 
    cv::pow(blur, 0.5, den); 
    enhanced = num/den; 
    cv::normalize(enhanced, enhanced, 0.0, 255.0, NORM_MINMAX, -1); 
    enhanced.convertTo(enhanced, CV_8UC1); 

    // Low-pass filter 
    Mat gaussian; 
    cv::GaussianBlur(enhanced, gaussian, Size(0,0), 3); 

    // High-pass filter on computed low-pass image 
    Mat laplace; 
    Laplacian(gaussian, laplace, CV_32F, 19); 
    double lapmin, lapmax; 
    minMaxLoc(laplace, &lapmin, &lapmax); 
    double scale = 127/ max(-lapmin, lapmax); 
    laplace.convertTo(laplace, CV_8U, scale, 128); 

    // Thresholding using empirical value of 150 to create a vein mask 
    Mat mask; 
    cv::threshold(laplace, mask, THRESHOLD, 255, CV_THRESH_BINARY); 

    // Clean-up the mask using open morphological operation 
    morphologyEx(mask,mask,cv::MORPH_OPEN, 
     getStructuringElement(cv::MORPH_ELLIPSE, cv::Size(5,5))); 

    // Connect the neighboring areas using close morphological operation 
    Mat connected; 
    morphologyEx(mask,mask,cv::MORPH_CLOSE, 
     getStructuringElement(cv::MORPH_ELLIPSE, cv::Size(11,11))); 

    // Blurry the mask for a smoother enhancement 
    cv::GaussianBlur(mask, mask, Size(15,15), 0); 

    // Blurry a little bit the image as well to remove noise 
    cv::GaussianBlur(enhanced, enhanced, Size(3,3), 0); 

    // The mask is used to amplify the veins 
    Mat result(enhanced); 
    ushort new_pixel; 
    double coeff; 
    for(int i=0;i<mask.rows;i++){ 
     for(int j=0;j<mask.cols;j++){ 
      coeff = (1.0-(mask.at<uchar>(i,j)/255.0))*BRIGHT + (1-DARK); 
      new_pixel = coeff * enhanced.at<uchar>(i,j); 
      result.at<uchar>(i,j) = (new_pixel>255) ? 255 : new_pixel; 
     } 
    } 

    // Show results 
    imshow("frame", img); 
    waitKey(); 

    imshow("frame", result); 
    waitKey(); 

    return 0; 
} 

所以本文的主要步驟如下。對於某些部分，我已經啓發了我自己找到的代碼。我申請的第一個處理就是這種情況，我找到了here。此外，對於高通濾波器（拉普拉斯算子），我已經啓發了我自己在OpenCV 2計算機視覺應用程序編程指南中給出的代碼。

最後，我通過允許修改背景的亮度和靜脈的黑暗度（參見定義BRIGHT和DARK）做了一些小的改進。我也決定模糊一下面具，以獲得更「自然」的增強效果。

---

這裏的結果（來源/紙結果/我的結果）：