(57)【要約】【課題】  病変部のコントラストを向上可能な画像補正方法を提供する。【解決手段】  本発明の画像補正方法は、（ａ）ＣＴ画像の各画素のＣＴ値を対応の線吸収係数に置換してなる吸収補正マップであって、骨領域の線吸収係数を水又は生体軟部組織の線吸収係数に置換してなる吸収補正マップを生成するステップと、（ｂ）上記のＣＴ画像が取得された断層面と実質的に同一の断層面に対して取得された核医学画像の各画素と、上記のステップにて生成された前記吸収補正マップの対応の画素とを乗算してなる補正画像を生成するステップと、を含む。【選択図】  図３
2006047152.tif 000002
【特許請求の範囲】
【請求項１】
ＣＴ画像の各画素のＣＴ値を対応の線吸収係数に置換してなる吸収補正マップであって、骨領域の線吸収係数を水又は生体軟部組織の線吸収係数に置換してなる吸収補正マップを生成するステップと、
前記ＣＴ画像が取得された断層面と実質的に同一の断層面に対して取得された核医学画像の各画素と、前記ステップにて生成された前記吸収補正マップの対応の画素とを乗算してなる補正画像を生成するステップと、
を含む、画像補正方法。
【請求項２】
前記吸収補正マップを生成するステップにおいて、
予め定められたＣＴ値と線吸収係数との対応に基づき、前記ＣＴ画像の各画素におけるＣＴ値を対応の線吸収係数に置換することによって、吸収補正マップを生成し、
該吸収補正マップにおける線吸収係数のうち、骨領域の線吸収係数として予め定められた範囲に含まれる線吸収係数を、水又は生体軟部組織の線吸収係数として予め定められた線吸収係数に置換する、
請求項１に記載の画像補正方法。
【請求項３】
前記吸収補正マップを生成するステップにおいて、
前記ＣＴ画像におけるＣＴ値のうち、骨領域におけるＣＴ値として予め定められた範囲に含まれているＣＴ値を、水又は生体軟部組織のＣＴ値として予め定められたＣＴ値に置換することによって、第２のＣＴ画像を生成し、
予め定められたＣＴ値と線吸収係数との対応に基づき、前記第２のＣＴ画像における各画素のＣＴ値を対応の線吸収係数に置換することによって、前記吸収補正マップを生成する、
請求項１に記載の画像補正方法。
【請求項４】
前記吸収補正マップを生成するステップにおいて、前記吸収補正マップの解像度を前記核医学画像の解像度に変換するためのフィルタ処理を、該吸収補正マップに適用し、
前記補正画像を生成するステップにおいて、前記核医学画像の各画素と、前記フィルタ処理が適用された前記吸収補正マップの対応の画素とを乗算することによって、前記補正画像を生成する、
請求項１〜３の何れか一項に記載の画像補正方法。
【請求項５】
ＣＴ画像の各画素のＣＴ値を対応の線吸収係数に置換してなる吸収補正マップであって、骨領域の線吸収係数を水又は生体軟部組織の線吸収係数に置換してなる吸収補正マップを生成する手段と、
前記ＣＴ画像が取得された断層面と実質的に同一の断層面に対して取得された核医学画像の各画素と、前記手段にて生成された前記吸収補正マップの対応の画素とを乗算してなる補正画像を生成する手段と、
を備える、画像補正装置。
【請求項６】
前記吸収補正マップを生成する手段は、
予め定められたＣＴ値と線吸収係数との対応に基づき、前記ＣＴ画像の各画素におけるＣＴ値を対応の線吸収係数に置換することによって、吸収補正マップを生成し、
該吸収補正マップにおける線吸収係数のうち、骨領域の線吸収係数として予め定められた範囲に含まれる線吸収係数を、水又は生体軟部組織の線吸収係数として予め定められた線吸収係数に置換する、
請求項５に記載の画像補正装置。
【請求項７】
前記吸収補正マップを生成する手段は、
前記ＣＴ画像におけるＣＴ値のうち、骨領域におけるＣＴ値として予め定められた範囲に含まれているＣＴ値を、水又は生体軟部組織のＣＴ値として予め定められたＣＴ値に置換することによって、第２のＣＴ画像を生成し、
予め定められたＣＴ値と線吸収係数との対応に基づき、前記第２のＣＴ画像における各画素のＣＴ値を対応の線吸収係数に置換することによって、前記吸収補正マップを生成する、
請求項５に記載の画像補正装置。
【請求項８】
前記吸収補正マップを生成する手段は、前記吸収補正マップの解像度を前記核医学画像の解像度に変換するためのフィルタ処理を、該吸収補正マップに適用し、
前記補正画像を生成する手段は、前記核医学画像の各画素と、前記フィルタが適用された前記吸収補正マップの対応の画素とを乗算することによって、前記補正画像を生成する、
請求項５〜７の何れか一項に記載の画像補正装置。
【請求項９】
コンピュータに、
ＣＴ画像の各画素のＣＴ値を対応の線吸収係数に置換してなる吸収補正マップであって、骨領域の線吸収係数を水又は生体軟部組織の線吸収係数に置換してなる吸収補正マップを生成するステップと、
前記ＣＴ画像が取得された断層面と実質的に同一の断層面に対して取得された核医学画像の各画素と、前記ステップにて生成された前記吸収補正マップの対応の画素とを乗算してなる補正画像を生成するステップと、
を実行させるための画像補正プログラム。
【請求項１０】
前記吸収補正マップを生成するステップにおいて、コンピュータに、
予め定められたＣＴ値と線吸収係数との対応に基づき、前記ＣＴ画像の各画素におけるＣＴ値を対応の線吸収係数に置換することによって、吸収補正マップを生成させ、
該吸収補正マップにおける線吸収係数のうち、骨領域の線吸収係数として予め定められた範囲に含まれる線吸収係数を、水又は生体軟部組織の線吸収係数として予め定められた線吸収係数に置換させる、
請求項９に記載の画像補正プログラム。
【請求項１１】
前記吸収補正マップを生成するステップにおいて、コンピュータに、
前記ＣＴ画像におけるＣＴ値のうち、骨領域におけるＣＴ値として予め定められた範囲に含まれているＣＴ値を、水又は生体軟部組織のＣＴ値として予め定められたＣＴ値に置換することによって、第２のＣＴ画像を生成させ、
予め定められたＣＴ値と線吸収係数との対応に基づき、前記第２のＣＴ画像における各画素のＣＴ値を対応の線吸収係数に置換することによって、前記吸収補正マップを生成させる、
請求項９に記載の画像補正プログラム。
【請求項１２】
前記吸収補正マップを生成するステップにおいて、更にコンピュータに、前記吸収補正マップの解像度を前記核医学画像の解像度に変換するためのフィルタ処理を該吸収補正マップに適用させ、
前記補正画像を生成するステップにおいて、前記核医学画像の各画素と、前記フィルタ処理が適用された前記吸収補正マップの対応の画素とを乗算することによって、前記補正画像をコンピュータに生成させる、
請求項９〜１１の何れか一項に記載の画像補正プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像補正方法、画像補正装置、及び画像補正プログラムに関する。
【背景技術】
【０００２】
陽電子放出断層撮影（以下、「ＰＥＴ」という）、及び単一光子放射断層撮影（以下、「ＳＰＥＣＴ」という）といった核医学検査によって得られる核医学画像は、癌を始めとする種々の疾患を診断するために、広く用いられている。
【０００３】
核医学画像は、被検者に放射性同位元素で標識した薬剤（以下、「放射性医薬品」という）を投与し、その医薬品の投与により直接又は二次的に放出されるγ線を、γカメラ等によって検出することによって得られる。すなわち、核医学画像は、投与した放射性医薬品の分布を画像化したものであるので、特定の機能を有した疾患部の検出に適している。例えば、放射性医薬品の一つである２−［１８F］フルオロ−２−デオキシ−D−グルコースは、糖代謝の盛んな腫瘍に集積する性質があるため、糖代謝が盛んな腫瘍を特異的に検出することが可能となる。
【０００４】
このように、核医学検査は、放射性医薬品の種類によって種々の疾患の診断に適用可能な核医学画像を得ることができるという利点を有するが、核医学画像の画質は、ＣＴ画像といった他の方法によって得られる生体断層像に比して低いという問題がある。例えば、放射源からの距離に依存するγ線の吸収を要因とする画質劣化、或いは、病変部のコントラストの低さといった問題がある。
【０００５】
前者、すなわち、γ線の吸収を要因とする画質劣化に対しては、ＣＴ画像より作成した吸収補正マップを用い、線吸収係数および体表面からの深さを考慮して核医学画像を補正する方法（Ordered Subset Expectation Maximization法（以下、「OSEM法」という）が提案されている（例えば、非特許文献１、非特許文献２を参照）。
【非特許文献１】Patton, J.A., Delbeke, D. and Sandler, M.P., “Image Fusion Using an Integrated, Dual-Head Coincidence Camera with X-Ray Tube-Based Attenuation Maps.”, The Journal of Nuclear Medicine, (USA), Society of Nuclear Medicine, 2000, 41, 8, p.1364-1368.
【非特許文献２】Hudson, H. M. and Larkin, R. S., “Accelerated Image Reconstruction Using Ordered Subsets of Projection Data.”, IEEE Transactions on Medical Imaging, (USA), IEEE Engineering in Medicine and Biology Society et al., 1994, 13, 4, p.601-609.
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかしながら、核医学画像における病変部のコントラストを向上するという課題は、未だ解決されていない。
【０００７】
本発明は上記問題点に鑑みてなされたもので、病変部のコントラストを向上可能な画像補正方法、画像補正装置、及び画像補正プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００８】
本発明の一側面に係る画像補正方法は、（ａ）ＣＴ画像の各画素のＣＴ値を対応の線吸収係数に置換してなる吸収補正マップであって、骨領域の線吸収係数を水又は生体軟部組織の線吸収係数に置換してなる吸収補正マップを生成するステップと、（ｂ）上記のＣＴ画像が取得された断層面と実質的に同一の断層面に対して取得された核医学画像の各画素と、上記のステップにて生成された前記吸収補正マップの対応の画素とを乗算してなる補正画像を生成するステップと、を含む。
【０００９】
本発明の別の一側面に係る画像補正装置は、（ａ）ＣＴ画像の各画素のＣＴ値を対応の線吸収係数に置換してなる吸収補正マップであって、骨領域の線吸収係数を水又は生体軟部組織の線吸収係数に置換してなる吸収補正マップを生成する手段と、（ｂ）上記のＣＴ画像が取得された断層面と実質的に同一の断層面に対して取得された核医学画像の各画素と、上記の手段にて生成された吸収補正マップの対応の画素とを乗算してなる補正画像を生成する手段と、
を備える。
【００１０】
本発明の更に別の一側面に係る画像補正処理プログラムは、コンピュータに、（ａ）ＣＴ画像の各画素のＣＴ値を対応の線吸収係数に置換してなる吸収補正マップであって、骨領域の線吸収係数を水又は生体軟部組織の線吸収係数に置換してなる吸収補正マップを生成するステップと、（ｂ）上記のＣＴ画像が取得された断層面と実質的に同一の断層面に対して取得された核医学画像の各画素と、上記のステップにて生成された吸収補正マップの対応の画素とを乗算してなる補正画像を生成するステップと、を実行させる。
【００１１】
本発明によれば、相対的に他の領域よりその値が大きい骨領域の線吸収係数が、水又は生体軟部組織の線吸収係数に置換された吸収補正マップが用いられる。この吸収補正マップが核医学画像に乗算されるので、その結果得られる補正画像では、バックグランドノイズが低減されて、病変部のコントラストが高くなる。また、吸収補正マップを生成するための置換処理、及び補正画像を生成するための乗算処理によって、補正画像が生成される。したがって、本発明によれば、術者の経験によらず短時間で補正画像を生成することが可能である。
【００１２】
上述の本発明において、予め定められたＣＴ値と線吸収係数との対応に基づき、（ｃ１）ＣＴ画像の各画素におけるＣＴ値を対応の線吸収係数に置換することによって、吸収補正マップを生成し、（ｄ１）骨領域の線吸収係数として予め定められた範囲に含まれる線吸収係数を、水又は生体軟部組織の線吸収係数として予め定められた線吸収係数に置換してもよい。
【００１３】
また、上述の本発明において、（ｃ２）ＣＴ画像におけるＣＴ値のうち、骨領域におけるＣＴ値として予め定められた範囲に含まれているＣＴ値を、水又は生体軟部組織のＣＴ値として予め定められたＣＴ値に置換することによって、第２のＣＴ画像を生成し、（ｄ２）予め定められたＣＴ値と線吸収係数との対応に基づき、第２のＣＴ画像における各画素のＣＴ値を対応の線吸収係数に置換することによって、吸収補正マップを生成してもよい。
【００１４】
更に、本発明では、（ｅ）吸収補正マップの解像度を前記核医学画像の解像度に変換するためのフィルタ処理を、当該吸収補正マップに適用し、（ｆ）核医学画像の各画素と、フィルタ処理が適用された吸収補正マップの対応の画素とを乗算することによって、補正画像を生成してもよい。
【発明の効果】
【００１５】
本発明によれば、ＰＥＴ及びＳＰＥＣＴによって得られる核医学画像における病変部のコントラストを向上させた補正画像を得ることが可能であり、かつ、当該補正画像を術者の経験によらずに短時間に得ることが可能な画像補正方法、画像補正装置、及び画像補正プログラムが提供される。
【発明を実施するための最良の形態】
【００１６】
以下、図面を参照して本発明の好適な実施形態について詳細に説明する。なお、各図面において同一又は相当の部分に対しては同一の符号を附すこととする。
【００１７】
図１は、本発明の実施の形態にかかる画像補正装置の構成を示す図である。図１に示す画像補正装置１０は、物理的には、ＣＰＵ（中央処理装置）、ハードディスク、メモリ、ディスプレイといった表示デバイス、キーボード・マウスといった入力デバイス等を備えるコンピュータである。
【００１８】
画像補正装置１０は、機能的に、入力部１２、吸収補正マップ生成部１４、補正画像生成部１６、及び出力部１８を備えている。以下、これらの構成要素の機能を、図１を参照しつつ説明し、併せて、本発明の実施の形態に係る画像補正方法について説明する。図２は、本発明の実施の形態に係る画像補正方法のフローチャートである。
【００１９】
入力部１２は、実質的に同一の生体断層面に対して取得された核医学画像及びＣＴ画像（Ｘ線ＣＴ画像）を入力する（ステップＳ０１）。入力部１２は、例えば、公知のＳＰＥＣＴ−ＣＴ又はＰＥＴ−ＣＴ装置によって取得された核医学画像、及びＣＴ画像を入力する。
【００２０】
ＳＰＥＣＴ−ＣＴ装置、ＰＥＴ−ＣＴ装置は、Ｘ線ＣＴ装置とＳＰＥＣＴ装置、Ｘ線ＣＴ装置とＰＥＴ装置が、各々一体化された装置である。ＳＰＥＣＴ−ＣＴ装置、ＰＥＴ−ＣＴ装置は、生体の同一断層面におけるＸ線ＣＴ画像とＳＰＥＣＴ画像、Ｘ線ＣＴ画像とＰＥＴ画像を、自動的に各々取得することが可能である。ＳＰＥＣＴ−ＣＴ又はＰＥＴ−ＣＴ装置の詳細については、Keidar, Z., Israel, O. and Krausz, Y., “SPECT/CT in Tumor Imaging: Technical Aspect and Clinical Applications.”, Seminars in Nuclear Medicine, (USA), Elsevier, 2003, 33,3, p.205-218を参照されたい。以下、上述の核医学画像としてＳＰＥＣＴ画像を用いる場合について説明する。
【００２１】
吸収補正マップ生成部１４は、吸収補正マップを生成する（ステップＳ０２）。吸収補正マップは、Ｘ線ＣＴ画像におけるＣＴ値を対応の線吸収係数に置換したマップであって、骨領域における線吸収係数が水又は生体軟部組織の線吸収係数に置換されたマップである。また、吸収補正マップは、ＳＰＥＣＴ画像の解像度に対応する解像度に変換されたものである。吸収補正マップの生成方法についての詳細は後述する。
【００２２】
補正画像生成部１６は、核医学画像の各画素と吸収補正マップの対応の画素との乗算を実行することによって、補正画像を生成する（ステップＳ０３）。
【００２３】
出力部１８は、上記の補正画像をディスプレイといった表示装置に出力する（ステップＳ０４）。出力部１８による表示装置への表示は、補正画像の信号強度に応じた輝度又は色彩等によるものある。なお、出力部１８による表示装置への表示は、補正画像の二次元表示のみならず、三次元表示であってもよい。三次元表示の場合には、複数の生体断面各々に対して生成された複数の補正画像をフォーマット処理（座標変換処理）し、矢状断像及び冠状断像を生成する。これらの画像を二次元空間上に表示することによって、三次元表示が可能となる。
【００２４】
以下、吸収補正マップ生成部１４による吸収補正マップの生成方法について詳細に説明する。図３は、吸収補正マップを生成する処理のフローチャートである。
【００２５】
図３に示す実施の形態の吸収補正マップの生成処理では、吸収補正マップ生成部１４は、まず、ＣＴ画像のＣＴ値を対応の線吸収係数に置換することによって、吸収補正マップを生成する（ステップＳ１１）。
【００２６】
ここで、ＣＴ画像では、骨のＣＴ値を１０００Ｈｕ（ハンスフィールド値）、水のＣＴ値を０Ｈｕ、空気のＣＴ値を−１０００Ｈｕとした相対表示によって画像が構成されている。一方、例えば、エネルギー５１１ｋｅＶでは、水（ＣＴ値＝０Ｈｕ）の線吸収係数は、０．０９６／ｃｍ、骨（ＣＴ＝１０００Ｈｕ）の線吸収係数は、０．１７２／ｃｍ、空気（ＣＴ値＝−１０００Ｈｕ）の線吸収係数は０／ｃｍである。
【００２７】
吸収補正マップ生成部１４は、これらの値を用い、空気領域（−１０００≦ＣＴ値≦０）に対しては次式（１）を用い、水又は固体領域に対しては次式（２）を用いて、ＣＴ値を線吸収係数に置換する。なお、次式（１）及び（２）において、ＸはＣＴ値、Ｙは線吸収係数である。
Ｙ＝０．００００９６Ｘ＋０．０９６（ＣＴ値＜０）        …（１）
Ｙ＝０．００００７６Ｘ＋０．０９６（ＣＴ値＞０）        …（２）
【００２８】
次に、吸収補正マップ生成部１４は、ステップＳ１１にて生成された吸収補正マップにおける線吸収係数のうち、骨領域の線吸収係数を水又は生体軟部組織の線吸収係数に置換する（ステップＳ１２）。例えば、エネルギー５１１ｋｅＶでは、骨の線吸収係数は０．１７２／ｃｍである。したがって、吸収補正マップ生成部１４は、０．１７２／ｃｍを閾値とし、ステップＳ１１にて生成された吸収補正マップにおける線吸収係数のうち当該閾値以上の線吸収係数を、水の線吸収係数である０．０９６／ｃｍに置換する。これによって、吸収補正マップが生成される。
【００２９】
次に、吸収補正マップ生成部１４は、吸収補正マップの解像度をＳＰＥＣＴ画像の解像度に一致させるためのフィルタ処理を実行する（ステップＳ１３）。このフィルタ処理において、吸収補正マップ生成部１４は、ＳＰＥＣＴ画像の再構成に用いられたフィルタと同様のフィルタを用いる。当該フィルタには、例えば、バターワースフィルタが用いられる。
【００３０】
より詳細には、吸収補正マップ生成部１４は、吸収補正マップを投影したプロファイルｐ（ｓ，θ）を複数生成する。すなわち、吸収補正マップ生成部１４は、θを所定角度のピッチで変更することによって、吸収補正マップのプロファイルｐ（ｓ，θ）を複数生成する。
【００３１】
次いで、吸収補正マップ生成部１４は、複数のプロファイルｐ（ｓ，θ）を各々フーリエ変換することによって、プロファイルの周波数空間表現Ｐθ（ｆ）を複数生成する。吸収補正マップ生成部１４は、Ｐθ（ｆ）の各々に、次式（３）に示すバターワースフィルタＢ（ｆ）を、次式（４）に示すようにＰθ（ｆ）に適用することによって、Ｇθ（ｆ）とする。
Ｂ（ｆ）＝１／（１＋（ｆ／ｆｃ））２ｎ      …（３）
Ｇθ（ｆ）＝Ｐθ（ｆ）・Ｂ（ｆ）              …（４）
【００３２】
なお、式（３）において、ｆｃは遮断周波数、ｎはオーダである。ｎの値には、任意に選択された値を用いることができるが、ｎが小さすぎると高周波成分が多くなり、統計雑音が増加する。一方、ｎの値が大きすぎるとリング状のアーチファクトが発生するので、好ましくない。ｎの値には、好ましくは２〜８が用いられ、特に好ましくは５が用いられる。
【００３３】
ｆｃの値は、用いた放射性医薬品の種類や測定対象部位等に応じて任意に選択されるが、０．５ｃｙｃｌｅ／ｐｉｘｅｌ以下である必要がある。ここで、ｆｃの値が大きすぎると統計雑音が増加する。一方、ｆｃの値が小さすぎると必要な情報までもが除去され、画像の分解能が劣化する。したがって、例えば、Ｇａにて腫瘍のイメージング（ピクセルサイズ：０．７７ｃｍ）を行う場合には、０．２〜０．３ｃｙｃｌｅ／ｃｍのｆｃが好適に用いられ、特に好適には０．２５ｃｙｃｌｅ／ｃｍのｆｃが用いられる。また、ＦＤＧによる腫瘍イメージングの場合には、０．３５〜０．４５ｃｙｃｌｅ／ｃｍのｆｃが好適に用いられ、特に好適には０．４ｃｙｃｌｅ／ｃｍのｆｃが用いられる。
【００３４】
次いで、吸収補正マップ生成部１４は、上記のＧθ（ｆ）を各々、角度θ方向に逆投影することによって吸収補正マップを得て、フィルタ処理（ステップＳ１３）を完了する。
【００３５】
次に、吸収補正マップ生成部１４は、吸収補正マップにスケーリング処理を施す（ステップＳ１４）。スケーリング処理は、補正画像の信号強度が上限値を超えないようにする処理である。具体的に、吸収補正マップ生成部１４は、次式（５）の関係式を満足する数値Ｎを求め、吸収補正マップの各線吸収係数を当該Ｎで除す。なお、次式（５）において、ａは吸収補正マップの各画素の信号強度、ｂは補正対象画像における対応画素の信号強度、Ｍｍａｘは上記の上限値である。以上の処理によって、補正画像の生成に用いられる吸収補正マップが生成される。
ａ×ｂ／Ｎ＜Ｍｍａｘ      …（５）
【００３６】
以下、別の実施の形態の吸収補正マップの生成方法について、説明する。図４は、吸収補正マップを生成する処理の他の実施の形態を示すフローチャートである。
【００３７】
図４に示す別の実施の形態の吸収補正マップの生成処理では、吸収補正マップ生成部１４は、まず、ＣＴ画像のＣＴ値のうち骨領域のＣＴを、水又は生体軟部組織のＣＴ値に置換する（ステップＳ３１）。具体的に、吸収補正マップ生成部１４は、１０００Ｈｕ以上のＣＴ値を０Ｈｕに置換する。
【００３８】
次いで、吸収補正マップ生成部１４は、ステップＳ３１にて置換されてなるＣＴ画像におけるＣＴ値を対応の線吸収係数に置換する（ステップＳ３２）。この置換処理においても、上述した式（１）及び（２）を用いることができる。
【００３９】
次いで、吸収補正マップ生成部１４は、上述したフィルタ処理と同様の処理によって、ステップＳ３２にて生成された吸収補正マップの解像度を、ＳＰＥＣＴ画像の解像度に一致させる（ステップＳ３３）。
【００４０】
次いで、吸収補正マップ生成部１４は、上述の式（５）の関係を満足するように、ステップＳ３３にて生成された吸収補正マップにスケーリング処理を施すことによって、補正画像の生成に用いる吸収補正マップを生成する（ステップＳ３４）。
【００４１】
以上のように、吸収補正マップ生成部１４は、種々の方法によって、吸収補正マップを生成することができる。
【００４２】
以下、コンピュータを画像補正装置１０として動作させる画像補正プログラムについて説明する。図５は、本発明の実施の形態に係る画像補正プログラムの構成を、記録媒体と共に示す図である。
【００４３】
図５に示す画像補正プログラム４０は、記録媒体１００に格納されて提供される。記録媒体１００としては、フロッピーディスク、ＣＤ−ＲＯＭ、ＤＶＤ、あるいはＲＯＭ等の記録媒体あるいは半導体メモリ等が例示される。
【００４４】
図６は、記録媒体に記憶されたプログラムを実行するためのコンピュータのハードウェア構成を示す図であり、図７は、記録媒体に記憶されたプログラムを実行するためのコンピュータの斜視図である。図６に示すように、コンピュータ１１０は、フロッピーディスクドライブ装置、ＣＤ−ＲＯＭドライブ装置、ＤＶＤドライブ装置等の読取装置１１２と、オペレーティングシステムを常駐させた作業用メモリ（ＲＡＭ）１１４と、記録媒体１００に記憶されたプログラムを記憶するメモリ１１６と、ディスプレイといった表示装置１１８と、入力装置であるマウス１２０及びキーボード１２２と、データ等の送受を行うための通信装置１２４と、プログラムの実行を制御するＣＰＵ１２６とを備えている。コンピュータ１１０は、記録媒体１００が読取装置１１２に挿入されると、読取装置１１２から記録媒体１００に格納された画像補正プログラム４０にアクセス可能になり、当該画像補正プログラム４０によって、画像補正装置１０として動作することが可能になる。
【００４５】
図７に示すように、画像補正プログラム４０は、搬送波に重畳されたコンピュータデータ信号１３０としてネットワークを介して提供されるものであってもよい。この場合、コンピュータ１１０は、通信装置１２４によって受信した画像補正プログラム４０をメモリ１１６に格納し、当該画像補正プログラム４０を実行することができる。
【００４６】
図５に示すように、画像補正プログラム４０は、処理を統括するメインモジュール４１と、入力モジュール４２と、吸収補正マップ生成モジュール４４と、補正画像生成モジュール４６と、出力モジュール４８とを備えている。
【００４７】
入力モジュール４２は上記のステップＳ０１、吸収補正マップ生成モジュール４４は上記のステップＳ０２、補正画像生成モジュール４６上記のステップＳ０３、出力モジュール４８は上記のステップＳ０４を、コンピュータに実行させる。
【実施例１】
【００４８】
以下、本発明の好適な実施例についてさらに詳細に説明するが、本発明はこれらの実施例に限定されるものではない。
【００４９】
（実施例１）
悪性リンパ腫の患者（75歳、男性）に対し、クエン酸ガリウム（67Ga）注ＮＭＰ（商品名、日本メジフィジックス株式会社製）を１４８ＭBq静脈内投与し、投与７２時間後にＳＰＥＣＴ−ＣＴ装置（形式Millennium VG＋Hawkeye、ＧＥ社製）を用いて胸部のＳＰＥＣＴ画像（横断像）（撮像条件、64×64マトリックス、60方向／360度、30秒／方向）を得た。
【００５０】
次いで、上記ＳＰＥＣＴ画像と同一の断層面におけるＣＴ画像の撮像を行い（撮像条件140kV、2.5mA、15秒／スライス、128x128マトリックス）、当該ＣＴ画像を用いて吸収補正マップを作成した。当該吸収補正マップにおける線吸収係数のヒストグラム上で、骨の線吸収係数の値（0.214/cm）に閾値を設定し、骨領域と骨以外の領域とに区分わけを行った。この吸収補正マップにおける骨領域の線吸収係数を、水における線吸収係数（0.119/cm）で置換した。
【００５１】
置換後の吸収補正マップにフィルタ処理を施し、吸収補正マップの空間分解能を、ＳＰＥＣＴ画像と同一の空間分解能に変換した。
【００５２】
上記フィルタ処理後の吸収補正マップとＳＰＥＣＴ画像との間で、同一の位置における画素同士の信号強度を乗算し、得られた結果を画面上に表し、補正画像を得た。
【００５３】
（比較例１）
実施例１にて取得したＣＴ画像及びＳＰＥＣＴ画像を用いて、ＯＳＥＭ法によるＳＰＥＣＴ画像の補正を行った。ＯＳＥＭ法による画像処理は、撮像に用いたＳＰＥＣＴ−ＣＴ装置（形式Millennium VG＋Hawkeye、ＧＥ社製）に内蔵されているプログラムを用い、以下のパラメータを入力して行った。
サブセット数：１０
繰り返し計算回数：４
使用フィルタ：３次元バターワースフィルタ
【００５４】
図８及び図９は、実施例１及び比較例１に用いたＳＰＥＣＴ画像であり、図９は、複数の断面各々のＳＰＥＣＴ画像を並べて表示したものである。図１０は、図８のＳＰＥＣＴ画像の実施例１における補正画像であり、図１１は、図９のＳＰＥＣＴ画像の実施例１における補正画像であり、図１２は、図９のＳＰＥＣＴ画像の比較例１における補正画像である。
【００５５】
図８〜図１２参照すれば明確なように、比較例１に対して、バックグランドノイズが低減され、病変部のコントラストが増加し、病変部の輪郭が鮮明な補正画像が、実施例１では得られた。
【図面の簡単な説明】
【００５６】
【図１】図１は、本発明の実施の形態に係る画像補正装置の構成を示す図である。
【図２】図２は、本発明の実施の形態に係る画像補正方法のフローチャートである。
【図３】図３は、吸収補正マップを生成する処理のフローチャートである。
【図４】図４は、吸収補正マップを生成する処理の他の実施の形態のフローチャートである。
【図５】図５は、本発明の実施の形態に係る画像補正プログラムの構成を、記録媒体と共に示す図である。
【図６】図６は、記録媒体に記憶されたプログラムを実行するためのコンピュータのハードウェア構成を示す図である。
【図７】図７は、記録媒体に記憶されたプログラムを実行するためのコンピュータの斜視図である。
【図８】図８は、実施例１及び比較例１に用いたＳＰＥＣＴ画像である。
【図９】図９は、実施例１及び比較例１に用いたＳＰＥＣＴ画像である。
【図１０】図１０は、図８のＳＰＥＣＴ画像の実施例１における補正画像である。
【図１１】図１１は、図９のＳＰＥＣＴ画像の実施例１における補正画像である。
【図１２】図１２は、図９のＳＰＥＣＴ画像の比較例１における補正画像である。
【符号の説明】
【００５７】
１０…画像補正装置、１２…入力部、１４…吸収補正マップ生成部、１６…補正画像生成部、１８…出力部。
【図１】
2006047152.tif 000003
【図２】
2006047152.tif 000004
【図３】
2006047152.tif 000005
【図４】
2006047152.tif 000006
【図５】
2006047152.tif 000007
【図６】
2006047152.tif 000008
【図７】
2006047152.tif 000009
【図８】
2006047152.tif 000010
【図９】
2006047152.tif 000011
【図１０】
2006047152.tif 000012
【図１１】
2006047152.tif 000013
【図１２】
2006047152.tif 000014