---
abstract: 'Light field image quality assessment (LF-IQA) plays a significant role due to its guidance to Light Field (LF) contents acquisition, processing and application. The LF can be represented as 4-D signal, and its quality depends on both angular consistency and spatial quality. However, few existing LF-IQA methods concentrate on effects caused by angular inconsistency. Especially, no-reference methods lack effective utilization of 2-D angular information. In this paper, we focus on measuring the 2-D angular consistency for LF-IQA. The Micro-Lens Image (MLI) refers to the angular domain of the LF image, which can simultaneously record the angular information in both horizontal and vertical directions. Since the MLI contains 2-D angular information, we propose a No-Reference Light Field image Quality assessment model based on MLI (LF-QMLI). Specifically, we first utilize Global Entropy Distribution (GED) and Uniform Local Binary Pattern descriptor (ULBP) to extract features from the MLI, and then pool them together to measure angular consistency. In addition, the information entropy of Sub-Aperture Image (SAI) is adopted to measure spatial quality. Extensive experimental results show that LF-QMLI achieves the state-of-the-art performance.'
author:
- 
bibliography:
- 'bibfile.bib'
title: |
    No-Reference Light Field Image Quality Assessment Based on Micro-Lens Image\
    [^1] 
---

Light field, Image quality assessment, Objective model, Micro-lens image, Angular consistency

**Introduction** {#sec:intro}
================

As a representative of the attractive technique for immersive multimedia data, Light Field (LF) images have attracted widespread attention [@wu2017light]. Unlike traditional 2D images, LF images can record radiance information in both spatial and angular dimensions [@levoy2009recording], leading to better immersive experience. In order to provide satisfactory viewing quality of experience (QoE), light field image quality assessment (LF-IQA) plays a crucial role in LF contents acquisition, processing and application.

The LF image is a 4-D signal containing spatial and angular information. In Fig. \[fig:intro\], we show the LF image with different formats. Fig. \[fig:intro\](a) shows the LF image captured by a lenslet LF camera called Lytro Illum [@ng2005light]. The parameters $u$ and $v$ refer to angular dimensions while $s$ and $t$ represent spatial dimensions. We can obtain the Sub-Aperture Image (SAI) by fixing $u$ and $v$ [@van2018light], while the Micro-Lens Image (MLI) is through fixing $s$ and $t$ [@cho2013modeling], as shown in Fig. \[fig:intro\](b-c). The bottom and right of Fig. \[fig:intro\](b) are Epipolar-Plane Images (EPIs) [@wu2017light], which are produced by fixing ($u$,$s$) and ($v$,$t$). The SAI only contains spatial information of the LF image, while the EPI include both spatial and angular dimensions. However, the EPI only contains horizontal or vertical angular direction. Unlike the SAI and EPI, the MLI includes 2-D angular information. On account of the 4-D characteristic, the perceptual quality of LF image mainly depends on spatio-angular resolution, angular consistency and spatial quality [@wu2017light]. Concretely, spatio-angular resolution refers to the LF image resolution (i.e. the values of $u$,$v$,$s$ and $t$). Angular consistency measures the visual coherence of LF images while spatial quality indicates the SAI quality. Since spatio-angular resolution is an inherent factor of the LF image, we consider the effects of angular consistency and spatial quality in this paper.

![LF image with different formats. (a) Lenset image captured by Lytro Illum; (b) SAI array with EPIs in the bottom and right; (c) 9 MLIs in the red bounding box of (a) at high magnification.[]{data-label="fig:intro"}](intro.png){width="\linewidth"}

Although the subjective evaluation of LF-IQA  [@viola2016objective; @kiran2017towards; @shi2018perceptual] is precise and reliable, it is resource and time-consuming. Therefore, an effective objective LF-IQA model is urgently required. In general, image quality assessment (IQA) methods can be classified into three categories: full-reference (FR), reduced-reference (RR) and no-reference (NR) [@zhou2016binocular]. Among FR methods, intact information of original images is needed. Structure similarity between original and distorted images is measured in structural similarity index (SSIM) [@SSIM], with several variants, e.g. MS-SSIM [@MSSIM] and FSIM [@FSIM]. MP-PSNR Full [@MP_F] and MP-PSNR Reduc [@MP_R] based on Morphological pyramid decomposition are proposed to evaluate the multi-view image quality. RR methods only require part of information from original images. NR methods only utilize distorted images, which can be applied to most applications where reference images are hardly available, e.g. Mittal $et \ al.$ [@BRI] uses scene statistics in spatial domain and binocular fusion and rivalry are concerned in BSVQE [@BSVQE].

There exist only a few objective LF-IQA models. Fang $et \ al.$ [@fang2018light] proposes a FR LF-IQA method to compute the gradient magnitude similarity between original and distorted EPIs. Paudyal $et \ al.$ [@LFIQM] predicts the LF image quality with the structure similarity between original depth map and distorted depth map. However, neither of the methods consider the spatial quality degradation on the SAI. In addition, the EPIs only contain horizontal or vertical angular dimension, leading to insufficient measurement of angular consistency for the LF image applications. Therefore, a LF-IQA method that considers spatial quality and 2-D angular consistency is necessary for practical application.

In this paper, we propose a novel NR Light Field image Quality assessment model based on Micro-Lens Image (LF-QMLI) to evaluate both of angular consistency and spatial quality for LF images. As shown in Fig. \[fig:intro\](c), each pixel in the MLI comes from the same point in spatial domain, but is captured in various directions. Hence, there exists quite strong dependence between MLI pixels for 2-D angular consistency. To our best knowledge, we are the first to utilize MLI to evaluate angular consistency of LF images.

In this work, we first obtain the MLI by fixing $s$ and $t$, while the SAI is generated through fixing $u$ and $v$. Second, the Global Entropy Distribution (GED) and Uniform Local Binary Pattern descriptor (ULBP) are proposed to measure the angular consistency on each MLI, and then we utilize content pooling. Third, the information entropy of SAI is utilized to evaluate spatial quality. Finally, we train a regression model that predicts the perceptual quality of distorted LF images. Our experimental results demonstrate that proposed model LF-QMLI achieves the state-of-the-art performance.

The rest of the paper is organized as follows: Section \[method\] introduces our LF-QMLI method. Experimental results are shown in Section \[experiment\] and we conclude the paper in Section \[conclusion\].

![Flow diagram of proposed LF-QMLI model. (a) LF image (lenslet format); (b) MLI Array; (c) SAI Array.[]{data-label="fig:flowdiagram"}](flowdiagram.png){width="0.9\linewidth"}

**Proposed Method** {#method}
===================

The flow diagram of our proposed LF-QMLI model is shown in Fig. \[fig:flowdiagram\]. We first convert LF images from lenslet format into MLI and SAI arrays. Since the entropy can measure the angular dependence between adjacent pixels [@liu2014no], the GED is adopted on MLI to measure angular consistency. As the textural variation shown in Fig. \[fig:MLIS\], the ULBP is selected to measure the textural features for original and distorted MLIs. In addition, the information entropy of SAI is utilized to measure spatial quality [@hu2008constructing]. After content pooling, a regression model is used to predict the perceptual quality of LF images.

**Angular Consistency Based on MLI**
------------------------------------

The distortion caused by angular inconsistency affects LF image quality. A LF camera captures the same object in spatial domain with various angles of view, and the MLI is composed of light rays from both horizontal and vertical directions, leading to a 2-D angular domain.

### **Global Entropy Distribution (GED)** {#sec:GED}

In previous works, information entropy is proved as an efficient method to measure spatial image quality [@sponring1996entropy]. However, without considering angular dimension, it cannot work well in LF images.

The entropy of undistorted images possesses certain statistical properties, owing to the dependence between adjacent pixels [@liu2014no]. As shown in Fig. \[subfig:Ori\], since the variation between pixels is piecewise linear [@heber2013variational], the MLI without angular distortion is regular and gradually varied. With increased distortion, the dependence between adjacent pixels is destroyed, leading to the change of global entropy. In Fig. \[fig:MLIS\](b)-(e), the distortion caused by angular inconsistency primarily affects the MLI entropy.

Our proposed GED includes global image entropy distribution and global frequency entropy distribution of MLI.

The image entropy is $$E_{I}=-\sum_{x}P_{x}log_{2}P_{x},$$ where $x$ is the pixel value within a MLI, ranging from 0 to 255, with empirical probability density $P_{x}$.

The frequency entropy is $$E_{F}=-\sum_{i}\sum_{j}P_{i,j}log_{2}P_{i,j},$$ where $P_{i,j}$ is the value of probability map locating at $(i,j)$ in the DCT domain of MLI.

Finally, the MLI global entropy $E_{MLI}$ includes image entropy and frequency entropy. $$E_{MLI}=\{E_{I},E_{F}\}.$$

  --------------------------------------------------------------------------- -- -- -- -- --
   **E$_{MLI}$ & **Ori. Img & **LN lv.2 & **LN lv.5 & **NN lv.2 & **NN lv.5\              
              **IE & 5.5772 & 5.4971 & 5.3812 & 5.1807 & 4.3343\                          
              **FE & 2.0344 & 2.0543 & 3.0031 & 2.2356 & 3.3460\                          
                               ****************                                           
  --------------------------------------------------------------------------- -- -- -- -- --

  : E$_{MLI}$ on various angular distortion[]{data-label="tab:EMLI"}

We conducted validation experiments on various levels of distortion caused by angular inconsistency. We considered 2 kinds of angular distortion {linear interpolation (LN), nearest neighbor interpolation (NN)} and 2 levels of distortion {Level2, Level5} [@shi2018perceptual]. The higher level represents higher distortion. The image entropy (IE) and frequency entropy (FE) of 5 MLIs (Fig. \[fig:MLIS\]) were computed to demonstrate the validity of $E_{MLI}$. It is shown in Table \[tab:EMLI\] that the angular distortion can affect $E_{MLI}$ in a conspicuous and predictable way. The loss of image details is caused by the distortion, leading to image entropy reduction. Generally, the angular distortion of higher level obtains smaller IE and greater FE. On the same distortion level, NN destroys angular consistency more acutely than LN, which is verified in subjective experiments [@shi2018perceptual].

However, there exist a large number of MLIs in the LF image, so we utilize the GED considering all MLIs in the LF image. We will propose our Content Pooling method in section \[sec:pooling\].

### **Uniform Local Binary Pattern (ULBP)**

Although image entropy and frequency entropy can measure the distortion caused by angular inconsistency, $E_{MLI}$ is the global characteristics of the whole MLI. As shown in Fig. \[fig:MLIS\], the increased angular distortion changes the local texture of the MLI. Thus, we utilize ULBP to measure local textural features in MLI.

Local binary patterns (LBP) has been proved as an efficient operator to extract local distribution information [@ojala1994performance; @ojala1996comparative]. LBP is a very simple but efficient texture descriptor, with rotation invariance, position invariance and robustness under various illumination. Since LBP can efficiently represent local distributions of joint pixels, we adopt modified ULBP descriptor to describe inconsistency of local adjacent pixels.

The LBP operator can be given as [@ojala2002multiresolution] $$LBP_{P,R}=\sum_{p=0}^{P-1}s(g_{p}-g_{c})2^{p}.$$

Here we set a neighbor of $P=4$ members on a circle of radius $R=1$. $g_{c}$ is the gray value of central pixel, while $g_{p}$ is the gray value of the neighbor pixel. $s$ represents the sign function.

To reduce the number of pattern types, we utilize modified Uniform LBP operator on each MLI$_{i}$. There exist $P+2$ types of uniform pattern class without regard to upright property. Then we combine the Probability of every Type of binary patterns ($PT$) into vector $\bm{PT}_{i}$. $$\bm{PT}_{i}=\{\underbrace{PT_{1},PT_{2},\cdots,PT_{n}}_{P+2}\}.$$

### **Content Pooling** {#sec:pooling}

In order to evaluate LF image quality, we pool the characteristics $E_{MLI}$ and $\bm{PT}_{i}$ of each MLI together into LF image angular features.

The method of pooling used on $E_{MLI}$ is percentile pooling [@moorthy2009visual]: central elements of $E_{MLI}$ are extracted while extremely large or small elements are neglected. Our experimental results verify that percentile pooling improves our proposed model.

We reserve 60% of central elements of $E_{MLI}$ here and show the global distribution histogram of IE and FE in Fig. \[fig:GEDhistogram\]. In accord with our analysis of Table \[tab:EMLI\], with the increase of angular distortion, the histogram of IE has a left shift, while the histogram of FE has a right shift. The higher level the distortion is, the larger the shift will become and the steeper the curve will appear.

![Histograms of GED for different types and levels of angular distortion. (a) IE histogram; (b) FE histogram.[]{data-label="fig:GEDhistogram"}](GEDhistogram.png){width="\linewidth"}

The mean values and skew values of IE and FE distribution histogram are selected as the GED features. $$f_{GED}=\{ mean(IE), skew(IE), mean(FE), skew(FE) \}.$$

The angular inconsistency of distorted LF images appears obviously at the edge between the foreground and background, while it appears inconspicuously in the background or mild areas. In these mild MLIs, less information is contained and the ULBP features might be misleading. Therefore, in order to exclude the MLIs in mild areas, we introduce a selector on $\bm{PT}_{i}$. The ULBP features of a whole LF image are exacted as: $$f_{ULBP}=avg\{\bm{PT}_{i}\} \quad \forall. i \quad if \quad R(i)>threshold,$$ where $R(i)=Max(MLI_{i})-Min(MLI_{i})$, the Range of a single MLI$_{i}$. $threshold$ is selected as gray value 20 here.

**Spatial Quality**
-------------------

The spatial quality plays an important part in LF image perceptual quality. Specifically, we utilize information entropy distribution of SAI to measure the changes in spatial quality. In the undistorted SAI, there exists the spatial dependence between adjacent pixels [@liu2014no]. With increased spatial distortion such as compression, the dependence is destroyed, resulting in changes of information entropy.

We divide the SAI into $8 \times 8$ blocks, then compute SAI Image Entropy (SIE) and SAI Frequency Entropy (SFE) of each block. Finally, pooling all the blocks together like what is shown in section \[sec:pooling\] and obtain Spatial Quality features. $$\begin{aligned}
f_{SQ}=\{ \ mean(SIE),skew(SIE), \\
    mean(SFE),skew(SFE) \ \}. 
\end{aligned}$$

**Experiment results** {#experiment}
======================

**Light Field Image Databases**
-------------------------------

To test the performance of our proposed LF-QMLI model, comparison experiments were conducted on Win5-LID [@shi2018perceptual] and VALID [@viola2018valid] databases. In Win5-LID database, there exist 220 distorted LF images with various distortion types and levels. The distortion types consist of JPEG2000, HEVC, linear interpolation (LN), nearest neighbor interpolation (NN) and two CNN models. The overall Mean Opinion Score (MOS) value is provided for each LF image.

VALID database includes 5 reference LF images captured by Lytro Illum with a number of distorted LF images caused by several types of compression methods. In this experiment, we utilize several distorted LF images obtained through the interactive methodology, including 40 LF images with two types of distortion: HEVC and VP9.

**Comparison with Previous Objective Metrics**
----------------------------------------------

We conducted comparison experiments between our proposed model and several FR, RR and NR metrics, including nine 2D-FR metrics [@SSIM; @MSSIM; @FSIM; @IW; @IFC; @VIF; @NQM; @VSNR], three 2D-NR metrics [@BRI; @NIQE; @FRI], one 3D-FR metric [@Chen], two 3D-NR metrics [@SINQ; @BSVQE], five Multi-views FR metrics [@MP_F; @MP_R; @MW; @DSwIM], one Multi-views NR metric [@APT] and one LF-RR metric [@LFIQM]. Three evaluation criteria are selected to measure the correlation between MOS and predicted results, consisting of Spearman Rank Order Correlation Coefficient (SROCC), Linear Correlation Coefficient (LCC) and Root Mean Squared Error (RMSE). The SROCC measures the monotonicity while LCC evaluates the linear relationship between predicted score and MOS. The RMSE computes the deviation of prediction. The better consistency with human perception is reflected in SROCC and LCC closing to 1 as well as RMSE closing to 0.

Then, we use SVR for regression [@smola2004tutorial]. LIBSVM package [@chang2011libsvm] is utilized to implement the SVR, which uses a radial basis function (RBF) kernel. We randomly select 80% of the database as the training set while the remaining 20% constitute the test set. The median of correlation coefficients across 1000 random trails were regarded as the final results.

The results of all metrics are shown in Table \[tab:performance\]. Here we can find that almost all 2D-FR metrics perform well on VALID database and LF-QMLI is competitive in NR metrics. This situation may be caused by its limited types of distortion. VALID only introduces two compression distortions, which destroy the spatial quality of LF images but do not take angular consistency into consideration, so previous 2D-FR metrics can almost excellently measure the distortion. The results show that VALID is not challenging for quality assessment of LF images. Therefore, we mainly analyze how our proposed model LF-QMLI performs on Win5-LID database.

On Win5-LID database, LF-QMLI outperforms all previous metrics. In general, the existing 2D and 3D metrics only consider the degradation of spatial quality, ignoring the degradation of angular consistency. Although multi-view metrics can measure the angular distortion, they do not take the compression distortion and similar spatial distortions into account. Therefore, we can reach a conclusion that the proposed model LF-QMLI can evaluate both angular consistency and spatial quality.

**Ablation Study**
------------------

In order to verify the validity of our proposed MLI-based model, we conducted an ablation study on Win5-LID database and the results are demonstrated in TABLE \[tab:ablation\]. The features extracted from MLI $f_{MLI}$ can obviously improve the model performance. One possible reason is that the $f_{MLI}$ provides the measurement of the 2-D angular consistency.

  -------------------------------------------------------------------------
                                                        SROCC   LCC   RMSE
  ---------------------------------------------------- ------- ----- ------
   **Model-f$ \bf _{MLI}$ & 0.6927 & 0.7890 & 0.5504\                
          **Model & 0.8802 & 0.9038 & 0.4147\                        
                          ****                                       
  -------------------------------------------------------------------------

  : Ablation Study[]{data-label="tab:ablation"}

**Conclusion**
==============

In this paper, we proposed a No-Reference Light Field image Quality assessment model based on Micro-Lens Image (LF-QMLI). We theoretically analyze the significance of MLI in LF-IQA and extract features for our evaluator. The model can effectively measure the angular consistency and spatial quality. The results show that LF-QMLI achieves state-of-the-art performance. In the future, we will consider more advanced features on the MLI to improve our model.

acknowledgment {#acknowledgment .unnumbered}
==============

This work was supported in part by NSFC under Grant 61571413, 61632001.

[^1]: \* Equal contribution
