最大熵直接迭代法,维纳滤波和盲恢复对水平方向匀速直线运动造成的模糊图像的恢复研究 - matlab图像处理 - 谷速源码
下载频道> 资源分类> matlab源码> 图像处理> 最大熵直接迭代法,维纳滤波和盲恢复对水平方向匀速直线运动造成的模糊图像的恢复研究

标题：最大熵直接迭代法,维纳滤波和盲恢复对水平方向匀速直线运动造成的模糊图像的恢复研究
分享到：

所属分类: 图像处理 资源类型：程序源码 文件大小: 2.29 KB 上传时间: 2019-06-16 15:28:38 下载次数: 568 资源积分：1分 提 供 者: zhangsan456 max_extropy_main
内容：
最大熵直接迭代法,维纳滤波和盲恢复对水平方向匀速直线运动造成的模糊图像的恢复研究
function deblurred_image=max_entropy(blurred_image,PSF,n,mask);%最大熵直接迭代算法程序
%blurred_image　模糊化图像，即算法中的Yk；
%PSF 点扩展函数；
%deblurred_image 去模糊化后的恢复图像，即算法中的Xi；
%mask 对残损图像进行掩模操作的矩阵；
%n 迭代次数
if nargin==3
    mask=ones(size(blurred_image));
end
if nargin==2
    n=100;%迭代次数有些讲究
    mask=ones(size(blurred_image));%掩模图像支持
end
if nargin<2
    error('too little input arguments');
end
classin=class(blurred_image);
classChanged=0;
if ~isa(blurred_image,'double') 
    blurred_image=im2double(blurred_image);%完成图像数据类型转换
    classChanged=1;
end
lambda=zeros(size(blurred_image));%lambda参数初始化为0，此时熵最大，但不满足约束方程
blurred_image=blurred_image+mean2(blurred_image)*ones(size(blurred_image))/1000;%对于有些点可能为0，在后面的递推式中会出现log0，故必须保证大于0
totd=sum(sum(blurred_image));
for i=1:n
    deblurred=exp(conv2(lambda,PSF,'same'));%卷积和指数运算
    deblurred=deblurred/sum(sum(deblurred));%归一化
    reblurred=conv2(deblurred,PSF,'same');%Rk，重模糊化图像
    scale=sum(sum(reblurred.*mask))/totd;%动态归一化因子
    err=abs(blurred_image.*scale-reblurred.*mask);%误差
    entropy=-sum(sum(deblurred.*log(deblurred)));%图像熵
    if i<50
        lambda=lambda+(blurred_image*scale-reblurred).*mask;%前面50次迭代选择Gk作为修正项
    else
        reblurred=0.00001*(reblurred==0)+reblurred;
        lambda=lambda+log(blurred_image*scale./reblurred).*mask;%后面选择lnYk-lnRk作为修正项
    end
end
deblurred_image=deblurred/scale-mean2(blurred_image)/1000;%减去前面的增量
figure,imshow(deblurred_image); ('max extropy deblurred image');%显示去模糊化后的恢复图像

文件列表（点击上边下载按钮，如果是垃圾文件请在下面评价差评或者投诉）:

max_extropy_main/
max_extropy_main/max_entropy.m
max_extropy_main/max_entropy_main.m
max_extropy_main/max_extropy_readme.txt

关键词： 最大熵 模糊图像

相关推荐

matlab源码下载排行

Top_arrow
回到顶部

友情链接

matlab源码 网站源码 编程语言 网络编程 手机通讯编程 KUNAG 官网 KUNAG 电商 更多
联系方式| 版权声明| 招聘信息| 广告服务| 银行汇款| 法律顾问| 兼职技术| 付款方式| 关于我们|
网站客服网站客服 程序员兼职招聘 程序员兼职招聘
沪ICP备19040327号-3
公安备案号：沪公网安备 31011802003874号
库纳格流体控制系统（上海）有限公司 版权所有
Copyright © 1999-2014, GUSUCODE.COM, All Rights Reserved