【图像重建】正则化图像超分辨重建【含Matlab源码 1882期】

news2024/11/24 16:24:05

⛄一、正则化图像超分辨重建简介

图像超分辨率重建的非局部正则化模型与算法研究

利用图像非局部不连续性测度的概念,建立了面向图像超分辨的非局部正则化能量泛函和相应的变分框架.理论分析了该框架与目前关于双边滤波等一类广义邻域滤波器和经典的变分偏微分方程模型之间的联系.推导了该非局部泛函约束的变分模型最优解满足的积分形式欧拉-拉格朗日方程,并研究了其最速下降流满足的若干重要性质.基于图理论,设计了图像超分辨的自适应加权迭代算法.最后通过选择不同势函数的非局部正则化泛函进行图像去噪、去马赛克效应和图像超分辨处理,性能分析表明:相同势函数下,非局部正则化方法优于同类局部正则化方法,峰值信噪比提高0.5~1.0dB.

⛄二、部分源代码

clear
close all
clc

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Demo file for deconvtv
% Image ‘salt and pepper’ noise removal
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Prepare images
f_orig = im2double(imread(‘./data/lena_color_512.png’));
[rows,cols,frames] = size(f_orig);
H = fspecial(‘gaussian’, [9 9], 2);
g = imfilter(f_orig, H, ‘circular’);
g = imnoise(g, ‘salt & pepper’, 0.05);

% Setup parameters (for example)
opts.rho_r = 5;
opts.rho_o = 100;
opts.beta = [1 1 0];
opts.print = true;
opts.alpha = 0.7;
opts.method = ‘l1’;

% Setup mu
mu = 20;

% Main routine
tic
out = deconvtv(g, H, mu, opts);
toc

PSNR_output = psnr(out.f,f_orig);
fprintf(‘\nPSNR = %3.2f dB \n’, PSNR_output);
% Display results
figure(1);
imshow(g);
set(get(gca, ‘Title’), ‘String’, ‘input’);

figure(2);
imshow(out.f);
set(get(gca, ‘Title’), ‘String’, ['PSNR = ', num2str(PSNR_output),‘dB’]);