Add a download module (#1057)

* Add a download module

* [pre-commit.ci] auto fixes from pre-commit.com hooks

for more information, see https://pre-commit.ci

---------

Co-authored-by: pre-commit-ci[bot] <66853113+pre-commit-ci[bot]@users.noreply.github.com>

Author: giswqs

docs/download.md

@@ -0,0 +1,3 @@
+# download module
+
+::: leafmap.download

docs/maplibre/choropleth.ipynb

@@ -126,14 +126,6 @@
    "id": "8",
    "metadata": {},
    "outputs": [],
-   "source": []
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "9",
-   "metadata": {},
-   "outputs": [],
    "source": [
     "m = leafmap.Map(style=\"liberty\")\n",
     "m.add_data(\n",
@@ -152,7 +144,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "10",
+   "id": "9",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -173,7 +165,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "11",
+   "id": "10",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -194,7 +186,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "12",
+   "id": "11",
    "metadata": {},
    "outputs": [],
    "source": [

docs/maplibre/overview.md

@@ -186,7 +186,7 @@ Use the upcase and downcase expressions to change the case of labels.
 
 Create choropleth maps with a variety of classification schemes.
 
-[![](https://github.com/user-attachments/assets/35ca3e17-956e-4d88-9e24-7e94c9e322f8)](https://leafmap.org/maplibre/change_case_of_labels)
+[![](https://github.com/user-attachments/assets/35ca3e17-956e-4d88-9e24-7e94c9e322f8)](https://leafmap.org/maplibre/choropleth)
 
 ## Cloud Optimized GeoTIFF (COG)
 
@@ -212,7 +212,6 @@ Utilize and refine data from the MapTiler Countries to create a Choropleth map o
 
 [![](https://i.imgur.com/k1d6k9I.png)](https://leafmap.org/maplibre/countries_filter)
 
-
 ## Create vector data
 
 Create vector data interactively on a map.

leafmap/common.py

@@ -8689,6 +8689,7 @@ def zonal_stats(
 
     import geopandas as gpd
     import rasterio
+    import xarray as xr
 
     try:
         import rasterstats
@@ -8708,6 +8709,9 @@ def zonal_stats(
             nodata = raster.nodata
             array = raster.read(band_num, masked=True)
             raster_crs = raster.crs
+        elif isinstance(raster, xr.DataArray):
+            array = raster
+            raster_crs = raster.rio.crs
         else:
             array = raster