CN107369137A - 一种基于局部线性近似的不完备图像补齐方法 - Google Patents
一种基于局部线性近似的不完备图像补齐方法 Download PDFInfo
- Publication number
- CN107369137A CN107369137A CN201710495129.8A CN201710495129A CN107369137A CN 107369137 A CN107369137 A CN 107369137A CN 201710495129 A CN201710495129 A CN 201710495129A CN 107369137 A CN107369137 A CN 107369137A
- Authority
- CN
- China
- Prior art keywords
- new
- local linear
- model
- repairing model
- incomplete image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000005498 polishing Methods 0.000 title claims abstract description 13
- 239000011159 matrix material Substances 0.000 claims abstract description 20
- 238000012986 modification Methods 0.000 claims description 2
- 230000004048 modification Effects 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 claims 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
- G06T5/94—Dynamic range modification of images or parts thereof based on local image properties, e.g. for local contrast enhancement
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明提供一种基于局部线性近似对不完备图像进行补齐方法,步骤一,根据不完备图像中待恢复数据建立修复模型X;所述修复模型X为具有缺失值的矩阵;步骤二,对所述修复模型X进行初始化,对所述修复模型X缺失值进行索引,给出近邻数目k和最大迭代次数M;步骤三,判断修复模型运算次数N是否小于迭代次数M;若满足判断条件,程序进入下一步,否则令X=Xnew终止,输出等步骤;该方法将邻域近似的思想引入矩阵填充,定义了局部线性近似及重建错误率的公式,提出了局部线性近似的矩阵填充算法。
Description
技术领域
本发明属于数据计算方法,尤其涉及矩阵填充方法,具体是指一种基于局部线性近似的不完备图像补齐方法。
背景技术
在很多领域中,获取的数据是不完整的,即可能存在部分缺失。而矩阵填充是一个处理不完备数据的有效方法,它可以从较少的观测样本恢复矩阵。其广泛出现在诸多科研和工程领域,例如协同过滤、图像修复、视频去噪、遥感技术等等。
通常,矩阵填充的方法主要围绕几个假设:低秩或近似低秩的假设;缺失值是均匀分布的假设。问题的关键是最小化矩阵的秩,在实际应用中可将其转化为最小化核范数。然而,基于低秩的方法,不仅复杂,而且没有很好的适用性。
局部线性嵌入可以由局部线性拟合发现全局的结构,每个局部小段可以线性地近似,每个点可以表示成其近邻的线性加权和。受此思想的启发,本发明从邻域的局部线性重建的角度来近似缺失数据,提出了不完备数据局部线性近似的定义和公式。在此之上,给出了一种不完备数据的局部线性近似方法,实现矩阵填充。
发明内容
针对现有技术存在的技术问题,本发明提供一种基于局部线性近似对不完备图像进行补齐方法,该方法将邻域近似的思想引入矩阵填充,定义了局部线性近似及重建错误率的公式,提出了局部线性近似的矩阵填充算法。
为了解决现有技术中存在技术问题,本发明采用如下技术方案:
一种基于局部线性近似对不完备图像进行补齐方法,
步骤一,根据不完备图像中待恢复数据建立修复模型X;所述修复模型X为具有缺失值的矩阵;
步骤二,对所述修复模型X进行初始化,对所述修复模型X缺失值进行索引,给出近邻数目k和最大迭代次数M;
步骤三,判断修复模型运算次数N是否小于迭代次数M;若满足判断条件,程序进入下一步,否则令X=Xnew终止,输出;
步骤四,对修复模型X采用局部线性近似方法估算缺失值建立新模型Xnew;
步骤五,对新模型Xnew进行||X-Xnew||2<tolerance*||X||2判断,若满足判断条件,则令X=Xnew终止,输出;否则令N=N+1返回步骤三。
所述步骤一中修复模型X通过公式(1)获得初次迭代的邻域;
其中,m是属性个数,δ(xik)表示狄垃克函数,若xik不是缺失值,δ(xik)=1;若xik是缺失值,δ(xik)=0;xik表示第k个属性下x样本的值。
所述步骤四中局部线性近似方法通过公式(2)获得修改模型重建系数Wij;
其中,XMiss表示缺失值集合,xi和xj表示矩阵X的行,Wij是xj对xi的最优线性重建系数,Nk表示xi的k近邻集合。
所述步骤四中通过公式(3)对存在的缺失值xip进行估计,以估计得到的xip构成新模型Xnew;
有益效果
本发明提出的一种基于局部线性近似的不完备图像补齐方法同时在行列维度(样本空间和属性空间)借用邻域的思想来局部近似矩阵当中的缺失值,计算简单,更加灵活,为矩阵填充提供了一个新的思路和角度。和其他低秩的方法相比,该方法通过邻域使用局部线性重建来近似缺失值,并不需要过多的假设,除此之外,局部线性近似的方法还可以处理高噪声的数据。
附图说明
图1是本发明流程图。
图2是利用本发明进行图像补齐的结果图。
具体实施方式
下面结合附图对本发明作出详细说明。
如图1所述,本发明提供一种基于局部线性近似对不完备图像进行补齐方法,包括如下步骤:
步骤一101,获取不完备图像中待恢复数据建立修复模型X;所述修复模型X为具有缺失值的矩阵;
步骤二102,对所述修复模型X进行初始化,给出近邻数目k和最大迭代次数M;给出迭代数M;所述步骤一中修复模型X通过公式(1)获得初次迭代的邻域;
其中,m是属性个数,δ(xik)表示狄垃克函数,若xik不是缺失值,δ(xik)=1;若xik是缺失值,δ(xik)=0;xik表示第k个属性下x样本的值。
步骤三103,判断修复模型运算次数N是否小于迭代次数M;若满足判断条件,程序进入下一步,否则令终止;
步骤四(104,105)对修复模型X采用局部线性近似方法估算缺失值建立新模型Xnew;所述步骤四中局部线性近似方法通过公式(2)获得修复模型重建系数Wij;
其中,XMiss表示缺失值集合,xi和xj表示矩阵X的行,Wij是xj对xi的最优线性重建系数,Nk表示xi的k近邻集合。
所述步骤四中通过公式(3)生成缺失值xip估算,以此构成新模型Xnew;
本发明中利用局部线性近似Locally linear approximation,即LLA,分别在行空间(样本空间)和列空间(属性空间)进行近似。即通过线性重建系数Wij生成修复模型X缺失值xip的近似估算进而建立新模型Xnew。
步骤五(106,107,108,109)对新模型Xnew进行||X-Xnew||2<tolerance*||X||2判断,若满足判断条件,则令X=Xnew终止,输出;否则令N=N+1返回步骤三。本发明通过重复步骤三、步骤四达到最大迭代次数,X=Xnew,算法结束,输出当前还原矩阵X。
实验例1:
通过将本发明方法在实际数据集LENA上的运行,显示其有效性。数据
集MATLAB版本来源于:
http://www.cad.zju.edu.cn/home/dengcai/Data/data.html
对比算法包括:OPT、SVT、TNN、FANCL、IRNN、LMAFIT
缺失方式分别为:CHUNK、LINE、RND
缺失率为0.1、0.3、0.5、0.7、0.9
运行的结果如图2所示.
上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护范围的情况下,还可以做出很多变形,这些均属于本发明的保护之列。
Claims (4)
1.一种基于局部线性近似的不完备图像补齐方法,其特征在于,
包括如下步骤:
步骤一,根据不完备图像中待恢复数据建立修复模型X;所述修复模型X为具有缺失值的矩阵;
步骤二,对所述修复模型X进行初始化,对所述修复模型X缺失值进行索引,给出近邻数目k和最大迭代次数M;
步骤三,判断修复模型运算次数N是否小于迭代次数M;若满足判断条件,程序进入下一步,否则令X=Xnew终止,输出;
步骤四,对修复模型X采用局部线性近似方法估算缺失值建立新模型Xnew;
步骤五,对新模型Xnew进行||X-Xnew||2<tolerance*||X||2判断,若满足判断条件,则令X=Xnew终止,输出;否则令N=N+1返回步骤三。
2.根据权利要求1所述的一种基于局部线性近似的不完备图像补齐方法,其特征在于,所述步骤一中修复模型X通过公式(1)获得初次迭代的邻域;
其中,m是属性个数,δ(xik)表示狄垃克函数,若xik不是缺失值,δ(xik)=1;若xik是缺失值,δ(xik)=0;xik表示第k个属性下x样本的值。
3.根据权利要求1所述的一种基于局部线性近似的不完备图像补齐方法,其特征在于,所述步骤四中局部线性近似方法通过公式(2)获得修改模型重建系数Wij;
其中,XMiss表示缺失值集合,xi和xj表示矩阵X的行,Wij是xj对xi的最优线性重建系数,Nk表示xi的k近邻集合。
4.根据权利要求2所述的一种基于局部线性近似的不完备图像补齐方法,其特征在于,所述步骤四中通过公式(3)生成缺失值xip的近似估计,以此构成新模型Xnew;
。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710495129.8A CN107369137A (zh) | 2017-06-26 | 2017-06-26 | 一种基于局部线性近似的不完备图像补齐方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710495129.8A CN107369137A (zh) | 2017-06-26 | 2017-06-26 | 一种基于局部线性近似的不完备图像补齐方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107369137A true CN107369137A (zh) | 2017-11-21 |
Family
ID=60305677
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710495129.8A Pending CN107369137A (zh) | 2017-06-26 | 2017-06-26 | 一种基于局部线性近似的不完备图像补齐方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107369137A (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030053707A1 (en) * | 2001-09-13 | 2003-03-20 | Bhattacharjya Anoop K. | Techniques for scratch and date removal from scanned film |
CN102938142A (zh) * | 2012-09-20 | 2013-02-20 | 武汉大学 | 基于Kinect的室内LiDAR缺失数据填补方法 |
CN104050631A (zh) * | 2013-11-25 | 2014-09-17 | 中国科学院上海应用物理研究所 | 一种低剂量ct图像重建方法 |
CN104461772A (zh) * | 2014-11-07 | 2015-03-25 | 沈阳化工大学 | 一种对缺失数据进行恢复处理的方法 |
-
2017
- 2017-06-26 CN CN201710495129.8A patent/CN107369137A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030053707A1 (en) * | 2001-09-13 | 2003-03-20 | Bhattacharjya Anoop K. | Techniques for scratch and date removal from scanned film |
CN102938142A (zh) * | 2012-09-20 | 2013-02-20 | 武汉大学 | 基于Kinect的室内LiDAR缺失数据填补方法 |
CN104050631A (zh) * | 2013-11-25 | 2014-09-17 | 中国科学院上海应用物理研究所 | 一种低剂量ct图像重建方法 |
CN104461772A (zh) * | 2014-11-07 | 2015-03-25 | 沈阳化工大学 | 一种对缺失数据进行恢复处理的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yao et al. | Dr2-net: Deep residual reconstruction network for image compressive sensing | |
Zhang et al. | Self-supervised image denoising for real-world images with context-aware transformer | |
Chen et al. | Denoising hyperspectral image with non-iid noise structure | |
CN109410146A (zh) | 一种基于Bi-Skip-Net的图像去模糊算法 | |
CN110349087B (zh) | 基于适应性卷积的rgb-d图像高质量网格生成方法 | |
CN101895297A (zh) | 一种面向压缩感知的块稀疏信号重构方法 | |
CN104199627B (zh) | 基于多尺度在线字典学习的可分级视频编码系统 | |
CN112291570B (zh) | 一种基于轻量级可形变卷积神经网络的实时视频增强方法 | |
CN110598673A (zh) | 基于残差网络的遥感图像道路提取方法 | |
CN112712488A (zh) | 基于自注意力融合的遥感图像超分辨率重建方法 | |
CN105590296B (zh) | 一种基于双字典学习的单帧图像超分辨率复原方法 | |
Yue et al. | IENet: Internal and external patch matching ConvNet for web image guided denoising | |
CN113506224A (zh) | 基于多尺度生成对抗网络的图像修复方法 | |
CN111861886A (zh) | 一种基于多尺度反馈网络的图像超分辨率重建方法 | |
Shi et al. | Multi-scale deep networks for image compressed sensing | |
Yin et al. | CSformer: Cross-scale features fusion based transformer for image denoising | |
Wang et al. | SCA-Net: multiscale contextual information network for building extraction based on high-resolution remote sensing images | |
Wu et al. | Multi-scale non-local attention network for image super-resolution | |
CN107369137A (zh) | 一种基于局部线性近似的不完备图像补齐方法 | |
CN105279742A (zh) | 一种快速的基于分块噪声能量估计的图像去噪方法 | |
CN117115605A (zh) | 一种彩色可见光与深度图像显著性目标检测方法 | |
Wang et al. | A novel method for point cloud completion: Adaptive region shape fusion network | |
CN102156967B (zh) | 基于多尺度的局部的图像内插方法 | |
CN112991181B (zh) | 一种基于反应扩散方程的图像超分辨率重建方法 | |
CN103236041A (zh) | 一种基于Contourlet变换的图像超分辨率重建方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20171121 |
|
WD01 | Invention patent application deemed withdrawn after publication |