| 42 | return _ssim(img1, img2, window, window_size, channel, size_average) |
| 43 | |
| 44 | def _ssim(img1, img2, window, window_size, channel, size_average=True): |
| 45 | mu1 = F.conv2d(img1, window, padding=window_size // 2, groups=channel) |
| 46 | mu2 = F.conv2d(img2, window, padding=window_size // 2, groups=channel) |
| 47 | |
| 48 | mu1_sq = mu1.pow(2) |
| 49 | mu2_sq = mu2.pow(2) |
| 50 | mu1_mu2 = mu1 * mu2 |
| 51 | |
| 52 | sigma1_sq = F.conv2d(img1 * img1, window, padding=window_size // 2, groups=channel) - mu1_sq |
| 53 | sigma2_sq = F.conv2d(img2 * img2, window, padding=window_size // 2, groups=channel) - mu2_sq |
| 54 | sigma12 = F.conv2d(img1 * img2, window, padding=window_size // 2, groups=channel) - mu1_mu2 |
| 55 | |
| 56 | C1 = 0.01 ** 2 |
| 57 | C2 = 0.03 ** 2 |
| 58 | |
| 59 | ssim_map = ((2 * mu1_mu2 + C1) * (2 * sigma12 + C2)) / ((mu1_sq + mu2_sq + C1) * (sigma1_sq + sigma2_sq + C2)) |
| 60 | |
| 61 | if size_average: |
| 62 | return ssim_map.mean() |
| 63 | else: |
| 64 | return ssim_map.mean(1).mean(1).mean(1) |
| 65 | |
| 66 | ms_ssim = MultiScaleStructuralSimilarityIndexMeasure(data_range=1.0) |
| 67 | |