JP2022513847A - 高精度マルチアパーチャスペクトルイメージングのためのシステムおよび方法 - Google Patents

高精度マルチアパーチャスペクトルイメージングのためのシステムおよび方法 Download PDF

Info

Publication number
JP2022513847A
JP2022513847A JP2021533780A JP2021533780A JP2022513847A JP 2022513847 A JP2022513847 A JP 2022513847A JP 2021533780 A JP2021533780 A JP 2021533780A JP 2021533780 A JP2021533780 A JP 2021533780A JP 2022513847 A JP2022513847 A JP 2022513847A
Authority
JP
Japan
Prior art keywords
multispectral
light
imaging system
frequency band
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.)
Granted
Application number
JP2021533780A
Other languages
English (en)
Other versions
JP7186298B2 (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.)
Spectral MD Inc
Original Assignee
Spectral MD Inc
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 Spectral MD Inc filed Critical Spectral MD Inc
Publication of JP2022513847A publication Critical patent/JP2022513847A/ja
Priority to JP2022188817A priority Critical patent/JP7529752B2/ja
Priority to JP2022188833A priority patent/JP7529753B2/ja
Application granted granted Critical
Publication of JP7186298B2 publication Critical patent/JP7186298B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J3/2823Imaging spectrometer
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/02Details
    • G01J3/0202Mechanical elements; Supports for optical elements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/02Details
    • G01J3/0205Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows
    • G01J3/0229Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows using masks, aperture plates, spatial light modulators or spatial filters, e.g. reflective filters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/02Details
    • G01J3/10Arrangements of light sources specially adapted for spectrometry or colorimetry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J3/2803Investigating the spectrum using photoelectric array detector
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J3/30Measuring the intensity of spectral lines directly on the spectrum itself
    • G01J3/36Investigating two or more bands of a spectrum by separate detectors
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • G06T7/55Depth or shape recovery from multiple images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/10Cameras or camera modules comprising electronic image sensors; Control thereof for generating image signals from different wavelengths
    • H04N23/13Cameras or camera modules comprising electronic image sensors; Control thereof for generating image signals from different wavelengths with multiple sensors
    • H04N23/15Image signal generation with circuitry for avoiding or correcting image misregistration
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/50Constructional details
    • H04N23/55Optical parts specially adapted for electronic image sensors; Mounting thereof
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/56Cameras or camera modules comprising electronic image sensors; Control thereof provided with illuminating means
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/70Circuitry for compensating brightness variation in the scene
    • H04N23/743Bracketing, i.e. taking a series of images with varying exposure conditions
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N25/00Circuitry of solid-state image sensors [SSIS]; Control thereof
    • H04N25/10Circuitry of solid-state image sensors [SSIS]; Control thereof for transforming different wavelengths into image signals
    • H04N25/11Arrangement of colour filter arrays [CFA]; Filter mosaics
    • H04N25/13Arrangement of colour filter arrays [CFA]; Filter mosaics characterised by the spectral characteristics of the filter elements
    • H04N25/134Arrangement of colour filter arrays [CFA]; Filter mosaics characterised by the spectral characteristics of the filter elements based on three different wavelength filter elements
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N25/00Circuitry of solid-state image sensors [SSIS]; Control thereof
    • H04N25/10Circuitry of solid-state image sensors [SSIS]; Control thereof for transforming different wavelengths into image signals
    • H04N25/11Arrangement of colour filter arrays [CFA]; Filter mosaics
    • H04N25/13Arrangement of colour filter arrays [CFA]; Filter mosaics characterised by the spectral characteristics of the filter elements
    • H04N25/135Arrangement of colour filter arrays [CFA]; Filter mosaics characterised by the spectral characteristics of the filter elements based on four or more different wavelength filter elements
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N25/00Circuitry of solid-state image sensors [SSIS]; Control thereof
    • H04N25/40Extracting pixel data from image sensors by controlling scanning circuits, e.g. by modifying the number of pixels sampled or to be sampled
    • H04N25/41Extracting pixel data from a plurality of image sensors simultaneously picking up an image, e.g. for increasing the field of view by combining the outputs of a plurality of sensors
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N5/00Details of television systems
    • H04N5/222Studio circuitry; Studio devices; Studio equipment
    • H04N5/262Studio circuits, e.g. for mixing, switching-over, change of character of image, other special effects ; Cameras specially adapted for the electronic generation of special effects
    • H04N5/265Mixing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J3/2823Imaging spectrometer
    • G01J2003/2826Multispectral imaging, e.g. filter imaging
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10032Satellite or aerial image; Remote sensing
    • G06T2207/10036Multispectral image; Hyperspectral image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10141Special mode during image acquisition
    • G06T2207/10144Varying exposure
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20221Image fusion; Image merging

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Geometry (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Spectrometry And Color Measurement (AREA)
  • Blocking Light For Cameras (AREA)
  • Image Processing (AREA)
  • Image Input (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

本明細書で概説するように、本願の1つ以上の態様は、各アパーチャの上方に配置された湾曲したマルチバンドパスフィルタを備えたマルチアパーチャシステムを使用したスペクトルイメージング用システムおよびスペクトルイメージング技術に関する。また、本開示は、そのようなイメージングシステムから受信した画像情報を使用して、スペクトルデータキューブを生成するために、スペクトルアンミキシングおよび画像の位置合わせを実装するための技術に関する。さらに、本開示の態様は、そのようなデータキューブを使用して、画像化された物体の分析、たとえば、臨床での組織分析、生体認証または材料分析を行うことに関する。

Description

関連出願の相互参照
本願は、「高精度マルチアパーチャスペクトルイメージングのためのシステムおよび方法」という名称で2018年12月14日に出願された米国仮出願第62/780,121号および「高精度マルチアパーチャスペクトルイメージングのためのシステムおよび方法」という名称で2019年3月14日に出願された米国仮出願第62/818,375号の利益を主張するものであり、これらの文献は引用によりその全体があらゆる目的で本明細書に明示的に援用される。
連邦政府の助成による研究開発に関する陳述
本開示に記載の発明の一部は、契約番号HHSO100201300022Cの下、米国保健福祉省の事前準備・対応担当次官補局内の米生物医学先端研究開発局(BARDA)により付与された米国政府の支援を受けてなされたものである。また、本開示に記載の発明の一部は、契約番号W81XWH-17-C-0170および/または契約番号W81XWH-18-C-0114の下、米国国防保健局(DHA)により付与された米国政府の支援を受けてなされたものである。米国政府は、本発明に関し一定の権利を保有する場合がある。
本明細書で開示するシステムおよび方法は、スペクトルイメージングに関し、より具体的には、マルチアパーチャイメージングシステムによる高精度スペクトルイメージングに関する。
電磁スペクトルとは、電磁放射線(たとえば光)の波長または周波数の帯域である。電磁スペクトルは、長い波長から短い波長の順に、電波、マイクロ波、赤外線光(IR)、可視光(すなわちヒトの眼の構造により検出可能な光)、紫外線光(UV)、X線、およびガンマ線を含む。また、スペクトルイメージングは、像平面の様々な位置においてスペクトル情報の一部または完全なスペクトルを収集する分光学および写真撮影術の一分野を指す。マルチスペクトルイメージングシステムは、複数のスペクトル帯域(通常、それぞれ異なるスペクトル領域にある十数個以下のスペクトル帯域)を撮影して、画素ごとにそのスペクトル帯域での測定値を収集することができ、1個のスペクトルチャネルあたり数十nmの帯域幅を参照することができる。一方、ハイパースペクトルイメージングシステムは、さらに多くのスペクトル帯域を測定することができ、たとえば200個以上ものスペクトル帯域の測定が可能であり、いくつかの種類のハイパースペクトルイメージングシステムは、電磁スペクトルの一部に沿って狭帯域(たとえば数nm以下のスペクトル帯域幅)を連続的にサンプリングすることができる。
本明細書で開示するマルチスペクトルイメージングシステムおよびマルチスペクトルイメージング技術は、いくつかの特徴を有するが、このうちのどれが欠けても所望の特性を達成することはできない。本明細書で開示するスペクトルイメージングの特定の特徴を簡潔に説明するが、以下の説明は、後述の請求項で示す本発明の範囲を限定するものではない。当業者であれば、本明細書で開示するスペクトルイメージングの特徴が、従来のシステムや方法と比べてどのような利点を有するのかを十分に理解できるであろう。
一態様は、マルチスペクトル画像システムであって、
少なくとも1つのイメージセンサ;
前記少なくとも1つのイメージセンサの第1のセンサ領域へと入射光を通過させるように配置された第1のアパーチャ;
前記少なくとも1つのイメージセンサの第2のセンサ領域へと入射光を通過させるように配置された第2のアパーチャ;
第1のアパーチャの上方に配置され、少なくとも共通周波数帯の光および第1の特有周波数帯の光を通過させるように構成された第1のマルチバンドパスフィルタ;
第2のアパーチャの上方に配置され、少なくとも前記共通周波数帯の光および第2の特有周波数帯の光を通過させるように構成された第2のマルチバンドパスフィルタ;
マルチスペクトルデータキューブを生成させる命令を格納するメモリ;ならびに
前記命令に従って、少なくとも、第1のセンサ領域および第2のセンサ領域からシグナルを受信し、受信したこのシグナルを処理して、前記マルチスペクトルデータキューブを生成するように構成された少なくとも1つのプロセッサ
を含み、
前記マルチスペクトルデータキューブが、前記共通周波数帯に対応する第1のスペクトルチャネル、第1の特有周波数帯に対応する第2のスペクトルチャネル、および第2の特有周波数帯に対応する第3のスペクトルチャネルを含み;
前記シグナルの処理が、
少なくとも、前記共通周波数帯に対応する第1のスペクトルチャネルからの情報を利用して、第1のセンサ領域により撮影された第1の画像データと第2のセンサ領域により撮影された第2の画像データの間の視差を概算すること、ならびに
前記視差を利用して、第1のスペクトルチャネルの画像データ、第2のスペクトルチャネルの画像データおよび第3のスペクトルチャネルの画像データの位置調整を行い、マルチスペクトルデータキューブを生成することを含む、
マルチスペクトル画像システムに関する。
いくつかの実施形態において、第1のマルチバンドパスフィルタは、前記マルチスペクトル画像システムのレンズの前面に配置されるか、前記マルチスペクトル画像システムのレンズの後ろに配置されるか、前記マルチスペクトル画像システムのレンズの内部に配置されるか、または前記マルチスペクトル画像システムのレンズの素子の表面上に直接配置される。いくつかの実施形態において、前記少なくとも1つのプロセッサは、前記命令に従って、少なくともスペクトルアンミキシングを行うことによって、第1のスペクトルチャネルの画像データ、第2のスペクトルチャネルの画像データおよび第3のスペクトルチャネルの画像データを生成するように構成されている。いくつかの実施形態において、前記マルチスペクトル画像システムは、第1のセンサ領域および第2のセンサ領域を覆うカラーフィルタアレイをさらに含み、該カラーフィルタアレイが、カラーフィルタの繰り返しパターンを含み;前記繰り返しパターンの第1のカラーフィルタが、前記共通周波数帯を含む第1の重み付け光波長分布を通過させるように構成されており;前記繰り返しパターンの第2のカラーフィルタが、第1の特有周波数帯を含む第2の重み付け光波長分布を通過させるように構成されており;前記繰り返しパターンの第3のカラーフィルタが、第2の特有周波数帯を含む第3の重み付け光波長分布を通過させるように構成されている。いくつかの実施形態において、前記共通周波数帯は、紫色光、青紫色光、青色光および青緑色光の少なくとも1つを含む。いくつかの実施形態において、前記共通周波数帯は、380nmよりも長く450nmよりも短い波長を有する光を含む。いくつかの実施形態において、前記共通周波数帯は緑色光を含む。いくつかの実施形態において、前記共通周波数帯は赤色光を含む。いくつかの実施形態において、第1のマルチバンドパスフィルタおよび第2のマルチバンドパスフィルタのそれぞれが、緑色光、青緑色光および青色光の少なくとも1つを含む第2の共通周波数帯を通過させるように構成されている。いくつかの実施形態において、第1のマルチバンドパスフィルタおよび第2のマルチバンドパスフィルタのそれぞれが、赤色光を含む第2の共通周波数帯を通過させるように構成されている。いくつかの実施形態において、第1のマルチバンドパスフィルタおよび第2のマルチバンドパスフィルタのそれぞれが、前記カラーフィルタアレイのパターンを構成するカラーフィルタのうち、第1の共通周波数帯を遮断しない少なくとも1枚のカラーフィルタにより遮断されない第2の共通周波数帯を通過させるように構成されている。いくつかの実施形態において、前記少なくとも1つのプロセッサは、前記命令に従って、少なくともスペクトルアンミキシングを行うことによって、撮影画像データの、前記共通周波数帯に対応する部分と第2の共通周波数帯に対応する部分を分離するように構成されている。いくつかの実施形態において、前記少なくとも1つのプロセッサは、前記命令に従って、前記マルチスペクトルデータキューブを使用して画像化組織を少なくとも分析するように構成されている。
いくつかの実施形態において、前記マルチスペクトル画像システムは、
前記少なくとも1つのイメージセンサの第3のセンサ領域へと入射光を通過させるように配置された第3のアパーチャ;
前記少なくとも1つのイメージセンサの第4のセンサ領域へと入射光を通過させるように配置された第4のアパーチャ;
第3のアパーチャの上方に配置され、少なくとも共通周波数帯の光および第3の特有周波数帯の光を通過させるように構成された第3のマルチバンドパスフィルタ;ならびに
第4のアパーチャの上方に配置され、少なくとも前記共通周波数帯の光および第4の特有周波数帯の光を通過させるように構成された第4のマルチバンドパスフィルタ
をさらに含む。
いくつかの実施形態において、前記マルチスペクトル画像システムは、
前記少なくとも1つのイメージセンサの第5のセンサ領域へと入射光を通過させるように配置された第5のアパーチャ;および
第3のアパーチャの上方に配置された第5のマルチバンドパスフィルタ
をさらに含み、
該第5のマルチバンドパスフィルタが、第1のマルチバンドパスフィルタ、第2のマルチバンドパスフィルタ、第3のマルチバンドパスフィルタおよび第4のマルチバンドパスフィルタの少なくとも1つを通過する少なくとも1つの周波数帯の光を通過させるように構成されている。
いくつかの実施形態において、前記マルチスペクトル画像システムは、
前記少なくとも1つのイメージセンサの最大総数25個のセンサ領域へと入射光を通過させるように配置された最大総数25個のアパーチャ;および
前記最大総数25個のアパーチャのそれぞれの上方に配置された最大総数25枚のマルチバンドパスフィルタ
を含み、
各マルチバンドパスフィルタが、少なくとも共通周波数帯の光と、前記画像システム内の他のマルチバンドパスフィルタを通過しない特有周波数帯の光とを通過させるように構成されており;
通過する周波数帯の総数が、4個を超え、かつ51個未満である。
いくつかの実施形態において、前記マルチスペクトル画像システムは、
前記少なくとも1つのイメージセンサの第3のセンサ領域へと入射光を通過させるように配置された第3のアパーチャ;および
第3のアパーチャの上方に配置された第3のマルチバンドパスフィルタ
をさらに含み、
該第3のマルチバンドパスフィルタが、第1のマルチバンドパスフィルタおよび第2のマルチバンドパスフィルタの少なくとも一方を通過する少なくとも1つの周波数帯の光を通過させるように構成されている。
いくつかの実施形態において、第1のマルチバンドパスフィルタ、第2のマルチバンドパスフィルタ、第3のマルチバンドパスフィルタおよび第4のマルチバンドパスフィルタは、これら全体で7個、8個、9個または10個の異なる周波数帯を通過させ、これらの異なる周波数帯に、共通周波数帯、第1の特有周波数帯、第2の特有周波数帯、第3の特有周波数帯および第4の特有周波数帯が含まれている。いくつかの実施形態において、第1のマルチバンドパスフィルタ、第2のマルチバンドパスフィルタ、第3のマルチバンドパスフィルタおよび第4のマルチバンドパスフィルタのそれぞれは、青色光を含む共通周波数帯、緑色光を含む追加の共通周波数帯、および2個の追加の周波数帯を通過させるように構成されている。いくつかの実施形態において、前記少なくとも1つのプロセッサは、前記命令に従って、前記7個、8個、9個または10個の異なる周波数帯に対応する10個のスペクトルチャネルを含むマルチスペクトルデータキューブを少なくとも生成するように構成されている。
いくつかの実施形態において、前記マルチスペクトル画像システムは、前記7個、8個、9個または10個の異なる周波数帯の光を発するように構成された複数の発光ダイオード(LED)を含む照明基板をさらに含み、個々のLEDが、前記7個、8個、9個または10個の異なる周波数帯のいずれかの光を発する。いくつかの実施形態において、前記少なくとも1つのイメージセンサは、該少なくとも1つのイメージセンサのフォトダイオードの上方に配置されたカラーフィルタアレイを含み、該カラーフィルタアレイが、少なくとも2枚の異なるカラーフィルタを有する。いくつかの実施形態において、前記カラーフィルタアレイは、3~25枚の異なるカラーフィルタを有する。いくつかの実施形態において、前記7個、8個、9個または10個の異なる周波数帯は、前記3枚の異なるカラーフィルタを通過する周波数帯よりも狭い周波数帯の光を含む。いくつかの実施形態において、前記マルチスペクトル画像システムは、前記共通周波数帯の光、第1の特有周波数帯の光および第2の特有周波数帯の光を発するように構成された少なくとも1つの光源をさらに含む。いくつかの実施形態において、前記マルチスペクトル画像システムは、前記少なくとも1つの光源の上方に配置された拡散光学素子をさらに含む。いくつかの実施形態において、少なくとも第1のマルチバンドパスフィルタは、湾曲したフィルタである。
別の一態様は、本明細書に記載のマルチスペクトル画像システムを使用してマルチスペクトルデータキューブを生成する方法であって、
第1のスペクトルチャネル、第2のスペクトルチャネルおよび第3のスペクトルチャネルを含む第1の露光を撮影する工程;
第1のスペクトルチャネル、第2のスペクトルチャネルおよび第3のスペクトルチャネルを含む第2の露光を撮影する工程;ならびに
第1の露光のシグナル/ノイズ比を低下させるため、第2の露光の使用に少なくとも部分的に基づいてマルチスペクトルデータキューブを生成する工程
を含む方法に関する。
いくつかの実施形態において、前記方法は、第1の露光中および第2の露光中に、各光源が全体として、共通周波数帯の光、第1の特有周波数帯の光および第2の特有周波数帯の光を発する光源群を起動する工程をさらに含む。いくつかの実施形態において、第1の露光および第2の露光のそれぞれは少なくとも第4のスペクトルチャネルを含み、第1の露光および第2の露光のそれぞれに含まれるスペクトルチャネルの総数は10個以下である。
いくつかの実施形態において、前記方法は、
環境光を使用して第3の露光を撮影する工程、ならびに
第1の露光および第2の露光から第3の露光を差し引く工程
をさらに含む。
いくつかの実施形態において、前記方法は、第1の露光および第2の露光から第3の露光を差し引いた結果に基づいて、前記マルチスペクトルデータキューブを生成する工程をさらに含む。いくつかの実施形態において、前記方法は、第1のスペクトルチャネルの情報に基づいて、第1のセンサ領域により撮影された画像と第2のセンサ領域により撮影された画像の位置合わせを行う工程をさらに含む。いくつかの実施形態において、前記方法は、第1のセンサ領域により撮影された画像と第2のセンサ領域により撮影された画像の位置合わせの結果に基づいて、前記マルチスペクトルデータキューブを生成する工程をさらに含む。
いくつかの実施形態において、前記画像のそれぞれは、第1のスペクトルチャネルにおいて第1のセンサ領域により撮影された情報および第1のスペクトルチャネルにおいて第2のセンサ領域により撮影された情報を示し、
前記方法は、
前記共通周波数帯に対応する第1のスペクトルチャネルにおいて、第1のセンサ領域により撮影された画像と第2のセンサ領域により撮影された画像の間の視差を特定する工程;
前記視差を使用して、第2のスペクトルチャネルにおいて第1のセンサ領域により撮影された情報を示す第1の追加の画像を前記複数の画像に位置合わせする工程;ならびに
前記視差を使用して、第3のスペクトルチャネルにおいて第2のセンサ領域により撮影された情報を示す第2の追加の画像を前記複数の画像に位置合わせする工程
をさらに含む。
いくつかの実施形態において、各露光は、4個以上のスペクトルチャネルを含み;
前記方法は、前記視差を使用して、第4のスペクトルチャネルにおいて第3のセンサ領域により撮影された情報を示す少なくとも第3の追加の画像を前記複数の画像に位置合わせする工程をさらに含み、
前記複数の画像に位置合わせされる追加の画像の総数が、スペクトルチャネルの総数よりも1つ少ない。
いくつかの実施形態において、前記マルチスペクトル画像システムは、1つ以上の点、縞模様、グリッド、ランダムな斑点またはその他の空間パターンを投影するように構成された光源をさらに含み、
前記方法は、
前記光源によって、1つ以上の点、縞模様、グリッド、ランダムな斑点またはその他の空間パターンを含む光を、前記少なくとも1つのイメージセンサの視野内に投影する工程;および
前記光源により投影された光を利用して、追加の露光を撮影する工程
をさらに含み、
前記位置合わせが、この追加の露光に少なくとも部分的に基づく。
いくつかの実施形態において、前記光源により投影された前記光は、前記共通周波数帯の波長を有する光、広帯域照明の波長を有する光および累積的に可視可能な照明の波長を有する光のうちの少なくとも1つを含む。いくつかの実施形態において、前記追加の露光は、前記画像の前記位置合わせの正確度の確認に使用される。いくつかの実施形態において、前記追加の露光は、前記位置合わせの正確度を向上させるために前記視差の計算に含まれる。いくつかの実施形態において、前記位置合わせの正確度をさらに向上させるため、前記光源は、前記共通周波数帯、広帯域照明または累積的に可視可能な照明からの、特有の構成の、複数の点、縞模様、グリッド、ランダムな斑点またはその他の空間パターンを投影する。いくつかの実施形態において、前記マルチスペクトル画像システムは、前記少なくとも1つのイメージセンサの視野内の物体の形状を検出するように構成された追加のセンサをさらに含み、前記位置合わせは、この追加のセンサにより検出された形状に少なくとも部分的に基づく。この追加センサは、ライトフィールド情報を検知可能な1個のアパーチャセンサまたはマルチアパーチャセンサであってもよく、超音波やパルスレーザなどのその他のシグナルを検知可能な1個のアパーチャまたはマルチアパーチャセンサであってもよい。いくつかの実施形態において、前記方法は、第1のスペクトルチャネル、第2のスペクトルチャネルおよび第3のスペクトルチャネルを分離するために、スペクトルアンミキシングを行う工程をさらに含む。いくつかの実施形態において、各露光は、4個以上のスペクトルチャネルを含み;前記方法は、これらのスペクトルチャネルのすべてを分離するために、スペクトルアンミキシングを行う工程をさらに含む。いくつかの実施形態において、前記方法は、前記視差を使用して、トポグラフィー情報および深さ情報の少なくとも一方を測定する工程をさらに含む。いくつかの実施形態において、前記視差を使用して、個々の画素の深さ情報を測定する。いくつかの実施形態において、前記方法は、前記深さ情報を使用して、画像化した創傷の体積を測定する工程をさらに含む。いくつかの実施形態において、前記方法は、非可視光を示す第4のスペクトルチャネルを少なくとも含む追加の露光を撮影する工程;および前記追加の露光に少なくとも部分的に基づいて前記マルチスペクトルデータキューブを生成する工程をさらに含む。いくつかの実施形態において、前記方法は、前記追加の露光中に前記非可視光を発する第3の光源群を起動する工程をさらに含む。
別の一態様は、自動組織分析方法であって、
本明細書に記載のマルチスペクトル画像システムのいずれかを使用して、画像化された組織を示すマルチスペクトルデータキューブを生成する工程;
組織の少なくとも1つの臨床的特徴を予測するように訓練された機械学習システムへの入力として前記マルチスペクトルデータキューブを提供する工程;および
前記機械学習システムの出力に基づいて、前記画像化された組織の少なくとも1つの臨床的特徴を予測する工程
を含む方法に関する。
いくつかの実施形態において、前記方法は、本明細書に記載のマルチスペクトルデータキューブ生成方法を使用して、前記マルチスペクトルデータキューブを生成する工程をさらに含む。いくつかの実施形態において、前記少なくとも1つの臨床的特徴は、複数の創傷状態または熱傷状態を含み、前記方法は、前記機械学習システムの前記出力を使用して、前記画像化された組織を、前記複数の創傷状態または熱傷状態の様々な種類ごとに分類した領域として示した画像を生成する工程をさらに含む。いくつかの実施形態において、たとえば、前記複数の熱傷状態は、少なくともI度熱傷、II度熱傷、III度熱傷または正常皮膚という種類を含み、前記方法は、前記画像化された組織を示す複数の画素のそれぞれを、I度熱傷、II度熱傷、III度熱傷または正常皮膚という種類に分類する工程をさらに含む。いくつかの実施形態において、たとえば、前記複数の創傷状態は、少なくとも止血、炎症、増殖、リモデリングまたは正常皮膚という種類を含み、前記方法は、前記画像化された組織を示す複数の画素のそれぞれを、止血、炎症、増殖、リモデリングまたは正常皮膚という種類に分類する工程をさらに含む。いくつかの実施形態において、前記少なくとも1つの臨床的特徴は治癒可能性を含み、前記方法は、前記機械学習システムの前記出力を使用して、前記画像化された組織の様々な領域において判定された治癒可能性を視覚的に示した、前記画像化された組織の画像を生成する工程をさらに含む。いくつかの実施形態において、前記少なくとも1つの臨床的特徴は、前記画像化された組織のがん状態または非がん状態を含み、前記方法は、前記画像化された組織を示す複数の画素のそれぞれを、がん状態または非がん状態として分類する工程をさらに含む。いくつかの実施形態において、前記少なくとも1つの臨床的特徴は、創傷の深さを含み、前記方法は、前記画像化された組織における創傷の深さを特定する情報を出力する工程をさらに含む。いくつかの実施形態において、前記少なくとも1つの臨床的特徴は、デブリードマンを行う際のマージンを含み、前記方法は、前記画像化された組織における創傷またはその周囲のデブリードマンを行う際のマージンを特定する情報を出力する工程をさらに含む。いくつかの実施形態において、前記少なくとも1つの臨床的特徴は、糖尿病性潰瘍、非糖尿病性潰瘍または慢性潰瘍を含む。いくつかの実施形態において、前記少なくとも1つの臨床的特徴は、糖尿病性足部潰瘍を含む。
別の一態様は、自動生体認証方法であって、
本明細書に記載のマルチスペクトル画像システムのいずれかを使用して、ユーザの生体的特徴を示すマルチスペクトルデータキューブを生成する工程;
特定の生体的特徴を認証するように訓練された機械学習システムへの入力として前記マルチスペクトルデータキューブを提供する工程;および
前記機械学習システムの出力に基づいて、前記ユーザを認証または拒否する工程
を含む方法に関する。
いくつかの実施形態において、前記方法は、本明細書に記載のマルチスペクトルデータキューブ生成方法のいずれかを使用して、前記マルチスペクトルデータキューブを生成する工程をさらに含む。
別の一態様は、自動材料分析方法であって、
本明細書に記載のマルチスペクトル画像システムのいずれかを使用して、画像化された材料を示すマルチスペクトルデータキューブを生成する工程;
1つ以上の材料特性を測定するように訓練された機械学習システムへの入力として前記マルチスペクトルデータキューブを提供する工程;および
前記機械学習システムの出力に基づいて、前記画像化された材料の1つ以上の材料特性を測定する工程
を含む方法に関する。
いくつかの実施形態において、前記方法は、本明細書に記載のマルチスペクトルデータキューブ生成方法のいずれかを使用して、前記マルチスペクトルデータキューブを生成する工程をさらに含む。
別の一態様は、
カラーフィルタアレイを有する少なくとも1つのイメージセンサ領域を含む携帯型のマルチアパーチャイメージング装置を使用して、該カラーフィルタアレイ中のカラーフィルタの色の数よりもチャネルの数が多い画像情報を撮影する工程;および
前記画像情報を処理して、各チャネルが0nmを超えかつ約40nm以下の範囲のそれぞれ異なる周波数帯の光を示す、前記チャネル数を含むマルチスペクトルデータキューブを作成する工程
を含む方法に関する。
いくつかの実施形態において、前記携帯型のマルチアパーチャイメージング装置は、本明細書に記載のマルチスペクトル画像システムのいずれかを含む。いくつかの実施形態において、前記方法は、本明細書に記載のマルチスペクトルデータキューブ生成方法のいずれかを使用して、前記マルチスペクトルデータキューブを生成する工程をさらに含む。
別の一態様は、
本明細書に記載のマルチスペクトル画像システムのいずれかを使用して、マルチスペクトルデータキューブを生成する工程;および
前記マルチスペクトル画像システムの少なくとも1つのプロセッサを使用して、前記マルチスペクトルデータキューブを処理することによって、該マルチスペクトルデータキューブにより示される物体の少なくとも1つの特性を測定する工程
を含む方法に関する。
いくつかの実施形態において、前記方法は、本明細書に記載のマルチスペクトルデータキューブ生成方法のいずれかを使用して、前記マルチスペクトルデータキューブを生成する工程をさらに含む。
別の一態様は、
本明細書に記載のマルチスペクトル画像システムのいずれかを使用して、マルチスペクトルデータキューブを生成する工程;
ネットワーク上の少なくとも1つの遠隔プロセッサに前記マルチスペクトルデータキューブを送信する工程;および
前記少なくとも1つの遠隔プロセッサを使用して、前記マルチスペクトルデータキューブを処理することによって、該マルチスペクトルデータキューブにより示される物体の少なくとも1つの特性を測定する工程
を含む方法に関する。
いくつかの実施形態において、前記方法は、本明細書に記載のマルチスペクトルデータキューブ生成方法のいずれかを使用して、前記マルチスペクトルデータキューブを生成する工程をさらに含む。
様々な主光線入射角でフィルタに入射する光の一例を示す。
様々な主光線入射角の光が透過した図1Aのフィルタから得られた透過効率の一例を示すグラフである。
マルチスペクトル画像データキューブの一例を示す。
特定のマルチスペクトルイメージング技術が、どのようにして図2Aのデータキューブを生成するのかを示した一例を示す。
図2Aのデータキューブを生成することができるスナップショットイメージングシステムの一例を示す。
本開示による、湾曲したマルチバンドパスフィルタを備えたマルチアパーチャイメージングシステムの一例の光学設計を示した断面概略図を示す。
図3Aに示したマルチアパーチャイメージングシステムの1つの光路を構成する光学部品の光学設計の一例を示す。
図3Aおよび図3Bに関連して述べた光学設計を備えたマルチスペクトル・マルチアパーチャイメージングシステムの一実施形態を示す。
図3Aおよび図3Bに関連して説明した光学設計を備えたマルチスペクトル・マルチアパーチャイメージングシステムの別の一実施形態を示す。
図3Aおよび図3Bに関連して説明した光学設計を備えたマルチスペクトル・マルチアパーチャイメージングシステムの別の一実施形態を示す。
図3Aおよび図3Bに関連して説明した光学設計を備えたマルチスペクトル・マルチアパーチャイメージングシステムの別の一実施形態を示す。
図3Aおよび図3Bに関連して説明した光学設計を備えたマルチスペクトル・マルチアパーチャイメージングシステムの別の一実施形態を示す。
図3Aおよび図3Bに関連して説明した光学設計を備えたマルチスペクトル・マルチアパーチャイメージングシステムの別の一実施形態を示す。
図3Aおよび図3Bに関連して説明した光学設計を備えたマルチスペクトル・マルチアパーチャイメージングシステムの別の一実施形態を示す。
図3A~10Bに示したマルチスペクトル・マルチアパーチャイメージングシステムのフィルタを通過することができる1組の周波数帯の一例を示す。
図3A~10Bに示したマルチスペクトル・マルチアパーチャイメージングシステムに使用することができるイメージングシステムの略ブロック図を示す。
図3A~10Bに示したマルチスペクトル・マルチアパーチャイメージングシステムを使用して、画像データを撮影するプロセスの一例を示したフローチャートである。
たとえば、図3A~10Bに示したマルチスペクトル・マルチアパーチャイメージングシステムを使用し、かつ/または図13に示したプロセスを使用して撮影された画像データなどの画像データを処理するためのワークフローの略ブロック図を示す。
たとえば、図3A~10Bに示したマルチスペクトル・マルチアパーチャイメージングシステムを使用し、かつ/または図13に示したプロセスを使用して撮影された画像データなどの画像データを処理するための視差および視差補正を図示する。
たとえば、図3A~10Bに示したマルチスペクトル・マルチアパーチャイメージングシステムを使用して、かつ/または図13に示したプロセスを使用して撮影した後、図14および図15に従って処理された画像データなどのマルチスペクトル画像データを画素単位で分類するためのワークフローを図示する。
図3A~10Bに示したマルチスペクトル・マルチアパーチャイメージングシステムを含む分散コンピューティングシステムの一例の略ブロック図を示す。
携帯型マルチスペクトル・マルチアパーチャイメージングシステムの実施形態の一例を示す。
携帯型マルチスペクトル・マルチアパーチャイメージングシステムの実施形態の一例を示す。
一般的なカメラのハウジングに収容された、小型USB3.0用のマルチスペクトル・マルチアパーチャイメージングシステムの一例を示す。
画像の位置合わせを向上させるための追加の光源を含むマルチスペクトル・マルチアパーチャイメージングシステムの一例を示す。
本明細書で概説するように、本開示は、各アパーチャの上方に配置された湾曲したマルチバンドパスフィルタを備えたマルチアパーチャシステムを使用したスペクトルイメージングに関する。また、本開示は、そのようなイメージングシステムから受信した画像情報を使用して、スペクトルデータキューブを生成するために、スペクトルアンミキシングおよび画像の位置合わせを実装するための技術に関する。本開示の技術は、後述するように、画像化される物体から反射された波長帯域に関する正確な情報を示す画像データを得る際に、スペクトルイメージングにおいて通常見られる様々な課題に対処するものである。いくつかの実施形態において、本明細書に記載のシステムおよび方法は、短時間で(たとえば6秒以内に)広範囲な組織(たとえば5.9×7.9インチ)の画像を取得することができ、イメージング造影剤の注射を必要とせずに、これを達成することができる。いくつかの態様において、たとえば、本明細書に記載のマルチスペクトル画像システムは、たとえば5.9×7.9インチといった広範囲な組織から6秒以内に画像を取得するように構成されており、かつ複数の熱傷状態、複数の創傷状態もしくは治癒可能性の識別;または画像化された組織のがん状態もしくは非がん状態、創傷の深さ、デブリードマンの際のマージン、糖尿病性潰瘍、非糖尿病性潰瘍もしくは慢性潰瘍の存在を含む臨床的特徴などの組織分析情報を、イメージング造影剤の非存在下で提供するように構成されている。同様に、本明細書に記載の方法のいくつかにおいて、本発明のマルチスペクトル画像システムは、たとえば5.9×7.9インチといった広範囲な組織から6秒以内に画像を取得し、かつ複数の熱傷状態、複数の創傷状態もしくは治癒可能性の識別;または画像化された組織のがん状態または非がん状態、創傷の深さ、デブリードマンの際のマージン、糖尿病性潰瘍、非糖尿病性潰瘍もしくは慢性潰瘍の存在を含む臨床的特徴などの組織分析情報を、イメージング造影剤の非存在下で出力することができる。
既存の解決策における課題の1つとして、撮影した画像の色が歪んでしまい、画像データの品質が損なわれることがある。これは、特に、光学フィルタを使用して特定の波長の光を正確に検出し分析することに依存した用途において問題となりうる。より具体的には、色の濃淡は、イメージセンサ領域を通過する位置に応じた光の波長の変化であり、これは、カラーフィルタへの光の入射角が大きくなるほど、カラーフィルタを透過する光の波長が短くなることによって起こる。通常、このような効果は、干渉系フィルタにおいて観察され、干渉系フィルタは、透明な基材に様々な屈折率の薄層を複数重ねて製造される。したがって、長い波長(赤色光など)であるほど、光線入射角が大きくなることから、イメージセンサの端部において阻止されやすくなり、イメージセンサを通過する際に、同じ波長の入射光が空間的には不均一な色として検出される。画像補正をしない場合、色の濃淡は、撮影された画像の端部付近の色の変化として現れる。
本開示の技術は、レンズおよび/またはイメージセンサの構成と、これらの視野またはアパーチャサイズが限定されないため、上市されているその他のマルチスペクトルイメージングシステムと比べて有利である。当業者であれば熟知しているように、本明細書で開示するイメージングシステムのレンズ、イメージセンサ、アパーチャサイズ、またはその他の部品の変更は、該イメージングシステムへのその他の調整を含んでいてもよいことを理解できるであろう。また、本開示の技術は、波長を分割する機能を有する部品、またはイメージングシステム全体を利用して波長を分割する機能を有する部品(たとえば光学フィルタなど)を、光エネルギーをデジタル出力に変換する部品(たとえばイメージセンサなど)から分離可能であるという点で、その他のマルチスペクトルイメージングシステムと比べて改良されている。これによって、コストを削減し、簡略化し、かつ/または異なるマルチスペクトル波長用のイメージングシステムを再構成するための開発時間を削減することができる。また、本開示の技術は、より小さくより軽量のフォームファクタにおいて、上市されているその他のマルチスペクトルイメージングシステムと同等の撮像特性を達成可能であるという点で、その他のマルチスペクトルイメージングシステムよりも堅牢であると考えられる。また、本開示の技術は、スナップショット、ビデオレートまたは高速ビデオレートでマルチスペクトル画像を取得できるという点でも、その他のマルチスペクトルイメージングシステムと比べて有利である。また、本開示の技術は、いくつかのスペクトル帯域を各アパーチャに多重化することができるため、イメージングデータセット中の特定の数のスペクトル帯域を取得するのに必要とされるアパーチャの数を減らすことができ、アパーチャの数の削減と改良された光収集によりコストが削減できることから(たとえば、市販のセンサアレイでは、固定されたサイズおよび寸法に収められる程度のさらに大きなアパーチャが使用されうる)、マルチアパーチャ技術に基づいたマルチスペクトルイメージングシステムのさらに堅牢な実装を提供することができる。さらに、本開示の技術は、解像度または画質に関して妥協することなく、これらの利点のすべてを提供することができる。
図1Aは、イメージセンサ110に向かう光路に沿って配置されたフィルタ108の一例と、様々な光線入射角でフィルタ108に入射する光を示す。光線102A、光線104A、光線106Aは、フィルタ108を通過した後、レンズ112により屈折されてセンサ110上に入射する線として示す。レンズ112は、ミラーおよび/またはアパーチャなどの(ただしこれらに限定されない)別の画像形成光学系と置き換えてもよい。図1Aにおいて、各光線の光は広帯域であると見なされ、たとえば、広い波長領域のスペクトル組成を有する光であり、フィルタ108により特定波長の光のみが選択的に透過される。3つの光線102A、104A、106Aは、互いに異なる入射角でフィルタ108に入射する。説明を目的として、光線102Aは、フィルタ108に対して実質的に垂直な入射として示され、光線104Aは、光線102Aよりも大きな入射角を有し、光線106Aは、光線104Aよりも大きい入射角を有する。フィルタを通過した各光線102B、104B、106Bは、入射角に依存したフィルタ108の透過特性により、特有のスペクトラムを示し、これらの特有のスペクトラムは、センサ110で検出される。この入射角度依存性という効果から、入射角が増加するほど波長が短くなるという、フィルタ108の帯域通過特性のシフトが起こる。さらに、入射角度依存性によって、フィルタ108の透過効率が低下することがあり、フィルタ108の帯域通過特性を示すスペクトル形状が変化することがある。このような複合効果は、入射角度依存性のスペクトル透過と呼ばれる。図1Bは、センサ110の位置に存在すると仮定した分光器で検出された図1Aの各光線のスペクトラムを示し、入射角が増加するほど、フィルタ108の帯域通過特性を示すスペクトルがシフトすることを説明している。曲線102C、曲線104Cおよび曲線106Cは、帯域通過特性の中心波長が短くなっていることを示し、したがって、この例では、光学系から照射された光の波長が短縮していることが示されている。また、この図に示すように、帯域通過特性のスペクトル形状および透過率のピークも入射角により変化する。特定のコンシューマー向けの用途では、このような入射角度依存性のスペクトル透過による目に見える効果を除去するため、画像処理を行うことができる。しかし、このような後加工技術では、どの波長の光が実際にフィルタ108に入射したのかという正確な情報を回復することはできない。したがって、得られる画像データは、特定の高精密な用途には使用できない可能性がある。
既存の特定のスペクトルイメージングシステムが直面している別の課題として、完全なスペクトル画像データセットの撮影に必要とされる時間が挙げられる。これについては、図2Aおよび図2Bに関連して検討する。スペクトルイメージセンサは、特定のシーンのスペクトル放射照度I(x,y,λ)を取得し、データキューブと一般に呼ばれる三次元(3D)データセットを収集する。図2Aは、スペクトル画像のデータキューブ120の一例を示す。この図に示すように、データキューブ120は、三次元の画像データを示し、このうち、2つの空間次元(xおよびy)は、イメージセンサの二次元(2D)表面に対応し、スペクトル次元(λ)は、特定の波長域に対応している。データキューブ120のサイズは、NxNyNλで求めることができ、NxおよびNyは、各空間次元(x,y)方向の標本点の数であり、Nλは、スペクトル軸λ方向の標本点の数である。データキューブは、現在利用可能な2D検出器アレイ(たとえばイメージセンサ)よりも次元が高いため、典型的なスペクトルイメージングシステムでは、データキューブ120の時系列2Dスライス(すなわち平面)(本明細書において「走査型」イメージングシステムと呼ぶ)を撮影するか、あるいは、演算過程でデータキューブ120を分割して、このデータキューブ120に再統合可能な複数の二次元データを得ることにより、データキューブ120のすべての標本点を同時に測定する(本明細書において「スナップショット型」イメージングシステムと呼ぶ)。
図2Bは、特定の走査スペクトルイメージング技術が、どのようにしてデータキューブ120を生成するのかを示した一例を示す。具体的には、図2Bは、1回の検出器の積分時間に収集することができるデータキューブ120の部分132、134および136を示す。ポイントスキャン分光器は、たとえば、1つの空間位置(x,y)においてスペクトル平面λの全体に延びる部分132を撮影することができる。ポイントスキャン分光器は、空間次元の各(x,y)位置に対して積分を複数回行うことによってデータキューブ120を構築することができる。フィルタホイールイメージングシステムは、たとえば、空間次元xおよびyの全体に延びるが、単一のスペクトル平面λのみに位置する部分134を撮影することができる。フィルタホイールイメージングシステムなどの波長走査イメージングシステムは、スペクトル平面λに対して積分を複数回行うことによってデータキューブ120を構築することができる。ラインスキャン分光器は、たとえば、スペクトル次元λの全体と空間次元の一方(xまたはy)の全体に延びるが、もう一方の空間次元(yまたはx)の1つの点のみに位置する部分136を撮影することができる。ラインスキャン分光器は、この一方の空間次元(yまたはx)位置に対して積分を複数回行うことによってデータキューブ120を構築することができる。
対象物体とイメージングシステムの両方が静止している(または露光時間において比較的静止した状態の)用途では、前述した走査型イメージングシステムは、高解像度のデータキューブ120が得られるという利点がある。ラインスキャンイメージングシステムおよび波長走査型イメージングシステムは、イメージセンサの全域を使用して各スペクトルまたは各空間画像が撮影されることにより利点が得られる場合がある。しかし、最初の露光と次の露光の間にイメージングシステムおよび/または物体に動きがあると、得られる画像データにアーチファクトが発生する場合がある。たとえば、データキューブ120において同じ(x,y)位置であっても、実際には、スペクトル次元λにおいて、画像化した物体上の別の物理的位置が示される場合がある。このようなアーチファクトによって、後続の分析において誤差が生じることがあり、かつ/または位置合わせを行う必要性が出てくることがある(たとえば、特定の(x,y)位置が物体上の同じ物理的位置に対応するように、スペクトル次元λの位置調整を行う必要性が出てくる)。
これに対して、スナップショットイメージングシステム140では、1つの積分時間または1回の露光でデータキューブ120の全体を撮影することができることから、前述のような動きが画質に与える影響を回避することができる。図2Cは、スナップショットイメージングシステムの作製に使用可能なイメージセンサ142および光学フィルタアレイ144(カラーフィルタアレイ(CFA)など)の一例を示す。この一例でのカラーフィルタアレイ(CFA)144は、イメージセンサ142の表面にカラーフィルタユニット146の繰り返しパターンが配置されたものである。スペクトル情報を得るためのこの方法は、マルチスペクトルフィルタアレイ(MSFA)またはスペクトル分解検出器アレイ(SRDA)と呼ぶこともできる。図に示した一例では、カラーフィルタユニット146は、5×5に配置された様々なカラーフィルタを含むことから、得られる画像データにおいて25個のスペクトルチャネルが生成される。CFAは、様々なカラーフィルタを使用することによって入射光を各カラーフィルタの帯域に分割し、分割した光をイメージセンサ上の各色の感光体へと導くことができる。このようにして実際には、所定の色148に対して、25個の感光体のうちの1個のみが、その色の波長の光を示すシグナルを検出する。したがって、このようなスナップショットイメージングシステム140を使用した場合、1回の露光で25種のカラーチャネルが生成されうるが、各カラーチャネルの測定データ量は、イメージセンサ142の総出力よりも小さくなる。いくつかの実施形態において、CFAは、フィルタアレイ(MSFA)およびスペクトル分解検出器アレイ(SRDA)の一方または両方を含んでいてもよく、かつ/または従来のべイヤーフィルタ、CMYKフィルタもしくはその他の吸収系フィルタもしくは干渉系フィルタを含んでいてもよい。干渉系フィルタの一種として、薄膜フィルタアレイがグリッド状に配置され、各グリッドの素子が1つ以上のセンサ素子に対応しているものがある。別の種類の干渉系フィルタとして、ファブリ・ペローフィルタがある。ナノエッチングされたファブリ・ペロー干渉フィルタは、通常、半値全幅(FWHM)が約20~50nmの帯域通過特性を示すが、フィルタの通過帯域の中心波長から阻止帯域への移行部においてロールオフが小さいことから、いくつかの実施形態において使用することできるという点で有利である。さらに、このようなフィルタは、阻止帯域において低い光学濃度(OD)を示すことから、通過帯域外の光に対する感度を向上させることができる。このような複合効果を有することから、これらの特定のフィルタは、蒸着やイオンビームスパッタリングなどのコーティング蒸着法により多数の薄膜層で構成された同様のFWHMの高光学濃度干渉フィルタの高いロールオフにより阻止されてしまうスペクトル領域に対して感度を示す。色素系CMYKフィルタまたはRGB系(べイヤー)フィルタを使用する構成の実施形態では、スペクトルのロールオフが小さく、かつ個々のフィルタの通過帯域のFWHMが大きいことが好ましく、これにより、観察されるスペクトラム全体において各波長に対するスペクトル透過率を特有のものとすることができる。
したがって、スナップショットイメージングシステムにより得られるデータキューブ120は、高精度なイメージング用途において問題となりうる以下の2つの特性の一方を有する。1つ目の問題点として、スナップショットイメージングシステムにより得られるデータキューブ120は、検出器アレイの(x,y)の大きさよりもNxおよびNyの大きさが小さくなりうることから、同じイメージセンサを有する走査型イメージングシステムにより得られる解像度と同程度までデータキューブ120の解像度が低くなる。2つ目の問題点として、スナップショットイメージングシステムにより得られるデータキューブ120では、特定の(x,y)位置に対する数値補間により、NxおよびNyの大きさが検出器アレイの(x,y)の大きさと同じになることがある。しかし、データキューブの生成において補間が行われた場合、データキューブの特定の数値は、センサ上の入射光の波長の実測値ではなく、周辺値に基づいた実測の推定値であることを意味する。
露光を1回のみ行うマルチスペクトルイメージングにおいて使用される別の既存の部材として、マルチスペクトルイメージング用ビームスプリッターがある。このようなイメージングシステムにおいて、ビームスプリッターキューブは、入射光を異なる色帯域へと分割し、各帯域は独立した複数のイメージセンサにより観測される。ビームスプリッターの設計を変更して、測定するスペクトル帯域を調整することができるが、システムのパフォーマンスを損なうことなく、4つを超えるビームに入射光を分割することは容易ではない。したがって、4個のスペクトルチャネルが、このアプローチでの実用限界であると見られる。これと密接に関連した方法では、かさ高いビームスプリッターキューブ/プリズムの代わりに、薄膜フィルタを使用して光を分割するが、連続した複数のフィルタによる累積的な透過率の減衰と空間的制限により、このアプローチでも、約6個のスペクトルチャネルに制限される。
いくつかの実施形態において、とりわけ前述の課題は、各アパーチャを通って入射する光を選択的に透過する湾曲したマルチバンドパスフィルタを備えた、本明細書で開示するマルチアパーチャスペクトルイメージングシステムと、これに関連する画像データ処理技術により対処することができる。この特定の構成によって、速いイメージング速度、高解像度画像、および検出波長の高い忠実度という設計目標のすべてを達成することができる。したがって、本明細書で開示する光学設計および関連する画像データ処理技術は、ポータブルスペクトルイメージングシステムにおいて使用することができ、かつ/または移動している対象のイメージングに使用することができるとともに、高精度用途(たとえば、臨床組織分析、生体認証、一時的な臨床徴候)に適したデータキューブを得ることができる。このような高精度用途として、転移前の悪性黒色腫の進行度(0~3)の診断、皮膚組織における熱傷の重症度の分類、または糖尿病性足部潰瘍の重症度の組織診断を挙げることができる。したがって、いくつかの実施形態で説明するスナップショットスペクトルの取得とスモールフォームファクタによって、様々な種類の網膜症(たとえば、非増殖性糖尿病性網膜症、増殖性糖尿病性網膜症および加齢黄斑変性)の診断、ならびに動きが多い小児患者のイメージングを含む、一時的事象を扱う臨床環境において本発明を使用することができる。したがって、本明細書で開示する平坦または湾曲したマルチバンドパスフィルタを備えたマルチアパーチャシステムの使用は、従来のスペクトルイメージングの実装と比べて大幅な技術的進歩となることを当業者であれば十分に理解できるであろう。より具体的には、本発明のマルチアパーチャシステムは、各アパーチャ間における算出された遠近差に基づいて、物体の弯曲、深さ、体積および/または面積の3D空間画像またはこれらに関する3D空間画像の収集を可能にするものであってもよい。しかし、本明細書で開示するマルチアパーチャを用いた方法は、特定のフィルタには限定されず、干渉フィルタリングまたは吸収フィルタリングに基づいた平坦なフィルタおよび/または薄型フィルタを含んでいてもよい。本明細書で開示する本発明は、適切なレンズまたはアパーチャに対する入射角が小さいか、許容域である場合に、イメージングシステムの像空間に平坦なフィルタを備えるように変更することができる。結像レンズの入射瞳および/もしくは射出瞳またはアパーチャ絞りにフィルタを配置してもよく、光工学分野の当業者であれば、このような構成を容易に採用できるであろう。
特定の例および実施形態に関連して、本開示の様々な態様について以下で述べる。これらの例および実施形態は説明を目的としたものであり、本開示を限定するものではない。本明細書で述べる例および実施形態は、説明を目的として、特定の計算およびアルゴリズムに焦点を当てているが、当業者であれば、これらの例が説明のみを目的としたものであり、本発明を限定するものではないことを十分に理解できるであろう。たとえば、マルチスペクトルイメージングに関していくつかの例を挙げているが、本明細書で開示するマルチアパーチャイメージングシステムおよびこれに関連するフィルタは、その他の実装においてハイパースペクトルイメージングを実施できるように構成することができる。さらに、特定の例は、携帯型用途および/または移動している対象物体に対する用途において利点を達成するものとして提示されているが、本明細書で開示するイメージングシステムの設計およびこれに関連する処理技術により、固定されたイメージングシステムおよび/または比較的静止した状態の対象の分析に適した高精度なデータキューブを得ることができることは十分に理解できるであろう。
電磁スペクトル範囲およびイメージセンサの概要
本明細書において、電磁スペクトルの特定の色または特定の部分に言及し、以下、ISO21348「放射照度スペクトルの種類の定義」による定義に従って、これらの波長を説明する。以下で詳述するように、特定のイメージング用途において、特定の色の波長領域を一括して特定のフィルタに通過させることができる。
760nm~380nmまたは約760nm~約380nmの波長の範囲の電磁放射線は、通常、「可視」スペクトラム、すなわち、ヒトの眼の色受容体により認識可能なスペクトラム部分であると考えられている。可視スペクトルにおいて、赤色光は、通常、700nmもしくは約700nmの波長を有するか、または760nm~610nmもしくは約760nm~約610nmの範囲にあると考えられている。橙色光は、通常、600nmもしくは約600nmの波長を有するか、または610nm~591nmもしくは約610nm~約591nmの範囲にあると考えられている。黄色光は、通常、580nmもしくは約580nmの波長を有するか、または591nm~570nmもしくは約591nm~約570nmの範囲にあると考えられている。緑色光は、通常、550nmもしくは約550nmの波長を有するか、または570nm~500nmもしくは約570nm~約500nmの範囲にあると考えられている。青色光は、通常、475nmもしくは約475nmの波長を有するか、または500nm~450nmもしくは約500nm~約450nmの範囲にあると考えられている。紫色(赤紫色)光は、通常、400nmもしくは約400nmの波長を有するか、または450nm~360nmもしくは約450nm~約360nmの範囲にあると考えられている。
可視スペクトル外の範囲に目を向けると、赤外線(IR)は、可視光よりも長い波長を有する電磁放射線を指し、通常、ヒトの眼には見えない。IRは、約760nmまたは760nmにある赤色の可視スペクトルの公称下限値から約1mmまたは1mmに及ぶ波長を有する。この範囲内において、近赤外線(NIR)は、赤色の範囲に隣接しているスペクトラムの一部を指し、約760nm~約1400nmまたは760nm~1400nmの波長範囲である。
紫外線(UV)光は、可視光よりも短い波長の電磁放射線の一部を指し、通常、ヒトの眼には見えない。UVは、約40nmまたは40nmにある紫色の可視スペクトルの公称値下限値から約400mmに及ぶ波長を有する。この範囲内において、近紫外線(NUV)は、紫色の範囲に隣接しているスペクトラムの一部を指し、約400nm~約300nmまたは400nm~300nmの波長範囲であり、中間紫外線(MUV)は、約300nm~約200nmまたは300nm~200nmの波長範囲であり、遠紫外線(FUV)は、約200nm~約122nmまたは200nm~122nmの波長範囲である。
本明細書に記載のイメージセンサは、特定の用途に適した特定の波長領域に応じて、上述の領域のうち、どの電磁放射線であっても検出できるように構成することができる。通常のシリコン電荷結合素子(CCD)または相補型金属酸化膜半導体(CMOS)センサのスペクトル感度は、可視スペクトルの全体に及ぶが、近赤外線(IR)スペクトラムの大半およびUVスペクトラムの一部にも及ぶ。いくつかの実装では、裏面入射型または表面入射型のCCDアレイまたはCMOSアレイを付加的または代替的に使用することができる。高いSN比および科学グレードの測定が必要とされる用途では、いくつかの実装において、科学用途の相補型金属酸化膜半導体(sCMOS)カメラまたは電子増幅CCDカメラ(EMCCD)を付加的または代替的に使用することができる。別の実装では、目的とする用途に応じて、特定の色の範囲(たとえば、短波長赤外線(SWIR)、中波長赤外線(MWIR)、または長波長赤外線(LWIR))において動作することが知られているセンサと、これに対応する光学フィルタアレイを付加的または代替的に使用することができる。これらのセンサおよび光学フィルタアレイは、インジウム・ガリウム・ヒ素(InGaAs)やインジウム・アンチモン(InSb)などの検出器の材質に応じたカメラ、またはマイクロボロメータアレイに応じたカメラを付加的または代替的に備えていてもよい。
本明細書で開示するマルチスペクトルイメージング技術において使用されるイメージセンサは、カラーフィルタアレイ(CFA)などの光学フィルタアレイに併用して使用してもよい。いくつかの種類のCFAは、入射光を赤色(R)、緑色(G)および青色(B)の可視領域に分割し、分割した可視光をイメージセンサ上のそれぞれ専用の赤色、緑色または青色のフォトダイオード感光体に導くことができる。CFAの一般的な例として、ベイヤーパターンがあり、これは、光検出器の長方形のグリッド上にRGBカラーフィルタを特定のパターンで配置したカラーフィルタアレイである。ベイヤーパターンは、緑色50%、赤色25%および青色25%で構成されており、赤色のカラーフィルタと緑色のカラーフィルタが繰り返された列と、青色のカラーフィルタと緑色のカラーフィルタが繰り返された列とが交互に並んでいる。また、いくつかの種類のCFA(たとえばRGB-NIRセンサ用のCFA)も、NIR光を分割し、分割したNIR光をイメージセンサ上のそれぞれ専用のフォトダイオード感光体に導くことができる。
したがって、CFAのフィルタ部品の波長領域により、撮影された画像において各画像チャネルが示す波長領域が決まることがある。したがって、様々な実施形態において、画像の赤色のチャネルは、カラーフィルタの赤色の波長領域に対応してもよく、黄色光および橙色光の一部を含むことができ、約570nm~約760nmまたは570nm~760nmの範囲である。様々な実施形態において、画像の緑色のチャネルは、カラーフィルタの緑色の波長領域に対応してもよく、黄色光の一部を含むことができ、約570nm~約480nmまたは570nm~480nmの範囲である。様々な実施形態において、画像の青色のチャネルは、カラーフィルタの青色の波長領域に対応してもよく、紫色光の一部を含むことができ、約490nm~約400nmまたは490nm~400nmの範囲である。当業者であれば、CFAの色(たとえば、赤色、緑色および青色)を定義する正確な開始波長と終了波長(または電磁スペクトルの一部)が、CFAの実装に応じて様々に異なることがあることを十分に理解できるであろう。
さらに、通常の可視光用CFAは、可視スペクトル外の光を透過する。したがって、多くのイメージセンサでは、赤外線波長を阻止するが、可視光を通過させる薄膜赤外線反射フィルタをセンサの前面に配置することによって赤外線に対する感度を制限している。しかし、本明細書で開示するイメージングシステムのいくつかでは、薄膜赤外線反射フィルタを省略して、赤外線光を通過させてもよい。したがって、赤色のチャネル、緑色のチャネルおよび/または青色のチャネルを使用して、赤外線の波長帯域を収集してもよい。いくつかの実装において、青色のチャネルを使用して、特定のNUV波長帯域を収集してもよい。スタックしたスペクトル画像の各波長において、赤色のチャネル、緑色のチャネルおよび青色のチャネルは、各チャネルが特有の透過効率を有することに関連して異なるスペクトル応答を示すことから、公知の透過プロファイルを使用して、アンミキシングを行う前の各スペクトル帯域に、特有の重み付け応答を提供してもよい。たとえば、この重み付け応答には、赤外線波長領域および紫外線波長領域における赤色のチャネル、青色のチャネルおよび緑色のチャネルの公知の透過応答が含まれていてもよく、これによって、これらの波長領域からの帯域の収集に各チャネルを使用することが可能となる。
以下でさらに詳述するように、イメージセンサに向かう光路に沿って、追加のカラーフィルタをCFAの前に配置することにより、イメージセンサ上に入射する光の特定の帯域を選択的に抽出することができる。本明細書で開示するカラーフィルタは、(薄膜の)ダイクロイックフィルタおよび/もしくは吸収フィルタの組み合わせ、または単独のダイクロイックフィルタ、および/もしくは単独の吸収フィルタである。本明細書で開示するカラーフィルタの一部は、(通過帯域の)特定の領域内の周波数を通すが、(阻止領域において)その領域外の周波数を遮断(減衰)するバンドパスフィルタであってもよい。本明細書で開示するカラーフィルタの一部は、不連続な複数の波長領域を通すマルチバンドパスフィルタであってもよい。これらの「周波数帯」は、CFAフィルタの広い色範囲と比べて、通過帯域が狭く、阻止領域における減衰量が大きく、かつスペクトルのロールオフが急峻であってもよい。なお、ロールオフは、フィルタの通過帯域から阻止領域への移行部におけるスペクトル応答の急峻さとして定義される。たとえば、本明細書で開示するカラーフィルタは、約20nm~約40nmまたは20nm~40nmの通過帯域をカバーすることができる。カラーフィルタのこのような特定の構成によって、実際にセンサに入射する波長帯域が決まることがあり、これによって、本明細書で開示するイメージング技術の精度が向上する場合がある。本明細書に記載のカラーフィルタは、特定の用途に適した特定の波長帯域に応じて、上述の領域のうち、電磁放射線の特定の帯域を選択的に阻止または通過させるように構成することができる。
本明細書において、「画素」という用語は、二次元検出器アレイの素子により生成された出力を説明する場合に使用される。これに対して、二次元検出器アレイのフォトダイオード、すなわち、1つの感光素子は、光電効果を介して光子を電子に変換し、次いでこの電子を画素値の決定に使用可能なシグナルに変換することができるトランスデューサとして挙動する。データキューブの1つの要素は、「ボクセル」(たとえば、体積のある要素)と呼ぶことができる。「スペクトルのベクトル」は、データキューブの特定の(x,y)位置におけるスペクトルデータ(たとえば、物体空間の特定の点から受光した光のスペクトラム)を示すベクトルを指す。本明細書において、データキューブの1つの水平面(たとえば、1つのスペクトル次元を表す画像)は、「画像チャネル」と呼ばれる。本明細書に記載の特定の実施形態では、スペクトルのビデオ情報を撮影してもよく、得られるデータ次元は、NxNyNλNtで示される「ハイパーキューブ」の形態であると想定することができる(ここで、Ntは、ビデオシーケンス中に撮影されたフレームの数である)。
湾曲したマルチバンドパスフィルタを備えるマルチアパーチャイメージングシステムの一例の概要
図3Aは、本開示による、湾曲したマルチバンドパスフィルタを備えるマルチアパーチャイメージングシステム200の一例の概略図を示す。ここに示した概略図は、第1のイメージセンサ領域225A(フォトダイオードPD1~PD3)および第2のイメージセンサ領域225B(フォトダイオードPD4~PD6)を含む。フォトダイオードPD1~PD6は、たとえば、半導体基板(たとえばCMOSイメージセンサ)に形成されたフォトダイオードであってもよい。通常、各フォトダイオードPD1~PD6は、何らかの材質、半導体、センサ素子または入射光を電流に変換できるその他の装置からなる単一ユニットであってもよい。この図では、マルチアパーチャイメージングシステムの構造およびその作動を説明することを目的として、マルチアパーチャイメージングシステム全体のごく一部のみを示しており、実装では、イメージセンサ領域は、何百個または何千個ものフォトダイオード(およびこれに対応するカラーフィルタ)を備えることができることは十分に理解できるであろう。第1のイメージセンサ領域225Aと第2のイメージセンサ領域225Bは、実装に応じて、別々のセンサとして実装してもよく、同じイメージセンサ上の別々の領域として実装してもよい。図3Aでは、2個のアパーチャとこれに対応する光路およびセンサ領域が示されているが、図3Aに示した光学設計原理は、実装に応じて、3個以上のアパーチャとこれに対応する光路およびセンサ領域を含む設計に拡大することができることは十分に理解できるであろう。
マルチアパーチャイメージングシステム200は、第1のセンサ領域225Aに向かう第1の光路を提供する第1の開口部210A、および第2のセンサ領域225Bに向かう第1の光路を提供する第2の開口部210Bを含む。これらのアパーチャは、画像に反映される光の明るさを増加さたり減少させたりするために調整可能であってもよく、あるいは、これらのアパーチャを調整することによって、特定の画像の露光時間を変更し、イメージセンサ領域に入射する光の明るさが変化しないようにしてもよい。これらのアパーチャは、光学設計分野の当業者が妥当であると判断可能であれば、このマルチアパーチャシステムの光軸に沿ったどの位置に配置してもよい。第1の光路に沿って配置した光学部品の光軸を破線230Aで示し、第2の光路に沿って配置した光学部品の光軸を破線230Bで示すが、これらの破線が、マルチアパーチャイメージングシステム200の物理構造を示すわけではないことは十分に理解できるであろう。光軸230Aと光軸230Bは距離Dにより隔てられており、この距離Dは、第1のセンサ領域225Aにより撮影される画像と第2のセンサ領域225Bにより撮影される画像の間の視差となりうる。「視差」は、立体写真画像の左側と右側(または上側と下側)において対応する2つの点の間の距離を指し、物体空間の同じ物理点が各画像の異なる位置に現われる。この視差を補正し、かつ利用した処理技術を以下でさらに詳細に述べる。
光軸230Aおよび光軸230Bは、それぞれに対応するアパーチャの中心Cを通過する。その他の光学部品も、これらの光軸を中心として配置することができる(たとえば、光学部品の回転対称点をこの光軸に沿って配置することができる)。たとえば、第1の光軸230Aを中心として、第1の湾曲マルチバンドパスフィルタ205Aおよび第1の結像レンズ215Aを配置することができ、第2の光軸230Bを中心として、第2の湾曲マルチバンドパスフィルタ205Bおよび第2の結像レンズ215Bを配置することができる。
本明細書において、光学素子の配置に関して述べる際に使用される「上」および「上方」という用語は、物体空間からイメージングシステム200に入った光が、別の構造に到達する(または入射する)前に通り抜ける特定の構造(たとえばカラーフィルタまたはレンズ)の位置を指す。これを説明すると、第1の光路に沿って、湾曲マルチバンドパスフィルタ205Aが、アパーチャ210Aの上方に位置し、アパーチャ210Aが、結像レンズ215Aの上方に位置し、結像レンズ215Aが、CFA 220Aの上方に位置し、CFA 220Aが、第1のイメージセンサ領域225Aの上方に位置する。したがって、まず、物体空間(たとえば、画像化の対象となる物理的空間)から入った光が、第1の湾曲マルチバンドパスフィルタ205Aを通過し、次にアパーチャ210Aを通過し、次に結像レンズ215Aを通過し、次にCFA 220Aを通過し、最後に、第1のイメージセンサ領域225Aに入射する。第2の光路(たとえば、湾曲マルチバンドパスフィルタ205B、アパーチャ210B、結像レンズ215B、CFA 220B、第2のイメージセンサ領域225Bを含む光路)も同様の配置となる。別の実装では、アパーチャ210A,210Bおよび/または結像レンズ215A,215Bを、湾曲マルチバンドパスフィルタ205A,205Bの上方に配置することができる。さらに、別の実装では、物理的なアパーチャを使用しなくてもよく、光学部品のクリアアパーチャに依存して、センサ領域225A上およびセンサ領域225B上に結像される光の明るさを制御してもよい。したがって、レンズ215A,215Bは、アパーチャ210A,210Bと湾曲マルチバンドパスフィルタ205A,205Bの上方に配置してもよい。この実装では、光学設計分野の当業者が必要であると判断可能であれば、アパーチャ210A,210Bをレンズ215A,215Bの上方に配置してもよく、あるいはアパーチャ210A,210Bをレンズ215A,215Bの下方に配置してもよい。
第1のセンサ領域225Aの上方に配置した第1のCFA 220Aと、第2のセンサ領域225Bの上方に配置した第2のCFA 220Bは、特定の波長を選択的に通過させるフィルタとして機能することができ、入射してきた可視領域の光を、赤色領域、緑色領域および青色領域(それぞれR、GおよびBの記号で示す)に分割することができる。選択された特定の波長のみを、第1のCFA 220Aおよび第2のCFA 220Bの各カラーフィルタに通過させることによって、光が「分割」される。分割された光は、イメージセンサ上のそれぞれ専用の赤色ダイオード、緑色ダイオードまたは青色ダイオードにより受光される。赤色フィルタ、青色フィルタおよび緑色フィルタが一般に使用されているが、別の実施形態では、撮影された画像データに必要とされるカラーチャネルに応じて様々なカラーフィルタを使用することができ、たとえば、RGB-IR CFAと同様に、紫外線、赤外線または近赤外線を選択的に通過させるフィルタを使用することができる。
図に示すように、単一のフォトダイオードPD1~PD6の上方にCFAの各フィルタが配置されている。図3Aでは、マイクロレンズ(MLで示す)の一例も示されており、これらのマイクロレンズは、各カラーフィルタ上に形成するか、それ以外の方法で各カラーフィルタ上に配置することができ、アクティブ検出器の領域上に入射光を集光することができる。別の実装では、単一のフィルタの下方に複数個のフォトダイオードを有していてもよい(たとえば、2個、4個またはそれ以上の個数の互いに隣接したフォトダイオードの集合体であってもよい)。図に示した一例において、フォトダイオードPD1とフォトダイオードPD4は、赤色フィルタの下方に位置することから、赤色のチャネルの画素情報を出力し、フォトダイオードPD2とフォトダイオードPD5は、緑色フィルタの下方に位置することから、緑色のチャネルの画素情報を出力し、フォトダイオードPD3とフォトダイオードPD6は、青色フィルタの下方に位置することから、青色のチャネルの画素情報を出力する。さらに、以下で詳述するように、所定のフォトダイオードによる特定の色のチャネル出力は、アクティブな光源および/またはマルチバンドパスフィルタ205A,205Bを通過した特定の周波数帯に基づいて、さらに狭い周波数帯とすることができ、これによって、所定のフォトダイオードから、様々な露光により様々な画像チャネル情報を出力できるようになる。
結像レンズ215A,215Bは、センサ領域225A,225B上で物体シーンの画像の焦点が合うように成形することができる。各結像レンズ215A,215Bは、図3Aに示すような単一の凸面レンズに限定されず、画像形成に必要とされる数の光学素子および表面で構成されていてもよく、市販または特注設計による様々な種類の結像レンズまたはレンズアセンブリを使用することができる。各光学素子またはレンズアセンブリは、止め輪またはベゼルを備えた光学機械式鏡筒を使用して、直列に収容または積み重なるように形成または互いに結合されていてもよい。いくつかの実施形態において、光学素子またはレンズアセンブリは、1つ以上の結合レンズ群を含んでいてもよく、結合レンズ群として、2つ以上の光学部品が接合されたものや、それ以外の方法で2つ以上の光学部品が結合されたものを挙げることができる。様々な実施形態において、本明細書に記載のマルチバンドパスフィルタは、マルチスペクトル画像システムのレンズアセンブリの前面に配置してもよく、マルチスペクトル画像システムの単レンズの前面に配置してもよく、マルチスペクトル画像システムのレンズアセンブリの後ろに配置してもよく、マルチスペクトル画像システムの単レンズの後ろに配置してもよく、マルチスペクトル画像システムのレンズアセンブリの内部に配置してもよく、マルチスペクトル画像システムの結合レンズ群の内部に配置してもよく、マルチスペクトル画像システムの単レンズの表面上に直接配置してもよく、マルチスペクトル画像システムのレンズアセンブリの素子の表面上に直接配置してもよい。さらに、アパーチャ210Aおよびアパーチャ210Bを取り除き、デジタル一眼レフ(DSLR)カメラまたはミラーレスカメラによる写真撮影で通常使用される種類のレンズをレンズ215A,215Bに使用してもよい。さらに、レンズ215A,215Bは、CマウントスレッドまたはSマウントスレッドを装着に使用したマシンビジョンに使用される種類のレンズであってもよい。ピント調整は、たとえば、手動フォーカス、コントラスト方式オートフォーカスまたはその他の適切なオートフォーカス技術に基づいて、センサ領域225A,225Bに対して結像レンズ215A,215Bを移動させることによって、あるいは結像レンズ215A,215Bに対してセンサ領域225A,225Bを移動させることによって行うことができる。
各マルチバンドパスフィルタ205A,205Bは、狭い周波数帯の複数の光を選択的に通過するように構成することができ、たとえば、いくつかの実施形態では10~50nmの周波数帯の複数の光(別の実施形態では、これよりも広い周波数帯または狭い周波数帯の複数の光)を選択的に通過するように構成することができる。図3Aに示すように、マルチバンドパスフィルタ205A,205Bはいずれも周波数帯λc(「共通周波数帯」)を通過させることができる。3つ以上の光路を備えた実装では、各マルチバンドパスフィルタは、この共通周波数帯を通過させることができる。このような方法によって、各センサ領域は、同じ周波数帯(「共通チャネル」)の画像情報を撮影する。以下でさらに詳細に述べるように、この共通チャネルで得られたこの画像情報を使用して、各センサ領域により撮影された一組の画像の位置合わせを行うことができる。いくつかの実装では、1つの共通周波数帯と、これに対応する共通チャネルを有していてもよく、あるいは、複数の共通周波数帯と、これらに対応する共通チャネルを有していてもよい。
各マルチバンドパスフィルタ205A,205Bは、共通周波数帯λcに加えて、1つ以上の特有周波数帯を選択的に通過させるように構成することができる。このような方法によって、イメージングシステム200は、センサ領域205A,205Bにより一括して撮影される異なるスペクトルチャネルの数を、単一のセンサ領域により撮影可能なスペクトルチャネル数よりも多くすることができる。この態様は、特有周波数帯λu1を通過させるマルチバンドパスフィルタ205Aと、特有周波数帯λu2を通過させるマルチバンドパスフィルタ205Bとして、図3Aに示されており、ここでλu1およびλu2は、互いに異なる周波数帯を示す。この図では、2個の周波数帯の通過として示したが、本開示のマルチバンドパスフィルタは、それぞれ2個以上の周波数帯を一度に通過させることができる。たとえば、図11Aおよび図11Bに関連して後述するように、いくつかの実装では、各マルチバンドパスフィルタは、4個の周波数帯を通過させることができる。様々な実施形態において、これよりも多くの周波数帯を通過させてもよい。たとえば、4台のカメラを使用した実装では、8個の周波数帯を通過させるように構成されたマルチバンドパスフィルタを備えていてもよい。いくつかの実施形態において、周波数帯の数は、たとえば、4個、5個、6個、7個、8個、9個、10個、12個、15個、16個またはそれ以上であってもよい。
マルチバンドパスフィルタ205A,205Bは、それぞれに対応するセンサ領域225A,225Bにおける入射角度依存性のスペクトル透過が低下するように選択された弯曲を有する。この結果、物体空間から狭帯域の照明を受光する場合、その波長に感度を有するセンサ領域225A,225Bの表面全体に配置された各フォトダイオード(たとえば、センサ領域を覆い、その波長を通過させるカラーフィルタ)は、図1Aに関連して前述したような、センサの端部近傍に配置されたフォトダイオードにおいて見られる波長のシフトは発生させず、実質的に同じ波長の光を受光すると考えられる。このような構成により、平坦なフィルタを使用した場合よりも正確なスペクトル画像データを生成することができる。
図3Bは、図3Aに示したマルチアパーチャイメージングシステムの1つの光路を構成する光学部品の光学設計の一例を示す。より具体的には、図3Bは、マルチバンドパスフィルタ205A,205Bとして使用することができる特注の色消しダブレット240を示す。この特注の色消しダブレット240は、ハウジング250を通して光を通過させ、イメージセンサ225へと送る。ハウジング250は、前述したような、開口部210A,210Bおよび結像レンズ215A,215Bを備えていてもよい。
色消しダブレット240は、マルチバンドパスフィルタコーティング205A,205Bに必要とされる表面を組み込んで導入されることから、光学収差を補正するように構成される。図に示した色消しダブレット240は、分散度と屈折率が異なるガラスまたはその他の光学材料から作製することが可能な2枚のレンズを含む。別の実装では、3枚以上のレンズを使用してもよい。これらの色消しダブレットレンズの設計は、湾曲した前面242上にマルチバンドパスフィルタコーティング205A,205Bが組み込まれ、フィルタコーティング205A,205Bを蒸着した湾曲単レンズの光学面を組み込むことにより光学収差を打ち消し、色消しダブレット240の湾曲した前面242と湾曲した裏面244の組み合わせ効果により屈折力または集光力を制限し、ハウジング250に収容されたこれらのレンズのみで、集光のための主要な素子を構成することができる。したがって、色消しダブレット240は、イメージングシステム200により撮影された画像データの精度の向上に寄与することができる。これらの個々のレンズは互いに隣接して装着することができ、たとえば、個々のレンズを結合または接合することにより、一方のレンズの収差が他方のレンズの収差によって相殺されるように成形することができる。色消しダブレット240の湾曲した前面242または湾曲した後面244は、マルチバンドパスフィルタコーティング205A,205Bでコーティングすることができる。別のダブレット設計を本明細書に記載のシステムに実装してもよい。
本明細書に記載の光学設計のさらなる変形例を実装してもよい。たとえば、いくつかの実施形態において、図3Bに示すダブレット240の代わりに、図3Aに示す凸レンズ(正レンズ)または凹レンズ(負レンズ)のような、単レンズまたはその他の光学単レンズが光路に含まれていてもよい。図3Cは、平坦なフィルタ252をレンズのハウジング250とセンサ225の間に配置した実装の一例を示す。複数帯域透過特性を備えた平坦なフィルタ252を組み込んで導入した図3Cの色消しダブレット240は、ハウジング250に収容された各レンズの屈折力に顕著に寄与することなく、光学収差を補正する。図3Dは、ハウジング250内に収容されたレンズアセンブリの前面にマルチバンドパスコーティング254を施すことにより、マルチバンドパスコーティングが実装された実装の別の一例を示す。このように、このマルチバンドパスコーティング254は、ハウジング250内に収容されたどの光学素子のどの曲面に施してもよい。
図4A~4Eは、図3Aおよび図3Bに関連して説明した光学設計を備えたマルチスペクトル・マルチアパーチャイメージングシステム300の一実施形態を示す。図4Aは、より具体的には、内部の部品を見せるために半透明で示したハウジング305を備えるイメージングシステム300の透視図を示す。実際のハウジング305は、たとえば、所望とする組み込み用コンピュータリソースの大きさに応じて、図に示したハウジング305よりも大きくてもよく、小さくてもよい。図4Bは、イメージングシステム300の正面図を示す。図4Cは、図4BのC-C線に沿って切断した、イメージングシステム300の破断側面図を示す。図4Dは、処理基板335を示した、イメージングシステム300の底面図を示す。以下、図4A~4Dについて述べる。
イメージングシステム300のハウジング305は、別のハウジングに収容されていてもよい。たとえば、携帯型の実装では、イメージングシステム300を安定に保持できるように成形された1本以上のハンドルを備えていてもよいハウジングの内部に該イメージングシステムが収容されていてもよい。携帯型の実装の一例は、図18A~18Cおよび図19A~19Bにより詳細に示している。ハウジング305の上面には、4つの開口部320A~320Dを備えている。各開口部320A~320Dの上方に互いに異なるマルチバンドパスフィルタ325A~325Dが配置されており、フィルタキャップ330A~330Bにより適所に収容されている。マルチバンドパスフィルタ325A~325Dは湾曲していてもよく、本明細書で述べるように、各マルチバンドパスフィルタは、共通周波数帯と少なくとも1つの特有周波数帯を通過させることによって、カラーフィルタアレイで覆われたイメージセンサで通常撮影されるスペクトルチャネルよりも多くのスペクトルチャネルにおいて高精度マルチスペクトルイメージングを実施することができる。前述したイメージセンサ、結像レンズおよびカラーフィルタは、カメラハウジング345A~345Dの内部に配置される。いくつかの実施形態において、たとえば、図20A~20Bに示すように、前述したイメージセンサ、結像レンズおよびカラーフィルタは、単一のカメラハウジングに収容されていてもよい。このように、図に示した実装では、別々の複数のセンサを使用する(たとえば、カメラハウジング345A~345Dのそれぞれの内部に1個のセンサを収容する)が、別の実装では、開口部320A~320Dを介して露出した領域全体に配置された単一のイメージセンサを使用できることも十分に理解できるであろう。この実施形態では、支持材340を使用して、カメラハウジング345A~345Dが、前記イメージングシステムのハウジング305に固定されているが、様々な実装において、その他の適切な手段を使用して固定することもできる。
ハウジング305の上面は、任意で設けられる照明基板310を支持しており、この照明基板310は光学拡散素子315で覆われている。照明基板310は、以下の図4Eに関連してより詳細に説明する。拡散素子315は、ガラス、プラスチックまたはその他の光学材料で構成することができ、照明基板310から発せられた光を拡散させて、物体空間が空間的に実質的に均一な照明を受けられるようにする。対象物体に均一な照明を当てると、各波長において、物体の表面全体に実質的に均一な量の照明を当てることができるため、たとえば画像撮影による組織の臨床分析などの特定のイメージング用途において有益な場合がある。いくつかの実施形態において、本明細書で開示するイメージングシステムは、任意で設けられる照明基板からの光の代わりに、またはこの照明基板からの光に加えて、環境光を利用してもよい。
使用中の照明基板310は熱を生じることから、イメージングシステム300は、放熱フィン355を多数備えたヒートシンク350を備える。放熱フィン355は、各カメラハウジング345A~345Dの間の空間内へと延伸することができ、ヒートシンク350の上部は、照明基板310からフィン355へと熱を逃がすことができる。ヒートシンク350は、適切な熱伝導性材料から作製することができる。ヒートシンク350は、その他の部品からの熱の放散を補助してもよく、これによって、イメージングシステムの実装のいくつかでは、ファンを設けなくてもよい。
ハウジング305内の多数の支持材365によって、カメラ345A~345Dと通信を行う処理基板335が固定されている。処理基板335は、イメージングシステム300の作動を制御することができる。図には示していないが、イメージングシステム300は、たとえば該イメージングシステムを使用して作成されたデータの格納用メモリなどの、1つ以上のメモリ、および/またはコンピュータが実行可能な命令を出すシステム制御用モジュールを備えるように構成することもできる。処理基板335は、システム設計の目的に応じて様々に構成することができる。たとえば、(たとえば、コンピュータが実行可能な命令を出すモジュールにより)、照明基板310の特定のLEDの起動を制御するように処理基板を構成することができる。いくつかの実装では、高度に安定な同期整流降圧型LEDドライバを使用することができ、このドライバによって、ソフトウェアによるLEDのアナログ電流制御と、LEDの故障の検出が可能となる。いくつかの実装では、(たとえば、コンピュータが実行可能な命令を出すモジュールにより)、処理基板335または別の処理基板に画像データ分析機能をさらに提供することができる。図には示していないが、イメージングシステム300は、処理基板335がセンサからデータを受信し、そのデータを処理できるように、センサと処理基板335の間にデータインターコネクトを備えていてもよく、処理基板が照明基板310の特定のLEDを起動させることができるように、照明基板310と処理基板335の間にデータインターコネクトを備えていてもよい。
図4Eは、イメージングシステム300に設けられていてもよい照明基板310の一例を示し、他の部品から照明基板310を分離して示している。照明基板310は、中央部から延びた4つのアームを備え、各アームに沿ってLEDが3列に配置されている。隣り合う列のLED間のスペースは、隣り合うLED同士が離れるように、横方向の中心が互いにずれるように配置されている。各LED列は、様々な色のLEDを有する行を多数含む。中央部の各角部に1個の緑色LEDが配置されるように、計4個の緑色LED 371が中央部に配置されている。最も内側の行(たとえば、中心部に最も近い行)から説明すると、深紅色LED 372を2個含む行が各列に含まれている(合計8個の深紅色LED)。径方向外側に説明を続けると、各アームは、中心の列に、1個の琥珀色LED 374を含む行を有し、両側の列に、2個の短く点滅する青色LED 376を含む行を有し(合計8個の短く点滅する青色LEDを有し)、中心の列に、1個の琥珀色LED 374を含む行をもう1つ有し(合計8個の琥珀色LEDを有し)、両側の列に、1個の非PPG用NIR LED 373と1個の赤色LED 375を有する行を有し(合計4個の非PPG用NIR LEDと合計4個の赤色LEDを有し)、中心の列に、1個のPPG用NIR LED 377を有する(合計4個のPPG用NIR LEDを有する)。「PPG用」LEDとは、生体組織の拍動血流を表す光電容積脈波(PPG)情報を撮影するための連続した多数回の露光において起動されるLEDを指す。別の実施形態では、様々なその他の色および/または配置のLEDを照明基板において使用してもよいことを理解できるであろう。
図5は、図3Aおよび図3Bに関連して説明した光学設計を備えたマルチスペクトル・マルチアパーチャイメージングシステム300の別の一実施形態を示す。イメージングシステム300の設計と同様に、イメージングシステム400は、4つの光路を含み、この図では、止め輪430A~430Dでハウジング405に固定されたマルチバンドパスフィルタレンズ群425A~425Dを有する開口部420A~420Dとして示している。イメージングシステム400は、対象物体上への空間的に均一な光の照射を補助するため、止め輪430A~430Dの間の、ハウジング405の正面部に固定された照明基板410と、照明基板410上に配置された拡散板415とをさらに含む。
イメージングシステム400の照明基板410は、十字形の4つの腕にLEDが配置されており、各腕には密集したLEDが2列に配置されている。このようにすれば、前述の照明基板310よりも照明基板410がコンパクトになり、より小さいフォームファクタ要件を有するイメージングシステムとの使用に適する場合がある。この一例としての構成において、各腕は、最外列に1個の緑色LEDと1個の青色LEDを有し、内側に向かって、黄色LEDを2行含み、橙色LEDを1行含み、1個の赤色LEDと1個の深紅色LEDを1行に含み、1個の琥珀色LEDと1個のNIR LEDを1行に含む。したがって、この実装では、より長い波長の光を発するLEDが照明基板410の中心部に位置し、より短い波長の光を発するLEDが照明基板410の端部に位置するように、各LEDが配置される。
図6A~6Cは、図3Aおよび図3Bに関連して説明した光学設計を備えたマルチスペクトル・マルチアパーチャイメージングシステム500の別の一実施形態を示す。より具体的には、図6Aは、イメージングシステム500の透視図を示し、図6Bは、イメージングシステム500の正面図を示し、図6Cは、図6BのC-C線に沿って切断した、イメージングシステム500の破断側面図を示す。イメージングシステム500は、イメージングシステム300に関連して前述したものと同様の部品(たとえば、ハウジング505、照明基板510、拡散板515、止め輪530A~530Dにより開口部の上方に固定されたマルチバンドパスフィルタ525A~525Dなど)を含むが、短いフォームファクタとなっている(たとえば、一実施形態において、組み込み用コンピュータ部品の数が少なくなっているか、かつ/または組み込み用コンピュータ部品の寸法が小さくなっている)。イメージングシステム500は、カメラアライメントの位置を強固かつ堅牢に固定するための直接カメラ-フレームマウント540をさらに含む。
図7A~7Bは、マルチスペクトル・マルチアパーチャイメージングシステム600の別の一実施形態を示す。図7A~7Bは、マルチアパーチャイメージングシステム600の周辺光源610A~610Cの別の可能な配置を示す。図に示すように、図3A~3Dに関連して説明した光学設計を備えたマルチバンドパスフィルタ625A~625Dを備えた4つのレンズアセンブリは、長方形または正方形に配置して、(イメージセンサを含む)4台のカメラ630A~630Dに受光させることができる。3つの長方形発光素子610A~610Cは、マルチバンドパスフィルタ625A~625Dを備えた各レンズアセンブリの外部かつ該レンズアセンブリ間に、互いに平行になるように配置することができる。これらの長方形発光素子は、広域スペクトルの発光パネルであってもよく、それぞれ異なる周波数帯の光を発するLEDが配置されたものであってもよい。
図8A~8Bは、マルチスペクトル・マルチアパーチャイメージングシステム700の別の一実施形態を示す。図8A~8Bは、マルチアパーチャイメージングシステム700の周辺光源710A~710Dの別の可能な配置を示す。図に示すように、図3A~3Dに関連して説明した光学設計を採用したマルチバンドパスフィルタ725A~725Dを備えた4つのレンズアセンブリは、長方形または正方形に配置して、(イメージセンサを含む)4台のカメラ730A~730Dに受光させることができる。4台のカメラ730A~730Dは、互いにより接近した構成の一例として示されており、これによって、レンズ間の遠近差を最小限に抑えてもよい。4つの長方形発光素子710A~710Dは、マルチバンドパスフィルタ725A~725Dを備えたレンズアセンブリを取り囲んだ正方形に配置することができる。これらの長方形発光素子は、広域スペクトルの発光パネルであってもよく、それぞれ異なる周波数帯の光を発するLEDが配置されたものであってもよい。
図9A~9Cは、マルチスペクトル・マルチアパーチャイメージングシステム800の別の一実施形態を示す。イメージングシステム800は、レンズ集合体前面フレーム830に接続されたフレーム805を含み、レンズ集合体前面フレームは、開口部820とマイクロビデオレンズ825用の支持構造を含み、マイクロビデオレンズ825は、図3A~3Dに関して説明した光学設計を使用したマルチバンドパスフィルタを備えることができる。マイクロビデオレンズ825は、レンズ集合体後面フレーム840に装着された4台のカメラ845(結像レンズおよびイメージセンサ領域を含む)に受光させる。直線状に4列に配置されたLED 811は、レンズ集合体前面フレーム830の4つの辺に沿って配置されており、4つの辺のそれぞれには各LED用の拡散素子815が備えられている。図9Bおよび図9Cは、マルチアパーチャイメージングシステム800の可能な寸法の一例をインチ単位で示したものである。
図10Aは、図3A~3Dに関連して説明した光学設計を備えたマルチスペクトル・マルチアパーチャイメージングシステム900の別の一実施形態を示す。イメージングシステム900は、モバイル端末910のマルチアパーチャカメラ915に取り付け可能な1組のマルチバンドパスフィルタ905として実装することができる。たとえば、スマートフォンなどの特定のモバイル端末910は、2つのイメージセンサ領域に通じる2つの開口部を有する立体イメージングシステムを装備することができる。本明細書で開示するマルチアパーチャスペクトルイメージング技術は、狭い周波数帯の複数の光をセンサ領域へと通過させる適切な1組のマルチバンドパスフィルタ905を提供することによって、このような装置に実装することができる。この1組のマルチバンドパスフィルタ905には、このような周波数帯の光を物体空間に照射する光源(LEDアレイと拡散板など)を装備してもよい。
イメージングシステム900は、マルチスペクトルデータキューブを生成する処理と、(たとえば、臨床組織分類用アプリケーション、生体認証用アプリケーション、材料分析用アプリケーションまたはその他のアプリケーションにおける)マルチスペクトルデータキューブの処理とを行うようにモバイル端末を構成するモバイルアプリケーションをさらに含むことができる。あるいは、このモバイルアプリケーションは、ネットワークを介して、離れた場所にある処理システムにマルチスペクトルデータキューブを送信し、その後、分析結果を受信して表示するようにモバイル端末910を構成するものであってもよい。このような用途のユーザインタフェース910の一例を図10Bに示す。
図11A~11Bは、図3A~10Bに示したマルチスペクトル・マルチアパーチャイメージングシステムに4種のフィルタを実装した場合に、これらのフィルタを通過することができる1組の周波数帯の一例を示し、これらの1組の周波数帯は、たとえば、べイヤーCFA(または別のRGB CFAまたはRGB-IR CFA)を備えたイメージセンサに入射する。
Figure 2022513847000002
これらの透過曲線には、この例で使用したセンサによる量子効率の効果も含んでいる。図に示すように、この4台1組のカメラ全体で、8個の特有チャネルすなわち8個の特有周波数帯が撮影される。各フィルタは、各カメラに対して2つの共通周波数帯(左端の2つのピーク)と、別の2つの周波数帯を通過させる。この実装では、第1のカメラおよび第3のカメラは、第1の共有NIR周波数帯の光(右端のピーク)を受光し、第2のカメラおよび第4のカメラは、第2の共有NIR周波数帯の光(右から2番目のピーク)を受光する。各カメラは、約550nm~約800nmまたは550nm~800nmの範囲の特有周波数帯をそれぞれ1つずつ受光する。したがって、これらのカメラは、コンパクトな構成を使用して、8個の特有のスペクトルチャネルを撮影することができる。図11Bのグラフ1010は、図11Aに示した4台のカメラの照明として使用してもよい、図4Eに示したLEDボードの分光放射照度を示す。
この実装では、臨床組織分類に適したスペクトルチャネルの生成に基づき、8個の周波数帯を選択したが、(イメージングシステム内に熱を伝える)LEDの数を制限しながら、シグナル/ノイズ比(SN比)およびフレームレートに関して最適化してもよい。組織(たとえば、ヒト組織などの動物組織)は、緑色の波長や赤色の波長よりも青色の波長において高いコントラストを示すことから、この8個の周波数帯には、4枚のフィルタすべてを通過する青色光の共通周波数帯(グラフ1000の左端のピーク)を含めている。より具体的には、グラフ1000に示すように、ヒト組織は、約420nmを中心とする周波数帯で画像化した場合にコントラストが最も高くなる。この共通周波数帯に対応するチャネルは視差補正に使用されることから、コントラストを高くするほど、より正確な補正行うことができる。たとえば、視差補正を行う際、画像処理プロセッサは、ローカル法またはグローバル法を使用して、ローカル画像パッチ間またはローカル画像間の類似性を示す性能指数が最大となるように、1組の視差を見出すことができる。あるいは、画像処理プロセッサは、同様の方法を使用して、相違性を示す性能指数を最小にすることができる。これらの性能指数は、エントロピー、相関性、絶対差または深層学習法に基づくものであってもよい。グローバル法による視差の計算は、反復して実行することができ、性能指数が安定になった時点で終了することができる。ローカル法は、画素ごとに視差を計算することができ、1つの画像中の固定されたパッチを性能指数の入力として使用し、試験中に、別の画像から切り出した多数のパッチに対してそれぞれ異なる視差の数値を求める。このような方法はいずれも、考慮される視差の範囲が拘束される場合がある。この拘束は、たとえば、物体の深さおよび距離に関する知識に基づく場合がある。また、物体において予想される傾きの範囲に基づいて拘束される場合もある。計算された視差は、エピポーラ拘束などの射影幾何学により拘束される場合がある。低い解像度で計算された視差の出力を、次のレベルの解像度での視差の計算に対する初期値または拘束として使用して、様々な解像度で視差を計算することができる。たとえば、1回目の計算において4画素の解像度で計算された視差を利用して、より高い解像度での次の視差の計算において±4画素の拘束を設定することができる。視差から計算されるあらゆるアルゴリズムは、コントラストが高いほど恩恵を受け、特に、そのコントラスト源がすべての視点と関連している場合に、高コントラストによる恩恵は顕著である。通常、共通周波数帯は、特定の用途において画像化されることが予想される材料の最も高いコントラストでの画像化に対応できるように選択することができる。
画像の撮影後の、隣接したチャネル間の色分解は完全ではない場合があることから、この実装では、青色の周波数帯に隣接した緑色の周波数帯として、グラフ1000に示したすべてのフィルタを通過する追加の共通周波数帯をさらに設ける。このような周波数帯を利用する理由として、青色カラーフィルタの画素は、スペクトル通過帯域が広いことから、緑色のスペクトラム領域に感度を示すことが挙げられる。これは、通常、スペクトルの重複として示され、隣接したRGB画素間での意図的なクロストークと見なしてもよい。このスペクトルの重複によって、カラーカメラのスペクトル感度をヒト網膜のスペクトル感度に近づけることができ、最終的に得られる色空間をヒトの視覚に質的に似たものとすることができる。したがって、各フィルタに共通する緑色のチャネルを設けることにより、受光した青色光に真に相当する、青色のフォトダイオードによって生成されたシグナルの一部を分離することが可能になる。この分離は、マルチバンドパスフィルタの透過率(Tの凡例で示した黒色の実線で示す)と、これに対応するCFAカラーフィルタの透過率(それぞれQの凡例で示した、赤色の破線、緑色の破線および青色の破線で示す)を要因に入れたスペクトルアンミキシングアルゴリズムを使用して達成することができる。いくつかの実装では、赤色光を共通周波数帯として使用してもよく、この場合、第2の共通チャネルは必要とされない場合もあることは十分に理解できるであろう。
図12は、高解像度スペクトルイメージング能を備えたコンパクト型イメージングシステム1100の一例のハイレベルブロック図を示し、このイメージングシステム1100は、マルチアパーチャスペクトルカメラ1160と光源1165に連結されたプロセッサ1120を含む1組の部品を備える。ワーキングメモリ1105、記憶装置1110、電子ディスプレイ1125およびメモリ1130もプロセッサ1120と通信可能な状態である。本明細書で述べるように、このイメージングシステム1100は、マルチアパーチャスペクトルカメラ1160の開口部の上方にそれぞれに異なるマルチバンドパスフィルタを配置することによって、イメージセンサのCFAに異なるカラーフィルタを配置した場合よりも多くの画像チャネルを撮影可能であってもよい。
コンパクト型イメージングシステム1100は、携帯電話、デジタルカメラ、タブレット型コンピュータ、パーソナルデジタルアシスタントなどの装置であってもよい。コンパクト型イメージングシステム1100は、画像撮影用の内部カメラまたは外部カメラを使用したデスクトップパソコンやテレビ会議ステーションなどの固定型デバイスであってもよい。コンパクト型イメージングシステム1100は、画像撮影装置と、これとは別個に設けられた、画像撮影装置からの画像データを受信する処理装置の組み合わせであってもよい。ユーザは、このコンパクト型イメージングシステム1100上で複数のアプリケーションを利用可能であってもよい。これらのアプリケーションには、従来の写真用アプリケーション、静止画像および動画の撮影用アプリケーション、動的色補正用アプリケーション、明度・濃淡の補正用アプリケーションなどが含まれていてもよい。
この画像撮影システム1100は、画像撮影用のマルチアパーチャスペクトルカメラ1160を含む。マルチアパーチャスペクトルカメラ1160は、たとえば、図3A~10Bに示す装置のいずれであってもよい。様々なスペクトルチャネルにおいて撮影された画像を、様々なセンサ領域から画像処理プロセッサ1120に送信できるように、マルチアパーチャスペクトルカメラ1160がプロセッサ1120に接続されていてもよい。また、以下で詳細に説明するように、単一または複数の光源1165をプロセッサで制御することによって、特定の露光中に特定の波長の光を発することができる。画像処理プロセッサ1120は、視差が補正された高品質なマルチスペクトルデータキューブを出力できるように、受信した撮影画像に対して様々な処理を行うように構成されていてもよい。
プロセッサ1120は、汎用処理ユニットまたはイメージング用途に特別に設計されたプロセッサであってもよい。図に示すように、プロセッサ1120は、メモリ1130とワーキングメモリ1105に接続されている。この図に示した実施形態では、メモリ1130は、撮影制御モジュール1135、データキューブ作成モジュール1140、データキューブ分析モジュール1145およびオペレーティングシステム1150を格納している。これらのモジュールは、様々な画像処理および装置管理タスクを実行できるようにプロセッサを構成する命令を含む。プロセッサ1120は、ワーキングメモリ1105を使用して、メモリ1130の各モジュールに含まれるプロセッサ命令のワーキングセットを格納する。あるいは、プロセッサ1120は、ワーキングメモリ1105を使用して、画像撮影システム1100の作動中に作成された動的データを格納する。
前述したように、プロセッサ1120は、メモリ1130に格納された数個のモジュールによって構成されている。いくつかの実装において、撮影制御モジュール1135は、マルチアパーチャスペクトルカメラ1160の焦点位置を調整できるようにプロセッサ1120を構成する命令を含む。撮影制御モジュール1135は、たとえば、様々なスペクトルチャネルで撮影されるマルチスペクトル画像や、同じスペクトルチャネル(たとえばNIRチャネル)で撮影されるPPG画像などの画像をマルチアパーチャスペクトルカメラ1160で撮影できるようにプロセッサ1120を構成する命令をさらに含む。非接触PPGイメージングでは、通常、近赤外線(NIR)波長を照明として使用し、この波長において組織への光子の透過が増加することを利用している。したがって、撮影制御モジュール1135、マルチアパーチャスペクトルカメラ1160およびワーキングメモリ1105を備えたプロセッサ1120は、一連のスペクトル画像および/または連続した画像を撮影する手段の1つである。
データキューブ作成モジュール1140は、異なるセンサ領域にある各フォトダイオードから受信した強度シグナルに基づいてマルチスペクトルデータキューブを生成するようにプロセッサ1120を構成する命令を含む。たとえば、データキューブ作成モジュール1140は、すべてのマルチバンドパスフィルタを通過した共通周波数帯に対応するスペクトルチャネルに基づいて、画像化された物体の同じ領域間の視差を推定することができ、この視差を使用して、撮影されたすべてのチャネルおけるすべてのスペクトル画像を互いに位置合わせすることができる(たとえば、すべてのスペクトルチャネルにおいて、特定の物体上の同じ点が実質的に同じ画素位置(x,y)となるように位置合わせを行うことができる)。位置合わせした画像を統合してマルチスペクトルデータキューブを形成し、視差の情報を使用して、別々に画像化された物体の深さを測定してもよく、たとえば、正常組織の深さと創傷部位の最も深い位置の深さとの差を測定してもよい。いくつかの実施形態において、データキューブ作成モジュール1140は、さらにスペクトルアンミキシングを行うことにより、たとえば、フィルタの透過率およびセンサの量子効率を要因に入れたスペクトルアンミキシングアルゴリズムに基づいて、フォトダイオードの強度シグナルのどの部分が、フィルタを通過したどの周波数帯に対応するのかを特定してもよい。
データキューブ分析モジュール1145は、データキューブ作成モジュール1140によって生成されたマルチスペクトルデータキューブを用途に応じて分析するために、様々な技術を実装することができる。たとえば、データキューブ分析モジュール1145の実装のいくつかでは、特定の状態に応じて各画素を分類するように訓練された機械学習モデルに、(深さ情報が含まれていてもよい)マルチスペクトルデータキューブを提供することができる。この特定の状態は、組織イメージングを行う際の臨床状態であってもよく、たとえば、熱傷の状態(たとえば、I度熱傷、II度熱傷、III度熱傷または正常組織という種類)、創傷の状態(たとえば、止血、炎症、増殖、リモデリングまたは正常な皮膚という種類)、治癒可能性(たとえば、特定の治療法を行った場合あるいは特に治療法を行わなかった場合に、組織が創傷状態から治癒する可能性を反映したスコア)、潅流状態、がん状態、または創傷に関連したその他の組織状態であってもよい。また、データキューブ分析モジュール1145は、生体認証および/または材料分析を目的としてマルチスペクトルデータキューブを分析することもできる。
オペレーティングシステムモジュール1150は、画像撮影システム1100のメモリおよび処理リソースを管理するようにプロセッサ1120を構成する。たとえば、オペレーティングシステムモジュール1150は、電子ディスプレイ1125、記憶装置1110、マルチアパーチャスペクトルカメラ1160、光源1165などのハードウェアリソースを管理するデバイスドライバを含んでいてもよい。したがって、いくつかの実施形態において、前述の画像処理モジュールに含まれる命令は、これらのハードウェアリソースと直接的に相互作用しなくてもよく、その代わりに、オペレーティングシステムコンポーネント1150に含まれる標準サブルーチンまたはAPIを介して相互作用してもよい。次に、オペレーティングシステム1150内の命令は、これらのハードウェアコンポーネントと直接相互作用してもよい。
さらに、プロセッサ1120は、ディスプレイ1125を制御して、撮影された画像および/またはマルチスペクトルデータキューブの分析結果(たとえば、分類された画像)をユーザに対して表示するように構成されていてもよい。ディスプレイ1125は、マルチアパーチャスペクトルカメラ1160を含むイメージング装置の外部に設けられてもよく、イメージング装置の一部であってもよい。また、ディスプレイ1125は、画像の撮影前にビューファインダーをユーザに提供するように構成されていてもよい。さらに、ディスプレイ1125は、LCDスクリーンまたはLEDスクリーンを含んでいてもよく、タッチ感応技術を実装していてもよい。
プロセッサ1120は、たとえば、撮影された画像、マルチスペクトルデータキューブおよびデータキューブの分析結果を表すデータなどのデータを記憶装置モジュール1110に書き込んでもよい。この図において、記憶装置モジュール1110は、従来のディスク装置として示されているが、当業者であれば、記憶装置モジュール1110が、どのような種類の記憶媒体装置として構成されていてもよいことを理解できるであろう。たとえば、記憶装置モジュール1110は、フロッピーディスクドライブ、ハードディスクドライブ、光ディスクドライブ、光磁気ディスクドライブなどのディスクドライブ、またはFLASHメモリ、RAM、ROMおよび/もしくはEEPROMなどのソリッドステートメモリを含んでいてもよい。記憶装置モジュール1110は、複数のメモリユニットをさらに含むことができ、メモリユニットのいずれか1つは、画像撮影装置1100の内部に配置されるように構成されていてもよく、画像撮影システム1100の外部に配置されていてもよい。たとえば、記憶装置モジュール1110は、画像撮影システム1100内に格納されたシステムプログラム命令を含むROMメモリを含んでいてもよい。さらに、記憶装置モジュール1110は、カメラから取り外し可能な、撮影された画像を格納するように構成されたメモリカードまたは高速メモリを含んでいてもよい。
図12では、プロセッサ、イメージセンサおよびメモリを備えるように、別々の部品を含むシステムが示されているが、当業者であれば、これらの別々の部品を様々に組み合わせて特定の設計目標を達成できることを認識できるであろう。たとえば、別の一実施形態において、コストの削減および性能の向上のために、メモリ部品をプロセッサ部品と組み合わせてもよい。
さらに、図12では、2つのメモリ部品として、いくつかのモジュールを含むメモリ部品1130と、ワーキングメモリを含む別個のメモリ1105が示されているが、当業者であれば、いくつかの実施形態において、別のメモリ構造を利用できることを認識できるであろう。たとえば、メモリ1130に含まれる各モジュールを実装したプロセッサ命令の格納にROMまたはスタティックRAMメモリを利用した設計としてもよい。あるいは、画像撮影システム1100に一体化されたディスク記憶装置、または外部デバイスポートを介して接続されたディスク記憶装置からシステム起動する際に、プロセッサ命令を読み込んでもよい。次に、このプロセッサ命令をRAMにロードして、プロセッサによる命令の実行を促してもよい。たとえば、ワーキングメモリ1105はRAMメモリであってもよく、プロセッサ1120により命令を実行する前に、ワーキングメモリ1105に命令をロードしてもよい。
画像処理技術の一例の概要
図13は、図3A~10Bおよび図12に示したマルチスペクトル・マルチアパーチャイメージングシステムを使用して、画像データを撮影するプロセス1200の一例を示したフローチャートである。図13は、本明細書に記載のマルチスペクトルデータキューブの生成に使用可能な4種の露光、すなわち、可視光露光1205、追加の可視光露光1210、非可視光露光1215および環境光露光1220の一例を示す。これら露光はどのような順序で撮影してもよく、これらの露光のうちのいくつかは、後述する特定のワークフローから省いてもよく、後述する特定のワークフローに追加してもよいことは十分に理解できるであろう。以下、図11Aおよび図11Bに示した周波数帯を参照しながらプロセス1200を説明するが、別の1組の周波数帯に基づいて生成した画像データを使用して同様のワークフローを実装することもできる。さらに、様々な実施形態において、公知の様々なフラットフィールド補正技術に従って、フラットフィールド補正をさらに実装することによって、画像収集および/または視差補正を向上させてもよい。
可視光露光1205を得るため、照明基板に制御信号を送ることにより、最初の5つのピーク(図11Aのグラフ1000の可視光に対応する左側の5つのピーク)に対応する複数のLEDを点灯させることができる。特定の複数のLEDにおいて同時に光の出力波形を安定化させる必要がある場合がある(たとえば10ミリ秒間)。その後、撮影制御モジュール1135により、4台のカメラの露光を開始させ、たとえば約30ミリ秒間にわたって、この露光を継続することができる。その後、撮影制御モジュール1135により露光を停止させ、(たとえば、フォトダイオードからの生の強度シグナルをワーキングメモリ1105および/またはデータストア1110に転送することにより)センサ領域からデータを取得する。このデータは、本明細書に記載の、視差補正に使用される共通スペクトルチャネルを含んでいてもよい。
SN比を向上させるため、いくつかの実装では、可視光露光1205について述べたのと同じプロセスを使用して、追加の可視光露光1210を撮影することができる。同一またはほぼ同一の2つの露光を撮影してSN比を向上させることにより、画像データをより正確に分析することができる。しかし、単一の画像のSN比のみで許容可能な実装では、この2つ目の露光を省略してもよい。また、いくつかの実装では、共通スペクトルチャネルで2つの露光を撮影することにより、より正確な視差補正を行ってもよい。
いくつかの実装では、NIR光またはIR光に対応する非可視光露光1215をさらに撮影することができる。たとえば、撮影制御モジュール1135は、図11Aに示した2つのNIRチャネルに対応する2種のNIR LEDを起動することができる。特定の複数のLEDにおいて同時に光の出力波形を安定化させる必要がある場合がある(たとえば10ミリ秒間)。その後、撮影制御モジュール1135により、4台のカメラの露光を開始させ、たとえば約30ミリ秒間にわたって、この露光を継続することができる。その後、撮影制御モジュール1135により露光を停止させ、(たとえば、フォトダイオードからの生の強度シグナルをワーキングメモリ1105および/またはデータストア1110に転送することにより)センサ領域からデータを取得する。この露光では、露光1205および露光1210に対する物体の形状または位置に変化はないと見なしても問題はないことから、前もって算出した視差値を使用して各NIRチャネルの位置合わせを行うことができるため、すべてのセンサ領域を通過する共通周波数帯は設けなくてもよい。
いくつかの実装では、拍動血流による組織部位の変形を示すPPGデータを生成するために、多重露光を使用して連続撮影することができる。いくつかの実装では、このようなPPG用露光を非可視光波長で撮影してもよい。PPGデータとマルチスペクトルデータを組み合わせることによって、特定の医用イメージング分析の正確度を向上させてもよいが、PPGデータの撮影によって、画像撮影プロセスに要する時間が長くなる場合がある。いくつかの実装では、画像撮影プロセスに要する時間が長くなったことによって、その間に携帯型撮像装置および/または物体が移動し、誤差が生じることがある。したがって、特定の実装では、PPGデータの撮影を省略してもよい。
いくつかの実装では、環境光露光1220をさらに撮影することができる。この露光では、すべてのLEDを消灯し、周囲の照明(たとえば、日光やその他の光源からの光)を利用して画像を撮影することができる。その後、撮影制御モジュール1135により、4台のカメラの露光を開始させ、所望の時間(たとえば約30ミリ秒間)にわたって、この露光を継続することができる。その後、撮影制御モジュール1135により露光を停止させ、(たとえば、フォトダイオードからの生の強度シグナルをワーキングメモリ1105および/またはデータストア1110に転送することにより)センサ領域からデータを取得する。可視光露光1205(または第2の露光1210によりSN比を補正した可視光露光1205)の数値から環境光露光1220の強度値を差し引き、さらに非可視光露光1215の数値からも環境光露光1220の強度値を差し引くことによって、マルチスペクトルデータキューブから環境光の影響を除去することができる。このような処理により、光源により発せられた光および物体/組織部位から反射した光を示す生成信号の一部を分離して、下流の分析の正確度を向上させることができる。いくつかの実装では、可視光露光1205,1210と非可視光露光1215のみを使用するだけで十分な分析正確度が得られる場合、この工程を省略してもよい。
上に挙げた特定の露光時間は、一実装における一例であり、別の実装では、イメージセンサ、光源の強度および画像化される物体に応じて、露光時間は様々に異なりうることを十分に理解できるであろう。
図14は、たとえば、図3A~10Bおよび図12に示したマルチスペクトル・マルチアパーチャイメージングシステムを使用し、かつ/または図13に示したプロセス1200を使用して撮影された画像データなどの画像データを処理するためのワークフロー1300の略ブロック図を示す。ワークフロー1300は、2つのRGBセンサ領域1301A,1301Bの出力を示しているが、ワークフロー1300のセンサ領域の数はさらに増やすことができ、異なるCFAカラーチャネルに対応するセンサ領域を設けることもできる。
2つのセンサ領域1301A,1301Bから出力されたRGBセンサ出力は、2Dセンサ出力モジュール1305A,1305Bにそれぞれ格納される。これらのセンサ領域の数値は、非線形マッピングモジュール1310A,1310Bに送信され、これらの非線形マッピングモジュール1310A,1310Bにおいて、共通チャネルを使用して撮影された画像間の視差を特定し、次に、この決定された視差をすべてのチャネルに対して適用してすべてのスペクトル画像を互いに位置合わせすることにより視差補正を行うことができる。
次に、これらの非線形マッピングモジュール1310A,1310Bの出力は、深さ計算モジュール1335に提供され、深さ計算モジュール1335は、画像データ中の特定の目的領域の深さを算出することができる。たとえば、深さは、物体とイメージセンサの間の距離を表してもよい。いくつかの実装では、複数の深さの値を算出し、これらを比較することによって、イメージセンサ以外の何らかの対象に対する物体の深さを測定することができる。たとえば、創床の最も深い深さを測定することができ、この創床を取り囲む正常組織の深さ(最も深い深さ、最も浅い深さまたは平均の深さ)も測定することができる。創床の深さから正常組織の深さを差し引くことによって、創傷の最も深い深さを測定することができる。この深さの比較は、創床中のその他の箇所(たとえば、所定のサンプリング点のすべてまたはその一部)でも行うことができ、これによって、様々な点(図14においてz(x,y)として示す(ここでzは深さの値である))において創傷の深さの3Dマップを構築することができる。視差が大きくなるほど、前述のような深さの計算用のアルゴリズムによる計算の負荷が高くなるが、いくつかの実施形態では、視差が大きいほど深さの計算が向上する場合がある。
次に、これらの非線形マッピングモジュール1310A,1310Bの出力は、線形方程式モジュール1320にも提供され、線形方程式モジュール1320は、検知された数値を、スペクトルアンミキシング用の1組の線形方程式として処理することができる。一実装では、少なくともセンサの量子効率およびフィルタの透過率値の関数としてムーア・ペンローズの擬似逆行式を使用することにより、実際のスペクトルの数値(たとえば、各像点(x,y)に入射する特定の波長の光の強さ)を算出することができる。スペクトルアンミキシングは、臨床診断やその他の生物学的応用などの、高い正確度が要求される実装において使用することができる。スペクトルアンミキシングを適用することにより、光量子束およびSN比の推定値を提供することもできる。
ワークフロー1300は、視差を補正したスペクトルチャネル画像およびスペクトルアンミキシングに基づいて、たとえば図に示したF(x,y,λ)という形式(ここでFは、特定波長または周波数帯λにおける特定の結像位置(x,y)における光の強さを示す)などでスペクトルデータキューブ1325を生成することができる。
図15は、たとえば、図3A~10Bおよび図12に示したマルチスペクトル・マルチアパーチャイメージングシステムを使用し、かつ/または図13に示したプロセスを使用して撮影された画像データなどの画像データを処理するための視差および視差補正を図示する。第1の組の画像1410は、4種のセンサ領域で撮影した同じ物体上の同じ物理的位置の画像データを示す。図に示すように、イメージセンサ領域のフォトダイオードグリッドの(x,y)座標フレームに基づくと、これらの生の画像において、物体位置は同じ位置にはない。第2の組の画像1420は、視差補正後に同じ物体位置にあることを示しており、これらの画像では、位置合わせした画像の座標フレームにおいて同じ(x,y)位置になっている。このような位置合わせは、互いに完全には重複していない画像の端部領域から特定のデータをトリミングすることを含んでいてもよいことは十分に理解できるであろう。
図16は、たとえば、図3A~10Bおよび図12に示したマルチスペクトル・マルチアパーチャイメージングシステムを使用して、かつ/または図13に示したプロセスを使用して撮影した後、図14および図15に従って処理された画像データなどのマルチスペクトル画像データを画素単位で分類するためのワークフロー1500を図示する。
ブロック1510において、マルチスペクトル・マルチアパーチャイメージングシステム1513により、物体1511上の物理的な点1512を示す画像データを撮影することができる。この一例において、物体1511は、創傷を有する患者の組織を含む。創傷は、熱傷、糖尿病性潰瘍(たとえば、糖尿病性足部潰瘍)、非糖尿病性潰瘍(たとえば褥瘡または治癒の遅い創傷)、慢性潰瘍、術後切創、(切断術の前またはその後の)切断部位、がん性病変、または損傷を受けた組織を包含しうる。PPG情報が含まれている場合、本明細書で開示するイメージングシステムは、組織灌流;心血管の健康状態;潰瘍などの損傷;末梢動脈疾患;および呼吸器の健康状態を含む、組織の血流の変化および脈拍数の変化を伴う病態を評価する方法を提供する。
ブロック1520において、マルチスペクトル・マルチアパーチャイメージングシステム1513により撮影されたデータを処理して、多数の異なる波長1523を有するマルチスペクトルデータキューブ1525を生成することができ、このマルチスペクトルデータキューブ1525は、異なる時間に対応する同じ波長の多数の異なる画像(PPGデータ1522)を含んでいてもよい。たとえば、画像処理プロセッサ1120は、データキューブ作成モジュール1140を介して、ワークフロー1300に従ってマルチスペクトルデータキューブ1525を生成するように構成することができる。いくつかの実装は、前述したように、空間次元の様々な点における深さの値と関連していてもよい。
ブロック1530において、マルチスペクトルデータキューブ1525は、機械学習モデル1532への入力データ1525として分析して、画像化された組織の分類マッピング1535を生成することができる。分類マッピングにより、画像データ中の各画素(各画素は、位置合わせ後の、画像化された物体1511上の特定の点を示す)を特定の組織分類または特定の治癒可能性スコアに割り当てることができる。分類された画像の出力において、視覚的に異なる色またはパターンを使用することにより、様々な分類または様々なスコアを示すことができる。したがって、物体1511の画像を数多く撮影したとしても、その出力は、画素単位で分類された複数の視覚表示を重ねた物体の単一の画像(たとえば、通常のRGB画像)になりうる。
いくつかの実装では、機械学習モデル1532は人工ニューラルネットワークであってもよい。人工ニューラルネットワークは、生物学的神経回路網に着想を得たものであるが、コンピュータ装置による実装を目的として改良された計算実体であるという意味で人工物である。入力と出力の間の依存性を見出すことは容易ではないが、人工ニューラルネットワークを使用することにより、入力と出力の間の複雑な関係をモデル化したり、データ中において一定のパターンを見つけたりすることができる。ニューラルネットワークは、通常、入力層、1層以上の中間層(「隠れ層」)および出力層を含み、各層には数多くのノードが含まれる。ノードの数は各層によって異なっていてもよい。2層以上の隠れ層を含むニューラルネットワークは、「深層」であると見なされる。各層のノードは、次の層のノードのすべてまたはその一部と結合しており、これらの結合の重みは、通常、訓練プロセス中のデータから学習され、たとえば、ニューラルネットワークのパラメータを調整して、ラベルを付けた訓練データ中の対応する入力から予想される出力を生成するバックプロパゲーション(誤差逆伝播法)を使用して行われる。したがって、人工ニューラルネットワークは、人工ニューラルネットワークを通して流れる情報に基づいて、訓練中にその構造(たとえば結合構成および/または重み)を変更するように構成された適応システムであり、隠れ層の重みは、データ中の意味のあるパターンをコードするものであると見なすことができる。
全結合ニューラルネットワークとは、入力層の各ノードが次の層(第1の隠れ層)の各ノードに結合され、第1の隠れ層の各ノードが次の隠れ層の各ノードに結合され、このような結合が、最終的な隠れ層の各ノードが出力層の各ノードに結合されるまで継続するニューラルネットワークである。
CNNは、人工ニューラルネットワークの1種であり、前述した人工ニューラルネットワークと同様に、複数のノードから構成され、学習可能な重みを有する。前述の人工ニューラルネットワークとは異なる点として、CNNの各層は、各ビデオフレームの2×2アレイの画素値(たとえば幅と高さ)と、シーケンス中のビデオフレームの数(たとえば深さ)に対応した、幅と高さと深さからなる三次元に配置されたノードを有することができる。各層のノードは、幅と高さからなるその前の層の小領域(受容野と呼ぶ)に局所結合されていてもよい。隠れ層の重みは、受容野に適用される畳み込みフィルタの形態を取ることができる。いくつかの実施形態において、畳み込みフィルタは二次元であってもよく、したがって、入力ボリュームの各フレームまたは指定されたフレームサブセットに対して同じフィルタを使用して畳み込み(または画像の畳み込み変換)を繰り返すことができる。別の実施形態では、畳み込みフィルタは三次元であってもよく、したがって、入力ボリュームの複数のノードの深さ全体に及ぶものであってもよい。CNNの各畳み込み層の各ノードは重みを共有して、入力ボリュームの幅と高さの全体に対して(たとえば全フレームに対して)所定の層の畳み込みフィルタを複製することにより、訓練可能な重みの総数を減らし、訓練データ外のデータセットに対するCNNの適用性を向上させることができる。特定の層の数値をプールして、次の層の計算の回数を減らしてもよく(たとえば、特定の画素を示す数値を次の層へと伝達してもよく、その他の数値は破棄する)、さらに、CNNの深さが進むに従って、廃棄された数値にプールマスクを再導入し、元の大きさに対してデータポイントの数を返してもよい。一部が全結合されていてもよい多数の層を積み重ねて、CNN構造を形成することができる。
訓練中、人工ニューラルネットワークをその訓練データ中のペアデータに適用して、パラメータを変更し、入力が提供された場合に特定のペアの出力を予測することができる。たとえば、訓練データは、マルチスペクトルデータキューブ(入力)と、たとえば、特定の臨床状態に対応する創傷領域を臨床医が指定することなどによりラベル付けされた分類マッピング(予想される出力)とを含むことができ、このラベル付けされた分類マッピングには、最初に創傷のイメージングを行ってからしばらく経過した後に実際の治癒が確認されてから付与された治癒ラベル(1)または未治癒ラベル(0)が含まれていてもよい。機械学習モデル1532の別の実装では、その他の種類の予測を行うように訓練することもでき、たとえば、特定の期間中に創傷が治癒して、創傷領域の割合が特定のパーセンテージまで減少する(たとえば、30日以内に創傷領域が少なくとも50%減少する)可能性を予測したり、止血、炎症、増殖、リモデリング、正常な皮膚などに分類される創傷の状態を予測するように訓練することもできる。いくつかの実装では、分類の正確度をさらに向上させるため、患者の測定値を入力データにさらに組み込んでもよく、あるいは、機械学習モデル1532に別の実例を訓練させるため、患者の測定値に基づいて訓練データをセグメント化し、同じ測定値を有する別の患者に使用してもよい。患者の測定値には、患者の特性またはその患者の健康状態の特性を説明する文字情報または病歴またはそれらの態様が含まれていてもよく、患者の特性またはその患者の健康状態の特性としては、たとえば、創傷、病変または潰瘍の面積;患者のBMI;患者の糖尿病の状態;患者における末梢血管疾患または慢性炎症の存在;患者が現在罹患しているその他の創傷の数または患者がこれまでに罹患した創傷の数;免疫抑制剤(たとえば化学療法)または創傷治癒率に正もしくは負の影響を及ぼすその他の薬剤が患者に投与されているか否か、または最近投与されたか否か;HbA1c;慢性腎不全ステージIV;1型糖尿病または2型糖尿病のいずれであるか;慢性貧血;喘息;薬剤の使用;喫煙の有無;糖尿病性神経障害;深部静脈血栓症;過去の心筋梗塞症;一過性脳乏血発作;および睡眠時無呼吸;ならびにこれらの任意の組み合わせが挙げられる。これらの測定値は、適切な処理によりベクトル表示に変換することができ、たとえば、単語をベクトル表現化した埋め込み、患者が特定の測定値を満たすか否か(たとえば、1型糖尿病に罹患しているか否か)を示す2進値のベクトル、または患者における各測定値の程度を示す数値などによって、ベクトル表示に変換することができる。
ブロック1540では、ユーザに対して分類マッピング1535を出力することができる。この一例では、分類マッピング1535では、第1の色1541を使用して、第1の状態に従って分類された画素が示され、第2の色1542を使用して、第2の状態に従って分類された画素が示される。分類とその結果得られた分類マッピング1535は、たとえば、物体認識、背景色の識別および/または深さの値に基づいたバックグラウンド画素が除外されていてもよい。図に示したように、マルチスペクトル・マルチアパーチャイメージングシステム1513のいくつかの実装では、組織部位上に分類マッピング1535を投影することができる。このような態様は、推奨される切除マージンおよび/または切除する深さの視覚表示が分類マッピングに含まれる場合に、特に有益となる場合がある。
このような方法およびシステムによって、熱傷の切除、切断レベル、病変除去および創傷のトリアージにおける判断などの、皮膚創傷の管理プロセスにおいて臨床医および外科医を補助してもよい。本明細書に記載の実施形態は、褥瘡性潰瘍、充血、四肢の悪化、レイノー現象、強皮症、慢性創傷、擦過創、裂創、出血、破裂損傷、穿刺創、穿通創、皮膚がん(基底細胞癌、扁平上皮癌、悪性黒色腫もしくは光線角化症)、または組織の特性および性質が正常状態とは異なる何らか種類の組織変化の重症度の特定および/または分類するために使用することができる。また、本明細書に記載の装置を使用して、正常組織のモニター;(たとえば、デブリードマンを行う際のマージンを決定するための、より速く、より洗練されたアプローチが可能となるような)創傷治療法の補助および改善;ならびに(特に治療が適用された後の)創傷または疾患からの回復の進行の評価を行うことができる。本明細書に記載のいくつかの実施形態では、損傷を受けた組織に隣接した正常組織の識別、切除の際のマージンおよび/もしくは深さの決定、左室補助人工心臓などの人工装具を移植した後の回復過程のモニタリング、組織移植片または再生細胞移植片の生存能力の評価、または(特に再建術を行った後の)術後回復のモニタリングが可能な装置を提供する。さらに、本明細書に記載の実施形態を利用して、特に、ステロイド、肝細胞成長因子、線維芽細胞成長因子、抗生物質、再生細胞(幹細胞、内皮細胞および/または内皮前駆細胞を含む、単離または濃縮された細胞集団など)などの治療剤の投与後の、創傷の変化または創傷後の正常組織の発生を評価することができる。
分散コンピューティング環境の一例の概要
図17は、図3A~10Bおよび図12に示したマルチスペクトル・マルチアパーチャイメージングシステムのいずれであってもよいマルチスペクトル・マルチアパーチャイメージングシステム1605を含む、分散コンピューティングシステム1600の一例の略ブロック図を示す。図に示すように、データキューブ分析サーバ1615は、恐らくはサーバクラスタまたはサーバファームとして配置される1つ以上のコンピュータを含んでいてもよい。これらのコンピュータを構成するメモリおよびプロセッサは、1つのコンピュータ内に配置されていてもよく、(互いに離れて設置されたコンピュータを含む)数多くのコンピュータに分散されていてもよい。
マルチスペクトル・マルチアパーチャイメージングシステム1605は、ネットワーク1610を介して、ユーザ装置1620およびデータキューブ分析サーバ1615と通信するためのネットワーキングハードウェア(たとえば、無線のインターネット、衛星通信、ブルートゥースまたはその他の通信機)を含むことができる。たとえば、いくつかの実装において、マルチスペクトル・マルチアパーチャイメージングシステム1605のプロセッサは、画像撮影を制御し、その後、生データをデータキューブ分析サーバ1615に送信するように構成されていてもよい。別の実装でのマルチスペクトル・マルチアパーチャイメージングシステム1605のプロセッサは、画像撮影を制御し、スペクトルアンミキシングおよび視差補正を行ってマルチスペクトルデータキューブを生成し、その後、マルチスペクトルデータキューブをデータキューブ分析サーバ1615に送信するように構成されていてもよい。いくつかの実装では、マルチスペクトル・マルチアパーチャイメージングシステム1605上においてすべての処理および分析を局所的に行うことができ、集計分析の実施および/または機械学習モデルの訓練もしくは再訓練での使用を目的として、マルチスペクトルデータキューブおよび得られた分析結果をデータキューブ分析サーバ1615に送信してもよい。したがって、データキューブ分析サーバ1615は、マルチスペクトル・マルチアパーチャイメージングシステム1605に最新の機械学習モデルを提供してもよい。マルチスペクトルデータキューブを分析した最終結果を得るための処理負荷は、マルチアパーチャイメージングシステム1605の処理能力に応じて、マルチアパーチャイメージングシステム1605とデータキューブ分析サーバ1615の間で様々な方法により分担してもよい。
ネットワーク1610は、イントラネット、インターネット、セルラーネットワーク、ローカルエリアネットワーク、同様のその他のネットワーク、これらの組み合わせなどの、適切なネットワークを含むことができる。ユーザ装置1620には、ネットワークを装備したあらゆるコンピュータ装置が含まれ、たとえば、デスクトップコンピュータ、ノートパソコン、スマートフォン、タブレット、電子書籍リーダー、テレビゲーム機などが挙げられる。たとえば、組織分類が必要となった際に、マルチアパーチャイメージングシステム1605およびデータキューブ分析サーバ1615により測定された結果(たとえば、分類された画像)を、患者もしくは医師の指定のユーザ装置、患者の電子的な医療記録を格納する病院情報システム、および/または(たとえば、米国疾病対策予防センターなどが保有する)集中型公衆衛生データベースに送信してもよい。
実装の結果の一例
背景:
負傷した戦闘員やその医療従事者にとって、熱傷による罹患および死亡は課題となっている。過去の戦闘死傷者における熱傷の受傷率は5~20%であり、これらの死傷者の約20%は、米国陸軍外科研究所(ISR)熱傷センターなどでの複雑な熱傷外科手術が必要であった。熱傷外科手術は専門的な訓練を必要とすることから、米国陸軍病院のスタッフではなく、ISRのスタッフにより熱傷外科手術が行われる。熱傷専門医の数が限られていることから、熱傷を負った兵士に医療を提供するための兵站業務が非常に複雑となる。したがって、術前および術中に熱傷の深さを検出する新規な客観的方法を利用することによって、戦闘継続中の熱傷患者に対する医療の提供への従事が可能な、(ISR以外の組織に所属する人員を含む)医療スタッフの数をより多く確保することができる。医療従事者の人員確保が増強されることによって、熱傷を負った戦闘員への医療的役割の進展に向けて、より複雑な熱傷治療を拡充させることができる。
このようなニーズへの対応を開始するため、術前における熱傷の治癒可能性の判定を補助することを目的として、マルチスペクトルイメージング(MSI)および人工知能(AI)アルゴリズムを使用した新規なカート用イメージング装置を開発した。この装置は、短時間で(たとえば、6秒以内、5秒以内、4秒以内、3秒以内、2秒以内または1秒以内に)広範囲な組織(たとえば、5.9×7.9平方インチ)の画像を取得することができ、イメージング造影剤の注射を必要としない。一般市民を対象としたこの研究では、この装置による熱傷の治癒可能性の判定の正確度は、熱傷の専門医師による臨床判断よりも優れていること(たとえば70~80%に達する)が示されている。
方法:
様々な重症度の熱傷を負った一般市民の患者において、熱傷受傷の72時間以内に画像化を実施し、熱傷受傷後の最長で7日間にわたりいくつかの時点でも画像化を実施した。各画像における正確な熱傷の重症度は、3週間にわたる治癒評価またはパンチ生検を使用して判定した。この装置での、I度熱傷、II度熱傷およびIII度熱傷における治癒過程の熱傷組織と非治癒熱傷組織の特定および判別の正確度は、画素単位で分析した。
結果:
データは、38人の一般市民の患者から収集し、合計で58個の熱傷および393個の画像を取得した。AIアルゴリズムにより、非治癒熱傷組織の予測において87.5%の感度と90.7%の特異度が達成された。
結論:
本発明の新規な装置およびそのAIアルゴリズムによる熱傷の治癒可能性の判定の正確度は、熱傷の専門医の臨床判断よりも優れていることが示された。将来的には、携帯を可能にするための本装置の再設計と術中の状況下での使用の評価に焦点を当てて検討を行う予定である。携帯を可能にするための設計変更としては、本装置のサイズを携帯型システムのサイズまで小さくすること、視野を広くすること、1回のスナップショットの取得時間を短くすること、およびブタモデルを使用して、術中の状況下での使用において本装置を評価することが挙げられる。これらの開発は、基本的な画像検査において同等の機能を示すベンチトップ型のマルチスペクトルイメージング(MSI)サブシステムを使用して実施した。
画像の位置合わせ用の追加の光源
様々な実施形態において、本明細書で開示した実施形態のいずれかと1個以上の追加の光源とを併用して、画像の位置合わせの正確度を向上させてもよい。図21は、プロジェクター2105を含むマルチアパーチャスペクトル撮像装置2100の実施形態の一例を示す。いくつかの実施形態において、プロジェクター2105またはその他の適切な光源は、たとえば、図12に関して前述した複数の光源1165のうちの1個であってもよい。位置合わせ用のプロジェクター2105などの追加の光源を含む実施形態において、前記方法は、追加の露光をさらに含んでいてもよい。プロジェクター2105などの追加の光源は、撮像装置2100のすべてのカメラを通して個別にまたは累積的に可視化することが可能な、単一のスペクトル帯域、複数のスペクトル帯域または広帯域の、1つ以上の点、縞模様、グリッド、ランダムな斑点またはその他の適切な空間パターンを、撮像装置2100の視野内に投影することができる。たとえば、プロジェクター2105は、前述の共通帯域を利用したアプローチに基づいて計算された画像の位置合わせの正確度の確認に使用することが可能な、共有チャネルもしくは共通チャネルの光、広帯域照明または累積的に可視可能な照明を投影してもよい。本明細書において、「累積的に可視可能な照明」は、選択された複数の波長であり、マルチスペクトルイメージングシステムの各イメージセンサによりそのパターンが変換される波長を指す。たとえば、累積的に可視可能な照明として、複数の波長であって、すべてのチャネルに共通の波長がこの複数の波長に含まれていなくても、この複数の波長のうちの少なくとも1つが各チャネルにより変換される波長を挙げることができる。いくつかの実施形態において、プロジェクター2105により投影されるパターンの種類は、そのパターンが画像化されるアパーチャの数に基づいて選択してもよい。たとえば、1個のアパーチャのみでしかパターンを見ることができない場合、このパターンは比較的高密度であることが好ましい場合があり(たとえば、約1~10画素、約20画素、約50画素未満、約100画素未満などの、比較的狭い自己相関を有していてもよい)、一方、複数個のアパーチャによりパターンが画像化される場合は、より低密度のパターンまたは自己相関がそれほど狭くないパターンが有用である場合がある。いくつかの実施形態において、投影された空間パターンとともに撮影される追加の露光は、投影された空間パターンを含めずに露光が撮影される実施形態よりも位置合わせの正確度が向上されることを意図として、視差の計算に含められる。いくつかの実施形態において、追加の光源は、すべてのカメラを通して個別にまたは累積的に可視化することが可能な、単一のスペクトル帯域、複数のスペクトル帯域または広帯域(たとえば、共有チャネルもしくは共通チャネル、または広帯域照明など)の縞模様を撮像装置の視野内に投影し、これを利用することによって、縞模様の位相に基づいて画像の位置合わせを向上させることができる。いくつかの実施形態において、追加の光源は、すべてのカメラを通して個別にまたは累積的に可視化することが可能な、単一のスペクトル帯域、複数のスペクトル帯域または広帯域(たとえば、共有チャネルもしくは共通チャネル、または広帯域照明など)のドット、グリッドおよび/または斑点からなる複数の特有の空間配置を撮像装置の視野内に投影し、これを利用することによって、画像の位置合わせを向上させることができる。いくつかの実施形態において、前記方法は、1個のアパーチャまたは複数個のアパーチャを備えた追加のセンサをさらに含み、この追加のセンサによって、視野内の単一または複数の物体の形状を検出することができる。たとえば、この追加のセンサは、LIDAR技術、ライトフィールド技術または超音波技術を使用して、前述の共通帯域を利用したアプローチによる画像の位置合わせの正確度をさらに向上させてもよい。この追加センサは、ライトフィールド情報を検知可能な1個のアパーチャまたはマルチアパーチャセンサであってもよく、超音波やパルスレーザなどのその他のシグナルを検知可能であってもよい。
用語
本明細書に記載の方法およびタスクはすべてコンピュータシステムにより実行されてもよく、完全に自動化されていてもよい。場合によっては、このコンピュータシステムは、ネットワーク上で通信し相互運用して、本明細書で述べる機能を実行する別々の複数のコンピュータまたはコンピュータ装置(たとえば、物理サーバ、ワークステーション、ストレージアレイ、クラウドコンピューティングリソースなど)を含む。このようなコンピュータ装置は、いずれも、通常、メモリまたはその他の非一時的なコンピュータ可読記憶媒体もしくは非一時的なコンピュータ可読記憶装置(たとえば、ソリッドステート記憶装置、ディスクドライブなど)に格納されたプログラム命令またはプログラムモジュールを実行するプロセッサ(または複数のプロセッサ)を含む。本明細書で開示した様々な機能は、そのようなプログラム命令で具体化してもよく、コンピュータシステムの特定用途向け回路(たとえば、ASICやFPGA)の形態で実装してもよい。コンピュータシステムが複数のコンピュータ装置を含む場合、これらの装置は、同じ場所に配置してもよいし、別の場所に配置してもよい。本明細書で開示した方法およびタスクによる結果は、ソリッドステートメモリチップや磁気ディスクなどの物理的記憶装置を様々な形態に変換することによって、永続的に保存してもよい。いくつかの実施形態において、コンピュータシステムは、複数の別個の企業体またはその他のユーザにより処理リソースが共有されるクラウドベースのコンピューティングシステムであってもよい。
本明細書で開示したプロセスは、ユーザまたはシステムアドミニストレータによって要求が開始された際に、所定のスケジュールや動的に決定したスケジュールなどのイベントに応答して開始してもよく、その他のいくつかのイベントに応答して開始してもよい。前記プロセスが開始されたら、1つ以上の非一時的なコンピュータ可読媒体(たとえば、ハードドライブ、フラッシュメモリ、取り外し可能な媒体など)に格納された実行可能プログラム命令群を、サーバまたはその他のコンピュータ装置のメモリ(たとえばRAM)にロードしてもよい。次に、実行可能命令は、コンピュータ装置が備えるハードウェアベースのコンピュータプロセッサにより、実行してもよい。いくつかの実施形態において、本明細書で開示したプロセスまたはその一部は、直列または並列に、複数のコンピュータ装置および/または複数のプロセッサ上に実装してもよい。
実施形態に応じて、本明細書に記載のプロセスまたはアルゴリズムのいずれかによる特定の行為、イベントまたは機能は、別の順序で行うことができ、互いに付加することができ、互いに合体することができ、すべてを省略することができる(たとえば、本明細書に記載の作動およびイベントのすべてが、本明細書に記載のアルゴリズムの実行に必要とは限らない)。さらに、特定の実施形態において、複数の作動またはイベントは、たとえば、マルチスレッド処理、割り込み処理、もしくは複数のプロセッサもしくは複数のプロセッサコアを使用して、またはその他の並列構造上で、順次ではなく同時に行うことができる。
本明細書で開示した実施形態に関連して記載した、例示的な様々な論理ブロック、モジュール、ルーチン、およびアルゴリズムの工程は、電子機器(たとえばASICまたはFPGA装置)、コンピュータハードウェア上で作動するコンピュータソフトウェア、またはこれらの組み合わせとして実装することができる。さらに、本明細書で開示した実施形態に関連して記載した、例示的な様々な論理ブロックおよびモジュールは、プロセッサ装置、デジタルシグナルプロセッサ(「DSP」)、特定用途向け集積回路(「ASIC」)、フィールドプログラマブルゲートアレイ(「FPGA」)もしくはその他のプログラマブルロジックデバイス、ディスクリートゲートロジックもしくはトランスファーロジック、ディスクリートハードウェアコンポーネント、または本明細書に記載の機能を実行するように設計された、これらの部品の組み合わせなどの機器により実装することができる。プロセッサ装置はマイクロプロセッサであってもよいが、別の態様において、プロセッサ装置は、制御装置、マイクロコントローラまたは状態機械、これらの組み合わせなどであってもよい。プロセッサ装置は、コンピュータ実行可能命令を処理するように構成された電気回路を含んでいてもよい。別の一実施形態において、プロセッサ装置は、FPGAを含むか、またはコンピュータ実行可能命令を処理することなく論理演算を行うその他のプログラマブルデバイスを含む。また、プロセッサ装置は、コンピュータ装置の組み合わせとして実装することができ、たとえば、DSPとマイクロプロセッサの組み合わせ、複数のマイクロプロセッサの組み合わせ、DSPコアと1個以上のマイクロプロセッサの組み合わせ、またはこのようなその他の構成として実装することができる。本明細書では、主にデジタル技術に関して説明をしてきたが、プロセッサ装置は、主にアナログ部品を含んでいてもよい。たとえば、本明細書に記載のレンダリング技術のすべてまたはそのうちの一部は、アナログ回路を使用して、またはアナログ回路とデジタル回路を併用して実装してもよい。コンピューティング環境は、何らかの種類のコンピュータシステムを含むことができ、このコンピュータシステムとして、いくつか例を挙げれば、マイクロプロセッサに基づくコンピュータシステム、メインフレームコンピュータ、デジタルシグナルプロセッサ、携帯型コンピュータ装置、デバイスコントローラ、電気器具内のコンピュータエンジンなどがあるが、これらに限定されない。
本明細書で開示した実施形態に関して述べた方法、プロセス、ルーチンまたはアルゴリズムの構成要素は、ハードウェアで直接具現化することができ、プロセッサ装置により実行されるソフトウェアモジュールで直接具現化することもでき、またはこれらのハードウェアとソフトウェアモジュールの組み合わせで直接具現化することもできる。ソフトウェアモジュールは、RAMメモリ、フラッシュメモリ、ROMメモリ、EPROMメモリ、EEPROMメモリ、レジスタ、ハードディスク、リムーバブルディスク、CD-ROM、またはその他の形態の非一時的なコンピュータ可読記憶媒体に含まれていてもよい。典型的な記憶媒体は、プロセッサ装置がこのような記憶媒体から情報を読み取ったり、このような記憶媒体に情報を書き込んだりできるように、プロセッサ装置に接続することができる。この態様において、記憶媒体はプロセッサ装置にとって不可欠でありうる。このプロセッサ装置および記憶媒体は、ASICに含まれていてもよい。このASICは、ユーザ端末に含まれていてもよい。この態様において、このプロセッサ装置および記憶媒体は、ユーザ端末のディスクリート部品として存在していてもよい。
本明細書で使用されている、「できる」、「場合がある」、「可能性がある」、「してもよい」、「たとえば」などの、条件を示す用語は、別段の記載がない限り、あるいはこれらの用語が使用されている文脈上で別の意味に解釈されない限り、通常、特定の実施形態が特定の特徴、構成要素または工程を含み、別の実施形態には、これらの特徴、構成要素または工程が含まれないことを意味する。したがって、条件を示すこのような用語は、通常、特定の実施形態において特定の特徴、構成要素または工程が含まれるのか、あるいは実施されるのかにかかわらず、その他の入力やプロンプティングの存在下または非存在下において、1つ以上の実施形態が、これらの特徴、構成要素または工程を何らかの形で必要とすることや、1つ以上の実施形態が、決定のロジックを必然的に含むことを含意しない。また、「含む」、「備える」、「有する」などの用語は同義語であり、オープンエンド形式で包括的な意味で使用され、追加の構成要素、特徴、行為、作動などを除外するものではない。さらに、「または」という用語も、(排他的な意味ではなく)包括的な意味で使用され、たとえば、「または」という用語が、構成要素を列記するために使用された場合、列記された構成要素のうちの1つ、その一部またはすべてを意味する。
別段の記載がない限り、「X、YまたはZの少なくとも1つ」などの選言的な用語は、本明細書中において、通常、特定の項目や用語などが、X、YもしくはZ、またはこれらの任意の組み合わせ(たとえば、X、YまたはZ)であってもよいことを示すために使用されると理解される。したがって、このような選言的な用語は、通常、特定の実施形態において、Xの少なくとも1つと、Yの少なくとも1つと、Zの少なくとも1つのそれぞれの存在が必要であることを意味することを意図せず、そのように意図されるべきではない。
前述の詳細な説明において、様々な実施形態として適用される新規な特徴を示し、説明し、指摘してきたが、前述の装置またはアルゴリズムの形態およびその細部は、本開示の範囲から逸脱することなく、様々に省略、置換および変更することができることを理解できるであろう。本明細書に示した特徴のいくつかは、その他の特徴とは別々に使用または実施できることから、本明細書に記載の特定の実施形態は、本明細書に示された特徴および利点のすべてを提供するわけではない形態として具体化することができることは容易に理解できるであろう。請求項の記載と等価の意味および範囲内のあらゆる変更は、これらの請求項の範囲内に包含される。

Claims (73)

  1. マルチスペクトル画像システムであって、
    少なくとも1つのイメージセンサ;
    前記少なくとも1つのイメージセンサの第1のセンサ領域へと入射光を通過させるように配置された第1のアパーチャ;
    前記少なくとも1つのイメージセンサの第2のセンサ領域へと入射光を通過させるように配置された第2のアパーチャ;
    第1のアパーチャの上方に配置され、少なくとも共通周波数帯の光および第1の特有周波数帯の光を通過させるように構成された第1のマルチバンドパスフィルタ;
    第2のアパーチャの上方に配置され、少なくとも前記共通周波数帯の光および第2の特有周波数帯の光を通過させるように構成された第2のマルチバンドパスフィルタ;
    マルチスペクトルデータキューブを生成させる命令を格納するメモリ;ならびに
    前記命令に従って、少なくとも、第1のセンサ領域および第2のセンサ領域からシグナルを受信し、受信したこのシグナルを処理して、前記マルチスペクトルデータキューブを生成するように構成された少なくとも1つのプロセッサ
    を含み、
    前記マルチスペクトルデータキューブが、前記共通周波数帯に対応する第1のスペクトルチャネル、第1の特有周波数帯に対応する第2のスペクトルチャネル、および第2の特有周波数帯に対応する第3のスペクトルチャネルを含み;
    前記シグナルの処理が、
    少なくとも、前記共通周波数帯に対応する第1のスペクトルチャネルからの情報を利用して、第1のセンサ領域により撮影された第1の画像データと第2のセンサ領域により撮影された第2の画像データの間の視差を概算すること、ならびに
    前記視差を利用して、第1のスペクトルチャネルの画像データ、第2のスペクトルチャネルの画像データおよび第3のスペクトルチャネルの画像データの位置調整を行い、マルチスペクトルデータキューブを生成することを含む、
    マルチスペクトル画像システム。
  2. 第1のマルチバンドパスフィルタが、前記マルチスペクトル画像システムのレンズアセンブリの前面に配置されるか、前記マルチスペクトル画像システムの単レンズの前面に配置されるか、前記マルチスペクトル画像システムのレンズアセンブリの後ろに配置されるか、前記マルチスペクトル画像システムの単レンズの後ろに配置されるか、前記マルチスペクトル画像システムのレンズアセンブリの内部に配置されるか、前記マルチスペクトル画像システムの結合レンズ群の内部に配置されるか、前記マルチスペクトル画像システムの単レンズの表面上に直接配置されるか、または前記マルチスペクトル画像システムのレンズアセンブリの素子の表面上に直接配置される、請求項1に記載のマルチスペクトル画像システム。
  3. 前記少なくとも1つのプロセッサが、前記命令に従って、少なくともスペクトルアンミキシングを行うことによって、第1のスペクトルチャネルの画像データ、第2のスペクトルチャネルの画像データおよび第3のスペクトルチャネルの画像データを生成するように構成されている、請求項1または2に記載のマルチスペクトル画像システム。
  4. 第1のセンサ領域および第2のセンサ領域を覆うカラーフィルタアレイをさらに含み、該カラーフィルタアレイが、カラーフィルタの繰り返しパターンを含み;
    前記繰り返しパターンの第1のカラーフィルタが、前記共通周波数帯を含む第1の重み付け光波長分布を通過させるように構成されており;
    前記繰り返しパターンの第2のカラーフィルタが、第1の特有周波数帯を含む第2の重み付け光波長分布を通過させるように構成されており;
    前記繰り返しパターンの第3のカラーフィルタが、第2の特有周波数帯を含む第3の重み付け光波長分布を通過させるように構成されている、
    請求項1~3のいずれか1項に記載のマルチスペクトル画像システム。
  5. 前記共通周波数帯が、紫色光、青紫色光、青色光および青緑色光の少なくとも1つを含む、請求項1~4のいずれか1項に記載のマルチスペクトル画像システム。
  6. 前記共通周波数帯が、380nmよりも長く450nmよりも短い波長を有する光を含む、請求項1~4のいずれか1項に記載のマルチスペクトル画像システム。
  7. 前記共通周波数帯が緑色光を含む、請求項1~4のいずれか1項に記載のマルチスペクトル画像システム。
  8. 前記共通周波数帯が赤色光を含む、請求項1~4のいずれか1項に記載のマルチスペクトル画像システム。
  9. 第1のマルチバンドパスフィルタおよび第2のマルチバンドパスフィルタのそれぞれが、緑色光、青緑色光および青色光の少なくとも1つを含む第2の共通周波数帯を通過させるように構成されている、請求項5~8のいずれか1項に記載のマルチスペクトル画像システム。
  10. 第1のマルチバンドパスフィルタおよび第2のマルチバンドパスフィルタのそれぞれが、赤色光を含む第2の共通周波数帯を通過させるように構成されている、請求項5~8のいずれか1項に記載のマルチスペクトル画像システム。
  11. 第1のマルチバンドパスフィルタおよび第2のマルチバンドパスフィルタのそれぞれが、前記カラーフィルタアレイのパターンを構成するカラーフィルタのうち、第1の共通周波数帯を遮断しない少なくとも1枚のカラーフィルタにより遮断されない第2の共通周波数帯を通過させるように構成されている、請求項5~8のいずれか1項に記載のマルチスペクトル画像システム。
  12. 前記少なくとも1つのプロセッサが、前記命令に従って、少なくともスペクトルアンミキシングを行うことによって、撮影画像データの、前記共通周波数帯に対応する部分と第2の共通周波数帯に対応する部分を分離するように構成されている、請求項9または10に記載のマルチスペクトル画像システム。
  13. 前記少なくとも1つのプロセッサが、前記命令に従って、前記マルチスペクトルデータキューブを使用して画像化組織を少なくとも分析するように構成されている、請求項1~12のいずれか1項に記載のマルチスペクトル画像システム。
  14. 前記少なくとも1つのイメージセンサの第3のセンサ領域へと入射光を通過させるように配置された第3のアパーチャ;
    前記少なくとも1つのイメージセンサの第4のセンサ領域へと入射光を通過させるように配置された第4のアパーチャ;
    第3のアパーチャの上方に配置され、少なくとも共通周波数帯の光および第3の特有周波数帯の光を通過させるように構成された第3のマルチバンドパスフィルタ;ならびに
    第4のアパーチャの上方に配置され、少なくとも前記共通周波数帯の光および第4の特有周波数帯の光を通過させるように構成された第4のマルチバンドパスフィルタ
    をさらに含む、請求項1~13のいずれか1項に記載のマルチスペクトル画像システム。
  15. 前記少なくとも1つのイメージセンサの第5のセンサ領域へと入射光を通過させるように配置された第5のアパーチャ;および
    第3のアパーチャの上方に配置された第5のマルチバンドパスフィルタ
    をさらに含み、
    該第5のマルチバンドパスフィルタが、第1のマルチバンドパスフィルタ、第2のマルチバンドパスフィルタ、第3のマルチバンドパスフィルタおよび第4のマルチバンドパスフィルタの少なくとも1つを通過する少なくとも1つの周波数帯の光を通過させるように構成されている、
    請求項14に記載のマルチスペクトル画像システム。
  16. 前記少なくとも1つのイメージセンサの最大総数25個のセンサ領域へと入射光を通過させるように配置された最大総数25個のアパーチャ;および
    前記最大総数25個のアパーチャのそれぞれの上方に配置された最大総数25枚のマルチバンドパスフィルタ
    を含み、
    各マルチバンドパスフィルタが、少なくとも共通周波数帯の光と、前記画像システム内の他のマルチバンドパスフィルタを通過しない特有周波数帯の光とを通過させるように構成されており;
    通過する周波数帯の総数が、4個を超え、かつ51個未満である、
    請求項1~14のいずれか1項に記載のマルチスペクトル画像システム。
  17. 前記少なくとも1つのイメージセンサの第3のセンサ領域へと入射光を通過させるように配置された第3のアパーチャ;および
    第3のアパーチャの上方に配置された第3のマルチバンドパスフィルタ
    をさらに含み、
    該第3のマルチバンドパスフィルタが、第1のマルチバンドパスフィルタおよび第2のマルチバンドパスフィルタの少なくとも一方を通過する少なくとも1つの周波数帯の光を通過させるように構成されている、
    請求項1~13のいずれか1項に記載のマルチスペクトル画像システム。
  18. 第1のマルチバンドパスフィルタ、第2のマルチバンドパスフィルタ、第3のマルチバンドパスフィルタおよび第4のマルチバンドパスフィルタが、これら全体で7個、8個、9個または10個の異なる周波数帯を通過させ、これらの異なる周波数帯に、共通周波数帯、第1の特有周波数帯、第2の特有周波数帯、第3の特有周波数帯および第4の特有周波数帯が含まれている、請求項14~16のいずれか1項に記載のマルチスペクトル画像システム。
  19. 第1のマルチバンドパスフィルタ、第2のマルチバンドパスフィルタ、第3のマルチバンドパスフィルタおよび第4のマルチバンドパスフィルタのそれぞれが、青色光を含む共通周波数帯、緑色光を含む追加の共通周波数帯、および2個の追加の周波数帯を通過させるように構成されている、請求項18に記載のマルチスペクトル画像システム。
  20. 前記少なくとも1つのプロセッサが、前記命令に従って、前記7個、8個、9個または10個の異なる周波数帯に対応する10個のスペクトルチャネルを含むマルチスペクトルデータキューブを少なくとも生成するように構成されている、請求項18または19に記載のマルチスペクトル画像システム。
  21. 前記7個、8個、9個または10個の異なる周波数帯の光を発するように構成された複数の発光ダイオード(LED)を含む照明基板をさらに含み、個々のLEDが、前記7個、8個、9個または10個の異なる周波数帯のいずれかの光を発する、請求項20に記載のマルチスペクトル画像システム。
  22. 前記少なくとも1つのイメージセンサが、該少なくとも1つのイメージセンサのフォトダイオードの上方に配置されたカラーフィルタアレイを含み、該カラーフィルタアレイが、少なくとも2枚の異なるカラーフィルタを有する、請求項18~21のいずれか1項に記載のマルチスペクトル画像システム。
  23. 前記カラーフィルタアレイが、3~25枚の異なるカラーフィルタを有する、請求項22に記載のマルチスペクトル画像システム。
  24. 前記7個、8個、9個または10個の異なる周波数帯が、前記3枚の異なるカラーフィルタを通過する周波数帯よりも狭い周波数帯の光を含む、請求項22または23に記載のマルチスペクトル画像システム。
  25. 前記共通周波数帯の光、第1の特有周波数帯の光および第2の特有周波数帯の光を発するように構成された少なくとも1つの光源をさらに含む、請求項1~24のいずれか1項に記載のマルチスペクトル画像システム。
  26. 前記少なくとも1つの光源の上方に配置された拡散光学素子をさらに含む、請求項25に記載のマルチスペクトル画像システム。
  27. 少なくとも第1のマルチバンドパスフィルタが湾曲したフィルタである、請求項1~26のいずれか1項に記載のマルチスペクトル画像システム。
  28. 請求項1~27のいずれか1項に記載のマルチスペクトル画像システムを使用してマルチスペクトルデータキューブを生成する方法であって、
    第1のスペクトルチャネル、第2のスペクトルチャネルおよび第3のスペクトルチャネルを含む第1の露光を撮影する工程;
    第1のスペクトルチャネル、第2のスペクトルチャネルおよび第3のスペクトルチャネルを含む第2の露光を撮影する工程;ならびに
    第1の露光のシグナル/ノイズ比を低下させるため、第2の露光の使用に少なくとも部分的に基づいてマルチスペクトルデータキューブを生成する工程
    を含む方法。
  29. 第1の露光中および第2の露光中に、各光源が全体として、共通周波数帯の光、第1の特有周波数帯の光および第2の特有周波数帯の光を発する光源群を起動する工程をさらに含む、請求項28に記載の方法。
  30. 第1の露光および第2の露光のそれぞれが少なくとも第4のスペクトルチャネルを含み、第1の露光および第2の露光のそれぞれに含まれるスペクトルチャネルの総数が10個以下である、請求項28または29に記載の方法。
  31. 環境光を使用して第3の露光を撮影する工程、ならびに
    第1の露光および第2の露光から第3の露光を差し引く工程
    をさらに含む、請求項28~30のいずれか1項に記載の方法。
  32. 第1の露光および第2の露光から第3の露光を差し引いた結果に基づいて、前記マルチスペクトルデータキューブを生成する工程をさらに含む、請求項31に記載の方法。
  33. 第1のスペクトルチャネルの情報に基づいて、第1のセンサ領域により撮影された画像と第2のセンサ領域により撮影された画像の位置合わせを行う工程をさらに含む、請求項28~32のいずれか1項に記載の方法。
  34. 第1のセンサ領域により撮影された画像と第2のセンサ領域により撮影された画像の位置合わせの結果に基づいて、前記マルチスペクトルデータキューブを生成する工程をさらに含む、請求項33に記載の方法。
  35. 前記画像のそれぞれが、第1のスペクトルチャネルにおいて第1のセンサ領域により撮影された情報および第1のスペクトルチャネルにおいて第2のセンサ領域により撮影された情報を示し、
    前記方法が、
    前記共通周波数帯に対応する第1のスペクトルチャネルにおいて、第1のセンサ領域により撮影された画像と第2のセンサ領域により撮影された画像の間の視差を特定する工程;
    前記視差を使用して、第2のスペクトルチャネルにおいて第1のセンサ領域により撮影された情報を示す第1の追加の画像を前記複数の画像に位置合わせする工程;ならびに
    前記視差を使用して、第3のスペクトルチャネルにおいて第2のセンサ領域により撮影された情報を示す第2の追加の画像を前記複数の画像に位置合わせする工程
    をさらに含む、請求項33または34に記載の方法。
  36. 各露光が、4個以上のスペクトルチャネルを含み;
    前記方法が、前記視差を使用して、第4のスペクトルチャネルにおいて第3のセンサ領域により撮影された情報を示す少なくとも第3の追加の画像を前記複数の画像に位置合わせする工程をさらに含み、
    前記複数の画像に位置合わせされる追加の画像の総数が、スペクトルチャネルの総数よりも1つ少ない、
    請求項35に記載の方法。
  37. 前記マルチスペクトル画像システムが、1つ以上の点、縞模様、グリッド、ランダムな斑点またはその他の空間パターンを投影するように構成された光源をさらに含み、
    前記方法が、
    前記光源によって、1つ以上の点、縞模様、グリッド、ランダムな斑点またはその他の空間パターンを含む光を、少なくとも第1のアパーチャの視野内に投影する工程;および
    前記光源により投影された光を利用して、追加の露光を撮影する工程
    をさらに含み、
    前記位置合わせが、この追加の露光に少なくとも部分的に基づく、
    請求項33~36のいずれか1項に記載の方法。
  38. 前記光源により投影された前記光が、前記共通周波数帯の波長を有する光、広帯域照明の波長を有する光および累積的に可視可能な照明の波長を有する光のうちの少なくとも1つを含む、請求項37に記載の方法。
  39. 前記追加の露光が、前記画像の前記位置合わせの正確度の確認に使用される、請求項37に記載の方法。
  40. 前記追加の露光が、前記位置合わせの正確度を向上させるために前記視差の計算に含まれる、請求項37に記載の方法。
  41. 前記光源が、前記共通周波数帯、広帯域照明または累積的に可視可能な照明からの、特有の構成の、複数の点、縞模様、グリッド、斑点またはその他の空間パターンを投影する、請求項37に記載の方法。
  42. 前記マルチスペクトル画像システムが、前記少なくとも1つのイメージセンサの視野内の物体の形状を検出するように構成された追加のセンサをさらに含み、前記位置合わせが、この追加のセンサにより検出された形状に少なくとも部分的に基づく、請求項33~41のいずれか1項に記載の方法。
  43. 前記光源が、前記マルチスペクトル画像システムの複数のアパーチャの視野内に前記光を投影する、請求項33~42のいずれか1項に記載の方法。
  44. 第1のスペクトルチャネル、第2のスペクトルチャネルおよび第3のスペクトルチャネルを分離するために、スペクトルアンミキシングを行う工程をさらに含む、請求項28~43のいずれか1項に記載の方法。
  45. 各露光が、4個以上のスペクトルチャネルを含み;前記方法が、これらのスペクトルチャネルのすべてを分離するために、スペクトルアンミキシングを行う工程をさらに含む、請求項37に記載の方法。
  46. 前記視差を使用して、トポグラフィー情報および深さ情報の少なくとも一方を測定する工程をさらに含む、請求項35~45のいずれか1項に記載の方法。
  47. 前記視差を使用して、個々の画素の深さ情報を測定する、請求項46に記載の方法。
  48. 前記深さ情報を使用して、画像化した創傷の体積を測定する工程をさらに含む、請求項47に記載の方法。
  49. 非可視光を示す第4のスペクトルチャネルを少なくとも含む追加の露光を撮影する工程;および
    前記追加の露光に少なくとも部分的に基づいて前記マルチスペクトルデータキューブを生成する工程
    をさらに含む、請求項44または45に記載の方法。
  50. 前記追加の露光中に前記非可視光を発する第3の光源群を起動する工程をさらに含む、請求項49に記載の方法。
  51. 自動組織分析方法であって、
    請求項1~27のいずれか1項に記載のマルチスペクトル画像システムを使用して、画像化された組織を示すマルチスペクトルデータキューブを生成する工程;
    創傷または熱傷などの組織の少なくとも1つの臨床的特徴を予測するように訓練された機械学習システムへの入力として前記マルチスペクトルデータキューブを提供する工程;および
    前記機械学習システムの出力に基づいて、前記画像化された組織の少なくとも1つの臨床的特徴を予測する工程
    を含む方法。
  52. 請求項28~50のいずれか1項に記載の方法を使用して、前記マルチスペクトルデータキューブを生成する工程をさらに含む、請求項51に記載の方法。
  53. 前記少なくとも1つの臨床的特徴が、複数の創傷状態または熱傷状態を含み、前記方法が、前記機械学習システムの前記出力を使用して、前記画像化された組織を、前記複数の創傷状態または熱傷状態の様々な種類ごとに分類した領域として示した画像を生成する工程をさらに含む、請求項51に記載の方法。
  54. 前記複数の熱傷状態が、少なくともI度熱傷、II度熱傷、III度熱傷または正常皮膚という種類を含むか、または前記複数の創傷状態が、少なくとも止血、炎症、増殖、リモデリングまたは正常皮膚という種類を含み、
    前記方法が、前記画像化された組織を示す複数の画素のそれぞれを、I度熱傷、II度熱傷、III度熱傷または正常皮膚という種類に分類する工程をさらに含むか、または前記方法が、前記画像化された組織を示す複数の画素のそれぞれを、止血、炎症、増殖、リモデリングまたは正常皮膚という種類に分類する工程をさらに含む、
    請求項53に記載の方法。
  55. 前記少なくとも1つの臨床的特徴が治癒可能性を含み、前記方法が、前記機械学習システムの前記出力を使用して、前記画像化された組織の様々な領域において判定された治癒可能性を視覚的に示した、前記画像化された組織の画像を生成する工程をさらに含む、請求項54に記載の方法。
  56. 前記少なくとも1つの臨床的特徴が、前記画像化された組織のがん状態または非がん状態を含み、前記方法が、前記画像化された組織を示す複数の画素のそれぞれを、がん状態または非がん状態として分類する工程をさらに含む、請求項51~55のいずれか1項に記載の方法。
  57. 前記少なくとも1つの臨床的特徴が、創傷の深さまたは体積を含み、前記方法が、前記画像化された組織における創傷の深さまたは体積を特定する情報を出力する工程をさらに含む、請求項51~56のいずれか1項に記載の方法。
  58. 前記少なくとも1つの臨床的特徴が、デブリードマンを行う際のマージンを含み、前記方法が、前記画像化された組織における創傷またはその周囲のデブリードマンを行う際のマージンを特定する情報を出力する工程をさらに含む、請求項51~57のいずれか1項に記載の方法。
  59. 前記少なくとも1つの臨床的特徴が、糖尿病性潰瘍、非糖尿病性潰瘍または慢性潰瘍を含む、請求項51~58のいずれか1項に記載の方法。
  60. 前記少なくとも1つの臨床的特徴が、糖尿病性足部潰瘍を含む、請求項59に記載の方法。
  61. 自動生体認証方法であって、
    請求項1~27のいずれか1項に記載のマルチスペクトル画像システムを使用して、ユーザの生体的特徴を示すマルチスペクトルデータキューブを生成する工程;
    特定の生体的特徴を認証するように訓練された機械学習システムへの入力として前記マルチスペクトルデータキューブを提供する工程;および
    前記機械学習システムの出力に基づいて、前記ユーザを認証または拒否する工程
    を含む方法。
  62. 請求項28~50のいずれか1項に記載の方法を使用して、前記マルチスペクトルデータキューブを生成する工程をさらに含む、請求項61に記載の方法。
  63. 自動材料分析方法であって、
    請求項1~27のいずれか1項に記載のマルチスペクトル画像システムを使用して、画像化された材料を示すマルチスペクトルデータキューブを生成する工程;
    1つ以上の材料特性を測定するように訓練された機械学習システムへの入力として前記マルチスペクトルデータキューブを提供する工程;および
    前記機械学習システムの出力に基づいて、前記画像化された材料の1つ以上の材料特性を測定する工程
    を含む方法。
  64. 請求項28~50のいずれか1項に記載の方法を使用して、前記マルチスペクトルデータキューブを生成する工程をさらに含む、請求項63に記載の方法。
  65. カラーフィルタアレイを有する少なくとも1つのイメージセンサ領域を含む携帯型のマルチアパーチャイメージング装置を使用して、該カラーフィルタアレイ中のカラーフィルタの色の数よりもチャネルの数が多い画像情報を撮影する工程;および
    前記画像情報を処理して、各チャネルが0nmを超えかつ約40nm以下の範囲のそれぞれ異なる周波数帯の光を示す、前記チャネル数を含むマルチスペクトルデータキューブを作成する工程
    を含む方法。
  66. 前記携帯型のマルチアパーチャイメージング装置が、請求項1~27のいずれか1項に記載のマルチスペクトル画像システムを含む、請求項65に記載の方法。
  67. 請求項28~50のいずれか1項に記載の方法を使用して、前記マルチスペクトルデータキューブを生成する工程をさらに含む、請求項65または66に記載の方法。
  68. 請求項1~27のいずれか1項に記載のマルチスペクトル画像システムを使用して、マルチスペクトルデータキューブを生成する工程;および
    前記マルチスペクトル画像システムの少なくとも1つのプロセッサを使用して、前記マルチスペクトルデータキューブを処理することによって、該マルチスペクトルデータキューブにより示される物体の少なくとも1つの特性を測定する工程
    を含む方法。
  69. 請求項28~50のいずれか1項に記載の方法を使用して、前記マルチスペクトルデータキューブを生成する工程をさらに含む、請求項68に記載の方法。
  70. 請求項1~27のいずれか1項に記載のマルチスペクトル画像システムを使用して、マルチスペクトルデータキューブを生成する工程;
    ネットワーク上の少なくとも1つの遠隔プロセッサに前記マルチスペクトルデータキューブを送信する工程;および
    前記少なくとも1つの遠隔プロセッサを使用して、前記マルチスペクトルデータキューブを処理することによって、該マルチスペクトルデータキューブにより示される物体の少なくとも1つの特性を測定する工程
    を含む方法。
  71. 請求項28~50のいずれか1項に記載の方法を使用して、前記マルチスペクトルデータキューブを生成する工程をさらに含む、請求項70に記載の方法。
  72. 前記マルチスペクトル画像システムが、たとえば5.9×7.9インチといった広範囲な組織から6秒以内に画像を取得するように構成されており、かつ複数の熱傷状態、複数の創傷状態もしくは治癒可能性の識別;または前記画像化された組織のがん状態もしくは非がん状態、創傷の深さ、デブリードマンの際のマージン、糖尿病性潰瘍、非糖尿病性潰瘍もしくは慢性潰瘍の存在を含む臨床的特徴などの組織分析情報を、イメージング造影剤の非存在下で提供するように構成されている、請求項1~27のいずれか1項に記載のマルチスペクトル画像システム。
  73. 前記マルチスペクトル画像システムが、たとえば5.9×7.9インチといった広範囲な組織から6秒以内に画像を取得し、かつ複数の熱傷状態、複数の創傷状態もしくは治癒可能性の識別;または前記画像化された組織のがん状態または非がん状態、創傷の深さ、創傷の体積、デブリードマンの際のマージン、糖尿病性潰瘍、非糖尿病性潰瘍もしくは慢性潰瘍の存在を含む臨床的特徴などの組織分析情報を、イメージング造影剤の非存在下で出力する、請求項43~60のいずれか1項に記載の方法。
JP2021533780A 2018-12-14 2019-12-11 高精度マルチアパーチャスペクトルイメージングのためのシステムおよび方法 Active JP7186298B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2022188817A JP7529752B2 (ja) 2018-12-14 2022-11-28 高精度マルチアパーチャスペクトルイメージングのためのシステムおよび方法
JP2022188833A JP7529753B2 (ja) 2018-12-14 2022-11-28 高精度マルチアパーチャスペクトルイメージングのためのシステムおよび方法

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US201862780121P 2018-12-14 2018-12-14
US62/780,121 2018-12-14
US201962818375P 2019-03-14 2019-03-14
US62/818,375 2019-03-14
PCT/US2019/065818 WO2020123722A1 (en) 2018-12-14 2019-12-11 System and method for high precision multi-aperture spectral imaging

Related Child Applications (2)

Application Number Title Priority Date Filing Date
JP2022188833A Division JP7529753B2 (ja) 2018-12-14 2022-11-28 高精度マルチアパーチャスペクトルイメージングのためのシステムおよび方法
JP2022188817A Division JP7529752B2 (ja) 2018-12-14 2022-11-28 高精度マルチアパーチャスペクトルイメージングのためのシステムおよび方法

Publications (2)

Publication Number Publication Date
JP2022513847A true JP2022513847A (ja) 2022-02-09
JP7186298B2 JP7186298B2 (ja) 2022-12-08

Family

ID=71077492

Family Applications (3)

Application Number Title Priority Date Filing Date
JP2021533780A Active JP7186298B2 (ja) 2018-12-14 2019-12-11 高精度マルチアパーチャスペクトルイメージングのためのシステムおよび方法
JP2022188833A Active JP7529753B2 (ja) 2018-12-14 2022-11-28 高精度マルチアパーチャスペクトルイメージングのためのシステムおよび方法
JP2022188817A Active JP7529752B2 (ja) 2018-12-14 2022-11-28 高精度マルチアパーチャスペクトルイメージングのためのシステムおよび方法

Family Applications After (2)

Application Number Title Priority Date Filing Date
JP2022188833A Active JP7529753B2 (ja) 2018-12-14 2022-11-28 高精度マルチアパーチャスペクトルイメージングのためのシステムおよび方法
JP2022188817A Active JP7529752B2 (ja) 2018-12-14 2022-11-28 高精度マルチアパーチャスペクトルイメージングのためのシステムおよび方法

Country Status (7)

Country Link
US (2) US11631164B2 (ja)
EP (1) EP3899463A4 (ja)
JP (3) JP7186298B2 (ja)
KR (1) KR20210099126A (ja)
CN (1) CN113260835A (ja)
BR (1) BR112021011113A2 (ja)
WO (1) WO2020123722A1 (ja)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3212057B1 (en) 2014-10-29 2021-12-01 Spectral MD, Inc. Reflective mode multi-spectral time-resolved optical imaging methods and apparatuses for tissue classification
CN110573066A (zh) 2017-03-02 2019-12-13 光谱Md公司 用于多光谱截肢部位分析的机器学习系统和技术
US10783632B2 (en) 2018-12-14 2020-09-22 Spectral Md, Inc. Machine learning systems and method for assessment, healing prediction, and treatment of wounds
EP3899463A4 (en) 2018-12-14 2022-12-21 Spectral MD, Inc. SYSTEM AND METHOD FOR HIGH-PRECISION MULTIAPERTURE SPECTRAL IMAGING
US10740884B2 (en) 2018-12-14 2020-08-11 Spectral Md, Inc. System and method for high precision multi-aperture spectral imaging
WO2020123724A1 (en) 2018-12-14 2020-06-18 Spectral Md, Inc. Machine learning systems and methods for assessment, healing prediction, and treatment of wounds
JP6945195B2 (ja) * 2019-01-16 2021-10-06 パナソニックIpマネジメント株式会社 光学フィルタ、光検出装置、および光検出システム
CN113820012B (zh) * 2020-06-19 2022-07-22 华为技术有限公司 一种多光谱图像处理方法及装置
CN116568214A (zh) * 2020-10-07 2023-08-08 光谱公司 使用光谱传感器系统进行健康分析
WO2022266145A1 (en) * 2021-06-18 2022-12-22 Canfield Scientific, Incorporated Multi-modal skin imaging
US20240335119A1 (en) * 2021-08-03 2024-10-10 Jubaan Ltd System and method for using multispectral imaging and deep learning to detect gastrointestinal pathologies by capturing images of a human tongue
KR102613052B1 (ko) * 2021-09-29 2023-12-12 삼성전자주식회사 멀티 스펙트럼 이미지 센서 및 이를 포함하는 전자 장치
CN113827234B (zh) * 2021-11-25 2022-03-18 之江实验室 一种基于高光谱人脸视频的远程脉搏波重建方法
CN114845061A (zh) * 2022-04-02 2022-08-02 维沃移动通信有限公司 图像传感器的控制方法、装置、图像传感器及电子设备
WO2024030297A1 (en) * 2022-08-05 2024-02-08 Battelle Memorial Institute Broadband radio frequency imaging surface
CN116760962B (zh) * 2023-08-17 2023-10-27 鸥朗创想(北京)科技有限公司 一种获取植物4d表型的无人车及方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014007869A2 (en) * 2012-06-05 2014-01-09 Hypermed Imaging, Inc. Methods and apparatus for coaxial imaging of multiple wavelengths
EP2955496A2 (en) * 2014-05-01 2015-12-16 Rebellion Photonics, Inc. Dual band divided aperture infrared spectral imager (daisi) for chemical detection
WO2016088483A1 (ja) * 2014-12-01 2016-06-09 ソニー株式会社 画像処理装置と画像処理方法
US20160349228A1 (en) * 2015-05-29 2016-12-01 Rebellion Photonics, Inc. Hydrogen sulfide imaging system
WO2017026296A1 (ja) * 2015-08-10 2017-02-16 株式会社リコー 試料測定装置
WO2017223206A1 (en) * 2016-06-21 2017-12-28 Sri International Hyperspectral imaging methods and apparatuses
WO2018160963A1 (en) * 2017-03-02 2018-09-07 Spectral Md, Inc. Machine learning systems and techniques for multispectral amputation site analysis

Family Cites Families (167)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4170987A (en) 1977-11-28 1979-10-16 California Institute Of Technology Medical diagnosis system and method with multispectral imaging
US4693255A (en) 1985-04-22 1987-09-15 Beall Harry C Medical apparatus method for assessing the severity of certain skin traumas
ATE145563T1 (de) 1989-12-14 1996-12-15 Elof Eriksson Behandlungssystem für wunden und andere fehler
US5074306A (en) 1990-02-22 1991-12-24 The General Hospital Corporation Measurement of burn depth in skin
US5701902A (en) 1994-09-14 1997-12-30 Cedars-Sinai Medical Center Spectroscopic burn injury evaluation apparatus and method
GB9521784D0 (en) 1995-10-24 1996-01-03 Rosslyn Medical Ltd Diagnostic apparatus
US6081612A (en) 1997-02-28 2000-06-27 Electro Optical Sciences Inc. Systems and methods for the multispectral imaging and characterization of skin tissue
US6008889A (en) 1997-04-16 1999-12-28 Zeng; Haishan Spectrometer system for diagnosis of skin disease
US6058352A (en) 1997-07-25 2000-05-02 Physical Optics Corporation Accurate tissue injury assessment using hybrid neural network analysis
US6353753B1 (en) 1998-05-05 2002-03-05 Stephen Thomas Flock Optical imaging of deep anatomic structures
US6352517B1 (en) 1998-06-02 2002-03-05 Stephen Thomas Flock Optical monitor of anatomical movement and uses thereof
US5982497A (en) 1998-07-09 1999-11-09 Optical Insights, Llc Multi-spectral two-dimensional imaging spectrometer
DE19850350C1 (de) 1998-11-02 2000-09-28 Jena Optronik Gmbh Verfahren und Vorrichtung zur Erzeugung von Daten für die Diagnose des Schädigungsgrades von Hautgewebe eines Patienten
US6381488B1 (en) 1999-06-15 2002-04-30 Sandia Corporation Method and apparatus to measure the depth of skin burns
US6640132B1 (en) 1999-11-17 2003-10-28 Hypermed, Inc. Forensic hyperspectral apparatus and method
AU2001229916A1 (en) 2000-01-27 2001-08-07 National Research Council Of Canada Visible-near infrared spectroscopy in burn injury assessment
US6889075B2 (en) 2000-05-03 2005-05-03 Rocky Mountain Biosystems, Inc. Optical imaging of subsurface anatomical structures and biomolecules
US6638668B2 (en) 2000-05-12 2003-10-28 Ocean Optics, Inc. Method for making monolithic patterned dichroic filter detector arrays for spectroscopic imaging
US20050187478A1 (en) 2001-07-16 2005-08-25 Art, Advanced Research Technologies Inc. Multi-wavelength imaging of highly turbid media
JP2005502895A (ja) * 2001-09-19 2005-01-27 ジュール マイクロシステムズ カナダ インコーポレイテッド 信号整合フィルタリングを組み込んだ分光計
US20060072109A1 (en) 2004-09-03 2006-04-06 Andrew Bodkin Hyperspectral imaging systems
TW512058B (en) 2001-12-24 2002-12-01 Yung-Jian Yang Spectral analysis system for burning and scalding injuries and device used in the system
EP1532431A4 (en) 2002-07-09 2010-03-31 Medispectra Inc METHODS AND APPARATUSES FOR CHARACTERIZING TISSUE SAMPLES
WO2004051242A1 (en) 2002-12-02 2004-06-17 Erasmus Universiteit Rotterdam Use of high wavenumber raman spectroscopy for measuring tissue
US20050033145A1 (en) 2003-07-02 2005-02-10 Graham John S. Wearable tissue viability diagnostic unit
US7693069B2 (en) 2003-07-28 2010-04-06 Alcatel-Lucent Usa Inc. Method, apparatus and system for improved inter-domain routing convergence
RU2372117C2 (ru) 2003-09-18 2009-11-10 Аркюо Медикал, Инк. Способ опто-термо-механического воздействия на биологическую ткань и устройство для его осуществления
US8634607B2 (en) 2003-09-23 2014-01-21 Cambridge Research & Instrumentation, Inc. Spectral imaging of biological samples
US7920908B2 (en) 2003-10-16 2011-04-05 David Hattery Multispectral imaging for quantitative contrast of functional and structural features of layers inside optically dense media such as tissue
WO2005046248A1 (ja) 2003-11-11 2005-05-19 Olympus Corporation マルチスペクトル画像撮影装置
US7433042B1 (en) 2003-12-05 2008-10-07 Surface Optics Corporation Spatially corrected full-cubed hyperspectral imager
US7648808B2 (en) 2004-01-12 2010-01-19 Ocean Thin Films, Inc. Patterned coated dichroic filter
US8583216B2 (en) 2004-04-27 2013-11-12 University Of Utah Research Foundation Skin type assessment and nevi screening for skin cancer with a noninvasive, portable reflectance spectrophotometer
US20070179482A1 (en) 2004-05-07 2007-08-02 Anderson Robert S Apparatuses and methods to treat biological external tissue
SE0402576D0 (sv) 2004-10-25 2004-10-25 Forskarpatent I Uppsala Ab Multispectral and hyperspectral imaging
US8374682B2 (en) 2005-04-04 2013-02-12 Hypermed Imaging, Inc. Hyperspectral imaging in diabetes and peripheral vascular disease
US8224425B2 (en) 2005-04-04 2012-07-17 Hypermed Imaging, Inc. Hyperspectral imaging in diabetes and peripheral vascular disease
US8548570B2 (en) 2004-11-29 2013-10-01 Hypermed Imaging, Inc. Hyperspectral imaging of angiogenesis
EP1835845A4 (en) 2004-12-28 2008-09-10 Hypermed Inc HYPERSPECTRAL / MULTISPECTRAL IMAGERY FOR THE DETERMINATION, EVALUATION AND MONITORING OF SYSTEMIC PHYSIOLOGY OR IMPACT CONDITIONS
US7729750B2 (en) 2005-01-20 2010-06-01 The Regents Of The University Of California Method and apparatus for high resolution spatially modulated fluorescence imaging and tomography
US20060241495A1 (en) 2005-03-23 2006-10-26 Eastman Kodak Company Wound healing monitoring and treatment
US8971984B2 (en) 2005-04-04 2015-03-03 Hypermed Imaging, Inc. Hyperspectral technology for assessing and treating diabetic foot and tissue disease
US8106956B2 (en) 2005-06-27 2012-01-31 Nokia Corporation Digital camera devices and methods for implementing digital zoom in digital camera devices and corresponding program products
EP1952106B1 (en) 2005-11-04 2014-08-27 George Themelis System for multispectral imaging
CA2636174C (en) 2006-01-05 2022-07-05 University Of Virginia Patent Foundation Method, system and computer program product for evaluation of blood glucose variability in diabetes from self-monitoring data
JP4655991B2 (ja) 2006-04-21 2011-03-23 カシオ計算機株式会社 撮像装置及び電子ズーム方法と、プログラム
EP2077753B1 (en) 2006-09-12 2012-12-19 The General Hospital Corporation Apparatus, probe and method for providing depth assessment in an anatomical structure
US20120321759A1 (en) 2007-01-05 2012-12-20 Myskin, Inc. Characterization of food materials by optomagnetic fingerprinting
BRPI0806109A2 (pt) 2007-01-05 2011-08-30 Myskin Inc sistema, dispositivo e método para imagem dérmica
ATE539553T1 (de) 2007-02-26 2012-01-15 Olympus Medical Systems Corp Kapselendoskop
US8639309B2 (en) 2007-07-31 2014-01-28 J&M Shuler, Inc. Method and system for monitoring oxygenation levels of compartments and tissue
US7879020B1 (en) 2007-03-19 2011-02-01 The United States Of America As Represented By The Secretary Of The Army Decision-assist method for resuscitation of patients
JP2008283442A (ja) 2007-05-10 2008-11-20 Olympus Corp 撮像装置
WO2009008745A2 (en) 2007-07-06 2009-01-15 Industrial Research Limited Laser speckle imaging systems and methods
US20090072142A1 (en) 2007-09-14 2009-03-19 Forensicare Incorporated Scanning system and techniques for medical and/or forensic assessment using the same
US20090118600A1 (en) 2007-11-02 2009-05-07 Ortiz Joseph L Method and apparatus for skin documentation and analysis
US8509879B2 (en) 2007-11-06 2013-08-13 The Regents Of The University Of California Apparatus and method for widefield functional imaging (WiFI) using integrated structured illumination and laser speckle imaging
US20100210931A1 (en) 2008-04-04 2010-08-19 Modulate Imaging Inc. Method for performing qualitative and quantitative analysis of wounds using spatially structured illumination
US8812083B2 (en) 2008-04-21 2014-08-19 Philadelphia Health & Education Corporation Methods of optically monitoring wound healing
US8838211B2 (en) 2008-04-30 2014-09-16 Board Of Regents, The University Of Texas System Multi-wavelength diagnostic imager
WO2009135081A2 (en) 2008-04-30 2009-11-05 Board Of Regents, The University Of Texas System Integrated patient bed system
CA2724973C (en) 2008-05-20 2015-08-11 University Health Network Device and method for fluorescence-based imaging and monitoring
WO2009142758A1 (en) 2008-05-23 2009-11-26 Spectral Image, Inc. Systems and methods for hyperspectral medical imaging
US8694266B2 (en) 2008-06-05 2014-04-08 The Regents Of The University Of Michigan Multimodal spectroscopic systems and methods for classifying biological tissue
US9572494B2 (en) 2008-08-12 2017-02-21 New Jersy Institute of Technology Method and apparatus for multi-spectral imaging and analysis of skin lesions and biological tissues
JP2010043979A (ja) 2008-08-13 2010-02-25 Yuichi Kamata 分光画像計測装置
US8488863B2 (en) 2008-11-06 2013-07-16 Los Alamos National Security, Llc Combinational pixel-by-pixel and object-level classifying, segmenting, and agglomerating in performing quantitative image analysis that distinguishes between healthy non-cancerous and cancerous cell nuclei and delineates nuclear, cytoplasm, and stromal material objects from stained biological tissue materials
US9206482B2 (en) 2009-04-29 2015-12-08 Genomedx Biosciences Inc. Systems and methods for expression-based classification of thyroid tissue
US8330832B2 (en) 2009-06-23 2012-12-11 Research In Motion Limited Adjustment of sharpness during digital zoom in digital photography
CN101627902B (zh) 2009-07-15 2011-12-28 深圳先进技术研究院 基于环境光的低功耗、高精度光电容积描记信号前端处理模块
JP5715132B2 (ja) 2009-08-20 2015-05-07 コーニンクレッカ フィリップス エヌ ヴェ 画像解析に関する方法及びシステム
US8666116B2 (en) 2009-10-06 2014-03-04 Koninklijke Philips N.V. Method and system for obtaining a first signal for analysis to characterize at least one periodic component thereof
US9220412B2 (en) 2009-11-19 2015-12-29 Modulated Imaging Inc. Method and apparatus for analysis of turbid media via single-element detection using structured illumination
US9823127B2 (en) 2010-01-22 2017-11-21 Duke University Systems and methods for deep spectroscopic imaging of biological samples with use of an interferometer and spectrometer
US8692912B2 (en) 2010-01-25 2014-04-08 Pixelteq, Inc. Multi-spectral camera
JP2011188481A (ja) 2010-02-10 2011-09-22 Panasonic Corp 撮像装置
US8294570B2 (en) 2010-02-24 2012-10-23 Clawson Jeffrey J Burn diagnostic and intervention tool for emergency dispatch
US20130051651A1 (en) 2010-05-07 2013-02-28 Purdue Research Foundation Quantitative image analysis for wound healing assay
JP5891424B2 (ja) * 2010-09-17 2016-03-23 パナソニックIpマネジメント株式会社 立体映像作成装置および立体映像作成方法
US20120078088A1 (en) 2010-09-28 2012-03-29 Point of Contact, LLC. Medical image projection and tracking system
CN103327894B (zh) 2011-01-19 2016-05-04 加利福尼亚大学董事会 组织的血氧定量测定和灌注成像的设备、系统及方法
US8761853B2 (en) 2011-01-20 2014-06-24 Nitto Denko Corporation Devices and methods for non-invasive optical physiological measurements
US8994819B2 (en) 2011-02-04 2015-03-31 Raytheon Company Integrated optical detection system
US8694067B2 (en) 2011-02-15 2014-04-08 General Electric Company Sensor, apparatus and method for non-invasively monitoring blood characteristics of a subject
US9372118B1 (en) 2011-03-07 2016-06-21 Fluxdata, Inc. Apparatus and method for multispectral based detection
KR101983402B1 (ko) 2011-03-07 2019-05-28 웨이크 포리스트 유니버시티 헬스 사이언시즈 전달 시스템
JP5704984B2 (ja) * 2011-03-22 2015-04-22 キヤノン株式会社 撮像装置
US20120288230A1 (en) 2011-05-13 2012-11-15 Kestrel Labs, Inc. Non-Reflective Optical Connections in Laser-Based Photoplethysmography
US8891087B2 (en) 2011-06-01 2014-11-18 Digital Light Innovations System and method for hyperspectral imaging
JP6131250B2 (ja) * 2011-06-24 2017-05-17 ケーエルエー−テンカー コーポレイション 光ルミネセンス画像化を使用する発光半導体デバイスの検査の方法および装置
US9270875B2 (en) 2011-07-20 2016-02-23 Broadcom Corporation Dual image capture processing
KR101134770B1 (ko) 2011-07-22 2012-04-13 양완석 휴대용 인체외시경 영상장치
US8942459B2 (en) 2011-09-12 2015-01-27 Perkinelmer Cellular Technologies Germany Gmbh Methods and apparatus for fast identification of relevant features for classification or regression
US9295402B1 (en) 2012-01-17 2016-03-29 University Of Washington Through Its Center For Commercialization Methods and systems for assessing a burn injury
US20150011892A1 (en) 2012-01-20 2015-01-08 Harvard Apparatus Regenerative Technology, Inc. Methods for evaluating tissue injuries
WO2013116316A1 (en) 2012-01-30 2013-08-08 Scanadu Incorporated Hyperspectral imaging systems, units, and methods
WO2013160861A1 (en) 2012-04-27 2013-10-31 Aïmago S.A. Optical coherent imaging medical device
US9031306B2 (en) 2012-05-02 2015-05-12 The Regents Of The University Of California Diagnostic and prognostic histopathology system using morphometric indices
US9593982B2 (en) * 2012-05-21 2017-03-14 Digimarc Corporation Sensor-synchronized spectrally-structured-light imaging
US9766382B2 (en) * 2012-06-05 2017-09-19 Hypermed Imaging, Inc. Single-sensor hyperspectral imaging device
US20130331669A1 (en) 2012-06-11 2013-12-12 Raytheon Company Multi-spectral imaging system and method for remote biometric measurement of human physiological parameters
US9036877B2 (en) 2012-06-20 2015-05-19 Xerox Corporation Continuous cardiac pulse rate estimation from multi-channel source video data with mid-point stitching
US9547178B2 (en) 2012-08-15 2017-01-17 Semrock. Inc. Dichroic image splitter
MX359713B (es) 2012-09-13 2018-10-08 Polyheal Ltd Composiciones mejoradas para la cicatrizacion de heridas que comprenden microesferas.
KR102251749B1 (ko) * 2012-11-07 2021-05-13 모듈레이티드 이미징, 아이엔씨. 효율적인 변조 이미지 진단
WO2014110025A1 (en) 2013-01-10 2014-07-17 Caliper Life Sciences, Inc. Whole slide multispectral imaging systems and methods
WO2014110027A1 (en) 2013-01-10 2014-07-17 Caliper Life Sciences, Inc. Multispectral imaging system and methods
US20140213910A1 (en) 2013-01-25 2014-07-31 The Regents Of The University Of California Method and apparatus for performing qualitative and quantitative analysis of burn extent and severity using spatially structured illumination
CN103973957B (zh) * 2013-01-29 2018-07-06 上海八运水科技发展有限公司 双目3d相机自动调焦系统及方法
US10012537B2 (en) 2013-01-31 2018-07-03 Ventana Medical Systems, Inc. Systems and methods for calibrating, configuring and validating an imaging device or system for multiplex tissue assays
US11653874B2 (en) 2013-02-01 2023-05-23 Acceleritas Corporation Method and system for characterizing tissue in three dimensions using multimode optical measurements
WO2014121152A1 (en) 2013-02-01 2014-08-07 Farkas Daniel Method and system for characterizing tissue in three dimensions using multimode optical measurements
GB201302451D0 (en) 2013-02-12 2013-03-27 Isis Innovation Method and system for signal analysis
AU2014224230A1 (en) 2013-03-06 2015-09-24 Marika Pty Ltd Assessing optical density gradients and variations
WO2014143235A1 (en) 2013-03-14 2014-09-18 Integrated Plasmonics Corporation Ambient light assisted spectroscopy
MX2016001866A (es) 2013-08-14 2016-06-14 Univ Nanyang Tech Sistemas para valoración de revascularización.
US9823128B2 (en) * 2013-10-16 2017-11-21 Arizona Board Of Regents On Behalf Of The University Of Arizona Multispectral imaging based on computational imaging and a narrow-band absorptive filter array
CN103815875B (zh) 2013-10-28 2015-06-03 重庆西南医院 一种用于烧伤皮肤坏死深度和面积诊断的近红外光谱成像系统
WO2015066297A1 (en) 2013-10-30 2015-05-07 Worcester Polytechnic Institute System and method for assessing wound
CN105979853A (zh) 2013-12-13 2016-09-28 莱文尼尔研究有限公司 医学成像
US10258242B2 (en) * 2014-01-06 2019-04-16 The Florida International University Board Of Trustees Near-infrared optical imaging system for hemodynamic imaging, pulse monitoring, and mapping spatio-temporal features
US20150208950A1 (en) 2014-01-28 2015-07-30 The Texas A&M University System Arterial and Venous Oxygenation Method and Apparatus
WO2015116823A1 (en) 2014-01-29 2015-08-06 The Johns Hopkins University System and method for wound imaging and debridement
WO2016178055A2 (en) * 2014-02-07 2016-11-10 Santec Corporation Self-operated health monitoring device
DE102014002514B4 (de) * 2014-02-21 2015-10-29 Universität Stuttgart Vorrichtung und Verfahren zur multi- oder hyperspektralen Bildgebung und / oder zur Distanz- und / oder 2-D oder 3-D Profilmessung eines Objekts mittels Spektrometrie
EP3120122B1 (en) 2014-03-21 2019-11-13 HyperMed Imaging, Inc. Compact light sensor
DE102014106974A1 (de) 2014-05-16 2015-11-19 Cubert GmbH Räumlich und spektral auflösende Hyperspektralkamera und Verfahren
WO2015187489A1 (en) 2014-06-04 2015-12-10 University Of Massachusetts Medical School Hyperspectral imaging for prediction of skin injury after exposure to thermal energy or ionizing radiation
FI3940371T3 (fi) 2014-06-05 2023-11-20 Univ Heidelberg Menetelmä ja kuvantamislaitteisto fluoresenssi- ja heijastuskuvien saamiseksi
JP6769949B2 (ja) 2014-07-24 2020-10-14 ユニバーシティー ヘルス ネットワーク 診断目的のためのデータの収集および解析
US10117617B2 (en) 2014-10-08 2018-11-06 Revealix, Inc. Automated systems and methods for skin assessment and early detection of a latent pathogenic bio-signal anomaly
US20220142484A1 (en) 2014-10-29 2022-05-12 Spectral Md, Inc. Reflective mode multi-spectral time-resolved optical imaging methods and apparatuses for tissue classification
US9717417B2 (en) 2014-10-29 2017-08-01 Spectral Md, Inc. Reflective mode multi-spectral time-resolved optical imaging methods and apparatuses for tissue classification
EP3212057B1 (en) 2014-10-29 2021-12-01 Spectral MD, Inc. Reflective mode multi-spectral time-resolved optical imaging methods and apparatuses for tissue classification
CN107209934B (zh) * 2014-12-03 2021-03-30 文塔纳医疗系统公司 用于定量分析异质生物标志物分布的方法、系统和装置
US9990472B2 (en) 2015-03-23 2018-06-05 Ohio State Innovation Foundation System and method for segmentation and automated measurement of chronic wound images
US11037070B2 (en) 2015-04-29 2021-06-15 Siemens Healthcare Gmbh Diagnostic test planning using machine learning techniques
BR112017027043B1 (pt) * 2015-06-15 2022-12-06 Agrowing Ltd Composto de lente e método para capturar imagens multiespectrais
WO2017027839A1 (en) 2015-08-13 2017-02-16 Somagenics, Inc. Methods and compositions of short small hairpin rnas and micrornas for wound healing
US10921186B2 (en) 2015-09-22 2021-02-16 Hypermed Imaging, Inc. Methods and apparatus for imaging discrete wavelength bands using a mobile device
CA2998699A1 (en) 2015-09-23 2017-03-30 Novadaq Technologies ULC Methods and systems for assessing healing of tissue
CN108471949B (zh) 2015-10-28 2021-11-12 光谱Md公司 组织分类用反射模式多光谱时间分辨光学成像方法和设备
NL2015804B1 (en) 2015-11-17 2017-06-02 Quest Photonic Devices B V Hyperspectral 2D imaging device.
US10572997B2 (en) * 2015-12-18 2020-02-25 Given Imaging Ltd. System and method for detecting anomalies in an image captured in-vivo using color histogram association
US11660012B2 (en) 2016-04-15 2023-05-30 The Regents Of The University Of California Assessment of wound status and tissue viability via analysis of spatially resolved THz reflectometry maps
GB2550582B (en) 2016-05-23 2020-07-15 Bluedrop Medical Ltd A skin inspection device identifying abnormalities
JP6585006B2 (ja) * 2016-06-07 2019-10-02 株式会社東芝 撮影装置および車両
US10013811B2 (en) 2016-06-13 2018-07-03 Umm-Al-Qura University Hyperspectral image visualization in patients with medical conditions
CA3021481A1 (en) 2016-07-29 2018-02-01 Novadaq Technologies ULC Methods and systems for characterizing tissue of a subject utilizing machine learning
US10223788B2 (en) 2016-08-31 2019-03-05 International Business Machines Corporation Skin lesion segmentation using deep convolution networks guided by local unsupervised learning
US10248713B2 (en) 2016-11-30 2019-04-02 Business Objects Software Ltd. Time series analysis using a clustering based symbolic representation
EP3555580B1 (en) * 2016-12-27 2023-12-06 Urugus S.A. Hyper-spectral imaging when observed object is still
US10145740B2 (en) * 2017-02-27 2018-12-04 Visera Technologies Company Limited Sensing multiple peak wavelengths using combination of dual-band filters
US10806334B2 (en) 2017-02-28 2020-10-20 Verily Life Sciences Llc System and method for multiclass classification of images using a programmable light source
US11003933B2 (en) * 2017-08-15 2021-05-11 Noblis, Inc. Multispectral anomaly detection
AU2017415626B2 (en) 2017-10-17 2020-08-06 Kronikare Pte Ltd System and method for facilitating analysis of a wound in a target subject
CA3076483A1 (en) * 2017-11-16 2019-05-23 MultiSensor Scientific, Inc. Systems and methods for multispectral imaging and gas detection using a scanning illuminator and optical sensor
CN108053482A (zh) * 2018-02-05 2018-05-18 喻强 一种基于手机扫描的人体3d建模方法
US20190290117A1 (en) 2018-03-22 2019-09-26 Kabushiki Kaisha Topcon Interferometric fundus imaging method
US10740884B2 (en) 2018-12-14 2020-08-11 Spectral Md, Inc. System and method for high precision multi-aperture spectral imaging
WO2020123724A1 (en) 2018-12-14 2020-06-18 Spectral Md, Inc. Machine learning systems and methods for assessment, healing prediction, and treatment of wounds
EP3899463A4 (en) 2018-12-14 2022-12-21 Spectral MD, Inc. SYSTEM AND METHOD FOR HIGH-PRECISION MULTIAPERTURE SPECTRAL IMAGING
US10783632B2 (en) 2018-12-14 2020-09-22 Spectral Md, Inc. Machine learning systems and method for assessment, healing prediction, and treatment of wounds
WO2021173763A1 (en) 2020-02-28 2021-09-02 Spectral Md, Inc. Machine learning systems and methods for assessment, healing prediction, and treatment of wounds
WO2022015597A1 (en) 2020-07-13 2022-01-20 Spectral Md, Inc. Spectral imaging systems and methods for histological assessment of wounds
CN117321394A (zh) 2021-03-30 2023-12-29 光谱Md公司 基于多路照明的高精度快照多光谱成像系统和方法
CN118475953A (zh) 2022-01-21 2024-08-09 光谱Md公司 用于组织的拓扑表征的系统和方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014007869A2 (en) * 2012-06-05 2014-01-09 Hypermed Imaging, Inc. Methods and apparatus for coaxial imaging of multiple wavelengths
US20150051498A1 (en) * 2012-06-05 2015-02-19 Hypermed Imaging, Inc. Methods and apparatus for coaxial imaging of multiple wavelengths
EP2955496A2 (en) * 2014-05-01 2015-12-16 Rebellion Photonics, Inc. Dual band divided aperture infrared spectral imager (daisi) for chemical detection
WO2016088483A1 (ja) * 2014-12-01 2016-06-09 ソニー株式会社 画像処理装置と画像処理方法
US20160349228A1 (en) * 2015-05-29 2016-12-01 Rebellion Photonics, Inc. Hydrogen sulfide imaging system
WO2017026296A1 (ja) * 2015-08-10 2017-02-16 株式会社リコー 試料測定装置
WO2017223206A1 (en) * 2016-06-21 2017-12-28 Sri International Hyperspectral imaging methods and apparatuses
WO2018160963A1 (en) * 2017-03-02 2018-09-07 Spectral Md, Inc. Machine learning systems and techniques for multispectral amputation site analysis

Also Published As

Publication number Publication date
KR20210099126A (ko) 2021-08-11
BR112021011113A2 (pt) 2021-08-31
WO2020123722A1 (en) 2020-06-18
US20230206413A1 (en) 2023-06-29
JP7186298B2 (ja) 2022-12-08
JP2023041063A (ja) 2023-03-23
CN113260835A (zh) 2021-08-13
US20220156903A1 (en) 2022-05-19
JP2023037629A (ja) 2023-03-15
US11631164B2 (en) 2023-04-18
EP3899463A4 (en) 2022-12-21
US11989860B2 (en) 2024-05-21
JP7529753B2 (ja) 2024-08-06
JP7529752B2 (ja) 2024-08-06
EP3899463A1 (en) 2021-10-27

Similar Documents

Publication Publication Date Title
JP7529752B2 (ja) 高精度マルチアパーチャスペクトルイメージングのためのシステムおよび方法
JP7261883B2 (ja) 創傷の評価、治癒予測および治療のための機械学習システム
US11182888B2 (en) System and method for high precision multi-aperture spectral imaging
US11599998B2 (en) Machine learning systems and methods for assessment, healing prediction, and treatment of wounds
US20230148951A1 (en) Spectral imaging systems and methods for histological assessment of wounds
US20230181042A1 (en) Machine learning systems and methods for assessment, healing prediction, and treatment of wounds
WO2014007869A9 (en) Methods and apparatus for coaxial imaging of multiple wavelengths
US20240048823A1 (en) System and method for high precision snapshot multi-spectral imaging based on multiplexed illumination
WO2023141216A2 (en) System and method for topological characterization of tissue
JP7574354B2 (ja) 創傷の評価、治癒予測および治療のための機械学習システム

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210811

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20220630

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220802

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20221018

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20221128

R150 Certificate of patent or registration of utility model

Ref document number: 7186298

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150