|
| 1 | +using SixLabors.ImageSharp.PixelFormats; |
| 2 | + |
| 3 | +namespace VerifyTests; |
| 4 | + |
| 5 | +public static class SsimComparer |
| 6 | +{ |
| 7 | + const double k1 = 0.01; |
| 8 | + const double k2 = 0.03; |
| 9 | + const double l = 255.0; |
| 10 | + const double c1 = k1 * l * k1 * l; |
| 11 | + const double c2 = k2 * l * k2 * l; |
| 12 | + |
| 13 | + public static double Calculate(Stream received, Stream verified) |
| 14 | + { |
| 15 | + using var imgReceived = Image.Load<Rgba32>(received); |
| 16 | + using var imgVerified = Image.Load<Rgba32>(verified); |
| 17 | + |
| 18 | + if (imgReceived.Width != imgVerified.Width || |
| 19 | + imgReceived.Height != imgVerified.Height) |
| 20 | + { |
| 21 | + return 0; |
| 22 | + } |
| 23 | + |
| 24 | + var width = imgReceived.Width; |
| 25 | + var height = imgReceived.Height; |
| 26 | + var pixelCount = (double) (width * height); |
| 27 | + |
| 28 | + double sumR1 = 0, sumG1 = 0, sumB1 = 0; |
| 29 | + double sumR2 = 0, sumG2 = 0, sumB2 = 0; |
| 30 | + double sumR1Sq = 0, sumG1Sq = 0, sumB1Sq = 0; |
| 31 | + double sumR2Sq = 0, sumG2Sq = 0, sumB2Sq = 0; |
| 32 | + double sumR12 = 0, sumG12 = 0, sumB12 = 0; |
| 33 | + |
| 34 | + imgReceived.ProcessPixelRows(imgVerified, (accessorReceived, accessorVerified) => |
| 35 | + { |
| 36 | + for (var y = 0; y < height; y++) |
| 37 | + { |
| 38 | + var row1 = accessorReceived.GetRowSpan(y); |
| 39 | + var row2 = accessorVerified.GetRowSpan(y); |
| 40 | + |
| 41 | + for (var x = 0; x < width; x++) |
| 42 | + { |
| 43 | + var p1 = row1[x]; |
| 44 | + var p2 = row2[x]; |
| 45 | + |
| 46 | + double r1 = p1.R, g1 = p1.G, b1 = p1.B; |
| 47 | + double r2 = p2.R, g2 = p2.G, b2 = p2.B; |
| 48 | + |
| 49 | + sumR1 += r1; |
| 50 | + sumG1 += g1; |
| 51 | + sumB1 += b1; |
| 52 | + sumR2 += r2; |
| 53 | + sumG2 += g2; |
| 54 | + sumB2 += b2; |
| 55 | + |
| 56 | + sumR1Sq += r1 * r1; |
| 57 | + sumG1Sq += g1 * g1; |
| 58 | + sumB1Sq += b1 * b1; |
| 59 | + sumR2Sq += r2 * r2; |
| 60 | + sumG2Sq += g2 * g2; |
| 61 | + sumB2Sq += b2 * b2; |
| 62 | + |
| 63 | + sumR12 += r1 * r2; |
| 64 | + sumG12 += g1 * g2; |
| 65 | + sumB12 += b1 * b2; |
| 66 | + } |
| 67 | + } |
| 68 | + }); |
| 69 | + |
| 70 | + var ssimR = CalculateChannel(pixelCount, sumR1, sumR2, sumR1Sq, sumR2Sq, sumR12); |
| 71 | + var ssimG = CalculateChannel(pixelCount, sumG1, sumG2, sumG1Sq, sumG2Sq, sumG12); |
| 72 | + var ssimB = CalculateChannel(pixelCount, sumB1, sumB2, sumB1Sq, sumB2Sq, sumB12); |
| 73 | + |
| 74 | + return (ssimR + ssimG + ssimB) / 3.0; |
| 75 | + } |
| 76 | + |
| 77 | + static double CalculateChannel(double n, double sum1, double sum2, double sum1Sq, double sum2Sq, double sum12) |
| 78 | + { |
| 79 | + var mu1 = sum1 / n; |
| 80 | + var mu2 = sum2 / n; |
| 81 | + var sigma1Sq = sum1Sq / n - mu1 * mu1; |
| 82 | + var sigma2Sq = sum2Sq / n - mu2 * mu2; |
| 83 | + var sigma12 = sum12 / n - mu1 * mu2; |
| 84 | + |
| 85 | + return (2 * mu1 * mu2 + c1) * (2 * sigma12 + c2) / |
| 86 | + ((mu1 * mu1 + mu2 * mu2 + c1) * (sigma1Sq + sigma2Sq + c2)); |
| 87 | + } |
| 88 | +} |
0 commit comments