低阶秩序图像修复模型
低阶秩图像修复模型,假设图像的数据维度秩(rank)较低,那么,可以利用该性质实现图像修复,具体模型如下:
通过将迹范数转化为核范数,将非凸优化问题转化为凸优化问题,如下:
通过迭代求解上述凸优化问题可以实现图像修复:
% demo2.m
% Date: 2015/6/12
% Author: HSW
% HARBIN INSTITUTE OF TECHNOLOGY
% set matlab
close all;
clear all;
clc;
% add path
addpath(genpath('DamagedImage\'));
addpath(genpath('ResultsImage\'));
addpath(genpath('Solvers\'));
addpath(genpath('TestImage\'));
% Set True rank of Image
Rank = 40;
% read image
Img = imread('TestImage2\Lowrank.bmp');
Img = imresize(Img,[256,256],'bicubic');
% Img = imread('TestImage2\mat.bmp');
% Img = imread('TestImage2\rattan.tif');
% Img = imread('TestImage2\Build.bmp'); %可以达到95%
% Img = imread('TestImage\barba256.png');
% Img = imread('D82.gif');
% Img = imresize(Img,[256,256],'bicubic');
if size(Img,3) == 3
OrigImg = rgb2gray(Img);
elseif size(Img,3) == 4
OrigImg = mean(Img,3);
else
OrigImg = Img;
end
OrigImg = OrigImg';
OrigImg = double(OrigImg);
[U,S,V] = svd(OrigImg);
OrigRank = length(find(S ~= 0));
S(Rank+1:end,:) = 0; %注意这里
LowRankImg = U*S*V';
figure;
subplot(1,3,1)
imshow(OrigImg/255,[]);
title(['Image of Rank = ', num2str(OrigRank)]);
subplot(1,3,2);
imshow(LowRankImg/255,[]);
title(['Image of Rank = ',num2str(Rank)]);
subplot(1,3,3);