1. PSNR(峰值信噪比)
PSNR Peak Signal-to-Noise Ratio
衡量图像质量的指标。
维基百科解释
表示信号最大可能功率和影响它的表示精度的破坏性噪声功率的比值的工程术语。由于许多信号都有非常宽的动态范围,峰值信噪比常用对数分贝单位来表示。
计算方法:
- 先计算MSE
M S E = 1 m n ∑ i = 0 m − 1 ∑ j = 0 n − 1 [ I ( i , j ) − K ( i , j ) ] MSE=\frac{1}{mn}\sum^{m-1}_{i=0}\sum^{n-1}_{j=0}[I(i,j)-K(i,j)] MSE=mn1i=0∑m−1j=0∑n−1[I(i,j)−K(i,j)] - 计算PSNR
P S N R = 10 ⋅ l o g 10 M A X I 2 M S E = 20 ⋅ l o g 10 ( M A X I M S E ) PSNR = 10 \cdot log_{10}\frac{MAX^2_I}{MSE}=20\cdot log_{10}(\frac{MAX_I}{\sqrt MSE}) PSNR=10⋅log10MSEMAXI2=20⋅log10(MSEMAXI)
M A X I MAXI MAXI图像中颜色最大的数值,比如对于8位就是255
PSNR 大于40dB,图片质量很好
30——40dB,察觉到失真但是可以接收
20——30dB,图片质量差
<20dB,不可接受
2.SSIM(Structural SIMilarity)结构相似性
SSIM对两张图片比较了三个方面:亮度(Iuminance),对比度(Contrast)和结构(Structure)
l
(
x
,
y
)
=
2
μ
x
μ
y
+
c
1
μ
x
2
+
μ
y
2
+
c
1
l(x,y) = \frac {2\mu_x\mu_y+c_1}{\mu_x^2+\mu_y^2+c_1}
l(x,y)=μx2+μy2+c12μxμy+c1
c
(
x
,
y
)
=
2
σ
x
σ
y
+
c
2
σ
x
2
+
σ
y
2
+
c
2
c(x, y) = \frac{2\sigma_x\sigma_y+c_2}{\sigma_x^2+\sigma_y^2+c_2}
c(x,y)=σx2+σy2+c22σxσy+c2
s
(
x
,
y
)
=
σ
x
y
+
c
3
σ
x
σ
y
+
c
3
s(x, y ) = \frac{\sigma_{xy} + c_3}{\sigma_x\sigma_y + c_3}
s(x,y)=σxσy+c3σxy+c3
一般,取
c
3
=
c
2
/
2
c_3 = c_2 /2
c3=c2/2
- μ x \mu_x μx为x的均值
- μ y \mu_y μy为y的均值
- σ x 2 \sigma_x^2 σx2为x的方差
- σ y 2 \sigma_y^2 σy2为y的方差
- σ x y \sigma_{xy} σxy为x和y的协方差
- c 1 = ( k 1 L ) 2 , c 2 = ( k 2 L ) 2 c_1 =(k_1L)^2,c2=(k_2L)^2 c1=(k1L)2,c2=(k2L)2两个常数,避免分母为0
- L L L像素范围, 2 B − 1 2^B-1 2B−1
- k 1 = 0.01 k1 = 0.01 k1=0.01, k 2 = 0.03 k2=0.03 k2=0.03默认值
S
S
I
(
x
,
y
)
=
[
l
(
x
,
y
)
α
⋅
c
(
x
,
y
)
β
⋅
s
(
x
,
y
)
γ
]
SSI(x,y)=[l(x,y)^{\alpha}\cdot c(x,y)^{\beta} \cdot s(x,y)^{\gamma} ]
SSI(x,y)=[l(x,y)α⋅c(x,y)β⋅s(x,y)γ]
α
β
,
γ
\alpha \beta ,\gamma
αβ,γ设为1
S
S
I
M
(
x
,
y
)
=
(
2
μ
x
μ
y
+
c
1
)
(
2
σ
x
y
+
c
2
)
(
μ
x
2
+
μ
y
2
+
c
1
)
(
σ
x
2
+
σ
y
2
+
c
2
)
SSIM(x,y) = \frac{(2\mu_x\mu_y + c_1)(2\sigma_{xy}+c_2)}{(\mu_x^2 + \mu_y^2 + c_1)(\sigma_x^2 + \sigma_y^2 + c_2)}
SSIM(x,y)=(μx2+μy2+c1)(σx2+σy2+c2)(2μxμy+c1)(2σxy+c2)
计算的时候取 N ∗ N N * N N∗N的窗口计算,最后取平均值
ssim = skimage.measure.compare_ssim(im1, im2, data_range=255)
3、IE(interpolation error)
这个在2020年的cvpr Space-Time-Aware Multi-Resolution Video Enhancement Muhammad看到
计算代码在https://github.com/alterzero/STARnet/blob/master/inter_error.m
function ie=inter_error(im1,im2)
diff_rgb = 128.0 + im2 - im1;
ie = mean(mean(mean(abs(diff_rgb - 128.0))));
4、FR-IQA(full reference image quality assessment)
5、NR-IQA(no reference image quality assessment)
又叫做NQSR,越高越好
6、NIQE
越低越好