train_shapes.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Mask R-CNN - Train on Shapes Dataset\n",
    "\n",
    "\n",
    "This notebook shows how to train Mask R-CNN on your own dataset. To keep things simple we use a synthetic dataset of shapes (squares, triangles, and circles) which enables fast training. You'd still need a GPU, though, because the network backbone is a Resnet101, which would be too slow to train on a CPU. On a GPU, you can start to get okay-ish results in a few minutes, and good results in less than an hour.\n",
    "\n",
    "The code of the *Shapes* dataset is included below. It generates images on the fly, so it doesn't require downloading any data. And it can generate images of any size, so we pick a small image size to train faster. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/usr/local/lib/python3.5/dist-packages/h5py/__init__.py:36: FutureWarning: Conversion of the second argument of issubdtype from `float` to `np.floating` is deprecated. In future, it will be treated as `np.float64 == np.dtype(float).type`.\n",
      "  from ._conv import register_converters as _register_converters\n",
      "Using TensorFlow backend.\n"
     ]
    }
   ],
   "source": [
    "import os\n",
    "import sys\n",
    "import random\n",
    "import math\n",
    "import re\n",
    "import time\n",
    "import numpy as np\n",
    "import cv2\n",
    "import matplotlib\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "from config import Config\n",
    "import utils\n",
    "import model as modellib\n",
    "import visualize\n",
    "from model import log\n",
    "\n",
    "%matplotlib inline \n",
    "\n",
    "# Root directory of the project\n",
    "ROOT_DIR = os.getcwd()\n",
    "\n",
    "# Directory to save logs and trained model\n",
    "MODEL_DIR = os.path.join(ROOT_DIR, \"logs\")\n",
    "\n",
    "# Local path to trained weights file\n",
    "COCO_MODEL_PATH = os.path.join(ROOT_DIR, \"mask_rcnn_coco.h5\")\n",
    "# Download COCO trained weights from Releases if needed\n",
    "if not os.path.exists(COCO_MODEL_PATH):\n",
    "    utils.download_trained_weights(COCO_MODEL_PATH)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Configurations"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "Configurations:\n",
      "BACKBONE                       resnet101\n",
      "BACKBONE_SHAPES                [[32 32]\n",
      " [16 16]\n",
      " [ 8  8]\n",
      " [ 4  4]\n",
      " [ 2  2]]\n",
      "BACKBONE_STRIDES               [4, 8, 16, 32, 64]\n",
      "BATCH_SIZE                     8\n",
      "BBOX_STD_DEV                   [0.1 0.1 0.2 0.2]\n",
      "DETECTION_MAX_INSTANCES        100\n",
      "DETECTION_MIN_CONFIDENCE       0.7\n",
      "DETECTION_NMS_THRESHOLD        0.3\n",
      "GPU_COUNT                      1\n",
      "IMAGES_PER_GPU                 8\n",
      "IMAGE_MAX_DIM                  128\n",
      "IMAGE_MIN_DIM                  128\n",
      "IMAGE_PADDING                  True\n",
      "IMAGE_SHAPE                    [128 128   3]\n",
      "LEARNING_MOMENTUM              0.9\n",
      "LEARNING_RATE                  0.001\n",
      "MASK_POOL_SIZE                 14\n",
      "MASK_SHAPE                     [28, 28]\n",
      "MAX_GT_INSTANCES               100\n",
      "MEAN_PIXEL                     [123.7 116.8 103.9]\n",
      "MINI_MASK_SHAPE                (56, 56)\n",
      "NAME                           shapes\n",
      "NUM_CLASSES                    4\n",
      "POOL_SIZE                      7\n",
      "POST_NMS_ROIS_INFERENCE        1000\n",
      "POST_NMS_ROIS_TRAINING         2000\n",
      "ROI_POSITIVE_RATIO             0.33\n",
      "RPN_ANCHOR_RATIOS              [0.5, 1, 2]\n",
      "RPN_ANCHOR_SCALES              (8, 16, 32, 64, 128)\n",
      "RPN_ANCHOR_STRIDE              1\n",
      "RPN_BBOX_STD_DEV               [0.1 0.1 0.2 0.2]\n",
      "RPN_NMS_THRESHOLD              0.7\n",
      "RPN_TRAIN_ANCHORS_PER_IMAGE    256\n",
      "STEPS_PER_EPOCH                100\n",
      "TRAIN_ROIS_PER_IMAGE           32\n",
      "USE_MINI_MASK                  True\n",
      "USE_RPN_ROIS                   True\n",
      "VALIDATION_STEPS               5\n",
      "WEIGHT_DECAY                   0.0001\n",
      "\n",
      "\n"
     ]
    }
   ],
   "source": [
    "class ShapesConfig(Config):\n",
    "    \"\"\"Configuration for training on the toy shapes dataset.\n",
    "    Derives from the base Config class and overrides values specific\n",
    "    to the toy shapes dataset.\n",
    "    \"\"\"\n",
    "    # Give the configuration a recognizable name\n",
    "    NAME = \"shapes\"\n",
    "\n",
    "    # Train on 1 GPU and 8 images per GPU. We can put multiple images on each\n",
    "    # GPU because the images are small. Batch size is 8 (GPUs * images/GPU).\n",
    "    GPU_COUNT = 1\n",
    "    IMAGES_PER_GPU = 8\n",
    "\n",
    "    # Number of classes (including background)\n",
    "    NUM_CLASSES = 1 + 3  # background + 3 shapes\n",
    "\n",
    "    # Use small images for faster training. Set the limits of the small side\n",
    "    # the large side, and that determines the image shape.\n",
    "    IMAGE_MIN_DIM = 128\n",
    "    IMAGE_MAX_DIM = 128\n",
    "\n",
    "    # Use smaller anchors because our image and objects are small\n",
    "    RPN_ANCHOR_SCALES = (8, 16, 32, 64, 128)  # anchor side in pixels\n",
    "\n",
    "    # Reduce training ROIs per image because the images are small and have\n",
    "    # few objects. Aim to allow ROI sampling to pick 33% positive ROIs.\n",
    "    TRAIN_ROIS_PER_IMAGE = 32\n",
    "\n",
    "    # Use a small epoch since the data is simple\n",
    "    STEPS_PER_EPOCH = 100\n",
    "\n",
    "    # use small validation steps since the epoch is small\n",
    "    VALIDATION_STEPS = 5\n",
    "    \n",
    "config = ShapesConfig()\n",
    "config.display()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Notebook Preferences"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "def get_ax(rows=1, cols=1, size=8):\n",
    "    \"\"\"Return a Matplotlib Axes array to be used in\n",
    "    all visualizations in the notebook. Provide a\n",
    "    central point to control graph sizes.\n",
    "    \n",
    "    Change the default size attribute to control the size\n",
    "    of rendered images\n",
    "    \"\"\"\n",
    "    _, ax = plt.subplots(rows, cols, figsize=(size*cols, size*rows))\n",
    "    return ax"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Dataset\n",
    "\n",
    "Create a synthetic dataset\n",
    "\n",
    "Extend the Dataset class and add a method to load the shapes dataset, `load_shapes()`, and override the following methods:\n",
    "\n",
    "* load_image()\n",
    "* load_mask()\n",
    "* image_reference()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "class ShapesDataset(utils.Dataset):\n",
    "    \"\"\"Generates the shapes synthetic dataset. The dataset consists of simple\n",
    "    shapes (triangles, squares, circles) placed randomly on a blank surface.\n",
    "    The images are generated on the fly. No file access required.\n",
    "    \"\"\"\n",
    "\n",
    "    def load_shapes(self, count, height, width):\n",
    "        \"\"\"Generate the requested number of synthetic images.\n",
    "        count: number of images to generate.\n",
    "        height, width: the size of the generated images.\n",
    "        \"\"\"\n",
    "        # Add classes\n",
    "        self.add_class(\"shapes\", 1, \"square\")\n",
    "        self.add_class(\"shapes\", 2, \"circle\")\n",
    "        self.add_class(\"shapes\", 3, \"triangle\")\n",
    "\n",
    "        # Add images\n",
    "        # Generate random specifications of images (i.e. color and\n",
    "        # list of shapes sizes and locations). This is more compact than\n",
    "        # actual images. Images are generated on the fly in load_image().\n",
    "        for i in range(count):\n",
    "            bg_color, shapes = self.random_image(height, width)\n",
    "            self.add_image(\"shapes\", image_id=i, path=None,\n",
    "                           width=width, height=height,\n",
    "                           bg_color=bg_color, shapes=shapes)\n",
    "\n",
    "    def load_image(self, image_id):\n",
    "        \"\"\"Generate an image from the specs of the given image ID.\n",
    "        Typically this function loads the image from a file, but\n",
    "        in this case it generates the image on the fly from the\n",
    "        specs in image_info.\n",
    "        \"\"\"\n",
    "        info = self.image_info[image_id]\n",
    "        bg_color = np.array(info['bg_color']).reshape([1, 1, 3])\n",
    "        image = np.ones([info['height'], info['width'], 3], dtype=np.uint8)\n",
    "        image = image * bg_color.astype(np.uint8)\n",
    "        for shape, color, dims in info['shapes']:\n",
    "            image = self.draw_shape(image, shape, dims, color)\n",
    "        return image\n",
    "\n",
    "    def image_reference(self, image_id):\n",
    "        \"\"\"Return the shapes data of the image.\"\"\"\n",
    "        info = self.image_info[image_id]\n",
    "        if info[\"source\"] == \"shapes\":\n",
    "            return info[\"shapes\"]\n",
    "        else:\n",
    "            super(self.__class__).image_reference(self, image_id)\n",
    "\n",
    "    def load_mask(self, image_id):\n",
    "        \"\"\"Generate instance masks for shapes of the given image ID.\n",
    "        \"\"\"\n",
    "        info = self.image_info[image_id]\n",
    "        shapes = info['shapes']\n",
    "        count = len(shapes)\n",
    "        mask = np.zeros([info['height'], info['width'], count], dtype=np.uint8)\n",
    "        for i, (shape, _, dims) in enumerate(info['shapes']):\n",
    "            mask[:, :, i:i+1] = self.draw_shape(mask[:, :, i:i+1].copy(),\n",
    "                                                shape, dims, 1)\n",
    "        # Handle occlusions\n",
    "        occlusion = np.logical_not(mask[:, :, -1]).astype(np.uint8)\n",
    "        for i in range(count-2, -1, -1):\n",
    "            mask[:, :, i] = mask[:, :, i] * occlusion\n",
    "            occlusion = np.logical_and(occlusion, np.logical_not(mask[:, :, i]))\n",
    "        # Map class names to class IDs.\n",
    "        class_ids = np.array([self.class_names.index(s[0]) for s in shapes])\n",
    "        return mask, class_ids.astype(np.int32)\n",
    "\n",
    "    def draw_shape(self, image, shape, dims, color):\n",
    "        \"\"\"Draws a shape from the given specs.\"\"\"\n",
    "        # Get the center x, y and the size s\n",
    "        x, y, s = dims\n",
    "        if shape == 'square':\n",
    "            cv2.rectangle(image, (x-s, y-s), (x+s, y+s), color, -1)\n",
    "        elif shape == \"circle\":\n",
    "            cv2.circle(image, (x, y), s, color, -1)\n",
    "        elif shape == \"triangle\":\n",
    "            points = np.array([[(x, y-s),\n",
    "                                (x-s/math.sin(math.radians(60)), y+s),\n",
    "                                (x+s/math.sin(math.radians(60)), y+s),\n",
    "                                ]], dtype=np.int32)\n",
    "            cv2.fillPoly(image, points, color)\n",
    "        return image\n",
    "\n",
    "    def random_shape(self, height, width):\n",
    "        \"\"\"Generates specifications of a random shape that lies within\n",
    "        the given height and width boundaries.\n",
    "        Returns a tuple of three valus:\n",
    "        * The shape name (square, circle, ...)\n",
    "        * Shape color: a tuple of 3 values, RGB.\n",
    "        * Shape dimensions: A tuple of values that define the shape size\n",
    "                            and location. Differs per shape type.\n",
    "        \"\"\"\n",
    "        # Shape\n",
    "        shape = random.choice([\"square\", \"circle\", \"triangle\"])\n",
    "        # Color\n",
    "        color = tuple([random.randint(0, 255) for _ in range(3)])\n",
    "        # Center x, y\n",
    "        buffer = 20\n",
    "        y = random.randint(buffer, height - buffer - 1)\n",
    "        x = random.randint(buffer, width - buffer - 1)\n",
    "        # Size\n",
    "        s = random.randint(buffer, height//4)\n",
    "        return shape, color, (x, y, s)\n",
    "\n",
    "    def random_image(self, height, width):\n",
    "        \"\"\"Creates random specifications of an image with multiple shapes.\n",
    "        Returns the background color of the image and a list of shape\n",
    "        specifications that can be used to draw the image.\n",
    "        \"\"\"\n",
    "        # Pick random background color\n",
    "        bg_color = np.array([random.randint(0, 255) for _ in range(3)])\n",
    "        # Generate a few random shapes and record their\n",
    "        # bounding boxes\n",
    "        shapes = []\n",
    "        boxes = []\n",
    "        N = random.randint(1, 4)\n",
    "        for _ in range(N):\n",
    "            shape, color, dims = self.random_shape(height, width)\n",
    "            shapes.append((shape, color, dims))\n",
    "            x, y, s = dims\n",
    "            boxes.append([y-s, x-s, y+s, x+s])\n",
    "        # Apply non-max suppression wit 0.3 threshold to avoid\n",
    "        # shapes covering each other\n",
    "        keep_ixs = utils.non_max_suppression(np.array(boxes), np.arange(N), 0.3)\n",
    "        shapes = [s for i, s in enumerate(shapes) if i in keep_ixs]\n",
    "        return bg_color, shapes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Training dataset\n",
    "dataset_train = ShapesDataset()\n",
    "dataset_train.load_shapes(500, config.IMAGE_SHAPE[0], config.IMAGE_SHAPE[1])\n",
    "dataset_train.prepare()\n",
    "\n",
    "# Validation dataset\n",
    "dataset_val = ShapesDataset()\n",
    "dataset_val.load_shapes(50, config.IMAGE_SHAPE[0], config.IMAGE_SHAPE[1])\n",
    "dataset_val.prepare()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAy4AAACnCAYAAAD35AgmAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAACPdJREFUeJzt3WmopmUdx/Hf3xSRCjIik3oR064vSspWK4ukPWkvWrAFijJaiTbIdoqgAtuXKSqoiFIpQzLbxjTFgjaIxOpFWZMlZmVT6r8Xzz10mKaZozU+fzmfDxzOc1/nPvdzPcP14nzPdT9nqrsDAAAw2UHrngAAAMD+CBcAAGA84QIAAIwnXAAAgPGECwAAMJ5wAQAAxtsy4VJVt6+qs/cYu/h6XOfMqjpmefyoqrq8qmo5fldVPWsT13hLVf1643yq6piqOreqvlNV51TVtmV82zL2rar6ZlXdbh/XvUNVXVRVf6mq4zaMv7eqzl8+XrNh/LVVdWFVXVBVr7iu/xYAAHBD2TLh8n+0I8kDlscPSHJRkqM3HH93E9f4QJKH7DF2aZJHdPeDkrw7yZuW8Rcl+Xh3H5/kU0leso/rXprkhCRf3GP8/d193yT3T3LiEjg3T/LcJLvHX1hVN93E3NmCquom654DALC1CZc9VNUHqurZVXVQVZ1VVffZ45QdSXbvZtw9yQeTHFdVhyY5ort/tb/n6O5Lk1y7x9jvuvvK5XBXkquXxz9Ncovl8eFJdlbVoVW1o6ruWlW3WXZMDu/uv3X3n/byfL9YPl+7XPeaJFcl+W2Sw5aPq5L8c39zZ6aqOrqqzlt25b5WVUct6+KrVfWFqjplOe/iDd/zsao6fnl81rKrd0FV3W8ZO6WqPllVZyR5SlU9uKq+vZz3od07jQAAN4SD1z2BG9g9q+pb+znnFUnOyWr35Bvd/f09vn5Bkk9U1SFJOqsdlncn+UmSC5Nk+cHvHXu59pu7+5x9Pfmy6/HWJM9bhs5OclZVPS/JoUnu3d27luPtSa5I8rLuvnw/rytV9Ywkl+yOq6o6M8nPswrYt3b3P/Z3DcZ6eJLt3f2RqjooyZeTvLS7z6uqj27i+5/Q3X+tqrsleX+Shy7ju7r7cUuk/CDJ8d19RVW9J8mjk3zlALwWAID/sNXC5aLuftjug729x6W7/15V25O8K8mR/+XrO5M8IckPu3tnVd0mq12YHcs55yU5/rpObomhzyd5Z3f/bBl+Z5I3dPeXqurpSd6e5MXd/fOq+mWSW3b39zZx7YcleU6Sxy7Hd07yxCTbsgqXb1fVad39m+s6b0bYnuT1VfXZJD9KcqesIjtJvp9kb++N2v3erMOSvK+q7pLVbtxtN5yze23dKsntk5y+bLTcLKvohf9JVZ2c5ElJLu7u5697PmxN1iHrZg1uzlYLl/2qqiOz2u14S1aRsLc3re9I8uokr1uOf5vkyVmFwfXacVl+S/6ZJKd192kbv5TksuXxziS3XM4/IckhSS6rqsd19xn7eE33WV7PI7v7qg3XvbK7dy3n7Mrqh1FunHZ196uSZPmjD79Pcq+souXYrN7/lCRXLKH9hyT3SPLpJI9Ick13P7CqjkqycS1ds3y+LMklSR7T3X9ZnueQA/uS2Aq6+9Qkp657Hmxt1iHrZg1ujnDZYImH7VndenV+VX2uqh7V3WfuceqOJK9Mcv5yfG6SE7O6XWy/Oy5LVT8tyd2WHzJfkOSYrG69OaKqnpnkx939kqxuG/twVV2dVai8oKpuneRtWd0edHWSs6vqB0n+nORLSY5KcnRVndndb0zy8eWpT1t+W/7K7r5oeT/D+VlFzDe722/Qb7yeXlUnZXX74u+yWjcfq6o/5t/hm6x2Er+e1Xundi5j5yV57bIWz93bxbu7l788d8Zy29i1SV6e1e4OAMABV9297jkAB9ASwnfs7lPWPRcAgOvLXxUDAADGs+MCAACMZ8cFAAAYT7gAAADjjfirYk898ST3q20hnz/9kyP/x/XDjjnZOtxCrvrhqdYhazdxHVqDW8vENZhYh1vNZtehHRcAAGA84QIAAIwnXAAAgPGECwAAMJ5wAQAAxhMuAADAeMIFAAAYT7gAAADjCRcAAGA84QIAAIwnXAAAgPGECwAAMJ5wAQAAxhMuAADAeMIFAAAY7+B1T4C9u+TrX173FDZt2wmPX/cUOEAuv/DUdU9h0w4/9uR1TwEAOIDsuAAAAOMJFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYDzhAgAAjHfwuifA3m074fHrngLk8GNPXvcUAACS2HEBAABuBIQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYLzq7nXPAQAAYJ/suAAAAOMJFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy4AAMB4/wLdza0OzcvZJQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f075e81ca20>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAy4AAACnCAYAAAD35AgmAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAADZdJREFUeJzt3H2sZHddx/HPF1oL1IdulYcmkCAkqFQ0DVkRaNqiJUJVmlQ0GoFEUGtkm0BrBI1GBBSLqCQuGEIomigRo01DQhNIaYu0trBW/rCrVgmiUSgF3UDFdsvDzz/mDFyHu/dh79w7vzPn9Uo2uzN3euZ7m9Pbec/3zFZrLQAAAD172KoHAAAA2I5wAQAAuidcAACA7gkXAACge8IFAADonnABAAC6N5lwqaonVtVNC/d97DSOc2NVXTD8+bKqOlFVNdx+Y1W9eAfHeF1V/dvGearqgqq6var+uqpurqonDfc/abjv1qq6paoev8Vxn1xVd1XV/1TVhRvuf3NV3Tn8evWG+3+lqo5V1Ueq6urd/rtgHKrqcVX1e7t4/K1bnWewUVWdU1UvOcXX3lxVj17S83zdz3AApmUy4bJEtyV59vDnZye5K8n5G25/aAfHeGuS5yzc96kkz2utXZTkTUl+c7j/F5O8o7V2SZI/SXLVFsf9VJLnJvnLhfvf0lr7/iTPSnL5EDjflOSlSeb3/0JVnb2D2RmZ1tq9rbVrFu+vqoevYh7WzjlJvi5cqurhrbVXtNY+s4KZAFhDwmVBVb21ql5SVQ+rqvdV1TMWHnJbkvk243uT/FGSC6vqrCSPba19YrvnaK19KslXFu67t7V2/3DzZJIvDX8+ntkLgyQ5lOS+qjqrqm6rqu8c3k3/SFUdaq39b2vtvzd5vn8Zfv/KcNwvJ3kgySeTPHL49UCSL243O+NQVddW1R3Dlu7K+TvVVfWaqvrjqnpPkp+oqucMm75bq+oPNjnOG6rqg8OxfuTAvxHG4OokTx/OoWML59etVfX4qvq2qvrAcPv2qnpKkgyPfXtVvXfYCD9muP/qqvrbqvqz4ZhP3PiEVfWE4Z+5efh9KVsdAPp2xqoHOGBPr6pbt3nM1Uluzmx78oHW2ocXvv6RJNdV1ZlJWmYbljcluTvJsSSpqmcmecMmx35ta+3mrZ582Hq8PsnLhrtuSvK+qnpZkrOSfF9r7eRw+51JPpfkFa21E9t8X6mqn07y8XlcVdWNSe7JLGBf31p7aLtj0L+quizJE5I8q7XWqurJSX58w0NOttZeMFzi+I9JLm6tfXpxA1NVz0tyqLV2cVU9KskdVfXe1lo7qO+FUfj9JE9trV1aVa9Jcl5r7QVJUlVXDo/5XJLnt9YeqqrnJ3l1ZhvfJDneWvu5qvrVzGLnL5K8OMnhJI9K8vFNnvN3k7yutXZnVV2e5FVJfmmfvj8AOjG1cLmrtXbp/MZmn3FprT1YVe9M8sYk553i6/cluSLJR1tr91XV4zLbwtw2POaOJJfsdrghht6d5NrW2j8Md1+b5Ndaa9dX1U8l+e0kL2+t3VNV/5rk3Nba3+zg2Jcm+ZkkPzrcfkqSH0vypMzC5YNVdUNr7T93Ozfd+e4kt2wIjC8vfH1+vjw6yX+11j6dJK21xcc9LcnFG2L/rCTfmuSzS5+YdbLZz6Nzkrxl+Fn5DUnu3/C1u4bf/z3Jk5N8e5K7W2tfSvL5qvqnTY73tCS/M2vvnJFk159XhI2q6kiSFyb5WGvtZ1c9D9PjHNwZl4otqKrzMtt2vC6zSNjMbUl+Ocntw+1PZvaO9oeGYzxzuCRi8dcPbPG8D0vyp0luaK3dsPFL+doLxfuSnDs8/rlJzkzy2ap6wTbf0zOG7+eFrbUHNhz3/tbayeG+k0m+cavjMBp3J7l4w+3F/87ngfKZJOfOL7MZzsGNjid5f2vtkuEzVt/TWhMtLHoo//9NsMUATpIXZfZGz0VJXpvZz5+5jRu8SvKJJOdX1RnDZ/G+Y5PjHU/yyuHcvDDJz+9hfkhr7ehwPnnByEo4B3dmahuXLQ0v3N6Z2aVXd1bVn1fVZa21GxceeluSa5LcOdy+Pcnlmb1g3HbjMlT1Tyb5ruGzB1cmuSDJDyd5bFW9KMnft9auyuyysbdV1ZcyC5Urh+vAfyvJD2X2mZWbqurvknw+yfVJnprZ//hvbK39RpJ3DE99w/AO5TWttbuGz8bcmdmLhVtaa/ecxr82OtNau7GqLqmqOzL77NK7T/G4VlUvT/KeqjqZ5KNJXrlwnGcNG5eW5D8yu4QHNro3yQNV9VdJHpPNtx/vT/Kuqroos+g4peGyxXcl+XCSf87svHsos03N3DWZbXDmb7Zcl9kbPwCssXK5OgA9qaozW2tfrKpvziyon7LJpYwATIyNCwC9eXVV/WCSb0ny66IFgMTGBQAAGAEfzgcAALonXAAAgO518RmXP3zbta5Xm5CrrnxVbf+og/fIC444DyfkgY8edR6ycj2eh87BaenxHEych1Oz0/PQxgUAAOiecAEAALonXAAAgO4JFwAAoHvCBQAA6J5wAQAAuidcAACA7gkXAACge8IFAADonnABAAC6J1wAAIDuCRcAAKB7wgUAAOiecAEAALonXAAAgO4JFwAAoHvCBQAA6J5w2cTx689e9QiQE8eOrnoEAIBunLHqAVZluzjZ6uvnX/GFZY/DRG0XJ1t9/dDhI8seBwCgW5MLl2VsU+bHEDCcrmVsU+bHEDAAwBRM5lKx49efvfRLwFxSxm6dOHZ06ZeAuaQMAJiCtd+47Hdc2L6wE/sdF7YvAMC6W9twOehtiIBhMwe9DREwAMC6WstLxVZ5CZfLx5hb5SVcLh8DANbN2oVLD+HQwwysVg/h0MMMAADLslbh0lMw9DQLB6unYOhpFgCAvVibcOkxFHqcif3VYyj0OBMAwG6tRbj0HAg9z8Zy9RwIPc8GALATaxEuAADAeht9uIxhozGGGdmbMWw0xjAjAMCpjD5cAACA9SdcDoitCz2wdQEAxkq4AAAA3Rt1uIxtizG2edmZsW0xxjYvAEAy8nABAACmYbThMtbtxVjnZnNj3V6MdW4AYLpGGy4AAMB0CBcAAKB7wgUAAOiecAEAALonXAAAgO4JFwAAoHujDJex/5XCY5+fmbH/lcJjnx8AmJZRhsv5V3xh1SPsydjnZ+bQ4SOrHmFPxj4/ADAtowwXAABgWoQLAADQPeECAAB0T7gAAADdEy4AAED3RhsuY/2bucY6N5sb69/MNda5AYDpGm24AAAA0zHqcBnb9mJs87IzY9tejG1eAIBk5OECAABMg3A5ILYt9MC2BQAYK+ECAAB0b/ThMoZNxhhmZG/GsMkYw4wAAKcy+nABAADW31qES88bjZ5nY7l63mj0PBsAwE6sRbgkfQZCjzOxv3oMhB5nAgDYrbUJl6SvUOhpFg5WT6HQ0ywAAHuxVuGS9BEMPczAavUQDD3MAACwLGsXLslqw0G0MLfKcBAtAMC6OWPVA+yXeUAcv/7sA30+2GgeECeOHT3Q5wMAWDdrGy5z+x0wgoWd2O+AESwAwLpby0vFNnP+FV9YemSIFnbr0OEjS48M0QIATMHab1wWLWMDI1jYq2VsYAQLADAlkwuXua3i4/j1Z4sTDsRW8XHi2FFxAgAwmMylYrshWuiBaAEA+BrhAgAAdE+4AAAA3RMuAABA94QLAADQPeECAAB0T7gAAADdEy4AAED3hAsAANA94QIAAHRPuAAAAN0TLkvw0gcvXfUIAF04cezoqkcAYE0Jlz2aR4t4AaZuHi3iBYD9IFwAAIDuCZc9WNyy2LoAU7W4ZbF1AWDZhAsAANA94bJkti4AM7YuACyTcDlNAgVgRqAAcBCEyz4QNQAzogaAZREup2EnYSJegCnYSZiIFwCWQbgAAADdEy67tJtNiq0LsM52s0mxdQFgr4TLLggRgBkhAsBBEy77TOwAzIgdAPZCuOzQXgJEvADrZC8BIl4AOF3CBQAA6J5w2YFlbExsXYB1sIyNia0LAKdDuGxDcADMCA4AVkm4HCARBDAjggDYLeGyBaEBMCM0AFg14XLAxBDAjBgCYDeEyynsZ2CIF2BM9jMwxAsAOyVcNiEsAGaEBQC9EC4rIo4AZsQRADshXBYICoAZQQFAT4TLCokkgBmRBMB2hMsGQgJgRkgA0BvhMlhVtIgloDerihaxBMBWhEsHxAvAjHgB4FSES4QDwJxwAKBXkw+XXqKllzmA6eolWnqZA4C+TDpcXvrgpbnuETeteoyvEi/Aqpw4djSHDh9Z9RhfJV4AWDTpcEnEAsCcWACgZ5MNl16Dpde5gPXVa7D0OhcAqzHZcOmZeAGYES8AzE0yXIQBwIwwAGAsJhcuvX0g/1TEFbDfevtA/qmIKwCSCYZLIgoA5kQBAGMxqXAZW7CMbV5gPMYWLGObF4Dlm1S4jM11j7hJvAAkOXT4iHgBmLjJhMsYA2CMMwP9G2MAjHFmAJZrEuEylg/kb8bWBVimsXwgfzO2LgDTNolwSca7vRjr3EC/xvrif6xzA7Acax8uY962bCRggL0a87ZlIwEDME1rHy6JF/0Ac170AzBWax0u6xYs6/b9AAdn3YJl3b4fALZXrbVVzwAAALCltd64AAAA60G4AAAA3RMuAABA94QLAADQPeECAAB0T7gAAADdEy4AAED3hAsAANA94QIAAHRPuAAAAN0TLgAAQPeECwAA0D3hAgAAdE+4AAAA3RMuAABA94QLAADQPeECAAB0T7gAAADdEy4AAED3hAsAANA94QIAAHRPuAAAAN0TLgAAQPf+Dz+E4Sju1ZzFAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f075e6cde48>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAy4AAACnCAYAAAD35AgmAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAACotJREFUeJzt3X+IZfdZx/HPE5OG+jMp2jZQoaYQtUUklFo1pVRM0URs/I1iW9AqEU0hJFZjUfyRamxTJH+kFZE2ClpU+iMUDFTStOpGt40xf9iIqaVW0abGYmgjrknTPv5xz5bJZnZ3djI795m5rxcMc8+5Z8/5znAG5r3PvbvV3QEAAJjsnHUvAAAA4HSECwAAMJ5wAQAAxhMuAADAeMIFAAAYT7gAAADjbUy4VNVzq+rOE/Z9bBfnuaOqLl0eX1lVD1dVLdtvqqpX7eAcN1bVv25dT1VdWlV3V9VfVdVdVXXxsv/iZd8Hq+oDVfWcU5z3eVV1b1X9T1W9ZMv+W6rq6PJxw5b9v1RV91TVh6vqujP9XrBeVXVBVb36JM/dUlVfs0fXedLPDgDAftuYcNlDR5Jctjy+LMm9SV6wZfuvd3COtyb5jhP2PZjku7v7pUnenOTXl/0/m+Rt3f2yJH+Y5LWnOO+DSV6e5J0n7H9Ld39rkm9PctUSOF+R5CeTHN//M1X1ZTtYO3NckORJ4VJVX9Ld13b3f61hTQAAZ4VwOUFVvbWqXl1V51TV+6rqxSccciTJ8WnGNyf53SQvqarzkzyruz9xumt094NJvnDCvk919yPL5qNJHl8e35/VL6hJcmGSh6rq/Ko6UlXfUFXPXiYmF3b3/3b3f29zvX9ePn9hOe/nkxxL8skkT18+jiX53OnWzijXJXnhMo27p6r+oKrem+RHln3Pqaqvrqr3L9t3V9UlSbIc+/tV9efLJO6Zy/7rqurvquqPl3M+d+sFq+prlz9z1/J5T6Y6AACnc+66F7DPXlhVHzzNMdcluSur6cn7u/tDJzz/4SRvr6rzknRWE5Y3J/lIknuSpKq+LclN25z7N7r7rlNdfJl6vCHJa5ZddyZ5X1W9Jsn5Sb6lux9dtm9L8pkk13b3w6f5ulJVP57k48fjqqruSPJAVgH7hu5+7HTnYJTfSfL87r68qn4tyUXd/Yokqaqrl2M+k+SK7n6sqq5IckNWk7Ykub+7f7qqXp9V7PxZklcleVGSL03y8W2ueXOSG7v7aFVdleQXk/z8Wfr6AAC+aNPC5d7uvvz4xnbvcenu/6uq25K8KclFJ3n+oSQ/kOS+7n6oqp6d1RTmyHLM3yZ52ZkubomhP03yxu7+x2X3G5P8cne/u6p+LMlvJfm57n6gqv4lyTO6+292cO7Lk/xEku9dti9J8oNJLs4qXP6yqm7v7v8403Uzxnb3wQVJ3rLco09L8siW5+5dPv9bkucl+bokH+nux5N8tqr+aZvzfVOS317e1nVukjN+nxhsVVXXJPmhJB/r7p9a93rYTO5D1s09uDObFi6nVVUXZTXtuDGrSNjuTetHkvxCktcv259M8sNZhcGuJi5VdU6SP0pye3ffvvWpJJ9eHj+U5BnL8S9Pcl6ST1fVK7r7vaf4ml68fD1XdPexLed9pLsfXY55NMmXn+wcjPRYnvgz/PltjnllVoF9U1VdmSfez73lcSX5RJIXVNW5Wb188Ou3Od/9SW7q7vuSpKqetvvlQ9Ldtya5dd3rYLO5D1k39+DOCJctlni4LauXXh2tqj+pqiu7+44TDj2S5PokR5ftu5NcldXLxU47cVmq+keTfOPyrzVdneTSJN+T5FlV9cok/9Ddr83qZWO/V1WPZxUqVy/vR/jNJN+V1XtW7qyqv0/y2STvTvL8rH4BvaO7fzXJ25ZL3778Tfn13X3v8t6Yo1n90vqB7n5gF9821udTSY5V1buSPDPbTz/+Isk7quqlWUXHSXX3f1bVO5J8KMlHk/x7VnG0NU6uz2qCczxy355VcAMAnFXV3ac/CtgIVXVed3+uqr4yyX1JLunu7SY5AAD7ysQF2OqGqvrOJF+V5FdECwAwhYkLAAAwnv/HBQAAGE+4AAAA4414j8v3ff+1Xq+2QW5/zy217jVs5+mXXuM+3CDH7rvVfcjaTbwP3YObZeI9mLgPN81O70MTFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy6c1Ote99F1LwHy8D23rnsJAMAAwoVtHY8W8cI6HY8W8QIACBcAAGA84cKTnDhlMXVhHU6cspi6AMBmEy4AAMB4woUnONl0xdSF/XSy6YqpCwBsLuHCF4kTJhAnAMB2hAs7JmyYQNgAwGYSLiTZeZSIF86mnUaJeAGAzSNcECOMIEYAgFMRLgAAwHjCZcPtZtpiQsNe2820xYQGADaLcNlgAoQJBAgAsBPChV0RPUwgegBgcwiXDSU8mEB4AAA7JVzYNfHDBOIHADaDcNlAexkc4oXd2svgEC8AcPgJlw0jNJhAaAAAZ0q48JSJISYQQwBwuAmXDSIwmEBgAAC7IVzYE6KICUQRABxewmVDCAsmEBYAwG4Jlw2wX9EijjiV/YoWcQQAh5NwYU+JFyYQLwBw+AiXQ05IMIGQAACeKuHCnhNLTCCWAOBwES6HmIBgAgEBAOwF4XJIrTta1n19Zlh3tKz7+gDA3hEuAADAeMLlEJoy7ZiyDtZjyrRjyjoAgKdGuAAAAOMJl0Nm2pRj2nrYH9OmHNPWAwCcOeHCWSdemEC8AMDBJlwOEYHABAIBADgbhMshMT1apq+PvTE9WqavDwA4OeECAACMJ1wOgYMyzTgo62R3Dso046CsEwB4IuECAACMJ1wOuIM2xTho62VnDtoU46CtFwAQLgeaCGACEQAA7AfhckCJFiYQLQDAfhEuAADAeMIFAAAYT7gcQF4mxgReJgYA7CfhAgAAjCdcDhjTFiYwbQEA9ptwAQAAxhMuB4hpCxOYtgAA6yBcAACA8YTLAWHawgSmLQDAugiXA0C0MIFoAQDWSbgAAADjCZfhTFuYwLQFAFg34QIAAIwnXAYzbWEC0xYAYIJz170ATu7mmy9Z9xLOisves+4VcCYufNE1617CWXHsPkEGAAeJiQsAADCecAEAAMYTLgAAwHjCBQAAGE+4AAAA4wkXAABgPOECAACMJ1wAAIDxhAsAADCecAEAAMYTLgAAwHjCBQAAGE+4AAAA4wkXAABgPOECAACMJ1wAAIDxhAsAADCecAEAAMYTLgAAwHjCBQAAGE+4AAAA4wkXAABgPOECAACMJ1wAAIDxhAsAADCecAEAAMYTLgAAwHjCBQAAGE+4AAAA4wkXAABgPOECAACMJ1wAAIDxhAsAADCecAEAAMYTLgAAwHjCBQAAGE+4AAAA4wkXAABgPOECAACMJ1wAAIDxhAsAADCecAEAAMYTLgAAwHjCBQAAGE+4AAAA4wkXAABgPOECAACMJ1wAAIDxhAsAADCecAEAAMYTLgAAwHjCBQAAGE+4AAAA4wkXAABgPOECAACMJ1wAAIDxhAsAADCecAEAAMYTLgAAwHjCBQAAGE+4AAAA4wkXAABgvOruda8BAADglExcAACA8YQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYDzhAgAAjCdcAACA8YQLAAAwnnABAADGEy4AAMB4wgUAABhPuAAAAOMJFwAAYLz/B0txlpbqWU3LAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f075e61dc88>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAy4AAACnCAYAAAD35AgmAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAACktJREFUeJzt3F2oZXd5x/HfE5VgbaETtc5ACyWC+C6hqDVKE4tSiVVBbbH4AmohRiO0BqSFQl+0lYilvZgaemGjpRcVRILggBKTqJMmZoi5UFtCpSqoE2PaobVtnFb99+Ks3Ww2M+fNc85+lvvzgcOcvfaatZ89rGHOd/5r7RpjBAAAoLNL1j0AAADAToQLAADQnnABAADaEy4AAEB7wgUAAGhPuAAAAO1tTLhU1S9W1a0r2766j+Ocqqorpu+vqapzVVXT4/dX1Rt3cYz3VNU3luepqiuq6s6q+lxV3VZVl0/bL5+23VFVt1fVz29z3CdX1b1V9Z9V9aKl7X9ZVXdPX7+3tP33q+pMVd1TVe/a658F81BVx6vqz/ew/x3bnWcAAOuwMeFygE4neeH0/QuT3JvkGUuPP7+LY3wwyYtXtp1N8rIxxq8k+UCSP562vz3Jh8YYVyf5SJJ3bnPcs0lemuRjK9v/aozxy0muTPKqKXB+Jslbkiy2v62qHreL2ZmZMcYDY4wbVrdX1aPWMQ8AwH4IlxVV9cGqelNVXVJVn6qq56/scjrJYjXjOUluSvKiqro0yZPGGF/f6TXGGGeT/Ghl2wNjjO9ND88n+cH0/VeS/Oz0/bEkD1bVpVV1uqqeOv1v+j1VdWyM8d9jjH+7wOv98/Trj6bj/jDJw0m+neSx09fDSf53p9mZh6q6sarumlbprl2s7lXVH1XVh6vqE0l+s6pePK303VFVf3GB47yvqj47HevXj/yNAABMHr3uAY7YL1XVHTvs864kt2Vr9eQzY4wvrDx/T5K/qarHJBnZWmH5QJIvJzmTJFX1giTvu8Cx/2SMcdt2Lz6terw3yVunTbcm+VRVvTXJpUmeN8Y4Pz2+Ocm/J/mdMca5Hd5Xqur1Sf5lEVdVdSrJ/dkK2PeOMf5np2PQX1Vdk+QXklw5xhhV9eQkv7G0y/kxxiunSxz/KclVY4zvrK7AVNXLkhwbY1xVVT+V5K6q+uQYYxzVewEAWNi0cLl3jPGSxYML3eMyxvh+Vd2c5P1JTlzk+QeTvDrJfWOMB6vqeLZWYU5P+9yV5Oq9DjfF0EeT3DjG+Mdp841J/mCM8fGq+q0kf5bkHWOM+6vqa0kuG2P8wy6O/ZIkb07yiunxU5K8Jsnl2QqXz1bVLWOMb+11btp5ZpLblwLjhyvPL86XJyb51zHGd5JkjLG637OSXLUU+5cmeXyShw58YjZWVV2f5LVJvjrG+O11z8Nmch6ybs7B3XGp2IqqOpGt1Y73ZCsSLuR0kncnuXN6/O1s/Y/256djvGC69Gb161e3ed1LkvxdklvGGLcsP5VHflB8MMll0/4vTfKYJA9V1St3eE/Pn97Pa8cYDy8d93tjjPPTtvNJfnq74zAbX05y1dLj1b/ni0D5bpLLquqJyf+fg8u+kuTTY4yrp3usnj3GEC0cqDHGyekc8w81a+M8ZN2cg7uzaSsu25p+cLs5W5de3V1Vf19V14wxTq3sejrJDUnunh7fmeRV2fqBcccVl6mqX5fkadO9B9cmuSLJy5M8qarekORLY4x3Zuuysb+uqh9kK1SuraqfS/KnSX4tW/es3FpVX0zyH0k+nuTpSZ5RVafGGH+Y5EPTS98yfQDaDWOMe6d7Y+7OVsTcPsa4fx9/bDQzxjhVVVdX1V3ZunfpoxfZb1TVO5J8oqrOJ7kvye+uHOfKacVlJPlmkh0/NQ8A4DCUy9UBAIDuXCoGAAC0J1wAAID2hAsAANCecAEAANpr8alib3rC43xCwAb524f+q9Y9w4U89orrnYcb5OH7TjoPWbuO56FzcLN0PAcT5+Gm2e15aMUFAABoT7gAAADtCRcAAKC9jQqXE8evW/cIkHNnTq57BACA2Wlxc/5B2y5QLvbc2QduOqxx2FDbBcrFnjv23OsPaxwAgFn7iQmXH3c1Zfn3ixj268ddTVn+/SIGAOARsw+Xw7j8a3FMAcNuHcblX4tjChgAgA27x2Wv3BNDB+6JAQCY8YrLUUWF1Re2c1RRYfUFANh0s1xxWcdKiNUXVq1jJcTqCwCwqWYXLusMCPHCwjoDQrwAAJtoVuHSIRw6zMB6dQiHDjMAAByl2YRLp2DoNAtHq1MwdJoFAOCwzSZcAACAzdU+XE4cv67lCkfHmTg8586cbLnC0XEmAIDD0D5cAAAAWodL91WN7vNxMLqvanSfDwDgILQOFwAAgES4AAAAM9A2XOZyGdZc5mR/5nIZ1lzmBADYr7bhAgAAsCBcAACA9oQLAADQnnABAADaEy4AAEB7LcNlbp/UNbd52Z25fVLX3OYFANiLluFy9oGb1j3CnsxtXnbn2HOvX/cIezK3eQEA9qJluAAAACwTLgAAQHvCBQAAaE+4AAAA7QkXAACgvbbhMpdP6prLnOzPXD6pay5zAgDsV9twAQAAWBAuAABAe63DpftlWN3n42B0vwyr+3wAAAehdbgAAAAkMwgXqxp0YFUDAGC92odL0jNeOs7E4eoYLx1nAgA4DLMIFwAAYLPNJlw6rXB0moWj1WmFo9MsAACHbTbhkvQIhg4zsF4dgqHDDAAAR2lW4ZKsNxxECwvrDAfRAgBsotmFS7KegBAtrFpHQIgWAGBTPXrdA+zXIiROHL/uSF4HLmQREufOnDyS1wEA2FSzDZeFwwgYscJeHUbAiBUAgEfMPlwWlmNjPxEjVjgIy7Gxn4gRKwAAF/YTEy7LLhYhJ45fJ1A4MheLkHNnTgoUAIA9muXN+fslWuhAtAAA7N1GhQsAADBPwgUAAGhPuAAAAO0JFwAAoD3hAgAAtCdcAACA9oQLAADQnnABAADaEy4AAEB7wgUAAGhPuAAAAO0JFwAAoD3hAgAAtCdcAACA9oQLAADQnnABAADaEy4AAEB7wgUAAGhPuAAAAO0JFwAAoD3hAgAAtCdcAACA9oQLAADQnnABAADaEy4AAEB7wgUAAGhPuAAAAO0JFwAAoD3hAgAAtCdcAACA9oQLAADQnnABAADaEy4AAEB7wgUAAGhPuAAAAO0JFwAAoD3hAgAAtCdcAACA9oQLAADQnnABAADaEy4AAEB7wgUAAGhPuAAAAO0JFwAAoD3hAgAAtCdcAACA9oQLAADQnnABAADaEy4AAEB7wgUAAGhPuAAAAO0JFwAAoD3hAgAAtCdcAACA9oQLAADQnnABAADaEy4AAEB7wgUAAGhPuAAAAO0JFwAAoD3hAgAAtCdcAACA9oQLAADQnnABAADaEy4AAEB7wgUAAGhPuAAAAO0JFwAAoD3hAgAAtCdcAACA9oQLAADQnnABAADaEy4AAEB7wgUAAGivxhjrngEAAGBbVlwAAID2hAsAANCecAEAANoTLgAAQHvCBQAAaE+4AAAA7QkXAACgPeECAAC0J1wAAID2hAsAANCecAEAANoTLgAAQHvCBQAAaE+4AAAA7QkXAACgPeECAAC0J1wAAID2hAsAANCecAEAANoTLgAAQHvCBQAAaE+4AAAA7QkXAACgvf8DtzFsma1CmoQAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f075e4cf278>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Load and display random samples\n",
    "image_ids = np.random.choice(dataset_train.image_ids, 4)\n",
    "for image_id in image_ids:\n",
    "    image = dataset_train.load_image(image_id)\n",
    "    mask, class_ids = dataset_train.load_mask(image_id)\n",
    "    visualize.display_top_masks(image, mask, class_ids, dataset_train.class_names)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Ceate Model"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create model in training mode\n",
    "model = modellib.MaskRCNN(mode=\"training\", config=config,\n",
    "                          model_dir=MODEL_DIR)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "scrolled": false
   },
   "outputs": [],
   "source": [
    "# Which weights to start with?\n",
    "init_with = \"coco\"  # imagenet, coco, or last\n",
    "\n",
    "if init_with == \"imagenet\":\n",
    "    model.load_weights(model.get_imagenet_weights(), by_name=True)\n",
    "elif init_with == \"coco\":\n",
    "    # Load weights trained on MS COCO, but skip layers that\n",
    "    # are different due to the different number of classes\n",
    "    # See README for instructions to download the COCO weights\n",
    "    model.load_weights(COCO_MODEL_PATH, by_name=True,\n",
    "                       exclude=[\"mrcnn_class_logits\", \"mrcnn_bbox_fc\", \n",
    "                                \"mrcnn_bbox\", \"mrcnn_mask\"])\n",
    "elif init_with == \"last\":\n",
    "    # Load the last model you trained and continue training\n",
    "    model.load_weights(model.find_last()[1], by_name=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Training\n",
    "\n",
    "Train in two stages:\n",
    "1. Only the heads. Here we're freezing all the backbone layers and training only the randomly initialized layers (i.e. the ones that we didn't use pre-trained weights from MS COCO). To train only the head layers, pass `layers='heads'` to the `train()` function.\n",
    "\n",
    "2. Fine-tune all layers. For this simple example it's not necessary, but we're including it to show the process. Simply pass `layers=\"all` to train all layers."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "Starting at epoch 0. LR=0.001\n",
      "\n",
      "Checkpoint Path: /home/chuqh/Documents/idrid-2/logs/shapes20180402T1637/mask_rcnn_shapes_{epoch:04d}.h5\n",
      "Selecting layers to train\n",
      "fpn_c5p5               (Conv2D)\n",
      "fpn_c4p4               (Conv2D)\n",
      "fpn_c3p3               (Conv2D)\n",
      "fpn_c2p2               (Conv2D)\n",
      "fpn_p5                 (Conv2D)\n",
      "fpn_p2                 (Conv2D)\n",
      "fpn_p3                 (Conv2D)\n",
      "fpn_p4                 (Conv2D)\n",
      "In model:  rpn_model\n",
      "    rpn_conv_shared        (Conv2D)\n",
      "    rpn_class_raw          (Conv2D)\n",
      "    rpn_bbox_pred          (Conv2D)\n",
      "mrcnn_mask_conv1       (TimeDistributed)\n",
      "mrcnn_mask_bn1         (TimeDistributed)\n",
      "mrcnn_mask_conv2       (TimeDistributed)\n",
      "mrcnn_mask_bn2         (TimeDistributed)\n",
      "mrcnn_class_conv1      (TimeDistributed)\n",
      "mrcnn_class_bn1        (TimeDistributed)\n",
      "mrcnn_mask_conv3       (TimeDistributed)\n",
      "mrcnn_mask_bn3         (TimeDistributed)\n",
      "mrcnn_class_conv2      (TimeDistributed)\n",
      "mrcnn_class_bn2        (TimeDistributed)\n",
      "mrcnn_mask_conv4       (TimeDistributed)\n",
      "mrcnn_mask_bn4         (TimeDistributed)\n",
      "mrcnn_bbox_fc          (TimeDistributed)\n",
      "mrcnn_mask_deconv      (TimeDistributed)\n",
      "mrcnn_class_logits     (TimeDistributed)\n",
      "mrcnn_mask             (TimeDistributed)\n",
      "WARNING:tensorflow:From /home/chuqh/Documents/idrid-2/model.py:2078: calling reduce_mean (from tensorflow.python.ops.math_ops) with keep_dims is deprecated and will be removed in a future version.\n",
      "Instructions for updating:\n",
      "keep_dims is deprecated, use keepdims instead\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/usr/local/lib/python3.5/dist-packages/tensorflow/python/ops/gradients_impl.py:97: UserWarning: Converting sparse IndexedSlices to a dense Tensor of unknown shape. This may consume a large amount of memory.\n",
      "  \"Converting sparse IndexedSlices to a dense Tensor of unknown shape. \"\n",
      "/usr/local/lib/python3.5/dist-packages/keras/engine/training.py:2095: UserWarning: Using a generator with `use_multiprocessing=True` and multiple workers may duplicate your data. Please consider using the`keras.utils.Sequence class.\n",
      "  UserWarning('Using a generator with `use_multiprocessing=True`'\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Epoch 1/1\n",
      " 99/100 [============================>.] - ETA: 0s - loss: 1.8654 - rpn_class_loss: 0.0306 - rpn_bbox_loss: 0.6004 - mrcnn_class_loss: 0.4111 - mrcnn_bbox_loss: 0.3997 - mrcnn_mask_loss: 0.4236"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Process Process-5:\n",
      "Process Process-2:\n",
      "Process Process-4:\n",
      "Process Process-3:\n",
      "Traceback (most recent call last):\n",
      "Traceback (most recent call last):\n",
      "Traceback (most recent call last):\n",
      "Traceback (most recent call last):\n",
      "  File \"/usr/lib/python3.5/multiprocessing/process.py\", line 249, in _bootstrap\n",
      "    self.run()\n",
      "  File \"/usr/lib/python3.5/multiprocessing/process.py\", line 249, in _bootstrap\n",
      "    self.run()\n",
      "  File \"/usr/lib/python3.5/multiprocessing/process.py\", line 249, in _bootstrap\n",
      "    self.run()\n",
      "  File \"/usr/lib/python3.5/multiprocessing/process.py\", line 93, in run\n",
      "    self._target(*self._args, **self._kwargs)\n",
      "  File \"/usr/lib/python3.5/multiprocessing/process.py\", line 93, in run\n",
      "    self._target(*self._args, **self._kwargs)\n",
      "  File \"/usr/lib/python3.5/multiprocessing/process.py\", line 93, in run\n",
      "    self._target(*self._args, **self._kwargs)\n",
      "  File \"/usr/lib/python3.5/multiprocessing/process.py\", line 249, in _bootstrap\n",
      "    self.run()\n",
      "  File \"/usr/local/lib/python3.5/dist-packages/keras/utils/data_utils.py\", line 664, in _data_generator_task\n",
      "    self.queue.put((True, generator_output))\n",
      "  File \"/usr/lib/python3.5/multiprocessing/process.py\", line 93, in run\n",
      "    self._target(*self._args, **self._kwargs)\n",
      "  File \"/usr/local/lib/python3.5/dist-packages/keras/utils/data_utils.py\", line 664, in _data_generator_task\n",
      "    self.queue.put((True, generator_output))\n",
      "  File \"/usr/local/lib/python3.5/dist-packages/keras/utils/data_utils.py\", line 664, in _data_generator_task\n",
      "    self.queue.put((True, generator_output))\n",
      "  File \"<string>\", line 2, in put\n",
      "  File \"<string>\", line 2, in put\n",
      "  File \"<string>\", line 2, in put\n",
      "  File \"/usr/local/lib/python3.5/dist-packages/keras/utils/data_utils.py\", line 662, in _data_generator_task\n",
      "    self.queue.qsize() < self.max_queue_size):\n",
      "  File \"/usr/lib/python3.5/multiprocessing/managers.py\", line 717, in _callmethod\n",
      "    kind, result = conn.recv()\n",
      "  File \"/usr/lib/python3.5/multiprocessing/managers.py\", line 717, in _callmethod\n",
      "    kind, result = conn.recv()\n",
      "  File \"/usr/lib/python3.5/multiprocessing/connection.py\", line 250, in recv\n",
      "    buf = self._recv_bytes()\n",
      "  File \"/usr/lib/python3.5/multiprocessing/connection.py\", line 407, in _recv_bytes\n",
      "    buf = self._recv(4)\n",
      "  File \"<string>\", line 2, in qsize\n",
      "  File \"/usr/lib/python3.5/multiprocessing/connection.py\", line 250, in recv\n",
      "    buf = self._recv_bytes()\n",
      "  File \"/usr/lib/python3.5/multiprocessing/managers.py\", line 717, in _callmethod\n",
      "    kind, result = conn.recv()\n",
      "  File \"/usr/lib/python3.5/multiprocessing/managers.py\", line 717, in _callmethod\n",
      "    kind, result = conn.recv()\n",
      "  File \"/usr/lib/python3.5/multiprocessing/connection.py\", line 250, in recv\n",
      "    buf = self._recv_bytes()\n",
      "  File \"/usr/lib/python3.5/multiprocessing/connection.py\", line 379, in _recv\n",
      "    chunk = read(handle, remaining)\n",
      "  File \"/usr/lib/python3.5/multiprocessing/connection.py\", line 250, in recv\n",
      "    buf = self._recv_bytes()\n",
      "  File \"/usr/lib/python3.5/multiprocessing/connection.py\", line 407, in _recv_bytes\n",
      "    buf = self._recv(4)\n",
      "  File \"/usr/lib/python3.5/multiprocessing/connection.py\", line 379, in _recv\n",
      "    chunk = read(handle, remaining)\n",
      "KeyboardInterrupt\n",
      "  File \"/usr/lib/python3.5/multiprocessing/connection.py\", line 407, in _recv_bytes\n",
      "    buf = self._recv(4)\n",
      "  File \"/usr/lib/python3.5/multiprocessing/connection.py\", line 379, in _recv\n",
      "    chunk = read(handle, remaining)\n",
      "  File \"/usr/lib/python3.5/multiprocessing/connection.py\", line 407, in _recv_bytes\n",
      "    buf = self._recv(4)\n",
      "KeyboardInterrupt\n",
      "  File \"/usr/lib/python3.5/multiprocessing/connection.py\", line 379, in _recv\n",
      "    chunk = read(handle, remaining)\n",
      "KeyboardInterrupt\n",
      "KeyboardInterrupt\n"
     ]
    },
    {
     "ename": "KeyboardInterrupt",
     "evalue": "",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mKeyboardInterrupt\u001b[0m                         Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-9-83fb3ae74319>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m      6\u001b[0m             \u001b[0mlearning_rate\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mconfig\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mLEARNING_RATE\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m      7\u001b[0m             \u001b[0mepochs\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;36m1\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 8\u001b[0;31m             layers='heads')\n\u001b[0m",
      "\u001b[0;32m~/Documents/idrid-2/model.py\u001b[0m in \u001b[0;36mtrain\u001b[0;34m(self, train_dataset, val_dataset, learning_rate, epochs, layers)\u001b[0m\n\u001b[1;32m   2241\u001b[0m             \u001b[0mmax_queue_size\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;36m100\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   2242\u001b[0m             \u001b[0mworkers\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mworkers\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 2243\u001b[0;31m             \u001b[0muse_multiprocessing\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mTrue\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   2244\u001b[0m         )\n\u001b[1;32m   2245\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mepoch\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mmax\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mepoch\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mepochs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/local/lib/python3.5/dist-packages/keras/legacy/interfaces.py\u001b[0m in \u001b[0;36mwrapper\u001b[0;34m(*args, **kwargs)\u001b[0m\n\u001b[1;32m     89\u001b[0m                 warnings.warn('Update your `' + object_name +\n\u001b[1;32m     90\u001b[0m                               '` call to the Keras 2 API: ' + signature, stacklevel=2)\n\u001b[0;32m---> 91\u001b[0;31m             \u001b[0;32mreturn\u001b[0m \u001b[0mfunc\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     92\u001b[0m         \u001b[0mwrapper\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_original_function\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfunc\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     93\u001b[0m         \u001b[0;32mreturn\u001b[0m \u001b[0mwrapper\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/local/lib/python3.5/dist-packages/keras/engine/training.py\u001b[0m in \u001b[0;36mfit_generator\u001b[0;34m(self, generator, steps_per_epoch, epochs, verbose, callbacks, validation_data, validation_steps, class_weight, max_queue_size, workers, use_multiprocessing, shuffle, initial_epoch)\u001b[0m\n\u001b[1;32m   2268\u001b[0m                                 \u001b[0mbatch_size\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mbatch_size\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   2269\u001b[0m                                 \u001b[0msample_weight\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mval_sample_weights\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 2270\u001b[0;31m                                 verbose=0)\n\u001b[0m\u001b[1;32m   2271\u001b[0m                         \u001b[0;32mif\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0misinstance\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mval_outs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mlist\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   2272\u001b[0m                             \u001b[0mval_outs\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0mval_outs\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/local/lib/python3.5/dist-packages/keras/engine/training.py\u001b[0m in \u001b[0;36mevaluate\u001b[0;34m(self, x, y, batch_size, verbose, sample_weight, steps)\u001b[0m\n\u001b[1;32m   1784\u001b[0m                                \u001b[0mbatch_size\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mbatch_size\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1785\u001b[0m                                \u001b[0mverbose\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mverbose\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1786\u001b[0;31m                                steps=steps)\n\u001b[0m\u001b[1;32m   1787\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1788\u001b[0m     def predict(self, x,\n",
      "\u001b[0;32m/usr/local/lib/python3.5/dist-packages/keras/engine/training.py\u001b[0m in \u001b[0;36m_test_loop\u001b[0;34m(self, f, ins, batch_size, verbose, steps)\u001b[0m\n\u001b[1;32m   1431\u001b[0m                     \u001b[0mins_batch\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mi\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mins_batch\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mi\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mtoarray\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1432\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1433\u001b[0;31m                 \u001b[0mbatch_outs\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mf\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mins_batch\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1434\u001b[0m                 \u001b[0;32mif\u001b[0m \u001b[0misinstance\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mbatch_outs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mlist\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1435\u001b[0m                     \u001b[0;32mif\u001b[0m \u001b[0mbatch_index\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0;36m0\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/local/lib/python3.5/dist-packages/keras/backend/tensorflow_backend.py\u001b[0m in \u001b[0;36m__call__\u001b[0;34m(self, inputs)\u001b[0m\n\u001b[1;32m   2473\u001b[0m         \u001b[0msession\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mget_session\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   2474\u001b[0m         updated = session.run(fetches=fetches, feed_dict=feed_dict,\n\u001b[0;32m-> 2475\u001b[0;31m                               **self.session_kwargs)\n\u001b[0m\u001b[1;32m   2476\u001b[0m         \u001b[0;32mreturn\u001b[0m \u001b[0mupdated\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0mlen\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0moutputs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   2477\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/local/lib/python3.5/dist-packages/tensorflow/python/client/session.py\u001b[0m in \u001b[0;36mrun\u001b[0;34m(self, fetches, feed_dict, options, run_metadata)\u001b[0m\n\u001b[1;32m    893\u001b[0m     \u001b[0;32mtry\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    894\u001b[0m       result = self._run(None, fetches, feed_dict, options_ptr,\n\u001b[0;32m--> 895\u001b[0;31m                          run_metadata_ptr)\n\u001b[0m\u001b[1;32m    896\u001b[0m       \u001b[0;32mif\u001b[0m \u001b[0mrun_metadata\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    897\u001b[0m         \u001b[0mproto_data\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mtf_session\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mTF_GetBuffer\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mrun_metadata_ptr\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/local/lib/python3.5/dist-packages/tensorflow/python/client/session.py\u001b[0m in \u001b[0;36m_run\u001b[0;34m(self, handle, fetches, feed_dict, options, run_metadata)\u001b[0m\n\u001b[1;32m   1126\u001b[0m     \u001b[0;32mif\u001b[0m \u001b[0mfinal_fetches\u001b[0m \u001b[0;32mor\u001b[0m \u001b[0mfinal_targets\u001b[0m \u001b[0;32mor\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0mhandle\u001b[0m \u001b[0;32mand\u001b[0m \u001b[0mfeed_dict_tensor\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1127\u001b[0m       results = self._do_run(handle, final_targets, final_fetches,\n\u001b[0;32m-> 1128\u001b[0;31m                              feed_dict_tensor, options, run_metadata)\n\u001b[0m\u001b[1;32m   1129\u001b[0m     \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1130\u001b[0m       \u001b[0mresults\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/local/lib/python3.5/dist-packages/tensorflow/python/client/session.py\u001b[0m in \u001b[0;36m_do_run\u001b[0;34m(self, handle, target_list, fetch_list, feed_dict, options, run_metadata)\u001b[0m\n\u001b[1;32m   1342\u001b[0m     \u001b[0;32mif\u001b[0m \u001b[0mhandle\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1343\u001b[0m       return self._do_call(_run_fn, self._session, feeds, fetches, targets,\n\u001b[0;32m-> 1344\u001b[0;31m                            options, run_metadata)\n\u001b[0m\u001b[1;32m   1345\u001b[0m     \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1346\u001b[0m       \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_do_call\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0m_prun_fn\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_session\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mhandle\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfeeds\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfetches\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/local/lib/python3.5/dist-packages/tensorflow/python/client/session.py\u001b[0m in \u001b[0;36m_do_call\u001b[0;34m(self, fn, *args)\u001b[0m\n\u001b[1;32m   1348\u001b[0m   \u001b[0;32mdef\u001b[0m \u001b[0m_do_call\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfn\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m*\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1349\u001b[0m     \u001b[0;32mtry\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1350\u001b[0;31m       \u001b[0;32mreturn\u001b[0m \u001b[0mfn\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1351\u001b[0m     \u001b[0;32mexcept\u001b[0m \u001b[0merrors\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mOpError\u001b[0m \u001b[0;32mas\u001b[0m \u001b[0me\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1352\u001b[0m       \u001b[0mmessage\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mcompat\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mas_text\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0me\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmessage\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/local/lib/python3.5/dist-packages/tensorflow/python/client/session.py\u001b[0m in \u001b[0;36m_run_fn\u001b[0;34m(session, feed_dict, fetch_list, target_list, options, run_metadata)\u001b[0m\n\u001b[1;32m   1327\u001b[0m           return tf_session.TF_Run(session, options,\n\u001b[1;32m   1328\u001b[0m                                    \u001b[0mfeed_dict\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfetch_list\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mtarget_list\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1329\u001b[0;31m                                    status, run_metadata)\n\u001b[0m\u001b[1;32m   1330\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1331\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m_prun_fn\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0msession\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mhandle\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfeed_dict\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfetch_list\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;31mKeyboardInterrupt\u001b[0m: "
     ]
    }
   ],
   "source": [
    "# Train the head branches\n",
    "# Passing layers=\"heads\" freezes all layers except the head\n",
    "# layers. You can also pass a regular expression to select\n",
    "# which layers to train by name pattern.\n",
    "model.train(dataset_train, dataset_val, \n",
    "            learning_rate=config.LEARNING_RATE, \n",
    "            epochs=1, \n",
    "            layers='heads')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "scrolled": false
   },
   "outputs": [],
   "source": [
    "# Fine tune all layers\n",
    "# Passing layers=\"all\" trains all layers. You can also \n",
    "# pass a regular expression to select which layers to\n",
    "# train by name pattern.\n",
    "model.train(dataset_train, dataset_val, \n",
    "            learning_rate=config.LEARNING_RATE / 10,\n",
    "            epochs=2, \n",
    "            layers=\"all\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Save weights\n",
    "# Typically not needed because callbacks save after every epoch\n",
    "# Uncomment to save manually\n",
    "# model_path = os.path.join(MODEL_DIR, \"mask_rcnn_shapes.h5\")\n",
    "# model.keras_model.save_weights(model_path)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Detection"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "class InferenceConfig(ShapesConfig):\n",
    "    GPU_COUNT = 1\n",
    "    IMAGES_PER_GPU = 1\n",
    "\n",
    "inference_config = InferenceConfig()\n",
    "\n",
    "# Recreate the model in inference mode\n",
    "model = modellib.MaskRCNN(mode=\"inference\", \n",
    "                          config=inference_config,\n",
    "                          model_dir=MODEL_DIR)\n",
    "\n",
    "# Get path to saved weights\n",
    "# Either set a specific path or find last trained weights\n",
    "# model_path = os.path.join(ROOT_DIR, \".h5 file name here\")\n",
    "model_path = model.find_last()[1]\n",
    "\n",
    "# Load trained weights (fill in path to trained weights here)\n",
    "assert model_path != \"\", \"Provide path to trained weights\"\n",
    "print(\"Loading weights from \", model_path)\n",
    "model.load_weights(model_path, by_name=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Test on a random image\n",
    "image_id = random.choice(dataset_val.image_ids)\n",
    "original_image, image_meta, gt_class_id, gt_bbox, gt_mask =\\\n",
    "    modellib.load_image_gt(dataset_val, inference_config, \n",
    "                           image_id, use_mini_mask=False)\n",
    "\n",
    "log(\"original_image\", original_image)\n",
    "log(\"image_meta\", image_meta)\n",
    "log(\"gt_class_id\", gt_class_id)\n",
    "log(\"gt_bbox\", gt_bbox)\n",
    "log(\"gt_mask\", gt_mask)\n",
    "\n",
    "visualize.display_instances(original_image, gt_bbox, gt_mask, gt_class_id, \n",
    "                            dataset_train.class_names, figsize=(8, 8))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "results = model.detect([original_image], verbose=1)\n",
    "\n",
    "r = results[0]\n",
    "visualize.display_instances(original_image, r['rois'], r['masks'], r['class_ids'], \n",
    "                            dataset_val.class_names, r['scores'], ax=get_ax())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Evaluation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Compute VOC-Style mAP @ IoU=0.5\n",
    "# Running on 10 images. Increase for better accuracy.\n",
    "image_ids = np.random.choice(dataset_val.image_ids, 10)\n",
    "APs = []\n",
    "for image_id in image_ids:\n",
    "    # Load image and ground truth data\n",
    "    image, image_meta, gt_class_id, gt_bbox, gt_mask =\\\n",
    "        modellib.load_image_gt(dataset_val, inference_config,\n",
    "                               image_id, use_mini_mask=False)\n",
    "    molded_images = np.expand_dims(modellib.mold_image(image, inference_config), 0)\n",
    "    # Run object detection\n",
    "    results = model.detect([image], verbose=0)\n",
    "    r = results[0]\n",
    "    # Compute AP\n",
    "    AP, precisions, recalls, overlaps =\\\n",
    "        utils.compute_ap(gt_bbox, gt_class_id, gt_mask,\n",
    "                         r[\"rois\"], r[\"class_ids\"], r[\"scores\"], r['masks'])\n",
    "    APs.append(AP)\n",
    "    \n",
    "print(\"mAP: \", np.mean(APs))"
   ]
  },
  {
   "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.5.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}