lab_6/6-morphology-1.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 6. morphology-1\n",
    "\n",
    "http://scikit-image.org/docs/dev/api/skimage.morphology.html"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from skimage.io import imread, imshow\n",
    "from skimage.morphology import (closing, opening, \n",
    "                                binary_closing, binary_opening,\n",
    "                                binary_erosion, binary_dilation, \n",
    "                                erosion, dilation, \n",
    "                                diamond, disk, rectangle, square, star)\n",
    "\n",
    "import matplotlib.pyplot as plt\n",
    "from skimage import data\n",
    "from sklearn.datasets import load_digits\n",
    "from skimage import img_as_float\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "from plot import plot2, plot3, plot4"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "digits = load_digits()\n",
    "imgs = digits.images[0:10]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAARsAAAEYCAYAAABsuVKPAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAADwtJREFUeJzt3WGoHWedx/HfzzTdCiJt7WkJuWVvF4K0L3ZbuJRCfbFEC7GKyYsutCuSF4G8caGi4MZdWBD2hX1jfbOsBFu8L8TUrUJKcVlCTBFhSXtjq7YGTSy4ZhuaKxq0b1yj/31xZuVuci5n7pmZ/zwz5/uBwz0zmTnzPOfc+8tz/uc5M44IAUDX3tV3AwAsB8IGQArCBkAKwgZACsIGQArCBkAKwgZACsIGQIpGYWP7gO2f2L5o+1hbjQIwPl50BrHtXZJ+KukRSZckvSLpiYj48Xb73HHHHbG6urrQ8VCuc+fO9d0E9OuXETGZt9FNDQ7woKSLEfGmJNk+IemgpG3DZnV1VRsbGw0OiRLZ7rsJ6NfP62zU5G3UXkm/2LJ8qVr3/9g+anvD9sbm5maDwwEYsiZhM+u/sxvek0XE8YhYi4i1yWTuSAvASDUJm0uS7t6yvCLprWbNATBWTWo2r0jaZ/seSf8t6XFJf7vTBxnD+/26RfYmfZ11jD6euyanJBnD6UzqPudN+jrrGH08d23/fi0cNhFxzfbfSfoPSbskPRsRb7TWMgCj0mRko4j4tqRvt9QWACPGDGIAKQgbACkavY1qw9CKhn0VtMdQSB9aH/r63Rza30RdjGwApCBsAKQgbACkIGwApOi9QFzKLNihFS/blvGclDILdqwF2Lr6ek4Y2QBIQdgASEHYAEhB2ABI0XuBuJRiXUY7SulrXW23t5QifEY7SulrXRntZWQDIAVhAyAFYQMgBWEDIAVhAyAFYQMgBWEDIAVhAyAFYQMgRe8ziJto86vyfX3tvo/jDvG0C22eAqOvU4z0cdySTqfCyAZACsIGQArCBkCKwdRs6tYZ6mxX0rWf6tZKFt23bl/rbtfXt+MXff1LuvZTk+e4zr51+9rkd64JRjYAUhA2AFIQNgBSEDYAUgymQNxm8aukiU6lKOk5abOoPcQJjF3julEARo2wAZCCsAGQgrABkGIwBeKulXTdqDYLeHULv8tUSC3pulFtFubrvoZ9fRjAyAZACsIGQArCBkCKuWFj+1nbV2y/vmXd7bZP2b5Q/byt22YCGLo6I5uvSjpw3bpjkk5HxD5Jp6tltCQibrh1/fhdHxOz2b7h1vXjd33M7cwNm4j4rqRfXbf6oKT16v66pEMttwvAyCxas7krIi5LUvXzzvaaBGCMOi8Q2z5qe8P2xubmZteHA1CoRcPmbdt7JKn6eWW7DSPieESsRcTaZDJZ8HAAhm7RsHlB0uHq/mFJJ9tpDoCxqvPR99cl/aek99u+ZPuIpC9IesT2BUmPVMsAsK25342KiCe2+acPttwWACPGDGIAKQgbACk4xcTINZkdyiziYWnyenGKCQCjQdgASEHYAEhB2ABIsRQF4kWLX30VSJucv7brdpReNF60fX1dkK/Jeam7bkfbx2RkAyAFYQMgBWEDIAVhAyBF7wXiPoqhfV2ka5ZSisEZx+yjGFpSkbuUYnBfx2RkAyAFYQMgBWEDIAVhAyBF7wXiUpQ+M7bNQmfdAnlJxdU29fVhQF1tfoBR9zVkBjGA0SBsAKQgbACkIGwApCiyQNxHEbKvYmjXx6j7+KWfTqNNfc0g7/oYdR+/rwI5IxsAKQgbACkIGwApCBsAKYosELddrB3DrNexartYW/rs4GXGyAZACsIGQArCBkAKwgZAit4LxBmzOa9/vGUqGDc5x3HbhfqMWdrXP94yFYybnOOYU0wAGA3CBkAKwgZAit5rNsv0nnqWrq8bVVJ9qqS29KHr60aV/rfEyAZACsIGQArCBkAKwgZAirlhY/tu22dsn7f9hu0nq/W32z5l+0L187bum7scIuKGm+0bbnUsuh9yzHp9Zr3+dSy6X5Y6I5trkj4TEfdKekjSJ23fJ+mYpNMRsU/S6WoZAGaaGzYRcTkivl/d/62k85L2Sjooab3abF3Soa4aCWD4dlSzsb0q6QFJZyXdFRGXpWkgSbpzm32O2t6wvbG5udmstQAGq3bY2H6PpG9K+lRE/KbufhFxPCLWImJtMpks0kYAI1BrBrHt3ZoGzdci4lvV6rdt74mIy7b3SLrSVSObWuZvAs/S1zWyMizzN/xn6esaWbPU+TTKkp6RdD4ivrjln16QdLi6f1jSyfabB2As6oxsHpb0CUk/sv1ate4fJH1B0jdsH5H0X5L+ppsmAhiDuWETEd+TtN2464PtNgfAWDGDGECK3k8x0UcBr6Si4U5mh3b5+G3vO0sfhcmSPgzYyazvLh+/7X3rYmQDIAVhAyAFYQMgBWEDIAVhAyAFYQMgBWEDIAVhAyAFYQMgBWEDIAVhAyAFYQMgBWEDIAVhAyAFYQMgBWEDIAVhAyAFYQMgBWEDIAVhAyAFYQMgBWEDIAVhAyAFYQMgBWEDIAVhAyAFYQMgBWEDIAVhAyDFTX03wHbfTWgsImptt0x9bXvfUtR9DZepr3UxsgGQgrABkIKwAZCCsAGQovcC8SwlF9fGUOQtScnPZ8m/h0PEyAZACsIGQArCBkCKuWFj+xbbL9v+ge03bH++Wn+P7bO2L9h+zvbN3TcXwFDVKRD/TtL+iHjH9m5J37P975I+LenpiDhh+8uSjkj61zYaVXLREO2iCLs85o5sYuqdanF3dQtJ+yU9X61fl3SokxYCGIVaNRvbu2y/JumKpFOSfibpakRcqza5JGnvNvsetb1he2Nzc7ONNgMYoFphExF/iIj7Ja1IelDSvbM222bf4xGxFhFrk8lk8ZYCGLQdfRoVEVclvSTpIUm32v6/ms+KpLfabRqAMZlbILY9kfT7iLhq+92SPiTpKUlnJD0m6YSkw5JOLtKAZSoQLlNfZ1mmwv8y9bWuOp9G7ZG0bnuXpiOhb0TEi7Z/LOmE7X+W9KqkZzpsJ4CBmxs2EfFDSQ/MWP+mpvUbAJiLGcQAUhA2AFIQNgBSEDYAUhA2AFIQNgBSFHla0LpmTZxadOJcm4/V9LiztNmWIV77aFZbFp041+ZjNT3uLG22paRrmjGyAZCCsAGQgrABkIKwAZBiMAXiugXcOtvVLYa1XTRu8niL7tuk8NdX0XzR13XWdnXb23bRuMnjLbpvk9cmo2jOyAZACsIGQArCBkAKwgZAisEUiJsUDbt8LLSvyYcBXT4WmmFkAyAFYQMgBWEDIAVhAyDFYArEs7RZwC1pZmzdfUs9ZhfaLOBmFIPbni1e6jF3gpENgBSEDYAUhA2AFIQNgBSDLhCPQR9F2JIKv8ukj1nKJc2MZmQDIAVhAyAFYQMgBWEDIAVhAyAFYQMgBWEDIAVhAyAFYQMgBTOIR67tGaTMPi5X268NF6kDMEiEDYAUhA2AFLXDxvYu26/afrFavsf2WdsXbD9n++bumglg6HYysnlS0vkty09Jejoi9kn6taQjbTasTbYXupXctpJOHVCSiFjoVnLbxlKUrxU2tlckfUTSV6plS9ov6flqk3VJh7poIIBxqDuy+ZKkz0r6Y7X8PklXI+JatXxJ0t5ZO9o+anvD9sbm5majxgIYrrlhY/ujkq5ExLmtq2dsOnOsFxHHI2ItItYmk8mCzQQwdHUm9T0s6WO2H5V0i6T3ajrSudX2TdXoZkXSW901E8DQzR3ZRMTnImIlIlYlPS7pOxHxcUlnJD1WbXZY0slFGtBHMXSsBbjS9VEMpbBejibzbP5e0qdtX9S0hvNMO00CMEY7+m5URLwk6aXq/puSHmy/SQDGiBnEAFIQNgBSEDYAUhA2AFIQNgBSEDYAUhA2AFJwDuKedT17mdnR5eh69nLps6MZ2QBIQdgASEHYAEhB2ABIQdgASEHYAEhB2ABIQdgASMGkvpHr4pSqKFMXp1RtEyMbACkIGwApCBsAKQgbACkoEPesSRFu0YJg3f1K/xbx0DQp4C76WtTdL6Pwz8gGQArCBkAKwgZACsIGQAoKxAWaVayjWDtOs17Xsc7SZmQDIAVhAyAFYQMgBWEDIMVSFohLL7Yu+/WF2lR6sXWZrhvGyAZACsIGQArCBkAKwgZAit4LxBkFrJKKZNdbpgJhRmG65OL3shf+GdkASEHYAEhB2ABIQdgASOHMAqLtTUk/l3SHpF+mHbgbY+iDNI5+0Id+/XlETOZtlBo2fzqovRERa+kHbtEY+iCNox/0YRh4GwUgBWEDIEVfYXO8p+O2aQx9kMbRD/owAL3UbAAsH95GAUhB2ABIkR42tg/Y/onti7aPZR9/EbaftX3F9utb1t1u+5TtC9XP2/ps4zy277Z9xvZ522/YfrJaP5h+2L7F9su2f1D14fPV+ntsn6368Jztm/tu6zy2d9l+1faL1fLg+rBTqWFje5ekf5H0YUn3SXrC9n2ZbVjQVyUduG7dMUmnI2KfpNPVcsmuSfpMRNwr6SFJn6ye+yH143eS9kfEX0m6X9IB2w9JekrS01Uffi3pSI9trOtJSee3LA+xDzuSPbJ5UNLFiHgzIv5H0glJB5PbsGMR8V1Jv7pu9UFJ69X9dUmHUhu1QxFxOSK+X93/raa/6Hs1oH7E1DvV4u7qFpL2S3q+Wl90HyTJ9oqkj0j6SrVsDawPi8gOm72SfrFl+VK1bojuiojL0vQPWdKdPbenNturkh6QdFYD60f19uM1SVcknZL0M0lXI+JatckQfqe+JOmzkv5YLb9Pw+vDjmWHzayz+/DZeyLb75H0TUmfiojf9N2enYqIP0TE/ZJWNB0p3ztrs9xW1Wf7o5KuRMS5ratnbFpsHxaVfaa+S5Lu3rK8Iumt5Da05W3beyLisu09mv5PWzTbuzUNmq9FxLeq1YPrhyRFxFXbL2laf7rV9k3VyKD036mHJX3M9qOSbpH0Xk1HOkPqw0KyRzavSNpXVd5vlvS4pBeS29CWFyQdru4flnSyx7bMVdUFnpF0PiK+uOWfBtMP2xPbt1b33y3pQ5rWns5IeqzarOg+RMTnImIlIlY1/f3/TkR8XAPqw8IiIvUm6VFJP9X0vfY/Zh9/wTZ/XdJlSb/XdHR2RNP32aclXah+3t53O+f04QOaDs1/KOm16vbokPoh6S8lvVr14XVJ/1St/wtJL0u6KOnfJP1Z322t2Z+/lvTikPuwkxtfVwCQghnEAFIQNgBSEDYAUhA2AFIQNgBSEDYAUhA2AFL8L65jWONWsMHoAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "img0 = imread('images/0.png', 1)\n",
    "imshow(img0)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPgAAAD8CAYAAABaQGkdAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAACtRJREFUeJzt3V2IXPUZx/Hfr6ul9Q2ltVV2Q2NEAlKo0RCQgNCYlFhFe1EhAcVKIblRlBYk9q53uRJ7UWRD1AqmSo0KIlabRcUKrXUT09a4saSDJdtoo3TFl0JC4tOLnZQ03TJnM//zMo/fDwR3N8P+nyF8PWdn55y/I0IAcvpC2wMAqA+BA4kROJAYgQOJETiQGIEDiRE4kBiBA4kROJDYGXV8U9sp3x530UUXNbre+Ph4Y2sdOXKksbVmZmYaW+v48eONrdW0iPCgx9QSeFa33XZbo+tt3bq1sbV6vV5ja61cubKxtebm5hpbq4s4RQcSI3AgMQIHEiNwIDECBxIjcCAxAgcSI3AgsUqB215v+23bB2xvqXsoAGUMDNz2mKSfS7pO0uWSNtq+vO7BAAyvyhF8laQDEdGLiKOSHpd0U71jASihSuDjkg6e9Pls/2sAOq7KxSYLXbHyP1eL2d4kadPQEwEopkrgs5KWnPT5hKRDpz4oIrZJ2iblvVwUGDVVTtFfl3SZ7Utsf1HSBknP1DsWgBIGHsEj4pjtOyS9IGlM0kMRsa/2yQAMrdINHyLiOUnP1TwLgMJ4JxuQGIEDiRE4kBiBA4kROJAYgQOJETiQGIEDiY38ziZN7v5x8803N7aWJG3evLmxtSYnJxtb66qrrmpsrampqcbW6iKO4EBiBA4kRuBAYgQOJEbgQGIEDiRG4EBiBA4kRuBAYlV2NnnI9mHbbzYxEIByqhzBfyFpfc1zAKjBwMAj4hVJ/2xgFgCF8TM4kFixq8nYugjonmKBs3UR0D2cogOJVfk12WOSfidpue1Z2z+sfywAJVTZm2xjE4MAKI9TdCAxAgcSI3AgMQIHEiNwIDECBxIjcCAxAgcSc0T5t403+V70ZcuWNbWU5ubmGltLkqanpxtdrymXXnpp2yOkEBEe9BiO4EBiBA4kRuBAYgQOJEbgQGIEDiRG4EBiBA4kRuBAYgQOJFblpotLbL9ke8b2Ptt3NTEYgOFVuS/6MUk/jog9ts+VtNv2roh4q+bZAAypyt5k70bEnv7HH0uakTRe92AAhreonU1sL5W0QtJrC/wdWxcBHVM5cNvnSHpS0t0R8dGpf8/WRUD3VHoV3faZmo97R0Q8Ve9IAEqp8iq6JT0oaSYi7qt/JAClVDmCr5Z0q6Q1tvf2/3y35rkAFFBlb7JXJQ28NQyA7uGdbEBiBA4kRuBAYgQOJEbgQGIEDiRG4EBiBA4ktqirybqo1+s1tlaT+6A1vd7U1FRja11wwQWNrdX0fnJdwxEcSIzAgcQIHEiMwIHECBxIjMCBxAgcSIzAgcQIHEisyk0Xv2T7D7b/2N+66KdNDAZgeFXeqnpE0pqI+KR/++RXbf86In5f82wAhlTlposh6ZP+p2f2/7CxATACqm58MGZ7r6TDknZFxIJbF9metj1dekgAp6dS4BFxPCKukDQhaZXtby7wmG0RsTIiVpYeEsDpWdSr6BHxoaSXJa2vZRoARVV5Ff1C2+f3P/6ypLWS9tc9GIDhVXkV/WJJj9ge0/z/EH4VEc/WOxaAEqq8iv4nze8JDmDE8E42IDECBxIjcCAxAgcSI3AgMQIHEiNwIDECBxLz/NWghb+pzeWkBTS5xc+uXbsaW6tJ69ata3S9JrdKiggPegxHcCAxAgcSI3AgMQIHEiNwIDECBxIjcCAxAgcSI3AgscqB9++N/oZt7scGjIjFHMHvkjRT1yAAyqu6s8mEpOslba93HAAlVT2C3y/pHkmf1TgLgMKqbHxwg6TDEbF7wOPYmwzomCpH8NWSbrT9jqTHJa2x/eipD2JvMqB7BgYeEfdGxERELJW0QdKLEXFL7ZMBGBq/BwcSq7I32X9ExMua310UwAjgCA4kRuBAYgQOJEbgQGIEDiRG4EBiBA4kRuBAYmxdBEnNbpM0OTnZ2Fq9Xq+xtSRpy5Ytja3F1kXA5xyBA4kROJAYgQOJETiQGIEDiRE4kBiBA4kROJBYpVs29e+o+rGk45KOcedUYDQs5p5s346ID2qbBEBxnKIDiVUNPCT9xvZu25vqHAhAOVVP0VdHxCHbX5O0y/b+iHjl5Af0wyd+oEMqHcEj4lD/v4clPS1p1QKPYesioGOqbD54tu1zT3ws6TuS3qx7MADDq3KK/nVJT9s+8fhfRsTztU4FoIiBgUdET9K3GpgFQGH8mgxIjMCBxAgcSIzAgcQIHEiMwIHECBxIjMCBxBZzPfjn3tatWxtdb2pqqrG1mty6aO3atY2t9cQTTzS2VhdxBAcSI3AgMQIHEiNwIDECBxIjcCAxAgcSI3AgMQIHEqsUuO3zbe+0vd/2jO2r6x4MwPCqvlX1Z5Kej4jv2/6ipLNqnAlAIQMDt32epGsk/UCSIuKopKP1jgWghCqn6MskvS/pYdtv2N7evz86gI6rEvgZkq6U9EBErJD0qaQtpz7I9ibb07anC88I4DRVCXxW0mxEvNb/fKfmg/8vbF0EdM/AwCPiPUkHbS/vf+laSW/VOhWAIqq+in6npB39V9B7km6vbyQApVQKPCL2SuLUGxgxvJMNSIzAgcQIHEiMwIHECBxIjMCBxAgcSIzAgcQIHEiMvckWYW5urtH1JicnG12vKU3uF7Z58+bG1uoijuBAYgQOJEbgQGIEDiRG4EBiBA4kRuBAYgQOJEbgQGIDA7e93Pbek/58ZPvuJoYDMJyBb1WNiLclXSFJtsck/V3S0zXPBaCAxZ6iXyvprxHxtzqGAVDWYi822SDpsYX+wvYmSZuGnghAMZWP4P1ND26UtOClQGxdBHTPYk7Rr5O0JyL+UdcwAMpaTOAb9X9OzwF0U6XAbZ8laZ2kp+odB0BJVfcm+5ekr9Q8C4DCeCcbkBiBA4kROJAYgQOJETiQGIEDiRE4kBiBA4k5Isp/U/t9SYu9pPSrkj4oPkw3ZH1uPK/2fCMiLhz0oFoCPx22p7NeiZb1ufG8uo9TdCAxAgcS61Lg29oeoEZZnxvPq+M68zM4gPK6dAQHUFgnAre93vbbtg/Y3tL2PCXYXmL7JdsztvfZvqvtmUqyPWb7DdvPtj1LSbbPt73T9v7+v93Vbc80jNZP0fv3Wv+L5u8YMyvpdUkbI+KtVgcbku2LJV0cEXtsnytpt6TvjfrzOsH2jyStlHReRNzQ9jyl2H5E0m8jYnv/RqNnRcSHbc91urpwBF8l6UBE9CLiqKTHJd3U8kxDi4h3I2JP/+OPJc1IGm93qjJsT0i6XtL2tmcpyfZ5kq6R9KAkRcTRUY5b6kbg45IOnvT5rJKEcILtpZJWSHqt3UmKuV/SPZI+a3uQwpZJel/Sw/0fP7bbPrvtoYbRhcC9wNfSvLRv+xxJT0q6OyI+anueYdm+QdLhiNjd9iw1OEPSlZIeiIgVkj6VNNKvCXUh8FlJS076fELSoZZmKcr2mZqPe0dEZLkj7WpJN9p+R/M/Tq2x/Wi7IxUzK2k2Ik6cae3UfPAjqwuBvy7pMtuX9F/U2CDpmZZnGppta/5nuZmIuK/teUqJiHsjYiIilmr+3+rFiLil5bGKiIj3JB20vbz/pWsljfSLoovdm6y4iDhm+w5JL0gak/RQROxreawSVku6VdKfbe/tf+0nEfFcizNhsDsl7egfbHqSbm95nqG0/msyAPXpwik6gJoQOJAYgQOJETiQGIEDiRE4kBiBA4kROJDYvwFBuZfzuATFqgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.imshow(imgs[3])\n",
    "plt.gray()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA00AAAEtCAYAAADUaItfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAADIVJREFUeJzt3UmI5GcZx/Hfo41Go3HFGHHDDdGI5KCXqBHEgyiKuBvUiAoaxIuI8RCQoOagggTEg3gQBTdUFLcgQhCXk6JRg4lkkaBxCSHBjMElvh6qJqlkZp7qdFdX/av784GB6e6a+v+re/rt+tbbVU+NMQIAAMDJ3W/TJwAAADBlogkAAKAhmgAAABqiCQAAoCGaAAAAGqIJAACgIZpYiap6clWNqtqZv/39qnrbps8L2F7WFWC/rCOsimjiQIwxXjbG+Pwqr7Oq3lhVV1fVbVX1t6r6fFWdscpjANN1QOvKBVV1Z1XdvvDnxas8BjAdB7SOnF1Vl1fVzVV1wgDUqnpkVX2zqo5V1R+r6s2rPD7rIZo4wfFHYybop0nOHWM8LMlTkuwk+chmTwnYjQmvK0ny8zHGQxb+XLHpEwJONOF15D9JvprkHaf4+KeT/DvJmUnOT/KZqnr2ms6NFRFNJEmq6oaq+mBVXZnkWFXtVNVFVXVtVf2jqq6qqlcvXP7+VfWJ+aMq1yV5+b2u74qqeuf87x+uqi8ufOzeW+UXVNV18+NcX1Xnn+wcxxg3jjFuXnjXnUmetrrPArBK27CuANO2DevIGOPqMcbnkvzuJOd/epLXJLl4jHH7GOMnSb6d5C37/uSwVlMtdjbjTZktLjePMf5bVdcmeWGSvyR5XZIvVtXTxhg3JXlXklckOSfJsSRf38sB54vJZUmeN8a4uqrOSvLI5vIvSPLdJGck+WeSV5/qssAkTH5dSXJOVd2c5JYkX0hy6Rjjv3s5NnAgtmEdOZVnJLlzjHHNwvt+neS8vZwXm2OniUWXzXdz7kiSMcbXxhh/HmP8b4zxlSR/SPL8+WVfn+RT88vfkuTSfRz3f0nOrqoHjTFuGmOc8EjNcWOMn8x/Pe/xST6e5IZ9HBc4eFNfV36c5Owkj8ns0eA3JfnAPo4LrN7U15HOQ5Lcdq/33Zbkofs4LzZANLHoxsU3quqtVfWrqrq1qm7N7I7Fo+cffty9Lv/HvRxwjHEsyRuSvDvJTVX13ap65i7+3Z+S/CDJl/dyXGBtJr2ujDGuG2NcP7/z9ZsklyR57V6OCxyYSa8jS9ye2W/HLDojyT/2cl5sjmhi0V2v+FJVT0ry2STvTfKoMcbDk/w2Sc0vclOSJyz82yc213ssyYMX3n7sPQ46xuVjjJcmOSvJ7+fH3Y2dJE/d5WWBzdi2dWUsnA8wDdu2jiy6JslOVT194X3PzUme/8S0iSZO5fTMFqm/J0lVvT2zR3KO+2qS91XV46vqEUkuaq7rV0leVFVPrKqHJfnQ8Q9U1ZlV9cr57w7/K7NHZO482ZVU1fnz66j5ovnRJD/a+00E1myK68rLqurM+d+fmeTiJN/a6w0EDtwU15GqqtOSPGD+9mlV9cDkrh2rbyS5pKpOr6pzk7wqs+dPskVEEyc1xrgqySeT/DzJX5M8J7OX/D7us0kuz+zJjL/MbEE41XX9MMlXklyZ5BdJvrPw4fsleX+SP2f2JOzzklx4iqt6VpKfZbZw/TTJ1Zk94RPYAhNdV16S5MqqOpbke/Njfuw+3jRgTSa6jjwpyR25e/fojszuoxx3YZIHJflbki8lec8enx/FBtUYJ8zgAgAAYM5OEwAAQEM0AQAANEQTAABAQzQBAAA0RBMAAEBjZ50Hq6qNv1TfslcLrNr8TMMxxuZPArbEFNaVg7abVzmdwtp10KyNrMtRWFeY2Yb7hazHsp8xdpoAAAAaogkAAKAhmgAAABqiCQAAoCGaAAAAGqIJAACgIZoAAAAaa53TtMxuZpEchnPwmv/AolWsO2aNAJxov+urtZXj7DQBAAA0RBMAAEBDNAEAADREEwAAQEM0AQAANEQTAABAQzQBAAA0JjWnCWDdpjAfbh2mcDvNMwFWaQrrmjlOR4edJgAAgIZoAgAAaIgmAACAhmgCAABoiCYAAICGaAIAAGiIJgAAgIZoAgAAaBhuCwDA5ExheO1+7eY2GIC7Hew0AQAANEQTAABAQzQBAAA0RBMAAEBDNAEAADREEwAAQEM0AQAANNY6p+kwvN7+Kvg8AABHnftDM8s+D+Y4TYOdJgAAgIZoAgAAaIgmAACAhmgCAABoiCYAAICGaAIAAGiIJgAAgMZa5zQte535o/J6/T4PAMB+uK9wdByGr/VhmDVlpwkAAKAhmgAAABqiCQAAoCGaAAAAGqIJAACgIZoAAAAaogkAAKCx1jlNAOt2GOZbHBbLvhaHYY4HAIeTnSYAAICGaAIAAGiIJgAAgIZoAgAAaIgmAACAhmgCAABoiCYAAICGOU3AVjOH6fDwtQRgquw0AQAANEQTAABAQzQBAAA0RBMAAEBDNAEAADREEwAAQEM0AQAANEQTAABAY1LDbavqwI+xbHjiOs4BWJ1l37MGpm4PX0sApspOEwAAQEM0AQAANEQTAABAQzQBAAA0RBMAAEBDNAEAADREEwAAQGOtc5rGGJMfgmQOCBwuZv9Mhzl4AGwrO00AAAAN0QQAANAQTQAAAA3RBAAA0BBNAAAADdEEAADQEE0AAAAN0QQAANAQTQAAAA3RBAAA0BBNAAAADdEEAADQEE0AAAAN0QQAANAQTQAAAI2dTZ8AAAD3NMbY9CnAyiz7/1xVazqTvbPTBAAA0BBNAAAADdEEAADQEE0AAAAN0QQAANAQTQAAAA3RBAAA0DCnCQBgzcxhgrttwxwnO00AAAAN0QQAANAQTQAAAA3RBAAA0BBNAAAADdEEAADQEE0AAAAN0QQAANAw3BY40qYwMG8dQy6ncDvhqDC4FlZrN99TB/1zzk4TAABAQzQBAAA0RBMAAEBDNAEAADREEwAAQEM0AQAANEQTAABAw5wmgA1bNltiCvMpgLuZwwTTs+z7cr8/J+00AQAANEQTAABAQzQBAAA0RBMAAEBDNAEAADREEwAAQEM0AQAANMqsAQCAe3DnCI6edpCTnSYAAICGaAIAAGiIJgAAgIZoAgAAaIgmAACAhmgCAABoiCYAAIDGzqZPAABgSqracS0x4xK2z36/r+00AQAANEQTAABAQzQBAAA0RBMAAEBDNAEAADREEwAAQEM0AQAANMxpAgC4D8xxgulZ9n25X3aaAAAAGqIJAACgIZoAAAAaogkAAKAhmgAAABqiCQAAoCGaAAAAGuY0AQCs0G7mxZjlBLt30DOYdsNOEwAAQEM0AQAANEQTAABAQzQBAAA0RBMAAEBDNAEAADREEwAAQEM0AQAANAy3BQBYs2XDOg2/5SiZwvDaZew0AQAANEQTAABAQzQBAAA0RBMAAEBDNAEAADREEwAAQEM0AQAANMxpAgCYGHOcOEy2YQ7TMnaaAAAAGqIJAACgIZoAAAAaogkAAKAhmgAAABqiCQAAoCGaAAAAGqIJAACgIZoAAAAaogkAAKAhmgAAABqiCQAAoCGaAAAAGqIJAACgIZoAAAAaogkAAKAhmgAAABqiCQAAoCGaAAAAGqIJAACgIZoAAAAaogkAAKAhmgAAABqiCQAAoCGaAAAAGqIJAACgIZoAAAAaogkAAKAhmgAAABqiCQAAoCGaAAAAGqIJAACgIZoAAAAaogkAAKAhmgAAABqiCQAAoCGaAAAAGqIJAACgIZoAAAAaogkAAKAhmgAAABqiCQAAoCGaAAAAGqIJAACgIZoAAAAaogkAAKAhmgAAABqiCQAAoLGz6RMAAOC+qapNn8K+jTE2fQpb4TB8rQ8DO00AAAAN0QQAANAQTQAAAA3RBAAA0BBNAAAADdEEAADQEE0AAAANc5oAAFi7ZfOHjsocJ3OYtoOdJgAAgIZoAgAAaIgmAACAhmgCAABoiCYAAICGaAIAAGiIJgAAgIY5TQAATM5hmONkBtPhYacJAACgIZoAAAAaogkAAKAhmgAAABqiCQAAoCGaAAAAGqIJAACgIZoAAAAahtsCALB1djM49qAH4Bpee3TYaQIAAGiIJgAAgIZoAgAAaIgmAACAhmgCAABoiCYAAICGaAIAAGiY0wQAwKG0bI7SsjlO5jBxnJ0mAACAhmgCAABoiCYAAICGaAIAAGiIJgAAgIZoAgAAaIgmAACARi17fXoAAICjzE4TAABAQzQBAAA0RBMAAEBDNAEAADREEwAAQEM0AQAANEQTAABAQzQBAAA0RBMAAEBDNAEAADREEwAAQEM0AQAANEQTAABAQzQBAAA0RBMAAEBDNAEAADREEwAAQEM0AQAANEQTAABAQzQBAAA0RBMAAEBDNAEAADREEwAAQOP/LAwKo/MsMB4AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 864x288 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot3(star(3), star(5), star(10), 'radius 3', 'radius 5', 'radius 10')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA08AAAEtCAYAAADQnVtiAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAADIdJREFUeJzt3V+I/Xldx/HX2w6mra1l0WaYRlmIbYQXdmNlEF1IUUh/bakUE0yimxD1QhCpvKgghOhCupCE0igpMl0ikND2qrDVpDVcFcn1zyIu7q8ldX17Medn429/v5/vc+ac+X7nN48H/OA3M9853/eZmfOZ8zyfOTPV3QEAAODmHrf0AAAAABeBeAIAABgQTwAAAAPiCQAAYEA8AQAADIgnAACAAfHEuaqq76mqrqrN9uV3VtVvLD0XcLFYS4BdWTc4BPHEorr7Bd395kNeZlX9SlXdV1UPVdWnq+rNVXX7Ic8BrMuR1pIXV9WjVfXwqX8/cchzAMs50rpxZ1XdXVUPVtVj/phqVT2lqt5eVVeq6mNV9auHPD/HJ57Y29VHblbovUme191PTvK9STZJfm/ZkYAbWfFakiT3dPeTTv1799IDAateN76Y5G1JXnqDt/9pki8kuSPJXUn+rKp+8Jxm4wDEEzupqo9W1auq6t4kV6pqU1WvrqoPV9Xnq+qDVfXCU8d/Q1X90fYRmPuT/PQ1l/fuqvrN7f9fV1VvOfW2a7fXX1xV92/P85Gquut6M3b3x7v7wVOvejTJMw/3UQDO6iKsJcC6XIR1o7vv6+4/T/Kf15n/tiQ/n+S13f1wd78nyd8n+bUzf3A4N2utdtbtRTlZgB7s7i9V1YeT/FiSTyb5xSRvqapndvcDSV6W5GeSPCfJlSR/s88JtwvOG5M8t7vvq6qnJnnKTY7/0STvSHJ7kv9N8sIbHQssZvVrSZLnVNWDST6b5C+SvKG7v7TPuYGDuAjrxo38QJJHu/tDp173H0mev89cLMPOE/t443Z355Ek6e6/7u5PdPeXu/utSf47yY9sj/2lJH+yPf6zSd5whvN+OcmdVfXE7n6gux/zqM5V3f2e7Y/tPS3JHyb56BnOCxzH2teSf0lyZ5LvyMmjxS9K8soznBc4u7WvGzfzpCQPXfO6h5J88xnm4pyJJ/bx8dMvVNWvV9X7qupzVfW5nNzZ+Pbtm7/rmuM/ts8Ju/tKkl9O8vIkD1TVO6rqWYP3+58k70ryV/ucFziqVa8l3X1/d39ke6fs/Ulen+QX9jkvcDCrXje+jodz8hMxp92e5PP7zMUyxBP7+Opvj6mqZyR5U5LfTvJt3f0tST6QpLaHPJDku0+979NvcrlXknzTqZe/82tO2n13d/9Ukqcm+a/teSc2Sb5veCxwfi7aWtKn5gGWcdHWjdM+lGRTVd9/6nU/nOs8P4r1Ek+c1W05Wcg+kyRV9ZKcPOpz1duS/E5VPa2qvjXJq29yWe9L8uNV9fSqenKS11x9Q1XdUVU/u/254//LyaM3j17vQqrqru1l1HZh/f0k/7z/VQTOwRrXkhdU1R3b/z8ryWuT/N2+VxA4uDWuG1VVT0jy+O3LT6iqb0y+uoP1t0leX1W3VdXzkvxcTp5PyQUhnjiT7v5gkj9Ock+STyX5oZz8qvCr3pTk7pw8IfLfc7Jo3Oiy/inJW5Pcm+TfkvzDqTc/LsnvJvlETp64/fwkr7jBRT07yb/mZHF7b5L7cvKkUWClVrqW/GSSe6vqSpJ/3J7zD3a8asCRrHTdeEaSR/L/u0mP5OR+yFWvSPLEJJ9O8pdJfmvP50+xkOp+zN/vAgAA4Bp2ngAAAAbEEwAAwIB4AgAAGBBPAAAAA+IJAABgYLP0ANdTVav6FYD7/EbCqnX9HcXuXtdAcGRrW0eO7VZYp47NOsg+Lttachntun5etrXzMrrZ9ws7TwAAAAPiCQAAYEA8AQAADIgnAACAAfEEAAAwIJ4AAAAGxBMAAMCAeAIAABgQTwAAAAPiCQAAYEA8AQAADGyWHuAQunvpER7j2DNV1VEvH7i5Na47u1rbdbCuwcW0trXk2C769bXWno2dJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABjZLDwCwBt299AiX3q6fg6o60iRw67C2ca3z+Jq4lddnO08AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAxslh7gerp76RFWz8cIAC4f3/+5CHb9Oq2qI01yeHaeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAY2Cw9AADAZdXdS48Ai9v1dlBVR5rk61tlPO36AbmMC4+PEQAAnC8/tgcAADAgngAAAAbEEwAAwIB4AgAAGBBPAAAAA+IJAABgQDwBAAAMiCcAAIAB8QQAADAgngAAAAY2Sw8AcAzdvfQIHJnPMQDnzc4TAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBgs/QAh1BVR7387t75fY49E3Bzu94G97mdsyyfY9bI1xkc3663s0PeL7fzBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAwGbpAa6nu2vpGc6qu5ceAQAAOCA7TwAAAAPiCQAAYEA8AQAADIgnAACAAfEEAAAwIJ4AAAAGxBMAAMCAeAIAABgQTwAAAAPiCQAAYEA8AQAADIgnAACAAfEEAAAwIJ4AAAAGxBMAAMCAeAIAABgQTwAAAAPiCQAAYEA8AQAADIgnAACAAfEEAAAwIJ4AAAAGxBMAAMDAZukBANagqnY6vruPNMnltevnAADOm50nAACAAfEEAAAwIJ4AAAAGxBMAAMCAeAIAABgQTwAAAAPiCQAAYEA8AQAADIgnAACAAfEEAAAwIJ4AAAAGNksPAHARVdXSI3yN7t75fdZ2HWCN9rltAetyyNuxnScAAIAB8QQAADAgngAAAAbEEwAAwIB4AgAAGBBPAAAAA+IJAABgQDwBAAAMiCcAAIAB8QQAADAgngAAAAaqu5eeAQBglapqpztK7lfB+lTVTsd39w3fwc4TAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMLBZegAAgFtFVe10fHcfaRK4de16OzskO08AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAY2Cw9AADAZVVVOx3f3UeaBJaz6+1gSXaeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAY2Cw9AAAAM1W10/HdfaRJ4MZ2/Tq9SOw8AQAADIgnAACAAfEEAAAwIJ4AAAAGxBMAAMCAeAIAABgQTwAAAAPiCQAAYEA8AQAADIgnAACAAfEEAAAwsFl6AAAAjqOqjn6O7j76OTic8/iauJXZeQIAABgQTwAAAAPiCQAAYEA8AQAADIgnAACAAfEEAAAwIJ4AAAAGxBMAAMCAeAIAABgQTwAAAAPiCQAAYGCz9AAAAFxcVbX0CGfS3Tsdf9GvL2dj5wkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgIHq7qVnAAAAWD07TwAAAAPiCQAAYEA8AQAADIgnAACAAfEEAAAwIJ4AAAAGxBMAAMCAeAIAABgQTwAAAAPiCQAAYEA8AQAADIgnAACAAfEEAAAwIJ4AAAAGxBMAAMCAeAIAABgQTwAAAAPiCQAAYEA8AQAADIgnAACAAfEEAAAwIJ4AAAAGxBMAAMDAVwCiTuN+WYYNJAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 864x288 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot3(disk(3), disk(5), disk(10), 'radius 3', 'radius 5', 'radius 10')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA08AAAEtCAYAAADQnVtiAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAADQRJREFUeJzt3V+I/Xldx/HX2w6lra1m0WaYRllIbYQXdWNlEF1IUUh/bamMClKim4jsQgipvKgghOhCupCE0qgosloikNC8MmwzaQ1XRXJLF3Fxfy2V67uLOT8df/v7/fb9PTNnvt/5zeMBP/jNzHfO9z3nzPnMec5nzkx1dwAAALi9p609AAAAwGUgngAAAAbEEwAAwIB4AgAAGBBPAAAAA+IJAABgQDxxoarqa6qqq2q3f/lvquqn1p4LuFysJcBS1g3Og3hiVd39su5+03leZlX9WFU9WFWPVtXHqupNVXX3eZ4D2JYjrSWvrKonquqxU/++6zzPAaznSOvGvVV1f1U9UlVP+mOqVfWcqvrzqrpWVR+uqh8/z/NzfOKJg13/zs0GvTPJS7r7WUm+Nskuya+vOxJwKxteS5LkXd39zFP/3r72QMCm143/S/LWJD9zi7f/XpL/TXJPkvuS/H5VfdMFzcY5EE8sUlUfqqpfqaoHklyrql1VvaaqPlBVn6qq91XVy08d/wVV9dv778A8lOR7b7i8t1fVz+7//2tV9eZTb7txe/2VVfXQ/jwfrKr7bjZjd3+kux859aonkrzw/K4F4Kwuw1oCbMtlWDe6+8Hu/oMk/3qT+e9K8oNJXtvdj3X3O5L8ZZKfOPOVw4XZarWzba/IyQL0SHd/uqo+kOQ7kvxnkh9O8uaqemF3P5zk55J8X5IXJ7mW5E8POeF+wXlDkm/t7ger6rlJnnOb4789yduS3J3kv5O8/FbHAqvZ/FqS5MVV9UiSTyT5wySv7+5PH3Ju4FxchnXjVr4hyRPd/f5Tr/vnJC89ZC7WYeeJQ7xhv7vzeJJ0959090e7+zPd/ZYk/57k2/bH/kiS390f/4kkrz/DeT+T5N6qekZ3P9zdT/quznXd/Y79j+09L8lvJfnQGc4LHMfW15J/SHJvkq/IyXeLX5Hkl89wXuDstr5u3M4zkzx6w+seTfIlZ5iLCyaeOMRHTr9QVT9ZVe+pqk9W1Sdz8mDjy/dv/qobjv/wISfs7mtJfjTJzyd5uKreVlUvGrzffyT52yR/fMh5gaPa9FrS3Q919wf3D8r+JcnrkvzQIecFzs2m142n8FhOfiLmtLuTfOqQuViHeOIQn/3tMVX1giRvTPILSb6su5+d5L1Jan/Iw0m++tT7Pv82l3styRefevkrP++k3fd39/ckeW6Sf9ufd2KX5OuGxwIX57KtJX1qHmAdl23dOO39SXZV9fWnXvctucnzo9gu8cRZ3ZWThezjSVJVP52T7/pc99Ykv1hVz6uqL03ymttc1nuSfGdVPb+qnpXkV6+/oaruqarv3//c8f/k5Ls3T9zsQqrqvv1l1H5h/Y0kf3/4hwhcgC2uJS+rqnv2/39Rktcm+YtDP0Dg3G1x3aiqenqSL9y//PSq+qLksztYf5bkdVV1V1W9JMkP5OT5lFwS4okz6e73JfmdJO9K8l9Jvjknvyr8ujcmuT8nT4j8p5wsGre6rL9L8pYkDyR5d5K/OvXmpyX5pSQfzckTt1+a5NW3uKhvTPKPOVnc3pnkwZw8aRTYqI2uJd+d5IGqupbkr/fn/M2FHxpwJBtdN16Q5PF8bjfp8Zw8Drnu1UmekeRjSf4oyasOfP4UK6nuJ/39LgAAAG5g5wkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABnZrD3AzVbWpXwF4yG8krNrW31Hs7m0NBEe2tXXk2O6EderYrIMc4qqtJVfR0vXzqq2dV9Htvl7YeQIAABgQTwAAAAPiCQAAYEA8AQAADIgnAACAAfEEAAAwIJ4AAAAGxBMAAMCAeAIAABgQTwAAAAPiCQAAYGC39gBr6O7NnaOqjjQJcBlZpwAOc+z109p5tdl5AgAAGBBPAAAAA+IJAABgQDwBAAAMiCcAAIAB8QQAADAgngAAAAbEEwAAwIB4AgAAGBBPAAAAA+IJAABgYLf2AOehu9ce4cyWfgxVdaRJgGOwTgEc5rKvn9bOO4udJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABnZrD3Az3b32CJvnOoJ1uQ8+taXXUVUdaRJgS6yft2ft3DY7TwAAAAPiCQAAYEA8AQAADIgnAACAAfEEAAAwIJ4AAAAGxBMAAMCAeAIAABgQTwAAAAPiCQAAYEA8AQAADOzWHgBgC7p77RGuvKW3QVUdaRJgCevnuqydF2uT8bT0Rr2Kd1rXEQAAXCw/tgcAADAgngAAAAbEEwAAwIB4AgAAGBBPAAAAA+IJAABgQDwBAAAMiCcAAIAB8QQAADAgngAAAAZ2aw8AcAzdvfYIHJnbGM6f+9Wdb+ltXFVHmuRyuiPiaemNusWFwScmAABsmx/bAwAAGBBPAAAAA+IJAABgQDwBAAAMiCcAAIAB8QQAADAgngAAAAbEEwAAwIB4AgAAGBBPAAAAA+IJAABgYLf2AGuoqkXHd/fRzwGcr4u4n7MutzE8NZ/3nJXHwZ/PzhMAAMCAeAIAABgQTwAAAAPiCQAAYEA8AQAADIgnAACAAfEEAAAwIJ4AAAAGxBMAAMCAeAIAABgQTwAAAAO7tQe4me6utWc4q+5eewRggaply477+PlbehvAVWTt4TJY+nl6mdZ/O08AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAzs1h4A4DKqqkXHd/eRJtmupdcRXEVXcW2AGy29H6z59cXOEwAAwIB4AgAAGBBPAAAAA+IJAABgQDwBAAAMiCcAAIAB8QQAADAgngAAAAbEEwAAwIB4AgAAGBBPAAAAA7u1BwC4Cqpq0fHdfaRJDrf0Y4CraIv3XbjTLL2fnefXLztPAAAAA+IJAABgQDwBAAAMiCcAAIAB8QQAADAgngAAAAbEEwAAwIB4AgAAGBBPAAAAA+IJAABgQDwBAAAM7NYeAIAnq6pFx3f30c8BV9Eh9y1gW87zfmznCQAAYEA8AQAADIgnAACAAfEEAAAwIJ4AAAAGxBMAAMCAeAIAABgQTwAAAAPiCQAAYEA8AQAADIgnAACAgerutWcAANikqlr0QMnjKtieqlp0fHff8h3sPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMLBbewAAgDtFVS06vruPNAncuZbez86TnScAAIAB8QQAADAgngAAAAbEEwAAwIB4AgAAGBBPAAAAA+IJAABgQDwBAAAMiCcAAIAB8QQAADAgngAAAAZ2aw8AAHBVVdWi47v7SJPAepbeD9Zk5wkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgIHd2gMAADBTVYuO7+4jTQK3tvTz9DKx8wQAADAgngAAAAbEEwAAwIB4AgAAGBBPAAAAA+IJAABgQDwBAAAMiCcAAIAB8QQAADAgngAAAAbEEwAAwMBu7QEAADiOqlp0fHcfaRIuq6WfQ3c6O08AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgIHd2gMAALANVbX4fbr7CJNwLIfcxnyOnScAAIAB8QQAADAgngAAAAbEEwAAwIB4AgAAGBBPAAAAA+IJAABgQDwBAAAMiCcAAIAB8QQAADAgngAAAAZ2aw8AAMDlVVWLju/uI01yNS29/jkbO08AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAzs1h4AAICro6oWHd/dR5pkm5ZeP1wsO08AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAzs1h4AAABupaoWHd/dR5rkMEvnZ9vsPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMLBbewAAADgvVbXo+O4+6uVzZ7HzBAAAMCCeAAAABsQTAADAgHgCAAAYEE8AAAAD4gkAAGBAPAEAAAyIJwAAgAHxBAAAMCCeAAAABsQTAADAQHX32jMAAABsnp0nAACAAfEEAAAwIJ4AAAAGxBMAAMCAeAIAABgQTwAAAAPiCQAAYEA8AQAADIgnAACAAfEEAAAwIJ4AAAAGxBMAAMCAeAIAABgQTwAAAAPiCQAAYEA8AQAADIgnAACAAfEEAAAwIJ4AAAAGxBMAAMCAeAIAABgQTwAAAAPiCQAAYOD/AUWHIpkS4oC4AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 864x288 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot3(diamond(3), diamond(5), diamond(10), 'radius 3', 'radius 5', 'radius 10')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA1gAAACiCAYAAABcfmxKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAABlVJREFUeJzt3THI7XUdx/HP13shpGi41lW4lYNBQltTDsbVxUEocqpBbEjEJRyamoKWiihokJaGxFWkwWgyXaolIqyUCFriyh00UUtE89dwr3ApH33O7XOO9//4esEZnsPvefhOX877+Z//ObPWCgAAAP+/697vAQAAAE4KgQUAAFAisAAAAEoEFgAAQInAAgAAKBFYAAAAJQILAACgRGBRMzOPzszzM/PyzPxlZr7+LmfvmJlnZualmXlhZh6fmXOHnBfYjh33y/mZeWtmXr3icd8h5wW2Ycfd8q3/2iuvXd41HzvkzFz7xhcN0zIzn03y17XW6zNza5Knkty91vrdO5y9McmptdaFmflQku8kuXWt9cWDDg1swo775XySR9danzjslMDW7LJb3uF3v53kC2utO/c7JVvjChY1a60/rbVef/vHy49bjjh7ca114Yqn/p3k03seEdioXfYLwHFd7W6ZmUlyb5Kf7XE8NkpgUTUzD8/Mv5I8l+T5JL94l7OfmpmXkryW5JtJvn+YKYEt2mW/JDk7Mxdn5m8z86OZ+fBhpgS2Zsfd8rbbk9yY5LF9zsY2eYsgdTNzKsltSc4n+d5a6433OH8myf1Jnl5r/Xb/EwJbdZz9MjM3JTmTSy+Wbs6l/zA/u9Z64ICjAhtyFa9dfppLtzp8bf/TsTUCi72ZmZ8k+fNa68fHOHtTkj8kObfWenPvwwGbtuN++XySJ9ZaN+x/MmDLjrNbZub6JBeTfGmt9auDDcdmeIsg+3Q6x79H4nSSs0k+ur9xgBNkl/2yksweZwFOjuPslnuSvJhLH4gB/0NgUTEzZ2fmKzPzkZk5NTN3JflqkiePOH/PzHxmZq6bmY8n+WGS36+1Xjzk3MC17yr2y/nL93jOzHwyyXeT/PyQMwPXvl13yxXuS/LI8jYwjiCwaFlJHkzy9yT/SPKDJA+ttY56UXMuyS+TvJLkmSRvJfnyAeYEtmfX/fK5JL9J8s8kv07yxyTfOMCcwLbsulty+Ts770zyyEEmZJPcgwUAAFDiChYAAECJwAIAACgRWAAAACUCCwAAoOT0Pv7ozPjkjBNkreX7Y7hm2C980NnJ+2G3ALs6ah+7ggUAAFAisAAAAEoEFgAAQInAAgAAKBFYAAAAJQILAACgRGABAACUCCwAAIASgQUAAFAisAAAAEoEFgAAQInAAgAAKBFYAAAAJQILAACgRGABAACUCCwAAIASgQUAAFAisAAAAEoEFgAAQInAAgAAKBFYAAAAJQILAACgRGABAACUCCwAAIASgQUAAFAisAAAAEoEFgAAQInAAgAAKBFYAAAAJQILAACgRGABAACUCCwAAIASgQUAAFAisAAAAEoEFgAAQInAAgAAKBFYAAAAJQILAACgRGABAACUCCwAAIASgQUAAFAisAAAAEoEFgAAQInAAgAAKBFYAAAAJQILAACgRGABAACUCCwAAIASgQUAAFAisAAAAEoEFgAAQInAAgAAKBFYAAAAJQILAACgRGABAACUCCwAAIASgQUAAFAisAAAAEoEFgAAQInAAgAAKBFYAAAAJQILAACgRGABAACUCCwAAIASgQUAAFAisAAAAEoEFgAAQInAAgAAKBFYAAAAJQILAACgRGABAACUCCwAAIASgQUAAFAisAAAAEoEFgAAQInAAgAAKBFYAAAAJQILAACgRGABAACUCCwAAIASgQUAAFAisAAAAEoEFgAAQInAAgAAKBFYAAAAJQILAACgRGABAACUCCwAAIASgQUAAFAisAAAAEoEFgAAQInAAgAAKBFYAAAAJQILAACgRGABAACUCCwAAIASgQUAAFAisAAAAEoEFgAAQInAAgAAKBFYAAAAJQILAACgRGABAACUCCwAAIASgQUAAFAisAAAAEoEFgAAQInAAgAAKBFYAAAAJQILAACgRGABAACUCCwAAIASgQUAAFAisAAAAEoEFgAAQMmstd7vGQAAAE4EV7AAAABKBBYAAECJwAIAACgRWAAAACUCCwAAoERgAQAAlAgsAACAEoEFAABQIrAAAABKBBYAAECJwAIAACgRWAAAACUCCwAAoERgAQAAlAgsAACAEoEFAABQIrAAAABKBBYAAECJwAIAACgRWAAAACUCCwAAoERgAQAAlAgsAACAkv8A10sSPuPQYgYAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 864x288 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot3(rectangle(3, 3), rectangle(3, 5), rectangle(3,7), '3 3', '3 5', '3 7')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPgAAAD8CAYAAABaQGkdAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAACstJREFUeJzt3V+IXOUZx/Hfr6vS+g9Da4vshsYVCUihxoSABITGtMQq2osaElCoFNYbRWlBY+9655XYiyKEqBVMlW5UELHaBBUrtNbdJG2NG0u6WLKJNoqRqIWGxKcXO4E0XTtnM+e858zj9wPB/TPs+0zWb87Z2ZnzOiIEIKcvtT0AgOYQOJAYgQOJETiQGIEDiRE4kBiBA4kROJAYgQOJndXEF7Wd8ulxS5YsKbre6OhosbWOHj1abK2DBw8WW+vEiRPF1iotItzvNo0EntW6deuKrnf//fcXW2vnzp3F1tq8eXOxtY4cOVJsrS7iFB1IjMCBxAgcSIzAgcQIHEiMwIHECBxIjMCBxCoFbnu97bdt77dd7lkKAAbSN3DbI5J+Kek6SVdI2mT7iqYHAzC4Kkfw1ZL2R8RsRByT9KSkm5odC0AdqgQ+KunAKe/P9T4GoOOqvNhkoVes/M+rxWxPSJoYeCIAtakS+Jykpae8Pybp0Ok3iogtkrZIeV8uCgybKqfob0i63Palts+RtFHSs82OBaAOfY/gEXHc9h2SXpQ0IumRiNjb+GQABlbpgg8R8byk5xueBUDNeCYbkBiBA4kROJAYgQOJETiQGIEDiRE4kBiBA4mxs8kilNxpRJLGx8eLrVVyW6YPP/yw2FobNmwotpYkTU5OFl2vH47gQGIEDiRG4EBiBA4kRuBAYgQOJEbgQGIEDiRG4EBiVXY2ecT2YdtvlhgIQH2qHMF/JWl9w3MAaEDfwCPiVUnlnjwMoDb8DA4kVturydi6COie2gJn6yKgezhFBxKr8muyJyT9QdJy23O2f9z8WADqUGVvsk0lBgFQP07RgcQIHEiMwIHECBxIjMCBxAgcSIzAgcQIHEhs6LcuWrlyZbG1Sm4lJEmXXXZZsbVmZ2eLrbVjx45ia5X8/0Ni6yIABRE4kBiBA4kROJAYgQOJETiQGIEDiRE4kBiBA4kROJBYlYsuLrX9su0Z23tt31ViMACDq/Jc9OOSfhoRu2xfIGna9o6IeKvh2QAMqMreZO9GxK7e2x9LmpE02vRgAAa3qFeT2V4maYWk1xf4HFsXAR1TOXDb50t6StLdEXH09M+zdRHQPZUeRbd9tubj3hYRTzc7EoC6VHkU3ZIeljQTEQ80PxKAulQ5gq+RdKuktbb39P58v+G5ANSgyt5kr0lygVkA1IxnsgGJETiQGIEDiRE4kBiBA4kROJAYgQOJETiQ2NDvTbZkyZJia01PTxdbSyq7X1hJpf8ev8g4ggOJETiQGIEDiRE4kBiBA4kROJAYgQOJETiQGIEDiVW56OKXbf/J9p97Wxf9vMRgAAZX5amq/5a0NiI+6V0++TXbv42IPzY8G4ABVbnoYkj6pPfu2b0/bGwADIGqGx+M2N4j6bCkHRGx4NZFtqdsT9U9JIAzUynwiDgREVdKGpO02va3FrjNlohYFRGr6h4SwJlZ1KPoEfGRpFckrW9kGgC1qvIo+sW2L+q9/RVJ6yTta3owAIOr8ij6JZIesz2i+X8QfhMRzzU7FoA6VHkU/S+a3xMcwJDhmWxAYgQOJEbgQGIEDiRG4EBiBA4kRuBAYgQOJMbWRYuwc+fOYmtlVvJ7duTIkWJrdRFHcCAxAgcSI3AgMQIHEiNwIDECBxIjcCAxAgcSI3AgscqB966Nvts212MDhsRijuB3SZppahAA9au6s8mYpOslbW12HAB1qnoEf1DSPZI+a3AWADWrsvHBDZIOR8R0n9uxNxnQMVWO4Gsk3Wj7HUlPSlpr+/HTb8TeZED39A08Iu6LiLGIWCZpo6SXIuKWxicDMDB+Dw4ktqgrukTEK5rfXRTAEOAIDiRG4EBiBA4kRuBAYgQOJEbgQGIEDiRG4EBiQ791UcmtaVauXFlsrdJKbidU8u9xcnKy2FpdxBEcSIzAgcQIHEiMwIHECBxIjMCBxAgcSIzAgcQIHEis0jPZeldU/VjSCUnHuXIqMBwW81TV70TEB41NAqB2nKIDiVUNPCT9zva07YkmBwJQn6qn6Gsi4pDtr0vaYXtfRLx66g164RM/0CGVjuARcaj338OSnpG0eoHbsHUR0DFVNh88z/YFJ9+W9D1JbzY9GIDBVTlF/4akZ2yfvP2vI+KFRqcCUIu+gUfErKRvF5gFQM34NRmQGIEDiRE4kBiBA4kROJAYgQOJETiQGIEDiTki6v+idv1f9HOMj4+XWkpTU1PF1pKk22+/vdhaN998c7G1Sn7PVq3K+9KIiHC/23AEBxIjcCAxAgcSI3AgMQIHEiNwIDECBxIjcCAxAgcSqxS47Ytsb7e9z/aM7aubHgzA4KpeF/0Xkl6IiB/aPkfSuQ3OBKAmfQO3faGkayT9SJIi4pikY82OBaAOVU7RxyW9L+lR27ttb+1dHx1Ax1UJ/CxJV0l6KCJWSPpU0ubTb2R7wvaU7bIvuQLwuaoEPidpLiJe772/XfPB/xe2LgK6p2/gEfGepAO2l/c+dK2ktxqdCkAtqj6Kfqekbb1H0Gcl3dbcSADqUinwiNgjiVNvYMjwTDYgMQIHEiNwIDECBxIjcCAxAgcSI3AgMQIHEiNwILGh35uspImJiaLr3XvvvcXWmp6eLrbWhg0biq2VGXuTAV9wBA4kRuBAYgQOJEbgQGIEDiRG4EBiBA4kRuBAYn0Dt73c9p5T/hy1fXeJ4QAMpu9FFyPibUlXSpLtEUkHJT3T8FwAarDYU/RrJf09Iv7RxDAA6lX1uugnbZT0xEKfsD0hqeyrMQD8X5WP4L1ND26UNLnQ59m6COiexZyiXydpV0T8s6lhANRrMYFv0uecngPopkqB2z5X0nclPd3sOADqVHVvsn9J+mrDswCoGc9kAxIjcCAxAgcSI3AgMQIHEiNwIDECBxIjcCCxprYuel/SYl9S+jVJH9Q+TDdkvW/cr/Z8MyIu7nejRgI/E7ansr4SLet94351H6foQGIEDiTWpcC3tD1Ag7LeN+5Xx3XmZ3AA9evSERxAzToRuO31tt+2vd/25rbnqYPtpbZftj1je6/tu9qeqU62R2zvtv1c27PUyfZFtrfb3tf73l3d9kyDaP0UvXet9b9p/ooxc5LekLQpIt5qdbAB2b5E0iURscv2BZKmJf1g2O/XSbZ/ImmVpAsj4oa256mL7cck/T4itvYuNHpuRHzU9lxnqgtH8NWS9kfEbEQck/SkpJtanmlgEfFuROzqvf2xpBlJo+1OVQ/bY5Kul7S17VnqZPtCSddIeliSIuLYMMctdSPwUUkHTnl/TklCOMn2MkkrJL3e7iS1eVDSPZI+a3uQmo1Lel/So70fP7baPq/toQbRhcC9wMfSPLRv+3xJT0m6OyKOtj3PoGzfIOlwREy3PUsDzpJ0laSHImKFpE8lDfVjQl0IfE7S0lPeH5N0qKVZamX7bM3HvS0islyRdo2kG22/o/kfp9bafrzdkWozJ2kuIk6eaW3XfPBDqwuBvyHpctuX9h7U2Cjp2ZZnGphta/5nuZmIeKDteeoSEfdFxFhELNP89+qliLil5bFqERHvSTpge3nvQ9dKGuoHRRe7N1ntIuK47TskvShpRNIjEbG35bHqsEbSrZL+antP72M/i4jnW5wJ/d0paVvvYDMr6baW5xlI678mA9CcLpyiA2gIgQOJETiQGIEDiRE4kBiBA4kROJAYgQOJ/Qcpuo92pLZ1pQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "img1 = imgs[0]\n",
    "img2 = imgs[1]\n",
    "img3 = imgs[2]\n",
    "plt.imshow(img1)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## erosion/dilation"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### disk(1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0wAAAEtCAYAAAA7quBhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAFtdJREFUeJzt3XvMbXlZH/Dvg4MMyk2Klw43KaREqdTaUm2DrUhTa1Pa/tE2iheI0qZN29QG1FSpzlgUayoSNTEmrVpABFN7UUovGKN0ICZIIzbaWJEOjAw64DByp87w6x9rn/DOOfv3nr3PWnuttff+fJKTvOd9917rt9Ze61nr2b/9PLtaawEAAOBaD1p6AAAAAGslYQIAAOiQMAEAAHRImAAAADokTAAAAB0SJgAAgA4JEwAPUFV3VNVf6vztS6vqN+ce05Kq6r9U1fOWHgccs2ONK1V1a1W9avPzE6rqQ1X1KZv//2JVveAGl/uAZbFuNy09AACOR2vtfyR56tLjmFNr7SuXHgOcsmOJK621dyV52I08t6ruSPKC1trPj10W8zPDxGK8qwJcUYODXZOqyhuEcGYOHVc4Hw4iuqrqW6vq3VX1war6zap6dlU9pKpeXlV3bf69vKoesnn886vq9quW0arqKZuff6KqfqSqXl9VH07yrKp6aFV9f1W9s6r+oKpur6qHbh7/JVX15qq6t6reVlVfNvc+gDP2jKr6jap6f1X9eFXdnCRV9WVV9TtXHrT5mM2LqurXNufway889jOq6nVV9d7Ncl5XVY+78NxfrKrvrqo3JflIkhdW1VsvDqKqXlhV/3HbAKvqkVX1b6rqPZtY9ZILH5V5flW9qap+oKruSXJrVT2oql68iTd3V9UrquqRm8ffXFWvqqrf38Sct1TVZ18Y5ws2P1+2jM/dxLznVdW7qup9VfXtk70icPyOIa48qap+aXPv84Ykj7nwtyvn+DVvwFTVk6vqFzYx5H1V9ZNV9ajN316Z5AlJfq6Gj+F9y9XLqqpbqupnq+qeqnp7Vf3dC8u+tap+ehNvPlhVv15Vf+YG9j83SMLEVlX11CT/KMkzWmsPT/IVSe5I8u1JviTJFyb5k0n+bJIX77Ho5yb57iQPT3J7kn+V5E8n+fNJHp3kW5J8oqoem+Q/J3nJ5vcvSvIzVfWZY7cN2MnXZDjvn5zkj+fy8/zvJPkrSZ6U5OlJnr/5/YOS/HiSJ2a4Wfhokh++6rlfl+TvZYgJP5jkSVX1eRf+/rVJXtlZ779Ncl+SpyT5U0n+cpKL9QRfnOQdST4rQ9x5/ubfs5L8sQwfh7kynucleWSSxyf5I0n+/ma8V7tsGVc8M8PHi56d5Duu2h44Z8cQV16d5K0ZEqV/kSE27KKSvDTJLUk+L0MsuTVJWmtfl+RdSZ7TWntYa+37tjz/p5L8zub5fyvJ91TVsy/8/a8neU2SRyX52Vy7zRyQhIme+5M8JMnnV9WDW2t3tNZ+O0Ow+67W2t2ttfcmuS1DYNrVf2qtvam19okk/y/JNyT5J621d7fW7m+tvbm19vEMwez1rbXXt9Y+0Vp7Q5JfSfJXp9xIoOuHW2t3ttbuyZBsfPUlj/3B1tpdm8f+XIY3VNJa+/3W2s+01j7SWvvgZjl/8arn/kRr7ddba/dtzv3XZjj/U1VPS/K5SV539Qo3sz9fmeSbWmsfbq3dneQHknzVhYfd1Vr7oc2yP5ohfr2stfaO1tqHkvyzJF+1eYf3DzMkSk/ZxKK3ttY+sGVbL1vGFbe11j7aWntbkrdleHMJWH9ceUKSZyT55621j7fW3rhZ93W11t7eWnvD5nnvTfKyLePaqqoen+GNlm9trX2stfarSf51Hnh/dfvmnuj+DMmeuDIjCRNbtdbenuSbMrw7cndVvaaqbsnwzsc7Lzz0nZvf7erOCz8/JsnNSX57y+OemORvbz4ac29V3ZshmPzRPdYF3LiL5+r1zvPfvfDzR7IpZK6qT6uqH918fO0DSd6Y5FH1wPrFi+tJhlmj51ZVZbhZ+OnNDc/VnpjkwUnecyFG/GiG2aTesrfFr5uSfHaGG5D/luQ1NXzc+Puq6sFb1nvZMq7Yuj+A1ceVW5K8v7X24avGeV1V9Vmbe6V3b8b1qlz4ON913JLknk0CeHG9j73w/6v3x83bPhrIYUiY6Gqtvbq19swMNyYtyb9Mctfm/1c8YfO7JPlwkk+78oeq+pxti73w8/uSfCzD1PzV7kzyytbaoy78+/TW2vfe8AYB+3j8hZ8vnuf7eGGGj6Z9cWvtEUn+wub3deExF2NCWmu/nGH2+UszfIS397GZO5N8PMljLsSIR7TWntZbdrbHr/uS/F5r7Q9ba7e11j4/w0eE/1qSr9+y3u4yOuMEPmntceU9ST6jqj79qnHu4qWb9T59M66vvWxMV7kryaOr6uFXrffdO66bA5MwsVVVPbWqvryGhg4fy/AZ4fszfMb2xVX1mVX1mCTfkeFdlGT46MnTquoLN8WZt162js3H8n4sycs2xY6fUlV/brPOVyV5TlV9xeb3N2+KQh932TKByfzDqnpcVT06ybdl+EjLvh6eIXbcu1nOd+74vFdk+Hz+fa2127c9oLX2niT/Pcn3V9UjamjG8OSquuwjMD+V5J9uirofluR7kry2tXZfVT2rqr5g8y71BzJ8RO/+fZax47bBOVt7XHlnho//31ZVn1pVz0zynD3G9aHNuB6b5Juv+vvvZah73LbeO5O8OclLN/c7T0/yjUl+csd1c2ASJnoekuR7M8wC/W6Gj7l8W4YmDL+S5NeS/K8k/3Pzu7TW/k+S70ry80l+K0NTh+t50WY5b0lyT4ZZrAdtgsff2KzzvRneTf7mOGZhLq/OkJC8Y/PvJTewjJcneWiGOPLLSf7rjs97ZZI/kf67wFd8fZJPTfIbSd6f5N/l8o/t/thmmW9M8n8zvBn0jzd/+5zN8z+Q5H8n+aV88s2gXZcBXO4Y4spzMzSMuSdDMvaKHZd/W5IvSvIHGZpW/fur/v7SDG8431tVL9ry/K/OUFt1V5L/kOQ7N/XbrEC1dtkMIQDMq4avFrg7yRe11n5r6fEAx09cYQzv1gOwNv8gyVvc1AATEle4YbprALAaVXVHhkLpv7nwUIATIa4wlo/kAQAAdPhIHgAAQIeECQAAoGPuGiaf/zsDwxdpw6C1dugD4pq4ckrH4K4fmx6zzdvWcUr7kNOzRFzhPGyLfadUvrJrbB+zzWvbh3tsc/eBZpgAAAA6JEwAAAAdEiYAAIAOCRMAAEDH4l9ce0qFxedYnD2miE8R5fh1HEMR5aGtZRyHMsf2nfo+BNjVKd2bbDPH9q1tH05x72yGCQAAoEPCBAAA0CFhAgAA6JAwAQAAdCze9GFthWG7WqpI+pSKs491W5Y6Zo/1XDm0tTVCgTUY04ToELHGOXk6xjQ5OsRxMPXxurYGS6yDGSYAAIAOCRMAAECHhAkAAKBDwgQAANCxeNOHJQtBFYv3zbFvliyiVNTZd2z7xjkL1xpzXjinuMyY68GaryVXHMMYmZ8ZJgAAgA4JEwAAQIeECQAAoEPCBAAA0FEzF7eppDsDCoa5qLV26APimrhySs1kdo3Rx3DeabTDVJaIK5yHtTU+mno8u8bcY2h+seu+2WObuw80wwQAANAhYQIAAOiQMAEAAHRImAAAADpuWnoAAFObuiD2EMWvUy9z10LXuQp5xxQWawQBrMnUjXsOEeOmXubUDRXGWroBkhkmAACADgkTAABAh4QJAACgQ8IEAADQUTN/k++kKztEQfNSRdJrK86e41ukxzhE4eFSxYxrK6Kcet2XfXP2VKu40SdOXch7qGXu4gCv22LrhutZc1zhuK0thi91XZj6HnSfsSy17sviihkmAACADgkTAABAh4QJAACgQ8IEAADQIWECAADouGnpAexq105u+3R8W6qz05iudGM72u26zXt0FNl53TdqzOt8iK4zY4zpSje2o92u2zzmWFyzqce75PbPse5je31ZjyU7f3LaluwofOj1LrUdc617bdu3b0wywwQAANAhYQIAAOiQMAEAAHRImAAAADqOpunD1AXrvccqTB2sad9M3aii99glCxLX5FT3zdTNQ9bWTAbWwjnAoazpWjT19WNtTap4IDNMAAAAHRImAACADgkTAABAh4QJAACgo2YuJlO5toe5iu93LdDddd0Kfsebq+nGHK9pa+3QB8Tq48oc5/LUx8c+43POMzdxhUM5hng2xz3CkteopdZ9WVwxwwQAANAhYQIAAOiQMAEAAHRImAAAADo0fWByx1AwyXwUZ6/LMTSTEUO4HnGFQxF/9nMuTarMMAEAAHRImAAAADokTAAAAB0SJgAAgI6blh4AAPNZssEDAKdlyQYPczLDBAAA0CFhAgAA6JAwAQAAdEiYAAAAOjR9ODMjvwF5wpEwtzGvn6L+a43dJ0udT7uOe65vb3dsAYy/JiwVS3cd97bxHeI6eKhrlxkmAACADgkTAABAh4QJAACgQ8IEAADQcRZNH86xOPsQ5irYm9I5FlEewlwNAKYy9bkzV7ODqc+nMeNe8+vLuq0tfnE6lrovmvqYPqZmB9dbxxzPHWuKfWOGCQAAoEPCBAAA0CFhAgAA6JAwAQAAdCze9OHUC+R2tbYGD3Os+9QL7nd1DAXSay7OnrqxwdobmfTMUdy75uOA9XCccChzxOepGxsc6/kwRxOiJZuW7csMEwAAQIeECQAAoEPCBAAA0CFhAgAA6Fi86cParL3g+xBFc7sWKZ5Sgfw2ay/MPEQx6a6v6SkVsq7F2s8dry/A8Vh7zF5b07J995cZJgAAgA4JEwAAQIeECQAAoEPCBAAA0LHKpg9LFoYt1dhgrm3edT1LvQZLFi0u1dhgrm3edT1rLxyd0inFmqmP333Gsuu2nNOxBZyfU7qHmfoatc9Ydt2WXcczxetihgkAAKBDwgQAANAhYQIAAOiQMAEAAHSssunDHMXQrNMcRYtwxa6x5lgbGBxDU4tj2I8Au9g1xi3VYGysY2hqcaj9aIYJAACgQ8IEAADQIWECAADokDABAAB0LN70YY4i4N7yjqHAbmq77ts5iuvmKNbrLe8cC83HfCP2mov1l4whazfmfB+zvEM9H+DQlrw3WbuprxVj98Ouz5/iPsEMEwAAQIeECQAAoEPCBAAA0CFhAgAA6JAwAQAAdCzeJU/XpMHaOorMYU1jWZIOZeOc6nbta8x+2LUrlK6TwKlzbzIYsx927Up3TF2szTABAAB0SJgAAAA6JEwAAAAdEiYAAICOxZs+0Ldr0dw+xXFjn888dn2d9im2H/t8zotjA4CpHPu9phkmAACADgkTAABAh4QJAACgQ8IEAADQcRZNH3w7/eVOuRHEMX2L9BI0gtiN/QTAVE75vutUmWECAADokDABAAB0SJgAAAA6JEwAAAAdNXORmYq2M6AYnotaa4c+IK6JK7seg4eIf0uu+2q7FhaPLUAes83iBTdiibjCeRgTDw8Rz5Zc99XGXD/2Gd+YbR6z7sviihkmAACADgkTAABAh4QJAACgQ8IEAADQoekDk1PEzUWKs9mXGML1iCscivhzvjR9AAAAuAESJgAAgA4JEwAAQIeECQAAoEPCBAAA0CFhAgAA6JAwAQAAdEiYAAAAOiRMAAAAHRImAACADgkTAABAh4QJAACgQ8IEAADQIWECAADokDABAAB0SJgAAAA6JEwAAAAdEiYAAIAOCRMAAEDHTUsPAGBqVbX0ECbTWtvpcae0zQBrsmscPga7XivOcZsvY4YJAACgQ8IEAADQIWECAADokDABAAB0LN704ZQKlc+xOHtMUeA5FhSe4zYf2lrGcSinvn0Aa3JK1+ltTn37ttm2zfteW80wAQAAdEiYAAAAOiRMAAAAHRImAACAjsWbPmxzDAVpCrEP4xj26zEcn8B5O4Y4dQzxnmsdw+t2DMc/x8UMEwAAQIeECQAAoEPCBAAA0CFhAgAA6Fhl04djKCjkMBRqAoznOsqhuE5zjswwAQAAdEiYAAAAOiRMAAAAHRImAACAjsWbPpxj8eA5bvM251iUfI7bfGjbzqdT2s+7xotT2maApWyLpad037brteIct/kyZpgAAAA6JEwAAAAdEiYAAIAOCRMAAECHhAkAAKBj8S55AFMb093nEB2S5ui6NKZb4CG6IY1Zt45/wJqMiUmHiHFzxM0x161DxPClO8aaYQIAAOiQMAEAAHRImAAAADokTAAAAB11iGLfS0y6smMtzt51vdusrTh7zPLGONYiyl3Xu82xFlG21g69E2/4pJi6kPdQy9zF1MfHPmPRpIG5rTmucNzWFsOXui5Mfb+5z1iWWvdlccUMEwAAQIeECQAAoEPCBAAA0CFhAgAA6Diapg+7NmPYp2nDmgr7xhTfL1mcPUejhDGv8zEUUY5punHsRZQT2WnDjrUxwVLNSJjGOb5+595MZi5LNala0hwNspZswjWHczxutunsB00fAAAA9iVhAgAA6JAwAQAAdEiYAAAAOm5aegC7GtvgYa5lMr2xDR7mWibrNvXru7ZmMqyTGMKhiBfLmvp+cW1NqnggM0wAAAAdEiYAAIAOCRMAAECHhAkAAKCjZi4aVKHYseQ3S0+9bkWG+1lyv86x7su+OXsik54Uh2j8MkczmamPj33G55xnbscWVzgec8TSseuYo0nVkteopdZ9WVwxwwQAANAhYQIAAOiQMAEAAHRImAAAADo0fWByCsC5SHH2co61mYwYwvWIKxyK+NN3zk2qzDABAAB0SJgAAAA6JEwAAAAdEiYAAICOm5YeAACHMVdTnzFFtgqsAY7DXPF6zLXrUI2EzDABAAB0SJgAAAA6JEwAAAAdEiYAAIAOTR/OzNQFe3MVlTPe1K/VuRfrj93+Oc6dkd94PunyAOgbe02YIz6PGeO28c11DznF9cwMEwAAQIeECQAAoEPCBAAA0CFhAgAA6DiLpg/nWJzN4ByLKNl9v+y67w+xn6cugJ16jI4tbtSS5xWnberYvqupj+lD3LdN3ahnqX24VmaYAAAAOiRMAAAAHRImAACADgkTAABAx+JNH06pQI79nFLBPcsZc3w4DjgnrnkcyinF0jHjdo6dLjNMAAAAHRImAACADgkTAABAh4QJAACgY/GmDwDA4c1VfK7I/fzM1eDhWBtJcPzMMAEAAHRImAAAADokTAAAAB0SJgAAgA5NH1ZCweT5UojNFKb+hvl9YoVjC+C0bIvrY+4h97lOTH2vOsU1ygwTAABAh4QJAACgQ8IEAADQIWECAADokDABAAB06JIHcKJ0xQRgKufcEdUMEwAAQIeECQAAoEPCBAAA0CFhAgAA6ND04cxMXbCnqPx4TP1anXPx57EY8xptO1685gDna8x9xLbrx1z3kFNcz8wwAQAAdEiYAAAAOiRMAAAAHRImAACADk0fVmLq4uxDPF/B92FMXUR5iOdr7rF+U5+fzneA8zX1df/Y7yPMMAEAAHRImAAAADokTAAAAB0SJgAAgA5NH1Zsim8m5jgt+Y3YAAB8khkmAACADgkTAABAh4QJAACgQ8IEAADQcdZNH9beQGGu8a19PxzC2hsozDW+te8HAIClmWECAADokDABAAB0SJgAAAA6JEwAAAAdNXPRtwrzM3COTSToa60d+oC4Jq7segyOiX9jj/M5Yu+YMe4zPuc8c1sirnAetsWzXePhXDF36nXvaq5r5tTXxz3uCboPNMMEAADQIWECAADokDABAAB0SJgAAAA65m76AAAAcDTMMAEAAHRImAAAADokTAAAAB0SJgAAgA4JEwAAQIeECQAAoEPCBAAA0CFhAgAA6JAwAQAAdEiYAAAAOiRMAAAAHRImAACADgkTAABAh4QJAACgQ8IEAADQIWECAADokDABAAB0SJgAAAA6JEwAAAAdEiYAAIAOCRMAAECHhAkAAKBDwgQAANDx/wFuemoSm0Z+kwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 864x288 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot3(img0, binary_erosion(img0, disk(1)), binary_dilation(img0, disk(1)), \\\n",
    "      'source', 'binary erosion', 'binary dilation')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0kAAAEtCAYAAADdgyslAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAEEJJREFUeJzt3XuMZnddx/HPtyy0XApdLBcLFNgSCaCItAgaUBAjYkT9w90AcmkUi0aNGApFQGgR6ErkEiQxJApIy3XFCyBeujGAhZDQNYIBg8DaUtpCC7ulUC7S8vOPc1a+DLOzl+7s2XZer2SSmefy+51nZp7fPu9zzjNbY4wAAAAwOW7pDQAAADiWiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBLABlZVl1bVz+7nukdV1aeO9jYtqar+saqevvR2wC3BzXV9qapzq+rC+fNTq+prVXWr+ev3V9UzDnPc7xmLY9umpTcAgGPTGOPfktx/6e04msYYj196G2AjuLmsL2OMzyW5w+Hct6ouTfKMMcbOmzoWR58jSaw7e0yAribr9u9PVdkBCBvUeq8vbBx+iUhVnVNVV1TVV6vqU1X12Ko6vqpeU1VXzh+vqarj59ufWVUXrxhjVNX95s/fVFV/XlXvq6rrkzymqm5bVa+sqsuq6itVdXFV3Xa+/SOq6sNVdW1VfayqHn20vwewwT2sqj5ZVXur6o1VdUKSVNWjq+rz+240nzpzdlV9fH4ev6PddnNVvbeqrpnHeW9V3bPd9/1V9bKq+lCSryd5dlXt6htRVc+uqr9bbQOr6k5V9ZdVddW8Xr20nf5yZlV9qKpeXVV7kpxbVcdV1QvnNefqqnpzVd1pvv0JVXVhVX15Xnc+WlV3a9v5jPnztca4z7zuPb2qPldVX6qqFxyxnwjcctwc1pf7VtUH5tdBFyU5uV2377n+fTtfquq0qvrXeS35UlW9papOmq+7IMmpSd5T0yl2z105VlWdUlXvrqo9VfWZqvrNNva5VfXOed35alV9oqrOOIzvP4dJJG1wVXX/JL+b5GFjjBOTPC7JpUlekOQRSR6S5EeT/HiSFx7C0E9O8rIkJya5OMmfJjk9yU8muXOS5yb5TlXdI8k/JHnpfPnZSd5VVXe5qY8NOGi/lum5f1qSH8raz/VtSX4+yX2TPDjJmfPlxyV5Y5J7Z3ph8I0kr1tx36cmOSvTuvDaJPetqge065+S5IL9zPtXSW5Icr8kP5bk55L09wU8PMnuJHfNtPacOX88JsmWTKe47Nuepye5U5J7JfmBJL81b+9Ka42xzyMznTL02CQvWvF4gJvH+vLWJLsyxdEfZ1ojDkYlOT/JKUkekGlNOTdJxhhPTfK5JE8YY9xhjPGKVe7/tiSfn+//q0leXlWPbdf/UpK3Jzkpybvz/Y+ZdSSSuDHJ8UkeWFW3HmNcOsb4bKZF7SVjjKvHGNckOS/TAnSw/n6M8aExxneS/G+SX0/y+2OMK8YYN44xPjzG+FamRet9Y4z3jTG+M8a4KMklSX7hSD5IYE2vG2NcPsbYkykwnrTGbV87xrhyvu17Mu1IyRjjy2OMd40xvj7G+Oo8zk+vuO+bxhifGGPcMD//35FpDUhVPSjJfZK8d+WE81Gexyd51hjj+jHG1UleneSJ7WZXjjH+bB77G5nWsFeNMXaPMb6W5A+TPHHeg/vtTHF0v3k92jXGuG6Vx7rWGPucN8b4xhjjY0k+lmmnEvBdx/r6cmqShyX5ozHGt8YYH5znPqAxxmfGGBfN97smyatW2a5VVdW9Mu1kOWeM8c0xxn8k+Yt872uti+fXRzdmCjzry1Ekkja4McZnkjwr056Pq6vq7VV1Sqa9Gpe1m142X3awLm+fn5zkhCSfXeV2906ydT7l5dqqujbTovGDhzAXcNP05+uBnutfaJ9/PfObkKvqdlX1+vnUtOuSfDDJSfW970ns8yTT0aEnV1VlemHwzvnFzUr3TnLrJFe1deL1mY4a7W/s1dawTUnulunFxj8neXtNpxO/oqpuvcq8a42xz6rfD+D/HevryylJ9o4xrl+xnQdUVXedXzddMW/XhWmn6h3AKUn2zNHX571H+3rl9+OE1U77Y32IJDLGeOsY45GZXoiMJH+S5Mr5631OnS9LkuuT3G7fFVV199WGbZ9/Kck3Mx1qX+nyJBeMMU5qH7cfY2w/7AcEHKp7tc/7c/1QPDvTaWcPH2PcMclPzZdXu01fFzLG+EimI82PynSK7v5Ohbk8ybeSnNzWiTuOMR60v7Gz+hp2Q5IvjjG+PcY4b4zxwEynAP9ikqetMu9+x9jPdgLf71hfX65Ksrmqbr9iOw/G+fO8D5636ylrbdMKVya5c1WduGLeKw5ybtaZSNrgqur+VfUzNf1Rhm9mOs/3xkznyb6wqu5SVScneVGmPSTJdErJg6rqIfObKs9da475lLs3JHnV/CbFW1XVT8xzXpjkCVX1uPnyE+Y3c95zrTGBI+p3quqeVXXnJM/PdJrKoTox0/px7TzOiw/yfm/OdJ79DWOMi1e7wRjjqiT/kuSVVXXHmv6gwmlVtdZpLW9L8gfzG7LvkOTlSd4xxrihqh5TVT8y74W+LtPpdzceyhgH+diAY399uSzTaf7nVdVtquqRSZ5wCNv1tXm77pHkOSuu/2Km9zOuNu/lST6c5Pz5tc+Dk/xGkrcc5NysM5HE8Um2Zzra84VMp688P9MfUrgkyceT/GeSf58vyxjjv5O8JMnOJJ/O9IcZDuTseZyPJtmT6WjVcfMi8cvznNdk2mP8nPjdhKPprZkiZPf88dLDGOM1SW6baS35SJJ/Osj7XZDkh7P/vbz7PC3JbZJ8MsneJH+dtU/LfcM85geT/E+mnUC/N1939/n+1yX5ryQfyHd3Ah3sGMDBuTmsL0/O9Mdf9mQKsDcf5PjnJXlokq9k+iNUf7Pi+vMz7XC+tqrOXuX+T8r0Xqkrk/xtkhfP783mGFBjrHUkEADWT03/FcDVSR46xvj00tsD3HJYX7gp7K0HYEm/neSjXsAA68D6wmHzFzIAWERVXZrpTc6/svCmALcw1hduKqfbAQAANE63AwAAaEQSAABAsy7vSaqqDXsO39atWxedf/v2Zf8P1p07dy46//Oe97zF5t67d+9icx8Lxhh14FvddBt5fVnaRj89u+qo/IqziqOxviy9tmz05xcb09Lr6lpriyNJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0m5begFua7du3Lzr/li1bFp1/8+bNi86/Z8+exebetm3bYnMnyY4dOxadH9ZbVS06/xhjsbmXfuwA62HJdfVAHEkCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0m5begPVw+umnLzb3li1bFps7SU477bRF59+9e/ei81900UWLzb3k712S7NixY9H5WX9jjEXnr6pF51/aRn/8t3RLP7+AY4sjSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANJuW3oD1sHnz5sXm3rVr12JzJ8nu3bsXnX9pS3//YT1V1dKbAAAbgiNJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0m5begPWwefPmxebeuXPnYnOz7M9+7969i80NAMCR40gSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAADNpqU3YD3s3bt3sblPP/30xeY+FmzevHnR+Zf8/u/YsWOxuYH1N8ZYbO6qWmxugI3IkSQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEBTY4wjP2jVkR/0EGzZsmWxuS+55JLF5k6SZz7zmYvOv3Xr1kXnX/Jnf8YZZyw297FgjFFHY56l15clrcd6fSiqjsqPeL82+uPfyI7S+rJh1xZYytLr6lpriyNJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0NcY48oNWHflBbybOOuusRec/55xzFp1/165di86/bdu2ReffyMYYdTTm2cjry9LW49+LQ1F1VH7FOAYdjfVl6bVl6ecXG9NGX1fXWlscSQIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoKkxxtLbAAAAcMxwJAkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACA5v8ASYCkl43lzBkAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 864x288 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot3(img1, binary_erosion(img1, disk(1)), binary_dilation(img1, disk(1)), \\\n",
    "      'source', 'binary erosion', 'binary dilation')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0kAAAEtCAYAAADdgyslAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAADwtJREFUeJzt3X2MpWdZx/HfBVsLSpU1rWKrkMyCqJhQM1tSDCraJkQjmGB28T2oyRKNIgrSphJTCZo1voQ/BOJqDISitAMEI/QPujFEEAQ6JNWQgNKhTWlBUKZCEWwpt3+cs+Ris6/tztyznc8n2WT2vDzXc2Zn7nm+55lztsYYAQAAYOFRs3cAAABgJxFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAM5LVfWRqnr27P0A5qqq11fVq6rqh6rqY+3yO6rq6oe4za/bFrvPntk7AAAPxRjjabP3Adg5xhjvSfLUh3LfqhpJnjLG+PjD3RaPDM4kseWq6tGz9wHYmarKk3UA7DgiiVTVNVV1d1V9oao+VlVXVdWFVfXqqrpn+efVVXXh8vYvrKr3HreNUVVPXn78+qp6XVXdXFVfTPKjVfXYqvqzqrqzqv6nqt5bVY9d3v7KqnpfVd1bVbf59Rk4/1XVpVX11qr6bFV9oqpevLz8+qp6S1XdUFWfT/LC06w3F1fVO5brw+eq6j1V9ajldV/7VZrTbOPZVfXJqnppVX2mqj5VVb886VMDPExV9QNV9eHlccuNSR6zvPzZVfXJk9znGVX1/uVa8qmq+ouq+obldf+0vNltVXVfVb3g+G1V1fdW1buX9/9IVT2vXff6qnpNVb1zuU8fqKp9W/cZYDuIpF2uqp6a5DeSXDHGuCjJc5LckeT3klyZ5PIkT0/yjCSvOItN/1ySP0xyUZL3JvnTJKtJfjDJtyZ5eZKvVtVlSd6Z5FXLy1+W5K1VdcnDfWzAHMuI+YcktyW5LMlVSV5SVc9Z3uSnkrwlyeOTvCmnXm9emuSTSS5J8u1JrksyTjD2dGvWE5J8y3J/fjXJa6pq78N/tMB2WobN25O8MYvjhrUkP30Gd30wyW8nuTjJM7NYl349ScYYP7y8zdPHGI8bY9x43MwLsljT3pXk25L8ZpI3LY+hjvnZJH+QZG+Sj2dxDMR5TCTxYJILk3xfVV0wxrhjjHF7kp9P8soxxmfGGJ/N4hv/F89iu38/xvjnMcZXk9yf5FeS/NYY4+4xxoNjjPeNMf4vyS8kuXmMcfMY46tjjFuS3JrkJ87lgwS21RVJLhljvHKMcf8YYyPJXyX5meX17x9jvH35Pf+lnHq9eSDJdyR50hjjgTHGe8YYJ4qk061ZDyyvf2CMcXOS++L1BnA+ujLJBUlevfx+fkuSD53uTmOM9THGv4wxvjLGuCPJXyb5kbOY+bgkh5dr2j8meUcWYXTM28YYHxxjfCWLJ38uP/OHxE4kkna55QsUX5Lk+iSfqao3V9WlSS5Ncme76Z3Ly87UXe3ji7M4FX77CW73pCQHlqev762qe5M8K4uDIuD89KQklx73fX1dFmeCkq9fH5JTrzd/ksWzsu+qqo2quvYkM0+3Zv338uDlmP/N4qAHOL9cmuTu454sufNkNz6mqr57+au7n17+qu8fZXF8cqYz71o+8dtnXtb+/un2sfXlEUAkkTHG344xnpXFgc1I8sdJ7ln+/ZgnLi9Lki8m+cZjV1TVE0602fbxfyX5cpIT/X7uXUneOMZ4fPvzTWOMww/5AQGz3ZXkE8d9X180xjh2hvj4M0EnXW/GGF8YY7x0jLGS5LlJfqeqrjrBzFOtWcAjx6eSXFZV1S574hnc73VJPprFO9h9cxZP3NSp7/I19yT5rmOvh2wz7z7D+3MeEkm7XFU9tap+bPkC5y8n+VIWv4L3d0leUVWXVNXFSX4/yQ3Lu92W5GlVdXlVPSaLs1AntXzm5W+S/PnyxdyPrqpnLmfekOS5VfWc5eWPWb5Y8ju35AED2+GDST5fizeFeezye/v7q+qKk9z+pOtNVf1kVT15eUD0+SzWpwfPZhvAI8r7k3wlyYurak9VPT+L1yCezkVZrCH3VdX3JPm1467/zyQrJ7nvB7J4gvjlVXXB8g2mnpvkzQ9h/zlPiCQuTHI4i7M9n87iBYnXZfFGCrcm+dck/5bkw8vLMsb49ySvTHI0yX9k8cYMp/Oy5XY+lORzWZytetQY464sXsR9XZLPZvEM9O/G1yact8YYD2ZxAHF5kk9ksb78dRZvnHAiJ11vkjwli7XmviwOjl47xnj3WW4DeIQYY9yf5PlJXphkM8kLkrztDO76sizeVOoLWbxG8sbjrr8+yRuWvyJ88AQzn5fkx7NYz16b5JfGGB99yA+EHa9O/PpXAACA3cmz9QAAAI1IAgAAaEQSAABAI5IAAAAakQQAANDs2YqNVtWufcu8AwcOTJ1/+PDc/4P16NGjU+dfe+2102Zvbm5Om70TjDHO9D/le1h28/oy2+x3Q92/f//U+evr61Pn72bbsb7MXltmHj/MPnZYWTnZf0+0PTY2NqbOn/n4Dx48ePobbaG1tbWp80+1tjiTBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAU2OMc7/RqnO/0fPE7bffPnX+ysrK1Plra2tT5x84cGDa7IMHD06bncz/3I8xajvm7Ob1ZbZDhw5NnX/kyJGp81dXV6fNXl9fnzZ7J9iO9WX22jLz+GH2sQO71+xjp5tuuumka4szSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgGbP7B3YCqurq9Nmr6ysTJudJPv27Zs6f2NjY+r8W265ZdrsmV93SbK2tjZ1Plvv0KFDU+cfOXJk6vzZ1tfXZ+8CW2j2Gj77+GE3O3r06NT5V1999bTZsx/77O+7U3EmCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACAZs/sHdgKe/funTZ7fX192uwk2djYmDp/ttmff9hKR44cmb0L8Ig189ghmfvza3V1ddrsZP7P7tnzZ9rNj/10nEkCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKDZM3sHtsLevXunzT569Oi02cz9t9/c3Jw2G4CHZ+bPj2R3Hz/Mfuyz/+1nzp/92HfysZMzSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANHtm78BW2NzcnDZ7dXV12uydYO/evVPnz/z8r62tTZsNbL2Z68v6+vq02bvFzGOHZO7X18rKyrTZyfxjp9nHLjPN/tzv5GMnZ5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAABNjTHO/Uarzv1Gz8LKysq02bfeeuu02Unyohe9aOr8AwcOTJ0/899+//7902bvBGOM2o45s9eXmQ4dOjR1/pEjR6bOX11dnTp/fX196vzdbDvWl3379k1dW2YfP+xmGxsbs3dh15p97HSqtcWZJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAmhpjnPuNVp37jZ4nDh06NHX+NddcM3X++vr61PkHDx6cOn83G2PUdszZzevLbFvx8+JsVG3Llxg70HasL7PXlpnHD4cPH542O0k2Nzenzp9t5uPfv3//tNk7wanWFmeSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoaowxex8AAAB2DGeSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGj+Hy8G4sRXUPq7AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 864x288 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot3(img1, erosion(img1, disk(1)), dilation(img1, disk(1)), \\\n",
    "      'source', 'erosion', 'dilation')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0kAAAEtCAYAAADdgyslAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAD69JREFUeJzt3XuMpXddx/HPtyy0XEoXLBdLaUGIBFDcoggasFtqRIxY/1ADyKVRNBo11lAgAtItAkUi0CCJIVGotFwj6gLiBWIWLIQEMAsGDHKxpbSFFtqlUArS8vOP51n5Mp2dzm539tndeb2SSWbnnPP8fjPd89vzfp7fmdYYIwAAAEyOWXoCAAAAhxORBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAtjEquqyqvrZfdz2uKr69KGe05Kq6p+q6plLzwOOBkfq+lJVO6rqkvnzU6rqG1V1h/nPu6rqWQd43O87Foe3LUtPAIDD0xjj35M8ZOl5HEpjjCcuPQfYDI6U9WWM8YUkdzuQx1bVZUmeNcZ43+09FoeeK0lsOGdMgK4mG/bvT1U5AQib1EavL2we/hKRqnpeVV1ZVV+vqk9X1ZlVdWxVXVhVV80fF1bVsfP9z66qS1ccY1TVg+fPL6qqv6yq91TVjUnOqKo7V9Urq+ryqvpaVV1aVXee7/+YqvpQVe2pqo9X1fZD/TOATe5RVfWpqrq+qt5QVcclSVVtr6ov7r3TvHXm3Kr6xPw8flu77z2q6t1Vde18nHdX1cntsbuq6qVV9cEk30zy7Kr6WJ9EVT27qv5htQlW1QlV9ddVdfW8Xr2kbX85u6o+WFWvrqrrkuyoqmOq6oXzmnNNVb2xqk6Y739cVV1SVV+d152PVNV92jyfNX++1jEeMK97z6yqL1TVV6rqBQftvwgcPY6E9eWBVfX++XXQe5Oc2G7b+1y/1cmXqnpQVf3bvJZ8pareVFVb59suTnJKknfVtMXuuSuPVVUnVdU7q+q6qvpsVf1WO/aOqnr7vO58vao+WVU/cQA/fw6QSNrkquohSX4/yaPGGMcneUKSy5K8IMljkmxL8mNJfjLJC/fj0E9N8tIkxye5NMmfJ/nxJD+d5J5Jnpvku1V1vyT/mOQl89fPTfKOqrrX7f3egHX79UzP/Qcl+eGs/Vz/tSQ/n+SBSR6R5Oz568ckeUOSUzO9MLgpyWtXPPbpSX4707rwmiQPrKqHttufluTifYz7N0luTvLgJKcl+bkk/X0Bj07y+ST3zrT2nD1/nJHkhzJtcdk7n2cmOSHJ/ZP8QJLfmee70lrH2OuxmbYMnZnkRSu+H+DIWF/enORjmeLoTzOtEetRSS5IclKSh2ZaU3YkyRjj6Um+kORJY4y7jTFescrj35Lki/PjfyXJy6rqzHb7LyV5a5KtSd6ZW3/PbCCRxC1Jjk3ysKq64xjjsjHG5zItai8eY1wzxrg2yfmZFqD12jnG+OAY47tJ/jfJbyT5wzHGlWOMW8YYHxpjfDvTovWeMcZ7xhjfHWO8N8lHk/zCwfwmgTW9doxxxRjjukyB8ZQ17vuaMcZV833flelESsYYXx1jvGOM8c0xxtfn45y+4rEXjTE+Oca4eX7+vy3TGpCqeniSByR598oB56s8T0xyzhjjxjHGNUleneTJ7W5XjTH+Yj72TZnWsFeNMT4/xvhGkj9O8uT5DO53MsXRg+f16GNjjBtW+V7XOsZe548xbhpjfDzJxzOdVAK+53BfX05J8qgkfzLG+PYY4wPz2LdpjPHZMcZ758ddm+RVq8xrVVV1/0wnWZ43xvjWGGN3kr/K97/WunR+fXRLpsCzvhxCImmTG2N8Nsk5mc58XFNVb62qkzKd1bi83fXy+WvrdUX7/MQkxyX53Cr3OzXJr85bXvZU1Z5Mi8YP7sdYwO3Tn6+39Vz/Uvv8m5nfhFxVd6mq181b025I8oEkW+v735PYx0mmq0NPrarK9MLg7fOLm5VOTXLHJFe3deJ1ma4a7evYq61hW5LcJ9OLjX9J8taathO/oqruuMq4ax1jr1V/HsD/O9zXl5OSXD/GuHHFPG9TVd17ft105TyvS9K26t2Gk5JcN0dfH/d+7c8rfx7Hrbbtj40hksgY481jjMdmeiEykvxZkqvmP+91yvy1JLkxyV323lBV913tsO3zryT5VqZL7StdkeTiMcbW9nHXMcbLD/gbAvbX/dvn/bm+P56dadvZo8cYd0/yM/PXq92nrwsZY3w405Xmx2XaoruvrTBXJPl2khPbOnH3McbD93XsrL6G3Zzky2OM74wxzh9jPCzTFuBfTPKMVcbd5zH2MU/g1g739eXqJPeoqruumOd6XDCP+4h5Xk9ba04rXJXknlV1/Ipxr1zn2GwwkbTJVdVDqurxNf1Shm9l2ud7S6Z9si+sqntV1YlJXpTpDEkybSl5eFVtm99UuWOtMeYtd69P8qr5TYp3qKqfmse8JMmTquoJ89ePm9/MefJaxwQOqt+rqpOr6p5Jnp9pm8r+Oj7T+rFnPs5563zcGzPts795jHHpancYY1yd5F+TvLKq7l7TL1R4UFWtta3lLUn+aH5D9t2SvCzJ28YYN1fVGVX1o/NZ6Bsybb+7ZX+Osc7vDTj815fLM23zP7+q7lRVj03ypP2Y1zfmed0vyXNW3P7lTO9nXG3cK5J8KMkF82ufRyT5zSRvWufYbDCRxLFJXp7pas+XMm1feX6mX6Tw0SSfSPKfSf5j/lrGGP+d5MVJ3pfkM5l+McNtOXc+zkeSXJfpatUx8yJx1jzmtZnOGD8n/m7CofTmTBHy+fnjJQdwjAuT3DnTWvLhJP+8zsddnORHsu+zvHs9I8mdknwqyfVJ/jZrb8t9/XzMDyT5n0wngf5gvu2+8+NvSPJfSd6f750EWu8xgPU5EtaXp2b65S/XZQqwN67z+OcneWSSr2X6JVR/t+L2CzKdcN5TVeeu8vinZHqv1FVJ/j7JefN7szkM1BhrXQkEgI1T0/8K4JokjxxjfGbp+QBHD+sLt4ez9QAs6XeTfMQLGGADWF84YH5DBgCLqKrLMr3J+ZcXngpwlLG+cHvZbgcAANDYbgcAANCIJAAAgGZD3pNUVfbwLWTXrl2Ljr9169ZFxz/vvPX+rxMOvp07dy429uFgjFG3fa/bz/qyeS29PbzqkPwVZxWHYn3ZzGvL0s8tNq+l19W11hZXkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaLYsPQEOrj179iw6/umnn77o+GecccZiY+/cuXOxsQEAOHhcSQIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADRblp7A0Wbbtm2Ljr99+/ZFx1/a7t27l54CHLXGGEtPAY5anl9weHElCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACAZsvSE9gI55xzzmJj79ixY7Gxk+SEE05YdPyl7dq1a+kpwIYZYyw9BThqeX4BnStJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0NcY4+AetOvgHPUJs3bp10fGvv/76Rcdf2mmnnbbY2Lt3715s7MPBGKMOxTibeX1Z2kb8e3EkqTokf8VZxaFYX5ZeWzb784vNael1da21xZUkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAACaLUtPAA6mbdu2LTb27t27FxsbAICDx5UkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQLNl6Qkcbfbs2bPo+Dt37lx0/LPOOmvR8bdv377Y2BdddNFiY8OhUFWLjj/GWHR82EhLPr88t+DWXEkCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKCpMcbScwAAADhsuJIEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQPN/nJeAXFJttK4AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 864x288 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot3(img2, binary_erosion(img2, disk(1)), binary_dilation(img2, disk(1)), \\\n",
    "      'source', 'binary erosion', 'binary dilation')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0kAAAEtCAYAAADdgyslAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAADiVJREFUeJzt3W3M3fVdx/HPlxUBHdIS0Ak6EthEnUlrGhZmphRpQjQykhmc90F9pNFZBafBRSqZBuNN+8BtEY25zJiOlC1r3HgwjWkEh9vo0mqWbLoNCLeCrhWYm9z9fHD+Xb5rerWl0P7b63q9kiZXz933fxquX37v8z/nUGOMAAAAsHDa3AcAAABwMhFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAE5JVfWZqto093EA86qqpap6V1X9QFV9rl3+QFVtPsbH/LrHYvVZM/cBAMCxGGO8Ye5jAE4eY4y7k1x6LPetqpHk9WOMz7/cx2JlcCaJ466qXjX3MQAnp6ryYh0AJx2RRKrqt6rqkap6uqo+V1VXVdUZVbW9qh6d/myvqjOm219fVfcc9Bijql43/bxUVe+tqruq6stJrqyqs6rqT6rqwar6n6q6p6rOmm5/eVV9vKr2V9Veb5+BU19VXVBVH6yqJ6vq/qp6+3T51qq6s6pur6qnklx/hPXmvKr6yLQ+fKmq7q6q06brvvZWmiM8xqaqeriqbqiqJ6rqsar6+Zn+aYCXqaq+r6o+Pe1b7khy5nT5pqp6eJn7vLGq7p3Wkseq6s+q6hum6/5putneqnqmqt528GNV1XdX1a7p/p+pqre065aq6t1V9dHpmD5RVZccv38BTgSRtMpV1aVJfiXJZWOMs5NcneSBJL+T5PIkG5KsT/LGJO98CQ/9U0l+P8nZSe5J8sdJNib5/iTnJnlHkher6sIkH03yrunyG5N8sKrOf7nPDZjHFDF/l2RvkguTXJVkS1VdPd3k2iR3Jlmb5P05/HpzQ5KHk5yf5FuT3JRkHGLskdas1yQ5ZzqeX0zy7qpa9/KfLXAiTWHz4STvy2LfsCPJjx3FXV9I8utJzkvypizWpV9OkjHGD063WT/GePUY446DZp6exZr2sSTfkuRXk7x/2kMd8JNJfi/JuiSfz2IPxClMJPFCkjOSfE9VnT7GeGCM8YUkP53kljHGE2OMJ7P4xf/Zl/C4O8cY/zzGeDHJs0l+IcmvjTEeGWO8MMb4+Bjj/5L8TJK7xhh3jTFeHGP8fZL7kvzIK/kkgRPqsiTnjzFuGWM8O8b4YpK/SPIT0/X3jjE+PP3OfyWHX2+eS/JtSS4aYzw3xrh7jHGoSDrSmvXcdP1zY4y7kjwTnzeAU9HlSU5Psn36fb4zyaeOdKcxxu4xxr+MMZ4fYzyQ5M+TXPESZr46ya3TmvaPST6SRRgd8KExxifHGM9n8eLPhqN/SpyMRNIqN31AcUuSrUmeqKoPVNUFSS5I8mC76YPTZUfrofbzeVmcCv/CIW53UZLrptPX+6tqf5I3Z7EpAk5NFyW54KDf65uyOBOUfP36kBx+vfmjLF6V/VhVfbGqfnuZmUdas/572rwc8L9ZbHqAU8sFSR456MWSB5e78QFV9Z3TW3cfn97q+wdZ7E+OduZD0wu/feaF7e+Pt5+tLyuASCJjjL8ZY7w5i43NSPKHSR6d/n7Aa6fLkuTLSb7xwBVV9ZpDPWz7+b+SfDXJod6f+1CS940x1rY/3zTGuPWYnxAwt4eS3H/Q7/XZY4wDZ4gPPhO07Hozxnh6jHHDGOPiJNck+Y2quuoQMw+3ZgErx2NJLqyqape99iju994kn83iG+y+OYsXburwd/maR5N8x4HPQ7aZjxzl/TkFiaRVrqouraofmj7g/NUkX8niLXh/m+SdVXV+VZ2X5HeT3D7dbW+SN1TVhqo6M4uzUMuaXnn5qyR/On2Y+1VV9aZp5u1Jrqmqq6fLz5w+LPntx+UJAyfCJ5M8VYsvhTlr+t3+3qq6bJnbL7veVNWPVtXrpg3RU1msTy+8lMcAVpR7kzyf5O1Vtaaq3prFZxCP5Ows1pBnquq7kvzSQdf/Z5KLl7nvJ7J4gfgdVXX69AVT1yT5wDEcP6cIkcQZSW7N4mzP41l8IPGmLL5I4b4k/5rk35J8erosY4x/T3JLkn9I8h9ZfDHDkdw4Pc6nknwpi7NVp40xHsriQ9w3JXkyi1egfzP+24RT1hjjhSw2EBuS3J/F+vKXWXxxwqEsu94keX0Wa80zWWyO3jPG2PUSHwNYIcYYzyZ5a5Lrk+xL8rYkHzqKu96YxZdKPZ3FZyTvOOj6rUn+enqL8I8fYuZbkvxwFuvZe5L83Bjjs8f8RDjp1aE//woAALA6ebUeAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAmjXH40GrylfmzWTXrl2zzl+7du2s82+++ebZZu/cuXO22SeDMcbR/k/5Xhbry+p13XXXzTp/x44ds85fzU7E+rKa15bVvndYv379rPP37t072+w5903J/Hunw60tziQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCsmfsAeGXt379/1vlXXHHFrPOvvPLK2Wbv3LlzttmwGmzevHnW+Tt27Jh1Phwvq33vMLf169fPNnvOfVNycu+dnEkCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0a+Y+gJVmw4YNs87ftGnTrPPntmfPnrkPAVasjRs3zjp/9+7ds86H42nLli2zzV7te4fVzL5pec4kAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQrJn7AI6HLVu2zDZ769ats81OknPOOWfW+XPbtWvX3IcAx826detmnX/ffffNOv+SSy6ZdT4r25x7h2Te/cNq3zusZvZNy3MmCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACAZs3cB3A8bN++fbbZS0tLs81Okn379s06f25r166d+xDguJn79/vcc8+ddf7FF18863xWtjn3DnPbtm3b3IfATOybludMEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAzZq5DwBeSRs2bJht9p49e2abDavBxo0bZ52/e/fuWeezsi0tLc02e9u2bbPNZl5z7puSk3vv5EwSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoFkz9wGsNPv37591/s6dO2edf+211846f9OmTbPNXlpamm02nAj79u2bdf7mzZtnnX/bbbfNOp+Vbc79w2rfO6xmc+6bkpN77+RMEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAATY0x5j4GAACAk4YzSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0/w94bYxnHd1E5AAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 864x288 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot3(img2, erosion(img2, disk(1)), dilation(img2, disk(1)), \\\n",
    "      'source', 'erosion', 'dilation')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### square(2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0wAAAEtCAYAAAA7quBhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAFbVJREFUeJzt3XusNHlZJ/Dvg4MMyk3E23ATIUsUZb3hZYO7IhtRs6z+4W4ULxDFjUaNGlCjos64KK5ZkaiJMVlRARGM7nph2VWMURyICWJEoxtXxIGRQQccRu7oDD//qHpjz5n+nbfP6e6q6u7PJ3mT857TXfWr6qqn6+lfPU9Xay0AAADc3T3mHgAAAMBSSZgAAAA6JEwAAAAdEiYAAIAOCRMAAECHhAkAAKBDwgTAXVTVTVX17zt/+5yq+oupxzSnqvo/VfXUuccBh+xQ40pVXV9VLxp/flhVvauqPmj8/+9W1dMvudy7LItlu2buAQBwOFprv5/k0XOPY0qttS+cewxwzA4lrrTW3pTkPpd5blXdlOTprbXf3nZZTM8ME7PxqQpwRQ329p5UVT4ghBOz77jC6XAQ0VVV31lVb66qd1bVX1TVE6vqXlX1vKq6Zfz3vKq61/j4p1XVjWeW0arqUePPP1dVP1VVL6+qdyd5QlXdu6p+tKreWFX/UFU3VtW9x8d/VlW9uqpur6rXVdXnTr0P4IQ9rqr+vKreXlU/W1XXJklVfW5V/c2VB4232Tyzqv5kPIdfuvLYD6uql1XVW8flvKyqHrLy3N+tqh+sqlcleU+SZ1TVa1cHUVXPqKpfXTfAqrp/Vf1MVb1ljFXPXrlV5mlV9aqq+rGqui3J9VV1j6p61hhvbq2qF1TV/cfHX1tVL6qqvx9jzmuq6qNWxvn08efzlvGxY8x7alW9qareVlXfs7NXBA7fIcSVR1TV743XPq9I8qCVv105x+/2AUxVPbKqfmeMIW+rql+oqgeMf3thkocl+Y0absP7jrPLqqrrqurXq+q2qnp9VX3dyrKvr6pfGuPNO6vqz6rq0y+x/7kkCRNrVdWjk3xTkse11u6b5ElJbkryPUk+K8knJ/nXST4jybMusOinJPnBJPdNcmOS/57k05L8myQPTPIdST5QVQ9O8r+TPHv8/TOT/EpVfcS22wZs5CsynPePTPKvcv55/p+TfEGSRyR5bJKnjb+/R5KfTfLwDBcL703yk2ee+1VJ/kuGmPDjSR5RVR+/8vevTPLCznp/PskdSR6V5FOSfH6S1XqCz0zyhiQfmSHuPG3894QkH5fhdpgr43lqkvsneWiSD0/y9eN4zzpvGVc8PsPtRU9M8n1ntgdO2SHElRcneW2GROm/ZogNm6gkz0lyXZKPzxBLrk+S1tpXJXlTkie31u7TWvuRNc//xSR/Mz7/S5P8UFU9ceXv/zHJS5I8IMmv5+7bzB5JmOi5M8m9knxCVd2ztXZTa+2vMgS7H2it3dpae2uSGzIEpk39WmvtVa21DyT5xyRfk+RbWmtvbq3d2Vp7dWvt/RmC2ctbay9vrX2gtfaKJH+Y5It2uZFA10+21m5urd2WIdn48nMe++OttVvGx/5Ghg9U0lr7+9bar7TW3tNae+e4nH935rk/11r7s9baHeO5/9IM53+q6jFJPjbJy86ucJz9+cIk39pae3dr7dYkP5bky1Yedktr7SfGZb83Q/x6bmvtDa21dyX5riRfNn7C+08ZEqVHjbHota21d6zZ1vOWccUNrbX3ttZel+R1GT5cApYfVx6W5HFJvre19v7W2ivHdV9Va+31rbVXjM97a5LnrhnXWlX10AwftHxna+19rbU/TvI/ctfrqxvHa6I7MyR74sqEJEys1Vp7fZJvzfDpyK1V9ZKqui7DJx9vXHnoG8ffbermlZ8flOTaJH+15nEPT/Kfxltjbq+q2zMEk4+5wLqAy1s9V692nv/tys/vyVjIXFUfUlU/Pd6+9o4kr0zygLpr/eLqepJh1ugpVVUZLhZ+abzgOevhSe6Z5C0rMeKnM8wm9Za9Ln5dk+SjMlyA/GaSl9Rwu/GPVNU916z3vGVcsXZ/AIuPK9cleXtr7d1nxnlVVfWR47XSm8dxvSgrt/NdxXVJbhsTwNX1Pnjl/2f3x7Xrbg1kPyRMdLXWXtxae3yGC5OW5L8luWX8/xUPG3+XJO9O8iFX/lBVH71usSs/vy3J+zJMzZ91c5IXttYesPLvQ1trP3zpDQIu4qErP6+e5xfxjAy3pn1ma+1+Sf7t+PtaecxqTEhr7Q8yzD5/ToZbeHu3zdyc5P1JHrQSI+7XWntMb9lZH7/uSPJ3rbV/aq3d0Fr7hAy3CP+HJF+9Zr3dZXTGCfyLpceVtyT5sKr60DPj3MRzxvU+dhzXV543pjNuSfLAqrrvmfW+ecN1s2cSJtaqqkdX1efV0NDhfRnuEb4zwz22z6qqj6iqByX5vgyfoiTDrSePqapPHoszrz9vHeNtec9P8tyx2PGDquqzx3W+KMmTq+pJ4++vHYtCH3LeMoGd+caqekhVPTDJd2e4peWi7pshdtw+Luf7N3zeCzLcn39Ha+3GdQ9orb0lyW8l+dGqul8NzRgeWVXn3QLzi0m+bSzqvk+SH0ry0tbaHVX1hKr6pPFT6ndkuEXvzossY8Ntg1O29Ljyxgy3/99QVR9cVY9P8uQLjOtd47genOTbz/z97zLUPa5b781JXp3kOeP1zmOTfG2SX9hw3eyZhImeeyX54QyzQH+b4TaX787QhOEPk/xJkj9N8kfj79Ja+/9JfiDJbyf5ywxNHa7mmeNyXpPktgyzWPcYg8cXj+t8a4ZPk789jlmYyoszJCRvGP89+xLLeF6Se2eII3+Q5P9u+LwXJvnE9D8FvuKrk3xwkj9P8vYkv5zzb9t9/rjMVyb56wwfBn3z+LePHp//jiT/L8nv5V8+DNp0GcD5DiGuPCVDw5jbMiRjL9hw+Tck+dQk/5ChadX/PPP352T4wPn2qnrmmud/eYbaqluS/K8k3z/Wb7MA1dp5M4QAMK0avlrg1iSf2lr7y7nHAxw+cYVt+LQegKX5hiSvcVED7JC4wqXprgHAYlTVTRkKpb9k5qEAR0JcYVtuyQMAAOhwSx4AAECHhAkAAKBj6hom9/+dgOGLtGHQWtv3ASGunABxhVXiCruwLq4ce6mKbe47L66YYQIAAOiQMAEAAHRImAAAADokTAAAAB2zf3HtMRXyblo0t802r1vHnPtwm0LBYyoyvEBB4U7XMec+XPK5u+Sx7cLS4gBwnHb93jZVnDqm6wuWwQwTAABAh4QJAACgQ8IEAADQIWECAADomL3pw6EW5s1VYH1Mhd2Hui1zHbOHeq4A09t1Ef5U8edQ3xdOnWYOHDszTAAAAB0SJgAAgA4JEwAAQIeECQAAoGP2pg9zFniuKx5UcDqYYt/MWby5blsUkw7sGzgdmjlwHrEfBmaYAAAAOiRMAAAAHRImAACADgkTAABAR01c0Kd68AQo7mVVa23fB8Td4sqxH4Ma1nDq5ogrHB9xk1XnxRUzTAAAAB0SJgAAgA4JEwAAQIeECQAAoOOauQcAsGun+O30x77NirOBpds0Dl8knu1jmUte79zr7jHDBAAA0CFhAgAA6JAwAQAAdEiYAAAAOg666cO6Yq9tC5/3scwlr3fudV/WuvFtW/y3j2Uueb1zr3upNt3+i5wj+1jmktc797oB9m2b2LW0uDfFtuzjWnPK5hBmmAAAADokTAAAAB0SJgAAgA4JEwAAQIeECQAAoONguuRt2l3jIl045uoGtk2nkG27jGy6zUvqcLXN69zbjrk61GzTlW7bjnabbvOUXWfmts02LG37p9iWfcTSpe1HYBkOsYPvsTvl18QMEwAAQIeECQAAoEPCBAAA0CFhAgAA6DiYpg+7LljvPVYB8mBJ+2bXjSp6jz2VwsWrsW+O25LObY6TY4xdWPr7zlTH9JL2w5LGMjUzTAAAAB0SJgAAgA4JEwAAQIeECQAAoKMmLuA63WqxS5iq+H7TwsVN1624d3tTFU1P8Zq21vZ9QCw+ruz6HNtmHdvojc85z9TEFXZB7LqYqa4D53ovPC+umGECAADokDABAAB0SJgAAAA6JEwAAAAdmj6wc4ooWaU4e1mmaEBxkfVsum5xhVXiCrsgrrBK0wcAAIBLkDABAAB0SJgAAAA6JEwAAAAd18w9AACmM1Wjn4kbCgFwYPbxPrGvRh5mmAAAADokTAAAAB0SJgAAgA4JEwAAQIemDydmm2I4RdyHbZvXz7ehc55dHx9iDbBpXDnEeHFIzQ5WHcJ1xLox7mLdZpgAAAA6JEwAAAAdEiYAAIAOCRMAAEDHSTR92LbYa66CwqUVPK4bz9KLLbcd31zNDjYd96EXUR6yfWz/FOfTITR+OfVj69AtLX7BRRxCY4MpLP36bmpmmAAAADokTAAAAB0SJgAAgA4JEwAAQMfsTR/mbGywpEL2pTV4mGLdcxYGL6mBxSEUSB9aIeshNDaYwqG9brBkS3rfmNNU++EU9y3LZYYJAACgQ8IEAADQIWECAADokDABAAB0zN70YWmWXmS4j2LLTZtfHHvB69IL5PfRpGTT13RJDVLgohy/g6n2w7Hu22N6v9uG/cApMsMEAADQIWECAADokDABAAB0SJgAAAA6Ftn0Yc6CwrkaG0y1zZuuZ67XYM5i4bkKw6fa5k3Xc6wF25wux/TAfgA2dUzXP7tghgkAAKBDwgQAANAhYQIAAOiQMAEAAHQssunDrhsv+Fbqw7HrxguHVFAIu7a0ZjLbcC4DHJep3qN28f5hhgkAAKBDwgQAANAhYQIAAOiQMAEAAHTM3vRh10X+6/SWd4rNIDbdt5u+Ltvsw10vb53e8k6xgHzTfbvp63KK+5D1HAsAx2XX13eHzgwTAABAh4QJAACgQ8IEAADQIWECAADokDABAAB0zN4l7xg7aVzGNt3rLmJJnQGXNJY5bdO97iKca6dnm9fc+Ql3NUVnVzgG++isO/f5Z4YJAACgQ8IEAADQIWECAADokDABAAB0zN70gb5Ni+YuUvQ2d9Ecm9n0dbpIEeU+ijA5XvuIFXMdg459gMM297WqGSYAAIAOCRMAAECHhAkAAKBDwgQAANBxEk0feoViin4Hx9wIovcaH8v2bUsxPFNybHHIvG+wK8d83XWszDABAAB0SJgAAAA6JEwAAAAdEiYAAICOmrjITEXbCVDYzarW2r4PiLvFlW2OwW1j4pzrPmvX56IGOizFHHGF4yN2seq8uGKGCQAAoEPCBAAA0CFhAgAA6JAwAQAAdFwz9wAAdm3Ob0xf0re1TzWWKdajOBuAuZhhAgAA6JAwAQAAdEiYAAAAOiRMAAAAHRImAACADgkTAABAh4QJAACgQ8IEAADQIWECAADokDABAAB0SJgAAAA6JEwAAAAdEiYAAIAOCRMAAECHhAkAAKBDwgQAANAhYQIAAOiQMAEAAHRImAAAADqumXsAAACwBK21uYewV1V1t9+d4jZflBkmAACADgkTAABAh4QJAACgQ8IEAADQMXvTh10UYi3FpkVzp7jNu37u0mz6mp7iNs9hyWPbhXXH0bFvMzC9Q40rx/ReyzKYYQIAAOiQMAEAAHRImAAAADokTAAAAB2zN31Y5xCK9Q61EHLpDmG/HsLxCZy2Q41Th/AeAJweM0wAAAAdEiYAAIAOCRMAAECHhAkAAKBjkU0fFH2erkMtVAaAY+M9GQZmmAAAADokTAAAAB0SJgAAgA4JEwAAQMfsTR9OsaDwFLd5nVNs7nGK2wxMT6yByznFc+cUt/mizDABAAB0SJgAAAA6JEwAAAAdEiYAAIAOCRMAAEDH7F3yAHbtFDtRHvs26+IELN2mcfgi8Wwfy9zletfZdixzrrvHDBMAAECHhAkAAKBDwgQAANAhYQIAAOg46KYP6wq7ti183scyL7vedfYxljnXfVnrxrKPIsMpCs3nKuice91LtY/zYa5zbJvXbR+xdKp1A8dprmu0dbZZ7z7GvM01zK7Hs+lYlrYfzmOGCQAAoEPCBAAA0CFhAgAA6JAwAQAAdBxM04dNi8UuUlS2aRHYrgvVtlnetmPZpvBtrmLLbV7n3vZuOu5dFw/uuihzH98Wvo91L9Wuz4dt7To2TDWWKdbNYTvWGML+zHWNRt8p72szTAAAAB0SJgAAgA4JEwAAQIeECQAAoONgmj7s41uD5/wmYja3bYOHqZYJh8Bxzr7N1VSH47L0a7SpmvTM1fhnnVO+HjbDBAAA0CFhAgAA6JAwAQAAdEiYAAAAOmriAq7TrRa7ik2L9fbxeu163Yp2L2bO/TrFultr+z4gdnpS7ONcnOL8nrPg1znP1ey6wcOhxRWWSezqm/P9aK51nxdXzDABAAB0SJgAAAA6JEwAAAAdEiYAAIAOTR/YOUWUrFKcPZ8ti18XtW5xhVXiCrsgrrBK0wcAAIBLkDABAAB0SJgAAAA6JEwAAAAd18w9AAD2Y+KmPotZNwDLN9X7xC6ae5hhAgAA6JAwAQAAdEiYAAAAOiRMAAAAHZo+nJhdf6u1wu7DsevXyjekc4VjAdi1beLK0q9N9jG+pV/fHfr7hBkmAACADgkTAABAh4QJAACgQ8IEAADQcRJNH7YtNJuiePCYixvntO2+maJIcZsxHnoR5SHbx75fepGtWMMmxLTjtulrdKjxwvHLOmaYAAAAOiRMAAAAHRImAACADgkTAABAx+xNH+YsHly3TAV709n0Nd3Ha7JumYdaoMrdaaICdzdnzD0G3jfgdJlhAgAA6JAwAQAAdEiYAAAAOiRMAAAAHbM3fQCAfdLgB3ZnqkYXGmqwJGaYAAAAOiRMAAAAHRImAACADgkTAABAh6YPC6GI8nRNVXyuyB1Om1gDbGqK8/iQYoUZJgAAgA4JEwAAQIeECQAAoEPCBAAA0CFhAgAA6NAlD+BIHVP3zUPqpgTA1U31HrWL9w8zTAAAAB0SJgAAgA4JEwAAQIeECQAAoEPThxOz68LpqQr22N6uXytF+FzhWAA4Lq4Z7soMEwAAQIeECQAAoEPCBAAA0CFhAgAA6ND0YSG2KYbbtjBv0+cfesHeUm3z+m37mmz6fM09DpMmLwAcg3XvZ1O+J5lhAgAA6JAwAQAAdEiYAAAAOiRMAAAAHZo+LNi6YjaNF07D3MWNsA9iGgCXMfc1kBkmAACADgkTAABAh4QJAACgQ8IEAADQcdJNH5ZebDzV+Ja+H/Zh7uLBq5lqfEvfDxyXU4w1HA/xEk6XGSYAAIAOCRMAAECHhAkAAKBDwgQAANAxe9OHKYooD6FQ8xSL/KcoAD+EInPNPWC/HPvAplybDOyHuzLDBAAA0CFhAgAA6JAwAQAAdEiYAAAAOmpJTQAAAACWxAwTAABAh4QJAACgQ8IEAADQIWECAADokDABAAB0SJgAAAA6JEwAAAAdEiYAAIAOCRMAAECHhAkAAKBDwgQAANAhYQIAAOiQMAEAAHRImAAAADokTAAAAB0SJgAAgA4JEwAAQIeECQAAoEPCBAAA0CFhAgAA6JAwAQAAdEiYAAAAOiRMAAAAHf8M3rRBDDTPMj8AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 864x288 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot3(img0, binary_erosion(img0, square(2)), binary_dilation(img0, square(2)), \\\n",
    "      'source', 'binary erosion', 'binary dilation')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### rectangle(1,2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0wAAAEtCAYAAAA7quBhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAFrJJREFUeJzt3XvQfXs9B/DPJ7+ck66S2+muRkPkGjGhZIQR/sCQSw0xDAZTGCTnENFQDWZMM4RKyrhLLhlDTo2ZZJSJQXLq6KRTTkf3OKevP9b6jX1+Z3+es/ez1l6XvV+vmd/M83uevdd9ffZ677U/352ttQAAAODWbjf3AgAAACyVwAQAAFAQmAAAAAoCEwAAQEFgAgAAKAhMAAAABYEJgFvIzGsy8wuKv31OZv7z1Ms0p8z848x83NzLAWu21rqSmVdm5vP7n++Tme/KzA/q//+XmfmEc073FtNi2S7MvQAArEdr7a8j4kFzL8eUWmtfPPcywDFbS11prb0xIu50nudm5jUR8YTW2p8PnRbTc4eJ2XhXBbgoOwd7TcpMbxDCiTl0XeF0OIgoZeYPZOabMvOdmfnPmfmozLwsM5+Vmdf1/56VmZf1j398Zl59yTRaZj6w//lXM/MXM/MlmfnuiHhkZt4hM382M9+Qmf+dmVdn5h36xz8sM1+RmTdm5qsz8xFTbwM4YQ/NzH/MzLdn5q9k5uUREZn5iMz8j4sP6j9m86TMfE1/Dr9o47Efmpkvzsy39tN5cWbea+O5f5mZP5GZL4+I90TEEzPzVZsLkZlPzMzf27aAmXnXzPzlzHxzX6ueuvFRmcdn5ssz85mZeUNEXJmZt8vMJ/f15vrMfG5m3rV//OWZ+fzM/K++5rwyMz9yYzmf0P981jTu19e8x2XmGzPzbZn5w6PtEVi/NdSV+2fmX/XXPi+NiHts/O3iOX6rN2Ay8wGZ+Rd9DXlbZv56Zt6t/9vzIuI+EfGH2X0M7/svnVZmXpGZf5CZN2Tm6zLzWzamfWVm/mZfb96Zma/NzE8/x/bnnAQmtsrMB0XEd0bEQ1trd46IR0fENRHxwxHxsIj45Ij4pIj4jIh48h6TfmxE/ERE3Dkiro6In4mIT4uIz46Iu0fE90fEBzLznhHxRxHx1P73T4qI387MDx+6bsBOvi668/4BEfGxcfZ5/tUR8UURcf+IeEhEPL7//e0i4lci4r7RXSy8NyJ+4ZLnfkNEfGt0NeHnIuL+mflxG3//+oh4XjHfX4uImyLigRHxKRHxhRGx2U/wmRHx+oj4iOjqzuP7f4+MiI+J7uMwF5fncRFx14i4d0R8WER8W7+8lzprGhc9PLqPFz0qIp5yyfrAKVtDXXlBRLwquqD049HVhl1kRDwtIq6IiI+LrpZcGRHRWvuGiHhjRDymtXan1trTtzz/NyLiP/rnf2VE/GRmPmrj718WES+MiLtFxB/ErdeZAxKYqNwcEZdFxMdn5u1ba9e01v4tumL3Y62161trb42Iq6IrTLv6/dbay1trH4iI/4mIb4qI726tvam1dnNr7RWttfdHV8xe0lp7SWvtA621l0bE30bEl4y5kkDpF1pr17bWbogubHztGY/9udbadf1j/zC6N1SitfZfrbXfbq29p7X2zn46n3fJc3+1tfba1tpN/bn/oujO/8jMB0fE/SLixZfOsL/788UR8T2ttXe31q6PiGdGxNdsPOy61trP99N+b3T16xmttde31t4VET8YEV/Tv8P7v9EFpQf2tehVrbV3bFnXs6Zx0VWttfe21l4dEa+O7s0lYPl15T4R8dCI+JHW2vtbay/r532bWmuva629tH/eWyPiGVuWa6vMvHd0b7T8QGvtfa21v4+IX4pbXl9d3V8T3Rxd2FNXJiQwsVVr7XUR8T3RvTtyfWa+MDOviO6djzdsPPQN/e92de3Gz/eIiMsj4t+2PO6+EfFV/UdjbszMG6MrJh+9x7yA89s8V2/rPP/PjZ/fE30jc2Z+SGY+u//42jsi4mURcbe8Zf/i5nwiurtGj83MjO5i4Tf7C55L3Tcibh8Rb96oEc+O7m5SNe1t9etCRHxkdBcgfxoRL8zu48ZPz8zbb5nvWdO4aOv2ABZfV66IiLe31t59yXLepsz8iP5a6U39cj0/Nj7OdxuuiIgb+gC4Od97bvz/0u1x+baPBnIYAhOl1toLWmsPj+7CpEXET0fEdf3/L7pP/7uIiHdHxIdc/ENmftS2yW78/LaIeF90t+YvdW1EPK+1dreNf3dsrf3UuVcI2Me9N37ePM/38cToPpr2ma21u0TE5/a/z43HbNaEaK39TXR3nz8nuo/wVh+buTYi3h8R99ioEXdprT24mnZsr183RcRbWmv/21q7qrX28dF9RPhLI+Ibt8y3nEaxnMD/W3pdeXNEfGhm3vGS5dzF0/r5PqRfrq8/a5kucV1E3D0z73zJfN+047w5MIGJrTLzQZn5+dkN6PC+6D4jfHN0n7F9cmZ+eGbeIyKeEt27KBHdR08enJmf3DdnXnnWPPqP5T0nIp7RNzt+UGZ+Vj/P50fEYzLz0f3vL++bQu911jSB0XxHZt4rM+8eET8U3Uda9nXn6GrHjf10fnTH5z03us/n39Rau3rbA1prb46IP4uIn83Mu2Q3GMMDMvOsj8D8RkR8b9/UfaeI+MmIeFFr7abMfGRmfmL/LvU7ovuI3s37TGPHdYNTtvS68oboPv5/VWZ+cGY+PCIes8dyvatfrntGxPdd8ve3RNf3uG2+10bEKyLiaf31zkMi4psj4td3nDcHJjBRuSwifiq6u0D/Gd3HXH4oukEY/jYiXhMR/xARf9f/Llpr/xIRPxYRfx4R/xrdoA635Un9dF4ZETdEdxfrdn3x+PJ+nm+N7t3k7wvHLEzlBdEFktf3/556jmk8KyLuEF0d+ZuI+JMdn/e8iPiEqN8FvugbI+KDI+IfI+LtEfFbcfbHdp/TT/NlEfHv0b0Z9F393z6qf/47IuKfIuKv4v/fDNp1GsDZ1lBXHhvdgDE3RBfGnrvj9K+KiE+NiP+ObtCq37nk70+L7g3nGzPzSVue/7XR9VZdFxG/GxE/2vdvswDZ2ll3CAFgWtl9tcD1EfGprbV/nXt5gPVTVxjCu/UALM23R8QrXdQAI1JXODejawCwGJl5TXSN0l8x86IAR0JdYSgfyQMAACj4SB4AAEBBYAIAAChM3cPk838noPsibei01g59QKgrJ0BdYZO6whjUFTadVVfcYQIAACgITAAAAAWBCQAAoCAwAQAAFGb/4tpjarjb9TuthqzztnnMuQ2HfI/XMX0H2K77YMg6b5vHnNtwyefukpdtDEurA0u31loz9mvFkHlMtQ2XfBwvednOMvbrzpB5TLUNXZvUlnYtMYUxjjt3mAAAAAoCEwAAQEFgAgAAKAhMAAAAhdkHfVhro9lczZ9rbTrdZq3rMtcxu9ZzZQ5r3VZrPSeWznatrfVcmYNtVTumc2yt6+L4PCx3mAAAAAoCEwAAQEFgAgAAKAhMAAAAhdkHfZizuW5bg9xam/3GNsW2mbNB8RS/6XpXx7BtnMdsmur43fW423V5juk4PoZ1mWId9jlWxz7e1mCKdTmm7cV43GECAAAoCEwAAAAFgQkAAKAgMAEAABRy4uY2nXQn4BiaexlPa+3QB4S6cgLUFTapK4xBXWHTWXXFHSYAAICCwAQAAFAQmAAAAAoCEwAAQEFgAgAAKAhMAAAABYEJAACgIDABAAAUBCYAAIDChbkXYIht39Dc2rAv5z7ENJc837nnfV7blm/oN3YfYppLnu/c816qXdd/n3PkENNc8nznnjfAHHatZ/u8zh5imkue79zzrrjDBAAAUBCYAAAACgITAABAQWACAAAoCEwAAACF1YySt+tIbvuM+DbXaGBDRqUbOqLdruu8pBGuhuznaj3mGplryKh0Q0e023Wdlzg6zaEMWYelrf8Ux9EhaunStuPYTnEkyiH19di3zXmtcTTboYYcC3Ndm0zx3EOY4hr0EMfwlNcr7jABAAAUBCYAAICCwAQAAFAQmAAAAAqrGfRh7Ib16rEaTjtL2jZjD1RRPXZpTZhzOaVtM9V6qSuny77v2A7DHGsN3teStsNUx/SS1vmUucMEAABQEJgAAAAKAhMAAEBBYAIAACjkxM1kOtf2MFXz/a6Ni1N+o/Kpm2rQjSn2aWvt0AeEurKHsc/3ofPZdd7qyrqNvU/VldM09nWRurKfqa4D59qnZ9UVd5gAAAAKAhMAAEBBYAIAACgITAAAAAWDPjA6TZRs0pzNGNQVNqkrjEFdYZNBHwAAAM5BYAIAACgITAAAAAWBCQAAoCAwAQAAFAQmAACAgsAEAABQEJgAAAAKAhMAAEDhwtwLwLSGfKt1a774fM2G7D/fhs5Zxj4+1Jp123X/qSucZdfjQ71Yr0Psu0PVFXeYAAAACgITAABAQWACAAAoCEwAAACFkxj0YWgD2FwNhUtreNy2PEtvthy6fHM1JS+taXrb8px6w/Yh1n+K82kNA7+s9dha2nm7iznPbXXl1vZZ/yW9/s55fXCK1ybbTHHurGHwqEPVFXeYAAAACgITAABAQWACAAAoCEwAAACF2Qd9mHNggyU1nC5tgIcp5j1ng/SSmkTX0Ci+tkbsNQxsMIW1boexa/NUtX7p58mSXvPmnvd5TPG6sYb6s6TXz0PMew0DG0xhrdvhUPN2hwkAAKAgMAEAABQEJgAAgILABAAAUJh90IelWXrD5SGaLXdtBF5ao+fYlt6seYiG7V336dKaxddmjefJnOe7Y2s6U9WVqebNPKa6Nplq3kuyxnNiznP7UPN2hwkAAKAgMAEAABQEJgAAgILABAAAUMiJG+NuNbMpmvWqZq9dG8OOqXlwCktq+K323a77eY3NlkvTWjv0Rtyprmx94sK+yVytqTkX2TRHXeH4DLnmG1KTDlHr1cjhzqor7jABAAAUBCYAAICCwAQAAFAQmAAAAAoX5l6AbcYeeEEj9XqMPfCCJkjOMqTW7FNXHIfzMqALYzAoFJuG1JV96o9jbBncYQIAACgITAAAAAWBCQAAoCAwAQAAFGYf9GGKZtxqeqfYSLfrtt11vwzZhlM00FbTO8WG7yHfXr62pvk568opWvq2WPryHZu11YtdneI1w5zmGmRjzmuTUzPVdhij/rjDBAAAUBCYAAAACgITAABAQWACAAAoCEwAAACF2UfJO4aRc8YwZPS6fSxpZJYlLcuchoxet49TOtdOaV3PMmQ77HN+jj1a5tiOddS2fU21HY512841atvSTLUd5tq2p7hPt5mirg+tFVPuK3eYAAAACgITAABAQWACAAAoCEwAAACF2Qd9oLZrg+4+TW+aVtdh1/20T8Okxnf2sU+tmGrQGoZxvjMG5zG3Zarr1ym5wwQAAFAQmAAAAAoCEwAAQEFgAgAAKJzEoA9DG5WP3TEPBFHt42NZv6EMBHF+6sr+TnHbnOI5dorrfGjH/DpdOeZ1PuZ1O5S5t487TAAAAAWBCQAAoCAwAQAAFAQmAACAQk7cRKWj7QRo7mVTa+3QB4S6cgLUFTapK4xBXWHTWXXFHSYAAICCwAQAAFAQmAAAAAoCEwAAQEFgAgAAKAhMAAAABYEJAACgIDABAAAUBCYAAICCwAQAAFAQmAAAAAoCEwAAQEFgAgAAKAhMAAAABYEJAACgIDABAAAUBCYAAICCwAQAAFAQmAAAAAoCEwAAQEFgAgAAKAhMAAAABYEJAACgIDABAAAUBCYAAICCwAQAAFAQmAAAAAoCEwAAQEFgAgAAKFyYewEyc+5FGE1rbafHneI6j/3cpdl1n57iOs9hycs2hm3H0bGv8xBrPe+G7NNTXOdDW/KynWXIsXCK67zWc2dX2/bpKa7zvtxhAgAAKAhMAAAABYEJAACgIDABAAAUZh/0YZs1NJ+ttRFy6dawXddwfLKbpe3LNRz/a2S7Moal1QvmtbS64vg8LHeYAAAACgITAABAQWACAAAoCEwAAACFRQ76sLRGOqajaZEpqTWnYa11xfG5LGvdH2s9/pfOdj0t7jABAAAUBCYAAICCwAQAAFAQmAAAAAo5cdOaDrkTsNbGWA6jtXboA0JdOQHqCpvUFcagrrDprLriDhMAAEBBYAIAACgITAAAAAWBCQAAoCAwAQAAFAQmAACAgsAEAABQEJgAAAAKAhMAAEBBYAIAAChcmHsBhsjMW/2utba4aZ53vtscYlnmnPd5bVuWXddjymmed77bHGJZ5pz3Uh3ifJjrHBuy3w5RS6eaN9vNVePmnjfjmus6ae557+IQr6lzvU4P2a6HuB6bat4Vd5gAAAAKAhMAAEBBYAIAACgITAAAAIWcuFnu3DPbtdFvn4bAuZqShzQtDm14nKIpcK557LMsczUUDtleQ7f1XAMKtNYO3dk9al1ZmrEbfsc+Bvd5/ik6xYFVpqivS64rh7DGwZmGmuIa7dgHpxn7uJnz+nWIMeqKO0wAAAAFgQkAAKAgMAEAABQEJgAAgMKFuRdgV0MHeJhqmoxv6AAPU02T9TnEAAjbnGLDNh01hDGoDcsz9gAIlVMcOGaJ3GECAAAoCEwAAAAFgQkAAKAgMAEAABRy4kZCXYuFOZvCp/gmaGpzbteJvg390AeEulKY81vsx563ujLcnNt17HmrK/OZc3CsseetrtTm3K5zzfusuuIOEwAAQEFgAgAAKAhMAAAABYEJAACgYNAHRqeJkk2asxmDusImdYUxqCtsMugDAADAOQhMAAAABYEJAACgIDABAAAUBCYAAICCwAQAAFAQmAAAAAoCEwAAQEFgAgAAKFyYewGY1tjfat2aL0Nfi7H3lW9I5yLHApuG1BrHEhcNORZcm6zDVPtpjLriDhMAAEBBYAIAACgITAAAAAWBCQAAoHASgz4MbfaaoilNc+NhDN02UzQga5Bep0Ns+6UPzKHWnG3X7TPXebv05WO/bT/X+bjrMqoXtUNsm6XX+7XXFXeYAAAACgITAABAQWACAAAoCEwAAACF2Qd9mLN5cNs0196UtiZzNiBvm6YG1eNhEBU2TVXrvX4ctyleN9Sf+RmIiW3cYQIAACgITAAAAAWBCQAAoCAwAQAAFAQmAACAgsAEAABQEJgAAAAKAhMAAEBBYAIAAChcmHsB6Ez17d6+RXx5pvpmcN9AfnqOqa44fodTaxjDMdUValOcx2uqFe4wAQAAFAQmAACAgsAEAABQEJgAAAAKAhMAAEBBYAIAACgITAAAAAWBCQAAoCAwAQAAFAQmAACAwoW5F4BpZeao02utjTo9DmfsfTX2scR6ORbYNKTWOJa4aMix4NpkHabaT2PUFXeYAAAACgITAABAQWACAAAoCEwAAAAFgz4sxJzNjbs+XzPuYczZIL3r8zXQrtPSBnnZ9nx1ZTpT1Xr79LhN9brhdWdeSxsoatvzpzxG3GECAAAoCEwAAAAFgQkAAKAgMAEAABQM+rBgGqRP19zNjcDx8frBGLwWcYrcYQIAACgITAAAAAWBCQAAoCAwAQAAFE560IelN8BOtXxL3w6HsPSm1amWb+nbgeNy7LVm1/Npru0w1fINqSvHfozcln3Wf676vesyDl2+IceC17bjM/c+dYcJAACgIDABAAAUBCYAAICCwAQAAFDIiZuodOGdgFNv2uWWWmuHPiDUlROgrrBJXWEM6gqbzqor7jABAAAUBCYAAICCwAQAAFAQmAAAAApTD/oAAACwGu4wAQAAFAQmAACAgsAEAABQEJgAAAAKAhMAAEBBYAIAACgITAAAAAWBCQAAoCAwAQAAFAQmAACAgsAEAABQEJgAAAAKAhMAAEBBYAIAACgITAAAAAWBCQAAoCAwAQAAFAQmAACAgsAEAABQEJgAAAAKAhMAAEBBYAIAACgITAAAAIX/A3STEAuG/x44AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 864x288 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot3(img0, binary_erosion(img0, rectangle(1,2)), binary_dilation(img0, rectangle(1,2)), \\\n",
    "      'source', 'binary erosion', 'binary dilation')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## opening/closing"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0wAAAEtCAYAAAA7quBhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAFwBJREFUeJzt3X/wbXtZF/D3gwcvkFeRkdILSE5MjOEIo1HY+IOCGcwZIcl+qFn2UysrJghnzElitGimjGEcq8mCRBFQCoWwqcairoShFBaNJsplrlwcMLjBSAhcPv2x1sF9z92f79n7rL33Wmvv12vmzpx7zt5rfdaP/ezPsz/7eXa11gIAAMADPWjuAQAAACyVhAkAAKBDwgQAANAhYQIAAOiQMAEAAHRImAAAADokTAAkSarqrqp6euffvqyqfuHUY1qSqnp7VT117nHAWq09xlTVU6vqVyZuY/HHyQNdm3sAACxfa+0/J3n83OOYU2vtCXOPAc7VpcSYSznOc2OFidlU1afMPQZgXjXwXgQchRjDIbiB6Kqqb6uqd1fVh6rqF6rqaVV1W1W9uKruGf97cVXdNj7+m6rqzhu20arqceOfX1ZV/6iq3lBVv57k91fVQ6vqH1TVu6rq/1bVnVX10PHxT6mqN1XVvVX1Nl+FgZN4clX9r6r6QFW9tKoekjzwqyjjV2ueV1U/N752X7Xx2M+sqtdX1fvG7by+qh698dz/WFXfXVU/leTDSZ5bVT+7OYiqem5VvXbbAKvqjqr68ap6f1W9o6r+/Ma/vaCqfnQcz4eq6q1V9cQbnvuacWzvrKq/esNzX11VPzA+9+1V9btvOOan7/jYL6qq/zb+24+M4/muW7kgcGbWEGMeMY7tnnH7vcd9/rive8cY8MyNf/uq8Tg/NM6lnrfvcY7//vyqes84lj9XG/MqTkfCxFZV9fgk35rkya2125M8I8ldSf5mkqckeVKSJyb5PUm+Y49Nf32S705ye5I7k/z9JF+c5PcleUSS5yf5RFU9Ksm/TvJd498/L8lrquqRU48NuNI3ZHi9/44kvzNXv77/aJKvTPJ5Sb4wyTeNf/+gJC9N8tgkn5vk/yX53hue+41J/kKGWPCSJJ9XVZ+/8e9/IsnLO/v94SS/kuSOJF+b5O9U1dM2/v1ZSX4kQ+x4RZLXVtWDa/iU+XVJ3pbkUUmeluQ5VfWMjec+M8krkzw8yY9vGXdu9tiq+tQk/yrJy8Yx/HCSr7liO3BJ1hBjXp7kYUmekOS3JvmHNz6gqh6cIZ782/ExfyXJD43zpyT5Z0m+eZxDfUGSn9z3OKvqK5P89SRPT/K4JF9xxTY4IgkTPfcluS3J76qqB7fW7mqt/VKGQPfC1tp7W2vvS/K3MwSlXf1Ya+2nWmufSPLRJH8myV9rrb27tXZfa+1NrbXfyBDI3tBae0Nr7ROttX+X5GeSfNUhDxJ4gO9trd3dWnt/hg83vu6Kx76ktXbP+NjXZfggJa21/9Nae01r7cOttQ+N27nxjf5lrbW3t9Y+Pr7mX5XhdZ+qekKS357k9TfusKoek+RLk3xba+0jrbX/nuT7c/849LOttR9trX0syfckeUiGD3qenOSRrbUXttY+2lr75ST/NMkf33junWPcuS/DpOmJ6es99ikZaoRf0lr7WGvtXyb5r1dsBy7J0mPM5yT5g0m+pbX2gfE1/MYtY3tKkk9L8qIxnvzkuL3rx/OxDHOoTx+389Z9jzNDIvXS8Tg+nGHOxQwkTGzVWntHkuckeUGS91bVK6vqjgyf6L5r46HvGv9uV3dv/PmzMkxkfmnL4x6b5I+My9z3VtW9GSZJn7PHvoD9bb5Gb/b6/tWNP384w+QhVfWwqvonNXzV9oNJ/lOSh9f96xY395Mk/yLJ11dVZUh+Xj1Ocm50R5L3j5OkzXE+atu2xw9nrq9GPTbJHTfElW9P8tuuOKaHVFWvQVLvsXckeXdrrV1xvHCplh5jHpMhxnzgJsdxR5K7xxizeTzXY9EfzvAh77uq6o1V9SVXbGvrcV7fx8a/iSMzkTDR1Vp7RWvtSzNMMlqSv5fknvH/r/vc8e+S5NczLGEnSarqs7dtduPPv5bkIxmW5W90d5KXt9YevvHfb2mtveiWDwjYxWM2/rz5+t7HczN0gfq9rbVPT/Ll49/XxmM2Y0Faa2/OsOr8ZRm+utv7qsw9SR5RVbffMM53b/z/J49h/Breo8fn3Z3knTfEldtba4deuX5PkkeNE7MHjAku3NJjzN0ZYszDbzKGe5I8pu7fUOKTsai19pbW2rMyfF3vtUlefbOD2uI9GeLXdeLITCRMbFVVj6+qP1BDQ4ePZPh+8H0Zvov/HVX1yKr6rCR/K8kPjk97W5InVNWTxoLFF1y1j/FTmX+e5HtqKMT+lKr6knGfP5jkq6vqGePfP2QslHz0VdsEJvvLVfXoqnpEhtWXV93CNm7PEDPuHbfznTs+7wcy1CF8vLV257YHtNbuTvKmJH93jAtfmOTPJvmhjYd9cVU9e1zteU6S30jy5gxfi/tgDQ1tHjrGli+oqiffwjFe5b9kiJffWlXXqupZGeo9geXHmPck+Ykk31dDc4kHV9WXb3noT2f4oPj542OemuSrk7yyqj61qr6hqj5j/GrwBzPEhH29OsmfHptLPCzDnIsZSJjouS3JizKsAv1qhk9Ivj1DE4afSfJzSf5HkreOf5fW2v9O8sIk/z7JL2Zo6nAzzxu385Yk78+wivWgcVL0rHGf78vwic/fiHsWju0VGYqYf3n871Y6u704yUMzxI83J/k3Oz7v5RmKo3uf/F73dRnqD+7J0FzhO8c6x+t+LMkfS/KBDF+9efZYh3BfhgnNk5K8cxzf9yf5jB3Ht5PW2keTPDtDIndvhrqJ12dI3ODSrSHGfGOGGqSfT/LeDB+83M/4On9mhnqnX0vyfUn+ZGvt5ze2cdf4lcFvyVg/tY/W2k9kaFjxH5K8I8OHMYlYcnJ1/69YA8A8avhJgfcm+aLW2i/e4jZekORxrbW9JyfHVFU/neQft9ZeOvdY4FIdIsbMaezy9z+T3NZa+/jc47kkPq0HYCn+YpK3rHEic6Oq+oqq+uzxK3l/KkOr4F0/BQeOY3Uxpqq+ZvyK32dm+BbO6yRLp9fr/AMAJ1NVd2Uo2P5DMw/lUB6fof7g0zJ0Av3asTYCmMGKY8w3Z/hNt/uSvDHJX5p1NBfKV/IAAAA6fCUPAACgQ8IEAADQceoaJt//uwD3/61GLl1r7dg3xE5xZQ335aG/Ir3tmKfsYw3nkMuwlLjCuk2JaWsoaTl0zN52zEs7hxPH032yFSYAAIAOCRMAAECHhAkAAKBDwgQAANAx+w/XnlMR8a7Fa4cukJvzHE4p2FtDweSudr0Ghy64n/McLvm1u+SxXWXKNd71mHfdx1rPIYMp983SCrFZlinvd8e4Dw49D1nae+0UU67Brse8tHN4rPvOChMAAECHhAkAAKBDwgQAANAhYQIAAOioExeyrbNqbotTFGevwSUWi5+iiHIN9jjmY98QO53Uc78vuVxrjSsTGyAtIq6cu7XGzbmacLFMh5ivWGECAADokDABAAB0SJgAAAA6JEwAAAAd1+YewJwFhZfYsGBXpzg3cxZMntMveR/aOZybc4or59RMRsw9jkM3ITrG6911XqdDNzk6xn0w5X7ddYznNF85pyZVp7xWVpgAAAA6JEwAAAAdEiYAAIAOCRMAAEBHnbioa/kVZEymuJdNV/1y9qF2cciNLa3gd4q1HosYws2sLa6wTKeINefU7OacjmWbq+KKFSYAAIAOCRMAAECHhAkAAKBDwgQAANBxbe4BACzJGpoi7Gqtx3LuhcXA5Tin2HVOx7IvK0wAAAAdEiYAAIAOCRMAAECHhAkAAKBj1U0fjvEr9sfY5pL3O/e+b9UxisLnKjSfs8Bdcf0DiSvz7/vS78G1E1c4FvOV421zyfude9+JFSYAAIAuCRMAAECHhAkAAKBDwgQAANAhYQIAAOhYTZe8XTuj7NNBZa6uPYfuHrVPh5hdj3nXx52iO82U69w7jrm66kzp8jK1Q8yux3xJnczOPa5MedyccYX1mBJX3A9cZUpcWfp7/FrnK1Met/b5ihUmAACADgkTAABAh4QJAACgQ8IEAADQsZqmD4cuAOs9VhHqYEnn5hgF5UsqEl2aSzo34gpMM7WgHXqmNuCZg/nK+bLCBAAA0CFhAgAA6JAwAQAAdEiYAAAAOurEhWOq1PZwqkK/XYsPT/mLypfuVEXTp7imrbVj3xDiyh7WEFfEEG5GXOEQxJrpLmW+YoUJAACgQ8IEAADQIWECAADokDABAAB0aPrAwSmiZJPibA5BXGGTuMIhiCts0vQBAADgFkiYAAAAOiRMAAAAHRImAACAjmtzDwCA87OtmPrETYYAOBPb3j9O2bTDChMAAECHhAkAAKBDwgQAANAhYQIAAOjQ9OHCTCmQU7C9blOun19D5yq73h+7NoJwv63brrHGdWZfu94z5ivrtuv1O+X7hxUmAACADgkTAABAh4QJAACgQ8IEAADQcRFNH6YWgM1VPLi04sZdC7aXZOr45ipKXlrRtML8B7rEuHKMa37p99HaTbmPxRWuMuVeMF85nSnzlWNck2PFFStMAAAAHRImAACADgkTAABAh4QJAACgY/amD3M2NlhSwenSGjycYt9zNjZYUkHo0ho8LG3ft2INr6dTOPR56D1ubfcHhzMllrpvuMqUOdqUmOae7jv0eeg9bsp787GugRUmAACADgkTAABAh4QJAACgQ8IEAADQMXvTh6VZehH4MZoV7Fr0uKRGCcew9GLNYzQp2fWaLqlByo2WMo6rLP11Mmdc4fy4zhzLKWLp0uN1svzX2DnOV6wwAQAAdEiYAAAAOiRMAAAAHRImAACAjkU2fZiz4G6uxganOuYpv359CnMWMs5VpH6qY576y9trJ67Mvx+NIIBzYb4y/35OeQ2sMAEAAHRImAAAADokTAAAAB0SJgAAgI5FNn04dIH0Gn61mcGhCxkVlJ+/Xe+ZU8UV99x+NIIAlmTX94pTzVfMYZfBChMAAECHhAkAAKBDwgQAANAhYQIAAOiYvenDKQp+FdL9pl3P7a7XZco5PPT2tlGY/5t2Pbe7Xpcln8M548olOnRcATiFXePPnPOVS7TE+YoVJgAAgA4JEwAAQIeECQAAoEPCBAAA0CFhAgAA6Ji9S54OSYMpXab2saQuLEsay5ymdIPZxyW91i7pWK8yZ1zZtm/XBViSU3QEPidzzle27fuU18UKEwAAQIeECQAAoEPCBAAA0CFhAgAA6Ji96QN9xyhGVMy4Drtep30KK6c+f6nO4Rjmdoy4ArBG5kTTHWO+MjcrTAAAAB0SJgAAgA4JEwAAQIeECQAAoOMimj70iscUKg/OuRFE7xqfy/FNdQ6NIOY6hrXGlTmv+Tncb8B5m2tOtNb5ypxzyFPu2woTAABAh4QJAACgQ8IEAADQIWECAADoqBMXky27co2DUMTNptbasW8IceXMiCHcjLjCIYg1bLoqrlhhAgAA6JAwAQAAdEiYAAAAOiRMAAAAHdfmHgAAbNrWjEhxNgBzscIEAADQIWECAADokDABAAB0SJgAAAA6JEwAAAAdEiYAAIAOCRMAAECHhAkAAKBDwgQAANAhYQIAAOiQMAEAAHRImAAAADokTAAAAB0SJgAAgA4JEwAAQIeECQAAoEPCBAAA0CFhAgAA6JAwAQAAdFybewAAsKmq5h4CAHySFSYAAIAOCRMAAECHhAkAAKBDwgQAANAxe9OHcyruba3t9LhLPOZDP3dpdr2ml3jMc1jy2PZ1iXGF6dYQa9yz67SG67br/b/tcduObw2vp12Zr9waK0wAAAAdEiYAAIAOCRMAAECHhAkAAKBj9qYP26yh0GwNRY9rtIbzuob785Ks4Z6Z4tyPj+nWGpN2LbhnPmu9HlMaN6z19bSrcz++Y7HCBAAA0CFhAgAA6JAwAQAAdEiYAAAAOhbZ9GGtRYZMpxiRfZ174fiur4lzOmb2M6XAfU7u2eVba3ydcv+v9fW0q12v3zkd8yFYYQIAAOiQMAEAAHRImAAAADokTAAAAB114qIuFWQXYA0FoZxOa+3YN8RB48oxCn7nKiKes7h3yr7FEG5mbXGFZTpFrDlGjJsrbs7ZhOgU+74qrlhhAgAA6JAwAQAAdEiYAAAAOiRMAAAAHRImAACADl3yODgdrtikmxWHIK6wSVzhEMQVNumSBwAAcAskTAAAAB0SJgAAgA4JEwAAQMe1uQcwxbZivalNLI6xzVvd7zbHGMuc+75V28YytXjzGNu81f1uc4yxzLnvpRJX5t/3Od1v535824grpzFXXJmT+crxtnmr+93mHOcrVpgAAAA6JEwAAAAdEiYAAIAOCRMAAEBHnbgw7pZ3tmtx4z5FkFMKw6actymFmlOLPA9dDHeKwsMp17k3lkNfv11NOV9Tz/WhX+t7jPvYlagHjSs773QFcWWKOePKWi3p+p2KuHIac8WVOZ3LfGVXa5ivTLH2+YoVJgAAgA4JEwAAQIeECQAAoEPCBAAA0HFt7gHsamqDh1Ntk8Ob2uDhVNtkfcQVDklc4VjElXUwXzlfVpgAAAA6JEwAAAAdEiYAAIAOCRMAAEBHnbhAUDVixx6/Qrz4fSs83M+c5/UU+77ql7MPRFzpWGtcEUPWTVzhWMxX5nXJ8xUrTAAAAB0SJgAAgA4JEwAAQIeECQAAoEPTBw5OESWbFGdzCOIKm8QVDkFcYZOmDwAAALdAwgQAANAhYQIAAOiQMAEAAHRcm3sAAJyfXYuptzUeUogNwKZdm9Qd6/3DChMAAECHhAkAAKBDwgQAANAhYQIAAOjQ9OHCHLoYbtciPOZ36GulMJ/rptwL7qPzM3dxNudrSjMZ1mPK9TtWIyErTAAAAB0SJgAAgA4JEwAAQIeECQAAoOMimj5MLfY6RfHglDEqbuybem5OUZQ8ZYyKpudziXHF/bZux7imSyzO5vS2Xbc533+PMZ5ju8T5ytKvySYrTAAAAB0SJgAAgA4JEwAAQIeECQAAoGP2pg9z/mqzgtN5zflr8GsvPuRqfg1+OrFw3Xa9t70Psq9d7481vM/uev/POV85F0u79vuywgQAANAhYQIAAOiQMAEAAHRImAAAADpmb/oAMIUi2+PYp0DXNVieXQvXp147BfKX59CNEqbuZ9fnTtkv0+3zel/iNbDCBAAA0CFhAgAA6JAwAQAAdEiYAAAAOjR9WIhTFbgtsZDu0p2q8FmB9eWZM66435ZHrOFYzGEuw5wxZO5mMlaYAAAAOiRMAAAAHRImAACADgkTAABAh4QJAACgQ8IEAADQIWECAADokDABAAB0SJgAAAA6JEwAAAAd1+YeAKdVVQfdXmvtoNvjeA59rQ59L7Fe7gUADmXKfGXbcw/xHmWFCQAAoEPCBAAA0CFhAgAA6JAwAQAAdGj6sBBTCtKmFvPv+nyF3ccx5fpNvSa7Pl9zj3U6dFwRAwAu16HnK2uaW1hhAgAA6JAwAQAAdEiYAAAAOiRMAAAAHZo+LJii68u19uJI1k+sAeBQ1j6HscIEAADQIWECAADokDABAAB0SJgAAAA6Lrrpw9KLmk81vqWfh2NYevHhqca39POwC81RjmOfe8P5Bs6FpkvHsc/7xBLPtxUmAACADgkTAABAh4QJAACgQ8IEAADQUScurFpeFRcHpwCcTa21Y98Qi4orSy8YPtX4xAGO6dLiCsexhji19MZGpxrfru9TU/Z9VVyxwgQAANAhYQIAAOiQMAEAAHRImAAAADpO3fQBAABgNawwAQAAdEiYAAAAOiRMAAAAHRImAACADgkTAABAh4QJAACgQ8IEAADQIWECAADokDABAAB0SJgAAAA6JEwAAAAdEiYAAIAOCRMAAECHhAkAAKBDwgQAANAhYQIAAOiQMAEAAHRImAAAADokTAAAAB0SJgAAgA4JEwAAQIeECQAAoEPCBAAA0PH/ATS+5rIZbk+MAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 864x288 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot3(img0, binary_opening(img0, disk(1)), binary_closing(img0, disk(1)), 'source', 'binary opening', 'binary closing')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0kAAAEtCAYAAADdgyslAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAEKBJREFUeJzt3XusJnddx/HPF7a0RQrdBhR6AdkSGywBwrZaDDeFpEgCFWSrgCjeWlRUIpUSJFIJ6JooEkJQI1qkXLugXCoY26BoQZAuCooBgaVN6ZYU7FYaKrfy84+ZxW+XvbNn5+ye1yvZ5Ow588xvZs95fjvvZ2aeU2OMAAAAMLnT0hsAAACwmogkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSwBpUVddW1eP28LVHVtUnD/c2rSZV9fGqeszS2wFHuiN9rqmqx1TV577Ddaz6/eTbrVt6AwBYXcYY/5TkjKW3Y0ljjDOX3gY42q2VuWat7OfRxpkkVlxV3XnpbQCWVxP/7wArylzDoeAHiFTVxVV1Q1XdWlWfrKrHVtWxVfWKqto+/3lFVR07L/+sqrp6l3WMqnrA/PFrq+qPq+rdVfXlJD9cVcdX1R9W1XVV9T9VdXVVHT8vf05VfaCqbqmqj7rEBQ6bs6vqP6tqR1VdWlXHJd9+ecl8ucxFVfWx+fn7lrbs+qq6oqq+MK/niqo6tT32H6rqZVX1/iS3JXleVW3tG1FVz6uqt+9uA6vq5Kp6Z1XdXFWfrqpfbF+7pKreOm/PrVX1kap6yC6Pfdu8bZ+tql/b5bGXV9Xr5sd+vKrO2mWfH7efyz6sqv51/tqWeXteejDfEDhKHQlzzUnztm2f17+n5R44j3XLPBc8qX3tCfN+3jofV110oPs5f/35VXXjvC2/UO0Yi8NHJK1xVXVGkuckOXuMcUKSc5Ncm+S3kpyT5KFJHpLkB5K86ABW/fQkL0tyQpKrk/xBko1JfijJSUmen+SbVXVKkr9J8tL58xcleVtV3es73Tdgn56R6Tl/epLvy96f4+cneXyS+yd5cJJnzZ+/U5JLk9wvyX2T/G+SV+3y2GcmuSDTfPDKJPevqge2r/9Uksv2MO6bknwuyclJnprkd6vqse3r5yXZkmn+eGOSt1fVMTW9ivyuJB9NckqSxyZ5blWd2x77pCRvTnJiknfuZruzr2Wr6i5J/jrJa+dteFOSJ+9lPbAWHQlzzWVJ7prkzCTfneSPdl2gqo7JNK/83bzMryZ5w3wslSR/nuTC+XjqQUnee6D7WVWPT/IbSR6X5AFJHr2XdbCCRBK3Jzk2yfdX1TFjjGvHGJ/JNKG9ZIxx0xjjC0l+J9Pks7/eMcZ4/xjjm0m+luTnkvz6GOOGMcbtY4wPjDG+mmnCevcY491jjG+OMa5Mck2SJxzKnQR261VjjOvHGDdnelHjaXtZ9pVjjO3zsu/K9AJKxhj/PcZ42xjjtjHGrfN6dv1P/bVjjI+PMb4xP+/fkum5n6o6M8n3Jrli1wGr6rQkj0hy8RjjK2OMf0vymtxxLto6xnjrGOPrSV6e5LhML/CcneReY4yXjDG+NsbYluTPkvxke+zV89xze6YDpIdkz/a07DmZ7u995Rjj62OMv0ryL3tZD6xFq32uuU+SH03y7DHGjvm5/L7dbNs5Se6WZPM8r7x3Xt/O/fl6puOpu8/r+ciB7memeLp03o/bMh1/sQCRtMaNMT6d5LlJLklyU1W9uapOzvSq7XVt0evmz+2v69vH98x04PKZ3Sx3vySb5tPWt1TVLZkOiu5zAGMBB6c/T/f1HP98+/i2TAcKqaq7VtWf1nQp7ZeS/GOSE+uO9yL2cZLkL5M8vaoqU/BcPh/Q7OrkJDfPB0R9O0/Z3brnF2V2nnW6X5KTd5lbXpjke/ayT8dV1Z7e0GhPy56c5IYxxtjL/sJat9rnmtMyzTU79rEfJye5fp5r+v7snJN+PNOLvNdV1fuq6uF7Wddu93PnGO1r5pOFiCQyxnjjGOMRmQ4qRpLfT7J9/vtO950/lyRfznRKOklSVffe3Wrbx19M8pVMp9l3dX2Sy8YYJ7Y/3zXG2HzQOwTsr9Pax/05fiCel+ldm35wjHH3JI+aP19tmT4fZIzxwUxnmB+Z6dLcPV3+sj3JSVV1wi7beUP7+7f2Yb7E7tT5cdcn+ewuc8sJY4xDfZb6xiSnzAdh37ZNQJLVP9dcn2muOXEf27A9yWl1xzeF+NacNMb48BjjvEyX4r09yeX72qnduDHTPLaT+WQhImmNq6ozqupHanpThq9kusb39kzX1b+oqu5VVfdM8ttJXj8/7KNJzqyqh843Gl6ytzHmV1z+IsnLa7qR+s5V9fB5zNcneWJVnTt//rj5BsdT97ZO4JD4lao6tapOynSW5S0HsY4TMs0bt8zrefF+Pu51me4n+MYY4+rdLTDGuD7JB5L83jw3PDjJzyd5Q1tsY1U9ZT6r89wkX03ywUyXvH2ppjemOX6eXx5UVWcfxD7uzT9nmjOfU1Xrquq8TPdwAv9vtc81NyZ5T5JX1/QGEcdU1aN2s+iHMr1Q/Px5mcckeWKSN1fVXarqGVV1j/ny3y9lmhsO1OVJfnZ+g4i7Zjr+YgEiiWOTbM50tufzmV79eGGmN1K4JsnHkvx7ko/Mn8sY47+SvCTJVUk+lemNGfblonk9H05yc6azVXeaD4LOm8f8QqZXc34zfjbhcHhjphuQt81/DuYd2V6R5PhMc8gHk/ztfj7uskw3Nu/pld2dnpbpPoLtmd4g4cXzvYs7vSPJTyTZkelymqfM9xPcnung5aFJPjtv32uS3GM/t2+/jDG+luQpmeLtlkz3P1yRKdaAyZEw1zwz0z1Fn0hyU6YXXe5gfr4/KdP9S19M8uokPz3G+ERbx7Xz5YDPznw/1IEYY7wn05tO/H2ST2d6ISYxpxx2dcfLqAFg5dX0KwBuSvKwMcanDnIdlyR5wBjjgA9EVlJVfSjJn4wxLl16W2CtOxRzzZLmd+f7jyTHjjG+sfT2rCVerQdgCb+U5MNH4kHLrqrq0VV17/lyu5/J9Ha++/sqN7Cyjri5pqqePF++tz7TlTfvEkiH357exQcAVkRVXZvpZusfW3hTDpUzMt1HcLdM7+L51PkeB2BBR/Bcc2Gm3712e5L3JfnlRbdmjXK5HQAAQONyOwAAgEYkAQAANCtyT1JVrdlr+DZt2rTo+Js3L/s7WK+66qpFx3/BC16w2Ng7duzrF3Uf3cYYte+lvnNreX5xefTadcffVbv2HI75xdzCWmRu2fPc4kwSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAADNuqU34GizefPmRcffsGHDouOvX79+0fFvvvnmxcY+//zzFxs7SbZs2bLo+MDKGWMsOn5VLTo+sDLMLXvmTBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKBZt/QGrISNGzcuNvaGDRsWGztJTj/99EXH37Zt26LjX3nllYuNveTPXZJs2bJl0fHXijHG0puwZlXVouMv+b1fet9ZeeYW1qLV/HPvTBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAM26pTdgJaxfv36xsbdu3brY2Emybdu2Rcdf2tL//hz9qmqxsccYi42dLLvvq8Fa3/+j3dLPL2B1cSYJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBm3dIbsBLWr1+/2NhXXXXVYmOz7Pd+x44di43N2lBVS2/CmjbGWGxs33uAw8uZJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAmnVLb8BK2LFjx2Jjb9y4cbGxV4P169cvOv6S//5btmxZbGzWhjHGouNX1aLjL73/S1p635f+3gNHp6Xnlr3Nrc4kAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAmhpjHPqVVh36lR6ADRs2LDb2Nddcs9jYSXLhhRcuOv6mTZsWHX/J7/1ZZ5212NirwRijDsc4S88vKzFnwr5UHZan16p1OOYXcwtrkbllz3OLM0kAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQ1xjj0K6069Cs9QlxwwQWLjn/xxRcvOv7WrVsXHf/8889fdPy1bIxRh2OctTy/rMR8zf6rOiw/4uzG4ZhfzC0sxdyynL3NLc4kAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQ1Bhj6W0AAABYNZxJAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKD5P4sz4iT7BZreAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 864x288 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot3(img1, binary_opening(img1, disk(1)), binary_closing(img1, disk(1)), \\\n",
    "      'source', 'binary opening', 'binary closing')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0kAAAEtCAYAAADdgyslAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAD8tJREFUeJzt3X2MpWdZx/HfBa0FKQnLm9paMFsibxqJsyBUElERQgXkD3cRKoYYnaJiMBEoAvEFNdk1SpAowpJoE17jCsZIm0g3IqYSeRkUDBGErpSWQqAwJVApBHr7x3kqF2t3acvs3DudzydpMnvOmXM9Z7rn3vt7zjMzNcYIAAAAK3eZfQAAAACnE5EEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCYEerqi9V1d7ZxwHMVVXPrqorv837uKiq3r5Vx8TOdcbsAwCAb8cY4+zZxwDcOYwx3pDkDbOPg/m8k8QpV1V3nX0MAABwW4kkUlWXVNUnq+qLVfWRqvrJqjqrql5RVdct/72iqs5abv//3s6uqlFVD1o+vrSq/qKqLq+qG5P8eFXdvar+pKqurqovVNWVVXX35faPrqp3VdUNVfWBqnrcdn8NgK1RVQ+tqn9ans8fqqqnLpdfWlWvrqorlrXmnVX1wPZ5D1mu+/yyDh1o111aVX9eVZctn/vuqjq/XX/8+nOy2z5huf8vVNWrluP4pe356gBbparOq6q3VtVnq+pzVfVnt3KbC6rqvcvz/b1VdUG77tlVdWxZJ/67qi5ql1/Zbjeq6jlV9dGq2lzWl1quu+uyt7l+uY/nLrd3ptadgEja5arqwUmem+SRY4x7Jnliko8neUmSRyd5RJIfSvKoJC+9HXf9zCR/mOSeSa5M8sdJ1pJckOTeSV6Y5OaqOjfJZUn+YLn8+UneUlX3+3YfG7C9qurMJH+f5O1J7p/k15O8YVlnkuSiJL+f5L5J/j3LKS1VdY8kVyR54/J5z0jyqqp6eLv7ZyT5vSR7knwsq/XlRG71tlV13yR/k+S3ktwnyUeyWpOAHWQ5Q+VtSa5O8n1Jzk3y5uNuc++s9hevzOr5/vIkl1XVfZY155VJnrTsfS7Iak06kScneWRW+6EDWe2VkuSXkzwpq73SDyd52hY8PE4TIomvJzkrycOq6swxxsfHGFdltZl52RjjM2OMz2a14XjW7bjfvxtj/MsY4+YkX03yi0meN8b45Bjj62OMd40xvpLk55NcPsa4fIxx8xjjiiTvS3LhVj5IYFs8OsnZSQ6OMb46xvjHrDYyz1iuv2yM8c/Lc/8lSR5TVedltQH5+Bjjr8YYXxtjvD/JW5L8bLvvt44x3jPG+FpWcfWIkxzHiW57YZIPjTHeulz3yiSf3pJHDmynRyU5J8kLxhg3jjFuGmMc/wMbfjrJR8cYr1vWlTcl+XCSpyzX35zkB6rq7mOMT40xPnSSeQfHGDeMMT6R5B35xppyIMmfjjGuHWNsJjm4VQ+Q+UTSLjfG+FiS30jyu0k+U1Vvrqpzslp8rm43vXq57La6pn183yR3S3LVrdzugUn2L6fm3FBVNyR5bJLvuR2zgNPDOUmuWV4cucXVWb3Km7R1YYzxpSSfXz7ngUl+5Lh14KIk393up8fM/2QVYydyotuec9wxjCTX3obHBZxezkty9fJix4kcv4/J8udzxxg3Jnl6kuck+dRyeu5DTnJft2lNOe5jdjiRRMYYbxxjPDarjcpIcijJdcufb/GA5bIkuTHJd95yRVX1jcz/3W37+PokNyU5/1Zud02S140x7tX+u8cYw6sxsPNcl+S8qur/tjwgySeXj8+75cKqOjurU2yvy2odeOdx68DZY4xf2eLj+1SS723HUP3PwI5xTZIHfIvv/Tl+H5O09WiM8Q9jjJ/K6kXZDyd57R04jm9aU9LWOHY+kbTLVdWDq+onlh/KcFOSL2d1Ct6bkry0qu63nMf/20lev3zaB5I8vKoeUVV3y+pdqBNaXlX+yyQvr6pzlm90fMwy8/VJnlJVT1wuv1tVPa6qbFxg53l3Vi+ivLCqzlx+CMtT8o3vFbiwqh5bVd+R1fcmvXuMcU1Wp+R9f1U9a/m8M6vqkVX10C0+vsuS/GBVPW3ZXP1avvndKmBneE9WgXKwqu6x7B1+9LjbXJ7VuvLMqjqjqp6e5GFJ3lZV31VVT12+N+krSb6U1d7n9vrrJM+rqnOr6l5JLrnjD4nTjUjirKzOob0+q7eT75/kxVn9IIX3Jflgkv9I8v7lsowx/ivJy5IcTfLRrH4ww7fy/OV+3pvVKTaHktxl2SD9zDLzs1m9OvSC+LsJO84Y46tJnprVNzJfn+RVSX5hjPHh5SZvTPI7Wa0Ba1mdUpcxxheTPCHJz2X16u+ns1ojztri47s+yf4kf5Tkc1ltmN6X1SYJ2CHGGF/P6gWYByX5RFanzT79uNt8Lqvvd/zNrJ7vL0zy5GUduMty+XVZrUc/luRX78ChvDarH1TzwST/llWYfS13LLg4zdTqlGwAOHWq6tIk144xbs9PyTylltMCr01y0RjjHbOPB9jZqupJSV49xjj+ND92IK/WA7BrLKf23ms53ffFSSrJv04+LGAHqtXvgLxwOZ3v3KzeKf/b2cfF1hBJAOwmj8nqJ21en9XpOk8bY3x57iEBO1Rl9StSNrM63e4/s/oebu4EnG4HAADQeCcJAACgEUkAAADNyX4J1x1WVbv2HL79+/dPnX/w4NzfwXr06NGp81/0ohdNm725uTlt9ulgjFHbMWc3ry9Oj55r375902ZvbGxMm3062I71ZTevLbt97zLbzL3LkSNHps0+HZxsbfFOEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAATY0xtv5Oq7b+TneIq666aur8vXv3Tp1/5MiRqfP3798/bfaBAwemzU7mf+3HGLUdc3bz+rK+vj51/mte85qp8y+++OKp8zc2NqbOn2n2Y9+O9WU3ry27fe9y9OjRqfMf//jHT5tt73LitcU7SQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgOaM2QdwKqytrU2bvXfv3mmzk+T888+fOv/YsWNT519xxRXTZs/8e5ckR44cmTp/t1hfX599CNNcfPHFU+cfPnx46vyZz/GNjY1ps9keu3nvcvTo0anzd/Pza/be5XTmnSQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAAJozZh/AqbBnz55pszc2NqbNTpJjx45NnT/b7K8/d36HDx+eNnt9fX3a7GTuYz8dWF/u3NbW1qbO3817l9nzZ9vNj3/28+5kvJMEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEBzxuwDOBX27NkzbfbRo0enzWbu//vNzc1ps9kdDh8+PPsQdrW1tbVpszc2NqbN3i1m/vsxe769y1y7ee8ye/7JeCcJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIDmjNkHcCpsbm5Om722tjZt9ulgz549U+fP/PofOXJk2mx2h/X19anzDx8+PHX+bl5fZz/2jY2NqfO3w7Fjx6bO37t377TZu/3vl73LPIcOHZo6/+DBgye8zjtJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACA5ozZB3AqHDt2bNrstbW1abOTZP/+/bt6/kyHDh2afQhsg/X19dmHMM0YY+r8ffv2TZ0/08bGxuxDuNObuXeYbfbeZbY9e/bMPoRp7F1OzDtJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0NcbY+jut2vo73SHW19enzr/kkkumzt/Y2Jg6/8CBA1Pn72ZjjNqOObt5fTkV6zW3XdW2/BXnVmzH+rKb15bdvnfZ3NycOn/fvn1T5+9mJ1tbvJMEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEBTY4zZxwAAAHDa8E4SAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAM3/Ar9DCEmqQvMWAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 864x288 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot3(img1, opening(img1, disk(1)), closing(img1, disk(1)), \\\n",
    "      'source', 'opening', 'closing')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0kAAAEtCAYAAADdgyslAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAD+VJREFUeJzt3XvsX3ddx/HXGzq2IWNlAYVuA4nEBUeggNNhgHVCApJABfECiOIVVNRFJiRIXCWgmCguhKBGdJNxHaAU5jBqsNOBIBcLigG5bSnryMC1sjC5jY9/nFN877fe+bWn7e/xSJr8+vud7zmfk+376Xmey/dXY4wAAAAwudPSAwAAADiWiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBLAGlRV11XVY/fxs0dV1ceP9piOJVX10aratPQ44Hh3vM81VbWpqj77La7jmN9P7mjd0gMA4NgyxvjnJOcsPY4ljTHOXXoMcKJbK3PNWtnPE40rSRxxVXXnpccALK8m/t0BjihzDavB/0Ckql5QVTdU1S1V9fGqekxVnVxVl1bVzvnPpVV18rz8s6rq2hXrGFX1gPnry6vqj6vq6qr6UpILq+rUqvrDqrq+qv6nqq6tqlPn5c+vqvdU1e6q+rBbXOCoOa+q/rOqdlXVZVV1SnLH20vm22UurqqPzO/fN7Vl71FVV1XV5+f1XFVVZ7XXbquql1bVu5PcmuR5VfXBPoiqel5VvW1vA6yqDVX19qq6uao+WVW/0H62pareMo/nlqr6UFU9ZMVr3zqP7TNV9WsrXntlVb1mfu1Hq+p7V+zzYw9y2YdV1b/NP3vzPJ6XHM5/EDhBHQ9zzRnz2HbO69/Xcg+ct7V7ngue1H72hHk/b5mPqy4+1P2cf/78qrpxHsvPVzvG4ugRSWtcVZ2T5LlJzhtjnJbkcUmuS/JbSc5PsjHJQ5J8X5IXHcKqn57kpUlOS3Jtkj9I8vAkP5DkjCTPT/KNqjozyd8kecn8/YuTvLWq7vWt7htwQM/I9J7/riTfnf2/x38syeOT3D/Jg5M8a/7+nZJcluR+Se6b5H+TvHLFa5+Z5BczzQevSHL/qnpg+/lPJrliH9t9Q5LPJtmQ5KlJfreqHtN+vjnJmzPNH69P8raqOqmms8jvSPLhJGcmeUySi6rqce21T0ryxiTrk7x9L+POgZatqrsk+eskl89jeEOSJ+9nPbAWHQ9zzRVJ7prk3CTfnuSPVi5QVSdlmlf+bl7mV5O8bj6WSpI/T/Ls+XjqQUnedaj7WVWPT/IbSR6b5AFJLtjPOjiCRBK3JTk5yfdU1UljjOvGGJ/KNKG9eIxx0xjj80l+J9Pkc7C2jjHePcb4RpKvJvnZJL8+xrhhjHHbGOM9Y4yvZJqwrh5jXD3G+MYY4++TfCDJE1ZzJ4G9euUYY8cY4+ZMJzWetp9lXzHG2Dkv+45MJ1AyxvjvMcZbxxi3jjFumdez8h/1y8cYHx1jfH1+378p03s/VXVuku9MctXKDVbV2UkemeQFY4wvjzG2J3l1bj8XfXCM8ZYxxteSvDzJKZlO8JyX5F5jjBePMb46xvh0kj9L8hPttdfOc89tmQ6QHpJ929ey52d6vvcVY4yvjTH+Ksm/7mc9sBYd63PNfZL8UJLnjDF2ze/la/YytvOT3C3Jy+Z55V3z+vbsz9cyHU/dfV7Phw51PzPF02Xzftya6fiLBYikNW6M8ckkFyXZkuSmqnpjVW3IdNb2+rbo9fP3DtaO9vU9Mx24fGovy90vyY/Ol613V9XuTAdF9zmEbQGHp79PD/Qe/1z7+tZMBwqpqrtW1Z/WdCvtF5P8U5L1dftnEft2kuQvkzy9qipT8Fw5H9CstCHJzfMBUR/nmXtb93xSZs9Vp/sl2bBibnlhku/Yzz6dUlX7+kCjfS27IckNY4yxn/2Fte5Yn2vOzjTX7DrAfmxIsmOea/r+7JmTfiTTSd7rq+qaqnrEfta11/3cs432M/PJQkQSGWO8fozxyEwHFSPJ7yfZOf99j/vO30uSL2W6JJ0kqap772217esvJPlypsvsK+1IcsUYY337821jjJcd9g4BB+vs9nV/jx+K52X61KbvH2PcPcmj5+9XW6bPBxljvDfTFeZHZbo1d1+3v+xMckZVnbZinDe0v39zH+Zb7M6aX7cjyWdWzC2njTFW+yr1jUnOnA/C7jAmIMmxP9fsyDTXrD/AGHYmObtu/6EQ35yTxhjvH2NsznQr3tuSXHmgndqLGzPNY3uYTxYikta4qjqnqn6wpg9l+HKme3xvy3Rf/Yuq6l5Vdc8kv53ktfPLPpzk3KraOD9ouGV/25jPuPxFkpfX9CD1navqEfM2X5vkiVX1uPn7p8wPOJ61v3UCq+JXquqsqjoj01WWNx3GOk7LNG/sntdzyUG+7jWZnif4+hjj2r0tMMbYkeQ9SX5vnhsenOTnkryuLfbwqnrKfFXnoiRfSfLeTLe8fbGmD6Y5dZ5fHlRV5x3GPu7Pv2SaM59bVeuqanOmZziB/3eszzU3JnlnklfV9AERJ1XVo/ey6PsynSh+/rzMpiRPTPLGqrpLVT2jqk6fb//9Yqa54VBdmeRn5g+IuGum4y8WIJI4OcnLMl3t+Vymsx8vzPRBCh9I8pEk/57kQ/P3Msb4ryQvTvIPST6R6YMZDuTieT3vT3JzpqtVd5oPgjbP2/x8prM5vxn/b8LR8PpMDyB/ev5zOJ/IdmmSUzPNIe9N8rcH+borMj3YvK8zu3s8LdNzBDszfUDCJfOzi3tsTfLjSXZlup3mKfPzBLdlOnjZmOQz8/heneT0gxzfQRljfDXJUzLF2+5Mzz9clSnWgMnxMNc8M9MzRR9LclOmky63M7/fn5Tp+aUvJHlVkp8aY3ysreO6+XbA52R+HupQjDHemelDJ/4xyScznYhJzClHXd3+NmoAOPJq+hUANyV52BjjE4e5ji1JHjDGOOQDkSOpqt6X5E/GGJctPRZY61ZjrlnS/Ol8/5Hk5DHG15cez1ribD0AS/ilJO8/Hg9aVqqqC6rq3vPtdj+d6eN8D/YsN3BkHXdzTVU9eb597x6Z7rx5h0A6+vb1KT4AcERU1XWZHrb+4YWHslrOyfQcwd0yfYrnU+dnHIAFHcdzzbMz/e6125Jck+SXFx3NGuV2OwAAgMbtdgAAAI1IAgAAaI7IM0lV5R6+hWzbtm3R7a9ff6Dfw3ZkXXLJwf7ahNW3devWxbZ9LBhj1IGX+taZX5bj9uzl3P531a49R2N+Mbcsx9yyHHPLvucWV5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGjWLT0AVtfu3bsX3f4FF1yw6PYvvPDCxba9devWxbYNAMDqcSUJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQrFt6ACeajRs3Lrr9TZs2Lbr9pW3fvn3pIcARM8ZYegjACcjcAnfkShIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAM26pQdwJFx00UWLbXvLli2LbTtJTj/99EW3v7Rt27YtPQROcGOMpYcAnIDMLXBscSUJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAICmxhirv9Kq1V/pcWL9+vWLbn/Xrl2Lbn9pD33oQxfb9vbt2xfb9rFgjFFHYztreX45EvM1x4eqo/L2OmYdjfnF3MJaZG7Z99ziShIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAM26pQcAq2njxo2LbXv79u2LbRsAgNXjShIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgWbf0AE40u3fvXnT7W7duXXT7mzdvXnT7mzZtWmzbl19++WLbZm2oqkW3P8ZYdPvAkWFugTtyJQkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgKbGGEuPAQAA4JjhShIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAzf8B2GKl75blp6YAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 864x288 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot3(img2, binary_opening(img2), binary_closing(img2), \\\n",
    "      'source', 'binary opening', 'binary closing')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0kAAAEtCAYAAADdgyslAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAADz1JREFUeJzt3X/M7nVdx/HXWw8hifMc/FFBYIbLX7XOOtOU3IA66SAV/whMycZagZXNs+WPpS7IasNWQq7Q4VbHKeggzbOELT1LbORSwI42FoYSyA+d4uE4JdEJn/74fsm3J84R8D73977PeTy2s933dX3v6/O5zrg/fJ7X93tdp8YYAQAAYPKIpScAAACwlogkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAOtaVX2jqn5y6XkAy6qqs6rq6h/wMc6sqg+v1JxYvzYsPQEA+EGMMY5ceg7AwWGMcUmSS5aeB8tzJokDrqoeufQcAADgwRJJpKpeX1W3V9XXq+qzVfVLVXV4VV1YVXfMfy6sqsPn4//f6eyqGlX1lPnr7VX19qq6sqruTnJyVR1RVX9ZVbdU1deq6uqqOmI+/jlV9fGq2lNVn66qk1b77wBYGVX19Kq6av59vr6qXjzfvr2q3lFVH5nXmo9V1ZPazz1tvm/3vA6d0e7bXlV/U1VXzD/7iao6vt2/9/qzv2OfPz/+16rqonkev7U6fzvASqmqY6vqA1X1lar6alX99QMcc0JVXTP/vl9TVSe0+86qqpvmdeK/q+rMdvvV7bhRVa+sqhur6q55fan5vkfOe5s758d41Xy8K7UOAiLpEFdVT03yqiTPGmM8JskLktyc5I1JnpNkc5KfTfLsJG96CA/98iR/luQxSa5O8hdJtiQ5IclRSV6X5L6qOibJFUn+dL79NUneX1VP+EGfG7C6quqwJP+Y5MNJnpjk95NcMq8zSXJmkj9J8vgkuzJf0lJVj07ykSSXzj/3siQXVdUz28O/LMkfJ9mU5HOZ1pd9ecBjq+rxSf4+yR8meVySz2Zak4B1ZL5C5UNJbknyE0mOSfK+vY45KtP+4m2Zft/fmuSKqnrcvOa8Lckp897nhExr0r68MMmzMu2Hzsi0V0qS305ySqa90s8leckKPD3WCJHEvUkOT/KMqjpsjHHzGOPzmTYzbx5jfHmM8ZVMG45XPITH3THG+Ncxxn1Jvp3kN5O8eoxx+xjj3jHGx8cY30ry60muHGNcOca4b4zxkSTXJjl1JZ8ksCqek+TIJOePMb49xvjnTBuZl833XzHG+Jf5d/+NSZ5bVcdm2oDcPMb4uzHGd8YYn0ry/iS/2h77A2OMT44xvpMprjbvZx77OvbUJNePMT4w3/e2JF9akWcOrKZnJzk6yWvHGHePMe4ZY+z9gQ2/kuTGMca753XlvUluSPKi+f77kvx0VR0xxvjiGOP6/Yx3/hhjzxjjC0k+mu+uKWck+asxxm1jjLuSnL9ST5DliaRD3Bjjc0m2JTkvyZer6n1VdXSmxeeWdugt820P1q3t68cneVSSzz/AcU9Kcvp8ac6eqtqT5HlJfuwhjAWsDUcnuXV+ceR+t2R6lTdp68IY4xtJds8/86QkP7/XOnBmkh9tj9Nj5n8yxdi+7OvYo/eaw0hy24N4XsDacmySW+YXO/Zl731M5u+PGWPcneSlSV6Z5Ivz5blP289jPag1Za+vWedEEhljXDrGeF6mjcpI8pYkd8zf3++4+bYkuTvJD99/R1X1jcz/PWz7+s4k9yQ5/gGOuzXJu8cYG9ufR48xvBoD688dSY6tqv7/luOS3D5/fez9N1bVkZkusb0j0zrwsb3WgSPHGL+zwvP7YpIfb3Oo/j2wbtya5Ljv896fvfcxSVuPxhj/NMb45Uwvyt6Q5J0PYx7fs6akrXGsfyLpEFdVT62qX5w/lOGeJN/MdAnee5O8qaqeMF/H/0dJ3jP/2KeTPLOqNlfVozKdhdqn+VXlv03y1qo6en6j43PnMd+T5EVV9YL59kdV1UlVZeMC688nMr2I8rqqOmz+EJYX5bvvFTi1qp5XVT+U6b1Jnxhj3JrpkryfqqpXzD93WFU9q6qevsLzuyLJz1TVS+bN1e/le89WAevDJzMFyvlV9eh57/ALex1zZaZ15eVVtaGqXprkGUk+VFU/UlUvnt+b9K0k38i093moLkvy6qo6pqo2Jnn9w39KrDUiicMzXUN7Z6bTyU9M8oZMH6RwbZLPJPmPJJ+ab8sY47+SvDnJziQ3Zvpghu/nNfPjXJPpEpu3JHnEvEE6bR7zK5leHXpt/LcJ684Y49tJXpzpjcx3JrkoyW+MMW6YD7k0ybmZ1oAtmS6pyxjj60men+TXMr36+6VMa8ThKzy/O5OcnuTPk3w104bp2kybJGCdGGPcm+kFmKck+UKmy2ZfutcxX830fsc/yPT7/rokL5zXgUfMt9+RaT06McnvPoypvDPTB9V8Jsm/Zwqz7+ThBRdrTE2XZAPAgVNV25PcNsZ4KJ+SeUDNlwXeluTMMcZHl54PsL5V1SlJ3jHG2PsyP9Yhr9YDcMiYL+3dOF/u+4YkleTfFp4WsA7V9G9AnjpfzndMpjPl/7D0vFgZIgmAQ8lzM33S5p2ZLtd5yRjjm8tOCVinKtM/kXJXpsvt/jPTe7g5CLjcDgAAoHEmCQAAoBFJAAAAzf7+Ea6Hrapcw7eQq666atHxN27cuOj455577mJj79ixY7Gx14IxRq3GONaX5Zx++umLjr9169ZFx9+5c+diY19++eWLjb0WrMb6Ym1Zjr2LvctS9re2OJMEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAECzYekJsLL27Nmz6PgnnnjiouOffPLJi429Y8eOxcaG1bB169alp7CoJZ//5ZdfvtjYcKDZu9i7rEXOJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAAJoNS0/gYLN58+ZFxz/ppJMWHX9pu3btWnoKcMBs2bJl0fHPPvvsRce/+OKLFx3/uuuuW3R8OFC2bdu26Pj2LvYua5EzSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANDXGWPkHrVr5B30Itm3bttjY55133mJjJ8ljH/vYRcdf2pOf/OTFxr755psXG3stGGPUaoyz9PqyadOmxcbevXv3YmOTHH/88YuNfdNNNy029lqwGuvL0muLvcuhy95lOftbW5xJAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACg2bD0BA6ECy+8cLGxt2/fvtjYSXLXXXctOv7SNm7cuPQUOMgt+Tt21FFHLTZ2kuzevXvR8c8555xFx9+0adOi43NwW3LvsrQLLrhg6Sksyt5lbXImCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACAZsPSE4CVtHnz5sXG3rVr12Jjw2o455xzlp7CorZs2bLY2Nddd91iY3Pw2759+6LjX3DBBYuOvzR7l7XJmSQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAECzYekJAMCDsXXr1kXH37lz56LjA7B6nEkCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKCpMcbKP2jVyj8oD8oHP/jBRcc/7bTTFh3/Xe9612Jjn3XWWYuNvRaMMWo1xrG+LOeyyy5begqHrDPOOGPpKSxqNdYXa8ty7F3sXZayv7XFmSQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAAJoaYyw9BwAAgDXDmSQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAmv8FXIrQ5HJAdjIAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 864x288 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot3(img2, opening(img2), closing(img2), \\\n",
    "      'source', 'opening', 'closing')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0kAAAEtCAYAAADdgyslAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAD8pJREFUeJzt3X3M3Wddx/HPFzYHMpLKk7qtrBEiT5o1EpBOklWNECaw/mFBmJjFaIKK2ag8RCBuoibDaJlE0YTENeExm2Ab2RLBiJhJ5SFswxBBQDtaBoEB9wKTh8Au/zi/yZe6lW20ve579+uVNLl7zrnP93ea3lev9zm/c1pjjAAAALByv9kHAAAAsJ6IJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAbWlV9tap+bPZxAHNV1UVVdd33eR8XVtW7jtcxsXGdMvsAAOD7McY4ffYxAPcNY4w3J3nz7ONgPq8kccJV1f1nHwMAANxdIolU1cur6jNV9ZWq+nhV/XxVnVZVV1TVzcuvK6rqtOX2/+/l7KoaVfXo5et9VfVXVXVtVd2W5Ger6oFV9WdVdVNV3VpV11XVA5fbP6Wq3ldVa1V1Y1XtPNl/BsDxUVWPq6p/Xn6eP1pVz14u31dVf11V717WmvdW1dnt+x67XPelZR16TrtuX1X9ZVVds3zv+6vqUe36o9efY932acv931pVr1+O49dPzp8OcLxU1daqekdVfaGqvlhVf3Entzm3qj64/Lx/sKrObdddVFX/tawT/11VF7bLr2u3G1X1wqr6RFV9eVlfarnu/sve5pblPl603N6ZWvcBImmTq6rHJHlRkieNMR6c5OlJDiV5ZZKnJNme5JwkT07yqntw189P8sdJHpzkuiR/muSJSc5N8pAkL0tye1WdmeSaJH+0XP6SJG+vqod/v48NOLmq6tQkf5/kXUkekeR3krx5WWeS5MIkf5jkYUluyHJKS1U9KMm7k7xl+b7nJXl9VT2h3f3zkvxBkh9K8sms1pe7cqe3raqHJfnbJL+X5KFJPp7VmgRsIMsZKu9MclOSbUnOTPK2o27zkKz2F6/L6ud9b5Jrquqhy5rzuiTPWPY+52a1Jt2VZyZ5Ulb7oedktVdKkt9I8oys9ko/lWTXcXh4rBMiiW8nOS3J46vq1DHGoTHGp7LazLx6jPH5McYXstpwvOAe3O+BMca/jjFuT/LNJL+W5OIxxmfGGN8eY7xvjPGNJL+S5NoxxrVjjNvHGO9O8qEk5x/PBwmcFE9JcnqSy8cY3xxj/FNWG5nnLddfM8b4l+Vn/5VJdlTV1qw2IIfGGFeOMb41xvhwkrcn+aV23+8YY3xgjPGtrOJq+zGO465ue36Sj44x3rFc97oknzsujxw4mZ6c5IwkLx1j3DbG+PoY4+gPbPjFJJ8YY7xxWVfemuRjSZ61XH97kp+oqgeOMT47xvjoMeZdPsZYG2N8Osl78p015TlJ/nyMcWSM8eUklx+vB8h8ImmTG2N8MsklSS5L8vmqeltVnZHV4nNTu+lNy2V31+H29cOSPCDJp+7kdmcn2b2cmrNWVWtJnprkR+/BLGB9OCPJ4eXJkTvclNWzvElbF8YYX03ypeV7zk7y00etAxcm+ZF2Pz1m/ierGLsrd3XbM446hpHkyN14XMD6sjXJTcuTHXfl6H1Mlt+fOca4Lclzk7wwyWeX03Mfe4z7ultrylFfs8GJJDLGeMsY46lZbVRGktckuXn5/R0euVyWJLcl+cE7rqiqvpH5v7ttX9+S5OtJHnUntzuc5I1jjC3t14PGGJ6NgY3n5iRbq6r/2/LIJJ9Zvt56x4VVdXpWp9jenNU68N6j1oHTxxi/eZyP77NJzmrHUP33wIZxOMkjv8d7f47exyRtPRpj/MMY4xeyelL2Y0necC+O47vWlLQ1jo1PJG1yVfWYqvq55UMZvp7ka1mdgvfWJK+qqocv5/H/fpI3Ld92Y5InVNX2qnpAVq9C3aXlWeW/SbK3qs5Y3ui4Y5n5piTPqqqnL5c/oKp2VpWNC2w878/qSZSXVdWpy4ewPCvfea/A+VX11Kr6gazem/T+McbhrE7J+/GqesHyfadW1ZOq6nHH+fiuSfKTVbVr2Vz9dr771SpgY/hAVoFyeVU9aNk7/MxRt7k2q3Xl+VV1SlU9N8njk7yzqn64qp69vDfpG0m+mtXe5566KsnFVXVmVW1J8vJ7/5BYb0QSp2V1Du0tWb2c/Igkr8jqgxQ+lOQjSf49yYeXyzLG+M8kr07yj0k+kdUHM3wvL1nu54NZnWLzmiT3WzZIFywzv5DVs0Mvjb+bsOGMMb6Z5NlZvZH5liSvT/KrY4yPLTd5S5JLs1oDnpjVKXUZY3wlydOS/HJWz/5+Lqs14rTjfHy3JNmd5E+SfDGrDdOHstokARvEGOPbWT0B8+gkn87qtNnnHnWbL2b1fsffzern/WVJnrmsA/dbLr85q/XovCS/dS8O5Q1ZfVDNR5Jcn1WYfSv3LrhYZ2p1SjYAnDhVtS/JkTHGPfmUzBNqOS3wSJILxxjvmX08wMZWVc9I8tdjjKNP82MD8mw9AJvGcmrvluV031ckqST/NvmwgA2oVv8H5PnL6XxnZvVK+d/NPi6OD5EEwGayI6tP2rwlq9N1do0xvjb3kIANqrL6L1K+nNXpdv+R1Xu4uQ9wuh0AAEDjlSQAAIBGJAEAADTH+k+47rWqcg7fJDt37pw6f//+/VPn33DDDdNmz/6zn22MUSdjjvVlnhe/+MVT5+/Zs2fq/L17906b/drXvnba7PXgZKwv1pZ5Zv/7ae+yeR1rbfFKEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAATY0xjv+dVh3/O90gtm/fPnX+9ddfP3X+rbfeOnX+2tratNnbtm2bNns9GGPUyZizmdeX3bt3T52/d+/eqfM3s61bt84+hKlOxvqymdcWexd7l83qWGuLV5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAADNKbMP4L5m165dU+ffeOONU+fv379/6vxLL7106nw4kXbv3j37EID7oM2+d9m2bdvU+WefffbU+dw5ryQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCcMvsA7muuuOKKqfMPHTo0df7sx3/gwIGp8+FE2rNnz9T5e/funTp/x44dU+cfPHhw6nw4UWb/2z3bJZdcMnW+vcv65JUkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAACaGmMc/zutOv53eg9s2bJl2uxLLrlk2uwk2bVr19T527Zt27Tz19bWps1eD8YYdTLmzF5fzjrrrGmzDx48OG32epg/2549e6bNPnLkyLTZ68HJWF9mry2bee8ye/5s9i7zHGtt8UoSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAADNKbMP4ES47LLLps2++OKLp81eD3bt2jV1/tra2tT53PcdPHhw9iFMs2PHjqnzr7766qnzjxw5MnU+J9aWLVumzrd3mefAgQNT55933nnTZs9+7OuZV5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAADNKbMP4ETYt2/ftNk7d+6cNjtJzjnnnKnz9+/fP3X+gQMHps2+8sorp81O5j72zWTr1q3TZh8+fHja7CS5+uqrp86fbYwxbXZVTZu9WaytrU2db+8yzwUXXDB1/ma2nvcuXkkCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKA5ZfYBnAg33HDDtNnbt2+fNns9zL/sssumzr/gggumzT506NC02Uly4MCBqfM58bZu3Tp1/u7du6fOv+qqq6bOn2mMMXV+VU2dvxnM3LtcdNFF02Ynyb59+6bOP+ecc6bO38x7l/XMK0kAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaEQSAABAI5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQ1xph9DAAAAOuGV5IAAAAakQQAANCIJAAAgEYkAQAANCIJAACgEUkAAACNSAIAAGhEEgAAQCOSAAAAGpEEAADQiCQAAIBGJAEAADQiCQAAoBFJAAAAjUgCAABoRBIAAEAjkgAAABqRBAAA0IgkAACARiQBAAA0IgkAAKARSQAAAI1IAgAAaP4XCsQFXlFMCTkAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 864x288 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot3(img3, opening(img3), closing(img3), 'source', 'opening', 'closing')"
   ]
  }
 ],
 "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.7.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}