Gridded GEDI Vegetation Structure Metrics and Biomass Density, 1KM pixel size

LARSE/GEDI/GRIDDEDVEG_002/V1/1KM
データセットの可用性
2019-04-17T00:00:00Z–2023-03-16T00:00:00Z
データセット プロバイダ
Earth Engine スニペット
ee.ImageCollection("LARSE/GEDI/GRIDDEDVEG_002/V1/1KM")
タグ
バイオマス キャノピー 森林 森林バイオマス GEDI LARSE LIDAR NASA 植生

説明

このデータセットは、NASA の地球規模生態系ダイナミクス調査(GEDI)のレベル 2 および 4A のプロダクトから導出された、ほぼ地球規模の分析可能な多重解像度グリッド化された植生構造指標で構成されています。これらの指標は、直径 25 m の LiDAR フットプリントに関連付けられています。このデータセットは、GEDI ライダーのみに基づいて、垂直方向のプロファイル全体を含む、ほぼ地球規模の植生構造を包括的に表現したもので、独立したデータで検証されています。

国際宇宙ステーション(ISS)に搭載された GEDI センサーは、地球表面で軌道に沿って 60 m、軌道に直交して 600 m の間隔で配置された 8 本のレーザービームを使用して、北緯約 52 度から南緯約 52 度の間の地盤の標高と植生構造を測定します。2019 年 4 月 17 日から 2023 年 3 月 16 日までの間に、GEDI はそれぞれ地表の標高と植生の構造の測定に適した 110 億個と 77 億個の高品質な波形を取得しました。

標準的な L2 および L4A のショット指標に加えて、地球システム モデルの炭素循環プロセスや水循環プロセス、森林管理、生物多様性モデリング、生息地の評価などのアプリケーションに特に役立つ可能性のある、いくつかの追加の指標が導出されています。変数には、樹冠の高さ、樹冠被覆率、植物面積指数、葉の高さの多様性、5 m 階層での植物面積体積密度が含まれます。詳しくは、グリッド化された GEDI 植生構造指標とバイオマス密度をご覧ください。

GEDI ショット指標ごとに、平均値、平均値のブートストラップ標準誤差、中央値、標準偏差、四分位範囲、95 パーセンタイル、シャノンの多様性指数、ショット数という 8 つの統計情報が含まれています。GEDI L4B Gridded Aboveground Biomass Density, Version 2.1 に準拠した品質ショット フィルタリング手法が使用されました。対応する GEDI L3 データセットと比較して、このデータセットは複数の空間解像度と複数の時間期間(年単位とミッション全体の期間)で追加のグリッド化された指標を提供します。

このデータセットは、3 つの空間解像度(1 km、6 km、12 km)でラスター グリッドに集計された GEDI ショット指標を提供します。このデータセットでは、1 km のピクセルサイズが使用されます。

バンド

Pixel Size
1000 メートル

帯域

名前 ピクセルサイズ 説明
mean メートル

ピクセル内の GEDI ショット指標値の平均

meanbse メートル

ブートストラップ再サンプリングを使用して計算された平均の標準誤差。100 個のブートストラップ サンプルが作成されました。各サンプルには、ランダムに選択されたショットの 70% が含まれています。標準誤差は、ブートストラップ サンプルの平均値を使用して計算されました。(グリッド セルに 10 個以上の GEDI ショットがある場合にのみ計算されます)。

median メートル

ピクセル内の GEDI ショット指標値の中央値(50 パーセンタイル)。

sd メートル

ピクセル内の GEDI ショット指標値の標準偏差。

iqr メートル

ピクセル内の GEDI ショット指標値の四分位範囲(75 パーセンタイル値から 25 パーセンタイル値を引いた値)。

p95 メートル

ピクセル内の GEDI ショット指標値の 95 パーセンタイル値。

shan メートル

ピクセル内の GEDI ショット指標値のシャノン多様性指数(H)。-1(sum(plog(p))) として計算されます。ここで、p はビンごとの GEDI ショット値の割合です。

countf メートル

ピクセル内の GEDI ショット指標値の数。30 m のサブグリッドを使用して、各 30 m サブグリッド セルで取得された(時間的に)最初の GEDI ショットを選択しました。

利用規約

利用規約

このデータセットはパブリック ドメインにあり、使用と配布に制限はありません。詳しくは、NASA の地球科学データと情報に関するポリシーをご覧ください。

引用

引用:
  • Burns, P.、Hakkenberg, C.R. & Goetz, S.J. GEDI による植生構造のマルチ解像度グリッド マップ。Sci Data 11, 881 (2024). doi:10.1038/s41597-024-03668-4

Earth Engine で探索する

コードエディタ(JavaScript)

var palettes = require('users/gena/packages:palettes');

// GEDI image collections at different spatial resolutions
var ic_1k = ee.ImageCollection('LARSE/GEDI/GRIDDEDVEG_002/V1/1KM')
var ic_6k = ee.ImageCollection('LARSE/GEDI/GRIDDEDVEG_002/V1/6KM')
var ic_12k = ee.ImageCollection('LARSE/GEDI/GRIDDEDVEG_002/V1/12KM')

// slopeshade basemap
var elev = ee.Image('MERIT/DEM/v1_0_3').select('dem')
var slope = ee.Terrain.slope(elev)
Map.setCenter(92.319, 27.129, 8)
Map.addLayer(
    slope, {min: 0, max: 40, palette: ['ffffff', '000000']}, 'Slopeshade')

var opac = 0.7
// View various measurement count metrics from 2019 to 2023
// "va" in the asset name corresponds to the count of high-quality, leaf-on
// vegetation measurements by GEDI
var i_counts_1k_19to23 = ee.Image(
    'LARSE/GEDI/GRIDDEDVEG_002/COUNTS/V1/1KM/gediv002_counts_va_20190417_20230316')
// Number of GEDI laser shots (i.e. footprints) per 1km pixel
Map.addLayer(
    i_counts_1k_19to23.select('shots_count'),
    {min: 0, max: 200, palette: palettes.matplotlib.magma[7]},
    'Shot count per 1km pixel, 2019 to 2023', 0, opac)
// Number of ISS orbits (with valid GEDI shots) per 1km pixel
Map.addLayer(
    i_counts_1k_19to23.select('orbits_uniq'),
    {min: 0, max: 10, palette: palettes.matplotlib.plasma[7]},
    'Unique orbits per 1km pixel, 2019 to 2023', 0, opac)
// The Nearest Neighbor Index (Evans et al. 2023), a proxy for quantifying
// spatial clustering/dispersion of GEDI shots
Map.addLayer(
    i_counts_1k_19to23.select('shots_nni'),
    {min: 0.5, max: 1.5, palette: palettes.colorbrewer.RdBu[7]},
    'Shot nearest neighbor index per 1km pixel, 2019 to 2023', 0, opac)

// View several GEDI vegetation structure metrics at 1km spatial res.
// For GEDI metric descriptions see Table 1 at
// https://daac.ornl.gov/GEDI/guides/GEDI_HighQuality_Shots_Rasters.html
// Relative height of the 98th percentile of returned energy (RH98), a proxy for
// tree canopy height
var i_rh98_1k_19to23 = ee.Image(
    'LARSE/GEDI/GRIDDEDVEG_002/V1/1KM/gediv002_rh-98-a0_vf_20190417_20230316')
// display the median value of GEDI RH98 measurements per 1km pixel, masking out
// values less than 3 consider using a threshold of 10 shots per 1km pixel. More
// shots generally yield more accurate estimates of the aggregation statistics
// (different bands)
var i_rh98_1k_19to23_med = i_rh98_1k_19to23.select('median')
var rh98_pal = palettes.crameri.bamako[10].reverse()
Map.addLayer(
    i_rh98_1k_19to23_med.updateMask(i_rh98_1k_19to23_med.gte(3).and(
        i_rh98_1k_19to23.select('countf').gte(10))),
    {min: 3, max: 40, palette: rh98_pal},
    'Median RH98 per 1km pixel, 2019 to 2023', 1, opac)
// Standard deviation of RH98 per 1km pixel. Captures variability of GEDI
// measurements and vegetation heterogeneity
Map.addLayer(
    i_rh98_1k_19to23.select('sd').updateMask(i_rh98_1k_19to23_med.gte(3).and(
        i_rh98_1k_19to23.select('countf').gte(10))),
    {min: 2, max: 20, palette: palettes.cmocean.Curl[7]},
    'SD of RH98 per 1km pixel, 2019 to 2023', 0, opac)
// Foliage height diversity of the 1m vertical Plant Area Index (PAI) profile
var i_fhd_1k_19to23 = ee.Image(
    'LARSE/GEDI/GRIDDEDVEG_002/V1/1KM/gediv002_fhd-pai-1m-a0_vf_20190417_20230316')
Map.addLayer(
    i_fhd_1k_19to23.select('median').updateMask(i_rh98_1k_19to23_med.gte(3).and(
        i_rh98_1k_19to23.select('countf').gte(10))),
    {min: 1.2, max: 3.2, palette: palettes.matplotlib.viridis[7].reverse()},
    'Median FHD per 1km pixel, 2019 to 2023', 0, opac)
// The height above ground associated with the peak of the vertical Plant Area
// Volume Density (PAVD) profile
var i_pavdmaxh_19to23 = ee.Image(
    'LARSE/GEDI/GRIDDEDVEG_002/V1/1KM/gediv002_pavd-max-h_vf_20190417_20230316')
Map.addLayer(
    i_pavdmaxh_19to23.select('mean').updateMask(i_rh98_1k_19to23_med.gte(3).and(
        i_pavdmaxh_19to23.select('countf').gte(10))),
    {min: 0, max: 25, palette: palettes.cmocean.Haline[7].reverse()},
    'Mean Height of Max. PAVD, 2019 to 2023', 0, opac)

// 1km coverage is not great in the low latitudes, try 6 or 12km for a more
// continuous depiction
var i_rh98_6k_19to23 = ee.Image(
    'LARSE/GEDI/GRIDDEDVEG_002/V1/6KM/gediv002_rh-98-a0_vf_20190417_20230316')
// display the median value of GEDI RH98 measurements per 6km pixel, masking out
// values less than 3 consider using a threshold of 100 shots per 6km pixel.
// More shots generally yield more accurate estimates of the aggregation
// statistics (different bands)
var i_rh98_6k_19to23_med = i_rh98_6k_19to23.select('median')
Map.addLayer(
    i_rh98_6k_19to23_med.updateMask(i_rh98_6k_19to23_med.gte(3).and(
        i_rh98_6k_19to23.select('countf').gte(100))),
    {min: 3, max: 40, palette: rh98_pal},
    'Median RH98 per 6km pixel, 2019 to 2023', 0, opac)
コードエディタで開く