JP6927046B2 - 情報処理装置、スペックルイメージングシステム、及び情報処理方法 - Google Patents
情報処理装置、スペックルイメージングシステム、及び情報処理方法 Download PDFInfo
- Publication number
- JP6927046B2 JP6927046B2 JP2017553697A JP2017553697A JP6927046B2 JP 6927046 B2 JP6927046 B2 JP 6927046B2 JP 2017553697 A JP2017553697 A JP 2017553697A JP 2017553697 A JP2017553697 A JP 2017553697A JP 6927046 B2 JP6927046 B2 JP 6927046B2
- Authority
- JP
- Japan
- Prior art keywords
- speckle
- contrast
- information processing
- brightness
- image
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims description 77
- 230000010365 information processing Effects 0.000 title claims description 65
- 238000003672 processing method Methods 0.000 title claims description 9
- 239000012530 fluid Substances 0.000 claims description 68
- 230000017531 blood circulation Effects 0.000 claims description 36
- 230000001427 coherent effect Effects 0.000 claims description 25
- 230000010354 integration Effects 0.000 claims description 24
- 239000006185 dispersion Substances 0.000 claims description 12
- 210000004204 blood vessel Anatomy 0.000 claims description 6
- 238000000034 method Methods 0.000 description 18
- 238000004891 communication Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 9
- 239000008280 blood Substances 0.000 description 5
- 210000004369 blood Anatomy 0.000 description 5
- 230000007423 decrease Effects 0.000 description 4
- 208000031481 Pathologic Constriction Diseases 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000005286 illumination Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 3
- 208000037804 stenosis Diseases 0.000 description 3
- 230000036262 stenosis Effects 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- CPBQJMYROZQQJC-UHFFFAOYSA-N helium neon Chemical compound [He].[Ne] CPBQJMYROZQQJC-UHFFFAOYSA-N 0.000 description 2
- 230000001771 impaired effect Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000001356 surgical procedure Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000005653 Brownian motion process Effects 0.000 description 1
- XKRFYHLGVUSROY-UHFFFAOYSA-N argon Substances [Ar] XKRFYHLGVUSROY-UHFFFAOYSA-N 0.000 description 1
- 229910052786 argon Inorganic materials 0.000 description 1
- 238000005537 brownian motion Methods 0.000 description 1
- 238000011088 calibration curve Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 229910052743 krypton Inorganic materials 0.000 description 1
- DNNSSWSSYDEUBZ-UHFFFAOYSA-N krypton atom Chemical compound [Kr] DNNSSWSSYDEUBZ-UHFFFAOYSA-N 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 229910044991 metal oxide Inorganic materials 0.000 description 1
- 150000004706 metal oxides Chemical class 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
- G01P5/26—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting optical wave
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
- A61B5/0261—Measuring blood flow using optical means, e.g. infrared light
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/1455—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/74—Details of notification to user or communication with user or patient ; user input means
- A61B5/742—Details of notification to user or communication with user or patient ; user input means using visual displays
- A61B5/7445—Display arrangements, e.g. multiple display units
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/704—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow using marked regions or existing inhomogeneities within the fluid stream, e.g. statistically occurring variations in a fluid parameter
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/704—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow using marked regions or existing inhomogeneities within the fluid stream, e.g. statistically occurring variations in a fluid parameter
- G01F1/708—Measuring the time taken to traverse a fixed distance
- G01F1/7086—Measuring the time taken to traverse a fixed distance using optical detecting arrangements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/704—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow using marked regions or existing inhomogeneities within the fluid stream, e.g. statistically occurring variations in a fluid parameter
- G01F1/708—Measuring the time taken to traverse a fixed distance
- G01F1/712—Measuring the time taken to traverse a fixed distance using auto-correlation or cross-correlation detection means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
- G01P5/18—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the time taken to traverse a fixed distance
- G01P5/22—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the time taken to traverse a fixed distance using auto-correlation or cross-correlation detection means
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
- G06T2207/30104—Vascular flow; Blood flow; Perfusion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Physics & Mathematics (AREA)
- Physiology (AREA)
- Hematology (AREA)
- Cardiology (AREA)
- Fluid Mechanics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Multimedia (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Optics & Photonics (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Description
本開示は、情報処理装置、スペックルイメージングシステム、及び情報処理方法に関する。
従来、スペックル血流イメージによって流速を表示する方法としてMulti-Exposure Speckle Imaging法などが提案されている（非特許文献１参照）。
"Ashwin B. Patharathyet al.,Robust flow measurement with multi-exposure speckle imaging",Optics Express (2008)
非特許文献１に開示された技術では、流速を測定する前提としてスペックルのコントラストを取得する際に、様々な露光時間による撮影が必要であり、それに応じたレーザー強度のコヒーレントな光源による照明が必要である。露光時間を変化させる場合、照射するレーザーの強度を調節する必要があり、照射するレーザー光強度を露光時間に対応させて調節することは容易ではない。
イメージャーの露光時間のダイナミックレンジは、レーザーの強度のダイナミックレンジより広く設定可能である。仮に、強力なレーザーを用いて、ＮＤフィルターなどを用いて照射強度を制御することが行われているが、レーザー光の一部しか使わない構成になるため効率的とはいえない。
本開示はかかる問題点に鑑みてなされたもので、その目的は、流体速度を測定する前提としてスペックルパターンのコントラストを簡便かつ効率的に取得することのできる情報処理装置、スペックルイメージングシステム、及び情報処理方法を提供することにある。
本願発明者は、前記の目的を解決するために鋭意研究を行った結果、撮像素子を用いて複数回撮像することにより得られた複数枚のスペックル画像の輝度を積算することに着目し、本開示を完成させるに至った。
即ち、本開示では、まず、コヒーレント光が照射された撮像対象から得られる散乱光を、撮像素子を用いて複数回撮像することにより得られた複数枚のスペックル画像の輝度を積算する輝度積算部と、
前記輝度積算部で積算したスペックル積算画像に基づいて、スペックルパターンのコントラストを算出するコントラスト算出部とを備える、情報処理装置を提供する。
前記輝度積算部で積算したスペックル積算画像に基づいて、スペックルパターンのコントラストを算出するコントラスト算出部とを備える、情報処理装置を提供する。
本開示では、次に、請求項１乃至１１のいずれか一項に記載の情報処理装置と、
撮像対象にコヒーレント光を照射する光源と、
コヒーレント光が照射された撮像対象から得られる散乱光を、撮像素子を用いて複数回撮像して前記複数枚のスペックル画像を出力する撮像装置と、
画像を表示する表示装置と、を備えるスペックルイメージングシステムを提供する。
撮像対象にコヒーレント光を照射する光源と、
コヒーレント光が照射された撮像対象から得られる散乱光を、撮像素子を用いて複数回撮像して前記複数枚のスペックル画像を出力する撮像装置と、
画像を表示する表示装置と、を備えるスペックルイメージングシステムを提供する。
本開示では、次に、コヒーレント光が照射された撮像対象から得られる散乱光を、撮像素子を用いて複数回撮像することにより得られた複数枚のスペックル画像の輝度を積算する輝度積算ステップと、
前記輝度積算ステップで積算したスペックル積算画像に基づいて、スペックルパターンのコントラストを算出するコントラスト算出ステップと、を少なくとも行う、情報処理方法を提供する。
前記輝度積算ステップで積算したスペックル積算画像に基づいて、スペックルパターンのコントラストを算出するコントラスト算出ステップと、を少なくとも行う、情報処理方法を提供する。
本開示によれば、流体速度を測定する前提としてスペックルパターンのコントラストを簡便かつ効率的に取得することができる。
なお、ここに記載された効果は、必ずしも限定されるものではなく、本開示中に記載されたいずれかの効果であってもよい。
なお、ここに記載された効果は、必ずしも限定されるものではなく、本開示中に記載されたいずれかの効果であってもよい。
以下、本開示を実施するための好適な形態について、図面を参照しながら説明する。以下に説明する実施形態は、本開示の代表的な実施形態の一例を示したものであり、これにより本開示の範囲が狭く解釈されることはない。なお、説明は以下の順序で行う。
１．第１実施形態（情報処理装置１）
（１）輝度積算部１１
（２）コントラスト算出部１２
（３）流体速度算出部１３
（４）表示制御部
（５）演算のフロー例
２．第２実施形態（スペックルイメージングシステム１０）
（１）光源１４
（２）撮像装置１５
（３）表示装置１６
（４）記憶装置１７
（５）撮像対象Ｏ
（６）スペックルイメージングの第１フロー例
（７）スペックルイメージングの第２フロー例
３．第３実施形態（情報処理方法）
４．ハードウェア構成
１．第１実施形態（情報処理装置１）
（１）輝度積算部１１
（２）コントラスト算出部１２
（３）流体速度算出部１３
（４）表示制御部
（５）演算のフロー例
２．第２実施形態（スペックルイメージングシステム１０）
（１）光源１４
（２）撮像装置１５
（３）表示装置１６
（４）記憶装置１７
（５）撮像対象Ｏ
（６）スペックルイメージングの第１フロー例
（７）スペックルイメージングの第２フロー例
３．第３実施形態（情報処理方法）
４．ハードウェア構成
＜１．第１実施形態（情報処理装置１）＞
図１は、本開示の第１実施形態に係る情報処理装置１を模式的に示す模式概念図である。本開示に係る情報処理装置１は、大別して、輝度積算部１１、コントラスト算出部１２、及び流体速度算出部１３を備える。必要に応じて、表示制御部などを更に備えることも可能である。以下、各部について詳細に説明する。
図１は、本開示の第１実施形態に係る情報処理装置１を模式的に示す模式概念図である。本開示に係る情報処理装置１は、大別して、輝度積算部１１、コントラスト算出部１２、及び流体速度算出部１３を備える。必要に応じて、表示制御部などを更に備えることも可能である。以下、各部について詳細に説明する。
（１）輝度積算部１１
輝度積算部１１は、コヒーレント光が照射された撮像対象から得られる散乱光を、撮像素子を用いて複数回撮像することにより得られた複数枚のスペックル画像の輝度を積算する。
輝度積算部１１は、コヒーレント光が照射された撮像対象から得られる散乱光を、撮像素子を用いて複数回撮像することにより得られた複数枚のスペックル画像の輝度を積算する。
（２）コントラスト算出部１２
コントラスト算出部１２は、前記輝度積算部１１で積算したスペックル積算画像に基づいて、スペックルパターンのコントラストを算出する。本開示では、非侵襲、非接触で血流を測定する前提としてスペックルパターンのコントラストが算出される。血液のような散乱体が移動する部分では、スペックルが変化してしまい、輝度の濃淡は低下する。この濃淡の低下具合の指標の一つとしてスペックルパターンのコントラストがある。スペックルパターンのコントラストＫは、輝度値をＩとすると、以下の式で定義される。
コントラスト算出部１２は、前記輝度積算部１１で積算したスペックル積算画像に基づいて、スペックルパターンのコントラストを算出する。本開示では、非侵襲、非接触で血流を測定する前提としてスペックルパターンのコントラストが算出される。血液のような散乱体が移動する部分では、スペックルが変化してしまい、輝度の濃淡は低下する。この濃淡の低下具合の指標の一つとしてスペックルパターンのコントラストがある。スペックルパターンのコントラストＫは、輝度値をＩとすると、以下の式で定義される。
ここでは、濃淡の低下具合の指標として、スペックルコントラストを挙げたが、この指標はスペックルコントラストに限定するものではない。
スペックルコントラストＫは、相関時間τcと露光時間Tに依存する。従来、流速を測定する方法として、Multi Exposure Speckle Contrastなどが提案されている（図２参照）。図２は、露光時間とスペックルコントラストと流速の関係を示すグラフである。図２に示されるように、従来の方法（Multi-exposure speckle imaging法）では、露光時間を変えて様々な値に振って測定する必要があった。
相関時間τcは、流速や粘性と相関がある物理量であり、流速が速いと小さい値になり、流速が遅いと大きい値になる。上述のように、従来のMulti-exposure speckle imaging法において、スペックルコントラストから流速を求めるには、露光時間を変えて様々な値に振って測定する必要があった。これに対し、本開示では、従来のMulti-exposure speckle imaging法より少ない露光時間のパターンによって流速を算出できるようにした。
本開示に係るコントラスト算出部１２は、スペックル積算画像の局所的な画像領域における輝度の分散値及び平均値を求め、求めた輝度の分散値及び平均値に基づいて、スペックルコントラストを算出する。コントラスト算出部１２は、前記求めた輝度の分散値の平方根をとることで輝度の標準偏差を算出し、上記の式（１）に前記求めた輝度の標準偏差及び平均値を代入することによって、スペックルコントラストを算出する。
（３）流体速度算出部１３
流体速度算出部１３は、前記複数枚のスペックル画像の露光時間を積算した積算露光時間と、前記コントラスト算出部１２で算出したスペックルパターンのコントラストとに基づいて、前記撮像対象の流体速度を算出する。前記流体速度は、例えば、血管内における血流速度である。前記複数枚のスペックル画像は、１０ｍｓ以下の露光時間で撮像された画像である。
流体速度算出部１３は、前記複数枚のスペックル画像の露光時間を積算した積算露光時間と、前記コントラスト算出部１２で算出したスペックルパターンのコントラストとに基づいて、前記撮像対象の流体速度を算出する。前記流体速度は、例えば、血管内における血流速度である。前記複数枚のスペックル画像は、１０ｍｓ以下の露光時間で撮像された画像である。
露光時間の上限値を１０ｍｓに設定するメリットについては、血流の場合、露光時間１５ｍｓ程度で撮影すると、１ｍｍ／ｓｅｃ程度でスペックルコントラスト（スペックルの分散)が小さくなる。そのため、スペックルコントラストを積算しても値は変わらない。測定したい血流程度の流速では、露光時間１０ｍｓの場合、既にスペックルの分散が下がっているため、積算のメリットが出にくくなるからである。
図２のグラフより１０−３ｓｅｃ,１０−４ｓｅｃの場合、それぞれの流速のスペックルバリアンスの値が十分下がりきっていないので、スペックルコントラストを積算することの効果が現れる。一方、１０−２ｓｅｃよりも大きい値になると、スペックルバリアンスの値が下がりきっていて、積算する効果が得られない。そのため、露光時間を１０−３ｓｅｃ以下,や、１０−４ｓｅｃ程度に設定されることが好ましい。
上述のように、流体速度算出部１３は、積算露光時間を用いて血流速度を算出する。一般に、露光時間の長さに依存して血流の流速感度は異なる。露光中に血液が移動していないように見えるほど短い露光時間で撮像したスペックル画像を平均した画像から、その領域の血流速度を求めることができる。たとえば、２枚、…、１０枚、…、１００枚、…、１０００枚、…の平均画像からは、その露光時間の２倍、…、１０倍、…、１００倍、…、１０００倍、…で撮像したのとほぼ等価の画像が得られる。
流体速度算出部１３は、得られた画像の局所的なスペックルパターンのコントラストＫと、積算露光時間Ｔとを下記の式（２）に代入することによって、相関時間τcを求める。
流体速度算出部１３は、前記積算露光時間と前記スペックルコントラストに基づいて相関時間を求め、求めた相関時間と、予め定められた相関時間とを比較することによって、血流速度を算出する。なお、予め定められた相関時間とは、検量線から得られた相関時間や、予め取得しておいた相関時間である。なお、相関時間を予め測定していない場合は、プロットで得られた相関時間を直接比較することで相対的な流速を求めることができる。
上記式（２）以外にも、露光時間、スペックルコントラスト、及び、相関時間の関係を示す下記の式（３）を用いることができる。式（３）のβは、既知の流速の血流をプロットした際に決定した値を用いることができる。
（４）表示制御部
表示制御部は、表示部に画像を表示させる。表示制御部は、前記流体速度算出部１３で算出された流体速度をマッピングし、流体速度の分布を前記表示部に更に表示させることも可能である。
表示制御部は、表示部に画像を表示させる。表示制御部は、前記流体速度算出部１３で算出された流体速度をマッピングし、流体速度の分布を前記表示部に更に表示させることも可能である。
（５）演算のフロー例
図３は、本開示に係る情報処理装置１を用いて行う演算の一例を示すフローチャートである。以下、時系列に沿って、フロー例を説明する。
図３は、本開示に係る情報処理装置１を用いて行う演算の一例を示すフローチャートである。以下、時系列に沿って、フロー例を説明する。
（ａ）輝度の積算
まず、ステップＳＴ１０１において、輝度積算部１１は、複数枚のスペックル画像の輝度を積算する。
まず、ステップＳＴ１０１において、輝度積算部１１は、複数枚のスペックル画像の輝度を積算する。
（ｂ）スペックルパターンのコントラストの算出
次に、ステップＳＴ１０２において、コントラスト算出部１２は、輝度積算部１１で積算したスペックル積算画像に基づいて、スペックルパターンのコントラストＫを算出する。
次に、ステップＳＴ１０２において、コントラスト算出部１２は、輝度積算部１１で積算したスペックル積算画像に基づいて、スペックルパターンのコントラストＫを算出する。
（ｃ）相関時間の算出
続いて、流体速度算出部１３は、積算露光時間Ｔと、スペックルパターンのコントラストＫとに基づいて、撮像対象の流体速度を算出する。ステップＳＴ１０３において、流体速度算出部１３は、積算露光時間ＴとコントラストＫとに基づいて、相関時間τcを算出する。具体的には、積算露光時間ＴとコントラストＫとを上記の式（２）に代入し、相関時間τcを算出する。
続いて、流体速度算出部１３は、積算露光時間Ｔと、スペックルパターンのコントラストＫとに基づいて、撮像対象の流体速度を算出する。ステップＳＴ１０３において、流体速度算出部１３は、積算露光時間ＴとコントラストＫとに基づいて、相関時間τcを算出する。具体的には、積算露光時間ＴとコントラストＫとを上記の式（２）に代入し、相関時間τcを算出する。
（ｄ）流体速度の算出
そして、ステップＳＴ１０４において、流体速度算出部１３は、算出した相関時間τcと、予め定められた相関時間τcとを比較することによって、流体速度を算出する。
そして、ステップＳＴ１０４において、流体速度算出部１３は、算出した相関時間τcと、予め定められた相関時間τcとを比較することによって、流体速度を算出する。
＜２．第２実施形態（スペックルイメージングシステム１０）＞
図４（Ａ）は、本開示の第２実施形態に係るスペックルイメージングシステム１０の配線例を示すブロック図である。図４（Ｂ）は、本開示の第２実施形態に係るスペックルイメージングシステム１０の内部構成を示すブロック図である。なお、図４（Ａ）、（Ｂ）においては、前述した本開示に係る情報処理装置１の構成要素と同じものには、同じ符号を付し、その詳細な説明は省略する。本開示に係るスペックルイメージングシステム１０は、大別して、情報処理装置１、光源１４、撮像装置１５及び表示装置１６を備える。また、必要に応じて、記憶装置１７などを更に備えることも可能である。
図４（Ａ）は、本開示の第２実施形態に係るスペックルイメージングシステム１０の配線例を示すブロック図である。図４（Ｂ）は、本開示の第２実施形態に係るスペックルイメージングシステム１０の内部構成を示すブロック図である。なお、図４（Ａ）、（Ｂ）においては、前述した本開示に係る情報処理装置１の構成要素と同じものには、同じ符号を付し、その詳細な説明は省略する。本開示に係るスペックルイメージングシステム１０は、大別して、情報処理装置１、光源１４、撮像装置１５及び表示装置１６を備える。また、必要に応じて、記憶装置１７などを更に備えることも可能である。
（１）光源１４
光源１４は、撮像対象Ｏにコヒーレント光を照射する。コヒーレント光とは、光束内の任意の二点における光波の位相関係が時間的に不変で一定であり、任意の方法で光束を分割した後、大きな光路差を与えて再び重ね合わせても完全な干渉性を示す光をいう。
光源１４は、撮像対象Ｏにコヒーレント光を照射する。コヒーレント光とは、光束内の任意の二点における光波の位相関係が時間的に不変で一定であり、任意の方法で光束を分割した後、大きな光路差を与えて再び重ね合わせても完全な干渉性を示す光をいう。
光源１４が照射するコヒーレント光の光源の種類は、本技術の効果を損なわない限り特に限定されない。一例としては、レーザー光等を挙げることができる。レーザー光を発する光源１４としては、例えば、アルゴンイオン（Ar）レーザー、ヘリウム−ネオン（He-Ne）レーザー、ダイ（dye）レーザー、クリプトン（Cr）レーザー、半導体レーザー、または、半導体レーザーと波長変換光学素子を組み合わせた固体レーザー等を、１種又は２種以上、自由に組み合わせて用いることができる。
（２）撮像装置１５
撮像装置１５は、コヒーレント光が照射された撮像対象Ｏから得られる散乱光を、撮像素子を用いて複数回撮像して前記複数枚のスペックル画像を出力する。
撮像装置１５は、コヒーレント光が照射された撮像対象Ｏから得られる散乱光を、撮像素子を用いて複数回撮像して前記複数枚のスペックル画像を出力する。
撮像装置１５が行う撮像方法は、本技術の効果が損なわれない限り特に限定されず、公知の撮像方法を１種又は２種以上選択して、自由に組み合わせて用いることができる。例えば、ＣＣＤ（Charge Coupled Device）やＣＭＯＳ（Complementary Metal Oxide Semiconductor）センサー等の撮像素子を用いた撮像方法を挙げることができる。
本開示では、上記撮像装置１５の露光時間は１０ｍｓ以下に設定されるが、１ｍｓ以下に設定されることが好ましく、１００μｓ程度に設定されることがより好ましい。
具体的な方法として、本開示では、レーザーなどコヒーレントな光源を用いて被写体（生体など）を照明し、ＣＣＤやＣＭＯＳなどの撮像素子を用いて１０ｍｓ以下の短い露光時間で複数枚（２枚以上）の撮像が行われる。露光時間は、散乱体の動きによりスペックルコントラストが低下しないような露光時間を用いてもよい。この場合、複数枚の画像を平均化することで、露光時間を延長させる効果が得られる。露光時間の最適値が未知の場合、十分短くした露光時間で撮影すれば、後処理によって露光時間を調節することも可能である。例えば、一定の露光時間（１０ｍｓ）で撮影された画像の場合、２フレーム足し合わせることで２０ｍｓ相当の画像とし、３フレーム足し合わせることで３０ｍｓ相当の画像とすることができる（図５参照）。図５は、スペックルコントラストの算出概念を示した図面代用写真である。このように、一定の露光時間（１０ｍｓ）を２、３フレーム分だけ足し合わせる数学的な処理によって得られる積算露光時間（２０ｍｓ、３０ｍｓ）は、前記流体速度算出部１３における血流速度の算出に用いられる。
（３）表示装置１６
表示装置１６は、前記輝度積算部１１で積算したスペックル積算画像などの画像を表示する。また、この表示装置１６においては、前記流体速度算出部１３で算出された血流速度をマッピングし、血流速度の分布を更に表示することも可能である。
表示装置１６は、前記輝度積算部１１で積算したスペックル積算画像などの画像を表示する。また、この表示装置１６においては、前記流体速度算出部１３で算出された血流速度をマッピングし、血流速度の分布を更に表示することも可能である。
（４）記憶装置１７
記憶装置１７は、前記輝度積算部１１で積算したスペックル積算画像、前記コントラスト算出部１２によって算出されたスペックルコントラストＫなどを記憶する。また、この記憶装置１７においては、血流速度の分布を更に記憶することも可能である。
記憶装置１７は、前記輝度積算部１１で積算したスペックル積算画像、前記コントラスト算出部１２によって算出されたスペックルコントラストＫなどを記憶する。また、この記憶装置１７においては、血流速度の分布を更に記憶することも可能である。
（５）撮像対象Ｏ
本開示に係るスペックルイメージングシステム１０は、様々なものを撮像対象とすることができるが、例えば、流体を含むものを撮像対象とするイメージングに、好適に用いることができる。より具体的には、撮像対象Ｏを生体とし、流体としては血液を挙げることができる。例えば、本開示に係るスペックルイメージングシステム１０を、手術用顕微鏡や手術用内視鏡などに適用すれば、血管の位置を確認しながら手術を行うことが可能である。そのため、より安全で高精度な手術を行うことができ、医療技術の更なる発展にも貢献することができる。
本開示に係るスペックルイメージングシステム１０は、様々なものを撮像対象とすることができるが、例えば、流体を含むものを撮像対象とするイメージングに、好適に用いることができる。より具体的には、撮像対象Ｏを生体とし、流体としては血液を挙げることができる。例えば、本開示に係るスペックルイメージングシステム１０を、手術用顕微鏡や手術用内視鏡などに適用すれば、血管の位置を確認しながら手術を行うことが可能である。そのため、より安全で高精度な手術を行うことができ、医療技術の更なる発展にも貢献することができる。
（６）スペックルイメージングの第１フロー例
図６は、本開示の第２実施形態に係るスペックルイメージングシステム１０を用いて行うスペックルイメージングの第１フロー例を示すフローチャートである。この第１フロー例では、予め測定しておいた相関時間τcと、流体速度算出部１３で求めた相関時間τcとを比較することにより血流速度が算出される。以下、時系列に沿って、第１フロー例を説明する。
図６は、本開示の第２実施形態に係るスペックルイメージングシステム１０を用いて行うスペックルイメージングの第１フロー例を示すフローチャートである。この第１フロー例では、予め測定しておいた相関時間τcと、流体速度算出部１３で求めた相関時間τcとを比較することにより血流速度が算出される。以下、時系列に沿って、第１フロー例を説明する。
（ａ）露光時間及びフレームレートの設定
まず、ステップＳＴ２０１では、撮像装置１５において、露光時間及びフレームレートが設定される。
まず、ステップＳＴ２０１では、撮像装置１５において、露光時間及びフレームレートが設定される。
（ｂ）積算枚数の決定
次に、ステップＳＴ２０２では、撮像装置１５において、ステップＳＴ２０１で設定された露光時間及びフレームレートに基づいて、輝度を積算すべき画像の積算枚数が決定される。
次に、ステップＳＴ２０２では、撮像装置１５において、ステップＳＴ２０１で設定された露光時間及びフレームレートに基づいて、輝度を積算すべき画像の積算枚数が決定される。
（ｃ）スペックル画像の撮像
次に、ステップＳＴ２０３において、撮像装置１５は、撮像素子を用いて、ステップＳＴ２０２で決定された積算枚数に相当する回数分を、撮像素子を用いて撮像して複数枚のスペックル画像を情報処理装置１に出力する。
次に、ステップＳＴ２０３において、撮像装置１５は、撮像素子を用いて、ステップＳＴ２０２で決定された積算枚数に相当する回数分を、撮像素子を用いて撮像して複数枚のスペックル画像を情報処理装置１に出力する。
（ｄ）スペックル画像の記憶
次に、ステップＳＴ２０４では、記憶装置１７において複数枚のスペックル画像が記録される。
次に、ステップＳＴ２０４では、記憶装置１７において複数枚のスペックル画像が記録される。
（ｅ）スペックル画像の表示
次に、ステップＳＴ２０５では、表示装置１６において複数枚のスペックル画像が表示される。
次に、ステップＳＴ２０５では、表示装置１６において複数枚のスペックル画像が表示される。
なお、上記ステップＳＴ２０４、ＳＴ２０５は、いずれか一方のみが実施されてもよい。また、スペックル画像の記憶処理（ステップＳＴ２０４）は、スペックル画像の表示処理（ステップＳＴ２０５）の後に行われてもよい。なお、スペックル画像の記憶は、スペックル画像の解析を後に行うために行われる。また、スペックル画像の表示は、撮影したスペックル画像を確認するために行われる。
（ｆ）輝度の積算
次に、ステップＳＴ２０６において、情報処理装置１の輝度積算部１１は、複数枚のスペックル画像の輝度を積算する。
次に、ステップＳＴ２０６において、情報処理装置１の輝度積算部１１は、複数枚のスペックル画像の輝度を積算する。
（ｇ）スペックルパターンのコントラストの算出
次に、ステップＳＴ２０７において、情報処理装置１のコントラスト算出部１２は、輝度積算部１１で積算したスペックル積算画像に基づいて、スペックルパターンのコントラストＫを算出する。なお、コントラストＫは、記憶装置１７に記憶されてもよい。
次に、ステップＳＴ２０７において、情報処理装置１のコントラスト算出部１２は、輝度積算部１１で積算したスペックル積算画像に基づいて、スペックルパターンのコントラストＫを算出する。なお、コントラストＫは、記憶装置１７に記憶されてもよい。
（ｈ）相関時間の算出
次に、情報処理装置１の流体速度算出部１３は、積算露光時間Ｔと、スペックルパターンのコントラストＫとに基づいて、撮像対象の流体速度を算出する。ステップＳＴ２０８において、流体速度算出部１３は、積算露光時間ＴとコントラストＫとに基づいて、相関時間τcを算出する。具体的には、積算露光時間ＴとコントラストＫとを上記の式（２）に代入し、相関時間τcを算出する。なお、相関時間τcは、記憶装置１７に記憶されてもよい。
次に、情報処理装置１の流体速度算出部１３は、積算露光時間Ｔと、スペックルパターンのコントラストＫとに基づいて、撮像対象の流体速度を算出する。ステップＳＴ２０８において、流体速度算出部１３は、積算露光時間ＴとコントラストＫとに基づいて、相関時間τcを算出する。具体的には、積算露光時間ＴとコントラストＫとを上記の式（２）に代入し、相関時間τcを算出する。なお、相関時間τcは、記憶装置１７に記憶されてもよい。
（ｉ）流体速度の算出
次に、ステップＳＴ２０９において、流体速度算出部１３は、算出した相関時間τcと、予め定められた相関時間τcとを比較することによって、流体速度を算出する。なお、コントラストＫは、記憶装置１７に記憶されてもよく、表示装置１６に表示されてもよい。
次に、ステップＳＴ２０９において、流体速度算出部１３は、算出した相関時間τcと、予め定められた相関時間τcとを比較することによって、流体速度を算出する。なお、コントラストＫは、記憶装置１７に記憶されてもよく、表示装置１６に表示されてもよい。
（ｊ）流体速度の分布の記憶
次に、ステップＳＴ２１０では、記憶装置１７において、前記流体速度算出部１３で算出された流体速度をマッピングすることにより得られた流体速度の分布が記憶される。
次に、ステップＳＴ２１０では、記憶装置１７において、前記流体速度算出部１３で算出された流体速度をマッピングすることにより得られた流体速度の分布が記憶される。
そして、ステップＳＴ２１１では、表示装置１６において、マッピングすることにより得られた流体速度の分布が表示される。
（７）スペックルイメージングの第２フロー例
図７は、本開示の第２実施形態に係るスペックルイメージングシステム１０を用いて行うスペックルイメージングの第２フロー例を示すフローチャートである。この第２フロー例では、輝度の積算を行うべき目的の枚数に到達するまで、スペックル画像の撮像が反復して実行される。以下、時系列に沿って、第２フロー例を説明する。
図７は、本開示の第２実施形態に係るスペックルイメージングシステム１０を用いて行うスペックルイメージングの第２フロー例を示すフローチャートである。この第２フロー例では、輝度の積算を行うべき目的の枚数に到達するまで、スペックル画像の撮像が反復して実行される。以下、時系列に沿って、第２フロー例を説明する。
（ａ）露光時間の設定
まず、ステップＳＴ３０１では、撮像装置１５において、露光時間が設定される。
まず、ステップＳＴ３０１では、撮像装置１５において、露光時間が設定される。
（ｂ）積算枚数の決定
次に、ステップＳＴ３０２では、撮像装置１５において、ステップＳＴ３０１で設定された露光時間に基づいて、輝度を積算すべき画像の積算枚数が決定される。
次に、ステップＳＴ３０２では、撮像装置１５において、ステップＳＴ３０１で設定された露光時間に基づいて、輝度を積算すべき画像の積算枚数が決定される。
（ｃ）スペックル画像の撮像
次に、ステップＳＴ３０３では、撮像装置１５は、撮像素子を用いて撮像してスペックル画像を情報処理装置１に出力する。
次に、ステップＳＴ３０３では、撮像装置１５は、撮像素子を用いて撮像してスペックル画像を情報処理装置１に出力する。
（ｄ）スペックル画像の記憶
次に、ステップＳＴ３０４では、記憶装置１７において複数枚のスペックル画像が記憶される。なお、記憶装置１７には、撮像回数（積算枚数）が記憶されてもよい。
次に、ステップＳＴ３０４では、記憶装置１７において複数枚のスペックル画像が記憶される。なお、記憶装置１７には、撮像回数（積算枚数）が記憶されてもよい。
（ｅ）スペックル画像の表示
次に、ステップＳＴ３０５では、表示装置１６においてスペックル画像が表示される。
次に、ステップＳＴ３０５では、表示装置１６においてスペックル画像が表示される。
なお、上記ステップＳＴ３０４、ＳＴ３０５は、いずれか一方のみが実施されてもよい。また、スペックル画像の記憶処理（ステップＳＴ３０４）は、スペックル画像の表示処理（ステップＳＴ３０５）の後に行われてもよい。なお、スペックル画像の記憶は、スペックル画像の解析を後に行うために行われる。また、スペックル画像の表示は、撮影したスペックル画像を確認するために行われる。
（ｆ）積算枚数のカウント
次に、ステップＳＴ３０６では、撮像装置１５において、撮像回数（積算枚数）のカウントが行われる。目的の撮像回数に到達していない場合は、ステップＳＴ３０３へ戻り、撮像が繰り返される。撮像回数が目的の回数に到達した場合は、ステップＳＴ３０７へ進む。
次に、ステップＳＴ３０６では、撮像装置１５において、撮像回数（積算枚数）のカウントが行われる。目的の撮像回数に到達していない場合は、ステップＳＴ３０３へ戻り、撮像が繰り返される。撮像回数が目的の回数に到達した場合は、ステップＳＴ３０７へ進む。
（ｇ）輝度の積算
次に、ステップＳＴ３０７において、情報処理装置１の輝度積算部１１は、複数枚のスペックル画像の輝度を積算する。
次に、ステップＳＴ３０７において、情報処理装置１の輝度積算部１１は、複数枚のスペックル画像の輝度を積算する。
（ｈ）スペックルパターンのコントラストの算出
次に、ステップＳＴ３０８において、情報処理装置１のコントラスト算出部１２は、輝度積算部１１で積算したスペックル積算画像に基づいて、スペックルパターンのコントラストＫを算出する。なお、コントラストＫは、記憶装置１７に記憶されてもよい。
次に、ステップＳＴ３０８において、情報処理装置１のコントラスト算出部１２は、輝度積算部１１で積算したスペックル積算画像に基づいて、スペックルパターンのコントラストＫを算出する。なお、コントラストＫは、記憶装置１７に記憶されてもよい。
（ｉ）スペックル積算画像の記憶
次に、ステップＳＴ３０９では、記憶装置１７においてスペックル積算画像が記憶される。
次に、ステップＳＴ３０９では、記憶装置１７においてスペックル積算画像が記憶される。
（ｊ）スペックル積算画像の表示
そして、ステップＳＴ３１０では、表示装置１６においてスペックル積算画像が表示される。
そして、ステップＳＴ３１０では、表示装置１６においてスペックル積算画像が表示される。
なお、上記ステップＳＴ３０９、ＳＴ３１０は、いずれか一方のみが実施されてもよい。また、スペックル画像の記憶処理（ステップＳＴ３０９）は、スペックル画像の表示処理（ステップＳＴ３１０）の後に行われてもよい。
＜３．第３実施形態（情報処理方法）＞
本開示に係る情報処理方法は、大別して、輝度積算ステップ、及び、スペックルコントラスト算出ステップを少なくとも行う方法である。また、必要に応じて、流体速度算出ステップ、記憶ステップ及び表示ステップなどを更に行うことも可能である。なお、輝度積算ステップ、スペックルコントラスト算出ステップ、流体速度算出ステップ、表示ステップ及び記憶ステップは、前述した本開示に係るスペックルイメージングシステム１０の輝度積算部１１、コントラスト算出部１２、流体速度算出部１３、表示装置１６及び記憶装置１７が、それぞれ行う方法と同一であるため、ここでは説明を省略する。
本開示に係る情報処理方法は、大別して、輝度積算ステップ、及び、スペックルコントラスト算出ステップを少なくとも行う方法である。また、必要に応じて、流体速度算出ステップ、記憶ステップ及び表示ステップなどを更に行うことも可能である。なお、輝度積算ステップ、スペックルコントラスト算出ステップ、流体速度算出ステップ、表示ステップ及び記憶ステップは、前述した本開示に係るスペックルイメージングシステム１０の輝度積算部１１、コントラスト算出部１２、流体速度算出部１３、表示装置１６及び記憶装置１７が、それぞれ行う方法と同一であるため、ここでは説明を省略する。
＜４．ハードウェア構成＞
上述した第１実施形態に係る情報処理装置１の処理は、ソフトウェアと、以下に説明するハードウェアとの協働により実現される。
上述した第１実施形態に係る情報処理装置１の処理は、ソフトウェアと、以下に説明するハードウェアとの協働により実現される。
図８は、本開示の第１実施形態に係る情報処理装置１のハードウェア構成を示した説明図である。図８に示したように、情報処理装置１は、ＣＰＵ（Central Processing Unit）１０１と、ＲＯＭ（Read Only Memory）１０２と、ＲＡＭ（Random Access Memory）１０３と、ブリッジ１０４と、バス１０５と、インターフェース１０６と、入力装置１０７と、出力装置１０８と、ストレージ装置１０９と、接続ポート１１０と、通信装置１１１とを備える。
ＣＰＵ１０１は、情報処理装置として機能し、各種プログラムと協働して情報処理装置１内の輝度積算部１１、コントラスト算出部１２及び流体速度算出部１３の動作を実現する。また、ＣＰＵ１０１は、マイクロプロセッサであってもよい。ＲＯＭ１０２は、ＣＰＵ１０１が使用するプログラムまたは演算パラメータ等を記憶する。ＲＡＭ１０３は、ＣＰＵ１０１の実行において使用するプログラムまたは実行において適宜変化するパラメータ等を一時記憶する。ＲＯＭ１０２およびＲＡＭ１０３により、情報処理装置１内のメモリの一部を実現する。ＣＰＵ１０１、ＲＯＭ１０２およびＲＡＭ１０３は、ＣＰＵバスなどから構成される内部バスにより相互に接続されている。
入力装置１０７は、タッチパネル、ボタン、マイクロフォン、スイッチおよびレバーなどユーザが情報を入力するための入力手段、およびユーザによる入力に基づいて入力信号を生成し、ＣＰＵ１０１に出力する入力制御回路などから構成されている。情報処理装置１のユーザは、入力装置１０７を操作することにより、情報処理装置１に対して各種のデータを入力したり処理動作を指示したりすることができる。
出力装置１０８は、例えば、液晶ディスプレイ（ＬＣＤ）装置、ＯＬＥＤ（Organic Light Emitting Diode）装置、ランプなどの装置への出力を行う。さらに、ユーザフレンドリーの観点から、出力装置１０８は、スピーカおよびヘッドフォンなどの音声出力を行ってもよい。
ストレージ装置１０９は、データ格納用の装置である。ユーザフレンドリーの観点から、ストレージ装置１０９は、記憶媒体、記憶媒体にデータを記録する記録装置、記憶媒体からデータを読み出す読出し装置および記憶媒体に記録されたデータを削除する削除装置等を含んでもよい。ストレージ装置１０９は、ＣＰＵ１０１が実行するプログラムや各種データを格納する。
接続ポート１１０は、例えば、情報処理装置１の外部の装置または周辺機器と接続するためのバスである。また、ユーザフレンドリーの観点から、接続ポート１１０は、ＵＳＢ（Universal Serial Bus）であってもよい。
通信装置１１１は、例えば、ネットワークに接続するための通信デバイスで構成された通信インターフェースである。また、ユーザフレンドリーの観点から、通信装置１１１は、赤外線通信対応装置であっても、無線ＬＡＮ（Local Area Network）対応通信装置であっても、ＬＴＥ（Long Term Evolution）対応通信装置であっても、有線による通信を行うワイヤー通信装置であってもよい。
以上、上記第１〜第３実施形態では、流体速度を測定する前提としてスペックルパターンのコントラストを簡便かつ効率的に取得することができる。
また、上記第１〜第３実施形態では、コヒーレントな光源の強度を変えることなく（一定の強度）で照明し、一定の露光時間で撮影することが可能である。これにより、血流速度を測定する際に、レーザーを減衰させるための外部制御を必要とせず、装置を簡略化することができ、コストや技術的な難易度を低減させることができる。
また、上記第１〜第３実施形態では、血流速度をマッピングした画像を取得することができる。従来、スペックルコントラストによって流速を表示する方法としてMulti-Exposure Speckle Imaging法などが提案されているが、この方法では、様々な露光時間による撮像が必要であり、それに応じた強度のコヒーレントな光源による照明が必要であった。露光時間の調節については容易であるが、コヒーレントな光源の照明強度を幅広く変えることは困難であった。上記第１〜第３実施形態では、従来のMulti-exposure speckle imaging法より少ない露光時間のパターンによって血流速度を算出できるようしたので、血流速度を測定する際に、従来のMulti-exposure speckle imaging法のようにコヒーレントな光源の照明強度を幅広く変える外部制御を必要とせず、装置を簡略化することができ、コストや技術的な難易度を低減させることができる。
散乱体の移動（流れ）をスペックルコントラストで観測する場合、ある程度、露光時間を長くする必要がある。例えばデジタルイメージャを用いる場合、散乱体が露光時間中に１画素（スペックルパターンが変わる程度）以上、移動する必要がある。従来手法では、この露光時間を測定対象物に最適な露光時間を設定しないと、スペックルコントラストを算出できなかった。本実施形態のように、後処理によって露光時間を調節できるが、後処理によって露光時間を短くすることは不可能である。後処理で最適な露光時間を設定するためには、短い露光時間を実現する必要がある。
また、上記第１〜第３実施形態では、スペックルコントラストを算出するために撮影後の露光時間の最適化をおこなうことが可能である。さらに、スペックル血流イメージでは描画できないような遅いあるいは止まった（但しブラウン運動はしている）などの流速についても描画が可能である（流速感度のダイナミックレンジを広げることができる）。
ここでの最適化とは、例えば、以下の用途（Ａ）〜（Ｃ）に最適な血流イメージングをすることをいう。つまり、（Ａ）血流が正常に流れている部分と、流れていない部分とを検出したり、（Ｂ）血管に狭窄（血流が流れていない箇所）が無くても何かが詰まっていて流れていない部分を検出したり、（Ｃ）狭窄のレベル評価を実現可能（ex.画像を何枚重ねたら狭窄を発見できるか等）を行うに当たり、血管の血流に応じて最適な血流イメージが得られる。
また、本開示は、以下のような構成をとることもできる。
（１）
コヒーレント光が照射された撮像対象から得られる散乱光を、撮像素子を用いて複数回撮像することにより得られた複数枚のスペックル画像の輝度を積算する輝度積算部と、
前記輝度積算部で積算したスペックル積算画像に基づいて、スペックルパターンのコントラストを算出するコントラスト算出部と、を備える情報処理装置。
（２）
前記複数枚のスペックル画像の露光時間を積算した積算露光時間と、前記コントラスト算出部で算出したスペックルパターンのコントラストとに基づいて、前記撮像対象の流体速度を算出する流体速度算出部を更に備える上記（１）に記載の情報処理装置。
（３）
前記複数枚のスペックル画像は、１０ｍｓ以下の露光時間で撮像された画像である上記（１）又は（２）に記載の情報処理装置。
（４）
前記コントラスト算出部は、前記スペックル積算画像の局所的な画像領域における輝度の分散値及び平均値を求め、求めた輝度の分散値及び平均値に基づいて、スペックルパターンのコントラストを算出する上記（１）から（３）のいずれかに記載の情報処理装置。
（５）
前記コントラスト算出部は、前記求めた輝度の分散値の平方根をとることで輝度の標準偏差を算出し、スペックルパターンのコントラスト、輝度の標準偏差、及び、輝度の平均値の関係を示す式（１）に前記求めた輝度の標準偏差及び平均値を代入することによって、スペックルパターンのコントラストを算出する上記（４）に記載の情報処理装置。
ただし、Ｋはスペックルパターンのコントラスト、σは輝度Ｉの標準偏差、＜Ｉ＞は輝度Ｉの平均値である。
（６）
前記流体速度算出部は、前記積算露光時間と前記スペックルパターンのコントラストに基づいて相関時間を求め、求めた相関時間と、予め定められた相関時間とを比較することによって、前記流体速度を算出する上記（２）に記載の情報処理装置。
（７）
前記流体速度算出部は、露光時間、スペックルパターンのコントラスト、及び、相関時間の関係を示す式（２）に前記積算露光時間と前記スペックルパターンのコントラストとを代入することによって、相関時間を求める上記（６）に記載の情報処理装置。
（８）
前記流体速度算出部は、露光時間、スペックルパターンのコントラスト、及び、相関時間の関係を示す式（３）に前記積算露光時間と前記スペックルパターンのコントラストとを代入することによって、相関時間を求める上記（６）に記載の情報処理装置。
（９）
表示部に画像を表示させる表示制御部を更に備える上記（２）に記載の情報処理装置。
（１０）
前記表示制御部は、前記流体速度算出部で算出された流体速度をマッピングし、流体速度の分布を前記表示部に更に表示させる上記（９）に記載の情報処理装置。
（１１）
前記流体速度は、血管内における血流速度である上記（２）に記載の情報処理装置。
（１２）
上記（１）乃至（１１）のいずれかに記載の情報処理装置と、
撮像対象にコヒーレント光を照射する光源と、
コヒーレント光が照射された撮像対象から得られる散乱光を、撮像素子を用いて複数回撮像して複数枚のスペックル画像を出力する撮像装置と、
画像を表示する表示装置と、を備えるスペックルイメージングシステム。
（１３）
前記表示装置は、前記流体速度算出部で算出された流体速度をマッピングし、流体速度の分布を更に表示する上記（１２）に記載のスペックルイメージングシステム。
（１４）
コヒーレント光が照射された撮像対象から得られる散乱光を、撮像素子を用いて複数回撮像することにより得られた複数枚のスペックル画像の輝度を積算する輝度積算ステップと、
前記輝度積算ステップで積算したスペックル積算画像に基づいて、スペックルパターンのコントラストを算出するコントラスト算出ステップと、を少なくとも行う、情報処理方法。
（１）
コヒーレント光が照射された撮像対象から得られる散乱光を、撮像素子を用いて複数回撮像することにより得られた複数枚のスペックル画像の輝度を積算する輝度積算部と、
前記輝度積算部で積算したスペックル積算画像に基づいて、スペックルパターンのコントラストを算出するコントラスト算出部と、を備える情報処理装置。
（２）
前記複数枚のスペックル画像の露光時間を積算した積算露光時間と、前記コントラスト算出部で算出したスペックルパターンのコントラストとに基づいて、前記撮像対象の流体速度を算出する流体速度算出部を更に備える上記（１）に記載の情報処理装置。
（３）
前記複数枚のスペックル画像は、１０ｍｓ以下の露光時間で撮像された画像である上記（１）又は（２）に記載の情報処理装置。
（４）
前記コントラスト算出部は、前記スペックル積算画像の局所的な画像領域における輝度の分散値及び平均値を求め、求めた輝度の分散値及び平均値に基づいて、スペックルパターンのコントラストを算出する上記（１）から（３）のいずれかに記載の情報処理装置。
（５）
前記コントラスト算出部は、前記求めた輝度の分散値の平方根をとることで輝度の標準偏差を算出し、スペックルパターンのコントラスト、輝度の標準偏差、及び、輝度の平均値の関係を示す式（１）に前記求めた輝度の標準偏差及び平均値を代入することによって、スペックルパターンのコントラストを算出する上記（４）に記載の情報処理装置。
ただし、Ｋはスペックルパターンのコントラスト、σは輝度Ｉの標準偏差、＜Ｉ＞は輝度Ｉの平均値である。
（６）
前記流体速度算出部は、前記積算露光時間と前記スペックルパターンのコントラストに基づいて相関時間を求め、求めた相関時間と、予め定められた相関時間とを比較することによって、前記流体速度を算出する上記（２）に記載の情報処理装置。
（７）
前記流体速度算出部は、露光時間、スペックルパターンのコントラスト、及び、相関時間の関係を示す式（２）に前記積算露光時間と前記スペックルパターンのコントラストとを代入することによって、相関時間を求める上記（６）に記載の情報処理装置。
（８）
前記流体速度算出部は、露光時間、スペックルパターンのコントラスト、及び、相関時間の関係を示す式（３）に前記積算露光時間と前記スペックルパターンのコントラストとを代入することによって、相関時間を求める上記（６）に記載の情報処理装置。
（９）
表示部に画像を表示させる表示制御部を更に備える上記（２）に記載の情報処理装置。
（１０）
前記表示制御部は、前記流体速度算出部で算出された流体速度をマッピングし、流体速度の分布を前記表示部に更に表示させる上記（９）に記載の情報処理装置。
（１１）
前記流体速度は、血管内における血流速度である上記（２）に記載の情報処理装置。
（１２）
上記（１）乃至（１１）のいずれかに記載の情報処理装置と、
撮像対象にコヒーレント光を照射する光源と、
コヒーレント光が照射された撮像対象から得られる散乱光を、撮像素子を用いて複数回撮像して複数枚のスペックル画像を出力する撮像装置と、
画像を表示する表示装置と、を備えるスペックルイメージングシステム。
（１３）
前記表示装置は、前記流体速度算出部で算出された流体速度をマッピングし、流体速度の分布を更に表示する上記（１２）に記載のスペックルイメージングシステム。
（１４）
コヒーレント光が照射された撮像対象から得られる散乱光を、撮像素子を用いて複数回撮像することにより得られた複数枚のスペックル画像の輝度を積算する輝度積算ステップと、
前記輝度積算ステップで積算したスペックル積算画像に基づいて、スペックルパターンのコントラストを算出するコントラスト算出ステップと、を少なくとも行う、情報処理方法。
なお、本明細書に記載された効果はあくまでも例示であって限定されるものではなく、また他の効果があってもよい。
以下、本開示の実験例により、本開示の効果について具体的に説明する。
＜第１実験例＞
第１実験例では、波長８２０ｎｍのコヒーレントなレーザー光源により均一に照明された生体を、ソニー製のグローバルシャッターのＣＭＯＳイメージャーを用いて、フレームレート１２０ｆｐｓ、露光時間１００ｕｓで撮像した。なお、輝度Ｉの取得はソニー製のグローバルシャッターのＣＭＯＳイメージャーによって実現される。１、１０、１００、１０００枚の時系列のスペックル画像を積算することにより、１００ｕｓ、１ｍｓ、１０ｍｓ、１００ｍｓの露光時間に相当するスペックル積算画像を求めた。それぞれのスペックル積算画像に基づき、上記式（１）を用い、局所的なスペックルコントラストＫを求めた。当該スペックルコントラストＫと、積算露光時間Ｔとでプロットし、上記式（２）でフィッティングを行い、相関時間τcを算出した。当該算出により得られた相関時間τcを、既知の流速によって算出した相関時間τcと比較することにより、局所的な部位の血流速度を求め、求めた血流速度をマッピングした（図９参照）。図９は、血流速度をマッピングしたマップ図を示す図面代用写真である。マップ図の色を見れば、血流が正常に流れている部分と、流れていない部分とを検出したり、血流速度を把握したりすることが可能である。
第１実験例では、波長８２０ｎｍのコヒーレントなレーザー光源により均一に照明された生体を、ソニー製のグローバルシャッターのＣＭＯＳイメージャーを用いて、フレームレート１２０ｆｐｓ、露光時間１００ｕｓで撮像した。なお、輝度Ｉの取得はソニー製のグローバルシャッターのＣＭＯＳイメージャーによって実現される。１、１０、１００、１０００枚の時系列のスペックル画像を積算することにより、１００ｕｓ、１ｍｓ、１０ｍｓ、１００ｍｓの露光時間に相当するスペックル積算画像を求めた。それぞれのスペックル積算画像に基づき、上記式（１）を用い、局所的なスペックルコントラストＫを求めた。当該スペックルコントラストＫと、積算露光時間Ｔとでプロットし、上記式（２）でフィッティングを行い、相関時間τcを算出した。当該算出により得られた相関時間τcを、既知の流速によって算出した相関時間τcと比較することにより、局所的な部位の血流速度を求め、求めた血流速度をマッピングした（図９参照）。図９は、血流速度をマッピングしたマップ図を示す図面代用写真である。マップ図の色を見れば、血流が正常に流れている部分と、流れていない部分とを検出したり、血流速度を把握したりすることが可能である。
＜第２実験例＞
第２実験例では、波長８２０ｎｍのコヒーレントなレーザー光源により均一に照明された生体を模したファントムに血液を送液し、ソニー製のグローバルシャッターのＣＭＯＳイメージャーを用いて、フレームレート３０ｆｐｓ、露光時間５００ｕｓで撮像した。１枚（比較例）、２枚、・・・、１０枚の時系列のスペックル画像を積算することにより、５００ｕｓ（比較例）、１ｍｓ、・・・、５ｍｓの積算露光時間に相当する各スペックル積算画像を求めた（図１０参照）。そして、図１０に示す各スペックル積算画像の中から、血流速度の測定対象物に応じて最適なスペックル積算画像を選択することができた。
第２実験例では、波長８２０ｎｍのコヒーレントなレーザー光源により均一に照明された生体を模したファントムに血液を送液し、ソニー製のグローバルシャッターのＣＭＯＳイメージャーを用いて、フレームレート３０ｆｐｓ、露光時間５００ｕｓで撮像した。１枚（比較例）、２枚、・・・、１０枚の時系列のスペックル画像を積算することにより、５００ｕｓ（比較例）、１ｍｓ、・・・、５ｍｓの積算露光時間に相当する各スペックル積算画像を求めた（図１０参照）。そして、図１０に示す各スペックル積算画像の中から、血流速度の測定対象物に応じて最適なスペックル積算画像を選択することができた。
スペックルコントラストＫは、撮像対象のエリアごとの移動速度に依存する。第２実施例では、血流速度が１ｍｍ/ｓｅｃの場合、５００ｕｓ、１ｍｓ、５ｍｓの積算露光時間を得るに当たり、スペックルコントラストＫの値を、それぞれ、０．６、０．５、０．４程度に設定した。図１１は、５００ｕｓ、１ｍｓ、５ｍｓの積算露光時間に相当する各スペックル積算画像と、実際の露光時間１ｍｓ、５ｍｓで撮像した画像と、を示している。図１１に示されるように、１ｍｓ、５ｍｓの積算露光時間に相当する各スペックル積算画像は、実際の露光時間１ｍｓ、５ｍｓで撮像した画像とほぼ等価の画像が得られることが分かった。
１ 情報処理装置
１０ スペックルイメージングシステム
１１ 輝度積算部
１２ コントラスト算出部
１３ 流体速度算出部
１４ 光源
１５ 撮像装置
１６ 表示装置
１７ 記憶装置
１０１ ＣＰＵ
１０２ ＲＯＭ
１０３ ＲＡＭ
１０４ ブリッジ
１０５ バス
１０６ インターフェース
１０７ 入力装置
１０８ 出力装置
１０９ ストレージ装置
１１０ 接続ポート
１１１ 通信装置
１０ スペックルイメージングシステム
１１ 輝度積算部
１２ コントラスト算出部
１３ 流体速度算出部
１４ 光源
１５ 撮像装置
１６ 表示装置
１７ 記憶装置
１０１ ＣＰＵ
１０２ ＲＯＭ
１０３ ＲＡＭ
１０４ ブリッジ
１０５ バス
１０６ インターフェース
１０７ 入力装置
１０８ 出力装置
１０９ ストレージ装置
１１０ 接続ポート
１１１ 通信装置
Claims (13)
- コヒーレント光が照射された撮像対象から得られる散乱光を、撮像素子を用いて複数回撮像することにより得られた複数枚のスペックル画像の輝度を積算する輝度積算部と、
前記輝度積算部で積算したスペックル積算画像に基づいて、スペックルパターンのコントラストを算出するコントラスト算出部と、を備え、
前記コントラスト算出部は、前記スペックル積算画像の局所的な画像領域における輝度の分散値及び平均値を求め、求めた輝度の分散値及び平均値に基づいて、スペックルパターンのコントラストを算出する、情報処理装置。 - 前記複数枚のスペックル画像の露光時間を積算した積算露光時間と、前記コントラスト算出部で算出したスペックルパターンのコントラストとに基づいて、前記撮像対象の流体速度を算出する流体速度算出部を更に備える請求項１に記載の情報処理装置。
- 前記複数枚のスペックル画像は、１０ｍｓ以下の露光時間で撮像された画像である請求項１又は２に記載の情報処理装置。
- 前記コントラスト算出部は、前記求めた輝度の分散値の平方根をとることで輝度の標準偏差を算出し、スペックルパターンのコントラスト、輝度の標準偏差、及び、輝度の平均値の関係を示す式（１）に前記求めた輝度の標準偏差及び平均値を代入することによって、スペックルパターンのコントラストを算出する請求項１から３のいずれか一項に記載の情報処理装置。
ただし、Ｋはスペックルパターンのコントラスト、σは輝度Ｉの標準偏差、＜Ｉ＞は輝度Ｉの平均値である。 - 前記流体速度算出部は、前記積算露光時間と前記スペックルパターンのコントラストに基づいて相関時間を求め、求めた相関時間と、予め定められた相関時間とを比較することによって、前記流体速度を算出する請求項２に記載の情報処理装置。
- 前記流体速度算出部は、露光時間、スペックルパターンのコントラスト、及び、相関時間の関係を示す式（２）に前記積算露光時間と前記スペックルパターンのコントラストとを代入することによって、相関時間を求める請求項５に記載の情報処理装置。
- 前記流体速度算出部は、露光時間、スペックルパターンのコントラスト、及び、相関時間の関係を示す式（３）に前記積算露光時間と前記スペックルパターンのコントラストとを代入することによって、相関時間を求める請求項５に記載の情報処理装置。
- 表示部に画像を表示させる表示制御部を更に備える請求項２に記載の情報処理装置。
- 前記表示制御部は、前記流体速度算出部で算出された流体速度をマッピングし、流体速度の分布を前記表示部に更に表示させる請求項８に記載の情報処理装置。
- 前記流体速度は、血管内における血流速度である請求項２に記載の情報処理装置。
- 請求項１乃至１０のいずれか一項に記載の情報処理装置と、
撮像対象にコヒーレント光を照射する光源と、
コヒーレント光が照射された撮像対象から得られる散乱光を、撮像素子を用いて複数回撮像して複数枚のスペックル画像を出力する撮像装置と、
画像を表示する表示装置と、を備えるスペックルイメージングシステム。 - 前記表示装置は、前記流体速度算出部で算出された流体速度をマッピングし、流体速度の分布を更に表示する請求項１１に記載のスペックルイメージングシステム。
- コヒーレント光が照射された撮像対象から得られる散乱光を、撮像素子を用いて複数回撮像することにより得られた複数枚のスペックル画像の輝度を積算する輝度積算ステップと、
前記輝度積算ステップで積算したスペックル積算画像に基づいて、スペックルパターンのコントラストを算出するコントラスト算出ステップと、を少なくとも行い、
前記コントラスト算出ステップにおいて、前記スペックル積算画像の局所的な画像領域における輝度の分散値及び平均値を求め、求めた輝度の分散値及び平均値に基づいて、スペックルパターンのコントラストを算出する、情報処理方法。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2015238047 | 2015-12-04 | ||
| JP2015238047 | 2015-12-04 | ||
| PCT/JP2016/080906 WO2017094380A1 (ja) | 2015-12-04 | 2016-10-19 | 情報処理装置、スペックルイメージングシステム、及び情報処理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPWO2017094380A1 JPWO2017094380A1 (ja) | 2018-09-20 |
| JP6927046B2 true JP6927046B2 (ja) | 2021-08-25 |
Family
ID=58797028
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2017553697A Active JP6927046B2 (ja) | 2015-12-04 | 2016-10-19 | 情報処理装置、スペックルイメージングシステム、及び情報処理方法 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US20180344176A1 (ja) |
| EP (1) | EP3385723B1 (ja) |
| JP (1) | JP6927046B2 (ja) |
| CN (1) | CN108291925B (ja) |
| WO (1) | WO2017094380A1 (ja) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109314764B (zh) * | 2016-04-20 | 2021-02-05 | 雷瑟联合科技有限公司 | 用于校准和修正散斑对比流量计的系统和方法 |
| WO2018100603A1 (ja) * | 2016-11-29 | 2018-06-07 | パイオニア株式会社 | 計測装置 |
| US20210235968A1 (en) * | 2018-08-28 | 2021-08-05 | Sony Corporation | Medical system, information processing apparatus, and information processing method |
| JP2022508237A (ja) * | 2018-11-27 | 2022-01-19 | コンティンユーズ バイオメトリクス リミテッド | 生物医学的パラメータを遠隔監視するためのシステムおよび方法 |
| WO2020153194A1 (ja) * | 2019-01-23 | 2020-07-30 | ソニー株式会社 | 情報処理装置、情報処理方法及びプログラム |
| CN110072061B (zh) * | 2019-05-31 | 2021-02-09 | 努比亚技术有限公司 | 一种交互式拍摄方法、移动终端及存储介质 |
| JP2021058478A (ja) * | 2019-10-08 | 2021-04-15 | ソニー株式会社 | 画像処理装置、画像処理システム及び画像処理方法 |
| CN111402313B (zh) * | 2020-03-13 | 2022-11-04 | 合肥的卢深视科技有限公司 | 图像深度恢复方法和装置 |
| CN111795800B (zh) * | 2020-06-03 | 2023-06-06 | 杭州三泰检测技术有限公司 | 一种视觉散斑对比度测量装置及方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH10290791A (ja) * | 1997-04-18 | 1998-11-04 | Advance Co Ltd | レーザー血流計 |
| WO2010096447A2 (en) * | 2009-02-17 | 2010-08-26 | Board Of Regents, The University Of Texas System | Quantitative imaging with multi-exposure speckle imaging (mesi) |
| US9028421B2 (en) * | 2009-05-13 | 2015-05-12 | Kyushu Institute Of Technology | Blood flow image diagnosing device |
| WO2013049123A1 (en) * | 2011-09-26 | 2013-04-04 | The Johns Hopkins University | Anisotropic processing of laser speckle images |
| JP6076111B2 (ja) * | 2013-02-07 | 2017-02-08 | 浜松ホトニクス株式会社 | 塊状細胞評価方法および塊状細胞評価装置 |
| US9282905B2 (en) * | 2013-05-29 | 2016-03-15 | University Of Washington Through Its Center For Commercialization | Methods for laser speckle contrast imaging of blood perfusion |
-
2016
- 2016-10-19 WO PCT/JP2016/080906 patent/WO2017094380A1/ja active Application Filing
- 2016-10-19 CN CN201680069354.9A patent/CN108291925B/zh not_active Expired - Fee Related
- 2016-10-19 US US15/778,758 patent/US20180344176A1/en not_active Abandoned
- 2016-10-19 EP EP16870320.5A patent/EP3385723B1/en active Active
- 2016-10-19 JP JP2017553697A patent/JP6927046B2/ja active Active
Also Published As
| Publication number | Publication date |
|---|---|
| EP3385723B1 (en) | 2020-12-02 |
| WO2017094380A1 (ja) | 2017-06-08 |
| EP3385723A1 (en) | 2018-10-10 |
| EP3385723A4 (en) | 2018-12-05 |
| JPWO2017094380A1 (ja) | 2018-09-20 |
| CN108291925B (zh) | 2020-10-09 |
| US20180344176A1 (en) | 2018-12-06 |
| CN108291925A (zh) | 2018-07-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6927046B2 (ja) | 情報処理装置、スペックルイメージングシステム、及び情報処理方法 | |
| US8648918B2 (en) | Method and system for obtaining a point spread function using motion information | |
| US10375330B2 (en) | Systems and methods for surface topography acquisition using laser speckle | |
| CN103857335B (zh) | 激光散斑图像的各向异性处理 | |
| JP6431535B2 (ja) | 歯牙表面撮像装置のためのビデオベースの自動取得 | |
| JP4420909B2 (ja) | 撮像装置 | |
| JP5873619B2 (ja) | 三角測量法を用いた歯科用３ｄカメラによる３次元物体の光学的走査方法および装置 | |
| JP7178423B2 (ja) | 遠隔歯科医療画像のためのガイダンス方法及びシステム | |
| JP5714875B2 (ja) | 蛍光内視鏡装置 | |
| JPWO2011158498A1 (ja) | 撮像装置及び撮像方法 | |
| CN107517346B (zh) | 基于结构光的拍照方法、装置及移动设备 | |
| WO2017163542A1 (ja) | 画像解析装置、画像解析方法 | |
| WO2018230098A1 (ja) | 内視鏡システム、内視鏡システムの作動方法 | |
| US20160157736A1 (en) | System and Methods for Assessment of Relative Fluid Flow Using Laser Speckle Imaging | |
| US20200267305A1 (en) | Imaging apparatus, imaging method, and imaging system | |
| JP2014153890A5 (ja) | ||
| JP2019501720A (ja) | 視線追跡のためのシステムおよび装置 | |
| JP2023099679A (ja) | 内視鏡システム及び、データ作成方法 | |
| JP2012050617A (ja) | 画像取得方法および装置 | |
| US11050931B2 (en) | Control device and control method | |
| WO2019226664A1 (en) | Noninvasive three-dimensional fluorescence microscopy for skin disease detection | |
| US9808150B2 (en) | Image processing apparatus, control method therefor, and storage medium storing program | |
| WO2018203429A1 (ja) | 流速測定方法、流速測定装置及びプログラム | |
| US20220196394A1 (en) | Imaging device, control method therefor, measuring device, and storage medium | |
| JP2014230172A5 (ja) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190827 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190827 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20201013 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20201203 |
|
| 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: 20210706 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210719 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 6927046 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |