warpedvrt_mask.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import rasterio\n",
    "from rasterio.plot import show\n",
    "from rasterio.vrt import WarpedVRT"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "([<MaskFlags.per_dataset: 2>], [<MaskFlags.per_dataset: 2>], [<MaskFlags.per_dataset: 2>])\n",
      "0.0\n"
     ]
    }
   ],
   "source": [
    "with rasterio.open(\"masked.tif\") as src:\n",
    "    print(src.mask_flag_enums)\n",
    "    print(src.nodata)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVEAAAD5CAYAAACeTgZ1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAPBUlEQVR4nO3db4ylZXnH8e+vdFgUNLiKmy1LClVaQowMyXSL0Re4FrulpmBijLQxvCBdm0iiiWmLNqmatAkmFcqLxmQVyppQkYIGQmh1u5JsTBpkwGVZ/hSQYmS7shogQJsuu3D1xXlIRpjZOTP3OXP+zPeTnMx57nNmn+sOZ37c57nOc55UFZKk1fm1URcgSZPMEJWkBoaoJDUwRCWpgSEqSQ0MUUlq8Ostv5xkO3AdcALwjaq6+njPPzEb6iRObtmlJtBvv/d/R13CcT22/82jLkFj7v/4H16uI1nssaz2c6JJTgAeAy4CngbuBS6rqoeX+p23ZmP9Xj60qv1pcn3vv/eNuoTj+oPfmB11CRpz99QeXqhnFw3RlrfzW4EnqurJqnoZuBm4pOHfk6SJ0xKipwM/W7D9dDf2K5LsSDKfZP4oRxp2J0njZ+iNparaWVVzVTU3w4Zh706S1lRLiB4EzliwvaUbk6R1o6U7fy9wdpKz6IXnJ4A/GUhVE2zcmyh6o6X+m9lwUj9WHaJVdSzJlcD36H3E6YaqemhglUnSBGj6nGhV3QXcNaBaJGnieMaSJDUwRCWpQdPb+fXOJtJ0s+GkfrgSlaQGhqgkNTBEJamBISpJDQxRSWpgd15aIbv2WsiVqCQ1MEQlqYEhKkkNDFFJamCISlIDu/N98Bx59WMlrxM7+dPDlagkNTBEJamBISpJDQxRSWpgiEpSg6bufJKngBeBV4BjVTU3iKKkaef599NjEB9x+mBV/XIA/44kTRzfzktSg9YQLeD7Se5LsmOxJyTZkWQ+yfxRjjTuTpLGS+vb+Q9U1cEk7wR2J3m0qvYufEJV7QR2Arw1G6txf5I0VppCtKoOdj8PJ/kusBXYe/zfGl+e3qlRW+w1aLNpvK367XySk5O85bX7wIeBA4MqTJImQctKdBPw3SSv/Tv/XFX/NpCqJGlCrDpEq+pJ4LwB1iJJE8ePOElSA0NUkhr4pczSmPMU0fHmSlSSGhiiktTAEJWkBoaoJDUwRCWpwbrsznuOvKaBXfvx4EpUkhoYopLUwBCVpAaGqCQ1MEQlqcG67M5L08yu/dpyJSpJDQxRSWpgiEpSA0NUkhos21hKcgPwEeBwVb2nG9sIfBs4E3gK+HhVPTe8MlfPUzylHi/HPBz9rERvBLa/buwqYE9VnQ3s6bYlad1ZNkSrai/w7OuGLwF2dfd3AZcOtixJmgyr/Zzopqo61N3/Ob1r0C8qyQ5gB8BJvHmVu5Ok8dTcWKqqAuo4j++sqrmqmpthQ+vuJGmsrDZEn0myGaD7eXhwJUnS5Fjt2/k7gMuBq7uftw+sIklrxlNE2y27Ek3yLeA/gN9J8nSSK+iF50VJHgd+v9uWpHVn2ZVoVV22xEMfGnAtkjRxPGNJkhoYopLUwBCVpAZT86XMniMvrZxd+HauRCWpgSEqSQ0MUUlqYIhKUgNDVJIaGKKS1MAQlaQGhqgkNTBEJamBISpJDSbutE9P75QGZ6V/T54m+kauRCWpgSEqSQ0MUUlqYIhKUgNDVJIa9HO1zxuSHE5yYMHYl5IcTLKvu1083DIlaTz1sxK9Edi+yPi1VTXb3e4abFmSNBmWDdGq2gs8uwa1SNLEaTkmemWS/d3b/bct9aQkO5LMJ5k/ypGG3UnS+FltiH4NeBcwCxwCvrrUE6tqZ1XNVdXcDBtWuTtJGk+rCtGqeqaqXqmqV4GvA1sHW5YkTYZVnTufZHNVHeo2PwocON7zV8Nz5KXxs5K/y/Vynv2yIZrkW8CFwDuSPA18EbgwySxQwFPAp4ZXoiSNr2VDtKouW2T4+iHUIkkTxzOWJKmBISpJDQxRSWowcd9sL2kyDOoTNuPe5XclKkkNDFFJamCISlIDQ1SSGoxFY8lTPCUtZbF8GKdmkytRSWpgiEpSA0NUkhoYopLUwBCVpAZj0Z2XpJVY6hM9o+jauxKVpAaGqCQ1MEQlqYEhKkkNDFFJatDP1T7PAL4JbKJ3dc+dVXVdko3At4Ez6V3x8+NV9dxqiliso+b59JJWahRd+35WoseAz1XVucAFwKeTnAtcBeypqrOBPd22JK0ry4ZoVR2qqvu7+y8CjwCnA5cAu7qn7QIuHVKNkjS2VvRh+yRnAucD9wCbqupQ99DP6b3dX+x3dgA7AE7izasuVJLGUd+NpSSnALcBn62qFxY+VlVF73jpG1TVzqqaq6q5GTY0FStJ46avlWiSGXoBelNVfacbfibJ5qo6lGQzcHiQhS11INiGk6SVGmbDadmVaJIA1wOPVNU1Cx66A7i8u385cHtzNZI0YfpZib4f+CTwYJJ93dgXgKuBW5JcAfwU+PhQKpSkMbZsiFbVD4Es8fCHBluOJE0Wz1iSpAaGqCQ1mLgvZbZrL2lQBnE5ZleiktTAEJWkBoaoJDUwRCWpgSEqSQ0mrju/FLv2kkbBlagkNTBEJamBISpJDQxRSWpgiEpSg6npzi/FyzFLGiZXopLUwBCVpAaGqCQ1MEQlqcGyjaUkZwDfBDbRu7b8zqq6LsmXgD8DftE99QtVddewCh0kTxGVNCj9dOePAZ+rqvuTvAW4L8nu7rFrq+rvh1eeJI23fq72eQg41N1/MckjwOnDLkySJsGKjokmORM4H7inG7oyyf4kNyR52xK/syPJfJL5oxxpq1aSxkzfIZrkFOA24LNV9QLwNeBdwCy9lepXF/u9qtpZVXNVNTfDhvaKJWmM9BWiSWboBehNVfUdgKp6pqpeqapXga8DW4dXpiSNp3668wGuBx6pqmsWjG/ujpcCfBQ4MJwS145de0kr1U93/v3AJ4EHk+zrxr4AXJZklt7Hnp4CPjWE+iRprPXTnf8hkEUemojPhErSMHnGkiQ1MEQlqYEhKkkNpv5LmQfBrr2kpbgSlaQGhqgkNTBEJamBISpJDQxRSWpgd76Bl2OWJttSn7xZCVeiktTAEJWkBoaoJDUwRCWpgY2lAfMUUWl9cSUqSQ0MUUlqYIhKUgNDVJIaGKKS1KCfSyafBOwFNnTPv7WqvpjkLOBm4O3AfcAnq+rlYRY7yezaS9Opn5XoEWBbVZ0HzALbk1wAfAW4tqreDTwHXDG0KiVpTC0botXzUrc5090K2Abc2o3vAi4dRoGSNM76Oiaa5IQk+4DDwG7gJ8DzVXWse8rTwOlL/O6OJPNJ5o9yZAAlS9L46CtEq+qVqpoFtgBbgXP63UFV7ayquaqam2HD6qqUpDG1ou58VT0P3A28Dzg1yWuNqS3AwcGWJknjr5/u/GnA0ap6PsmbgIvoNZXuBj5Gr0N/OXD7MAudVnbtpcnWzxeQbAZ2JTmB3sr1lqq6M8nDwM1J/hb4MXD9EOuUpLG0bIhW1X7g/EXGn6R3fFSS1i3PWJKkBoaoJDUwRCWpgd9sP6a8HLM0GVyJSlIDQ1SSGhiiktTAEJWkBjaWJoiniEqrs9TfziC4EpWkBoaoJDUwRCWpgSEqSQ0MUUlqYHd+Cti1l0bHlagkNTBEJamBISpJDQxRSWqwbIgmOSnJj5I8kOShJF/uxm9M8l9J9nW32aFXK0ljpp/u/BFgW1W9lGQG+GGSf+0e+4uqunV45amFXXuN0jDPVx8n/Vzts4CXus2Z7lbDLEqSJkVfx0STnJBkH3AY2F1V93QP/V2S/UmuTbJhid/dkWQ+yfxRjgymakkaE32FaFW9UlWzwBZga5L3AJ8HzgF+F9gI/NUSv7uzquaqam6GRXNWkibWirrzVfU8cDewvaoOVc8R4J+ArUOoT5LGWj/d+dOSnNrdfxNwEfBoks3dWIBLgQPDK1OSxlM/3fnNwK4kJ9AL3Vuq6s4kP0hyGhBgH/DnwytTg7SSrum4d/LXSwdY46uf7vx+4PxFxrcNpSJJmiCesSRJDQxRSWpgiEpSA7+UWcdl40Y6PleiktTAEJWkBoaoJDUwRCWpgSEqSQ3S+7rQNdpZ8gvgp93mO4BfrtnOR8M5Tof1MEdYH/Nc7Rx/s6pOW+yBNQ3RX9lxMl9VcyPZ+RpxjtNhPcwR1sc8hzFH385LUgNDVJIajDJEd45w32vFOU6H9TBHWB/zHPgcR3ZMVJKmgW/nJamBISpJDdY8RJNsT/KfSZ5IctVa739YktyQ5HCSAwvGNibZneTx7ufbRlljqyRnJLk7ycNJHkrymW58auaZ5KQkP0ryQDfHL3fjZyW5p3vdfjvJiaOutVV3KfQfJ7mz256qOSZ5KsmDSfYlme/GBv5aXdMQ7a7T9I/AHwLnApclOXctaxiiG4Htrxu7CthTVWcDe7rtSXYM+FxVnQtcAHy6++83TfM8AmyrqvOAWWB7kguArwDXVtW7geeAK0ZX4sB8BnhkwfY0zvGDVTW74LOhA3+trvVKdCvwRFU9WVUvAzcDl6xxDUNRVXuBZ183fAmwq7u/i95VUSdWd5ns+7v7L9L7AzydKZpndxnwl7rNme5WwDbg1m58oucIkGQL8EfAN7rtMGVzXMLAX6trHaKnAz9bsP10NzatNlXVoe7+z4FNoyxmkJKcSe8ChvcwZfPs3ubuAw4Du4GfAM9X1bHuKdPwuv0H4C+BV7vttzN9cyzg+0nuS7KjGxv4a9Vvtl8jVVVJpuLzZElOAW4DPltVL/QWMT3TMM+qegWYTXIq8F3gnNFWNFhJPgIcrqr7klw44nKG6QNVdTDJO4HdSR5d+OCgXqtrvRI9CJyxYHtLNzatnkmyGaD7eXjE9TRLMkMvQG+qqu90w1M3T4Cqeh64G3gfcGqS1xYdk/66fT/wx0meondIbRtwHdM1R6rqYPfzML3/GW5lCK/VtQ7Re4Gzuy7gicAngDvWuIa1dAdweXf/cuD2EdbSrDtudj3wSFVds+ChqZlnktO6FShJ3gRcRO/Y793Ax7qnTfQcq+rzVbWlqs6k9zf4g6r6U6ZojklOTvKW1+4DHwYOMIzXalWt6Q24GHiM3nGmv17r/Q9xXt8CDgFH6R1PuoLecaY9wOPAvwMbR11n4xw/QO84035gX3e7eJrmCbwX+HE3xwPA33TjvwX8CHgC+Bdgw6hrHdB8LwTunLY5dnN5oLs99FrWDOO16mmfktTAM5YkqYEhKkkNDFFJamCISlIDQ1SSGhiiktTAEJWkBv8PoZqtwMRcftwAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "with rasterio.open(\"masked.tif\") as src:\n",
    "    show(src.dataset_mask())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVEAAAD5CAYAAACeTgZ1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAObElEQVR4nO3db4hl9X3H8fenk020UVA3G9m6SzdNpEFCHWG6GMwDs6nJ1oZqIITYEnwgbAoJKEhbk0KbQAsGGm0ehMAmWveBVdP8QRFbs90sSKCsGXVdd10bjTXE7epuVkVtqXbHbx/cI0xkxrk7v3tn7r3zfsHlnvM75875/vDsx3PP755zUlVIkpbnN1a7AEkaZ4aoJDUwRCWpgSEqSQ0MUUlqYIhKUoN3tHw4yXbgG8AU8J2quvHt1n/POVO1ZfO6lk1KGnM/O/Cbq13CKftf/pvX67UstGzZIZpkCvgmcBnwLPDTJPdU1eOLfWbL5nU8eP/m5W5S0gT4xG9Nr3YJp2xf7Vl0WcvX+a3AU1X1dFW9DtwJXNHw9yRp7LSE6HnAL+fNP9u1/ZokO5LMJpk9fmKuYXOSNHqGPrBUVTuraqaqZjasnxr25iRpRbWE6BFg/gnOTV2bJK0ZLaPzPwXOT/I+euH5WeBPBlKVpIl1/3/tX7B9HAecoCFEq+pkki8C99P7idOtVXVoYJVJ0hho+p1oVd0H3DegWiRp7HjFkiQ1MEQlqUHT13lJOlWnOoC00EDUKA1CeSQqSQ0MUUlqYIhKUgNDVJIaGKKS1MDReUkjbZRG4hfikagkNTBEJamBISpJDQxRSWpgiEpSA0NUkhoYopLUwBCVpAaGqCQ1MEQlqYEhKkkNmq6dT/IM8AowB5ysqplBFCVp7RnXRykP4gYkH62qXw3g70jS2PHrvCQ1aA3RAn6U5KEkOxZaIcmOJLNJZo+fmGvcnCSNltav8x+pqiNJ3gvsTvJEVT0wf4Wq2gnsBJi58LRq3J4kjZSmEK2qI937sSQ/BLYCD7z9pyRNmkEM/oz6ANJilv11Psm7k5z55jTwceDgoAqTpHHQciR6LvDDJG/+nX+qqn8dSFWSNCaWHaJV9TRw4QBrkaSx40+cJKmBISpJDQxRSWpgiEpSA0NUkhoYopLUwBCVpAaGqCQ1GMT9RCWtcYvdUHkhi10jP643ZfZIVJIaGKKS1MAQlaQGhqgkNTBEJamBo/OShuJUR9VHfRR+MR6JSlIDQ1SSGhiiktTAEJWkBksOLCW5FfgkcKyqPtS1nQPcBWwBngE+U1UvDq9MSaNgXAd/hqmfI9HbgO1vabsB2FNV5wN7unlJWnOWDNGqegB44S3NVwC7uuldwJWDLUuSxsNyz4meW1VHu+nn6D2DfkFJdiSZTTJ7/MTcMjcnSaOpeWCpqgqot1m+s6pmqmpmw/qp1s1J0khZbog+n2QjQPd+bHAlSdL4WO5ln/cAVwM3du93D6wiSSNrXG+cPExLHokmuQP4d+B3kzyb5Bp64XlZkieBP+jmJWnNWfJItKquWmTRxwZciySNHa9YkqQGhqgkNTBEJamBN2WW1Le1PAq/GI9EJamBISpJDQxRSWpgiEpSA0NUkhoYopLUwBCVpAaGqCQ1MEQlqYEhKkkNvOxTUt8WuynzqZi0S0c9EpWkBoaoJDUwRCWpgSEqSQ0MUUlq0M/TPm9NcizJwXltX0lyJMn+7nX5cMuUpNHUz5HobcD2Bdpvrqrp7nXfYMuSpPGwZIhW1QPACytQiySNnZZzol9McqD7un/2Yisl2ZFkNsns8RNzDZuTpNGz3BD9FvB+YBo4Cnx9sRWramdVzVTVzIb1U8vcnCSNpmWFaFU9X1VzVfUG8G1g62DLkqTxsKxr55NsrKqj3eyngINvt74kvWmx6+/H9Zr6JUM0yR3ApcB7kjwL/A1waZJpoIBngM8Pr0RJGl1LhmhVXbVA8y1DqEWSxo5XLElSA0NUkhoYopLUwDvbSxoJg7hr/mKGOfLvkagkNTBEJamBISpJDQxRSWrgwJKkibEal456JCpJDQxRSWpgiEpSA0NUkhoYopLUwNF5SRNjWJeObv3E/yy6zCNRSWpgiEpSA0NUkhoYopLUwBCVpAZLhmiSzUn2Jnk8yaEk13bt5yTZneTJ7v3s4ZcrSaOlnyPRk8D1VXUBcDHwhSQXADcAe6rqfGBPNy9Ja8qSIVpVR6vq4W76FeAwcB5wBbCrW20XcOWQapSkkXVK50STbAEuAvYB51bV0W7Rc8C5i3xmR5LZJLPHT8y11CpJI6fvEE1yBvB94Lqqenn+sqoqoBb6XFXtrKqZqprZsH6qqVhJGjV9hWiSdfQC9Paq+kHX/HySjd3yjcCx4ZQoSaOrn9H5ALcAh6vqpnmL7gGu7qavBu4efHmSNNr6uQHJJcDngMeS7O/avgzcCHw3yTXAL4DPDKVCSRphS4ZoVf0EyCKLPzbYciRpvHjFkiQ1MEQlqYEhKkkNDFFJamCISlIDQ1SSGhiiktTAEJWkBoaoJDUwRCWpgSEqSQ0MUUlqYIhKUgNDVJIaGKKS1MAQlaQGhqgkNTBEJamBISpJDfp52ufmJHuTPJ7kUJJru/avJDmSZH/3unz45UrSaOnnaZ8ngeur6uEkZwIPJdndLbu5qv5+eOVJ0mjr52mfR4Gj3fQrSQ4D5w27MEkaB6d0TjTJFuAiYF/X9MUkB5LcmuTsRT6zI8lsktnjJ+baqpWkEdN3iCY5A/g+cF1VvQx8C3g/ME3vSPXrC32uqnZW1UxVzWxYP9VesSSNkL5CNMk6egF6e1X9AKCqnq+quap6A/g2sHV4ZUrSaOpndD7ALcDhqrppXvvGeat9Cjg4+PIkabT1Mzp/CfA54LEk+7u2LwNXJZkGCngG+PwQ6pOkkdbP6PxPgCyw6L7BlyNJ48UrliSpgSEqSQ0MUUlqYIhKUgNDVJIaGKKS1MAQlaQGhqgkNTBEJamBISpJDQxRSWpgiEpSA0NUkhoYopLUwBCVpAaGqCQ1MEQlqYEhKkkNDFFJamCISlKDfh6ZfFqSB5M8muRQkq927e9Lsi/JU0nuSvLO4ZcrSaOlnyPR14BtVXUhMA1sT3Ix8DXg5qr6APAicM3QqpSkEbVkiFbPq93suu5VwDbge137LuDKYRQoSaOsr3OiSaaS7AeOAbuBnwMvVdXJbpVngfMW+eyOJLNJZo+fmBtAyZI0OvoK0aqaq6ppYBOwFfhgvxuoqp1VNVNVMxvWTy2vSkkaUac0Ol9VLwF7gQ8DZyV5R7doE3BksKVJ0ujrZ3R+Q5KzuunTgcuAw/TC9NPdalcDdw+pRkkaWe9YehU2AruSTNEL3e9W1b1JHgfuTPK3wCPALUOsU5JG0pIhWlUHgIsWaH+a3vlRSVqzvGJJkhoYopLUwBCVpAaGqCQ1MEQlqYEhKkkNDFFJamCISlIDQ1SSGhiiktTAEJWkBoaoJDUwRCWpgSEqSQ0MUUlqYIhKUgNDVJIaGKKS1MAQlaQG/Tzt87QkDyZ5NMmhJF/t2m9L8p9J9nev6aFXK0kjpp+nfb4GbKuqV5OsA36S5F+6ZX9eVd8bXnmSNNr6edpnAa92s+u6Vw2zKEkaF32dE00ylWQ/cAzYXVX7ukV/l+RAkpuTvGuRz+5IMptk9viJucFULUkjoq8Qraq5qpoGNgFbk3wI+BLwQeD3gXOAv1zkszuraqaqZjasnxpM1ZI0Ik5pdL6qXgL2Atur6mj1vAb8I7B1CPVJ0kjrZ3R+Q5KzuunTgcuAJ5Js7NoCXAkcHF6ZkjSa+hmd3wjsSjJFL3S/W1X3Jvlxkg1AgP3Anw2vTEkaTf2Mzh8ALlqgfdtQKpKkMeIVS5LUwBCVpAaGqCQ1MEQlqYEhKkkNDFFJamCISlIDQ1SSGhiiktQgvduFrtDGkuPAL7rZ9wC/WrGNrw77OBnWQh9hbfRzuX387arasNCCFQ3RX9twMltVM6uy8RViHyfDWugjrI1+DqOPfp2XpAaGqCQ1WM0Q3bmK214p9nEyrIU+wtro58D7uGrnRCVpEvh1XpIaGKKS1GDFQzTJ9iT/keSpJDes9PaHJcmtSY4lOTiv7Zwku5M82b2fvZo1tkqyOcneJI8nOZTk2q59YvqZ5LQkDyZ5tOvjV7v29yXZ1+23dyV552rX2qp7FPojSe7t5ieqj0meSfJYkv1JZru2ge+rKxqi3XOavgn8IXABcFWSC1ayhiG6Ddj+lrYbgD1VdT6wp5sfZyeB66vqAuBi4Avdf79J6udrwLaquhCYBrYnuRj4GnBzVX0AeBG4ZvVKHJhrgcPz5iexjx+tqul5vw0d+L660keiW4GnqurpqnoduBO4YoVrGIqqegB44S3NVwC7uuld9J6KOra6x2Q/3E2/Qu8f4HlMUD+7x4C/2s2u614FbAO+17WPdR8BkmwC/gj4TjcfJqyPixj4vrrSIXoe8Mt58892bZPq3Ko62k0/B5y7msUMUpIt9B5guI8J62f3NXc/cAzYDfwceKmqTnarTMJ++w/AXwBvdPPrmbw+FvCjJA8l2dG1DXxf7eeRyRqAqqokE/F7siRnAN8Hrquql3sHMT2T0M+qmgOmk5wF/BD44OpWNFhJPgkcq6qHkly6yuUM00eq6kiS9wK7kzwxf+Gg9tWVPhI9AmyeN7+pa5tUzyfZCNC9H1vlepolWUcvQG+vqh90zRPXT4CqegnYC3wYOCvJmwcd477fXgL8cZJn6J1S2wZ8g8nqI1V1pHs/Ru9/hlsZwr660iH6U+D8bhTwncBngXtWuIaVdA9wdTd9NXD3KtbSrDtvdgtwuKpumrdoYvqZZEN3BEqS04HL6J373Qt8ulttrPtYVV+qqk1VtYXev8EfV9WfMkF9TPLuJGe+OQ18HDjIMPbVqlrRF3A58DN655n+aqW3P8R+3QEcBf6P3vmka+idZ9oDPAn8G3DOatfZ2MeP0DvPdADY370un6R+Ar8HPNL18SDw11377wAPAk8B/wy8a7VrHVB/LwXunbQ+dn15tHsdejNrhrGvetmnJDXwiiVJamCISlIDQ1SSGhiiktTAEJWkBoaoJDUwRCWpwf8DCuUwyzBK934AAAAASUVORK5CYII=\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 rasterio.open(\"masked.tif\") as src:\n",
    "        with WarpedVRT(src, add_alpha=True) as vrt:\n",
    "            show(vrt.dataset_mask())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Creating output file that is 53P x 38L.\r\n",
      "Processing input file masked.tif.\r\n",
      "Using internal nodata values (e.g. 0) for image masked.tif.\r\n",
      "Copying nodata values from source masked.tif to destination warped.tif.\r\n",
      "0...10...20...30...40...50...60...70...80...90...100 - done.\r\n"
     ]
    }
   ],
   "source": [
    "!gdalwarp -dstalpha masked.tif warped.tif"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/denis/satellogic/env/lib/python3.6/site-packages/ipykernel_launcher.py:2: NodataShadowWarning: The dataset's nodata attribute is shadowing the alpha band. All masks will be determined by the nodata attribute\n",
      "  \n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVEAAAD5CAYAAACeTgZ1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAPBUlEQVR4nO3db4ylZXnH8e+vdFgUNLiKmy1LClVaQowMyXSL0Re4FrulpmBijLQxvCBdm0iiiWmLNqmatAkmFcqLxmQVyppQkYIGQmh1u5JsTBpkwGVZ/hSQYmS7shogQJsuu3D1xXlIRpjZOTP3OXP+zPeTnMx57nNmn+sOZ37c57nOc55UFZKk1fm1URcgSZPMEJWkBoaoJDUwRCWpgSEqSQ0MUUlq8Ostv5xkO3AdcALwjaq6+njPPzEb6iRObtmlJtBvv/d/R13CcT22/82jLkFj7v/4H16uI1nssaz2c6JJTgAeAy4CngbuBS6rqoeX+p23ZmP9Xj60qv1pcn3vv/eNuoTj+oPfmB11CRpz99QeXqhnFw3RlrfzW4EnqurJqnoZuBm4pOHfk6SJ0xKipwM/W7D9dDf2K5LsSDKfZP4oRxp2J0njZ+iNparaWVVzVTU3w4Zh706S1lRLiB4EzliwvaUbk6R1o6U7fy9wdpKz6IXnJ4A/GUhVE2zcmyh6o6X+m9lwUj9WHaJVdSzJlcD36H3E6YaqemhglUnSBGj6nGhV3QXcNaBaJGnieMaSJDUwRCWpQdPb+fXOJtJ0s+GkfrgSlaQGhqgkNTBEJamBISpJDQxRSWpgd15aIbv2WsiVqCQ1MEQlqYEhKkkNDFFJamCISlIDu/N98Bx59WMlrxM7+dPDlagkNTBEJamBISpJDQxRSWpgiEpSg6bufJKngBeBV4BjVTU3iKKkaef599NjEB9x+mBV/XIA/44kTRzfzktSg9YQLeD7Se5LsmOxJyTZkWQ+yfxRjjTuTpLGS+vb+Q9U1cEk7wR2J3m0qvYufEJV7QR2Arw1G6txf5I0VppCtKoOdj8PJ/kusBXYe/zfGl+e3qlRW+w1aLNpvK367XySk5O85bX7wIeBA4MqTJImQctKdBPw3SSv/Tv/XFX/NpCqJGlCrDpEq+pJ4LwB1iJJE8ePOElSA0NUkhr4pczSmPMU0fHmSlSSGhiiktTAEJWkBoaoJDUwRCWpwbrsznuOvKaBXfvx4EpUkhoYopLUwBCVpAaGqCQ1MEQlqcG67M5L08yu/dpyJSpJDQxRSWpgiEpSA0NUkhos21hKcgPwEeBwVb2nG9sIfBs4E3gK+HhVPTe8MlfPUzylHi/HPBz9rERvBLa/buwqYE9VnQ3s6bYlad1ZNkSrai/w7OuGLwF2dfd3AZcOtixJmgyr/Zzopqo61N3/Ob1r0C8qyQ5gB8BJvHmVu5Ok8dTcWKqqAuo4j++sqrmqmpthQ+vuJGmsrDZEn0myGaD7eXhwJUnS5Fjt2/k7gMuBq7uftw+sIklrxlNE2y27Ek3yLeA/gN9J8nSSK+iF50VJHgd+v9uWpHVn2ZVoVV22xEMfGnAtkjRxPGNJkhoYopLUwBCVpAZT86XMniMvrZxd+HauRCWpgSEqSQ0MUUlqYIhKUgNDVJIaGKKS1MAQlaQGhqgkNTBEJamBISpJDSbutE9P75QGZ6V/T54m+kauRCWpgSEqSQ0MUUlqYIhKUgNDVJIa9HO1zxuSHE5yYMHYl5IcTLKvu1083DIlaTz1sxK9Edi+yPi1VTXb3e4abFmSNBmWDdGq2gs8uwa1SNLEaTkmemWS/d3b/bct9aQkO5LMJ5k/ypGG3UnS+FltiH4NeBcwCxwCvrrUE6tqZ1XNVdXcDBtWuTtJGk+rCtGqeqaqXqmqV4GvA1sHW5YkTYZVnTufZHNVHeo2PwocON7zV8Nz5KXxs5K/y/Vynv2yIZrkW8CFwDuSPA18EbgwySxQwFPAp4ZXoiSNr2VDtKouW2T4+iHUIkkTxzOWJKmBISpJDQxRSWowcd9sL2kyDOoTNuPe5XclKkkNDFFJamCISlIDQ1SSGoxFY8lTPCUtZbF8GKdmkytRSWpgiEpSA0NUkhoYopLUwBCVpAZj0Z2XpJVY6hM9o+jauxKVpAaGqCQ1MEQlqYEhKkkNDFFJatDP1T7PAL4JbKJ3dc+dVXVdko3At4Ez6V3x8+NV9dxqiliso+b59JJWahRd+35WoseAz1XVucAFwKeTnAtcBeypqrOBPd22JK0ry4ZoVR2qqvu7+y8CjwCnA5cAu7qn7QIuHVKNkjS2VvRh+yRnAucD9wCbqupQ99DP6b3dX+x3dgA7AE7izasuVJLGUd+NpSSnALcBn62qFxY+VlVF73jpG1TVzqqaq6q5GTY0FStJ46avlWiSGXoBelNVfacbfibJ5qo6lGQzcHiQhS11INiGk6SVGmbDadmVaJIA1wOPVNU1Cx66A7i8u385cHtzNZI0YfpZib4f+CTwYJJ93dgXgKuBW5JcAfwU+PhQKpSkMbZsiFbVD4Es8fCHBluOJE0Wz1iSpAaGqCQ1mLgvZbZrL2lQBnE5ZleiktTAEJWkBoaoJDUwRCWpgSEqSQ0mrju/FLv2kkbBlagkNTBEJamBISpJDQxRSWpgiEpSg6npzi/FyzFLGiZXopLUwBCVpAaGqCQ1MEQlqcGyjaUkZwDfBDbRu7b8zqq6LsmXgD8DftE99QtVddewCh0kTxGVNCj9dOePAZ+rqvuTvAW4L8nu7rFrq+rvh1eeJI23fq72eQg41N1/MckjwOnDLkySJsGKjokmORM4H7inG7oyyf4kNyR52xK/syPJfJL5oxxpq1aSxkzfIZrkFOA24LNV9QLwNeBdwCy9lepXF/u9qtpZVXNVNTfDhvaKJWmM9BWiSWboBehNVfUdgKp6pqpeqapXga8DW4dXpiSNp3668wGuBx6pqmsWjG/ujpcCfBQ4MJwS145de0kr1U93/v3AJ4EHk+zrxr4AXJZklt7Hnp4CPjWE+iRprPXTnf8hkEUemojPhErSMHnGkiQ1MEQlqYEhKkkNpv5LmQfBrr2kpbgSlaQGhqgkNTBEJamBISpJDQxRSWpgd76Bl2OWJttSn7xZCVeiktTAEJWkBoaoJDUwRCWpgY2lAfMUUWl9cSUqSQ0MUUlqYIhKUgNDVJIaGKKS1KCfSyafBOwFNnTPv7WqvpjkLOBm4O3AfcAnq+rlYRY7yezaS9Opn5XoEWBbVZ0HzALbk1wAfAW4tqreDTwHXDG0KiVpTC0botXzUrc5090K2Abc2o3vAi4dRoGSNM76Oiaa5IQk+4DDwG7gJ8DzVXWse8rTwOlL/O6OJPNJ5o9yZAAlS9L46CtEq+qVqpoFtgBbgXP63UFV7ayquaqam2HD6qqUpDG1ou58VT0P3A28Dzg1yWuNqS3AwcGWJknjr5/u/GnA0ap6PsmbgIvoNZXuBj5Gr0N/OXD7MAudVnbtpcnWzxeQbAZ2JTmB3sr1lqq6M8nDwM1J/hb4MXD9EOuUpLG0bIhW1X7g/EXGn6R3fFSS1i3PWJKkBoaoJDUwRCWpgd9sP6a8HLM0GVyJSlIDQ1SSGhiiktTAEJWkBjaWJoiniEqrs9TfziC4EpWkBoaoJDUwRCWpgSEqSQ0MUUlqYHd+Cti1l0bHlagkNTBEJamBISpJDQxRSWqwbIgmOSnJj5I8kOShJF/uxm9M8l9J9nW32aFXK0ljpp/u/BFgW1W9lGQG+GGSf+0e+4uqunV45amFXXuN0jDPVx8n/Vzts4CXus2Z7lbDLEqSJkVfx0STnJBkH3AY2F1V93QP/V2S/UmuTbJhid/dkWQ+yfxRjgymakkaE32FaFW9UlWzwBZga5L3AJ8HzgF+F9gI/NUSv7uzquaqam6GRXNWkibWirrzVfU8cDewvaoOVc8R4J+ArUOoT5LGWj/d+dOSnNrdfxNwEfBoks3dWIBLgQPDK1OSxlM/3fnNwK4kJ9AL3Vuq6s4kP0hyGhBgH/DnwytTg7SSrum4d/LXSwdY46uf7vx+4PxFxrcNpSJJmiCesSRJDQxRSWpgiEpSA7+UWcdl40Y6PleiktTAEJWkBoaoJDUwRCWpgSEqSQ3S+7rQNdpZ8gvgp93mO4BfrtnOR8M5Tof1MEdYH/Nc7Rx/s6pOW+yBNQ3RX9lxMl9VcyPZ+RpxjtNhPcwR1sc8hzFH385LUgNDVJIajDJEd45w32vFOU6H9TBHWB/zHPgcR3ZMVJKmgW/nJamBISpJDdY8RJNsT/KfSZ5IctVa739YktyQ5HCSAwvGNibZneTx7ufbRlljqyRnJLk7ycNJHkrymW58auaZ5KQkP0ryQDfHL3fjZyW5p3vdfjvJiaOutVV3KfQfJ7mz256qOSZ5KsmDSfYlme/GBv5aXdMQ7a7T9I/AHwLnApclOXctaxiiG4Htrxu7CthTVWcDe7rtSXYM+FxVnQtcAHy6++83TfM8AmyrqvOAWWB7kguArwDXVtW7geeAK0ZX4sB8BnhkwfY0zvGDVTW74LOhA3+trvVKdCvwRFU9WVUvAzcDl6xxDUNRVXuBZ183fAmwq7u/i95VUSdWd5ns+7v7L9L7AzydKZpndxnwl7rNme5WwDbg1m58oucIkGQL8EfAN7rtMGVzXMLAX6trHaKnAz9bsP10NzatNlXVoe7+z4FNoyxmkJKcSe8ChvcwZfPs3ubuAw4Du4GfAM9X1bHuKdPwuv0H4C+BV7vttzN9cyzg+0nuS7KjGxv4a9Vvtl8jVVVJpuLzZElOAW4DPltVL/QWMT3TMM+qegWYTXIq8F3gnNFWNFhJPgIcrqr7klw44nKG6QNVdTDJO4HdSR5d+OCgXqtrvRI9CJyxYHtLNzatnkmyGaD7eXjE9TRLMkMvQG+qqu90w1M3T4Cqeh64G3gfcGqS1xYdk/66fT/wx0meondIbRtwHdM1R6rqYPfzML3/GW5lCK/VtQ7Re4Gzuy7gicAngDvWuIa1dAdweXf/cuD2EdbSrDtudj3wSFVds+ChqZlnktO6FShJ3gRcRO/Y793Ax7qnTfQcq+rzVbWlqs6k9zf4g6r6U6ZojklOTvKW1+4DHwYOMIzXalWt6Q24GHiM3nGmv17r/Q9xXt8CDgFH6R1PuoLecaY9wOPAvwMbR11n4xw/QO84035gX3e7eJrmCbwX+HE3xwPA33TjvwX8CHgC+Bdgw6hrHdB8LwTunLY5dnN5oLs99FrWDOO16mmfktTAM5YkqYEhKkkNDFFJamCISlIDQ1SSGhiiktTAEJWkBv8PoZqtwMRcftwAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "with rasterio.open(\"warped.tif\") as src:\n",
    "    show(src.dataset_mask())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "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.6.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}