1470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com/* 2470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com * Copyright (c) 2011 The WebRTC project authors. All Rights Reserved. 3470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com * 4470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com * Use of this source code is governed by a BSD-style license 5470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com * that can be found in the LICENSE file in the root of the source 6470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com * tree. An additional intellectual property rights grant can be found 7470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com * in the file PATENTS. All contributing project authors may 8470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com * be found in the AUTHORS file in the root of the source tree. 9470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com */ 10470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 110f59a88b32d95430bf12bd8c4d9b2963e018ead2Henrik Kjellander#include "webrtc/modules/video_processing/deflickering.h" 12470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 13470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com#include <math.h> 14470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com#include <stdlib.h> 15470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 16854e84c7fb61972bd8e39c1482e5f3e31f796b5fpbos#include "webrtc/base/logging.h" 176f3d8fcfc08686a9e8852ad225d54cca6e197875pbos@webrtc.org#include "webrtc/common_audio/signal_processing/include/signal_processing_library.h" 1898f53510b222f71fdd8b799b2f33737ceeb28c61Henrik Kjellander#include "webrtc/system_wrappers/include/sort.h" 19470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 20470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.comnamespace webrtc { 21470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 22470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com// Detection constants 23b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org// (Q4) Maximum allowed deviation for detection. 24b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.orgenum { kFrequencyDeviation = 39 }; 25b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org// (Q4) Minimum frequency that can be detected. 26b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.orgenum { kMinFrequencyToDetect = 32 }; 27b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org// Number of flickers before we accept detection 28b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.orgenum { kNumFlickerBeforeDetect = 2 }; 29b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.orgenum { kmean_valueScaling = 4 }; // (Q4) In power of 2 30b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org// Dead-zone region in terms of pixel values 31b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.orgenum { kZeroCrossingDeadzone = 10 }; 32b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org// Deflickering constants. 33470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com// Compute the quantiles over 1 / DownsamplingFactor of the image. 34470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.comenum { kDownsamplingFactor = 8 }; 35470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.comenum { kLog2OfDownsamplingFactor = 3 }; 36470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 37470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com// To generate in Matlab: 38b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org// >> probUW16 = round(2^11 * 39b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org// [0.05,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,0.95,0.97]); 40470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com// >> fprintf('%d, ', probUW16) 41b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org// Resolution reduced to avoid overflow when multiplying with the 42b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org// (potentially) large number of pixels. 4399ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodmanconst uint16_t VPMDeflickering::prob_uw16_[kNumProbs] = { 4499ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman 102, 205, 410, 614, 819, 1024, 4599ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman 1229, 1434, 1638, 1843, 1946, 1987}; // <Q11> 46470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 47470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com// To generate in Matlab: 48470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com// >> numQuants = 14; maxOnlyLength = 5; 49b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org// >> weightUW16 = round(2^15 * 50b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org// [linspace(0.5, 1.0, numQuants - maxOnlyLength)]); 51470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com// >> fprintf('%d, %d,\n ', weightUW16); 5299ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodmanconst uint16_t VPMDeflickering::weight_uw16_[kNumQuants - kMaxOnlyLength] = { 5399ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman 16384, 18432, 20480, 22528, 24576, 26624, 28672, 30720, 32768}; // <Q15> 54470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 554161715e3f7e744bc9ef3d3ae437da1e8e4de38dtommi@webrtc.orgVPMDeflickering::VPMDeflickering() { 56b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org Reset(); 57470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com} 58470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 59b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.orgVPMDeflickering::~VPMDeflickering() {} 60470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 61b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.orgvoid VPMDeflickering::Reset() { 62b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org mean_buffer_length_ = 0; 63b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org detection_state_ = 0; 64b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org frame_rate_ = 0; 65b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 66b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org memset(mean_buffer_, 0, sizeof(int32_t) * kMeanBufferLength); 67b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org memset(timestamp_buffer_, 0, sizeof(int32_t) * kMeanBufferLength); 68b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 69b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // Initialize the history with a uniformly distributed histogram. 70b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org quant_hist_uw8_[0][0] = 0; 71b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org quant_hist_uw8_[0][kNumQuants - 1] = 255; 72b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org for (int32_t i = 0; i < kNumProbs; i++) { 730d394f3609e039db3c4168bf3f8273b3a1e02acbbjornv@webrtc.org // Unsigned round. <Q0> 7499ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman quant_hist_uw8_[0][i + 1] = 7599ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman static_cast<uint8_t>((prob_uw16_[i] * 255 + (1 << 10)) >> 11); 76b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 77b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 78b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org for (int32_t i = 1; i < kFrameHistory_size; i++) { 79b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org memcpy(quant_hist_uw8_[i], quant_hist_uw8_[0], 80b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org sizeof(uint8_t) * kNumQuants); 81b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 82b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org} 83470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 8499ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodmanint32_t VPMDeflickering::ProcessFrame(VideoFrame* frame, 8599ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman VideoProcessing::FrameStats* stats) { 86b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org assert(frame); 87b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org uint32_t frame_memory; 88b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org uint8_t quant_uw8[kNumQuants]; 89b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org uint8_t maxquant_uw8[kNumQuants]; 90b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org uint8_t minquant_uw8[kNumQuants]; 91b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org uint16_t target_quant_uw16[kNumQuants]; 92b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org uint16_t increment_uw16; 93b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org uint8_t map_uw8[256]; 94b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 95b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org uint16_t tmp_uw16; 96b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org uint32_t tmp_uw32; 97b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org int width = frame->width(); 98b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org int height = frame->height(); 99b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 100b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org if (frame->IsZeroSize()) { 101b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org return VPM_GENERAL_ERROR; 102b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 103b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 104b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // Stricter height check due to subsampling size calculation below. 105b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org if (height < 2) { 1062a770828d8594584b1275db06ad43fbe0d081279asapersson@webrtc.org LOG(LS_ERROR) << "Invalid frame size."; 107b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org return VPM_GENERAL_ERROR; 108b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 109b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 110a8565425bc6b1d3d7b5a0b0e0a49bc4e57c64f21mflodman if (!VideoProcessing::ValidFrameStats(*stats)) { 111b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org return VPM_GENERAL_ERROR; 112b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 113b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 11499ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman if (PreDetection(frame->timestamp(), *stats) == -1) 11599ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman return VPM_GENERAL_ERROR; 116b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 117b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // Flicker detection 118b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org int32_t det_flicker = DetectFlicker(); 119b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org if (det_flicker < 0) { 120b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org return VPM_GENERAL_ERROR; 121b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } else if (det_flicker != 1) { 122b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org return 0; 123b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 124b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 125b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // Size of luminance component. 126b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org const uint32_t y_size = height * width; 127b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 12899ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman const uint32_t y_sub_size = 12999ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman width * (((height - 1) >> kLog2OfDownsamplingFactor) + 1); 130b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org uint8_t* y_sorted = new uint8_t[y_sub_size]; 131b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org uint32_t sort_row_idx = 0; 132b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org for (int i = 0; i < height; i += kDownsamplingFactor) { 13399ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman memcpy(y_sorted + sort_row_idx * width, frame->buffer(kYPlane) + i * width, 13499ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman width); 135b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org sort_row_idx++; 136b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 137b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 138b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org webrtc::Sort(y_sorted, y_sub_size, webrtc::TYPE_UWord8); 139b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 140b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org uint32_t prob_idx_uw32 = 0; 141b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org quant_uw8[0] = 0; 142b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org quant_uw8[kNumQuants - 1] = 255; 143b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 144b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // Ensure we won't get an overflow below. 145b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // In practice, the number of subsampled pixels will not become this large. 146b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org if (y_sub_size > (1 << 21) - 1) { 1472a770828d8594584b1275db06ad43fbe0d081279asapersson@webrtc.org LOG(LS_ERROR) << "Subsampled number of pixels too large."; 148b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org return -1; 149b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 150b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 151b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org for (int32_t i = 0; i < kNumProbs; i++) { 152b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // <Q0>. 153b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org prob_idx_uw32 = WEBRTC_SPL_UMUL_32_16(y_sub_size, prob_uw16_[i]) >> 11; 154b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org quant_uw8[i + 1] = y_sorted[prob_idx_uw32]; 155b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 156b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 15799ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman delete[] y_sorted; 158b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org y_sorted = NULL; 159b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 160b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // Shift history for new frame. 161b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org memmove(quant_hist_uw8_[1], quant_hist_uw8_[0], 16299ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman (kFrameHistory_size - 1) * kNumQuants * sizeof(uint8_t)); 163b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // Store current frame in history. 164b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org memcpy(quant_hist_uw8_[0], quant_uw8, kNumQuants * sizeof(uint8_t)); 165b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 166b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // We use a frame memory equal to the ceiling of half the frame rate to 167b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // ensure we capture an entire period of flicker. 168b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org frame_memory = (frame_rate_ + (1 << 5)) >> 5; // Unsigned ceiling. <Q0> 169b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // frame_rate_ in Q4. 170b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org if (frame_memory > kFrameHistory_size) { 171b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org frame_memory = kFrameHistory_size; 172b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 173b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 174b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // Get maximum and minimum. 175b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org for (int32_t i = 0; i < kNumQuants; i++) { 176b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org maxquant_uw8[i] = 0; 177b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org minquant_uw8[i] = 255; 178b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org for (uint32_t j = 0; j < frame_memory; j++) { 179b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org if (quant_hist_uw8_[j][i] > maxquant_uw8[i]) { 180b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org maxquant_uw8[i] = quant_hist_uw8_[j][i]; 181b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 182b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 183b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org if (quant_hist_uw8_[j][i] < minquant_uw8[i]) { 184b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org minquant_uw8[i] = quant_hist_uw8_[j][i]; 185b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 186470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com } 187b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 188b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 189b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // Get target quantiles. 190b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org for (int32_t i = 0; i < kNumQuants - kMaxOnlyLength; i++) { 1910d394f3609e039db3c4168bf3f8273b3a1e02acbbjornv@webrtc.org // target = w * maxquant_uw8 + (1 - w) * minquant_uw8 1920d394f3609e039db3c4168bf3f8273b3a1e02acbbjornv@webrtc.org // Weights w = |weight_uw16_| are in Q15, hence the final output has to be 1930d394f3609e039db3c4168bf3f8273b3a1e02acbbjornv@webrtc.org // right shifted by 8 to end up in Q7. 19499ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman target_quant_uw16[i] = static_cast<uint16_t>( 19599ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman (weight_uw16_[i] * maxquant_uw8[i] + 19699ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman ((1 << 15) - weight_uw16_[i]) * minquant_uw8[i]) >> 19799ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman 8); // <Q7> 198b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 199b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 200b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org for (int32_t i = kNumQuants - kMaxOnlyLength; i < kNumQuants; i++) { 201b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org target_quant_uw16[i] = ((uint16_t)maxquant_uw8[i]) << 7; 202b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 203b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 204b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // Compute the map from input to output pixels. 205b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org uint16_t mapUW16; // <Q7> 206b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org for (int32_t i = 1; i < kNumQuants; i++) { 207b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // As quant and targetQuant are limited to UWord8, it's safe to use Q7 here. 20899ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman tmp_uw32 = 20999ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman static_cast<uint32_t>(target_quant_uw16[i] - target_quant_uw16[i - 1]); 210b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org tmp_uw16 = static_cast<uint16_t>(quant_uw8[i] - quant_uw8[i - 1]); // <Q0> 211b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 212b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org if (tmp_uw16 > 0) { 21399ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman increment_uw16 = 21499ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman static_cast<uint16_t>(WebRtcSpl_DivU32U16(tmp_uw32, 21599ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman tmp_uw16)); // <Q7> 216b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } else { 217b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // The value is irrelevant; the loop below will only iterate once. 218b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org increment_uw16 = 0; 219470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com } 220470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 221b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org mapUW16 = target_quant_uw16[i - 1]; 222b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org for (uint32_t j = quant_uw8[i - 1]; j < (uint32_t)(quant_uw8[i] + 1); j++) { 223b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // Unsigned round. <Q0> 224b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org map_uw8[j] = (uint8_t)((mapUW16 + (1 << 6)) >> 7); 225b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org mapUW16 += increment_uw16; 226470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com } 227b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 228470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 229b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // Map to the output frame. 230b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org uint8_t* buffer = frame->buffer(kYPlane); 231b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org for (uint32_t i = 0; i < y_size; i++) { 232b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org buffer[i] = map_uw8[buffer[i]]; 233b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 234470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 235b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // Frame was altered, so reset stats. 236a8565425bc6b1d3d7b5a0b0e0a49bc4e57c64f21mflodman VideoProcessing::ClearFrameStats(stats); 237470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 238b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org return VPM_OK; 239470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com} 240470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 241470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com/** 242b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org Performs some pre-detection operations. Must be called before 243470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com DetectFlicker(). 244470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 245470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com \param[in] timestamp Timestamp of the current frame. 246470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com \param[in] stats Statistics of the current frame. 247b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 248470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com \return 0: Success\n 249470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 2: Detection not possible due to flickering frequency too close to 250470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com zero.\n 251470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com -1: Error 252470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com*/ 25399ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodmanint32_t VPMDeflickering::PreDetection( 25499ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman const uint32_t timestamp, 25599ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman const VideoProcessing::FrameStats& stats) { 256b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org int32_t mean_val; // Mean value of frame (Q4) 257b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org uint32_t frame_rate = 0; 258b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org int32_t meanBufferLength; // Temp variable. 259b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 260b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org mean_val = ((stats.sum << kmean_valueScaling) / stats.num_pixels); 261b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // Update mean value buffer. 262b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // This should be done even though we might end up in an unreliable detection. 263b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org memmove(mean_buffer_ + 1, mean_buffer_, 26499ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman (kMeanBufferLength - 1) * sizeof(int32_t)); 265b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org mean_buffer_[0] = mean_val; 266b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 267b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // Update timestamp buffer. 268b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // This should be done even though we might end up in an unreliable detection. 26999ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman memmove(timestamp_buffer_ + 1, timestamp_buffer_, 27099ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman (kMeanBufferLength - 1) * sizeof(uint32_t)); 271b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org timestamp_buffer_[0] = timestamp; 272b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 27399ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman /* Compute current frame rate (Q4) */ 274b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org if (timestamp_buffer_[kMeanBufferLength - 1] != 0) { 275b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org frame_rate = ((90000 << 4) * (kMeanBufferLength - 1)); 276b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org frame_rate /= 277b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org (timestamp_buffer_[0] - timestamp_buffer_[kMeanBufferLength - 1]); 278b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } else if (timestamp_buffer_[1] != 0) { 279b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org frame_rate = (90000 << 4) / (timestamp_buffer_[0] - timestamp_buffer_[1]); 280b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 281b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 282b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org /* Determine required size of mean value buffer (mean_buffer_length_) */ 283b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org if (frame_rate == 0) { 284b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org meanBufferLength = 1; 285b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } else { 286b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org meanBufferLength = 287b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org (kNumFlickerBeforeDetect * frame_rate) / kMinFrequencyToDetect; 288b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 289b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org /* Sanity check of buffer length */ 290b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org if (meanBufferLength >= kMeanBufferLength) { 291b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org /* Too long buffer. The flickering frequency is too close to zero, which 292b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org * makes the estimation unreliable. 293470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com */ 294b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org mean_buffer_length_ = 0; 295b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org return 2; 296b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 297b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org mean_buffer_length_ = meanBufferLength; 298b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 299b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org if ((timestamp_buffer_[mean_buffer_length_ - 1] != 0) && 300b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org (mean_buffer_length_ != 1)) { 301b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org frame_rate = ((90000 << 4) * (mean_buffer_length_ - 1)); 302b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org frame_rate /= 303b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org (timestamp_buffer_[0] - timestamp_buffer_[mean_buffer_length_ - 1]); 304b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } else if (timestamp_buffer_[1] != 0) { 305b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org frame_rate = (90000 << 4) / (timestamp_buffer_[0] - timestamp_buffer_[1]); 306b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 307b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org frame_rate_ = frame_rate; 308b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 309b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org return VPM_OK; 310470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com} 311470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 312470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com/** 313b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org This function detects flicker in the video stream. As a side effect the 314b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org mean value buffer is updated with the new mean value. 315b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 316470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com \return 0: No flickering detected\n 317470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 1: Flickering detected\n 318470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 2: Detection not possible due to unreliable frequency interval 319470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com -1: Error 320470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com*/ 321b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.orgint32_t VPMDeflickering::DetectFlicker() { 32299ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman uint32_t i; 32399ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman int32_t freqEst; // (Q4) Frequency estimate to base detection upon 32499ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman int32_t ret_val = -1; 325b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 326b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org /* Sanity check for mean_buffer_length_ */ 327b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org if (mean_buffer_length_ < 2) { 328b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org /* Not possible to estimate frequency */ 32999ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman return 2; 330b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 331b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // Count zero crossings with a dead zone to be robust against noise. If the 332b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // noise std is 2 pixel this corresponds to about 95% confidence interval. 333b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org int32_t deadzone = (kZeroCrossingDeadzone << kmean_valueScaling); // Q4 334b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org int32_t meanOfBuffer = 0; // Mean value of mean value buffer. 33599ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman int32_t numZeros = 0; // Number of zeros that cross the dead-zone. 33699ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman int32_t cntState = 0; // State variable for zero crossing regions. 33799ab9447d1b48043adf6bfdf72cf820bbda3ee3fmflodman int32_t cntStateOld = 0; // Previous state for zero crossing regions. 338b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 339b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org for (i = 0; i < mean_buffer_length_; i++) { 340b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org meanOfBuffer += mean_buffer_[i]; 341b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 342b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org meanOfBuffer += (mean_buffer_length_ >> 1); // Rounding, not truncation. 343b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org meanOfBuffer /= mean_buffer_length_; 344b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 345b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // Count zero crossings. 346b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org cntStateOld = (mean_buffer_[0] >= (meanOfBuffer + deadzone)); 347b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org cntStateOld -= (mean_buffer_[0] <= (meanOfBuffer - deadzone)); 348b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org for (i = 1; i < mean_buffer_length_; i++) { 349b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org cntState = (mean_buffer_[i] >= (meanOfBuffer + deadzone)); 350b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org cntState -= (mean_buffer_[i] <= (meanOfBuffer - deadzone)); 351b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org if (cntStateOld == 0) { 352b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org cntStateOld = -cntState; 353470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com } 354b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org if (((cntState + cntStateOld) == 0) && (cntState != 0)) { 355b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org numZeros++; 356b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org cntStateOld = cntState; 357470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com } 358b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 359b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // END count zero crossings. 360b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 361b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org /* Frequency estimation according to: 362b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org * freqEst = numZeros * frame_rate / 2 / mean_buffer_length_; 363b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org * 364b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org * Resolution is set to Q4 365b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org */ 366b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org freqEst = ((numZeros * 90000) << 3); 367b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org freqEst /= 368b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org (timestamp_buffer_[0] - timestamp_buffer_[mean_buffer_length_ - 1]); 369b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org 370b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org /* Translate frequency estimate to regions close to 100 and 120 Hz */ 371b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org uint8_t freqState = 0; // Current translation state; 372b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // (0) Not in interval, 373b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // (1) Within valid interval, 374b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org // (2) Out of range 375b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org int32_t freqAlias = freqEst; 376b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org if (freqEst > kMinFrequencyToDetect) { 377b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org uint8_t aliasState = 1; 378a8565425bc6b1d3d7b5a0b0e0a49bc4e57c64f21mflodman while (freqState == 0) { 379b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org /* Increase frequency */ 380b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org freqAlias += (aliasState * frame_rate_); 381b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org freqAlias += ((freqEst << 1) * (1 - (aliasState << 1))); 382b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org /* Compute state */ 383b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org freqState = (abs(freqAlias - (100 << 4)) <= kFrequencyDeviation); 384b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org freqState += (abs(freqAlias - (120 << 4)) <= kFrequencyDeviation); 385b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org freqState += 2 * (freqAlias > ((120 << 4) + kFrequencyDeviation)); 386b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org /* Switch alias state */ 387b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org aliasState++; 388b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org aliasState &= 0x01; 389470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com } 390b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 391b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org /* Is frequency estimate within detection region? */ 392b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org if (freqState == 1) { 393b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org ret_val = 1; 394b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } else if (freqState == 0) { 395b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org ret_val = 2; 396b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } else { 397b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org ret_val = 0; 398b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org } 399b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org return ret_val; 400470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com} 401470e71d3649f6cac4688e83819640b012b5d38bbniklase@google.com 402b43d8078a152d91037b25175712a57a9a377220amikhal@webrtc.org} // namespace webrtc 403