7. 稀疏表示之OMP,SOMP算法及openCV实现
一、前言
稀疏表示是自上世纪90年代开始,从人眼的视觉感受野获得启示,逐渐被人们所研究。现在已经发展为一种重要的信息表示方法。所谓稀疏表示是指,一个信号在过完备字典中,可以由少数个原子线性表达,
其数学模型可以表达如下:
这个数学模型解算是一个NP-hard问题,也就是说只能通过穷举去获得最优解,其时间复杂度很大,几乎无法获得其精确的解算。在这种情况下,我们常用贪婪算法去获得该模型的次最优解。本文介绍一种主流的贪婪算法——
正交匹配追踪(OMP)。
二、OMP算法
贪婪算法的核心是每次从字典的原子中选择一个最优原子来表示原始的信号。贪婪算法最大的缺点是,在贪婪算法的思想里,认为全局最优是每个局部最优得到的,这实际上很容易进入局部最优解,无法得到数据的全局最优解。
OMP作为贪婪算法中比较具有代表性的算法,其主要思想在于以下两点:
1 认为字典原子在信号投影中的越大,对信号的描述越好;
2 每一次选择的原子都与之前的原子正交。
介于以上两点,OMP算法的描述如下:
上图是从网上摘抄下来的。大概就是那样。但是值得注意的是:这样的OMP算法在解算的时候其效果往往没有ORMP算法好,目前好多人说的OMP算法其实质往往是ORMP算法。
比如:开源的工具箱SPAMS上的OMP算法解算部分,其核心就是ORMP算法。ORMP算法相比如OMP算法的不同在于,在计算完残差后对字典原子进行了另一个的拉伸(具体拉伸见后面代码部分),如下图:
三、SOMP算法
SOMP算法又叫同步OMP算法,主要思想为:相似的原子具有相同的稀疏特性。因此在对相似原子进行稀疏表示时,假设稀疏原子位于相同的位置,及其在过完备字典的选择的原子相同,OMP算法是SOMP算法在原始信号为一个原子
时的特殊情况。OMP算法可以统一到SOMP算法当中,其解算流程几乎同OMP算法部分。
四、代码实现
代码如下:
- cv::Mat ormpSparseRepresentation::ompSparseL2(const cv::Mat Dict, const cv::Mat Y, const int K)
- {
- int r = Dict.rows;
- int c = Dict.cols;
- int n = Y.cols;
- cv::Mat ERR(r,,CV_32F);
- ERR = Y;
- int size[] = {c,n};
- cv::Mat A = cv::Mat(,size,CV_32F,cv::Scalar(.f));
- QVector<int> index;
- cv::Mat U = cv::Mat::ones(,c,CV_32F);
- cv::Mat tmpA;
- for(int i = ;i<K;i++)
- {
- cv::Mat S = ERR.t()*Dict;
- cv::pow(S,,S);
- if(S.rows != )
- cv::reduce(S,S,,CV_REDUCE_SUM);
- cv::sqrt(S,S);
- S = S/U;
- if(i!=)
- {
- for(int j = ;j<index.size();j++)
- {
- S.at<float>(,index[j]) = .f;
- }
- }
- cv::Point maxLoc;
- cv::minMaxLoc(S,NULL,NULL,NULL,&maxLoc);
- int pos = maxLoc.x;
- index.append(pos);
- cv::Mat subDict;
- getColDictFormIndex(Dict,index,subDict);
- cv::Mat invSubDict;
- cv::invert(subDict,invSubDict,cv::DECOMP_SVD);
- tmpA = invSubDict*Y;
- ERR = Y - subDict*tmpA;
- cv::Mat Dict_T_Dict;
- cv::mulTransposed(subDict,Dict_T_Dict,);
- cv::Mat invDict_T_Dict;
- cv::invert(Dict_T_Dict,invDict_T_Dict,cv::DECOMP_SVD);
- cv::Mat P = (cv::Mat::eye(r,r,CV_32F) - subDict*invDict_T_Dict*subDict.t())*Dict;
- cv::pow(P,,P);
- cv::reduce(P,P,,CV_REDUCE_SUM);
- cv::sqrt(P,U);
- }
- for(int i = ;i<K;i++)
- {
- int tmpC = index[i];
- const float *inP=tmpA.ptr<float>(i);
- float *outP=A.ptr<float>(tmpC);
- for(int j = ;j<n;j++)
- {
- outP[j] = inP[j];
- }
- }
- return A;
- }
- void ormpSparseRepresentation::getColDictFormIndex(const cv::Mat Dict, const QVector<int> index, cv::Mat &res)
- {
- if(index.size() == )
- return;
- if(!Dict.data)
- return;
- int r = Dict.rows;
- int c = index.size();
- cv::Mat Dict_T;
- cv::transpose(Dict,Dict_T);
- cv::Mat res_T = cv::Mat(c,r,Dict.type());
- for(int i = ;i<index.size();i++)
- {
- int tmpC = index[i];
- const float *inP=Dict_T.ptr<float>(tmpC);
- float *outP=res_T.ptr<float>(i);
- for(int j = ;j<r;j++)
- {
- outP[j] = inP[j];
- }
- }
- cv::transpose(res_T,res);
- res_T.release();
- Dict_T.release();
- }
7. 稀疏表示之OMP,SOMP算法及openCV实现的更多相关文章
- 稀疏自动编码之反向传播算法(BP)
假设给定m个训练样本的训练集,用梯度下降法训练一个神经网络,对于单个训练样本(x,y),定义该样本的损失函数: 那么整个训练集的损失函数定义如下: 第一项是所有样本的方差的均值.第二项是一个归一化项( ...
- (转载)利用SIFT和RANSAC算法(openCV框架)实现物体的检测与定位,并求出变换矩阵(findFundamentalMat和findHomography的比较) 置顶
原文链接:https://blog.csdn.net/qq_25352981/article/details/46914837#commentsedit 本文目标是通过使用SIFT和RANSAC算法, ...
- SLAM: 图像角点检测的Fast算法(OpenCV文档)
官方链接:http://docs.opencv.org/trunk/doc/py_tutorials/py_feature2d/py_fast/py_fast.html#fast-algorithm- ...
- SAD算法在opencv上的实现代码(c++)
#include <opencv2/opencv.hpp>#include <opencv2/core/core.hpp>#include <opencv2/highgu ...
- 稀疏分解中的MP与OMP算法
MP:matching pursuit匹配追踪 OMP:正交匹配追踪 主要介绍MP与OMP算法的思想与流程,解释为什么需要引入正交? !!今天发现一个重大问题,是在读了博主的正交匹配追踪(OMP)在稀 ...
- OMP算法代码学习
正交匹配追踪(OMP)算法的MATLAB函数代码并给出单次测试例程代码 测量数M与重构成功概率关系曲线绘制例程代码 信号稀疏度K与重构成功概率关系曲线绘制例程代码 参考来源:http://blog ...
- 浅谈压缩感知(二十一):压缩感知重构算法之正交匹配追踪(OMP)
主要内容: OMP的算法流程 OMP的MATLAB实现 一维信号的实验与结果 测量数M与重构成功概率关系的实验与结果 稀疏度K与重构成功概率关系的实验与结果 一.OMP的算法流程 二.OMP的MATL ...
- MP算法、OMP算法及其在人脸识别的应用
主要内容: 1.MP算法 2.OMP算法 3.OMP算法的matlab实现 4.OMP在压缩感知和人脸识别的应用 一.MP(Matching Pursuits)与OMP(Orthogonal Matc ...
- scikit-learn 线性回归算法库小结
scikit-learn对于线性回归提供了比较多的类库,这些类库都可以用来做线性回归分析,本文就对这些类库的使用做一个总结,重点讲述这些线性回归算法库的不同和各自的使用场景. 线性回归的目的是要得到输 ...
随机推荐
- C# 创建验证码图片
using System; using System.Drawing; using System.Drawing.Drawing2D; using System.Drawing.Imaging; us ...
- js获取get值
//获取get值 function getPar(par) { //获取当前URL var local_url = document.location.href; //获取要取得的get参数位置 va ...
- 我永远的 dell 15r
陪伴我三年多的15r 让我疯狂过.努力过.更记录了我很多的成长,与很多个瞬间.看到它是有感情的.还记得第一次失去她.我好像失去了自己. 我是一个程序员.每个程序员都有自己的环境配置.不同的软件.就是 ...
- Objective-C 字典、可变字典
字典相当于c++ stl中的map 字典NSDictionary #import <UIKit/UIKit.h> #import "AppDelegate.h" int ...
- 3、HelloKhala示例说明
最简单的dome程序只需3行代码 int main() { //设置端口号 InetAddress listenAddr(USER_PORT); //将端口号绑定到Server NodeServer ...
- 初学QML之qmlRegisterType
qmlRegisterType 是一个可以将C++实现的类在QML中调用的,连接C++和QML的一个工具 首先来看QtAssistant的介绍 int qmlRegisterType(const ch ...
- www
dddd int vec_rotate(char *vec,int rotdist, int length) { int i,j,k,times; char t; times = gcd(rotdis ...
- Windows下Mysql解压缩版配置安装与卸载
安装: ①解压Mysql到合适的位置! ②以管理员身份运行命令提示符(cmd),cd C:\Documents and Settings\Administrator\桌面\mysql-5.6.24-w ...
- 混入模式(max-in)实现继承
混入模式并不是一种复制完整的对象,而是从多个对象中复制出任意的成员并将这些成员组合成一个新的对象. 实现如下: function mix(){ var arg,prop,child = {}; for ...
- [译]终极塔防——运用HTML5从头创建一个塔防游戏
翻译共享一篇CodeProject的高星力作,原文地址:http://www.codeproject.com/Articles/737238/Ultimate-Tower-Defense 下载演示项目 ...