HLC_map_download.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import ee\n",
    "import geemap\n",
    "from ipywidgets import interactive, interactive_output\n",
    "from IPython.display import display\n",
    "\n",
    "Map = geemap.Map()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "tags": []
   },
   "source": [
    "# Hydrological Landscape Classes (HLC) generator app\n",
    "### Application to download HLC map in any region via Google Earth Engine (GEE)\n",
    "#### Developed by Rafael Barbedo Fontana (2021)\n",
    "For more information, access the GitHub page: https://github.com/rbfont/hlc-generator\n",
    "\n",
    "This application has some limitations. If your region of study is too large or if you want more flexibility with the definitions, please use the JavaScript code directly in GEE platform: https://code.earthengine.google.com/f3c135517e2387fcd4f9abc5baf54f72\n",
    "\n",
    "For questions and suggestions, please contact: rbarbedofontana@gmail.com\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "from ipywidgets import widgets, Dropdown, Layout\n",
    "from ipyleaflet import WidgetControl\n",
    "\n",
    "def linked_maps_wlegend(\n",
    "    rows=2,\n",
    "    cols=2,\n",
    "    height=\"400px\",\n",
    "    ee_objects=[],\n",
    "    vis_params=[],\n",
    "    labels=[],\n",
    "    legends=[],\n",
    "    label_position=\"topright\",\n",
    "    **kwargs,\n",
    "                ):\n",
    "    \"\"\"Create linked maps of Earth Engine data layers.\n",
    "\n",
    "    Args:\n",
    "        rows (int, optional): The number of rows of maps to create. Defaults to 2.\n",
    "        cols (int, optional): The number of columns of maps to ceate. Defaults to 2.\n",
    "        height (str, optional): The height of each map in pixels. Defaults to \"400px\".\n",
    "        ee_objects (list, optional): The list of Earth Engine objects to use for each map. Defaults to [].\n",
    "        vis_params (list, optional): The list of visualization parameters to use for each map. Defaults to [].\n",
    "        labels (list, optional): The list of labels to show on the map. Defaults to [].\n",
    "        label_position (str, optional): The position of the label, can be [topleft, topright, bottomleft, bottomright]. Defaults to \"topright\".\n",
    "\n",
    "    Raises:\n",
    "        ValueError: If the length of ee_objects is not equal to rows*cols.\n",
    "        ValueError: If the length of vis_params is not equal to rows*cols.\n",
    "        ValueError: If the length of labels is not equal to rows*cols.\n",
    "\n",
    "    Returns:\n",
    "        ipywidget: A GridspecLayout widget.\n",
    "    \"\"\"\n",
    "    grid = widgets.GridspecLayout(rows, cols, grid_gap=\"0px\")\n",
    "    count = rows * cols\n",
    "\n",
    "    maps = []\n",
    "\n",
    "    if len(ee_objects) > 0:\n",
    "        if len(ee_objects) == 1:\n",
    "            ee_objects = ee_objects * count\n",
    "        elif len(ee_objects) < count:\n",
    "            raise ValueError(f\"The length of ee_objects must be equal to {count}.\")\n",
    "\n",
    "    if len(vis_params) > 0:\n",
    "        if len(vis_params) == 1:\n",
    "            vis_params = vis_params * count\n",
    "        elif len(vis_params) < count:\n",
    "            raise ValueError(f\"The length of vis_params must be equal to {count}.\")\n",
    "\n",
    "    if len(labels) > 0:\n",
    "        if len(labels) == 1:\n",
    "            labels = labels * count\n",
    "        elif len(labels) < count:\n",
    "            raise ValueError(f\"The length of labels must be equal to {count}.\")\n",
    "\n",
    "    for i in range(rows):\n",
    "        for j in range(cols):\n",
    "            index = i * rows + j\n",
    "            m = geemap.Map(\n",
    "                height=height,\n",
    "                lite_mode=True,\n",
    "                add_google_map=False,\n",
    "                layout=widgets.Layout(margin=\"0px\", padding=\"0px\"),\n",
    "                **kwargs,\n",
    "            )\n",
    "\n",
    "            if len(ee_objects) > 0:\n",
    "                m.addLayer(ee_objects[index], vis_params[index], labels[index])\n",
    "                m.add_legend(legend_title=\" \", legend_dict=legends[index])\n",
    "\n",
    "            if len(labels) > 0:\n",
    "                label = widgets.Label(\n",
    "                    labels[index], layout=widgets.Layout(padding=\"0px 5px 0px 5px\")\n",
    "                )\n",
    "                control = WidgetControl(widget=label, position=label_position)\n",
    "                m.add_control(control)\n",
    "\n",
    "            maps.append(m)\n",
    "            widgets.jslink((maps[0], \"center\"), (m, \"center\"))\n",
    "            widgets.jslink((maps[0], \"zoom\"), (m, \"zoom\"))\n",
    "\n",
    "            output = widgets.Output()\n",
    "            with output:\n",
    "                display(m)\n",
    "            grid[i, j] = output\n",
    "\n",
    "    return grid\n",
    "\n",
    "def ee_export_image_zip(\n",
    "    ee_object, filename, scale=None, crs=None, region=None, file_per_band=False\n",
    "):\n",
    "    \"\"\"Exports an ee.Image as a GeoTIFF.\n",
    "\n",
    "    Args:\n",
    "        ee_object (object): The ee.Image to download.\n",
    "        filename (str): Output filename for the exported image.\n",
    "        scale (float, optional): A default scale to use for any bands that do not specify one; ignored if crs and crs_transform is specified. Defaults to None.\n",
    "        crs (str, optional): A default CRS string to use for any bands that do not explicitly specify one. Defaults to None.\n",
    "        region (object, optional): A polygon specifying a region to download; ignored if crs and crs_transform is specified. Defaults to None.\n",
    "        file_per_band (bool, optional): Whether to produce a different GeoTIFF per band. Defaults to False.\n",
    "    \"\"\"\n",
    "    import os\n",
    "    import requests\n",
    "\n",
    "    if not isinstance(ee_object, ee.Image):\n",
    "        print(\"The ee_object must be an ee.Image.\")\n",
    "        return\n",
    "\n",
    "    filename = os.path.abspath(filename)\n",
    "    basename = os.path.basename(filename)\n",
    "    name = os.path.splitext(basename)[0]\n",
    "    filetype = os.path.splitext(basename)[1][1:].lower()\n",
    "    filename_zip = filename.replace(\".tif\", \".zip\")\n",
    "\n",
    "    if filetype != \"tif\":\n",
    "        print(\"The filename must end with .tif\")\n",
    "        return\n",
    "\n",
    "    try:\n",
    "        print(\"Generating URL ...\")\n",
    "        params = {\"name\": name, \"filePerBand\": file_per_band}\n",
    "        if scale is None:\n",
    "            scale = ee_object.projection().nominalScale().multiply(10)\n",
    "        params[\"scale\"] = scale\n",
    "        if region is None:\n",
    "            region = ee_object.geometry()\n",
    "        params[\"region\"] = region\n",
    "        if crs is not None:\n",
    "            params[\"crs\"] = crs\n",
    "\n",
    "        url = ee_object.getDownloadURL(params)\n",
    "        print(\"Click link to download data: {}\".format(url))\n",
    "        r = requests.get(url, stream=True)\n",
    "\n",
    "        if r.status_code != 200:\n",
    "            print(\"An error occurred while downloading.\")\n",
    "            return\n",
    "\n",
    "        with open(filename_zip, \"wb\") as fd:\n",
    "            for chunk in r.iter_content(chunk_size=1024):\n",
    "                fd.write(chunk)\n",
    "\n",
    "    except Exception as e:\n",
    "        print(\"An error occurred while requesting.\")\n",
    "        print(e)\n",
    "        return\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Select map atributes:\n"
     ]
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "f78cb0da2a0f4365a2c1786407d6a239",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "interactive(children=(Dropdown(description='DEM: ', options={'MERIT-DEM 90m': <ee.image.Image object at 0x0000…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import numpy as np\n",
    "'''\n",
    "The app is organized in a nested functions scheme\n",
    "Perhaps not the most elegant way, but it is a lot easier to understand the process (and to perform changes)\n",
    "Funtion 1 (inputs) receives input maps and the scale parameter\n",
    "    and returns the first year available of Land Cover database\n",
    "    and the region of interest (ROI) extension selection\n",
    "Function 2 (roi_from_coords) receives extension of ROI coordinates\n",
    "    and returns widgets for parameters selection\n",
    "Function 3 (params) receives the parameter definition and\n",
    "    performs computations on map visualization and downloading\n",
    "'''\n",
    "\n",
    "def inputs(DEM, HAND, LandCover, Scale):\n",
    "    def roi_from_coords(west, south, east, north):\n",
    "        def params(Slope_th,\n",
    "                   HAND_th,\n",
    "                   LC_year):\n",
    "            \n",
    "            # Image manipulation\n",
    "            slpTh = Slope_th\n",
    "            handTh = HAND_th\n",
    "            lcY = LC_year\n",
    "\n",
    "            tc_img = ee.Image(1).where(\n",
    "                        slp_img.gt(slpTh).And(hand_img.gt(handTh)), 2).where(\n",
    "                        slp_img.lte(slpTh).And(hand_img.gt(handTh)), 3)#.reproject(crs='EPSG:4326', scale=scale)\n",
    "            \n",
    "            if lc_imgcol == mapbiomas:\n",
    "                old_classes = [1, 2, 3, 4, 5, 9, 10, 11, 12, 13, 14, 15, 18, 19, 39, 20, 41, 36, 21, 22, 23, 24, 25, 26, 29, 30, 31, 32, 33]\n",
    "                new_classes = [1, 1, 1, 2, 1, 1,  2,  2,  2,  2,  3,  2,  3,  3,  3,  3,  3,  3,  3,  4,  4,  4,  4,  5,  4,  4,  5,  4,  5]\n",
    "            elif lc_imgcol == igbp:\n",
    "                old_classes = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17]\n",
    "                new_classes = [1, 1, 1, 1, 1, 2, 2, 2, 2,  2,  2,  3,  4,  2,  4,  4,  5]\n",
    "            elif lc_imgcol == umd:\n",
    "                old_classes = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]\n",
    "                new_classes = [5, 1, 1, 1, 1, 1, 2, 2, 2, 2,  2,  2,  3,  4,  2,  4]\n",
    "            elif lc_imgcol == lai:\n",
    "                old_classes = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10]\n",
    "                new_classes = [5, 2, 2, 3, 2, 1, 1, 1, 1, 4,  4]\n",
    "            elif lc_imgcol == pft:\n",
    "                old_classes = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11]\n",
    "                new_classes = [5, 1, 1, 1, 1, 2, 2, 3, 3, 4,  4,  4]\n",
    "            \n",
    "            try:\n",
    "                lcr_img = lc_imgcol.filter(ee.Filter.calendarRange(lcY, lcY, 'year')).first()\n",
    "            except:\n",
    "                lcr_img = lc_imgcol.select('classification_'+ str(lcY))\n",
    "            \n",
    "            lc_img = lcr_img.remap(old_classes, new_classes)\n",
    "            lc_img = ee.Image(5).where(lc_img.lte(5), lc_img)\n",
    "            lc_img = ee.Image(1).where(lc_img.gt(1), lc_img)\n",
    "            \n",
    "            # HAND + Slope + Land cover\n",
    "            #\n",
    "            #  '1': Wetland Forest\n",
    "            #  '2': Wetland non-Forest natural\n",
    "            #  '3': Wetland Farmland\n",
    "            #  '4': Hillslope Forest\n",
    "            #  '5': Hillslope non-Forest natural\n",
    "            #  '6': Hillslope Farmland\n",
    "            #  '7': Plateau Forest\n",
    "            #  '8': Plateau non-Forest natural\n",
    "            #  '9': Plateau Farmland\n",
    "            #  '10': Semi-permeable\n",
    "            #  '11': Water\n",
    "            #\n",
    "            hlc_img = tc_img.multiply(10).addBands(lc_img).reduce(ee.Reducer.sum())\n",
    "            hlc_img = hlc_img.where(\n",
    "                hlc_img.remap([14, 24, 34], [1, 1, 1]), 40).where(\n",
    "                hlc_img.remap([15, 25, 35], [1, 1, 1]), 50)\n",
    "            hlc_img = hlc_img.remap([11, 12, 13, 21, 22, 23, 31, 32, 33, 40, 50],\n",
    "                                    [ 1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11])\n",
    "            hlc_img = ee.Image(1).where(hlc_img.gt(1), hlc_img)#.reproject(crs='EPSG:4326', scale=scale)\n",
    "            \n",
    "            # Render maps -------------------------------------------------------\n",
    "            tc_legend = {\n",
    "                        'Wetland': \"#aec3d4\", \n",
    "                        'Hillslope': \"#cc0013\", \n",
    "                        'Plateau': \"#369b47\"\n",
    "                        }\n",
    "            parvis_tc = {'palette': list(tc_legend.values())}\n",
    "            #Map.addLayer(tc_img, parvis, 'Terrain Classes')\n",
    "\n",
    "            lc_legend = {\n",
    "                        'Forest': \"#152106\",\n",
    "                        'Savanna-Grassland': \"#387242\",\n",
    "                        'Farmland': \"#cdb33b\",\n",
    "                        'Non-vegetated': \"#cc0013\",\n",
    "                        'Water': \"#aec3d4\"\n",
    "                        }        \n",
    "            parvis_lc = {'palette': list(lc_legend.values())}\n",
    "            #Map.addLayer(lc_img, parvis, 'Land Cover Classes')\n",
    "\n",
    "\n",
    "            hlc_legend = {\n",
    "                        \"01 Wetland Forest\": \"#152106\",\n",
    "                        \"02 Wetland Savanna-Grassland\": \"#387242\",\n",
    "                        \"03 Wetland Farmland\": \"#d9903d\",\n",
    "                        \"04 Hillslope Forest\": \"#225129\",\n",
    "                        \"05 Hillslope Savanna-Grassland\": \"#c3aa69\",\n",
    "                        \"06 Hillslope Farmland\": \"#cdb33b\",\n",
    "                        \"07 Plateau Forest\": \"#6a2325\",\n",
    "                        \"08 Plateau Savanna-Grassland\": \"#369b47\",\n",
    "                        \"09 Plateau Farmland\": \"#91af40\",\n",
    "                        \"10 Non-vegetated\": \"#cc0013\",\n",
    "                        \"11 Water\": \"#aec3d4\"\n",
    "                        }\n",
    "            parvis_hlc = {'palette': list(hlc_legend.values())}\n",
    "\n",
    "            imgs = [tc_img, lc_img, hlc_img]\n",
    "            parvis_w = [parvis_tc, parvis_lc, parvis_hlc]\n",
    "            labels = ['Terrain Classes', 'Land Cover Classes', 'Hydrological Landscape Classes']\n",
    "            legends = [tc_legend, lc_legend, hlc_legend]\n",
    "\n",
    "            Maps = linked_maps_wlegend(rows=1, cols=3, height=\"600px\", center=centroid, zoom=7,\n",
    "                                ee_objects=imgs, vis_params=parvis_w, labels=labels, legends=legends, label_position=\"topright\")\n",
    "            \n",
    "            display(Maps)\n",
    "            \n",
    "            # Buttons to download data --------------------------------------------------------\n",
    "            btn_layout = Layout(width='500px')\n",
    "            dem_button = widgets.Button(description='Download DEM map', layout=btn_layout)\n",
    "            slp_button = widgets.Button(description='Download Slope map', layout=btn_layout)\n",
    "            hand_button = widgets.Button(description='Download HAND map', layout=btn_layout)\n",
    "            lcr_button = widgets.Button(description='Download Land Cover map', layout=btn_layout)\n",
    "            buttons_raw = widgets.VBox([dem_button, slp_button, hand_button, lcr_button], layout=Layout(grid_gap='3px'))\n",
    "            \n",
    "            tc_button = widgets.Button(description='Download Terrain Classes map', layout=btn_layout)\n",
    "            lc_button = widgets.Button(description='Download Land Cover Classes map (reclassified)', layout=btn_layout)\n",
    "            hlc_button = widgets.Button(description='Download Hydrological Landscape Classes map', layout=btn_layout)\n",
    "            buttons = widgets.VBox([tc_button, lc_button, hlc_button], layout=Layout(grid_gap='3px'))\n",
    "            \n",
    "            print('\\n\\n\\nDownload raw maps:')\n",
    "            display(buttons_raw)            \n",
    "            \n",
    "            print('\\n\\n\\nDownload processed maps:')\n",
    "            display(buttons)\n",
    "\n",
    "            out = widgets.Output()\n",
    "\n",
    "            def dwn_dem(_):\n",
    "                with out:\n",
    "                    ee_export_image_zip(dem_img, 'dem_map.tif', scale=scale, region=roi)            \n",
    "            def dwn_slp(_):\n",
    "                with out:\n",
    "                    ee_export_image_zip(slp_img, 'slp_map.tif', scale=scale, region=roi)            \n",
    "            def dwn_hand(_):\n",
    "                with out:\n",
    "                    ee_export_image_zip(hand_img, 'hand_map.tif', scale=scale, region=roi)\n",
    "            def dwn_lcr(_):\n",
    "                with out:\n",
    "                    ee_export_image_zip(lcr_img, 'lc_raw_map.tif', scale=scale, region=roi)\n",
    "            \n",
    "            def dwn_tc(_):\n",
    "                with out:\n",
    "                    ee_export_image_zip(tc_img, 'tc_map.tif', scale=scale, region=roi)\n",
    "            def dwn_lc(_):\n",
    "                with out:\n",
    "                    ee_export_image_zip(lc_img, 'lc_map.tif', scale=scale, region=roi)\n",
    "            def dwn_hlc(_):\n",
    "                with out:\n",
    "                    ee_export_image_zip(hlc_img, 'hlc_map.tif', scale=scale, region=roi)\n",
    "\n",
    "            dem_button.on_click(dwn_dem)\n",
    "            slp_button.on_click(dwn_slp)\n",
    "            hand_button.on_click(dwn_hand)\n",
    "            lcr_button.on_click(dwn_lcr)\n",
    "            \n",
    "            tc_button.on_click(dwn_tc)\n",
    "            lc_button.on_click(dwn_lc)\n",
    "            hlc_button.on_click(dwn_hlc)\n",
    "            \n",
    "            display(out)\n",
    "            \n",
    "            return\n",
    "        \n",
    "        roi = ee.Geometry.BBox(float(west), float(south), float(east), float(north))\n",
    "        centroid = list(reversed(roi.centroid(1).getInfo()['coordinates']))    \n",
    "        \n",
    "        # Parameters\n",
    "        lcy_range = widgets.IntSlider(min=start_y, max=2019, value=2010, description='LC classification year: ', style={'description_width': 'initial'})\n",
    "        slope_range = widgets.IntSlider(min=1, max=30, value=10,         description='Slope threshold:        ', style={'description_width': 'initial'})\n",
    "        hand_range = widgets.IntSlider(min=1, max=30, value=10,          description='HAND threshold:         ', style={'description_width': 'initial'})\n",
    "        param_ops = interactive(params, #{'manual': True},\n",
    "                                Slope_th = slope_range,\n",
    "                                HAND_th = hand_range,\n",
    "                                LC_year = lcy_range)\n",
    "        print('\\nSet parameters for map computation:')\n",
    "        display(param_ops)        \n",
    "        \n",
    "        return\n",
    "    \n",
    "    scale = float(Scale)\n",
    "    dem_img = DEM\n",
    "    slp_img = ee.Terrain.slope(dem_img).divide(180).multiply(ee.Number(np.pi)).tan().multiply(100)\n",
    "    hand_img = HAND#.reproject(crs='EPSG:4326', scale=scale)\n",
    "    lc_imgcol = LandCover\n",
    "\n",
    "    if lc_imgcol == mapbiomas:\n",
    "        start_y = 1985\n",
    "    else:\n",
    "        start_y = 2001\n",
    "    \n",
    "    # Region Of Interest    \n",
    "    west = widgets.Text(description='West limit:', value='-55', layout=Layout(width='150px', grid_area='west'))\n",
    "    south = widgets.Text(description='South limit:', value='-15', layout=Layout(width='150px', grid_area='south'))\n",
    "    east = widgets.Text(description='East limit:', value='-51', layout=Layout(width='150px', grid_area='east'))\n",
    "    north = widgets.Text(description='North limit:', value='-11', layout=Layout(width='150px', grid_area='north'));\n",
    "    bbox_coords = widgets.GridBox(children=[north, west, east, south],\n",
    "                layout=Layout(\n",
    "                width='50%',\n",
    "                grid_template_rows='auto auto auto',\n",
    "                grid_template_columns='auto auto auto',\n",
    "                grid_template_areas='''\n",
    "                \". north .\"\n",
    "                \"west . east\"\n",
    "                \". south .\"\n",
    "                '''))    \n",
    "    output = interactive_output(roi_from_coords, {'west': west, 'south': south, 'east': east, 'north': north})\n",
    "    \n",
    "    print('\\nDelimit extension of region to download data:')\n",
    "    display(bbox_coords, output)\n",
    "    \n",
    "    return\n",
    "\n",
    "#folder = widgets.Text(description='Folder', value='D:/')#get_download_folder())\n",
    "\n",
    "# Maps\n",
    "demList = {'MERIT-DEM 90m': ee.Image(\"MERIT/DEM/v1_0_3\").select('dem'),\n",
    "           'SRTM-DEM V4 90m': ee.Image(\"CGIAR/SRTM90_V4\").select('elevation'),\n",
    "           'SRTM-DEM  V3 30m': ee.Image(\"USGS/SRTMGL1_003\").select('elevation'),\n",
    "           'NASA-DEM 30m': ee.Image(\"NASA/NASADEM_HGT/001\").select('elevation')}          \n",
    "handList ={'TPS 90m (South America only)': ee.Image(\"projects/ee-rbfontana/assets/HAND_SA_merit_90m\"),\n",
    "           'FA=8km² 90m': ee.Image(\"users/gena/GlobalHAND/90m-global/hand-1000\"),\n",
    "           'FA=5km² 30m': ee.Image(\"users/gena/GlobalHAND/30m/hand-5000\"),\n",
    "           'FA=1km² 30m': ee.Image(\"users/gena/GlobalHAND/30m/hand-1000\")}\n",
    "lcList = {'Mapbiomas 30m (Brazil only)': ee.Image('projects/mapbiomas-workspace/public/collection5/mapbiomas_collection50_integration_v1'),\n",
    "          'MODIS IGBP 500m': ee.ImageCollection(\"MODIS/006/MCD12Q1\").select('LC_Type1'),\n",
    "          'MODIS UMD 500m': ee.ImageCollection(\"MODIS/006/MCD12Q1\").select('LC_Type2'),\n",
    "          'MODIS LAI 500m': ee.ImageCollection(\"MODIS/006/MCD12Q1\").select('LC_Type3'),\n",
    "          'MODIS PTT 500m': ee.ImageCollection(\"MODIS/006/MCD12Q1\").select('LC_Type5')}\n",
    "\n",
    "mapbiomas = list(lcList.values())[0]\n",
    "igbp = list(lcList.values())[1]\n",
    "umd = list(lcList.values())[2]\n",
    "lai = list(lcList.values())[3]\n",
    "pft = list(lcList.values())[4]\n",
    "\n",
    "\n",
    "demDdwn = widgets.Dropdown(options=demList, description='DEM: ')\n",
    "handDdwn = widgets.Dropdown(options=handList, description='HAND: ')\n",
    "lcDdwn = widgets.Dropdown(options=lcList, description='Land Cover: ')\n",
    "scaleTxt = widgets.Text(description='Scale (m): ', value='90')\n",
    "input_ops = interactive(inputs,\n",
    "                        #folder=folder,\n",
    "                        DEM=demDdwn,\n",
    "                        HAND=handDdwn,\n",
    "                        LandCover=lcDdwn,\n",
    "                        Scale=scaleTxt)\n",
    "\n",
    "print('Select map atributes:')\n",
    "display(input_ops)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}