Nothing Special   »   [go: up one dir, main page]

JP2010521272A - 標準化された撮像データを実現するために撮像デバイスに自動品質フィードバックを与える方法 - Google Patents

標準化された撮像データを実現するために撮像デバイスに自動品質フィードバックを与える方法 Download PDF

Info

Publication number
JP2010521272A
JP2010521272A JP2009554535A JP2009554535A JP2010521272A JP 2010521272 A JP2010521272 A JP 2010521272A JP 2009554535 A JP2009554535 A JP 2009554535A JP 2009554535 A JP2009554535 A JP 2009554535A JP 2010521272 A JP2010521272 A JP 2010521272A
Authority
JP
Japan
Prior art keywords
image
cross
polarized
polarized image
region
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.)
Pending
Application number
JP2009554535A
Other languages
English (en)
Inventor
ハーグローブ,ジョン
グ,ジア
リ,ウェンジン
ウォルタース,ロルフ・ホルガー
Original Assignee
エスティーアイ・メディカル・システムズ・エルエルシー
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 エスティーアイ・メディカル・システムズ・エルエルシー filed Critical エスティーアイ・メディカル・システムズ・エルエルシー
Publication of JP2010521272A publication Critical patent/JP2010521272A/ja
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/40Image enhancement or restoration using histogram techniques
    • 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
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/90Dynamic range modification of images or parts thereof
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/56Extraction of image or video features relating to colour
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0082Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes
    • A61B5/0084Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes for introduction into the body, e.g. by catheters
    • 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/10016Video; Image sequence
    • 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/10024Color 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/10152Varying illumination
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30168Image quality inspection

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Quality & Reliability (AREA)
  • Multimedia (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Testing, Inspecting, Measuring Of Stereoscopic Televisions And Televisions (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Studio Devices (AREA)
  • Endoscopes (AREA)

Abstract

対象領域の位置を特定する機能、画像コントラストを最大にする機能、及び対象領域が正確に画像の中心に位置することを確実にする機能を実行する、自動画像品質評価アルゴリズム。その対象領域は、画像自体のサンプルから得られる目標スペクトルを用いて、スペクトル整合フィルタによって特定される。
【選択図】図1

Description

本発明は包括的には医用画像に関し、より具体的には、アーカイブ品質診療記録及びコンピュータ支援診断(CAD)システムにおいて使用するために、高品質の標準化デジタル画像を実現する画像処理に関する。本発明は、アルゴリズムを用いて、デジタル撮像装置から得られる画像の品質を解析し、画像品質を能動的に改善するためのフィードバックを生成する処理システムを提供する。
本発明は子宮頸癌との関連で開示されているが、本発明は数多くの他の医学分野にも適用することができる。子宮頸癌は世界中の女性において2番目に多い癌であり、一年間に実質的に500000人以上の新たな症例が発生し、270000人を超える女性が死亡している(非特許文献1、参照により本明細書に援用される)。コルポスコピーは、子宮頸部の前癌病変及び癌を検出するために用いられる診断方法である(参照により本明細書に援用される非特許文献2)。リード指数(Reid’s index)のようなスコアリング方式は、アセトホワイトニング(acetowhitening)、血管パターン及び病変周縁部のような種々の特徴に基づいて、コルポスコピー診断を行なうための補助手段である(参照により本明細書に援用される非特許文献2)。これらの特徴を個別に評価して、スコア化した後に、全ての特徴のスコアが組み合わされて、疾患重症度を格付けする複合スコアが生成される。しかしながら、信頼性のあるスコアリングを保証するために、さらに解析を進める前に、画像の品質を評価しなければならない。
子宮頸部画像の品質が制限されるのは、機器の設定が不正確であること、機器の配置が不正確であること、グリント、合焦不良に起因するぶれ、物理的な汚れを含む、いくつかの要因によるものと考えることができる。グリント(鏡面反射)は、影響を及ぼされるピクセル内の色情報をなくし、それゆえ、特徴抽出アルゴリズムにアーティファクトを導入する可能性がある。鏡面反射は、鏡のように表面から光を完全に反射し、単一の入射方向からの光(光線)が単一の出射方向に反射される。ピクセルはグラフィック画像内の1つの点であり、画像の最も小さな単一要素である。1つの画像内の各ピクセルは、その輝度(brightness)又は強度(intensity)と相関する独自の値を有する。カラー画像では、各ピクセルは、色相(H(hue))、彩度(S(saturation))及び明度(V(value))(HSV)、又は色相(H)、彩度(S)、明度(L(lightness))(HSL)を用いて示すことができるが、通常は、代わりに赤色、緑色及び青色(RGB)強度(intensities)として表される。色相(H)、彩度(S)及び強度(I(intensity))(HSI)並びに色相(H)、彩度(S)及び輝度(B(brightness))(HSB)は、HSV及びHSLの代わりに用いられる名称である。HSL及びHSVは、その中心軸が最も下の黒色から最も上の白色までの範囲を有し、その間に中間の色を有する円柱形内の点として色を表すために用いることができ、その軸を中心にした角度は「色相」に対応し、その軸からの距離は「彩度」に対応し、その軸に沿った距離は「明度(lightness)」、「明度(value)」及び「輝度」に対応する。
コルポスコピーのためのCADは、医用画像処理の新たな応用を表す。本発明者らは、子宮頸部異常の重症度を評価するために、リードインデックスのような診断方法に類似のCADシステムを開発した(参照により本明細書に援用される非特許文献4)。コルポスコピーの特徴を正確に評価するために、CADが処理する画像は高品質でなければならない。それゆえ、本発明は、アクティブ画像品質評価(Active Image Quality Assessment)(AIQA)と呼ばれる、アルゴリズムの体系的な枠組みを含み、それらのアルゴリズムは、対象領域の位置を特定すること、コントラスト最大化、及び構図解析を含む。これらのアルゴリズムは、デジタルコルポスコープから得られる子宮頸部画像を自動的に評価する。フィルタリングされたデータセットは、その後、CADアルゴリズム及びアーカイブ品質診療記録のために用いることができ、遠隔医療子宮頸癌診断においても用いることができる。
適当なコントラスト及び構図の要件が満たされる時点を判断するために、本発明の評価アルゴリズムは、最初に、デジタルコルポスコープからのライブビデオデータストリームからサンプリングされる低解像度の生画像に適用される。評価アルゴリズムからのフィードバックに応答して、必要に応じてシステム調整が行われる。これらの調整は、自動システム調整、及び機器の指示を受けながらオペレータによって行なわれる手動調整の双方から成る場合がある。任意の調整後に、その評価アルゴリズムは繰り返される。そのアルゴリズムが、システムに対するそれ以上の調整が不要であることを指示した後、そのシステムは、CADシステムにおいて用いるための高解像度の静止画像を収集できる状態にある。
多くの評価アルゴリズムが当該技術分野において既知であるが、本発明者らは、対象領域の位置を特定するために、グリッドボックスサンプリング法を用いて、スペクトル整合フィルタによって用いられる目標スペクトルを推定する他のいかなる画像処理方法も知らない。以下の特許及び特許出願は、本発明の分野に関連すると考えることができる。
参照により本明細書に援用される、Galに対する特許文献1は、撮像センサと、広視野光学系と、照明手段と、電子回路と、通信手段とを備える、膣検査を実行するための医療デバイスを開示する。さらに、特許文献1は、プロセッサ及び専用医用画像理解・分析(MIUA)アルゴリズムを含む、MIUAユニットを開示する。撮像センサ、光学系及び照明手段が、種々の波長において画像を収集し、マルチスペクトル画像を提供する。その撮像デバイスは、少なくとも1つのフィルタから成る。MIUAアルゴリズムは、それらの画像に対象領域が存在するかを調べて、得られた画像を所定の品質基準と比較し、得られた画像が所定の基準を満たさない場合には、照明強度、露出時間、角度及びカメラズームを調整することができる。最後に、MIUAアルゴリズムは、所定の対象領域の十分な画像が得られたことを検証し、保存し、その後、それをユーザに通知する。
参照により本明細書に援用される、Collins、Jeffery他に対する特許文献2は、医用画像分析及び悪性病変検出に関するコンピュータ支援システム並びに方法を開示する。多様なモダリティにより取得された医用画像を分析する。形態学的な特徴並びに時間的な、すなわち動力学的な特性を組み合わせて、医用画像中に検出した病変可能性の総合評価を算定する。システムは少なくとも1個の動力学的モジュールを有し、この動力学的モジュールは、患者に造影剤を投与した後に取得される時系列MRI画像又はMRSデータから動力学的特性を抽出することができるものとする。総合評価は確認又は改善のためにユーザに提示する。
参照により本明細書に援用される、Utzinger他に対する特許文献3は、組織のマルチスペクトル画像を生じさせるための方法及び装置を開示する。このマルチスペクトル画像は、頸部癌のような状態の検出及び診断のための、診断ツールとして使用され得る。一次放射線が、照射源を用いて発生される。この一次放射線は、第1の波長及び第1の偏光を選択するために、フィルタリングされる。組織が、このフィルタリングされた一次放射線で照射されて、二次放射線を発生させ、これは、第2の波長及び第2の偏光を選択するために、フィルタリングされる。このフィルタリングされた二次放射線は、検出器を用いて収集され、そして組織の複数のマルチスペクトル画像が、第1及び第2の波長並びに第1及び第2の偏光の異なる組合せに従って、検出器と作動可能な関係にある分析ユニットを用いて、生成される。
参照により本明細書に援用される、Hsieh他に対する特許文献4は、コンピュータ支援の診断アルゴリズムによって最初の画像の分析に基づいて画像データを生成及び処理するための方法及びシステムを開示する。そのCADアルゴリズムは、画像データ評価、データ又は画像セグメント化、特徴選択又は抽出、分類、トレーニング及び視覚化のようなモジュールを含むことができる。特許文献4は、セグメント化又は構造識別によって、関心のある特徴の位置を特定することを開示する。また、特許文献4は、コントラスト、空間解像度(たとえば、ズーム)及び色のような種々のパラメータに重点を置く後処理技法も開示する。
Hsieh他に対する特許文献5は、コンピュータ支援診断アルゴリズムによる最初の画像データの解析に基づいて、医療診断との関連で後続の画像データを収集し、付加的な画像データ収集が適切であるか否かを判断するための技法を開示する。
参照により本明細書に援用される、Jiangに対する特許文献6は、入力画像から明るさ(luminance)の値をフィルタリングし、フィルタリングされた値を変換して、画質を高めた画像を生成することによって、診断のための組織サンプルの画像の画質を高める方法を開示する。その発明は、画像輝度、コントラストを調整する変換アルゴリズムを使用し、カラーバランシングを提供することができる。また、特許文献6は、遮られるか、又は対象の診断ゾーンの外側にあるサンプルの領域のためのマスキング技法も開示する。
参照により本明細書に援用される、Blairに対する特許文献7は、子宮頸部のデジタル画像を生成することができる低コストのハンドヘルドコルポスコピーアセンブリを開示しており、それによれば、癌性組織を検出するために子宮頸部の画像をリアルタイムに撮像し、記録できるようになる。作成されるデジタル画像を瞬時に処理して、反射グレアを除去するか、又は任意のデジタル画像画質向上処理を実行して、組織の外見、組織及び病変の境界、並びに組織血管分布を確定することができる。
参照により本明細書に援用される、Chen、Ying-Ling Ann他に対する特許文献8は、適応的フォトスクリーニングシステムを開示する。
参照により本明細書に援用される、Cresens及びMarcに対する特許文献9は、直列番号付与済み品質管理(QC)目標を使う品質管理(QC)解析時に得られた、デジタル画像形成システムの画像品質性能への物理的特性変動の影響を、物理的特性偏差を制御される動作モデル修正を原画像品質性能に適用することにより、打ち消す方法。該QC解析用に使われる該直列番号付与済みQC目標はQC目標の内側又は外側にエンコードされた、目標固有の、測定された物理的特性データを有する。
参照により本明細書に援用される、Reiner及びBruceに対する特許文献10は、放射線医学に基づく医用画像研究においてデジタルデバイスを使用する科学技術者のための品質保証(QA)スコアカードを生成する品質保証システム及び方法を開示する。
参照により本明細書に援用される、Selse、Emil他に対する特許文献11は、患者の胸部の測位に関するデジタルマンモグラフィ画像の自動品質評価を開示する。
参照により本明細書に援用される、Krauter他に対する特許文献12は、電子緑色フィルタを備えるビデオコルポスコープのための装置及び方法を開示する。ビデオカメラが、対象物の対象電子画像を取得し、アルゴリズム駆動デジタル信号処理回路(DSP)を用いて、DSP基準フィルタアルゴリズム及び格納されている基準色バランスレベルに従って、対象電子画像の彩度、色相及び強度レベルが修正され、それにより、対象電子画像に対応する修正済み電子画像が生成される。対応する対象画像がビデオカメラから得られるのに応じて、修正済み電子画像が、連続してリアルタイムにディスプレイに出力される。この修正済み電子画像は、光学緑色フィルタを通じて得られた画像をエミュレートし、シミュレートされたホワイトバランスを組み込む。
参照により本明細書に援用される、Manianに対する特許文献13は、放射線検査画像の画像品質の迅速な自動定量評価を実行するためのシステム及び方法を開示する。
国際公開第2008/012813号明細書 米国特許出願公開第2007/0133852号明細書 米国特許第6,766,184号明細書 米国特許出願公開第2004/0068167号明細書 米国特許第6,687,329号明細書 米国特許出願公開第2004/0208385号明細書 米国特許第6,277,067号明細書 米国特許出願公開第2007/0171363号明細書 米国特許出願公開第2007/0009096号明細書 米国特許出願公開第2007/0237308号明細書 米国特許出願公開第2007/0248210号明細書 米国特許第6,147,705号明細書 米国特許第5,565,678号明細書
http://wwwdepdb.iarc.fr/globocan2002.htm B. S. Apgar、Brotzman, G, L.及びSpitzer, M.著「Colposcopy: Principle and Practice」(W. B. Saunders Company: Philadelphia, 2002) Reid R及びScalzi P.著「Genital warts and cervical cancer. VII. An improved colposcopic index for differentiating benign papillomaviral infection from high-grade cervical intraepithelial neoplasia」(Am J Obstet Gynecol 1985; 153: 611-618) Lange H.及びFerris, Daron G.著「Computer-Aided-Diagnosis (CAD) for colposcopy」(SPIE Medical Imaging 2005; SPIE Proc. 5747, 2005)
本発明は、標準化された撮像データを記録するために、撮像デバイスに自動品質フィードバックを与える方法を開示する。本発明は、光源で視野を照明して、交差偏光(XP)画像及び単一偏光(PP)画像(平行偏光画像と呼ばれることもある)を得るステップを含む。XP画像は、第1の偏光方向を有する光源で視野を照明し、第1の方向に対して実質的に垂直である第2の偏光方向を有する偏光フィルタを通じて検出器で視野を視認することによって作成されることが好ましい。PP画像は、第1の偏光方向が第2の偏光方向に実質的に平行であるときに作成されることが好ましい。また、PP画像は、光源又は検出器のいずれかに(双方ではない)単一偏光フィルタを置くことによって得ることもできる。XP画像及びPP画像の双方は、互いに重ね合わせられる(位置合わせされる)。その後、交差偏光画像の空間的に離散したサンプルの平均スペクトルを求め、最も明るいサンプルの約0%〜20%及び最も暗いサンプルの10%〜50%を廃棄して、残りのスペクトルを平均することによって目標スペクトルを規定する画像分類アルゴリズムを用いて、XP画像上で対象領域の位置が特定される。その後、目標スペクトルに対応するXP画像の部分を照合することによって、対象領域の大きさが求められる。次に、信号評価及び露出時間調整を用いて適当なコントラストを得るために、対象領域及びPP画像が解析される。その後、XP画像に構図解析を適用して、対象領域が画像境界内の正確に中央に配置されるようにする。最後に、コントラスト最大化及び構図解析のステータスが評価され、そのシステムは、「良好」又は「不良」いずれかのブール変数を報告する。そのステータスが「不良」である場合には、「良好」の評価ステータスが達成されるまで、全過程(2つのバージョンの生画像を収集し、対象領域の位置を特定し、コントラスト最大化を実行し、構図解析を適用すること)が繰り返される。ステータスが「良好」である場合には、これはコントラスト最大化及び構図解析が十分であることを意味しており、そのシステムは、対象領域の高解像度静止画像を収集できる状態にある。
アクティブ画像品質評価(AIQA)の流れ図である。 子宮頸部であるが、膣壁部、光沢のない黒い膣鏡及び外部組織の一部も含む、赤色-緑色-青色(RGB)画像の例である。 スペクトル整合フィルタのための目標スペクトルを推定するために用いられる9個のボックスを有する子宮頸部のRGB画像である。 AIQA ROIアルゴリズムを用いて計算される、図4及び図5の画像のための子宮頸部の対象領域(ROI)を示す図である。 ROIによってマスクされる子宮頸部のRGB画像である。 ROI内の赤色ピクセルのデジタル数(DN)値のヒストグラムである信号DNの最小許容値及び最大許容値が2つの垂直な破線として示されている。860DNにおける線が、コントラストアルゴリズムによって推定される信号DNの値である。この画像は許容可能なコントラストを有する。 ROIによってマスクされる子宮頸部のRGB画像である。 ROI内の赤色ピクセルのDN値のヒストグラムである。信号DNのための最小許容値及び最大許容値が2つの垂直な破線として示されている。約900DNにおける線が、コントラストアルゴリズムによって推定される信号DNの値である。この画像は露出過度である。 ROIによってマスクされる子宮頸部のRGB画像である。 ROI内の赤色ピクセルのDN値のヒストグラムである。信号DNの最小許容値及び最大許容値が2つの垂直な破線として示されている。約750DNにおける線が、コントラストアルゴリズムによって推定される信号DNの値である。この画像は露出不足である。
本発明の評価アルゴリズムが実行される前に、デジタル撮像装置を用いて、調査中のライブビデオから、視野の低解像度の生画像が収集されなければならない。2つのバージョンの視野、XP画像及びPP画像が捕捉される。これは、第1の偏光方向を有する偏光光源で視野を照明し、第1の偏光方向に対して実質的に垂直である第2の偏光方向を有する偏光フィルタを通じて検出器で視野を視認することによって達成されることが好ましい。これが、視野のXP画像を生成する。視野のPP画像を生成するために、第2のバージョンの視野は、偏光光源が、検出器フィルタの第2の偏光方向に対して実質的に平行である第1の偏光方向を有するときに、又は光源若しくは検出器のいずれかにおいてただ1つの偏光フィルタが存在するときに達成されることが好ましい。
(ほぼ無限の数の波長を有する)広帯域白色光源、並びに赤色、緑色及び青色チャネルを有する(すなわち、これら3つの波長を検出することができる)標準色センサの使用によるマルチスペクトル検出が用いられることが好ましい。
1.評価アルゴリズムシーケンス
1.1 対象領域(ROI)の位置を特定する
本発明の現時点で好ましい実施形態の流れ図が図1において与えられる。本発明は、システムが焦点サブシステムによって正確に焦点を合わせられているものと仮定する。好ましくは、3つのアルゴリズム、すなわち(1)対象領域位置特定、(2)コントラスト最大化、及び(3)構図解析が存在する。これらのアルゴリズムが適用された後に、総合ステータス評価が実行される。第1のアルゴリズムは、子宮頸部に対応するはずである対象領域(ROI)の計算である。ROIは、XP画像を用いて推定される。なぜなら、XP画像は、グリントアーティファクト(鏡面反射)の影響を受けにくいためである。鏡面反射は、鏡のような、表面からの光の完全な反射であり、単一の入射方向からの光(光線)が単一の出射方向に反射される。その後、推定されたROIは、残りの2つのアルゴリズム、すなわちコントラスト最大化及び構図解析のための入力として用いられる。
1.2. コントラスト最大化
このアルゴリズムの場合、視野のXP画像及びPP画像のためのコントラストが個別に最大にされる。偏光効果に起因して、2つの画像は異なる量の光を受け取るので、PP画像においてコントラスト最大化を実行する必要がある。結果として、各画像のコントラストを最大にするために、露出毎に個別に露出時間が調整されなければならない。XP画像及びPP画像のコントラスト最大化は、XP画像からのROI推定値を入力として利用する。
本発明のこの部分の目標は、カメラのダイナミックレンジを完全に利用するのを確実にし、それにより、画像コントラストを最大にすることである。コントラスト最大化は、信号評価及び露出時間調整を通じて達成される。これは、最も明るい子宮頸部ピクセル(生画像において最も大きなデジタル数(DN)を生成するピクセル)が、そうしなくても、カメラの飽和にできる限り近いように、カメラ露出時間が調整されることを意味する。ピクセルは、グラフィック画像内の1つの点であり、1つの画像の最も小さな単一要素である。1つの画像内の各ピクセルは、その輝度又は強度と相関する独自の値を有する。0の数値表現は通常黒色を表し、取り得る最大値は通常、白色を表す。DNは、デジタル画像内のピクセルの相対的な輝度を表す正の整数である。カラー画像では、各ピクセルは、その色相、彩度及び明度(HSV)を用いて表すことができるが、通常は、代わりに赤色、緑色及び青色(RBG)強度として表される。
コントラスト最大化アルゴリズムからの一次出力は、そのシステムが満足のいくコントラストを達成したか否かを指示するメッセージである。そのシステムが満足のいくコントラストを達成していない場合には、テスト画像に適用する新たな露出時間が与えられ、評価アルゴリズムの残りは実行されず、その画像のためのシーケンスは終了し、新たな露出時間による新たなテスト画像が得られる。一方で、許容可能なコントラストが達成されている場合には、出力される露出時間は変更されず、評価アルゴリズムは構図解析に続く。
1.3 構図解析
構図解析は、本発明の現時点で好ましい実施形態の第3のステージである。構図解析の目的は、画像境界と交差することなく、ROI(子宮頸部であることが好ましい)が各画像内の中心に位置し、できる限り視野の多くを占めるのを確実にすることである。XP画像からのROIが入力として用いられる。ROIの中心が計算され、画像の中心と比較される。ROI境界が画像境界に近いかが評価される。最後に、カメラ視野内で子宮頸部が遮られる場合を検出するために、ROIの大きさが計算される。ROIが中心に位置し、が画像境界と交差せず、且つ十分に大きい(XP画像の全面積の少なくとも25%を構成する)場合には、出力構図メッセージは、構図が良好であることを指示する。これが当てはまらない場合には、その構図は不良であり、3つの品質フラグ出力(後述する)において指示される。この情報を用いて、人間オペレータが構図に対する調整を行なうことができる。
1.4 状態評価
本発明の現時点で好ましい実施形態の最終ステージは、評価アルゴリズムのための総合状態の判断であり、これは「良好」又は「不良」のいずれかのブール変数として報告される。総合ステータスは、2つの評価アルゴリズム(コントラスト最大化及び構図解析)がいずれも満足のいくときにのみ、値「良好」(真)を割り当てられる。総合ステータスが「良好」であるとき、そのシステムは、臨床的に関連する特徴を観測するために、ROIの高解像度静止画像を収集できる状態にある。2つの評価アルゴリズム(コントラスト最大化及び構図解析)のいずれかが不満足である場合には、総合ステータスは「不良」であり、コントラストフラグ及び構図フラグに基づいてシステム調整が行われた後に、全過程(ROIの位置を特定すること、コントラスト最大化、構図解析、及びステータス評価)が繰り返される。
2 対象領域の位置を特定すること
2.1 問題の記述
ROIアルゴリズムの目的は、子宮頸部に対応する生画像内のピクセルを特定することである。その後、これらのピクセルを用いて、画像コントラスト及び構図が解析される。一例が図2として示されており、ここで、生画像が赤色-緑色-青色(RGB)画像に変換されている。子宮頸部は、その画像の中心において明らかに視認することができるが、膣壁部、光沢のない黒い膣鏡及び外部組織の一部のようないくつかの他の特徴も、その画像において視認することができる。子宮頸部ピクセルのみが、画像品質を解析するための対象である。
2.2 アルゴリズムの記述
XP画像及びPP画像が互いに重ね合わせられる(位置合わせされる)ので、単一の画像からROIが十分に計算される。1つのカメラが用いられる(画像捕捉間に時間遅延がある)場合には、その時間遅延は、画像が互いに重ね合わせられるように、組織の動きを最小限に抑えるほど十分に短くしなければならないか、又はそうしない場合には2つの画像を同時に撮影するための半透鏡シャッタビーム又は他の機構が用いられなければならない。上述したように、XP画像はグリントアーティファクトの影響を受けにくいので、ROIはXP画像を用いて推定される。生XP画像は、通常はバイエルパターンフォーマットであり、最初にRGB画像に変換されなければならない。バイエルパターンフォーマットは、カラー画像を作成するためにデジタルカメラにおいて用いられる大部分の単一チップデジタル画像センサにおいて用いられる色フィルタの特定の配列を指している。フィルタパターンは50%緑色、25%赤色及び25%青色であり、それゆえ、RGBGとも呼ばれる。赤色又は青色の2倍の数の緑色要素を用いて、人の目が緑色光に対してより高い感度を有するのを模倣する。これをRGBに変換するために、複雑なデモゼイキング過程は、計算負荷及び関連する時間遅延に起因して望ましくない。代わりに、簡単なビンニングが用いられ、それは迅速であり、且つここでこの目的を果たすのに十分である。各RGBピクセルは、生からの単一赤色ピクセル、単一青色ピクセル、及び2つの隣接する緑色ピクセルの平均から構成される。
ROIを計算するために利用される基本的な数学関数はスペクトル整合フィルタリングである。スペクトル整合フィルタリングは非常に良好に機能し、計算が速い。電気通信の場合、既知の信号、すなわちテンプレートと、未知の信号との相関をとり、未知の信号内のテンプレートの存在を検出することによって、整合フィルタが得られる。画像処理では、2次元整合フィルタが一般的に用いられる。本発明において用いられる整合フィルタは、子宮頸部に対応する目標スペクトルと、子宮頸部上に存在する可能性がある各ピクセルとの間のスペクトル角を推定する。スペクトルは、特定の1組の値に限定されるのではなく、連続体内で無限に変化することができる1つの条件である。
目標スペクトル及びスペクトルの全てのサンプルを多次元空間内で表すことができる。この空間では、スペクトルの赤色ピクセル値が赤色次元内の値に対応し、スペクトルの緑色ピクセル値が緑色次元内の値に対応し、スペクトルの青色ピクセル値が青色次元内の値に対応する。赤色値、緑色値及び青色値は合わせて、この空間内で1つのベクトルを定義する。目標スペクトル及びスペクトルの全てのサンプルは、この空間内でベクトルとして表すことができる。目標スペクトルベクトルと各サンプルスペクトルベクトルとの間のスペクトル角を評価することができる。この角度が或る最小しきい値未満であるピクセルは、子宮頸部目標ピクセルに十分に類似していると見なされ、フィルタリングに合格する。目標スペクトルとサンプルスペクトルとの間のスペクトル角は、好ましくは、約65度〜85度未満にすべきである。最適には、その角度は、約75度未満にすべきである。スペクトル角は、「白色化された」RGB空間、すなわち、ピクセルの色共分散行列(後に定義される)によって回転及び正規化されている色空間において計算されることが好ましい。
整合フィルタを推定するために、最初に不良ピクセルの組が、0に極めて近いか、又は飽和に極めて近い値を有する全てのピクセルと定義される。残りのピクセルは、これ以降、背景ピクセルと呼ばれる。最初に、不良ピクセルが廃棄され、その後、背景ピクセルの平均(平均値)スペクトルμbkg及び共分散行列Σが計算される。統計学において、共分散行列は、ベクトルの要素間の共分散の行列である。共分散は、2つの確率変数が共に変化する程度に関する指標である(単一の変数が変化する程度を評価する分散とは異なる)。2つの変数が共に変化する傾向がある場合には(すなわち、それらの変数のうちの一方が、その期待値より大きいときに、他方の変数も、その期待値より大きい傾向がある)、それらの変数間の共分散は正である。一方、それらの変数のうちの一方が、その期待値より大きいときに、他方の変数が、その期待値より小さい傾向がある場合には、2つの変数間の共分散は負になる。
次に、子宮頸部に対応する目標スペクトルが生成される。これは、XP画像におけるROI内のスペクトルの空間的に離散したサンプルを取ることを含む、グリッドボックスサンプリング法を用いることによって達成される。スペクトル整合フィルタについての従来からの問題は、画像ピクセルと比較する、又は「照合する」ために、「目標スペクトル」を得ることであるので、ROIの位置を特定するこの方法は望ましい。子宮頸部を特徴付ける固有スペクトルは、自然変動に起因して、女性によって常に異なる。画像そのものから目標スペクトルを得ることが好ましいが、ROIアルゴリズムの目的がその画像において子宮頸部を見つけることであるときに、画像そのものから子宮頸部目標スペクトルを如何に得るかという問題が残る。グリッドボックス法は、サンプリングを通じて、この問題に対処する。
図3は、離隔して配置される9個のピクセルボックス(サンプル)を示す。ボックスの大きさ及びボックス間の間隔は調整可能である。好ましくは、ボックスは均一な大きさ及び均一な間隔を有するが、必ずしもそうである必要はなく、ボックス長はXP画像の境界長の約10%〜30%であり、ボックスはXP画像の全面積の約20%〜60%を覆う。ボックスは正方形である必要はなく、任意の所望の形状にすることができる。各ボックスの平均スペクトルが計算される。最適には、その後、最も明るいボックス(最も高いDN)からスペクトルの約12%、最も暗いボックス(最も低いDN)からスペクトルの30%が廃棄され、残りのスペクトルを平均して、目標スペクトルμtrgを生成する。しかしながら、最も明るいサンプルの約0%〜20%及び最も暗いサンプルの約10%〜50%が廃棄される場合にも、目標スペクトルを使用することができる。
その後、フィルタ値、好ましくはスペクトル角の余弦(白色化された空間内にあることが好ましい)が、好ましくは以下の式を用いて、背景ピクセルA毎に計算される。
このフィルタ値を用いて、ピクセルがフィルタリング除去されるか、保持されるかを判断する。
RGB画像と同じ数のピクセルを有する2次元行列が0で初期化される。背景ピクセル毎のフィルタ値が、この行列に書き込まれる。不良ピクセルの値は0のままにされる。75度の角しきい値が選択される場合には、これは、フィルタ値f(A)=cos(75)=0.259に対応する。0の第2の2次元行列が初期化される。フィルタ値f(A)>0.259を有するピクセルの場合に、その行列は1の値を割り当てられる。上記の結果は2つの行列であり、一方の行列は、背景ピクセル毎にフィルタ値を、そして不良ピクセル毎に0を有し、第2の行列は、十分に高いフィルタ値を有するピクセル毎に1の値を、そして全ての他のピクセルの場合に0を有する。
その後、この行列は、モルフォロジー演算を受けて、子宮頸部ROIがさらに精緻化される。隣接しているピクセルの「ブロブ」又はクラスタが特定される。隣接しているピクセルは、辺を共有するピクセルと定義され、一方、角で接触しているピクセルは隣接しているとは定義されない。各ブロブの面積が計算され、その後、最も大きなブロブを除く全てのブロブが除去される。その結果がROIである。一例が図4として示される。最後に、ROI内のピクセルの数がカウントされる。
3.コントラスト最大化
3.1 問題の記述
コントラストは、情報を与える明るさ(luminance)(輝度(brightness))の階調(gradation)(差(difference))の指標(measure)である。それは、視野の隣接している部分における比(明るさの差)/(平均の明るさ)として表される。一般的に言うと、コントラストは、2つの直ぐ近くの信号が有することができ、且つ依然として区別可能である最も小さな差である。最適な条件下では、人の目は、2%のコントラストの存在を検出することができる。それゆえ、デジタルコルポスコープによって生成される画像において、膣鏡検査熟練者によって観測されるのと同じ臨床的に関連する特徴を見分けるために、そのデバイスは、少なくとも2%のコントラストを検出できなければならない。目標として、そのデバイスは、2%未満のコントラストを検出することができることになる。
コントラストは以下のように定義される。
分子は暗雑音の標準偏差である。これは、デバイス特徴抽出手順中に「暗ファイル」から得られる(その手順は、雑音、ダイナミックレンジ、信号対雑音比、出力電力などに関して機器の特徴を抽出するか、又はその機器を測定する)。分母は、暗ファイルの平均値(基本的にカメラの0信号点)が減算される信号強度である。上記のようにコントラストを定義する正当性は、システムによって検出されることができる明るさ(luminosity)の最小差がシステム雑音(暗雑音から推定される)に等しく、一方、ほとんど同じ振幅の2つの信号の場合の平均の明るさ(luminosity)が、近似的に、それらの信号のうちの一方の信号の振幅になることである。所与の信号の明るさの場合に記録される信号DNが、アナログ/デジタルコンバータの利得を高めることによって増加する場合には、暗雑音も同様に増幅され、結果として、コントラストは改善されない。上述したように、デジタル数(DN)は、デジタル画像内のピクセルの輝度(brightness)を表す正の整数である。コントラスト不良は、光源強度が不十分であるか、又はカメラ露出時間が短すぎる結果である可能性が高い。光源強度を高めると、暗雑音を変更することなく、信号DNが増加する。この場合、式(2)において、分子は一定のままであるが、分母が増加し、それにより、コントラストが下がる(改善される)。しかしながら、デジタルコルポスコープの設計によって、一般的には、光源の強度は一定のままに制限される。結果として、通常は、露出時間を長くすることが、コントラストを改善するために唯一利用可能な方法である。
当然、露出時間が長すぎる場合には、カメラの光子井戸(photon well)、又はそのアナログ/デジタルコンバータをオーバーフローすることによって、子宮頸部ピクセルが飽和した状態になる場合がある。これが発生すると、それらのピクセルのための情報は破壊され、子宮頸部の特有の明るさ(luminance)を測定することはできない。飽和に近づくと、カメラの応答が非線形になるので、実際には、飽和点前に情報が破損する場合がある。それゆえ、最も明るい信号を、飽和よりもわずかに小さな或る最大DNに制限することが望ましい。この値も、デバイス特徴抽出手順から得られるべきである。
グリントピクセルは、光源の鏡面反射のために、その信号DNが異常に(その画像の残りの部分に対して)高いピクセルである。上述したように、鏡面反射は、鏡のような表面からの光の完全な反射であり、単一の入射方向からの光(光線)が単一の出射方向に反射される。PPカメラを用いて得られる偏光によって、少なくともいくつかのグリントピクセルが、ほぼ全ての画像において見られる可能性がある(XPカメラでは、グリントピクセルが抑圧される)。その数がほとんどないものと仮定すると、グリントピクセルは子宮頸部ピクセルの信号DNを表さないので、コントラストを特徴付けるために用いられるべきではない。
コントラスト最大化アルゴリズムの目的は、臨床的に関連する特徴を観測するために収集される画像が、実現し得る最良のコントラストを示すことを確実にすることである。カメラ雑音は基本的には一定であるので、実現し得る最良のコントラストは、最も明るい子宮頸部ピクセルの場合に最も高い値の信号DNを生成するが、それでも非線形範囲又は飽和の開始よりも短い露出時間を、好ましくは所定のリストから決定することによって達成される。
3.2 アルゴリズムの記述
コントラスト最大化アルゴリズムは2つの主要部分を有する。第1の部分は信号評価である。信号評価の目的は、最も明るい子宮頸部ピクセルの信号DNを推定することである。第2の部分は露出時間調整であり、それは、必要に応じて、信号DNの値に基づいて露出時間を調整する。コントラスト最大化アルゴリズムは、1つの画像内で測定されるデジタル数範囲をさらに広い範囲に拡張することを伴うコントラスト伸張(正規化と呼ばれることが多い)を利用して、画像及びその構成部分のコントラストを改善する。「伸張」は、より高い(より飽和点に近い)DN値において、より多くの数のピクセルを得るように露出時間を調整することによって達成される。
子宮頸部の反射率は、赤色に対応する光学スペクトル範囲において最も高い。結果として、子宮頸部のための信号DNは常に、画像の赤色帯において最も高い。それゆえ、そのアルゴリズムは、赤色帯のみを用いる(たとえば、青色帯内の信号DNが最大にされたなら、赤色帯はおそらく飽和するであろう)。その後、このアルゴリズムは、対象領域内にない、赤色帯内の全てのピクセルを除去する。その後、残りのピクセルのヒストグラムが計算される。当然、他の用途の場合には、異なるスペクトル帯を用いることができる。
このアルゴリズムは、入力として、信号DNの最大許容値及び最小許容値を受け入れる。XP画像及びPP画像の双方が評価される。好ましくは、最小許容値は、その画像の最大DN値の約80%〜90%に相当することが好ましい。8ビット画像の場合、DN値は0〜255の範囲内にある。本発明の場合、0〜1023の範囲内のDN値を有する10ビット画像を用いる。
最大許容値において開始して、このアルゴリズムは、ヒストグラムの形状を解析し、「エッジ」を探索する。エッジでは、或る特定のDN範囲内のピクセルの数が急激に増加し始める。PP画像の一例が、図5(a)として与えられる。
図5(b)において、そのヒストグラムによれば、この画像(DN=1023)にはほぼ200の飽和ピクセルがある。しかしながら、これらのピクセルは、ヒストグラムエッジを探索する際に無視される。飽和値(又は「飽和ビン」)の左隣のヒストグラム値(又は「ヒストグラムビン」)は、非常に少ない数のピクセルを有する。その後、ヒストグラムは、上に向かって曲がり始め、最も明るい子宮頸部ピクセルを特徴付ける。この画像の最大信号DN、860は、800の最小許容値と920の最大許容値との中間にある。ここでも、最小許容値及び最大許容値は、この画像に利用することができるDNの最大値(1023)の約80%(800)、及びこの画像に利用することができるDNの最大値(1023)の約90%(920)に相当する。結果として、図5(a)の画像は満足のいくコントラストを有し、露出時間の調整は不要である。
露出過度の画像の一例が図6(a)として与えられる。飽和したピクセルの数(DN=1023)は多く、ヒストグラム上の任意の他の値よりも多い。図6(b)を参照されたい。しかしながら、ヒストグラムエッジの探索は、最大許容DN値において開始する。このアルゴリズムは、そのヒストグラムが平坦域にあり、エッジにはないことを直ちに検出する。飽和したピクセルの場合、信号DNのために、最大許容値として入力される値が報告されるが、これは真の値ではない。なぜなら、「真の」値は飽和に起因してわからないためである。この場合、コントラストアルゴリズムは、その画像が露出過度であり、不満足なコントラストを有することを報告する。その後、この情報はコントラスト最大化アルゴリズムの第2のステージに渡される。最も明るい子宮頸部ピクセルを特徴付ける値はわからないので、最適な露出時間の推定値を定量化することはできない。露出時間は、離散した1組の許容値に制限されることが好ましい。このため、画像が露出過度である場合には、露出時間は約5%〜10%だけ短くされることが好ましく、新たな画像が収集される。この新たな露出時間も長すぎる場合には、そのアルゴリズムは再び、約5%〜10%だけ露出時間を短くする。
最後に、露出不足の画像の一例が図7(a)として与えられる。ここでも、多数の飽和したピクセルが存在する(DN=1023)。図7(b)を参照されたい。しかしながら、ほぼ空のヒストグラム値が、飽和値の左側に長く延在する。そのヒストグラムは、約DN=750まで、上に向かって曲がり始めない。これは、800の最小許容値未満である。結果として、この画像のコントラストは不満足であり、露出不足と呼ばれる。
露出不足の場合、最適な露出時間を推定することができる。信号DNは、露出時間とともに線形に増減する。新たな露出時間は式(3)を用いて推定される。
ここで、860は信号DNの目標値であり、それは最小許容値と最大許容値との平均として求められる。適当に増減させるために、0信号に対応するDN値が、目標値及び現在値から減算される。その後、そのアルゴリズムは、リストにおいて、新たな推定値に最も近い露出時間を選択する。
露出不足の場合には、新たな露出時間を定量可能に推定することができるので、露出過度よりも露出不足の方が好ましい。それゆえ、相対的に短い露出時間で評価アルゴリズムを開始することが好ましい。コントラスト最大化アルゴリズムは、時間を短くする方に1つずつ戻らなければならない場合よりも、迅速に正確な露出時間を求めることができる。
4.構図解析
4.1 問題の記述
上述したように、子宮頸部は各画像内の中心に配置されるべきであり、画像境界と交差することなく、できる限り視野の多くを占めるべきである。構図アルゴリズムの目的は、これが当てはまるか否かを評価し、当てはまらない場合にはオペレータにフィードバックを与えることである。
4.2 アルゴリズムの記述
XP画像からのROIが、このアルゴリズムの入力として用いられる。最初に、ROI境界が画像境界に近いかが評価される。好ましくは、ROI境界は、画像境界と交差することなく、画像境界に極めて近くすべきである。これを表す1つの方法は、全ROI面積を全画像面積の約90%〜99%にし、95%を最適値にすることである。しかしながら、本発明は、ROIの面積が全画像面積の80%〜99%の範囲内にある場合にも実施することができる。こうして、ROIの全面積が画像面積の約90%未満であるか、又は画像面積の100%よりも大きい場合には、BooleanBadBorders品質フラグが真に設定される。次に、ROIの中心が、最大行位置及び最小行位置、並びに最大列位置及び最小列位置の中点であると見なされ、画像中心からのROI中心の距離が計算される。ROI中点と画像中心との間の距離の差は、5%の最大値しきい値を有し、それを超える場合には、BooleanBadCenter品質フラグが真に設定される。最後に、カメラ視野内で子宮頸部が遮られる場合を検出するために、ROIの大きさが、画像の全面積の少なくとも約25%の最小しきい値と比較される。ROIの大きさがこの最小しきい値よりも小さい場合には、BooleanBadArea品質フラグが真に設定される。上記の品質フラグのいずれかが真である場合には、その構図は不良であり、テスト(複数可)は失敗である。しかしながら、ROIが中心にあり、画像境界と交差せず、且つ十分に大きい場合には、出力構図診断メッセージは、その構図が良好であることを指示する。
5.ステータス評価
ステータス評価は別のブール変数である。それは、XP及びPPのコントラスト最大化及びXPの場合の構図がいずれも満たされる場合には、良好(真)の値を割り当てられる。その際、そのデバイスは、臨床的に関連する特徴を観測するために高解像度のXP画像及びPP画像を収集できる状態にある。
ステータス評価が偽である場合には、或る構成要素が失敗しており、露出時間メッセージ又は構図品質フラグにおいて推奨される調整が利用されるべきである。これらの調整が行なわれた後に、評価アルゴリズムが繰り返されるべきである。
詳細な説明において記述され、図面において例示される実施形態を参照しながら本発明を具体的に図示及び説明してきたが、特許請求の範囲によって規定されるような、本発明の精神及び範囲から逸脱することなく、細部において種々の変更が加えられる場合があることは、当業者には理解されよう。したがって、特許請求の範囲において明示される場合を除いて、限定が暗示、又は推論されるべきではない。
本発明は、組織診断システムが、捕捉された画像の品質を自動的に解析し、調整するための方法を提供する。それは、背景とは異なるスペクトルを有する対象領域の任意の撮像及び解析に適用することができる。

Claims (14)

  1. 電子診療記録及びコンピュータ支援診断(CAD)システムにおいて用いるためのアーカイブ品質標準化画像データを達成するために、撮像デバイスに自動品質フィードバックを与える方法であって、
    デジタル撮像装置を用いて調査中のライブビデオから2つのバージョンの生画像を収集することであって、該生画像の第1のバージョンは第1の偏光方向を有する光源を用いて視野を照明し、検出器を用いて第2の偏光方向を有する偏光フィルタを通じて該視野を視認することによって作成されて、該視野の交差偏光画像が生成され、前記第1の偏光方向は前記第2の偏光方向に対して実質的に垂直であり、前記生画像の第2のバージョンは、前記第2の偏光方向が前記第1の偏光方向に対して実質的に平行であるときに生成されて、前記視野の平行偏光画像が生成される、収集することであり、ここで、前記視野の前記画像はいずれも実質的に互いに重ね合わせられ、それにより該画像は実質的に同じ境界であって、該境界はそれぞれ、該境界内の面積を画定する境界長及び境界幅を有する、境界、及び1つの中心を有するものであり、及び、
    画像分類アルゴリズムを用いて前記交差偏光画像内の対象領域の位置を特定することであって、該画像分類アルゴリズムは、前記交差偏光画像の空間的に離散したサンプルの輝度平均スペクトルを求めること、最も明るいサンプルの約0%〜20%及び最も暗いサンプルの約10%〜50%を廃棄すること、残りのサンプルのスペクトルを平均して、輝度目標スペクトルを求めることによって、前記交差偏光画像から該目標スペクトルを規定すること、及び、該目標スペクトルと実質的に一致する前記交差偏光画像の部分を特定することによって該対象領域の大きさを求めることを含む、位置を特定すること、
    前記交差偏光画像及び前記平行偏光画像内の前記対象領域に対する信号評価及び露出時間調整によってコントラスト最大化を実行すること、
    前記交差偏光画像内の前記対象領域に構図解析を適用することであって、該対象領域が前記交差偏光画像内の中心に位置するのを確実にする、適用すること、
    前記コントラスト最大化及び前記構図解析の評価ステータスを判断することであって、該評価ステータスは良好又は不良いずれかのブール変数として報告される、判断すること、
    前記ステータスが不良である場合には、前記2つのバージョンの生画像を収集すること、前記対象領域の位置を特定すること、前記コントラスト最大化及び前記構図解析を実行することを繰り返すこと、並びに、
    前記ステータスが良好である場合には、高解像度交差偏光静止画像及び高解像度平行偏光静止画像を収集することを含む、方法。
  2. 電子診療記録及びコンピュータ支援診断(CAD)システムにおいて用いるためのアーカイブ品質標準化画像データを達成するために、撮像デバイスに自動品質フィードバックを与える方法であって、
    デジタル撮像装置を用いて調査中のライブビデオから2つのバージョンの生画像を収集することであって、該生画像の第1のバージョンは第1の偏光方向を有する光源を用いて視野を照明し、検出器を用いて第2の偏光方向を有する偏光フィルタを通じて該視野を視認することによって作成されて、該視野の交差偏光画像が生成され、前記第1の偏光方向は前記第2の偏光方向に対して実質的に垂直であり、前記生画像の第2のバージョンは、前記光源及び前記検出器から成る群から選択されるロケーションに偏光フィルタを配置することによって得られ、該視野の単一偏光画像が生成される、収集することであり、ここで、前記視野の前記画像はいずれも実質的に互いに重ね合わせられ、それにより該画像は実質的に同じ境界であって、該境界はそれぞれ、該境界内の面積を画定する境界長及び境界幅を有する、境界、及び1つの中心を有するものであり、及び、
    画像分類アルゴリズムを用いて前記交差偏光画像の対象領域の位置を特定することであって、該画像分類アルゴリズムは、前記交差偏光画像の空間的に離散したサンプルの輝度平均スペクトルを求めること、最も明るいサンプルの約0%〜20%及び最も暗いサンプルの約10%〜50%を廃棄すること、残りのサンプルのスペクトルを平均して、輝度目標スペクトルを求めることによって、前記交差偏光画像から該目標スペクトルを規定すること、及び、該目標スペクトルと実質的に一致する前記交差偏光画像の部分を特定することによって該対象領域の大きさを求めることを含む、位置を特定すること、
    前記交差偏光画像及び前記単一偏光画像内の前記対象領域に対する信号評価及び露出時間調整によってコントラスト最大化を実行すること、
    前記交差偏光画像内の前記対象領域に構図解析を適用することであって、該対象領域が前記交差偏光画像内の中心に位置するのを確実にする、適用すること、
    前記コントラスト最大化及び前記構図解析の評価ステータスを判断することであって、該評価ステータスは良好又は不良いずれかのブール変数として報告される、判断すること、
    前記状態が不良である場合には、前記2つのバージョンの生画像を収集すること、前記対象領域の位置を特定すること、前記コントラスト最大化及び前記構図解析を実行することを繰り返すこと、並びに、
    前記ステータスが良好である場合には、高解像度交差偏光静止画像及び高解像度単一偏光静止画像を収集することを含む、方法。
  3. 前記画像分類アルゴリズムは、スペクトル整合フィルタを用いて前記対象領域の大きさを求めて、残りのサンプルの前記スペクトルと前記目標スペクトルとの間のスペクトル角を推定し、前記交差偏光画像は、前記スペクトル角が約65度〜85度未満であるときに、前記目標スペクトルと実質的に一致する、請求項1又は2に記載の方法。
  4. 前記画像分類アルゴリズムは、前記交差偏光画像の前記中心付近に分散する複数のボックスを含む空間的に離散したスペクトルサンプルを用いて実行され、前記ボックスはそれぞれ、ボックス長及びボックス幅を有し、前記ボックス長は前記交差偏光画像の前記境界長の約10%〜30%であり、前記ボックスは、前記交差偏光画像の前記境界内にある前記面積の約20%〜60%を覆う、請求項1又は2に記載の方法。
  5. 前記画像分類アルゴリズムは、最も明るいサンプルのスペクトルの約10〜15%、及び最も暗いサンプルのスペクトルの約20〜40%を廃棄することによって実行される、請求項1又は2に記載の方法。
  6. 前記画像分類アルゴリズムは、最も明るいサンプルのスペクトルの12%、及び最も暗いサンプルのスペクトルの30%を廃棄することによって実行される、請求項1若しくは2に記載の方法。
  7. 前記コントラスト最大化ステップは、前記対象領域の露出時間を調整して、その輝度(brightness)がデジタル数によって決定される、前記交差偏光画像及び前記平行偏光画像内の最も明るいピクセルを生成することによって実行され、前記最も明るいピクセルは、満足のいくコントラストを生成するために、所望の範囲に入るデジタル数を有し、該所望の範囲の最小値は前記交差偏光画像の最大デジタル数値の約80%であり、且つ前記平行偏光画像の最大デジタル数値の80%であり、前記所望の範囲の最大値は、前記交差偏光画像の最大デジタル数値の約90%に対応し、且つ、前記平行偏光画像の最大デジタル数値の約90%であることが好ましい、請求項1に記載の方法。
  8. 前記コントラスト最大化ステップは、前記対象領域の露出時間を調整して、その輝度(brightness)がデジタル数によって決定される、前記交差偏光画像及び前記単一偏光画像内の最も明るいピクセルを生成することによって実行され、前記最も明るいピクセルは、満足のいくコントラストを生成するために、所望の範囲に入るデジタル数を有し、該所望の範囲の最小値は前記交差偏光画像の最大デジタル数値の約80%であり、且つ前記単一偏光画像の最大デジタル数値の80%であり、前記所望の範囲の最大値は、前記交差偏光画像の最大デジタル数値の約90%に対応し、且つ、前記単一偏光画像の最大デジタル数値の約90%であることが好ましい、請求項2に記載の方法。
  9. 前記構図解析ステップは、前記交差偏光画像内で、前記対象領域の境界が前記交差偏光画像の前記境界が近いかを評価し、前記対象領域の中心が前記交差偏光画像の前記中心に近いかを判断し、前記対象領域の大きさを計算し、構図診断メッセージを報告することによって実行される、請求項1又は2に記載の方法。
  10. 前記構図診断メッセージは、前記対象領域が前記交差偏光画像内の中心に位置し、前記交差偏光画像境界内に包含され、前記対象領域の前記大きさが前記交差偏光画像の全面積の少なくとも約25%であるときに、満足のいくものとして報告される、請求項9に記載の方法。
  11. 前記コントラスト最大化及び前記構図解析が満足であるときにのみ、良好の前記ブール変数が当てはまる、請求項1又は2に記載の方法。
  12. 前記コントラスト最大化若しくは前記構図解析のいずれかが不満足であるときに、又は前記コントラスト最大化及び前記構図解析の双方が不満足であるときに、不良の前記ブール変数が当てはまる、請求項1又は2に記載の方法。
  13. 電子診療記録及びコンピュータ支援診断(CAD)システムにおいて用いるためのアーカイブ品質標準化画像データを達成するために、撮像デバイスに自動品質フィードバックを与える方法であって、
    デジタル撮像装置を用いて調査中のライブビデオから2つのバージョンの生画像を収集することであって、該生画像の第1のバージョンは第1の偏光方向を有する光源を用いて視野を照明し、検出器を用いて第2の偏光方向を有する偏光フィルタを通じて該視野を視認することによって作成され、該視野の交差偏光画像が生成され、前記第1の偏光方向は前記第2の偏光方向に対して実質的に垂直であり、前記生画像の第2のバージョンは、前記第2の偏光方向が前記第1の偏光方向に対して実質的に平行であるときに生成されて、前記視野の平行偏光画像が生成される、収集することであり、ここで、前記視野の前記画像はいずれも実質的に互いに重ね合わせられ、それにより該画像は実質的に同じ境界であって、該境界はそれぞれ、該境界内の面積を画定する境界長及び境界幅を有する、境界、及び1つの中心を有するものであり、及び、
    画像分類アルゴリズムを用いて前記交差偏光画像内の対象領域の位置を特定(locating)することであって、該画像分類アルゴリズムは、前記交差偏光画像の空間的に離散したサンプルの輝度平均スペクトルを求めること、最も明るいサンプルの約0%〜20%及び最も暗いサンプルの約10%〜50%を廃棄すること、残りのサンプルのスペクトルを平均して、輝度目標スペクトルを求めることによって、前記交差偏光画像から該目標スペクトルを規定すること、及び、該目標スペクトルと実質的に一致する前記交差偏光画像の部分を特定することによって該対象領域の大きさを求めることを含む、位置を特定すること、
    前記対象領域の露出時間を調整して、その輝度(brightness)がデジタル数によって決定される、前記交差偏光画像及び前記平行偏光画像内の最も明るいピクセルを生成することによって、前記交差偏光画像及び前記平行偏光画像に対する信号評価及び露出時間調整によりコントラスト最大化を実行することであって、該最も明るいピクセルは、満足のいくコントラストを生成するために、所望の範囲に入るデジタル数を有し、該所望の範囲の最小値は前記交差偏光画像の最大デジタル数値の約80%であり、且つ前記平行偏光画像の最大デジタル数値の80%であり、該所望の範囲の最大値は、前記交差偏光画像の最大デジタル数値の約90%に対応し、且つ前記平行偏光画像の最大デジタル数値の約90%であることが好ましい、実行すること、
    前記交差偏光画像内で、前記対象領域の境界が前記交差偏光画像の前記境界が近いかを評価し、前記対象領域の中心が前記交差偏光画像の前記中心に近いかを判断し、前記対象領域の大きさを計算し、構図診断メッセージを報告することによって、前記交差偏光画像に構図解析を適用することであって、該構図診断メッセージは、前記対象領域が前記交差偏光画像内の中心に位置し、前記交差偏光画像境界内に包含され、前記対象領域の該大きさが前記交差偏光画像の全面積の少なくとも約25%であるときに、満足のいくものとして報告される、適用すること、
    前記コントラスト最大化及び前記構図解析の評価ステータスを判断することであって、該評価ステータスは良好又は不良いずれかのブール変数として報告され、該良好のブール変数は、前記コントラスト最大化及び前記構図解析が満足のいくときにのみ当てはまり、該不良のブール変数は、前記コントラスト最大化若しくは前記構図解析のいずれかが不満足であるときに、又は前記コントラスト最大化及び前記構図解析の双方が不満足であるときに当てはまる、判断すること、
    前記ステータスが不良である場合には、前記2つのバージョンの生画像を収集すること、前記対象領域の位置を特定すること、前記コントラスト最大化及び前記構図解析を実行することを繰り返すこと、並びに、
    前記ステータスが良好である場合には、高解像度交差偏光静止画像及び高解像度平行偏光静止画像を収集することを含む、方法。
  14. 電子診療記録及びコンピュータ支援診断(CAD)システムにおいて用いるためのアーカイブ品質標準化画像データを達成するために、撮像デバイスに自動品質フィードバックを与える方法であって、
    デジタル撮像装置を用いて調査中のライブビデオから2つのバージョンの生画像を収集することであって、該生画像の第1のバージョンは第1の偏光方向を有する光源を用いて視野を照明し、検出器を用いて第2の偏光方向を有する偏光フィルタを通じて該視野を視認することによって作成され、該視野の交差偏光画像が生成され、前記第1の偏光方向は該第2の偏光方向に対して実質的に垂直であり、前記生画像の第2のバージョンは、前記光源及び前記検出器から成る群から選択されるロケーションに偏光フィルタを配置することによって得られ、前記視野の単一偏光画像が生成される、収集することであり、ここで、前記視野の該画像はいずれも実質的に互いに重ね合わせられ、それにより該画像は実質的に同じ境界であって、該境界はそれぞれ、該境界内の面積を画定する境界長及び境界幅を有する、境界、及び1つの中心を有するものであり、及び、
    画像分類アルゴリズムを用いて前記交差偏光画像内の対象領域の位置を特定することであって、該画像分類アルゴリズムは、前記交差偏光画像の空間的に離散したサンプルの輝度平均スペクトルを求めること、最も明るいサンプルの約0%〜20%及び最も暗いサンプルの約10%〜50%を廃棄すること、残りのサンプルのスペクトルを平均して、輝度目標スペクトルを求めることによって、前記交差偏光画像から該目標スペクトルを規定すること、及び、該目標スペクトルと実質的に一致する前記交差偏光画像の部分を特定することによって該対象領域の大きさを求めることを含む、位置を特定すること、
    前記対象領域の露出時間を調整して、その輝度(brightness)がデジタル数によって決定される、前記交差偏光画像及び前記単一偏光画像内の最も明るいピクセルを生成することによって、前記交差偏光画像及び前記単一偏光画像に対する信号評価及び露出時間調整によりコントラスト最大化を実行することであって、該最も明るいピクセルは、満足のいくコントラストを生成するために、所望の範囲に入るデジタル数を有し、該所望の範囲の最小値は前記交差偏光画像の最大デジタル数値の約80%であり、且つ前記単一偏光画像の最大デジタル数値の80%であり、該所望の範囲の最大値は、前記交差偏光画像の最大デジタル数値の約90%に対応し、且つ、前記単一偏光画像の最大デジタル数値の約90%であることが好ましい、実行すること、
    前記交差偏光画像内で、前記対象領域の境界が前記交差偏光画像の前記境界が近いかを評価し、前記対象領域の中心が前記交差偏光画像の前記中心に近いかを判断し、前記対象領域の大きさを計算し、構図診断メッセージを報告することによって、前記交差偏光画像に構図解析を適用することであって、該構図診断メッセージは、前記対象領域が前記交差偏光画像内の中心に位置し、前記交差偏光画像境界内に包含され、前記対象領域の該大きさが前記交差偏光画像の全面積の少なくとも約25%であるときに、満足のいくものとして報告される、適用すること、
    前記コントラスト最大化及び前記構図解析の評価ステータスを判断することであって、該評価ステータスは良好又は不良いずれかのブール変数として報告され、該良好のブール変数は、前記コントラスト最大化及び前記構図解析が満足のいくときにのみ当てはまり、該不良のブール変数は、前記コントラスト最大化若しくは前記構図解析のいずれかが不満足であるときに、又は前記コントラスト最大化及び前記構図解析の双方が不満足であるときに当てはまる、判断すること、
    前記ステータスが不良である場合には、前記2つのバージョンの生画像を収集すること、前記対象領域の位置を特定すること、前記コントラスト最大化及び前記構図解析を実行することを繰り返すこと、並びに、
    前記ステータスが良好である場合には、高解像度交差偏光静止画像及び高解像度単一偏光静止画像を収集することを含む、方法。
JP2009554535A 2007-03-16 2008-03-14 標準化された撮像データを実現するために撮像デバイスに自動品質フィードバックを与える方法 Pending JP2010521272A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US91852707P 2007-03-16 2007-03-16
PCT/US2008/003378 WO2008115410A2 (en) 2007-03-16 2008-03-14 A method to provide automated quality feedback to imaging devices to achieve standardized imaging data

Publications (1)

Publication Number Publication Date
JP2010521272A true JP2010521272A (ja) 2010-06-24

Family

ID=39762742

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009554535A Pending JP2010521272A (ja) 2007-03-16 2008-03-14 標準化された撮像データを実現するために撮像デバイスに自動品質フィードバックを与える方法

Country Status (4)

Country Link
US (2) US8401258B2 (ja)
EP (1) EP2137696A2 (ja)
JP (1) JP2010521272A (ja)
WO (2) WO2008115410A2 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013066166A (ja) * 2011-08-30 2013-04-11 Ricoh Co Ltd 撮像装置、並びに、これを用いる画像解析装置及び移動装置
WO2022014258A1 (ja) * 2020-07-17 2022-01-20 富士フイルム株式会社 プロセッサ装置、プロセッサ装置の作動方法
JP7517325B2 (ja) 2019-03-25 2024-07-17 ソニーグループ株式会社 医療システム、信号処理装置、及び、信号処理方法
JP7547856B2 (ja) 2020-08-21 2024-09-10 カシオ計算機株式会社 撮像制御装置、撮像装置、撮像制御方法及びプログラム

Families Citing this family (111)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ES2691646T3 (es) 2005-07-15 2018-11-28 Micell Technologies, Inc. Revestimientos poliméricos que contienen polvo de fármaco de morfología controlada
WO2007011708A2 (en) 2005-07-15 2007-01-25 Micell Technologies, Inc. Stent with polymer coating containing amorphous rapamycin
EP2412301B1 (en) * 2005-12-28 2014-04-23 Olympus Medical Systems Corp. Image processing device and image processing method in image processing device
AU2007243268B2 (en) 2006-04-26 2013-06-13 Micell Technologies, Inc. Coatings containing multiple drugs
US11426494B2 (en) 2007-01-08 2022-08-30 MT Acquisition Holdings LLC Stents having biodegradable layers
EP2111184B1 (en) 2007-01-08 2018-07-25 Micell Technologies, Inc. Stents having biodegradable layers
US8086007B2 (en) * 2007-10-18 2011-12-27 Siemens Aktiengesellschaft Method and system for human vision model guided medical image quality assessment
BRPI0910969B8 (pt) 2008-04-17 2021-06-22 Micell Technologies Inc dispositivo
US9486431B2 (en) 2008-07-17 2016-11-08 Micell Technologies, Inc. Drug delivery medical device
US8270695B2 (en) * 2008-10-07 2012-09-18 Carestream Health, Inc. Diagnostic image processing with automatic self image quality validation
US8483454B2 (en) * 2008-10-10 2013-07-09 Sti Medical Systems, Llc Methods for tissue classification in cervical imagery
CN101448174A (zh) * 2008-12-26 2009-06-03 深圳华为通信技术有限公司 一种图像质量评价装置和方法
US8290295B2 (en) * 2009-03-03 2012-10-16 Microsoft Corporation Multi-modal tone-mapping of images
JP2012522589A (ja) 2009-04-01 2012-09-27 ミシェル テクノロジーズ,インコーポレイテッド 被覆ステント
EP2419058B1 (en) 2009-04-17 2018-02-28 Micell Technologies, Inc. Stents having controlled elution
EP2453834A4 (en) 2009-07-16 2014-04-16 Micell Technologies Inc MEDICAL DEVICE DISPENSING MEDICINE
US8825131B2 (en) 2009-10-14 2014-09-02 Nocimed, Llc MR spectroscopy system and method for diagnosing painful and non-painful intervertebral discs
US8761860B2 (en) 2009-10-14 2014-06-24 Nocimed, Llc MR spectroscopy system and method for diagnosing painful and non-painful intervertebral discs
US20110255762A1 (en) * 2010-04-15 2011-10-20 Harald Deischinger Method and system for determining a region of interest in ultrasound data
EP2560576B1 (en) 2010-04-22 2018-07-18 Micell Technologies, Inc. Stents and other devices having extracellular matrix coating
WO2011139362A1 (en) * 2010-05-03 2011-11-10 Sti Medical Systems, Llc Image analysis for cervical neoplasia detection and diagnosis
US8805112B2 (en) 2010-05-06 2014-08-12 Nikon Corporation Image sharpness classification system
WO2011151821A1 (en) * 2010-05-31 2011-12-08 Dvp Technologies Ltd. Inspection of region of interest
US8666148B2 (en) 2010-06-03 2014-03-04 Adobe Systems Incorporated Image adjustment
JP5761601B2 (ja) * 2010-07-01 2015-08-12 株式会社リコー 物体識別装置
CA2805631C (en) 2010-07-16 2018-07-31 Micell Technologies, Inc. Drug delivery medical device
US8773577B2 (en) 2010-10-27 2014-07-08 Qualcomm Incorporated Region of interest extraction
JP5725194B2 (ja) 2010-11-03 2015-05-27 株式会社ニコン 夜景画像ボケ検出システム
US9280718B2 (en) * 2010-11-24 2016-03-08 Nocimed, Llc Systems and methods for automated voxelation of regions of interest for magnetic resonance spectroscopy
FR2968811A1 (fr) * 2010-12-13 2012-06-15 Univ Paris Diderot Paris 7 Methode de detection et de quantification du flou dans une image numerique.
IT1404810B1 (it) * 2011-01-26 2013-11-29 St Microelectronics Srl Riconoscimento di texture nell'elaborazione di immagini
JP4831259B1 (ja) * 2011-03-10 2011-12-07 オムロン株式会社 画像処理装置、画像処理方法、および制御プログラム
US8798342B2 (en) * 2011-05-10 2014-08-05 General Electric Company Method and system for ultrasound imaging with cross-plane images
US8854041B2 (en) * 2011-05-20 2014-10-07 Kabushiki Kaisha Toshiba Spatially shaped pre-saturation profile for enhanced non-contrast MRA
WO2013012689A1 (en) 2011-07-15 2013-01-24 Micell Technologies, Inc. Drug delivery medical device
US9251439B2 (en) 2011-08-18 2016-02-02 Nikon Corporation Image sharpness classification system
CN102955947B (zh) * 2011-08-19 2017-09-22 北京百度网讯科技有限公司 一种用于确定图像清晰度的设备及其方法
US8903169B1 (en) 2011-09-02 2014-12-02 Adobe Systems Incorporated Automatic adaptation to image processing pipeline
US9008415B2 (en) * 2011-09-02 2015-04-14 Adobe Systems Incorporated Automatic image adjustment parameter correction
CN103020947B (zh) 2011-09-23 2016-04-06 阿里巴巴集团控股有限公司 一种图像的质量分析方法及装置
CA2851152A1 (en) * 2011-10-05 2013-04-11 Cireca Theranostics, Llc Method and system for analyzing biological specimens by spectral imaging
US10188772B2 (en) 2011-10-18 2019-01-29 Micell Technologies, Inc. Drug delivery medical device
DE102011090047A1 (de) 2011-12-28 2013-07-25 Klinikum der Universität München - Campus Innenstadt Kontrollverfahren und Kontrollsystem
WO2013126568A1 (en) 2012-02-21 2013-08-29 Massachusetts Eye & Ear Infirmary Calculating conjunctival redness
US8965094B2 (en) 2012-04-14 2015-02-24 Nocimed, Llc Magnetic resonance spectroscopy pulse sequence, acquisition, and processing system and method
CN102708568B (zh) * 2012-05-11 2014-11-05 宁波大学 一种基于结构失真的立体图像客观质量评价方法
CN102708567B (zh) * 2012-05-11 2014-12-10 宁波大学 一种基于视觉感知的立体图像质量客观评价方法
JP5789341B2 (ja) * 2012-09-18 2015-10-07 富士フイルム株式会社 静止画表示装置及びシステム並びに撮像装置
US9798918B2 (en) * 2012-10-05 2017-10-24 Cireca Theranostics, Llc Method and system for analyzing biological specimens by spectral imaging
EP2920965A1 (en) * 2012-11-16 2015-09-23 VID SCALE, Inc. Perceptual preprocessing filter for viewing-conditions-aware video coding
US9936858B2 (en) * 2013-02-04 2018-04-10 Orpheus Medical Ltd Color reduction in images of an interior of a human body
WO2014165264A1 (en) 2013-03-12 2014-10-09 Micell Technologies, Inc. Bioabsorbable biomedical implants
CN105307554B (zh) 2013-04-18 2018-01-09 皇家飞利浦有限公司 采集宫颈图像
WO2014186532A1 (en) 2013-05-15 2014-11-20 Micell Technologies, Inc. Bioabsorbable biomedical implants
CN103325128B (zh) * 2013-05-16 2016-08-17 深圳市理邦精密仪器股份有限公司 一种智能识别阴道镜所采集的图像特征的方法及装置
KR20140137715A (ko) * 2013-05-23 2014-12-03 삼성디스플레이 주식회사 엑스레이 검출 장치 및 방법
RU2016119373A (ru) * 2013-10-22 2017-11-28 Конинклейке Филипс Н.В. Визуализация изображений
US20150244946A1 (en) * 2013-11-04 2015-08-27 Sos Agaian Method and systems for thermal image / video measurements and processing
EP2887266A1 (en) * 2013-12-23 2015-06-24 Koninklijke Philips N.V. Automatic extraction of a region of interest based on speculum detection
US10045050B2 (en) * 2014-04-25 2018-08-07 Vid Scale, Inc. Perceptual preprocessing filter for viewing-conditions-aware video coding
AU2015252852B2 (en) 2014-05-02 2019-03-28 Massachusetts Eye And Ear Infirmary Grading corneal fluorescein staining
US10026010B2 (en) * 2014-05-14 2018-07-17 At&T Intellectual Property I, L.P. Image quality estimation using a reference image portion
CN104504676B (zh) * 2014-11-07 2017-03-08 嘉兴学院 一种多视觉敏感特征相似度的全参考图像质量评价方法
US10516865B2 (en) * 2015-05-17 2019-12-24 Endochoice, Inc. Endoscopic image enhancement using contrast limited adaptive histogram equalization (CLAHE) implemented in a processor
US9832484B2 (en) * 2015-05-20 2017-11-28 Texas Instruments Incorporated Still block detection in a video sequence
US9741107B2 (en) * 2015-06-05 2017-08-22 Sony Corporation Full reference image quality assessment based on convolutional neural network
US10360675B2 (en) 2015-06-12 2019-07-23 International Business Machines Corporation Methods and systems for automatically analyzing clinical images using rules and image analytics
US10289940B2 (en) * 2015-06-26 2019-05-14 Here Global B.V. Method and apparatus for providing classification of quality characteristics of images
US10043075B2 (en) * 2015-11-19 2018-08-07 Microsoft Technology Licensing, Llc Eye feature identification
US10402967B2 (en) * 2015-12-21 2019-09-03 Koninklijke Philips N.V. Device, system and method for quality assessment of medical images
CN106971386B (zh) * 2016-01-14 2020-03-20 广州市动景计算机科技有限公司 判断图像完整度及页面加载度的方法、装置和客户端设备
CN106022354B (zh) * 2016-05-07 2018-09-18 浙江大学 基于svm的图像mtf测量方法
CN107451959B (zh) * 2016-05-31 2020-10-16 宇龙计算机通信科技(深圳)有限公司 图片处理方法及系统
AU2017282665B2 (en) 2016-06-19 2022-05-19 Aclarion, Inc. Magnetic resonance spectroscopy system and method for diagnosing pain or infection associated with propionic acid
KR101789513B1 (ko) * 2016-07-11 2017-10-26 주식회사 인피니트헬스케어 디지털 병리 시스템의 이미지 품질 평가 방법
CN106210446B (zh) * 2016-07-18 2019-08-30 北京小米移动软件有限公司 饱和度增强方法及装置
US10482633B2 (en) * 2016-09-12 2019-11-19 Zebra Medical Vision Ltd. Systems and methods for automated detection of an indication of malignancy in a mammographic image
JP6801434B2 (ja) * 2016-12-20 2020-12-16 富士通株式会社 生体画像処理装置、生体画像処理方法および生体画像処理プログラム
CN107203991B (zh) * 2017-04-19 2020-05-19 山西农业大学 一种基于光谱残差的半参考图像质量评价方法
US10713537B2 (en) 2017-07-01 2020-07-14 Algolux Inc. Method and apparatus for joint image processing and perception
EP3699865B1 (en) * 2017-10-20 2022-11-30 Nec Corporation Three-dimensional face shape derivation device, three-dimensional face shape deriving method, and non-transitory computer readable medium
EP3489892B1 (en) * 2017-11-24 2022-01-05 Ficosa Adas, S.L.U. Determining clean or dirty captured images
US10832808B2 (en) 2017-12-13 2020-11-10 International Business Machines Corporation Automated selection, arrangement, and processing of key images
CN108230314B (zh) * 2018-01-03 2022-01-28 天津师范大学 一种基于深度激活池化的图像质量评估方法
KR20200116107A (ko) * 2018-01-23 2020-10-08 모바일오디티 엘티디. 의료적 이미지화 절차의 자동화된 모니터링
DE102018201794B3 (de) * 2018-02-06 2019-04-11 Heidelberger Druckmaschinen Ag Adaptive Bildglättung
CN108830823B (zh) * 2018-03-14 2021-10-26 西安理工大学 基于空域结合频域分析的全参考型图像质量评价方法
CN108447059B (zh) * 2018-04-09 2021-06-29 华侨大学 一种全参考光场图像质量评价方法
US11288800B1 (en) * 2018-08-24 2022-03-29 Google Llc Attribution methodologies for neural networks designed for computer-aided diagnostic processes
CN109473149B (zh) * 2018-11-09 2021-01-15 天津开心生活科技有限公司 数据质量评估方法、装置、电子设备及计算机可读介质
US11024023B2 (en) * 2018-11-19 2021-06-01 Vision Guided Robotics LLC Inspection system
US11436720B2 (en) * 2018-12-28 2022-09-06 Shanghai United Imaging Intelligence Co., Ltd. Systems and methods for generating image metric
EP3933387B1 (en) * 2019-02-27 2023-10-18 Daikin Industries, Ltd. Information providing system
US11783917B2 (en) 2019-03-21 2023-10-10 Illumina, Inc. Artificial intelligence-based base calling
US11210554B2 (en) 2019-03-21 2021-12-28 Illumina, Inc. Artificial intelligence-based generation of sequencing metadata
US11423306B2 (en) 2019-05-16 2022-08-23 Illumina, Inc. Systems and devices for characterization and performance analysis of pixel-based sequencing
US11593649B2 (en) 2019-05-16 2023-02-28 Illumina, Inc. Base calling using convolutions
WO2021035098A2 (en) * 2019-08-21 2021-02-25 The Regents Of The University Of California Systems and methods for machine learning-based identification of sepsis
US11669949B2 (en) * 2019-10-22 2023-06-06 Canon U.S.A., Inc. Apparatus and method for inferring contrast score of an image
KR102268335B1 (ko) 2019-10-28 2021-06-24 주식회사 아이도트 자궁경부 영상 획득기
US11823378B2 (en) * 2019-12-20 2023-11-21 PAIGE.AI, Inc. Systems and methods for processing electronic images to detect contamination in specimen preparations
US11048971B1 (en) * 2019-12-24 2021-06-29 Ping An Technology (Shenzhen) Co., Ltd. Method for training image generation model and computer device
BR112022016415A2 (pt) 2020-02-20 2022-10-11 Illumina Inc Múltiplas chamadas de base à base de inteligência artificial
RU2734575C1 (ru) * 2020-04-17 2020-10-20 Общество с ограниченной ответственностью "АЙРИМ" (ООО "АЙРИМ") Способ и система идентификации новообразований на рентгеновских изображениях
US11288771B2 (en) * 2020-04-29 2022-03-29 Adobe Inc. Texture hallucination for large-scale image super-resolution
US20220000339A1 (en) * 2020-07-06 2022-01-06 Maine Medical Center Diagnostic cervical scanning and treatment device
US11688041B2 (en) 2021-03-02 2023-06-27 International Business Machines Corporation System and method of automatic image enhancement using system generated feedback mechanism
US20220336054A1 (en) 2021-04-15 2022-10-20 Illumina, Inc. Deep Convolutional Neural Networks to Predict Variant Pathogenicity using Three-Dimensional (3D) Protein Structures
CN113393461B (zh) * 2021-08-16 2021-12-07 北京大学第三医院(北京大学第三临床医学院) 基于深度学习的分裂中期染色体图像质量筛选方法和系统
US11798151B1 (en) * 2022-04-25 2023-10-24 Rivian Ip Holdings, Llc Systems and methods for determining image capture degradation of a camera sensor
CN116309379B (zh) * 2023-02-24 2023-11-03 飞燕航空遥感技术有限公司 基于多数据融合的航摄质量自动检查方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0236836A (ja) * 1988-04-22 1990-02-06 Olympus Optical Co Ltd 内視鏡画像処理装置
JP2002094979A (ja) * 2000-07-11 2002-03-29 Asahi Optical Co Ltd 電子内視鏡
JP2003047588A (ja) * 2001-08-03 2003-02-18 Olympus Optical Co Ltd 内視鏡装置
JP2003527916A (ja) * 2000-03-28 2003-09-24 ボード・オブ・リージェンツ,ザ・ユニヴァーシティ・オヴ・テキサス・システム 診断用マルチスペクトルデジタル画像化のための方法および装置
JP2004283421A (ja) * 2003-03-24 2004-10-14 Pentax Corp 電子内視鏡装置
JP2007506485A (ja) * 2003-09-26 2007-03-22 タイダール フォトニクス,インク. カラー画像内視鏡システムに関する装置と方法

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9408923D0 (en) * 1994-05-05 1994-06-22 Eaton Corp Power synchronizer for a compound transmission
US5565678A (en) 1995-06-06 1996-10-15 Lumisys, Inc. Radiographic image quality assessment utilizing a stepped calibration target
US6147705A (en) 1996-08-20 2000-11-14 Welch Allyn Inc. Apparatus and method for video colposcope with electronic green filter
US6081612A (en) * 1997-02-28 2000-06-27 Electro Optical Sciences Inc. Systems and methods for the multispectral imaging and characterization of skin tissue
US6738494B1 (en) * 2000-06-23 2004-05-18 Eastman Kodak Company Method for varying an image processing path based on image emphasis and appeal
US6713978B2 (en) * 2001-07-18 2004-03-30 Texas A&M University System Method and system for determining induction motor speed
US7133070B2 (en) * 2001-09-20 2006-11-07 Eastman Kodak Company System and method for deciding when to correct image-specific defects based on camera, scene, display and demographic data
US7038820B1 (en) * 2002-04-03 2006-05-02 Eastman Kodak Company Automatic exposure control for an image sensor
US6674885B2 (en) * 2002-06-04 2004-01-06 Amersham Biosciences Corp Systems and methods for analyzing target contrast features in images of biological samples
JP2006507775A (ja) * 2002-11-25 2006-03-02 サーノフ・コーポレーション 基準なしで圧縮ビデオシーケンスの品質を測定する方法及び装置
JP2007528137A (ja) * 2003-06-27 2007-10-04 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 非対称に強調される画像及びビデオのシャープネスメトリック
US7711174B2 (en) * 2004-05-13 2010-05-04 The Charles Stark Draper Laboratory, Inc. Methods and systems for imaging cells
US20060159325A1 (en) * 2005-01-18 2006-07-20 Trestle Corporation System and method for review in studies including toxicity and risk assessment studies
US7272207B1 (en) * 2006-03-24 2007-09-18 Richard Aufrichtig Processes and apparatus for variable binning of data in non-destructive imaging
US7940970B2 (en) * 2006-10-25 2011-05-10 Rcadia Medical Imaging, Ltd Method and system for automatic quality control used in computerized analysis of CT angiography

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0236836A (ja) * 1988-04-22 1990-02-06 Olympus Optical Co Ltd 内視鏡画像処理装置
JP2003527916A (ja) * 2000-03-28 2003-09-24 ボード・オブ・リージェンツ,ザ・ユニヴァーシティ・オヴ・テキサス・システム 診断用マルチスペクトルデジタル画像化のための方法および装置
JP2002094979A (ja) * 2000-07-11 2002-03-29 Asahi Optical Co Ltd 電子内視鏡
JP2003047588A (ja) * 2001-08-03 2003-02-18 Olympus Optical Co Ltd 内視鏡装置
JP2004283421A (ja) * 2003-03-24 2004-10-14 Pentax Corp 電子内視鏡装置
JP2007506485A (ja) * 2003-09-26 2007-03-22 タイダール フォトニクス,インク. カラー画像内視鏡システムに関する装置と方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JPN6013011713; Jia Gu and Wenjing Li: 'Automatic Image Quality Assessment for Uterine Cervical Imagery' Proceedings of SPIE - The International Society for Optical Engineering Vol. 6146, 2006, pp. 398-409, SPIE *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013066166A (ja) * 2011-08-30 2013-04-11 Ricoh Co Ltd 撮像装置、並びに、これを用いる画像解析装置及び移動装置
JP7517325B2 (ja) 2019-03-25 2024-07-17 ソニーグループ株式会社 医療システム、信号処理装置、及び、信号処理方法
WO2022014258A1 (ja) * 2020-07-17 2022-01-20 富士フイルム株式会社 プロセッサ装置、プロセッサ装置の作動方法
JP7547856B2 (ja) 2020-08-21 2024-09-10 カシオ計算機株式会社 撮像制御装置、撮像装置、撮像制御方法及びプログラム

Also Published As

Publication number Publication date
WO2008115405A3 (en) 2009-08-13
EP2137696A2 (en) 2009-12-30
WO2008115410A3 (en) 2009-08-13
US8401258B2 (en) 2013-03-19
WO2008115410A4 (en) 2009-10-15
US20080226147A1 (en) 2008-09-18
WO2008115410A2 (en) 2008-09-25
WO2008115405A2 (en) 2008-09-25
US20080226148A1 (en) 2008-09-18
US8295565B2 (en) 2012-10-23

Similar Documents

Publication Publication Date Title
US8401258B2 (en) Method to provide automated quality feedback to imaging devices to achieve standardized imaging data
US7454046B2 (en) Method and system for analyzing skin conditions using digital images
US20230190180A1 (en) Enhancing pigmentation in dermoscopy images
US8472682B2 (en) Medical image processing
US6902935B2 (en) Methods of monitoring effects of chemical agents on a sample
JP6917183B2 (ja) 撮像装置
US10467747B2 (en) Image analysis apparatus, imaging system, surgery support system, image analysis method, storage medium, and detection system
US20230363697A1 (en) Acne severity grading methods and apparatuses
WO2018161078A1 (en) Image adjustment and standardization
WO2022014258A1 (ja) プロセッサ装置、プロセッサ装置の作動方法
CN117314754B (zh) 一种双摄超光谱图像成像方法、系统及双摄超光谱内窥镜
CN114216867A (zh) 高光谱图像采集识别装置及方法
CN113784658A (zh) 用于生物组织的增强成像的系统和方法
JP7142754B1 (ja) 畳み込みニューラルネットワークを使用して食道癌の画像を検出する方法
TW202336707A (zh) 於超頻譜影像之物件偵測方法
CN116645565A (zh) 一种基于拉曼光谱和支持向量机分类器构建方法
TW202336690A (zh) 以波段用於超頻譜檢測物件影像之方法
Long et al. A Practical Dental Image Enhancement Network for Early Diagnosis of Oral Dental Disease
Vega García et al. Development of a Hyperspectral Colposcope for Early Detection and Assessment of Cervical Dysplasia
CN114627045A (zh) 医用图像处理系统、医用图像处理系统的工作方法
Voipio et al. Practical Imaging in Dermatology
CN114270399A (zh) 信息处理装置、脉搏波测量系统以及脉搏波测量程序

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20110118

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20121101

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130312

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20130903