根据运动模型生成运动模糊图像

Myth丶恋晨 2023-07-18 06:20 42阅读 0赞

## 根据运动模型生成运动模糊图像 ##

**具体的步骤：先用程序二生成txt文件——执行主函数即可（参数选择）**  
1.确定图像退化模型  
2.生成退化图像  
（ps:以下是用到的函数）  
一、函数

%%%%%%%motionblur函数
    function h=motionblur(dirangle,winsize)
    if nargin<2
        winsize=9;
    end
    h=zeros(winsize);
    ext=(winsize-1)/2;
    if (abs(abs(dirangle)-90) >=45) & (abs(abs(dirangle)-270)>=45),
        slope=tan(dirangle*pi/180); 
        rloc=round(slope*[-ext:ext]);
        for i=1:winsize,
            h(ext-rloc(i)+1,i)=1;
        end
    else
        slope=cot(dirangle*pi/180);
        cloc=round(slope*[-ext:ext]);
        for i=1:winsize,
            h(i,ext-cloc(i)+1)=1;
        end
    end

二、图片生成txt文件格式的程序

%%%%%图像数据%%%
    rgb = imread('cdl.jpg'); 
    % 打开文件用于保存结果
    subplot(121)
    image(rgb)%原始图
    %%%fid = fopen('JPEG2TXT.txt', 'wt+');%%如果没有找到txt文件则用下面这句
    fid = fopen('C:\Users\huawei\Desktop\JPEG2TXT.txt', 'wt+');%%这是我的电脑桌面位置，仿改即可
    if fid == -1, return, end % 写文件，按照红绿蓝顺序，每种颜色写成矩阵形式
    for i = 1 : 3   
        for r = 1 : size(rgb, 1)       
            fprintf(fid, '%4i', double(rgb(r, :, i)));       
            fprintf(fid, '\n');    
        end
        fprintf(fid, '\n');
    end
    fclose(fid);
    %%打开txt并转成图片
    rgb1=load('JPEG2TXT.txt','-ascii');
    [r,c]=size(rgb1);
    rgb2=permute(reshape(rgb1,[r/3, 3, c]),[1 3 2]);
    subplot(122)
    image(rgb2/255)%重新用txt 打开的图

三、主程序

%%根据运动模型生成运动模糊图像
    load JPEG2TXT.txt;
    [m,n]=size(JPEG2TXT);
    winsize=input('Blur operator window size (an odd number, default= 9):');
    if isempty (winsize), winsize=9;
        elseif rem (winsize,2)==0,
            winsize=winsize+1; 
            disp(['Use odd number for window size =' int2str(winsize)])
    end
    disp(['1.Linear motion blur;'])
    chos=input('Enter a number to choose type of blur applied(default=1):');
    if isempty (chos),chos=1;end
    if chos==1,
        dirangle=input('Bluring direction(an angle in degrees, default=45)=');
        if isempty(dirangle)
            dirangle=45;
        end
        h=motionblur(dirangle,winsize);%%%调用运动模糊函数生成模糊模型
    end
    %%根据模糊模型生成模糊图像
    F=fft2(JPEG2TXT);
    Hmat=fft2(h,372,200);%%%%%%%%%%%%%%特别注意，这里的维度要统一，即图像的像素大小
    Gmat=F.*Hmat;
    g=ifft2(Gmat);
    figure(1)
    subplot(121),imagesc(JPEG2TXT),colormap('gray'),title('originalimage')
    subplot(122),imagesc(abs(g)),colormap('gray'),title('blurredimage')
    figure(2),
    subplot(212),imagesc(log(1+abs(Gmat))),colormap('gray'),title('blurring filter')
    subplot(211),imagesc(h),colormap('gray'),title('blurring filtermask')

五、采用顺序统计滤波器对图像进行滤波主程序

%%采用顺序统计滤波器对图像进行滤波
    img=rgb2gray(imread('cdl.png'));
    figure;imshow(img);
    img=double(imnoise(img,'salt & pepper',0.1));
    figure,imshow(img,[]);
    %%采用均值滤波
    N=5;
    h=fspecial('average',N);
    I=filter2(h,img);
    figure,imshow(I,[])
    %%中值滤波
    I=medfilt2(img,[N N]);
    figure,imshow(I,[])
    %%最大值滤波
    I=ordfilt2(img,N*N,true(N));
    figure,imshow(I,[])
    %%最小值滤波
    I=ordfilt2(img,1,true(N));
    figure,imshow(I,[])