計算 GLCM 紋理

灰度共生矩陣 （Haralick 等，1973）紋理是用於影象分析的強大影象特徵。glcm 包提供了一個易於使用的功能來計算 R 中 RasterLayer 物件的這種中性特徵。

library(glcm)
library(raster)

r <- raster("C:/Program Files/R/R-3.2.3/doc/html/logo.jpg")
plot(r)

在一個方向上計算 GLCM 紋理

rglcm <- glcm(r, 
              window = c(9,9), 
              shift = c(1,1), 
              statistics = c("mean", "variance", "homogeneity", "contrast", 
                             "dissimilarity", "entropy", "second_moment")
              )

plot(rglcm)

計算旋轉不變紋理特徵

紋理特徵也可以在所有 4 個方向（0°，45°，90°和 135°）中計算，然後組合成一個旋轉不變紋理。關鍵是 shift 引數：

rglcm1 <- glcm(r, 
              window = c(9,9), 
              shift=list(c(0,1), c(1,1), c(1,0), c(1,-1)), 
              statistics = c("mean", "variance", "homogeneity", "contrast", 
                             "dissimilarity", "entropy", "second_moment")
              )

plot(rglcm1)