-
Notifications
You must be signed in to change notification settings - Fork 1
Expand file tree
/
Copy pathtest_only.py
More file actions
191 lines (157 loc) · 8.15 KB
/
test_only.py
File metadata and controls
191 lines (157 loc) · 8.15 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
import os
import argparse
import torch.nn.functional as F
import torch.optim as optim
import torch.utils.data as data
from ignite.contrib.handlers import ProgressBar
from ignite.engine import Engine, Events
from ignite.handlers import ModelCheckpoint, Timer
from ignite.metrics import RunningAverage
from tensorboardX import SummaryWriter
from imgaug import augmenters as iaa
from misc.train_ultils_all_iter import *
import importlib
import dataset as dataset
from config import Config
from loss.ceo_loss import count_pred
from tqdm import tqdm
from sklearn.metrics import precision_score, recall_score, f1_score, accuracy_score, confusion_matrix, cohen_kappa_score
####
class Tester(Config):
def __init__(self, _args=None):
super(Tester, self).__init__(_args=_args)
if _args is not None:
self.__dict__.update(_args.__dict__)
print(self.run_info)
####
def infer_step(self, net, batch, device):
net.eval() # infer mode
imgs, true = batch
imgs = imgs.permute(0, 3, 1, 2) # to NCHW
# push data to GPUs and convert to float32
imgs = imgs.to(device).float()
true = true.to(device).long() # not one-hot
# -----------------------------------------------------------
with torch.no_grad(): # dont compute gradient
out_net = net(imgs, tax=True) # a list contains all the out put of the network
if "CLASS" in self.task_type:
logit_class = out_net
prob = nn.functional.softmax(logit_class, dim=-1)
return dict(logit_c=prob.cpu().numpy(), # from now prob of class task is called by logit_c
true=true.cpu().numpy())
if "REGRESS" in self.task_type:
if "rank_ordinal" in self.loss_type:
logits, probas = out_net[0], out_net[1]
predict_levels = probas > 0.5
pred = torch.sum(predict_levels, dim=1)
return dict(logit_r=pred.cpu().numpy(),
true=true.cpu().numpy())
if "rank_dorn" in self.loss_type:
pred, softmax = net(imgs)
return dict(logit_r=pred.cpu().numpy(),
true=true.cpu().numpy())
if "soft_label" in self.loss_type:
logit_regress = (self.nr_classes - 1) * out_net
return dict(logit_r=logit_regress.cpu().numpy(),
true=true.cpu().numpy())
if "FocalOrdinal" in self.loss_type:
logit_regress = out_net
pred = count_pred(logit_regress)
return dict(logit_r=pred.cpu().numpy(),
true=true.cpu().numpy())
else:
logit_regress = out_net
return dict(logit_r=logit_regress.cpu().numpy(),
true=true.cpu().numpy())
if "MULTI" in self.task_type:
logit_class, logit_regress = out_net[0], out_net[1]
prob = nn.functional.softmax(logit_class, dim=-1)
return dict(logit_c=prob.cpu().numpy(),
logit_r=logit_regress.cpu().numpy(),
true=true.cpu().numpy())
####
def run_once(self, data_root_dir, dataset_=None, fold_idx=None):
log_dir = self.log_dir
check_manual_seed(self.seed)
if dataset_ == 'colon_tma_test_1':
self.dataset = "colon_tma"
_, _, test_pairs = getattr(dataset, ('prepare_colon_tma_data_test_1'))()
elif dataset_ == 'colon_tma_test_2':
self.dataset = "colon_tma"
_, _, test_pairs = getattr(dataset, ('prepare_colon_tma_data_test_2'))()
elif dataset_ == 'prostate_uhu':
self.dataset = "prostate_uhu"
_, _, test_pairs = getattr(dataset, ('prepare_%s_data' % self.dataset))()
elif dataset_ == 'prostate_ubc':
self.dataset = "prostate_ubc"
test_pairs = getattr(dataset, ('prepare_%s_data' % self.dataset))()
elif dataset_ == 'aggc2022':
self.dataset = "aggc2022"
test_pairs = getattr(dataset, ('prepare_%s_data' % self.dataset))()
# --------------------------- Dataloader
infer_augmentors = self.infer_augmentors() # HACK at has_aux
test_dataset = dataset.DatasetSerial(test_pairs[:], has_aux=False,
shape_augs=iaa.Sequential(infer_augmentors[0]))
test_loader = data.DataLoader(test_dataset,
num_workers=self.nr_procs_valid,
batch_size=self.infer_batch_size,
shuffle=False, drop_last=False)
device = 'cuda'
# Define network
net_def = importlib.import_module('model_lib.efficientnet_pytorch.model') # dynamic import
net = net_def.jl_efficientnet(task_mode=self.task_type.lower(), pretrained=True)
PATH_model = args.checkpoint
net = torch.nn.DataParallel(net).to(device)
checkpoint = torch.load(PATH_model)
pytorch_total_params = sum(p.numel() for p in net.parameters() if p.requires_grad)
print('Num params:', pytorch_total_params)
net.load_state_dict(checkpoint)
net.eval()
logits_c = []
trues = []
# Evaluating
with tqdm(desc='Epoch %d - evaluation', unit='it', total=len(test_loader)) as pbar:
for it, (images, gts) in enumerate(iter(test_loader)):
results = self.infer_step(net, (images, gts), device)
logits_c.append(results['logit_c'])
trues.append(results['true'])
pbar.update()
logits_c = np.concatenate(logits_c, axis=0)
trues = np.concatenate(trues)
preds_c = np.argmax(logits_c, axis=-1)
# For class indices 1, 2, 3, 4
if max(trues) == 4:
trues -= 1
print('----------------------------- Predictions by classification head -----------------------------')
print('Precision: ', precision_score(trues, preds_c, average='macro'))
print('Recall: ', recall_score(trues, preds_c, average='macro'))
print('F1: ', f1_score(trues, preds_c, average='macro'))
print('Accuracy: ', accuracy_score(trues, preds_c))
print('Kw:', cohen_kappa_score(trues, preds_c, weights='quadratic'))
print('Confusion matrix: ')
print(confusion_matrix(trues, preds_c))
return
####
def run(self, data_root_dir=None, dataset=None):
self.run_once(data_root_dir, dataset, self.fold_idx)
return
####
if __name__ == '__main__':
parser = argparse.ArgumentParser()
parser.add_argument('--gpu', help='comma separated list of GPU(s) to use.')
parser.add_argument('--view', help='view dataset', action='store_true')
parser.add_argument('--run_info', type=str, default='REGRESS_rank_dorn',
help='CLASS, REGRESS, MULTI + loss, '
'loss ex: MULTI_mtmr, REGRESS_rank_ordinal, REGRESS_rank_dorn'
'REGRESS_FocalOrdinalLoss, REGRESS_soft_ordinal')
parser.add_argument('--dataset', type=str, default='colon_tma', help='colon_tma_test_1, colon_tma_test_2, prostate_uhu, prostate_ubc, panda_512')
parser.add_argument('--data_root_dir', type=str, default='../../anhnguyen/dataset/KBSMC_512_test2/KBSMC_test_2/')
parser.add_argument('--seed', type=int, default=5, help='number')
parser.add_argument('--alpha', type=int, default=5, help='number')
parser.add_argument('--checkpoint', type=str, default='/home/compu/doanhbc/JCO_Learning-pytorch/experiments_dir/log_prostate_uhu_20230327_HUBER_SEESAW_350x350/MULTI_ce_mse_cancer_Effi_seed5_BS64/_net_1550.pth')
parser.add_argument('--log_path', type=str, default='', help='log path')
args = parser.parse_args()
tester = Tester(_args=args)
if args.gpu:
os.environ['CUDA_VISIBLE_DEVICES'] = args.gpu
tester.run(data_root_dir=args.data_root_dir, dataset=args.dataset)