我们经常提取影像的纹理信息,而提取纹理信息,我们常用灰度共生矩阵,下面就是利用skimage计算图像的GLCM

那么,计算GLCM的是什么包呢:

from skimage.feature import greycomatrix, greycoprops

下面导入的包有些是不用的我也懒得删掉,不影响使用

import math

import numpy as np

import rasterio

from rasterio.mask import mask

import geopandas as gpd

from shapely.geometry import mapping

import pandas as pd

from sklearn.metrics import classification_report

from sklearn.metrics import confusion_matrix

from sklearn.model_selection import train_test_split

import time

import matplotlib.pyplot as plt

import cv2

from PIL import Image

import os

from skimage.feature import greycomatrix, greycoprops

def get_img(s): # s为图像路径

values_temp = []

input = cv2.imread(s, cv2.IMREAD_GRAYSCALE) # 读取图像,灰度模式

# 得到共生矩阵,参数:图像矩阵,距离,方向,灰度级别,是否对称,是否标准化

# [0, np.pi / 4, np.pi / 2, np.pi * 3 / 4] 一共计算了四个方向,你也可以选择一个方向

# 统计得到glcm

glcm = greycomatrix(input, [2, 8, 16], [0, np.pi / 4, np.pi / 2, np.pi * 3 / 4], 256, symmetric=True, normed=True) # , np.pi / 4, np.pi / 2, np.pi * 3 / 4

print(glcm.shape)

# 循环计算表征纹理的参数

for prop in {'contrast', 'dissimilarity','homogeneity', 'energy', 'correlation', 'ASM'}:

temp = greycoprops(glcm, prop)

# temp=np.array(temp).reshape(-1)

values_temp.append(temp)

print(prop, temp)

print('len:',len(temp))

print('')

return (values_temp)

values = []

if __name__ == '__main__':

for n in range(318):

values.append([])

temp_ = get_img(r'C:\Users\Desktop/grey/'+str(n)+'.png')

values[n].append(np.array(temp_).ravel())

# 因为选择了四个方向,所以最终每一个波段有四个值

python灰度共生矩阵_四:利用python的skimage计算灰度共生矩阵相关推荐

  1. python epub 精品_如何利用Python打包HTML页面为epub?

    最近沉迷于将各种博客和官方文档html转化成pdf,结果用手机看还是不太方便,所以想到将html转化成epub格式的电子书,要用os,re,requests,lxm,zipfile,五个库,在这里分享 ...

  2. 怎么用python找论文_如何利用Python绘制学术论文图表方法

    论文中图表形式多样,常用的处理工具有excel.MATLAB以及Python等,excel自处理的方法有两个缺陷: 1.当数据较多时,容易出现excel"翻白眼"的现象: 2.需要 ...

  3. 利用python进行数据分析_资料 | 利用Python进行数据分析

    下载地址:https://www.yanxishe.com/resourceDetail/1443?from=leiphonecolumn_res0518 以下内容节选自豆瓣: 内容简介 · · · ...

  4. python求加速度_如何利用Python 为自然语言处理加速度

    自去年发布 Python 的指代消解包(coreference resolution package)之后,很多用户开始用它来构建许多应用程序,而这些应用与我们最初的对话应用完全不同. 利用 spaC ...

  5. python数据论文_如何利用Python绘制学术论文图表

    论文中图表形式多样,常用的处理工具有excel.MATLAB以及Python等,excel自处理的方法有两个缺陷:1.当数据较多时,容易出现excel"翻白眼"的现象:2.需要使用 ...

  6. 用python做头像_如何利用python制作微信好友头像照片墙?

    这个不难,主要用到itchat和pillow这2个库,其中itchat用于获取微信好友头像照片,pillow用于拼接头像生成一个照片墙,下面我简单介绍一下实现过程,代码量不多,也很好理解,实验环境wi ...

  7. python绘制太阳系_如何利用Python动态模拟太阳系运转

    前言 提到太阳系,大家可能会想到哥白尼和他的日心说,或是捍卫.发展日心说的斗士布鲁诺,他们像一缕光一样照亮了那个时代的夜空,对历史感兴趣的小伙伴可以深入了解一下,这里就不多说了. 太阳以巨大的引力使周 ...

  8. python画四角星_如何利用Python|对多角星的绘制

    本文的文字及图片来源于网络,仅供学习.交流使用,不具有任何商业用途,版权归原作者所有,如有问题请及时联系我们以作处理 以下文章来源于腾讯云,作者:编程与算法之美. 1 前言 在初步学习python的过 ...

  9. ios安装python的步骤_如何利用 Python 爬虫实现给微信群发新闻早报?(详细)

    点击上方"AirPython",选择"加为星标" 第一时间关注 Python 技术干货! 1. 场景 经常有小伙伴在交流群问我,每天的早报新闻是怎么获取的? 其 ...

  10. python刷新页面_小伙利用Python制作浏览器,网友点评这小伙将来要进腾讯

    导语 好久以前看到有人做过(CSDN上?)..今天心血来潮..就想自己试试..仅此而已..用的PyQt版本是5.10...网上关于PyQt5.10的中文教程=None..所以最后做出来的效果一般(毕竟 ...

最新文章

  1. Eclipse配置注释模板
  2. java之filter原理与责任链模式
  3. Django实现一个简单的中间件,不熟悉中间件的爬坑之路
  4. Java实训项目:GUI学生信息管理系统(2019)【下】
  5. 宇宙是否可以了解,宇宙和计算机科学有关系吗?
  6. Java-JDK:JDK百科
  7. 来自大数据的反思:需要你读懂的10个小故事
  8. 什么是scrum敏捷项目管理
  9. 2018年上半年阅读总结(系统架构师)
  10. C#针式打印机自定义纸张连续打印爬坑过程
  11. NFT和游戏化的革命-通过 Play Earn 尝试世界上第一个国际象棋游戏
  12. 四川c语言二级成绩查询,求历年四川省计算机二级考试c语言真题.
  13. 【转载】TCP/IP协议详解
  14. SpringBoot+Vue实现前后端分离旅游资源信息系统
  15. MySQL 综合分析题_MySQL练习题:综合面试题
  16. 彻底删除柚子新闻、柚子壁纸软件
  17. 二本考生计算机考研建议哪个学校,一般二本建议考研学校都有哪些
  18. Vue报错:Root file specified for compilation Vetur(1261)
  19. android 6.0默认壁纸,我伙呆!安卓6.0壁纸竟然是这样得来的
  20. 考研数二第十二讲 复合函数、反函数、隐函数及参数方程所确定的函数的微分法与一阶微分形式的不变性

热门文章

  1. js中设置元素class的三种方法小结
  2. 用正则表达式去除行首行尾空格
  3. vi命令,印在杯子上,还是印在脑子里。
  4. BS-XX-020基于SSM实现停车位租赁系统
  5. hadoop配置文件聚集功能(设置日志可以查看)
  6. linux下初步实现Keepalived+Nginx高可用
  7. 深入理解空间搜索算法 ——数百万数据中的瞬时搜索
  8. java与javascript之间json格式数据互转
  9. SQL Server中EXISTS结构
  10. 如何从开发人员走向架构师