{"metadata":{"kernelspec":{"language":"python","display_name":"Python 3","name":"python3"},"language_info":{"name":"python","version":"3.10.12","mimetype":"text/x-python","codemirror_mode":{"name":"ipython","version":3},"pygments_lexer":"ipython3","nbconvert_exporter":"python","file_extension":".py"},"kaggle":{"accelerator":"gpu","dataSources":[{"sourceId":45867,"databundleVersionId":6924515,"sourceType":"competition"},{"sourceId":7253378,"sourceType":"datasetVersion","datasetId":4202787},{"sourceId":7253829,"sourceType":"datasetVersion","datasetId":4203111},{"sourceId":7258084,"sourceType":"datasetVersion","datasetId":4206058},{"sourceId":7260435,"sourceType":"datasetVersion","datasetId":4207669},{"sourceId":7261045,"sourceType":"datasetVersion","datasetId":4208092},{"sourceId":7263797,"sourceType":"datasetVersion","datasetId":4210032},{"sourceId":7263977,"sourceType":"datasetVersion","datasetId":4210149},{"sourceId":3729,"sourceType":"modelInstanceVersion","modelInstanceId":2656},{"sourceId":3731,"sourceType":"modelInstanceVersion","modelInstanceId":2658},{"sourceId":3736,"sourceType":"modelInstanceVersion","modelInstanceId":2663},{"sourceId":3737,"sourceType":"modelInstanceVersion","modelInstanceId":2664}],"dockerImageVersionId":30627,"isInternetEnabled":false,"language":"python","sourceType":"notebook","isGpuEnabled":true}},"nbformat_minor":4,"nbformat":4,"cells":[{"cell_type":"code","source":"import os\nimport gc\nimport cv2\nimport math\nimport copy\nimport time\nimport random\nimport glob\nfrom PIL import Image\nfrom matplotlib import pyplot as plt\n\n# For data manipulation\nimport numpy as np\nimport pandas as pd\n\n# Pytorch Imports\nimport torch\nimport torch.nn as nn\nimport torch.optim as optim\nimport torch.nn.functional as F\nfrom torch.optim import lr_scheduler\nfrom torch.utils.data import Dataset, DataLoader\nfrom torch.cuda import amp\nimport torchvision\n\n# Utils\nimport joblib\nfrom tqdm import tqdm\nfrom collections import defaultdict\n\n# Sklearn Imports\nfrom sklearn.preprocessing import LabelEncoder\nfrom sklearn.model_selection import StratifiedKFold\n\n# For Image Models\nimport timm\n\n# Albumentations for augmentations\nimport albumentations as A\nfrom albumentations.pytorch import ToTensorV2\n\n# For colored terminal text\nfrom colorama import Fore, Back, Style\nb_ = Fore.BLUE\nsr_ = Style.RESET_ALL\n\nimport warnings\nwarnings.filterwarnings(\"ignore\")\n\n# For descriptive error messages\nos.environ['CUDA_LAUNCH_BLOCKING'] = \"1\"\nimport os\nimport gc\nimport cv2\nimport math\nimport copy\nimport time\nimport random\nimport glob\nfrom matplotlib import pyplot as plt\n\n# For data manipulation\nimport numpy as np\nimport pandas as pd\n\n# Pytorch Imports\nimport torch\nimport torch.nn as nn\nimport torch.optim as optim\nimport torch.nn.functional as F\nfrom torch.optim import lr_scheduler\nfrom torch.utils.data import Dataset, DataLoader\nfrom torch.cuda import amp\nimport torchvision\n\n# Utils\nimport joblib\nfrom tqdm import tqdm\nfrom collections import defaultdict\n\n# Sklearn Imports\nfrom sklearn.preprocessing import LabelEncoder\nfrom sklearn.model_selection import StratifiedKFold\n\n# For Image Models\nimport timm\n\n# Albumentations for augmentations\nimport albumentations as A\nfrom albumentations.pytorch import ToTensorV2\n\n# For colored terminal text\nfrom colorama import Fore, Back, Style\nb_ = Fore.BLUE\nsr_ = Style.RESET_ALL\n\nimport warnings\nwarnings.filterwarnings(\"ignore\")\n\n# For descriptive error messages\nos.environ['CUDA_LAUNCH_BLOCKING'] = \"1\"","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:32.489461Z","iopub.execute_input":"2023-12-23T02:30:32.490272Z","iopub.status.idle":"2023-12-23T02:30:39.445006Z","shell.execute_reply.started":"2023-12-23T02:30:32.490234Z","shell.execute_reply":"2023-12-23T02:30:39.44402Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"CONFIG = {\n    \"seed\": 42,\n    \"img_size\": 2048,\n    \"model_name\": \"tf_efficientnet_b0_ns\",\n    \"num_classes\": 5,\n    \"valid_batch_size\": 1,\n    \"device\": torch.device(\"cuda:0\" if torch.cuda.is_available() else \"cpu\"),\n}","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:39.446641Z","iopub.execute_input":"2023-12-23T02:30:39.447093Z","iopub.status.idle":"2023-12-23T02:30:39.474116Z","shell.execute_reply.started":"2023-12-23T02:30:39.447065Z","shell.execute_reply":"2023-12-23T02:30:39.47327Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"CONFIG","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:39.475424Z","iopub.execute_input":"2023-12-23T02:30:39.475812Z","iopub.status.idle":"2023-12-23T02:30:39.496424Z","shell.execute_reply.started":"2023-12-23T02:30:39.475777Z","shell.execute_reply":"2023-12-23T02:30:39.49558Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"def set_seed(seed=42):\n    '''Sets the seed of the entire notebook so results are the same every time we run.\n    This is for REPRODUCIBILITY.'''\n    np.random.seed(seed)\n    torch.manual_seed(seed)\n    torch.cuda.manual_seed(seed)\n    # When running on the CuDNN backend, two further options must be set\n    torch.backends.cudnn.deterministic = True\n    torch.backends.cudnn.benchmark = False\n    # Set a fixed value for the hash seed\n    os.environ['PYTHONHASHSEED'] = str(seed)\n    \nset_seed(CONFIG['seed'])","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:39.499059Z","iopub.execute_input":"2023-12-23T02:30:39.499476Z","iopub.status.idle":"2023-12-23T02:30:39.509325Z","shell.execute_reply.started":"2023-12-23T02:30:39.499442Z","shell.execute_reply":"2023-12-23T02:30:39.508526Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"ROOT_DIR = '/kaggle/input/UBC-OCEAN'\nTEST_DIR = '/kaggle/input/UBC-OCEAN/test_thumbnails'\n\nLABEL_ENCODER_BIN = \"/kaggle/input/label-encoder-pkl/label_encoder.pkl\"\n#BEST_WEIGHT = \"/kaggle/input/acc0-66-loss0-8905-epoch9-bin/Acc0.66_Loss0.8905_epoch9.bin\"\n#BEST_WEIGHT = \"/kaggle/input/loss0-0000-epoch5-bin/Acc1.00_Loss0.0000_epoch5.bin\"\nBEST_WEIGHT='/kaggle/input/acc0-92-loss0-2957-epoch10-bin/Acc0.92_Loss0.2957_epoch10.bin'","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:39.510317Z","iopub.execute_input":"2023-12-23T02:30:39.510589Z","iopub.status.idle":"2023-12-23T02:30:39.515572Z","shell.execute_reply.started":"2023-12-23T02:30:39.510566Z","shell.execute_reply":"2023-12-23T02:30:39.514641Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"def get_test_file_path(image_id):\n    return f\"{TEST_DIR}/{image_id}_thumbnail.png\"","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:39.516597Z","iopub.execute_input":"2023-12-23T02:30:39.516931Z","iopub.status.idle":"2023-12-23T02:30:39.524995Z","shell.execute_reply.started":"2023-12-23T02:30:39.516906Z","shell.execute_reply":"2023-12-23T02:30:39.524136Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"df = pd.read_csv(f\"{ROOT_DIR}/test.csv\")\ndf['file_path'] = df['image_id'].apply(get_test_file_path)\ndf['label'] = 0 # dummy\ndf","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:39.526109Z","iopub.execute_input":"2023-12-23T02:30:39.526443Z","iopub.status.idle":"2023-12-23T02:30:39.555835Z","shell.execute_reply.started":"2023-12-23T02:30:39.526401Z","shell.execute_reply":"2023-12-23T02:30:39.554779Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"df_sub = pd.read_csv(f\"{ROOT_DIR}/sample_submission.csv\")\ndf_sub","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:39.557211Z","iopub.execute_input":"2023-12-23T02:30:39.557583Z","iopub.status.idle":"2023-12-23T02:30:39.568954Z","shell.execute_reply.started":"2023-12-23T02:30:39.557551Z","shell.execute_reply":"2023-12-23T02:30:39.568079Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"encoder = joblib.load( LABEL_ENCODER_BIN )","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:39.570116Z","iopub.execute_input":"2023-12-23T02:30:39.570438Z","iopub.status.idle":"2023-12-23T02:30:39.580862Z","shell.execute_reply.started":"2023-12-23T02:30:39.570406Z","shell.execute_reply":"2023-12-23T02:30:39.579945Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"def get_cropped_images(file_path, image_id, th_area = 1000):\n    image = Image.open(file_path)\n    # Aspect ratio\n    as_ratio = image.size[0] / image.size[1]\n    \n    sxs, exs, sys, eys = [],[],[],[]\n    if as_ratio >= 1.5:\n        # Crop\n        mask = np.max( np.array(image) > 0, axis=-1 ).astype(np.uint8)\n        retval, labels = cv2.connectedComponents(mask)\n        if retval >= as_ratio:\n            x, y = np.meshgrid( np.arange(image.size[0]), np.arange(image.size[1]) )\n            for label in range(1, retval):\n                area = np.sum(labels == label)\n                if area < th_area:\n                    continue\n                xs, ys= x[ labels == label ], y[ labels == label ]\n                sx, ex = np.min(xs), np.max(xs)\n                cx = (sx + ex) // 2\n                crop_size = image.size[1]\n                sx = max(0, cx-crop_size//2)\n                ex = min(sx + crop_size - 1, image.size[0]-1)\n                sx = ex - crop_size + 1\n                sy, ey = 0, image.size[1]-1\n                sxs.append(sx)\n                exs.append(ex)\n                sys.append(sy)\n                eys.append(ey)\n        else:\n            crop_size = image.size[1]\n            for i in range(int(as_ratio)):\n                sxs.append( i * crop_size )\n                exs.append( (i+1) * crop_size - 1 )\n                sys.append( 0 )\n                eys.append( crop_size - 1 )\n    else:\n        # Not Crop (entire image)\n        sxs, exs, sys, eys = [0,],[image.size[0]-1],[0,],[image.size[1]-1]\n\n    df_crop = pd.DataFrame()\n    df_crop[\"image_id\"] = [image_id] * len(sxs)\n    df_crop[\"file_path\"] = [file_path] * len(sxs)\n    df_crop[\"sx\"] = sxs\n    df_crop[\"ex\"] = exs\n    df_crop[\"sy\"] = sys\n    df_crop[\"ey\"] = eys\n    return df_crop","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:39.584689Z","iopub.execute_input":"2023-12-23T02:30:39.585Z","iopub.status.idle":"2023-12-23T02:30:39.59769Z","shell.execute_reply.started":"2023-12-23T02:30:39.584976Z","shell.execute_reply":"2023-12-23T02:30:39.596727Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"dfs = []\nfor (file_path, image_id) in zip(df[\"file_path\"], df[\"image_id\"]):\n    dfs.append( get_cropped_images(file_path, image_id) )\n\ndf_crop = pd.concat(dfs)\ndf_crop[\"label\"] = 0 # dummy\ndf_crop","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:39.598849Z","iopub.execute_input":"2023-12-23T02:30:39.599122Z","iopub.status.idle":"2023-12-23T02:30:40.1902Z","shell.execute_reply.started":"2023-12-23T02:30:39.599099Z","shell.execute_reply":"2023-12-23T02:30:40.189193Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"df_crop = df_crop.drop_duplicates(subset=[\"image_id\", \"sx\", \"ex\", \"sy\", \"ey\"]).reset_index(drop=True)\ndf_crop","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:40.191328Z","iopub.execute_input":"2023-12-23T02:30:40.191656Z","iopub.status.idle":"2023-12-23T02:30:40.209118Z","shell.execute_reply.started":"2023-12-23T02:30:40.191628Z","shell.execute_reply":"2023-12-23T02:30:40.208216Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"class UBCDataset(Dataset):\n    def __init__(self, df, transforms=None):\n        self.df = df\n        self.file_names = df['file_path'].values\n        self.labels = df['label'].values\n        self.transforms = transforms\n        self.sxs = df[\"sx\"].values\n        self.exs = df[\"ex\"].values\n        self.sys = df[\"sy\"].values\n        self.eys = df[\"ey\"].values\n        \n    def __len__(self):\n        return len(self.df)\n    \n    def __getitem__(self, index):\n        img_path = self.file_names[index]\n        sx = self.sxs[index]\n        ex = self.exs[index]\n        sy = self.sys[index]\n        ey = self.eys[index]\n        img = cv2.imread(img_path)\n        img = cv2.cvtColor(img, cv2.COLOR_BGR2RGB)\n        label = self.labels[index]\n        \n        img = img[ sy:ey, sx:ex, : ]\n        \n        if self.transforms:\n            img = self.transforms(image=img)[\"image\"]\n            \n        return {\n            'image': img,\n            'label': torch.tensor(label, dtype=torch.long)\n        }","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:40.210282Z","iopub.execute_input":"2023-12-23T02:30:40.210618Z","iopub.status.idle":"2023-12-23T02:30:40.220816Z","shell.execute_reply.started":"2023-12-23T02:30:40.210587Z","shell.execute_reply":"2023-12-23T02:30:40.219951Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"data_transforms = {\n    \"valid\": A.Compose([\n        A.Resize(CONFIG['img_size'], CONFIG['img_size']),\n        A.Normalize(\n                mean=[0.485, 0.456, 0.406], \n                std=[0.229, 0.224, 0.225], \n                max_pixel_value=255.0, \n                p=1.0\n            ),\n        ToTensorV2()], p=1.)\n}","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:40.221992Z","iopub.execute_input":"2023-12-23T02:30:40.222623Z","iopub.status.idle":"2023-12-23T02:30:40.231865Z","shell.execute_reply.started":"2023-12-23T02:30:40.222589Z","shell.execute_reply":"2023-12-23T02:30:40.231138Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"class GeM(nn.Module):\n    def __init__(self, p=3, eps=1e-6):\n        super(GeM, self).__init__()\n        self.p = nn.Parameter(torch.ones(1)*p)\n        self.eps = eps\n\n    def forward(self, x):\n        return self.gem(x, p=self.p, eps=self.eps)\n        \n    def gem(self, x, p=3, eps=1e-6):\n        return F.avg_pool2d(x.clamp(min=eps).pow(p), (x.size(-2), x.size(-1))).pow(1./p)\n        \n    def __repr__(self):\n        return self.__class__.__name__ + \\\n                '(' + 'p=' + '{:.4f}'.format(self.p.data.tolist()[0]) + \\\n                ', ' + 'eps=' + str(self.eps) + ')'","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:40.232777Z","iopub.execute_input":"2023-12-23T02:30:40.233035Z","iopub.status.idle":"2023-12-23T02:30:40.242903Z","shell.execute_reply.started":"2023-12-23T02:30:40.233013Z","shell.execute_reply":"2023-12-23T02:30:40.24208Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"class UBCModel(nn.Module):\n    def __init__(self, model_name, num_classes, pretrained=False, checkpoint_path=None):\n        super(UBCModel, self).__init__()\n        self.model = timm.create_model(model_name, pretrained=pretrained)\n\n        in_features = self.model.classifier.in_features\n        self.model.classifier = nn.Identity()\n        self.model.global_pool = nn.Identity()\n        self.pooling = GeM()\n        # 追加のLinearレイヤーとDropout\n        self.linear1 = nn.Linear(in_features, 512)\n        self.dropout = nn.Dropout(p=0.3)\n        self.linear2 = nn.Linear(512, num_classes)\n        self.softmax = nn.Softmax(dim=1)\n\n    def forward(self, images):\n        features = self.model(images)\n        pooled_features = self.pooling(features).flatten(1)\n        # 新しいLinearレイヤーを使用\n        x = self.linear1(pooled_features)\n        x = self.dropout(x)\n        output = self.linear2(x)\n\n        return output\n\n    \n\nmodel = UBCModel(CONFIG['model_name'], CONFIG['num_classes'])\nmodel.load_state_dict(torch.load( BEST_WEIGHT ))\nmodel.to(CONFIG['device']);","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:40.243813Z","iopub.execute_input":"2023-12-23T02:30:40.244055Z","iopub.status.idle":"2023-12-23T02:30:43.588876Z","shell.execute_reply.started":"2023-12-23T02:30:40.244033Z","shell.execute_reply":"2023-12-23T02:30:43.587865Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"test_dataset = UBCDataset(df_crop, transforms=data_transforms[\"valid\"])\ntest_loader = DataLoader(test_dataset, batch_size=CONFIG['valid_batch_size'], \n                          num_workers=2, shuffle=False, pin_memory=True)","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:43.590063Z","iopub.execute_input":"2023-12-23T02:30:43.590354Z","iopub.status.idle":"2023-12-23T02:30:43.596098Z","shell.execute_reply.started":"2023-12-23T02:30:43.590328Z","shell.execute_reply":"2023-12-23T02:30:43.595037Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"preds = []\nwith torch.no_grad():\n    bar = tqdm(enumerate(test_loader), total=len(test_loader))\n    for step, data in bar:        \n        images = data['image'].to(CONFIG[\"device\"], dtype=torch.float)        \n        batch_size = images.size(0)\n        outputs = model(images)\n        outputs = model.softmax(outputs)\n        preds.append( outputs.detach().cpu().numpy() )\n\npreds = np.vstack(preds)\nprint(preds.shape)","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:43.597312Z","iopub.execute_input":"2023-12-23T02:30:43.597562Z","iopub.status.idle":"2023-12-23T02:30:48.843754Z","shell.execute_reply.started":"2023-12-23T02:30:43.597539Z","shell.execute_reply":"2023-12-23T02:30:48.84272Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"for i in range(preds.shape[-1]):\n    df_crop[f\"cat{i}\"] = preds[:, i]\n\ndict_label = {}\nfor image_id, gdf in df_crop.groupby(\"image_id\"):\n    dict_label[image_id] = np.argmax( gdf[ [f\"cat{i}\" for i in range(preds.shape[-1])] ].values.max(axis=0) )\n    #dict_label[image_id] = np.argmax( gdf[ [f\"cat{i}\" for i in range(preds.shape[-1])] ].values.mean(axis=0) )\npreds = np.array( [ dict_label[image_id] for image_id in df[\"image_id\"].values ] )","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:48.845201Z","iopub.execute_input":"2023-12-23T02:30:48.845495Z","iopub.status.idle":"2023-12-23T02:30:48.861503Z","shell.execute_reply.started":"2023-12-23T02:30:48.845466Z","shell.execute_reply":"2023-12-23T02:30:48.860665Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"pred_labels = encoder.inverse_transform( preds )\ndf_sub[\"label\"] = pred_labels\ndf_sub.to_csv(\"submission.csv\", index=False)","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:48.862664Z","iopub.execute_input":"2023-12-23T02:30:48.863014Z","iopub.status.idle":"2023-12-23T02:30:48.880536Z","shell.execute_reply.started":"2023-12-23T02:30:48.862986Z","shell.execute_reply":"2023-12-23T02:30:48.879751Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"df_sub","metadata":{"execution":{"iopub.status.busy":"2023-12-23T02:30:48.88174Z","iopub.execute_input":"2023-12-23T02:30:48.882045Z","iopub.status.idle":"2023-12-23T02:30:48.892706Z","shell.execute_reply.started":"2023-12-23T02:30:48.88202Z","shell.execute_reply":"2023-12-23T02:30:48.891879Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"","metadata":{},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"","metadata":{},"execution_count":null,"outputs":[]}]}