Created
February 21, 2019 03:31
-
-
Save prl900/f132738ce6f89b941a0e652f14d37827 to your computer and use it in GitHub Desktop.
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
{ | |
"cells": [ | |
{ | |
"cell_type": "code", | |
"execution_count": 8, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"from netCDF4 import Dataset\n", | |
"\n", | |
"ds= Dataset(\"/home/prl900/Downloads/WOFS_3577_16_-33_2016_summary.nc\", \"r\", format=\"NETCDF4\")\n", | |
"\n", | |
"wet_2016 = ds[\"count_wet\"][0,:]\n", | |
"\n", | |
"ds= Dataset(\"/home/prl900/Downloads/WOFS_3577_16_-33_2000_summary.nc\", \"r\", format=\"NETCDF4\")\n", | |
"\n", | |
"wet_2000 = ds[\"count_wet\"][0,:]\n" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 9, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"<matplotlib.image.AxesImage at 0x7fc07b833518>" | |
] | |
}, | |
"execution_count": 9, | |
"metadata": {}, | |
"output_type": "execute_result" | |
}, | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAARQAAAD8CAYAAAC2EFsiAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztnX+wFdWV77/r/LpX4IBwr/rgashVCMQQMAwK98bKmF9odGYMNTFlfPXCc7R8Y2KZmZd6CU6scZxxpoxViVPWS5wxJhmSimPUhAoxJkpMfJMUSiAoiD+QKwQViMAFuQeR++Oc9f7o3Yc+fbv79I/d3bv77E/VrXtOn+7du7v3Xr322muvRcwMjUajkUEh7QpoNJr8oAWKRqORhhYoGo1GGlqgaDQaaWiBotFopKEFikajkUbiAoWILiOinUQ0RERrkj6/RqOJD0rSD4WIigBeAfBxAG8A2AzgM8z8YmKV0Gg0sZG0hnIRgCFm3s3MYwAeBHBlwnXQaDQxUUr4fH0AXrd8fwPAcusORHQDgBsAoIjin0zB9ORqlxREQA49lKlQADcaaVdDE5Aajh5m5jNklJW0QCGHbS09i5nvA3AfAEynWbycPppEvZKHACpXwONj0YuSVE5kGM5PWKM0v+RH9soqK+khzxsAzrF8PxvA/oTroAyyhIASwsQHVK6kXQVNzCQtUDYDmE9E/URUAXA1gPUJ1yExCtVq2lVQiqwIvqTY8y8DaVdBOokOeZh5gohuAvA4gCKA7zDzC0nWIUkatVraVdAoTP/fPZ12FaSTtA0FzPwYgMeSPq9KKGPz0Ggkoz1lE8S0IWhhoskrWqAkiClItHEyGCreLxXrpAJaoKRAXBpKXhu5ihqdinWyQuVKKu1BC5ScoO0yGhNzdjGN9qAFSk7oVGGSV60sCo1aLbX2oAWKJtPE3XFUFVhO9SpUq6nXVwsUQVoPIu0GoPFGVc3PqV5OmknS7UsLFEFaDUfVBps3ZHSsLAr/pNuXFiiajkBGx0pT+GdFmGmBotFopKEFiqaJ21swK2/HPJOVobEWKJombo02K405z2RFqCe+ODAIVCjYwi9lg9ErLkS9UkD1N0OoHx5OuzqajODlnJgVoa60QOFGI1MRwEr9czE6tweFcUZhvK6FSQehPZUNcjfkSUs1LPXPRePwEUxMLaJRJpSOj6dSD006RBUmUdttoVrF8auWt99R0vncUFpDCQOPj6Xytqj3VFEEMD7FkNHjU7oxLdEaaLJM1PbaqNUw7eFNiZ3PjdxpKEDy480Tq5ZjbFY33pl/BqY9vAmnHTqlnWTFmKZJnzy0ldxpKEliakJT1hlvhvGVywAAxae2NrUTN+Gmx9zZxu/z67TnHElDIaI/ENHzRPQcEW0R22YR0QYi2iX+zxTbiYjuESlItxPRUh8niFK92LE3lPITWybt4/bW4fGxTAWxDvP2zMMb1w2/QsK+X/0S92bfrky/99N8sckoKygyhjwfZuYLmNm8ijUAnmTm+QCeFN8B4BMA5ou/GwDc27bkHCTD8mokfHI0wZpEI8xbNs9v5nYd0u334lNbQ5/T7/10erGFLSsocdhQrgSwVnxeC+CTlu3fY4NnAJxORLNjOL8jxQXzQh0X51vW+lDz/DbPE8XeHgDtO2TYDpv1dhBVoDCAJ4jo9yKFKACcxcwHAED8P1Nsd0pD2udZusQhT33nUOhjozxkr2Otv1kboNloZYXxy9LQSlXM59A4Fm9qFCdBFIeQUXXa+IPMvJ+IzgSwgYhe9ti3bRpSoDW3cTempO7YZk5DRzm+sHghaP8hNI61xqtwNdhWpwESneJ0fqDomM8qjWFcHOdUcsjDzPvF/4MA1gG4CMCb5lBG/D8odveVhpSZ72PmZcy8rIyu5va0gu4CEnwEtr+M+uFhzxmflv0PH2nODuTZDqE5hamVWvEy4AJAYfFC7L19MK4qhSK0QCGiqURUNT8DWAlgB4zUoqvFbqsB/ER8Xg/gs2K2ZwWAY+bQyA957lxegibrY2qNM8XenpZn67RMo1Hx7p57/2IWGmW1Ji6iDHnOArCODDtHCcADzPwLItoM4CEiug7AawCuEvs/BuByAEMATgC4NsK5c02jVtOCRCFk+5JQuTJp+OtEu9mac+7Y6PucxQXzMDZ7OgpjddDGbb6PC0pogcLMuwEscdg+DOCjDtsZwOeDnicLjkGFalW6nSLpa87CfU4L+33Z8y8DkfISU3dX4nat+s4hFHeKtgpDQ/Ij1IKivOt9Fhp5HoyeKt7nNLU0r3NHTXKeZnsxz+1l04uC8gIlLKpPlfqdTu5k0hRyKgpYP6TddnIrUFTXGjw9aDPamPNI2h00KGm3ndwKFI07523uTrsKmUUFAeNnrU5aaIEimVL/XN/7xuWB245XLzzZdp92PhB5xX5frW/8QrUaSgOQLYT8rNVJCy1QJDOxZ6/vfdNWT72Isogty3g9kyDT+W7LKlzPOzhpwjSTZEqgaEevU6gsjLLI+MplvrSysGEL/EDLFgU+xotDNw5ILc8PmRIoefaW1aRL79/vwcTUorTygr74aOM28JYdk7Z7rZI/8MVBvHmzu+v9Gfe2n96WPRuaKYHihNOD01rMZPQ98ebYxcPo+tlm199fvzXYmpmwLz77mh6vVfKzv7YRZ93j31vWifE/mR/peDuZFihu3p3mNpmdKOvDLa3ZtfKXLx1sv5OF8encoi2U+uZIq4tViCSdeqV0vMM8Zb2IK8iNW1l+Vwtr5DJ8vVxbwLU79+J4vXXq3Gm1r5Vzv9Q6fJjYd2qhvGmrMIcPQYWNXyHSro5Ow5d26WWdhllRyLRAUQX99o+P8lOzce7/fEVqmd9dMBePL5re/F6oVn116vrOIccOeuK/EfatGWw6U1qFjUzMOroJCSdnzqQzEWZaoATRDLQWkU3eNfUoerveDjxECUIQr2qnjtgoM/ru3JhYGyue2RtoBqeweKH0GSQ3iBUOBD2dZvFyMhYuh1kNq9oKWtXqo5FLsbcnERvI67cOojABjM5iTH2DfBlmadki1+HNL/mR31uCzEciMxqKV0f0SlXhh6TeLFkXJlnX8uKuf/3wcCgP45FrVgSyu5z9q7dx5tYxFN/xJ0yAVltJnPchMwLFi6gdNemOrvpKaDdk3idZTleHbhzAmzcP+vI0TeI52z2M23XeE6uWo+tY3bfdhQeXoHh8FIWxBs594GAoARbnfciFQMnalK7qK6Hjxpy1iSpUhq8fQOkEUDrBeGv+FF/HhJnudWpbfl8K7TrvlHWbPP1fnBifNQWFsToa1W7llkhkxobiRju7hNfv2qaRHkdXD2BiCnD6q+OpLnY7sWo5pr7+Nnb+dTfec334ehQWL0Rju1fSh+ThwSWu4R6pXGlGjpNpQ1E7t7FHXh6/wkDHHVFTcM5ca/h1DF8/AG/vinjpOjqOt8+ZigU3Pjs5p4sPzFAChbEG5DnuR4eWLULh+CgaaA1Ras2m0EwRI7FptB3yENF3iOggEe2wbAucv5iIVov9dxHRaqdzTcJBe2o65IiboVpHcSLt4ZjK9+hkL+HAF8Olghjc5n1dpf65GL3iQs/4L8WntmLKuk2+75F9qFN+YgvKT2xRbujBW3Y0NSbrENvsN4VqNZZ26ceG8h8ALrNtC5S/mIhmAbgNwHIYuXtuM4VQUOyJstLurH5QuUOnzZQDjMoxd93A6/luXOL8m2mbmdizF10/24wPTd+JwW1jOPDFQZxYtTxSmzE75+gVF4YuI02aL2SRVzvxINXM/F8Ajtg2B81ffCmADcx8hJmPAtiAyUIqFJ3YWQuLF6ZdBWnMXPs0eu53XxUb5vnaV9l+d8FcPL1sGirHGO/0GE3eHKpYhUsQQRPUkNqOJF0XVJzlCZq/2HdeYyK6gYi2ENGWcYxO/j0DGkncBDX+6XtmdKTThhvoXft7AED31j3N7aaA7rSg2HFo+LKNsm75i33lNQaMVKQA7gOMWZ5Jv3egRhKVPN+zo6sHMDbjlIOXl11t6qPPonDuuzA2e3qLzcNNQMu00Zn+IubqXtmL8sIQhw0yrIYSNH+xr7zGdqjgr3qypGyYcnQ6jHQ5MZtw0jJN5NRBrIb8+s4h3wZUr84W1KGs+NRWvHNGGRPTKkoIk7iGPmEFStD8xY8DWElEM4UxdqXY5gk3Gr4qY70xUTqx7LVCedYM0sL+fPvu3IhuH8tnZDg/mmXQskWhZnWmPbxJidmgOB1B/Uwb/yeApwEsIKI3RM7iOwF8nIh2Afi4+A4Y+Yt3w8hf/C0AnwMAZj4C4J8AbBZ//yi2hSbo+p04bqCOj+Kfo6v9ecVSuYJLd4y4/m6/59fu3Ive7SdRqFZxYtVyz2NkLNHg8TFXDUOlmZ+9tw/i+FXu9yOul10mPWW9Uh0kgemL0Oku9FFx0vD+8qWD+NF7z3Q5YjIj16zArP/3WnMtTKlvDuoHD+PtP/sApqzbJLW+blC5gj/cugxzb4sWjlEmb948iDO3vN30lD1044BrjNmOXG0MoOmMY/X0SxoqV9Co1bQwkYD1+Z23uRuLt5KjMDFfIE7rcE7/6QuY2Le/mQ9pYt9+8PhYU5iU+ubErjHy+BjGZtUDx52Nitd1jU0H/jgwtfndT8BqGWRSQ7GSFW9ZjTuX7hjBgbEZ2L7U2TOaurtwy7bfoK94HJ+bezEA4+XyzocWBvIHkdFWrG7sQZBxbh5cgolp5Uhrn6z1MDXtJ0a+25kaihNxBKS20i4mZ1T8rlrNasgDP1jDMTrRqNVw5yc+hc/NvbhpUGzUauh+wnnhmxthO7T1WQcVJtYZpqiUht/GaS9ECy9prUccmnamBUoSxk/7DJLsaWK/2ejyPsTavpQd74N5/810EtahrlcnHblmhWtw6913BQubEEUYyNSeG9Vu1GfPCnWs3SPY+icTtVcbt8FpXU8Sw5+oEeJkHddpBHm+0x94xnH73ofej8Y+dYf5XpizS2HaedyavInyGkqYXLJx1cMuwMKeO8trcW4eiifmR6lvTmxpUUauWdGcQp376ecx72+dhY2VOIaYXmUGaUNRNaY4JzSUFyjAqZvtlSUwiVkft/FwYIe4/Yek1Slp7pkXjzD0EwJx7+3GLMo39/7Wd7kzftuD0juM0w6NB6qP1xAz7MvLq8ykNAinc8pEeYFivdFp3HS3+kTBKTJ6np3hZFzb8PUDTT8Pc6an3bnO29yNlw6dhSnr5Hqoern3x1F2O6JMHKiylicx7DclLae2uDWgPNlR7LFig16bGXDp6OqB5p9XiAO3c7164UnMWfVioHOHwXzZWYMWWYc3ceaPsg/Frec165QkahtlLSEgZXY4WcZb7QPjTOmE8/YZv+3BsYtbtTOne9h92DCammEiVccaBqA5A3VytMUJE/A3hRzW2GpiH1bxydHQvjNhUF5DAeSPL2WplVqYOOMkCArVKlacvmfSdvs9PP6Lc9Hz4/RX44bBjIJmxTo1awqXJLQGe4hH85yFarXlt31r5Hr3qi1QmCfZUFqC6yZEnu0bUWh3X0avuBC71hrL/Bu1WtOBze24kWtWYNplu1sCKkc5f5LYh8TWttpOi4gT6u5qfjaFCp8cbdap706564/UFigWos6saOTT7hl0/Wwzpm0/FSC6ncpv9x3xKj/p4SaVKyj2Bo/P7yT0TA0hiEuEH0c06z6NWg08PhYogboMMiFQ7DcyzngObmgBFo7ZX/N+A4bVQpIWJtTdFSpvsVM9rUOjdonH7Jq5H/uLk0bUqNWA/j7wgrmgcgX71gyi1D9X+pAnM4sD7W+kJN9Q2vgqF7f7qfJ9jlK3Ut8cNEZqk4ZyQcqL8950ZPiCtISJ/dxuqDSeV5HDP31P246kqjABotVtYt/+Sblx/JRnDb3A42OhhlxJo7xAoWWLJm1L2ijrB5U7gwr0/vkr+h7ZMAVEqW+Oo7CoHzzcaugdVf/+KS9Q3MLt2RtnXHYV1QRXksi6dqdyoiZKzzqFarVpk5nYt9/RPmN3l8jCivOwqUj/gYj2EdFz4u9yy2+3iFSkO4noUsv2y8S2ISJaYz9PGKxrfOKaTu7kt2rUa/ca4jhFEFNReHutH4uCH+EgM5ZKUoRNRQoAdzPzBeLvMQAgovMBXA3gfeKYbxJRkYiKAL4BI1Xp+QA+I/aNRJQVv5r42XVXsFQTKj5Ppzoludwja4RNRerGlQAeZOZRZt4DI/r9ReJviJl3M/MYgAfFvlJwioviNgSKsq7CPnXdKXgZA71WgvsJE+BGmp1JhanqrBLFhnITEW0XQyIz8XnsqUgn7e/iEu9mSQ/SKLyc6TqpcdUPD7vaPKy+D0n7B8Vxvk56UcRBWIFyL4DzAFwA4ACAr4ntUlKRMvMyZl5WRpfTLo6ag5P9JM7G1mkN760Vo814JG4kEZMmrvMlGVcnz4RabczMb5qfiehbAB4VX71SjgZORep6fofVnfbPTt+tRAmj167sPDJ/9eR4Iio7ovnB+lLI+nWoUv9QGoqZ11iwCoA5A7QewNVE1EVE/QDmA/gdjGyB84mon4gqMAy368NW2kkTCar+qvIAssbEL9+FwW3tg0RnAT/u7FlApfq31VBEKtJLAPQS0RsAbgNwCRFdAGPY8gcA/wsAmPkFInoIwIsAJgB8npnropybYOQzLgL4DjO/ELSyTQNgd1dzPYTbcMd6jEo3POuUPvYaNqKzhnsa/2RmLY9JoVptWVyVhku+FlLRiHr/9P2XS0eu5bHGcjCxh75LupF1mmG2HX7ux6EbBzDxwcnLKdqVmxd7R97JjEAxE2IVZrTG6rTGf7DjNivTbpFVu5gTbsvE8855m7s9f7cG87FinR1qlClUwOhOu9dZJZNDHqvwsHZwrQrHx6U7RtqmDNWEJ82225FDHkDcdMuQx65J6DAD8aGiMGmnSWaJoMJE1evLlEAJsq7Cy31a1YehEl73SJXE7TKjx2cNVa8vUwIFQNOGUpjhbYRV9YZnBa/7F3UZvRbo+SVzAqVxrNby34l2DVYLm2CE1Uisx+35l4FmBHyZ99+tblpopUPmBArQ3oAlK46HbpQGYTUS63GVEWqJgB8Vc6bOrW76pZEOamcO9MDNhhLFWm4N1ORUdlyoPDslq25u+V9K/XMxsWdv4PKcIpwlmSFP40zmNBTq7nLs8DKC36TVqe2h/tLE7rsT98xZGGHidk4tTNIncwLFdHALsrLYDyp05rQilrl5ofq9p2kIYlU1uqwhu51lTqAAcufsVY3bmVR9Sv1z255LBWGbRnK3TkB2O8ucQAnTqLx8VVQTJECyHdhryFGoVrF4K2H2b5xd6uPGrjmp+KySRnWhmjmBIhNVG2iSxmAvGrUati9l9FaOJ1IfTXtUbbMmmRMoSfgwdALF3h7f93L70nTWe6neeTSTyZRA8aPu+VUJszTFGFTN9eM275X4e/j6AZxYtTzUuTWtdNr9y5RA8fPG8rNPsbdHeWFidnyZNiMAk1ZqA0BxwbzmtuHrB9AoA1PWbWpbVpxkIY+vHzpNy8qUQImCdXjj9XZWBT45in1rBpsNMmwH8xrWmWXXdw41t/Xc/7RjVr+kycIz0kzGTyrSc4jo10T0EhG9QERfENtnEdEGItol/s8U24mI7hEpR7cT0VJLWavF/ruIaHUcF+TW8bzCHsRNUGFA5Qr+eOOp8BQ8Pha6g6muiWnyhR8NZQLAF5n5vQBWAPi8SCO6BsCTzDwfwJPiO2CkG50v/m6AkcMHRDQLRoDr5TAyCd5mSRDmCz+CwK3jOWUXTAo/wmDkmhXNzzw+hsI4T3JXDxPdP6zhudPG/llD1efjJxXpAWbeKj7XALwEI+vflQDWit3WAvik+HwlgO+xwTMAThdpNy4FsIGZjzDzUQAb4JwzORFUG9tOf6A1bafTsKOZEN4l1KKVvQ+9H4CzhmIVMpfuGHE8XrX7o2lF1ecTyIZCRO8G8AEAmwCcxcwHAEPoADhT7BYpHak9FWmQ4MTt3t6qSnU33OrrZxgz99PPu/5mDv9uHnrZNRJb0veqk6fw84Tv1cZENA3AjwD8DTOPEDllFzV2ddjmOx0pM98H4D7AiCnrVxL78XqNQ6rH6W0bZ7kj16zAPfO890kSbevJB740FCIqwxAmP2DmH4vNb5oZBMX/g2K7WzpSrzSlkUjThV7Geb1yMpv+IEGxD2XqlyxtKfudXu9HH0ZDyZoGqJGPn1keAvBtAC8x89ctP60HYM7UrAbwE8v2z4rZnhUAjokh0eMAVhLRTGGMXSm2hebAF72Td2cFr5XTU9ZtwtHVA83vpb45vsq0D2XM1BVm2Wfd02rwPXTjqXOEFdBOx+ihTGfRNo0GEV0M4DcAngfQEJv/DoYd5SEA7wLwGoCrmPmIEED/F4bB9QSAa5l5iyjrr8SxAPDPzPxdr3M7pdFo1svW6Iu9Pdp3wYFS3xzUDx4G4K5NfWX3c/jpWx/A88vLjvuouohSIweZaTQyl5cHSFZ4xNmZvMoeuWZFc+bnvM3dmFocxfPLywDUtfBbKfXPBeC8mlkLKLWQKVAyGQIySU0kzobvVbZ1GvnVC08G6oRxJD8r9c3BxL5TJi/TWc/+LMzzeYVF0MIkv+TW9T5vBsIg9gnrvq/cv6jFBhOW+sHDk5YvOAl2LSw6m0xqKH6Q1bCTUM/DnsPPVOv81VuN8sNUzIJTgKMw9dbDnXyjNRSXY5MMDRk2SPVX92wKVL5swtRbC5N8k1uBErbhFqrV1MINBj3nD46uaL+TIM4hoBYSGpNcCRQZncYpXohsaNkiKeXYI6k5+ajM+G26cUXyZsvSeJMpG0rcGQPjLq9Z7pYdbaeMS+8wJk6jltmedtdvnYUxOXaxYTjNsiexJjtkSqComsnPvn9Uw6N95bGJV5lu59y/7nzMWfVi6LpoNEHIxJAnabU5qDBwcp1v53JuPebNm6MvIZhUh8ElANAiTPatycdShTTRQzhv1BYoRIlOM8pKJlXqmwPqMsoZ3DbWtG28cr+zM6J9XQ1geMc61c8vYzMrzQWBALD7rgFMnDZ5v+HrB8CDS3DoxgHf64Q6mSwN4dIQfsoPefzEQJH1kGWVM7FvP3hwCejwMDYuqcBcVH3aH/w/4FcvPBmpft1PbGvZ/9wvGQGb7Per5/6n8fqtgzjnjo2Y8F26JgukIfzU1lAs64zcpK1qOY3NNSy0cdukss+5Y7Im4oXfutkTnAOT74upfTndL7NeSUWaz/KwIS/R+ONCbYFiIazzl99yo1Dqm4O9D70fc56pYur3jxvaSbniK1SjjLpZY+San53uU7vy9Grt9uh75E1mBIqJCmNYu8H11Rvejbmffh63zv45bprzJGijMdwIEoVs7+2DeP3WQWl2DPM+mXW13zf7NSQZt0SFZ2glyxqTamQufEGaa0HandsMM9AudadXOeZvpf65nit2g2CWeemOkZbAS+b2wz99D3r//BUp50oLUyPUoSSD09HhC9J8u/H4GIoL5oGOn2hxIjOX9jsZUt3KAZwFC48bs0KyhIlZ5uC2MQy9cyaAU6EQzHNnUZjYbUVpLZfQtJK5IU+aFHt7wPvfRP3g4WaDpnLF0UPVD27+KmHL82Ljkgp2Hjuzed6sYwoQ6/0Lk7dIIxflBYoKjcOaYLxRq6F49my89uVlzYWEUWjUaon5f5Q+9loi55GB3+duDnFMrYu6u3IhMLNKlFSk/0BE+4joOfF3ueWYW0Qq0p1EdKll+2Vi2xARrXE6nx0VGod9XD6xZy/OuWPjpO3WThDEyOlHI/GKjO8XFYSzX6IsfdCkR5RUpABwNzNfIP4eAwDx29UA3gcjUPU3iahIREUA34CRqvR8AJ+xlJMLrI3a+ua0d+QoGomZsjRLHej4VeFSgfjBLkysOaw1yRMlFakbVwJ4kJlHmXkPgCEYuYwvAjDEzLuZeQzAg2Lf3GHVTpyMhef/9EDgMs0yxqcQFm81/oDkghtF0W5GZ8QzsnYzamvSI0oqUgC4iYi2E9F3LInPpaYizQqD28aasUfaTV22m1Z2g8oV9Nz/NLYvZez402kAnDvQyDUrlHIA7Ll/cp5mGWjhoR6+BYo9FSmAewGcB+ACAAcAfM3c1eHwQKlImXkZMy8rI5qnqWy8OukLtdkY+XAyPhBUrrQILbu9ZvoDz7jm1zGXBvghD0m6smQ3ygOhU5Ey85vMXGfmBoBvwRjSACmkIk0Kp/UxJscuHnaNTC+rUZtu9VY/Fqey3Ww07dJb2GnUapGmYa3ZCNNCazHJEjoVqZnXWLAKwA7xeT2Aq4moi4j6AcwH8DsAmwHMJ6J+IqrAMNyul3MZ6eDl7WrSqNXA42PYt2YwcjwSKldQmOFtnwHk+rFEcRg74954hjoy0RqMXPxoKB8E8D8AfMQ2RXwXET1PRNsBfBjA3wIAM78AI0XpiwB+AeDzQpOZAHATjHzGLwF4SOyrLGFVfnv0NgDouzPYSmO3cp0WpznZbLxWxQbpRHnvcFqDkUvm1vKoTqFaneRslRZpn1+TDWSu5VHeUzZNgmoodmNpGp25JTZKxPAJUQkS2tL0r9FR47KNFigOmIIk6MpVuwAZ3DbmOxlXELyGIU7OdXGf0w2n0JZumIG5w9h/nAJMadJBCxQHZHXEn339T7GvPgOAe6MP0wGykIQsSVo8ZRWuZyegBUqMzFx7apbDrdHrDqDJE1qgWIhj/H7PvIVG2QEcylQlD45umnjRAsUCj8p39TeHNDIDJsXNoRsHmkZSKzoamqYdmYvYFieyAhBbp2tlD2mSmArOgkOaiZ4aVwutoViQ6SKvKmZGwbyg8r3uRDpSoLgFK0q7cfqxUUSt4ztnRctIqJmMvn+nyI1ACRsXpNjb4+oAFldDKVSrTdd4q4t8EjaKqY8+O2lb2oI062Tp/sUt/HJjQwn7UN3sJn5SXYSlUasBNe/zx0USjT9tu0ba51eZuO9LbjQU2cQlTNLg9VujrXIOStr3x+v8b948qIcoMZJbgRImL7AfrGk/2xHWr2XX2qWhjnNi+PoBTPmjugtAk+b0V8cx8cFFaVcjt+RWoPh5SxYWLwz1NvV7jJ91KSdWnQrgvPuuAey+awDzV28NXCc3eu5/GoXBLF4fAAARL0lEQVRxacUlQmHxwtjKPu2/XkbxKeP+xqGpdLr2k1uBYsXtITe2v5xwTU4xcs0KjF5xYcu2c7/0NM79krMPSDttxyuyvHUJQBaoT+sCLYtHi5CxGtz6ErCT9nAvCHEkResIgZLEQ/YKaORGvVJAcazha9+JffsxvvJUyAp7QyhMBD69stDGbeAtO059t1yrXQinwZR18leQp4XsvpGbWZ4gFHt7mrMrZl5iN/waYYPO1pjL9YNQfmJL87O1TsevWp4rgWKnmUJk5TKU3q6jsHhhqtplXojjRduRAsXs/H7yEvu96WnO/kx7OD9vTC9MgepPp8smWZxFtOInSHU3Ef2OiLaJVKS3i+39RLSJiHYR0Q9F4GmI4NQ/FOlGN4lcPmZZjilK48DP2DBouktZZWniIQ8GUbdZxKwsmfBjQxkF8BFmXgIjB89lRLQCwFdhpCKdD+AogOvE/tcBOMrM8wDcLfZzTVEq82KsyO7gWmCoT1JOe3HjmFNp47bYzysDP6lImZmPi69l8ccAPgLgEbF9LYBPis9Xiu8Qv39UpOJwS1GqPHl483UChcULY5sdMtEvFm/8JvoqEtFzAA4C2ADgVQBvidQYQGta0WbKUfH7MQA9UCQVqTm92k5IWKdp3bLwAckGHTp+1XLUL1mK/etylWPekTBCvD6tC41KbEpvLkll2ljk1bkARra/iwC812k38V/pVKTTHt6EQzcOtH3T+DXWJh10qFEp4O2RySuG80YYTYA2bmsODYavH1Aic2GnEWiWh5nfIqKnAKwAcDoRlYQWYk0raqYcfYOISgBmADgChVKRnnHv04ZTT3dXpqKQmbM5859IuSIZIK4E7XlD9hDOzyzPGUR0uvh8GoCPwcj892sAnxK7rQbwE/F5vfgO8fuv2Mgm5paiNBWchIm2lWjql8hbR9WJ+NFQZgNYK2ZkCjBSiD5KRC8CeJCI7gDwLIz8xxD/v09EQzA0k6sBI0UpEZkpSicgUpTKvBj7HL7Vgc2KNbufSTsHt7jIut9B3jDX+WjCketUpKZAKfXPRf2NA7rjdgjHr1qORplCeSN3EubLTKciDcjEnr0twsQ6gxNmDU7SuC38G7lmBU6sWu65MLATmfbwJhTG431R8uAS0LJFjulW4566lkFcmnGmXe/dboq53Wm4Y3e3TzpimiyOX7UcXcfqqFcKaJQJxd4eNI7VQjeStIZ8cRH3cgTauA0M4Kwtk38rDtfAC+ahvnMo1jpEIS5tPbdDnk6zTZT65oBHRzMrIDXx4dQXzG1UrmDD2AN6yONFoVpF8ezZaVcjcXi0cwSoxj9OL9a4XraZHvK40ajVMuVf4oe2GleljML0KvjkaEdpZppoJO6HkhXqlyzNhIE1LG0f/Ng4GiM115QgUelEH52woSjdlmN0wj3MhUAp9c+NNU5oJrBoKHEQ9U2Wxfit47OmhDrOTTvuBM0xF0Oe+hsHACDSLEfW4dpxpQ2ypgFQ5vOJ+1lrJ7fg5EJDiSsxeVagckVpYWJiPp8kV2jHTVYCHyVFLgQK0FlDHbNDmv+zJkj55Ghunpc98FGUtUC0bFEmnOK8yLxAKfXPBZUrKJ49O3BaAOu+WWrg5hjd2jFLfXMy8+bn8TFlhWBQw7693UxMDR+PhbfsAI1NoLhgXugy0ibzNpSJPXtb/gfBbNRZdIIz61yoVsHjY5nwcjU7n8r3OujQ0X4t41MKkaL4ZD2af+Y1FBlTxSo38HaY2oqXhlXs7YlNewmi2VkFeJY0wiB0SgYCNzItUArVKhrH8uXA1g63zmjvrFSutIaxFMMj2YLFjzA262tqVVkW4BpvlBcobm+yYm8PqCt7Q5WomB3S7bp5fAyFGYbQMIdB9cPDzWPi9iB2Enj2Wbi8aicmr986eQWyFWsGSB5cEmsu56RR3obi1HHcAidpDGTcGypXUDyzF/WDhwMJbb/75vlFcM4dGz1/t2aALB4fzbzdxIryGooTYTrMeZvzH9g5CvahkGnojaPju5WZd83FiTwJEyCjAsXa8PwaZV+98GRc1cksVnuL36EQlSvGcDOGzt8pQyK/ZNEvJUoq0v8goj1E9Jz4u0BsJyK6R6Qc3U5ESy1lrRapS3cR0Wq3cwYhT0OfpDtSO3uMHdOoGtYQHodROM/wlh0AsiVY/NhQzFSkx4moDOC3RPRz8dv/YeZHbPt/AkZE+/kAlgO4F8ByIpoF4DYAy2Dk4/k9Ea1n5qNBK53X8bfq1+W2xMGvH0+W/X7SgrfsMATK2MSkJFZFBaPCRUlF6saVAL4njnsGRv6e2QAuBbCBmY8IIbIBRo5jX2g1OBnC3OegRtt2QxstbFrhLTscbS10/IRy/SJUKlJmNr13/lkMa+4mItNB0C3lqK9UpG7oRpYMUe9z0CFN0h0izQ7Ig0uk5v3hUfWCaYVKRUpEiwDcAmAhgAsBzALwZbF7pFSkbrmNVZPEeSXqfQ7i52JqK7KerZ9yrOdLsk3x4BI0KkUUxuSlolLRfhholoeZ3wLwFIDLmPmAGNaMAvgujJzHgHvKUV+pSN1yG6smiWWQNyEZ9nrcMhcELS+oD4wMbcxvHWnjNhSf2orCWF26gVWldhQ2FenLwi4CIiIAnwSwQxyyHsBnxWzPCgDHmPkAgMcBrCSimUQ0E8BKsa0tKt0wmeRNSMq+Hh4fUzqsZ6MWPKAXb9lxavZGUrtWqR1FSUX6KyI6A8ZQ5jkAfy32fwzA5QCGAJwAcC0AMPMRIvonAJvFfv/IzEf8VJK6u0LPLKiKNY2BjOuIUo71WNXuaf3wsJLP2qtOfj25VbsmGeQ2L49GXWjZouZbOm3SCKlQWLwQ47OmOIaY5MElk4I2xY1ORaqRRju12/SMVaUuMoe/aQ2luVJC6fhkAVa/ZGniwkQ2WqCkTNr2oXZvZreUrpHOuWVHS2iFIHUBwgsW++xOUE9hWZh2FHs82jwExdYCJWVkTpuGJY3zm6EVnPxW7FqR9bNVAPjxebHHYjHLUIGsayNOKC9QrLEj8kqSDdwrOJMqmOuFTE2kfngYpb45k4RDo1Zrq62oJkTcSPulIgvlBYo1dkTWSarR2M9j/e5HI0rU4atNBPymIKiUm3U3A0iV+uY0MyXKqHNc1+3HBqW6wPOL8gLFStaleBKNxskL1H5ev7aKOGk3u2L3ojWTuZk2HSpXMLFvv6GldHdNSsEa95qkIFBXPOlhVSRTAiUvUlwmTuEW/eQ3Tls4W5+ll03EOmSxC0kz+HajVpuk6agUWyXJjARpJx7LlEDpNPyuTTH3NfcP6sGpQqcD2q8Dsl9T/fAwCtOrKPXNcV0XpILRO0lKw2+nmtdH+ZiynUzQsABhUlrYj01bCwzqFWt9+7t5qMq4JlMTsqOaF29951CqAlRrKAkS94OOGqu1XcdIoqFG8TWxChNz4Z6Tv0sY3LSnuITJiVXLQ0fDT1PAaYGSIGk96CDR1KLuIwtzCFOoVgOt6jUxh32NkfaJ0FRkyrpNeOec7IXL1AIlA/jpDGE7jOodrVGrGTM57z23RdvwK2SamRXbzAJFdeunckVq8CQA6PrZ5vY7KYbyAkX1Bp8EfoYibvu08yb1KlsFf5XignnG9b20u8VewidHPY6aTGF6qwCyXrc1o2Gpf+6kY/0E1+bxMRSf2pr6LEvaKC9QVDJ4qYrXPYqSKTCsv4pMQVPfOeS63iZI22iM1FwN11a7zcSevS0BkKhcCeQ8Z7rTq/YiTCpqvvICBVDbC1IFglxbkobVuLDGkvGLXaMZveLCZlnNfUyfFxFawdRKzGFX0Fk31SgumBd7GpNMCBQZD0fFB+yHIL4o7faPS5gkLazbrc8xO4197Y/1e/cT2zzLAE4ZdmV0wrRfaLxlB07OPR3o9x0XPhSZEChupP2QksApUl2Q/e2/xSFYVRPW5jDPLnjs/01DrNUgaxUepuHXa9joN+B1mvfIdHQrjDViT32aaYEi+yE5NQp7g0sb1TpvEsi8907ajfWzuYLZ/OynPNWc2+w0dr8GIJl4K5kTKHGq7W6Gv6iBfTTRSMImY/0PtGoqKmsffkg0vKXfHUWyr2eJ6FHxvZ+INok8xT8koorY3iW+D4nf320p4xaxfScRXRqmwk43R0ZH9zOjkUZ0LzuyPD+zjizB7vRMzVgrbr9r3AmioXwBwEuW718FcDczzwdwFMB1Yvt1AI4y8zwAd4v9QETnA7gawPtgpCD9poikL4WgtgYnrIJJVU0kyZWrbqhwb2Qb6rXQkIPfVKRnA7gCwP3iOwH4CAAzUfpaGLl5ACO38Vrx+REAHxX7XwngQWYeZeY9MNJsmMnBImN3Pgq72jYrEb7SxM+9UVWTsj5r68sjDTuICoJZNn5XG/8rgC8BMHtsD4C3mHlCfLfmKW7mMGbmCSI6JvbvA/CMpUzH3MZEdAOAG8TX0V/yI+3zLYzZ/gcl6vHh6AVwONEzJoNxXW+kXQ0Xxhw++3/ucp+ZOu+sBbIKaitQiOjPABxk5t8T0SXmZodduc1vvnIbM/N9AO4T594iK1+IauT12vJ6XUB+r42IpMVZ9aOhfBDAXxDR5QC6AUyHobGcTkQloaVY8xSbOYzfIKISgBkAjsBnbmONRpNd2tpQmPkWZj6bmd8Nw6j6K2b+7wB+DeBTYrfVAH4iPq8X3yF+/xUb6QnXA7hazAL1A5gP4HfSrkSj0aROlIhtXwbwIBHdAeBZAN8W278N4PtENARDM7kaAJj5BSJ6CMCLACYAfJ6Z623OcV+E+qlOXq8tr9cF5PfapF2X0rmNNRpNtsicp6xGo1EXLVA0Go00lBUoRHSZcNEfIqI1adenHUT0HSI6SEQ7LNtmEdEGsTxhAxHNFNuJiO4R17adiJZajlkt9t9FRKudzpUkRHQOEf2aiF4ioheI6Atiex6urZuIfkdE28S13S62p7KsRDapLJdhZuX+ABQBvArgXAAVANsAnJ92vdrU+UMAlgLYYdl2F4A14vMaAF8Vny8H8HMYvjkrAGwS22cB2C3+zxSfZ6Z8XbMBLBWfqwBeAXB+Tq6NAEwTn8sANok6PwTgarH93wDcKD5/DsC/ic9XA/ih+Hy+aKNdAPpF2y0q0Cb/N4AHADwqvsd+XalesMeNGADwuOX7LQBuSbtePur9bptA2Qlgtvg8G8BO8fnfAXzGvh+AzwD4d8v2lv1U+IPhHvDxvF0bgCkAtgJYDsMbtmRviwAeBzAgPpfEfmRvn9b9UryeswE8CWOJzKOinrFfl6pDnqb7vsDRTT8DnMXMBwBA/D9TbHe7PqWvW6jCH4DxJs/FtYlhwXMADgLYAOMt7GtZCQDrshLVrs1cLtMQ330vl0GE61JVoPhy088wkZYnpAERTQPwIwB/w8wjXrs6bFP22pi5zswXwHijXwTgvU67if+ZuDbrchnrZoddpV+XqgIlL276bxLRbAAQ/w+K7W7Xp+R1E1EZhjD5ATP/WGzOxbWZMPNbAJ6CYUM5XSwbAZyXlUDxZSXmcpk/AHgQxrCnuVxG7BPLdakqUDYDmC+s0hUYhqL1KdcpDNZlCPblCZ8VMyIrABwTw4bHAawkopli1mSl2JYaIvTEtwG8xMxft/yUh2s7g4hOF59PA/AxGDF/Mr2shNNcLpO2MczDqHQ5jBmFVwF8Je36+KjvfwI4AGAchmS/DsY49EkAu8T/WWJfAvANcW3PA1hmKeevYMSKGQJwrQLXdTEMNXc7gOfE3+U5ubbFMJaNbAewA8Dfi+3nio4zBOBhAF1ie7f4PiR+P9dS1lfENe8E8Im0r81Sr0twapYn9uvSrvcajUYaqg55NBpNBtECRaPRSEMLFI1GIw0tUDQajTS0QNFoNNLQAkWj0UhDCxSNRiON/w8IUl+6bH7iUwAAAABJRU5ErkJggg==\n", | |
"text/plain": [ | |
"<Figure size 432x288 with 1 Axes>" | |
] | |
}, | |
"metadata": { | |
"needs_background": "light" | |
}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"%matplotlib inline\n", | |
"\n", | |
"import matplotlib.pyplot as plt\n", | |
"\n", | |
"plt.imshow(wet_2016)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 24, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"mask_2000 = (wet_2000 >= 1)\n", | |
"mask_2016 = (wet_2016 >= 1)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 25, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"<matplotlib.image.AxesImage at 0x7fc07b4480b8>" | |
] | |
}, | |
"execution_count": 25, | |
"metadata": {}, | |
"output_type": "execute_result" | |
}, | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAARQAAAD8CAYAAAC2EFsiAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAHgRJREFUeJztnX/sXlV9x1+ffcsPmSgtCilFQ9WyiQuryGiNZkHUFplZMdGsbJEGSbpMyHQum60mQzf5QxNlISoKA4VFLQx1NATtKmK2/cGPorUUaulXYA7a0LgW1JAUwc/+uOept7f3Ps/9ce6959z7eSVPnnvPc597zzn3nPf5fT6iqhiGYfjgd/r2gGEYw8EExTAMb5igGIbhDRMUwzC8YYJiGIY3TFAMw/BG54IiIheKyG4RmReRDV0/3zCM9pAu56GIyBzwKPBO4EngAeASVX2kM08YhtEaXddQzgPmVfUxVX0e2ASs6dgPhmG0xIKOn7cE+N/U+ZPAivQFIrIeWA8wx9ybTuBl3fmuI848+zke3XFC394wDAB+ycGfq+orfdyra0GRHLcj2lyqej1wPcDLZJGukLd34a9ueQhWCGzZu53Vpy1vfDtf9zHGyff09v/xda+umzxPAq9KnZ8O7O3YD8HgSwRiEZMte7f37QWjZboWlAeAZSKyVESOBdYCmzv2Q2dYBjqSWISvK4aYPjpt8qjqCyJyJbAFmANuUtWHu/RDl1gGMqYxxPTR+TwUVb1LVc9U1deq6tVdPz8EhlgyGQbYTNlOmQjJEEsmwwATlE6ZCInVUKoRYnyF6KcQMEHpgbZqKENN5CHW6EL0U5ote7f3kh5MUAaCzUUxJvTZtDZBGQhjFZOh1sqasPq05b2lBxMUI2razjihClaev/pq5qQxQXH09SL6TgDGdEKt+eX5K69m0nX6MkFx9JVwQk2wQ8NHxopR/LtOXyYoxijwkbH6FP9YxMwExTAMb5igGIcpKgVjKR2HTCxNYxMU4zBFiTaWxDxkYhH1oAXlzLOf69sLtQhh+M6Ij2lpJhZRD1pQYtsm0YRkvNh7TwhaUOpg80mMPmhag2iafqr+v6302qkZjarU3VO2j3UttjWBESvf09sfVNVzfdxrcDUU6D5TW+3E8MEQ0tEgBaUr0rWSqiI2hMQzZsq+v7G950aCIiJPiMhDIrJdRLY5t0UislVE9rjvhc5dRORaZ4J0h4icM+v+oY/yZEUkT1SKEtTq05ZHldjq+DWm8FWlbAFSpaCZda3P+Gzr3fioobxNVZen2mAbgLtVdRlwtzsHeBewzH3WA9fNunFsozx5DKVPpU44hhL2OrSRYX3GZ1vvpo0mzxrgZnd8M3Bxyv0WTbgXOElEFrfwfK+0WcqmX+qQS/MhUfY91c2wsaeDpoKiwH+IyIPOhCjAqaq6D8B9n+Lc88yQLpl281CaPE1e8rT/pn8rai6NdZVsaHQVh1WazU1oKzxN7fK8RVX3isgpwFYR+cmUa2eaIYUjbRsfzwm8Ne9fHdK0ryP7/3SCmVWK+epnGXPTwxeTOOwjLtt4ZpBNHlXd6773A98GzgOenjRl3Pd+d3kpM6Sqer2qnquq5x7DcYfd+5yF2jTyJ6NARfcp2n3Lx7ONOJiWBqb9J7TaZ21BEZHfFZETJ8fAKmAniWnRde6ydcAd7ngzcKkb7VkJPDtpGpWhz30y22ZauEJLMIYfsmIwlLTdpMlzKvBtEZnc5+uq+l0ReQC4TUQuB34GvM9dfxdwETAPPAdc1uDZgye2YeUh43vmddnaZ9Pf+yD4qffP7Ds5yIhLMwQTFkMIQyz0Gdd5z/Y59T54QamzlscYBqFlvKFia3lKEHpzoexw8pjpM0PHKiZ9px2roRjGFMZQU7EaimH0RN81gNAxQfFMlQTX1gxcH4w142TDXWYBaNV7DhkTFM/4XF3aJyH7rU18hbvqHJOhiE5UghLizMC+GGuGbwvfaSuEmkwfeaXpWp5OsUxkxELVzty+1uts2budOY9r/qOqoeRRZw3EGLE4mY7vpR2xbF/gW8iiFpSiUmDi5nuHq5gzpdXumtPW+x/Smp6oBaXLtQ5VVwsb46VqeiibTsusPi77HzOjYRgtUaW/I4SJbr79YBPbHF3N+TCGja/O0z53dZtGl831aAQlL0Kq7BLed6kCJmpGP3SZ9qMRlDolg6+2qS9CELUmxC6Ibfu/yR42Xfiti2dZH0oPhNAOHwqh1UKzxPCurQ8lQ2xDuqEnsFjo6p03mes0tncdfQ1lVgkw7fcYSg+jfca4IXg6zKOpoUyzy+OjhBhLAoqp9tYHQ08H6feftcftO23MFBQRuUlE9ovIzpRbZfvFIrLOXb9HRNblPStLninSdITEUsPoO0PHEEdtUKUpHFoHvk+ytqAm8dJGWMrUUL4KXJhxq2S/WEQWAVcBK0hs91w1EaGq5EVO6Iw1Q/ugqZE135b4Yn+Xs/Z7acpMQVHV/wQOZJyr2i9eDWxV1QOqehDYytEiVYvYX7AxHZ/7k6RL5ex39rhrupy60GaeqduHUtV+cWm7xiKyXkS2ici2X3PoqN9jqJGEhsVZsTCla7lj2xS7lz6UihTZLy5l1xiKTZFOsBpJdcYUZ7MySJ6N4rEu+myjD7KuoFS1X1zKrnGWaaM8aXy9+Dr3MXMYYTGtz6SqwfMQ1u20RVtNn7qCUtV+8RZglYgsdJ2xq5zbVPJGefLwNa24zqIrG5buliZT25uKQLofJuZ32+ZE0JkT20TkG8D5wCuAp0lGa/4duA14Nc5+saoekMTQ8edJOlyfAy5T1W3uPh8APuZue7WqfmWW56ZNbKv6UrtMBLEnuNgoM/3e3kkxozdF2vbQ1yxC6MQbAr4yefo+fbyb2MVqNDNls6Srm20Pf03zQ1/PHho+JpNN3kd2Bmj695BXGTchxH6cKGsoaWIvHYzpFNU46tREfKSVuvcIJZ3m1ebmFs+Ps4aSRxsbUqdpe0/OsvcJsTTqirwaR53M6WMn+jod902e7ZvscLlvf0VdQ8mqbdsvbdY6othLr77xGQ+hZeS+yauZTLAaiqOvdT1Nd4jz9b+xUeX9Drmfq046b7smPyF4QQmlqp8tPZvMfQklTCFRpnbSpMkSwjvyNQmyiVC2PaARvKBA/kKu7G/QfklfVI0eaknYJaHFYRuzZMtMguyysAlppmxn5PWR9FnC+3gJvpfUh06XYct7lu8Suc33V7cProp70+dNI3hByUZKX5Pa2m6Th1ZC+6TuyEgXz/JBei5M020R6jTNiprffcyPCVpQ0osDfSaUPhcTGkcyhDgs2gIh61amL6dqOi/T/O4yjoMWlAm+mzq+qpVDrlV0xdDjMC0ik1puFxk8K17pmlOviwP7JD0PJa9q1+ViPxh+4q+K73k3Vd+xzd8ppkpNaPSLA7vGBMUPvieudb0AEPxscVF1ZLJsv2HRfWf936egLPBxk7bJa4umz7vAxKQfioQjBjEp8x+ftbCiGdxdxlVUfShF50Yc1OkLC+Fd+9piIe++ZSfzlY2zvuMryiZPiG3nEP0UEtZsrEcX69VGtR9KkbKHNtxoGWU6Q15bU5dZaTjGQip4QSmK0LyOrq7XXwwdm6/TLlWbOzGIS11TpJ8QkadEZLv7XJT6baMzRbpbRFan3C90bvMisiH7nDqkx9bbGuOP4SW2RdOwV80IIQrPrPVjXTw7pjRY1xQpwDWqutx97gIQkbOAtcAb3H++KCJzIjIHfIHEVOlZwCXu2kY0WfFrtE+TWZ+hvM8+R5hiEpIJdU2RFrEG2KSqh1T1cWCexJbxecC8qj6mqs8Dm9y1XsjbF6WoCdRkXUXd9RmxU2bZ/ZBqhmN6t75pMg/lShG5FNgG/K2zWbwEuDd1TdrkaNYU6Yq8m4rIehJD6xzPbLs8VcfdqyTUafeNsfSoS5n4LFrP0iZtPM/EpBl1O2WvA14LLAf2AZ917q2bIoX8mkNe/4nvxDHWGkpZuh7J8fm87Hobox61aiiq+vTkWERuAO50p9NMjlY2RVpEUWlYpaZSZ0hurDWUImIc1kzT5eZcbRLSe6glKCKy2JkYBXgPMBkB2gx8XUQ+B5wGLAPuJ6mhLBORpcBTJB23f17X09NqIlWmKRvNiD0OY/f/hJDCMVNQ0qZIReRJElOk54vIcpJmyxPAXwKo6sMichvwCPACcIWqvujucyWwBZgDblLVh6t6tqiTddpwcUjqbRhDJ7qp99NWTnYlHiZSzWgafxb/fhnV1PsJszY46iORWcfskZSJjzpxlp0CYGISLtEIStEozuQza+u7KiM0s+ZdpIdKx0SfO+aNLa5jJcomT1HNxKrCRqz0mXZH2eQBPzvgWzNlOJSZwRsLdXaCC5GoBCVvrUeZ7fCy9wj1ZYREDJnV14zoGAk1fFEJSpYxJ6i2aTNuQxEkwz9RC0oRVXb8NmZTVwCm1RJ9MW3+kdE9UQrKrA4sXyWoJcqEuvHZpnDPavJaodEPUQoKFO/Q1kQEsom061WzIdK23+ref9Y0AaMfohSUvAzvYyl7X6Vand3g2yI7X6ftkbO6cT5roqPRD9HNQ4F2xuxDm8PSpX9CC7vRHVv2bmdu8fw456FM8DlmH+q+nSGJSSg1pxD8MTR8p7PoBKVOopo2VyU0IYFuM/CsPWP6zMTZ9TshvivjSKITFJ+EmkBD6Qy2TGxUJTpBsf1D/VCldtZ3Z7URD1EJis/l8aE2d/KoKnxN+4zSTZ0xi64PxhZ/UQlKGQEoO8wZupg0ydBlJ/1Ne0bfWzQMJSOGns58E5WgNCHWDXqazlGZ9r8i0QghfkLwg1GdMqZIXyUi94jILhF5WEQ+5NwXichWEdnjvhc6dxGRa53J0R0ick7qXuvc9XtEZF0bAWpr17Am1GmyZHfYD3H6u2FkKVNDeYHEkNfrgZXAFc6M6AbgblVdBtztziExN7rMfdaT2PBBRBaRbHC9gsSS4FUTESpLmYxZZm1H11sYNDHXkSa9Q11Z2qjZGP0T6vspY4p0n6r+0B3/EthFYg1wDXCzu+xm4GJ3vAa4RRPuBU4SkcXAamCrqh5wVga3km8zuRNiLLnr1FTqrnmJMX7GRKjvp1IfioicAbwRuA84dWKbx32f4i5bwtFmR5dMcc8+Y72IbBORbQuX/LJS38es0jtUVS+ii+X/VZ89lOcZ7VBaUETkpcA3gQ+r6i+mXZrjVtocadoU6cGnTiydeSb9Dl1vutRmRmhLOMo0+bouAUMtcY1qlBIUETmGREy+pqrfcs5Pu6YM7nu/cy8yRzrNTGkj+hwG9vHcWXvl+qTs2qW65i6McVNmlEeAG4Fdqvq51E+bgclIzTrgjpT7pW60ZyXwrGsSbQFWichC1xm7yrmNnirDtk0zbZvzdGyPEmPm9gUi8lbgv4CHgN8454+R9KPcBrwa+BnwPlU94ATo8yQdrs8Bl6nqNnevD7j/Alytql+Z9uyi7QsgjslpIVC1D6pIFCyuh4tPMxq2H8oM2nxW1XuHutVCHtP8agIVFqMXFGM6bRg/yzOwBkcLholFfIzW0FcVxtB2LzOs3Nbq7KIRNROTcTNYQfGVsLsQJp8bNfu8f/ZZebWRqoxB6MfMYAWl6e73XfZXtL1JdZvzWaCav60GM2ysDyXDUPsAhhouoznWh1KAr6q9z/vl0eV9+25i9P18o1uiqqEMqZStE5Yhhd8Ih9HWUELZvHnW9X30hVhNwAiBKASl68zSdIuAPhbfWc2lG0y4pxO0oJx59nOdz4r1kWBi63vp6v5DICbh7uN9Bi0o0M6q2GnP8pFg2h6mLUNRvJTZY9YYBn28z6AF5dEdJxw+bmuzoTb7O5reu4pJkOxx3iS0WbW9rkq0mGtCMfu9C6Ia5YEwRzomfppmjqIrP2SP8343jDSjHeWBMKrlWeHI65QNxYzntBrLtPMxMeaw+yY6QQnh5ZfZWX8aZfa9rRrOogWB6ZpTdnFfrLaK8mjSoR572EMiuibPkOhyQ6M2tjTom6GEo29G3eQJgXRp2KTG1GW/S1tbGvRJXvz5Gvo36hF8DeWZfScPJgMY5amzm52lk3p0WkOZYor0EyLylIhsd5+LUv/Z6EyR7haR1Sn3C53bvIhsyHtelpgSSd7wbRv3rnv/mEpuX9YWjW5pYooU4BpVXe4+dwG439YCbyDZqPqLIjInInPAF0hMlZ4FXJK6zyDIa1bkVcGb7tWSfVbodGEWxAiDJqZIi1gDbFLVQ6r6ODBPYsv4PGBeVR9T1eeBTe7awTFrq0RfJiqyz6p7j7YJZZNvo32amCIFuFJEdojITSnD595Mkf6aQ1W8FwxdZKBZNZUxlNwmJuHRxBTpdcBrgeXAPuCzk0tz/l7LFOkxHFfWe50wa/5Il1PXZ20EZUa3EsYY5j5ZUOaiPFOkqvp06vcbgDvd6TSTo62YIu2KadPZu9p7Nq+WkqWoKTDGEn2MYe6T2qZIJ3aNHe8BdrrjzcBaETlORJYCy4D7gQeAZSKyVESOJem43ewnGP1QpbnhowaTFQozZdEcq8H4pUyT5y3A+4ELMkPEnxGRh0RkB/A24G8AVPVhEhOljwDfBa5Q1RdV9QXgShJ7xruA29y1weJjKrfPUZkqU/7LTO83THx9E/zEttim3oc0xb3v5xtxYFPvO6LO3rJ9T3Ef0oI/Iz5MUHKo20wJYRNtExGjT0xQcmgjU/reqtIwQsQEpSNMBIwxYIJiGIY3TFAMw/CGCYphGN4wQUkRw4SvGPxojJdRCkrR/iQxdJyGPgJljJvBCEob+4J0kdFCyMwxCKnhh7bTm029r4FNaZ+OxU9c2NT7nrHMMp3Q4yeEWuFQGayg1LELbIyD0AUvZgYrKGUSjVXNwyRmkY/Z7z4YrKCkKXrJfW8t0Ob1sRNreGMqoNrYtnQUghLiS66zkjnWTFaV7E50Ywl3H/jOG6MQlCaElJhDFMYuGGu426aNeB21oJQRi5DmrBjDJ/Z0VGaT6uNF5H4R+bEzRfpJ575URO4TkT0icqvbeBq3OfWtztzofc6Wz+ReuSZK+6Kq7Vxf9zLaIfbMCPE3bcvUUA4BF6jqH5LY4LlQRFYCnyYxRboMOAhc7q6/HDioqq8DrnHXFZoo9RmYNjHBCJ8u3lEXmT3mtFbGFKmq6q/c6THuo8AFwO3O/WbgYne8xp3jfn+7M8VRZKLUMKIh5szeBaX6UJyx8+3AfmAr8FPgGWcaA440K3rY5Kj7/VngZAIzRRpztdIwfNHLsLGzq7OcxNrfecDr8y5z31GYIrWSxjD8U2mUR1WfAX4ArAROEpGJKdO0WdHDpkjd7y8HDjDdRKlhGD3Q+TwUEXmliJzkjl8CvIPE8t89wHvdZeuAO9zxZneO+/37mixpLjJRahjGQChTQ1kM3ONMjj4AbFXVO4GPAh8RkXmSPpIb3fU3Aic7948AG6DYRKnPwMTYLxKjnw2jiNHsh2ILAQ3jSCZ5wud+KAtmXzIMhigmZnbUqEtbBWzUU+/rNBdibGLk+Tm9D+5kdmWMYTP6oa0CaDRNnjFgzTqjLJO0smXvduYWz9sWkMbRmJgYZWkrrZigRII1Z4w2sP1QRorVPuJhzOI/CEEZ8wvMYnHRP2MW/0EIivFbxpyYjf4ZhKBYJjKMMBiEoBiGEQaDE5QqfQjpa63vwaiDpZsjsYlt2IQwY9yYbWPPmJgYhh9GszhwLGSr4Hm7qJuAGm0RfJPnmX0n52YAa6YUk144aBizGFWTx8SkOllTnk2wTkejCsELSh4mJv4J0aC8ER9RCooN9/qhzgZNFt/GNJqYIv2qiDwuItvdZ7lzFxG51pkc3SEi56Tutc6ZLt0jIuuKnjmLdOIfUgnadWatE3exm8o02qXMKM/EFOmvROQY4L9F5Dvut79T1dsz17+LZEf7ZcAK4DpghYgsAq4CziWxx/OgiGxW1YM+AjIEYhHHWPxpdE8TU6RFrAFucf+7l8R+z2JgNcmO+QeciGwlsXFcCisVDeNoQssXtUyRqup97qerXbPmGhGZmPkrMjlayhRpEVYqGsaRhDjaWcsUqYj8AbAR+H3gj4BFJHZ6oKEp0iLbxqEpsWH0TWhiAvVNkV6oqvtcs+YQ8BUSm8dQbHK0lCnSItvGIUZeU0wkjaFR1xTpT1y/CCIiwMXATveXzcClbrRnJfCsqu4DtgCrRGShiCwEVjm3mQw14w1RJI1x08QU6ddE5CHgIeAVwKfc9XcBjwHzwA3ABwFU9QDwT+4eDwD/6NxqEbvITPzvKxyxx4cxDIJfy2N2eQyjXUa1lsdol1k1G7NIaFTBBKVn+s6ss/pxfC40NIaPCUrPhDCVve/nG8PBBCUAuqwB5IlHyDWQsYjdUMJpgtIhISSaMuIRgj8nhCx2PhlKOE1QOiSWRNOFP0MSLcMfJiiRUzdj9p2hYxFXoxomKJGTzpixiosxHExQRkp2kyqbb2L4wASlQ9rOsEXNiLLPndUMMcExZmFT7w1j5NjUe+MorP/ECIHgBcUS/Gya7Nw17X9l1vmMgbGE0wfW5DG8E+LWhEYxo2vy+CghhlzKVAlbF/FgYjJeohAUHwk01kReRgDKzkVpS0yGLNZGNaIQlCLGkJCzQjgrzNOEs62tCGIVax/43nkvdqwPJUVe27+OuU5jXMTeZ+SzDyU6Qenr5ZmwGEPFp6CUMUUKJMa+gG3AU6r6bhFZCmwiscnzQ+D9qvq8M/h1C/Am4P+AP1PVJ9w9NgKXAy8Cf62qpXa9T5OXmSeZvc2MbiJiGLOp0ofyIWBX6vzTwDWqugw4SCIUuO+Dqvo64Bp3HSJyFrAWeAOJCdIvOpHyQtW+hjzS/7E2cTEWN0YRZU2Rng78CfAv7lyAC4CJofSbSWzzQGLb+GZ3fDvwdnf9GmCTqh5S1cdJzGxMjIM1Jru4rUqNIrZdzPrG4sYoomyT55+BvwdOdOcnA8+o6gvuPG2n+LANY1V9QUSeddcvAe5N3TPXtrGIrAfWu9ND39Pbd2avyTK3OH02XyI4ef+fP/zdEa8Aft7VwzpkqOGC4Ybt93zdaKagiMi7gf2q+qCInD9xzrlUZ/xWyraxql4PXO+evc1XZ1FoDDVsQw0XDDdsIrLN173K1FDeAvypiFwEHA+8jKTGcpKILHC1lLSd4okN4ydFZAHwcuAAJW0bG4YRLzP7UFR1o6qerqpnkHSqfl9V/wK4B3ivu2wdcIc73uzOcb9/X5Ox6c3AWhE5zo0QLQPu9xYSwzB6p/SwcQ4fBTaJyKeAHwE3OvcbgX8VkXmSmslaAFV9WERuAx4BXgCuUNUXZzzj+gb+C52hhm2o4YLhhs1buIKe2GYYRlxEvZbHMIywMEExDMMbwQqKiFwoIrtFZF5ENvTtn1mIyE0isl9EdqbcFonIVhHZ474XOncRkWtd2HaIyDmp/6xz1+8RkXV5z+oSEXmViNwjIrtE5GER+ZBzH0LYjheR+0Xkxy5sn3TuS0XkPufPW0XkWOd+nDufd7+fkbrXRue+W0RW9xOiIxGRORH5kYjc6c7bD5eqBvcB5oCfAq8BjgV+DJzVt79m+PmPgXOAnSm3zwAb3PEG4NPu+CLgOyRzc1YC9zn3RcBj7nuhO17Yc7gWA+e44xOBR4GzBhI2AV7qjo8B7nN+vg1Y69y/BPyVO/4g8CV3vBa41R2f5dLoccBSl3bnAkiTHwG+DtzpzlsPV68BnhIRbwa2pM43Ahv79lcJf5+REZTdwGJ3vBjY7Y6/DFySvQ64BPhyyv2I60L4kEwPeOfQwgacQLLIdQXJbNgF2bQIbAHe7I4XuOskmz7T1/UYntOBu0mWyNzp/Nl6uEJt8hyevu/InaYfAaeq6j4A932Kcy8KX9DhdlXhN5KU5IMIm2sWbAf2A1tJSuFSy0qA9LKS0MI2WS7zG3deerkMDcIVqqCUmqYfMY2WJ/SBiLwU+CbwYVX9xbRLc9yCDZuqvqiqy0lK9POA1+dd5r6jCFt6uUzaOedS7+EKVVCGMk3/aRFZDOC+9zv3ovAFGW4ROYZETL6mqt9yzoMI2wRVfQb4AUkfyklu2QjkLysh8GUlk+UyT5DsWXQBqeUy7ppWwhWqoDwALHO90seSdBRt7tlPdUgvQ8guT7jUjYisBJ51zYYtwCoRWehGTVY5t95wW0/cCOxS1c+lfhpC2F4pIie545cA7yDZ8yfqZSXa53KZvjvDpnQqXUQyovBT4ON9+6eEf78B7AN+TaLsl5O0Q+8G9rjvRe5aAb7gwvYQcG7qPh8g2UNhHrgsgHC9laSauwPY7j4XDSRsZ5MsG9kB7AT+wbm/xmWceeDfgOOc+/HufN79/prUvT7uwrwbeFffYUv563x+O8rTerhs6r1hGN4ItcljGEaEmKAYhuENExTDMLxhgmIYhjdMUAzD8IYJimEY3jBBMQzDG/8PgR/h5cGp2woAAAAASUVORK5CYII=\n", | |
"text/plain": [ | |
"<Figure size 432x288 with 1 Axes>" | |
] | |
}, | |
"metadata": { | |
"needs_background": "light" | |
}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"import numpy as np\n", | |
"\n", | |
"plt.imshow(mask_2016 & np.logical_not(mask_2000))" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 18, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"<matplotlib.image.AxesImage at 0x7fc07b6b4668>" | |
] | |
}, | |
"execution_count": 18, | |
"metadata": {}, | |
"output_type": "execute_result" | |
}, | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAARQAAAD8CAYAAAC2EFsiAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAFhdJREFUeJzt3X+s3XV9x/Hny/JLJ9hWgRQxoWh14qJX1tEal4WJtLVbVkw0KVlGw0hwExLNls0yk+GvJbpMWUgUxdFZjVoYamgIris/jFkyfmotLYi9AnPYhsa1oIaEWXzvj+/n4OFy7r3nx+d8f74eyc0953u+95zv59zzfZ3P9/P98VZEYGaWw0uqXgAzaw8Hipll40Axs2wcKGaWjQPFzLJxoJhZNqUHiqQNkh6RNCtpa9mvb2bTozKPQ5G0BPgRcCHwBHAfcHFEPFTaQpjZ1JTdQzkPmI2IRyPi/4AdwKaSl8HMpuS4kl/v1cD/9N1/AljTP4Oky4HLAZaw5HdfxinlLZ1N5PVvfoYf7X1Z1YthI/oFR38WEafmeK6yA0UDpr1gmysirgeuBzhFy2ONLihjuSyHB2HNoP+w1drtcfN/53qusjd5ngBe03f/TOBgyctgZlNSdqDcB6yStFLSCcBmYGfJy2BmU1LqJk9EHJN0JbALWAJsi4j9ZS6DmU1P2WMoRMRtwG1lv66ZTZ+PlDWzbBwoZpaNA8XMsnGgmFk2DhQzy8aBYmbZOFDMLBsHipll40Axs2wcKGaWjQPFzLJxoJhZNg6UKdp1cE/Vi2BWKgfKFK0/Y6bqRTArlQPFrCXq0CMu/Xoobdf/T3UPxcpUh8+beyiZ1eGfalYV91CmwKFiXTVRD0XS45IelLRH0v1p2nJJuyUdSL+XpemSdG0qQbpX0rk5GmBmL1TlWEqOTZ4/jIiZiFid7m8F7oiIVcAd6T7Au4BV6edy4LoMr21mc1TZQ57GGMomYHu6vR24qG/6l6NwN7BU0oopvL6ZVWTSQAngPyQ9kEqIApweEYcA0u/T0vRBZUhfPeHrt1Iddv+ZjWPSQdm3R8RBSacBuyX9cIF5Fy1DCi+sbXwSrpNr1iQT9VAi4mD6fRj4FnAe8GRvUyb9PpxmH6oMaURcHxGrI2L18Zw4yeKZWcnGDhRJvyXp5N5tYB2wj6K06JY02xbglnR7J3BJ2tuzFni6t2lkL+TdztZUk2zynA58S1Lveb4WEf8u6T7gJkmXAT8B3pvmvw3YCMwCzwCXTvDaZlZDYwdKRDwKvGXA9P8FLhgwPYArxn09s7rYdXCPe5Hz8KH3ZiNymMzPgWJm2ThQzCwbB0oH+cC5Zqvz/8+B0kEeA2i2Ov//HChmlo0DxcyycaCY1dhC4yW7Du6p3XiKr9hmVmMLjZfkGEvZdXAPSzJeRMQ9FOukUb/Z585ft55BXThQrPZyr7zjPN/c3kD//SaHS+49Rg4Uq7VprKzrz5jJviI1OVRycqBYI9R9hS3z2JA6Dsb2eFDWaq3OB3H1VLWMvVCp03vkHopZxSbpbdQpTMCBYhWpa5e9bG3b2+RAsdL1VoK6rQxVmHSA2D0U67y27HLtV1Y7egE0TBBV8d56UNYq0VsZmh4odV3+qpZr0R6KpG2SDkva1zdt5PrFkrak+Q9I2jLotaybxv3w12FlHra3ULaqlmmYTZ4vARvmTBupfrGk5cDVwBqK2j1X90LIbFzzrTBlBE0dwqyOFg2UiPgucGTO5FHrF68HdkfEkYg4CuzmxSFlQ2rTh3ka36RlfDPXrUdSF+MOyo5av3jousaSLpd0v6T7f8WzYy5eu/nDbHWVe1B2vvrFQ9U1hqIUKXA9wClaPnAes546Hi06yNxeZd2Xd1zj9lBGrV88VF1js7bqD5C2hgmMHyij1i/eBayTtCwNxq5L08wm0qSVs457g3IbZrfx14H/At4g6YlUs/iTwIWSDgAXpvtQ1C9+lKJ+8ReB9wNExBHg48B96edjaZp1xCgDyaPOW5eVtE6D5VWdkbzoGEpEXDzPQyPVL46IbcC2kZbOWmPYlT5XQJQdNHUKtkHKWj4fem+10PtGXehDP+gbd775y165158xU6seylxlvR8+9N4qN8yemjqvrD3TXmlH2aNVVW/JPRSr3LArSJ03KazgQLFacFi0gzd5zBqiCaHrHoqNpAljGVYdB4qNpMpvyXGu9OYALJcDxRqhf5fyKMe0WLk8hmLZ5Tw4rWec52vCmEPbOFCsNKMGjQOhebzJY9lNEgTeTClX7vdbxek39XSKlscaveiUIWuIplyrpOtuj5sfiIjVOZ7LPRSbGgdJ93gMxabKoVIf/T3G/k2dJSvyvYYDxawj5r9q3Gy21/Amj5ll40Axq1hVV1ebBgeKVaItK1AObbo0w7ilSD8i6aeS9qSfjX2PXZVKkT4iaX3f9A1p2qykrXNfx7qlLSuQvdC4pUgBromImfRzG4Ckc4DNwJvS33xO0hJJS4DPUpQqPQe4OM1rLeZeSPcMc5Hq70o6a8jn2wTsiIhngcckzVLUMgaYjYhHASTtSPM+NPISW2O4F9I9k4yhXClpb9ok6hU+dynSlnJvw4YxbqBcB7wWmAEOAZ9O07OUIo2I1RGx+nhOHHPxbBocKraYsQ5si4gne7clfRG4Nd1dqOSoS5E2mDdfbBhj9VB6dY2TdwO9PUA7gc2STpS0ElgF3EtRLXCVpJWSTqAYuN05/mJbVdp0zITlt2gPJZUiPR94laQngKuB8yXNUGy2PA68DyAi9ku6iWKw9RhwRUQ8l57nSop6xkuAbRGxP3trbOrcU7GF+PIFZh3nyxeYWS05UGzqPObSHQ4UG8okoeBxl+5woNhQHAo2DAeKmWXjQDGzbBwoVgofENcNDhQzy8YXqbZSeFC3G9xDMbNsHChmlo0DxcyycaCYWTYOFDPLxoFiZtk4UMwsGweKmWXjQKkxH6puTTNMKdLXSLpL0sOS9kv6QJq+XNJuSQfS72VpuiRdm0qO7pV0bt9zbUnzH5C0ZXrNaj6HiTXRMD2UY8BfR8QbgbXAFamM6FbgjohYBdyR7kNRbnRV+rmcooYPkpZTXOB6DUU1wav7CoR13qAA8eHq1jSLBkpEHIqI76XbvwAepqj6twnYnmbbDlyUbm8CvhyFu4GlqezGemB3RByJiKPAbgbXTO6kueExaZgM28NxT8hyGmkMJdU4fitwD3B6RByCInSA09JsE5UjdSnSPIYJpF0H97gXZFkNHSiSXg58A/hgRPx8oVkHTBu6HKlLkZbHYWK5DRUoko6nCJOvRsQ30+QnexUE0+/Dafp85UgXKlNq8xh3k2Tu3y123yyHYfbyCLgBeDgiPtP30E6gt6dmC3BL3/RL0t6etcDTaZNoF7BO0rI0GLsuTRtbF1aK9WfMjNXOxcZk5t7vwntp0zfMBZbeDvwZ8KCk3qfu74BPAjdJugz4CfDe9NhtwEZgFngGuBQgIo5I+jhFnWOAj0XEkUkWvitd9mm302FiubgUaU31D5j2VviuBKiVK2cpUl8Csqb6w8NBYk3hQ+9bzlebtzK5h9Jy7t1YmdxDaTD3PKxuHChmlo0DpcG8OWN14zGUGurflHFoWJM4UGooZ4j4BEArkzd5amAag6uDDobzIK5NmwNlyvqPA5lvhR7Ug8ix8vc/x3zP11s+h43l4EPvS1TXzY+6LpeVI+eh9+6hlKiuK21dl8uax4GSiTcZzLyXJ4scmwzeVWxt4EDJIEcAVBkiHkOxXLzJk0FTNndG2ctkNg4HSgZeIc0KDpQx5eyVlNXDcfDZtE1SivQjkn4qaU/62dj3N1elUqSPSFrfN31DmjYraeug12uKnCunV3Rri2EGZXulSL8n6WTgAUm702PXRMQ/9c+cypRuBt4EnAHcLun16eHPAhdSlNS4T9LOiHgoR0PMrHqLBkoqgdGrEPgLSb1SpPPZBOyIiGeBxyTNUtQyBpiNiEcBJO1I8zpQRuSLVltdTVKKFOBKSXslbesrfN7JUqRlng+z/owZh4nV0iSlSK8DXgvMUPRgPt2bdcCft74UqVdyswlKkUbEkxHxXET8Gvgiv9mscSlSA5pzfI7lM3Yp0l5d4+TdwL50eyewWdKJklYCq4B7KSoGrpK0UtIJFAO3O/M0w+rIPbbumaQU6cWSZig2Wx4H3gcQEfsl3UQx2HoMuCIingOQdCVFPeMlwLaI2J+xLWZWMV8PxazjfD0Ua4RRxlA83tIODpQp6vqlFUcZQ/F4Szs4UMwsGwfKFPlb17rGgTJlbQqVLm++2XAcKPYi8439tCkcbTp8CUh7EQeHjcs9FDPLxoFSIx6jsKbrfKBUtRJ7jMLaqPOBUtVK7PCwNup8oFi3eLNyuhwoNeQP/XT5/Z0e7zauWH/VPl8rdvr83k6XeygVmvtN6ctINoN7OPNzoGTiD1k56vA+O/Tn1/lAyfUBHedD5h7JaHr/qzqEig3W+TEUr9DN4f9V/Q1zkeqTJN0r6QepFOlH0/SVku6RdEDSjenC06SLU9+Yyo3ek2r59J5rYIlSM1tYU3plw2zyPAu8IyLeQlGDZ4OktcCnKEqRrgKOApel+S8DjkbE64Br0nxzS5RuAD4naUnOxpi1VVN6Z4sGShR+me4en34CeAdwc5q+Hbgo3d6U7pMevyCV4ni+RGlEPAb0lyg1m1jXL7lZB8MW+lqSSmgcBnYDPwaeiohjaZb+sqLPlxxNjz8NvJKalCId9gOXe74ceiuMVxqrq6ECJVUInKGo9nce8MZBs6XftS5Fuv6MmaFWyGG7mE3pinZB/14zB281RtrLExFPSfoOsBZYKum41AvpLyvaKzn6hKTjgFcAR6hRKdKmhkBTl7sKfq+qMcxenlMlLU23Xwq8E3gYuAt4T5ptC3BLur0z3Sc9fmcU1cTmK1FqVhvu1UxmmB7KCmB72iPzEuCmiLhV0kPADkmfAL5PUf+Y9PsrkmYpeiabYeESpXXQf06NdZc/A5NxKVKzjnMp0gn1d2vdxbVxLLTHrcufKfdQDPAmX05NuwyFeyj2IpN+Kzblw98EXT7ps7WB0uVup1lVWhsoXfuG6Fp7rZ5aEyhdPzLSh+Xn5/dydK0JFH9DF/w+WJW8l8es47yXxxqtTZsSbWpLDg6UBmr6tVXbtFk2ty2T/E/aMAbWmkBp+j9iFL0Pcf+HuUvtn6Yq38e59ZmaqDWB0qZvvVE07ajMuqv6fWz6QXGtCZSuavKHr426/v9woJhZNo0OlCZva1p7Lfa5nHu2e5s+x40OlK53L62eFvtczn28TZ/jRgfKKNr0LWDt0aYwgQ4FStv+cdZ+TdwcmqQU6ZckPSZpT/qZSdMl6dpUcnSvpHP7nmtLKl16QNKW+V7TzF54XEpTgmWYi1T3SpH+UtLxwH9K+nZ67G8i4uY587+L4or2q4A1wHXAGknLgauB1RT1eB6QtDMijuZoiFkbDVtHqi4WDZRUAmNQKdL5bAK+nP7ubklLJa0Azgd2R8QRAEm7KWocf338xTdrvyZtro9VijQi7kkP/UParLlGUq/M33wlR4cqRWrN1qRv07I1adNlXGOVIpX0O8BVwG8DvwcsBz6UZp+oFOm0axvbdDXp27RMbQ+SnpH28kTEU8B3gA0RcSgKzwL/SlHzGOYvOTpUKdJp1zY2q8LcusttNW4p0h+mcREkCbgI2Jf+ZCdwSdrbsxZ4OiIOAbuAdZKWSVoGrEvTzDqj6Sf/LWaSUqR3SjqVYlNmD/AXaf7bgI3ALPAMcClARByR9HHgvjTfx3oDtGbWDr4EpJWu60XFFrrkRBXvjS8BadZCbQha91DMKlKXAHEPxawF6hAmudU+UNq8i82sbWofKG1McRtdk79Ymrzso6p9oJiBv1iawoHSQl36RmyCMsOw6v+9A6WF/G3ebVWGyjBHyppZQ1T9ZeIeilkNNfVSBw4Uq70mrliTqrqnMS4HitVe7zKIdQ+W3MvXxFBxoFit9Z9IV8YKNkkoNO36r9PgQLFaK/tbetJQqGuvoqygc6BYJy20grX1ymplbDY6UKyThulJtClYyuo5OVDMFlHXzZhh9QfitNviQDFruV6IlBGMDhQzy2boQEnFvr4v6dZ0f6Wke1Kd4hslnZCmn5juz6bHz+p7jqvS9Eckrc/dGDOr1ig9lA8AD/fd/xRwTUSsAo4Cl6XplwFHI+J1wDVpPiSdA2wG3kRRgvRz6Ur6ZtYSw5YiPRP4I+Bf0n0B7wB6hdK3U9TmgaK28fZ0+2bggjT/JmBHRDwbEY9RlNnoFQczsxYY9mzjfwb+Fjg53X8l8FREHEv3++sUP1/DOCKOSXo6zf9q4O6+5xxY21jS5cDl6e6zt8fN++bO0xKvAn5W9UJMQVvbBe1t2xtyPdGigSLpj4HDEfGApPN7kwfMGos8NlRt44i4Hrg+vfb9ua7GXTdtbVtb2wXtbZuk+3M91zA9lLcDfyJpI3AScApFj2WppONSL6W/TnGvhvETko4DXgEcYcjaxmbWXIuOoUTEVRFxZkScRTGoemdE/ClwF/CeNNsW4JZ0e2e6T3r8ziiK/+wENqe9QCuBVcC92VpiZpWb5IptHwJ2SPoE8H3ghjT9BuArkmYpeiabASJiv6SbgIeAY8AVEfHcIq9x/QTLV3dtbVtb2wXtbVu2dtW6cqCZNYuPlDWzbBwoZpZNbQNF0oZ0iP6spK1VL89iJG2TdFjSvr5pyyXtTqcn7Ja0LE2XpGtT2/ZKOrfvb7ak+Q9I2jLotcok6TWS7pL0sKT9kj6QprehbSdJulfSD1LbPpqmt+K0kkpOl4mI2v0AS4AfA2cDJwA/AM6perkWWeY/AM4F9vVN+0dga7q9FfhUur0R+DbFsTlrgXvS9OXAo+n3snR7WcXtWgGcm26fDPwIOKclbRPw8nT7eOCetMw3AZvT9M8Df5luvx/4fLq9Gbgx3T4nfUZPBFamz+6SGnwm/wr4GnBruj/1dlXa4AXeiLcBu/ruXwVcVfVyDbHcZ80JlEeAFen2CuCRdPsLwMVz5wMuBr7QN/0F89Xhh+LwgAvb1jbgZcD3gDUUR8MeN/ezCOwC3pZuH5fm09zPZ/98FbbnTOAOilNkbk3LOfV21XWT5/nD95OBh+k3wOkRcQgg/T4tTZ+vfbVud+oKv5Xim7wVbUubBXuAw8Buim/hoU4rAfpPK6lb23qny/w63R/6dBkmaFddA2Wow/QbbKLTE6og6eXAN4APRsTPF5p1wLTati0inouIGYpv9POANw6aLf1uRNv6T5fpnzxg1uztqmugtOUw/SclrQBIvw+n6fO1r5btlnQ8RZh8NSK+mSa3om09EfEU8B2KMZSl6bQRGHxaCTU/raR3uszjwA6KzZ7nT5dJ80ylXXUNlPuAVWlU+gSKgaKdFS/TOPpPQ5h7esIlaY/IWuDptNmwC1gnaVnaa7IuTatMuvTEDcDDEfGZvofa0LZTJS1Nt18KvJPimj+NPq0kqjxdpurBsAUGlTZS7FH4MfDhqpdniOX9OnAI+BVFsl9GsR16B3Ag/V6e5hXw2dS2B4HVfc/z5xTXipkFLq1Bu36fopu7F9iTfja2pG1vpjhtZC+wD/j7NP3stOLMAv8GnJimn5Tuz6bHz+57rg+nNj8CvKvqtvUt1/n8Zi/P1NvlQ+/NLJu6bvKYWQM5UMwsGweKmWXjQDGzbBwoZpaNA8XMsnGgmFk2/w//1fDjmwcOSAAAAABJRU5ErkJggg==\n", | |
"text/plain": [ | |
"<Figure size 432x288 with 1 Axes>" | |
] | |
}, | |
"metadata": { | |
"needs_background": "light" | |
}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"plt.imshow(mask_2016 & mask_2000)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 52, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"(4000, 4000, 4)" | |
] | |
}, | |
"execution_count": 52, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"\n", | |
"\n", | |
"b = mask_2016 & mask_2000\n", | |
"r = mask_2016 & np.logical_not(mask_2000)\n", | |
"g = np.logical_not(mask_2016) & mask_2000\n", | |
"a = np.logical_not(np.logical_not(mask_2016) & np.logical_not(mask_2000))\n", | |
"\n", | |
"rgba = np.dstack((r,g,b,a))\n", | |
"rgba.shape" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 50, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"<matplotlib.image.AxesImage at 0x7fc07b3ec860>" | |
] | |
}, | |
"execution_count": 50, | |
"metadata": {}, | |
"output_type": "execute_result" | |
}, | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAARQAAAD8CAYAAAC2EFsiAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztnXvwZ2V931+fLrcQIezGSxalA+qakTh2VcbFidPSGliyTYvOJHHpxfUyQwzaam2nLjUt1OAkbS7VTC0GCwkk6kqJjjtUu12RjM1MUC6uXCTAT6QJLoWxi0THDgn46R/n+cLh7HM75zzn8v1+P6+Z3/y+3/M957mcc5738/l8nuecR1QVwzCMEvyNqQtgGMbqYIJiGEYxTFAMwyiGCYphGMUwQTEMoxgmKIZhFGN0QRGR80XkXhHZEJG9Y+dvGMZwyJjzUERkE3AfcC7wEHALcKGqfmO0QhiGMRhjWyivBTZU9QFV/StgH3DByGUwDGMgjhk5vxcCf1H7/hCwo76DiFwEXASwiU2vOZGTxyudYawh3+Ox76jq80qkNbagiGfbs3wuVb0SuBLgZNmiO+QNY5TLMNaWL+r1/7tUWmO7PA8Bp9W+vwg4PHIZDMMYiLEF5RZgm4icISLHAbuB/SOXwTCMgRjV5VHVJ0Xk3cABYBNwtarePWYZDMMYjtHnoajq51X1Zar6ElX90Nj5t+HA4UNTF8EwlgqbKRth56nbpy6CYSwVJiiGYRTDBMUwjGKYoBiGUQwTFMMwimGCYhhGMUxQDMMohgmKYRjFMEExDKMYJijGymEznKfDBMVYOWyG83SYoBiGUQwTFMNYY0q7hyYohrHGlHYPTVAMwyiGCYphGMUwQSmMDVka64wJimEYxeglKCLyoIjcKSKHRORWt22LiBwUkfvd/81uu4jI77glSO8QkVeXqMDcsDkQxjpTwkL5u6q6XVXPct/3Ajeq6jbgRvcd4GeBbe7vIuCKAnkbhjEjhnB5LgCucZ+vAd5Y236tVtwMnCIiWwfI3zBmw9QxtbHz7ysoCvxPEbnNLSEK8AJVfRjA/X++2+5bhvSFPfM3jFkztQs8dv591+X5aVU9LCLPBw6KyJ9F9k0uQwrPXtv4BE7sWTzDMMakl4Wiqofd/0eBzwKvBR5ZuDLu/6Nu96xlSFX1SlU9S1XPOpbj+xTPMCani8tRyk2ppzOW69NZUETkR0XkpMVn4DzgLqqlRfe43fYAn3Of9wNvcaM9ZwOPL1wjw1hVurgcqWO6iMNYrk8fl+cFwGdFZJHOJ1X1f4jILcB1IvIO4M+BX3D7fx7YBWwAPwDe1iNvw5gFBw4fGj1OkZvfFPEbUT0qjDEbTpYtukPeMHUxDGM0mgJVUrBCaX1Rr7+tNu2jFzZT1jAGoGvMotngm+LShzEsFhMUwxiALo03JRhTD0HnYIJiGEvKgcOHJp8416TvPBTDMArR1gKZo8ViFoqxtozZu5fIa27WiA8b5TGMNcdGeQzDmCUmKIZhFMMExVg7liEWsayYoBiGUQwTlMJY7zd/+gy32vWNY4JSmDnODTDKYdc3jgmKYTiGtj7WwboxQTEMx9DWxzpYNyYohmEUwwTFMIximKAYa8siprEOsY2xMEEx1pZFTGMdYhtjYYJiGEYxkoIiIleLyKMicldtW+v1i0Vkj9v/fhHZ48vLMIzlJsdC+X3g/Ma2VusXi8gW4FJgB9XaPZcuRMgwUowR47A4ShmSgqKqXwaONDa3Xb94J3BQVY+o6mPAQY4WKWMNyWnIY8Q4LI5Shq4xlLbrF2evaywiF4nIrSJy61/zRMfiGcuCNeTVonRQNrR+cda6xmBLkRrLg7lJR9NVUNquX5y1rrFhLNPcELOujqaroLRdv/gAcJ6IbHbB2PPcNsN4mrqIWGNdTnKGjT8F/CnwkyLykFuz+NeBc0XkfuBc9x2q9YsfoFq/+OPAxQCqegT4VeAW9/dBt80wnsWqCckyWFolsbfeG0vJFIuUryr21ntj7VlVMVl2i8YExTBmxLILpQmKYRjFMEEx1pJldy3migmKsZYsu2uRYirBXCpBOXD4kPUshpHBVIJ5zCS5dmTVexXDWHaWykIxDGPemKAYhlEMExTDMIphgmKMigXVVxsTFGNUhgism0jNBxMUY+mx0b/5YIJiGEYxlkpQ2pi2ZgYvF3a9VoOlEpQ2pq2ZwcvDXF/7OLfyLANLJSjGatBsqHMV/7mWa86YoBijYw11dem6FOllIvJtETnk/nbVfrvELUV6r4jsrG0/323bEJG9zXyM9WbnqdtNaFaArkuRAvwnVd3u/j4PICJnAruBn3LH/BcR2SQim4CPUi1VeiZwodvXMIwVIvm0sap+WUROz0zvAmCfqj4BfEtENqjWMgbYUNUHAERkn9v3G61LbBjGbOkTQ3m3iNzhXKLFwue2FKmxNtgo0NF0FZQrgJcA24GHgd9y220pUmNtsJjP0XQSFFV9RFWfUtUfUi3otXBrbClSw5iAuVhLnQRlsa6x403AYgRoP7BbRI4XkTOAbcBXqVYL3CYiZ4jIcVSB2/3di20YRp25WEvJoKxbivQc4Lki8hBwKXCOiGynclseBH4JQFXvFpHrqIKtTwLvUtWnXDrvplrPeBNwtareXbw2hmFMii1FahgFWOalUW0pUsOYCYvYxbKKSWlMUAyjByYkz8YExTCMYpigGJORGuqcy1CokY8JijEZKXfB3InlwwTFMIximKAYQZoux1AuiLk2q4MJihFs0E2XYygXxFyb1cEExbAGbRTDBGXNaetumHtixDBBWXMW1kmuUKyyNXPg8KHZC+bcy2eCssZMeXPOsWEsw3tt514+E5Q1pfkw2xQ36hxFZU6E1isqed5KXwN72tjwssxPz86NUg8QDrWekT1tbAyOicnRdOnNF8Jcj1XlpLMM8RwfJiiGgd+9aDbqLiLrO2bnqduz3JjUPKA5Co65PEaUIVyfublTY5Qn1vh9eafKVPI9LCVdHhMUI8jcGv6QLOo65rrLsfMbs4xyrJk2jCooInIacC3wE8APgStV9SMisgX4NHA61Xtlf1FVHxMRAT4C7AJ+ALxVVW93ae0BfsUlfbmqXhPL2wTFgPGEzZfPkHn3FYZS4jd2UPZJ4F+q6suBs4F3uWVE9wI3quo24Eb3HarlRre5v4uo1vDBCdClwA6qZTcurS0QZhhBxrKSls0am+O8maSgqOrDCwtDVb8H3EO16t8FwMLCuAZ4o/t8AXCtVtwMnOKW3dgJHFTVI6r6GHAQ/5rJhmOOQbd1Y8wGW8+rzbWfk7C0GuVxaxy/CvgK8AJVfRgq0QGe73brtRypLUX6DHO5SVaFOQm0ryx9R5TmQLagiMhzgD8C3quqfxnb1bMtezlSW4p0nsypMXZlzFjImOnO6dpkCYqIHEslJp9Q1c+4zY8sVhB0/x9122050hVkWXtMH6UbYMlzkzvNfk4iUicpKG7U5irgHlX97dpP+4E97vMe4HO17W+RirOBx51LdAA4T0Q2u2DseW6bEWGuN84yM/UEsZw5KakyxYabp7xncoaNXw/8L+BOqmFjgH9DFUe5DvibwJ8Dv6CqR5wA/WeqgOsPgLep6q0urbe7YwE+pKq/F8vbho2NVaXU6yLqE9y6DnHbxDbDyMA3m3Ruk/VyhGXo8pYUlORi6YaxjISEw8RkWOzhQGMShvTz21ghU74CM6eMoXeizBVzeZaIVZinMDRzc2n6MFZdzOVZc1alwZSiHiupj5Is43la9k7DXJ4lYk5TrOdE7gN9y+A+LMq9rNfZBMWYhOaLjPqS837cZWmsOa80mCsmKMZoDGXOL0NDK8HchRBMUIwByF3aNLU9l9jM13URm7lggjIT5rw0whyJPfMyxfIg63DOczBBmQklb/ypTeNS+ccejJvbC5ynzn8umKAYk1N/oK3+v4+LFBrlGauhTy3quZQ+HzYPxZiE+kuhm3NHhppDsiyNfJkxC2XNmerR/dgb5mNzSJaBZShryWU46pigrAExU3/sXrtpjbR5PL8v6+7upF47WQJ7lscoSmqtmeYLhNpM4pprQ11mDhw+xKatG7a2sVExN/Pa9+Z2n3ik3ky2sF7qf6tA6evVJdCcE/juiglKD+bQmOfc0FJT3RdCMfQzNl3SHaosvthR3/Ta3gOly1DHBKUHQ16YZaFN/X2T0eq95FDiWGqR85A10PYemEMnMFQZcl5SfZqI3CQi94jI3SLyHrf9MhH5togccn+7asdcIiIbInKviOysbT/fbdsQkb2+/JaNPktHLgttF/oOpdF1Iaspac7AHeMtcMscPM55SfVWYKuq3i4iJwG3Ua0S+IvA91X1Nxv7nwl8imq50VOBLwIvcz/fB5xLtaTGLcCFqvqNUN4WlJ2WIRcOT4nU2O8zyV3XeFnfsxJj1BcsuSUwFisEfk9EFkuRhrgA2KeqTwDfEpENKnEB2FDVBwBEZJ/bNygoxnT4xKTt8g+xhhfbf4pGO/f3zy4LfZYiBXi3iNwhIlfXFj63pUiXnC5xgthoQ67LVI+pzNUlWtYXN41Fn6VIrwBeAmynsmB+a7Gr53BbinTm+CyFvj107pBm6MG6xfFdn8EZ+kXYC0Lna8xnh+ZC56VIVfURVX1KVX8IfJxn3BpbinSJaN70zQf1YsOSvnklobkjvnzq/0P0CeQOEfNpM3N3zvNnhhK6nKCsANcAR1T1vbXtW118BRH5F8AOVd0tIj8FfJJngrI3AtuoLJT7gDcA36YKyv4jVb07lPc6BmVLxA/aplF6XkRO+rFneXKtmqZ71BSfOTXmuZWnzqgrB0aWIr2Qyt1R4EHgl2oC8wHg7cCTVC7SF9z2XcCHgU3A1ar6oVje6ygoMYa6KXNGXErSZqg5Vo6mmOQGhOcwgjQnbCnSNWDMm7CUYLR5yVBO3VJiEbN85tyA50ZJQbGZsjNljAYxRNCwjZWQm79vUlwolhGKzTQDvUMxx0DsmOUxQenJ3G6eXPoGOHMCtX3OTUgQmvk0xcYnPKHZrm1FLYc5BmLHnKFsb2zrQD0YOLebB9LuUsiVSKW5OGZB8/iSApJLc8g2lE6uOxayfHJfuzB3hi63WSgdGPphtr6EXiEQGrqdmi5zVZrb6pPiQnn4rJVmGVJzSdqMRK0jFpRtyRBWyZCWTijtLg1iyBmsKQsgZwZtTrA2Z79YnnPqRErdNxaUnZAhJkt1GUYtQd3fb05I6yJCXc9NV8HzuTupGazN2EssDd/xc7FM5lKOJiYoLUldyDYxibYNsEuDXfSwsZhB1ynjpQOQOdZHyIJIPWMTG/3xWT6hIHAzvea2UD4laYph6nqlrn9JltblmUNAdMwAXU6gNVWOUOPp2gPnuClt0wodN9QEu1CAOiRQuedqyNnOMfcvp9xNbGKbYw6iUprmDe6bYl4nddPF8qgf33cYOVT+Zh4+CyOnjCXEpG06qTr4rlWo/CWJdQwhizN0DcAEZRaUELMhA7yxm61v40zFPNrWKdQIQun1HWL2iUJo/xxraYhOrc81yrVUF1hQdgaUuIlK3Yg5/ntJ376NpZRD7NjcdH2BV9/vzcYUiwOl4g5d65wTh2vGjFIB8wVTW+wmKB2ZU5Q9dRONdZNNmU/9eoQENWYRdjmHXUQ6ZY0188opc11ocgLqQ14nE5SOTNETdBGxocoZunG7xmR86XTpiYfqvX0NOze9UDwoZenEzsFch7QthtKCOQaBc4d3h/DJc+IcJeIpqbzbEBvZaXtciXuh9D3lO+9jjvKYhdKCnGDXEDR7tGbMJFW2PjGPNr1wzNVo03CHJOUaxY7rYpGktpeMo+WUYWhMUAoyZGNojtCkXILUPjl+dsqcDqUfcg9Cwlg/Licu4RPS0ue+Xrd6fCInH995S9WrqwDEzlnu+SyJuTwTk3OTpkz0UnGLelopYsPTof1yxS3X5ekqlrnEArF9yzgEOZ2E71zbPJSZMfSNEwrqNekrLKkGmprL0vzdV56Q4MTy61pW3/YuhBpis+y+76FjQvVuE/sI5ROrh++YTVs3Rn2n7AnAl4Hjqd6fcr2qXioiZwD7gC3A7cA/VdW/EpHjgWuB1wD/F3izqj7o0roEeAfwFPDPVfVALO9lEZQhiDVYH13EJMeMj6WXe3yO0NTTC+WVQ0xQYhZVTChjpILRITHu0gHFBCdlucYEeGxBEeBHVfX7bjmNPwHeA7wP+Iyq7hORjwFfV9UrRORi4JWq+k4R2Q28SVXfHFqiVFWfCuW9boLSZ8RmCEEJWQo5Izq+mz00whIbQckVrND3VDlj5JznqUb9UuLURrQmc3lE5EQqQfll4L8DP6GqT4rI64DLVHWniBxwn/9URI4B/g/wPGAvgKr+mkvr6f1C+c1NUHyB0RJpNtMK9XJ1Ur15yhLIKVNs35RrEXIJfOn6rIPQ7zkuU4rc89HW+ipF23R990ub8owuKCKyiWqR9JcCHwV+A7hZVV/qfj8N+IKqvkJE7gLOV9WH3G/fBHYAl7lj/tBtv8odc30jr4uAiwBO4MTXvF52lahnZ3KDiqVuqJzG7Ns3RE4sIVec+jTePudniIYbE74u6SzS8m3vQk5wu829EmP0eShuhcDtVKv9vRZ4uW83938llyKNXbA+wVBfb7n461uunDKF4gr1MrRxEZp/qXIstjd/9zWWvqJU/5yqX5tAZ1NI6mKVKnOzHCmrsmkhx/KfglbzUFT1u8AfA2cDpziXBp69rOjTS466338MOMISLUVaSvmHIFS2lGkeutGa6XURqJQVUm8kOQ0l1SBiv+cEdXOttjaERCm3Lr7yhY5v07mNLS5JQRGR54nIKe7zjwA/A9wD3AT8vNttD/A593m/+477/Uta+VX7gd0icrwbIdoGfLVURUrSVURSjdq3bQzByu0pF2VKiYZPFFI9s084QmWIWXz1hhKzgmIWTyjfEtSthhwrr37umnXqU64hYjs55IzyvJJqbeNNVAJ0nap+UERezDPDxl8D/omqPuGGmf8AeBWVZbJbVR9waXmXKA0x16Bs1/19/m/fCx7KIxbs7UJOrx7av6/fHwvu+vId4jyXpFm+ZpxqQeg65qRZTyOFTWxbAYbuPVJB5BAxC2LxeyqN3H1LxURCaeaO5oTSHuL6+EZkmnn5BMf3my/tNiK0wATFyCbV8Lr+XoIuDTbUIH1phwQlJeZthaiLAMViIz4rJTTqE0s3t74mKC2Ywo8ckr5DnLFRhDaNcwhy6tW2HDF3ok0+Q91DKTFq676EXKixBGUlnzYeapRmyEYVGxFIBepSozeL47v2yrn1zgn++tJvlrOZZ9cRC9/IUdcy+oLAJWhTnjZuZu4IYGlWUlCGOoFD91Il9232TqFRlxJlCuWfk08o76ZF0aU8OQIdOiY2KtO1PDllC+XXNQ2fqA7JSro8y+bm9A0ChnzsRZq57k0J16Zv2r5RmlD8wCc8XcrY9rx3PTbnvpzi3rUYygrRZVgw97eSPVJfscnt2WOCtGCIWE+XIeyx4iyh/Lvu37znTFBWgNBwYXPbYntOjx3LJ0TTNRrKLC6dzxCjUL7h2dwh9lAgtEuQNEUzzb6UFJRj0rsYfenqL/tiCKVuyma+sYZRgnqaJdJvGz9ok2auWDWvT/Pc1WMXQ1y3Ehw4fIhNW8ulZxbKgLSJyvdJx9eztj12DMskld9Yw9NtyY25tImtlIyV9E3LXJ6J6OK/Qtzc75peqFwl4hzN/PoQc+m6BFNjx5WI8+SkXUosfFbhFJjLMwFdL3iJeQFtGknJRte1gaaGXGPli5Gqk0+wcoPVvuvkE5Cuwducc9KWOY5mruQ8lFL0GVpsm36M5jySXFeqzcjP0DemLx60yLdv3vWG1RTUhdDUv7fJN0cUfMfkiJAvppRzbWOW0tQu41q7PEMqfO5FzzWNc1yllAkdGuaMDX92GeZNlS1FCasoFPjNLVto25D3Sxfxqh+/oO2xFkNZYnKEpsuIQJubsG2wuE8Djw2bpkY/uo6O+Y6NnZtY2cakjyiYoGSQIyhdAqVtG+kUw31d8kzVLdYT98F3jmLppkZ6upYrVa8cC6B04DWXXAulyzVOYYJSkJBoDBU/CaUbG6npMhIUOi7UGzcpYaF09fX7jtbE0g3RtKRyCbmWXe6Zvvdc147PnjYemFiQry+LgF3zooeCeF1jDqlAYPPYerl8Qlff7hsR8eUdOneh8g5pAabSTtWpSVc3Lcdyy83bd/zYblqTtbVQQj10TtBuDsN1oXiLz7rJ6UXbWCzNY3z75lpZvnI0haiNxRLbd+prlsvQLnbzHM9lKdLfB/4O8Ljb9a2qesitNPgRYBfwA7f9dpfWHuBX3P6Xq+o1sbxjglLqpMcCcTkjHFOT26hz/O4UbXtjnzCEjm07+tKFUtcrNoLWJ4jextLM7dRy9hs1hhJZivSdwA2ehbp2Af+MSlB2AB9R1R0isgW4FTiLaj2e24DXqOpjobxLWihdLtZcBCUnENf199y6xmiO2IT+L/at593c1ixb6Jh63qHfQmWN0SVo34VQELbNfdqH+nUpaaEkYyha8X339Vj3F1OhC4Br3XE3U63fsxXYCRxU1SNORA4C5/crfh7NmEiOP506ps9Fbtt4Yz1Zrs/u26/+W9sYQiqPuoj4XJeuoxWh+FKbHrtZ1lD6KUJWVc41aRMny0mzLUMFvjstRaqq73cuz+uAJ4Abgb1uGY0bgF9X1T9xx94IvB84BzhBVS932/8t8P9U9TdD+foslC69QpcI/lCjPn3jL6HYSZtjhyRWrpwRruZ+zbSb6YW2pVzBNuevr4XYhVA9hrBWRn+WR1WfAra7Bb8+KyKvAC6hWgj9OOBKKtH4ID2XIm2sbXzUAW1OaNueIJRXF0FKpdmG5khM3/Ri5DTGVN6hnttntaTybJtPCp8rlmq0sfp0LUeKetpNi7RrhzoGXZciPV9VH3ZuzRPA71GteQzhJUezliItvbZxXz8Xju7VhurpF2mnbvLF720skz7noSmsTWGLnZOmFZAbfA3l0TWu1WygMREL5ROzsHJdzzaUSnMsMYHuS5H+mYuLLIK2bwTucofsB94iFWcDj6vqw8AB4DwR2Swim4Hz3LbBiN28Ocf6eodQuiWoN9ZUEDXHTG+mu/hc355TntD+McspN4/YfrkB46YQp1ym1P9mur4yhs6v77hQ/XJpXr9QB5GKD41BjoWyFbhJRO4AbqEKrN4AfEJE7gTuBJ4LXO72/zzwALABfBy4GEBVjwC/6tK4Bfig2zYobVyE+gWp3xRjmow+usRI2gTdfPGM1HkrYfHk/BbLp03D8V3Ten5tBTd0fkveK23d2yFdsFzWdmJbDF/QbmpRSeELbKZ6/VCjGKKevnxDMZk2QpjjAuXEZ3KC8LFzNbT12qRkfjb1fiBCVkmOmExpZkL/3izlavWhq6Xnq1PI7UjlnZN+LI2ceEsbq7ALQ6RdOs21EZScExfySXPIuTmHDOh2ZWjLKxTI9Vl9zXMT+p4bswmJbD0+1rwm9fhEl/MSSjdE23us9P1T+tqby5Mgx7ScuzsE8TkbQ+XTdiQKwvNTQkHOXJcn5V7FAq6xIG+fOFUfSl67tX19Qd+T2Pb4kM+cE8Ef26du5h2LAY1Rtj7XKhWIbdZtsT2HtkHfmMD5fmsbe+lCaOSpKxZD6UBqSC120evmb2y4rnlsH/Hq6h6Fgq1NV2BooSuVftMdCbk0fUecfCN8bTqKWCwq9FtX96VtHMm331Cu91IJSt+bNNf3Dm0P9fB9Yi++fEoMyc7dBUtRb9ShGEdMdH3XuotApLb1oW28ZUEq9hTKK/a9FEslKH0odQJjZm6sZ+uaRz2vZsPykeOOzZGQGV8/x77tPvenTiro29dlCtEmvTZlaJ6P5jnoWo5SLFUMpQRjx2FKpNfGZ16GAHGdWN1y4hE+qzEkQs396tty0h/6vLaNbbWJ8cRYu6DsmI0kFPQbi1iwsW1MZpmEBfJNed8+uecrNdJTOuDZl9R9GBOV0H7NfdZOUMZmisaYa6nM4SYfk5wh3FQQOse1aVo2zTRKkbqGXQQyJrChPOvb1m6UZ2w/cI5iAvPoMcema5wgFOj0jcIs9h/TCg6RG/hf7BOzXHzfh7a8zUKZIetoifhIBWqb+8QIBXab6Q8V1A7lleuW+X4PxYliZfCdx7WzUGKMbb30yb/LEGHfPOdI2xGQUkOeoRhJn8Br7siM73sqr+bQuO+4NuUdY+h46QVl6p68S/45N9KqExLW5tB4bJ/60Gluns1GWs8np8dvHpuyMBb7+FyslFtWTz82+pQr0E1hGoKsV0DOmTm6B7nDlV2ZW33bEmokTaGo7+s7PhZXyS1DbKSnrVvVTNsnAn0ac+xcLPKL3Rsh8du0tXORjsJiKAPSdwhyjmLZldK9os998cUW2lgescBnSBy6DOX7YkC5QeJUfbrcczZsvET0FYXcG2gM4fE1Uh85cyZiDbW+/2JbKO/QuckZHo7V01ePrgHbVJlzrmHMguqLCcqM8d08pdL0ba/TpbfsUxbfTR4jRwhy8s9JN6dcIXGLWTuhNNrWITTa5Eu/jbXV5ZyaoPSg602cOsZnag7psvh6+1i5xrJgQoREqM15baYXskxy4yAxl8L3OVW3NrSxkvp0HDlMMmwsIptE5GtuIS9E5AwR+YqI3C8inxaR49z24933Dff76bU0LnHb7xWRnSUq0JYxG9aQedXTzu2Jp8A3whELuIaOr//VSY1yxEZS6vuEjqmPAMXyaEPKZVyUL7RPqr5Tkm2hiMj7qNYlPllVf05ErgM+o6r7RORjwNdV9QoRuRh4paq+U0R2A29S1TeLyJnAp6jW7zkV+CLwMreImJez/tYJuvnOf9CvhiOQ8u8Xn1eR1A1c0jSPpZOK5+SISp9RmBxLKyRisYDsGPfP6BaKiLwI+PvAf3XfBfh7wGKh9Guo1uaBam3ja9zn64E3uP0vAPap6hOq+i2qZTYWi4N5ue+Oo1cO7MNQ6t3sceu9i69XzSnX1D1NLqH65VgSoW0+UseGrJhmgwxZR33EJFSmZvrN76H7oml95gac50Du2sbXA78GnAT8K+CtwM2q+lL3+2mBFsBUAAAEQUlEQVTAF1T1FSJyF9XKgg+5374J7AAuc8f8odt+lTvm+kZeTy9FCryCZxYQWzWeC3xn6kIMwKrWC1a3bj+pqieVSCg5sU1Efg54VFVvE5FzFps9u2rit6y1jVX1Sqq1khGRW0uZYnNjVeu2qvWC1a2biNxaKq2cmbI/DfxDEdkFnACcDHwYOEVEjlHVJ3n2OsWLNYwfEpFjgB8DjpC5trFhGMtLMoaiqpeo6otU9XRgN/AlVf3HwE3Az7vd9gCfc5/3u++437+klV+1H9jtRoHOALYBXy1WE8MwJqfPszzvB/aJyOXA14Cr3PargD8QkQ0qy2Q3gKre7UaGvgE8CbwrNsLjuLJH+ebOqtZtVesFq1u3YvWa9cQ2wzCWi6V/fYFhGPPBBMUwjGLMVlBE5Hw3RX9DRPZOXZ4UInK1iDzq5uEstm0RkYPu8YSDIrLZbRcR+R1XtztE5NW1Y/a4/e8XkT2+vMZERE4TkZtE5B4RuVtE3uO2r0LdThCRr4rI113d/r3bvpSPlTSZ5HEZVZ3dH7AJ+CbwYuA44OvAmVOXK1Hmvw28Grirtu0/Anvd573Af3CfdwFfoJqbczbwFbd9C/CA+7/Zfd48cb22Aq92n08C7gPOXJG6CfAc9/lY4CuuzNcBu932jwG/7D5fDHzMfd4NfNp9PtPdo8cDZ7h7d9MM7sn3AZ8EbnDfB6/XpBWOnIjXAQdq3y8BLpm6XBnlPr0hKPcCW93nrcC97vPvAhc29wMuBH63tv1Z+83hj2p6wLmrVjfgROB2qlnd3wGOad6LwAHgde7zMW4/ad6f9f0mrM+LgBupHpG5wZVz8HrN1eV5IfAXte8PuW3LxgtU9WEA9//5bnuofrOutzOFX0XVk69E3ZxbcAh4FDhI1Qt/V6sJm/Dscj5dB/f748CPM8+6fRj418AP3fcfZ4R6zVVQsqbpLzG9Hk+YAhF5DvBHwHtV9S9ju3q2zbZuqvqUqm6n6tFfC7zct5v7vxR1qz8uU9/s2bV4veYqKKsyTf8REdkK4P4/6raH6jfLeovIsVRi8glV/YzbvBJ1W6Cq3wX+mCqGcop7bAT8j5Uw88dKFo/LPAjso3J7nn5cxu0zSL3mKii3ANtcVPo4qkDR/onL1IX6YwjNxxPe4kZEzgYed27DAeA8EdnsRk3Oc9smQ0SEavbzPar627WfVqFuzxORU9znHwF+BriHJX+sRKd8XGbqYFgkqLSLakThm8AHpi5PRnk/BTwM/DWVsr+Dyg+9Ebjf/d/i9hXgo65udwJn1dJ5O9W7YjaAt82gXq+nMnPvAA65v10rUrdXUj02cgfVazL+ndv+YtdwNoD/Bhzvtp/gvm+4319cS+sDrs73Aj87dd1q5TqHZ0Z5Bq+XTb03DKMYc3V5DMNYQkxQDMMohgmKYRjFMEExDKMYJiiGYRTDBMUwjGKYoBiGUYz/D3kes4bhoSodAAAAAElFTkSuQmCC\n", | |
"text/plain": [ | |
"<Figure size 432x288 with 1 Axes>" | |
] | |
}, | |
"metadata": { | |
"needs_background": "light" | |
}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"plt.imshow(b)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 53, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"<matplotlib.image.AxesImage at 0x7fc07b0b9400>" | |
] | |
}, | |
"execution_count": 53, | |
"metadata": {}, | |
"output_type": "execute_result" | |
}, | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAARQAAAD8CAYAAAC2EFsiAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztfX3sLkd13nNyjU0KNL5ODHL5ECa5LZiqcXxvbUdUFYVgG7eqiZRIF1Xlilpy2hqJqFUbu5FKAqkUqiaukIDEKW5MlMS4JBFXllP31hBV/QPb94Ix/ojjy0dTgoWdXuMkQnJqevrHzu47O3tm5szs7O7s+9vn5+v3fWd3Z87MzjxzzpkvYmZs2LBhQwl8z9ICbNiwYX+wEcqGDRuKYSOUDRs2FMNGKBs2bCiGjVA2bNhQDBuhbNiwoRhmJxQiuo6IniSis0R0y9zpb9iwYTrQnPNQiOgQgD8C8A4A3wDwEIB3M/PjswmxYcOGyTC3hnIlgLPM/FVm/ksAdwG4YWYZNmzYMBHOmzm9VwP439bvbwC4yr6BiG4CcBMAvOxlLzv6xje+cT7p5sKZM8DRo0tLURx7mq29x5kzZ/6UmS8uEdfchEJCWM/mYubbAdwOAMeOHePTp0/PIdcyIAJKmJyl4tlwIEFE/6tUXHObPN8A8Frr92sAfHNmGepBKRJYC5mQ1J9s2CfMTSgPAThCRJcS0fkAjgM4ObMM82FrQH2shfjmwh7Wj1lNHmZ+kYjeB+A+AIcA3MHMj80pw6zYGtCGEPawfsztQwEz3wvg3rnTzQGBwJjgpW8+jw17im2mbADFyaRVcTcy2bCn2AhlTrREsoe286SosbxqlKkCbISyBKbSUPa1kteo0dUokw2iRerDRij7gs0vs6HFgqb1Rij7goNKJvuqlY0B82L1YSOUDevG1A2nVsKS5FrIzLGxEUqLpV5ErRV2Q4NaNT9JLkkzmbl+bYTSYqmKU2uF3TeUaFhrJP+Z69dGKBsOBko0rCXJfyVkthHKhr0DiYvaN8yBjVA27ODrBVfSO7aYZLnE0liJabwRyoYdfJV2JZV5r7ESUq+bUM6cWVqCPJjhu5XUgQ21IFRhJiL10uZh3YSytv0EnXkAW8d+gLDS3qO0eVg3oeRgm0+yYQmM7T3G1p9UjXii+jr7fiiTg3nRdS1kMf6moGxQY2x95URdY6L2sX8aCjA/mdhsz+S9tGFDEHtQWfaTUOaCvarTIjH7p5fb9qDyHGho398Be8+jCIWIvk5EXyaih4notAm7iIhOEdFT5vOwCSci+og5gvQRIroimkDtozwuWzAPnVy+CtWaZmtBjqxryl8qtFqwc1+wSGJxKstTddtE76aEhvL3mPlyZj5mft8C4H5mPgLgfvMbAN4J4Ij5dxOAj0djXtsoj4R9GerJyce+5F1AtD16bhhVJMqHVbetyIdyA4A7zfc7AbzLCv8kN/g8gAuJ6JIJ0heRTchT9rL2S93n3nyfYN5TtD2aG5LneRSuB3MvQxhLKAzgvxHRGXOEKAC8ipmfBgDz+UoTLh1D+upg7AVNnlGEPOYlh57tOXN5GFxqf4uNrMYjswyT53lIFXXE+/OmX+mw8VuY+ZtE9EoAp4joDwP3Ro8hBfpnG7/uda8bKV4BjPV1uM/bFSbCcs0hHgVe/B6bHrMh6mWfIe3a48RIDYWZv2k+nwHwewCuBPCt1pQxn8+Y21XHkDLz7cx8jJmPXXyxdX7zkrtRFZgjENyWz7P7FnOBtDcsCrXJYdWB7plYfY+0CTvtuUyfbEIhopcR0Sva7wCuAfAomqNFT5jbTgD4jPl+EsB7zGjP1QCeb00jFRbcJ3NyePJFhM1cWTmCJodHc42ZSTnkMNcK7DEmz6sA/B41hXIegN9i5v9KRA8BuJuIbgTwxwB+0tx/L4DrAZwF8B0A7x2R9l6jsZIKmTsbxiMw8zrrdEntrvSe6116kedtuXxWd2lkEwozfxXADwvh/wfA24VwBnBzckJrOB5iAhmb6GbM9xrKeSkIc0k6l0ruO1ISVKnjcJlNXAQAPNlJG/XPlF1DJV+DjDHUmIclzb1A2pqiCpolgQhc8nDJZQwYO7fBVB6E+gklF7X7HrTDyQcZS5LcyLRztIoYYaji3I7RmAg19rg2QvLVLvtBQsXkTuavh4Xrzv4SygYvKm4j9WPCwtNqNZ3/A8LasYWxEUpppFS4qWbgRqDqxA4A64gmhptvdwZzjgZAgiaRgTaOmhXYjVBKI+Vt11wzapatEMTePeQwTZk07YzTRjUJRcS1aSMS1kUoFZzdWg0OQIOfFcq6pS72TE2mKBZoK+silH2eLbthv5DamD11OzYYOGp/lVgCGVgXoUjwnUK/oY+tTDqI/oxIZ5VcfLkdn5NQbDBw9N7Y2673FnyzO9uwko1o7ebWptn1kPUq7QV8JavCXPPiZ8C6CSVzLUR2WorVwtu5uuVRnMcpzbQArD5KGGlpn+3tY5MCbT3VrD7WPqPdKCoR6yaUWpCwUnRDGqZQCgd9Q7tVhOJBDrTPLo6ptIyY5h3S1rXhI7FuQplrzseGxTHp60vZrtG3ApgRFbKU9krW4j7dA/OZ6+shFKlAUnYJr8E23UgtGbG9qYqmlaldqjfAL6y9qnlixrq/HkIJFUqKCpjyfGnUQGpjsHZCnFr+zO1Cx2yklTUvZsJyWA+hhFBii8Y5sdaGWbCcihVB201rIpxFzXHSiJlB7WWlbG50WeU4YTnsB6GsbUh37ZrKSAxGRcZGlPlI69OI+jbGzHWK7ao20pyrrSqt/7D02E5joevbLmWLwLYMRr2CjAd7hw4Y52brkk0aEQns4jYlUtIIyjTRlm11ayihc3lK9BAHhUwq1N5qWkUxZrZorcpxT6aeWmYRSabPJ4QooRDRHUT0DBE9aoUln19MRCfM/U8R0QmVdNJRpHaBrEXDWLrGVV5GucUTfc609tB9HbFlOvDnHIVKQU+mnlrGaX6nRGg0lF8HcJ0TlnR+MRFdBOADAK5Cc3bPB1oSSoZUOLWjttq2JgT2EolOXLZaVa8NCXVGPUfEim+VCOz3UgJRQmHm/wHgnBOcen7xtQBOMfM5Zn4OwCkMSSoPB7CxrrYyS+DILFXNXiLDRwa/uZkOtjuaxPEMM9I6p+LVbs6pCxWO8qSeX6w+15iIbiKi00R0+tlnn5VuyBR5f5BcHyous9mWKoQcrhM5KJOwRNpL+FASIUnHgfBhoO8o0t0NowQ8kNjjMhu4AqLT301ZuKazL/JgXPrGaJtc1fD7BD7IXEJJPb9Yda7xAKFRHhul3lBOPLEdcDYkQT03xAepgVgaSJIDNjIdIUW70vDXrJjI9MkllNTzi+8DcA0RHTbO2GtMWBjSKI+EUtOKUws4xvBV1Jz1wCaR4JnAFtRFXEI1sFWMzDkwVVSJCdWk6MQ2IvptAG8F8ANE9A00ozW/iITzi5n5HBF9CMBD5r4PMrPr6E2D76XOOYltmzCnhrpIqHGbeu+VJpSBdka1r06UeCeR56UkSh0lmoqgW2jCuklcccU/duwYnz59enhh4qGvKGpw4u0DhBaobfdtQ+0IpX2IWpOJ53s9RPOmp4BbRcMTxukMMx8rkW7dM2Vd2OrmUvqjnf6GcbDKMGhJtCRh9SNtr9+RiTMDtItnDi8oMxjhCXRaJPmOlhzm9mBda3lqaMQ1yLBnCCp8ApmID2nW4ZQwewKmdjBmZdqhWAblFIgvuIbH0eZKYl0aioQJCycY79wjS3s8YhRtZ+YGZoxyimaTSU81ynDcj0m7NNzhpkqGjevAHI3M7eFKDxNrJxfVUiEnQkcW4gWnV45pJUDwXRElmhZjyr6W9+YuEJxoUsy6CWWpdT2+dHIrTy2VrnakvF9P79spDCvcTLzLkqIcBoQ5tSZvUD+hWAVA8C8Um5xLXDV7zNyXFZsvk4muMWPGmCztmh2tlj9FRktpt2N8MRMPaNRPKMCuMpBQUFZFmbyn99nDB0jDmCyrioid9XwqZPNCzJwqHedMGoSYZkHUTyjSQPqSPXyJlxCaHr6PKJA3+/UHX4Gt0dqT3Uo2ninfX46cYwYODpxT1i2UpSa1TT33ZI+0nEE9zhwZyZrpbt0423ShVkDbyek6QVPiyklben5Ov6JB3YRiLw4sWTPmHvLdAMBTXIEyXHL+YhLciXUtpOmqsTqTmlmN+T1jPa2bUFqUNnVKqZXV1/RlEKrT3tEH5769gU0iLTvOkUmXvGzNybpWWpS6CeXo0aEPZc6X0mLvankhRMpFakvBndGI+s71WLm3jaKGA+pdVcquq7U48YUJeqVFWc/U+1peyoYdIu9gwPt2p6CJLzbS0jaKOeaU5M54lRxAqTNvtX5DX7wz+h3XQSiSLWr/ngMbgWUhVmxRh2tg7cxsGFPfxhZAike6ZXAPMXe8MqFCV7fJ06ItEOhfbBVq8IY+BFOmC/Kt1q2ByEfKQATZfNN4nFP9h4H4XCtsCof3OggFANja/VzB2kXV4By1dEMfrkPFvtR1FHWVYa9TGtHyBpMulS25R0QLDAHnoH5C8TF7bYVbQ09aM6xG5GqPta6rmVwuZ77NMP2UCTh1oH5C0TrwJlg52cV7UDHRfJ12p7WDDCL0RlrEar6kzzATuUeR/hwR/QkRPWz+XW9du9UcRfokEV1rhV9nws4S0S1uOlmwx9anGk5ewUucDGPzHmgIA78huDaLp4FUnwrUMVXRrohIWuQeRQoAtzHz5ebfvQBARJcBOA7gzeaZjxHRISI6BOCjaI4qvQzAu8294xAaGtuwOJIPIa/xfVoy9dYGzZz2WpB7FKkPNwC4i5lfYOavodn9/krz7ywzf5WZ/xLAXebeMnAdXs5swB7GrKvIXZ+xdoTyGlgCPKo9LNmYPPldYfueHWN8KO8jokeMSdQefD79UaTDB4ZhoTGxlFoRmmh1kGoXB3webTnYRD4XpkgvIb5tasIQuYTycQA/COByAE8D+CUTLs4kCIQPA2NHkQKy5iD5T6asbAdJQzGIZnnulXwl03PXCGiSr3R0aklkzZRl5m+134no1wDcY36GjhxNP4rUL4DssMqcup2UribuPYSY3ZwyrAmpU+BrBRGI6yC4LA2lPdfY4McBtCNAJwEcJ6ILiOhSAEcAPIjmxMAjRHQpEZ2PxnF7MltqSRNJVX/XXIEWRK+Y116GU00XnRv2pM+FkXsU6VuJ6HI0ZsvXAfwUADDzY0R0N4DHAbwI4GZm/q6J531ozjM+BOAOZn4sWVqfkzU0XLz2XrQybEW5IYT1HUUaWjk5F3lsJDUOY8uvwvJf6gzjEjiYR5HGNjhaopIdQMdsEJryyCkz15ytiEzakZ61kklprIdQfKM47T+pkvlGZZQb93iv2UOlBwmZBDp4bIptABbCRiR9rIdQgH5Dtp1pviXevlEZ7ZLx1Gt7DHvtiRee673gPSKTotgTbXddhBKbZ7JtMzAZpmjTsYlh0Yljmhm8a0FqAVeav3URirTWw/cifAVe49YHNWKGxhozF6LmxEHWJCvN37oIxcVBrlBTY8qy3Qh9b7FuQvEhVmE3slGhMznajYBS165Y7yE0IW7UmpjQ/KMNs2OdhBIbIi7Vgx6QSulr0J3J0W4ENGY7gtBtOSMlMZN36zQWwToJBSizPYHv2bmHhRcmrmCDLuUv8c1wz40/Nk1gwyJYJ6FIDb7E7lZL9Wq+Ye8ZMNBO3Pk6GvNkjNy5ZR6b6LhhEayPUNpRmpSVxRrU0LstsGNZe5KfeESJtHWjpM0s0ZAnSJNA1e9xUly+wvVsfYQClB2zr3XfzhnkIQRmGfduLN/Ishy8E5Mso55Vuz4Ul69wPVsfoeRUqtBcldqIBJhFO+kWs4UcrhM34qSJa/uwzUAi2vKxy4morJZSWuNZH6GURK0VdAa5goeWw1S0CRvxGrSBUsg1pdryaT+bd5I5hB9JoxTWRyglK3gNfpOlENHOehVtaWd1JchpxK0maJOCJp41+HMkrItQSi6Pr9XckZBKfGN9Rrapc0BIVzQvnEad05tLZw0xeEAWEnm46bm/ayScdRGKhgC0CwRrJ5MxDVo76U9Iw53NSgssCnRlmhr25kh2oy1ilvFQM7E/U/dTqd1UXBehjEGlG/SE0NsyILeBhZ4T5vMwo9cIpqi8qjhnfEe2xkDOX8k03L823AefDCEtamlojiJ9LRF9joieIKLHiOj9JvwiIjpFRE+Zz8MmnIjoI+bI0UeI6AorrhPm/qeI6MQkOZpq17AxyDBZeqfujXGOjmiYtfSEczSYllTcRj912j4Cc4nHlbVW80ejobwI4F8y85sAXA3gZnOM6C0A7mfmIwDuN7+B5rjRI+bfTWjO8AERXYRmg+ur0Jwk+AHrgDAdNA1Ts7Zj7i0MFI26J047lOs+lrO7/xSazcyYi9hqIVAViLyEsyQ0R5E+zcxfMN//HMATaE79uwHAnea2OwG8y3y/AcAnucHnAVxojt24FsApZj7HzM8BOAX5zOR5UJnZI038HfQ6CZpK8KgLyWcSE+iAYs4Ga6cV1ThsM7UiYkk66IuIXg/gRwA8AOBVzPw00JAOEb3S3DbqOFIiugmNZoPXNQH6xXqaTZdW1FCYALGeKPIw2jc9d1nNkF5NO9P7fCMa30rNUDtliejlAH4HwE8z85+FbhXCOBDeD7CPIj16VF/J2go596ZLU5oGUzUwY/IFoy+5Z4lWpokxVSOdqmy08dbiPwGUhEJEL0FDJr/JzL9rgr/VniBoPp8x4b7jSEPHlI7DkppHiXQDe+WWcoF0v7vNktLW76gaY0V+lxBKN8CSRKWZnxIKXxqaUR4C8AkATzDzL1uXTgJoR2pOAPiMFf4eM9pzNYDnjWl0H4BriOiwccZeY8KysZL6G0do5TRTP5/KTHuj9Dh8B2mU2p2+wpe09AhJKL3eNPsAfCS29DCyxofyFgD/GMCXiehhE/ZvAPwigLuJ6EYAfwzgJ821ewFcD+AsgO8AeC8AMPM5IvoQmnOOAeCDzHwuW3IisOVbWZl7RI1mAZ8dkJhJxfybnvYimY1jCncFL2Vuf0VoONoNj8nmzklZ2veyvqNIgXlNnCnTCsRtX+p8zW3lqfiddQg5yGd6f9Is1BoanQ2NNjG1vAfzKFIbczaoKdMKLc5j4bt22Njdda0EJKdMaAvOkKwzkYk0nLqRybRYJ6FoUKHtPgZJ+7Gy3SNzuaLQ7E+inSczoZ2fooWkylF6Or42vVqdsC6S5qGsCqV6wTnUc2Uag5Wv2vko5JvQkoBSpot5Zsoh3JS4U+UoLXdsLc/a5qVsGorv2Tm3hkxcANiNBGizOOV8liRBppFFWrW7lh7dhb2uZy0kYmN/CSW34momyE2FSJqjplivcQKeNnlhGFgqpzWYD2vTSFzsF6GUaDTSfiGlUShecfHgRGllw5O+O9w5Fu4+JqF7am+soTkmtWNdhBJrHKV7yomnvvvgXVSck3/tOqip0HMQOz6gQlhDQyuB2okQWBuhzDn3ZMz9I7WC1uIaZDe2baM+eFL4GnhMc8hFaObrQSGbWrAOQpm7VaQSl7j3gF6bKHI0gqMJSD7lpa2fkgitedGYP/kJr2ttzdyom1DOnJl/VmyJVucs7ottD8ss9LLSvSNGgXyP9syrUr6dQo04ND19sQ2cPXVx6fVBIhboQeomFCBOJiULrdTojrNHa6l1dikRMfrbE3RyuKuIbRfLQn4W3wbOY0wk3yjPXA29Cn9H4tymEqibUI4e3X1XzArNQmEW7+2U5mm86ZEl3OebP9OqIYIQXdBcPRrRgDjsPV2n3Bx7dLz7ZDdOgLoJxUbOJKqUeMeAqDMZGL4hmgxoZbPJy3wXkxfi6/VQM62xIYtAbISGdaswIYBFev0pkHp8hxbrIZQWNay0dVorgftieYdpwlGW4qFWhDZiAg3P1yGnQs3U87IpK9/iPRclK763oRfOexXmjoDYkRwlsD5CWVDl7JkzFjrlSVuRhDx0HNS+4NR8isuTm+9MppK7i/va3dsIxYg6VEG1u5G112yzKFrxFWzsbeg1dFIzYNLRL4P1EcqCL79tg75pJ2rHZsh8axMpuWq3NYNcbcSkUbJIpZ3bJS0jtubGPR9HHAXTrH6uDKU1ghxHs8bxnYv1EcqC6G101DoWKbLZcwi++SpJozkJfpaM+McgNtXddsYCaY2tI29niF5i/JxGPJUfJLRbW258OSump8pf/YRSgVfd3TOocb0yQDz+7F8uuF9JPKniSKmYksnjXaKvLRRL22vfzcAkzeiFU4adUxtnDT6WxUyewFGkP0dEf0JED5t/11vP3GqOIn2SiK61wq8zYWeJ6BYpvQEqUGOlibA7n4enEQgNwqvea7JYYnp/JnNpNlXWxJF7kFU07oI+oEHclpw+baB041zzXBnNBkvtUaRfIKJXADhDRKfMtduY+T/YN5tjSo8DeDOAvwbgvxPRXzeXPwrgHWiO1HiIiE4y8+MlMlIFJLPCUm96Iyor2F6h2AbIZG0oHkkH2KnkqvQicWshpacNK40aNJhcRAnFHIHRnhD450TUHkXqww0A7mLmFwB8jYjOojnLGADOMvNXAYCI7jL37g+htLAJQ5z7kVFlnLU/XVAKOaX0+m4vSQziuKYiramR8hu+X08mvhG3VMyheRwEJPlQnKNIAeB9RPQIEd1hHXw++ihSIjpNRKefffbZFPEWRc8XGNsoKbeettsn2sPLEmG5IyCpyUgqd4BM2me8qrqCiOw41E7DBczhtW7cNBfGHEX6cQA/COByNBrML7W3Co9zILwfYB9FevHFWvHmQaCRsnHVziZHYBmxd/jaN4MWsqYw1uOsHdL0Lawj5y8n/ak83lJ5Sfk4aCSj2qRaOoqUmb9lXf81APeYn6EjR6c5inQuSOtjIhsY+ZbWZKdvpxlaQiwlyMP+NTjZjGlHKoKWYWsSkt+jcZbKc1G063bcYeUUM4TBQpeVB1fe3AWKtWAqX1D2UaTtucYGPw7gUfP9JIDjRHQBEV0K4AiAB9GcGHiEiC4lovPROG5PlsnGQvCxhNXQe6t8Fb1lsEezJ6f1InduS6gn0ZEXJq/J4jptB4+6o2CQzRqbbGLbALgzaN2wYF7ETOjui82piSZTmaYyFdmNOYr03UR0OZo+4OsAfgoAmPkxIrobjbP1RQA3M/N3AYCI3ofmPONDAO5g5scK5qU8ckdjXHOkm8CSrv73xPF1uNK8i5DsKfkijvpPercnmDix7Qlcs0IadXI1JDcOn4zd9ZlGy2rWVkpinUeR1gzXJMks3yKVUJl+qd4zpl2E7pUQIwtRK1KYURv62I4inQupDj23AY8g6zFzPrqvkTimcBr24gvYXq6zNbSrnTQpzvXJuOn7TKCpHaU1OmLnlGcjFAkJntS+E9Op3FMNMoQitWQIiZ9cyUgwq6Tkez4efRqubzvmE+lpKbZsNCQe32zXlBEkbXnlrK2ZGkkzlEdif48iHYMIkWhHKPoP7bSX0VaRUr7QdWBoSoTTHGoEkl8jF7nP7kZyTNl6SExrjvk0H9fMqo00tJha7k1DyYDW4z/YdEkIn8KFJfVIUm9ci2qeM1fFDbNHj3xpSNqKK0NsLknIf7Nh01B60GgLyXMhFJrIHE7EFIepF5Hp92MQ0wB8Q81SHL74Q6NGdhrSXJnQ0PZSqNH5vGkoiUh9gSF3TKr6XLpXtBsHW38glh2qETLJrdyh6eya51yfiO+6e6+bfzcO6flaNJNa5HCxEYoFjfkRe5EEGmwk5rsvtQEy0vdOcXtYlRyByWxu3CV7SM0w86Bxk9+XJY3s+LQSrRPYjc8NC2lQpSBpUqH0Yu+/JFZLKJO8sMSJZ9KLYrBuF8hUTafrNdMrRlv5WrILVfqx5drTdJTwmRNSIx3Ix+kN2Ecqbj4k2PnzjRKVMEW0I09an5hqVK4AVksoSSMULoKbwqbJULSHJrk37dnyHHcw2s+68ophJDeMEFwzItRbek0ryGUo/S7RgYTKyKehuCaQKAcN8+Yjm1xZtc9oNOgpsWqnbK/QM/cFIYL/1QXiHNML7bY5oKaD7W0WPbT3taNKrlyyE1JoVInZiPk8fPL6fmvKspQfoy0HqYxCCJlE7ntr0xmDHGLwPdMjf2qG1qdy5q5WQxkgc/zVt9Q/uK0gjXshve0jC+02FrPfu7ACozRR7SfR0RMqS20599P3DzG75BXSMmN+h+wORUEWUmeiMSejMk281GZ/CKUwgguJebpdw3MQrUSjd9JOlGPiShvTkCTS7AiXWCS8FA3JTTPVfxMzT11tLOb/kEytEKbSToB9JhRtL5nam3JrJinMj0Iuni48h8QmIhNfxU01JXzI7Yk1z6gPZLOQ4riVnrWfcYko5gD2XXOfraGT219CUY0BJ/hdYFWAdFdNUAT7e0gkXy8ZrEgFTBxX3Y45YW3ZUtFoEboRjtQ0bbPPHvHS9uoSAQwavUL7scmgJAmkmm5TYH8JxYavgmaM6gSTSXlxwnaMvs3s20bmjjz00gtpIp7d1jRI6YV9vpuU0Y4uvakOK7LLInGfl45EQ48FzGHtyFsOQk5j/0M68k7BwSCUOTbRUZpB3f2JFand56irKBQ2CYA8h6PvWV9j8GlNoaFfaUi295udA9Sc0wE1I0pjYefN9k8wOFqdpHKL+VtytYiQMz7mqG+ELVtuqx42zkXSal+lWZT6XvrL9FNsKFttb2Usp+pqfRNt4wr1unYj9I1c+OTspWet0m6+NOyaEl8qOvJwzJuUER4pzz7kEmKqD6XvvE6eMRDFgSSUrvFrDolS7olizydJkiVJq/E0YmdrgdwGFdNmXGKQRidcAnE1Ex8xhdKzEkiSVQpPgVcOj+yxERwfwfaet+pRaITHl0YM7rsr7V/RbFL9UiJ6kIi+ZI4i/XkTfikRPUBETxHRp8zG0zCbU3/KHDf6gDnLp41LPKJ0CqhMwxQCCB6h4XF+jIDkewg58XLIxFXpfXK4391G65oHrjw+J6Rk+kjX7TRCssZeuvScJPsgWg+ZSHHZeZG0Gx+ZdPdbew/3yIqGaca0pVDeguU4AhofygsA3sbMP4zmDJ7riOhqAB9GcxTpEQD1tIsoAAAX0UlEQVTPAbjR3H8jgOeY+YcA3Gbuc48ovQ7Ax4joUMnM2CjuNplpM+OQA1PbY5WUBwhXRrcRSc/GHIY+8tGYC30Nxn+/a65I5ozvWi+eiK9MitdHtv5IBPl5aN748iS9r5C8JRElFG7wF+bnS8w/BvA2AJ824XcCeJf5foP5DXP97eYoju6IUmb+GgD7iNKq0e13KvSoo+N2ejRfZbTvb6+F4vE974ZL10Nx2Y1fkrsNlyq7RCRuj+0zo6RrA7kJwdEuO70YebvPdHkUSKu02aDJb4jQ5+54bKhGeYjokDlC4xkApwB8BcC3mflFc4t9rGh35Ki5/jyA70clR5HuHHuKygnzclhwsNHOUTimQrmN3qdKd2n69iqxZXcqpI8MXFPD18CkCq4lK4l4fD2ofV07mpnTcFzSjpmRoXh8jd8t45Q/KX2fSegzr9RYYtiYmb/LzJejOe3vSgBvkm4zn5KEHAh305r0KNLORE08f3jQCHj36etdNfD12tpKEbtP0wPHZNDKEmocITNOvM7DhhVt4JYry51iIfXgMS1NA6lD8Gmcsedd0nfvaz9DBF1Se85B0jwUZv42gD8AcDWAC4moHSWyjxXtjiI1178PwDmEjyidFS2pJJ3zWwGahYy0IzNlr6TRPto4cgjKZ4a5skkySA0latoErneyOH7yUFn5tMIUxExHH0Iahu/5kGwuIUXJpbBvUDPKczERXWi+fy+AHwPwBIDPAfgJc9sJAJ8x30+a3zDXP8vN2KzviNLFIB1VDOSTSKxRS2FzEJa2p2xlipGGRAqxnlkiDp8MIY0vZKKFNKJQw+r18CM7d9fnEtOu7LJz8zSmbrTxzt0hauahXALgTjMi8z0A7mbme4jocQB3EdEvAPgimvOPYT5/g4jOotFMjgMIHlFaDM4kNN+cNCl8ig2q7fvdF9xV8MT1RFLvKqWRWiljDVDrd7BlsvPa+UYsQtD0tO69UtquliWVs9qEnKD9SUTpfnffmw1NvXPvWcrs2eujSNu22kyLT2u4U2Pq3iMUf6y31vTmmgobu3dM/kNx2htY5aRXQkPQxi2lJRGOdE2KO6kzMQ2k5FGkB2KmrLujVqHjh8fJlKLp+GQM7Iatbcy+3j52XSX3BL2kSj5GcCRMS+ahhtkzjdyZyoE6FZJb0jBcjctHNPa9Uj4GeZio4q+bUDyF0ga7w73tRVsrq0hp8aLZmNqpSNSaDqaiefbFDfV8PbSrD6177TimJAdJ1tRnB+Bhg9Tmwb1XbLzcf2IgX6S9hjQUiVhcGWLmYtQpO1HFX/dqY283QH4H64iCnNIuDY0ItBtTA1Z+DGO2Iz/u6mW3QUR75YB5kNIQ1b4KqRFaZSA5WFPLX2pYuTJKTmDvc6yvZinypJiZ/U7U40i0Pwth3YTiw0S7lE3l8yhlRtik4trSvlGXEjJJiI1uxNKWtISceHzhoWtAeFQmVx6NbL70cuNg7DTP3v2bhpIAa67GmhCqxOzTuBy4ywTa724PK1XGsYQZ0iqk0SkJ0mgSO3/2fakyu2ZfzOyS0s7VmFwzMnaPJENIVh+CZTX3PJS1wJ6oNuXoSWloPfKpSleoJ431zmME0IwQaeMJNVzpWq7p4C0PIdgmOFvOXhzeGZNx+Uprjr6ym8p83w9CIeqN2qwBtl8j1GO24VpV2z5ZUNPAtI6+sXvT2muectR4t+E2jmmd1qOSz0MyPu1OcpJ213lIeITdmrAxcqocrs79UthUne66R3kE1Gjq5NrLkgkRmmcgHTilTdeNv/jITnPCFNoNosfGT6BoVlPTGJRrhETd9+OWXY9MAu9tSZTWVPaDUHon79WDUo5Nn3Oy+07c21Iyt8KOdYJK6DVqs20jEUCejbNLVvDcuCS/is+06s1pcgjEhW8IPxWuT2RsXCXLfD9MHsxj6uQWvNqsCDzvDgO2Kr/txbfl1KrWrhnhcwCOgeTQlPwJcXMuLNdYmUOOUq/D3AqWzLBQWgMzivLqcU0az+o1FPdcm6Qp9glTZnN7AcnDnjMy0YvT5hYCwJZDmll8xk47xQTLIdHYkCs4j5xjeXJNCt87k+Lp7nO0vd41CHnzTa4EReUYxM27wQXVmU4jNJOpsHpCGRZ8QgFHyMR+YVO8OG2FEIcrTcX3TJAVzSXpt3v/VGhlbhIbmg9j/Sm2/8KOzzVVJNKxESpH7bCrdnhX8l2BdSM9Pm1taZJZvckzxtTZ2fbx3iM7bs+1HGehhHaDJ/sUPPuZTg33HCAeHP5UCZdC4Ja545wrVMLUijVkqSFqiVerrSUfp2L+bNk0pp+LsaNHpbBqDaXE3hVTIWSP+677RnJsLcmb5/acmmZj1a6ZgJvY2p/AjgN8FdDVzIIV1XMwea+hEjeKo71eiPxbLaRqLbGyDt1jX/fJoY1jDEo7R5dC/RqKd5LQbuQgKboMzaAk+4d6YU0P3a4TkVRzgmO69fap5E5R0WgDKRqL6zBuycIe3el9WlpNdARLAY2cobglItOYHd5rEdGl86xDZowmXdHxvQDq11ACOyS11SApOo8a7LOtS9uk3uFHt1cXwmOw7fDmt+e6lYa3Z6S2dClq1fSe52G8wzLf3SfFlVveYzSaVH9ZkHzYeMp9z7q3cxNjMnk7CGm/c6F+DUVCxoSTFPNIqlylNZSQOu377kNw1EIIl3rr9ryZ7s9oQZJGE2t40bLzHd4+0aJOW5bYdW1j9JOlXlMCHIUyMW0XY31QJbBOQnF1RgUGWz6aP6mHjjWIOVRKrUnik0ly0NpxA6ZdtyYKD008caSEWjdNfBvIWMUfEFTr26F0Uo3dW/qd5ZqtWmIjkHdeyljSyDUtNRhzFOmvE9HXiOhh8+9yE05E9BFz5OgjRHSFFdcJc3TpU0R0IldogmWjZ06PbV98zEnqFr6vkRZBghrlyu/KFNOsGI0K0u6nIsXRE81U8J6T1wPJoew2EC9hZTbUmPwlIaWR6mdTxdFuOyAQizatmPZaWqPR+FDao0j/goheAuB/EtHvm2v/ipk/7dz/TjQ72h8BcBWAjwO4ioguAvABAMfQGJhniOgkMz+XKnTOmhVfwUqFW9pvooaHHEPytKQiXdeYaz4tTU7M3MXNk7tIuBtlcuPr/TbreRoXWJm1La4fKgStszVFi0g1lUJpiu+hXQMlyKCZABfS2NTvPQFjjiL14QYAnzTPfR7N+T2XALgWwClmPmdI5BSaM451GDFG7L54TYWJPTO2AqUgZKdrbXZfj9h+dmnk+DGEDXxcUmnD2r1qCHJDAbC7jmHefb2rthOIaZhjzSuNpuLKGk2TyevEzm0XU5AJkHkUKTM/YC79O2PW3EZEF5gw35GjqqNIvYhMKU+FtvJJvXyJ3mjM8+73mM1u936uf2QQ78htCuwhY3b+fAQW8lFJ8kq9ubqBOvNfNPVJS9q2bOJ9NPQPuZB8Sz45xp7kMIUTN+soUiL6mwBuBfBGAH8bwEUAfsbcLpU+B8J78J5tTOnq8fAFpBVeadU293m74vvMtGKIHDbu/h78sVxRXSLtTKDWDPIQnZ1OTB5VZ22l151PHSBaKZ1O/sD1wf2JPK0iFSWZzOFXapF7FOl1zPy0MWteAPCf0Zx5DPiPHFUdReo92ziTicdoBHYjTu3VdAnIIzNSGm5jTMlX6v0uXGJ1G3gnsVAsrhYlysBCZ+EsFwg1+C4aDjt1u892VnFgLkzIRGzjsSelReuEtS2pllzU5mws6dKdTgC5R5H+ofGLgIgIwLsAPGoeOQngPWa052oAzzPz0wDuA3ANER0mosMArjFhcYywEwdRKclAUqs1JkaagOEeX4LdqEPwqeEp2k1MG3K1JonztWXlJZrI8HQrR0zLsNOIfbrxSjISaDdr2UOwYtqMZGJx35/bQUj7sSyFMUeRfpaILkZjyjwM4J+a++8FcD2AswC+A+C9AMDM54joQwAeMvd9kJnPZUuuXOMtqeneKCF7vn3ho9A71lA3qhCLr52E1j6j1krsNTYsPxMzYVIQKsvBNWunt4HYiVqam4ZEtHb6wXdODBLKqjOlFL6o1BEaTQeguW9K7PVRpLlwCSSpcS4EyUfh09BcE87FFPnspecMM7vppmiRsXuj9zhbKkgyhORyr6X4U7TnbotxFNSWSx5Fus6ZshPBVlXtBqkhkyXVTGBEb0a736IplbuNmB2FUKbifU4yQ1mGfhVt2vLFNB9U7NpYLTZGJlPUsdJxHhhC0RRc0KOe+KxPhkGcC2/Tb5tJIlKOwZPix9CRy4ze8O0uLf/wce+6Z38XIOIj6vZtMLK0o0zOTGGtn8qXftSXwztrt6cJKUyg0gRQWhs9MISS4hyMOUZdpDp6+4E8C6lIKnLXoAum7zqxfWVobwwV0kQa+RzSEbZB8KVvk4i7jYJLPvb70YyMxZy4oTAgPCqlSbNGE7z67QvcXcXHFGLq874RHp8HX3pWlV5hP5bkA9pdpJ3boA0vmH6KmSg91/wgi1Saw+3J9nU4GzV14WLEfhIR03by4MuT75rP9xLtdLTOEyvtGv161Wsopep6bEgt5HBzey6NaeT1R0Tks3/nqLdtr9wjQ3LIgwTZlHJqkUTcWjeHu67F9z0BvUZvzSkJObXt59x4RO2W5Gv90aw88tXWEZ+WVBrVE0oPI/fLCJkxmnBf76T2vQQqjWiK5MBW8T3JRuvuHCN/7cxnX1LOLFopPLzGZfiuYwQR25uktDag8bdIkEb0NGmFfpfCqghlTD0vVYBSw7e1nxJqqK/iSybYcGgkbo6Jz80N7hNeqyEAQ/Onc6C6vpRuT8uhT6WLPjAqJJkwUbGV71brbPVqLIF42zpnazxdHXR9/hNpIj6silBKYGwBi6MWgmlURB7yVzTX/9FsemQ5+1IcfguSS2/Y1aMhdGXuZscmEltzsZ9jdNqKrRHYafjMWJHAlUjVCLSkYtc1119jk0qSD68gVkEoJVk2VsCutpH6/Fh53AYmOfnEeBxHa7TM3M2sC8w3yYHby3r3JHeIoCNvQySMvrZik4Q4gUzhU5D8IIN4ChaZ1CGFzCHRCexxPo8hxxSsglDmZFlXnSwKhVPWTVdSa73g/nOadH3Py0JOSzi2FtJuM+leB/qO0zac4PhTqL94T4IUh9eX5imaVDM89g67nDn8nkQCwnwet04daKfs3HbgZASWsSNb7/FcuTS1vtQ9hdDNp2v9JIFFgmK4O2HN8s/4NIG5Oq5Qfd6ZsmHHfxuPT2afBja1KbQKQpnbDpwbsVEilWKQqz0s7ZwNoNM6LF9JV1aB2bI9ePbLHWS71WgyTYOYxejz3QQRMF/cAYHuEYEwbQ03dXQoFfUTSqTCz629jEk/Z4gQGCoGg+cDk6KifBFcPBJ5OJOMkkZU7Fmzjp8oPWFrsp/ttG4PSPfs3woK+3a6+FguEp+2oHHU2vVFei6lLFIdxTmon1BKLO2fEDnpaz3++gjjztosRBeX5A9J+4h1oCHwbt+RwT2c6Otydoezwwm0c+raI0UEdAe70lDuQSM1jnS7CELOVp8z2DV1JE3EjicGW/s68E7ZlBmmOXGURqihlCDAQRxJp5jNUA6KTsDnE7A/23uHN3JvgWFWA3EXCsLSWqRp/aFJdA5a5zI5f2PfvfS8XZbRcqDh2dilyaX6tTwAwj2wunead6QI0A/5+kAgc6plgqYQWhMyFZko1qEkOVQR6U3Z/I/75dvda2sa7md7vX1OEHvnGEVfcxFkDA7l2/4KQrchk/2cj2zsMNfv4RuREoeRrXxKEwdLj/isQ0PxoWKHIhD20Guf79pA11iG9nhffQ+kVWArAveTQN2pgyFtwS6LULn4nIrup8+fYA89Dw5pdz67YWVCRzK99TwmXJok2DVokt/L4F7uPxdz0KY4Yt3n3PIjZfUogXVoKD6ULh2pp018G74Kk4tBb2TJ4PobUtLKUcG1PWOIVOxy0fbw7cFgUvxSeUvfd5FapOKsYG58NthpJFa4z0RjxTaekj8kZPK5z0pE4iu3UJnOoaSvT0MptHuYGK2o/+60hNDYYMxZNwbeyhzp7SW5xsonaQSa+HLTax2zvvJ1NR/XnFA5bQVTyP0eNM0SW6pEDu6fmA7SzKwubD5rX08o5rCvLxLRPeb3pUT0gDmn+FNEdL4Jv8D8Pmuuv96K41YT/iQRXZslsW8jTu2h6bns3S2Q0WkpobTGIMcOnkIOLaQGkmIK+hqbdEayhBihdsPFvYuWA9ae2h9IIwUppkvIwV9CltJI0VDeD+AJ6/eHAdzGzEcAPAfgRhN+I4DnmPmHANxm7gMRXQbgOIA3ozmC9GNmJ30vzuCMXrqQqaKFPd+A9NqG5EBz/Qyl4J0C7iPKKcjEKpuY78Qnz5gyaf0krk9FcuSGevzevJbICM7Yd+jKFHI6+wgwpMHY9y0J7VGkrwHw9wH8J/ObALwNQHtQ+p1ozuYBmrON7zTfPw3g7eb+GwDcxcwvMPPX0Byz0R4OJuIojupz4mopQsvzVgqbPNj9EkdM/dYQklrWGuCsGRLz50wQk0jHFyYh9qyvwQ20RWuzo54/VXn0RQgap3Tvt3JfXPfZmuuG1in7HwH8awCvML+/H8C3mflF89s+p7g7w5iZXySi5839rwbweStO8WxjIroJwE3m5wtE9Kh7TxC5M2sV3vocBF7+DwD409B9NVecAEy+yPwXJ4xcUomFy9f7cpFzKYLeO/OlkfTehDLSPVa0bvyNUhFFCYWI/gGAZ5j5DBG9tQ0WbuXItdAzuwDm2wHcbtI+Xeq8kNqwr3nb13wB+5s3Iip2+JVGQ3kLgH9IRNcDeCmAv4pGY7mQiM4zWop9TnF7hvE3iOg8AN8H4ByUZxtv2LBhvYj6UJj5VmZ+DTO/Ho1T9bPM/I8AfA7AT5jbTgD4jPl+0vyGuf5ZbiYRnARw3IwCXQrgCIAHi+Vkw4YNi2PMxLafAXAXEf0CgC8C+IQJ/wSA3yCis2g0k+MAwMyPEdHdAB4H8CKAm5n5u5E0bh8hX+3Y17zta76A/c1bsXxVfbbxhg0b1oX1zZTdsGFDtdgIZcOGDcVQLaEQ0XVmiv5ZIrplaXliIKI7iOgZe94MEV1ERKfM8oRTRHTYhBMRfcTk7REiusJ65oS5/ykiOiGlNSeI6LVE9DkieoKIHiOi95vwfcjbS4noQSL6ksnbz5vwZZaVFMYiy2WYubp/AA4B+AqANwA4H8CXAFy2tFwRmf8ugCsAPGqF/XsAt5jvtwD4sPl+PYDfRzM352oAD5jwiwB81XweNt8PL5yvSwBcYb6/AsAfAbhsT/JGAF5uvr8EwANG5rsBHDfhvwLgn5nv/xzAr5jvxwF8yny/zNTRCwBcauruoQrq5L8A8FsA7jG/J8/XohkOFMSPArjP+n0rgFuXlksh9+sdQnkSwCXm+yUAnjTffxXAu937ALwbwK9a4b37aviHZnrAO/YtbwD+CoAvALgKzWzY89y6COA+AD9qvp9n7iO3ftr3LZif1wC4H80SmXuMnJPnq1aTp5u+byBO018BXsXMTwOA+XylCfflr+p8G1X4R9D05HuRN2MWPAzgGQCn0PTCqmUlAOxlJbXlrV0u8//Mb/VyGYzIV62Eopqmv2KMWp6wBIjo5QB+B8BPM/OfhW4VwqrNGzN/l5kvR9OjXwngTdJt5nMVebOXy9jBwq3F81UroezLNP1vEdElAGA+nzHhvvxVmW8iegkaMvlNZv5dE7wXeWvBzN8G8AdofCgXmmUjgLysBJUvK2mXy3wdwF1ozJ5uuYy5Z5J81UooDwE4YrzS56NxFJ1cWKYc2MsQ3OUJ7zEjIlcDeN6YDfcBuIaIDptRk2tM2GIgIkIz+/kJZv5l69I+5O1iIrrQfP9eAD+GZs+fVS8r4SWXyyztDAs4la5HM6LwFQA/u7Q8Cnl/G8DTAP4vGma/EY0dej+Ap8znReZeAvBRk7cvAzhmxfNP0OwVcxbAeyvI199Bo+Y+AuBh8+/6Pcnb30KzbOQRAI8C+Lcm/A2m4ZwF8F8AXGDCX2p+nzXX32DF9bMmz08CeOfSebPkeit2ozyT52uber9hw4ZiqNXk2bBhwwqxEcqGDRuKYSOUDRs2FMNGKBs2bCiGjVA2bNhQDBuhbNiwoRg2QtmwYUMx/H/eV9fktuOHxgAAAABJRU5ErkJggg==\n", | |
"text/plain": [ | |
"<Figure size 432x288 with 1 Axes>" | |
] | |
}, | |
"metadata": { | |
"needs_background": "light" | |
}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"plt.imshow(rgba.astype(np.uint8)*255)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 48, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"masked_array(\n", | |
" data=[[[ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" ...,\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255]],\n", | |
"\n", | |
" [[ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" ...,\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255]],\n", | |
"\n", | |
" [[ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" ...,\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255]],\n", | |
"\n", | |
" ...,\n", | |
"\n", | |
" [[ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" ...,\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255]],\n", | |
"\n", | |
" [[ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" ...,\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255]],\n", | |
"\n", | |
" [[ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" ...,\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255],\n", | |
" [ 0, 0, 0, 255]]],\n", | |
" mask=False,\n", | |
" fill_value=999999,\n", | |
" dtype=uint8)" | |
] | |
}, | |
"execution_count": 48, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"rgba.astype(np.uint8)*255" | |
] | |
}, | |
{ | |
"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.8" | |
} | |
}, | |
"nbformat": 4, | |
"nbformat_minor": 2 | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment