{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "0",
   "metadata": {},
   "source": [
    "<a href=\"https://githubtocolab.com/gee-community/geemap/blob/master/examples/notebooks/usda_naip_imagery.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open in Colab\"/></a>\n",
    "\n",
    "Uncomment the following line to install [geemap](https://geemap.org) if needed."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "1",
   "metadata": {},
   "outputs": [],
   "source": [
    "# !pip install geemap"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "2",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Installs geemap package\n",
    "import subprocess\n",
    "\n",
    "try:\n",
    "    import geemap\n",
    "except ImportError:\n",
    "    print(\"geemap package not installed. Installing ...\")\n",
    "    subprocess.check_call([\"python\", \"-m\", \"pip\", \"install\", \"geemap\"])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3",
   "metadata": {},
   "outputs": [],
   "source": [
    "import ee"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4",
   "metadata": {},
   "outputs": [],
   "source": [
    "import geemap\n",
    "from ipyleaflet import *\n",
    "from ipywidgets import *\n",
    "from geemap.utils import *\n",
    "\n",
    "Map = geemap.Map(center=(-100, 40), zoom=4)\n",
    "Map.default_style = {\"cursor\": \"pointer\"}\n",
    "Map\n",
    "\n",
    "Map.setOptions(\"ROADMAP\")\n",
    "# Load National Hydrography Dataset (NHD)\n",
    "HUC10 = ee.FeatureCollection(\"USGS/WBD/2017/HUC10\")\n",
    "# 18,487 HUC10 watersheds in the U.S.\n",
    "\n",
    "# Add HUC layer to the map\n",
    "Map.setCenter(-99.00, 47.01, 8)\n",
    "Map.addLayer(\n",
    "    ee.Image().paint(HUC10, 0, 1), {}, \"HUC-10 Watershed\"\n",
    ")  # HUC10 for the entire U.S.\n",
    "\n",
    "label = Label(\"Click on the map to select a watershed\")\n",
    "widget_control = WidgetControl(widget=label, position=\"bottomright\")\n",
    "Map.add_control(widget_control)\n",
    "\n",
    "layer = None\n",
    "\n",
    "\n",
    "def handle_interaction(**kwargs):\n",
    "    latlon = kwargs.get(\"coordinates\")\n",
    "    if kwargs.get(\"type\") == \"click\":\n",
    "        Map.default_style = {\"cursor\": \"wait\"}\n",
    "        xy = ee.Geometry.Point(latlon[::-1])\n",
    "        watershed = HUC10.filterBounds(xy)\n",
    "        huc10_id = watershed.first().get(\"huc10\").getInfo()\n",
    "        Map.layers = Map.layers[:3]\n",
    "        Map.addLayer(\n",
    "            ee.Image().paint(watershed, 0, 2), {\"palette\": \"red\"}, \"HUC ID: \" + huc10_id\n",
    "        )\n",
    "\n",
    "        NAIP_images = find_NAIP(watershed)\n",
    "        first_image = ee.Image(NAIP_images.toList(5).get(0))\n",
    "        Map.addLayer(first_image, {\"bands\": [\"N\", \"R\", \"G\"]}, \"first image\")\n",
    "        count = NAIP_images.size().getInfo()\n",
    "        for i in range(0, count):\n",
    "            image = ee.Image(NAIP_images.toList(count).get(i))\n",
    "            Map.addLayer(image, {\"bands\": [\"N\", \"R\", \"G\"]}, str(i))\n",
    "\n",
    "        Map.default_style = {\"cursor\": \"pointer\"}\n",
    "\n",
    "\n",
    "Map.on_interaction(handle_interaction)\n",
    "\n",
    "Map"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "5",
   "metadata": {},
   "outputs": [],
   "source": [
    "Map.layers"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}