4_NCP_extract_regional.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Zonal Statistics As Table operation completed for all rasters.\n",
      "csvs exported\n",
      "dbfs deleted\n"
     ]
    }
   ],
   "source": [
    "### Script to calculate PA and OECM coverage on NCP layers (produced by Vignesh) by region\n",
    "## count number of grid cells in each NCP layer by region and save to a csv\n",
    "\n",
    "## stop outputs being added to the map\n",
    "arcpy.env.addOutputsToMap = False\n",
    "\n",
    "rasters = [\n",
    "    \"O:/PP2024/NCP/global/potential_cbc_5km_0.05g_1t.tif\",\n",
    "    \"O:/PP2024/NCP/global/realized_cbc_5km_0.05g_1t.tif\"\n",
    "]\n",
    "\n",
    "\n",
    "output_tables = ['O:/PP2024/NCP/global/potential_table.dbf',\n",
    "                'O:/PP2024/NCP/global/realised_table.dbf']\n",
    "\n",
    "zones = 'O:/PP2024/basemap.gdb/EEZv8_WVS_DIS_V3_ALL_final_v7dis_UNSD_Join_moll'\n",
    "zone_field = \"GEOandUNEP\"\n",
    "\n",
    "# Loop through each raster and perform Zonal Statistics As Table\n",
    "for rasters, output_tables in zip(rasters, output_tables):\n",
    "    arcpy.sa.ZonalStatisticsAsTable(\n",
    "        in_zone_data=zones,\n",
    "        zone_field=zone_field,\n",
    "        in_value_raster=rasters,\n",
    "        out_table=output_tables,\n",
    "        statistics_type=\"SUM\"  \n",
    "    )\n",
    "\n",
    "print(\"Zonal Statistics As Table operation completed for all rasters.\")\n",
    "\n",
    "### export as csvs\n",
    "arcpy.conversion.ExportTable('O:/PP2024/NCP/global/potential_table.dbf', \"O:/PP2024/NCP/output_tables/potential_ncp_cell_count_regional.csv\")\n",
    "arcpy.conversion.ExportTable('O:/PP2024/NCP/global/realised_table.dbf', \"O:/PP2024/NCP/output_tables/realized_ncp_cell_count_regional.csv\")\n",
    "\n",
    "print(\"csvs exported\")\n",
    "\n",
    "### delete dbfs\n",
    "arcpy.Delete_management(\"O:/PP2024/NCP/global/potential_table.dbf\")\n",
    "arcpy.Delete_management(\"O:/PP2024/NCP/global/realised_table.dbf\")\n",
    "\n",
    "print(\"dbfs deleted\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "### Extraxt NCP cells within PAs \n",
    "## if this crashes dice the PA layers (OECMs should be ok)\n",
    "\n",
    "import arcpy\n",
    "import os\n",
    "\n",
    "raster_folder = \"O:/PP2024/NCP/global\"\n",
    "feature_class_gdb = \"O:/PP2024/ecoregion_test/wdpa_aug24_gte_5km_dice.gdb\"  # Update with the path to your feature class GDB\n",
    "output_workspace = \"O:/PP2024/NCP/extract/global\"  # 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",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Zonal Statistics As Table operation completed for all rasters.\n"
     ]
    }
   ],
   "source": [
    "## count number of grid cells in each NCP PA layer by region\n",
    "\n",
    "rasters = [\n",
    "    \"O:/PP2024/NCP/extract/global/potential_cbc_5km_0.05g_1t_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_5km.tif\",\n",
    "    \"O:/PP2024/NCP/extract/global/potential_cbc_5km_0.05g_1t_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_5km.tif\",\n",
    "    \"O:/PP2024/NCP/extract/global/potential_cbc_5km_0.05g_1t_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_5km.tif\",\n",
    "    \"O:/PP2024/NCP/extract/global/potential_cbc_5km_0.05g_1t_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_5km.tif\",\n",
    "    \"O:/PP2024/NCP/extract/global/potential_cbc_5km_0.05g_1t_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_5km.tif\",\n",
    "    \"O:/PP2024/NCP/extract/global/potential_cbc_5km_0.05g_1t_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_5km.tif\",\n",
    "    \"O:/PP2024/NCP/extract/global/realized_cbc_5km_0.05g_1t_wdoecm_poly_point_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_5km.tif\",\n",
    "    \"O:/PP2024/NCP/extract/global/realized_cbc_5km_0.05g_1t_wdoecm_poly_point_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_5km.tif\",\n",
    "    \"O:/PP2024/NCP/extract/global/realized_cbc_5km_0.05g_1t_wdoecm_poly_point_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_5km.tif\",\n",
    "    \"O:/PP2024/NCP/extract/global/realized_cbc_5km_0.05g_1t_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_August2024_PA_DEF_moll_gte_5km.tif\",\n",
    "    \"O:/PP2024/NCP/extract/global/realized_cbc_5km_0.05g_1t_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_August2024_PA_DEF_moll_gte_5km.tif\",\n",
    "    \"O:/PP2024/NCP/extract/global/realized_cbc_5km_0.05g_1t_wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_August2024_PA_DEF_moll_gte_5km.tif\",\n",
    "]\n",
    "\n",
    "output_tables = [\"O:/PP2024/NCP/extract/global/potential_oecm_2020.dbf\",\n",
    "                \"O:/PP2024/NCP/extract/global/potential_oecm_2022.dbf\",\n",
    "                \"O:/PP2024/NCP/extract/global/potential_oecm_2024.dbf\",\n",
    "                \"O:/PP2024/NCP/extract/global/potential_wdpa_2020.dbf\",\n",
    "                \"O:/PP2024/NCP/extract/global/potential_wdpa_2022.dbf\",\n",
    "                \"O:/PP2024/NCP/extract/global/potential_wdpa_2024.dbf\",\n",
    "                \"O:/PP2024/NCP/extract/global/realized_oecm_2020.dbf\",\n",
    "                \"O:/PP2024/NCP/extract/global/realized_oecm_2022.dbf\",\n",
    "                \"O:/PP2024/NCP/extract/global/realized_oecm_2024.dbf\",\n",
    "                \"O:/PP2024/NCP/extract/global/realized_wdpa_2020.dbf\",\n",
    "                \"O:/PP2024/NCP/extract/global/realized_wdpa_2022.dbf\",\n",
    "                \"O:/PP2024/NCP/extract/global/realized_wdpa_2024.dbf\"]\n",
    "\n",
    "zones = 'O:/PP2024/basemap.gdb/EEZv8_WVS_DIS_V3_ALL_final_v7dis_UNSD_Join_moll'\n",
    "zone_field = \"GEOandUNEP\"\n",
    "\n",
    "# Loop through each raster and perform Zonal Statistics As Table\n",
    "for rasters, output_tables in zip(rasters, output_tables):\n",
    "    arcpy.sa.ZonalStatisticsAsTable(\n",
    "        in_zone_data=zones,\n",
    "        zone_field=zone_field,\n",
    "        in_value_raster=rasters,\n",
    "        out_table=output_tables,\n",
    "        statistics_type=\"SUM\"  \n",
    "    )\n",
    "\n",
    "print(\"Zonal Statistics As Table operation completed for all rasters.\")\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
}