ddist.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Distance Transform"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-11-18T18:54:58.907524Z",
     "start_time": "2018-11-18T18:54:58.751934Z"
    }
   },
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import ia870 as mm\n",
    "from PIL import Image"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-11-18T18:54:58.942080Z",
     "start_time": "2018-11-18T18:54:58.910395Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "a:\n",
      " [[ True  True  True  True  True  True  True  True  True]\n",
      " [ True  True  True  True  True  True  True  True  True]\n",
      " [ True  True  True  True False  True  True  True  True]\n",
      " [ True  True  True  True  True  True  True  True  True]\n",
      " [ True  True  True  True  True  True  True  True  True]]\n",
      "f4:\n",
      " [[6 5 4 3 2 3 4 5 6]\n",
      " [5 4 3 2 1 2 3 4 5]\n",
      " [4 3 2 1 0 1 2 3 4]\n",
      " [5 4 3 2 1 2 3 4 5]\n",
      " [6 5 4 3 2 3 4 5 6]]\n",
      "f8:\n",
      " [[4 3 2 2 2 2 2 3 4]\n",
      " [4 3 2 1 1 1 2 3 4]\n",
      " [4 3 2 1 0 1 2 3 4]\n",
      " [4 3 2 1 1 1 2 3 4]\n",
      " [4 3 2 2 2 2 2 3 4]]\n",
      "fe:\n",
      " [[4 4 3 2 2 2 3 4 4]\n",
      " [4 3 2 1 1 1 2 3 4]\n",
      " [4 3 2 1 0 1 2 3 4]\n",
      " [4 3 2 1 1 1 2 3 4]\n",
      " [4 4 3 2 2 2 3 4 4]]\n"
     ]
    }
   ],
   "source": [
    "a = mm.iaframe( np.ones((5,9),bool),WT=4,HT=2)\n",
    "f4 = mm.iadist(a)\n",
    "f8 = mm.iadist(a, mm.iasebox())\n",
    "fe = mm.iadist(a, mm.iasebox(),'EUCLIDEAN')\n",
    "print('a:\\n',   a)\n",
    "print('f4:\\n', f4)\n",
    "print('f8:\\n', f8)\n",
    "print('fe:\\n', fe)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-11-18T18:59:55.947384Z",
     "start_time": "2018-11-18T18:59:55.926478Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAIsAAACHCAAAAAADsPVYAAAM3mlDQ1BJQ0MgUHJvZmlsZQAAeJylVwdYU8kWnluS0BJ6lRI60gwoXUqkBpBeBFGJIZBACDEFAbEhiyu4dhHBsqKiKIsdgcWGBQtrB7sLuigo6+IqNixvEopYdt/7vnfzzb3/nXPOnDpnbgBQ5TAFAh4KAMjki4WBUfSEKQmJVNJdIAe0gTKwB8pMlkhAj4gIhSyAn8Vng2+uV+0AkT6v2UnX+pb+rxchhS1iwedxOHJTRKxMAJCJAJC6WQKhGAB5MzhvOlsskOIgiDUyYqJ8IU4CQE5pSFZ6GQWy+Wwhl0UNFDJzqYHMzEwm1dHekRohzErl8r5j9f97ZfIkI7rhUBJlRIfApz20vzCF6SfFrhDvZzH9o4fwk2xuXBjEPgCgJgLxpCiIgyGeKcmIpUNsC3FNqjAgFmIviG9yJEFSPAEATCuPExMPsSHEwfyZYeEQu0PMYYl8EyG2griSw2ZI8wRjhp3nihkxEEN92DNhVpSU3xoAfGIK289/cB5PzcgKkdpgAvFBUXa0/7DNeRzfsEFdeHs6MzgCYguIX7J5gVGD6xD0BOII6ZrwneDH54WFDvpFKGWLZP7Cd0K7mBMjzZkjAEQTsTAmatA2YkwqN4ABcQDEORxhUNSgv8SjAp6szmBMiO+FkqjYQR9JAWx+rHRNaV0sYAr9AwdjRWoCcQgTsEEWmAnvLMAHnYAKRIALsmUoDTBBJhxUaIEtHIGQiw+HEHKIQIaMQwi6RujDElIZAaRkgVTIyYNyw7NUkALlB+Wka2TBIX2Trtstm2MN6bOHw9dwO5DANw7ohXQORJNBh2wmB9qXCZ++cFYCaWkQj9YyKB8hs3XQBuqQ9T1DWrJktjBH5D7b5gvpfJAHZ0TDnuHaOA0fD4cHHop74jSZlBBy5AI72fxE2dyw1s+eS33rGdE6C9o62vvR8RqO4QkoJYbvPOghfyg+ImjNeyiTMST9lZ8r9CRWAkHpskjG9KoRK6gFwhlc1qXlff+StW+1232Rt/DRdSGrlJRv6gLqIlwlXCbcJ1wHVPj8g9BG6IboDuEu/N36blSyRmKQBkfWSAWwv8gNG3LyZFq+tfNzzgbX+WoFBBvhpMtWkVIz4eDKeEQj+ZNALIb3VJm03Ve5C/xab0t+kw6gti89fg5Qa1Qazn6Odhten3RNqSU/lb9CTyCYXpU/wBZ8pkrzwF4c9ioMFNjS9tJ6adtoNbQXtPufOWg3aH/S2mhbIOUptho7hB3BGrBGrBVQ4VsjdgJrkKEarAn+9v1Dhad9p8KlFcMaqmgpVTxUU6Nrf3Rk6aOiJeUfjnD6P9Tr6IqRZux/s2j0Ol92BPbnXUcxpThQSBRrihOFTkEoxvDnSPGByJRiQgmlaENqEMWS4kcZMxKP4VrnDWWY+8X+HrQ4AVKHK4Ev6y5MyCnlYA75+7WP1C+8lHrGHb2rEDLcVdxRPeF7vYj6xc6KhbJcMFsmL5Ltdr5MTvBF/YlkXQjOIFNlOfyObbgh7oAzYAcKB1ScjjvhPkN4sCsN9yVZpnBvSPXC/XBXaR/7oi+w/qv1o3cGm+hOtCT6Ey0/04l+xCBiAHw6SOeJ44jBELtJucTsHLH0kPfNEuQKuWkcMZUOv3LYVAafZW9LdaQ5wNNN+s00+CnwIlL2LYRotbIkwuzBOVx6IwAF+D2lAXThqWoKT2s7qNUFeMAz0x+ed+EgBuZ1OvSDA+0Wwsjmg4WgCJSAFWAtKAebwTZQDWrBfnAYNMEeewZcAJdBG7gDz5Mu8BT0gVdgAEEQEkJG1BFdxAgxR2wQR8QV8UL8kVAkCklAkpE0hI9IkHxkEVKCrELKkS1INbIPaUBOIOeQK8gtpBPpQf5G3qEYqoRqoAaoBToOdUXpaAgag05D09BZaB5aiC5Dy9BKtAatQ0+gF9A2tAN9ivZjAFPEtDBjzA5zxXyxcCwRS8WE2DysGCvFKrFa2ANasGtYB9aLvcWJuDpOxe1gFoPwWJyFz8Ln4UvxcnwnXoefwq/hnXgf/pFAJugTbAjuBAZhCiGNMJtQRCglVBEOEU7DDt1FeEUkErVgflxg3hKI6cQ5xKXEjcQ9xOPEK8SHxH4SiaRLsiF5ksJJTJKYVERaT6ohHSNdJXWR3sgpyhnJOcoFyCXK8eUK5Erldskdlbsq91huQF5F3lzeXT5cPkU+V365/Db5RvlL8l3yAwqqCpYKngoxCukKCxXKFGoVTivcVXihqKhoouimGKnIVVygWKa4V/GsYqfiWyU1JWslX6UkJYnSMqUdSseVbim9IJPJFmQfciJZTF5GriafJN8nv6GoU+wpDEoKZT6lglJHuUp5piyvbK5MV56unKdcqnxA+ZJyr4q8ioWKrwpTZZ5KhUqDyg2VflV1VQfVcNVM1aWqu1TPqXarkdQs1PzVUtQK1baqnVR7qI6pm6r7qrPUF6lvUz+t3qVB1LDUYGika5Ro/KJxUaNPU01zgmacZo5mheYRzQ4tTMtCi6HF01qutV+rXeudtoE2XZutvUS7Vvuq9mudMTo+OmydYp09Om0673Spuv66GbordQ/r3tPD9az1IvVm623SO63XO0ZjjMcY1pjiMfvH3NZH9a31o/Tn6G/Vb9XvNzA0CDQQGKw3OGnQa6hl6GOYbrjG8Khhj5G6kZcR12iN0TGjJ1RNKp3Ko5ZRT1H7jPWNg4wlxluMLxoPmFiaxJoUmOwxuWeqYOpqmmq6xrTZtM/MyGyyWb7ZbrPb5vLmruYc83XmLeavLSwt4i0WWxy26LbUsWRY5lnutrxrRbbytpplVWl1fSxxrOvYjLEbx162Rq2drDnWFdaXbFAbZxuuzUabK7YEWzdbvm2l7Q07JTu6XbbdbrtOey37UPsC+8P2z8aZjUsct3Jcy7iPNCcaD55udxzUHIIdChwaHf52tHZkOVY4Xh9PHh8wfv74+vHPJ9hMYE/YNOGmk7rTZKfFTs1OH5xdnIXOtc49LmYuyS4bXG64arhGuC51PetGcJvkNt+tye2tu7O72H2/+18edh4ZHrs8uidaTmRP3DbxoaeJJ9Nzi2eHF9Ur2etnrw5vY2+md6X3Ax9TnxSfKp/H9LH0dHoN/dkk2iThpEOTXvu6+871Pe6H+QX6Fftd9Ffzj/Uv978fYBKQFrA7oC/QKXBO4PEgQlBI0MqgGwwDBotRzegLdgmeG3wqRCkkOqQ85EGodagwtHEyOjl48urJd8PMw/hhh8NBOCN8dfi9CMuIWRG/RhIjIyIrIh9FOUTlR7VEq0fPiN4V/SpmUszymDuxVrGS2OY45bikuOq41/F+8aviO6aMmzJ3yoUEvQRuQn0iKTEusSqxf6r/1LVTu5KckoqS2qdZTsuZdm663nTe9CMzlGcwZxxIJiTHJ+9Kfs8MZ1Yy+2cyZm6Y2cfyZa1jPU3xSVmT0sP2ZK9iP071TF2V2p3mmbY6rYfjzSnl9HJ9ueXc5+lB6ZvTX2eEZ+zI+MSL5+3JlMtMzmzgq/Ez+KeyDLNysq4IbARFgo5Z7rPWzuoThgirRIhomqherAH/YLZKrCQ/SDqzvbIrst/Mjpt9IEc1h5/TmmuduyT3cV5A3vY5+BzWnOZ84/yF+Z1z6XO3zEPmzZzXPN90fuH8rgWBC3YuVFiYsfC3AlrBqoKXi+IXNRYaFC4ofPhD4A+7iyhFwqIbiz0Wb/4R/5H748Ul45esX/KxOKX4fAmtpLTk/VLW0vM/OfxU9tOnZanLLi53Xr5pBXEFf0X7Su+VO1eprspb9XD15NV1a6hrite8XDtj7bnSCaWb1ymsk6zrKAstq19vtn7F+vflnPK2ikkVezbob1iy4fXGlI1XN/lsqt1ssLlk87ufuT/f3BK4pa7SorJ0K3Fr9tZH2+K2tWx33V5dpVdVUvVhB39Hx86onaeqXaqrd+nvWr4b3S3Z3VOTVHP5F79f6mvtarfs0dpTshfslex9si95X/v+kP3NB1wP1B40P7jhkPqh4jqkLreu7zDncEd9Qv2VhuCG5kaPxkO/2v+6o8m4qeKI5pHlRxWOFh79dCzvWP9xwfHeE2knHjbPaL5zcsrJ66ciT108HXL67JmAMydb6C3HznqebTrnfq7hvOv5wxecL9S1OrUe+s3pt0MXnS/WXXK5VH/Z7XLjlYlXjl71vnrimt+1M9cZ1y+0hbVdaY9tv3kj6UbHzZSb3bd4t57fzr49cGcB/Igvvqdyr/S+/v3K38f+vqfDueNIp19n64PoB3cesh4+/UP0x/uuwkfkR6WPjR5Xdzt2N/UE9Fx+MvVJ11PB04Heoj9V/9zwzOrZwb98/mrtm9LX9Vz4/NPfS1/ovtjxcsLL5v6I/vuvMl8NvC5+o/tm51vXty3v4t89Hpj9nvS+7MPYD40fQz7e/ZT56dN/AC1d8ByJG09FAAACiElEQVR4nO2a2XKDMAxF5U7//5fdh9LQgHZZcjLRfUhhCtaxFtsYAFqtVqvVarU+T3NRO18rGpnEsVELWCZ5YpSTZRKnM0Lj9ctEjqJyx2hiZ5S7Mlkmbv7CZ/VYIHdvOftMZ49dpI7+2b//z5FGBpZ768/xCY8y34ZrJwwWJkZijZHO894qN7AMgIPGbkzVCUuMjmZvkVLcopKdxeqW36s1/B4WiyzcuSw2Dy5ZvyySzHKWgDlnT6lu1cRoxjjUUtg4/eJfqizzy5+iqyapujXA61ZuvMH3qqO14uKgYamoIYByv7DdUrEscwzf0CvlLkO6tJQlYwCsX6pSVmew2DPFLKxRPhBZMLhdvo7SUgZtWKjpJBi82S3jC9FDiSXDMVSbG/xCdq++ph3ruuKxBYBc7+YNLIwozEQYe75kesY+B5TTcHWUmjSI4S3zNGFaGFdLU1ga4yuHmex9KVy4B3aspSib9XXk2WdIyhSm8/VzgNUvxQPLIXwtlfuYRvV0y9xor+lymm3zkXVu/OD5aIdfPHNA9aYH+6yW4BWuSY6l+sGEYclBYVqlWbK8QrdLsmx4ctzxTOKraYAUIudeEMAo2+fdskdGuNryvjEo0dQrvZtQvJ9eZUpsSGQpHGcqYyR1S2KpHH61z/ajAEte143zMNma4JfSd2t6YwG/jMftrLmyvaCTh71EpUi6KI281xxwUeISwsjiI1n8rZ/nm6lhTDJ9HXk8YqPRxsgRnPH4WcxCGhMv0tP4WZ6MDPTQSONmGcTJcGzyxFjGFYVwkkk+loH9OY+dMB4WNAFuMbMDOViuRrDwuDwTnxtx8x6YGMvg6qmYJWo9T1s+PGi1Wq1Wq9VqvZx+AN55dcQIaGKBAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<PIL.TiffImagePlugin.TiffImageFile image mode=L size=139x135 at 0x11158CEF0>"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "f_pil = Image.open('data/gear.tif')\n",
    "f_pil"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-11-18T19:01:53.450151Z",
     "start_time": "2018-11-18T19:01:53.440651Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAIsAAACHCAAAAAADsPVYAAADmUlEQVR4nO2b23LrIAxFpc75/1/WeXB8AesOIp0OeqmTxnixBQKBQYJfYz/fBnjYZuFts/C2WXjbLLxtFt42C29/jQUnlAEwhQUbmAGwOT7Cx1UeJsmCz0tiv1/F0tWezo8IRHmetI/w+ivPmINUWZa2+p9PLRdG206KBYG4x3co4WLz/ahRhgBPlON7RKBouhNg6Sv6eD7323ji5WdBwI7mqQw9sTKihFg+D2lh7udTI0QqG/0X+O3bKew14ZtE6/n3j/w1QKCz0ThvOgnQd09El6PEaF89QDx3RVkAKKLlKYnLR3GWUMMM9e3SeR1BqG/burjklWECZvoe71LzVL47He2FEAYnSU4zWQiRJrFYscnbjwhCcVHgALUUR/kjjbdBsdpcJr4kzFWbxXmjqrGDZWBmH7O1uuhNz8MyTRijFyhtd5FnbpN1afOfGWYFByW+zIgrXYF6iVqsq4ABAJFoUaw7GQ4aQSB1DJguzFUuC6P26aooR7wuenwpghH0/sY6puR6g6VCGLEVqv1opYNAYTFiQQGKxBJKmyehCCyppRwvjVgy33brNvPptYhzm0RZqIwQdRXFSt3E06ycM1g0+pyhlqYvXIu7VBTtBBQrbyzyE99g9LlUjSxSR9JYxM5XQqKx1JBoHUJiKSLRt5v4CFhEAqoubJ/G2PpjyJRAwbJUBhZlR1KIL4VOgnh+VErDw6hjYynNu2hlPKLSZsMU/Y3xSJpMW+PR/DQABBIzPypBCcW6MhR9F0dhqZxjsqat11WgaAm6yFKligIjseAX1lTVMQAACkKvqLiRk5TQSDBafDk33Uf2sCK2aC/rKOe25HrdMpSla81WnUxdpsUZexnSYlk5ECxc3zWFMdrLJUvp+oePBR7tvzzM6D66RyVhOyFklpN0XWiOa5wZhVv4gQ51D7N6IWti3Se7GetHlw3HGcfLPGFdCl9ZCsY6rJzvRVgQY90Kr3jgw3f6CD9vDkb8g1GHunVBAIouVgUXB5wsbL6nPweBgEjeoUmzcCTPV8ub8wLNfcp+UZblZdh+QObygPH7NXtuonmLujtF0b+KTc4WnD038Sy+Oy3x7sJOaTIspyjSEYHsBntOl3dFpdMCtSzUNUfu5EROmYQubX35x6fGz6E8gAC700d0/yuOM5yTPJ6O7PerWPraj823BnWh5npw6jczbxzNGP7aWctZtll42yy8bRbeNgtvm4W3zcLbZuHtP4R8EC4cGxaVAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<PIL.Image.Image image mode=L size=139x135 at 0x11155CE48>"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "f = np.array(f_pil)\n",
    "fn = mm.ianeg( mm.iagradm(f>0))\n",
    "Image.fromarray(fn.astype(np.uint8)*255)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## City-block distance"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-11-18T19:13:28.017188Z",
     "start_time": "2018-11-18T19:13:27.887107Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAIsAAACHCAAAAAADsPVYAAAGrUlEQVR4nO2b0Y7kKAxF7wr//+/uC5bmAWxMgo1Jd9VIu0GqUXoqgRvn2AaH+uffyrXmPgUE+YuoEMl3VIioUD+LiAoRaoUcEypIjoNPXkoFAPkLAFQYWhtntS+hx+YoaJSXUgFA/jKDq5ZJGLgea8lL6VrM3U/C5Dsd3mphJBrlpTBweSyTMPnOaBkKknZ5eVnb5eVlreXlZW2Xl5e1lpeXtV3+W7zISH+Vl8ttT8K+zUsfvVgpOvjXeGl3X3onhbl/hwKWJ/ZxXubem1locdf1C7zQbIreSTOG0VUrKhtzfIQXQcQOL+MzikidjIL6GV5mW5AcgxmMQgSAGagMlOFSH+bFjH/npTIMRpN3h1ryUlhEGDH6X2xkMVD0PHD6GeWQ1W4nw3SzNNfmQqTjF4KYpf4mL3puAWPyHgBEJMAaVRiqmpbf4kU6Zi6iZlijrZ31UFQSiSlqbeZMqEnwYs+2nTa7jCdjVIGKRN8eoTOR5oyXNr6WEiYLwajS7CC5ooNFITgnvAxjkJpjxQuKgstAAYZ3B0pyvKiGoUXN4fEiUoxvb5k54wUtN/faDeDzov2D1Bw7x97zorfT02HpkOR4oUYBCsWstHO3rVxckjK83PvfBzwqmxPAXHCJuGleTIKs+xCTsIuepL58wstwKFyPbsNEzizxViKvPJsjXoaAZt+ojpmbr0y+/IgXGI9y7JJ05vFACuXjy5iey1iRWaL4IiKY0SB/wktkiKtdssHf+vI5L4AuGvy0lOFF5vz0E15YpbielOVluuYBL5klZMCLOWsNG7DnZe4xnjSQDxebp77kXlUOwaLlbo7eZdQCLWAuxvz3j6ocgrsqx304FrPzuR5wlw9JVQ7BcpHXLZt/72NFYpj7Cc958RSt5HgXaSoChhOf8+KKWd2dF1/AshoEvJ4Jh7y05mcCz6ETK72m5YiXKfpm3zeiAqhzmWUl5oyXGkrx8hEqEms98ngJJ9CelICXMYL7cXhJSFnlApcXoK8/Ail046VmpHg5yeGl64AXctWrZ152zxQI57tetBYCIyl0jS8pKUH+dy7h5iZnvCS0BGrc6/sVcfd04SWlxVUTXB+GFhVj/s3x4quJ5wzZriU1pbUs1cR5Oldng6amTImworZqVTYfAVRHbsyJCdd9c0jXml8mHw312UbpEnpTsgqibmp01XvunZTSlSyCqHMzCNSvPvVESgGBVkE0tmOcAoaUOt6K7qR0JQe8BOpXUmrskKb1gL58+LEdU7y0K9Ji/DXsMu4yuFSqza93rcGeF+MGCjdPH7hz5VYH+KmYXW5MdV2lJlHHodfI7zxcTyfV6DywFzZCV+p6VqYJ83xSjQwh6MQZQGaY9643zziVj4SXHmMQSRnLi/U80R0jroDKlfXKyy4DAE5Ef7aqKjTq+2aIPS9NB07ykcTrMAVNkIzDzZSoLbYP8tGQ4hbYjJYLL5swQ01LOh+JlGCJr7ObGy+hmF7TOZi/1G2hD+zy4uXsVjvhvvY64GW7snd5Eft4LuS2IAdsxPi8wPFqwA25TUsQd8NCX8SL49V9xNM8PSlaX/qAl01L5PmLlPGu+Tkvh3MGGVGLN9fbOubFXOvEz3g9PnY0Lr7a8nK5YgzlhPIwLo3oe3/YCV7mYuwo6TurC1/LHH3vWhK83A1kSp6LT2CWMeoqKSV4GTekhWtTLV/ko3Dmw+jQPuTFRDlpQUFyaxW+Hcife16GgUy+JJeXjVnGxopuWL3HY1404PqTop1RZqe3O1FPeMm0/UOiPnLXZuyS5qXf1ebtQlbLLfAneemtpPaS7eugNN6zTHtyU7xMN7RZVux4wbzF05bNUrzM+6w2bxcSz8hYwlr6MB/1GnAoJV+T7YsbVZDjhcgG3FjKjheNv5c3j7zn5QJrwbYUuZGiY7dde8YYAS8yOoNNxKUf8yJSCghmI2yNeBHPu0y0d6XINC/TC027tXzBi/EeO13elCLT71nmHYTGt9v4DDUH22swb837YXwZ/ZpIwyPumADbIedrWOGNhENe7ObTaWt5n+QYLi4/HQC2yD7gRVor/wqY0y8UFFg52+x+3ks54GV4a7m+mbeLwSZGt2kBe07OeRGzFJDZ5n795YT+XsFc9vu8jO7HVd7ws5bP8UI602C22xig+wNgfKdp+W1expEpE5nReTpTBgE+yAvMTxH63U/CdFPrgPc7vNj7mR4L7L1+gxe5+0lY//KLvMhxMX3Yx/J9XmxAmx/LX+ZleiwvLy8vqurlZa3l5WW0l5e9lpeX0f5fvPwBLzPCQTK6J6AAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<PIL.Image.Image image mode=L size=139x135 at 0x1116000B8>"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "d4 = mm.iadist(fn)\n",
    "Image.fromarray((d4 * 255 / d4.max()).astype(np.uint8))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-11-18T19:14:27.758018Z",
     "start_time": "2018-11-18T19:14:27.752527Z"
    }
   },
   "source": [
    "## Chessboard distance"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-11-18T19:15:04.812704Z",
     "start_time": "2018-11-18T19:15:04.704288Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAIsAAACHCAAAAAADsPVYAAAFv0lEQVR4nO2buZIcMQhAscsB/P/HQuhAB0iNENLMlgM3wZb62OHt09nqWYA33njjjTfeeOON/y/wS5/z+xsfgovyYXyBBZcHh3HJMqQkPcRPaG69oC2Rnidybv5hFgNDeoSre36QBbVesBz10+am09q69EIKQ72VVMIOeVpd9+2Fan5AaC2GtPHgOcoJC05H1CiqDTTdCG8a8YkXnEYSAsTWj6jUSSHAdu0s/hzdjSCtSIzSagaglG1LuRhmDlgECRhBah7ilg6llamzARD7f8RXWEoeRhDN2WF6CQyVkqTiuI4I2DbSRzJDpdeIE1pOvRQaqM31GTwrQP++L7AI9lbizoKieVlJMCPllIWtDceMsSL1lkLMkIk9i3YBqQNtvzCZcawAlA6VMpPxgmA/jGzJmnGtnMSeRZBqHlG0XjL5fStJJzkWACYgMH+3a2ZpRdsKbtASLIKllszYZks1Q2ylD4DRQJPrR6Yvtw/V0uDMt+JIumVhAh5b4mQmZQVA4saTYRFs9aQxmIGUFZbNOJOrIyaY+6g1k7EiDMBxn0qxSFus2bBmTCUsrPC2hg7ngJ2ZwIqyXbE4q47YTGhF6W5YxFnKR2Y2VvYR1ZE464KoN0FsZbvSjNsL00PO0ozAzgq1IXFBFLIIMj3krMxAwko9gz5O7EWwDLl+xhkta4UW48ymT5chN2OGc1YYZDnkbccX1geyOb1Fk4wVCVG2LNIfyB7pLVrKCscoB+NuZCZpBUKUmAXXR9O4k7JSzl2s63BIAY8j/yE1tjKcTbOgebhxUs4HfdFrjjwrcSxYSpPFtBlrIrISzgOrOhJkoE405XQOytihiCsrE1mKBQSQqRLZ/Vt7z2gmaWVNE64ZmASQx02mqDdBzsqKJl4z9HrKmFG4rRWfJp6nixVNMhe9uYntRRl2ykDawtd93N/PjSbNs/jsTTJaEZtTIpKTOWDVnbTcrHC7UKwYMxyQ7FiGGtBxngS0bwkAAJLdiqiLcEZz2G4ONmMilnmoq2SCbDZwuWHatsJjW5nNHLMgAHlmmMRMVlJvncZ4WR6ez0dKMpupViyCQmnw6jDYSPRZLMnYa7sVTW8rS2NhhtZTwIJFEJj0p/lEtdKSj43YUHuHKNqkkyw1af3ZYcqpMbUh3ZshRhBawYRzI5AZeUMrj87rsHGzsqqnaG7MW3E6r8PWfmdVT9n5KGHFye61mU50xmLno4QVJ7tnhp0LGRazqN5amdYLLtuSY8/SNxkSViZQGG9eHOZZ+ivxlJXFU8fIFr2/iefG+nY1ZUVzBmZoupZl6ShHVubsAxvK+XykKMAAZ1bG7AMbMSwH3YgFyzwIAKkeFDfRNhkBSlBLwRwAWB++Uj1ouUiwv0gc1VL4TFISZ8YVyZnpPMcs7UMz4wonzUDQlw724IPR1t3tX48zp2sG85vDR3tW2jSTNbOIjBfbNhdW6rmDAf+KBW2fXVqpZyMz29f3Wxa0Dz6BlXp+bSYc/nMsQDrAh1bqIL0yEw//KZY6eOhz0NIKCQRmiPdi9l64feTGStszXZnZf9tjv64bl22hlXp92WY2OLv1rikP045vpd7hmqFtJZ28+5wlOVYqoWMGGXZizt7D7q1UxqcBAW+r8Y4FQTdPIis1ryvgwz6tQb1zx1YABjOPV7nfYOkbXDsr5VZd/fwEC4Ox0QuLZGxXyskvkmVZcNxpEzNDKda0F9a2KvYz0RnL1PDabiQOWJawakTON5gsy2i57dNSQalYWmRqv2K2977G8kDjtglR3mdC+epQKYppRJJvvlffsRYBBikry/5qVWqDfbzelOzy7u773tzSSa+kygHyhMmuNa9YmpTSVJsVaSfCF79fZumZhtWTeS4RYGdG+hEW2wB6NxZTvjRzweJ071rSazdV9On//DDb9OMjwXH8+oxlfJt4p6PHp/8LJW7x38dm2fbGG2+88cYbb7zxxn8SfwHrMNnXu/PwOgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<PIL.Image.Image image mode=L size=139x135 at 0x111600DA0>"
      ]
     },
     "execution_count": 27,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "d8 = mm.iadist(f, mm.iasebox())\n",
    "Image.fromarray((d8 * 255 / d8.max()).astype(np.uint8))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Euclidean distance"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-11-18T19:15:35.123203Z",
     "start_time": "2018-11-18T19:15:34.997624Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAIsAAACHCAAAAAADsPVYAAAHRUlEQVR4nO2by5LjNgxFb1Ip8P8/l1jNAm+KlCjZrlQlwqJHbcviER4XoDwNvPbaa6+99tprr/3/jL50nb+/cRFaHN+0L7DQ8peb9pBlWJ/S689pnvqF8lFanv6NGCWYoKHVOT9koYgLIZxBxS23XfQ0XyhgiAwJfoiSRZv2zzMWgMDEujABHFisL93OnBt+qdc2z5CGorqB/PUbdscvJLccaCzLtS7HFHnzqJ7uxajQqD9IYQAitA7N7Ar+bRYmAhMsSQwAjYAuSds8gQ7VvUF21y9CQ8kbjQASmEQlHoqP7djNOiJYljTyl8gSJ1HVjxDvBOwBC7DSkUhnj8md2r6tL54kNGWxIw4SU5/vsjBAzTMDHph5O3Ld4a0I7bBECWii6AquIA2asZ7OYlZnehffYNGELb+J51t6wQ+FgD2LS0F9ysJkKRnpCFJdObIAECD+SQ9gCYdfXGqjUUqYgQU9DgePnoRrg4XTrOSrtxhWkpbYcgmP0+sEnKXxZh15OBoRkZCI5k/1gy1ZwPmEC8nbYhlqWcIjRDOWLkLcc/akS33EwgSgl1omCZOLbOXpKW+ioHClM3sxYpNXqeUcJlTHVCpObaEz87nO7MUotT6tZQ1TqF85m130rKC46z19ymJGVsvUdCsyTV73DSOOLiN0yjKpD/JalnG35W6k1oeQ9eiNFzAnLJw7sB5ZBalTaMICdwxrVu82gbMYJRgTFhlqQRIn3QqM+dIkZ6yss/B+kLucHYKoIItTGnHTDZQjXz1GwhXRKYvuxhrykNB0zG9EEqbMwu6jboO6YqQ84sV0deEXMmUJYTF/EIHaRGE6sa1tPiqKB2ChMxc1LZJrVSPCov4gWGUfPmBQvqQqnqTPUvIu9SX2zCYsCJS0xS8k4ptuHDFC9BOUS5a4zxbCAkcxuZveg3EYCJ+j3NtPu7BEmKgdkkVS04tIKryDmc9RzlkOItbI4kRaTi1njDmxBIjz9H7aHpcsXh9Ddk7idGwWHiDGqLv35zqKfU9IrlojOpRTse5bh7a7TTtjYR9WsuQ23yvKXk1RvBk30xYbLScgJ31gFSN/ANbKvgxEjXRk0C7QhiDlAjqCINR4mwUse2KS2TLtnxXJ57v68FKl5Lha6gMrmtOZQWD8Jock9vnOHTN0wwQxTr5TmvOZQWEOHLaMdgHJ3l4SJDUffyHtDmY0530a/th0ZVZOdd21demMPDt3pzdemGdvuTpzZ+YeQy4zg/ua5INnzRVnMZEkKADgDl6SXLFseEXt8skGez961BtlhNrj2ZuvH/YjG+YOLWC2xBbI5bkrFieZtYDx6g0M2xpeGd31SyKZtoA8J6Sr89XekM5awIKFCSxTUV+MbQYCChhmAH2BEz2AV65Z9Wl/QtuL3tXdjryVPCO1y33SBlIPWMFc9cbcj1CKnNEcJn+OuY8bedMXe+i6iNNZb8SyH3UA5DDpI+IRZpbnLUfZ1WtNXfOgH3GMSY1B6M2v20VmkeI0yC4szhOYu/0oDY7EQOPxXVlf4pSK3GTXTnrQjw6hkUhHfA4wFhdR++FdXnJcsxRt0c0FaewaB1hi6fA4debqnBOOSxZBkfFSarshCqgxQExlre77duHxwPglH/SAhOI3Q15A8rMd71ImAvOFxYlDER76RQPk1+lxNepxh0NVG44WUrcJqt7JTRb/XtcWFhWXyznMeGmWHbwchwyzb1WWHeB8DxuCn/YCXZKFwIdkISUBlzB16DSnX6ov7/6kB6R7du1NcdFqLkWtPmGUlE2Sax++x4IikFLY3cEovxnmU74MtgfJBUC3+3Qh8iatWULmtDZelOFhYndMktwZ/x2WvJCNEK2GsLJYmMIXGBr33ZkhfTJQ7HcLfs1dOafLpM++AclvntrW97Bp+6MyI7l8uHpH3hbWlL22ne+nbRWXGY7arsbB0hH7RPiHP2MR9Z2MeLMFtENLPR1U7mqpDb+Q35uGxrqC/ZtERxyj3igqd+C+z0J2ffsN8A5gMEN++v79MFherLWjL/oz15PWtlgKmgibaoyc3q+Ds8eS7ubwiMV3LWPdii8YQFY5ugzS1bybuaiQ9dlJ6sRcPqFyn+ZLtZbSNXTuwBLlE69OJ4xPWLo3gsibY74Ag+I6z/nVt1m0hILAXu8Zy6coK5+yHf9ajIjibodHuDEOhLC5C/ftXp/mGCFyRsbUFJuAVD6/YcmreoeasKAOCZed6B4L1crQjiOSx9E4rS9r7hDjxv/M3PdLuaAJCHWQbjk6yHYfvalbJq38CyyH6cDT1FKDug4KpaSYeDtlnj1rlmFWv4e21KB8mE7dds1jFvOUfp9J6ZBr9Ww75hFLeXjN/ljd5HA67v2Mxa4uvRgWL3u+ciNHPmUJlPhysfxkOrSjHXvy9wHDM0oDTG8Mffp3LHXazoNKvPfAKx//zU9GmbXLW/bXZyxZ33e2Haf26d9CfeiKn9mTOn7ttddee+2111577b9nfwAWeRidwOGzywAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<PIL.Image.Image image mode=L size=139x135 at 0x111600A20>"
      ]
     },
     "execution_count": 28,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "de = mm.iadist(f, mm.iasebox(),'EUCLIDEAN')\n",
    "Image.fromarray((de * 255 / de.max()).astype(np.uint8))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.5"
  },
  "toc": {
   "base_numbering": 1,
   "nav_menu": {},
   "number_sections": true,
   "sideBar": true,
   "skip_h1_title": false,
   "title_cell": "Table of Contents",
   "title_sidebar": "Contents",
   "toc_cell": false,
   "toc_position": {},
   "toc_section_display": true,
   "toc_window_display": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}