<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta>
      <journal-title-group>
        <journal-title>Asia-PacificWorkshop on Mixed and Aug-
mented Reality, Dec.</journal-title>
      </journal-title-group>
    </journal-meta>
    <article-meta>
      <title-group>
        <article-title>Manipulation of anisotropic reflections based on optical models using multiple projectors</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Shogo Ohsumi</string-name>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Toshiyuki Amano</string-name>
        </contrib>
      </contrib-group>
      <pub-date>
        <year>2022</year>
      </pub-date>
      <volume>0</volume>
      <fpage>2</fpage>
      <lpage>03</lpage>
      <abstract>
        <p>This paper proposes a novel appearance-manipulation technique that parametrically manipulates the visible anisotropic reflection property with illumination projection from multiple projectors. This method obtains a reflectance matrix corresponding to the bidirectional reflectance distribution function (BRDF) from images captured using multiple cameras. The reflectance matrix was then fitted to the Ashikhmin BRDF model to estimate its parameters of the BRDF model. The reflectance matrix corresponding to the target BRDF was then calculated by manipulating the estimated parameters. The anisotropic reflection was manipulated based on the optical model by projecting images from multiple projectors that changed the texture represented by the reflectance matrix calculated in this manner.</p>
      </abstract>
      <kwd-group>
        <kwd>1 Anisotropy</kwd>
        <kwd>Light-field projection</kwd>
        <kwd>BRDF</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>
        The angular light intensity distribution on the
surface is formed by its properties (e.g.,
bidirectional reflectance distribution function (BRDF),
bidirectional transmittance distribution function
(BTDF)) and represents rich materiality, such as
glossy metallic reflection, clear glass caustic, and
beautiful structural color. Meanwhile, precisely
designed light-field projection, instead of normal
environmental illumination, has the potential to
manipulate light angular distribution and alter our
perception of materiality [
        <xref ref-type="bibr" rid="ref1 ref2 ref3 ref4">1-4</xref>
        ]. Such material
appearance manipulation is a key challenge in
spatial augmented reality (SAR), known as projection
mapping. This paper proposes an anisotropic
reflection property manipulation, which is an
angular distribution manipulation of the reflected light
ray on an anisotropic reflection surface using light
field projection as a novel SAR technique.
      </p>
    </sec>
    <sec id="sec-2">
      <title>2. Related work</title>
    </sec>
    <sec id="sec-3">
      <title>2.1. Auto-stereoscopic display</title>
      <p>
        Horizontally aligned projectors and a screen
composed of lenticular lenses with a diffusing
screen can achieve a projection-based
autostereoscopic display [
        <xref ref-type="bibr" rid="ref5 ref6">5,6</xref>
        ]. Jones et al. [
        <xref ref-type="bibr" rid="ref7">7</xref>
        ] demonstrated
a wide-viewing and high-angular-resolution
autostereoscopic 3D display using 216 projectors.
Nagano et al. [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ] proposed an autostereoscopic
projection display with 72 overlay images projected
onto a vertically oriented lenticular screen with
black back. Such front-projection
auto-stereoscopic displays can be used to show complex
materiality on an object using retroreflection paint on
a 3D object [
        <xref ref-type="bibr" rid="ref9">9</xref>
        ]. However, they only displayed a
BRDF and did not realize the alternation or
manipulation of the BRDF that the object originally
had.
2.2.
      </p>
    </sec>
    <sec id="sec-4">
      <title>VDDAM</title>
      <p>
        Amano et al. [
        <xref ref-type="bibr" rid="ref10">10</xref>
        ] demonstrated view-direction
dependent appearance manipulation (VDDAM),
using multiple projector-camera feedback
systems. Murakami et al. [
        <xref ref-type="bibr" rid="ref11">11</xref>
        ] proposed another
method for VDDAM based on reflectance
measurement among multiple projectors and cameras,
which was equivalent to roughly sampled BRDF.
Amano and Yoshioka [
        <xref ref-type="bibr" rid="ref12">12</xref>
        ] combined reflectance
analysis with multiple projector-camera feedback
and expanded the applicable reflection property to
retroreflection and improved robustness against
environmental lighting changes. However, the
parameters for anisotropy are unknown for Amano
et al. [
        <xref ref-type="bibr" rid="ref10">10</xref>
        ], whereas Murakami et al. [
        <xref ref-type="bibr" rid="ref11">11</xref>
        ] and
Amano and Yoshioka [
        <xref ref-type="bibr" rid="ref12">12</xref>
        ] must consider in
advance what type of anisotropy was used to create
the target image.
      </p>
      <p>In this paper, we propose a method to
parameterize the VDDAM applying the Ashikhmin
BRDF model to the previously acquired reflection
characteristics. With this parameterization, we
enhanced or reduced the anisotropy and then
recreated the reflectance. Finally, we calculated the
manipulation references for each viewing
direction using the recreated reflectance.
2.3.</p>
    </sec>
    <sec id="sec-5">
      <title>Reflectance matrix</title>
      <p>
        Murakami and Amano proposed a response
model for multiple projectors and cameras that
considers color [
        <xref ref-type="bibr" rid="ref11">11</xref>
        ]. The RGB values at a point A
in the captured and projected images are defined
as follows:
"! = $!", !#, !$'%, !" ≥ 0, !# ≥ 0, !$ ≥ 0,
ℎ  = 1,2, … , ,
7&amp; = 8&amp;", &amp;#, &amp;$:%, &amp;" ≥ 0, &amp;# ≥ 0, &amp;$ ≥ 0,
ℎ  = 1,2, … , ,
where  denotes the number of cameras, and 
denotes the number of projectors. In this case, by
expressing the reflection at an object surface as a
matrix 7 ∈ ℛ'×', it can be described as follows:
"! = 7!&amp;!&amp;7&amp; ,
where the matrix !&amp; represents the color-mixing
matrix [
        <xref ref-type="bibr" rid="ref13">13</xref>
        ] that calibrates the color. Furthermore,
when multiple projectors or cameras are used,
they are represented as follows:
      </p>
      <p>B = CCC ,
where
" = $%!", %#", … , %$"(" , + = $,!", ,#", … , ,$"(" , (5)
!!
* = , ⋮"!
!" ⋯ !# 6!! 6!" ⋯ 6!#
"" ⋯ ⋮"#0 , * = ⎛6"! 6"" ⋯ 6"#⎞ . (6)</p>
      <p>⋮ ⋱ ⋮ ⋮ ⋱ ⋮
$! $" ⋯ $# ⎝6$! 6$" ⋯ 6$#⎠
Hereafter, we regard the color spaces as calibrated,
and we write CC as C in the following sections.
2.4.</p>
    </sec>
    <sec id="sec-6">
      <title>Ashikhmin BRDF model[14]</title>
      <p>(1)
(2)
(3)
(4)</p>
      <p>In this section, we introduce the Ashikhmin
BRDF model to fit the reflectance matrix. The
Ashikhmin BRDF model is described by the sum
of the specular and diffuse components, with the
specular component = defined as follows:
!(", #) ('!()"+'#()")
= ( ($ + 1)(% + 1) () ("-()")</p>
      <p>8 () 3("), (#)4 5! + (1 − !) 31 − ()4/8, (7)
where &gt; and ? represent normalized vector to
the light and viewer.  represent normalized
halfvector between &gt; and ? .  represent surface
normal to macroscopic surface. @ and A
represent two phong-like exponents that control the
specular lobe shape. The larger the value of @,
the higher the directivity of reflection in the u
direction. Similarly, the larger the value of A, the
higher the directivity of reflection in the v
direction.</p>
    </sec>
    <sec id="sec-7">
      <title>3. Proposed method</title>
      <p>Our proposed method obtained a reflectance
matrix representing the optical response between
projectors and cameras corresponding to a
roughly sampled BRDF on every single point on
the object's surface with the experimental devices
shown in Figure 1. Subsequently, we fitted the
reflectance matrix with the Ashikhmin BRDF
model and parameterized the reflectance
relationship. Then, the parameters were manipulated to
design a desired anisotropic reflection, yielding a
recreated reflection matrix. Finally, the VDDAM
based on the optical model achieved the desired
appearance that the reflectance matrix represented
by projecting images from multiple projectors.</p>
    </sec>
    <sec id="sec-8">
      <title>Multiple projector-camera sys</title>
      <p>In this paper, we employed 7 cameras (Ximea，
MQ013CGE2, resolution: 1280 × 1024 ) and 7
projectors (EPSON,EB-W05, resolution:1280 ×
800) in order to achieve high quality perceptual
BRDF manipulation with complex reflection
characteristics. The cameras and projectors were
placed in front of the target object and the other
projectors are placed radially at 15. intervals
around the projector 4 (Figure 1). The cameras are
placed close to each projector. This arrangement
takes into account the measurement and
manipulation of anisotropic reflections.</p>
      <p>To obtain the reflectance matrix, first capture
the red projection from one projector with all
cameras. Similarly, the green and blue projection from
one projector is captured by all cameras. This
process is repeated with seven projectors. The
reflectance is obtained by dividing the image thus
acquired by the RGB of the projection image.
3.2.</p>
    </sec>
    <sec id="sec-9">
      <title>Manipulation Target Object</title>
      <p>We used a drawing foil of Nishijin silk textile,
which contains patterns of birds, flowers, clouds,
and a mountain with rivers, as the manipulation
target. Various threads, including gold thread and
dyed thread, are used in this Nishijin silk textile,
and differences in gloss can be seen, such as the
gold thread being more reflective than the dyed
thread. The weaving method also causes
differences in reflectance characteristics. In the case of
twill weave, the ratio of warp to weft threads on
the surface is close, resulting in isotropic
reflections. On the other hand, a satin weave has a
higher ratio of warp threads than weft threads on
the surface, resulting in anisotropic reflections. In
this study, we regard the target object as a plane.
3.3.</p>
    </sec>
    <sec id="sec-10">
      <title>Parameter estimation</title>
      <p>We employed the Levenberg-Marquardt
method, a nonlinear optimization scheme, and
obtained anisotropic reflection parameters by
minimizing the error function as follows:
(!"#, $"#, %"#, ""#)
, *
where  ∈ (, , ) ,  ∈ (, , ) , and &amp;"'#
represents the reflectance of the  color component of
the  color projection of projector  captured by
camera . @=B, A=B represent anisotropic scattering
for each direction, and u, v, C=B, ==B represent the
intensities of the diffuse and specular components,
respectively. Because the four parameters @, A,
=, and C must be positive, we applied a
nonnegative condition to the Levenberg-Marquardt
method. This allows us to obtain the four
parameters of the Ashikhmin model for a single pixel
from the reflectance matrix at a single pixel.</p>
      <p>Figure 2 shows the estimated @ and A.
Because there is no significant difference the color
channel, this figure shows only the R channel.
Brightness expresses the value of each parameter,
and the brighter area has a sharp specular
reflection along each direction. A small difference
between the @ and A values indicates isotropy,
whereas a large difference indicates anisotropy.
Area (a) in the figure shows a twill weave using
gold threads and has almost isotropic reflections.
Area (b) is a satin weave that uses gold threads
and exhibits strong anisotropy. Area (c) is a satin
weave using dyed threads, and it has weak
anisotropy. Area (d) is a twill weave using dyed threads
and exhibits diffuse reflection.
3.4.</p>
    </sec>
    <sec id="sec-11">
      <title>Anisotropy manipulation</title>
      <p>Our anisotropic manipulation aims to enhance
or reduce its reflection of the anisotropic
reflection optically while maintaining the glossiness of
the isotropic reflection. Based on this, we updated
the parameters as follows:
:; = : + (: − &lt;), &lt;; = &lt; (: ≥ &lt;) . (9)
= )&amp;+( )'+((&amp;"'# − ((', )&amp;; !"#, $"#, %"#, ""#))) , (8) !:; = :, &lt;; = &lt; + (&lt; − :) (: &lt; &lt;)
(a) Anisotropic reduction (b) Anisotropic enhancement</p>
      <p>Figure 3:Target images.</p>
      <p>(a) Anisotropic reduction (b) Anisotropic enhancement</p>
      <p>Figure 4: Projection images.
(a) Anisotropic reduction
(c) Anisotropic enhancement
(b) Original appearance</p>
      <p>Figure 5: Projection results.</p>
      <p>If @was greater than A, the difference, which
is an anisotropy, was added to @with a
multiplication of the scaling factor α. Otherwise, the
difference was added to A as well. When  &gt; 1, the
difference was expanded, and anisotropy was
enhanced. On the contrary, the anisotropy was
completely removed when  = −1.</p>
      <p>We applied this manipulation to the estimated
parameters and obtained the desired target images
for the entire manipulation area.
3.5.</p>
    </sec>
    <sec id="sec-12">
      <title>Calculation of projection image</title>
      <p>From the reflectance matrix and the target
images BB = (BB%&gt;, BB%?, … , BB%@)% , we obtained the
projection images CB = (CB%&gt;, CB%?, … , CB%@)% for
each projector. However, the projection images
B must be positive. Therefore, in this paper, a
C
non-negative conditional optimization problem
min ∥ CCB − BB ∥?,</p>
      <p>D</p>
      <p>
        #
ℎ B"&gt; ≥ 0, B&gt; ≥ 0, B$&gt; ≥ 0, … , B$A ≥ 0. (10)
used the Lawson-Hanson algorithm [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ] to obtain
the projection value of each projector at a certain
point. This calculation was performed for all the
points in the operating range to obtain the
projection image for each projector.
(a) Satin weave
(b) Twill weave
      </p>
    </sec>
    <sec id="sec-13">
      <title>4. Result</title>
      <p>The target images were created by updating @
and A, and by manipulating the value of  in Eq.
(9), the projection images were obtained using Eq.
(10), and then projected.</p>
      <p>Figures 3 and 4 shows examples of the target
and projection images, respectively. In Figure 5,
the projection results are arranged to correspond
to each camera position. All of these images are
shown in identical aspects by geometrical
transformation to the common coordinate (cam 4).</p>
    </sec>
    <sec id="sec-14">
      <title>Anisotropic enhancement</title>
      <p>Figure 3 shows the target images created using
 = 8. We solved the non-negative optimization
problem described in Section 3.3 and obtained the
projection images shown in Figure 4(b). The
manipulation results from the projection are shown
in Figure 5(c). It should be noted that a significant
difference in glossiness between viewing
directions along u and v was observed. Figure 6 shows
the brightness change from cam1 to cam5 and
cam3 to cam7 when the viewpoint is moved from
left to right. Figure 6(a) shows the average
brightness of the 3 × 3 pixels in the area shown in
Figure 2(b). When the viewpoint is moved
horizontally in the range of cam1 to cam5 and cam3 to
cam7, the gloss change of the anisotropy
enhanced image is sharper than that of the original
appearance. These results confirmed the enhanced
anisotropy.</p>
      <p>Figure 6(b) shows the average brightness
values of the 3 × 3 pixels in the area shown in Figure
2(a). When the viewpoint is moved horizontally
in the range of cam1 to cam5 and cam3 to cam7,
it can be confirmed that the anisotropic
enhancement follows the brightness change of the original
appearance, although the brightness of the
anisotropic enhancement is reduced compared to the
original appearance. This confirms that the
isotropic reflection was maintained in the region of
the gold thread twill weave.
4.2.</p>
    </sec>
    <sec id="sec-15">
      <title>Anisotropic reduction</title>
      <p>even with anisotropic reduction in the region of
the gold thread twill weave.</p>
    </sec>
    <sec id="sec-16">
      <title>5. Conclusion</title>
      <p>In this paper, we propose a method for
estimating the parameters of the Ashikhmin model from
the reflectance matrix. Furthermore, we propose a
parameter manipulation method that can enhance
and reduce anisotropy. Projection images were
calculated using non-negative conditional
optimization. The projection results showed that
anisotropy could be enhanced and reduced in areas with
anisotropic reflections. In areas with isotropic
reflections, although the brightness was slightly
reduced, the isotropic reflections were maintained.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>S.</given-names>
            <surname>Shimazu</surname>
          </string-name>
          et al.,
          <article-title>“3d high dynamic range display system”</article-title>
          ,
          <source>ISMAR</source>
          , pp.
          <fpage>235</fpage>
          -
          <lpage>236</lpage>
          , (
          <year>2011</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>T.</given-names>
            <surname>Amano</surname>
          </string-name>
          , H. Kato, “
          <article-title>Appearance control using projection with model predictive control”</article-title>
          ,
          <source>ICPR</source>
          , pp.
          <fpage>2832</fpage>
          -
          <lpage>2835</lpage>
          , (
          <year>2010</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>T.</given-names>
            <surname>Okazaki</surname>
          </string-name>
          et al.,
          <article-title>“A projector-camera system for high-quality synthesis of virtual reflectance on real object surfaces”</article-title>
          , IPSJ, pp.
          <fpage>71</fpage>
          -
          <lpage>83</lpage>
          , (
          <year>2010</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>A. J.</given-names>
            <surname>Law</surname>
          </string-name>
          et al.,
          <article-title>“Perceptually based appearance modification for compliant appearance editing”</article-title>
          ,
          <source>CG Forum</source>
          ,
          <volume>30</volume>
          (
          <issue>8</issue>
          ), pp.
          <fpage>2288</fpage>
          -
          <lpage>2300</lpage>
          , (
          <year>2011</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>W.</given-names>
            <surname>Matusik</surname>
          </string-name>
          , H. Pfister, “
          <article-title>3D TV: a scalable system for real-time acquisition, transmission, and autostereoscopic display of dynamic scenes”</article-title>
          ,
          <source>SIGGRAPH</source>
          ,
          <volume>23</volume>
          (
          <issue>3</issue>
          ), pp.
          <fpage>814</fpage>
          -
          <lpage>824</lpage>
          , (
          <year>2004</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>R.</given-names>
            <surname>Yang</surname>
          </string-name>
          et al.,
          <article-title>“Toward the light field display: Autostereoscopic rendering via a cluster of projectors”</article-title>
          ,
          <source>IEEE TVCG</source>
          ,
          <volume>14</volume>
          (
          <issue>1</issue>
          ), pp.
          <fpage>84</fpage>
          -
          <lpage>96</lpage>
          , (
          <year>2008</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>A.</given-names>
            <surname>Jones</surname>
          </string-name>
          et al.,
          <article-title>“An automultiscopic projector array for interactive digital humans</article-title>
          ” SIGGRAPH, (
          <year>2015</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [8]
          <string-name>
            <given-names>K.</given-names>
            <surname>Nagano</surname>
          </string-name>
          et al.,
          <article-title>“An autostereoscopic projector array optimized for 3D facial display”</article-title>
          , SIGGRAPH, (
          <year>2013</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          [9]
          <string-name>
            <given-names>T.</given-names>
            <surname>Amano</surname>
          </string-name>
          , K. Minami, “
          <article-title>Structural color display on retroreflective objects”</article-title>
          ,
          <source>ICAT - EGVE</source>
          , pp.
          <fpage>37</fpage>
          -
          <lpage>44</lpage>
          , (
          <year>2015</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          [10]
          <string-name>
            <given-names>T.</given-names>
            <surname>Amano</surname>
          </string-name>
          et al., “
          <article-title>Viewpoint-Dependent Appearance-Manipulation with Multiple Projector-Camera Systems"</article-title>
          , ICAT-EGVE, pp.
          <fpage>101</fpage>
          -
          <lpage>107</lpage>
          , (
          <year>2017</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          [11]
          <string-name>
            <given-names>K.</given-names>
            <surname>Murakami</surname>
          </string-name>
          , T. Amano, “
          <article-title>Materiality Manipulation by Light-Field Projection from Reflectance Analysis”</article-title>
          ,
          <source>ICAT-EGVE</source>
          , pp.
          <fpage>99</fpage>
          -
          <lpage>105</lpage>
          , (
          <year>2018</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          [12]
          <string-name>
            <given-names>T.</given-names>
            <surname>Amano</surname>
          </string-name>
          , H. Yoshioka, “
          <article-title>Viewing-Direction Dependent Appearance Manipulation Based on Light-Field Feedback”</article-title>
          ,
          <source>EuroVR</source>
          , pp.
          <fpage>192</fpage>
          -
          <lpage>205</lpage>
          , (
          <year>2020</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          [13]
          <string-name>
            <given-names>S. K.</given-names>
            <surname>Nayar</surname>
          </string-name>
          et al. “
          <article-title>A projection system with radiometric compensation for screen imperfections”</article-title>
          , PROCAMS, (
          <year>2003</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          [14]
          <string-name>
            <given-names>M.</given-names>
            <surname>Ashikhmin</surname>
          </string-name>
          , P. Shirley, “
          <article-title>An anisotropic phong BRDF model”</article-title>
          ,
          <source>Journal of Graphics Tools</source>
          , (
          <year>2000</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref15">
        <mixed-citation>
          [15]
          <string-name>
            <given-names>C. L.</given-names>
            <surname>Lawson</surname>
          </string-name>
          ,
          <string-name>
            <given-names>R. J.</given-names>
            <surname>Hanson</surname>
          </string-name>
          , “Solving least squares problems”, (
          <year>1974</year>
          ).
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>