python 使用numpy计算混淆矩阵

假如一个模型要预测的类别有三个，分别为A、B、C，使用模型预测测试集得到以下结果:

我们一列一列来看，先看第一列：30、15、5

这里我们的测试集有且只有三个分类A、B、C；也就是真实分类A、B、C就对应着测试集的总体，对于一个样本的预测也只可能是这三者之一。

模型预测值为A的，实际标签不一定就是A，但它一定是A、B、C三者之一，这里预测为A的前提下：真实值为A的有30个、真实值为B的有15个、真实值为C的有5个。

上述表格用numpy表示如下

import numpy as np
# 混淆矩阵
c_matrix = np.array([[30,  7,   3],
                    [ 15, 22,  3],
                    [ 5,   1,  14]])
print(c_matrix.shape)  # (3,3)
print(c_matrix[0][1])  # 7

也可对混淆矩阵进行标准化，使其值在0到1之间

# 混淆矩阵标准化（这里使用L1规范化，是对每一行来说规范化)
print(c_matrix.sum(axis=1))
print(c_matrix.sum(axis=1)[:, np.newaxis])  
c_matrix = c_matrix / c_matrix.sum(axis=1)[:, np.newaxis]
print(c_matrix)

[40 40 20]
[[40]
 [40]
 [20]]
[[0.75  0.175 0.075]
 [0.375 0.55  0.075]
 [0.25  0.05  0.7  ]]

在测试模型时计算

confusion_matrix = np.zeros( (len(class_names), len(class_names)) )  # 混淆矩阵
for images, labels in test_ds.take(total_batch):
        labels = labels.numpy()
        predictions = model.predict(images)
        score = tf.nn.softmax(predictions)
        for index, elem in enumerate(score):
            r, c = np.argmax(elem), labels[index]
            confusion_matrix[r][c] += 1