模糊噪声图像的建模

%checkerboard产生测试板图像，第一个参数为每个正方形一边的像素数，第二个为参数行数，第三个为列数
f=checkerboard(8);%产生一个一边为8个正方形的测试板
subplot(2,2,1),imshow(f),title('原图像')

PSF=fspecial('motion',7,45)%用fspecial处理图像的运动模糊，PSF为空间滤波器
gb=imfilter(f,PSF,'circular');%减少边界效应
subplot(2,2,2),imshow(gb);title('运动模糊图像')

noise=imnoise(zeros(size(f)),'gaussian',0,0.001);%高斯噪声
subplot(2,2,3),imshow(noise);title('运动模糊图像')

g=gb+noise;%添加运动模糊的高斯噪声
subplot(2,2,4),imshow(pixeldup(g,8),[]);title('运动模糊+高斯噪声图像')%大图像运算过程过慢，选用小图像来节省时间