我正在尝试使用 Haralick 描述的 GLCM(能量、均匀性等)为我拥有的一系列 4 波段(R、G、B、NIR)航空照片计算一些纹理测量值。我已经在一个子集上尝试过这个,但最终得到的图像大部分是空白的。我目前的理解是它与灰度和levels
参数有关,但我无法弄清楚。
我的日期非常大(几 GB),所以我试图通过使用模块 RIOS 来提高效率(以 400×400×nbands numpy 数组的形式读取图像,处理数据并写入输出图像)。
我的输入场景可以在这里找到(200 MB)。
我的输出图像看起来像(这可能很难看到,因为黑色像素非常小):
我的代码是:
#Set up input and output filenames
infiles = applier.FilenameAssociations()
infiles.image1 = "infile.tif"
outfiles = applier.FilenameAssociations()
outfiles.outimage = "outfile.tif"
controls = applier.ApplierControls()
controls.progress = cuiprogress.CUIProgressBar()
# I ultimately want to use a window here
# which RIOS easily allows you to set up.
# For a 3x3 the overlap is 1, 5x5 overlap is 2 etc
#controls.setOverlap(4)
def doFilter(info, infiles, outfiles, controls=controls):
grayImg = img_as_ubyte(color.rgb2gray(infiles.image1[3]))
g = greycomatrix(grayImg, distances=[1], angles=[0, np.pi/4, np.pi/2, 3*np.pi/4], symmetric=True, normed=True)
filtered = greycoprops(g, 'energy')
# create 3d image from 2d array
outfiles.outimage = numpy.expand_dims(filtered, axis=0)
applier.apply(doFilter, infiles, outfiles, controls=controls)
显然这里有问题,因为我的输出与我期望的不一样。我猜这与'levels'参数有关。我在这里被指出了一个解释:GLCM 结果中的黑线很好地解释了参数,但我无法改进我的结果。
有人可以向我解释为什么我的结果如图所示以及我该如何补救?