(57)【特許請求の範囲】
【請求項１】
被検体を捉えた複数枚のコンピュータ断層像に係る画像データにアクセスするステップと、
それらコンピュータ断層像のうち高密度部分が１個又は複数個写っている高密度部分含有像を識別する像識別ステップと、
高密度部分含有像それぞれから１個又は複数個の高密度部分を検出する検出ステップと、
それら高密度部分を一種類又は複数種類の構成部分特性に従い分類する分類ステップと、
典型的な組織内画素値を分類済の高密度部分上に置換することで１枚又は複数枚の補償像を生成する補償像生成ステップと、
前記１枚又は複数枚の補償像に基づいて補償サイノグラムＰを演算するステップと、
分類済の高密度部分についてのサイノグラムＭを補正前のコンピュータ断層像についてのサイノグラムＷから減ずることで高密度部分含有像毎に差分サイノグラムＷｓＭを生成するステップと、
前記補償サイノグラムＰを差分サイノグラムＷｓＭに加えることで、高密度部分含有像毎に最終的なサイノグラムＷｓＭａＰを生成するステップと、
生成した最終的なサイノグラムＷｓＭａＰに従い高密度部分含有像毎に補正版の再構成コンピュータ断層像を生成するステップと、
を有する補正版再構成コンピュータ断層像生成方法。
【請求項２】
請求項１記載の補正版再構成コンピュータ断層像生成方法であって、その像識別ステップにて、補正前の再構成コンピュータ断層像を表示させ、それに従いオペレータが入力する１個又は複数個の指令を受け付ける補正版再構成コンピュータ断層像生成方法。
【請求項３】
請求項１記載の補正版再構成コンピュータ断層像生成方法であって、その像識別ステップにて、被検体を捉えた複数枚のコンピュータ断層像のなかから１枚又は複数枚を候補像に指定した上で、
ａ）指定された候補像内で画素値を並べて順序統計量を導出するサブステップと、
ｂ）その順序統計量から求まる一種類又は複数種類の数値に従い、指定された候補像のうち高密度部分が１個又は複数個写っている像を識別するサブステップと、
ｃ）指定された候補像それぞれについてサブステップａ）及びｂ）を実行させるサブステップと、
を含む手順を実行する補正版再構成コンピュータ断層像生成方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明はＸ線ＣＴ（コンピュータ断層撮像）時画像処理、特に金属等の高密度物質で形成された高密度物体が原因で生じる密度性偽像を像の再構成時に抑圧する技術に関する。
【背景技術】
【０００２】
ＣＢＣＴ（コーンビームＣＴ）、スパイラルＸ線ＣＴ等のＣＴで発生しがちな偽像は多様であるが、非特許文献１の記載によれば物理性、患者性、スキャナ性の三種類に分類することができる。物理性偽像はビームハードニング、フォトンスターベーション、アンダーサンプリング等が原因で生じるもの、患者性偽像は高密度物体、不測体動等が原因で生じるもの、またスキャナ性偽像は検出器感度、機械的不安定性等が原因で生じるもののことである。
【０００３】
なかでもＣＴシステムの性能上大きな問題となるのは、金属等の高密度物体が原因で生じる密度性偽像である。そうした高密度物体で形成された部分を有する被検体をスキャンした場合、高密度部分の像から放散する髭といった形態の高周波ストリーク付ＣＴ像が得られることがある。これは、高密度部分における強い減衰がもとで、ＣＴシステム内検出器への光子到達個数が減少するために生じるものであり、信号対雑音比の低下にもつながっている。加えて、Ｘ線ビームの透過減衰率には波長依存性があるため、高密度部分を透過したＸ線ビームはハードニングしていることが多い。それによる非線形的なプロジェクション情報内変動は、高密度部分を取り巻き又は高密度部分同士を結ぶ低周波テイル偽像となって現れることがある。これらの密度性偽像が医用診断撮像で問題になるのは、被検体内の高密度部分によって高周波ストリーク等が引き起こされると、その周囲にある軟組織や骨格構造を調べるのが難しくなるためである。例えば歯牙ＣＢＣＴ撮像で得られたＣＴ像に、歯牙充填材（詰め物；フィクスチャ（人工歯根）等も含む）による偽像が生じていると、その周囲にある歯牙、骨、組織構造等を適正に調べるのが不可能又は困難となる。撮像したい被検体例えば患者の身体に備わる高密度部分でＸ線ビームがこうして様々に減衰するので、撮像で得た一群の断層像に基づき立体的なＣＴ像を形成する処理（いわゆる再構成）が複雑になってしまう。それでもなお不要な偽像が残ることがある。
【０００４】
ＣＴ像再構成の際に高密度物体の影響を抑える策は既に幾種類か検討されている。第１に、Ｘ線減衰率が低い充填／修復材や断面積が小さいインプラント／装具を開発して使用する、という単純な予防策がある。第２に、高密度部分によるプロジェクション情報損失を抑えるため浸透性が高い強エネルギＸ線ビームを使用する、という策がある。これらの策は、金属等で形成されている高密度部分の影響を軽減又は抑圧するのに有効だが、かといって万能なものではない。特に、Ｘ線ビームのエネルギを強める策は、患者にとってのリスクが増す点で貧弱な策である。
【０００５】
そのため、密度性偽像に対し画像処理技術で対処する策が開発されている。その代表は、体軸横断ＣＴ像のデータを形態変換で補間し新たな体軸横断像、冠状／矢状断像又は体軸斜交像を合成する、という策である。スキャンした平面以外の面へと画像データを形態変換すると、元々の体軸横断ＣＴ像が隣り合うスライス（断層像）同士で結合して平均化されるので、概ねランダムに散在している偽像信号に比べ、真正な画像信号が強調されることとなる。このほか、雑音を含む像に対しポスト再構成フィルタリングを直接適用して画質を高める、という策もある。
【０００６】
偽像の抑圧及び失われたプロジェクション情報の再生を目的とする画像処理手法には、大別して、プロジェクション情報補間法と反復的再構成法がある。反復的再構成法の一例は、特許文献１（名称：コンピュータ断層像内偽像抑圧方法及び装置(Method and Apparatus for Reduction of Artifacts in Computed Tomography Images)，発明者：Ｍａｎ）に記載のＣＴ像内偽像抑圧方法である。この方法は、個々のサイノグラム要素に係る荷重に基づく補正でサイノグラムを反復的に再構成することで、より秀逸なＣＴ像を再構成しようというものである。補正版のサイノグラムは、例えば、ＣＴスキャナで得られた当初のサイノグラムに対する補間によって生成される。
【０００７】
他方のプロジェクション情報補間法は、反復的でないサイノグラム再構成で同じ問題に対処する手法である。その一例は、特許文献２（名称：Ｘ線走査装置で生成された画像における金属性偽像を抑圧する方法及びシステム(Method and System for Reducing Metal Artifacts in Images Generated by X-ray Scanning Devices)、発明者：Ｎａｉｄｕ  ｅｔ  ａｌ．）に記載の方法である。この方法では、所与の像に含まれる高密度物体を識別し、識別した高密度物体からその物体のプロジェクションを生成し、生成したプロジェクションを所与のプロジェクションから減じて補正版のプロジェクションを生成し、そしてその補正版のプロジェクションから最終的な補正版の像を再構成する。同文献には、薄いシート状高密度物体を保持することも記載されている。
【０００８】
これら、従来型の画像処理手法では、密度性偽像の問題に対処するため、識別した高密度物体の近隣から密度性偽像抑圧に役立つ情報を抽出して利用している。この手法には確かに一理があるが、歯牙ＣＢＣＴ像に現れる偽像を正確に補償するには十分なものではない。ことに、歯牙撮像分野には、他種撮像分野では生じないその分野特有の問題が数多くある。例えば、他種ＣＴ撮像分野と違い歯牙ＣＴ撮像では、骨、象牙質、エナメル質、歯肉組織等、様々な種類の硬組織及び軟組織が写り、また充填材、インプラント、歯冠その他の修復材を初めとする様々な高密度物体も写る。金属等で形成された高密度物体が概ね均質な領域内に埋入している場合であれば、特許文献２に記載のプロジェクション情報補間法等も使用可能であるが、周囲の組織に対し大きく異なる特性を有する諸領域内に充填材等の高密度物体が存在している歯牙撮像環境では、そうした補間法では好適な結果を得ることができない。そして、自動金属検出処理は総じて情報処理依存的であり、その動作を最適化するため患者関連の情報を利用することもできない。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】米国特許第７０２３９５１号明細書
【特許文献２】米国特許第６７２１３８７号明細書
【非特許文献】
【００１０】
【非特許文献１】Mehran Yazdi, “Artifacts in Spiral X-ray CT Scanners: Problems and Solutions", Proceedings of World Academy of Science, Engineering and Technology, Volume 26, December 2007, pp. 376-380
【発明の概要】
【発明が解決しようとする課題】
【００１１】
このように、発明者の知見によれば、歯牙ＣＴスキャンを初めとするＣＢＣＴ等のＣＴ撮像時に使用でき、金属等で形成された高密度物体が原因で生じる偽像をより好適に補償可能な方法が求められている。本発明の目的は、ＣＴ像再構成における密度性偽像抑圧をより好適に行えるようにすることにある。
【課題を解決するための手段】
【００１２】
こうした目的に鑑み、本願では、少なくとも部分的にコンピュータ上のプロセッサで実行可能な再構成コンピュータ断層像補正方法であって、被検体を捉えた複数枚のＣＴ像に係る画像データにアクセスするステップと、それらのＣＴ像のうち高密度部分が１個又は複数個写っている高密度部分含有像を識別するステップと、高密度部分含有像それぞれから１個又は複数個の高密度部分を検出するステップと、それら高密度部分を一種類又は複数種類の構成部分特性に従い分類するステップと、典型的な組織内画素値を分類済の高密度部分上に割り振ることで１枚又は複数枚の補償像を生成するステップと、分類済の高密度部分についてのサイノグラムを補正前のＣＴ像についてのサイノグラムから減ずることで高密度部分含有像毎に差分サイノグラムを生成するステップと、生成済の補償像に従い生成したサイノグラムを差分サイノグラムに加えることで高密度部分含有像毎に最終的なサイノグラムを生成するステップと、高密度部分含有像毎に最終的なサイノグラムに従い補正版の再構成ＣＴ像を生成するステップと、を有する方法を提案する。
【００１３】
本方法は、偽像補償の対象となる個別の領域についての情報をオペレータが撮像システムに供給する形態で実施することができる。
【発明の効果】
【００１４】
本発明によれば、再構成ＣＴ像上の所要部位のみに直に密度性偽像補償処理を施すことができる。
【００１５】
本発明によれば、更に、システムを操作するオペレータの人智・技量を活用すると共に患者の受診歴に関する情報を利用し、コンピュータの密度性偽像抑圧性能を高めることができる。これにより、人間の長所である創造性、自発性、柔軟性及び判断力を、コンピュータの長所である情報処理速度、継続的精密処理能力及びデータ授受能力と結合させることができる。
【図面の簡単な説明】
【００１６】
【図１】本発明の一実施形態に係る処理を示す論理フローチャートである。
【図２】高密度部分含有像識別処理の一例を示す論理フローチャートである。
【図３】密度性偽像が生じている再構成ＣＢＣＴ像の一例を示す図である。
【図４】再構成ＣＢＣＴボリューム内のある候補像における最大画素値対順序統計量比曲面を示す図である。
【図５】高密度部分像の一例を示す図である。
【図６】高密度部分が写っている歯牙ボリュームの一例を示す図である。
【図７】図５に示した高密度部分像に基づき生成される補償像の一例を示す図である。
【図８】本発明の一実施形態に係る再構成方法を示す論理フローチャートである。
【図９Ａ】高密度部分補正が施されていない初期的なＣＴ像についてのサイノグラムを示す図である。
【図９Ｂ】高密度部分像についてのサイノグラムを示す図である。
【図９Ｃ】補償像についてのサイノグラムを示す図である。
【図９Ｄ】新たな像の再構成に用いられる補正版のサイノグラムを示す図である。
【図１０Ａ】密度性偽像が写っている歯牙ＣＴ像を示す図である。
【図１０Ｂ】図１０Ａの歯牙ＣＴ像に対し偽像補償を適用した結果を示す図である。
【図１１】本発明の他の実施形態に係るオペレータ操作利用型処理を示す論理フローチャートである。
【図１２Ａ】オペレータコンソールに備わるタッチスクリーン上で像の部分指定機能や処理対象患者情報入力機能に係る表示を行っている状態を示す平面図である。
【図１２Ｂ】オペレータコンソールに備わるタッチスクリーン上で分類用情報入力機能に係る表示を行っている状態を示す平面図である。
【発明を実施するための形態】
【００１７】
以下、本発明の好適な実施形態について別紙図面を参照しつつ説明する。以下の説明を図面と併せ参照することで、別紙請求の範囲に記載されている発明の目的、構成、効果及び諸実施形態について、上掲のものを含めてより深く理解することができよう。それらの説明で参照する図面では、同一の部材に対し同一の参照符号を付してある。本発明の実施形態に係る構成の機能や相互関係を総括的且つ例示的に示すための図であるので、実際の寸法又はその比に対し忠実なものでないことを了解されたい。
【００１８】
また、以下の説明中、「高密度部分」とは、周囲の組織に比べ高密度な充填材で埋められた部分のように、金属等の高密度物質で形成された物体、その物体が占める領域又はその領域の像のことである。「高密度物質」とは、その密度が約２．７ｇ／ｃｍ3を超えるアルミニウムのように、本件技術分野で習熟を積まれた方々（いわゆる当業者）が高密度と認める物質のことである。一般的な使用量に違いがあるため、高密度か否かを決める絶対的なしきい値を定めることは現実的ではない。像の種類によっては、アルミニウム等の金属に比べかなり低密度の物質が高密度物質として扱われる場合がある。以下の説明では金属質の高密度物質を念頭に置いているが、高密度部分の検出、その結果に基づく偽像の補償等といった処理は非金属質の高密度物質に対しても適用しうるものである。更に、「集合」とは、周知の集合論でいうところの元即ち集合構成要素を１個以上集めたもののことである。「部分集合」とは、特記がない限り、１個以上の元を有するより大きな集合（母集合）の一部分を形成していて１個以上の元を有する集合のことである。集合Ｓの部分集合がその集合Ｓそのものであることもある。「真部分集合」とは部分集合のうち母集合と等しくないもの、即ち母集合の元のうち１個以上がその集合に属していない部分集合のことである。説明中で使用している「第１」「第２」等の語は、その語が冠されている要素間の順序や優先度を示すためのものではなく、それらの要素を区別するためのものである。
【００１９】
そして、本発明の主題はディジタル画像処理技術及びコンピュータビジョン技術に関わっている。それらの技術は、画像データをディジタル的に処理することで、人間が理解可能な物体、属性又は状態を認識してそれに有用な意味合いを与え、その画像データで得られた結果を利用に供する技術である。ディジタル画像処理における画像形成プロセスは、一般に、画像データを生成する処理及びその画像データを（将来的な利用例えば更なる処理や表示に備え）電子メモリに保存する処理を含むプロセスとして理解されている。
【００２０】
ＣＴ撮像に当たり金属等の高密度部分で生じる諸問題への対策として、本発明では、従来と違う策を採用している。高密度部分が存する領域を分類し、その結果に基づき元々のプロジェクション像を補正し、補正版のプロジェクション像に基づきＣＴ像を再構成することで、高密度部分による偽像を抑圧する、という策である。
【００２１】
図１に、本発明の一実施形態に係る歯牙ＣＢＣＴボリューム内密度性偽像抑圧手順で実行される一連のステップをフローチャートにより示す。この手順では、まず、ＣＴスキャナで得られた複数枚のプロジェクション乃至サイノグラムを利用する標準的な再構成アルゴリズムに従い、且つ電子的に保存されている画像データに基づき、初期的な再構成ＣＴボリューム１０２を生成する。ボリュームとは、１枚又は複数枚の断層像からなる画像データのことである。通常、ボリューム１０２を構成している像のなかで、金属等で形成された高密度物体が写っている像は、そのボリューム１０２を構成している像の部分集合に過ぎない。寧ろ、ボリューム１０２を構成している像の多くには、高密度部分が写っておらず歯牙や軟組織が正確に捉えられている。従って、その実行に若干の時間や情報処理リソースが必要になる密度性偽像抑圧手順を、ボリューム１０２を構成している像全てについて実行する必要はない。そこで、像識別ステップ１０４ではそうした部分集合に属する像を識別・選出する。検出ステップ１０６では、ステップ１０４で選出された個々の像から高密度部分を検出する。次いで、従来手法と違い、検出された高密度部分によって占められている領域を分類ステップ１０８にて分類し、その結果に基づき偽像抑圧ステップ１１０にて偽像を抑圧する。この手順で最終的に得られるのは、画像処理によって偽像が抑圧された再構成ボリューム１１２である。
【００２２】
［高密度部分含有像の識別］
本実施形態では、偽像抑圧の効率及び実効性を全般的に高めるため、その対象を、ＣＴ像再構成用の像のうち高密度部分が写っている像（高密度部分含有像）に絞っている。高密度部分含有像の識別は、次の通り自動処理で実行することも、また後述の通りオペレータからシステムに与えられる指令に従い実行することも可能である。
【００２３】
図２に、上述した像識別ステップ１０４を自動処理で実行する際に実行される諸ステップの詳細をフローチャートにより示す。この手順では、初期的な再構成ＣＴボリューム１０２に基づき順序統計量を導出・利用している。順序統計量は人口構成解析等のデータ解析問題を解くのに役立つ非パラメトリック統計ツールであり、標本値をその値順に並べ直すことで導出することができる。本実施形態では、ボリューム１０２を構成している個々の像乃至スライスを抜取標本値の集合と見なし、そのスライスを構成している諸画素をその強度値ｋ（但しｋ＝０，…，ｎ）即ち画素値の順に並べ直すことで、そのスライスにおける順序統計量を導出する。更に、こうして導出した順序統計量を処理することで像を解析し、金属等で形成されている高密度物体がその像に写っているか否かを判別する。
【００２４】
図２に示した手順は、初期的な再構成ＣＴボリューム１０２を構成している個々の像を高密度部分含有像かもしれない候補像として扱い、それら候補像それぞれについて実行される。その手順のうち順序統計量導出ステップ２０２では、手始めに、個々の候補像を構成している諸画素をその画素値に従いソートして（並べ替えて）順序統計量の形態へと変換する。同ステップ２０２では、後続の処理に備え、そうした順序統計量を複数種類導出する。本実施形態で順序統計量として導出されるのは、例えばその候補像における画素の最大値（像内最大画素値）やメディアン値（像内メディアン画素値）である。その次の比算出ステップ２０４ではそれら順序統計量同士の比を算出する。そして、評価ステップ２０６では、順序統計量同士の比、例えば像内最大値に対する他の順序統計量（像内最大画素値自体も含む）の比を並べ、像のコンテンツに関する有用な情報が得られる候補像を選別する。即ち、高密度部分を含むと評価された候補像を高密度部分含有像として選出する。
【００２５】
図３〜図５に、初期的な再構成ＣＴボリューム１０２から像を選出する自動処理がどのようにして実行されどのような像が選出されるかを例により示す。まず、図３に示した再構成ＣＴ像７０２はボリューム１０２に属する候補像の例であり、この歯牙像には歯牙充填材によって生じた偽像が写っている。図２に示した手順では、こうした像７０２における順序統計量がステップ２０２で導出される。例えば、ボリューム１０２を構成している個々のスライス内でその強度値順に画素を並べ替えることで、像内最大画素値、像内メディアン画素値等の順序統計量が導出される。その次のステップ２０４では、導出された個々の順序統計量が像内最大画素値に対し有している比が算出される。
【００２６】
図４として示した立体的なグラフ１００２は、再構成ＣＴボリューム１０２中の第１〜第４７０スライスに関し、個々のスライスにおける順序統計量（但し像内メディアン画素値以上の値を有するものに限る）が像内最大画素値に対し有している比を求め、それらの比で形成される面の一部を描いたものである。このグラフ１００２の縦軸は比の算出結果に係る軸、０〜５００なる数字が振られている軸はボリューム１０２内スライス番号に係る軸である。縦軸に直交していて１０００、２０００及び３０００なる数字が振られている軸は、縦軸沿いに示した比をもたらす順序統計量の番号（順序）ｋに係る軸である。図示の通り、算出された比はスライス毎に曲線を形成している。この図でその曲線のうち比が大きな部分のみを示したのは、像内メディアン画素値より小さな順序統計量が像内最大画素値に対し有している比に、あまり意味がないからである。また、それらの曲線は、図示の通りいずれも１から始まっており、画素を増やすにつれ比の値が減少する傾向の曲率を呈している。留意すべきことに、ボリューム１０２内にその像とそっくりな像があるのでない限り、それらスライス間では一般にその曲線が異なった形状となる。
【００２７】
図４を用いビジュアル的に示した解析手法であれば、再構成ＣＴボリューム１０２を構成しているスライスのなかから、金属等の高密度物質で形成された部分が写っているものを好適に検出することができる。即ち、比で形成される曲面の窪みを以て高密度部分付のスライスを検出することができる。例えば、図４中で第９０スライス付近から第３００スライス付近に亘っている窪みは、ボリューム１０２のうちＸ線が透過しにくい歯牙充填材存在部分に相当している。また、こうした解析法を用い高密度部分含有像を自動検出する場合、像内最大画素値と像内メディアン画素値の比をしきい値と比較することで、図２中の評価ステップ２０６を実行することができる。即ち、像内最大画素値と像内メディアン画素値の比がしきい値未満（例えば０．２未満）の像を、金属等で形成された高密度部分を捉えている蓋然性が非常に高い像、即ち高密度部分含有像というに相応しい候補像であると評価し、後続の諸ステップで処理させることができる。図３に示した再構成ＣＴ像７０２は、歯牙充填材を捉えていてそうした評価が下される歯牙像の一例である。なお、自明な通り、ある候補像を高密度部分含有像と認定すべきか否かを迅速に判別可能な手法は数多くある。一種類又は複数種類の順序統計量について比を算出する手法はその一種に過ぎない。また、比に基づき高密度部分の存否を調べる際のしきい値は、経験的には０．２程度とするのが望ましいが、それとは違う値にすることもできる。上述した手法以外の手法で、ボリューム１０２内の個々のスライスにて画素値の概要を捉え評価することもできる。
【００２８】
［高密度部分の検出］
また、前掲の図１では、高密度部分含有像が選出された後に幾つかのステップが実行されている。最初に実行されているのは、ステップ１０４で選出された像を処理し、金属等で形成されている高密度部分を検出するステップ１０６である。検出される高密度部分は、当該高密度部分に属する画素からなる画素の集合（但し空集合でないもの）である。その集合には、その高密度部分内の全画素が同じ高密度部分内の他の画素から所定距離（例えば１．０画素）以内にある、という特性がある。
【００２９】
本実施形態では、金属等で形成された高密度部分を検出する手法としてしきい値判別を採用している。これは、像を構成している画素、例えば図３にて再構成ＣＴ像７０２を構成している画素の強度値が所定の強度しきい値を上回っているか否かを調べ、上回っている画素を高密度部分像の一部として保存する手法である。使用する強度しきい値は、例えば、像内最大画素値に対し所定の比を呈するよう像毎に設定する。その比を例えば９０％にすると、像内最大画素値に対し０．９倍の強度を有する画素が高密度部分像の一部として保存されることとなる。図５に示した高密度部分像７０４は、図３に例示した像７０２に対し、そうした検出ステップ１０６を適用して得られたものである。一見して理解できるように、この像７０４には金属等（ここでは歯牙充填材）で形成された高密度部分が複数個写っている。それらの高密度部分は、後続の処理ステップにて個別に処理することができる。
【００３０】
［分類］
歯牙ＣＢＣＴに従来型金属性偽像抑圧手法を適用しづらかった事情の一つは、歯牙撮像に関わる組織が多種多様なことである。背景技術の欄で述べた通り、金属等で形成された高密度物体によって引き起こされる偽像を抑圧するため、従来手法ではその高密度物体の近隣にある画素から抽出された情報を利用している。即ち、ミスプロジェクション（高密度物体によって妨げられているＸ線伝搬路に沿ったインテグレーション）を、その隣にある別のプロジェクション（高密度物体によって妨げられていないＸ線伝搬路に沿ったインテグレーション）によって補間している。こうした手法が好適に機能するのは、組織を被覆又は置換している高密度物体が、その高密度物体の直近を囲んでいる別の組織と同じ特性を有している状況下である。これに対し、歯牙充填の場合には、背景技術の欄で既に述べた通りそうした状況にならないことがある。その一例は歯根管充填である。その場合、組織（例えば歯根管）を被覆又は置換している充填材とその周囲にある組織（例えば象牙質）とが全く異質であるため、近隣のプロジェクションによる補間処理を実行すると、そのプロジェクション内にある密度的に不整合な要素が補間に使用されてしまうことになりかねない。
【００３１】
図６に、そうした問題を解決するのに分類が役立つことを示すため、歯牙ＣＴボリュームの例を立体的に示す。この立体表示９０２には、金属で形成された二種類の高密度部分が現れている。第１の高密度部分９０６は例えば象牙質充填材の階級、第２の高密度部分９０４は例えば歯根管充填材の階級に分類されうる部分である。そうした分類が可能なのは、高密度部分間に寸法、位置、形状等といった構成部分特性の差があるからである。即ち、高密度部分に備わる種々の構成部分特性を計測し、その結果を分類条件として利用することで、高密度部分を自動的に分類することができる。分類条件として利用できる特性の例としては、その部分の直径、長さ及び向きや、そのボリューム乃至歯牙を構成している他の部分乃至領域に対するその部分の位置関係がある。
【００３２】
なお、検出済の高密度部分を個別に調べる処理は、位置、形状等に関わる特性の計測結果に基づく自動処理以外に、オペレータから与えられる情報に基づく処理でも実行することができる。いずれにしても、個々の高密度部分を分類することで、その部分に相応しい偽像補償処理をその部分に適用することが可能になる。
【００３３】
［補償像の生成］
高密度部分含有像の選出、高密度部分像７０４の検出及び高密度部分の分類が終わったら、次は、偽像抑圧に使用される補償像を生成する。図７に、像７０４内該当部分に対する画素値割り振りで形成された補償像８０４を示す。割り振られる画素値は典型的な組織内画素値、例えば象牙質でよく現れる画素値である。象牙質によく現れる画素値を示す標本を予め複数枚の歯牙像から収集しておき、その標本を利用し像７０４内該当部分を置換することで、こうした補償像８０４を生成することができる。例えば、歯牙像標本保存ライブラリから歯牙像画素の標本をパッチとして切り出し、像７０４内該当部分に対応する領域内にその標本をランダムに埋め込むことで、像８０４を生成することができる。補償像内で使用される典型的な置換画素値としては、こうした例に限らず、統計的に習得した値を利用することができる。
【００３４】
［偽像補償を伴う再構成］
図８に、本実施形態におけるＣＴ像再構成時偽像補償手順をフローチャートにより示す。この手順は、図１に示した偽像抑圧ステップ１１０における処理動作に概ね対応している。その実行に際しては、前述の初期的な画像処理を実行しておくことが必要である。それには、金属等からなる高密度部分が写っている像（例えば初期的な再構成ＣＴ像７０２）を識別する処理や、識別した個々の高密度部分含有像毎に高密度部分像（例えば７０４）ひいては補償像（例えば８０４）を生成する処理が含まれる。
【００３５】
この図に示す手順中、ステップ３０２では図９Ａに示すサイノグラムＷ（１２０２）を生成する。この処理は、与えられている一群の像乃至スライスから選出された像（例えば初期的な再構成ＣＴ像７０２）毎に、相応の画像処理ホスト上でのラドン変換によって実行する。サイノグラム像形成のためのラドン変換及びその逆変換に係る情報処理は、医用画像処理の分野ではいわゆる当業者にとり周知なものである。ステップ３０４では、図９Ｂに示すサイノグラムＭ（１２０４）を生成する。この処理は、与えられている一群の像乃至スライスから選出された像（例えば図３及び図１０Ａ中の像７０２）に係る画像データから検出された個々の高密度部分像（例えば７０４）毎に実行する。同様に、ステップ３０６では、図９Ｃに示すサイノグラムＰ（１２０６）を生成する。この処理も、一群の所与像から選出済の像（例えば７０２）に基づき生成された補償像（例えば８０４）毎に実行する。減算ステップ３０８では、サイノグラムＷ（１２０２）からサイノグラムＭ（１２０４）を減じて差分サイノグラムＷｓＭを生成する。加算ステップ３１０では、補償像についてのサイノグラムＰ（１２０６）をそのサイノグラムＷｓＭに加算し、図９Ｄに示す最終的なサイノグラムＷｓＭａＰ（１２０８）を生成する。処理ステップ３１２では、補正版のＣＴ像を再構成するため、そのサイノグラムＷｓＭａＰに対しあらゆる角度から逆ラドン変換を適用する。このようにすると、サイノグラムＷｓＭａＰの構成要素のうちサイノグラムＷ又はＰ内の帯に対応する帯内にあるものがそのサイノグラムＷｓＭａＰ内で帯境界直近に存する構成要素で補間され、高密度部分像７０４の構成要素が同ステップ３１２の終了前処理にて最終再構成像１１０２（図１０Ｂ参照）内の対応位置に再移植され、その像１１０２内の画素のうち境界にあるものが平滑されるので、処理ステップ３１２の実行により補正版の再構成ＣＴ像を生成することができる。
【００３６】
図１０Ａ及び図１０Ｂに、本実施形態における密度性偽像抑圧効果を一例により示す。まず、初期的な再構成ＣＴ像７０２では、歯牙組織の像に係る領域１１０４内に、歯牙内高密度部分による偽像が複数現れている。この像７０２を図８に示した手順で処理すると、図１０Ｂに示すように、本実施形態の方法に従い対応する領域１１０６内で高密度部分の影響が補償され、偽像が抑圧された最終再構成像１１０２が得られる。即ち、この図の像１１０２では、元々の像７０２で高密度物体間の暗点になってしまっていた画素群が領域１１０６内で少なくとも部分的に復元されている。
【００３７】
［患者についての予備知識の利用］
構成部分特性に基づく分類に加え、患者の状態や受診歴から得られる情報も、偽像を検出及び抑圧するのに役立てることができる。即ち、その患者についての予備知識を参照することで、補償像（群）の生成ひいては失われた情報の復元をより好適に行うことができる。その予備知識例えば充填材の個数、位置及び組成は患者の受診歴等から知ることができる。
【００３８】
例えば、図３に示した再構成ＣＴ像７０２に対し検出ステップ１０６を適用すると、図５に示すような高密度部分像７０４が得られる。この像７０４には歯牙充填材等の高密度部分が複数個所に白く現れている。こうした高密度部分例えば金属質部分は、患者についての予備知識を利用することで、より詳細に分類することができる。例えば、図５及び図６に示した情報に基づく分類によって、或いは該当する部分の幾何学的特徴に基づく分類によって、ある部分９０４を歯根管充填材の階級、別の部分９０６を象牙質充填材の階級、といった具合に分類することができる。図５に示した像７０４であれば、その高密度部分が象牙質充填材の階級に分類され、その結果が後続のステップで利用されることとなる。
【００３９】
［オペレータの介在］
本発明は、オペレータの介在を通じ患者関連の有用な情報を取得及び適用し、高密度部分含有像の識別、高密度部分像の分類等を初め、諸機能に係る動作の効率を全体として向上させる形態にて実施することができる。
【００４０】
これは、金属等で形成されている高密度部分、即ちかなりの情報処理能力及びリソースを必要とする偽像補償処理を適用すべき部分が、一般にＣＢＣＴボリュームの一部を占めるに過ぎないからである。本実施形態では、患者についての予備知識を参照し、像のうち金属等の高密度部分を捉えている部分に偽像補償処理の対象を限定しているので、分類をより好適に実行し、また情報処理リソースをより効率的に利用することができる。
【００４１】
図１１にその際の処理をフローチャートにより示す。患者の歯牙に充填が施されていることをオペレータ例えば技能士が知っている場合、この図に示す手順を使用することで、偽像補償をより簡略に実行することができる。図中、ステップ５０２では、１枚又は複数枚の再構成ＣＴ像が入力される。その次の像識別ステップ５０４では、処理対象となる高密度含有像を識別する。この識別は、患者の受診歴から得られる予備知識に基づき効率的に実行することができる。更にその次の検出ステップ５０６では、識別された像内の高密度部分例えば歯牙充填材部分を検出する。
【００４２】
図１２Ａに、像識別ステップ５０４における高密度含有像識別や検出ステップ５０６における歯牙充填材部分検出に当たり利用可能なオペレータインタフェースを示す。これらの処理では、受診歴データベースにアクセスすることも、図中のインタフェースを介しオペレータから与えられる指令を利用することもできる。図示例では、オペレータコンソール１４０に設けられているタッチスクリーンディスプレイの画面上に、密度性偽像抑圧処理が施されていない再構成ＣＴボリューム１４２を表示させると共に、オペレータが手又はポインティングデバイスで操作可能な操作子１４６を複数個表示させている。ポインティングデバイスとしてはマウスその他のカーソル操作型デバイス等を使用することができる。オペレータは、そうした操作を通じ、表示されているボリューム１４２内の注目領域を括り、密度性偽像抑圧処理の対象に指定することができる。即ち、オペレータによる注目部分乃至領域指定操作を通じ、図１１中の像識別ステップ５０４における高密度部分含有像識別を効率的に実行することができる。オペレータが図中のメニュー型インタフェースを利用し、高密度部分例えば歯牙充填材で形成されている部分をリスト表示させ個別に指定する、といった操作で、同図中の検出ステップ５０６も効率的に実行することができる。
【００４３】
図１２Ｂに、高密度部分及びその周囲の組織を分類する処理に当たり、オペレータが表示を参照して指令を与えることができるよう、オペレータコンソール１４０上に再構成像１４８を表示させた状態を示す。なお、図１２Ａ及び図１２Ｂに示したインタフェースは単なる例であり、像識別ステップ５０４における高密度部分含有像識別、検出ステップ５０６における高密度部分例えば歯牙充填材の検出、並びに分類ステップ５０８における高密度部分の分類等、図１１中の諸ステップでは様々な種類のオペレータインタフェースユーティリティを利用することができる。分類、検出等といった処理にオペレータを介在させることで、再構成ステップ５１０にて生成される再構成像を、高密度部分による偽像がより好適に抑圧された像とすることができる。
【００４４】
このように、本発明の諸実施形態によれば、そのシステムを操作するオペレータの人智・技量をコンピュータの密度性偽像抑圧能力と協働的に結びつける実用的な密度性偽像抑圧システムが得られる。これにより、人間の長所である創造性、自発性、柔軟性及び判断力を、コンピュータの長所である情報処理速度、継続的精密処理能力、データ授受能力及び多彩表示機能と結合させることができる。
【００４５】
また、本発明は、一群のストアド命令からなるコンピュータプログラムに従い電子メモリ上の画像データに密度性偽像抑圧処理を施す形態、即ち一群のストアド命令に従いコンピュータを密度性偽像抑圧処理用プロセッサとして動作させる形態にて実施することができる。画像処理分野ではいわゆる当業者とり自明な通り、本発明を実施するためのコンピュータプログラムは相応の汎用コンピュータシステム、例えばパーソナルコンピュータやワークステーションで実行できるほか、プロセッサネットワークを含め様々なコンピュータシステム乃至プログラム式論理プロセッサで実行することができる。本発明に係る方法を実施するためのコンピュータプログラムは、また、様々な種類のコンピュータ可読記録媒体に格納しておくことができる。使用できる媒体としては、例えば、磁気ディスク（例えばハードディスクドライブ）、磁気テープ等の磁気記録媒体や、光ディスク、光テープ、機械可読バーコード等の光記録媒体や、ＲＡＭ（ランダムアクセスメモリ）、ＲＯＭ（リードオンリメモリ）等の固体電子記憶デバイスがある。コンピュータプログラムを格納可能な有形デバイス乃至媒体であればよい。更に、インターネットその他の通信媒体を介し画像プロセッサからコンピュータ可読記録媒体にアクセスし、その媒体上に格納されているコンピュータプログラムを実行することで、本発明に係る方法を実施する形態にしてもよい。いわゆる当業者にはご理解頂けるように、これらのコンピュータプログラムで実現される機能と同等の機能をハードウェアで実現することもできる。
【００４６】
また、ご理解頂けるように、本発明を実施するためのコンピュータプログラムを作成する際に使用する画像操作アルゴリズム乃至処理のうち幾つかは、周知のものでかまわない。本願中で説明していないものを含め、また画像処理分野にていわゆる当業者にとり既知のユーティリティを含め、その実現に役立つ様々なアルゴリズム乃至処理を利用することができる。本発明を実施するためのコンピュータプログラムと連携し像の生成乃至処理を担うアルゴリズム、システム、ハードウェア乃至ソフトウェアのうち、本願中で具体的に説明していないものは、本件技術分野で既知の要素で形成されたアルゴリズム、システム、ハードウェア乃至ソフトウェアとすることができる。条件を指定し検索することが可能な情報ライブラリ乃至データベースをコンピュータハードウェア上に構築し利用すること、例えば適切な組織像群でライブラリを構築しておき図７に示した補償像８０４の生成に利用することも可能である。
【００４７】
念のために述べると、本発明に係る補正版再構成ＣＴ像生成方法は、その分類ステップにて、補正前の再構成ＣＴ像を表示させ、それに従いオペレータが入力する１個又は複数個の指令を受け付ける構成にすることが可能である。
【００４８】
本発明に係る補正版再構成ＣＴ像生成方法は、その像識別ステップにて、画像データを並べ替えて順序統計量を導出する構成にすることが可能である。
【００４９】
本発明に係る補正版再構成ＣＴ像生成方法は、金属で形成されている高密度部分等に適用可能である。
【００５０】
本発明に係る補正版再構成ＣＴ像生成方法は、その補償像生成ステップにて、１枚又は複数枚の保存像から典型的な組織内画素値を取得する構成にすることが可能である。
【００５１】
本発明に係る補正版再構成ＣＴ像生成方法は、コーンビームＸ線スキャナで得られたＣＴ像を処理する構成にすることが可能である。
【００５２】
本発明に係る補正版再構成ＣＴ像生成方法は、その検出ステップにてしきい値を使用する構成にすることが可能である。
【００５３】
本発明に係る補正版再構成ＣＴ像生成方法は、その分類ステップにて、高密度部分の位置、寸法又はその双方を参照する構成にすることが可能である。
【００５４】
本発明に係る補正版再構成ＣＴ像生成方法は、その分類ステップにて、高密度部分の周囲にある組織がエナメル質、象牙質、歯根及び骨のいずれであるかを識別する構成にすることが可能である。
【００５５】
本発明に係る補正版再構成ＣＴ像生成方法は、補正版の再構成ＣＴ像を表示するステップを有する構成にすることが可能である。
【００５６】
本発明に係る補正版再構成ＣＴ生成方法は、その像識別ステップにて、被検体についての予備知識を参照する構成にすることが可能である。
【００５７】
本発明に係る補正版再構成ＣＴ像生成方法は、順序統計量の比を求め利用する構成にすることが可能である。
【００５８】
本発明に係る補正版再構成ＣＴ像生成方法は、候補内最大画素値に対する一種類又は複数種類の順序統計量の比を求め利用する構成にすることが可能である。
【００５９】
そして、本発明に係る補正版再構成ＣＴ像生成方法は、被検体を捉えた複数枚のＣＴ像に係る画像データにアクセスするステップと、それらＣＴ像のうち高密度部分が１個又は複数個写っている高密度部分含有像を順序統計量たる画像値の並びに従い識別するステップと、高密度部分含有像それぞれから１個又は複数個の高密度部分を検出するステップと、それら高密度部分を一種類又は複数種類の構成部分特性に従い分類するステップと、典型的な組織内画素値を分類済の高密度部分上に割り振ることで１枚又は複数枚の補償像を生成するステップと、分類済の高密度部分についてのサイノグラムを補正前のＣＴ像についてのサイノグラムから減ずることで高密度部分含有像毎に差分サイノグラムを生成するステップと、生成済の補償像に従い生成したサイノグラムを差分サイノグラムに加えることで高密度部分含有像毎に最終的なサイノグラムを生成するステップと、生成した最終的なサイノグラムに従い高密度部分含有像毎に補正版の再構成ＣＴ像を生成するステップと、補正版の再構成ＣＴ像を表示するステップと、を有する構成にすることが可能である。
【符号の説明】
【００６０】
１０２，１４２  再構成ＣＴボリューム、１０４，５０４  像識別ステップ、１０６，５０６  検出ステップ、１０８，５０８  分類ステップ、１１０，５１０  偽像抑圧ステップ、１１２  偽像抑圧済再構成ＣＴボリューム、１４０  オペレータコンソール、１４６  操作子、１４８，７０２  初期的な再構成ＣＴ像、２０２  順序統計量導出ステップ、２０４  比算出ステップ、２０６  評価ステップ、３０２〜３０６  サイノグラム生成ステップ、３０８  減算ステップ、３１０  加算ステップ、３１２  処理ステップ、５０２  再構成ＣＴ像入力ステップ、７０４  高密度部分像、８０４  補償像、９０２  歯牙ボリュームの立体的な表示、９０４，９０６  高密度部分、１００２  最大強度値対順序統計量比曲面の立体的なグラフ、１１０２  最終的な再構成ＣＴ像、１１０４，１１０６  像内領域、１２０２〜１２０８  サイノグラムＷ，Ｍ，Ｐ，ＷｓＭａＰ。
【図１】
0005858625.tif 000002
【図２】
0005858625.tif 000003
【図８】
0005858625.tif 000004
【図１１】
0005858625.tif 000005
【図３】
0005858625.tif 000006
【図４】
0005858625.tif 000007
【図５】
0005858625.tif 000008
【図６】
0005858625.tif 000009
【図７】
0005858625.tif 000010
【図９Ａ】
0005858625.tif 000011
【図９Ｂ】
0005858625.tif 000012
【図９Ｃ】
0005858625.tif 000013
【図９Ｄ】
0005858625.tif 000014
【図１０Ａ】
0005858625.tif 000015
【図１０Ｂ】
0005858625.tif 000016
【図１２Ａ】
0005858625.tif 000017
【図１２Ｂ】
0005858625.tif 000018