Spaces:
Sleeping
Sleeping
# ''' | |
# https://github.com/One-sixth/ms_ssim_pytorch/blob/master/ssim.py | |
# ''' | |
# | |
# import torch | |
# import torch.jit | |
# import torch.nn.functional as F | |
# | |
# | |
# @torch.jit.script | |
# def create_window(window_size: int, sigma: float, channel: int): | |
# ''' | |
# Create 1-D gauss kernel | |
# :param window_size: the size of gauss kernel | |
# :param sigma: sigma of normal distribution | |
# :param channel: input channel | |
# :return: 1D kernel | |
# ''' | |
# coords = torch.arange(window_size, dtype=torch.float) | |
# coords -= window_size // 2 | |
# | |
# g = torch.exp(-(coords ** 2) / (2 * sigma ** 2)) | |
# g /= g.sum() | |
# | |
# g = g.reshape(1, 1, 1, -1).repeat(channel, 1, 1, 1) | |
# return g | |
# | |
# | |
# @torch.jit.script | |
# def _gaussian_filter(x, window_1d, use_padding: bool): | |
# ''' | |
# Blur input with 1-D kernel | |
# :param x: batch of tensors to be blured | |
# :param window_1d: 1-D gauss kernel | |
# :param use_padding: padding image before conv | |
# :return: blured tensors | |
# ''' | |
# C = x.shape[1] | |
# padding = 0 | |
# if use_padding: | |
# window_size = window_1d.shape[3] | |
# padding = window_size // 2 | |
# out = F.conv2d(x, window_1d, stride=1, padding=(0, padding), groups=C) | |
# out = F.conv2d(out, window_1d.transpose(2, 3), stride=1, padding=(padding, 0), groups=C) | |
# return out | |
# | |
# | |
# @torch.jit.script | |
# def ssim(X, Y, window, data_range: float, use_padding: bool = False): | |
# ''' | |
# Calculate ssim index for X and Y | |
# :param X: images [B, C, H, N_bins] | |
# :param Y: images [B, C, H, N_bins] | |
# :param window: 1-D gauss kernel | |
# :param data_range: value range of input images. (usually 1.0 or 255) | |
# :param use_padding: padding image before conv | |
# :return: | |
# ''' | |
# | |
# K1 = 0.01 | |
# K2 = 0.03 | |
# compensation = 1.0 | |
# | |
# C1 = (K1 * data_range) ** 2 | |
# C2 = (K2 * data_range) ** 2 | |
# | |
# mu1 = _gaussian_filter(X, window, use_padding) | |
# mu2 = _gaussian_filter(Y, window, use_padding) | |
# sigma1_sq = _gaussian_filter(X * X, window, use_padding) | |
# sigma2_sq = _gaussian_filter(Y * Y, window, use_padding) | |
# sigma12 = _gaussian_filter(X * Y, window, use_padding) | |
# | |
# mu1_sq = mu1.pow(2) | |
# mu2_sq = mu2.pow(2) | |
# mu1_mu2 = mu1 * mu2 | |
# | |
# sigma1_sq = compensation * (sigma1_sq - mu1_sq) | |
# sigma2_sq = compensation * (sigma2_sq - mu2_sq) | |
# sigma12 = compensation * (sigma12 - mu1_mu2) | |
# | |
# cs_map = (2 * sigma12 + C2) / (sigma1_sq + sigma2_sq + C2) | |
# # Fixed the issue that the negative value of cs_map caused ms_ssim to output Nan. | |
# cs_map = cs_map.clamp_min(0.) | |
# ssim_map = ((2 * mu1_mu2 + C1) / (mu1_sq + mu2_sq + C1)) * cs_map | |
# | |
# ssim_val = ssim_map.mean(dim=(1, 2, 3)) # reduce along CHW | |
# cs = cs_map.mean(dim=(1, 2, 3)) | |
# | |
# return ssim_val, cs | |
# | |
# | |
# @torch.jit.script | |
# def ms_ssim(X, Y, window, data_range: float, weights, use_padding: bool = False, eps: float = 1e-8): | |
# ''' | |
# interface of ms-ssim | |
# :param X: a batch of images, (N,C,H,W) | |
# :param Y: a batch of images, (N,C,H,W) | |
# :param window: 1-D gauss kernel | |
# :param data_range: value range of input images. (usually 1.0 or 255) | |
# :param weights: weights for different levels | |
# :param use_padding: padding image before conv | |
# :param eps: use for avoid grad nan. | |
# :return: | |
# ''' | |
# levels = weights.shape[0] | |
# cs_vals = [] | |
# ssim_vals = [] | |
# for _ in range(levels): | |
# ssim_val, cs = ssim(X, Y, window=window, data_range=data_range, use_padding=use_padding) | |
# # Use for fix a issue. When c = a ** b and a is 0, c.backward() will cause the a.grad become inf. | |
# ssim_val = ssim_val.clamp_min(eps) | |
# cs = cs.clamp_min(eps) | |
# cs_vals.append(cs) | |
# | |
# ssim_vals.append(ssim_val) | |
# padding = (X.shape[2] % 2, X.shape[3] % 2) | |
# X = F.avg_pool2d(X, kernel_size=2, stride=2, padding=padding) | |
# Y = F.avg_pool2d(Y, kernel_size=2, stride=2, padding=padding) | |
# | |
# cs_vals = torch.stack(cs_vals, dim=0) | |
# ms_ssim_val = torch.prod((cs_vals[:-1] ** weights[:-1].unsqueeze(1)) * (ssim_vals[-1] ** weights[-1]), dim=0) | |
# return ms_ssim_val | |
# | |
# | |
# class SSIM(torch.jit.ScriptModule): | |
# __constants__ = ['data_range', 'use_padding'] | |
# | |
# def __init__(self, window_size=11, window_sigma=1.5, data_range=255., channel=3, use_padding=False): | |
# ''' | |
# :param window_size: the size of gauss kernel | |
# :param window_sigma: sigma of normal distribution | |
# :param data_range: value range of input images. (usually 1.0 or 255) | |
# :param channel: input channels (default: 3) | |
# :param use_padding: padding image before conv | |
# ''' | |
# super().__init__() | |
# assert window_size % 2 == 1, 'Window size must be odd.' | |
# window = create_window(window_size, window_sigma, channel) | |
# self.register_buffer('window', window) | |
# self.data_range = data_range | |
# self.use_padding = use_padding | |
# | |
# @torch.jit.script_method | |
# def forward(self, X, Y): | |
# r = ssim(X, Y, window=self.window, data_range=self.data_range, use_padding=self.use_padding) | |
# return r[0] | |
# | |
# | |
# class MS_SSIM(torch.jit.ScriptModule): | |
# __constants__ = ['data_range', 'use_padding', 'eps'] | |
# | |
# def __init__(self, window_size=11, window_sigma=1.5, data_range=255., channel=3, use_padding=False, weights=None, | |
# levels=None, eps=1e-8): | |
# ''' | |
# class for ms-ssim | |
# :param window_size: the size of gauss kernel | |
# :param window_sigma: sigma of normal distribution | |
# :param data_range: value range of input images. (usually 1.0 or 255) | |
# :param channel: input channels | |
# :param use_padding: padding image before conv | |
# :param weights: weights for different levels. (default [0.0448, 0.2856, 0.3001, 0.2363, 0.1333]) | |
# :param levels: number of downsampling | |
# :param eps: Use for fix a issue. When c = a ** b and a is 0, c.backward() will cause the a.grad become inf. | |
# ''' | |
# super().__init__() | |
# assert window_size % 2 == 1, 'Window size must be odd.' | |
# self.data_range = data_range | |
# self.use_padding = use_padding | |
# self.eps = eps | |
# | |
# window = create_window(window_size, window_sigma, channel) | |
# self.register_buffer('window', window) | |
# | |
# if weights is None: | |
# weights = [0.0448, 0.2856, 0.3001, 0.2363, 0.1333] | |
# weights = torch.tensor(weights, dtype=torch.float) | |
# | |
# if levels is not None: | |
# weights = weights[:levels] | |
# weights = weights / weights.sum() | |
# | |
# self.register_buffer('weights', weights) | |
# | |
# @torch.jit.script_method | |
# def forward(self, X, Y): | |
# return ms_ssim(X, Y, window=self.window, data_range=self.data_range, weights=self.weights, | |
# use_padding=self.use_padding, eps=self.eps) | |
# | |
# | |
# if __name__ == '__main__': | |
# print('Simple Test') | |
# im = torch.randint(0, 255, (5, 3, 256, 256), dtype=torch.float, device='cuda') | |
# img1 = im / 255 | |
# img2 = img1 * 0.5 | |
# | |
# losser = SSIM(data_range=1.).cuda() | |
# loss = losser(img1, img2).mean() | |
# | |
# losser2 = MS_SSIM(data_range=1.).cuda() | |
# loss2 = losser2(img1, img2).mean() | |
# | |
# print(loss.item()) | |
# print(loss2.item()) | |
# | |
# if __name__ == '__main__': | |
# print('Training Test') | |
# import cv2 | |
# import torch.optim | |
# import numpy as np | |
# import imageio | |
# import time | |
# | |
# out_test_video = False | |
# # 最好不要直接输出gif图,会非常大,最好先输出mkv文件后用ffmpeg转换到GIF | |
# video_use_gif = False | |
# | |
# im = cv2.imread('test_img1.jpg', 1) | |
# t_im = torch.from_numpy(im).cuda().permute(2, 0, 1).float()[None] / 255. | |
# | |
# if out_test_video: | |
# if video_use_gif: | |
# fps = 0.5 | |
# out_wh = (im.shape[1] // 2, im.shape[0] // 2) | |
# suffix = '.gif' | |
# else: | |
# fps = 5 | |
# out_wh = (im.shape[1], im.shape[0]) | |
# suffix = '.mkv' | |
# video_last_time = time.perf_counter() | |
# video = imageio.get_writer('ssim_test' + suffix, fps=fps) | |
# | |
# # 测试ssim | |
# print('Training SSIM') | |
# rand_im = torch.randint_like(t_im, 0, 255, dtype=torch.float32) / 255. | |
# rand_im.requires_grad = True | |
# optim = torch.optim.Adam([rand_im], 0.003, eps=1e-8) | |
# losser = SSIM(data_range=1., channel=t_im.shape[1]).cuda() | |
# ssim_score = 0 | |
# while ssim_score < 0.999: | |
# optim.zero_grad() | |
# loss = losser(rand_im, t_im) | |
# (-loss).sum().backward() | |
# ssim_score = loss.item() | |
# optim.step() | |
# r_im = np.transpose(rand_im.detach().cpu().numpy().clip(0, 1) * 255, [0, 2, 3, 1]).astype(np.uint8)[0] | |
# r_im = cv2.putText(r_im, 'ssim %f' % ssim_score, (10, 30), cv2.FONT_HERSHEY_PLAIN, 2, (255, 0, 0), 2) | |
# | |
# if out_test_video: | |
# if time.perf_counter() - video_last_time > 1. / fps: | |
# video_last_time = time.perf_counter() | |
# out_frame = cv2.cvtColor(r_im, cv2.COLOR_BGR2RGB) | |
# out_frame = cv2.resize(out_frame, out_wh, interpolation=cv2.INTER_AREA) | |
# if isinstance(out_frame, cv2.UMat): | |
# out_frame = out_frame.get() | |
# video.append_data(out_frame) | |
# | |
# cv2.imshow('ssim', r_im) | |
# cv2.setWindowTitle('ssim', 'ssim %f' % ssim_score) | |
# cv2.waitKey(1) | |
# | |
# if out_test_video: | |
# video.close() | |
# | |
# # 测试ms_ssim | |
# if out_test_video: | |
# if video_use_gif: | |
# fps = 0.5 | |
# out_wh = (im.shape[1] // 2, im.shape[0] // 2) | |
# suffix = '.gif' | |
# else: | |
# fps = 5 | |
# out_wh = (im.shape[1], im.shape[0]) | |
# suffix = '.mkv' | |
# video_last_time = time.perf_counter() | |
# video = imageio.get_writer('ms_ssim_test' + suffix, fps=fps) | |
# | |
# print('Training MS_SSIM') | |
# rand_im = torch.randint_like(t_im, 0, 255, dtype=torch.float32) / 255. | |
# rand_im.requires_grad = True | |
# optim = torch.optim.Adam([rand_im], 0.003, eps=1e-8) | |
# losser = MS_SSIM(data_range=1., channel=t_im.shape[1]).cuda() | |
# ssim_score = 0 | |
# while ssim_score < 0.999: | |
# optim.zero_grad() | |
# loss = losser(rand_im, t_im) | |
# (-loss).sum().backward() | |
# ssim_score = loss.item() | |
# optim.step() | |
# r_im = np.transpose(rand_im.detach().cpu().numpy().clip(0, 1) * 255, [0, 2, 3, 1]).astype(np.uint8)[0] | |
# r_im = cv2.putText(r_im, 'ms_ssim %f' % ssim_score, (10, 30), cv2.FONT_HERSHEY_PLAIN, 2, (255, 0, 0), 2) | |
# | |
# if out_test_video: | |
# if time.perf_counter() - video_last_time > 1. / fps: | |
# video_last_time = time.perf_counter() | |
# out_frame = cv2.cvtColor(r_im, cv2.COLOR_BGR2RGB) | |
# out_frame = cv2.resize(out_frame, out_wh, interpolation=cv2.INTER_AREA) | |
# if isinstance(out_frame, cv2.UMat): | |
# out_frame = out_frame.get() | |
# video.append_data(out_frame) | |
# | |
# cv2.imshow('ms_ssim', r_im) | |
# cv2.setWindowTitle('ms_ssim', 'ms_ssim %f' % ssim_score) | |
# cv2.waitKey(1) | |
# | |
# if out_test_video: | |
# video.close() | |
""" | |
Adapted from https://github.com/Po-Hsun-Su/pytorch-ssim | |
""" | |
import torch | |
import torch.nn.functional as F | |
from torch.autograd import Variable | |
import numpy as np | |
from math import exp | |
def gaussian(window_size, sigma): | |
gauss = torch.Tensor([exp(-(x - window_size // 2) ** 2 / float(2 * sigma ** 2)) for x in range(window_size)]) | |
return gauss / gauss.sum() | |
def create_window(window_size, channel): | |
_1D_window = gaussian(window_size, 1.5).unsqueeze(1) | |
_2D_window = _1D_window.mm(_1D_window.t()).float().unsqueeze(0).unsqueeze(0) | |
window = Variable(_2D_window.expand(channel, 1, window_size, window_size).contiguous()) | |
return window | |
def _ssim(img1, img2, window, window_size, channel, size_average=True): | |
mu1 = F.conv2d(img1, window, padding=window_size // 2, groups=channel) | |
mu2 = F.conv2d(img2, window, padding=window_size // 2, groups=channel) | |
mu1_sq = mu1.pow(2) | |
mu2_sq = mu2.pow(2) | |
mu1_mu2 = mu1 * mu2 | |
sigma1_sq = F.conv2d(img1 * img1, window, padding=window_size // 2, groups=channel) - mu1_sq | |
sigma2_sq = F.conv2d(img2 * img2, window, padding=window_size // 2, groups=channel) - mu2_sq | |
sigma12 = F.conv2d(img1 * img2, window, padding=window_size // 2, groups=channel) - mu1_mu2 | |
C1 = 0.01 ** 2 | |
C2 = 0.03 ** 2 | |
ssim_map = ((2 * mu1_mu2 + C1) * (2 * sigma12 + C2)) / ((mu1_sq + mu2_sq + C1) * (sigma1_sq + sigma2_sq + C2)) | |
if size_average: | |
return ssim_map.mean() | |
else: | |
return ssim_map.mean(1) | |
class SSIM(torch.nn.Module): | |
def __init__(self, window_size=11, size_average=True): | |
super(SSIM, self).__init__() | |
self.window_size = window_size | |
self.size_average = size_average | |
self.channel = 1 | |
self.window = create_window(window_size, self.channel) | |
def forward(self, img1, img2): | |
(_, channel, _, _) = img1.size() | |
if channel == self.channel and self.window.data.type() == img1.data.type(): | |
window = self.window | |
else: | |
window = create_window(self.window_size, channel) | |
if img1.is_cuda: | |
window = window.cuda(img1.get_device()) | |
window = window.type_as(img1) | |
self.window = window | |
self.channel = channel | |
return _ssim(img1, img2, window, self.window_size, channel, self.size_average) | |
window = None | |
def ssim(img1, img2, window_size=11, size_average=True): | |
(_, channel, _, _) = img1.size() | |
global window | |
if window is None: | |
window = create_window(window_size, channel) | |
if img1.is_cuda: | |
window = window.cuda(img1.get_device()) | |
window = window.type_as(img1) | |
return _ssim(img1, img2, window, window_size, channel, size_average) | |