|
10 | 10 | "\n", |
11 | 11 | "## Introduction\n", |
12 | 12 | "\n", |
13 | | - "This Colab notebook is designed to process and analyze geospatial data using Google Earth Engine (GEE). It focuses on deriving topographic, climatic, and spectral information for a specified region of interest (ROI). By leveraging GEE's powerful cloud computing capabilities, this workflow enables efficient data handling and export of key environmental variables for further analysis.\n", |
| 13 | + "This Colab notebook is designed to process and analyze geospatial data using Google Earth Engine (GEE). It focuses on deriving topographic, climatic, vegetation structure and spectral information for a specified region of interest (ROI). By leveraging GEE's powerful cloud computing capabilities, this workflow enables efficient data handling and export of key environmental variables for further analysis.\n", |
14 | 14 | "\n", |
15 | 15 | "\n", |
16 | 16 | "### What Does This Notebook Do?:\n", |
|
24 | 24 | "\n", |
25 | 25 | "- Topographic Variables: Elevation, slope, aspect, landforms, and topographic position index (TPI).\n", |
26 | 26 | "- Climate Variables: Total precipitation and average temperature over a defined period.\n", |
| 27 | + "- Structural Variables: Vegetation Height derived from the GPW dataset (30 m resolution). This dataset provides wall-to-wall coverage of all terrestrial ecosystems, used here to represent pre-fire standing fuel from the prior calendar year.\n", |
27 | 28 | "- Spectral Indices: NDVI, NDBI, NDWI, and EVI derived from Sentinel-2 imagery.\n", |
28 | 29 | "- Sentinel-2 Variables: Red, Green, Blue, R1, R2, R3, NIR, R4, SWIR1, and SWIR2 bands with cloud correction included.\n", |
29 | 30 | "- Radar Data: VV and VH polarizations from Sentinel-1 with slope correction applied.\n", |
|
646 | 647 | "metadata": {}, |
647 | 648 | "outputs": [], |
648 | 649 | "source": [ |
649 | | - "# Vegetation Height data from Global Pasture Watch (using season_year - 1 for pre-fire state)\n", |
650 | | - "canopy_collection = ee.ImageCollection(\"projects/global-pasture-watch/assets/gsvh-30m/v1/short-veg-height_m\")\n", |
651 | | - "\n", |
| 650 | + "# GPW provides median vegetation height for all terrestrial ecosystems (not limited to short veg)\n", |
| 651 | + "veg_height_collection = ee.ImageCollection(\"projects/global-pasture-watch/assets/gsvh-30m/v1/short-veg-height_m\")\n", |
652 | 652 | "veg_height = apply_region_mask(\n", |
653 | | - " canopy_collection.filterDate(f'{season_year-1}-01-01', f'{season_year}-01-01') \\\n", |
| 653 | + " veg_height_collection.filterDate(f'{season_year-1}-01-01', f'{season_year}-01-01') \\\n", |
654 | 654 | " .median() \\\n", |
| 655 | + " # The dataset stores values as scaled integers; multiply by 0.1 to convert to meters\n", |
655 | 656 | " .multiply(0.1) \\\n", |
656 | 657 | " .rename('veg_height') \\\n", |
657 | 658 | " .toFloat(),\n", |
|
1092 | 1093 | "with open('/content/processed_data.pkl', 'wb') as f:\n", |
1093 | 1094 | " pickle.dump(data, f)\n", |
1094 | 1095 | "\n", |
1095 | | - "print(\"Data Serilized and saved as processed_data.pkl\")" |
| 1096 | + "print(\"Data Serialized and saved as processed_data.pkl\")" |
1096 | 1097 | ] |
1097 | 1098 | }, |
1098 | 1099 | { |
|
1168 | 1169 | "# Access the variables\n", |
1169 | 1170 | "landform = data['landform']\n", |
1170 | 1171 | "tpi = data['tpi']\n", |
1171 | | - "nasadem_variables = data['elevation_variables']\n", |
| 1172 | + "nasadem_variables = data['nasadem_variables']\n", |
1172 | 1173 | "climate_variables = data['climate_variables']\n", |
1173 | 1174 | "veg_height = data['veg_height'] \n", |
1174 | 1175 | "sentinel2_collections = data['sentinel2_collections']\n", |
|
1260 | 1261 | "batch_2 = [\n", |
1261 | 1262 | " (tpi, 'TPI', 'tpi', region), # 05 minutes aprox\n", |
1262 | 1263 | " (landform, 'Landform', 'landform', region), # 03 minutes aprox\n", |
1263 | | - " (climate_variables, 'Climate Variables', 'climate_variables', region),\n", |
1264 | | - " (veg_height, 'Vegetation Height', 'veg_height', region), # 06 minutes aprox\n", |
| 1264 | + " (climate_variables, 'Climate Variables', 'climate_variables', region), # 06 minutes aprox\n", |
| 1265 | + " (veg_height, 'Vegetation Height', 'veg_height', region), # 06 minutes aprox\n", |
1265 | 1266 | "]\n", |
1266 | 1267 | "\n", |
1267 | 1268 | "batch_3 = [\n", |
|
0 commit comments