思路:
1转换为二进制图像并覆盖。
2生成名为的翻转图像。
3使用 ICM 对翻转的图像进行降噪。
代码:

结果:

注意知识点:

1原图:

quzao3.jpg

我试着用其他图做了这个实验,但是有灰色区域就会出现问题。

2最终成像时:给图片加标题,如语句:prop = font_manager.FontProperties(fname=r'E:\学习代码\ttf字体\华康娃娃体.ttf'),其中的华康娃娃体.ttf需要下载,用的时候指定路径即可。

3 ICM我认为是这个程序的重点!!ICM 算法是计算最小能量函数的经典算法。

4接下来就是用matlab去实现降噪

close all;
clear all;
I=imread('quzao3.jpg');
I1=im2bw(I);
subplot(2,2,1)
imshow(I1);
title('原图像');
J = imnoise(I,'salt & pepper', 0.2);
J1=im2bw(J);
subplot(2,2,2)
imshow(J1);
title('噪声图')
 
Y=ones(size(J1));
Y(J1==0)=-1;
[m,n]=size(Y);
X=Y;h=0;beta=3.5;eta=.1;
while 1
    tot=0
    for i=2:1:m-1
        for j=2:1:n-1
            temp=X(i,j);
            X(i,j)=-1;%根据定义计算势函数
        E1=h*X(i,j)-beta*X(i,j)*(X(i-1,j)+X(i+1,j)+X(i,j-1)+X(i,j+1))-eta*X(i,j)*Y(i,j);
           X(i,j)=1;
           E2=h*X(i,j)-beta*X(i,j)*(X(i-1,j)+X(i+1,j)+X(i,j-1)+X(i,j+1))-eta*X(i,j)*Y(i,j);
           if E1<E2
                X(i,j)=-1;
            else
                X(i,j)=1;
            end   
            if temp~=X(i,j)
                tot=tot+1;
            end
        end
    end
    if tot<1
        break;
    end
end
 
J2=X;
J2(X==-1)=0;
subplot(2,2,3)
imshow(J2);
title('mrf降噪结果')

成像结果:

这个代码重点就是势函数那一部分。

!!!马尔科夫随机场:重点在于变量的联合密度函数以及可以与吉普斯随机场对照学习。

MRF马尔科夫链在图像去噪中的应用(python+matlab)相关推荐

  1. 第十四课.马尔科夫链

    目录 马尔科夫链 马尔科夫链的平稳分布 马尔科夫链进入稳态的转移过程 马尔科夫链稳态的意义 稳态判定:细致平稳条件 马尔科夫链 随机过程指的是一个随机变量序列:{Xt}\left\{X_{t}\rig ...

  2. 马尔科夫链与MCMC方法

    马尔科夫链概述 基本思想: 过去所有的信息都已经被保存到了现在的状态,基于现在就可以预测未来. Example: 假如每天的天气是一个状态的话,那今天是不是晴天只依赖于昨天的天气,而和前天的天气没有任 ...

  3. 求解平稳分布matlab,随机过程课程设计--应用马尔科夫链的平稳分布预测空调市场的占有率.doc...

    课程名称: <随机过程> 课程设计(论文) 题 目: 应用马尔科夫链的平稳分 布预测市场占有率 学 院: 理学院 专 业: 数学与应用数学 班 级: 15-1 学 生 姓 名: 邹光睿 学 ...

  4. 马尔科夫链和马尔科夫链蒙特卡洛方法

    前言 译自:<Training Restricted Boltzmann Machines: An Introduction > 马尔科夫链在RBM的训练中占据重要地位,因为它提供了从复杂 ...

  5. 结对编程-马尔科夫链作业成绩

    结对编程-马尔科夫链作业成绩 作业博客:http://www.cnblogs.com/vertextao/p/6881960.html 问题 本次作业仅有三位同学按时提交,两位同学延迟提交,出人意料, ...

  6. 第十五课.马尔科夫链蒙特卡洛方法

    目录 M-H采样 Metropolis-Hastings采样原理 M-H采样步骤 Gibbs方法 Gibbs核心流程 Gibbs采样的合理性证明 Gibbs采样实验 在 第十四课中讲述了马尔科夫链与其 ...

  7. MCMC+马尔科夫链蒙特卡罗

    MCMC+马尔科夫链蒙特卡罗 为了解决什么问题,所以出现了这一种方法? 后来又因为出现了什么新情况,所以产生了对应的变种?

  8. 从马尔科夫链到吉布斯采样与PageRank

    马尔科夫链表示state的链式关系,下一个state只跟上一个state有关. 吉布斯采样通过采样条件概率分布得到的样本点,近似估计概率分布P(z)P(z).PageRank通过节点间的连接,估计节点 ...

  9. 13张动图助你彻底看懂马尔科夫链、PCA和条件概率!

    来源:新智元 本文共2100字,建议阅读9分钟. 本文用可视化的方式来解释抽象的理论概念,使这些抽象概念变得生动而立体! [ 导读 ]马尔科夫链.主成分分析以及条件概率等概念,是计算机学生必学的知识点 ...

最新文章

  1. 查看Flink的Job Graph时的问题
  2. Unity重置Animator到初始状态和重复播放同一个Animation
  3. jquery获取html代码怎么写,jQuery从html代码中获取对应标签的写法
  4. 使用DotNetty编写跨平台网络通信程序
  5. 【OpenCV 例程200篇】91. 高斯噪声、瑞利噪声、爱尔兰噪声
  6. js designMode contentEditable 编辑在线网页
  7. c# MODBUS协议 上位机 (转载)
  8. Android中的EditText默认时不弹出软键盘的方法
  9. python scrapy框架爬虫当当图书网
  10. OEM XP 让你的盗版xp像品牌机一样华丽
  11. 强化学习推荐系统_推荐人系统:价值调整,强化学习和道德规范
  12. 【2019牛客暑期多校训练营(第八场)E】Explorer【线段树分治维护图连通】
  13. win10磁盘管理 磁盘分区和合并
  14. python hist2d_具有二维数组的Python Matplotlib Hist2d
  15. 使用开源项目DCRM搭建Cydia源服务器
  16. 尚硅谷大厂面试题第二季周阳主讲整理笔记【涨薪必看】
  17. ffmpeg每隔1秒抽取视频的1帧画面
  18. 资源 步金国 的高手
  19. I.MX6ULL ARM驱动开发---设备树下的LED驱动实验
  20. 浙大MBA复试上岸经历分享——万事俱备复试必得

热门文章

  1. 兄弟连Linux学习笔记之文件处理命令
  2. STM32代码空间不够问题的解决
  3. 完美国际地图文件在服务器那个文件夹,完美国际------【大地图安全区】------修改教程...
  4. 《完美世界国际版》客户端下载
  5. 教学管理软件——教学评价系统
  6. 汇编学习教程:CPU内的小秘密
  7. 1106-材料 聚氨酯,FR4,PC+GF,PC+ABS,钕铁硼,mylar,3M,
  8. OneThink 忘记密码 OneThink重置密码
  9. INJECT_EVENTS权限与输入法的全屏手写
  10. kaldi学习 - egs/yesno —— 数据准备(二)