JP4138558B2 - 画像再構成装置、画像再構成方法および放射線断層像撮影装置 - Google Patents
画像再構成装置、画像再構成方法および放射線断層像撮影装置 Download PDFInfo
- Publication number
- JP4138558B2 JP4138558B2 JP2003100668A JP2003100668A JP4138558B2 JP 4138558 B2 JP4138558 B2 JP 4138558B2 JP 2003100668 A JP2003100668 A JP 2003100668A JP 2003100668 A JP2003100668 A JP 2003100668A JP 4138558 B2 JP4138558 B2 JP 4138558B2
- Authority
- JP
- Japan
- Prior art keywords
- data
- continuous
- region
- radiation
- image reconstruction
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 73
- 230000005855 radiation Effects 0.000 title claims description 62
- 238000003325 tomography Methods 0.000 title claims description 14
- 238000001914 filtration Methods 0.000 claims description 71
- 238000001514 detection method Methods 0.000 claims description 28
- 238000000926 separation method Methods 0.000 claims description 7
- 238000000605 extraction Methods 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 claims description 3
- 230000001678 irradiating effect Effects 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 230000002194 synthesizing effect Effects 0.000 claims description 2
- 239000000126 substance Substances 0.000 description 12
- 238000013480 data collection Methods 0.000 description 9
- 229910052751 metal Inorganic materials 0.000 description 9
- 239000002184 metal Substances 0.000 description 9
- 238000002591 computed tomography Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 8
- 238000010521 absorption reaction Methods 0.000 description 6
- 238000007781 pre-processing Methods 0.000 description 6
- 238000003745 diagnosis Methods 0.000 description 4
- RPPBZEBXAAZZJH-UHFFFAOYSA-N cadmium telluride Chemical compound [Te]=[Cd] RPPBZEBXAAZZJH-UHFFFAOYSA-N 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/448—Computed tomography involving metal artefacts, streaking artefacts, beam hardening or photon starvation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S378/00—X-ray or gamma ray systems or devices
- Y10S378/901—Computer tomography program or processor
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、たとえば、Ｘ線ＣＴ(Computed Tomography) 装置等の放射線断層像撮影装置に関する。
【０００２】
【従来の技術】
Ｘ線ＣＴ装置は、放射線出射系・検出系、すなわち、Ｘ線出射・検出系を被検体を挟んで相対向させた状態でＸ線を被検体に照射しながら被検体の回りで回転（スキャン(scan)）させて、被検体の周囲の複数のビュー(view)方向で被検体を透過したＸ線を検出系により検出し、それぞれＸ線による被検体の投影データを測定し、この投影データに基づいて、被検体の断層像を再構成するように構成されている。
【０００３】
また、被検体の断層像は、被検体の投影データに対して、所定の前処理を施し、この前処理を施した投影データを、たとえば、フィルタード・バックプロジェクション(filtered back projection)法により再構成される。
フィルタード・バックプロジェクションでは、被検体の投影データに対して画像再構成のために必要なフィルタリング処理を行ったのち、フィルタリング処理後の投影データを逆投影することにより、画像を再構成する。
上記のフィルタリング処理では、診断目的や対象部位に応じて各種のフィルタが選択される。このフィルタリング処理は、投影データを高速フーリエ変換(FFT:Fast Fourier Transfer) したのち、このデータに選択されたフィルタが乗算され、再び逆高速フーリエ変換(IFFT: Inverse FFT) される。フィルタリング処理では、投影データに含まれる高周波のノイズ成分も同時に除去される。なお、投影データに含まれる高周波のノイズ成分を除去する技術は、たとえば、特許文献１に開示されている。
【０００４】
【特許文献１】
特開２０００−８３９４６号公報
【０００５】
【発明が解決しようとする課題】
ところで、金属等の人体とは異なる異物質が被検体に存在すると、検出された投影データのうち、この異物質に相当する部分の値が急峻に変化し、この部分は他の部分とは非連続な部分となる。金属は、Ｘ線の吸収率が非常に大きいからである。
金属等の異物質の像を含む投影データに上記のフィルタリング処理を施すと、異物質に対応する部分のデータは値が急峻に変化する高周波成分をもつので、この部分のデータ値が変化し、逆高速フーリエ変換したのちの投影データにはアンダーシュートやオーバシュートと呼ばれる歪みが発生することがある。
この歪みが発生した投影データを逆投影して画像再構成すると、金属等の異物質に対応する画像にいわゆるストリーク・アーチファクト(streak artifact) が発生するという不利益が存在した。
【０００６】
本発明は、上述の問題に鑑みて成されたものであって、その目的は、金属等の放射線吸収率の高い物質が被検体に存在することにより再構成画像に発生するストリーク・アーチファクトを低減可能な放射線断層像撮影装置、画像再構成装置および画像再構成方法を提供することにある。
【０００７】
【課題を解決するための手段】
本発明の画像再構成装置は、被検体を透過した放射線を複数の放射線検出素子によって検出して得られた投影データに基づいて、前記被検体の断層像を再構成する画像再構成装置であって、前記放射線検出素子の配列順に並ぶ前記投影データのうち、非連続的に値が異なるデータを含む非連続領域を除くデータに対して所定のフィルタリング処理を施こし、前記非連続領域のデータと前記フィルタリング処理が施されたデータとを前記断層像を再構成するための投影データとする非連続データ処理手段を有する。
【０００８】
好適には、前記非連続データ処理手段は、前記非連続領域のデータを前記投影データから分離し、当該投影データを値が連続的に変化する連続データへ変換し、当該連続データのみを前記フィルタリング処理の対象にする分離手段と、前記フィルタリング処理された前記連続データと前記非連続領域のデータとを合成し、前記画像再構成手段に与える投影データとする合成手段とを有する。
【０００９】
前記非連続データ処理手段が、前記フィルタリング処理されたのちの投影データと、前記フィルタリング処理されていない投影データとの差に基づいて、前記非連続領域を抽出する抽出手段と、前記フィルタリング処理された投影データのうち、抽出された前記非連続領域のデータをフィルタリング処理前の投影データに置き換え、前記画像再構成手段に与える投影データとする合成手段と
を有する構成とすることも可能である。
【００１０】
本発明の画像再構成方法は、被検体を透過した放射線を複数の放射線検出素子によって検出して得られた投影データに基づいて、前記被検体の断層像を再構成する画像再構成方法であって、前記放射線検出素子の配列順に並ぶ前記投影データのうち、非連続的に値が異なるデータを含む非連続領域を除くデータに対して所定のフィルタリング処理を施こし、前記非連続領域のデータと前記フィルタリング処理が施されたデータとを前記断層像を再構成するための投影データとする。
【００１１】
本発明の放射線断層像撮影装置は、被検体を挟むようにして相対向して配される放射線照射手段および複数の放射線検出素子からなる放射線検出器と、前記被検体を透過した放射線を前記複数の放射線検出素子によって検出して得られた投影データに基づいて、前記被検体の断層像を再構成する画像再構成装置とを有し、前記画像再構成装置は、前記放射線検出素子の配列順に並ぶ前記投影データのうち、非連続的に値が異なるデータを含む非連続領域を除くデータに対して所定のフィルタリング処理を施こし、前記非連続領域のデータと前記フィルタリング処理が施されたデータとを前記断層像を再構成するための投影データとする非連続データ処理手段を有する。
【００１２】
本発明では、投影データから断層像を再構成する前に行うフィルタリング処理において、投影データに含まれる非連続データにはフィルタリング処理を施さずにそのまま画像再構成に用いる。これにより、非連続データのもつ高周波成分がフィルタリング処理によって歪むことがなく、この非連続部分の画像が正確に再構成される。
【００１３】
【発明の実施の形態】
以下、本発明の実施の形態について図面を参照して説明する。
第１実施形態
図１は、本発明の放射線断層像撮影装置の一実施形態としてのＸ線ＣＴ装置の全体構成を示すブロック図である。
【００１４】
図１に示すように、Ｘ線ＣＴ装置１は、走査ガントリ(gantry)２、操作コンソール(console) ３、および撮影テーブル（クレードル）４を有している。
【００１５】
走査ガントリ２は、Ｘ線管２１を有する。Ｘ線管２１から放射された図示しないＸ線は、コリメータ(collimator)２２により、たとえば、扇状のＸ線ビーム(beam)となるように成形され、Ｘ線検出器２３に照射される。
【００１６】
Ｘ線検出器２３は、ファンビームＸ線の広がり方向にアレイ状に配列された複数のＸ線検出素子を有する。Ｘ線検出器２３は、複数のＸ線検出素子をアレイ状に配列した、多チャネル(channel) の検出器となっている。
【００１７】
Ｘ線検出器２３は、全体として、円筒凹面状に湾曲したＸ線入射面を形成する。Ｘ線検出器２３は、たとえば、シンチレータ(scintillator)とフォトダイオード(photo diode) の組み合わせによって構成される。なお、これに限られず、たとえば、カドミウム・テルル（ＣｄＴｅ）等を利用した半導体Ｘ線検出素子またはＸｅガス(gas) を用いる電離箱型のＸ線検出素子であっても良い。
【００１８】
Ｘ線検出器２３にはデータ収集部２４が接続されている。データ収集部２４は、Ｘ線検出器２３の個々のＸ線検出素子の検出データを収集する。Ｘ線管２１からのＸ線の照射は、Ｘ線コントローラ(controller)２５によって制御される。なお、Ｘ線管２１とＸ線コントローラ２５との接続関係およびコリメータ２２とコリメータコントローラ２６との接続関係については図示を省略する。コリメータ２２は、コリメータコントローラ２６によって制御される。
【００１９】
Ｘ線管２１、コリメータ２２、Ｘ線検出器２３、データ収集部２４、Ｘ線コントローラ２５およびコリメータコントローラ２６が、走査ガントリ２の回転部２７に搭載されている。ここで、被検体は、回転部２７の中心に位置するボア(bore)２９内のクレードル(cradle)上に載置される。回転部２７は、回転コントローラ２８により制御されつつ回転し、Ｘ線管２１からＸ線を爆射し、Ｘ線検出器２３において被検体の透過Ｘ線を、各ビューごとの投影情報として検出する。なお、回転部２７と回転コントローラ２８との接続関係については図示を省略する。
【００２０】
操作コンソール６は、中央処理装置３１、入力装置３２、表示装置３３および記憶装置３４を有する。
なお、中央処理装置３１および記憶装置３４が、本発明の画像再構成装置を構成している。
【００２１】
中央処理装置３１は、たとえば、マイクロプロセッサ、メモリ等によって構成される。
この中央処理装置３１は、記憶装置３４に記憶されたソフトウエアにしたがって、走査ガントリ２の動作を制御し、被検体を透過したＸ線をＸ線検出器２３で検出して得られる投影データを収集する機能と、収集したＸ線投影データに基づいて、被検体の断層像を再構成する機能とを少なくとも有する。
【００２２】
中央処理装置３１には、データ収集部２４で収集されたデータが図示しないデータ収集バッファを介して入力され、中央処理装置３１は、収集した投影データを用いて画像再構成を行う。この中央処理装置３１による画像再構成処理については後述する。
【００２３】
中央処理装置３１には、また、表示装置３３と入力装置３２がそれぞれ接続されている。表示装置３３は、中央処理装置３１から出力される断層画像情報やその他の情報を表示する。入力装置３２は、オペレータによって操作され、各種の指示や情報等を中央処理装置３１に入力する。オペレータは表示装置３３および入力装置３２を使用してインタラクティブ(interactive) に本装置を操作する。
【００２４】
上記構成のＸ線ＣＴ装置１では、走査ガントリ２において、Ｘ線管２１およびＸ線検出器２３とを相対的に移動させ、被検体の被検対象領域を走査することにり、データ収集部２４に被検体の投影データが収集され、この投影データが中央処理装置３１へ送られ、中央処理装置３１において被検体の断層画像が再構成される。
【００２５】
図２は、中央処理装置３１による画像再構成処理手順の概略を示す図である。中央処理装置３１における画像再構成は、フィルタード・バックプロジェクション(filtered back projection)法によって行われる。なお、フィルタード・バックプロジェクション法については周知の方法であるので詳細な説明については省略する。
【００２６】
図２に示すように、中央処理装置３１による画像再構成処理は、大きく分けて、前処理（ステップＳＴ１）と、分離処理（ステップＳＴ２）、フィルタリング処理（ステップＳＴ３）および合成処理（ステップＳＴ４）からなる非連続データ処理と、バックプロジェクション処理（ステップＳＴ５）とからなる。
【００２７】
前処理（ステップＳＴ１）では、正常な画像が再構成されるように、データ収集部２４から得られた被検体の投影データに、物理的および電気的な補正処理を施す。
【００２８】
ここで、図３は分離処理（ステップＳＴ２）の内容を説明するための図であり、（ａ）は、データ収集部２４から得られた被検体の投影データの一例を示しており、（ｂ）は分離処理後の投影データを示しており、（ｃ）は分離された非連続データを示している。なお、これらのＸ線検出器２３の各検出器アレイのチャンネル順（配列順）に並んでいる。
非連続データとは、他の連続的に変化するデータとは複数チャンネルにわたって値が異なるデータを含む部分的なデータである。
【００２９】
たとえば、金属等の人体とは異なりＸ線の吸収率が高い異物質が被検体に存在すると、図３（ａ）に示すように、投影データｄａｔ＿Ｒ中において、異物質に対応する部分のデータ値が急峻に変化する。投影データｄａｔ＿Ｒの異物質が存在していない部分は、連続的に滑らかにデータ値が変化する。
【００３０】
図３（ａ）に示す投影データｄａｔ＿Ｒ中の、値が急峻に変化するデータを含む非連続領域のデータを非連続データｄａｔ＿ＥＸとすると、この非連続データｄａｔ＿ＥＸは、データ値を周波数領域でみたときに高周波成分を含んでいる。非連続データｄａｔ＿ＥＸ以外の部分には、データ値が滑らかに変化するため、このような高周波成分は存在していない。
【００３１】
投影データｄａｔ＿Ｒに含まれる高周波のノイズ成分を除去するために、投影データｄａｔ＿Ｒの全てにフィルタリングを施すと、上記の非連続データｄａｔ＿ＥＸに含まれる高周波成分も除去される可能性がある。
このため、分離処理（ステップＳＴ２）においては、まず、図３（ｃ）に示すように、投影データｄａｔ＿Ｒから非連続データｄａｔ＿ＥＸを分離する。
【００３２】
また、非連続データｄａｔ＿ＥＸを分離したのちの投影データｄａｔ＿Ｒは、図３（ｂ）に示すように、非連続データｄａｔ＿ＥＸが存在していた部分をデータｄａｔ＿Ｃにより補完する。具体的には、非連続データｄａｔ＿ＥＸの周辺のデータを用いてフィッティングを行う。
これにより、すべての領域でデータ値が滑らかに変化する連続データｄａｔ＿Ｓが得られる。
この連続データｄａｔ＿Ｓは、次にフィルタリング処理される。
【００３３】
フィルタリング処理（ステップＳＴ３）では、連続データｄａｔ＿Ｓを、まず、高速フーリエ変換（ＦＦＴ）する。
次いで、診断目的や対象部位に応じて選択されたフィルタによってフィルタリング処理される。
次いで、フィルタリング処理された連続データｄａｔ＿Ｓを逆高速フーリエ変換（ＩＦＦＴ）して復元する。
【００３４】
フィルタリング処理（ステップＳＴ３）におけるフィルタは、診断目的や対象部位に応じて適宜選択されるが、連続データｄａｔ＿Ｓに含まれる高周波のノイズ成分も同時に除去する。
ここで、図４は、フィルタリング処理（ステップＳＴ３）におけるフィルタの特性を説明するためのグラフであり、（ａ）はＦＦＴ前のデータを示しており、（ｂ）はＩＦＦＴ後のデータを示している。
【００３５】
図４（ａ）からわかるように、ＦＦＴ前のデータにはノイズ成分が存在しているが、フィルタリング処理によってこれらのノイズは除去される。
また、データ中に値が急峻に変化する部分が存在すると、ＩＦＦＴ後のデータには、図４（ｂ）の円Ａ内に示すように、この部分周辺においてオーバシュートやアンダーシュートと呼ばれる歪みが発生する。
【００３６】
本実施形態では、上記した分離処理（ステップＳＴ２）において、高周波成分をもつ非連続データｄａｔ＿ＥＸを分離してフィルタリング処理の処理対象から除外している。このため、連続データｄａｔ＿Ｓは、フィルタリング処理によってノイズ成分だけが除去され、歪みが発生することがない。
【００３７】
合成処理（ステップＳＴ４）では、分離された非連続データｄａｔ＿ＥＸと、フィルタリング処理された連続データｄａｔ＿Ｓとを合成し、この合成されたデータをバックプロジェクション処理（ステップＳＴ５）の画像再構成に用いるデータとする。
【００３８】
以上のように、本実施形態では、投影データｄａｔ＿Ｒに含まれる非連続データｄａｔ＿ＥＸを分離して得られる連続データｄａｔ＿Ｓのみにフィルタリング処理を施す。そして、フィルタリング処理を施した連続データｄａｔ＿Ｓと非連続データｄａｔ＿ＥＸとを合成し、この合成したデータをバックプロジェクション処理に用いる。これにより、再構成された画像に、フィルタリング処理による高周波成分の除去に起因するストリーク・アーチファクトが発生しにくくなり、金属等のＸ線の吸収率の高い異物質の断層像を正確に再現することが可能となる。
【００３９】
第２実施形態
図５は、本発明の放射線断層像撮影装置の他の実施形態に係る画像再構成処理手順の概略を示す図である。
なお、本実施形態は、上述した中央処理装置３１における画像再構成処理手順以外は第１の実施形態の構成と同じである。
【００４０】
図５に示すように、本実施形態に係る中央処理装置３１による画像再構成処理は、大きく分けて、前処理（ステップＳＴ１１）と、フィルタリング処理（ステップＳＴ１２）、差分処理（ステップＳＴ１３）、抽出処理（ステップＳＴ１４）および合成処理（ステップＳＴ１５）からなる非連続データ処理と、バックプロジェクション処理（ステップＳＴ１６）とからなる。
【００４１】
前処理（ステップＳＴ１１）は、上述した第１の実施形態と同様の処理を行う。
【００４２】
フィルタリング処理（ステップＳＴ１２）は、上述した第１の実施形態と同様のフィルタリング処理を行う。すなわち、データ収集部２４から得られた被検体の投影データに、まず、高速フーリエ変換（ＦＦＴ）し、診断目的や対象部位に応じて選択されたフィルタによってフィルタリング処理し、フィルタリング処理されたデータを逆高速フーリエ変換（ＩＦＦＴ）して復元する。
【００４３】
ここで、図６（ａ）は、フィルタリング処理前の投影データの一例であり、（ｂ）はフィルタリング処理後の投影データの一例であり、（ｃ）はフィルタリング処理前後の投影データの差分値を示している。
図６（ａ）に示すように、フィルタリング処理前の投影データｄａｔ＿Ｂにはノイズが存在するとともに、非連続的に値が異なるデータを含む非連続領域Ｒａが存在している。この非連続領域Ｒａは、上述したように、金属等のＸ線吸収率の高い異物質の存在により発生する。
【００４４】
図６（ａ）に示すような投影データｄａｔ＿Ｂに対してフィルタリング処理を施すと、たとえば、図６（ｂ）に示すデータｄａｔ＿Ａとなる。
図６（ｂ）に示すフィルタリング後のデータｄａｔ＿Ａには、非連続領域Ｒａにおけるデータ値が急峻に変化する位置で歪みｓｔが発生しているのが分かる。
【００４５】
次いで、差分処理（ステップＳＴ１３）では、フィルタリング処理前の投影データｄａｔ＿Ｂとフィルタリング処理後の投影データｄａｔ＿Ａとの差分を算出する。この差分値データは、図６（ｃ）に示すように、投影データｄａｔ＿Ｂに含まれるノイズと、フィルタリングによって発生した歪みｓｔを含んでいる。
【００４６】
次いで、抽出処理（ステップＳＴ１４）では、図６（ｃ）に示したような差分値データから、非連続領域Ｒｂ，Ｒｃを抽出する。この非連続領域Ｒｂ，Ｒｃは、投影データｄａｔ＿Ｂのうち、フィルタリングによって歪みｓｔが発生する領域である。
この非連続領域Ｒｂ，Ｒｃは、たとえば、所定のしきい値を設けて、差分値データとこのしきい値とを比較することにより抽出することができる。なお、差分値データとしきい値とを比較する際に、ノイズ成分と歪みｓｔとを区別する必要がある。歪みｓｔは複数のチャンネルにわたって存在し、ノイズ成分は単一のチャンネルに存在する場合が多く、また、歪みｓｔとノイズ成分のレベルは異なるので、両者の区別は容易である。
【００４７】
合成処理（ＳＴ１５）は、フィルタリング処理後の投影データｄａｔ＿Ａのうち、抽出処理（ステップＳＴ１４）によって抽出された非連続領域Ｒｂ，Ｒｃにおけるデータをフィルタリング処理前の投影データｄａｔ＿Ｂに置き換える。すなわち、非連続領域Ｒｂ，Ｒｃに歪みｓｔの発生していないデータを用いる。この合成されたデータをバックプロジェクション処理（ステップＳＴ１６）の画像再構成に用いるデータとする。
【００４８】
以上のように、本実施形態では、投影データｄａｔ＿Ｂに実際にフィルタリング処理を行って歪みｓｔを発生させることにより、非連続領域Ｒｂ，Ｒｃを特定し、この非連続領域Ｒｂ，Ｒｃにフィルタリング処理を行っていないデータを用いて画像再構成に用いるデータとする。これにより、バックプロジェクション処理によって再構成された画像にフィルタリング処理に起因するストリーク・アーチファクトが発生するのを確実に防ぐことができる。
【００４９】
なお以上に説明した各実施形態においては、放射線としてＸ線を用いた例について説明したが、放射線はＸ線に限るものではなく、たとえばγ線等の他の種類の放射線であっても良い。ただし、現時点では、Ｘ線がその発生、検出および制御等に関し実用的な手段が最も充実している点で好ましい。
【００５０】
【発明の効果】
本発明によれば、金属等の放射線吸収率の高い物質が被検体に存在することにより再構成画像に発生するストリーク・アーチファクトを低減可能な放射線断層像撮影装置、画像再構成装置および画像再構成方法が提供される。
【図面の簡単な説明】
【図１】本発明が適用される放射線断層像撮影装置としてのＸ線ＣＴ装置の全体構成を示すブロック図である。
【図２】本発明の一実施形態に係る中央処理装置の画像再構成処理手順を示す図である。
【図３】分離処理の内容を説明するための図である。
【図４】フィルターの特性を説明するためのグラフである。
【図５】本発明の画像再構成装置の他の実施形態に係る中央処理装置の処理手順を示す図である。
【図６】フィルタリング処理前、フィルタリング処理後および差分処理後の投影データの一例を示すグラフである。
【符号の説明】
１…Ｘ線ＣＴ装置
２…走査ガントリ
３…操作コンソール
４…撮影テーブル
２１…Ｘ線管
２２…コリメータ
２３…Ｘ線検出器
２４…データ収集部
２５…Ｘ線コントローラ
２６…コリメータコントローラ
３１…中央処理装置
３２…入力装置
３３…表示装置
３４…記憶装置
【発明の属する技術分野】
本発明は、たとえば、Ｘ線ＣＴ(Computed Tomography) 装置等の放射線断層像撮影装置に関する。
【０００２】
【従来の技術】
Ｘ線ＣＴ装置は、放射線出射系・検出系、すなわち、Ｘ線出射・検出系を被検体を挟んで相対向させた状態でＸ線を被検体に照射しながら被検体の回りで回転（スキャン(scan)）させて、被検体の周囲の複数のビュー(view)方向で被検体を透過したＸ線を検出系により検出し、それぞれＸ線による被検体の投影データを測定し、この投影データに基づいて、被検体の断層像を再構成するように構成されている。
【０００３】
また、被検体の断層像は、被検体の投影データに対して、所定の前処理を施し、この前処理を施した投影データを、たとえば、フィルタード・バックプロジェクション(filtered back projection)法により再構成される。
フィルタード・バックプロジェクションでは、被検体の投影データに対して画像再構成のために必要なフィルタリング処理を行ったのち、フィルタリング処理後の投影データを逆投影することにより、画像を再構成する。
上記のフィルタリング処理では、診断目的や対象部位に応じて各種のフィルタが選択される。このフィルタリング処理は、投影データを高速フーリエ変換(FFT:Fast Fourier Transfer) したのち、このデータに選択されたフィルタが乗算され、再び逆高速フーリエ変換(IFFT: Inverse FFT) される。フィルタリング処理では、投影データに含まれる高周波のノイズ成分も同時に除去される。なお、投影データに含まれる高周波のノイズ成分を除去する技術は、たとえば、特許文献１に開示されている。
【０００４】
【特許文献１】
特開２０００−８３９４６号公報
【０００５】
【発明が解決しようとする課題】
ところで、金属等の人体とは異なる異物質が被検体に存在すると、検出された投影データのうち、この異物質に相当する部分の値が急峻に変化し、この部分は他の部分とは非連続な部分となる。金属は、Ｘ線の吸収率が非常に大きいからである。
金属等の異物質の像を含む投影データに上記のフィルタリング処理を施すと、異物質に対応する部分のデータは値が急峻に変化する高周波成分をもつので、この部分のデータ値が変化し、逆高速フーリエ変換したのちの投影データにはアンダーシュートやオーバシュートと呼ばれる歪みが発生することがある。
この歪みが発生した投影データを逆投影して画像再構成すると、金属等の異物質に対応する画像にいわゆるストリーク・アーチファクト(streak artifact) が発生するという不利益が存在した。
【０００６】
本発明は、上述の問題に鑑みて成されたものであって、その目的は、金属等の放射線吸収率の高い物質が被検体に存在することにより再構成画像に発生するストリーク・アーチファクトを低減可能な放射線断層像撮影装置、画像再構成装置および画像再構成方法を提供することにある。
【０００７】
【課題を解決するための手段】
本発明の画像再構成装置は、被検体を透過した放射線を複数の放射線検出素子によって検出して得られた投影データに基づいて、前記被検体の断層像を再構成する画像再構成装置であって、前記放射線検出素子の配列順に並ぶ前記投影データのうち、非連続的に値が異なるデータを含む非連続領域を除くデータに対して所定のフィルタリング処理を施こし、前記非連続領域のデータと前記フィルタリング処理が施されたデータとを前記断層像を再構成するための投影データとする非連続データ処理手段を有する。
【０００８】
好適には、前記非連続データ処理手段は、前記非連続領域のデータを前記投影データから分離し、当該投影データを値が連続的に変化する連続データへ変換し、当該連続データのみを前記フィルタリング処理の対象にする分離手段と、前記フィルタリング処理された前記連続データと前記非連続領域のデータとを合成し、前記画像再構成手段に与える投影データとする合成手段とを有する。
【０００９】
前記非連続データ処理手段が、前記フィルタリング処理されたのちの投影データと、前記フィルタリング処理されていない投影データとの差に基づいて、前記非連続領域を抽出する抽出手段と、前記フィルタリング処理された投影データのうち、抽出された前記非連続領域のデータをフィルタリング処理前の投影データに置き換え、前記画像再構成手段に与える投影データとする合成手段と
を有する構成とすることも可能である。
【００１０】
本発明の画像再構成方法は、被検体を透過した放射線を複数の放射線検出素子によって検出して得られた投影データに基づいて、前記被検体の断層像を再構成する画像再構成方法であって、前記放射線検出素子の配列順に並ぶ前記投影データのうち、非連続的に値が異なるデータを含む非連続領域を除くデータに対して所定のフィルタリング処理を施こし、前記非連続領域のデータと前記フィルタリング処理が施されたデータとを前記断層像を再構成するための投影データとする。
【００１１】
本発明の放射線断層像撮影装置は、被検体を挟むようにして相対向して配される放射線照射手段および複数の放射線検出素子からなる放射線検出器と、前記被検体を透過した放射線を前記複数の放射線検出素子によって検出して得られた投影データに基づいて、前記被検体の断層像を再構成する画像再構成装置とを有し、前記画像再構成装置は、前記放射線検出素子の配列順に並ぶ前記投影データのうち、非連続的に値が異なるデータを含む非連続領域を除くデータに対して所定のフィルタリング処理を施こし、前記非連続領域のデータと前記フィルタリング処理が施されたデータとを前記断層像を再構成するための投影データとする非連続データ処理手段を有する。
【００１２】
本発明では、投影データから断層像を再構成する前に行うフィルタリング処理において、投影データに含まれる非連続データにはフィルタリング処理を施さずにそのまま画像再構成に用いる。これにより、非連続データのもつ高周波成分がフィルタリング処理によって歪むことがなく、この非連続部分の画像が正確に再構成される。
【００１３】
【発明の実施の形態】
以下、本発明の実施の形態について図面を参照して説明する。
第１実施形態
図１は、本発明の放射線断層像撮影装置の一実施形態としてのＸ線ＣＴ装置の全体構成を示すブロック図である。
【００１４】
図１に示すように、Ｘ線ＣＴ装置１は、走査ガントリ(gantry)２、操作コンソール(console) ３、および撮影テーブル（クレードル）４を有している。
【００１５】
走査ガントリ２は、Ｘ線管２１を有する。Ｘ線管２１から放射された図示しないＸ線は、コリメータ(collimator)２２により、たとえば、扇状のＸ線ビーム(beam)となるように成形され、Ｘ線検出器２３に照射される。
【００１６】
Ｘ線検出器２３は、ファンビームＸ線の広がり方向にアレイ状に配列された複数のＸ線検出素子を有する。Ｘ線検出器２３は、複数のＸ線検出素子をアレイ状に配列した、多チャネル(channel) の検出器となっている。
【００１７】
Ｘ線検出器２３は、全体として、円筒凹面状に湾曲したＸ線入射面を形成する。Ｘ線検出器２３は、たとえば、シンチレータ(scintillator)とフォトダイオード(photo diode) の組み合わせによって構成される。なお、これに限られず、たとえば、カドミウム・テルル（ＣｄＴｅ）等を利用した半導体Ｘ線検出素子またはＸｅガス(gas) を用いる電離箱型のＸ線検出素子であっても良い。
【００１８】
Ｘ線検出器２３にはデータ収集部２４が接続されている。データ収集部２４は、Ｘ線検出器２３の個々のＸ線検出素子の検出データを収集する。Ｘ線管２１からのＸ線の照射は、Ｘ線コントローラ(controller)２５によって制御される。なお、Ｘ線管２１とＸ線コントローラ２５との接続関係およびコリメータ２２とコリメータコントローラ２６との接続関係については図示を省略する。コリメータ２２は、コリメータコントローラ２６によって制御される。
【００１９】
Ｘ線管２１、コリメータ２２、Ｘ線検出器２３、データ収集部２４、Ｘ線コントローラ２５およびコリメータコントローラ２６が、走査ガントリ２の回転部２７に搭載されている。ここで、被検体は、回転部２７の中心に位置するボア(bore)２９内のクレードル(cradle)上に載置される。回転部２７は、回転コントローラ２８により制御されつつ回転し、Ｘ線管２１からＸ線を爆射し、Ｘ線検出器２３において被検体の透過Ｘ線を、各ビューごとの投影情報として検出する。なお、回転部２７と回転コントローラ２８との接続関係については図示を省略する。
【００２０】
操作コンソール６は、中央処理装置３１、入力装置３２、表示装置３３および記憶装置３４を有する。
なお、中央処理装置３１および記憶装置３４が、本発明の画像再構成装置を構成している。
【００２１】
中央処理装置３１は、たとえば、マイクロプロセッサ、メモリ等によって構成される。
この中央処理装置３１は、記憶装置３４に記憶されたソフトウエアにしたがって、走査ガントリ２の動作を制御し、被検体を透過したＸ線をＸ線検出器２３で検出して得られる投影データを収集する機能と、収集したＸ線投影データに基づいて、被検体の断層像を再構成する機能とを少なくとも有する。
【００２２】
中央処理装置３１には、データ収集部２４で収集されたデータが図示しないデータ収集バッファを介して入力され、中央処理装置３１は、収集した投影データを用いて画像再構成を行う。この中央処理装置３１による画像再構成処理については後述する。
【００２３】
中央処理装置３１には、また、表示装置３３と入力装置３２がそれぞれ接続されている。表示装置３３は、中央処理装置３１から出力される断層画像情報やその他の情報を表示する。入力装置３２は、オペレータによって操作され、各種の指示や情報等を中央処理装置３１に入力する。オペレータは表示装置３３および入力装置３２を使用してインタラクティブ(interactive) に本装置を操作する。
【００２４】
上記構成のＸ線ＣＴ装置１では、走査ガントリ２において、Ｘ線管２１およびＸ線検出器２３とを相対的に移動させ、被検体の被検対象領域を走査することにり、データ収集部２４に被検体の投影データが収集され、この投影データが中央処理装置３１へ送られ、中央処理装置３１において被検体の断層画像が再構成される。
【００２５】
図２は、中央処理装置３１による画像再構成処理手順の概略を示す図である。中央処理装置３１における画像再構成は、フィルタード・バックプロジェクション(filtered back projection)法によって行われる。なお、フィルタード・バックプロジェクション法については周知の方法であるので詳細な説明については省略する。
【００２６】
図２に示すように、中央処理装置３１による画像再構成処理は、大きく分けて、前処理（ステップＳＴ１）と、分離処理（ステップＳＴ２）、フィルタリング処理（ステップＳＴ３）および合成処理（ステップＳＴ４）からなる非連続データ処理と、バックプロジェクション処理（ステップＳＴ５）とからなる。
【００２７】
前処理（ステップＳＴ１）では、正常な画像が再構成されるように、データ収集部２４から得られた被検体の投影データに、物理的および電気的な補正処理を施す。
【００２８】
ここで、図３は分離処理（ステップＳＴ２）の内容を説明するための図であり、（ａ）は、データ収集部２４から得られた被検体の投影データの一例を示しており、（ｂ）は分離処理後の投影データを示しており、（ｃ）は分離された非連続データを示している。なお、これらのＸ線検出器２３の各検出器アレイのチャンネル順（配列順）に並んでいる。
非連続データとは、他の連続的に変化するデータとは複数チャンネルにわたって値が異なるデータを含む部分的なデータである。
【００２９】
たとえば、金属等の人体とは異なりＸ線の吸収率が高い異物質が被検体に存在すると、図３（ａ）に示すように、投影データｄａｔ＿Ｒ中において、異物質に対応する部分のデータ値が急峻に変化する。投影データｄａｔ＿Ｒの異物質が存在していない部分は、連続的に滑らかにデータ値が変化する。
【００３０】
図３（ａ）に示す投影データｄａｔ＿Ｒ中の、値が急峻に変化するデータを含む非連続領域のデータを非連続データｄａｔ＿ＥＸとすると、この非連続データｄａｔ＿ＥＸは、データ値を周波数領域でみたときに高周波成分を含んでいる。非連続データｄａｔ＿ＥＸ以外の部分には、データ値が滑らかに変化するため、このような高周波成分は存在していない。
【００３１】
投影データｄａｔ＿Ｒに含まれる高周波のノイズ成分を除去するために、投影データｄａｔ＿Ｒの全てにフィルタリングを施すと、上記の非連続データｄａｔ＿ＥＸに含まれる高周波成分も除去される可能性がある。
このため、分離処理（ステップＳＴ２）においては、まず、図３（ｃ）に示すように、投影データｄａｔ＿Ｒから非連続データｄａｔ＿ＥＸを分離する。
【００３２】
また、非連続データｄａｔ＿ＥＸを分離したのちの投影データｄａｔ＿Ｒは、図３（ｂ）に示すように、非連続データｄａｔ＿ＥＸが存在していた部分をデータｄａｔ＿Ｃにより補完する。具体的には、非連続データｄａｔ＿ＥＸの周辺のデータを用いてフィッティングを行う。
これにより、すべての領域でデータ値が滑らかに変化する連続データｄａｔ＿Ｓが得られる。
この連続データｄａｔ＿Ｓは、次にフィルタリング処理される。
【００３３】
フィルタリング処理（ステップＳＴ３）では、連続データｄａｔ＿Ｓを、まず、高速フーリエ変換（ＦＦＴ）する。
次いで、診断目的や対象部位に応じて選択されたフィルタによってフィルタリング処理される。
次いで、フィルタリング処理された連続データｄａｔ＿Ｓを逆高速フーリエ変換（ＩＦＦＴ）して復元する。
【００３４】
フィルタリング処理（ステップＳＴ３）におけるフィルタは、診断目的や対象部位に応じて適宜選択されるが、連続データｄａｔ＿Ｓに含まれる高周波のノイズ成分も同時に除去する。
ここで、図４は、フィルタリング処理（ステップＳＴ３）におけるフィルタの特性を説明するためのグラフであり、（ａ）はＦＦＴ前のデータを示しており、（ｂ）はＩＦＦＴ後のデータを示している。
【００３５】
図４（ａ）からわかるように、ＦＦＴ前のデータにはノイズ成分が存在しているが、フィルタリング処理によってこれらのノイズは除去される。
また、データ中に値が急峻に変化する部分が存在すると、ＩＦＦＴ後のデータには、図４（ｂ）の円Ａ内に示すように、この部分周辺においてオーバシュートやアンダーシュートと呼ばれる歪みが発生する。
【００３６】
本実施形態では、上記した分離処理（ステップＳＴ２）において、高周波成分をもつ非連続データｄａｔ＿ＥＸを分離してフィルタリング処理の処理対象から除外している。このため、連続データｄａｔ＿Ｓは、フィルタリング処理によってノイズ成分だけが除去され、歪みが発生することがない。
【００３７】
合成処理（ステップＳＴ４）では、分離された非連続データｄａｔ＿ＥＸと、フィルタリング処理された連続データｄａｔ＿Ｓとを合成し、この合成されたデータをバックプロジェクション処理（ステップＳＴ５）の画像再構成に用いるデータとする。
【００３８】
以上のように、本実施形態では、投影データｄａｔ＿Ｒに含まれる非連続データｄａｔ＿ＥＸを分離して得られる連続データｄａｔ＿Ｓのみにフィルタリング処理を施す。そして、フィルタリング処理を施した連続データｄａｔ＿Ｓと非連続データｄａｔ＿ＥＸとを合成し、この合成したデータをバックプロジェクション処理に用いる。これにより、再構成された画像に、フィルタリング処理による高周波成分の除去に起因するストリーク・アーチファクトが発生しにくくなり、金属等のＸ線の吸収率の高い異物質の断層像を正確に再現することが可能となる。
【００３９】
第２実施形態
図５は、本発明の放射線断層像撮影装置の他の実施形態に係る画像再構成処理手順の概略を示す図である。
なお、本実施形態は、上述した中央処理装置３１における画像再構成処理手順以外は第１の実施形態の構成と同じである。
【００４０】
図５に示すように、本実施形態に係る中央処理装置３１による画像再構成処理は、大きく分けて、前処理（ステップＳＴ１１）と、フィルタリング処理（ステップＳＴ１２）、差分処理（ステップＳＴ１３）、抽出処理（ステップＳＴ１４）および合成処理（ステップＳＴ１５）からなる非連続データ処理と、バックプロジェクション処理（ステップＳＴ１６）とからなる。
【００４１】
前処理（ステップＳＴ１１）は、上述した第１の実施形態と同様の処理を行う。
【００４２】
フィルタリング処理（ステップＳＴ１２）は、上述した第１の実施形態と同様のフィルタリング処理を行う。すなわち、データ収集部２４から得られた被検体の投影データに、まず、高速フーリエ変換（ＦＦＴ）し、診断目的や対象部位に応じて選択されたフィルタによってフィルタリング処理し、フィルタリング処理されたデータを逆高速フーリエ変換（ＩＦＦＴ）して復元する。
【００４３】
ここで、図６（ａ）は、フィルタリング処理前の投影データの一例であり、（ｂ）はフィルタリング処理後の投影データの一例であり、（ｃ）はフィルタリング処理前後の投影データの差分値を示している。
図６（ａ）に示すように、フィルタリング処理前の投影データｄａｔ＿Ｂにはノイズが存在するとともに、非連続的に値が異なるデータを含む非連続領域Ｒａが存在している。この非連続領域Ｒａは、上述したように、金属等のＸ線吸収率の高い異物質の存在により発生する。
【００４４】
図６（ａ）に示すような投影データｄａｔ＿Ｂに対してフィルタリング処理を施すと、たとえば、図６（ｂ）に示すデータｄａｔ＿Ａとなる。
図６（ｂ）に示すフィルタリング後のデータｄａｔ＿Ａには、非連続領域Ｒａにおけるデータ値が急峻に変化する位置で歪みｓｔが発生しているのが分かる。
【００４５】
次いで、差分処理（ステップＳＴ１３）では、フィルタリング処理前の投影データｄａｔ＿Ｂとフィルタリング処理後の投影データｄａｔ＿Ａとの差分を算出する。この差分値データは、図６（ｃ）に示すように、投影データｄａｔ＿Ｂに含まれるノイズと、フィルタリングによって発生した歪みｓｔを含んでいる。
【００４６】
次いで、抽出処理（ステップＳＴ１４）では、図６（ｃ）に示したような差分値データから、非連続領域Ｒｂ，Ｒｃを抽出する。この非連続領域Ｒｂ，Ｒｃは、投影データｄａｔ＿Ｂのうち、フィルタリングによって歪みｓｔが発生する領域である。
この非連続領域Ｒｂ，Ｒｃは、たとえば、所定のしきい値を設けて、差分値データとこのしきい値とを比較することにより抽出することができる。なお、差分値データとしきい値とを比較する際に、ノイズ成分と歪みｓｔとを区別する必要がある。歪みｓｔは複数のチャンネルにわたって存在し、ノイズ成分は単一のチャンネルに存在する場合が多く、また、歪みｓｔとノイズ成分のレベルは異なるので、両者の区別は容易である。
【００４７】
合成処理（ＳＴ１５）は、フィルタリング処理後の投影データｄａｔ＿Ａのうち、抽出処理（ステップＳＴ１４）によって抽出された非連続領域Ｒｂ，Ｒｃにおけるデータをフィルタリング処理前の投影データｄａｔ＿Ｂに置き換える。すなわち、非連続領域Ｒｂ，Ｒｃに歪みｓｔの発生していないデータを用いる。この合成されたデータをバックプロジェクション処理（ステップＳＴ１６）の画像再構成に用いるデータとする。
【００４８】
以上のように、本実施形態では、投影データｄａｔ＿Ｂに実際にフィルタリング処理を行って歪みｓｔを発生させることにより、非連続領域Ｒｂ，Ｒｃを特定し、この非連続領域Ｒｂ，Ｒｃにフィルタリング処理を行っていないデータを用いて画像再構成に用いるデータとする。これにより、バックプロジェクション処理によって再構成された画像にフィルタリング処理に起因するストリーク・アーチファクトが発生するのを確実に防ぐことができる。
【００４９】
なお以上に説明した各実施形態においては、放射線としてＸ線を用いた例について説明したが、放射線はＸ線に限るものではなく、たとえばγ線等の他の種類の放射線であっても良い。ただし、現時点では、Ｘ線がその発生、検出および制御等に関し実用的な手段が最も充実している点で好ましい。
【００５０】
【発明の効果】
本発明によれば、金属等の放射線吸収率の高い物質が被検体に存在することにより再構成画像に発生するストリーク・アーチファクトを低減可能な放射線断層像撮影装置、画像再構成装置および画像再構成方法が提供される。
【図面の簡単な説明】
【図１】本発明が適用される放射線断層像撮影装置としてのＸ線ＣＴ装置の全体構成を示すブロック図である。
【図２】本発明の一実施形態に係る中央処理装置の画像再構成処理手順を示す図である。
【図３】分離処理の内容を説明するための図である。
【図４】フィルターの特性を説明するためのグラフである。
【図５】本発明の画像再構成装置の他の実施形態に係る中央処理装置の処理手順を示す図である。
【図６】フィルタリング処理前、フィルタリング処理後および差分処理後の投影データの一例を示すグラフである。
【符号の説明】
１…Ｘ線ＣＴ装置
２…走査ガントリ
３…操作コンソール
４…撮影テーブル
２１…Ｘ線管
２２…コリメータ
２３…Ｘ線検出器
２４…データ収集部
２５…Ｘ線コントローラ
２６…コリメータコントローラ
３１…中央処理装置
３２…入力装置
３３…表示装置
３４…記憶装置
Claims (13)
- 被検体を透過した放射線を複数の放射線検出素子によって検出して得られた投影データに基づいて、前記被検体の断層像を再構成する画像再構成装置であって、
前記放射線検出素子の配列順に並ぶ前記投影データのうち、非連続的に値が異なるデータを含む非連続領域のフィルタリング処理が施されていないデータと、前記非連続領域を除く領域の所定のフィルタリング処理が施されたデータとを前記断層像を再構成するためのデータとする非連続データ処理手段であって、前記非連続領域のデータを分離し、前記非連続領域が存在した領域を補間した連続データに対しフィルタリング処理を行う分離手段と、前記フィルタリング処理された連続データとフィルタリング処理が施されていない前記非連続領域のデータとを合成する合成手段とを有する非連続データ処理手段
を有することを特徴とする画像再構成装置。 - 被検体を透過した放射線を複数の放射線検出素子によって検出して得られた投影データに基づいて、前記被検体の断層像を再構成する画像再構成装置であって、
前記放射線検出素子の配列順に並ぶ前記投影データのうち、非連続的に値が異なるデータを含む非連続領域のフィルタリング処理が施されていないデータと、前記非連続領域を除く領域の所定のフィルタリング処理が施されたデータとを前記断層像を再構成するためのデータとする非連続データ処理手段であって、前記非連続領域及び前記非連続領域を除く領域を含む領域におけるフィルタリング処理されたデータと、前記非連続領域及び前記非連続領域を除く領域を含む領域におけるフィルタリング処理されていないデータとの差に基づいて、前記非連続領域を抽出する抽出手段と、前記非連続領域及び前記非連続領域を除く領域を含む領域のフィルタリング処理されたデータのうち、抽出された前記非連続領域のデータをフィルタリング処理されていない前記非連続領域のデータに置き換える合成手段とを有する非連続データ処理手段
を有することを特徴とする画像再構成装置。 - 前記放射線は、Ｘ線であることを特徴とする請求項１または２に記載の画像再構成装置。
- 被検体を透過した放射線を複数の放射線検出素子によって検出して得られた投影データに基づいて、前記被検体の断層像を再構成する画像再構成方法であって、
前記放射線検出素子の配列順に並ぶ前記投影データのうち、非連続的に値が異なるデータを含む非連続領域のフィルタリング処理が施されていないデータと、前記非連続領域を除く領域の所定のフィルタリング処理が施されたデータとを前記断層像を再構成するためのデータとする非連続データ処理工程であって、前記非連続領域のデータを分離し、前記非連続領域が存在した領域を補間した連続データに対しフィルタリング処理を行う工程と、前記フィルタリング処理された連続データとフィルタリング処理が施されていない前記非連続領域のデータとを合成する工程とを含む非連続データ処理工程
を含むことを特徴とする画像再構成方法。 - 被検体を透過した放射線を複数の放射線検出素子によって検出して得られた投影データに基づいて、前記被検体の断層像を再構成する画像再構成方法であって、
前記放射線検出素子の配列順に並ぶ前記投影データのうち、非連続的に値が異なるデータを含む非連続領域のフィルタリング処理が施されていないデータと、前記非連続領域を除く領域の所定のフィルタリング処理が施されたデータとを前記断層像を再構成するためのデータとする非連続データ処理工程であって、前記非連続領域及び前記非連続領域を除く領域を含む領域におけるフィルタリング処理されたデータと、前記非連続領域及び前記非連続領域を除く領域を含む領域におけるフィルタリング処理されていないデータとの差に基づいて、前記非連続領域を抽出する工程と、前記非連続領域及び前記非連続領域を除 く領域を含む領域のフィルタリング処理されたデータのうち、抽出された前記非連続領域のデータをフィルタリング処理されていない前記非連続領域のデータに置き換える工程とを含む非連続データ処理工程
を含むことを特徴とする画像再構成方法。 - 前記放射線に、Ｘ線を用いることを特徴とする請求項４または５に記載の画像再構成方法。
- 被検体を挟むようにして相対向して配置される放射線照射手段および複数の放射線検出素子からなる放射線検出器と、
前記被検体を透過した放射線を前記複数の放射線検出素子によって検出して得られた投影データに基づいて、前記被検体の断層像を再構成する画像再構成装置とを有し、
前記画像再構成装置は、前記放射線検出素子の配列順に並ぶ前記投影データのうち、非連続的に値が異なるデータを含む非連続領域のフィルタリング処理が施されていないデータと、前記非連続領域を除く領域の所定の前記フィルタリング処理が施されたデータとを前記断層像を再構成するためのデータとする非連続データ処理手段であって、前記非連続領域のデータを分離し、前記非連続領域が存在した領域を補間した連続データに対しフィルタリング処理を行う分離手段と、前記フィルタリング処理された連続データとフィルタリング処理が施されていない前記非連続領域のデータとを合成する合成手段とを有する非連続データ処理手段
を有することを特徴とする放射線断層像撮影装置。 - 被検体を挟むようにして相対向して配置される放射線照射手段および複数の放射線検出素子からなる放射線検出器と、
前記被検体を透過した放射線を前記複数の放射線検出素子によって検出して得られた投影データに基づいて、前記被検体の断層像を再構成する画像再構成装置とを有し、
前記画像再構成装置は、前記放射線検出素子の配列順に並ぶ前記投影データのうち、非連続的に値が異なるデータを含む非連続領域のフィルタリング処理が施されていないデータと、前記非連続領域を除く領域の所定の前記フィルタリング処理が施されたデータとを前記断層像を再構成するためのデータとする非連続データ処理手段であって、前記非連続領域及び前記非連続領域を除く領域を含む領域におけるフィルタリング処理されたデータと、前記非連続領域及び前記非連続領域を除く領域を含む領域におけるフィルタリング処理されていないデータとの差に基づいて、前記非連続領域を抽出する抽出手段と、前記非連続領域及び前記非連続領域を除く領域を含む領域のフィルタリング処理されたデータのうち、抽出された前記非連続領域のデータをフィルタリング処理されていない前記非連続領域のデータに置き換える合成手段とを有する非連続データ処理手段
を有することを特徴とする放射線断層像撮影装置。 - 前記放射線照射手段は、Ｘ線管を有することを特徴とする請求項７または８に記載の放射線断層像撮影装置。
- 前記Ｘ線管から放射されたＸ線を、扇状のＸ線ビームに成形するコリメータをさらに有することを特徴とする請求項９に記載の放射線断層像撮影装置。
- 前記放射線検出器は、複数の放射線検出素子をアレイ状に配列して構成されることを特徴とする請求項７〜１０の何れかに記載の放射線断層像撮影装置。
- 前記画像再構成装置は、フィルタード・バックプロジェクション（ｆｉｌｔｅｒｅｄ ｂａｃｋ ｐｒｏｊｅｃｔｉｏｎ）法を用いることを特徴とする請求項７〜１１のいずれかに記載の放射線断層像撮影装置。
- 前記放射線検出器は、シンチレータ（ｓｃｉｎｔｉｌｌａｔｏｒ）とフォトダイオード（ｐｈｏｔｏ ｄｉｏｄｅ） の組み合わせによって構成されることを特徴とする請求項１１に記載の放射線断層像撮影装置。
Priority Applications (6)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003100668A JP4138558B2 (ja) | 2003-04-03 | 2003-04-03 | 画像再構成装置、画像再構成方法および放射線断層像撮影装置 |
| DE200460000481 DE602004000481T2 (de) | 2003-04-03 | 2004-03-30 | Verfahren und Vorrichtung zur Bildrekonstruktion und Vorrichtung zur Röntgenstrahlbildgebung |
| EP20040251882 EP1475750B1 (en) | 2003-04-03 | 2004-03-30 | Method and apparatus for image reconstruction and X-Ray CT imaging apparatus |
| KR1020040022857A KR20040086814A (ko) | 2003-04-03 | 2004-04-02 | 화상 재구성 방법 및 장치와 방사선 토모그라피 촬상 장치 |
| US10/817,222 US6980622B2 (en) | 2003-04-03 | 2004-04-02 | Method and apparatus for image reconstruction and X-ray CT imaging apparatus |
| CNB2004100342908A CN1299647C (zh) | 2003-04-03 | 2004-04-05 | 图像重构的方法与装置和x射线计算层析照像成像装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003100668A JP4138558B2 (ja) | 2003-04-03 | 2003-04-03 | 画像再構成装置、画像再構成方法および放射線断層像撮影装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2004305331A JP2004305331A (ja) | 2004-11-04 |
| JP4138558B2 true JP4138558B2 (ja) | 2008-08-27 |
Family
ID=32985526
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2003100668A Expired - Fee Related JP4138558B2 (ja) | 2003-04-03 | 2003-04-03 | 画像再構成装置、画像再構成方法および放射線断層像撮影装置 |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US6980622B2 (ja) |
| EP (1) | EP1475750B1 (ja) |
| JP (1) | JP4138558B2 (ja) |
| KR (1) | KR20040086814A (ja) |
| CN (1) | CN1299647C (ja) |
| DE (1) | DE602004000481T2 (ja) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006167161A (ja) * | 2004-12-16 | 2006-06-29 | Tokyo Institute Of Technology | Ｘ線ｃｔ画像再構成方法、装置及びｘ線ｃｔ画像再構成プログラムを記憶した記憶媒体 |
| JP4268976B2 (ja) * | 2006-06-15 | 2009-05-27 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | イメージング装置 |
| CN101495038A (zh) | 2006-08-03 | 2009-07-29 | 加州大学董事会 | 层析中剂量减少和图像增强的迭代方法 |
| KR100923098B1 (ko) * | 2008-01-15 | 2009-10-22 | (주)이우테크놀로지 | 엑스 레이 씨티 촬영 영상의 메탈 아티팩트를 제거하는방법 |
| JP5601675B2 (ja) * | 2008-02-29 | 2014-10-08 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ｘ線ｃｔ装置およびプログラム |
| JP5324880B2 (ja) * | 2008-10-21 | 2013-10-23 | アズビル株式会社 | Ｃｔ装置および金属形状抽出方法 |
| JP5324883B2 (ja) * | 2008-10-24 | 2013-10-23 | アズビル株式会社 | Ｃｔ装置および金属形状抽出方法 |
| WO2012164921A1 (ja) * | 2011-05-31 | 2012-12-06 | 株式会社島津製作所 | 放射線断層画像生成方法および放射線断層画像生成プログラム |
| JP6122269B2 (ja) * | 2011-12-16 | 2017-04-26 | キヤノン株式会社 | 画像処理装置、画像処理方法、及びプログラム |
| CN105374012B (zh) * | 2014-08-27 | 2018-11-27 | 通用电气公司 | 用于消除由性能差异的探测器单元所致的条状伪影的方法 |
| CN106994021B (zh) * | 2016-01-22 | 2022-10-14 | 通用电气公司 | 一种计算ct影像上的噪声的方法及装置 |
Family Cites Families (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4590558A (en) * | 1981-12-30 | 1986-05-20 | General Electric Company | Method and apparatus for removing objects from CT images |
| US4920491A (en) * | 1988-05-16 | 1990-04-24 | General Electric Company | Enhancement of image quality by utilization of a priori information |
| GB8900165D0 (en) * | 1989-01-05 | 1989-03-01 | Glaverbel | Glass coating |
| US5243664A (en) * | 1991-09-16 | 1993-09-07 | Picker International, Inc. | Post-processing technique for reducing metallic clip artifacts in CT images |
| JP3673041B2 (ja) * | 1996-11-21 | 2005-07-20 | ジーイー横河メディカルシステム株式会社 | Ｃｔイメージング方法およびｘ線ｃｔ装置 |
| US5933471A (en) * | 1997-02-26 | 1999-08-03 | International Business Machines Corporation | System and method for reducing reconstruction artifacts in computed tomography images |
| IL121773A0 (en) * | 1997-09-15 | 1998-02-22 | Elscint Ltd | Method for improving CT images having high attenuation objects |
| US6201888B1 (en) * | 1998-02-18 | 2001-03-13 | International Business Machines Corporation | System and method for restoring, describing and graphically displaying noise-corrupted boundaries in tomography images |
| US6108575A (en) * | 1998-02-20 | 2000-08-22 | General Electric Company | Helical weighting algorithms for fast reconstruction |
| US6125193A (en) * | 1998-06-01 | 2000-09-26 | Kabushiki Kaisha Toshiba | Method and system for high absorption object artifacts reduction and superposition |
| JP2000083946A (ja) | 1998-09-14 | 2000-03-28 | Ge Yokogawa Medical Systems Ltd | プロジェクション補正方法および装置並びに放射線断層撮影装置 |
| US6351514B1 (en) * | 2000-06-22 | 2002-02-26 | Ge Medical Systems Global Technology Company, Llc | Slice-adaptive multislice helical weighting for computed tomography imaging |
| DE10123797B4 (de) * | 2001-05-16 | 2006-10-19 | Siemens Ag | Verfahren zur Herstellung eines Bildes mit einem Computertomographen |
| JP3963086B2 (ja) * | 2001-07-13 | 2007-08-22 | 株式会社島津製作所 | 断層再構成ソフトウエアとその記録媒体と断層撮影装置 |
-
2003
- 2003-04-03 JP JP2003100668A patent/JP4138558B2/ja not_active Expired - Fee Related
-
2004
- 2004-03-30 DE DE200460000481 patent/DE602004000481T2/de not_active Expired - Lifetime
- 2004-03-30 EP EP20040251882 patent/EP1475750B1/en not_active Expired - Lifetime
- 2004-04-02 KR KR1020040022857A patent/KR20040086814A/ko not_active Application Discontinuation
- 2004-04-02 US US10/817,222 patent/US6980622B2/en not_active Expired - Fee Related
- 2004-04-05 CN CNB2004100342908A patent/CN1299647C/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| EP1475750A2 (en) | 2004-11-10 |
| US6980622B2 (en) | 2005-12-27 |
| US20040234022A1 (en) | 2004-11-25 |
| CN1299647C (zh) | 2007-02-14 |
| CN1536533A (zh) | 2004-10-13 |
| JP2004305331A (ja) | 2004-11-04 |
| EP1475750A3 (en) | 2004-12-15 |
| KR20040086814A (ko) | 2004-10-12 |
| DE602004000481D1 (de) | 2006-05-11 |
| EP1475750B1 (en) | 2006-03-15 |
| DE602004000481T2 (de) | 2006-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5180181B2 (ja) | コンピュータ断層撮影データ収集装置及び方法 | |
| JP4575909B2 (ja) | Ｘ線断層撮影装置 | |
| EP2446821B1 (en) | Dynamic collimator for wide coverage and low dose cardiac CT imaging | |
| JP2004174249A (ja) | ボリューム灌流を計算するための方法及び装置 | |
| KR20230153347A (ko) | 좁은 시야 ｘ선 이미징 시스템 및 방법 | |
| JP4138558B2 (ja) | 画像再構成装置、画像再構成方法および放射線断層像撮影装置 | |
| JP2004174264A (ja) | 計算機式断層写真法（ｃｔ）スカウト画像を形成する方法及び装置 | |
| JP3987024B2 (ja) | 横方向のフィルタリング処理を用いたトモシンセシス画像を強調する方法及びシステム | |
| US20200226800A1 (en) | Tomographic imaging apparatus and method of generating tomographic image | |
| CN102074276A (zh) | 用于过滤x射线的过滤器和x射线计算机断层成像仪 | |
| JP5097355B2 (ja) | 放射線断層撮影装置 | |
| EP1433137B1 (en) | A method of reducing artifacts in object images | |
| JP5144723B2 (ja) | Ｘ線断層撮影装置 | |
| JP4090970B2 (ja) | 放射線断層画像撮影装置と放射線断層画像撮影方法および画像生成装置と画像生成方法 | |
| KR20180063753A (ko) | 의료 영상 장치 및 동작 방법 | |
| US6647084B1 (en) | Method and apparatus for filtering projection data of a helical scan | |
| JP4559723B2 (ja) | 放射線ｃｔ装置、画像処理装置、及び画像処理方法 | |
| JP5596100B2 (ja) | Ｘ線断層撮影装置 | |
| JPH1075947A (ja) | 画像再構成処理装置におけるアーチファクト低減方法 | |
| JP3678375B2 (ja) | 放射線断層撮影装置 | |
| JP4282110B2 (ja) | Ｘ線コンピュータ断層撮影装置 | |
| JP2000083946A (ja) | プロジェクション補正方法および装置並びに放射線断層撮影装置 | |
| JP2009106759A (ja) | コンピュータ断層撮影装置、その処理方法及び記録媒体 | |
| US8005183B2 (en) | Computed tomography device with active adaptation of the measuring electronics | |
| JP2000229076A (ja) | Ｘ線コンピュータ断層撮影装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A625 | Written request for application examination (by other person) |
Free format text: JAPANESE INTERMEDIATE CODE: A625 Effective date: 20041020 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20060220 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060228 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060530 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20061121 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20070219 |
|
| A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20070222 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20071113 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080213 |
|
| 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: 20080527 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080605 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110613 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110613 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110613 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120613 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120613 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130613 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130613 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130613 Year of fee payment: 5 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |