Next Article in Journal
Towards a Connected Mobile Cataract Screening System: A Future Approach
Next Article in Special Issue
Atmospheric Correction for High-Resolution Shape from Shading on Mars
Previous Article in Journal
X-ray Tomography Unveils the Construction Technique of Un-Montu’s Egyptian Coffin (Early 26th Dynasty)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

DRM-Based Colour Photometric Stereo Using Diffuse-Specular Separation for Non-Lambertian Surfaces

by
Boren Li
1,† and
Tomonari Furukawa
2,*,†,‡
1
Beijing Institute for General Artificial Intelligence, Beijing 100124, China
2
School of Engineering and Applied Science, University of Virginia, Charlottesville, VA 22904, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Current address: 675 Old Reservoir Rd, Charlottesville, VA 22903, USA.
J. Imaging 2022, 8(2), 40; https://doi.org/10.3390/jimaging8020040
Submission received: 29 November 2021 / Revised: 3 January 2022 / Accepted: 5 January 2022 / Published: 8 February 2022
(This article belongs to the Special Issue Photometric Stereo)

Abstract

:
This paper presents a photometric stereo (PS) method based on the dichromatic reflectance model (DRM) using colour images. The proposed method estimates surface orientations for surfaces with non-Lambertian reflectance using diffuse-specular separation and contains two steps. The first step, referred to as diffuse-specular separation, initialises surface orientations in a specular invariant colour subspace and further separates the diffuse and specular components in the RGB space. In the second step, the surface orientations are refined by first initialising specular parameters via solving a log-linear regression problem owing to the separation and then fitting the DRM using Levenburg-Marquardt algorithm. Since reliable information from diffuse reflection free from specularities is adopted in the initialisations, the proposed method is robust and feasible with less observations. At pixels where dense non-Lambertian reflectances appear, signals from specularities are exploited to refine the surface orientations and the additionally acquired specular parameters are potentially valuable for more applications, such as digital relighting. The effectiveness of the newly proposed surface normal refinement step was evaluated and the accuracy in estimating surface orientations was enhanced around 30 % on average by including this step. The proposed method was also proven effective in an experiment using synthetic input images comprised of twenty-four different reflectances of dielectric materals. A comparison with nine other PS methods on five representative datasets further prove the validity of the proposed method.

1. Introduction

Photometric stereo (PS) estimates surface orientations using images captured from a fixed viewpoint under various illuminations and is especially powerful in acquisition of fine surface details at pixel level [1]. Surface orientation is important in a variety of fields, such as geometry segmentation for three-dimensional (3D) object recognition [2] and digital re-rendering in computer graphics. Surface geometries, which can be obtained via integrating surface orientations, have also been proven useful for applications, such as industrial quality control and reverse engineering. Due to the strength of PS and the importance of surface orientation acquisition, PS has drawn increasing interests since its debut [3]. However, making PS for a general real scene remains challenging due to the diverse reflectance properties of different materials that appear non-Lambertian [4]. This has given rise to the need for reliable estimation of surface orientations for a wide range of non-Lambertian reflectance, which essentially requires a proper imaging photometry model characterizing the forward problem and a subsequent PS method that inversely derives surface orientations.
Existing PS methods dealing with non-Lambertian reflectance can be classified into three categories. The first approximates surface reflectance using analytical bidirectional reflectance distribution function (BRDF) [5] and formulates the estimation of surface orientations as a nonlinear fitting problem. Nayar et al. [6] derived surface orientations using the Torrance-Sparrow BRDF [7] and incorporated extended sources to ensure sufficient information from specularities. Georghiades [8] inverted the same BRDF to simultaneously estimate surface orientations and resolve the generalized bas-relief ambiguity. Goldman et al. [9] assumed that general material reflectance can be represented by a convex combination of fundamental materials characterized by the Ward BRDF [10] and recovered surface orientations, fundamental material BRDFs and weight maps simultaneously for further scene editing purpose. Methods in this category exploiting information from surface reflectance are capable to derive not only surface orientations but also the other parameters in the analytical BRDFs, allowing more functionalities, such as digital relighting and material classification. However, due to the nonlinearity of analytical BRDFs and larger number of parameters to be estimated, these methods are sensitive to initialisations, numerically unstable under heavily corrupted outliers (e.g., shadows), and inapplicable when the number of observations is limited.
The second category of methods infers surface orientations through adopting the general properties of BRDF, such as isotropy, monotonicity and reciprocity. Alldrin et al. [11] developed a non-parametric PS method using bi-variate approximation of the isotropy property. Higo et al. [12] analyzed the general BRDF constraints and employed the BRDF properties of monotonicity, isotropy and visibility to vote for the most possible surface orientations for single-lobed reflectance. Shi et al. [13] proposed a bi-polynomial representation for low-frequency reflectance that was especially adaptable to the inverse problem as PS, while Ikehata et al. [14] developed another general isotropic BRDF as sum of lobes with unknown center directions. Methods in this category capitalize on the most fundamental properties of BRDF and therefore have the potential to deal with a broader range of reflectance. However, these methods are only capable to derive surface orientations with limited other functionalities and require an even larger set of observations compared with the first category. To widen the applicability, recent years have seen non-parametric BRDFs based on machine learning. Santo et al. [15,16] used a deep neural network for the first time whereas Taniai and Maehara [17] estimated surface normals and BRDFs by unsupervised learning. Ikehata [18] estimated surface normals more straightforward by deriving the so-called observation maps and using convolutional neural networks. Such data-driven approaches clearly improves the accuracy as observations directly create the BRDF or estimate surface normals. The use of a significantly large set of observations, which is key to the accuracy, and the expensive training effort using the observations make the approaches outside the scope of this paper.
Methods in the third category assume that non-Lambertian effects appear sparsely among observations and treat them as outliers. A substantial corpus of methods rely on robust statistical techniques for outlier rejection. Wu et al. [19] formulated the PS with outlier rejection as a global rank minimization problem, while Ikehata et al. [20] employed sparse Bayesian regression instead. Barsky et al. [21] initiated another line of PS researches using colour images where they first identified the significance of using specular colour in the dichromatic reflectance model (DRM) [22] for specularity rejection. Using the cue from known specular colour under the same theoretical foundation, Zickler et al. [23] derived a PS method in a novel two-dimensional (2D) specular invariant colour subspace. The major advantages of methods in the third category are their robustness and requirement for less images, while they are inefficient in the presence of dense non-Lambertian effects.
This paper presents a colour PS method to estimate surface orientations using diffuse-specular separation dealing with non-Lambertian reflectance. The proposed method models the colour imaging photometry using DRM as the forward problem and inverts it following a two-step procedure, the diffuse-specular separation and the surface normal refinement. The first step, initialising surface orientations using known specular colour similarly to [23], separates the diffuse and specular components in the RGB space and identifies outliers to reject in the UV space. In the second step, the parameters characterizing the specularity are initialised via solving a log-linear regression problem and the surface orientations are finally refined by fitting the nonlinear DRM using Levenburg-Marquardt algorithm (LMA). The proposed method robustly initialises surface orientations in a specular-free colour subspace and further fits for the nonlinear DRM with a newly derived parameter initialisation strategy. The proposed method preserves the advantage of robustness as the third category of methods, while it also benefits from the separated specular component to tackle with dense non-Lambertian reflectance as the first category. Furthermore, DRM parameters besides the surface orientations can be additionally obtained where dense non-Lambertian reflectance appears, making more potential applications feasible, such as digital relighting and material classification.
The rest of the paper is organized as follows. The next section presents the problem formulation of colour PS incorporating DRM and the relevant PS methods. Section 3 first overviews the flow and original contribution of the proposed colour DRM-based PS method and then, elucidates the specific two steps in the subsequent subsections. It follows the comprehensive evaluations on surface orientation estimation using both synthetic and real images in Section 4, while the last section summarizes the conclusions and proposes the future works.

2. Colour PS Incorporating Dichromatic Reflectance Model

2.1. Generic Colour PS Problem Formulation

Figure 1 shows the schematic diagram of the hardware components of the sensor and the necessary data for the generic colour PS. The hardware components are an RGB camera and a set of directional point sources whose positions are priorly known. Assume that ambient light is completely blocked, the colour PS requires that the RGB camera captures a colour image for the target surface with one source lit at a time. Given the N colour images illuminated under N different sources, the objective of the colour PS is to identify surface orientations pixel-wise. More specifically, given the colour image irradiance measurements at pixel i , j , { e k i , j * } k 1 , N , under the unit illumination directions, { l ¯ k i , j } k 1 , N , estimated from the light positions, the colour PS aims to derive the unit surface normal, n ¯ i , j , by minimizing the sum of N corresponding square residuals:
ϵ i , j = k = 1 N e k i , j * T w f n ¯ i , j , l ¯ k i , j 2 min n ¯ i , j .
where w is a vector of weights for colour channels whereas f is a forward function of n ¯ i , j having l ¯ k i , j as parameters, which corresponds to the weighted measurements. Note that ( ) ¯ and ( ) * mean the unit version and the measurement of ( ) respectively throughout the paper.

2.2. Colour Image Formation Model

Figure 2 shows the coordinate setups that are necessary for the colour PS to model the imaging geometry and the light configuration. The coordinate frames are the 3D camera frame, { C } , the 2D pixel frame, { P } , and the world frame, { W } , where its XY plane is known and referred to as the reference plane. These coordinate frames follow the convention given in [24]. A 3D surface point, { C } X ˜ , is projected to be a 2D image point, { P } x , in the RGB camera using the perspective projection as:
{ P } x = 1 { C } Z M p K ^ { C } X ˜ = 1 { C } Z m x 0 o x 0 m y o y 0 0 1 f 0 0 0 f 0 0 0 1 { C } X ˜ ,
where m x and m y are each the number of pixels per unit distance along x and y direction respectively. o x and o y are the principal point location in { P } , and f is the effective focal length. Note that { P } x = { P } x , { P } y , 1 and { C } X ˜ = { C } X , { C } Y , { C } Z . { P } x is rasterized into a regular grid in the camera as:
j i = ceil { P } x + 0.5 { P } y + 0.5 ,
where ceil ( · ) denotes the ceiling operator. { W } and { C } are related with a 3D rigid transformation: { C } X ˜ = R W { W } X ˜ + t W .
Let the kth light position be { C } S ˜ k = { C } X k , { C } Y k , { C } Z k . { C } S ˜ k is associated in the spherical space as:
{ C } X k = { C } r k sin { C } θ k cos { C } ϕ k { C } Y k = { C } r k sin { C } θ k sin { C } ϕ k { C } Z k = { C } r k cos { C } θ k ,
where { C } r k , { C } θ k and { C } ϕ k are the radius, zenith and azimuth angle, respectively. The N lights are configured such that they have the same zenith angle and radius, and uniformly distributed azimuth angles around the optical axis, i.e., { C } θ k = θ , { C } r k = r and { C } ϕ k = 360 ( k 1 ) N + ϕ 0 .
Having established the imaging geometry and light configuration, the imaging photometry is then modeled using the physically-based DRM [22]. The DRM properly characterizes reflectance for dielectric materials [25] and is represented by a linear combination of the diffuse and specular components. Assume that pixel i , j is shadow-free, the image irradiance for colour channel c after light strength normalisation and vignetting correction is written as:
e c i , j = e d , c i , j + e s , c i , j = d c i , j f d i , j + s c f s i , j ,
where e d , c i , j and e s , c i , j are the diffuse and specular components, and d c i , j and s c are diffuse and specular colours which are wavelength-dependent and represented by:
d c i , j = 0 Φ λ R i , j λ C c λ d λ , S c = 0 Φ λ C c λ d λ .
In the equation, λ is the wavelength, Φ λ is the source spectral power density (SPD), R i , j λ is the spectral body reflectance, and C c λ is the camera spectral sensitivity for the colour channel. f d i , j and f s i , j represent the diffuse and specular geometrical scaling factors that are characterized by BRDF.
Since the scaling factors are illumination-dependent, let the image irradiance Equation (5) generalized in the RGB colour space be with the kth illumination as:
e k i , j = e d , k i , j + e s , k i , j = d ¯ R G B i , j f d , k i , j + s ¯ R G B f s , k i , j = d ¯ R i , j d ¯ G i , j d ¯ B i , j f d , k i , j + s ¯ R s ¯ G s ¯ B f s , k i , j ,
where d ¯ R G B i , j = { d c i , j } c = R , G , B and s ¯ R G B = { s c } c = R , G , B are the unit diffuse and specular colour, respectively. The diffuse geometrical scaling factor is characterized with the kth unit illumination direction l ¯ k i , j by Lambertian model as:
f d , k i , j = k d i , j n ¯ i , j l ¯ k i , j ,
where k d i , j is the diffuse reflectance factor. Since the surface deviation is small compared with the light-surface distance, l ¯ k i , j can be approximately derived using a vector from the point intersected by ray i , j and the reference plane to { C } S ˜ k . One of the common options of the specular geometrical scaling factor, on the other hand, is a model using the Blinn-Phong BRDF [26] because of its computational efficiency and high performance in characterizing a wide range of isotropic reflectance [27]:
f s , k i , j = k s i , j n ¯ i , j T h ¯ k i , j β i , j ,
where k s i , j is the specular reflectance factor, h ¯ k i , j is the unit half vector, and β i , j is the shininess coefficient. The unit half vector is written as:
h ¯ k i , j = l ¯ k i , j + v ¯ i , j l ¯ k i , j + v ¯ i , j ,
where v ¯ i , j is the unit viewer direction and determined by:
v i , j = j o x m x , i o y m y , f 𐊗 , v ¯ i , j = v i , j / v i , j .

2.3. Colour PS Methods Using DRM

Existing colour PS methods using DRM contain two major variants, representative works of which were developed by Barsky et al. [21] and by Zickler et al. [23]. The PS method of Barsky finds the unit surface normal n ¯ i , j and the diffuse reflectance factor k d i , j such that the residual between the colour image irradiances weighted by d ¯ R G B i , j and the diffuse geometrical scaling factor is minimised:
ϵ i , j = k = 1 N b e k i , j * d ¯ R G B i , j f d , k i , j 2 = k = 1 N b e k i , j * d ¯ R G B i , j k d i , j n ¯ i , j l ¯ k i , j 2 min n ¯ i , j , k d i , j ,
where N b ( N ) is the number of specular-free and shadow-free images. The reason for the use of N b instead of N is that the Barsky’s PS method assumes that non-Lambertian reflectance appears sparsely among observations, which should be rejected as outliers. As the first process for outlier rejection, the direct principal component analysis (DPCA) [28] on E R G B i , j * E R G B i , j * is performed to estimate d ¯ R G B i , j , where E R G B i , j * = e 1 i , j * , e 2 i , j * , , e N i , j * . Given s ¯ R G B and d ¯ R G B i , j , the method then determines f s , k i , j as:
f s , k i , j = e k i , j * T s ¯ R G B e k i , j * d ¯ R G B i , j d ¯ R G B i , j s ¯ R G B 1 d ¯ R G B i , j T s ¯ R G B 2 .
A threshold on f s , k i , j is finally used to determine pixels to be rejected as outliers, and this subsequently determines N b .
The objective function (12) of this PS method is not different from that of the conventional Lambertian-based PS method [3]. Given { l ¯ k i , j } k 1 , M a priori, n ¯ i , j together with k d i , j , which minimizes the objective function (12), is analytically determined as:
k d i , j = L i , j L i , j 1 L i , j E R G B i , j * d ¯ R G B i , j , n ¯ i , j = L i , j L i , j 1 L i , j E R G B i , j * d ¯ R G B i , j k d i , j ,
where L i , j = l ¯ 1 i , j , l ¯ 2 i , j , , l ¯ N b i , j and rank ( L i , j ) = 3 . While the effect of specularity rejection through the use of s ¯ R G B is exhibited, the estimation of d ¯ R G B i , j using DPCA is erroneous due to the existence of specularities. Furthermore, all the specularities are rejected as outliers though useful information may be contained.
Unlike Barsky’s PS method, which performs PS in the RGB colour space, the PS method of Zickler estimates surface orientations in the SUV colour space, or the UV colour space as S is not considered:
ϵ i , j = k = 1 N e U V , k i , j * T d ¯ U V i , j α i , j n ¯ i , j T l ¯ k i , j 2 min n ¯ i , j ,
where α i , j is a scaling factor. The UV colour subspace includes more observations than the RGB colour space in surface orientation estimation. To operate in the SUV space, the colour image irradiance is transformed from RGB to SUV as:
e S U V , k i , j e k i , j = R s e k i , j = d ¯ S U V i , j f d , k i , j + s ¯ S U V f s , k i , j ,
where d ¯ S U V i , j = d ¯ U i , j , d ¯ V i , j , d ¯ S i , j T , s ¯ S U V = 0 , 0 , 1 T , and R s S O ( 3 ) is any transformation that yields s ¯ S U V = R s s ¯ R G B = 0 , 0 , 1 T . UV forms a specular-free colour subspace that is invariant to f s , k i , j . The Zickler’s PS method first estimates d ¯ U V i , j of unit length using DPCA on E U V i , j * E U V i , j * . It then becomes the conventional Lambertian-based PS problem and the solution is obtained in the similar form to Equation (14). Compared with the Barsky’s PS, the Zickler’s PS utilizes more observations and advantageously performs PS in the specular-free UV colour subspace. Shape information from specularities is however neglected again, and the method can only provide diffuse component up to the normalised RGB space. The diffuse-specular separation in the RGB space remains open.

3. DRM-Based Colour PS Method Using Diffuse-Specular Separation

3.1. Overview of the Proposed Color PS Method

Figure 3 shows the flow and original contributions of the proposed DRM-based colour PS method using the diffuse-specular separation. The proposed method consists of two processes: Step 1, which is diffuse-specular separation using a PS, and Step 2, which is the surface normal refinement. Step 1 identifies good approximate surface normals, whereas Step 2 further refines and finalizes them with maximum accuracy. The steps are marked by red boxes.
Step 1 is composed of four sub-processes: diffuse colour estimation, diffuse-specular separation in RGB space, PS in UV space and the outlier estimation. After completing shadow rejection as a preprocess, and generating the shadow-free image irradiance matrix E ̲ R G B i , j * , Step 1 begins with the diffuse colour estimation using the robust principal component analysis (RPCA) proposed in this paper and derives d ¯ R G B i , j . Note that ( ) ̲ is a quantity of ( ) after shadow rejection. The diffuse color estimation also outputs a specularity map, which describes the distribution of specular reflection. The diffuse-specular separability is then checked for each pixel by deriving the diffuse-specular chromatic angle ψ i , j in RGB space. For separable pixels, the proposed method performs PS in the specular-free UV colour subspace similarly to Zickler’s PS method and derives the initial guess of the unit surface normal n ¯ i , j 1 and the diffuse reflectance factor k d i , j 1 similarly to Barsky’s PS. ( ) 1 means the initial guess of ( ) . In parallel, the specular geometrical scaling factors, f ̲ s i , j 1 , are derived through the diffuse-specular separation process using Equation (13) for the subsequent surface normal refinement. Observations not in the specularity map are clamped to zero in f ̲ s i , j 1 . The surface normals and the other parameters are then refined by fitting the nonlinear DRM using LMA with regularisation.
In Step 2, the initial guess of the specular parameters, k s i , j 1 and β i , j 1 , are drived from n ¯ i , j 1 and f ̲ s i , j 1 . Finally, PS for the surface normal refinement derives n ¯ i , j and k d i , j with enhanced accuracy while additionally updating k s i , j and β i , j . The strength of the proposed method lies in this PS with the diffuse-specular separation. For pixels where sparse non-Lambertian reflectances appear, the first step is sufficient to reliably estimate the surface orientations. To tackle with more general problems with dense non-Lambertian reflectances, the second step is additionally required to exploit surface normal information from specularities in the DRM with the designed parameter initialisation strategy. The rest of this section elucidates the proposed method in detail.

3.2. Diffuse-Specular Separation

3.2.1. Diffuse Color Estimation

The proposed RPCA derives d ¯ R G B i , j and the specularity map as follows:
1.
Perform principle component analysis (PCA) for E ̲ R G B i , j * E ̲ R G B i , j * to estimate d ¯ R G B i , j ;
2.
Compute residual matrix: R d i , j = E ̲ R G B i , j * d ¯ R G B i , j d ¯ R G B i , j E ̲ R G B i , j * .
3.
Compute residual vector: r d i , j = r d , 1 i , j 2 + r d , 2 i , j 2 + r d , 3 i , j 2 , where R d i , j = r d , 1 i , j , r d , 2 i , j , r d , 3 i , j and 2 denotes the hadamard square.
4.
If the mean of r d i , j , mean ( r d i , j ) , is smaller than the threshold T d , terminate RPCA and output the current estimate of d ¯ R G B i , j and the specularity map. Otherwise, find the element that provides the maximal value of r d i , j mean ( r d i , j ) std ( r d i , j ) and register ( i , j ) in the specularity map. Remove the corresponding row vector in E ̲ R G B i , j * and repeat from step 1.
The RPCA eliminates image irradiances that are inappropriate for diffuse colour estimation.
The validity of the RPCA can be explained as follows. The shadow-free image irradiance matrix E ̲ R G B i , j * is expanded from Equation (7) and decomposed as
E ̲ R G B i , j * = f ̲ d i , j d ¯ R G B i , j + f ̲ s i , j s ¯ R G B = f d , 1 i , j f d , 2 i , j f d , N p i , j d ¯ R i , j d ¯ G i , j d ¯ B i , j + f s , 1 i , j f s , 2 i , j f s , N p i , j s ¯ R s ¯ G s ¯ B
where most entries in f ̲ s i , j are near-zero since images are often near specular-free. In such a case, a non-zero residual occurs if f ̲ s i , j contains non-negligible entries. The RPCA algorithms robustly and accurately estimate d ¯ R G B i , j by keeping rejecting specularities until the residual becomes smaller than T d . Unlike Barsky’s PS method, which identifies specularities after d ¯ R G B i , j estimation, the proposed algorithms derive d ¯ R G B i , j and the specularity map simultaneously.
While the originality and superiority of the proposed RPCA algorithms have been identified, the performance of the RPCA algorithms depends on the value of T d . If T d is too large, the RPCA algorithms are tolerant to specularities and perform similarly to the DPCA algorithms. If T d is too small, the RPCA algorithms reject innocent pixels as specularities and make images sensitive to noise.

3.2.2. Diffuse-Specular Separability Check

Having d ¯ R G B i , j estimated, the diffuse-specular chromatic angle, used to perform the separability check, is derived as
cos ψ i , j = d ¯ R G B i , j T s ¯ R G B ,
Asymptotically, the chromatic angle gives
lim ψ i , j 0 d ¯ R G B i , j = s ¯ R G B
This means that ψ i , j should be sufficiently larger than 0 in order for d ¯ R G B i , j and s ¯ R G B to be distinct and separable. The threshold of the chromatic angle to determine the separability is T c . If ψ i , j T c , d ¯ R G B i , j and s ¯ R G B are sufficiently distinct and pixel i , j is considered as separable. Similarly to T d , T c should be chosen carefully. A small T c results in yielding false separation. Once d ¯ R G B i , j and s ¯ R G B have been identified, f ̲ s i , j 1 is then initialised using Equation (13).

3.2.3. PS in UV Space

Since the PS is performed in the UV colour space, the unit diffuse colour d ¯ R G B i , j should be converted into that in UV colour space d ¯ U V i , j . The PS in UV color space allows more observations as described in the Zickler’s PS method. With R s , the estimated d ¯ R G B i , j is transformed to the SUV space as d ^ i , j = d ^ U i , j , d ^ V i , j , d ^ S i , j . The diffuse color in UV space d ^ U V i , j is then given by
d ^ U V i , j = d ^ U i , j , d ^ V i , j = κ z i , j d ¯ U V i , j
where d ¯ U V i , j = d ¯ U i , j , d ¯ V i , j is a unit vector, and
κ z i , j = d ^ U i , j 2 + d ^ V i , j 2 > 0 .
Image irradiances of separable pixels are also transformed to the SUV colour space as E ̲ S U V i , j using the same R s . E ̲ U V i , j is then formed by picking the first two columns of E ̲ S U V i , j . Due to the presence of image noise, the PS in UV space is modified from Equation (15) to reject noise-corrupted image irradiances as outliers using studentised residuals [29]:
ϵ i , j = k = 1 N p e U V , k i , j * d ¯ U V i , j κ z i , j 2 f d , k i , j 2 = k = 1 N p e U V , k i , j * d ¯ U V i , j κ z i , j 2 k d i , j 1 n ¯ i , j 1 l ¯ k i , j 2 min n ¯ i , j 1 , k d i , j 1 ,
where N p is the number of shadow-free images, so N b N p N . The derivation of k d i , j 1 in addition to n ¯ i , j 1 similarly to the Barsky’s PS method enables full recovery of the colour diffuse component. N p is used instead of N and N b because shadow-free images are effective for reliability.
After the PS, n ¯ i , j and k d i , j are initialised as n ¯ i , j 1 and k d i , j 1 respectively. The diffuse component matrix, E d , R G B i , j , is then re-rendered as:
E d , R G B i , j = max k d i , j 1 L i , j n ¯ i , j 1 , 0 N d ¯ R G B i , j T ,
where 0 N is an N × 1 vector with all zero entries. The specular component matrix, E s , R G B i , j , is obtained by max f ̲ s i , j 1 s ¯ R G B T , 0 N . With the separated colour diffuse and specular components, the proposed method allows the functionality of specular removal and intrinsic image decomposition.

3.2.4. Outlier Estimation

Let the residual vector r ̲ i , j be written as:
r ̲ i , j = E ̲ U V i , j d ¯ U V i , j L ̲ i , j n i , j = E ̲ U V i , j d ¯ U V i , j L ̲ i , j ρ i , j n ¯ i , j ,
where n i , j = ρ i , j n ¯ i , j is the scaled normal, L ̲ i , j = l ¯ 1 i , j , l ¯ 2 i , j , , l ¯ M i , j T , and ρ i , j = k d i , j κ z i , j . The hat matrix H a i , j [29] is represented by:
H ̲ a i , j = L ̲ i , j L ̲ i , j T L ̲ i , j 1 L ̲ i , j T .
Diagonal entries in H ̲ a i , j are leverages and the kth leverage is denoted as h a , k i , j . These leverages quantify the influence that the observed E ̲ U V i , j d ¯ U V i , j on their predicted values L ̲ i , j n i , j . Each entry in the studentized residual vector r ̲ ˜ i , j is then approximated using:
r ˜ k i , j = r k i , j ( M S E ) i , j 1 h a , k i , j ,
where ( MSE ) i , j represents the mean squared error of r ̲ i , j . If any entry in | r ̲ ˜ i , j | exceeds T o , this entry is estimated as outliers. The proposed PS rejects outliers until no more can be detected or ( MSE ) i , j is smaller than a tolerance, T m . With the general rule of thumb in detecting outliers using studentised residuals, T o and T m can be chosen as 2.5 and 9 σ n 2 , where σ n is the standard deviation of additive white Gaussian noise that can be estimated using the method given by [30].

3.3. Surface Normal Refinement

3.3.1. Specular Parameter Initialisation

If dense non-Lambertian reflectance appears at pixel i , j , which implies that more than one entry in f ̲ s i , j 1 are in the specularity map, surface normal refinement is necessary since the sparse non-Lambertian reflectance assumption in the diffuse-specular separation is violated. With the estimated f ̲ s i , j 1 , the next objective is to initialise the specular parameters, k s i , j 1 and β i , j 1 . The specific cost functional is given by:
ϵ s i , j = k = 1 N q f s , k i , j 1 k s i , j 1 h ¯ k i , j n ¯ i , j 1 β i , j 1 2 min k s i , j 1 , β i , j 1 ,
where N q is the number of specularity maps showing pixel i , j and N q 2 . Equation (26) suggests a nonlinear least-squares problem, while it can be manipulated to the natural logarithm domain as a linear problem given by:
ϵ l s i , j = k = 1 N q ln f s , k i , j 1 ln k s i , j 1 β i , j 1 ln h ¯ k i , j n ¯ i , j 1 2 min k s i , j 1 , β i , j 1 .
Let f ̲ ̲ s i , j 1 with N q rows consist of entries in f ̲ s i , j 1 in the specularity map and H ̲ ̲ i , j = h ¯ 1 i , j , h ¯ 2 i , j , . . . , h ¯ N q i , j T be the matrix comprised of the corresponding unit half vectors. Then, the solution of ln k s i , j 1 and β i , j 1 to Equation (27) is given analytically by:
β i , j 1 ln k s i , j 1 = S ̲ ̲ i , j S ̲ ̲ i , j 1 S ̲ ̲ i , j ln f ̲ ̲ s i , j 1 ,
where
S ̲ ̲ i , j = ln H ̲ ̲ i , j n ¯ i , j 1 , 1 N q
and 1 N q is a N q × 1 vector with all entries equal to 1.

3.3.2. Surface Normal Refinement in DRM

All parameters in the DRM have been initialised up to this point. The LMA is then employed to iteratively refine these parameters by solving the optimization problem:
ϵ i , j = k = 1 N e k i , j s ¯ R G B k d i , j n ¯ i , j l ¯ k i , j d ¯ R G B i , j s ¯ R G B k s i , j n ¯ i , j h ¯ k i , j β i , j + T α 1 n ¯ i , j n ¯ i , j 2 min n ¯ i , j , k d i , j , k s i , j , β i , j ,
where the first three terms are from the generalised formulation, and, the fourth regularization term with constant T α is added. The regularisation term reduces the change of surface normal refinement from the initial guess of n ¯ i , j 1 and prevents overfitting. Larger value of T α strengthens the robustness but limits the capability of surface normal refinement.
The strength of the proposed PS method lies in the determination of the surface normals while separating the diffuse and specular components in RGB space, eliminating outliers and then tuning the diffuse and specular reflectance factors simultaneously. The determination considering diffuse and specular reflections simultaneously makes the surface normal estimation more accurately. In addition, the proposed method can be used for specular removal [31] and intrinsic image decomposition [32]. Parameters characterizing f s , k are determined at pixels, so the proposed method could also be used for digital relighting and material classification where dense non-Lambertian reflectance appears.

4. Performance Evaluation on Surface Orientation Estimation

4.1. Evaluations Using Synthetic Images

The first experiment aims to evaluate the effectiveness of surface normal refinement for dense non-Lambertian reflectance using synthetic input images. A scene with six different-coloured spheres were rendered under 32 illuminants using the Blinn-Phong model. The sphere was adopted as the scene geometry since it samples all the possibilities of surface orientations for the visible surface. Different colours of the spheres were introduced for more comprehensive evaluations of the proposed PS method on surfaces with various spectral reflectances. The 32 illuminants were chosen to provide a sufficient number of pixels with dense non-Lambertian reflectance for evaluation. s ¯ R G B was set at 0.5774 , 0.5774 , 0.5774 T . Three colours of the spheres were red, green and blue with the same ψ i , j of 57 . 74 , while the other three spheres were yellow, cyan and magenta with ψ i , j of 35 . 26 . The 32 light positions were configured with r = 442 mm and θ = 20 . The reference plane was located where R W = d i a g ( 1 , 1 , 1 ) and t W = 0 , 0 , 678 T . k d i , j , k s i , j and β i , j were set the same across the field of view as 0.4 , 0.2 and 100, respectively. Additive Gaussian noise was introduced with σ n = 0.02 . The image irradiance under the first illuminant is shown by Figure 4a. The parameters of the proposed PS method were configured as: T d = 0.01 , T c = 5 , T o = 2.5 , T m = 0.0036 , T α = 3 , T x = 0.02 and T y = 0.02 . 100 repeated tests were conducted and the mean value was adopted for evaluation.
Figure 4b shows the angular error of surface orientation estimation without surface normal refinement. As is shown, the central regions of the spheres where dense non-Lambertian reflectances appear have larger error. By including the proposed surface normal refinement step, the angular errors of surface orientations around the detected dense non-Lambertian reflectance region, shown by a ring, are significantly reduced as shown by Figure 4c. This is where the reflectance parameters are most tuned through the surface normal refinement. Errors in such regions are even smaller than those with sparse non-Lambertian reflectance, which further enhances that making specularities as meaningful signals is beneficial for surface orientation estimation. Figure 4d demonstrates the improvement percentage for quantitative evaluation. The improvement percentage is defined reduction of the error divided by the error without surface normal refinement. The accuracy of surface orientations was enhanced by 34.33 % in median and 32.25 % on average. The first and third quantile values of the improvement were 15.76 % and 54.23 % , respectively. The performance improvements onto the six different-coloured spheres are shown in Table 1. The improvements were obvious for all spheres is significant, implying that the method was applicable to a wide range of surfaces with different spectral reflectances. For spheres that had larger chromatic angle, ψ i , j , the d ¯ R G B i , j estimations were more accurate. The better estimated d ¯ R G B i , j led to the more accurate recovery of f ̲ s i , j 1 , resulting in the more reliable initialisation of n ¯ i , j 1 . Such initialisations of n ¯ i , j ultimately affected the different improvement performances. Overall, the accuracy was enhancement by including the additional surface normal refinement step and for spectral reflectances with larger chromatic angles, the improvement was more phenomenal.
In the second experiment, the performance improvement due to surface normal refinement was verified onto a wide range of reflectances. Experimental settings and parameters of the proposed method were the same as the first experiment except that different values of k d i , j , k s i , j and β i , j were applied to represent various surface reflectances when generating synthetic images. The ratio of κ d s i , j = k d i , j / k s i , j represents the relative strength between the diffuse and specular components, whereas β i , j indicates the width of the specular lobe. The mean, median, first and third quantile values of the improvement were adopted for evaluation and the angular error of d ¯ R G B i , j estimation was used for analyses. As shown from Table 2, the improvements on all the different reflectances were significant and more obvious for larger value of β i , j and κ d s i , j . Larger value of β i , j suggesting narrower specular lobe made the d ¯ R G B i , j estimation more accurate since the sparse non-Lambertian reflectance assumption was valid. Similarly, larger value of κ d s i , j indicating stronger diffuse components made better d ¯ R G B i , j estimation due to stronger inliers. More accurate d ¯ R G B i , j estimation led to bigger improvement from the surface normal refinement. In summary, with the additional surface normal refinement step, the accuracy of surface orientation estimation is enhanced by around 30 % on average and the improvement is more obvious for reflectance that has narrower specular lobe and stronger diffuse component.
The third experiment aims to evaluate the overall performance of the proposed colour PS method onto a wide variety of dielectric material reflectances. Twenty-four material BRDFs divided into six categories were evaluated, where the different BRDFs were provided from the MERL database [33]. Experimental settings and the method parameters were exactly the same as the first experiment except that the BRDFs were replaced. The image irradiance for the first six materials under the first illuminant is shown in Figure 5a. The angular errors of the surface orientation estimation are given in Figure 5b. As is shown, the error is not only larger at the near central regions where specularities overlap but also at the boundaries. Boundary pixels in the spheres have shallow grazing angles and the error is due to the lack of modeling of the fresnel reflection in the Blinn-Phong model. Figure 5c shows the performance of surface normal estimation on the 24 materials using box-and-whisker plot. The red and black dot represent the mean and median value, respectively. The lower and upper bound of the box indicates the first and third quantile values. The six colours of the boxes suggest the six material categories. As is shown, the proposed method can accurately estimate n ¯ i , j for most dielectric materials within 5 degrees on average. The estimation performances on phenolic, plastic and rubber were consistent, whereas those on the wood stain, fabric and acrylic were also acceptable with two exceptions, the violet acrylic and the green fabric. The degradation for violet acrylic was due to its wide specular lobe with small κ d s i , j , while the reason for green fabric was because of its small ψ i , j .

4.2. Evaluations Using Real Images

The proposed method was evaluated on five different datasets comprised of real images in the DiLigenT database [4]. Parameters of the proposed method were set the same as the first experiment. s ¯ R G B was estimated as 0.5774 , 0.5774 , 0.5774 T . From the first to fifth column in Figure 6, the results on BUDDHA, BEAR, POT2, READING and GOBLET are respectively shown. Figure 6a shows the image irradiance under the first illuminant, while Figure 6b,c demonstrate the estimated normal map and the angular error of surface orientation estimation, respectively. As is shown from READING, the proposed method is only feasible for surfaces whose d ¯ R G B i , j and s ¯ R G B are distinct. Figure 6c shows that the dense non-Lambertian problem where the specularity overlaps is mostly solved for BUDDHA, BEAR and POT2, while deficits still exist on GOBLET, which is because the DRM is not proper to characterize metallic reflectance.
The results were compared with those of nine other PS methods which are listed in Table 3. The results from the nine methods were provided from the DiLigenT database and these methods were representative and well-recognized. It is to be noted that the results of the recent machine learning based techniques reviewed in Introduction are not shown for comparison since the effect of the physics based models is investigated in this paper. Due to the data availability, refs. [21,23] were not involved in the comparison, while the essence of their works in using DRM with known s ¯ R G B was well-inherited by the proposed method.
Figure 6d shows the results from ten PS methods on the five datasets. The X-axis represents the method index and the Y-axis indicates the angular error of n ¯ i , j . The five colours of the bar represent the different datasets. The angular error of n ¯ i , j is also demonstrated using the box-and-whisker plot similar to Figure 5. As is shown, the proposed method performs well in estimating surface orientations for all the five datasets except GOBLET. The results suggest that modeling reflectance using DRM is proper for surfaces made of dielectric materials. The strength of the proposed method lies on the reliable estimation of d ¯ R G B i , j with simultaneous specularity detection. The advantage is also accredited by exploiting n ¯ i , j not only from the diffuse components but also from the specularities owing to the diffuse-specular separation. In spite of its advantages, the proposed method is not effective for metallic surfaces compared with methods, such as [14]. The limitation originates from the fundamental assumption of DRM that s ¯ R G B are the same across the surface. This assumption is violated for metallic surfaces since the wavelength and geometry exhibit inter-dependency.

5. Conclusions and Future Works

A PS method using colour images dealing with non-Lambertian reflectance has been proposed. The method formulates the imaging photometry using DRM with known specular colour. It extracts surface orientations not only from the diffuse components but also specularities owing to the diffuse-specular separation. Introducing the additional surface refinement step using information from specularities, the proposed method particularly improves the accuracy for surface orientation estimation at pixels where dense non-Lambertian reflectance appear. The simultaneously acquired specular parameters can be applied for more potential functionalities, such as digital relighting and material classification.
From the experiment of validating the effectiveness of the newly proposed surface normal refinement step on surface where dense non-Lambertian reflectance appears, the results indicate that with the additional step, the accuracy is enhanced by around 30 % on average and the improvement is more phenominal for surfaces with larger chromatic angles, stronger diffuse components and narrower specular lobes due to the better estimated unit diffuse colour. The result investigating the proposed method applicability on 24 different reflectances of dielectric materials suggests that the average angular error of surface orientation estimation for most materials are within 5 and bigger errors occur at surface patches with shallower grazing angles due to the limitation of the Blinn-Phong model that does not incorporate the fresnel reflection. From the systematic comparison on five datasets with nine other representative methods, the proposed method shows its descent performance on reflectances of dielectric materials and degradation on metallic surface due to the limitation of DRM. The result has also demonstrated that the proposed method is only feasible for surfaces whose diffuse and specular colour are distinct.
This paper presented the first set of results to estimate surface orientations using colour images for non-Lambertian reflectance, while it can be extended in a variety of ways. First, other analytical BRDFs can be employed instead of the Blinn-Phong model to better account for the fresnel reflection. Second, machine learning based models can be incorporated to learn and reduce errors of analytical models. Further, an alternative method can be proposed specifically for surfaces whose diffuse and specular colours are close to extend the method’s applicability. Making the method adaptive is one possible solution. Last but not least, it is also of particular interest to infer surface properties from the estimated specular parameters for the purpose of material classification.

Author Contributions

Conceptualization, B.L. and T.F.; methodology, B.L. and T.F.; software, B.L.; validation, B.L.; formal analysis, B.L.; investigation, B.L.; resources, B.L.; data curation, B.L.; writing—original draft preparation, B.L.; writing—review and editing, T.F.; visualization, B.L.; supervision, T.F.; project administration, T.F.; funding acquisition, T.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work is sponsored in part by US Office of Naval Research (N00014-20-1-2468).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to express their sincere gratitude to Tom McKenna at Office of Naval Research for his invaluable advice.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ackermann, J.; Goesele, M. A survey of photometric stereo techniques. J. Found. Trends Comput. Graph. Vis. 2015, 9, 149–254. [Google Scholar] [CrossRef]
  2. Klasing, K.; Althoff, D.; Wollherr, D.; Buss, M. Comparison of surface normal estimation methods for range sensing applications. In Proceedings of the 2009 IEEE International Conference on Robotics and Automation, Kobe, Japan, 12–17 May 2009; pp. 3206–3211. [Google Scholar]
  3. Woodham, R.J. Photometric method for determining surface orientation from multiple images. Opt. Eng. 1980, 19, 139–144. [Google Scholar] [CrossRef]
  4. Shi, B.; Zhe, W.; Mo, Z.; Duan, D.; Yeung, S.; Tan, P. A benchmark dataset and evaluation for non-Lambertian and uncalibrated photometric stereo. In Proceedings of the Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Las Vegas, NV, USA, 27–30 June 2016; pp. 3707–3716. [Google Scholar]
  5. Shell, J.R. Bidirectional Reflectance: An Overview with Remote Sensing Applications & Measurement Recommendations. 2004, pp. 37–50. Available online: http://web.gps.caltech.edu/~vijay/Papers/BRDF/shell-04.pdf (accessed on 4 January 2022).
  6. Nayar, S.K.; Ikeuchi, K.; Kanade, T. Determining shape and reflectance of hybrid surfaces by photometric sampling. IEEE Trans. Robot. Autom. 1990, 6, 418–431. [Google Scholar] [CrossRef]
  7. Torrance, K.E.; Sparrow, E.M. Theory for off-specular reflection from roughened surfaces. J. Opt. Soc. Am. 1967, 57, 1105–1114. [Google Scholar] [CrossRef]
  8. Georghiades, A. Incorporating the Torrance and Sparrow model of reflectance in uncalibrated photometric stereo. In Proceedings of the Proceedings Ninth IEEE International Conference on Computer Vision, Nice, France, 13–16 October 2003; pp. 816–823. [Google Scholar]
  9. Goldman, D.B.; Curless, B.; Hertzmann, A.; Seitz, S.M. Shape and spatially-varying BRDFs from photometric stereo. IEEE Trans. PAMI 2010, 32, 1060–1071. [Google Scholar] [CrossRef] [PubMed]
  10. Ward, G.J. Measuring and modeling anisotropic reflection. In SIGGRAPH92: Proceedings of the 19th Annual Conference and Exhibition on Computer Graphics and Interactive Techniques; Association for Computing Machinery: New York, NY, USA, 1992; pp. 265–272. [Google Scholar]
  11. Alldrin, N.G.; Zickler, T.; Kriegman, D.J. Photometric stereo with non-parametric and spatially-varying reflectance. In Proceedings of the 2008 IEEE Conference on Computer Vision and Pattern Recognition, Anchorage, AK, USA, 23–28 June 2008; pp. 1–8. [Google Scholar]
  12. Higo, T.; Matsushita, Y.; Ikeuchi, K. Consensus photometric stereo. In Proceedings of the 2010 IEEE Computer Society Conference on Computer Vision and Pattern Recognition, San Francisco, CA, USA, 13–18 June 2010; pp. 1157–1164. [Google Scholar]
  13. Shi, B.; Tan, P.; Matsushita, Y.; Ikeuchi, K. Bipolynomial modeling of low-frequency reflectances. IEEE Trans. PAMI 2014, 36, 1078–1091. [Google Scholar] [CrossRef] [PubMed]
  14. Ikehata, S.; Aizawa, K. Photometric stereo using constrained bivariate regression for general isotropic surfaces. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Columbus, OH, USA, 23–28 June 2014; pp. 23–28. [Google Scholar]
  15. Santo, H.; Samejima, M.; Sugano, Y.; Shi, B.; Matsushita, Y. Deep photometric stereo network. In Proceedings of the International Workshop on PBDL in Conjunction with IEEE ICCV, Venice, Italy, 22–29 October 2017; pp. 1–9. [Google Scholar]
  16. Santo, H.; Samejima, M.; Sugano, Y.; Shi, B.; Matsushita, Y. Deep photometric stereo networks for determining surface normal and reflectances. IEEE Trans. Pattern Anal. Mach. Intell. 2022, 44, 114–128. [Google Scholar] [CrossRef] [PubMed]
  17. Taniai, T.; Maehara, T. Neural inverse rendering for general reflectance photometric stereo. In Proceedings of the 35th International Conference on Machine Learning, Stockholm, Sweden, 10–15 July 2018. [Google Scholar]
  18. Ikehata, S. CNN-PS: CNN-based photometric stereo for general non-convex surfaces. In Proceedings of the European Conference on Computer Vision (ECCV), Munich, Germany, 8–14 September 2018; pp. 3–18. [Google Scholar]
  19. Wu, L.; Ganesh, A.; Shi, B.; Matsushita, Y.; Wang, Y.; Ma, Y. Robust photometric stereo via low-rank matrix completion and recovery. In Proceedings of the 10th Asian Conference on Computer Vision, Queenstown, New Zealand, 8–12 November 2010; pp. 703–717. [Google Scholar]
  20. Ikehata, S.; Wipf, D.; Matsushita, Y.; Aizawa, K. Robust photometric stereo using sparse regression. In Proceedings of the 2012 IEEE Conference on Computer Vision and Pattern Recognition, Providence, RI, USA, 16–21 June 2012; pp. 318–325. [Google Scholar]
  21. Barsky, S.; Petrou, M. The 4-source photometric stereo technique for three-dimensional surfaces in the presence of highlights and shadows. IEEE Trans. PAMI 2003, 25, 1239–1252. [Google Scholar] [CrossRef] [Green Version]
  22. Shafer, S.A. Using colour to separate reflection components. Color Res. Appl. 1985, 10, 210–218. [Google Scholar] [CrossRef] [Green Version]
  23. Zickler, T.; Mallick, S.P.; Kriegman, D.J.; Belhumeur, P.N. Color subspace as photometric invariants. Int. J. Comput. Vis. 2008, 79, 13–30. [Google Scholar] [CrossRef] [Green Version]
  24. Hartley, R.; Zisserman, A. Chapter 6: Camera Models. In Multiple View Geometry; Cambridge University Press: Cambridge, UK, 2003; pp. 153–177. [Google Scholar]
  25. Angelopoulou, E. Specular highlight detection based on the fresnel reflection coefficient. In Proceedings of the 2007 IEEE 11th International Conference on Computer Vision, Rio de Janeiro, Brazil, 14–21 October 2007; pp. 1–8. [Google Scholar]
  26. Blinn, J.F. Models of light reflection for computer synthesized pictures. In Proceedings of the 4th Annual Conference on Computer Graphics and Interactive Techniques, San Jose, CA, USA, 20–22 July 1977; pp. 192–198. [Google Scholar]
  27. Ngan, A.; Dur, F.; Matusik, W. Experimental validation of analytical BRDF models. In Proceedings of the ACM SIGGRAPH 2004 Sketches, Los Angeles, CA, USA, 8–12 August 2004; p. 90. [Google Scholar]
  28. Kao, Y.H.; Roy, B.V. Directed Principal Component Analysis. Oper. Res. 2014, 62, 957–972. [Google Scholar] [CrossRef] [Green Version]
  29. Pope, A.J. The Statics of Residuals and the Detection of Outliers; NOAA Technical Report; 1976; Volume 1, p. 65. Available online: https://www.ngs.noaa.gov/PUBS_LIB/TRNOS65NGS1.pdf (accessed on 4 January 2022).
  30. Liu, C.; Freeman, W.T.; Szeliski, R.; Kang, S.B. Noise estimation from a single image. In Proceedings of the 2006 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR’06), New York, NY, USA, 17–22 June 2006; pp. 901–908. [Google Scholar]
  31. Artusi, A.; Banterle, F.; Chetverikov, D. A survey of specularity removal methods. Comput. Graph. Forum 2011, 30, 2208–2230. [Google Scholar] [CrossRef]
  32. Barrow, H.; Tanenbaum, J. Recovering intrinsic scene characteristic from images. Proc. Comput. Vis. Syst. 1978, 2, 3–26. [Google Scholar]
  33. Matusik, W.; Pfister, H.; Br, M.; McMillan, L. A Data-Driven Reflectance Model. ACM Trans. Graph. 2003, 22, 759–769. [Google Scholar] [CrossRef]
  34. Shi, B.; Tan, P.; Matsushita, Y.; Ikeuchi, K. Elevation angle from reflectance monotonicity: Photometric stereo for general isotropic reflectances. In Proceedings of the European Conference on Computer Vision, Florence, Italy, 7–13 October 2012; pp. 455–468. [Google Scholar]
Figure 1. Schematic diagram for generic colour PS.
Figure 1. Schematic diagram for generic colour PS.
Jimaging 08 00040 g001
Figure 2. Coordinate setups for modeling imaging geometry and light configuration.
Figure 2. Coordinate setups for modeling imaging geometry and light configuration.
Jimaging 08 00040 g002
Figure 3. Flows and the original contribution of the DRM-based colour PS.
Figure 3. Flows and the original contribution of the DRM-based colour PS.
Jimaging 08 00040 g003
Figure 4. Evaluation of surface normal refinement: (a) Image irradiance under the first illuminant; (b) Angular error of surface orientations without surface normal refinement in degrees; (c) Angular error of surface orientations with surface normal refinement in degrees; (d) Improvement of surface orientation estimation by including surface normal refinement in percentage.
Figure 4. Evaluation of surface normal refinement: (a) Image irradiance under the first illuminant; (b) Angular error of surface orientations without surface normal refinement in degrees; (c) Angular error of surface orientations with surface normal refinement in degrees; (d) Improvement of surface orientation estimation by including surface normal refinement in percentage.
Jimaging 08 00040 g004
Figure 5. (a) Image irradiance for the first six material BRDFs under the first illuminant; (b) Angular error of surface orientation estimation for the first six material BRDFs; (c) Angular error of surface orientation for the twenty-four dielectric materals.
Figure 5. (a) Image irradiance for the first six material BRDFs under the first illuminant; (b) Angular error of surface orientation estimation for the first six material BRDFs; (c) Angular error of surface orientation for the twenty-four dielectric materals.
Jimaging 08 00040 g005
Figure 6. Evaluations on surface oreintation estimation: (a) Image irradiance under the first illuminant; (b) Estimated normal map; (c) Angular error of surface orientations; (d) Angular error of surface orientations on five datasets using ten different PS methods.
Figure 6. Evaluations on surface oreintation estimation: (a) Image irradiance under the first illuminant; (b) Estimated normal map; (c) Angular error of surface orientations; (d) Angular error of surface orientations on five datasets using ten different PS methods.
Jimaging 08 00040 g006
Table 1. Improvement on different-coloured spheres.
Table 1. Improvement on different-coloured spheres.
IndexColour ψ i , j Error of d ¯ RGB i , j Mean Improvement
1red 57 . 74 1 . 15 23.39 %
2yellow 35 . 26 1 . 33 32.85 %
3green 57 . 74 1 . 15 23.42 %
4cyan 35 . 26 1 . 32 32.04 %
5blue 57 . 74 1 . 15 23.15 %
6magenta 35 . 26 1 . 33 32.24 %
Table 2. Evaluation of the improvement for different reflectances.
Table 2. Evaluation of the improvement for different reflectances.
k d / k s β MeanMedianFirst QuantileThird QuantileError of d ¯ RGB
0.4 / 0.2 100 32.25 % 34.33 % 15.76 % 54.23 % 1 . 23
0.4 / 0.4 100 31.14 % 32.97 % 15.21 % 53.75 % 1 . 58
0.4 / 0.8 100 30.48 % 32.77 % 14.27 % 52.90 % 2 . 34
0.4 / 0.2 20 27.49 % 26.70 % 6.80 % 51.66 % 2 . 99
0.4 / 0.4 20 25.92 % 25.98 % 6.67 % 51.25 % 5 . 06
0.4 / 0.8 20 25.02 % 24.07 % 6.57 % 50.81 % 8 . 30
Table 3. Comparative studies of PS methods.
Table 3. Comparative studies of PS methods.
Method Index12345678910
Reference[3][9][11][12][19][20][13][34][14]proposed
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, B.; Furukawa, T. DRM-Based Colour Photometric Stereo Using Diffuse-Specular Separation for Non-Lambertian Surfaces. J. Imaging 2022, 8, 40. https://doi.org/10.3390/jimaging8020040

AMA Style

Li B, Furukawa T. DRM-Based Colour Photometric Stereo Using Diffuse-Specular Separation for Non-Lambertian Surfaces. Journal of Imaging. 2022; 8(2):40. https://doi.org/10.3390/jimaging8020040

Chicago/Turabian Style

Li, Boren, and Tomonari Furukawa. 2022. "DRM-Based Colour Photometric Stereo Using Diffuse-Specular Separation for Non-Lambertian Surfaces" Journal of Imaging 8, no. 2: 40. https://doi.org/10.3390/jimaging8020040

APA Style

Li, B., & Furukawa, T. (2022). DRM-Based Colour Photometric Stereo Using Diffuse-Specular Separation for Non-Lambertian Surfaces. Journal of Imaging, 8(2), 40. https://doi.org/10.3390/jimaging8020040

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop