vp9_aq_complexity.c 5.9 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164
  1. /*
  2. * Copyright (c) 2014 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. #include <limits.h>
  11. #include <math.h>
  12. #include "vp9/encoder/vp9_aq_complexity.h"
  13. #include "vp9/encoder/vp9_aq_variance.h"
  14. #include "vp9/encoder/vp9_encodeframe.h"
  15. #include "vp9/common/vp9_seg_common.h"
  16. #include "vp9/encoder/vp9_segmentation.h"
  17. #define AQ_C_SEGMENTS 5
  18. #define DEFAULT_AQ2_SEG 3 // Neutral Q segment
  19. #define AQ_C_STRENGTHS 3
  20. static const double aq_c_q_adj_factor[AQ_C_STRENGTHS][AQ_C_SEGMENTS] =
  21. { {1.75, 1.25, 1.05, 1.00, 0.90},
  22. {2.00, 1.50, 1.15, 1.00, 0.85},
  23. {2.50, 1.75, 1.25, 1.00, 0.80} };
  24. static const double aq_c_transitions[AQ_C_STRENGTHS][AQ_C_SEGMENTS] =
  25. { {0.15, 0.30, 0.55, 2.00, 100.0},
  26. {0.20, 0.40, 0.65, 2.00, 100.0},
  27. {0.25, 0.50, 0.75, 2.00, 100.0} };
  28. static const double aq_c_var_thresholds[AQ_C_STRENGTHS][AQ_C_SEGMENTS] =
  29. { {-4.0, -3.0, -2.0, 100.00, 100.0},
  30. {-3.5, -2.5, -1.5, 100.00, 100.0},
  31. {-3.0, -2.0, -1.0, 100.00, 100.0} };
  32. #define DEFAULT_COMPLEXITY 64
  33. static int get_aq_c_strength(int q_index, vpx_bit_depth_t bit_depth) {
  34. // Approximate base quatizer (truncated to int)
  35. const int base_quant = vp9_ac_quant(q_index, 0, bit_depth) / 4;
  36. return (base_quant > 10) + (base_quant > 25);
  37. }
  38. void vp9_setup_in_frame_q_adj(VP9_COMP *cpi) {
  39. VP9_COMMON *const cm = &cpi->common;
  40. struct segmentation *const seg = &cm->seg;
  41. // Make SURE use of floating point in this function is safe.
  42. vp9_clear_system_state();
  43. if (cm->frame_type == KEY_FRAME ||
  44. cpi->refresh_alt_ref_frame ||
  45. (cpi->refresh_golden_frame && !cpi->rc.is_src_frame_alt_ref)) {
  46. int segment;
  47. const int aq_strength = get_aq_c_strength(cm->base_qindex, cm->bit_depth);
  48. // Clear down the segment map.
  49. memset(cpi->segmentation_map, DEFAULT_AQ2_SEG, cm->mi_rows * cm->mi_cols);
  50. vp9_clearall_segfeatures(seg);
  51. // Segmentation only makes sense if the target bits per SB is above a
  52. // threshold. Below this the overheads will usually outweigh any benefit.
  53. if (cpi->rc.sb64_target_rate < 256) {
  54. vp9_disable_segmentation(seg);
  55. return;
  56. }
  57. vp9_enable_segmentation(seg);
  58. // Select delta coding method.
  59. seg->abs_delta = SEGMENT_DELTADATA;
  60. // Default segment "Q" feature is disabled so it defaults to the baseline Q.
  61. vp9_disable_segfeature(seg, DEFAULT_AQ2_SEG, SEG_LVL_ALT_Q);
  62. // Use some of the segments for in frame Q adjustment.
  63. for (segment = 0; segment < AQ_C_SEGMENTS; ++segment) {
  64. int qindex_delta;
  65. if (segment == DEFAULT_AQ2_SEG)
  66. continue;
  67. qindex_delta =
  68. vp9_compute_qdelta_by_rate(&cpi->rc, cm->frame_type, cm->base_qindex,
  69. aq_c_q_adj_factor[aq_strength][segment],
  70. cm->bit_depth);
  71. // For AQ complexity mode, we dont allow Q0 in a segment if the base
  72. // Q is not 0. Q0 (lossless) implies 4x4 only and in AQ mode 2 a segment
  73. // Q delta is sometimes applied without going back around the rd loop.
  74. // This could lead to an illegal combination of partition size and q.
  75. if ((cm->base_qindex != 0) && ((cm->base_qindex + qindex_delta) == 0)) {
  76. qindex_delta = -cm->base_qindex + 1;
  77. }
  78. if ((cm->base_qindex + qindex_delta) > 0) {
  79. vp9_enable_segfeature(seg, segment, SEG_LVL_ALT_Q);
  80. vp9_set_segdata(seg, segment, SEG_LVL_ALT_Q, qindex_delta);
  81. }
  82. }
  83. }
  84. }
  85. #define DEFAULT_LV_THRESH 10.0
  86. #define MIN_DEFAULT_LV_THRESH 8.0
  87. #define VAR_STRENGTH_STEP 0.25
  88. // Select a segment for the current block.
  89. // The choice of segment for a block depends on the ratio of the projected
  90. // bits for the block vs a target average and its spatial complexity.
  91. void vp9_caq_select_segment(VP9_COMP *cpi, MACROBLOCK *mb, BLOCK_SIZE bs,
  92. int mi_row, int mi_col, int projected_rate) {
  93. VP9_COMMON *const cm = &cpi->common;
  94. const int mi_offset = mi_row * cm->mi_cols + mi_col;
  95. const int bw = num_8x8_blocks_wide_lookup[BLOCK_64X64];
  96. const int bh = num_8x8_blocks_high_lookup[BLOCK_64X64];
  97. const int xmis = MIN(cm->mi_cols - mi_col, num_8x8_blocks_wide_lookup[bs]);
  98. const int ymis = MIN(cm->mi_rows - mi_row, num_8x8_blocks_high_lookup[bs]);
  99. int x, y;
  100. int i;
  101. unsigned char segment;
  102. if (0) {
  103. segment = DEFAULT_AQ2_SEG;
  104. } else {
  105. // Rate depends on fraction of a SB64 in frame (xmis * ymis / bw * bh).
  106. // It is converted to bits * 256 units.
  107. const int target_rate = (cpi->rc.sb64_target_rate * xmis * ymis * 256) /
  108. (bw * bh);
  109. double logvar;
  110. double low_var_thresh;
  111. const int aq_strength = get_aq_c_strength(cm->base_qindex, cm->bit_depth);
  112. vp9_clear_system_state();
  113. low_var_thresh = (cpi->oxcf.pass == 2)
  114. ? MAX(cpi->twopass.mb_av_energy, MIN_DEFAULT_LV_THRESH)
  115. : DEFAULT_LV_THRESH;
  116. vp9_setup_src_planes(mb, cpi->Source, mi_row, mi_col);
  117. logvar = vp9_log_block_var(cpi, mb, bs);
  118. segment = AQ_C_SEGMENTS - 1; // Just in case no break out below.
  119. for (i = 0; i < AQ_C_SEGMENTS; ++i) {
  120. // Test rate against a threshold value and variance against a threshold.
  121. // Increasing segment number (higher variance and complexity) = higher Q.
  122. if ((projected_rate <
  123. target_rate * aq_c_transitions[aq_strength][i]) &&
  124. (logvar < (low_var_thresh + aq_c_var_thresholds[aq_strength][i]))) {
  125. segment = i;
  126. break;
  127. }
  128. }
  129. }
  130. // Fill in the entires in the segment map corresponding to this SB64.
  131. for (y = 0; y < ymis; y++) {
  132. for (x = 0; x < xmis; x++) {
  133. cpi->segmentation_map[mi_offset + y * cm->mi_cols + x] = segment;
  134. }
  135. }
  136. }