data_logging_utils.py 9.03 KB
Newer Older
1
2
3
4
"""Data Logging Functions

Description:

5
    This folder contains several functions which, either on their own or included in larger pieces of software, perform data logging tasks.
6

7
8
9
10
11
Usage:

    To use content from this folder, import the functions and instantiate them as you wish to use them:

        from utils.data_logging_utils import function_name
12
13
14

"""

15
import os
16
import matplotlib
17
18
19
import matplotlib.pyplot as plt
import shutil
import logging
20
import numpy as np
21
import re
22
from textwrap import wrap
23
import torch
24
25
26

# The SummaryWriter class provides a high-level API to create an event file in a given directory and add summaries and events to it.
# More here: https://tensorboardx.readthedocs.io/en/latest/tensorboard.html
27

28
from tensorboardX import SummaryWriter
29

30
31
import utils.data_evaluation_utils as evaluation

32
33
plt.axis('scaled')

34

35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
class LogWriter():

    """Log Writer class for the BrainMapper U-net.

    This class contains the pytorch implementation of the several logging functions required for the BrainMapper project.
    These functions are designed to keep track of progress during training, and also aid debugging.

    Args:
        number_of_classes (int): Number of classes
        logs_directory (str): Directory for outputing training logs
        experiment_name (str): Name of the experiment
        use_last_checkpoint (bool): Flag for loading the previous checkpoint
        labels (arr): Vector/Array of labels (if applicable)
        confusion_matrix_cmap (class): Colour Map to be used for the Conusion Matrix
    """

51
52
    def __init__(self, number_of_classes, logs_directory, experiment_name, use_last_checkpoint=False, labels=None, confusion_matrix_cmap=plt.cm.Blues):

53
        self.number_of_classes = number_of_classes
54
55
        training_logs_directory = os.path.join(
            logs_directory, experiment_name, "train")
56
57
        validation_logs_directory = os.path.join(
            logs_directory, experiment_name, "validation")
58
59
60
61
62

        # If the logs directory exist, we clear their contents to allow new logs to be created
        if not use_last_checkpoint:
            if os.path.exists(training_logs_directory):
                shutil.rmtree(training_logs_directory)
63
64
            if os.path.exists(validation_logs_directory):
                shutil.rmtree(validation_logs_directory)
65
66

        self.log_writer = {
67
68
            'train': SummaryWriter(logdir=training_logs_directory),
            'validation': SummaryWriter(logdir=validation_logs_directory)
69
        }
70
71
72
73
74

        self.confusion_matrix_color_map = confusion_matrix_cmap

        self.current_iteration = 1

75
76
77
        if labels is not None:
            self.labels = self.labels_generator(labels)
        else:
78
            self.labels = ['rsfMRI']
79
80

        self.logger = logging.getLogger()
81
82
        file_handler = logging.FileHandler(
            "{}/{}.log".format(os.path.join(logs_directory, experiment_name), "console_logs"))
83
84
        self.logger.addHandler(file_handler)

85
86
87
88
89
90
91
92
93
    def log(self, message):
        """Log function

        This function logs a message in the logger.

        Args:
            message (str): Message to be logged
        """

94
        self.logger.info(msg=message)
95

96
    def loss_per_iteration(self, loss_per_iteration, batch_index, iteration):
97
        """Log of loss / iteration
98

99
100
101
102
103
        This function records the loss for every iteration.

        Args:
            loss_per_iteration (torch.tensor): Value of loss for every iteration step
            batch_index (int): Index of current batch
104
            iteration (int): Current iteration value
105
106
        """

107
108
109
        print("Loss for Iteration {} is: {}".format(
            batch_index, loss_per_iteration))
        self.log_writer['train'].add_scalar(
Andrei-Claudiu Roibu's avatar
Andrei-Claudiu Roibu committed
110
            'loss/iteration', loss_per_iteration, iteration)
111

112
    def loss_per_epoch(self, losses, phase, epoch, previous_loss=None):
113
114
115
116
117
118
119
120
        """Log function

        This function records the loss for every epoch.

        Args:
            losses (list): Values of all the losses recorded during the training epoch
            phase (str): Current run mode or phase
            epoch (int): Current epoch value
121
            previous_loss(float): Value of the previous epoch's loss
122
123
        """

124
125
        loss = np.mean(losses)

126
        if phase == 'train':
127
128
            # loss = losses[-1]
            print("Loss for Epoch {} of {} is: {}".format(epoch, phase, loss))
129
        else:
130
131
132
133
134
            # loss = np.mean(losses)
            if previous_loss == None:
                print("Loss for Epoch {} of {} is: {}".format(epoch, phase, loss))
            else:
                print("Loss for Epoch {} of {} is {} and Absolute Change is {}".format(epoch, phase, loss, previous_loss - loss))
135

Andrei-Claudiu Roibu's avatar
Andrei-Claudiu Roibu committed
136
        self.log_writer[phase].add_scalar('loss/epoch', loss, epoch)
137

Andrei-Claudiu Roibu's avatar
Andrei-Claudiu Roibu committed
138
139
    def close(self):
        """Close the log writer
140

Andrei-Claudiu Roibu's avatar
Andrei-Claudiu Roibu committed
141
142
        This function closes the two log writers.
        """
143

Andrei-Claudiu Roibu's avatar
Andrei-Claudiu Roibu committed
144
145
146
147
148
149
150
151
152
        self.log_writer['train'].close()
        self.log_writer['validation'].close()

    def add_graph(self, model):
        """Produces network graph

        This function produces the network graph

        NOTE: Currently, the function suffers from bugs and is not implemented.
153

154
        Args:
Andrei-Claudiu Roibu's avatar
Andrei-Claudiu Roibu committed
155
            model (torch.nn.Module): Model to draw.
156
157
        """

Andrei-Claudiu Roibu's avatar
Andrei-Claudiu Roibu committed
158
159
160
        self.log_writer['train'].add_graph(model)

    # DEPRECATED / UNDEBUGGED FUNCTIONS
161

162
    def plot_dice_score(self, dice_score, phase, plot_name, title, epochs=None):
163
        """Function plotting dice score for multiple epochs
164

165
        This function plots the dice score for each epoch.
166

167
168
169
170
171
172
173
174
        Args:
            dice_score (torch.tensor): Dice score value for each class
            phase (str): Current run mode or phase
            plot_name (str): Caption name for later refference
            title (str): Plot title
            epoch (int): Current epoch value
        """

175
        figure = matplotlib.figure.Figure()  # Might add some arguments here later
176
177
178
179
180
        ax = figure.add_subplot(1, 1, 1)
        ax.set_xlabel(title)
        ax.xaxis.set_label_position('top')
        ax.bar(np.arange(self.number_of_classes), dice_score)
        ax.set_xticks(np.arange(self.number_of_classes))
181
182
183
184
185
186

        if self.labels is None:
            pass
        else:
            ax.set_xticklabels(self.labels)
            ax.xaxis.tick_bottom()
187

188
        if epochs:
189
190
            self.log_writer[phase].add_figure(
                plot_name + '/' + phase, figure, global_step=epochs)
191
192
        else:
            self.log_writer[phase].add_figure(plot_name + '/' + phase, figure)
193

Andrei-Claudiu Roibu's avatar
Andrei-Claudiu Roibu committed
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
    def dice_score_per_epoch(self, phase, outputs, correct_labels, epoch):
        """Function calculating dice score for each epoch

        This function computes the dice score for each epoch.

        Args:
            phase (str): Current run mode or phase
            outputs (torch.tensor): Tensor of all the network outputs (Y-hat)
            correct_labels (torch.tensor): Output ground-truth labelled data (Y)
            epoch (int): Current epoch value
        """

        print("Dice Score is being calculated...", end='', flush=True)
        dice_score = evaluation.dice_score_calculator(
            outputs, correct_labels, self.number_of_classes)
        mean_dice_score = torch.mean(dice_score)
        self.plot_dice_score(
            dice_score, phase, plot_name='dice_score_per_epoch', title='Dice Score', epochs=epoch)
        print("Dice Score calculated successfully")
        return mean_dice_score.item()
214

215
216
217
218
219
220
221
222
223
224
225
226
    def sample_image_per_epoch(self, prediction, ground_truth, phase, epoch):
        """Function plotting mirrored images

        This function plots a predicted and a grond truth images side-by-side.

        Args:
            prediction (torch.tensor): Predicted image after passing throught the network
            ground_truth (torch.tensor): Labelled ground truth image
            phase (str): Current run mode or phase
            epoch (int): Current epoch value
        """

227
228
        print("Sample Image is being loaded...", end='', flush=True)
        figure, ax = plt.subplots(nrows=len(prediction), ncols=2)
229
230
231
232
233
234
235
236
237
238
239

        for i in range(len(prediction)):
            ax[i][0].imshow(prediction[i])
            ax[i][0].set_title("Predicted Image")
            ax[i][0].axis('off')

            ax[i][1].imshow(ground_truth[i])
            ax[i][1].set_title('Ground Truth Image')
            ax[i][1].axis('off')

        figure.set_tight_layout()
240
241
        self.log_writer[phase].add_figure(
            'sample_prediction/'+phase, figure, epoch)
242
243

        print("Sample Image successfully loaded!")
244

Andrei-Claudiu Roibu's avatar
Andrei-Claudiu Roibu committed
245
246
    def labels_generator(self, labels):
        """ Label Generator Function
Andrei-Claudiu Roibu's avatar
Andrei-Claudiu Roibu committed
247

Andrei-Claudiu Roibu's avatar
Andrei-Claudiu Roibu committed
248
249
250
251
252
253
254
        This function processess an input array of labels.

        Args:
            labels (arr): Vector/Array of labels (if applicable)

        Returns:
            label_classes (list): List of processed labels
Andrei-Claudiu Roibu's avatar
Andrei-Claudiu Roibu committed
255
        """
256

Andrei-Claudiu Roibu's avatar
Andrei-Claudiu Roibu committed
257
258
259
260
261
262
263
264
265
266
267
        label_classes = []

        for label in labels:

            label_class = re.sub(
                r'([a-z](?=[A-Z])|[A-Z](?=[A-Z][a-z]))', r'\1 ', label)
            label_class = ['\n'.join(wrap(element, 40))
                           for element in label_class]
            label_classes.append(label_class)

        return label_classes