5_EFG_to_biome.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F1.1.web.mix_v2.0.tif to O:/PP2024/ecoregion_test/GET\\F1_1_web_mix_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F1.2.web.map_v1.0.tif to O:/PP2024/ecoregion_test/GET\\F1_2_web_map_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F1.3.web.map_v1.0.tif to O:/PP2024/ecoregion_test/GET\\F1_3_web_map_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F1.4.web.map_v1.0.tif to O:/PP2024/ecoregion_test/GET\\F1_4_web_map_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F1.5.web.map_v1.0.tif to O:/PP2024/ecoregion_test/GET\\F1_5_web_map_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F1.6.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\F1_6_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F1.7.web.orig_v2.0.tif to O:/PP2024/ecoregion_test/GET\\F1_7_web_orig_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F2.1.web.alt_v4.0.tif to O:/PP2024/ecoregion_test/GET\\F2_1_web_alt_v4_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F2.10.web.mix_v2.0.tif to O:/PP2024/ecoregion_test/GET\\F2_10_web_mix_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F2.2.web.mix_v2.0.tif to O:/PP2024/ecoregion_test/GET\\F2_2_web_mix_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F2.3.web.mix_v2.0.tif to O:/PP2024/ecoregion_test/GET\\F2_3_web_mix_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F2.4.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\F2_4_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F2.5.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\F2_5_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F2.6.web.mix_v4.0.tif to O:/PP2024/ecoregion_test/GET\\F2_6_web_mix_v4_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F2.7.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\F2_7_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F2.8.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\F2_8_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F2.9.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\F2_9_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F3.1.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\F3_1_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F3.2.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\F3_2_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F3.3.web.alt_v2.0.tif to O:/PP2024/ecoregion_test/GET\\F3_3_web_alt_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F3.4.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\F3_4_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\F3.5.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\F3_5_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\FM1.1.web.map_v1.0.tif to O:/PP2024/ecoregion_test/GET\\FM1_1_web_map_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\FM1.2.web.map_v1.0.tif to O:/PP2024/ecoregion_test/GET\\FM1_2_web_map_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\FM1.3.web.orig_v2.0.tif to O:/PP2024/ecoregion_test/GET\\FM1_3_web_orig_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M1.1.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M1_1_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M1.10.WM.nwx_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M1_10_WM_nwx_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M1.2.web.orig_v2.2.tif to O:/PP2024/ecoregion_test/GET\\M1_2_web_orig_v2_2.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M1.3.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M1_3_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M1.4.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M1_4_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M1.5.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M1_5_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M1.6.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M1_6_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M1.7.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M1_7_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M1.8.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M1_8_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M1.9.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M1_9_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M2.1.web.orig_v2.1.tif to O:/PP2024/ecoregion_test/GET\\M2_1_web_orig_v2_1.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M2.2.web.orig_v2.0.tif to O:/PP2024/ecoregion_test/GET\\M2_2_web_orig_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M2.3.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M2_3_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M2.4.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M2_4_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M2.5.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M2_5_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M3.1.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M3_1_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M3.2.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M3_2_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M3.3.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M3_3_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M3.4.web.orig_v2.0.tif to O:/PP2024/ecoregion_test/GET\\M3_4_web_orig_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M3.5.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M3_5_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M3.6.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M3_6_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M3.7.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M3_7_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M4.1.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M4_1_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\M4.2.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\M4_2_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\MFT1.1.map.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\MFT1_1_map_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\MFT1.2.web.orig_v2.0.tif to O:/PP2024/ecoregion_test/GET\\MFT1_2_web_orig_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\MFT1.3.web.orig_v2.0.tif to O:/PP2024/ecoregion_test/GET\\MFT1_3_web_orig_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\MT1.1.web.map_v1.0.tif to O:/PP2024/ecoregion_test/GET\\MT1_1_web_map_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\MT1.2.IM.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\MT1_2_IM_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\MT1.3.web.map_v1.0.tif to O:/PP2024/ecoregion_test/GET\\MT1_3_web_map_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\MT1.4.web.map_v1.0.tif to O:/PP2024/ecoregion_test/GET\\MT1_4_web_map_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\MT2.1.web.map_v1.0.tif to O:/PP2024/ecoregion_test/GET\\MT2_1_web_map_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\MT2.2.IM.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\MT2_2_IM_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\MT3.1.web.map_v1.0.tif to O:/PP2024/ecoregion_test/GET\\MT3_1_web_map_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\S1.1.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\S1_1_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\S2.1.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\S2_1_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\SF1.1.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\SF1_1_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\SF1.2.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\SF1_2_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\SF2.1.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\SF2_1_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\SF2.2.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\SF2_2_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\SM1.1.web.grid_v3.0.tif to O:/PP2024/ecoregion_test/GET\\SM1_1_web_grid_v3_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\SM1.2.web.grid_v3.0.tif to O:/PP2024/ecoregion_test/GET\\SM1_2_web_grid_v3_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\SM1.3.web.grid_v3.0.tif to O:/PP2024/ecoregion_test/GET\\SM1_3_web_grid_v3_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T1.1.web.mix_v2.0.tif to O:/PP2024/ecoregion_test/GET\\T1_1_web_mix_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T1.2.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T1_2_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T1.3.WM.nwx_v1.1.tif to O:/PP2024/ecoregion_test/GET\\T1_3_WM_nwx_v1_1.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T1.4.web.orig_v2.0.tif to O:/PP2024/ecoregion_test/GET\\T1_4_web_orig_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T2.1.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T2_1_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T2.2.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T2_2_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T2.3.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T2_3_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T2.4.web.orig_v2.0.tif to O:/PP2024/ecoregion_test/GET\\T2_4_web_orig_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T2.5.web.alt_v2.0.tif to O:/PP2024/ecoregion_test/GET\\T2_5_web_alt_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T2.6.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T2_6_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T3.1.web.orig_v2.0.tif to O:/PP2024/ecoregion_test/GET\\T3_1_web_orig_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T3.2.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T3_2_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T3.3.web.alt_v4.0.tif to O:/PP2024/ecoregion_test/GET\\T3_3_web_alt_v4_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T3.4.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T3_4_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T4.1.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T4_1_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T4.2.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T4_2_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T4.3.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T4_3_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T4.4.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T4_4_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T4.5.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T4_5_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T5.1.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T5_1_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T5.2.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T5_2_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T5.3.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T5_3_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T5.4.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T5_4_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T5.5.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T5_5_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T6.1.web.map_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T6_1_web_map_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T6.2.web.alt_v2.0.tif to O:/PP2024/ecoregion_test/GET\\T6_2_web_alt_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T6.3.web.map_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T6_3_web_map_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T6.4.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T6_4_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T6.5.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T6_5_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T7.1.web.alt_v2.0.tif to O:/PP2024/ecoregion_test/GET\\T7_1_web_alt_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T7.2.web.map_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T7_2_web_map_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T7.3.web.alt_v2.0.tif to O:/PP2024/ecoregion_test/GET\\T7_3_web_alt_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T7.4.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T7_4_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\T7.5.web.mix_v1.0.tif to O:/PP2024/ecoregion_test/GET\\T7_5_web_mix_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\TF1.1.web.map_v1.0.tif to O:/PP2024/ecoregion_test/GET\\TF1_1_web_map_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\TF1.2.web.orig_v2.0.tif to O:/PP2024/ecoregion_test/GET\\TF1_2_web_orig_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\TF1.3.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\TF1_3_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\TF1.4.web.map_v1.0.tif to O:/PP2024/ecoregion_test/GET\\TF1_4_web_map_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\TF1.5.web.map_v2.0.tif to O:/PP2024/ecoregion_test/GET\\TF1_5_web_map_v2_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\TF1.6.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\TF1_6_web_orig_v1_0.tif\n",
      "Renamed file: O:/PP2024/ecoregion_test/GET\\TF1.7.web.orig_v1.0.tif to O:/PP2024/ecoregion_test/GET\\TF1_7_web_orig_v1_0.tif\n",
      "All .tif files have been processed.\n"
     ]
    }
   ],
   "source": [
    "#### This script is to aggregate functional group layers up to biome level using data from Global Ecosystem Typology https://global-ecosystems.org/explore\n",
    "## rename files to replace . with _ as the code below fails otherwise\n",
    "\n",
    "import os\n",
    "\n",
    "# Define the path to the target folder\n",
    "folder_path = \"path_to_GET_layers\"\n",
    "\n",
    "# Iterate over all files in the folder\n",
    "for filename in os.listdir(folder_path):\n",
    "    # Check if the file is a .tif file\n",
    "    if filename.lower().endswith('.tif'):\n",
    "        # Generate the new filename by replacing periods with underscores\n",
    "        # Only replace periods in the filename part, not the extension\n",
    "        name, ext = os.path.splitext(filename)\n",
    "        new_name = name.replace('.', '_') + ext\n",
    "        \n",
    "        # Construct the full file paths\n",
    "        original_file_path = os.path.join(folder_path, filename)\n",
    "        new_file_path = os.path.join(folder_path, new_name)\n",
    "        \n",
    "        # Rename the file\n",
    "        os.rename(original_file_path, new_file_path)\n",
    "        \n",
    "        print(f\"Renamed file: {original_file_path} to {new_file_path}\")\n",
    "\n",
    "print(\"All .tif files have been processed.\")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F1_1_web_mix_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F1_2_web_map_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F1_3_web_map_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F1_4_web_map_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F1_5_web_map_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F1_6_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F1_7_web_orig_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F2_10_web_mix_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F2_1_web_alt_v4_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F2_2_web_mix_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F2_3_web_mix_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F2_4_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F2_5_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F2_6_web_mix_v4_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F2_7_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F2_8_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F2_9_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F3_1_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F3_2_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F3_3_web_alt_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F3_4_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/F3_5_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/FM1_1_web_map_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/FM1_2_web_map_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/FM1_3_web_orig_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M1_10_WM_nwx_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M1_1_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M1_2_web_orig_v2_2.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M1_3_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M1_4_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M1_5_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M1_6_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M1_7_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M1_8_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M1_9_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M2_1_web_orig_v2_1.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M2_2_web_orig_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M2_3_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M2_4_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M2_5_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M3_1_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M3_2_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M3_3_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M3_4_web_orig_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M3_5_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M3_6_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M3_7_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M4_1_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/M4_2_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/MFT1_1_map_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/MFT1_2_web_orig_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/MFT1_3_web_orig_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/MT1_1_web_map_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/MT1_2_IM_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/MT1_3_web_map_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/MT1_4_web_map_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/MT2_1_web_map_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/MT2_2_IM_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/MT3_1_web_map_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/S1_1_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/S2_1_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/SF1_1_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/SF1_2_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/SF2_1_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/SF2_2_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/SM1_1_web_grid_v3_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/SM1_2_web_grid_v3_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/SM1_3_web_grid_v3_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T1_1_web_mix_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T1_2_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T1_3_WM_nwx_v1_1.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T1_4_web_orig_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T2_1_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T2_2_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T2_3_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T2_4_web_orig_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T2_5_web_alt_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T2_6_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T3_1_web_orig_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T3_2_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T3_3_web_alt_v4_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T3_4_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T4_1_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T4_2_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T4_3_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T4_4_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T4_5_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T5_1_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T5_2_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T5_3_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T5_4_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T5_5_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T6_1_web_map_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T6_2_web_alt_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T6_3_web_map_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T6_4_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T6_5_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T7_1_web_alt_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T7_2_web_map_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T7_3_web_alt_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T7_4_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/T7_5_web_mix_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/TF1_1_web_map_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/TF1_2_web_orig_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/TF1_3_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/TF1_4_web_map_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/TF1_5_web_map_v2_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/TF1_6_web_orig_v1_0.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/TF1_7_web_orig_v1_0.tif\n",
      "Reclassification completed for all rasters.\n"
     ]
    }
   ],
   "source": [
    "### reclassify GET rasters - keeping major occurances and removing everything else\n",
    "\n",
    "import arcpy\n",
    "from arcpy.sa import *\n",
    "\n",
    "# Set the environment settings\n",
    "arcpy.env.workspace = \"path_to_GET_layers\"\n",
    "arcpy.env.overwriteOutput = True\n",
    "arcpy.env.addOutputsToMap = False\n",
    "\n",
    "# Define the output directory\n",
    "output_directory = \"output_path/reclass\"\n",
    "\n",
    "# List all rasters in the workspace\n",
    "rasters = arcpy.ListRasters()\n",
    "\n",
    "# Check out the Spatial Analyst extension\n",
    "arcpy.CheckOutExtension(\"Spatial\")\n",
    "\n",
    "# Define the reclassification rules\n",
    "remap = RemapValue([[1, 1], [2, 0], [3, 0], [\"NODATA\", 0]])\n",
    "\n",
    "# Process each raster\n",
    "for raster in rasters:\n",
    "    # Perform the reclassification\n",
    "    reclassified_raster = Reclassify(raster, \"Value\", remap, \"NODATA\")\n",
    "    \n",
    "    # Define the output path\n",
    "    output_raster = f\"{output_directory}/{raster}\"\n",
    "    \n",
    "    # Save the reclassified raster\n",
    "    reclassified_raster.save(output_raster)\n",
    "    \n",
    "    print(f\"Reclassified raster saved as {output_raster}\")\n",
    "\n",
    "# Release the Spatial Analyst extension\n",
    "arcpy.CheckInExtension(\"Spatial\")\n",
    "\n",
    "print(\"Reclassification completed for all rasters.\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div class='gpresult'><h2>Messages</h2><div id='messages'>Start Time: 25 July 2024 11:28:03<br>Building Pyramids...<br>Succeeded at 25 July 2024 11:28:54 (Elapsed Time: 50.98 seconds)<br></div><div id = 'default' /></div></div>"
      ],
      "text/plain": [
       "<Result 'O:\\\\PP2024\\\\ecoregion_test\\\\GET\\\\reclass\\\\moll\\\\moll_F1_1_web_mix_v2_0.tif'>"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "## project a single raster so all others can be matched to it\n",
    "\n",
    "# raster to project\n",
    "target = \"O:/PP2024/ecoregion_test/GET/reclass/F1_1_web_mix_v2_0.tif\"\n",
    "\n",
    "# coordinate system to project to\n",
    "sr = arcpy.SpatialReference(\"Mollweide (world)\")\n",
    "\n",
    "# project raster\n",
    "arcpy.management.ProjectRaster(in_raster = target, \n",
    "                               out_raster = \"O:/PP2024/ecoregion_test/GET/reclass/moll/moll_F1_1_web_mix_v2_0.tif\",\n",
    "                               out_coor_system = sr, \n",
    "                               resampling_type = \"NEAREST\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "All rasters have been projected and aligned.\n",
      "snap layers deleted\n"
     ]
    }
   ],
   "source": [
    "### ## project other rasters to match the one above\n",
    "import arcpy\n",
    "import os\n",
    "\n",
    "# Set workspace and output directory\n",
    "workspace = r\"O:/PP2024/ecoregion_test/GET/reclass\"\n",
    "output_folder = r\"O:/PP2024/ecoregion_test/GET/reclass/moll\"\n",
    "arcpy.env.workspace = workspace\n",
    "\n",
    "\n",
    "# List all rasters in the workspace\n",
    "rasters = arcpy.ListRasters()\n",
    "\n",
    "# Define the target projection (Mollweide)\n",
    "mollweide_proj = arcpy.SpatialReference(\"Mollweide (world)\")\n",
    "\n",
    "# Determine the target cell size (use the cell size of the first raster as reference)\n",
    "reference_raster = \"O:/PP2024/ecoregion_test/GET/reclass/moll/moll_F1_1_web_mix_v2_0.tif\"\n",
    "desc = arcpy.Describe(reference_raster)\n",
    "target_cell_size = desc.meanCellWidth  # or you can specify your own cell size\n",
    "\n",
    "# Align rasters by snapping to the lower-left corner of the reference raster\n",
    "snap_raster = os.path.join(output_folder, \"snap_raster.tif\")\n",
    "arcpy.management.ProjectRaster(reference_raster, snap_raster, mollweide_proj, \n",
    "                               cell_size=target_cell_size)\n",
    "arcpy.env.snapRaster = snap_raster\n",
    "\n",
    "# Reproject all rasters\n",
    "for raster in rasters:\n",
    "    # Define the output path\n",
    "    output_raster = os.path.join(output_folder, os.path.basename(raster))\n",
    "    \n",
    "    # Reproject the raster\n",
    "    arcpy.management.ProjectRaster(\n",
    "        in_raster=raster,\n",
    "        out_raster=output_raster,\n",
    "        out_coor_system=mollweide_proj,\n",
    "        resampling_type=\"NEAREST\",  # or BILINEAR/CUBIC if preferred\n",
    "        cell_size=target_cell_size,\n",
    "        geographic_transform=\"\"\n",
    "    )\n",
    "\n",
    "print(\"All rasters have been projected and aligned.\")\n",
    "\n",
    "arcpy.env.snapRaster = None\n",
    "arcpy.management.Delete(\"O:/PP2024/ecoregion_test/GET/reclass/moll/moll_F1_1_web_mix_v2_0.tif\")\n",
    "arcpy.management.Delete(\"O:/PP2024/ecoregion_test/GET/reclass/moll/snap_raster.tif\")\n",
    "print(\"snap layers deleted\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Summed rasters for prefix F1 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/F1_sum.tif\n",
      "Summed rasters for prefix F2 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/F2_sum.tif\n",
      "Summed rasters for prefix F3 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/F3_sum.tif\n",
      "Summed rasters for prefix FM1 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/FM1_sum.tif\n",
      "Summed rasters for prefix M1 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/M1_sum.tif\n",
      "Summed rasters for prefix M2 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/M2_sum.tif\n",
      "Summed rasters for prefix M3 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/M3_sum.tif\n",
      "Summed rasters for prefix M4 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/M4_sum.tif\n",
      "Summed rasters for prefix MFT1 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/MFT1_sum.tif\n",
      "Summed rasters for prefix MT1 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/MT1_sum.tif\n",
      "Summed rasters for prefix MT2 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/MT2_sum.tif\n",
      "Summed rasters for prefix MT3 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/MT3_sum.tif\n",
      "Summed rasters for prefix S1 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/S1_sum.tif\n",
      "Summed rasters for prefix S2 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/S2_sum.tif\n",
      "Summed rasters for prefix SF1 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/SF1_sum.tif\n",
      "Summed rasters for prefix SF2 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/SF2_sum.tif\n",
      "Summed rasters for prefix SM1 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/SM1_sum.tif\n",
      "Summed rasters for prefix T1 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/T1_sum.tif\n",
      "Summed rasters for prefix T2 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/T2_sum.tif\n",
      "Summed rasters for prefix T3 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/T3_sum.tif\n",
      "Summed rasters for prefix T4 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/T4_sum.tif\n",
      "Summed rasters for prefix T5 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/T5_sum.tif\n",
      "Summed rasters for prefix T6 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/T6_sum.tif\n",
      "Summed rasters for prefix T7 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/T7_sum.tif\n",
      "Summed rasters for prefix TF1 saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/TF1_sum.tif\n",
      "All raster summations completed.\n"
     ]
    }
   ],
   "source": [
    "## sum rasters together based on prefix which identifies the biome \n",
    "import arcpy\n",
    "from arcpy.sa import *\n",
    "\n",
    "# Set the environment settings\n",
    "arcpy.env.workspace = r\"O:/PP2024/ecoregion_test/GET/reclass/moll\"\n",
    "arcpy.env.addOutputsToMap = False\n",
    "\n",
    "# List all rasters in the workspace\n",
    "rasters = arcpy.ListRasters()\n",
    "\n",
    "# Create a dictionary to store the lists of rasters by their prefixes\n",
    "raster_groups = {}\n",
    "\n",
    "# Group rasters by their prefixes\n",
    "for raster in rasters:\n",
    "    prefix = raster.split('_')[0]\n",
    "    if prefix not in raster_groups:\n",
    "        raster_groups[prefix] = []\n",
    "    raster_groups[prefix].append(raster)\n",
    "\n",
    "# Sum the rasters in each group\n",
    "for prefix, raster_list in raster_groups.items():\n",
    "    # Initialize the sum with the first raster\n",
    "    raster_sum = Raster(raster_list[0])\n",
    "    \n",
    "    # Add the rest of the rasters in the group\n",
    "    for raster in raster_list[1:]:\n",
    "        raster_sum += Raster(raster)\n",
    "    \n",
    "    # Save the resulting summed raster\n",
    "    output_raster = f\"O:/PP2024/ecoregion_test/GET/reclass/moll/sum/{prefix}_sum.tif\"\n",
    "    raster_sum.save(output_raster)\n",
    "    \n",
    "    print(f\"Summed rasters for prefix {prefix} saved as {output_raster}\")\n",
    "\n",
    "print(\"All raster summations completed.\")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/F1_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/F2_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/F3_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/FM1_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/M1_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/M2_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/M3_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/M4_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/MFT1_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/MT1_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/MT2_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/MT3_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/S1_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/S2_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/SF1_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/SF2_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/SM1_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/T1_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/T2_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/T3_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/T4_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/T5_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/T6_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/T7_sum.tif\n",
      "Reclassified raster saved as O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass/TF1_sum.tif\n",
      "Reclassification completed for all rasters.\n"
     ]
    }
   ],
   "source": [
    "### reclassify GET rasters again to 1 and 0\n",
    "\n",
    "import arcpy\n",
    "from arcpy.sa import *\n",
    "\n",
    "# Set the environment settings\n",
    "arcpy.env.workspace = \"O:/PP2024/ecoregion_test/GET/reclass/moll/sum\"\n",
    "arcpy.env.overwriteOutput = True\n",
    "\n",
    "# Define the output directory\n",
    "output_directory = \"O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass\"\n",
    "\n",
    "# List all rasters in the workspace\n",
    "rasters = arcpy.ListRasters()\n",
    "\n",
    "# Check out the Spatial Analyst extension\n",
    "arcpy.CheckOutExtension(\"Spatial\")\n",
    "\n",
    "# Define the reclassification rules\n",
    "remap = RemapValue([[0, 0], [1, 1], [2, 1], [3, 1], [4, 1],[5, 1],[6, 1],[7, 1],[8, 1], [\"NODATA\", 0]])\n",
    "\n",
    "# Process each raster\n",
    "for raster in rasters:\n",
    "    # Perform the reclassification\n",
    "    reclassified_raster = Reclassify(raster, \"Value\", remap, \"NODATA\")\n",
    "    \n",
    "    # Define the output path\n",
    "    output_raster = f\"{output_directory}/{raster}\"\n",
    "    \n",
    "    # Save the reclassified raster\n",
    "    reclassified_raster.save(output_raster)\n",
    "    \n",
    "    print(f\"Reclassified raster saved as {output_raster}\")\n",
    "\n",
    "# Release the Spatial Analyst extension\n",
    "arcpy.CheckInExtension(\"Spatial\")\n",
    "\n",
    "print(\"Reclassification completed for all rasters.\")"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "ArcGISPro",
   "language": "Python",
   "name": "python3"
  },
  "language_info": {
   "file_extension": ".py",
   "name": "python",
   "version": "3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}