Earth Engine 提供多種特殊方法,可用於估算空間紋理。如果圖片為離散值 (非浮點值),您可以使用 image.entropy()
計算鄰域中的熵:
程式碼編輯器 (JavaScript)
// Load a high-resolution NAIP image. var image = ee.Image('USDA/NAIP/DOQQ/m_3712213_sw_10_1_20140613'); // Zoom to San Francisco, display. Map.setCenter(-122.466123, 37.769833, 17); Map.addLayer(image, {max: 255}, 'image'); // Get the NIR band. var nir = image.select('N'); // Define a neighborhood with a kernel. var square = ee.Kernel.square({radius: 4}); // Compute entropy and display. var entropy = nir.entropy(square); Map.addLayer(entropy, {min: 1, max: 5, palette: ['0000CC', 'CC0000']}, 'entropy');
請注意,由於熵值計算需要離散值輸入,因此 NIR 頻帶會在呼叫 entropy()
之前縮放至 8 位元。核中的非零元素會指定鄰域。
另一種評估紋理的做法是使用灰階共現矩陣 (GLCM)。使用上一個範例中的圖像和核,計算以 GLCM 為基礎的對比度,如下所示:
程式碼編輯器 (JavaScript)
// Compute the gray-level co-occurrence matrix (GLCM), get contrast. var glcm = nir.glcmTexture({size: 4}); var contrast = glcm.select('N_contrast'); Map.addLayer(contrast, {min: 0, max: 1500, palette: ['0000CC', 'CC0000']}, 'contrast');
image.glcm()
會輸出許多紋理測量值。如需輸出內容的完整參考資料,請參閱 Haralick 等人 (1973) 和 Conners 等人 (1984)。
您可以使用 image.neighborhoodToBands()
在 Earth Engine 中計算空間關聯的區域指標,例如 Geary 的 C (Anselin 1995)。使用上一個範例中的圖片:
程式碼編輯器 (JavaScript)
// Create a list of weights for a 9x9 kernel. var row = [1, 1, 1, 1, 1, 1, 1, 1, 1]; // The center of the kernel is zero. var centerRow = [1, 1, 1, 1, 0, 1, 1, 1, 1]; // Assemble a list of lists: the 9x9 kernel weights as a 2-D matrix. var rows = [row, row, row, row, centerRow, row, row, row, row]; // Create the kernel from the weights. // Non-zero weights represent the spatial neighborhood. var kernel = ee.Kernel.fixed(9, 9, rows, -4, -4, false); // Convert the neighborhood into multiple bands. var neighs = nir.neighborhoodToBands(kernel); // Compute local Geary's C, a measure of spatial association. var gearys = nir.subtract(neighs).pow(2).reduce(ee.Reducer.sum()) .divide(Math.pow(9, 2)); Map.addLayer(gearys, {min: 20, max: 2500, palette: ['0000CC', 'CC0000']}, "Geary's C");
如需使用鄰近區域標準差來計算圖片紋理的範例,請參閱「圖片鄰近區域的統計資料」一文。