python open cv 图像纹理计算_使用scikit-image和Python实现GLCM纹理特征

我有同样的问题,不同的数据.这是我编写的一个脚本,它使用并行处理和滑动窗口方法：

import gdal, osr

import numpy as np

from scipy.interpolate import RectBivariateSpline

from numpy.lib.stride_tricks import as_strided as ast

import dask.array as da

from joblib import Parallel, delayed, cpu_count

import os

from skimage.feature import greycomatrix, greycoprops

def im_resize(im,Nx,Ny):

'''

resize array by bivariate spline interpolation

'''

ny, nx = np.shape(im)

xx = np.linspace(0,nx,Nx)

yy = np.linspace(0,ny,Ny)

try:

im = da.from_array(im, chunks=1000) #dask implementation

except:

pass

newKernel = RectBivariateSpline(np.r_[:ny],np.r_[:nx],im)

return newKernel(yy,xx)

def p_me(Z, win):

'''

loop to calculate greycoprops

'''

try:

glcm = greycomatrix(Z, [5], [0], 256, symmetric=True, normed=True)

cont = greycoprops(glcm, 'contrast')

diss = greycoprops(glcm, 'dissimilarity')

homo = greycoprops(glcm, 'homogeneity')

eng = greycoprops(glcm, 'energy')

corr = greycoprops(glcm, 'correlation')

ASM = greycoprops(glcm, 'ASM')

return (cont, diss, homo, eng, corr, ASM)

except:

return (0,0,0,0,0,0)

def read_raster(in_raster):

in_raster=in_raster

ds = gdal.Open(in_raster)

data = ds.GetRasterBand(1).ReadAsArray()

data[data<=0] = np.nan

gt = ds.GetGeoTransform()

xres = gt[1]

yres =