{"metadata":{"kernelspec":{"language":"python","display_name":"Python 3","name":"python3"},"language_info":{"pygments_lexer":"ipython3","nbconvert_exporter":"python","version":"3.6.4","file_extension":".py","codemirror_mode":{"name":"ipython","version":3},"name":"python","mimetype":"text/x-python"}},"nbformat_minor":4,"nbformat":4,"cells":[{"cell_type":"code","source":"# This Python 3 environment comes with many helpful analytics libraries installed\n# It is defined by the kaggle/python Docker image: https://github.com/kaggle/docker-python\n# For example, here's several helpful packages to load\n\nimport numpy as np # linear algebra\nimport pandas as pd # data processing, CSV file I/O (e.g. pd.read_csv)\n\n# Input data files are available in the read-only \"../input/\" directory\n# For example, running this (by clicking run or pressing Shift+Enter) will list all files under the input directory\n\nimport os\nfor dirname, _, filenames in os.walk('../input/concrete-crack-images'):\n    for filename in filenames:\n        print(os.path.join(dirname, filename))\n\n# You can write up to 20GB to the current directory (/kaggle/working/) that gets preserved as output when you create a version using \"Save & Run All\" \n# You can also write temporary files to /kaggle/temp/, but they won't be saved outside of the current session","metadata":{"_uuid":"8f2839f25d086af736a60e9eeb907d3b93b6e0e5","_cell_guid":"b1076dfc-b9ad-4769-8c92-a6c4dae69d19","execution":{"iopub.status.busy":"2022-11-03T10:38:01.998308Z","iopub.execute_input":"2022-11-03T10:38:01.998661Z","iopub.status.idle":"2022-11-03T10:38:02.029069Z","shell.execute_reply.started":"2022-11-03T10:38:01.998629Z","shell.execute_reply":"2022-11-03T10:38:02.028130Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# U-Net Original Research Paper\n### Paper Name: U-Net: Convolutional Networks for Biomedical Image Segmentation\n### Paper Link: https://arxiv.org/pdf/1505.04597.pdf","metadata":{}},{"cell_type":"markdown","source":"# Configs","metadata":{}},{"cell_type":"code","source":"class ROOTDIR:\n    train = \"../input/concrete-crack-images/images\"\n    train_mask = \"../input/concrete-crack-images/masks\"","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:02.030665Z","iopub.execute_input":"2022-11-03T10:38:02.031495Z","iopub.status.idle":"2022-11-03T10:38:02.036629Z","shell.execute_reply.started":"2022-11-03T10:38:02.031461Z","shell.execute_reply":"2022-11-03T10:38:02.035413Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# General Imports","metadata":{}},{"cell_type":"code","source":"import os\nimport cv2\nimport zipfile\nimport numpy as np\nimport pandas as pd\nfrom PIL import Image\nfrom tqdm.auto import tqdm\nimport matplotlib.pyplot as plt","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:02.038523Z","iopub.execute_input":"2022-11-03T10:38:02.038862Z","iopub.status.idle":"2022-11-03T10:38:02.052720Z","shell.execute_reply.started":"2022-11-03T10:38:02.038829Z","shell.execute_reply":"2022-11-03T10:38:02.051870Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# Extracting Files","metadata":{}},{"cell_type":"markdown","source":"dirs = [\"../input/carvana-image-masking-challenge/train.zip\",\n        \"../input/carvana-image-masking-challenge/train_masks.zip\",\n        \"../input/carvana-image-masking-challenge/metadata.csv.zip\"]\n\nfor i in tqdm(dirs):\n    with zipfile.ZipFile(i) as z:\n        z.extractall()\n    ","metadata":{"execution":{"iopub.status.busy":"2022-11-02T09:12:58.429001Z","iopub.execute_input":"2022-11-02T09:12:58.429390Z","iopub.status.idle":"2022-11-02T09:13:09.189807Z","shell.execute_reply.started":"2022-11-02T09:12:58.429352Z","shell.execute_reply":"2022-11-02T09:13:09.188835Z"}}},{"cell_type":"markdown","source":"# Working with Train Images and Masks","metadata":{}},{"cell_type":"code","source":"train_img_lst = os.listdir('../input/concrete-crack-images/images') # \"./train\"\ntrain_mask_lst = os.listdir('../input/concrete-crack-images/masks') # \"./train_masks\"","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:45:35.537471Z","iopub.execute_input":"2022-11-03T10:45:35.537828Z","iopub.status.idle":"2022-11-03T10:45:35.553209Z","shell.execute_reply.started":"2022-11-03T10:45:35.537797Z","shell.execute_reply":"2022-11-03T10:45:35.552335Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"print(train_mask_lst[:5])\nprint(train_img_lst[:5])","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:45:38.645272Z","iopub.execute_input":"2022-11-03T10:45:38.645637Z","iopub.status.idle":"2022-11-03T10:45:38.651921Z","shell.execute_reply.started":"2022-11-03T10:45:38.645603Z","shell.execute_reply":"2022-11-03T10:45:38.650950Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"print(len(train_mask_lst))\nprint(len(train_img_lst))","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:45:41.265105Z","iopub.execute_input":"2022-11-03T10:45:41.265681Z","iopub.status.idle":"2022-11-03T10:45:41.270562Z","shell.execute_reply.started":"2022-11-03T10:45:41.265646Z","shell.execute_reply":"2022-11-03T10:45:41.269577Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"#### Sorting to make sure we get right image and right mask","metadata":{}},{"cell_type":"code","source":"sorted_train_mask_lst = sorted(train_mask_lst)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:45:47.464196Z","iopub.execute_input":"2022-11-03T10:45:47.464862Z","iopub.status.idle":"2022-11-03T10:45:47.470133Z","shell.execute_reply.started":"2022-11-03T10:45:47.464816Z","shell.execute_reply":"2022-11-03T10:45:47.469108Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"sorted_train_img_lst = sorted(train_img_lst)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:45:49.055589Z","iopub.execute_input":"2022-11-03T10:45:49.055945Z","iopub.status.idle":"2022-11-03T10:45:49.060521Z","shell.execute_reply.started":"2022-11-03T10:45:49.055906Z","shell.execute_reply":"2022-11-03T10:45:49.059630Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"print(sorted_train_mask_lst[:16])\nprint(sorted_train_img_lst[:16])","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:45:51.253342Z","iopub.execute_input":"2022-11-03T10:45:51.253709Z","iopub.status.idle":"2022-11-03T10:45:51.259372Z","shell.execute_reply.started":"2022-11-03T10:45:51.253677Z","shell.execute_reply":"2022-11-03T10:45:51.257541Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# Visualizing Images with their Mask\n### Making sure images and mask are paired correctly.","metadata":{}},{"cell_type":"code","source":"def show_images(imgs_lst,masks_lst,loops=2):\n    for i in range(loops):\n        img_path = os.path.join(ROOTDIR.train,imgs_lst[i])\n        mask_path = os.path.join(ROOTDIR.train_mask,masks_lst[i])\n        img = Image.open(img_path)\n        mask = Image.open(mask_path)\n        print(img_path)\n        print(img.size)\n        print(type(img))\n        plt.imshow(img)\n        plt.show()\n        print(mask_path)\n        print(mask.size)\n        plt.imshow(mask)\n        plt.show()\n        print(\"----------------------------------------------------\")\n\nshow_images(sorted_train_img_lst, sorted_train_mask_lst)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:45:53.083173Z","iopub.execute_input":"2022-11-03T10:45:53.083539Z","iopub.status.idle":"2022-11-03T10:45:53.840056Z","shell.execute_reply.started":"2022-11-03T10:45:53.083506Z","shell.execute_reply":"2022-11-03T10:45:53.839162Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# PyTorch Imports","metadata":{}},{"cell_type":"code","source":"pwd","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:02.940582Z","iopub.execute_input":"2022-11-03T10:38:02.941245Z","iopub.status.idle":"2022-11-03T10:38:02.948256Z","shell.execute_reply.started":"2022-11-03T10:38:02.941208Z","shell.execute_reply":"2022-11-03T10:38:02.947039Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"import torch\nimport torchvision\nimport torch.nn as nn\nimport albumentations as A\nimport torch.optim as optim\nfrom torchvision import models\nimport torch.nn.functional as F\nfrom torch.optim import lr_scheduler\nimport torchvision.datasets as datasets\nimport torchvision.transforms as transforms\nfrom torchvision.datasets import ImageFolder\nfrom albumentations.pytorch import ToTensorV2 \nfrom torch.utils.data import DataLoader, Dataset","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:02.949807Z","iopub.execute_input":"2022-11-03T10:38:02.950205Z","iopub.status.idle":"2022-11-03T10:38:02.958431Z","shell.execute_reply.started":"2022-11-03T10:38:02.950170Z","shell.execute_reply":"2022-11-03T10:38:02.957453Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# PyTorch Configuration","metadata":{}},{"cell_type":"code","source":"class CFG:\n    device = torch.device(\"cuda\" if torch.cuda.is_available() else \"cpu\")\n    split_pct = 0.2\n    learning_rate = 3e-4\n    batch_size = 4\n    epochs = 3","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:02.961669Z","iopub.execute_input":"2022-11-03T10:38:02.962622Z","iopub.status.idle":"2022-11-03T10:38:02.967784Z","shell.execute_reply.started":"2022-11-03T10:38:02.962587Z","shell.execute_reply":"2022-11-03T10:38:02.967141Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"seed = 123\nnp.random.seed(seed)\ntorch.manual_seed(seed)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:02.969065Z","iopub.execute_input":"2022-11-03T10:38:02.969820Z","iopub.status.idle":"2022-11-03T10:38:02.980942Z","shell.execute_reply.started":"2022-11-03T10:38:02.969782Z","shell.execute_reply":"2022-11-03T10:38:02.980226Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"CFG.device","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:02.982112Z","iopub.execute_input":"2022-11-03T10:38:02.983137Z","iopub.status.idle":"2022-11-03T10:38:02.991689Z","shell.execute_reply.started":"2022-11-03T10:38:02.983103Z","shell.execute_reply":"2022-11-03T10:38:02.991000Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# Working with data","metadata":{}},{"cell_type":"markdown","source":"### Shuffling the data.","metadata":{}},{"cell_type":"code","source":"#permuted_train_img_lst = np.random.permutation(np.array(sorted_train_img_lst))\n#permuted_train_mask_lst = [x.replace(\".jpg\", \"_mask.gif\") for x in permuted_train_img_lst]\npermuted_train_img_lst=sorted_train_img_lst\npermuted_train_mask_lst=sorted_train_mask_lst\nprint(permuted_train_img_lst[:5])\nprint(permuted_train_mask_lst[:5])","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:02.992840Z","iopub.execute_input":"2022-11-03T10:38:02.993825Z","iopub.status.idle":"2022-11-03T10:38:03.003544Z","shell.execute_reply.started":"2022-11-03T10:38:02.993792Z","shell.execute_reply":"2022-11-03T10:38:03.002502Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"show_images(permuted_train_img_lst,permuted_train_mask_lst)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:03.004672Z","iopub.execute_input":"2022-11-03T10:38:03.005502Z","iopub.status.idle":"2022-11-03T10:38:03.767215Z","shell.execute_reply.started":"2022-11-03T10:38:03.005466Z","shell.execute_reply":"2022-11-03T10:38:03.766231Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"### Splitting into Training and Validation","metadata":{}},{"cell_type":"code","source":"length = len(permuted_train_img_lst)\nprint(length*0.2) # convert this to int","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:03.768808Z","iopub.execute_input":"2022-11-03T10:38:03.769461Z","iopub.status.idle":"2022-11-03T10:38:03.774866Z","shell.execute_reply.started":"2022-11-03T10:38:03.769424Z","shell.execute_reply":"2022-11-03T10:38:03.773965Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"train_images_list = permuted_train_img_lst[int(CFG.split_pct*len(permuted_train_img_lst)) :]\ntrain_masks_list = permuted_train_mask_lst[int(CFG.split_pct*len(permuted_train_mask_lst)) :]\nprint(len(train_masks_list))\n\nval_images_list = permuted_train_img_lst[: int(CFG.split_pct*len(permuted_train_img_lst))]\nval_masks_list = permuted_train_mask_lst[: int(CFG.split_pct*len(permuted_train_mask_lst))]\nprint(len(val_masks_list))\n\n# 4071+1017=5088 (split includes all items)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:03.776132Z","iopub.execute_input":"2022-11-03T10:38:03.777066Z","iopub.status.idle":"2022-11-03T10:38:03.788466Z","shell.execute_reply.started":"2022-11-03T10:38:03.777030Z","shell.execute_reply":"2022-11-03T10:38:03.787542Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"### Visualizing Train Dataset","metadata":{}},{"cell_type":"code","source":"show_images(train_images_list,train_masks_list)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:03.789739Z","iopub.execute_input":"2022-11-03T10:38:03.790022Z","iopub.status.idle":"2022-11-03T10:38:04.540983Z","shell.execute_reply.started":"2022-11-03T10:38:03.789989Z","shell.execute_reply":"2022-11-03T10:38:04.540020Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"### Visualizing Validation Dataset","metadata":{}},{"cell_type":"code","source":"show_images(val_images_list,val_masks_list)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:04.542567Z","iopub.execute_input":"2022-11-03T10:38:04.543253Z","iopub.status.idle":"2022-11-03T10:38:05.310000Z","shell.execute_reply.started":"2022-11-03T10:38:04.543214Z","shell.execute_reply":"2022-11-03T10:38:05.309104Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# Dataset Class","metadata":{}},{"cell_type":"code","source":"class CarvanaDataset(Dataset):\n    def __init__(self,img_list,mask_list,transform=None):\n        self.img_list = img_list\n        self.mask_list = mask_list\n        self.transform = transform\n        \n    def __len__(self):\n        return len(self.img_list)\n    \n    def __getitem__(self,index):\n        img_path = os.path.join(ROOTDIR.train,self.img_list[index])\n        mask_path = os.path.join(ROOTDIR.train_mask,self.mask_list[index])\n        img = Image.open(img_path)\n        mask = Image.open(mask_path)\n        img = np.array(img)\n        mask = np.array(mask)\n        mask[mask==255.0] = 1.0\n        #img_mask_dict = {\"image\": img, \"mask\": mask}\n        \n        if self.transform:\n            augmentation = self.transform(image=img, mask=mask)\n            img = augmentation[\"image\"]\n            mask = augmentation[\"mask\"]\n            mask = torch.unsqueeze(mask,0)\n            #transformations = self.transform(image=img, mask=mask)\n            #img = transformations[\"image\"]\n            #mask = transformations[\"mask\"]\n            \n        return img,mask","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:46:11.197518Z","iopub.execute_input":"2022-11-03T10:46:11.197888Z","iopub.status.idle":"2022-11-03T10:46:11.206410Z","shell.execute_reply.started":"2022-11-03T10:46:11.197855Z","shell.execute_reply":"2022-11-03T10:46:11.205134Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"train_transform = A.Compose([A.Resize(572,572), \n                             A.Rotate(limit=15,p=0.1),\n                             A.HorizontalFlip(p=0.5),\n                             A.Normalize(mean=(0,0,0),std=(1,1,1),max_pixel_value=255),\n                             ToTensorV2()])\n\nval_transform = A.Compose([A.Resize(572,572),\n                           A.Normalize(mean=(0,0,0),std=(1,1,1),max_pixel_value=255),\n                           ToTensorV2()])","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:05.325610Z","iopub.execute_input":"2022-11-03T10:38:05.326052Z","iopub.status.idle":"2022-11-03T10:38:05.339601Z","shell.execute_reply.started":"2022-11-03T10:38:05.326015Z","shell.execute_reply":"2022-11-03T10:38:05.338651Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"train_dataset = CarvanaDataset(train_images_list, train_masks_list, transform = train_transform)\nval_dataset = CarvanaDataset(val_images_list, val_masks_list, transform = train_transform)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:05.341034Z","iopub.execute_input":"2022-11-03T10:38:05.341389Z","iopub.status.idle":"2022-11-03T10:38:05.350609Z","shell.execute_reply.started":"2022-11-03T10:38:05.341356Z","shell.execute_reply":"2022-11-03T10:38:05.349747Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"idx = 200\nimg,mask = train_dataset[idx]","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:05.352277Z","iopub.execute_input":"2022-11-03T10:38:05.352625Z","iopub.status.idle":"2022-11-03T10:38:05.373952Z","shell.execute_reply.started":"2022-11-03T10:38:05.352592Z","shell.execute_reply":"2022-11-03T10:38:05.373080Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"mask.shape","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:05.375809Z","iopub.execute_input":"2022-11-03T10:38:05.376439Z","iopub.status.idle":"2022-11-03T10:38:05.382600Z","shell.execute_reply.started":"2022-11-03T10:38:05.376405Z","shell.execute_reply":"2022-11-03T10:38:05.381678Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"img.max()","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:05.384441Z","iopub.execute_input":"2022-11-03T10:38:05.384804Z","iopub.status.idle":"2022-11-03T10:38:05.394748Z","shell.execute_reply.started":"2022-11-03T10:38:05.384771Z","shell.execute_reply":"2022-11-03T10:38:05.393548Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"def show_single_img(img,mask,index=None,train=True):\n    if index:\n        if train:\n            img,mask = train_dataset[index]\n        else:\n            img,mask = val_dataset[index]\n    plt.imshow(img.permute(1,2,0),cmap=\"gray\")  # Convert (3, 572, 572) -> (572, 572, 3)\n    plt.show()\n    plt.imshow(mask.permute(1,2,0), cmap=\"gray\")  # Convert (1, 572, 572) -> (572, 572, 1)\n    print(mask.shape)\n    plt.show()\n    ","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:05.396574Z","iopub.execute_input":"2022-11-03T10:38:05.397068Z","iopub.status.idle":"2022-11-03T10:38:05.404563Z","shell.execute_reply.started":"2022-11-03T10:38:05.396990Z","shell.execute_reply":"2022-11-03T10:38:05.403731Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"print(\"---------------Train---------------\")\nshow_single_img(img,mask,index=15,train=False)\nprint(\"---------------Validation---------------\")\nshow_single_img(img,mask,index=15,train=True)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:05.405684Z","iopub.execute_input":"2022-11-03T10:38:05.406530Z","iopub.status.idle":"2022-11-03T10:38:06.457844Z","shell.execute_reply.started":"2022-11-03T10:38:05.406471Z","shell.execute_reply":"2022-11-03T10:38:06.456902Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# Dataloader","metadata":{}},{"cell_type":"code","source":"train_dataloader = DataLoader(train_dataset,batch_size=CFG.batch_size,shuffle=True)\nval_dataloader = DataLoader(val_dataset,batch_size=CFG.batch_size,shuffle=False)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:06.459178Z","iopub.execute_input":"2022-11-03T10:38:06.459923Z","iopub.status.idle":"2022-11-03T10:38:06.466486Z","shell.execute_reply.started":"2022-11-03T10:38:06.459887Z","shell.execute_reply":"2022-11-03T10:38:06.465018Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"a = iter(train_dataloader)\nimg,mask = a.next()\nprint(img.shape,mask.shape)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:38:06.468479Z","iopub.execute_input":"2022-11-03T10:38:06.468834Z","iopub.status.idle":"2022-11-03T10:38:06.519834Z","shell.execute_reply.started":"2022-11-03T10:38:06.468800Z","shell.execute_reply":"2022-11-03T10:38:06.518875Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"### Utility Functions","metadata":{}},{"cell_type":"code","source":"def double_conv(in_ch, out_ch):\n    conv = nn.Sequential(\n        nn.Conv2d(in_channels=in_ch,out_channels=out_ch,kernel_size=3,stride=1,padding=1),\n        nn.BatchNorm2d(out_ch),                                                            \n        nn.ReLU(inplace=True),\n        nn.Conv2d(in_channels=out_ch,out_channels=out_ch,kernel_size=3,stride=1,padding=1), \n        nn.BatchNorm2d(out_ch),                                                            \n        nn.ReLU(inplace=True)\n    )\n    \n    return conv\n\n#def cropper(og_tensor, target_tensor):\n#    og_shape = og_tensor.shape[2]\n#    target_shape = target_tensor.shape[2]\n#    delta = (og_shape - target_shape) // 2\n#    cropped_og_tensor = og_tensor[:,:,delta:og_shape-delta,delta:og_shape-delta]\n#    return cropped_og_tensor\n \n    \ndef padder(left_tensor, right_tensor): \n    # left_tensor is the tensor on the encoder side of UNET\n    # right_tensor is the tensor on the decoder side  of the UNET\n    \n    if left_tensor.shape != right_tensor.shape:\n        padded = torch.zeros(left_tensor.shape)\n        padded[:, :, :right_tensor.shape[2], :right_tensor.shape[3]] = right_tensor\n        return padded.to(CFG.device)\n    \n    return right_tensor.to(CFG.device)\n    ","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:46:19.610022Z","iopub.execute_input":"2022-11-03T10:46:19.610414Z","iopub.status.idle":"2022-11-03T10:46:19.619173Z","shell.execute_reply.started":"2022-11-03T10:46:19.610378Z","shell.execute_reply":"2022-11-03T10:46:19.618018Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# UNET MODEL FROM SCRATCH","metadata":{}},{"cell_type":"code","source":"class UNET(nn.Module):\n    def __init__(self,in_chnls, n_classes):\n        super(UNET,self).__init__()\n        \n        self.in_chnls = in_chnls\n        self.n_classes = n_classes\n        \n        self.max_pool = nn.MaxPool2d(kernel_size=2,stride=2)\n        \n        self.down_conv_1 = double_conv(in_ch=self.in_chnls,out_ch=64)\n        self.down_conv_2 = double_conv(in_ch=64,out_ch=128)\n        self.down_conv_3 = double_conv(in_ch=128,out_ch=256)\n        self.down_conv_4 = double_conv(in_ch=256,out_ch=512)\n        self.down_conv_5 = double_conv(in_ch=512,out_ch=1024)\n        #print(self.down_conv_1)\n        \n        self.up_conv_trans_1 = nn.ConvTranspose2d(in_channels=1024,out_channels=512,kernel_size=2,stride=2)\n        self.up_conv_trans_2 = nn.ConvTranspose2d(in_channels=512,out_channels=256,kernel_size=2,stride=2)\n        self.up_conv_trans_3 = nn.ConvTranspose2d(in_channels=256,out_channels=128,kernel_size=2,stride=2)\n        self.up_conv_trans_4 = nn.ConvTranspose2d(in_channels=128,out_channels=64,kernel_size=2,stride=2)\n        \n        self.up_conv_1 = double_conv(in_ch=1024,out_ch=512)\n        self.up_conv_2 = double_conv(in_ch=512,out_ch=256)\n        self.up_conv_3 = double_conv(in_ch=256,out_ch=128)\n        self.up_conv_4 = double_conv(in_ch=128,out_ch=64)\n        \n        self.conv_1x1 = nn.Conv2d(in_channels=64,out_channels=self.n_classes,kernel_size=1,stride=1)\n        \n    def forward(self,x):\n        \n        # encoding\n        x1 = self.down_conv_1(x)\n        #print(\"X1\", x1.shape)\n        p1 = self.max_pool(x1)\n        #print(\"p1\", p1.shape)\n        x2 = self.down_conv_2(p1)\n        #print(\"X2\", x2.shape)\n        p2 = self.max_pool(x2)\n        #print(\"p2\", p2.shape)\n        x3 = self.down_conv_3(p2)\n        #print(\"X2\", x3.shape)\n        p3 = self.max_pool(x3)\n        #print(\"p3\", p3.shape)\n        x4 = self.down_conv_4(p3)\n        #print(\"X4\", x4.shape)\n        p4 = self.max_pool(x4)\n        #print(\"p4\", p4.shape)\n        x5 = self.down_conv_5(p4)\n        #print(\"X5\", x5.shape)\n        \n        # decoding\n        d1 = self.up_conv_trans_1(x5)  # up transpose convolution (\"up sampling\" as called in UNET paper)\n        pad1 = padder(x4,d1) # padding d1 to match x4 shape\n        cat1 = torch.cat([x4,pad1],dim=1) # concatenating padded d1 and x4 on channel dimension(dim 1) [batch(dim 0),channel(dim 1),height(dim 2),width(dim 3)]\n        uc1 = self.up_conv_1(cat1) # 1st up double convolution\n        \n        d2 = self.up_conv_trans_2(uc1)\n        pad2 = padder(x3,d2)\n        cat2 = torch.cat([x3,pad2],dim=1)\n        uc2 = self.up_conv_2(cat2)\n        \n        d3 = self.up_conv_trans_3(uc2)\n        pad3 = padder(x2,d3)\n        cat3 = torch.cat([x2,pad3],dim=1)\n        uc3 = self.up_conv_3(cat3)\n        \n        d4 = self.up_conv_trans_4(uc3)\n        pad4 = padder(x1,d4)\n        cat4 = torch.cat([x1,pad4],dim=1)\n        uc4 = self.up_conv_4(cat4)\n        \n        conv_1x1 = self.conv_1x1(uc4)\n        return conv_1x1\n        #print(conv_1x1.shape)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:46:22.483534Z","iopub.execute_input":"2022-11-03T10:46:22.483885Z","iopub.status.idle":"2022-11-03T10:46:22.499868Z","shell.execute_reply.started":"2022-11-03T10:46:22.483853Z","shell.execute_reply":"2022-11-03T10:46:22.498991Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# Training and Validation","metadata":{}},{"cell_type":"markdown","source":"### Train Function\n","metadata":{}},{"cell_type":"code","source":"def train_model(model,dataloader,criterion,optimizer):\n    model.train()\n    train_running_loss = 0.0\n    for j,img_mask in enumerate(tqdm(dataloader)):\n        img = img_mask[0].float().to(CFG.device)\n        #print(\" ----- IMAGE -----\")\n        #print(img)\n        mask = img_mask[1].float().to(CFG.device)\n        #print(\" ----- MASK -----\")\n        #print(mask)\n        \n        y_pred = model(img)\n        #print(\" ----- Y PRED -----\")\n        #print(y_pred)\n        #print(\" ----- Y PRED SHAPE -----\")#\n        #print(y_pred.shape)\n        optimizer.zero_grad()\n        \n        loss = criterion(y_pred,mask)\n        \n        train_running_loss += loss.item() * CFG.batch_size\n        \n        loss.backward()\n        optimizer.step()\n        \n    train_loss = train_running_loss / (j+1)\n    return train_loss","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:46:27.792581Z","iopub.execute_input":"2022-11-03T10:46:27.792943Z","iopub.status.idle":"2022-11-03T10:46:27.800336Z","shell.execute_reply.started":"2022-11-03T10:46:27.792911Z","shell.execute_reply":"2022-11-03T10:46:27.799004Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"### Validation Function","metadata":{}},{"cell_type":"code","source":"def val_model(model,dataloader,criterion,optimizer):\n    model.eval()\n    val_running_loss = 0\n    with torch.no_grad():\n        for j,img_mask in enumerate(tqdm(dataloader)):\n            img = img_mask[0].float().to(CFG.device)\n            mask = img_mask[1].float().to(CFG.device)\n            y_pred = model(img)\n            loss = criterion(y_pred,mask)\n            \n            val_running_loss += loss.item() * CFG.batch_size\n            \n        val_loss = val_running_loss / (j+1)\n    return val_loss","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:46:31.230148Z","iopub.execute_input":"2022-11-03T10:46:31.230517Z","iopub.status.idle":"2022-11-03T10:46:31.237877Z","shell.execute_reply.started":"2022-11-03T10:46:31.230484Z","shell.execute_reply":"2022-11-03T10:46:31.236610Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"model = UNET(in_chnls = 3, n_classes = 1).to(CFG.device)\noptimizer = optim.Adam(model.parameters(), lr = CFG.learning_rate)\ncriterion = nn.BCEWithLogitsLoss()\ntrain_loss_lst = []\nval_loss_lst = []  ","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:46:33.770002Z","iopub.execute_input":"2022-11-03T10:46:33.771488Z","iopub.status.idle":"2022-11-03T10:46:34.065225Z","shell.execute_reply.started":"2022-11-03T10:46:33.771437Z","shell.execute_reply":"2022-11-03T10:46:34.064289Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"### Train and Validation Loop","metadata":{}},{"cell_type":"code","source":"for i in tqdm(range(CFG.epochs)):\n    train_loss = train_model(model=model,dataloader=train_dataloader,criterion=criterion,optimizer=optimizer)\n    val_loss = val_model(model=model,dataloader=val_dataloader,criterion=criterion,optimizer=optimizer)\n    train_loss_lst.append(train_loss)\n    val_loss_lst.append(val_loss)\n    print(f\" Train Loss : {train_loss:.4f}\")\n    print(f\" Validation Loss : {val_loss:.4f}\")","metadata":{"execution":{"iopub.status.busy":"2022-11-03T10:52:20.926997Z","iopub.execute_input":"2022-11-03T10:52:20.927662Z","iopub.status.idle":"2022-11-03T11:10:59.159947Z","shell.execute_reply.started":"2022-11-03T10:52:20.927627Z","shell.execute_reply":"2022-11-03T11:10:59.158996Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"### Training and Validation Loss Plot","metadata":{}},{"cell_type":"code","source":"plt.plot(train_loss_lst, color=\"green\", label='train loss')\nplt.plot(val_loss_lst, color=\"red\", label='validation loss')\nplt.xlabel(\"epochs\")\nplt.ylabel(\"loss\")\nplt.legend()\nplt.show()","metadata":{"execution":{"iopub.status.busy":"2022-11-03T11:15:49.684471Z","iopub.execute_input":"2022-11-03T11:15:49.684832Z","iopub.status.idle":"2022-11-03T11:15:49.893846Z","shell.execute_reply.started":"2022-11-03T11:15:49.684800Z","shell.execute_reply":"2022-11-03T11:15:49.892934Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# Saving Model","metadata":{}},{"cell_type":"code","source":"TRAINED_FILE = \"./unet_scratch.pth\"","metadata":{"execution":{"iopub.status.busy":"2022-11-03T11:16:20.205260Z","iopub.execute_input":"2022-11-03T11:16:20.205617Z","iopub.status.idle":"2022-11-03T11:16:20.212085Z","shell.execute_reply.started":"2022-11-03T11:16:20.205586Z","shell.execute_reply":"2022-11-03T11:16:20.210812Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"torch.save(model.state_dict(), TRAINED_FILE)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T11:16:21.541993Z","iopub.execute_input":"2022-11-03T11:16:21.542680Z","iopub.status.idle":"2022-11-03T11:16:21.837208Z","shell.execute_reply.started":"2022-11-03T11:16:21.542643Z","shell.execute_reply":"2022-11-03T11:16:21.835657Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"from IPython.display import FileLink\nFileLink(TRAINED_FILE)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T11:16:22.974065Z","iopub.execute_input":"2022-11-03T11:16:22.974684Z","iopub.status.idle":"2022-11-03T11:16:22.981207Z","shell.execute_reply.started":"2022-11-03T11:16:22.974648Z","shell.execute_reply":"2022-11-03T11:16:22.980217Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# Testing","metadata":{}},{"cell_type":"code","source":"trained_model = UNET(in_chnls = 3, n_classes = 1)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T11:16:25.519428Z","iopub.execute_input":"2022-11-03T11:16:25.519846Z","iopub.status.idle":"2022-11-03T11:16:25.774336Z","shell.execute_reply.started":"2022-11-03T11:16:25.519805Z","shell.execute_reply":"2022-11-03T11:16:25.773386Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"UNET_TRAINED = \"./unet_scratch.pth\"","metadata":{"execution":{"iopub.status.busy":"2022-11-03T11:16:48.403200Z","iopub.execute_input":"2022-11-03T11:16:48.404242Z","iopub.status.idle":"2022-11-03T11:16:48.408773Z","shell.execute_reply.started":"2022-11-03T11:16:48.404198Z","shell.execute_reply":"2022-11-03T11:16:48.407820Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"trained_model.load_state_dict(torch.load(UNET_TRAINED))","metadata":{"execution":{"iopub.status.busy":"2022-11-03T11:16:50.287456Z","iopub.execute_input":"2022-11-03T11:16:50.287819Z","iopub.status.idle":"2022-11-03T11:16:50.406814Z","shell.execute_reply.started":"2022-11-03T11:16:50.287787Z","shell.execute_reply":"2022-11-03T11:16:50.405431Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"trained_model = trained_model.to(\"cuda\")\ntrained_model.eval()\n","metadata":{"execution":{"iopub.status.busy":"2022-11-03T11:17:09.611101Z","iopub.execute_input":"2022-11-03T11:17:09.611450Z","iopub.status.idle":"2022-11-03T11:17:09.621748Z","shell.execute_reply.started":"2022-11-03T11:17:09.611419Z","shell.execute_reply":"2022-11-03T11:17:09.620652Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"img_path = \"../input/carvana-image-masking-challenge/29bb3ece3180_11.jpg\"\n\nimg = cv2.imread(img_path)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T11:19:24.390913Z","iopub.execute_input":"2022-11-03T11:19:24.391643Z","iopub.status.idle":"2022-11-03T11:19:24.491307Z","shell.execute_reply.started":"2022-11-03T11:19:24.391581Z","shell.execute_reply":"2022-11-03T11:19:24.490029Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"plt.imshow(img)\nplt.show()","metadata":{"execution":{"iopub.status.busy":"2022-11-03T11:19:26.768800Z","iopub.execute_input":"2022-11-03T11:19:26.769177Z","iopub.status.idle":"2022-11-03T11:19:27.217778Z","shell.execute_reply.started":"2022-11-03T11:19:26.769144Z","shell.execute_reply":"2022-11-03T11:19:27.216863Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"test_transform = A.Compose([A.Resize(572,572),\n                           A.Normalize(mean=(0,0,0),std=(1,1,1),max_pixel_value=255),\n                           ToTensorV2()])","metadata":{"execution":{"iopub.status.busy":"2022-11-03T11:19:29.981583Z","iopub.execute_input":"2022-11-03T11:19:29.982057Z","iopub.status.idle":"2022-11-03T11:19:29.989977Z","shell.execute_reply.started":"2022-11-03T11:19:29.982009Z","shell.execute_reply":"2022-11-03T11:19:29.988930Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"test_image = test_transform(image = img)\n\nprint(test_image)\n\nprint(test_image[\"image\"].dtype)\nprint(test_image[\"image\"].shape)\n\nimg = test_image[\"image\"].unsqueeze(0)\nprint(img.shape)\n\nimg = img.to(\"cuda\")","metadata":{"execution":{"iopub.status.busy":"2022-11-03T11:19:31.586082Z","iopub.execute_input":"2022-11-03T11:19:31.586997Z","iopub.status.idle":"2022-11-03T11:19:31.604742Z","shell.execute_reply.started":"2022-11-03T11:19:31.586913Z","shell.execute_reply":"2022-11-03T11:19:31.603493Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"pred = trained_model(img)\npred.shape","metadata":{"execution":{"iopub.status.busy":"2022-11-03T11:19:34.465350Z","iopub.execute_input":"2022-11-03T11:19:34.465721Z","iopub.status.idle":"2022-11-03T11:19:34.585617Z","shell.execute_reply.started":"2022-11-03T11:19:34.465690Z","shell.execute_reply":"2022-11-03T11:19:34.584562Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"mask = pred.squeeze(0).cpu().detach().numpy()\nprint(mask.shape)\nmask = mask.transpose(1,2,0)\nprint(mask.shape)","metadata":{"execution":{"iopub.status.busy":"2022-11-03T11:19:36.072462Z","iopub.execute_input":"2022-11-03T11:19:36.073050Z","iopub.status.idle":"2022-11-03T11:19:36.084252Z","shell.execute_reply.started":"2022-11-03T11:19:36.073005Z","shell.execute_reply":"2022-11-03T11:19:36.083183Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"display_test_img = test_image[\"image\"].cpu().detach().numpy()\nprint(display_test_img.shape)\ndisplay_test_img = display_test_img.transpose(1,2,0)\ndisplay_test_img.shape","metadata":{"execution":{"iopub.status.busy":"2022-11-03T11:19:37.756805Z","iopub.execute_input":"2022-11-03T11:19:37.757989Z","iopub.status.idle":"2022-11-03T11:19:37.768243Z","shell.execute_reply.started":"2022-11-03T11:19:37.757919Z","shell.execute_reply":"2022-11-03T11:19:37.767223Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"mask[mask < 0]=0\nmask[mask > 0]=1","metadata":{"execution":{"iopub.status.busy":"2022-11-03T11:19:41.605816Z","iopub.execute_input":"2022-11-03T11:19:41.606506Z","iopub.status.idle":"2022-11-03T11:19:41.612243Z","shell.execute_reply.started":"2022-11-03T11:19:41.606469Z","shell.execute_reply":"2022-11-03T11:19:41.611187Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"print(\"-------Original Image-------\")\nplt.imshow(display_test_img, cmap=\"gray\")\nplt.show()\nprint(\"-------Image Mask-------\")\nplt.imshow(mask,cmap=\"gray\")\nplt.show()","metadata":{"execution":{"iopub.status.busy":"2022-11-03T11:19:43.181597Z","iopub.execute_input":"2022-11-03T11:19:43.181972Z","iopub.status.idle":"2022-11-03T11:19:43.632332Z","shell.execute_reply.started":"2022-11-03T11:19:43.181919Z","shell.execute_reply":"2022-11-03T11:19:43.631316Z"},"trusted":true},"execution_count":null,"outputs":[]}]}