(57)【要約】  組織に生じた変化の程度を精度良く評価する。組織を撮像した１つの撮像画像に対して、二値化の基準値が異なる複数の二値化画像を生成し、複数の二値化画像のそれぞれについて、二値化後の一方の画素値の画素が連結して成る連結領域の数を示す第１特徴数を算出するベッチ数算出部（４２）と、二値化の基準値と第１特徴数との関係性を示すパターンを所定の基準パターンと比較した結果に基づいて、組織に生じた変化に関する情報を決定する予測スコア決定部（４４）とを備える。
2019102829.tif 000002
【特許請求の範囲】
【請求項１】
組織を撮像した撮像画像を解析する方法であって、
前記撮像画像から、二値化の基準値が異なる複数の二値化画像を生成する二値化ステップと、
複数の前記二値化画像のそれぞれについて、二値化後の一方の画素値の画素が連結して成る連結領域の数を示す第１特徴数を算出する算出ステップと、
前記二値化の基準値と前記第１特徴数との関係性を示す第１パターンを、所定の基準パターンと比較した結果に基づいて、前記組織に生じた変化に関する情報を決定する決定ステップと、を含むことを特徴とする画像解析方法。
【請求項２】
前記二値化ステップにて生成した複数の二値化画像のそれぞれについて、前記一方の画素値の画素に囲まれた、他方の画素値の画素からなる穴形状の領域の数を示す第２特徴数を算出し、
前記決定ステップにおいて、前記二値化の基準値の変化に応じて前記第２特徴数が変化する第２パターンを、所定の基準パターンと比較した結果に基づいて、前記組織に生じた変化に関する情報を決定することを特徴とする請求項１に記載の画像解析方法。
【請求項３】
前記組織に生じた変化に関する確定情報が予め対応付けられた前記撮像画像を取得して、該撮像画像の前記第１パターンと、前記確定情報との組み合わせを含む学習データを用いた機械学習によって、前記基準パターンを生成する学習ステップをさらに含むことを特徴とする請求項１または２に記載の画像解析方法。
【請求項４】
前記確定情報は、前記組織に生じた変化の程度に応じた複数の段階に関する情報を含み、
前記基準パターンは、前記段階ごとに生成されており、
前記決定ステップにおいて、前記組織に生じた変化が前記段階のいずれに該当するかを決定することを特徴とする請求項３に記載の画像解析方法。
【請求項５】
前記撮像画像から解析対象となる組織に対応する領域を抽出する抽出ステップをさらに含むことを特徴とする請求項１から４のいずれか１項に記載の画像解析方法。
【請求項６】
組織を撮像した撮像画像を解析する画像解析装置であって、
前記撮像画像から、二値化の基準値が異なる複数の二値化画像を生成する二値化部と、
複数の二値化画像のそれぞれについて、二値化後の一方の画素値の画素が連結して成る連結領域の数を示す第１特徴数を算出する算出部と、
前記二値化の基準値と前記第１特徴数との関係性を示す第１パターンを、所定の基準パターンと比較する比較部と、
前記比較の結果に基づいて、前記組織に生じた変化に関する情報を決定する決定部と、を備えることを特徴とする画像解析装置。
【請求項７】
前記算出部は、前記二値化部によって生成された複数の二値化画像のそれぞれについて、前記一方の画素値の画素に囲まれた、他方の画素値の画素からなる穴形状の領域の数を示す第２特徴数を算出し、
前記比較部は、前記二値化の基準値の変化に応じて前記第２特徴数が変化する第２パターンを、所定の基準パターンと比較し、
前記決定部は、比較結果に基づいて、前記組織に生じた変化に関する情報を決定することを特徴とする請求項６に記載の画像解析装置。
【請求項８】
前記組織に生じた変化に関する確定情報が予め対応付けられた前記撮像画像を取得して、該撮像画像の前記第１パターンと、前記確定情報との組み合わせを含む学習データを用いた機械学習によって、前記比較部に前記基準パターンを生成させる学習制御部をさらに備えることを特徴とする請求項６または７に記載の画像解析装置。
【請求項９】
前記確定情報は、前記組織に生じた変化の程度を複数の段階に区分する情報を含み、
前記基準パターンは、前記段階ごとに生成されており、
前記決定部は、前記組織に生じた変化が前記段階のいずれに該当するかを決定することを特徴とする請求項８に記載の画像解析装置。
【請求項１０】
前記撮像画像から解析対象となる組織に対応する領域を抽出する領域抽出部をさらに備えることを特徴とする請求項６から９のいずれか１項に記載の画像解析装置。
【請求項１１】
請求項６から１０のいずれか１項に記載の画像解析装置と、
前記撮像画像の画像データを前記画像解析装置へ送信する外部機器と、
前記画像解析装置において決定された情報を取得して該情報を提示する提示装置と、を含むことを特徴とする画像解析システム。
【請求項１２】
請求項６に記載の画像解析装置としてコンピュータを機能させるための画像解析プログラムであって、前記二値化部、前記算出部、前記比較部、および前記決定部としてコンピュータを機能させるための画像解析プログラム。
【請求項１３】
請求項１２に記載の画像解析プログラムを記録したコンピュータ読み取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は組織を撮像した画像を解析して、組織に生じた変化の程度を評価するための情報を出力する画像解析方法および画像解析装置等に関する。
【背景技術】
【０００２】
ＣＴ（computed tomography）は、医療画像診断において必要不可欠な技術である。ＣＴ画像を用いた画像診断が主に用いられている臓器の一つとして肺が挙げられる。肺に生じた疾患を診断する場合、被検体の胸部を撮像した胸部ＣＴ画像が用いられる。【０００３】
従来、医師などが視覚的に診断を行っていたが、近年では、胸部ＣＴ画像を解析した結果に基づいた診断が行われている。例えば、肺のＣＴ画像のスクリーニングは、慢性閉塞性肺疾患（Chronic obstructive pulmonary disease (COPD)、以後、「ＣＯＰＤ」と略記する）などを早期に発見・診断し、治療を行うことを可能にしており、日本のみならず世界全体において健康寿命を維持する上で重要な技術である。【０００４】
例えば、胸部ＣＴ画像の画像解析によって肺に生じたＣＯＰＤ、肺気腫、および肺癌などの病変について診断する場合に用いられている定量指標の１つにＬＡＡ％が挙げられる。ＬＡＡ％を用いた画像診断では、胸部ＣＴ画像に含まれる肺の上部、中部、および下部の３領域において、Ｘ線の吸収率が低い（ＣＴ値が−９５０ＨＵ以下）領域である低吸収肺領域（low-attenuation lung area）が占める割合が指標として用いられる。【０００５】
肺気腫は、肺胞が拡大して、Ｘ線を吸収する肺胞壁が萎縮消失し、肺が持続的に拡張する状態である。そのため、肺のＣＴ画像において、肺気腫が生じた領域は正常な肺の画像に比べてより黒い領域となる。ＬＡＡ％を適用する場合、例えば、胸部ＣＴ画像に含まれる肺の領域が抽出され、ＣＴ値が−９５０ＨＵなどを閾値として二値化される。そして、肺の領域の総画素数に占める低吸収肺領域（すなわち、二値化の結果、画素値が０の領域）のピクセル数の割合（％）を算出する。この割合は、肺気腫の重症度に相関性が認められる値として知られている。【０００６】
肺のＣＴ画像を診断する放射線医および呼吸器内科医などは、肺気腫が肺癌のリスク要因であることを指摘している。非特許文献１および２に記載されているように、肺気腫検出のための定量指標としてＬＡＡ％を適用することは有効である。しかし、ＬＡＡ％は肺気腫の重症度を正確に評価できていない場合も多い。それゆえ、ＬＡＡ％を用いた肺気腫の重症度判定に基づいて肺癌に進行するリスクを評価することは難しい、との指摘がなされている。
【先行技術文献】
【非特許文献】
【０００７】
【非特許文献１】Smith BM et al., 「Emphysema detected on computed tomographyand risk of lung cancer: a systematic review and meta-analysis」 Lung Cancer, Vol.77(1), p58-63, 2012
【非特許文献２】Wille MMW et al., 「Visual assessment of early emphysema andinterstitial abnormalities on CT is useful in lung cancer risk analysis」 Eur Radiol. Vol.26(2), p487-494, 2016
【発明の概要】
【発明が解決しようとする課題】
【０００８】
医師が肺のＣＴ画像を視覚的に診断して肺気腫から肺癌へ進行するリスクを評価できるのと同様に、肺気腫から肺癌に進行するリスクを評価することが可能な画像解析技術が求められている。【０００９】
一般に、画像診断では、多くの場合、癌化のリスクに関する評価が重要である。そのため、画像解析技術としては、組織に生じた病変について構造上の特性を捉えて、その病変がさらに重症度の疾患へと進行する可能性を判定するための情報を提供することが重要である。肺気腫から肺癌に進行するリスクを評価する場合に限らず、さまざまな臓器および組織においても同様である。【００１０】
本発明の一態様は、上記の問題点を解決するためになされたものであり、その目的は、組織の画像を解析することにより、組織に生じた変化の程度を医師のように精度良く評価するための情報を出力する画像解析方法および画像解析装置等を実現することを目的とする。
【課題を解決するための手段】
【００１１】
上記の課題を解決するために、本発明の一態様に係る画像解析方法は、組織を撮像した撮像画像を解析する方法であって、前記撮像画像から、二値化の基準値が異なる複数の二値化画像を生成する二値化ステップと、複数の前記二値化画像のそれぞれについて、二値化後の一方の画素値の画素が連結して成る連結領域の数を示す第１特徴数を算出する算出ステップと、前記二値化の基準値と前記第１特徴数との関係性を示す第１パターンを、所定の基準パターンと比較した結果に基づいて、前記組織に生じた変化に関する情報を決定する決定ステップと、を含む。【００１２】
また、本発明の一態様に係る画像解析装置は、組織を撮像した撮像画像を解析する画像解析装置であって、前記撮像画像から、二値化の基準値が異なる複数の二値化画像を生成する二値化部と、複数の二値化画像のそれぞれについて、二値化後の一方の画素値の画素が連結して成る連結領域の数を示す第１特徴数を算出する算出部と、前記二値化の基準値と前記第１特徴数との関係性を示す第１パターンを、所定の基準パターンと比較する比較部と、前記比較の結果に基づいて、前記組織に生じた変化に関する情報を決定する決定部と、を備える。【００１３】
また、本発明の一態様に係る画像解析システムは、前記の画像解析装置と、前記撮像画像の画像データを前記画像解析装置へ送信する外部機器と、前記画像解析装置において決定された情報を取得して該情報を提示する提示装置と、を含む。
【発明の効果】
【００１４】
本発明の一態様によれば、組織の画像を解析することにより、組織に生じた変化の程度を精度良く評価するための情報を出力することができる。
【図面の簡単な説明】
【００１５】
【図１】本発明の実施形態１に係る画像解析装置の構成例を示す機能ブロック図である。
【図２】本発明の一態様に係る画像解析方法の処理の流れの一例を示すフローチャートである。
【図３】図２のステップＳ１における処理の一例を示すフローチャートである。
【図４】学習に用いる胸部ＣＴ画像の例を示す図である。
【図５】胸部ＣＴ画像から肺のＣＴ画像を切り出して、異なる二値化の基準値を用いて生成された二値化画像の例を示す図である。
【図６】二値化の基準値に対して、０次元ベッチ数および１次元ベッチ数が変化する様子をプロットしたグラフである。
【図７】解析対象の胸部ＣＴ画像から予測スコアを決定する処理の流れの一例を示すフローチャートである。
【図８】本発明の実施形態２に係る画像解析装置の構成例を示す機能ブロック図である。
【図９】本発明に係る画像解析装置を含む画像解析システムの概略構成の例を示す図である。
【図１０】本発明の一態様に係る画像解析方法を胸部ＣＴ画像に適用して、肺癌に罹患している患者を評価した場合の精度を、他の指標を適用した場合の精度と比較した表である。
【図１１】年齢・性別・喫煙履歴（ブリンクマン指数）という基本の指標を用いた場合、基本の指標に加えてＬＡＡ％を指標として用いた場合、および基本の指標に加えてベッチ数を指標として用いた場合における、肺癌の予測精度を比較した結果を示す表である。
【発明を実施するための形態】
【００１６】
〔実施形態１〕
以下、本発明の一実施形態について、詳細に説明する。
【００１７】
（本発明の技術思想）
まず、本発明の技術思想について以下に説明する。本発明の目的は、組織を撮像した撮像画像を解析することにより、当該撮像画像に写っている組織に生じた病変の有無、および癌化の可能性の程度を精度良く判別するために有用な情報を出力することである。
【００１８】
本発明の発明者は、被検体の胸部を撮像したＣＴ（computed tomography）画像を調べ、肺気腫の重症度と癌化の可能性との関係を詳細に比較・検討した。具体的には、胸部ＣＴ画像のうち、肺が写っている領域に対して、二値化の基準値が異なる複数の二値化画像を生成した。そして、二値化の基準値に応じて、二値化後の各画像に含まれる肺のＣＴ画像について算出された０次元ベッチ数（ｂ０、第１特徴数）、および１次元ベッチ数（ｂ１、第２特徴数）がどのように変化するかを調べた。その結果、肺のＣＴ画像において、二値化の基準値を変化させたことに伴う、０次元ベッチ数および１次元ベッチ数の変化はいずれも、肺気腫の重症度を判別したり、癌化の可能性を評価したりするための指標として利用することが可能であることを見出した。
【００１９】
ところで、０次元ベッチ数と１次元ベッチ数のいずれか一方を利用する構成であってもよいが、以下では、０次元ベッチ数と１次元ベッチ数の双方を利用する構成を採用した場合について説明する。０次元ベッチ数と１次元ベッチ数の双方を利用することにより、より精度の高い情報を出力することができる。
【００２０】
なお、本明細書では、二値化の基準値は各画素のＣＴ値に基づいて設定されるものとする。すなわち、二値化の基準値は、肺を撮像対象とする胸部ＣＴ画像において一般的に設定されるＣＴ値−１０００〜−７００ＨＵの範囲において設定される。例えば、二値化の基準値を−９００ＨＵに設定した場合、ＣＴ値が−９００ＨＵ以下の画素の画素値は０となり、ＣＴ値が−９００より高い画素の画素値を１となる。
【００２１】
以下では、Ｘ線ＣＴ画像を解析して、肺に生じた肺気腫の重症度を判別するために有用な情報を出力する画像解析装置１を例に挙げて説明するが、本発明の解析対象は、肺に限定されない。また、Ｘ線ＣＴ画像に限定されるものでもない。本発明の一態様に係る画像解析装置１は、任意の医用画像に適用可能であり、例えば、（１）レントゲン（マンモグラフィーを含む）検査にて撮像される肝臓および乳房などの画像、（２）ＭＲＩ（磁気共鳴像、magnetic resonance imaging）検査にて撮像される脳および膵臓などの画像、（３）ＰＥＴ（陽電子放射断層撮影、positron emission tomography）、ＲＩ検査（アイソトープ検査とも呼称される）にて撮像される画像の解析にも好適に適用され得る。
【００２２】
＜肺のＣＴ画像＞
例えば、肺のＣＴ画像において正常な肺は淡い乳白色の靄がかかった黒っぽい領域として表示されるのに対し、肺気腫を患っている部位では肺胞壁が萎縮消失しているため靄のうすい黒い領域として表示される。
【００２３】
一般に肺のＣＴ画像において、肺胞の破壊および融合は、その初期において、正常な肺に囲まれた低吸収肺領域として検出される。正常な肺の領域と低吸収肺領域との境界には明瞭な壁が無いことが多いが、肺気腫の初期段階では、直径が約１ｃｍ程度の低吸収領域が散在するようになり、重症化するにつれて、低吸収領域が融合して増加する一方、正常な肺は減少する。
【００２４】
胸部ＣＴ画像の画像解析によって肺に生じた肺気腫などの病変について診断する場合、ＬＡＡ％等の定量指標が適用され得る。ＬＡＡ％は、肺のＣＴ画像において、低吸収肺領域（low-attenuation lung area）が全肺領域に占める割合（％）を算出し、この割合に基づいて、肺気腫の重症度を判定する。ここで、重症度とは肺に生じた変化の程度を意図し、健常な肺には無い病変の進行具合を指す。
【００２５】
医師が、胸部ＣＴ画像に基づいて肺気腫の有無および重症度を判断する場合、撮像された肺に生じた肺気腫の重症度に応じて複数の段階に区分された視覚的スコア（visual score）が確定情報として胸部ＣＴ画像に対応付けられる。
【００２６】
以下では、肺気腫の所見が無い０、および肺気腫の重症度が１（軽度）、２、３（中度）、４、５（重度）の計６段階に肺気腫の重症度を区分する視覚的スコアを適用した場合を例に挙げて説明する。
【００２７】
なお、視覚的スコアには幾つかの種類があり、これに限定されるものではない。例えば、肺のＣＴ画像において、６つの部位における視覚的評価（０〜４点の５段階）を合計して肺気腫の重症度を段階別に区分するGoddard分類の肺気腫スコア（満点は２４点）を用いてもよい。
【００２８】
＜肺のＣＴ画像を解析するための数学的表現＞
組織に生じた変化を定量化するために、本発明の発明者らは、ホモロジーの概念、特にパーシステントホモロジーの適用を試みた。ホモロジーとは、図形の形態上の性質を代数的に置き換えて、図形の結合などの解析を容易にする数学の一分野である。本発明の発明者らは、ホモロジーの概念の中でも、二次元画像における０次元ベッチ数および１次元ベッチ数の利用に着目した。
【００２９】
ホモロジーの概念は、構成要素の連結および接触を表す数学的な概念である。病理画像において、適切な二値化の基準値（二値化パラメータともいう）を設定して二値化する。そして、二値化された画像から単位面積当たりの０次元ベッチ数および１次元ベッチ数を算出する。算出された０次元ベッチ数および１次元ベッチ数を用いて、その組織の構成要素同士の連結の程度、および構成要素同士の接触の程度の評価が可能である。
【００３０】
ベッチ数とは、図形（例えば、組織の構成要素に該当）の形状には無関係であり、図形同士の接触と分離とにのみ関係するトポロジカルな示唆数である、ｑ次元特異ホモロジー群が有限生成のとき、このｑ次特異ホモロジー群は、自由アーベル群と有限アーベル群との直和に分けられる。この自由アーベル群の階級をベッチ数という。２次元画像の場合、０次元ベッチ数は連結成分の数であり、１次元ベッチ数は当該連結成分が外縁となって囲んだ空間（以下、「穴形状の領域」と記す場合がある）の数である。穴形状の領域の数とは、すなわち連結成分中に存在する「穴」の数である。
【００３１】
＜０次元ベッチ数＞
０次元ベッチ数は、数学的には以下のように定義される。一般に有限個の線分を繋ぎ合わせて成る図形（１次元複体とも呼称される）Ｋの連結成分の個数を０次元ベッチ数という。「有限個の点を有限個の線分で結んだ図形が連結である」とは、この図形の任意の頂点から他の任意の頂点に、この図形の辺を辿って到達し得ることを意図する。
【００３２】
肺のＣＴ画像に対して、異なる二値化の基準値を用いて生成した複数の二値化画像のそれぞれにおいて、二値化後の一方の画素値（例えば、二値化の結果、画素値が０）の画素が連結して成る連結領域の数が、０次元ベッチ数である。それゆえ、０次元ベッチ数は、肺のＣＴ画像における低吸収肺領域の数に相当するといえる。
【００３３】
＜１次元ベッチ数＞
１次元ベッチ数は、数学的には以下のように定義される。以下の（１）および（２）の条件が満たされる場合、図形Kの１次元ベッチ数はｒである。（１）有限個の線分を繋ぎ合わせて成る図形（連結な１次元複体）Ｋに対して、適当なｒ個の、開いた（両端を含まない）１次元単体（例えば、線分）を図形Ｋから取り去っても図形Ｋの連結成分の個数は増加しない。（２）任意の（ｒ＋１）個の、開いた１次元単体をＫから取り去った場合にはＫは連結でなくなる（すなわち、Ｋの連結成分の個数が１つ増加する）。
【００３４】
肺のＣＴ画像に対して、異なる二値化の基準値を用いて生成した複数の二値化画像のそれぞれにおいて、二値化後の一方の画素値（例えば、二値化の結果、画素値が０）の画素に囲まれた穴形状の領域の数が、１次元ベッチ数である。それゆえ、１次元ベッチ数は、肺のＣＴ画像における、低吸収肺領域に囲まれた正常な肺の領域の数に相当するといえる。
【００３５】
（画像解析方法の概要）
まず、本発明の一態様に係る画像解析方法は、被検体を撮像した撮像画像を解析し、撮像されている組織に生じた変化に関する情報を決定・出力する方法である。本発明の一態様に係る画像解析方法の概要について図２を用いて説明する。図２は、本発明の一態様に係る画像解析方法における処理の流れの一例を示すフローチャートである。なお、以下では、被検体の胸部ＣＴ画像を解析して、肺気腫などの重症度に関する情報を出力する画像解析装置１が行う処理を例に挙げて説明する。
【００３６】
なお、以下では、０次元ベッチ数および１次元ベッチ数を算出する構成について説明するが、いずれか一方のみを算出して利用する構成であってもよい。ただし、組織に生じた変化の特性は様々であり、撮像画像に表れる特徴も様々であるため、０次元ベッチ数および１次元ベッチ数を算出する構成が望ましい。なぜなら、組織に生じた変化を、撮像画像中の連結領域の数の変化として検出することに加え、穴形状の数の変化としても検出する構成である方が、精度の良い予測スコアを決定することができるからである。
【００３７】
まず、医師などによって過去に診察された被検体を撮像した胸部ＣＴ画像であって、肺気腫の重症度および肺癌への進行に関するリスクが既に医師などによって評価されている胸部ＣＴ画像を用いた機械学習が行われる（ステップＳ１）。ステップＳ１では、過去に診察された被検体を撮像した胸部ＣＴ画像が、ステップＳ３において解析対象のＣＴ画像に対して行われる解析と同様に解析される。そして、肺のＣＴ画像を解析することで得られる特徴量と、既に決定されている視覚的スコアとの相関関係が学習される。ここで、特徴量とは、ＣＴ画像についての二値化の基準値に応じた０次元ベッチ数および１次元ベッチ数の変化を示す数値、該変化をプロットしたグラフにおいて表される増減パターン、および該グラフを近似的に表す関数の係数値などであってもよい。
【００３８】
なお、本明細書では、肺気腫の視覚的スコアを予測スコアとして決定する場合を例に挙げて説明するが、本発明の適用範囲は視覚的スコアの決定に限定されない。例えば、ＣＴ、ＭＲＩ、ＰＥＴ、ＲＩなどを用いて撮像された断層撮像を解析することで得られる特徴量と、既に決定されている該病変に関する臨床的意義のある任意の指標との相関関係が学習され得る。臨床的意義のある任意の指標としては、例えば、癌化のリスク、癌の重症度（例えば、ＴＮＭステージングなどの病期分類）、および癌の予後予測などが挙げられる。このような学習を行うことによって、画像解析装置１は、被検体の画像を解析して、臨床的に意義のある任意の指標に関する情報を出力することができる。
【００３９】
ここで、「予測スコア」とは、医師によって決定される視覚的スコアに倣った指標（スコア）であり、肺のＣＴ画像を解析した結果として得られた、肺気腫の重症度および肺癌への進行に関するリスク評価に利用可能なスコアである。ステップＳ１の詳細な工程については後に具体例を挙げて説明する（図３参照）。
【００４０】
次に、被検体の胸部ＣＴ画像が入力される（ステップＳ２）。
【００４１】
続いて、入力された胸部ＣＴ画像を解析して、肺気腫の重症度および肺癌への進行に関するリスク評価に利用可能な予測スコアが決定される（ステップＳ３）。ステップＳ３の詳細な工程については後に具体例を挙げて説明する（図７参照）。
【００４２】
最後に、ステップＳ３にて決定された予測スコアが出力される（ステップＳ４）。予測スコアを出力する態様は、予測スコアの利用を希望するユーザに提示されるものであればよく、特に限定されない。例えば、表示装置５に表示されてもよいし、プリンター（図示せず）、スピーカー（図示せず）から出力されてもよい。
【００４３】
（画像解析装置１の構成）
次に、図２に示す画像解析方法を実現する画像解析装置１の構成について、図１を用いて説明する。図１は、本発明の実施例１に係る画像解析装置１の構成の一例を示すブロック図である。なお、図１において、図２のステップＳ１に関連するデータの流れは破線で示し、図２のステップＳ２〜４に関連するデータの流れは実線で示している。
【００４４】
図１に示すように、画像解析装置１は、画像取得部２（受信部）、記憶部３、制御部４を備えている。記憶部３には少なくとも、ＣＴ画像３１および基準パターン３２が記憶されている。
【００４５】
なお、画像解析装置１が別体として設けられた表示装置５（提示部）と接続されている例を図示しているが、これに限定されない。例えば、画像解析装置１が表示装置５（提示部）を内蔵する構成であってもよいし、複数の表示装置５と無線／有線で接続された構成であってもよい。
【００４６】
画像取得部２は、外部機器８から組織を撮像した撮像画像（以下、ＣＴ画像と称する）を取得し、取得したＣＴ画像を記憶部３に格納する。あるいは、画像取得部２は、取得したＣＴ画像を二値化部４１に出力してもよい。なお、外部機器８は、例えば、ＣＴ装置に接続されたパーソナルコンピュータ、画像データを保存可能な記憶装置などである。上記ＣＴ画像は、解析対象の組織が肺である場合、被検体の胸部を適当なコントラストで撮像して得られたＣＴ画像である。なお、ＣＴ画像のコントラストなどは、解析の対象となる組織のＸ線吸収率に応じて適宜設定されればよい。
【００４７】
画像解析装置１が肺のＣＴ画像に基づいてＣＯＰＤに関する診断を行う場合、画像取得部２は、肺の上部、中部、および下部の断面をそれぞれ含む少なくとも３枚の胸部ＣＴ画像を一組として取得することが望ましい。
【００４８】
また、画像取得部２は、医師によって決定された、組織に関する視覚的スコアが確定情報として予め対応付けられたＣＴ画像とその視覚的スコアとの組み合わせを取得し、記憶部３に格納する。視覚的スコア（確定情報）が予め対応付けられたＣＴ画像は、予測スコア決定部４４の機械学習において入力される学習データとして用いられる。
【００４９】
記憶部３は、画像取得部２が取得したＣＴ画像３１、図２のステップＳ１で生成した基準パターン３２に加え、制御部４が実行する（１）各部の制御プログラム、（２）ＯＳプログラム、（３）アプリケーションプログラム、および、（４）これらプログラムを実行するときに読み出す各種データを記録するものである。記憶部３は、ハードディスク、フラッシュメモリなどの不揮発性の記憶装置によって構成される。なお、記憶部３の他に、上述の各種プログラムを実行する過程でデータを一時的に保持するための作業領域として使用される記憶装置であり、ＲＡＭ（Random Access Memory）などの揮発性の記憶装置を備えていてもよい。
【００５０】
表示装置５は、制御部４から出力される出力情報等を表示する表示装置であり、例えば液晶ディスプレイである。なお、画像解析装置１が専用の表示装置５を備える構成であってもよい。なお、表示装置５の表示画面にタッチセンサが重畳されており、ユーザによる表示面へのタッチ操作を表示装置５が検出できる構成であってもよい。
【００５１】
（制御部４の構成）
制御部４は、画像取得部２が取得したＣＴ画像を解析することにより、当該ＣＴ画像に含まれる肺気腫などＣＯＰＤの重症度を示す情報を出力する。この制御部４は、さらに二値化部４１（二値化部、領域抽出部）、ベッチ数算出部４２（算出部）、パターン特定部４３、予測スコア決定部４４（比較部、決定部）、表示制御部４５、および学習制御部４６を備えている。
【００５２】
二値化部４１は、画像取得部２が取得した胸部ＣＴ画像から、被検体の肺に対応する領域を肺のＣＴ画像として切り出し（抽出し）、切り出された肺のＣＴ画像に対して二値化処理を行う。すなわち、二値化部４１が領域抽出部（図示せず）としての機能も備えている。被検体の肺に対応する領域は、肺のＸ線吸収率が低いことを利用した公知の画像処理の手法を用いて胸部ＣＴ画像から切り出すことが可能である。なお、被検体の肺に対応する領域を肺のＣＴ画像として切り出して、切り出された肺のＣＴ画像を二値化部４１に出力する領域抽出部を二値化部４１とは別に備える構成としてもよい。
【００５３】
図１は、二値化部４１が胸部ＣＴ画像を記憶部３から取り出す構成について図示しているが、これに限定されない。例えば、画像取得部２が取得した解析対象となる胸部ＣＴ画像を記憶部３に記憶させることなく、二値化部４１に出力する構成であってもよい。
【００５４】
二値化部４１は、二値化処理において、二値化の基準値よりも大きいＣＴ値を有する画素を白色の画素に変換し、二値化の基準値以下のＣＴ値を有する画素を黒色の画素に変換する。このとき二値化部４１は、肺を撮像した１つのＣＴ画像に対して、二値化の基準値を変更する度に二値化処理を行い、複数の二値化画像を生成する。すなわち、二値化部４１は、肺のＣＴ画像について、二値化の基準値が異なる複数の二値化画像を生成する。
【００５５】
ベッチ数算出部４２は、二値化部４１が生成した複数の二値化画像について、当該二値化画像に含まれる、（１）二値化後の画素値が０である画素が連結して成る連結領域の数（２）二値化後の画素値が０である画素に囲まれた穴形状の領域の数を算出する。すなわち、ベッチ数算出部４２は、複数の二値化画像にそれぞれ含まれる組織の像について、０次元ベッチ数および１次元ベッチ数を算出する。
【００５６】
上記連結領域は、画素値が０である画素が隣接しつつ集まった領域である。それぞれの連結領域は、画素値が１である画素によって囲まれており、互いに独立している。
【００５７】
上記穴は、１つ以上の構成成分の外縁の少なくとも一部（構成成分が１つの場合には、その構成成分の外縁の全部）をその外縁として有している開口部である。換言すれば、構成成分の外縁によって囲まれた穴には、１つの構成成分がその内部に有する穴と、互いに連結した複数の構成成分がそれぞれ有する外縁によって囲まれた穴とが含まれる。
【００５８】
ベッチ数算出部４２として、既存のプログラムを用いてもよい。このようなプログラムの一例としてＣＨｏｍＰを挙げることができる。ＣＨｏｍＰは、ＧＮＵ（General PublicLicense）に準拠したフリーウェアであり、ホームページ（http://chomp.rutgers.edu/）からダウンロードできる。なお、これに限定されるものではなく、画像に関する０次元ベッチ数および１次元ベッチ数を算出できるプログラムであれば、ＣＨｏｍＰ以外のプログラムを用いてもよい。
【００５９】
パターン特定部４３は、二値化の基準値に応じた０次元ベッチ数の変化のパターン（第１パターン）、および二値化の基準値に応じた１次元ベッチ数の変化のパターン（第２パターン）のうちの少なくとも一方を特定する。例えば、パターン特定部４３は、二値化の基準値に対して、該二値化の基準値を設定して生成した二値化画像において算出された０次元ベッチ数および１次元ベッチ数の少なくとも何れかをプロットしたグラフを生成する。この場合、０次元ベッチ数の変化について生成されたグラフは第１パターンであり、１次元ベッチ数の変化について生成されたグラフは第２パターンである。
【００６０】
なお、パターン特定部４３は、二値化の基準値の変化に伴う０次元ベッチ数および１次元ベッチ数の変化を近似的に示す関数を特定してもよい。この場合、パターン特定部４３は、特定した関数の係数などを、０次元ベッチ数および１次元ベッチ数の変化のパターンとして出力してもよい。ここで、０次元ベッチ数および１次元ベッチ数の変化を近似的に示す関数（近似関数）としては、任意の関係式で表記可能な関数であってもよく、例えば、ｎ次関数（ｎ≧１）（線形近似および対数近似を含む）、指数関数（累乗近似）、および三角関数などであってもよい。
【００６１】
予測スコア決定部４４は、パターン特定部４３によって特定された第１パターンを、記憶部３に記憶されている基準パターン３２と比較し、比較結果に基づいて、肺に生じた変化に関する情報を決定する。
【００６２】
ここで、基準パターン３２とは、過去に医師などによって診察され、既に視覚的スコアが決定されている胸部ＣＴ画像を、その視覚的スコア毎に解析することでパターン生成部によって生成される特徴量である。基準パターン３２は、図１のステップＳ１の学習によって生成される。具体的には、基準パターン３２は、パターン特定部４３によって生成された、二値化の基準値に応じた０次元ベッチ数および１次元ベッチ数の変化をプロットしたグラフである。基準パターン３２は、視覚的スコアが複数の段階に区分されている場合、その段階ごとに生成され得る。
【００６３】
すなわち、予測スコア決定部４４は、パターン特定部４３によって特定された第１パターンを、視覚的スコアの段階ごとに生成された基準パターン３２と比較し、最も類似した基準パターン３２と同じスコアを予測スコアとして決定する。
【００６４】
（学習ステップＳ１の処理の流れ）
次に、図２のステップＳ１における処理の流れの一例について、図４〜６を参照しながら、図３を用いて説明する。図３は、学習ステップＳ１における処理の流れの一例を示すフローチャートである。
【００６５】
画像取得部２は、外部機器８から視覚的スコアなどの診断結果が確定されている胸部ＣＴ画像を取得する（ステップＳ１１）。画像取得部２は、取得した胸部ＣＴ画像を記憶部３に格納する。
【００６６】
図４に、学習に用いる胸部ＣＴ画像の例を示す。図４の（Ａ）は肺気腫が疑われる所見が無く、正常の肺（視覚的スコア＝０）と診断された画像の例であり、（Ｂ）は、中度の肺気腫（視覚的スコア＝３）が認められる画像の例であり、（Ｃ）は、重度の肺気腫（視覚的スコア＝５）が認められる画像の例である。
【００６７】
次に、学習制御部４６は、ＣＴ画像３１から胸部ＣＴ画像を１つ選択し、二値化部４１に送る。この胸部ＣＴ画像には、予め決定された視覚的スコアが対応付けられている。二値化部４１は、肺に対応する肺のＣＴ画像を切り出す（抽出する）（ステップＳ１２：抽出ステップ）。さらに、二値化部４１は、肺のＣＴ画像に対して、二値化の基準値を変化させながら、肺のＣＴ画像から複数の二値化画像を生成する（ステップＳ１３）。
【００６８】
図５は、胸部ＣＴ画像から肺のＣＴ画像を切り出して、異なる二値化の基準値を用いて生成された二値化画像の例を示す図である。なお、図５では、肺の領域を分かりやすく示すために、白と黒を反転させた画像を示している。図５の（Ａ）は、図４の（Ｃ）に示した胸部ＣＴ画像と同じものであり、（Ｂ）〜（Ｅ）は、（Ａ）の画像から、二値化の基準値を−９７５ＨＵ〜−９００ＨＵに設定した場合に生成される二値化画像の例を示している。
【００６９】
次に、ベッチ数算出部４２は、各二値化画像について０次元ベッチ数（ｂ０）および１次元ベッチ数（ｂ１）を算出する（ステップＳ１４）。パターン特定部４３は、二値化の基準値と０次元ベッチ数との関係性を示すグラフ、および二値化の基準値と１次元ベッチ数との関係性を示すグラフを生成する（ステップＳ１５）。すなわち、二値化の基準値の変化に対する０次元ベッチ数の値の変化を示すグラフと、二値化の基準値の変化に対する１次元ベッチ数の値の変化を示すグラフとを生成する。
【００７０】
図６は、二値化の基準値に対して、０次元ベッチ数および１次元ベッチ数が変化する様子をプロットしたグラフである。図６の左側に示すグラフはいずれも、横軸に二値化の基準値（すなわち、二値化の閾値として設定した値）としたＣＴ値（ＨＵ）、縦軸に０次元ベッチ数をプロットしたグラフである。また、図６の右側に示すグラフはいずれも、横軸に二値化の基準値としたＣＴ値（ＨＵ）、縦軸に１次元ベッチ数をプロットしたグラフである。図６の（Ａ）〜（Ｃ）はそれぞれ、図４の（Ａ）〜（Ｃ）に示した胸部ＣＴ画像から生成された二値化画像について生成されたグラフである。なお、図６では、０次元ベッチ数または１次元ベッチ数の値を肺領域の総画素数で正規化した値を縦軸にプロットしているが、これに限定されない。例えば、肺領域の総画素数で正規化を施さなくても、同様の学習効果が得られる場合、正規化は必須の処理ではない。なお、学習において、敢えて正規化を施さないようにすることで予測精度を向上させることが可能であることを、本発明の発明者は確認している。
【００７１】
学習制御部４６は、パターン特定部４３から生成されたグラフを取得し、胸部ＣＴ画像に対応付けられていた診断結果（例えば、視覚的スコア）とともに予測スコア決定部４４に送る。学習制御部４６は、予測スコア決定部４４に、グラフを入力変数、視覚的スコアを目的変数とする学習を行わせる（ステップＳ１６）。なお、この学習のアルゴリズムとしては任意の機械学習用のアルゴリズムを適用することができる。このようなアルゴリズムの一例として、「ランダムフォレスト」が挙げられる。
【００７２】
図６に示す６つのグラフに示すように、二値化の基準値に応じた０次元ベッチ数または１次元ベッチ数の値の変化のパターン（すなわち、グラフの形態上の特徴）は、肺気腫の重症度によって異なっている。例えば、０次元ベッチ数のグラフの極大値を与える二値化の基準値を見れば、肺気腫が重症化するにつれて低下する傾向があることが分かる。一方、グラフのピーク幅を見れば、肺気腫が重症化するにつれて広くなる傾向が分かる。また、肺気腫が重症化すると、グラフに変曲点が出現する傾向があることも分かる。学習制御部４６は、予測スコア決定部４４に、このようなグラフの特徴を入力変数とし、視覚的スコアを目的変数とする学習を行わせる。
【００７３】
学習制御部４６は、学習用に格納されたＣＴ画像３１のすべてが選択されるまでステップＳ１２〜ステップＳ１６を繰り返す（ステップＳ１７）。予測スコア決定部４４にすべてのＣＴ画像３１についての学習が終わった場合、学習制御部４６は、パターン特定部４３によって生成されたグラフを基準パターンのグラフとして、視覚的スコア毎に基準パターン３２に格納する。
【００７４】
（予測スコアを決定する処理の流れ）
続いて、図２のステップＳ３における処理の流れの一例について、図７を用いて説明する。図７は、解析対象の胸部ＣＴ画像から予測スコアを決定する処理の流れの一例を示すフローチャートである。
【００７５】
ステップＳ３１において、二値化部４１は、記憶部３から胸部ＣＴ画像を１つ取得し、被検体の肺に対応する肺のＣＴ画像を切り出す（ステップＳ３１）。さらに、二値化部４１は、肺のＣＴ画像に対して、二値化の基準値を変化させながら、肺のＣＴ画像から複数の二値化画像を生成する（ステップＳ３２：二値化ステップ）。
【００７６】
次に、ベッチ数算出部４２は、各二値化画像について０次元ベッチ数（ｂ０）および１次元ベッチ数（ｂ１）を算出する（ステップＳ３３：算出ステップ）。パターン特定部４３は、二値化の基準値に応じて０次元ベッチ数および１次元ベッチ数が変化する様子をプロットしたグラフを生成する（ステップＳ３４）。
【００７７】
予測スコア決定部４４は、パターン特定部４３によって生成されたグラフと、予め学習で得た、０次元ベッチ数および１次元ベッチ数が変化する様子をプロットした基準パターン３２のグラフとを比較し、比較結果に基づいて予測スコアを決定する。予測スコア決定部４４は、例えば、グラフの極大値を与える二値化の基準値、ピーク幅、変曲点の有無および位置などを基準パターン３２のグラフと比較することにより、グラフの特徴が最も類似している基準パターン３２のグラフに対応付けられている視覚的スコアを、予測スコアとして決定する（ステップＳ３５：決定ステップ）。
【００７８】
医師などが、肺を撮像した撮像画像（例えば、ＣＴ画像）に基づいて肺気腫などの疾患の有無および重症度を決定する場合、例えば、以下の（１）および（２）のような特徴が手掛かりとして利用される。（１）軽度の肺気腫では、直径約１ｃｍ程度の低吸収領域の散在が認められる。（２）重度の肺気腫では、低吸収領域は融合して増加する一方、正常な肺の領域は減少する。
【００７９】
図２および７に示す方法により、組織に生じた変化の有無および程度に関する情報を精度良く決定することができる。例えば、決定した情報を医師に提示することによって、組織に生じた変化の有無および程度の判断を適切にサポートすることができる。
【００８０】
〔実施形態２〕
本発明の他の実施形態について、以下に説明する。なお、説明の便宜上、上記実施形態にて説明した部材と同じ機能を有する部材については、同じ符号を付記し、その説明を繰り返さない。
【００８１】
（画像解析装置１ａの構成）
次に、図８を用いて、画像解析装置１ａの構成について説明する。図８は、本発明の実施例２に係る画像解析装置１ａの構成の一例を示すブロック図である。画像解析装置１ａは通信部６を備えている点で、画像解析装置１と異なっている。通信部６は、外部機器８から受信した画像データが示す撮像画像を解析した結果に基づいて決定した予測スコアを予測スコア決定部４４から取得し、提示装置７に送信する。
【００８２】
なお、外部機器８および提示装置７の数は１つに限定されるものではなく、複数であってもよい。
【００８３】
（画像解析システム）
ここでは、画像解析装置１ａを含む画像解析システム１００および１００ａの構成例について、図９を用いて説明する。図９は、本発明に係る画像解析装置１ａを含む画像解析システム１００および１００ａの構成例を示す概略図である。図９の（ａ）は、外部機器８と提示装置７とが離れた場所に設置されている例を示しており、（ｂ）は、提示装置７が外部機器８ａに接続されている例を示している。
【００８４】
画像解析システム１００は、外部機器８、画像解析装置１ａ、および提示装置７を有している。外部機器８、画像解析装置１ａ、および提示装置７は、インターネットなどの情報通信ネットワーク５０に接続されており、相互にデータの送受信が可能である。
【００８５】
外部機器８は、例えば、ＣＴ装置に接続されたパーソナルコンピュータであってもよいし、組織を撮像した画像を集約して管理するサーバ（電子カルテサーバ、ＣＴ画像データサーバなど）であってもよい。
【００８６】
提示装置７は、ユーザに画像解析の結果を提示する機能を有する装置であればよく、例えば、提示装置７はディスプレイを備える表示装置である。あるいは、提示装置７は、医療関係者などが携帯するタブレット端末などの通信端末機器であってもよい。
【００８７】
組織を撮像した撮像画像の画像データは、外部機器８から画像解析装置１ａへと送信される。画像データを受信した画像解析装置１ａは該画像を解析し、予測スコア決定部４４によって決定された予測スコアを、通信部６を介して提示装置７に送信する。
【００８８】
画像解析システム１００ａは、外部機器８ａ、画像解析装置１ａ、および提示装置７を有している。外部機器８ａおよび画像解析装置１ａは、インターネットなどの情報通信ネットワーク５０に接続されており、相互にデータの送受信が可能である。提示装置７は、外部機器８ａに接続されている。
【００８９】
すなわち、画像解析装置１ａは、遠隔地で撮像された画像を外部機器８または８ａから受信して、画像解析を行い、解析結果である予測スコアを、提示装置７に送信することができる。なお、予測スコアを、解析に供したＣＴ画像と対応付けて提示装置７に送信してもよい。提示装置７は、外部機器８ａに接続された装置であってもよいし、画像解析装置１ａおよび外部機器８と独立している装置であってもよい。
【００９０】
この構成を採用することにより、画像解析システム１００および１００ａは、遠隔地の外部機器８および８ａからの撮像画像を受信して、当該画像を解析することができ、かつ撮像画像および決定された情報を遠隔地のユーザに提示することができる。よって、医師が不在であったり、不足していたりする遠隔地のユーザに対しても、高い精度の画像診断の結果を提供することができる。
【００９１】
〔ソフトウェアによる実現例〕
画像解析装置１、１ａの制御ブロック（特に制御部４）は、集積回路（ＩＣチップ）等に形成された論理回路（ハードウェア）によって実現してもよいし、ソフトウェアによって実現してもよい。
【００９２】
後者の場合、画像解析装置１、１ａは、各機能を実現するソフトウェアであるプログラムの命令を実行するコンピュータを備えている。このコンピュータは、例えば１つ以上のプロセッサを備えていると共に、上記プログラムを記憶したコンピュータ読み取り可能な記録媒体を備えている。そして、上記コンピュータにおいて、上記プロセッサが上記プログラムを上記記録媒体から読み取って実行することにより、本発明の目的が達成される。上記プロセッサとしては、例えばＣＰＵ（Central Processing Unit）を用いることができる。上記記録媒体としては、「一時的でない有形の媒体」、例えば、ＲＯＭ（Read Only Memory）等の他、テープ、ディスク、カード、半導体メモリ、プログラマブルな論理回路などを用いることができる。また、上記プログラムを展開するＲＡＭ（Random Access Memory）などをさらに備えていてもよい。また、上記プログラムは、該プログラムを伝送可能な任意の伝送媒体（通信ネットワークや放送波等）を介して上記コンピュータに供給されてもよい。なお、本発明の一態様は、上記プログラムが電子的な伝送によって具現化された、搬送波に埋め込まれたデータ信号の形態でも実現され得る。
【００９３】
本発明は上述した各実施形態に限定されるものではなく、請求項に示した範囲で種々の変更が可能であり、異なる実施形態にそれぞれ開示された技術的手段を適宜組み合わせて得られる実施形態についても本発明の技術的範囲に含まれる。
【実施例１】
【００９４】
本発明の一実施例について、図１０を用いて説明する。図１０は、本発明の一態様に係る画像解析方法を胸部ＣＴ画像に適用して、肺癌に罹患している患者を評価した場合の精度を、他の指標を適用した予測方法の精度と比較した表である。【００９５】
図１０において比較されている予測方法は、（１）ブリンクマン指標を用いた方法、（２）肺領域、（３）ＬＡＡ％、および（４）本発明に係る画像解析方法の４種類である。ブリンクマン指標は、喫煙指標とも呼ばれ、患者の喫煙履歴等を数値化した指標である。【００９６】
なお、ＬＡＡ％は、二値化の基準値をＣＴ値＝−９５０ＨＵとした場合、ＣＴ値＝−９１０ＨＵとした場合、およびＣＴ＝−８８０ＨＵとした場合の３つの場合について調べた。【００９７】
また、本発明に係る画像解析方法に関しては、二値化の基準値をＣＴ値＝−８８０ＨＵと設定して生成された二値化画像における０次元ベッチ数、および１次元ベッチ数を指標として用いた場合について調べた。このように、本発明に係る画像解析方法は、二値化の基準値に応じた０次元ベッチ数の変化のパターン、および１次元ベッチ数の変化のパターンではなく、所定の二値化の基準値にて生成された二値化画像における０次元ベッチ数、および１次元ベッチ数の値を指標として採用してもよい。【００９８】
図１０に示すように、肺癌の有無予測において、０次元ベッチ数、および１次元ベッチ数を指標として用いた予測は他の予測方法に比べて、精度が高く、最も信頼性の高い予測方法であることが示された。
【実施例２】
【００９９】
図１０では、単一の指標を用いて肺癌の有無を予測した場合の結果を示した。しかし、肺癌のリスク評価を行う場合、年齢・性別・ブリンクマン指数などの各指標（図１１では「予測変数」と表記されている）を組み合わせて評価することが有効であるとされている（非特許文献１参照）。【０１００】
そこで、複数の指標を組み合わせて肺癌のリスクを評価した場合の結果を図１１に示す。図１１は、（１）年齢・性別・喫煙履歴（ブリンクマン指数）という基本の指標を用いた予測モデル、（２）基本の指標に加えてＬＡＡ％を指標として用いた予測モデル、および（３）基本の指標に加えてベッチ数を指標として用いた予測モデルにおける、肺癌の予測精度を比較した結果を示す表である。【０１０１】
図１１に示すように、年齢・性別・ブリンクマン指数に加えて、ＬＡＡ％を指標として用いた予測モデルは、肺癌のリスク評価の精度の向上は見られなかった。これに対し、年齢・性別・ブリンクマン指数に加えて、ベッチ数（ここでは、二値化の基準値をＣＴ値＝−８８０ＨＵと設定して生成された二値化画像における１次元ベッチ数を用いた）を指標として用いた予測モデルでは、肺癌のリスク評価の精度が有意に向上した。図１１に示す例では、上記（１）〜（３）の各予測モデルの肺癌のリスク評価の精度を比較するために、赤池情報量基準（ＡＩＣ）を用いた。ここで、ＡＩＣはその値が低いほど精度の良い予測モデルであるとされる。また、ＡＩＣの差が１．５〜２程度あれば、統計的に有意な差であると扱われる。基本の指標のみを用いた予測モデルのＡＩＣは７８５．９であった。基本の指標に加えてＬＡＡ％を指標として用いた予測モデルのＡＩＣは７８６．２であり、ＬＡＡ％を指標に追加した効果は見られなかった。これに対して、基本の指標に加えて１次元ベッチ数ｂ１を指標として用いた場合のＡＩＣは７７９．８であり、基本の指標のみを用いた場合に比べて６．１低かった。すなわち、年齢・性別・ブリンクマン指数に加えて、１次元ベッチ数ｂ１を指標として用いた予測モデルは、比較した他の予測モデルに比べてＡＩＣが低く、肺癌のリスク評価の予測モデルとして、最も精度の良い予測モデルであることが示された。【０１０２】
一般に肺癌の大多数では、肺のＣＴ画像上に原発巣となる肺結節が確認される。本実施例では、肺癌の候補となるような肺結節の有無に依存することなく、肺癌のリスク評価をすることが可能であることが示された。すなわち、本発明の一実施例は、肺結節の有無によらず実施されたこれが肺癌リスクの層別化に有用と期待される。【０１０３】
〔まとめ〕
本発明の一態様に係る画像解析方法は、組織を撮像した撮像画像を解析する方法であって、前記撮像画像から、二値化の基準値が異なる複数の二値化画像を生成する二値化ステップ（Ｓ３２）と、複数の前記二値化画像のそれぞれについて、二値化後の一方の画素値の画素が連結して成る連結領域の数を示す第１特徴数を算出する算出ステップ（Ｓ３３）と、前記二値化の基準値と前記第１特徴数との関係性を示す第１パターンを、所定の基準パターンと比較した結果に基づいて、前記組織に生じた変化に関する情報を決定する決定ステップ（Ｓ３５）と、を含む。【０１０４】
上記の構成によれば、被検体を撮像した撮像画像に対して、二値化の基準値が異なる複数の二値化画像を生成し、二値化後の一方の画素値の画素が連結して成る連結領域の数と二値化の基準値との関係性を示す第１パターンと所定の基準パターンとの比較に基づいて、組織に生じた変化に関する情報を決定する。【０１０５】
「連結領域」とは、二次元平面である撮像画像における連結成分である。連結領域の数を算出する処理としては、上記複数の二値化画像の各々について第１特徴数として０次元ベッチ数を求める処理を行ってもよい。なお、画像からベッチ数を算出するためのプログラムは公知であり、そのプログラムを実行する装置を算出部として用いればよい。【０１０６】
また「組織」とは、１または複数種類の構成要素が一定のパターンで集合した構造単位を意味しており、細胞などの集まりから構成される生体組織である臓器などが含まれる。【０１０７】
組織を撮像した撮像画像を用いた検査・診断では、組織に生じた変化が、正常な組織とは異なる所見の有無およびその程度が利用される。画像診断を自動化する場合、組織に生じた変化の程度を評価するための手掛かりとなる特徴を如何に正確に検出して定量化するか、が重要である。【０１０８】
医師などが肺気腫などの疾患の有無および重症度を、肺を撮像した撮像画像（例えば、ＣＴ画像）に基づいて決定する場合、例えば、以下の（１）および（２）のような特徴が手掛かりとされる。（１）軽度の肺気腫では、直径約１ｃｍ程度の低吸収領域の散在が認められる。（２）重度の肺気腫では、低吸収領域は融合して増加する一方、正常な肺の領域は減少する。【０１０９】
本発明の発明者は、上記の方法を適用することによって、（１）および（２）のような組織に生じた変化が示す特徴を組織画像から検出・定量化し、該変化の程度を高精度に評価できる可能性があることを見出した。【０１１０】
上記の方法を採用することにより、組織に生じた変化の特性に合致する特徴に基づいて、該変化の有無および程度に関する情報を精度良く決定することができる。例えば、決定した情報を医師に提示することによって、組織に生じた変化の有無および程度の判断を適切にサポートすることができる。【０１１１】
前記二値化ステップにて生成した複数の二値化画像のそれぞれについて、前記一方の画素値の画素に囲まれた、他方の画素値の画素からなる穴形状の領域の数を示す第２特徴数を算出し、前記決定ステップにおいて、前記二値化の基準値の変化に応じて前記第２特徴数が変化する第２パターンを、所定の基準パターンと比較した結果に基づいて、前記組織に生じた変化に関する情報を決定してもよい。【０１１２】
組織に生じた変化の特性は様々であり、撮像画像に表れる特徴も様々である。それゆえ、組織に生じた変化を、撮像画像中の連結領域の数の変化として検出するよりも、穴形状の数の変化として検出した方が適切である場合がある。【０１１３】
上記の構成によれば、被検体を撮像した撮像画像に対して、二値化の基準値が異なる複数の二値化画像を生成し、二値化後の一方の画素値の画素に囲まれた穴形状の領域の数の変化パターンと所定の基準パターンとの比較に基づいて、組織に生じた変化に関する情報を決定する。これにより、二次元平面である撮像画像において、組織に生じたさまざまな変化を高精度に検出・定量化することができる。よって、多様な組織において生じた変化を撮像した撮像画像に適用することができる。【０１１４】
「穴形状の領域」とは、二次元平面である撮像画像における空間（穴）である。穴形状の領域の数を算出する処理としては、上記複数の二値化画像の各々について第２特徴数として１次元ベッチ数を求める処理を行ってもよい。なお、画像からベッチ数を算出するためのプログラムは前述のように公知であり、そのプログラムを実行する装置を算出部として用いればよい。【０１１５】
前記組織に生じた変化に関する確定情報が予め対応付けられた前記撮像画像を取得して、該撮像画像の前記第１パターンと、前記確定情報との組み合わせを含む学習データを用いた機械学習によって、前記基準パターンを生成する学習ステップ（Ｓ１）をさらに含んでいてもよい。【０１１６】
「確定情報」とは、医師によって予め評価・決定された、組織に生じた変化に関する情報であり、例えば、組織に生じた疾患の有無および重症度などを示す情報である。【０１１７】
上記の構成によれば、組織に関する確定情報が予め対応付けられた撮像画像の第１パターンと確定情報との相関関係について学習する。これにより、学習された相関関係に倣って、解析対象となる撮像画像の第１パターンから組織に生じた変化に関する情報を決定することができる。【０１１８】
前記確定情報は、前記組織に生じた変化の程度に応じた複数の段階に関する情報を含み、前記基準パターンは、前記段階ごとに生成されており、前記決定ステップにおいて、前記組織に生じた変化が前記段階のいずれに該当するかを決定してもよい。【０１１９】
これにより、学習された相関関係に倣って、解析対象となる撮像画像の第１パターンから、組織に生じた変化の程度を決定することができる。【０１２０】
前記撮像画像から解析対象となる組織に対応する領域を抽出する抽出ステップをさらに含んでいてもよい。【０１２１】
組織を撮像した画像には、解析対象ではない部位に対応する領域も含まれている場合がある。上記のように構成すれば、解析対象の組織に対応していない不要な領域であり、組織に生じた変化の検出・定量化を阻害する可能性がある領域を抽出することができる。【０１２２】
本発明の一態様に係る画像解析装置（１、１ａ）は、組織を撮像した撮像画像を解析する画像解析装置であって、前記撮像画像から、二値化の基準値が異なる複数の二値化画像を生成する二値化部（４１）と、複数の二値化画像のそれぞれについて、二値化後の一方の画素値の画素が連結して成る連結領域の数を示す第１特徴数を算出する算出部（ベッチ数算出部４２）と、前記二値化の基準値と前記第１特徴数との関係性を示す第１パターンを、所定の基準パターンと比較する比較部（パターン特定部４３）と、前記比較の結果に基づいて、前記組織に生じた変化に関する情報を決定する決定部（予測スコア決定部４４）と、を備える。【０１２３】
上記の構成によれば、組織に生じた変化の特性に合致する特徴に基づいて、該変化の有無および程度に関する情報を精度良く決定することができる、という上記方法と同様の効果を奏する。【０１２４】
前記算出部は、前記二値化部によって生成された複数の二値化画像のそれぞれについて、前記一方の画素値の画素に囲まれた、他方の画素値の画素からなる穴形状の領域の数を示す第２特徴数を算出し、前記比較部は、前記二値化の基準値の変化に応じて前記第２特徴数が変化する第２パターンを、所定の基準パターンと比較し、前記決定部は、比較結果に基づいて、前記組織に生じた変化に関する情報を決定する構成であってもよい。【０１２５】
上記の構成によれば、二次元平面である撮像画像において、組織に生じたさまざまな変化を高精度に検出・定量化することができる。よって、多様な組織を撮像した撮像画像に適用することができる。【０１２６】
前記組織に生じた変化に関する確定情報が予め対応付けられた前記撮像画像を取得して、該撮像画像の前記第１パターンと、前記確定情報との組み合わせを含む学習データを用いた機械学習によって、前記比較部に前記基準パターンを生成させる学習制御部（４６）をさらに備えていてもよい。【０１２７】
上記の構成によれば、組織に関する確定情報が予め対応付けられた撮像画像の第１パターンと確定情報との相関関係について学習する。これにより、学習された相関関係に倣って、解析対象となる撮像画像の第１パターンから組織に生じた変化に関する情報を決定することができる。【０１２８】
前記確定情報は、前記組織に生じた変化の程度を複数の段階に区分する情報を含み、前記基準パターンは、前記段階ごとに生成されており、前記決定部は、前記組織に生じた変化が前記段階のいずれに該当するかを決定する構成であってもよい。【０１２９】
これにより、解析対象となる撮像画像の第１パターンから、組織に生じた変化の程度を決定することができる。【０１３０】
前記撮像画像から解析対象となる組織に対応する領域を抽出する領域抽出部（二値化部４１）をさらに備える構成であってもよい。【０１３１】
上記の構成によれば、解析対象の組織に対応していない不要な領域であり、組織に生じた変化の検出・定量化を阻害する可能性がある領域を抽出することができる。【０１３２】
本発明の一態様に係る画像解析システム（１００、１００ａ）は、前記の画像解析装置（１、１ａ）と、前記撮像画像の画像データを前記画像解析装置へ送信する外部機器（８、８ａ）と、前記画像解析装置において決定された情報を取得して該情報を提示する提示装置（７）と、を含む。【０１３３】
これにより、例えば遠隔地の外部機器を用いて撮像された撮像画像を受信して、当該画像を解析することができ、かつ撮像画像および決定された情報を遠隔地のユーザに提示することができる。【０１３４】
本発明の各態様に係る画像解析装置は、コンピュータによって実現してもよく、この場合には、コンピュータを上記画像解析装置が備える各部（ソフトウェア要素）として動作させることにより上記画像解析装置をコンピュータにて実現させる画像解析プログラム、およびそれを記録したコンピュータ読み取り可能な記録媒体も、本発明の範疇に入る。
【符号の説明】
【０１３５】
１、１ａ  画像解析装置
２  画像取得部
４  制御部
７  提示装置
８、８ａ  外部機器
４１  二値化部（二値化部、領域抽出部）
４２  ベッチ数算出部（算出部）
４３  パターン特定部（比較部）
４４  予測スコア決定部（決定部）
４６  学習制御部
１００、１００ａ  画像解析システム
Ｓ１  学習ステップ
Ｓ３２  二値化ステップ
Ｓ３３  算出ステップ
Ｓ３５  決定ステップ
【図１】
2019102829.tif 000003
【図２】
2019102829.tif 000004
【図３】
2019102829.tif 000005
【図４】
2019102829.tif 000006
【図５】
2019102829.tif 000007
【図６】
2019102829.tif 000008
【図７】
2019102829.tif 000009
【図８】
2019102829.tif 000010
【図９】
2019102829.tif 000011
【図１０】
2019102829.tif 000012
【図１１】
2019102829.tif 000013
【国際調査報告】
2019102829.tif 000014
2019102829.tif 000015
2019102829.tif 000016
2019102829.tif 000017