5_biome_extract.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Results saved to O:/PP2024/ecoregion_test/GET/output_tables/biome_cell_count.csv\n"
     ]
    }
   ],
   "source": [
    "### script to calculate PA and OECM coverage of Global Ecosystem Typology biomes\n",
    "## count number of grid cells in each biome and save to a csv (values of 1 are cells in biome). This, along with PA cell count, is used to calculate % coverage\n",
    "\n",
    "import arcpy\n",
    "import os\n",
    "import csv\n",
    "import numpy\n",
    "\n",
    "# Set the workspace and get the list of rasters\n",
    "workspace = r\"O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass\"\n",
    "arcpy.env.workspace = workspace\n",
    "rasters = arcpy.ListRasters()\n",
    "\n",
    "# Output CSV file\n",
    "output_csv = r\"O:/PP2024/ecoregion_test/GET/output_tables/biome_cell_count.csv\"\n",
    "\n",
    "# Initialize a dictionary to hold the results\n",
    "results = {}\n",
    "all_unique_values = set()\n",
    "\n",
    "# Loop through each raster\n",
    "for raster in rasters:\n",
    "    # Get the unique values and their counts using NumPy\n",
    "    raster_path = os.path.join(workspace, raster)\n",
    "    arr = arcpy.RasterToNumPyArray(raster_path)\n",
    "    unique, counts = numpy.unique(arr, return_counts=True)\n",
    "    \n",
    "    # Create a dictionary for the current raster's value counts\n",
    "    value_count_dict = dict(zip(unique, counts))\n",
    "    results[raster] = value_count_dict\n",
    "    \n",
    "    # Update the set of all unique values across all rasters\n",
    "    all_unique_values.update(unique)\n",
    "\n",
    "# Convert the set of all unique values to a sorted list\n",
    "all_unique_values = sorted(all_unique_values)\n",
    "\n",
    "# Write the results to a CSV file\n",
    "with open(output_csv, 'w', newline='') as csvfile:\n",
    "    csvwriter = csv.writer(csvfile)\n",
    "    \n",
    "    # Write the header\n",
    "    header = [\"Raster\"] + [f\"Value_{int(val)}\" for val in all_unique_values]\n",
    "    csvwriter.writerow(header)\n",
    "    \n",
    "    # Write the results for each raster\n",
    "    for raster, value_count_dict in results.items():\n",
    "        row = [raster] + [value_count_dict.get(val, 0) for val in all_unique_values]\n",
    "        csvwriter.writerow(row)\n",
    "\n",
    "print(\"Results saved to\", output_csv)\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Extracted F1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted F1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted F1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted F1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F1_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted F1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F1_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted F1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F1_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted F2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted F2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted F2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted F2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F2_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted F2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F2_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted F2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F2_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted F3_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F3_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted F3_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F3_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted F3_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F3_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted F3_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F3_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted F3_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F3_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted F3_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\F3_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted FM1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\FM1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted FM1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\FM1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted FM1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\FM1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted FM1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\FM1_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted FM1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\FM1_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted FM1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\FM1_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M1_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M1_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M1_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M2_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M2_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M2_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M3_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M3_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M3_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M3_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M3_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M3_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M3_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M3_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M3_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M3_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M3_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M3_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M4_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M4_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M4_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M4_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M4_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M4_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M4_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M4_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M4_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M4_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted M4_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\M4_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MFT1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MFT1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MFT1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MFT1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MFT1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MFT1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MFT1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MFT1_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MFT1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MFT1_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MFT1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MFT1_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MT1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MT1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MT1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MT1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT1_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MT1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT1_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MT1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT1_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Extracted MT2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MT2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MT2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MT2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT2_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MT2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT2_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MT2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT2_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MT3_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT3_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MT3_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT3_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MT3_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT3_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MT3_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT3_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MT3_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT3_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted MT3_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\MT3_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted S1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\S1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted S1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\S1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted S1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\S1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted S1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\S1_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted S1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\S1_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted S1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\S1_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted S2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\S2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted S2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\S2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted S2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\S2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted S2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\S2_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted S2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\S2_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted S2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\S2_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SF1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SF1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SF1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SF1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SF1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SF1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SF1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SF1_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SF1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SF1_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SF1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SF1_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SF2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SF2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SF2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SF2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SF2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SF2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SF2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SF2_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SF2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SF2_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SF2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SF2_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SM1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SM1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SM1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SM1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SM1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SM1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SM1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SM1_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SM1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SM1_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted SM1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\SM1_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T1_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T1_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T1_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T2_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T2_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T2_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T2_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T2_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T2_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T3_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T3_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T3_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T3_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T3_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T3_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T3_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T3_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T3_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T3_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T3_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T3_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Extracted T4_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T4_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T4_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T4_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T4_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T4_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T4_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T4_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T4_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T4_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T4_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T4_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T5_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T5_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T5_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T5_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T5_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T5_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T5_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T5_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T5_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T5_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T5_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T5_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T6_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T6_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T6_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T6_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T6_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T6_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T6_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T6_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T6_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T6_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T6_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T6_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T7_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T7_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T7_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T7_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T7_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T7_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T7_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T7_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T7_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T7_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted T7_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\T7_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted TF1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\TF1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted TF1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\TF1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted TF1_sum.tif by wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\TF1_sum_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted TF1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\TF1_sum_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted TF1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\TF1_sum_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extracted TF1_sum.tif by wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m and saved as O:/PP2024/ecoregion_test/GET/PA_extract\\TF1_sum_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_835m.tif\n",
      "Extraction process completed.\n"
     ]
    }
   ],
   "source": [
    "### Extraxt biome cells within PAs\n",
    "\n",
    "import arcpy\n",
    "import os\n",
    "\n",
    "raster_folder = \"O:/PP2024/ecoregion_test/GET/reclass/moll/sum/sum_reclass\"\n",
    "feature_class_gdb = \"O:/PP2024/ecoregion_test/wdpa_aug24_gte_835m_dice.gdb\"  # Update with the path to your feature class GDB\n",
    "output_workspace = \"O:/PP2024/ecoregion_test/GET/PA_extract\"  # Update with the path to your output folder\n",
    "\n",
    "\n",
    "# Set environment settings\n",
    "arcpy.env.workspace = raster_folder\n",
    "\n",
    "# List all rasters in the raster folder\n",
    "rasters = arcpy.ListRasters()\n",
    "\n",
    "\n",
    "# Change workspace to the geodatabase and list all feature classes\n",
    "arcpy.env.workspace = feature_class_gdb\n",
    "feature_classes = arcpy.ListFeatureClasses()\n",
    "\n",
    "\n",
    "# Loop through each raster and feature class to perform Extract by Mask\n",
    "for raster in rasters:\n",
    "    raster_path = os.path.join(raster_folder, raster)\n",
    "    \n",
    "    for feature_class in feature_classes:\n",
    "        feature_class_path = os.path.join(feature_class_gdb, feature_class)\n",
    "        \n",
    "        # Define the output raster path\n",
    "        output_raster = os.path.join(output_workspace, f\"{os.path.splitext(raster)[0]}_{feature_class}.tif\")\n",
    "        \n",
    "\n",
    "        # Perform Extract by Mask\n",
    "        out_extract_by_mask = arcpy.sa.ExtractByMask(raster_path, feature_class_path)\n",
    "        out_extract_by_mask.save(output_raster)\n",
    "        print(f\"Extracted {raster} by {feature_class} and saved as {output_raster}\")\n",
    "\n",
    "\n",
    "print(\"Extraction process completed.\")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Results saved to O:/PP2024/ecoregion_test/GET/output_tables/biome_cell_count_PAs.csv\n"
     ]
    }
   ],
   "source": [
    "## count number of grid cells in each biome winthin PAs and save to a csv (values of 1 are cells in biome)\n",
    "## may get a 3 column in the output csv, this is the nodata value\n",
    "\n",
    "import arcpy\n",
    "import os\n",
    "import csv\n",
    "import numpy\n",
    "\n",
    "# Set the workspace and get the list of rasters\n",
    "workspace = r\"O:/PP2024/ecoregion_test/GET/PA_extract\"\n",
    "arcpy.env.workspace = workspace\n",
    "rasters = arcpy.ListRasters()\n",
    "\n",
    "# Output CSV file\n",
    "output_csv = r\"O:/PP2024/ecoregion_test/GET/output_tables/biome_cell_count_PAs.csv\"\n",
    "\n",
    "# Initialize a dictionary to hold the results\n",
    "results = {}\n",
    "all_unique_values = set()\n",
    "\n",
    "# Loop through each raster\n",
    "for raster in rasters:\n",
    "    # Get the unique values and their counts using NumPy\n",
    "    raster_path = os.path.join(workspace, raster)\n",
    "    arr = arcpy.RasterToNumPyArray(raster_path)\n",
    "    unique, counts = numpy.unique(arr, return_counts=True)\n",
    "    \n",
    "    # Create a dictionary for the current raster's value counts\n",
    "    value_count_dict = dict(zip(unique, counts))\n",
    "    results[raster] = value_count_dict\n",
    "    \n",
    "    # Update the set of all unique values across all rasters\n",
    "    all_unique_values.update(unique)\n",
    "\n",
    "# Convert the set of all unique values to a sorted list\n",
    "all_unique_values = sorted(all_unique_values)\n",
    "\n",
    "# Write the results to a CSV file\n",
    "with open(output_csv, 'w', newline='') as csvfile:\n",
    "    csvwriter = csv.writer(csvfile)\n",
    "    \n",
    "    # Write the header\n",
    "    header = [\"Raster\"] + [f\"Value_{int(val)}\" for val in all_unique_values]\n",
    "    csvwriter.writerow(header)\n",
    "    \n",
    "    # Write the results for each raster\n",
    "    for raster, value_count_dict in results.items():\n",
    "        row = [raster] + [value_count_dict.get(val, 0) for val in all_unique_values]\n",
    "        csvwriter.writerow(row)\n",
    "\n",
    "print(\"Results saved to\", output_csv)\n",
    "\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "ArcGISPro",
   "language": "Python",
   "name": "python3"
  },
  "language_info": {
   "file_extension": ".py",
   "name": "python",
   "version": "3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}