util.h 1.4 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546
  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. #ifndef TEST_UTIL_H_
  11. #define TEST_UTIL_H_
  12. #include <stdio.h>
  13. #include <math.h>
  14. #include "third_party/googletest/src/include/gtest/gtest.h"
  15. #include "vpx/vpx_image.h"
  16. // Macros
  17. #define GET_PARAM(k) std::tr1::get<k>(GetParam())
  18. inline double compute_psnr(const vpx_image_t *img1, const vpx_image_t *img2) {
  19. assert((img1->fmt == img2->fmt) && (img1->d_w == img2->d_w) &&
  20. (img1->d_h == img2->d_h));
  21. const unsigned int width_y = img1->d_w;
  22. const unsigned int height_y = img1->d_h;
  23. unsigned int i, j;
  24. int64_t sqrerr = 0;
  25. for (i = 0; i < height_y; ++i) {
  26. for (j = 0; j < width_y; ++j) {
  27. int64_t d = img1->planes[VPX_PLANE_Y][i * img1->stride[VPX_PLANE_Y] + j] -
  28. img2->planes[VPX_PLANE_Y][i * img2->stride[VPX_PLANE_Y] + j];
  29. sqrerr += d * d;
  30. }
  31. }
  32. double mse = static_cast<double>(sqrerr) / (width_y * height_y);
  33. double psnr = 100.0;
  34. if (mse > 0.0) {
  35. psnr = 10 * log10(255.0 * 255.0 / mse);
  36. }
  37. return psnr;
  38. }
  39. #endif // TEST_UTIL_H_