dea-tools 0.3.7.dev17__tar.gz → 0.3.7.dev19__tar.gz

{dea_tools-0.3.7.dev17 → dea_tools-0.3.7.dev19}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: dea-tools
-Version: 0.3.7.dev17
+Version: 0.3.7.dev19
 Summary: Functions and algorithms for analysing Digital Earth Australia data.
 Project-URL: Homepage, https://github.com/GeoscienceAustralia/dea-notebooks
 Project-URL: Repository, https://github.com/GeoscienceAustralia/dea-notebooks
@@ -22,6 +22,7 @@ Requires-Python: <4.0,>=3.10
 Requires-Dist: aiohttp>=3.5.0
 Requires-Dist: dask-ml>=2023.3.24
 Requires-Dist: dask>=2023.1.0
+Requires-Dist: eo-tides>=0.8.0
 Requires-Dist: fiona>=1.10.0
 Requires-Dist: folium>=0.16.0
 Requires-Dist: geopandas>=0.10.0
@@ -38,7 +39,7 @@ Requires-Dist: owslib>=0.26.0
 Requires-Dist: packaging>=22.0
 Requires-Dist: pandas>=2.2.0
 Requires-Dist: pyproj>=3.7.0
-Requires-Dist: pytmd==2.1.6
+Requires-Dist: python-dateutil>=2.8.2
 Requires-Dist: pytz>=2022.1
 Requires-Dist: rasterio>=1.3.11
 Requires-Dist: rasterstats>=0.16.0
@@ -54,10 +55,11 @@ Requires-Dist: traitlets>=5.0.0
 Requires-Dist: xarray>=2022.3.0
 Provides-Extra: all
 Requires-Dist: affine>=2.3.1; extra == 'all'
-Requires-Dist: bs4>=0.0.1; extra == 'all'
+Requires-Dist: beautifulsoup4>=4.12.0; extra == 'all'
 Requires-Dist: cmocean>=3.0.0; extra == 'all'
 Requires-Dist: contextily>=1.3.0; extra == 'all'
 Requires-Dist: dask-gateway>=2023.1.0; extra == 'all'
+Requires-Dist: datacube-ows>=1.9.3; extra == 'all'
 Requires-Dist: datacube>=1.8.5; extra == 'all'
 Requires-Dist: gcsfs>=2022.1.0; extra == 'all'
 Requires-Dist: hdstats>=0.1.8.post1; extra == 'all'
@@ -73,6 +75,7 @@ Requires-Dist: statsmodels>=0.14.0; extra == 'all'
 Provides-Extra: dask-gateway
 Requires-Dist: dask-gateway>=2023.1.0; extra == 'dask-gateway'
 Provides-Extra: datacube
+Requires-Dist: datacube-ows>=1.9.3; extra == 'datacube'
 Requires-Dist: datacube>=1.8.5; extra == 'datacube'
 Requires-Dist: odc-ui>=0.2.1; extra == 'datacube'
 Provides-Extra: hdstats
@@ -82,10 +85,12 @@ Requires-Dist: ipyleaflet>=0.17.0; extra == 'jupyter'
 Requires-Dist: jupyter>=1.0.0; extra == 'jupyter'
 Provides-Extra: notebooks
 Requires-Dist: affine>=2.3.1; extra == 'notebooks'
-Requires-Dist: bs4>=0.0.1; extra == 'notebooks'
+Requires-Dist: beautifulsoup4>=4.12.0; extra == 'notebooks'
 Requires-Dist: cmocean>=3.0.0; extra == 'notebooks'
 Requires-Dist: contextily>=1.3.0; extra == 'notebooks'
 Requires-Dist: gcsfs>=2022.1.0; extra == 'notebooks'
+Requires-Dist: ipyleaflet>=0.17.0; extra == 'notebooks'
+Requires-Dist: jupyter>=1.0.0; extra == 'notebooks'
 Requires-Dist: opencv-python>=4.6.0.66; extra == 'notebooks'
 Requires-Dist: pydotplus>=2.0.0; extra == 'notebooks'
 Requires-Dist: statsmodels>=0.14.0; extra == 'notebooks'
@@ -95,9 +100,9 @@ Requires-Dist: planetary-computer>=1.0.0; extra == 'stac'
 Requires-Dist: pystac-client>=0.8.3; extra == 'stac'
 Description-Content-Type: text/markdown
 
-# dea-tools Python package
+# DEA Tools Python package
 
-The `dea-tools` Python package provides a collection of functions and algorithms for geospatial analysis using Open Data Cube, Xarray, and Digital Earth Australia.
+`dea_tools`` is an open-source Python package providing functions and algorithms for geospatial analysis using Open Data Cube, Xarray, and Digital Earth Australia.
 It includes utilities for loading data, plotting, spatial and temporal analysis, and applying machine learning to satellite data.
 
 The package is organised into the following modules:
@@ -110,7 +115,7 @@ The package is organised into the following modules:
 * `dea_tools.classification`: Machine learning utilities (e.g. training and applying models on satellite data)
 * `dea_tools.dask`: Utilities for parallel processing with Dask (e.g. creating scalable Dask clusters)
 * `dea_tools.landcover`: Tools for accessing and visualising DEA Land Cover data
-* `dea_tools.coastal`: Coastal and intertidal analysis tools (e.g. tidal tagging, coastal change time series)
+* `dea_tools.coastal`: Coastal and intertidal analysis tools (e.g. coastal change time series, sunglint mapping)
 * `dea_tools.bom`: Accessing Bureau of Meteorology water data (e.g. gauge and discharge data)
 * `dea_tools.waterbodies`: Accessing and analysing DEA Waterbodies data (e.g. loading waterbody time series)
 * `dea_tools.maps`: Tools for interactive mapping (e.g. folium and ipyleaflet maps)

{dea_tools-0.3.7.dev17 → dea_tools-0.3.7.dev19}/README_tools.md

@@ -1,6 +1,6 @@
-# dea-tools Python package
+# DEA Tools Python package
 
-The `dea-tools` Python package provides a collection of functions and algorithms for geospatial analysis using Open Data Cube, Xarray, and Digital Earth Australia.
+`dea_tools`` is an open-source Python package providing functions and algorithms for geospatial analysis using Open Data Cube, Xarray, and Digital Earth Australia.
 It includes utilities for loading data, plotting, spatial and temporal analysis, and applying machine learning to satellite data.
 
 The package is organised into the following modules:
@@ -13,7 +13,7 @@ The package is organised into the following modules:
 * `dea_tools.classification`: Machine learning utilities (e.g. training and applying models on satellite data)
 * `dea_tools.dask`: Utilities for parallel processing with Dask (e.g. creating scalable Dask clusters)
 * `dea_tools.landcover`: Tools for accessing and visualising DEA Land Cover data
-* `dea_tools.coastal`: Coastal and intertidal analysis tools (e.g. tidal tagging, coastal change time series)
+* `dea_tools.coastal`: Coastal and intertidal analysis tools (e.g. coastal change time series, sunglint mapping)
 * `dea_tools.bom`: Accessing Bureau of Meteorology water data (e.g. gauge and discharge data)
 * `dea_tools.waterbodies`: Accessing and analysing DEA Waterbodies data (e.g. loading waterbody time series)
 * `dea_tools.maps`: Tools for interactive mapping (e.g. folium and ipyleaflet maps)

{dea_tools-0.3.7.dev17/dea_tools → dea_tools-0.3.7.dev19}/__init__.py

@@ -1,4 +1,4 @@
 try:
     from .__version__ import version as __version__
 except ImportError:
-    __version__ = None
+    __version__ = None

{dea_tools-0.3.7.dev17/dea_tools → dea_tools-0.3.7.dev19}/app/animations.py

@@ -1,57 +1,55 @@
 # -*- coding: utf-8 -*-
 """
-Satellite imagery animation widget, which can be used to interactively 
+Satellite imagery animation widget, which can be used to interactively
 produce animations for multiple DEA products.
 """
 
 # Import required packages
-import fiona
-import sys
-import datacube
+import datetime
+import itertools
+import json
 import warnings
+from io import BytesIO
+
+import datacube
+import geopandas as gpd
+import ipywidgets as widgets
+import matplotlib as mpl
 import matplotlib.pyplot as plt
-from datacube.utils.geometry import CRS
+import numpy as np
+import pandas as pd
+from datacube.utils.geometry import Geometry
+from datacube.utils.masking import mask_invalid_data
 from ipyleaflet import (
-    WMSLayer,
-    basemaps,
-    basemap_to_tiles,
-    Map,
-    DrawControl,
-    WidgetControl,
-    SearchControl,
-    Marker,
     LayerGroup,
-    LayersControl,
-    GeoData,
+    Marker,
+    SearchControl,
+    basemap_to_tiles,
+    basemaps,
 )
-from traitlets import Unicode
 from ipywidgets import (
-    GridspecLayout,
+    HTML,
     Button,
-    Layout,
+    GridspecLayout,
     HBox,
-    VBox,
-    HTML,
+    Layout,
     Output,
+    VBox,
 )
-import json
-import itertools
-import numpy as np
-import geopandas as gpd
-from io import BytesIO
-import ipywidgets as widgets
-import datetime
-from skimage import exposure
 from skimage.filters import unsharp_mask
 
-from datacube.utils import masking
-from datacube.utils.geometry import Geometry
-from datacube.utils.masking import mask_invalid_data
-import dea_tools.app.widgetconstructors as deawidgets
+from dea_tools.app.widgetconstructors import (
+    create_checkbox,
+    create_datepicker,
+    create_drawcontrol,
+    create_dropdown,
+    create_html,
+    create_map,
+)
+from dea_tools.coastal import get_coastlines
 from dea_tools.dask import create_local_dask_cluster
 from dea_tools.spatial import reverse_geocode
 
-import warnings
 warnings.filterwarnings("ignore")
 
 # WMS params and satellite style bands
@@ -137,14 +135,12 @@ def update_map_layers(self, update_basemap=False):
     self.query_params = None
 
     if update_basemap:
-
         # Clear all layers and add basemap
         self.map_layers.clear_layers()
         self.map_layers.add_layer(self.basemap)
 
 
 def extract_data(self):
-
     # Connect to datacube database
     dc = datacube.Datacube(app="Exporting satellite images")
 
@@ -166,15 +162,11 @@ def extract_data(self):
     }
 
     # Find matching datasets
-    dss = [
-        dc.find_datasets(product=i, **self.query_params)
-        for i in sat_params[self.dealayer]["products"]
-    ]
+    dss = [dc.find_datasets(product=i, **self.query_params) for i in sat_params[self.dealayer]["products"]]
     dss = list(itertools.chain.from_iterable(dss))
 
     # If data is found
     if len(dss) > 0:
-
         # Get CRS
         crs = str(dss[0].crs)
 
@@ -224,26 +216,21 @@ def extract_data(self):
 
 
 def plot_data(self, fname):
-
     # Data to plot
     to_plot = self.timeseries_ds
 
     # If rolling median specified
     if self.rolling_median:
         with self.status_info:
-            print(
-                f"\nApplying rolling median ({self.rolling_median_window} timesteps window)"
-            )
-        to_plot = to_plot.rolling(
-            time=int(self.rolling_median_window), center=True, min_periods=1
-        ).median()
+            print(f"\nApplying rolling median ({self.rolling_median_window} timesteps window)")
+        to_plot = to_plot.rolling(time=int(self.rolling_median_window), center=True, min_periods=1).median()
 
     # Raise by power to dampen bright features and enhance dark.
     # Raise vmin and vmax by same amount to ensure proper stretch
     if self.power < 1.0:
         with self.status_info:
             print(f"\nApplying power transformation ({self.power})")
-        to_plot = to_plot ** self.power
+        to_plot = to_plot**self.power
 
     # Apply unsharp masking to enhance overall dynamic range,
     # and improve fine scale detail
@@ -257,9 +244,7 @@ def plot_data(self, fname):
 
         funcs_list = [
             rescale_intensity,
-            lambda x: unsharp_mask(
-                x, radius=self.unsharp_mask_radius, amount=self.unsharp_mask_amount
-            ),
+            lambda x: unsharp_mask(x, radius=self.unsharp_mask_radius, amount=self.unsharp_mask_amount),
         ]
     else:
         funcs_list = None
@@ -288,13 +273,6 @@ def plot_data(self, fname):
 
 
 def deacoastlines_overlay(ds):
-
-    import geopandas as gpd
-    import pandas as pd
-    import matplotlib
-    from shapely.geometry import box, Point
-    from dea_tools.coastal import get_coastlines
-
     # Get bounding box of data
     xmin, ymin, xmax, ymax = ds.geobox.geographic_extent.boundingbox
     bounds = [xmin, ymin, xmax, ymax]
@@ -308,8 +286,8 @@ def deacoastlines_overlay(ds):
     deacl_gdf = deacl_gdf.dissolve("year")  # values("year", ascending=True)
 
     # Apply colours
-    norm = matplotlib.colors.Normalize(vmin=0, vmax=len(deacl_gdf.index))
-    cmap = matplotlib.cm.get_cmap("inferno")
+    norm = mpl.colors.Normalize(vmin=0, vmax=len(deacl_gdf.index))
+    cmap = mpl.cm.get_cmap("inferno")
     rgba = cmap(norm(deacl_gdf.reset_index().index))
     deacl_gdf["color"] = list(rgba)
     deacl_gdf["start_time"] = pd.to_datetime(deacl_gdf.index) + pd.DateOffset(months=0)
@@ -317,8 +295,7 @@ def deacoastlines_overlay(ds):
 
     if len(deacl_gdf.index) > 0:
         return deacl_gdf
-    else:
-        return None
+    return None
 
 
 class animation_app(HBox):
@@ -338,15 +315,13 @@ class animation_app(HBox):
 
         # Satellite data
         end_date = datetime.datetime.today()
-        start_date = datetime.datetime(
-            year=end_date.year - 3, month=end_date.month, day=end_date.day
-        )
+        start_date = datetime.datetime(year=end_date.year - 3, month=end_date.month, day=end_date.day)
         self.start_date = start_date.strftime("%Y-%m-%d")
         self.end_date = end_date.strftime("%Y-%m-%d")
         self.dealayer_list = [
             ("Landsat", "Landsat"),
             ("Sentinel-2", "Sentinel-2"),
-            ("Sentinel-2 and Landsat", "Sentinel-2 and Landsat")
+            ("Sentinel-2 and Landsat", "Sentinel-2 and Landsat"),
         ]
         self.dealayer = self.dealayer_list[0][1]
 
@@ -395,16 +370,14 @@ class animation_app(HBox):
         ##################
 
         # Create the Header widget
-        header_title_text = (
-            "<h3>Digital Earth Australia satellite imagery animations</h3>"
-        )
+        header_title_text = "<h3>Digital Earth Australia satellite imagery animations</h3>"
         instruction_text = (
             "<p>Select the desired satellite data, imagery date range "
             "and image style, then zoom in and draw a rectangle to "
             "select an area export as a satellite imagery time-series "
             "animation.</p>"
         )
-        self.header = deawidgets.create_html(f"{header_title_text}{instruction_text}")
+        self.header = create_html(f"{header_title_text}{instruction_text}")
         self.header.layout = make_box_layout()
 
         #####################################
@@ -413,7 +386,6 @@ class animation_app(HBox):
 
         # Define the action to take once something is drawn on the map
         def update_geojson(target, action, geo_json):
-
             # Get data from action
             self.action = action
 
@@ -443,12 +415,7 @@ class animation_app(HBox):
                     '<span style="color: #33cc33"> '
                     "<b>(Overriding maximum size limit; use with caution as may lead to memory issues)</b></span>"
                 )
-                self.header.value = (
-                    header_title_text
-                    + instruction_text
-                    + polyarea_text
-                    + confirmation_text
-                )
+                self.header.value = header_title_text + instruction_text + polyarea_text + confirmation_text
                 self.gdf_drawn = gdf
             elif area <= 50000:
                 confirmation_text = (
@@ -456,12 +423,7 @@ class animation_app(HBox):
                     "<b>(Area to extract falls within "
                     "recommended 50000 ha limit)</b></span>"
                 )
-                self.header.value = (
-                    header_title_text
-                    + instruction_text
-                    + polyarea_text
-                    + confirmation_text
-                )
+                self.header.value = header_title_text + instruction_text + polyarea_text + confirmation_text
                 self.gdf_drawn = gdf
             else:
                 warning_text = (
@@ -470,9 +432,7 @@ class animation_app(HBox):
                     "please select an area less than 50000 "
                     "ha)</b></span>"
                 )
-                self.header.value = (
-                    header_title_text + instruction_text + polyarea_text + warning_text
-                )
+                self.header.value = header_title_text + instruction_text + polyarea_text + warning_text
                 self.gdf_drawn = None
 
         ###########################
@@ -488,24 +448,26 @@ class animation_app(HBox):
 
         # Create drawing tools
         desired_drawtools = ["rectangle"]
-        draw_control = deawidgets.create_drawcontrol(desired_drawtools)
+        draw_control = create_drawcontrol(desired_drawtools)
 
         # Begin by displaying an empty layer group, and update the group with desired WMS on interaction.
         self.map_layers = LayerGroup(layers=())
         self.map_layers.name = "Map Overlays"
 
         # Create map widget
-        self.m = deawidgets.create_map(map_center=(-33.96, 151.20), zoom_level=13)
+        self.m = create_map(map_center=(-33.96, 151.20), zoom_level=13)
         self.m.layout = make_box_layout()
 
         # Add tools to map widget
         self.m.add_control(draw_control)
-        self.m.add_control(SearchControl(
-        position="topleft",
-        url='https://nominatim.openstreetmap.org/search?format=json&q={s}',
-        zoom=13, # 'Village / Suburb' level zoom
-        marker=Marker(draggable=False)
-        ))
+        self.m.add_control(
+            SearchControl(
+                position="topleft",
+                url="https://nominatim.openstreetmap.org/search?format=json&q={s}",
+                zoom=13,  # 'Village / Suburb' level zoom
+                marker=Marker(draggable=False),
+            )
+        )
         self.m.add_layer(self.map_layers)
 
         # Update all maps to starting defaults
@@ -516,24 +478,16 @@ class animation_app(HBox):
         ############################
 
         # Create parameter widgets
-        dropdown_basemap = deawidgets.create_dropdown(
-            self.basemap_list, self.basemap_list[0][1]
-        )
-        dropdown_dealayer = deawidgets.create_dropdown(
-            self.dealayer_list, self.dealayer_list[0][1]
-        )
-        dropdown_output = deawidgets.create_dropdown(
-            self.output_list, self.output_list[0][1]
-        )
-        date_picker_start = deawidgets.create_datepicker(
+        dropdown_basemap = create_dropdown(self.basemap_list, self.basemap_list[0][1])
+        dropdown_dealayer = create_dropdown(self.dealayer_list, self.dealayer_list[0][1])
+        dropdown_output = create_dropdown(self.output_list, self.output_list[0][1])
+        date_picker_start = create_datepicker(
             value=start_date,
         )
-        date_picker_end = deawidgets.create_datepicker(
+        date_picker_end = create_datepicker(
             value=end_date,
         )
-        dropdown_styles = deawidgets.create_dropdown(
-            self.styles_list, self.styles_list[0]
-        )
+        dropdown_styles = create_dropdown(self.styles_list, self.styles_list[0])
         slider_percentile = widgets.FloatRangeSlider(
             value=[0.01, 0.99],
             min=0,
@@ -553,7 +507,7 @@ class animation_app(HBox):
             layout={"width": "85%"},
         )
 
-        checkbox_rolling_median = deawidgets.create_checkbox(
+        checkbox_rolling_median = create_checkbox(
             self.rolling_median,
             "Apply rolling median to produce smooth, cloud-free animations",
             layout={"width": "85%"},
@@ -571,26 +525,20 @@ class animation_app(HBox):
         )
 
         # Expandable advanced section
-        text_interval = widgets.IntText(
-            value=100, description="", step=50, layout={"width": "95%"}
-        )
+        text_interval = widgets.IntText(value=100, description="", step=50, layout={"width": "95%"})
         text_resolution = widgets.FloatText(
             value=30,
             description="",
             layout={"width": "95%", "margin": "0px", "padding": "0px"},
         )
-        text_width = widgets.IntText(
-            value=900, description="", step=50, layout={"width": "95%"}
-        )
-        dropdown_resampling = deawidgets.create_dropdown(
+        text_width = widgets.IntText(value=900, description="", step=50, layout={"width": "95%"})
+        dropdown_resampling = create_dropdown(
             self.resample_list,
             self.resample_freq,
             description="",
             layout={"width": "95%"},
         )
-        checkbox_cloud_mask = deawidgets.create_checkbox(
-            self.cloud_mask, "Mask out cloudy pixels", layout={"width": "95%"}
-        )
+        checkbox_cloud_mask = create_checkbox(self.cloud_mask, "Mask out cloudy pixels", layout={"width": "95%"})
         slider_power = widgets.FloatSlider(
             value=1.0,
             min=0.01,
@@ -599,9 +547,7 @@ class animation_app(HBox):
             description="",
             layout={"width": "95%"},
         )
-        checkbox_unsharp_mask = deawidgets.create_checkbox(
-            self.unsharp_mask, "Enable", layout={"width": "95%"}
-        )
+        checkbox_unsharp_mask = create_checkbox(self.unsharp_mask, "Enable", layout={"width": "95%"})
         text_unsharp_mask_radius = widgets.FloatText(
             value=20,
             step=1,
@@ -624,12 +570,10 @@ class animation_app(HBox):
                 "display": "none",
             },
         )
-        checkbox_deacoastlines = deawidgets.create_checkbox(
+        checkbox_deacoastlines = create_checkbox(
             self.deacoastlines, "Add DEA Coastlines overlay", layout={"width": "95%"}
         )
-        checkbox_max_size = deawidgets.create_checkbox(
-            self.max_size, "Enable", layout={"width": "95%"}
-        )
+        checkbox_max_size = create_checkbox(self.max_size, "Enable", layout={"width": "95%"})
         expand_box = widgets.VBox(
             [
                 HTML("Frame interval (milliseconds):"),
@@ -649,9 +593,7 @@ class animation_app(HBox):
                 checkbox_unsharp_mask,
                 text_unsharp_mask_radius,
                 text_unsharp_mask_amount,
-                HTML(
-                    "</br>Override maximum size limit: (use with caution; may cause memory issues/crashes)"
-                ),
+                HTML("</br>Override maximum size limit: (use with caution; may cause memory issues/crashes)"),
                 checkbox_max_size,
             ],
         )
@@ -680,13 +622,9 @@ class animation_app(HBox):
         dropdown_styles.observe(self.update_styles, "value")
 
         slider_percentile.observe(self.update_slider_percentile, "value")
-        floatslider_max_cloud_cover.observe(
-            self.update_floatslider_max_cloud_cover, "value"
-        )
+        floatslider_max_cloud_cover.observe(self.update_floatslider_max_cloud_cover, "value")
         checkbox_rolling_median.observe(self.update_checkbox_rolling_median, "value")
-        text_rolling_median_window.observe(
-            self.update_text_rolling_median_window, "value"
-        )
+        text_rolling_median_window.observe(self.update_text_rolling_median_window, "value")
         dropdown_output.observe(self.update_output, "value")
         run_button.on_click(self.run_app)
         draw_control.on_draw(update_geojson)
@@ -708,34 +646,30 @@ class animation_app(HBox):
         # COLLECTION OF ALL APP CONTROLS #
         ##################################
 
-        parameter_selection = VBox(
-            [
-                HTML("<b>Start date:</b>"),
-                date_picker_start,
-                HTML("<b>End date:</b>"),
-                date_picker_end,
-                HTML("<b>Satellite imagery:</b>"),
-                dropdown_dealayer,
-                HTML("<b>Style:</b>"),
-                dropdown_styles,
-                HTML("<b>Colour percentile stretch:</b>"),
-                slider_percentile,
-                HTML("<b>Maximum cloud cover (%):</b>"),
-                floatslider_max_cloud_cover,
-                checkbox_rolling_median,
-                text_rolling_median_window,
-                HTML("</br><b>Output file format:</b>"),
-                dropdown_output,
-                HTML("</br>"),
-                expand,
-            ]
-        )
-        map_selection = VBox(
-            [
-                HTML("</br><b>Map overlay:</b>"),
-                dropdown_basemap,
-            ]
-        )
+        parameter_selection = VBox([
+            HTML("<b>Start date:</b>"),
+            date_picker_start,
+            HTML("<b>End date:</b>"),
+            date_picker_end,
+            HTML("<b>Satellite imagery:</b>"),
+            dropdown_dealayer,
+            HTML("<b>Style:</b>"),
+            dropdown_styles,
+            HTML("<b>Colour percentile stretch:</b>"),
+            slider_percentile,
+            HTML("<b>Maximum cloud cover (%):</b>"),
+            floatslider_max_cloud_cover,
+            checkbox_rolling_median,
+            text_rolling_median_window,
+            HTML("</br><b>Output file format:</b>"),
+            dropdown_output,
+            HTML("</br>"),
+            expand,
+        ])
+        map_selection = VBox([
+            HTML("</br><b>Map overlay:</b>"),
+            dropdown_basemap,
+        ])
         parameter_selection.layout = make_box_layout()
         map_selection.layout = make_box_layout()
 
@@ -810,10 +744,10 @@ class animation_app(HBox):
 
         elif change.new == "Sentinel-2":
             self.text_resolution.value = 10
-        
+
         elif change.new == "Sentinel-2 and Landsat":
             self.text_resolution.value = 30
-            
+
         # Clear data load params to trigger data re-load
         update_map_layers(self)
 
@@ -915,16 +849,13 @@ class animation_app(HBox):
         update_map_layers(self)
 
     def run_app(self, change):
-
         # Clear progress bar and output areas before running
         self.status_info.clear_output()
         self.output_plot.clear_output()
 
         # Verify that polygon was drawn
         if self.gdf_drawn is not None:
-
             with self.status_info:
-
                 # Load data and add to attribute
                 if self.timeseries_ds is None:
                     self.timeseries_ds = extract_data(self)
@@ -933,9 +864,7 @@ class animation_app(HBox):
                     print("Using previously loaded data")
 
             if self.timeseries_ds is not None:
-
                 with self.status_info:
-
                     # Create unique file name
                     centre_coords = self.gdf_drawn.geometry[0].centroid.coords[0][::-1]
                     site = reverse_geocode(coords=centre_coords)
@@ -947,9 +876,7 @@ class animation_app(HBox):
                         .lower()
                     )
 
-                    print(
-                        f"\nExporting animation for {site}.\nThis may take several minutes..."
-                    )
+                    print(f"\nExporting animation for {site}.\nThis may take several minutes...")
 
                 ############
                 # Plotting #
@@ -968,6 +895,4 @@ class animation_app(HBox):
 
         else:
             with self.status_info:
-                print(
-                    'Please draw a valid rectangle on the map, then press "Generate animation".'
-                )
+                print('Please draw a valid rectangle on the map, then press "Generate animation".')