JP2017506384A - 少なくとも二画像をレジストレーション及び視覚化するための方法 - Google Patents

少なくとも二画像をレジストレーション及び視覚化するための方法 Download PDF

Info

Publication number
JP2017506384A
JP2017506384A JP2016549165A JP2016549165A JP2017506384A JP 2017506384 A JP2017506384 A JP 2017506384A JP 2016549165 A JP2016549165 A JP 2016549165A JP 2016549165 A JP2016549165 A JP 2016549165A JP 2017506384 A JP2017506384 A JP 2017506384A
Authority
JP
Japan
Prior art keywords
image
navigation speed
deformation
images
solving
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
JP2016549165A
Other languages
English (en)
Other versions
JP6573894B2 (ja
JP2017506384A5 (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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
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 Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of JP2017506384A publication Critical patent/JP2017506384A/ja
Publication of JP2017506384A5 publication Critical patent/JP2017506384A5/ja
Application granted granted Critical
Publication of JP6573894B2 publication Critical patent/JP6573894B2/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
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/14Transformations for image registration, e.g. adjusting or mapping for alignment of images
    • G06T3/147Transformations for image registration, e.g. adjusting or mapping for alignment of images using affine transformations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/14Transformations for image registration, e.g. adjusting or mapping for alignment of images
    • G06T3/153Transformations for image registration, e.g. adjusting or mapping for alignment of images using elastic snapping
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F3/00Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
    • G06F3/01Input arrangements or combined input and output arrangements for interaction between user and computer
    • G06F3/048Interaction techniques based on graphical user interfaces [GUI]
    • G06F3/0484Interaction techniques based on graphical user interfaces [GUI] for the control of specific functions or operations, e.g. selecting or manipulating an object, an image or a displayed text element, setting a parameter value or selecting a range
    • G06F3/0485Scrolling or panning
    • 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/10004Still image; Photographic image
    • G06T2207/10008Still image; Photographic image from scanner, fax or copier
    • 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/10056Microscopic image
    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

本発明は少なくとも二画像をレジストレーション及び視覚化するための方法に関し、方法は、少なくとも二画像の第1画像のスクロール可能及び/又はズーム可能な視覚化中にユーザによって使用されるナビゲーション速度を測定し、ナビゲーション速度をナビゲーション速度閾値と比較するステップS1と;ナビゲーション速度がナビゲーション速度閾値を下回る場合、第1レベルの計算負荷を用いて第1画像の変形を解くことによって第1画像を視覚化するステップS2と;ナビゲーション速度がナビゲーション速度閾値を上回る場合、第2レベルの計算負荷を用いて第1画像の変形を解くことによって第1画像を視覚化するステップS3とを有する。

Description

本発明は少なくとも二画像をレジストレーション及び視覚化するための方法、並びにデジタル画像レジストレーション及び視覚化のためのイメージングシステムに関する。
US2013/0089249 A1はデジタル病理学における画像処理法を記載する。記載の処理法は領域内のサンプルをあらわす第1デジタル画像を処理し、当該画像は顕微鏡イメージングシステムを用いてサンプルから収集されている。画像は多重解像度画像データ構造に保存され、第1解像度で第1デジタル画像のサブ領域を読み出し、読み出されたサブ領域に対して変換関数を実行し、変換関数は第1デジタル画像の第2解像度表現から導出される少なくとも一つのメトリックに従ってサブ領域のコンテンツを修正する。
US2008/013852 A1は画質を改良するためのシステム、方法及びコンピュータプログラム製品を記載する。記載の方法によれば、二つの異なるスケール空間値に基づいてピクセルの二つのリストを生成するために入力画像のピクセルが使用される。その後、二つのピクセルリストが処理され、改良された画質を持つ出力画像を取得する。US2010/157079 A1は画像を選択的に組み合わせるシステムと方法を記載する。装置は第1画像に対応する第1画像データと第2画像に対応する第2画像データに基づいて動きベクトルデータのセットを生成するように構成されるレジストレーション回路を含む。
US2013016886 A1は組織サンプルの隣接スライスの高解像度デジタルイメージングのための方法を記載する。各タイルにおいて検出される画像オブジェクトに関連する値が計算される。隣接画像中のタイルがコレジストレーションされる。第1ハイパースペクトル画像が第1画像を用いて生成され、第2ハイパースペクトル画像が第2画像を用いて生成される。
US2012/068928 A1はオペレーションのセットを含むデジタル病理学環境において標本スライドを検査及び操作するための装置と方法を記載する。いずれかの方向にスクロールホイールを素早く動かすことが、最小若しくは最大倍率へのジャンプ、又は次の順次インクリメントを超える前進への指示として構成され得、中間倍率の負荷を未然に防ぐ。
高速若しくはリアルタイムの複雑な画像レジストレーション及び視覚化のためにデジタル画像処理を改良する必要性があり得る。
これらの必要性は独立クレームの主題によって満たされる。さらなる実施形態例は従属クレーム及び以下の記載から明らかとなる。
本発明の一態様は少なくとも二画像をレジストレーション及び視覚化するための方法に関し、当該方法は、少なくとも二画像の第1画像のスクロール可能及び/又はズーム可能な視覚化中にユーザによって使用されるナビゲーション速度を測定し、ナビゲーション速度をナビゲーション速度閾値と比較するステップ;ナビゲーション速度がナビゲーション速度閾値を下回る場合、第1レベルの計算負荷を用いて第1画像の変形を解くことによって第1画像を視覚化するステップ;並びにナビゲーション速度がナビゲーション速度閾値を上回る場合、第2レベルの計算負荷を用いて第1画像の変形を解くことによって第1画像を視覚化するステップを有する。
本発明のさらなる態様はデジタル画像レジストレーション及び視覚化のためのイメージングシステムに関し、当該イメージングシステムは、ナビゲーション速度をナビゲーション速度閾値と比較するように構成される画像プロセッサであって、当該ナビゲーション速度は少なくとも二画像の第1画像のスクロール可能及び/又はズーム可能な視覚化中にユーザによって使用され、当該画像プロセッサはさらに、ナビゲーション速度がナビゲーション速度閾値を下回る場合、第1レベルの計算負荷を用いて第1画像の変形を解くことによって第1画像を視覚化するように構成され、当該画像プロセッサはさらに、ナビゲーション速度がナビゲーション速度閾値を上回る場合、第2レベルの計算負荷を用いて第1画像の変形を解くことによって第1画像を視覚化するように構成される、画像プロセッサと、少なくとも二画像を表示するように構成されるディスプレイとを有する。
本発明はリアルタイム若しくは高速の複雑な画像レジストレーション及び視覚化を実行することを可能にする。二画像がレジストレーションのために使用される。両画像間の粗い変換を特定するアルゴリズムが実行される。両画像間の精密な変換を特定するアルゴリズムが実行される。連続画像間の変化量を示すアルゴリズムが実行され、これは粗雑から精密へ及びその逆の変換をトリガするために使用される。
本発明は、ユーザへのフィードバックタイムを削減するためにレジストレーションが可能な限り速く起こるように、デジタル病理画像などサイズの大きい二画像をレジストレーションするための方法を有利に提供する。方法は粗雑レジストレーション及び視覚化の後精密レジストレーション及び視覚化が続くことに基づく。
サンプル関心領域(ROI)及び/又は残りのサンプルの画像が、サンプルROIの分離後に生成される。この画像はサンプル全体の画像も生成したサンプル画像生成ユニットで、又は別の装置で生成され得る。サンプルROIの(若しくはサンプルの残りの)画像はサンプル全体の画像に、特に選択された画像ROIに比較されることができ、従って実際のサンプルROIが所望の関心領域に対応するか否かの検証を可能にする。
本発明の一実施形態例によれば、第1画像の変形を解くステップは第1画像の画像並進、第1画像の画像回転、第1画像のスケール変動、第1画像の追加線形変形若しくは第1画像の非線形変形を有する。
本発明の一実施形態例によれば、第1画像のスクロール可能及び/又はズーム可能な視覚化中にユーザによって使用されるナビゲーション速度はユーザインターフェース装置から取得される。
二つのデジタル画像のレジストレーションの場合、ナビゲーション速度はユーザの最終更新要求からの時間によって有利に決定され得る。ユーザが絶えずパン及び/又はズームしている場合、精密変換は省略される。ユーザが視点を変えない場合、例えばパン及び/又はズームしない場合、精密変換が実行される。
本発明の一実施形態例によれば、第1画像のスクロール可能及び/又はズーム可能な視覚化中にユーザによって使用されるナビゲーション速度はSIFT特徴マッチングを用いて測定される。
本発明の一実施形態例によれば、方法は少なくとも二画像の局所レジストレーションを実行するステップをさらに有する。
本発明の一実施形態例によれば、ナビゲーション速度がナビゲーション速度閾値を上回る場合、第1画像の変形を解くことは並進及び回転変換を適用することを有し、ナビゲーション速度がナビゲーション速度閾値を下回る場合、第1画像の変形を解くことは非線形精緻化を適用することを有する。
本発明の一実施形態例によれば、ナビゲーション速度がナビゲーション速度閾値を上回る場合、第1画像の変形を解くことは第1解像度レベルで変換を適用することを有し、ナビゲーション速度がナビゲーション速度閾値を下回る場合、第1画像の変形を解くことは第2解像度レベルで変換を適用することを有し、第1解像度レベルは第2解像度レベルよりも低い。
本発明はさらに分子診断、分子病理学のために、特にオンコロジーアプリケーション、生物試料分析、化学試料分析、食品分析、及び/又は法医学的分析のために使用されるデジタル画像解析に関する。分子診断は例えば標的分子に直接若しくは間接的に付着される磁気ビーズ若しくは蛍光粒子を利用して達成され得る。
本発明及びそれに付随する利点のより完全な評価は、原寸通りではない以下の略図を参照してより明確に理解される。
本発明にかかるサンプルの検査を概略的に示す。 本発明を説明するための二つの免疫組織化学染色スライドの画像レジストレーションの例を示す。 本発明の一実施形態例にかかる構成要素の概観を示す。 本発明の一実施形態例にかかる少なくとも二画像をレジストレーション及び視覚化するための方法の略フローチャート図を示す。
図中の例示は概略であって原寸通りではない。異なる図において、同様の若しくは同一の要素は同じ参照数字を与えられる。本特許出願で使用される"ナビゲート"という語は、全体が見えていないより大きな画像のように一見思われるものにわたる、ユーザのビューのスクロール、パン、ズーム若しくは傾斜を含む、任意の動きをあらわす。言い換えれば本特許出願にかかるスクロールとは、コンピュータのディスプレイ上でデータのビューを変更することを意味する。スクロールは離散的単位で、おそらく画像の一つ若しくはいくつかの部分を一度に若しくは連続的に、起こり得る。
概して、同一の部分、ユニット、エンティティ若しくはステップは図において同じ参照符号を与えられる。
患者の材料(組織及び細胞)の病理学診断検査は、特にオンコロジーにおいて、多くの治療診断の基礎である。標準的な、生検由来の薄いスライスが顕微鏡スライド上に示され、組織の形態学を視覚化するために特定プロトコルに従って染色される。より最近は疾患特異的バイオマーカについてのin‐situ染色が標的薬の比較診断のために開発されている。評価は明視野顕微鏡でなされ得る。診断が再評価される必要がある場合のバックアップとして、スライドは検査後に長期間保存される必要がある。
デジタル病理学は、顕微鏡がデジタルスキャナによって置き換えられる新しい技術であり、デジタルスキャナは被染色組織切片を自動的にスキャンし、十分な解像度とカラーレンダリングを伴うデジタル形式で画像を保存するので、病理学者は顕微鏡において直接し得るのと同じ診断をデジタル画像からすることができる。後者はデジタル画像が物理的スライドを置き換えることができることを意味する。デジタル画像がスライドの代わりに保存され得る。
上記形式の分析は別にして、組織及び細胞生検はq‐PCR及びシークエンシングのような分子法("分子診断"若しくは"MDx"と略される)でも検査され得る。このいわゆる分子病理学は新たな分子バイオマーカの出現により重要性が高まっている。病理学者は正しい治療選択のために癌組織の生物学特徴を同定する分子検査を行うことを形態学的情報に基づいて決定することが多い。多くの分子バイオマーカは組織上でin situに定量化されることができず、又は少なくとも要求精度では不可能であり、PCR若しくはシークエンシングのような別の分子診断テストが、一般に生検から既に採取されているクーペからの、生検から採取されるサンプル上で実行される。この組織切片はDNA若しくはmRNAマーカの測定前に細胞溶解によって処理される。結果として空間情報が失われる。
腫瘍組織は一般に癌細胞だけでなく、異なる細胞型から成り、癌細胞でも腫瘍の異なるエリアにおいて分子構造が大いに異なり得る。従って分子解析の結果は分子検査用のサンプルとして使用される組織切片の厳密な組成に依存する。癌細胞がより希釈されるほど、検査結果は低感度で不確定になる。加えて、癌細胞集団内の不均質性はMDxアッセイにおけるノイズも生じ、感度と特異性だけでなく再現性も削減する。再現性は単に各サンプルが固有であるという理由で一般に主要な問題である。
デジタル病理学では、現在生検クーペのサブ切片上で分子検査を行う可能性はない。全クーペ上での検査は、異なる出所(例えば内皮、繊維芽細胞及び免疫細胞)の良性細胞と癌細胞不均質性による標的細胞の希釈のために、準最適な精度と感度を持つ。
マニュアル選択は正確さを欠くことが多く、従って非常に低濃度の汚染さえも増幅するPCRの場合特に、汚染を生じやすい。
関心細胞の純粋画分を取得することは時間がかかる可能性があり、除去により高解像度を要する。上述のアプローチは選択の要件を緩和することを可能にし、どれ位容易な若しくは信頼できる切片が除去されることができるかを考慮する追加パラメータとのバランスをとり、選択の総サイズを制御する。好都合な選択のために、より大きな連続切片が好適である。
画像解析は、最終的に解剖された組織における同定された細胞型の各々の数及び割合のような、表形式パラメータを提供することができる。エリアの形状は、総面積、許容曲率及び結合性のようなパラメータを含む設計基準によって制限され得る。アルゴリズムに基づいて、最適条件が各MDxに特異的であり得る所定選択アルゴリズムを所与として決定され得る。
このアプローチの一実施形態例において、組織スライドはHER2免疫組織化学若しくは免疫蛍光染色、又は染色アッセイの組み合わせで所定の臨床的適応に従って染色され得る。
そしてスライドはデジタルスキャナによってスキャンされ得、得られる画像は共通の特徴のエリアを同定し示すためにコンピュータプログラムによって解析され得る。こうしたエリアは必要な場合若しくはその他望ましい場合に、確認及び適応のために病理学者に提示され得る。これらのエリアから"サンプルROI"とよばれる関心領域がソフトウェアプログラムによって定義され得、これはMDx検査のために選択されるサンプルの一部をあらわす。
所与の実施例におけるHER2の平均発現、及び細胞全体の発現の統計的分布のような、典型的なパラメータがそのサンプルROIに注釈付けされる。選択されたサンプルROIの座標はMDx用のサンプルの物理的選択を処理するサンプル分離ユニットへ送信され得る。選択された"サンプルROI"は転送装置へ転送され得るか、又は直接、MDx検査のために使用される消耗品であり得る。
MDx検査は、qRT‐PCR、qrt‐PCR、シークエンシング若しくは次世代シークエンシング、又はマイクロアレイハイブリダイゼーション、又はこれらの組み合わせなど、サンプル調製と分子解析を有し得る。その分析の結果は最終的に組織選択アルゴリズムからの情報に関連付けられ、オプションとして解釈され、MDx用に選択されたサンプルROIも示される組織のデジタル画像と一緒に病理学者へ提示される。
図1は本発明にかかるサンプルの検査を概略的に示す。
検査はサンプル調製ユニット100で開始し、ここで生体組織のスライス11がキャリアとして機能する顕微鏡スライド10上に準備される。典型的に、組織サンプルはパラフィン包埋生検から約4‐8ミクロンの薄い切片を切り出すことによって取得される。いわゆるクーペが水膜上の顕微鏡ガラススライド10上に置かれ、ミクロトーム負荷から緩和し、そして乾燥させる。さらに、サンプル11はオプションとして適切な染色12によって、例えばヘマトキシリン‐エオジン(H&E)若しくはIHCによって染色され得る。異なるアプリケーションに利用可能な複数の染色プロトコルがある。染色プロトコルは試薬を含む異なる溶液中にクーペとともにスライドを浸すことによって手動でベンチ上で実行され得るが、自動化された方法で実行されることもできる。
一つ以上のマーカMが顕微鏡スライド10上にプリントされるか若しくはその中に彫刻され、これは後にスライド10上の実際の座標に画像座標を関連付けるための基準点として機能し得る。さらに、後のスキャン中に高画質を可能にするため、及び長期保管期間にわたってサンプルの保存を保証するために、スライド10はサンプル11とともに好適にはカバースリップ(不図示)でカバーされる。
調製ステップ後、スライド10はサンプル11とともに特にデジタルスキャン顕微鏡であり得る画像生成ユニット200へ移される。
サンプルのデジタル画像Iは顕微鏡200によって生成される。顕微鏡はイメージングシステム300に接続され、これは顕微鏡200をビデオ画像源として使用し、ディスプレイ301とビデオ画像源データのデジタル画像処理用の画像プロセッサ302とを有する。
イメージングシステム300は例えばマウス若しくはキーボードなどのユーザインターフェース装置を有し得、イメージングシステム300は第1画像及び/又は第2画像のスクロール可能な及び/又はズーム可能な視覚化を提供するように構成され得る。スクロール可能な及び/又はズーム可能な視覚化はユーザインターフェース装置を用いることによって入力コマンドに従って実行され得る。ユーザは第1画像及び/又は第2画像のスクロール可能及び/又はズーム可能な視覚化中に使用されるナビゲーション速度を変更し制御し得る。イメージングシステム300はデジタル画像処理において使用される様々な画像変換に基づいて第1画像及び/又は第2画像の視覚化を提供するように構成され得る。画像プロセッサ302を用いて、デジタル画像に対し画像処理を実行するアルゴリズムが使用される。アルゴリズムは入力データに適用され得、処理中のノイズ及び信号ひずみの蓄積などの問題を回避することができる。
異なる計算負荷能力が画像プロセッサ302による画像変換のために使用され得、使用される能力の最大閾値もイメージングシステム300によって定義され得、これは同じ処理ユニットのマルチユーザがいる場合に有利である。
図2は二つの免疫組織化学染色スライドの画像レジストレーションの例を示す。図2の左上部分において、両スライド間のミスアライメント誤差の例が示される。図2の左下部分において、完全な画像のレジストレーション結果の例が提示される。図2の右側において、両スライドにおけるマッチする視野の同時視覚化が提示される。
組織病理学画像の解析に加えて、分子法、例えばPCR若しくはシークエンシングも(腫瘍)組織の生物学的特徴を同定するために使用される。ポリメラーゼ連鎖反応(PCR)は単一若しくは数コピーのDNA断片を数桁を超えて増幅し、特定DNAシーケンスの数千から数百万コピーを生成する分子生物学における生化学技術である。
分子法では、組織の一部が追加の分子解析のために選択される。分子解析は選択された組織領域の厳密な組成に大きく依存し、関心領域(ROI)の選択の正確さが重要である。現在、ROIはH&E染色サンプルから手動で選択される。
MDxが次のステップとして選択される場合、一つ以上の連続サンプルが使用される。新たなサンプルを使用する理由の一つは、H&Eサンプルがガラスプレート間に密封されるという事実である。他の理由はH&Eサンプルを将来の参考と長期保管のために保全しておきたいという願望(法的な理由)及びH&E染色のMDxとの不適合性である。
典型的に、医療専門家は第2のサンプルを(カバースリップなしで)顕微鏡下に置き、先に検出されたROIとマッチする組織部分を見つけるためにサンプルをブラウズする。H&Eサンプルと第2のサンプル間の外観の違いに依存して、第2のサンプルにおけるROIの検索は煩雑になり得る。ROI選択ステップの正確さの重要性のために、このステップにおいてオペレータを支援する手段が有益である。
画像解析アルゴリズムは顕微鏡で取得された第2のサンプルのライブビデオストリームに対してH&Eサンプル中の選択されたROIの画像データをマッチさせるように開発され得る。この画像マッチングのプロセスは計算時間を要し、結果として画像パイプラインにおいて遅延を導入する。画像マッチングに加えて、画像マッチングの視覚化も計算時間を消費することになり、従って遅延を追加する。総遅延が小さいまま、例えば10Hzを上回る限りは、ビデオストリームに導入される遅延は許容可能である。サンプル間のより複雑な変換は典型的にはより長い計算時間をもたらす。サンプル間の変換の相違は大きくなる可能性があり、単純な並進と回転から成るが、組織の引き伸ばし、裂け目、及び折り目によって生じる非線形局所変形及びスケールのわずかな変動の可能性もある。
図3は本発明の一実施形態例にかかる構成要素の外観を示す。
図3は本発明にかかるアルゴリズム構造を示す。アルゴリズムの入力は組織サンプル1の画像("基準画像")と組織サンプル2の画像("標的画像")である。時間tにおける標的画像の画像コンテンツが、以前にキャプチャされ画像バッファに保存されている標的画像と比較される。本発明の一実施形態では、画像中のピクセルのサブセットに対して差分画像が計算される。
動き若しくは既定閾値を超える動き若しくは倍率変化がある場合、時間tにおける標的画像と時間t−1における標的画像は実質的に異なり、(完全)基準画像との画像の複雑なマッチングにつながり、非線形精緻化(非線形変形の推定だけでなく視覚化の両方)にとって不十分な時間を残す。従って、並進{Tx,Ty}及び回転{R}並びにスケールのみに関して基準画像と標的画像間の単純変換を記述する剛性変換パラメータを推定する。当初は、並進における変化のみについてレジストレーションを補正することが我々の意図である。しかしながら、変換パラメータの推定は回転が無視される場合不安定若しくは不正確になりやすい。
動きがない場合又は限られた計算時間のみが動き推定において失われる場合、両画像間の変換を記述するより正確な非線形変換パラメータを推定する時間が残る。標的画像の視覚化は並進及び回転で更新される。次のステップにおいて変換もその非線形性成分で図においてさらに精緻化される。
小さな動きの場合、小さな動きは既定閾値を下回る動きをあらわすが、変位は比較的単純な/効率的な動き推定法(例えば投影強度プラス1D相関に基づく)で推定されることができる。加えて、視野において依然観察される部分について結果を再利用することが可能である。これは先と同様、より詳細な変換の推定と視覚化の余地を残す。
ナビゲーション速度閾値若しくは動き閾値は、例えば毎秒0から5画像の範囲の間で予め定義され得るか、又は毎秒スクリーン高さ若しくはスクリーン幅単位で、例えばそれぞれ毎秒0から5スクリーン高さ若しくは幅の範囲の間で定義され得る。ナビゲーション速度閾値は動的に設定されることができる、すなわち閾値は固定値をあらわさなくてもよい。ナビゲーション速度閾値は利用可能な計算リソースに基づき得る。
構成要素B1‐1は基準画像tをあらわし、構成要素B1‐2は特徴抽出サブルーチンを象徴する。構成要素B2‐1は推定剛性特徴変換に対応し、構成要素B2‐2はレジストレーション{Tx,Ty,R}とスケールをあらわす。構成要素B3‐1は標的画像tをあらわし、構成要素B3‐2は特徴抽出を象徴する。構成要素B3‐3は画像バッファに対応し、構成要素B3‐4は計算されたSIFTベース動き推定若しくは任意の他の動き推定をあらわす。構成要素B3‐5は動きが存在するかどうか若しくはナビゲーション速度閾値が超えられるかどうかに関するクエリのシンボルである。
構成要素B4‐1は非線形変換の推定をあらわし、構成要素B4‐2は{Tx,Ty,R}の形のレジストレーション及びスケールを象徴する。構成要素B4‐3は{Tx,Ty,非線形}の形のレジストレーションに対応する。
図4は本発明の一実施形態例にかかる少なくとも二画像をレジストレーション及び視覚化するための方法の略フローチャート図を示す。
少なくとも二画像をレジストレーション及び視覚化する方法は以下のステップを有する。
方法の第1ステップとして、少なくとも二画像の第1画像のスクロール可能及び/又はズーム可能な視覚化中にユーザによって使用されるナビゲーション速度を測定し、ナビゲーション速度をナビゲーション速度閾値と比較するステップS1が実行される。
方法の第2のステップとして、ナビゲーション速度がナビゲーション速度閾値を上回る場合、第1レベルの計算負荷を用いて第1画像の変形を解くことによって第1画像を視覚化するステップS2が実行される。
方法の第3のステップとして、ナビゲーション速度がナビゲーション速度閾値を下回る場合、第2レベルの計算負荷を用いて第1画像の変形を解くことによって第1画像を視覚化するステップS3が実行される。
本発明の一実施形態例において、デジタルスキャンを基準として使用するカメラフィードのライブレジストレーションが、少なくとも二画像をレジストレーション及び視覚化するための方法によって実行される。
本発明の別の実施形態例において、二つ以上のデジタルスキャンの画像レジストレーションが、少なくとも二画像をレジストレーション及び視覚化するための方法によって実行される。
二つ以上のデジタルスキャンの場合、変位情報はユーザインタラクション装置(例えばデジタル環境におけるナビゲーションのために使用されるマウス、キーボード、若しくは他の装置)から直接利用可能であることがある。この場合は動き推定のための計算時間を要しないが、粗雑視覚化({Tx,Ty,R}のみ)と精密視覚化との差が、サーバの計算負荷を節約するために依然望ましい。
高解像度画像であり得る二画像間のレジストレーションを実行するために、以下のステップのサブセットが実行される。
画像中の局所特徴を検出して記述するコンピュータビジョンにおけるアルゴリズムであるSIFTを用いて特徴の計算が実行される。
変換の計算が、例えばSIFT特徴マッチングに基づいて、又は投影強度の相関に基づいて実行され、これは比較的小さな変位、例えば閾値未満の変位に適用可能であり得る。局所変形の計算は例えば局所エッジ構造のマッチングによって実行される。
画像の変換は並進と回転のみを用いて実行される。これはスケールによっても実行され得る。追加ステップとして、少なくとも二画像のうち一つの局所変形が実行され得る。
ユーザ経験が依然快適であるようなフレームレート若しくは要求視点のビューイング若しくはローディング遅延によって制約され得る、利用可能な処理される計算負荷及び/又は時間量に依存して、上述のステップを用いる異なるストラテジが適用可能である。これらのステップは混合されるか若しくは同時に処理され得る。
本発明の別の実施形態例において、大きな変位の場合、局所変形も適用するのではなく粗い概観が提示される。従って、局所変形の計算及び少なくとも二画像のうち一つの局所変形は無視されることができ、粗いが許容可能なレジストレーション、限られた遅延で許容可能な高フレームレートをもたらす。
本発明の別の実施形態例において、小さな変位、例えば変位閾値未満の場合、例えば投影強度の相関に基づく、動き推定のためにより効率的なアルゴリズムを用いることによって、又は前のフレームと重なるエリアからの情報を再利用することによって、時間若しくは計算負荷が節約され得る。そして、節約された時間は局所変形を計算するため及び局所変形の視覚化のために使用されることができる。
本発明の別の実施形態例において、アクティブなパン及びズームが適用される。ユーザがアクティブにパン及びズームしている場合、並進と回転の補正が実行される。パン及びズームにアクティブでありながら、レジストレーションの非線形及び/又は精密な詳細は必要なく、省略され得る。
本発明の別の実施形態例において、パン及び/又はズームは起こらず、例えばユーザはパン及び/又はズームしておらず、画像は並進と回転を用いて表示される。そして、局所変形の計算が開始される。一旦、局所変形画像が利用可能になると、視点が更新される。
本発明の別の実施形態例において、適切なシステム上で先の実施形態の一つにかかる方法の方法ステップを実行するのに適していることを特徴とするコンピュータプログラム若しくはコンピュータプログラム素子が提供される。
従ってコンピュータプログラム素子は本発明の一実施形態の一部でもあり得るコンピュータユニット上に保存され得る。このコンピュータユニットは上記方法のステップを実行するか若しくはその実行を誘導するように適応され得る。
さらに、これは上記装置の構成要素を操作するように適応され得る。コンピュータユニットは自動的に作動する、及び/又はユーザの命令を実行するように適応され得る。コンピュータプログラムはデータプロセッサのワーキングメモリへロードされ得る。従ってデータプロセッサは本発明の方法を実行するように装備され得る。
本発明のこの実施形態例は初めから本発明を使用するコンピュータプログラムと、アップデートにより既存のプログラムを本発明を使用するプログラムへ変えるコンピュータプログラムの両方をカバーする。
さらに、コンピュータプログラム素子は上記方法の一実施形態例の手順を満たす全所要ステップを提供することができてもよい。
本発明のさらなる実施形態例によれば、前章によって記載されるコンピュータプログラム素子を記憶した、CD‐ROMなどのコンピュータ可読媒体が提示される。
コンピュータプログラムは、他のハードウェアと一緒に若しくはその一部として供給される光学記憶媒体若しくはソリッドステート媒体などの適切な媒体上に記憶及び/又は分散され得るが、インターネット又は他の有線若しくは無線通信システムなどを介して他の形式で分散されてもよい。
しかしながら、コンピュータプログラムはワールドワイドウェブのようなネットワークを介して提示されてもよく、かかるネットワークからデータプロセッサのワーキングメモリへダウンロードされることができる。
本発明のさらなる実施形態例によれば、コンピュータプログラム素子をダウンロードに利用可能にするための媒体が提供され、かかるコンピュータプログラム素子は本発明の前記実施形態の一つにかかる方法を実行するように構成される。
本発明の実施形態は異なる主題に関して記載されていることが留意されなければならない。特に、一部の実施形態は方法タイプクレームに関して記載されるが、他の実施形態は装置タイプクレームに関して記載される。
しかしながら、当業者は上記及び下記から、他に明記されない限り、一タイプの主題に属する特徴の任意の組み合わせに加えて、異なる主題に関する特徴間の任意の組み合わせも本願とともに開示されるとみなされることがわかるだろう。しかしながら、全特徴は組み合わされて特徴の単なる総和以上の相乗効果をもたらし得る。
本発明は図面と先の説明において詳細に図示され記載されているが、かかる図示と記載は例示若しくは説明であって限定ではないとみなされるものとする。本発明は開示の実施形態に限定されない。開示の実施形態への他の変更は図面、開示及び添付のクレームの考察から、請求される発明を実施する当業者によって理解されもたらされることができる。
クレーム中、"有する"という語は他の要素若しくはステップを除外せず、不定冠詞"a"若しくは"an"は複数を除外しない。単一プロセッサ若しくはコントローラ若しくは他のユニットはクレームに列挙される複数の項目の機能を満たし得る。特定の手段が相互に異なる従属クレームに列挙されるという単なる事実はこれら手段の組み合わせが有利に使用されることができないことを示さない。クレーム中の任意の参照符号は範囲を限定するものと解釈されてはならない。

Claims (10)

  1. 少なくとも二画像をレジストレーション及び視覚化するための方法であって、
    前記少なくとも二画像からの第1画像のスクロール可能及び/又はズーム可能な視覚化中にユーザによって使用されるナビゲーション速度を測定し、前記ナビゲーション速度をナビゲーション速度閾値と比較するステップと、
    前記ナビゲーション速度が前記ナビゲーション速度閾値を下回る場合、第1レベルの計算負荷を用いて前記第1画像の変形を解くことによって前記第1画像を視覚化するステップと、
    前記ナビゲーション速度が前記ナビゲーション速度閾値を上回る場合、第2レベルの計算負荷を用いて前記第1画像の変形を解くことによって前記第1画像を視覚化するステップと
    を有し、前記方法が前記少なくとも二画像の局所レジストレーションを実行するステップをさらに有する、方法。
  2. 前記第1画像の変形を解くステップが、前記第1画像の画像並進、前記第1画像の画像回転、前記第1画像のスケール変動、前記第1画像の追加線形変形若しくは前記第1画像の非線形変形を有する、請求項1に記載の方法。
  3. 前記第1画像のスクロール可能及び/又はズーム可能な視覚化中に前記ユーザによって使用される前記ナビゲーション速度がユーザインターフェース装置から取得される、請求項1又は2に記載の方法。
  4. 前記第1画像のスクロール可能及び/又はズーム可能な視覚化中に前記ユーザによって使用される前記ナビゲーション速度がスケール不変特徴変換特徴マッチングを用いて測定される、請求項1又は2に記載の方法。
  5. 前記ナビゲーション速度が前記ナビゲーション速度閾値を上回る場合、前記第1画像の変形を解くステップが並進及び回転変換を適用するステップを有し、前記ナビゲーション速度が前記ナビゲーション速度閾値を下回る場合、前記第1画像の変形を解くステップが非線形精緻化を適用するステップを有する、請求項1から4のいずれか一項に記載の方法。
  6. 前記ナビゲーション速度が前記ナビゲーション速度閾値を上回る場合、前記第1画像の変形を解くステップが第1解像度レベルで変換を適用するステップを有し、前記ナビゲーション速度が前記ナビゲーション速度閾値を下回る場合、前記第1画像の変形を解くステップが第2解像度レベルで変換を適用するステップを有し、前記第1解像度レベルが前記第2解像度レベルより低い、請求項1から5のいずれか一項に記載の方法。
  7. プロセッサによって実行されるとき、請求項1から6のいずれか一項に記載の方法のステップを実行するコンピュータプログラム。
  8. 請求項7に記載のコンピュータプログラムが格納された、コンピュータ可読媒体。
  9. 請求項1から6のいずれか一項に記載の方法を実行するように構成されるプロセッサ。
  10. デジタル画像レジストレーション及び視覚化のためのイメージングシステムであって、
    ナビゲーション速度をナビゲーション速度閾値と比較するように構成される画像プロセッサを有し、前記ナビゲーション速度は少なくとも二画像の第1画像のスクロール可能及び/又はズーム可能な視覚化中にユーザによって使用され、
    前記画像プロセッサがさらに、前記ナビゲーション速度が前記ナビゲーション速度閾値を下回る場合、第1レベルの計算負荷を用いて前記第1画像の変形を解くことによって前記第1画像を視覚化するように構成され、
    前記画像プロセッサがさらに、前記ナビゲーション速度が前記ナビゲーション速度閾値を上回る場合、第2レベルの計算負荷を用いて前記第1画像の変形を解くことによって前記第1画像を視覚化するように構成され、
    前記画像プロセッサがさらに、前記少なくとも二画像の局所レジストレーションを実行するように構成され、
    前記少なくとも二画像を表示するように構成されるディスプレイを有する、
    イメージングシステム。
JP2016549165A 2014-02-04 2015-01-29 少なくとも二画像をレジストレーション及び視覚化するための方法 Active JP6573894B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP14153744 2014-02-04
EP14153744.9 2014-02-04
PCT/EP2015/051745 WO2015117875A1 (en) 2014-02-04 2015-01-29 Method for registering and visualizing at least two images

Publications (3)

Publication Number Publication Date
JP2017506384A true JP2017506384A (ja) 2017-03-02
JP2017506384A5 JP2017506384A5 (ja) 2018-04-19
JP6573894B2 JP6573894B2 (ja) 2019-09-11

Family

ID=50064450

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016549165A Active JP6573894B2 (ja) 2014-02-04 2015-01-29 少なくとも二画像をレジストレーション及び視覚化するための方法

Country Status (5)

Country Link
US (1) US9972068B2 (ja)
EP (1) EP3103100B1 (ja)
JP (1) JP6573894B2 (ja)
CN (1) CN106133787B (ja)
WO (1) WO2015117875A1 (ja)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10198590B2 (en) 2015-11-11 2019-02-05 Adobe Inc. Content sharing collections and navigation
US10783431B2 (en) 2015-11-11 2020-09-22 Adobe Inc. Image search using emotions
US10249061B2 (en) 2015-11-11 2019-04-02 Adobe Inc. Integration of content creation and sharing
US20170132252A1 (en) * 2015-11-11 2017-05-11 Adobe Systems Incorporated Content Sharing Collections and Navigation
US10389804B2 (en) 2015-11-11 2019-08-20 Adobe Inc. Integration of content creation and sharing
WO2020037570A1 (zh) * 2018-08-22 2020-02-27 深圳市真迈生物科技有限公司 图像配准方法、装置和计算机程序产品

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000270169A (ja) * 1999-03-05 2000-09-29 Hewlett Packard Co <Hp> 運動識別方法
JP2013530467A (ja) * 2010-06-15 2013-07-25 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 顕微鏡検査法の画像処理方法
JP2013152452A (ja) * 2011-12-27 2013-08-08 Canon Inc 画像処理装置、画像表示システム、画像処理方法および画像処理プログラム
JP2013200640A (ja) * 2012-03-23 2013-10-03 Canon Inc 画像処理装置、画像処理システム、画像処理方法、およびプログラム

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6711283B1 (en) * 2000-05-03 2004-03-23 Aperio Technologies, Inc. Fully automatic rapid microscope slide scanner
US6859199B2 (en) * 2001-11-06 2005-02-22 Omnivision Technologies, Inc. Method and apparatus for determining relative movement in an optical mouse using feature extraction
US20040047497A1 (en) * 2002-09-10 2004-03-11 Confirma, Inc. User interface for viewing medical images
US7266249B2 (en) 2003-05-22 2007-09-04 Ge Medical Systems Global Technology Company, Llc Optimized region growing algorithm for scale space analysis
JP2007501470A (ja) * 2003-05-27 2007-01-25 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 多数の制御機能を備える診断撮像システムの制御部
JP3924558B2 (ja) * 2003-11-17 2007-06-06 富士通株式会社 生体情報採取装置
CN101006733B (zh) * 2004-08-18 2010-05-05 夏普株式会社 图像数据显示装置
JP2008538080A (ja) * 2005-02-03 2008-10-09 ブラッコ イメージング ソチエタ ペル アチオニ 対象物の移動による画像アーチファクトの減少を伴う生物医学画像の登録のための方法と、コンピュータプログラムプロダクト
US8681104B2 (en) * 2007-06-13 2014-03-25 Apple Inc. Pinch-throw and translation gestures
US8063901B2 (en) * 2007-06-19 2011-11-22 Siemens Aktiengesellschaft Method and apparatus for efficient client-server visualization of multi-dimensional data
JP4752900B2 (ja) * 2008-11-19 2011-08-17 ソニー株式会社 画像処理装置、画像表示方法および画像表示プログラム
US20100157079A1 (en) 2008-12-19 2010-06-24 Qualcomm Incorporated System and method to selectively combine images
US8437577B2 (en) * 2009-03-05 2013-05-07 Tektronix, Inc. Methods and systems for image registration
US8274592B2 (en) * 2009-12-22 2012-09-25 Eastman Kodak Company Variable rate browsing of an image collection
WO2012037417A1 (en) * 2010-09-16 2012-03-22 Omnyx, LLC Control configuration for digital image system
US8699769B2 (en) 2011-07-12 2014-04-15 Definiens Ag Generating artificial hyperspectral images using correlated analysis of co-registered images
JP5726706B2 (ja) * 2011-10-14 2015-06-03 クラリオン株式会社 ナビゲーション装置
JP5792607B2 (ja) * 2011-12-09 2015-10-14 株式会社ソニー・コンピュータエンタテインメント 画像処理装置および画像処理方法
US8755633B2 (en) * 2012-03-14 2014-06-17 Sony Corporation Automated synchronized navigation system for digital pathology imaging

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000270169A (ja) * 1999-03-05 2000-09-29 Hewlett Packard Co <Hp> 運動識別方法
JP2013530467A (ja) * 2010-06-15 2013-07-25 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 顕微鏡検査法の画像処理方法
JP2013152452A (ja) * 2011-12-27 2013-08-08 Canon Inc 画像処理装置、画像表示システム、画像処理方法および画像処理プログラム
JP2013200640A (ja) * 2012-03-23 2013-10-03 Canon Inc 画像処理装置、画像処理システム、画像処理方法、およびプログラム

Also Published As

Publication number Publication date
JP6573894B2 (ja) 2019-09-11
US20170011489A1 (en) 2017-01-12
US9972068B2 (en) 2018-05-15
CN106133787B (zh) 2020-09-01
CN106133787A (zh) 2016-11-16
EP3103100A1 (en) 2016-12-14
EP3103100B1 (en) 2019-11-06
WO2015117875A1 (en) 2015-08-13

Similar Documents

Publication Publication Date Title
JP6573894B2 (ja) 少なくとも二画像をレジストレーション及び視覚化するための方法
US11232565B2 (en) Examining device for processing and analyzing an image
US8229194B2 (en) Feature-based registration of sectional images
Glaser et al. Stereology, morphometry, and mapping: the whole is greater than the sum of its parts
JP4071186B2 (ja) 生物標本内の関心対象を識別するための方法及びシステム
EP1978486B1 (en) Method of image analysis based on chromogen separation
KR20150095646A (ko) 바이오마커 발현의 선택 및 디스플레이
JP6960935B2 (ja) コントロールスライドを使用する改善された画像解析アルゴリズム
US10685210B2 (en) Tissue microarray registration and analysis
US9424640B2 (en) Pathological diagnosis support apparatus and pathological diagnosis support method
JP6646552B2 (ja) 画像診断支援装置、画像診断支援方法及び試料分析システム
JP2024112965A (ja) 画像解析装置
JP6624369B2 (ja) 画像処理装置、画像処理方法、画像表示装置およびコンピュータプログラム
JP4920507B2 (ja) 画像処理装置および画像処理プログラム
Perdigão et al. Okapi-EM: A napari plugin for processing and analyzing cryogenic serial focused ion beam/scanning electron microscopy images
US20180025211A1 (en) Cell tracking correction method, cell tracking correction device, and storage medium which stores non-transitory computer-readable cell tracking correction program
JP2008188012A (ja) 分子組織学
Arlt et al. One size fits all: evaluation of the transferability of a new “learning” histologic image analysis application
JP5910637B2 (ja) 生体画像解析システム、生体画像解析方法および生体画像解析プログラム
WO2015089564A1 (en) Thickness estimation for microscopy
Aziz et al. Effectiveness of color correction on the quantitative analysis of histopathological images acquired by different whole-slide scanners
Perdigão et al. Okapi-EM–a napari plugin for processing and analysing cryogenic serial FIB/SEM images
EP4207071A1 (en) Systems and methods for measuring proximity of objects across serial sections of overlayed stained slide pathology images

Legal Events

Date Code Title Description
RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20170214

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20171108

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20171108

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180307

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20180307

A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20180607

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20180607

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180612

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180907

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180918

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20181217

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20190402

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190717

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

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20190725

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190814

R150 Certificate of patent or registration of utility model

Ref document number: 6573894

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