灰度共生矩阵 Gray-Level Co-occurrence Matrix(GLCM)原理(一)

2023-10-09 01:59

本文主要是介绍灰度共生矩阵 Gray-Level Co-occurrence Matrix(GLCM)原理(一),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

目录

重要概念

GLCM属性

1. 该矩阵是方形的,即N*N大小,一般为8*8, 16*16,32*32

2.它的行数和列数与图像的量化级别相同。

3. 该矩阵沿对角线对称。

4. offset位移,距离 ,定义为[row_offset, col_offset]

如何产生对称的共生矩阵GLCM图例说明

水平方向共生矩阵

垂直方向共生矩阵

归一化公式

产生对称的共生矩阵GLCM总结

官方说明


重要概念

neighbour pixel: 共生点
ref pixel:参考点,基准点

灰度共生矩阵(GLCM):是一种统计表格,不是一幅图片

GLCM属性

1. 该矩阵是方形的,即N*N大小,一般为8*8, 16*16,32*32

2.它的行数和列数与图像的量化级别相同。

如果图片是8位深度,即灰度范围[0,255]共2^8=256个级别(level),那么行数与列数都是256;但是这样的矩阵计算量太大,一般会缩放至2^3=82^4=16个级别,也有2^5=32个级别的。

3. 该矩阵沿对角线对称。

如果共生点对(0,2)在图片中出现了一次【往东边计算】,那么反过来,往西边计算,共生点对(2,0)在图片中也出现一次。

这样做的目的是:最东边的元素没有右边的共生点可以计算,而最西边的元素没有左边的共生点,让共生矩阵对称后就不会有此问题。

换句话说,就是从东往西计算一次,再从西往东计算一次。【参考offset位移,距离 ,定义为[row_offset, col_offset]】

以下面的原始图为例,以共生点对(0,2)为例,从东往西(0,2)计算出现一次,从西往东(2,0)计算也出现一次,这时不用管谁在左谁在右,只要它们之间的距离是1就行。即水平方向它们出现两次。

问:如何让共生矩阵对称?

答:把产生的统计矩阵 与 它的转置矩阵(transpose)相加,即得到对称的共生矩阵。

4. offset位移,距离 ,定义为[row_offset, col_offset]

此处借用MATLAB的图。

坐标系如下,注意方括号里的值代表的是[row, col]【即(y, x)】上对应的位移,与OpenCV的图片坐标表示是一样的。

如果要产生对称的GLCM,需要计算八个方向【东&西、南&北、东北&西南、西北&东南】;

不对称的GLCM,只需计算四个方向【东、东北、北、西北】

如何产生对称的共生矩阵GLCM图例说明

原始图:

水平方向共生矩阵

Framework matrix 框架矩阵(组合矩阵):

东边(右边)方向的共生点对(0,0) (0,1) (2,2)【(参考点,共生点)】【注意:(2,2)出现了三次】

参考点东边E(1,0)的共生统计结果如下,记为原矩阵 O :

将上面矩阵转置得到下面的转置矩阵T

原矩阵 O 与转置矩阵 T 相加,得到水平方向的共生矩阵,其中共生点对之间的距离为1。

这个共生矩阵有效点对的次数为4+2+1+2+4+1+6+1+1+2=24,而共生点对(2,2)共出现了6次。

那么共生点对(2,2)出现的概率为6/24=1/4 即25%, 共生点对(2,3)的出现概率为1/24=4.2%.

最后一步:归一化

每个格子的值(Value)除以总数(sum)就得到概率矩阵。

上面水平方向的共生矩阵中总数为24,归一化后即得到如下矩阵

垂直方向共生矩阵

Count(south) matrix + Transpose(north) matrix = Symmetrical(vertical) matrix

参考点南边(0,1)的矩阵 + 转置矩阵 = 垂直方向上对称的共生矩阵

概率矩阵:垂直方向上的对称共生矩阵归一化之后的结果如下:

 

归一化公式

产生对称的共生矩阵GLCM总结

产生对称的共生矩阵的步骤:

1. Create a framework matrix taking into account the bit depth【考虑图像位深、灰度级别,建立一个框架矩阵】
2. Decide on the spatial relation between the reference and neighbour pixel【确定参考像素与共生像素的空间关系,如水平方向--东、垂直方向 -- 北、对角线方向--东北、对角线方向--东南等】【注意:只考虑一个方向】
3. Count the occurrences and fill in the framework matrix【计算点对的出现次数,并填入矩阵中,得到原矩阵 O】
4. Add the matrix to its transpose to make it symmetrical【将原矩阵 O 的转置矩阵 T 加到 原矩阵 O 中,产生对称的共生矩阵】
5. Normalize the matrix to conceptually turn it into probabilities.【归一化,转为概率矩阵】

官方说明


1、https://www.mathworks.com/help/images/ref/graycomatrix.html   

【灰度共生矩阵原理,步距、方向原理等说明,及MATLAB代码】
2、https://prism.ucalgary.ca/handle/1880/51900   

【一本入门教程,可以下载为pdf小册子,里面有详细的范例《GLCM Texture: A Tutorial v. 3.0 March 2017》】

3、Haralick, R.M., K. Shanmugan, and I. Dinstein, "Textural Features for Image Classification", IEEE Transactions on Systems, Man, and Cybernetics, Vol. SMC-3, 1973, pp. 610-621.  【方法提出者】
 

这篇关于灰度共生矩阵 Gray-Level Co-occurrence Matrix(GLCM)原理(一)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/169551

相关文章

Java中流式并行操作parallelStream的原理和使用方法

《Java中流式并行操作parallelStream的原理和使用方法》本文详细介绍了Java中的并行流(parallelStream)的原理、正确使用方法以及在实际业务中的应用案例,并指出在使用并行流... 目录Java中流式并行操作parallelStream0. 问题的产生1. 什么是parallelS

Java中Redisson 的原理深度解析

《Java中Redisson的原理深度解析》Redisson是一个高性能的Redis客户端,它通过将Redis数据结构映射为Java对象和分布式对象,实现了在Java应用中方便地使用Redis,本文... 目录前言一、核心设计理念二、核心架构与通信层1. 基于 Netty 的异步非阻塞通信2. 编解码器三、

Java HashMap的底层实现原理深度解析

《JavaHashMap的底层实现原理深度解析》HashMap基于数组+链表+红黑树结构,通过哈希算法和扩容机制优化性能,负载因子与树化阈值平衡效率,是Java开发必备的高效数据结构,本文给大家介绍... 目录一、概述:HashMap的宏观结构二、核心数据结构解析1. 数组(桶数组)2. 链表节点(Node

Redis中Hash从使用过程到原理说明

《Redis中Hash从使用过程到原理说明》RedisHash结构用于存储字段-值对,适合对象数据,支持HSET、HGET等命令,采用ziplist或hashtable编码,通过渐进式rehash优化... 目录一、开篇:Hash就像超市的货架二、Hash的基本使用1. 常用命令示例2. Java操作示例三

Redis中Set结构使用过程与原理说明

《Redis中Set结构使用过程与原理说明》本文解析了RedisSet数据结构,涵盖其基本操作(如添加、查找)、集合运算(交并差)、底层实现(intset与hashtable自动切换机制)、典型应用场... 目录开篇:从购物车到Redis Set一、Redis Set的基本操作1.1 编程常用命令1.2 集

Redis中的有序集合zset从使用到原理分析

《Redis中的有序集合zset从使用到原理分析》Redis有序集合(zset)是字符串与分值的有序映射,通过跳跃表和哈希表结合实现高效有序性管理,适用于排行榜、延迟队列等场景,其时间复杂度低,内存占... 目录开篇:排行榜背后的秘密一、zset的基本使用1.1 常用命令1.2 Java客户端示例二、zse

Redis中的AOF原理及分析

《Redis中的AOF原理及分析》Redis的AOF通过记录所有写操作命令实现持久化,支持always/everysec/no三种同步策略,重写机制优化文件体积,与RDB结合可平衡数据安全与恢复效率... 目录开篇:从日记本到AOF一、AOF的基本执行流程1. 命令执行与记录2. AOF重写机制二、AOF的

java程序远程debug原理与配置全过程

《java程序远程debug原理与配置全过程》文章介绍了Java远程调试的JPDA体系,包含JVMTI监控JVM、JDWP传输调试命令、JDI提供调试接口,通过-Xdebug、-Xrunjdwp参数配... 目录背景组成模块间联系IBM对三个模块的详细介绍编程使用总结背景日常工作中,每个程序员都会遇到bu

Python中isinstance()函数原理解释及详细用法示例

《Python中isinstance()函数原理解释及详细用法示例》isinstance()是Python内置的一个非常有用的函数,用于检查一个对象是否属于指定的类型或类型元组中的某一个类型,它是Py... 目录python中isinstance()函数原理解释及详细用法指南一、isinstance()函数

java 恺撒加密/解密实现原理(附带源码)

《java恺撒加密/解密实现原理(附带源码)》本文介绍Java实现恺撒加密与解密,通过固定位移量对字母进行循环替换,保留大小写及非字母字符,由于其实现简单、易于理解,恺撒加密常被用作学习加密算法的入... 目录Java 恺撒加密/解密实现1. 项目背景与介绍2. 相关知识2.1 恺撒加密算法原理2.2 Ja