2.基礎的な画像特徴量の計算

【概要】


画像の特徴を表す量として,「特徴量」がさまざまに定義できる.

特徴量には,画素値の分布に基づくもの(演習2-1)と,対象の形に基づくもの(演習2-2,2-3)などがある.

病変の画素値や形などの特徴量を数多く求めて診断に役立てる手法はRadiomics(ラジオミクス)と呼ばれ,近年研究が進んでいる.


【画素値に基づく特徴量】


画素値に基づく基礎的な特徴量として,画素値の平均(Mean),分散(Variance),ヒストグラムなどの統計的な特徴量を求めることが多い.


平均は,画像全体あるいは一定領域内の画素値の和を,画素数で割ることによって求めることができる.

分散は,同じ範囲の画素値の2乗平均と平均の2乗の差で求めることができる.

これは,画素値のバラツキの指標となる.この指標は,分散の平方根である標準偏差(Standard Deviation)で示すことも多い.

例えば,CT画像の画質評価では,均一なファントムに設定したROI内の平均値(CT値)と標準偏差(ノイズ)が用いられる.


ヒストグラムは,横軸に画素値,縦軸にその画素値の頻度をとった分布グラフである.


 演習2-1:基礎的な特徴量として,平均,標準偏差,ヒストグラムを求めてみよう!

使用する画像: 胸部X線画像, 同心円チャート画像, シーメンスチャート画像

*3つの画像はいずれも反転LUT(inverting LUT)で表示される(第1回 演習1-1を参照).特徴量やヒストグラムは表示の明暗ではなく画素値から計算されるため,暗く表示されている部分ほど画素値は大きい.


方法:

1.ImageJを起動後,対象画像を開く.

⇒「File」→「Open」を押し,対象画像を選択.


2.計測項目を設定する.

⇒「Analyze」→「Set Measurements」を選択.

「Area」「Mean gray value」「Standard deviation」「Min & max gray value」「Display label」にチェックを入れる.
*「Display label」にチェックを入れると,Resultsに画像名が表示される.


3.特徴量を測定し,ヒストグラムを表示する.

⇒「Analyze」→「Measure」を選択.手順2で設定した特徴量が「Results」ウィンドウに表示される.

⇒「Analyze」→「Histogram」を選択.ヒストグラムが表示される.

*領域を選択していないときは,画像全体が対象となる.


4.領域を選択して,特徴量の相違を比較する.

⇒矩形選択ツールで任意の領域を選択し,「Edit」→「Selection」→「Specify」でWidthとHeightを指定したうえで,手順3を実施する.


考察:


・3つの画像のヒストグラムの形,平均値,標準偏差を比較し,それぞれの画像の濃淡の特徴(どのような画素値がどの程度含まれているか)がどのように数値に表れているかを説明してみよう.


・胸部X線画像の肺野内と肺野外(縦隔など)に同じ大きさの領域を設定し,特徴量の違いを,画像の見え方(明るさ,濃淡のばらつき)と関連づけて説明してみよう.



【形状特徴量】


画像中の対象(病変や臓器など)の特徴は,画素値だけでなく,その「形」によっても表すことができる.

例えば,肺結節では,良性のものは円形で辺縁が滑らかなことが多く,悪性のものは分葉状や棘状(スピキュラ)の不整な形を示すことが多い.

形を数値で表したものを形状特徴量と呼び,ImageJでは次のような特徴量を求めることができる.


・面積(Area),周囲長(Perimeter)

・円形度(Circularity)= 4π×面積/周囲長2 (真円で1となり,細長い形や凹凸の多い形ほど0に近づく)

・凸性(Solidity)= 面積/凸包(対象を囲む最小の凸多角形)の面積 (凹みがあるほど小さくなる)


 演習2-2:図形の形状特徴量を求め,形の違いを数値で比較してみよう!

使用する画像: 図形画像(shapes.tif,120×600画素,8ビット2値画像,Tif形式)

*上から円,楕円,四角形,三角形,棘状(スピキュラ)の5つの図形が並んでいる.


方法:

1.ImageJを起動後,対象画像を開く.

⇒「File」→「Open」を押し,対象画像を選択.


2.2値画像の設定を確認する.

⇒「Process」→「Binary」→「Options」を選択し,「Black background」のチェックを外してOKを押す.

*この画像は図形が黒,背景が白である.チェックが入っていると,背景が測定対象になってしまう.


3.計測項目を設定する.

⇒「Analyze」→「Set Measurements」を選択.

「Area」「Perimeter」「Shape descriptors」「Display label」にチェックを入れる.
*「Shape descriptors」を選ぶと,円形度(Circ.)と凸性(Solidity)のほかにAR,Roundも表示されるが,今回は主に円形度と凸性に注目する.


4.図形ごとの特徴量を一括で測定する.

⇒「Analyze」→「Analyze Particles」を選択し,「Display results」にチェックを入れ,「Show」を「Outlines」にしてOKを押す.

*5個の図形の特徴量が「Results」ウィンドウに表示される.Resultsの番号は,Outlines画像に表示される番号と対応している.
番号1が円,2が楕円,3が四角形,4が三角形,5が棘状の図形である.


考察:


・円,楕円,四角形,三角形,棘状の図形で,円形度と凸性の値はどのように異なったか.2つの特徴量がそれぞれ形のどのような違いを表しているかを,自分の結果に基づいて説明してみよう.


・円形度と凸性だけでは区別しにくい図形の組はあったか.その図形を区別するには,どのような特徴量が役立つか,Resultsに表示された他の特徴量も参考に考えてみよう.



 演習2-3:肺結節の形状特徴量を求め,良性と悪性で比較してみよう!

使用する画像:

良性結節:CT画像(nodule_benign.tif),病変マスク(nodule_benign_mask.tif)

悪性結節:CT画像(nodule_malignant.tif),病変マスク(nodule_malignant_mask.tif)

*いずれも胸部CT画像の1スライスから結節の周囲(96×96画素)を切り出したものである.
*良性結節は2年間の経過観察で変化がないことにより良性と確定したもの,悪性結節は生検により原発性肺癌と確定したものである.
*病変マスクは,放射線科医1名が描いた結節の輪郭から作成した2値画像である(結節が黒,背景が白).
*画像には画素の寸法(mm)が設定されているため,面積はmm2で表示される.


方法:

1.ImageJを起動後,良性結節のCT画像と病変マスクを開く.

⇒「File」→「Open」を押し,対象画像を選択.CT画像で結節の形を確認する.


2.計測項目を確認する.

⇒演習2-2と同じ設定(「Black background」のチェックなし,「Area」「Perimeter」「Shape descriptors」「Display label」にチェック)であることを確認する.


3.病変マスクの形状特徴量を測定する.

⇒病変マスクの画像を選択し,「Analyze」→「Analyze Particles」を選択.「Display results」にチェックを入れてOKを押す.


4.悪性結節についても,手順1~3を行う.

*「Results」ウィンドウには,測定結果が1行ずつ追加される.「Label」の列に表示される画像名で,良性と悪性の結果を見分けられる.


考察:


・良性と悪性の結節で,円形度と凸性の値はどのように異なったか.演習2-2の図形の結果と比較し,それぞれの結節の形の特徴と関連づけて説明してみよう.


・CT画像で見た結節の形と,病変マスクの形を見比べてみよう.悪性結節の周囲に見られる細い線状の影は,マスクにどの程度含まれているか.それが円形度や凸性の値にどのように影響していると考えられるか.



[演習2-3で使用した画像の出典]
The Lung Image Database Consortium(LIDC)and Image Database Resource Initiative(IDRI)のデータ(症例 LIDC-IDRI-0234,LIDC-IDRI-0265)を使用した.ライセンス:CC BY 3.0.
改変内容:CT画像の1スライスから結節の周囲を切り出してCT値に変換し,TIFF形式で保存した.病変マスクは,同データに含まれる放射線科医の輪郭から作成した.
・Armato SG III, et al. Data From LIDC-IDRI [Data set]. The Cancer Imaging Archive, 2015. https://doi.org/10.7937/K9/TCIA.2015.LO9QL9SX
・Armato SG III, et al. The Lung Image Database Consortium (LIDC) and Image Database Resource Initiative (IDRI): A completed reference database of lung nodules on CT scans. Medical Physics, 38(2): 915-931, 2011.
・Clark K, et al. The Cancer Imaging Archive (TCIA): Maintaining and Operating a Public Information Repository. Journal of Digital Imaging, 26(6): 1045-1057, 2013.


 実験項目へ戻る