clf;
clc,clear,close all;I=imread('peppers.png'); % 讀入圖像
I1=rgb2gray(I); % 圖像灰度轉換
figure;imshow(I);title('Original Image');
figure;imshow(I1);title('Original Grayimage');
PSF=fspecial('gaussian',7,10); % 設置高斯低通濾波器
Blurred=imfilter(I1,PSF,'symmetric','conv'); % 圖像卷積計算,symmetric圖像大小通過鏡像反射其邊界來擴展
figure;imshow(Blurred); title('Blurred Image');
UNDERPSF = ones(size(PSF)-4); % 欠尺寸PSF，比原PSF小4
[J1,P1] = deconvblind(Blurred,UNDERPSF); % J1表示圖形，P1表示PSF模型
figure;imshow(J1);title('Deblurring with Undersized PSF'); % 欠尺寸PSF,圖像復原,比原PSF小4
OVERPSF = padarray(UNDERPSF,[4 4],'replicate','both'); % 過尺寸PSF,比原PSF大4
[J2,P2] = deconvblind(Blurred,OVERPSF); % padarray填充圖像或填充數組
figure;imshow(J2);title('Deblurring with Oversized PSF'); % 過尺寸圖像復原，比原PSF大4
INITPSF = padarray(UNDERPSF,[2 2],'replicate','both'); % 全尺寸PSF
[J3,P3] = deconvblind(Blurred,INITPSF); % 圖像復原，原尺寸PSF
figure;imshow(J3);title('Deblurring with INITPSF'); % 設置圖像標題
figure;
subplot(221);imshow(PSF,[],'InitialMagnification','fit');title('True PSF'); % 圖像標題顯示
subplot(222);imshow(P1,[],'InitialMagnification','fit');title('Reconstructed Undersized PSF'); % 圖像標題顯示
subplot(223);imshow(P2,[],'InitialMagnification','fit');title('Reconstructed Oversized PSF'); % 圖像標題顯示
subplot(224);imshow(P3,[],'InitialMagnification','fit');title('Reconstructed true PSF'); % 圖像標題顯示
WEIGHT = edge(I1,'sobel',0.1); % 圖像邊緣顯示
se = strel('disk',2);
WEIGHT = 1-double(imdilate(WEIGHT,se)); % imdilate函數用於對圖像實現膨脹操作
WEIGHT([1:3 end-[0:2]],:) = 0; % 上下邊緣
WEIGHT(:,[1:3 end-[0:2]]) = 0; % 左右邊緣
figure;imshow(WEIGHT);title('Weight array');
[J,P] = deconvblind(Blurred,INITPSF,30,[],WEIGHT); % 圖像復原
figure;
imshow(J);
title('Deblurred Image');
P1 = 2;
P2 = 2;
FUN = @(PSF) padarray(PSF(P1+1:end-P1,P2+1:end-P2),[P1 P2]); % padarray填充圖像或填充數組。
[JF,PF] = deconvblind(Blurred,OVERPSF,30,[],WEIGHT,FUN); % 圖像復原
figure; % 新建窗口
imshow(JF); % 顯示處理後的圖像
title('Deblurred Image'); % 去模糊圖像