cougar_sightings.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "742735b4-8b0d-47ea-bbb3-ea94a2f65cda",
   "metadata": {},
   "source": [
    "ArcGis Story Map"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1fbb9ed3-fb1f-4eef-9c4b-6d0a8aa7acd6",
   "metadata": {},
   "source": [
    "https://arcg.is/11Sbyi1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "70f13e6c-4bdf-4db1-9022-6fa51038a258",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "41a65f94-5667-4753-a494-feb63ea1e40f",
   "metadata": {},
   "source": [
    "Import Libraries and Excel Table"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "ba896fd3-a000-41fa-a7c9-9bf4746997ae",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import arcpy\n",
    "import os\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "c0be61e1-69b5-4577-8963-6ed7cfa41471",
   "metadata": {},
   "outputs": [],
   "source": [
    "xcel = 'https://raw.githubusercontent.com/ecorey/Eastern_Cougar_Arcpy_/main/sightings_table.xlsx'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "a8abd66c-2190-4f20-b4e0-8f12298f7f7e",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "ea4639aa-272f-4bf2-a1d2-dcc87fda5899",
   "metadata": {},
   "outputs": [],
   "source": [
    "data=pd.read_excel(xcel)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0a3a5c61-aae8-4bd3-bc6e-256cd219a8b4",
   "metadata": {},
   "outputs": [],
   "source": [
    "data.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "8cb700f0-0363-4666-ba4c-6d80bf344b3e",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "bae10b39-66f0-4e71-a368-2c2f4f68d5aa",
   "metadata": {},
   "source": [
    "Query the Table for each of the Sighting Classes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d222478b-a585-4067-bc26-d2f1959cbedd",
   "metadata": {},
   "outputs": [],
   "source": [
    "filtered_data1 = data[data['SightingClass'] == 'Class 1a']\n",
    "filtered_data1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "859d87ba-0f8c-4d20-be8d-5070f99b5c16",
   "metadata": {},
   "outputs": [],
   "source": [
    "filtered_data2 = data[data['SightingClass'] == 'Class  1d']\n",
    "filtered_data2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "dafa92f7-43ed-4f25-9bcf-dea5097d3657",
   "metadata": {},
   "outputs": [],
   "source": [
    "filtered_data3 = data[data['SightingClass'] == 'Class 2a']\n",
    "filtered_data3"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "cfe5ca4e-eea9-4c94-88b2-4ace75078fda",
   "metadata": {},
   "outputs": [],
   "source": [
    "filtered_data4 = data[data['SightingClass'] == 'Class 3b']\n",
    "filtered_data4"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "9dc5b103-971e-4afd-ab55-6914d92c8e26",
   "metadata": {},
   "outputs": [],
   "source": [
    "filtered_data5 = data[data['SightingClass'] == 'Class 3c']\n",
    "filtered_data5"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "a500ad78-af61-4d4e-9a81-26abee3c1219",
   "metadata": {},
   "outputs": [],
   "source": [
    "filtered_data6 = data[data['SightingClass'] == 'Class 4a']\n",
    "filtered_data6"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "5044ace7-8c0f-4e3a-b102-5a81c530c946",
   "metadata": {},
   "outputs": [],
   "source": [
    "filtered_data7 = data[data['SightingClass'] == 'Class 4b']\n",
    "filtered_data7"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "b8d0d650-b4f5-49b4-a45e-8642d9f63c59",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "fe1bd7c5-ab16-450f-885e-eda941745f35",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "783c9c20-79a6-48d2-86d8-8528b8a09d8d",
   "metadata": {},
   "source": [
    "Plot the Frequency of the Sighting Classes with a Bar Graph"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "b21a80d0-d0a9-45c6-af42-941d5a6b1790",
   "metadata": {},
   "outputs": [],
   "source": [
    "value_counts = data['SightingClass'].value_counts()\n",
    "\n",
    "# Plot the bar graph\n",
    "value_counts.plot(kind='bar')\n",
    "plt.xlabel('Attribute Values')\n",
    "plt.ylabel('Frequency')\n",
    "plt.title('Frequency of Attribute Values')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3e4621e8-1661-421c-8911-d6dadb4f8aaf",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Plot the pie chart\n",
    "fig, ax = plt.subplots(figsize=(10, 6))\n",
    "wedges, texts, autotexts = ax.pie(value_counts, startangle=90, autopct='')\n",
    "\n",
    "# Remove the ylabel to make it cleaner\n",
    "ax.set_ylabel('')\n",
    "\n",
    "# Set the title\n",
    "ax.set_title('Frequency of Attribute Values')\n",
    "ax.axis('equal')\n",
    "\n",
    "# Combine the labels with their corresponding percentages\n",
    "total = value_counts.sum()\n",
    "percentage_labels = [f'{label}: {count / total * 100:.1f}%' for label, count in zip(value_counts.index, value_counts)]\n",
    "\n",
    "# Add a legend to the side with the percentage labels\n",
    "ax.legend(wedges, percentage_labels, loc=\"best\", bbox_to_anchor=(1, 0, 0.5, 1))\n",
    "\n",
    "# Show the plot\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4e350de4-d4aa-4c71-b62e-9955c318c133",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "2ddd49d1-c1fe-41f5-995d-581a79db20ff",
   "metadata": {},
   "source": [
    "FEATURE CLASSES FOR SIGHTINGS"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "141f6044-b379-4a6b-bb86-e80dd01f458b",
   "metadata": {},
   "source": [
    "Created FC for Class 1a "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "1ec357db-602f-452f-bb4d-4b4675389943",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Set Excel File and set Sighting Class for Query \n",
    "csv_file = csv_file = 'https://raw.githubusercontent.com/ecorey/Eastern_Cougar_Arcpy_/main/sightingTable_with_degrees.csv'\n",
    "\n",
    "data = pd.read_csv(csv_file)\n",
    "\n",
    "filtered_data = data[data['SightingClass'] == 'Class 1a']\n",
    "\n",
    "# NEED TO SET FOT LOCAL SYSTEM\n",
    "# Set a parent folder for the new file geodatabase\n",
    "parent_folder = r\"C:\\Users\\xyz\\arcGIS_test\\test\"\n",
    "\n",
    "# Set the name of the new file geodatabase\n",
    "gdb_name = \"cougar_sighting_points.gdb\"\n",
    "\n",
    "# Create the new file geodatabase\n",
    "arcpy.CreateFileGDB_management(parent_folder, gdb_name)\n",
    "\n",
    "# NEED TO SET FOT LOCAL SYSTEM\n",
    "# Set the output workspace\n",
    "arcpy.env.workspace = r'C:\\Users\\xyz\\arcGIS_test\\test\\cougar_sighting_points.gdb'\n",
    "arcpy.env.overwriteOutput = True\n",
    "\n",
    "# set the spatial reference system\n",
    "sr = arcpy.SpatialReference(4269)  # WGS 84\n",
    "\n",
    "# NEED TO SET FOT LOCAL SYSTEM\n",
    "# output path\n",
    "outFCPt = r'C:\\Users\\xyz\\arcGIS_test\\test\\cougar_sighting_points.gdb\\fc_pt_1a'\n",
    "geometry_type = \"POINT\"\n",
    "\n",
    "arcpy.CreateFeatureclass_management(os.path.dirname(outFCPt),\n",
    "                                    os.path.basename(outFCPt),\n",
    "                                    \"Point\", spatial_reference=sr)\n",
    "\n",
    "# Fields in the Output Feature Class \n",
    "fields = arcpy.ListFields(outFCPt)\n",
    "field_names = [f.name for f in arcpy.ListFields(outFCPt)]\n",
    "print(field_names)\n",
    "\n",
    "# Creating Point Features from Extracted Coordinates \n",
    "parray = arcpy.Array()  # An array to hold point objects\n",
    "point = arcpy.Point()  # Initialize a point object\n",
    "\n",
    "# Iterate over the filtered_data DataFrame\n",
    "for _, row in filtered_data.iterrows():\n",
    "    lat = row['Y_Cord']  \n",
    "    lon = row['X_Cord']  \n",
    "\n",
    "    point.X = lon\n",
    "    point.Y = lat\n",
    "    parray.append(point)\n",
    "\n",
    "# Creating Point Features from Extracted Coordinates \n",
    "with arcpy.da.InsertCursor(outFCPt, ('SHAPE@XY')) as icur2:  \n",
    "    for p in parray:\n",
    "        icur2.insertRow([p])\n",
    "\n",
    "# Retreive Values From FC for POINT Values        \n",
    "fields = ['SHAPE@XY']\n",
    "\n",
    "with arcpy.da.SearchCursor(outFCPt, fields) as cursor:\n",
    "    for row in cursor:\n",
    "        print(f'{row[0][0]:.2f}')\n",
    "        print(f'{row[0][1]:.2f}')\n",
    "\n",
    "# List Feature Class\n",
    "featureclasses = arcpy.ListFeatureClasses()\n",
    "\n",
    "for fc in featureclasses:\n",
    "    print(fc)\n",
    "\n",
    "# Print number of points\n",
    "num_points = len(parray)\n",
    "print(f\"Total number of points created: {num_points}\")\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5025e432-3d44-4b43-b476-507246c15439",
   "metadata": {},
   "source": [
    "Created FC for Class 1d"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "1cccc6c1-6876-4e74-baed-1ee762508940",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Set Excel File and set Sighting Class for Query \n",
    "csv_file = 'https://raw.githubusercontent.com/ecorey/Eastern_Cougar_Arcpy_/main/sightingTable_with_degrees.csv'\n",
    "data = pd.read_csv(csv_file)\n",
    "\n",
    "filtered_data = data[data['SightingClass'] == 'Class  1d']\n",
    "\n",
    "# NEED TO SET FOT LOCAL SYSTEM\n",
    "# Set a parent folder for the new file geodatabase\n",
    "parent_folder = r\"C:\\Users\\xyz\\arcGIS_test\\test\"\n",
    "\n",
    "# NEED TO SET FOT LOCAL SYSTEM\n",
    "# Set the output workspace\n",
    "arcpy.env.workspace = r'C:\\Users\\xyz\\arcGIS_test\\test\\cougar_sighting_points.gdb'\n",
    "arcpy.env.overwriteOutput = True\n",
    "\n",
    "# set the spatial reference system\n",
    "sr = arcpy.SpatialReference(4269)  # WGS 84\n",
    "\n",
    "# output path\n",
    "outFCPt = r'C:\\Users\\xyz\\arcGIS_test\\test\\cougar_sighting_points.gdb\\fc_pt_1d'\n",
    "geometry_type = \"POINT\"\n",
    "\n",
    "arcpy.CreateFeatureclass_management(os.path.dirname(outFCPt),\n",
    "                                    os.path.basename(outFCPt),\n",
    "                                    \"Point\", spatial_reference=sr)\n",
    "\n",
    "# Fields in the Output Feature Class \n",
    "fields = arcpy.ListFields(outFCPt)\n",
    "field_names = [f.name for f in arcpy.ListFields(outFCPt)]\n",
    "print(field_names)\n",
    "\n",
    "# Creating Point Features from Extracted Coordinates \n",
    "parray = arcpy.Array()  # An array to hold point objects\n",
    "point = arcpy.Point()  # Initialize a point object\n",
    "\n",
    "# Iterate over the filtered_data DataFrame\n",
    "for _, row in filtered_data.iterrows():\n",
    "    lat = row['Y_Cord']  \n",
    "    lon = row['X_Cord']  \n",
    "\n",
    "    point.X = lon\n",
    "    point.Y = lat\n",
    "    parray.append(point)\n",
    "\n",
    "# Creating Point Features from Extracted Coordinates \n",
    "with arcpy.da.InsertCursor(outFCPt, ('SHAPE@XY')) as icur2:  \n",
    "    for p in parray:\n",
    "        icur2.insertRow([p])\n",
    "\n",
    "# Retreive Values From FC for POINT Values        \n",
    "fields = ['SHAPE@XY']\n",
    "\n",
    "with arcpy.da.SearchCursor(outFCPt, fields) as cursor:\n",
    "    for row in cursor:\n",
    "        print(f'{row[0][0]:.2f}')\n",
    "        print(f'{row[0][1]:.2f}')\n",
    "\n",
    "# List Feature Class\n",
    "featureclasses = arcpy.ListFeatureClasses()\n",
    "\n",
    "for fc in featureclasses:\n",
    "    print(fc)\n",
    "\n",
    "# Print number of points\n",
    "num_points = len(parray)\n",
    "print(f\"Total number of points created: {num_points}\")\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e76333a1-73db-497d-8290-5bb15f7e573b",
   "metadata": {},
   "source": [
    "Create FC for Class 2a"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d1f4d851-4d27-48be-9bfa-a5a30192e37f",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Set Excel File and set Sighting Class for Query \n",
    "csv_file = 'https://raw.githubusercontent.com/ecorey/Eastern_Cougar_Arcpy_/main/sightingTable_with_degrees.csv'\n",
    "data = pd.read_csv(csv_file)\n",
    "\n",
    "filtered_data = data[data['SightingClass'] == 'Class 2a']\n",
    "\n",
    "# NEED TO SET FOT LOCAL SYSTEM\n",
    "# Set a parent folder for the new file geodatabase\n",
    "parent_folder = r\"C:\\Users\\xyz\\arcGIS_test\\test\"\n",
    "\n",
    "# NEED TO SET FOT LOCAL SYSTEM\n",
    "# Set the output workspace\n",
    "arcpy.env.workspace = r'C:\\Users\\xyz\\arcGIS_test\\test\\cougar_sighting_points.gdb'\n",
    "arcpy.env.overwriteOutput = True\n",
    "\n",
    "# set the spatial reference system\n",
    "sr = arcpy.SpatialReference(4269)  # WGS 84\n",
    "\n",
    "# output path\n",
    "outFCPt = r'C:\\Users\\xyz\\arcGIS_test\\test\\cougar_sighting_points.gdb\\fc_pt_2a'\n",
    "geometry_type = \"POINT\"\n",
    "\n",
    "arcpy.CreateFeatureclass_management(os.path.dirname(outFCPt),\n",
    "                                    os.path.basename(outFCPt),\n",
    "                                    \"Point\", spatial_reference=sr)\n",
    "\n",
    "# Fields in the Output Feature Class \n",
    "fields = arcpy.ListFields(outFCPt)\n",
    "field_names = [f.name for f in arcpy.ListFields(outFCPt)]\n",
    "print(field_names)\n",
    "\n",
    "# Creating Point Features from Extracted Coordinates \n",
    "parray = arcpy.Array()  # An array to hold point objects\n",
    "point = arcpy.Point()  # Initialize a point object\n",
    "\n",
    "# Iterate over the filtered_data DataFrame\n",
    "for _, row in filtered_data.iterrows():\n",
    "    lat = row['Y_Cord']  \n",
    "    lon = row['X_Cord']  \n",
    "\n",
    "    point.X = lon\n",
    "    point.Y = lat\n",
    "    parray.append(point)\n",
    "\n",
    "# Creating Point Features from Extracted Coordinates \n",
    "with arcpy.da.InsertCursor(outFCPt, ('SHAPE@XY')) as icur2:  \n",
    "    for p in parray:\n",
    "        icur2.insertRow([p])\n",
    "\n",
    "# Retreive Values From FC for POINT Values        \n",
    "fields = ['SHAPE@XY']\n",
    "\n",
    "with arcpy.da.SearchCursor(outFCPt, fields) as cursor:\n",
    "    for row in cursor:\n",
    "        print(f'{row[0][0]:.2f}')\n",
    "        print(f'{row[0][1]:.2f}')\n",
    "\n",
    "# List Feature Class\n",
    "featureclasses = arcpy.ListFeatureClasses()\n",
    "\n",
    "for fc in featureclasses:\n",
    "    print(fc)\n",
    "\n",
    "# Print number of points\n",
    "num_points = len(parray)\n",
    "print(f\"Total number of points created: {num_points}\")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "57b3698e-1145-4680-b886-847a7fa62a14",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "cce3ddbd-30b3-4f20-991c-58fd914d8dbd",
   "metadata": {},
   "source": [
    "Create FC for Class 3b"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "86f2e68a-927d-4b8e-9d7c-5d31059a64b2",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Set Excel File and set Sighting Class for Query \n",
    "csv_file = 'https://raw.githubusercontent.com/ecorey/Eastern_Cougar_Arcpy_/main/sightingTable_with_degrees.csv'\n",
    "data = pd.read_csv(csv_file)\n",
    "\n",
    "filtered_data = data[data['SightingClass'] == 'Class 3b']\n",
    "\n",
    "# NEED TO SET FOT LOCAL SYSTEM\n",
    "# Set a parent folder for the new file geodatabase\n",
    "parent_folder = r\"C:\\Users\\xyz\\arcGIS_test\\test\"\n",
    "\n",
    "# NEED TO SET FOT LOCAL SYSTEM\n",
    "# Set the output workspace\n",
    "arcpy.env.workspace = r'C:\\Users\\xyz\\arcGIS_test\\test\\cougar_sighting_points.gdb'\n",
    "arcpy.env.overwriteOutput = True\n",
    "\n",
    "# set the spatial reference system\n",
    "sr = arcpy.SpatialReference(4269)  # WGS 84\n",
    "\n",
    "# output path\n",
    "outFCPt = r'C:\\Users\\xyz\\arcGIS_test\\test\\cougar_sighting_points.gdb\\fc_pt_3b'\n",
    "geometry_type = \"POINT\"\n",
    "\n",
    "arcpy.CreateFeatureclass_management(os.path.dirname(outFCPt),\n",
    "                                    os.path.basename(outFCPt),\n",
    "                                    \"Point\", spatial_reference=sr)\n",
    "\n",
    "# Fields in the Output Feature Class \n",
    "fields = arcpy.ListFields(outFCPt)\n",
    "field_names = [f.name for f in arcpy.ListFields(outFCPt)]\n",
    "print(field_names)\n",
    "\n",
    "# Creating Point Features from Extracted Coordinates \n",
    "parray = arcpy.Array()  # An array to hold point objects\n",
    "point = arcpy.Point()  # Initialize a point object\n",
    "\n",
    "# Iterate over the filtered_data DataFrame\n",
    "for _, row in filtered_data.iterrows():\n",
    "    lat = row['Y_Cord']  \n",
    "    lon = row['X_Cord']  \n",
    "\n",
    "    point.X = lon\n",
    "    point.Y = lat\n",
    "    parray.append(point)\n",
    "\n",
    "# Creating Point Features from Extracted Coordinates \n",
    "with arcpy.da.InsertCursor(outFCPt, ('SHAPE@XY')) as icur2:  \n",
    "    for p in parray:\n",
    "        icur2.insertRow([p])\n",
    "\n",
    "# Retreive Values From FC for POINT Values        \n",
    "fields = ['SHAPE@XY']\n",
    "\n",
    "with arcpy.da.SearchCursor(outFCPt, fields) as cursor:\n",
    "    for row in cursor:\n",
    "        print(f'{row[0][0]:.2f}')\n",
    "        print(f'{row[0][1]:.2f}')\n",
    "\n",
    "# List Feature Class\n",
    "featureclasses = arcpy.ListFeatureClasses()\n",
    "\n",
    "for fc in featureclasses:\n",
    "    print(fc)\n",
    "\n",
    "# Print number of points\n",
    "num_points = len(parray)\n",
    "print(f\"Total number of points created: {num_points}\")\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "7f4df86b-edca-49a4-ae73-0a0054534c78",
   "metadata": {},
   "source": [
    "Create FC for Class 3c"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "06cf3ee5-d77e-4cbd-bafe-b54906482568",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Set Excel File and set Sighting Class for Query \n",
    "csv_file = 'https://raw.githubusercontent.com/ecorey/Eastern_Cougar_Arcpy_/main/sightingTable_with_degrees.csv'\n",
    "data = pd.read_csv(csv_file)\n",
    "\n",
    "filtered_data = data[data['SightingClass'] == 'Class 3c']\n",
    "\n",
    "# NEED TO SET FOT LOCAL SYSTEM\n",
    "# Set a parent folder for the new file geodatabase\n",
    "parent_folder = r\"C:\\Users\\xyz\\arcGIS_test\\test\"\n",
    "\n",
    "# NEED TO SET FOT LOCAL SYSTEM\n",
    "# Set the output workspace\n",
    "arcpy.env.workspace = r'C:\\Users\\xyz\\arcGIS_test\\test\\cougar_sighting_points.gdb'\n",
    "arcpy.env.overwriteOutput = True\n",
    "\n",
    "# set the spatial reference system\n",
    "sr = arcpy.SpatialReference(4269)  # WGS 84\n",
    "\n",
    "# output path\n",
    "outFCPt = r'C:\\Users\\xyz\\arcGIS_test\\test\\cougar_sighting_points.gdb\\fc_pt_3c'\n",
    "geometry_type = \"POINT\"\n",
    "\n",
    "arcpy.CreateFeatureclass_management(os.path.dirname(outFCPt),\n",
    "                                    os.path.basename(outFCPt),\n",
    "                                    \"Point\", spatial_reference=sr)\n",
    "\n",
    "# Fields in the Output Feature Class \n",
    "fields = arcpy.ListFields(outFCPt)\n",
    "field_names = [f.name for f in arcpy.ListFields(outFCPt)]\n",
    "print(field_names)\n",
    "\n",
    "# Creating Point Features from Extracted Coordinates \n",
    "parray = arcpy.Array()  # An array to hold point objects\n",
    "point = arcpy.Point()  # Initialize a point object\n",
    "\n",
    "# Iterate over the filtered_data DataFrame\n",
    "for _, row in filtered_data.iterrows():\n",
    "    lat = row['Y_Cord']  \n",
    "    lon = row['X_Cord']  \n",
    "\n",
    "    point.X = lon\n",
    "    point.Y = lat\n",
    "    parray.append(point)\n",
    "\n",
    "# Creating Point Features from Extracted Coordinates \n",
    "with arcpy.da.InsertCursor(outFCPt, ('SHAPE@XY')) as icur2:  \n",
    "    for p in parray:\n",
    "        icur2.insertRow([p])\n",
    "\n",
    "# Retreive Values From FC for POINT Values        \n",
    "fields = ['SHAPE@XY']\n",
    "\n",
    "with arcpy.da.SearchCursor(outFCPt, fields) as cursor:\n",
    "    for row in cursor:\n",
    "        print(f'{row[0][0]:.2f}')\n",
    "        print(f'{row[0][1]:.2f}')\n",
    "\n",
    "# List Feature Class\n",
    "featureclasses = arcpy.ListFeatureClasses()\n",
    "\n",
    "for fc in featureclasses:\n",
    "    print(fc)\n",
    "\n",
    "# Print number of points\n",
    "num_points = len(parray)\n",
    "print(f\"Total number of points created: {num_points}\")\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f49e15c4-f892-4c7c-8ffb-cf1a6a21b44d",
   "metadata": {},
   "source": [
    "Create FC for Class 4a"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "c57c1f62-28f3-4ccb-bc12-0cb4b0cdfa4e",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Set Excel File and set Sighting Class for Query \n",
    "csv_file = 'https://raw.githubusercontent.com/ecorey/Eastern_Cougar_Arcpy_/main/sightingTable_with_degrees.csv'\n",
    "data = pd.read_csv(csv_file)\n",
    "\n",
    "filtered_data = data[data['SightingClass'] == 'Class 4a']\n",
    "\n",
    "# NEED TO SET FOT LOCAL SYSTEM\n",
    "# Set a parent folder for the new file geodatabase\n",
    "parent_folder = r\"C:\\Users\\xyz\\arcGIS_test\\test\"\n",
    "\n",
    "# NEED TO SET FOT LOCAL SYSTEM\n",
    "# Set the output workspace\n",
    "arcpy.env.workspace = r'C:\\Users\\xyz\\arcGIS_test\\test\\cougar_sighting_points.gdb'\n",
    "arcpy.env.overwriteOutput = True\n",
    "\n",
    "# set the spatial reference system\n",
    "sr = arcpy.SpatialReference(4269)  # WGS 84\n",
    "\n",
    "# output path\n",
    "outFCPt = r'C:\\Users\\xyz\\arcGIS_test\\test\\cougar_sighting_points.gdb\\fc_pt_4a'\n",
    "geometry_type = \"POINT\"\n",
    "\n",
    "arcpy.CreateFeatureclass_management(os.path.dirname(outFCPt),\n",
    "                                    os.path.basename(outFCPt),\n",
    "                                    \"Point\", spatial_reference=sr)\n",
    "\n",
    "# Fields in the Output Feature Class \n",
    "fields = arcpy.ListFields(outFCPt)\n",
    "field_names = [f.name for f in arcpy.ListFields(outFCPt)]\n",
    "print(field_names)\n",
    "\n",
    "# Creating Point Features from Extracted Coordinates \n",
    "parray = arcpy.Array()  # An array to hold point objects\n",
    "point = arcpy.Point()  # Initialize a point object\n",
    "\n",
    "# Iterate over the filtered_data DataFrame\n",
    "for _, row in filtered_data.iterrows():\n",
    "    lat = row['Y_Cord']  \n",
    "    lon = row['X_Cord']  \n",
    "\n",
    "    point.X = lon\n",
    "    point.Y = lat\n",
    "    parray.append(point)\n",
    "\n",
    "# Creating Point Features from Extracted Coordinates \n",
    "with arcpy.da.InsertCursor(outFCPt, ('SHAPE@XY')) as icur2:  \n",
    "    for p in parray:\n",
    "        icur2.insertRow([p])\n",
    "\n",
    "# Retreive Values From FC for POINT Values        \n",
    "fields = ['SHAPE@XY']\n",
    "\n",
    "with arcpy.da.SearchCursor(outFCPt, fields) as cursor:\n",
    "    for row in cursor:\n",
    "        print(f'{row[0][0]:.2f}')\n",
    "        print(f'{row[0][1]:.2f}')\n",
    "\n",
    "# List Feature Class\n",
    "featureclasses = arcpy.ListFeatureClasses()\n",
    "\n",
    "for fc in featureclasses:\n",
    "    print(fc)\n",
    "\n",
    "# Print number of points\n",
    "num_points = len(parray)\n",
    "print(f\"Total number of points created: {num_points}\")\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "325e7592-db0f-43df-b095-5d680eb482b3",
   "metadata": {},
   "source": [
    "Create FC for Class 4b"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "09fb9dbd-8b19-4707-8bd3-4a4f73ab2c7f",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Set Excel File and set Sighting Class for Query \n",
    "csv_file = 'https://raw.githubusercontent.com/ecorey/Eastern_Cougar_Arcpy_/main/sightingTable_with_degrees.csv'\n",
    "data = pd.read_csv(csv_file)\n",
    "\n",
    "filtered_data = data[data['SightingClass'] == 'Class 4b']\n",
    "\n",
    "# NEED TO SET FOT LOCAL SYSTEM\n",
    "# Set a parent folder for the new file geodatabase\n",
    "parent_folder = r\"C:\\Users\\xyz\\arcGIS_test\\test\"\n",
    "\n",
    "# NEED TO SET FOT LOCAL SYSTEM\n",
    "# Set the output workspace\n",
    "arcpy.env.workspace = r'C:\\Users\\xyz\\arcGIS_test\\test\\cougar_sighting_points.gdb'\n",
    "arcpy.env.overwriteOutput = True\n",
    "\n",
    "# set the spatial reference system\n",
    "sr = arcpy.SpatialReference(4269)  # WGS 84\n",
    "\n",
    "# output path\n",
    "outFCPt = r'C:\\Users\\xyz\\arcGIS_test\\test\\cougar_sighting_points.gdb\\fc_pt_4b'\n",
    "geometry_type = \"POINT\"\n",
    "\n",
    "arcpy.CreateFeatureclass_management(os.path.dirname(outFCPt),\n",
    "                                    os.path.basename(outFCPt),\n",
    "                                    \"Point\", spatial_reference=sr)\n",
    "\n",
    "# Fields in the Output Feature Class \n",
    "fields = arcpy.ListFields(outFCPt)\n",
    "field_names = [f.name for f in arcpy.ListFields(outFCPt)]\n",
    "print(field_names)\n",
    "\n",
    "# Creating Point Features from Extracted Coordinates \n",
    "parray = arcpy.Array()  # An array to hold point objects\n",
    "point = arcpy.Point()  # Initialize a point object\n",
    "\n",
    "# Iterate over the filtered_data DataFrame\n",
    "for _, row in filtered_data.iterrows():\n",
    "    lat = row['Y_Cord']  \n",
    "    lon = row['X_Cord']  \n",
    "\n",
    "    point.X = lon\n",
    "    point.Y = lat\n",
    "    parray.append(point)\n",
    "\n",
    "# Creating Point Features from Extracted Coordinates \n",
    "with arcpy.da.InsertCursor(outFCPt, ('SHAPE@XY')) as icur2:  \n",
    "    for p in parray:\n",
    "        icur2.insertRow([p])\n",
    "\n",
    "# Retreive Values From FC for POINT Values        \n",
    "fields = ['SHAPE@XY']\n",
    "\n",
    "with arcpy.da.SearchCursor(outFCPt, fields) as cursor:\n",
    "    for row in cursor:\n",
    "        print(f'{row[0][0]:.2f}')\n",
    "        print(f'{row[0][1]:.2f}')\n",
    "\n",
    "# List Feature Class\n",
    "featureclasses = arcpy.ListFeatureClasses()\n",
    "\n",
    "for fc in featureclasses:\n",
    "    print(fc)\n",
    "\n",
    "# Print number of points\n",
    "num_points = len(parray)\n",
    "print(f\"Total number of points created: {num_points}\")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "1b9efb28-7bcd-4fdd-b955-d8888fa464c9",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "fda76ceb-1fe5-4fc0-85a6-8ffc59ae71df",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "c6ff5531-9f0e-41ad-9df8-df96a695e32b",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "e434c1d3-f454-429d-821c-7c0f3966d722",
   "metadata": {},
   "source": [
    "CREATE MAPS AND SHOW IMAGE"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 50,
   "id": "7e792dab-b3c4-40c7-b321-6df00d4fc730",
   "metadata": {},
   "outputs": [],
   "source": [
    "from PIL import Image\n",
    "\n",
    "Image.MAX_IMAGE_PIXELS = None \n",
    "image_path = 'https://raw.githubusercontent.com/ecorey/Eastern_Cougar_Arcpy_/main/Layout2.png'\n",
    "image = Image.open(image_path)\n",
    "image.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "ddc269b6-106c-4be3-b713-45c6b82d5f36",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "089f1b0f-9da8-48d3-adc1-8b683e85da0f",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "57ee775c-c894-42a2-8003-ae1d4a60d20a",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0481643f-7020-4f3b-8b82-a7080cca7634",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.16"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}