c++画出分割图像,水平线和垂直线

1、pca 找到图像某个区域的垂直线,并画出来

复制代码
    // 1、 斑块的框+ 血管二值化图,pca 找到垂直血管壁的直线,  还是根据斑块找主轴方向吧
    // Step 1: 提取斑块左右范围内的血管像素点坐标,
    std::vector<cv::Point> points;
    for (int y = 0; y < binaryVessel.rows; y++) {
        for (int x = x1; x <= x2; x++) {
            if (binaryVessel.at<uchar>(y, x) == 255) {
                points.push_back(cv::Point(x, y));
            }
        }
    }

    if (points.empty()) {
        std::cout << "====================>error No valid pixels in the specified range.  斑块附近的血管分割区域不存在,结束函数功能 " << std::endl;
        return;

    }

    // 转换为 cv::Mat 格式
    cv::Mat pointsMat(points.size(), 2, CV_32F);
    for (int i = 0; i < points.size(); i++) {
        pointsMat.at<float>(i, 0) = points[i].x;
        pointsMat.at<float>(i, 1) = points[i].y;
    }

    // Step 2: 使用 PCA  获取垂直线段斜率,感觉有点异常,不要了,改成用最小外接矩形
    // binaryPlaque    pointsMat
    cv::PCA pca(pointsMat, cv::noArray(), cv::PCA::DATA_AS_ROW);

    // Step 3: 获取主成分(直线的方向)和主成分对应的直线上的点
    cv::Mat directions = pca.eigenvectors;
    cv::Mat linePoints = pca.mean;

    // Step 4: 计算垂直线的斜率和截距
    float kVertical = -directions.at<float>(0,0) / directions.at<float>(0,1);
    float interceptHorizontal = linePoints.at<float>(0, 1) - kVertical * linePoints.at<float>(0, 0);


    // Step 5: 绘制垂直线
    float x;
    cv::Mat image22;  // 存储处理后的图像
    cv::cvtColor(binaryVessel, image22, cv::COLOR_GRAY2BGR);  // binaryVessel 转换为彩色图像  image22,存储
    for (int i = 0; i < image22.cols; i++) {
        x = static_cast<float>(i);
        float yVertical = kVertical * x + interceptHorizontal;
        cv::circle(image22, cv::Point(static_cast<int>(x), static_cast<int>(yVertical)), 1, cv::Scalar(0, 255, 0), -1);
    }
    cv::imwrite("/home/hebin/Desktop/shiyuanyin/image2.png", image22);  // 保存画图的直线

2、使用最小外接矩形的框,来找到水平线和垂直线

相关推荐
nbsaas-boot14 分钟前
Go语言生态成熟度分析:为何Go还无法像Java那样实现注解式框架?
java·开发语言·golang
hi0_618 分钟前
03 数组 VS 链表
java·数据结构·c++·笔记·算法·链表
xiaocainiao88124 分钟前
Python 实战:构建可扩展的命令行插件引擎
开发语言·python
碧海蓝天20221 小时前
C++法则21:避免将#include放在命名空间内部。
开发语言·c++
兮动人1 小时前
Java应用全链路故障排查实战指南:从系统资源到JVM深度诊断
java·开发语言·jvm
R-sz1 小时前
导出word并且插入图片
开发语言·c#·word
CodeWithMe1 小时前
【读书笔记】《C++ Software Design》第一章《The Art of Software Design》
开发语言·c++
kyle~1 小时前
Opencv---深度学习开发
人工智能·深度学习·opencv·计算机视觉·机器人
脑袋大大的1 小时前
判断当前是否为钉钉环境
开发语言·前端·javascript·钉钉·企业应用开发
Wy. Lsy2 小时前
Kotlin基础学习记录
开发语言·学习·kotlin