JP2022133479A - 予後推定装置及び予後推定方法 - Google Patents

Abstract

【課題】肺がんの予後を精度よく推定する。【解決手段】予後推定装置１は、分析対象の肺がん組織を撮像した対象画像について、所定の閾値を用いて閾値処理を行った画像を準備する閾値処理部２１と、閾値処理後の画像について、所定の個別関心領域毎に組織の穴及び連結部分の数をそれぞれカウントすることにより、ベッチマップを作成するベッチマップ作成部２２と、ベッチマップから、予後に関連する特徴量を算出する、特徴量算出部２３と、を有する。【選択図】図１

Description

本発明は、予後推定装置及び予後推定方法に関する。

がん患者の患部のＣＴ画像を撮像し、画像処理技術を用いてがんに係る判定を行う方法が検討されている。例えば、特許文献１では、組織の画像を解析することで、画像にがん組織の像が含まれているかを判定する方法が開示されている。また、非特許文献１では、テクスチャ特徴およびウェーブレット分解等を使用してがん組織の画像に係る特徴量を求め、肺がんの予後を推定することが検討されている。

特許第５５２２４８１号公報

Arimura H, Soufi M, Kamezawa H,Ninomiya K and Yamada M "Radiomics with artificial intelligence forprecision medicine in radiation therapy" Journal of Radiation Research, September2018, 60 (1), Pages 150-157

しかしながら、肺がんの予後の推定については精度に関して改善の余地があった。

本発明は上記を鑑みてなされたものであり、肺がんの予後を精度よく推定することが可能な予後推定装置及び予後推定方法を提供することを目的とする。

上記目的を達成するため、本発明の一形態に係る予後推定装置は、分析対象の肺がん組織を撮像した対象画像について、所定の閾値を用いて閾値処理を行った画像を準備する閾値処理部と、閾値処理後の画像について、所定の個別関心領域毎に組織の穴及び連結部分の数をそれぞれカウントすることにより、ベッチマップを作成するベッチマップ作成部と、前記ベッチマップから、予後に関連する特徴量を算出する、特徴量算出部と、を有する。

また、本発明の一形態に係る予後推定方法は、分析対象の肺がん組織を撮像した対象画像について、所定の閾値を用いて閾値処理を行った画像を準備する画像準備ステップと、閾値処理後の画像について、所定の個別関心領域毎に組織の穴及び連結部分の数をそれぞれカウントすることにより、ベッチマップを作成するベッチマップ作成ステップと、前記ベッチマップから、予後に関連する特徴量を算出する、特徴量算出ステップと、を含む。

上記の予後推定装置及び予後推定方法によれば、分析対象の肺がん組織を撮像した対象画像について、閾値に基づく閾値処理を行った画像を準備し、当該画像について、所定の関心領域毎に組織の穴及び連結部分の数をそれぞれカウントすることにより、ベッチマップを作成し、このベッチマップに基づいて予後に関連する特徴量が算出される。このようにがん組織の穴及び連結部分に係るベッチマップを作成した上で予後に関連する特徴量を算出する構成とすることで、肺がんの予後を精度よく推定することが可能となる。

前記閾値処理部において、１の前記対象画像から、閾値を変更した複数種類の閾値処理を行った画像を準備し、前記ベッチマップ作成部において、前記複数種類の閾値処理を行った画像から、複数の前記ベッチマップを作成し、前記特徴量算出部において、前記複数のベッチマップに基づいて特徴量を算出する態様とすることができる。

上記のように、１の対象画像から閾値を変更した複数種類の閾値処理を行った画像を準備し、複数のベッチマップを作成した上で特徴量を算出する構成とすることで、分析対象の肺がん組織に含まれる種々の穴等に係る情報を反映したベッチマップを作成することができると共に、これらの情報を考慮して予後に関連する特徴量を算出することができる。したがって、肺がんの予後をより精度よく推定することが可能となる。

前記特徴量算出部において、予後の予測に係る数理モデルに対して所定のアルゴリズムを複数回適用することで、予後に関連する特徴量を選択する態様とすることができる。

上記のように数理モデルに対して所定のアルゴリズムを適用して特徴量を選択する際に、これを複数回適用して特徴量を選択することにより、特徴量の選択をより精度良く行うことができるため、肺がんの予後をより精度よく推定することが可能となる。

本発明によれば、肺がんの予後を精度よく推定することが可能な予後推定装置及び予後推定方法が提供される。

図１は、一実施形態に係る予後推定装置の構成を説明するブロック図である。 図２は、予後推定装置において行われる予後の推定に係る処理を説明するフロー図である。 図３（ａ）、図３（ｂ）、図３（ｃ）は、予後推定装置が取得する画像の例を説明する図である。 図４は、閾値処理を行った画像を説明する図である。 図５は、ベッチマップ画像の作成手順について説明するフロー図である。 図６は、ベッチマップ画像の作成手順について説明する図である。 図７は、閾値毎に作成されたベッチマップ画像を説明する図である。 図８（ａ）、図８（ｂ）、図８（ｃ）、図８（ｄ）は、異なる手順で算出した特徴量に基づいてがん患者を２群に分けた際の、２群の分離具合について説明する図である。

以下、添付図面を参照して、本発明を実施するための形態を詳細に説明する。なお、図面の説明においては同一要素には同一符号を付し、重複する説明を省略する。

図１は、本発明の一形態に係る予後推定装置の構成を説明するブロック図である。本実施形態に係る予後推定装置１は、肺がん患者の患部の画像から肺がんの予後を推定する装置である。対象となる肺がんには、非小細胞肺がん（腺がん、扁平上皮がん、大細胞がん、腺扁平上皮がん等）、及び、小細胞肺がんの両方が含まれる。以下の実施形態では、これらを総称して肺がんという。予後推定装置１では、肺がん患者の患部組織（がん組織）を撮像した画像に基づいて肺がんの予後に関連性の高い特徴量（radiomic features）を算出して出力する機能を有する。予後推定装置１では、肺がんの予後と関連性が高い特徴量を算出することができるため、画像に基づいて特徴量の算出結果が、肺がんの予後の推定結果に対応するといえる。

予後推定装置１における肺がんの予後の推定に利用する肺がん患者の患部組織を撮像した画像（対象画像）としては、例えば、ＣＴ（Computed Tomography：コンピュータ断層撮影）により得られるＣＴ画像が挙げられる。予後推定装置１では、肺がんの患部組織及びその周囲を撮像したＣＴ画像に対して所定の処理を行うことで、予後に関係する特徴量を算出することができる。

図１に示すように、予後推定装置１は、画像取得部１１、画像解析部１２、記憶部１３、及び、出力部１４を有する。

予後推定装置１は、例えば、ＣＰＵ（Central Processing Unit）、主記憶装置であるＲＡＭ（Random Access Memory）及びＲＯＭ（Read Only Memory）、他の機器との間の通信を行う通信モジュール、並びにハードディスク等の補助記憶装置等のハードウェアを備えるコンピュータとして構成される。そして、これらの構成要素が動作することにより、予後推定装置１としての機能が発揮され得る。なお、予後推定装置１としての機能は、複数のコンピュータの組み合わせによって実現されていてもよい。

予後推定装置１の画像取得部１１は、外部装置等から患者の患部組織を撮像したＣＴ画像である対象画像を取得する機能を有する。また、画像取得部１１は、ＣＴ画像に含まれる患部組織（肺がん組織）の領域を特定する情報も併せて取得してもよい。一般的に、患者のＣＴ画像として、肺がん組織及びその周辺の組織を撮像した画像が取得される。したがって、画像取得部１１では、ＣＴ画像と共に当該ＣＴ画像のうち肺がん組織を撮像した領域を特定する情報を取得する。これにより、画像取得部１１は、肺がん組織を撮像した画像（肺がん組織として特定された像が含まれる画像）を取得することができる。なお、ＣＴ画像自体から肺がん組織の領域を特定することができる場合には、肺がん組織の領域を特定する情報を別途取得しなくてもよい。

画像解析部１２は、肺がん組織を撮像したＣＴ画像についての解析を行い、肺がんの予後に係る特徴量を算出する機能を有する。画像解析部１２で行われる解析の詳細については後述するが、ＣＴ画像を所定の閾値に基づいた二値化処理を行った後、肺がん組織に含まれる孔の数に係るベッチ数（Betti numbers）の分布に基づくベッチマップを作成する。そして、このベッチ数の分布に基づいてがんの予後に関連する特徴量を算出する。そのため、画像解析部１２は、閾値処理部２１、ベッチマップ作成部２２、及び、特徴量算出部２３を有する。

記憶部１３は、画像取得部１１において取得された画像及びがん組織の領域を特定する情報、画像解析部１２による解析結果等を記憶する機能を有する。

出力部１４は、画像解析部１２による解析結果を出力する機能を有する。出力先としては、予後推定装置１に設けられたモニタ、外部装置等が挙げられる。また、出力内容等は特に限定されず、例えば、画像解析部１２により算出された特徴量をそのまま出力する態様としてもよいし、特徴量に基づいた予後の判定を行いその結果を出力する態様としてもよい。予後の判定を行う場合には、判定を行うためのロジック等を記憶部１３において保持し、特徴量を算出した後に当該ロジックを用いて判定を行うこととしてもよい。

次に、図２～図７を参照しながら、予後推定装置１による画像に基づく予後推定の手順について説明する。図２は、予後推定装置１において行われる予後の推定に係る処理を説明するフロー図である。この図２を参照しながら、予後の推定に係る処理について説明する。

まず、予後推定装置１では、画像取得部１１により、肺がん組織を含むＣＴ画像と、当該ＣＴ画像における肺がん組織（腫瘍組織）を撮像した領域を特定する情報を取得する（ステップＳ０１）。画像取得部１１が取得する情報の例を図３に示す。図３（ａ）は、肺がん組織を含む領域を撮像したＣＴ画像の例である。図３（ｂ）は、ＣＴ画像のうち、肺がん組織の像が含まれる領域を特定した情報の例であり、白色部分が肺がん組織であることを示している。図３（ｃ）は、図３（ａ）で示すＣＴ画像に対して図３（ｂ）で示す肺がん組織の場所を特定する情報を重ね合わせたものである。画像取得部１１において図３（ｃ）で示す情報を取得することで、予後推定装置１では、ＣＴ画像のうち肺がん組織を撮像した領域を特定することができ、当該領域に含まれる情報を予後の推定に利用することができる。

次に、予後推定装置１の画像解析部１２の閾値処理部２１では、前処理としてＣＴ画像を８ビットのグレースケール画像に変換する（ステップＳ０２）。この処理を行うことにより、２５６段階の光の強さで表示されたグレースケールのＣＴ画像を得ることができる。その後、複数の閾値を用いて二値化処理を行い、閾値処理画像を作成する（ステップＳ０３：画像準備ステップ）。上記の通り、ＣＴ画像は２５６段階（０～２５５）の光の強度の組み合わせで構成されている。したがって、そのうちいずれかの段階を閾値として設定し、光の強度が閾値またはそれ以下の領域を黒とし、光の強度が閾値よりも大きい領域を白とすることで、二値化された画像を得ることができる。０～２５５のそれぞれを閾値として二値化処理を行うことで、２５６枚の二値化画像を得ることができる。図４は上記の処理に得られた二値化画像を閾値（threshold value）に並べたものを模式的に示している。図４に示すように、閾値が０の場合には、画像全体が白くなる。一方、閾値を大きくすると徐々に画像に含まれる黒色の領域が大きくなる。また、図示していないが閾値が２５５の場合には、画像全体が黒くなる。このように１枚のＣＴ画像から２５６枚の閾値処理を行った画像（閾値処理画像）を生成する。

次に、予後推定装置１のベッチマップ作成部２２では、閾値処理画像毎にベッチマップを作成する（ステップＳ０４：ベッチマップ作成ステップ）。ベッチマップとは、閾値処理画像に含まれるがん組織の穴の数に係るベッチ数及び連結部分の数に係るベッチ数に基づいて作成された画像である。ベッチマップの作成については、図５及び図６を参照しながら説明する。

まず、上述の通り、閾値処理画像を準備する（ステップＳ１１）。次に、準備した閾値処理画像について、ベッチ数を算出するための領域毎に画像を分割する（ステップＳ１２）。図６では、閾値を７０として二値化した画像Ｄ１を示している。この画像Ｄ１を例えばカーネルサイズを縦横共に７ピクセルとして分割する。なお、以降の処理を行う対象となる領域は、肺がん組織を撮像した領域である。この領域を関心領域（ＲＯＩ：regions of interests）という場合があり、分割された領域の一つずつを個別関心領域（local ROI）という場合がある。

次に、ベッチマップ作成部２２では、分割された領域毎にベッチ数（Bettinumbers）ｂ０及びｂ１を算出した画像（ｂ０画像及びｂ１画像）を生成する（ステップＳ１３、ステップＳ１４）。ベッチ数とは、位相空間に対する不変量であり、自然数で記述される。本実施形態では、ベッチ数ｂ０（０次のベッチ数）は分割された領域における連結成分の数を示す。ベッチ数ｂ１（１次のベッチ数）は分割された領域における穴の数を示す。すなわち、ベッチ数ｂ１，ｂ０を算出することは、所定の関心領域毎に組織の穴及び連結部分の数をそれぞれカウントすることに対応するといえる。

図６では、一部の関心領域ＲＯＩに対応する画像Ｄ２及び画像Ｄ３を反転した画像を示している。画像Ｄ２のうち黒色となっている領域は、がん組織が存在しないため暗くなっている、すなわち、穴となっている領域を示している。画像Ｄ２では、白色の領域である連結部分（白色の画素同士が連結されている部分）が１つとなっている。すなわち、画像Ｄ２では１つの連結部分のみで構成されている。したがって、この画像Ｄ２に基づくベッチ数ｂ０は１となる。

また、画像Ｄ３では、穴となっている領域が白色で示されているが、この領域が３つ独立した状態（穴となっている領域に対して隣接する４辺の画素が全て黒色である状態）で存在している。したがって、この画像Ｄ３（画像Ｄ２と同じ）に基づくベッチ数ｂ１は３となる。このように、分割した領域毎にベッチ数ｂ０，ｂ１をそれぞれ算出する。分割した領域毎にベッチ数は算出されるので、各領域でのベッチ数ｂ０の算出結果をまとめたｂ０画像（図６の画像Ｄ４）及び各領域でのベッチ数ｂ１の算出結果をまとめたｂ１画像（図６の画像Ｄ５）を得ることができる。ｂ０画像の作成（Ｓ１３）及びｂ１画像の作成（Ｓ１４）は、独立した工程であるので、同時に行ってもよいし、いずれかを先に行ってもよい。

次に、ベッチマップ作成部２２では、ｂ０画像及びｂ１画像から、ｂ１／ｂ０画像を作成する（ステップＳ１５）。上述のようにベッチ数ｂ０及びｂ１は、いずれも自然数であるので、ｂ１／ｂ０を算出することができる。上記のように分割された領域毎にベッチ数がｂ０，ｂ１が算出されるので、これらの数値に基づいて領域毎にｂ１／ｂ０を算出する。算出されたｂ１／ｂ０の分布を図示したのがｂ１／ｂ０画像である。図６では、ｂ０画像Ｄ４及びｂ１画像Ｄ５に基づいて作成されたｂ１／ｂ０画像Ｄ６を示している。このｂ１／ｂ０画像Ｄ６が、ベッチマップ画像に対応する。

ベッチマップ画像は、１つの二値化画像から１つ生成することができる。すなわち、本実施形態で説明したように２５６個の二値化画像（閾値処理画像）を作成した場合、２５６個のベッチマップを作成することができる。

図２に戻り、予後推定装置１の特徴量算出部２３では、上記の手順で得られた複数のベッチマップ画像から、予後推定用の特徴量（radiomic features）を算出する（ステップＳ０５：特徴量算出ステップ）。予後推定用の特徴量とは、ベッチマップ画像から算出される特徴量であり肺がんの予後と関連性が高いと考えられる特徴量である。予後推定用の特徴量を算出するための手法としては、従来、Soufi M, Arimura H, Nakamoto T, Hirose T aki, Ohga S, Umezu Y, HondaH and Sasaki T "Exploration of temporal stability and prognostic power ofradiomic features based on electronic portal imaging device images" PhysicaMedica, February 2018, 46, Pages 32-44等で開示されている特徴量の算出手法を適用することができる。従来開示されている特徴量の算出手法とは、８ビットグレースケールに変換した１枚のＣＴ画像からウェーブレット関数に基づく４種類の画像を準備し、この４種類の画像に基づいて５４個の特徴量を算出するものである。特徴量を算出する手法自体は上記の文献に記載されている。また、本実施形態でも上記と同様の方法により１枚のｂ１／ｂ０画像から５４種類の特徴量を得ることができる。本実施形態の場合、例えば、２５６個のベッチマップ画像に対して上記の特徴量の算出を行った結果、合計１３８２４の特徴量（homological radiomic features）を得ることができる。

なお、本実施形態では、予後の予測に係る数理モデルを用いて、上記の１３８２４の特徴量から予後の推定に適した特徴量として７種類の特徴量を選択している。具体的には、ＣＰＨＭ（コックス比例ハザードモデル：Cox Proportional-Hazards Model）を用いている。このモデルはある患者が他の患者が長生きしているかを評価する関数であり、患者の生存時間を予測するものである。このモデルに対して、ｅｌａｓｔｉｃ－ｎｅｔというスパースモデリングを適用したＣｏｘ－ｎｅｔアルゴリズムを適用することで、特徴量を選択する処理を行った。Ｃｏｘ－ｎｅｔアルゴリズム自体は予後の推定に用いられている公知のものであり、例えば、Simon N, Friedman J, Hastie T and Tibshirani R "RegularizationPaths for Cox’s Proportional Hazards Model via Coordinate Descent" Journalof Statistical Software, Mar 2011, 39(5), Pages.1-13に具体的に記載されている。このＣｏｘ－ｎｅｔアルゴリズムを用いることで、肺がんの予後の推定に適した特徴量を５４種類の特徴量から選択することができる。

また、本実施形態では、１３８２４の特徴量に対してＣｏｘ－ｎｅｔアルゴリズムを１００回適用することで、患者の生存時間の予測（すなわち、予後の推定）との相関がより強いと考えられる７個の特徴量を選択している。Ｃｏｘ－ｎｅｔアルゴリズムを１００回適用することで、より多様な条件において予後との相関が高い特徴量を特定することができる。

なお、上記の手順で算出された特徴量が予後の予測に有用かどうかを評価した結果を図８に示す。具体的には、上記実施形態で説明した特徴量に基づいて肺がん患者を２群に分けた場合、その予後に差があるかを評価した。また、従来から行われている特徴量の算出方法等、上記実施形態で説明した方法とは異なる予後の推定方法に基づいて肺がん患者を２群に分けた場合と比較して、群間の差が大きくなるかを評価した。

図８（ａ）では、本実施形態で説明したベッチマップを利用して算出した特徴量を用いて２群に分けた場合の結果を示している。また、図８（ｂ）では、従来から行われているウェーブレット関数に基づく４種類の画像から得られる特徴量を用いて２群に分けた場合の結果を示している。また、図８（ｃ）では、本実施形態で説明したベッチマップを利用して算出した特徴量と従来からのウェーブレット関数を用いて得られる特徴量との組み合わせた場合の結果を示している。また、図８（ｄ）では、深層学習（Deep Learning）を用いて患者を２群に分けた場合の結果を示している。図８（ａ）～図８（ｄ）では、同一の肺がん患者の画像を利用してそれぞれ２群に分けている。

肺がん患者を特徴量に基づいて２群に分ける方法について説明する。本実施形態に係る手法に基づいてがん患者を２群に分ける場合、まず、ＣＰＨＭを作成すると、７個の特徴量それぞれに対する重み付けに係る重み係数が得られる。患者毎にＣＴ画像から得られる７つの特徴量それぞれに関する重み係数と特徴量との積の総和を算出する。このように算出されるradiomic scores（rad-scores）の中央値を用いて、肺がん患者を２群に分けた。なお、図８（ｂ）、（ｃ）に示す結果に対応する手法についても同様の方法を用いて肺がん患者を２群に分けた。また、図８（ｄ）に示す結果に対応する手法については、深層学習を行う際に利用した特徴量を利用して患者を２群に分けた。

図８（ａ）～図８（ｄ）では、２群それぞれの患者について、経過年に対する生存率をプロットしている。各図では、２群のうち、生存年が長くなった患者群をＬｏｗとし、生存年が短くなった患者群をＨｉｇｈとしている。また、各図では、ｐ値を示している。ｐ値は、２群の分離度合いを示す値であって、値が小さいほど２群が精度良く分離されていることを示している。具体的には、図８（ａ）に示す結果に係るｐ値は１．５×１０^－３であり、図８（ｂ）に示す結果に係るｐ値は０．３１であり、図８（ｃ）に示す結果に係るｐ値は２．９×１０^－３であり、図８（ｄ）に示す結果に係るｐ値は０．０７９であった。

図８（ａ）～図８（ｄ）に示すように、本実施形態で説明したベッチマップを利用して算出した特徴量による２つのがん患者群は、他の手法に基づいて分けられた２つのがん患者群よりも経過年に対する生存率という点で精度よく分離できていることが示されている。すなわち、本実施形態で説明した手法で算出される特徴量及びこの特徴量に基づくradiomic scores（rad-scores）は、予後の推定に対して有用であるといえる。

画像解析部１２の特徴量算出部２３では、ベッチマップ画像を用いて特徴量を算出するだけでなく、当該画像の元となるＣＴ画像を撮像したがん患者がＨｉｇｈ群またはＬｏｗ群のいずれに属するかを判定するような処理を行ってもよい。

なお、上記の７つの特徴量に対して重み係数を反映してradiomic scores（rad-scores）を算出する方法は、上記の特徴量（radiomic features）が予後の推定に関して有用であるかを評価する際に用いた手法の一例である。すなわち、実際に予後が不明ながん患者について予後を推定する場合には、上記とは異なる手法を用いてもよい。すなわち、ベッチマップ画像に基づいた特徴量を算出した後、上記とは異なる方法（例えば、特徴量に関する何らかの計算を行う等）で特徴量に基づいてがん患者の予後を推定してもよい。

図２に戻り、予後推定装置１の出力部１４では、推定結果を出力する（ステップＳ０６）。出力する結果とは、予後の推定時に利用した特徴量自体であってもよいし、上記の図８等で示したように、がん患者を２群に分けた場合にどちらに属するかを特定する情報であってもよい。

上記のように、本実施形態に係る予後推定装置及び予後推定方法によれば、分析対象の肺がん組織を撮像した対象画像について、閾値に基づく閾値処理を行った画像を準備し、当該画像について、所定の関心領域毎に組織の穴及び連結部分の数をそれぞれカウントすることにより、ベッチマップを作成し、このベッチマップに基づいて予後に関連する特徴量が算出される。このようにがん組織の穴及び連結部分に係るベッチマップを作成した上で予後に関連する特徴量を算出する構成とすることで、肺がんの予後を精度よく推定することが可能となる。

従来から、肺がん組織を撮像した画像（特に、ＣＴ画像）から予後を推定する方法は検討されている。また、肺がんの場合、従来から肺がん組織に形成された穴と患者の予後に相関がある可能性が指摘されていた。したがって、予後の推定に適した特徴量を算出できれば、予後の推定を精度良く行うことができると考えられていた。しかしながら、従来から行われた方法では、撮像した画像から予後の推定に対応すると考えられる特徴量を算出することができるものの、推定精度の向上が望まれていた。

これに対して、本実施形態に係る予後推定装置及び予後推定方法では、組織の穴及び連結部分の数に着目したベッチマップを作成している。このベッチマップに用いられるベッチ数とは、位相幾何学における位相空間に対する不変量であり、これを利用することで、肺がん組織の画像に含まれる組織の穴及び連結部分に係る不変量を利用した特徴量の算出が可能となる。したがって、従来の手法（例えば、ベッチマップを利用しない方法）と比較して、予後の推定に係る精度の高い特徴量を算出することができる。したがって、予後の推定を精度よく行うことが可能となる。

また、上記の予後推定装置１では、１の対象画像から閾値を変更した複数種類の閾値処理を行った画像を準備し、複数のベッチマップを作成した上で特徴量を算出する構成とされている。このような構成とすることで、分析対象の肺がん組織に含まれる種々の穴等に係る情報を反映したベッチマップを作成することができると共に、これらの情報を考慮して予後に関連する特徴量を算出することができる。したがって、肺がんの予後をより精度よく推定することが可能となる。

なお、上記実施形態では、８ビットのグレースケール画像における２５６段階全てを利用して、２５６種類の閾値処理を行った画像を作成しているが、この数を変更したとしても予後の推定精度を高くすることができると考えられる。例えば、光の強度が０～２５５である２５６段階全てではなく、１，３，５，７・・・というように一部の光の強度のみを閾値として用いて複数種類の閾値処理を行った画像を作成してもよい。この場合、複数種類の画像の閾値は、特定の光の強度範囲に集中して選択してもよいし、ある程度分散させて選択してもよい。

また、ベッチマップの作成に使用する閾値処理を行った画像は１種類であってもよい。すなわち、１つの二値化画像のみを用いてベッチマップを作成して特徴量を算出してもよい。

また、上記の予後推定装置１では、特徴量算出部２３において、予後の予測に係る数理モデル（ＣＰＨＭ）に対して所定のアルゴリズム（Ｃｏｘ－ｎｅｔアルゴリズム）を複数回適用することで、予後に関連する特徴量を選択している。数理モデルに対して所定のアルゴリズムを適用して特徴量を選択する際に、これを複数回適用して特徴量を選択することにより、特徴量の選択をより精度良く行うことができるため、肺がんの予後をより精度よく推定することが可能となる。なお、上記の実施形態では、ＣＰＨＭに対してＣｏｘ－ｎｅｔアルゴリズムを複数回適用する場合について説明したが、数理モデル及びアルゴリズムは、上記で説明したものに限定されない。すなわち、予後の予測に係る数理モデルであれば、ＣＰＨＭ以外の数理モデルを利用してもよい。また、数理モデルの選択に対応させて生存時間の分析に係るアルゴリズムから適切なアルゴズムを選択して組み合わせる構成としてもよい。

本発明は、上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合せにより種々の発明を形成できる。

例えば、上記実施形態では、患者の肺がん組織を撮像した画像（対象画像）がＣＴ画像である場合について説明したが、画像の種類は限定されない。

また、ベッチマップを作成する際の分割する領域の大きさ（画素数等）は、画像の大きさ及び細かさ等を考慮して適宜変更してもよい。また、ベッチマップの作成方法は上記に限定されない。

また、１の患者について１の画像を用いて特徴量を算出する構成について説明したが、例えば、複数の画像を用いて特徴量を算出する構成としてもよい。

また、予後推定装置１は、少なくとも予後の推定に関連する特徴量を算出する機能を有していればよい。すなわち、予後推定装置１とは異なる装置で、特徴量から具体的な予後を算出する処理等を行ってもよい。

１…予後推定装置、１１…画像取得部、１２…画像解析部、１３…記憶部、１４…出力部、２１…閾値処理部、２２…ベッチマップ作成部、２３…特徴量算出部。

Claims (4)

1. 分析対象の肺がん組織を撮像した対象画像について、所定の閾値を用いて閾値処理を行った画像を準備する閾値処理部と、
   閾値処理後の画像について、所定の個別関心領域毎に組織の穴及び連結部分の数をそれぞれカウントすることにより、ベッチマップを作成するベッチマップ作成部と、
   前記ベッチマップから、予後に関連する特徴量を算出する、特徴量算出部と、
   を有する、予後推定装置。
2. 前記閾値処理部において、１の前記対象画像から、閾値を変更した複数種類の閾値処理を行った画像を準備し、
   前記ベッチマップ作成部において、前記複数種類の閾値処理を行った画像から、複数の前記ベッチマップを作成し、
   前記特徴量算出部において、前記複数のベッチマップに基づいて特徴量を算出する、請求項１に記載の予後推定装置。
3. 前記特徴量算出部において、予後の予測に係る数理モデルに対して所定のアルゴリズムを複数回適用することで、予後に関連する特徴量を選択する、請求項２に記載の予後推定装置。
4. 分析対象の肺がん組織を撮像した対象画像について、所定の閾値を用いて閾値処理を行った画像を準備する画像準備ステップと、
   閾値処理後の画像について、所定の個別関心領域毎に組織の穴及び連結部分の数をそれぞれカウントすることにより、ベッチマップを作成するベッチマップ作成ステップと、
   前記ベッチマップから、予後に関連する特徴量を算出する、特徴量算出ステップと、
   を含む、予後推定方法。

Images