5_feow.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "649499\n",
      "641030\n",
      "628791\n",
      "54799\n",
      "52869\n",
      "40837\n"
     ]
    }
   ],
   "source": [
    "## script to calculate PA and OECM coverage of marine ecoregions, realms, biomes and provences\n",
    "## stop outputs being added to the map\n",
    "arcpy.env.addOutputsToMap = False\n",
    "\n",
    "# read in ecorgions and wdpa\n",
    "## these should be projected to Mollweide and geom repaired\n",
    "eco = 'O:/PP2024/feow.gdb/feow_hydrosheds_moll' ### need to add ID field\n",
    "wdpa24 = 'O:/PP2024/wdpa_Aug2024_v2.gdb/wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_Aug2024_PA_DEF_moll'\n",
    "wdpa22 = 'O:/PP2024/wdpa_Aug2024_v2.gdb/wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_Aug2024_PA_DEF_moll'\n",
    "wdpa20 = 'O:/PP2024/wdpa_Aug2024_v2.gdb/wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_Aug2024_PA_DEF_moll'\n",
    "oecm24 = 'O:/PP2024/wdoecm_Aug2024_v2.gdb/wdoecm_poly_point_basemap_intersect_statusyr2024_Aug2024_PA_DEF_moll'\n",
    "oecm22 = 'O:/PP2024/wdoecm_Aug2024_v2.gdb/wdoecm_poly_point_basemap_intersect_statusyr2022_Aug2024_PA_DEF_moll'\n",
    "oecm20 = 'O:/PP2024/wdoecm_Aug2024_v2.gdb/wdoecm_poly_point_basemap_intersect_statusyr2020_Aug2024_PA_DEF_moll'\n",
    "\n",
    "wdpa_timeserise = [wdpa24, wdpa22, wdpa20, oecm24, oecm22, oecm20]\n",
    "\n",
    "## count features in each file read in to make sure correct files are being used\n",
    "for i in wdpa_timeserise:\n",
    "    count = arcpy.management.GetCount(i)\n",
    "    print(count)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "intersect complete\n"
     ]
    }
   ],
   "source": [
    "## intersect wdpas with ecoregion\n",
    "\n",
    "output_intersects = ['O:/PP2024/chapter5_final/ecoregions_freshwater/ecoregions_freshwater.gdb/wdpa24_feow_intersect',\n",
    "                    'O:/PP2024/chapter5_final/ecoregions_freshwater/ecoregions_freshwater.gdb/wdpa22_feow_intersect',\n",
    "                    'O:/PP2024/chapter5_final/ecoregions_freshwater/ecoregions_freshwater.gdb/wdpa20_feow_intersect',\n",
    "                     'O:/PP2024/chapter5_final/ecoregions_freshwater/ecoregions_freshwater.gdb/oecm24_feow_intersect',\n",
    "                     'O:/PP2024/chapter5_final/ecoregions_freshwater/ecoregions_freshwater.gdb/oecm22_feow_intersect',\n",
    "                     'O:/PP2024/chapter5_final/ecoregions_freshwater/ecoregions_freshwater.gdb/oecm20_feow_intersect'\n",
    "                    ]\n",
    "\n",
    "for wdpa_timeserise, output_intersects in zip(wdpa_timeserise, output_intersects):\n",
    "    input_features = [eco, wdpa_timeserise]\n",
    "    eco_intersect = arcpy.analysis.PairwiseIntersect(input_features, output_intersects)\n",
    "    \n",
    "print(\"intersect complete\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "repair complete\n",
      "repair complete\n",
      "repair complete\n",
      "repair complete\n",
      "repair complete\n",
      "repair complete\n"
     ]
    }
   ],
   "source": [
    "## repair geometry of intersected layers \n",
    "wdpa24_intersect = 'O:/PP2024/chapter5_final/ecoregions_freshwater/ecoregions_freshwater.gdb/wdpa24_feow_intersect'\n",
    "wdpa22_intersect = 'O:/PP2024/chapter5_final/ecoregions_freshwater/ecoregions_freshwater.gdb/wdpa22_feow_intersect'\n",
    "wdpa20_intersect = 'O:/PP2024/chapter5_final/ecoregions_freshwater/ecoregions_freshwater.gdb/wdpa20_feow_intersect'\n",
    "oecm24_intersect = 'O:/PP2024/chapter5_final/ecoregions_freshwater/ecoregions_freshwater.gdb/oecm24_feow_intersect'\n",
    "oecm22_intersect = 'O:/PP2024/chapter5_final/ecoregions_freshwater/ecoregions_freshwater.gdb/oecm22_feow_intersect'\n",
    "oecm20_intersect = 'O:/PP2024/chapter5_final/ecoregions_freshwater/ecoregions_freshwater.gdb/oecm20_feow_intersect'\n",
    "\n",
    "eco_intersects = [wdpa24_intersect, wdpa22_intersect, wdpa20_intersect, oecm24_intersect, oecm22_intersect, oecm20_intersect]\n",
    "\n",
    "for i in eco_intersects:\n",
    "    arcpy.management.RepairGeometry(i, \"DELETE_NULL\", \"OGC\")\n",
    "    print(\"repair complete\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "## calculate area\n",
    "for i in eco_intersects:\n",
    "     arcpy.management.CalculateGeometryAttributes(\n",
    "        in_features = i,\n",
    "        geometry_property = [[\"area_km2\", \"AREA_GEODESIC\"]],\n",
    "        area_unit = \"SQUARE_KILOMETERS\",\n",
    "        coordinate_system = 54009\n",
    "     )\n",
    "\n",
    "print(\"Geometry calculated\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "statistics calculated\n"
     ]
    }
   ],
   "source": [
    "## calculate PA coverage of each ecoregion\n",
    "output_tables = ['O:/PP2024/chapter5_final/ecoregions_freshwater/eco_fresh_output_tables/wdpa24_feow_intersect.csv',\n",
    "                    'O:/PP2024/chapter5_final/ecoregions_freshwater/eco_fresh_output_tables/wdpa22_feow_intersect.csv',\n",
    "                    'O:/PP2024/chapter5_final/ecoregions_freshwater/eco_fresh_output_tables/wdpa20_feow_intersect.csv',\n",
    "                 'O:/PP2024/chapter5_final/ecoregions_freshwater/eco_fresh_output_tables/oecm24_feow_intersect.csv',\n",
    "                    'O:/PP2024/chapter5_final/ecoregions_freshwater/eco_fresh_output_tables/oecm22_feow_intersect.csv',\n",
    "                    'O:/PP2024/chapter5_final/ecoregions_freshwater/eco_fresh_output_tables/oecm20_feow_intersect.csv'\n",
    "                ]\n",
    "\n",
    "for eco_intersects, output_tables in zip(eco_intersects, output_tables):\n",
    "    arcpy.analysis.Statistics(\n",
    "        in_table = eco_intersects, \n",
    "        out_table = output_tables, \n",
    "        statistics_fields = [[\"area_km2\", \"SUM\"]], \n",
    "        case_field = [\"FEOW_ID\"]\n",
    "    )\n",
    "\n",
    "print(\"statistics calculated\")"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "ArcGISPro",
   "language": "Python",
   "name": "python3"
  },
  "language_info": {
   "file_extension": ".py",
   "name": "python",
   "version": "3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}