Python OpenCV 直方图 (五) - 范文中心

Python OpenCV 直方图 (五)

01/26

直方图

Python 调用 calcHist 返回 hist (直方图)

calcHist 函数:

[python] view plain copy

print?

cv2.calcHist([images], channels, mask, histSize, ranges[, hist[, accumulate ]]) #返回hist

参数说明: images  --   图像对像

channels  --  用于计算直方图的通道

Mask  --

histSize  --   表示直方图分成多少份(多少个直方柱)

ranges  --  表示直方图中各个像素的值,[0.0, 256.0]表示直方图能表示像素值从0.0到256的像素

hist   --

accumulate  --  是一个布尔值,表示直方图是否叠加

灰度图计算

[python] view plain copy

print?

#!/usr/bin/env python

# encoding: utf-8

import cv2

import numpy as np

image = cv2.imread("D:/histTest.jpg", 0)

hist = cv2.calcHist([image],

[0], #使用的通道

None, #没有使用mask

[256], #HistSize

[0.0,255.0]) #直方图柱的范围

彩色图像(多通道)直方图

使用 OpenCV 方法

步骤: 1. 读取并分离各通道

2.计算每个通道的直方图

[python] view plain copy

print?

#!/usr/bin/env python

# encoding: utf-8

import cv2

import numpy as np

def calcAndDrawHist(image, color):

hist= cv2.calcHist([image], [0], None, [256], [0.0,255.0])

'''''

minMaxLoc寻找矩阵(一维数组当作向量,用Mat定义) 中最小值和最大值的位置

注意:多通道图像在使用minMaxLoc()函数是不能给出其最大最小值坐标的,因为每个像

素点其实有多个坐标,所以是不会给出的。因此在编程时,这2个位置应该给NULL。

'''

minVal, maxVal, minLoc, maxLoc = cv2.minMaxLoc(hist)

#创建绘制直方图的的图像,由于值全为0 所以它是黑色的

histImg = np.zeros([256,256,3], np.uint8)

hpt = int(0.9* 256) #直方图的范围限定在0-255×0.9之间

for h in range(256):

intensity = int(hist[h]*hpt/maxVal)  # 计算直方图的最大值再乘以一个系数

'''''

绘制线

histImg --   图像

(h,256)  --  线段的第一个端点

(h,256-intensity)  --  线段的第二个端点

color  --  线段的颜色

'''

cv2.line(histImg,(h,256), (h,256-intensity), color)

return histImg;

if __name__ == '__main__':

img = cv2.imread("mini.jpg")

b, g, r = cv2.split(img)     # 使用Opencv 自带的分离函数 split

histImgB = calcAndDrawHist(b, [255, 0, 0])

histImgG = calcAndDrawHist(g, [0, 255, 0])

histImgR = calcAndDrawHist(r, [0, 0, 255])

cv2.imshow("histImgB", histImgB)

cv2.imshow("histImgG", histImgG)

cv2.imshow("histImgR", histImgR)

cv2.imshow("Img", img)

cv2.waitKey(0)

cv2.destroyAllWindows()

在一张图上绘制,这样不用再分离通道,用折线来描绘直方图的边界即可

[python] view plain copy

print?

#!/usr/bin/env python

# encoding: utf-8

import cv2

import numpy as np

img = cv2.imread('mini.jpg')

h = np.zeros((256,256,3)) #创建用于绘制直方图的全0图像

bins = np.arange(256).reshape(256,1) #直方图中各bin的顶点位置

color = [ (255,0,0),(0,255,0),(0,0,255) ] #BGR三种颜色

'''''

对三个通道遍历一次,每次绘制相应通道的直方图的折线

'''

for ch, col in enumerate(color):

#计算对应通道的直方图

originHist = cv2.calcHist([img],[ch],None,[256],[0,256])

'''''

OpenCV的归一化函数。该函数将直方图的范围限定在0-255×0.9之间

'''

cv2.normalize(originHist, originHist,0,255*0.9,cv2.NORM_MINMAX)

'''''

先将生成的原始直方图中的每个元素四舍六入五凑偶取整

(cv2.calcHist函数得到的是float32类型的数组)

注意,这里必须使用np.int32(...)进行转换,

numpy的转换函数可以对数组中的每个元素都进行转换,

而Python的int(...)只能转换一个元素,如果使用int(...),

将导致only length-1 arrays can be converted to Python scalars错误。

'''

hist=np.int32(np.around(originHist))

'''''

将直方图中每个bin的值转成相应的坐标。

如hist[0] =3,...,hist[126] = 178,...,hist[255] = 5;

而bins的值为[[0],[1],[2]...,[255]]

使用np.column_stack将其组合成[0, 3]、[126, 178]、[255, 5]这样的坐标

作为元素组成的数组。

'''

pts = np.column_stack((bins,hist))

'''''

polylines 根据这些点绘制出折线

False  --  指出这个折线不需要闭合

col   --  指定了折线的颜色

'''

cv2.polylines(h,[pts],False,col)

'''''

反转绘制好的直方图,因为绘制时,[0,0]在图像的左上角

'''

h=np.flipud(h)

cv2.imshow('colorhist',h)

cv2.waitKey(0)

使用Numpy 直方图计算

NumPy中histogram函数应用到一个数组返回一对变量:直方图数组和箱式向量。注意:matplotlib也有一个用来建立直方图的函数(叫作hist,正如matlab中一样)与NumPy中的不同。主要的差别是pylab.hist自动绘制直方图,而numpy.histogram仅仅产生数据。

[python] view plain copy

print?

#!/usr/bin/env python

# encoding: utf-8

import cv2

import numpy as np

img = cv2.imread('mini.jpg')

h = np.zeros((300,256,3))

bins = np.arange(257)

bin = bins[0:-1]

color = [ (255,0,0),(0,255,0),(0,0,255) ]

for ch,col in enumerate(color):

item = img[:,:,ch]

N,bins = np.histogram(item,bins)

v=N.max()

N = np.int32(np.around((N*255)/v))

N=N.reshape(256,1)

pts = np.column_stack((bin,N))

cv2.polylines(h,[pts],False,col)

h=np.flipud(h)

cv2.imshow('img',h)

cv2.waitKey(0)


相关内容

  • 大学生创新项目个人总结
    个人总结与思考 时光荏苒,一年的创新项目如期就要结束.能够参加这次项目,我感到非常荣幸,过程中的点点滴滴依然历历在目.因为是第一次参加创新项目,我也从中学习到了很多在日常学习中很难学习到的"知识",也让我感受到实践带来的 ...
  • 锐捷认证掉线问题探究
    锐捷认证掉线问题探究 胡 磊,戴 瑜,李 婧 (空军预警学院训练部,湖北武汉430019) 目前,国内大部分高校校园网均部署锐捷安全管理平台,而锐捷认证掉线问题导致用户无法正常上网则是大摘 要: 多数高校亟需解决的问题.结合锐捷认证系统原理 ...
  • 基于RaspberryPi实现的智能家居
    摘 要 本文以树莓派为载体进行智能家居控制器的设计.在对现在三种市面上的网络电视的比较中引出毕设的主角--树莓派.首先是LINUX 系统常识性的介绍:然后是在ARM/LINUX平台上面开发应用的两种方法:调用系统已有资源.做程序(C 经典法 ...
  • 23. 描述对象的特征
    Python实例属性 类的属性 类的属性包括: 实例属性 类属性 私有属性 特殊属性 实例属性 1.类被实例化后才会具有的属性 2.一般在_init_()方法中创建并初始化 3.直接使用即定义:self. 4.引用方法:self. 5.se ...
  • Adept机械手
    Adept 工业机器人(SCARA机器人.六轴.联轴.线性模块) 在过去的 25 年,Adept 一直是在高速精密制造,包装创新,成为强大的工业机器人发展的行业领先者.Adpet 的工业机器人 创新有 SCARA 机器人,六轴机器人,联轴机 ...
  • 从此乱码是路人
    相信大家在刚学 Qt 的时候一定遇到过 #include #include #include int main(int argc, char *argv[]) { QApplication a(argc, argv); QLabel lb; ...
  • 生物医学图像处理实验指导书 20**年
    实验一 直方图Matlab 运算及C 之间转换 一.实验目的 1. 熟悉利用Matlab 进行图像处理的基本操作,了解图像数据的存储形式及进行图像处理编程的步骤方法. 2. 巩固图像处理编程的步骤格式,理解图像直方图的原理,掌握图像直方图的 ...
  • 动作识别中局部时空特征的运动表示方法研究
    ComputerEngineering andApplications计算机工程与应用 2010.46(34) 7 动作识别巾局部时空特征的运动表示方法研究 雷 LEI 庆1,2,3李绍滋h2Qin91.2.3 LI Shao-zil・2 ...
  • 七年级数学下册统计复习
    一.统计调查 1. 数据处理的过程 (1).数据处理一般包括收集数据.整理数据.描述数据和分析数据等过程. 收集数据的方法:a .民意调查:如投票选举 b.实地调查:如现场进行观察.收集.统计数据 c.媒体调查:报纸.电视.电话.网络等调查 ...
  • MATLAB实现频数直方图
    "hist" is short for "Histogram(直方图.柱状图)". 1.N = hist(Y) bins the elements of Y into 10 equally space ...