Python实现直方图均衡基本原理解析


Posted in Python onAugust 08, 2019

1. 基本原理

通过一个变换,将输入图像的灰度级转换为`均匀分布`,变换后的灰度级的概率密度函数为

$$P_s(s) = \frac{1}{L-1}$$

直方图均衡的变换为

$$s = T(r) = (L-1)\int_0^r {P_r(c)} \,{\rm d}c $$

$s$为变换后的灰度级,$r$为变换前的灰度级$P_r(r)$为变换前的概率密度函数2. 测试结果

Python实现直方图均衡基本原理解析

图源自skimage

3.代码

import numpy as np
def hist_equalization(input_image):
  '''
  直方图均衡(适用于灰度图)
  :param input_image: 原图像
  :return: 均衡后的图像
  '''
  output_imgae = np.copy(input_image) # 输出图像,初始化为输入
  input_image_cp = np.copy(input_image) # 输入图像的副本
  m, n = input_image_cp.shape # 输入图像的尺寸(行、列)
  pixels_total_num = m * n # 输入图像的像素点总数
  input_image_grayscale_P = [] # 输入图像中各灰度级出现的概率,亦即输入图像直方图
  # 求输入图像中各灰度级出现的概率,亦即输入图像直方图
  for i in range(256):
    input_image_grayscale_P.append(np.sum(input_image_cp == i) / pixels_total_num)
  # 求解输出图像
  t = 0        # 输入图像的灰度级分布函数F
  for i in range(256):
    t = t + input_image_grayscale_P[i]
    output_imgae[np.where(input_image_cp == i)] = 255 * t
  return output_imgae

4. 数学证明目标变换

  • $$S = T(r) = (L-1)\int_0^rp_r(w)dw$$
  • $T(r)$为严格单调函数,可保证反映射时,消除二义性$p_r(w)$为源图像归一化后的直方图

4.1 假定

  • 图像灰度级为:$[0, L-1]$
  • 源图像中,$k$灰度级的像素个数:$n_k$
  • 源图像像素总数:$n$原图像直方
  • 图$h(r_k) = n$4.2 归一化后的直方图

$$p(r_k) = n_k / n$$

$p(r_k)$即为灰度级$r_k$在源图像中出现的概率估计

4.3 证明

概率密度函数的积分为分布函数,即对分布函数的导数为概率密度函数。

因为$p_r(r)$与$T(r)$已知,则由

$$\frac{{\rm d}r}{{\rm d}S} = \frac{p_s(s)}{p_r(r)}$$

又因为

$$S = T(r)$$


$$\frac{{\rm d}S}{{\rm d}r} = \frac{T(r)}{r}$$

联立上三式及目标变换

$$S = T(r) = (L-1)\int_0^rp_r(w)dw$$

可得

$$p_s(s) = \frac{1}{L-1}$$

故,这意味着变换之后的图像的灰度级为均匀分布,证毕。

总结

以上所述是小编给大家介绍的Python实现直方图均衡基本原理解析,希望对大家有所帮助,如果大家有任何疑问请给我留言,小编会及时回复大家的。在此也非常感谢大家对三水点靠木网站的支持!
如果你觉得本文对你有帮助,欢迎转载,烦请注明出处,谢谢!

Python 相关文章推荐
python正则表达式修复网站文章字体不统一的解决方法
Feb 21 Python
用Python代码来绘制彭罗斯点阵的教程
Apr 03 Python
使用Python的Bottle框架写一个简单的服务接口的示例
Aug 25 Python
Django 导出 Excel 代码的实例详解
Aug 11 Python
python 匹配url中是否存在IP地址的方法
Jun 04 Python
python 正确保留多位小数的实例
Jul 16 Python
对pandas的行列名更改与数据选择详解
Nov 12 Python
Numpy截取指定范围内的数据方法
Nov 14 Python
python opencv 读取本地视频文件 修改ffmpeg的方法
Jan 26 Python
django 使用全局搜索功能的实例详解
Jul 18 Python
python协程gevent案例 爬取斗鱼图片过程解析
Aug 27 Python
浅谈Pytorch中的自动求导函数backward()所需参数的含义
Feb 29 Python
python获取指定日期范围内的每一天,每个月,每季度的方法
Aug 08 #Python
Python 调用 Outlook 发送邮件过程解析
Aug 08 #Python
python判断自身是否正在运行的方法
Aug 08 #Python
Python 日期区间处理 (本周本月上周上月...)
Aug 08 #Python
python各类经纬度转换的实例代码
Aug 08 #Python
Python3.0 实现决策树算法的流程
Aug 08 #Python
Python使用ffmpy将amr格式的音频转化为mp3格式的例子
Aug 08 #Python
You might like
php中使用explode查找某个字符是否存在的方法
2011/07/12 PHP
功能强大的php分页函数
2016/07/20 PHP
ajax 文件上传应用简单实现
2009/03/03 Javascript
javascript类继承机制的原理分析
2009/09/12 Javascript
对setInterval在火狐和chrome切换标签产生奇怪的效果之探索,与解决方案!
2011/10/29 Javascript
目前流行的JavaScript库的介绍及对比
2013/09/29 Javascript
下拉列表select 由左边框移动到右边示例
2013/12/04 Javascript
javascript动态创建及删除元素的方法
2014/12/22 Javascript
jQuery中clearQueue()方法用法实例
2014/12/29 Javascript
jQuery基于ajax实现页面加载后检查用户登录状态的方法
2017/02/10 Javascript
微信小程序实战之自定义抽屉菜单(7)
2017/04/18 Javascript
在JS中如何把毫秒转换成规定的日期时间格式实例
2017/05/11 Javascript
AngularJS常见过滤器用法实例总结
2017/07/06 Javascript
JavaScript基于面向对象实现的猜拳游戏
2018/01/03 Javascript
解决vue页面DOM操作不生效的问题
2018/03/17 Javascript
基于AngularJS拖拽插件ngDraggable.js实现拖拽排序功能
2019/04/02 Javascript
你了解vue3.0响应式数据怎么实现吗
2019/06/07 Javascript
[01:28:56]2014 DOTA2华西杯精英邀请赛 5 24 CIS VS DK
2014/05/26 DOTA
[00:52]DOTA2齐天大圣预告片
2016/08/13 DOTA
使用python实现链表操作
2018/01/26 Python
在Python中实现shuffle给列表洗牌
2018/11/08 Python
在Python中关于使用os模块遍历目录的实现方法
2019/01/03 Python
PyQt5实现简易计算器
2020/05/30 Python
Python 爬虫实现增加播客访问量的方法实现
2019/10/31 Python
Python调用shell cmd方法代码示例解析
2020/06/18 Python
详解Python Celery和RabbitMQ实战教程
2021/01/20 Python
轻松制作精彩视频:Animoto
2018/09/19 全球购物
英国钻石公司:British Diamond Company
2020/02/16 全球购物
加拿大大码女装购物网站:Penningtons
2020/12/26 全球购物
合作协议书怎么写
2014/04/18 职场文书
企业趣味活动方案
2014/08/21 职场文书
户籍证明模板
2014/09/28 职场文书
房屋租赁意向书范本
2015/05/09 职场文书
七年级英语教学反思
2016/02/15 职场文书
Python面向对象之成员相关知识总结
2021/06/24 Python
mysql全面解析json/数组
2022/07/07 MySQL