Course 2 - Part 2 - Lesson 4 - Notebook.ipynb

{"nbformat":4,"nbformat_minor":0,"metadata":{"colab":{"name":"Horse-or-Human-WithValidation-150.ipynb","version":"0.3.2","provenance":[{"file_id":"1bNhoDuxgvvq_cxjQCNw4nVKL3UPWZ0dE","timestamp":1549908220201},{"file_id":"1Pjy7fKwEhiqYGXTIgERTLFN1BpjYpJRc","timestamp":1549904042321},{"file_id":"13vKUPc7zaldx3RoMrDtHX08F8KJHYncx","timestamp":1549861587694}],"collapsed_sections":[]},"kernelspec":{"name":"python3","display_name":"Python 3"},"accelerator":"GPU"},"cells":[{"metadata":{"colab_type":"code","id":"RXZT2UsyIVe_","outputId":"acb843be-97b7-42aa-c675-f915afcaeb6c","executionInfo":{"status":"ok","timestamp":1549908367232,"user_tz":480,"elapsed":2920,"user":{"displayName":"Laurence Moroney","photoUrl":"https://lh4.googleusercontent.com/-wUzpekukCVw/AAAAAAAAAAI/AAAAAAAAAHw/pQPstOOJqqE/s64/photo.jpg","userId":"17858265307580721507"}},"colab":{"base_uri":"https://localhost:8080/","height":206}},"cell_type":"code","source":["!wget --no-check-certificate \\\n","    https://storage.googleapis.com/laurencemoroney-blog.appspot.com/horse-or-human.zip \\\n","    -O /tmp/horse-or-human.zip"],"execution_count":1,"outputs":[{"output_type":"stream","text":["--2019-02-11 18:06:05--  https://storage.googleapis.com/laurencemoroney-blog.appspot.com/horse-or-human.zip\n","Resolving storage.googleapis.com (storage.googleapis.com)... 74.125.20.128, 2607:f8b0:400e:c09::80\n","Connecting to storage.googleapis.com (storage.googleapis.com)|74.125.20.128|:443... connected.\n","HTTP request sent, awaiting response... 200 OK\n","Length: 149593353 (143M) [application/zip]\n","Saving to: ‘/tmp/horse-or-human.zip’\n","\n","/tmp/horse-or-human 100%[===================>] 142.66M   123MB/s    in 1.2s    \n","\n","2019-02-11 18:06:06 (123 MB/s) - ‘/tmp/horse-or-human.zip’ saved [149593353/149593353]\n","\n"],"name":"stdout"}]},{"metadata":{"id":"0mLij6qde6Ox","colab_type":"code","outputId":"52891eae-998f-4d85-9983-8133e9fee36d","executionInfo":{"status":"ok","timestamp":1549908370052,"user_tz":480,"elapsed":1424,"user":{"displayName":"Laurence Moroney","photoUrl":"https://lh4.googleusercontent.com/-wUzpekukCVw/AAAAAAAAAAI/AAAAAAAAAHw/pQPstOOJqqE/s64/photo.jpg","userId":"17858265307580721507"}},"colab":{"base_uri":"https://localhost:8080/","height":206}},"cell_type":"code","source":["!wget --no-check-certificate \\\n","    https://storage.googleapis.com/laurencemoroney-blog.appspot.com/validation-horse-or-human.zip \\\n","    -O /tmp/validation-horse-or-human.zip"],"execution_count":2,"outputs":[{"output_type":"stream","text":["--2019-02-11 18:06:09--  https://storage.googleapis.com/laurencemoroney-blog.appspot.com/validation-horse-or-human.zip\n","Resolving storage.googleapis.com (storage.googleapis.com)... 74.125.20.128, 2607:f8b0:400e:c08::80\n","Connecting to storage.googleapis.com (storage.googleapis.com)|74.125.20.128|:443... connected.\n","HTTP request sent, awaiting response... 200 OK\n","Length: 11485819 (11M) [application/zip]\n","Saving to: ‘/tmp/validation-horse-or-human.zip’\n","\n","/tmp/validation-hor 100%[===================>]  10.95M  65.1MB/s    in 0.2s    \n","\n","2019-02-11 18:06:09 (65.1 MB/s) - ‘/tmp/validation-horse-or-human.zip’ saved [11485819/11485819]\n","\n"],"name":"stdout"}]},{"metadata":{"id":"9brUxyTpYZHy","colab_type":"text"},"cell_type":"markdown","source":["The following python code will use the OS library to use Operating System libraries, giving you access to the file system, and the zipfile library allowing you to unzip the data. "]},{"metadata":{"colab_type":"code","id":"PLy3pthUS0D2","colab":{}},"cell_type":"code","source":["import os\n","import zipfile\n","\n","local_zip = '/tmp/horse-or-human.zip'\n","zip_ref = zipfile.ZipFile(local_zip, 'r')\n","zip_ref.extractall('/tmp/horse-or-human')\n","local_zip = '/tmp/validation-horse-or-human.zip'\n","zip_ref = zipfile.ZipFile(local_zip, 'r')\n","zip_ref.extractall('/tmp/validation-horse-or-human')\n","zip_ref.close()"],"execution_count":0,"outputs":[]},{"metadata":{"colab_type":"text","id":"o-qUPyfO7Qr8"},"cell_type":"markdown","source":["The contents of the .zip are extracted to the base directory `/tmp/horse-or-human`, which in turn each contain `filtered-horses` and `filtered-humans` subdirectories.\n","\n","In short: The training set is the data that is used to tell the neural network model that 'this is what a horse looks like', 'this is what a human looks like' etc. \n","\n","One thing to pay attention to in this sample: We do not explicitly label the images as horses or humans. If you remember with the handwriting example earlier, we had labelled 'this is a 1', 'this is a 7' etc.  Later you'll see something called an ImageGenerator being used -- and this is coded to read images from subdirectories, and automatically label them from the name of that subdirectory. So, for example, you will have a 'training' directory containing a 'filtered-horses' directory and a 'filtered-humans' one. ImageGenerator will label the images appropriately for you, reducing a coding step. \n","\n","Let's define each of these directories:"]},{"metadata":{"colab_type":"code","id":"NR_M9nWN-K8B","colab":{}},"cell_type":"code","source":["# Directory with our training horse pictures\n","train_horse_dir = os.path.join('/tmp/horse-or-human/filtered-horses')\n","\n","# Directory with our training human pictures\n","train_human_dir = os.path.join('/tmp/horse-or-human/filtered-humans')\n","\n","# Directory with our training horse pictures\n","validation_horse_dir = os.path.join('/tmp/validation-horse-or-human/validation-horses')\n","\n","# Directory with our training human pictures\n","validation_human_dir = os.path.join('/tmp/validation-horse-or-human/validation-humans')"],"execution_count":0,"outputs":[]},{"metadata":{"colab_type":"text","id":"5oqBkNBJmtUv"},"cell_type":"markdown","source":["## Building a Small Model from Scratch\n","\n","But before we continue, let's start defining the model:\n","\n","Step 1 will be to import tensorflow."]},{"metadata":{"id":"qvfZg3LQbD-5","colab_type":"code","colab":{}},"cell_type":"code","source":["import tensorflow as tf"],"execution_count":0,"outputs":[]},{"metadata":{"colab_type":"text","id":"BnhYCP4tdqjC"},"cell_type":"markdown","source":["We then add convolutional layers as in the previous example, and flatten the final result to feed into the densely connected layers."]},{"metadata":{"id":"gokG5HKpdtzm","colab_type":"text"},"cell_type":"markdown","source":["Finally we add the densely connected layers. \n","\n","Note that because we are facing a two-class classification problem, i.e. a *binary classification problem*, we will end our network with a [*sigmoid* activation](https://wikipedia.org/wiki/Sigmoid_function), so that the output of our network will be a single scalar between 0 and 1, encoding the probability that the current image is class 1 (as opposed to class 0)."]},{"metadata":{"id":"PixZ2s5QbYQ3","colab_type":"code","outputId":"17e65df8-01a6-4df5-d2d9-273b56f2808e","executionInfo":{"status":"ok","timestamp":1549908388497,"user_tz":480,"elapsed":441,"user":{"displayName":"Laurence Moroney","photoUrl":"https://lh4.googleusercontent.com/-wUzpekukCVw/AAAAAAAAAAI/AAAAAAAAAHw/pQPstOOJqqE/s64/photo.jpg","userId":"17858265307580721507"}},"colab":{"base_uri":"https://localhost:8080/","height":88}},"cell_type":"code","source":["model = tf.keras.models.Sequential([\n","    # Note the input shape is the desired size of the image 150x150 with 3 bytes color\n","    # This is the first convolution\n","    tf.keras.layers.Conv2D(16, (3,3), activation='relu', input_shape=(150, 150, 3)),\n","    tf.keras.layers.MaxPooling2D(2, 2),\n","    # The second convolution\n","    tf.keras.layers.Conv2D(32, (3,3), activation='relu'),\n","    tf.keras.layers.MaxPooling2D(2,2),\n","    # The third convolution\n","    tf.keras.layers.Conv2D(64, (3,3), activation='relu'),\n","    tf.keras.layers.MaxPooling2D(2,2),\n","    # The fourth convolution\n","    #tf.keras.layers.Conv2D(64, (3,3), activation='relu'),\n","    #tf.keras.layers.MaxPooling2D(2,2),\n","    # The fifth convolution\n","    #tf.keras.layers.Conv2D(64, (3,3), activation='relu'),\n","    #tf.keras.layers.MaxPooling2D(2,2),\n","    # Flatten the results to feed into a DNN\n","    tf.keras.layers.Flatten(),\n","    # 512 neuron hidden layer\n","    tf.keras.layers.Dense(512, activation='relu'),\n","    # Only 1 output neuron. It will contain a value from 0-1 where 0 for 1 class ('horses') and 1 for the other ('humans')\n","    tf.keras.layers.Dense(1, activation='sigmoid')\n","])"],"execution_count":6,"outputs":[{"output_type":"stream","text":["WARNING:tensorflow:From /usr/local/lib/python3.6/dist-packages/tensorflow/python/ops/resource_variable_ops.py:435: colocate_with (from tensorflow.python.framework.ops) is deprecated and will be removed in a future version.\n","Instructions for updating:\n","Colocations handled automatically by placer.\n"],"name":"stdout"}]},{"metadata":{"colab_type":"text","id":"s9EaFDP5srBa"},"cell_type":"markdown","source":["The model.summary() method call prints a summary of the NN "]},{"metadata":{"colab_type":"code","id":"7ZKj8392nbgP","outputId":"f7f2a9d3-6551-4209-9da1-5229c8d7f5d1","executionInfo":{"status":"ok","timestamp":1549908392560,"user_tz":480,"elapsed":318,"user":{"displayName":"Laurence Moroney","photoUrl":"https://lh4.googleusercontent.com/-wUzpekukCVw/AAAAAAAAAAI/AAAAAAAAAHw/pQPstOOJqqE/s64/photo.jpg","userId":"17858265307580721507"}},"colab":{"base_uri":"https://localhost:8080/","height":446}},"cell_type":"code","source":["model.summary()"],"execution_count":7,"outputs":[{"output_type":"stream","text":["_________________________________________________________________\n","Layer (type)                 Output Shape              Param #   \n","=================================================================\n","conv2d (Conv2D)              (None, 148, 148, 16)      448       \n","_________________________________________________________________\n","max_pooling2d (MaxPooling2D) (None, 74, 74, 16)        0         \n","_________________________________________________________________\n","conv2d_1 (Conv2D)            (None, 72, 72, 32)        4640      \n","_________________________________________________________________\n","max_pooling2d_1 (MaxPooling2 (None, 36, 36, 32)        0         \n","_________________________________________________________________\n","conv2d_2 (Conv2D)            (None, 34, 34, 64)        18496     \n","_________________________________________________________________\n","max_pooling2d_2 (MaxPooling2 (None, 17, 17, 64)        0         \n","_________________________________________________________________\n","flatten (Flatten)            (None, 18496)             0         \n","_________________________________________________________________\n","dense (Dense)                (None, 512)               9470464   \n","_________________________________________________________________\n","dense_1 (Dense)              (None, 1)                 513       \n","=================================================================\n","Total params: 9,494,561\n","Trainable params: 9,494,561\n","Non-trainable params: 0\n","_________________________________________________________________\n"],"name":"stdout"}]},{"metadata":{"colab_type":"text","id":"DmtkTn06pKxF"},"cell_type":"markdown","source":["The \"output shape\" column shows how the size of your feature map evolves in each successive layer. The convolution layers reduce the size of the feature maps by a bit due to padding, and each pooling layer halves the dimensions."]},{"metadata":{"colab_type":"text","id":"PEkKSpZlvJXA"},"cell_type":"markdown","source":["Next, we'll configure the specifications for model training. We will train our model with the `binary_crossentropy` loss, because it's a binary classification problem and our final activation is a sigmoid. (For a refresher on loss metrics, see the [Machine Learning Crash Course](https://developers.google.com/machine-learning/crash-course/descending-into-ml/video-lecture).) We will use the `rmsprop` optimizer with a learning rate of `0.001`. During training, we will want to monitor classification accuracy.\n","\n","**NOTE**: In this case, using the [RMSprop optimization algorithm](https://wikipedia.org/wiki/Stochastic_gradient_descent#RMSProp) is preferable to [stochastic gradient descent](https://developers.google.com/machine-learning/glossary/#SGD) (SGD), because RMSprop automates learning-rate tuning for us. (Other optimizers, such as [Adam](https://wikipedia.org/wiki/Stochastic_gradient_descent#Adam) and [Adagrad](https://developers.google.com/machine-learning/glossary/#AdaGrad), also automatically adapt the learning rate during training, and would work equally well here.)"]},{"metadata":{"colab_type":"code","id":"8DHWhFP_uhq3","colab":{}},"cell_type":"code","source":["from tensorflow.keras.optimizers import RMSprop\n","\n","model.compile(loss='binary_crossentropy',\n","              optimizer=RMSprop(lr=0.001),\n","              metrics=['acc'])"],"execution_count":0,"outputs":[]},{"metadata":{"colab_type":"text","id":"Sn9m9D3UimHM"},"cell_type":"markdown","source":["### Data Preprocessing\n","\n","Let's set up data generators that will read pictures in our source folders, convert them to `float32` tensors, and feed them (with their labels) to our network. We'll have one generator for the training images and one for the validation images. Our generators will yield batches of images of size 300x300 and their labels (binary).\n","\n","As you may already know, data that goes into neural networks should usually be normalized in some way to make it more amenable to processing by the network. (It is uncommon to feed raw pixels into a convnet.) In our case, we will preprocess our images by normalizing the pixel values to be in the `[0, 1]` range (originally all values are in the `[0, 255]` range).\n","\n","In Keras this can be done via the `keras.preprocessing.image.ImageDataGenerator` class using the `rescale` parameter. This `ImageDataGenerator` class allows you to instantiate generators of augmented image batches (and their labels) via `.flow(data, labels)` or `.flow_from_directory(directory)`. These generators can then be used with the Keras model methods that accept data generators as inputs: `fit_generator`, `evaluate_generator`, and `predict_generator`."]},{"metadata":{"colab_type":"code","id":"ClebU9NJg99G","outputId":"556c1d09-6125-460f-917a-b4d4d08dcebd","executionInfo":{"status":"ok","timestamp":1549908409693,"user_tz":480,"elapsed":454,"user":{"displayName":"Laurence Moroney","photoUrl":"https://lh4.googleusercontent.com/-wUzpekukCVw/AAAAAAAAAAI/AAAAAAAAAHw/pQPstOOJqqE/s64/photo.jpg","userId":"17858265307580721507"}},"colab":{"base_uri":"https://localhost:8080/","height":51}},"cell_type":"code","source":["from tensorflow.keras.preprocessing.image import ImageDataGenerator\n","\n","# All images will be rescaled by 1./255\n","train_datagen = ImageDataGenerator(rescale=1/255)\n","validation_datagen = ImageDataGenerator(rescale=1/255)\n","\n","# Flow training images in batches of 128 using train_datagen generator\n","train_generator = train_datagen.flow_from_directory(\n","        '/tmp/horse-or-human/',  # This is the source directory for training images\n","        target_size=(150, 150),  # All images will be resized to 150x150\n","        batch_size=128,\n","        # Since we use binary_crossentropy loss, we need binary labels\n","        class_mode='binary')\n","\n","# Flow training images in batches of 128 using train_datagen generator\n","validation_generator = validation_datagen.flow_from_directory(\n","        '/tmp/validation-horse-or-human/',  # This is the source directory for training images\n","        target_size=(150, 150),  # All images will be resized to 150x150\n","        batch_size=32,\n","        # Since we use binary_crossentropy loss, we need binary labels\n","        class_mode='binary')"],"execution_count":9,"outputs":[{"output_type":"stream","text":["Found 1027 images belonging to 2 classes.\n","Found 256 images belonging to 2 classes.\n"],"name":"stdout"}]},{"metadata":{"colab_type":"text","id":"mu3Jdwkjwax4"},"cell_type":"markdown","source":["### Training\n","Let's train for 15 epochs -- this may take a few minutes to run.\n","\n","Do note the values per epoch.\n","\n","The Loss and Accuracy are a great indication of progress of training. It's making a guess as to the classification of the training data, and then measuring it against the known label, calculating the result. Accuracy is the portion of correct guesses. "]},{"metadata":{"colab_type":"code","id":"Fb1_lgobv81m","outputId":"898d7e34-0d9d-4b78-90b2-45638bc061b5","executionInfo":{"status":"ok","timestamp":1549908525339,"user_tz":480,"elapsed":112000,"user":{"displayName":"Laurence Moroney","photoUrl":"https://lh4.googleusercontent.com/-wUzpekukCVw/AAAAAAAAAAI/AAAAAAAAAHw/pQPstOOJqqE/s64/photo.jpg","userId":"17858265307580721507"}},"colab":{"base_uri":"https://localhost:8080/","height":860}},"cell_type":"code","source":["history = model.fit_generator(\n","      train_generator,\n","      steps_per_epoch=8,  \n","      epochs=15,\n","      verbose=1,\n","      validation_data = validation_generator,\n","      validation_steps=8)"],"execution_count":10,"outputs":[{"output_type":"stream","text":["WARNING:tensorflow:From /usr/local/lib/python3.6/dist-packages/tensorflow/python/ops/math_ops.py:3066: to_int32 (from tensorflow.python.ops.math_ops) is deprecated and will be removed in a future version.\n","Instructions for updating:\n","Use tf.cast instead.\n","Epoch 1/15\n","8/8 [==============================] - 1s 181ms/step - loss: 0.6093 - acc: 0.6523\n","9/9 [==============================] - 8s 876ms/step - loss: 1.5314 - acc: 0.5998 - val_loss: 0.6093 - val_acc: 0.6523\n","Epoch 2/15\n","8/8 [==============================] - 1s 178ms/step - loss: 0.8493 - acc: 0.7578\n","9/9 [==============================] - 7s 810ms/step - loss: 0.4189 - acc: 0.8306 - val_loss: 0.8493 - val_acc: 0.7578\n","Epoch 3/15\n","8/8 [==============================] - 1s 172ms/step - loss: 2.0065 - acc: 0.7305\n","9/9 [==============================] - 7s 811ms/step - loss: 0.1664 - acc: 0.9426 - val_loss: 2.0065 - val_acc: 0.7305\n","Epoch 4/15\n","8/8 [==============================] - 1s 180ms/step - loss: 1.5915 - acc: 0.8008\n","9/9 [==============================] - 7s 827ms/step - loss: 0.0786 - acc: 0.9669 - val_loss: 1.5915 - val_acc: 0.8008\n","Epoch 5/15\n","8/8 [==============================] - 1s 176ms/step - loss: 2.3525 - acc: 0.7344\n","9/9 [==============================] - 7s 804ms/step - loss: 0.0396 - acc: 0.9883 - val_loss: 2.3525 - val_acc: 0.7344\n","Epoch 6/15\n","8/8 [==============================] - 1s 177ms/step - loss: 1.1861 - acc: 0.8125\n","9/9 [==============================] - 7s 810ms/step - loss: 0.3273 - acc: 0.9007 - val_loss: 1.1861 - val_acc: 0.8125\n","Epoch 7/15\n","8/8 [==============================] - 1s 176ms/step - loss: 1.3978 - acc: 0.8438\n","9/9 [==============================] - 7s 801ms/step - loss: 0.0430 - acc: 0.9864 - val_loss: 1.3978 - val_acc: 0.8438\n","Epoch 8/15\n","8/8 [==============================] - 1s 176ms/step - loss: 1.5442 - acc: 0.8438\n","9/9 [==============================] - 7s 809ms/step - loss: 0.0198 - acc: 0.9942 - val_loss: 1.5442 - val_acc: 0.8438\n","Epoch 9/15\n","8/8 [==============================] - 1s 168ms/step - loss: 1.5385 - acc: 0.8516\n","9/9 [==============================] - 7s 806ms/step - loss: 0.0144 - acc: 0.9981 - val_loss: 1.5385 - val_acc: 0.8516\n","Epoch 10/15\n","8/8 [==============================] - 1s 176ms/step - loss: 1.5228 - acc: 0.8594\n","9/9 [==============================] - 7s 810ms/step - loss: 0.0052 - acc: 1.0000 - val_loss: 1.5228 - val_acc: 0.8594\n","Epoch 11/15\n","8/8 [==============================] - 1s 179ms/step - loss: 1.1527 - acc: 0.8750\n","9/9 [==============================] - 7s 812ms/step - loss: 0.0033 - acc: 1.0000 - val_loss: 1.1527 - val_acc: 0.8750\n","Epoch 12/15\n","8/8 [==============================] - 1s 173ms/step - loss: 1.6821 - acc: 0.8359\n","9/9 [==============================] - 7s 784ms/step - loss: 0.1890 - acc: 0.9435 - val_loss: 1.6821 - val_acc: 0.8359\n","Epoch 13/15\n","8/8 [==============================] - 1s 178ms/step - loss: 1.8283 - acc: 0.8359\n","9/9 [==============================] - 7s 820ms/step - loss: 0.0039 - acc: 1.0000 - val_loss: 1.8283 - val_acc: 0.8359\n","Epoch 14/15\n","8/8 [==============================] - 1s 179ms/step - loss: 2.0000 - acc: 0.8203\n","9/9 [==============================] - 7s 799ms/step - loss: 0.0027 - acc: 1.0000 - val_loss: 2.0000 - val_acc: 0.8203\n","Epoch 15/15\n","8/8 [==============================] - 1s 176ms/step - loss: 0.7265 - acc: 0.8398\n","9/9 [==============================] - 7s 786ms/step - loss: 0.0069 - acc: 0.9990 - val_loss: 0.7265 - val_acc: 0.8398\n"],"name":"stdout"}]},{"metadata":{"id":"o6vSHzPR2ghH","colab_type":"text"},"cell_type":"markdown","source":["###Running the Model\n","\n","Let's now take a look at actually running a prediction using the model. This code will allow you to choose 1 or more files from your file system, it will then upload them, and run them through the model, giving an indication of whether the object is a horse or a human."]},{"metadata":{"id":"DoWp43WxJDNT","colab_type":"code","outputId":"49b79085-3189-43c5-a13f-6058abebbe06","executionInfo":{"status":"ok","timestamp":1549908583241,"user_tz":480,"elapsed":7478,"user":{"displayName":"Laurence Moroney","photoUrl":"https://lh4.googleusercontent.com/-wUzpekukCVw/AAAAAAAAAAI/AAAAAAAAAHw/pQPstOOJqqE/s64/photo.jpg","userId":"17858265307580721507"}},"colab":{"resources":{"http://localhost:8080/nbextensions/google.colab/files.js":{"data":"Ly8gQ29weXJpZ2h0IDIwMTcgR29vZ2xlIExMQwovLwovLyBMaWNlbnNlZCB1bmRlciB0aGUgQXBhY2hlIExpY2Vuc2UsIFZlcnNpb24gMi4wICh0aGUgIkxpY2Vuc2UiKTsKLy8geW91IG1heSBub3QgdXNlIHRoaXMgZmlsZSBleGNlcHQgaW4gY29tcGxpYW5jZSB3aXRoIHRoZSBMaWNlbnNlLgovLyBZb3UgbWF5IG9idGFpbiBhIGNvcHkgb2YgdGhlIExpY2Vuc2UgYXQKLy8KLy8gICAgICBodHRwOi8vd3d3LmFwYWNoZS5vcmcvbGljZW5zZXMvTElDRU5TRS0yLjAKLy8KLy8gVW5sZXNzIHJlcXVpcmVkIGJ5IGFwcGxpY2FibGUgbGF3IG9yIGFncmVlZCB0byBpbiB3cml0aW5nLCBzb2Z0d2FyZQovLyBkaXN0cmlidXRlZCB1bmRlciB0aGUgTGljZW5zZSBpcyBkaXN0cmlidXRlZCBvbiBhbiAiQVMgSVMiIEJBU0lTLAovLyBXSVRIT1VUIFdBUlJBTlRJRVMgT1IgQ09ORElUSU9OUyBPRiBBTlkgS0lORCwgZWl0aGVyIGV4cHJlc3Mgb3IgaW1wbGllZC4KLy8gU2VlIHRoZSBMaWNlbnNlIGZvciB0aGUgc3BlY2lmaWMgbGFuZ3VhZ2UgZ292ZXJuaW5nIHBlcm1pc3Npb25zIGFuZAovLyBsaW1pdGF0aW9ucyB1bmRlciB0aGUgTGljZW5zZS4KCi8qKgogKiBAZmlsZW92ZXJ2aWV3IEhlbHBlcnMgZm9yIGdvb2dsZS5jb2xhYiBQeXRob24gbW9kdWxlLgogKi8KKGZ1bmN0aW9uKHNjb3BlKSB7CmZ1bmN0aW9uIHNwYW4odGV4dCwgc3R5bGVBdHRyaWJ1dGVzID0ge30pIHsKICBjb25zdCBlbGVtZW50ID0gZG9jdW1lbnQuY3JlYXRlRWxlbWVudCgnc3BhbicpOwogIGVsZW1lbnQudGV4dENvbnRlbnQgPSB0ZXh0OwogIGZvciAoY29uc3Qga2V5IG9mIE9iamVjdC5rZXlzKHN0eWxlQXR0cmlidXRlcykpIHsKICAgIGVsZW1lbnQuc3R5bGVba2V5XSA9IHN0eWxlQXR0cmlidXRlc1trZXldOwogIH0KICByZXR1cm4gZWxlbWVudDsKfQoKLy8gTWF4IG51bWJlciBvZiBieXRlcyB3aGljaCB3aWxsIGJlIHVwbG9hZGVkIGF0IGEgdGltZS4KY29uc3QgTUFYX1BBWUxPQURfU0laRSA9IDEwMCAqIDEwMjQ7Ci8vIE1heCBhbW91bnQgb2YgdGltZSB0byBibG9jayB3YWl0aW5nIGZvciB0aGUgdXNlci4KY29uc3QgRklMRV9DSEFOR0VfVElNRU9VVF9NUyA9IDMwICogMTAwMDsKCmZ1bmN0aW9uIF91cGxvYWRGaWxlcyhpbnB1dElkLCBvdXRwdXRJZCkgewogIGNvbnN0IHN0ZXBzID0gdXBsb2FkRmlsZXNTdGVwKGlucHV0SWQsIG91dHB1dElkKTsKICBjb25zdCBvdXRwdXRFbGVtZW50ID0gZG9jdW1lbnQuZ2V0RWxlbWVudEJ5SWQob3V0cHV0SWQpOwogIC8vIENhY2hlIHN0ZXBzIG9uIHRoZSBvdXRwdXRFbGVtZW50IHRvIG1ha2UgaXQgYXZhaWxhYmxlIGZvciB0aGUgbmV4dCBjYWxsCiAgLy8gdG8gdXBsb2FkRmlsZXNDb250aW51ZSBmcm9tIFB5dGhvbi4KICBvdXRwdXRFbGVtZW50LnN0ZXBzID0gc3RlcHM7CgogIHJldHVybiBfdXBsb2FkRmlsZXNDb250aW51ZShvdXRwdXRJZCk7Cn0KCi8vIFRoaXMgaXMgcm91Z2hseSBhbiBhc3luYyBnZW5lcmF0b3IgKG5vdCBzdXBwb3J0ZWQgaW4gdGhlIGJyb3dzZXIgeWV0KSwKLy8gd2hlcmUgdGhlcmUgYXJlIG11bHRpcGxlIGFzeW5jaHJvbm91cyBzdGVwcyBhbmQgdGhlIFB5dGhvbiBzaWRlIGlzIGdvaW5nCi8vIHRvIHBvbGwgZm9yIGNvbXBsZXRpb24gb2YgZWFjaCBzdGVwLgovLyBUaGlzIHVzZXMgYSBQcm9taXNlIHRvIGJsb2NrIHRoZSBweXRob24gc2lkZSBvbiBjb21wbGV0aW9uIG9mIGVhY2ggc3RlcCwKLy8gdGhlbiBwYXNzZXMgdGhlIHJlc3VsdCBvZiB0aGUgcHJldmlvdXMgc3RlcCBhcyB0aGUgaW5wdXQgdG8gdGhlIG5leHQgc3RlcC4KZnVuY3Rpb24gX3VwbG9hZEZpbGVzQ29udGludWUob3V0cHV0SWQpIHsKICBjb25zdCBvdXRwdXRFbGVtZW50ID0gZG9jdW1lbnQuZ2V0RWxlbWVudEJ5SWQob3V0cHV0SWQpOwogIGNvbnN0IHN0ZXBzID0gb3V0cHV0RWxlbWVudC5zdGVwczsKCiAgY29uc3QgbmV4dCA9IHN0ZXBzLm5leHQob3V0cHV0RWxlbWVudC5sYXN0UHJvbWlzZVZhbHVlKTsKICByZXR1cm4gUHJvbWlzZS5yZXNvbHZlKG5leHQudmFsdWUucHJvbWlzZSkudGhlbigodmFsdWUpID0+IHsKICAgIC8vIENhY2hlIHRoZSBsYXN0IHByb21pc2UgdmFsdWUgdG8gbWFrZSBpdCBhdmFpbGFibGUgdG8gdGhlIG5leHQKICAgIC8vIHN0ZXAgb2YgdGhlIGdlbmVyYXRvci4KICAgIG91dHB1dEVsZW1lbnQubGFzdFByb21pc2VWYWx1ZSA9IHZhbHVlOwogICAgcmV0dXJuIG5leHQudmFsdWUucmVzcG9uc2U7CiAgfSk7Cn0KCi8qKgogKiBHZW5lcmF0b3IgZnVuY3Rpb24gd2hpY2ggaXMgY2FsbGVkIGJldHdlZW4gZWFjaCBhc3luYyBzdGVwIG9mIHRoZSB1cGxvYWQKICogcHJvY2Vzcy4KICogQHBhcmFtIHtzdHJpbmd9IGlucHV0SWQgRWxlbWVudCBJRCBvZiB0aGUgaW5wdXQgZmlsZSBwaWNrZXIgZWxlbWVudC4KICogQHBhcmFtIHtzdHJpbmd9IG91dHB1dElkIEVsZW1lbnQgSUQgb2YgdGhlIG91dHB1dCBkaXNwbGF5LgogKiBAcmV0dXJuIHshSXRlcmFibGU8IU9iamVjdD59IEl0ZXJhYmxlIG9mIG5leHQgc3RlcHMuCiAqLwpmdW5jdGlvbiogdXBsb2FkRmlsZXNTdGVwKGlucHV0SWQsIG91dHB1dElkKSB7CiAgY29uc3QgaW5wdXRFbGVtZW50ID0gZG9jdW1lbnQuZ2V0RWxlbWVudEJ5SWQoaW5wdXRJZCk7CiAgaW5wdXRFbGVtZW50LmRpc2FibGVkID0gZmFsc2U7CgogIGNvbnN0IG91dHB1dEVsZW1lbnQgPSBkb2N1bWVudC5nZXRFbGVtZW50QnlJZChvdXRwdXRJZCk7CiAgb3V0cHV0RWxlbWVudC5pbm5lckhUTUwgPSAnJzsKCiAgY29uc3QgcGlja2VkUHJvbWlzZSA9IG5ldyBQcm9taXNlKChyZXNvbHZlKSA9PiB7CiAgICBpbnB1dEVsZW1lbnQuYWRkRXZlbnRMaXN0ZW5lcignY2hhbmdlJywgKGUpID0+IHsKICAgICAgcmVzb2x2ZShlLnRhcmdldC5maWxlcyk7CiAgICB9KTsKICB9KTsKCiAgY29uc3QgY2FuY2VsID0gZG9jdW1lbnQuY3JlYXRlRWxlbWVudCgnYnV0dG9uJyk7CiAgaW5wdXRFbGVtZW50LnBhcmVudEVsZW1lbnQuYXBwZW5kQ2hpbGQoY2FuY2VsKTsKICBjYW5jZWwudGV4dENvbnRlbnQgPSAnQ2FuY2VsIHVwbG9hZCc7CiAgY29uc3QgY2FuY2VsUHJvbWlzZSA9IG5ldyBQcm9taXNlKChyZXNvbHZlKSA9PiB7CiAgICBjYW5jZWwub25jbGljayA9ICgpID0+IHsKICAgICAgcmVzb2x2ZShudWxsKTsKICAgIH07CiAgfSk7CgogIC8vIENhbmNlbCB1cGxvYWQgaWYgdXNlciBoYXNuJ3QgcGlja2VkIGFueXRoaW5nIGluIHRpbWVvdXQuCiAgY29uc3QgdGltZW91dFByb21pc2UgPSBuZXcgUHJvbWlzZSgocmVzb2x2ZSkgPT4gewogICAgc2V0VGltZW91dCgoKSA9PiB7CiAgICAgIHJlc29sdmUobnVsbCk7CiAgICB9LCBGSUxFX0NIQU5HRV9USU1FT1VUX01TKTsKICB9KTsKCiAgLy8gV2FpdCBmb3IgdGhlIHVzZXIgdG8gcGljayB0aGUgZmlsZXMuCiAgY29uc3QgZmlsZXMgPSB5aWVsZCB7CiAgICBwcm9taXNlOiBQcm9taXNlLnJhY2UoW3BpY2tlZFByb21pc2UsIHRpbWVvdXRQcm9taXNlLCBjYW5jZWxQcm9taXNlXSksCiAgICByZXNwb25zZTogewogICAgICBhY3Rpb246ICdzdGFydGluZycsCiAgICB9CiAgfTsKCiAgaWYgKCFmaWxlcykgewogICAgcmV0dXJuIHsKICAgICAgcmVzcG9uc2U6IHsKICAgICAgICBhY3Rpb246ICdjb21wbGV0ZScsCiAgICAgIH0KICAgIH07CiAgfQoKICBjYW5jZWwucmVtb3ZlKCk7CgogIC8vIERpc2FibGUgdGhlIGlucHV0IGVsZW1lbnQgc2luY2UgZnVydGhlciBwaWNrcyBhcmUgbm90IGFsbG93ZWQuCiAgaW5wdXRFbGVtZW50LmRpc2FibGVkID0gdHJ1ZTsKCiAgZm9yIChjb25zdCBmaWxlIG9mIGZpbGVzKSB7CiAgICBjb25zdCBsaSA9IGRvY3VtZW50LmNyZWF0ZUVsZW1lbnQoJ2xpJyk7CiAgICBsaS5hcHBlbmQoc3BhbihmaWxlLm5hbWUsIHtmb250V2VpZ2h0OiAnYm9sZCd9KSk7CiAgICBsaS5hcHBlbmQoc3BhbigKICAgICAgICBgKCR7ZmlsZS50eXBlIHx8ICduL2EnfSkgLSAke2ZpbGUuc2l6ZX0gYnl0ZXMsIGAgKwogICAgICAgIGBsYXN0IG1vZGlmaWVkOiAkewogICAgICAgICAgICBmaWxlLmxhc3RNb2RpZmllZERhdGUgPyBmaWxlLmxhc3RNb2RpZmllZERhdGUudG9Mb2NhbGVEYXRlU3RyaW5nKCkgOgogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAnbi9hJ30gLSBgKSk7CiAgICBjb25zdCBwZXJjZW50ID0gc3BhbignMCUgZG9uZScpOwogICAgbGkuYXBwZW5kQ2hpbGQocGVyY2VudCk7CgogICAgb3V0cHV0RWxlbWVudC5hcHBlbmRDaGlsZChsaSk7CgogICAgY29uc3QgZmlsZURhdGFQcm9taXNlID0gbmV3IFByb21pc2UoKHJlc29sdmUpID0+IHsKICAgICAgY29uc3QgcmVhZGVyID0gbmV3IEZpbGVSZWFkZXIoKTsKICAgICAgcmVhZGVyLm9ubG9hZCA9IChlKSA9PiB7CiAgICAgICAgcmVzb2x2ZShlLnRhcmdldC5yZXN1bHQpOwogICAgICB9OwogICAgICByZWFkZXIucmVhZEFzQXJyYXlCdWZmZXIoZmlsZSk7CiAgICB9KTsKICAgIC8vIFdhaXQgZm9yIHRoZSBkYXRhIHRvIGJlIHJlYWR5LgogICAgbGV0IGZpbGVEYXRhID0geWllbGQgewogICAgICBwcm9taXNlOiBmaWxlRGF0YVByb21pc2UsCiAgICAgIHJlc3BvbnNlOiB7CiAgICAgICAgYWN0aW9uOiAnY29udGludWUnLAogICAgICB9CiAgICB9OwoKICAgIC8vIFVzZSBhIGNodW5rZWQgc2VuZGluZyB0byBhdm9pZCBtZXNzYWdlIHNpemUgbGltaXRzLiBTZWUgYi82MjExNTY2MC4KICAgIGxldCBwb3NpdGlvbiA9IDA7CiAgICB3aGlsZSAocG9zaXRpb24gPCBmaWxlRGF0YS5ieXRlTGVuZ3RoKSB7CiAgICAgIGNvbnN0IGxlbmd0aCA9IE1hdGgubWluKGZpbGVEYXRhLmJ5dGVMZW5ndGggLSBwb3NpdGlvbiwgTUFYX1BBWUxPQURfU0laRSk7CiAgICAgIGNvbnN0IGNodW5rID0gbmV3IFVpbnQ4QXJyYXkoZmlsZURhdGEsIHBvc2l0aW9uLCBsZW5ndGgpOwogICAgICBwb3NpdGlvbiArPSBsZW5ndGg7CgogICAgICBjb25zdCBiYXNlNjQgPSBidG9hKFN0cmluZy5mcm9tQ2hhckNvZGUuYXBwbHkobnVsbCwgY2h1bmspKTsKICAgICAgeWllbGQgewogICAgICAgIHJlc3BvbnNlOiB7CiAgICAgICAgICBhY3Rpb246ICdhcHBlbmQnLAogICAgICAgICAgZmlsZTogZmlsZS5uYW1lLAogICAgICAgICAgZGF0YTogYmFzZTY0LAogICAgICAgIH0sCiAgICAgIH07CiAgICAgIHBlcmNlbnQudGV4dENvbnRlbnQgPQogICAgICAgICAgYCR7TWF0aC5yb3VuZCgocG9zaXRpb24gLyBmaWxlRGF0YS5ieXRlTGVuZ3RoKSAqIDEwMCl9JSBkb25lYDsKICAgIH0KICB9CgogIC8vIEFsbCBkb25lLgogIHlpZWxkIHsKICAgIHJlc3BvbnNlOiB7CiAgICAgIGFjdGlvbjogJ2NvbXBsZXRlJywKICAgIH0KICB9Owp9CgpzY29wZS5nb29nbGUgPSBzY29wZS5nb29nbGUgfHwge307CnNjb3BlLmdvb2dsZS5jb2xhYiA9IHNjb3BlLmdvb2dsZS5jb2xhYiB8fCB7fTsKc2NvcGUuZ29vZ2xlLmNvbGFiLl9maWxlcyA9IHsKICBfdXBsb2FkRmlsZXMsCiAgX3VwbG9hZEZpbGVzQ29udGludWUsCn07Cn0pKHNlbGYpOwo=","ok":true,"headers":[["content-type","application/javascript"]],"status":200,"status_text":""}},"base_uri":"https://localhost:8080/","height":104}},"cell_type":"code","source":["import numpy as np\n","from google.colab import files\n","from keras.preprocessing import image\n","\n","uploaded = files.upload()\n","\n","for fn in uploaded.keys():\n"," \n","  # predicting images\n","  path = '/content/' + fn\n","  img = image.load_img(path, target_size=(150, 150))\n","  x = image.img_to_array(img)\n","  x = np.expand_dims(x, axis=0)\n","\n","  images = np.vstack([x])\n","  classes = model.predict(images, batch_size=10)\n","  print(classes[0])\n","  if classes[0]>0.5:\n","    print(fn + \" is a human\")\n","  else:\n","    print(fn + \" is a horse\")\n"," "],"execution_count":15,"outputs":[{"output_type":"display_data","data":{"text/html":["\n","     <input type=\"file\" id=\"files-1c60183f-5f52-4b2e-80a9-830d1e834dab\" name=\"files[]\" multiple disabled />\n","     <output id=\"result-1c60183f-5f52-4b2e-80a9-830d1e834dab\">\n","      Upload widget is only available when the cell has been executed in the\n","      current browser session. Please rerun this cell to enable.\n","      </output>\n","      <script src=\"/nbextensions/google.colab/files.js\"></script> "],"text/plain":["<IPython.core.display.HTML object>"]},"metadata":{"tags":[]}},{"output_type":"stream","text":["Saving woman-918784_640.jpg to woman-918784_640.jpg\n","[0.]\n","woman-918784_640.jpg is a horse\n"],"name":"stdout"}]},{"metadata":{"colab_type":"text","id":"-8EHQyWGDvWz"},"cell_type":"markdown","source":["### Visualizing Intermediate Representations\n","\n","To get a feel for what kind of features our convnet has learned, one fun thing to do is to visualize how an input gets transformed as it goes through the convnet.\n","\n","Let's pick a random image from the training set, and then generate a figure where each row is the output of a layer, and each image in the row is a specific filter in that output feature map. Rerun this cell to generate intermediate representations for a variety of training images."]},{"metadata":{"colab_type":"code","id":"-5tES8rXFjux","colab":{}},"cell_type":"code","source":["import numpy as np\n","import random\n","from tensorflow.keras.preprocessing.image import img_to_array, load_img\n","\n","# Let's define a new Model that will take an image as input, and will output\n","# intermediate representations for all layers in the previous model after\n","# the first.\n","successive_outputs = [layer.output for layer in model.layers[1:]]\n","#visualization_model = Model(img_input, successive_outputs)\n","visualization_model = tf.keras.models.Model(inputs = model.input, outputs = successive_outputs)\n","# Let's prepare a random input image from the training set.\n","horse_img_files = [os.path.join(train_horse_dir, f) for f in train_horse_names]\n","human_img_files = [os.path.join(train_human_dir, f) for f in train_human_names]\n","img_path = random.choice(horse_img_files + human_img_files)\n","\n","img = load_img(img_path, target_size=(300, 300))  # this is a PIL image\n","x = img_to_array(img)  # Numpy array with shape (150, 150, 3)\n","x = x.reshape((1,) + x.shape)  # Numpy array with shape (1, 150, 150, 3)\n","\n","# Rescale by 1/255\n","x /= 255\n","\n","# Let's run our image through our network, thus obtaining all\n","# intermediate representations for this image.\n","successive_feature_maps = visualization_model.predict(x)\n","\n","# These are the names of the layers, so can have them as part of our plot\n","layer_names = [layer.name for layer in model.layers]\n","\n","# Now let's display our representations\n","for layer_name, feature_map in zip(layer_names, successive_feature_maps):\n","  if len(feature_map.shape) == 4:\n","    # Just do this for the conv / maxpool layers, not the fully-connected layers\n","    n_features = feature_map.shape[-1]  # number of features in feature map\n","    # The feature map has shape (1, size, size, n_features)\n","    size = feature_map.shape[1]\n","    # We will tile our images in this matrix\n","    display_grid = np.zeros((size, size * n_features))\n","    for i in range(n_features):\n","      # Postprocess the feature to make it visually palatable\n","      x = feature_map[0, :, :, i]\n","      x -= x.mean()\n","      x /= x.std()\n","      x *= 64\n","      x += 128\n","      x = np.clip(x, 0, 255).astype('uint8')\n","      # We'll tile each filter into this big horizontal grid\n","      display_grid[:, i * size : (i + 1) * size] = x\n","    # Display the grid\n","    scale = 20. / n_features\n","    plt.figure(figsize=(scale * n_features, scale))\n","    plt.title(layer_name)\n","    plt.grid(False)\n","    plt.imshow(display_grid, aspect='auto', cmap='viridis')"],"execution_count":0,"outputs":[]},{"metadata":{"colab_type":"text","id":"tuqK2arJL0wo"},"cell_type":"markdown","source":["As you can see we go from the raw pixels of the images to increasingly abstract and compact representations. The representations downstream start highlighting what the network pays attention to, and they show fewer and fewer features being \"activated\"; most are set to zero. This is called \"sparsity.\" Representation sparsity is a key feature of deep learning.\n","\n","\n","These representations carry increasingly less information about the original pixels of the image, but increasingly refined information about the class of the image. You can think of a convnet (or a deep network in general) as an information distillation pipeline."]},{"metadata":{"colab_type":"text","id":"j4IBgYCYooGD"},"cell_type":"markdown","source":["## Clean Up\n","\n","Before running the next exercise, run the following cell to terminate the kernel and free memory resources:"]},{"metadata":{"colab_type":"code","id":"651IgjLyo-Jx","colab":{}},"cell_type":"code","source":["import os, signal\n","os.kill(os.getpid(), signal.SIGKILL)"],"execution_count":0,"outputs":[]}]}