Opencv中的内存管理

时间:2014-08-26 00:12:02

标签: c opencv

我试图找到字节数 cvCreateImage(size,d,nChan); 函数调用分配给返回的指针。 (说大小:宽度= 1200,Hieght = 600,d = 32,nChan = 3)

opencv是分配对齐的内存位置还是只是随机内存位置?

使用 cvReleaseImage 功能释放的内存是否可立即供正在运行的应用程序使用? (假设一个程序,其中图像是在循环中创建和释放的,循环大约30000次,是否会导致内存不足错误或任何碎片错误)

1 个答案:

答案 0 :(得分:4)

您可以在the OpenCV source中查找问题的答案。

cvCreateImage()分配标头sizeof(IplImage)(当前为112字节)。数据有点复杂;如果您想在alloc.cpp中为image->imageSize获得精确答案grep,但它大约为size.height * size.width * d * nChan / 8(例如,8,640,000)。

内存最终由alloc.cpp中的fastMalloc()分配并与16字节边界对齐 - 但返回的指针是指向该边界上方的指针的大小。

cvReleaseImage递减引用计数,仅在引用计数达到零时释放数据。 如果您的程序只是在循环中创建和发布图像,那么除了性能之外,您应该没有问题。在下一次迭代中再次需要时,解除图像内存的设计很糟糕。

确切的分配/解除分配行为将取决于您的OpenCV是否使用CV_USE_SYSTEM_MALLOC构建。

如果使用cv::Mat而不是IplImage,OpenCV会处理所有内存管理问题automaticly。我强烈建议您重写代码以使用C ++ API(即使用Mat的)

例如,以下是您在上述评论中链接到的代码的翻译:

// Based on http://mehdi.rabah.free.fr/SSIM/SSIM.cpp
// Converted to OpenCV C++ API by B...

/*
 * The equivalent of Zhou Wang's SSIM matlab code using OpenCV.
 * from http://www.cns.nyu.edu/~zwang/files/research/ssim/index.html
 * The measure is described in :
 * "Image quality assessment: From error measurement to structural similarity"
 * C++ code by Rabah Mehdi. http://mehdi.rabah.free.fr/SSIM
 *
 * This implementation is under the public domain.
 * @see http://creativecommons.org/licenses/publicdomain/
 * The original work may be under copyrights. 
 */

//#include <cv.h>
#include <opencv2/core/core.hpp>
#include <opencv2/imgproc/imgproc.hpp>

//#include <highgui.h>
#include "opencv2/highgui/highgui.hpp"

//#include <iostream.h>
#include <iostream>
using namespace std;


/*
 * Parameters : complete path to the two image to be compared
 * The file format must be supported by your OpenCV build
 */
int main(int argc, char** argv)
{
    if(argc!=3)
        return -1;

    // default settings
    double C1 = 6.5025, C2 = 58.5225;

    /*
    IplImage
        *img1=NULL, *img2=NULL, *img1_img2=NULL,
        *img1_temp=NULL, *img2_temp=NULL,
        *img1_sq=NULL, *img2_sq=NULL,
        *mu1=NULL, *mu2=NULL,
        *mu1_sq=NULL, *mu2_sq=NULL, *mu1_mu2=NULL,
        *sigma1_sq=NULL, *sigma2_sq=NULL, *sigma12=NULL,
        *ssim_map=NULL, *temp1=NULL, *temp2=NULL, *temp3=NULL;
    */

    /***************************** INITS **********************************/
    //img1_temp = cvLoadImage(argv[1]);
    cv::Mat img1 = cv::imread(argv[1]);
    //img2_temp = cvLoadImage(argv[2]);
    cv::Mat img2 = cv::imread(argv[2]);


    //if(img1_temp==NULL || img2_temp==NULL)
    if(img1.empty() || img2.empty())
        return -1;

    //int x=img1_temp->width, y=img1_temp->height;
    //int nChan=img1_temp->nChannels, d=IPL_DEPTH_32F;
    //CvSize size = cvSize(x, y);

    //img1 = cvCreateImage( size, d, nChan);
    //img2 = cvCreateImage( size, d, nChan);

    //cvConvert(img1_temp, img1);
    img1.convertTo(img1, CV_32F);
    //cvConvert(img2_temp, img2);
    img2.convertTo(img2, CV_32F);
    //cvReleaseImage(&img1_temp);
    //cvReleaseImage(&img2_temp);


    //img1_sq = cvCreateImage( size, d, nChan);
    cv::Mat img1_sq;
    //img2_sq = cvCreateImage( size, d, nChan);
    cv::Mat img2_sq;
    //img1_img2 = cvCreateImage( size, d, nChan);
    cv::Mat img1_img2;

    //cvPow( img1, img1_sq, 2 );
    cv::pow(img1, 2, img1_sq);
    //cvPow( img2, img2_sq, 2 );
    cv::pow(img1, 2, img1_sq);
    //cvMul( img1, img2, img1_img2, 1 );
    cv::multiply(img1, img2, img1_img2);

    //mu1 = cvCreateImage( size, d, nChan);
    cv::Mat mu1;
    //mu2 = cvCreateImage( size, d, nChan);
    cv::Mat mu2;

    //mu1_sq = cvCreateImage( size, d, nChan);
    cv::Mat mu1_sq;
    //mu2_sq = cvCreateImage( size, d, nChan);
    cv::Mat mu2_sq;
    //mu1_mu2 = cvCreateImage( size, d, nChan);
    cv::Mat mu1_mu2;


    //sigma1_sq = cvCreateImage( size, d, nChan);
    cv::Mat sigma1_sq;
    //sigma2_sq = cvCreateImage( size, d, nChan);
    cv::Mat sigma2_sq;
    //sigma12 = cvCreateImage( size, d, nChan);
    cv::Mat sigma12;

    //temp1 = cvCreateImage( size, d, nChan);
    //temp2 = cvCreateImage( size, d, nChan);
    //temp3 = cvCreateImage( size, d, nChan);

    //ssim_map = cvCreateImage( size, d, nChan);
    cv::Mat ssim_map;
    /*************************** END INITS **********************************/


    //////////////////////////////////////////////////////////////////////////
    // PRELIMINARY COMPUTING
    //cvSmooth( img1, mu1, CV_GAUSSIAN, 11, 11, 1.5 );
    cv::GaussianBlur(img1, mu1, cv::Size(11, 11), 1.5);
    //cvSmooth( img2, mu2, CV_GAUSSIAN, 11, 11, 1.5 );
    cv::GaussianBlur(img2, mu2, cv::Size(11, 11), 1.5);

    //cvPow( mu1, mu1_sq, 2 );
    cv::pow(mu1, 2, mu1_sq);
    //cvPow( mu2, mu2_sq, 2 );
    cv::pow(mu2, 2, mu2_sq);
    //cvMul( mu1, mu2, mu1_mu2, 1 );
    cv::multiply(mu1, mu2, mu1_mu2);

    //cvSmooth( img1_sq, sigma1_sq, CV_GAUSSIAN, 11, 11, 1.5 );
    cv::GaussianBlur(img1_sq, sigma1_sq, cv::Size(11, 11), 1.5);
    //cvAddWeighted( sigma1_sq, 1, mu1_sq, -1, 0, sigma1_sq );
    cv::addWeighted(sigma1_sq, 1.0, mu1_sq, -1.0, 0.0, sigma1_sq);

    //cvSmooth( img2_sq, sigma2_sq, CV_GAUSSIAN, 11, 11, 1.5 );
    cv::GaussianBlur(img2_sq, sigma2_sq, cv::Size(11, 11), 1.5);
    //cvAddWeighted( sigma2_sq, 1, mu2_sq, -1, 0, sigma2_sq );
    cv::addWeighted(sigma2_sq, 1.0, mu2_sq, -1.0, 0.0, sigma2_sq);

    //cvSmooth( img1_img2, sigma12, CV_GAUSSIAN, 11, 11, 1.5 );
    cv::GaussianBlur(img1_img2, sigma12, cv::Size(11, 11), 1.5);
    //cvAddWeighted( sigma12, 1, mu1_mu2, -1, 0, sigma12 );
    cv::addWeighted(sigma12, 1.0, mu1_mu2, -1.0, 0.0, sigma12);


    //////////////////////////////////////////////////////////////////////////
    // FORMULA

    // (2*mu1_mu2 + C1)
    //cvScale( mu1_mu2, temp1, 2 );
    //cvAddS( temp1, cvScalarAll(C1), temp1 );
    cv::Mat temp3 = 2 * mu1_mu2 + C1;

    // (2*sigma12 + C2)
    //cvScale( sigma12, temp2, 2 );
    //cvAddS( temp2, cvScalarAll(C2), temp2 );

    // ((2*mu1_mu2 + C1).*(2*sigma12 + C2))
    //cvMul( temp1, temp2, temp3, 1 );
    temp3 = temp3.mul(2 * sigma12 + C2);

    // (mu1_sq + mu2_sq + C1)
    //cvAdd( mu1_sq, mu2_sq, temp1 );
    //cvAddS( temp1, cvScalarAll(C1), temp1 );
    cv::Mat temp1 = mu1_sq + mu2_sq + C1;

    // (sigma1_sq + sigma2_sq + C2)
    //cvAdd( sigma1_sq, sigma2_sq, temp2 );
    //cvAddS( temp2, cvScalarAll(C2), temp2 );

    // ((mu1_sq + mu2_sq + C1).*(sigma1_sq + sigma2_sq + C2))
    //cvMul( temp1, temp2, temp1, 1 );
    temp1 = temp1.mul(sigma1_sq + sigma2_sq + C2);

    // ((2*mu1_mu2 + C1).*(2*sigma12 + C2))./((mu1_sq + mu2_sq + C1).*(sigma1_sq + sigma2_sq + C2))
    //cvDiv( temp3, temp1, ssim_map, 1 );
    cv::divide(temp3, temp1, ssim_map);


    //CvScalar index_scalar = cvAvg( ssim_map );
    cv::Scalar index_scalar = cv::mean(ssim_map);

    // through observation, there is approximately 
    // 1% error max with the original matlab program

    cout << "(R, G & B SSIM index)" << endl ;
    cout << index_scalar.val[2] * 100 << "%" << endl ;
    cout << index_scalar.val[1] * 100 << "%" << endl ;
    cout << index_scalar.val[0] * 100 << "%" << endl ;

    /*
    // if you use this code within a program
    // don't forget to release the IplImages
    */
    return 0;
}