OpenCV - 将向量的向量转换为Mat

Question

我想将一个 vector<vector<double> > 转换为 Mat，因为我需要对这些值应用自定义平滑滤波器。

下面的图片展示了正确的数值：

我尝试过以下方法：

std::vector<std::vector<double> > angles;
calculateAngles(angles);
Mat matAngles(angles.size(), angles.at(0).size(), CV_64FC1, angles.data());

但是第一列中的值被错误地转换为2.12566e-314。

生成的图像：

我也尝试直接将这些值放入Mat中。

void calculateAngles(cv::Mat& im, cv::Mat& angles, int blockSize, int(*f)(int x, int y), int(*g)(int x, int y)){

static int ySobel[3][3] = {{-1, 0, 1}, {-2, 0, 2}, {-1, 0, 1}};
static int xSobel[3][3] = {{-1, -2, -1}, {0, 0, 0}, {1, 2, 1}};

angles.create(cv::Size(im.cols/blockSize+1, im.rows/blockSize+1), CV_64FC1);

int nominator;
int denominator;
int GX, GY;
double angle;
for(int i = 1; i < im.cols; i += blockSize){
    for(int j = 1; j < im.rows; j += blockSize){
        nominator = 0;
        denominator = 0;
        for(int k = i; k < std::min(i + blockSize, im.cols-1); k++){
            for(int l = j; l < std::min(j + blockSize, im.rows-1); l++){
                GX = applyKernelAt(im, xSobel, k, l);
                GY = applyKernelAt(im, ySobel, k, l);

                nominator += f(GX, GY);
                denominator += g(GX, GY);
            }
        }

        angle = (PI + std::atan2(nominator, denominator)) / 2;
        angles.at<double>((i-1)/blockSize, (j-1)/blockSize) = angle;
    }
}
}

但是在一定点之后，数值会重复。 以下是数值的打印方式。

void drawLines(cv::Mat& src, cv::Mat& dst, cv::Mat& angles, int blockSize){
cv::Size size = src.size();
dst.create(src.size(), src.type());

// for a white background
dst += 255;

cv::cvtColor(dst, dst, CV_GRAY2BGR);

int x1,y1,x2,y2;

for(int i = 1; i < size.width; i += blockSize){
    for(int j = 1; j < size.height; j += blockSize){
        double tang = std::tan(angles.at<double>((i-1)/blockSize,(j-1)/blockSize));

        if( tang >= -1 && tang <= 1 ){
            x1 = i;
            y1 = (-blockSize/2) * tang + j + blockSize/2;
            x2 = i + blockSize;
            y2 = (blockSize/2) * tang + j + blockSize/2;
        }else{
            x1 = i + blockSize/2 + blockSize/(2*tang);
            y1 = j + blockSize/2;
            x2 = i + blockSize/2 - blockSize/(2*tang);
            y2 = j -blockSize/2;
        }

        cv::line(dst, cv::Point(x1,y1), cv::Point(x2,y2), cv::Scalar(0,0,255));
    }
}
}

- Rayron Victor

Answer 1

这是一种从向量的向量中填充cv::Mat的方法。

std::vector<std::vector<double> > angles;
cv::Mat matAngles(angles.size(), angles.at(0).size(), CV_64FC1);
for(int i=0; i<matAngles.rows; ++i)
     for(int j=0; j<matAngles.cols; ++j)
          matAngles.at<double>(i, j) = angles.at(i).at(j);

- bruvel

---

谢谢，我把这个放在calculateAngles(angles)之后，它就起作用了。但是，如果我在这个函数内部进行转换并将Mat&作为参数传递，问题仍然存在。 - Rayron Victor

Answer 2

如果你想将一个大的std::vector of std::vector数据临时复制到cv::Mat中而不进行复制，则最有效的方法是：

std:vector<std::vector<float> > OrigSamples;
...
... // Here OrigSamples is filled with some rows with the same row size
...

// Create a new, _empty_ cv::Mat with the row size of OrigSamples
cv::Mat NewSamples(0, OrigSamples[0].size(), cv::DataType<float>::type);

for (unsigned int i = 0; i < OrigSamples.size(); ++i)
{
  // Make a temporary cv::Mat row and add to NewSamples _without_ data copy
  cv::Mat Sample(1, OrigSamples[0].size(), cv::DataType<float>::type, OrigSamples[i].data());

  NewSamples.push_back(Sample);
}

请注意：

• 不要操纵NewSamples中的数据，只使用它。
• OrigSamples必须在作用域内（不能被销毁）或修改，直到使用NewSamples。

否则会出现内存损坏/崩溃。

- kecsap