issue401.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "advised-religious",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Populating the interactive namespace from numpy and matplotlib\n"
     ]
    }
   ],
   "source": [
    "import rasterio\n",
    "from rio_tiler.io import COGReader\n",
    "%pylab inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "respected-suggestion",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{\"bounds\": [-77.92959037999877, 25.148614663322203, -77.76954571373244, 25.254492610438696], \"center\": [-77.84956804686561, 25.201553636880448, 9], \"minzoom\": 9, \"maxzoom\": 9, \"band_metadata\": [[\"1\", {}], [\"2\", {}], [\"3\", {}]], \"band_descriptions\": [[\"1\", \"\"], [\"2\", \"\"], [\"3\", \"\"]], \"dtype\": \"uint8\", \"nodata_type\": \"Mask\", \"colorinterp\": [\"red\", \"green\", \"blue\"], \"scale\": null, \"offset\": null, \"colormap\": null, \"height\": 38, \"width\": 53, \"driver\": \"GTiff\", \"count\": 3, \"overviews\": []}\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVEAAAD5CAYAAACeTgZ1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAfb0lEQVR4nO3dfXCddZUH8O9J2jSh76Wl1lJsgUJBlCKxguAsgsWqKDoqK7qKDmvdXVl11vF9RnfH3VnxBcTd9aUKax1fAEFXZVGotQ6yo2gLBVpKoWCFlrSh0DckTZrk7B/36UzlOSe5J7/fvbm5/X5mMk1+efI8v+feJ6c39zzn/ERVQUREI9My2hMgIhrLGESJiBIwiBIRJWAQJSJKwCBKRJSAQZSIKMG4lB8WkWUArgHQCuBbqvq5obZva2nXjpbJ5W9Yt1m1OPHduyUrx61a3j5E0vcdFTkfb36NdD4R1rwbfc5ROc4nes0bx9TBQW/noeGIk848Pn0ndbZu3bpdqjrL+p6M9D5REWkF8BCApQC2AfgDgEtV9QHvZ6aOm6VnT7m4NK4H+8v7P+ooeycH+8xhax8uL0AfPGiPt7ZWvx/vovSO6W0/MGCPW7z5efsYP776fbu/ZIboOXqseUeeg6hcz1ngepDxSa9fAAxxzTvzltby+GDPAXsfzrWjg4F4MWjvY9XgD6vfR4MQkXWq2ml9L+UKXAJgi6o+qqp9AK4HUI6QRERNLCWIzgXw+GFfbyvG/oKILBeRtSKytk97Eg5HRNR4ap5YUtUVqtqpqp1t0lHrwxER1VVKEN0OYN5hXx9bjBERHTFS3t3+A4CFIrIAleD5NgBvH/InVKtOAOkB5w3vaJIiwkteeKy5RBJCOedi8RJI1rydZIT228+XjAtcOsFEnrRPqHrXoYSiN5ccySkg9rg68w4lnIKJrxythqTFvqvAu04sS8e9zRxf1X/9iOY02kYcRFW1X0SuAHAbKrc4XaeqG7PNjIhoDEi6z0JVbwVwa6a5EBGNOaxYIiJKwCBKRJQgvWwiBysZ471pnqMqxqtM8ngJGms/ORJCQCzZ4Z175DHxKmva2mq2b+9x1QFj+1xVRRl4yR9z3kHWPqxKoyEFrh3v+dUe+57uUMWSI5pAshJRjZSE4itRIqIEDKJERAkYRImIEjCIEhElYBAlIkpQ/+x8tZnDXBl0K0sbybYPJVcm3hLILrvZ4mC/yZBcpZJVCp9jRDDznyMLHzkfLx/uld66JZjG9R3tJ+yWfcL4XXD6iUY1UibewleiREQJGESJiBIwiBIRJWAQJSJKwCBKRJSgcWvng5lvr77YzD162VjvmJEVOaOrUmZYOdLNFkdXNTVIh72kS6QJr7tvryZcyuPaZ6/yGl7R1Dr3yOqnQbWse3cfE28f1rXpPI9ujbyTcbfuFNAa9k9vJHwlSkSUgEGUiCgBgygRUQIGUSKiBAyiREQJkrLzIrIVwH4AAwD6VbVz2B+yMsxW1jC6HOyB3mEPPeTxhpKhTjx09wCQZzlmb96BbHSOLHx02eDQssbR58Y692jX/IAcdfZhzt0X0Tp5U4v9u5PjOhmrSynnuMXplaq6K8N+iIjGHP45T0SUIDWIKoDbRWSdiCy3NhCR5SKyVkTW9umBxMMRETWW1D/nz1XV7SJyDIBVIvKgqt5x+AaqugLACgCY2jozw5syRESNIymIqur24t9uEfkxgCUA7hj6pwzRhIklR4NkZx5y1FHmuD77bNW7DicYIsmVGi4P7LKO6czPLX30km2RRss1bA7tJgMjz2UwaWUd0z2e93vj/C5Y7ZS1107Ius2XnQTSqsEf2nMJaPQEkmfEV6CITBSRyYc+B3AhgA25JkZENBakvBKdDeDHInJoP99X1V9kmRUR0Rgx4iCqqo8COD3jXIiIxhze4kRElIBBlIgoQWMsmRxoQJwlGxvMZkey8OH5ReaSqfQx0ig4y7LLRpPlys4zPMfROxPqvNSzyytfDtz1EF7628jau8sueyWiXrPmI1iDXFFERGMTgygRUQIGUSKiBAyiREQJGESJiBLUPzsfyT7W6niZSFtbaUx7euyNczSC9rKuwSV/zVrsHHdDeBnnYMPeUP14hjsTzIz4UMeMPFY5lsqOcq41r04+B7OhsrO8sldnv7TlraHtGwVfiRIRJWAQJSJKwCBKRJSAQZSIKAGDKBFRgrpm51XVzNS2zJ9X3vaJnfZOcmQ1c3TSB6B9feXBYKbcY2aREdx3vbPIme6yyNI53mP1Dcg071BPAu8bkbtXAjXy7jyitfBOxt0Szao3ehbew1eiREQJGESJiBIwiBIRJWAQJSJKMGxiSUSuA3ARgG5VPa0YmwHgBgDzAWwFcImq7h52X5WfLY/3ld8gD7d+jSRAooklb3sriRRNfEWa80ZlSJjI+Opzj1nKJD2Zmi+bTaajj5O37wNGWaWX5Ik02PbOPdqU2RJIFAFjN/lTS9VcPd8GsOw5Yx8HsFpVFwJYXXxNRHTEGTaIquodAJ5+zvDFAFYWn68E8Ma80yIiGhtGep/obFXtKj7fgcoa9CYRWQ5gOQC0y8QRHo6IqDElv2mmlRWt3LcwVXWFqnaqamcbJqQejoiooYw0iO4UkTkAUPzbnW9KRERjx0j/nP8pgMsAfK749ydV/ZSIma0c7N5V3ja6PHAke+tkTL1MdDjrHBE5n8hdAkMJlBa6SyYHRDL8gPN4R5/3GjbjdkXu1oiInrt3R4DVHLvF2dbJ2o/Vxsm1NOwrURH5AYDfAjhZRLaJyOWoBM+lIvIwgFcVXxMRHXGGfYmgqpc637og81yIiMYcViwRESVgECUiSsAgSkSUoP5LJlty1FBnyNqHmgFH5aj9zrX8buTxzlB/n6Wm3qkHl44O+5iBZZq9ZsoNcz1En0fWztcVX4kSESVgECUiSsAgSkSUgEGUiCgBgygRUYL6ZudVsy1XXJJh+diwSBdy77y9On5ryWSvjj3DnQnS1maOm8tCe8fM1H3ezJa32h3A3PkFjpmjPwBg9wgILY0MZLlTpdJYzRi3lkcO1s5TGV+JEhElYBAlIkrAIEpElIBBlIgoQX0TS05T5ixJisjStJGyuFyCS+eGGhNnoD099jechs9m4svbeYbGxDUtwazh45rlmNEkVIbkrYyzQ4PXlDmi2UpH+UqUiCgBgygRUQIGUSKiBAyiREQJGESJiBIMm50XkesAXASgW1VPK8b+GcB7ATxZbPZJVb11xLMYjeVtIyJLFeco5/P2kyuLbJ1PIAsP1DhbHhF8TLzzifDKRN27Eyw1vOa9RtWIlshavGbXXvnoEaCaK+rbAJYZ41er6uLiY+QBlIhoDBs2iKrqHQCersNciIjGnJS/ba4QkftE5DoRme5tJCLLRWStiKzt0wMJhyMiajwjDaJfA3ACgMUAugB8ydtQVVeoaqeqdrZJ+wgPR0TUmEYURFV1p6oOqOoggG8CWJJ3WkREY8OIaudFZI6qdhVfvgnAhnxTKgSbGLusLGh0H07mOpRhzZW1D7CaBAMArObBXsY5w7zdeTiyZP69hsWRvgne857jOQs0pHabRgefG68e3uItOe3tw2z47DR29urvx2pNfTW3OP0AwHkAZorINgCfAXCeiCxG5a6OrQDeV7spEhE1rmGDqKpeagxfW4O5EBGNOaxYIiJKwCBKRJSAQZSIKEF9O9t7rEx8NIOeYdngsBwd+SOC2djQEss5Ms5eRjx9z3He+dSyi33kuXfuPtHBKld+GOp43h0I1t0GzrZuFt7eMzBoXGvBevql497m7Nt4rIJLPdcy889XokRECRhEiYgSMIgSESVgECUiStAYiaVoEskSeVN/NJJQngZpSB1uvhxoGh3dd6Rxsh7otffRPqH6Y+Z6DiIJUq+k1JJr+XBrPzlKmuGU9gZLt0XEHFetft+399e/dJSvRImIEjCIEhElYBAlIkrAIEpElIBBlIgoQWNk5wPlk9FmsSavLC6S1fRkWsI30pjY3Udg3tFGyFY2NstSwsG5uFl4r+TVkutuDeP68RpShx7vHMttZ+L+nhnZ8kGrUTNgl4jCv04iTb0vbLM6dzpNowMmY/qZ3vf4SpSIKAGDKBFRAgZRIqIEDKJERAkYRImIElSz2uc8AN8BMBuVBNoKVb1GRGYAuAHAfFRW/LxEVXePaBZWNtGpjfXqa0MZVi8LX8uGvY4cWWS3Bj1HZjhiFHoSuDnXyHOZazlrYz+h3gOI3fUQFmkiHl2y3LwzwX52vPNxs/DeXKx9u3cEVL+PlokTS2PyrH8tVHOV9AP4sKqeCuAsAO8XkVMBfBzAalVdCGB18TUR0RFl2CCqql2qenfx+X4AmwDMBXAxgJXFZisBvLFGcyQialihm+1FZD6AMwDcBWC2qnYV39qByp/71s8sB7AcANql/DKZiGgsq/pNHxGZBOBmAB9S1X2Hf09VFc7bU6q6QlU7VbWzTdqTJktE1GiqeiUqIuNRCaDfU9UfFcM7RWSOqnaJyBwA3cPuSLX6xrWRprVDqeWKnDUsB80hlEAKPiahIrocZYu1LH0MJC7CcwkmbrRWK7F6ovsOJJykrc3c1G++XP1V5ZWftnR0mOODzz5r78g4n5bJk8rb9SYklqRyxtcC2KSqVx32rZ8CuKz4/DIAPxluX0REzaaaV6LnAHgngPtFZH0x9kkAnwNwo4hcDuBPAC6pyQyJiBrYsEFUVe8E4NyciQvyToeIaGxhxRIRUQIGUSKiBPVtyixiZ+KNDKH29dm7iDRfdvYd5TY9tga9cr7IPjzRuwoi517LxylY+hjOlltG4W6I0J0g3rLBxmPoln1magBuCZfTGufpHm+cnUFvmTbFnotVDuo8rvtfdIw5PunRfeZ4y75y1n6we1d5wyHudOErUSKiBAyiREQJGESJiBIwiBIRJWAQJSJKUN/svFc7b4y5zZcziDYrDmeXrX1EG9Fa+6hlzb+TdY1k3MPNnr1MrzXuLXOdq8dCtfPAEI+J9RxneFxrugx3oDk0ELxTwLmTpufMBeb4zk671n7C3vLY3lPsOzgW/MieX89cox4ewFNLp5fGnn9H+S4B3WAvzQ3wlSgRURIGUSKiBAyiREQJGESJiBIwiBIRJah/7byVTa3hkrqW8BK0kSxydC61Wr4YCD2uoYwzUNvadGveXhY+eu1Yd4cEM/yRxyTaTyBSO+91jveIVa4fqKcf6pgydXJprOs1c81tBybYd96IU7C/7/jyY3XCoifMbXecOc8cb3m5vZr79I4DpbHHO+aUxvoeS1symYiIHAyiREQJGESJiBIwiBIRJRg2sSQi8wB8B8BsVHq1rlDVa0TknwG8F8CTxaafVNVbh9yZql2+V2Wj5iF5CYYMy+/mKKOLLDULoO4NlRsmgZRLDZOBkX1kSRxG56zB8mWDd833LzrOHH/yzImlsfbX7TS33fHYDHP81ldfY46f0nZUaewln/17c9uPfPBGc/yA2gmxzc8+rzT2PydNK41pu//YVZOd7wfwYVW9W0QmA1gnIquK712tql+sYh9ERE2pmtU+uwB0FZ/vF5FNAOx7F4iIjjChvxNEZD6AMwDcVQxdISL3ich1IlJuh1L5meUislZE1vahN222REQNpuogKiKTANwM4EOqug/A1wCcAGAxKq9Uv2T9nKquUNVOVe1sg99OiohoLKoqiIrIeFQC6PdU9UcAoKo7VXVAVQcBfBPAktpNk4ioMVWTnRcA1wLYpKpXHTY+p3i/FADeBGDDsEcLLJnsZm5rmXV1hJoyj0IWPrIUrsdt+Bx9HiIi+4jOo4bzDi09XMOyXu0Pli8bj4kumm9u+tSL7OWLnzzbPuY7Xvab0tjqrpPMbaXXfvw+s+315viZUx8rje0+027S/a4pxnLHAE773TvM8beecE9prGNjeUnnlh7/Oa8mO38OgHcCuF9E1hdjnwRwqYgsRuW2p60A3lfFvoiImko12fk7AVgdA4a+J5SI6AjAiiUiogQMokRECRhEiYgS1Lcpcw45srS1rAev5RK+zrl7mfXQ0tA5lvatpWi/g8BS1J7wEtqWYH8Apy9xaN8DLz7R3nxC+Q6R3un29do3xW6cfPTv7cf15+vOLY1NfPMOc9uZJzxtjl9z3E/M8Qu+8dHS2LgpoUcK15x+vTm+o39aaewFN2wvjW1/2lmyG3wlSkSUhEGUiCgBgygRUQIGUSKiBAyiREQJGiM7H8lgerXptZyHl7k2MsChOvuhWMcM1omH5uLsO0smOsMqBaE7DQDoAbvtomo5q9t73ovMbTvu/pM53vtiu7v7hG17y4M7nyyPRYlz/U0sd3wHgD0LyrXfADB18/7S2ON/b9fCt2yc5E3G3r6v/Lj++eZy13gAOOXdm8zxv3no7eb4rPvKc5x0TzmDDqBSoG64oMOOG2955IzS2NHfLy+vPO7dfp8CvhIlIkrAIEpElIBBlIgoAYMoEVGCxkgsWckir4mxV1YZSWrUcqleT7R5sMEtwfR+oN7nmauc1thPtMxU2u2laFomlRMmHev+6OzEPp8J95WbBAMADvZVvQ+P9vSUxvZfdLq57fj9drJkYpcxDwAtPeXSxY7/m2puu+9kb9/2+Yx/pnwVqnPq635xqjm+6e++ao6f8Mq/K409r2Oeue0/bD/LHP/lI3aD6PtfcW1p7JxPf6A01rPrdvPnAb4SJSJKwiBKRJSAQZSIKAGDKBFRAgZRIqIE1SyZ3A7gDgATiu1vUtXPiMgCANcDOBrAOgDvVFU7LTicapdRBuKlmUZGO9bO1T9mjga6ke3DWfha3hGQoylzYN4yzr5UvSy8d3fH4N595X1kWHLadfQ0e3yX3ZjYMuXXW8zxk24zykwBtMJ+XB98Szmj/fzbus1tn/9L+zF56D1Hm+MDU8rZ/EVffcbctvtCuyz1w10vMcenbyyXmj5xod0k+a1H2Y2gf9sx3xz/2I6zS2O7X1j+TRv4hfnjAKp7JdoL4HxVPR3AYgDLROQsAFcCuFpVTwSwG8DlVeyLiKipDBtEteLQfynjiw8FcD6Am4rxlQDeWIsJEhE1sqr+jhGRVhFZD6AbwCoAjwDYo6qHWptsAzDX+dnlIrJWRNb26YEMUyYiahxVBVFVHVDVxQCOBbAEwKJqD6CqK1S1U1U726R9ZLMkImpQoXfUVXUPgDUAzgYwTUQOvdt/LACnwR8RUfOqJjs/C8BBVd0jIh0AlqKSVFoD4C2oZOgvA2Cvd1qFyPK20cx6lqbCkUbQXkbc20ekF0B0fhlq2XM8ftrvNP49cb45fnBGudnwuL3lmnIAwOM77WMaNeiAnc2XyZND+4D3mKgx/tQee9sI53ib32MvjbznhdPM8RnoKo3tOnuWuW3vVLv5cpt9QwDmf36r/Q3DxA0LzfEvLb3bHH/X35b7Hdz5h1PMbb+1+Rxz/Nh/t38XNkx6cWlswrnG3SFDhIBqotccACtFpBWVV643quotIvIAgOtF5F8B3AOgXMlPRNTkhg2iqnofgFIPfVV9FJX3R4mIjlisWCIiSsAgSkSUgEGUiChB/TvbW13LDxrZWyez3OIsEzu4367TrXYOQ4p20zdYS/UC3gK0iHXk91YByLHvDBn+lg77/uDBo9qq3sdgu/0ctM6xs8vS5SxVbGTQveWV937X7vr+xNaZ5vipV5bvFOg5wd52wl0P2fMzrjXv7obNl08zx3WKXVc+Y0252KVnpn0FHnezfcfiM6fNNscx1Vhi2ekPMPeXdor/7O3lDvYA0DulPMdFa+znd89i+/Fu2bLZHG8zrof5T5Svqa7dXDKZiKgmGESJiBIwiBIRJWAQJSJKIF7Coxamts7Usya9oW7HA2AnTKLJkkZZjjnXksQRkVJTJ8ElbXYCyWuoPPC8cuPfF137gLntxgtnmOObv3ycOX78nF2lseMm7ja33faBBeb4k2caSRQA/R3lBMjH3neDue33ziqXG7qchtR6wO6Ktve1LzTHWwbKv+uDrXZiadrqh+25eCWvVmNrZ7loGWdfJ16Cz7pOdPoUc9vBDjsB2X5V+XkHgKevekFp7KjH/lwa+92DK7Dvz0+YDxZfiRIRJWAQJSJKwCBKRJSAQZSIKAGDKBFRgvqXfdaIu7SvNZgrg27sx2sw7WUevZJSc6lnqzwW8LP2B+3yv8gS1V5m3dLzMrtJ8NMn2+d4zN120+PWA+Xz/MisO81tp99rl5R+fc/x5viKh8pNe1ed8jNz2zOW2Bn0PruHMzb+41ftbxh+doudFX/9zHtLY1/8yl+b2z7/5kfM8Ynb7az9ziXlkuk/v8R+DvYsPNkcn/+fG81x9FdfTjvo3BHUeoxdwnvgpHKp6eB4+3rdu8D+/RvYbf/Od0wq72fHJeXMf99/+KXVfCVKRJSAQZSIKAGDKBFRAgZRIqIEwwZREWkXkd+LyL0islFE/qUY/7aI/FFE1hcfi2s+WyKiBlNNdr4XwPmq+oyIjAdwp4j8vPjeR1T1ptARU2vZvWV5I5noKG9JYmPf7hLDwcbOGsj8exTOMSP7cBoCW5n/9t/Y9e1TOk4zx9+64jZz/PP/W+6v8OhBOwv/2b+62Bzf+k67dn7umnJd9IW4zNx28GN2TX3//dPM8YW/fnd52/32c3DlK280xz+19k2lsQkX2E2M9z0x3xxX55If11POiutTdv+C476wzhzf9BX7joXWyeXr4eR/Ki/RDAB69DRzvPtldh+E1t7yvNveZS+V3Tbg9G/4rt2seeeyvtLYhC3GtTbEDT3VrPapAA61jR9ffNSvawkRUQOr6iWgiLSKyHoA3QBWqepdxbf+TUTuE5GrRcT8L01ElovIWhFZ26f2/WtERGNVVUFUVQdUdTGAYwEsEZHTAHwCwCIALwUwA8DHnJ9doaqdqtrZJvafZEREY1UoO6+qewCsAbBMVbu0ohfAfwNYUoP5ERE1tGqy87NEZFrxeQeApQAeFJE5xZgAeCOADbWbJhFRY6om3TsHwEoRaUUl6N6oqreIyK9EZBYqq/6uB2Cvd/pcRiberBP3stxeJn80ur7XUrVLSzvbAn4/AYu3b+11av6NOxM2X32quekfL/qGOX7jM/aSxFve/nXrgOa2r72tXGsOAFets+uwb79iZWnskkcvsLc9frU5/ovT7Iz2lVe8qzT2+Kvsea9435vN8TlTyr+Su0+ylwkfaLN/RyY9ZtfDT364nIlu6Z9mbvvUO15ijnc8btf8T9hdvguh/7hjzG33nTjRHG9xbrDZ+Vflu2Oky87kn7bAXup542uchgdaPp/emeXj6RCRsprs/H0AzjDGzx/uZ4mIml2TvXwjIqovBlEiogQMokRECRqiKbOZRIosUxzkNnCONj3OcczkPQ+RQHKWrIUGSkqdktf+l55SGjv104+Z237h5SeY4x+Y/qA5/uXdC0tj33jgFea202+ykxRnXLHFHLfc6CSQrtj+MnP89kfshsVzjbEpW+xEzLhn7CzKhK1PlcZae+0ETc/RdtKq9dlyAgkA9p1cTuR5yZzxz9pX5jHftRM36C0fc+dF9pLTT51lH3R8t10i29JR/r2cN9suyW1rsX+HX71okzn+W6N0du8+I+FpP42V+fnfIiKi4TCIEhElYBAlIkrAIEpElIBBlIgogaizfGlNDibyJIA/FV/OBLCrbgcfHTzH5nAknCNwZJznSM/xBapq1hLXNYj+xYFF1qpq56gcvE54js3hSDhH4Mg4z1qcI/+cJyJKwCBKRJRgNIPoilE8dr3wHJvDkXCOwJFxntnPcdTeEyUiagb8c56IKAGDKBFRgroHURFZJiKbRWSLiHy83sevFRG5TkS6RWTDYWMzRGSViDxc/Dt9NOeYSkTmicgaEXlARDaKyAeL8aY5TxFpF5Hfi8i9xTn+SzG+QETuKq7bG0SkbbTnmqpYCv0eEbml+LqpzlFEtorI/SKyXkTWFmPZr9W6BtFinab/AvAaAKcCuFRE7IV5xp5vA1j2nLGPA1itqgsBrC6+Hsv6AXxYVU8FcBaA9xfPXzOdZy+A81X1dACLASwTkbMAXAngalU9EcBuAJeP3hSz+SCAw3vENeM5vlJVFx92b2j2a7Xer0SXANiiqo+qah+A6wFcXOc51ISq3gHg6ecMXwzg0MpoK1FZFXXMKpbJvrv4fD8qv4Bz0UTnWSwD/kzx5fjiQwGcD+CmYnxMnyMAiMixAF4H4FvF14ImO0dH9mu13kF0LoDHD/t6G+xets1itqp2FZ/vADB7NCeTk4jMR2UBw7vQZOdZ/Jm7HkA3gFUAHgGwR1UPdfxthuv2ywA+CuBQh+6j0XznqABuF5F1IrK8GMt+rTZEZ/sjgaqqiDTF/WQiMgnAzQA+pKr7Ki9iKprhPFV1AMBiEZkG4McAFo3ujPISkYsAdKvqOhE5b5SnU0vnqup2ETkGwCoR+YulFHJdq/V+JbodwLzDvj62GGtWO0VkDgAU/3aP8nySich4VALo91T1R8Vw050nAKjqHgBrAJwNYJqIHHrRMdav23MAvEFEtqLyltr5AK5Bc50jVHV78W83Kv8ZLkENrtV6B9E/AFhYZAHbALwNwE/rPId6+imAy4rPLwPwk1GcS7LifbNrAWxS1asO+1bTnKeIzCpegUJEOgAsReW93zUA3lJsNqbPUVU/oarHqup8VH4Hf6Wq70ATnaOITBSRyYc+B3AhgA2oxbWqqnX9APBaAA+h8j7Tp+p9/Bqe1w8AdAE4iMr7SZej8j7TagAPA/glgBmjPc/EczwXlfeZ7gOwvvh4bTOdJ4AXA7inOMcNAD5djB8P4PcAtgD4IYAJoz3XTOd7HoBbmu0ci3O5t/jYeCjW1OJaZdknEVECViwRESVgECUiSsAgSkSUgEGUiCgBgygRUQIGUSKiBAyiREQJ/h+bh87b5wZp3AAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "with COGReader(\"masked.tif\") as cog:\n",
    "    print(cog.info().json())\n",
    "    imshow(cog.preview().as_masked()[0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "fuzzy-mounting",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAASsAAAD4CAYAAAC0ecCBAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAMX0lEQVR4nO3dX4xc9XnG8efZXW+9YARJoFVrW7UvEImVNjjdolCqSIUSkSYiF20lkIjUqJUvmj9QRYpIe9X7KiIXUSQLSCuFgBonSFFEHahCFEVqaWzjtGAblbpJsENiUEX5kyYbe59e7EB3vbves/Wc39l3+X6kFTs7o/k9Z3d4fM6cmXmdRACw0U0MHQAAuqCsAJRAWQEogbICUAJlBaCEqT7udHpia2YmtvVx18vNNzybOeF2a0lSyxO1jTetmU38O7z62l1tF2zg8OHDLya5aqXreimrmYltun7bh/q462UyN9dkHUny9HSztSSp5ctK7M3ZVpv5d3jw0H1N12vB9g9Wu47DQAAlUFYASqCsAJRAWQEogbICUAJlBaAEygpACZQVgBIoKwAldCor27fYfsb2s7bv7jsUAJxvzbKyPSnpc5LeL2mPpNtt7+k7GAAs1mXP6jpJzyY5mWRO0kOS2rzxDwBGupTVdknPLbp8avSzJWzvs33I9qG5+Z+NKx8ASBrjE+xJ9ieZTTI7PbF1XHcLAJK6ldVpSTsXXd4x+hkANNOlrL4r6Wrbu21PS7pN0tf6jQUAS6354XtJztr+mKRvSJqUdH+Sp3tPBgCLdPqk0CSPSHqk5ywAsCpewQ6gBMoKQAmUFYASKCsAJVBWAEqgrACUQFkBKKGXicxK20m4rbTeJk/18+cZWs6eHTpCLzbjY34jYc8KQAmUFYASKCsAJVBWAEqgrACUQFkBKIGyAlACZQWgBMoKQAmUFYASukxkvt/2GdtPtQgEACvpsmf1t5Ju6TkHAFzQmmWV5NuS/qtBFgBY1dje1m97n6R9krTVl47rbgFAUl/j4834eADjxdlAACVQVgBK6PLShQcl/ZOka2yfsv2n/ccCgKXWfII9ye0tggDAhXAYCKAEygpACZQVgBIoKwAlUFYASqCsAJRAWQEooZ/55JZs93LXy8zMtFlHA4w9P3eu3VqTk82WavbYaGz+pz8dOsKmxp4VgBIoKwAlUFYASqCsAJRAWQEogbICUAJlBaAEygpACZQVgBIoKwAldPkM9p22H7d9zPbTtu9sEQwAFuvy3sCzkj6Z5IjtyyQdtv1YkmM9ZwOAN3QZH/98kiOj71+RdFzS9r6DAcBi6/rUBdu7JO2V9MQK1zE+HkBvOj/BbnubpK9IuivJy+dfv2R8/ATj4wGMV6eysr1FC0X1QJKv9hsJAJbrcjbQku6TdDzJZ/qPBADLddmzukHShyXdaPvo6OsPes4FAEt0GR//HUmb83NoAZTBK9gBlEBZASiBsgJQAmUFoATKCkAJlBWAEigrACVQVgBKWNenLnQ2H2Vurpe7XmbLljbrDCBJs7Vavuq32WNDkqen2601M9NsLUm6eeKPm6312PyXm621GvasAJRAWQEogbICUAJlBaAEygpACZQVgBIoKwAlUFYASqCsAJTQZWDEVtv/Yvt7o/Hxf90iGAAs1uXtNj+XdGOSV0cjub5j+x+S/HPP2QDgDV0GRkTSq6OLW0Zf7d60BgDqPuR00vZRSWckPZZkxfHxtg/ZPjSnn485JoA3u05lleRckmsl7ZB0ne13rnCb/xsfr18ac0wAb3brOhuY5CVJj0u6pZc0ALCKLmcDr7J9xej7GUk3SzrRcy4AWKLL2cBflfR3tie1UG5/n+Tr/cYCgKW6nA38V0l7G2QBgFXxCnYAJVBWAEqgrACUQFkBKIGyAlACZQWgBMoKQAmUFYAS+hkfP+GmY7tbaTnOvblz54ZO0IvNOqpekia2bm221vtm7miyzmV+62+tdh17VgBKoKwAlEBZASiBsgJQAmUFoATKCkAJlBWAEigrACVQVgBKoKwAlNC5rEaDTp+0zbAIAM2tZ8/qTknH+woCABfSdXz8DkkfkHRvv3EAYGVd96zukfQpSfOr3cD2PtuHbB+am//ZOLIBwBu6TGT+oKQzSQ5f6HZJ9ieZTTI7PdHuoysAvDl02bO6QdKttr8v6SFJN9r+Yq+pAOA8a5ZVkk8n2ZFkl6TbJH0zSZtP4gKAEV5nBaCEdX2scZJvSfpWL0kA4ALYswJQAmUFoATKCkAJlBWAEigrACVQVgBKoKwAlNDP+Pi0G7Vuu8k6rdeSNu+4+tZj1ltp/ff69/ve0Wyt39x5qsk6E382ufp1TRIAwEWirACUQFkBKIGyAlACZQWgBMoKQAmUFYASKCsAJVBWAEqgrACU0OntNqPJNq9IOifpbJLZPkMBwPnW897A30vyYm9JAOACOAwEUELXsoqkR20ftr1vpRssGR8fxscDGK+uh4G/m+S07V+W9JjtE0m+vfgGSfZL2i9Jl09euTk/2wTAYDrtWSU5PfrvGUkPS7quz1AAcL41y8r2pbYve/17Se+T9FTfwQBgsS6Hgb8i6eHRp2ROSfpSkoO9pgKA86xZVklOSnpXgywAsCpeugCgBMoKQAmUFYASKCsAJVBWAEqgrACUQFkBKKGf8fFuP2q9hdbjwVv+Dltum6f6edgNzVe+pel6b73i1WZrXTL1iybrTHr1xyF7VgBKoKwAlEBZASiBsgJQAmUFoATKCkAJlBWAEigrACVQVgBKoKwAlNCprGxfYfuA7RO2j9u+vu9gALBY1zdpfVbSwSR/ZHta0iU9ZgKAZdYsK9uXS3qvpD+RpCRzkub6jQUAS3U5DNwt6QVJX7D9pO17R/MDl1gyPn6e8fEAxqtLWU1JerekzyfZK+k1SXeff6Mk+5PMJpmdntg65pgA3uy6lNUpSaeSPDG6fEAL5QUAzaxZVkl+LOk529eMfnSTpGO9pgKA83Q9G/hxSQ+MzgSelPSR/iIBwHKdyirJUUmz/UYBgNXxCnYAJVBWAEqgrACUQFkBKIGyAlACZQWgBMoKQAmUFYASur6CfcNKMnSETeGh4482W+sP7/hos7X+8Uv3N1vrN+7582ZrSdKOv3qp2VonfvsdTdb5nxcPrnode1YASqCsAJRAWQEogbICUAJlBaAEygpACZQVgBIoKwAlUFYASlizrGxfY/vooq+Xbd/VIBsAvGHNt9skeUbStZJke1LSaUkP9xsLAJZa72HgTZL+I8kP+ggDAKtZ7xuZb5P04EpX2N4naZ8kbV0+XR4ALkrnPavRzMBbJX15pesZHw+gT+s5DHy/pCNJftJXGABYzXrK6natcggIAH3rVFa2L5V0s6Sv9hsHAFbWdXz8a5Le1nMWAFgVr2AHUAJlBaAEygpACZQVgBIoKwAlUFYASqCsAJRAWQEowX2MX7f9gqT1fozMlZJeHHuYjWGzbhvbVc9G37ZfT3LVSlf0Ulb/H7YPJZkdOkcfNuu2sV31VN42DgMBlEBZAShhI5XV/qED9GizbhvbVU/Zbdswz1kBwIVspD0rAFgVZQWghA1RVrZvsf2M7Wdt3z10nnGwvdP247aP2X7a9p1DZxon25O2n7T99aGzjJPtK2wfsH3C9nHb1w+daRxs/8XocfiU7Qdtl5vqMnhZjQanfk4LAyn2SLrd9p5hU43FWUmfTLJH0nskfXSTbNfr7pR0fOgQPfispINJ3i7pXdoE22h7u6RPSJpN8k5Jk1oYq1fK4GUl6TpJzyY5mWRO0kOSPjRwpouW5PkkR0bfv6KFB/32YVONh+0dkj4g6d6hs4yT7cslvVfSfZKUZC7JS4OGGp8pSTO2pyRdIulHA+dZt41QVtslPbfo8iltkv+pX2d7l6S9kp4YOMq43CPpU5LmB84xbrslvSDpC6ND3HtHw1JKS3Ja0t9I+qGk5yX9d5JHh021fhuhrDY129skfUXSXUleHjrPxbL9QUlnkhweOksPpiS9W9Lnk+yV9Jqk8s+h2n6LFo5Wdkv6NUmX2r5j2FTrtxHK6rSknYsu7xj9rDzbW7RQVA8k2SxjzG6QdKvt72vhkP1G218cNtLYnJJ0Ksnre8AHtFBe1f2+pP9M8kKSX2hhpN7vDJxp3TZCWX1X0tW2d49G1N8m6WsDZ7potq2F5z6OJ/nM0HnGJcmnk+xIsksLf6tvJin3r/RKkvxY0nO2rxn96CZJxwaMNC4/lPQe25eMHpc3qeCJg05zA/uU5Kztj0n6hhbOUtyf5OmBY43DDZI+LOnfbB8d/ewvkzwyXCR08HFJD4z+4Twp6SMD57loSZ6wfUDSES2cpX5SBd92w9ttAJSwEQ4DAWBNlBWAEigrACVQVgBKoKwAlEBZASiBsgJQwv8CNvjZcz/WaG4AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "with COGReader(\"masked.tif\") as cog:\n",
    "    imshow(cog.preview(max_size=10).as_masked()[0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "dental-hughes",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{\"bounds\": [-77.92959037999877, 25.148614663322203, -77.76954571373244, 25.254492610438696], \"center\": [-77.84956804686561, 25.201553636880448, 9], \"minzoom\": 9, \"maxzoom\": 9, \"band_metadata\": [[\"1\", {}], [\"2\", {}], [\"3\", {}]], \"band_descriptions\": [[\"1\", \"\"], [\"2\", \"\"], [\"3\", \"\"]], \"dtype\": \"uint8\", \"nodata_type\": \"Mask\", \"colorinterp\": [\"red\", \"green\", \"blue\"], \"scale\": null, \"offset\": null, \"colormap\": null, \"height\": 38, \"width\": 53, \"driver\": \"VRT\", \"count\": 3, \"overviews\": [2, 4, 8, 13]}\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAASsAAAD4CAYAAAC0ecCBAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAMTUlEQVR4nO3db4xddZ3H8c+nM20K2AAKa7QlS02wpjFKcUJAjMmW1RQxmBgf0EQTjcn4QLQYEoM+2+cbgg+MpiLuJrKQ3SoJGkRJhKCJotNSFVowWFFa0ZYYBarstNOPD+ZCZqZzO2fsPefMd3i/kglz77k5v++5c/n0d+7583USAcBKt6bvAgCgCcIKQAmEFYASCCsAJRBWAEoYb2Ol69aszzljG9pY9elOdXg0090NJUlarQdqu3wf3eFgHR9Zz8ypzsZ667ve0sk4e/fufT7JxYstayWszhnboKsv+HAbqz5Npk90Mo4kucsPvqTVelpJp+/j2lY+4os7cbK7sSTNvHS8s7EenLqnk3Fs/27YMnYDAZRAWAEogbACUAJhBaAEwgpACYQVgBIIKwAlEFYASiCsAJTQKKxs77D9lO2nbd/adlEAsNCSYWV7TNKXJV0naauknba3tl0YAMzVZGZ1paSnkxxKMi3pHkkfarcsAJivSVhtlPTsnMeHB8/NY3vS9pTtqelTL4+qPgCQNMIv2JPsTjKRZGLdmvWjWi0ASGoWVkckXTLn8abBcwDQmSZh9XNJl9nebHudpBsl3dduWQAw35J3Jkty0vZNkr4vaUzSnUmeaL0yAJij0W0Uk9wv6f6WawGAoTiDHUAJhBWAEggrACUQVgBKIKwAlEBYASiBsAJQQjvtapPOu9OuRnQuxhmdmum7gk4xswJQAmEFoATCCkAJhBWAEggrACUQVgBKIKwAlEBYASiBsAJQAmEFoIQmHZnvtH3U9uNdFAQAi2kys/ovSTtargMAzmjJsEryiKQ/d1ALAAw1skvtbU9KmpSk9T5vVKsFAEm0jwdQBEcDAZRAWAEoocmpC3dL+omkLbYP2/5k+2UBwHxLfsGeZGcXhQDAmbAbCKAEwgpACYQVgBIIKwAlEFYASiCsAJRAWAEooaX28VKSVla9UKct1jESXX02JMkdtqrvcrtei5hZASiBsAJQAmEFoATCCkAJhBWAEggrACUQVgBKIKwAlEBYASiBsAJQQpN7sF9i+yHbB2w/YXtXF4UBwFxNrg08KemWJPtsb5C01/aDSQ60XBsAvKpJ+/jnkuwb/P6ipIOSNrZdGADMtay7Lti+VNI2SY8usoz28QBa0/gLdtuvk/QtSTcneWHh8nnt4037eACj1SisbK/VbFDdleTb7ZYEAKdrcjTQkr4u6WCS29ovCQBO12RmdY2kj0nabnv/4OcDLdcFAPM0aR//Y0ncOxhArziDHUAJhBWAEggrACUQVgBKIKwAlEBYASiBsAJQAmEFoIRl3XVhJcrJk52N5fHyb9dwJzp8H93dOcZJOhury+16LWJmBaAEwgpACYQVgBIIKwAlEFYASiCsAJRAWAEogbACUAJhBaCEJg0j1tv+me1fDNrH/0cXhQHAXE2uH/l/SduTvDRoyfVj299L8tOWawOAVzVpGBFJLw0erh38dHfBFQCoeZPTMdv7JR2V9GCSRdvH256yPTWdl0dcJoDXukZhlWQmyeWSNkm60vbbF3kN7eMBtGZZRwOT/EXSQ5J2tFINAAzR5GjgxbYvGPx+jqT3SXqy5boAYJ4mRwPfJOm/bY9pNtz+N8l32y0LAOZrcjTwl5K2dVALAAzFGewASiCsAJRAWAEogbACUAJhBaAEwgpACYQVgBIIKwAl1O+HPjbWdwXtWdvhn6fD9vGrtaV7l9slSWs2bOhsrPev29nJOBt84buGLWNmBaAEwgpACYQVgBIIKwAlEFYASiCsAJRAWAEogbACUAJhBaAEwgpACY3DatDo9DHbNIsA0LnlzKx2STrYViEAcCZN28dvknS9pDvaLQcAFtd0ZnW7pM9LOjXsBbYnbU/ZnprOy6OoDQBe1aQj8wclHU2y90yvS7I7yUSSiXVeP7ICAUBqNrO6RtINtp+RdI+k7ba/2WpVALDAkmGV5AtJNiW5VNKNkn6Y5KOtVwYAc3CeFYASlnXf3CQPS3q4lUoA4AyYWQEogbACUAJhBaAEwgpACYQVgBIIKwAlEFYASmivP/nMTGurnsvj3bVY77o9uKZPdDZUl23WuxxLa7v7fHS4VZKkdPj5+Nv1V3QyzqmHfzR0GTMrACUQVgBKIKwAlEBYASiBsAJQAmEFoATCCkAJhBWAEggrACUQVgBKaHQtwqCzzYuSZiSdTDLRZlEAsNByLpz6tyTPt1YJAJwBu4EASmgaVpH0A9t7bU8u9gLaxwNoU9PdwPckOWL7XyQ9aPvJJI/MfUGS3ZJ2S9L5Yxd1fC8VAKtdo5lVkiOD/x6VdK+kK9ssCgAWWjKsbJ9ne8Mrv0t6v6TH2y4MAOZqshv4Rkn3Du7uOC7pf5I80GpVALDAkmGV5JCkd3ZQCwAMxakLAEogrACUQFgBKIGwAlACYQWgBMIKQAmEFYAS2uutPTbW2qrn6ryl+yrV5fu45qLXdzZWXnixs7Fu3/edzsaSpOseuamzsX69/audjHPVjmNDlzGzAlACYQWgBMIKQAmEFYASCCsAJRBWAEogrACUQFgBKIGwAlACYQWghEZhZfsC23tsP2n7oO2r2y4MAOZqem3glyQ9kOQjttdJOrfFmgDgNEuGle3zJb1X0sclKcm0pOl2ywKA+ZrsBm6WdEzSN2w/ZvuOQf/AeWgfD6BNTcJqXNIVkr6SZJuk45JuXfiiJLuTTCSZWOf1Iy4TwGtdk7A6LOlwkkcHj/doNrwAoDNLhlWSP0p61vaWwVPXSjrQalUAsEDTo4GfkXTX4EjgIUmfaK8kADhdo7BKsl/SRLulAMBwnMEOoATCCkAJhBWAEggrACUQVgBKIKwAlEBYASiBsAJQQtMz2JfHku1WVr1Qkk7GkSTNzHQ3liSPt/PnWczf371l6ReNyEN3fq2zse473t2t1z41uauzsSRpy9ShzsZ6x+RNnYzzzNHbhi5jZgWgBMIKQAmEFYASCCsAJRBWAEogrACUQFgBKIGwAlACYQWghCXDyvYW2/vn/Lxg++YOagOAVy15PUeSpyRdLkm2xyQdkXRvu2UBwHzL3Q28VtJvkvyujWIAYJjlXil7o6S7F1tge1LSpCStP727PACclcYzq0HPwBsk/d9iy+e1j19D+3gAo7Wc3cDrJO1L8qe2igGAYZYTVjs1ZBcQANrWKKxsnyfpfZK+3W45ALC4pu3jj0t6Q8u1AMBQnMEOoATCCkAJhBWAEggrACUQVgBKIKwAlEBYASiBsAJQgttov277mKTl3kbmIknPj7yYlWG1bhvbVc9K37Z/TXLxYgtaCat/hu2pJBN919GG1bptbFc9lbeN3UAAJRBWAEpYSWG1u+8CWrRat43tqqfstq2Y76wA4ExW0swKAIYirACUsCLCyvYO20/Zftr2rX3XMwq2L7H9kO0Dtp+wvavvmkbJ9pjtx2x/t+9aRsn2Bbb32H7S9kHbV/dd0yjY/tzgc/i47bttl+vq0ntYDRqnflmzDSm2Stppe2u/VY3ESUm3JNkq6SpJn14l2/WKXZIO9l1EC74k6YEkb5P0Tq2CbbS9UdJnJU0kebukMc221Sul97CSdKWkp5McSjIt6R5JH+q5prOW5Lkk+wa/v6jZD/3GfqsaDdubJF0v6Y6+axkl2+dLeq+kr0tSkukkf+m1qNEZl3SO7XFJ50r6Q8/1LNtKCKuNkp6d8/iwVsn/1K+wfamkbZIe7bmUUbld0uclneq5jlHbLOmYpG8MdnHvGDRLKS3JEUn/Ken3kp6T9NckP+i3quVbCWG1qtl+naRvSbo5yQt913O2bH9Q0tEke/uupQXjkq6Q9JUk2yQdl1T+O1TbF2p2b2WzpDdLOs/2R/utavlWQlgdkXTJnMebBs+VZ3utZoPqriSrpY3ZNZJusP2MZnfZt9v+Zr8ljcxhSYeTvDID3qPZ8Kru3yX9NsmxJCc021Lv3T3XtGwrIax+Luky25sHLepvlHRfzzWdNdvW7HcfB5Pc1nc9o5LkC0k2JblUs3+rHyYp96/0YpL8UdKztrcMnrpW0oEeSxqV30u6yva5g8/ltSp44KBR38A2JTlp+yZJ39fsUYo7kzzRc1mjcI2kj0n6le39g+e+mOT+/kpCA5+RdNfgH85Dkj7Rcz1nLcmjtvdI2qfZo9SPqeBlN1xuA6CElbAbCABLIqwAlEBYASiBsAJQAmEFoATCCkAJhBWAEv4BNUrh8GNdj68AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "with rasterio.Env(GDAL_DISABLE_READDIR_ON_OPEN=\"FALSE\"):\n",
    "    with COGReader(\"data.vrt\") as cog:\n",
    "        print(cog.info().json())\n",
    "        imshow(cog.preview(max_size=10).as_masked()[0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "possible-wireless",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQYAAAD8CAYAAACVSwr3AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAPFElEQVR4nO3dXYxc5X3H8e9vFzAtUMDxYrnGCQS5F1AlBK0oUmiUCgUIUmVyQ4mqxEqQnAuQkiqVShKpiXqVviSRIrWojoLiVGkoVRJhtbQNQZGiXpCwjojBUILDi7Bj8PIioE1i8Prfiz1rBj+73rV3ZmcWvh9pNGee85yZ/x6tf/uc58w5TlUhSb3Ghl2ApNFjMEhqGAySGgaDpIbBIKlhMEhqDCwYklyb5NEke5PcOqjPkdR/GcT3GJKMAz8HPgDsA+4HPlxVD/f9wyT13aBGDJcDe6vq8ap6FbgD2DKgz5LUZ6cM6H03Ak/3vN4H/MFCndetW1cXXHDBgEqRBLBr167nqmpiKX0HFQyLSrIN2Abw9re/nampqWGVIr0lJHlqqX0HdSixH9jU8/r8ru2oqtpeVZNVNTkxsaQQk7RCBhUM9wObk1yY5DTgRmDngD5LUp8N5FCiqg4nuQX4L2AcuL2q9gzisyT138DmGKrqbuDuQb2/pMHxm4+SGgaDpIbBIKlhMEhqGAySGgaDpIbBIKlhMEhqGAySGgaDpIbBIKlhMEhqGAySGgaDpIbBIKlhMEhqGAySGgaDpIbBIKlhMEhqGAySGgaDpIbBIKlhMEhqGAySGgaDpIbBIKlhMEhqGAySGgaDpIbBIKlhMEhqGAySGgaDpMYpy9k4yZPAK8AMcLiqJpOsBf4FuAB4Erihql5cXpmSVlI/Rgx/VFWXVtVk9/pW4N6q2gzc272WtIoM4lBiC7CjW94BXD+Az5A0QMsNhgK+n2RXkm1d2/qqOtAtPwOsn2/DJNuSTCWZmp6eXmYZkvppWXMMwJVVtT/JecA9Sf6nd2VVVZKab8Oq2g5sB5icnJy3j6ThWNaIoar2d88Hge8BlwPPJtkA0D0fXG6RklbWSQdDkjOSnDW3DFwNPATsBLZ23bYCdy23SEkrazmHEuuB7yWZe59/rqr/THI/cGeSm4CngBuWX6aklXTSwVBVjwPvnqf9eeCq5RQlabj85qOkhsEgqWEwSGoYDJIaBoOkhsEgqWEwSGoYDJIaBoOkhsEgqWEwSGoYDJIaBoOkhsEgqWEwSGoYDJIaBoOkhsEgqWEwSGoYDJIaBoOkhsEgqWEwSGoYDJIaBoOkhsEgqWEwSGoYDJIaBoOkhsEgqWEwSGoYDJIaBoOkxqLBkOT2JAeTPNTTtjbJPUke657P7dqT5KtJ9ibZneSyQRYvaTCWMmL4BnDtMW23AvdW1Wbg3u41wAeBzd1jG3Bbf8qUtJIWDYaq+hHwwjHNW4Ad3fIO4Pqe9m/WrPuAc5Js6FOtklbIyc4xrK+qA93yM8D6bnkj8HRPv31dm6RVZNmTj1VVQJ3odkm2JZlKMjU9Pb3cMiT10ckGw7Nzhwjd88GufT+wqaff+V1bo6q2V9VkVU1OTEycZBmSBuFkg2EnsLVb3grc1dP+0e7sxBXASz2HHJJWiVMW65Dk28D7gXVJ9gGfB74I3JnkJuAp4Iau+93AdcBe4FfAxwZQs6QBWzQYqurDC6y6ap6+Bdy83KIkDZfffJTUMBgkNQwGSY1F5xikUXLN2R+HI0dgbGz2GWaXAV57jSOHDnHPzJ3DK/BNwmDQUF1z9sdnF44cgSpIXl9Z3ffmett6ZHx8ttvMTNe9IA6C+8G9qOHq/upnfLwNgGNeZ3z89bBgNhBqZmZ2xNCFwz2H7xhsvW8RBoOG78iRo3/1c/qaN67rwiHj47N95sJi7jCiV/ceWj6DQcPVMwKgivr1b2b/gVe9vm5sjJqZOXrosOBbGQx94xyDRk9vAIyNkYSim0MY6/lbNjdqmG/0oGVxxKDhWmRegZmZ2UAAssAk5OvbjvGB8RuO30dL4ohBw7XQmYe5kUHP5GRVvfHsRdcnCVVFTj2Feu3wSlX+puaIQcM3dygw34ig9/sKx5McHVFcveZP+1jcW5PBoOHrnTfoHSkcEwhJZkcLc3MQXZ/qmaTMqQ6C+8Fg0OiaO1ToTlXW4cOzo4r5AmNuuQuNayc+wZUf+tuVq/VNxnjV8B05cnQUMHcG4mj72NjrI4KefsAbRhpHT2WOj1OXXMTL7/htXnrn8U9vamEGg4YrecM/9pr7DkPP5CLwxuV5XlcV2bSBp7as4/CZxfrJZ3j3WS8OsvI3NQ8lNFrGxtpJyGMOHcbOW8eRSy4kp58+exgxPs7Y2nN46ffXcmQNnPmu5zl43wb+cuO/r2Dhby4Gg0bHQmcgeuYaxs46E379G8afeIac8Vtk7TkcuWgjj95yPmfveZGxV+HFF87k9OfhrlfetcI/wJuHhxIart4w6L2MuvesxNgYY2f/Dq9eeB6MwfjLr7LvmrM5/fni1P8tzvzlIc5+DDg8wzv+9RlqzWk8d/ka/u0zV/EXd837qVqEwaDhm/va89yFVHMTkD2XYT/6Z5sYOxRu+5PtPPnqBP/4+B9y/abdnH/a8/zVrj/mgq8d4pdXn8fb9hzi9CeeY+zwuUy/21/vk5WqE/6/YvpucnKypqamhl2G9KaWZFdVTS6lr3MMkhoGg6SGwSCpYTBIahgMkhoGg6SGwSCpYTBIahgMkhoGg6SGwSCpYTBIahgMkhqLBkOS25McTPJQT9sXkuxP8kD3uK5n3WeS7E3yaJJrBlW4pMFZyojhG8C187R/paou7R53AyS5GLgRuKTb5h+SeEdOaZVZNBiq6kfAC0t8vy3AHVV1qKqeAPYCly+jPklDsJw5hluS7O4ONc7t2jYCT/f02de1NZJsSzKVZGp6enoZZUjqt5MNhtuAi4BLgQPAl070Dapqe1VNVtXkxMTESZYhaRBOKhiq6tmqmqmqI8DXeP1wYT+wqafr+V2bpFXkpIIhyYaelx8C5s5Y7ARuTLImyYXAZuAnyytR0kpb9Da6Sb4NvB9Yl2Qf8Hng/UkuBQp4EvgEQFXtSXIn8DBwGLi5qmYGUrmkgfEu0dJbhHeJlrQsBoOkhsEgqWEwSGoYDJIaBoOkhsEgqWEwSGoYDJIaBoOkhsEgqWEwSGoYDJIaBoOkhsEgqWEwSGoYDJIaBoOkhsEgqWEwSGoYDJIaBoOkhsEgqWEwSGoYDJIaBoOkhsEgqWEwSGoYDJIaBoOkhsEgqWEwSGoYDJIaiwZDkk1Jfpjk4SR7knyya1+b5J4kj3XP53btSfLVJHuT7E5y2aB/CEn9tZQRw2Hg01V1MXAFcHOSi4FbgXurajNwb/ca4IPA5u6xDbit71VLGqhFg6GqDlTVT7vlV4BHgI3AFmBH120HcH23vAX4Zs26DzgnyYZ+Fy5pcE5ojiHJBcB7gB8D66vqQLfqGWB9t7wReLpns31dm6RVYsnBkORM4DvAp6rq5d51VVVAncgHJ9mWZCrJ1PT09IlsKmnAlhQMSU5lNhS+VVXf7ZqfnTtE6J4Pdu37gU09m5/ftb1BVW2vqsmqmpyYmDjZ+iUNwFLOSgT4OvBIVX25Z9VOYGu3vBW4q6f9o93ZiSuAl3oOOSStAqcsoc97gY8ADyZ5oGv7LPBF4M4kNwFPATd06+4GrgP2Ar8CPtbPgiUN3qLBUFX/DWSB1VfN07+Am5dZl6Qh8puPkhoGg6SGwSCpYTBIahgMkhoGg6SGwSCpYTBIahgMkhoGg6SGwSCpYTBIahgMkhoGg6SGwSCpYTBIahgMkhoGg6SGwSCpYTBIahgMkhoGg6SGwSCpYTBIahgMkhoGg6SGwSCpYTBIahgMkhoGg6SGwSCpYTBIahgMkhoGg6TGosGQZFOSHyZ5OMmeJJ/s2r+QZH+SB7rHdT3bfCbJ3iSPJrlmkD+ApP47ZQl9DgOfrqqfJjkL2JXknm7dV6rq73o7J7kYuBG4BPhd4AdJfq+qZvpZuKTBWXTEUFUHquqn3fIrwCPAxuNssgW4o6oOVdUTwF7g8n4UK2llnNAcQ5ILgPcAP+6abkmyO8ntSc7t2jYCT/dsto95giTJtiRTSaamp6dPvHJJA7PkYEhyJvAd4FNV9TJwG3ARcClwAPjSiXxwVW2vqsmqmpyYmDiRTSUN2JKCIcmpzIbCt6rquwBV9WxVzVTVEeBrvH64sB/Y1LP5+V2bpFViKWclAnwdeKSqvtzTvqGn24eAh7rlncCNSdYkuRDYDPykfyVLGrSlnJV4L/AR4MEkD3RtnwU+nORSoIAngU8AVNWeJHcCDzN7RuNmz0hIq0uqatg1kGQa+D/guWHXsgTrWB11wuqp1Tr7b75a31FVS5rQG4lgAEgyVVWTw65jMaulTlg9tVpn/y23Vr8SLalhMEhqjFIwbB92AUu0WuqE1VOrdfbfsmodmTkGSaNjlEYMkkbE0IMhybXd5dl7k9w67HqOleTJJA92l5ZPdW1rk9yT5LHu+dzF3mcAdd2e5GCSh3ra5q0rs77a7ePdSS4bgVpH7rL949xiYKT264rcCqGqhvYAxoFfAO8ETgN+Blw8zJrmqfFJYN0xbX8D3Not3wr89RDqeh9wGfDQYnUB1wH/AQS4AvjxCNT6BeDP5+l7cfd7sAa4sPv9GF+hOjcAl3XLZwE/7+oZqf16nDr7tk+HPWK4HNhbVY9X1avAHcxetj3qtgA7uuUdwPUrXUBV/Qh44ZjmheraAnyzZt0HnHPMV9oHaoFaFzK0y/Zr4VsMjNR+PU6dCznhfTrsYFjSJdpDVsD3k+xKsq1rW19VB7rlZ4D1wymtsVBdo7qfT/qy/UE75hYDI7tf+3krhF7DDobV4Mqqugz4IHBzkvf1rqzZsdrIndoZ1bp6LOuy/UGa5xYDR43Sfu33rRB6DTsYRv4S7ara3z0fBL7H7BDs2bkhY/d8cHgVvsFCdY3cfq4RvWx/vlsMMIL7ddC3Qhh2MNwPbE5yYZLTmL1X5M4h13RUkjO6+1yS5AzgamYvL98JbO26bQXuGk6FjYXq2gl8tJtFvwJ4qWdoPBSjeNn+QrcYYMT260J19nWfrsQs6iIzrNcxO6v6C+Bzw67nmNreyexs7s+APXP1AW8D7gUeA34ArB1Cbd9mdrj4GrPHjDctVBezs+Z/3+3jB4HJEaj1n7padne/uBt6+n+uq/VR4IMrWOeVzB4m7AYe6B7Xjdp+PU6dfdunfvNRUmPYhxKSRpDBIKlhMEhqGAySGgaDpIbBIKlhMEhqGAySGv8PGMdkXqTjyOMAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "with rasterio.Env(GDAL_DISABLE_READDIR_ON_OPEN=\"FALSE\"):\n",
    "    with COGReader(\"data.vrt\") as cog:\n",
    "        imshow(cog.tile(72, 109, 8, tilesize=256).as_masked()[0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "iraqi-shelter",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "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.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}