=Paper= {{Paper |id=Vol-3271/Paper11_CVCS2022 |storemode=property |title=Characterising Appearance of Gold Foils and Gilding in Conservation and Restoration |pdfUrl=https://ceur-ws.org/Vol-3271/Paper11_CVCS2022.pdf |volume=Vol-3271 |authors=Yoko Arteaga,Aditya Sole,Jon Yngve Hardeberg,Clotilde Boust |dblpUrl=https://dblp.org/rec/conf/cvcs/ArteagaSHB22 }} ==Characterising Appearance of Gold Foils and Gilding in Conservation and Restoration== https://ceur-ws.org/Vol-3271/Paper11_CVCS2022.pdf
Characterising appearance of gold foils and gilding in
conservation and restoration
Yoko Arteaga1,2,* , Aditya Sole2 , Jon Yngve Hardeberg2 and Clotilde Boust1,3
1
  Centre for Research and Restoration of the Museums of France, Palais du Louvre, Porte des Lions, 14 quai Frano̧is
Mitterand, 75001 Paris, France
2
  Colourlab, Norwegian University of Science and Technology, Teknologivegen 22, 2815 Gjøvik, Norway
3
  PSL-PCMTH UMR8247, CNRS, 11 Rue Pierre et Marie Curie, 75005 Paris, France


                                         Abstract
                                         Materials with complex optical properties like gold foils are commonly used for gilding in the field of
                                         conservation and restoration to render the visual appearance of the restoration object. Characterising the
                                         appearance of these materials is important to render an accurate visual appearance that can be obtained
                                         with different gilding techniques. In this paper, the bidirectional reflectance of gold foils commonly used
                                         for gilding is characterised and as well as gilding mock-ups. The dataset of gold foils is publicly available.

                                         Keywords
                                         BRDF measurement, gold, material appearance




1. Introduction
Gold material is commonly found in artworks in the form of gilding. In the past, artists used
gold not only because of its cultural and religious importance but also because of its interesting
appearance properties. Gilding is a form of polychromy found in medieval paintings and
altarpieces [1]. It consists of a ground layer where a gold foil is attached to by different means.
Different gilding techniques render different ap pearances. For example, oil gilding is known to
have a matte and rough appearance whereas water gilding is smooth and glossy due to the final
burnishing process using an agate stone [2].
   In the field of conservation, it is important to have a deep understanding of the properties of
different materials and how they contribute to the original appearance of s urfaces. Gilding is
a composite material in which its elements are arranged in layers, gold being one of them. It
is thus, vital to characterise the appearance of the gold layer in the gilding to carry out better
conservation campaigns which can render a more faithful appearance of the gilded object.
   Gold material shows complex reflectance properties as its appearance changes with the change
in the illumination direction [3]. Characterising the appearance of such materials is therefore
difficult using the reflectance measurement geometries (like 0:diffuse, and 0:45) defined by the
The 11th Colour and Visual Computing Symposium, September 08–09, 2022, Gjøvik, Norway
*
 Corresponding author.
$ yoko.arteaga@culture.gouv.fr (Y. Arteaga); aditya.sole@ntnu.no (A. Sole); jon.hardeberg@ntnu.no
(J. Y. Hardeberg); clotilde.boust@culture.gouv.fr (C. Boust)
 0000-0002-6787-789X (Y. Arteaga); 0000-0002-7916-5363 (A. Sole); 0000-0003-1150-2498 (J. Y. Hardeberg);
0000-0001-6019-1118 (C. Boust)
                                       © 2022 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0).
    CEUR
    Workshop
    Proceedings
                  http://ceur-ws.org
                  ISSN 1613-0073
                                       CEUR Workshop Proceedings (CEUR-WS.org)
standardising bodies like ASTM and ISO. Bidirectional reflectance measurements are therefore
necessary to characterise these materials appearance. Bidirectional reflectance measurements
are performed using multi-angle spectrophotometers and goniospectrophotometers. Moreover,
various types of gold foils are used in gilding that vary in karat, composition, and provenance
which affects the gilding appearance visually. Famous examples besides traditional gold are
white gold and rose gold [1].
   In this paper, ten different gold foils commonly used for restoration, and two gilding mock-ups
made by professional restorers have been bidirectionally measured and characterised using a
bidirectional reflectance measurement device and a well established reflectance model. The
gilding mock-ups are produced using one of the gold foils presented in this paper.
   The structure of the paper is as follows: the next section of this paper presents a brief literature
review of appearance modelling of gold and gilding. The third section presents the gold foils
used and describes the fabrication method of the gilding mock-ups, the setup used for the
bidirectional measurements, and the reflectance model fitting. Finally, the results and discussion
are presented followed by the conclusion and final remarks.


2. Background
2.1. Appearance modelling of gold and gilding in restoration
Several attempts have been made to model and render the appearance of gold surfaces in cultural
heritage objects. Callet et al. [4] proposed an empirical model based on the metallic structure of
gold and used Kubelka Munk’s two flux theory to render the gilding of polychrome statues.
   The influence of the colour of the substrate on the appearance of gilding is actively studied
within the field of cultural heritage and has been tested empirically by Dumazet et al. [5] by
changing the colour of the substrate in their model and thus rendering different appearances of
gilded surface. A similar study has been carried out by Mounier et al. [6] on medieval mural
paintings. Mounier et al. [6] observed that the final appearance of the gilding is affected by the
substrate colour.
   Wu et al. [7] studied the influence of the substrate colour on the gild appearance by fabricating
gilding mock-ups using different methods such as oil gilding, water gilding and ground gilding.
Colorimetric and interferometric microscopic measurements were made to characterise the
mock-ups. They concluded that the colour of the substrate does not alter the visual appearance
of gilding. However, it was observed that burnishing the gild shows a significant visual change
in its appearance.
   MacDonald et al. [8] conducted a round-robin test to evaluate the quality of different multi-
and hyper-spectral imaging devices both in terms of spectral and spatial quality. They studied
objects with metal like finish and golden appearance. It was observed that the specularity of the
surface and the golden appearance of the object increase the difficulty of accurately reproducing
both colour and spectral characteristics of the object.
   In [9] an HDR hyperspectral imaging methodology was proposed to accurately capture the
reflectance spectra of an illuminated image with gold areas. The diffuse reflectance was then
rendered within an uniform dynamic range.
   A different type of gilding is found in medieval sculptures in Norway where instead of using a
gold foil, a silver foil is used and then a yellow resin is applied to it, which gives the appearance
of gold. In [10], bidirectional reflectance measurements of these samples are performed. It was
observed that the samples showed two specular peaks due to its composite structure. One being
from the interface air-resin and one from the interface resin-silver.

2.2. Bidirectional reflectance modelling
Bidirectional reflectance of a material can be modelled using a distribution function called
bidirectional reflectance distribution function (BRDF) defined by Nicodemus in [11]. A BRDF
describes how light is reflected off a surface. It is defined by the ratio between the differential
irradiance along l at point 𝑥 and the differential outgoing radiance along v:
                                                      𝑑𝐿𝑟 (𝑥, v)
                                      𝑓 (𝑥, l, v) =                                               (1)
                                                       𝑑𝐸(𝑥, l)
   For a BRDF model to be physically valid it must obey two properties: be energy conserving
and reciprocal. Due to the law of conservation of energy, the radiant energy reflected at a
surface must not be greater than the energy incident on the surface. The second property is
known as Helmholtz reciprocity, which dictates that the BRDF must be symmetrical. This means
that reversing incoming and outgoing light does not change the BRDF outcome.
   BRDF of a material can be further divided into isotropic and anisotropic. Regular surfaces with
a rotational symmetry are called isotropic and the BRDF is expressed as 𝑓 (𝑥, 𝜃𝑖 , 𝜃𝑟 ), following
the notation in Figure 4. Anisotropic surfaces exhibit a change in BRDF when rotated around
the surface normal, giving 𝑓 (𝑥, 𝜃𝑖 , 𝜑𝑖 , 𝜃𝑟 , 𝜑𝑟 ).
   A number of empirical and physical reflectance models have been defined till now to model
different types of materials and for applications [12]. Cook-Torrance (CT) is a well established
physical model and is used extensively to model specular materials like gold [3, 13, 14].

2.2.1. Cook-Torrance BRDF:
Cook-Torrance (CT) is a physical model based on the micro-facet theory [13]. CT model consists
of specular and diffuse reflection components and mathematically is defined as follows:

                                        𝑓 = 𝑘𝑑 𝑅𝑑 + 𝑘𝑠 𝑅𝑠 ,                                       (2)

In Eq. 2, 𝑘𝑑 and 𝑘𝑠 are weighting parameters and 𝑘𝑑 + 𝑘𝑠 = 1. The specular term is based on
micro-facet theory, which dictates that only the micro-facets on the surface with orientations in
between the viewing vector and the incident vector, called half-vector h = l.v, will contribute
to the reflected light.
   The general form of the CT specular term is:
                                              1 𝐹 𝐷𝐺
                                       𝑅𝑠 =                ,                                      (3)
                                              𝜋 (n.l)(n.v)
where 𝐹 is the Fresnel term [15], 𝐺 is the geometrical attenuation factor accounting for shad-
owing and masking of micro-facets, and 𝐷 is the distribution of normals facing h.
GGX distribution Walter et al. [16] propose a Bidirectional Scattering Distribution Function
(BSDF) called GGX. This GGX BRDF is inspired by the micro-facet model proposed by Cook
and Torrance. The general form of the GGX is given by:

                                             1 𝐹 𝐷𝐺
                                      𝑅𝑠 =                ,                                    (4)
                                             4 (n.l)(n.v)
which is the same as the general form of the Cook-Torrance in Eq. 3, except that the normalisation
factor for the Fresnel term is changed to 4. The Fresnel term in the GGX distribution is the
same as the one proposed by Cook and Torrance. The GGX normal distribution function of
micro-facets is given by:

                                             𝛼𝑔2 𝜒+ (h.n)
                                𝐷=                                ,                            (5)
                                      𝜋 cos4 𝜃ℎ (𝛼𝑔2 + tan2 𝛼ℎ )2
where 𝛼𝑔 is a width parameter for the specular lobe, 𝜃ℎ is the angle between the half-vector, h,
and the surface normal, n, and 𝜒+ is a positive characteristic function which equals to one if its
parameter is greater than zero and zero if its parameter is lesser or equal to zero.
  The geometrical attenuation factor which describes the shadowing and masking effects in
the GGX is derived from D. This is given by the monodirectional shadowing term defined as
                                       (︂     )︂
                                    + v.h                 2
                         𝐺1 (v) = 𝜒                 √︁                  ,                      (6)
                                          v.n 1 + 1 + 𝛼2 tan2 𝜃
                                                             𝑔        𝑣

where 𝜃𝑣 is the angle between the viewing vector v and the surface normal n.

2.3. BRDF metrics
When optimising a BRDF model, different metrics can be used to find the best approximation to
the surface reflectance. There is no consensus in the literature as to which metric should be
used. A common metric to use is the Root Mean Square error (RMSE). In [17], the BRDF of a set
of materials is estimated using different error metrics. The estimated BRDFs are then tested in a
psychophysical experiment to correlate the choice of error metric and perceived appearance.
The authors find that while metrics such as RMSE can provide a good optimisation when
evaluating the BRDF plot, this does not always correlate to a faithful appearance reproduction
if the rendered images are compared.


3. Method
3.1. Measurement samples: Gold foils
Ten gold foils that are commonly used in restoration were measured bidirectionally using a
multiangle spectrophotometer. Table 1 gives the name and specification of these 10 gold foils.
Since an extremely thin layer (approximately 0.2 microns in thickness) of the gold foils is used
for gilding, the foil samples were prepared by a professional restorer to allow for production
Table 1
Specifications of the gold foils used. ©C2RMF/Stéphanie Courtier.
                     Name            Colour                  Karat   Manufacturer
                     Or Alienor      Deep yellow/orange      23.5    Dauvet (FR)
                     Or Blanc        Light grey              12      Eytzinger (DE)
                     Or Citron       Light yellow            18      Eytzinger (DE)
                     Or Jaune        Deep yellow             21      Eytzinger (DE)
                     Or Noble        Deep yellow/orange      23.75   Eytzinger (DE)
                     Or Orange       Deep orange             22      Eytzinger (DE)
                     Or Rose         Rose                    22      Eytzinger (DE)
                     Or Rouge        Deep orange/red         23      Eytzinger (DE)
                     Or Versailles   Deep yellow             22      Dauvet (FR)
                     Or Vert         Light yellow/green      16.7    Eytzinger (DE)


accuracy, layer uniformity, sample transportation, and measurement possibility. The gold foils
have been prepared following the steps of oil gilding as follows:

    • A ground layer of glue is applied to an acrylic sheet,
    • When the glue is sticky, the gold foil is placed on top of the glue with the aid of a trowel
      and a gilder’s tip, illustrated in Figure 2a,
    • A gilder’s duster is used to tap on any air bubbles which may have formed between the
      gold sheet and the acrylic.

  Figure 1 shows the 10 gold foils used in this study.




(a) or Orange, or Noble, or Vert, or Rouge, or Rose, or Jaune, or Blanc, or (b) or Versailles used for the
    Citron, or Alienor.                                                         gilding mock-ups

Figure 1: Gold foil samples
3.1.1. Gilding samples:
Gilding using the gold foils can be carried out using 2 well established techniques called water
gilding, and oil gilding. In both the techniques, a wood piece surface is first prepared using
a ground layer made of chalk diluted with rabbit skin glue. In the water gilding technique, a
substrate layer of red coloured clay, called bole, is applied to the surface and is made wet using
ethanol and water before the gold foil is applied to it. In the next step, excess water (if any) is
removed using a gilder’s duster. The gold foil is then burnished manually with an agate stone
in a single direction of the surface plane. This process gives the surface a very high gloss. The
water gilding technique is therefore also known as "glossy gilding" technique.
   In the oil gilding technique, a layer of oil-based adhesive called ’gold-size’ is used as both,
the substrate and adhesive. When the gold-size is tacky but not completely dry, the gold foil is
applied to it using a gilder’s tip. Due to the nature of the gold size it is not possible to burnish
this surface, giving it a matte finish.




(a) Gold foil on a gilders cushion, cut with
    a trowel before being placed on the          (b) Gilding                           mock-ups.
    substrate using the gilder’s tip (pink)          Oil gilding on the left and water gilding
    ©C2RMF/Stéphanie Courtier                        on the right. ©Inp/Angèle Dequier



3.2. Bidirectional reflectance measurements
Bidirectional reflectance from the gold foil sample surface was measured using an X-Rite MA-
T12 multi-angle spectrophotometer. The MA-T12 instrument measures at six illumination
sources and two pick-up points according to the specifications of industry standards for metallic
colours [18, 19, 20]. It measures surface reflectance in the ranges of 400 - 700nm at 10nm
intervals. The illumination source is a polychromatic white LED with blue enhancement and its
spot size is of 9 mm by 12 mm. Figure 3 illustrates the combination of illumination and viewing
directions obtained by the MA-T12.




Figure 3: Measurement geometry of the MA-T12. Left: Illumination directions for viewing at -15∘ .
Right: Illumination directions for viewing at -45∘ .


   For each gold foil the measurements were repeated at five different locations, illustrated in
Figure 4. For each location, the MA-T12 was rotated clockwise around the normal, four times.
Since each measurement of the MA-T12 gives 12 spectra, 48 reflectance spectra were measured
per area at five different locations, giving a total of 240 reflectance spectra per sample. Table 2
shows all the combinations of incidence and viewing directions. The convention for the angles
is illustrated in Figure 4. The elevation angle (𝜃) ranges from -90∘ to 90∘ and is measured from
the normal of the surface and the azimuth angle (𝜑) ranges from from 0∘ to 360∘ and is measured
clockwise from the x-axis. Illumination directions are positive clockwise from the normal and
viewing directions are positive anti-clockwise from the normal.




Figure 4: Left: Diagram representing five areas measured for each gold foil: the four corners and the
centre. Right: Proposed parametrisation of the BRDF as a function of (𝜃𝑖 ) and (𝜃𝑟 ). For each area in the
gold foil 48 reflectance spectra are taken.


  The gilding mock-ups were measured in the same way as the gold foils with the exception
that the measurement was not repeated at different areas due to limited sample surface size.
Table 2
Elevation and azimuth angles for incident and viewing positions. The elevation angles are given from
the normal of the surface and the azimuth angles are given clockwise from the x-axis. All angles are in
degrees.
                        Illumination direction (𝜃𝑖 )     Viewing direction (𝜃𝑟 )
                                   -15 ∘
                                                        -60∘ -45∘ -30∘ 0∘ 30∘ 65∘
                                   -45∘                 -60∘ -30∘ -20∘ 0∘ 30∘ 65∘
                        Illumination azimuth (𝜑𝑖 )       Viewing azimuth (𝜑𝑟 )
                             0∘ 90∘ 180∘ 270∘               180∘ 270∘ 0∘ 90∘


3.3. BRDF fitting
Spectral reflectance measurements obtained from the MA-T12 were used to train the isotropic
CT model with the GGX distribution. The model was trained in the CIEXYZ colorimetric space.
CIEXYZ values were calculated using Equation (7). In Equation (7), 𝑟(𝜆) is the measured spectral
reflectance at the sample surface, 𝑆(𝜆) is the spectral power distribution of the D65 illuminant,
and 𝑥(𝜆), 𝑦(𝜆), 𝑧(𝜆) are the CIE 1964 colour matching functions,
                             ⎡ ⎤                        ⎡      ⎤
                               𝑋          ∫︁             𝑥(𝜆)
                             ⎣𝑌 ⎦ = 1        𝑆(𝜆)𝑟(𝜆) ⎣𝑦(𝜆)⎦ 𝑑𝜆,                              (7)
                                      𝑁 𝜆
                               𝑍                          𝑧(𝜆)

where                                           ∫︁
                                           𝑁=        𝑟(𝜆)𝑦(𝜆)𝑑𝜆.                                    (8)
                                                 𝜆
   The diffuse reflection component coefficient (𝑅𝑑 ) in the CT model was optimised individually
(𝑅𝑑𝑋 , 𝑅𝑑𝑌 , and 𝑅𝑑𝑍 ) in the CIEXYZ colour space while a single coefficient for the specular
reflection (𝑅𝑠 ) and roughness (𝑚) component was optimised across CIEXYZ as shown in
Equation (9).
                         ⎡      ⎤                  ⎛                   ⎡     ⎤⎞
                           𝐼𝑃 𝑥                                         𝑅𝑑𝑋
                𝐼𝑃 = 𝑠 ⎣𝐼𝑃 𝑦 ⎦ = 𝐼𝑎 𝑅𝑎 + 𝐼𝑖 cos 𝜃𝑖 ⎝𝑘𝑠 𝑅𝑠 + (1 − 𝑘𝑠 ) ⎣𝑅𝑑𝑌 ⎦⎠                (9)
                           𝐼𝑃 𝑧                                          𝑅𝑑𝑍

  In Equation (9), 𝐼𝑃 is the CIE tristimulus value at point P with illumination direction 𝜃𝑖
and viewing direction 𝜃𝑟 . 𝐼𝑎 𝑅𝑎 is the ambient light term which is assumed to be zero as the
experiment is performed in a dark environment. 𝐼𝑖 is the incident light intensity, 𝑅𝑑𝑋𝑌 𝑍 are
the diffuse reflectance components to be optimised. The specular components 𝑅𝑠 are given by
Equation 4 following the GGX distribution of micro-facets normals.
  An extra parameter 𝑠 is optimised as a scaling factor due to the fact the CIE Y values go well
above one. This is because the domain is limited to that of a white Lambertian surface and in
this case the samples have reflectance values above one. The scaling parameter violates the
conservation of energy constraint of the BRDF model but it is necessary to add it in order to fit
the model. To simplify the optimisation problem, the diffuse component 𝑅𝑑 is not estimated
and the value at 0/45 measurement geometry is used. Using a genetic algorithm, the RMSE
formula is used as an objective function and only the CIE Y value is estimated. The RMSE is
given by                             √︂
                                        Σ((𝑀 (𝜃𝑖 , 𝜃𝑟 ) − 𝐸(𝜃𝑖 , 𝜃𝑟 , 𝑝))2 )
                         𝑅𝑀 𝑆𝐸 =                                             ,         (10)
                                                        𝑛
where 𝑀 are the measured CIE XYZ values (Eq. 7) and 𝐸 are the estimated BRDF values
obtained using the Cook-Torrance GGX BRDF model (Eq. 9), with the parameters 𝑝, calculated
for the 𝑛 pairs of incident and reflected directions.


4. Results and Discussion
The gold foils and gilding mock-ups were measured to obtain their bi-directional reflectance
according to the geometries previously described. CIE L*a*b* values are calculated and to
evaluate the luminance of the samples, the L* values are analysed as a function of angle of
incidence. A BRDF model is estimated for the gilding mock-ups and the gold foil used for their
fabrication, namely Or Versailles.
   All the gold foils have a different appearance due to their karat composition, which gives
them a different colour. However, due to the metallic properties of gold and the fact that the gold
foils have been produced in a similar way, they should have a similar shine. On the contrary,
the gilding mock-ups have been manufactured using techniques which alter the surface of the
gold foil and thus, the shape of their BRDF is expected to be different even if their diffuse colour
is similar to that of the original gold leaf.

4.1. Gold foils: colour and reflectance spectra
The reflectance spectra of the gold foils is plotted as a function of wavelength in the visible
 spectrum range from 400 to 700 nm in Figure 5 and Figure 6. The measurement geometries
 are specified in Table 2. All the measurements have been averaged, regardless of their location
 in the gold foil and the orientation of measurement, assuming the samples are isotropic. The
 curves are smooth and for all samples at both viewing angles the reflectance factors go above
 one at illumination directions -30∘ and 0∘ (yellow and purple curves), and -60∘ and -45∘ (blue
 and orange curves) respectively. This is because the white reference used is a Lambertian
white reference which assumes a value of one for all directions and wavelengths. The different
 gold foils reflectance curves have more or less contribution on the longer end of the spectrum
 depending on their colour.
    Figure 7, shows the sRGB visualisation of the gold foils at all measurement geometries.
The rendering was done using CIE standard observer 10∘ colour matching functions and D65
 illuminant. As described previously, each gold foil has a different hue due to its different karat
 composition. The colour is rendered as white in the case of all gold foils except Or Versailles at
 angles of incidence -30∘ and 0∘ for viewing direction -15∘ , and -60∘ and -30∘ for viewing direction
-45∘ . This is because the reflectance factors exceed one for these measurement geometries and
 thus the luminance of the samples cannot be rendered in a standard dynamic range display.
                         or Alienor                                           or Alienor
     10                                                   10

      9                                                    9

      8                                                    8

      7                                                    7

      6                                                    6

      5                                                    5

      4                                                    4

      3                                                    3

      2                                                    2

      1                                                    1

      0                                                    0
      400   450    500      550       600   650   700      400   450    500      550       600   650   700



                         or Blanc                                             or Blanc
     10                                                   10

      9                                                    9

      8                                                    8

      7                                                    7

      6                                                    6

      5                                                    5

      4                                                    4

      3                                                    3

      2                                                    2

      1                                                    1

      0                                                    0
      400   450    500      550       600   650   700      400   450    500      550       600   650   700



                         or Citron                                            or Citron
     10                                                   10

      9                                                    9

      8                                                    8

      7                                                    7

      6                                                    6

      5                                                    5

      4                                                    4

      3                                                    3

      2                                                    2

      1                                                    1

      0                                                    0
      400   450    500      550       600   650   700      400   450    500      550       600   650   700



                         or Jaune                                             or Jaune
     10                                                   10

      9                                                    9

      8                                                    8

      7                                                    7

      6                                                    6

      5                                                    5

      4                                                    4

      3                                                    3

      2                                                    2

      1                                                    1

      0                                                    0
      400   450    500      550       600   650   700      400   450    500      550       600   650   700



                         or Noble                                             or Noble
     10                                                   10

      9                                                    9

      8                                                    8

      7                                                    7

      6                                                    6

      5                                                    5

      4                                                    4

      3                                                    3

      2                                                    2

      1                                                    1

      0                                                    0
      400   450    500      550       600   650   700      400   450    500      550       600   650   700




Figure 5: Reflectance curve of gold foils (1/2). The viewing angles are -15∘ and -45∘ . For each viewing
direction the illumination directions are defined in Table 2.


   Figure 8 shows the mean CIE L* value for all gold foils as a function of the illumination
direction. The error bars represent two standard deviations. For all gold foils the value of L*
greatly increases as the illumination angle approaches the mirror angle. The values of L* go
much higher than 100 because the spectrophotometer is taking the reference white to be a
Lambertian white surface.
                           or Orange                                              or Orange
     10                                                      10

      9                                                      9

      8                                                      8

      7                                                      7

      6                                                      6

      5                                                      5

      4                                                      4

      3                                                      3

      2                                                      2

      1                                                      1

      0                                                      0
      400   450    500        550          600   650   700   400    450   500         550         600   650   700



                            or Rose                                                or Rose
     10                                                      10

      9                                                      9

      8                                                      8

      7                                                      7

      6                                                      6

      5                                                      5

      4                                                      4

      3                                                      3

      2                                                      2

      1                                                      1

      0                                                      0
      400   450    500        550          600   650   700   400    450   500        550          600   650   700



                            or Rouge                                               or Rouge
     10                                                      10

      9                                                       9

      8                                                       8

      7                                                       7

      6                                                       6

      5                                                       5

      4                                                       4

      3                                                       3

      2                                                       2

      1                                                       1

      0                                                       0
      400   450    500         550         600   650   700    400   450   500         550         600   650   700



                          or Versailles                                          or Versailles
     10                                                      10

      9                                                       9

      8                                                       8

      7                                                       7

      6                                                       6

      5                                                       5

      4                                                       4

      3                                                       3

      2                                                       2

      1                                                       1

      0                                                       0
      400   450    500        550          600   650   700    400   450   500        550          600   650   700



                             or Vert                                                or Vert
     10                                                      10

      9                                                       9

      8                                                       8

      7                                                       7

      6                                                       6

      5                                                       5

      4                                                       4

      3                                                       3

      2                                                       2

      1                                                       1

      0                                                       0
      400   450    500         550         600   650   700    400   450   500         550         600   650   700
                         Wavelength (nm)                                        Wavelength (nm)




Figure 6: Reflectance curve of gold foils (2/2). The viewing angles are -15∘ and -45∘ . For each viewing
direction the illumination directions are defined in Table 2.


4.2. Gold foils: bidirectional measurements
4.3. Comparing gold foil, water gilding, and oil gilding
Different manual processes and techniques are expected to alter the surface of the gold foil
and change its appearance when used for gilding. While oil gilding is known to have a more
matte appearance, water gilding has a higher gloss due to the final burnishing step which is not
performed in the former.
  Figure 9, shows the CIE Y values as a function of the illumination direction for the oil gilding
mock-up, water gilding mock-up and Or Versailles gold foil at viewing directions -15∘ and -45∘ .
These samples have also been used to train a BRDF model using the GGX distribution and a
genetic algorithm to optimise the RMSE as objective function. The BRDF model was trained for
Figure 7: sRGB visualisation of the gold foils at all measurement geometries using CIE standard observer
10∘ colour matching functions and illuminant D65.




Figure 8: Measured CIE L* of gold foils at viewing directions -15∘ (left) and -45∘ (right).


each viewing direction rather than training a single model with both viewing directions. This
is because the data is sparsely sampled and an acceptable model that satisfied both viewing
directions was not achieved.
Figure 9: Measured CIE Y at viewing angle -15∘ (left) and -45∘ (right) respectively. BRDF model with
estimated parameters in Table 3


Table 3
BRDF model estimated parameters using GGX distribution and genetic algorithm to optimise RMSE.
                                    Water gilding     Oil gilding    Or Versailles
                                    -15∘    -45∘     -15∘     -45∘   -15∘    -45∘
                   Scaling factor   1.506   2.086   3.092    3.379   2.951   3.156
                   𝑘𝑠               0.694   0.667   0.819    0.988   0.706   0.672
                   𝛼                0.137   0.090   0.096    0.099   0.071   0.092


    The BRDF parameters estimated and used to plot the values in Figure 9 are shows in Table 3.
The scaling factor was limited to the range [0 5] and both 𝑘𝑠 and 𝛼 defined in the range [0 1].
    The three samples show a similar behaviour at both viewing directions. At viewing direction
-15∘ , oil gilding shows the higher the values of CIE Y for all angles of illumination. The gold foil
 has lower values of CIE Y for angles further from the mirror angle, but for angles of illumination
-30∘ and 0∘ , the water gilding has the lowest values of CIE Y. The same trend is seen at viewing
 direction -45∘ . However, the BRDF models show a different relative difference between the
 samples. At viewing direction -15∘ , water gilding has a much lower peak than oil gilding and
 the gold foil has the highest peak. On the other hand, for -45∘ viewing direction, oil gilding has
 the highest peak. The shape of the curves is more similar at this viewing direction as well.
    The gloss of the water gilding and oil gilding mock-ups was measured using a Rhopoint IQ
 glossmeter. Table 4 shows the gloss values obtained at 20∘ , 60∘ and 85∘ geometry. It can be seen
 that water gilding has a much higher gloss than oil gilding.

4.4. Discussion and future work
Complex materials such as gold foils and gilding present particular appearance properties
which cannot be simply described by colorimetry. Moreover, the metallic nature of the samples
requires the specular peak to be sampled with a fine resolution to properly characterise these
Table 4
Gloss measurements of water and oil gilding. The measurements were done using a Rhopoint IQ
glossmeter.
                                          Gloss 20   Gloss 60   Gloss 85
                          Water gilding     315        360         77
                          Oil gilding       25         85          22


surfaces.
   In the field of heritage science, gilded or golden surfaces are frequently characterised by
colour measurements performed at single measurement geometries. Normally, measurements
are done in a diffuse domain for example 0/45, which is insufficient to properly describe the
appearance of gold. In Figure 7, the sRGB render of the reflectances at different measurement
geometries shows that the colour of gold is very different from the visual perception that
observers have of these materials. For example at 0/-45, column ten in Figure 7, the colour of
the gold foils is very dark.
   Evidently, pure colour measurements are not sufficient to properly characterise the appearance
of golden surfaces. As demonstrated, a BRDF model certainly can help to render the appearance
of gold and gilding in more realistic way which is related to visual perception. An attempt to
estimate the BRDF model of the gilding mock ups and its gold foil was presented in Figure 9.
   The BRDF models obtained fail to properly characterise the differences in appearance of the
gold foil, water gilding, and oil gilding. This is probably due to the fact the BRDF model is a
data hungry model which requires more number of measurements to train coefficients that
can accurately represent the surface at many combinations of angles of incidence and viewing.
In this paper, the models are trained with six measurements for each viewing direction and
moreover, the reflectances are not sampled at angles close to the specular peak. The fact that
the BRDF model is trained on each individual viewing direction independently rather than just
having a single model does not allow for the models to then be directly comparable. This is
because the specular peak is unknown and thus the model cannot be properly trained.
   Finally, a remark must be made on the fact that the gilding mock-ups have been made with
the same gold foil. The only difference between the pure gold foil (Or Versailles), the oil gilding
mock-up, and the water gilding mock-up is the orientation of the micro-facets of each surface.
The Fresnel coefficient of the three surfaces will be the same since the chemical composition of
the gold foil is unchanged by either gilding techniques. The amount of radiance being absorbed
by the three surfaces is the same. Thus, the distribution of the scattered rays can be normalised
over all incident directions to a constant that will be the same for the three surfaces.
   Future work should include improvement in the fitting of the BRDF model. Although the
Cook-Torrance BRDF model using a GGX distribution seems to be the appropriate model for
the samples, some limitations have to be overcome to obtain more faithful models. First of all,
the high specularity of the surfaces must be accounted in the BRDF model so the scaling factor
is not violating any rule such as conservation of energy. Additionally, the lack of data close
to the specular peak requires that some assumptions are to be made to characterise these. For
example, the fact that the Fresnel coefficient of the surfaces is the same should be exploited.
   Moreover, once the BRDF models are obtained, they will have to be evaluated. It would
be valuable to perform a psychophysical experiment to evaluate how good the models are at
rendering the subtleties in appearance difference between these surfaces. Also, the models can
be tested for their ability to characterise the surfaces with respect to each other by performing
a classification task based on the coefficients.
   In the field of conservation science, it would be valuable to understand how the different
gilding techniques change the BRDF model and thus the appearance of the surfaces. This could
give insights into restoration techniques and procedures that render an appearance closest to
that of the original gilding. Another interesting possibility would be to characterise the effects
of ageing on gilded surfaces in order to create a model linking the original appearance of the
gilding and its aged appearance for virtual restoration purposes.


5. Conclusion
In this paper, ten different gold foils commonly used for restoration have been measured and
characterised using a bidirectional reflectance measurement device and an well established
reflectance model. Two types of gilding mock-ups, water gilding and oil gilding, fabricated with
one of the gold foils studied, have also been measured. It is found that the gold foils and gilded
mock-ups are very specular and present reflectance factors higher than one at angles as close
as 15∘ from the mirror angle.
   It is found that these measurement angles are not enough to fit a BRDF model that accurately
describes the appearance of the surfaces. It is suggested that the model is modified so the high
specularity of the surface is better described. Also, given the surfaces should have the same
transmittance and absorption, the BRDF models can be normalised to have the same value over
all angles of incidence.


Acknowledgments
This work has received funding from the European Union’s Horizon 2020 research and innova-
tion program under grant agreement No. 813789, through the the MSCA ITN project CHANGE
(https://change-itn.eu/).
   We would like to thank Steṕhanie Courtier from the Centre of Research and Restoration of
the Museums of France, for her help in the preparation of samples.


References
 [1] K. Kollandsrud, U. Plahter, Twelfth and early thirteenth century polychromy at the
     northernmost edge of europe: past analyses and future research, Medievalista online (2019).
     URL: https://dx.doi.org/10.4000/medievalista.2303. doi:10.4000/medievalista.2303.
 [2] G. P. Mastrotheodoros, K. G. Beltsios, Y. Bassiakos, V. Papadopoulou, On the metal-
     leaf decorations of post-byzantine greek icons, Archaeometry 60 (2018) 269–289. URL:
     https://onlinelibrary.wiley.com/doi/abs/10.1111/arcm.12287. doi:https://doi.org/10.
     1111/arcm.12287.
 [3] A. Sole, G. C. Guarnera, I. Farup, P. Nussbaum, Measurement and rendering of complex
     non-diffuse and goniochromatic packaging materials, The Visual Computer 37 (2021)
     2207–2220.
 [4] P. Callet, A. Zymla, A. Mofakhami, Virtual metallurgy and archaeology. application to the
     visual simulation of a work of art, in: ICCVG, 2002, pp. 25–29.
 [5] S. Dumazet, A. Genty, A. Zymla, F. De Contencin, A. Texier, N. Ruscassier, B. Bonnet,
     P. Callet, Influence of the substrate colour on the visual appearance of gilded sculptures,
     in: XXI International CIPA Symposium, 2007, pp. 01–06.
 [6] A. Mounier, F. Daniel, The role of the under-layer in the coloured perception of gildings
     in mediaeval mural paintings, Open Journal of Archaeometry 1 (2013). doi:10.4081/arc.
     2013.e16.
 [7] Q. Wu, M. Hauldenschild, B. Rösner, T. Lombardo, K. Schmidt-Ott, B. Watts, F. Nolting,
     D. Ganz, Does substrate colour affect the visual appearance of gilded medieval sculptures?
     part I: colorimetry and interferometric microscopy of gilded models, Heritage Science 8
     (2020). doi:10.1186/s40494-020-00463-3.
 [8] L. W. MacDonald, T. Vitorino, M. Picollo, R. Pillay, M. Obarzanowski, J. Sobczyk, S. Nasci-
     mento, J. Linhares, Assessment of multispectral and hyperspectral imaging systems for digi-
     tisation of a Russian icon, Heritage Science 5 (2017). doi:10.1186/s40494-017-0154-1.
 [9] M. A. Martinez, E. M. Valero, J. L. Nieves, R. Blanc, E. Manzano, J. L. Vilchez, Multifocus
     hdr vis/nir hyperspectral imaging and its application to works of art, Opt Express 27 (2019)
     11323–11338. doi:10.1364/OE.27.011323.
[10] O. Sidorov, J. Y. Hardeberg, S. George, J. S. Harvey, H. E. Smithson, Changes in the visual
     appearance of polychrome wood caused by (accelerated) aging, in: Electronic Imaging,
     Society for Imaging Science and Technology, 2020, pp. 060–1–060–7. doi:doi.org/10.
     2352/ISSN.2470-1173.2020.5.MAAP-060.
[11] F. E. Nicodemus, R. JC, et al., Geometrical considerations and nomenclature for reflectance.
     (1977).
[12] D. Guarnera, G. C. Guarnera, A. Ghosh, C. Denk, M. Glencross, Brdf representation and
     acquisition, in: Computer Graphics Forum, volume 35, 2016, pp. 625–650.
[13] R. L. Cook, K. E. Torrance, A reflectance model for computer graphics, ACM Transactions
     on Graphics 1 (1982) 7–24.
[14] A. Sole, I. Farup, P. Nussbaum, S. Tominaga, Bidirectional reflectance measurement
     and reflection model fitting of complex materials using an image-based measurement
     setup, Journal of Imaging 4 (2018). URL: https://www.mdpi.com/2313-433X/4/11/136.
     doi:10.3390/jimaging4110136.
[15] K. E. Torrance, E. M. Sparrow, Theory for off-specular reflection from roughened surfaces,
     Journal of the Optical Society of America 57 (1967) 1105–1112.
[16] B. Walter, S. R. Marschner, H. Li, K. E. Torrance, Microfacet models for refraction through
     rough surfaces, in: Eurographics conference on Rendering Techniques, Eurographics
     Association, 2007, pp. 195–206.
[17] A. Fores, J. Ferwerda, J. Gu, Toward a perceptually based metric for brdf modeling, in:
     Color and imaging conference, Society for Imaging Science and Technology, 2012, pp.
     142–148.
[18] Standard Practice for Specifying the Geometry of Multiangle Spectrophotometers, ASTM
     E2175-01, Standard, American Society for Testing and Materials, Pennsylvania, USA, 2001.
[19] Standard Practice for Multiangle Color Measurement of Metal Flake Pigmented Materials,
     ASTM E2194-12, Standard, American Society for Testing and Materials, Pennsylvania,
     USA, 2012.
[20] Tolerances for Automotive Paint—Part 2: Goniochromatic Paints, DIN 6175-2, Standard,
     Deutsches Institut für Normung, e.V, Berlin, DE, 2001.



6. Online Resources
The dataset of gold foils reflectances is available in the CHANGE-ITN repository (https://zenodo.
org/communities/change/?page=1&size=20).