JP2020185223A

JP2020185223A - 単色ｃｔ画像作成方法、単色ｃｔ画像作成装置および単色ｃｔ画像作成プログラム

Info

Publication number: JP2020185223A
Application number: JP2019092373A
Authority: JP
Inventors: 創大田中; Sodai Tanaka; 博樹白土; Hiroki Shirato
Original assignee: Hokkaido University NUC
Current assignee: Hokkaido University NUC
Priority date: 2019-05-15
Filing date: 2019-05-15
Publication date: 2020-11-19
Also published as: US20200359980A1; US11324466B2

Abstract

【課題】適切な単色ＣＴ画像を短時間で自動的に得る。【解決手段】画像処理部２２は、複数のエネルギーレベルの放射線を照射して得た複数のＣＴ画像について演算処理して、単色ＣＴ画像を得る。第１エネルギーレベルの放射線を照射して得た第１エネルギーレベルＣＴ画像と、第２エネルギーレベルの放射線を照射して得た第２エネルギーレベルＣＴ画像を取得し、取得した第１および第２エネルギーレベルＣＴ画像について複数種類の重み付け演算し、それぞれの重み付け演算結果における複数の単色ＣＴ画像を算出し、複数の単色ＣＴ画像のそれぞれについて、高吸収体の周辺領域を周方向において所定の面積ごとに複数の注目領域に分割し、各注目領域の画像データの平均値を用い周辺領域についての標準偏差を算出し、標準偏差が小さい単色ＣＴ画像を選択する。【選択図】図１

Description

本発明は、複数のエネルギーレベルの放射線を照射して得た複数のＣＴ画像について演算処理し単色ＣＴ画像を得る、単色ＣＴ画像の作成に関する。

従来、放射線診断・治療などにおいて、コンピュータ断層撮影（ＣＴ：Computed Tomography）装置が広く利用されている。ＣＴ装置では、対象（患者）にＸ線を照射して、その投影画像から断層画像（ＣＴ画像）を作成する。

ここで、ＣＴ画像においては、各種のアーチファクトにより画像が劣化することが知られている。例えば、生体中に金属などの高吸収体が存在すると、Ｘ線のエネルギースペクトルが変化し、または検出器におけるＸ線強度が十分でなくなり、正しい投影画像が得られなくなる、いわゆる金属（高吸収体）アーチファクトが発生する。この金属アーチファクトは、２種類の管電圧で得たＣＴ画像について重み付け演算を行って仮想単色ＣＴ画像（以下、単に単色ＣＴ画像という）を得ることで除去することができる。そこで、２つのＸ線ＣＴ画像から単色ＣＴ画像を得る機能が、商用Ｘ線ＣＴ装置に導入されてきている（特許文献１参照）。

ＷＯ２０１４／０３４６１８号公報

ここで、単色ＣＴ画像を得るためには、２つのＣＴ画像を適切な重みで重み付け演算する必要がある。現状のＸ線ＣＴ装置では、単色ＣＴ画像を算出する際の重みは手動で設定する。すなわち、重みを変更しながら、操作者が得られた単色ＣＴ画像を視覚的に判断し、適切な重みを選択する。このため、操作者の判断が必要であるとともに、その調整に時間がかかるという問題もあった。

本発明は、適切な単色ＣＴ画像を短時間で自動的に得ることを目的とする。

本発明は、複数のエネルギーレベルの放射線を照射して得た複数のＣＴ画像について演算処理して単色ＣＴ画像を得る単色ＣＴ画像作成方法であって、第１エネルギーレベルの放射線を照射して第１エネルギーレベルＣＴ画像を得、第２エネルギーレベルの放射線を照射して第２エネルギーレベルＣＴ画像を得、得られた第１および第２エネルギーレベルＣＴ画像について複数種類の重み付け演算し、それぞれの重み付け演算結果における複数の単色ＣＴ画像を得、複数の単色ＣＴ画像のそれぞれについて、高吸収体領域の周辺領域を周方向において所定の面積ごとに複数の注目領域に分割し、各注目領域の画像データの平均値を用い周辺領域についての標準偏差を算出し、標準偏差が小さい単色ＣＴ画像を選択する。

第１および第２エネルギーレベルＣＴ画像について複数種類の重み付け演算は、重み付け加算であるとよい。

高吸収体の周辺領域は、放射線を照射して得た複数のＣＴ画像において、高吸収体アーチファクトの影響が出る領域に基づいて決定するとよい。

高吸収体の周辺領域は、放射線を照射して得た複数のＣＴ画像において、高吸収体アーチファクトが存在すると推定される領域の周辺における近接領域を除外した領域であるとよい。

本発明は、複数のエネルギーレベルの放射線を照射して得た複数のＣＴ画像について演算処理して、単色ＣＴ画像を得る単色ＣＴ画像作成装置であって、第１エネルギーレベルの放射線を照射して得た第１エネルギーレベルＣＴ画像と、第２エネルギーレベルの放射線を照射して得た第２エネルギーレベルＣＴ画像を取得し、取得した第１および第２エネルギーレベルＣＴ画像について複数種類の重み付け演算し、それぞれの重み付け演算結果における複数の単色ＣＴ画像を算出し、複数の単色ＣＴ画像のそれぞれについて、高吸収体の周辺領域を周方向において所定の面積ごとに複数の注目領域に分割し、各注目領域の画像データの平均値を用い周辺領域についての標準偏差を算出し、標準偏差が小さい単色ＣＴ画像を選択する。

また、本発明は、複数のエネルギーレベルの放射線を照射して得た複数のＣＴ画像について演算処理して、単色ＣＴ画像を得る単色ＣＴ画像作成プログラムであって、単色作成プログラムは、コンピュータに、第１エネルギーレベルの放射線を照射して得た第１エネルギーレベルＣＴ画像と、第２エネルギーレベルの放射線を照射して得た第２エネルギーレベルＣＴ画像を取得させ、取得した第１および第２エネルギーレベルＣＴ画像について複数種類の重み付け演算し、それぞれの重み付け演算結果における複数の単色ＣＴ画像を算出させ、複数の単色ＣＴ画像のそれぞれについて、高吸収体の周辺領域を周方向において所定の面積ごとに複数の注目領域に分割し、各注目領域の画像データの平均値についての標準偏差を算出し、標準偏差が小さい単色ＣＴ画像を選択させる。

本発明によれば、適切な重みを短時間で自動的に決定でき、高吸収体アーチファクトを除去した単色ＣＴ画像を効果的に得ることができる。

本実施形態に係るＣＴ画像作成装置を搭載するＣＴ装置１００の全体構成を示すブロック図である。金属領域Ｍｅｔａｌの周辺を示す図であり、（ａ）は、注目領域ＲＯＩ、（ｂ）は周辺注目領域ＳＲＯＩを示す。最適な単色ＣＴ画像を選択する手順を示すフローチャートである。ＣＴ画像の例を示す図であり、（ａ）は管電圧８０−ｋＶでの撮影ＣＴ画像、（ｂ）は管電圧１１０−ｋＶでの撮影ＣＴ画像、（ａ）は管電圧１４０−ｋＶでの撮影ＣＴ画像、（ａ）は最適化した単色ＣＴ画像、（ａ）はホワイトノイズを付加した場合の最適化した単色ＣＴ画像を示す。注目領域ＲＯＩの他の例を示す図である。

以下、本発明の実施形態について、図面に基づいて説明する。なお、本発明は、ここに記載される実施形態に限定されるものではない。

図１は、本実施形態に係るＣＴ画像作成装置を搭載するＣＴ装置１００の全体構成を示すブロック図である。

例えば、診断、治療の対象となる患者などの対象物１０は、診察台（図示せず）に固定される。対象物１０の上方には線源１２が配置される。線源１２は、放射線を対象物１０に照射するものであり、この例ではＸ線管である。なお、線源１２としては、Ｘ線の他、γ線などの高エネルギー電磁波を発生するものも採用される。

線源１２には、電圧発生部１４が接続されており、線源１２において放射線を発生するために必要な高電圧が発生され、これが線源１２に供給される。この例では、電圧発生部１４から供給される電圧により線源１２のＸ線管の管電圧が変更され、放射されるＸ線のエネルギーレベルが変更される。また、対象物１０について所定範囲の断層画像を得るために、対象物１０を固定する診察台や線源１２は移動（回転および軸方向の移動）できるようになっている。移動制御部１６は、これらの移動を制御する。

また、対象物１０の下方（線源１２と反対側）には、検出器１８が配置されており、対象物１０を透過した放射線が入射することで、対象物１０の透過像が画素ごとの電気信号として得られる。

検出器１８には、信号処理部２０が接続され、検出器１８からの電気信号に基づいて撮影された２次元画像データを順次生成する。すなわち、線源１２および検出器１８は、対象物１０の回りを回転して、３６０度方向から画像データを取得する。取得された画像データは画像処理部２２に供給され、ここで各種処理がなされ、逐次変化する３次元データが構築され、これが記憶部２４に記憶される。そして、記憶部２４に記憶されたデータに基づき、所望のＣＴ画像を得ることができる。

ここで、本実施形態では、観察したい１つのＣＴ画像について、２つの管電圧（第１エネルギーレベルおよび第２エネルギーレベル）の撮影ＣＴ画像を選択し、それらを重み付け加算することによって単色ＣＴ画像を作成する。特に、重みを変更しながら複数の単色ＣＴ画像を作成し、その中から後述するようにして適切な単色ＣＴ画像を選択する。これらの処理は、記憶部２４に記憶されているデータを利用して、画像処理部２２が行うとよい。なお、少なくとも画像処理部２２はコンピュータで構成され、単色ＣＴ画像作成処理プログラムを実行することで、記憶部２４に記憶されている撮影ＣＴが像を取得しこれらの演算処理を行う。

ＣＴ装置１００は、制御部３０を有しており、電圧発生部１４、移動制御部１６、信号処理部２０、画像処理部２２などの動作を制御する。制御部３０には、入力部３２、出力部３４が接続されている。入力部３２は、キーボードやタッチパネルからなり、操作者が各種情報を入力する。出力部３４は、ディスプレイやプリンタからなり、所望の画像や、文字データを出力する。

「単色ＣＴ画像の作成」
線源１２からの放射線のエネルギーレベルは、制御部３０からの指令に応じた電圧が電圧発生部１４から出力されることによって制御される。本実施形態においては、対象物１０を移動（線源１２および検出器１８の回転軸の方向に移動）して得る撮影を２種類のエネルギーレベル（第１エネルギーレベルおよび第２エネルギーレベル）の放射線によって行う。この２種類の撮影は、対象物１０を移動しながら順次行ってもよいし、１種類のエネルギーレベルでの撮影を終了した後、再度対象物１０を移動しながら他の種類のエネルギーレベルでの撮影を行ってもよい。

２種類のエネルギーレベルでの撮影が終了した場合には、それぞれの撮影画像データが記憶部２４に記憶される。そこで、これら撮影画像データから所望の撮影ＣＴ画像を作成することができ、この際にエネルギーレベルの異なるだけで他は同一の撮影ＣＴ画像を作成する。これらは、記憶部２４に記憶される。

そして、入力部３２からの指令に基づき、１つのＣＴ画像を特定し、この特定されたＣＴ画像と、対応するエネルギーレベルの異なるＣＴ画像を記憶部２４から読み出し、画像処理部２２において重み付け加算する。

「高吸収体領域（金属領域）の周辺領域についての処理」
金属は放射線の高吸収体であり、金属の周囲に高吸収体アーチファクト（金属アーチファクト）が生じる。そこで、この金属アーチファクトを除去するために、２つのＣＴ画像の重み付け加算を行う。そして、この重みは、金属周囲を取り囲む領域を複数に分割した注目領域ＲＯＩの平均値の分散度合いである標準偏差を評価値として、標準偏差が最も小さくなるように選択する。なお、金属領域は、ＣＴ画像において容易に認識できる。

図２（ａ）は、金属領域Ｍｅｔａｌとその周辺の扇状の注目領域ＲＯＩを示す。この例では、周辺注目領域ＳＲＯＩが２４の注目領域ＲＯＩに分割されている。従って、注目領域ＲＯＩは、金属領域Ｍｅｔａｌの中心から円周方向３０度ずつの幅を有する。図２（ｂ）には、金属領域Ｍｅｔａｌとその周囲の周辺注目領域ＳＲＯＩを示す。

「標準偏差を用いた重み付け演算」
図３は、最適な単色ＣＴ画像を選択する手順を示すフローチャートである。まず、撮影ＣＴ画像の中から、金属領域Ｍｅｔａｌを抽出する（Ｓ１１）。これは、１つの撮影ＣＴ画像における画素毎のＣＴ値の変化からエッジを検出することなどによって行えばよい。なお、複数の撮影ＣＴ画像から金属領域Ｍｅｔａｌを検出してもよく、１つであればエネルギーレベル（管電圧）の高い方の画像から検出するとよい。

金属領域Ｍｅｔａｌが抽出された場合、この周縁から数画素分を近接領域として抽出し（Ｓ１２）、近接領域から外側に向けて拡大し周辺注目領域ＳＲＯＩを形成する（Ｓ１３）。

ここで、必ずしも近接領域を設定しなくてもよいが、金属領域Ｍｅｔａｌが周辺領域に入ったり、入らなかったりすると、後述する標準偏差の誤差が大きくなるため、近接領域抽出して周辺注目領域ＳＲＯＩから除外するのがよい。なお、近接領域をあまり大きくすると金属アーチファクトが存在する領域が少なくなり、評価の精度が落ちるため、周辺の放射方向数画素を近接領域に設定するとよい。また、金属領域Ｍｅｔａｌから周辺領域に拡大してから、近接領域を削除して周辺注目領域ＳＲＯＩを形成してもよい。

周辺注目領域ＳＲＯＩは、金属アーチファクトが存在する領域に対応する領域（金属アーチファクトが存在すると推定される領域）であることが好ましいが、その領域より若干小さくても、若干大きくてもよい。例えば、金属領域Ｍｅｔａｌの直径が１ｃｍ程度の場合に、その周囲数ｃｍの厚みのドーナツ状の領域を周辺注目領域ＳＲＯＩにするとよい。

このように、周辺注目領域ＳＲＯＩが設定された場合、これを円周方向において等分割して扇台形状の複数の注目領域ＲＯＩを形成する（Ｓ１４）。１つの金属領域Ｍｅｔａｌについてその周辺の分割領域である注目領域ＲＯＩが得られた場合には、重みについての最適演算処理を行う（Ｓ１５）。

Ｓ１５の最適演算では、選択した２つの撮影ＣＴ画像について、両者の和が１となるようにしながら、一方の重みを−９→１０、他方の重みを１０→−９の範囲で段階的に変化させ、それぞれの重みでの重み付け加算した単色ＣＴ画像を作成する（Ｓ１５−１）。そして、得られた重み付けの異なる単色ＣＴ画像のそれぞれについて注目領域ＲＯＩごとに画像データ（ＣＴ値）の平均値を算出し（Ｓ１５−２）、得られた各ＲＯＩの平均値についての標準偏差を算出する（Ｓ１５−３）。そして、最も標準偏差の小さい単色ＣＴ画像を選択し（Ｓ１５-４）、選択された単色ＣＴ画像を出力部３４に表示する（Ｓ１６）。

ここで、各種の撮影ＣＴ画像を用い、検証したところ、このような重み付け演算結果に従って選択された、標準偏差が最小の単色ＣＴ画像が金属アーチファクトの影響を除去した最適な単色ＣＴ画像であることが確認された。

このように、本実施形態では、多数の異なる重み付け加算により得られた単色ＣＴ画像について、注目領域ＲＯＩを用いて評価をして、自動的に適切な単色ＣＴ画像を得ることができる。

「実施例」
図４は、対象物１０としてA SUN NUCLEAR COMPANY のGAMMEXファントムを採用し、General Electric CompanyのＣＴ装置Optima CT 580W を用いて得られたＣＴ画像である。

図４（ａ）は管電圧８０−ｋＶ、（ｂ）は管電圧１２０−ｋＶ、（ａ）は管電圧１４０−ｋＶでの撮影ＣＴ画像である。なお、GAMMEXファントムには、金属として、１つだけチタンを埋め込んでいる。図４（ａ）−図４（ｃ）において、輝いて見える（周辺にアーチファクト）が見える丸がチタンの画像である。なお、対象物１０に利用したGAMMEXファントムは、水を模擬した基体に骨等を模擬した挿入物を埋め込んだもので、挿入物が基体とは異なる明度になっている。この画像では、吸収が大きいものが白くなっている。

図４（ｄ）は、本実施形態による最適化演算、すなわち８０−ｋＶ（図４（ａ））の撮影ＣＴ画像と１４０−ｋＶ（図４（ｃ））の撮影ＣＴ画像を重み付け加算して得られた複数の単色ＣＴ画像の中で図２（ａ）に示された２４の注目領域ＲＯＩの平均値の標準偏差ＳＤＭＡＤが最小になるように決定された最適な単色ＣＴ画像である。また、図４（ｅ）は、処理対象となるＣＴ画像にランダムなホワイトノイズを付加した場合において同様にして得た最適な単色ＣＴ画像である。

図４（ａ）−図４（ｅ）の各ＣＴ画像における、注目領域ＲＯＩの平均値の標準偏差ＳＤＭＡＤ、周辺注目領域ＳＲＯＩの標準偏差ＳＤＳＲＯＩは、（ａ）：９３．５８８，１２８．２８６、（ｂ）：３５．３５７，５０．３６３、（ｃ）：２６．５７９，３７．６４３、（ｄ）：７．３４３，１８．６８７、（ｅ）：７．５９７，３０．９６１である。

このように、撮影ＣＴ画像では、標準偏差ＳＤＭＡＤ，ＳＤＳＲＯＩはともに管電圧が高いものほど小さくなる。また、図４（ｄ）の本実施形態の手法による標準偏差ＳＤＭＡＤを最小とした単色ＣＴ画像によれば、標準偏差ＳＤＳＲＯＩも撮影ＣＴ画像よりも十分小さくなっており、金属アーチファクトが除去されていることがわかる。

そして、ホワイトノイズを付加した図４（ｄ）における標準偏差ＳＤＭＡＤ，ＳＤＲＯＩをみると、標準偏差ＳＤＳＲＯＩは、ホワイトノイズの影響で６５％程度上昇しているが、標準偏差ＳＤＭＡＤは３−４％しか上昇していない。これより、本実施形態で利用する、注目領域ＲＯＩの平均値を用いる標準偏差ＳＤＭＡＤは、ホワイトノイズの影響が小さく、ホワイトノイズがあっても適切な重みを算出できることがわかる。すなわち、標準偏差ＳＤＭＡＤが最小の単色ＣＴ画像を選択することによって、ホワイトノイズの影響に依存せず、金属アーチファクトを除去した画像を選択できる。

このように、本実施形態によれば、２つの撮影ＣＴ画像から適切な重み付け加算を選択して、金属アーチファクトの影響が抑制された単色ＣＴ画像を自動的に得ることができ、これに要する時間は３秒程度である。従来は、手動で最適な画像を選ぶ必要があり、少なくとも１５秒程度は必要であり、５倍程度の時間短縮が可能である。なお、本実施形態の方法は、高性能のコンピュータを使用することで、さらに時間の短縮は可能であり、操作者は実質的に待つことなく即時に適切な単色ＣＴ画像をみることができるといえる。

本実施形態によれば、放射線診断・治療において、金属アーチファクトの影響が最小化されるような画像が即座に自動的に提示される。このため、アダプティブ放射線治療（ART：Adaptive Radiotherapy、即時適合型放射線治療）に特に有効と考えられる。すなわち、治療時点での適切なＣＴ画像が即時に得られるため、治療期間中の体形変化や病巣の大きさの変化、臓器の位置変化などにより放射線の照射が初期の計画通りにならない場合でも、その時のＣＴ画像を利用して治療期間の途中から再計画を行い、より適切な治療を続けることができる。

「注目領域への分割」
上述の例では、金属領域Ｍｅｔａｌについて円形とした。しかし、金属は円形には限らない。このような場合には、例えば金属領域の接線に垂直な方向に一定距離伸ばして、金属領域の周辺に沿って周辺注目領域ＳＲＯＩを設定し、得らえた周辺注目領域ＳＲＯＩを注目領域ＲＯＩに分割すればよい。例えば、周辺注目領域ＳＲＯＩの面積を所望の分割数で除算してＲＯＩ面積を設定し、任意の基準位置（放射方向の線）から決定したＲＯＩ面積毎に個別の注目領域ＲＯＩを設定すればよい。

この場合、金属領域Ｍｅｔａｌの重心を基準として、幅（中心角）を広げていき、設定されたＲＯＩ面積毎に注目領域ＲＯＩを設定することができる。また、図５に示すように、金属の周縁に対し直角な方向に１つの境界線を設定し、ＲＯＩ面積毎にもう一方の境界線を設定してもよい。

「最適化演算」
上述の例では、重み付け演算において、重みを所定の間隔で分割し、それぞれの注目領域ＲＯＩの平均値に基づいた標準偏差ＳＤＭＡＤを計算して、標準偏差ＳＤＭＡＤが最小になる重み付けの単色ＣＴ画像を選択する最適演算を行った。しかし、１度の演算ではなく、分割間隔を順次狭くする方法を採用してもよい。例えば、最も標準偏差の小さい方から２つを選択し、その２つの間をさらに分割し、最小の標準偏差の重みを決定するなどの手段を採用することができる。また、重み付け加算ではなく、重み付け乗算などを採用としてもよい。

「標準偏差」
上述の例では、評価値として「標準偏差」を用いた。しかし、ばらつきが最小となる評価値として用いるのであれば、分散を用いても同一の結果が得られるため、分散を用いてもよい。

「分野」
本実施形態によるＣＴ画像作成は、医療分野だけでなく、材料分析など各種の分野においても使用できる。

「重み付け演算の対象」
上述の実施形態では、２つの撮影ＣＴ画像を対象として、重み付け演算を行った。しかし、３以上の撮影ＣＴ画像を対象として、重みの組み合わせを変更して最適化演算を行い、最適な組み合わせを決定して、単色ＣＴ画像を得てもよい。また、３以上の撮影ＣＴ画像について、複数の２つの組み合わせについて、それぞれ最適化演算を行い、得られた最適な単色ＣＴ画像の中で一番標準偏差が小さいものを選択してもよい。

１０対象物、１２線源、１４電圧発生部、１６移動制御部、１８検出器、２０信号処理部、２２画像処理部、２４記憶部、３０制御部、３２入力部、３４出力部。

Claims

複数のエネルギーレベルの放射線を照射して得た複数のＣＴ画像について演算処理して単色ＣＴ画像を得る単色ＣＴ画像作成方法であって、
第１エネルギーレベルの放射線を照射して第１エネルギーレベルＣＴ画像を得、
第２エネルギーレベルの放射線を照射して第２エネルギーレベルＣＴ画像を得、
得られた第１および第２エネルギーレベルＣＴ画像について複数種類の重み付け演算し、それぞれの重み付け演算結果における複数の単色ＣＴ画像を得、
複数の単色ＣＴ画像のそれぞれについて、高吸収体領域の周辺領域を周方向において所定の面積ごとに複数の注目領域に分割し、各注目領域の画像データの平均値を用い周辺領域についての標準偏差を算出し、標準偏差が小さい単色ＣＴ画像を選択する、
単色ＣＴ画像作成方法。
請求項１に記載の単色ＣＴ画像作成方法であって、
第１および第２エネルギーレベルＣＴ画像について複数種類の重み付け演算は、重み付け加算である、
単色ＣＴ画像作成方法。
請求項１または２に記載の単色ＣＴ画像作成方法であって、
高吸収体の周辺領域は、放射線を照射して得た複数のＣＴ画像において、高吸収体アーチファクトの影響が出る領域に基づいて決定する、
単色ＣＴ画像作成方法。
請求項１〜３のいずれか１つに記載の単色ＣＴ画像作成方法であって、
高吸収体の周辺領域は、放射線を照射して得た複数のＣＴ画像において、高吸収体アーチファクトが存在すると推定される領域の周辺における近接領域を除外した領域である、
単色ＣＴ画像作成方法。
複数のエネルギーレベルの放射線を照射して得た複数のＣＴ画像について演算処理して、単色ＣＴ画像を得る単色ＣＴ画像作成装置であって、
第１エネルギーレベルの放射線を照射して得た第１エネルギーレベルＣＴ画像と、第２エネルギーレベルの放射線を照射して得た第２エネルギーレベルＣＴ画像を取得し、
取得した第１および第２エネルギーレベルＣＴ画像について複数種類の重み付け演算し、それぞれの重み付け演算結果における複数の単色ＣＴ画像を算出し、
複数の単色ＣＴ画像のそれぞれについて、高吸収体の周辺領域を周方向において所定の面積ごとに複数の注目領域に分割し、各注目領域の画像データの平均値を用い周辺領域についての標準偏差を算出し、標準偏差が小さい単色ＣＴ画像を選択する、
単色ＣＴ画像作成装置。
複数のエネルギーレベルの放射線を照射して得た複数のＣＴ画像について演算処理して、単色ＣＴ画像を得る単色ＣＴ画像作成プログラムであって、
単色作成プログラムは、コンピュータに、
第１エネルギーレベルの放射線を照射して得た第１エネルギーレベルＣＴ画像と、第２エネルギーレベルの放射線を照射して得た第２エネルギーレベルＣＴ画像を取得させ、
取得した第１および第２エネルギーレベルＣＴ画像について複数種類の重み付け演算し、それぞれの重み付け演算結果における複数の単色ＣＴ画像を算出させ、
複数の単色ＣＴ画像のそれぞれについて、高吸収体の周辺領域を周方向において所定の面積ごとに複数の注目領域に分割し、各注目領域の画像データの平均値を用い周辺領域についての標準偏差を算出し、標準偏差が小さい単色ＣＴ画像を選択させる、
単色ＣＴ画像作成プログラム。