2016-03-15 175 views
0

我使用R中的glcm包在Landsat波段上運行灰度共生矩陣(GLCM)的圖像紋理度量,我想知道mean'和'mean_ENVI'/'方差'和'variance_ENVI'。圖像紋理GLCM(mean_Envi)

這似乎是一個基本問題,但我找不到一個很好的解釋。

library(glcm) 

Feb2014_B2 <-raster ("Feb2014_band2x.tif") 

Feb2014.B2.textures3x3 <- glcm(Feb2014_B2, window = c(3, 3), shift = c(1, 1), 
statistics = c("mean", "mean_ENVI", "variance", "variance_ENVI", 
"homogeneity", "contrast", "dissimilarity", "entropy", "second_moment", 
"correlation"), na_opt="any", na_val=NA) 
+0

您能否提供一個小例子(數據和代碼)來顯示我們在說什麼? –

+0

RomanLuštrik我剛剛用一個例子更新了我的問題。謝謝 – Nery

+0

'glcm'函數的引用部分引用了4個引用。你有沒有讀過它們,看看它們是否列出了有用的東西? –

回答

1

我是R glcm軟件包的作者。均值和均值_ENVI之間以及方差和方差_ENVI之間的差異與均值如何計算有關。 glcm包中的平均值和方差按照原始的Haralick publication執行。 ENVI中的均值和方差的實現方式不同。 ENVI將平均值計算爲給定窗口大小內像素值的簡單平均值 - 因此ENVI平均值根本不是紋理度量,而只是原始圖像的平滑版本。我在我的包中實現了它,以便可以在R中複製ENVI結果。

我使用的公式與Haralick中的一樣。現在看到它們的唯一方法是使用C++代碼(https://github.com/azvoleff/glcm/blob/master/src/calc_texture.cpp)。我將最終將這些內容加入到R文檔中,但還沒有時間將它們插入到文檔中。

+0

非常感謝。現在兩者之間的區別很明顯。如果這些信息寫在R文檔中,這將非常有用。 – Nery