Validation_and_Classification_Maps.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Define Dependencies"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Using TensorFlow backend.\n",
      "/anaconda3/envs/benchmarking/lib/python3.6/site-packages/tensorflow/python/framework/dtypes.py:526: FutureWarning: Passing (type, 1) or '1type' as a synonym of type is deprecated; in a future version of numpy, it will be understood as (type, (1,)) / '(1,)type'.\n",
      "  _np_qint8 = np.dtype([(\"qint8\", np.int8, 1)])\n",
      "/anaconda3/envs/benchmarking/lib/python3.6/site-packages/tensorflow/python/framework/dtypes.py:527: FutureWarning: Passing (type, 1) or '1type' as a synonym of type is deprecated; in a future version of numpy, it will be understood as (type, (1,)) / '(1,)type'.\n",
      "  _np_quint8 = np.dtype([(\"quint8\", np.uint8, 1)])\n",
      "/anaconda3/envs/benchmarking/lib/python3.6/site-packages/tensorflow/python/framework/dtypes.py:528: FutureWarning: Passing (type, 1) or '1type' as a synonym of type is deprecated; in a future version of numpy, it will be understood as (type, (1,)) / '(1,)type'.\n",
      "  _np_qint16 = np.dtype([(\"qint16\", np.int16, 1)])\n",
      "/anaconda3/envs/benchmarking/lib/python3.6/site-packages/tensorflow/python/framework/dtypes.py:529: FutureWarning: Passing (type, 1) or '1type' as a synonym of type is deprecated; in a future version of numpy, it will be understood as (type, (1,)) / '(1,)type'.\n",
      "  _np_quint16 = np.dtype([(\"quint16\", np.uint16, 1)])\n",
      "/anaconda3/envs/benchmarking/lib/python3.6/site-packages/tensorflow/python/framework/dtypes.py:530: FutureWarning: Passing (type, 1) or '1type' as a synonym of type is deprecated; in a future version of numpy, it will be understood as (type, (1,)) / '(1,)type'.\n",
      "  _np_qint32 = np.dtype([(\"qint32\", np.int32, 1)])\n",
      "/anaconda3/envs/benchmarking/lib/python3.6/site-packages/tensorflow/python/framework/dtypes.py:535: FutureWarning: Passing (type, 1) or '1type' as a synonym of type is deprecated; in a future version of numpy, it will be understood as (type, (1,)) / '(1,)type'.\n",
      "  np_resource = np.dtype([(\"resource\", np.ubyte, 1)])\n"
     ]
    }
   ],
   "source": [
    "# Import the necessary libraries\n",
    "from sklearn.decomposition import PCA\n",
    "import os\n",
    "import scipy.io as sio\n",
    "import numpy as np\n",
    "from keras.models import load_model\n",
    "from keras.utils import np_utils\n",
    "from sklearn.metrics import classification_report, confusion_matrix\n",
    "import itertools\n",
    "import spectral\n",
    "import matplotlib\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Load the Global values (windowSize, numPCAcomponents, testRatio) from the text file global_variables.txt\n",
    "myFile = open('global_variables.txt', 'r') \n",
    "file = myFile.readlines()[:]\n",
    "\n",
    "\n",
    "for line in file:\n",
    "\n",
    "    if line[0:3] == \"win\":\n",
    "\n",
    "        ds = line.find('=')\n",
    "        windowSize = int(line[ds+1:-1],10)\n",
    "\n",
    "    elif line[0:3] == \"num\":\n",
    "\n",
    "        ds = line.find('=')\n",
    "        numPCAcomponents = int(line[ds+2:-1],10)\n",
    "\n",
    "    else:\n",
    "\n",
    "        ds = line.find('=')\n",
    "        testRatio = float(line[ds+1:])\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Define the neccesary functions for later use"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [],
   "source": [
    "def loadIndianPinesData():\n",
    "    data_path = os.path.join(os.getcwd(),'Data')\n",
    "    data = sio.loadmat(os.path.join(data_path, 'Indian_pines_corrected.mat'))['indian_pines_corrected']\n",
    "    labels = sio.loadmat(os.path.join(data_path, 'Indian_pines_gt.mat'))['indian_pines_gt']\n",
    "    \n",
    "    return data, labels\n",
    "\n",
    "def reports (X_test,y_test):\n",
    "    Y_pred = model.predict(X_test)\n",
    "    y_pred = np.argmax(Y_pred, axis=1)\n",
    "    target_names = ['Alfalfa', 'Corn-notill', 'Corn-mintill', 'Corn'\n",
    "               ,'Grass-pasture', 'Grass-trees', 'Grass-pasture-mowed', \n",
    "                'Hay-windrowed', 'Oats', 'Soybean-notill', 'Soybean-mintill',\n",
    "               'Soybean-clean', 'Wheat', 'Woods', 'Buildings-Grass-Trees-Drives',\n",
    "               'Stone-Steel-Towers']\n",
    "\n",
    "    \n",
    "    classification = classification_report(np.argmax(y_test, axis=1), y_pred, target_names=target_names)\n",
    "    confusion = confusion_matrix(np.argmax(y_test, axis=1), y_pred)\n",
    "    score = model.evaluate(X_test, y_test, batch_size=32)\n",
    "    Test_Loss =  score[0]*100\n",
    "    Test_accuracy = score[1]*100\n",
    "    \n",
    "    return classification, confusion, Test_Loss, Test_accuracy\n",
    "\n",
    "\n",
    "def applyPCA(X, numPCAcomponents=75):\n",
    "    newX = np.reshape(X, (-1, X.shape[2]))\n",
    "    pca = PCA(n_components=numPCAcomponents, whiten=True)\n",
    "    newX = pca.fit_transform(newX)\n",
    "    newX = np.reshape(newX, (X.shape[0],X.shape[1], numPCAcomponents))\n",
    "    return newX, pca\n",
    "\n",
    "def Patch(data,height_index,width_index):\n",
    "    #transpose_array = data.transpose((2,0,1))\n",
    "    #print transpose_array.shape\n",
    "    height_slice = slice(height_index, height_index+PATCH_SIZE)\n",
    "    width_slice = slice(width_index, width_index+PATCH_SIZE)\n",
    "    patch = data[height_slice, width_slice, :]\n",
    "    \n",
    "    return patch"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "X_test = np.load(\"X_testPatches_\" + str(windowSize) + \"PCA\" + str(numPCAcomponents) + \"testRatio\" + str(testRatio) + \".npy\")\n",
    "\n",
    "\n",
    "y_test = np.load(\"y_testPatches_\" + str(windowSize) + \"PCA\" + str(numPCAcomponents) + \"testRatio\" + str(testRatio) + \".npy\")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [],
   "source": [
    "X_test  = np.reshape(X_test, (X_test.shape[0], X_test.shape[3], X_test.shape[1], X_test.shape[2]))\n",
    "y_test = np_utils.to_categorical(y_test)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [],
   "source": [
    "# load the model architecture and weights\n",
    "model = load_model('my_model' + str(windowSize) + 'PCA' + str(numPCAcomponents) + \"testRatio\" + str(testRatio) + '.h5')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2563/2563 [==============================] - 0s 112us/step\n"
     ]
    }
   ],
   "source": [
    "# Using the pretrained model make predictions and print the results into a report\n",
    "classification, confusion, Test_loss, Test_accuracy = reports(X_test,y_test)\n",
    "classification = str(classification)\n",
    "confusion = str(confusion)\n",
    "file_name = 'report' + \"WindowSize\" + str(windowSize) + \"PCA\" + str(numPCAcomponents) + \"testRatio\" + str(testRatio) +\".txt\"\n",
    "with open(file_name, 'w') as x_file:\n",
    "    x_file.write('{} Test loss (%)'.format(Test_loss))\n",
    "    x_file.write('\\n')\n",
    "    x_file.write('{} Test accuracy (%)'.format(Test_accuracy))\n",
    "    x_file.write('\\n')\n",
    "    x_file.write('\\n')\n",
    "    x_file.write('{}'.format(classification))\n",
    "    x_file.write('\\n')\n",
    "    x_file.write('{}'.format(confusion))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Build the Classification Maps"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [],
   "source": [
    "# load the original image\n",
    "X, y = loadIndianPinesData()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [],
   "source": [
    "X,pca = applyPCA(X,numPCAcomponents)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [],
   "source": [
    "height = y.shape[0]\n",
    "width = y.shape[1]\n",
    "PATCH_SIZE = windowSize\n",
    "numPCAcomponents = numPCAcomponents"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [],
   "source": [
    "# calculate the predicted image\n",
    "outputs = np.zeros((height,width))\n",
    "for i in range(height-PATCH_SIZE+1):\n",
    "    for j in range(width-PATCH_SIZE+1):\n",
    "        target = y[int(i+PATCH_SIZE/2), int(j+PATCH_SIZE/2)]\n",
    "        if target == 0 :\n",
    "            continue\n",
    "        else :\n",
    "            image_patch=Patch(X,i,j)\n",
    "            #print (image_patch.shape)\n",
    "            X_test_image = image_patch.reshape(1,image_patch.shape[2],image_patch.shape[0],image_patch.shape[1]).astype('float32')                                   \n",
    "            prediction = (model.predict_classes(X_test_image))                         \n",
    "            outputs[int(i+PATCH_SIZE/2)][int(j+PATCH_SIZE/2)] = prediction+1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAATkAAAEyCAYAAACBJqcyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAWd0lEQVR4nO3df+xddX3H8edrRVBUpB0rq7SuuDCXuszpGlbUGGP9wSyhLBlJyVjqbNMsOIdui7SSjOwPkjqN02Tplm8KUidCGsTRgL861JglK9qiaEvBFsHylUphZrjMBa177497vvT29txf55x7zzmf+3qQcu895/54f3+98zqfz/mhiMDMLFW/UncBZmaT5CZnZklzkzOzpLnJmVnS3OTMLGlucmaWtIk1OUmXS3pU0lFJWyf1OWZmg2gS+8lJWgR8H3gHMA98C7gmIh6u/MPMzAaYVJK7FDgaET+IiJ8DdwLrJ/RZZmZ9nTWh970IeLLr8TzwB/2eLF0QsHJCpZRxAPj9AetsUl7xilfVXcLYnnvu2GmPL/jtmgoZwbOP1F1BpZ6NiF/rt3JSTU45y07bLpa0BdjSefQqYP+ESilD9K8r70u0qrz5ze0bxr3vvutOe3zVbfXUMYqda+quoFI/HLRyUpur88CKrsfLgae6nxARcxGxOiJWQ98mXBPhJmaWhkk1uW8Bl0i6WNLZwAZgz4Q+y8ysr4lsrkbESUl/AXwZWATcGhGHJvFZ1epNbyVmnjfvO3NZYtsIZm0wqTE5IuILwBcm9f5mZqOYWJNrlwoT3CB56Q6c8MwmaMabXEMmF7xp+4J163bUXYIlxseumlnSZjTJ9UtwFW6mFkli3YmuN93NaLIzK8tJzsyS1pAkd4DT09WkLq5TfYKLntcq+4x9bB75Pdaws3NnIa3ljdF50sKsECc5M0taQ5Jcr4XEVVWim8IYXF08M2s2kJOcmSWtoUluQRWJriH7wk1TC8fvvH+cTYqTnJklreFJboEYnubGSWwJjMUV0eB973rPxdbLSW+y1u1YV9l73XfdfZW9VxWc5MwsaS1JclDN2NqMJrgEDEp6Tnk2iJOcmSWtRUmuDCe40ho8Y+vxPBvESc7MkjYjSa7qIygarsp01S/BDVrfgHTXrV/Sc8KbDTPS5GZTbM4/SYB27pzsBzd407abJzNmgzdXzSxpM5bkZmyztY/ehJeb7IZtphbRkoQHnsxIiZOcmSVtxpJc9dS0EwDsXNN3LK6f7udPfLwuT3fCa2CqyzMs6VlzOMmZWdJmNMlN/iSaL5zSvG0mMRZX5PNbkuis+ZzkzCxp7Uly05gQVV7CG+2Dy8y2NWV8Rzt3wpjjedbfoDBad2CeJU5yZpa09iS5aegNbXrhf82XxYZxZ1ZP43gxNf1Snn8E1Suc5CStkPQ1SYclHZJ0fbZ8iaS9ko5kt4urK9fMbDxlktxJ4K8j4kFJLwcOSNoLvAe4PyK2S9oKbAVuKF9qDbqTXUsCXREv7Bvn8bjaDZtUdtIbX+EkFxHHI+LB7P5/A4eBi4D1wK7sabuAq8oWaWZWVCVjcpJWAq8HHgAujIjj0GmEkpb2ec0WYEsVnz8VueN17XYqwTketEVe0qvixzfOxWeqvOjNNJRucpJeBnwO+EBE/FS5u2GcKSLmgLnsPWb7iPkySkw4LLxGC385bnatNO1N3IWG2JZmV2oXEkkvotPgbo+Iu7PFT0talq1fBpwoV6KZWXGFk5w6ke0W4HBEfLxr1R5gI7A9u72nVIVmVsqsHyFXZnP1TcCfAt+T9J1s2YfpNLfdkjYBx4Cry5VoZlZc4SYXEf9O/+H3tUXf18ZTaudfsxLGmayokw/rMrOk+bCutsoGWnJPKdC2dDfrg0Y2UU5yZpY0J7lhmrbT7wipp/cU5q1LdmYVcpIzs6Q5yTVRxWNUC8nOic5mkZOcmSWtPUmuaWNjmaQuMlw2QfrYV2sgJzkzS5oi6j8BiM9CYmYlHIiI1f1Wtmdz1WbaPvpPmgy7xm3va7uf3+99i1w3ty2nHsrTlkO0ivDmqpklLaEkN+0t3p6ZkFnf4G7oxFCeQamwjLadTHJWOMmZWdLc5MwsaW5yZpa0Ro/JjbJ3y6nr5vQOCs36IJlNQvd4XpEZWJs+JzkzS1qjk9woetNefrJzqrPqLaS6mUt0+4bMTq9p1vfDSc7Mktb6JNcrP9l5vM7GM6l96VptWIIb9Lwa052TnJklLbkk16s72Xkm1qygUVNc0ddPMOk5yZlZ0pJPct1Gm4kFJzuzTNkEV/ZzKkh4M9Xkei00PZ1xcLl3P7EZN63mNgXeXDWzpM10khuNN2Wbzrt7VCihBLegdJKTtEjStyXdmz1eImmvpCPZ7eLyZZqZFVNFkrseOAyclz3eCtwfEdslbc0e31DB5zREi84OabVI+VTibVQqyUlaDqyD0w7eWw/syu7vAq4q8xlmZmWUTXKfAD4EvLxr2YURcRwgIo5LWlr0zc+c9Rzt9Etmjdb9O6ycZXmy5+WNPy6cIGCcy96WvcRumxROcpKuAE5ExIGCr98iab+k/UVrMDMbpkySexNwpaR3Ay8GzpP0GeBpScuyFLcMOJH34oiYA+ZgvOuuLqS7xiW6vKG6ptU4SVV/rR76tIoUTnIRsS0ilkfESmAD8NWIuBbYA2zMnrYRuKd0lWZmBU1iP7ntwG5Jm4BjwNUT+Ix28C52s8U/30aqpMlFxNeBr2f3/xNYW8X7mpmV5SMeKDa+d8bMb/+j/7uW9Stg/M83s9H42FUzS5qT3KTkn60zn8fuzCbGTW4aRtmUPW395Eo5jZupzQBvrppZ0pzk6jBuspsUbybbDHCSM7OkOck1QROTnVOdJcJJzsyS5iTXUL2nzZn6qXE8XmeJcJIzs6Q5ybWEk51ZMU5yZpa01ia5uiYgm6JxyW4UTn9WAyc5M0tas5PcOFfmmLaGXQmk+1vVsNJOGfWiLWYVcpIzs6Q1O8lZIbWP15k1iJvcDGhc05vxSSObLm+umlnSnORmUOOSndkEOcmZWdKc5Kwdu5+YFeQkZ2ZJc5Kz03i8zlLjJGdmSXOSs4Gc7KztnOTMLGnNTnKODY3jZGdtUyrJSTpf0l2SHpF0WNJlkpZI2ivpSHa7uKpizczGVTbJfRL4UkT8saSzgXOBDwP3R8R2SVuBrcANxd6+Cefk8YGWg3gfu+ab9Z9L4SQn6TzgLcAtABHx84j4L2A9sCt72i7gqrJFmpkVVSbJvRp4BviUpNcBB4DrgQsj4jhARByXtDTvxZK2AFtKfL41jMfrErBmZ90VVK5MkzsLeAPw/oh4QNIn6WyajiQi5oA5AElN2C4daDOn/wXvxH/Bw7jpWROUmXiYB+Yj4oHs8V10mt7TkpYBZLcnypVoZlZc4SQXET+W9KSk10TEo8Ba4OHs30Zge3Z7TyWVWus52ZW3hv6bk/5+5is7u/p+4PZsZvUHwJ/RSYe7JW0CjgFXl/wMM7PCFFH/cFj/Mbn6a1vYhWTomFyV30ep2i892wsmxYuf7WNztYUUMChdjWzYz9t7Mg1yICJW91vpw7rMLGnNPqzLZkbR8bruFNWEVAenQpnDVzM4yZlZ0pzkLBm9Y2N1J7uxEp1j38Q4yZlZ0pzkrBX6zQwPGrtrSrLzGF29nOTMLGlOctZqeQmvX7qrO9lNK9HFgJ3uNIN50knOzJLmJGczq6597AKPz02Tm5wZ09+U9WTE9Hhz1cyS5iRnlqPuSQqrjpOcmSWt4UmuOSMWUz/d+QS+9CpPqtjk0zZNgpNdeznJmVnSGp7kWkTNSZ0patqpvfNOlDmNU7zuGxAgE7zQViWc5MwsaY1KcoMOR2mKcQ6Lacq4VdNSUCFN+dUY8OPvt6oppc8qJzkzS1qjklxqFhJUUxKdVSAvlg0J972rneymy0nOzJLmJGdjcSotz/Pw0+UkZ2ZJc5Mzs6S5yZlZ0tzkzCxpnniYIZ40sFlUKslJ+qCkQ5IOSrpD0oslLZG0V9KR7HZxVcWamY2rcJKTdBHwl8CqiPhfSbuBDcAq4P6I2C5pK7AVuKGSas1mnA/CH1/ZMbmzgJdIOgs4F3gKWA/sytbvAq4q+RlmZoUVbnIR8SPgY8Ax4DjwXER8BbgwIo5nzzkOLK2iUDMbLLL/7HSFm1w21rYeuBh4JfBSSdeO8fotkvZL2l+0BjOzYcpsrr4deDwinomIXwB3A28Enpa0DCC7PZH34oiYi4jVEbG6RA1mllGf/2ZdmSZ3DFgj6VxJAtYCh4E9wMbsORuBe8qVaGZWXOHZ1Yh4QNJdwIPASeDbwBzwMmC3pE10GuHVVRRqVpkZCzeznuZK7QwcETcBN/Usfp5OqjMzq50P65qCnWsSOQW5WQu5yZlZ0nzs6ohmfVzDrK3c5IZwczNrN2+umlnSnORyOL2ZpcNJzsyS5iTXxQnOLD1OcmaWtEYlOScpM6uak5yZJc1NzsyS5iZnZklzkzOzpLnJmVnS3OTMLGlucmaWNDc5M0tao3YGNuvL+4lbQU5yZpY0J7mSNu+ruwIzG8RJzsyS5iRn7RBR3XvJA3yzxEnOzJLmJGcj8XVjra2c5MwsaW5yZpY0NzkzS5rH5Gw8Vc5yDuIZUKvI0CQn6VZJJyQd7Fq2RNJeSUey28Vd67ZJOirpUUnvmlThZmajGGVz9Tbg8p5lW4H7I+IS4P7sMZJWARuA12av2SFpUWXVmlUpopp/1mhDm1xEfAP4Sc/i9cCu7P4u4Kqu5XdGxPMR8ThwFLi0olrNzMZWdOLhwog4DpDdLs2WXwQ82fW8+WzZGSRtkbRf0v6CNZiZDVX1xEPeaHFuno+IOWAOQJIzv03ewmSGNzFnStEk97SkZQDZ7Yls+Tywout5y4GnipdnZlZO0Sa3B9iY3d8I3NO1fIOkcyRdDFwCfLNciWYVk079s+QN3VyVdAfwVuACSfPATcB2YLekTcAx4GqAiDgkaTfwMHASeF9E/HJCtZuZDTW0yUXENX1Wre3z/JuBm8sUZTY1vWnO43XJ8WFdZpY0Nzmzbh6rS46bnJklzU3OLI8TXTLc5MwsaT7VkqXPiWymOcmZWdKc5Kz9nNRsACc5M0uam5yZJc1NztrNm6o2hJucmSXNEw/WTk5wNiInOTNLmpucmSXNTc7MkuYxOWsXj8XZmJzkzCxpTnLWDnUlOCfH1nOSM7OkucmZWdLc5MwsaYoGXIJNUv1FmFlbHYiI1f1WOsmZWdI8uzpDguKBWXiW0drJTa6gHevW1V3CyK677766S2isdTva83OcpvuuS+d3xpurZpY0NzkzO8O6HeuSSblDm5ykWyWdkHSwa9lHJT0i6buSPi/p/K512yQdlfSopHdNqnAzs1GMkuRuAy7vWbYX+J2I+F3g+8A2AEmrgA3Aa7PX7JC0qLJqG2DHunWtGo8zKyOFRDe0yUXEN4Cf9Cz7SkSczB7uA5Zn99cDd0bE8xHxOHAUuLTCes3MxlLFmNx7gS9m9y8CnuxaN58taz0nOJtlbU5zpXYhkXQjcBK4fWFRztNyd86StAXYUubzzcyGKdzkJG0ErgDWxqljw+aBFV1PWw48lff6iJgD5rL3auxhXU5vZh0Laa5t+9AV2lyVdDlwA3BlRPysa9UeYIOkcyRdDFwCfLN8mWZmxQxNcpLuAN4KXCBpHriJzmzqOcBedU4quC8i/jwiDknaDTxMZzP2fRHxy0kVP2lOcWZnaluiG9rkIuKanMW3DHj+zcDNZYoyM6uKj13N4QRnlg4f1mVmSXOS6+IEZ5YeN7kuPiXRbGnzDq42Om+umlnSnOSAPgdlJMRn9bXZ5SRnZklzkpshvk6DzSInOTNLmpucmSXNTc7MkuYmZ2ZJc5Mzs6S5yZlZ0tzkzCxpbnJmljQ3OTNLmo94sJnjs48U05bTnfdykzOzvtra2Lp5c9XMkuYmZ2ZJc5Mzs6R5TA7wSSXN0uUkZ2ZJc5Mzs6S5yZlZ0jwmNyH72Fx3CS9Yw866SzCrjZOcmSXNTc7Mkja0yUm6VdIJSQdz1v2NpJB0QdeybZKOSnpU0ruqLtjMbByjjMndBvwj8OnuhZJWAO8AjnUtWwVsAF4LvBL4N0m/FRG/rKpgs7JSOB7TRjc0yUXEN4Cf5Kz6B+BDnH75+fXAnRHxfEQ8DhwFLq2iUDOzIgqNyUm6EvhRRDzUs+oi4Mmux/PZsrz32CJpv6T9RWowMxvF2LuQSDoXuBF4Z97qnGWRs4yImAPmsvfMfY6ZWVlF9pP7TeBi4CFJAMuBByVdSie5reh67nLgqbJFmpkVNXaTi4jvAUsXHkt6AlgdEc9K2gN8VtLH6Uw8XAJ8s6Jak+OddM0mb5RdSO4A/gN4jaR5SZv6PTciDgG7gYeBLwHv88yqmdVpaJKLiGuGrF/Z8/hm4OZyZZmZVcNHPJhZ0hRR/8SmZ1fNrIQDEbG630onOTNLWlNOtfQs8D/ZbdNcQPPqamJN0My6mlgTNLOuJtYEw+v6jUEvbsTmKoCk/YMiZ12aWFcTa4Jm1tXEmqCZdTWxJihflzdXzSxpbnJmlrQmNbm5ugvoo4l1NbEmaGZdTawJmllXE2uCknU1ZkzOzGwSmpTkzMwq14gmJ+ny7HTpRyVtramGFZK+JumwpEOSrs+WL5G0V9KR7HZxDbUtkvRtSfc2qKbzJd0l6ZHse3ZZ3XVJ+mD2szso6Q5JL66jprxLBgyqY1qXDOhT10ezn+F3JX1e0vnTrGsql1eIiFr/AYuAx4BXA2cDDwGraqhjGfCG7P7Lge8Dq4C/B7Zmy7cCH6mhtr8CPgvcmz1uQk27gM3Z/bOB8+usi87JWR8HXpI93g28p46agLcAbwAOdi3LrSP7HXsIOIfOKcweAxZNsa53Amdl9z8y7bryasqWrwC+DPwQuKBMTVP9w+jzRV4GfLnr8TZgWwPquofONSweBZZly5YBj065juXA/cDbuppc3TWdlzUU9SyvrS5OnZV6CZ2d3O/N/oBrqQlY2dNMcuvo/X3P/rAvm1ZdPev+CLh92nXl1QTcBbwOeKKryRWqqQmbqyOfMn1aJK0EXg88AFwYEccBstul/V85EZ+gcy2N/+taVndNrwaeAT6VbUbvlPTSOuuKiB8BH6NzYaXjwHMR8ZU6a+rRr44m/f6/F/hidr+2uqq4vEK3JjS5kU+ZPg2SXgZ8DvhARPy0rjqyWq4ATkTEgTrryHEWnU2Mf4qI19M5JK+WsdQF2RjXejqbMa8EXirp2jprGlEjfv8l3QicBG5fWJTztInX1XV5hb/NW52zbGhNTWhyjTlluqQX0Wlwt0fE3dnipyUty9YvA05MsaQ3AVdmZ1++E3ibpM/UXBN0fmbzEfFA9vguOk2vzrreDjweEc9ExC+Au4E31lxTt3511P77L2kjcAXwJ5FtB9ZYV/flFZ7g1OUVfr1oTU1oct8CLpF0saSz6Vy3dc+0i1DnghW3AIcj4uNdq/YAG7P7G+mM1U1FRGyLiOXROTHpBuCrEXFtnTVldf0YeFLSa7JFa+mcDbrOuo4BaySdm/0s1wKHa66pW7869gAbJJ0j6WKmfMkASZcDNwBXRsTPeuqdel0R8b2IWBoRK7Pf+3k6E4I/LlzTpAY4xxx4fDed2czHgBtrquHNdKLvd4HvZP/eDfwqnYH/I9ntkprqeyunJh5qrwn4PWB/9v36V2Bx3XUBfwc8AhwE/oXOLNzUawLuoDMu+Ivsj3TToDrobJ49Rmdy4g+nXNdROuNcC7/z/zzNuvJq6ln/BNnEQ9GafMSDmSWtCZurZmYT4yZnZklzkzOzpLnJmVnS3OTMLGlucmaWNDc5M0uam5yZJe3/ASVWZpuJNQ9wAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Plot the Ground Truth Image\n",
    "ground_truth = spectral.imshow(classes = y,figsize =(5,5))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAATkAAAEyCAYAAACBJqcyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3dfaxk9X3f8fe3EOzgBwGl0DWLCo6oIxw1NaUwTirLCklMc1deKgUJK442DbOret3UTVuFJRa1qgiVNJGb/LOudsdOti0BrQgJaG8Sm5C4VqXO2oBjhwdj1sGBNWvWrps+JXKC++0fc87cc8+cM3PmPJ/ffF7Sau6cMzPnt3Pn/uZzfk/H3B0RkVD9ta4LICLSJFVyIhI0VXIiEjRVciISNFVyIhI0VXIiErTGKjkzu9XMnjezM2Z2pKnjiIgsY02MkzOzC4AvAz8CnAU+B7zP3Z+t/WAiIks0leRuAs64+5+4+18CDwL7GzqWiEiuCxt63auAlxP3zwI35z3YzDTtQkTK+qa7/428nU1VcpaxbVdFZmaHgEMNHV9ENsefLtvZVCV3Frg6cX8v8EryAe5+DDgGSnIi0pym2uQ+B1xnZtea2UXAHcCjDR1LRCRXI0nO3V8zs38KfBK4APiEuz+zxissbJlyEIDReLx7xyRq6vPZGfLUZvtHHJ8/ZDyd7ZuMVgXG9Fn2ssdnnZEnjKeL2yajFceX2NbW0QpP3s7evr1V/jUL2N4+vOt+/BGY5LZGVzM+vXz/ZJLxGZzv3JzPYlOnq7j77wC/09Tri4gU0VglV7d5MpvMbubJLu60jZr14u+nnf3Ha05wFWSlO9iob9VWNZzc+ipOcOPxaNf9TdXIYOC1C7HQ8bDkdDVxGlrDkTO3Tlk85S363Lm8Cq2oDa34Kp2mdix9utrUd2VMp6tzT7r7jXk7NXdVRII2vNPVlOKdCkl5KWz3KW8lZb4pk+kvnQTD+uaVBqw6LR2PRxt56qokJyJBG0ySS4sTXPr+8kS3PMGVEj3Vo2ZFi44xZcxkOuslGY/GmU+dD4eJU1p8m9GeN28nXBhCo4TXWwdTv8fj3fyu4g6ITaUkJyJB62mSixNXfsKKE9s8wS0dcNlAglthMp3kJrjYdDJLekW+Z9MJLn4uccLz6P7BqdJdT+QO4Tg4XT/Vxakw43mr2tniv5UpB5ks/0gGSUlORILW03FyaVXKuGw8W/XX9dRrJNvkihrFI5xjS8bYzdNfum2uyGN6nPCGPD4u1sQ4uayxcPEZQtzmW/i1RuPFduAwaJyciGyugSQ5yPtaHDP7qpusNbqtjv9zlBCj9y9aH6DxJLeOnXa76DgDSnZpQ0h6jUzQz+ihLTvTIfl5SH/mto7mTIFLTo3LW/gg/ZTDxR5XIyU5EdlcPe1dzZLdtrZeq0T9qdVXTGFt23QyyW2vW0h2bRSoJgvtXQl9TXlxgpsv/+XrfVqBeW9qnN4mB6dQctxbkR7/BXF6297aSXUFE11fKMmJSNAGlOSq6L7dsS1ZKS5/jF1xCzMujsdLXHUfZZelPGgo6WWlmdSmOH3FCa5Kopu36908gpw2uby2uHi83topLmlg6S1JSU5EgrYhSW71DIqgRL2m85abKm1x8y7CKI3ErxXd/kKUkk6dOsy+fbOfVyWrtuWVp+m2vDh9zdvTEuPa5tsm+TMZko+T8jakkmvOvq2PAbM/cmBen47G41KnhXXynA4IK1GuvFPe7T03zF7z4BT2RDsHsgJypc6M7a3djfKzH3Y/Jhr+McmaihWffq74VWQNP1m34qt0mhoAna6KSNA2LMnVf9oaJ5khSSe8zGSXk8YqpdOBJDwo0JlxdKv4NSQaXnKpNwth5iXajinJiUjQBjStq04VDjdvMV7xbZxILavST63Tuiaj3La4IuaprqapZWvrYaorLf0e1pDglrXHFZvOFV9GYPfQn9xpXeuIklwHHU+a1iUim2vD2uRiFRbRLDrbOpFIhpRN4hTY1RDfuFezb8NQapFumyNjWawKaS9vkc7kMvw7l/bcHEpyIhK04bTJlW0jWmudm6z8kle0VJtGYlAsMB8YW8RCaumoTc4mk/LHPj5aXKI7I7kUFlLbXJZV73MNg4PTiW46mTTbJhfRUksiIi0aTJvcygUg84xPr5Hm0qnNWLd1ap0EV6so+VTpWfXxuNa2OI/SiK1KdIkUuNOuVGNB+ignqc4XW8254sxoPC58icH04yZjYJR9kfaQlU5yZna1mf2hmT1nZs+Y2Yei7ZeZ2WNm9kJ0e2l9xRURWU+VJPca8C/d/SkzexPwpJk9BvwU8Li732dmR4AjwF3Vi1rOLAGmJ6jnfJulGzwmnlzXPFi1jI3LSGvztyzdvpR+bJW2u8AsjJlc0NI81PRsjk1casndz7n7U9HP/xt4DrgK2A+ciB52AritaiFFRMqqpU3OzK4B3sFsOb8r3f0czCpCM7si5zmHgEOlD7pqPFEiHcSLFM6DWrpdImcN8+lkDLb7m3O+4GEAyS5OcF2PjQPmv8/Qm+Iqm4wW2yvHu9szy1jsEd3OvVtnT2wbKg8hMbM3Av8FuNfdHzazP3P3SxL7/4e7L22XKzKEJO/qV3mrrE4tP9bvXHsy6pBIVXLJs9a8ZWoWrnZUYW2ySkNIauhwyDxdzfsSWXZqmR46ss7A1vRzAvgS6Y2sz1MNQ3TmlV3q2g9BDSExs+8CfhO4390fjja/amZ7ov17gPNVjiEiUkXp01UzM+DjwHPu/tHErkeBA8B90e0jlUpYUnpd/aSdqTS7t6f7HZIpLjk1ZiOU6QzI61AokuhqXn5IEpoaWJ3unOjpUktV2uR+EPhJ4I/N7I+ibT/PrHI7aWZ3Ai8Bt1croohIeaUrOXf/r+S3nNxS9nXrsrRNbp7yZvfjcZdxSotTW/I6lfHtTnteP5rIq7TFde5ggTZA6a2hLKKgaV0iErTBTOvKsyyxwe4LysQJLf2cdIJLbu+qLW7lNLaonSVzSYG+pbtV7XuTkcaOSGOU5EQkaINZailvnNw60mloYQJz6jqZULxjqrVxciV6ylYlu2XTuhYWdayTxsRJPbTUkohsrsG1ycWT6+NlnAs9JyeFLCwqOJosHGPSRWNRzeOa0tO31tFIgovN2+rUsyrNUZITkaANJsntzBWN09YaCo5p23nN4seo0hbXO1UT5JpLNXn0PqtJTpqkJCciQRtM76qISI6lvauDOV2VzRY3C9yz/dTCvlWr6Y6ZnUZPiK/1kThBjr/ko01bR/cB5ZYLmjddpFbRPXX41Mrn7ouO25UOlkdqjU5XRSRoASW5nDPe5IKYOWfFnlo0c1lDePwKln7QeI1rjq5akDLrGqbrWHX2n7MSciUN9x5kJbiixlHSG887lHben6lFV5SPhspsP3VPtGf9ZBMP6t7ayl45NyutFUl5Uo2SnIgELaAkt1ucznYFjHSCiRLPvN8jtd8Tz1/ZM5JOXGUSWPI5VZYeml9hbFbqzPdCdg0ob3T6WiQrtXXdFrcJlOREJGi9TnJFRrfstI2l2tWy2sRWpKOsw609tqVKW1rNlOCyrTMlcPG5icVTCw6eTqc1tcO1S0lORILW6yRXRDrtzZNd8ls27vkseHGWwu1wWXqQ4GJKcOubt80VelyU6tZcxGHw7XDTFe2Wo36tgKokJyJBG0ySK5quMpNduu0kTnbx6pjxapl16FObXHSrRLe+OhZpDc6qBLfscR2mOyU5EQnaYJJc2Rn8yWSX2V6XePX0OgFWZmZADxJcLKQEF49fS1/gZ63XKLHgahlx7+ng296Siqa4ss9vMOkpyYlI0AaR5JKJpMqaTLk9senMEz/u4HQhmamdS1ZRgqvxODUkvEFUck0tNjdfZSddY83vj2CSOoWdX8qrP6elUkzTp6lBaatya4FOV0UkaINIcma+sBxSi0ffdc+joQVtl6aJ0S6h0HCPCrbjZaGipaUCSnCxyknOzC4ws8+b2ano/mVm9piZvRDdXlq9mCIi5dSR5D4EPAe8Obp/BHjc3e8zsyPR/buqHKC7FLfI8oYwpKeM1TyURAluOIa1lPiQylpOpSRnZnuBLdg1eW8/cCL6+QRwW5VjiIhUUTXJ/Qrwc8CbEtuudPdzAO5+zsyuKPviC72eFFt+ScJTZRBw32wd3eLUB2YJKv6MxzMN48S+tW/3RXHidLh1dHFp9fS+7Q8USGf9OTlqXOkkZ2b7gPPu/mTJ5x8ysyfM7ImyZRARWaX0dVfN7N8CPwm8BryeWZvcw8DfB94dpbg9wKfd/W0rXmvtQtSZ6LISYy02NHXWMmA650JBfUx06y61tOtzUXTliVTiS5rPUlzn85b3y+lb72qxwcBLr7taOsm5+93uvtfdrwHuAP7A3d8PPAociB52AHik7DFERKpqYpzcfcBJM7sTeAm4vYFjDEP623JDkl0jwTjurZ7MksZkuvgNPx51m0Li8Xpx2druEddYymy1VHLu/mng09HP/x24pY7XFRGpahAzHppWpn0v3Y63tTXrDds+9YHsB0B+xImOr2/iauIE1VWii5dyYhTdVhjfueyzMLk5/YHdoK7SEjR3VUSCpiRXk3mCi2Wu1pkj2j1vZRovLvE0f9nsp679GGlQDTN0lqX58enZ68eJNf1J0RnBbqrk1lTq45u/kF22yYi8UQlFjl+4jAPsCMnqcNhU8/dinRmE8+mHm7NUmE5XRSRoSnJrqmWg67rJrimBDXHpegiJ9JOSnIgETUmuD/qY7AaY6roeQlJHh0PavBOhynOTHRDzQdVVSjUsSnIiEjQlud6aRanxdJYOFi4V27SetdfNrx+07DE9aZOrcwiHhoFUpyQnIkEbbJKzaLyPr7HMeNFVbXqp6xG+6UvTeubmxsSJZmr9SGvLKH31i5KciARtsEkubqSyEr1EQ5jmFLfFze+n2qQmHbfZleoAHmSELq7OtrjF3/cazx3lL0e1iZTkRCRo/U5yWWs990Xr3Z27pZPerveq47LlWtIoumomyRDa4uqQToNFepXnz13Ruzw+vZnthUpyIhK0fie5AuILm4zGYX/Tr/yWjlptJtMJk0nXXbErxCthrNEz3vlshjVUaZsrk+AWXiOnLW5yc7UZFEM12EoufdWmUCu79B91ocbkeHrR/CJoPav0JotTi3pWDQ9S3mcjWenqdFVEJDCDTXKhJba0+Tr+6Q6GAuanOgWf29d+imWGdPoq3VKSE5Gg9TvJ1dEKKyvF1wuFEleDb0k6sfUxwQ1hEO4mdjwoyYlI0Pqd5GLJLqFNSXXxFbxyrsla59uQTB5xPtq5EtQmfefXY57ovPx710Qv6Kb86aQpyYlI0IaR5OpWMA6Nx4vdjpNJO1PN5tO2UkVc9W2cHPC58hipNqRdg0XjdJdOFBrQJgOjJCciQet5ksuIDXU2Ea14rUnm/n4OKis1hSjVC5iVABe2pYLsEMfYbbJzT21FP20XeryPdj4jNu1fj3YRlZKcmV1iZg+Z2ZfM7Dkze6eZXWZmj5nZC9HtpXUVVkRkXVWT3K8Cv+fuP25mFwEXAz8PPO7u95nZEeAIcFe5l+9mlcUpB+c/9713MSvBtTm8MF7haTwa9/692lS703axBBcbanpLKp3kzOzNwLuAjwO4+1+6+58B+4ET0cNOALdVLaSISFlVktxbgW8Av2Zm3w88CXwIuNLdzwG4+zkzuyLryWZ2CDhU4fiNGXE8ca/f6aTOJbermEwnGmMX2fn8dPP/r5TkR+H9zqpUchcCNwA/4+6nzexXmZ2aFuLux4BjAGbW+9X/x6kW90lPOyCg+0Gf88nz0f1Nr/SKyPuyqjINa91OhlBV6Xg4C5x19/hP6iFmld6rZrYHILo9X62IIiLllU5y7v51M3vZzN7m7s8DtwDPRv8OAPdFt4/UUlIppOsUlyVOdtP04OP+huFSdjdz1COd7Ja9Z+l92xue4GJVe1d/Brg/6ln9E+AfM0uHJ83sTuAl4PaKxxARKc3cu28Oy2+T675s8YDklW1ydb6PZpC6rmpRWdO64m/4EC9+llwmqit1tDWmfzcLHUmaTrfMk+5+Y95OTesSkaD1fFqXbIqFJFMw2SVTVB9SHay+hqy0S0lORIKmJLeM22C+jrseDNwH6baxrpPdOokutJ7mPlGSE5GgKcll8YHEt4Q+jo+rU17P8LIE1HWyiye3x8sVDe9TFQYlOREJmpJcll3D9vT922fJhLdqjmxWsotnKSSX16qLpya7t9Xr6kvGl9oGfp6V5EQkaEpyPbXuTIdNVvZC07NkN4l+nqWfJhJdmtPO+cG+o/sAOHX4VLRlMz9TquRk8OZLO5Ws7GZmFUDcj9F0J0Ubp647ldtm0+mqiAQt2CQXn3Y0sfxNH2UNBl6cqF9u0n8fJf+/C9eKrWFgbdfDT5rR1olyvyjJiUjQer7U0oDUvtRSvzW9bNPCIN913t4W3r5kskunvq6LuoFDSLTUkohsrmDb5Fo3gPQ1aD17e7MGHLdxOrLsMqgBXmirFkpyIhK0Xia5ZW0KXVurTSNup4ueku7VjHs7Y031egaxjE8TH4llTcHxIg0Lj8n/HeXt6e+neTMoyYlI0HqZ5IIxb6db77s8mfBWpb+8x0kBpZbUynr/l79Oeq+SXbuU5EQkaEpyLchLWfH2ZW1z49QFmeOZHOn5mm11rPX5soaNyEt75ol2u/VeUpm7XUpyIhI0JbkWpNPYwv5U0staTaPaChuyrpV9qgNcIn9TqZJr0crKbo2KLP3YVa/dd30rv6qwcOh0VUSCpiTXga5OPdOdBjvXRDieeEw3GaYvCS6tresySHMqJTkz+1kze8bMnjazB8zs9WZ2mZk9ZmYvRLeX1lVYEZF1lU5yZnYV8M+A6939L8zsJHAHcD3wuLvfZ2ZHgCPAXbWUVuaqJJ94yEr6FZqbVtbO9LUm9K2kmoS/vqptchcC321mFwIXA68A+4ET0f4TwG0VjyEiUlrpJOfuXzOzXwZeAv4C+JS7f8rMrnT3c9FjzpnZFTWVdaOlk1uVpczTg5CbTlZDSm5Dtmphi31bH4PtlgrTI6WTXNTWth+4FngL8AYze/8azz9kZk+Y2RNlyyAiskqV3tUfBl50928AmNnDwA8Ar5rZnijF7QHOZz3Z3Y8Bx6Lnas5yyqo2tzrSUWc9qYFcUKdvPa95y4BtbR1tuST9UqVN7iVgZGYXm5kBtwDPAY8CB6LHHAAeqVZEEZHyqrTJnTazh4CngNeAzzNLZm8ETprZncwqwtvrKOgm6OtYsbp1nuBqOnxfEtwq29uHuy5CpyoNBnb3jwAfSW3+NrNUJyLSOc14aEF80ZP4MnabkthE+kBzV0UkaEpyBdVxUd6FhS6V6JoVd9prWaSNpkpuhSauOK7KrWVrXHFLwqPTVREJmpJchibSWx+FMig3V+71GdothnRLSU5EgqYkl7ApCS4WbIJLUwfERlOSE5Gg9TLJhZqobviFLQCeumcD17sR6YiSnIgErZdJLlSH79kDQFdXUA2+NzWP2uI2mpKciARNlVyL/OYJfnP3sx0mI991cZn0/dbKcfPsn0iTVMmJSNDMvfuVxzdm+fP4f3lwdpXn8XiU+bDJzTA+Xf1w85Q0b5Lqpk0uLsc6/6dJ+q1Z4xMSH2eeTtNtcmqiC82T7n5j3k4lOREJmnpXWzS/8MnxWUw5Ht0/GCePRPtUmfSTfo2+qCOV5jLPT6ZRghtnvL+yOVTJtSj+U4z/6KzAkI6ylV1wct8iWzy1Tem+q0e6pNNVEQmaklxF8TUt99xT4IpIqTQ2v5J9jadRWYlPSUY2mZKciARNSa6i7T03ADtTtcpcv6GthvF029+gpnnVOdTJBvD/ldooyYlI0JTkKppOZoltEkW5KhepaasHtcwUrlU9mCJ9pSQnIkFTkttAg2iDE6mJkpyIBE1JrqKNu1B0Wws6qAdUarIyyZnZJ8zsvJk9ndh2mZk9ZmYvRLeXJvbdbWZnzOx5M3tPUwUXESmiyOnqrwO3prYdAR539+uAx6P7mNn1wB3A26PnHDWzC2orbQ+NR+P52DgZGPd6/kmvrazk3P0zwLdSm/cDJ6KfTwC3JbY/6O7fdvcXgTPATTWVVURkbWXb5K5093MA7n7OzK6Itl8FTBOPOxttW2Bmh4BDJY/fGxvXJicyMHV3PGS1FmfmeXc/BhyDYa8MPKipUWsoMz2t9+LODJ1ibpSyQ0heNbM9ANHt+Wj7WeDqxOP2Aq+UL56ISDVlk9yjwAHgvuj2kcT23zCzjwJvAa4DPlu1kNK+oBJcWnJ4ilJd8FZWcmb2APBu4HIzOwt8hFnldtLM7gReAm4HcPdnzOwk8CzwGvBBd/9OQ2UXEVlpZSXn7u/L2XVLzuPvBe6tUqhBiWbVB9mGtQnSg46V7IKjaV0iEjRdd7Wi8XT1Y7q2s0xS+Z7g+Wt0Na2rjuOuM1VsneNpClrXdN1VEdlcmqC/QUIby9cojakLhpKciARNSW6DhDo7I8/Wvo8BsL1d4HKREiwlOREJmpKcDF9O7+Z2y8WQflKSE5GgKckFrOmLVYsMgSq5AE0m0Qjl+QyzgC+aqoG4soJOV0UkaEpyIToeJbeDA5hzVlZLCW4+DOXUB1o5ntRPSU5EgqYJ+hVN2X2lrjJLLTW9TFN6gv78uGsMCu7LBP22B/hubR3ddT8z0aldsGuaoC8im0ttchWNxrMUNp2UT2FtL7Q55eDsuK0etRpN0ZKylOREJGhqk6sqWjUzK8n1ZSn0IBbNFMmnNjkR2Vxqk2vQur2myx6f3lfkteePSbW+9fmiO/O2t47LIeFQkhORoKlNTkSGTm1yIrK51Ca3QZzygdlQb6cMkyq5ko5ubXVdhMIOb6sZP8/W0eH8Htu0fTicz4xOV0UkaKrkRGTB1tGtYFLuykrOzD5hZufN7OnEtl8ysy+Z2RfN7LfM7JLEvrvN7IyZPW9m72mq4CIiRRRJcr8O3Jra9hjwfe7+d4AvA3cDmNn1wB3A26PnHDWzC2orbQ8c3doaVHucSBUhJLqVlZy7fwb4Vmrbp9z9tejuFNgb/bwfeNDdv+3uLwJngJtqLK+IyFrqaJP7aeB3o5+vAl5O7DsbbRs8JTjZZENOc5WGkJjZh4HXgPvjTRkPyxycZWaHgENVji8iskrpSs7MDgD7gFt8Z27YWeDqxMP2Aq9kPd/djwHHotfq7bQupTepy6nDpwDYd3RfxyUpJ05zQxtDV+p01cxuBe4C3uvuf57Y9Shwh5m9zsyuBa4DPlu9mCIi5axMcmb2APBu4HIzOwt8hFlv6uuAx2y2uOHU3f+Juz9jZieBZ5mdxn7Q3b/TVOGbphQndRpqgksbWqJbWcm5+/syNn98yePvBe6tUigRkbpo7moGJTiRcGhal4gETUkuQQlOJDyq5BK0JNFmGfIAVylOp6siEjQlOSBnUkZAtKqvbC4lOREJmpJchvgK85NR/xPeOmXVdRpkEynJiUjQlOQGakhpU6RLSnIiEjQluQx9Tkdxgsu73+eyi3RBSU5EgqYkl2FZe1fRtrBkwiqTrtIJbRUlOJFsSnIiEjQluTXVkcqKvEb8mFWJTglOZDklOREJmpJcRcuSVh0pq2iik+K0+kg5Q1nuPE2VXIZ1Kqcqp6/pCqzMaayGkEiThlqxJel0VUSCpiQHOystpc4Ix1NrJBmVec30c5TYRIpRkhORoCnJAVh2o/5kBG0sODk7TjvHEtk0SnIiEjRVciISNFVyIhI0tck1ZMq46yLMjZh0XQSRzijJiUjQVMkFYjKdMJlOFn4W2XQrKzkz+4SZnTezpzP2/SszczO7PLHtbjM7Y2bPm9l76i6wiMg6zH35yHkzexfwf4D/6O7fl9h+NTABvhf4e+7+TTO7HngAuAl4C/D7wN929++sOEZww/fHUxiPlrfLxWkr/bgyKWzZsdQmJ4F70t1vzNu5Msm5+2eAb2Xs+vfAz7H78vP7gQfd/dvu/iJwhlmFJyLSiVK9q2b2XuBr7v4F2z1b4Cpgmrh/NtqW9RqHgENljj8EyWSVl9jS+8u8vtreRJZbu5Izs4uBDwM/mrU7Y1vmqai7HwOORa8Z3OlqsvJZdTq67FQzrxJT5SZSTJkk9z3AtUCc4vYCT5nZTcyS29WJx+4FXqlaSBGRslZ2PACY2TXAqWTHQ2LfV4Ebo46HtwO/wU7Hw+PAdZvY8QBe4PoMLRVFJGzVOh7M7AHgvwFvM7OzZnZn3mPd/RngJPAs8HvAB1dVcCIiTSqU5BovRJBJTkRaUi3JiYgMmSo5EQmaKjkRCVpfllr6JvB/o9u+uZz+lauPZYJ+lquPZYJ+lquPZYLV5fpby57ci44HADN7YlnjYVf6WK4+lgn6Wa4+lgn6Wa4+lgmql0unqyISNFVyIhK0PlVyx7ouQI4+lquPZYJ+lquPZYJ+lquPZYKK5epNm5yISBP6lORERGrXi0rOzG6Nlks/Y2ZHOirD1Wb2h2b2nJk9Y2YfirZfZmaPmdkL0e2lHZTtAjP7vJmd6lGZLjGzh8zsS9F79s6uy2VmPxv97p42swfM7PVdlCnrkgHLytHWJQNyyvVL0e/wi2b2W2Z2SZvlauXyCu7e6T/gAuArwFuBi4AvANd3UI49wA3Rz28CvgxcD/w74Ei0/Qjwix2U7V8wW93lVHS/D2U6AYyjny8CLumyXMwWZ30R+O7o/kngp7ooE/Au4Abg6cS2zHJEn7EvAK9jtoTZV4ALWizXjwIXRj//YtvlyipTtP1q4JPAnwKXVylTq38YOf/JdwKfTNy/G7i7B+V6BPgR4HlgT7RtD/B8y+XYy2zJqh9KVHJdl+nNUYViqe2dlSuq5F4GLmM2yP1U9AfcSZmAa1KVSWY50p/36A/7nW2VK7XvHwH3t12urDIBDwHfD3w1UcmVKlMfTlfjD2csd8n0tkTr570DOA1c6e7nAKLbK1ouzq8wu5bG/0ts67pMbwW+AfxadBo9MbM3dFkud/8a8MvAS8A54H+6+6e6LFNKXjn69Pn/aeB3o587K1fy8gqpXaXK1IdKrvCS6W0wszcCvwn8c3f/X12VIyrLPuC8uz/ZZTkyXMjsFONj7v4OZti3dcIAAAHiSURBVFPyOmlLjUVtXPuZnca8BXiDmb2/yzIV1IvPv5l9GHgNuD/elPGwxsuVuLzCv87anbFtZZn6UMn1Zsl0M/suZhXc/e7+cLT5VTPbE+3fA5xvsUg/CLw3Wn35QeCHzOw/d1wmmP3Ozrr76ej+Q8wqvS7L9cPAi+7+DXf/K+Bh4Ac6LlNSXjk6//yb2QFgH/ATHp0Hdliu5OUVvsrO5RX+Ztky9aGS+xxwnZlda2YXAXcAj7ZdCJtdsOLjwHPu/tHErkeBA9HPB5i11bXC3e92973ufg2z9+UP3P39XZYpKtfXgZfN7G3RpluYrQbdZbleAkZmdnH0u7wFeK7jMiXlleNR4A4ze52ZXQtcB3y2rUKZ2a3AXcB73f3PU+VtvVzu/sfufoW7XxN97s8y6xD8eukyNdXAuWbD448x6838CvDhjsrwD5hF3y8CfxT9+zHgrzNr+H8hur2so/K9m52Oh87LBPxd4Ino/fpt4NKuywX8G+BLwNPAf2LWC9d6mZhdYP0c8FfRH+mdy8rB7PTsK8w6J/5hy+U6w6ydK/7M/4c2y5VVptT+rxJ1PJQtk2Y8iEjQ+nC6KiLSGFVyIhI0VXIiEjRVciISNFVyIhI0VXIiEjRVciISNFVyIhK0/w8ZJq89/EIwFAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Plot the Predicted image\n",
    "predict_image = spectral.imshow(classes = outputs.astype(int),figsize =(5,5))"
   ]
  }
 ],
 "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.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}