Skip to content

Instantly share code, notes, and snippets.

@soxofaan
Created October 2, 2023 14:51
Show Gist options
  • Save soxofaan/aaf001c8ee6aca4198d1475b9ca8741e to your computer and use it in GitHub Desktop.
Save soxofaan/aaf001c8ee6aca4198d1475b9ca8741e to your computer and use it in GitHub Desktop.
Open-EO/openeo-geopyspark-driver#527 Preserve nan in comparisons
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"id": "0fd0c8cc-dfb7-487c-a7ce-cf2c1e1e55f8",
"metadata": {},
"source": [
"https://discuss.eodc.eu/t/boolean-mask-and-mask-polygon-issue/626/4"
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "e4866fab-914b-433f-9251-67f948e490dc",
"metadata": {},
"outputs": [],
"source": [
"import openeo\n",
"import shapely.geometry\n",
"import xarray\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "c512637e-b1e0-4fbb-adf0-fc2d2a880e85",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Authenticated using refresh token.\n"
]
},
{
"data": {
"text/plain": [
"<Connection to 'https://openeo.vito.be/openeo/1.1/' with OidcBearerAuth>"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"connection = openeo.connect(\n",
" \"openeo.vito.be\"\n",
")\n",
"connection.authenticate_oidc()"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "e60d5247-b111-4edf-8a13-ab0cac0c3f97",
"metadata": {},
"outputs": [
{
"data": {
"image/svg+xml": [
"<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"100.0\" height=\"100.0\" viewBox=\"5.084 51.184000000000005 0.4320000000000004 0.1319999999999908\" preserveAspectRatio=\"xMinYMin meet\"><g transform=\"matrix(1,0,0,-1,0,102.5)\"><path fill-rule=\"evenodd\" fill=\"#66cc99\" stroke=\"#555555\" stroke-width=\"0.008640000000000007\" opacity=\"0.6\" d=\"M 5.1,51.3 L 5.1,51.2 L 5.5,51.2 L 5.14,51.22 L 5.1,51.3 z\" /></g></svg>"
],
"text/plain": [
"<shapely.geometry.polygon.Polygon at 0x7f9076204a30>"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"geometry = shapely.geometry.shape(\n",
" # {\"type\": \"Polygon\", \"coordinates\": [[[5.1, 51.23], [5.08, 51.21], [5.13, 51.2], [5.1, 51.23]]]}\n",
" {\"type\": \"Polygon\", \"coordinates\": [[[5.1, 51.3], [5.1, 51.2], [5.5, 51.2], [5.14, 51.22], [5.1, 51.3]]]}\n",
")\n",
"geometry"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "3c38c4c9-64b0-4677-8e74-21de7d90c3e5",
"metadata": {},
"outputs": [],
"source": [
"collection = \"SENTINEL2_L2A\"\n",
"# spatial_extent = dict(zip([\"west\", \"south\", \"east\", \"north\"], geometry.bounds))\n",
"temporal_extent = [\"2018-09-15\", \"2018-09-20\"]\n",
"bands = ['B03', 'B11']"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "f90bbbed-a467-44c4-9588-11b97d39d1ca",
"metadata": {},
"outputs": [],
"source": [
"s2cube = connection.load_collection(\n",
" collection,\n",
" temporal_extent=temporal_extent,\n",
" # spatial_extent=spatial_extent,\n",
" bands=bands,\n",
")"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "f2ea33c3-d22f-46cb-8e9c-8c223e987196",
"metadata": {},
"outputs": [],
"source": [
"s2cube = s2cube.mask_polygon(geometry)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "4c36fa64-9065-4ddb-b708-04e3fd3f166a",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"CPU times: user 27.7 ms, sys: 16.2 ms, total: 43.9 ms\n",
"Wall time: 24.1 s\n"
]
}
],
"source": [
"%%time\n",
"s2cube.download(\"forum626-mask-polygon.nc\")"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "cf2ec21e-fa6a-4986-9529-7990014ed8a3",
"metadata": {},
"outputs": [],
"source": [
"plt.rcParams[\"figure.dpi\"] = 50\n",
"\n",
"def show(nc_path, var=\"var\"):\n",
" data = xarray.load_dataset(nc_path)\n",
" data[var].isel(t=0).plot(figsize=(5, 3))\n",
"\n",
"# show(\"forum626-mask-polygon.nc\", \"B03\")"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "3feeac09-eb79-405c-858d-5e9ee1ab10af",
"metadata": {},
"outputs": [],
"source": [
"green = s2cube.band(\"B03\")\n",
"swir = s2cube.band(\"B11\")\n",
"ndsi = (green - swir) / (green + swir)"
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "f7cbdec2-8bfb-4376-a38b-8abd7df382f8",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"CPU times: user 35.1 ms, sys: 15.5 ms, total: 50.7 ms\n",
"Wall time: 28.5 s\n"
]
}
],
"source": [
"%%time\n",
"ndsi.download(\"forum626-ndsi.nc\")"
]
},
{
"cell_type": "code",
"execution_count": 11,
"id": "894cc117-2907-4c4e-a861-614da1895566",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAPoAAACdCAYAAACdHWXfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/av/WaAAAACXBIWXMAAAexAAAHsQEGxWGGAABDZklEQVR4nO2dd3gUdf74XzvbN9vSE5IQCD00pQnSEUSxnCIKllNUlCI2/OJhQT3RU1BPLIgK2Auiwul5ehSlCJ5KkZpQA2mkl+195/fHkv0ZsouBAAk4r+eZ50km79nPeyb7nk97F5koiiISEhLnNUJzKyAhIXHmkQxdQuJPgGToEhJ/AiRDl5D4EyAZ+hngyJEjrFq16pSvnzx5MhdffDH9+/dn9erVADgcDsaOHcugQYOYN28eABaLhX79+qHX69m9e3f4+gULFtCvXz/69evHl19+2eDzS0tLufTSSxk4cCAfffQRAFarlauvvprhw4czc+bMiHqNGDECs9nMN998Ez735Zdf0rdvXy666CJef/31U75niTOMKHHaWbt2rfjQQw+d8vWHDh0SRVEUq6urxb59+4qiKIrz588XFy1aJIqiKI4ePVosKioSvV6vWF5eLt52223irl27wtdnZ2eLPp9PdDgcYq9evRp8/gMPPCCuXr1a9Pl8Yv/+/UWXyyW+8MIL4uLFi0VRFMWpU6eKv/zyS4PriouLxSeffFL897//HT7Xp08f0WKxiIFAQOzevbsYCARO+b4lzhxSj34GWLhwIZ999hnDhg2jurr6pK/PysoCQK1WI5PJAPjpp5+49NJLARg1ahT/+9//UCqVJCYmRrze5XJhs9kwm80N/v7rr78yYsQIFAoFffr0Yffu3Rw6dIgLLrgAgF69erFhw4YG17Vq1arBuU6dOmG1WnG73Wi1WgRB+kq1RBTNrcD5yNSpU8nIyODFF1+sd37Xrl3ce++9DeTXrVsX8XMeeeSRsHxNTQ1GoxEAk8l0whfIFVdcQZcuXQgEAixZsqTB330+X9gg6z4rOzubH374gd69e7NmzRratWvXqHsdP348/fr1Qy6X8/jjjzfqGomzT4t//ebl5XHnnXcybty4iH8vKSnhvvvu495772XTpk1nWbuTo3v37qxbt67BEYl33nkHv9/PLbfcAoDZbMZqtQKhuXlcXFzE66xWKwsXLuTAgQPs3buX2bNnI4oil156KcOGDWPXrl0olUqCwWC9z5o0aRK5ubmMHDkSvV5PSkoKn3zyCcOGDeOhhx6Kek+zZs1i165dHDx4kA8//JCampomPCGJM0WL79GzsrJYsmRJ2NAXLlzI/v37qa2tZc6cObz88ssYDAbsdjvp6enNrG0IpVJJIBBocL6xPfqaNWv48ssv+eqrr8LnLr74YtasWcMdd9zBmjVrWLRoUcS2BUFAq9Wi0WhQKpV4vV5EUay3ONi3b1/WrVvHkCFD2Lp1K/PmzUOr1fLOO+8AMGnSJK688kratm3LTTfddMJ7ValUGAwGVCoVCoUCt9t9QnmJZqK5Fwkay3XXXSfabDaxf//+4pNPPik+8MAD4ldffSWOHj1a3Lt3r1hdXS3eeeedza2mKIqiaLFYxMGDB4vjxo0Ta2pqTvr6rKwssVevXuLQoUPFyy67TBRFUbTZbOJf/vIXceDAgeJzzz0Xlr388svF1NRUsX///uK7774riqIozp07V7zooovEvn37igsXLmzw+UePHhVHjhwpDhgwQHz//fdFURTF3377TRw6dKg4fPhw8Z133omo1+233y62bdtW7NmzZ1iHDz74QOzbt6/Yv39/cfbs2Sd9rxJnh3PK0K1Wq3jdddfVO3/nnXeKZWVlotfrFW+++eZm0k5ComUjE8WWHdRSVVXFY489xurVq5k0aRJyuZzCwkJEUeSOO+5Ao9Ewb948ZDIZkyZNYuDAgc2tsoREi6PFG7qEhETTadKq+7p16xg8eDBTpkypt6CUm5vLPffcw3333UdOTg5Wq5XrrruOu+66ixkzZkSUkZCQOHM0adVdJpOh1+txu931VrxffPFFkpOT8Xg8pKSkUFhYyIUXXsjjjz/ODTfcEFEmGmPHjqVNmzZNUVNColEcOXKE5cuXn1Cml8yEn4aD4KsevZdnn332TKnWZJpk6IMHD2bo0KGUlZUxY8YMPv74YwC2bt3K+vXrKSoqYv78+TzyyCNs3ryZyy+/nIsuuiiizNNPPx2xjTZt2vDPf/6zKWpKSDSKutHmCZHBJUJ8g9Mul+sMaHT6aNLQvc67KjY2Fo/HEz6flZVFTEwMsbGx2Gw2/vOf/3Ddddfx3XffcfToUaqrqxvIHM/KlSuZMWMGR44caYqKEhKnFbkMVIKswdHSaVKPvnz5clauXEltbS3Tp0/nr3/9Kx9++CEPPPAAkydPxuv18vjjj2MymZg+fTo///wzPp+P2NjYBjLHM3r0aEaPHt24t6yExFlCLots2L5m0OVkaJKhjx07lrFjx4Z/HzZsGABDhgxhyJAh9WS/+OKLer9HkpGQaOkIgFLW8nvw42nxLrASEi0J4djQ/VxDMnQJiZNAIZOhlbf4WLAGRDT0UaNGERsbW++cKIqkp6fz8ssvnxXFJCRaItHm6C2diIY+adIkxo8f3+D8Z599dsYVkpBoyYRW3Ztbi5MnoqHXGbnFYglvmyUlJUU0fgmJPxPn3Rz91ltvRavVolKpkMlkvPrqq2dTLwmJFsl5NUcHSE9P5x//+MfZ1EVCosVzXs3RAeLj4xk3blw4T1ld9hEJiT8zMgHkisb16A6Hg2nTpqFSqRg2bBg333wzAEuWLGHz5s0UFhbSo0cPnnvuOdq1a8eoUaPo1asXd99992nXO6qh79q1i3feeQe1Wn3aGz0Zlm0tRPJ0l2gpCHIBuUreKNnly5czbtw4rrrqKsaPHx829DvvvJM777yTBx98kIkTJwKg1+txuVxkZGScEb2jGrrZbCYnJyfco2dnZ58RBf4IvVbZLO1KSERCJsiQN3LZvaioiO7duwMgl9d/Objdbg4fPkynTp0A+O233xBFkSuuuILLL7/89CrNCYJaYmNjWblyJZ9//nkD99WzidPTMMmihERzIShkKDSKBseRI0eYMWMGK1euDMump6dTVFQEEM66W8cXX3xRz31cEATkcjkajaaB7OkgYo9eVFTEk08+GfH82c60Ggic/puWkDhVZELkoXukcOqxY8cyffp0/vOf/3DVVVeFg74APv/8c5YuXQrAvn37mDt3LhCKFzkTRTAiGvrkyZO54oor6p0TRZGdO3fy1ltvnXYlToRwDq5wSpy/CHIZcmXj5ugxMTG8++674d/r5uhAvVTenTp1OuOL3REN/eGHH44oPGjQoHq/r1u3jtmzZ9O1a1cmTJgQjl7Lzc3l9ddfRy6XM2XKFNLT07n99tuJi4vDYDCE33wrV67knnvu4eDBg1EVTDJpTuW+JCTODIIM4Rx0jYto6EOHDm3UxU1JJVVbW8u6devC9b6iUVbrYtjL61n3YON0kpA4kwhyWaNX3VsSTXo1DR48mO+++465c+fWm9Nv3bqVv/3tb9xxxx3Mnz+frKyscCqputX7Z599Nmp53t8TCIgYYlRNUVNC4rQhKASUGkWDo6UT1dA9Hg+//PILGzZsiFhZE049lVRxcTEHDx7k6aefZseOHeEa3b+nLpWU0l2F3dnS83dI/FmQCTIElbzB0dKJ+iq64447GDBgAEplaB87UjaYU00l1apVK1asWAGEVvLrCgn+nrpUUl+MGE+PzNgGf5eQaA5kggy58jyZowP06dOH6dOnn/DipqSS+qPzdXicPnYX1p5QRkLibHEynnEtiaiG/u6777JmzRp0Oh0ymYxly5adTb3CqLVKembGMvGTbbx3U69m0UFCog6ZQkChOffWjKIa+s6dO8+mHlEJiiL5VQ7aJeqbWxUJCWSCgKBq+YtvxxNV40WLFrF69WpkMhnDhw9nypQpZ1OvMAq5gFoh4PRKrrASzY9MLiAoI8RftHAHzhP26HXD9QceeOBs6dMAg1aB0xvA5nY0mw4SEnXIBAG5MoLZeBqeiham+tRTT5Gbm0tsbCxPPPEEqampTJkyBYVCQVZWFg899NBp1zvq8qHFYmH9+vVs2LCB6urq095wY3FLPblEC6Ju6H78EYm6MNVFixbx9ddfh88rFApUKhVKpRKz2czGjRvp1q0bCxYsYPv27Xi93tOud1RDnz9/Prt372bXrl3Mnz//tDfcWEQRysod2Bxe/vrRlmbTQ0ICQFDIkWtUDY5I0WtFRUXh+PLfh6k++uijfPjhh4waNYrFixfXk0tKSqKqqur06x3p5LJly8KRNUC9n8823kCQGKManVZJokHye5doXkL76IoGR1302ujRo8Oy0cJU6xzNkpKSsNvt9eQqKiqIj29YxLGpRBxztGvXjuTk5P9/c81cgkalEPD6g1hckoecRPMik0dbdW9YSjlamOo//vEPCgsLqays5NVXXyUlJYVPPvmE+++/n549e6JSnf7tu4iG3rt3b/7xj3/w6KOPAvDkk082W500xbEwVZVCoKja2Sw6SEjUEVqMi5T1qOG8OlqYap1d/Z6FCxeeNh0jEdHQr7/+enJycti+fTvAGXnDNBa1Uo5KIWBz+qQkFBLNTt0cvSGnfwHtdBLR0D///HM2bdrEwIEDz7Y+DXB4/Li8AdQqOXJBweg3NrFyWvPrJfEn5Rx1mIm66r5q1arwz88+++xZUSYSgkxGIPj/5z9+v9SrSzQfMkFAUCgbHC2dqIZeU1MT/vlMLPc3FoNWQcAfpEOygZoaN/Zad7PpIiGBIIBC2fBo4UQdg7Ru3ZrbbrsNQRCaLdUzQDAIKpWcaoeH1ql6jlZJC3ISzYdMkCNTnkdBLTNmzKC4uJjY2Nioi3FNyRk3c+ZMrFYrtbW1LF68GIPBELENh8ePXqMg1aRld7EFv1cauks0I3I5MtW5588Rdej+wAMP8NRTT6HX67n33nsjypwoZ5zJZEKpVNbLGbdo0aKwY8ALL7zAW2+9Rf/+/cOr+5HITNARCIps2lOGy+0nOVHHJa/+eIq3KyHRNGSCHJlC2eBo6UTt0eVyOZmZmQCYTKaIMoMHD2bo0KGUlZUxY8YMPv74YyCUM279+vUUFRUxf/58HnnkkXDOuIsuuih8fWlpKVu2bOG+++5r8NkrV65k5cqVHDp4mAsvEclKNxIIiqgUcsql4btEcyETQHHuDd2j9uhqtTo8BP/9wly9i08xZ1x1dTXFxcXMnDmTBQsWNChXA6FUUv/85z9BnxBylql0IhdkuLx+vC5/U+9bQuLUkMuRKZUNjkg4HA5uu+027rrrrnAnCPD8888zefJkrrvuOoqKijhy5AgXXHABU6ZM4csvvzwjakft0Z999llWr16NKIpMmzYtosyp5oyLjY1l1KhRdOzYkVmzZjF16lR69uwZVUm7249ep6RHhpnf8mvQm8+9OZLE+YFMJiBrZI8ercjirFmzAFixYgVr165l8ODB6PV6nE4nrVu3PiN6RzT0+fPnc/ToUUQxtH+9fv16Bg0axJVXXllPrik547Zu3dooBVVKOXqNgiqLmzV7ygCorZBi0yWaCbki4mJcXfRaXVJTOHGRRbvdzrJly3j77bfR6/Vs3LgRp9PJhAkT6oW0ni4iGvr48ePxer3hYJZgMMjcuXMbGPrZQCkX0KnkeGNU4eCW7l2SuHPpbyyZcOFZ10fiT44gRByqR6q9VheVdsEFF9SLXrNardxzzz3Mmzev3m6TTqc7Y2pHNPTU1FS+/vprPvroIwRB4MYbbzzjTvfRCIoi23eW0a5jPNVWD33bxWPSKSmRHGckmgNB3ujFuGjRaxMnTsTn8/Hss89yww03oFQq+eCDD3A6ndx0001nRO2oc/RvvvkmnEpq6tSp/OUvfzkjCvwRaqUcU7wOu9uPzxPA4w9SZffSNc3YLPpI/LmRCUKjt9OiRa8tX768gezgwYNPj4JRiLrq7nK5KCgooKCgAIej+ebEDo+fC9rHo1IIyATIr3KQHqelwhYhSZeExBlGlMkR5aoGR0snao9+//3389prrwEwe/bss6bQ8Zh1KnYX1NAm2YDD4eXKHqlYPX72FFuaTSeJPzEyWWj4fo4R1dB/+OEHXnjhhbOpS0ScHj+ZBjW1zlC870ebjvDGjReyt8TG49/l8szlXZpZQ4k/FTIBUX7uhalG1fjrr7/m0KFDYa+4efPmnTWlfo/F6sHm9OF2+tDoQnOjl9YdYtuuUtIzzc2ik8SfGJkMhPPI0D/++OPw9lrdfnpzoNIoGNAhgX99f4g+AzMJBEW2H6zClKBDpTj3it1JnOPIBMTzydA3bNjAN998g9/v59prrw37vZ9tkkwaurUyYhuYyd4SK3qNArVWgcPqwaVt+cEEEucZMgHk5973LmqXmJeXx2effcaXX37J//73v7OpUz30GgVLNxdSYXMjF2ToVHLapRjITDeSatY2m14Sf05EITRHP/5o6UQ1dKvVyk8//cSGDRuQyWTk5OScTb3CtIsNGXOiQYNckOH0BnB5Q0EtWqWcmz7Y3Cx6SfxJkQmhOfrxRzNw8ODBRstGNXSj0cjq1atZu3YtiYmJf1jH/Eyxp8xOklFDtcMT7sGPVrvQa5TsLbFSJnnISZxNjs3Rjz8iES16bffu3dx8883cfPPN7N69GwgFutx3333hgJfGEMnxJhpRX0VPPvlkoz/kTFK3reb1B3lvsIy7/6fC4vSxv7AWlVqBQXfuzZckzmUav48eLXrtlVdeYcGCBchkMh5++GEee+wxfD4fr776KjNnzqSwsDBcoulEfPDBB3z//fcYjUZkMlnYkzUSTRpzNCWV1Nq1a3nvvffw+/288MILtGrVKmIbNbVuvP4g/qDI9M2gVyvwefyIQUgwqtm7t7IptyAhcVKIgrzRnnDRotcsFgtmsxkAm81GcXFx2LBbt25drxbbiagbDTSGiEP3xYsX89///vcPL25KKqk333yTd999l0ceeYQlS5ZEbSMlIYat20u4umcrSmpd7MyvQaGUk5yow+MP0qNbUqNvVkKiqYhAEFmDI1KRxWi110wmExaLBavVisFgIC0tLSxXWFhYz5ZOxNdff8348eMZM2YM11xzzQllI/boK1aswGKx8MMPPyCKIjKZLKLDTFNSSYmiiCAIZGZmhm/y99SlkrKUF5M5PIYvthZh0CjIahuHWadkfW454/plYPf4efb7/Tx2ScdGPRwJiaYgihCI4FcSKUw1WvTa/fffH87D+PDDD9O6dWuUSiUzZsxArVY3qjcH+O6778jOzuaRRx5p0PbxRDT0pUuX8sEHH3DVVVed0FnmZFNJ3Xrrrdx9991UV1cjCALBYJCCgoKIb7C6AP7PBl9PpwwzE/tl8Mq6Q+wsrEWrkpOREMO3O0vQqeRcmBnbqAcjIdFURFGsV1DkRESLXuvWrRsffPBBPdnnnnvupHWJj4/H4XDw66+/snfv3hPKRjR0g8FAnz59ePzxxwGYNm1aRIeZpqSSuvvuu5k0aRI+n4+5c+dGVVCQh/LEdU/SMaJLMtd2TmDipzsoqXWh14TUd3kDjX44EhJNIQD4WkjG8c2bNzNq1Cg2b97M1KlTTygbdTFu4cKFLF68OJwzrn///g1kmpJKasSIEYwYMeIPb0Z2rJrqgl+KGNU+kUVbionXqzhYYMGrU6JSyWklOc5InCWiDd2bg//+97/88MMPfPnll2zfvr1ehuXjibqPXjc3FwShWX3dtUo5VXYv1XYvGwuq+flQFfvzahjcPYXUOC0qhcCbXzePM4/Enw9RhECw4dEc1NbWcuDAAWpqakhNTT2hbNQefdq0adx+++0ATJ8+/fRqeBI4PH66pZkosbhIi9PidPvp1jGB/WU2vP4ggaBIty6JUpVVibNCELHF9OiPPvoo48ePZ/LkyeEAtGhENfSLLrrohEOBs8UNvdLZfKiKp6/uyqxlO0hP0SMXZBg0CmqdPlQKgW3bSkjNkhbkJM48oR69ZRj6yeRxPGfiPF/64QAqjYIUkxazVonXHyQzPgavP8jooW3xeaSiDhJnnqAI/mDDo6UT1dAdDgdbtmwBoLi4+KwpFAmvO2TEep2SVLOGOL2KHhlmdCo5DocXrz+I0+5tVh0l/hyIhBbjjj9aOicssli3B/iPf/zjrCkUiYu6JFFy1EbftnF4/UHyq5x4/UFqXT4EuUB+lYPEZH2z6ijx56AlLcadDFEN3WAwkJycDIBG07wlkPYUW0htZWBbfg3f7yyhqNrJpgOVPDGyPcKx+bpBo2Dk6xubVU+J859ghN78nO7RExIS+Omnn/i///u/iEUQzybpcToGdkhArRBQqhVkxsfg9wV4/L/7GdE1mf2Hqqmye6XUUhJnnEBQxOMPNjhOhqNHj3LzzTdz6623snbt2np/mzlzJpMnT2b8+PHYbDbWrVvH4MGDmTJlCuvWrTtlvaOuut92222MHTuWYDAYtWzy2SI+RhWuuzawQwIZsVqSjGo6JxsosriIMarRaxQ8PLJDs+opcf4TRMQXbNpYfcmSJcyaNYuuXbtyyy23MHz48PDf6jIvv/zyy2zfvj1q4NjJErULXLJkCZ07dyY7O5tXXnnllBs4HVzaOYlb+mdy1YWtSDSqyU7S88PuUt74/gByQUZKvI5Us5Z5aw40q54S5z/BIPgCYoMjUvQawK5du7jyyivrHXl5eWRkZIRjRY6ntLSULVu2cPHFFzN48GC+++475s6d26QcERF79Ouvv56cnJxwvGvXrl1PuYHTwff7K/j1UBVX907jsg4JeAMi79/amydX7ueLn/IRg5AZH4M89swVqZOQgNAc3RdhHz1S9BpA9+7d+eabb+qdmzNnDkVFRRiNDcuKFRcXM2vWLBYsWFBvynx84NjJEtHQP//8c4qLi0lLSzvlDz6d5Fc50B0LYLln2U6u6Z1Gql6NXJAxaUR7lIIMX1BErRDo+/dVbH7y0mbWWOJ8JSCCr4nL7HfeeSezZs1CoVAwadIkgHBA2NVXX03Hjh2ZNWsWU6dO5dChQ/UCx06VqHP0Dz/8kBUrVhATE0NlZSU7d+485UaaytEqJ72y4vn3b0fx+wL8cqiKi9snkB6rJSlGRYndw6aDVagUAmqNlFpK4swRCIq4m+gh06pVqwZhqh9++CEQyuXwe3r27FkvcOxUiTpHLy0t5bLLLuOHH37gtttuiygTbUUwNzeXe+65h/vuu4+cnBxycnKYMmUKU6ZMoV27dkBoDeCWW25h7NixbN++/YRK3jq4LdvyqojXq1Ao5Zi0KiwuHyadkv/sKeNAmR2vP0B6rFbykJM4o9QN3Y8/WjpRe/SEhASCwSDPPPMMO3bsiChzolRSycnJeDweUlJSiIuL480332T79u3ExcUBsH79et599122bNnCunXruOCCC6Iq+cnPBbRJNlBU7aRrmokt+yooSYrB4vSRHqejb5tYym0eAkERQ5wUsipx5giKIv5zwUPmOKIael3Sid9++437778/okxjU0k9/fTTQCgX3YMPPgjAhAkTGD58OH6/ny+//PKESrodPlRJAr0yY9lRWEuH1mbGdEthZU5oy82oVpCVGINckOF2+E7yEUhINJ6gyDnRgx9PVEN///33WblyZTgu/ZNPPmkg09hUUgBOp5Pi4uLw0P3NN99k48aNFBcX88ILLzB//vx6n12XM+7IkSMYO6rZub8STYySfu0TOFRhZ8nGw3RNM2Fx+fjnd3uZOrIDBbUuendMaPJDkZCIRkAUm7wY1xxENfS9e/dGNO7f09hUUgCfffYZ119/ffjaoUOHctddd2G1WsMrj7+nLmfcjBkz0GXGUmTUsH13KbkxKpxuPxkJMVTY3AzvnMTBYgvvbDpCVpKePXnV3L9iF69c2/1Un4mERFQCQRH3+WToPp+PzZs3ExMTA0B2dnYDmZNJJVWXxKKOhx56qNFK/u9AJa3itIhBMGgU2J0+ap1eMuNj2LC/kknD2/PjwUp6Z8aiVyuwSwtyEmeIgCjiPRfiUo8jqqEbDAa+/fZbILTo9sQTT5w1pY6nVZwWuSCQ0MpAYZmdYEDEnKQnyahmb4mV7/eWY3P7qLR72FNsIcmobjZdJc5vAsHzyNAPHz7M9ddfH84V90dpas40cTFqzDolBwstBAJBendK5LcDlcTFqNCq5KgUAg+NaI9epcDi9HF9z8hVXyQkmkrwfDL0VatWUVZWVu9cc/boBo2CjfsqcDu9GGK1GDQKumaFtumqalzE69W892shl2cnM7ZHKjM+286PM4f/wadKSJw8/oDY5PTiR48eZebMmcjlcm6//fZ6QS0TJ05EoVCgUCh45ZVX8Pv9TJs2DZVKxbBhw8K54U+WiA4z6enp9OnTh969e9O7d2/69Olzand0GqmtcKA3a3BaPeSWWBnZOQmLy8ekEe3RKuUcOGol1aDmzU1HAHh69b7mVVjivCQgNj1MtS567b333mPRokX1/qbVapHJZJjNZpRKZbhQ46JFi/j6669PWe+IPXplZSVfffUVAwYMQKVS8eOPPzJmzJhTbqSp1Dp99O6eQl65HYVKjs3pY2VOGXcNyKRbopZtBbXIFQJOX4DOqQb+0iOVBF3jCuFJSJwM0ebo0b5tu3bt4pFHHql3LjExMWr02oIFCxAEgVdffZVvvvkmaqHGkyWiod92223k5uYyc+ZMAAKB5q2EYvf4Kap2olII2LwBMlqbuX9IFoeqnSz4MY/HRnVk6qfbeWvjYWxuPxv2VfCvCe2bVWeJ85OAKOKNsL129FiYat22cB0nG71WZ/xJSUnY7fZwocYLLrigXqHGkyXqqrtMJuOOO+5ApVL9YXL4M41Jq8SUZmJ7XhXmBB0Pj2iPwxsgTqvk3REGArn/ZtWAGObZs/n8+4MM6p3OiLf2sPFvUqVVidNLMEjEHj1amGokThS99tBDD+FyuaipqWHx4sUA9Qo1nioRDV0URfr27cvIkSMBIr55ziajOicx+6NtPH9rb7JitdS6/Sz6OZ+J/TIIHt5B0FrFtXvb0r+dyJuT+9NHVUVwpOQwI3H68QeCuLxN89M4UfTaSy+91ED+94UaT5WIi3EymYzdu3fjcDhQKBQ4nc4mN9QUhm5eiEIl56PNBXx7oIJql4/embE8umI3CHI+TvkLXn+Qcd1T6bnrUyg9hKJsH+Pf+7VZ9ZY4/zgdOeOag6hD98OHD/PYY4+Ff3/nnXfOikKR2DZwOq39h3h3mJZXD4fCAtvFhbLJ/CtmAN/8WsB/B9qRuQ7iU2nwHtyJsnVH4vRtm01nifOT884z7u2332bVqlWIoshll112NnVqwEDrFv6TmMHt61wkGkR+OVTF9CGh4JhrMgQu2/ot3jwjVxT249v2bkS/F5lCKZVTljjtnKuecVETT0yZMoXq6mpqamqYPHny2dSpAUNXqnhmSDIHjtSSV2Gna5qJTfnVGGJUXPNZHjK5gDw2iVmXd8Z+4CA1e/O5ZF0M+eV2rnvnl2bVXeL8QgyKBPzBBkdLJ2qPbjQa+etf/woQNfHE2eIf1/eAnWtYN+MqXvpfMXJBhl6joLzUzqfT+qP4+Tec+3Lo8csvxHRsyyudr2BCrJZ/bz+KWafiqrf/h9XmISVeR8HRUNhs23QjJdUu1lwXS566Da9vOiJFvEn8IcGgiN937o0Uoxp6bGwsd955JwCZmZkRZdatW8fs2bPp2rUrEyZMCEev5ebm8vrrryOXy5kyZQoAr776KgCrV6/m0KFDDWQiRcfV8eL3B3C601iZVUD31FhKbB7eXX2A+KQYMgKV+DwuFDEaTHEmAg47a3KL6dExgeFdkthZWItZp8Tm8FJUaqdVip7yKidltW4WTLgAKjezrExNICjyzpYCdhdbsbh8HCqxIgZF2qYa6ZxqIDvZwIJ1h0g1awkEg3x4S/N7C0qcfcSgSDBwHiWeeOKJJ8JJIwwGQ0SZpqSSOl7mRCz/SwpCVT5iaR4F3mxMaiUje4cy1MptZfgARVwi+L1o2/dAflSgpNaFXJBxuNhG93ZxmA1qqvwuzDoVFqcPl93Lc2v2s/jaflQdygdgW34t+VUOUs1axKCIQimnqNqJXq0gyaAmyahGpRDommbmrx9tYdu2EmSCjFuu6ERehQO720dxhYOYGBXVlU5+fnxk4/8TEucEosg5MVQ/nqiG/swzz5CfHzKAjIyMiEEtTUklFU0mEnJHFaLWiBifSUKpinaxOj79tYD+7eK5eKmduTdOZKB1C4IpgZLYLng9m5ALMsxaJYFAkPxKB4GgiNmkoaTWhUmnpGuaiXmXtmHBtvLw4kqFzU3ftnGs3VOG3xfE4/ITY1RTYXOzYlsxz16ZzYNf7KTK7qGoyEpmx3gK82pYvauUjq2M7NpdjjlRR2ZCDPF6FUNfWofRoCYrUU9RjZNembFU2b3YPX7yyu34PH5aJ+ul0cE5hBgUCZxPiSeqq6vDDvd1rrDH05RUUpFkfs/vU0mJMoEn92o5UFbBsM5JfHW0lA7JBuxuPx3amJnz3V66pqWx40g1fu+vpKUZiderSTFriI3TolXJqbZ6kAsyqiud9OiYgFmnRMhZS7qpDxv2VVBR6USQyygqtYd1UKrlJBjV9GkbR6xWyWPf5FBVasOpV9OzcyIWpxeVWo5MkPHjz4WYE2Nw2b1s31+JXC4QY1ChU8mxe/xU1rj4ptROdlYcVXYPnVON7C2xUmX3MvqNTaFn5PCSkawnLkbN2J6pdEvUYfcFmfGvPSQZNcgFGQvH9Typf7DE6SUYDOI/B3dzohp6bW0t77//PjKZjKqqqogyTUklFUnm9/w+lVT/DyppnWGi9lgN9BsvSGNXuY1X/rOX5KQYasrs/FztIilFT2FFDYFAkCGdEvl6azEZCTGUW93UlNvRt4ll+IWtyE41crjKwRvebvx2oBSvP4hCKaDRKkMOES4fuhgVXm+AKruXaruXNXvKMOuUmOJ1ZGfGUlLrwmL3ckHXZPYcrMKUoMPt9NKtUyLlVjdWi4e8PeX0vSgdlULggqx48qscmLRKbu2bwVubDiMGRTzHvjQyQUbGsdLPmQk6kmPUPPtDHhU2N4GgSO/MWDQKgXKLg9xKFy+vPYhBo0AuCASCQV64qguvbAyNwJ67Ivp6h0TTEIM0eY5+ojDVefPmkZeXR05ODjfffDOjR4/mmmuuoX///owaNYrrrrvulNqUiWLkmq8+ny+8jz569GiUyuYpjDBjxgzs/W8lyaBmd7GFqhoXuhgVHo8ft8NH56w48svtxBnVlJY7SEzQ4fIGSDKqeWRkR17beJi88tCWXLxeRVGNi3Krm9JyB+kpespqXKjVCizVLvRmDWJQxF7rRmtQcW2/DG7omszD3+SiUsg5Um7nxv6t2VpQy/6jVnQaBU63ny5pJvKrQtODzqlG9hRb6JpmosLmxukNYHH66Jxq5JY+6SzbfpRSi4tyq4fMhBgOlVhRqxVoj0Xl1TGgQwI3XZhGhcOLzRsgTqtg4cbDPDKyI+9tLiTFrGFtbjk3X9Sa9Qcq2b6/EkGQodEpaZdiYE9eNZ3bxJJXGhot2S1u+nRNRiUXUCsEln+3j8suac/r1/Volv9rS2TGjBl/6K+e0OcaUoY1zHF4afCHRvu6z5kzh2uuuSZcZDFSbsYbb7yRhQsXUltbyy233EJWVhb33nsvffv2bdzNHEfUHv2NN94Ip3l+4403mDZt2ik1cDrokKznP9uKUasV6GJU2K0eBLmMIT1T2bSnjMxWBrQqBUntQ3XckwxqXN4Ai38pIDNeh1wmIz1OS36lk2q7B4vThyCXUWn1EG/SkJ9fi96soWuaiReScvnrgTaM6JKMLxjkzqU7CARFtCo5bZL05JTY6JFuIie/hvJaN4JCYPOeMjpmxZJ7oAqrxUO/zolU2NzUOn0YNAoyE2LISorBqFYwN+Ugf1d14miVk7xSG+ZjuiYZNRg0CgorHKjVCjYfrmZMl6RwauF0oxqb288bm47w4uXt+PvafFLNWuweP7mFtUwY2pZfDlWRYtJSYXMz8sJWlNs8+H0BurWNo0f/1hyucFBicXFNz1b4R3ekwuZmxle7cfkCBIIickGG3e0j1aRFpRD4W8VnFK3dRodn5qJI69Js//+WRDDKHP1IfuTotZMNU4VQj6/VajGbzZhMJjZu3IjT6WTChAmnHJMe1dALCgrCP9ctyjUXeRUOZDIZWUl6SmpdOGQyzCYNhyrsdGhtpsruCfsbp5q01Dp9DGkfzy9HathZWMvQTknolHK2HK6m0urB4/IjE8Dr8VPu8SMTQqmy0uO0lHcZg+rwXvaV2riyazJ5qcZQ7TeVnB4ZJlzeAHaPn+zM2PDIwHvMUPp2T+FAqY28CjvxejUGTaiH31FYy/RuMRz0QHnHUfBzAf+OX8vcxPEEgiLb8mvo1yaWz7cUodMqcbp8tIrXsa3ESoe4GAwqOan+SkZ2TWbNnjLGfrADlULAoFHw6c8FmA1q8sodqBSheGW9RklRjQu9WoHBoEYlF9hbYiO7lZHuaSae/ng7F/VqhV6j5MIMM/vKbBwos5GXX8tF3VKocnjRqxXMib2ev18bR8UHr/Nmr3sx6ULehi5vAK1Kzs7CWuL0auJjVOSWWPnyjoua7TtythCDIsGTiF472TBVCLmbT5w4Efj/adx0uqYVEI1q6HK5nNmzZyMIQpPiYE8HOfk1KJRydh+uRq1VIlfIUCkEioqstGltwusPkpWoDw2fRZG8MjsuX4CiaidDOiWy8Ks96E1aYgwqrDUuBJkMpVpB93Zx7DxYhVwu0DEtZMTbyxwUlIUW5F5dd4huaSaOltoRBBlefxB/UMSgUZIZr6Oo2knAH6RTh3h2bC8lNSuWeL2K8lo3Uwdn8WNeFV5/kHaJeg64dWw5asGoVpBXYWecfRTe6gr8vgBqrZJPfy4gwaimtMpJrEnD0SonOzUWxnRI4P2txTyVW06cUU16rI6t+ypwO7zEpxrwuPz0zIzF7vGjVysIBEUsTi8mnYoKm5vR3VNY/nMBglwIv7AMsVqq7B5q7V5c3tBzSo/ToVDKyUqKocLqoajGSVyMmud1l1LUdjDx3gAqhYDd7cfjD5JXYadzqhG9RsHeEhsDOySw86iFv766idjkGOJNGtQKgc6pRi5IM7Fg3SE6JBuQCzIsLh/v3dSrWb9Tp4oYzWHmJGa2JwpTFUWRjRs3htetfvzxRz744AOcTic33XTTKesddY4OIccXURRP6MxyppkxYwb/S7iCVil6DBol5VY3LrefJLMGjz9IVY2LdmmmY19CHyOzkym3hXYAVu8qpV+7eGxuP95AkLxyO2UFFjKyYqmpdjFvwgU88fUeumaY2ZJbzqdTLmLKsp2kx+kot7rp3y6ebfk1GDRKTFolpRYXcXo1Zq2S/CoHNrefeL2a4Z0S+flwNbmFtSTHaik4aqN9axOZ8TFsy69BpRDQqeTcN7QdPxXUsOlAJR2TDfy4uwyjSY3HG+Dq3mnY3X7yq5wUVTuxWjzEGELTlK5ZccgFGTn5NbgdPvRmDcGgiCiKDM1O5kCZDbkgIz1WRyAo8lNOOTEGFalmLfuO1OD3BXj0uu4sWHeI9DgdO/aUcdnFmWw+XE2vzFi2Hq6mVZwWrUqBxemla5qJViYNOSU2Si0uemSY0ark5Fc6SY/T4vUHWbOnjJFdkymqdoXb75lhptblo6jaSXmpHb1ZQ7JZQ5Xdi9cbIDNJT5JBjVYlp5VZy/p95VzZPZWDlQ4OlNnw+oPY3H46p4Z6uqykGKYPOHuBSY2Zoxu7XUFc/1sbnB9r/KXRc/TmIGqPDtClS8uYl7VK0RMIilTZPWEjvKF3OkVWNyt3lXJJ5yTe/+kIWpWcr7cfJRAUuePiNnxUYGGD00dWupGqYyv2/Xq14ki5neSkGA5UO+iaYWZ3QQ09Oibw4L/2AKBSCEwenMV/9pQSr1ezr7CW1sl69h+o4t5ru5JX4aBzqpFap4+tByvZva+C9AwTf7+6K0//J5fkpBj0GiW1x2rD/fRrEa3bhxyFxnVNpvqYLpoYZWgEYPWwNrccm9OH7Jh7b6dMczju+cBRa/hZCAqBYFCktsKBJkaF1x8k0aDh19xySmrd4VTXY3un4wsGKbe6yUyIYdXecrQqOSW1Lvpf2Ioh7eLZXxZaqFOr5KTH6sJz9B2FtewtEdCqFOGXx55iK0kGNfmVTrKSYuiQYmBviQ2tSs6QTonoNQqW/XgEjU6JvdaNzqjm6t5prPi1kEGdk9Cq5HzxwyFSM0x0TjXy7uoDmOK0fLXjKJnxMQSCIk5vgKxEPdUOL+VWNwOz4pj+5U62H6zCaFJTUWYnLc2Ixe6lb7t4Lkg38dLy3Wydc/aCrkQxQNDvPWvtnS5OaOgthSSjBrvbh0ohJ7/SQXqcjiKrmz3FVlxuPzplaFFjSKfE0Dmvn+9yypArBPRGNU5vALNOyYWZsSxZvofEdCOZGWbW5JTx2KhOfKQJDXkPHLVSU24nEBR5Zm8FYy7KYFDbOJ7Mr+FggYXxI9vj8QepdfqosLkpr3XTOlmPMzbUy81dvZ9AIEiXVCObD1fTI8NMfpWDCy9IobDCwbLtRxndOYmnOznYqmjHz4equDAzltHtE1HIIafcwXc5Zbi8fo6Uh6YPBp0SuSDD6fLhdvgwJ+jomWHml0AQp92LzR16GXRvH09RtZP8ozZMZg37Sm3sLbGSHhcy1NRjzkJZSXriY1R8sqUQs07FoQo7/dvFsz63nC5pJjz+ICOzk0mKUTFvxR5GHPNAzIzXYXH6kAsy9pbYwqm9hndK5NtdJZTXuumaFUd+pQNzgg5rrZvlPxeQlhjD7mILFruXzLaxZMbHUO3wYorTMrp7ChVWD7uLLRg0Csy60EgpyahhSKdE5q3YQ3q6EVEUKS+xodWruK5XOsu2FFJu82DzBjDG6+h0778YOrgNRdVO4vUq4mLUeANBzFol+8tsVFnc9GwTh04lb/LWoxgIEvSde4YeNXrt1ltv5bvvvjubukTF4vSyP68Gl9dPklFN51QDxdUuJg9oQ3qCjm/3lGHSKam2eymqduL0BujXJpa27eMY2CGBXpmxYQeb4YPb0DMrHrkgY2inJKpdPmpdPnYXWxjYKZGHr+9B1zQTKSl6yq0enlyxm/jYUIXWTon60IvEF0AuyBhzYWhBq2OygYry0AvIVu1i97EiEi5fgCGdEikos9Mu1ciz8vXolHLcv60jp8IeKh6pEPhXTilPfLuX1384yI69FZRbPSQY1Zj1KmxOHzXVLjwuP067l9pKJ+t/OxraWmwXGiXIBRkur58OyQa0ehVZiXrSj1WVHdYhgfgYFetzypjQNwO5TEa3VqERTpJBTbxezfc7SzDrVXgDQZKMGrYcqeG5L3aRnm7EpAv5Fhwos6NVyfH6g+wtshCvV5MZH8OK34opLnPgsnvJK7fj8fjx+4OotQoCfhGb249ckNE53UQgKJJf5aDK7sHvD7Kz0MLuYgvd0kzoNaF2buybgdcfYMO+Cu66rCM6lRy/N0DHdnEYDGp2FVsY2TWZWqeXjQcqUKnk9LsoHbkgIxAUKbd62LT9KNV2DyUWN0lGDZZKJ9UODzsKa7l72fYmfRfFoJ+A19XgaOlE7dEXL17M0qVLmTBhAgMGDGDSpEnh8kxnmyu7p3Iw1ciGfRVkJenZW2JjdHYyxcf2qRMNmtDiWI2LXpmxpMVpySmxkX+4JtSTJOnJjNeRV2GnQ7KB2mP71RqFQJxWybSBbXj310I27as45tSiIhAUqXZ4QltT6SYAlm8vpqTWzQOXdOCllfvw+INU1ri4qGMiQ3qmsmFHCYY4Ldf2SmNviY1f9lcwtEMCl17Qij3FFq460pMF/ZR8knE9Pl8ApVrO2txyDJpQr33DRa1Z9ksBSUY1RZVOvG4/fm+AXtlJHCi1IQZFHFYPCa0MmPQhHV2+AFqlnFqnj5x9lai0Sg4INnYdqSbWpKHY6qaoxklthQOnL4DLF+CHfRVkJsTg8Qc5VGqjU4aZy7OTef/nfOSCjHi9mo7t4jDpVFTbvcgFGZnxOmqdoRdTj8xYimqcjOmWQn6Vg9aperQqBeVWN6jkuLwB4vWasINTu2MuwC5vgA4phtD6R9dkVu8sISMxhg7Jejb/XIBeo2DR+jx6ZMZSUuvm058LcNq99O2a3GC///6BWY3/Ap1Gb0IxGCBwPg3dq6qqyMvLw2g0kpKSwp133snSpUvPpm5hPvoltNWnUghkxuv48OvcsANNqllLepyWLzflM7xnKrVOH//aXASAXCFQcriG9ITQS6BDsoEDZbbwFOD+QW3YVmJje7ElNJw81htu3FcBgPlYbza6fTxPrzpAudXD01d2YUepDZNeRVmFE78vQLnNg8vrJy5Bh8cbYO3ecnIPVBFjVLP+QCVDOySwemcJ8SYNJrWcFIOGD3/JD1fCyUqMYdW2o6wHAoEgT2x6nhkXPoRcr8Ji95KTX4PH5cdlD7ncup0+VCo5mfEx7MyvwaRXUV3tok3bWGqPLXylJcZwtMrJz4eqiNeriU3S8873h2idqqfW6WN09xQ6J+jZkV+DxenlYJUDg0aB0xugR4aJH/dVhFx85aFBn83tR3esKk6Vw4vF6eO9n47g9QeprXDQrVMiXn+QzIQYSmpDwUPVVg//vntA1P/rY5d0DP98R5/WZ+Krc9oRg+fZHP3FF1/knnvuISsr9ObMyMg4a0odj9cfxO8LMH5AJks3HeHFKf15be1BPrrlQj7ZVUpeuYOuWXFsy6+hc6qROKOaIZ0SWbOnDFdiTHiRSq9RkGjQ0CnFwMR+rUnRiIxuF0uyXsWRGhff7Cph04FKEoxq5ILsmHupyJP/3YfN7cda6+b1DXmUVLsQBBkKlcDlfdJYs6uUDq2M4e2mnYW1qLVKqkrs5B7box/QMZH7B2VSavcxKMPAjqNmDqgUHCq20PfiNvTNMDP/h4O0itfRZuItZBbFsKewFr8vgFwh4LJ7MSeGtsA8rtBw2O7xM3FwW1b8VkxcnJZauxe1Sk6vzFjW7y5FF6NiZHYyB8rsZCbpOVxipaTahUolZ+n6w/h9AXp0SqSo2klehYPyY/EAv+XX0iHZwDOXt4zF2JaEGDw35+hRDf34bJQXX3zxGVcmGknHvMKW/VzAbUOzeGrpdsyJMbyw/jBxehUuX4ByqxuDRsFPu0rp0yUJlzeAzeElOVaL0xsIOZFUu3D5AhRWO5mYbUKwV5Lrj8V9zAGiuMIR9mMPeYoFqLV5cNm8dGkfT5UyFPoqCDICgSBuh4/Nh6sZ2CkRm9uPSiGwNrcclUJAkMswJejwe0PRc1cMbYdeKfBDXhXf7ijh2wlteD1Xjsvr5/GPf2Purb15/+YLePvXIu44IKNrmpbtB6tQquX4fQFSWptINWs5WGjBZfegVIeuffVfe8jINDOwQwKzhnf4/w9NMtIzwnnXo7ck+rSNY+WuUiYObsua3DJ6dU3G6w+GF3Lkgoz+7eKxOH2U1biwuLyY3UqSY7XE6dX0SYxhT7EVO6Ge0OkN8Oq2Sga3iSOnwsrwNrHMXbUfQS5QbfcQr1ezc18FKrWC9BQ92mQD+eV2MpL12NyhRbFLjrmY7smrJv/Yqv3Pu0rRGlQ4AUGQodYq8biOucGaNazYW8XQNvH0zzAz8ds8IFSF5tL+rRnZIRGAJ0Z1Ct93PcOVaBGIwQABr7u51ThpzglDVykEbumfyeK1B7ljWDu0Sjl7SqzsO1KDJ9mAVimnpNaNNxCkW+tYNu8pY3jnUPGG9btLmXBhyBll474KMhJi+GlLMWMGt+HTbcU8NiKLH46Ehqo6lZyiGmcoGi1eR48MM7/klpOYoEMURYorHCTEamndKjTX/2xiv5O6j+6ppvDPF7eJP63PSOLsEAwECEQcurdsU2rZ2h3jl2P7zSnxOvIqHWTEatl8uJrsrDi8gSB2h4+sRD15FXZsbj+PjeuOIJOxu9jK36/tRpnDy9XZKWw+XE1cjIodz9evIzehZwwTeqZFbvwvZ+EGJc4dojrMtGxTatnaHSM9VsemA5WUVznDveikfpHz2P2eK7sk1/v92ynNt84gcX4gBoNRDL1pQSdnmhZv6EeOHKHN94vJBDKBGTP+3bTPatPmdKnWotqT7u30tPNHzLrrL7hcDR1ktFrtGdDoNCL+iXjwwQfP2/ake5M4EVFdYM9Hfp8Q4HxrT7o3iRNxwjBVCQmJ84M/VY8uIfFnpcUvxkUjGAwye/ZsrFYrffr0YenSpWRmZqLX63nxxRcjZtp86aWXOHz4MD6fjzfffJOSkpI/lKlL5fNH7b333nssW7aM1q1bc88999C9e3dmzZqF0+lEp9Px/PPPs3v3bp577jkAHnnkEbp169ZAJlJb7dq1Y+nSpcjlcmbNmoXRaGTatGmoVCqGDRvGzTffzCeffMLatWvxeDwsXLgQ4A9l6oKU/qi9ffv2NajIc6rP8vi2cnJysFgs/PLLLzz//PMMGjTotN6bRIhztkf/6quvKCoqQqlUkp6ejk6nIxgMkpwc2lJbsmQJs2bN4r333mPRokV4vV62bdvG66+/Tvfu3dm4cWOjZBrbniAI6HQ6AoEAycnJFBQU4PP5ePXVVwkEAhQWFvLKK6+wYMEC3njjDV577bWIMpHamj9/PjExMcTExBAXF8fy5csZN24cixYtCicLXLFiBYsWLeKGG25g+fLljZKJdm/Ht3d8RZ6mPMvj25o7dy5vvvkmaWlpjBw58rTfm0SIc7ZH37dvHxdffDGTJ09m3LhxfP755wiCwIwZM9i5cydFRUX1Mm1WVVWRmBhyM83MzKSoqKhRMo1t75ZbbuHWW29l586dPP/881x//fXhQKDWrVtTVFSExWLBbDYDYLPZKC4ubiCTkZHRoK0dO3awdOlSVq5cyccff0xZWRndu4cKQsrloYSQdSOPzMxMdu3aBdAomUj3dnx7EydOrFeR58UXXzzlZ3l8W5dccgm//vorvXr1Qi6XU1RUdFrvTSLEOdujp6enExsbC4T+2XVfsKSkJOx2O+np6RQVFYUTW8bHx1NZWQmEMtymp6c3Sqax7R3/e1paWvjLXVhYSHp6OiaTCYvFgtVqxWAwRJSJ1FaXLl1QKBThqjZ1egMNEncef28nkol2b8e3d3xFnqY8y+PbglDugzvuuCP899N5bxIhztlVd6fTyb333otOp6Nz5878+uuv6HQ6/H4/b731FqWlpeFMm7fccgsjRozgn//8J/n5+eF5XElJyR/K1PUUf9Te4sWL2bZtG1VVVTz55JN069aNRx55BI/Hg1qt5rnnnmP37t3MmzcPgIcffjiiTKS2EhMTWbt2LQ6Hg5deegmdTsf06dPRaDQMGjQoPEf98ccfcblcLFiwAOAPZermsX/U3o8//hiuyDN16lSGDRt2ys/y+LZuvfVWbrvttvBw2+FwnNZ7kwhxzhq6hIRE4zlnh+4SEhKNRzJ0CYk/AZKhS0j8CZAM/Qzx+9Xg9957j2+++YY333yTQ4cONer6cePGnXSbVquVMWPG8PbbbzdK/vjlmaNHj/Lyyy+flH5btmzhs88+Oyk9161bx5gxY+rV94vE0qVLGTly5El9tkRkztl99NPJihUrKCwsJDs7m19//ZVHH300/LeZM2fi9/vDHm+TJ0/GZDIRFxfHE088wd///ndqamqwWCy89dZb3H333bRp04bu3buzefNmRFFk3759TJo0idLSUlwuF0899RQ2mw2FQkHnzp25/fbbmTNnDpWVlaSkpPDXv/6VnJwcnnrqKSZNmsT7779PZWUlNpuNl19+GYPBAMDPP//MG2+8AYQ8xbZs2UJtbW297aWnnnoqvJ3Xrl07hgwZws0338xVV13FJZdcwiuvvIJGo+Gqq66iZ8+eFBcX4/f7eeyxx/B4PAQCAebPn88///lPCgoKMJvNjB49mpycHJ555hkGDRpERUUFVVVVPPDAAxiNRnr06MHkyZPp3r07EydOZOvWrSxZsqReKOeYMWNo3bo1Tz31FDU1NdTU1NCjRw+qqqooKiriww8/ZMKECXzxxRdn6VtwfiP16MC1117L/v37eeutt/jb3/4WPp+Tk4NKpeLll1/mwQcfZNWqVQwZMoT58+eTn5+P1WrlyJEjzJ8/n2HDhrFq1SoA7rrrLkaOHElpaSlz585lwICGKY9vuOEG5s6dy6pVqwgGgwSDQYxGI8uXLyc9PZ3s7OzwC2HDhg2YTCaUSiW5ubnhz3jttdd4++23efvtt1mwYAFXXnklF198MWPG1M+gM3bsWP7+97+Hq3pmZ2cza9YsPvvsM+bMmcPixYv56KOPwvKrV6/myJEjmM1m7HY7xcXF7Nq1i9dee405c+YwaNAgsrOzw4UAAT799FPuuOMOFixYwJo1a4DQnvhDDz3EgAED2L59e9TnP378eF577TV++uknnnvuObRaLTU1NSfxH5T4I6QendAQtra2FkEQ8Pv9YUcOURTr1bAWRTG8r348MpksPBQ2mUwEAgFUKhUAarW6gXzdPq8oivz222/IZDLmzJkTdhX9vY99165deeqppyLqXSd3ol1Svz9Utsnn84X1O/76399XMBhk4MCB3HfffUDIi+/4+z7+90jPpu4elUolHo8nqn5GoxG1Wh0uI6xSqU4oL3HySD068Morr3DjjTfy8MMPM3v27PD5rl274nK5mDlzJq+//jqXXnop69ev56GHHiIjIwOj0UhmZiYPPfQQa9eurRc3bTKZSE1N5aWXXmLDhg0nbL99+/bs2LGDF198kYqKinDb//d//4fRaAy72k6dOjXsDw8hp5EpU6YwZcoUpk2bFvXzly1bxtSpU7n22mvrnZ88eTKzZ89m8uTJ3HjjjeEX2+jRo9m5cyczZ87krrvuQqlUkp2dzYMPPsicOXOAkPH+/lndeOONLFmyhPvuu4/hw4c34qlLnE0kh5nznKeeeopx48bRrVu3P5SdMmUKY8aM4eqrrz6jOq1bt47du3czffr0P5QdN26cNE8/DUhD9/OcSEP+aLz55ptnTpHfYTab2bFjBwUFBbRuHb0U09KlS0lLi5KdV+KkkHp0CYk/AdIcXULiT4Bk6BISfwIkQ5eQ+BPw/wCMNDic3029PwAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 250x150 with 2 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"show(\"forum626-ndsi.nc\")"
]
},
{
"cell_type": "code",
"execution_count": 12,
"id": "e89bc9bd-66b2-49f5-ac59-5cd26c999edf",
"metadata": {},
"outputs": [],
"source": [
"mask = ndsi > 0"
]
},
{
"cell_type": "code",
"execution_count": 13,
"id": "8da5963c-2d89-4b4f-865d-e72a3a68b52c",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"CPU times: user 8.74 ms, sys: 1.72 ms, total: 10.5 ms\n",
"Wall time: 50.2 s\n"
]
}
],
"source": [
"%%time\n",
"mask.download(\"forum626-ndsi-mask-binary.nc\")"
]
},
{
"cell_type": "code",
"execution_count": 14,
"id": "56121fa8-c725-403c-bf56-5bfb2d1fece5",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAPIAAACdCAYAAACOyiUrAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/av/WaAAAACXBIWXMAAAexAAAHsQEGxWGGAAAil0lEQVR4nO2deXyM1/7H35NIELJRtcVS7qu3DVpXLbEloQhKqxFC7USTWIKQNkpIi2qsUVuIfauluFyuBiXoaqkleQW9qGVi+VVWIrLN+f2R5qlkZmSyTSbJeb9ez+uVeeb7POf7TOY755zvOedzVEIIgUQiKdOYlbYDEomk6MhAlkjKATKQJZJygAxkiaQcIAO5BLh9+zZHjx4t9PXe3t506NABJycnjh07BkBKSgru7u506tSJBQsWAJCUlETbtm2pXr060dHRyvUrV66kbdu2tG3blr1792rd/+HDh/To0YOOHTuybds2AJKTk3n//ffp0qULAQEBOv3q2rUrdnZ2HDp0SDm3d+9e2rRpQ7t27VixYkWhn1lSRISk2Dl58qSYOnVqoa+/efOmEEKI+Ph40aZNGyGEEKGhoSI8PFwIIYSbm5tQq9UiPT1d/N///Z8YMWKEiIqKUq53dHQUGRkZIiUlRbRq1Urr/pMnTxbHjh0TGRkZwsnJSaSmpoqFCxeKdevWCSGE8PX1Fb/++qvWdbGxsWL27NniP//5j3KudevWIikpSWRlZYkWLVqIrKysQj+3pPDIGrkEWL16Nbt27cLV1ZX4+PgCX9+kSRMAKleujEqlAuCnn36iR48eAHTv3p2ff/4ZCwsLatWqpfP61NRUnjx5gp2dndb7Z8+epWvXrlSqVInWrVsTHR3NzZs3admyJQCtWrXi9OnTWtfVq1dP69w///lPkpOTef78OVWrVsXMTH6lSoNKpe1AecTX15cGDRqwaNGiXOejoqKYOHGiln1kZKTO+0yfPl2xT0hIwMbGBgBbW9uX/kC89957vPnmm2RlZbF+/Xqt9zMyMpSAy7mXo6MjJ06c4J133uH48eM0bdrUoGf19PSkbdu2mJubM3PmTIOukRQ/Jv/zeevWLcaMGYOHh4fO9x88eICfnx8TJ07kxx9/NLJ3BaNFixZERkZqHbrYsGEDmZmZDB06FAA7OzuSk5OB7L5xjRo1dF6XnJzM6tWr+d///se1a9cICgpCCEGPHj1wdXUlKioKCwsLNBpNrnt5eXlx9epVunXrRvXq1alTpw47duzA1dWVqVOn6n2mwMBAoqKiuHHjBlu3biUhIaEIn5CksJh8jdykSRPWr1+vBPLq1av5/fffSUxMZM6cOSxduhRra2uePn2Kg4NDKXubjYWFBVlZWVrnDa2Rjx8/zt69ezlw4IByrkOHDhw/fpzRo0dz/PhxwsPDdZZtZmZG1apVqVKlChYWFqSnpyOEyJV8a9OmDZGRkTg7O3PhwgUWLFhA1apV2bBhAwBeXl706dOH1157jY8++uilz2ppaYm1tTWWlpZUqlSJ58+fv9ReUkKUdifdUPr37y+ePHkinJycxOzZs8XkyZPFgQMHhJubm7h27ZqIj48XY8aMKW03hRBCJCUlic6dOwsPDw+RkJBQ4OubNGkiWrVqJVxcXETPnj2FEEI8efJEfPDBB6Jjx45i/vz5im2vXr1E3bp1hZOTk9i4caMQQoiQkBDRrl070aZNG7F69Wqt+9+/f19069ZNtG/fXmzevFkIIcTFixeFi4uL6NKli9iwYYNOv0aNGiVee+018fbbbys+bNmyRbRp00Y4OTmJoKCgAj+rpHgoU4GcnJws+vfvn+v8mDFjxKNHj0R6eroYMmRIKXknkZQuKiFMe9FEXFwcM2bM4NixY3h5eWFubs69e/cQQjB69GiqVKnCggULUKlUeHl50bFjx9J2WSIxOiYfyBKJJH+KlLWOjIykc+fO+Pj45ErYXL16lfHjx+Pn50dMTAzJycn079+fsWPH4u/vr9NGIpEUniJlrVUqFdWrV+f58+e5MsaLFi2idu3apKWlUadOHe7du8e//vUvZs6cycCBA3Xa6ONVVX2qUq0obhYI94//NFpZEtPizp9d2Ldv30ttPptUg9Tn2o1Yq1fGMW/evJJyLV+KFMidO3fGxcWFR48e4e/vz/bt2wG4cOECp06dQq1WExoayvTp0zl37hy9evWiXbt2Om2++OILnWVUpRqvq94uipsFYvHnl41WlsS0mLagcb42T59nERKsPYY/Y2FqCXhkOEVqWufMDrK3tyctLU0536RJE6pVq4a9vT1Pnjzh8OHD9O/fnyNHjnD//n3i4+O1bPISERGBv78/qaQUxUWJpFjRCA0ZIkvrKG2KVCPv27ePiIgIEhMTmTBhAsOGDWPr1q1MnjwZb29v0tPTmTlzJra2tkyYMIFffvmFjIwM7O3ttWzy4ubmhpubG3uXHtJRskRSOmQhyEBT2m5oUaRAdnd3x93dXXnt6uoKgLOzM87Ozrlsv/3221yvddlIJKaOBsgQ5SyQJZKKhgZBBqY3YisDWSIpAJlCxXOhMtj+1q1bzJs3j6SkJL799ltGjRqFpaUl6enprFu3jkePHhEQEIC5uTmjRo2iS5cuhfJLZyB3794de3v7XOeEEDg4OLB06dJCFSSRlAeym9aGB3LeRT8bN24EYNKkSdy/f59NmzYRGBhIs2bNGDp0aPEGspeXF56enlrnd+3aVahCJJLygkaoyNChx3H79m38/f2VJO3LuHbtGmlpaTRo0AC1Wk2DBg2KLMigM5BzgjgpKUkZVnr11Vd1BrdEUpHIQkU65lrnGzduzJIlS/K9Pjo6mtDQUFatWgWAg4MDarVaEY0oLHr7yMOHD6dq1apYWlqiUqn4+uuvi1SQRFIeyMKM5xrDU0s5i34uXrzIvHnzWLFiBb169cLPz4+ZM2cyZswYAgMDqVSpEl5eXoX2S69HDg4OfPnll4W+sURSHslCRYaOGlkfNWvWJCwsTHk9Y8YMLZstW7YU2S+9gVyzZk08PDyUKj9HPUIiqchohBkZwvBANhZ6AzkqKooNGzZQuXJlY/ojkZg0GlSkC9MbtdXrkZ2dHTExMUqN7OjoaDSnJBJTJaus1cj29vZEREQA2csVZ82aZTSnJBJTJQszngvL0nZDC52BrFarmT17ts7zpqJUKZGUBmWqj+zt7c17772X65wQgitXrrBmzRqjOCaRmCJZqMpOIH/yySc6jTt16pTrdWRkJEFBQTRr1oxBgwYpq5+uXr3KihUrMDc3x8fHBwcHB0aNGkWNGjWwtrZWBs4jIiIYP348N27cKMZHkkhKDo0wKzvJLhcXF4MuLorUT2JiIpGRkcp+QxJJWcBUa+QiTfDs3LkzR44cISQkJFef+sKFC3z66aeMHj2a0NBQmjRpokj95GS/582bp3f7TonEVMkSZqRpLLSO0kZvIKelpfHrr79y+vRpnTvzQeGlfmJjY7lx4wZffPEFly9fVvbofREp9SMxRTRkJ7vyHqWN3sb+6NGjad++PRYW2b82utQ8Civ1U69ePfbv3w9kZ8JzNip7ESn1IzFFslc/lX7g5kVvILdu3ZoJEya89OKiSP3kd14iMUU0QkWGpgwF8saNGzl+/DhWVlaoVCp2795tTL8kEpMkC3PShOF94rwKITt27ODkyZOkpaWxevVqAMaNG4elpSWurq4MGTKkUH7pDeQrV64U6oYSSXkmS6jI0GinlvTV0XkVQvbv38+ePXs4dOiQIobv4eFB37598fT0LHQg6012hYeHM3DgQDw9PXMtw5JIKjKav4af8h45CiE505r1oVJlywQ1atQItVqtKIQAmJsXvsn+0ho5pzk9efLkQhcgkZQnNEJFpo4+sqEKITncvXtXmXuhVqtp2bIlGk3hZXb1BnJSUhKnTp1CpVIRHx9f6AIkkvKEPs0ufbyoEDJ//nz69euHr68vqamprFy5EoAJEyZw+PBh+vbtW2i/9AZyaGgo33zzjfK3RCLJnhCSXgCpn7wKIQAfffRRrtc5yppFQedPy+7du9m5c6fy+sW/JZKKjPiraZ33KG10/rQ0bdqU2rVrK69zOugSSUVHg4rMAjStjYVOj9555x1+/PFHXFxccHFx4fvvvze2XxKJSZKd7DLTOkobnTXygAEDiImJ4dKlSwBYWpqeIoJEUhpkCRXpJtCUzovOQN6zZw8//vgjHTt2NLY/EolJI/6qkU0NvR4dPXpU+XvevHlGcUYiMXXKVNMaICEhQfk7Li7OKM5IJKaOQEWWzmRXltF9eRG9gdywYUNGjBiBmZmZlMKVSP5CI1Rk6ayBTTSQ/f39iY2Nxd7eXm+yqyiaXQEBASQnJ5OYmMi6deuwtrYukQeUSIoTjVCRkaUr2ZVhdF9eRG/jfvLkyQQHB1O9enUmTpyo0+Zlml22trZYWFjk0uwKDw9HrVYDsHDhQtasWYOTk5OSHZdITB0hVGTpOEobvTWyubk5jRo1AsDW1lanTefOnXFxceHRo0f4+/uzfft2IFuz69SpU6jVakJDQ5k+fbqi2dWuXTvl+ocPH3L+/Hn8/Py07h0REUFERISU+pGYFPqb1qWLXo8qV66sNJFfTHzluriQml3x8fHExsYSEBDAypUrdS7fcnNzY8mSJVSlWlGfUSIpNgSg0ai0jtJGb408b948jh07hhCCcePG6bQprGaXvb093bt35/XXXycwMBBfX1/efvvtEntIiaS4EAWoke/evYufnx81atTg9ddfp2HDhrnUQapVK75KSiWEEHlPhoaGcv/+fXLeMjc3p1OnTvTp06fYCjaURqrXeV1lvCCPuH/ZaGVJTItpC/rku6a4Rt8OvDqqp9b53j8kal17+PBhEhISGDp0KJ6enmg0GkUdJCEhgWHDhhWb7zprZE9PT9LT05XFEhqNhpCQkFIJZInElBBk95PzkqMQkqP+CuDk5ISHhwcbNmxg2LBhHDlyBMhWB4mKiipWv3QGct26dTl48CDbtm3DzMyMwYMHK0JhEklFRggVQkefWJdCyMaNG/n8889xdnbGw8NDySm9qA5SXOjtIx86dEiR+vH19eWDDz4o1oIlkjKJwODkVs+ePQkODmbHjh00btyYVq1aaamDFBd6Azk1NZW7d+8CkJIih4AkEvirRs4yLJCbN2+updueVx2kuNAbyJMmTWL58uUABAUFlUjhEkmZQ6CzaV3a6A3kEydOsHDhQmP6YhK41Svfw2AyK180RFkL5IMHD3Lz5k1lVteCBQuM5pREYrqowASmZOZFbyBv375dGX7SMdQskVRMBFCWauTTp09z6NAhMjMz+fDDD5V51xJJhUZjeLLLmOida3br1i127drF3r17+fnnn43pk0RiughQabSP0kZvICcnJ/PTTz9x+vRpVCoVMTExxvRLIjFdNCrtowS4ceOGwbZ6A9nGxoZjx45x8uRJatWqJfcxlkjgrz6yjqMEyNmt0RD09pFnz55dLM4UKyozECbQjpFUWFRChcpIya4tW7bw/fffY2Njk+8e5YZvYqODokj9nDx5kk2bNpGZmcnChQupV69e/gXKIJaUNsJ4feLo6GiDbXU2rdetW8d3332X78VFkfoJCwtj48aNTJ8+nfXr1xvssERSquQMPxmhj3zw4EE8PT3p3bs3/fr1e6mtzhp5//79JCUlceLECYQQqFQqnRNCiiL1I4TAzMxM2fA5L1LqR2KKqApQI2s0GoKCgkhOTqZ169ZYWFgUSFjgyJEjODo6Mn369HzXSesM5J07d7Jlyxb69u370skgBZX6GT58OB9//DHx8fGYmZmh0Wj0LunKWde5d+mhlz6ARGJUcpJdBnDgwAHUajU1a9bEwcGBsLAwRVhg3759+QoL1KxZk5SUFM6ePcu1a9deaqszkK2trWndujUzZ84EYNy4cTonhBRF6ufjjz/Gy8uLjIwMQkJCDPtkJJLSRqhQGTgh5Pr163To0AFvb+9c65ENFRY4d+4c3bt359y5c/j6+r7UVm+ya/Xq1axbt07R7HJyctKycXd3x93dXXmdk+xydnbG2dk5l23e4auuXbvStWvXfB9GIjEp9DStdSmEODg4KJrw5ubmSuvWUGGB7777jhMnTrB3714uXbqUS4E2L3oDOadvrFKp5FxrieQv9PWRdSmEuLu7M3HiRM6cOYOzszP29vYFEhZITEzkf//7HwkJCflOkdYbyOPGjWPUqFEATJgwId9CJZIKQQGSXVZWVlojMgURFvjss8/w9PTE29tbWcCkD72B3K5du5dW5RJJhcVIDdSC6OQVaUKIRFLRUGlAVbr7telE71zrlJQUzp8/D0BsbKzRHJJITJ0ytfpp8uTJbNy4EYAvv/zSaA5JJKaMykSXMeptWltbW2NnZwdAlSpVjOWPRGLaGHGudUHQG8ivvPIKp0+fZtq0aTo3WZNIKiKm2kfWG8gjRozA3d0djUajd1tViaQiYoo1st4+8vr163njjTdwdHRk2bJlxvRJIjFdjCgsUBB01sgDBgwgJiZGWQ/ZrFkzozolkZgsIjvhpUUp6/HpDOQ9e/YQGxtL/fr1je2PRGLS6F3GWMppJL195K1bt7J//36qVavG48ePuXLlijH9kkhME6En2VXKgay3j/zw4UN69uzJiRMnGDFihE6byMhIOnfujI+PD5GRkcr5q1evMn78ePz8/IiJiSEmJgYfHx98fHxo2rQpkN0HHzp0KO7u7ly6dKlYH0oiKSnK3DjyK6+8gkajYe7cuVy+rHu/oJdJ/dSuXZu0tDTq1KlDjRo1CAsL49KlS9SoUQOAU6dOsXHjRs6fP09kZCQtW7Ys3ieTSEqCAo4jp6Sk4OLiQnBwMNevX+ePP/4gIyODsLCwfBdCFAS9gZwjKnDx4kUmTZqk08ZQqZ8vvvgCyNYCmzJlCgCDBg2iS5cuZGZmsnfv3mJ7IImkJCmI1A9ASEgIAwcORKPR8Ntvv7F9+3ZWrFjBDz/8QOfOnYvNL72BvHnzZiIiIpR1yTt27NCyMVTqB+DZs2fExsYqTeuwsDB++OEHYmNjWbhwIaGhobnuLTW7JCaJEKg0hi1/OnbsGI6Ojjx//pykpCRq1aoFoFenrijoDeRr167pDN4XMVTqB2DXrl0MGDBAudbFxYWxY8eSnJyMl5eX1r2lZpfEFFHpSXbpUgiJjIwkJSWFmJgYzM3Nsbe3B7IVQt56661i9kuP/Me0adPw9PRUlP4cHR2LtWBDaaR6nddV5XvPYmMi90fWz7QFffJVq6zj6ELDNh9one9UU6332k2bNvHKK6/w+++/c+fOHUVF0yh9ZGtra/773/8C2UmtWbNmFVuhEkmZpRCLJkaOHFkirryIzkD+448/GDBggKLVVZy/HBJJWSY72WV6GnY6A/no0aM8evQo1zlZI0sk2YFsllnaXmijM5AdHBxo0KCBrJElkryUpfXIjx8/5sCBA7Rv3x5LS0vOnDlD7969je2bRGJyqAow/GRMdAbyiBEjuHr1KgEBAQBkZZngSur8kFuwSkqCslQjQ3ZzevTo0VhaWlK3bl1j+lQwVC9MF38xcGUQS0qAMpXsEkLQpk0bunXrBoCNjY1RnSoULwZuTnDLYJYUNxqBKtP0Alnn6ieVSkV0dDQpKSlUqlSJZ8+eGduvPA7pXaSVja4glkhKgrK2+umPP/5gxowZyusNGzYYxSGd5ASqrn7vi+9JJCWMijLUtAZYu3YtR48eRQhBz549jemTNjkBrKupnLePLANaUpKYaNZa77fex8eH+Ph4EhIS8Pb2NqZP2rwsQPX1g2X/WFICZMvhCq2jtNFbI9vY2Cg7qusTFjA6+Q0pFSbJJYepJAVBYBKBmxe9gWxvb8+YMWMA9O7NGhkZSVBQEM2aNWPQoEHKRudXr15lxYoVmJub4+PjA8DXX38NZK/RvHnzppZNvqur9NXKsiktMSIqYXgN/O9//5vDhw+TnJzMmDFjiIqKMr5CyKxZsxRRAGtra502RZH6yWvzUmSwSkwFARjYR+7Xrx/9+vUjISEBf39/0tPTja8QMnfuXO7cuQNAgwYNdC6aKIrUjz6bApM3ay2b1ZKSRAhUmoJ9Z+bOnYuXlxd79uwBSkYhRG9VFx8fT3h4OOHh4UrNrHVxEaR+dNm8SEREBP7+/n9L/ejLWueQ3/sSSTGg0ghUmRqtI0chJCIiQrEVQvDpp5/Sq1cv2rRpw+PHj4FshZAXW7DFgd4aOTExkc2bN6NSqYiLi9NpUxSpH102L5JL6qekArQCBr5bPam2oo8Wkw0wEgJ09JEbN26spRCyfPlyjh8/TlJSEjdu3KBVq1ZMmjSJtLQ0xo0bVzxO/4VeqZ+MjAxlHNnNzQ0LC4tiLdhQpNSPxFi0mNwgX6mfxrWceKNed63zju+m5HttSaK3ab1q1Sree+89+vTpQ3h4uDF9Kl1KOrFWHPc3VvJPZfb3IclGCNBotI9SRm/T+u7du8rfOUmvCkFJNbfzzk4raKLtRXtjdQnkTDltCjD8ZEz0BrK5uTlBQUGYmZmhMYFfnDKPvjnihb3eWFTAPMJL0QjINL31+XoDecGCBVy9ehUhRKlJ4UokJkdO09rE0BvIAG+++aax/JBIygZCgAkq5rw0kCUlhJyIUnbRCMjS9b8r3VyC3tKHDx/OkSNHjOlL+Ucmjso+QpNdI+c9Shm936x169bx559/MmjQIJYtW0ZKitxMrcgYO+ssKX5ymtZlJZDj4uK4desWNjY21KlTR1kJJTECsuY2XXKa1nmPUkZvH3nRokWMHz+eJk2aANkLJ8odsq8qKShCIEygBs6L3kBevHhxrtcdOnQocWeMSmnXei/7EZE/LqaLzFqbGKUdLKVdvqRwaDSITNPb/KniBrJEUhhystYGkJKSwrhx47C0tMTV1ZUhQ4aUmFsyqyKRFAChye4j5z10sW/fPjw8PAgPD+fgwYMl6pcMZImkIAjDA1mtVitJYnNz8xJ1y+Sb1u982JzGjYsnY3779m0aN25cLPcytfLksxVPOfkx9LMBpKamap1/8OAB/v7+iiAGZG9PrFaradmyZckvPBIViClTppTb8uSzmR5Pnz4VI0eOFD4+PmLbtm0lWpbJ18jFSc4vZXksTz6b6VGtWjU2btxolLL0Sv1IJJKyg0x2SSTlgDLbtNZoNAQFBZGcnEzr1q3ZuXMnjRo1onr16ixatIj79+8TEBCAubk5o0aNokuXLixevDiX0v+DBw/ytcnZDSC/8jZt2sTu3btp2LAh48ePp0WLFgQGBvLs2TOsrKz46quviI6OZv78+QBMnz6d5s2ba9noKqtp06bs3LkTc3NzAgMDsbGx0Rqf3LFjBydPniQtLY3Vq1cD5GtTrVo1g8q7fv261o4ihf0s85YVExNDUlISv/76K1999RWdOnUq1merKJTZGvnAgQOo1WosLCxwcHDAysoKjUZD7dq1AVi/fj2BgYFs2rSJ8PBw0tPT+e2331ixYgUtWrTghx9+MMjG0PLMzMywsrIiKyuL2rVrc/fuXTIyMvj666/Jysri3r17LFu2jJUrV7Jq1SqWL1+u00ZXWaGhoVSrVo1q1apRo0YNneOT+/fvJzw8nIEDB7Jv3z6DbPQ9W97y8u4oUpTPMm9ZISEhhIWFUb9+fbp161bsz1ZRKLM18vXr1+nQoQPe3t54eHiwZ88ezMzM8Pf358qVK8oYXo6IflxcHLVq1QL+Vvo3xMbQ8oYOHcrw4cO5cuUKX331FQMGDFDGEBs2bIharSYpKQk7OzsAnjx5QmxsrJZNgwYNtMq6fPkyO3fuJCIigu3bt/Po0SNatGgB/D0+mdNyaNSoEVFRUQAG2eh6trzljRw5MteOIosWLSr0Z5m3rHfffZezZ8/SqlUrzM3NUavVxfpsFYUyWyM7ODhgb28PZP8zc75Ar776Kk+fPlXG8HLG72rWrKml9G+IjaHl5X1dv3595ct77949HBwcsLW1JSkpieTkZKytrXXa6CrrzTffpFKlSsquHDl+A1rjk3mf7WU2+p4tb3l5dxQpymeZtyzIXvs+evRo5f3ifLaKQpnNWj979oyJEydiZWXFG2+8wdmzZ7GysiIzM5M1a9bw8OFDAgMDqVSpEkOHDqVr164sWbKEO3fuKP2oBw8e5GuT80ufX3nr1q3jt99+Iy4ujtmzZ9O8eXOmT59OWloalStXZv78+URHR7NgwQIAPvnkE502usqqVasWJ0+eJCUlhcWLF2NlZcWECROoUqUKnTp1UvqIZ86cITU1lZUrVwLka5PTj8yvvDNnzig7ivj6+uLq6lrozzJvWcOHD2fEiBFKczglJaVYn62iUGYDWSKR/E2ZbVpLJJK/kYEskZQDZCBLJOUAGcglxIvZ1E2bNnHo0CHCwsK4efOmQdd7eHgUuMzk5GR69+7N2rVrDbLPmx65f/8+S5cuLZB/58+fZ9euXQXyMzIykt69e+faX0wXO3fupFu3bgW6d0WlzI4jFyf79+/n3r17ODo6cvbsWT777DPlvYCAADIzM5UZW97e3tja2lKjRg1mzZrF559/TkJCAklJSaxZs4aPP/6Yxo0b06JFC86dO4cQguvXr+Pl5cXDhw9JTU0lODiYJ0+eUKlSJd544w1GjRrFnDlzePz4MXXq1GHYsGHExMQQHByMl5cXmzdv5vHjxzx58oSlS5dibW0NwC+//MKqVauA7JlO58+fJzExMdfwS3BwsDLc1bRpU5ydnRkyZAh9+/bl3XffZdmyZVSpUoW+ffvy9ttvExsbS2ZmJjNmzCAtLY2srCxCQ0NZsmQJd+/exc7ODjc3N2JiYpg7dy6dOnXizz//JC4ujsmTJ2NjY8Nbb72Ft7c3LVq0YOTIkVy4cIH169dTtWpVxa/evXvTsGFDgoODSUhIICEhgbfeeou4uDjUajVbt25l0KBBfPvtt0b6FpRtZI0MfPjhh/z++++sWbOGTz/9VDkfExODpaUlS5cuZcqUKRw9ehRnZ2dCQ0O5c+cOycnJ3L59m9DQUFxdXTl69CgAY8eOpVu3bjx8+JCQkBDat2+vVebAgQMJCQnh6NGjaDQaNBoNNjY27Nu3DwcHBxwdHZWAP336NLa2tlhYWHD16lXlHsuXL2ft2rWsXbuWlStX0qdPHzp06EDv3r1zleXu7s7nn3/OoUOHAHB0dCQwMJBdu3YxZ84c1q1bx7Zt2xT7Y8eOcfv2bezs7Hj69CmxsbFERUWxfPly5syZQ6dOnXB0dMy1Qf0333zD6NGjWblyJcePHweyx4SnTp1K+/btuXTpkt7P39PTk+XLl/PTTz8xf/58qlatSkJCQgH+gxJZI5PdxExMTMTMzIzMzExlooIQQpkMkfM6Z1w5LyqVSmmq2trakpWVhaWlJQCVK1fWss8Z5xRCcPHiRVQqFXPmzFGmMr44x7tZs2YEBwfr9DvH7mWjiJl/icVlZGQo/uW9/sXn0mg0dOzYET8/PyB7Flre5877Wtdnk/OMFhYWpKWl6fXPxsaGypUrY2NjA4ClpeVL7SXayBoZWLZsGYMHD+aTTz4hKChIOd+sWTNSU1MJCAhgxYoV9OjRg1OnTjF16lQaNGiAjY0NjRo1YurUqZw8eTLXullbW1vq1q3L4sWLOX369EvL/8c//sHly5dZtGgRf/75p1L2tGnTsLGxUaaC+vr6KvOxIXtShI+PDz4+PowbN07v/Xfv3o2vry8ffvhhrvPe3t4EBQXh7e3N4MGDlR8uNzc3rly5QkBAAGPHjsXCwgJHR0emTJnCnDlzgOzgfPGzGjx4MOvXr8fPz48uXboY8KlLihM5IaScExwcjIeHB82bN8/X1sfHh969e/P++++XqE+RkZFER0czYcKEfG09PDxkP9kAZNO6nKOrSa6PsLCwknPkBezs7Lh8+TJ3796lYcOGeu127txJ/fr1jeJTWUfWyBJJOUD2kSWScoAMZImkHCADWSIpB/w/A/1dIwBpEd0AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 250x150 with 2 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"show(\"forum626-ndsi-mask-binary.nc\")"
]
},
{
"cell_type": "code",
"execution_count": 15,
"id": "12d9d067-e5e7-47e4-a7b5-c3b17a4df40a",
"metadata": {},
"outputs": [],
"source": [
"mask = (ndsi > 0) * 2.0"
]
},
{
"cell_type": "code",
"execution_count": 16,
"id": "b4af4255-c4e5-4d05-9ae2-0116ad1d4207",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"CPU times: user 5.78 ms, sys: 0 ns, total: 5.78 ms\n",
"Wall time: 24.4 s\n"
]
}
],
"source": [
"%%time\n",
"mask.download(\"forum626-ndsi-mask-float.nc\")"
]
},
{
"cell_type": "code",
"execution_count": 17,
"id": "246ab779-702d-429d-af71-16496d4197df",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAPQAAACdCAYAAACD1FVsAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/av/WaAAAACXBIWXMAAAexAAAHsQEGxWGGAAAl0ElEQVR4nO2daVRUV7aAvyoGFRk1DkQUG3vFBCX6jMZ5SjQkDukEcUichQg4i42NRpSOI2oiThEHnGM7RH36zEujPkFi0h017YAPY7dTBIw+BQRBZKrzfpC6LVRdhgKLAs+31l2LOrXv3fuU7nvGvY9GCCGQSCS1Am11GyCRSKoO6dASSS1COrREUouQDi2R1CKkQz8Hbt++zfHjx02+PyAggG7dutGlSxdOnDgBQHZ2Nj4+PvTo0YPly5cDkJGRwZtvvom9vT1XrlxR7l+/fj1vvvkmb775JgcPHjR4/r1793jnnXfo3r07u3fvBiAzM5P333+fvn37EhISYtSut956C2dnZ44dO6aUHTx4kE6dOtG5c2fWrVtncp0lVYSQVDmxsbFi1qxZJt9/48YNIYQQaWlpolOnTkIIISIjI8XmzZuFEEJ4e3uL5ORkkZeXJ/7v//5PjB07ViQkJCj3e3p6ivz8fJGdnS06dOhg8PwZM2aIEydOiPz8fNGlSxeRk5MjVqxYIbZs2SKEECIoKEj8+OOPBvelpKSIBQsWiP/6r/9Syjp27CgyMjJEYWGh8PLyEoWFhSbXW1J5ZAv9HNiwYQP79u2jT58+pKWlVfh+Dw8PAOrUqYNGowHghx9+4J133gGgf//+/O1vf8PGxoZGjRoZvT8nJ4fHjx/j7Oxs8P3Zs2d56623sLa2pmPHjly5coUbN27Qvn17ADp06EB8fLzBfS+//LJBWevWrcnMzOTp06fUq1cPrVb+l6pOrKvbgNpIUFAQzZs3Z+XKlcXKExISmDp1qoF8XFyc0efMmTNHkU9PT8fR0REAJyenUl8UAwcO5LXXXqOwsJDo6GiD7/Pz8xXH0z/L09OTU6dO8cYbb3Dy5ElatWpVrroOHz6cN998EysrK+bNm1eueyTPD4t/nd68eRM/Pz98fX2Nfv/rr78ybdo0pk6dyvfff29m6yqGl5cXcXFxBpcxtm7dSkFBAaNGjQLA2dmZzMxMoGjs3KBBA6P3ZWZmsmHDBv71r3/x888/ExYWhhCCd955hz59+pCQkICNjQ06na7Ys/z9/bl69Sr9+vXD3t6epk2bsmfPHvr06cOsWbNU6xQaGkpCQgLXr19n165dpKenV+IXklQWi2+hPTw8iI6OVhx6w4YN/POf/+TRo0csXLiQVatW4eDgQFZWFm5ubtVsbRE2NjYUFhYalJe3hT558iQHDx7kyJEjSlm3bt04efIkEyZM4OTJk2zevNmobq1WS7169ahbty42Njbk5eUhhCg2SdepUyfi4uLo1asXP/30E8uXL6devXps3boVAH9/fwYNGsTvfvc7Pv7441Lramtri4ODA7a2tlhbW/P06dNS5SXPmeoexJeXIUOGiMePH4suXbqIBQsWiBkzZogjR44Ib29v8fPPP4u0tDTh5+dX3WYKIYTIyMgQPXv2FL6+viI9Pb3C93t4eIgOHTqI3r17i3fffVcIIcTjx4/FH/7wB9G9e3exdOlSRfa9994Trq6uokuXLmLbtm1CCCEiIiJE586dRadOncSGDRsMnn/37l3Rr18/0bVrV7Fjxw4hhBAXLlwQvXv3Fn379hVbt241atf48ePF7373O9GuXTvFhp07d4pOnTqJLl26iLCwsArXVVK11CiHzszMFEOGDClW7ufnJ+7fvy/y8vLEyJEjq8k6icQy0Ahh2cEZqampfPrpp5w4cQJ/f3+srKxISkpCCMGECROoW7cuy5cvR6PR4O/vT/fu3avbZImk2rB4h5ZIair/+Z//yTfffENmZiZ+fn7KsmNsbCzbt2+noKCAFStW4OrqSmBgINbW1nh4eJQ6CVkmlWneY2NjRY8ePURAQICIjY1VyhMTE8WkSZPE1KlTxf/+7/+KjIwM4ePjI/z9/cXMmTONykgktZW0tDQxYcIE5fOwYcNEYWGhSEhIEJ999pmIj48Xa9asEUIIMWrUKJGbm2uyrkrNcms0Guzt7Xn69GmxGeaVK1fSpEkTcnNzadq0KUlJSfzHf/wH8+bNY9iwYUZl1GisaUY96lfGzAoxZOYgs+mSWBa3b9/m0KFDpcrMnd6AnKeGndp7me/i6uqKt7c33t7exb5btGgRkydPVj4LIdBqtbi7u5OcnExycjLNmzcHoHHjxqSmpuLq6mpSHSrl0D179qR3797cv3+f4OBgvvrqKwB++uknTp8+TXJyMpGRkcyZM4dz587x3nvv0blzZ6Myn332mVEd9ajPK5p2lTGzQnzxxRdm0yWxLIKDg8uUyXpaSES44R6AT1e4GvzfEUIQGhrKe++9R4cOHZRyrVaLTqfjzp07uLm54ebmxqVLlwB48OABDRs2NLkOlXJo/W4jFxcXcnNzlXIPDw/q16+Pi4sLjx8/5ptvvmHIkCGMGTOGiRMnkpaWZiBTkpiYGGJiYsghuzImSiRVik7oyBeGewyMsXbtWk6ePElGRgbXr1/n+++/Z9euXUycOBF/f3/y8/OJiIjA1dWVPXv2MH36dNq1a4etra3J9lVqUuzQoUPExMTw6NEjgoKCiI6OZteuXcTHx7Njxw7y8vKYN28eTk5OTJkyhcaNG5OTk8PWrVv57rvvism0bt3aqA53zStmbaFP6A6YTZfEsggODi6zhxY00Z5F4c4G5YtXDrOI3l2lWmgfHx98fHyUz3369AGgV69e9OrVq5js119/XeyzMRmJxNLRAflCV91mqGLxWz8lEktChyAfy13plQ4tkVSAAqHhqdBUtxmqGHXo/v374+LiUqxMCIGbmxurVq0yi2ESiSVS1OWuYQ7t7+/P8OHDDcr37dv33A2SSCwZndCQb8F5QYw6tN6ZMzIylOWoxo0bG3VyieRFohANeVhVtxmqqI6hx4wZQ7169bC1tUWj0bBmzRpz2iWRWCSFaHmqs9ypJ1XL3NzcWLJkiTltkUgsnkI05NfEFrphw4b4+voqeaz02SwkkhcZndCSL8rv0Ddv3mTx4sVkZGQU24sRHR3NuXPnSEpK4vXXX2fp0qW0atWK/v3706FDByZOnGiSfaoOnZCQwNatW6lTp45JD5ZIaiM6NORVIKapZAotPX5+fvj5+TFz5kzGjRsHgL29PTk5OUqghimoWubs7ExiYqLSQnt6epqsRCKpLRRWsIUujadPn3Lr1i1l2/OFCxcQQjBw4EDee+89k56pOv/u4uJCTEwMBw4cMNi2KZG8qBSi5amwNbhu375NcHAwMTEx5X7W119/XWzrtFarxcrKirp16ypZWSuK0RY6OTmZBQsWGC23lMyaEkl1oDaGbtmypdHgDH0KrQsXLrB06VISExPZtWsXAAcOHGDv3r0AXLt2jYiICKAoJsLUAwuMOnRAQAADBw4sViaE4PLly2zcuNEkRRJJbaAQTYW63A0bNiQqKsrod8+maW7dunWVTDwbdejZs2cbFe7Ro0exz3FxcYSFhdGmTRtGjBihRFtdvXqVdevWYWVlRWBgIG5ubowfP54GDRrg4OCgvMliYmKYPHky169fr3RFJBJzoBPaCk2KmRujlvXu3btcN1cmBdGjR4+Ii4tTzlOSSGoCFW2hzU2lNqX27NmTb7/9loiIiGJj7p9++ok//elPTJgwgcjISDw8PJQURPrZ8sWLF6seWyqRWCqFQkuuzsbgshRUHTo3N5cff/yR+Ph4oycRQsVTEH377bfcvXuXlJQUrl+/zmeffcalS5eUM4qfJSYmhuDgYJmCSGJR6CiaFCt5WQqqg4EJEybQtWtXbGyK3j7Gsos8m4JoypQpjB49ml27djFjxgwCAgIMUhD9/e9/Jz8/n5dffpnDhw8DRTPn+gPZnkWfPfHgqmMG30kk1UVRtJXlOHBJVB26Y8eOTJkypdSbK5OCqKxyicQS0QkN+boa6NDbtm3j5MmT2NnZodFo2L9/vzntkkgskkKsyBWWM2YuiapDX7582Zx2SCQ1gkKhIV9nOPVkKW22qkNv3ryZEydOoNFo6Nu3L4GBgea0SyKxSHQqy1YW79CXL19WutkzZswwlz0SiUWjExoKKjCGVgufDA8P5+rVq7i4uDB//vwqO7BOddkqIyOD06dPEx8fT1pamkkPl0hqG/qcYiUvNfThkyWxtrbG1tYWGxsbnJ2dOXPmDG3btmX9+vVcvHiRvLw8k+xTbaEjIyP5y1/+ovwtkUiKNpbkGUlBpI+2MnZYnTHmzp2LVqvl6NGjbNmyhUaNGlXJgXVGXy379+9XokCAYn9LJC8y4rcud8lLH21VHmeGf2/Katy4MVlZWbi5uZGcnAxU7sA6oy10q1ataNKkifJZo7HcPMQSiTnRoaGgAml81cInlyxZQlJSEg8fPmTNmjU0bdq0Sg6sM+rQb7zxBkuWLGHu3LkALFiwQJ5DJZGgnxQrv0OrhU/qfetZNmzYUCnbQMWhhw4dSmJiIhcvXgSo1PGWEkltolBoyKtpO8UOHDjA999/T/fu3c1tj0Ri0YgKttDmRtWy48ePK38vXrzYLMZIJJaOvstd8rIUVJet0tPTlb9TU1PNYoxEYukINBQanRQrNLstxlB16BYtWjB27Fi0Wq1M4SuR/IZOaCg02iJbuEMHBweTkpKCi4uL6qRYZXKKhYSEkJmZyaNHj9iyZQsODg7PpYISSVWiExryC41NiuWb3RZjqHb+Z8yYQXh4OPb29kydOtWoTGk5xZycnLCxsSmWU2zz5s3K4vmKFSvYuHEjXbp0UWbTJRJLRwgNhUYuS0G1hbayssLd3R0AJycnozI9e/akd+/e3L9/n+DgYL766iugKKfY6dOnSU5OJjIykjlz5ig5xTp37qzcf+/ePc6fP8+0adMMnh0TE0NMTIxMQSSxKNS73JaBqmV16tRRus7PTpAVu9nEnGJpaWmkpKQQEhLC+vXrsbIy7MJ4e3vzxRdfUI/6la2jRFJlCECn0xhcaty8eRM/Pz+Ds62WLVtGQEAAQ4YMITk5mdu3b9O+fXsCAwM5ePCgyfapttCLFy/mxIkTCCGYNGmSURlTc4q5uLjQv39/XnnlFUJDQwkKCqJdu3YmV0IiMReigi202mF1oaGhABw+fJjY2Fh69uyJvb09T548oUWLFibbpxFCiJKFkZGR3L17F/1XVlZW9OjRg0GDBpmsyFTcNa/wisZ8zn5Cd8BsuiSWRXBwsNHjbJ6lweBuNB7/rkG55+6LtGzZUjXaytfX1yB/XlZWFp988gmbNm3C3t4ejUbDkydPGDFiBEePHjWpDkZb6OHDh5OXl6cEZeh0OiIiIqrFoSUSS0JQNI4uidrZVmpkZmYyefJkli9fXmyFx87OrlL2GXVoV1dXjh49yu7du9FqtXz00UdVsnFcIqnpCKFBlDJmLolatNW4cePIz89n8eLFDBs2DBsbG3bu3MmTJ0/4+OOPTbZPdQx97NgxJQVRUFAQf/jDH0xWIpHUGgSlToKVRC3a6tChQwZlPXv2rJRpUIpD5+TkcOfOHQCys+XSkUQCv7XQhZaz7lwSVYeePn06a9euBSAsLMxsBkkkFo2gQl1uc6Pq0KdOnWLFihXmtMUi6K8dWt0mPFfkLH7lEDXVoY8ePcqNGzeUXWLLly83m1ESieWiAQva6lkSVYf+6quvlGUrI0vVEsmLiQBqYgsdHx/PsWPHKCgo4MMPP1T2dUskLzQ6y54UU93DdvPmTfbt28fBgwf529/+Zk6bJBLLRYBGZ3hZCqoOnZmZyQ8//EB8fDwajYbExERz2iWRWC46jeH1nLl+/Xq55FQd2tHRkRMnThAbG0ujRo3kOc4SCfw2hjZyqaAWbXXlyhVGjhzJyJEjuXLlClAUsDFt2jQlcONZjG1EMYbqGHrBggXleoA5ibl7Ce+XZVSWpPrQCA2aCrTIatFWq1evZv369Wg0GmbPns2nn35Kfn4+a9asISQkhKSkJOVoHICdO3fyP//zPzg6OpZ6XruqQ5eHyqQgio2NZfv27RQUFLBixQpefvnlMvVJZ5ZUO6JqxswZGRk4OzsD8PjxY1JSUhQHbtGiBcnJycUcWt+Kl4XRLveWLVv461//WubNlUlBFBUVxbZt25gzZ47R0/kkEotEv2xV4tIfVhcTE1Ouxzg5OZGRkUFmZiYODg40a9ZM8Y2kpKRi/gRF+0KGDx/OgAED+OCDD1Sfa7SFPnz4MBkZGZw6dQohBBqNxujGksqkIBJCoNVqcXd3VyryLDIFkcQS0ai00Grhk2rRVtOnT1dy9c2ePZsWLVpgY2NDcHAwderUKdY6A3z77bd4enoyZ86cUsM0jTr03r172blzJ4MHDy51U0lFUxCNGTOGiRMnkpaWhlarRafTcefOHYO3EaAEih9cdUxVv0RidvSTYuVELdqqbdu27Ny5s1jZ0qVLS31OdnY2Z8+e5eeff1aVM+rQDg4OdOzYkXnz5gEwadIkoxtLKpOCaOLEifj7+5Ofn09ERISqgRKJRSE0aKphY8m5c+fo378/586dIygoSFXOaAoigHHjxrFp0yYlp1h1jXPNnYKotiODM9QpTwoi5x59aDTQMDfA4AdJFcpYUlGEEJw6dYqDBw+SnZ3Njh07jMqprkPrx85arVbu5ZZIfkNTTTvFHj16xL/+9S/S09NxdXVVlVNdtpo0aRLjx48HYMqUKVVvoURSE6miZauKMnfuXIYPH05AQIASNGUMVYfu3LlzsaT4EonkN6qhw1renH6V2lgikbxoaHSgsYxz6YyiOobOzs7m/PnzAKSkpJjNIInE0qmR0VYzZsxg27ZtACxZssRsBkkklkx1TYqVF9Uut4ODg7LXtG7duuayRyKxbKppUqy8qDr0Sy+9RHx8PH/84x+NHiYnkbyIVGQMnZ2dzaRJk7C1taVPnz6MHDkSgOjoaM6dO0dSUhKvv/46S5cupVWrVvTv358OHTowceJEk+1TdeixY8fi4+ODTqdTPU5WInkRKW8LfejQIXx9fRk8eDDDhw9XHNrPzw8/Pz9mzpzJuHHjALC3tycnJ8dgD3dFUR1DR0dH8+qrr+Lp6cnq1asrpUQiqTWoJDgwFm31bAhkyV7u06dPuXXrFq1btwbgwoULbN26tdK+ZrSFHjp0KImJiUoMZps2bSqlRCKpNYiiibGSGIu2cnNzIzk5mfbt26PTFW/Wv/76a3x8fJTP+kCnunXrotPplM8VxahDHzhwgJSUFJo1a2bSQyWS2opa+CRGppl8fHyYMmUK33zzDYMHD1aCl6DIx/bu3QvAtWvXlAClPn36mOzMUMoYeteuXRw+fJj69evz8OFDLl++bLISiaTWIFQmxYw4dP369ZWlX0AZQwMcOXJE+bt169Zs3bq1SsxTfRXcu3ePd999l1OnTjF27FijMnFxcfTs2ZPAwEDi4uKU8qtXrzJ58mSmTZtGYmIiiYmJBAYGEhgYSKtWrYCiMfqoUaPw8fHh4sWLVVIZieR5U2PXoV966SV0Oh2LFi3i0qVLRmVKS0HUpEkTcnNzadq0KQ0aNCAqKoqLFy/SoEEDAE6fPs22bds4f/48cXFxtG/fvmprJpE8D2rqOrQ+ucGFCxeYPn26UZnypiD67LPPgKJcZTNnzgRgxIgR9O3bl4KCAg4ePFillZJInheqY2gLQdWhd+zYQUxMjBIXvWfPHgOZ8qYgAnjy5AkpKSlKlzsqKoozZ86QkpLCihUriIyMLPZsmVNMYpEIgUZnufkBVB36559/NurEz1LeFEQA+/btY+jQfx/V2rt3bz755BMyMzPx9/c3eLbMKSaxRDRqk2IWgqpD5+fnc+7cOerXrw+Ap6engYyPj0+xtTR9Xu5evXrRq1evYrL6ZAl6Zs2aZbLREkm1UVO73A4ODvz3f/83UDT5NX/+fLMZJZFYLDXRoW/dusXQoUOVXGKlpTyRSF4kiibFatgY+vjx49y/f79YmWyhJZIih9YWlE9WLdoqPDycq1ev4uLiwvz583F1dSUwMBBra2s8PDwqNRw16tBubm40b95cttASSUkq0OVWi7aytrbG1tYWGxsbnJ2dOXPmDG3btmXq1KmMHj2avLw8bG1tTTLPqEM/fPiQI0eO0LVrV2xtbfnuu+8YMGCASQokktqEpgLLVsnJyXh5eQHFo63mzp2LVqvl6NGjbNmyhUaNGilRWY0bNyY1NbXUVL2lYXTr59ixY3nllVcICQlh+vTpdOvWzaSHVycxd43vbpNIKoXK1k9j4ZP6aCugWLSVfv9G48aNycrKKib34MEDGjZsaLJ5qrPcGo2GCRMmYGtra/Lbwhw867jPHjcrj56VPA/UJsVatjIMn1SLtlqyZAlJSUk8fPiQNWvW0LRpU/bs2cP06dNp166dyd1tUHFoIQSdOnWiX79+ADg6OpqswFw868B6J5dOLalydAJNQfm63GrRVnPnzjWQLW/e7bIw2uXWaDRcuXKF7OxsrK2tefLkSZUoM5Wyus/GnFkieS7U1GirW7du8emnnyqfqype0xT0Dhtz95JBq/vsdxLJ80ZDDVyHBti0aRPHjx9HCMG7775rTpsM0DuysS50yTG0dGzJc8XCgzNUExwEBgaSlpZGeno6AQEB5rTJgNIcVW2cLMfPkudBURpfYXBZCqoO7ejoyOjRoxk9ejQuLi7mtEmVslpfUybDZIsuqRDCsh1atcvt4uKCn58fAO7u7kZl4uLiCAsLo02bNowYMUKJtrp69Srr1q3DysqKwMBAANasWQPAiRMnuHHjhoGMsWiuZ1FrpaVDSsyJRliWA5dE1aHnz5+vJCdwcHAwKlOZFEQlZUpDOq3EYhCABY+hVR160aJF/PLLLwA0b97caHBGZVIQqclUlJKz3BXtbsuxtqRCCIFGZ0HrVCVQdei0tDQ2b94MQEhIiFGZyqQgMibzLMVTENUv0/GkY0rMgUYn0BTUQId+9OgRO3bsQKPRkJqaalSmMimIjMk8y7MpiJ6Xs76IL4H+2qFlC72geM0ox7lSQkA5x9Bq4ZPLli3j1q1bPHz4kNWrV1NQUMAHH3xAly5d6N+/P0OGDDG5Dhqhj5EsQX5+vrIO7e3tjY2NjclKKoO75hVe0bx4jicxP14zmhvsxy5Jy0ZdePXl/gblnm9nG9y7a9cunJ2dlfDJffv2Ffv+8OHDZGVl0bNnT0aNGoWHhwdTp06lU6dOJtdBddnqyy+/ZODAgQwaNEjper8IPO8JuKp4vrkmCWPuXlIuyW8IATqdwVXRw+qysrLYv38/H3zwAe7u7pw5c4aoqCgWLlxYKfNUHfrOnTvK3/rJsReB59UNLzlpV1EnUYsqe568iEOSMhGGa9CaQqEcVuft7a2IqoVPZmZmEhQUxPLly3FwcFASiNjZ2VXaPNUxtJWVFWFhYWi1WoOT8yQVR20Puqn3mwvp1CXQCSgoXx5ftfDJcePGkZ+fz+LFixk2bBg2Njbs3LmTJ0+e8PHHH1fKPNUxNBRtEBFClLnp43kix9ASc1GuMbRzR157qa9B+WvvF5Z5rzlQbaEBXnvtNXPZIZHUDISAQsvNtG/6QbQSk5GTTDUYnYBCneFlIag69JgxY/j222/NaUutRzpyLUDoilrokpeFoOrQW7Zs4cGDB4wYMYLVq1eTnS0Pjass+gkmOdFUg9F3uWuaQ6empnLz5k0cHR1p2rSpEnklef7IltyCsfAut+qk2MqVK5k8eTIeHh4AygJ5bUIGZ0gqjBAIC2qRS6LaQn/++eeKMwM1Mjd3aVR3K1iafvmSsWAsvMtd6rJVbaa6naa69UtMRKdDFJTzcKtq4IV1aInEJPSz3OVALdrqypUrLF26FIA5c+bQtm1bQkNDefLkCXZ2dixbtsxk8+Q6tERSAYSuaAxd8jKG/rC6zZs3c/ToUaV89erVrF+/ni+//JK1a9dy584d8vPzWbNmDYWFhSQlJZlsn3RoiaQiiPI7tFq0VUZGBs7Ozjg5OfH48WNSUlIUuRYtWigBHaZg8V3uNz5sS8uWVTPDfvv2bVq2bFklz7I0fbJuVaOnLEbNHUpOTo5B+a+//kpwcLCSmAP+HW3Vvn37YgFOTk5OZGRkoNFocHBwoFmzZooTJyUl8cEHH5heCfECMXPmzFqrT9bN8sjKyhLjxo0TgYGBYvfu3WLUqFFCCCESEhLE6NGjxejRo0VCQoIQQojQ0FAxc+ZMERoaWimdL5RD//Wvf621+mTdJEIIUWr4pEQiqVnISTGJpBZh8ZNiauh0OsLCwsjMzKRjx47s3bsXd3d37O3tWblyJXfv3iUkJAQrKyvGjx9P3759+fzzz7l16xb5+flERUXx66+/limjTw9Tlr7t27ezf/9+WrRoweTJk/Hy8jJYWyzv+mNJXa1atWLv3r1YWVkRGhqKo6Ojwfrmnj17iI2NJTc3VzlruCyZ+vXrl0vftWvXDE5IMfW3LKkrMTGRjIwMfvzxR5YtW0aPHj2qtG4vGjW2hT5y5AjJycnY2Njg5uaGnZ0dOp2OJk2aABAdHU1oaCjbt29n8+bN5OXl8Y9//IN169bh5eXFmTNnyiVTXn1arRY7OzsKCwtp0qSJ0bXF8q4/ltQVGRlJ/fr1qV+/Pg0aNDC6vnn48GE2b97MsGHDOHToULlk1OpWUl/JE1Iq81uW1BUREUFUVBTNmjWjX79+VV63F40a20Jfu3aNbt26ERAQgK+vLwcOHECr1RIcHMzly5eVNUD9YQCpqak0atQIKDqrKzk5uVwy5dU3atQoxowZw+XLl1m2bBlDhw41WFvUrz8CquuPzZs3N9B16dIl9u7dS0xMDF999RX379/Hy8sL+Pf6pr4n4e7uTkJCAkC5ZIzVraS+cePGFTshZeXKlSb/liV1vf3225w9e5YOHTpgZWVFcnJyldbtRaPGttBubm7KqZhWVlbKf6TGjRuTlZWlrAHq1/8aNmzIw4cPgaKMpm5ubuWSKa++kp9Lri26ubkp64+ZmZlG1x/1+krqeu2117C2tlZOGVHLJmmsbqXJqNWtpL6SJ6RU5rcsqQuKYu8nTJigfF+VdXvRqLGz3E+ePGHq1KnY2dnx6quvcvbsWezs7CgoKGDjxo3cu3eP0NBQrK2tGTVqFG+99RZffPEFv/zyizLO+vXXX8uU0b/5y9K3ZcsW/vGPf5CamsqCBQto27Ytc+bMITc3lzp16rB06VKuXLnC8uXLAZg9e7ZRGWO6GjVqRGxsLNnZ2Xz++efY2dkxZcoU6tatS48ePZQx5HfffUdOTg7r168HKFNGP84sS993332nnJASFBREnz59TP4tS+oaM2YMY8eOVbrJ2dnZVVq3F40a69ASicSQGtvllkgkhkiHlkhqEdKhJZJahHTo58Szs6/bt2/n2LFjREVFcePGjXLd7+vrW2GdmZmZDBgwgE2bNpVLvuT0yd27d1m1alWF7Dt//rzBqYplERcXx4ABA4qdn2aMvXv30q9fvwo9+0Wnxq5DVyWHDx8mKSkJT09Pzp49y9y5c5XvQkJCKCgoUHaABQQE4OTkRIMGDZg/fz5//vOfSU9PJyMjg40bNzJx4kRatmyJl5cX586dQwjBtWvX8Pf35969e+Tk5BAeHs7jx4+xtrbm1VdfZfz48SxcuJCHDx/StGlTRo8eTWJiIuHh4fj7+7Njxw4ePnzI48ePWbVqFQ4ODgD8/e9/58svvwSKdk6dP3+eR48eFVu2CQ8PV5bJWrVqRa9evRg5ciSDBw/m7bffZvXq1dStW5fBgwfTrl07UlJSKCgo4NNPPyU3N5fCwkIiIyP54osvuHPnDs7Oznh7e5OYmMiiRYvo0aMHDx48IDU1lRkzZuDo6Mjrr79OQEAAXl5ejBs3jp9++ono6Gjq1aun2DVgwABatGhBeHg46enppKen8/rrr5OamkpycjK7du1ixIgRfP3112b6X1A7kC008OGHH/LPf/6TjRs38qc//UkpT0xMxNbWllWrVjFz5kyOHz9Or169iIyM5JdffiEzM5Pbt28TGRlJnz59OH78OACffPIJ/fr14969e0RERNC1a1cDncOGDSMiIoLjx4+j0+nQ6XQ4Ojpy6NAh3Nzc8PT0VBw/Pj4eJycnbGxsuHr1qvKMtWvXsmnTJjZt2sT69esZNGgQ3bp1Y8CAAcV0+fj48Oc//5ljx44B4OnpSWhoKPv27WPhwoVs2bKF3bt3K/InTpzg9u3bODs7k5WVRUpKCgkJCaxdu5aFCxfSo0cPPD09mTdvnnLPX/7yFyZMmMD69es5efIkULSmPGvWLLp27crFixdVf//hw4ezdu1afvjhB5YuXUq9evVIT0+vwL+gRI9soSnqej569AitVktBQYGy4UEIoWyq0H/Wr0uXRKPRKF1YJycnCgsLsbW1BaBOnToG8vp1UiEEFy5cQKPRsHDhQmWL5LN7yNu0aUN4eLhRu/Vypa0+FvyW1C4/P1+xr+T9z9ZLp9PRvXt3pk2bBhTtaitZ75Kfjf02+jra2NiQm5urap+joyN16tTB0dERAFtb21LlJerIFpqiHE8fffQRs2fPJiwsTClv06YNOTk5hISEsG7dOt555x1Onz7NrFmzaN68OY6Ojri7uzNr1ixiY2OLnQ3s5OSEq6srn3/+OfHx8aXq//3vf8+lS5dYuXIlDx48UHT/8Y9/xNHRUdliGhQUVCzf1JQpUwgMDCQwMJBJkyapPn///v0EBQXx4YcfFisPCAggLCyMgIAAPvroI+UF5u3tzeXLlwkJCeGTTz7BxsYGT09PZs6cqRxIbmNjU+y3+uijj4iOjmbatGn07Wt4OqPEPMiNJbWc8PBwfH19adu2bZmygYGBDBgwgPfff/+52hQXF8eVK1eYMmVKmbK+vr5yHF0BZJe7lmOsq65GVFTU8zPkGZydnbl06RJ37tyhRYsWqnJ79+6lWbNmZrGptiBbaImkFiHH0BJJLUI6tERSi5AOLZHUIv4fiqYubPbhSc4AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 250x150 with 2 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"show(\"forum626-ndsi-mask-float.nc\")"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "d7f96673-e02b-498b-bb16-da8f0f972f6f",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "7bb9a7f5-ebf2-4524-ae99-7b0f0497bada",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "195619c2-7686-465b-98a8-2e30f7dab216",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "54998979-8881-4b5f-9d84-10529d4678c4",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "openeo-client-venv3.9",
"language": "python",
"name": "openeo-3.9"
},
"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.9.15"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment