首页 文章

从X射线图像中提取手骨

提问于
浏览
7

我有一只手的X射线图像 . 我需要自动提取骨骼 . 我可以使用不同的技术轻松分割手 . 但我需要得到骨头并使用这些技术无济于事 . 有些骨头比其他骨头更亮,所以如果我使用阈值处理,其中一些骨头会消失,而其他骨骼则会变得更清晰 . 我想也许我应该只限制一个手的区域?是否可以限制不是正方形的投资回报率?哦,也许你有任何其他解决方案,建议?也许有一些像OpenCV这样的图书馆?任何帮助都会非常棒!

Extended:

Raw Image

Expected Output

Raw Image Expected Output

1 回答

  • 8

    一种方法可能是从图像中 segment the hand and fingers

    enter image description here

    然后使用 hand silhouette 创建另一个图像:

    enter image description here

    一旦你有了轮廓,你可以侵蚀图像,使它变得更小 . 这用于 subtract the hand from the hand & fingers image ,导致手指:

    enter image description here

    下面的代码显示了执行此方法:

    void detect_hand_and_fingers(cv::Mat& src);
    void detect_hand_silhoutte(cv::Mat& src);
    
    int main(int argc, char* argv[])
    {
        cv::Mat img = cv::imread(argv[1]);
        if (img.empty())
        {
            std::cout << "!!! imread() failed to open target image" << std::endl;
            return -1;        
        }
    
        // Convert RGB Mat to GRAY
        cv::Mat gray;
        cv::cvtColor(img, gray, CV_BGR2GRAY);
        cv::Mat gray_silhouette = gray.clone();
    
        /* Isolate Hand + Fingers */
    
        detect_hand_and_fingers(gray);
        cv::imshow("Hand+Fingers", gray);
        cv::imwrite("hand_fingers.png", gray);
    
        /* Isolate Hand Sillhoute and subtract it from the other image (Hand+Fingers) */
    
        detect_hand_silhoutte(gray_silhouette);
        cv::imshow("Hand", gray_silhouette);
        cv::imwrite("hand_silhoutte.png", gray_silhouette);
    
        /* Subtract Hand Silhoutte from Hand+Fingers so we get only Fingers */
    
        cv::Mat fingers =  gray - gray_silhouette;
        cv::imshow("Fingers", fingers);
        cv::imwrite("fingers_only.png", fingers);
        cv::waitKey(0);
    
        return 0;
    }
    
    void detect_hand_and_fingers(cv::Mat& src)
    {        
        cv::Mat kernel = cv::getStructuringElement(cv::MORPH_ELLIPSE, cv::Size(3,3), cv::Point(1,1));
        cv::morphologyEx(src, src, cv::MORPH_ELLIPSE, kernel);    
    
        int adaptiveMethod = CV_ADAPTIVE_THRESH_GAUSSIAN_C; // CV_ADAPTIVE_THRESH_MEAN_C, CV_ADAPTIVE_THRESH_GAUSSIAN_C
        cv::adaptiveThreshold(src, src, 255, 
                              adaptiveMethod, CV_THRESH_BINARY, 
                              9, -5);
    
        int dilate_sz = 1;
        cv::Mat element = cv::getStructuringElement(cv::MORPH_ELLIPSE,
                                           cv::Size(2*dilate_sz, 2*dilate_sz),
                                           cv::Point(dilate_sz, dilate_sz) );
        cv::dilate(src, src, element);
    }
    
    void detect_hand_silhoutte(cv::Mat& src)
    {
        cv::Mat kernel = cv::getStructuringElement(cv::MORPH_ELLIPSE, cv::Size(7, 7), cv::Point(3, 3));
        cv::morphologyEx(src, src, cv::MORPH_ELLIPSE, kernel);        
    
        int adaptiveMethod = CV_ADAPTIVE_THRESH_MEAN_C; // CV_ADAPTIVE_THRESH_MEAN_C, CV_ADAPTIVE_THRESH_GAUSSIAN_C
        cv::adaptiveThreshold(src, src, 255, 
                              adaptiveMethod, CV_THRESH_BINARY, 
                              251, 5); // 251, 5
    
        int erode_sz = 5;
        cv::Mat element = cv::getStructuringElement(cv::MORPH_ELLIPSE,
                                           cv::Size(2*erode_sz + 1, 2*erode_sz+1),
                                           cv::Point(erode_sz, erode_sz) );
        cv::erode(src, src, element);
    
        int dilate_sz = 1;
        element = cv::getStructuringElement(cv::MORPH_ELLIPSE,
                                           cv::Size(2*dilate_sz + 1, 2*dilate_sz+1),
                                           cv::Point(dilate_sz, dilate_sz) );
        cv::dilate(src, src, element);
    
        cv::bitwise_not(src, src);
    }
    

相关问题