Question

这是我在C ++学习的第一个课程，我在大学里做的是它的缝纫。起初它工作正常，但后来我在this图像上尝试了它，我得到了一个分段错误。我试过调试，似乎出现在

for(int j=0; j<n; j++){
   C[0][j]=sobel(0,j, im);
}

我正在使用

 g++ -g -std=c++11 myprogram.cpp -o app `pkg-config --cflags --libs opencv`

在Ubuntu上编译。也许我的数组C有问题？我似乎也无法获得初始化程序列表。如果我试试

my_args = { left, center, right }

而不是

      min_args[0]=left;
      min_args[1]=center;
      min_args[2]=right;

我收到错误，即使我专门启用了c ++ 11。

编辑：我现在用malloc初始化我的数组（我之后也使用了免费）。

int *C; C=(int *)malloc(n*m*sizeof(int));

现在它有效！怎么会这样？在C ++中它不应该是一样的吗？这可能会使我的程序变慢吗？

这是我的代码：

＆＃13;

#include <opencv2/opencv.hpp>
#include <iostream>
#include <algorithm>

using namespace cv;
using namespace std;

struct MinArgmin{
  float min;
  int argmin;
};

struct twoMat{
  Mat gray;
  Mat color;
};


float sobel(int i, int j, const Mat &im);
int m; int n;
void find_seam(int arr[], const Mat &im);
MinArgmin getMinArgmin(const float* arr, const int length);
Mat seam_carving(const Mat &imcol, Mat &imgray, int m_neu, int n_new);
twoMat delete_seam(const int *arr, const Mat &im_gray, const Mat &im_col);


int main( int argc, char** argv )
{

  Mat original_image;
  Mat im_gray;
  original_image = imread("Camera_obscura.jpg", 1 );
  im_gray = imread("Camera_obscura.jpg", 0 );

  if( argc == 2 || !original_image.data )
    {
      printf( "No image data \n" );
      return -1;
    }

  m=im_gray.rows;
  n=im_gray.cols;
  Mat final_result=seam_carving(original_image, im_gray, m-70, n-80);
  namedWindow( "Result", WINDOW_AUTOSIZE );
  imshow( "Result", final_result);
  namedWindow( "Original", WINDOW_AUTOSIZE );
  imshow( "Original", original_image);  


  waitKey(0);
  return 0;
}

float sobel(int i, int j, const Mat &im){
  int pos_j=std::max(j-1,0);
  int pos_i=std::max(i-1,0);
  int m_i=std::min(i+1,im.rows);
  int n_j=std::min(j+1,im.cols);
  float ix=-2*im.at<uchar>(i,pos_j) + 2* im.at<uchar>(i,n_j) - im.at<uchar>(m_i,pos_j) + im.at<uchar>(m_i,n_j) - im.at<uchar>(pos_i,pos_j) + im.at<uchar>(pos_i,n_j);
  float iy=im.at<uchar>(m_i,pos_j) + 2* im.at<uchar>(m_i,j) + im.at<uchar>(m_i,n_j) - ( im.at<uchar>(pos_i,pos_j) + 2* im.at<uchar>(pos_i,j) + im.at<uchar>(pos_i,n_j));
  return abs(ix) + abs(iy);
}



void find_seam(int arr[], const Mat &im){
  int m=im.rows;
  int n=im.cols;
  float C[m][n];
  int trace[m][n];
  float min_args[3];
  int j_min;
  MinArgmin cj;
  for(int j=0; j<n; j++){
    C[0][j]=sobel(0,j, im);
  }
  for(int i=1; i<m; i++){
    for(int j=0; j<n; j++){
      float left;
      float right;
      float center;
      C[i-1][j-1];
      center=C[i-1][j];
      right=C[i-1][j+1];
      if(j==0){
        right=C[i-1][j+1];
        center=C[i-1][j];
        left=right+1; }
        else if(j==n-1) {
          center=C[i-1][j];
          left=C[i-1][j-1];
          right=left+1;
        }
        else {
          left=C[i-1][j-1];
          center=C[i-1][j];
          right=C[i-1][j+1];
        }
      min_args[0]=left;
      min_args[1]=center;
      min_args[2]=right;
      cj=getMinArgmin(min_args, 3);
      C[i][j]=cj.min+sobel(i,j, im);

      j_min=cj.argmin-1;
      trace[i][j]=j+j_min;
    }
  }

  int optcol=getMinArgmin(C[m-1],n).argmin;

  arr[m-1]=optcol;
  for(int k=m-2; k>=0; k--){
    arr[k]=trace[k+1][arr[k+1]];
  }
}

MinArgmin getMinArgmin(const float* arr, const int length){
  int argmin=0;
  float min=arr[0];
  for(int i=1; i<length; i++){
    if(arr[i]<min){min=arr[i]; argmin=i;}
  }
  return {min, argmin};
}

Mat seam_carving(const Mat &imcol, Mat &imgray, int m_new, int n_new){
  Mat result=imcol.clone();
  int horz_diff=imcol.rows-m_new;
  int vert_diff=imcol.cols-n_new;
  int m_n=max(imcol.rows, imcol.cols);
  int seam[m_n];
  for(int i=0; i<vert_diff; i++){
      find_seam(seam,imgray);
      twoMat both=delete_seam(seam,imgray,result);
      imgray=both.gray;
      result=both.color;
  }
  imgray=imgray.t();
  result=result.t();
  for(int i=0; i<horz_diff; i++){
      find_seam(seam,imgray);
      twoMat both=delete_seam(seam,imgray,result);
      imgray=both.gray;
      result=both.color;
  }
  result=result.t();
  return result;
}

twoMat delete_seam(const int *arr, const Mat &imgray, const Mat &im_col){
  if(imgray.cols!=im_col.cols || imgray.rows!=im_col.rows){
    printf("Matrices don't have the same dimension \n");
  }
  int m=imgray.rows;
  int n=imgray.cols;
  Mat temp_gray(m,n-1,CV_8UC1);
  Mat temp_col(m,n-1,CV_8UC3);
  for(int i=0; i<m; i++){
    for(int j=0; j<n-1; j++){
      if(j<arr[i]){
        temp_gray.at<uchar>(i,j)=imgray.at<uchar>(i,j);
        temp_col.at<Vec3b>(i,j)=im_col.at<Vec3b>(i,j);
      }
      else {
        temp_gray.at<uchar>(i,j)=imgray.at<uchar>(i,j+1);
        temp_col.at<Vec3b>(i,j)=im_col.at<Vec3b>(i,j+1);
      }
    }
  }
  return {temp_gray, temp_col};
}

＆＃13;

C ++ OpenCV - 程序适用于某些图像，为其他人获取分段错误

0 个答案: