JP7247431B2 - 減弱係数マップ作成装置、減弱係数マップ作成方法、およびプログラム - Google Patents
減弱係数マップ作成装置、減弱係数マップ作成方法、およびプログラム Download PDFInfo
- Publication number
- JP7247431B2 JP7247431B2 JP2018092351A JP2018092351A JP7247431B2 JP 7247431 B2 JP7247431 B2 JP 7247431B2 JP 2018092351 A JP2018092351 A JP 2018092351A JP 2018092351 A JP2018092351 A JP 2018092351A JP 7247431 B2 JP7247431 B2 JP 7247431B2
- Authority
- JP
- Japan
- Prior art keywords
- region
- image
- bone
- coefficient map
- attenuation coefficient
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims description 13
- 210000000988 bone and bone Anatomy 0.000 claims description 68
- 238000002603 single-photon emission computed tomography Methods 0.000 claims description 64
- 238000012937 correction Methods 0.000 claims description 38
- 210000004072 lung Anatomy 0.000 claims description 13
- 210000004872 soft tissue Anatomy 0.000 claims description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 10
- 210000004185 liver Anatomy 0.000 claims description 5
- 238000012545 processing Methods 0.000 description 41
- 238000009206 nuclear medicine Methods 0.000 description 30
- 238000013500 data storage Methods 0.000 description 10
- 238000003384 imaging method Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 230000005855 radiation Effects 0.000 description 7
- 238000006243 chemical reaction Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 230000002238 attenuated effect Effects 0.000 description 4
- 238000010521 absorption reaction Methods 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000000877 morphologic effect Effects 0.000 description 3
- 238000011158 quantitative evaluation Methods 0.000 description 3
- 229940121896 radiopharmaceutical Drugs 0.000 description 3
- 239000012217 radiopharmaceutical Substances 0.000 description 3
- 230000002799 radiopharmaceutical effect Effects 0.000 description 3
- 208000001132 Osteoporosis Diseases 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 229940079593 drug Drugs 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 201000010235 heart cancer Diseases 0.000 description 1
- 208000019622 heart disease Diseases 0.000 description 1
- 208000024348 heart neoplasm Diseases 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 230000015654 memory Effects 0.000 description 1
- 230000002503 metabolic effect Effects 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 230000006335 response to radiation Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Images
Landscapes
- Nuclear Medicine (AREA)
Description
本発明は、陽電子放射断層撮像（以下、「ＰＥＴ」と称す）画像及び単一光子放射断層撮像（以下、「ＳＰＥＣＴ」と称す）画像の減弱補正に用いる減弱係数マップを作成する技術に関し、特に骨ＳＰＥＣＴの減弱補正に用いられる減弱係数マップを作成する技術に関する。
ＰＥＴ画像およびＳＰＥＣＴ画像に代表される核医学画像は、心臓疾患や癌をはじめとする種々の疾患の診断に有効である。これらの画像は、特定の放射性同位元素でラベルされた薬剤（以下、「放射性医薬品」と称す）を投与し、該薬剤より直接的または間接的に放出されたγ線を専用のカメラによって検出し、再構築することによって得られる。核医学画像は、疾患に対する特異度や感度が高いという優れた性質を有しているばかりでなく、病変部の機能に関する情報を得ることができるという他の診断画像にはない特徴を有している。
核医学画像は被験者に投与された放射性医薬品から直接的または間接的に放出されたγ線を検出することにより得られる。γ線は生体内において構成物質と吸収や散乱といった相互作用を起こし、減弱を受けるので、検出されたγ線のカウントは体表からの深さに応じて減弱され、過小評価されたものとなる。従って、検出されたγ線カウントを用いて核医学画像を正確に再構築するためには、生体におけるγ線の減弱を考慮した減弱補正を行う必要がある。
減弱補正とは、撮像された核医学画像空間のある断面内において、任意の線源（例えば、Tc-99m線源）から放出されたγ線が断面内で放射状に空間内を進行し、最外側に位置する検出器に到達するまでに生じる光子減弱を計算することにより、その減弱分を検出データに上乗せする方法である。
このような減弱補正を行う手法としては、外部線源を用いた減弱補正と、Ｘ線ＣＴを用いた減弱補正法がある（非特許文献１）。外部線源を用いた減弱補正法では、減弱体がある場合のトランスミッションデータと減弱体がない場合のブランクデータとを用いて、減弱係数のマップを求める。Ｘ線ＣＴを用いた減弱補正法では、被験者のＸ線ＣＴ画像を撮像し、そのＣＴ値より線減弱係数を見積もって減弱係数マップを求める。
氏家盛通、宇野和也、廣野圭司、森田立美「ＰＥＴ検査における吸収補正法（外部線源とＸ線ＣＴ）の比較～吸収補正法と被ばくに関する臨床現場での対応～」診療放射線技術科学科論集第６号（２００８）
近年、ＳＰＥＣＴ／ＣＴ装置の開発により骨検査領域での定量評価が注目されている。この定量評価には、減弱補正の実施が必須である。全身検査を行う骨ＳＰＥＣＴに対応するＣＴデータを撮像するにはＳＰＥＣＴ／ＣＴ装置が必要である。しかし、ＳＰＥＣＴ／ＣＴ装置を保有していない施設も少なからず存在し、施設によっては、骨ＳＰＥＣＴデータの減弱補正を行い、骨検査領域での定量評価を行うことが困難な場合があった。
従来、ＳＰＥＣＴ／ＣＴ装置を有しない施設では、Ｃｈａｎｇ法を用いた減弱補正が行われていたが、骨ＳＰＥＣＴの減弱補正の方法は体幹部の構造が複雑であるため従来からあるＣｈａｎｇ法の使用が困難であった。そこで、本発明は、減弱係数マップの新しい作成技術を提供することを目的とする。
本発明の減弱係数マップ作成装置は、被験者の核医学画像データを入力する入力部と、前記核医学画像の中で減弱の程度が均一になる領域に関心領域を設定する関心領域設定部と、前記関心領域の平均カウント値の定数倍を閾値として、被験者の核医学画像を、平均カウント値が前記閾値以上の第１の領域、前記閾値未満の第２の領域に分類する領域判定部と、テンプレートのＣＴ画像において、前記第１の領域の画素値を骨領域に対応するＣＴ値に変換し、前記第２の領域の画素値を骨以外の領域に対応するＣＴ値に変換して、仮想のＣＴ画像を生成し、前記仮想のＣＴ画像を用いて前記核医学画像の減弱補正を行う減弱係数マップを作成する減弱係数マップ作成部とを備える。
本発明の減弱係数マップ作成装置は、被験者の核医学画像データを入力する入力部と、前記核医学画像の中で減弱の程度が均一になる領域に関心領域を設定する関心領域設定部と、前記関心領域の平均カウント値の定数倍を閾値として、被験者の核医学画像を、平均カウント値が前記閾値以上の第１の領域、前記閾値未満の第２の領域に分類する領域判定部と、前記第１の領域の画素値を骨領域に対応するＣＴ値に変換し、前記第２の領域の画素値を骨以外の領域に対応するＣＴ値に変換すると共に、ＣＴ画像を構成するＣＴ情報を付加して、仮想のＣＴ画像を生成し、前記仮想のＣＴ画像を用いて前記核医学画像の減弱補正を行う減弱係数マップを作成する減弱係数マップ作成部とを備える。
本発明によれば、減弱の程度が均一になる領域に設定した関心領域の平均カウント値の定数倍を閾値として、被験者の核医学画像の領域を分類することにより、適切な減弱係数マップを作成できる。また、本発明は、分類した領域に基づいてＣＴ値を設定した仮想のＣＴ画像を生成し、仮想のＣＴ画像から減弱係数マップを生成するので、被験者のＸ線ＣＴ画像から減弱係数マップを作成する従来法（例えば、非特許文献１）を利用することができる。なお、テンプレートのＣＴ画像としては、被験者のＣＴ画像を用いてもよいし、被験者以外の平均的なＣＴ画像を用いてもよい。
ここで、前記関心領域設定部は、被験者の核医学画像の中で減弱の程度が均一になる領域に関心領域を設定できればよい。減弱の程度が均一になる領域とは、放射性医薬品が均一に分布し、そこから発せられる放射線が同程度に減弱される領域である。具体例として、関心領域設定部は、肋骨の領域を避けて肺に関心領域を設定してもよく、又は肝臓に設定してもよい。
骨や軟部組織が写っていない肺の領域は、放射線の減弱に関しては均一で、かつ、空気とほぼ同じとみなすことができるので、肺に関心領域を設定することにより、実質的に減弱のないカウント値を基準とすることができる。また、肝臓の領域は均一に放射能が分布するため減弱の程度が均一であるので、関心領域を設定することで適切な閾値を得られる。なお、肺や肝臓以外であっても、前記核医学画像の中に減弱の程度が均一と評価される領域がある場合には、その領域内に関心領域を設定し、閾値の基準となる平均カウント値を求めてもよい。
本発明の減弱係数マップ作成装置において、前記領域判定部は、前記第２の領域のうち、カウント値が前記平均カウント値以上の領域を第３の領域、前記平均カウント値未満の領域を第４の領域に分類し、前記減弱係数マップ作成部は、前記第３の領域の画素値を軟部組織又は水領域に対応するＣＴ値に変換し、前記第４の領域の画素値を空気領域に対応するＣＴ値に変換して、仮想のＣＴ画像を生成し、前記核医学画像の減弱補正を行う減弱係数マップを作成する減弱係数マップ作成部とを備えてもよい。この構成により、骨以外の領域を軟部組織又は水領域と、空気領域とに分けて、各領域に対応するＣＴ値に変換して、仮想のＣＴ画像を生成し、きめ細かく減弱係数を設定した減弱係数マップを作成することができる。
本発明の減弱係数マップ作成装置において、前記核医学画像は、骨ＳＰＥＣＴ画像データであってもよい。本発明は、体幹部の構造が複雑である骨ＳＰＥＣＴ画像データに好適に適用することができる。
本発明は、被験者の核医学画像に基づいて、適切に減弱係数マップを作成することができる。
以下、本発明の実施の形態の減弱係数マップ作成装置について図面を参照して説明する。以下に説明する実施の形態では、核医学画像として骨ＳＰＥＣＴ画像を用いる例を挙げるが、本発明は、骨以外のＳＰＥＣＴ画像やＰＥＴ画像等の核医学画像についても適用することができる。
（第１の実施の形態）
図１は、第１の実施の形態の減弱係数マップ作成装置を備えた画像処理装置１０の構成を示す図である。
画像処理装置１０は、核医学画像撮像装置３０にて撮像された被験者の骨ＳＰＥＣＴ画像データの減弱補正を行う機能を有している。画像処理装置１０は、核医学画像撮像装置３０にて撮像された骨ＳＰＥＣＴ画像データを入力する画像入力部１１と、骨ＳＰＥＣＴ画像データの減弱補正を行う演算処理部１２と、減弱補正された画像を出力する出力部１３とを有している。画像入力部１１は、被験者のＣＴ画像を入力する機能も有している。なお、ここで入力される被験者のＣＴ画像は、ＳＰＥＣＴ／ＣＴ装置で撮像されたＣＴ画像のようにＳＰＥＣＴ画像との位置合わせがなされたものでなくてもよい。入力されたＣＴ画像は、ＣＴ形態画像データ記憶部１６に記憶される。
図１は、第１の実施の形態の減弱係数マップ作成装置を備えた画像処理装置１０の構成を示す図である。
画像処理装置１０は、核医学画像撮像装置３０にて撮像された被験者の骨ＳＰＥＣＴ画像データの減弱補正を行う機能を有している。画像処理装置１０は、核医学画像撮像装置３０にて撮像された骨ＳＰＥＣＴ画像データを入力する画像入力部１１と、骨ＳＰＥＣＴ画像データの減弱補正を行う演算処理部１２と、減弱補正された画像を出力する出力部１３とを有している。画像入力部１１は、被験者のＣＴ画像を入力する機能も有している。なお、ここで入力される被験者のＣＴ画像は、ＳＰＥＣＴ／ＣＴ装置で撮像されたＣＴ画像のようにＳＰＥＣＴ画像との位置合わせがなされたものでなくてもよい。入力されたＣＴ画像は、ＣＴ形態画像データ記憶部１６に記憶される。
演算処理部１２は、減弱補正を行うために用いる減弱係数マップを作成する減弱係数マップ作成装置に相当する処理部２０と、作成された減弱係数マップを用いて骨ＳＰＥＣＴ画像の減弱補正を行う減弱補正部２８と、減弱補正された画像を再構成する画像再構成部２９とを有している。
減弱係数マップ作成装置に相当する処理部２０について説明する。関心領域設定部２１は、被験者の骨ＳＰＥＣＴ画像に対して関心領域（以下、「ＶＯＩ」とも称す）を設定する機能を有する。関心領域設定部２１は、被験者の骨ＳＰＥＣＴデータ画像の体軸断面（axial）像、冠状断面（coronal）像を参照して、被験者の右肺上葉の肋骨のない領域内にＶＯＩを設定する。ＶＯＩは、一例として、直径６０ｍｍ程度の球形の形状を有する。設定するＶＯＩの位置や形状のデータは、ＶＯＩデータ記憶部１４に記憶されており、関心領域設定部２１は、ＶＯＩデータ記憶部１４からＶＯＩの形状のデータを読み出して用いる。
ＶＯＩ平均値算出部２２は、設定したＶＯＩの平均カウント値を算出する機能を有する。ここで、平均カウント値を「Ｌ」とする。数値比較部２３は、被験者の骨ＳＰＥＣＴ画像データの各ピクセルのカウント値と閾値とを比較する。第１の閾値として平均カウント値Ｌの５倍である５Ｌを用い、第２の閾値として平均カウント値Ｌを用いる。第１の閾値、第２の閾値としてどのような値を用いるかは、領域区分閾値データ記憶部１５に記憶されている。数値比較部２３は、領域区分閾値データ記憶部１５から閾値のデータを取得して用いる。ここでは、第１の閾値として、平均カウント値の５倍を用いているが、５倍というのは一例であり、例えば、平均カウント値Ｌの他の定数倍の数値（４倍とか６倍等）を用いてもよい。
領域判定部２４は、数値比較部２３による比較結果に基づいて、平均カウント値が第１の閾値以上の領域をＡ領域、平均カウント値が第２の閾値以上、第１の閾値未満の領域をＢ領域、平均カウント値が第２の閾値未満の領域をＣ領域と判定する。
数値変換部２５は、ＣＴ形態画像データ記憶部１６に記憶されたＣＴ画像を読み出し、読み出したＣＴ画像のＣＴ値を領域の種類（Ａ領域、Ｂ領域、Ｃ領域の種類）に応じて変換する。つまり、読み出したＣＴ画像をテンプレートとし、領域の種類に応じたＣＴ値に置き換えていく。Ａ領域と判定された領域はピクセル値を３００または４００に変換し、Ｂ領域と判定された領域はピクセル値を５０とし、Ｃ領域と判定された領域は、ピクセル値を－１０００とする。なお、Ａ領域においてピクセル値を３００と４００のいずれにするかは、被験者の年齢や症例によって判定する。被験者が骨粗鬆症である場合、または５０歳以上の場合にはピクセル値を３００とし、それ以外では４００とする。これにより、被験者の仮想ＣＴ画像が生成される。なお、ここで変換に用いた数値は一例であり、他の値に変換してもよい。
減弱係数マップ作成部２６は、数値変換部２５で生成された被験者の仮想ＣＴ画像を用いて、減弱係数マップを作成する。減弱係数マップ作成部２６による減弱係数マップの生成には、従来の方法を適用することができる。
図２は、第１の実施の形態の画像処理装置１０の動作を示す図である。まず、核医学画像撮像装置３０にて被験者の骨ＳＰＥＣＴ画像を撮像し（Ｓ１０）、続いて被験者のＣＴ画像データを読み込んで入力する（Ｓ１２）。画像処理装置１０は、核医学画像撮像装置３０にて撮像された骨ＳＰＥＣＴ画像データを読み込んで入力する（Ｓ１４）。
画像処理装置１０は、関心領域設定部２１にて骨ＳＰＥＣＴ画像データの右肺上葉に球形のＶＯＩを設定し（Ｓ１６）、ＶＯＩの平均カウント値Ｌを算出する（Ｓ１８）。続いて、画像処理装置１０は、平均カウント値Ｌから求めた第１の閾値（５Ｌ）、第２の閾値（Ｌ）を用いて、骨ＳＰＥＣＴ画像データのＡ領域、Ｂ領域、Ｃ領域に分類する（Ｓ２０）。次に、画像処理装置１０は、Ａ領域は骨、Ｂ領域は軟部組織又は水、Ｃ領域は空気として、ＣＴ画像データ記憶部１６から読み出したＣＴ画像のデータを、領域の種類に応じたＣＴ値に変換し、仮想のＣＴ画像を生成する（Ｓ２２）。画像処理装置１０は、仮想のＣＴ画像に基づいて減弱係数マップを作成する（Ｓ２４）。
次に、画像処理装置１０は、作成された減弱係数マップを用いて、骨ＳＰＥＣＴ画像データを減弱補正し、画像再構成を行う（Ｓ２６）。そして、画像処理装置１０は、再構成された画像を出力する（Ｓ２８）。
以上、本実施の形態の減弱係数マップ作成装置を備えた画像処理装置１０の構成について説明したが、上記した画像処理装置１０のハードウェアの例は、ＣＰＵ、ＲＡＭ、ＲＯＭ、ハードディスク、ディスプレイ、キーボード、マウス、通信インターフェース等を備えたコンピュータである。上記した各機能を実現するモジュールを有するプログラムをＲＡＭまたはＲＯＭに格納しておき、ＣＰＵによって当該プログラムを実行することによって、上記した画像処理装置１０が実現される。このようなプログラムも本発明の範囲に含まれる。
第１の実施の形態の減弱係数マップ作成装置を備えた画像処理装置１０は、被験者のＸ線ＣＴ画像データをテンプレートとし、被験者の骨ＳＰＥＣＴ画像データに基づいて仮想のＣＴ画像を生成し、仮想のＣＴ画像データに基づいて減弱係数マップを作成することができる。本実施の形態では、被験者のＣＴ画像をテンプレートとして用いて仮想のＣＴ画像を生成する例を挙げたが、テンプレートとして用いるＣＴ画像は被験者のＣＴ画像でなくてもよく、第三者のＣＴ画像をテンプレートとして用いてもよい。
第１の実施の形態では、被験者の骨ＳＰＥＣＴ画像を、骨、軟部組織又は水、空気のそれぞれに対応すると考えられる領域に分類することによって、体幹部の構造が複雑である骨ＳＰＥＣＴ画像の減弱係数マップを適切に作成することができる。
また、本実施の形態で用いるＣＴ画像は、ＳＰＥＣＴ／ＣＴ画像で撮像される画像のように、ＳＰＥＣＴ画像との位置合わせがなされたものでなくてもよいので、ＳＰＥＣＴ／ＣＴ撮像装置を保有していない施設でも本実施の形態を用いることができる。
（第２の実施の形態）
図３は、第２の実施の形態の減弱係数マップ作成装置を備えた画像処理装置１０の構成を示す図である。第２の実施の形態の画像処理装置１０は、第１の実施の形態の画像処理装置１０と同様に、被験者の骨ＳＰＥＣＴ画像データをＡ領域、Ｂ領域、Ｃ領域に分類する。第２の実施の形態の画像処理装置１０では、仮想のＣＴ画像から減弱係数マップを作成するために、ＣＴ画像のテンプレートを用いず、ＣＴ画像のＣＴ情報を用いる。以下、第１の実施の形態と相違する点を中心に説明する。
図３は、第２の実施の形態の減弱係数マップ作成装置を備えた画像処理装置１０の構成を示す図である。第２の実施の形態の画像処理装置１０は、第１の実施の形態の画像処理装置１０と同様に、被験者の骨ＳＰＥＣＴ画像データをＡ領域、Ｂ領域、Ｃ領域に分類する。第２の実施の形態の画像処理装置１０では、仮想のＣＴ画像から減弱係数マップを作成するために、ＣＴ画像のテンプレートを用いず、ＣＴ画像のＣＴ情報を用いる。以下、第１の実施の形態と相違する点を中心に説明する。
第２の実施の形態の減弱係数マップ作成装置は、ＣＴ画像を構成するＣＴ情報を記憶したＣＴ情報記憶部１７を有している。ＣＴ情報は、ＣＴ画像を構成する情報、例えば撮像条件、再構成条件等である。ＣＴ画像は、被験者の形態を表すＣＴ値とＣＴ情報を含んで構成される。なお、本実施の形態では、ＣＴ情報記憶部１７に記憶されたＣＴ情報は、ＤＩＣＯＭ等のコンピュータで利用可能な形式にて、ＣＤ、ＤＶＤ、ハードディスク、あるいは半導体メモリ等のコンピュータで読取り可能な媒体に保存された状態で提供されるデータを用いてもよい。
数値変換部２５は、前記核医学画像のピクセル値を領域の種類（Ａ領域、Ｂ領域、Ｃ領域の種類）に応じたＣＴ値に変換する。Ａ領域と判定された領域はピクセル値を３００または４００に変換し、Ｂ領域と判定された領域はピクセル値を５０とし、Ｃ領域と判定された領域は、ピクセル値を－１０００とする。なお、Ａ領域においてピクセル値を３００と４００のいずれにするかは、被験者の年齢や症例によって判定する。被験者が骨粗鬆症である場合、または５０歳以上の場合にはピクセル値を３００とし、それ以外では４００とする。これにより、被験者の仮想のＣＴ画像が生成される。なお、ここで変換に用いた数値は一例であり、他の値に変換してもよい。変換されたＣＴ値と、ＣＴ情報記憶部１７から読み出したＣＴ情報とに基づいて仮想のＣＴ画像を生成する。
減弱係数マップ作成部２６は、仮想のＣＴ画像に基づいて減弱係数マップを作成する。減弱係数マップ出力部２７は、作成した減弱係数マップを減弱補正部２８に対して出力する。これにより、減弱補正部２８は、減弱係数マップを用いて減弱補正を行うことができる。
図４は、第２の実施の形態の画像処理装置１０の動作を示す図である。まず、核医学画像撮像装置３０にて被験者の骨ＳＰＥＣＴ画像を撮像する（Ｓ１０）。画像処理装置１０は、核医学画像撮像装置３０にて撮像された骨ＳＰＥＣＴ画像データを読み込んで入力する（Ｓ１４）。
画像処理装置１０は、ＶＯＩ設定部２１にて骨ＳＰＥＣＴ画像データの右肺上葉に球形のＶＯＩを設定し（Ｓ１６）、ＶＯＩの平均カウント値Ｌを算出する（Ｓ１８）。続いて、画像処理装置１０は、平均カウント値Ｌから求めた第１の閾値（５Ｌ）、第２の閾値（Ｌ）を用いて、骨ＳＰＥＣＴ画像データのＡ領域、Ｂ領域、Ｃ領域に分類する（Ｓ２０）。画像処理装置１０は、Ａ領域は骨、Ｂ領域は軟部組織又は水、Ｃ領域は空気として、各領域に対応するＣＴ値に変換し、変換によって得られたＣＴ値とＣＴ情報記憶部１７から読み出したＣＴ情報とにより仮想のＣＴ画像を生成する。画像処理装置１０は、仮想のＣＴ画像に基づいて、各領域に対応する減弱係数を設定して減弱係数マップを作成する（Ｓ２４）。
次に、画像処理装置１０は、作成された減弱係数マップを用いて、骨ＳＰＥＣＴ画像データを減弱補正し、画像再構成を行う（Ｓ２６）。そして、画像処理装置１０は、再構成された画像を出力する（Ｓ２８）。
以上、第２の実施の形態の画像処理装置１０の構成および動作について説明した。第２の実施の形態の画像処理装置１０は、第１の実施の形態の画像処理装置１０と同様に、被験者の骨ＳＰＥＣＴ画像に基づいて減弱係数マップを作成することができる。また、第２の実施の形態では、ＣＴ情報と被験者の骨ＳＰＥＣＴ画像データに基づいて、仮想のＣＴ画像を生成することにより、仮想のＣＴ画像に対して、Ｘ線ＣＴ画像から減弱係数マップを作成する従来の手法を用いることができる。また、本実施の形態の画像処理装置１０は、被験者のＣＴ画像なしに減弱係数マップを生成できる。
ＳＰＥＣＴ画像のカウント値は骨組織の代謝活性だけではなく、組織血流や血中クリアランスなどに依存するため、肺の平均カウント値を基準として変換された仮想ＣＴ画像データにおける骨・軟部組織・空気の範囲は、真のＣＴデータと比較すると多少異なる。しかし、仮想ＣＴデータを用いた減弱補正によるＳＰＥＣＴ画像のカウント値と、真のＣＴデータを用いた減弱補正によるＳＰＥＣＴ画像のカウント値の誤差は少ないことを臨床データによって確認した。したがって、本実施の形態のように仮想のＣＴ画像を用いた方法によって、適切な減弱係数マップを作成することができる。
以上、本発明の減弱係数マップ作成装置の構成について実施の形態を挙げて詳細に説明したが、本発明の減弱係数マップ作成装置は、上記した実施の形態に限定されるものではない。
上記した実施の形態では、骨ＳＰＥＣＴ画像データの領域をＡ～Ｃの３つの領域に分類する例を挙げたが、必ずしも３つの領域に分類しなくてもよく、例えば、第１の閾値（５Ｌ）を境界として２つの領域に分類してもよい。骨であると想定される領域か否かによって領域を分類することにより、減弱係数マップを作成することができる。
現段階では、ＳＰＥＣＴ画像に基づいて軟部組織と水を分離することは困難である。その一方で、軟部組織と水を分離することによる減弱補正の精度向上は僅かであると考えられるので、本実施の形態では分離していない。ただし、将来に、軟部組織または水を分離可能な技術が現れたときには、Ｂ領域を、軟部組織と水の領域に分類してもよい。
上記した実施の形態では、右肺上葉にＶＯＩを設定する例を挙げたが、ＶＯＩの設定領域は必ずしも右肺上葉に限られず、左肺上葉や、右肺中葉等でもよい。ただし、右肺上葉が、肋骨と心臓を避ける意味で好適である。また、被験者の骨ＳＰＥＣＴ画像に肺が含まれていない場合には、減弱の程度が均一であると評価される領域内にＶＯＩを設定してもよい。減弱の程度が均一であると評価される領域として、肺以外に、肝臓にＶＯＩを設定してもよい。
本発明は、核医学画像の減弱補正に用いる減弱係数マップを作成する装置として有用である。
１０ 画像処理装置
１１ 画像入力部
１２ 演算処理部
１３ 出力部
１４ ＶＯＩデータ記憶部
１５ 領域区分閾値データ記憶部
１６ ＣＴ形態画像データ記憶部
１７ ＣＴ情報記憶部
２０ 減弱係数マップ作成装置の処理部
２１ 関心領域設定部
２２ ＶＯＩ平均値算出部
２３ 数値比較部
２４ 領域判定部
２５ 数値変換部
２６ 減弱係数マップ作成部
２７ 減弱係数マップ出力部
２８ 減弱補正部
２９ 画像再構成部
３０ 核医学画像撮像装置
１１ 画像入力部
１２ 演算処理部
１３ 出力部
１４ ＶＯＩデータ記憶部
１５ 領域区分閾値データ記憶部
１６ ＣＴ形態画像データ記憶部
１７ ＣＴ情報記憶部
２０ 減弱係数マップ作成装置の処理部
２１ 関心領域設定部
２２ ＶＯＩ平均値算出部
２３ 数値比較部
２４ 領域判定部
２５ 数値変換部
２６ 減弱係数マップ作成部
２７ 減弱係数マップ出力部
２８ 減弱補正部
２９ 画像再構成部
３０ 核医学画像撮像装置
Claims (6)
- 被験者の骨ＳＰＥＣＴ画像データを入力する入力部と、
前記骨ＳＰＥＣＴ画像の中で減弱の程度が均一であると評価される領域内に関心領域を設定する関心領域設定部と、
前記関心領域の平均カウント値の定数倍を閾値として、被験者の骨ＳＰＥＣＴ画像を、カウント値が前記閾値以上の第１の領域、前記閾値未満の第２の領域に分類する領域判定部と、
前記第１の領域の画素値を骨領域に対応するＣＴ値に変換し、前記第２の領域の画素値を骨以外の領域に対応するＣＴ値に変換すると共に、ＣＴ画像を構成するＣＴ情報を付加して、仮想のＣＴ画像を生成し、前記仮想のＣＴ画像を用いて前記骨ＳＰＥＣＴ画像の減弱補正を行う減弱係数マップを作成する減弱係数マップ作成部と、
を備える減弱係数マップ作成装置。 - 前記関心領域設定部は、肋骨の領域を避けて肺領域に関心領域を設定する請求項１に記載の減弱係数マップ作成装置。
- 前記関心領域設定部は、肝臓領域に関心領域を設定する請求項１に記載の減弱係数マップ作成装置。
- 前記領域判定部は、前記第２の領域のうち、カウント値が前記平均カウント値以上の領域を第３の領域、前記平均カウント値未満の領域を第４の領域に分類し、
前記減弱係数マップ作成部は、前記第３の領域の画素値を軟部組織又は水領域に対応するＣＴ値に変換し、前記第４の領域の画素値を空気領域に対応するＣＴ値に変換して、仮想のＣＴ画像を生成し、前記骨ＳＰＥＣＴ画像の減弱補正を行う減弱係数マップを作成する減弱係数マップ作成部と、
を備える請求項１乃至３のいずれかに記載の減弱係数マップ作成装置。 - 被験者の骨ＳＰＥＣＴ画像から当該骨ＳＰＥＣＴ画像の減弱補正に用いる減弱係数マップを生成する方法であって、
被験者の骨ＳＰＥＣＴ画像データを入力するステップと、
前記骨ＳＰＥＣＴ画像の中で減弱の程度が均一であると評価される領域内に関心領域を設定するステップと、
前記関心領域の平均カウント値の定数倍を閾値として、被験者の骨ＳＰＥＣＴ画像を、カウント値が前記閾値以上の第１の領域、前記閾値未満の第２の領域に分類するステップと、
前記第１の領域の画素値を骨領域に対応するＣＴ値に変換し、前記第２の領域の画素値を骨以外の領域に対応するＣＴ値に変換すると共に、ＣＴ画像を構成するＣＴ情報を付加して、仮想のＣＴ画像を生成するステップと、
前記仮想のＣＴ画像を用いて前記骨ＳＰＥＣＴ画像の減弱補正を行う減弱係数マップを作成するステップと、
を備える減弱係数マップ作成方法。 - 被験者の骨ＳＰＥＣＴ画像から当該骨ＳＰＥＣＴ画像の減弱補正に用いる減弱係数マップを生成するためのプログラムであって、コンピュータに、
被験者の骨ＳＰＥＣＴ画像データを入力するステップと、
前記骨ＳＰＥＣＴ画像の中で減弱の程度が均一であると評価される領域内に関心領域を設定するステップと、
前記関心領域の平均カウント値の定数倍を閾値として、被験者の骨ＳＰＥＣＴ画像を、カウント値が前記閾値以上の第１の領域、前記閾値未満の第２の領域に分類するステップと、
前記第１の領域の画素値を骨領域に対応するＣＴ値に変換し、前記第２の領域の画素値を骨以外の領域に対応するＣＴ値に変換すると共に、ＣＴ画像を構成するＣＴ情報を付加して、仮想のＣＴ画像を生成するステップと、
前記仮想のＣＴ画像を用いて前記骨ＳＰＥＣＴ画像の減弱補正を行う減弱係数マップを作成するステップと、
を実行させるプログラム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2018092351A JP7247431B2 (ja) | 2018-05-11 | 2018-05-11 | 減弱係数マップ作成装置、減弱係数マップ作成方法、およびプログラム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2018092351A JP7247431B2 (ja) | 2018-05-11 | 2018-05-11 | 減弱係数マップ作成装置、減弱係数マップ作成方法、およびプログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2019197026A JP2019197026A (ja) | 2019-11-14 |
| JP7247431B2 true JP7247431B2 (ja) | 2023-03-29 |
Family
ID=68538359
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2018092351A Active JP7247431B2 (ja) | 2018-05-11 | 2018-05-11 | 減弱係数マップ作成装置、減弱係数マップ作成方法、およびプログラム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP7247431B2 (ja) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007014759A (ja) | 2005-06-17 | 2007-01-25 | Cti Molecular Imaging Inc | Ｐｅｔ／ｃｔ撮像における画像ベースのアーティファクト低減方法 |
| JP2009047602A (ja) | 2007-08-21 | 2009-03-05 | Toshiba Corp | 陽電子放出コンピュータ断層撮影装置、減弱マップ作成装置および減弱マップ作成プログラム |
| JP2009525780A (ja) | 2006-02-03 | 2009-07-16 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Ｃｔベースの減衰マップを作成するときの異質対象物の明示 |
| CN107115119A (zh) | 2017-04-25 | 2017-09-01 | 上海联影医疗科技有限公司 | Pet图像衰减系数的获取方法、衰减校正的方法及系统 |
-
2018
- 2018-05-11 JP JP2018092351A patent/JP7247431B2/ja active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007014759A (ja) | 2005-06-17 | 2007-01-25 | Cti Molecular Imaging Inc | Ｐｅｔ／ｃｔ撮像における画像ベースのアーティファクト低減方法 |
| JP2009525780A (ja) | 2006-02-03 | 2009-07-16 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Ｃｔベースの減衰マップを作成するときの異質対象物の明示 |
| JP2009047602A (ja) | 2007-08-21 | 2009-03-05 | Toshiba Corp | 陽電子放出コンピュータ断層撮影装置、減弱マップ作成装置および減弱マップ作成プログラム |
| CN107115119A (zh) | 2017-04-25 | 2017-09-01 | 上海联影医疗科技有限公司 | Pet图像衰减系数的获取方法、衰减校正的方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| Y. Takahashi et al.，Segmented attenuation correction for myocardial SPECT，Annals of Nuclear Medicine，2004年04月，Vol. 18, No. 2，p.137-143 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2019197026A (ja) | 2019-11-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Shi et al. | Deep learning-based attenuation map generation for myocardial perfusion SPECT | |
| CN109741284B (zh) | 用于校正pet成像中由呼吸运动引起的失配的系统和方法 | |
| Ashrafinia | Quantitative nuclear medicine imaging using advanced image reconstruction and radiomics | |
| US10395353B2 (en) | Model-based scatter in multi-modality multi-energy SPECT reconstruction | |
| US10628972B2 (en) | Diagnostic imaging method and apparatus, and recording medium thereof | |
| US9619905B2 (en) | Apparatus and method for generation of attenuation map | |
| US20220207791A1 (en) | Method and system for generating attenuation map from spect emission data | |
| JP2013200306A (ja) | エミッション・データに基づいた核医学イメージングにおける減衰補償のためのシステム及び方法 | |
| Zaidi et al. | Strategies for attenuation compensation in neurological PET studies | |
| US11574154B2 (en) | Patient-adaptive nuclear imaging | |
| Clark et al. | Photon-counting cine-cardiac CT in the mouse | |
| Grosser et al. | Optimization of SPECT-CT hybrid imaging using iterative image reconstruction for low-dose CT: a phantom study | |
| KR102283934B1 (ko) | Pet 영상에 기반한 가상 ct 영상 및 감쇠보정 pet 영상 생성 방법 및 시스템 | |
| US20190114772A1 (en) | Image processing apparatus and medical image taking apparatus | |
| CN107886554B (zh) | 流数据的重构 | |
| US7242004B2 (en) | Image correction method, image correction apparatus, and image correction program | |
| JP4997575B2 (ja) | 減弱係数マップ作成装置、減弱係数マップ作成方法およびプログラム | |
| JP7247431B2 (ja) | 減弱係数マップ作成装置、減弱係数マップ作成方法、およびプログラム | |
| JP6143533B2 (ja) | 核医学画像再構成装置、核医学画像再構成方法、及びプログラム | |
| Zaidi et al. | Attenuation correction strategies in emission tomography | |
| CN115774281A (zh) | 用于断层摄影不一致性检测的基于衰减校正的加权 | |
| JP6052425B2 (ja) | 輪郭画像生成装置および核医学診断装置 | |
| CN114730476A (zh) | 有限角度重建的网络确定 | |
| US11911193B2 (en) | Method and apparatus for generating a resultant image dataset of a patient | |
| JP2019124672A (ja) | 画像処理装置、画像処理方法およびプログラム |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20210325 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20220405 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220525 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20220726 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220922 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20221129 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20221223 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7247431 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |