JP6501800B2 - 信頼度マッチング付き生体内マルチカメラカプセルからの画像の再構築 - Google Patents

信頼度マッチング付き生体内マルチカメラカプセルからの画像の再構築 Download PDF

Info

Publication number
JP6501800B2
JP6501800B2 JP2016570804A JP2016570804A JP6501800B2 JP 6501800 B2 JP6501800 B2 JP 6501800B2 JP 2016570804 A JP2016570804 A JP 2016570804A JP 2016570804 A JP2016570804 A JP 2016570804A JP 6501800 B2 JP6501800 B2 JP 6501800B2
Authority
JP
Japan
Prior art keywords
image
images
matching
displayed
target image
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2016570804A
Other languages
English (en)
Other versions
JP2017522072A (ja
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 JP2017522072A publication Critical patent/JP2017522072A/ja
Application granted granted Critical
Publication of JP6501800B2 publication Critical patent/JP6501800B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B1/00Instruments for performing medical examinations of the interior of cavities or tubes of the body by visual or photographical inspection, e.g. endoscopes; Illuminating arrangements therefor
    • A61B1/04Instruments for performing medical examinations of the interior of cavities or tubes of the body by visual or photographical inspection, e.g. endoscopes; Illuminating arrangements therefor combined with photographic or television appliances
    • A61B1/041Capsule endoscopes for imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B1/00Instruments for performing medical examinations of the interior of cavities or tubes of the body by visual or photographical inspection, e.g. endoscopes; Illuminating arrangements therefor
    • A61B1/00002Operational features of endoscopes
    • A61B1/00004Operational features of endoscopes characterised by electronic signal processing
    • A61B1/00009Operational features of endoscopes characterised by electronic signal processing of image signals during a use of endoscope
    • A61B1/000095Operational features of endoscopes characterised by electronic signal processing of image signals during a use of endoscope for image enhancement
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B1/00Instruments for performing medical examinations of the interior of cavities or tubes of the body by visual or photographical inspection, e.g. endoscopes; Illuminating arrangements therefor
    • A61B1/00002Operational features of endoscopes
    • A61B1/00043Operational features of endoscopes provided with output arrangements
    • A61B1/00045Display arrangement
    • A61B1/0005Display arrangement combining images e.g. side-by-side, superimposed or tiled
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/22Matching criteria, e.g. proximity measures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/60Editing figures and text; Combining figures or text
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4038Image mosaicing, e.g. composing plane images from plane sub-images
    • 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
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/74Image or video pattern matching; Proximity measures in feature spaces
    • G06V10/75Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
    • 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/10068Endoscopic image
    • 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/20076Probabilistic 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/30004Biomedical image processing
    • G06T2207/30092Stomach; Gastric

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Theoretical Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • General Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Biomedical Technology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Signal Processing (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Multimedia (AREA)
  • Databases & Information Systems (AREA)
  • Computing Systems (AREA)
  • Quality & Reliability (AREA)
  • Endoscopes (AREA)
  • Image Processing (AREA)
  • Instruments For Viewing The Inside Of Hollow Bodies (AREA)

Description

本発明は2013年5月29日に出願され、発明の名称が「Reconstruction of Images from an in vivo Multi−Cameras Capsule」である米国仮特許出願第61/828,653号に関する。ここでは上記米国仮特許出願の全体内容は参照として引用されている。
本発明は、生体内カプセルカメラで撮像された画像からの画像ステッチおよびその表示に関する。
カプセル内視鏡は従来の内視鏡の様々な問題を対処する生体内イメージングデバイスである。カメラは、データをベースステーション受信器または送受信器に伝送する無線(radio)伝送器と共に、嚥下可能なカプセル内に収納される。体外のデータ記録器は伝送されたデータの受信および記録に使用されてもよい。そのデータは主にデジカルカメラに記録された画像を含む。カプセルは、ベースステーションから指令またはその他のデータを受信するための無線受信機を含んでもよい。無線周波数の送信を使用する代わりに、周波数の低い電磁信号を使用してもよい。電力は、誘導方式により外部インダクタからカプセルの内部インダクタに供給されてもよく、またはカプセル内のバッテリから供給されてもよい。内蔵(on−board)記憶装置付きの他のタイプのカプセルカメラにおいて、撮像された画像は外部デバイスへの伝送ではなく、内蔵記憶装置に保存される。内蔵記憶装置付きカプセルはカプセルの排出後に回収される。内蔵記憶装置付きカプセルは、データ記録器を取付けたり、無線データ受信器の近傍に制限されたりすることなく、患者には快適さと自由を与える。
進行方向に向ける(forward−looking)カプセルカメラはカメラを1つ含むが、サイドビューまたはパノラマビューを提供するマルチカメラを使用する他のタイプのカプセルカメラもある。生体の組織の表面を適切に観察するため、サイドアングルまたはリバースアングルが必要である。ポリープまたは他の異常は正確に診断するために徹底的に観察することが必要であるので、医師または診断医がこれらの器官の全てのエリアを見ることが重要である。2006年12月19日に出願され、発明の名称が「In vivo sensor with panoramic camera」である米国特許出願第11/642,275号に、カメラ周辺環境のパノラマ画像を撮影するカメラが開示されている。
自律(autonomous)カプセルシステムにおいて、カプセルカメラが胃腸(GI)管を通って移動する間に多数の画像が他のデータとともに収集される。通常、取得されてプロセスされた画像とデータは表示デバイスに表示されて診断医または医療専門家に検査される。しかし、各画像はGI管の小さいセクションの限られた視野(view)のみを提供する。多数のカプセル画像から単一の複合ビューを表する大きいピクチャを形成することが望ましい。例えば、多数のカプセル画像を内部GI管の表面の断面図の形成に使用することができる。大きいピクチャは、高い解像度の大きいスクリーン表示デバイスを利用することで、ユーザがさらに多くの情報を同時に見えるようにすることが可能である。画像ステッチ(stitching)プロセスは、画像の間の余分の重畳エリアを取り除くことを含めて、内部GI管の表面のより大きいエリアは単一複合ピクチャとして同時に見ることを可能にするようにしてもよい。さらに、大きいピクチャは内部GI管の表面の全景または重要な部分を提供することができる。そのため、診断医または医療専門家がより容易かつ迅速に関心のあるエリア(例えばポリープ)を特定することができなければならない。
コンピュテーショナル・フォトグラフィーの分野において、画像モザイク処理技術はより小さい画像から大きいピクチャにステッチするために開発された。画像アライメントおよびステッチに関する一般的な技術アプローチのレビューは題名が「Image Alignment and Stitching: A Tutorial」であり、2006年12月10日に公開され、著者がSzeliskiであるマイクロソフト・リサーチ・テクニカル・レポートMSR−TR−2004−92に記載されている。
画像モザイク処理では、まず、画像間の対応する部分、オブジェクトまたはエリアを識別する。対応する部分、オブジェクトまたはエリアを識別した後、すなわち、2つの画像が登録された後、対応する部分、オブジェクトまたはエリアを揃えることで、画像をステッチすることが可能である。2つの画像は直接に画素ドメインに登録されてもよく、または画像から抽出された特徴に基いてマッチングされてもよい。画素に基づく画像マッチングは直接マッチングと称してもよい。画像マッチングの品質を評価するための類似度測定はいくつかあり、例えば距離自乗和(sum of squared distance;SSD)、正規化相互相関(normalized cross correlation;NCC)、相互情報(mutual information;MI)などが挙げられる。SSDに対してノーマライゼーションを行えば、NCCはSSDと同様である。
特に、画像Aと画像Bの相互情報は以下のように定義し、2つ異なるモダリティからの画像をマッチングする。
Figure 0006501800
相互情報は、画像強度値の同時分布p(a,b)と、それらの画像が独立である時の画像強度値の同時分布p(a)p(b)との間の距離を測定する。これは、2つの画像の間の依存関係の測定である。前提は、それらの画像が正確にアラインする時に、画像の強度値の間に最大依存関係があることである。誤った登録は相互情報の測定値が減少する結果につながる。そのため、より大きい相互情報はより信頼性の高い登録を意味する。
特徴に基づくマッチングはまず各画像における特徴点の集合を決めてから対応の特徴記述子を比較する。2つの異なる視野角から取得された2つの画像パッチまたは特徴をマッチングするために、対応度に基いて、スケールやローテーションなどを含むリジッドモデルを評価すること。変形のオブジェクトからの2つの画像をマッチングするために、ローカルディフォーメーションを含む非リジッドモデルを計算してもよい。
米国特許出願第11/642,275号
通常、特徴点の数は対応画像の画素の数よりもはるかに少ない。そのため、特徴に基づく画像マッチングの計算負荷は実質的に画素に基づく画像マッチングの計算負荷よりも少ない。しかし、ペアワイズマッチングは依然時間がかかる。通常、この分野において公知技術であるk−d treeはこの手順のスピードアップに活用される。したがって、特徴に基づく画像マッチングはこの分野によく使用されている。それにもかかわらず、特徴に基づくマッチングはいくつかの状況において画像に対してうまく動作しない可能性がある。この場合、直接画像マッチングは予備(fall back)モードとして常に使用可能であり、または上記2つのアプローチの組合せが好ましい。
通常、画像マッチング技術は一定のモーションモデルを仮定する。カメラで撮像されたシーンがリジッドオブジェクトからなる時、特徴マッチングまたは画素ドメインマッチングに基づく画像マッチングは十分うまく動作する。しかし、もしシーンにおけるオブジェクトが変形したり、区別できる特徴を欠いたりする場合、画像マッチングタスクは非常に困難になる。GI管を通って移動間に撮像されたカプセル画像は状況がさらに厳しくなる。カメラの移動時におけるGI管の壁に対応するシーンが変形し、よく区別できる特徴を欠くだけでなく、シーンがカメラから短い距離で撮像されることもある。オブジェクトとカメラの間の距離が短いため、よく使用されていたリニアカメラモデルは異なるシーンの間のマッチングをよく行えない可能性がある。さらに、近いオブジェクトからの光の反射はオブジェクトのいくつかの部分の過度な露光を引き起こす可能性がある。そのため、上記の課題を克服できることが期待されている。
本開示は、画像マッチングの品質によりカメラで撮像された画像を適応的に表示する方法を開示している。画像ステッチは効率的なビューまたは莫大な数の画像の検査を提供するが、画像ステッチは、特にカメラモデルに適しない画像に対し、不自然な結果をもたらす。本発明は画像マッチングの品質を活用することで、基礎画像をステッチするか否かをガイドする。したがって、改善した画像の再構築およびより視覚的に見やすいビューを達成した。一実施形態において、カメラで撮像された複数の画像は受信され、画像の各ペアのための画像マッチングの品質が決定される。画像マッチングの品質はマッチングされる2つの画像の間で測定される。画像ペアは画像マッチングの品質により、高信頼度画像ペアまたは低信頼度画像ペアとされる。高信頼度画像ペアは1つまたはそれ以上のより大きい複合ピクチャにステッチされ、表示デバイスに表示される。一方、低信頼度画像ペアはステッチせず、個別の画像として表示される。
高信頼度マッチング画像と低信頼度アンマッチング画像は交互に表示デバイスに表示されてもよく、すなわち高信頼度マッチング画像と低信頼度アンマッチング画像は順番に表示される。高信頼度マッチング画像と低信頼度アンマッチング画像は同時に表示デバイスに表示されてもよい。この場合、高信頼度マッチング画像はある表示エリアに表示され、低信頼度アンマッチング画像は他の表示エリアに表示される。
本発明の一側面は画像マッチングの品質の様々な測定を対処する。例えば、画像マッチングの品質は、特徴が抽出された正確な画像マッチングに対応する事後確率に基づいてもよく、各マッチングされた特徴ペアはインライア(inlier)または外れ値(outlier)となるバイナリー確率変数としてモデル化される。この場合、画像マッチングの品質は単にインライアに属するマッチングされた特徴ペアの数をカウントすることで測定することが可能である。インライアに属する特徴の数が閾値を超えたら、画像ペアを高信頼度マッチング画像ペアとし、そうでなければ、画像ペアを低信頼度アンマッチング画像ペアとする。閾値はインライアに属する特徴の確率および外れ値に属する特徴の確率に依存する。
画像マッチングの品質は画素ドメイン測定に基いてもよい。例えば、画像マッチングの品質はマッチングされる2つの画像の間の距離自乗和(SSD)に基いてもよい。もしSSDは閾値よりも小さければ、画像ペアを高信頼度マッチング画像ペアとし、そうでなければ、画像ペアを低信頼度アンマッチング画像ペアとする。また、画像マッチングの品質はマッチングされる2つの画像の間の正規化相互相関(NCC)または相互情報(MI) に基いてもよい。もしNCCまたはMIは閾値を超えたら、画像ペアを高信頼度マッチング画像ペアとし、そうでなければ、画像ペアを低信頼度アンマッチング画像ペアとする。
グローバルトランスフォーメーションに加えてローカルディフォーメーションをも含むより複雑なケースについて、画像マッチングの品質は、画像ピラミッドと、フリーフォーム・ディフォーメーション、例えばB−splineディフォーメーションとの組合せで評価してもよい。よいグローバルマッチングのための検索は、オリジナル画像をダウンサンプリングすることで画像ピラミッドの粗いレベルから開始してもよい。各ピラミッドレベルにおいて、グローバルトランスフォーメーションは、ローカルシフトを平均し、または徹底的な検索を行うことで評価してもよい。このようなグローバルトランスフォーメーションは、オリジナル画像である最終レベルまで、画像ピラミッドの次のレベルにおいてリファインされる。この手順を繰り返すことで、外れ値効果を補償してもよい。グローバルトランスフォーメーションの評価後、重畳エリアを制御点/グリッドの集合に分割することで、重畳エリアにおいてフリーフォーム・ディフォーメーションを評価してもよい。最適化関数のアウトプットは、画像ペアを高信頼度マッチング画像ペアまたは低信頼度アンマッチング画像ペアに分類する信頼度スコアとして使用してもよい。
図1は本発明の一実施例による、画像マッチングに基づいた高信頼度と低信頼度画像ペアの典型的な決定、および同じ表示エリアにおいて交互に高信頼度マッチング画像とアンマッチングされた低信頼度画像を表示することを示す。 図2は本発明の一実施例による、画像マッチングに基づいた高信頼度と低信頼度画像ペアの典型的な決定、および高信頼度マッチング画像と低信頼度アンマッチング画像をそれぞれの表示エリアに表示することを示す。 図3は本発明の一実施例による、画像マッチングの品質にガイドされ、画像ステッチを含む、典型的な画像を表示するためのシステムのフローチャートを示す。
概して記載されており、図面に示されている本発明の構成要素は、幅広い種類の異なる形態にアレンジおよび設計されてもよいことは明らかにに分かる。したがって、下方に図面に示す本発明のシステムおよび方法の実施形態のより詳しい記載は、本発明の本発明のスコープを意図的に制限するものではなく、本発明の選択された実施形態の代表例に過ぎない。本明細書の全ての「1つの実施形態」、「一実施形態」または似たような文言は、実施例に関して記載されている特定の特徴、構造または特性は少なくとも本発明の1つの実施例に含まれてもよいことを意味する。したがって、本明細書における異なる箇所に記載の「1つの実施形態において」または「一実施形態において」というフレーズは全て同じ実施形態に属する必要がない。
さらに、記載の特徴、構造または特性を、1つまたはそれ以上の実施形態に適切な任意形式で組合せてもよい。しかし、関係領域における通常の知識を有する者は、1つまたはそれ以上の特定細部または他の方法や構成要素などがなくても、本発明は実行可能であると理解している。本発明の側面を分かりやすくするため、他の例において、公知構造または操作を詳細に図示しない、または記載しない。本発明の示された実施形態は図面を参照すれば、よく理解でき、図面において、類似部分は全て類似数字を振る。以下の記載はあくまでも例であり、請求項のように本発明と合致する装置と方法の特定の選択させた実施形態を簡単に示している。
前述のように、画像マッチングはいくつかの状況において画像(特に人間の胃腸(GI)管を通るカプセルで撮像した画像)に対してうまく動作しない可能性がある。本発明による実施形態は画像マッチングの品質測定を使用する。測定されたマッチング品質によれば、マッチング信頼度レベルは決定される。マッチング信頼度レベルが十分である時、基礎画像はステッチされる。そうでなければ、基礎画像はステッチされない。例えば、特徴に基づく画像マッチングが使用される場合、画像マッチングは実行され、様々な対応箇所をマッチングする。マッチング後、ランダム・サンプル・コンセンサス(RANdom Sample Consensus;RANSAC)プロセスは画像の間のトランスフォーメーション・モデルに適合するインライアの集合の選択に使用される。この分野において、RANSACは2つの画像の間の特徴点における最適なトランスフォームを探すために用いられる公知技術である。マッチングを検証するために、信頼度スコアは確率モデルに基いて計算される。
マッチング画像の各ペアについて、特徴対応箇所の部分集合は幾何学的に一致し(すなわち、RANSAC インライア)、残りの特徴は一致しない(すなわち、RANSAC外れ値)。対応箇所の集合を検証するためには、特徴の集合が正確な画像マッチング(すなわちインライア)または正しくない画像マッチング(すなわち外れ値)によって生成される確率は評価される。ある画像について、特徴の総数はnと記され、インライアの数はnと記され、画像マッチングの正確/不正確はバイナリー変数m∈{0,1}で記され、m=1は正確なマッチングを示し、m=0は不正確なマッチングを示す。i番目の特徴対応箇所がインライア/外れ値であり、バイナリー変数f(i)∈{0,1}で記され、式中、f=1はインライアを示し、f=0は外れ値を示すイベントは、独立ベルヌーイ分布であると仮定される。したがって、全ての特徴がインライアである確率は二項分布である:
Figure 0006501800
Figure 0006501800
式中、pは、特徴が正確な画像マッチングが得られるインライアである確率であり、pは、特徴が正しくない画像マッチングが得られるインライアである確率である。インライアの総数(n)は
Figure 0006501800
によって計算される。画像マッチングが正確である事後確率はベイズ定理で評価可能である:
Figure 0006501800
正確/不正確な画像マッチングのイベントをユニフォーム・プライア(すなわち、プライア確率分布)とすると、p(m=0)=p(m=1)である。画像マッチングを受け入れる基準は
Figure 0006501800
か否かに基いており、式中、pminは画像マッチングを受け入れる基準として使用される最小確率閾値である。数学的な展開のシーケンスを介し、この条件は尤度比検定(likelihood ratio test)となる:
Figure 0006501800
Figure 0006501800
min、pとpの値は画像モデルにより、またはテストデータに基いて、適切に選択される。上記決定プロセスはされに下記テストのように単純化してもよい:
Figure 0006501800
式中、gはpmin、pとpの関数である。すなわち、pmin、pとpの値の決定後、gは決定できる。決定プロセスは簡単にインライアの数(n)をカウントすることになる。式(7)の条件が満たされたら、画像マッチングは検証され、登録は信頼な登録とされる。そうでなければ、画像マッチング検証されず、登録は低信頼度を有する。上記実施形態において、画像マッチングの品質は、式(4)に示すように、特徴が抽出された画像マッチングが正確である事後確率として測定される。もし画像マッチングの品質は閾値(すなわち、Pmin)を超えたら、画像マッチングは検証される。さらなる展開後、本発明の1つの実施例による決定プロセスは簡単にインライアの数(n)をカウントし、その結果を閾値と比較することになる。インライアの数をカウントすることで画像マッチングの品質を測定することができるが、外れ値の数をカウントすることで画像マッチングの品質を測定してもよい。この場合、外れ値の数は第二閾値より少なかったら、画像マッチングは検証される。そうでなければ、画像マッチングは検証されない。
他の実施形態において、このシステムは非・特徴に基づく直接マッチングを利用し、画像マッチングの品質の測定として距離自乗和(SSD)を計算する。画像Aと画像Bの間のSSDは以下のように定義される:
Figure 0006501800
式中、(x,y)は重畳エリアにおける画素であり、Tは画像Aから画像Bへのトランスフォーメーションである。注意深く閾値Dmaxを選び、もしDssd(A,B|T)<Dmaxであれば、画像マッチングは検証され、登録は高信頼度を有する。そうでなければ、登録は検証されず、登録は信頼的でない。他の実施例において、正規化相互相関(NCC)または相互情報(MI)は、マッチングの品質を評価して信頼度スコアを計算する基準として使用してもよい。
他の実施形態において、連続する2つの画像をステッチするため、まず、そのペア画像の各画像はダウンサンプルされ、画像ピラミッドを作る。粗いレベルから、グローバルトランスフォーメーションは予め定めた範囲内に徹底的な検索で評価される。得られたグローバルトランスフォーメーションは、オリジナル画像である最終レベルまで、次のレベルにおいてリファインされる。グローバルトランスフォーメーションの評価後、フリーフォーム・ディフォーメーション・トランスフォーメーションは重畳エリアに適用され、ローカルディフォーメーションを評価してもよい。最適化オブジェクト関数のアウトプットは、マッチングの品質を評価して信頼度スコアを計算する基準として使用されてもよい。
他の実施形態において、連続する2つの画像をステッチするために、まず、そのペア画像の各画像はダウンサンプルされ、画像ピラミッドを作る。粗いレベルから、グローバルトランスフォーメーションは、フリーフォーム・ディフォーメーションを全体画像に適用させて計算されたローカルトランスフォーメーションを平均することで評価される。得られたグローバルトランスフォーメーションは、オリジナル画像である最終レベルまで、次のレベルにおいてリファインされる。プロセスが収束し、外れ値効果を除去するまでにこのような手順を繰り返す。グローバルトランスフォーメーションの評価後、フリーフォーム・ディフォーメーション・トランスフォーメーションは重畳エリアに適用され、ローカルディフォーメーションを評価してもよい。最適化オブジェクト関数のアウトプットは、マッチングの品質を評価して信頼度スコアを計算する基準として使用されてもよい。
他の実施形態において、以下の条件が真実である場合のみ、3つ以上の画像を高信頼度で一緒にステッチしてよい。画像il,i2, ... ,iNの集合の場合、各画像ij(j=l,2,...N)について、他の画像から高信頼度でijにマッチする画像を少なくとも1つ見つけることができる。高信頼度でijにマッチする画像は複数ある可能性がある。そうでなければ、ijは他の画像とステッチできず、この画像集合から取り除かれることを意味する。上記プロセスは、画像集合から取り除かれる画像がなくなるまで繰り返されてもよい。次に、この集合にある全ての画像は一緒にステッチし、大きい複合画像を形成してもよい。取り除かれた全ての画像は個別に表示される。
1つの実施例において、il, i2,...,iNは時間ドメインに沿う画像のシーケンスであり、il、i2、i3、i5、i6、i7、i8、il2はマッチングが高信頼度であるので、一緒にステッチされ、複合画像I1として表示されるが、i4、i9、i1Oとillはマッチングが高信頼度でなく、単一画像として表示される。もしi4とi9とillはマッチングが信頼であるが、i1Oはそうでなければ、映像において、i4、i9とillは複合画像I2として一緒にステッチされ、表示されるが、i1O単一画像として表示される。
あまり少ない画像をステッチし、1つの複合画像に表示することは、デメリットがメリットを上回る場合もある。例えば、単一画像の寸法とアスペクト比が固定される一方、ステッチされた画像は任意のサイズを有するので、一定のフレームレートで表示された映像では、2つの画像がステッチされた複合画像を見ることよりも、2つの画像を見たほうが時間的に効率的であることがある。画像がステッチされ、複合画像として表示される前に、閾値Tを選んで、高信頼度でマッチングされる画像の最小値を設定してもよい。
上記の画像マッチングの品質は画像ステッチをガイドするのに使用されることが可能である。画像マッチングの品質が高い時、登録は信頼であるとされてもよい。一実施形態において、画像の通過時間が不連続であっても、画像はステッチされ、画像のためのより大きい複合ピクチャを高信頼度で形成される。それら低信頼度だとされている画像について、画像はステッチされない。ステッチされない画像は個別の画像または画像シーケンスとして扱われ、映像として閲覧される。図1は本発明による1つの実施例を示し、Al、A2とA3は、映像全体にわたって、3つの高信頼度の画像グループを表す。Al、A2とA3はそれぞれ高信頼度を有する期間tA1、tA2とtA3の画像に対応する。各グループ(すなわち、Al、A2またはA3)の画像は1つまたはそれ以上より大きい複合ピクチャにステッチされる。Bl、B2とB3は低信頼度を有する期間tB1、tB2とtB3の画像に対応する。一実施形態において、Al、A2、A3と関連付けられた画像は表示エリア110に表示され、その次にBl、B2、B3と関連付けられた画像が表示されてもよい。図1は、グループA1と対応する複合ピクチャが表示される場合を示している。表示順番はAl、A2、A3、そしてBl、B2、B3であってもよい。表示はAl、B1、A2、B2、A3、B3との順番でもよい。Al、A2、A3と関連付けられた画像が表示される時、ステッチされたより大きい画像は同時に多数の画像の検査のための観察に使用されることができる。Bl、B2、B3と関連付けられた画像が表示される時、画像は個別の画像として扱われ、手動で1つずつ表示してもよく、または好ましい再生レートで映像シーケンスとして表示してもよい。図1において、画像は一定の速度で撮像される。他の実施形態において、画像は一定でないフレームレートで撮像されてもよい。
図2は本発明による他の実施例を示している。同様に、Al、A2とA3は映像全体にわたって高信頼度の画像を表す。Bl、B2とB3は低信頼度を有する画像に対応する。表示エリアが2つあり、1つはA1、A2とA3を表示し、他の表示エリアはB1、B2とB3を表示する。2つの表示エリア(210と220)は独立してA1/A2/A3とB1/B2/B3を表示する。Al、A2、A3と関連付けられた画像は表示エリア210にステッチされたより大きい複合画像として表示されることができる。Bl、B2、B3と関連付けられた画像は個別の画像として表示されることができる。これらの画像は、表示エリア220において、手動で1つずつ表示されてもよく、または好ましい再生レートで映像として表示されてもよい。
図3は、本発明の実施例によるシステムの典型的なフローチャートを示しており、このシステムは画像マッチングの品質でガイドされた画像ステッチを含む画像を表示する。カメラで撮像された複数の画像はステップ310に示すように受信される。画像はメモリから回収されてもよく、プロセッサから受信されてもよい。各画像の画像マッチングの品質はステップ320に示すように決定される。画像マッチングの品質はマッチングされる2つの画像の間で測定される。そして、画像ペアはステップ330に示すように画像マッチングの品質により、高信頼度マッチング画像ペアまたは低信頼度アンマッチング画像ペアとされる。高信頼度画像はステップ340に示すように1つまたはそれ以上のより大きい複合ピクチャにステッチされる。ステッチされたより大きい複合ピクチャはステップ350に示すように表示デバイスに表示される。
特定の例はカプセル画像に向けているが、本発明による画像マッチングの品質に基づく画像ステッチは異なる視野角で撮像された自然なシーンの画像に適用してもよい。
本発明は、本発明の精神または本質的な特性から逸脱しない限り、他の特定の形態に具体化してもよい。上記の例は全ての側面において単に説明だと考えられべき、制限ではない。そのため、本発明のスコープは上記の説明ではなく、特許請求の範囲に示す通りである。特許請求の範囲と均等する意味および範囲に該当する全ての変更は本発明のスコープに含まれる。
110 表示エリア
210 表示エリア
220 表示エリア

Claims (8)

  1. カプセルカメラで撮像された複数の画像を受信するステップと、
    前記複数の画像から選択された対象画像集合における対象画像毎に、
    選択された対象画像に対応する第一の画像と、前記対象画像集合におけるもう一つの画像に対応する第二の画像とにより構成された各画像ペアと関連付けられた画像マッチングの品質を決定するステップと、
    前記画像マッチングの品質の決定において、各画像ペアについて、インライア特徴の数をカウントする、または前記第一の画像と前記第二の画像の間の一つ又はそれ以上の重畳エリアについて距離自乗和(SSD)を決定するステップと、
    前記インライア特徴の數が第一の閾値より大きいこと、または前記距離自乗和(SSD)が第二の閾値より小さいことに対応する前記画像マッチングの品質により、前記選択された対象画像と関連付けられた少なくとも一つの対応の画像ペアが高信頼度マッチングを有する場合、前記選択された対象画像をマッチング画像とするステップと、
    前記選択された対象画像と関連付けられた全ての対応の画像ペアは前記第一の閾値より小さいまたは等しいインライア特徴の數、もしくは前記第二の閾値より大きいまたは等しい距離自乗和(SSD)有する場合、前記選択された対象画像をアンマッチング画像とするステップと、
    を実行するステップと、
    前記対象画像集合におけるマッチングした画像のみを1つまたはそれ以上の複合ピクチャにステッチするステップと、
    前記1つまたはそれ以上の複合ピクチャを表示デバイスに表示するステップと、
    を含む、人間の胃腸(GI)管を通って移動するカプセルカメラで撮像されたGI管の画像を表示する方法。
  2. 前記1つまたはそれ以上の複合ピクチャと前記対象画像集合におけるアンマッチング画像を前記表示デバイスに交互に表示し、前記1つまたはそれ以上の複合ピクチャを第一期間に表示し、前記アンマッチング画像を第二期間に表示し、前記第一期間と前記第二期間は重畳しない、
    請求項1に記載の方法。
  3. 前記1つまたはそれ以上の複合ピクチャを前記表示デバイスの第一表示エリアに表示し、前記対象画像集合におけるアンマッチング画像を前記表示デバイスの第二表示エリアに表示する、
    請求項1に記載の方法。
  4. 前記対象画像集合におけるマッチングした画像のみを1つまたはそれ以上の複合ピクチャにステッチするステップは、1つの複合ピクチャについて対応の画像ペアの数が閾値以上である場合のみに実行される、
    請求項1に記載の方法。
  5. 表示デバイスと、前記表示デバイスに接続されたプロセッサと、を含む、人間の胃腸(GI)管を通って移動するカプセルカメラで撮像されたGI管の画像を表示するシステムであって、
    前記プロセッサは、
    前記カプセルカメラで撮像された複数の画像を受信し、
    選択された対象画像に対応する第一の画像と、前記対象画像集合におけるもう一つの画像に対応する第二の画像とにより構成された各画像ペアと関連付けられた画像マッチングの品質を決定し、
    前記画像マッチングの品質の決定において、各画像ペアについて、インライア特徴の数をカウントする、または前記第一の画像と前記第二の画像の間の一つ又はそれ以上の重畳エリアについて距離自乗和(SSD)を決定し、
    前記インライア特徴の數が第一の閾値より大きいこと、または前記距離自乗和(SSD)が第二の閾値より小さいことに対応する前記画像マッチングの品質により、前記選択された対象画像と関連付けられた少なくとも一つの対応の画像ペアが高信頼度マッチングを有する場合、前記選択された対象画像をマッチング画像とし、
    前記選択された対象画像と関連付けられた全ての対応の画像ペアは前記第一の閾値より小さいまたは等しいインライア特徴の數、もしくは前記第二の閾値より大きいまたは等しい距離自乗和(SSD)有する場合、前記選択された対象画像をアンマッチング画像とし、
    前記対象画像集合におけるマッチングした画像のみを1つまたはそれ以上の複合ピクチャにステッチし、
    前記1つまたはそれ以上の複合ピクチャを表示デバイスに表示するようにされている、
    システム。
  6. 前記対象画像集合におけるアンマッチング画像がステッチせず、個別の画像として前記表示デバイスに表示される、
    請求項に記載のシステム。
  7. 前記1つまたはそれ以上の複合ピクチャと前記アンマッチング画像が前記表示デバイスに交互に表示され、前記1つまたはそれ以上の複合ピクチャが第一期間に表示され、前記アンマッチング画像が第二期間に表示され、前記第一期間と前記第二期間は重畳しない、
    請求項に記載のシステム。
  8. 前記1つまたはそれ以上の複合ピクチャが前記表示デバイスの第一表示エリアに表示され、前記アンマッチング画像が前記表示デバイスの第二表示エリアに表示される、
    請求項に記載のシステム。
JP2016570804A 2014-06-01 2015-04-27 信頼度マッチング付き生体内マルチカメラカプセルからの画像の再構築 Active JP6501800B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201462006257P 2014-06-01 2014-06-01
US62/006,257 2014-06-01
PCT/US2015/027813 WO2015187261A1 (en) 2014-06-01 2015-04-27 Reconstruction of images from an in vivo multi-camera capsule with confidence matching

Publications (2)

Publication Number Publication Date
JP2017522072A JP2017522072A (ja) 2017-08-10
JP6501800B2 true JP6501800B2 (ja) 2019-04-17

Family

ID=54767145

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016570804A Active JP6501800B2 (ja) 2014-06-01 2015-04-27 信頼度マッチング付き生体内マルチカメラカプセルからの画像の再構築

Country Status (5)

Country Link
US (1) US20170055814A1 (ja)
EP (1) EP3148399B1 (ja)
JP (1) JP6501800B2 (ja)
CN (1) CN106455956A (ja)
WO (1) WO2015187261A1 (ja)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018101936A1 (en) * 2016-11-30 2018-06-07 CapsoVision, Inc. Method and apparatus for image stitching of images captured using a capsule camera
WO2018194599A1 (en) * 2017-04-19 2018-10-25 CapsoVision, Inc. Method of image processing and display for images captured by a capsule camera
CN108022211A (zh) * 2017-10-27 2018-05-11 平安科技(深圳)有限公司 电子装置、车祸现场全景图像展示方法和存储介质
CN108024061A (zh) * 2017-12-08 2018-05-11 合肥工业大学 医用内窥镜人工智能系统的硬件架构及图像处理方法
US10510145B2 (en) * 2017-12-27 2019-12-17 Industrial Technology Research Institute Medical image comparison method and system thereof
CN108335319A (zh) * 2018-02-06 2018-07-27 中南林业科技大学 一种基于自适应阈值及ransac的图像角点匹配方法
CN113409198B (zh) * 2021-08-19 2021-11-05 北京欧应信息技术有限公司 用于图像拼接的方法、计算设备和存储介质

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6075905A (en) * 1996-07-17 2000-06-13 Sarnoff Corporation Method and apparatus for mosaic image construction
US6044181A (en) * 1997-08-01 2000-03-28 Microsoft Corporation Focal length estimation method and apparatus for construction of panoramic mosaic images
JP4631057B2 (ja) * 2004-02-18 2011-02-16 国立大学法人大阪大学 内視鏡システム
CN1284505C (zh) * 2004-02-28 2006-11-15 重庆金山科技(集团)有限公司 医用无线电胶囊式内窥系统
CN101065052B (zh) * 2004-12-27 2010-12-22 奥林巴斯株式会社 医疗图像处理装置和医疗图像处理方法
JP2007236700A (ja) * 2006-03-09 2007-09-20 Fujifilm Corp カプセル内視鏡システム
US20100149183A1 (en) * 2006-12-15 2010-06-17 Loewke Kevin E Image mosaicing systems and methods
US7978932B2 (en) * 2007-08-02 2011-07-12 Mauna Kea Technologies Robust mosaicing method, notably with correction of motion distortions and tissue deformations for in vivo fibered microscopy
US20090074265A1 (en) * 2007-09-17 2009-03-19 Capsovision Inc. Imaging review and navigation workstation system
WO2011003168A1 (en) * 2009-07-10 2011-01-13 Isee3D Inc. Method and apparatus for generating three dimensional image information using a single imaging path
DE102009039251A1 (de) * 2009-08-28 2011-03-17 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Verfahren und Vorrichtung zum Zusammenfügen von mehreren digitalen Einzelbildern zu einem Gesamtbild
US8724868B2 (en) * 2009-10-12 2014-05-13 Capso Vision, Inc. System and method for display of panoramic capsule images
CN103347436B (zh) * 2011-01-11 2016-06-22 爱德芳世株式会社 口腔内拍摄显示系统
US8712151B2 (en) * 2011-02-14 2014-04-29 Intuitive Surgical Operations, Inc. Method and structure for image local contrast enhancement
CN102436665A (zh) * 2011-08-25 2012-05-02 清华大学 消化道图像二维平面表示方法
CN103190881A (zh) * 2012-01-04 2013-07-10 清华大学 胶囊式内窥镜及其图像处理方法

Also Published As

Publication number Publication date
EP3148399B1 (en) 2021-09-22
EP3148399A1 (en) 2017-04-05
US20170055814A1 (en) 2017-03-02
CN106455956A (zh) 2017-02-22
WO2015187261A1 (en) 2015-12-10
EP3148399A4 (en) 2018-01-17
JP2017522072A (ja) 2017-08-10

Similar Documents

Publication Publication Date Title
JP6501800B2 (ja) 信頼度マッチング付き生体内マルチカメラカプセルからの画像の再構築
CN107920722B (zh) 针对自胶囊相机所撷取的图像通过对象检测进行重建
Ozyoruk et al. EndoSLAM dataset and an unsupervised monocular visual odometry and depth estimation approach for endoscopic videos
CN111557665B (zh) 用于医学成像中的患者定位的根据相机数据的密集身体标记估计
Bergen et al. Stitching and surface reconstruction from endoscopic image sequences: a review of applications and methods
JP4631057B2 (ja) 内視鏡システム
US10198872B2 (en) 3D reconstruction and registration of endoscopic data
KR102237441B1 (ko) 뉴럴 네트워크를 이용하여 캡슐 내시경 영상으로부터 병변 판독 방법 및 장치
US20170084036A1 (en) Registration of video camera with medical imaging
CN104000655B (zh) 用于腹腔镜外科手术的组合的表面重构和配准
CN111035351B (zh) 用于胃肠道中的胶囊相机的行进距离测量的方法及装置
CN102065744A (zh) 图像处理装置、图像处理程序及图像处理方法
US20160295126A1 (en) Image Stitching with Local Deformation for in vivo Capsule Images
Seshamani et al. Real-time endoscopic mosaicking
CN111127461B (zh) 胸部图像处理方法、装置、存储介质及医疗设备
JP6824078B2 (ja) 内視鏡位置特定装置、方法およびプログラム
US11120547B2 (en) Reconstruction of images from an in vivo multi-camera capsule with two-stage confidence matching
Ahmedt-Aristizabal et al. Monitoring of pigmented skin lesions using 3D whole body imaging
Safavian et al. Endoscopic measurement of the size of gastrointestinal polyps using an electromagnetic tracking system and computer vision-based algorithm
WO2016160862A1 (en) Method of overlap-dependent image stitching for images captured using a capsule camera
Ozyoruk et al. Quantitative evaluation of endoscopic slam methods: Endoslam dataset
Hu et al. Homographic patch feature transform: a robustness registration for gastroscopic surgery
Hackner et al. Deep-learning based reconstruction of the stomach from monoscopic video data
Hackner et al. 3D Reconstruction of the Colon from Monocular Sequences Evaluation by 3D-printed Phantom Data
Hu et al. Image Registration

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20171027

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20171113

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20171227

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180205

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180507

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20181029

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20181207

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20181217

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190319

R150 Certificate of patent or registration of utility model

Ref document number: 6501800

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250