JP2022534157A - 組織学的画像および術後腫瘍辺縁評価における腫瘍のコンピュータ支援レビュー - Google Patents

組織学的画像および術後腫瘍辺縁評価における腫瘍のコンピュータ支援レビュー Download PDF

Info

Publication number
JP2022534157A
JP2022534157A JP2021538332A JP2021538332A JP2022534157A JP 2022534157 A JP2022534157 A JP 2022534157A JP 2021538332 A JP2021538332 A JP 2021538332A JP 2021538332 A JP2021538332 A JP 2021538332A JP 2022534157 A JP2022534157 A JP 2022534157A
Authority
JP
Japan
Prior art keywords
tumor
image
histological
tissue
dye
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
JP2021538332A
Other languages
English (en)
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.)
Leica Biosystems Imaging Inc
Original Assignee
Leica Biosystems Imaging Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Leica Biosystems Imaging Inc filed Critical Leica Biosystems Imaging Inc
Publication of JP2022534157A publication Critical patent/JP2022534157A/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • 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
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • G06F18/2148Generating training patterns; Bootstrap methods, e.g. bagging or boosting characterised by the process organisation or structure, e.g. boosting cascade
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/94Hardware or software architectures specially adapted for image or video understanding
    • G06V10/95Hardware or software architectures specially adapted for image or video understanding structured as a network, e.g. client-server architectures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/60Type of objects
    • G06V20/69Microscopic objects, e.g. biological cells or cellular parts
    • G06V20/693Acquisition
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/60Type of objects
    • G06V20/69Microscopic objects, e.g. biological cells or cellular parts
    • G06V20/695Preprocessing, e.g. image segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/60Type of objects
    • G06V20/69Microscopic objects, e.g. biological cells or cellular parts
    • G06V20/698Matching; Classification
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/40ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
    • 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/04842Selection of displayed objects or displayed text elements
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00Indexing scheme for image data processing or generation, in general
    • G06T2200/24Indexing scheme for image data processing or generation, in general involving graphical user interfaces [GUIs]
    • 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/20Special algorithmic details
    • G06T2207/20081Training; Learning
    • 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/20084Artificial neural networks [ANN]
    • 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/30024Cell structures in vitro; Tissue sections in vitro
    • 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/30096Tumor; Lesion
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V2201/00Indexing scheme relating to image or video recognition or understanding
    • G06V2201/03Recognition of patterns in medical or anatomical images

Abstract

組織学的画像において腫瘍を同定および可視化し、腫瘍辺縁を測定するためのコンピュータ装置および方法が提供される。CNNを使用して、画像内のピクセルは、非腫瘍組織に関連すると判別されたのか、または腫瘍組織の1つ以上のクラスに関連すると判別されたのかに従って分類される。CNNの結果に基づいてセグメンテーションが実行され、マスクが生成され、当該マスクにより、個々の腫瘍によって占有された領域がマーキングされる。各腫瘍の要約統計量が計算され、フィルタに提供される。フィルタにより、重要でないとみなされた腫瘍を除外することによってセグメンテーションマスクが編集される。任意に、フィルタを通過する腫瘍は、要約統計量に従って、例えば臨床的関連性の順序で、または病理医にとって常識的なレビューの順序でランク付けすることができる。次に、ビジュアライゼーションアプリケーションにより、セグメンテーションマスク、要約統計量、および/またはランク付けを考慮して組織学的画像を表示することができる。切除によって抽出された腫瘍塊は、その表面領域を強調表示するために染料で着色される。CNNは、染料組織および非染料組織、ならびに腫瘍組織および非腫瘍組織を区別するように訓練されている。CNNは組織学的画像に適用され、ピクセルが組織クラスに割り当てられている出力画像を生成する。組織切片の腫瘍辺縁状態は、腫瘍染料分類ピクセルの有無によって判別される。腫瘍の辺縁病変および腫瘍辺縁距離は、分類で指定されたピクセル間距離パラメータに基づいて追加のパラメータを計算することによって求められる。

Description

本開示は、切除によって患者から除去された腫瘍塊の辺縁を評価する病理医によるレビューおよび診断を支援するために、畳み込みニューラルネットワーク(CNN)を用いた、腫瘍を含む組織学的画像の画像処理に関する。
デジタルパソロジーにより、病理医がスライドを観察して診断する手法は変化を続けている。従来の手法では、病理医は、スライドガラスを顕微鏡下で観察することによりスライドを検査する。病理医はまず、低倍率の対物レンズでスライドを観察する。診断価値がある可能性が高い領域が観察された場合、病理医は、高倍率の対物レンズに切り替えて、その領域をより詳細に観察する。その後、病理医は低倍率に戻し、スライド上の他の領域の検査を続ける。スライドに対する確定的かつ完全な診断が可能となるまで、低倍率-高倍率-低倍率の観察シーケンスがスライド上で数回繰り返される場合がある。
過去20年間で、デジタルスキャナの導入によりこのワークフローが変化している。デジタルスキャナは、スライドガラス全体の画像、いわゆる全スライド画像(WSI)を取得し、病理医を必要としない、大部分が自動化されたプロセスにおいて、デジタル画像データファイルとしてWSIを保存することができる。得られた画像データファイルは、通常、スライドデータベースに保存される。高解像度ディスプレイを備えた観察ワークステーションにいる病理医は、スライドデータベースからクリニカルネットワークを介して、これを利用することができる。当該ワークステーションは、この目的のためのビジュアライゼーションアプリケーションを有している。
病理学におけるより最近の進展として、CNN法が研究上の関心を高めていることが挙げられる。組織学的画像から腫瘍を同定し診断する際に、CNN法が、病理医と同等のまたはそれ以上の性能を発揮していることが報告されるようになってきている。
Wang et al.(2016)では、乳癌のリンパ節への転移を検出するためのCNNアプローチについて説明されている。
米国特許出願公開第2015/213302号明細書では、癌性組織の領域で細胞の有糸分裂がどのように検出されるかが説明されている。CNNを訓練した後、有糸分裂カウントを実行する自動核検出システムに基づいて分類が実行され、有糸分裂カウントは、腫瘍の等級付けに使用されている。
Hou et al.(2016)では、脳癌および肺癌の画像処理が行われている。WSIの画像パッチが、パッチレベルCNNによって与えられるパッチレベルの予測を行うために使用されている。
Liu et al.(2017)では、ギガピクセルの乳癌の組織学的画像から抽出された画像パッチをCNNで処理し、画像内の全てのピクセルに腫瘍確率を割り当てることで、腫瘍が検出され、その位置が特定されている。
Bejnordi et al.(2017)では、2つの積層CNNを適用して、ヘマトキシリンおよびエオシン(H&E)染色で染色された乳房組織のWSIから抽出された画像パッチの腫瘍が分類されている。その性能は、これらの病理画像での物体の検出およびセグメンテーションに適していることを示している。さらに、Bejnordi et al.は、乳癌サンプルに適用される他のCNNベースの腫瘍分類方法の概要も提供していることに留意されたい(参考文献10~13を参照)。
Esteva et al.(2017)では、深層CNNを適用して皮膚病変が分析され、木構造分類に従って、病変が種々の悪性型、非悪性型および非腫瘍性型に分類されている。これには、悪性型の末端黒子型黒色腫、無色素性黒色腫および黒子黒色腫、ならびに非悪性型の青色母斑、ハロー母斑および蒙古斑が含まれている。皮膚病変(黒色腫など)の画像が、臨床クラス全体の確率分布に順次ワーピングされて、分類が実行される。
Mobadersany et al.(2017)では、脳腫瘍と診断された患者の全生存率を予測する生存CNNに基づく計算方法が開示されている。患者の転帰を予測するために、組織生検による病理画像データ(組織学的画像データ)がモデルおよび患者固有のゲノムバイオマーカーに入力されている。この方法では、適応フィードバックを使用して、患者の転帰に関連する視覚パターンおよび分子バイオマーカーが同時に学習される。
腫瘍切除術は、通常、臓器または腺の一部または全部を切除することによって腫瘍組織を切除する手術であり、部分的な切除術は、部分切除術と称される。腫瘍切除の目的は、腫瘍の辺縁の周囲の少量の正常組織とともに全ての腫瘍を含む組織塊を除去することである。
切除により腫瘍が完全に切除されたことを確認するために、身体から組織塊を切除した後、以下の病理検査が行われる。切除された組織塊の表面に着色染料で着色される。微視的な組織構造を保存するために組織塊が種々の方法で準備される。次に、組織塊がミクロトームで薄片にスライスされる。スライスされた後、微視的組織構造のさらなる保存措置を実施することができる。組織切片がスライドにマウントされる。組織切片は、概して、特殊な染色またはマーカーで染色され、例えばコントラストによって、細胞の視認性が高められ、あるいは特定の遺伝子発現または関連する受容体に固有のマーカーで、特定の種類の細胞が同定される。乳癌腫瘍の場合、標的となる典型的な受容体は、ER+、ER-、PR+、PR-、HER2+、およびHER2-である。ここで、ERはエストロゲン受容体を表し、PRはプロゲステロン受容体を表し、HER2はヒト上皮成長因子受容体2を表す。
従来の辺縁の評価方法では、病理医が染色されたスライドを顕微鏡下で観察し、いずれかの腫瘍細胞が染料表面上にあるか、または染料表面に過度に近接しているかどうかを判別している。染料上に腫瘍細胞が見つかった場合、病理医は、その外科的切除辺縁が陽性であると報告する(Lester et al.(2009))。陽性の辺縁は、通常、より多くの組織塊を除去するためにさらなる手術が必要であることを示す。陽性の辺縁の場合、病理医は、辺縁への湿潤の量に関連する、巣性分類、最小/中程度分類、および広範囲分類などの下位分類を使用して、辺縁病変の程度を報告することができる。癌細胞が染料表面に近接しているが、完全に染料表面の上にはない場合、例えば2mm未満で近接している場合、病理医は腫瘍の辺縁を「クロース」と報告することができる。染料表面に癌細胞がない(「陽性」ではない)場合、またはその近傍にない(「クロース」ではない)場合、病理医は辺縁を「クリア」として報告する。この場合、例えば、辺縁は2mm以上である。クリアな辺縁またはクロースな辺縁の場合、病理医は、染料表面と染料表面に最も近接した腫瘍細胞との間の距離を示す辺縁距離を報告することができる。辺縁が陽性の場合、残った癌性組織を除去するためにさらに手術を行う必要がある場合がある。辺縁状態を正確に報告することはその後の癌治療にとって重要であるが、辺縁状態の報告の質には大きなばらつきがある(Persing et al.(2011))。
このため、上記のような従来のシステムに見られる重大な問題を克服したシステムおよび方法が求められている。
本開示の一態様によれば、組織学的画像(またはそのセット)内の腫瘍を同定する方法が提供され、この方法は、
ピクセルの2次元アレイを含む組織学的画像を受信することと、
畳み込みニューラルネットワーク(CNN)を適用して、組織学的画像の2次元アレイに対するマッピングを有するピクセルの2次元アレイを有する出力画像を生成することであって、出力画像が、複数の組織クラスのうちの1つを各ピクセルに割り当てることによって生成され、複数の組織クラスが、非腫瘍組織を表す少なくとも1つのクラスおよび腫瘍組織を表す少なくとも1つのクラスを含む、ことと、
個々の腫瘍によって占有された領域がマーキングされている出力画像から、セグメンテーションマスクを生成することと、
各腫瘍の要約統計量を計算することと、
各腫瘍の要約統計量にフィルタを適用して、フィルタに従って腫瘍を選択および選択解除することにより、セグメンテーションマスクを編集することと、
を含む。
提案しているアプローチにより、画像内の潜在的なまたは潜在的に重要な腫瘍がどこにあるかを病理医に示すことができ、これにより、病理医に、これらの領域への直接のナビゲートを誘導することが可能となる。このようにして、腫瘍を含む領域が病理医によって見落とされる可能性が低くなるため、分析の信頼性を向上させることができる。さらに、病理医が腫瘍領域を発見するためにスライド領域全体を検査する必要がなくなるため、スループットを向上させることができる。病理医はまた、CNN分析結果に基づいて行われたフィルタリングおよびランク付けを考慮に入れて、限られた分析時間を配分するように暗黙的に支援される。これにより、特定のスライド画像の調査における病理医の時間がより良好に全体的に配分される。
この方法は、スコアリングアルゴリズムに従って各腫瘍をスコアリングして、各腫瘍にスコアを割り当てることによって拡張することができる。このスコアに基づいて、腫瘍をランク付けすることができる。ランク付け順は、多種多様な有用な順序を表すことができ、例えば、知覚された臨床的関連性を示してよく、臨床医が腫瘍のレビューを希望する順序を予測してもよい。ランク付けは、フィルタを適用した後に残った腫瘍に限定されうる。ランク付けは、組織学的画像に関連付けられた追加のメタデータとして、レコードに保存されうる。組織学的画像は、仮想スライドライブラリもしくは他のデータベースに保存されたレコードから受信することができ、またはデジタルスキャナによる取得の直後に受信することができ、この場合、この方法は、デジタルスキャナによって、または例えばデジタルスキャナに接続されているコンピュータ装置によって、または別のデータソースから、実行可能である。
次に、編集されたセグメンテーションマスクを考慮して、組織学的画像のビジュアライゼーションを作成することができる。次に、ビジュアライゼーションは、コンピュータ装置上で実行されているビジュアライゼーションアプリケーションによって、ディスプレイ上でユーザに表示されうる。例えば、ビジュアライゼーションには、セグメンテーションマスクが組織学的画像上にオーバレイされる概観表示ペインが含まれうる。ビジュアライゼーションの別の選択肢は、セグメンテーションマスクおよび組織学的画像が1対1の比較のために、例えば並べて、互いに隣接して提示されるそれぞれの概観表示ペインを含むことである。概観表示ペインには、ランク付けに基づいた各腫瘍のランク付けラベルが含まれうる。ユーザがビジュアライゼーションとのインタラクションを行って、編集されたセグメンテーションマスク内の腫瘍を選択することを許可するために、ユーザインタフェース腫瘍選択制御が提供されうる。ビジュアライゼーションには、現在選択されている腫瘍にズームインされた拡大表示ペインが含まれうる。腫瘍選択制御は、編集されたセグメンテーションマスク内の腫瘍をランク付け順でスイープするためのスクロール機能を組み込むことができる。腫瘍選択制御はまた、ユーザが現在選択されている腫瘍に対して追加の計算診断プロセスを実行することを可能にするための診断機能を組み込むことができる。
最も単純な形式では、CNNには2つの組織クラスしか存在しない。1つは腫瘍組織のクラスであり、もう1つは非腫瘍組織のクラスである。ただし、2つ以上の腫瘍組織のクラスが提供されてもよい。例えば、腫瘍組織の組織クラスでは、浸潤性腫瘍と非浸潤性(in situ)腫瘍とが区別されてもよい。非腫瘍組織の組織クラスには、1つの組織クラスしかない場合もあれば、例えば骨と他の組織とを区別するために複数のクラスがある場合もある。また、非組織のクラス、すなわちスライド画像全体の中でサンプルが存在しない領域のクラスがあってもよく、これは、特に組織マイクロアレイサンプルに有効でありうる。
この方法は、割り当てステップの後に、出力画像パッチを組織学的画像(またはそのセット)の確率マップにアセンブルするステップをさらに含みうる。
この方法は、アセンブルするステップの後に、確率マップを組織学的画像(またはそのセット)にリンクするように、例えばデータリポジトリ内のレコードに確率マップを保存するステップをさらに含みうる。
本実装形態において、連続する各畳み込み段階では、次元が削減されるにつれて深度が増大するため、畳み込み層の深度は増大し、次元は削減される。また、連続する転置畳み込み段階では、次元が増大するにつれて深度が減少するため、逆畳み込み層の深度は減少し、次元は増大する。最終畳み込み層は、最大の深度および最小の次元を有する。畳み込み段階および逆畳み込み段階でそれぞれ深度が増大および減少するアプローチの代わりに、入力層と出力層とを除く全ての層が同じ深度を有するニューラルネットワークを設計することもできる。
この方法はさらに、組織学的画像(またはそのセット)を確率マップとともにディスプレイ上に表示すること、例えば、その上にオーバレイして、または互いに並べて表示することを含みうる。確率マップを使用して、使用する免疫組織化学(IHC)スコアリングアルゴリズムがどんなものであれ、これによってスコアリングすべき領域を判別できる。確率マップを使用して、例えば、病理医がCNNによって生成された結果を評価するために、ディスプレイに表示されうる腫瘍細胞の周囲の等高線のセットを生成することもできる。
特定の実施形態では、畳み込みニューラルネットワークは、1つ以上のスキップ接続を有する。各スキップ接続は、畳み込み層のうちの最終畳み込み層よりも大きな次元の少なくとも1つの畳み込み層から中間結果を取得し、当該結果を、なし、1つ、または2つ以上でありうる必要な数の転置畳み込みに供して、入力画像パッチにサイズが一致する少なくとも1つのさらなる復元された層を取得する。次に、これらは、各ピクセルに組織クラスを割り当てる該ステップの前に、上記の復元された層と組み合わされる。さらなる処理ステップでは、確率を再計算するために、復元された層をさらなる復元された層の各々と組み合わせ、これによって、スキップ接続から得られた結果が考慮される。
特定の実施形態では、ソフトマックス演算を使用して確率を生成する。
組織学的画像から抽出された画像パッチは、画像の全領域をカバーすることができる。パッチは、重複しない画像タイル、またはマージン部分で重複して確率マップのつなぎ合わせを支援する画像タイルでありうる。CNNは固定サイズのピクセルアレイのみを受け入れるように設計されているため、各画像パッチは、CNNに一致する幅および高さの固定数のピクセルを有する必要があるが、このことは、各画像パッチが組織学的画像上の同じ物理的領域に対応する必要があることを意味するものではない。なぜなら、組織学的画像のピクセルは、より広い領域をカバーする低解像度のパッチに結合される可能性があるからである。例えば、隣接するピクセルの各2×2アレイを組み合わせて、1つの「スーパー」ピクセルとし、組織学的画像のネイティブ解像度で抽出されたパッチの物理的領域の4倍のパッチを形成することができる。
本方法は、CNNが訓練された後、予測のために実行することができる。この訓練の目的は、層間接続に適切な重み値を割り当てることである。訓練のために、組織学的画像(またはそのセット)の各ピクセルを組織クラスの1つに割り当てるグラウンドトゥルースデータが提供される。グラウンドトゥルースデータは、専門の臨床医が十分な数の画像に注釈を付けることにより得られる。訓練は、CNNを繰り返し適用することによって実行され、各反復では、グラウンドトゥルースデータと出力画像パッチとの比較に基づいて重み値が調整される。本実装形態では、重みは、勾配降下法による訓練中に調整される。
組織クラスの設定には種々の選択肢があるが、全ての実施形態でないにしても、大部分の実施形態では、非腫瘍組織と腫瘍組織との間でクラスが区別されることが共通している。非腫瘍組織クラスは、1つ、2つ、またはそれ以上のクラスを含みうる。腫瘍組織クラスはまた、1つ、2つ、またはそれ以上のクラスを含みうる。例えば、本実装形態では、3つの組織クラスがある。1つは非腫瘍組織クラス、2つは腫瘍組織クラスである。2つの腫瘍組織クラスは、浸潤性腫瘍クラスおよび非浸潤性腫瘍クラスである。
幾つかの実施形態では、CNNは、一度に1つの組織学的画像に適用される。他の実施形態では、CNNは、組織の領域の異なって染色された隣接するセクションから取得された組織学的画像のセットを組み合わせることによって形成された合成組織学的画像に適用されうる。さらに別の実施形態では、CNNは、組織の領域の異なって染色された隣接するセクションから取得された画像のセットの画像の各々と並行して適用されうる。
CNNの結果を用いて、本方法は、ピクセル分類および確率マップを参照してその分類から定義される腫瘍に基づくスコアリングプロセスを含むように拡張されうる。例えば、本方法は、確率マップに従って、腫瘍に対応する組織学的画像内の領域を画定することと、スコアリングアルゴリズムに従って各腫瘍をスコアリングして、各腫瘍にスコアを割り当てることと、スコアを、例えばデータリポジトリのレコードに保存することと、をさらに含みうる。したがって、スコアリングは、組織学的画像上で行われるが、確率マップによって腫瘍組織を含むと特定された領域に限定される。
結果は、臨床医のディスプレイに表示されうる。すなわち、組織学的画像は、関連する確率マップとともに表示することができ、例えば、その上にオーバレイして、または互いに並べて表示することができる。腫瘍スコアは、例えば、腫瘍にテキストラベルを付けたり、腫瘍を指し示したり、画像の横に表示したりするなど、便利な方法で表示することもできる。
畳み込みニューラルネットワークは、完全な畳み込みニューラルネットワークでありうる。
本発明のさらなる一態様は、組織学的画像(またはそのセット)内の腫瘍を同定するためのコンピュータプログラム製品、上記の方法を実行するための機械可読命令を有するコンピュータプログラム製品に関する。
本発明のさらに別の態様は、組織学的画像(またはそのセット)内の腫瘍を同定するためのコンピュータ装置に関し、当該装置は、
例えば、保存されたレコードから、組織学的画像(またはそのセット)を受信するように動作可能な入力部であって、組織学的画像が、ピクセルの2次元アレイを含む、入力部と、
処理モジュールであって、
畳み込みニューラルネットワークを適用して、組織学的画像の2次元アレイに対するマッピングを有するピクセルの2次元アレイを有する出力画像を生成することであって、出力画像が、複数の組織クラスのうちの1つを各ピクセルに割り当てることによって生成され、複数の組織クラスが、非腫瘍組織を表す少なくとも1つのクラスおよび腫瘍組織を表す少なくとも1つのクラスを含む、ことと、
個々の腫瘍によって占有された領域がマーキングされている出力画像から、セグメンテーションマスクを生成することと、
各腫瘍の要約統計量を計算することと、
各腫瘍の要約統計量にフィルタを適用して、フィルタに従って腫瘍を選択および選択解除することにより、セグメンテーションマスクを編集することと、を行うように構成された、処理モジュールと、を含む。
コンピュータ装置は、出力画像パッチを組織学的画像(またはそのセット)の確率マップにアセンブルするように構成された後処理モジュールをさらに含みうる。さらに、コンピュータ装置は、確率マップを、例えば、データリポジトリ内のレコードに保存するように動作可能な出力部をさらに含んでもよく、その結果、確率マップは、組織学的画像(またはそのセット)にリンクされる。装置はさらに、組織学的画像(またはそのセット)および確率マップをディスプレイに送信するように動作可能なディスプレイおよびディスプレイ出力部を含んでもよく、その結果、組織学的画像は、確率マップとともに、例えば、その上にオーバレイして、または確率マップとともに表示される。
本発明の別の態様は、1つ以上の他の要素と組み合わせて上記で指定されたコンピュータ装置を含むシステムに関する。例えば、システムは、組織学的画像またはそのセットを取得し、それらをデータリポジトリ内のレコードに保存するように動作可能な顕微鏡など、組織学的画像を取得するように動作可能な画像取得装置を含みうる。かかる顕微鏡は、デジタルスキャナに組み込むことができる。システムはまた、組織学的画像を含む患者データのレコードを保存するように構成されたデータリポジトリと、コンピュータ装置とデータリポジトリとの間の患者データレコードまたはその一部の転送を可能にするネットワーク接続と、を含みうる。
要約すると、組織学的画像において腫瘍を同定および可視化するための方法、コンピュータ装置、コンピュータプログラム製品およびシステムが提供されうる。畳み込みニューラルネットワークCNNを使用して、画像内のピクセルは、非腫瘍組織に関連すると判別されたのか、または腫瘍組織を表す1つ以上のクラスに関連すると判別されたのかに従って分類される。CNNの結果に基づいてセグメンテーションが実行され、マスクが生成され、当該マスクにより、個々の腫瘍によって占有された領域がマーキングされる。次に、各腫瘍の要約統計量が計算され、フィルタに提供される。フィルタにより、重要でないとみなされた腫瘍を除外することによってセグメンテーションマスクが編集される。任意に、フィルタを通過する腫瘍は、要約統計量に従って、例えば臨床的関連性の順序で、または病理医にとって常識的なレビューの順序でランク付けすることができる。次に、ビジュアライゼーションアプリケーションにより、セグメンテーションマスク、要約統計量、および/またはランク付けを考慮して組織学的画像を表示することができる。
組織学的画像は、1つ以上のさらなる2次元ピクセルアレイを含んでもよく、したがって、幾つかの場合では、例えば、複数の染色の各々に1つ、または有限の深度の透明または半透明のサンプルを通して顕微鏡の焦点面をステップすることによって得られたサンプル(いわゆるzスタック)の異なる深度の各々に1つなど、ピクセルの複数の2次元アレイを含みうる。CNNによって生成された出力画像はまた、ピクセルの1つ以上の2次元アレイを含む。ここで、(入力)組織学的画像と出力画像との間に定義されたマッピングが存在しており、これは1対1のマッピング、多対1のマッピング、または1対多のマッピングでありうる。少なくとも幾つかの実施形態では、組織学的画像は、顕微鏡、特に光学顕微鏡によってセクション化された組織サンプルから撮影された2次元画像のデジタル表現であることが理解されよう。これは、従来の光学顕微鏡、共焦点顕微鏡、または未染色もしくは染色された組織サンプルの組織学的画像を取得するのに適した他の種類の顕微鏡でありうる。組織学的画像のセットの場合、これらは、組織の領域の隣接するセクション(すなわちスライス)から撮影された一連の顕微鏡画像であってよく、各セクションは、異なって染色されていてよい。単一の組織学的画像はまた、単純なオーバレイまたはワーピング変換によって支援されたオーバレイを使用して、組織の領域の異なって染色された隣接するセクションから取得された複数の組織学的画像をマージすることによって計算された合成物でありうる。
一実施形態では、CNNは、外科的切除後の腫瘍辺縁の状態を自動的に求めるためのコンピュータ自動化方法に適用されうる。
本開示の一態様によれば、腫瘍塊の辺縁状態を求めるために組織学的画像を画像処理する方法が提供され、当該方法は、
切除によって抽出され、その表面を強調表示するために染料で着色された腫瘍塊から得られた組織サンプル切片の組織学的画像を受信することであって、組織学的画像が、ピクセルの2次元アレイを含む、ことと、
組織クラスを区別するように訓練された少なくとも1つのニューラルネットワークを使用してCNNプロセスを適用して、組織学的画像の2次元アレイに対するマッピングを有するピクセルの2次元アレイを有する出力画像を生成することであって、出力画像が、組織クラスを各ピクセルに割り当てて、染料および非染料ならびに腫瘍および非腫瘍を区別することによって生成される、ことと、
腫瘍染料ピクセルの有無に応じて、辺縁状態を陽性または陰性としてそれぞれ計算することと、
辺縁状態を出力することと、
を含む。
陽性または陰性の辺縁状態は、腫瘍染料ピクセルの数が閾値数を上回っているか、または下回っているかに少なくとも部分的に基づいて計算されうる。陽性状態および陰性状態は細分化できる。辺縁状態の計算により、陰性状態をクロースサブ状態とクリアサブ状態との間で細分化することができ、クロースとクリアとは、腫瘍非染料ピクセルと染料非腫瘍ピクセルとの間の距離が閾値を下回っているかまたは上回っていることによって区別される。辺縁状態の計算により、当技術分野でまとめて辺縁病変と称される複数のサブ状態間で陽性状態を細分化することができる。辺縁病変は、巣性、最小/中程度、および広範囲であると言われている。辺縁病変は、腫瘍染料ピクセルの有病率に基づいて求めることができ、有病率は、腫瘍染料ピクセル間の最大距離の計算に少なくとも部分的に基づいて定量化することができる。
発明者らは、深層学習を使用して腫瘍の辺縁状態を検出し、辺縁距離および辺縁病変をも報告しうる画像分析アルゴリズムを使用することを提案する。アルゴリズムは、病理医がレビューできる仮想スライド画像のバッチに対して教師なしで実行できる。このコンピュータ自動化方法を通じてこの辺縁情報を病理医に提供することにより、病理医のワークフローの速度を向上させ、より高品質でより一貫性のある腫瘍辺縁レポートを作成できる。
実施形態の1つのグループでは、単一の畳み込みニューラルネットワークを使用して、完全な分類が実行される。つまり、畳み込みニューラルネットワークプロセスは、腫瘍染料、腫瘍非染料、非腫瘍染料、および非腫瘍非染料の組織クラスを区別するように訓練された1つのニューラルネットワークを適用することを含む。
実施形態の別のグループでは、別個の畳み込みニューラルネットワークが適用され、一方は腫瘍型を分類し、他方は染料と非染料とを分類し、これらの2つのニューラルネットワークの結果を論理的に組み合わせてクラスの完全なセットが生成される。すなわち、畳み込みニューラルネットワークプロセスは、腫瘍および非腫瘍の組織クラスを区別し、第1の出力画像を出力するように訓練された第1のニューラルネットワークと、染料および非染料の組織クラスを区別して第2の出力画像を出力するように訓練された第2のニューラルネットワークとを別々に適用することを含み、当該出力画像は、第1の出力画像および第2の出力画像を組み合わせて、組織クラスを腫瘍染料、腫瘍非染料、非腫瘍染料、および非腫瘍非染料の組織クラスに割り当てることによって生成される。
分類スキームには、腫瘍組織型が1つしかない場合もあれば、複数ある場合もある。一実施形態では、浸潤性腫瘍用および非浸潤性腫瘍用に別個の腫瘍組織型があり、非腫瘍組織用に1つのクラスがあり、そのため、腫瘍組織型と並行して、染料と非染料との間の違いを考慮した6つのクラスのセットが提供される。
処理のための方法によって受信された組織学的画像は、仮想スライドライブラリまたは他のデータベースに保存されたレコードからロードされうるか、またはデジタルスキャナによる組織学的画像の取得の直後に取得されうる。その場合、この方法は、デジタルスキャナ自体によって実行可能であり、あるいは例えば、デジタルスキャナに接続されているか、または別のデータソースから接続されたコンピュータ装置によって実行されうる。
この方法は、割り当てステップの後に、出力画像を組織学的画像(またはそのセット)の確率マップにアセンブルするステップをさらに含みうる。この方法は、アセンブルするステップの後、組織学的画像またはそのセットにリンクされた確率マップを、例えば、組織学的画像データのメタデータを構成する確率マップとの共通レコードに保存するステップをさらに含みうる。
本実装形態において、連続する各畳み込み段階では、次元が削減されるにつれて深度が増大するため、畳み込み層の深度は増大し、次元は削減される。また、連続する転置畳み込み段階では、次元が増大するにつれて深度が減少するため、逆畳み込み層の深度は減少し、次元は増大する。最終畳み込み層は、最大の深度および最小の次元を有する。畳み込み段階および逆畳み込み段階でそれぞれ深度が増大および減少するアプローチの代わりに、入力層と出力層とを除く全ての層が同じ深度を有するニューラルネットワークを設計することもできる。
この方法はさらに、組織学的画像またはそのセットを確率マップとともにディスプレイ上に表示すること、例えば、その上にオーバレイして、または互いに並べて表示することを含みうる。
特定の実施形態では、畳み込みニューラルネットワークは、1つ以上のスキップ接続を有する。各スキップ接続は、畳み込み層のうちの最終畳み込み層よりも大きな次元の少なくとも1つの畳み込み層から中間結果を取得し、当該結果を、なし、1つ、または2つ以上でありうる必要な数の転置畳み込みに供して、入力画像パッチにサイズが一致する少なくとも1つのさらなる復元された層を取得する。次に、これらは、各ピクセルに組織クラスを割り当てる該ステップの前に、上記の復元された層と組み合わされる。さらなる処理ステップでは、確率を再計算するために、復元された層をさらなる復元された層の各々と組み合わせ、これによって、スキップ接続から得られた結果が考慮される。
特定の実施形態では、ソフトマックス演算を使用して確率を生成する。
組織学的画像から抽出された画像パッチは、画像の全領域をカバーすることができる。パッチは、重複しない画像タイル、またはマージン部分で重複して確率マップのつなぎ合わせを支援する画像タイルでありうる。CNNは固定サイズのピクセルアレイのみを受け入れるように設計されているため、各画像パッチは、CNNに一致する幅および高さの固定数のピクセルを有する必要があるが、このことは、各画像パッチが組織学的画像上の同じ物理的領域に対応する必要があることを意味するものではない。なぜなら、組織学的画像のピクセルは、より広い領域をカバーする低解像度のパッチに結合される可能性があるからである。例えば、隣接するピクセルの各2×2アレイを組み合わせて、1つの「スーパー」ピクセルとし、組織学的画像のネイティブ解像度で抽出されたパッチの物理的領域の4倍のパッチを形成することができる。
本方法は、CNNが訓練された後、予測のために実行することができる。この訓練の目的は、層間接続に適切な重み値を割り当てることである。訓練の場合、使用されるグラウンドトゥルースデータは、訓練データの組織学的画像の各ピクセルを組織クラスの1つに割り当てる。グラウンドトゥルースデータは、専門の臨床医が十分な数の画像に注釈を付けることにより得られる。訓練は、CNNを繰り返し適用することによって実行され、各反復では、グラウンドトゥルースデータと出力画像パッチとの比較に基づいて重み値が調整される。本実装形態では、重みは、勾配降下法による訓練中に調整される。
一方が腫瘍組織型を分類し、他方が染料および非染料を分類する別々の畳み込みニューラルネットワークを使用する実施形態では、一方は腫瘍のセグメンテーションのみに基づき、他方は染料境界のセグメンテーションのみに基づく2つの別々の訓練データセットを使用できるため、訓練データの提供はより容易でありうる。したがって、染料訓練データは、腫瘍型、または腫瘍組織型の下位分類の種類に固有ではない可能性がある。さらに、腫瘍訓練データには、具体的に染料を使用した切除スライドからの画像ではなく、任意の適切な病理画像を使用することができる。一方、同じ畳み込みニューラルネットワークを4方向分類、6方向分類、またはそれ以上の偶数の分類で使用する場合、訓練データセットは、適切に専門的に分析され、グラウンドトゥルースデータで補強された、着色された切除組織塊の画像である必要がある。
腫瘍組織クラスの設定には種々の選択肢があるが、全ての実施形態でないにしても、大部分の実施形態では、非腫瘍組織と腫瘍組織との間でクラスが区別されることが共通している。非腫瘍組織クラスは、1つ、2つ、またはそれ以上のクラスを含みうる。腫瘍組織クラスはまた、1つ、2つ、またはそれ以上のクラスを含みうる。例えば、本実装形態では、3つの組織クラスがある。1つは非腫瘍組織クラス、2つは腫瘍組織クラスである。2つの腫瘍組織クラスは、浸潤性腫瘍クラスおよび非浸潤性腫瘍クラスである。
幾つかの実施形態では、CNNは、一度に1つの組織学的画像に適用される。他の実施形態では、CNNは、組織の領域の異なって染色された隣接するセクションから取得された組織学的画像のセットを組み合わせることによって形成された合成組織学的画像に適用されうる。さらに別の実施形態では、CNNは、組織の領域の異なって染色された隣接するセクションから取得された画像のセットの画像の各々と並行して適用されうる。
結果は、臨床医のディスプレイに表示されうる。すなわち、組織学的画像は、関連する確率マップとともに表示することができ、例えば、その上にオーバレイして、または互いに並べて表示することができる。腫瘍スコアは、例えば、腫瘍にテキストラベルを付けたり、腫瘍を指し示したり、画像の横に表示したりするなど、便利な方法で表示することもできる。
畳み込みニューラルネットワークは、完全な畳み込みニューラルネットワークでありうる。
本発明のさらなる一態様は、組織学的画像またはそのセット内の腫瘍を同定するためのコンピュータプログラム製品、上記の方法を実行するための機械可読命令を有するコンピュータプログラム製品に関する。
本発明のさらに別の態様は、組織学的画像を画像処理して腫瘍塊の辺縁状態を求めるためのコンピュータ装置に関し、当該装置は、
切除によって抽出され、その表面を強調表示するために染料で着色された腫瘍塊から得られた組織サンプル切片の組織学的画像を受信するように動作可能な入力部であって、組織学的画像が、ピクセルの2次元アレイを含む、入力部と、
組織クラスを区別し、組織学的画像の2次元アレイに対するマッピングを有するピクセルの2次元アレイを有する出力画像を生成するように訓練された少なくとも1つのニューラルネットワークを含む畳み込みニューラルネットワーク処理モジュールであって、出力画像が、組織クラスを各ピクセルに割り当てて、染料および非染料ならびに腫瘍および非腫瘍を区別することによって生成される、畳み込みニューラルネットワーク処理モジュールと、
出力画像内の腫瘍染料ピクセルの有無に応じて、辺縁状態を陽性または陰性としてそれぞれ計算するように動作可能な辺縁評価モジュールと、
辺縁状態を出力するように動作可能な出力部と、
を含む。
出力をデータリポジトリに接続して、辺縁状態を保存し、任意に出力画像をレコードに保存することもできる。保存された出力画像は、各ピクセルの組織クラスおよび/または各ピクセルの確率、すなわち確率マップを含みうる。出力画像を保存することは、ビジュアライゼーションを構築するのに役立ちうる。これは、例えば、組織学的画像は、組織クラスマップまたは確率マップをオーバレイして表示できるためである。レコードに保存されているかどうかにかかわらず、出力部は、組織学的画像が出力画像とともに表示されるように、組織学的画像および出力画像をディスプレイに送信するように動作可能でありうる。
本発明の別の態様は、以下の要素、すなわち、組織学的画像またはそのセットを取得するように動作可能な顕微鏡などの画像取得装置、および組織学的画像(またはそのセット)を含む患者データのレコードを保存するように構成されたデータリポジトリのうちの1つ以上と組み合わせた上記のようなコンピュータ装置を含む、クリニカルネットワークなどのシステムに関する。コンピュータ装置とデータリポジトリとの間の患者データレコードまたはその一部の転送など、データ転送を可能にするこれらの要素のうちのいずれかの要素間に適切なネットワーク接続が提供されることが理解されよう。
少なくとも幾つかの実施形態では、組織学的画像は、顕微鏡、特に光学顕微鏡によってセクション化された組織サンプルから撮影された2次元画像のデジタル表現であることが理解されよう。これは、従来の光学顕微鏡、共焦点顕微鏡、または未染色もしくは染色された組織サンプルの組織学的画像を取得するのに適した他の種類の顕微鏡でありうる。組織学的画像のセットの場合、これらは、組織の領域の隣接するセクション(すなわちスライス)から撮影された一連の顕微鏡画像であってよく、各セクションは、異なって染色されていてよい。
本発明の他の特徴および利点は、以下の詳細な説明および添付の図面を検討した後、当業者にはより容易に明らかになるであろう。
本発明の構造および動作は、以下の詳細な説明および添付図面を検討することで理解され、ここで、同様の参照番号は、同様の部品を指す。
本発明の一実施形態で使用されるニューラルネットワークアーキテクチャの概略図である。 図1Aのニューラルネットワークアーキテクチャ内で、大域特徴および局所特徴マップがどのように組み合わされて、入力画像パッチ内の各ピクセルの個々のクラスを予測する特徴マップを生成するかを示す図である。 生の組織パッチ画像を示す図であり、当該画像は、動作中、カラー画像である。 図2AのCNN予測を示す図であり、当該図面は、動作中、カラー画像である。CNN予測は、予測された腫瘍領域を示しており、これは、カラー画像では暗赤色で示されている。 生の組織学的画像を示す図であり、当該画像は、動作中、カラー画像である。マーキングされた生の組織学的画像は、病理医の手による注釈(赤色)とCNN予測(ピンク色および黄色)とを示している。 CNN予測を示す図であり、当該図面は、動作中、カラー画像である。CNN予測画像は、非腫瘍領域(緑色)、浸潤性腫瘍領域(赤色、図2Cのピンク色に対応)、および非浸潤性腫瘍(青色、図2Cの黄色に対応)を示している。 入力RGB画像パッチの例を示す図であり、当該図面は、動作中、カラー画像である。画像パッチは、病理医による浸潤性腫瘍の手動での輪郭付け(赤色)を示し、さらにニューラルネットワークの予測(ピンク色および黄色)のオーバレイを示している。 最終出力腫瘍確率ヒートマップを示す図であり、当該図面は、動作中、カラー画像である。ヒートマップは、ニューラルネットワークの予測のオーバレイを示している(それぞれ赤褐色および青色)。 CNNの訓練に含まれるステップを示すフロー図である。 CNNを使用した予測に含まれるステップを示すフロー図である。 訓練済みCNNを使用して辺縁を評価することに含まれるステップを示すフロー図である。 本開示の一実施形態による方法のフロー図である。 図1Aおよび図1Bのニューラルネットワークアーキテクチャの実装に含まれる計算を実行するために使用されうるTPUのブロック図である。 本発明の実施形態と併せて使用することができる例示的なコンピュータネットワークを示す図である。 例えば、図8のTPUのホストコンピュータとして使用することができるコンピューティング装置のブロック図である。 本明細書に記載の種々の実施形態に関連して使用されうる例示的なプロセッサ対応デバイス550を示すブロック図である。 単一の線形アレイを有する例示的なラインスキャンカメラを示すブロック図である。 3つの線形アレイを有する例示的なラインスキャンカメラを示すブロック図である。 複数の線形アレイを有する例示的なラインスキャンカメラを示すブロック図である。
以下の詳細な説明では、限定ではなく説明の目的で、本開示のより良い理解を提供するために特定の詳細が記載されている。本開示は、これらの特定の詳細から逸脱する他の実施形態で実施されうることが当業者には明らかであろう。
癌は、北米女性の第2に多い死因である。全種類の女性の癌の中で、乳癌は最も一般的なものであり、癌による死亡原因の第2位となっている。そのため、女性は、人生のどこかで乳癌に罹患する割合が多く、乳癌治療の精度は、女性の寿命および生活の質に大きな影響を及ぼしている。
乳癌は、特定の遺伝子の発現に基づいて、種々の分子サブタイプに分類できる。概して使用される分類スキームは、
1.ルミナルA型:ER+、PR+、HER2-
2.ルミナルB型:ER+、PR-、HER2+
3.トリプルネガティブ乳癌(TNBC)型:ER-、PR-、HER2-
4.HER2濃縮型:HER2+、ER-、PR-
5.ノーマル-like型
の通りである。
ERはエストロゲン受容体の略である。PRはプロゲステロン受容体の略である。HER2はヒト上皮成長因子受容体2の略である。
浸潤性および非浸潤性乳癌細胞の核を自動的に検出して輪郭を描画するコンピュータ自動腫瘍発見法について説明する。この方法は、WSIなどの単一の入力画像、またはWSIのセットなどの入力画像のセットに適用される。各入力画像は、WSIなどのデジタル化された組織学的画像である。入力画像のセットの場合、これらは隣接する組織切片の異なって染色された画像でありうる。染色という用語は、広く使用されており、バイオマーカーによる染色だけでなく、従来のコントラスト増強染色による染色も含む。
コンピュータで自動化された腫瘍の輪郭付けは手動の輪郭付けよりもはるかに高速であるため、画像から選択した抽出タイルに手動で注釈を付すだけでなく、画像全体を処理することができる。したがって、提案している自動腫瘍輪郭付けにより、病理医は画像内の全ての腫瘍細胞の陽性(または陰性)割合を計算できるようになり、これにより、より正確で再現性のある結果が得られる。
腫瘍の発見、輪郭付け、および分類のために提案されたコンピュータ自動化方法は、畳み込みニューラルネットワーク(CNN)を使用して、WSI上の各核ピクセルを発見し、かかる各ピクセルを、本実装形態の乳房腫瘍クラスにおいて、非腫瘍クラスの1つと複数の腫瘍クラスの1つとに分類する。
発明者らの実装形態におけるニューラルネットワークは、設計の点で、<http://www.robots.ox.ac.uk/vgg/research/very_deep/>で入手可能なVGG-16アーキテクチャと類似しており、Simonyan and Zisserman(2014)により説明されている。その内容は全て、参照により本明細書に組み込まれている。
入力画像は、病理画像であり、本明細書の他の箇所でより詳細に説明されるように、複数の従来の染色のうちのいずれかで染色された画像である。CNNの場合、画像パッチは特定のピクセル寸法、例えば128×128、256×256、512×512、または1024×1024ピクセルで抽出される。画像パッチは任意のサイズであってよく、正方形である必要はないが、パッチの行および列のピクセル数は2nに一致し、nは正の整数であることが理解されよう。これは、かかる数値が、一般に、適切な単一CPU(中央処理装置)、GPU(グラフィックス処理装置)、TPU(テンソル処理装置)、またはそのアレイによる直接のデジタル処理により適しているためである。
「パッチ」は、WSIから取得された画像部分を指すために使用される専門用語であり、通常、正方形または長方形の形状であることに留意されたい。この点で、WSIには10億以上のピクセル(ギガピクセル画像)が含まれうるため、画像処理は通常、CNNで処理するための管理可能なサイズ(例えば、約500×500ピクセル)のパッチに適用される。したがって、WSIは、パッチに分割し、CNNでパッチを分析し、出力(画像)パッチをWSIと同じサイズの確率マップに再アセンブルすることに基づいて処理される。次に、確率マップを、例えば半透明に、WSIまたはその一部にオーバレイすることができ、その結果、病理画像および確率マップの両方を一緒に観察することができる。その意味で、確率マップは病理画像のオーバレイ画像として使用される。CNNによって分析されるパッチは、全て同じ倍率のものでありうるか、または異なる倍率、例えば5倍、20倍、50倍などの混合物を有しうるため、サンプル組織の異なるサイズの物理的領域に対応する。異なる倍率によって、これらは、WSIが取得された物理的倍率に対応することができ、あるいはより高倍率(すなわち、より高解像度)の物理的画像をデジタル的に縮小することで得られる有効倍率に対応することができる。
図1Aは、発明者らのニューラルネットワークアーキテクチャの概略図である。層C1,C2,...,C10は畳み込み層である。層D1、D2、D3、D4、D5およびD6は、転置畳み込み(つまり逆畳み込み)層である。特定の層を相互接続する線は、畳み込みC層と逆畳み込みD層との間のスキップ接続を示している。スキップ接続により、より次元が大きく深度が浅い層(「大きい」および「浅い」は、より低いインデックスの畳み込み層を意味する)の局所特徴を、最後の(つまり、最小、最深の)畳み込み層からの大域特徴と組み合わせることができる。これらのスキップ接続により、より正確なアウトラインが提供される。最大プール層は、その各々がパッチの幅および高さを2分の1に減らすために使用されており、層C2、C4およびC7の後に存在しているが、概略図には直接に示されておらず、パッチのサイズを結果的に小さくすることで暗示的に示されている。ニューラルネットワークの幾つかの実現形態では、最大プール層が1×1畳み込みに置き換えられ、完全な畳み込みネットワークになる。
ニューラルネットワークの畳み込み部分は、順番に、入力層(RGB入力画像パッチ)と、2つの畳み込み層C1、C2と、第1の最大プール層(図示せず)と、2つの畳み込み層C3、C4と、第2の最大プール層(図示せず)と、3つの畳み込み層C5、C6、C7と、第3の最大プール層(図示せず)とを有する。第2および第3の最大プール層からの出力は、それぞれ、層C5およびC8への通常の接続に加えて、スキップ接続を使用して逆畳み込み層に直接に接続される。
次に、最後の畳み込み層C10、第2の最大プール層(つまり、層C4の後の層)からの出力、および第3の最大プール層(つまり、層C7の後の層)からの出力は、それぞれ「逆畳み込み層」の個別のシーケンスに接続され、これらは、入力(画像)パッチと同じサイズにアップスケールされる。つまり、畳み込み特徴マップを、入力画像パッチと同じ幅および高さで、検出される組織クラスの数に等しい、すなわち非腫瘍型および1つ以上の腫瘍型に等しいチャネル数(つまり、特徴マップの数)を有する特徴マップに変換する。第2の最大プール層では、逆畳み込みの1つのステージのみが必要なため、層D6への直接のリンクが表示される。第3の最大プール層の場合、層D5に到達するには、中間逆畳み込み層D4を介して2段階の逆畳み込みが必要である。最も深い畳み込み層C10の場合、D1およびD2を経由して層D3に至る3段階の逆畳み込みが必要である。結果は、入力パッチと同じサイズの3つのアレイD3、D5、D6である。
図1Aに示されているものを単純化したバージョンでは、おそらく性能は劣るものの、スキップ接続を省略することができ、この場合、層D4、D5およびD6は存在せず、出力パッチは層D3からのみ計算される。
図1Bは、図1Aのニューラルネットワークアーキテクチャにおける最終ステップがどのように実行されるかをより詳細に示している。すなわち、大域特徴マップ層D3および局所特徴マップ層D5、D6が組み合わされて、入力画像パッチの各ピクセルの個々のクラスを予測する特徴マップを生成している。具体的には、図1Bは、最後の3つの転置畳み込み層D3、D5、D6が腫瘍クラスの出力パッチにどのように処理されるかを示している。
ここで、上記のアプローチが、デジタルパソロジーで現在使用されている既知のCNNとどのように異なるかについて説明する。この既知のCNNでは、複数の利用可能なクラスから選択された1つのクラスを各画像パッチに割り当てている。かかるタイプのCNNの例として、Wang et al.(2016)、Liu et al.(2017)、Cruz-Roa et al.(2017)、Vandenberghe et al.(2017)の論文が挙げられる。ただし、ここで説明されているのは、所与の画像パッチ内で、複数の使用可能なクラスから選択された1つのクラスが全てのピクセルに割り当てられるということである。したがって、画像パッチごとに1つのクラスラベルを生成する代わりに、発明者らのニューラルネットワークは、所与のパッチの個々のピクセルごとにクラスラベルを出力している。発明者らの出力パッチは、入力パッチと1対1のピクセル間で対応しているため、出力パッチの各ピクセルには、複数の利用可能なクラスのうちの1つ(非腫瘍、腫瘍1、腫瘍2、腫瘍3など)が割り当てられている。
かかる既知のCNNでは、各パッチに単一のクラスを割り当てるために、一連の畳み込み層が使用され、その後に1つ以上の全結合層が続き、検出するクラスと同じ数の値を有する出力ベクトルが続いている。予測されるクラスは、出力ベクトルの最大値の位置によって求められる。
訓練済みCNNは、入力としてデジタルスライド画像からピクセルを取得し、各ピクセルの確率のベクトルを返す(Goodfellow,Bengio,and Courville(2016))。ベクトル長はNであり、Nは、CNNが検出するように訓練されたクラスの数である。例えば、CNNが、浸潤性腫瘍、非浸潤性腫瘍、非腫瘍の3つのクラスを区別するように訓練されている場合、ベクトル長vは3になる。ベクトルの各座標は、ピクセルが特定のクラスに属する確率を示している。したがって、v[0]は、ピクセルが浸潤性腫瘍クラスに属する確率、v[1]はピクセルが非浸潤性クラスに属する確率、v[2]はピクセルが非腫瘍クラスに属する確率を示すことができる。各ピクセルのクラスは、確率ベクトルから求められる。ピクセルをクラスに割り当てる単純な方法は、ピクセルを最も確率の高いクラスに割り当てることである。
個々のピクセルのクラスを予測するために、発明者らのCNNでは、畳み込み層に続く異なるアーキテクチャが使用される。一連の全結合層の代わりに、一連の転置畳み込み層により畳み込み層が追跡される。全結合層は、このアーキテクチャから削除される。各転置層は、特徴マップの幅および高さを2倍にすると同時に、チャネル数を半分にする。このようにして、特徴マップは入力パッチのサイズにアップスケールされる。
さらに、予測を改善するために、Long et al.(2015)に記載されているスキップ接続が使用され、その内容は全て、参照により本明細書に組み込まれている。
スキップ接続により、より浅い特徴を使用して、最終畳み込み層C10からのアップスケーリングによって行われる粗予測が改善される。図1Aの層D5およびD6に含まれるスキップ接続からの局所特徴は、図1Aの層D3に含まれる大域特徴を最終畳み込み層からアップスケーリングすることによって生成された特徴と連結される。次に、大域特徴層および局所特徴層D3、D5およびD6は、図1Bに示されるように、組み合わされた層に連結される。
図1Bの連結された層から(またはスキップ接続が使用されない場合は最後の逆畳み込み層D3から直接に)、チャネルの数は、組み合わされた層の1×1畳み込みによってクラスの数と一致するように減少される。次に、この分類層でのソフトマックス演算により、結合された層の値が確率に変換される。出力パッチ層のサイズはN×N×Kであり、Nは入力パッチのピクセル単位の幅および高さ、Kは検出されているクラスの数である。したがって、画像パッチ内の任意のピクセルPに対して、サイズKの出力ベクトルVが存在する。次に、対応するベクトルVの最大値の位置によって、一意のクラスを各ピクセルPに割り当てることができる。
したがって、CNNは、各ピクセルを、非癌性または複数の異なる癌(腫瘍)タイプのうちの1つ以上に属するものとしてラベル付けする。特に関心のある癌は乳癌であるが、本方法は、膀胱癌、結腸癌、直腸癌、腎臓癌、血液(白血病)、子宮内膜癌、肺癌、肝臓癌、皮膚癌、膵臓癌、前立腺癌、脳癌、脊椎癌、甲状腺癌などの他の癌の組織学的画像にも適用可能である。
発明者らの特定のニューラルネットワークの実装形態は、特定の固定ピクセル寸法を有する入力画像で動作するように構成されている。したがって、訓練および予測の両方の前処理ステップとして、N×N×nピクセルなどの目的のピクセル寸法を有するパッチがWSIから抽出され、ここで、WSIが従来の可視光顕微鏡によって取得されたカラー画像であるとき、各物理的位置に3つの原色(通常はRGB)に関連付けられた3つのピクセルがある場合はn=3となる(以下でさらに説明するように、2つ以上のカラーWSIが組み合わされている場合、「n」は合成WSIの数の3倍になりうる)。さらに、単一のモノクロWSIの場合、「n」の値は1になる。訓練を高速化するために、入力パッチもこの段階でセンタリングされ、正規化される。
発明者らの好ましいアプローチは、WSI全体、または少なくとも組織を含むWSIの領域全体を処理することである。したがって、この場合のパッチは、少なくともWSIの組織領域全体をカバーするタイルである。タイルは、オーバラップすることなく隣接している場合があり、あるいは例えば1、2、3、4、5、6、7、8、9、または10ピクセル幅の重複するエッジマージン領域を有しているため、CNNの出力パッチを、不一致を考慮してつなぎ合わせることができる。しかしながら、発明者らのアプローチは、必要に応じて、従来技術のように、または病理医によって実行されるように、同じまたは異なる倍率であるWSI上のパッチのランダムサンプルに適用することもできる。
発明者らのニューラルネットワークは、設計の点で、Simonyan and Zisserman(2014)のVGG-16アーキテクチャに類似している。全ての畳み込みフィルタにおいて、非常に小さい3×3カーネルが使用される。最大プーリングは、2×2の小さなウィンドウおよび2のストライドで実行される。畳み込み層の後に一連の全結合層を有するVGG-16アーキテクチャとは対照的に、発明者らは、畳み込み層の後に一連の「逆畳み込み」(より正確には転置畳み込み)のシーケンスを行い、セグメンテーションマスクを生成している。セマンティックセグメンテーションのためのこのタイプのアップサンプリングは、過去に、Long et al.(2015)によって自然画像処理に使用されており、その内容は全て、参照により本明細書に組み込まれている。
各逆畳み込み層は、入力特徴マップを幅および高さの次元で2倍に拡大する。これにより、最大プール層の縮小効果が打ち消され、入力画像と同じサイズのクラス特徴マップが作成される。各畳み込み層および逆畳み込み層からの出力は、非線形アクティベーション層によって変換される。ここでは、非線形アクティベーション層は、整流関数ReLU(x)=max(0,x)を使用している。必要に応じて、ReLU、リーキーReLU、eLUなどの種々のアクティベーション関数を使用できる。
提案している方法は、任意の数の組織クラスに変更を加えることなく適用できる。その制約は、ニューラルネットワークで複製することが望ましい方法で分類された適切な訓練データの可用性にすぎない。さらなる乳房病変の例として、浸潤性小葉癌または浸潤性乳管癌が挙げられる。すなわち、前述の例の単一の浸潤性腫瘍クラスは、複数の浸潤性腫瘍クラスで置き換えることができる。ニューラルネットワークの精度は、ほとんどの場合、各クラスで使用可能な画像の数、クラスの類似性、およびメモリ制限に陥る前にニューラルネットワークをどれだけ深くすることができるかによって決まる。概して、クラスごとの画像の数が多い、ネットワークが深い、クラスが異なるなどの条件が揃うと、ネットワークの精度が高くなる。
ソフトマックス回帰層(つまり、多項ロジスティック回帰層)がチャネルパッチの各々に適用され、特徴マップの値が確率に変換される。
最終特徴マップのチャネルCの位置(x,y)の値には、確率P(x,y)が含まれ、確率P(x,y)では、入力画像パッチの位置(x,y)のピクセルがチャネルCによって検出された腫瘍型に属している。
畳み込み層および逆畳み込み層の数は、必要に応じて増減され、ニューラルネットワークを実行するハードウェアのメモリ制限を受けうることが理解されよう。
ミニバッチ勾配降下法を使用してニューラルネットワークが訓練される。学習率は、指数関数的減衰を用いて、初期率0.1から減少する。Srivastava et al.(2014)[2017]によって説明されている「ドロップアウト」手順を使用することにより、ニューラルネットワークの過剰適合が防止されており、その内容は全て、参照により本明細書に組み込まれている。ネットワークの訓練は、利用可能な複数の深層学習フレームワークのいずれかを使用して、GPU、CPU、またはFPGAで実行できる。本実装形態では、Google Tensorflowを使用しているが、同じニューラルネットワークをMicrosoft CNTKなどの別の深層学習フレームワークに実装することもできる。
ニューラルネットワークは、サイズN×N×Kの確率マップを出力する。ここで、Nは入力パッチのピクセル単位の幅および高さ、Kは検出されているクラスの数である。これらの出力パッチは、サイズW×H×Kの確率マップにステッチバックされ、WおよびHは、パッチに分割される前の元のWSIの幅および高さである。
次に、ラベル画像の各位置(x,y)で最大確率でクラスインデックスを記録することにより、確率マップをW×Hラベル画像に折りたたむことができる。
本実装形態では、発明者らのニューラルネットワークにより、全てのピクセルが3つのクラス(非腫瘍、浸潤性腫瘍、非浸潤性腫瘍)のいずれかに割り当てられる。
複数の腫瘍クラスを使用する場合、出力画像を後処理して、非腫瘍および腫瘍のより単純なバイナリ分類にすることができる。つまり、複数の腫瘍クラスを組み合わせることができる。バイナリ分類は、ベースデータから画像を作成する際の選択肢として使用できるが、マルチクラスの腫瘍分類は保存データに保持される。
本発明の特定の実装形態に関する上記の説明は、CNNを使用する特定のアプローチに集中しているが、発明者らのアプローチは、多種多様な異なるタイプの畳み込みニューラルネットワークで実装できることが理解されよう。概して、畳み込みを使用してますます複雑となる特徴を検出し、その後転置畳み込み(「逆畳み込み」)を使用して特徴マップを入力画像の幅および高さにアップスケールするニューラルネットワークが適している。
実施例1
図2A~図2Dは、動作中のカラー画像であり、Camelyonコンペティション用のGoogleのCNNソリューション(Liu et al.(2017))によって生成されたパッチレベルの予測と、本明細書で説明されているCNNによって生成された腫瘍クラスのピクセルレベルの予測との違いを示している。
図2Aおよび図2Bは、Liu et al.(2017)の図7からコピーされたものであり、一方、図2Cおよび図2Dは、腫瘍組織クラス(染料/非染料を含まない)によって分類するために適用される、本明細書で説明されているCNNによる例の比較可能なタイルである。
図2Aは、H&E染色されたWSIからのパッチであり、右下の象限にある大きな濃い紫色の細胞のクラスタは腫瘍であり、小さな濃い紫色の細胞はリンパ球である。
図2Bは、Liu et al.(2017)のCNNによって生成された腫瘍確率ヒートマップであり、著者は、結合組織およびリンパ球を無視して、腫瘍細胞を正確に同定していると述べている。
図2Cは、本発明を具現化したCNN法が適用された例示的なWSIの生の画像パッチである。生の画像と同様に、図2Cは病理医によって手動で描画された輪郭(赤い実線の周囲線)を示している。さらに、図2Dを参照すると、図2Cはまた、本明細書で説明されているCNN法の結果を示している(ピンク色の周囲線でピンク色に網掛けされた第1の領域は、第1の腫瘍型、すなわち、図2Dに赤色で示されている腫瘍型に対応し、ピンク色の周囲線で黄色に網掛けされた第2の領域は、第2の腫瘍型、すなわち、図2Dで青色に網掛けされた腫瘍型に対応する)。
図2Dは、発明者らのCNNによって生成された腫瘍確率ヒートマップである。発明者らのピクセルレベルの予測のアプローチにより、滑らかな周囲の輪郭を有する領域が得られることがわかる。発明者らのヒートマップの場合、異なる(任意に選択された)色は異なるクラスを示す。つまり、非腫瘍の場合は緑色、第1の腫瘍型の場合は赤色、第2の腫瘍型の場合は青色で示されている。
実施例2
図3A~図3Bは、動作中のカラー画像であり、入力RGB画像パッチ(図3A)および最終出力腫瘍確率ヒートマップ(図3B)の例を示している。
図3Aは、さらに、図3Bに別々に示されているように、浸潤性腫瘍の病理医の手動の輪郭付け(赤色の輪郭)、および発明者らのニューラルネットワークの予測のオーバレイ(影付きのピンク色および黄色の領域)をさらに示している。
図3Bは、発明者らのCNNによって生成された腫瘍確率ヒートマップである。発明者らのヒートマップでは、異なる(任意に選択された)色は異なるクラスを示す。つまり、非腫瘍の場合は緑色、浸潤性腫瘍の場合は赤褐色(図3Aにピンク色で表示)、非浸潤性腫瘍の場合は青色(図3Aに黄色で表示)で示されている。繰り返しになるが、発明者らのピクセルレベルの予測のアプローチにより、滑らかな周囲の輪郭を有する領域が得られることがわかる。さらに、CNN予測が病理医の手動マーキングといかに適合しているかがわかる。さらに、CNNは、病理医によって実行されなかった浸潤性組織および非浸潤性(in situ)組織をさらに区別する。これは、本質的には、組織を必要に応じて臨床的に関連する様々なタイプに分類するようにプログラムおよび訓練されうるマルチチャネルCNN設計の一部である。
取得および画像処理
本方法では、組織サンプルがセクション化されていること、すなわちスライスされており、隣接するセクションが異なる染色で染色されていることが出発点となる。隣接するセクションは、セクションが薄いため非常に類似した組織構造を有しているが、異なる層であるため、同一ではない。
例えば、5つの隣接するセクションがあり、それぞれがER、PR、p53、HER2、H&EおよびKi-67などの異なる染色でありうる。次に、各セクションの顕微鏡画像が取得される。隣接するセクションは非常に類似した組織形状を有するが、染色によって、様々な特徴、例えば、核、細胞質、一般的なコントラスト強調による全ての特徴が強調表示される。次に、異なる画像が位置合わせ、ワーピング、または前処理されて、ある画像の所与の特徴の座標が他の画像の同じ特徴にマッピングされる。マッピングにより、わずかに異なる倍率、顕微鏡でのスライドの位置合わせの違いによる方向の違い、またはスライドへの組織スライスの取り付けなどの要因によって引き起こされる画像間の差異が処理される。
異なって染色された隣接セクションを含むセットの異なるWSI間の座標マッピングにより、WSIを単一の合成WSIにマージすることができ、そこから合成パッチを抽出して、CNNによる処理が可能であることに留意されたい。かかる合成パッチの寸法はN×N×3mであり、「m」は、セットを形成する合成WSIの数である。
次に、標準的な画像処理が実行される。これらの画像処理ステップは、WSIレベルまたは個々の画像パッチのレベルで実行できる。CNNがカラー画像ではなくモノクロ画像で動作するように構成されている場合、画像はカラーからグレースケールに変換されうる。画像は、コントラスト強調フィルタを適用することによって変更できる。次に、画像のセット内の一般的な組織領域を特定するため、または単に組織に関係のない背景を除外するために、幾つかのセグメンテーションが実行されうる。セグメンテーションには、次の画像処理技術、すなわち
1.シード組織領域を特定するための分散ベースの分析
2.適応閾値処理
3.モルフォロジー演算(例えば、ブロブ分析)
4.輪郭特定
5.近接ヒューリスティックルールに基づく輪郭のマージ
6.不変画像モーメントの計算
7.エッジ抽出(例えば、ソーベル(Sobel)エッジ検出)
8.曲率流フィルタリング
9.連続セクション間の強度変動を排除するためのヒストグラムマッチング
10.多重解像度のリジッド/アフィン画像レジストレーション(勾配降下オプティマイザ)
11.非剛体変位/変形
12.スーパーピクセルクラスタリング
のうちのいずれかまたは全てが含まれうる。
上記の種類の画像処理ステップがWSI上またはパッチ抽出後の個々のパッチ上で実行できることも理解されよう。幾つかの場合では、パッチ抽出の前後の両方で、つまりそれぞれCNN前処理およびCNN後処理として、同じタイプの画像処理を実行することが有効でありうる。つまり、一部の画像処理はパッチ抽出の前にWSIで実行されてもよく、他の画像処理がWSIからのパッチ抽出後に実行されてもよい。
これらの画像処理ステップは例として説明されており、本発明の範囲を制限するものとして解釈されるべきではない。例えば、十分な処理能力が利用できる場合、CNNはカラー画像を直接に操作できる。
訓練および予測
図4は、CNNの訓練に含まれるステップを示すフロー図である。
ステップS40では、腫瘍を発見し、輪郭を描画し、分類するために臨床医によって注釈が付された処理用のWSIを含む訓練データが取得される。臨床医の注釈は、グラウンドトゥルースデータを表している。
ステップS41では、WSIが、CNNの入力画像パッチである画像パッチに分解される。つまり、画像パッチが、WSIから抽出される。
ステップS42では、画像パッチが、上記のように前処理される(代替的にもしくは付加的に、WSIは、ステップS41の前に上記のように前処理されうる)。
ステップS43では、CNNの重みすなわち層間の重みに対して初期値が設定される。
図1Aおよび図1Bを参照してさらに上述したように、ステップS44では、入力画像パッチのバッチの各々がCNNに入力され、ピクセルごとにパッチを発見し、輪郭を描画し、分類するために処理される。発明者らの方法は各腫瘍(または腫瘍型)ピクセルを特定するためのものであるため、本明細書の用語の概要は必ずしも厳密に技術的に正しい用語として使用されるとは限らず、したがって、CNNにより各腫瘍型の腫瘍領域が求められるというほうがおそらくより正確である。
ステップS45では、CNN出力画像パッチがグラウンドトゥルースデータと比較される。このことは、パッチごとに実行されうる。代替的にWSI全体をカバーするパッチが抽出されている場合、このことは、WSIレベルで、またはパッチの連続したバッチで構成されるWSIのサブエリア(例えば、WSIの1つの象限)で実行されうる。かかる変形例では、出力画像パッチは、WSI全体またはその隣接部分の確率マップに再アセンブルでき、確率マップは、例えば確率マップがWSIへの半透明のオーバレイとしてディスプレイに表示される場合、コンピュータおよびユーザの両方によって、グラウンドトゥルースデータと視覚的に比較できる。
次に、ステップS46で、CNNは、この比較から学習を行い、例えば勾配降下法を用いてCNNの重みを更新する。このようにして、学習は、図4にプロセスフローのリターンループによって示すように、訓練データの繰り返し処理にフィードバックされ、CNNの重みを最適化することができる。
訓練後、CNNは、グラウンドトゥルースデータとは関係なくWSIに適用できる。つまり、CNNは、予測のためにライブで使用できる。
図5は、CNNを使用した予測に含まれるステップを示すフロー図である。
ステップS50では、例えば、検査情報システム(LIS)または他の組織学的データリポジトリから、1つ以上のWSIが処理のために取得される。WSIは、例えば上記のように前処理される。
ステップS51では、画像パッチが、WSIから、または各WSIから抽出される。パッチは、WSI全体をカバーする場合もあれば、ランダムまたは非ランダムの選択である場合もある。
ステップS52では、画像パッチが例えば上記のように前処理される。
図1Aおよび図1Bを参照してさらに上述したように、ステップS53では、入力画像パッチのバッチの各々はCNNに入力され、ピクセルごとにパッチを発見し、輪郭を描画し、分類するために処理される。次に、出力パッチは、入力画像パッチが抽出されたWSIの確率マップとして再アセンブルできる。確率マップは、例えば、確率マップがWSIへの半透明のオーバレイとしてディスプレイに表示されている場合、またはWSIと並んでディスプレイに表示されている場合に、コンピュータ装置によるデジタル処理だけでなく、ユーザが視覚的にWSIと比較することもできる。
ステップS54では、偽陽性である可能性が高い腫瘍、例えば、小さすぎる領域またはエッジアーチファクトである可能性がある領域を除外して、腫瘍領域がフィルタリングされる。
ステップS55では、スコアリングアルゴリズムが実行される。スコアリングは細胞特異的であり、スコアは腫瘍ごとに集計され、かつ/またはWSI(あるいはWSIのサブエリア)についてさらに集計することができる。
ステップS56では、結果が、例えば適切な高解像度モニタ上に注釈付きWSIを表示することによって、診断のために、病理医または他の関連する熟練した臨床医に提示される。
ステップS57では、CNNの結果、すなわち確率マップデータ、および任意にCNNパラメータに関連するメタデータ、ならびに病理医によって追加された追加の診断情報が、CNNによって処理されたWSIまたはWSIのセットを含む患者データファイルにリンクされた方法で保存される。したがって、LISまたはその他の組織学的データリポジトリ内の患者データファイルは、CNNの結果により補完される。
辺縁の評価
図6は、CNNを使用した予測に含まれるステップを示すフロー図である。
ステップS1では、処理のために、染料で着色された組織サンプルと染色した組織サンプルのうちの1つ以上のWSIとを取得し、例えば、検査情報システム(LIS)を使用して、仮想スライドライブラリまたはその他の組織学的データリポジトリからデータを取得する。
ステップS2では、WSIは、例えば上記のように前処理される。
ステップS3では、第1のCNNを適用して、3方向の組織型分類(非腫瘍、浸潤性腫瘍、非浸潤性腫瘍)を実行する。画像パッチは、1つ以上のWSIから抽出される。パッチは、WSI全体をカバーする場合もあれば、ランダムまたは非ランダムの選択である場合もある。画像パッチは、例えば上記のように前処理することができる。入力画像パッチのバッチの各々は、CNNに入力され、図1Aおよび図1Bを参照してさらに上述したように、ピクセルごとにパッチを分類するために処理される。出力パッチは、腫瘍の位置の確率マップとして再アセンブルされる。
ステップS4では、第2のCNNが適用されて、バイナリ分類を実行して、染料ピクセル(染料、非染料)が認識される。ステップS3について説明したように確率マップが生成され、ここで、確率マップにより、染料の位置、すなわち組織塊の境界が特定される。
ステップS5では、2つの確率マップの後処理を実行して、単純な閾値処理またはより複雑な数値処理に従って確率値がフィルタリングされる。例えば、セグメンテーションを確率値に適用して染料境界を特定し、これによって染料カテゴリのピクセルの確率値を調整することができる。後処理では、セグメンテーションを使用して、偽陽性の可能性が高い腫瘍領域、例えば小さすぎる領域またはエッジアーチファクトである可能性がある領域を除外することもできる。必要に応じて、ステップS3およびS4から出力された2つの確率マップを単一の確率マップにマージまたは結合し、マージされた確率マップに後処理を適用することができる。
上記でさらに説明したように、ステップS3およびS4は、単一のCNNが6方向分類を実行し、単一の確率マップを生成する他の実施形態で効果的にマージされる。この場合、ステップS5の後処理は、単一の確率マップで実行される(確率マップはまた、例えば、図6のプロセスフローから出力された辺縁評価データとともに病理医に提示されうる組織学的画像のオーバレイとして、観察アプリケーションにおいてビジュアライゼーションを生成するために使用されうることに留意されたい)。
ステップS6では、各ピクセルを6つのクラスのうちの1つに割り当てることによってマスクが作成される(マスクはまた、例えば、染料境界オーバレイまたは腫瘍ピクセルオーバレイなどの組織学的画像の1つ以上のオーバレイとして、観察アプリケーションでビジュアライゼーションを生成するために使用されうることにも留意されたい。これらのオーバレイは、図6のプロセスフローから出力された辺縁評価データと一緒に病理医に提示することができる)。
ステップS7において、腫瘍染料クラスに閾値数εを超えるピクセルが存在するかどうかが判別される。「ノー」の場合、辺縁状態はステップS8で「陰性」に設定され、プロセスフローは辺縁距離の計算に進む。一方、「イエス」の場合、辺縁状態はステップS12で「陽性」に設定され、プロセスフローは辺縁病変の計算、すなわち腫瘍組織が切除手術によって切除された組織塊の境界に浸潤した程度の計算に進む。閾値テキストは、全ての腫瘍染料クラスのピクセルの合計に単一の閾値を適用することも、各腫瘍染料クラスの個別の合計に個別の閾値を適用することもできる。例えば、後者の場合、「浸潤性腫瘍染料」または「染料上の非浸潤性腫瘍」クラスのいずれかがそれぞれの(同じか、または異なりうる)閾値よりも高いレベルで検出された場合、閾値テスト結果は「イエス」であり、辺縁状態は「陽性」となり、一方、両方の閾値の比較でレベルが閾値以下である場合、テスト結果は「ノー」であり、辺縁状態は「陰性」となる。閾値は、経験的に定義された値またはその他の適切な値に設定できる。
陰性の分岐に続いて、ステップS9で辺縁距離が計算され、ステップS10で適切な閾値に対してテストされる。染料上腫瘍クラスのピクセルが閾値ε未満のレベルで検出された場合、辺縁が存在するため、辺縁状態は、辺縁距離に応じて「クリア」または「クロース」のいずれかである必要がある。ステップS11で、辺縁距離が閾値を下回っている場合、状態は陰性/クロースに設定される。ステップS14において、辺縁距離が閾値を下回っている場合、状態は陰性/クリアに設定される。辺縁距離は、腫瘍組織と染料表面との間の最も近い距離として定義される。集合S1を全ての腫瘍ピクセルの座標(xt,yt)∈Rを含むものとして定義し、集合S2を全ての辺縁ピクセルの座標(xm,ym)∈Rを含むものとして定義できる。したがって、辺縁距離を計算する問題は、2つの集合の平面点間の最小距離を計算する問題となる。この測定値はO(nLogn)時間で計算できるため、大きな点の集合でもすばやく計算できることが示されている(Toussaint and Bhattacharya(1983))。辺縁距離Dmが、ユーザによって設定されうる指定の閾値θを下回っている場合、辺縁状態は「クロース」として返される。Dm>θの場合、辺縁状態は「クリア」として返される。もちろん、実際に計算された辺縁距離が出力されうる。これは、最小距離を表す単純な数値、または最小辺縁距離および最大辺縁距離の個別の距離、平均辺縁距離などの統計値、ならびに標準偏差などの確率分布などの追加データを提供するより高度な出力でありうる。陽性の分岐に続いて、ステップS13で辺縁病変が計算される。腫瘍染料クラスのピクセルは、閾値ε、またはそのうちの1つ以上のレベルで検出されているため、腫瘍のない辺縁は存在せず、陽性の辺縁が存在することになる。陽性の辺縁の場合、辺縁病変は、染料表面の腫瘍細胞の全長として計算できる。これは、各腫瘍に属する2つのピクセル間の最大距離を見出し、結果を合計することと同じである。これは最小距離の問題と同等の問題であるため、O(nLogn)時間で解決することもできる(Vaidya(1989))。現在の米国病理学会のガイドラインでは、おおよその範囲を報告することしか求められていない。全長に基づいて辺縁病変の程度を分類するための1つのスキームは、次の通り、すなわち
-単巣性:最大距離が閾値未満、例えば4mm未満の辺縁にある癌腫の1つの巣性領域
-多巣性:辺縁にある2つ以上の癌腫の病巣で、最大距離の合計が単巣性の閾値未満
-最小/中程度:単巣性または多巣性の場合と同様に、最大距離の合計が、単巣性または多巣性の分類に採用された閾値よりも大きい第2の閾値(例えば5mm未満)よりも小さい
-広範囲:広範な前部の辺縁に癌腫が存在する(最大距離の合計>5mm)
である。
必要に応じて、単巣性および多巣性を、「巣性」というラベルの付いた単一のクラスにマージできる。
最後に、プロセスフローはステップS15で終了し、辺縁状態が報告される。つまり、辺縁状態がディスプレイなどに出力され、かつ/またはファイルに保存される。
要約すると、切除によって抽出され、その表面領域を強調表示するために染料で着色された腫瘍塊から得られた組織切片の組織学的画像を画像処理するためのコンピュータ自動化方法が説明された。1つまたは2つの畳み込みニューラルネットワークは、染料および非染料ならびに腫瘍および非腫瘍を区別する組織クラスを区別するように訓練されている。訓練されたニューラルネットワークは、組織学的画像に適用され、ピクセルが組織クラスに割り当てられている出力画像を生成する。次に、組織切片の陽性または陰性の辺縁状態は、腫瘍染料分類ピクセルの有無に応じてそれぞれ判別することができる。さらに、陽性および陰性の辺縁状態をさらに調べて、分類で指定されたピクセル間距離に基づいて追加のパラメータを計算することにより、辺縁病変(陽性のサブ状態)および辺縁距離(陰性のサブ状態)のサブ状態を生成できる。
腫瘍領域のフィルタリングおよびランク付け
図7は、本開示の一実施形態によるフロー図である。
ステップS71では、スライドスキャナによって生成された可能性があるWSIの画像データを含む画像データファイルが提供される。画像データファイルには、複数の画像が含まれていてよく、例えば、複数の染色の各々に1つ、または有限の深度の透明または半透明のサンプルを通して顕微鏡の焦点面をステップすることによって得られたサンプル(いわゆるzスタック)の異なる深度の各々に1つの画像が含まれうることが理解されよう。
ステップS72は、上記の例としてさらに説明したように、分散ベースの分析、適応閾値処理、モルフォロジー演算など、幾つかのCNN前処理を実行しうる任意のステップである。
ステップS73では、特に図5のステップS51~S54を参照して説明したように、上記のCNNが実行される。組織型のピクセルごとの分類を実行して腫瘍ピクセルをマーキングし、続いて、腫瘍(すなわち腫瘍領域)の輪郭を描画するためにセグメンテーションが実行される。組織型は、癌の種類による分類である。セグメンテーションの場合、概して、隣接する腫瘍ピクセル、すなわち互いに接触しているかまたは互いに近接しているピクセルが、共通の腫瘍に属していると考えられる。ただし、信頼性を向上させるために、例えば2つの異なる癌細胞分類に関連付けられた、異なるピクセル分類を有する2つの接触する腫瘍を同定するために、通常、より複雑なセグメンテーション基準が含まれている。CNNにより、各ピクセルに確率が割り当てられる。これは、CNNが検出するように訓練されたN個のクラスの各々に属するピクセルの確率を表す確率ベクトルである。例えば、浸潤性領域、非浸潤性領域、および非腫瘍領域を区別するように訓練されたCNNの場合、各ピクセルには長さ3のベクトルが割り当てられる。位置kのピクセルは、確率ベクトル[0.1,0.2,0.7]を有してもよく、これは、ピクセルが浸潤性領域にある確率が10%、非浸潤性領域にある確率が20%、非腫瘍領域にある確率が70%であることを示す。
ステップS74では、ステップS73で腫瘍を発見するCNNによって生成されたデータ、すなわち腫瘍固有のデータを使用して、各腫瘍の要約統計量のセットを計算する。例えば、各腫瘍について、スコアは、その腫瘍(領域)に含まれる全てのピクセルについての上記の確率値の数学的平均として計算されうる。中央値、加重平均などの他の要約統計量もスコアの計算に使用できる。他の測定値、例えば、腫瘍内のピクセル数または腫瘍領域の形状によって測定される腫瘍領域、または浸潤性腫瘍および非浸潤性腫瘍などの特定のピクセル分類の有病率などの腫瘍の寸法属性またはモルフォロジー属性を使用することができる。通常、各腫瘍について、腫瘍の確率、腫瘍の面積、および腫瘍の最大寸法の長さの平均および標準偏差が含まれる。腫瘍領域は、必ずしも単一のスライドからの領域ではなく、別々のスライドに属していてもよい。例えば、2つのスライドの組織サンプルは異なる染色で染色され、異なるクラスの腫瘍細胞が強調表示されうるため、幾つかの腫瘍は第1のスライドで同定され、他の腫瘍は第2のスライドで同定される。任意に、より高度なスコアを計算できる。例えば、組織学的画像データ(すなわち画像データで同定された腫瘍)と患者固有の遺伝子(すなわちゲノム)データとの組み合わせで訓練されたCNNを使用して、患者のリスクを予測することができる。この種のCNNは、Mobadersany et al.(2018)によって説明されている。他の実装形態では、スコアは、CNNによって同定された腫瘍に適用される従来の画像処理技術(例えば、セグメンテーションに関連して上記にさらに列挙されている)を使用して計算することができる。例えば、形状およびテクスチャの測定値を遺伝子データと組み合わせて、要約統計量に含める一連の統計測定値を作成できる。サポートベクターマシンまたはランダムフォレストアルゴリズムは、これらの特徴を使用して転移のリスクを予測することができる。転移のリスクは、この腫瘍の細胞が身体の他の部分に転移する確率として定義される。いずれにせよ、転移リスクスコアが計算され、各腫瘍領域に関連付けられる。
ステップS75では、ステップS74で計算および編集された要約統計量に基づいて腫瘍のフィルタリングが実行される。フィルタは、要約統計量のパラメータの1つの値、あるいは要約統計量に存在するパラメータの複数のものを含む論理的な組み合わせおよび/または式から導出された複合値と、閾値とを比較し、腫瘍のリストを単純に通過することで動作することができ、比較結果に基づいて腫瘍物体を削除または保持する。例えば、病理医によって構成されうるか、または事前設定されうるフィルタにより、閾値、例えば、閾値50%を超えた平均確率を有する、100マイクロメートルを超える最大寸法を有する腫瘍のみを通過させることが選択されてもよい。別の例のフィルタは、浸潤性として分類された腫瘍のみを通過させ、浸潤性腫瘍の中で、平均確率が例えば80%を超えるもののみを通過させる。
ステップS76では、ステップS75でフィルタを通過した腫瘍のセットが取得され、順番にランク付けされる。次に、要約統計量に基づいて、一連の基準を適用して、腫瘍のランク付けを実行できる。標準的なランク付けアプローチのプリセットがユーザに提供され、ユーザがプリセットの1つを選択できるようにされうる。さらに、どの基準が適用されるかを定義するためのユーザインタフェースがユーザに提供されうる。例えば、長さの寸法または面積、または整数カウント(例えば、セル数)に基づくスカラー値を基準とする場合、ユーザは、これらの基準の閾値または値の範囲を設定できる。ランク付け順は、患者の生存の重要性を示す複合スコア(5年生存率など)に基づく場合もあれば、単純な単一パラメータのランク付けに基づく場合もある。例えば、面積または最大寸法などの腫瘍のサイズパラメータ、または真円度などのモルフォロジーパラメータに基づいてもよい。
ステップS77では、スライド画像のビジュアライゼーションが生成され、それが表示デバイスのGUIウィンドウでユーザに表示される。ビジュアライゼーションでは、フィルタリングおよびランク付けも考慮される。特に、WSI内の最も重要な腫瘍とみなされる領域は、臨床的に関連性のある方式でその要約統計量とともに表示され、当該要約統計量は、全ての場合に全体的な表示の一部として、またはユーザからのGUIコマンド、例えば、WSIの関心腫瘍のカーソルでの「クリック」に選択的に応答して利用でき、その腫瘍の統計を表形式またはその他の適切な形式で表示するポップアップウィンドウが生成される。このアプローチにより、潜在的に重要な腫瘍を強調表示して、潜在的に重要な腫瘍間のランク付け情報およびこれらの腫瘍の統計的要約をユーザに提供する方法で、ユーザに提示することができる。統計的要約により、各腫瘍の個々のパラメータ値、特にフィルタリングおよび/またはランク付け基準として使用されるもの、ならびに複数のフィルタリングおよび/またはランク付け基準の公式および/またはブール論理積から計算されるランク付け番号または有意性スコアなどの合成パラメータが提示されうる。
通常、表示される画像は、オーバレイビューまたはマルチタイルビューの組み合わせのいずれかの形態になる。オーバレイビューでは、(処理済みの可能性がある)生データが表示され、セグメンテーションデータがその上にオーバレイされる。セグメンテーションデータは、ビジュアライゼーションのために各腫瘍のシェーディングおよび/または輪郭に変換される。シェーディングまたは輪郭付けは、例えば腫瘍分類によって色分けすることができる。組織の非腫瘍領域は、まったくマーキングされていなくてもよく、または透明度の高いカラーウォッシュ(ブルーウォッシュなど)でシェーディングされてもよい。マルチタイルビューでは、オーバレイビューの種々の層がタイルとして並べて表示されるため、フィルタリングされた腫瘍領域の(処理済みの可能性がある)生の画像データおよびセグメンテーションデータを示すタイルが表示される。必要に応じて、腫瘍分類タイプごとに個別のセグメンテーションデータタイルを表示できる。ディスプレイ内の腫瘍領域の提示では、ステップS75で実行されるフィルタリング、および好ましくはステップS76で実行されるランク付けも考慮に入れられる。腫瘍スコア、分類された腫瘍型、および腫瘍領域に関連する要約統計量の他のパラメータなどの要素を単独で、または組み合わせて使用して、表示を構成することができる。
検出され、フィルタリングされ、ランク付けされた腫瘍をユーザに表示するための複数の選択肢がある。
WSIビューにおいて、腫瘍情報を表示する1つの方法は、サムネイル画像および/またはテキスト情報を含みうる一連のマーカーをWSIの低解像度の画像にオーバレイすることである。マーカーは、ランク付けで並べ替えることができ、またはランク付けを明示的に示すこともできる。病理医には、マーカーを選択することで腫瘍を選択するための、また任意に、ランクの高い腫瘍からランクの低い腫瘍へまたは逆にランクの低い腫瘍からランクの高い腫瘍への両方でのランク付けによって腫瘍をナビゲートするための、適切なGUIツールが提供される。適切なGUIツールとして、キーボードショートカット、キーボードの上下もしくは左右の矢印キー、キーボードのページアップキーおよびページダウンキー、マウスナビゲーション(スクロールホイールで上下にスクロールするなど)またはその他の入力デバイス(音声ナビゲーション、タッチセンサでのマルチタッチジェスチャなど)が挙げられる。GUIでランク付けマーカーを選択すると、その腫瘍の要約統計量もしくはそのサブセットの表示、および/または低解像度の画像での腫瘍の表示にリンクされたポップアップ、あるいは別の高解像度の表示ペインでの腫瘍の高解像度ビューの表示が促されうる。ユーザは、キーボードコマンド、マウスコマンド(例えば、ダブルクリック)または他の適切な入力などの適切なGUIコマンドによって、腫瘍の高解像度の、例えば完全なネイティブの解像度の画像のディスプレイを引き起こすことができる。低解像度のWSIビューに戻るための対応するGUIコマンドが用意されている。ビジュアライゼーションアプリケーションにより、好ましくは、GUI制御が提供され、これにより、ユーザは、低解像度および高解像度の両方のビューで順序をランク付けして、腫瘍を上下にステップスルーさせることができる。
WSIビューの一例は、ステップS76のフィルタを通過した全ての腫瘍を、ランク付けマーカーラベル1、2、3などとともに表示する(つまり、そのセグメンテーションデータを表示する)ことである。次に、ランク付けマーカー(または腫瘍)のラベルをクリックすると、選択された要約統計量のセット、特にフィルタで使用される統計量および/または腫瘍のサムネイルビューをWSIビューより高い解像度で一覧表示するポップアップが生成されうる。代替的に、ビューは、画面の一部にオーバレイ画像または画像タイルが表示され、画面の別の部分にフィルタリングされた腫瘍のテーブルが表示される、分割画面ビューであってもよい。テーブルは最初にランク付けして並べ替えて表示できるが、GUIは、他の列または複数の列の組み合わせでユーザが並べ替えることができる機能も有しうる。ここで、他の列とは、腫瘍領域、腫瘍分類、転移確率など、要約統計量またはフィルタ基準からの任意の基準でありうる。例えば、腫瘍分類による並べ替えの後、腫瘍領域による並べ替えを行うことができる。マルチスライド画像ファイルの場合、スライド番号による並べ替えの後、他のパラメータで並べ替えを行うことができる。
多重解像度ビューには、通常はWSIを再現する低解像度(10倍の倍率など)の1つの表示ペインと、高解像度(例えば、倍率60倍、つまり10倍ビューに対して6倍または600%ズーム)の別の表示ペインとが含まれる。例えば、表示される初期ビューは、WSIの低解像度の画像ペインおよび上位の腫瘍領域を中心とした高解像度の画像ペインでありうる。次に、ステップスルー下矢印(もしくは下矢印および上矢印のペア)または他の適切なGUIボタンもしくはボタンの組み合わせ(物理的もしくは仮想スクロールホイールなど)により、ユーザは、フィルタリングされた腫瘍領域を1つずつランク付けしてステップスルーすることができる。GUIを使用することで、ユーザは、ユーザ入力を介して高解像度の画像ペインの解像度を調整できる。また、GUIによって、腫瘍が高解像度の表示ペインを実質的に満たすようなサイズ設定が行われるように、高解像度の画像ペインの初期解像度を選択することができる。
このように、ビジュアライゼーションアプリケーションにより、ユーザに表示されるもの、セグメンテーションデータおよび要約統計量で強調表示される腫瘍、ならびに選択ビューの場合、腫瘍がユーザに表示される時系列、つまり順序が求められる。
スライドセットに特に適した腫瘍領域を表示する別の方法は、腫瘍の低解像度の画像およびその要約統計量が並べ替えられたタイルとして表示されるモンタージュを作成することである。並べ替えは、1次元(1D)リストまたは2次元(2D)グリッドを表示する方法で行うことができる。バーチャルリアリティゴーグルを使用して、腫瘍を3次元で配置することも可能である。ユーザは、キーストローク、マウススクロール、またはその他の入力モダリティ(音声、タッチなど)を使用して、1Dリストまたは2Dタイルアレイをナビゲートできる。タイルが選択されると、病理医は、適切なスライドで腫瘍の高解像度バージョンにすばやく移動し、要約統計量およびセグメンテーションデータによって提示されたCNN分析の結果と比較して、任意に腫瘍のより詳細な検証を行うことができ、また、任意に、診断に役立ちうる1つ以上の追加の分析アルゴリズムを腫瘍に適用することにより、腫瘍に対して追加の処理を実行することができる。
ステップS78において、フィルタリングされランク付けされた腫瘍を観察した後、病理医は、さらなる調査のためにこれらの腫瘍のいずれか(またはランク付けされた腫瘍のサブセット、もしくは実際にフィルタリングされた腫瘍の完全なセット)を選択して追加の調査を行う選択肢を有する。
ステップS79では、この追加の調査が実行される。追加の調査とは、1つ以上の追加アルゴリズムを適用して、選択した腫瘍に関する詳細情報を提供し、それらの診断を支援することを意味する。これは、暫定的な診断を確認するのに役立ち、または病理医が出発点として検討するための初期診断を行うのに役立ちうる。
例えば、乳癌の場合、病理医が、ステップS78で浸潤性乳癌を含む腫瘍を選択することができ、ステップS79で腫瘍の有糸分裂数を計算するアルゴリズムを適用する。つまり、ビジュアライゼーションアプリケーションには、かかるアルゴリズムが含まれている。IHC染色スライドからの画像データに関連する別の乳癌の例では、アルゴリズムが適用されて、選択した腫瘍におけるER、PR、またはHER2などの1つ以上の診断に関連する遺伝子の発現のスコアが計算される。
要約すると、ビジュアライゼーションは、潜在的に重要な腫瘍を含む画像内の領域を調べるようにユーザを誘導するように構成されている。特に、従来のビジュアライゼーションアプリケーションと比較して、ユーザは、潜在的な腫瘍領域に遭遇するたびに、低解像度でスライド領域全体を手動で視覚的に走査し、高解像度のビューにズームインおよびズームアウトする必要がない。したがって、上記のようにCNN支援フィルタリングとランク付けとに基づいて病理医に提示されたビジュアライゼーションの自動前処理により、病理医がスライドをレビューして適切な診断を行うために必要な時間が短縮されるだけでなく、人為的ミスによって重要な腫瘍を見落とす可能性を減らすことができる。
例示的な実施形態
一実施形態では、組織学的画像内の腫瘍を同定するためのシステムは、組織学的画像内の腫瘍を同定する方法を実行するように構成される。この方法は、ピクセルの2次元アレイを含む組織学的画像を受信することと、畳み込みニューラルネットワークを適用して、ピクセルの2次元アレイを有する出力画像を生成し、出力画像内のピクセルを組織学的画像のピクセルにマッピングすることと、を含む。出力画像は、複数の組織クラスのうちの1つを各ピクセルに割り当てることによって生成され、複数の組織クラスは、非腫瘍組織を表す少なくとも1つのクラスおよび腫瘍組織を表す少なくとも1つのクラスを含む。
次に、システムにより、出力画像からセグメンテーションマスクが生成される。セグメンテーションマスクにより、個々の腫瘍によって占有された領域が特定される。システムはまた、各腫瘍の要約統計量を計算し、システムは、セグメンテーションマスクをフィルタリングして、各腫瘍の要約統計量に従って個々の腫瘍を選択および選択解除するように構成されている。システムはまた、編集されたセグメンテーションマスクに従って組織学的画像のビジュアライゼーションを作成し、フィルタリングされたセグメンテーションマスク内の個々の腫瘍の選択を受信し、次に、現在選択されている腫瘍に対して追加の計算診断プロセスの実行を調整するように構成されている。
一実施形態では、追加の計算診断プロセスは、処理タスクを追加の計算診断プロセスに送信することと、処理タスクの結果を追加の計算診断から受信することと、を含む。追加の計算診断プロセスは、ローカルで実行することも、データ通信ネットワークを介して通信可能に結合されたリモートプロセッサ対応デバイスで実行することもできる。処理タスクの一部には、組織学的画像が含まれうる。処理タスクはまた、組織学的画像に関連するメタデータを含んでもよく、状況によっては、組織学的画像に関連するメタデータは、セグメンテーションマスクを含みうる。さらに、組織学的画像に関連するメタデータはまた、患者情報、染色情報、染色プロトコル情報、および走査プロトコル情報のうちの1つ以上を含みうる。
CNNコンピューティングプラットフォーム
提案している画像処理は、多種多様なコンピューティングアーキテクチャ、特に、CPU、GPU、TPU、FPGA、および/またはASICに基づくニューラルネットワーク用に最適化されたコンピューティングアーキテクチャで実行されうる。幾つかの実施形態では、ニューラルネットワークは、テスラK80 GPUなどの、カリフォルニア州サンタクララ在のNvidia CorporationのNvidia GPU上で実行されるGoogleのTensorflowソフトウェアライブラリを使用して実装される。他の実施形態では、ニューラルネットワークは、汎用CPU上で実行することができる。より高速な処理は、CNN計算を実行するための目的に合わせて設計されたプロセッサ、例えば、Jouppi et al.(2017)に開示されたTPUによって得ることができ、その内容は全て、参照により本明細書に組み込まれている。
図8は、Jouppi et al.(2017)のTPUを示しており、Jouppiの図1を簡略化したものである。TPU100は、256×256MACを含む収縮行列乗算ユニット(MMU)102を有しており、この256×256MACにより、符号付きまたは符号なし整数に対する8ビットの乗算および加算が実行されうる。MMUの重みは、重みFIFOバッファ104を介して供給され、次に、適切なメモリインタフェース108を介して、オフチップ8GB DRAMの形態で、メモリ106から重みを読み取る。中間結果を保存するために、統合バッファ(UB)110が提供される。MMU102は、重みFIFOインタフェース104およびUB110からの入力を(収縮期データセットアップユニット112を介して)受信し、MMU処理の16ビット積をアキュムレータユニット114に出力するように接続されている。アクティベーションユニット116は、アキュムレータユニット114に保持されているデータに対して非線形な機能を実行する。正規化ユニット118およびプーリングユニット120によるさらなる処理の後、中間結果は、データセットアップユニット112を介してMMU102に再供給するためにUB110に送信される。プーリングユニット120は、必要に応じて、最大プーリング(すなわち、マックスプーリング)または平均プーリングを実行することができる。プログラム可能なDMAコントローラ122は、TPUのホストコンピュータとUB110との間でデータを転送する。TPU命令は、ホストインタフェース124および命令バッファ126を介して、ホストコンピュータからコントローラ122に送信される。
ニューラルネットワークを実行するために使用される計算能力は、CPU、GPU、またはTPUのいずれに基づくものであっても、例えば後述するクリニカルネットワークでローカルにホストされうるものであり、またはデータセンタでリモートでホストされうるものであることが理解されよう。
ネットワークおよびコンピューティングおよびスキャン環境
提案しているコンピュータ自動化方法は、より大きなクリニカルネットワーク環境の一部である検査情報システム(LIS)、例えば、病院情報システム(HIS)または画像アーカイブおよび通信システム(PACS)などのコンテキストで動作する。LISでは、WSIは、データベース、通常は個々の患者の電子医療記録を含む患者情報データベースに保持される。WSIは、スライドにマウントされた染色組織サンプルから取得される。WSIを取得する顕微鏡にはバーコードリーダが装備されているため、スライドには、WSIに適切なメタデータがタグ付けされたバーコードラベルが印刷される。ハードウェアの観点からは、LISは、必要に応じて有線および無線接続を備えたローカルエリアネットワーク(LAN)などの従来のコンピュータネットワークとなる。
図9は、本発明の実施形態と併せて使用することができる例示的なコンピュータネットワークを示す図である。ネットワーク150は、病院152内のLANを含む。病院152は、複数のワークステーション154を備えており、複数のワークステーション154は、それぞれがローカルエリアネットワークを介して、関連するストレージデバイス158を有する病院コンピュータサーバ156にアクセスすることができる。LIS、HISまたはPACSアーカイブは、ストレージデバイス158に保存され、これにより、アーカイブ内のデータは、任意のワークステーション154からアクセスすることができる。ワークステーション154のうちの1つ以上は、グラフィックカードおよびソフトウェアにアクセスして、前述の画像生成方法をコンピュータで実装することができる。ソフトウェアは、またはいずれかのワークステーション154または各ワークステーション154にローカルに保存可能であるか、またはリモートにて保存可能であり、必要に応じてネットワーク150を介してワークステーション154にダウンロードされうる。他の例では、本発明を具現化する方法は、端末として動作するワークステーション154を備えたコンピュータサーバ上で実行されうる。例えば、ワークステーションは、所望の組織学的画像データセットを定義するユーザ入力を受信し、CNN分析がシステムの他の場所で実行されている間、結果の画像を表示するように構成されうる。また、複数の組織学的および他の医療撮像デバイス160、162、164、166が、病院のコンピュータサーバ156に接続されている。デバイス160、162、164、166で収集された画像データは、ストレージデバイス156上のLIS、HIS、またはPACSアーカイブに直接に保存することができる。したがって、組織学的画像は、対応する組織学的画像データが記録された直後に観察および処理することができる。ローカルエリアネットワークは、病院インターネットサーバ170によってインターネット168に接続されており、これにより、LIS、HIS、またはPACSアーカイブへのリモートアクセスが可能になる。これは、データへのリモートアクセスや、例えば患者が移動した場合などに病院間でデータを転送したり、外部調査を実施したりする場合に役立つ。
図10は、本明細書に記載の種々の実施形態に関連して使用されうる例示的なコンピューティング装置500を示すブロック図である。例えば、コンピューティング装置500は、上記のLISまたはPACSシステムにおけるコンピューティングノード、例えば、適切なGPU、または図8に示されるTPUと組み合わせてCNN処理が実行されるホストコンピュータとして使用されてもよい。
コンピューティング装置500は、サーバまたは任意の従来のパーソナルコンピュータ、あるいは有線または無線のデータ通信が可能な他の任意のプロセッサ対応デバイスでありうる。当業者には明らかであるように、有線または無線のデータ通信が不可能なデバイスを含む、他のコンピューティング装置、システム、および/またはアーキテクチャも使用することができる。
コンピューティング装置500は、好ましくは、プロセッサ510などの1つ以上のプロセッサを含む。プロセッサ510は、例えば、CPU、GPU、TPU、またはアレイ、あるいはこれらの組み合わせであってよく、例えば、CPUとTPUとの組み合わせ、またはCPUとGPUとの組み合わせであってよい。入出力を管理するための補助プロセッサ、浮動小数点数学演算を実行するための補助プロセッサ(例えば、TPU)、信号処理アルゴリズム(デジタルシグナルプロセッサ、イメージプロセッサなど)の高速実行に適したアーキテクチャを備えた専用マイクロプロセッサ、メイン処理システムに従属するスレーブプロセッサ(バックエンドプロセッサなど)、デュアルもしくはマルチプロセッサシステム用の追加のマイクロプロセッサもしくはコントローラ、またはコプロセッサなどの追加のプロセッサを提供することができる。かかる補助プロセッサは、個別のプロセッサでありうるか、またはプロセッサ510と統合されうる。コンピューティング装置500で使用できるCPUの例は、Pentiumプロセッサ、Core i7プロセッサ、およびXeonプロセッサであり、これらは全て、カリフォルニア州サンタクララ在のインテルコーポレーションから市販されている。コンピューティング装置500で使用されうるGPUの例は、カリフォルニア州サンタクララ在のNvidia CorporationのテスラK80 GPUである。
プロセッサ510は、通信バス505に接続されている。通信バス505には、ストレージとコンピューティング装置500の他の周辺コンポーネントとの間の情報転送を容易にするためのデータチャネルが含まれうる。通信バス505はさらに、プロセッサ510との通信に使用される信号のセットを提供することができ、これには、データバス、アドレスバス、および制御バス(図示せず)が含まれる。通信バス505には、任意の標準または非標準バスアーキテクチャが含まれていてよく、例えば業界標準アーキテクチャ(ISA)、拡張業界標準アーキテクチャ(EISA)に準拠するバスアーキテクチャなど、マイクロチャネルアーキテクチャ(MCA)、周辺コンポーネントインタコネクト(PCI)ローカルバス、またはIEEE488汎用インタフェースバス(GPIB)、IEEE696/S-100などを含む、米国電気電子学会(IEEE)によって公布された規格が含まれていてよい。
コンピューティング装置500は、好ましくはメインメモリ515を含み、セカンダリメモリ520も含みうる。メインメモリ515は、上記で論じた機能および/またはモジュールのうちの1つ以上など、プロセッサ510上で実行されるプログラムのための命令およびデータの記憶を提供する。メモリに保存され、プロセッサ510によって実行されるコンピュータ可読プログラム命令は、アセンブラ命令、命令セットアーキテクチャ(ISA)命令、機械命令、機械依存命令、マイクロコード、ファームウェア命令、状態設定データ、集積回路の構成データ、もしくはSmalltalk、C/C++、Java、JavaScript、Perl、VisualBasic.NETなどを含むがこれらに限定されない1つ以上のプログラミング言語の任意の組み合わせで記述および/またはコンパイルされたソースコードまたはオブジェクトコードのいずれかでありうることを理解されたい。メインメモリ515は、典型的には、ダイナミックランダムアクセスメモリ(DRAM)および/または静的ランダムアクセスメモリ(SRAM)などの半導体ベースのメモリである。他の半導体ベースのメモリタイプには、例えば、同期動的ランダムアクセスメモリ(SDRAM)、ランバス動的ランダムアクセスメモリ(RDRAM)、強誘電体ランダムアクセスメモリ(FRAM)などが含まれ、読み取り専用メモリ(ROM)が含まれる。
コンピュータ可読プログラム命令は、スタンドアロンソフトウェアパッケージとして、完全にユーザのコンピュータ上で、一部はユーザのコンピュータ上で実行可能であり、一部はユーザのコンピュータ上で、一部はリモートコンピュータ上で、または完全にリモートコンピュータもしくはサーバ上で実行されうる。後者のシナリオでは、リモートコンピュータは、ローカルエリアネットワーク(LAN)またはワイドエリアネットワーク(WAN)を含む任意のタイプのネットワークを介してユーザのコンピュータに接続可能であり、または(例えば、インターネットサービスプロバイダを使用したインターネット経由で)外部コンピュータに接続可能である。
セカンダリメモリ520は、任意に、内部メモリ525および/またはリムーバブル媒体530を含みうる。リムーバブル媒体530は、任意の周知の方法で読み取られ、かつ/または書き込まれる。リムーバブル記憶媒体530は、例えば、磁気テープドライブ、コンパクトディスク(CD)ドライブ、デジタル多用途ディスク(DVD)ドライブ、他の光学ドライブ、フラッシュメモリドライブなどであってもよい。
リムーバブル記憶媒体530は、非一時的コンピュータ可読媒体であり、コンピュータ実行可能コード(すなわち、ソフトウェア)および/またはデータがその上に保存されている。リムーバブル記憶媒体530に保存されたコンピュータソフトウェアまたはデータは、プロセッサ510による実行のためにコンピューティング装置500に読み込まれる。
セカンダリメモリ520には、コンピュータプログラムまたは他のデータもしくは命令をコンピューティング装置500にロードすることを可能にするための他の同様の要素が含まれうる。かかる手段には、例えば、外部記憶媒体545および通信インタフェース540が含まれていてよく、これにより、ソフトウェアおよびデータが外部記憶媒体545からコンピューティング装置500に転送されることが可能となる。外部記憶媒体545の例として、外付けハードディスクドライブ、外付け光学ドライブ、外付け光磁気ドライブなどが挙げられる。セカンダリメモリ520の他の例として、プログラム可能な読み取り専用メモリ(PROM)、消去可能なプログラム可能な読み取り専用メモリ(EPROM)、電気的に消去可能な読み取り専用メモリ(EEPROM)、またはフラッシュメモリ(EEPROMと類似のブロック指向のメモリ)などの半導体ベースのメモリが挙げられる。
上記のように、コンピューティング装置500は、通信インタフェース540を含みうる。通信インタフェース540により、ソフトウェアおよびデータが、コンピューティング装置500と外部デバイス(例えば、プリンタ)、ネットワーク、または他の情報源との間で転送されることが可能となる。例えば、コンピュータソフトウェアまたは実行可能コードは、通信インタフェース540を介してネットワークサーバからコンピューティング装置500に転送されうる。通信インタフェース540の例として、内蔵ネットワークアダプタ、ネットワークインタフェースカード(NIC)、パーソナルコンピュータメモリカード国際協会(PCMCIA)ネットワークカード、カードバスネットワークアダプタ、ワイヤレスネットワークアダプタ、ユニバーサルシリアルバス(USB)ネットワークアダプタ、モデム、ネットワークインタフェースカード(NIC)、ワイヤレスデータカード、通信ポート、赤外線インタフェース、IEEE1394ファイアワイヤ、またはネットワークもしくは別のコンピューティングデバイスによりシステム550とインタフェースを有しうるその他のデバイスが挙げられる。好ましくは、通信インタフェース540により、業界で公布されたプロトコル標準が実装され、好ましくは、イーサネットIEEE802標準、ファイバチャネル、デジタル加入者線(DSL)、非同期デジタル加入者線(ADSL)、フレームリレー、非同期転送モード(ATM)、統合デジタルサービスネットワーク(ISDN)、パーソナル通信サービス(PCS)、伝送制御プロトコル/インターネットプロトコル(TCP/IP)、シリアルラインインターネットプロトコル/ポイントツーポイントプロトコル(SLIP/PPP)などが実装されるが、カスタマイズされた、または非標準のインタフェースプロトコルも実装されうる。
通信インタフェース540を介して転送されるソフトウェアおよびデータは、概して、電気通信信号555の形態である。これらの信号555は、通信チャネル550を介して通信インタフェース540に提供されうる。一実施形態では、通信チャネル550は、有線または無線ネットワーク、あるいは任意の多種多様な他の通信リンクでありうる。通信チャネル550は、信号555を伝送しており、幾つかの例を挙げると、ワイヤもしくはケーブル、光ファイバ、従来の電話回線、携帯電話リンク、無線データ通信リンク、無線周波数(「RF」)リンク、または赤外線リンクを含む多種多様な有線または無線通信手段を使用して実装することができる。
コンピュータ実行可能コード(すなわち、コンピュータプログラムまたはソフトウェア)は、メインメモリ515および/またはセカンダリメモリ520に保存される。コンピュータプログラムはまた、通信インタフェース540を介して受信され、メインメモリ515および/またはセカンダリメモリ520に保存されうる。本明細書の他の箇所で説明されるように、かかるコンピュータプログラムは、実行される際に、コンピューティング装置500に、開示される実施形態の種々の機能を実行させることを可能にする。
本明細書において、「コンピュータ可読媒体」という用語は、コンピュータ実行可能コード(例えば、ソフトウェアおよびコンピュータプログラム)をコンピューティング装置500に提供するために使用される任意の非一時的コンピュータ可読記憶媒体を指すために使用される。かかる媒体の例として、メインメモリ515、セカンダリメモリ520(内部メモリ525、リムーバブル媒体530、および外部記憶媒体545を含む)、および通信インタフェース540と通信可能に結合された任意の周辺デバイス(ネットワーク情報サーバまたは他のネットワークデバイスを含む)が挙げられる。これらの非一時的コンピュータ可読媒体は、実行可能コード、プログラミング命令、およびソフトウェアをコンピューティング装置500に提供するための手段である。ソフトウェアを使用して実装される一実施形態では、ソフトウェアは、コンピュータ可読媒体に保存され、リムーバブル媒体530、I/Oインタフェース535、または通信インタフェース540を介してコンピューティング装置500にロードされうる。かかる一実施形態では、ソフトウェアは、電気通信信号555の形態でコンピューティング装置500にロードされる。ソフトウェアは、プロセッサ510によって実行される際に、好ましくはプロセッサ510に、本明細書の他の箇所で説明されている特徴および機能を実行させる。
I/Oインタフェース535により、コンピューティング装置500の1つ以上のコンポーネントと1つ以上の入力および/または出力デバイスとの間のインタフェースが提供される。入力デバイスの例として、キーボード、タッチスクリーンまたは他のタッチセンシティブデバイス、生体認証デバイス、コンピュータマウス、トラックボール、ペンベースのポインティングデバイスなどが挙げられるが、これらに限定されない。出力デバイスの例として、ブラウン管(CRT)、プラズマディスプレイ、発光ダイオード(LED)ディスプレイ、液晶ディスプレイ(LCD)、プリンタ、真空蛍光表示管(VFD)、表面伝導型電子放出ディスプレイ(SED)、電界放出ディスプレイ(FED)などが挙げられるが、これらに限定されない。
コンピューティング装置500にはまた、音声ネットワークおよび/またはデータネットワークを介した無線通信を容易にする任意の無線通信コンポーネントが含まれる。無線通信コンポーネントには、アンテナシステム570、無線システム565、およびベースバンドシステム560が含まれる。コンピューティング装置500では、無線周波数(RF)信号は、無線システム565の管理下にあるアンテナシステム570によって無線で送受信される。
アンテナシステム570は、アンテナシステム570に送信および受信信号経路を提供するためにスイッチング機能を実行する1つ以上のアンテナおよび1つ以上のマルチプレクサ(図示せず)を含みうる。受信経路において、受信されたRF信号は、低ノイズ増幅器(図示せず)に結合可能であり、低ノイズ増幅器は、マルチプレクサから、受信されたRF信号を増幅し、増幅された信号を無線システム565に送信する。
無線システム565は、種々の周波数で通信するように構成された1つ以上の無線を含みうる。一実施形態では、無線システム565は、復調器(図示せず)および変調器(図示せず)を1つの集積回路(IC)に組み合わせることができる。復調器および変調器は、別々のコンポーネントにすることもできる。着信経路において、復調器は、RF搬送波信号を除去し、ベースバンド受信オーディオ信号を残す。ベースバンド受信オーディオ信号は、無線システム565からベースバンドシステム560に送信される。
受信信号が音声情報を含む場合、ベースバンドシステム560は、信号を復号し、これをアナログ信号に変換する。次に、信号が増幅されてスピーカに送信される。ベースバンドシステム560はまた、マイクロフォンからアナログオーディオ信号を受信する。これらのアナログオーディオ信号は、デジタル信号に変換され、ベースバンドシステム560によって符号化される。ベースバンドシステム560はまた、送信用のデジタル信号をコード化し、無線システム565の変調器部分にルーティングされるベースバンド送信オーディオ信号を生成する。変調器は、ベースバンド送信オーディオ信号をRF搬送波信号と混合して、RF送信信号を生成し、RF送信信号は、アンテナシステム570にルーティングされ、電力増幅器(図示せず)を通過しうる。電力増幅器は、RF送信信号を増幅し、これをアンテナシステム570にルーティングし、そこで、信号は、送信のためにアンテナポートに切り替えられる。
ベースバンドシステム560はまた、プロセッサ510と通信可能に結合されており、プロセッサ510は、中央処理装置(CPU)でありうる。プロセッサ510は、データ記憶領域515および520にアクセスすることができる。プロセッサ510は、好ましくは、メインメモリ515またはセカンダリメモリ520に保存されうる命令(すなわち、コンピュータプログラムまたはソフトウェア)を実行するように構成される。コンピュータプログラムはまた、ベースバンドプロセッサ560から受信され、メインメモリ510またはセカンダリメモリ520に保存されるか、または受信時に実行されうる。かかるコンピュータプログラムは、実行される際に、コンピューティング装置500に、開示される実施形態の種々の機能を実行させることを可能にするためのものである。例えば、データ記憶領域515または520は、種々のソフトウェアモジュールを含んでもよい。
コンピューティング装置は、通信バス505に直接に接続されたディスプレイ575をさらに含み、これは、上記のI/Oインタフェース535に接続された任意のディスプレイに代えて、またはこれに加えて提供されうる。
種々の実施形態はまた、例えば、特定用途向け集積回路(ASIC)、プログラマブルロジックアレイ(PLA)、またはフィールドプログラマブルゲートアレイ(FPGA)などのコンポーネントを使用して、主にハードウェアに実装されうる。本明細書に記載の機能を実行しうるハードウェアステートマシンの実装形態もまた、関連技術分野の当業者には明らかであろう。ハードウェアおよびソフトウェアの両方の組み合わせを使用して、種々の実施形態を実装することもできる。
さらに、当業者は、上記の図および本明細書に開示される実施形態に関連して説明される種々の例示的な論理ブロック、モジュール、回路、および方法ステップが、多くの場合、電子ハードウェア、コンピュータソフトウェア、または両者の組み合わせとして実装されうることを理解するであろう。このハードウェアおよびソフトウェアの互換性を明確に説明するために、種々の例示的なコンポーネント、ブロック、モジュール、回路、およびステップを、概してその機能の観点から上述した。かかる機能がハードウェアとして実装されるか、またはソフトウェアとして実装されるかは、システム全体に課せられる特定のアプリケーションおよび設計上の制約によって異なる。当業者は、特定のアプリケーションごとに種々の方法で説明した機能を実装できるが、かかる実装形態の決定は、本発明の範囲からの逸脱を引き起こすと解釈されるべきではない。さらに、モジュール、ブロック、回路、またはステップ内の機能のグループ化は、説明を容易にするためのものである。特定の機能またはステップは、本発明から逸脱することなく、あるモジュール、ブロック、または回路から別のモジュール、ブロック、または回路に移動することができる。
さらに、本明細書に開示される実施形態に関連して説明される種々の例示的な論理ブロック、モジュール、機能、および方法は、汎用プロセッサ、デジタル信号プロセッサ(DSP)、ASIC、FPGA、または他のプログラマブルロジックデバイス、ディスクリートゲートもしくはトランジスタロジック、ディスクリートハードウェアコンポーネント、あるいは本明細書に記載の機能を実行するように設計されたこれらの任意の組み合わせによって実装または実行することができる。汎用プロセッサはマイクロプロセッサであってよく、代替的に、プロセッサが、任意のプロセッサ、コントローラ、マイクロコントローラ、またはステートマシンであってもよい。プロセッサはまた、コンピューティングデバイスの組み合わせとして実装可能であり、例えば、DSPおよびマイクロプロセッサの組み合わせ、複数のマイクロプロセッサ、DSPコアと組み合わせた1つ以上のマイクロプロセッサ、または任意の他のかかる構成として実装可能である。
付加的に、本明細書に開示される実施形態に関連して説明される方法またはアルゴリズムのステップは、ハードウェア、プロセッサによって実行されるソフトウェアモジュール、または両者の組み合わせで直接に具現化することができる。ソフトウェアモジュールは、RAMメモリ、フラッシュメモリ、ROMメモリ、EPROMメモリ、EEPROMメモリ、レジスタ、ハードディスク、リムーバブルディスク、CD-ROM、またはネットワーク記憶媒体を含むその他の形態の記憶媒体に常駐しうる。例示的な記憶媒体は、プロセッサが記憶媒体から情報を読み取り、記憶媒体に情報を書き込むことができるように、プロセッサに結合することができる。別の方法として、記憶媒体をプロセッサに統合することもできる。プロセッサおよび記憶媒体はASICに常駐することもできる。
本明細書で言及されるコンピュータ可読記憶媒体は、それ自体が、電波もしくは他の自由に伝搬する電磁波、導波管もしくは他の伝送媒体を伝搬する電磁波(例えば、光ファイバケーブルを通過する光パルス)、またはワイヤを介して伝送される電気信号などの一時的な信号であると解釈されるべきではない。
本明細書に記載されているソフトウェアコンポーネントのいずれも、多種多様な形態をとりうる。例えば、コンポーネントは、スタンドアロンのソフトウェアパッケージであってもよく、またはより大きなソフトウェア製品に「ツール」として組み込まれているソフトウェアパッケージであってもよい。これは、Webサイトなどのネットワークからスタンドアロン製品として、または既存のソフトウェアアプリケーションにインストールするためのアドインパッケージとしてダウンロード可能であってもよい。また、クライアントサーバソフトウェアアプリケーション、Web対応ソフトウェアアプリケーション、および/またはモバイルアプリケーションとして利用されてもよい。
本発明の実施形態は、本発明の実施形態による方法、装置(システム)、およびコンピュータプログラム製品のフローチャート図および/またはブロック図を参照して本明細書に記載されている。フローチャート図および/またはブロック図の各ブロック、ならびにフローチャート図および/またはブロック図のブロックの組み合わせは、コンピュータ可読プログラム命令によって実装できることが理解されよう。
コンピュータ可読プログラム命令は、汎用コンピュータ、専用コンピュータ、または他のプログラム可能なデータ処理装置のプロセッサに提供され、コンピュータまたは他のプログラム可能なデータ処理装置のプロセッサを介して実行される命令が、フローチャートおよび/またはブロック図の単一もしくは複数のブロックで指定された機能/動作を実施するための手段を作成するような機械を作り出すことができる。これらのコンピュータ可読プログラム命令はまた、コンピュータ、プログラム可能なデータ処理装置、および/または他のデバイスに特定の方法で機能するように指示しうるコンピュータ可読記憶媒体に保存可能であり、命令が保存されたコンピュータ可読記憶媒体には、フローチャートおよび/またはブロック図の単一もしくは複数のブロックで指定された機能/動作の態様を実施する命令を含む製造品が含まれる。
コンピュータ可読プログラム命令はまた、コンピュータ、他のプログラム可能なデータ処理装置、または他のデバイスにロードされて、コンピュータ、他のプログラム可能な装置または他のデバイス上で一連の操作ステップを実行させ、コンピュータ実装プロセスを生成可能であり、その結果、コンピュータ、他のプログラム可能な装置、または他のデバイス上で実行される命令により、フローチャートおよび/またはブロック図の単一もしくは複数のブロックで指定された機能/動作が実装される。
図示のフローチャートおよびブロック図は、本発明の種々の実施形態による、システム、方法、およびコンピュータプログラム製品の可能な実装形態のアーキテクチャ、機能、および動作を示している。これに関して、フローチャートまたはブロック図の各ブロックは、モジュール、セグメント、または命令の一部を表すことができ、これには、指定された論理機能を実装するための1つ以上の実行可能命令が含まれる。幾つかの代替の実装形態では、ブロックに示されている機能は、図に示されている順序とは異なる場合がある。例えば、連続して表示される2つのブロックは、実際には実質的に同時に実行される場合があり、または関連する機能に応じてブロックが逆の順序で実行される場合がある。ブロック図および/またはフローチャート図の各ブロック、およびブロック図および/またはフローチャート図のブロックの組み合わせは、指定された機能もしくは動作を実行するか、または専用ハードウェアおよびコンピュータ命令の組み合わせを実行する専用ハードウェアベースのシステムによって実装されうることにも留意されたい。
本発明を具現化する装置および方法は、クラウドコンピューティング環境でホストされ、クラウドコンピューティング環境によって配信されうる。クラウドコンピューティングは、構成可能なコンピューティングリソース(例えば、ネットワーク、ネットワーク帯域幅、サーバ、処理、メモリ、ストレージ、アプリケーション、仮想マシン、およびサービス)の共有プールへの便利なオンデマンドネットワークアクセスを可能にするサービス提供のモデルであり、最小限の管理作業またはサービスのプロバイダとのインタラクションにより、迅速にプロビジョニングおよびリリースされうる。このクラウドモデルには、少なくとも5つの特徴、少なくとも3つのサービスモデル、および少なくとも4つの展開モデルが含まれうる。
その特徴は以下の通りである。
オンデマンドセルフサービス:クラウドコンシューマは、人間のサービスプロバイダとのインタラクションを必要とせずに、必要に応じてサーバ時間およびネットワークストレージなどのコンピューティング機能を一方的にプロビジョニングできる。
ブロードネットワークアクセス:ネットワーク経由で機能が利用可能であり、異種のシンクライアントまたはシッククライアントプラットフォーム(携帯電話、ラップトップ、PDAなど)による使用を促進する標準メカニズムを介してアクセスされる。
リソースプーリング:プロバイダのコンピューティングリソースは、マルチテナントモデルを使用して複数のコンシューマにサービスを提供するためにプールされ、種々の物理リソースおよび仮想リソースが需要に応じて動的に割り当てられ、かつ再割り当てされる。コンシューマは通常、提供されたリソースの正確な場所を制御したり、知ったりすることはできないが、より高い抽象度で場所(国、州、データセンタなど)を指定する可能性があるという点で、場所の独立性がある。
迅速な弾力性(Rapid elasticity):迅速にかつ弾力的に、幾つかの場合には自動的に、機能がプロビジョニングされ、迅速にスケールアウトされ、迅速にリリースされて迅速にスケールインされる。コンシューマにとっては、プロビジョニングに利用できる機能は無制限で、いつでも好きなだけ購入できるように見えることが多い。
計測可能なサービス(Measured service):クラウドシステムにより、サービスのタイプ(ストレージ、処理、帯域幅、アクティブなユーザアカウントなど)に適した抽象化レベルで計測機能を活用することで、リソースの使用が自動的に制御および最適化される。リソースの使用状況を監視、制御、および報告できるため、利用するサービスのプロバイダおよびコンシューマの両方に透明性が提供される。
サービスモデルは以下の通りである。
サービスとしてのソフトウェア(SaaS):コンシューマに提供される機能は、クラウドインフラストラクチャで実行されているプロバイダのアプリケーションを使用することである。アプリケーションには、Webブラウザ(Webベースの電子メールなど)などのシンクライアントインタフェースを介して、種々のクライアントデバイスからアクセスできる。コンシューマは、限られたユーザ固有のアプリケーション構成設定を除いて、ネットワーク、サーバ、オペレーティングシステム、ストレージ、さらには個々のアプリケーション機能を含む基盤となるクラウドインフラストラクチャを管理または制御しない。
サービスとしてのプラットフォーム(PaaS):コンシューマに提供される機能は、プロバイダがサポートするプログラミング言語およびツールを使用して作成された、コンシューマが作成または取得したアプリケーションをクラウドインフラストラクチャにデプロイすることである。コンシューマは、ネットワーク、サーバ、オペレーティングシステム、ストレージなどの基盤となるクラウドインフラストラクチャを管理または制御しないが、デプロイされたアプリケーションと、幾つかの場合にはアプリケーションホスティング環境の構成とを制御する。
サービスとしてのインフラストラクチャ(IaaS):コンシューマに提供される機能は、処理、ストレージ、ネットワーク、およびその他の基本的なコンピューティングリソースをプロビジョニングすることであり、コンシューマはオペレーティングシステムおよびアプリケーションを含みうる任意のソフトウェアをデプロイし、かつ実行できる。コンシューマは、基盤となるクラウドインフラストラクチャを管理または制御しないが、オペレーティングシステム、ストレージ、デプロイされたアプリケーションを制御し、幾つかの場合には選択したネットワークコンポーネント(ホストファイアウォールなど)の制御を制限する。
展開モデルは以下の通りである。
プライベートクラウド:クラウドインフラストラクチャは、組織のためのみに運用される。組織またはサードパーティによって管理可能であり、オンプレミスまたはオフプレミスに存在しうる。
コミュニティクラウド:クラウドインフラストラクチャは複数の組織で共有されており、問題(ミッション、セキュリティ要件、ポリシー、コンプライアンスの考慮事項など)を共有している特定のコミュニティをサポートする。組織またはサードパーティによって管理可能であり、オンプレミスまたはオフプレミスに存在しうる。
パブリッククラウド:クラウドインフラストラクチャは、一般の人々または大規模な業界団体が利用できるようになっており、クラウドサービスを販売する組織が所有している。
ハイブリッドクラウド:クラウドインフラストラクチャは、2つ以上のクラウド(プライベート、コミュニティ、またはパブリック)の構成であり、当該クラウドは、一意のエンティティのままであるが、データおよびアプリケーションの移植性(例えば、クラウド間の負荷分散のためのクラウドバースト)を可能にする標準化されたテクノロジーまたは独自のテクノロジーによって結び付けられている。
クラウドコンピューティング環境は、サービス指向であり、ステートレス、低結合、モジュール性、およびセマンティック相互運用性に重点を置いている。クラウドコンピューティングの中心は、相互接続されたノードのネットワークを含むインフラストラクチャである。
本開示の範囲から逸脱することなく、前述の例示的な実施形態に対して多くの改善および修正を行うことができることは、当技術分野の当業者には明らかであろう。
図11Aは、本明細書に記載の種々の実施形態に関連して使用されうる例示的なプロセッサ対応デバイス551を示すブロック図である。デバイス551の代替形態もまた、当業者によって理解されるように使用されうる。図示の実施形態では、デバイス551は、デジタル撮像デバイス(本明細書では、スキャナシステムまたはスキャンシステムとも称される)として提示されており、当該デジタル撮像デバイスには、1つ以上のプロセッサ556、1つ以上のメモリ566、1つ以上のモーションコントローラ571、1つ以上のインタフェースシステム576、各々が1つ以上のサンプル590を備えた1つ以上のスライドガラス585を支持する1つ以上の可動ステージ580、サンプルを照明する1つ以上の照明システム595、各々が光軸に沿って移動する光路605を確定する1つ以上の対物レンズ600、1つ以上の対物レンズポジショナ630、(例えば、蛍光スキャナシステムに含まれる)1つ以上の任意の落射照明システム635、1つ以上の集束光学系610、1つ以上のラインスキャンカメラ615および/または1つ以上のエリアスキャンカメラ620が含まれ、各々がサンプル590および/またはスライドガラス585上に別個の視野625を定めている。スキャナシステム551の種々の要素は、1つ以上の通信バス560を介して通信可能に結合されている。スキャナシステム551の種々の各要素のうちの1つ以上が存在しうるが、以下の説明を容易にするために、これらの要素は、適切な情報を伝えるために複数形で説明する必要がある場合を除いて、単数形で説明する。
1つ以上のプロセッサ556は、例えば、命令を並列に処理しうる中央処理装置(「CPU」)および別個のグラフィックス処理装置(「GPU」)を含みうるか、あるいは1つ以上のプロセッサ556は、命令を並行して処理しうるマルチコアプロセッサを含みうる。特定のコンポーネントを制御したり、画像処理などの特定の機能を実行したりするために、追加の個別のプロセッサを提供することもできる。例えば、追加のプロセッサには、データ入力を管理するための補助プロセッサ、浮動小数点数学演算を実行するための補助プロセッサ、信号処理アルゴリズムの高速実行に適したアーキテクチャを有する専用プロセッサ(例えば、デジタル信号プロセッサ)、メインプロセッサに従属するスレーブプロセッサ(例えば、バックエンドプロセッサ)、ラインスキャンカメラ615、ステージ580、対物レンズ600、および/またはディスプレイ(図示せず)を制御するための追加のプロセッサが含まれてもよい。かかる追加のプロセッサは、別個の離散的なプロセッサであってもよいし、またはプロセッサ556と統合されてもよい。
メモリ566は、プロセッサ556によって実行されうるプログラムのためのデータおよび命令のストレージを提供する。メモリ566は、データおよび命令、例えば、ランダムアクセスメモリ、読み取り専用メモリ、ハードディスクドライブ、リムーバブルストレージドライブなどを保存する1つ以上の揮発性かつ永続性のコンピュータ可読記憶媒体を含みうる。プロセッサ556は、メモリ566に保存され、通信バス560を介してスキャナシステム551の種々の要素と通信して、スキャナシステム551の全体的な機能を実行する命令を実行するように構成される。
1つ以上の通信バス560は、アナログ電気信号を伝達するように構成された通信バス560を含むことができ、デジタルデータを伝達するように構成された通信バス560を含むことができる。したがって、1つ以上の通信バス560を介したプロセッサ556、モーションコントローラ571、および/またはインタフェースシステム576からの通信には、電気信号およびデジタルデータの両方が含まれうる。プロセッサ556、モーションコントローラ571、および/またはインタフェースシステム576はまた、無線通信リンクを介してスキャンシステム551の種々の要素のうちの1つ以上と通信するように構成されうる。
モーション制御システム571は、ステージ580および対物レンズ600のXYZ移動を(例えば、対物レンズポジショナ630を介して)正確に制御および調整するように構成される。モーション制御システム571はまた、スキャナシステム551内の他の可動部分の移動を制御するように構成される。例えば、蛍光スキャナの実施形態では、モーション制御システム571は、落射照明システム635内の光学フィルタなどの移動を調整するように構成される。
インタフェースシステム576により、スキャナシステム551は、他のシステムおよび人間のオペレータとのインタフェースを有することが可能となる。例えば、インタフェースシステム576は、情報をオペレータに直接に提供するため、かつ/またはオペレータからの直接の入力を可能にするためのユーザインタフェースを含みうる。インタフェースシステム576はまた、スキャンシステム551と、直接に接続されている1つ以上の外部デバイス(例えば、プリンタ、リムーバブル記憶媒体)または外部デバイス、例えばネットワーク(図示せず)を介してスキャナシステム551に接続されている画像サーバシステム、オペレータステーション、ユーザステーション、および管理サーバシステムとの間の通信およびデータ転送を容易にするように構成される。
照明システム595は、サンプル590の一部を照明するように構成される。照明システムには、例えば、光源および照明光学系が含まれてもよい。光源は、光出力を最大化するための凹面反射鏡と熱を抑制するためのKG-1フィルタとを備えた可変強度ハロゲン光源でありうる。光源はまた、任意のタイプのアークランプ、レーザ、または他の光源でありうる。一実施形態では、照明システム595は、ラインスキャンカメラ615および/またはエリアスキャンカメラ620がサンプル590を介して送信される光エネルギをセンシングするように、透過モードでサンプル590を照明する。代替的に、または組み合わせて、照明システム595はまた、反射モードでサンプル590を照明するように構成可能であり、これにより、ラインスキャンカメラ615および/またはエリアスキャンカメラ620は、サンプル590から反射される光エネルギをセンシングする。全体として、照明システム595は、任意の既知の光学顕微鏡モードでの顕微鏡サンプル590の調査に適しているように構成される。
一実施形態では、スキャナシステム551は、任意に落射照明システム635を含み、落射照明システム635は、スキャナシステム551を最適化して蛍光スキャンを行う。蛍光スキャンは、蛍光分子を含むサンプル590のスキャンであり、この蛍光分子は、特定の波長(励起)の光を吸収しうる光子感受性分子である。これらの光子感受性分子はまた、より高い波長(発光)で発光する。このフォトルミネッセンス現象の効率は非常に低いため、多くの場合、放出される光の量は非常に少ない。この放出光が少量であるため、従来の技術(例えば、透過モード顕微鏡)では、通常、サンプル590をスキャンおよびデジタル化することは困難である。有利には、スキャナシステム551の任意の蛍光スキャナシステムの実施形態では、複数の線形センサアレイ(例えば、時間遅延積分(「TDI」)ラインスキャンカメラ)を含むラインスキャンカメラ615を使用して、サンプル590の同じ領域をラインスキャンカメラ615の複数の線形センサアレイの各々に露光することによって、ラインスキャンカメラの光に対する感度が高められる。これは、放出光が少ない微弱な蛍光サンプルを走査する場合に特に便利である。
したがって、蛍光スキャナシステムの実施形態では、ラインスキャンカメラ615は、好ましくは、モノクロTDIラインスキャンカメラである。有利には、モノクロ画像は、サンプル上に存在する種々のチャネルからの実際の信号のより正確な表現を提供するため、蛍光顕微鏡に理想的な画像である。当業者によって理解されるように、蛍光サンプル590は、異なる波長で光を放出する複数の蛍光色素で標識することができ、これは「チャネル」とも称される。
さらに、種々の蛍光サンプルのローエンドおよびハイエンドの信号レベルは、ラインスキャンカメラ615がセンシングしうる広範なスペクトルの波長を提示するので、ラインスキャンカメラ615がセンシングしうるローエンドおよびハイエンドの信号レベルは、同様に広範であることが望ましい。したがって、蛍光スキャナの実施形態では、蛍光スキャンシステム551で使用されるラインスキャンカメラ615は、モノクロの10ビット64線形アレイTDIラインスキャンカメラである。スキャンシステム551の蛍光スキャナの実施形態で使用するために、ラインスキャンカメラ615の多種多様なビット深度が利用されうることに留意されたい。
可動ステージ580は、プロセッサ556またはモーションコントローラ571の制御下の正確なXY移動のために構成される。可動ステージはまた、プロセッサ556またはモーションコントローラ571の制御下でZ方向内で移動するように構成されうる。可動ステージは、ラインスキャンカメラ615および/またはエリアスキャンカメラによる画像データの取り込み中に、所望の位置にサンプルを位置決めするように構成される。可動ステージはまた、サンプル590を走査方向に実質的に一定の速度に加速し、ついで、ラインスキャンカメラ615による画像データの取り込み中に実質的に一定の速度を維持するように構成される。一実施形態では、スキャナシステム551は、可動ステージ580上のサンプル590の位置を補助するために、高精度で密に調整されたXYグリッドを使用することができる。一実施形態では、可動ステージ580は、高精度エンコーダがX軸およびY軸の両方に採用された、リニアモータベースのXYステージである。例えば、非常に精密なナノメートルエンコーダを、走査方向の軸上で使用することができ、かつ走査方向に垂直な方向にあり、走査方向と同じ平面上にある軸上で使用することができる。ステージはまた、サンプル590が配設されたスライドガラス585を支持するように構成される。
サンプル590は、光学顕微鏡によって調査されうるものであればいかなるものであってもよい。例えば、顕微鏡スライドガラス585は、標本の観察基板として頻繁に使用されており、当該観察基板としては、組織および細胞、染色体、DNA、タンパク質、血液、骨髄、尿、細菌、ビーズ、生検材料、または死んでいるもしくは生きている、染色もしくは非染色の、標識を有するもしくは非標識であるその他の種類の生物学的材料または物質が挙げられる。サンプル590はまた、任意のタイプのスライドまたは他の基板上に堆積される、任意のタイプのDNAまたはcDNAもしくはRNAもしくはタンパク質などのDNA関連材料のアレイであってよく、これには、マイクロアレイとして一般に既知である全てのサンプルが含まれる。サンプル590は、マイクロタイタープレート、例えば、96ウェルプレートでありうる。サンプル590の他の例として、集積回路基板、電気泳動記録、ペトリ皿、フィルム、半導体材料、法医学材料、または機械加工部品が挙げられる。
対物レンズ600は、一実施形態では、対物レンズポジショナ630に取り付けられており、対物レンズポジショナ630は、非常に精密なリニアモータを使用して、対物レンズ600によって定められる光軸に沿って対物レンズ600を移動させうる。例えば、対物レンズポジショナ630のリニアモータは、50ナノメートルのエンコーダを含んでもよい。XYZ軸におけるステージ580および対物レンズ600の相対位置は、情報および命令を保存するためにメモリ566を使用するプロセッサ556の制御下でモーションコントローラ571を使用して、閉ループ方式で調整および制御される。ここで、当該情報および命令には、全体的なスキャンシステム551の動作のためのコンピュータ実行可能なプログラムされたステップが含まれる。
一実施形態では、対物レンズ600は、望ましい最高の空間分解能に対応する開口数を備えた平面アポクロマティック(「APO」)無限遠補正対物レンズであり、ここで、対物レンズ600は、透過モード照明顕微鏡、反射モード照明顕微鏡、および/または落射照明モード蛍光顕微鏡(例えば、オリンパス社40X、0.75NAまたは20X、0.75NA)に適している。有利には、対物レンズ600は、色収差および球面収差を補正することができる。対物レンズ600は無限遠補正されているので、集束光学系610は、対物レンズ600の上の光路605に配置することができ、当該光路において、対物レンズを通過する光ビームは、コリメートされた光ビームとなる。集束光学系610は、対物レンズ600によって捕捉された光信号を、ラインスキャンカメラ615および/またはエリアスキャンカメラ620の光応答要素に集束させる。また、集束光学系610は、フィルタ、倍率変更レンズなどの光学部品を含みうる。集束光学系610と組み合わされた対物レンズ600により、スキャンシステム551の総合倍率が提供される。一実施形態では、集束光学系610は、管レンズおよび任意の2倍倍率チェンジャを含みうる。有利には、2倍倍率チェンジャにより、ネイティブ20倍の対物レンズ600は、40倍の倍率でサンプル590を走査することが可能となる。
ラインスキャンカメラ615は、画像要素(「ピクセル」)の少なくとも1つの線形アレイを含む。ラインスキャンカメラは、モノクロカメラまたはカラーカメラでありうる。カラーラインスキャンカメラは通常、少なくとも3つの線形アレイを有するが、モノクロラインスキャンカメラは、単一の線形アレイまたは複数の線形アレイを有しうる。また、任意のタイプの単数または複数の線形アレイは、カメラの一部としてパッケージ化されているか、撮像電子モジュールにカスタム統合されているかにかかわらず使用できる。例えば、3線形アレイ(「赤色-緑色-青色」または「RGB」)カラーラインスキャンカメラまたは96線形アレイモノクロTDIも使用できる。TDIラインスキャンカメラは、通常、過去に撮像された標本の領域の強度データを合計することにより、実質的により優れた出力信号の信号対雑音比(「SNR」)をもたらす。ここで、SNRの増大は、積分段数の平方根に比例する。TDIラインスキャンカメラは複数の線形アレイで構成され、例えば、TDIラインスキャンカメラは24、32、48、64、96個、またはそれ以上の線形アレイで利用できる。スキャナシステム551はまた、512ピクセルのもの、1024ピクセルのもの、および4096ピクセルものものを含む多種多様なフォーマットで製造される線形アレイをサポートしている。同様に、多種多様なピクセルサイズの線形アレイもスキャナシステム551で使用することができる。任意のタイプのラインスキャンカメラ615を選択するための重要な要件は、ステージ580のモーションが、ラインスキャンカメラ615のラインレートと同期されうることにより、サンプル590のデジタル画像の取り込み中に、ステージ580がラインスキャンカメラ615に対して移動しうることである。
ラインスキャンカメラ615によって生成された画像データは、メモリ566の一部に保存され、プロセッサ556によって処理されて、サンプル590の少なくとも一部の連続デジタル画像が生成される。連続デジタル画像は、プロセッサ556によってさらに処理することができ、修正された連続デジタル画像はまた、メモリ566に保存することができる。
2つ以上のラインスキャンカメラ615を有する一実施形態では、ラインスキャンカメラ615の少なくとも1つは集束センサとして機能するように構成可能であり、当該集束センサは、イメージセンサとして機能するように構成されたラインスキャンカメラ615のうちの少なくとも1つと組み合わされて動作する。集束センサは、イメージセンサと同じ光軸上に論理的に位置決め可能であり、あるいは集束センサは、スキャナシステム551の走査方向に関してイメージセンサの前または後に論理的に位置決め可能である。かかる集束センサとして機能する少なくとも1つのラインスキャンカメラ615を有する一実施形態では、集束センサによって生成された画像データは、メモリ566の一部に保存され、1つ以上のプロセッサ556によって処理されて、焦点情報を生成し、スキャナシステム551は、サンプル590と対物レンズ600との間の相対距離を調整して、走査中にサンプルに焦点を合わせ続ける。さらに、一実施形態では、集束センサとして機能する少なくとも1つのラインスキャンカメラ615は、集束センサの複数の個々のピクセルの各々が、光路605に沿って異なる論理高さに位置決めされるように配向されうる。
動作中、スキャナシステム551の種々の構成要素およびメモリ566に保存されたプログラムされたモジュールにより、スライドガラス585上に配設されたサンプル590の自動走査およびデジタル化が可能となる。スライドガラス585は、サンプル590を走査するために、スキャナシステム551の可動ステージ580上に堅固に配置されている。プロセッサ556の制御下で、可動ステージ580は、ラインスキャンカメラ615によるセンシングのためにサンプル590を実質的に一定の速度に加速する。ここで、ステージの速度は、ラインスキャンカメラ615のライン速度と同期されている。画像データのストリップを走査した後、可動ステージ580は減速して、サンプル590を実質的に完全に停止させる。次に、可動ステージ580は、走査方向に直交して移動して、画像データの後続のストリップ、例えば、隣接するストリップを走査するために、サンプル590を位置決めする。続いて、サンプル590の全体またはサンプル590全体が走査されるまで、追加のストリップが走査される。
例えば、サンプル590のデジタル走査中、複数の連続する視野としてサンプル590の連続デジタル画像が取得され、これらは一緒に組み合わされて画像ストリップが形成される。複数の隣接する画像ストリップが同様に一緒に組み合わされて、サンプル590の一部または全体の連続したデジタル画像が形成される。サンプル590のスキャンには、垂直画像ストリップまたは水平画像ストリップの取得が含まれうる。サンプル590のスキャン方向は、上から下、下から上、またはその両方(双方向)のいずれかであってよく、サンプル上の任意の点で開始されうる。代替的に、サンプル590のスキャン方向は、左から右、右から左、またはその両方(双方向)のいずれかであってもよく、サンプル上の任意の点で開始されうる。さらに、画像ストリップを隣接または連続して取得する必要はない。さらに、結果として得られるサンプル590の画像は、サンプル590全体の画像であってもよく、サンプル590の一部のみの画像であってもよい。
一実施形態では、コンピュータ実行可能命令(例えば、プログラムされたモジュールおよびソフトウェア)は、メモリ566に保存され、実行される際に、スキャンシステム551は、本明細書に記載の種々の機能を実施することができる。この説明では、「コンピュータ可読記憶媒体」という用語は、プロセッサ556によって実行するために、コンピュータ実行可能命令を保存し、スキャンシステム551に提供するために使用される任意の媒体を指すために使用される。これらの媒体の例には、メモリ566、および例えばネットワーク(図示せず)を介して直接的にまたは間接的にスキャンシステム551と通信可能に結合された任意のリムーバブル媒体または外部記憶媒体(図示せず)が含まれる。
図11Bは、電荷結合デバイス(「CCD」)アレイとして実装されうる単一の線形アレイ640を有するラインスキャンカメラを示している。単一の線形アレイ640は、複数の個別のピクセル645を含む。図示の実施形態では、単一の線形アレイ640は、4096ピクセルを有する。代替的な実施形態では、線形アレイ640は、より多くのピクセルまたはより少ないピクセルを有しうる。例えば、線形アレイの一般的なフォーマットには、512、1024、および4096ピクセルが含まれる。ピクセル645は、線形アレイ640の視野625を定めるために線形に配置されている。視野の大きさは、スキャナシステム551の倍率に応じて変化する。
図11Cは、3つの線形アレイを有するラインスキャンカメラを示しており、その各々は、CCDアレイとして実装されうる。3つの線形アレイが組み合わさって、カラーアレイ650を形成している。一実施形態では、カラーアレイ650内の個々の線形アレイは、様々な色強度、例えば、赤色、緑色、または青色を検出する。カラーアレイ650内の個々の線形アレイからのカラー画像データは、組み合わされて、カラー画像データの単一の視野625を形成している。
図11Dは、複数の線形アレイを有するラインスキャンカメラを示しており、その各々は、CCDアレイとして実装されうる。複数の線形アレイが組み合わさって、TDIアレイ655を形成している。有利には、TDIラインスキャンカメラは、過去に撮像された標本の領域の強度データを合計することにより、実質的により優れた出力信号のSNRをもたらしうる。ここで、SNRの増大は、(積分段とも称される)線形アレイの数の平方根に比例する。TDIラインスキャンカメラは、より多くの多種多様な線形アレイを含みうる。例えば、TDIラインスキャンカメラの一般的なフォーマットには、24、32、48、64、96、120個、およびさらに多くの線形アレイが含まれる。
開示した実施形態の上記の説明は、当業者が本発明を作製または使用することを可能にするために提供されている。これらの実施形態に対する種々の修正例は、当業者には容易に明らかとなり、本明細書に記載の全般的な基本方式は、本発明の精神または範囲から逸脱することなく、他の実施形態に適用することができる。したがって、本明細書に提示される説明および図面は、本発明の現在の好ましい実施形態を表し、したがって、本発明によって広く企図される主題を代表するものであることが理解されるべきである。さらに、本発明の範囲は、当技術分野の当業者に明らかになりうる他の実施形態を完全に包含しており、これによって本発明の範囲は限定されないことが理解されよう。

Claims (71)

  1. 組織学的画像内の腫瘍を同定する方法であって、前記方法が、
    ピクセルの2次元アレイを含む組織学的画像を受信することと、
    畳み込みニューラルネットワークを適用して、組織学的画像の2次元アレイに対するマッピングを有するピクセルの2次元アレイを有する出力画像を生成することであって、前記出力画像が、複数の組織クラスのうちの1つを各ピクセルに割り当てることによって生成され、前記複数の組織クラスは、非腫瘍組織を表す少なくとも1つのクラスおよび腫瘍組織を表す少なくとも1つのクラスを含む、ことと、
    個々の腫瘍によって占有された領域がマーキングされている出力画像から、セグメンテーションマスクを生成することと、
    各腫瘍の要約統計量を計算することと、
    前記各腫瘍の要約統計量にフィルタを適用して、前記フィルタに従って腫瘍を選択および選択解除することにより、前記セグメンテーションマスクを編集することと、
    を含む、方法。
  2. 前記方法が、
    スコアリングアルゴリズムに従って各腫瘍をスコアリングして、各腫瘍にスコアを割り当てることと、
    そのスコアに従って前記腫瘍をランク付けすることと、
    をさらに含む、
    請求項1記載の方法。
  3. 前記ランク付けは、前記フィルタを適用した後に残った腫瘍に限定される、
    請求項2記載の方法。
  4. 前記方法が、前記ランク付けを保存することをさらに含む、
    請求項2または3記載の方法。
  5. 前記方法は、編集された前記セグメンテーションマスクを考慮して前記組織学的画像のビジュアライゼーションを作成することをさらに含む、
    請求項1、2、3または4記載の方法。
  6. 前記ビジュアライゼーションは、前記セグメンテーションマスクが前記組織学的画像上にオーバレイされる概観表示ペインを含む、
    請求項5記載の方法。
  7. 前記概観表示ペインは、各腫瘍のランク付けラベルを含む、
    請求項2に付加される請求項6記載の方法。
  8. 前記ビジュアライゼーションは、前記セグメンテーションマスクおよび前記組織学的画像が1対1の比較のために互いに隣接して提示されるそれぞれの概観表示ペインを含む、
    請求項5記載の方法。
  9. 前記セグメンテーションマスクの前記概観表示ペインは、各腫瘍のランク付けラベルを含む、
    請求項2に付加される請求項8記載の方法。
  10. 前記方法が、編集された前記セグメンテーションマスク内の腫瘍を選択するためにユーザが前記ビジュアライゼーションとのインタラクションを行うことを可能にするように動作可能なユーザインタフェース腫瘍選択制御をさらに含む、
    請求項5から9までのいずれか1項記載の方法。
  11. 前記ビジュアライゼーションは、現在選択されている腫瘍にズームインされた拡大表示ペインを含む、
    請求項10記載の方法。
  12. 前記腫瘍選択制御は、編集された前記セグメンテーションマスク内の前記腫瘍をランク付け順でスイープするためのスクロール機能を有する、
    請求項10または11記載の方法。
  13. 前記腫瘍選択制御は、ユーザが前記現在選択されている腫瘍に対して追加の計算診断プロセスを実行することを可能にするための診断機能を有する、
    請求項10、11または12記載の方法。
  14. 前記組織クラスが、腫瘍組織について少なくとも2つのクラスを含む、
    請求項1から13までのいずれか1項記載の方法。
  15. 腫瘍組織の前記組織クラスが、少なくとも浸潤性腫瘍の第1のクラスおよび非浸潤性腫瘍の第2のクラスを含む、
    請求項14記載の方法。
  16. 非腫瘍組織に対して1つの組織クラスが存在している、
    請求項1から15までのいずれか1項記載の方法。
  17. 前記方法が、前記ビジュアライゼーションをディスプレイ上に表示することをさらに含む、
    請求項1から16までのいずれか1項記載の方法。
  18. 前記畳み込みニューラルネットワークを適用することが、
    前記組織学的画像から画像パッチを抽出することであって、前記画像パッチが、幅および高さのピクセル数によって定められるサイズを有する前記組織学的画像の領域部分である、ことと、
    前記畳み込みニューラルネットワークに重みのセットおよび複数のチャネルを提供することであって、各チャネルが、特定される複数の前記組織クラスのうちの1つに対応する、ことと、
    各画像パッチを入力画像パッチとして前記畳み込みニューラルネットワークに入力することと、
    多段畳み込みを実行して、最小次元の最終畳み込み層までさらに次元削減された畳み込み層であって、最小次元の最終畳み込み層を含む、畳み込み層を生成し、続いて、前記入力画像パッチとサイズが一致する層が復元されるまで、さらに次元増大された逆畳み込み層を生成することにより、前記畳み込みを反転する多段転置畳み込みを実行することであって、前記復元された層の各ピクセルが、前記組織クラスの各々に属する確率を含む、ことと、
    出力画像パッチに到達する前記確率に基づいて、前記復元された層の各ピクセルに前記組織クラスを割り当てることと、
    を含む、
    請求項1から17までのいずれか1項記載の方法。
  19. 前記方法が、前記出力画像パッチを前記組織学的画像の確率マップにアセンブルすることをさらに含む、
    請求項18記載の方法。
  20. 前記方法が、前記確率マップに従って腫瘍に対応する前記組織学的画像内の領域を画定することをさらに含む、
    請求項18または19記載の方法。
  21. 前記方法が、前記組織学的画像にリンクされた前記確率マップを保存することをさらに含む、
    請求項18記載の方法。
  22. 前記方法が、
    前記畳み込みニューラルネットワークに少なくとも1つのスキップ接続を提供することであって、前記スキップ接続の各々が、前記畳み込み層のうちの前記最終畳み込み層よりも大きな次元の少なくとも1つの畳み込み層から中間結果を取得し、前記結果を、なし、1つ、または2つ以上でありうる必要な数の転置畳み込みに供して、前記入力画像パッチにサイズが一致する少なくとも1つのさらなる復元された層を取得する、ことと、
    各ピクセルに組織クラスを割り当てる前記ステップの前に、前記復元された層をさらに処理して、前記少なくとも1つのスキップ接続を考慮して前記確率を再計算するために、前記少なくとも1つのさらなる復元された層と組み合わせることと、
    をさらに含む、
    請求項18記載の方法。
  23. ソフトマックス演算を使用して前記確率を生成する、
    請求項18記載の方法。
  24. 前記方法が、予測のために実行され、前記畳み込みニューラルネットワークに、事前の訓練中、重み値が割り当てられている、
    請求項18記載の方法。
  25. 前記方法が、訓練のために実行され、前記組織学的画像内の各ピクセルを前記組織クラスのうちの1つに割り当てるグラウンドトゥルースデータを提供することをさらに含み、
    請求項18記載の方法が反復的に実行され、ここで、各反復が、前記グラウンドトゥルースデータと前記出力画像パッチとの比較に基づいて前記畳み込みニューラルネットワークの重み値を調整することを含む、
    請求項18記載の方法。
  26. 訓練中に前記重みを調整することが、勾配降下法によって実行される、
    請求項25記載の方法。
  27. 前記組織学的画像は、組織の領域の異なって染色された隣接するセクションから取得された複数の組織学的画像を含む合成物である、
    請求項1から26までのいずれか1項記載の方法。
  28. 組織学的画像内の腫瘍を同定するためのコンピュータプログラム製品であって、請求項1から27までのいずれか1項記載の方法を実行するための機械可読命令を有する、コンピュータプログラム製品。
  29. 組織学的画像内の腫瘍を同定するためのコンピュータ装置であって、前記装置が、
    ピクセルの2次元アレイを含む組織学的画像を受信するように動作可能な入力部と、
    処理モジュールであって、
    畳み込みニューラルネットワークを適用して、組織学的画像の2次元アレイに対するマッピングを有するピクセルの2次元アレイを有する出力画像を生成することであって、前記出力画像が、複数の組織クラスのうちの1つを各ピクセルに割り当てることによって生成され、前記複数の組織クラスが、非腫瘍組織を表す少なくとも1つのクラスおよび腫瘍組織を表す少なくとも1つのクラスを含む、ことと、
    個々の腫瘍によって占有された領域がマーキングされている出力画像から、セグメンテーションマスクを生成することと、
    各腫瘍の要約統計量を計算することと、
    前記各腫瘍の要約統計量にフィルタを適用して、前記フィルタに従って腫瘍を選択および選択解除することにより、前記セグメンテーションマスクを編集することと、
    を行うように構成された、処理モジュールと、
    を含む、装置。
  30. 前記処理モジュールが、前記組織学的画像および編集された前記セグメンテーションマスクのビジュアライゼーションを作成するようにさらに構成されている、
    請求項29記載の装置。
  31. 前記装置が、前記セグメンテーションマスクを前記組織学的画像とともにメタデータとして保存するように動作可能な出力部をさらに含む、
    請求項29または30記載の装置。
  32. 前記装置が、前記要約統計量を前記組織学的画像とともにメタデータとして保存するように動作可能な出力部をさらに含む、
    請求項29、30、または31記載の装置。
  33. 前記装置が、
    ディスプレイと、
    前記組織学的画像が編集された前記セグメンテーションマスクで表示されるように、前記組織学的画像および前記セグメンテーションマスクを前記ディスプレイに送信するように動作可能なディスプレイ出力部と、
    をさらに含む、
    請求項29、30、31または32記載の装置。
  34. システムであって、
    請求項29から33までのいずれか1項記載のコンピュータ装置と、
    組織学的画像を取得するように動作可能な画像取得装置と、
    を含む、システム。
  35. システムであって、
    請求項29から33までのいずれか1項記載のコンピュータ装置と、
    組織学的画像を含む患者データのレコードを保存するように構成されたデータリポジトリと、
    前記コンピュータ装置と前記データリポジトリとの間の患者データレコードまたはその一部の転送を可能にするネットワーク接続と、
    を含む、システム。
  36. システムであって、
    請求項29から33までのいずれか1項記載のコンピュータ装置と、
    組織学的画像を取得するように動作可能な画像取得装置と、
    組織学的画像を含む患者データのレコードを保存するように構成されたデータリポジトリと、
    前記コンピュータ装置と前記データリポジトリとの間の患者データレコードまたはその一部の転送を可能にするネットワーク接続と、
    を含む、システム。
  37. 腫瘍塊の辺縁状態を求めるために組織学的画像を画像処理する方法であって、前記方法が、
    切除によって抽出され、その表面を強調表示するために染料で着色された腫瘍塊から得られた組織サンプル切片の組織学的画像を受信することであって、前記組織学的画像が、ピクセルの2次元アレイを含む、ことと、
    組織クラスを区別するように訓練された少なくとも1つのニューラルネットワークを使用した畳み込みニューラルネットワークプロセスを適用して、前記組織学的画像の2次元アレイに対するマッピングを有するピクセルの2次元アレイを有する出力画像を生成することであって、前記出力画像が、組織クラスを各ピクセルに割り当てて、染料および非染料ならびに腫瘍および非腫瘍を区別することによって生成される、ことと、
    腫瘍染料ピクセルの有無に応じて、辺縁状態を陽性または陰性としてそれぞれ計算することと、
    前記辺縁状態を出力することと、
    を含む、方法。
  38. 前記畳み込みニューラルネットワークプロセスが、腫瘍染料、腫瘍非染料、非腫瘍染料、および非腫瘍非染料の組織クラスを区別するように訓練された1つのニューラルネットワークを適用することを含む、
    請求項37記載の方法。
  39. 前記畳み込みニューラルネットワークプロセスが、腫瘍および非腫瘍の組織クラスを区別し、第1の出力画像を出力するように訓練された第1のニューラルネットワークと、染料および非染料の組織クラスを区別して第2の出力画像を出力するように訓練された第2のニューラルネットワークとを別々に適用することを含み、前記出力画像は、前記第1の出力画像および前記第2の出力画像を組み合わせて、前記組織クラスを腫瘍染料、腫瘍非染料、非腫瘍染料、および非腫瘍非染料の組織クラスに割り当てることによって生成される、
    請求項37記載の方法。
  40. 複数の腫瘍組織型を同定するためのクラスが存在している、
    請求項37から39までのいずれか1項記載の方法。
  41. 前記複数の腫瘍組織型は、浸潤性腫瘍および非浸潤性腫瘍を含む、
    請求項40記載の方法。
  42. 辺縁状態の前記計算により、陰性状態がクロースサブ状態とクリアサブ状態との間で細分化され、腫瘍非染料ピクセルと染料非腫瘍ピクセルとの間の距離が閾値を下回っているかまたは上回っていることによってクロースとクリアとが区別される、
    請求項37から41までのいずれか1項記載の方法。
  43. 辺縁状態の前記計算により、腫瘍染料ピクセルの有病率に基づいて複数のサブ状態間で前記陽性状態が細分化される、
    請求項1から42までのいずれか1項記載の方法。
  44. 有病率が、腫瘍染料ピクセル間の最大距離を計算することに少なくとも部分的に基づいて定量化される、
    請求項43記載の方法。
  45. 前記陽性または陰性の辺縁状態が、腫瘍染料ピクセルの数が閾値数を上回っているか、または下回っているかに少なくとも部分的に基づいて計算される、
    請求項1から44までのいずれか1項記載の方法。
  46. 前記畳み込みニューラルネットワークまたは各畳み込みニューラルネットワークを適用することが、
    前記組織学的画像から画像パッチを抽出することであって、前記画像パッチが、幅および高さのピクセル数によって定められるサイズを有する前記組織学的画像の領域部分である、ことと、
    重みのセットおよび複数のチャネルを提供することであって、各チャネルが、特定される複数の組織クラスのうちの1つに対応する、ことと、
    各画像パッチを入力画像パッチとして前記畳み込みニューラルネットワークに入力することと、
    多段畳み込みを実行して、最小次元の最終畳み込み層までさらに次元削減された畳み込み層であって、最小次元の最終畳み込み層を含む、畳み込み層を生成し、続いて、前記入力画像パッチとサイズが一致する層が復元されるまで、さらに次元増大された逆畳み込み層を生成することにより、前記畳み込みを反転する多段転置畳み込みを実行することであって、前記復元された層の各ピクセルが、前記組織クラスの各々に属する確率を含む、ことと、
    出力画像パッチに到達する前記確率に基づいて、前記復元された層の各ピクセルに前記組織クラスを割り当てることと、
    を含む、
    請求項37から45までのいずれか1項記載の方法。
  47. 前記組織学的画像が、組織の領域の異なって染色された隣接するセクションから取得された組織学的画像のセットである、
    請求項46記載の方法。
  48. 前記方法が、前記出力画像パッチを前記組織学的画像の確率マップにアセンブルすることをさらに含む、
    請求項46または47記載の方法。
  49. 前記方法が、前記組織学的画像にリンクされた前記確率マップを保存することをさらに含む、
    請求項46、47、または48記載の方法。
  50. 前記方法が、前記確率マップを用いて前記組織学的画像をディスプレイ上に表示することをさらに含む、
    請求項46、47、48、または49記載の方法。
  51. 前記方法が、
    前記畳み込みニューラルネットワークまたは各畳み込みニューラルネットワークに少なくとも1つのスキップ接続を提供することであって、前記スキップ接続の各々が、前記畳み込み層のうちの前記最終畳み込み層よりも大きな次元の少なくとも1つの畳み込み層から中間結果を取得し、前記結果を、なし、1つ、または2つ以上でありうる必要な数の転置畳み込みに供して、前記入力画像パッチにサイズが一致する少なくとも1つのさらなる復元された層を取得する、ことと、
    各ピクセルに組織クラスを割り当てる前記ステップの前に、前記復元された層をさらに処理して、前記少なくとも1つのスキップ接続を考慮して前記確率を再計算するために、前記少なくとも1つのさらなる復元された層と組み合わせることと、
    をさらに含む、
    請求項46から50までのいずれか1項記載の方法。
  52. ソフトマックス演算を使用して前記確率を生成する、
    請求項46から51までのいずれか1項記載の方法。
  53. 前記方法が、訓練のために実行され、前記組織学的画像内の各ピクセルを前記組織クラスのうちの1つに割り当てるグラウンドトゥルースデータを提供することをさらに含み、
    請求項46から52までのいずれか1項記載の方法が反復的に実行され、ここで、各反復が、前記グラウンドトゥルースデータと前記出力画像パッチとの比較に基づいて前記畳み込みニューラルネットワークの重み値を調整することを含む、
    請求項46から52までのいずれか1項記載の方法。
  54. 訓練中に前記重みを調整することが、勾配降下法によって実行される、
    請求項53記載の方法。
  55. 組織学的画像内の腫瘍を同定するためのコンピュータプログラム製品であって、
    請求項37から54までのいずれか1項記載の方法を実行するための機械可読命令を有する、コンピュータプログラム製品。
  56. 腫瘍塊の辺縁状態を求めるために組織学的画像を画像処理するためのコンピュータ装置であって、前記装置が、
    切除によって抽出され、その表面を強調表示するために染料で着色された腫瘍塊から得られた組織サンプル切片の組織学的画像を受信するように動作可能な入力部であって、前記組織学的画像が、ピクセルの2次元アレイを含む、入力部と、
    組織クラスを区別し、前記組織学的画像の2次元アレイに対するマッピングを有するピクセルの2次元アレイを有する出力画像を生成するように訓練された少なくとも1つのニューラルネットワークを含む畳み込みニューラルネットワーク処理モジュールであって、前記出力画像が、組織クラスを各ピクセルに割り当てて、染料および非染料ならびに腫瘍および非腫瘍を区別することによって生成される、畳み込みニューラルネットワーク処理モジュールと、
    前記出力画像内の腫瘍染料ピクセルの有無に応じて、辺縁状態を陽性または陰性としてそれぞれ計算するように動作可能な辺縁評価モジュールと、
    前記辺縁状態を出力するように動作可能な出力部と、
    を含む、装置。
  57. 前記出力部が、データリポジトリに接続されて、前記辺縁状態を前記組織学的画像とともにメタデータとして保存する、
    請求項56記載の装置。
  58. 前記出力部が、データリポジトリに接続されて、前記出力画像を前記組織学的画像とともにメタデータとして保存する、
    請求項56または57記載の装置。
  59. 前記出力画像が、各ピクセルの組織クラスを含む、
    請求項56、57、または58記載の装置。
  60. 前記出力画像が、各ピクセルの確率を含み、これにより確率マップを構成する、
    請求項58記載の装置。
  61. 前記出力部は、前記組織学的画像が前記出力画像とともに表示されるように、前記組織学的画像および前記出力画像をディスプレイに送信するように動作可能である、
    請求項56から60までのいずれか1項記載の装置。
  62. システムであって、
    請求項56から61までのいずれか1項記載のコンピュータ装置と、
    組織学的画像またはそのセットを含む患者データのレコードを保存するように構成されたデータリポジトリと、
    前記コンピュータ装置と前記データリポジトリとの間の患者データレコードまたはその一部の転送を可能にするネットワーク接続と、
    を含む、システム。
  63. システムであって、
    請求項56から61までのいずれか1項記載のコンピュータ装置と、
    組織学的画像またはそのセットを取得し、前記組織学的画像またはそのセットを前記データリポジトリ内のレコードに保存するように動作可能な画像取得装置と、
    を含む、システム。
  64. 前記システムが、
    組織学的画像またはそのセットを含む患者データのレコードを保存するように構成されたデータリポジトリと、
    前記コンピュータ装置と前記データリポジトリとの間の患者データレコードまたはその一部の転送を可能にするネットワーク接続と、
    をさらに含む、
    請求項63記載のシステム。
  65. 組織学的画像内の腫瘍を同定する方法であって、前記方法が、
    ピクセルの2次元アレイを含む組織学的画像を受信することと、
    畳み込みニューラルネットワークを適用して、ピクセルの2次元アレイおよび組織学的画像の2次元アレイに対するマッピングを有する出力画像を生成することであって、前記出力画像が、複数の組織クラスのうちの1つを各ピクセルに割り当てることによって生成され、前記複数の組織クラスが、非腫瘍組織を表す少なくとも1つのクラスおよび腫瘍組織を表す少なくとも1つのクラスを含む、ことと、
    前記出力画像からセグメンテーションマスクを生成することであって、前記セグメンテーションマスクが、個々の腫瘍によって占有された領域を特定する、ことと、
    各腫瘍の要約統計量を計算することと、
    前記セグメンテーションマスクをフィルタリングして、各腫瘍の前記要約統計量に従って個々の腫瘍を選択および選択解除することと、
    編集されたセグメンテーションマスクに従って前記組織学的画像のビジュアライゼーションを作成することと、
    前記フィルタリングされたセグメンテーションマスクで個々の腫瘍の選択を受信することと、
    現在選択されている腫瘍に対して追加の計算診断プロセスを実行することと、
    を含む、方法。
  66. 前記追加の計算診断プロセスが、処理タスクを前記追加の計算診断プロセスに送信することと、前記処理タスクの結果を前記追加の計算診断から受信することと、を含む、
    請求項65記載の方法。
  67. 前記追加の計算診断プロセスが、データ通信ネットワークを介して通信可能に結合されたプロセッサ対応デバイス上で実行される、
    請求項66記載の方法。
  68. 前記処理タスクが、前記組織学的画像を含む、
    請求項66記載の方法。
  69. 前記処理タスクが、前記組織学的画像に関連するメタデータを含む、
    請求項68記載の方法。
  70. 前記組織学的画像に関連する前記メタデータが、前記セグメンテーションマスクを含む、
    請求項68記載の方法。
  71. 前記組織学的画像に関連する前記メタデータが、患者情報、染色情報、染色プロトコル情報、および走査プロトコル情報のうちの1つ以上を含む、
    請求項68記載の方法。
JP2021538332A 2019-05-29 2020-05-29 組織学的画像および術後腫瘍辺縁評価における腫瘍のコンピュータ支援レビュー Pending JP2022534157A (ja)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US201962854014P 2019-05-29 2019-05-29
US201962854055P 2019-05-29 2019-05-29
US62/854,014 2019-05-29
US62/854,055 2019-05-29
PCT/US2020/035286 WO2020243545A1 (en) 2019-05-29 2020-05-29 Computer supported review of tumors in histology images and post operative tumor margin assessment

Publications (1)

Publication Number Publication Date
JP2022534157A true JP2022534157A (ja) 2022-07-28

Family

ID=71787097

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021538332A Pending JP2022534157A (ja) 2019-05-29 2020-05-29 組織学的画像および術後腫瘍辺縁評価における腫瘍のコンピュータ支援レビュー

Country Status (6)

Country Link
US (2) US11893732B2 (ja)
EP (1) EP3884465A1 (ja)
JP (1) JP2022534157A (ja)
KR (1) KR20220015368A (ja)
CN (1) CN113366530A (ja)
WO (1) WO2020243545A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7312510B1 (ja) 2022-08-16 2023-07-21 之江実験室 腫瘍微小環境を考慮した全スライド病理画像分類システム及び構築方法

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11893732B2 (en) * 2019-05-29 2024-02-06 Leica Biosystems Imaging, Inc. Computer supported review of tumors in histology images and post operative tumor margin assessment
AU2020344517A1 (en) * 2019-09-09 2022-03-24 PAIGE.AI, Inc. Systems and methods for processing images of slides for digital pathology
CN111325220B (zh) * 2020-02-17 2023-04-07 腾讯科技(深圳)有限公司 图像生成方法、装置、设备及存储介质
KR102406080B1 (ko) * 2020-09-01 2022-06-07 성균관대학교산학협력단 경계특성 기반의 병변 분석 장치 및 방법
US20220164576A1 (en) * 2020-11-25 2022-05-26 Metal Industries Research & Development Centre Surgical instrument inventory system and surgical instrument inventory method
CN112967267B (zh) * 2021-03-23 2024-01-23 湖南珞佳智能科技有限公司 一种全卷积神经网络的激光定向能量沉积溅射计数方法
FR3126528B1 (fr) * 2021-08-30 2023-12-22 Totalenergies Se Procédé et dispositif électronique de surveillance de la fabrication d’un liquide avec un additif anti-mousse, système de surveillance et programme d’ordinateur associés
WO2023073405A1 (en) * 2021-10-28 2023-05-04 Cloud Pathology Group Srl Digital microscopy tissue image analysis method and system for digital pathology
CN115240000B (zh) * 2022-07-22 2023-05-02 司法鉴定科学研究院 一种用于法医鉴定的糖尿病检测装置及方法
CN115063425B (zh) * 2022-08-18 2022-11-11 合肥综合性国家科学中心人工智能研究院(安徽省人工智能实验室) 基于读片知识图谱的结构化检查所见生成方法及系统
CN115294126B (zh) * 2022-10-08 2022-12-16 南京诺源医疗器械有限公司 一种病理图像的癌细胞智能识别方法
CN116046647B (zh) * 2023-01-28 2023-06-09 深圳安侣医学科技有限公司 血液成像分析系统和方法
CN116230214B (zh) * 2023-05-08 2023-08-01 浙江大学滨江研究院 Hcc及vetc辅助诊断装置及设备
CN116664722B (zh) * 2023-07-27 2023-10-13 山东新卓诚智慧科技有限公司 一种智慧环保的大数据可视化处理方法及系统

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102043133B1 (ko) * 2012-11-16 2019-11-12 삼성전자주식회사 컴퓨터 보조 진단 지원 장치 및 방법
US9430829B2 (en) 2014-01-30 2016-08-30 Case Western Reserve University Automatic detection of mitosis using handcrafted and convolutional neural network features
US10169685B2 (en) * 2014-07-07 2019-01-01 The Regents Of The University Of California Automatic segmentation and quantitative parameterization of brain tumors in MRI
EP3185000B1 (en) * 2014-08-22 2021-12-01 National University Corporation Nagoya University Three-dimensional image construction method, image processing device, and electron microscope
WO2017070585A1 (en) * 2015-10-23 2017-04-27 Novartis Ag Method of scoring a sample comprising tumor tissue
EP3721373A1 (en) * 2017-12-07 2020-10-14 Ventana Medical Systems, Inc. Deep-learning systems and methods for joint cell and region classification in biological images
US11348239B2 (en) * 2018-05-14 2022-05-31 Tempus Labs, Inc. Predicting total nucleic acid yield and dissection boundaries for histology slides
US11348661B2 (en) * 2018-05-14 2022-05-31 Tempus Labs, Inc. Predicting total nucleic acid yield and dissection boundaries for histology slides
AU2020274091A1 (en) * 2019-05-14 2021-12-09 Tempus Ai, Inc. Systems and methods for multi-label cancer classification
US20210142904A1 (en) * 2019-05-14 2021-05-13 Tempus Labs, Inc. Systems and methods for multi-label cancer classification
US11893732B2 (en) * 2019-05-29 2024-02-06 Leica Biosystems Imaging, Inc. Computer supported review of tumors in histology images and post operative tumor margin assessment

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7312510B1 (ja) 2022-08-16 2023-07-21 之江実験室 腫瘍微小環境を考慮した全スライド病理画像分類システム及び構築方法

Also Published As

Publication number Publication date
KR20220015368A (ko) 2022-02-08
WO2020243545A1 (en) 2020-12-03
US11893732B2 (en) 2024-02-06
EP3884465A1 (en) 2021-09-29
US20220076410A1 (en) 2022-03-10
CN113366530A (zh) 2021-09-07
US20240119595A1 (en) 2024-04-11

Similar Documents

Publication Publication Date Title
US20220076411A1 (en) Neural netork based identification of areas of interest in digital pathology images
US11893732B2 (en) Computer supported review of tumors in histology images and post operative tumor margin assessment
US11164316B2 (en) Image processing systems and methods for displaying multiple images of a biological specimen
US20220084660A1 (en) Artificial intelligence processing system and automated pre-diagnostic workflow for digital pathology
US10650221B2 (en) Systems and methods for comprehensive multi-assay tissue analysis
AU2015265811B2 (en) An image processing method and system for analyzing a multi-channel image obtained from a biological tissue sample being stained by multiple stains
WO2017093524A1 (en) Devices, systems and methods for automated quantitative scoring of digitized tissue images
WO2022126010A1 (en) Systems and methods for generating histology image training datasets for machine learning models

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20211130

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230512

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20240313

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20240423