{"metadata":{"kernelspec":{"language":"python","display_name":"Python 3","name":"python3"},"language_info":{"name":"python","version":"3.7.6","mimetype":"text/x-python","codemirror_mode":{"name":"ipython","version":3},"pygments_lexer":"ipython3","nbconvert_exporter":"python","file_extension":".py"},"kaggle":{"accelerator":"gpu","dataSources":[{"sourceId":13836,"databundleVersionId":1718836,"sourceType":"competition"},{"sourceId":1721243,"sourceType":"datasetVersion","datasetId":1008017}],"dockerImageVersionId":30034,"isInternetEnabled":false,"language":"python","sourceType":"notebook","isGpuEnabled":true}},"nbformat_minor":4,"nbformat":4,"cells":[{"cell_type":"markdown","source":"## Loading packages","metadata":{}},{"cell_type":"code","source":"import numpy as np\nimport pandas as pd\nimport matplotlib.pyplot as plt\nimport seaborn as sns\nimport datetime\n\nfrom sklearn.model_selection import train_test_split\nfrom sklearn.metrics import accuracy_score\nimport tensorflow as tf\nfrom tensorflow.keras import models, layers\nfrom tensorflow.keras.preprocessing import image\nfrom tensorflow.keras.preprocessing.image import ImageDataGenerator\nfrom tensorflow.keras.callbacks import ModelCheckpoint, EarlyStopping, ReduceLROnPlateau\nfrom tensorflow.keras.applications import EfficientNetB0\nfrom tensorflow.keras.optimizers import Adam\n\n# ignoring warnings\nimport warnings\nwarnings.simplefilter(\"ignore\")\n\nimport os, cv2, json\nfrom PIL import Image","metadata":{"_uuid":"d629ff2d2480ee46fbb7e2d37f6b5fab8052498a","_cell_guid":"79c7e3d0-c299-4dcb-8224-4455121ee9b0","_kg_hide-input":true,"execution":{"iopub.status.busy":"2025-04-09T13:47:15.083398Z","iopub.execute_input":"2025-04-09T13:47:15.083705Z","iopub.status.idle":"2025-04-09T13:47:15.090002Z","shell.execute_reply.started":"2025-04-09T13:47:15.083679Z","shell.execute_reply":"2025-04-09T13:47:15.089083Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"## Work directory","metadata":{}},{"cell_type":"code","source":"WORK_DIR = '../input/cassava-leaf-disease-classification'\nos.listdir(WORK_DIR)","metadata":{"execution":{"iopub.status.busy":"2025-04-09T13:47:21.168682Z","iopub.execute_input":"2025-04-09T13:47:21.169005Z","iopub.status.idle":"2025-04-09T13:47:21.175805Z","shell.execute_reply.started":"2025-04-09T13:47:21.168976Z","shell.execute_reply":"2025-04-09T13:47:21.175040Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"# First look at the data","metadata":{}},{"cell_type":"code","source":"print('Train images: %d' %len(os.listdir(\n    os.path.join(WORK_DIR, \"train_images\"))))","metadata":{"_kg_hide-input":true,"execution":{"iopub.status.busy":"2025-04-09T13:47:24.384257Z","iopub.execute_input":"2025-04-09T13:47:24.384755Z","iopub.status.idle":"2025-04-09T13:47:24.396744Z","shell.execute_reply.started":"2025-04-09T13:47:24.384706Z","shell.execute_reply":"2025-04-09T13:47:24.396055Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"code","source":"with open(os.path.join(WORK_DIR, \"label_num_to_disease_map.json\")) as file:\n    print(json.dumps(json.loads(file.read()), indent=4))","metadata":{"_kg_hide-input":true,"execution":{"iopub.status.busy":"2025-04-09T13:47:27.864526Z","iopub.execute_input":"2025-04-09T13:47:27.864805Z","iopub.status.idle":"2025-04-09T13:47:27.871336Z","shell.execute_reply.started":"2025-04-09T13:47:27.864781Z","shell.execute_reply":"2025-04-09T13:47:27.870307Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"code","source":"train_labels = pd.read_csv(os.path.join(WORK_DIR, \"train.csv\"))\ntrain_labels.head()","metadata":{"execution":{"iopub.status.busy":"2025-04-09T13:47:31.482700Z","iopub.execute_input":"2025-04-09T13:47:31.483043Z","iopub.status.idle":"2025-04-09T13:47:31.505693Z","shell.execute_reply.started":"2025-04-09T13:47:31.483010Z","shell.execute_reply":"2025-04-09T13:47:31.504878Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"code","source":"sns.set_style(\"whitegrid\")\nfig, ax = plt.subplots(figsize = (6, 4))\n\nfor i in ['top', 'right', 'left']:\n    ax.spines[i].set_visible(False)\nax.spines['bottom'].set_color('black')\n\nsns.countplot(train_labels.label, edgecolor = 'black',\n              palette = reversed(sns.color_palette(\"viridis\", 5)))\nplt.xlabel('Classes', fontfamily = 'serif', size = 15)\nplt.ylabel('Count', fontfamily = 'serif', size = 15)\nplt.xticks(fontfamily = 'serif', size = 12)\nplt.yticks(fontfamily = 'serif', size = 12)\nax.grid(axis = 'y', linestyle = '--', alpha = 0.9)\nplt.show()","metadata":{"execution":{"iopub.status.busy":"2025-04-09T13:47:34.964566Z","iopub.execute_input":"2025-04-09T13:47:34.964879Z","iopub.status.idle":"2025-04-09T13:47:35.098577Z","shell.execute_reply.started":"2025-04-09T13:47:34.964854Z","shell.execute_reply":"2025-04-09T13:47:35.097741Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"## Some photos of \"0\": \"Cassava Bacterial Blight (CBB)\"","metadata":{}},{"cell_type":"code","source":"sample = train_labels[train_labels.label == 0].sample(3)\nplt.figure(figsize=(15, 5))\nfor ind, (image_id, label) in enumerate(zip(sample.image_id, sample.label)):\n    plt.subplot(1, 3, ind + 1)\n    img = cv2.imread(os.path.join(WORK_DIR, \"train_images\", image_id))\n    img = cv2.cvtColor(img, cv2.COLOR_BGR2RGB)\n    plt.imshow(img)\n    plt.axis(\"off\")\n    \nplt.show()","metadata":{"_kg_hide-input":true,"execution":{"iopub.status.busy":"2025-04-09T13:47:39.458407Z","iopub.execute_input":"2025-04-09T13:47:39.458693Z","iopub.status.idle":"2025-04-09T13:47:40.044915Z","shell.execute_reply.started":"2025-04-09T13:47:39.458669Z","shell.execute_reply":"2025-04-09T13:47:40.044152Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"## Some photos of \"1\": \"Cassava Brown Streak Disease (CBSD)\"","metadata":{}},{"cell_type":"code","source":"sample = train_labels[train_labels.label == 1].sample(3)\nplt.figure(figsize=(15, 5))\nfor ind, (image_id, label) in enumerate(zip(sample.image_id, sample.label)):\n    plt.subplot(1, 3, ind + 1)\n    img = cv2.imread(os.path.join(WORK_DIR, \"train_images\", image_id))\n    img = cv2.cvtColor(img, cv2.COLOR_BGR2RGB)\n    plt.imshow(img)\n    plt.axis(\"off\")\n    \nplt.show()","metadata":{"_kg_hide-input":true,"execution":{"iopub.status.busy":"2025-04-09T13:47:45.361596Z","iopub.execute_input":"2025-04-09T13:47:45.361901Z","iopub.status.idle":"2025-04-09T13:47:45.954672Z","shell.execute_reply.started":"2025-04-09T13:47:45.361873Z","shell.execute_reply":"2025-04-09T13:47:45.953746Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"## Some photos of \"2\": \"Cassava Green Mottle (CGM)\"","metadata":{}},{"cell_type":"code","source":"sample = train_labels[train_labels.label == 2].sample(3)\nplt.figure(figsize=(15, 5))\nfor ind, (image_id, label) in enumerate(zip(sample.image_id, sample.label)):\n    plt.subplot(1, 3, ind + 1)\n    img = cv2.imread(os.path.join(WORK_DIR, \"train_images\", image_id))\n    img = cv2.cvtColor(img, cv2.COLOR_BGR2RGB)\n    plt.imshow(img)\n    plt.axis(\"off\")\n    \nplt.show()","metadata":{"_kg_hide-input":true,"execution":{"iopub.status.busy":"2025-04-09T13:47:51.994789Z","iopub.execute_input":"2025-04-09T13:47:51.995095Z","iopub.status.idle":"2025-04-09T13:47:52.803929Z","shell.execute_reply.started":"2025-04-09T13:47:51.995068Z","shell.execute_reply":"2025-04-09T13:47:52.803170Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"## Some photos of \"3\": \"Cassava Mosaic Disease (CMD)\"","metadata":{}},{"cell_type":"code","source":"sample = train_labels[train_labels.label == 3].sample(3)\nplt.figure(figsize=(15, 5))\nfor ind, (image_id, label) in enumerate(zip(sample.image_id, sample.label)):\n    plt.subplot(1, 3, ind + 1)\n    img = cv2.imread(os.path.join(WORK_DIR, \"train_images\", image_id))\n    img = cv2.cvtColor(img, cv2.COLOR_BGR2RGB)\n    plt.imshow(img)\n    plt.axis(\"off\")\n    \nplt.show()","metadata":{"_kg_hide-input":true,"execution":{"iopub.status.busy":"2025-04-09T13:47:59.028598Z","iopub.execute_input":"2025-04-09T13:47:59.028944Z","iopub.status.idle":"2025-04-09T13:47:59.602153Z","shell.execute_reply.started":"2025-04-09T13:47:59.028903Z","shell.execute_reply":"2025-04-09T13:47:59.601313Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"## Some photos of \"4\": \"Healthy\"","metadata":{}},{"cell_type":"code","source":"sample = train_labels[train_labels.label == 4].sample(3)\nplt.figure(figsize=(15, 5))\nfor ind, (image_id, label) in enumerate(zip(sample.image_id, sample.label)):\n    plt.subplot(1, 3, ind + 1)\n    img = cv2.imread(os.path.join(WORK_DIR, \"train_images\", image_id))\n    img = cv2.cvtColor(img, cv2.COLOR_BGR2RGB)\n    plt.imshow(img)\n    plt.axis(\"off\")\n    \nplt.show()","metadata":{"_kg_hide-input":true,"execution":{"iopub.status.busy":"2025-04-09T13:48:05.913326Z","iopub.execute_input":"2025-04-09T13:48:05.913625Z","iopub.status.idle":"2025-04-09T13:48:06.495702Z","shell.execute_reply.started":"2025-04-09T13:48:05.913591Z","shell.execute_reply":"2025-04-09T13:48:06.494833Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"# Preparation for modeling","metadata":{}},{"cell_type":"code","source":"# Main parameters\nBATCH_SIZE = 8\nSTEPS_PER_EPOCH = len(train_labels)*0.8 / BATCH_SIZE\nVALIDATION_STEPS = len(train_labels)*0.2 / BATCH_SIZE\nEPOCHS = 3\nTARGET_SIZE = 512","metadata":{"execution":{"iopub.status.busy":"2025-04-09T13:48:16.962936Z","iopub.execute_input":"2025-04-09T13:48:16.963303Z","iopub.status.idle":"2025-04-09T13:48:16.967677Z","shell.execute_reply.started":"2025-04-09T13:48:16.963274Z","shell.execute_reply":"2025-04-09T13:48:16.966895Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"### ImageDataGenerator","metadata":{}},{"cell_type":"code","source":"train_labels.label = train_labels.label.astype('str')\n\ntrain_datagen = ImageDataGenerator(validation_split = 0.2,\n                                     preprocessing_function = None,\n                                     rotation_range = 45,\n                                     zoom_range = 0.2,\n                                     horizontal_flip = True,\n                                     vertical_flip = True,\n                                     fill_mode = 'nearest',\n                                     shear_range = 0.1,\n                                     height_shift_range = 0.1,\n                                     width_shift_range = 0.1)\n\ntrain_generator = train_datagen.flow_from_dataframe(train_labels,\n                         directory = os.path.join(WORK_DIR, \"train_images\"),\n                         subset = \"training\",\n                         x_col = \"image_id\",\n                         y_col = \"label\",\n                         target_size = (TARGET_SIZE, TARGET_SIZE),\n                         batch_size = BATCH_SIZE,\n                         class_mode = \"sparse\")\n\n\nvalidation_datagen = ImageDataGenerator(validation_split = 0.2)\n\nvalidation_generator = validation_datagen.flow_from_dataframe(train_labels,\n                         directory = os.path.join(WORK_DIR, \"train_images\"),\n                         subset = \"validation\",\n                         x_col = \"image_id\",\n                         y_col = \"label\",\n                         target_size = (TARGET_SIZE, TARGET_SIZE),\n                         batch_size = BATCH_SIZE,\n                         class_mode = \"sparse\")","metadata":{"execution":{"iopub.status.busy":"2025-04-09T13:48:24.408794Z","iopub.execute_input":"2025-04-09T13:48:24.409101Z","iopub.status.idle":"2025-04-09T13:48:44.755682Z","shell.execute_reply.started":"2025-04-09T13:48:24.409076Z","shell.execute_reply":"2025-04-09T13:48:44.754802Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"### Random photo before augmentation","metadata":{}},{"cell_type":"code","source":"img_path = os.path.join(WORK_DIR, \"train_images\", train_labels.image_id[20])\nimg = image.load_img(img_path, target_size = (TARGET_SIZE, TARGET_SIZE))\nimg_tensor = image.img_to_array(img)\nimg_tensor = np.expand_dims(img_tensor, axis = 0)\nimg_tensor /= 255.\n\nplt.imshow(img_tensor[0])\nplt.axis('off')\nplt.show()","metadata":{"execution":{"iopub.status.busy":"2025-04-09T13:49:16.190458Z","iopub.execute_input":"2025-04-09T13:49:16.190761Z","iopub.status.idle":"2025-04-09T13:49:16.387677Z","shell.execute_reply.started":"2025-04-09T13:49:16.190734Z","shell.execute_reply":"2025-04-09T13:49:16.386858Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"### This photo after augmentation","metadata":{}},{"cell_type":"code","source":"generator = train_datagen.flow_from_dataframe(train_labels.iloc[20:21],\n                         directory = os.path.join(WORK_DIR, \"train_images\"),\n                         x_col = \"image_id\",\n                         y_col = \"label\",\n                         target_size = (TARGET_SIZE, TARGET_SIZE),\n                         batch_size = BATCH_SIZE,\n                         class_mode = \"sparse\")\n\naug_images = [generator[0][0][0]/255 for i in range(10)]\nfig, axes = plt.subplots(2, 5, figsize = (20, 10))\naxes = axes.flatten()\nfor img, ax in zip(aug_images, axes):\n    ax.imshow(img)\n    ax.axis('off')\nplt.tight_layout()\nplt.show()","metadata":{"execution":{"iopub.status.busy":"2025-04-09T13:49:21.723130Z","iopub.execute_input":"2025-04-09T13:49:21.723422Z","iopub.status.idle":"2025-04-09T13:49:24.635235Z","shell.execute_reply.started":"2025-04-09T13:49:21.723397Z","shell.execute_reply":"2025-04-09T13:49:24.634045Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"# Modeling","metadata":{}},{"cell_type":"markdown","source":"### Creating the CNN","metadata":{}},{"cell_type":"code","source":"def create_model():\n    conv_base = EfficientNetB0(include_top = False, weights = None,\n                               input_shape = (TARGET_SIZE, TARGET_SIZE, 3))\n    model = conv_base.output\n    model = layers.GlobalAveragePooling2D()(model)\n    model = layers.Dense(5, activation = \"softmax\")(model)\n    model = models.Model(conv_base.input, model)\n\n    model.compile(optimizer = Adam(lr = 0.001),\n                  loss = \"sparse_categorical_crossentropy\",\n                  metrics = [\"acc\"])\n    return model","metadata":{"execution":{"iopub.status.busy":"2025-04-09T13:49:42.802436Z","iopub.execute_input":"2025-04-09T13:49:42.802758Z","iopub.status.idle":"2025-04-09T13:49:42.808467Z","shell.execute_reply.started":"2025-04-09T13:49:42.802725Z","shell.execute_reply":"2025-04-09T13:49:42.807609Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"code","source":"model = create_model()\nmodel.summary()","metadata":{"_kg_hide-output":true,"execution":{"iopub.status.busy":"2025-04-09T13:49:49.169913Z","iopub.execute_input":"2025-04-09T13:49:49.170259Z","iopub.status.idle":"2025-04-09T13:49:50.834007Z","shell.execute_reply.started":"2025-04-09T13:49:49.170230Z","shell.execute_reply":"2025-04-09T13:49:50.833153Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"code","source":"print('Our EfficientNet CNN has %d layers' %len(model.layers))","metadata":{"_kg_hide-input":true,"execution":{"iopub.status.busy":"2025-04-09T13:50:05.856693Z","iopub.execute_input":"2025-04-09T13:50:05.857027Z","iopub.status.idle":"2025-04-09T13:50:05.861239Z","shell.execute_reply.started":"2025-04-09T13:50:05.856991Z","shell.execute_reply":"2025-04-09T13:50:05.860584Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"### Loading weights","metadata":{}},{"cell_type":"code","source":"model.load_weights('../input/cassava-leaf-disease-models/basic_EfNetB0_imagenet_512.h5')","metadata":{"execution":{"iopub.status.busy":"2025-04-09T13:50:08.646828Z","iopub.execute_input":"2025-04-09T13:50:08.647132Z","iopub.status.idle":"2025-04-09T13:50:09.435714Z","shell.execute_reply.started":"2025-04-09T13:50:08.647107Z","shell.execute_reply":"2025-04-09T13:50:09.434811Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"### Training","metadata":{}},{"cell_type":"code","source":"model_save = ModelCheckpoint('./EffNetB0_512_8_best_weights.h5', \n                             save_best_only = True, \n                             save_weights_only = True,\n                             monitor = 'val_loss', \n                             mode = 'min', verbose = 1)\nearly_stop = EarlyStopping(monitor = 'val_loss', min_delta = 0.001, \n                           patience = 5, mode = 'min', verbose = 1,\n                           restore_best_weights = True)\nreduce_lr = ReduceLROnPlateau(monitor = 'val_loss', factor = 0.3, \n                              patience = 2, min_delta = 0.001, \n                              mode = 'min', verbose = 1)\n\n\nhistory = model.fit(\n    train_generator,\n    steps_per_epoch = STEPS_PER_EPOCH,\n    epochs = EPOCHS,\n    validation_data = validation_generator,\n    validation_steps = VALIDATION_STEPS,\n    callbacks = [model_save, early_stop, reduce_lr]\n)","metadata":{"_kg_hide-output":true,"execution":{"iopub.status.busy":"2025-04-09T13:50:12.616225Z","iopub.execute_input":"2025-04-09T13:50:12.616546Z","iopub.status.idle":"2025-04-09T14:46:17.692757Z","shell.execute_reply.started":"2025-04-09T13:50:12.616515Z","shell.execute_reply":"2025-04-09T14:46:17.692115Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"code","source":"acc = history.history['acc']\nval_acc = history.history['val_acc']\nloss = history.history['loss']\nval_loss = history.history['val_loss']\n\nepochs = range(1, len(acc) + 1)\n\nfig, (ax1, ax2) = plt.subplots(1, 2, figsize=(15, 5))\nsns.set_style(\"white\")\nplt.suptitle('Train history', size = 15)\n\nax1.plot(epochs, acc, \"bo\", label = \"Training acc\")\nax1.plot(epochs, val_acc, \"b\", label = \"Validation acc\")\nax1.set_title(\"Training and validation acc\")\nax1.legend()\n\nax2.plot(epochs, loss, \"bo\", label = \"Training loss\", color = 'red')\nax2.plot(epochs, val_loss, \"b\", label = \"Validation loss\", color = 'red')\nax2.set_title(\"Training and validation loss\")\nax2.legend()\n\nplt.show()","metadata":{"_kg_hide-input":true,"execution":{"iopub.status.busy":"2025-04-09T14:48:43.034133Z","iopub.execute_input":"2025-04-09T14:48:43.034466Z","iopub.status.idle":"2025-04-09T14:48:43.483725Z","shell.execute_reply.started":"2025-04-09T14:48:43.034434Z","shell.execute_reply":"2025-04-09T14:48:43.483044Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"code","source":"model.save('./EffNetB0_512_8.h5')","metadata":{"execution":{"iopub.status.busy":"2024-09-30T13:54:37.742922Z","iopub.status.idle":"2024-09-30T13:54:37.743724Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"# Visualization of CNN intermediate activations","metadata":{}},{"cell_type":"code","source":"def activation_layer_vis(img, activation_layer = 0, layers = 10):\n    layer_outputs = [layer.output for layer in model.layers[:layers]]\n    activation_model = models.Model(inputs = model.input, outputs = layer_outputs)\n    activations = activation_model.predict(img)\n    \n    rows = int(activations[activation_layer].shape[3] / 3)\n    cols = int(activations[activation_layer].shape[3] / rows)\n    fig, axes = plt.subplots(rows, cols, figsize = (15, 15 * cols))\n    axes = axes.flatten()\n    \n    for i, ax in zip(range(activations[activation_layer].shape[3]), axes):\n        ax.matshow(activations[activation_layer][0, :, :, i], cmap = 'viridis')\n        ax.axis('off')\n    plt.tight_layout()\n    plt.show()","metadata":{"execution":{"iopub.status.busy":"2025-04-09T14:48:57.172204Z","iopub.execute_input":"2025-04-09T14:48:57.172498Z","iopub.status.idle":"2025-04-09T14:48:57.179588Z","shell.execute_reply.started":"2025-04-09T14:48:57.172474Z","shell.execute_reply":"2025-04-09T14:48:57.178583Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"### Visualization of the first layer","metadata":{}},{"cell_type":"code","source":"activation_layer_vis(img_tensor, 0)","metadata":{"_kg_hide-input":true,"execution":{"iopub.status.busy":"2025-04-09T14:49:03.492862Z","iopub.execute_input":"2025-04-09T14:49:03.493222Z","iopub.status.idle":"2025-04-09T14:49:03.998661Z","shell.execute_reply.started":"2025-04-09T14:49:03.493193Z","shell.execute_reply":"2025-04-09T14:49:03.997643Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"code","source":"def all_activations_vis(img, layers = 10):\n    layer_outputs = [layer.output for layer in model.layers[:layers]]\n    activation_model = models.Model(inputs = model.input, outputs = layer_outputs)\n    activations = activation_model.predict(img)\n    \n    layer_names = []\n    for layer in model.layers[:layers]: \n        layer_names.append(layer.name) \n\n    images_per_row = 3\n    for layer_name, layer_activation in zip(layer_names, activations): \n        n_features = layer_activation.shape[-1] \n\n        size = layer_activation.shape[1] \n\n        n_cols = n_features // images_per_row \n        display_grid = np.zeros((size * n_cols, images_per_row * size)) \n\n        for col in range(n_cols): \n            for row in range(images_per_row): \n                channel_image = layer_activation[0, :, :, col * images_per_row + row] \n                channel_image -= channel_image.mean() \n                channel_image /= channel_image.std() \n                channel_image *= 64 \n                channel_image += 128 \n                channel_image = np.clip(channel_image, 0, 255).astype('uint8') \n                display_grid[col * size : (col + 1) * size, \n                             row * size : (row + 1) * size] = channel_image \n        scale = 1. / size \n        plt.figure(figsize=(scale * 5 * display_grid.shape[1], \n                            scale * 5 * display_grid.shape[0])) \n        plt.title(layer_name) \n        plt.grid(False)\n        plt.axis('off')\n        plt.imshow(display_grid, aspect = 'auto', cmap = 'viridis')","metadata":{"execution":{"iopub.status.busy":"2025-04-09T14:49:16.519824Z","iopub.execute_input":"2025-04-09T14:49:16.520152Z","iopub.status.idle":"2025-04-09T14:49:16.530275Z","shell.execute_reply.started":"2025-04-09T14:49:16.520122Z","shell.execute_reply":"2025-04-09T14:49:16.529325Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"### Visualization of the first 5 layers","metadata":{}},{"cell_type":"code","source":"all_activations_vis(img_tensor, 5)","metadata":{"_kg_hide-input":true,"execution":{"iopub.status.busy":"2025-04-09T14:49:21.781313Z","iopub.execute_input":"2025-04-09T14:49:21.781628Z","iopub.status.idle":"2025-04-09T14:49:24.708882Z","shell.execute_reply.started":"2025-04-09T14:49:21.781595Z","shell.execute_reply":"2025-04-09T14:49:24.707861Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"markdown","source":"Visualization of intermediate activations gives a rough step-by-step understanding of how CNN works.","metadata":{}},{"cell_type":"markdown","source":"# Prediction","metadata":{}},{"cell_type":"code","source":"ss = pd.read_csv(os.path.join(WORK_DIR, \"sample_submission.csv\"))\nss","metadata":{"execution":{"iopub.status.busy":"2025-04-09T14:49:48.816455Z","iopub.execute_input":"2025-04-09T14:49:48.816770Z","iopub.status.idle":"2025-04-09T14:49:48.833585Z","shell.execute_reply.started":"2025-04-09T14:49:48.816738Z","shell.execute_reply":"2025-04-09T14:49:48.832753Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"code","source":"preds = []\n\nfor image_id in ss.image_id:\n    image = Image.open(os.path.join(WORK_DIR,  \"test_images\", image_id))\n    image = image.resize((TARGET_SIZE, TARGET_SIZE))\n    image = np.expand_dims(image, axis = 0)\n    preds.append(np.argmax(model.predict(image)))\n\nss['label'] = preds\nss","metadata":{"execution":{"iopub.status.busy":"2025-04-09T14:50:00.194258Z","iopub.execute_input":"2025-04-09T14:50:00.194532Z","iopub.status.idle":"2025-04-09T14:50:01.424186Z","shell.execute_reply.started":"2025-04-09T14:50:00.194509Z","shell.execute_reply":"2025-04-09T14:50:01.423455Z"},"trusted":true},"outputs":[],"execution_count":null},{"cell_type":"code","source":"ss.to_csv('submission.csv', index = False)","metadata":{"execution":{"iopub.status.busy":"2025-04-09T14:50:08.174912Z","iopub.execute_input":"2025-04-09T14:50:08.175257Z","iopub.status.idle":"2025-04-09T14:50:08.575521Z","shell.execute_reply.started":"2025-04-09T14:50:08.175230Z","shell.execute_reply":"2025-04-09T14:50:08.574822Z"},"trusted":true},"outputs":[],"execution_count":null}]}