[go: up one dir, main page]

JP7609075B2 - 情報処理装置および情報処理方法 - Google Patents

情報処理装置および情報処理方法 Download PDF

Info

Publication number
JP7609075B2
JP7609075B2 JP2021556079A JP2021556079A JP7609075B2 JP 7609075 B2 JP7609075 B2 JP 7609075B2 JP 2021556079 A JP2021556079 A JP 2021556079A JP 2021556079 A JP2021556079 A JP 2021556079A JP 7609075 B2 JP7609075 B2 JP 7609075B2
Authority
JP
Japan
Prior art keywords
image
information processing
processing device
captured image
information
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
JP2021556079A
Other languages
English (en)
Other versions
JPWO2021095672A1 (ja
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.)
Sony Corp
Sony Group Corp
Original Assignee
Sony Corp
Sony Group 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 Sony Corp, Sony Group Corp filed Critical Sony Corp
Publication of JPWO2021095672A1 publication Critical patent/JPWO2021095672A1/ja
Application granted granted Critical
Publication of JP7609075B2 publication Critical patent/JP7609075B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/107Measuring physical dimensions, e.g. size of the entire body or parts thereof
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/30Measuring arrangements characterised by the use of optical techniques for measuring roughness or irregularity of surfaces
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/27Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands using photo-electric detection ; circuits for computing concentration
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/10Image acquisition
    • G06V10/12Details of acquisition arrangements; Constructional details thereof
    • G06V10/14Optical characteristics of the device performing the acquisition or on the illumination arrangements
    • G06V10/141Control of illumination
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/25Determination of region of interest [ROI] or a volume of interest [VOI]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/60Extraction of image or video features relating to illumination properties, e.g. using a reflectance or lighting model
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/74Image or video pattern matching; Proximity measures in feature spaces
    • G06V10/761Proximity, similarity or dissimilarity measures
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/82Arrangements for image or video recognition or understanding using pattern recognition or machine learning using neural networks
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/60Type of objects
    • G06V20/69Microscopic objects, e.g. biological cells or cellular parts
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10056Microscopic image
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V2201/00Indexing scheme relating to image or video recognition or understanding
    • G06V2201/03Recognition of patterns in medical or anatomical images

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Multimedia (AREA)
  • General Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Software Systems (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Artificial Intelligence (AREA)
  • Computing Systems (AREA)
  • Databases & Information Systems (AREA)
  • Molecular Biology (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Mathematical Physics (AREA)
  • Chemical & Material Sciences (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Image Analysis (AREA)

Description

本開示は、情報処理装置および情報処理方法に関する。
物体の微細な状態を観測する装置として、比較的導入コストが低く、また、簡易に計測を行うことができるマイクロスコープが広く利用されている。
マイクロスコープに関する技術として、照明ユニットの入射角度の相違を利用して、皮膚表面の色や汚れを解析する技術が知られている(例えば、特許文献1)。また、マイクロスコープの先端ドームからの所定距離に透明ガラスを配置することで、肌表面撮影時のピントボケやゆがみを低減させる技術が知られている(例えば、特許文献2)。
特開平10-333057号公報 特開2008-253498号公報
従来技術によれば、マイクロスコープによって撮影される画像の質を向上させることができる。
しかしながら、従来技術は、あくまで平面的な画像の質を向上させるものであり、物体の微小な形状(凹凸)を再現した3D画像を得ることは難しい。なお、物体の微小な形状を計測する装置として、非接触3D測定器や3Dスキャナ等が利用されるが、これらを導入する場合、コストが比較的高くなるという問題がある。なお、ToF(Time of Flight)方式による測距装置は、コストが比較的低廉であるものの、精度が充分でない場合がある。
そこで、本開示では、形状計測の精度を向上させることができる情報処理装置および情報処理方法を提供する。
本開示によれば、情報処理装置が提供される。情報処理装置は、制御部を備える。制御部は、対象をセンサが撮像した撮像画像を取得する。前記撮像画像は、それぞれ異なる位置に配置される複数の光源から前記対象に照射された光の反射光から得られる画像である。制御部は、前記撮像画像の輝度値に基づき、前記撮像画像から平坦領域を抽出する。制御部は、前記センサに関する情報および前記撮像画像の前記平坦領域に基づき、前記対象の表面の形状に関する形状情報を算出する。
撮像画像から深度を算出する算出方法について説明するための図である。 対象表面の大まかな形状について説明するための表である。 本開示の第1実施形態の概要について説明するための図である。 本開示の第1実施形態に係る情報処理システムの構成例を示すブロック図である。 本開示の第1実施形態に係るマイクロスコープの構成例を示す図である。 本開示の第1実施形態に係る情報処理装置の構成例を示すブロック図である。 撮像画像の輝度値の一例を示す図である。 領域取得部が生成するマスク画像の一例を示す図である。 領域取得部による平坦領域の取得例を説明するための図である。 法線情報の一例について説明するための図である。 本開示の実施形態に係る学習器について説明するための図である。 撮像画像の1pixelあたりの長さについて説明するための図である。 表示制御部が表示部に表示させる画像の一例を示す図である。 本開示の第1実施形態に係る深度算出処理の一例を説明するためのフローチャートである。 本開示の第2実施形態に係る取得部が取得する撮像画像について説明するための図である。 本開示の第2実施形態に係る取得部が取得する撮像画像について説明するための図である。 本開示の第3実施形態に係る領域取得部によるスムージング処理を説明するための図である。 本開示の第4実施形態に係る領域取得部による平坦領域の取得について説明するための図である。 本開示の第4実施形態に係る領域取得部による平坦領域の取得について説明するための図である。 センサの被写界深度を説明するための図である。 本開示の第5実施形態に係る領域取得部による平坦領域の取得について説明するための図である。 本開示の第5実施形態に係る領域取得部による平坦領域の取得について説明するための図である。 本開示の第6実施形態に係る取得部が取得する複数の撮像画像について説明するための図である。 本開示の第6実施形態に係る領域取得部による平坦領域の取得について説明するための図である。 本開示の第7実施形態に係る情報処理装置の構成例を示す図である。 本開示の第7実施形態に係る法線周波数分離部による周波数分離を説明するための図である。 本開示の第8実施形態に係る情報処理装置の構成例を示す図である。 本開示の第8実施形態に係る法線周波数分離部による周波数分離を説明するための図である。 第9実施形態に係る学習器の学習方法について説明するための図である。 本開示の第10実施形態に係る情報処理装置の制御部の構成例を示す図である。 本開示に係る技術(本技術)が適用され得る、カプセル型内視鏡を用いた患者の体内情報取得システムの概略的な構成の一例を示すブロック図である。 本開示に係る技術(本技術)が適用され得る内視鏡手術システムの概略的な構成の一例を示す図である。 図32に示すカメラヘッドおよびCCUの機能構成の一例を示すブロック図である。
以下に添付図面を参照しながら、本開示の好適な実施の形態について詳細に説明する。なお、本明細書および図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。
なお、説明は以下の順序で行うものとする。
1.背景
1.1.表面形状の算出方法
1.2.算出方法の課題
2.第1実施形態
2.1.第1実施形態の概要
2.2.システム構成例
2.3.マイクロスコープの構成例
2.4.情報処理装置の構成例
2.5.深度算出処理
3.第2実施形態
4.第3実施形態
5.第4実施形態
6.第5実施形態
7.第6実施形態
8.第7実施形態
9.第8実施形態
10.第9実施形態
11.第10実施形態
12.その他の実施形態
13.適応例
14.補足
<1.背景>
<1.1.表面形状の算出方法>
まず、本開示の実施形態の詳細を説明する前に、本発明者らが本開示の実施形態を創作するに至った背景について説明する。
撮像装置によって撮影される撮像画像(RGB画像)から物体の形状(凹凸)を再現した3D画像を得る方法として、例えばCNN(Convolutional Neural Network)等の機械学習を利用してRGB画像から直接表面形状(表面までの深度)を算出する方法が知られている。しかしながら、RGB画像からCNNで深度を直接算出する場合、不確定性があり精度を向上させることが難しかった。
撮像画像から物体表面までの深度を算出する他の方法として、撮像画像からCNN等の機械学習を利用して物体表面の法線情報を算出し、かかる法線情報を深度に式変換する方法がある。かかる方法について、図1を用いて説明する。図1は、撮像画像から深度を算出する算出方法について説明するための図である。図1を用いて説明する方法は、例えば情報処理装置200によって実行される。
かかる方法では、情報処理装置200は、まずマイクロスコープが撮像した撮像画像を取得する(ステップS1)。このとき、情報処理装置200は、例えばマイクロスコープに関する情報も取得するものとする。情報処理装置200は、撮像画像に基づいて、対象の表面における法線(法線情報)を算出する(ステップS2)。情報処理装置200は、例えば、CNN等を利用して学習された学習器(モデル)に撮像画像を入力することで、出力データとして法線情報を得る。
情報処理装置200は、法線情報に基づき、対象の表面までの距離(深度情報)を算出する(ステップS3)。マイクロスコープに関する情報の1つであるマイクロスコープのセンサと対象との距離が既知であれば、情報処理装置200は、法線情報に基づいて対象の表面までの測距を行える。ここで、ユーザがマイクロスコープで対象の撮像を行う場合、マイクロスコープのヘッドマウント部を対象に接触させて撮像する。そのため、マイクロスコープのセンサと対象との距離は、ヘッドマウント部の長さに相当する既知の値となる。
具体的に、情報処理装置200は、以下の式(1)のWを最小化することで対象表面までの深度を法線情報に基づいて算出する。
Figure 0007609075000001
なお、式(1)の各パラメータは以下の通りである。
p:算出した法線のx方向
q:算出した法線のy方向
:求める深度のx方向の偏微分(法線のx方向)
:求める深度のy方向の偏微分(法線のy方向)
xx:求める深度のx方向の2回偏微分(法線のx方向のx方向への偏微分)
yy:求める深度のy方向の2回偏微分(法線のy方向のy方向への偏微分)
xy:求める深度のx、y方向の2回偏微分
なお、x、yは、撮像画像上の座標を示しており、x方向は、例えば撮像画像の横方向、y方向は撮像画像の縦方向とする。
式(1)のフーリエ変換(周波数変換)を行うと式(2)になる。これを最小化することで、対象表面までの深度が算出される。なお、式(2)のZは求める深度をフーリエ変換したものであり、u、vは、周波数空間での座標を表している。
Figure 0007609075000002
式(2)を式展開したものが、以下の式(3)であり、式(3)を逆フーリエ変換することで、対象表面までの深度が算出される。
Figure 0007609075000003
上述した式(1)~式(3)の式変換を用いて、情報処理装置200は、対象の表面までの距離(深度)を算出する。
<1.2.算出方法の課題>
上述した算出方法において、式(1)の右辺の2段目および3段目は、式(1)におけるcost項であり、2段目が法線の絶対値の総和、3段目が法線の微分の絶対値の総和を表している。
ここで、上述した算出方法では、cost項の重み(λ、μ)が小さい、すなわち法線の絶対値の総和が少なく、法線の微分の絶対値の総和が少ないという仮定で深度を算出している。法線の絶対値の総和が少ないという仮定(以下、仮定1とも称する)は、深度を算出する対象表面が平らであることを意味する。また、法線の微分の絶対値の総和が少ないという仮定(以下、仮定2とも称する)は、深度を算出する対象表面の曲率が小さいことを意味する。すなわち、上述した算出方法では、深度を算出する対象の大まかな表面形状が平面であるという仮定で、対象までの深度を算出している。
ここで、図2を用いて、対象表面の大まかな形状について説明する。図2は、対象表面の大まかな形状について説明するための表である。なお、図2では、実際の対象表面の形状を点線で、対象表面の大まかな形状について直線で示している。また、対象表面の法線を矢印で示している。また、図2では、対象の横方向をx方向、縦方向をz方向(深度)としている。
対象表面の大まかな形状とは、対象表面の撮像画像全体における形状であり、大まかな形状が平面であるとは、対象表面の撮像画像全体における形状を見たときに変動が少ないことを意味する。なお、対象表面の局所的な凹凸の変化は、算出対象であるため大まかな形状には含まれない。上述した式(1)に示す最適化式のcost項の仮定を満たす、すなわち対象表面の撮像画像全体における表面形状が平面に近いほど、情報処理装置200は、対象表面までの深度を高精度に算出することができる。
例えば、図2の表の(1)に示す対象は、まっすぐな表面に細かな凹凸が含まれる形状となっている。そのため、(1)に示す対象は、上述した仮定1、仮定2の両方を満たす(表では「○」で示している)ため、かかる対象の深度を情報処理装置200は高精度に算出することができる。
図2の表の(2)に示す対象は、わずかに上向きにカーブした表面に細かな凹凸が含まれる形状となっている。そのため、(2)に示す対象は、上述した仮定1をやや満たし(表では「△」で示している)、仮定2を満たす。そのため、情報処理装置200は、(2)に示す対象の深度を(1)よりは精度が下がるが、精度よく算出することができる。
図2の表の(3)に示す対象は、左側が大きく上向きにカーブした表面に細かな凹凸が含まれる形状となっている。そのため、(3)に示す対象は、上述した仮定1、2ともに満たさない(表では「×」で示している)。そのため、情報処理装置200は、(3)に示す対象の深度を精度よく算出することができない。
図2の表の(4)に示す対象は、表面が段々になっている。このように、(4)に示す対象は、高さや方向が異なるまっすぐな面を複数有する表面に細かな凹凸が含まれる形状になっているため、上述した仮定1、2ともに満たさない。そのため、情報処理装置200は、(4)に示す対象の深度を精度よく算出することができない。
このように、上述した算出方法では、深度の算出対象の表面の大まかな形状が平らであるという仮定を満たさない場合、深度の算出精度が低下するという問題があった。
そこで、このような状況を鑑みて、本件の開示者は、深度の算出精度を向上させることで、形状計測の精度を向上させることができる情報処理装置200に係る本開示の各実施形態を創作するに至った。そこで、以下に、本開示に係る各実施形態の詳細について順次説明する。
<2.第1実施形態>
<2.1.第1実施形態の概要>
図3は、本開示の第1実施形態の概要について説明するための図である。第1実施形態に係る情報処理装置200A(図示省略)は、上述した仮定1、2を満たす領域(平坦領域)を撮像画像から抽出することで、対象の表面形状に関する形状情報(深度情報)の精度劣化を抑制し、形状計測の精度を向上させる。
情報処理装置200Aは、まずマイクロスコープ100のセンサ150が対象Sを撮像した撮像画像M11を取得する。撮像画像M11は、それぞれ異なる位置に配置される複数の光源160A、160Bから対象Sに照射された光IA、IBの反射光から得られる画像である。
ここで、マイクロスコープ100について簡単に説明する。図3の左図に示すように、マイクロスコープ100は、センサ150と、センサ150と撮像対象との間に設置される筒状の機構であるヘッドマウント部10と、複数の光源160A、160Bと、を有する。なお、センサ150は、レンズやカメラ等と読み替えてもよい。マイクロスコープ100は、ユーザが手に持って、対象Sにヘッドマウント部10を接触させることで、対象Sにセンサ150を向けて使用される撮像装置である。
マイクロスコープ100は、光源160A、160Bから対象Sに同時に照射された光IA、IBの反射光を露光させ、対象Sを撮像する。このとき、例えば対象Sの表面の大まかな形状が平面でない場合、対象Sの表面には、オクルージョン(光が当たらない領域)が発生する。
例えば、光源160Aから照射された光IAは、対象Sの表面のうち領域SAおよび領域SABには当たるが、領域SBには当たらない。一方、光源160Bから照射された光IBは、対象Sの表面のうち領域SBおよび領域SABには当たるが、領域SAには当たらない。複数の光源160A、160Bからの光IA、IBが両方当たる領域SABは、オクルージョンがない平坦な領域である。
マイクロスコープ100の撮像画像M11は、図3の右図に示すように、光源160A、160Bの一方の光IA、IBしか当たらない領域SA、SBが、両方の光IA、IBが当たる領域SABに比べて輝度値が低く暗い画像となる。
そこで、情報処理装置200Aは、撮像画像M11の輝度値に基づき、撮像画像M11から平坦領域SABを抽出する。例えば、情報処理装置200Aは、撮像画像M11の輝度値が閾値以上である領域SABを平坦領域として抽出する。
情報処理装置200Aは、センサ150に関する情報および撮像画像M11の平坦領域SABに基づき、対象Sの表面の形状に関する形状情報(深度情報)を算出する。例えば、情報処理装置200Aは、CNNを利用して学習された学習器に、撮像画像M11の平坦領域SABを入力することで平坦領域SABの法線情報を取得する。情報処理装置200Aは、取得した法線情報を上述した式(1)~式(3)を用いた式変換を行うことで、対象Sの表面までの深度情報を算出する。
このように、情報処理装置200Aは、撮像画像M11から式(1)の仮定1、2を満たす領域を抽出することで、形状計測の精度を向上させることができる。
以下、上述した情報処理装置200Aを有する情報処理システム1の詳細について説明する。
<2.2.システム構成例>
図4は、本開示の第1実施形態に係る情報処理システム1の構成例を示すブロック図である。図4に示すように、情報処理システム1は、マイクロスコープ100と、情報処理装置200Aと、を有する。
マイクロスコープ100は、ユーザが手に持って、撮像対象Sにセンサ150を向けて使用される撮像装置である。
情報処理装置200Aは、マイクロスコープ100が撮像した撮像画像M11に基づき、撮像対象Sの表面形状に関する形状情報を算出する。情報処理装置200Aの詳細は、図6を用いて後述する。
マイクロスコープ100と情報処理装置200Aとは例えばケーブルを用いて接続される。あるいは、マイクロスコープ100と情報処理装置200AとがBluetooth(登録商標)やNFC(Near Field Communication)等の無線通信によって直接接続されていてもよい。マイクロスコープ100と情報処理装置200Aとが、例えばネットワーク(図示省略)を介して有線または無線で接続されていてもよい。あるいは、マイクロスコープ100と情報処理装置200Aとが、例えばハードディスク、磁気ディスク、光磁気ディスク、光ディスク、USBメモリ、またはメモリカードなどの外部装着型の記憶媒体を介して撮像画像をやりとりしてもよい。また、マイクロスコープ100と情報処理装置200Aとを一体として構成してもよい。具体的には、情報処理装置200Aが、例えば、マイクロスコープ100の本体内部に配置される構成であってもよい。
<2.3.マイクロスコープの構成例>
図5は、本開示の第1実施形態に係るマイクロスコープ100の構成例を示す図である。マイクロスコープ100は、センサ150と、センサ150と撮像対象との間に設置される筒状の機構であるヘッドマウント部10と、を有する。また、マイクロスコープ100は、それぞれ異なる位置に配置される点光源160A、160B(図示省略)を有する。
ヘッドマウント部10は、マイクロスコープ100の先端にマウントされる機構である。ヘッドマウント部10は、例えば先端ヘッドや鏡筒とも称される。ヘッドマウント部10の内部が、例えば鏡になっており、点光源160A、160Bから照射された光がヘッドマウント部10の側面で全反射するようになっていてもよい。
ユーザは、ヘッドマウント部10を対象Sに接触させて、対象Sを撮像する。これにより、センサ150と対象Sとの距離が固定され、撮像においてピント(焦点距離)がずれることを防止することができる。
マイクロスコープ100は、ヘッドマウント部10の内部に点光源160A、160Bを有する。これにより、マイクロスコープ100は、点光源160A、160Bから対象Sに照射された光IA、IBの反射光を露光させ、対象Sを撮像する。なお、本明細書における点光源とは、理想的には点による光源を意味するが、現実的には点による光源は存在しえないため、きわめて微少な大きさ(数ミリメートル以内など)の光源を含むものである。
なお、本明細書では、光源の一例として点光源を用いた場合について説明するが、光源は点光源に限定されない。また、光源の数も、2つに限定されない。光源は複数あればよく、それぞれことなる位置に配置されていれば光源が3つ以上あってもよい。
<2.4.情報処理装置の構成例>
図6は、本開示の第1実施形態に係る情報処理装置200Aの構成例を示すブロック図である。情報処理装置200Aは、制御部220と、記憶部230と、を有する。
(制御部)
制御部220は、情報処理装置200Aの動作を制御する。制御部220は、取得部221と、領域取得部225と、法線算出部222と、深度算出部223と、表示制御部224と、を有する。取得部221と、領域取得部225と、法線算出部222と、深度算出部223と、表示制御部224と、の各機能部は、例えば、制御部220によって、制御部220内部に記憶されたプログラムがRAM等を作業領域として実行されることにより実現される。なお、制御部220の内部構造は、図6に示した構成に限られず、後述する情報処理を行う構成であれば他の構成であってもよい。また、制御部220が有する各処理部の接続関係は、図6に示した接続関係に限られず、他の接続関係であってもよい。
(取得部)
取得部221は、マイクロスコープ100が撮像した撮像画像M11や、マイクロスコープ100に関する情報を取得する。マイクロスコープ100に関する情報には、例えば焦点距離fやヘッドマウント部10の長さdなど、マイクロスコープ100の構造に関する情報が含まれる。
取得部221が、例えばマイクロスコープ100の点光源160A、160Bやセンサ150を制御するようにしてもよい。この場合、取得部221は、点光源160A、160Bから光IA、IBが同時に照射されるように、点光源160A、160Bを制御する。また、取得部221は、点光源160A、160Bから光IA、IBが同時に照射されている間にセンサ150が対象Sを撮像するように、センサ150を制御する。このように、情報処理装置200Aがマイクロスコープ100を制御するようにしてもよい。
あるいは、取得部221が、撮像画像M11に加え、撮像条件に関する情報を取得するようにしてもよい。撮像条件に関する情報とは、例えば、撮像画像M11が点光源160A、160Bから光IA、IBが同時に照射された間に撮像された画像であることを示す情報である。この場合、後述する領域取得部225は、例えば撮像条件に応じて撮像画像M11から平坦領域SABを抽出する。
(領域取得部)
領域取得部225は、撮像画像M11から平坦領域SABを抽出する。領域取得部225は、例えば、撮像画像M11の各画素の輝度値L(x、y)と閾値thとを比較する。図7に示すように、領域取得部225は、輝度値L(x、y)が閾値th以上の画素を含む領域を処理領域とし、輝度値L(x、y)が閾値th未満の画素を含む領域を除外領域とする。なお、図7は、撮像画像M11の輝度値の一例を示す図であり、図7では、y軸の値を所定値「β」とした場合の輝度値L(x、β)と画素(x、β)との関係を示している。
なお、撮像画像M11がRGB画像の場合、輝度値L(x、y)は、L(x、y)=(R(x、y)+2G(x、y)+B(x、y))/4で算出される。R(x、y)は画素(x、y)における画素値のR(Red)成分、G(x、y)は画素(x、y)における画素値のG(Green)成分、B(x、y)は画素(x、y)における画素値のB(Blue)成分である。
領域取得部225は、全ての画素において輝度値L(x、y)と閾値thとを比較することで、処理領域と判定した領域を撮像画像M11から抽出する抽出領域(平坦領域SAB)として取得する。
領域取得部225は、例えば、画素の輝度値L(x、y)が閾値th以上である場合にかかる画素を白色とし、閾値th未満である場合にかかる画素を黒色とすることで、例えば図8に示すマスク画像M12を生成する。なお、図8は、領域取得部225が生成するマスク画像M12の一例を示す図である。
領域取得部225は、図9に示すように、生成したマスク画像M12と撮像画像M11とを比較し、同座標においてマスク画像M12が白色である画素を撮像画像M11から抽出することで、平坦領域SABを取得する。なお、図9は、領域取得部225による平坦領域SABの取得例を説明するための図である。
あるいは、領域取得部255は、画素の輝度値L(x、y)が閾値th以上である場合にかかる画素の輝度値を「1」とし、閾値th未満である場合にかかる画素の輝度値を「0」とすることで、マスク画像M12を生成してもよい。この場合、領域取得部255は、撮像画像M11とマスク画像M12とを乗算することで、平坦領域SABを取得する。
領域取得部255は、取得した平坦領域SABを法線算出部222に出力する。なお、領域取得部255が生成したマスク画像M12を法線算出部222に出力し、法線算出部222が、マスク画像M12を用いて撮像画像M11から平坦領域SABを取得するようにしてもよい。
(法線算出部)
法線算出部222は、撮像画像M11の平坦領域SABに基づいて、対象Sの表面における法線情報を法線画像として算出する。法線算出部222は、例えば、CNNを用いて法線画像を生成する。法線算出部222は、CNNを利用して事前に学習された学習器を用いて法線画像を生成する。学習器300(図11参照)の入力チャネル数は、入力画像(ここでは平坦領域SAB)のRGBに相当する3チャネル、あるいは、入力画像のグレースケールの1チャネルとする。また、学習器300の出力チャネル数は、出力画像(ここでは法線画像)のRGBに相当する3チャネルとする。また、出力画像の解像度は、入力画像の解像度と等しいものとする。
ここで、法線情報と法線画像のRGBとの関係について、図10を用いて説明する。図10は、法線情報の一例について説明するための図である。図10に示す法線Nは、対象Sを撮像した撮像画像M11の所定画素における法線ベクトルであるとする。この場合、法線Nの天頂角をθ、方位角をφとすると、極座標系における法線N(θ、φ)と直交座標系における法線N(x、y、z)との間には、以下の式(4)の関係が成り立つ。なお、法線Nは、単位ベクトルであるものとする。
Figure 0007609075000004
法線情報は、上述した直交座標系における法線N(x、y、z)であり、例えば撮像画像M11の平坦領域SABの画素ごとに算出される。また、法線画像は、各画素の法線情報をRGBに置き換えることで得られる画像である。すなわち、法線画像は、法線NのxをR(レッド)、yをG(グリーン)、zをB(ブルー)の各値として置き換えることで得られる画像である。
図6に戻る。法線算出部222は、平坦領域SABを入力画像として学習器300に入力することで、出力画像である法線画像を生成する。学習器300は、学習の真値を法線として、CNN中の重みの学習を行うことで事前に生成されているものとする。
ここで、図11を用いて学習器300の生成方法について説明する。図11は、本開示の実施形態に係る学習器300について説明するための図である。ここでは、説明を簡略化するために、情報処理装置200Aの制御部220が学習器300を生成するものとして説明するが、学習器300は、例えば他の情報処理装置(図示省略)等によって生成されてもよい。
図11に示す例では、情報処理装置200Aは、入力画像M14を学習器300に入力した場合の出力画像M15と、学習の真値(Ground Truth)である法線画像M16との類似度を算出する。ここでは、情報処理装置200Aは、出力画像M15と法線画像M16との最小二乗誤差(L2)を算出し、L2をロスとして学習処理を行う。情報処理装置200Aは、例えば算出したロスをバックプロパゲーションすることにより、学習器300の重みを更新する。これにより、情報処理装置200Aは、撮像画像を入力すると法線画像を出力する学習器300を生成する。
なお、ここでは、機械学習の一例としてCNNを利用して学習器300を生成するとしたが、これに限定されない。機械学習としてCNN以外にも例えばRNN(Recurrent Neural Network)等、種々の手法を用いて学習器300を生成してもよい。また、上述した例では、算出したロスをバックプロパゲーションすることにより、学習器300の重みを更新するとしたが、これに限定されない。バックプロパゲーション以外にも例えば確率的勾配降下法等の任意の学習手法を用いて学習器300の重みを更新するようにしてもよい。また、上述した例では、ロスを最小二乗誤差としたがこれに限定されない。ロスを最小平均誤差としてもよい。
図6に戻る。法線算出部222は、生成した学習器300を利用して法線画像を生成する。法線算出部222は、算出した法線画像を深度算出部223に出力する。
(深度算出部)
深度算出部223は、法線画像の各画素値を入力として、上述した式(1)~(3)を用いた式変換を行ってセンサ150から対象Sまでの距離(深度)を算出する。具体的には、深度算出部223は、以下の通り、上述した式(3)を逆フーリエ変換してZ’を算出し、求める深度ZをZ=Z’×Pとして算出する。
Figure 0007609075000005
なお、上述した式(3)等の各パラメータは以下の通りである。
P(u、v):法線のx方向のフーリエ変換
Q(u、v):法線のy方向のフーリエ変換
u、v:各画素の周波数空間中での座標
P:撮像画像xyの1pixel(1画素)あたりの長さ[μm]
:求める深度のフーリエ変換
Z’:求める1pixelあたりの深度[pixel]
Z:求める深度[μm]
上記P(撮像画像xyの1pixel(1画素)あたりの長さ[μm])は、撮像装置(マイクロスコープ100)のセンサ150の構成から値が決まるパラメータである。
例えば、図12に示すように、センサ150の焦点距離がf[mm]、筐体(ヘッドマウント部10)の長さがd[mm]であるとする。また、センサ150の画素ピッチがp[μm]とすると、撮像画像の1pixelの実長Pは、P:d=p:fからP=(d×p)/f[μm]となる。なお、図12は、撮像画像の1pixelあたりの長さについて説明するための図である。
深度算出部223は、法線画像とセンサ150に関する情報とに基づき、例えば撮像画像M11の画素ごとに対象Sの深度を算出する。センサ150に関する情報とは、例えばセンサ150の構造に関する情報であり、具体的には、上述したセンサ150の焦点距離やセンサ150と対象Sとの間の距離に関する情報である。撮像装置がマイクロスコープ100の場合、センサ150と対象Sとの間の距離はヘッドマウント部10の長さdとなる。
(表示制御部)
表示制御部224は、各種の画像を液晶ディスプレイ等の表示部(図示省略)に表示させる。図13は、表示制御部224が表示部に表示させる画像M17の一例を示す図である。図13に示すように、表示制御部224は、例えば撮像画像M17a、法線画像M17b、深度画像M17c、深度グラフを示す画像M17dおよび対象Sの3D画像M17eを含む画像M17を表示部に表示させる。
撮像画像M17aは、例えば、マイクロスコープ100で撮像した撮像画像M11から平坦領域SABを切り取った画像である。法線画像M17bは、撮像画像M17aの法線情報をRGBで表示した画像である。深度画像M17cは、撮像画像M17aの各画素における深度を示す画像であり、例えば色が薄いほど深度が大きい(センサ150から対象Sまでの距離が大きい)ことを示している。深度グラフを示す画像M17dは、法線画像M17b、深度画像M17cおよび対象Sの3D画像M17eに示す直線における深度をグラフとして表示した画像である。対象Sの3D画像M17eは、深度画像M17cに基づいて対象Sを3次元表示した画像である。
このように、表示制御部224が、撮像画像や生成した法線画像および深度画像に加え、深度を示すグラフや3次元画像を表示部に表示させるようにしてもよい。
(記憶部)
記憶部230は、制御部220の処理に用いられるプログラムや演算パラメータ等を記憶するROM(Read Only Memory)、および適宜変化するパラメータ等を一時記憶するRAM(Random Access Memory)により実現される。
<2.5.深度算出処理>
図14は、本開示の第1実施形態に係る深度算出処理の一例を説明するためのフローチャートである。図14に示す深度算出処理は、情報処理装置200Aの制御部220がプログラムを実行することによって実現される。図14に示す深度算出処理は、マイクロスコープ100が撮像画像M11を撮像した後に実行される。あるいは、図14に示す深度算出処理は、ユーザからの指示に応じて実行されるようにしてもよい。
図14に示すように、制御部220は、例えば、マイクロスコープ100から撮像画像M11を取得する(ステップS101)。あるいは、制御部220が、例えばネットワーク(図示省略)を介して他の装置から撮像画像M11を取得するようにしてもよい。他の装置として、例えば他の情報処理装置やクラウドサーバ等が挙げられる。
制御部220は、センサ150に関する情報を取得する(ステップS102)。制御部220は、マイクロスコープ100からセンサ150に関する情報を取得してもよい。あるいは、記憶部230にセンサ150に関するセンサ情報を記憶している場合、制御部220は、記憶部230からセンサ情報を取得してもよい。また、制御部220は、例えばユーザからの入力によってセンサ情報を取得してもよい。
制御部220は、取得した撮像画像M11から平坦領域SABを取得する(ステップS103)。具体的に、制御部220は、撮像画像M11の各画素の輝度値が閾値th以上である領域を平坦領域SABとして抽出する。
制御部220は、取得した平坦領域SABの法線情報を算出する(ステップS104)。具体的に、制御部220は、平坦領域SABを入力画像として学習器300に入力することで、出力画像として法線情報を含む法線画像を生成する。
制御部220は、法線情報およびセンサ150に関する情報に基づき、深度情報を算出する(ステップS105)。具体的に、制御部220は、法線情報に対して、上述した式(3)等に基づく式変換を行うことで深度情報を算出する。
以上のように、第1実施形態に係る情報処理装置200Aは、制御部220を備える。制御部220は、対象Sをセンサ150が撮像した撮像画像M11を取得する。撮像画像M11は、それぞれ異なる位置に配置される複数の点光源160A、160Bから対象Sに照射された光IA、IBの反射光から得られる画像である。制御部220は、撮像画像M11の輝度値に基づき、撮像画像M11から平坦領域SABを抽出する。制御部220は、センサ150に関する情報および撮像画像M11の平坦領域SABに基づき、対象Sの表面の形状に関する形状情報(深度情報)を算出する。
特に、第1実施形態に係る情報処理装置200Aは、複数の点光源160A、160Bから同時に対象Sに照射された光IA、IBの反射光から得られる撮像画像M11から平坦領域SABを取得する。
このように、情報処理装置200Aが平坦領域SABを取得することで、深度算出部223における式変換において仮定1、2を満たすことができ、深度算出の精度低下を抑制することができる。これにより、情報処理装置200Aは、形状計測の精度を向上させることができる。
<3.第2実施形態>
上記第1実施形態では、複数の点光源160A、160Bから同時に光IA、IBが照射される間にマイクロスコープ100が撮像した撮像画像M11に基づき、情報処理装置200Aが深度を算出する場合を示した。上記例以外にも、情報処理装置200Aは、1つの点光源から光が照射される間にマイクロスコープ100が撮像した撮像画像に基づいて深度を算出してもよい。そこで、第2実施形態では、情報処理装置200Aが、複数の点光源160A、160Bの1つから対象Sに照射された光の反射光から得られる撮像画像であって、点光源160A、160Bごとに得られる撮像画像に基づき、深度を算出する例について説明する。なお、第2実施形態に係る情報処理装置200Aは、取得部221および領域取得部225の動作を除き、第1実施形態に係る情報処理装置200Aと同じ構成および動作であるため、同一符号を付し、一部の説明を省略する。
図15は、本開示の第2実施形態に係る取得部221が取得する撮像画像M21、M22について説明するための図である。
図15に示すように、マイクロスコープ100は、まず点光源160Aから光IAを照射して、対象Sを撮像する。このとき、例えば対象Sの表面の大まかな形状が平面出ない場合、対象Sの表面にはオクルージョン(光が当たらない領域)が発生する。
例えば、光源160Aから照射された光IAは、対象Sの表面のうち領域SA1には当たるが、領域SA2には当たらない。この場合、取得部221は、例えば図15に示すように、領域SA2が領域SA1に比べて暗い撮像画像M21を取得する。
続いて、マイクロスコープ100は、点光源160Bから光IBを照射して、対象Sを撮像する。このとき、点光源160Aから光IAを照射した場合と同様に、対象Sの表面にはオクルージョン(光が当たらない領域)が発生する。
例えば、光源160Bから照射された光IBは、対象Sの表面のうち領域SB1には当たるが、領域SB2には当たらない。この場合、取得部221は、例えば図15に示すように、領域SB2が領域SB1に比べて暗い撮像画像M22を取得する。
このように、取得部221は、マイクロスコープ100が2つの点光源160A、160Bを順次点灯させて撮像した撮像画像M21、M22を取得する。
領域取得部225は、取得部221が取得した撮像画像M21、M22から平坦領域SABを取得する。上述したように、平坦領域SABには、点光源160A、160Bの両方の光IA、IBが当たる。換言すると、光IA、IBの少なくとも一方が当たらないオクルージョン領域は平坦でない領域となる。このように、撮像画像M21、M22において、平坦領域SABには、点光源160A、160Bからの光IA、IBが当たるため、両方の撮像画像M21、M22において平坦領域SABの画素の輝度値はほぼ同じ値になる。
そこで、領域取得部225は、撮像画像M21、M22の輝度値に基づき、平坦領域SABを取得する。領域取得部225は、撮像画像M21、M22のそれぞれ対応する画素において輝度値の変化がない、あるいは変化が小さい(閾値th2以内である)領域を平坦領域SABとして抽出する。
領域取得部225による平坦領域SABの具体的な取得方法について、図16を用いて説明する。図16は、本開示の第2実施形態に係る領域取得部225による平坦領域の取得について説明するための図である。
領域取得部225は、例えば撮像画像M21、M22の対応する画素(x、y)の輝度値の差分D(x、y)=L1(x、y)-L2(x、y)を算出する。なお、L1(x、y)は、撮像画像M21の画素(x、y)の輝度値であり、L2(x、y)は、撮像画像M22の画素(x、y)の輝度値である。図16は、画素のyを所定の値に固定した場合の、差分Dとxとの関係を示すグラフである。なお、図15に示すように、撮像画像M21、M22の横方向をx方向、縦方向をy方向としている。
例えば図16に示す領域SB2では、撮像画像M21の輝度値が高く、撮像画像M22の輝度値が低い(図15参照)。そのため、領域SB2における差分Dの値は正方向に大きな値となる。一方、図16に示す領域SA2では、撮像画像M21の輝度値が低く、撮像画像M22の輝度値が高い(図15参照)。そのため、領域SA2における差分Dの値は負方向に大きな値となる。また、領域SB2、SA2を除く領域では、撮像画像M21、M22の両方とも同程度に輝度値が高いため、差分Dはゼロに近い値となる。
そこで、領域取得部225は、差分Dが閾値th2以上である領域SB2を、深度算出を行わない除外領域とする。また、領域取得部225は、差分Dが閾値-th2以下である領域SA2を、深度算出を行わない除外領域とする。一方、領域取得部225は、差分Dが閾値±th2の範囲内である領域を、深度算出を行う処理領域とする。領域取得部225は、処理領域と判定した領域を撮像画像M21、M22から抽出する抽出領域(平坦領域SAB)として取得する。
より具体的に、領域取得部225は、上述した差分Dの絶対値abs(D)を閾値判定し、差分Dの絶対値abs(D)が閾値th2以上の画素を含む領域を黒色の領域、閾値th2未満の画素を含む領域を白色の領域とすることで、例えば図8に示すマスク画像M12を生成する。領域取得部225は、例えば、マスク画像M12と撮像画像M21とを比較し、同座標においてマスク画像M12が白色である画素を撮像画像M21から抽出することで、平坦領域SABを取得する。あるいは、領域取得部225が、マスク画像M12と撮像画像M22とを比較して平坦領域SABを取得してもよい。なお、マスク画像M12を生成した後の処理は、第1実施形態と同じである。
第2実施形態では、点光源160A、160Bのいずれか一方を照射して撮像を行う。この場合、上述した第1実施形態のように同時に照射する場合に比べて、対象Sに当たる光量を低く抑えることができるため、撮像画像M21、M22の全体的な明るさを抑え、コントラストを高く保つことができる。そのため、撮像画像M21、M22の光が当たる領域の輝度値と光が当たらないオクルージョン領域の輝度値の差を大きくすることができる。これにより、領域取得部225は、より容易に撮像画像M21、M22から平坦領域SABを取得することができる。
一方、第1実施形態に係る情報処理装置200Aは、取得する撮像画像M11の枚数が1枚でよいため、第2実施形態に比べて画像取得時間を半分の時間に抑えることができる。また、撮像画像M21、M22を順次撮像する間に、対象Sまたはセンサ150が動いてしまうと、領域取得部225は、平坦領域SABを取得する前に、例えば撮像画像M21、M22の位置合わせ等の処理を行う必要がでてくる。そのため、リアルタイムで深度算出を行う場合は、画像取得時間が短く位置合わせ等の処理が不要である第1実施形態に係る情報処理装置200Aを用いて深度算出を行うことが望ましい。
以上のように、第2実施形態に係る情報処理装置200Aは、複数の点光源160A、160Bの1つから対象Sに照射された光IA、IBの反射光から得られる撮像画像M21、M22を、点光源160A、160Bごとに取得する。情報処理装置200Aは、点光源160A、160Bごとに取得した複数の撮像画像M21、M22間における輝度値の変化(差分Dの絶対値)が所定閾値th2未満である領域を平坦領域SABとして抽出する。
これにより、情報処理装置200Aは、撮像画像M21、M22のコントラストを高くすることができ、平坦領域SABの抽出精度を向上させることができ、形状計測の精度を向上させることができる。
なお、ここでは、マイクロスコープ100が2つの点光源160A、160Bを備える場合を例に説明したが、点光源の数は2つに限定されない。例えば、マイクロスコープ100が3つ以上の点光源を備えていてもよい。
この場合、マイクロスコープ100は、3つ以上の点光源を順次点灯させて対象Sを撮像する。情報処理装置200Aは、マイクロスコープ100が撮像した複数の撮像画像を取得する。情報処理装置200Aは、取得した複数の撮像画像間で輝度値の変化が少ない領域、換言すると複数の撮像画像の対応する画素の輝度値の差分の絶対値が所定閾値th2未満である領域を平坦領域SABとして取得する。
なお、マイクロスコープ100が3つ以上の点光源を有している場合であっても、全ての点光源を順次点灯させて対象Sを撮像する必要はない。例えば、マイクロスコープ100が3つの点光源を有している場合、3つの点光源のうち2つの点光源を順次点灯させて撮像を行ってもよい。この場合、マイクロスコープ100は、例えば複数の点光源のうち最も離れた点光源を選択して点灯させるようにしてもよい。
<4.第3実施形態>
上記第1、第2実施形態では、取得部221が取得した撮像画像から平坦領域SABを取得する場合を示した。上記例以外にも、情報処理装置200Aは、撮像画像に対してスムージング処理を行ってから平坦領域SABを取得してもよい。そこで、第3実施形態では、情報処理装置200Aが撮像画像から平坦領域SABを取得する前にスムージング処理を行う例について説明する。なお、第3実施形態に係る情報処理装置200Aは、領域取得部225の動作を除き、第1実施形態に係る情報処理装置200Aと同じ構成および動作であるため、同一符号を付し、一部の説明を省略する。
図17は、本開示の第3実施形態に係る領域取得部225によるスムージング処理を説明するための図である。
図17に示すように、取得部221が取得した撮像画像M31には、対象Sの局所の凹凸やノイズが含まれる。なお、対象Sの局所の凹凸は、深度算出部223における深度情報の算出対象である。
領域取得部225が、図17に示す撮像画像M31を画素単位で閾値判定をしようとすると、上述した凹凸やノイズの影響を受けてしまい、平坦領域SABの取得精度が低下してしまう。
そこで、領域取得部225は、閾値判定を行って平坦領域SABを取得する前に、撮像画像M31に対してスムージング処理を行い、図17の下図に示すようにスムージング画像M32を生成する。具体的に、領域取得部225は、撮像画像M31にローパスフィルターを適用することで、スムージング画像M32を生成する。領域取得部225は、スムージング画像M32に対して閾値判定を行って図8に示すマスク画像M12を生成する。
領域取得部225は、例えば、マスク画像M12と撮像画像M31とを比較し、同座標においてマスク画像M12が白色である画素を撮像画像M31から抽出することで、撮像画像M31における平坦領域SABを取得する。
このように、領域取得部225が、スムージング画像M32に基づき、マスク画像M12を生成することで、撮像画像M31の細かな凹凸やノイズの影響を低減しつつマスク画像M12を生成することができる。また、領域取得部225がマスク画像M12を用いて、撮像画像M31の平坦領域SABを取得することで、局所の凹凸を含んだ平坦領域SABを抽出することができ、対象Sの深度を精度よく算出することができる。
なお、第2実施形態で説明したように、取得部221が複数の撮像画像を取得する場合、領域取得部225は、取得した全ての撮像画像に対してスムージング処理を行う。これにより、情報処理装置200Aは、平坦領域SABの取得精度低下を抑制することができ、形状計測の精度を向上させることができる。
<5.第4実施形態>
上記第1、第2実施形態では、領域取得部225が画素の輝度値を閾値判定することで平坦領域SABを取得する場合を示した。上記例以外にも、情報処理装置200Aは、撮像画像を複数のブロックに分割することで平坦領域を取得してもよい。そこで、第4実施形態では、情報処理装置200Aが、撮像画像を複数のブロックに分割する例について説明する。なお、第4実施形態に係る情報処理装置200Aは、領域取得部225の動作を除き、第1実施形態に係る情報処理装置200Aと同じ構成および動作であるため、同一符号を付し、一部の説明を省略する。
図18は、本開示の第4実施形態に係る領域取得部225による平坦領域の取得について説明するための図である。
図18に示すように、領域取得部225は、取得部221が取得した撮像画像M41を、例えば複数のブロックに分割し、分割した1つのブロックM42を平坦領域として取得する。なお、領域取得部225は、マイクロスコープ100に関する情報に基づき、ブロックの大きさを決定し、決定した大きさのブロックで撮像画像M41を分割する。
図19に示すように、対象Sの形状の変動が撮像画像M41全体では大きかったとしても、局所領域SCでは、概ね平面とみなすことができる。なお、図19は、本開示の第4実施形態に係る領域取得部225による平坦領域の取得について説明するための図である。
例えば、対象Sが人間の皮膚である場合、経験的に5mm×5mm以下の面積ならば、上述した仮定1、仮定2を満たす平面としてよい。そこで、領域取得部225は、例えば撮像画像M41を5mm×5mm以下のブロックで分割する。撮像画像の1pixelの実長P[μm]は、上述したように、センサ150の焦点距離f[mm]、筐体(ヘッドマウント部10)の長さd[mm]、および、センサ150の画素ピッチp[μm]に基づいて算出される。
例えば上述したように、撮像画像M41を5mm×5mm以下のブロックで分割する場合、ブロックの幅w[pixel]は、w=5000/P=5000f/(d×p)となる。具体的に、センサ150の焦点距離f=16[mm]、ヘッドマウント部10の長さd=100[mm]、および、センサ150の画素ピッチp=5[μm]とすると、1pixelあたりの長さPはP=31.25[μm]となる。よって、撮像画像M41を5mm×5mm以下のブロックで分割する場合、ブロックサイズは、160×160[pixel]となる。
なお、法線算出部222および深度算出部223は、分割したブロック全てに対して法線情報および深度情報を算出する。算出対象とするブロック(平坦領域)の数が異なる点を除き、深度算出処理は図13に示す算出処理と同じである。
このように、第4実施形態では、情報処理装置200Aが、撮像画像M41を複数のブロックに分割して深度を算出する。これにより、情報処理装置200Aは、上述した式(1)~式(3)における仮定1、仮定2を満たす平坦領域を取得することができ、精度よく深度算出を行うことができる。このように、第4実施形態に係る情報処理装置200Aは、形状計測の精度を向上させることができる。
なお、ここでは、ブロックのサイズを5mm×5mmとしたが、これに限定されない。例えば、マイクロスコープ100の撮像対象は人間の皮膚以外であってもよい。撮像対象によって平坦とみなせるブロックのサイズは異なるため、ブロックのサイズは、撮像対象に応じて異なるサイズとしてもよい。
例えば、マイクロスコープ100が複数種類の対象を撮像する場合、記憶部230に撮像対象の種類と適切なブロックのサイズとを対応付けた表を記憶しておくものとする。この場合、情報処理装置200Aは、撮像対象の種類に応じてブロックのサイズを選択する。あるいは、ユーザがブロックのサイズを指定するようにしてもよい。このように、ブロックのサイズは、上述した例に限定されず、種々のサイズが選択可能である。
<6.第5実施形態>
上記第4実施形態では、領域取得部225が撮像画像をブロックに分割することで平坦領域を取得する場合を示した。上記例以外にも、情報処理装置200Aは、分割したブロックのコントラスト値に応じて平坦領域を取得してもよい。そこで、第5実施形態では、情報処理装置200Aが、分割したブロックのコントラスト値に応じて平坦領域を取得する例について説明する。なお、第5実施形態に係る情報処理装置200Aは、領域取得部225の動作を除き、第4実施形態に係る情報処理装置200Aと同じ構成および動作であるため、同一符号を付し、一部の説明を省略する。
図20は、センサ150の被写界深度を説明するための図である。センサ150(カメラ)には、レンズの絞りに応じた被写界深度(被写体がぼけずに撮影できる奥行きの範囲)が存在する。カメラの光学系には、合焦することで被写体を最も解像度高く鮮明に撮影することができる焦点距離fが存在する。例えば、図20に示す面M52は、合焦した面、すなわちピントが合っており被写体を最も解像度高く鮮明に撮影することができる面である。また、面M53から面M51の範囲が、被写体がぼけずに撮影できる奥行きの範囲(被写界深度)である。
マイクロスコープ100のようにヘッドマウント部10を対象Sに接触させて撮像する場合、対象Sまでの距離が短くなるためセンサ150に入射する光量が少なくなってしまう。そのため、センサ150に入射する光量をできるだけ増やすために、マイクロスコープ100では絞り値Fを小さく、すなわち絞りを開放して撮像を行う。そのため、図21に示すように被写界深度P1が浅くなる。なお、図21は、本開示の第5実施形態に係る領域取得部225による平坦領域の取得について説明するための図である。
図21に示すように、マイクロスコープ100の被写界深度P1が浅いため、被写界深度P1に含まれる対象S、すなわち大まかな形状においてセンサ150からの距離がほぼ同じである対象Sの表面SD1はぼけずに撮影することができる。一方、被写界深度P1に含まれない対象S、すなわちセンサ150からの距離が表面SD1とは異なる表面SD2、SD3はぼやけて撮像されてしまう。
そこで、領域取得部225は、被写界深度P1に含まれる撮像領域を平坦領域として取得する。領域取得部225は、ピントがあってぼけずに撮像された領域を撮像画像から取得することで、平坦領域を取得する。かかる点について図22を用いて説明する。図22は、本開示の第5実施形態に係る領域取得部225による平坦領域の取得について説明するための図である。
領域取得部225は、取得部221が取得した撮像画像M54を複数のブロックに分割する。なお、ここで分割するブロックのサイズは、上述した第4実施形態のブロックのサイズと同じであっても異なっていてもよい。
続いて、領域取得部225は、分割したブロックごとにコントラスト値を算出する。コントラスト値contrastは、contrast=(Lmax-Lmin)/(Lmax+Lmin)で求まる値である。なお、Lmaxは、ブロック内の輝度値Lの最大値であり、Lmax=max(L(x、y))である。また、Lminは、ブロック内の輝度値Lの最小値であり、Lmin=min(L(x、y))である。
ピントがあってぼけずに撮像された領域はコントラスト値contrastが高く、ピントがあっておらずぼやけて撮像された領域はコントラスト値contrastが低い。そこで、領域取得部225は、算出したコントラスト値contrastに基づいて平坦領域を取得する。具体的に、領域取得部225は、算出したコントラスト値contrastと閾値th3とをブロックごとに比較する。領域取得部225は、図22に示すように、コントラスト値contrastが閾値th3以上であるブロックを、深度算出を行う平坦領域として抽出する。
法線算出部222および深度算出部223は、領域取得部225が取得した平坦領域の法線情報および深度情報を算出する。算出方法は第1実施形態と同じである。なお、法線算出部222および深度算出部223は、ブロック単位で法線情報および深度情報を算出してもよい。あるいは、法線算出部222および深度算出部223は、領域取得部225が平坦領域として抽出した領域、すなわちコントラスト値contrastが閾値th3以上である複数のブロックをまとめて、法線情報および深度情報を算出してもよい。
以上のように、第5実施形態では、情報処理装置200Aは、撮像画像M54を複数のブロックに分割する。情報処理装置200Aは、分割したブロックごとに算出したコントラスト値contrastに応じて平坦領域を取得する。
このように、情報処理装置200Aは、撮像画像M54のコントラスト値contrastに応じて平坦領域を取得することで、上述した式(1)~式(3)における仮定1、仮定2を満たす平坦領域における深度を算出することができる。これにより、第5実施形態に係る情報処理装置200Aは、形状計測の精度を向上させることができる。
<7.第6実施形態>
上記第5実施形態では、領域取得部225が1枚の撮像画像を複数のブロックに分割し、分割したブロックごとのコントラスト値に応じて平坦領域を取得する場合を示した。上記例以外にも、情報処理装置200Aは、被写体深度が異なる複数の撮像画像をそれぞれ複数のブロックに分割し、分割したブロックごとのコントラスト値に応じて平坦領域を取得してもよい。そこで、第6実施形態では、情報処理装置200Aが複数の撮像画像それぞれのブロックごとのコントラスト値に応じて平坦領域を取得する例について説明する。なお、第6実施形態に係る情報処理装置200Aは、取得部221および領域取得部225の動作を除き、第5実施形態に係る情報処理装置200Aと同じ構成および動作であるため、同一符号を付し、一部の説明を省略する。
図23は、本開示の第6実施形態に係る取得部221が取得する複数の撮像画像について説明するための図である。取得部221は、被写体深度が異なる複数の撮像画像を取得する。図23に示すように、マイクロスコープ100のセンサ150を動かすことで、マイクロスコープ100は、被写体深度(フォーカス面)がそれぞれ異なる撮像画像を複数取得する。図23の例では、マイクロスコープ100が、センサ150を上下に動かしながら対象Sを3回撮像することで3枚の撮像画像を撮像する。マイクロスコープ100は、例えば、1回目に被写体深度がセンサ150に最も近い撮像画像を撮像し、3回目に被写体深度がセンサ150から最も遠い撮像画像を撮像する。マイクロスコープ100は、例えば2回目に被写体深度が1回目と2回目の中間に位置する撮像画像を取得する。
なお、図23に示す被写体深度と撮像枚数は一例であり、これに限定されない。例えば、マイクロスコープ100が2枚、あるいは、4枚以上の撮像画像を撮像するようにしてもよい。また、各撮像時の被写体深度の範囲が連続、あるいは、一部重複していてもよい。
取得部221は、例えばマイクロスコープ100を制御して、被写体深度がそれぞれ異なる撮像画像M61~M63を取得する。
図24は、本開示の第6実施形態に係る領域取得部225による平坦領域の取得について説明するための図である。ここでは、取得部221が、マイクロスコープ100を制御して、図23に示すように被写体深度が異なる撮像画像を撮像させることで、図24に示すような撮像画像M61~M63を取得するものとする。
領域取得部225は、撮像画像M61~M63それぞれをブロックに分割する。領域取得部225は、分割したブロックごとにコントラスト値contrastを算出する。領域取得部225は、算出したブロックのコントラスト値contrastと閾値th3とを比較し、コントラスト値contrastが閾値th3以上であるブロックを平坦領域として抽出する。
例えば、図24では、領域取得部225は、撮像画像M61のコントラスト値contrastが閾値th3以上であるブロックを含む領域M61Aを撮像画像M61から抽出する。同様に、領域取得部225は、撮像画像M62のコントラスト値contrastが閾値th3以上であるブロックを含む領域M62Aを撮像画像M62から抽出する。また、撮像画像M63のコントラスト値contrastが閾値th3以上であるブロックを含む領域M63Aを撮像画像M63から抽出する。
最後に、領域取得部225は、抽出した領域M61A~M63Aを合成して平坦領域を含む画像M64を生成する。換言すると、領域取得部225は、撮像画像M61~M63のピントが合っており、ぼけていない領域を合成して画像M64を生成する。
なお、画像M64を平坦領域として、法線算出部222および深度算出部223が実行する処理は第5実施形態と同じである。
以上のように、第6実施形態に係る情報処理装置200Aは、被写体深度が異なる複数の撮像画像を取得する。情報処理装置200Aは、取得した複数の撮像画像をそれぞれ複数のブロックに分割する。情報処理装置200Aは、分割したブロックごとのコントラスト値contrastに応じて平坦領域を取得する。
これにより、情報処理装置200Aは、コントラスト値contrastに応じて平坦領域を取得することで、上述した式(1)~式(3)における仮定1、仮定2を満たす平坦領域における深度を算出することができる。また、情報処理装置200Aは、複数の撮像画像から平坦領域を取得することで、コントラスト値contrastに応じた平坦領域を拡大することができる。
なお、ここでは、領域取得部225が領域M61A~M63Aを合成するとしたが、これに限定されない。例えば、法線算出部222および深度算出部223がそれぞれ領域M61A~M63Aごとに法線情報および深度情報を算出し、例えば表示制御部224が表示部に結果を表示させるときに、領域M61A~M63Aを合成した結果を表示するようにしてもよい。
<8.第7実施形態>
上記第1~6実施形態では、領域取得部225が撮像画像から平坦領域を取得する場合を示した。上記例以外にも、法線算出部222が算出した法線情報から高周波成分を抽出するようにしてもよい。そこで、第7実施形態では、法線情報から高周波成分を抽出して深度情報を算出する例について説明する。
図25は、本開示の第7実施形態に係る情報処理装置200Bの構成例を示す図である。図25に示す情報処理装置200Bは、制御部220Bが領域取得部225を備えておらず、法線周波数分離部226Bを備える点を除き、図6に示す情報処理装置200Aと同じ構成要素を有する。
法線算出部222は、取得部221が取得した撮像画像の各画素における法線情報を算出する。
法線周波数分離部226Bは、法線情報から高周波成分(以下、法線高周波ともいう)を分離する。ここで、上述したように、深度算出部223の算出対象である深度情報は、対象Sの局所的な凹凸形状に関する情報である。換言すると、深度算出部223が算出する深度情報は、対象Sの表面の形状の高周波成分であり、低周波成分を含んでいないものである。法線周波数分離部226は、法線情報から大まかな形状の法線情報である低周波成分を除去して法線高周波を生成する。法線高周波は、低周波成分(大まかな形状の法線情報)を含まないため、上述した式(1)~式(3)の仮定1、仮定2を満たすものとなる。
図26は、本開示の第7実施形態に係る法線周波数分離部226Bによる周波数分離を説明するための図である。
例えば、図26の上図に示すような法線情報Xが法線周波数分離部226Bに入力されたとする。なお、図26では、実際の対象Sの表面を点線で、対象Sの表面の大まかな形状を実線で、対象Sの表面の法線を矢印で示している。
この場合、法線周波数分離部226Bは、例えば畳み込みフィルタF(X)を用いて法線情報の高周波成分と低周波成分を分離し、図26の中図に示す法線高周波Yを抽出する。なお、畳み込みフィルタF(X)は、法線情報を入力としたときに高周波成分を抽出するように予め設計されたフィルタである。
深度算出部223は、法線周波数分離部226Bが分離した法線高周波を上述した式(1)~式(3)等を用いて式変換することで、仮定1、仮定2を満たす状態で深度情報Z(図26の下図参照)を算出することができる。
このように、情報処理装置200Bは、法線情報から高周波成分を抽出し、抽出した高周波成分に基づいて深度情報を算出する。
これにより、情報処理装置200Bは、仮定1、仮定2を満たす状態で上述した(1)~式(3)等を用いた式変換を行えるため、深度情報の算出精度の低下を抑制することができる。そのため、情報処理装置200Bは、周波数分離を行わない場合に比べて、形状計測の精度を向上させることができる。
<9.第8実施形態>
上記第7実施形態では、法線情報の高周波成分から深度情報を算出する場合を示した。上記例以外にも、法線情報の高周波成分に加え、低周波成分も用いて深度情報を算出するようにしてもよい。そこで、第8実施形態では、法線情報から高周波成分および低周波成分を分離して深度情報を算出する例について説明する。
図27は、本開示の第8実施形態に係る情報処理装置200Cの構成例を示す図である。図27に示す情報処理装置200Cは、制御部220Cの法線周波数分離部226Cおよび深度算出部223Cを除き、図25に示す情報処理装置200Bと同じ構成要素を有する。
法線周波数分離部226Cは、法線情報から高周波成分と低周波成分(以下、法線低周波ともいう)を分離して、それぞれ深度算出部223Cに出力する。上述したように、対象Sの局所的な凹凸形状が法線情報の高周波成分に寄与する。また、対象Sの大まかな形状が法線情報の低周波成分に寄与する。かかる点について図28を用いて説明する。
図28は、本開示の第8実施形態に係る法線周波数分離部226Cによる周波数分離を説明するための図である。
例えば、図28(a)に示すような法線情報Xが法線周波数分離部226Cに入力されたとする。なお、図28では、実際の対象Sの表面を点線で、対象Sの表面の大まかな形状を実線で、対象Sの表面の法線を矢印で示している。
この場合、法線周波数分離部226Cは、例えば畳み込みフィルタF(X)を用いて法線情報の高周波成分を抽出し、法線高周波YH=F(X)を抽出する。図28(b)の上図に示すように、法線高周波YHで表される表面形状は、大まかな形状(法線低周波YL)が除去され、平坦な面に局所的な凹凸がある形状となる。なお、畳み込みフィルタF(X)は、法線情報を入力としたときに高周波成分を抽出するように予め設計されたフィルタである。
また、法線周波数分離部226Cは、法線情報Xから法線高周波YHを減算することで法線低周波YLを抽出する。すなわち、法線周波数分離部226Cは、法線低周波YL=X-YH=X-F(X)を算出することで、法線低周波YLを抽出する。図28(b)の下図に示すように、法線低周波YLで表される表面形状は、局所的な凹凸形状が除去された大まかな形状となる。
深度算出部223Cは、法線高周波YHおよび法線低周波YLそれぞれに対して式変換を行い、図28(c)に示すように、深度情報の高周波成分ZH(以下、深度高周波ともいう)および低周波成分ZL(以下、深度低周波ともいう)を算出する。
ここで、法線高周波YHは、上述したように、式変換の仮定1、仮定2を満たす状態であるため、深度算出部223Cは、精度よく深度高周波ZHを算出できる。
また、深度算出部223Cは、上述した式(1)~式(3)のλおよびμをそれぞれゼロ(λ=0、μ=0)として、法線低周波YLの式変換を行うことで、深度低周波ZLを算出する。このように、cost項の係数をゼロにすることで、仮定1、仮定2を排除した状態で深度低周波ZLを算出することができる。なお、法線低周波YLは、法線高周波YHが除去された法線情報Xの低周波成分であり、対象Sの表面形状のうち局所的な凹凸形状を除去した状態を表す法線情報である。そのため、深度算出部223Cは、仮定1、仮定2を除去した状態で法線低周波YLの式変換を行っても、精度よく深度低周波ZLを算出することができる。
深度算出部223Cは、算出した深度高周波ZHおよび深度低周波ZLを合成することで、図28(d)に示す深度情報Zを算出する。
このように、情報処理装置200Cは、撮像画像から算出した法線情報を周波数分離する。情報処理装置200Cは、分離した周波数ごとに深度情報を算出し、算出した深度情報を合成して対象Sの深度情報を算出する。
これにより、情報処理装置200Cは、局所的な凹凸形状の深度に加え、大まかな形状の深度を算出することができる。また、情報処理装置200Cは、これらの深度を合成することで、対象Sの深度情報をより精度よく算出することができる。
なお、ここでは、深度算出部223Cが、上述した式(1)~式(3)のλおよびμをそれぞれゼロ(λ=0、μ=0)として、法線低周波YLの式変換を行うとしたが、これに限定されない。深度算出部223Cが仮定1および仮定2の影響が小さい状態で深度低周波ZLを算出できればよく、λおよびμの値を、仮定1および仮定2が深度低周波ZLの算出に影響を与えない程度に小さい値としてもよい。
<10.第9実施形態>
上記第1実施形態では、真値である画像M16を用いて学習器300の学習を行う場合を示した。上記例以外にも、情報処理装置200Aは、画像M16をシフトさせた画像を正解データとして学習器300の学習を行ってもよい。そこで、第9実施形態では、情報処理装置200Aが画像M16をシフトさせた画像を正解データとして学習器300の学習を行う例について図29を用いて説明する。図29は、第9実施形態に係る学習器300の学習方法について説明するための図である。なお、第9実施形態に係る情報処理装置200Aは、法線算出部222の動作を除き、第1実施形態に係る情報処理装置200Aと同じ構成および動作であるため、同一符号を付し、一部の説明を省略する。
上述したように、法線算出部222は、撮像画像M14を入力データとし、法線画像M16を正解データ(真値)として学習器300の重みを更新して学習を行う。ここで、撮像画像M14は、例えばマイクロスコープ100等のRGB画像を撮像するデバイスによって撮像される。一方、法線画像M16は、例えば非接触3D計測器のように直接対象Sの表面形状を計測するデバイスが計測した形状情報に基づいて生成される。
このように、入力画像を取得するデバイスの光学系と正解画像を取得するデバイスの光学系とが異なっていると、入力画像と正解画像の位置合わせが難しくなってしまう。そのため、入力画像を学習器300に入力することで得られる出力画像および正解画像が数画素ずれてしまうことがある。このように、出力画像と正解画像とがずれた状態で学習器300の学習を行うと、入力画像に対してぼけた出力画像を出力する学習器300が生成されてしまう。
そこで、図29に示すように、法線算出部222は、学習器(CNN)300のフィルタ係数(重み)を決める学習時に、法線画像M16をx方向およびy方向にシフトさせ、シフト後の法線画像M16と出力画像M15との類似度を算出する。なお、法線算出部222は、法線画像M16をx方向において0から±xの範囲でシフトさせ、y方向において0から±yの範囲でシフトさせるものとする。
法線算出部222は、例えばシフト後の法線画像M16と出力画像M15との最小二乗誤差をシフト量ごとに算出する。法線算出部222は、算出した最小二乗誤差のうち最小値を最終的なロスとして、学習器300の重みを更新する。なお、最小二乗誤差のかわりに、シフト後の法線画像M16と出力画像M15との最小平均誤差をシフト量ごとに算出し、算出した最小平均誤差のうち最小値を最終的なロスとしてもよい。
このように、学習器300は、真値画像(法線画像)をシフトしたシフト画像(シフト後の法線画像)と出力画像M15との類似度に応じて重みが更新されて生成される。より具体的には、学習器300は、真値画像(法線画像)を異なるシフト量(画素数)でシフトした複数のシフト画像(シフト後の法線画像、正解候補画像ともいう)それぞれと、出力画像M15との類似度に応じて重みが更新されて生成される。
これにより、学習器300の出力画像M15と法線画像M16との位置合わせを行わなくても、両方の画像の位置があった状態で学習器300の重みを更新することができる。これにより、学習器300の学習精度をより向上させることができる。また、法線算出部222がかかる学習器300を用いて法線情報を算出することで、より精度よく法線情報を算出することができ、情報処理装置200Aの形状計測の精度をより向上させることができる。
<11.第10実施形態>
上記第1実施形態では、撮像画像から平坦領域を取得する場合を示した。上記例以外にも、情報処理装置200Aは、撮像画像を複数の領域に分割し、分割した領域ごとに平坦領域を取得するようにしてもよい。そこで、第10実施形態では、情報処理装置200Aが撮像画像を複数の領域に分割し、分割した領域ごとに平坦領域を取得する例について図30を用いて説明する。図30は、本開示の第10実施形態に係る情報処理装置200Aの制御部220Dの構成例を示す図である。
図30に示すように、情報処理装置200Aの制御部220Dは、領域取得部225の代わりに領域取得部225Dを備える点、および、深度合成部227をさらに備える点を除き、図6に示す制御部220と同じである。
図30に示すように、領域取得部225Dは、取得部221が取得した撮像画像M71を3つの分割領域M71A~M71Cに分割する。なお、領域取得部225Dが分割する領域の数は、3つに限定されず、2つであっても4つ以上であってもよい。
領域取得部225Dは、分割領域M71A~M71Cごとに分割平坦領域M72A~M72Cを取得する。法線算出部222は、分割平坦領域M72A~M72Cそれぞれに基づき、法線領域画像M73A~M73Cを算出し、深度算出部223は、法線領域画像M73A~M73Cに基づき、深度領域画像M74A~M74Cを算出する。
深度合成部227は、深度領域画像M74A~M74Cを合成して深度画像M75を生成し、表示制御部224に出力する。
なお、深度合成部227は、深度領域画像M74A~M74のほかにも、分割平坦領域M72A~M72Cの合成や、法線領域画像M73A~M73Cの合成を行うようにしてもよい。また、各部における処理は、領域ごとに順次行われてもよく、領域ごとに並列に行われていてもよい。
このように、情報処理装置200Aは、撮像画像を複数の分割領域M71A~M71Cに分割し、分割領域M71A~M71Cごとに分割平坦領域M72A~M72Cを取得する。また、情報処理装置200Aは、分割平坦領域M72A~M72Cごとに法線情報および深度情報を算出する。
これにより、各部の処理対象である撮像画像(領域)のサイズを小さくすることができ、処理負担を軽減したり、処理速度を向上させたりすることができる。
<12.その他の実施形態>
上述した各実施形態に係る処理は、上記各実施形態以外にも種々の異なる形態にて実施されてよい。
上記各実施形態では、ヘッドマウント部10がマイクロスコープ100の先端にマウントされる筒部である例を示した。しかし、ヘッドマウント部10は、対象Sとマイクロスコープ100のセンサ150との距離を一定にするための構造物であれば、必ずしも筒状の形態でなくてもよい。
<13.適応例>
図31を用いて、第1~第10実施形態に係る情報処理装置200A~200Cの適用例について説明する。図31は、本開示に係る技術(本技術)が適用され得る、カプセル型内視鏡を用いた患者の体内情報取得システムの概略的な構成の一例を示すブロック図である。
体内情報取得システム10001は、カプセル型内視鏡10100と、外部制御装置10200とから構成される。
カプセル型内視鏡10100は、検査時に、患者によって飲み込まれる。カプセル型内視鏡10100は、撮像機能および無線通信機能を有し、患者から自然排出されるまでの間、胃や腸等の臓器の内部を蠕動運動等によって移動しつつ、当該臓器の内部の画像(以下、体内画像ともいう)を所定の間隔で順次撮像し、その体内画像についての情報を体外の外部制御装置10200に順次無線送信する。
外部制御装置10200は、体内情報取得システム10001の動作を統括的に制御する。また、外部制御装置10200は、カプセル型内視鏡10100から送信されてくる体内画像についての情報を受信し、受信した体内画像についての情報に基づいて、表示装置(図示しない)に当該体内画像を表示するための画像データを生成する。
体内情報取得システム10001では、このようにして、カプセル型内視鏡10100が飲み込まれてから排出されるまでの間、患者の体内の様子を撮像した体内画像を随時得ることができる。
カプセル型内視鏡10100と外部制御装置10200の構成および機能についてより詳細に説明する。
カプセル型内視鏡10100は、カプセル型の筐体10101を有し、その筐体10101内には、光源部10111、撮像部10112、画像処理部10113、無線通信部10114、給電部10115、電源部10116、および制御部10117が収納されている。
光源部10111は、例えばLED(Light Emitting Diode)等の光源から構成され、撮像部10112の撮像視野に対して光を照射する。
撮像部10112は、撮像素子、および当該撮像素子の前段に設けられる複数のレンズからなる光学系から構成される。観察対象である体組織に照射された光の反射光(以下、観察光という)は、当該光学系によって集光され、当該撮像素子に入射する。撮像部10112では、撮像素子において、そこに入射した観察光が光電変換され、その観察光に対応する画像信号が生成される。撮像部10112によって生成された画像信号は、画像処理部10113に提供される。
画像処理部10113は、CPUやGPU(Graphics Processing Unit)等のプロセッサによって構成され、撮像部10112によって生成された画像信号に対して各種の信号処理を行う。画像処理部10113は、信号処理を施した画像信号を、RAWデータとして無線通信部10114に提供する。
無線通信部10114は、画像処理部10113によって信号処理が施された画像信号に対して変調処理等の所定の処理を行い、その画像信号を、アンテナ10114Aを介して外部制御装置10200に送信する。また、無線通信部10114は、外部制御装置10200から、カプセル型内視鏡10100の駆動制御に関する制御信号を、アンテナ10114Aを介して受信する。無線通信部10114は、外部制御装置10200から受信した制御信号を制御部10117に提供する。
給電部10115は、受電用のアンテナコイル、当該アンテナコイルに発生した電流から電力を再生する電力再生回路、および昇圧回路等から構成される。給電部10115では、いわゆる非接触充電の原理を用いて電力が生成される。
電源部10116は、二次電池によって構成され、給電部10115によって生成された電力を蓄電する。図31では、図面が煩雑になることを避けるために、電源部10116からの電力の供給先を示す矢印等の図示を省略しているが、電源部10116に蓄電された電力は、光源部10111、撮像部10112、画像処理部10113、無線通信部10114、および制御部10117に供給され、これらの駆動に用いられ得る。
制御部10117は、CPU等のプロセッサによって構成され、光源部10111、撮像部10112、画像処理部10113、無線通信部10114、および、給電部10115の駆動を、外部制御装置10200から送信される制御信号に従って適宜制御する。
外部制御装置10200は、CPU、GPU等のプロセッサ、又はプロセッサとメモリ等の記憶素子が混載されたマイクロコンピュータ若しくは制御基板等で構成される。外部制御装置10200は、カプセル型内視鏡10100の制御部10117に対して制御信号を、アンテナ10200Aを介して送信することにより、カプセル型内視鏡10100の動作を制御する。カプセル型内視鏡10100では、例えば、外部制御装置10200からの制御信号により、光源部10111における観察対象に対する光の照射条件が変更され得る。また、外部制御装置10200からの制御信号により、撮像条件(例えば、撮像部10112におけるフレームレート、露出値等)が変更され得る。また、外部制御装置10200からの制御信号により、画像処理部10113における処理の内容や、無線通信部10114が画像信号を送信する条件(例えば、送信間隔、送信画像数等)が変更されてもよい。
また、外部制御装置10200は、カプセル型内視鏡10100から送信される画像信号に対して、各種の画像処理を施し、撮像された体内画像を表示装置に表示するための画像データを生成する。当該画像処理としては、例えば現像処理(デモザイク処理)、高画質化処理(帯域強調処理、超解像処理、ノイズリダクション処理、手ブレ補正処理等)、拡大処理(電子ズーム処理)等、それぞれ単独で、あるいは、組み合わせて、各種の信号処理を行うことができる。外部制御装置10200は、表示装置の駆動を制御して、生成した画像データに基づいて撮像された体内画像を表示させる。あるいは、外部制御装置10200は、生成した画像データを記録装置(図示しない)に記録させたり、印刷装置(図示しない)に印刷出力させたりしてもよい。
以上、本開示に係る技術が適用され得る体内情報取得システムの一例について説明した。本開示に係る技術は、以上説明した構成のうち、例えば、外部制御装置10200に適用され得る。外部制御装置10200に本開示に係る技術を適用することにより、カプセル型内視鏡10100が撮像した体内画像から体内の表面形状を計測することができる。
(内視鏡手術システムへの適用例)
本開示に係る技術は、さらに、内視鏡手術システムに適用されてもよい。図32は、本開示に係る技術(本技術)が適用され得る内視鏡手術システムの概略的な構成の一例を示す図である。
図32では、術者(医師)11131が、内視鏡手術システム11000を用いて、患者ベッド11133上の患者11132に手術を行っている様子が図示されている。図示するように、内視鏡手術システム11000は、内視鏡11100と、気腹チューブ11111やエネルギー処置具11112等の、その他の術具11110と、内視鏡11100を支持する支持アーム装置11120と、内視鏡下手術のための各種の装置が搭載されたカート11200と、から構成される。
内視鏡11100は、先端から所定の長さの領域が患者11132の体腔内に挿入される鏡筒11101と、鏡筒11101の基端に接続されるカメラヘッド11102と、から構成される。図示する例では、硬性の鏡筒11101を有するいわゆる硬性鏡として構成される内視鏡11100を図示しているが、内視鏡11100は、軟性の鏡筒を有するいわゆる軟性鏡として構成されてもよい。
鏡筒11101の先端には、対物レンズが嵌め込まれた開口部が設けられている。内視鏡11100には光源装置11203が接続されており、当該光源装置11203によって生成された光が、鏡筒11101の内部に延設されるライトガイドによって当該鏡筒の先端まで導光され、対物レンズを介して患者11132の体腔内の観察対象に向かって照射される。なお、内視鏡11100は、直視鏡であってもよいし、斜視鏡又は側視鏡であってもよい。
カメラヘッド11102の内部には光学系および撮像素子が設けられており、観察対象からの反射光(観察光)は当該光学系によって当該撮像素子に集光される。当該撮像素子によって観察光が光電変換され、観察光に対応する電気信号、すなわち観察像に対応する画像信号が生成される。当該画像信号は、RAWデータとしてCCU(カメラコントロールユニット)11201に送信される。
CCU11201は、CPUやGPU等によって構成され、内視鏡11100および表示装置11202の動作を統括的に制御する。さらに、CCU11201は、カメラヘッド11102から画像信号を受け取り、その画像信号に対して、例えば現像処理(デモザイク処理)等の、当該画像信号に基づく画像を表示するための各種の画像処理を施す。
表示装置11202は、CCU11201からの制御により、当該CCU11201によって画像処理が施された画像信号に基づく画像を表示する。
光源装置11203は、例えばLED(Light Emitting Diode)等の光源から構成され、術部等を撮影する際の照射光を内視鏡11100に供給する。
入力装置11204は、内視鏡手術システム11000に対する入力インタフェースである。ユーザは、入力装置11204を介して、内視鏡手術システム11000に対して各種の情報の入力や指示入力を行うことができる。例えば、ユーザは、内視鏡11100による撮像条件(照射光の種類、倍率および焦点距離等)を変更する旨の指示等を入力する。
処置具制御装置11205は、組織の焼灼、切開又は血管の封止等のためのエネルギー処置具11112の駆動を制御する。気腹装置11206は、内視鏡11100による視野の確保および術者の作業空間の確保の目的で、患者11132の体腔を膨らめるために、気腹チューブ11111を介して当該体腔内にガスを送り込む。レコーダ11207は、手術に関する各種の情報を記録可能な装置である。プリンタ11208は、手術に関する各種の情報を、テキスト、画像又はグラフ等各種の形式で印刷可能な装置である。
なお、内視鏡11100に術部を撮影する際の照射光を供給する光源装置11203は、例えばLED、レーザ光源又はこれらの組み合わせによって構成される白色光源から構成することができる。RGBレーザ光源の組み合わせにより白色光源が構成される場合には、各色(各波長)の出力強度および出力タイミングを高精度に制御することができるため、光源装置11203において撮像画像のホワイトバランスの調整を行うことができる。また、この場合には、RGBレーザ光源それぞれからのレーザ光を時分割で観察対象に照射し、その照射タイミングに同期してカメラヘッド11102の撮像素子の駆動を制御することにより、RGBそれぞれに対応した画像を時分割で撮像することも可能である。当該方法によれば、当該撮像素子にカラーフィルタを設けなくても、カラー画像を得ることができる。
また、光源装置11203は、出力する光の強度を所定の時間ごとに変更するようにその駆動が制御されてもよい。その光の強度の変更のタイミングに同期してカメラヘッド11102の撮像素子の駆動を制御して時分割で画像を取得し、その画像を合成することにより、いわゆる黒つぶれおよび白とびのない高ダイナミックレンジの画像を生成することができる。
また、光源装置11203は、特殊光観察に対応した所定の波長帯域の光を供給可能に構成されてもよい。特殊光観察では、例えば、体組織における光の吸収の波長依存性を利用して、通常の観察時における照射光(すなわち、白色光)に比べて狭帯域の光を照射することにより、粘膜表層の血管等の所定の組織を高コントラストで撮影する、いわゆる狭帯域光観察(Narrow Band Imaging)が行われる。あるいは、特殊光観察では、励起光を照射することにより発生する蛍光により画像を得る蛍光観察が行われてもよい。蛍光観察では、体組織に励起光を照射し当該体組織からの蛍光を観察すること(自家蛍光観察)、又はインドシアニングリーン(ICG)等の試薬を体組織に局注するとともに当該体組織にその試薬の蛍光波長に対応した励起光を照射し蛍光像を得ること等を行うことができる。光源装置11203は、このような特殊光観察に対応した狭帯域光および/又は励起光を供給可能に構成され得る。
図33は、図32に示すカメラヘッド11102およびCCU11201の機能構成の一例を示すブロック図である。
カメラヘッド11102は、レンズユニット11401と、撮像部11402と、駆動部11403と、通信部11404と、カメラヘッド制御部11405と、を有する。CCU11201は、通信部11411と、画像処理部11412と、制御部11413と、を有する。カメラヘッド11102とCCU11201とは、伝送ケーブル11400によって互いに通信可能に接続されている。
レンズユニット11401は、鏡筒11101との接続部に設けられる光学系である。鏡筒11101の先端から取り込まれた観察光は、カメラヘッド11102まで導光され、当該レンズユニット11401に入射する。レンズユニット11401は、ズームレンズおよびフォーカスレンズを含む複数のレンズが組み合わされて構成される。
撮像部11402は、撮像素子で構成される。撮像部11402を構成する撮像素子は、1つ(いわゆる単板式)であってもよいし、複数(いわゆる多板式)であってもよい。撮像部11402が多板式で構成される場合には、例えば各撮像素子によってRGBそれぞれに対応する画像信号が生成され、それらが合成されることによりカラー画像が得られてもよい。あるいは、撮像部11402は、3D(Dimensional)表示に対応する右目用および左目用の画像信号をそれぞれ取得するための1対の撮像素子を有するように構成されてもよい。3D表示が行われることにより、術者11131は術部における生体組織の奥行きをより正確に把握することが可能になる。なお、撮像部11402が多板式で構成される場合には、各撮像素子に対応して、レンズユニット11401も複数系統設けられ得る。
また、撮像部11402は、必ずしもカメラヘッド11102に設けられなくてもよい。例えば、撮像部11402は、鏡筒11101の内部に、対物レンズの直後に設けられてもよい。
駆動部11403は、アクチュエータによって構成され、カメラヘッド制御部11405からの制御により、レンズユニット11401のズームレンズおよびフォーカスレンズを光軸に沿って所定の距離だけ移動させる。これにより、撮像部11402による撮像画像の倍率および焦点が適宜調整され得る。
通信部11404は、CCU11201との間で各種の情報を送受信するための通信装置によって構成される。通信部11404は、撮像部11402から得た画像信号をRAWデータとして伝送ケーブル11400を介してCCU11201に送信する。
また、通信部11404は、CCU11201から、カメラヘッド11102の駆動を制御するための制御信号を受信し、カメラヘッド制御部11405に供給する。当該制御信号には、例えば、撮像画像のフレームレートを指定する旨の情報、撮像時の露出値を指定する旨の情報、並びに/又は撮像画像の倍率および焦点を指定する旨の情報等、撮像条件に関する情報が含まれる。
なお、上記のフレームレートや露出値、倍率、焦点等の撮像条件は、ユーザによって適宜指定されてもよいし、取得された画像信号に基づいてCCU11201の制御部11413によって自動的に設定されてもよい。後者の場合には、いわゆるAE(Auto Exposure)機能、AF(Auto Focus)機能およびAWB(Auto White Balance)機能が内視鏡11100に搭載されていることになる。
カメラヘッド制御部11405は、通信部11404を介して受信したCCU11201からの制御信号に基づいて、カメラヘッド11102の駆動を制御する。
通信部11411は、カメラヘッド11102との間で各種の情報を送受信するための通信装置によって構成される。通信部11411は、カメラヘッド11102から、伝送ケーブル11400を介して送信される画像信号を受信する。
また、通信部11411は、カメラヘッド11102に対して、カメラヘッド11102の駆動を制御するための制御信号を送信する。画像信号や制御信号は、電気通信や光通信等によって送信することができる。
画像処理部11412は、カメラヘッド11102から送信されたRAWデータである画像信号に対して各種の画像処理を施す。
制御部11413は、内視鏡11100による術部等の撮像、および、術部等の撮像により得られる撮像画像の表示に関する各種の制御を行う。例えば、制御部11413は、カメラヘッド11102の駆動を制御するための制御信号を生成する。
また、制御部11413は、画像処理部11412によって画像処理が施された画像信号に基づいて、術部等が映った撮像画像を表示装置11202に表示させる。この際、制御部11413は、各種の画像認識技術を用いて撮像画像内における各種の物体を認識してもよい。例えば、制御部11413は、撮像画像に含まれる物体のエッジの形状や色等を検出することにより、鉗子等の術具、特定の生体部位、出血、エネルギー処置具11112の使用時のミスト等を認識することができる。制御部11413は、表示装置11202に撮像画像を表示させる際に、その認識結果を用いて、各種の手術支援情報を当該術部の画像に重畳表示させてもよい。手術支援情報が重畳表示され、術者11131に提示されることにより、術者11131の負担を軽減することや、術者11131が確実に手術を進めることが可能になる。
カメラヘッド11102およびCCU11201を接続する伝送ケーブル11400は、電気信号の通信に対応した電気信号ケーブル、光通信に対応した光ファイバ、又はこれらの複合ケーブルである。
ここで、図33の例では、伝送ケーブル11400を用いて有線で通信が行われていたが、カメラヘッド11102とCCU11201との間の通信は無線で行われてもよい。
以上、本開示に係る技術が適用され得る内視鏡手術システムの一例について説明した。本開示に係る技術は、以上説明した構成のうち、例えば、CCU11201に適用され得る。具体的には、上述した制御部220、220B、220Cは、画像処理部11412に適用することができる。画像処理部11412に本開示に係る技術を適用することにより、内視鏡11100が撮像した体内画像から体内の表面形状を計測することができる。
なお、ここでは、一例として内視鏡手術システムについて説明したが、本開示に係る技術は、その他、例えば、顕微鏡手術システム等に適用されてもよい。
<14.補足>
以上、添付図面を参照しながら本開示の好適な実施形態について詳細に説明したが、本開示の技術的範囲はかかる例に限定されない。本開示の技術分野における通常の知識を有する者であれば、請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本開示の技術的範囲に属するものと了解される。
上記各実施形態において説明した各処理のうち、自動的に行われるものとして説明した処理の全部または一部を手動的に行うこともでき、あるいは、手動的に行われるものとして説明した処理の全部または一部を公知の方法で自動的に行うこともできる。この他、上記文書中や図面中で示した処理手順、具体的名称、各種のデータやパラメータを含む情報については、特記する場合を除いて任意に変更することができる。例えば、各図に示した各種情報は、図示した情報に限られない。
また、図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部または一部を、各種の負荷や使用状況などに応じて、任意の単位で機能的または物理的に分散・統合して構成することができる。
また、上述してきた各実施形態および変形例は、処理内容を矛盾させない範囲で適宜組み合わせることが可能である。また、上記実施形態では、画像処理装置としてマイクロスコープを例に挙げたが、本開示の画像処理は、マイクロスコープ以外の撮影機器にも適用可能である。
また、本明細書に記載された効果は、あくまで説明的または例示的なものであって限定的ではない。つまり、本開示にかかる技術は、上記の効果とともに、または上記の効果に代えて、本明細書の記載から当業者には明らかな他の効果を奏しうる。
なお、以下のような構成も本開示の技術的範囲に属する。
(1)
対象をセンサが撮像した撮像画像を取得し、
前記撮像画像は、それぞれ異なる位置に配置される複数の光源から前記対象に照射された光の反射光から得られる画像であり、
前記撮像画像の輝度値に基づき、前記撮像画像から平坦領域を抽出し、
前記センサに関する情報および前記撮像画像の前記平坦領域に基づき、前記対象の表面の形状に関する形状情報を算出する制御部、
を備える情報処理装置。
(2)
前記制御部は、
複数の前記光源から同時に前記対象に照射された前記光の前記反射光から得られる前記撮像画像を取得する、
(1)に記載の情報処理装置。
(3)
前記制御部は、
前記撮像画像の輝度値が、所定閾値以上である領域を前記平坦領域とする、
(1)または(2)に記載の情報処理装置。
(4)
前記制御部は、
複数の前記光源の1つから前記対象に照射された前記光の前記反射光から得られる前記撮像画像を、前記光源ごとに取得し、
前記光源ごとに取得した複数の前記撮像画像間における前記輝度値の変化が所定閾値未満である領域を前記平坦領域として抽出する、
(1)に記載の情報処理装置。
(5)
前記制御部は、
前記撮像画像にスムージング処理を行ってスムージング画像を生成し、
前記スムージング画像に基づき、前記撮像画像から前記平坦領域を抽出する、
(1)~(4)のいずれか1つに記載の情報処理装置。
(6)
前記制御部は、
前記撮像画像を複数の分割領域に分割し、
複数の前記分割領域ごとに前記平坦領域を抽出し、形状情報を算出する、
(1)~(5)のいずれか1つに記載の情報処理装置。
(7)
前記制御部は、
前記撮像画像における前記平坦領域の法線情報を算出し、
前記法線情報を機械学習により生成されたモデルに入力して前記形状情報を得る、
(1)~(6)のいずれか1つに記載の情報処理装置。
(8)
前記モデルは、
複数の正解候補画像と前記モデルの出力データとの比較結果に基づいて重みを更新することで生成され、
複数の前記正解候補画像は、正解画像をそれぞれ異なる画素数ずつずらして生成される、
(7)に記載の情報処理装置。
(9)
前記モデルは、
複数の正解候補画像と前記出力データとの最小二乗誤差をそれぞれ算出し、
複数の前記最小二乗誤差の最小値に基づいて前記重みを更新することで生成される、
(8)に記載の情報処理装置。
(10)
前記制御部は、
前記撮像画像の輝度値から算出されるコントラスト値に基づき、前記平坦領域を抽出する、
(1)に記載の情報処理装置。
(11)
前記制御部は、
前記撮像画像を複数の分割領域に分割し、
複数の前記分割領域ごとに前記コントラスト値を算出し、
前記コントラスト値に応じて前記分割領域が平坦であるか否かを判定することで、前記平坦領域を抽出する、
(1)に記載の情報処理装置。
(12)
前記制御部は、
前記センサの被写界深度面がそれぞれ異なる複数の前記撮像画像を取得し、
複数の前記撮像画像からそれぞれ前記平坦領域を抽出する、
(10)または(11)に記載の情報処理装置。
(13)
対象をセンサが撮像した撮像画像を取得し、
前記撮像画像は、それぞれ異なる位置に配置される複数の光源から前記対象に照射された光の反射光から得られる画像であり、
前記撮像画像の輝度値に基づき、前記撮像画像から平坦領域を抽出し、
前記センサに関する情報および前記撮像画像の前記平坦領域に基づき、前記対象の表面の形状に関する形状情報を算出する、
情報処理方法。
10 ヘッドマウント部
100 マイクロスコープ
150 センサ
160 点光源
200 情報処理装置
220 制御部
221 取得部
222 法線算出部
223 深度算出部
224 表示制御部
225 領域取得部
226 法線周波数分離部
227 深度合成部
230 記憶部

Claims (11)

  1. 対象をセンサが撮像した撮像画像であって、それぞれ異なる位置に配置される複数の光源から同時に前記対象に照射された光の反射光から得られる前記撮像画像を取得し
    前記撮像画像の輝度値に基づき、前記撮像画像から前記撮像画像の輝度値が、所定閾値以上である平坦領域を抽出し、
    前記センサに関する情報および前記撮像画像の前記平坦領域に基づき、前記対象の表面の形状に関する形状情報を算出する制御部、
    を備える情報処理装置。
  2. 前記制御部は、
    複数の前記光源の1つから前記対象に照射された前記光の前記反射光から得られる前記撮像画像を、前記光源ごとに取得し、
    前記光源ごとに取得した複数の前記撮像画像間における前記輝度値の変化が所定閾値未満である領域を前記平坦領域として抽出する、
    請求項1に記載の情報処理装置。
  3. 前記制御部は、
    前記撮像画像にスムージング処理を行ってスムージング画像を生成し、
    前記スムージング画像に基づき、前記撮像画像から前記平坦領域を抽出する、
    請求項に記載の情報処理装置。
  4. 前記制御部は、
    前記撮像画像を複数の分割領域に分割し、
    複数の前記分割領域ごとに前記平坦領域を抽出し、形状情報を算出する、
    請求項に記載の情報処理装置。
  5. 前記制御部は、
    前記撮像画像における前記平坦領域の法線情報を算出し、
    前記法線情報を機械学習により生成されたモデルに入力して前記形状情報を得る、
    請求項に記載の情報処理装置。
  6. 前記モデルは、
    複数の正解候補画像と前記モデルの出力データとの比較結果に基づいて重みを更新することで生成され、
    複数の前記正解候補画像は、正解画像をそれぞれ異なる画素数ずつずらして生成される、
    請求項に記載の情報処理装置。
  7. 前記モデルは、
    複数の正解候補画像と前記出力データとの最小二乗誤差をそれぞれ算出し、
    複数の前記最小二乗誤差の最小値に基づいて前記重みを更新することで生成される、
    請求項に記載の情報処理装置。
  8. 前記制御部は、
    前記撮像画像の輝度値から算出されるコントラスト値に基づき、前記平坦領域を抽出する、
    請求項1に記載の情報処理装置。
  9. 前記制御部は、
    前記撮像画像を複数の分割領域に分割し、
    複数の前記分割領域ごとに前記コントラスト値を算出し、
    前記コントラスト値に応じて前記分割領域が平坦であるか否かを判定することで、前記平坦領域を抽出する、
    請求項に記載の情報処理装置。
  10. 前記制御部は、
    前記センサの被写界深度面がそれぞれ異なる複数の前記撮像画像を取得し、
    複数の前記撮像画像からそれぞれ前記平坦領域を抽出する、
    請求項に記載の情報処理装置。
  11. 対象をセンサが撮像した撮像画像であって、それぞれ異なる位置に配置される複数の光源から同時に前記対象に照射された光の反射光から得られる前記撮像画像を取得し
    前記撮像画像の輝度値に基づき、前記撮像画像から前記撮像画像の輝度値が、所定閾値以上である平坦領域を抽出し、
    前記センサに関する情報および前記撮像画像の前記平坦領域に基づき、前記対象の表面の形状に関する形状情報を算出する、
    情報処理方法。
JP2021556079A 2019-11-15 2020-11-09 情報処理装置および情報処理方法 Active JP7609075B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2019207103 2019-11-15
JP2019207103 2019-11-15
PCT/JP2020/041676 WO2021095672A1 (ja) 2019-11-15 2020-11-09 情報処理装置および情報処理方法

Publications (2)

Publication Number Publication Date
JPWO2021095672A1 JPWO2021095672A1 (ja) 2021-05-20
JP7609075B2 true JP7609075B2 (ja) 2025-01-07

Family

ID=75912071

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021556079A Active JP7609075B2 (ja) 2019-11-15 2020-11-09 情報処理装置および情報処理方法

Country Status (4)

Country Link
US (1) US12361573B2 (ja)
JP (1) JP7609075B2 (ja)
CN (1) CN114667431B (ja)
WO (1) WO2021095672A1 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7799375B2 (ja) * 2020-07-02 2026-01-15 キヤノン株式会社 画像処理装置、画像処理方法及びプログラム
CN114359063B (zh) * 2021-12-03 2025-08-22 上海微创医疗机器人(集团)股份有限公司 目标亮度确定方法、光源控制方法和系统、设备和介质
JP7693585B2 (ja) * 2022-03-03 2025-06-17 株式会社日立製作所 撮影装置および認証装置
WO2025110014A1 (ja) * 2023-11-22 2025-05-30 株式会社Ihi 生成装置および生成方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007139067A1 (ja) 2006-05-29 2007-12-06 Panasonic Corporation 画像高解像度化装置、画像高解像度化方法、画像高解像度化プログラムおよび画像高解像度化システム
WO2018159738A1 (ja) 2017-03-03 2018-09-07 ソニーセミコンダクタソリューションズ株式会社 固体撮像装置および電子機器
JP2018200328A (ja) 2018-09-18 2018-12-20 株式会社キーエンス 検査装置、検査方法およびプログラム
WO2019139694A1 (en) 2018-01-11 2019-07-18 Exxonmobil Chemical Patents Inc. Polyethylene compositions and articles manufactured therefrom

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10307011A (ja) * 1995-12-05 1998-11-17 New Kurieishiyon:Kk 表面検査方法および表面検査装置
JP3980722B2 (ja) 1997-04-03 2007-09-26 株式会社モリテックス Ccdマイクロスコープ
JP3738291B2 (ja) * 2003-06-09 2006-01-25 住友大阪セメント株式会社 三次元形状測定装置
JP2008253498A (ja) 2007-04-04 2008-10-23 Shiseido Co Ltd マイクロスコープの先端ドーム
JP5152077B2 (ja) * 2009-04-01 2013-02-27 ソニー株式会社 生体像提示装置、生体像提示方法及びプログラム並びに生体像提示システム
JP5908668B2 (ja) * 2010-04-12 2016-04-26 富士フイルム株式会社 可搬型放射線撮影装置
KR101087180B1 (ko) * 2011-06-22 2011-11-28 동국대학교 경주캠퍼스 산학협력단 신뢰성 있는 금속 표면 3차원 형상 추출 기법 및 시스템
KR101087172B1 (ko) * 2011-07-13 2011-11-28 동국대학교 경주캠퍼스 산학협력단 고온 금속 표면 3차원 형상 추출 시스템 및 방법
JP6063315B2 (ja) * 2013-03-26 2017-01-18 富士フイルム株式会社 真贋判定システム,特徴点登録装置およびその動作制御方法,ならびに照合判定装置およびその動作制御方法
US9903710B2 (en) 2015-06-05 2018-02-27 Nippon Steel & Sumitomo Metal Corporation Shape inspection apparatus for metallic body and shape inspection method for metallic body
JP2017146938A (ja) * 2016-02-19 2017-08-24 富士通株式会社 書籍検出装置、書籍検出方法及び書籍検出用コンピュータプログラム
WO2017194276A1 (en) * 2016-05-13 2017-11-16 Basf Se System and method for detecting plant diseases
JPWO2018079255A1 (ja) * 2016-10-24 2019-09-12 パナソニックIpマネジメント株式会社 画像処理装置、画像処理方法、および画像処理プログラム
JP6820515B2 (ja) * 2017-03-06 2021-01-27 株式会社東京精密 表面形状測定装置及び表面形状測定方法
JP7286268B2 (ja) * 2018-02-15 2023-06-05 キヤノン株式会社 画像処理方法、画像処理装置、撮像装置、画像処理プログラム、および、記憶媒体
JP7163049B2 (ja) * 2018-04-03 2022-10-31 キヤノン株式会社 情報処理装置、情報処理方法及びプログラム
CN109272459B (zh) * 2018-08-20 2020-12-01 Oppo广东移动通信有限公司 图像处理方法、装置、存储介质及电子设备
KR102644946B1 (ko) * 2019-05-02 2024-03-08 삼성전자주식회사 영상 처리 방법 및 장치
US11651505B2 (en) * 2020-01-27 2023-05-16 Cognex Corporation Systems and method for vision inspection with multiple types of light

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007139067A1 (ja) 2006-05-29 2007-12-06 Panasonic Corporation 画像高解像度化装置、画像高解像度化方法、画像高解像度化プログラムおよび画像高解像度化システム
WO2018159738A1 (ja) 2017-03-03 2018-09-07 ソニーセミコンダクタソリューションズ株式会社 固体撮像装置および電子機器
WO2019139694A1 (en) 2018-01-11 2019-07-18 Exxonmobil Chemical Patents Inc. Polyethylene compositions and articles manufactured therefrom
JP2018200328A (ja) 2018-09-18 2018-12-20 株式会社キーエンス 検査装置、検査方法およびプログラム

Also Published As

Publication number Publication date
WO2021095672A1 (ja) 2021-05-20
CN114667431B (zh) 2025-10-28
US12361573B2 (en) 2025-07-15
US20230044097A1 (en) 2023-02-09
CN114667431A (zh) 2022-06-24
JPWO2021095672A1 (ja) 2021-05-20

Similar Documents

Publication Publication Date Title
JP7609075B2 (ja) 情報処理装置および情報処理方法
US11004197B2 (en) Medical image processing apparatus, medical image processing method, and program
WO2018084003A1 (ja) 医療用画像処理装置、医療用画像処理方法、プログラム
JP6983067B2 (ja) 医療用観察装置及び医療用観察方法
JP2013034753A (ja) 内視鏡装置
CA2627611A1 (en) Imaging system and method to improve depth perception
US11699215B2 (en) Imaging device, method and program for producing images of a scene having an extended depth of field with good contrast
WO2018159083A1 (ja) 内視鏡システム、プロセッサ装置、及び、内視鏡システムの作動方法
CN212326346U (zh) 一种内窥镜成像系统
JPWO2018221041A1 (ja) 医療用観察システム及び医療用観察装置
JP2023164610A (ja) 画像処理装置、画像処理方法および画像処理システム
WO2017145531A1 (ja) 医療用画像処理装置、システム、方法及びプログラム
CN115191915A (zh) 一种内窥镜成像系统
KR20200026135A (ko) 달팽이관에서의 미세 혈류 측정 방법 및 장치
US12364386B2 (en) Medical image generation apparatus, medical image generation method, and medical image generation program
WO2022019057A1 (ja) 医療用アーム制御システム、医療用アーム制御方法及びプログラム
JP7456385B2 (ja) 画像処理装置、および画像処理方法、並びにプログラム
JP5422180B2 (ja) 内視鏡装置およびその作動方法
CN118512135B (zh) 一种双光源荧光鼻颅底镜影像处理系统及装置
WO2020116067A1 (ja) 医療システム、情報処理装置及び情報処理方法
JPWO2019053804A1 (ja) 内視鏡装置、内視鏡装置の作動方法及びプログラム
US20240065525A1 (en) Method, computer program, and data processing unit for creating at least one correction value for correcting fluorescence intensities in a fluorescence image, and optical observation system
US20240366074A1 (en) Imaging systems with multiple fold optical path
JP5508140B2 (ja) 眼底撮像装置及びその処理方法
JP7140113B2 (ja) 内視鏡

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230925

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20240806

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240911

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: 20241119

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20241202

R150 Certificate of patent or registration of utility model

Ref document number: 7609075

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150