VTM中YUV-PSNR的计算

2023-10-10 02:10
文章标签 计算 psnr yuv vtm

本文主要是介绍VTM中YUV-PSNR的计算,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

名词解释

参见文章cnblogs-什么是SAD,SAE,SATD,SSD,SSE,MAD,MAE,MSD,MSE?

PSNR的计算

较为标准的MSE计算公式和PSNR计算公式如下:
在这里插入图片描述
引用链接:CSDN-图像的峰值信噪比(PSNR)的计算方法
这里需要注意的是,PSNR的计算公式中,应该使用的是MAX(更标准),而不是常见的255,因为255是针对8bit的数据。
对于10bit的数据,不同的工具可能有不同的计算方式,大概在 2 10 − 1 2^{10}-1 2101附近。
VTM中maxval的计算方式如下:

const uint32_t maxval                = /*useWPSNR ? (1 << maximumBitDepth) - 1 :*/ 255 << (maximumBitDepth - 8); 
// fix with WPSNR: 1023 (4095) instead of 1020 (4080) for bit depth 10 (12)

对应的10bit是 255 ∗ 4 = 1020 255*4=1020 2554=1020
对应的12bit是 255 ∗ 16 = = 4080 255*16==4080 25516==4080

YUV-PSRN的计算

VTM中YUV-PSNR的计算并不是三个分量通过加权计算得到的,而是对各个分量的MSE进行加权计算得到的,进而计算得到YUV-PSNR。具体的权重计算分析见代码注释

void calculateCombinedValues(const ChromaFormat chFmt, double &PSNRyuv, double &MSEyuv, const BitDepths &bitDepths){MSEyuv    = 0;int scale = 0;//计算maximumBitDepthint maximumBitDepth = bitDepths.recon[CHANNEL_TYPE_LUMA];for (uint32_t channelTypeIndex = 1; channelTypeIndex < MAX_NUM_CHANNEL_TYPE; channelTypeIndex++){if (bitDepths.recon[channelTypeIndex] > maximumBitDepth){maximumBitDepth = bitDepths.recon[channelTypeIndex];}}
//计算maxval
#if ENABLE_QPAconst uint32_t maxval                = /*useWPSNR ? (1 << maximumBitDepth) - 1 :*/ 255 << (maximumBitDepth - 8); // fix with WPSNR: 1023 (4095) instead of 1020 (4080) for bit depth 10 (12)
#elseconst uint32_t maxval                = 255 << (maximumBitDepth - 8);
#endif//计算分量数const uint32_t numberValidComponents = getNumberValidComponents(chFmt);//遍历各个分量,计算对应的MSEfor (uint32_t comp=0; comp<numberValidComponents; comp++){const ComponentID compID        = ComponentID(comp);//下面会对函数getComponentScaleX部分单独const uint32_t        csx           = getComponentScaleX(compID, chFmt);const uint32_t        csy           = getComponentScaleY(compID, chFmt);const int         scaleChan     = (4>>(csx+csy));	//权重的计算//注意8bit和10bit的计算方式不同const uint32_t        bitDepthShift = 2 * (maximumBitDepth - bitDepths.recon[toChannelType(compID)]); //*2 because this is a squared number//计算分量对应的MSEconst double      channelMSE    = (m_MSEyuvframe[compID] * double(1 << bitDepthShift)) / double(getNumPic());//累加权重scale  += scaleChan;//加权MSEyuv += scaleChan * channelMSE;}//平均MSEyuv /= double(scale);  // i.e. divide by 6 for 4:2:0, 8 for 4:2:2 etc.//计算YUVPSNRyuv = (MSEyuv == 0) ? 999.99 : 10.0 * log10((maxval * maxval) / MSEyuv);}

对于权重的计算函数,参见:

static inline uint32_t        getComponentScaleX        (const ComponentID id, const ChromaFormat fmt) { return getChannelTypeScaleX(toChannelType(id), fmt);               }
static inline uint32_t        getComponentScaleY        (const ComponentID id, const ChromaFormat fmt) { return getChannelTypeScaleY(toChannelType(id), fmt);               }

深入调用

static inline uint32_t        getChannelTypeScaleX      (const ChannelType id, const ChromaFormat fmt) { return (isLuma(id) || (fmt==CHROMA_444)) ? 0 : 1;                  }
static inline uint32_t        getChannelTypeScaleY      (const ChannelType id, const ChromaFormat fmt) { return (isLuma(id) || (fmt!=CHROMA_420)) ? 0 : 1;   

计算流程

SSD -> MSE -> PSNR
MSE的计算代码实现:

const uint32_t size   = width * height;MSEyuvframe[comp] = (double)uiSSDtemp / size;

SSD的计算代码实现

const uint64_t uiSSDtemp = xFindDistortionPlane(recPB, orgPB, useWPSNR ? bitDepth : 0, ::getComponentScaleX(compID, format), ::getComponentScaleY(compID, format));uint64_t EncGOP::xFindDistortionPlane(const CPelBuf& pic0, const CPelBuf& pic1, const uint32_t rshift
#if ENABLE_QPA, const uint32_t chromaShiftHor /*= 0*/, const uint32_t chromaShiftVer /*= 0*/
#endif)
{uint64_t uiTotalDiff;const  Pel*  pSrc0 = pic0.bufAt(0, 0);const  Pel*  pSrc1 = pic1.bufAt(0, 0);TCHECK(pic0.width  != pic1.width , "Unspecified error");TCHECK(pic0.height != pic1.height, "Unspecified error");if( rshift > 0 ){
#if ENABLE_QPAconst   uint32_t  BD = rshift;      // image bit-depthif (BD >= 8){const uint32_t   W = pic0.width;  // image widthconst uint32_t   H = pic0.height; // image heightconst double     R = double(W * H) / (1920.0 * 1080.0);const uint32_t   B = Clip3<uint32_t>(0, 128 >> chromaShiftVer, 4 * uint32_t(16.0 * sqrt(R) + 0.5)); // WPSNR block size in integer multiple of 4 (for SIMD, = 64 at full-HD)uint32_t x, y;if (B < 4) // image is too small to use WPSNR, resort to traditional PSNR{uiTotalDiff = 0;for (y = 0; y < H; y++){for (x = 0; x < W; x++){const           int64_t iDiff = (int64_t)pSrc0[x] - (int64_t)pSrc1[x];uiTotalDiff += uint64_t(iDiff * iDiff);}pSrc0 += pic0.stride;pSrc1 += pic1.stride;}return uiTotalDiff;}double wmse = 0.0, sumAct = 0.0; // compute activity normalized SNR valuefor (y = 0; y < H; y += B){for (x = 0; x < W; x += B){wmse += calcWeightedSquaredError(pic1,   pic0,sumAct, BD,W,      H,x,      y,B,      B);}}// integer weighted distortionsumAct = 16.0 * sqrt ((3840.0 * 2160.0) / double((W << chromaShiftHor) * (H << chromaShiftVer))) * double(1 << BD);return (wmse <= 0.0) ? 0 : uint64_t(wmse * pow(sumAct, BETA) + 0.5);}
#endif // ENABLE_QPAuiTotalDiff = 0;for (int y = 0; y < pic0.height; y++){for (int x = 0; x < pic0.width; x++){Intermediate_Int iTemp = pSrc0[x] - pSrc1[x];uiTotalDiff += uint64_t((iTemp * iTemp) >> rshift);}pSrc0 += pic0.stride;pSrc1 += pic1.stride;}}else{uiTotalDiff = 0;for (int y = 0; y < pic0.height; y++){for (int x = 0; x < pic0.width; x++){Intermediate_Int iTemp = pSrc0[x] - pSrc1[x];uiTotalDiff += uint64_t(iTemp * iTemp);}pSrc0 += pic0.stride;pSrc1 += pic1.stride;}}return uiTotalDiff;
}

这篇关于VTM中YUV-PSNR的计算的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/177183

相关文章

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu

uva 1342 欧拉定理(计算几何模板)

题意: 给几个点,把这几个点用直线连起来,求这些直线把平面分成了几个。 解析: 欧拉定理: 顶点数 + 面数 - 边数= 2。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc

uva 11178 计算集合模板题

题意: 求三角形行三个角三等分点射线交出的内三角形坐标。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <

XTU 1237 计算几何

题面: Magic Triangle Problem Description: Huangriq is a respectful acmer in ACM team of XTU because he brought the best place in regional contest in history of XTU. Huangriq works in a big compa

音视频入门基础:WAV专题(10)——FFmpeg源码中计算WAV音频文件每个packet的pts、dts的实现

一、引言 从文章《音视频入门基础:WAV专题(6)——通过FFprobe显示WAV音频文件每个数据包的信息》中我们可以知道,通过FFprobe命令可以打印WAV音频文件每个packet(也称为数据包或多媒体包)的信息,这些信息包含该packet的pts、dts: 打印出来的“pts”实际是AVPacket结构体中的成员变量pts,是以AVStream->time_base为单位的显

计算数组的斜率,偏移,R2

模拟Excel中的R2的计算。         public bool fnCheckRear_R2(List<double[]> lRear, int iMinRear, int iMaxRear, ref double dR2)         {             bool bResult = true;             int n = 0;             dou

GPU 计算 CMPS224 2021 学习笔记 02

并行类型 (1)任务并行 (2)数据并行 CPU & GPU CPU和GPU拥有相互独立的内存空间,需要在两者之间相互传输数据。 (1)分配GPU内存 (2)将CPU上的数据复制到GPU上 (3)在GPU上对数据进行计算操作 (4)将计算结果从GPU复制到CPU上 (5)释放GPU内存 CUDA内存管理API (1)分配内存 cudaErro

Java - BigDecimal 计算分位(百分位)

日常开发中,如果使用数据库来直接查询一组数据的分位数,就比较简单,直接使用对应的函数就可以了,例如:         PERCENT_RANK() OVER(PARTITION BY 分组列名 ORDER BY 目标列名) AS 目标列名_分位数         如果是需要在代码逻辑部分进行分位数的计算,就需要我们自己写一个工具类来支持计算了 import static ja

OpenStack离线Train版安装系列—2计算节点-环境准备

本系列文章包含从OpenStack离线源制作到完成OpenStack安装的全部过程。 在本系列教程中使用的OpenStack的安装版本为第20个版本Train(简称T版本),2020年5月13日,OpenStack社区发布了第21个版本Ussuri(简称U版本)。 OpenStack部署系列文章 OpenStack Victoria版 安装部署系列教程 OpenStack Ussuri版

新一代车载(E/E)架构下的中央计算载体---HPC软件架构简介

老规矩,分享一段喜欢的文字,避免自己成为高知识低文化的工程师: 屏蔽力是信息过载时代一个人的特殊竞争力,任何消耗你的人和事,多看一眼都是你的不对。非必要不费力证明自己,无利益不试图说服别人,是精神上的节能减排。 无人问津也好,技不如人也罢,你都要试着安静下来,去做自己该做的事.而不是让内心的烦躁、焦虑、毁掉你本就不多的热情和定力。 时间不知不觉中,快要来到夏末秋初。一年又过去了一大半,成