1 /* 2 * Copyright (c) 2012 The WebM project authors. All Rights Reserved. 3 * 4 * Use of this source code is governed by a BSD-style license 5 * that can be found in the LICENSE file in the root of the source 6 * tree. An additional intellectual property rights grant can be found 7 * in the file PATENTS. All contributing project authors may 8 * be found in the AUTHORS file in the root of the source tree. 9 */ 10 11 #ifndef VPX_TEST_UTIL_H_ 12 #define VPX_TEST_UTIL_H_ 13 14 #include <stdio.h> 15 #include <math.h> 16 #include <tuple> 17 18 #include "gtest/gtest.h" 19 #include "vpx/vpx_image.h" 20 21 // Macros 22 #define GET_PARAM(k) std::get<k>(GetParam()) 23 compute_psnr(const vpx_image_t * img1,const vpx_image_t * img2)24inline double compute_psnr(const vpx_image_t *img1, const vpx_image_t *img2) { 25 assert((img1->fmt == img2->fmt) && (img1->d_w == img2->d_w) && 26 (img1->d_h == img2->d_h)); 27 28 const unsigned int width_y = img1->d_w; 29 const unsigned int height_y = img1->d_h; 30 unsigned int i, j; 31 32 int64_t sqrerr = 0; 33 for (i = 0; i < height_y; ++i) { 34 for (j = 0; j < width_y; ++j) { 35 int64_t d = img1->planes[VPX_PLANE_Y][i * img1->stride[VPX_PLANE_Y] + j] - 36 img2->planes[VPX_PLANE_Y][i * img2->stride[VPX_PLANE_Y] + j]; 37 sqrerr += d * d; 38 } 39 } 40 double mse = static_cast<double>(sqrerr) / (width_y * height_y); 41 double psnr = 100.0; 42 if (mse > 0.0) { 43 psnr = 10 * log10(255.0 * 255.0 / mse); 44 } 45 return psnr; 46 } 47 48 #endif // VPX_TEST_UTIL_H_ 49