3

画像のスケルトンのパス/カーブで処理できるコードを開発しようとしています。2 点間のスケルトンからの点のベクトルが必要です。

ここに画像の説明を入力

このコードは、いくつかのポイントを追加した後に終了します。解決策が見つかりませんでした。

 #include "opencv2/highgui/highgui.hpp"
 #include "opencv2/imgproc/imgproc.hpp"
 #include "opencv2/imgproc/imgproc_c.h"
 #include <opencv2/legacy/compat.hpp>
 #include <opencv2/core/core.hpp>
 #include <iostream>
 #include <cmath>
 #include <list>
 using namespace cv;
 using namespace std;

  //this method to find the 8-neighbors of a point from image 
 vector<Point> search8Neighbor(cv::Mat mat,Point startPoint)
 {
  vector<Point> segment;
  Point p;
    //uchar p1 = mat.at<uchar>(startPoint.x, startPoint.y);
    uchar p2 = mat.at<uchar>(startPoint.x-1, startPoint.y);
    uchar p3 = mat.at<uchar>(startPoint.x-1, startPoint.y+1);
    uchar p4 = mat.at<uchar>(startPoint.x, startPoint.y+1);
    uchar p5 = mat.at<uchar>(startPoint.x+1, startPoint.y+1);
    uchar p6 = mat.at<uchar>(startPoint.x+1, startPoint.y);
    uchar p7 = mat.at<uchar>(startPoint.x+1, startPoint.y-1);
    uchar p8 = mat.at<uchar>(startPoint.x, startPoint.y-1);
    uchar p9 = mat.at<uchar>(startPoint.x-1, startPoint.y-1);

    //if(p1==255) segment.push_back(startPoint.x, startPoint.y);
    if (p2 == 255) {
            p.x=startPoint.x-1;
            p.y=startPoint.y;
            segment.push_back(p);
    }
    if(p3==255) {
            p.x=startPoint.x-1;
            p.y=startPoint.y+1;
            segment.push_back(p);
    }

    if(p4==255) {
            p.x=startPoint.x;
            p.y=startPoint.y+1;
            segment.push_back(p);
    }

    if(p5==255) {
            p.x=startPoint.x+1;
            p.y=startPoint.y+1;
            segment.push_back(p);
    }

    if(p6==255) {
            p.x=startPoint.x+1;
            p.y=startPoint.y;
            segment.push_back(p);
    }
    if(p7==255) {
            p.x=startPoint.x+1;
            p.y=startPoint.y-1;
            segment.push_back(p);
    }
    if(p8==255){
            p.x=startPoint.x;
            p.y=startPoint.y-1;
            segment.push_back(p);
    }

    if(p9==255) {
            p.x=startPoint.x-1;
            p.y=startPoint.y-1;
            segment.push_back(p);
    }

    return segment ;

 }
  //this method return a vector of points from a skeleton that contains all the points 
 //  between a start point "peak" and an end point 
 //this method use the idea of Breadth first search 
  vector<Point> traceLine(Mat img , Point peak)
{
    vector<Point> vect1;
    vector<Point> vect2;
    img.at<uchar>(peak.x,peak.y)=0;//
    vect1.push_back(peak);//add peak to vect1
    while(vect1.size() != 0)
        {
        Point p=vect1[0];
        vect1.erase(vect1.begin());
        vect2.push_back(p);
        vector<Point> vectN;
        vectN=search8Neighbor(img,p);
        vector<Point>::iterator it;
        it = vect1.begin();
        //cout<<vectN.size()<<endl;
        for(int i=0;i<vectN.size();i++)
            {
                 img.at<uchar>(vectN[i].x,vectN[i].y)=0;
                vect1.insert(it,vectN[i]);

            }
        }
return vect2;
}

int main( int argc, char** argv )
{
cv::Mat im = cv::imread("aa.jpg",0);
if (im.empty())
return -1;
 cv::Mat img;
cv::threshold(im, img, 155, 255, CV_THRESH_BINARY);
vector<Point> vect1;
Point p;
p.x=20;
p.y=30;
if(mat.at<uchar>(p.x-1, p.y)==255)
vect1=traceLine(img,p);
imshow("image",im);
cv::waitKey();
return 0 ;
}
4

1 に答える 1

2

Mat は、行列の通常の行/列の方法でインデックス付けされるため、必要な場合、必要mat.at<uchar>(y, x)でないmat.at<uchar>(x, y)場合、または混乱が生じます。

あなたが持っている場所：

vector<Point>::iterator it;
it = vect1.begin();
...
for(...)
    vect1.insert(it,vectN[i]);

挿入によってバッファが再割り当てされると、これは壊れます。これは、割り当てが解除された、またはおそらく別の場所に再割り当てされたメモリを指すためです。代わりに使用

 vect1.insert(vect1.begin(),vectN[i]);

(これにより、プログラムが実行しているように見える順序が少し異なります) またはpush_back()、実際に幅優先の動作を取得するために使用します。

最大の問題は search8Neighbor() にあります。uchar p2 = mat.at<uchar>(startPoint.x-1, startPoint.y);startPoint が画像の端にあるピクセルを参照するとどうなるでしょうか? 画像の外側のピクセルを参照することになり、プログラムは例外が発生するまでさまよってしまいます。ここでは、画像のどこにいるのかを確認し、画像の外側からのピクセルが近傍に含まれていないことを確認する必要があります。

編集 私が言ったように変更すると、コードが機能しないとは思いませんでした。修正を実装しましたが、動作します：

#include "opencv2/highgui/highgui.hpp"
#include "opencv2/imgproc/imgproc.hpp"
#include "opencv2/imgproc/imgproc_c.h"
#include <opencv2/legacy/compat.hpp>
#include <opencv2/core/core.hpp>
#include <iostream>
#include <cmath>
#include <list>
using namespace cv;
using namespace std;

 // Find the 8-neighbors of a point from image.
 vector<Point> search8Neighbor(cv::Mat mat,Point startPoint)
 {
    vector<Point> neighbors;
    vector<Point> offsets;
    offsets.push_back(Point(-1,  0)); // p2
    offsets.push_back(Point(-1,  1)); // p3
    offsets.push_back(Point( 0,  1)); // p4
    offsets.push_back(Point( 1,  1)); // p5
    offsets.push_back(Point( 1,  0)); // p6
    offsets.push_back(Point( 1, -1)); // p7
    offsets.push_back(Point( 0, -1)); // p8
    offsets.push_back(Point(-1, -1)); // p9

    vector<Point>::iterator it;
    for(it = offsets.begin(); it < offsets.end(); ++it)
    {
        if(!((it->x < 0 && startPoint.x == 0)
            || (it->y < 0 && startPoint.y == 0)
            || (it->x > 0 && startPoint.x == mat.cols - 1)
            || (it->y > 0 && startPoint.y == mat.rows - 1)))
        {
            Point p(startPoint + *it);
            if(mat.at<uchar>(p) == 255)
                neighbors.push_back(p);
        }
    }

    return neighbors;
 }

 //this method return a vector of points from a skeleton that contains all the points 
 //  between a start point "peak" and an end point 
 //this method use the idea of Breadth first search 
  vector<Point> traceLine(Mat img , Point peak)
{
    vector<Point> vect1;
    vector<Point> vect2;
    img.at<uchar>(peak.y, peak.x)=0;//
    vect1.push_back(peak);//add peak to vect1
    while(vect1.size() != 0)
        {
        Point p=vect1[0];
        vect1.erase(vect1.begin());
        vect2.push_back(p);
        vector<Point> vectN;
        vectN = search8Neighbor(img, p);

        //cout<< " vectN.size()=" << vectN.size()<<endl;

        for(int i = 0; i < int(vectN.size()); ++i)
            {
                img.at<uchar>(vectN[i].y, vectN[i].x)=0;
                vect1.push_back(vectN[i]);
            }
        }
    return vect2;
}

int main( int argc, char** argv )
{
    cv::Mat im = cv::imread("aa.jpg",0);
    if (im.empty())
        return -1;
    cv::Mat img;
    cv::threshold(im, img, 155, 255, CV_THRESH_BINARY);
    imshow("thresholded image",im);
    vector<Point> vect1;

    Point p(3, 32);
    uchar u = img.at<uchar>(p.x, p.y);
    if(img.at<uchar>(p) == 255)
        vect1 = traceLine(img, p);

    Mat output(im.rows, im.cols, CV_8UC3, Scalar(0, 0, 0));
    vector<Point>::iterator it;

    for(it = vect1.begin(); it < vect1.end(); ++it)
    {
        Vec3b green(0, 255, 0);
        output.at<Vec3b>(*it) = green;
    }
    imshow("output", output);
    cv::waitKey();
    return 0 ;
}
于 2013-05-19T01:30:15.240 に答える