关于precision、recall、ap、PRC曲线、ROC曲线及AUC等参数的说明

 原本的代码voc_eval.py给出了prec、rec以及ap的值,在此基础上绘制了PRC曲线、ROC曲线及计算AUC等参数。参考了Faster R-CNN/R-FCN里mAP的计算过程(voc_eval.py解析)

prec:它计算的是所有"正确被检索的item(TP)"占所有"实际被检索到的(TP+FP)"的比例

rec:它计算的是所有"正确被检索的item(TP)"占所有"应该检索到的item(TP+FN)"的比例

要想查看prec及rec的值,在voc_dataset_evaluator.py中_do_python_eval添加以下, 代码注意缩进

        logger.info('REC for {} = {:.4f}'.format(cls, rec[-1]))
        logger.info('PREC for {} = {:.4f}'.format(cls, prec[-1]))

PRC曲线:横坐标为rec,纵坐标为prec,积分面积为ap

同样在上述位置添加

        pl.plot(rec, prec, lw=2, 
              label='Precision-recall curve of class {} (area = {:.4f})'
                    ''.format(cls, ap))
    pl.xlabel('Recall') 
    pl.ylabel('Precision')  
    plt.grid(True) 
    pl.ylim([0.0, 1.05]) 
    pl.xlim([0.0, 1.0]) 
    pl.title('Precision-Recall') 
    pl.legend(loc="upper right") 
    plt.show()

这是我的PRC曲线:

 

ROC曲线:横坐标FPR,纵坐标TPR,其值与REC相同,曲线积分面积为AUC,这一点与PRC曲线类似,AUC越大说明训练的模型性能越好。其中, 曲线的横坐标FPR为假正例率,在voc_eval.py中,tn+fp可由nd-npos得到。

    an=nd - npos
    fpr=fp/an

要计算AUC的值,得编写添加一个函数voc_auc

def voc_auc(rec,fpr):
    mrec = np.concatenate(([0.], rec, [1.]))
    mfpr = np.concatenate(([0.], fpr, [0.]))
    for l in range(mfpr.size - 1, 0, -1):
         mfpr[l - 1] = np.maximum(mfpr[l - 1], mfpr[l])
    l = np.where(mrec[1:] != mrec[:-1])[0]
    auc = np.sum((mrec[l + 1] - mrec[l]) * mfpr[l + 1])
    return auc

从而调用这个函数就可得到AUC值。

绘图命令还是在voc_dataset_evaluator.py中添加

        pl.plot(fpr, rec, lw=2, 
              label='fpr-tpr curve of class {} (area = {:.4f})'
                    ''.format(cls, auc))
    pl.xlabel('FPR') 
    pl.ylabel('TPR')  
    plt.grid(True) 
    pl.ylim([0.0, 1.05]) 
    pl.xlim([0.0, 1.0]) 
    pl.title('FPR-TPR') 
    pl.legend(loc="upper right") 
    plt.show()

我的ROC曲线

 

 

  • 1
    点赞
  • 10
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
好的,以下是给出GPU上训练CIFAR10并计算指标的PyTorch代码: ```python import torch import torch.nn as nn import torch.optim as optim import torchvision import torchvision.transforms as transforms from sklearn.metrics import accuracy_score, precision_score, recall_score, f1_score, roc_curve, auc # 定义网络模型 class Net(nn.Module): def __init__(self): super(Net, self).__init__() self.conv1 = nn.Conv2d(3, 32, 5) self.pool = nn.MaxPool2d(2, 2) self.conv2 = nn.Conv2d(32, 64, 5) self.fc1 = nn.Linear(64 * 5 * 5, 256) self.fc2 = nn.Linear(256, 10) def forward(self, x): x = self.pool(nn.functional.relu(self.conv1(x))) x = self.pool(nn.functional.relu(self.conv2(x))) x = x.view(-1, 64 * 5 * 5) x = nn.functional.relu(self.fc1(x)) x = self.fc2(x) return x # 加载数据集并进行数据增强 transform_train = transforms.Compose([ transforms.RandomCrop(32, padding=4), transforms.RandomHorizontalFlip(), transforms.ToTensor(), transforms.Normalize((0.5, 0.5, 0.5), (0.5, 0.5, 0.5)) ]) transform_test = transforms.Compose([ transforms.ToTensor(), transforms.Normalize((0.5, 0.5, 0.5), (0.5, 0.5, 0.5)) ]) trainset = torchvision.datasets.CIFAR10(root='./data', train=True, download=True, transform=transform_train) trainloader = torch.utils.data.DataLoader(trainset, batch_size=128, shuffle=True, num_workers=2) testset = torchvision.datasets.CIFAR10(root='./data', train=False, download=True, transform=transform_test) testloader = torch.utils.data.DataLoader(testset, batch_size=128, shuffle=False, num_workers=2) # 定义损失函数和优化器 net = Net() criterion = nn.CrossEntropyLoss() optimizer = optim.SGD(net.parameters(), lr=0.001, momentum=0.9) # 将模型移动到GPU上进行训练 device = torch.device("cuda:0" if torch.cuda.is_available() else "cpu") net.to(device) # 训练模型 for epoch in range(20): running_loss = 0.0 for i, data in enumerate(trainloader, 0): inputs, labels = data inputs, labels = inputs.to(device), labels.to(device) optimizer.zero_grad() outputs = net(inputs) loss = criterion(outputs, labels) loss.backward() optimizer.step() running_loss += loss.item() if i % 100 == 99: print('[%d, %5d] loss: %.3f' % (epoch + 1, i + 1, running_loss / 100)) running_loss = 0.0 # 在测试集上进行预测并计算各项指标 y_true = [] y_pred = [] y_score = [] net.eval() with torch.no_grad(): for data in testloader: images, labels = data images, labels = images.to(device), labels.to(device) outputs = net(images) _, predicted = torch.max(outputs.data, 1) y_true += labels.cpu().numpy().tolist() y_pred += predicted.cpu().numpy().tolist() y_score += nn.functional.softmax(outputs, dim=1).cpu().numpy().tolist() acc = accuracy_score(y_true, y_pred) precision = precision_score(y_true, y_pred, average='macro') recall = recall_score(y_true, y_pred, average='macro') f1 = f1_score(y_true, y_pred, average='macro') fpr = {} tpr = {} roc_auc = {} for i in range(10): fpr[i], tpr[i], _ = roc_curve([1 if j == i else 0 for j in y_true], [score[i] for score in y_score]) roc_auc[i] = auc(fpr[i], tpr[i]) print("Accuracy:", acc) print("Precision:", precision) print("Recall:", recall) print("F1 score:", f1) for i in range(10): print("ROC AUC of class", i, ":", roc_auc[i]) ``` 这段代码会输出ACC PRECISION F1 RECALL AUC以及每一类ROC曲线的值。
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值