Aggregates_Number_of_Pixels_MODIS.ipynb

{
  "nbformat": 4,
  "nbformat_minor": 0,
  "metadata": {
    "colab": {
      "provenance": [],
      "authorship_tag": "ABX9TyMebBIqPRBUjf5pETCaeKkD"
    },
    "kernelspec": {
      "name": "python3",
      "display_name": "Python 3"
    },
    "language_info": {
      "name": "python"
    }
  },
  "cells": [
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "oXnkggvzD5HX"
      },
      "outputs": [],
      "source": [
        "\n",
        "\n",
        "# This code was used to aggregate number of pixels coresponding to the IFL_2000 global forest dataset.\n",
        "\n",
        "# \"\"\" Created on Mon July 2 2024\n",
        "\n",
        "# @author: Pius N.Nwachukwu \"\"\"\n"
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "pip install earthengine-api matplotlib geopandas folium\n"
      ],
      "metadata": {
        "id": "LA4SoMET1eYi"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "# Import necessary libraries\n",
        "import ee\n",
        "import geopy\n",
        "import folium\n",
        "import pandas as pd\n",
        "import matplotlib.pyplot as plt\n",
        "from datetime import datetime\n",
        "from google.colab import auth\n",
        "from sklearn.metrics.pairwise import rbf_kernel\n",
        "import numpy as np\n",
        "from google.colab import drive\n",
        "import gdown\n",
        "from shapely.geometry import Polygon\n",
        "import geopandas as gpd\n"
      ],
      "metadata": {
        "id": "BeOLe0HBKI14"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "# Initialize Earth Engine\n",
        "auth.authenticate_user()\n",
        "ee.Initialize(project='ee-')"
      ],
      "metadata": {
        "id": "dqRYCIK0KJqD"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "source": [
        "\n",
        "# Function to aggregate pixel count\n",
        "def aggregate_pixel_count(product, geometries, start_date, end_date):\n",
        "    # Load MODIS ImageCollection\n",
        "    collection = ee.ImageCollection(product)\\\n",
        "        .filterDate(start_date, end_date)\n",
        "\n",
        "    # Calculate pixel count for each geometry\n",
        "    count_images = []\n",
        "    for geometry in geometries:\n",
        "        count_image = collection.filterBounds(geometry) \\\n",
        "            .reduce(ee.Reducer.count()) \\\n",
        "            .select([0], ['pixel_count']) \\\n",
        "            .clip(geometry)\n",
        "        count_images.append(count_image)\n",
        "\n",
        "    return count_images\n",
        "\n",
        "# Example usage:\n",
        "product = 'MODIS/061/MOD13A2'  # MODIS product path\n",
        "# Define rectangles for different continents\n",
        "continents = [\n",
        "    ee.Geometry.Rectangle([19.85,-3.16, 22.39,-1.31]),  # Africa\n",
        "    ee.Geometry.Rectangle([101.23,29.02, 102.10,29.54]),    # Asia\n",
        "    ee.Geometry.Rectangle([122.37,-33.08, 122.85,-32.49]),   # Aust_Ocean\n",
        "    ee.Geometry.Rectangle([22.45,45.12, 22.97,45.47]),  # Europe\n",
        "    ee.Geometry.Rectangle([-121.68,48.02, -119.80,49.01]),    # Nth_America\n",
        "    ee.Geometry.Rectangle([-60.70,-9.18, -58.42,-7.56]),   # Sth_America\n",
        "]\n",
        "start_date = '2020-01-01'\n",
        "end_date = '2020-12-30'\n",
        "\n",
        "pixel_count_images = aggregate_pixel_count(product, continents, start_date, end_date)\n",
        "\n",
        "# Define a function to add Earth Engine layers to a folium map\n",
        "def add_ee_layer(self, ee_image_object, vis_params, name):\n",
        "    map_id_dict = ee.Image(ee_image_object).getMapId(vis_params)\n",
        "    folium.raster_layers.TileLayer(\n",
        "        tiles=map_id_dict['tile_fetcher'].url_format,\n",
        "        attr='Google Earth Engine',\n",
        "        name=name,\n",
        "        overlay=True,\n",
        "        control=True\n",
        "    ).add_to(self)\n",
        "\n",
        "# Add EE drawing method to folium\n",
        "folium.Map.add_ee_layer = add_ee_layer\n",
        "\n",
        "# Create a folium map\n",
        "map_center = [20, 20]\n",
        "my_map = folium.Map(location=map_center, zoom_start=3)\n",
        "\n",
        "# Visualization parameters\n",
        "vis_params = {\n",
        "    'min': 0,\n",
        "    'max': 10,\n",
        "    'palette': ['blue', 'green', 'yellow', 'orange', 'red']\n",
        "}\n",
        "\n",
        "# Add the pixel count layers to the map\n",
        "for i, count_image in enumerate(pixel_count_images):\n",
        "    my_map.add_ee_layer(count_image, vis_params, f'Pixel Count {i+1}')\n",
        "\n",
        "# Add a layer control panel to the map\n",
        "my_map.add_child(folium.LayerControl())\n",
        "\n",
        "# Display the map\n",
        "my_map  # Or my_map.save(\"pixel_count_map.html\") to save it"
      ],
      "cell_type": "code",
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 674
        },
        "id": "CNM_109MbYIm",
        "outputId": "16b7115c-ca94-4d1c-92ba-efc26e1a01c8"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "<folium.folium.Map at 0x7d17fc5e19c0>"
            ],
            "text/html": [
              "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><span style=\"color:#565656\">Make this Notebook Trusted to load map: File -> Trust Notebook</span><iframe srcdoc=\"&lt;!DOCTYPE html&gt;\n",
              "&lt;html&gt;\n",
              "&lt;head&gt;\n",
              "    \n",
              "    &lt;meta http-equiv=&quot;content-type&quot; content=&quot;text/html; charset=UTF-8&quot; /&gt;\n",
              "    \n",
              "        &lt;script&gt;\n",
              "            L_NO_TOUCH = false;\n",
              "            L_DISABLE_3D = false;\n",
              "        &lt;/script&gt;\n",
              "    \n",
              "    &lt;style&gt;html, body {width: 100%;height: 100%;margin: 0;padding: 0;}&lt;/style&gt;\n",
              "    &lt;style&gt;#map {position:absolute;top:0;bottom:0;right:0;left:0;}&lt;/style&gt;\n",
              "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://code.jquery.com/jquery-1.12.4.min.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/js/bootstrap.bundle.min.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/css/bootstrap.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://netdna.bootstrapcdn.com/bootstrap/3.0.0/css/bootstrap.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/@fortawesome/fontawesome-free@6.2.0/css/all.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/gh/python-visualization/folium/folium/templates/leaflet.awesome.rotate.min.css&quot;/&gt;\n",
              "    \n",
              "            &lt;meta name=&quot;viewport&quot; content=&quot;width=device-width,\n",
              "                initial-scale=1.0, maximum-scale=1.0, user-scalable=no&quot; /&gt;\n",
              "            &lt;style&gt;\n",
              "                #map_88506dba3e37d0c84d220899ed1c3f18 {\n",
              "                    position: relative;\n",
              "                    width: 100.0%;\n",
              "                    height: 100.0%;\n",
              "                    left: 0.0%;\n",
              "                    top: 0.0%;\n",
              "                }\n",
              "                .leaflet-container { font-size: 1rem; }\n",
              "            &lt;/style&gt;\n",
              "        \n",
              "&lt;/head&gt;\n",
              "&lt;body&gt;\n",
              "    \n",
              "    \n",
              "            &lt;div class=&quot;folium-map&quot; id=&quot;map_88506dba3e37d0c84d220899ed1c3f18&quot; &gt;&lt;/div&gt;\n",
              "        \n",
              "&lt;/body&gt;\n",
              "&lt;script&gt;\n",
              "    \n",
              "    \n",
              "            var map_88506dba3e37d0c84d220899ed1c3f18 = L.map(\n",
              "                &quot;map_88506dba3e37d0c84d220899ed1c3f18&quot;,\n",
              "                {\n",
              "                    center: [20.0, 20.0],\n",
              "                    crs: L.CRS.EPSG3857,\n",
              "                    zoom: 3,\n",
              "                    zoomControl: true,\n",
              "                    preferCanvas: false,\n",
              "                }\n",
              "            );\n",
              "\n",
              "            \n",
              "\n",
              "        \n",
              "    \n",
              "            var tile_layer_8c8b232a72bbe131daf16798bdd4bd18 = L.tileLayer(\n",
              "                &quot;https://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Data by \\u0026copy; \\u003ca target=\\&quot;_blank\\&quot; href=\\&quot;http://openstreetmap.org\\&quot;\\u003eOpenStreetMap\\u003c/a\\u003e, under \\u003ca target=\\&quot;_blank\\&quot; href=\\&quot;http://www.openstreetmap.org/copyright\\&quot;\\u003eODbL\\u003c/a\\u003e.&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_52b780d4d9aa48792dd5b9b012abbec7 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/734d5bcea396a99b6eef038ab5343d73-637ec5db63ba17d42ff93bd755d6ec8e/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_328393ad0698fe2ac1b68a308d0e5b76 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/70b5d265a68f02f054886bb35f40b47f-104a74040078397b9cf709d2407563ab/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_6feff248935aa2da82baa46d3001330a = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/3b0be940298b5cd43a7bf8655190b784-ac9b3a3b2ad95469df0dea481fc6998c/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_d6b50773807c73dd7f8c99a0d6ce60c6 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/2ffceb15e0d3aa3d6c2ecb607eca85ab-1488b8a63c5a3862c378912d1e878807/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_3deb4fe552fb1e6b385db989a3b76226 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/81b561aed477752cded11a2cc1746bd9-04c915c55303b8b4e7ae2d8600449553/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_a9b71e91c05830699d0102e256d92a32 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/e9bc0eb17dfbdf3df05a499f02159326-a4301103c274f7176f9547760e7d255b/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var layer_control_ada72f09852909c29e558597ef57af77 = {\n",
              "                base_layers : {\n",
              "                    &quot;openstreetmap&quot; : tile_layer_8c8b232a72bbe131daf16798bdd4bd18,\n",
              "                },\n",
              "                overlays :  {\n",
              "                    &quot;Pixel Count 1&quot; : tile_layer_52b780d4d9aa48792dd5b9b012abbec7,\n",
              "                    &quot;Pixel Count 2&quot; : tile_layer_328393ad0698fe2ac1b68a308d0e5b76,\n",
              "                    &quot;Pixel Count 3&quot; : tile_layer_6feff248935aa2da82baa46d3001330a,\n",
              "                    &quot;Pixel Count 4&quot; : tile_layer_d6b50773807c73dd7f8c99a0d6ce60c6,\n",
              "                    &quot;Pixel Count 5&quot; : tile_layer_3deb4fe552fb1e6b385db989a3b76226,\n",
              "                    &quot;Pixel Count 6&quot; : tile_layer_a9b71e91c05830699d0102e256d92a32,\n",
              "                },\n",
              "            };\n",
              "            L.control.layers(\n",
              "                layer_control_ada72f09852909c29e558597ef57af77.base_layers,\n",
              "                layer_control_ada72f09852909c29e558597ef57af77.overlays,\n",
              "                {&quot;autoZIndex&quot;: true, &quot;collapsed&quot;: true, &quot;position&quot;: &quot;topright&quot;}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "&lt;/script&gt;\n",
              "&lt;/html&gt;\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
            ]
          },
          "metadata": {},
          "execution_count": 51
        }
      ]
    },
    {
      "cell_type": "markdown",
      "source": [
        "# Or this"
      ],
      "metadata": {
        "id": "hlQb_Us2EUwl"
      }
    },
    {
      "source": [
        "\n",
        "# Function to aggregate pixel count\n",
        "def aggregate_pixel_count(product, geometries, start_date, end_date):\n",
        "    # Load MODIS ImageCollection\n",
        "    collection = ee.ImageCollection(product)\\\n",
        "        .filterDate(start_date, end_date)\n",
        "\n",
        "    # Calculate pixel count for each geometry\n",
        "    count_images = []\n",
        "    for geometry in geometries:\n",
        "        count_image = collection.filterBounds(geometry) \\\n",
        "            .reduce(ee.Reducer.count()) \\\n",
        "            .select([0], ['pixel_count']) \\\n",
        "            .clip(geometry)\n",
        "        count_images.append(count_image)\n",
        "\n",
        "    return count_images\n",
        "\n",
        "# Example usage:\n",
        "product = 'MODIS/061/MOD13A2'  # MODIS product path\n",
        "# Define rectangles for different continents\n",
        "continents = [\n",
        "    ee.Geometry.Rectangle([19.85,-3.16, 22.39,-1.31]),  # Africa\n",
        "    ee.Geometry.Rectangle([101.23,29.02, 102.10,29.54]),    # Asia\n",
        "    ee.Geometry.Rectangle([122.37,-33.08, 122.85,-32.49]),   # Aust_Ocean\n",
        "    ee.Geometry.Rectangle([22.45,45.12, 22.97,45.47]),  # Europe\n",
        "    ee.Geometry.Rectangle([-121.68,48.02, -119.80,49.01]),    # Nth_America\n",
        "    ee.Geometry.Rectangle([-60.70,-9.18, -58.42,-7.56]),   # Sth_America\n",
        "\n",
        "]\n",
        "start_date = '2020-02-01'\n",
        "end_date = '2020-04-30'\n",
        "\n",
        "pixel_count_images = aggregate_pixel_count(product, continents, start_date, end_date)\n",
        "\n",
        "# Define a function to add Earth Engine layers to a folium map\n",
        "def add_ee_layer(self, ee_image_object, vis_params, name, label): # Added label parameter\n",
        "    map_id_dict = ee.Image(ee_image_object).getMapId(vis_params)\n",
        "    folium.raster_layers.TileLayer(\n",
        "        tiles=map_id_dict['tile_fetcher'].url_format,\n",
        "        attr='Google Earth Engine',\n",
        "        name=name,\n",
        "        overlay=True,\n",
        "        control=True,\n",
        "        show=False # Initially hide the layers\n",
        "    ).add_to(self)\n",
        "\n",
        "    # Add a label for the layer\n",
        "    folium.Marker(\n",
        "        location=ee_image_object.geometry().centroid().coordinates().getInfo()[::-1], # Get the center of the geometry\n",
        "        popup=label,\n",
        "        icon=folium.Icon(color='blue')\n",
        "    ).add_to(self)\n",
        "\n",
        "# Add EE drawing method to folium\n",
        "folium.Map.add_ee_layer = add_ee_layer\n",
        "\n",
        "# Create a folium map\n",
        "map_center = [20, 20]  # Adjust for a global view\n",
        "my_map = folium.Map(location=map_center, zoom_start=3)\n",
        "\n",
        "# Visualization parameters\n",
        "vis_params = {\n",
        "    'min': 0,\n",
        "    'max': 10,  # Adjust this based on expected range of pixel counts\n",
        "    'palette': ['blue', 'green', 'yellow', 'orange', 'red']\n",
        "}\n",
        "\n",
        "# Labels for each continent\n",
        "continent_labels = [\n",
        "    'Africa', 'Asia', 'Australia/Oceania', 'Europe', 'North America', 'South America'\n",
        "]\n",
        "\n",
        "# Add the pixel count layers to the map with labels\n",
        "for i, count_image in enumerate(pixel_count_images):\n",
        "    my_map.add_ee_layer(count_image, vis_params, f'Pixel Count {i+1}', continent_labels[i]) # Pass the label from continent_labels\n",
        "\n",
        "# Add a layer control panel to the map\n",
        "my_map.add_child(folium.LayerControl())\n",
        "\n",
        "# Display the map\n",
        "my_map  # Or my_map.save(\"pixel_count_map.html\") to save it"
      ],
      "cell_type": "code",
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 674
        },
        "id": "YafmWDFr6rlh",
        "outputId": "47f63fbe-f153-469d-d344-82834850f749"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "<folium.folium.Map at 0x7d17fcadfc10>"
            ],
            "text/html": [
              "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><span style=\"color:#565656\">Make this Notebook Trusted to load map: File -> Trust Notebook</span><iframe srcdoc=\"&lt;!DOCTYPE html&gt;\n",
              "&lt;html&gt;\n",
              "&lt;head&gt;\n",
              "    \n",
              "    &lt;meta http-equiv=&quot;content-type&quot; content=&quot;text/html; charset=UTF-8&quot; /&gt;\n",
              "    \n",
              "        &lt;script&gt;\n",
              "            L_NO_TOUCH = false;\n",
              "            L_DISABLE_3D = false;\n",
              "        &lt;/script&gt;\n",
              "    \n",
              "    &lt;style&gt;html, body {width: 100%;height: 100%;margin: 0;padding: 0;}&lt;/style&gt;\n",
              "    &lt;style&gt;#map {position:absolute;top:0;bottom:0;right:0;left:0;}&lt;/style&gt;\n",
              "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://code.jquery.com/jquery-1.12.4.min.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/js/bootstrap.bundle.min.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/css/bootstrap.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://netdna.bootstrapcdn.com/bootstrap/3.0.0/css/bootstrap.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/@fortawesome/fontawesome-free@6.2.0/css/all.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/gh/python-visualization/folium/folium/templates/leaflet.awesome.rotate.min.css&quot;/&gt;\n",
              "    \n",
              "            &lt;meta name=&quot;viewport&quot; content=&quot;width=device-width,\n",
              "                initial-scale=1.0, maximum-scale=1.0, user-scalable=no&quot; /&gt;\n",
              "            &lt;style&gt;\n",
              "                #map_69d08b70ea0decdd9797229d0bfb624b {\n",
              "                    position: relative;\n",
              "                    width: 100.0%;\n",
              "                    height: 100.0%;\n",
              "                    left: 0.0%;\n",
              "                    top: 0.0%;\n",
              "                }\n",
              "                .leaflet-container { font-size: 1rem; }\n",
              "            &lt;/style&gt;\n",
              "        \n",
              "&lt;/head&gt;\n",
              "&lt;body&gt;\n",
              "    \n",
              "    \n",
              "            &lt;div class=&quot;folium-map&quot; id=&quot;map_69d08b70ea0decdd9797229d0bfb624b&quot; &gt;&lt;/div&gt;\n",
              "        \n",
              "&lt;/body&gt;\n",
              "&lt;script&gt;\n",
              "    \n",
              "    \n",
              "            var map_69d08b70ea0decdd9797229d0bfb624b = L.map(\n",
              "                &quot;map_69d08b70ea0decdd9797229d0bfb624b&quot;,\n",
              "                {\n",
              "                    center: [20.0, 20.0],\n",
              "                    crs: L.CRS.EPSG3857,\n",
              "                    zoom: 3,\n",
              "                    zoomControl: true,\n",
              "                    preferCanvas: false,\n",
              "                }\n",
              "            );\n",
              "\n",
              "            \n",
              "\n",
              "        \n",
              "    \n",
              "            var tile_layer_2374a722904e935627241e89a5e81510 = L.tileLayer(\n",
              "                &quot;https://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Data by \\u0026copy; \\u003ca target=\\&quot;_blank\\&quot; href=\\&quot;http://openstreetmap.org\\&quot;\\u003eOpenStreetMap\\u003c/a\\u003e, under \\u003ca target=\\&quot;_blank\\&quot; href=\\&quot;http://www.openstreetmap.org/copyright\\&quot;\\u003eODbL\\u003c/a\\u003e.&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var tile_layer_b2f2f840c82df8a1b7658dd09c21ecd6 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/316b0b920a12bb241dbbcbdec2c07829-94add5e7720c5e518fc8a26b47b33e05/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_187b9c01a122896ebc703302e29709f7 = L.marker(\n",
              "                [-2.235354019871469, 21.119999999999973],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_3be35ef941368229fe39fc0b28f291c7 = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_187b9c01a122896ebc703302e29709f7.setIcon(icon_3be35ef941368229fe39fc0b28f291c7);\n",
              "        \n",
              "    \n",
              "        var popup_3f383803add628cc1f55cacaf1a9f3d1 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_469532255e4128bf0fe45e9fd1a49ed1 = $(`&lt;div id=&quot;html_469532255e4128bf0fe45e9fd1a49ed1&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Africa&lt;/div&gt;`)[0];\n",
              "                popup_3f383803add628cc1f55cacaf1a9f3d1.setContent(html_469532255e4128bf0fe45e9fd1a49ed1);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_187b9c01a122896ebc703302e29709f7.bindPopup(popup_3f383803add628cc1f55cacaf1a9f3d1)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var tile_layer_193f5e7ca9cb17ca77178769bb92533d = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/005a562f519482986bd8eedd47dbac4f-e9f90ae193d82d64e8ca1068864c3c7b/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_5769fdc1253a610b828b32ab972d55d3 = L.marker(\n",
              "                [29.280483897479627, 101.66499999999978],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_e5dc8127e1c05b10f0e0eb5793806ffe = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_5769fdc1253a610b828b32ab972d55d3.setIcon(icon_e5dc8127e1c05b10f0e0eb5793806ffe);\n",
              "        \n",
              "    \n",
              "        var popup_19fbc22d8c630000f643626fb0ab124e = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_9a79bcf6e5e799f4759211fded558521 = $(`&lt;div id=&quot;html_9a79bcf6e5e799f4759211fded558521&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Asia&lt;/div&gt;`)[0];\n",
              "                popup_19fbc22d8c630000f643626fb0ab124e.setContent(html_9a79bcf6e5e799f4759211fded558521);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_5769fdc1253a610b828b32ab972d55d3.bindPopup(popup_19fbc22d8c630000f643626fb0ab124e)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var tile_layer_258559f0c534a78882dbd6b8f77e269b = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/345dfaca4bd38f08da3d7d44eb2ebc0a-7ead80d5e960ac00747818bd16f62339/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_1b5e7dda867e27af29eb84aff6c699fe = L.marker(\n",
              "                [-32.7849027186664, 122.61000000000037],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_54276d2725d83abd203d121f0a12b43f = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_1b5e7dda867e27af29eb84aff6c699fe.setIcon(icon_54276d2725d83abd203d121f0a12b43f);\n",
              "        \n",
              "    \n",
              "        var popup_b162901fe9faeefeb7f7e07ba1d7e5da = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_574140113edf6905777522015f9b9527 = $(`&lt;div id=&quot;html_574140113edf6905777522015f9b9527&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Australia/Oceania&lt;/div&gt;`)[0];\n",
              "                popup_b162901fe9faeefeb7f7e07ba1d7e5da.setContent(html_574140113edf6905777522015f9b9527);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_1b5e7dda867e27af29eb84aff6c699fe.bindPopup(popup_b162901fe9faeefeb7f7e07ba1d7e5da)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var tile_layer_199fa9e5bebdc9b29bd6e02fed913575 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/13ce548180f932639f88663c4822fd63-3ce10d84e5609784df0b07e2552ac5ef/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_f603ee0fbbee9d00d4ee49e8f855c60c = L.marker(\n",
              "                [45.295114927680466, 22.710000000000402],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_a6b432ca6fe5be9d826a90d0b651c3c2 = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_f603ee0fbbee9d00d4ee49e8f855c60c.setIcon(icon_a6b432ca6fe5be9d826a90d0b651c3c2);\n",
              "        \n",
              "    \n",
              "        var popup_8000889f0b304991ca6634f68c06d0fd = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_84dd324758be8caca7e8c0cb7f7ea19c = $(`&lt;div id=&quot;html_84dd324758be8caca7e8c0cb7f7ea19c&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Europe&lt;/div&gt;`)[0];\n",
              "                popup_8000889f0b304991ca6634f68c06d0fd.setContent(html_84dd324758be8caca7e8c0cb7f7ea19c);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_f603ee0fbbee9d00d4ee49e8f855c60c.bindPopup(popup_8000889f0b304991ca6634f68c06d0fd)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var tile_layer_26fad73cf6c2104eefe1648fa6d9009c = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/d97624f199747268b38ae7f56b11e99c-639cbf94b36befdd274a34d54f08b315/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_f8ade20a1b8ef636bfe0f09f275d5a97 = L.marker(\n",
              "                [48.51721400867622, -120.73999999999994],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_1af20e2c038078735d6bfa2e30c99481 = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_f8ade20a1b8ef636bfe0f09f275d5a97.setIcon(icon_1af20e2c038078735d6bfa2e30c99481);\n",
              "        \n",
              "    \n",
              "        var popup_c34d9def621c7a408f69698cff1f5d35 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_813f6d354fb6872a9f5b1700a8c495a6 = $(`&lt;div id=&quot;html_813f6d354fb6872a9f5b1700a8c495a6&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;North America&lt;/div&gt;`)[0];\n",
              "                popup_c34d9def621c7a408f69698cff1f5d35.setContent(html_813f6d354fb6872a9f5b1700a8c495a6);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_f8ade20a1b8ef636bfe0f09f275d5a97.bindPopup(popup_c34d9def621c7a408f69698cff1f5d35)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var tile_layer_827b458f30bd1bc3ee96822c7fbaaac2 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/2cea209f6deb604adb85c69bfe594257-2e01d8c92410b0340b3624772e559b8f/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_64d262330bfcc9ef9730e9cf42a8d2ab = L.marker(\n",
              "                [-8.371071242702605, -59.55999999999991],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_eecceb7d0e032353bfdfb178d8dafb79 = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_64d262330bfcc9ef9730e9cf42a8d2ab.setIcon(icon_eecceb7d0e032353bfdfb178d8dafb79);\n",
              "        \n",
              "    \n",
              "        var popup_9701b3437217e8fcd80f57489d3547a1 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_305b51c966b7f76eb15e3de83ba5411d = $(`&lt;div id=&quot;html_305b51c966b7f76eb15e3de83ba5411d&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;South America&lt;/div&gt;`)[0];\n",
              "                popup_9701b3437217e8fcd80f57489d3547a1.setContent(html_305b51c966b7f76eb15e3de83ba5411d);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_64d262330bfcc9ef9730e9cf42a8d2ab.bindPopup(popup_9701b3437217e8fcd80f57489d3547a1)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var layer_control_1417db652dc8eacbeb40497b221ef399 = {\n",
              "                base_layers : {\n",
              "                    &quot;openstreetmap&quot; : tile_layer_2374a722904e935627241e89a5e81510,\n",
              "                },\n",
              "                overlays :  {\n",
              "                    &quot;Pixel Count 1&quot; : tile_layer_b2f2f840c82df8a1b7658dd09c21ecd6,\n",
              "                    &quot;Pixel Count 2&quot; : tile_layer_193f5e7ca9cb17ca77178769bb92533d,\n",
              "                    &quot;Pixel Count 3&quot; : tile_layer_258559f0c534a78882dbd6b8f77e269b,\n",
              "                    &quot;Pixel Count 4&quot; : tile_layer_199fa9e5bebdc9b29bd6e02fed913575,\n",
              "                    &quot;Pixel Count 5&quot; : tile_layer_26fad73cf6c2104eefe1648fa6d9009c,\n",
              "                    &quot;Pixel Count 6&quot; : tile_layer_827b458f30bd1bc3ee96822c7fbaaac2,\n",
              "                },\n",
              "            };\n",
              "            L.control.layers(\n",
              "                layer_control_1417db652dc8eacbeb40497b221ef399.base_layers,\n",
              "                layer_control_1417db652dc8eacbeb40497b221ef399.overlays,\n",
              "                {&quot;autoZIndex&quot;: true, &quot;collapsed&quot;: true, &quot;position&quot;: &quot;topright&quot;}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "            tile_layer_b2f2f840c82df8a1b7658dd09c21ecd6.remove();\n",
              "            tile_layer_193f5e7ca9cb17ca77178769bb92533d.remove();\n",
              "            tile_layer_258559f0c534a78882dbd6b8f77e269b.remove();\n",
              "            tile_layer_199fa9e5bebdc9b29bd6e02fed913575.remove();\n",
              "            tile_layer_26fad73cf6c2104eefe1648fa6d9009c.remove();\n",
              "            tile_layer_827b458f30bd1bc3ee96822c7fbaaac2.remove();\n",
              "        \n",
              "&lt;/script&gt;\n",
              "&lt;/html&gt;\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
            ]
          },
          "metadata": {},
          "execution_count": 52
        }
      ]
    },
    {
      "source": [
        "# Use this to print the pixel count\n",
        "\n",
        "def aggregate_pixel_count(product, geometry, start_date, end_date):\n",
        "    collection = ee.ImageCollection(product) \\\n",
        "        .filterBounds(geometry) \\\n",
        "        .filterDate(start_date, end_date)\n",
        "\n",
        "    count_image = collection.reduce(ee.Reducer.count()).select([0], ['pixel_count'])\n",
        "    count_image = count_image.clip(geometry)\n",
        "    return count_image\n",
        "\n",
        "# Example usage:\n",
        "product = 'MODIS/061/MOD13A2'\n",
        "start_date = '2000-02-01'\n",
        "end_date = '2024-04-30'\n",
        "\n",
        "# Define geometries for each continent\n",
        "geometries = {\n",
        "    'Africa': ee.Geometry.Rectangle([19.85,-3.16, 22.39,-1.31]),\n",
        "    'Asia': ee.Geometry.Rectangle([101.23,29.02, 102.10,29.54]),\n",
        "    'Aust_Ocean': ee.Geometry.Rectangle([122.37,-33.08, 122.85,-32.49]),\n",
        "    'Europe': ee.Geometry.Rectangle([22.45,45.12, 22.97,45.47]),\n",
        "    'Nth_America': ee.Geometry.Rectangle([-121.68,48.02, -119.80,49.01]),\n",
        "    'Sth_America': ee.Geometry.Rectangle([-60.70,-9.18, -58.42,-7.56]),\n",
        "}\n",
        "\n",
        "pixel_counts = {}\n",
        "for continent, geometry in geometries.items():\n",
        "    pixel_count_image = aggregate_pixel_count(product, geometry, start_date, end_date)\n",
        "    # Get the total pixel count for the continent\n",
        "    pixel_count = pixel_count_image.reduceRegion(\n",
        "        reducer=ee.Reducer.sum(),\n",
        "        geometry=geometry,\n",
        "        scale=500,\n",
        "        maxPixels=1e13  # Increase the maximum number of pixels allowed\n",
        "        #bestEffort=True\n",
        "    ).get('pixel_count').getInfo()\n",
        "    pixel_counts[continent] = pixel_count\n",
        "\n",
        "# Print the pixel counts for each continent\n",
        "for continent, count in pixel_counts.items():\n",
        "    print(f'Pixel count for {continent}: {count}')"
      ],
      "cell_type": "code",
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "9JjM6Wz2cZek",
        "outputId": "d5be6c80-daac-4e57-dc2d-49fda02153b9"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Pixel count for Africa: 129715142.91764683\n",
            "Pixel count for Asia: 12479566.815686226\n",
            "Pixel count for Aust_Ocean: 7814286.243137279\n",
            "Pixel count for Europe: 5020308.470588274\n",
            "Pixel count for Nth_America: 51338842.070588194\n",
            "Pixel count for Sth_America: 101914020.35686333\n"
          ]
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "pip install earthengine-api matplotlib geopandas folium\n"
      ],
      "metadata": {
        "id": "HAqH0cxHGLxQ"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "# Import necessary libraries\n",
        "import ee\n",
        "import geopy\n",
        "import folium\n",
        "import pandas as pd\n",
        "import matplotlib.pyplot as plt\n",
        "from datetime import datetime\n",
        "from google.colab import auth\n",
        "from sklearn.metrics.pairwise import rbf_kernel\n",
        "import numpy as np\n",
        "from google.colab import drive\n",
        "import gdown\n",
        "# Import the Polygon class from shapely.geometry\n",
        "from shapely.geometry import Polygon\n",
        "import geopandas as gpd\n"
      ],
      "metadata": {
        "id": "5NpGn4mqGLxR"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "# Initialize Earth Engine\n",
        "auth.authenticate_user()\n",
        "ee.Initialize(project='ee-pnnwachukwu17')"
      ],
      "metadata": {
        "id": "51tSBUxaGLxS"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "source": [
        " # For single continent like Europe"
      ],
      "metadata": {
        "id": "qxicYoqe1upb"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "# For single continent like Europe\n",
        "\n",
        "# Initialize Earth Engine\n",
        "# auth.authenticate_user()  # Assuming you have already authenticated\n",
        "ee.Initialize(project='ee-pnnwachukwu17')\n",
        "\n",
        "# Function to aggregate pixel count\n",
        "def aggregate_pixel_count(product, geometry, start_date, end_date):\n",
        "    # Load MODIS ImageCollection\n",
        "    collection = ee.ImageCollection(product) \\\n",
        "        .filterBounds(geometry) \\\n",
        "        .filterDate(start_date, end_date)\n",
        "\n",
        "    # Reduce the collection to get the total pixel count in the region of interest\n",
        "    count_image = collection.reduce(ee.Reducer.count()).select([0], ['pixel_count']) # Select the first band and rename it\n",
        "\n",
        "    # Clip to the geometry\n",
        "    count_image = count_image.clip(geometry)\n",
        "\n",
        "    return count_image\n",
        "\n",
        "# Example usage:\n",
        "product = 'MODIS/061/MOD13A2'  # MODIS product path\n",
        "# Define a square geometry (replace with your desired coordinates)\n",
        "geometry = ee.Geometry.Rectangle([\n",
        "    [-122.1, 37.4],  # Bottom left corner\n",
        "    [-122.0, 37.5]   # Top right corner\n",
        "])\n",
        "start_date = '2020-01-01'\n",
        "end_date = '2020-12-31'\n",
        "\n",
        "pixel_count_image = aggregate_pixel_count(product, geometry, start_date, end_date)\n",
        "\n",
        "# Define a function to add Earth Engine layers to a folium map\n",
        "def add_ee_layer(self, ee_image_object, vis_params, name):\n",
        "    map_id_dict = ee.Image(ee_image_object).getMapId(vis_params)\n",
        "    folium.raster_layers.TileLayer(\n",
        "        tiles=map_id_dict['tile_fetcher'].url_format,\n",
        "        attr='Google Earth Engine',\n",
        "        name=name,\n",
        "        overlay=True,\n",
        "        control=True\n",
        "    ).add_to(self)\n",
        "\n",
        "# Add EE drawing method to folium\n",
        "folium.Map.add_ee_layer = add_ee_layer\n",
        "\n",
        "# Create a folium map\n",
        "map_center = geometry.centroid().coordinates().getInfo()[::-1]\n",
        "my_map = folium.Map(location=map_center, zoom_start=12)\n",
        "\n",
        "# Visualization parameters\n",
        "vis_params = {\n",
        "    'min': 0,\n",
        "    'max': 10,  # Adjust this based on expected range of pixel counts\n",
        "    'palette': ['blue', 'green', 'yellow', 'orange', 'red']\n",
        "}\n",
        "\n",
        "# Add the pixel count layer to the map\n",
        "my_map.add_ee_layer(pixel_count_image, vis_params, 'Pixel Count')\n",
        "\n",
        "# Add a layer control panel to the map\n",
        "my_map.add_child(folium.LayerControl())\n",
        "\n",
        "# Display the map\n",
        "my_map  # Or my_map.save(\"pixel_count_map.html\") to save it\n",
        "\n"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "cdi7lFymZROg",
        "outputId": "2b07f449-1536-494e-d5d3-9947c0e8be9e"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Total number of pixels: None\n"
          ]
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [],
      "metadata": {
        "id": "-_QGw8rvYsan"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "# This is used to print the number of pixels\n",
        "\n",
        "# Function to aggregate pixel count\n",
        "def aggregate_pixel_count(product, geometry, start_date, end_date):\n",
        "    # Load MODIS ImageCollection\n",
        "    collection = ee.ImageCollection(product) \\\n",
        "        .filterBounds(geometry) \\\n",
        "        .filterDate(start_date, end_date)\n",
        "\n",
        "    # Reduce the collection to get the total pixel count in the region of interest\n",
        "    count_image = collection.reduce(ee.Reducer.count()).select([0], ['pixel_count'])\n",
        "\n",
        "    # Clip to the geometry\n",
        "    count_image = count_image.clip(geometry)\n",
        "\n",
        "    return count_image\n",
        "\n",
        "# Example usage:\n",
        "product = 'MODIS/061/MOD13A2'  # MODIS product path\n",
        "\n",
        "# Define a square geometry (replace with your desired coordinates)\n",
        "geometry = ee.Geometry.Rectangle([\n",
        "    22.45,45.12,  # Lower left corner (longitude, latitude)\n",
        "    22.97,45.47   # Upper right corner (longitude, latitude)\n",
        "])\n",
        "\n",
        "start_date = '2020-01-01'\n",
        "end_date = '2020-12-31'\n",
        "\n",
        "pixel_count_image = aggregate_pixel_count(product, geometry, start_date, end_date)\n",
        "\n",
        "# Get the total number of pixels\n",
        "pixel_count = pixel_count_image.reduceRegion(\n",
        "    reducer=ee.Reducer.sum(),\n",
        "    geometry=geometry,\n",
        "    scale=500  # Adjust scale as needed\n",
        ").get('pixel_count').getInfo()\n",
        "\n",
        "print(\"Total number of pixels:\", pixel_count)\n",
        "\n",
        "# ... (rest of the code for creating and displaying the map)"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "Xy5IMW2NYBoo",
        "outputId": "d9b02b6f-10e9-4c87-9a3e-5b698b851044"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Total number of pixels: 207374.13333333423\n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "source": [
        " # For Multiple Continents"
      ],
      "metadata": {
        "id": "p0L8hOfJae1A"
      }
    },
    {
      "cell_type": "code",
      "source": [],
      "metadata": {
        "id": "dgo7rX6lakC8"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "source": [
        "import ee\n",
        "import folium\n",
        "import matplotlib.pyplot as plt\n",
        "\n",
        "# Initialize Earth Engine\n",
        "# auth.authenticate_user()  # Assuming you have already authenticated\n",
        "ee.Initialize(project='ee-pnnwachukwu17')\n",
        "\n",
        "# Function to aggregate pixel count\n",
        "def aggregate_pixel_count(product, geometries, start_date, end_date):\n",
        "    # Load MODIS ImageCollection\n",
        "    collection = ee.ImageCollection(product)\\\n",
        "        .filterDate(start_date, end_date)\n",
        "\n",
        "    # Calculate pixel count for each geometry\n",
        "    count_images = []\n",
        "    for geometry in geometries:\n",
        "        count_image = collection.filterBounds(geometry) \\\n",
        "            .reduce(ee.Reducer.count()) \\\n",
        "            .select([0], ['pixel_count']) \\\n",
        "            .clip(geometry)\n",
        "        count_images.append(count_image)\n",
        "\n",
        "    return count_images\n",
        "\n",
        "# Example usage:\n",
        "product = 'MODIS/061/MOD13A2'  # MODIS product path\n",
        "# Define rectangles for different continents (replace with your desired coordinates)\n",
        "continents = [\n",
        "    ee.Geometry.Rectangle([19.85,-3.16, 22.39,-1.31]),  # Africa (approx.)\n",
        "    ee.Geometry.Rectangle([101.23,29.02, 102.10,29.54]),    # Asia (approx.)\n",
        "    ee.Geometry.Rectangle([122.37,-33.08, 122.85,-32.49]),   # Aust_Ocean (approx.)\n",
        "    ee.Geometry.Rectangle([22.45,45.12, 22.97,45.47]),  # Europe (approx.)\n",
        "    ee.Geometry.Rectangle([-121.68,48.02, -119.80,49.01]),    # Nth_America (approx.)\n",
        "    ee.Geometry.Rectangle([-60.70,-9.18, -58.42,-7.56]),   # Sth_America (approx.)\n",
        "    # Add more rectangles for other continents\n",
        "]\n",
        "start_date = '2020-01-01'\n",
        "end_date = '2020-12-30'\n",
        "\n",
        "pixel_count_images = aggregate_pixel_count(product, continents, start_date, end_date)\n",
        "\n",
        "# Define a function to add Earth Engine layers to a folium map\n",
        "def add_ee_layer(self, ee_image_object, vis_params, name):\n",
        "    map_id_dict = ee.Image(ee_image_object).getMapId(vis_params)\n",
        "    folium.raster_layers.TileLayer(\n",
        "        tiles=map_id_dict['tile_fetcher'].url_format,\n",
        "        attr='Google Earth Engine',\n",
        "        name=name,\n",
        "        overlay=True,\n",
        "        control=True\n",
        "    ).add_to(self)\n",
        "\n",
        "# Add EE drawing method to folium\n",
        "folium.Map.add_ee_layer = add_ee_layer\n",
        "\n",
        "# Create a folium map (adjust map center as needed)\n",
        "map_center = [20, 20]  # Adjust for a global view\n",
        "my_map = folium.Map(location=map_center, zoom_start=3)\n",
        "\n",
        "# Visualization parameters\n",
        "vis_params = {\n",
        "    'min': 0,\n",
        "    'max': 10,  # Adjust this based on expected range of pixel counts\n",
        "    'palette': ['blue', 'green', 'yellow', 'orange', 'red']\n",
        "}\n",
        "\n",
        "# Add the pixel count layers to the map\n",
        "for i, count_image in enumerate(pixel_count_images):\n",
        "    my_map.add_ee_layer(count_image, vis_params, f'Pixel Count {i+1}')\n",
        "\n",
        "# Add a layer control panel to the map\n",
        "my_map.add_child(folium.LayerControl())\n",
        "\n",
        "# Display the map\n",
        "my_map  # Or my_map.save(\"pixel_count_map.html\") to save it"
      ],
      "cell_type": "code",
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 674
        },
        "outputId": "16b7115c-ca94-4d1c-92ba-efc26e1a01c8",
        "id": "Ibu65nfhGLxZ"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "<folium.folium.Map at 0x7d17fc5e19c0>"
            ],
            "text/html": [
              "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><span style=\"color:#565656\">Make this Notebook Trusted to load map: File -> Trust Notebook</span><iframe srcdoc=\"&lt;!DOCTYPE html&gt;\n",
              "&lt;html&gt;\n",
              "&lt;head&gt;\n",
              "    \n",
              "    &lt;meta http-equiv=&quot;content-type&quot; content=&quot;text/html; charset=UTF-8&quot; /&gt;\n",
              "    \n",
              "        &lt;script&gt;\n",
              "            L_NO_TOUCH = false;\n",
              "            L_DISABLE_3D = false;\n",
              "        &lt;/script&gt;\n",
              "    \n",
              "    &lt;style&gt;html, body {width: 100%;height: 100%;margin: 0;padding: 0;}&lt;/style&gt;\n",
              "    &lt;style&gt;#map {position:absolute;top:0;bottom:0;right:0;left:0;}&lt;/style&gt;\n",
              "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://code.jquery.com/jquery-1.12.4.min.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/js/bootstrap.bundle.min.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/css/bootstrap.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://netdna.bootstrapcdn.com/bootstrap/3.0.0/css/bootstrap.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/@fortawesome/fontawesome-free@6.2.0/css/all.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/gh/python-visualization/folium/folium/templates/leaflet.awesome.rotate.min.css&quot;/&gt;\n",
              "    \n",
              "            &lt;meta name=&quot;viewport&quot; content=&quot;width=device-width,\n",
              "                initial-scale=1.0, maximum-scale=1.0, user-scalable=no&quot; /&gt;\n",
              "            &lt;style&gt;\n",
              "                #map_88506dba3e37d0c84d220899ed1c3f18 {\n",
              "                    position: relative;\n",
              "                    width: 100.0%;\n",
              "                    height: 100.0%;\n",
              "                    left: 0.0%;\n",
              "                    top: 0.0%;\n",
              "                }\n",
              "                .leaflet-container { font-size: 1rem; }\n",
              "            &lt;/style&gt;\n",
              "        \n",
              "&lt;/head&gt;\n",
              "&lt;body&gt;\n",
              "    \n",
              "    \n",
              "            &lt;div class=&quot;folium-map&quot; id=&quot;map_88506dba3e37d0c84d220899ed1c3f18&quot; &gt;&lt;/div&gt;\n",
              "        \n",
              "&lt;/body&gt;\n",
              "&lt;script&gt;\n",
              "    \n",
              "    \n",
              "            var map_88506dba3e37d0c84d220899ed1c3f18 = L.map(\n",
              "                &quot;map_88506dba3e37d0c84d220899ed1c3f18&quot;,\n",
              "                {\n",
              "                    center: [20.0, 20.0],\n",
              "                    crs: L.CRS.EPSG3857,\n",
              "                    zoom: 3,\n",
              "                    zoomControl: true,\n",
              "                    preferCanvas: false,\n",
              "                }\n",
              "            );\n",
              "\n",
              "            \n",
              "\n",
              "        \n",
              "    \n",
              "            var tile_layer_8c8b232a72bbe131daf16798bdd4bd18 = L.tileLayer(\n",
              "                &quot;https://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Data by \\u0026copy; \\u003ca target=\\&quot;_blank\\&quot; href=\\&quot;http://openstreetmap.org\\&quot;\\u003eOpenStreetMap\\u003c/a\\u003e, under \\u003ca target=\\&quot;_blank\\&quot; href=\\&quot;http://www.openstreetmap.org/copyright\\&quot;\\u003eODbL\\u003c/a\\u003e.&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_52b780d4d9aa48792dd5b9b012abbec7 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/734d5bcea396a99b6eef038ab5343d73-637ec5db63ba17d42ff93bd755d6ec8e/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_328393ad0698fe2ac1b68a308d0e5b76 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/70b5d265a68f02f054886bb35f40b47f-104a74040078397b9cf709d2407563ab/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_6feff248935aa2da82baa46d3001330a = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/3b0be940298b5cd43a7bf8655190b784-ac9b3a3b2ad95469df0dea481fc6998c/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_d6b50773807c73dd7f8c99a0d6ce60c6 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/2ffceb15e0d3aa3d6c2ecb607eca85ab-1488b8a63c5a3862c378912d1e878807/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_3deb4fe552fb1e6b385db989a3b76226 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/81b561aed477752cded11a2cc1746bd9-04c915c55303b8b4e7ae2d8600449553/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_a9b71e91c05830699d0102e256d92a32 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/e9bc0eb17dfbdf3df05a499f02159326-a4301103c274f7176f9547760e7d255b/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var layer_control_ada72f09852909c29e558597ef57af77 = {\n",
              "                base_layers : {\n",
              "                    &quot;openstreetmap&quot; : tile_layer_8c8b232a72bbe131daf16798bdd4bd18,\n",
              "                },\n",
              "                overlays :  {\n",
              "                    &quot;Pixel Count 1&quot; : tile_layer_52b780d4d9aa48792dd5b9b012abbec7,\n",
              "                    &quot;Pixel Count 2&quot; : tile_layer_328393ad0698fe2ac1b68a308d0e5b76,\n",
              "                    &quot;Pixel Count 3&quot; : tile_layer_6feff248935aa2da82baa46d3001330a,\n",
              "                    &quot;Pixel Count 4&quot; : tile_layer_d6b50773807c73dd7f8c99a0d6ce60c6,\n",
              "                    &quot;Pixel Count 5&quot; : tile_layer_3deb4fe552fb1e6b385db989a3b76226,\n",
              "                    &quot;Pixel Count 6&quot; : tile_layer_a9b71e91c05830699d0102e256d92a32,\n",
              "                },\n",
              "            };\n",
              "            L.control.layers(\n",
              "                layer_control_ada72f09852909c29e558597ef57af77.base_layers,\n",
              "                layer_control_ada72f09852909c29e558597ef57af77.overlays,\n",
              "                {&quot;autoZIndex&quot;: true, &quot;collapsed&quot;: true, &quot;position&quot;: &quot;topright&quot;}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "&lt;/script&gt;\n",
              "&lt;/html&gt;\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
            ]
          },
          "metadata": {},
          "execution_count": 51
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [],
      "metadata": {
        "id": "x18HndNNcMmT"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [],
      "metadata": {
        "id": "GTuhjinmcYlE"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "source": [
        "import ee\n",
        "import folium\n",
        "import matplotlib.pyplot as plt\n",
        "\n",
        "# Initialize Earth Engine\n",
        "# auth.authenticate_user()  # Assuming you have already authenticated\n",
        "ee.Initialize(project='ee-pnnwachukwu17')\n",
        "\n",
        "# Function to aggregate pixel count\n",
        "def aggregate_pixel_count(product, geometries, start_date, end_date):\n",
        "    # Load MODIS ImageCollection\n",
        "    collection = ee.ImageCollection(product)\\\n",
        "        .filterDate(start_date, end_date)\n",
        "\n",
        "    # Calculate pixel count for each geometry\n",
        "    count_images = []\n",
        "    for geometry in geometries:\n",
        "        count_image = collection.filterBounds(geometry) \\\n",
        "            .reduce(ee.Reducer.count()) \\\n",
        "            .select([0], ['pixel_count']) \\\n",
        "            .clip(geometry)\n",
        "        count_images.append(count_image)\n",
        "\n",
        "    return count_images\n",
        "\n",
        "# Example usage:\n",
        "product = 'MODIS/061/MOD13A2'  # MODIS product path\n",
        "# Define rectangles for different continents (replace with your desired coordinates)\n",
        "continents = [\n",
        "    ee.Geometry.Rectangle([19.85,-3.16, 22.39,-1.31]),  # Africa (approx.)\n",
        "    ee.Geometry.Rectangle([101.23,29.02, 102.10,29.54]),    # Asia (approx.)\n",
        "    ee.Geometry.Rectangle([122.37,-33.08, 122.85,-32.49]),   # Aust_Ocean (approx.)\n",
        "    ee.Geometry.Rectangle([22.45,45.12, 22.97,45.47]),  # Europe (approx.)\n",
        "    ee.Geometry.Rectangle([-121.68,48.02, -119.80,49.01]),    # Nth_America (approx.)\n",
        "    ee.Geometry.Rectangle([-60.70,-9.18, -58.42,-7.56]),   # Sth_America (approx.)\n",
        "    # Add more rectangles for other continents\n",
        "]\n",
        "start_date = '2020-02-01'\n",
        "end_date = '2020-04-30' # Corrected the end date to a valid date\n",
        "\n",
        "pixel_count_images = aggregate_pixel_count(product, continents, start_date, end_date)\n",
        "\n",
        "# Define a function to add Earth Engine layers to a folium map\n",
        "def add_ee_layer(self, ee_image_object, vis_params, name, label): # Added label parameter\n",
        "    map_id_dict = ee.Image(ee_image_object).getMapId(vis_params)\n",
        "    folium.raster_layers.TileLayer(\n",
        "        tiles=map_id_dict['tile_fetcher'].url_format,\n",
        "        attr='Google Earth Engine',\n",
        "        name=name,\n",
        "        overlay=True,\n",
        "        control=True,\n",
        "        show=False # Initially hide the layers\n",
        "    ).add_to(self)\n",
        "\n",
        "    # Add a label for the layer\n",
        "    folium.Marker(\n",
        "        location=ee_image_object.geometry().centroid().coordinates().getInfo()[::-1], # Get the center of the geometry\n",
        "        popup=label,\n",
        "        icon=folium.Icon(color='blue')\n",
        "    ).add_to(self)\n",
        "\n",
        "# Add EE drawing method to folium\n",
        "folium.Map.add_ee_layer = add_ee_layer\n",
        "\n",
        "# Create a folium map (adjust map center as needed)\n",
        "map_center = [20, 20]  # Adjust for a global view\n",
        "my_map = folium.Map(location=map_center, zoom_start=3)\n",
        "\n",
        "# Visualization parameters\n",
        "vis_params = {\n",
        "    'min': 0,\n",
        "    'max': 10,  # Adjust this based on expected range of pixel counts\n",
        "    'palette': ['blue', 'green', 'yellow', 'orange', 'red']\n",
        "}\n",
        "\n",
        "# Labels for each continent\n",
        "continent_labels = [\n",
        "    'Africa', 'Asia', 'Australia/Oceania', 'Europe', 'North America', 'South America'\n",
        "]\n",
        "\n",
        "# Add the pixel count layers to the map with labels\n",
        "for i, count_image in enumerate(pixel_count_images):\n",
        "    my_map.add_ee_layer(count_image, vis_params, f'Pixel Count {i+1}', continent_labels[i]) # Pass the label from continent_labels\n",
        "\n",
        "# Add a layer control panel to the map\n",
        "my_map.add_child(folium.LayerControl())\n",
        "\n",
        "# Display the map\n",
        "my_map  # Or my_map.save(\"pixel_count_map.html\") to save it"
      ],
      "cell_type": "code",
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 674
        },
        "outputId": "47f63fbe-f153-469d-d344-82834850f749",
        "id": "ws73Et5uGLxc"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "<folium.folium.Map at 0x7d17fcadfc10>"
            ],
            "text/html": [
              "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><span style=\"color:#565656\">Make this Notebook Trusted to load map: File -> Trust Notebook</span><iframe srcdoc=\"&lt;!DOCTYPE html&gt;\n",
              "&lt;html&gt;\n",
              "&lt;head&gt;\n",
              "    \n",
              "    &lt;meta http-equiv=&quot;content-type&quot; content=&quot;text/html; charset=UTF-8&quot; /&gt;\n",
              "    \n",
              "        &lt;script&gt;\n",
              "            L_NO_TOUCH = false;\n",
              "            L_DISABLE_3D = false;\n",
              "        &lt;/script&gt;\n",
              "    \n",
              "    &lt;style&gt;html, body {width: 100%;height: 100%;margin: 0;padding: 0;}&lt;/style&gt;\n",
              "    &lt;style&gt;#map {position:absolute;top:0;bottom:0;right:0;left:0;}&lt;/style&gt;\n",
              "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://code.jquery.com/jquery-1.12.4.min.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/js/bootstrap.bundle.min.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/css/bootstrap.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://netdna.bootstrapcdn.com/bootstrap/3.0.0/css/bootstrap.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/@fortawesome/fontawesome-free@6.2.0/css/all.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/gh/python-visualization/folium/folium/templates/leaflet.awesome.rotate.min.css&quot;/&gt;\n",
              "    \n",
              "            &lt;meta name=&quot;viewport&quot; content=&quot;width=device-width,\n",
              "                initial-scale=1.0, maximum-scale=1.0, user-scalable=no&quot; /&gt;\n",
              "            &lt;style&gt;\n",
              "                #map_69d08b70ea0decdd9797229d0bfb624b {\n",
              "                    position: relative;\n",
              "                    width: 100.0%;\n",
              "                    height: 100.0%;\n",
              "                    left: 0.0%;\n",
              "                    top: 0.0%;\n",
              "                }\n",
              "                .leaflet-container { font-size: 1rem; }\n",
              "            &lt;/style&gt;\n",
              "        \n",
              "&lt;/head&gt;\n",
              "&lt;body&gt;\n",
              "    \n",
              "    \n",
              "            &lt;div class=&quot;folium-map&quot; id=&quot;map_69d08b70ea0decdd9797229d0bfb624b&quot; &gt;&lt;/div&gt;\n",
              "        \n",
              "&lt;/body&gt;\n",
              "&lt;script&gt;\n",
              "    \n",
              "    \n",
              "            var map_69d08b70ea0decdd9797229d0bfb624b = L.map(\n",
              "                &quot;map_69d08b70ea0decdd9797229d0bfb624b&quot;,\n",
              "                {\n",
              "                    center: [20.0, 20.0],\n",
              "                    crs: L.CRS.EPSG3857,\n",
              "                    zoom: 3,\n",
              "                    zoomControl: true,\n",
              "                    preferCanvas: false,\n",
              "                }\n",
              "            );\n",
              "\n",
              "            \n",
              "\n",
              "        \n",
              "    \n",
              "            var tile_layer_2374a722904e935627241e89a5e81510 = L.tileLayer(\n",
              "                &quot;https://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Data by \\u0026copy; \\u003ca target=\\&quot;_blank\\&quot; href=\\&quot;http://openstreetmap.org\\&quot;\\u003eOpenStreetMap\\u003c/a\\u003e, under \\u003ca target=\\&quot;_blank\\&quot; href=\\&quot;http://www.openstreetmap.org/copyright\\&quot;\\u003eODbL\\u003c/a\\u003e.&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var tile_layer_b2f2f840c82df8a1b7658dd09c21ecd6 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/316b0b920a12bb241dbbcbdec2c07829-94add5e7720c5e518fc8a26b47b33e05/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_187b9c01a122896ebc703302e29709f7 = L.marker(\n",
              "                [-2.235354019871469, 21.119999999999973],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_3be35ef941368229fe39fc0b28f291c7 = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_187b9c01a122896ebc703302e29709f7.setIcon(icon_3be35ef941368229fe39fc0b28f291c7);\n",
              "        \n",
              "    \n",
              "        var popup_3f383803add628cc1f55cacaf1a9f3d1 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_469532255e4128bf0fe45e9fd1a49ed1 = $(`&lt;div id=&quot;html_469532255e4128bf0fe45e9fd1a49ed1&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Africa&lt;/div&gt;`)[0];\n",
              "                popup_3f383803add628cc1f55cacaf1a9f3d1.setContent(html_469532255e4128bf0fe45e9fd1a49ed1);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_187b9c01a122896ebc703302e29709f7.bindPopup(popup_3f383803add628cc1f55cacaf1a9f3d1)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var tile_layer_193f5e7ca9cb17ca77178769bb92533d = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/005a562f519482986bd8eedd47dbac4f-e9f90ae193d82d64e8ca1068864c3c7b/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_5769fdc1253a610b828b32ab972d55d3 = L.marker(\n",
              "                [29.280483897479627, 101.66499999999978],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_e5dc8127e1c05b10f0e0eb5793806ffe = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_5769fdc1253a610b828b32ab972d55d3.setIcon(icon_e5dc8127e1c05b10f0e0eb5793806ffe);\n",
              "        \n",
              "    \n",
              "        var popup_19fbc22d8c630000f643626fb0ab124e = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_9a79bcf6e5e799f4759211fded558521 = $(`&lt;div id=&quot;html_9a79bcf6e5e799f4759211fded558521&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Asia&lt;/div&gt;`)[0];\n",
              "                popup_19fbc22d8c630000f643626fb0ab124e.setContent(html_9a79bcf6e5e799f4759211fded558521);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_5769fdc1253a610b828b32ab972d55d3.bindPopup(popup_19fbc22d8c630000f643626fb0ab124e)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var tile_layer_258559f0c534a78882dbd6b8f77e269b = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/345dfaca4bd38f08da3d7d44eb2ebc0a-7ead80d5e960ac00747818bd16f62339/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_1b5e7dda867e27af29eb84aff6c699fe = L.marker(\n",
              "                [-32.7849027186664, 122.61000000000037],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_54276d2725d83abd203d121f0a12b43f = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_1b5e7dda867e27af29eb84aff6c699fe.setIcon(icon_54276d2725d83abd203d121f0a12b43f);\n",
              "        \n",
              "    \n",
              "        var popup_b162901fe9faeefeb7f7e07ba1d7e5da = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_574140113edf6905777522015f9b9527 = $(`&lt;div id=&quot;html_574140113edf6905777522015f9b9527&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Australia/Oceania&lt;/div&gt;`)[0];\n",
              "                popup_b162901fe9faeefeb7f7e07ba1d7e5da.setContent(html_574140113edf6905777522015f9b9527);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_1b5e7dda867e27af29eb84aff6c699fe.bindPopup(popup_b162901fe9faeefeb7f7e07ba1d7e5da)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var tile_layer_199fa9e5bebdc9b29bd6e02fed913575 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/13ce548180f932639f88663c4822fd63-3ce10d84e5609784df0b07e2552ac5ef/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_f603ee0fbbee9d00d4ee49e8f855c60c = L.marker(\n",
              "                [45.295114927680466, 22.710000000000402],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_a6b432ca6fe5be9d826a90d0b651c3c2 = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_f603ee0fbbee9d00d4ee49e8f855c60c.setIcon(icon_a6b432ca6fe5be9d826a90d0b651c3c2);\n",
              "        \n",
              "    \n",
              "        var popup_8000889f0b304991ca6634f68c06d0fd = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_84dd324758be8caca7e8c0cb7f7ea19c = $(`&lt;div id=&quot;html_84dd324758be8caca7e8c0cb7f7ea19c&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Europe&lt;/div&gt;`)[0];\n",
              "                popup_8000889f0b304991ca6634f68c06d0fd.setContent(html_84dd324758be8caca7e8c0cb7f7ea19c);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_f603ee0fbbee9d00d4ee49e8f855c60c.bindPopup(popup_8000889f0b304991ca6634f68c06d0fd)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var tile_layer_26fad73cf6c2104eefe1648fa6d9009c = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/d97624f199747268b38ae7f56b11e99c-639cbf94b36befdd274a34d54f08b315/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_f8ade20a1b8ef636bfe0f09f275d5a97 = L.marker(\n",
              "                [48.51721400867622, -120.73999999999994],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_1af20e2c038078735d6bfa2e30c99481 = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_f8ade20a1b8ef636bfe0f09f275d5a97.setIcon(icon_1af20e2c038078735d6bfa2e30c99481);\n",
              "        \n",
              "    \n",
              "        var popup_c34d9def621c7a408f69698cff1f5d35 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_813f6d354fb6872a9f5b1700a8c495a6 = $(`&lt;div id=&quot;html_813f6d354fb6872a9f5b1700a8c495a6&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;North America&lt;/div&gt;`)[0];\n",
              "                popup_c34d9def621c7a408f69698cff1f5d35.setContent(html_813f6d354fb6872a9f5b1700a8c495a6);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_f8ade20a1b8ef636bfe0f09f275d5a97.bindPopup(popup_c34d9def621c7a408f69698cff1f5d35)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var tile_layer_827b458f30bd1bc3ee96822c7fbaaac2 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/2cea209f6deb604adb85c69bfe594257-2e01d8c92410b0340b3624772e559b8f/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_64d262330bfcc9ef9730e9cf42a8d2ab = L.marker(\n",
              "                [-8.371071242702605, -59.55999999999991],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_eecceb7d0e032353bfdfb178d8dafb79 = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_64d262330bfcc9ef9730e9cf42a8d2ab.setIcon(icon_eecceb7d0e032353bfdfb178d8dafb79);\n",
              "        \n",
              "    \n",
              "        var popup_9701b3437217e8fcd80f57489d3547a1 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_305b51c966b7f76eb15e3de83ba5411d = $(`&lt;div id=&quot;html_305b51c966b7f76eb15e3de83ba5411d&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;South America&lt;/div&gt;`)[0];\n",
              "                popup_9701b3437217e8fcd80f57489d3547a1.setContent(html_305b51c966b7f76eb15e3de83ba5411d);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_64d262330bfcc9ef9730e9cf42a8d2ab.bindPopup(popup_9701b3437217e8fcd80f57489d3547a1)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var layer_control_1417db652dc8eacbeb40497b221ef399 = {\n",
              "                base_layers : {\n",
              "                    &quot;openstreetmap&quot; : tile_layer_2374a722904e935627241e89a5e81510,\n",
              "                },\n",
              "                overlays :  {\n",
              "                    &quot;Pixel Count 1&quot; : tile_layer_b2f2f840c82df8a1b7658dd09c21ecd6,\n",
              "                    &quot;Pixel Count 2&quot; : tile_layer_193f5e7ca9cb17ca77178769bb92533d,\n",
              "                    &quot;Pixel Count 3&quot; : tile_layer_258559f0c534a78882dbd6b8f77e269b,\n",
              "                    &quot;Pixel Count 4&quot; : tile_layer_199fa9e5bebdc9b29bd6e02fed913575,\n",
              "                    &quot;Pixel Count 5&quot; : tile_layer_26fad73cf6c2104eefe1648fa6d9009c,\n",
              "                    &quot;Pixel Count 6&quot; : tile_layer_827b458f30bd1bc3ee96822c7fbaaac2,\n",
              "                },\n",
              "            };\n",
              "            L.control.layers(\n",
              "                layer_control_1417db652dc8eacbeb40497b221ef399.base_layers,\n",
              "                layer_control_1417db652dc8eacbeb40497b221ef399.overlays,\n",
              "                {&quot;autoZIndex&quot;: true, &quot;collapsed&quot;: true, &quot;position&quot;: &quot;topright&quot;}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "            tile_layer_b2f2f840c82df8a1b7658dd09c21ecd6.remove();\n",
              "            tile_layer_193f5e7ca9cb17ca77178769bb92533d.remove();\n",
              "            tile_layer_258559f0c534a78882dbd6b8f77e269b.remove();\n",
              "            tile_layer_199fa9e5bebdc9b29bd6e02fed913575.remove();\n",
              "            tile_layer_26fad73cf6c2104eefe1648fa6d9009c.remove();\n",
              "            tile_layer_827b458f30bd1bc3ee96822c7fbaaac2.remove();\n",
              "        \n",
              "&lt;/script&gt;\n",
              "&lt;/html&gt;\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
            ]
          },
          "metadata": {},
          "execution_count": 52
        }
      ]
    },
    {
      "source": [
        "import ee\n",
        "import folium\n",
        "\n",
        "# Initialize Earth Engine\n",
        "# auth.authenticate_user()\n",
        "ee.Initialize(project='ee-pnnwachukwu17')\n",
        "\n",
        "def aggregate_pixel_count(product, geometry, start_date, end_date):\n",
        "    collection = ee.ImageCollection(product) \\\n",
        "        .filterBounds(geometry) \\\n",
        "        .filterDate(start_date, end_date)\n",
        "\n",
        "    count_image = collection.reduce(ee.Reducer.count()).select([0], ['pixel_count'])\n",
        "    count_image = count_image.clip(geometry)\n",
        "    return count_image\n",
        "\n",
        "# Example usage:\n",
        "product = 'MODIS/061/MOD13A2'\n",
        "start_date = '2000-02-01'\n",
        "end_date = '2024-04-30'\n",
        "\n",
        "# Define geometries for each continent (replace with actual coordinates)\n",
        "geometries = {\n",
        "    'Africa': ee.Geometry.Rectangle([19.85,-3.16, 22.39,-1.31]),\n",
        "    'Asia': ee.Geometry.Rectangle([101.23,29.02, 102.10,29.54]),\n",
        "    'Aust_Ocean': ee.Geometry.Rectangle([122.37,-33.08, 122.85,-32.49]),\n",
        "    'Europe': ee.Geometry.Rectangle([22.45,45.12, 22.97,45.47]),\n",
        "    'Nth_America': ee.Geometry.Rectangle([-121.68,48.02, -119.80,49.01]),\n",
        "    'Sth_America': ee.Geometry.Rectangle([-60.70,-9.18, -58.42,-7.56]),\n",
        "    # ... add other continents\n",
        "}\n",
        "\n",
        "pixel_counts = {}\n",
        "for continent, geometry in geometries.items():\n",
        "    pixel_count_image = aggregate_pixel_count(product, geometry, start_date, end_date)\n",
        "    # Get the total pixel count for the continent\n",
        "    pixel_count = pixel_count_image.reduceRegion(\n",
        "        reducer=ee.Reducer.sum(),\n",
        "        geometry=geometry,\n",
        "        scale=500,  #\n",
        "        maxPixels=1e13  # Increase the maximum number of pixels allowed\n",
        "        #bestEffort=True # Or set bestEffort to True to aggregate at a coarser scale\n",
        "    ).get('pixel_count').getInfo()\n",
        "    pixel_counts[continent] = pixel_count\n",
        "\n",
        "# Print the pixel counts for each continent\n",
        "for continent, count in pixel_counts.items():\n",
        "    print(f'Pixel count for {continent}: {count}')"
      ],
      "cell_type": "code",
      "metadata": {
        "id": "-ixr7JLpGLxe"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "pip install earthengine-api matplotlib geopandas folium\n"
      ],
      "metadata": {
        "id": "aIxXGvQeGizY"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "# Import necessary libraries\n",
        "import ee\n",
        "import geopy\n",
        "import folium\n",
        "import pandas as pd\n",
        "import matplotlib.pyplot as plt\n",
        "from datetime import datetime\n",
        "from google.colab import auth\n",
        "from sklearn.metrics.pairwise import rbf_kernel\n",
        "import numpy as np\n",
        "from google.colab import drive\n",
        "import gdown\n",
        "# Import the Polygon class from shapely.geometry\n",
        "from shapely.geometry import Polygon\n",
        "import geopandas as gpd\n"
      ],
      "metadata": {
        "id": "nhwNaVv-GizZ"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "# Initialize Earth Engine\n",
        "auth.authenticate_user()\n",
        "ee.Initialize(project='ee-pnnwachukwu17')"
      ],
      "metadata": {
        "id": "iXVnKrTxGiza"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "source": [
        " # For single continent like Europe"
      ],
      "metadata": {
        "id": "XgtrTzLZGizb"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "# For single continent like Europe\n",
        "\n",
        "# Initialize Earth Engine\n",
        "# auth.authenticate_user()  # Assuming you have already authenticated\n",
        "ee.Initialize(project='ee-pnnwachukwu17')\n",
        "\n",
        "# Function to aggregate pixel count\n",
        "def aggregate_pixel_count(product, geometry, start_date, end_date):\n",
        "    # Load MODIS ImageCollection\n",
        "    collection = ee.ImageCollection(product) \\\n",
        "        .filterBounds(geometry) \\\n",
        "        .filterDate(start_date, end_date)\n",
        "\n",
        "    # Reduce the collection to get the total pixel count in the region of interest\n",
        "    count_image = collection.reduce(ee.Reducer.count()).select([0], ['pixel_count']) # Select the first band and rename it\n",
        "\n",
        "    # Clip to the geometry\n",
        "    count_image = count_image.clip(geometry)\n",
        "\n",
        "    return count_image\n",
        "\n",
        "# Example usage:\n",
        "product = 'MODIS/061/MOD13A2'  # MODIS product path\n",
        "# Define a square geometry (replace with your desired coordinates)\n",
        "geometry = ee.Geometry.Rectangle([\n",
        "    [-122.1, 37.4],  # Bottom left corner\n",
        "    [-122.0, 37.5]   # Top right corner\n",
        "])\n",
        "start_date = '2020-01-01'\n",
        "end_date = '2020-12-31'\n",
        "\n",
        "pixel_count_image = aggregate_pixel_count(product, geometry, start_date, end_date)\n",
        "\n",
        "# Define a function to add Earth Engine layers to a folium map\n",
        "def add_ee_layer(self, ee_image_object, vis_params, name):\n",
        "    map_id_dict = ee.Image(ee_image_object).getMapId(vis_params)\n",
        "    folium.raster_layers.TileLayer(\n",
        "        tiles=map_id_dict['tile_fetcher'].url_format,\n",
        "        attr='Google Earth Engine',\n",
        "        name=name,\n",
        "        overlay=True,\n",
        "        control=True\n",
        "    ).add_to(self)\n",
        "\n",
        "# Add EE drawing method to folium\n",
        "folium.Map.add_ee_layer = add_ee_layer\n",
        "\n",
        "# Create a folium map\n",
        "map_center = geometry.centroid().coordinates().getInfo()[::-1]\n",
        "my_map = folium.Map(location=map_center, zoom_start=12)\n",
        "\n",
        "# Visualization parameters\n",
        "vis_params = {\n",
        "    'min': 0,\n",
        "    'max': 10,  # Adjust this based on expected range of pixel counts\n",
        "    'palette': ['blue', 'green', 'yellow', 'orange', 'red']\n",
        "}\n",
        "\n",
        "# Add the pixel count layer to the map\n",
        "my_map.add_ee_layer(pixel_count_image, vis_params, 'Pixel Count')\n",
        "\n",
        "# Add a layer control panel to the map\n",
        "my_map.add_child(folium.LayerControl())\n",
        "\n",
        "# Display the map\n",
        "my_map  # Or my_map.save(\"pixel_count_map.html\") to save it\n",
        "\n"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "outputId": "2b07f449-1536-494e-d5d3-9947c0e8be9e",
        "id": "dsvtuoNcGizb"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Total number of pixels: None\n"
          ]
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [],
      "metadata": {
        "id": "xo3NaaFaGizd"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "# This is used to print the number of pixels\n",
        "\n",
        "# Function to aggregate pixel count\n",
        "def aggregate_pixel_count(product, geometry, start_date, end_date):\n",
        "    # Load MODIS ImageCollection\n",
        "    collection = ee.ImageCollection(product) \\\n",
        "        .filterBounds(geometry) \\\n",
        "        .filterDate(start_date, end_date)\n",
        "\n",
        "    # Reduce the collection to get the total pixel count in the region of interest\n",
        "    count_image = collection.reduce(ee.Reducer.count()).select([0], ['pixel_count'])\n",
        "\n",
        "    # Clip to the geometry\n",
        "    count_image = count_image.clip(geometry)\n",
        "\n",
        "    return count_image\n",
        "\n",
        "# Example usage:\n",
        "product = 'MODIS/061/MOD13A2'  # MODIS product path\n",
        "\n",
        "# Define a square geometry (replace with your desired coordinates)\n",
        "geometry = ee.Geometry.Rectangle([\n",
        "    22.45,45.12,  # Lower left corner (longitude, latitude)\n",
        "    22.97,45.47   # Upper right corner (longitude, latitude)\n",
        "])\n",
        "\n",
        "start_date = '2020-01-01'\n",
        "end_date = '2020-12-31'\n",
        "\n",
        "pixel_count_image = aggregate_pixel_count(product, geometry, start_date, end_date)\n",
        "\n",
        "# Get the total number of pixels\n",
        "pixel_count = pixel_count_image.reduceRegion(\n",
        "    reducer=ee.Reducer.sum(),\n",
        "    geometry=geometry,\n",
        "    scale=500  # Adjust scale as needed\n",
        ").get('pixel_count').getInfo()\n",
        "\n",
        "print(\"Total number of pixels:\", pixel_count)\n",
        "\n",
        "# ... (rest of the code for creating and displaying the map)"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "outputId": "d9b02b6f-10e9-4c87-9a3e-5b698b851044",
        "id": "6DfKBSS3Gize"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Total number of pixels: 207374.13333333423\n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "source": [
        " # For Multiple Continents"
      ],
      "metadata": {
        "id": "DMAsAq4EGizf"
      }
    },
    {
      "cell_type": "code",
      "source": [],
      "metadata": {
        "id": "IY9QjYmYGizf"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "source": [
        "import ee\n",
        "import folium\n",
        "import matplotlib.pyplot as plt\n",
        "\n",
        "# Initialize Earth Engine\n",
        "# auth.authenticate_user()  # Assuming you have already authenticated\n",
        "ee.Initialize(project='ee-pnnwachukwu17')\n",
        "\n",
        "# Function to aggregate pixel count\n",
        "def aggregate_pixel_count(product, geometries, start_date, end_date):\n",
        "    # Load MODIS ImageCollection\n",
        "    collection = ee.ImageCollection(product)\\\n",
        "        .filterDate(start_date, end_date)\n",
        "\n",
        "    # Calculate pixel count for each geometry\n",
        "    count_images = []\n",
        "    for geometry in geometries:\n",
        "        count_image = collection.filterBounds(geometry) \\\n",
        "            .reduce(ee.Reducer.count()) \\\n",
        "            .select([0], ['pixel_count']) \\\n",
        "            .clip(geometry)\n",
        "        count_images.append(count_image)\n",
        "\n",
        "    return count_images\n",
        "\n",
        "# Example usage:\n",
        "product = 'MODIS/061/MOD13A2'  # MODIS product path\n",
        "# Define rectangles for different continents (replace with your desired coordinates)\n",
        "continents = [\n",
        "    ee.Geometry.Rectangle([19.85,-3.16, 22.39,-1.31]),  # Africa (approx.)\n",
        "    ee.Geometry.Rectangle([101.23,29.02, 102.10,29.54]),    # Asia (approx.)\n",
        "    ee.Geometry.Rectangle([122.37,-33.08, 122.85,-32.49]),   # Aust_Ocean (approx.)\n",
        "    ee.Geometry.Rectangle([22.45,45.12, 22.97,45.47]),  # Europe (approx.)\n",
        "    ee.Geometry.Rectangle([-121.68,48.02, -119.80,49.01]),    # Nth_America (approx.)\n",
        "    ee.Geometry.Rectangle([-60.70,-9.18, -58.42,-7.56]),   # Sth_America (approx.)\n",
        "    # Add more rectangles for other continents\n",
        "]\n",
        "start_date = '2020-01-01'\n",
        "end_date = '2020-12-30'\n",
        "\n",
        "pixel_count_images = aggregate_pixel_count(product, continents, start_date, end_date)\n",
        "\n",
        "# Define a function to add Earth Engine layers to a folium map\n",
        "def add_ee_layer(self, ee_image_object, vis_params, name):\n",
        "    map_id_dict = ee.Image(ee_image_object).getMapId(vis_params)\n",
        "    folium.raster_layers.TileLayer(\n",
        "        tiles=map_id_dict['tile_fetcher'].url_format,\n",
        "        attr='Google Earth Engine',\n",
        "        name=name,\n",
        "        overlay=True,\n",
        "        control=True\n",
        "    ).add_to(self)\n",
        "\n",
        "# Add EE drawing method to folium\n",
        "folium.Map.add_ee_layer = add_ee_layer\n",
        "\n",
        "# Create a folium map (adjust map center as needed)\n",
        "map_center = [20, 20]  # Adjust for a global view\n",
        "my_map = folium.Map(location=map_center, zoom_start=3)\n",
        "\n",
        "# Visualization parameters\n",
        "vis_params = {\n",
        "    'min': 0,\n",
        "    'max': 10,  # Adjust this based on expected range of pixel counts\n",
        "    'palette': ['blue', 'green', 'yellow', 'orange', 'red']\n",
        "}\n",
        "\n",
        "# Add the pixel count layers to the map\n",
        "for i, count_image in enumerate(pixel_count_images):\n",
        "    my_map.add_ee_layer(count_image, vis_params, f'Pixel Count {i+1}')\n",
        "\n",
        "# Add a layer control panel to the map\n",
        "my_map.add_child(folium.LayerControl())\n",
        "\n",
        "# Display the map\n",
        "my_map  # Or my_map.save(\"pixel_count_map.html\") to save it"
      ],
      "cell_type": "code",
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 674
        },
        "outputId": "16b7115c-ca94-4d1c-92ba-efc26e1a01c8",
        "id": "2w0NxQNnGizg"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "<folium.folium.Map at 0x7d17fc5e19c0>"
            ],
            "text/html": [
              "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><span style=\"color:#565656\">Make this Notebook Trusted to load map: File -> Trust Notebook</span><iframe srcdoc=\"&lt;!DOCTYPE html&gt;\n",
              "&lt;html&gt;\n",
              "&lt;head&gt;\n",
              "    \n",
              "    &lt;meta http-equiv=&quot;content-type&quot; content=&quot;text/html; charset=UTF-8&quot; /&gt;\n",
              "    \n",
              "        &lt;script&gt;\n",
              "            L_NO_TOUCH = false;\n",
              "            L_DISABLE_3D = false;\n",
              "        &lt;/script&gt;\n",
              "    \n",
              "    &lt;style&gt;html, body {width: 100%;height: 100%;margin: 0;padding: 0;}&lt;/style&gt;\n",
              "    &lt;style&gt;#map {position:absolute;top:0;bottom:0;right:0;left:0;}&lt;/style&gt;\n",
              "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://code.jquery.com/jquery-1.12.4.min.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/js/bootstrap.bundle.min.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/css/bootstrap.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://netdna.bootstrapcdn.com/bootstrap/3.0.0/css/bootstrap.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/@fortawesome/fontawesome-free@6.2.0/css/all.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/gh/python-visualization/folium/folium/templates/leaflet.awesome.rotate.min.css&quot;/&gt;\n",
              "    \n",
              "            &lt;meta name=&quot;viewport&quot; content=&quot;width=device-width,\n",
              "                initial-scale=1.0, maximum-scale=1.0, user-scalable=no&quot; /&gt;\n",
              "            &lt;style&gt;\n",
              "                #map_88506dba3e37d0c84d220899ed1c3f18 {\n",
              "                    position: relative;\n",
              "                    width: 100.0%;\n",
              "                    height: 100.0%;\n",
              "                    left: 0.0%;\n",
              "                    top: 0.0%;\n",
              "                }\n",
              "                .leaflet-container { font-size: 1rem; }\n",
              "            &lt;/style&gt;\n",
              "        \n",
              "&lt;/head&gt;\n",
              "&lt;body&gt;\n",
              "    \n",
              "    \n",
              "            &lt;div class=&quot;folium-map&quot; id=&quot;map_88506dba3e37d0c84d220899ed1c3f18&quot; &gt;&lt;/div&gt;\n",
              "        \n",
              "&lt;/body&gt;\n",
              "&lt;script&gt;\n",
              "    \n",
              "    \n",
              "            var map_88506dba3e37d0c84d220899ed1c3f18 = L.map(\n",
              "                &quot;map_88506dba3e37d0c84d220899ed1c3f18&quot;,\n",
              "                {\n",
              "                    center: [20.0, 20.0],\n",
              "                    crs: L.CRS.EPSG3857,\n",
              "                    zoom: 3,\n",
              "                    zoomControl: true,\n",
              "                    preferCanvas: false,\n",
              "                }\n",
              "            );\n",
              "\n",
              "            \n",
              "\n",
              "        \n",
              "    \n",
              "            var tile_layer_8c8b232a72bbe131daf16798bdd4bd18 = L.tileLayer(\n",
              "                &quot;https://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Data by \\u0026copy; \\u003ca target=\\&quot;_blank\\&quot; href=\\&quot;http://openstreetmap.org\\&quot;\\u003eOpenStreetMap\\u003c/a\\u003e, under \\u003ca target=\\&quot;_blank\\&quot; href=\\&quot;http://www.openstreetmap.org/copyright\\&quot;\\u003eODbL\\u003c/a\\u003e.&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_52b780d4d9aa48792dd5b9b012abbec7 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/734d5bcea396a99b6eef038ab5343d73-637ec5db63ba17d42ff93bd755d6ec8e/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_328393ad0698fe2ac1b68a308d0e5b76 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/70b5d265a68f02f054886bb35f40b47f-104a74040078397b9cf709d2407563ab/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_6feff248935aa2da82baa46d3001330a = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/3b0be940298b5cd43a7bf8655190b784-ac9b3a3b2ad95469df0dea481fc6998c/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_d6b50773807c73dd7f8c99a0d6ce60c6 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/2ffceb15e0d3aa3d6c2ecb607eca85ab-1488b8a63c5a3862c378912d1e878807/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_3deb4fe552fb1e6b385db989a3b76226 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/81b561aed477752cded11a2cc1746bd9-04c915c55303b8b4e7ae2d8600449553/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var tile_layer_a9b71e91c05830699d0102e256d92a32 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/e9bc0eb17dfbdf3df05a499f02159326-a4301103c274f7176f9547760e7d255b/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "    \n",
              "            var layer_control_ada72f09852909c29e558597ef57af77 = {\n",
              "                base_layers : {\n",
              "                    &quot;openstreetmap&quot; : tile_layer_8c8b232a72bbe131daf16798bdd4bd18,\n",
              "                },\n",
              "                overlays :  {\n",
              "                    &quot;Pixel Count 1&quot; : tile_layer_52b780d4d9aa48792dd5b9b012abbec7,\n",
              "                    &quot;Pixel Count 2&quot; : tile_layer_328393ad0698fe2ac1b68a308d0e5b76,\n",
              "                    &quot;Pixel Count 3&quot; : tile_layer_6feff248935aa2da82baa46d3001330a,\n",
              "                    &quot;Pixel Count 4&quot; : tile_layer_d6b50773807c73dd7f8c99a0d6ce60c6,\n",
              "                    &quot;Pixel Count 5&quot; : tile_layer_3deb4fe552fb1e6b385db989a3b76226,\n",
              "                    &quot;Pixel Count 6&quot; : tile_layer_a9b71e91c05830699d0102e256d92a32,\n",
              "                },\n",
              "            };\n",
              "            L.control.layers(\n",
              "                layer_control_ada72f09852909c29e558597ef57af77.base_layers,\n",
              "                layer_control_ada72f09852909c29e558597ef57af77.overlays,\n",
              "                {&quot;autoZIndex&quot;: true, &quot;collapsed&quot;: true, &quot;position&quot;: &quot;topright&quot;}\n",
              "            ).addTo(map_88506dba3e37d0c84d220899ed1c3f18);\n",
              "        \n",
              "&lt;/script&gt;\n",
              "&lt;/html&gt;\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
            ]
          },
          "metadata": {},
          "execution_count": 51
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [],
      "metadata": {
        "id": "5pdII-FUGizi"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [],
      "metadata": {
        "id": "vKc5afO6Gizj"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "source": [
        "import ee\n",
        "import folium\n",
        "import matplotlib.pyplot as plt\n",
        "\n",
        "# Initialize Earth Engine\n",
        "# auth.authenticate_user()  # Assuming you have already authenticated\n",
        "ee.Initialize(project='ee-pnnwachukwu17')\n",
        "\n",
        "# Function to aggregate pixel count\n",
        "def aggregate_pixel_count(product, geometries, start_date, end_date):\n",
        "    # Load MODIS ImageCollection\n",
        "    collection = ee.ImageCollection(product)\\\n",
        "        .filterDate(start_date, end_date)\n",
        "\n",
        "    # Calculate pixel count for each geometry\n",
        "    count_images = []\n",
        "    for geometry in geometries:\n",
        "        count_image = collection.filterBounds(geometry) \\\n",
        "            .reduce(ee.Reducer.count()) \\\n",
        "            .select([0], ['pixel_count']) \\\n",
        "            .clip(geometry)\n",
        "        count_images.append(count_image)\n",
        "\n",
        "    return count_images\n",
        "\n",
        "# Example usage:\n",
        "product = 'MODIS/061/MOD13A2'  # MODIS product path\n",
        "# Define rectangles for different continents (replace with your desired coordinates)\n",
        "continents = [\n",
        "    ee.Geometry.Rectangle([19.85,-3.16, 22.39,-1.31]),  # Africa (approx.)\n",
        "    ee.Geometry.Rectangle([101.23,29.02, 102.10,29.54]),    # Asia (approx.)\n",
        "    ee.Geometry.Rectangle([122.37,-33.08, 122.85,-32.49]),   # Aust_Ocean (approx.)\n",
        "    ee.Geometry.Rectangle([22.45,45.12, 22.97,45.47]),  # Europe (approx.)\n",
        "    ee.Geometry.Rectangle([-121.68,48.02, -119.80,49.01]),    # Nth_America (approx.)\n",
        "    ee.Geometry.Rectangle([-60.70,-9.18, -58.42,-7.56]),   # Sth_America (approx.)\n",
        "    # Add more rectangles for other continents\n",
        "]\n",
        "start_date = '2020-02-01'\n",
        "end_date = '2020-04-30' # Corrected the end date to a valid date\n",
        "\n",
        "pixel_count_images = aggregate_pixel_count(product, continents, start_date, end_date)\n",
        "\n",
        "# Define a function to add Earth Engine layers to a folium map\n",
        "def add_ee_layer(self, ee_image_object, vis_params, name, label): # Added label parameter\n",
        "    map_id_dict = ee.Image(ee_image_object).getMapId(vis_params)\n",
        "    folium.raster_layers.TileLayer(\n",
        "        tiles=map_id_dict['tile_fetcher'].url_format,\n",
        "        attr='Google Earth Engine',\n",
        "        name=name,\n",
        "        overlay=True,\n",
        "        control=True,\n",
        "        show=False # Initially hide the layers\n",
        "    ).add_to(self)\n",
        "\n",
        "    # Add a label for the layer\n",
        "    folium.Marker(\n",
        "        location=ee_image_object.geometry().centroid().coordinates().getInfo()[::-1], # Get the center of the geometry\n",
        "        popup=label,\n",
        "        icon=folium.Icon(color='blue')\n",
        "    ).add_to(self)\n",
        "\n",
        "# Add EE drawing method to folium\n",
        "folium.Map.add_ee_layer = add_ee_layer\n",
        "\n",
        "# Create a folium map (adjust map center as needed)\n",
        "map_center = [20, 20]  # Adjust for a global view\n",
        "my_map = folium.Map(location=map_center, zoom_start=3)\n",
        "\n",
        "# Visualization parameters\n",
        "vis_params = {\n",
        "    'min': 0,\n",
        "    'max': 10,  # Adjust this based on expected range of pixel counts\n",
        "    'palette': ['blue', 'green', 'yellow', 'orange', 'red']\n",
        "}\n",
        "\n",
        "# Labels for each continent\n",
        "continent_labels = [\n",
        "    'Africa', 'Asia', 'Australia/Oceania', 'Europe', 'North America', 'South America'\n",
        "]\n",
        "\n",
        "# Add the pixel count layers to the map with labels\n",
        "for i, count_image in enumerate(pixel_count_images):\n",
        "    my_map.add_ee_layer(count_image, vis_params, f'Pixel Count {i+1}', continent_labels[i]) # Pass the label from continent_labels\n",
        "\n",
        "# Add a layer control panel to the map\n",
        "my_map.add_child(folium.LayerControl())\n",
        "\n",
        "# Display the map\n",
        "my_map  # Or my_map.save(\"pixel_count_map.html\") to save it"
      ],
      "cell_type": "code",
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 674
        },
        "outputId": "47f63fbe-f153-469d-d344-82834850f749",
        "id": "qdTOlkeTGizj"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "<folium.folium.Map at 0x7d17fcadfc10>"
            ],
            "text/html": [
              "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><span style=\"color:#565656\">Make this Notebook Trusted to load map: File -> Trust Notebook</span><iframe srcdoc=\"&lt;!DOCTYPE html&gt;\n",
              "&lt;html&gt;\n",
              "&lt;head&gt;\n",
              "    \n",
              "    &lt;meta http-equiv=&quot;content-type&quot; content=&quot;text/html; charset=UTF-8&quot; /&gt;\n",
              "    \n",
              "        &lt;script&gt;\n",
              "            L_NO_TOUCH = false;\n",
              "            L_DISABLE_3D = false;\n",
              "        &lt;/script&gt;\n",
              "    \n",
              "    &lt;style&gt;html, body {width: 100%;height: 100%;margin: 0;padding: 0;}&lt;/style&gt;\n",
              "    &lt;style&gt;#map {position:absolute;top:0;bottom:0;right:0;left:0;}&lt;/style&gt;\n",
              "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://code.jquery.com/jquery-1.12.4.min.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/js/bootstrap.bundle.min.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;script src=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.js&quot;&gt;&lt;/script&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/css/bootstrap.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://netdna.bootstrapcdn.com/bootstrap/3.0.0/css/bootstrap.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/@fortawesome/fontawesome-free@6.2.0/css/all.min.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.css&quot;/&gt;\n",
              "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/gh/python-visualization/folium/folium/templates/leaflet.awesome.rotate.min.css&quot;/&gt;\n",
              "    \n",
              "            &lt;meta name=&quot;viewport&quot; content=&quot;width=device-width,\n",
              "                initial-scale=1.0, maximum-scale=1.0, user-scalable=no&quot; /&gt;\n",
              "            &lt;style&gt;\n",
              "                #map_69d08b70ea0decdd9797229d0bfb624b {\n",
              "                    position: relative;\n",
              "                    width: 100.0%;\n",
              "                    height: 100.0%;\n",
              "                    left: 0.0%;\n",
              "                    top: 0.0%;\n",
              "                }\n",
              "                .leaflet-container { font-size: 1rem; }\n",
              "            &lt;/style&gt;\n",
              "        \n",
              "&lt;/head&gt;\n",
              "&lt;body&gt;\n",
              "    \n",
              "    \n",
              "            &lt;div class=&quot;folium-map&quot; id=&quot;map_69d08b70ea0decdd9797229d0bfb624b&quot; &gt;&lt;/div&gt;\n",
              "        \n",
              "&lt;/body&gt;\n",
              "&lt;script&gt;\n",
              "    \n",
              "    \n",
              "            var map_69d08b70ea0decdd9797229d0bfb624b = L.map(\n",
              "                &quot;map_69d08b70ea0decdd9797229d0bfb624b&quot;,\n",
              "                {\n",
              "                    center: [20.0, 20.0],\n",
              "                    crs: L.CRS.EPSG3857,\n",
              "                    zoom: 3,\n",
              "                    zoomControl: true,\n",
              "                    preferCanvas: false,\n",
              "                }\n",
              "            );\n",
              "\n",
              "            \n",
              "\n",
              "        \n",
              "    \n",
              "            var tile_layer_2374a722904e935627241e89a5e81510 = L.tileLayer(\n",
              "                &quot;https://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Data by \\u0026copy; \\u003ca target=\\&quot;_blank\\&quot; href=\\&quot;http://openstreetmap.org\\&quot;\\u003eOpenStreetMap\\u003c/a\\u003e, under \\u003ca target=\\&quot;_blank\\&quot; href=\\&quot;http://www.openstreetmap.org/copyright\\&quot;\\u003eODbL\\u003c/a\\u003e.&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var tile_layer_b2f2f840c82df8a1b7658dd09c21ecd6 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/316b0b920a12bb241dbbcbdec2c07829-94add5e7720c5e518fc8a26b47b33e05/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_187b9c01a122896ebc703302e29709f7 = L.marker(\n",
              "                [-2.235354019871469, 21.119999999999973],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_3be35ef941368229fe39fc0b28f291c7 = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_187b9c01a122896ebc703302e29709f7.setIcon(icon_3be35ef941368229fe39fc0b28f291c7);\n",
              "        \n",
              "    \n",
              "        var popup_3f383803add628cc1f55cacaf1a9f3d1 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_469532255e4128bf0fe45e9fd1a49ed1 = $(`&lt;div id=&quot;html_469532255e4128bf0fe45e9fd1a49ed1&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Africa&lt;/div&gt;`)[0];\n",
              "                popup_3f383803add628cc1f55cacaf1a9f3d1.setContent(html_469532255e4128bf0fe45e9fd1a49ed1);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_187b9c01a122896ebc703302e29709f7.bindPopup(popup_3f383803add628cc1f55cacaf1a9f3d1)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var tile_layer_193f5e7ca9cb17ca77178769bb92533d = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/005a562f519482986bd8eedd47dbac4f-e9f90ae193d82d64e8ca1068864c3c7b/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_5769fdc1253a610b828b32ab972d55d3 = L.marker(\n",
              "                [29.280483897479627, 101.66499999999978],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_e5dc8127e1c05b10f0e0eb5793806ffe = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_5769fdc1253a610b828b32ab972d55d3.setIcon(icon_e5dc8127e1c05b10f0e0eb5793806ffe);\n",
              "        \n",
              "    \n",
              "        var popup_19fbc22d8c630000f643626fb0ab124e = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_9a79bcf6e5e799f4759211fded558521 = $(`&lt;div id=&quot;html_9a79bcf6e5e799f4759211fded558521&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Asia&lt;/div&gt;`)[0];\n",
              "                popup_19fbc22d8c630000f643626fb0ab124e.setContent(html_9a79bcf6e5e799f4759211fded558521);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_5769fdc1253a610b828b32ab972d55d3.bindPopup(popup_19fbc22d8c630000f643626fb0ab124e)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var tile_layer_258559f0c534a78882dbd6b8f77e269b = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/345dfaca4bd38f08da3d7d44eb2ebc0a-7ead80d5e960ac00747818bd16f62339/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_1b5e7dda867e27af29eb84aff6c699fe = L.marker(\n",
              "                [-32.7849027186664, 122.61000000000037],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_54276d2725d83abd203d121f0a12b43f = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_1b5e7dda867e27af29eb84aff6c699fe.setIcon(icon_54276d2725d83abd203d121f0a12b43f);\n",
              "        \n",
              "    \n",
              "        var popup_b162901fe9faeefeb7f7e07ba1d7e5da = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_574140113edf6905777522015f9b9527 = $(`&lt;div id=&quot;html_574140113edf6905777522015f9b9527&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Australia/Oceania&lt;/div&gt;`)[0];\n",
              "                popup_b162901fe9faeefeb7f7e07ba1d7e5da.setContent(html_574140113edf6905777522015f9b9527);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_1b5e7dda867e27af29eb84aff6c699fe.bindPopup(popup_b162901fe9faeefeb7f7e07ba1d7e5da)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var tile_layer_199fa9e5bebdc9b29bd6e02fed913575 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/13ce548180f932639f88663c4822fd63-3ce10d84e5609784df0b07e2552ac5ef/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_f603ee0fbbee9d00d4ee49e8f855c60c = L.marker(\n",
              "                [45.295114927680466, 22.710000000000402],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_a6b432ca6fe5be9d826a90d0b651c3c2 = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_f603ee0fbbee9d00d4ee49e8f855c60c.setIcon(icon_a6b432ca6fe5be9d826a90d0b651c3c2);\n",
              "        \n",
              "    \n",
              "        var popup_8000889f0b304991ca6634f68c06d0fd = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_84dd324758be8caca7e8c0cb7f7ea19c = $(`&lt;div id=&quot;html_84dd324758be8caca7e8c0cb7f7ea19c&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Europe&lt;/div&gt;`)[0];\n",
              "                popup_8000889f0b304991ca6634f68c06d0fd.setContent(html_84dd324758be8caca7e8c0cb7f7ea19c);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_f603ee0fbbee9d00d4ee49e8f855c60c.bindPopup(popup_8000889f0b304991ca6634f68c06d0fd)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var tile_layer_26fad73cf6c2104eefe1648fa6d9009c = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/d97624f199747268b38ae7f56b11e99c-639cbf94b36befdd274a34d54f08b315/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_f8ade20a1b8ef636bfe0f09f275d5a97 = L.marker(\n",
              "                [48.51721400867622, -120.73999999999994],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_1af20e2c038078735d6bfa2e30c99481 = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_f8ade20a1b8ef636bfe0f09f275d5a97.setIcon(icon_1af20e2c038078735d6bfa2e30c99481);\n",
              "        \n",
              "    \n",
              "        var popup_c34d9def621c7a408f69698cff1f5d35 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_813f6d354fb6872a9f5b1700a8c495a6 = $(`&lt;div id=&quot;html_813f6d354fb6872a9f5b1700a8c495a6&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;North America&lt;/div&gt;`)[0];\n",
              "                popup_c34d9def621c7a408f69698cff1f5d35.setContent(html_813f6d354fb6872a9f5b1700a8c495a6);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_f8ade20a1b8ef636bfe0f09f275d5a97.bindPopup(popup_c34d9def621c7a408f69698cff1f5d35)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var tile_layer_827b458f30bd1bc3ee96822c7fbaaac2 = L.tileLayer(\n",
              "                &quot;https://earthengine.googleapis.com/v1/projects/ee-pnnwachukwu17/maps/2cea209f6deb604adb85c69bfe594257-2e01d8c92410b0340b3624772e559b8f/tiles/{z}/{x}/{y}&quot;,\n",
              "                {&quot;attribution&quot;: &quot;Google Earth Engine&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 18, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var marker_64d262330bfcc9ef9730e9cf42a8d2ab = L.marker(\n",
              "                [-8.371071242702605, -59.55999999999991],\n",
              "                {}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "        \n",
              "    \n",
              "            var icon_eecceb7d0e032353bfdfb178d8dafb79 = L.AwesomeMarkers.icon(\n",
              "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;blue&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
              "            );\n",
              "            marker_64d262330bfcc9ef9730e9cf42a8d2ab.setIcon(icon_eecceb7d0e032353bfdfb178d8dafb79);\n",
              "        \n",
              "    \n",
              "        var popup_9701b3437217e8fcd80f57489d3547a1 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
              "\n",
              "        \n",
              "            \n",
              "                var html_305b51c966b7f76eb15e3de83ba5411d = $(`&lt;div id=&quot;html_305b51c966b7f76eb15e3de83ba5411d&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;South America&lt;/div&gt;`)[0];\n",
              "                popup_9701b3437217e8fcd80f57489d3547a1.setContent(html_305b51c966b7f76eb15e3de83ba5411d);\n",
              "            \n",
              "        \n",
              "\n",
              "        marker_64d262330bfcc9ef9730e9cf42a8d2ab.bindPopup(popup_9701b3437217e8fcd80f57489d3547a1)\n",
              "        ;\n",
              "\n",
              "        \n",
              "    \n",
              "    \n",
              "            var layer_control_1417db652dc8eacbeb40497b221ef399 = {\n",
              "                base_layers : {\n",
              "                    &quot;openstreetmap&quot; : tile_layer_2374a722904e935627241e89a5e81510,\n",
              "                },\n",
              "                overlays :  {\n",
              "                    &quot;Pixel Count 1&quot; : tile_layer_b2f2f840c82df8a1b7658dd09c21ecd6,\n",
              "                    &quot;Pixel Count 2&quot; : tile_layer_193f5e7ca9cb17ca77178769bb92533d,\n",
              "                    &quot;Pixel Count 3&quot; : tile_layer_258559f0c534a78882dbd6b8f77e269b,\n",
              "                    &quot;Pixel Count 4&quot; : tile_layer_199fa9e5bebdc9b29bd6e02fed913575,\n",
              "                    &quot;Pixel Count 5&quot; : tile_layer_26fad73cf6c2104eefe1648fa6d9009c,\n",
              "                    &quot;Pixel Count 6&quot; : tile_layer_827b458f30bd1bc3ee96822c7fbaaac2,\n",
              "                },\n",
              "            };\n",
              "            L.control.layers(\n",
              "                layer_control_1417db652dc8eacbeb40497b221ef399.base_layers,\n",
              "                layer_control_1417db652dc8eacbeb40497b221ef399.overlays,\n",
              "                {&quot;autoZIndex&quot;: true, &quot;collapsed&quot;: true, &quot;position&quot;: &quot;topright&quot;}\n",
              "            ).addTo(map_69d08b70ea0decdd9797229d0bfb624b);\n",
              "            tile_layer_b2f2f840c82df8a1b7658dd09c21ecd6.remove();\n",
              "            tile_layer_193f5e7ca9cb17ca77178769bb92533d.remove();\n",
              "            tile_layer_258559f0c534a78882dbd6b8f77e269b.remove();\n",
              "            tile_layer_199fa9e5bebdc9b29bd6e02fed913575.remove();\n",
              "            tile_layer_26fad73cf6c2104eefe1648fa6d9009c.remove();\n",
              "            tile_layer_827b458f30bd1bc3ee96822c7fbaaac2.remove();\n",
              "        \n",
              "&lt;/script&gt;\n",
              "&lt;/html&gt;\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
            ]
          },
          "metadata": {},
          "execution_count": 52
        }
      ]
    },
    {
      "source": [
        "import ee\n",
        "import folium\n",
        "\n",
        "# Initialize Earth Engine\n",
        "# auth.authenticate_user()\n",
        "ee.Initialize(project='ee-pnnwachukwu17')\n",
        "\n",
        "def aggregate_pixel_count(product, geometry, start_date, end_date):\n",
        "    collection = ee.ImageCollection(product) \\\n",
        "        .filterBounds(geometry) \\\n",
        "        .filterDate(start_date, end_date)\n",
        "\n",
        "    count_image = collection.reduce(ee.Reducer.count()).select([0], ['pixel_count'])\n",
        "    count_image = count_image.clip(geometry)\n",
        "    return count_image\n",
        "\n",
        "# Example usage:\n",
        "product = 'MODIS/061/MOD13A2'\n",
        "start_date = '2000-02-01'\n",
        "end_date = '2024-04-30'\n",
        "\n",
        "# Define geometries for each continent (replace with actual coordinates)\n",
        "geometries = {\n",
        "    'Africa': ee.Geometry.Rectangle([19.85,-3.16, 22.39,-1.31]),\n",
        "    'Asia': ee.Geometry.Rectangle([101.23,29.02, 102.10,29.54]),\n",
        "    'Aust_Ocean': ee.Geometry.Rectangle([122.37,-33.08, 122.85,-32.49]),\n",
        "    'Europe': ee.Geometry.Rectangle([22.45,45.12, 22.97,45.47]),\n",
        "    'Nth_America': ee.Geometry.Rectangle([-121.68,48.02, -119.80,49.01]),\n",
        "    'Sth_America': ee.Geometry.Rectangle([-60.70,-9.18, -58.42,-7.56]),\n",
        "    # ... add other continents\n",
        "}\n",
        "\n",
        "pixel_counts = {}\n",
        "for continent, geometry in geometries.items():\n",
        "    pixel_count_image = aggregate_pixel_count(product, geometry, start_date, end_date)\n",
        "    # Get the total pixel count for the continent\n",
        "    pixel_count = pixel_count_image.reduceRegion(\n",
        "        reducer=ee.Reducer.sum(),\n",
        "        geometry=geometry,\n",
        "        scale=500,  # Adjust scale as needed\n",
        "        maxPixels=1e13  # Increase the maximum number of pixels allowed\n",
        "        #bestEffort=True # Or set bestEffort to True to aggregate at a coarser scale\n",
        "    ).get('pixel_count').getInfo()\n",
        "    pixel_counts[continent] = pixel_count\n",
        "\n",
        "# Print the pixel counts for each continent\n",
        "for continent, count in pixel_counts.items():\n",
        "    print(f'Pixel count for {continent}: {count}')"
      ],
      "cell_type": "code",
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "outputId": "d5be6c80-daac-4e57-dc2d-49fda02153b9",
        "id": "WHJiqg8-Gizl"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Pixel count for Africa: 129715142.91764683\n",
            "Pixel count for Asia: 12479566.815686226\n",
            "Pixel count for Aust_Ocean: 7814286.243137279\n",
            "Pixel count for Europe: 5020308.470588274\n",
            "Pixel count for Nth_America: 51338842.070588194\n",
            "Pixel count for Sth_America: 101914020.35686333\n"
          ]
        }
      ]
    }
  ]
}