| #!/usr/bin/env python |
| ## Copyright (c) 2021, Alliance for Open Media. All rights reserved |
| ## |
| ## This source code is subject to the terms of the BSD 3-Clause Clear License and the |
| ## Alliance for Open Media Patent License 1.0. If the BSD 3-Clause Clear License was |
| ## not distributed with this source code in the LICENSE file, you can obtain it |
| ## at aomedia.org/license/software-license/bsd-3-c-c/. If the Alliance for Open Media Patent |
| ## License 1.0 was not distributed with this source code in the PATENTS file, you |
| ## can obtain it at aomedia.org/license/patent-license/. |
| ## |
| __author__ = "maggie.sun@intel.com, ryanlei@meta.com" |
| |
| import os |
| import re |
| import logging |
| import math |
| from Config import BinPath, LoggerName, VMAF |
| from Utils import GetShortContentName, ExecuteCmd |
| from AV2CTCVideo import CTC_VERSION |
| |
| subloggername = "CalcQtyMetrics_VMAFTool" |
| loggername = LoggerName + '.' + '%s' % subloggername |
| logger = logging.getLogger(loggername) |
| |
| Model_Pkg_File = os.path.join(BinPath, 'vmaf_v0.6.1.pkl') |
| VMAFMetricsFullList = ['VMAF_Y','VMAF_Y-NEG','PSNR_Y','PSNR_U','PSNR_V','SSIM_Y(dB)', |
| 'MS-SSIM_Y(dB)','PSNR-HVS','CIEDE2000','APSNR_Y','APSNR_U','APSNR_V','CAMBI'] |
| |
| def ParseVMAFLogFile(vmaf_log): |
| floats = len(VMAFMetricsFullList) * [0.0] |
| per_frame_log = [] |
| flog = open(vmaf_log, 'r') |
| psnr_y='0'; psnr_cb='0'; psnr_cr='0'; ssim='0'; psnr_hvs='0'; ms_ssim='0';ciede2000='0' |
| frameNum='0'; vmaf='0'; vmaf_neg='0'; cambi='0' |
| for line in flog: |
| m = re.search(r"<frame\s+frameNum=\"(\d+)\"",line) |
| if m: |
| frameNum = m.group(1) |
| m = re.search(r"<frame\s+(.*)\s+psnr_y=\"(\d+\.?\d*)\"",line) |
| if m: |
| psnr_y = m.group(2) |
| m = re.search(r"<frame\s+(.*)\s+psnr_cb=\"(\d+\.?\d*)\"", line) |
| if m: |
| psnr_cb = m.group(2) |
| m = re.search(r"<frame\s+(.*)\s+psnr_cr=\"(\d+\.?\d*)\"", line) |
| if m: |
| psnr_cr = m.group(2) |
| m = re.search(r"<frame\s+(.*)\s+float_ssim=\"(\d+\.?\d*)\"", line) |
| if m: |
| ssim = m.group(2) |
| m = re.search(r"<frame\s+(.*)\s+psnr_hvs=\"(\d+\.?\d*)\"", line) |
| if m: |
| psnr_hvs = m.group(2) |
| m = re.search(r"<frame\s+(.*)\s+float_ms_ssim=\"(\d+\.?\d*)\"", line) |
| if m: |
| ms_ssim = m.group(2) |
| m = re.search(r"<frame\s+(.*)\s+ciede2000=\"(\d+\.?\d*)\"", line) |
| if m: |
| ciede2000 = m.group(2) |
| m = re.search(r"<frame\s+(.*)\s+vmaf=\"(\d+\.?\d*)\"", line) |
| if m: |
| vmaf = m.group(2) |
| m = re.search(r"<frame\s+(.*)\s+cambi=\"(\d+\.?\d*)\"", line) |
| if m: |
| cambi = m.group(2) |
| m = re.search(r"<frame\s+(.*)\s+vmaf_neg=\"(\d+\.?\d*)\"", line) |
| if m: |
| vmaf_neg = m.group(2) |
| per_frame_log.append("%s,%s,%s,%s,%s,%s,%s,%s,%s,%s"% |
| (psnr_y,psnr_cb,psnr_cr,ssim,ms_ssim,vmaf,vmaf_neg,psnr_hvs, |
| ciede2000,cambi)) |
| m = re.search(r"\"vmaf\".*\s+mean=\"(\d+\.?\d*)\"\s+",line) |
| if m: |
| floats[0] = m.group(1) |
| m = re.search(r"\"vmaf_neg\".*\s+mean=\"(\d+\.?\d*)\"\s+", line) |
| if m: |
| floats[1] = m.group(1) |
| m = re.search(r"\"psnr_y\".*\s+mean=\"(\d+\.?\d*)\"\s+", line) |
| if m: |
| floats[2] = m.group(1) |
| m = re.search(r"\"psnr_cb\".*\s+mean=\"(\d+\.?\d*)\"\s+", line) |
| if m: |
| floats[3] = m.group(1) |
| m = re.search(r"\"psnr_cr\".*\s+mean=\"(\d+\.?\d*)\"\s+", line) |
| if m: |
| floats[4] = m.group(1) |
| m = re.search(r"\"float_ssim\".*\s+mean=\"(\d+\.?\d*)\"\s+", line) |
| if m: |
| floats[5] = m.group(1) |
| m = re.search(r"\"float_ms_ssim\".*\s+mean=\"(\d+\.?\d*)\"\s+", line) |
| if m: |
| floats[6] = m.group(1) |
| m = re.search(r"\"psnr_hvs\".*\s+mean=\"(\d+\.?\d*)\"\s+", line) |
| if m: |
| floats[7] = m.group(1) |
| m = re.search(r"\"ciede2000\".*\s+mean=\"(\d+\.?\d*)\"\s+", line) |
| if m: |
| floats[8] = m.group(1) |
| m = re.search(r"\"cambi\".*\s+mean=\"(\d+\.?\d*)\"\s+", line) |
| if m: |
| floats[12] = m.group(1) |
| #<aggregate_metrics apsnr_y="46.817276" apsnr_cb="49.092538" apsnr_cr="50.014785" /> |
| m = re.search(r"aggregate_metrics\s+apsnr_y=\"(\d+\.?\d*)\"\s+apsnr_cb=\"(\d+\.?\d*)\"\s+apsnr_cr=\"(\d+\.?\d*)\"", line) |
| if m: |
| floats[9] = m.group(1) |
| floats[10] = m.group(2) |
| floats[11] = m.group(3) |
| flog.close() |
| floats = [float(i) for i in floats] |
| print_str = "VMAF quality metrics: " |
| for metrics, idx in zip(VMAFMetricsFullList, range(len(VMAFMetricsFullList))): |
| print_str += "%s = %2.5f, " % (metrics, floats[idx]) |
| logger.info(print_str) |
| |
| return floats[0:len(VMAFMetricsFullList)], per_frame_log |
| |
| |
| def GetVMAFLogFile(recfile, path): |
| filename = GetShortContentName(recfile, False) + '_vmaf.log' |
| file = os.path.join(path, filename) |
| return file |
| |
| ################################################################################ |
| ##################### Exposed Functions ######################################## |
| def VMAF_CalQualityMetrics(origfile, recfile, logfilePath, LogCmdOnly=False): |
| vmaf_log = GetVMAFLogFile(recfile, logfilePath) |
| args = " -r %s -d %s -q --threads 4 -o %s" \ |
| % (origfile, recfile, vmaf_log) |
| |
| if CTC_VERSION in ['3.0', '4.0', '5.0']: |
| args += " --aom_ctc v2.0 --feature cambi" |
| else: |
| args += " --aom_ctc v1.0" |
| |
| cmd = VMAF + args |
| ExecuteCmd(cmd, LogCmdOnly) |
| |
| def VMAF_GatherQualityMetrics(recfile, logfilePath): |
| vmaf_log = GetVMAFLogFile(recfile, logfilePath) |
| if not os.path.exists(vmaf_log): |
| return [], [] |
| results, per_frame_log = ParseVMAFLogFile(vmaf_log) |
| return results, per_frame_log |
