数字图像处理作业图像运动模糊约束最小二乘方滤波MATLAB源码及实验报告

数字图像处理作业图像运动模糊&约束最小二乘方滤波MATLAB源码及实验报告

%%%%%约束最小二乘方滤波%%%%%
clear;close all;clc;
gk=imread（‘E:\学习\A+数字图像处理\作业4\ImageBlur_2.bmp‘）;
[MNK]=size（gk）;
gk=double（gk）;
g=gk（::1）*0.299+gk（::2）*0.587+gk（::3）*0.114;%Gray = R*0.299 + G*0.587 + B*0.114
%% 加入高斯噪声
%%%构造一个高斯噪声矩阵
mu=0;sigma=1;
noise=normrnd（musigmaMN）;

%%% 构造一个矩阵temp随机选出20%的像素用于加噪声
temp=zeros（MN）;%%%默认每个元素都为0，即不加噪声
amount=round（M*N*0.2）; %%%加噪像素总量
for j=1:amount
    x=randi（M11）;
    y=randi（N11）;
    temp（xy）=1;%%%将选中的像素点的值置为1，对这些像素点加噪声
end

%%%加高斯噪声
gn=g+temp.*noise;
gnoise=gn-g;
gnn=uint8（gn）;
hold on
imshow（gnn）;
imwrite（gnn‘E:\学习\A+数字图像处理\作业4\ImageBlurN.bmp‘）;

G=fftshift（myft2mn（gn））;%% 二维离散傅立叶变换
%% 退化函数
T=1; a=0.1; b=0.1;
H=zeros（MN）;
for u=1:M
    for v=1:N
        u1=u-M/2-1;%将H（uv）的（00）点移到矩阵中心
        v1=v-N/2-1

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-04-15 21:44  图像运动模糊&约束最小二乘方滤波\
     文件        1913  2018-11-18 20:12  图像运动模糊&约束最小二乘方滤波\CLSF.m
     文件         612  2018-11-18 15:52  图像运动模糊&约束最小二乘方滤波\ImageBlur.m
     文件       91078  2018-11-18 18:57  图像运动模糊&约束最小二乘方滤波\ImageBlurN.bmp
     文件      270054  2018-11-18 15:55  图像运动模糊&约束最小二乘方滤波\ImageBlur_2.bmp
     文件      270054  2018-11-18 17:40  图像运动模糊&约束最小二乘方滤波\ImageDeblur.bmp
     文件         473  2018-11-18 16:11  图像运动模糊&约束最小二乘方滤波\myft2mn.m
     文件         533  2018-11-18 16:11  图像运动模糊&约束最小二乘方滤波\myift2mn.m
     文件      237170  2019-04-15 21:44  图像运动模糊&约束最小二乘方滤波\实验报告.pdf