OpenCV输出均值、标准差、协方差

1、概述

  案例:使用OpenCV输出图片的均值和标准差以及协方差矩阵

  相关定义:

  1.均值:平均值,指讲一个数据集合中的值相加再除以数据集中数据的个数得出的值

  2.方差:数据集中每个样本值与全体样本值的平均数之差的平方值的平均数

  3.标准差:标准差是方差的算数平方根(ps:定义参考方差)

  4.协方差:在概率论和统计学中用于衡量两个变量的总体误差

  函数定义:

均值及标准差:
void meanStdDev(InputArray src, OutputArray mean, OutputArray stddev,
                             InputArray mask=noArray());
1.src:输入图像
2.mean:输出参数--->均值
3.stddev:输出参数--->标准差
协方差矩阵、均值矩阵<br>void calcCovarMatrix( InputArray samples, OutputArray covar,
                                   InputOutputArray mean, int flags, int ctype = CV_64F);<br>samples:输入数据,一般channgle==1<br>covar:输出参数表示协方差矩阵<br>mean:输出参数表示,均值矩阵<br>flags:操作标志:CV_COVAR_SCRAMBLED,CV_COVAR_NORMAL,CV_COVAR_USE_AVG,CV_COVAR_SCALE,CV_COVAR_ROWS,CV_COVAR_COLS.

2、代码演示

Face_Means_Covar::Face_Means_Covar(QWidget *parent)
    : MyGraphicsView{parent}
{
    this->setWindowTitle("计算Mat均值、标准差、及协方差");
    QLabel *labelTitle = new QLabel(this);
    labelTitle->setText("输出结果如下:");
    //均值
    QLabel *labelMeansTitle = new QLabel(this);
    labelMeansTitle->setText("均值:");
    labelMeansTitle->move(0,labelTitle->y()+labelTitle->height()+10);
    labelMeansValue = new QLabel(this);
    labelMeansValue->setFixedWidth(150);
    labelMeansValue->move(labelMeansTitle->x()+labelMeansTitle->width()+10,labelMeansTitle->y());
    //标准差
    QLabel *labelStddevTitle = new QLabel(this);
    labelStddevTitle->setText("标准差:");
    labelStddevTitle->move(0,labelMeansTitle->y()+labelMeansTitle->height()+10);
    labelStddevValue = new QLabel(this);
    labelStddevValue->setFixedWidth(150);
    labelStddevValue->move(labelStddevTitle->x()+labelStddevTitle->width()+10,labelStddevTitle->y());
    //协方差
    QLabel *labelCovarTitle = new QLabel(this);
    labelCovarTitle->setText("协方差:");
    labelCovarTitle->move(0,labelStddevTitle->y()+labelStddevTitle->height()+10);
    labelCovarValue = new QLabel(this);
    labelCovarValue->setFixedWidth(150);
    labelCovarValue->move(labelCovarTitle->x()+labelCovarTitle->width()+10,labelCovarTitle->y());

}

void Face_Means_Covar::dropEvent(QDropEvent *event){
    path = event->mimeData()->urls().at(0).toLocalFile();
    qDebug()<<"文件路径:"<<path;
    showMeansAndCovar(path.toStdString().c_str());

}

void Face_Means_Covar::showMeansAndCovar(const char *filePath){
    Mat src = imread(filePath);
    if(src.empty()){
        qDebug()<<"图片文件为空";
        return;

    }
    imshow("src",src);

    //计算图像均值以及标准差
    Mat means,stddev;
    meanStdDev(src,means,stddev,noArray());
    //输出均值矩阵,图像有多少通道means就有多少个值。例如:三通道图像,则有第一通道的均值,第二通道的均值,第三通道的均值
    cout << means<<endl;
    //输出标准差矩阵,图像有多少个通道stddev就有多少个值。例如:单通道图像,则0的位置就是通道标准差
    cout<<stddev<<endl;

    //计算所有通道的平均值
    double resultMeans = 0.0;
    for(int i=0;i<means.rows;i++){
        resultMeans+=means.at<double>(i);
         qDebug()<<"均值:"<< means.at<double>(i);
    }
    QString rMeans=QString("%1").arg(resultMeans/3);
    labelMeansValue->setText(rMeans);


    //计算所有通道的标准差
    double resultStddev=0.0;
    for(int i=0;i<stddev.rows;i++){
        resultStddev+=stddev.at<double>(i);
        qDebug()<<"标准差:"<<stddev.at<double>(i);
    }
    QString rStddev=QString("%1").arg(resultStddev/3);
    labelStddevValue->setText(rStddev);


    //计算协方差
    Mat covar,means2;
    cvtColor(src,src,COLOR_BGR2GRAY);
    //这里的输入需要是单通道的channgle==1
    calcCovarMatrix(src,covar,means2,COVAR_NORMAL|cv::COVAR_ROWS);
    qDebug()<<"以下是协方差矩阵";
    //一定要搞一个小一点的图测试
    cout << covar<<endl;
    cout << means2<<endl;



}

3、示例图片

 本文福利,费领取Qt开发学习资料包、技术视频,内容包括(C++语言基础,Qt编程入门,QT信号与槽机制,QT界面开发-图像绘制,QT网络,QT数据库编程,QT项目实战,QSS,OpenCV,Quick模块,面试题等等)↓↓↓↓↓↓见下面↓↓文章底部点击费领取↓↓

  • 1
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
均值方差模型是指通过最小化投资组合的方差来实现风险最小化,同时保证预期收益率达到一定水平的投资模型。其数学公式为: $$\min_{w} \frac{1}{2}w^{T}\Sigma w$$ $$s.t. \quad w^{T}\mu = \bar{r}$$ $$\quad \quad w^{T}1 = 1$$ 其中,$w$为投资组合中每个资产的权重向量,$\Sigma$为资产收益率的协方差矩阵,$\mu$为资产收益率的期望向量,$\bar{r}$为预期收益率,$1$为全为 $1$ 的列向量。 对于均值方差模型的优化问题,可以使用Python中的SciPy库中的优化函数来解决,代码如下: ```python import numpy as np from scipy.optimize import minimize # 定义目标函数 def objective(w, cov): return 0.5 * w.dot(cov).dot(w) # 定义约束条件 def constraint1(w): return w.sum() - 1 def constraint2(w, rets, target_return): return rets.dot(w) - target_return # 定义初始权重向量 w0 = np.ones(cov.shape[0]) / cov.shape[0] # 定义约束条件 cons = [{'type': 'eq', 'fun': constraint1}, {'type': 'eq', 'fun': constraint2, 'args':(rets, target_return)}] # 最小化目标函数 opt = minimize(objective, w0, args=(cov,), method='SLSQP', constraints=cons) # 输出结果 weights = opt.x print('Optimal weights:', weights) print('Minimized volatility:', np.sqrt(opt.fun)) ``` 其中,`cov`为资产收益率的协方差矩阵,`rets`为资产收益率的期望向量,`target_return`为目标预期收益率。函数`objective()`为目标函数,函数`constraint1()`和`constraint2()`为约束条件。最小化函数使用`minimize()`,采用SLSQP算法进行求解。最终得到的`weights`即为最优投资组合的权重向量,`np.sqrt(opt.fun)`即为最小化的投资组合的方差

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值