核医学画像の画像処理方法、画像処理装置、及び画像処理プログラム
【課題】心筋核医学画像において、肝臓や胆嚢といった心筋以外の高集積による心筋画像の影響を除去する方法、当該方法をコンピュータに実行させるプログラム、および当該方法を実行するための装置を提供する。
【解決手段】心筋核医学画像における各投影画像上において、心筋全体に関心領域を設定し、当該関心領域内に心筋部以外の高集積部位が存在している投影画像を削除し、残りの投影画像につき心筋領域に設定した関心領域以外の部位をマスク処理して除去するといった処理を行う。このような処理により、従来は除去が困難であった、心筋部位と心筋以外の高集積部位とが接近している投影画像における当該高集積部位の影響を、効果的に取り除く事が可能となった。
【解決手段】心筋核医学画像における各投影画像上において、心筋全体に関心領域を設定し、当該関心領域内に心筋部以外の高集積部位が存在している投影画像を削除し、残りの投影画像につき心筋領域に設定した関心領域以外の部位をマスク処理して除去するといった処理を行う。このような処理により、従来は除去が困難であった、心筋部位と心筋以外の高集積部位とが接近している投影画像における当該高集積部位の影響を、効果的に取り除く事が可能となった。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、核医学画像の画像処理技術に関する。より詳しくは、心筋核医学画像において、アーチファクトを有効に除去し得る、画像処理技術に関する。
【背景技術】
【０００２】
心臓疾患の画像診断の一つとして、ＳＰＥＣＴに代表される、核医学画像診断法が用いられている。核医学画像診断法では、放射性同位元素で標識された薬剤（放射性医薬品）を被験者に投与し、放出されたγ線を検出して画像化処理を行う事により、診断画像を得る事ができる。心臓疾患の画像診断に用いる放射性医薬品としては、塩化タリウム−２０１、ヘキサキス（２−メトキシイソブチルイソニトリル）テクネチウム（９９ｍＴｃ）、テトロホスミンテクネチウム（９９ｍＴｃ）、３−ヨードベンジルグアニジン（１２３Ｉ）、１５−（４−ヨードフェニル）−３（Ｒ，Ｓ）−メチルデカン酸（１２３Ｉ）等が開発され、臨床において用いられている。
【０００３】
これらの放射性医薬品のうち、放射性同位元素として９９ｍＴｃを用いた製剤（以下、９９ｍＴｃ製剤という）は、キット化できるために汎用性に優れ、また、ＳＰＥＣＴ装置に適したエネルギーを有しているために良好な画像を得る事ができるといった利点を有している。しかし、９９ｍＴｃ製剤は、一般に分子サイズが大きく、肝胆道系に高い集積を有する。このような心筋以外における高集積は、再構成画像上でストリークアーチファクトの原因となり、診断の妨げになるといった問題がある（非特許文献１、非特許文献２、非特許文献３、非特許文献４）。また、放射性ヨード製剤においても、肝臓に高集積し、心筋ＳＰＥＣＴにおいて下壁部が欠損するといったことが報告されている（非特許文献５）。
【０００４】
心筋以外への高集積に起因する再構成画像上でのアーチファクトの発生を防ぐ方法として、マスク処理によって高集積部位を除去する方法がこれまでに検討されている（非特許文献１、非特許文献２、非特許文献３）。これまでの報告によれば、再構成処理前の投影データにマスク処理を施す事により、当該高集積による心筋画像への影響を減少させ得ることが開示されている（非特許文献１、非特許文献２、非特許文献３）。しかし、心筋以外の高集積部位が心筋に近接していた場合には、画像データ上で当該高集積部位を除去する事は難しく、心筋画像への当該高集積部位の影響を完全に取り除く事は困難である（非特許文献２、非特許文献３）。
【先行技術文献】
【０００５】
【非特許文献１】 久野晴丘 他、「９９ｍＴｃ−ｔｅｔｒｏｆｏｓｍｉｎによる胆嚢の高集積が心筋ＳＰＥＣＴ画像に及ぼす影響と改善法についての検討」、日本放射線技術学会誌、２０００年８月、５６、８、ｐ．１０４４−１０５１
【非特許文献２】 今井嘉門 他、「ＭＩＢＩの心筋シンチのアーティファクトの原因に関して」、核医学、１９９５、３２、３、ｐ．３０７−３１０
【非特許文献３】 高木昭浩 他、「９９ｍＴｃ心筋血流製剤を用いたＳＰＥＣＴ撮像における肝の高集積が心筋に及ぼす影響の軽減；マスク処理法の有用性と問題点」、核医学、１９９９、３６、５、ｐ．４５９−４６５
【非特許文献４】 大西英雄 他、「９９ｍＴｃ標識心筋血流製剤を用いた心筋ＳＰＥＣＴ画像への肝臓からの影響」、核医学、１９９８、３５、６、ｐ．３７５−３８３
【非特許文献５】 小林秀樹 他、「肝高集積が原因となって出現する１２３Ｉ−ＭＩＢＧ心筋ＳＰＥＣＴの下壁欠損像の特徴とその対策−ファントムを用いた検討−」、核医学、１９９４、３１、４、ｐ．３５９−３６６
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明は上記事実に鑑みてなされたものであり、心筋以外の高集積の心筋画像への影響を除去する方法、当該方法をコンピュータに実行させるプログラム、および当該方法を実行するための装置を提供する事を目的とした。
【課題を解決するための手段】
【０００７】
発明者は検討の結果、投影データセットにおいて、心筋以外の高集積部位が心筋部位に近接している投影画像を削除しても、残りのデータを用いて作成された再構成画像の画質に本質的な悪影響を与えない事を見出した。この知見に基づき、各投影画像上で心筋全体に関心領域を設定し、当該関心領域内に心筋部以外の高集積部位が存在している投影画像を削除し、残りの投影画像につき前記関心領域以外の部位（すなわち、心筋領域に設定した関心領域以外の部位）をマスク処理等を施して除去するといった処理を行う事により、上記の問題点が解消し、当該高集積部位の影響を効果的に除去しうる事を見出し、本発明を完成させた。
【０００８】
なお、本明細書において、投影データセットとは、一度の撮像処理により得られた各投影角度での投影画像のセットをいう。例えば、患者に対し、２°ステップで３６０°の投影を行った場合、投影角度２°、４°、６°、…、３５８°、３６０°において得られた１８０枚の投影画像のセットが、本明細書における投影データセットとなる。
【０００９】
本発明の一側面に係る画像処理方法は、コンピュータが心筋核医学画像の投影データセットを取得する投影データセット取得ステップと、コンピュータが、前記投影データセットに含まれる各投影画像上で心筋部位を含む領域に関心領域を設定する関心領域設定ステップと、コンピュータが、心筋部以外の高集積部位が前記関心領域と重なっている投影画像を削除する画像削除ステップと、コンピュータが、前記画像削除ステップにて心筋部以外の高集積部位が前記関心領域内と重なっている投影画像を削除した後の投影データセットに含まれる各投影画像につき、前記関心領域設定ステップにて設定した関心領域以外の部分を削除する非関心領域削除ステップと、コンピュータが、非関心領域削除ステップ後の投影データセットを用いて画像再構成を行う画像再構成ステップと、を含む。
【００１０】
本発明の別の一側面に係る画像処理プログラムは、コンピュータに、心筋核医学画像の投影データセットを取得する投影データ取得ステップと、前記投影データセットに含まれる各投影画像上で心筋部位を含む領域に関心領域を設定する関心領域設定ステップと、心筋部以外の高集積部位が前記関心領域と重なっている投影画像を削除する画像削除ステップと、前記画像削除ステップにて心筋部以外の高集積部位が前記関心領域内と重なっている投影画像を削除した後の投影データセットに含まれる各投影画像につき、前記関心領域設定ステップにて設定した関心領域以外の部分を削除する非関心領域削除ステップと、非関心領域削除ステップ後の投影データセットを用いて画像再構成を行う画像再構成ステップと、を実行させる。
【００１１】
本発明のさらに別の一側面に係る画像処理装置は、心筋核医学画像の投影データセットを取り込む投影データセット取得部と、前記投影データセットに含まれる各投影画像上で心筋部位を含む領域に関心領域を設定する関心領域設定部と、前記投影データセットに含まれる各投影画像のうち、心筋部位以外の高集積部位が前記関心領域と重なっている投影画像を削除する画像削除部と、前記関心領域設定部の機能により設定された前記関心領域以外の部分を各投影画像上において削除する非関心領域削除部と、投影データセットを用いて画像再構成を行う画像再構成部と、を備える。
【００１２】
従来技術において指摘されてきたように、投影画像上において心筋部位と心筋以外の高集積部位とが接近している場合には、当該高集積部位を完全に削除する事が難しく、このことが、再構成画像上で当該高集積部位の影響を完全に除去することを困難にしてきた（今井嘉門 他、核医学、１９９５、３２、３、ｐ．３０７−３１０；高木昭浩 他、核医学、１９９９、３６、５、ｐ．４５９−４６５）。本発明に係る画像処理方法および画像処理プログラムでは、心筋部位を包含するよう設定した関心領域と心筋部位以外の高集積とが重なった投影画像を画像ごと削除し、それ以外の投影画像では前記関心領域以外の領域を削除する処理を行っているので、再構成画像において、心筋以外の高集積部位の影響を、効果的に取り除く事ができる。
また同様に、本発明に係る画像処理装置も、心筋部位を包含するよう設定した関心領域と心筋部位以外の高集積とが重なった投影画像を画像ごと削除し、それ以外の投影画像では前記関心領域以外の領域を削除する機能を備えているので、心筋部以外の高集積部位からの影響を完全に取り除いた再構成画像を生成する事が可能となる。
【００１３】
本発明に係る画像処理方法および画像処理プログラムにおいて、心筋部位に設定した関心領域以外の部分を削除する非関心領域削除ステップは、当該関心領域以外の部分をマスク処理によって削除するといった方法や、当該関心領域以外の領域に対応するピクセルにヌルコードやバックグラウンドにおけるカウントの平均値を代入するといった方法にて行う事ができる。
同様に、本発明に係る画像処理装置における非関心領域削除部は、心筋部位に設定した関心領域以外の部分をマスク処理で削除するといった方法や、当該関心領域以外の領域に対応するピクセルにヌルコードやバックグラウンドにおけるカウントの平均値を代入するといった方法により、当該関心領域以外の部分を削除するものであって良い。
特に、バックグラウンドの平均値を代入するといった方法を用いれば、関心領域境界部における不連続性の影響を抑える事が可能となる。
【００１４】
本発明に係る画像処理方法は、関心領域設定ステップの前に、コンピュータが、投影データセットに含まれる各投影画像において、心筋部位の中心を画面の回転中心に移動させる、画像シフトステップを含み、関心領域設定ステップは、コンピュータが、画像シフトステップ処理後の投影データセットにおける各投影画像上で心筋部位を含む領域に関心領域を設定するといった構成としても良い。
このような構成とすることにより、画像シフトステップ後は、心筋部位が全ての投影画像上で回転中心に位置しているので、代表的な投影画像上で関心領域を設定し、それを全ての投影画像に当てはめる事ができる。これにより、それぞれの投影画像につき別々に関心領域を設定する必要がなく、処理負担を大幅に軽減する事が可能となる。
【００１５】
同様に、本発明に係る画像処理プログラムは、関心領域設定ステップの前に、コンピュータに、投影データセットに含まれる各投影画像において、心筋部位の中心を画面の回転中心に移動させる、画像シフトステップを実行させ、関心領域設定ステップにおいては、コンピュータに、画像シフトステップ処理後の投影データセットにおける各投影画像上で心筋部位を含む領域に関心領域を設定する処理を実行させる構成としても良い。
このような構成とすることにより、画像シフトステップの実行後においては、画像処理方法に係る発明と同様に、心筋部位が全ての投影画像上で回転中心に位置しているので、代表的な投影画像上で関心領域を設定し、それを全ての投影画像に当てはめる事ができる。これにより、それぞれの投影画像につき別々に関心領域を設定する必要がなく、処理負担を大幅に軽減する事が可能となる。
【００１６】
また、本発明に係る画像処理装置は、投影データセットに含まれる各投影画像について、心筋部位の中心を画面の回転中心に移動させる機能を有する、画像シフト部を備えるものであっても良い。
このような機能を具備させることにより、全ての投影画像上で心筋部位を回転中心に位置させることができ、その結果、関心領域設定部は、代表的な投影画像上で関心領域を設定し、それを全ての投影画像に当てはめるといった処理を行う事ができる。これにより、それぞれの投影画像につき別々に関心領域を設定する必要がなく、関心領域設定部における処理負担を大幅に軽減する事が可能となる。
【００１７】
本発明に係る画像処理方法は、コンピュータが、非関心領域削除ステップ後の画像上で、非関心領域削除ステップによって削除された領域を関心領域境界における値で補完する、データ補完ステップをさらに含み、画像再構成ステップは、データ補完ステップ後の投影データセットを用いて画像再構成を行うものであっても良い。
同様に、本発明に係る画像処理プログラムは、コンピュータに、非関心領域削除ステップ後の画像上で、非関心領域削除ステップによって削除された領域を関心領域境界における値で補完する、データ補完ステップをさらに実行させ、画像再構成ステップは、コンピュータに、データ補完ステップ後の投影データセットを用いて画像再構成を実行させるものであるものであっても良い。
同様に、本発明に係る画像処理装置は、非関心領域の削除処理後の画像上で、削除された領域を関心領域境界における値で補完する機能を有する、データ補完部をさらに備えるものであっても良い。
このような構成とすることにより、削除処理後における関心領域境界部の不連続性の影響を抑える事が可能となる。
補完は、例えば、非関心領域の削除処理後の画像につき、関心領域における境界部の値を、縦方向および横方向に延長させる方法を用いて行う事ができる。このとき、縦方向および横方向からともに値が延長されてくる座標におけるピクセルについては、縦または横方向からの何れか一方の値を埋めるか、縦および横方向からの値の平均値で埋める処理を行えばよい。また、縦および横方向の何れからも値が延長されてこない座標におけるピクセルについては、当該ピクセルの上下または左右方向に存在するピクセルに延長されてきた値を埋めるか、それらの値の平均値を埋める処理を行えばよい。
境界部の値を延長する別の方法としては、画像の中心から外に向けて引いた放射状の直線に沿って、境界部の値を延長させる方法を用いても良い。このとき、値の延長されてこない座標におけるピクセルについては、当該ピクセルの上下または左右方向に存在するピクセルに延長されてきた値を埋めるか、それらの値の平均値を埋める処理を行えばよい。なお、上記の構成は、そのまま本発明に係る画像処理装置に適用する事も可能である。
【発明の効果】
【００１８】
本発明により、心筋核医学画像において、心筋部位以外の高集積部位を実質的に除去する事ができ、再構成画像におけるストリークアーチファクトを効果的に減少させる事が可能となった。
【図面の簡単な説明】
【００１９】
【図１】本発明に係る画像処理方法の好ましい態様における処理の流れの一例を示す図
【図２】本発明に係る画像処理方法における画像シフトステップにおける処理の流れの一例を示す図
【図３】本発明に係る画像処理方法におけるシフトパラメータ算出処理の流れの一例を示す図
【図４】本発明に係る画像処理プログラムの好ましい態様における構成の一例を示す図
【図５】本発明に係る画像処理装置の好ましい態様における機能ブロック図の一例を示す図
【図６】心筋部位を画像の回転中心にシフトさせる前の画像（ａ）と、回転中心にシフトさせた後の画像（ｂ）とを表す図（回転中心を丸にて表示）
【図７】マスク処理を実施した例。（ａ）関心領域外に高集積部位が存在している投影画像の例（マスク処理前）、（ｂ）（ａ）の画像につきマスク処理を行った画像、（ｃ）関心領域内に高集積部位が存在している投影画像の例、（ｄ）（ｃ）の画像につき、画像全体にマスク処理を行った画像
【図８】肝臓部に高集積を呈していないファントムにおける再構成画像（ＦＢＰ処理）
【図９】肝臓部に高集積を呈しているファントムにおける再構成画像（ＦＢＰ処理）
【図１０】肝臓部に高集積を呈しているファントムにおける再構成画像（マスク処理後にＦＢＰ処理）
【図１１】心筋部に設定した関心領域を示す図
【発明を実施するための形態】
【００２０】
以下、本発明につき、図面を参照して説明する。なお、以下に示す例は、あくまでも好ましい形態につき説明するものであり、本発明の内容はこれらの記載により何ら限定されるものではない。
【００２１】
まず、本発明に係る画像処理方法について、説明する。図１〜図３は、本発明に係る画像処理方法における、処理の流れを示すフローチャートである。本発明に係る画像処理方法では、まず、投影データセット取得ステップをコンピュータに実行させ、処理に供する投影データセットを本発明を実施するコンピュータシステムに入力する（ステップＳ１）。投影データセットは、ＳＰＥＣＴ装置やＰＥＴ装置といった、通常の核医学画像撮像装置にて撮像されたものを用いる事ができる。投影データセットは、ＤＩＣＯＭ形式等の、コンピュータで読み取り可能な形式で保存されたものを、ハードディスク、ＣＤ−ＲＯＭ、ＤＶＤ等といった、コンピュータで読み取り可能な記憶媒体に格納された状態で取得することができる。取得された投影データセットは、コンピュータシステムに備え付けられた読取装置により読み込まれ、当該システムに入力される。なお、上記投影データセットは、核医学画像撮像装置から、ネットワークを通じて直接入力されるものであっても良い。
【００２２】
投影データの入力が完了したら、画像シフトステップをコンピュータに実行させ、各投影画像における心筋部位の中心を、各投影画像の回転中心に合わせる（ステップＳ２）。この処理は、別に算出したシフトパラメータを各投影画像に適用する事によって、行う事ができる（ステップＳ２１、Ｓ２２）。図３を参照し、このシフトパラメータの算出方法の好ましい例について、説明する。
【００２３】
好ましい態様において、シフトパラメータの算出は、代表スライスの再構成画像を用いて行われる。まず、入力された投影データセットを用い、心筋部位を含む、代表スライスの横断像を一枚作成する（ステップＳ２３）。代表スライスとしては、心筋の中心付近が含まれる横断像を、好ましく用いる事ができる。
【００２４】
代表スライスが得られたら、代表スライス上で、心筋の中心部を決定する（ステップＳ２４）。心筋中心部の決定は、ユーザがディスプレイ等を介して目視にて決定し、マウスポインタ等を用いて指定するといった方法にて行う事ができる。また、代表スライス上で心筋部位の重心を算出し、それを心筋中心部としても良い。
【００２５】
代表スライス上にて心筋中心部を決定したら、画像のシフトパラメータを算出する（ステップＳ２５）。シフトパラメータは、各投影画像において心筋部位を回転中心に移動させるためのシフトパラメータである。シフトパラメータは、上記で決定した心筋中心部の座標と代表スライスの回転中心との関係に基づき、各投影角度の投影画像について、容易に算出する事ができる。このステップにて算出されたシフトパラメータを、画像シフトステップにて用いるシフトパラメータとする。
【００２６】
図１に戻り、本発明に係る画像処理方法について説明する。画像シフトステップが完了したら、コンピュータに関心領域設定ステップを実行させ、心筋部位を包含する関心領域（図１ではＲＯＩと表示）を設定する（ステップＳ３）。関心領域の設定は、投影画像上で、心筋部位を包含する領域の関心領域を設定するといった公知の方法にて行う事ができる。関心領域の形状は特に限定する必要は無いが、画面の回転中心に中心を持つ円や楕円といった既定の形状の関心領域を用いるのがユーザ間の再現性を向上させる上で好ましい。その場合、回転中心に中心を持つ円や楕円形状の関心領域の大きさを、心筋領域を完全に包含するように調整するといった方法を用いる事ができる。また、閾値法等の公知の方法を用いて心筋部位の輪郭を予め抽出し、当該輪郭を構成する各点の回転中心からの距離の最大値を算出した上で、当該最大値より大きな半径の円として、関心領域を求める事も可能である。
また、本実施の形態においては、画像シフトステップを実行させる事により、予め各投影画像における心筋部位を、画像の回転中心にシフトさせる処理が行われている。従って、各投影画像において、心筋部位が、実質的に同じ位置に表示されている。従って、ある代表的な投影画像について、上記の様な方法によって関心領域を求めれば、全ての投影画像に同様の関心領域を適用する事によって、関心領域の設定が完了する。すなわち、各投影画像について別々に関心領域の設定を行う必要がなく、処理負担が大幅に軽減される。
なお、関心領域の大きさは、心筋部位全体を包含したできるだけ小さな領域とすることが好ましいことは、いうまでもない。
【００２７】
関心領域が設定されたら、コンピュータに画像削除ステップを実行させ、心筋部以外の高集積部位が前記関心領域と重なっている投影画像を削除する（ステップＳ４）。画像削除ステップは、それぞれの投影画像を目視にて確認して該当する画像を削除するといった方法により行う事ができるが、サイノグラム上で、心筋部位以外の高集積部位が前記関心領域の輪郭に対応する線（中心をはさんで２本存在）の内側に存在している投影角度を調べ、当該投影角度に対応する投影画像を削除するといった方法を用いても良い。
【００２８】
次に、画像削除ステップ後の投影データセットにつき、非関心領域削除ステップをコンピュータに実行させ、各投影画像上で、前記関心領域として選択された部位以外の部位を、削除する（ステップＳ５）。このステップは、各投影画像について公知のマスク処理を行って前記関心領域以外の領域を削除するといった方法や、当該関心領域以外の領域に対応するピクセルにヌルコードやバックグラウンドにおけるカウントの平均値を代入するといった方法にて行う事ができる。ここで、バックグラウンドにおける平均値は、例えば、投影画像上で信号を発しない部位（肺野等）に一定の領域を指定し、当該領域内にて各ピクセルのカウントの平均値を計算する事によって、容易に求める事ができる。
【００２９】
非関心領域削除ステップが完了したら、画像再構成ステップを実施して、再構成画像を得る（ステップＳ６）。画像再構成は、公知の方法にて行う事ができる。そして、得られた再構成画像を、ディスプレイといった出力機器に出力する（ステップＳ７）。
以上の処理を実行する事により、本発明に係る画像処理方法が完了する。
【００３０】
次に、本発明の別の一側面に係る、画像処理プログラムについて、説明する。図４は、本発明に係る画像処理プログラム２０の最も好ましい態様における構成を、記憶媒体１０とともに示す図である。好ましい態様において、本発明に係る画像処理プログラム２０は、処理を統括するメインモジュール２２と、投影データセット取得モジュール２４と、代表スライス作成モジュール２６と、心筋中心決定モジュール２８と、シフトパラメータ算出モジュール３０と、画像シフトモジュール３２と、関心領域情報入力モジュール３４と、関心領域設定モジュール３６と、画像削除情報入力モジュール３８と、画像削除モジュール４０と、非関心領域削除モジュール４２と、再構成モジュール４４と、出力モジュール４６とにより、構成されている。
好ましい実施態様において、本発明に係る画像処理プログラム２０は、ハードディスク、ＣＤ−ＲＯＭ、ＤＶＤ等といった記憶媒体１０に格納されて提供される。画像処理プログラム２０が格納された記憶媒体１０が、コンピュータに備え付けられた読取装置に挿入される事により、コンピュータが画像処理プログラム２０にアクセス可能となり、後述する画像処理装置２００として動作することが可能となる。なお、画像処理プログラム２０は、記憶媒体を介さず、ネットワークを介して直接提供されるものであっても良い。
【００３１】
投影データセット取得モジュール２４は、コンピュータに、ステップＳ１に係る処理を実行させる。代表スライス作成モジュール２６は、コンピュータに、ステップＳ２３に係る処理を実行させる。心筋中心情報入力モジュール２８は、コンピュータに、ステップＳ２４に係る処理を実行させる。シフトパラメータ算出モジュール３０は、コンピュータに、ステップＳ２５に係る処理を実行させる。画像シフトモジュール３２は、コンピュータに、ステップＳ２に係る処理を実行させる。
【００３２】
関心領域情報入力モジュール３４は、心筋領域を包含するような関心領域に関する情報（形状、大きさ等）に関する情報を受け付ける。関心領域の形状や大きさは、上述した画像処理方法において説明した方法と同様の方法にて決定する事ができる。すなわち、ユーザインターフェースを介し、ユーザが投影画像上で設定した、心筋部位を包含する領域の関心領域を入力するといった方法や、回転中心に中心を持つ円や楕円といった既定の形状を、心筋領域を包含するように大きさを調整して関心領域とするといった方法を用いる事ができる。大きさの調整は、（１）ユーザがディスプレイ等で目視にて確認しながら行うといった方法を用いてもよいが、（２）予め公知の方法を用いて心筋部位の輪郭を抽出し、当該輪郭を構成する各点の回転中心からの距離の最大値を算出した上で、当該最大値より大きな半径の円として求めても良い。また、上記（１）および（２）の方法を組み合わせ、（２）と同様の方法にて自動的に関心領域の大きさを決め、これをディスプレイ等に表示させてユーザが確認を行い、必要に応じて微調整を行うといった方法を用いる事もできる。
【００３３】
関心領域設定モジュール３６は、コンピュータに、ステップＳ３に係る処理を実行させる。画像削除情報入力モジュール３８は、後述する画像削除モジュール４０にて削除すべき投影画像、すなわち、心筋部以外の高集積部位が上記で設定した関心領域内に含まれている画像に関する情報を受け付ける。この処理は、ユーザがディスプレイ上等でそれぞれの投影画像を目視にて確認して該当する画像を指定し、コンピュータがその入力を受け付けるといった方法により行う事ができる。また、サイノグラム上で、心筋部位以外の高集積部位が前記関心領域の輪郭に対応する線（中心をはさんで２本存在）の内側に存在している投影角度を検出し、当該投影角度に対応する投影画像を削除すべき画像として入力するといった方法を用いても良い。画像削除モジュール４０は、コンピュータに、ステップＳ４に係る処理を実行させる。非関心領域削除モジュール４２は、コンピュータに、ステップＳ５に係る処理を実行させる。再構成モジュール４４は、コンピュータに、ステップＳ６に係る処理を実行させる。出力モジュール４６は、コンピュータに、ステップＳ７に係る処理を実行させる。
【００３４】
次に、本発明のさらに別の一側面に係る、画像処理装置について説明する。図５は、本発明に係る画像処理装置２００の、好ましい態様における機能ブロック図である。好ましい態様において、本発明に係る画像処理装置２００は、ＳＰＥＣＴやＰＥＴといった、核医学画像撮像装置１００と、有機的に結合されている。好ましい態様において、画像処理装置２００は、機能的に、投影データセット取得部２２０と、代表スライス作成部３２０と、パラメータ算出部３４０と、入力部３６０と、画像シフト部２４０と、関心領域設定部（図５では、ＲＯＩ設定部と表示）２６０と、判断部３８０と、画像削除部２８０と、非関心領域削除部３００と、再構成部４００と、出力部４２０とにより構成されている。
【００３５】
本発明に係る画像処理装置２００は、典型的には、前記画像処理プログラム２０を読み込んだコンピュータとして構成されるが、それぞれの機能を機械的に実行する機能を有した部品の組み合わせによって構成する事も可能である。
【００３６】
投影データセット取得部２２０は、ステップＳ１に係る処理を実行する。代表スライス作成部３２０は、ステップＳ２３に係る処理を実行する。パラメータ算出部３４０は、入力部３６０を介して入力された心筋中心に関する情報を用い、ステップＳ２５に係る処理を実行する。画像シフト部２４０は、ステップＳ２に関する処理を実行する。関心領域設定部（図５では、ＲＯＩ設定部と表示）２６０は、入力部３６０より入力された関心領域の形状や大きさに関する情報を用い、ステップＳ３に関する処理を実行する。判断部３８０は、入力部３６０を介して入力された情報を用い、各投影画像が画像削除部２８０において削除すべき投影画像か否かを判断する。例えば、サイノグラム上で心筋部位以外の高集績部位が前記関心領域の輪郭に対応する線（中心をはさんで２本存在）の内側に存在しているか否かを確認し、内側に存在している投影角度における投影画像を、削除すべき投影画像と判断する。画像削除部２８０は、ステップＳ４に係る処理を実行する。非関心領域削除部３００は、ステップＳ５に係る処理を実行する。再構成部４００は、ステップＳ６に係る処理を実行する。出力部４２０は、ステップＳ７に係る処理を実行する。
【実施例】
【００３７】
以下に、実施例を用いて本発明をさらに詳しく説明するが、本発明の内容はこの内容になんら限定されるものではない。
【００３８】
（心筋ファントムを用いた検証）
心筋ファントム（ＲＨ２型、株式会社京都科学製、心筋下壁部に人工欠損有り）を用い、心筋部にヘキサキス（２−メトキシイソブチルイソニトリル）テクネチウム（９９ｍＴｃ）を１５ＭＢｑとなるように注入した。
また、肝臓に高集積を有する症例を再現するため、上記の心筋ファントムに、さらに肝臓部に心筋部の１０倍の放射能濃度となるようにヘキサキス（２−メトキシイソブチルイソニトリル）テクネチウム（９９ｍＴｃ）を注入したファントムも用意した。
作成した各ファントムにつき、ＳＰＥＣＴ装置を用いた撮像を行った。ＳＰＥＣＴ装置並びに主な撮像条件は、表１の通りである。
【００３９】
【表１】
【００４０】
得られた投影データセットにつき、通常のＦＢＰ処理を用いて画像再構成を行った。また、肝臓部に高集積部位を設けたファントムについては、下記の方法による処理を行った。
（１）各投影画像につき、心筋部位が画像の回転中心にくるようにシフトさせた（図６）。
（２）心筋部位に楕円の関心領域を設定した。
（３）この関心領域内に肝臓部（高集積部）が含まれている画像については、画像全体にマスク処理を施し、画像全体の削除を行った。一方、それ以外の画像については、関心領域外にマスク処理を施し、肝臓部（高集積部）を含む心筋外の領域の削除を行った。（図７）
（４）マスク処理後の投影データセットを用い、表１記載の方法によるＦＢＰ処理を用いた画像再構成を行った。
【００４１】
得られた再構成画像の一例を、図８〜図１０に示す。肝臓部に高集積を有していないファントムでは、通常のＦＢＰ処理によって得られた再構成画像において、ストリークアーチファクトのない良好な画像が得られ、下壁部の人口欠損も良好に描出されていた（図８）。
肝臓部に高集積を呈したファントムにおいては、通常のＦＢＰ処理により得られた再構成画像においてストリークアーチファクトが発生し、その影響により下壁部の人口欠損も描出されていなかった（図９）。
一方、肝臓部に高集積を呈したファントムにおいて、通常のＦＢＰ処理に先立ちマスク処理を施した例では、ストリークアーチファクトが消失し、下壁部の人口欠損も明瞭に描出されていた（図１０）。
【００４２】
次に、肝臓部に高集積を呈したファントムにて得られた画像につき、同一断面における横断像上で心筋部位にドーナツ状の関心領域を設定し（図１１）、％アップテークを求めてばらつきの比較を行った。その結果、通常のＦＢＰ処理を施した再構成画像では分散が９７．４５であったのに対し、ＦＢＰ処理に先立ちマスク処理を施した再構成画像では分散が４８．５４と、減少していた。
【産業上の利用可能性】
【００４３】
本発明は、核医学画像の画像処理プログラム並びに核医学画像処理装置の分野において、利用する事ができる。
【符号の説明】
【００４４】
２０ 本発明に係る画像処理プログラム
１０ 記憶媒体
２２ メインモジュール
２４ 投影データセット取得モジュール
２６ 代表スライス作成モジュール
２８ 心筋中心決定モジュール
３０ シフトパラメータ算出モジュール
３２ 画像シフトモジュール
３４ 関心領域情報入力モジュール
３６ 関心領域設定モジュール
３８ 画像削除情報入力モジュール
４０ 画像削除モジュール
４２ 非関心領域削除モジュール
４４ 再構成モジュール
４６ 出力モジュール
１００ 核医学画像撮像装置
２００ 画像処理装置
２２０ 投影データセット取得部
２４０ 画像シフト部
２６０ 関心領域設定部（図５では、ＲＯＩ設定部と表示）
２８０ 画像削除部
３００ 非関心領域削除部
３２０ 代表スライス作成部
３４０ パラメータ算出部
３６０ 入力部
３８０ 判断部
４００ 再構成部
４２０ 出力部
【技術分野】
【０００１】
この発明は、核医学画像の画像処理技術に関する。より詳しくは、心筋核医学画像において、アーチファクトを有効に除去し得る、画像処理技術に関する。
【背景技術】
【０００２】
心臓疾患の画像診断の一つとして、ＳＰＥＣＴに代表される、核医学画像診断法が用いられている。核医学画像診断法では、放射性同位元素で標識された薬剤（放射性医薬品）を被験者に投与し、放出されたγ線を検出して画像化処理を行う事により、診断画像を得る事ができる。心臓疾患の画像診断に用いる放射性医薬品としては、塩化タリウム−２０１、ヘキサキス（２−メトキシイソブチルイソニトリル）テクネチウム（９９ｍＴｃ）、テトロホスミンテクネチウム（９９ｍＴｃ）、３−ヨードベンジルグアニジン（１２３Ｉ）、１５−（４−ヨードフェニル）−３（Ｒ，Ｓ）−メチルデカン酸（１２３Ｉ）等が開発され、臨床において用いられている。
【０００３】
これらの放射性医薬品のうち、放射性同位元素として９９ｍＴｃを用いた製剤（以下、９９ｍＴｃ製剤という）は、キット化できるために汎用性に優れ、また、ＳＰＥＣＴ装置に適したエネルギーを有しているために良好な画像を得る事ができるといった利点を有している。しかし、９９ｍＴｃ製剤は、一般に分子サイズが大きく、肝胆道系に高い集積を有する。このような心筋以外における高集積は、再構成画像上でストリークアーチファクトの原因となり、診断の妨げになるといった問題がある（非特許文献１、非特許文献２、非特許文献３、非特許文献４）。また、放射性ヨード製剤においても、肝臓に高集積し、心筋ＳＰＥＣＴにおいて下壁部が欠損するといったことが報告されている（非特許文献５）。
【０００４】
心筋以外への高集積に起因する再構成画像上でのアーチファクトの発生を防ぐ方法として、マスク処理によって高集積部位を除去する方法がこれまでに検討されている（非特許文献１、非特許文献２、非特許文献３）。これまでの報告によれば、再構成処理前の投影データにマスク処理を施す事により、当該高集積による心筋画像への影響を減少させ得ることが開示されている（非特許文献１、非特許文献２、非特許文献３）。しかし、心筋以外の高集積部位が心筋に近接していた場合には、画像データ上で当該高集積部位を除去する事は難しく、心筋画像への当該高集積部位の影響を完全に取り除く事は困難である（非特許文献２、非特許文献３）。
【先行技術文献】
【０００５】
【非特許文献１】 久野晴丘 他、「９９ｍＴｃ−ｔｅｔｒｏｆｏｓｍｉｎによる胆嚢の高集積が心筋ＳＰＥＣＴ画像に及ぼす影響と改善法についての検討」、日本放射線技術学会誌、２０００年８月、５６、８、ｐ．１０４４−１０５１
【非特許文献２】 今井嘉門 他、「ＭＩＢＩの心筋シンチのアーティファクトの原因に関して」、核医学、１９９５、３２、３、ｐ．３０７−３１０
【非特許文献３】 高木昭浩 他、「９９ｍＴｃ心筋血流製剤を用いたＳＰＥＣＴ撮像における肝の高集積が心筋に及ぼす影響の軽減；マスク処理法の有用性と問題点」、核医学、１９９９、３６、５、ｐ．４５９−４６５
【非特許文献４】 大西英雄 他、「９９ｍＴｃ標識心筋血流製剤を用いた心筋ＳＰＥＣＴ画像への肝臓からの影響」、核医学、１９９８、３５、６、ｐ．３７５−３８３
【非特許文献５】 小林秀樹 他、「肝高集積が原因となって出現する１２３Ｉ−ＭＩＢＧ心筋ＳＰＥＣＴの下壁欠損像の特徴とその対策−ファントムを用いた検討−」、核医学、１９９４、３１、４、ｐ．３５９−３６６
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明は上記事実に鑑みてなされたものであり、心筋以外の高集積の心筋画像への影響を除去する方法、当該方法をコンピュータに実行させるプログラム、および当該方法を実行するための装置を提供する事を目的とした。
【課題を解決するための手段】
【０００７】
発明者は検討の結果、投影データセットにおいて、心筋以外の高集積部位が心筋部位に近接している投影画像を削除しても、残りのデータを用いて作成された再構成画像の画質に本質的な悪影響を与えない事を見出した。この知見に基づき、各投影画像上で心筋全体に関心領域を設定し、当該関心領域内に心筋部以外の高集積部位が存在している投影画像を削除し、残りの投影画像につき前記関心領域以外の部位（すなわち、心筋領域に設定した関心領域以外の部位）をマスク処理等を施して除去するといった処理を行う事により、上記の問題点が解消し、当該高集積部位の影響を効果的に除去しうる事を見出し、本発明を完成させた。
【０００８】
なお、本明細書において、投影データセットとは、一度の撮像処理により得られた各投影角度での投影画像のセットをいう。例えば、患者に対し、２°ステップで３６０°の投影を行った場合、投影角度２°、４°、６°、…、３５８°、３６０°において得られた１８０枚の投影画像のセットが、本明細書における投影データセットとなる。
【０００９】
本発明の一側面に係る画像処理方法は、コンピュータが心筋核医学画像の投影データセットを取得する投影データセット取得ステップと、コンピュータが、前記投影データセットに含まれる各投影画像上で心筋部位を含む領域に関心領域を設定する関心領域設定ステップと、コンピュータが、心筋部以外の高集積部位が前記関心領域と重なっている投影画像を削除する画像削除ステップと、コンピュータが、前記画像削除ステップにて心筋部以外の高集積部位が前記関心領域内と重なっている投影画像を削除した後の投影データセットに含まれる各投影画像につき、前記関心領域設定ステップにて設定した関心領域以外の部分を削除する非関心領域削除ステップと、コンピュータが、非関心領域削除ステップ後の投影データセットを用いて画像再構成を行う画像再構成ステップと、を含む。
【００１０】
本発明の別の一側面に係る画像処理プログラムは、コンピュータに、心筋核医学画像の投影データセットを取得する投影データ取得ステップと、前記投影データセットに含まれる各投影画像上で心筋部位を含む領域に関心領域を設定する関心領域設定ステップと、心筋部以外の高集積部位が前記関心領域と重なっている投影画像を削除する画像削除ステップと、前記画像削除ステップにて心筋部以外の高集積部位が前記関心領域内と重なっている投影画像を削除した後の投影データセットに含まれる各投影画像につき、前記関心領域設定ステップにて設定した関心領域以外の部分を削除する非関心領域削除ステップと、非関心領域削除ステップ後の投影データセットを用いて画像再構成を行う画像再構成ステップと、を実行させる。
【００１１】
本発明のさらに別の一側面に係る画像処理装置は、心筋核医学画像の投影データセットを取り込む投影データセット取得部と、前記投影データセットに含まれる各投影画像上で心筋部位を含む領域に関心領域を設定する関心領域設定部と、前記投影データセットに含まれる各投影画像のうち、心筋部位以外の高集積部位が前記関心領域と重なっている投影画像を削除する画像削除部と、前記関心領域設定部の機能により設定された前記関心領域以外の部分を各投影画像上において削除する非関心領域削除部と、投影データセットを用いて画像再構成を行う画像再構成部と、を備える。
【００１２】
従来技術において指摘されてきたように、投影画像上において心筋部位と心筋以外の高集積部位とが接近している場合には、当該高集積部位を完全に削除する事が難しく、このことが、再構成画像上で当該高集積部位の影響を完全に除去することを困難にしてきた（今井嘉門 他、核医学、１９９５、３２、３、ｐ．３０７−３１０；高木昭浩 他、核医学、１９９９、３６、５、ｐ．４５９−４６５）。本発明に係る画像処理方法および画像処理プログラムでは、心筋部位を包含するよう設定した関心領域と心筋部位以外の高集積とが重なった投影画像を画像ごと削除し、それ以外の投影画像では前記関心領域以外の領域を削除する処理を行っているので、再構成画像において、心筋以外の高集積部位の影響を、効果的に取り除く事ができる。
また同様に、本発明に係る画像処理装置も、心筋部位を包含するよう設定した関心領域と心筋部位以外の高集積とが重なった投影画像を画像ごと削除し、それ以外の投影画像では前記関心領域以外の領域を削除する機能を備えているので、心筋部以外の高集積部位からの影響を完全に取り除いた再構成画像を生成する事が可能となる。
【００１３】
本発明に係る画像処理方法および画像処理プログラムにおいて、心筋部位に設定した関心領域以外の部分を削除する非関心領域削除ステップは、当該関心領域以外の部分をマスク処理によって削除するといった方法や、当該関心領域以外の領域に対応するピクセルにヌルコードやバックグラウンドにおけるカウントの平均値を代入するといった方法にて行う事ができる。
同様に、本発明に係る画像処理装置における非関心領域削除部は、心筋部位に設定した関心領域以外の部分をマスク処理で削除するといった方法や、当該関心領域以外の領域に対応するピクセルにヌルコードやバックグラウンドにおけるカウントの平均値を代入するといった方法により、当該関心領域以外の部分を削除するものであって良い。
特に、バックグラウンドの平均値を代入するといった方法を用いれば、関心領域境界部における不連続性の影響を抑える事が可能となる。
【００１４】
本発明に係る画像処理方法は、関心領域設定ステップの前に、コンピュータが、投影データセットに含まれる各投影画像において、心筋部位の中心を画面の回転中心に移動させる、画像シフトステップを含み、関心領域設定ステップは、コンピュータが、画像シフトステップ処理後の投影データセットにおける各投影画像上で心筋部位を含む領域に関心領域を設定するといった構成としても良い。
このような構成とすることにより、画像シフトステップ後は、心筋部位が全ての投影画像上で回転中心に位置しているので、代表的な投影画像上で関心領域を設定し、それを全ての投影画像に当てはめる事ができる。これにより、それぞれの投影画像につき別々に関心領域を設定する必要がなく、処理負担を大幅に軽減する事が可能となる。
【００１５】
同様に、本発明に係る画像処理プログラムは、関心領域設定ステップの前に、コンピュータに、投影データセットに含まれる各投影画像において、心筋部位の中心を画面の回転中心に移動させる、画像シフトステップを実行させ、関心領域設定ステップにおいては、コンピュータに、画像シフトステップ処理後の投影データセットにおける各投影画像上で心筋部位を含む領域に関心領域を設定する処理を実行させる構成としても良い。
このような構成とすることにより、画像シフトステップの実行後においては、画像処理方法に係る発明と同様に、心筋部位が全ての投影画像上で回転中心に位置しているので、代表的な投影画像上で関心領域を設定し、それを全ての投影画像に当てはめる事ができる。これにより、それぞれの投影画像につき別々に関心領域を設定する必要がなく、処理負担を大幅に軽減する事が可能となる。
【００１６】
また、本発明に係る画像処理装置は、投影データセットに含まれる各投影画像について、心筋部位の中心を画面の回転中心に移動させる機能を有する、画像シフト部を備えるものであっても良い。
このような機能を具備させることにより、全ての投影画像上で心筋部位を回転中心に位置させることができ、その結果、関心領域設定部は、代表的な投影画像上で関心領域を設定し、それを全ての投影画像に当てはめるといった処理を行う事ができる。これにより、それぞれの投影画像につき別々に関心領域を設定する必要がなく、関心領域設定部における処理負担を大幅に軽減する事が可能となる。
【００１７】
本発明に係る画像処理方法は、コンピュータが、非関心領域削除ステップ後の画像上で、非関心領域削除ステップによって削除された領域を関心領域境界における値で補完する、データ補完ステップをさらに含み、画像再構成ステップは、データ補完ステップ後の投影データセットを用いて画像再構成を行うものであっても良い。
同様に、本発明に係る画像処理プログラムは、コンピュータに、非関心領域削除ステップ後の画像上で、非関心領域削除ステップによって削除された領域を関心領域境界における値で補完する、データ補完ステップをさらに実行させ、画像再構成ステップは、コンピュータに、データ補完ステップ後の投影データセットを用いて画像再構成を実行させるものであるものであっても良い。
同様に、本発明に係る画像処理装置は、非関心領域の削除処理後の画像上で、削除された領域を関心領域境界における値で補完する機能を有する、データ補完部をさらに備えるものであっても良い。
このような構成とすることにより、削除処理後における関心領域境界部の不連続性の影響を抑える事が可能となる。
補完は、例えば、非関心領域の削除処理後の画像につき、関心領域における境界部の値を、縦方向および横方向に延長させる方法を用いて行う事ができる。このとき、縦方向および横方向からともに値が延長されてくる座標におけるピクセルについては、縦または横方向からの何れか一方の値を埋めるか、縦および横方向からの値の平均値で埋める処理を行えばよい。また、縦および横方向の何れからも値が延長されてこない座標におけるピクセルについては、当該ピクセルの上下または左右方向に存在するピクセルに延長されてきた値を埋めるか、それらの値の平均値を埋める処理を行えばよい。
境界部の値を延長する別の方法としては、画像の中心から外に向けて引いた放射状の直線に沿って、境界部の値を延長させる方法を用いても良い。このとき、値の延長されてこない座標におけるピクセルについては、当該ピクセルの上下または左右方向に存在するピクセルに延長されてきた値を埋めるか、それらの値の平均値を埋める処理を行えばよい。なお、上記の構成は、そのまま本発明に係る画像処理装置に適用する事も可能である。
【発明の効果】
【００１８】
本発明により、心筋核医学画像において、心筋部位以外の高集積部位を実質的に除去する事ができ、再構成画像におけるストリークアーチファクトを効果的に減少させる事が可能となった。
【図面の簡単な説明】
【００１９】
【図１】本発明に係る画像処理方法の好ましい態様における処理の流れの一例を示す図
【図２】本発明に係る画像処理方法における画像シフトステップにおける処理の流れの一例を示す図
【図３】本発明に係る画像処理方法におけるシフトパラメータ算出処理の流れの一例を示す図
【図４】本発明に係る画像処理プログラムの好ましい態様における構成の一例を示す図
【図５】本発明に係る画像処理装置の好ましい態様における機能ブロック図の一例を示す図
【図６】心筋部位を画像の回転中心にシフトさせる前の画像（ａ）と、回転中心にシフトさせた後の画像（ｂ）とを表す図（回転中心を丸にて表示）
【図７】マスク処理を実施した例。（ａ）関心領域外に高集積部位が存在している投影画像の例（マスク処理前）、（ｂ）（ａ）の画像につきマスク処理を行った画像、（ｃ）関心領域内に高集積部位が存在している投影画像の例、（ｄ）（ｃ）の画像につき、画像全体にマスク処理を行った画像
【図８】肝臓部に高集積を呈していないファントムにおける再構成画像（ＦＢＰ処理）
【図９】肝臓部に高集積を呈しているファントムにおける再構成画像（ＦＢＰ処理）
【図１０】肝臓部に高集積を呈しているファントムにおける再構成画像（マスク処理後にＦＢＰ処理）
【図１１】心筋部に設定した関心領域を示す図
【発明を実施するための形態】
【００２０】
以下、本発明につき、図面を参照して説明する。なお、以下に示す例は、あくまでも好ましい形態につき説明するものであり、本発明の内容はこれらの記載により何ら限定されるものではない。
【００２１】
まず、本発明に係る画像処理方法について、説明する。図１〜図３は、本発明に係る画像処理方法における、処理の流れを示すフローチャートである。本発明に係る画像処理方法では、まず、投影データセット取得ステップをコンピュータに実行させ、処理に供する投影データセットを本発明を実施するコンピュータシステムに入力する（ステップＳ１）。投影データセットは、ＳＰＥＣＴ装置やＰＥＴ装置といった、通常の核医学画像撮像装置にて撮像されたものを用いる事ができる。投影データセットは、ＤＩＣＯＭ形式等の、コンピュータで読み取り可能な形式で保存されたものを、ハードディスク、ＣＤ−ＲＯＭ、ＤＶＤ等といった、コンピュータで読み取り可能な記憶媒体に格納された状態で取得することができる。取得された投影データセットは、コンピュータシステムに備え付けられた読取装置により読み込まれ、当該システムに入力される。なお、上記投影データセットは、核医学画像撮像装置から、ネットワークを通じて直接入力されるものであっても良い。
【００２２】
投影データの入力が完了したら、画像シフトステップをコンピュータに実行させ、各投影画像における心筋部位の中心を、各投影画像の回転中心に合わせる（ステップＳ２）。この処理は、別に算出したシフトパラメータを各投影画像に適用する事によって、行う事ができる（ステップＳ２１、Ｓ２２）。図３を参照し、このシフトパラメータの算出方法の好ましい例について、説明する。
【００２３】
好ましい態様において、シフトパラメータの算出は、代表スライスの再構成画像を用いて行われる。まず、入力された投影データセットを用い、心筋部位を含む、代表スライスの横断像を一枚作成する（ステップＳ２３）。代表スライスとしては、心筋の中心付近が含まれる横断像を、好ましく用いる事ができる。
【００２４】
代表スライスが得られたら、代表スライス上で、心筋の中心部を決定する（ステップＳ２４）。心筋中心部の決定は、ユーザがディスプレイ等を介して目視にて決定し、マウスポインタ等を用いて指定するといった方法にて行う事ができる。また、代表スライス上で心筋部位の重心を算出し、それを心筋中心部としても良い。
【００２５】
代表スライス上にて心筋中心部を決定したら、画像のシフトパラメータを算出する（ステップＳ２５）。シフトパラメータは、各投影画像において心筋部位を回転中心に移動させるためのシフトパラメータである。シフトパラメータは、上記で決定した心筋中心部の座標と代表スライスの回転中心との関係に基づき、各投影角度の投影画像について、容易に算出する事ができる。このステップにて算出されたシフトパラメータを、画像シフトステップにて用いるシフトパラメータとする。
【００２６】
図１に戻り、本発明に係る画像処理方法について説明する。画像シフトステップが完了したら、コンピュータに関心領域設定ステップを実行させ、心筋部位を包含する関心領域（図１ではＲＯＩと表示）を設定する（ステップＳ３）。関心領域の設定は、投影画像上で、心筋部位を包含する領域の関心領域を設定するといった公知の方法にて行う事ができる。関心領域の形状は特に限定する必要は無いが、画面の回転中心に中心を持つ円や楕円といった既定の形状の関心領域を用いるのがユーザ間の再現性を向上させる上で好ましい。その場合、回転中心に中心を持つ円や楕円形状の関心領域の大きさを、心筋領域を完全に包含するように調整するといった方法を用いる事ができる。また、閾値法等の公知の方法を用いて心筋部位の輪郭を予め抽出し、当該輪郭を構成する各点の回転中心からの距離の最大値を算出した上で、当該最大値より大きな半径の円として、関心領域を求める事も可能である。
また、本実施の形態においては、画像シフトステップを実行させる事により、予め各投影画像における心筋部位を、画像の回転中心にシフトさせる処理が行われている。従って、各投影画像において、心筋部位が、実質的に同じ位置に表示されている。従って、ある代表的な投影画像について、上記の様な方法によって関心領域を求めれば、全ての投影画像に同様の関心領域を適用する事によって、関心領域の設定が完了する。すなわち、各投影画像について別々に関心領域の設定を行う必要がなく、処理負担が大幅に軽減される。
なお、関心領域の大きさは、心筋部位全体を包含したできるだけ小さな領域とすることが好ましいことは、いうまでもない。
【００２７】
関心領域が設定されたら、コンピュータに画像削除ステップを実行させ、心筋部以外の高集積部位が前記関心領域と重なっている投影画像を削除する（ステップＳ４）。画像削除ステップは、それぞれの投影画像を目視にて確認して該当する画像を削除するといった方法により行う事ができるが、サイノグラム上で、心筋部位以外の高集積部位が前記関心領域の輪郭に対応する線（中心をはさんで２本存在）の内側に存在している投影角度を調べ、当該投影角度に対応する投影画像を削除するといった方法を用いても良い。
【００２８】
次に、画像削除ステップ後の投影データセットにつき、非関心領域削除ステップをコンピュータに実行させ、各投影画像上で、前記関心領域として選択された部位以外の部位を、削除する（ステップＳ５）。このステップは、各投影画像について公知のマスク処理を行って前記関心領域以外の領域を削除するといった方法や、当該関心領域以外の領域に対応するピクセルにヌルコードやバックグラウンドにおけるカウントの平均値を代入するといった方法にて行う事ができる。ここで、バックグラウンドにおける平均値は、例えば、投影画像上で信号を発しない部位（肺野等）に一定の領域を指定し、当該領域内にて各ピクセルのカウントの平均値を計算する事によって、容易に求める事ができる。
【００２９】
非関心領域削除ステップが完了したら、画像再構成ステップを実施して、再構成画像を得る（ステップＳ６）。画像再構成は、公知の方法にて行う事ができる。そして、得られた再構成画像を、ディスプレイといった出力機器に出力する（ステップＳ７）。
以上の処理を実行する事により、本発明に係る画像処理方法が完了する。
【００３０】
次に、本発明の別の一側面に係る、画像処理プログラムについて、説明する。図４は、本発明に係る画像処理プログラム２０の最も好ましい態様における構成を、記憶媒体１０とともに示す図である。好ましい態様において、本発明に係る画像処理プログラム２０は、処理を統括するメインモジュール２２と、投影データセット取得モジュール２４と、代表スライス作成モジュール２６と、心筋中心決定モジュール２８と、シフトパラメータ算出モジュール３０と、画像シフトモジュール３２と、関心領域情報入力モジュール３４と、関心領域設定モジュール３６と、画像削除情報入力モジュール３８と、画像削除モジュール４０と、非関心領域削除モジュール４２と、再構成モジュール４４と、出力モジュール４６とにより、構成されている。
好ましい実施態様において、本発明に係る画像処理プログラム２０は、ハードディスク、ＣＤ−ＲＯＭ、ＤＶＤ等といった記憶媒体１０に格納されて提供される。画像処理プログラム２０が格納された記憶媒体１０が、コンピュータに備え付けられた読取装置に挿入される事により、コンピュータが画像処理プログラム２０にアクセス可能となり、後述する画像処理装置２００として動作することが可能となる。なお、画像処理プログラム２０は、記憶媒体を介さず、ネットワークを介して直接提供されるものであっても良い。
【００３１】
投影データセット取得モジュール２４は、コンピュータに、ステップＳ１に係る処理を実行させる。代表スライス作成モジュール２６は、コンピュータに、ステップＳ２３に係る処理を実行させる。心筋中心情報入力モジュール２８は、コンピュータに、ステップＳ２４に係る処理を実行させる。シフトパラメータ算出モジュール３０は、コンピュータに、ステップＳ２５に係る処理を実行させる。画像シフトモジュール３２は、コンピュータに、ステップＳ２に係る処理を実行させる。
【００３２】
関心領域情報入力モジュール３４は、心筋領域を包含するような関心領域に関する情報（形状、大きさ等）に関する情報を受け付ける。関心領域の形状や大きさは、上述した画像処理方法において説明した方法と同様の方法にて決定する事ができる。すなわち、ユーザインターフェースを介し、ユーザが投影画像上で設定した、心筋部位を包含する領域の関心領域を入力するといった方法や、回転中心に中心を持つ円や楕円といった既定の形状を、心筋領域を包含するように大きさを調整して関心領域とするといった方法を用いる事ができる。大きさの調整は、（１）ユーザがディスプレイ等で目視にて確認しながら行うといった方法を用いてもよいが、（２）予め公知の方法を用いて心筋部位の輪郭を抽出し、当該輪郭を構成する各点の回転中心からの距離の最大値を算出した上で、当該最大値より大きな半径の円として求めても良い。また、上記（１）および（２）の方法を組み合わせ、（２）と同様の方法にて自動的に関心領域の大きさを決め、これをディスプレイ等に表示させてユーザが確認を行い、必要に応じて微調整を行うといった方法を用いる事もできる。
【００３３】
関心領域設定モジュール３６は、コンピュータに、ステップＳ３に係る処理を実行させる。画像削除情報入力モジュール３８は、後述する画像削除モジュール４０にて削除すべき投影画像、すなわち、心筋部以外の高集積部位が上記で設定した関心領域内に含まれている画像に関する情報を受け付ける。この処理は、ユーザがディスプレイ上等でそれぞれの投影画像を目視にて確認して該当する画像を指定し、コンピュータがその入力を受け付けるといった方法により行う事ができる。また、サイノグラム上で、心筋部位以外の高集積部位が前記関心領域の輪郭に対応する線（中心をはさんで２本存在）の内側に存在している投影角度を検出し、当該投影角度に対応する投影画像を削除すべき画像として入力するといった方法を用いても良い。画像削除モジュール４０は、コンピュータに、ステップＳ４に係る処理を実行させる。非関心領域削除モジュール４２は、コンピュータに、ステップＳ５に係る処理を実行させる。再構成モジュール４４は、コンピュータに、ステップＳ６に係る処理を実行させる。出力モジュール４６は、コンピュータに、ステップＳ７に係る処理を実行させる。
【００３４】
次に、本発明のさらに別の一側面に係る、画像処理装置について説明する。図５は、本発明に係る画像処理装置２００の、好ましい態様における機能ブロック図である。好ましい態様において、本発明に係る画像処理装置２００は、ＳＰＥＣＴやＰＥＴといった、核医学画像撮像装置１００と、有機的に結合されている。好ましい態様において、画像処理装置２００は、機能的に、投影データセット取得部２２０と、代表スライス作成部３２０と、パラメータ算出部３４０と、入力部３６０と、画像シフト部２４０と、関心領域設定部（図５では、ＲＯＩ設定部と表示）２６０と、判断部３８０と、画像削除部２８０と、非関心領域削除部３００と、再構成部４００と、出力部４２０とにより構成されている。
【００３５】
本発明に係る画像処理装置２００は、典型的には、前記画像処理プログラム２０を読み込んだコンピュータとして構成されるが、それぞれの機能を機械的に実行する機能を有した部品の組み合わせによって構成する事も可能である。
【００３６】
投影データセット取得部２２０は、ステップＳ１に係る処理を実行する。代表スライス作成部３２０は、ステップＳ２３に係る処理を実行する。パラメータ算出部３４０は、入力部３６０を介して入力された心筋中心に関する情報を用い、ステップＳ２５に係る処理を実行する。画像シフト部２４０は、ステップＳ２に関する処理を実行する。関心領域設定部（図５では、ＲＯＩ設定部と表示）２６０は、入力部３６０より入力された関心領域の形状や大きさに関する情報を用い、ステップＳ３に関する処理を実行する。判断部３８０は、入力部３６０を介して入力された情報を用い、各投影画像が画像削除部２８０において削除すべき投影画像か否かを判断する。例えば、サイノグラム上で心筋部位以外の高集績部位が前記関心領域の輪郭に対応する線（中心をはさんで２本存在）の内側に存在しているか否かを確認し、内側に存在している投影角度における投影画像を、削除すべき投影画像と判断する。画像削除部２８０は、ステップＳ４に係る処理を実行する。非関心領域削除部３００は、ステップＳ５に係る処理を実行する。再構成部４００は、ステップＳ６に係る処理を実行する。出力部４２０は、ステップＳ７に係る処理を実行する。
【実施例】
【００３７】
以下に、実施例を用いて本発明をさらに詳しく説明するが、本発明の内容はこの内容になんら限定されるものではない。
【００３８】
（心筋ファントムを用いた検証）
心筋ファントム（ＲＨ２型、株式会社京都科学製、心筋下壁部に人工欠損有り）を用い、心筋部にヘキサキス（２−メトキシイソブチルイソニトリル）テクネチウム（９９ｍＴｃ）を１５ＭＢｑとなるように注入した。
また、肝臓に高集積を有する症例を再現するため、上記の心筋ファントムに、さらに肝臓部に心筋部の１０倍の放射能濃度となるようにヘキサキス（２−メトキシイソブチルイソニトリル）テクネチウム（９９ｍＴｃ）を注入したファントムも用意した。
作成した各ファントムにつき、ＳＰＥＣＴ装置を用いた撮像を行った。ＳＰＥＣＴ装置並びに主な撮像条件は、表１の通りである。
【００３９】
【表１】
【００４０】
得られた投影データセットにつき、通常のＦＢＰ処理を用いて画像再構成を行った。また、肝臓部に高集積部位を設けたファントムについては、下記の方法による処理を行った。
（１）各投影画像につき、心筋部位が画像の回転中心にくるようにシフトさせた（図６）。
（２）心筋部位に楕円の関心領域を設定した。
（３）この関心領域内に肝臓部（高集積部）が含まれている画像については、画像全体にマスク処理を施し、画像全体の削除を行った。一方、それ以外の画像については、関心領域外にマスク処理を施し、肝臓部（高集積部）を含む心筋外の領域の削除を行った。（図７）
（４）マスク処理後の投影データセットを用い、表１記載の方法によるＦＢＰ処理を用いた画像再構成を行った。
【００４１】
得られた再構成画像の一例を、図８〜図１０に示す。肝臓部に高集積を有していないファントムでは、通常のＦＢＰ処理によって得られた再構成画像において、ストリークアーチファクトのない良好な画像が得られ、下壁部の人口欠損も良好に描出されていた（図８）。
肝臓部に高集積を呈したファントムにおいては、通常のＦＢＰ処理により得られた再構成画像においてストリークアーチファクトが発生し、その影響により下壁部の人口欠損も描出されていなかった（図９）。
一方、肝臓部に高集積を呈したファントムにおいて、通常のＦＢＰ処理に先立ちマスク処理を施した例では、ストリークアーチファクトが消失し、下壁部の人口欠損も明瞭に描出されていた（図１０）。
【００４２】
次に、肝臓部に高集積を呈したファントムにて得られた画像につき、同一断面における横断像上で心筋部位にドーナツ状の関心領域を設定し（図１１）、％アップテークを求めてばらつきの比較を行った。その結果、通常のＦＢＰ処理を施した再構成画像では分散が９７．４５であったのに対し、ＦＢＰ処理に先立ちマスク処理を施した再構成画像では分散が４８．５４と、減少していた。
【産業上の利用可能性】
【００４３】
本発明は、核医学画像の画像処理プログラム並びに核医学画像処理装置の分野において、利用する事ができる。
【符号の説明】
【００４４】
２０ 本発明に係る画像処理プログラム
１０ 記憶媒体
２２ メインモジュール
２４ 投影データセット取得モジュール
２６ 代表スライス作成モジュール
２８ 心筋中心決定モジュール
３０ シフトパラメータ算出モジュール
３２ 画像シフトモジュール
３４ 関心領域情報入力モジュール
３６ 関心領域設定モジュール
３８ 画像削除情報入力モジュール
４０ 画像削除モジュール
４２ 非関心領域削除モジュール
４４ 再構成モジュール
４６ 出力モジュール
１００ 核医学画像撮像装置
２００ 画像処理装置
２２０ 投影データセット取得部
２４０ 画像シフト部
２６０ 関心領域設定部（図５では、ＲＯＩ設定部と表示）
２８０ 画像削除部
３００ 非関心領域削除部
３２０ 代表スライス作成部
３４０ パラメータ算出部
３６０ 入力部
３８０ 判断部
４００ 再構成部
４２０ 出力部
【特許請求の範囲】
【請求項１】
コンピュータが心筋核医学画像の投影データセットを取得する投影データセット取得ステップと、
コンピュータが、前記投影データセットに含まれる各投影画像上で心筋部位を含む領域に関心領域を設定する関心領域設定ステップと、
コンピュータが、心筋部以外の高集積部位が前記関心領域と重なっている投影画像を削除する画像削除ステップと、
コンピュータが、前記画像削除ステップにて心筋部以外の高集積部位が前記関心領域内と重なっている投影画像を削除した後の投影データセットに含まれる各投影画像につき、前記関心領域設定ステップにて設定した関心領域以外の部分を削除する非関心領域削除ステップと、
コンピュータが、非関心領域削除ステップ後の投影データセットを用いて画像再構成を行う画像再構成ステップと、を含む、画像処理方法。
【請求項２】
関心領域設定ステップの前に、コンピュータが、投影データセットに含まれる各投影画像において、心筋部位の中心を画面の回転中心に移動させる、画像シフトステップをさらに含み、
関心領域設定ステップは、コンピュータが、画像シフトステップ処理後の投影データセットにおける各投影画像上で心筋部位を含む領域に関心領域を設定するものである、請求項１に記載の画像処理方法。
【請求項３】
非関心領域削除ステップにおける関心領域以外の部分の削除が、マスク処理にて行われるものである、請求項１または請求項２に記載の、画像処理方法。
【請求項４】
非関心領域削除ステップにおける関心領域以外の部分の削除が、関心領域外の部分をバックグラウンドの平均値で補完する事により行われるものである、請求項１〜請求項３のいずれか１項に記載の、画像処理方法。
【請求項５】
非関心領域削除ステップ後の各投影画像において、非関心領域の削除処理後の画像上で、削除された領域を関心領域境界における値で補完するデータ補完ステップをさらに含み、
画像再構成ステップは、データ補完ステップ後の投影データセットを用いて画像再構成を行うものである、請求項１〜３のいずれか１項に記載の、画像処理方法。
【請求項６】
コンピュータに、心筋核医学画像の投影データセットを取得する投影データ取得ステップと、
前記投影データセットに含まれる各投影画像上で心筋部位を含む領域に関心領域を設定する関心領域設定ステップと、
心筋部以外の高集積部位が前記関心領域と重なっている投影画像を削除する画像削除ステップと、
前記画像削除ステップにて心筋部以外の高集積部位が前記関心領域内と重なっている投影画像を削除した後の投影データセットに含まれる各投影画像につき、前記関心領域設定ステップにて設定した関心領域以外の部分を削除する非関心領域削除ステップと、
非関心領域削除ステップ後の投影データセットを用いて画像再構成を行う画像再構成ステップと、を実行させる、画像処理プログラム。
【請求項７】
関心領域設定ステップの前に、コンピュータに、投影データセットに含まれる各投影画像において、心筋部位の中心を画面の回転中心に移動させる、画像シフトステップをさらに実行させ、
関心領域設定ステップにおいては、コンピュータに、画像シフトステップ処理後の投影データセットにおける各投影画像上で心筋部位を含む領域に関心領域を設定する処理を実行させるものである、請求項６に記載の画像処理プログラム。
【請求項８】
非関心領域削除ステップにおける関心領域以外の部分の削除が、コンピュータにマスク処理を実行させる事により行われるものである、請求項６または請求項７に記載の、画像処理プログラム。
【請求項９】
非関心領域削除ステップにおける関心領域以外の部分の削除が、コンピュータに関心領域外の部分をバックグラウンドの平均値で補完する処理を実行させる事により行われるものである、請求項６〜請求項８のいずれか１項に記載の、画像処理プログラム。
【請求項１０】
コンピュータに、非関心領域の削除処理後の画像上で、削除された領域を関心領域境界における値で補完する、データ補完ステップをさらに実行させ、
画像再構成ステップは、コンピュータに、データ補完ステップ後の投影データセットを用いて画像再構成を実行させるものである、請求項６〜８のいずれか１項に記載の、画像処理プログラム。
【請求項１１】
心筋核医学画像の投影データセットを取り込む投影データセット取得部と、
前記投影データセットに含まれる各投影画像上で心筋部位を含む領域に関心領域を設定する関心領域設定部と、
前記投影データセットに含まれる各投影画像のうち、心筋部位以外の高集積部位が前記関心領域と重なっている投影画像を削除する画像削除部と、
前記関心領域設定部の機能により設定された前記関心領域以外の部分を各投影画像上において削除する非関心領域削除部と、
投影データセットを用いて画像再構成を行う画像再構成部と、を備える、画像処理装置。
【請求項１２】
投影データセットに含まれる各投影画像について、心筋部位の中心を画面の回転中心に移動させる機能を有する、画像シフト部をさらに備える、請求項１１に記載の、画像処理装置。
【請求項１３】
非関心領域の削除処理後の画像上で、削除された領域を関心領域境界における値で補完する機能を有する、データ補完部をさらに備える、請求項１１または請求項１２に記載の、画像処理装置。
【請求項１】
コンピュータが心筋核医学画像の投影データセットを取得する投影データセット取得ステップと、
コンピュータが、前記投影データセットに含まれる各投影画像上で心筋部位を含む領域に関心領域を設定する関心領域設定ステップと、
コンピュータが、心筋部以外の高集積部位が前記関心領域と重なっている投影画像を削除する画像削除ステップと、
コンピュータが、前記画像削除ステップにて心筋部以外の高集積部位が前記関心領域内と重なっている投影画像を削除した後の投影データセットに含まれる各投影画像につき、前記関心領域設定ステップにて設定した関心領域以外の部分を削除する非関心領域削除ステップと、
コンピュータが、非関心領域削除ステップ後の投影データセットを用いて画像再構成を行う画像再構成ステップと、を含む、画像処理方法。
【請求項２】
関心領域設定ステップの前に、コンピュータが、投影データセットに含まれる各投影画像において、心筋部位の中心を画面の回転中心に移動させる、画像シフトステップをさらに含み、
関心領域設定ステップは、コンピュータが、画像シフトステップ処理後の投影データセットにおける各投影画像上で心筋部位を含む領域に関心領域を設定するものである、請求項１に記載の画像処理方法。
【請求項３】
非関心領域削除ステップにおける関心領域以外の部分の削除が、マスク処理にて行われるものである、請求項１または請求項２に記載の、画像処理方法。
【請求項４】
非関心領域削除ステップにおける関心領域以外の部分の削除が、関心領域外の部分をバックグラウンドの平均値で補完する事により行われるものである、請求項１〜請求項３のいずれか１項に記載の、画像処理方法。
【請求項５】
非関心領域削除ステップ後の各投影画像において、非関心領域の削除処理後の画像上で、削除された領域を関心領域境界における値で補完するデータ補完ステップをさらに含み、
画像再構成ステップは、データ補完ステップ後の投影データセットを用いて画像再構成を行うものである、請求項１〜３のいずれか１項に記載の、画像処理方法。
【請求項６】
コンピュータに、心筋核医学画像の投影データセットを取得する投影データ取得ステップと、
前記投影データセットに含まれる各投影画像上で心筋部位を含む領域に関心領域を設定する関心領域設定ステップと、
心筋部以外の高集積部位が前記関心領域と重なっている投影画像を削除する画像削除ステップと、
前記画像削除ステップにて心筋部以外の高集積部位が前記関心領域内と重なっている投影画像を削除した後の投影データセットに含まれる各投影画像につき、前記関心領域設定ステップにて設定した関心領域以外の部分を削除する非関心領域削除ステップと、
非関心領域削除ステップ後の投影データセットを用いて画像再構成を行う画像再構成ステップと、を実行させる、画像処理プログラム。
【請求項７】
関心領域設定ステップの前に、コンピュータに、投影データセットに含まれる各投影画像において、心筋部位の中心を画面の回転中心に移動させる、画像シフトステップをさらに実行させ、
関心領域設定ステップにおいては、コンピュータに、画像シフトステップ処理後の投影データセットにおける各投影画像上で心筋部位を含む領域に関心領域を設定する処理を実行させるものである、請求項６に記載の画像処理プログラム。
【請求項８】
非関心領域削除ステップにおける関心領域以外の部分の削除が、コンピュータにマスク処理を実行させる事により行われるものである、請求項６または請求項７に記載の、画像処理プログラム。
【請求項９】
非関心領域削除ステップにおける関心領域以外の部分の削除が、コンピュータに関心領域外の部分をバックグラウンドの平均値で補完する処理を実行させる事により行われるものである、請求項６〜請求項８のいずれか１項に記載の、画像処理プログラム。
【請求項１０】
コンピュータに、非関心領域の削除処理後の画像上で、削除された領域を関心領域境界における値で補完する、データ補完ステップをさらに実行させ、
画像再構成ステップは、コンピュータに、データ補完ステップ後の投影データセットを用いて画像再構成を実行させるものである、請求項６〜８のいずれか１項に記載の、画像処理プログラム。
【請求項１１】
心筋核医学画像の投影データセットを取り込む投影データセット取得部と、
前記投影データセットに含まれる各投影画像上で心筋部位を含む領域に関心領域を設定する関心領域設定部と、
前記投影データセットに含まれる各投影画像のうち、心筋部位以外の高集積部位が前記関心領域と重なっている投影画像を削除する画像削除部と、
前記関心領域設定部の機能により設定された前記関心領域以外の部分を各投影画像上において削除する非関心領域削除部と、
投影データセットを用いて画像再構成を行う画像再構成部と、を備える、画像処理装置。
【請求項１２】
投影データセットに含まれる各投影画像について、心筋部位の中心を画面の回転中心に移動させる機能を有する、画像シフト部をさらに備える、請求項１１に記載の、画像処理装置。
【請求項１３】
非関心領域の削除処理後の画像上で、削除された領域を関心領域境界における値で補完する機能を有する、データ補完部をさらに備える、請求項１１または請求項１２に記載の、画像処理装置。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【公開番号】特開２０１１−２２１１９（Ｐ２０１１−２２１１９Ａ）
【公開日】平成２３年２月３日（２０１１．２．３）
【国際特許分類】
【出願番号】特願２００９−１８５２０２（Ｐ２００９−１８５２０２）
【出願日】平成２１年７月１７日（２００９．７．１７）
【特許番号】特許第４４９４５１０号（Ｐ４４９４５１０）
【特許公報発行日】平成２２年６月３０日（２０１０．６．３０）
【出願人】（０００１４９８３７）富士フイルムＲＩファーマ株式会社 (54)
【Ｆターム（参考）】
【公開日】平成２３年２月３日（２０１１．２．３）
【国際特許分類】
【出願日】平成２１年７月１７日（２００９．７．１７）
【特許番号】特許第４４９４５１０号（Ｐ４４９４５１０）
【特許公報発行日】平成２２年６月３０日（２０１０．６．３０）
【出願人】（０００１４９８３７）富士フイルムＲＩファーマ株式会社 (54)
【Ｆターム（参考）】
[ Back to top ]