gov_type.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Clear cache\n",
    "arcpy.management.Delete(\"in_memory\")\n",
    "arcpy.management.Delete(\"memory\")\n",
    "# Set the workspace (working directory)\n",
    "arcpy.env.workspace = r\"F:\\LewisK\\PPR_GovType\\PPR_GovType.gdb\"\n",
    "arcpy.env.overwriteOutput = True\n",
    "# Import necessary Libraries\n",
    "from dateutil.relativedelta import relativedelta\n",
    "import datetime\n",
    "from arcgis.gis import GIS\n",
    "import arcpy\n",
    "import os\n",
    "from arcgis.features import FeatureLayer\n",
    "from arcgis.features import FeatureLayerCollection\n",
    "import pandas as pd\n",
    "import math\n",
    "#Coverage by governance type\n",
    "# Identify the month and year of update\n",
    "date = datetime.datetime.today()\n",
    "month = date.strftime(\"%b\")\n",
    "year = date.strftime(\"%Y\")\n",
    "# Define Mollweide projection\n",
    "mollweide_sr = arcpy.SpatialReference(54009)\n",
    "# Define dataset paths in a list of tuples, where each tuple has a dataset name and the dataset path\n",
    "datasets = [\n",
    "     ('wdpa2024', r'O:/PP2024/gov_type_flat.gdb/wdpa_poly_point_public_restricted_basemap_intersect_statusyr2024_Aug2024_GOV_TYPE_union'),\n",
    "     ('wdpa2022', r'O:/PP2024/gov_type_flat.gdb/wdpa_poly_point_public_restricted_basemap_intersect_statusyr2022_Aug2024_GOV_TYPE_union'),\n",
    "    ('wdpa2020', r'O:/PP2024/gov_type_flat.gdb/wdpa_poly_point_public_restricted_basemap_intersect_statusyr2020_Aug2024_GOV_TYPE_union)',\n",
    "    ('wdoecm2024', r'O:/PP2024/gov_type_flat.gdb/wdoecm_poly_point_basemap_intersect_statusyr2024_Aug2024_GOV_TYPE_union'),\n",
    "    ('wdoecm2022', r'O:/PP2024/gov_type_flat.gdb/wdoecm_poly_point_basemap_intersect_statusyr2022_Aug2024_GOV_TYPE_union'),\n",
    "    ('wdoecm2020', r'O:/PP2024/gov_type_flat.gdb/wdoecm_poly_point_basemap_intersect_statusyr2020_Aug2024_GOV_TYPE_union')\n",
    "]\n",
    "print(\"Dataset paths initialized\")\n",
    "# Load basemap and repair geometry\n",
    "basemap = r'O:/PP2024/basemap.gdb/ungeo_eez_gie_national_unsd49regions_v2'\n",
    "arcpy.management.MakeFeatureLayer(basemap, \"basemap2\")\n",
    "arcpy.management.RepairGeometry(\"basemap2\", \"DELETE_NULL\", \"OGC\")\n",
    "print(\"Basemap loaded and repaired\")\n",
    "# Loop through each dataset\n",
    "for dataset_name, dataset_path in datasets:\n",
    "    print(f\"Processing {dataset_name}\")\n",
    "    # Step 1: Load dataset into memory and convert multipart to singlepart\n",
    "    arcpy.management.MakeFeatureLayer(dataset_path, f\"{dataset_name}_mp\")\n",
    "    arcpy.management.MultipartToSinglepart(f\"{dataset_name}_mp\", f\"{dataset_name}_sp\")\n",
    "    print(f\"{dataset_name}: Multipart converted to Singlepart\")\n",
    "    # Step 2: Dice the dataset\n",
    "    arcpy.management.Dice(f\"{dataset_name}_sp\", f\"{dataset_name}_dice\", vertex_limit=\"5000\")\n",
    "    print(f\"{dataset_name}: Singlepart features diced\")\n",
    "    # Step 3: Repair geometry\n",
    "    arcpy.management.RepairGeometry(f\"{dataset_name}_dice\", \"DELETE_NULL\", \"OGC\")\n",
    "    print(f\"{dataset_name}: Geometry repaired\")\n",
    "    # Step 4: Intersect with basemap\n",
    "    arcpy.analysis.PairwiseIntersect(\n",
    "        in_features=[f\"{dataset_name}_dice\", \"basemap2\"],\n",
    "        out_feature_class=f\"{dataset_name}_dice_intersect\",\n",
    "        join_attributes=\"ALL\",\n",
    "        output_type=\"INPUT\"\n",
    "    )\n",
    "    print(f\"{dataset_name}: Intersected with basemap\")\n",
    "    # Step 5: Repair geometry again (after intersect)\n",
    "    arcpy.management.RepairGeometry(f\"{dataset_name}_dice_intersect\", \"DELETE_NULL\", \"OGC\")\n",
    "    print(f\"{dataset_name}: Post-intersect geometry repaired\")\n",
    "    # Step 6: Project to Mollweide\n",
    "    arcpy.management.Project(\n",
    "        in_dataset=f\"{dataset_name}_dice_intersect\",\n",
    "        out_dataset=f\"{dataset_name}_dice_intersect_mollweide\",\n",
    "        out_coor_system=mollweide_sr\n",
    "    )\n",
    "    print(f\"{dataset_name}: Projected to Mollweide\")\n",
    "    # Step 7: Repair geometry one last time (after projection)\n",
    "    arcpy.management.RepairGeometry(f\"{dataset_name}_dice_intersect_mollweide\", \"DELETE_NULL\", \"OGC\")\n",
    "    print(f\"{dataset_name}: Post-projection geometry repaired\")\n",
    "    # Step 8: Calculate geometry attributes\n",
    "    arcpy.management.CalculateGeometryAttributes(\n",
    "        f\"{dataset_name}_dice_intersect_mollweide\",\n",
    "        [[\"PA_AREA\", \"AREA_GEODESIC\"]],\n",
    "        area_unit=\"SQUARE_KILOMETERS\",\n",
    "        coordinate_system=mollweide_sr\n",
    "    )\n",
    "    print(f\"{dataset_name}: Geometry attributes calculated\")\n",
    "    # Step 9: Perform Summary Statistics\n",
    "    # 9.1 Summing PA_AREA by GOV_TYPE\n",
    "    govtype_sum_table = f'F:\\\\LewisK\\PPR_GovType\\\\PPR_GovType.gdb\\\\{dataset_name}_govtype_sum_table'\n",
    "    arcpy.analysis.Statistics(\n",
    "        in_table=f\"{dataset_name}_dice_intersect_mollweide\",\n",
    "        out_table=govtype_sum_table,\n",
    "        statistics_fields=[[\"PA_AREA\", \"SUM\"]],\n",
    "        case_field=[\"GOV_TYPE\", \"TYPE\"]\n",
    "        )\n",
    "    print(f\"{dataset_name}: PA_AREA summed by GOV_TYPE and exported to {govtype_sum_table}\")\n",
    "    # 9.2 Summing PA_AREA by GOV_TYPE and Sub Region Name\n",
    "    govtype_subregion_sum_table = f'F:\\\\LewisK\\PPR_GovType\\\\PPR_GovType.gdb\\\\{dataset_name}_govtype_subregion_sum_table'\n",
    "    arcpy.analysis.Statistics(\n",
    "        in_table=f\"{dataset_name}_dice_intersect_mollweide\",\n",
    "        out_table=govtype_subregion_sum_table,\n",
    "        statistics_fields=[[\"PA_AREA\", \"SUM\"]],\n",
    "        case_field=[\"GOV_TYPE\", \"Sub-region Name\"]\n",
    "    )\n",
    "    print(f\"{dataset_name}: PA_AREA summed by GOV_TYPE and Region Name, and exported to {govtype_subregion_sum_table}\")\n",
    "    # 9.3 Summing PA_AREA by GOV_TYPE and Region Name\n",
    "    govtype_region_sum_table = f'F:\\\\LewisK\\PPR_GovType\\\\PPR_GovType.gdb\\\\{dataset_name}_govtype_region_sum_table'\n",
    "    arcpy.analysis.Statistics(\n",
    "        in_table=f\"{dataset_name}_dice_intersect_mollweide\",\n",
    "        out_table=govtype_region_sum_table,\n",
    "        statistics_fields=[[\"PA_AREA\", \"SUM\"]],\n",
    "        case_field=[\"GOV_TYPE\", \"Region Name\"]\n",
    "        )\n",
    "    print(f\"{dataset_name}: PA_AREA summed by GOV_TYPE and Region Name, and exported to {govtype_region_sum_table}\")\n",
    "print(\"All datasets processed successfully\")"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "ArcGISPro",
   "language": "Python",
   "name": "python3"
  },
  "language_info": {
   "file_extension": ".py",
   "name": "python",
   "version": "3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}