图像质量指标：峰值信噪比PSNR和结构相似性SSIM

图像质量指标

PSNR（Peak Signal-to-Noise Ratio 峰值信噪比）

PSNR（Peak Signal-to-Noise Ratio），峰值信噪比，是一个表示信号最大可能功率和影响它的表示精度的破坏性噪声功率的比值的工程术语。由于许多信号都有非常宽的动态范围，所以峰值信噪比常用对数分贝单位来表示。

以灰度图为例，设大小均为$m\times n$原始图像和处理图像分别为$I_o$和$I_p$，均方误差（MSE）定义为
$$MSE=\frac{1}{mn}\sum _{i=0}^{m-1}\sum _{j=0}^{m-1}[I_o(i,j)-I_p(i,j){]}^2$$
PSNR通过MSE计算，定义为
$$PSNR=10\cdot {\log }_{10}(\frac{{MA{X}_I}^2}{MSE})$$
其中$MA{X}_I$为表示图像像素值的最大值，例如图像一个像素点用8bit来表示，则$MA{X}_I=2^8-1=255$。

处理图像与原始图像越接近，MSE越小，则PSNR越大。一般来说：

• $PSNR>40dB$时，说明图像质量极好（非常接近原始图像）
• $30dB<PSNR<40dB$时，说明图像质量好（存在可以接受的失真）
• $20dB<PSNR<30dB$时，说明图像质量差（存在不可接受的失真）
• $PSNR<20dB$时，说明图像质量极差

对于彩色图像，以RGB图像为例，有以下3种方式计算PSNR：

1. 计算RGB三通道每个通道的PSNR值，再求平均
2. 计算RGB三通道每个通道的MSE值，取平均，再计算PSNR
3. 将RGB转换成YUV颜色空间，仅仅计算Y（亮度）分量的PSNR

其中方法2和方法3较常用。

PSNR是最普遍，最广泛使用的评鉴画质的客观量测法，不过许多实验结果都显示，PSNR的分数无法和人眼看到的视觉品质完全一致，有可能PSNR较高者看起来反而比PSNR较低者差。这是因为人眼的视觉对于误差的敏感度并不是绝对的，其感知结果会受到许多因素的影响而产生变化（例如：人眼对空间频率较低的对比差异敏感度较高，人眼对亮度对比差异的敏感度较色度高，人眼对一个区域的感知结果会受到其周围邻近区域的影响）。

SSIM（Structural Similarity 结构相似性）

SSIM(Structural Similarity)，结构相似性，是一种衡量两幅图像相似度的指标。SSIM 主要考量图片的三个关键特征：亮度（Luminance）,对比度（Contrast）,结构 (Structure)
SSIM使用的两张图像中，一张为未经压缩的无失真图像x，另一张为失真后的图像y。SSIM定义为
$$SSIM(x,y)=l(x,y{)}^{\alpha }\cdot c(x,y{)}^{\beta }\cdot s(x,y{)}^{\gamma }$$
其中$\alpha >0,\beta >0,\gamma >0$

亮度（Luminance）

亮度通过平均所有像素值来衡量
$$u_x=\frac{1}{N}\sum _{i=1}^N{x}_i$$
$$u_y=\frac{1}{N}\sum _{i=1}^N{y}_i$$
$$l(x,y)=\frac{2u_x{u}_y+c_1}{u_x^2+u_y^2+c_1}$$
$c_1$为常数，为了避免分母为0

对比度（Contrast）

对比度通过像素的样本标准差、协方差来衡量
$${\sigma }_x=\sqrt{\frac{{\sum }_{i=1}^N(x_i-u_x{)}^2}{N-1}}$$
$${\sigma }_y=\sqrt{\frac{{\sum }_{i=1}^N(y_i-u_y{)}^2}{N-1}}$$
$${\sigma }_{xy}=\frac{{\sum }_{i=1}^N(x_i-u_x)(y_i-u_y)}{N-1}$$
$$c(x,y)=\frac{2{\sigma }_{xy}+c_2}{{\sigma }_x^2+{\sigma }_y^2+c_2}$$
$c_2$为常数，作用同$c_1$

结构 (Structure)

$$s(x,y)=\frac{{\sigma }_{xy}+c_3}{{\sigma }_x{\sigma }_y+c_3}$$
$c_3$为常数，作用同$c_1$

关于常数$c_1,c_2,c_3$

常数$c_1=(k_{1}L{)}^2$，$c_2=(K_{2}L{)}^2$，$c_3=\frac{c_2}{2}$，经验上常取$k_1=0.01,k_2=0.03$，$L$为像素动态取值范围，即$L=2^{bitsperpixel}-1$（与PSNR中的$MA{X}_I$是一个意思）

实际工程计算中，一般设定$\alpha =\beta =\gamma =1$，SSIM可以简化为公式
$$SSIM(x,y)=\frac{(2u_x{u}_y+c_1)(2{\sigma }_{xy}+c_2)}{(u_x^2+u_y^2+c_1)({\sigma }_x^2+{\sigma }_y^2+c_2)}$$
SSIM满足以下3个条件

1. 对称性：$SSIM(x,y)=SSIM(y,x)$
2. 有界性：$-1\le SSIM(x,y)\le 1$
3. 最大值唯一性：当且仅当$x=y$时，$SSIM=1$

MSSIM（Mean SSIM 平均结构相似性）

在实际应用当中，当需要衡量一整张图片的质量时，可以利用滑动窗将图像分块（blocks）。计算每个block的平均值、标准差、协方差，然后得到每个block的SSIM，再对每个block的SSIM求平均，缺点是会造成块效应（block artifacts）。考虑到滑动窗的分块影响，可以采用高斯加权（圆对称高斯加权）计算每个block的平均值、标准差、协方差，然后得到每个block的SSIM，再对每个block的SSIM求平均即得到MSSIM。

---

参考：

• PSNR-峰值信噪比(原理及python代码实现)
• 图像质量评价指标之 PSNR 和 SSIM
• YCbCr与YUV的区别
• SSIM—结构相似性算法（OpenCV+Python)
• SSIM_百度百科 - Baidu Baike
• SSIM (Structure Similarity Index Measure) 结构衡量指标+代码
• SSIM(结构相似性)-数学公式及python实现