在使用google earth engine導(dǎo)出圖像時，遇到輸出圖像全黑且數(shù)值為零的問題，通常是由于`export.image.todrive`函數(shù)中`region`參數(shù)設(shè)置不當(dāng)，導(dǎo)致僅導(dǎo)出了單個像素，或因gis軟件默認(rèn)可視化設(shè)置無法正確拉伸單像素值所致。本文將深入探討這一常見問題，并提供詳細(xì)的解決方案和最佳實(shí)踐，確保您能夠成功導(dǎo)出并可視化處理結(jié)果。

Google Earth Engine導(dǎo)出圖像全黑與零值問題解析

在使用Google Earth Engine (GEE) 進(jìn)行遙感數(shù)據(jù)處理后，用戶常會將處理結(jié)果導(dǎo)出到Google Drive或Google Cloud Storage。然而，有時導(dǎo)出的GeoTIFF圖像在ArcGIS、QGIS等GIS軟件中打開時，卻顯示為全黑，并且檢查像素值發(fā)現(xiàn)全部為零。這通常不是數(shù)據(jù)處理錯誤，而是導(dǎo)出參數(shù)配置不當(dāng)或后續(xù)可視化操作不當(dāng)引起的。

問題的核心原因

1. region參數(shù)設(shè)置不當(dāng)： 這是最常見的原因。在Export.image.toDrive函數(shù)中，region參數(shù)定義了要導(dǎo)出圖像的地理范圍。如果將region設(shè)置為一個單一的ee.Geometry.Point對象，GEE只會導(dǎo)出該點(diǎn)所在的單個像素。

   • 后果： 導(dǎo)出的GeoTIFF文件將只包含一個像素點(diǎn)的數(shù)據(jù)。即使這個像素點(diǎn)的值是有效的（非零），在大多數(shù)GIS軟件中，默認(rèn)的拉伸或符號化設(shè)置通常難以正確顯示單個像素。軟件可能會嘗試基于一個像素的“值范圍”進(jìn)行拉伸，導(dǎo)致顯示為全黑。
   • 誤解： 用戶可能會誤認(rèn)為導(dǎo)出的數(shù)據(jù)值就是0，因?yàn)榭梢暬Y(jié)果是黑色的。但實(shí)際上，單個像素可能具有有效的數(shù)值，只是無法被肉眼或默認(rèn)設(shè)置正確識別。

2. 可視化拉伸問題： 即使導(dǎo)出了正確的區(qū)域，如果數(shù)據(jù)的動態(tài)范圍很小或者值集中在某個特定區(qū)間，GIS軟件的默認(rèn)拉伸方式（如Min-Max拉伸）可能無法有效顯示圖像細(xì)節(jié)，導(dǎo)致圖像看起來很暗甚至全黑。對于浮點(diǎn)型數(shù)據(jù)，尤其需要注意。

解決方案與最佳實(shí)踐

解決此問題的關(guān)鍵在于正確配置Export.image.toDrive函數(shù)的region參數(shù)，并理解GIS軟件中的可視化設(shè)置。

1. 修正region參數(shù)

確保region參數(shù)定義了一個有意義的地理區(qū)域，而不是一個單一的點(diǎn)。

• 使用輸入圖像的幾何范圍： 如果您希望導(dǎo)出整個輸入圖像或處理結(jié)果的范圍，可以使用圖像本身的幾何范圍作為導(dǎo)出區(qū)域。

  Export.image.toDrive({
    image: constrained, // 待導(dǎo)出的圖像
    description: 'unmix_result', // 導(dǎo)出任務(wù)描述
    scale: 30, // 導(dǎo)出分辨率，例如Landsat為30米
    region: image.geometry(), // 使用原始輸入圖像的幾何范圍
    maxPixels: 1e9, // 最大像素?cái)?shù)，避免導(dǎo)出任務(wù)因像素過多而失敗
    fileFormat: 'GeoTIFF' // 文件格式
  });

  登錄后復(fù)制

  在上述代碼中，image.geometry()會獲取用于分析的原始影像的邊界框，從而導(dǎo)出整個影像區(qū)域。

• 使用自定義的感興趣區(qū)域 (AOI)： 如果您只想導(dǎo)出特定區(qū)域，應(yīng)定義一個多邊形（ee.Geometry.Polygon）或特征集合（ee.FeatureCollection）作為region。

  

  Imagen – Google Research

  Google Brain team推出的圖像生成模型。

  16

  查看詳情

  // 假設(shè)您已經(jīng)定義了一個名為 'aoi' 的多邊形幾何對象
  var aoi = ee.Geometry.Polygon([
    [[-122.41, 37.77],
     [-122.41, 37.78],
     [-122.42, 37.78],
     [-122.42, 37.77]]
  ]);
  
  Export.image.toDrive({
    image: constrained,
    description: 'unmix_aoi_result',
    scale: 30,
    region: aoi, // 使用自定義AOI
    maxPixels: 1e9,
    fileFormat: 'GeoTIFF'
  });

  登錄后復(fù)制

• 使用當(dāng)前地圖視圖的邊界： 在交互式開發(fā)中，有時您可能想導(dǎo)出當(dāng)前地圖視窗內(nèi)的數(shù)據(jù)。

  Export.image.toDrive({
    image: constrained,
    description: 'unmix_map_view',
    scale: 30,
    region: Map.getBounds(true), // 導(dǎo)出當(dāng)前地圖視圖的邊界
    maxPixels: 1e9,
    fileFormat: 'GeoTIFF'
  });

  登錄后復(fù)制

2. scale參數(shù)的重要性

scale參數(shù)定義了導(dǎo)出圖像的像素分辨率（以米為單位）。

• 過高的分辨率（過小的scale值）： 會導(dǎo)致導(dǎo)出文件過大，可能超出maxPixels限制，或?qū)е聦?dǎo)出時間過長。
• 過低的分辨率（過大的scale值）： 會丟失圖像細(xì)節(jié)。
• 建議： 對于Landsat數(shù)據(jù)，通常使用scale: 30。在測試階段，為了快速導(dǎo)出小文件，可以使用較大的scale值（例如300或3000），但請注意這會降低圖像質(zhì)量。在最終導(dǎo)出時，應(yīng)使用與原始數(shù)據(jù)或所需分析相匹配的合適分辨率。

3. 導(dǎo)出代碼示例（基于原問題優(yōu)化）

針對原問題中region: point導(dǎo)致的問題，修正后的導(dǎo)出代碼如下：

// ... (之前的圖像處理和解混代碼) ...

// Constrained: constraining the result to be non-negative and sum to one.
var constrained = image.unmix([bareMean, vegMean, waterMean, snowMean], true, true);
Map.addLayer(constrained, {}, 'constrained fractions');

// Export output to Google Drive
Export.image.toDrive({
  image: constrained,
  description: 'unmix_result_corrected', // 更改描述以區(qū)分
  scale: 30, // 建議使用原始數(shù)據(jù)分辨率，如30米
  region: image.geometry(), // 關(guān)鍵修正：使用原始圖像的幾何范圍
  maxPixels: 1e9,
  fileFormat: 'GeoTIFF'
});

在上述修正中，我們將region從point改為了image.geometry()，這將確保導(dǎo)出的是整個處理后的圖像區(qū)域，而不是單個像素。同時，將scale調(diào)整為30，以保持Landsat數(shù)據(jù)的原始分辨率。

4. GIS軟件中的可視化設(shè)置

即使導(dǎo)出了正確的圖像，在GIS軟件中打開時仍可能需要手動調(diào)整顯示設(shè)置：

• 檢查數(shù)據(jù)類型和值范圍： 在GIS軟件中查看圖層的屬性，了解其數(shù)據(jù)類型（例如，浮點(diǎn)型）和像素值范圍。
• 調(diào)整拉伸方式：
  • Min-Max拉伸： 適用于值分布較均勻的數(shù)據(jù)。
  • 標(biāo)準(zhǔn)差拉伸： 適用于值集中在均值附近的數(shù)據(jù)。
  • 直方圖均衡化/規(guī)定化： 可以增強(qiáng)圖像的對比度。
  • 自定義拉伸： 根據(jù)數(shù)據(jù)的實(shí)際值范圍手動設(shè)置顯示的最大最小值。
• 應(yīng)用合適的調(diào)色板： 對于單波段圖像，使用合適的漸變色方案（如灰度、光譜色）可以更好地顯示數(shù)據(jù)差異。對于解混結(jié)果，每個波段代表一個組分分?jǐn)?shù)，通常值在0到1之間，可以為每個組分波段單獨(dú)設(shè)置符號化，或?qū)⑺鼈兘M合成RGB圖像（例如，裸地、植被、水體分別映射到R、G、B通道）。

總結(jié)

當(dāng)Google Earth Engine導(dǎo)出的圖像顯示為全黑或零值時，首要排查的是Export.image.toDrive函數(shù)中region參數(shù)的設(shè)置。確保region定義了一個合適的地理范圍（例如，image.geometry()或自定義的AOI），而不是一個單一的像素點(diǎn)。同時，注意scale參數(shù)的設(shè)置，并在GIS軟件中通過調(diào)整可視化拉伸方式和調(diào)色板來正確顯示導(dǎo)出的數(shù)據(jù)。通過這些步驟，您將能夠有效地導(dǎo)出并可視化GEE處理結(jié)果。

以上就是Google Earth Engine導(dǎo)出圖像全黑與零值問題的解析與解決方案的詳細(xì)內(nèi)容，更多請關(guān)注php中文網(wǎng)其它相關(guān)文章！