From 3ca00b80b7446f2237ef530d135988db3ef7a2a0 Mon Sep 17 00:00:00 2001
From: Maxime Gervais <gervais.maxime@gmail.com>
Date: Mon, 3 Jun 2024 11:02:15 +0200
Subject: [PATCH] DeckLink Input: Add frame saturation statistics

Signed-off-by: Maxime Gervais <gervais.maxime@gmail.com>
---
 Project/GNU/CLI/Makefile.am           |   9 +-
 Project/MSVC2017/CLI/DVRescue.vcxproj |   2 +
 Project/MSVC2019/CLI/DVRescue.vcxproj |   2 +
 Project/MSVC2022/CLI/DVRescue.vcxproj |   2 +
 Source/Common/DecklinkWrapper.cpp     |  28 +++++
 Source/Common/Merge.h                 |   1 +
 Source/Common/Output_Xml.cpp          |  62 ++++++++++-
 Source/Common/Output_Xml.h            |   5 +
 Source/Common/ProcessFile.cpp         |   1 +
 Source/Common/ProcessFile.h           |  10 ++
 Source/Common/ProcessFileWrapper.cpp  |  20 +++-
 Source/Common/ProcessFileWrapper.h    |   5 +-
 Source/Common/SignalStats.cpp         | 151 ++++++++++++++++++++++++++
 Source/Common/SignalStats.h           |  37 +++++++
 Source/Common/SimulatorWrapper.cpp    |  24 ++++
 tools/dvrescue.xsd                    |  38 +++++++
 16 files changed, 387 insertions(+), 10 deletions(-)
 create mode 100644 Source/Common/SignalStats.cpp
 create mode 100644 Source/Common/SignalStats.h

diff --git a/Project/GNU/CLI/Makefile.am b/Project/GNU/CLI/Makefile.am
index 72cc22e7..dee9e8b0 100644
--- a/Project/GNU/CLI/Makefile.am
+++ b/Project/GNU/CLI/Makefile.am
@@ -53,14 +53,17 @@ if BUILD_AVFCTL
 endif
 
 if BUILD_DECKLINK
-    dvrescue_SOURCES += ../../../Source/Common/DecklinkWrapper.cpp
+    dvrescue_SOURCES += ../../../Source/Common/DecklinkWrapper.cpp \
+                        ../../../Source/Common/Output_Mkv.cpp \
+                        ../../../Source/Common/SignalStats.cpp
 endif
 
 if BUILD_DECKLINK
-    dvrescue_SOURCES += ../../../Source/Common/Output_Mkv.cpp
+
 else
 if BUILD_SIMULATOR
-    dvrescue_SOURCES += ../../../Source/Common/Output_Mkv.cpp
+    dvrescue_SOURCES += ../../../Source/Common/Output_Mkv.cpp \
+                        ../../../Source/Common/SignalStats.cpp
 endif
 endif
 
diff --git a/Project/MSVC2017/CLI/DVRescue.vcxproj b/Project/MSVC2017/CLI/DVRescue.vcxproj
index 7d5cb63b..aab516bc 100644
--- a/Project/MSVC2017/CLI/DVRescue.vcxproj
+++ b/Project/MSVC2017/CLI/DVRescue.vcxproj
@@ -161,6 +161,7 @@
     <ClCompile Include="..\..\..\Source\Common\Output_Webvtt.cpp" />
     <ClCompile Include="..\..\..\Source\Common\Output_Xml.cpp" />
     <ClCompile Include="..\..\..\Source\Common\Output_Mkv.cpp" />
+    <ClCompile Include="..\..\..\Source\Common\SignalStats.cpp" />
     <ClCompile Include="..\..\..\Source\Common\ProcessFile.cpp" />
     <ClCompile Include="..\..\..\Source\Common\ProcessFileWrapper.cpp" />
     <ClCompile Include="..\..\..\Source\Common\SimulatorWrapper.cpp" />
@@ -184,6 +185,7 @@
     <ClInclude Include="..\..\..\Source\Common\Output_Webvtt.h" />
     <ClInclude Include="..\..\..\Source\Common\Output_Xml.h" />
     <ClInclude Include="..\..\..\Source\Common\Output_Mkv.h" />
+    <ClInclude Include="..\..\..\Source\Common\SignalStats.h" />
     <ClInclude Include="..\..\..\Source\Common\ProcessFile.h" />
     <ClInclude Include="..\..\..\Source\Common\ProcessFileWrapper.h" />
     <ClInclude Include="..\..\..\Source\Common\SimulatorWrapper.h" />
diff --git a/Project/MSVC2019/CLI/DVRescue.vcxproj b/Project/MSVC2019/CLI/DVRescue.vcxproj
index 65106f7a..a5a792f4 100644
--- a/Project/MSVC2019/CLI/DVRescue.vcxproj
+++ b/Project/MSVC2019/CLI/DVRescue.vcxproj
@@ -162,6 +162,7 @@
     <ClCompile Include="..\..\..\Source\Common\Output_Webvtt.cpp" />
     <ClCompile Include="..\..\..\Source\Common\Output_Xml.cpp" />
     <ClCompile Include="..\..\..\Source\Common\Output_Mkv.cpp" />
+    <ClCompile Include="..\..\..\Source\Common\SignalStats.cpp" />
     <ClCompile Include="..\..\..\Source\Common\ProcessFile.cpp" />
     <ClCompile Include="..\..\..\Source\Common\ProcessFileWrapper.cpp" />
     <ClCompile Include="..\..\..\Source\Common\SimulatorWrapper.cpp" />
@@ -185,6 +186,7 @@
     <ClInclude Include="..\..\..\Source\Common\Output_Webvtt.h" />
     <ClInclude Include="..\..\..\Source\Common\Output_Xml.h" />
     <ClInclude Include="..\..\..\Source\Common\Output_Mkv.h" />
+    <ClInclude Include="..\..\..\Source\Common\SignalStats.h" />
     <ClInclude Include="..\..\..\Source\Common\ProcessFile.h" />
     <ClInclude Include="..\..\..\Source\Common\ProcessFileWrapper.h" />
     <ClInclude Include="..\..\..\Source\Common\SimulatorWrapper.h" />
diff --git a/Project/MSVC2022/CLI/DVRescue.vcxproj b/Project/MSVC2022/CLI/DVRescue.vcxproj
index beebcf95..94dbd8f6 100644
--- a/Project/MSVC2022/CLI/DVRescue.vcxproj
+++ b/Project/MSVC2022/CLI/DVRescue.vcxproj
@@ -162,6 +162,7 @@
     <ClCompile Include="..\..\..\Source\Common\Output_Webvtt.cpp" />
     <ClCompile Include="..\..\..\Source\Common\Output_Xml.cpp" />
     <ClCompile Include="..\..\..\Source\Common\Output_Mkv.cpp" />
+    <ClCompile Include="..\..\..\Source\Common\SignalStats.cpp" />
     <ClCompile Include="..\..\..\Source\Common\ProcessFile.cpp" />
     <ClCompile Include="..\..\..\Source\Common\ProcessFileWrapper.cpp" />
     <ClCompile Include="..\..\..\Source\Common\SimulatorWrapper.cpp" />
@@ -185,6 +186,7 @@
     <ClInclude Include="..\..\..\Source\Common\Output_Webvtt.h" />
     <ClInclude Include="..\..\..\Source\Common\Output_Xml.h" />
     <ClInclude Include="..\..\..\Source\Common\Output_Mkv.h" />
+    <ClInclude Include="..\..\..\Source\Common\SignalStats.h" />
     <ClInclude Include="..\..\..\Source\Common\ProcessFile.h" />
     <ClInclude Include="..\..\..\Source\Common\ProcessFileWrapper.h" />
     <ClInclude Include="..\..\..\Source\Common\SimulatorWrapper.h" />
diff --git a/Source/Common/DecklinkWrapper.cpp b/Source/Common/DecklinkWrapper.cpp
index dfd1a982..0981bcc3 100644
--- a/Source/Common/DecklinkWrapper.cpp
+++ b/Source/Common/DecklinkWrapper.cpp
@@ -43,6 +43,17 @@ static uint32_t decklink_audio_sources[Decklink_Audio_Source_Max] =
     bmdAudioConnectionMicrophone
 };
 
+//---------------------------------------------------------------------------
+static uint32_t decklink_pixel_formats[Decklink_Pixel_Format_Max] =
+{
+    bmdFormatUnspecified,
+    bmdFormat8BitYUV,
+    bmdFormat10BitYUV,
+    bmdFormat8BitARGB,
+    bmdFormat8BitBGRA,
+    bmdFormat10BitRGB
+};
+
 //---------------------------------------------------------------------------
 static uint32_t decklink_timecode_formats[Decklink_Timecode_Format_Max] =
 {
@@ -129,9 +140,26 @@ HRESULT DecklinkWrapper::CaptureDelegate::VideoInputFrameArrived(IDeckLinkVideoI
             }
         }
 
+        uint8_t PixelFormat = (uint8_t)Decklink_Pixel_Format_Unspecified;
+        switch (VideoFrame->GetPixelFormat())
+        {
+            case bmdFormat8BitYUV:
+                PixelFormat = (uint8_t)Decklink_Pixel_Format_8BitYUV; break;
+            case bmdFormat10BitYUV:
+                PixelFormat = (uint8_t)Decklink_Pixel_Format_10BitYUV; break;
+            case bmdFormat8BitARGB:
+                PixelFormat = (uint8_t)Decklink_Pixel_Format_8BitARGB; break;
+            case bmdFormat8BitBGRA:
+                PixelFormat = (uint8_t)Decklink_Pixel_Format_8BitBGRA; break;
+            case bmdFormat10BitRGB:
+                PixelFormat = (uint8_t)Decklink_Pixel_Format_10BitRGB; break;
+            default:;
+        }
+
         decklink_frame Buffer = {
             .Width = (uint32_t)VideoFrame->GetWidth(),
             .Height = (uint32_t)VideoFrame->GetHeight(),
+            .Pixel_Format = PixelFormat,
             .Video_Buffer = (uint8_t*)VideoBuffer,
             .Video_Buffer_Size = VideoBufferSize,
             .Audio_Buffer = (uint8_t*)AudioBuffer,
diff --git a/Source/Common/Merge.h b/Source/Common/Merge.h
index 0255a742..39012cd2 100644
--- a/Source/Common/Merge.h
+++ b/Source/Common/Merge.h
@@ -28,6 +28,7 @@ extern const char* Control_Port;
 extern uint8_t DeckLinkVideoMode;
 extern uint8_t DeckLinkVideoSource;
 extern uint8_t DeckLinkAudioSource;
+extern uint8_t DeckLinkPixelFormat;
 extern uint8_t DeckLinkTimecodeFormat;
 #endif
 extern size_t Merge_Rewind_Count;
diff --git a/Source/Common/Output_Xml.cpp b/Source/Common/Output_Xml.cpp
index b8ce44cc..aeef607e 100644
--- a/Source/Common/Output_Xml.cpp
+++ b/Source/Common/Output_Xml.cpp
@@ -136,7 +136,7 @@ static void Aud_Element(string& Text, size_t o, size_t n, vector<uint16_t> Audio
 
 #if defined(ENABLE_DECKLINK) || defined(ENABLE_SIMULATOR)
 //---------------------------------------------------------------------------
-static string decklink_videosource_to_string(uint8_t value)
+string decklink_videosource_to_string(uint8_t value)
 {
     string ToReturn;
     switch (value)
@@ -153,7 +153,7 @@ static string decklink_videosource_to_string(uint8_t value)
 }
 
 //---------------------------------------------------------------------------
-static string decklink_audiosource_to_string(uint8_t value)
+string decklink_audiosource_to_string(uint8_t value)
 {
     string ToReturn;
     switch (value)
@@ -170,7 +170,23 @@ static string decklink_audiosource_to_string(uint8_t value)
 }
 
 //---------------------------------------------------------------------------
-static string decklink_timecodeformat_to_string(uint8_t value)
+string decklink_pixelformat_to_string(uint8_t value)
+{
+    string ToReturn;
+    switch (value)
+    {
+    case (uint8_t)Decklink_Pixel_Format_8BitYUV: ToReturn = "uyvy"; break;
+    case (uint8_t)Decklink_Pixel_Format_10BitYUV: ToReturn = "v210"; break;
+    case (uint8_t)Decklink_Pixel_Format_8BitARGB: ToReturn = "argb"; break;
+    case (uint8_t)Decklink_Pixel_Format_8BitBGRA: ToReturn = "bgra"; break;
+    case (uint8_t)Decklink_Pixel_Format_10BitRGB: ToReturn = "r210"; break;
+    default:;
+    }
+    return ToReturn;
+}
+
+//---------------------------------------------------------------------------
+string decklink_timecodeformat_to_string(uint8_t value)
 {
     string ToReturn;
     switch (value)
@@ -359,6 +375,14 @@ return_value Output_Xml(ostream& Out, std::vector<file*>& PerFile, bitset<Option
                 Text += to_string(File->Wrapper->FramesInfo.video_height);
                 Text += '\"';
             }
+
+            if (File->Wrapper->FramesInfo.pixel_format)
+            {
+                Text += " pixel_format=\"";
+                Text += decklink_pixelformat_to_string(File->Wrapper->FramesInfo.pixel_format);
+                Text += '\"';
+            }
+
             if (File->Wrapper->FramesInfo.video_rate_num)
             {
                 Text += " video_rate=\"";
@@ -418,7 +442,37 @@ return_value Output_Xml(ostream& Out, std::vector<file*>& PerFile, bitset<Option
                         Text += " tc_nc=\"1\"";
                 }
 
-                Text += "/>\n";
+                if (File->Wrapper->FramesInfo.frames[Pos].st.HasValue())
+                {
+                    Text += ">\n";
+                    Text += "\t\t\t\t<signalstats";
+                    if (!isnan(File->Wrapper->FramesInfo.frames[Pos].st.SatAvg))
+                    {
+                        stringstream ss;
+                        ss << File->Wrapper->FramesInfo.frames[Pos].st.SatAvg;
+                        Text += " satavg=\"";
+                        Text += ss.str();
+                        Text += '\"';
+                    }
+
+                    if (File->Wrapper->FramesInfo.frames[Pos].st.SatHi != (uint16_t)-1)
+                    {
+                        Text += " sathi=\"";
+                        Text += to_string(File->Wrapper->FramesInfo.frames[Pos].st.SatHi);
+                        Text += '\"';
+                    }
+
+                    if (File->Wrapper->FramesInfo.frames[Pos].st.SatMax != (uint16_t)-1)
+                    {
+                        Text += " satmax=\"";
+                        Text += to_string(File->Wrapper->FramesInfo.frames[Pos].st.SatMax);
+                        Text += '\"';
+                    }
+                    Text += "/>\n";
+                    Text += "\t\t\t</frame>\n";
+                }
+                else
+                    Text += "/>\n";
             }
 
             // Media footer
diff --git a/Source/Common/Output_Xml.h b/Source/Common/Output_Xml.h
index ddebed85..775c8cb3 100644
--- a/Source/Common/Output_Xml.h
+++ b/Source/Common/Output_Xml.h
@@ -13,3 +13,8 @@ class file;
 //---------------------------------------------------------------------------
 
 return_value Output_Xml(ostream& Out, vector<file*>& PerFile, bitset<Option_Max> Options, ostream* Err = nullptr);
+
+//---------------------------------------------------------------------------
+#if defined(ENABLE_DECKLINK) || defined(ENABLE_SIMULATOR)
+string decklink_pixelformat_to_string(uint8_t value);
+#endif
diff --git a/Source/Common/ProcessFile.cpp b/Source/Common/ProcessFile.cpp
index 46c0b017..3a0241a6 100644
--- a/Source/Common/ProcessFile.cpp
+++ b/Source/Common/ProcessFile.cpp
@@ -43,6 +43,7 @@ bool DeckLinkNativeControl = false;
 uint8_t DeckLinkVideoMode = (uint8_t)Decklink_Video_Mode_NTSC;
 uint8_t DeckLinkVideoSource = (uint8_t)Decklink_Video_Source_SDI;
 uint8_t DeckLinkAudioSource = (uint8_t)Decklink_Audio_Source_Embedded;
+uint8_t DeckLinkPixelFormat = (uint8_t)Decklink_Pixel_Format_Unspecified;
 uint8_t DeckLinkTimecodeFormat = (uint8_t)Decklink_Timecode_Format_VITC;
 #endif
 bool InControl = false;
diff --git a/Source/Common/ProcessFile.h b/Source/Common/ProcessFile.h
index 67a6e6af..6fbe7e37 100644
--- a/Source/Common/ProcessFile.h
+++ b/Source/Common/ProcessFile.h
@@ -82,6 +82,16 @@ enum decklink_audio_source {
     Decklink_Audio_Source_Max
 };
 
+enum decklink_pixel_format {
+    Decklink_Pixel_Format_Unspecified,
+    Decklink_Pixel_Format_8BitYUV,
+    Decklink_Pixel_Format_10BitYUV,
+    Decklink_Pixel_Format_8BitARGB,
+    Decklink_Pixel_Format_8BitBGRA,
+    Decklink_Pixel_Format_10BitRGB,
+    Decklink_Pixel_Format_Max
+};
+
 enum decklink_timecode_format {
     Decklink_Timecode_Format_RP188_VITC,
     Decklink_Timecode_Format_RP188_VITC2,
diff --git a/Source/Common/ProcessFileWrapper.cpp b/Source/Common/ProcessFileWrapper.cpp
index dceb220a..bb8eb4a6 100644
--- a/Source/Common/ProcessFileWrapper.cpp
+++ b/Source/Common/ProcessFileWrapper.cpp
@@ -8,10 +8,13 @@
 #include "Common/ProcessFile.h"
 
 #if defined(ENABLE_DECKLINK) || defined(ENABLE_SIMULATOR)
+#include "Common/SignalStats.h"
 #include "Common/Output_Mkv.h"
+#include "Common/Output_Xml.h"
 #include "Common/Merge.h"
 
 #include <iostream>
+#include <cmath>
 #endif
 
 using namespace std;
@@ -78,11 +81,15 @@ void FileWrapper::Parse_Buffer(const uint8_t *Buffer, size_t Buffer_Size)
     if (IsMatroska)
     {
         decklink_frame* Frame=(decklink_frame*)Buffer;
+        if (!FramesInfo.pixel_format && Frame->Pixel_Format)
+            FramesInfo.pixel_format = Frame->Pixel_Format;
+
         for (matroska_output& Output : Outputs)
         {
             Output.Writer->write_frame((const char*)Frame->Video_Buffer, Frame->Video_Buffer_Size,
                                        (const char*)Frame->Audio_Buffer, Frame->Audio_Buffer_Size, Frame->TC);
         }
+        SignalStats::Stats ST = SignalStats::ComputeStats(Frame);
 
         bool TimeCode_Repeat = false;
         bool TimeCode_NonConsecutive = false;
@@ -108,6 +115,7 @@ void FileWrapper::Parse_Buffer(const uint8_t *Buffer, size_t Buffer_Size)
         double FrameRate = (double)FramesInfo.video_rate_num / FramesInfo.video_rate_den;
         double ElapsedTime = (double)FramesInfo.frames.size() / FrameRate;
         FramesInfo.frames.push_back(decklink_framesinfo::frame {
+            ST,
             Frame->TC,
             TimeCode_Repeat,
             (uint8_t)(TimeCode_NonConsecutive ? (TimeCode_NonConsecutive_IsLess ? 2 : 1) : 0),
@@ -119,12 +127,16 @@ void FileWrapper::Parse_Buffer(const uint8_t *Buffer, size_t Buffer_Size)
         {
             if (!FrameCount)
             {
-                cout << "FramePos,abst,abst_r,abst_nc,tc,tc_r,tc_nc,rdt,rdt_r,rdt_nc,rec_start,rec_end,Used,Status,Comments,BlockErrors,BlockErrors_Even,IssueFixed"
+                cout << "FramePos,abst,abst_r,abst_nc,tc,tc_r,tc_nc,rdt,rdt_r,rdt_nc,rec_start,rec_end,pix_fmt,satavg,sathi,satmax,Used,Status,Comments,BlockErrors,BlockErrors_Even,IssueFixed"
                      << (Verbosity > 5 ? ",SourceSpeed,FrameSpeed,InputPos,OutputPos" : "")
                      << (ShowFrames_Intermediate ? ",RewindStatus" : "")
                      << endl;
             }
 
+            stringstream ss;
+            if (!isnan(ST.SatAvg))
+                ss << ST.SatAvg;
+
             cout << FrameCount++ // framePos
                  << "," // abst
                     "," // abst_r
@@ -138,7 +150,11 @@ void FileWrapper::Parse_Buffer(const uint8_t *Buffer, size_t Buffer_Size)
                     "," // rec_start
                     "," // rec_end
                     "," // Used
-                    "," // Status
+                 << "," << decklink_pixelformat_to_string(Frame->Pixel_Format) // pix_fmt
+                 << "," << ss.str() // satavg
+                 << "," << (ST.SatHi != (uint16_t)-1 ? to_string(ST.SatHi) : "") // sathi
+                 << "," << (ST.SatMax != (uint16_t)-1 ? to_string(ST.SatMax) : "")  // satmax
+                 << "," // Status
                     "," // Comments
                     "," // BlockErrors
                     "," // BlocksErrors_Even
diff --git a/Source/Common/ProcessFileWrapper.h b/Source/Common/ProcessFileWrapper.h
index cd5b5ccb..bcbce919 100644
--- a/Source/Common/ProcessFileWrapper.h
+++ b/Source/Common/ProcessFileWrapper.h
@@ -13,6 +13,7 @@
 #include <vector>
 
 #include "ThirdParty/TimeCode/TimeCode.h"
+#include "Common/SignalStats.h"
 
 class file;
 #if defined(ENABLE_DECKLINK) || defined(ENABLE_SIMULATOR)
@@ -37,6 +38,7 @@ struct decklink_frame
 {
     uint32_t                    Width = 0;
     uint32_t                    Height = 0;
+    uint8_t                     Pixel_Format = 0;
     uint8_t*                    Video_Buffer = nullptr;
     size_t                      Video_Buffer_Size = 0;
     uint8_t*                    Audio_Buffer = nullptr;
@@ -46,6 +48,7 @@ struct decklink_frame
 
 struct decklink_framesinfo {
     struct frame {
+        SignalStats::Stats st;
         TimeCode tc;
         bool tc_r;
         uint8_t tc_nc;
@@ -57,9 +60,9 @@ struct decklink_framesinfo {
     uint32_t video_height = 0;
     uint32_t video_rate_num = 0;
     uint32_t video_rate_den = 0;
+    uint8_t pixel_format = 0;
     uint8_t audio_channels = 0;
     uint32_t audio_rate = 0;
-
     std::vector<frame> frames;
 };
 
diff --git a/Source/Common/SignalStats.cpp b/Source/Common/SignalStats.cpp
new file mode 100644
index 00000000..a838dae0
--- /dev/null
+++ b/Source/Common/SignalStats.cpp
@@ -0,0 +1,151 @@
+/*  Copyright (c) MIPoPS. All Rights Reserved.
+ *
+ *  Use of this source code is governed by a BSD-3-Clause license that can
+ *  be found in the LICENSE.txt file in the root of the source tree.
+ */
+
+//---------------------------------------------------------------------------
+
+//---------------------------------------------------------------------------
+#include "Common/SignalStats.h"
+#include "Common/ProcessFileWrapper.h"
+
+#include <ZenLib/Utils.h>
+
+#include <cstring>
+#include <cmath>
+#include <iostream>
+
+
+using namespace std;
+//---------------------------------------------------------------------------
+
+//***************************************************************************
+// Helpers
+//***************************************************************************
+
+//---------------------------------------------------------------------------
+#define BITDEPTH_8_DEPTH   8
+#define BITDEPTH_8_OFFSET  (1 << (BITDEPTH_8_DEPTH - 1))
+#define BITDEPTH_10_DEPTH  10
+#define BITDEPTH_10_OFFSET (1 << (BITDEPTH_10_DEPTH - 1))
+
+//---------------------------------------------------------------------------
+void ComputeSatUYVY422(const decklink_frame* Frame, size_t* Sats)
+{
+    size_t LineSize = Frame->Width * 16 / 8;
+    if (Frame->Video_Buffer_Size < Frame->Height * LineSize)
+        return;
+
+    const uint8_t* Buffer = (const uint8_t*)Frame->Video_Buffer;
+    for (size_t Height = 0; Height < Frame->Height; Height++)
+    {
+        for (size_t Width = 0; Width < Frame->Width; Width += 2)
+        {
+            size_t Buffer_Offset = (Height * LineSize) + (Width * 16 / 8);
+            int8_t U, V;
+
+            U = Buffer[Buffer_Offset + 2] - BITDEPTH_8_OFFSET;
+            V = Buffer[Buffer_Offset]     - BITDEPTH_8_OFFSET;
+            Sats[(uint8_t)sqrt(U * U + V * V)]++;
+        }
+    }
+}
+
+//---------------------------------------------------------------------------
+void ComputeSatV210(const decklink_frame* Frame, size_t* Sats)
+{
+    //TODO: Indianness
+    size_t LineSize = (((Frame->Width + 47) / 48) * 128) / sizeof(uint32_t);
+    if (Frame->Video_Buffer_Size < Frame->Height * (LineSize / 6 * 4) * 4)
+        return;
+
+    const uint32_t* Buffer = (const uint32_t*)Frame->Video_Buffer;
+    for (size_t Height = 0; Height < Frame->Height; Height++)
+    {
+        for (size_t Width = 0; Width < Frame->Width; Width += 6)
+        {
+            size_t Buffer_Offset = Height * LineSize + (Width / 6 * 4);
+            int32_t U, V;
+
+            U = (((Buffer[Buffer_Offset]          ) & 0x3FF) - BITDEPTH_10_OFFSET);
+            V = (((Buffer[Buffer_Offset]     >> 20) & 0x3FF) - BITDEPTH_10_OFFSET);
+            Sats[(uint16_t)sqrt(U * U + V * V)]++;
+
+            U = (((Buffer[Buffer_Offset + 1] >> 10) & 0x3FF) - BITDEPTH_10_OFFSET);
+            V = (((Buffer[Buffer_Offset + 2]      ) & 0x3FF) - BITDEPTH_10_OFFSET);
+            Sats[(uint16_t)sqrt(U * U + V * V)]++;
+
+            U = (((Buffer[Buffer_Offset + 2] >> 20) & 0x3FF) - BITDEPTH_10_OFFSET);
+            V = (((Buffer[Buffer_Offset + 3] >> 10) & 0x3FF) - BITDEPTH_10_OFFSET);
+            Sats[(uint16_t)sqrt(U * U + V * V)]++;
+        }
+    }
+}
+
+//***************************************************************************
+// ComputeStats
+//***************************************************************************
+
+//---------------------------------------------------------------------------
+SignalStats::Stats SignalStats::ComputeStats(const decklink_frame* Frame)
+{
+    Stats Result;
+
+    if (!Frame)
+        return Stats();
+
+    size_t Sats_Count = 0, Sats_Size = 0, *Sats = nullptr;
+    switch (Frame->Pixel_Format)
+    {
+    case 1://UYVY
+    {
+        Sats_Count = (Frame->Width / 2) * Frame->Height;
+        Sats_Size = 1 << BITDEPTH_8_DEPTH;
+        Sats = new size_t[Sats_Size];
+        memset(Sats, 0, Sats_Size * sizeof(size_t));
+
+        ComputeSatUYVY422(Frame, Sats);
+    }
+    break;
+    case 2: //v210
+    {
+        Sats_Count = (Frame->Width / 2) * Frame->Height;
+        Sats_Size = 1 << BITDEPTH_10_DEPTH;
+        Sats = new size_t[Sats_Size];
+        memset (Sats, 0, Sats_Size * sizeof(size_t));
+
+        ComputeSatV210(Frame, Sats);
+    }
+    break;
+    default:;
+    }
+
+    if (Sats)
+    {
+        int32_t Hi = -1;
+        int32_t Max = -1;
+        int32_t Count = 0;
+        int64_t Total = 0;
+        int32_t Sats_Count_90Percentile = lrint(Sats_Count * 90 / 100.0f);
+        for (uint16_t Sat = 0; Sat < Sats_Size; Sat++)
+        {
+            if (Sats[Sat] > 0)
+            {
+                Count += Sats[Sat];
+                Total += Sat * Sats[Sat];
+                if (Hi < 0 && Count >= Sats_Count_90Percentile)
+                    Hi = Sat;
+
+                Max = Sat;
+            }
+        }
+        Result.SatAvg = 1.0f * Total / Sats_Count;
+        Result.SatHi = Hi;
+        Result.SatMax = (uint16_t)Max;
+
+        delete[] Sats;
+    }
+
+    return Result;
+}
diff --git a/Source/Common/SignalStats.h b/Source/Common/SignalStats.h
new file mode 100644
index 00000000..fd34c7eb
--- /dev/null
+++ b/Source/Common/SignalStats.h
@@ -0,0 +1,37 @@
+/*  Copyright (c) MIPoPS. All Rights Reserved.
+ *
+ *  Use of this source code is governed by a BSD-3-Clause license that can
+ *  be found in the LICENSE.txt file in the root of the source tree.
+ */
+
+//---------------------------------------------------------------------------
+#pragma once
+
+#include <cstdint>
+#include <cmath>
+//***************************************************************************
+// Class SignalStats
+//***************************************************************************
+
+struct decklink_frame;
+class SignalStats {
+public:
+    struct Stats {
+        double SatAvg = -NAN;
+        uint16_t SatHi = -1;
+        uint16_t SatMax = -1;
+
+        bool HasValue() const
+        {
+            return !std::isnan(SatAvg) ||
+                   SatHi != (uint16_t)-1 ||
+                   SatMax != (uint16_t)-1;
+        };
+    };
+
+    static Stats ComputeStats(const decklink_frame* Frame);
+
+private:
+    int Width;
+    int Height;
+};
\ No newline at end of file
diff --git a/Source/Common/SimulatorWrapper.cpp b/Source/Common/SimulatorWrapper.cpp
index e4bc695c..76ef9f8e 100644
--- a/Source/Common/SimulatorWrapper.cpp
+++ b/Source/Common/SimulatorWrapper.cpp
@@ -5,6 +5,7 @@
  */
 
 #include "Common/SimulatorWrapper.h"
+#include "Common/ProcessFile.h"
 #include "ZenLib/File.h"
 #include "ZenLib/ZtringListList.h"
 #include "TimeCode.h"
@@ -26,6 +27,13 @@ using namespace ZenLib;
     #define SIM_PREFIX __T("/tmp/dvrescue_simulator.")
 #endif
 
+
+#define V_UNCOMPRESSED_ID "V_UNCOMPRESSED"
+#define V210_ID           "V_MS/VFW/FOURCC"
+
+#define V_UNCOMPRESSED_UYVY_CC 0x55595659
+#define V210_CC                0x76323130
+
 //---------------------------------------------------------------------------
 struct ctl
 {
@@ -532,6 +540,7 @@ SimulatorWrapper::SimulatorWrapper(std::size_t DeviceIndex)
         else
             P->Buffer_Size = 120000;
         P->Buffer = new int8u[P->Buffer_Size];
+        memset(P->Buffer, 0, P->Buffer_Size);
         if (P->IsMatroska)
         {
             P->F.back()->Read(P->Buffer, P->Buffer_Size);
@@ -753,6 +762,21 @@ bool SimulatorWrapper::WaitForSessionEnd(uint64_t Timeout)
                 P->Mutex.unlock();
                 break;
             }
+            if (!P->Tracks.empty())
+            {
+                if (P->Tracks.back().CodecID == V_UNCOMPRESSED_ID)
+                {
+                    switch (P->Tracks.back().FourCC)
+                    {
+                    case V_UNCOMPRESSED_UYVY_CC:
+                        P->Decklink_Sim.Pixel_Format = Decklink_Pixel_Format_8BitYUV; break;
+                    default:;
+                    }
+                }
+                else if (P->Tracks.back().CodecID == V210_ID && P->Tracks.back().FourCC == V210_CC)
+                    P->Decklink_Sim.Pixel_Format = Decklink_Pixel_Format_10BitYUV;
+            }
+
             P->Wrapper->Parse_Buffer((uint8_t*)&P->Decklink_Sim, sizeof(P->Decklink_Sim));
         }
         else
diff --git a/tools/dvrescue.xsd b/tools/dvrescue.xsd
index 69eefbe6..59004827 100644
--- a/tools/dvrescue.xsd
+++ b/tools/dvrescue.xsd
@@ -180,6 +180,7 @@
       </xsd:annotation>
     </xsd:attribute>
     <xsd:attribute name="size" type="sizeType"/>
+    <xsd:attribute name="pixel_format" type="pixel_formatType"/>
     <xsd:attribute name="video_rate" type="video_rateType"/>
     <xsd:attribute name="chroma_subsampling" type="chroma_subsamplingType"/>
     <xsd:attribute name="aspect_ratio" type="aspect_ratioType"/>
@@ -218,6 +219,18 @@
     </xsd:restriction>
   </xsd:simpleType>
 
+  <xsd:simpleType name="pixel_formatType">
+    <xsd:annotation>
+      <xsd:documentation xml:lang="en">
+        Pixel format.
+      </xsd:documentation>
+    </xsd:annotation>
+    <xsd:restriction base="xsd:string">
+      <xsd:enumeration value="uyvy"/>
+      <xsd:enumeration value="v210"/>
+    </xsd:restriction>
+  </xsd:simpleType>
+
   <xsd:simpleType name="video_rateType">
     <xsd:annotation>
       <xsd:documentation xml:lang="en">
@@ -309,6 +322,7 @@
       <xsd:element name="dseq" type="dseqType" minOccurs="0" maxOccurs="unbounded"/>
       <xsd:element name="sta" type="staType" minOccurs="0" maxOccurs="unbounded"/>
       <xsd:element name="aud" type="audType" minOccurs="0" maxOccurs="unbounded"/>
+      <xsd:element name="signalstats" type="signalstatsType" minOccurs="0" maxOccurs="unbounded"/>
     </xsd:sequence>
     <xsd:attribute name="n" type="xsd:integer">
       <xsd:annotation>
@@ -635,4 +649,28 @@
     </xsd:attribute>
   </xsd:complexType>
 
+  <xsd:complexType name="signalstatsType">
+    <xsd:attribute name="satavg" type="xsd:float">
+      <xsd:annotation>
+        <xsd:documentation xml:lang="en">
+          Average saturation of the frame.
+        </xsd:documentation>
+      </xsd:annotation>
+    </xsd:attribute>
+    <xsd:attribute name="sathi" type="xsd:integer">
+      <xsd:annotation>
+        <xsd:documentation xml:lang="en">
+          Highest saturation of the frame excluding the to 10% values.
+        </xsd:documentation>
+      </xsd:annotation>
+    </xsd:attribute>
+    <xsd:attribute name="satmax" type="xsd:integer">
+      <xsd:annotation>
+        <xsd:documentation xml:lang="en">
+          Highest saturation of the frame.
+        </xsd:documentation>
+      </xsd:annotation>
+    </xsd:attribute>
+  </xsd:complexType>
+
 </xsd:schema>