灰度直方图
图像处理6灰度直方图和直方图均衡化
图像处理6灰度直⽅图和直⽅图均衡化灰度直⽅图介绍灰度直⽅图（Gray histogram）是关于灰度级分布的函数，是对图像中灰度级分布的统计。
灰度直⽅图是将数字图像中的所有像素，按照灰度值的⼤⼩，统计其出现的频率。
灰度直⽅图是灰度级的函数，它表⽰图像中具有某种灰度级的像素的个数，反映了图像中某种灰度出现的频率。
如果将图像总像素亮度（灰度级别）看成是⼀个随机变量，则其分布情况就反映了图像的统计特性，这可⽤probability density function （PDF）来刻画和描述，表现为灰度直⽅图。
实现以下代码便于理解灰度直⽅图的计算，其中histogram函数是基于numpy简化的，运⾏结果如下。
# coding: utf8from skimage import dataimport matplotlib.pyplot as pltimport numpy as npdef histogram(a, bins=10, range=None):"""Compute the histogram of a set of data."""import numpy as npfrom numpy.core import linspacefrom numpy.core.numeric import (arange, asarray)# 转成⼀维数组a = asarray(a)a = a.ravel()mn, mx = [mi + 0.0 for mi in range]ntype = np.dtype(np.intp)n = np.zeros(bins, ntype)# 预计算直⽅图缩放因⼦norm = bins / (mx - mn)# 均分，计算边缘以进⾏潜在的校正bin_edges = linspace(mn, mx, bins + 1, endpoint=True)# 分块对于⼤数组可以降低运⾏内存，同时提⾼速度BLOCK = 65536for i in arange(0, len(a), BLOCK):tmp_a = a[i:i + BLOCK]tmp_a_data = tmp_a.astype(float)# 减去Range下限，乘以缩放因⼦,向下取整tmp_a = tmp_a_data - mntmp_a *= normindices = tmp_a.astype(np.intp)# 对indices标签分别计数，标签等于bins减⼀indices[indices == bins] -= 1n += np.bincount(indices, weights=None,minlength=bins).astype(ntype)return n, bin_edgesif__name__ =="__main__":img=data.coffee()fig = plt.figure()f1 = fig.add_subplot(141)f1.imshow(img)f1.set_title("image")f2 = fig.add_subplot(142)arr=img.flatten()n, bins, patches = f2.hist(arr, bins=256, facecolor='red')f2.set_title("plt_hist")f3 = fig.add_subplot(143)hist, others = np.histogram(arr, range=(0, arr.max()), bins=256)f3.plot(others[1:],hist)f3.set_title("np_hist1")f4 = fig.add_subplot(144)hist, others = histogram(arr, range=(0, arr.max()), bins=256)f4.plot(others[1:], hist)f4.set_title("np_hist2")plt.show()关于bincount函数，可以参考Xurtle的博⽂https:///xlinsist/article/details/51346523bin的数量⽐x中的最⼤值⼤1，每个bin给出了它的索引值在x中出现的次数。
灰度直方图的用途
灰度直方图的用途灰度直方图（grayscale histogram）是一种统计图表，用于表示图像中每个灰度级别的像素数量。
它将图像的灰度级别作为横坐标，像素数量作为纵坐标，能够直观地展示图像中各个灰度级别像素的分布情况。
灰度直方图在计算机视觉和图像处理领域有着广泛的应用，以下将详细介绍它的用途。
1. 图像增强与调整灰度直方图可以帮助我们快速了解图像的亮度分布情况。
通过观察图像的灰度直方图，我们可以判断图像的亮度范围，并根据需要进行图像增强和调整。
例如，如果图像的灰度分布集中在较低的亮度范围，我们可以通过直方图均衡化等算法来增强图像的对比度，使得图像更加清晰明亮。
2. 图像分割与阈值处理利用灰度直方图，我们可以选择适当的阈值来进行图像分割。
图像分割是将图像划分为几个具有独立意义的区域，常用于目标检测和图像识别等应用。
通过分析灰度直方图可以确定一个或多个阈值，将图像分割成不同的区域，从而提取出感兴趣的目标。
3. 图像质量评估灰度直方图可以用于图像质量的评估和比较。
不同的图像质量可能对应着不同的灰度分布情况，通过对比不同图像的灰度直方图，我们可以直观地了解图像的细节信息和对比度。
在图像压缩、图像传输和图像复原等应用中，可以通过灰度直方图的对比来评估图像的质量，并确定是否需要使用图像增强或去噪等算法来改善图像的质量。
4. 图像匹配与检索利用灰度直方图可以进行图像的匹配与检索。
在图像数据库中，通过计算图像的灰度直方图，我们可以建立索引来加速图像的匹配和搜索。
通过对比查询图像与数据库中图像的灰度直方图相似性，可以找到与查询图像最相似的图像，从而实现图像的检索和识别。
5. 图像分类与识别灰度直方图也可用于图像分类与识别。
通过将图像的灰度直方图作为特征向量，我们可以训练分类器或者使用统计模型来对图像进行分类和识别。
例如，基于灰度直方图的人脸识别算法可以通过计算人脸图像的灰度直方图来进行人脸匹配和识别。
6. 图像分析与检测利用灰度直方图还可以进行图像的分析与检测。
第五章 遥感图像处理—图像增强
特征；其余三个分量与地物特征没有明确的对应关系。
七、多元信息复合
遥感图像信息融合(Fusion)是将多源遥感数据在统一的 地理坐标系中，采用一定的算法生成一组新的信息或合
其中：
k ( g 'max g 'min ) /( gmax gmin ) 255/ 52 4.9
b g 'ij kgij 0 49 49
2、非线性拉伸
（1）指数变换
xb be
（2）对数变换
axa
c
xb b度进行分层，每一层所包含的亮度值范围可以不
同。
图像密度分割原理可以按如下步骤进行：
（1）求图像的极大值dmax和极小值dmin； （2）求图像的密度区间ΔD = dmax-dmin + 1； （3）求分割层的密度差Δd =ΔD／n ，其中 n为需分割的层数；
（4）求各层的密度区间；
（5）定出各密度层灰度值或颜色。
减法运算可以增加不同地物间光谱反射率以及在 两个波段上变化趋势相反时的反差。不同时相同 一波段图像相减时，可以提取波段间的变化信息。
T M 4 影 像
T M 3 影 像
TM4-TM3影像
87 年 影 像
92 年 影 像 变化监测结果影像
（二）加法运算
B= i /m
i=1 m
加法运算可以加宽波段，如绿色波段和红色波 段图像相加可以得到近似全色图像；而绿色波 段，红色波段和红外波段图像相加可以得到全 色红外图像。
-1 -2 -1 0 0 0 1 2 1 1 2 0 -2 1 0 -1
第三章.图像灰度直方图变换
第三章图像灰度直方图变换在数字图像处理中，灰度直方图是最简单且最有用的工具，可以说，对图像的分析与观察直到形成一个有效的处理方法，都离不开直方图。
直方图的定义：一个灰度级别在范围[0，L-1]的数字图象的直方图是一个离散函数p(rk)= nk/nn 是图象的像素总数，nk是图象中第k个灰度级的像素总数，rk 是第k个灰度级，k = 0,1,2,…,L-直方图的性质1）灰度直方图只能反映图像的灰度分布情况，而不能反映图像像素的位置，即丢失了像素的位置信息。
2）一幅图像对应唯一的灰度直方图，反之不成立。
不同的图像可对应相同的直方图。
直方图的应用：用来判断图像量化是否恰当灰度变换一、对比度展宽的目的：是一点对一点的灰度级的影射。
设新、旧图的灰度级分别为g 和f，g和f 均在[0，255]间变化。
目的：将人所关心的部分强调出来。
对比度展宽方法：二、灰级窗：只显示指定灰度级范围内的信息。
如: α=γ=0三、灰级窗切片：只保留感兴趣的部分，其余部分置为0。
直方图均衡化算法：设f、g分别为原图象和处理后的图像。
求出原图f的灰度直方图，设为h。
h为一个256维的向量。
求出图像f的总体像素个数Nf=m*n （m,n分别为图像的长和宽）计算每个灰度级的像素个数在整个图像中所占的百分比。
hs(i)=h(i)/Nf (i=0,1, (255)3）计算图像各灰度级的累计分布hp。
4）求出新图像g的灰度值。
作业1. 在图像灰度变换处理中，请总结出线性变换，非线性变换的适应性及各自的特点？. 已知一幅图像为：∑==ikkhihp)()(255,...,2,1=i⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=22789321227881112388712439881228291010636921001001073910101002552547120025520010022525551f请对其进行灰度直方图的均衡化处理。
图像增强—灰度变换及直方图均衡化试验目的试验原理及知识点
图像增强—灰度变换及直方图均衡化一、实验目的1、了解图像增强的目的及意义，加深对图像增强的感性认识，巩固所学理论知识。
2、掌握直接灰度变换的图像增强方法。
3、掌握灰度直方图的概念及其计算方法；4、掌握直方图均衡化的计算过程；二、实验原理及知识点1、图像增强是指按特定的需要突出一幅图像中的某些信息，同时，消弱或去除某些不需要的信息的处理方法。
其主要目的是处理后的图像对某些特定的应用比原来的图像更加有效。
图像增强可以在空间域中执行，也可以在变换域中执行。
2、空间域指的是图像平面本身，在空间域内处理图像是直接对图像的像素进行处理。
空间域处理方法分为两种：灰度级变换、空间滤波。
空间域技术直接对像素进行操作，其表达式为g(x,y)=T[f(x,y)]其中f(x,y)为输入图像，g(x,y)为输出图像，T是对图像f进行处理的操作符，定义在点(x,y)的指定邻域内。
定义点(x,y)的空间邻近区域的主要方法是，使用中心位于(x,y)的正方形或长方形区域。
此区域的中心从原点(如左上角)开始逐像素点移动，在移动的同时，该区域会包含不同的邻域。
T应用于每个位置(x,y)，以便在该位置得到输出图像g。
在计算(x,y)处的g值时，只使用该领域的像素。
2、灰度变换T的最简单形式是使用领域大小为1×1，此时,(x,y)处的g值仅由f在该点处的亮度决定，T也变为一个灰度变换函数。
由于灰度变换函数仅取决于亮度的值，而与(x,y)无关，所以亮度函数通常可写做如下所示的简单形式：s=T(r)其中，r表示图像f中相应点(x,y)的亮度，s表示图像g中相应点(x,y)的亮度。
灰度拉伸又叫对比度拉伸是最基本的一种灰度变换，使用简单的分段线性变换函数，可以提高灰度的动态范围，适用于低对比度图像的处理，增强对比度。
3、直方图是多种空间城处理技术的基础。
直方图操作能有效地用于图像增强。
除了提供有用的图像统计资料外，直方图固有的信息在其他图像处理应用中也是非常有用的，如图像压缩与分割。
灰度直方图
第二章
一、灰度直方图的定义
3. 离散图像灰度直方图表明该灰度级上有多少个象素
H ( D) A( D) A( D 1)
H ( D) lim
D
A( D) A( D D) d A( D) D dD
第二章
一、灰度直方图的定义
4. 彩色图像二维直方图
坐标（Dx，Dy）处的值是指在红光图像中具有灰度值Dx，同时在蓝 光图像中具有灰度值Dy的象素对个数。
255
3.通过除以图像的面积来归一化灰度图像可以得到概率密度函数 （PDF）－－与图像本身的象素个数无关 4.对面积函数进行同样的归一化处理可得到图象的累积分布函数 （CDF） 5.图像直方图为不连接区域直方图之和
第二章
三、灰度直方图的应用
1. 数字化参数
观察直方图可
以看出不合适的 数字化
第二章：灰度直方图
一、灰度直方图的定义 二、灰度直方图的性质 三、灰度直方图的应用
第二章
一、灰度直方图的定义
1.灰度直方图是灰度级的函数 2. 连续图像灰度直方图是阀值面积函数的导数的负值
A( D) A( D D) d H ( D) lim A( D) D D dD
第二章
三、灰度直方图的应用
2. 边界阈值选择(物体边界和面积的确定)
简单物体面积的 确定(见灰度直方 图的性质)
选取146为阀值点的二值化
第二章
三、灰度直方图的应用 3. 综合光密度、灰度级平均值
IOD是反映图像质量的一种有用度量
IOD可以通过图像的直方图求得
数字图像处理_期末考试及答案(三篇)
《数字图像处理》试卷答案（2009级）一、名词解释(每题4分，共20分)1．灰度直方图：灰度直方图（histogram）是灰度级的函数，它表示图象中具有每种灰度级的象素的个数，反映图象中每种灰度出现的频率。
它是多种空间域处理技术的基础。
直方图操作能够有效用于图像增强；提供有用的图像统计资料，其在软件中易于计算，适用于商用硬件设备。
灰度直方图性质：1）表征了图像的一维信息。
只反映图像中像素不同灰度值出现的次数（或频数）而未反映像素所在位置。
2）与图像之间的关系是多对一的映射关系。
一幅图像唯一确定出与之对应的直方图，但不同图像可能有相同的直方图。
3）子图直方图之和为整图的直方图。
2．线性移不变系统：一个系统，如果满足线性叠加原理，则称为线性系统，用数学语言可作如下描述：对于，若T[a+b]=aT[]+bT[]=a+b(2.15)则系统T[·]是线性的。
这里，、分别是系统输入，、分别是系统输出。
T[·]表示系统变换，描述了输入输出序列关系，反映出系统特征。
对T[·]加上不同的约束条件，可定义不同的系统。
一个系统，如果系统特征T[·]不受输入序列移位（序列到来的早晚）的影响，则系统称为移不变系统。
由于很多情况下序号对应于时间的顺序，这时也把“移不变”说成是“时不变”。
用数学式表示：对于y(n)= T[x(n)] 若y(n-)=T[x(n-)] (2.16)则系统是移不变的。
既满足线性，又满足移不变条件的系统是线性移不变系统。
这是一种最常用、也最容易理论分析的系统。
这里约定：此后如不加说明，所说的系统均指线性移不变/时不变系统，简称LSI/LTI系统。
3．图像分割：为后续工作有效进行而将图像划分为若干个有意义的区域的技术称为图像分割(Image Segmentation)而目前广为人们所接受的是通过集合所进行的定义：令集合R代表整个图像区域，对R的图像分割可以看做是将R分成N个满足以下条件的非空子集R1，R2，R3，…，RN；（1）在分割结果中，每个区域的像素有着相同的特性（2）在分割结果中，不同子区域具有不同的特性，并且它们没有公共特性（3）分割的所有子区域的并集就是原来的图像（4）各个子集是连通的区域4．数字图像处理：数字图像处理（Digital Image Processing）又称为计算机图像处理，它是指将图像信号转换成数字信号并利用计算机对其进行处理的过程。
灰度直方图
1.灰度直方图灰度直方图（histogram）是灰度级的函数，它表示图象中具有每种灰度级的象素的个数，反映图象中每种灰度出现的频率。
如下图所示，灰度直方图的横坐标是灰度级，纵坐标是该灰度级出现的频率，是图象的最基本的统计特征。
从概率的观点来理解，灰度出现的频率可看作其出现的概率，这样直方图就对应于概率密度函数pdf(probability density function)，而概率分布函数就是直方图的累积和，即概率密度函数的积分，如下图所示：若直接从代表每种灰度的象素数目的直方图来观察，常用如下的表示：drr dP r p dr r p r P r)()(,)()(0==⎰∑⎰⎰⎰===≈======ki ik k k k k rrnnr P n n r p n r n A dAr p dr drr H A r P A drr dA A r H r p drr H A dr r H r A 00000025500)()()(1)(1)(,/)()()()()(,)()(，而概率分布函数，则概率密度的象素数为，灰度为若记象素总数为，时，在离散情况下，取概率密度象素总数一幅图象的总面积，或灰度直方图的计算是很简单的，依据定义，若图象具有L（通常L=256,即8位灰度级）级灰度，则大小为MxN的灰度图象f(x,y)的灰度直方图hist[0…L-1]可用如下计算获得：1. 1.初始化 hist[k]=0 ; k=0,…,L-12. 2.统计 hist[f(x,y)]++ ; x, y =0,…,M-1, 0,…,N-13. 3.标准化 hist[f(x,y)]/=M*N2.直方图均衡化直方图均衡化是通过灰度变换将一幅图象转换为另一幅具有均衡直方图，即在每个灰度级上都具有相同的象素点数的过程。
设灰度变换s=f（r）为斜率有限的非减连续可微函数，它将输入图象A(x,y)转换为输出图象B(x,y)，输入图象的直方图为HA (r)，输出图象的直方图为HB(s)，则它们的关系可由如下过程导出：例如，下图是直方图均衡化后的飞机图片及其直方图，可见其直方图与原图的直方图相比是很均衡的，但必须说明的是，离散情况下不可能作到绝对的一致。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
边界阈值选择
轮廓线提供了一个确立图像中简单物体的边界的 有效的方法。使用轮廓线作为边界的技术被称为 阈值化。 假定一幅图背景是深色的，其中有一个浅色的物 体。物体中的浅色像素产生了直方图上的右峰， 而背景中大量的灰度级产生了直方图上的左峰。 物体边界附近具有两个峰值之间灰度级的像素数 目相对较少，从而产生了两峰之间的谷。选择谷 作为灰度阈值将得到合理的边界。
直方图均衡化
直方图均衡化方法的基本思想是，对在图 像中像素个数多的灰度级进行展宽，而对 像素个数少的灰度级进行缩减。从而达到 清晰图像的目的（增强图像的整体对比 度）。 直方图均衡化处理的中心思想是把原始图 像的灰度直方图从比较集中的某个灰度区 间变成在全部灰度范围内的均匀分布。
直方图匹配（规定化）
对于两峰之间的最低点的灰度级作为阈值来确定边界是最 适合的。直方图是面积函数的导数。在谷底的附近，直方 图的值相对较小，意味着面积函数随阈值灰度级的变化很 缓慢。如果我们选择谷底处的灰度作为阈值，将可以使其 对物体的边界的影响达到最小。如果我们试图测量物体的 面积，选择谷底处阈值将使测量对于阈值灰度变化的敏感 降低到最小。
直方图均衡化校正不具备交互作用特性， 而直方图规定化校正在运用均衡化原理的 基础上，向人们提供了根据给定直方图作 图像增强的手段。 指定希望处理的图像所具有的直方图形状， 这种用于产生处理后又特殊直方图的图像 方法，叫做直方图匹配或直方图规格化处 理。
特点：
直方图均衡化的优点是能自动的增强整体 图像的对比度，但是它的具体增强效果却 不易控制，处理的结果总是得到全局均衡 化的直方图。 直方图规定化可以有选择的增强某个灰度 值范围内的对比度或使图像灰度值的分布 满足特的事图 像中该灰度级的像素个数。即：横坐标表 示灰度级，纵坐标表示图像中该灰度级出 现的频率（对数字图像即指像素的个数）。
直方图的两个性质
图像被缩成直方图后，所有的空间信息都 丢失了。 具有可加性。如果一图像由两个不连接的 区域组成，并且每个区域的直方图已知， 则整幅图像的直方图是该两个区域的直方 图之和。
连续图像中被具有灰度级D的所有轮廓线所包围的面积— —阈值面积函数A(D),则直方图定义为：
综合光密度
综合光密度IOD是反映图像“质量”的一 种度量，可以由直方图计算出来。
其中，K1和K2是所划定的图像区 域的最小和最大灰度值。H（k）代 表灰度级为k时所对应的像素的个 数。该式为灰度级加权的直方图之 和。 连续域分析： 其中a,b为图像区域边界，D （x,y）是图像上不同点的灰度， 相当于“灰度*面积”，面积一定 时，灰度值高则IOD高。 如果图像中的物体被阈值灰度 级 为T的边界勾划出来，则物体边界 内的IOD为： 内部灰度的平均值MGL（Mean Gray Level）等于IOD与面积之比： 对数字图像来说，离散域分析： 是所有像素点的灰度值的总和。