# Patent

## Title
Objective assessment method for color image quality based on online manifold learning

## Abstract

An objective assessment method for a color image quality based on online manifold learning considers a relationship between a saliency and an image quality objective assessment. Through a visual saliency detection algorithm, saliency maps of a reference image and a distorted image are obtained for further obtaining a maximum fusion saliency map. Based on maximum saliencies of image blocks in the maximum fusion saliency map, a saliency difference between each reference image block and a corresponding distorted image block is measured through an absolute difference, and thus reference visual important image blocks and distorted visual important image blocks are screened and extracted. Through manifold eigenvectors of the reference visual important image blocks and the distorted visual important image blocks, an objective quality assessment value of the distorted image is calculated. The method has an increased assessment effect and a higher correlation between an objective assessment result and a subjective perception.

## Background
CROSS REFERENCE OF RELATED APPLICATION
The application claims priority under 35 U.S.C. 119(a-d) to CN 201610202181.5, filed Mar. 31, 2016.
BACKGROUND OF THE PRESENT INVENTION
Field of Invention
The present invention relates to an image quality assessment method, and more particularly to an objective assessment method for a color image quality based on online manifold learning.
Description of Related Arts
Because of the performance limitation of the image processing system, during image obtaining, transmission and coding, various kinds of the distortions are introduced. The introduction of the distortions decreases the image quality, and meanwhile has an adverse impact on obtaining the information from the image. The image quality is an important index for comparing the performances of the various image processing algorithms and the image processing system parameters. Thus, in the fields of image transmission, multi-media network communication and video analysis, it is important to construct an effective image quality assessment method. Conventionally, the image quality assessment methods are divided into the subjective assessment methods and the objective assessment methods. Because the final receiver of the image is human, the subjective assessment methods are the most reliable assessment methods. However, the subjective assessment methods are time-consuming, and difficult to be embedded into the image processing system, so that the subjective assessment methods are limited in the practical application. In comparison, the objective assessment methods have the advantages of the simple operation and the good practicability, and currently are the research emphasis in academia and even industry.
Conventionally, the easiest and the most widely applied objective assessment methods are the peak signal-to-noise ratio (PSNR) method and the mean square error (MSE) method. The two methods have the simple calculation and the definite physical significance. However, because the visual characteristics of the human eyes are not considered, the assessment result thereof generally does not conform to the subjective perception of the human eyes. In fact, the human eyes do not process the image signal point by point. Therefore, the researchers introduce the visual characteristics of the human eyes, so that the objective assessment result has a higher conformity with the visual perception of the human eyes. For example, a method based on the structural similarity (SSIM) represents the structural information of the image from the brightness, the contrast, and the structure of the image, so as to further assess the image quality. Subsequently, based on the SSIM method, the multi-scale SSIM assessment method, the complex wavelet SSIM assessment method and the SSIM assessment method based on the information content weighting are proposed, which improve the performance of the SSIM method. Besides the assessment methods based on the SSIM, Sheikh et al. see the full-reference image quality assessment as the information fidelity problem, and, according to the loss amount of the image information during the quantization distortion, propose an image quality assessment method based on the visual information fidelity (VIF). From the critical threshold and the above-threshold characteristics of the visual perception of the image, combined with the wavelet transform, Chandler et al. propose an image quality assessment method based on the wavelet visual signal-to-noise ratio (VSNR), which is able to relatively well adapt to the different visual conditions. Although the researchers have deeply researched the human visual system, because of the complexity of the human visual system, the researchers still have little cognition on the human visual system. Thus, it is still failed to propose an objective assessment method for the image quality, which is completely consistent with the subjective perception of the human eyes.
In order to better reflect the characteristics of the human visual system, the objective assessment methods for the image quality based on the sparse representation and the visual attention have gained more and more attentions. A number of the researches indicate that the sparse representation is able to well describe the neuronal activities in the primary visual cortex of the human brain. For example, Guha et al. disclose an image quality assessment method based on the sparse representation. The method is divided into two phases, respectively the dictionary learning phase and the assessment phase. The dictionary learning phase as the first phase is to randomly choose the image blocks from the reference image to serve as the training samples, and then obtain the over-complete dictionary through training with the k-means singular value decomposition (K-SVD) algorithm. The assessment phase as the second phase is to process the image blocks of the reference image and the image blocks of the corresponding distorted image with the sparse coding through the orthogonal matching pursuit (OMP) algorithm, then obtain the sparse coefficients of the reference image and the distorted image, and further obtain the objective image assessment value of the distorted image. However, all of the objective assessment methods for the image quality based on the sparse representation utilize the OMP algorithm to proceed the sparse coding, which needs the high computational cost. Moreover, the objective assessment methods for the image quality based on the sparse representation obtain the over-complete dictionary through the off-line operation, which needs a large number of the effective natural images to serve as the training samples.
For the digital image which is the high dimensional data, a large number of the information redundancies exist substantially, which requires to be processed through the dimensionality reduction technology. Moreover, while reducing the dimensionality, the essential structure thereof is expected to be maintained. Since the manifold learning method was firstly proposed in the famous science magazine, Science, in 2000, the manifold learning method has become the research hotspot in the information science field. If the data are assumed to be uniformly sampled from the low dimensional manifold in the high dimensional Euclidean space, the manifold learning is to recover the low dimensional structure from the high dimensional sampling data, namely to find the low dimensional manifold in the high dimensional space and obtain the corresponding embedded mapping, for realizing the dimensionality reduction. Some researches indicate that: the manifold is the basis of the perception, and the human brain perceives the objects in the manifold manner. In recent years, the manifold learning has been widely applied in image denoising, face recognition and human behavior detection, and achieved relatively good effects. For the problem that the column vector in the locality preserving projection (LPP) algorithm is not orthogonal, Deng et al. improve the LPP algorithm and obtain the orthogonal locality preserving projection (OLPP) algorithm which is able to find the manifold structure of the data, has the linear features, and achieves the good locality preserving ability and the judging ability. The manifold learning is able to simulate the description of the image signal in the primary visual cortex cells and further accurately extract the visual perception characteristics of the image. The low dimensional manifold characteristics of the image relatively well describe the nonlinear variation relationship among each distorted image, and the distorted images in the manifold space are arranged according to the variation type and intensity. Thus, it is necessary to research an objective assessment method for the image quality based on the manifold learning with the objective assessment result having the high conformity with the visual perception of the human eyes.
SUMMARY OF THE PRESENT INVENTION
An object of the present invention is to provide an objective assessment method for a color image quality based on online manifold learning, which is able to effectively increase a correlation between an objective assessment result and a subjective perception.
Technical solutions of the present invention are described as follows.
An objective assessment method for a color image quality based on online manifold learning comprises steps of:
{circle around (1)} representing an undistorted reference image having a width W and a height H by IR; and representing a distorted image to be assessed, which is corresponding to the IR, by ID;
{circle around (2)} through a visual saliency detection algorithm, respectively obtaining saliency maps of the IR and the ID, correspondingly denoted as MR and MD; then, according to the MR and the MD, calculating a maximum fusion saliency map, denoted as MF; and denoting a pixel value of a pixel having coordinates of (x, y) in the MF as MF(x, y), MF(x, y)=max(MR(x, y), MD(x, y)), wherein: 1≦x≦W, 1≦y≦H; the max( ) is a function to find a maximum; the MR(x, y) represents a pixel value of a pixel having coordinates of (x, y) in the MR; and the MD(x, y) represents a pixel value of a pixel having coordinates of (x, y) in the MD;
{circle around (3)} respectively dividing the IR, the ID, the MR, the MD, and the MF into




  
    W
    ×
    H
  
  
    8
    ×
    8
  




image blocks, which are not overlapping mutually and have a size of 8×8;

vectorizing color values of R, G, and B channels of all pixels in each image block of the IR and the ID; denoting a color vector obtained through vectorizing the color values of the R, the G, and the B channels of all the pixels in a jth image block of the IR as XjR; denoting a color vector obtained through vectorizing the color values of the R, the G, and the B channels of all the pixels in a jth image block of the ID as XjD; wherein: the j has an initial value of 1,




  
    1
    ≤
    j
    ≤
    
      
        W
        ×
        H
      
      
        8
        ×
        8
      
    
  
  ;




both of the XjR and the XjD have a dimensionality of 192×1; values of a 1st element to a 64th element in the XjR respectively correspond to the color value of the R channel of each pixel in the jth image block of the IR in a line-by-line scanning manner; values of a 65th element to a 128th element in the XjR respectively correspond to the color value of the G channel of each pixel in the jth image block of the IR in the line-by-line scanning manner; values of a 129th element to a 192nd element in the XjR respectively correspond to the color value of the B channel of each pixel in the jth image block of the IR in the line-by-line scanning manner; values of a 1st element to a 64th element in the XjD respectively correspond to the color value of the R channel of each pixel in the jth image block of the ID in the line-by-line scanning manner; values of a 65th element to a 128th element in the XjD respectively correspond to the color value of the G channel of each pixel in the jth image block of the ID in the line-by-line scanning manner; and values of a 129th element to a 192nd element in the XjD respectively correspond to the color value of the B channel of each pixel in the jth image block of the ID in the line-by-line scanning manner; and

vectorizing pixel values of all pixels in each image block of the MR, the MD, and the MF; denoting a pixel value vector obtained through vectorizing the pixel values of all the pixels in a jth image block of the MF as SjR; denoting a pixel value vector obtained through vectorizing the pixel values of all the pixels in a jth image block of the MD as SjD; and denoting a pixel value vector obtained through vectorizing the pixel values of all the pixels in a jth image block of the MF as SjF; wherein: the SjR, the SjD, and the SjF all have a dimensionality of 64×1; values of a 1st element to a 64th element in the SjR respectively correspond to the pixel value of each pixel in the jth image block of the MR in the line-by-line scanning manner; values of a 1st element to a 64th element in the SjD respectively correspond to the pixel value of each pixel in the jth image block of the MD in the line-by-line scanning manner; and values of a 1st element to a 64th element in the SjF respectively correspond to the pixel value of each pixel in the jth image block of the MF in the line-by-line scanning manner;
{circle around (4)} calculating a saliency of each image block in the MF; and denoting the saliency of the jth image block in the MF as dj,




  
    
      d
      j
    
    =
    
      
        ∑
        
          i
          =
          1
        
        64
      
      ⁢
      
        
          S
          j
          F
        
        ⁡
        
          (
          i
          )
        
      
    
  
  ,




wherein: 1≦i≦64 and the SjF(i) represents a value of an ith element in the SjF;

orderly arranging the saliencies of all the image blocks in the MF, from the biggest to the smallest; and, after arranging, determining sequence numbers of the image blocks corresponding to former t1 saliencies, wherein:




  
    
      t
      1
    
    =
    
      
        λ
        1
      
      ×
      
        
          W
          ×
          H
        
        
          8
          ×
          8
        
      
    
  
  ;




the λ1 represents an image block selection proportionality coefficient; and λ1ε(0,1]; and

finding the image blocks in the IR, which are corresponding to the determined t1 sequence numbers, and defining as reference image blocks; finding the image blocks in the ID, which are corresponding to the determined t1 sequence numbers, and defining as distorted image blocks; finding the image blocks in the MR, which are corresponding to the determined t1 sequence numbers, and defining as reference saliency image blocks; finding the image blocks in the MD, which are corresponding to the determined t1 sequence numbers, and defining as distorted saliency image blocks;
{circle around (5)} measuring a saliency difference between each reference image block in the IR and a corresponding distorted image block in the ID through an absolute difference; and denoting a saliency difference between a t′th reference image block in the IR and a t′th distorted image block in the ID as et′,




  
    
      e
      
        t
        ′
      
    
    =
    
      
        1
        64
      
      ⁢
      
        
          ∑
          
            i
            =
            1
          
          64
        
        ⁢
        
          
          
            
              
                
                  S
                  ~
                
                
                  t
                  ′
                
                R
              
              ⁡
              
                (
                i
                )
              
            
            -
            
              
                
                  S
                  ~
                
                
                  t
                  ′
                
                D
              
              ⁡
              
                (
                i
                )
              
            
          
          
        
      
    
  
  ,




wherein: the t′ has an initial value of 1, 1≦t′≦t1; the symbol “∥” is an absolute value symbol; the {tilde over (S)}t′R(i) represents a value of an ith element in a pixel value vector {tilde over (S)}t′R corresponding to a t′th reference saliency image block in the MR; and the {tilde over (S)}t′D(i) represents a value of an ith element in a pixel value vector {tilde over (S)}t′D corresponding to a t′th distorted saliency image block in the MD; and

orderly arranging the measured t1 saliency differences, from the biggest to the smallest; after arranging, determining the reference image blocks and the distorted image blocks corresponding to former t2 saliency differences; defining the determined t2 reference image blocks as reference visual important image blocks, and adopting a matrix formed by color vectors corresponding to all the reference visual important image blocks as a reference visual important image block matrix, denoted as YR; defining the determined t2 distorted image blocks as distorted visual important image blocks, and adopting a matrix formed by color vectors corresponding to all the distorted visual important image blocks as a distorted visual important image block matrix, denoted as YD; wherein: t2=λ2×t1, the λ2 represents a selection proportionality coefficient of the reference image blocks and the distorted image blocks, and λ2ε(0,1]; the YR and the YD have a dimensionality of 192×t2; a t′′th column vector in the YR is a color vector corresponding to a determined t′′th reference visual important image block; a t′′th column vector in the YD is a color vector corresponding to a determined t′′th distorted visual important image block; and the t has an initial value of 1, 1≦t′′≦t2;
{circle around (6)} centralizing the YR through subtracting a mean value of values of all elements in each column vector from a value of each element in the same column vector of the YR; and denoting an obtained centralized matrix as Y, wherein the Y has a dimensionality of 192×t2; and
processing the Y with dimensionality reduction and whitening through a principal component analysis; and, denoting an obtained matrix after the dimensionality reduction and the whitening as Yw, Yw=W×Y; wherein: the Yw has a dimensionality of M×t2; the W represents a whitened matrix and has a dimensionality of M×192, 1<M<<192; and the symbol “<<” is a much-less-than symbol;
{circle around (7)} online training the Yw through an orthogonal locality preserving projection algorithm; and obtaining a characteristic basis matrix of the Yw, denoted as D, wherein the D has a dimensionality of M×192;
{circle around (8)} according to the YR and the D, calculating a manifold eigenvector of each reference visual important image block; denoting the manifold eigenvector of a t′′th reference visual important image block as ut′′, ut′′=D×yt′′R, wherein: the ut′′ has a dimensionality of M×1, and the yt′′R, is the t′′th column vector in the YR; according to the YD and the D, calculating a manifold eigenvector of each distorted visual important image block; and denoting the manifold eigenvector of a t′′th distorted visual important image block as vt′′, vt′′=D×yt′′D, wherein: the vt′′ has a dimensionality of M×1, and the yt′′D is the t′′th column vector in the YD; and
{circle around (9)} according to the manifold eigenvectors of all the reference visual important image blocks and the manifold eigenvectors of all the distorted visual important image blocks, calculating an objective quality assessment value of the ID, denoted as Score,




  
    Score
    =
    
      
        1
        
          
            t
            2
          
          ×
          M
        
      
      ⁢
      
        
          ∑
          
            
              t
              ′′
            
            =
            1
          
          
            t
            2
          
        
        ⁢
        
          
            ∑
            
              m
              =
              1
            
            M
          
          ⁢
          
            
              
                2
                ×
                
                  
                    u
                    
                      t
                      ′′
                    
                  
                  ⁡
                  
                    (
                    m
                    )
                  
                
                ×
                
                  
                    v
                    
                      t
                      ′′
                    
                  
                  ⁡
                  
                    (
                    m
                    )
                  
                
              
              +
              c
            
            
              
                
                  (
                  
                    
                      u
                      
                        t
                        ′′
                      
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                  )
                
                2
              
              +
              
                
                  (
                  
                    
                      v
                      
                        t
                        ′′
                      
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                  )
                
                2
              
              +
              c
            
          
        
      
    
  
  ,




wherein: 1≦m≦M; the ut′′(m) represents a value of an mth element in the ut′′; the vt′′(m) represents a value of an mth element in the vt′′; and the c is a small constant for guaranteeing a result stability.

Preferably, the Yw in the step {circle around (6)} is obtained through following steps of:
{circle around (6)}_1, representing a covariance matrix of the Y by C,




  
    C
    =
    
      
        1
        
          t
          2
        
      
      ⁢
      
        (
        
          Y
          ×
          
            Y
            T
          
        
        )
      
    
  
  ,




wherein: the C has a dimensionality of 192×192, and the YT is a transposition of the Y;

{circle around (6)}_2, processing the C with eigenvalue decomposition, and obtaining all maximum eigenvalues and corresponding eigenvectors, wherein the eigenvectors have a dimensionality of 192×1;
{circle around (6)}_3, choosing M maximum eigenvalues and corresponding M eigenvectors;
{circle around (6)}_4, according to the chosen M maximum eigenvalues and the corresponding M eigenvectors, calculating the whitened matrix W, W=Ψ−1/2×ET, wherein: the Ψ has a dimensionality of M×M, Ψ=diag(ψ1, . . . , ψM), Ψ−1/2=diag(1/√{square root over (ψ1)}, . . . , 1/√{square root over (ψM)}); the E has a dimensionality of 192×M, E=[e1, . . . , eM]; the diag( ) is a main-diagonal matrix representation; the ψ1, . . . , ψM correspondingly represent a 1st to a Mth chosen maximum eigenvalue; and, the e1, . . . , eM correspondingly represent a 1st to a Mth chosen eigenvector; and
{circle around (6)}_5, according to the W, processing the Y with the whitening, and obtaining the Yw after the dimensionality reduction and the whitening, Yw=W×Y.
Preferably, in the step {circle around (4)}, λ1=0.7.
Preferably, in the {circle around (5)}, λ2=0.6.
Preferably, in the step {circle around (9)}, c=0.04.
Compared with prior arts, the present invention has following advantages.
1) The method provided by the present invention considers a relationship between the saliency and the objective assessment of the image quality. Through the visual saliency detection algorithm, the saliency maps of the reference image and the distorted image are obtained for further obtaining the maximum fusion saliency map. Moreover, based on maximum saliencies of the image blocks in the maximum fusion saliency map, the saliency difference between each reference image block and the corresponding distorted image block is measured through the absolute difference, and thus the reference visual important image blocks and the distorted visual important image blocks are screened and extracted. Furthermore, through the manifold eigenvectors of the reference visual important image blocks and the distorted visual important image blocks, the objective quality assessment value of the distorted image is calculated. The method provided by the present invention has an obviously increased assessment effect and a higher correlation between an objective assessment result and a subjective perception.
2) The method provided by the present invention finds an inner geometric structure of data from image data through the manifold learning, and obtains the characteristic basis matrix after training. Through the characteristic basis matrix, the method further processes the reference visual important image blocks and the distorted visual important image blocks with the dimensionality reduction and obtains the manifold eigenvector. The manifold eigenvector after the dimensionality reduction still maintains geometric characteristics of high dimensional image data, and thus redundant information is decreased. Thus, when calculating the objective quality assessment value of the distorted image, the method provided by the present invention is easier and more accurate.
3) For conventional objective assessment methods for the image quality based on sparse representation, off-line learning of a over-complete dictionary needs a large number of effective training samples, and has a limitation on processing images having real-time requirements. For the above-described problem, the present invention online learns and trains the extracted reference visual important image blocks through the orthogonal locality preserving projection algorithm for obtaining the characteristic basis matrix in real-time. Thus, the present invention has a higher robustness and a more stable assessment effect.
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 implementation block diagram of an objective assessment method for a color image quality based on online manifold learning according to a preferred embodiment of the present invention.
FIG. 2a is a (scattered point)-(fitting curve) graph of the objective assessment method on a LIVE image database according to the preferred embodiment of the present invention.
FIG. 2b is a (scattered point)-(fitting curve) graph of the objective assessment method on a CSIQ image database according to the preferred embodiment of the present invention.
FIG. 2c is a (scattered point)-(fitting curve) graph of the objective assessment method on a TID2008 image database according to the preferred embodiment of the present invention.



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT
The present invention is further described with accompanying drawings and the preferred embodiment.
According to a preferred embodiment of the present invention, the present invention provides an objective assessment method for a color image quality based on online manifold learning, wherein an implementation block diagram thereof is showed in FIG. 1, and the method comprises steps of:
{circle around (1)} representing an undistorted reference image having a width W and a height H by IR; and representing a distorted image to be assessed, which is corresponding to the IR, by ID;
{circle around (2)} through a conventional visual saliency detection algorithm, named Saliency Detection based on Simple Priors (SDSP) herein, respectively obtaining saliency maps of the IR and the ID, correspondingly denoted as MR and MD; then, according to the MR and the MD, calculating a maximum fusion saliency map, denoted as MF; and denoting a pixel value of a pixel having coordinates of (x, y) in the MF as MF(x, y), MF(x, y)=max (MR(x, y), MD(x, y)), wherein: 1≦x≦W, 1≦y≦H; the max( ) is a function to find a maximum; the MR(x, y) represents a pixel value of a pixel having coordinates of (x, y) in the MR; and the MD(x, y) represents a pixel value of a pixel having coordinates of (x, y) in the MD;
{circle around (3)} respectively dividing the IR, the ID, the MR, the MD, and the MF into




  
    W
    ×
    H
  
  
    8
    ×
    8
  




image blocks, which are not overlapping mutually and have a size of 8×8; if a size of the IR, the ID, the MR, the MD, and the MF is indivisible by 8×8, redundant pixels are not processed;

vectorizing color values of R, G, and B channels of all pixels in each image block of the IR and the ID; denoting a color vector obtained through vectorizing the color values of the R, the G, and the B channels of all the pixels in a jth image block of the IR as XjR; denoting a color vector obtained through vectorizing the color values of the R, the G, and the B channels of all the pixels in a jth image block of the ID as XjD; wherein: the j has an initial value of 1,




  
    1
    ≤
    j
    ≤
    
      
        W
        ×
        H
      
      
        8
        ×
        8
      
    
  
  ;




both of the XjR and the XjD have a dimensionality of 192×1; values of a 1st element to a 64th element in the XjR respectively correspond to the color value of the R channel of each pixel in the jth image block of the IR in a line-by-line scanning manner, namely the value of the 1st element in the XjR is the color value of the R channel of a pixel in a 1st row and a 1st column of the jth image block of the IR, the value of a 2nd element in the XjR is the color value of the R channel of a pixel in the 1st row and a 2nd column of the jth image block of the IR, and so on; values of a 65th element to a 128th element in the XjR respectively correspond to the color value of the G channel of each pixel in the jth image block of the IR in the line-by-line scanning manner, namely the value of the 65th element in the XjR is the color value of the G channel of a pixel in a 1st row and a 1st column of the jth image block of the IR, the value of a 66th element in the XjR is the color value of the G channel of a pixel in the 1st row and a 2nd column of the jth image block of the IR, and so on; values of a 129th element to a 192nd element in the XjR respectively correspond to the color value of the B channel of each pixel in the jth image block of the IR in the line-by-line scanning manner, namely the value of the 129th element in the XjR is the color value of the B channel of a pixel in a 1st row and a 1st column of the jth image block of the IR, the value of a 130th element in the XjR is the color value of the B channel of a pixel in the 1st row and a 2nd column of the jth image block of the IR, and so on; values of a 1st element to a 64th element in the XjD respectively correspond to the color value of the R channel of each pixel in the jth image block of the ID in the line-by-line scanning manner, namely the value of the 1st element in the XjD is the color value of the R channel of a pixel in a 1st row and a 1st column of the jth image block of the ID, the value of a 2nd element in the XjD is the color value of the R channel of a pixel in the 1st row and a 2nd column of the jth image block of the ID, and so on; values of a 65th element to a 128th element in the XjD respectively correspond to the color value of the G channel of each pixel in the jth image block of the ID in the line-by-line scanning manner, namely the value of the 65th element in the XjD is the color value of the G channel of a pixel in a 1st row and a 1st column of the jth image block of the ID, the value of a 66th element in the XjD is the color value of the G channel of a pixel in the 1st row and a 2nd column of the jth image block of the ID, and so on; and values of a 129th element to a 192nd element in the XjD respectively correspond to the color value of the B channel of each pixel in the jth image block of the ID in the line-by-line scanning manner, namely the value of the 129th element in the XjD is the color value of the B channel of a pixel in a 1st row and a 1st column of the jth image block of the ID, the value of a 130th element in the XjD is the color value of the B channel of a pixel in the 1st row and a 2nd column of the jth image block of the ID, and so on; and

vectorizing pixel values of all pixels in each image block of the MR, the MD, and the MF; denoting a pixel value vector obtained through vectorizing the pixel values of all the pixels in a jth image block of the MR as SjR; denoting a pixel value vector obtained through vectorizing the pixel values of all the pixels in a jth image block of the MD as SjD; and denoting a pixel value vector obtained through vectorizing the pixel values of all the pixels in a jth image block of the MF as SjF; wherein: the SjR, the SjD, and the SjF all have a dimensionality of 64×1; values of a 1st element to a 64th element in the SjR respectively correspond to the pixel value of each pixel in the jth image block of the MR in the line-by-line scanning manner, namely the value of the 1st element in the SjR is the pixel value of a pixel in a 1st row and a 1st column of the jth image block of the MR, the value of a 2nd element in the SjR is the pixel value of a pixel in the 1st row and a 2nd column of the jth image block of the MR, and so on; values of a 1st element to a 64th element in the SjD respectively correspond to the pixel value of each pixel in the jth image block of the MD in the line-by-line scanning manner, namely the value of the 1st element in the SjD is the pixel value of a pixel in a 1st row and a 1st column of the jth image block of the MD, the value of a 2nd element in the SD is the pixel value of a pixel in the 1st row and a 2nd column of the jth image block of the MD, and so on; and values of a 1st element to a 64th element in the SjF respectively correspond to the pixel value of each pixel in the jth image block of the MF in the line-by-line scanning manner, namely the value of the 1st element in the SjF is the pixel value of a pixel in a 1st row and a 1st column of the jth image block of the MF, the value of a 2nd element in the SjF is the pixel value of a pixel in the 1st row and a 2nd column of the jth image block of the MF, and so on;
{circle around (4)} calculating a saliency of each image block in the MF; and denoting the saliency of the jth image block in the MF as dj,




  
    
      d
      j
    
    =
    
      
        ∑
        
          i
          =
          1
        
        64
      
      ⁢
      
        
          S
          j
          F
        
        ⁡
        
          (
          i
          )
        
      
    
  
  ,




wherein: 1≦i≦64; and the SjF(i) represents a value of an ith element in the SjF, namely a pixel value of an ith pixel in the jth image block of the MF;

orderly arranging the saliencies of all the image blocks in the MF, from the biggest to the smallest; and, after arranging, determining sequence numbers of the image blocks corresponding to former t1 saliencies (namely maximum t1 saliencies), wherein:




  
    
      t
      1
    
    =
    
      
        λ
        1
      
      ×
      
        
          W
          ×
          H
        
        
          8
          ×
          8
        
      
    
  
  ;




the λ1 represents an image block selection proportionality coefficient, λ1ε(0,1]; and it is embodied that λ1=0.7 herein; and

finding the image blocks in the IR, which are corresponding to the determined t1 sequence numbers, and defining as reference image blocks; finding the image blocks in the ID; which are corresponding to the determined t1 sequence numbers, and defining as distorted image blocks; finding the image blocks in the MR, which are corresponding to the determined t1 sequence numbers, and defining as reference saliency image blocks; finding the image blocks in the MD, which are corresponding to the determined t1 sequence numbers, and defining as distorted saliency image blocks;
{circle around (5)} measuring a saliency difference between each reference image block in the IR and a corresponding distorted image block in the ID through an absolute difference; and denoting a saliency difference between a t′th reference image block in the IR and a t′th distorted image block in the ID as et′,




  
    
      e
      
        t
        ′
      
    
    =
    
      
        1
        64
      
      ⁢
      
        
          ∑
          
            i
            =
            1
          
          64
        
        ⁢
        
          
          
            
              
                
                  S
                  ~
                
                
                  t
                  ′
                
                R
              
              ⁡
              
                (
                i
                )
              
            
            -
            
              
                
                  S
                  ~
                
                
                  t
                  ′
                
                D
              
              ⁡
              
                (
                i
                )
              
            
          
          
        
      
    
  
  ,




wherein: the t′ has an initial value of 1, 1≦t′≦t1; the symbol “∥” is an absolute value symbol; the {tilde over (S)}t′R(i) represents a value of an ith element in a pixel value vector {tilde over (S)}t′R corresponding to a t′th reference saliency image block in the MR, namely a pixel value of an ith pixel in the t′th reference saliency image block of the MR; and the {tilde over (S)}t′D(i) represents a value of an ith element in a pixel value vector {tilde over (S)}t′D corresponding to a t′th distorted saliency image block in the MD, namely a pixel value of an ith pixel in the t′th distorted saliency image block of the MD; and

orderly arranging the measured t1 saliency differences, from the biggest to the smallest; after arranging, determining the reference image blocks and the distorted image blocks corresponding to former t2 saliency differences (namely maximum t2 saliency differences); defining the determined t2 reference image blocks as reference visual important image blocks, and adopting a matrix formed by color vectors corresponding to all the reference visual important image blocks as a reference visual important image block matrix, denoted as YR; defining the determined t2 distorted image blocks as distorted visual important image blocks, and adopting a matrix formed by color vectors corresponding to all the distorted visual important image blocks as a distorted visual important image block matrix, denoted as YD; wherein: t2=λ2×t1; the λ2 represents a selection proportionality coefficient of the reference image blocks and the distorted image blocks, λ2ε(0,1]; it is embodied that λ2=0.6 herein; the YR and the YD have a dimensionality of 192×t2; a t′′th column vector in the YR is a color vector corresponding to a determined t′′th reference visual important image block; a t′′th column vector in the YD is a color vector corresponding to a determined t′′th distorted visual important image block; and the t′′ has an initial value of 1, 1≦t′′≦t2;
{circle around (6)} centralizing the YR through subtracting a mean value of values of all elements in each column vector from a value of each element in the same column vector of the YR; and denoting an obtained centralized matrix as Y, wherein the Y has a dimensionality of 192×t2; and
processing the Y obtained after centralizing the YR with dimensionality reduction and whitening through a conventional principal component analysis (PCA); and, denoting an obtained matrix after the dimensionality reduction and the whitening as Yw, Yw=W×Y; wherein: the Yw has a dimensionality of M×t2; the W represents a whitened matrix and has a dimensionality of M×192, 1<M<<192; the symbol “<<” is a much-less-than symbol; and, it is embodied that the PCA is realized through processing a covariance matrix of sample data with eigenvalue decomposition, namely the Yw in the step {circle around (6)} is obtained through following steps of:
{circle around (6)}_1, representing a covariance matrix of the Y by C,




  
    C
    =
    
      
        1
        
          t
          2
        
      
      ⁢
      
        (
        
          Y
          ×
          
            Y
            T
          
        
        )
      
    
  
  ,




wherein: the C has a dimensionality of 192×192, and the YT is a transposition of the Y;

{circle around (6)}_2, processing the C with the eigenvalue decomposition, and obtaining all maximum eigenvalues and corresponding eigenvectors, wherein the eigenvectors have a dimensionality of 192×1;
{circle around (6)}_3, choosing M maximum eigenvalues and corresponding M eigenvectors, so as to realize the dimensionality reduction of the Y, wherein: it is embodied that M=8 herein, namely merely former eight principle components are chosen for training and thus the dimensionality is decreased from 192 to M=8;
{circle around (6)}_4, according to the chosen M maximum eigenvalues and the corresponding M eigenvectors, calculating the whitened matrix W, W=Ψ−1/2×ET, wherein: the Ψ has a dimensionality of M×M, Ψ=diag(ψ1, . . . , ψM), Ψ−1/2=diag(1/√{square root over (ψ1)}, . . . , 1/√{square root over (ψM)}); the E has a dimensionality of 192×M, E=[e1, . . . , eM]; the diag( ) is a main-diagonal matrix representation; the ψ1, . . . , ψM correspondingly represent a 1st to a Mth chosen maximum eigenvalue; and, the e1, . . . , eM correspondingly represent a 1st to a Mth chosen eigenvector; and
{circle around (6)}_5, according to the W, processing the Y with the whitening, and obtaining the Yw after the dimensionality reduction and the whitening, Yw=W×Y;
{circle around (7)} online training the Yw through an orthogonal locality preserving projection (OLPP) algorithm; and obtaining a characteristic basis matrix of the Yw, denoted as D, wherein the D has a dimensionality of M×192;
{circle around (8)} according to the YR and the D, calculating a manifold eigenvector of each reference visual important image block; denoting the manifold eigenvector of a t′′th reference visual important image block as ut′′, ut′′=D×yt′′R, wherein: the ut′′ has a dimensionality of M×1, and the yt′′R is the t′′th column vector in the YR; according to the YD and the D, calculating a manifold eigenvector of each distorted visual important image block; and denoting the manifold eigenvector of a t′′th distorted visual important image block as vt′′, vt′′=D×yt′′D, wherein: the vt′′ has a dimensionality of M×1, and the yt′′D is the t′′th column vector in the YD; and
{circle around (9)} according to the manifold eigenvectors of all the reference visual important image blocks and the manifold eigenvectors of all the distorted visual important image blocks, calculating an objective quality assessment value of the ID, denoted as Score,




  
    Score
    =
    
      
        1
        
          
            t
            2
          
          ×
          M
        
      
      ⁢
      
        
          ∑
          
            
              t
              ′′
            
            =
            1
          
          
            t
            2
          
        
        ⁢
        
          
            ∑
            
              m
              =
              1
            
            M
          
          ⁢
          
            
              
                2
                ×
                
                  
                    u
                    
                      t
                      ′′
                    
                  
                  ⁡
                  
                    (
                    m
                    )
                  
                
                ×
                
                  
                    v
                    
                      t
                      ′′
                    
                  
                  ⁡
                  
                    (
                    m
                    )
                  
                
              
              +
              c
            
            
              
                
                  (
                  
                    
                      u
                      
                        t
                        ′′
                      
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                  )
                
                2
              
              +
              
                
                  (
                  
                    
                      v
                      
                        t
                        ′′
                      
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                  )
                
                2
              
              +
              c
            
          
        
      
    
  
  ,




wherein: 1≦m≦M; the ut′′(m) represents a value of an mth element in the ut′′; the vt′′(m) represents a value of an mth element in the vt′′; the c is a small constant for guaranteeing a result stability, and it is embodied that c=0.04 herein.

In order to further illustrate effectiveness and feasibility of the method provided by the present invention, the method is tested.
According to the preferred embodiment of the present invention, three open authoritative image databases are chosen to be tested, respectively a LIVE image database, a CSIQ image database, a TID2008 image database. In Table 1, various indexes of the three image databases are described in detail. The various indexes comprise a reference image number, a distorted image number and a distorted type number. All of the three databases provide a mean subjective assessment difference of each distorted image.







TABLE 1







various indexes of authoritative image databases











Reference image
Distorted image
Distorted type


Image database
number
number
number













LIVE
29
779
5


CSIQ
30
866
6


TID2008
25
1700
17









Then, a correlation between the objective quality assessment value obtained by the method of the present invention and the mean subjective assessment difference of each distorted image is analyzed. Herein, three common objective parameters for assessing an image quality assessment method serve as assessment indexes. The three objective parameters are respectively a Pearson linear correlation coefficient (PLCC) which reflects a prediction accuracy, a Spearman rank order correlation coefficient (SROCC) which reflects a prediction monotonicity, and a root mean squared error (RMSE) which reflects a prediction consistency. A value range of the PLCC and the SROCC is [0, 1]. The nearer a value of the PLCC and the SROCC approximates to 1, the better an image quality objective assessment method is; otherwise, the image quality objective assessment method is worse. The smaller RMSE, the higher predication accuracy and the better performance of the image quality objective assessment method; otherwise, the predication accuracy is lower and the performance is worse.
For all distorted images in the above LIVE image database, CSIQ image database and TID2008 image database, the objective quality assessment value of each distorted image is calculated in a same manner through the steps {circle around (1)}-{circle around (9)} of the method provided by the present invention. The obtained correlation between the objective quality assessment value and the mean subjective assessment difference of the distorted image is analyzed. Firstly, the objective quality assessment value is obtained; then, the objective quality assessment value is processed with five-parameter Logistic function non-linear fitting; and finally, a performance index value between an objective assessment result and the mean subjective assessment difference is obtained. In order to verify the effectiveness of the present invention, on the three image databases listed in Table 1, the method provided by the present invention and other six conventional full-reference image quality objective assessment methods having a relatively advanced performance are comparatively analyzed. The PLCC, the SROCC, and the RMSE for representing an assessment performance of the methods are listed in Table 2. In Table 2, the six methods for comparing are respectively a classical peak signal-to-noise ratio (PSNR) method, an assessment method based on a structural similarity (SSIM) proposed by Z. Wang, a method based on a degradation model named Information Fidelity Criterion (IFC) proposed by N. Damera Venkata, a method based on a visual information fidelity (VIF) proposed by H. R. Sheikh, a method based on a wavelet visual signal-to-noise ratio (VSNR) proposed by D. M. Chandler, and an image quality assessment method based on a sparse representation denoted as Sparse Representation-based Quality (SPARQ) proposed by T. Guha. According to data listed in Table 2, on the LIVE image database, the method provided by the present invention has the second best performance after the VIF method; and, on the CSIQ image database and the TID image database, the method provided by the present invention has the best performance. Thus, for all of the above-described three image databases, the objective quality assessment value of the distorted image obtained by the method of the present invention has a good correlation with the mean subjective assessment difference. Moreover, the values of the PLCC and the SROCC of the LIVE image database and the CSIQ image database are all above 0.94; the values of the PLCC and the SROCC of the TID2008 image database having more complex distorted types also reach 0.82; and, after weighted averaging, compared with all of the six conventional methods, the performance of the method provided by the present invention has different degrees of improvement. Thus, the objective assessment result of the method provided by the present invention is relatively consistent with the subjective perception of the human eyes, and has a stable assessment effect, which fully illustrates the effectiveness of the method provided by the present invention.







TABLE 2







performance comparison between the method provided by the present invention


and the conventional image quality objective assessment methods























Method of


Image







present


database

PSNR
SSIM
IFC
VIF
VSNR
SPARQ
invention


















LIVE
SROCC
0.8756
0.9479
0.9259
0.9636
0.9274
0.9310
0.9523



PLCC
0.8723
0.9449
0.9268
0.9604
0.9231
0.9280
0.9506



RMSE
13.3600
8.9455
10.2641
7.6137
10.5060
10.1850
8.4433


CSIQ
SROCC
0.8057
0.8756
0.7671
0.9195
0.8106
0.9460
0.9465



PLCC
0.8000
0.8613
0.8384
0.9277
0.8002
0.9390
0.9433



RMSE
0.1575
0.1344
0.1431
0.0980
0.1575
0.0900
0.0871


TID2008
SROCC
0.5531
0.7749
0.5675
0.7491
0.7046
0.7920
0.8356



PLCC
0.5734
0.7732
0.7340
0.8084
0.6820
0.8200
0.8228



RMSE
1.0994
0.8511
0.9113
0.7899
0.9815
0.7680
0.5975


Average
SROCC
0.6936
0.8413
0.7026
0.8432
0.7839
0.8642
0.9115



PLCC
0.7017
0.8360
0.8059
0.8747
0.7687
0.8760
0.9056



RMSE
4.8723
3.3103
3.7728
2.8339
3.8817
3.6810
3.0426









FIG. 2a shows a (scattered point)-(fitting curve) graph of the method provided by the present invention on the LIVE image database. FIG. 2b shows a (scattered point)-(fitting curve) graph of the method provided by the present invention on the CSIQ image database. FIG. 2c shows a (scattered point)-(fitting curve) graph of the method provided by the present invention on the TID2008 image database. From FIG. 2a, FIG. 2b, and FIG. 2c, it can be clearly seen that the scattered points are uniformly distributed near the fitting curve and show a good monotonicity and consistency.
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. An objective assessment method for a color image quality based on online manifold learning, comprising steps of:
{circle around (1)} representing an undistorted reference image having a width W and a height H by IR; and representing a distorted image to be assessed, which is corresponding to the IR, by ID;
{circle around (2)} through a visual saliency detection algorithm, respectively obtaining saliency maps of the IR and the ID, correspondingly denoted as MR and MD; then, according to the MR and the MD, calculating a maximum fusion saliency map, denoted as MF; and denoting a pixel value of a pixel having coordinates of (x, y) in the MF as MF(x, y), MF(x, y)=max (MR(x, y), MD(x, y)), wherein: 1≦x≦W, 1≦y≦H; the max( ) is a function to find a maximum; the MR(x, y) represents a pixel value of a pixel having coordinates of (x, y) in the MR; and the MD(x, y) represents a pixel value of a pixel having coordinates of (x, y) in the MD;
{circle around (3)} respectively dividing the IR, the ID, the MR, the MD, and the MF into





  
    W
    ×
    H
  
  
    8
    ×
    8
  




image blocks, which are not overlapping mutually and have a size of 8×8;
vectorizing color values of R, G, and B channels of all pixels in each image block of the IR and the ID; denoting a color vector obtained through vectorizing the color values of the R, the G, and the B channels of all the pixels in a jth image block of the IR as XjR; denoting a color vector obtained through vectorizing the color values of the R, the G, and the B channels of all the pixels in a jth image block of the ID as XjD; wherein: the j has an initial value of 1,





  
    1
    ≤
    j
    ≤
    
      
        W
        ×
        H
      
      
        8
        ×
        8
      
    
  
  ;




both of the XjR and the XjD have a dimensionality of 192×1; values of a 1st element to a 64th element in the XjR respectively correspond to the color value of the R channel of each pixel in the jth image block of the IR in a line-by-line scanning manner; values of a 65th element to a 128th element in the XjR respectively correspond to the color value of the G channel of each pixel in the jth image block of the IR in the line-by-line scanning manner; values of a 129th element to a 192nd element in the XjR respectively correspond to the color value of the B channel of each pixel in the jth image block of the IR in the line-by-line scanning manner; values of a 1st element to a 64th element in the XjD respectively correspond to the color value of the R channel of each pixel in the jth image block of the ID in the line-by-line scanning manner; values of a 65th element to a 128th element in the XjD respectively correspond to the color value of the G channel of each pixel in the jth image block of the ID in the line-by-line scanning manner; and values of a 129th element to a 192nd element in the XjD respectively correspond to the color value of the B channel of each pixel in the jth image block of the ID in the line-by-line scanning manner; and
vectorizing pixel values of all pixels in each image block of the MR, the MD, and the MF; denoting a pixel value vector obtained through vectorizing the pixel values of all the pixels in a jth image block of the MR as SjR; denoting a pixel value vector obtained through vectorizing the pixel values of all the pixels in a jth image block of the MD as SjD; and denoting a pixel value vector obtained through vectorizing the pixel values of all the pixels in a jth image block of the MF as SjF; wherein: the SjR, the SjD, and the SjF all have a dimensionality of 64×1; values of a 1st element to a 64th element in the SjR respectively correspond to the pixel value of each pixel in the jth image block of the MR in the line-by-line scanning manner; values of a 1st element to a 64th element in the SjD respectively correspond to the pixel value of each pixel in the jth image block of the MD in the line-by-line scanning manner; and values of a 1st element to a 64th element in the SjF respectively correspond to the pixel value of each pixel in the jth image block of the MF in the line-by-line scanning manner;
{circle around (4)} calculating a saliency of each image block in the MF; and denoting the saliency of the jth image block in the MF as dj,





  
    
      d
      j
    
    =
    
      
        ∑
        
          i
          =
          1
        
        64
      
      ⁢
      
        
          S
          j
          F
        
        ⁡
        
          (
          i
          )
        
      
    
  
  ,




wherein: 1≦i≦64 and the SjF(i) represents a value of an ith element in the SjF;
orderly arranging the saliencies of all the image blocks in the MF, from the biggest to the smallest; and, after arranging, determining sequence numbers of the image blocks corresponding to former t1 saliencies, wherein:





  
    
      t
      1
    
    =
    
      
        λ
        1
      
      ×
      
        
          W
          ×
          H
        
        
          8
          ×
          8
        
      
    
  
  ;




the λ1 represents an image block selection proportionality coefficient; and λ1ε(0,1]; and
finding the image blocks in the IR, which are corresponding to the determined t1 sequence numbers, and defining as reference image blocks; finding the image blocks in the ID, which are corresponding to the determined t1 sequence numbers, and defining as distorted image blocks; finding the image blocks in the MR, which are corresponding to the determined t1 sequence numbers, and defining as reference saliency image blocks; finding the image blocks in the MD, which are corresponding to the determined t1 sequence numbers, and defining as distorted saliency image blocks;
{circle around (5)} measuring a saliency difference between each reference image block in the IR and a corresponding distorted image block in the ID through an absolute difference; and denoting a saliency difference between a t′th reference image block in the IR and a t′th distorted image block in the ID as et′,





  
    
      e
      
        t
        ′
      
    
    =
    
      
        1
        64
      
      ⁢
      
        
          ∑
          
            i
            =
            1
          
          64
        
        ⁢
        
          
          
            
              
                
                  S
                  ~
                
                
                  t
                  ′
                
                R
              
              ⁡
              
                (
                i
                )
              
            
            -
            
              
                
                  S
                  ~
                
                
                  t
                  ′
                
                D
              
              ⁡
              
                (
                i
                )
              
            
          
          
        
      
    
  
  ,




wherein: the t′ has an initial value of 1, 1≦t′≦t1; the symbol “∥” is an absolute value symbol; the {tilde over (S)}t′R(i) represents a value of an ith element in a pixel value vector {tilde over (S)}t′R corresponding to a t′th reference saliency image block in the MR; and the {tilde over (S)}t′D(i) represents a value of an ith element in a pixel value vector {tilde over (S)}t′D corresponding to a t′th distorted saliency image block in the MD; and
orderly arranging the measured t1 saliency differences, from the biggest to the smallest; after arranging, determining the reference image blocks and the distorted image blocks corresponding to former t2 saliency differences; defining the determined t2 reference image blocks as reference visual important image blocks, and adopting a matrix formed by color vectors corresponding to all the reference visual important image blocks as a reference visual important image block matrix, denoted as YR; defining the determined t2 distorted image blocks as distorted visual important image blocks, and adopting a matrix formed by color vectors corresponding to all the distorted visual important image blocks as a distorted visual important image block matrix, denoted as YD; wherein: t2=λ2×t1, the λ2 represents a selection proportionality coefficient of the reference image blocks and the distorted image blocks, and λ2ε(0,1]; the YR and the YD have a dimensionality of 192×t2; a t′′th column vector in the YR is a color vector corresponding to a determined t′′th reference visual important image block; a t′′th column vector in the YD is a color vector corresponding to a determined t′′th distorted visual important image block; and the t′′ has an initial value of 1, 1≦t′′≦t2;
{circle around (6)} centralizing the YR through subtracting a mean value of values of all elements in each column vector from a value of each element in the same column vector of the YR; and denoting an obtained centralized matrix as Y, wherein the Y has a dimensionality of 192×t2; and
processing the Y with dimensionality reduction and whitening through a principal component analysis; and, denoting an obtained matrix after the dimensionality reduction and the whitening as Yw, Yw=W×Y; wherein: the Yw has a dimensionality of M×t2; the W represents a whitened matrix and has a dimensionality of M×192, 1<M<<192; and the symbol “<<” is a much-less-than symbol;
{circle around (7)} online training the Yw through an orthogonal locality preserving projection algorithm; and obtaining a characteristic basis matrix of the Yw, denoted as D, wherein the D has a dimensionality of M×192;
{circle around (8)} according to the YR and the D, calculating a manifold eigenvector of each reference visual important image block; denoting the manifold eigenvector of a t′′th reference visual important image block as ut′′, ut′′=D×yt′′R, wherein: the ut′′ has a dimensionality of M×1, and the yt′′R is the t′′th column vector in the YR; according to the YD and the D, calculating a manifold eigenvector of each distorted visual important image block; and denoting the manifold eigenvector of a t′′th distorted visual important image block as vt′′, vt′′=D×yt′′D, wherein: the vt′′ has a dimensionality of M×1, and the yt′′D is the t′′th column vector in the YD; and
{circle around (9)} according to the manifold eigenvectors of all the reference visual important image blocks and the manifold eigenvectors of all the distorted visual important image blocks, calculating an objective quality assessment value of the ID, denoted as Score,





  
    Score
    =
    
      
        1
        
          
            t
            2
          
          ×
          M
        
      
      ⁢
      
        
          ∑
          
            
              t
              ′′
            
            =
            1
          
          
            t
            2
          
        
        ⁢
        
          
            ∑
            
              m
              =
              1
            
            M
          
          ⁢
          
            
              
                2
                ×
                
                  
                    u
                    
                      t
                      ′′
                    
                  
                  ⁡
                  
                    (
                    m
                    )
                  
                
                ×
                
                  
                    v
                    
                      t
                      ′′
                    
                  
                  ⁡
                  
                    (
                    m
                    )
                  
                
              
              +
              c
            
            
              
                
                  (
                  
                    
                      u
                      
                        t
                        ′′
                      
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                  )
                
                2
              
              +
              
                
                  (
                  
                    
                      v
                      
                        t
                        ′′
                      
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                  )
                
                2
              
              +
              c
            
          
        
      
    
  
  ,




wherein: 1≦m≦M; the ut′′(m) represents a value of an mth element in the ut′′; the vt′′(m) represents a value of an mth element in the vt′′; and the c is a small constant for guaranteeing a result stability.

2. The objective assessment method for the color image quality based on the online manifold learning, as recited in claim 1, wherein the Yw in the step {circle around (6)} is obtained through following steps of:
{circle around (6)}_1, representing a covariance matrix of the Y by C,





  
    C
    =
    
      
        1
        
          t
          2
        
      
      ⁢
      
        (
        
          Y
          ×
          
            Y
            T
          
        
        )
      
    
  
  ,




wherein: the C has a dimensionality of 192×192, and the YT is a transposition of the Y;
{circle around (6)}_2, processing the C with eigenvalue decomposition, and obtaining all maximum eigenvalues and corresponding eigenvectors, wherein the eigenvectors have a dimensionality of 192×1;
{circle around (6)}_3, choosing M maximum eigenvalues and corresponding M eigenvectors;
{circle around (6)}_4, according to the chosen M maximum eigenvalues and the corresponding M eigenvectors, calculating the whitened matrix W, W=Ψ−1/2×ET, wherein: the Ψ has a dimensionality of M×M, Ψ=diag(ψ1, . . . , ψM), Ψ−1/2=diag(1/√{square root over (ψ1)}, . . . , 1/√{square root over (ψM)}); the E has a dimensionality of 192×M, E=[e1, . . . , eM]; the diag( ) is a main-diagonal matrix representation; the ψ1, . . . , ψM correspondingly represent a 1st to a Mth chosen maximum eigenvalue; and, the e1, . . . , eM correspondingly represent a 1st to a Mth chosen eigenvector; and
{circle around (6)}_5, according to the W, processing the Y with the whitening, and obtaining the Yw after the dimensionality reduction and the whitening, Yw=W×Y.

3. The objective assessment method for the color image quality based on the online manifold learning, as recited in claim 1, wherein: in the step {circle around (4)}, λ1=0.7.

4. The objective assessment method for the color image quality based on the online manifold learning, as recited in claim 1, wherein: in the {circle around (5)}, λ2=0.6.

5. The objective assessment method for the color image quality based on the online manifold learning, as recited in claim 1, wherein: in the step {circle around (9)}, c=0.04.