Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7633790B2 - Texture acquisition system, texture acquisition device, and texture acquisition program - Google Patents
[go: Go Back, main page]

JP7633790B2 - Texture acquisition system, texture acquisition device, and texture acquisition program - Google Patents

Texture acquisition system, texture acquisition device, and texture acquisition program Download PDF

Info

Publication number
JP7633790B2
JP7633790B2 JP2020168205A JP2020168205A JP7633790B2 JP 7633790 B2 JP7633790 B2 JP 7633790B2 JP 2020168205 A JP2020168205 A JP 2020168205A JP 2020168205 A JP2020168205 A JP 2020168205A JP 7633790 B2 JP7633790 B2 JP 7633790B2
Authority
JP
Japan
Prior art keywords
subject
polarized
reflection coefficient
information
texture
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2020168205A
Other languages
Japanese (ja)
Other versions
JP2022060641A (en
Inventor
秀樹 三ツ峰
淳 洗井
俊枝 三須
寛史 盛岡
敦志 荒井
太一 杉之下
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Japan Broadcasting Corp
Original Assignee
Japan Broadcasting Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Japan Broadcasting Corp filed Critical Japan Broadcasting Corp
Priority to JP2020168205A priority Critical patent/JP7633790B2/en
Publication of JP2022060641A publication Critical patent/JP2022060641A/en
Application granted granted Critical
Publication of JP7633790B2 publication Critical patent/JP7633790B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Image Generation (AREA)
  • Studio Devices (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Description

本発明は、コンピュータグラフィックス(CG)において扱う素材データである被写体情報として、形状や質感などを取得する、質感取得システム、質感取得装置および質感取得プログラムに関する。 The present invention relates to a texture acquisition system, texture acquisition device, and texture acquisition program that acquire shape, texture, and other information about a subject, which is material data used in computer graphics (CG).

近年のAR(Augmented Reality)/VR(Virtual Reality)技術の普及により、実世界の被写体をCG技術などで扱える被写体情報(以下、「3次元モデル」と称する場合がある。)に変換するニーズが高まっている。これは一般的にヴォリュメトリックキャプチャと呼ばれる技術で、被写体を取り囲むように配置されたカメラや奥行き計測装置により被写体全体の表面模様や形状を取得する(非特許文献1,2参照)。 With the recent spread of AR (Augmented Reality)/VR (Virtual Reality) technology, there is an increasing need to convert real-world subjects into subject information (hereinafter sometimes referred to as "3D models") that can be handled using CG technology. This is a technology generally known as volumetric capture, which captures the surface pattern and shape of the entire subject using cameras and depth measurement devices arranged to surround the subject (see Non-Patent Documents 1 and 2).

3次元モデルをCG技術において希望する照明条件で描画する場合、従来のヴォリュメトリックキャプチャでは、フラットな照明環境下に被写体を置くことで、被写体表面に一様に照明を当て、陰影なく撮影している。これは被写体の反射率を直接計測できないことから、疑似的に照明が方向性を有しないような照明効果として撮影することにより、被写体の各部位への入射光量を正規化し、得られる撮影映像上の輝度を疑似的に拡散反射係数のように扱う。これにより、レンダリング時にシェーディング(陰影)を付与することで仮想空間における照明条件で見込まれるような3次元モデルの見え方を表現することができる。 When rendering a 3D model under the desired lighting conditions using CG technology, conventional volumetric capture involves placing the subject in a flat lighting environment, lighting the subject's surface evenly, and photographing the subject without shadows. Since the reflectance of the subject cannot be measured directly, the lighting is photographed as if it were a pseudo-non-directional lighting effect, normalizing the amount of light incident on each part of the subject, and treating the brightness in the resulting photographed image as if it were a pseudo diffuse reflectance coefficient. This makes it possible to express the appearance of the 3D model as expected under the lighting conditions in the virtual space by adding shading during rendering.

O. Schreer, et al., “Capture and 3D Video Processing of Volumetric Video,”2019 IEEE International Conference on Image Processing (ICIP), Taipei, Taiwan, 2019, pp. 4310-4314.O. Schreer, et al., “Capture and 3D Video Processing of Volumetric Video,” 2019 IEEE International Conference on Image Processing (ICIP), Taipei, Taiwan, 2019, pp. 4310-4314. O. Schreer, et al.,“ADVANCED VOLUMETRIC CAPTURE AND PROCESSING,” [online],IBC2018,[令和2年7月20日検索],インターネット<URL:https://www.ibc.org/download?ac=6559>O. Schreer, et al., “ADVANCED VOLUMETRIC CAPTURE AND PROCESSING,” [online], IBC2018, [Retrieved July 20, 2020], Internet <URL: https://www.ibc.org/download?ac=6559>

しかしながら、従来のヴォリュメトリックキャプチャ技術を用いた手法では、反射光の拡散反射成分を疑似的に求めている。また、鏡面反射成分については、扱われていないか、手動により被写体全体で一様あるいはエリアごとに係数(鏡面反射係数)を設定しており、非剛体が対象である場合には、エリアごとに手動で係数を付与するのは困難である。このことから、特にライブ映像に対応させることは難しかった。
また、被写体表面を微視的に観察すると細かな凹凸の有無などで、鏡面反射光の振る舞いが異なる。この被写体表面の凹凸の有無は、被写体を見た際のつやや手触り感などの質感について、見た目の印象に大きく影響を与える要素である。そのため、この質感を出すためには、エリア単位以上に細かい単位で表面反射の振る舞いを記述する必要があるが、現状では表現できない。
However, in conventional volumetric capture techniques, the diffuse reflection component of reflected light is calculated artificially. In addition, the specular reflection component is not handled, or the coefficient (specular reflection coefficient) is set manually either uniformly for the entire subject or for each area. When the subject is a non-rigid body, it is difficult to manually assign coefficients for each area. For this reason, it is particularly difficult to apply it to live video.
Furthermore, when observing the surface of a subject microscopically, the behavior of specularly reflected light differs depending on the presence or absence of fine irregularities. The presence or absence of irregularities on the subject's surface is a factor that greatly affects the visual impression of the subject when looking at it, such as its gloss and feel to the touch. Therefore, in order to express this texture, it is necessary to describe the behavior of surface reflection in units finer than area units, but this cannot be expressed at present.

本発明は、以上のような点を鑑みてなされたものであり、被写体の質感を表現するための被写体表面の細かな凹凸の状態を、3次元モデルとしてリアルタイムに取得することができる、質感取得システム、質感取得装置および質感取得プログラムを提供することを課題とする。 The present invention has been made in consideration of the above points, and aims to provide a texture acquisition system, texture acquisition device, and texture acquisition program that can acquire the fine unevenness of a subject's surface in real time as a three-dimensional model in order to express the subject's texture.

前記課題を解決するために、本発明の質感取得システムは、被写体を取り囲んで配置される、奥行き計測装置、偏光照明装置および偏光撮影装置を有する撮影ブースと、奥行き計測装置、偏光照明装置および偏光撮影装置に通信接続され、被写体表面の凹凸の状態で示される質感を取得する質感取得装置と、を備える構成とした。
また、質感取得装置は、記憶部と、被写体形状取得部と、凹凸状態算出部と、反射係数算出部と、表面形状取得部と、3次元モデル詳細情報出力部とを備える構成とした。
In order to solve the above problems, the texture acquisition system of the present invention is configured to include a photography booth having a depth measurement device, a polarized lighting device, and a polarized photography device arranged to surround a subject, and a texture acquisition device that is communicatively connected to the depth measurement device, the polarized lighting device, and the polarized photography device, and acquires the texture represented by the uneven state of the subject's surface.
The texture acquisition device includes a storage unit, a subject shape acquisition unit, a concave-convex state calculation unit, a reflection coefficient calculation unit, a surface shape acquisition unit, and a three-dimensional model detailed information output unit.

かかる構成によれば、質感取得装置は、被写体形状取得部によって、奥行き計測装置が計測した奥行きの情報を受信することにより、被写体の形状情報を取得し、計測制御部によって、複数色の偏向光を偏光照明装置が被写体に照射する制御を行い、凹凸状態算出部が、複数色の偏向光が照射された被写体からの反射光を偏光撮影装置が撮影した情報を取得し、被写体表面の凹凸により生じた偏向光の色ベクトルの分散から被写体表面の粗さを算出することができる。
これにより、質感取得システムの質感取得装置は、被写体表面の凹凸により生じた偏向光の色ベクトルの分散に基づき、被写体表面の粗さを算出することが可能となる。
With this configuration, the texture acquisition device acquires shape information of the subject by receiving depth information measured by the depth measurement device through the subject shape acquisition unit, and the measurement control unit controls the polarized lighting device to irradiate the subject with polarized light of multiple colors.The unevenness state calculation unit acquires information on the reflected light from the subject irradiated with polarized light of multiple colors captured by the polarized photography device, and can calculate the roughness of the subject surface from the variance of the color vectors of the polarized light generated by the unevenness of the subject surface.
This enables the texture acquisition device of the texture acquisition system to calculate the roughness of the object surface based on the variance of the color vector of the deflected light caused by the unevenness of the object surface.

また、質感取得装置は、反射係数算出部によって、反射光を拡散反射成分と鏡面反射成分とに分離し、拡散反射成分、被写体への入射光量、被写体形状取得部が取得した被写体の形状情報および照明位置を用いて、被写体表面の拡散反射係数を算出するともに、鏡面反射成分、被写体への入射光量、凹凸状態算出部が算出した被写体表面の粗さおよびカメラ位置を用いて、被写体表面の鏡面反射係数を算出することができる。
これにより、質感取得システムの質感取得装置は、反射光を拡散反射成分と鏡面反射成分に分離した上で、被写体表面の拡散反射係数および鏡面反射係数を算出することができる。
In addition, the texture acquisition device separates reflected light into a diffuse reflection component and a specular reflection component using the reflection coefficient calculation unit, and calculates the diffuse reflection coefficient of the subject surface using the diffuse reflection component, the amount of light incident on the subject, and the shape information of the subject acquired by the subject shape acquisition unit and the lighting position, and also calculates the specular reflection coefficient of the subject surface using the specular reflection component, the amount of light incident on the subject, the roughness of the subject surface calculated by the unevenness calculation unit, and the camera position.
This allows the texture acquisition device of the texture acquisition system to separate the reflected light into a diffuse reflection component and a specular reflection component, and then calculate the diffuse reflection coefficient and the specular reflection coefficient of the subject surface.

また、質感取得装置は、表面形状取得部によって、計測制御部を介して複数の異なる光源方向から複数色の偏向光を被写体に照射させ、偏光撮影装置が撮影した被写体表面の拡散反射成分についての輝度に基づき法線を推定することにより被写体表面の形状情報を被写体表面形状情報として取得することができる。
これにより、質感取得システムの質感取得装置は、鏡面反射成分を除いた拡散反射成分のみにより被写体表面の法線を推定できるため、被写体表面の詳細な形状情報をより正確に取得することができる。
In addition, the texture acquisition device can acquire shape information of the subject's surface as subject's surface shape information by irradiating the subject with polarized light of multiple colors from multiple different light source directions via the measurement control unit, and estimating normals based on the brightness of the diffuse reflection components of the subject's surface photographed by the polarization imaging device.
As a result, the texture acquisition device of the texture acquisition system can estimate the normal of the subject's surface using only the diffuse reflection components excluding the specular reflection components, thereby making it possible to more accurately acquire detailed shape information of the subject's surface.

また、質感取得装置は、3次元モデル詳細情報出力部によって、被写体表面の粗さ、拡散反射係数、鏡面反射係数、および、被写体表面形状情報を含む情報を、その被写体表面の粗さ、拡散反射係数、鏡面反射係数、および、被写体表面形状情報を含む情報を算出するために偏光撮影装置が撮影した時刻に合わせて3次元モデル詳細情報として出力することができる。
これにより、質感取得システムの質感取得装置は、3次元モデルのより詳細な情報を、例えばCG描画装置等にリアルタイムに出力することができる。
In addition, the texture acquisition device can output information including the subject's surface roughness, diffuse reflection coefficient, specular reflection coefficient, and subject's surface shape information as three-dimensional model detailed information via the three-dimensional model detailed information output unit in accordance with the time when the polarized imaging device captured the image in order to calculate the information including the subject's surface roughness, diffuse reflection coefficient, specular reflection coefficient, and subject's surface shape information .
This allows the texture acquisition device of the texture acquisition system to output more detailed information of the three-dimensional model to, for example, a CG drawing device in real time.

このように、本発明の質感取得システムによれば、被写体の質感を表現するための被写体表面の細かな凹凸の状態について、従来のヴォリュメトリックキャプチャ技術を用いた手法では取得できない、正確な被写体表面の拡散反射係数や鏡面反射係数、被写体表面の凹凸の状態(粗さ)、および、被写体表面の詳細な形状情報(被写体表面形状情報)をリアルタイムに取得することが可能となる。 In this way, the texture acquisition system of the present invention makes it possible to acquire in real time accurate diffuse reflection coefficients and specular reflection coefficients of the subject's surface, the state of the subject's surface unevenness (roughness), and detailed shape information of the subject's surface (subject's surface shape information) in order to express the subject's texture, which cannot be acquired using conventional volumetric capture techniques.

本発明によれば、被写体の質感を表現するための被写体表面の細かな凹凸の状態を、3次元モデルとしてリアルタイムに取得する、質感取得システム、質感取得装置および質感取得プログラムを提供することができる。 The present invention provides a texture acquisition system, texture acquisition device, and texture acquisition program that can acquire the fine unevenness of a subject's surface in real time as a three-dimensional model to express the subject's texture.

本実施形態に係る質感取得装置を含む質感取得システムの全体構成を示す機能ブロック図である。1 is a functional block diagram showing an overall configuration of a texture acquisition system including a texture acquisition device according to an embodiment of the present invention. 本実施形態における反射光の拡散反射成分と鏡面反射成分の分離に関する説明図である。4 is an explanatory diagram relating to separation of a diffuse reflection component and a specular reflection component of reflected light in the present embodiment. FIG. 本実施形態における撮影ブースを、正6面体で構成する例を示す説明図である。FIG. 13 is an explanatory diagram showing an example of a photography booth in this embodiment configured as a regular hexahedron. 本実施形態に係る偏光撮影装置のカメラ撮像板のRGBベイヤ配列を示す説明図である。2 is an explanatory diagram showing an RGB Bayer array of a camera imaging plate of the polarization imaging device according to the present embodiment. FIG. 被写体表面の凹凸の状態を示す粗さと、鏡面反射成分のぼける度合いとの関係を示す説明図である。1 is an explanatory diagram showing the relationship between roughness, which indicates the uneven state of the subject surface, and the degree of blurring of the specular reflection component; 本実施形態に係る偏光照明装置において、コード化照明を被写体に照射する例を示す説明図である。1 is an explanatory diagram showing an example in which coded illumination is irradiated onto an object in the polarized illumination device according to the present embodiment. FIG. 被写体表面の凹凸の状態を示す粗さと、コード化照明によるRGBの色ベクトルの分散との関係を示す説明図である。1 is an explanatory diagram showing the relationship between roughness, which indicates the uneven state of the subject surface, and the dispersion of RGB color vectors due to coded illumination; RGB空間上における鏡面反射光の色ベクトルの分布を示す説明図である。FIG. 10 is an explanatory diagram showing the distribution of color vectors of specularly reflected light in RGB space. 主成分分析に基づき算出した粗さσを、<A>滑らかな表面の被写体、<B>粗さ(小)の表面の被写体、<C>粗さ(大)の表面の被写体、のそれぞれについて算出したグラフである。13 is a graph showing roughness σ calculated based on principal component analysis for <A> a subject with a smooth surface, <B> a subject with a surface having low roughness, and <C> a subject with a surface having high roughness. 本実施形態における反射光の拡散反射成分と鏡面反射成分の分離に関する説明図である。4 is an explanatory diagram relating to separation of a diffuse reflection component and a specular reflection component of reflected light in the present embodiment. FIG. ランバートの余弦則に基づく、反射モデルの拡散反射成分に関する説明図である。FIG. 2 is an explanatory diagram of the diffuse reflection component of a reflection model based on Lambert's cosine law. Phong反射モデルに基づく鏡面反射に関する説明図である。FIG. 2 is an explanatory diagram of specular reflection based on the Phong reflection model. 本実施形態に係るコード化照明の動的パターンの例を示す説明図である。FIG. 11 is an explanatory diagram showing an example of a dynamic pattern of coded illumination according to the present embodiment. 本実施形態に係るコード化照明の動的パターンにおける他の例を示す説明図である。FIG. 13 is an explanatory diagram showing another example of a dynamic pattern of coded illumination according to the embodiment. 本実施形態における撮影ブースを、ジオディシックドーム形状で構成する例を示す説明図である。FIG. 1 is an explanatory diagram showing an example of a photography booth in this embodiment configured in a geodesic dome shape. 本実施形態における撮影ブースを、12面体で構成する例を示す説明図である。FIG. 13 is an explanatory diagram showing an example of a photography booth in this embodiment configured as a dodecahedron. 本実施形態に係る質感取得装置の処理の流れを示すフローチャートである。4 is a flowchart showing a processing flow of the texture acquisition device according to the present embodiment.

以下、本発明を実施するための形態(以下、「実施形態」と称する。)について図面を参照して説明する。
<質感取得システムの構成>
まず、図1を参照して、質感取得装置1を含む質感取得システム1000の概要を先に説明し、その後各構成について説明する。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Hereinafter, modes for carrying out the present invention (hereinafter, referred to as "embodiments") will be described with reference to the drawings.
<Configuration of the texture acquisition system>
First, with reference to FIG. 1, an overview of a texture acquisition system 1000 including a texture acquisition device 1 will be described, and then each component will be described.

質感取得システム1000は、図1に示すように、被写体を撮影する撮影ブースBHに配置される、奥行き計測装置20、偏光照明装置30および偏光撮影装置40と、その奥行き計測装置20、偏光照明装置30および偏光撮影装置40に通信接続される質感取得装置1とを備えている。 As shown in FIG. 1, the texture acquisition system 1000 includes a depth measurement device 20, a polarized lighting device 30, and a polarized photography device 40 that are arranged in a photography booth BH where a subject is photographed, and a texture acquisition device 1 that is communicatively connected to the depth measurement device 20, the polarized lighting device 30, and the polarized photography device 40.

本実施形態に係る質感取得装置1は、AR/VRに利用する被写体情報(3次元モデル)の計測や、3次元テレビにおける被写体情報のキャプチャ手法として有用なものであり、被写体の質感を取得する処理(以下、「質感取得処理」と称する。)として、(1)拡散反射成分・鏡面反射成分の分離処理、(2)被写体表面の凹凸状態の計測処理、(3)照度差ステレオによる被写体表面の詳細な形状情報の算出処理、を実現する。 The texture acquisition device 1 according to this embodiment is useful for measuring subject information (three-dimensional model) for use in AR/VR and as a method for capturing subject information in three-dimensional televisions, and performs the following processes for acquiring the texture of a subject (hereinafter referred to as "texture acquisition process"): (1) separation process of diffuse reflection components and specular reflection components, (2) measurement process of the unevenness of the subject surface, and (3) calculation process of detailed shape information of the subject surface using photometric stereo.

(1)拡散反射成分・鏡面反射成分の分離処理において、質感取得装置1は、図2で示すように、偏光照明装置30から偏光照明(例えば、赤照明、緑照明、青照明)による偏向光(照明光)を被写体Sに照射する。そして、照射された被写体Sの反射光を、例えば、偏光方向0度、45度、90度、135度それそれで選択的に透過する偏光撮影装置40で撮影する(後記する図4,図10参照)。撮影画像Pは、光源が映り込んだようなハイライトとして映る鏡面反射成分と、被写体内部に侵入した入射光が物体の色成分を吸収して拡散する拡散反射成分とを合わせた画像となる。質感取得装置1では、偏向光の被写体表面での反射挙動の違いを利用し、拡散反射成分Idと鏡面反射成分Isを分離することにより、拡散反射係数および鏡面反射係数を取得する(詳細は後記)。 (1) In the process of separating the diffuse reflection component and the specular reflection component, the texture acquisition device 1, as shown in Fig. 2, irradiates the subject S with polarized light (illumination light) by polarized illumination (e.g., red illumination, green illumination, and blue illumination) from the polarized illumination device 30. Then, the reflected light of the irradiated subject S is photographed by a polarized imaging device 40 that selectively transmits the light in the polarization directions of 0 degrees, 45 degrees, 90 degrees, and 135 degrees (see Figs. 4 and 10 described later). The photographed image P is an image that combines the specular reflection component that is reflected as a highlight as if a light source is reflected, and the diffuse reflection component that is the incident light that penetrates the inside of the subject and diffuses it by absorbing the color components of the object. The texture acquisition device 1 utilizes the difference in the reflection behavior of the polarized light on the surface of the subject to separate the diffuse reflection component Id and the specular reflection component Is , thereby acquiring the diffuse reflection coefficient and the specular reflection coefficient (details will be described later).

(2)被写体表面の凹凸状態の計測処理では、被写体に偏光照明(例えば、RGB(赤、緑、青)の3色)を照射する。鏡面反射の場合において、例えば誘電体の被写体表面では、入射角=反射角が成立するものしか撮影されないが、被写体表面に凹凸がある場合には、その凹凸により反射光の光ベクトルが広がり、その結果凹凸が混色として観測される。本実施形態では、被写体の凹凸の状態を示す粗さσを、注目する画素近傍の色ベクトルの分散状態から主成分分析の手法を用いて推定する(詳細は後記)。 (2) In the process of measuring the unevenness of the subject's surface, polarized lighting (e.g., three colors: RGB (red, green, and blue)) is irradiated onto the subject. In the case of specular reflection, for example, on the surface of a dielectric subject, only images are captured where the angle of incidence = angle of reflection. However, if the subject surface is uneven, the light vector of the reflected light spreads due to the unevenness, and as a result, the unevenness is observed as a mixed color. In this embodiment, the roughness σ, which indicates the unevenness of the subject, is estimated using a principal component analysis method from the variance state of the color vectors near the pixel of interest (details will be described later).

(3)照度差ステレオによる被写体表面の詳細な形状情報の算出処理では、被写体のサイズや計測する表面の凹凸のスケールに対して、撮影ブースBHの照明を動的に変化させる。照度差ステレオは、複数の異なる光源方向から撮影した被写体表面の輝度に基づき法線を推定することにより面の傾きを求め、その被写体の形状を計測する手法である。
質感取得装置1は、(1)拡散反射成分・鏡面反射成分の分離処理において算出した、拡散反射成分Idおよび拡散反射係数を用いて、撮影ブースBHの照明を動的に変化させて撮影することにより、従来の手法に比べ、被写体表面のより詳細な形状を取得することができる。
(3) In the process of calculating detailed shape information of the subject surface using photometric stereo, the lighting in the photography booth BH is dynamically changed according to the size of the subject and the scale of the unevenness of the surface to be measured. Photometric stereo is a method of measuring the shape of the subject by estimating the normal based on the luminance of the subject surface photographed from multiple different light source directions to determine the inclination of the surface.
The texture acquisition device 1 (1) uses the diffuse reflection component Id and the diffuse reflection coefficient calculated in the process of separating the diffuse reflection component and the specular reflection component to dynamically change the lighting in the photography booth BH while photographing, thereby making it possible to acquire a more detailed shape of the subject's surface than with conventional methods.

このように、質感取得装置1は、従来のヴォリュメトリックキャプチャ技術を用いた手法では取得できない、正確な被写体表面での拡散反射係数や鏡面反射係数、被写体表面の凹凸の状態(粗さσ)、および、照度差ステレオに基づく被写体表面の詳細な形状情報を取得することが可能となる。そして、質感取得装置1は、このようにして取得した3次元モデルの詳細な情報(後記する「3次元モデル詳細情報」)を、図示を省略したCG描画装置などに出力する。
以下、本実施形態に係る質感取得装置1を含む質感取得システム1000の詳細について説明する。
In this way, the texture acquisition device 1 is able to acquire accurate diffuse reflection coefficients and specular reflection coefficients on the subject's surface, the uneven state of the subject's surface (roughness σ), and detailed shape information of the subject's surface based on photometric stereo, which cannot be acquired by methods using conventional volumetric capture technology.The texture acquisition device 1 then outputs detailed information of the three-dimensional model acquired in this way ("detailed three-dimensional model information" described below) to a CG drawing device (not shown) or the like.
Hereinafter, the texture acquisition system 1000 including the texture acquisition device 1 according to this embodiment will be described in detail.

<質感取得装置の構成>
本実施形態に係る質感取得装置1は、一般的なコンピュータにより構成され、図1に示すように、制御部10と、入出力部11と、記憶部12とを備える。
入出力部11は、情報の送受信を行うための通信インタフェース、および、キーボード等の入力装置や、モニタ等の出力装置との間で情報の送受信を行うための入出力インタフェースからなる。
<Configuration of texture acquisition device>
The texture acquisition device 1 according to this embodiment is configured by a general computer, and includes a control unit 10, an input/output unit 11, and a storage unit 12, as shown in FIG.
The input/output unit 11 includes a communication interface for transmitting and receiving information, and an input/output interface for transmitting and receiving information to and from an input device such as a keyboard and an output device such as a monitor.

記憶部12は、フラッシュメモリやハードディスク、RAM(Random Access Memory)等により構成される。この記憶部12には、質感取得装置1による質感取得処理の初期設定情報200(詳細は後記)が格納される。
また、記憶部12には、制御部10の各機能を実行させるためのプログラム(質感取得プログラム)や、制御部10の処理に必要な情報が一時的に記憶される。
The storage unit 12 is configured with a flash memory, a hard disk, a RAM (Random Access Memory), etc. The storage unit 12 stores initial setting information 200 (details of which will be described later) for the texture acquisition process by the texture acquisition device 1.
The storage unit 12 also temporarily stores a program (texture acquisition program) for executing each function of the control unit 10 and information required for processing by the control unit 10.

制御部10は、質感取得装置1の制御全般を司り、図1に示すように、初期設定情報取得部110、被写体形状取得部120、計測制御部130、質感計算部140、3次元モデル詳細情報出力部150を含んで構成される。
この制御部10は、例えば、記憶部12に記憶されたプログラム(質感取得プログラム)を、図示を省略したCPU(Central Processing Unit)が、RAMに展開し実行することにより実現される。
The control unit 10 is responsible for overall control of the texture acquisition device 1, and as shown in FIG. 1, is composed of an initial setting information acquisition unit 110, a subject shape acquisition unit 120, a measurement control unit 130, a texture calculation unit 140, and a three-dimensional model detailed information output unit 150.
The control unit 10 is realized, for example, by a CPU (Central Processing Unit) (not shown) expanding a program (texture acquisition program) stored in the storage unit 12 into a RAM and executing it.

初期設定情報取得部110は、質感取得装置1が質感取得処理を実行するに際し、撮影ブースBHに設置した各装置の照明位置、配光特性、カメラ姿勢を含む情報を、初期設定情報200として予め取得し、記憶部12に記憶しておく。この初期設定情報取得部110が初期設定情報200を取得する一例を、図3を参照して説明する。
図3は、撮影ブースBHを、正6面体で構成する例を示している。奥行き計測装置20、偏光照明装置30および偏光撮影装置40は、被写体を取り囲んで配置されている。図3に示すように、例えば、正6面体の頂点および辺の中央に、奥行き計測装置(奥行カメラ)20および偏光撮影装置40(例えば、偏光RGBカメラ)を配置し、撮影ブースBHの内部に置かれた被写体を撮影する。偏光照明装置30は、正6面体の面ごとに配置される。ここでは、フラットパネルディスプレイ(面光源)で構成し、直線偏光フィルタを配置しているものとする。
なお、正6面体は一例であり、これに限定されない(後記する図15,図16参照)。
When the texture acquisition device 1 executes the texture acquisition process, the initial setting information acquisition unit 110 acquires information including the lighting positions, light distribution characteristics, and camera postures of the devices installed in the photography booth BH in advance as initial setting information 200, and stores the information in the storage unit 12. An example of how the initial setting information acquisition unit 110 acquires the initial setting information 200 will be described with reference to FIG.
Fig. 3 shows an example in which the photography booth BH is configured as a regular hexahedron. The depth measurement device 20, the polarized illumination device 30, and the polarized photography device 40 are arranged to surround the subject. As shown in Fig. 3, for example, the depth measurement device (depth camera) 20 and the polarized photography device 40 (e.g., a polarized RGB camera) are arranged at the vertices and the center of the sides of the regular hexahedron to photograph the subject placed inside the photography booth BH. The polarized illumination device 30 is arranged on each face of the regular hexahedron. Here, it is assumed that the photography booth BH is configured with a flat panel display (surface light source) and a linear polarizing filter is arranged.
Note that a regular hexahedron is just one example and is not limited to this (see FIGS. 15 and 16 described later).

初期設定情報取得部110は、偏光照明装置30が配置される照明位置、空間上の各方向への光度等を測定した配光特性、奥行き計測装置20および偏光撮影装置40が配置されるカメラ位置(X,Y,Z座標)やカメラ姿勢(ロール・ピッチ・ヨー)の情報を、初期設定情報200として外部装置等から取得し、記憶部12に記憶する。
なお、撮影ブースBHの1箇所に配置される奥行き計測装置20と偏光照明装置30との組は、同じカメラ位置であり同じカメラ姿勢であるものとする。
The initial setting information acquisition unit 110 acquires information on the lighting position where the polarized lighting device 30 is placed, the light distribution characteristics measuring the luminous intensity in each direction in space, and the camera positions (X, Y, Z coordinates) and camera attitudes (roll, pitch, yaw) where the depth measurement device 20 and the polarized shooting device 40 are placed from an external device, etc., as initial setting information 200, and stores it in the memory unit 12.
It is assumed that the pair of depth measurement device 20 and polarized illumination device 30 arranged at one location in the photography booth BH has the same camera position and the same camera attitude.

図1に戻り、被写体形状取得部120は、計測制御部130を介して、撮影ブースBHに配置された各奥行き計測装置20に計測を実行させる指示情報を送信する。そして、被写体形状取得部120は、各奥行き計測装置20からの計測情報を受信することにより、(大まかな)被写体形状を取得する。
なお、奥行き計測装置20は、例えば、カメラから近赤外線を発光し、物体に反射して戻ってくるまでの時間を測定することで奥行きを算出するToF(Time of Flight)方式のカメラを用いることができる。
1 , the subject shape acquisition unit 120 transmits instruction information for causing each depth measurement device 20 arranged in the photography booth BH to perform measurement via the measurement control unit 130. Then, the subject shape acquisition unit 120 receives measurement information from each depth measurement device 20 to acquire a (rough) subject shape.
The depth measurement device 20 may be, for example, a Time of Flight (ToF) camera that emits near-infrared light, reflects off an object, and measures the time it takes for the light to return to calculate depth.

計測制御部130は、撮影ブースBHに配置された、奥行き計測装置20、偏光照明装置30および偏光撮影装置40に対し、被写体形状取得部120や質感計算部140と連携して、指示情報を送信することにより各装置を制御する。
具体的には、計測制御部130は、奥行き計測装置20に対して指示情報を送信することにより被写体形状の計測を実行させる。
The measurement control unit 130 cooperates with the subject shape acquisition unit 120 and the texture calculation unit 140 to control each device by sending instruction information to the depth measurement device 20, polarized lighting device 30, and polarized photography device 40 arranged in the photography booth BH.
Specifically, the measurement control unit 130 transmits instruction information to the depth measurement device 20 to cause it to measure the subject shape.

また、計測制御部130は、質感計算部140が質感取得処理を実行する際に、偏光照明装置30に対し、後記する図13や図14で示すような所定パターンの照明状態になるように制御する指示情報を送信する。
図13に示す所定のパターンは、被写体サイズや計測する表面の凹凸のスケールに対し、照明の状態を動的に変更するものであり、ここでは、RGB(赤、緑、青)の照明を所定のサイズの区画ごとに分けてコード化して照明する。なお、各色の区画において斜線で示す線は、直線偏光フィルタの偏光方向を示しており、ここでは、偏光方向0度、45度、90度、135度の偏向光が照射されることを示している。
計測制御部130は、このコード化照明の動的パターンを制御する指示情報を偏光照明装置30に送信する。
In addition, when the texture calculation unit 140 executes the texture acquisition process, the measurement control unit 130 transmits instruction information to the polarized lighting device 30 to control the device to achieve a predetermined lighting pattern as shown in Figures 13 and 14 described below.
The predetermined pattern shown in Fig. 13 dynamically changes the lighting conditions according to the size of the object and the scale of the unevenness of the surface to be measured, and here, RGB (red, green, blue) lighting is divided into sections of a predetermined size and coded for illumination. Note that the hatched lines in each color section indicate the polarization direction of the linear polarizing filter, and here indicate that polarized light with polarization directions of 0 degrees, 45 degrees, 90 degrees, and 135 degrees is irradiated.
The measurement control unit 130 transmits instruction information for controlling the dynamic pattern of this coded illumination to the polarized illumination device 30 .

また、計測制御部130は、質感計算部140が質感取得処理を実行し、偏光照明装置30による被写体への照明光の照射が行われることに対応して、偏光撮影装置40に反射光の撮影を実行する指示情報を送信する。 In addition, the measurement control unit 130 transmits instruction information to the polarized photographing device 40 to photograph the reflected light in response to the texture calculation unit 140 executing the texture acquisition process and the polarized lighting device 30 irradiating the subject with illumination light.

偏光撮影装置40(ここでは、「偏光RGBカメラ」とする。)のカメラ撮像板は、例えば図4に示すように、RGBベイヤ配列となっている。さらに、RGBの各画素は、4分割の画素で構成され、それぞれ偏光方向0度、45度、90度、135度の成分を選択的に透過する画素となっている。なお、図4においては、斜線で示すハッチの方向が透過する偏向光成分の方向を示している。 The camera imaging plate of the polarization imaging device 40 (herein referred to as a "polarized RGB camera") has an RGB Bayer array, for example as shown in FIG. 4. Furthermore, each RGB pixel is composed of four divided pixels, and each pixel selectively transmits components with polarization directions of 0 degrees, 45 degrees, 90 degrees, and 135 degrees. Note that in FIG. 4, the direction of the hatched area indicated by diagonal lines indicates the direction of the polarized light components that are transmitted.

質感計算部140は、質感取得処理として、上記した、拡散反射成分・鏡面反射成分の分離処理、被写体表面の凹凸状態の計測処理、および、照度差ステレオによる被写体表面の詳細な形状情報の算出処理を実行する。この質感計算部140は、凹凸状態算出部141と、反射係数算出部142と、照度差ステレオ形状取得部(表面形状取得部)143とを備える。 As part of the texture acquisition process, the texture calculation unit 140 performs the above-mentioned separation process of diffuse reflection components and specular reflection components, measurement process of the unevenness of the subject's surface, and calculation process of detailed shape information of the subject's surface using photometric stereo. This texture calculation unit 140 includes a roughness calculation unit 141, a reflection coefficient calculation unit 142, and a photometric stereo shape acquisition unit (surface shape acquisition unit) 143.

凹凸状態算出部141は、被写体表面の凹凸により生じた偏向光の色ベクトルの分散から、被写体表面の凹凸の状態を示す粗さσを算出する。
被写体の凹凸の状態は、凹凸の細かさなどの多様なパラメータで表現できるが、光源などの映り込むものの状態が被写体表面でぼやける度合いとして扱うことができる。
例えば、図5に示すように、符号51で示すRGBの照明を、球状の被写体52に照射し、鏡面反射成分を撮影すると、符号53で示すように、被写体52の表面が粗いほど、RGB各色の成分のぼける程度が大きくなる。つまり、色ベクトルの分散が大きいほど、被写体表面の凹凸の粗さが大きいことが示される。
The unevenness calculating section 141 calculates roughness σ, which indicates the unevenness of the object surface, from the variance of the color vector of the deflected light caused by the unevenness of the object surface.
The unevenness of a subject can be expressed by various parameters such as the fineness of the unevenness, but it can also be treated as the degree to which reflected objects such as light sources are blurred on the subject's surface.
5, when RGB illumination indicated by reference numeral 51 is irradiated onto a spherical subject 52 and the specular reflection component is photographed, the degree to which each of the RGB color components is blurred increases as the surface of the subject 52 becomes rougher, as indicated by reference numeral 53. In other words, the greater the variance of the color vectors, the greater the roughness of the subject surface.

本実施形態では、偏光照明を面光源とし、RGBの3色を用いて撮影ブース壁面でコード化し、これを照明(コード化照明γl)として被写体に照射する例として説明する(図6参照)。
図6では、RGBの3色をそれぞれ0度,45度,90度のように異なる偏向光(γl1,γl2,γl3)として、被写体に照射する例を示している。これは、単一方向では、入射角αがブリュースター角となり計測不可能になるケースを避けるためである。なお、偏向光は、0度,45度,90度に限定されるものではない。
また、本発明は、偏光照明をRGBの3色に限定するものではなく、複数色の偏向光であればよく、5色分光や多色分光であっても適応可能である。
In this embodiment, an example will be described in which polarized illumination is used as a surface light source, and the three colors RGB are used to code on the wall of the photography booth, and this is irradiated onto the subject as illumination (coded illumination γ l ) (see FIG. 6).
6 shows an example in which the three colors of RGB are irradiated onto the subject as differently polarized light (γ l1 , γ l2 , γ l3 ) at angles of 0 degrees, 45 degrees, and 90 degrees, respectively. This is to avoid cases in which the incident angle α becomes Brewster's angle and measurement becomes impossible in a single direction. Note that the polarized light is not limited to 0 degrees, 45 degrees, and 90 degrees.
Furthermore, the present invention is not limited to using polarized illumination in three colors, RGB, but may use polarized light in a plurality of colors, and may also be applicable to five-color or polychromatic light.

被写体表面のある部位に着目した場合、鏡のようにフラットな表面であるならば、入射角=反射角が成立するものしかカメラ(偏光RGBカメラ40a)には映り込まない(図7の符号71参照)。しかしながら、被写体表面に凹凸がある場合には、その粗さの程度に応じて、着目した画素位置において、RGBが混色として観測される度合いが高まる、つまり、被写体表面の部位(微小面構造)の粗さの程度が大きい程、色ベクトルの分散が大きくなる(図7の符号72,73参照)。本実施形態では、この被写体表面の粗さと、色ベクトルの分散との関連性を用いて、色ベクトルの分散から被写体表面の粗さを推定する。 When focusing on a certain part of the subject's surface, if the surface is flat like a mirror, only those parts for which the angle of incidence = the angle of reflection are captured by the camera (polarized RGB camera 40a) (see reference number 71 in FIG. 7). However, if the subject's surface is uneven, the degree to which RGB is observed as a mixture of colors at the focused pixel position increases depending on the degree of roughness. In other words, the greater the degree of roughness of a part of the subject's surface (microsurface structure), the greater the variance of the color vectors (see reference numbers 72 and 73 in FIG. 7). In this embodiment, the relationship between the roughness of the subject's surface and the variance of the color vectors is used to estimate the roughness of the subject's surface from the variance of the color vectors.

具体的には、凹凸状態算出部141は、被写体の凹凸の状態を示す粗さσを、注目する画素近傍の色ベクトルの分散状態から主成分分析の手法を用いて推定する。
図8は、RGB空間上における鏡面反射光の色ベクトルの分布を示す図である。注目した被写体の部位周辺の鏡面反射光を、偏光RGBカメラ40a(図6参照)で撮影し、RGB空間上に鏡面反射光をプロットしたものである。図8(a)は、被写体が滑らかな表面の場合の色ベクトルの分布を示す。図8(b)は、被写体が粗い表面の場合の色ベクトルの分布を示す。
Specifically, the unevenness calculation section 141 estimates the roughness σ, which indicates the unevenness of the subject, from the variance state of color vectors in the vicinity of a pixel of interest by using a principal component analysis technique.
Fig. 8 shows the distribution of color vectors of specular reflected light in RGB space. The specular reflected light around a part of a subject of interest was photographed with a polarized RGB camera 40a (see Fig. 6), and the specular reflected light was plotted in RGB space. Fig. 8(a) shows the distribution of color vectors when the subject has a smooth surface. Fig. 8(b) shows the distribution of color vectors when the subject has a rough surface.

そして、凹凸状態算出部141は、RGB空間上の鏡面反射光色ベクトルに対し主成分分析を適用し、第1主成分寄与率と第3主成分寄与率を算出する。凹凸状態算出部141は、算出した第1主成分寄与率と第3主成分寄与率を用いて、図8に示す式(1)により、粗さσを算出する。 Then, the unevenness calculation unit 141 applies principal component analysis to the specular reflected light color vector in the RGB space to calculate the first principal component contribution rate and the third principal component contribution rate. The unevenness calculation unit 141 uses the calculated first principal component contribution rate and third principal component contribution rate to calculate the roughness σ according to formula (1) shown in FIG. 8.

ここで、式(1)によれば、第1主成分寄与率と第3主成分寄与率の値が近いほど、粗さσが「1」に近くなる。つまり、鏡面反射光の色ベクトルの分散が大きいほど、色ベクトルのばらつきが大きく、第1主成分寄与率と第3主成分寄与率の値の差が小さくなり、粗さσが「1」に近くなる。このようにして、粗さσの状態を、主成分分析により得られた第1主成分寄与率と第3主成分寄与率とを用いて定量化することができる。 Here, according to formula (1), the closer the values of the first principal component contribution rate and the third principal component contribution rate are, the closer the roughness σ is to "1". In other words, the greater the variance of the color vector of the specular reflected light, the greater the variation in the color vector, the smaller the difference between the values of the first principal component contribution rate and the third principal component contribution rate, and the closer the roughness σ is to "1". In this way, the state of roughness σ can be quantified using the first principal component contribution rate and the third principal component contribution rate obtained by principal component analysis.

図9は、主成分分析に基づき算出した粗さσを、<A>滑らかな表面の被写体、<B>粗さ(小)の表面の被写体、<C>粗さ(大)の表面の被写体、のそれぞれについて算出したグラフである。
このグラフにおいて、粗さσの平均値は、<A>滑らかな表面の被写体において「0.55」、<B>粗さ(小)の表面の被写体において「0.82」、<C>粗さ(大)の表面の被写体において「0.91」となった。このように、粗さσは、値が小さいほど凹凸が少なく滑らかな表面であることを示し、値が大きいほど凹凸が大きく粗い状態の表面であることを示している。なお、ここでは、注目画素の周囲±3画素のサイズで主成分分析を行っている。
FIG. 9 is a graph showing the roughness σ calculated based on the principal component analysis for <A> an object with a smooth surface, <B> an object with a surface having low roughness, and <C> an object with a surface having high roughness.
In this graph, the average roughness σ is 0.55 for the subject with a smooth surface (A), 0.82 for the subject with a small roughness (B), and 0.91 for the subject with a large roughness (C). Thus, the smaller the roughness σ value, the smoother the surface is, and the larger the value, the rougher the surface is. Note that here, principal component analysis is performed on a size of ±3 pixels around the pixel of interest.

このように、凹凸状態算出部141は、主成分分析により算出した寄与率(第1主成分寄与率と第3主成分寄与率)を用いて、被写体表面の注目部分の粗さσを算出することができる。 In this way, the unevenness calculation unit 141 can calculate the roughness σ of the target portion of the subject's surface using the contribution rates (the first principal component contribution rate and the third principal component contribution rate) calculated by principal component analysis.

図1に戻り、反射係数算出部142は、偏向光の被写体表面での反射挙動の違いを利用して、反射光における拡散反射成分Idと鏡面反射成分Isの分離を行う。そして、反射係数算出部142は、取得した拡散反射成分Idと、被写体への入射光量と、被写体の形状情報と、照明位置とを用いて、被写体の拡散反射係数Kdを算出する。また、反射係数算出部142は、取得した鏡面反射成分Isと、被写体への入射光量と、粗さσと、カメラ位置とに基づき、鏡面反射係数Ksを算出する。以下、具体的に説明する。 Returning to Fig. 1, the reflection coefficient calculation unit 142 separates the reflected light into a diffuse reflection component Id and a specular reflection component Is by utilizing the difference in the reflection behavior of the polarized light on the surface of the subject. The reflection coefficient calculation unit 142 then calculates a diffuse reflection coefficient Kd of the subject using the acquired diffuse reflection component Id , the amount of light incident on the subject, shape information of the subject, and the illumination position. The reflection coefficient calculation unit 142 also calculates a specular reflection coefficient Ks based on the acquired specular reflection component Is, the amount of light incident on the subject, roughness σ, and the camera position. This will be described in detail below.

≪拡散反射成分と鏡面反射成分の分離処理≫
まず、図10を参照して、反射光における拡散反射成分と鏡面反射成分の分離について説明する。
ここでは、偏光照明装置30を面光源とし、計測制御部130の制御により、RGBの3色を用いて撮影ブース壁面でコード化し、照明として被写体に照射する例を示す。その際、偏光撮影装置40(偏光RGBカメラ40a)により被写体を撮影する。
この偏光RGBカメラ40aのカメラ撮像板は、図10の符号101に示すように、RGBベイヤ配列となっており、さらに、RGBの各画素は、4分割の画素で構成され、それぞれ偏光方向0度,45度,90度,135度の成分を選択的に透過する画素となっている。なお、4分割の画素および偏光方向0度,45度,90度,135度は一例であり、これに限定されない。
<Separation of diffuse reflection and specular reflection components>
First, with reference to FIG. 10, separation of the diffuse reflection component and the specular reflection component in the reflected light will be described.
Here, the polarized lighting device 30 is used as a surface light source, and the three colors RGB are used to code the light on the wall of the photography booth under the control of the measurement control unit 130, and the light is irradiated onto the subject as illumination. At this time, the subject is photographed by the polarized photography device 40 (polarized RGB camera 40a).
The camera imaging plate of this polarization RGB camera 40a has an RGB Bayer array as shown by reference numeral 101 in Fig. 10, and each RGB pixel is composed of a pixel divided into four, and each pixel selectively transmits components with polarization directions of 0 degrees, 45 degrees, 90 degrees, and 135 degrees. Note that the pixel divided into four and the polarization directions of 0 degrees, 45 degrees, 90 degrees, and 135 degrees are merely examples, and the present invention is not limited to these.

観測される輝度Iは、二色性反射モデルに基づき、次の式(2)により表される。なお、二色性反射モデルは、反射光のスペクトルは拡散反射成分と鏡面反射成分のスペクトルの線形和で表すことができる、というものである。 The observed luminance I is expressed by the following equation (2) based on the dichromatic reflection model. The dichromatic reflection model states that the spectrum of reflected light can be expressed as a linear sum of the spectra of the diffuse reflection component and the specular reflection component.

I=Id +a(1+cos(2Θ-β)) ・・・式(2)
ここで、Id は、拡散反射成分、a(1+cos(2Θ-β))が、鏡面反射成分Isを表す。また、aは、振幅を示し、Θは画素の偏光方向を示す。βは位相である。
I=I d +a(1+cos(2Θ-β))...Formula (2)
Here, Id is the diffuse reflection component, a(1+cos(2Θ-β)) is the specular reflection component Is , a is the amplitude, Θ is the polarization direction of the pixel, and β is the phase.

偏光方向0度,45度,90度,135度の画素で得られた輝度値(それぞれ、I0 ,I45 ,I90,I135とする。)は、図10に示すように、その平均値Iaveに対して正弦波で表される挙動を示すが、偏光角度の倍角で正弦波となる。そこで、三角関数の公式sin(θ-β)+cos(θ-β)=1(βは位相を考慮)を用いて、例えば、I0 +I45 =a(aは三角関数の振幅)となる。これにより、振幅aを算出することができる。また、観測された輝度Iの最小値Iminが拡散反射成分Idであり、振幅aの2倍である2aが鏡面反射成分Isとなる。 The luminance values obtained from pixels with polarization directions of 0 degrees, 45 degrees, 90 degrees, and 135 degrees (represented by I 0 , I 45 , I 90 , and I 135 , respectively) behave as sine waves with respect to their average value I ave , as shown in FIG. 10, but become sine waves at double the polarization angle. Therefore, using the trigonometric formula sin(θ-β) 2 +cos(θ-β) 2 =1 (β takes into account the phase), for example, I 0 2 +I 45 2 =a (a is the amplitude of the trigonometric function). This makes it possible to calculate the amplitude a. Furthermore, the minimum value I min of the observed luminance I is the diffuse reflection component I d , and 2a, which is twice the amplitude a, is the specular reflection component I s .

次に、反射光の拡散反射成分においては、Lambert反射モデル(ランバートの余弦則に基づく拡散面の反射)を例にして、拡散反射係数を算出する処理を説明する。また、鏡面反射成分においては、Phong反射モデルを例にして、鏡面反射係数を算出する処理を説明する。 Next, for the diffuse reflection component of reflected light, we will explain the process of calculating the diffuse reflection coefficient using the Lambert reflection model (reflection of a diffuse surface based on Lambert's cosine law) as an example. Also, for the specular reflection component, we will explain the process of calculating the specular reflection coefficient using the Phong reflection model as an example.

≪拡散反射係数の算出≫
反射係数算出部142による、拡散反射係数Kdの算出処理について説明する。
反射係数算出部142は、反射光の強さが入射角の余弦に比例する、というランバートの余弦則に基づき、拡散反射係数Kdを算出する。
<Calculation of diffuse reflectance coefficient>
The calculation process of the diffuse reflection coefficient Kd by the reflection coefficient calculation unit 142 will be described.
The reflection coefficient calculation unit 142 calculates the diffuse reflection coefficient Kd based on Lambert's cosine law, which states that the intensity of reflected light is proportional to the cosine of the angle of incidence.

図11に示すように、被写体に対し角度θで入射した入射光量I0の拡散反射光は、ランバートの余弦則に基づき、入射光と面の法線Nとの間の角度θの余弦に正比例する。このときの拡散反射係数Kdは、上記した拡散反射成分と鏡面反射成分の分離処理において算出されたId、被写体の形状情報から求まる法線N、光源位置(照明位置)、入射光量I0が既知であれば、下記の式(3)により求めることができる。
d = Kd0 cosθ ・・・式(3)
11, the diffuse reflection light of the incident light amount I0 incident on the subject at an angle θ is directly proportional to the cosine of the angle θ between the incident light and the normal line N of the surface, based on Lambert's cosine law. The diffuse reflection coefficient Kd at this time can be calculated by the following formula (3) if Id calculated in the above-mentioned process of separating the diffuse reflection component and the specular reflection component, the normal line N calculated from the shape information of the subject, the light source position (illumination position), and the incident light amount I0 are known.
I d = K d I 0 cosθ...Formula (3)

≪鏡面反射係数の算出≫
反射係数算出部142による、鏡面反射係数Ksの算出処理について説明する。
図12は、Phong反射モデルに基づく鏡面反射についての説明図である。
ここで、入射光(入射光量I0)が正反射する方向に対し、観測(撮影)する視点がαずれている場合、鏡面反射成分であるその反射光の強度Isは、下記の式(4)により求めることができる。なお、観測(撮影)する視点のずれαは、カメラ位置から求まる。
s = Ks0 cosnα ・・・式(4)
<Calculation of specular reflection coefficient>
The calculation process of the specular reflection coefficient Ks by the reflection coefficient calculation unit 142 will be described.
FIG. 12 is an explanatory diagram of specular reflection based on the Phong reflection model.
If the observation (photography) viewpoint is shifted by α from the direction in which the incident light (incident light amount I0 ) is specularly reflected, the intensity I s of the reflected light, which is the specular reflection component, can be calculated by the following formula (4). Note that the shift α of the observation (photography) viewpoint is calculated from the camera position.
I s = K s I 0 cos n α ...Formula (4)

式(4)における三角関数cosのべき数nは、被写体の表面の粗さに応じて決まる値である。
本実施形態においては、べき数nと、上記のようにして計算した粗さσとの関係を、以下の手法により紐づける。ここでは、注目する被写体の部位を単一視点から得た情報により、べき数nおよび鏡面反射係数Ksを求める〔手法1〕と、複数視点から得た情報により、べき数nおよび鏡面反射係数Ksを求める〔手法2〕について説明する。
The exponent n of the trigonometric function cos in equation (4) is a value determined according to the roughness of the surface of the subject.
In this embodiment, the relationship between the exponent n and the roughness σ calculated as above is linked by the following methods. Here, we will explain [Method 1], which calculates the exponent n and the specular reflection coefficient Ks from information obtained from a single viewpoint of a part of a subject of interest, and [Method 2], which calculates the exponent n and the specular reflection coefficient Ks from information obtained from multiple viewpoints.

〔手法1〕
反射係数算出部142は、被写体の部位を単一視点から得られた粗さσ(ここではσ1とする。)を用いて、事前に様々な表面の粗さの被写体を観測することにより、例えば、以下の手法1-1~手法1-3により、べき数nおよび鏡面反射係数Ksを求める。
[Method 1]
The reflection coefficient calculation unit 142 uses the roughness σ (here, σ 1 ) obtained from a single viewpoint of a part of the subject to observe subjects having various surface roughnesses in advance, and calculates the exponent n and the specular reflection coefficient Ks , for example, by the following methods 1-1 to 1-3.

(手法1-1)反射係数算出部142は、事前に作成しておいたルックアップテーブルを用いて、べき数nを算出する。
事前に、α、σ1 、Is 、Ks 、nをルックアップテーブルとして作成しておく。
そして、凹凸状態算出部141が算出した粗さσ1と、事前の観測により得られたα、鏡面反射成分Isとから、ルックアップテーブルを参照し、最適なべき数nと鏡面反射係数Ksとを抽出する。
(Method 1-1) The reflection coefficient calculation unit 142 calculates the exponent n by using a lookup table created in advance.
α, σ 1 , I s , K s and n are created in advance as a look-up table.
Then, from the roughness σ 1 calculated by the unevenness calculation unit 141, α obtained by prior observation, and the specular reflection component I s , a lookup table is referenced to extract the optimal exponent n and specular reflection coefficient K s .

(手法1-2)3次元以上の高次方程式にパラメータフィッティングする。
上記した式(4)において、α、I0 、Is が計測条件および計測結果から既知となる。この式に、例えば以下の式(5)で求めたnを用いることで、鏡面反射係数Ksを算出する。
n=aσ1 2+bσ1+c ・・・式(5)
(Method 1-2) Parameter fitting to a three-dimensional or higher order equation.
In the above formula (4), α, I 0 , and I s are known from the measurement conditions and the measurement results. By using n obtained by the following formula (5) in this formula, for example, the specular reflection coefficient K s is calculated.
n=aσ 1 2 +bσ 1 +c ... Formula (5)

(手法1-3)学習済みニューラルネットワークを用いてKs 、nを推定する。
入力を、α、σ1 、Is とし、出力をKs 、nとしたDNN(Deep Neural Network)を用いて、事前に計測した学習用データセットにより重みを学習させる。この学習用データセットには、様々な被写体、様々な光源位置(照明位置)、カメラの計測位置で取得されたデータを用いる。
反射係数算出部142が実際にKs 、nを求める際には、学習済みのDNNにα、σ1 、Is を入力し、Ks 、nを出力させる。
(Method 1-3) Estimate K s and n using a trained neural network.
A deep neural network (DNN) with inputs of α, σ 1 , and I s and outputs of K s and n is used to learn weights from a training data set measured in advance. This training data set uses data acquired from various subjects, various light source positions (illumination positions), and camera measurement positions.
When the reflection coefficient calculation unit 142 actually finds K s and n, α, σ 1 and I s are input to a trained DNN, which outputs K s and n.

〔手法2〕
反射係数算出部142は、被写体の部位を複数の視点から得られた粗さσ(ここではσ1 ,σ2 とする。)を用いて、事前に様々な表面の粗さの被写体を観測することにより、例えば、以下の手法2-1,手法2-2により、べき数nおよび鏡面反射係数Ksを求める。
[Method 2]
The reflection coefficient calculation unit 142 uses roughness σ (here, σ 1 and σ 2 ) obtained from multiple viewpoints of a part of the subject to observe subjects with various surface roughnesses in advance, and calculates the exponent n and the specular reflection coefficient K s , for example, by the following methods 2-1 and 2-2.

(手法2-1)学習済みニューラルネットワークを用いてnを推定する。
非線形回帰問題とし、入力をα1 、α2 、σ1 、σ2 、出力をnとするDNNを作成し、事前に計測した学習用データセットを用いて重みを学習させる。
反射係数算出部142が実際にnを求める際には、学習済みのDNNにα1 、α2 、σ1 、σ2 を入力し、nを出力させる。出力したnを上記した式(4)に挿入して鏡面反射係数Ksを求める。
(Method 2-1) Estimate n using a trained neural network.
A nonlinear regression problem is considered, and a DNN with inputs α 1 , α 2 , σ 1 , σ 2 and output n is created, and weights are learned using a learning dataset measured in advance.
When the reflection coefficient calculation unit 142 actually calculates n, α 1 , α 2 , σ 1 , and σ 2 are input to the trained DNN, which outputs n. The output n is inserted into the above equation (4) to calculate the specular reflection coefficient K s .

(手法2-2)学習済みニューラルネットワークを用いてKs 、nを推定する。
入力を、α1 、α2 、σ1 、σ2 、Is1、Is2 とし、出力をKs 、nとしたDNNに対し、事前に計測した学習用データセットを用いて重みを学習させる。この学習用データセットには、様々な被写体、様々な光源位置(照明位置)、カメラの計測位置で取得されたデータを用いる。
反射係数算出部142が実際にKs 、nを求める際には、学習済みのDNNにα1 、α2 、σ1 、σ2 、Is1、Is2 を入力し、Ks 、nを出力させる。
(Method 2-2) Estimate K s and n using a trained neural network.
The inputs of a DNN are α1 , α2 , σ1 , σ2 , Is1 , and Is2 , and the outputs are Ks and n. Weights are trained using a training data set measured in advance. This training data set uses data acquired from various subjects, various light source positions (illumination positions), and camera measurement positions.
When the reflection coefficient calculation unit 142 actually finds K s ,n, α 1 , α 2 , σ 1 , σ 2 , I s1 , and I s2 are input to the trained DNN, which outputs K s ,n.

このようにして、反射係数算出部142は、被写体の拡散反射係数Kdおよび鏡面反射係数Ksを算出することができる。 In this manner, the reflection coefficient calculation section 142 can calculate the diffuse reflection coefficient K d and the specular reflection coefficient K s of the subject.

≪照明の動的制御≫
ここで、質感計算部140が計測制御部130を介して実行する、偏光照明装置30の動的制御について説明する。
本実施形態では、偏光照明装置30からの偏光照明を面光源とし、RGBの3色を用いてコード化した照明(コード化照明)を被写体に照射する(図6,図10参照)。その際に、計測制御部130は、被写体のサイズや計測する表面の凹凸のスケールに応じて、照明の状態を動的に変更する。これにより、照度差ステレオ形状取得部143において、詳細な被写体表面の形状が取得可能な照明条件を満たすようにする。
Dynamic lighting control
Here, the dynamic control of the polarized lighting device 30 that the texture calculation unit 140 executes via the measurement control unit 130 will be described.
In this embodiment, polarized light from the polarized lighting device 30 is used as a surface light source, and the subject is illuminated with coded light (coded light) using the three colors RGB (see Figs. 6 and 10). At this time, the measurement control unit 130 dynamically changes the lighting conditions according to the size of the subject and the scale of the unevenness of the surface to be measured. This satisfies the lighting conditions that allow the photometric stereo shape acquisition unit 143 to acquire the detailed shape of the subject surface.

図13を参照して、コード化照明の動的パターンの例を説明する。
計測制御部130は、図13で示すように、RGBの3色について、色、偏光方向、および、スケールのパターンをフレーム単位で変化させて照明を被写体に照射する。
図13に示す例では、全6フレームのうち、2フレームごとにスケールが、小・中・大と変更される。そして、その2フレームのうち、色と偏光方向との組み合わせが1フレームごとに変化する例を示している。なお、図13の「黒」は、照明が発光していない状態を示している。
このようにすることにより、サイズが異なる被写体を対象にする場合や、被写体の部位によって表面の凹凸のスケールが異なる場合であっても、より正確な計測が可能となる。
なお、動的パターンは、図13に示すような所定の規則性をもってパターンを変更する場合に限られない。異方性反射がない状態であれば、RGBの3色や、偏光方向、スケールについてランダムにコードを変更することも可能である。
With reference to FIG. 13, an example of a dynamic pattern of coded illumination is illustrated.
As shown in FIG. 13, the measurement control unit 130 changes the color, polarization direction, and scale pattern for the three colors RGB on a frame-by-frame basis to illuminate the subject.
In the example shown in Fig. 13, the scale is changed from small to medium to large every two frames out of a total of six frames. And, in the two frames, the combination of color and polarization direction changes every other frame. Note that "black" in Fig. 13 indicates a state where the illumination is not emitting light.
In this way, even when subjects of different sizes are targeted, or when the scale of the surface irregularities varies depending on the part of the subject, more accurate measurements are possible.
The dynamic pattern is not limited to changing the pattern with a predetermined regularity as shown in Fig. 13. As long as there is no anisotropic reflection, it is also possible to change the code randomly with respect to the three colors of RGB, the polarization direction, and the scale.

図14は、本実施形態に係るコード化照明の動的パターンにおける他の例を示す図である。
図14で示す例では、正六面体の撮影ブースの面光源において、基準頂点に隣接する3面でRGBの各色を照射し、6フレームで1シーケンスを構成する。また2フレーム単位で1セットとなっており、各色が次のフレームで正反対の面に移行する。なお、この1セットで同一色の光源方向を変えることにより、後記する照度差ステレオによる形状計測を行い、被写体表面の詳細な形状情報(後記する「被写体表面形状情報」)を取得する。
FIG. 14 is a diagram showing another example of a dynamic pattern of coded illumination according to the present embodiment.
In the example shown in Fig. 14, the surface light source of a regular hexahedral photography booth illuminates the three surfaces adjacent to the reference vertex with each of the RGB colors, and one sequence is made up of six frames. One set is made up of two frames, and each color moves to the opposite surface in the next frame. By changing the direction of the light source of the same color within this set, shape measurement is performed using photometric stereo, which will be described later, and detailed shape information of the subject's surface ("subject surface shape information", which will be described later) is obtained.

なお、本実施形態に係る撮影ブースは、正6面体に限定するものではなく、例えば、図15に示すようなジオディシックドーム形状や、図16に示すような12面体としてもよい。また、照明を、例えばフラットパネルディスプレイのような面光源で構成すれば、構成する一面を光源単位としてもよいし(図15(a)参照)、照明として提示するパターンを調整することにより、面単位で発行色を制御する必要はなく、被写体表面で任意に各色を映りこませることも可能である(図15(b)参照)。 The photography booth according to this embodiment is not limited to a regular hexahedron, and may be, for example, a geodesic dome shape as shown in FIG. 15 or a dodecahedron as shown in FIG. 16. If the lighting is configured as a surface light source such as a flat panel display, each surface may be treated as a light source unit (see FIG. 15(a)). By adjusting the pattern presented as lighting, it is not necessary to control the emitted color on a surface-by-surface basis, and it is also possible to reflect each color arbitrarily on the surface of the subject (see FIG. 15(b)).

図1に戻り、照度差ステレオ形状取得部(表面形状取得部)143は、上記のように、コード化照明を、被写体のサイズや表面の凹凸にスケールに応じて、動的に変化させて、偏光撮影装置(偏光RGBカメラ)40が撮影した情報を取得する。
その際、照度差ステレオ形状取得部143は、複数の異なる光源方向から撮影した被写体表面の輝度に基づき法線を推定することにより面の傾きを求め、その被写体の詳細な形状情報を被写体表面形状情報として取得する。このとき、照度差ステレオ形状取得部143は、鏡面反射成分Is を除いた拡散反射成分Idのみを用いることにより、被写体表面の凹凸の状態をより正確に反映させて面の傾きを求めることができる。
Returning to FIG. 1 , the photometric stereo shape acquisition unit (surface shape acquisition unit) 143 dynamically changes the coded illumination in accordance with the size of the subject and the scale of the surface irregularities, as described above, and acquires the information captured by the polarization imaging device (polarized RGB camera) 40.
In this case, the photometric stereo shape acquisition unit 143 obtains the inclination of the surface by estimating a normal based on the luminance of the object surface photographed from a plurality of different light source directions, and obtains detailed shape information of the object as object surface shape information. In this case, the photometric stereo shape acquisition unit 143 uses only the diffuse reflection component Id , excluding the specular reflection component Is , so that the photometric stereo shape acquisition unit 143 can obtain the inclination of the surface while more accurately reflecting the uneven state of the object surface.

3次元モデル詳細情報出力部150は、質感計算部140が計算した、拡散反射係数Kd 、鏡面反射係数Ks 、被写体表面の凹凸の状態を示す粗さσ、および、被写体表面の詳細な形状情報(被写体表面形状情報)を含む情報を、各情報を算出するために偏光撮影装置40が撮影した時刻に合わせて3次元モデル詳細情報としてCG描画装置等の外部装置へ出力する。この3次元モデル詳細情報を取得することにより、CG描画装置は、被写体表面の細かな凹凸の状態を質感として、リアルタイムに表現することが可能となる。 The three-dimensional model detailed information output section 150 outputs information calculated by the texture calculation section 140, including the diffuse reflection coefficient K d , the specular reflection coefficient K s , the roughness σ indicating the uneven state of the subject's surface, and detailed shape information of the subject's surface (subject's surface shape information), as three-dimensional model detailed information to an external device such as a CG drawing device in accordance with the time of image capture by the polarization imaging device 40 to calculate each piece of information. By acquiring this three-dimensional model detailed information, the CG drawing device can represent the fine uneven state of the subject's surface as a texture in real time.

<質感取得装置の動作>
次に、本実施形態に係る質感取得装置1が実行する動作について、図17を参照して説明する。
<Operation of the texture acquisition device>
Next, the operation performed by the texture acquisition device 1 according to the present embodiment will be described with reference to FIG.

まず、質感取得装置1の初期設定情報取得部110は、撮影ブースBHに設置した、各装置(奥行き計測装置20、偏光照明装置30および偏光撮影装置40)の、カメラ位置、カメラ姿勢、照明位置、配光特性等の情報を、初期設定情報200として外部装置等から取得し、記憶部12に記憶しておく(ステップS1)。 First, the initial setting information acquisition unit 110 of the texture acquisition device 1 acquires information such as the camera position, camera attitude, lighting position, and light distribution characteristics of each device (depth measurement device 20, polarized lighting device 30, and polarized shooting device 40) installed in the shooting booth BH from an external device, etc., as initial setting information 200, and stores it in the memory unit 12 (step S1).

次に、質感取得装置1の被写体形状取得部120は、計測制御部130を介して、撮影ブースBHに配置された奥行き計測装置20に計測を実行させる指示情報を送信する。そして、被写体形状取得部120は、奥行き計測装置20からの計測情報を受信することにより、(大まかな)被写体形状を取得する(ステップS2)。 Next, the subject shape acquisition unit 120 of the texture acquisition device 1 transmits instruction information to the depth measurement device 20 arranged in the photography booth BH via the measurement control unit 130 to perform measurement. The subject shape acquisition unit 120 then receives measurement information from the depth measurement device 20 to acquire a (rough) subject shape (step S2).

続いて、質感取得装置1の質感計算部140(凹凸状態算出部141)は、被写体表面の凹凸の状態を示す粗さσを算出する(ステップS3)。
ここでは、凹凸状態算出部141が、注目する画素近傍の色ベクトルの分散状態から主成分分析の手法を用いて粗さσを算出するものとする。
具体的には、凹凸状態算出部141は、図6に示すように、RGBの偏向光を計測制御部130を介して偏光照明装置30に照射させ、偏光RGBカメラ40aに撮影させる。そして、凹凸状態算出部141は、RGB空間上の鏡面反射光色ベクトルに対し主成分分析を適用し、第1主成分寄与率と第3主成分寄与率とを計算する。凹凸状態算出部141は、計算した第1主成分寄与率と第3主成分寄与率を用いて、図8の式(1)に基づき、粗さσを算出する。
Next, the texture calculation unit 140 (unevenness calculation unit 141) of the texture acquisition device 1 calculates the roughness σ that indicates the unevenness of the subject surface (step S3).
Here, the unevenness calculation unit 141 calculates the roughness σ from the variance of color vectors in the vicinity of a pixel of interest using a principal component analysis technique.
Specifically, as shown in Fig. 6, the unevenness calculation unit 141 causes the polarized lighting device 30 to irradiate RGB polarized light via the measurement control unit 130, and causes the polarized RGB camera 40a to capture the image. Then, the unevenness calculation unit 141 applies principal component analysis to the specular reflected light color vector in the RGB space, and calculates the first principal component contribution rate and the third principal component contribution rate. The unevenness calculation unit 141 calculates the roughness σ based on the formula (1) in Fig. 8 using the calculated first principal component contribution rate and third principal component contribution rate.

次に、質感取得装置1の質感計算部(反射係数算出部142)は、偏向光の被写体表面での反射光について、拡散反射成分Idと鏡面反射成分Isの分離を行う(ステップS4)。
具体的には、反射係数算出部142は、図10に示すように、偏光方向0度,45度,90度,135度の画素で得られた輝度値から、正弦波を生成する。そして、反射係数算出部142は、その正弦波の輝度Iの最小値Iminを拡散反射成分Idとし、振幅aの2倍である2aを鏡面反射成分Isとする。
この際、反射係数算出部142は、計測制御部130を介し、偏光照明装置30のコード化照明を、例えば、図13に示すように、RGBの3色について、色、偏光方向、および、スケールのパターンをフレーム単位で動的に変化させて被写体に照射する。
Next, the texture calculation unit (reflection coefficient calculation unit 142) of the texture acquisition device 1 separates the polarized light reflected from the subject surface into a diffuse reflection component Id and a specular reflection component Is (step S4).
Specifically, the reflection coefficient calculation unit 142 generates a sine wave from the luminance values obtained from pixels with polarization directions of 0 degrees, 45 degrees, 90 degrees, and 135 degrees, as shown in Fig. 10. Then, the reflection coefficient calculation unit 142 determines the minimum value Imin of the luminance I of the sine wave as the diffuse reflection component Id , and determines 2a, which is twice the amplitude a, as the specular reflection component Is .
At this time, the reflection coefficient calculation unit 142, via the measurement control unit 130, illuminates the subject with coded illumination from the polarized illumination device 30 by dynamically changing the color, polarization direction, and scale pattern for the three colors RGB on a frame-by-frame basis, for example, as shown in FIG. 13 .

続いて、反射係数算出部142は、拡散反射係数Kdを算出する(ステップS5)。
ここで、反射係数算出部142は、観測した角度θ、入射光量I0、拡散反射成分Idを、ランバートの余弦則に基づく上記した式(3)に挿入し、拡散反射係数Kdを算出する。
Next, the reflection coefficient calculation unit 142 calculates the diffuse reflection coefficient K d (step S5).
Here, the reflection coefficient calculation unit 142 inserts the observed angle θ, the amount of incident light I 0 , and the diffuse reflection component I d into the above-mentioned equation (3) based on Lambert's cosine law to calculate the diffuse reflection coefficient K d .

次に、反射係数算出部142は、鏡面反射係数Ksを算出する(ステップS6)。
反射係数算出部142は、Phong反射モデルに基づく上記した式(4)において、粗さσとべき数nとの関係を、上記した〔手法1〕(手法1-1~手法1-3)や〔手法2〕(手法2-1,手法2-2)を用いて算出する。
Next, the reflection coefficient calculation unit 142 calculates the specular reflection coefficient Ks (step S6).
The reflection coefficient calculation unit 142 calculates the relationship between the roughness σ and the exponent n in the above-mentioned equation (4) based on the Phong reflection model, using the above-mentioned [Method 1] (Methods 1-1 to 1-3) or [Method 2] (Method 2-1, Method 2-2).

そして、照度差ステレオ形状取得部(表面形状取得部)143は、計測制御部130を介して、偏光照明装置30のコード化照明を、例えば図14に示すように、動的に変化させて、偏光撮影装置(偏光RGBカメラ)40が撮影した情報を取得する。これにより、照度差ステレオ形状取得部143は、照度差ステレオの手法を用いて、複数の異なる光源方向から撮影した被写体表面の拡散反射成分についての輝度に基づき法線を推定することにより面の傾きを求め、その被写体表面の詳細な形状情報を被写体表面形状情報として取得する(ステップS7)。 The photometric stereo shape acquisition unit (surface shape acquisition unit) 143 dynamically changes the coded illumination of the polarized illumination device 30 via the measurement control unit 130, for example as shown in FIG. 14, and acquires information captured by the polarized imaging device (polarized RGB camera) 40. As a result, the photometric stereo shape acquisition unit 143 uses the photometric stereo technique to estimate the normal based on the luminance of the diffuse reflection components of the subject surface captured from multiple different light source directions, thereby determining the inclination of the surface, and acquires detailed shape information of the subject surface as subject surface shape information (step S7).

次に、3次元モデル詳細情報出力部150は、質感計算部140が生成した、拡散反射係数Kd 、鏡面反射係数Ks 、被写体表面の凹凸の状態を示す粗さσ、および、被写体表面の詳細な形状情報(被写体表面形状情報)を含む情報を、各情報を算出するために偏光撮影装置40が撮影した時刻に合わせて3次元モデル詳細情報としてCG描画装置等の外部装置へ出力し(ステップS8)、処理を終える。 Next, the 3D model detailed information output unit 150 outputs the information generated by the texture calculation unit 140, including the diffuse reflection coefficient K d , the specular reflection coefficient K s , the roughness σ indicating the uneven state of the subject's surface, and detailed shape information of the subject's surface (subject's surface shape information), to an external device such as a CG drawing device as 3D model detailed information in accordance with the time when the polarization imaging device 40 captured the image in order to calculate each piece of information (step S8), and then the processing is completed.

このようにすることにより、本実施形態に係る質感取得装置1を含む質感取得システム1000は、被写体の質感を表現するための被写体表面の細かな凹凸の状態について、正確な被写体表面の拡散反射係数や鏡面反射係数、被写体表面の凹凸の状態(粗さσ)、および、照度差ステレオに基づく被写体表面の詳細な形状情報(被写体表面形状情報)をリアルタイムに取得することが可能となる。 In this way, the texture acquisition system 1000 including the texture acquisition device 1 according to this embodiment can acquire in real time the accurate diffuse reflection coefficient and specular reflection coefficient of the subject surface, the state of the unevenness of the subject surface (roughness σ), and detailed shape information of the subject surface based on photometric stereo (subject surface shape information) for the fine unevenness of the subject surface to express the texture of the subject.

1 質感取得装置
10 制御部
11 入出力部
12 記憶部
20 奥行き計測装置
30 偏光照明装置
40 偏光撮影装置
40a 偏光RGBカメラ
110 初期設定情報取得部
120 被写体形状取得部
130 計測制御部
140 質感計算部
141 凹凸状態算出部
142 反射係数算出部
143 照度差ステレオ形状取得部(表面形状取得部)
150 3次元モデル詳細情報出力部
200 初期設定情報
1000 質感取得システム
BH 撮影ブース
REFERENCE SIGNS LIST 1 Material quality acquisition device 10 Control unit 11 Input/output unit 12 Memory unit 20 Depth measurement device 30 Polarized lighting device 40 Polarized photography device 40a Polarized RGB camera 110 Initial setting information acquisition unit 120 Subject shape acquisition unit 130 Measurement control unit 140 Material quality calculation unit 141 Roughness state calculation unit 142 Reflection coefficient calculation unit 143 Photometric stereo shape acquisition unit (surface shape acquisition unit)
150 3D model detailed information output unit 200 Initial setting information 1000 Texture acquisition system BH Photo booth

Claims (5)

被写体を取り囲んで配置される、奥行き計測装置、偏光照明装置および偏光撮影装置を有する撮影ブースと、
前記奥行き計測装置、前記偏光照明装置および前記偏光撮影装置に通信接続され、前記被写体表面の凹凸の状態で示される質感を取得する質感取得装置と、を備える質感取得システムであって、
前記質感取得装置は、
前記偏光照明装置の照明位置、前記奥行き計測装置および前記偏光撮影装置のカメラ位置およびカメラ姿勢を含む初期設定情報が記憶される記憶部と、
前記奥行き計測装置が計測した奥行きの情報を受信することにより、前記被写体の形状情報を取得する被写体形状取得部と、
複数色の偏向光を前記偏光照明装置が前記被写体に照射する制御を行う計測制御部と、
前記複数色の偏向光が照射された前記被写体からの反射光を前記偏光撮影装置が撮影した情報を取得し、前記被写体表面の凹凸により生じた前記偏向光の色ベクトルの分散から被写体表面の粗さを算出する凹凸状態算出部と、
前記反射光を拡散反射成分と鏡面反射成分とに分離し、前記拡散反射成分、前記被写体への入射光量、前記被写体形状取得部が取得した前記被写体の形状情報および前記照明位置を用いて、前記被写体表面の拡散反射係数を算出するともに、前記鏡面反射成分、前記被写体への入射光量、前記凹凸状態算出部が算出した前記被写体表面の粗さおよび前記カメラ位置を用いて、前記被写体表面の鏡面反射係数を算出する反射係数算出部と、
前記計測制御部を介して複数の異なる光源方向から前記複数色の偏向光を前記被写体に照射させ、前記偏光撮影装置が撮影した被写体表面の拡散反射成分についての輝度に基づき法線を推定することにより前記被写体表面の形状情報を被写体表面形状情報として取得する表面形状取得部と、
前記被写体表面の粗さ、前記拡散反射係数、前記鏡面反射係数、および、前記被写体表面形状情報を含む情報を、前記被写体表面の粗さ、前記拡散反射係数、前記鏡面反射係数、および、前記被写体表面形状情報を含む情報を算出するために前記偏光撮影装置が撮影した時刻に合わせて3次元モデル詳細情報として出力する3次元モデル詳細情報出力部と、
を備えることを特徴とする質感取得システム。
a photography booth having a depth measurement device, a polarized lighting device, and a polarized photography device, the photography booth being arranged to surround a subject;
a texture acquisition device that is communicatively connected to the depth measurement device, the polarized illumination device, and the polarized imaging device, and that acquires a texture represented by a state of unevenness on the surface of the subject,
The texture acquisition device includes:
a storage unit for storing initial setting information including an illumination position of the polarized illumination device, and camera positions and camera orientations of the depth measurement device and the polarized imaging device;
a subject shape acquisition unit that acquires shape information of the subject by receiving depth information measured by the depth measurement device;
a measurement control unit that controls the polarized illumination device to irradiate the subject with polarized light of a plurality of colors;
an unevenness calculation unit that acquires information obtained by capturing, by the polarization imaging device, reflected light from the object irradiated with the polarized light of the plurality of colors, and calculates the roughness of the object surface from the variance of the color vector of the polarized light generated by the unevenness of the object surface;
a reflection coefficient calculation unit that separates the reflected light into a diffuse reflection component and a specular reflection component, calculates a diffuse reflection coefficient of the subject surface using the diffuse reflection component, the amount of light incident on the subject, the shape information of the subject acquired by the subject shape acquisition unit, and the illumination position, and calculates a specular reflection coefficient of the subject surface using the specular reflection component, the amount of light incident on the subject, the roughness of the subject surface calculated by the unevenness calculation unit, and the camera position;
a surface shape acquisition unit that irradiates the subject with the polarized light of the multiple colors from multiple different light source directions via the measurement control unit, and acquires shape information of the subject surface as subject surface shape information by estimating a normal based on the luminance of a diffuse reflection component of the subject surface photographed by the polarization imaging device;
a three-dimensional model detailed information output unit that outputs information including the roughness of the object surface, the diffuse reflection coefficient, the specular reflection coefficient, and the object surface shape information as three-dimensional model detailed information in accordance with the time when the polarization imaging device captures the image in order to calculate information including the roughness of the object surface, the diffuse reflection coefficient, the specular reflection coefficient, and the object surface shape information;
A texture acquisition system comprising:
撮影ブースにおいて被写体を取り囲んで配置される、奥行き計測装置、偏光照明装置および偏光撮影装置に通信接続され、前記被写体表面の凹凸の状態で示される質感を取得する質感取得装置であって、
前記偏光照明装置の照明位置、前記奥行き計測装置および前記偏光撮影装置のカメラ位置およびカメラ姿勢を含む初期設定情報が記憶される記憶部と、
前記奥行き計測装置が計測した奥行きの情報を受信することにより、前記被写体の形状情報を取得する被写体形状取得部と、
複数色の偏向光を前記偏光照明装置が前記被写体に照射する制御を行う計測制御部と、
前記複数色の偏向光が照射された前記被写体からの反射光を前記偏光撮影装置が撮影した情報を取得し、前記被写体表面の凹凸により生じた前記偏向光の色ベクトルの分散から被写体表面の粗さを算出する凹凸状態算出部と、
前記反射光を拡散反射成分と鏡面反射成分とに分離し、前記拡散反射成分、前記被写体への入射光量、前記被写体形状取得部が取得した前記被写体の形状情報および前記照明位置を用いて、前記被写体表面の拡散反射係数を算出するともに、前記鏡面反射成分、前記被写体への入射光量、前記凹凸状態算出部が算出した前記被写体表面の粗さおよび前記カメラ位置を用いて、前記被写体表面の鏡面反射係数を算出する反射係数算出部と、
前記計測制御部を介して複数の異なる光源方向から前記複数色の偏向光を前記被写体に照射させ、前記偏光撮影装置が撮影した被写体表面の拡散反射成分についての輝度に基づき法線を推定することにより前記被写体表面の形状情報を被写体表面形状情報として取得する表面形状取得部と、
前記被写体表面の粗さ、前記拡散反射係数、前記鏡面反射係数、および、前記被写体表面形状情報を含む情報を、前記被写体表面の粗さ、前記拡散反射係数、前記鏡面反射係数、および、前記被写体表面形状情報を含む情報を算出するために前記偏光撮影装置が撮影した時刻に合わせて3次元モデル詳細情報として出力する3次元モデル詳細情報出力部と、
を備えることを特徴とする質感取得装置。
A texture acquisition device that is communicatively connected to a depth measurement device, a polarized lighting device, and a polarized photography device that are arranged to surround a subject in a photography booth, and that acquires a texture represented by a state of unevenness on the surface of the subject,
a storage unit for storing initial setting information including an illumination position of the polarized illumination device, and camera positions and camera orientations of the depth measurement device and the polarized imaging device;
a subject shape acquisition unit that acquires shape information of the subject by receiving depth information measured by the depth measurement device;
a measurement control unit that controls the polarized illumination device to irradiate the subject with polarized light of a plurality of colors;
an unevenness calculation unit that acquires information obtained by capturing, by the polarization imaging device, reflected light from the object irradiated with the polarized light of the plurality of colors, and calculates the roughness of the object surface from the variance of the color vector of the polarized light generated by the unevenness of the object surface;
a reflection coefficient calculation unit that separates the reflected light into a diffuse reflection component and a specular reflection component, calculates a diffuse reflection coefficient of the subject surface using the diffuse reflection component, the amount of light incident on the subject, the shape information of the subject acquired by the subject shape acquisition unit, and the illumination position, and calculates a specular reflection coefficient of the subject surface using the specular reflection component, the amount of light incident on the subject, the roughness of the subject surface calculated by the unevenness calculation unit, and the camera position;
a surface shape acquisition unit that irradiates the subject with the polarized light of the multiple colors from multiple different light source directions via the measurement control unit, and acquires shape information of the subject surface as subject surface shape information by estimating a normal based on the luminance of a diffuse reflection component of the subject surface photographed by the polarization imaging device;
a three-dimensional model detailed information output unit that outputs information including the roughness of the object surface, the diffuse reflection coefficient, the specular reflection coefficient, and the object surface shape information as three-dimensional model detailed information in accordance with the time when the polarization imaging device captures the image in order to calculate information including the roughness of the object surface, the diffuse reflection coefficient, the specular reflection coefficient, and the object surface shape information;
A texture acquisition device comprising:
前記凹凸状態算出部は、
前記偏向光の色ベクトルに対し主成分分析を行い、算出した第1主成分寄与率と第3主成分寄与率とを用いて前記被写体表面の粗さの算出すること
を特徴とする請求項2に記載の質感取得装置。
The unevenness state calculation unit
3. The texture acquisition device according to claim 2, further comprising: a principal component analysis being performed on the color vector of the polarized light; and a roughness of the surface of the object being calculated using a first principal component contribution rate and a third principal component contribution rate.
前記計測制御部は、
前記偏光照明装置が前記被写体に照射する前記複数色の偏向光について、光源を面光源とし、色、偏光方向およびスケールをフレーム単位で変化させて、前記被写体に照射させること
を特徴とする請求項2または請求項3に記載の質感取得装置。
The measurement control unit is
4. The texture acquisition device according to claim 2 or 3, wherein the polarized lighting device irradiates the subject with the polarized light of the plurality of colors, the polarized light being a surface light source, and the color, polarization direction, and scale of the polarized light being changed on a frame-by-frame basis.
コンピュータを、請求項2乃至請求項4のいずれか1項に記載の質感取得装置として機能させるための質感取得プログラム。 A texture acquisition program for causing a computer to function as the texture acquisition device according to any one of claims 2 to 4.
JP2020168205A 2020-10-05 2020-10-05 Texture acquisition system, texture acquisition device, and texture acquisition program Active JP7633790B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2020168205A JP7633790B2 (en) 2020-10-05 2020-10-05 Texture acquisition system, texture acquisition device, and texture acquisition program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020168205A JP7633790B2 (en) 2020-10-05 2020-10-05 Texture acquisition system, texture acquisition device, and texture acquisition program

Publications (2)

Publication Number Publication Date
JP2022060641A JP2022060641A (en) 2022-04-15
JP7633790B2 true JP7633790B2 (en) 2025-02-20

Family

ID=81125215

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020168205A Active JP7633790B2 (en) 2020-10-05 2020-10-05 Texture acquisition system, texture acquisition device, and texture acquisition program

Country Status (1)

Country Link
JP (1) JP7633790B2 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024176831A1 (en) * 2023-02-22 2024-08-29 ソニーグループ株式会社 Image processing apparatus and image processing system
CN117523112B (en) * 2024-01-05 2024-05-17 深圳市宗匠科技有限公司 Three-dimensional model building method and system, control equipment and storage medium thereof
WO2025205731A1 (en) * 2024-03-28 2025-10-02 富士フイルム株式会社 Three-dimensional model generation system, photographing condition setting device, and photographing system
JP2026036012A (en) * 2024-08-20 2026-03-05 富士フイルムビジネスイノベーション株式会社 Image capturing system and program

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012039086A1 (en) 2010-09-24 2012-03-29 パナソニック株式会社 Image processing device
JP2012141758A (en) 2010-12-28 2012-07-26 Toshiba Corp Three-dimensional data processing device, method and program
JP2012173916A (en) 2011-02-21 2012-09-10 Nippon Hoso Kyokai <Nhk> Imaging device and image processing information generating program

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012039086A1 (en) 2010-09-24 2012-03-29 パナソニック株式会社 Image processing device
JP2012141758A (en) 2010-12-28 2012-07-26 Toshiba Corp Three-dimensional data processing device, method and program
JP2012173916A (en) 2011-02-21 2012-09-10 Nippon Hoso Kyokai <Nhk> Imaging device and image processing information generating program

Also Published As

Publication number Publication date
JP2022060641A (en) 2022-04-15

Similar Documents

Publication Publication Date Title
TWI870399B (en) System and method for creation of topical agents with improved image capture
JP7633790B2 (en) Texture acquisition system, texture acquisition device, and texture acquisition program
US11210839B2 (en) Photometric image processing
US6639594B2 (en) View-dependent image synthesis
CN110148204B (en) Method and system for representing virtual objects in a view of a real environment
US8264490B2 (en) Practical modeling and acquisition of layered facial reflectance
US11989971B2 (en) Techniques for facial appearance capture
Pitard et al. Discrete modal decomposition: a new approach for the reflectance modeling and rendering of real surfaces
CN109785423A (en) Image light compensation method, device and computer equipment
EP3382645A2 (en) Method for generation of a 3d model based on structure from motion and photometric stereo of 2d sparse images
US20230043787A1 (en) Lighting assembly for producing realistic photo images
JP2023553259A (en) dark flash normal camera
Ciortan et al. A practical reflectance transformation imaging pipeline for surface characterization in cultural heritage
US12482179B2 (en) Freestyle acquisition method for high-dimensional material
Ono et al. Practical BRDF reconstruction using reliable geometric regions from multi-view stereo
Li et al. Dpcs: Path tracing-based differentiable projector-camera systems
KR20240169499A (en) Learning method of determining light patterns for three-dimensional scenes, method and apparatus of three-dimensional scenes
JP7681470B2 (en) Texture acquisition system, texture acquisition device, and texture acquisition program
Thul et al. Precomputed radiance transfer for reflectance and lighting estimation
Sitnik et al. Processing paths from integrating multimodal 3D measurement of shape, color and BRDF
Boll Nielsen et al. A variational study on BRDF reconstruction in a structured light scanner
JP2024068744A (en) Subsurface scattering analysis device, subsurface scattering analysis system, and subsurface scattering analysis program
Coutinho et al. Assisted color acquisition for 3d models
Feng et al. Learning Photometric Feature Transform for Free-Form Object Scan
Padhye Techniques to Create a Realistic Digital Collection of Near Planar Cultural Heritage Objects

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230904

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20241112

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20241216

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20250114

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20250207

R150 Certificate of patent or registration of utility model

Ref document number: 7633790

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150