# Patent

## Title
Video quality objective assessment method based on spatiotemporal domain structure

## Abstract

A video quality objective assessment method based on a spatiotemporal domain structure firstly combines a spatiotemporal domain gradient magnitude and color information for calculating a spatiotemporal domain local similarity, and then uses variance fusion for spatial domain fusion. The spatiotemporal domain local similarity is fused into frame-level objective quality value, and then a temporal domain fusion model is established by simulating three important global temporal effects, which are a smoothing effect, an asymmetric track effects and a recency effect, of a human visual system. Finally, the objective quality values of the distorted video sequence are obtained. By modeling the human visual temporal domain effect, the temporal domain weighting method of the present invention is able to accurately and efficiently evaluate the objective quality of the distorted video.

## Background
CROSS REFERENCE OF RELATED APPLICATION
The present invention claims priority under 35 U.S.C. 119(a-d) to CN 201610367702.2, filed May 27, 2016.
BACKGROUND OF THE PRESENT INVENTION
Field of Invention
The present invention relates to a video quality assessment method, and more particularly to a video quality objective assessment method based on a spatiotemporal domain structure.
Description of Related Arts
Video quality assessment (VQA for short) plays a key role in the rapid development of video-related applications, which is commonly used in quality control of a wide variety of video services and performance comparison of various video processing algorithms. Video quality objective assessment focuses on estimating the video quality by algorithms, so as to approximate the video quality perceived by human visual system. According to availability of original video without distortion, video quality objective assessment method is generally divided into three categories: full-reference (FR), reduced-reference and no-reference video quality assessment methods. Due to a higher accuracy, the full reference video quality assessment method (FR-VQA for short) is widely used in perceptual video processing such as encoding module selection, parameter quantization, rate control, and error concealment in video compression.
Full-reference video quality assessment method requires the ability to quickly and accurately evaluate the video objective quality, so real-time processing capability is very important, such as online source video quality monitoring and distortion metrics of rate-distortion optimized video encoder. Besides, in other applications, a low complexity is also a very important performance requirement for the full reference video quality assessment method. Pixel-level peak signal-to-noise ratio (PSNR for short) in the full-reference assessment method is conventionally the most widely used performance quantized indicator in video processing. Due to advantages such as convenient implement, fast evaluation, stable performance, and clear physical meaning, the peak signal-to-noise ratio is still the most widely used objective assessment method for most widely used video coding standard H.264/AVC and the latest H.265/HEVC. However, distortion of each pixel is treated equally without considering perceptual characteristics of the human visual system, which causes a low consistency between the peak signal-to-noise ratio and the subjective quality perception, thus impeding the progress of video processing technology, especially the progress of the video compression efficiency.
Conventional full-reference video quality assessment methods are divided into two categories. The first one, which is the most direct video quality assessment method, is to use effective image quality assessment method on independent frame, and then use average or weighted average for obtaining the video quality. However, such video quality assessment method lacks temporal information and assessment effect is poor. The second one respectively investigates spatial and temporal domains distortion for obtaining the video quality, or directly evaluates according to the spatial domain information. Although assessment effect of such method is better than that of the peak signal-to-noise ratio, complexity thereof is high and some need very time-consuming motion estimation. Therefore, disadvantages, such as difficult implement, impossible real-time processing and difficult integration, hinder the wide application of such video quality assessment method.
SUMMARY OF THE PRESENT INVENTION
An object of the present invention is to provide a video quality objective assessment method based on a spatiotemporal domain structure, which effectively improves relativity between objective assessment results and subjective perception, and lowers calculation complexity.
Accordingly, in order to accomplish the above object, the present invention provides a video quality objective assessment method based on a spatiotemporal domain structure, comprising steps of:
(1) marking a reference video sequence without distortion as Sr, and marking a distorted video sequence, which is obtained after the Sr is distorted, as Sd; wherein a total Sr frame quantity is F, a total Sd frame quantity is also F, and F>1; widths of images in both the Sr and the Sd are W, and heights of the images in both the Sr and the Sd are H; defining an image luminance component sequence of the images in the Sr as a luminance component sequence of the Sr and marking as Yr; defining a first image chrominance component sequence of the images in the Sr as a first chrominance component sequence of the Sr and marking as Ur; defining a second image chrominance component sequence of the images in the Sr as a second chrominance component sequence of the Sr and marking as Vr; defining an image luminance component sequence of the images in the Sd as a luminance component sequence of the Sd and marking as Yd; defining a first image chrominance component sequence of the images in the Sd as a first chrominance component sequence of the Sd and marking as Ud; defining a second image chrominance component sequence of the images in the Sd as a second chrominance component sequence of the Sd and marking as Vd; wherein widths of images in the Yr, the Ur, the Vr, the Yd, the Ud and the Vd are W, and heights of the images in the Yr, the Ur, the Vr, the Yd, the Ud and the Vd are H;
(2) calculating a spatiotemporal domain gradient magnitude sequence of the Yr with a three-dimensional Prewitt operator and marking as Gr, and marking a pixel value of a pixel at a position of (x,y) in a number t frame in the Gr as Gr(x,y,t), wherein




  
    
      
        G
        r
      
      ⁡
      
        (
        
          x
          ,
          y
          ,
          t
        
        )
      
    
    =
    
      
        
          
            (
            
              
                Y
                rx
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          
            (
            
              
                Y
                ry
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          
            (
            
              
                Y
                rt
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
      
    
  
  ,




Yrx=Yr{circle around (×)}Fx, Yry=Yr{circle around (×)}Fy, Yrt=Yr{circle around (×)}Ft;

similarly, calculating a spatiotemporal domain gradient magnitude sequence of the Yd with the three-dimensional Prewitt operator and marking as Gd, and marking a pixel value of a pixel at a position of (x,y) in a number t frame in the Gd as Gd(x,y,t), wherein




  
    
      
        G
        d
      
      ⁡
      
        (
        
          x
          ,
          y
          ,
          t
        
        )
      
    
    =
    
      
        
          
            (
            
              
                Y
                dx
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          
            (
            
              
                Y
                dy
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          
            (
            
              
                Y
                dt
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
      
    
  
  ,




Ydr=Yd{circle around (×)}Fx, Ydy=Yd{circle around (×)}Fy, Ydt=Yd{circle around (×)}Ft;

wherein an initial value of the t is 1, 1≦t≦F, 1≦x≦W, 1≦y≦H; Yrx(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a horizontal gradient magnitude sequence Yrx of the Yr, Yry(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a vertical gradient magnitude sequence Yry of the Yr, and Yrt(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a temporal gradient magnitude sequence Yrt of the Yr; Ydx(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a horizontal gradient magnitude sequence Ydx of the Yd, Ydy(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a vertical gradient magnitude sequence Ydy of the Yd, and Ydt(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a temporal gradient magnitude sequence Ydt of the Yd; a symbol {circle around (×)} is a zeros truncated convolution symbol; after convolution, dimensions of the Yrx, the Yry and the Yrt are same as a dimension of the Yr, and dimensions of the Ydx, Ydy and Ydt are same as a dimension of the Yd; Fx, Fy and Ft correspond to a horizontal mask, a vertical mask and a temporal mask of the three-dimensional Prewitt operator;
(3) calculating a spatiotemporal domain local gradient similarity between each pixel point in each frame in the Sr and a corresponding pixel point in a corresponding frame in the Sd; marking the spatiotemporal domain local gradient similarity between a pixel point at a position of (x,y) in a number t frame in the Sr and a pixel point at a position of (x,y) in a number t frame in the Sd as Gsim(x,y,t); wherein




  
    
      
        G
        sim
      
      ⁡
      
        (
        
          x
          ,
          y
          ,
          t
        
        )
      
    
    =
    
      
        
          2
          ⁢
          
            
          
          ⁢
          
            
              G
              r
            
            ⁡
            
              (
              
                x
                ,
                y
                ,
                t
              
              )
            
          
          ×
          
            
              G
              d
            
            ⁡
            
              (
              
                x
                ,
                y
                ,
                t
              
              )
            
          
        
        +
        
          c
          1
        
      
      
        
          
            (
            
              
                G
                r
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          
            (
            
              
                G
                d
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          c
          1
        
      
    
  
  ,




c1 is a positive constant preventing the fractional from being meaningless;

(4) calculating a spatiotemporal domain local color similarity between each pixel point in each frame in the Sr and the corresponding pixel point in the corresponding frame in the Sd; marking the spatiotemporal domain local color similarity between the pixel point at the position of (x,y) in the number t frame in the Sr and the pixel point at the position of (x,y) in the number t frame in the Sd as Csim(x,y,t); wherein




  
    
      
        C
        sim
      
      ⁡
      
        (
        
          x
          ,
          y
          ,
          t
        
        )
      
    
    =
    
      
        
          
            2
            ⁢
            
              
            
            ⁢
            
              
                U
                r
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            ×
            
              
                U
                d
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
          
          +
          
            c
            2
          
        
        
          
            
              (
              
                
                  U
                  r
                
                ⁡
                
                  (
                  
                    x
                    ,
                    y
                    ,
                    t
                  
                  )
                
              
              )
            
            2
          
          +
          
            
              (
              
                
                  U
                  d
                
                ⁡
                
                  (
                  
                    x
                    ,
                    y
                    ,
                    t
                  
                  )
                
              
              )
            
            2
          
          +
          
            c
            2
          
        
      
      ×
      
        
          
            2
            ⁢
            
              
            
            ⁢
            
              
                V
                r
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            ×
            
              
                V
                d
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
          
          +
          
            c
            3
          
        
        
          
            
              (
              
                
                  V
                  r
                
                ⁡
                
                  (
                  
                    x
                    ,
                    y
                    ,
                    t
                  
                  )
                
              
              )
            
            2
          
          +
          
            
              (
              
                
                  V
                  d
                
                ⁡
                
                  (
                  
                    x
                    ,
                    y
                    ,
                    t
                  
                  )
                
              
              )
            
            2
          
          +
          
            c
            3
          
        
      
    
  
  ,




Ur(x,y,t) refers to a pixel value of a pixel point at a position of (x,y) in a number t frame in the Ur, which is also a pixel value of a pixel point at a position of (x,y) in a first chrominance component in a number t frame in the Sr; Vr(x,y,t) refers to a pixel value of a pixel point at a position of (x,y) in a number t frame in the Vr, which is also a pixel value of a pixel point at a position of (x,y) in a second chrominance component in a number t frame in the Sr; Ud(x,y,t) refers to a pixel value of a pixel point at a position of (x,y) in a number t frame in the Ud, which is also a pixel value of a pixel point at a position of (x,y) in a first chrominance component in a number t frame in the Sd; Vd(x,y,t) refers to a pixel value of a pixel point at a position of (x,y) in a number t frame in the Vd, which is also a pixel value of a pixel point at a position of (x,y) in a second chrominance component in a number t frame in the Sd; c2 and c3 are positive constants preventing the fractional from being meaningless;

(5) calculating a spatiotemporal domain local similarity between each pixel point in each frame in the Sr and the corresponding pixel point in the corresponding frame in the Sd according to the spatiotemporal domain local gradient similarity between each pixel point in each frame in the Sr and the corresponding pixel point in the corresponding frame in the Sd, and the spatiotemporal domain local color similarity between each pixel point in each frame in the Sr and the corresponding pixel point in the corresponding frame in the Sd; marking the spatiotemporal domain local similarity between the pixel point at the position of (x,y) in the number t frame in the Sr and the pixel point at the position of (x,y) in the number t frame in the Sd as QLS(x,y,t), wherein QLS(x,y,t)=Gsim(x,y,t)×(Csim(x,y,t))λ, λ is used for adjusting weights of color components, λ>0;
(6) calculating an objective quality value of each frame in the Sd with a variance fusion method, and marking the objective quality value of the number t frame in the Sd as Qframe(t), wherein




  
    
      
        Q
        frame
      
      ⁡
      
        (
        t
        )
      
    
    =
    
      
        
          1
          
            W
            ×
            H
          
        
        ⁢
        
          
            ∑
            
              x
              =
              1
            
            W
          
          ⁢
          
            
              ∑
              
                y
                =
                1
              
              H
            
            ⁢
            
              
            
            ⁢
            
              
                (
                
                  
                    
                      Q
                      LS
                    
                    ⁡
                    
                      (
                      
                        x
                        ,
                        y
                        ,
                        t
                      
                      )
                    
                  
                  -
                  
                    
                      Q
                      mean
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                )
              
              2
            
          
        
      
    
  
  ,




Qmean(t) refers to an average value of the spatiotemporal domain local similarity between all pixel points in the number t frame in the Sr and all pixel points in the number t frame in the Sd,





  
    
      
        Q
        mean
      
      ⁡
      
        (
        t
        )
      
    
    =
    
      
        1
        
          W
          ×
          H
        
      
      ⁢
      
        
          ∑
          
            x
            =
            1
          
          W
        
        ⁢
        
          
            ∑
            
              y
              =
              1
            
            H
          
          ⁢
          
            
          
          ⁢
          
            
              Q
              LS
            
            ⁡
            
              (
              
                x
                ,
                y
                ,
                t
              
              )
            
          
        
      
    
  
  ;




and

(7) calculating an objective quality value of the Sd with a temporal domain weighting method and marking as Q, wherein




  
    Q
    =
    
      
        1
        F
      
      ⁢
      
        
          ∑
          
            t
            =
            1
          
          F
        
        ⁢
        
          
        
        ⁢
        
          (
          
            
              
                Q
                LP
              
              ⁡
              
                (
                t
                )
              
            
            ×
            
              log
              ⁡
              
                (
                
                  
                    γ
                    ×
                    t
                  
                  +
                  1
                
                )
              
            
          
          )
        
      
    
  
  ,
  
    

  
  ⁢
  
    
      
        Q
        LP
      
      ⁡
      
        (
        t
        )
      
    
    =
    
      {
      
        
          
            
              
                
                  
                    Q
                    frame
                  
                  ⁡
                  
                    (
                    t
                    )
                  
                
                ,
                
                  
                    if
                    ⁢
                    
                      
                    
                    ⁢
                    t
                  
                  =
                  1
                
              
            
          
          
            
              
                
                  
                    
                      Q
                      LP
                    
                    ⁡
                    
                      (
                      
                        t
                        -
                        1
                      
                      )
                    
                  
                  +
                  
                    α
                    ×
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        Q
                        frame
                      
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                  
                
                ,
                
                  
                    if
                    ⁢
                    
                      
                    
                    ⁢
                    2
                  
                  ≤
                  t
                  ≤
                  
                    F
                    ⁢
                    
                      
                    
                    ⁢
                    and
                    ⁢
                    
                      
                    
                    ⁢
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        Q
                        frame
                      
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                  
                  ≤
                  0
                
              
            
          
          
            
              
                
                  
                    
                      Q
                      LP
                    
                    ⁡
                    
                      (
                      
                        t
                        -
                        1
                      
                      )
                    
                  
                  +
                  
                    β
                    ×
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        Q
                        frame
                      
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                  
                
                ,
                
                  
                    if
                    ⁢
                    
                      
                    
                    ⁢
                    2
                  
                  ≤
                  t
                  ≤
                  
                    F
                    ⁢
                    
                      
                    
                    ⁢
                    and
                    ⁢
                    
                      
                    
                    ⁢
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        Q
                        frame
                      
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                  
                  >
                  0
                
              
            
          
        
        ,
      
    
  




ΔQframe(t)=Qframe(t)−QLP(t−1), α refers to a weight during quality increase and β refers to a weight during quality decrease, γ is used to adjust a strength of a recency effect.

In the step (3), λ=3.
In the step (7), α=0.03, β=0.2, γ=1000.
Compared with conventional technologies, the present invention has advantages as follows.
1) Spatiotemporal domain gradient is able to reflect both spatial and temporal distortions, and color information in the video quality assessment is very important. The method of the present invention combines the spatiotemporal domain gradient magnitude and the color information for calculating the spatiotemporal domain local similarity between each pixel point in each frame in the reference video sequence without distortion and the corresponding pixel point in the corresponding frame in the distorted video sequence, and reflects the spatial and the temporal distortions at the same time with such bottom features, avoiding a widely used method which is based on motion estimation. Therefore, the present invention is accurate and fast.
2) According to the present invention, spatial domain fusion is provided by variance fusion. The spatiotemporal domain local similarity is fused into frame-level objective quality value, and then a temporal domain fusion model is established to fuse the frame-level objective quality values by simulating three important global temporal effects, which are a smoothing effect, an asymmetric track effects and a recency effect, of a human visual system. Finally, the objective quality values of the distorted video sequence are obtained. By modeling the human visual temporal domain effect, the temporal domain weighting method of the present invention is able to accurately and efficiently evaluate the objective quality of the distorted video.
3) During measuring the temporal domain distortion, the present invention firstly reflects local temporal domain distortion through the spatiotemporal domain local gradient similarity, then processes local temporal domain distortions with global integration through weighting by simulating the global temporal effect during temporal domain fusion, so as to avoid widely used temporal domain distortion assessment with motion vectors, and avoid time-consuming motion estimation, which ensures efficiency and low computational complexity of the present invention.
These and other objectives, features, and advantages of the present invention will become apparent from the following detailed description, the accompanying drawings, and the appended claims.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is an overall flow chart of the present invention.
FIG. 2 illustrates a horizontal mask, a vertical mask and a temporal mask of a three-dimensional Prewitt operator.



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT
Referring to drawings and a preferred embodiment, the present invention is further illustrated.
Outstanding full-reference video quality assessment methods not only are highly accurate in prediction and sensitive to both spatial and temporal distortions, but also lower a computational complexity as much as possible and provide real-time processing of video sequences. The present invention combines a spatiotemporal domain gradient magnitude and color information for calculating a spatiotemporal domain local similarity, and then treats a whole video sequence as a dynamic temporal sequence. Firstly, frame-level objective quality value of each frame is obtained by variance fusion, then a temporal domain fusion model is established by simulating three important global temporal effects, which are a smoothing effect, an asymmetric track effects and a recency effect, of a human visual system, so as to obtain objective quality values of a distorted video sequence. Bottom feature calculation is simple and sensitive to both temporal and spatial distortions, and the temporal domain fusion model analogs a temporal domain effect for ensuring accuracy and efficiency of the present invention.
The present invention provides a video quality objective assessment method based on a spatiotemporal domain structure, whose overall flow chart is shown in FIG. 1. The video quality objective assessment method comprises steps of:
(1) marking a reference video sequence without distortion as Sr, and marking a distorted video sequence, which is obtained after the Sr is distorted, as Sd; wherein a total Sr frame quantity is F, a total Sd frame quantity is also F, and F>1; widths of images in both the Sr and the Sd are W, and heights of the images in both the Sr and the Sd are H; defining an image luminance component sequence of the images in the Sr as a luminance component sequence of the Sr and marking as Yr; defining a first image chrominance component sequence of the images in the Sr as a first chrominance component sequence of the Sr and marking as Ur; defining a second image chrominance component sequence of the images in the Sr as a second chrominance component sequence of the Sr and marking as Vr; defining an image luminance component sequence of the images in the Sd as a luminance component sequence of the Sd and marking as Yd; defining a first image chrominance component sequence of the images in the Sd as a first chrominance component sequence of the Sd and marking as Ud; defining a second image chrominance component sequence of the images in the Sd as a second chrominance component sequence of the Sd and marking as Vd; wherein widths of images in the Yr, the Ur, the Vr, the Yd, the Ud and the Vd are W, and heights of the images in the Yr, the Ur, the Vr, the Yd, the Ud and the Vd are H;
(2) calculating a spatiotemporal domain gradient magnitude sequence of the Yr with a three-dimensional Prewitt operator and marking as Gr, and marking a pixel value of a pixel at a position of (x,y) in a number t frame in the Gr as Gr(x,y,t), wherein




  
    
      
        G
        r
      
      ⁡
      
        (
        
          x
          ,
          y
          ,
          t
        
        )
      
    
    =
    
      
        
          
            (
            
              
                Y
                rx
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          
            (
            
              
                Y
                ry
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          
            (
            
              
                Y
                rt
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
      
    
  
  ,




Yrx=Yr{circle around (×)}Fx, Yry=Yr{circle around (×)}Fy, Yrt=Yr{circle around (×)}Ft;

similarly, calculating a spatiotemporal domain gradient magnitude sequence of the Yd with the three-dimensional Prewitt operator and marking as Gd, and marking a pixel value of a pixel at a position of (x,y) in a number t frame in the Gd as Gd(x,y,t), wherein




  
    
      
        G
        d
      
      ⁡
      
        (
        
          x
          ,
          y
          ,
          t
        
        )
      
    
    =
    
      
        
          
            (
            
              
                Y
                dx
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          
            (
            
              
                Y
                dy
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          
            (
            
              
                Y
                dt
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
      
    
  
  ,




Ydx=Yd{circle around (×)}Fx, Ydy=Yd{circle around (×)}Fy, Ydt=Yd{circle around (×)}Ft;

wherein an initial value of the t is 1, 1≦t≦F, 1≦x≦W, 1≦y≦H; Yrx(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a horizontal gradient magnitude sequence Yrx of the Yr, Yry(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a vertical gradient magnitude sequence Yry of the Yr, and Yrt(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a temporal gradient magnitude sequence Yrt of the Yr; Ydx(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a horizontal gradient magnitude sequence Ydx of the Yd, Ydy(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a vertical gradient magnitude sequence Ydy of the Yd, and Ydt(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a temporal gradient magnitude sequence Ydt of the Yd; a symbol {circumflex over (×)} is a zeros truncated convolution symbol; after convolution, dimensions of the Yrx, the Yry and the Yrt are same as a dimension of the Yr, and dimensions of the Ydx, Ydy and Ydt are same as a dimension of the Yd; Fx, Fy and Ft correspond to a horizontal mask, a vertical mask and a temporal mask of the three-dimensional Prewitt operator; the Fx, the Fy and the Ft are shown in FIG. 2, and the three-dimensional Prewitt operator is available based on a two-dimensional Prewitt operator through conventional methods;
(3) calculating a spatiotemporal domain local gradient similarity between each pixel point in each frame in the Sr and a corresponding pixel point in a corresponding frame in the Sd; marking the spatiotemporal domain local gradient similarity between a pixel point at a position of (x,y) in a number t frame in the Sr and a pixel point at a position of (x,y) in a number t frame in the Sd as Gsim(x,y,t); wherein




  
    
      
        G
        sim
      
      ⁡
      
        (
        
          x
          ,
          y
          ,
          t
        
        )
      
    
    =
    
      
        
          2
          ⁢
          
            
          
          ⁢
          
            
              G
              r
            
            ⁡
            
              (
              
                x
                ,
                y
                ,
                t
              
              )
            
          
          ×
          
            
              G
              d
            
            ⁡
            
              (
              
                x
                ,
                y
                ,
                t
              
              )
            
          
        
        +
        
          c
          1
        
      
      
        
          
            (
            
              
                G
                r
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          
            (
            
              
                G
                d
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          c
          1
        
      
    
  
  ,




c1 is a positive constant preventing the fractional from being meaningless; according the preferred embodiment, c1=90;

(4) calculating a spatiotemporal domain local color similarity between each pixel point in each frame in the Sr and the corresponding pixel point in the corresponding frame in the Sd; marking the spatiotemporal domain local color similarity between the pixel point at the position of (x,y) in the number t frame in the Sr and the pixel point at the position of (x,y) in the number t frame in the Sd as Csim(x,y,t); wherein




  
    
      
        C
        sim
      
      ⁡
      
        (
        
          x
          ,
          y
          ,
          t
        
        )
      
    
    =
    
      
        
          
            2
            ⁢
            
              
                U
                r
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            ×
            
              
                U
                d
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
          
          +
          
            c
            2
          
        
        
          
            
              (
              
                
                  U
                  r
                
                ⁡
                
                  (
                  
                    x
                    ,
                    y
                    ,
                    t
                  
                  )
                
              
              )
            
            2
          
          +
          
            
              (
              
                
                  U
                  d
                
                ⁡
                
                  (
                  
                    x
                    ,
                    y
                    ,
                    t
                  
                  )
                
              
              )
            
            2
          
          +
          
            c
            2
          
        
      
      ×
      
        
          
            2
            ⁢
            
              
                V
                r
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            ×
            
              
                V
                d
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
          
          +
          
            c
            3
          
        
        
          
            
              (
              
                
                  V
                  r
                
                ⁡
                
                  (
                  
                    x
                    ,
                    y
                    ,
                    t
                  
                  )
                
              
              )
            
            2
          
          +
          
            
              (
              
                
                  V
                  d
                
                ⁡
                
                  (
                  
                    x
                    ,
                    y
                    ,
                    t
                  
                  )
                
              
              )
            
            2
          
          +
          
            c
            3
          
        
      
    
  
  ,




Ur(x,y,t) refers to a pixel value of a pixel point at a position of (x,y) in a number t frame in the Ur, which is also a pixel value of a pixel point at a position of (x,y) in a first chrominance component in a number t frame in the Sr; Vr(x,y,t) refers to a pixel value of a pixel point at a position of (x,y) in a number t frame in the Vr, which is also a pixel value of a pixel point at a position of (x,y) in a second chrominance component in a number t frame in the Sr; Ud(x,y,t) refers to a pixel value of a pixel point at a position of (x,y) in a number t frame in the Ud, which is also a pixel value of a pixel point at a position of (x,y) in a first chrominance component in a number t frame in the Sd; Vd(x,y,t) refers to a pixel value of a pixel point at a position of (x,y) in a number t frame in the Vd, which is also a pixel value of a pixel point at a position of (x,y) in a second chrominance component in a number t frame in the Sd; c2 and c3 are positive constants preventing the fractional from being meaningless; according to the preferred embodiment, c2=c3=300;

(5) calculating a spatiotemporal domain local similarity between each pixel point in each frame in the Sr and the corresponding pixel point in the corresponding frame in the Sd according to the spatiotemporal domain local gradient similarity between each pixel point in each frame in the Sr and the corresponding pixel point in the corresponding frame in the Sd, and the spatiotemporal domain local color similarity between each pixel point in each frame in the Sr and the corresponding pixel point in the corresponding frame in the Sd; marking the spatiotemporal domain local similarity between the pixel point at the position of (x,y) in the number t frame in the Sr and the pixel point at the position of (x,y) in the number t frame in the Sd as QLS(x,y,t), wherein QLS(x,y,t)=Gsim(x,y,t)×(Csim(x,y,t))λ, λ is used for adjusting weights of color components, λ>0; according to the preferred embodiment, λ=3;
(6) calculating an objective quality value of each frame in the Sd with a variance fusion method, and marking the objective quality value of the number t frame in the Sd as Qframe(t), wherein




  
    
      
        Q
        frame
      
      ⁡
      
        (
        t
        )
      
    
    =
    
      
        
          1
          
            W
            ×
            H
          
        
        ⁢
        
          
            ∑
            
              x
              =
              1
            
            W
          
          ⁢
          
            
          
          ⁢
          
            
              ∑
              
                y
                =
                1
              
              H
            
            ⁢
            
              
            
            ⁢
            
              
                (
                
                  
                    
                      Q
                      LS
                    
                    ⁡
                    
                      (
                      
                        x
                        ,
                        y
                        ,
                        t
                      
                      )
                    
                  
                  -
                  
                    
                      Q
                      mean
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                )
              
              2
            
          
        
      
    
  
  ,




Qmean(t) refers to an average value of the spatiotemporal domain local similarity between all pixel points in the number t frame in the Sr and all pixel points in the number t frame in the Sd,





  
    
      
        Q
        mean
      
      ⁡
      
        (
        t
        )
      
    
    =
    
      
        1
        
          W
          ×
          H
        
      
      ⁢
      
        
          ∑
          
            x
            =
            1
          
          W
        
        ⁢
        
          
        
        ⁢
        
          
            ∑
            
              y
              =
              1
            
            H
          
          ⁢
          
            
          
          ⁢
          
            
              Q
              LS
            
            ⁡
            
              (
              
                x
                ,
                y
                ,
                t
              
              )
            
          
        
      
    
  
  ;




and

(7) calculating an objective quality value of the Sd with a temporal domain weighting method and marking as Q, wherein




  
    Q
    =
    
      
        1
        F
      
      ⁢
      
        
          ∑
          
            t
            =
            1
          
          F
        
        ⁢
        
          
        
        ⁢
        
          (
          
            
              
                Q
                LP
              
              ⁡
              
                (
                t
                )
              
            
            ×
            
              log
              ⁡
              
                (
                
                  
                    γ
                    ×
                    t
                  
                  +
                  1
                
                )
              
            
          
          )
        
      
    
  
  ,
  
    

  
  ⁢
  
    
      
        Q
        LP
      
      ⁡
      
        (
        t
        )
      
    
    =
    
      {
      
        
          
            
              
                
                  
                    Q
                    frame
                  
                  ⁡
                  
                    (
                    t
                    )
                  
                
                ,
              
            
            
              
                
                  if
                  ⁢
                  
                    
                  
                  ⁢
                  t
                
                =
                1
              
            
          
          
            
              
                
                  
                    
                      Q
                      LP
                    
                    ⁡
                    
                      (
                      
                        t
                        -
                        1
                      
                      )
                    
                  
                  +
                  
                    α
                    ×
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        Q
                        frame
                      
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                  
                
                ,
              
            
            
              
                
                  if
                  ⁢
                  
                    
                  
                  ⁢
                  2
                
                ≤
                t
                ≤
                
                  F
                  ⁢
                  
                    
                  
                  ⁢
                  and
                  ⁢
                  
                    
                  
                  ⁢
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      Q
                      frame
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                ≤
                0
              
            
          
          
            
              
                
                  
                    
                      Q
                      LP
                    
                    ⁡
                    
                      (
                      
                        t
                        -
                        1
                      
                      )
                    
                  
                  +
                  
                    β
                    ×
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        Q
                        frame
                      
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                  
                
                ,
              
            
            
              
                
                  if
                  ⁢
                  
                    
                  
                  ⁢
                  2
                
                ≤
                t
                ≤
                
                  F
                  ⁢
                  
                    
                  
                  ⁢
                  and
                  ⁢
                  
                    
                  
                  ⁢
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      Q
                      frame
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                >
                0
              
            
          
        
        ,
      
    
  




ΔQframe(t)=Qframe(t)−QLP(t−1), α refers to a weight during quality increase and β refers to a weight during quality decrease, γ is used to adjust a strength of a recency effect; according to the preferred embodiment, α=0.03, β=0.2, γ=1000.

For illustrating feasibility and effectiveness, the present invention is tested as follows.
A LIVE video database and a CSIQ video database are used. Referring to the LIVE video database, there are 10 video segments without distortion; four distortion types are involved, which are MPEG-2 compression distortion, H.264 compression distortion, transmission distortion of bit-stream compressed with H.264 standard through IP network, and transmission distortion through wireless network; there are totally 150 segments of distorted videos which are all YUV420 format with a resolution of 768×432 and a length of 10 seconds; and two frame rates, 25 frames per second and 50 frames per second, are used. Referring to the CSIQ video database, there are 12 video segments without distortion; six distortion types are involved, which are motion JPEG compression distortion, H.264 compression distortion, HEVC compression distortion, wavelet compression distortion, loss distortion of wireless transmission, and additive white Gaussian noise distortion; there are totally 216 segments of distorted videos which are all YUV420 format with a resolution of 832×480 and a length of 10 seconds; and five frame rates, 24 frames per second, 25 frames per second, 30 frames per second, 50 frames per second and 60 frames per second, are used. Both the LIVE video database and the CSIQ video database provide average subjective opinion scores of each distorted video. Since two chrominance components are both 1⁄4 of a luminance component in the YUV420 format, the luminance component is processed with 2×2 mean filter in the spatial domain and down-sampling by a factor of 2, so as to matching dimensions of luminance component and chrominance components. Then the steps (1) to (7) are executed, and objective quality values of all distorted videos are calculated in a same way. For other formats such as YUV444 and YUYV, since the present invention requires that the dimension of the luminance component equal to the dimensions of the first chrominance component and the second chrominance component, down-sampling or up-sampling is also needed, so as to matching dimensions of the luminance component and the first as well as the second chrominance components. According to the present invention, in each treated video, sizes of all frames are equal.
Three common objective parameters of the video quality assessment method are used as evaluation criteria, which are Pearson linear correlation coefficient (PLCC), Spearman rank order correlation coefficient (SROCC) and rooted mean squared error (RMSE). Valve ranges of PLCC and SROCC are [0,1], wherein if a value is closer to 1, the assessment method is better; otherwise, the assessment method is worse. If a RMSE value is smaller, the assessment method is better; otherwise, the assessment method is worse. The PLCC, the SROCC and the RMSE values indicating assessment performances of the LIVE video database are shown in a Table 1, and the PLCC, the SROCC and the RMSE values indicating assessment performances of the CSIQ video database are shown in a Table 2. Referring to Table 1, values of the PLCC and the SROCC are higher than 0.84. Referring to Table 2, values of the PLCC and the SROCC are higher than 0.80. That is to say, with the present invention, the objective quality values calculated is highly related to the average subjective opinion scores, indicating that objective assessment results are relatively consistent with subjective perception of human eyes, which illustrates the effectiveness the present invention.







TABLE 1







relativity between objective quality values of distorted


videos obtained by the present invention and average


subjective opinion scores for LIVE video database















transmission
trans-






distortion of
mission



MPEG-2
H.264
H.264 compres-
distortion
all



compres-
compres-
sed bit flow
through
dis-



sion dis-
sion dis-
through IP
wireless
torted



tortion
tortion
network
network
videos
















PLCC
0.9115
0.8260
0.8401
0.8816
0.8632


SROCC
0.9122
0.7660
0.8126
0.8332
0.8475


RMSE
4.6162
6.3169
5.4474
5.3174
5.5439
















TABLE 2







relativity between objective quality values of distorted videos obtained by the


present invention and average subjective opinion scores for CSIQ video database




















additive




motion



loss
white



JPEG
H.264
HEVC
wavelet
distortion of
Gaussian
all



compression
compression
compression
compression
wireless
noise
distorted



distortion
distortion
distortion
distortion
transmission
distortion
videos


















PLCC
0.9366
0.7976
0.8590
0.8748
0.9030
0.9586
0.8099


SROCC
0.9331
0.7985
0.8474
0.8440
0.8723
0.9416
0.8302


RMSE
6.9568
6.3461
7.6241
7.3420
4.6853
6.4334
9.7535









One skilled in the art will understand that the embodiment of the present invention as shown in the drawings and described above is exemplary only and not intended to be limiting.
It will thus be seen that the objects of the present invention have been fully and effectively accomplished. Its embodiments have been shown and described for the purposes of illustrating the functional and structural principles of the present invention and is subject to change without departure from such principles. Therefore, this invention includes all modifications encompassed within the spirit and scope of the following claims.

## Claims

1. A video quality objective assessment method based on a spatiotemporal domain structure, comprising steps of:
(1) marking a reference video sequence without distortion as Sr, and marking a distorted video sequence, which is obtained after the Sr is distorted, as Sd; wherein a total Sr frame quantity is F, a total Sd frame quantity is also F, and F>1; widths of images in both the Sr and the Sd are W, and heights of the images in both the Sr and the Sd are H; defining an image luminance component sequence of the images in the Sr as a luminance component sequence of the Sr and marking as Yr; defining a first image chrominance component sequence of the images in the Sr as a first chrominance component sequence of the Sr and marking as Ur; defining a second image chrominance component sequence of the images in the Sr as a second chrominance component sequence of the Sr and marking as Vr; defining an image luminance component sequence of the images in the Sd as a luminance component sequence of the Sd and marking as Yd; defining a first image chrominance component sequence of the images in the Sd as a first chrominance component sequence of the Sd and marking as Ud; defining a second image chrominance component sequence of the images in the Sd as a second chrominance component sequence of the Sd and marking as Vd; wherein widths of images in the Yr, the Ur, the Vr, the Yd, the Ud and the Vd are W, and heights of the images in the Yr, the Ur, the Vr, the Yd, the Ud and the Vd are H;
(2) calculating a spatiotemporal domain gradient magnitude sequence of the Yr with a three-dimensional Prewitt operator and marking as Gr, and marking a pixel value of a pixel at a position of (x,y) in a number t frame in the Gr as Gr(x,y,t), wherein





  
    
      
        G
        r
      
      ⁡
      
        (
        
          x
          ,
          y
          ,
          t
        
        )
      
    
    =
    
      
        
          
            (
            
              
                Y
                rx
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          
            (
            
              
                Y
                ry
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          
            (
            
              
                Y
                rt
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
      
    
  
  ,




Yrx=Yr{circle around (×)}Fx, Yry=Yr{circle around (×)}Fy, Yrt=Yr{circle around (×)}Ft,
similarly, calculating a spatiotemporal domain gradient magnitude sequence of the Yd with the three-dimensional Prewitt operator and marking as Gd, and marking a pixel value of a pixel at a position of (x,y) in a number t frame in the Gd as Gd(x,y,t), wherein





  
    
      
        G
        d
      
      ⁡
      
        (
        
          x
          ,
          y
          ,
          t
        
        )
      
    
    =
    
      
        
          
            (
            
              
                Y
                dx
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          
            (
            
              
                Y
                dy
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          
            (
            
              
                Y
                dt
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
      
    
  
  ,




Ydx=Yd{circle around (×)}Fx, Ydy=Yd{circle around (×)}Fy, Ydt=Yd{circle around (×)}Ft;
wherein an initial value of the t is 1, 1≦t≦F, 1≦x≦W, 1≦y≦H; Yrx(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a horizontal gradient magnitude sequence Yrx of the Yr, Yry(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a vertical gradient magnitude sequence Yry of the Yr, and Yrrt(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a temporal gradient magnitude sequence Yrt of the Yr; Ydx(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a horizontal gradient magnitude sequence Ydx of the Yd, Ydy(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a vertical gradient magnitude sequence Ydy of the Yd, and Ydt(x,y,t) refers to a pixel value of a pixel at a position of (x,y) in a number t frame in a temporal gradient magnitude sequence Ydt of the Yd; a symbol {circle around (×)} is a zeros truncated convolution symbol; after convolution, dimensions of the Yrx, the Yry and the Yrt are same as a dimension of the Yr, and dimensions of the Ydx, Ydy and Ydt are same as a dimension of the Yd; Fx, Fy and Ft correspond to a horizontal mask, a vertical mask and a temporal mask of the three-dimensional Prewitt operator;
(3) calculating a spatiotemporal domain local gradient similarity between each pixel point in each frame in the Sr and a corresponding pixel point in a corresponding frame in the Sd; marking the spatiotemporal domain local gradient similarity between a pixel point at a position of (x,y) in a number t frame in the Sr and a pixel point at a position of (x,y) in a number t frame in the Sd as Gsim(x,y,t); wherein





  
    
      
        G
        sim
      
      ⁡
      
        (
        
          x
          ,
          y
          ,
          t
        
        )
      
    
    =
    
      
        
          2
          ⁢
          
            
              G
              r
            
            ⁡
            
              (
              
                x
                ,
                y
                ,
                t
              
              )
            
          
          ×
          
            
              G
              d
            
            ⁡
            
              (
              
                x
                ,
                y
                ,
                t
              
              )
            
          
        
        +
        
          c
          1
        
      
      
        
          
            (
            
              
                G
                r
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          
            (
            
              
                G
                d
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            )
          
          2
        
        +
        
          c
          1
        
      
    
  
  ,




c1 is a positive constant preventing the fractional from being meaningless;
(4) calculating a spatiotemporal domain local color similarity between each pixel point in each frame in the Sr and the corresponding pixel point in the corresponding frame in the Sd; marking the spatiotemporal domain local color similarity between the pixel point at the position of (x,y) in the number t frame in the Sr and the pixel point at the position of (x,y) in the number t frame in the Sd as Csim(x,y,t); wherein





  
    
      
        C
        sim
      
      ⁡
      
        (
        
          x
          ,
          y
          ,
          t
        
        )
      
    
    =
    
      
        
          
            2
            ⁢
            
              
                U
                r
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            ×
            
              
                U
                d
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
          
          +
          
            c
            2
          
        
        
          
            
              (
              
                
                  U
                  r
                
                ⁡
                
                  (
                  
                    x
                    ,
                    y
                    ,
                    t
                  
                  )
                
              
              )
            
            2
          
          +
          
            
              (
              
                
                  U
                  d
                
                ⁡
                
                  (
                  
                    x
                    ,
                    y
                    ,
                    t
                  
                  )
                
              
              )
            
            2
          
          +
          
            c
            2
          
        
      
      ×
      
        
          
            2
            ⁢
            
              
                V
                r
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
            ×
            
              
                V
                d
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
          
          +
          
            c
            3
          
        
        
          
            
              (
              
                
                  V
                  r
                
                ⁡
                
                  (
                  
                    x
                    ,
                    y
                    ,
                    t
                  
                  )
                
              
              )
            
            2
          
          +
          
            
              (
              
                
                  V
                  d
                
                ⁡
                
                  (
                  
                    x
                    ,
                    y
                    ,
                    t
                  
                  )
                
              
              )
            
            2
          
          +
          
            c
            3
          
        
      
    
  
  ,




Ur(x,y,t) refers to a pixel value of a pixel point at a position of (x,y) in a number t frame in the Ur, which is also a pixel value of a pixel point at a position of (x,y) in a first chrominance component in a number t frame in the Sr; Vr(x,y,t) refers to a pixel value of a pixel point at a position of (x,y) in a number t frame in the Vr, which is also a pixel value of a pixel point at a position of (x,y) in a second chrominance component in a number t frame in the Sr; Ud(x,y,t) refers to a pixel value of a pixel point at a position of (x,y) in a number t frame in the Ud, which is also a pixel value of a pixel point at a position of (x,y) in a first chrominance component in a number t frame in the Sd; Vd(x,y,t) refers to a pixel value of a pixel point at a position of (x,y) in a number t frame in the Vd, which is also a pixel value of a pixel point at a position of (x,y) in a second chrominance component in a number t frame in the Sd; c2 and c3 are positive constants preventing the fractional from being meaningless;
(5) calculating a spatiotemporal domain local similarity between each pixel point in each frame in the Sr and the corresponding pixel point in the corresponding frame in the Sd according to the spatiotemporal domain local gradient similarity between each pixel point in each frame in the Sr and the corresponding pixel point in the corresponding frame in the Sd, and the spatiotemporal domain local color similarity between each pixel point in each frame in the Sr and the corresponding pixel point in the corresponding frame in the Sd; marking the spatiotemporal domain local similarity between the pixel point at the position of (x,y) in the number t frame in the Sr and the pixel point at the position of (x,y) in the number t frame in the Sd as QLS(x,y,t), wherein QLS(x,y,t)=Gsim(x,y,t)×(Csim(x,y,t))λ, λ is used for adjusting weights of color components, λ>0;
(6) calculating an objective quality value of each frame in the Sd with a variance fusion method, and marking the objective quality value of the number t frame in the Sd as Qframe(t), wherein





  
    
      
        Q
        frame
      
      ⁡
      
        (
        t
        )
      
    
    =
    
      
        
          1
          
            W
            ×
            H
          
        
        ⁢
        
          
            ∑
            
              x
              =
              1
            
            W
          
          ⁢
          
            
          
          ⁢
          
            
              ∑
              
                y
                =
                1
              
              H
            
            ⁢
            
              
            
            ⁢
            
              
                (
                
                  
                    
                      Q
                      LS
                    
                    ⁡
                    
                      (
                      
                        x
                        ,
                        y
                        ,
                        t
                      
                      )
                    
                  
                  -
                  
                    
                      Q
                      mean
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                )
              
              2
            
          
        
      
    
  
  ,




Qmean(t) refers to an average value of the spatiotemporal domain local similarity between all pixel points in the number t frame in the Sr and all pixel points in the number t frame in the




  
    S
    d
  
  ,
  
    
      
        
          Q
          mean
        
        ⁡
        
          (
          t
          )
        
      
      =
      
        
          1
          
            W
            ×
            H
          
        
        ⁢
        
          
            ∑
            
              x
              =
              1
            
            W
          
          ⁢
          
            
          
          ⁢
          
            
              ∑
              
                y
                =
                1
              
              H
            
            ⁢
            
              
            
            ⁢
            
              
                Q
                LS
              
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  t
                
                )
              
            
          
        
      
    
    ;
  




and
(7) calculating an objective quality value of the Sd with a temporal domain weighting method and marking as Q, wherein





  
    Q
    =
    
      
        1
        F
      
      ⁢
      
        
          ∑
          
            t
            =
            1
          
          F
        
        ⁢
        
          
        
        ⁢
        
          (
          
            
              
                Q
                LP
              
              ⁡
              
                (
                t
                )
              
            
            ×
            
              log
              ⁡
              
                (
                
                  γ
                  +
                  t
                  +
                  1
                
                )
              
            
          
          )
        
      
    
  
  ,
  
    

  
  ⁢
  
    
      
        Q
        LP
      
      ⁡
      
        (
        t
        )
      
    
    =
    
      {
      
        
          
            
              
                
                  
                    Q
                    frame
                  
                  ⁡
                  
                    (
                    t
                    )
                  
                
                ,
              
            
            
              
                
                  if
                  ⁢
                  
                    
                  
                  ⁢
                  t
                
                =
                1
              
            
          
          
            
              
                
                  
                    
                      Q
                      LP
                    
                    ⁡
                    
                      (
                      
                        t
                        -
                        1
                      
                      )
                    
                  
                  +
                  
                    α
                    ×
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        Q
                        frame
                      
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                  
                
                ,
              
            
            
              
                
                  if
                  ⁢
                  
                    
                  
                  ⁢
                  2
                
                ≤
                t
                ≤
                
                  F
                  ⁢
                  
                    
                  
                  ⁢
                  and
                  ⁢
                  
                    
                  
                  ⁢
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      Q
                      frame
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                ≤
                0
              
            
          
          
            
              
                
                  
                    
                      Q
                      LP
                    
                    ⁡
                    
                      (
                      
                        t
                        -
                        1
                      
                      )
                    
                  
                  +
                  
                    β
                    ×
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        Q
                        frame
                      
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                  
                
                ,
              
            
            
              
                
                  if
                  ⁢
                  
                    
                  
                  ⁢
                  2
                
                ≤
                t
                ≤
                
                  F
                  ⁢
                  
                    
                  
                  ⁢
                  and
                  ⁢
                  
                    
                  
                  ⁢
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      Q
                      frame
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                >
                0
              
            
          
        
        ,
      
    
  




ΔQframe(t)=Qframe(t)−QLP(t−1), α refers to a weight during quality increase and β refers to a weight during quality decrease, γ is used to adjust a strength of a recency effect.

2. The video quality objective assessment method, as recited in claim 1, wherein in the step (3), λ=3.

3. The video quality objective assessment method, as recited in claim 1, wherein in the step (7), α=0.03, β=0.2, γ=1000.