Question

（NB：对不起，我的英语不好...）

大家好，我需要代码（C ++）的帮助：我需要使用OpenCV从灰度图像中标记区域。

为此，我将一个新的矩阵初始化为零。

我搜索第一个0（未标记的像素），获取坐标并开始递归：我看看它的8个邻居。对于每个像素：如果像素存在（如果我当前使用的像素不在边缘上），那么如果当前像素与其相邻像素之间的“距离”小于阈值，则我将其放置在同一区域中（相同的标签）并将其保存到列表中。

只要列表不为空，我就会从同一区域中选择另一个像素，然后重新开始（列表中的像素）。

当一个区域不再有像素可以使用，并且如果没有标记所有像素时，我将开始一个新的区域，依此类推。

所以问题是，在几个像素之后，我遇到了分段错误，根据gdb，这是一个来自std :: list.back（）的问题，带有“无法访问内存地址...”。

我开始使用一对std：list。我现在仅使用带有2 push_front的std :: list来在列表中添加坐标，然后使用2 back（紧随pop_back之后）来将其恢复。

我尝试使用向量没有任何好的结果。

当我使用少量图片（64x64）时，它可以正常工作。但是当涉及到更大的图片/非正方形图片时，我总是会遇到分割错误（使用gdb，似乎在“ int b = listePix.back（）”处）。

谢谢。

#include <opencv2/core/core.hpp>
#include <opencv2/highgui/highgui.hpp>
#include <opencv2/contrib/contrib.hpp>
#include <iostream>
#include <stdio.h>
#include <utility>
#include <list>

using namespace cv;
using namespace std;

void agregation(Mat src, Mat dst);
void expendReg(Mat src, Mat dst, int compteur);
pair<int,int> premiereOcc(Mat dst);

std::list <int> listePix;

int seuil = 20;

int main( int argc, char** argv )
{
    Mat image;
    image = imread(argv[1], CV_LOAD_IMAGE_GRAYSCALE);

    Mat reg = Mat::zeros(image.rows, image.cols, 0);

    agregation(image,reg);

    cout << "WORKED" << endl;

    waitKey(0);
    return 0;
}

void agregation(Mat src, Mat dst){
    bool flag = 0;
    int compteur = 0;

    do{
      compteur++;

      pair<int,int> coord = premiereOcc(dst);

      if(coord.first == -1 || coord.second == -1){
        flag = 1;
      }
      else{
        int pm = coord.first;
        int pn = coord.second;

        dst.at<uchar>(pm,pn) = compteur;
        listePix.push_front(pn);
        listePix.push_front(pm);

        expendReg(src,dst,compteur);
        listePix.clear();
      }
    }while(!flag);
}

void expendReg(Mat src, Mat dst, int compteur){
    if(listePix.size()>=2){

      int b = listePix.back();
      listePix.pop_back();
      int a = listePix.back();
      listePix.pop_back();

      int dist;

      if((a+1)<dst.rows){
          dist = abs(src.at<uchar>(a,b) - src.at<uchar>(a+1,b));
          int temp = dst.at<uchar>(a+1,b);
          if((dist<=seuil) && (temp == 0)){
            dst.at<uchar>(a+1,b) = compteur;
            listePix.push_front(b);
            listePix.push_front(a+1);
          }
      }
      if((b+1)<dst.cols){
          dist = abs(src.at<uchar>(a,b) - src.at<uchar>(a,b+1));
          int temp = dst.at<uchar>(a,b+1);
          if((dist<=seuil) && (temp == 0)){
            dst.at<uchar>(a,b+1) = compteur;
            listePix.push_front(b+1);
            listePix.push_front(a);
          }
      }
      expendReg(src,dst,compteur);
    }
}

pair<int,int> premiereOcc(Mat dst){
  //Return the coord of the first pixel without region
  for(int i = 0; i<dst.rows; i++){
    for(int j = 0; j<dst.cols; j++){
      int tmp = dst.at<uchar>(i,j);
      if(tmp == 0){
        pair<int,int> coord = make_pair(i,j);
        return coord;
      }
    }
  }
  pair<int,int> coord = make_pair(-1,-1);
  return coord;
}

EDIT2：我正在使用的图片：https://www.bogotobogo.com/Matlab/images/MATLAB_DEMO_IMAGES/blobs.png

Answer 1

我怀疑这里可能发生堆栈溢出，因为：

它仅在大图像上发生
堆栈溢出会导致分段错误：Why stack overflow causes segmentation fault instead of stack overflow in Linux?
expendReg函数是递归的

可能的修复

避免递归并使用迭代方法
限制递归深度

在非空std :: list <int>上使用std :: list.back（）时出现“分段错误”

1 个答案: