python计算psnr和ssim (用于超分辨率)

python计算psnr和ssim (用于超分辨率)

from utils import utils_image as util
import os
import cv2

HR_path = 'dataset/benchmark/Set5/HR'
SR_path = 'experiments/results/Set5/x4'


n_channels = 3

def evulate():
    hr_paths = util.get_image_paths(HR_path)
    numbers = len(hr_paths)
    sum_psnr = 0
    max_psnr = 0
    min_psnr = 100
    sum_ssim = 0
    max_ssim = 0
    min_ssim = 1
    for hr_path in hr_paths:
        # img_name, ext = os.path.splitext(os.path.basename(img_path))
        img_name = os.path.basename(hr_path)
        sr_path = os.path.join(SR_path,img_name)
        print(img_name)
        # print(hr_path)
        # print(sr_path)
        img_Hr = util.imread_uint(hr_path, n_channels=n_channels)  # HR image, int8
        img_Sr = util.imread_uint(sr_path, n_channels=n_channels)  # HR image, int8
        psnr = util.calculate_psnr(img_Sr, img_Hr,)
        print(psnr)
        sum_psnr += psnr
        max_psnr = max(max_psnr,psnr)
        min_psnr = min(min_psnr, psnr)
        ssim = util.calculate_ssim(img_Sr, img_Hr,)
        # print(ssim)
        sum_ssim += ssim
        max_ssim = max(max_ssim,ssim)
        min_ssim = min(min_ssim, ssim)
    print('Average psnr = ', sum_psnr / numbers)
    print('min_psnr = ', min_psnr)
    print('Max_psnr = ', max_psnr)
    print('Average ssim = ', sum_ssim / numbers)
    print('min_ssim = ', min_ssim)
    print('Max_ssim = ', max_ssim)


def evulate_diff_name():
    hr_paths = util.get_image_paths(HR_path)
    numbers = len(hr_paths)
    sum_psnr = 0
    max_psnr = 0
    min_psnr = 100
    sum_ssim = 0
    max_ssim = 0
    min_ssim = 1
    for hr_path in hr_paths:
        name, ext = os.path.splitext(os.path.basename(hr_path))
        img_name = os.path.basename(hr_path)
        print(img_name,end="       ")
        temp = str(name) + '_x4_SR.png'
        # print(temp)
        sr_path = os.path.join(SR_path, temp)
        #print(sr_path)
        # print(hr_path)
        # print(sr_path)
        img_Hr = util.imread_uint(hr_path, n_channels=n_channels)  # HR image, int8
        img_Sr = util.imread_uint(sr_path, n_channels=n_channels)  # HR image, int8
        # img_Hr = cv2.imread(hr_path)
        # img_Sr = cv2.imread(sr_path)
        psnr = util.calculate_psnr(img_Sr, img_Hr,)
        print(psnr)
        sum_psnr += psnr
        max_psnr = max(max_psnr,psnr)
        min_psnr = min(min_psnr, psnr)
        ssim = util.calculate_ssim(img_Sr, img_Hr,)
        # print(ssim)
        sum_ssim += ssim
        max_ssim = max(max_ssim,ssim)
        min_ssim = min(min_ssim, ssim)
    print('Average psnr = ', sum_psnr / numbers)
    print('min_psnr = ', min_psnr)
    print('Max_psnr = ', max_psnr)
    print('Average ssim = ', sum_ssim / numbers)
    print('min_ssim = ', min_ssim)
    print('Max_ssim = ', max_ssim)


if __name__ == '__main__':
    print('-------------------------compute psnr and ssim for evulate sr model---------------------------------')
    # evulate()
    evulate_diff_name()
### 回答1: 可以使用Python的第三方库来实现PSNRSSIM的计算PSNR: 可以使用pytube库中的'psnr'函数进行计算。 ``` from pytube import psnr original_image = # 原始图像 compressed_image = # 压缩后图像 psnr_value = psnr(original_image, compressed_image) print(psnr_value) ``` SSIM: 可以使用scikit-image库中的'ssim'函数进行计算。 ``` from skimage.metrics import structural_similarity as ssim original_image = # 原始图像 compressed_image = # 压缩后图像 ssim_value = ssim(original_image, compressed_image) print(ssim_value) ``` 注意: - 图像数据需要是numpy array格式 - PSNRSSIM的计算都依赖于图像的分辨率,因此请确保两个图像的分辨率相同 - PSNR值越高,说明图像质量越高; SSIM值越接近1,说明图像相似度越高. ### 回答2: PSNRSSIM是评估图像质量的经典评价指标。PSNR是峰值信噪比(Peak Signal-to-Noise Ratio),用于衡量图像的失真程度,而SSIM是结构相似性指标(Structural Similarity Index),用于衡量两幅图像的结构相似度。在python中,可以使用相关的包来实现这两个指标的计算PSNRpython实现: 对于两个图像x和y,它们的PSNR计算公式为: ``` PSNR = 10 * log10((MAX^2) / MSE) ``` 其中,MAX是像素值的最大值,通常为255;MSE是图像差分的均方误差,计算公式为: ``` MSE = 1 / (m * n) * sum_{i=0}^{m-1}sum_{j=0}^{n-1}(x(i,j) - y(i,j))^2 ``` 其中,m和n分别是图像的高和宽,x和y分别是两幅图像的像素矩阵。 python中可以使用numpy和PIL库来实现PSNR计算。具体代码如下: ``` import numpy as np from PIL import Image def psnr(img1, img2): img1 = np.array(img1).astype('float32') img2 = np.array(img2).astype('float32') mse = np.mean((img1 - img2) ** 2) if mse == 0: return float('inf') else: return 10 * np.log10((255 ** 2) / mse) #使用PIL库读取图像 img1 = Image.open('img1.jpg') img2 = Image.open('img2.jpg') #计算PSNRpsnr_value = psnr(img1, img2) print("PSNR值为:", psnr_value) ``` SSIM的python实现: 对于两幅图像x和y,它们的SSIM计算公式为: ``` SSIM(x, y) = (2 * mu_x * mu_y + C1) * (2 * sigma_xy + C2) / ((mu_x^2 + mu_y^2 + C1) * (sigma_x^2 + sigma_y^2 + C2)) ``` 其中,mu_x和mu_y分别是x和y的像素均值;sigma_x和sigma_y分别是x和y的像素标准差;sigma_xy是x和y的像素协方差;C1和C2是常数,其中C1=(K1 * L)^2,C2=(K2 * L)^2,L是像素值范围(通常为255),K1和K2是用于稳定分母的常数。 python中可以使用skimage和scipy库来实现SSIM的计算。具体代码如下: ``` from skimage.measure import compare_ssim from scipy.misc import imread #读取图像 img1 = imread('img1.jpg', mode='L') img2 = imread('img2.jpg', mode='L') #计算SSIM值 ssim_value = compare_ssim(img1, img2) print("SSIM值为:", ssim_value) ``` 注意:这里的图像读取方式使用的是scipy.misc库,如果读取不成功,可以考虑使用PIL库中的Image类读取。 ### 回答3: PSNR(Peak Signal to Noise Ratio)和SSIM(Structural Similarity Index)是图像处理中常用的评价指标。PSNR主要用于评估一幅图像和原始图像在像素级别上的相似程度,而SSIM可以测量两幅图像的结构相似性,对于一些噪声较多,图片质量较差的场合,与PSNR相比,SSIM更能反应图像的视觉质量。下面介绍如何使用python实现PSNRSSIM的计算。 一、PSNRPython实现 PSNR计算公式为:PSNR=10*log10(MAX^2/MSE),其中MAX为像素点最大值,一般为255,MSE是均方误差,即MSE=1/N * sum((I1-I2)^2),其中I1和I2分别为两幅图像中的像素值,N为像素出现的总次数。 - 实现方法 下面是使用Python实现PSNR的代码: ```python import cv2 import numpy as np def psnr(img1, img2): mse = np.mean((img1 - img2) ** 2) if mse == 0: return "Same Image" PIXEL_MAX = 255.0 return 10 * np.log10((PIXEL_MAX ** 2) / mse) ``` 其中img1和img2表示两幅图像的像素矩阵。 二、SSIM的Python实现 SSIM的计算公式分为三部分:亮度相似性(luminance)、对比度相似性(contrast)和结构相似性(structure),公式为:SSIM(x,y)=l(x,y)^alpha * c(x,y)^beta * s(x,y)^gamma,其中x和y分别表示两幅图像中对应像素点,^表示次方运算,alpha、beta、gamma分别为权重参数。 - 实现方法 下面是使用Python实现SSIM的代码: ```python import cv2 import numpy as np def ssim(img1, img2): C1 = (0.01 * 255) ** 2 C2 = (0.03 * 255) ** 2 img1 = img1.astype(np.float64) img2 = img2.astype(np.float64) kernel = cv2.getGaussianKernel(11, 1.5) window = np.outer(kernel, kernel.transpose()) mu1 = cv2.filter2D(img1, -1, window)[5:-5, 5:-5] mu2 = cv2.filter2D(img2, -1, window)[5:-5, 5:-5] mu1_sq = mu1 ** 2 mu2_sq = mu2 ** 2 mu1_mu2 = mu1 * mu2 sigma1_sq = cv2.filter2D(img1 ** 2, -1, window)[5:-5, 5:-5] - mu1_sq sigma2_sq = cv2.filter2D(img2 ** 2, -1, window)[5:-5, 5:-5] - mu2_sq sigma12 = cv2.filter2D(img1 * img2, -1, window)[5:-5, 5:-5] - mu1_mu2 ssim_map = ((2 * mu1_mu2 + C1) * (2 * sigma12 + C2)) / ((mu1_sq + mu2_sq + C1) * (sigma1_sq + sigma2_sq + C2)) return np.mean(ssim_map) ``` 其中img1和img2表示两幅图像的像素矩阵。该函数会返回两幅图像的SSIM值。 以上就是使用python实现PSNRSSIM的计算方法,可以根据实际需求选择一种或同时使用两种评价指标来衡量图像的质量。
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值