JP5243160B2 - X-ray CT system - Google Patents
X-ray CT system Download PDFInfo
- Publication number
- JP5243160B2 JP5243160B2 JP2008235985A JP2008235985A JP5243160B2 JP 5243160 B2 JP5243160 B2 JP 5243160B2 JP 2008235985 A JP2008235985 A JP 2008235985A JP 2008235985 A JP2008235985 A JP 2008235985A JP 5243160 B2 JP5243160 B2 JP 5243160B2
- Authority
- JP
- Japan
- Prior art keywords
- ray
- reprojection
- projection data
- component
- image
- 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
- 238000012937 correction Methods 0.000 claims description 81
- 239000006096 absorbing agent Substances 0.000 claims description 69
- 238000004364 calculation method Methods 0.000 claims description 35
- 239000000126 substance Substances 0.000 claims description 29
- 210000004872 soft tissue Anatomy 0.000 claims description 19
- 230000002745 absorbent Effects 0.000 claims description 2
- 239000002250 absorbent Substances 0.000 claims description 2
- 238000005510 radiation hardening Methods 0.000 claims 1
- 238000003384 imaging method Methods 0.000 description 44
- 239000002184 metal Substances 0.000 description 44
- 229910052751 metal Inorganic materials 0.000 description 44
- 210000000988 bone and bone Anatomy 0.000 description 32
- 239000000463 material Substances 0.000 description 20
- 238000010521 absorption reaction Methods 0.000 description 16
- 238000000034 method Methods 0.000 description 14
- 238000012545 processing Methods 0.000 description 14
- 238000001914 filtration Methods 0.000 description 12
- 230000009466 transformation Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 230000009467 reduction Effects 0.000 description 7
- 230000005540 biological transmission Effects 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 6
- 238000001514 detection method Methods 0.000 description 6
- 230000003287 optical effect Effects 0.000 description 6
- 210000000588 acetabulum Anatomy 0.000 description 4
- 210000000527 greater trochanter Anatomy 0.000 description 4
- 210000004394 hip joint Anatomy 0.000 description 4
- 238000011540 hip replacement Methods 0.000 description 4
- 238000007781 pre-processing Methods 0.000 description 4
- 239000000919 ceramic Substances 0.000 description 3
- 239000002872 contrast media Substances 0.000 description 3
- 210000004513 dentition Anatomy 0.000 description 3
- 239000000835 fiber Substances 0.000 description 3
- 230000000399 orthopedic effect Effects 0.000 description 3
- 238000011002 quantification Methods 0.000 description 3
- 230000036346 tooth eruption Effects 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 150000002739 metals Chemical class 0.000 description 2
- 238000004904 shortening Methods 0.000 description 2
- 210000001015 abdomen Anatomy 0.000 description 1
- 238000012884 algebraic function Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000012888 cubic function Methods 0.000 description 1
- 210000003298 dental enamel Anatomy 0.000 description 1
- 210000003414 extremity Anatomy 0.000 description 1
- 210000001624 hip Anatomy 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 210000004705 lumbosacral region Anatomy 0.000 description 1
- 210000004197 pelvis Anatomy 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 210000000689 upper leg Anatomy 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Description
本発明はX線CT装置に係り、特に、X線の高吸収体成分から生じるアーチファクトを補正するために有効な技術に関する。 The present invention relates to an X-ray CT apparatus, and more particularly to a technique effective for correcting an artifact caused by a high-absorber component of X-rays.
X線CT装置を用いて被検体のX線CT像を生成する場合、該被検体に金属等の高吸収体成分が存在すると、所謂高吸収体成分アーチファクトが発生することが知られている。 When an X-ray CT image of a subject is generated using an X-ray CT apparatus, it is known that a so-called high-absorber component artifact occurs when a high-absorber component such as a metal exists in the subject.
このため、一度再構成CT像を生成し、前記金属等を含む高吸収体成分領域を抽出し再投影データ上で非線形処理を施し、該高吸収体成分アーチファクトを低減させる手法が、例えば、下記特許文献1、あるいは非特許文献1において開示されている。
なお、2次元検出器の回転撮影データから3次元的CT像を生成するための再構成演算に、下記非特許文献2に記載のフェルドカンプの方法に基づいた再構成アルゴリズムが広く用いられている。
Note that a reconstruction algorithm based on the Feldkamp method described in Non-Patent Document 2 below is widely used for the reconstruction calculation for generating a three-dimensional CT image from the rotational imaging data of the two-dimensional detector. .
しかし、2つ以上の異なるX線吸収係数の高吸収体成分を含む被検体を再構成しようとする場合、上述したX線CT装置では、1つの再投影閾値を設定して再投影データを生成し補正演算を行うようになっていた。そのため、X線吸収係数の高い方の閾値で再投影データを生成し補正演算を行うと、図6−2に示すようにアーチファクトの低減が不十分なものとなり、低い再投影閾値で補正演算を行うと、図6−3のようにX線吸収係数の高い方の2つの物体を結ぶ方向のアーチファクトが信号として認識され、逆に偽像として表示されるものとなっていた。またこのような場合に、まず第1段階において高い再投影閾値を設定しX線吸収係数の高い方の吸収体から発生するアーチファクトを消去し、次に第2段階において低い再投影閾値を設定する事により2つのX線吸収係数の異なる高吸収体成分から生じるアーチファクトを低減するというように、2段階のステップを踏む手法も考えられるが、このような手法は再投影、再構成演算を再投影閾値の数相当の回数分繰り返す事になり、演算時間を要するものとなるという問題がある。 However, when reconstructing a subject that includes a high-absorber component having two or more different X-ray absorption coefficients, the above-described X-ray CT apparatus sets one reprojection threshold and generates reprojection data. Then, the correction calculation was performed. Therefore, if reprojection data is generated with a threshold having a higher X-ray absorption coefficient and correction calculation is performed, artifact reduction is insufficient as shown in FIG. 6-2, and correction calculation is performed with a low reprojection threshold. As a result, as shown in FIG. 6-3, the artifact in the direction connecting the two objects having the higher X-ray absorption coefficient is recognized as a signal and is displayed as a false image. In such a case, first, a high reprojection threshold is set in the first stage to eliminate artifacts generated from the absorber having the higher X-ray absorption coefficient, and then a low reprojection threshold is set in the second stage. A technique that takes two steps, such as reducing artifacts caused by two high-absorber components with different X-ray absorption coefficients, can be considered, but this technique reprojects the reprojection and reconstruction operations. There is a problem in that it is repeated for the number of times corresponding to the number of thresholds, and calculation time is required.
本発明の目的は、高吸収体成分アーチファクトの低減に要する演算時間を短くし、かつ高精度な高吸収体補正を行い得るX線CT装置を提供することにある。 An object of the present invention is to provide an X-ray CT apparatus capable of shortening the calculation time required for reducing the high-absorber component artifact and performing high-accuracy high-absorber correction.
本発明は、複数の再投影係数を有する再投影データを一度に生成し、被検体の物質の成分ごとに、成分ごとに異なる非線形変換を施した補正投影データを生成する事により、高吸収体補正の高度化と、演算時間の増大を抑えたX線CT装置を実現する。 The present invention generates reprojection data having a plurality of reprojection coefficients at a time, and generates corrected projection data that is subjected to different nonlinear transformation for each component of the substance of the subject, thereby obtaining a high absorber. An X-ray CT apparatus that realizes advanced correction and an increase in calculation time is realized.
より詳しくは、本発明に係るX線CT装置は、X線を発生させるX線源と、前記X線源に対向して配置され、高吸収体成分を含む被検体を透過した前記X線を検出して前記被検体の投影データを出力するX線検出器と、前記X線源と前記X線検出器とを回転移動させる回転手段と、前記被検体の投影データに基づいて再構成演算を行い、前記被検体のX線CT像を生成する再構成手段と、前記再構成手段が生成した画像を表示する表示手段と、前記X線CT像を前記高吸収体成分に対応させた再投影係数に基づいて再投影演算をすることにより前記高吸収体成分の再投影データを生成し、その高吸収体成分の再投影データと前記被検体の投影データとに基づいて、前記高吸収体成分の投影データと前記被検体に含まれる前記高吸収体成分を除く物質成分の投影データとを生成し、それら前記高吸収体成分の投影データと前記物質成分の投影データとの各々に対して前記X線の線質硬化に起因するアーチファクトを低減するための補正を施す高吸収体補正手段と、を備え、前記再構成手段は、前記補正後の前記高吸収体成分の投影データ及び前記物質成分の投影データに基づいて再構成演算を行い、前記アーチファクトを補正したX線CT像を生成し、前記表示手段は、前記アーチファクトを補正したX線CT像を表示する、ことを特徴とする。 More specifically, an X-ray CT apparatus according to the present invention includes an X-ray source that generates X-rays, and the X-ray that is disposed to face the X-ray source and that has passed through a subject containing a high-absorber component. An X-ray detector that detects and outputs projection data of the subject, a rotating means that rotationally moves the X-ray source and the X-ray detector, and a reconstruction calculation based on the projection data of the subject Reconstructing means for generating an X-ray CT image of the subject, display means for displaying an image generated by the reconstructing means, and reprojection in which the X-ray CT image is made to correspond to the superabsorbent component Reprojection data of the high absorber component is generated by performing a reprojection operation based on a coefficient, and the high absorber component is generated based on the reprojection data of the high absorber component and the projection data of the subject. Projection data and the superabsorbent component contained in the subject. The projection data of the material component is generated, and correction for reducing artifacts caused by X-ray quality hardening is performed on each of the projection data of the superabsorber component and the projection data of the material component. A high-absorber correction unit that performs correction, and the reconstruction unit corrects the artifact by performing a reconstruction calculation based on the corrected projection data of the high-absorber component and the projection data of the substance component. An X-ray CT image is generated, and the display means displays an X-ray CT image in which the artifact is corrected.
上記「X線の線質硬化に起因するアーチファクトを低減するための補正」は、ビームハードニング補正と呼ばれることもあるが、その名称を問わず、「X線の線質硬化に起因するアーチファクトを低減する」ことを目的とする補正は、上記の補正に含まれるものとする。 The above “correction for reducing artifacts caused by X-ray quality hardening” is sometimes called beam hardening correction, but regardless of its name, “correction for artifacts caused by X-ray quality hardening” Corrections intended to “reduce” are included in the above corrections.
また、本明細書において、「高吸収体成分を除く物質成分」ではなく、これを略して「物質成分」と短く記載することがある。 In the present specification, instead of “substance component excluding superabsorbent component”, this may be abbreviated as “substance component” for short.
また「成分」は、高吸収体成分と、高吸収体成分を除く物質成分(「物質成分」と略して記載することもある)と、の両方について言及している用語として用いる。 The “component” is used as a term referring to both the superabsorber component and the substance component excluding the superabsorber component (sometimes abbreviated as “substance component”).
本発明によれば、物質成分毎に補正後の投影データを生成した後、再構成演算を行ってX線CT像を生成するため、物質の成分の数が増えても、補正後のX線CT像を得るために必要な再構成演算は、被検体の投影データを再構成する時と補正後のX線CT像を再構成する時との合計2回になる。そのため、高吸収体成分アーチファクトの低減に要する演算時間を短くし、かつ高精度な高吸収体補正を行うことで、CT値定量性に優れた高度な高吸収体成分アーチファクト低減機能を有するX線CT装置を提供することができる。
特に、被検体に2成分以上の高吸収体成分が含まれる時、本発明によるアーチファクト低減機能はより顕著に有効である。
According to the present invention, the corrected projection data is generated for each substance component, and then the reconstruction calculation is performed to generate an X-ray CT image. Therefore, even if the number of substance components increases, the corrected X-ray The reconstruction calculation necessary for obtaining a CT image is a total of two times when reconstructing the projection data of the subject and when reconstructing the corrected X-ray CT image. Therefore, X-rays with advanced superabsorber component artifact reduction function with excellent CT value quantification by shortening the calculation time required to reduce superabsorber component artifacts and performing high-accuracy high-absorber correction A CT apparatus can be provided.
In particular, when the subject contains two or more superabsorbent components, the artifact reduction function according to the present invention is significantly more effective.
以下、添付図面を用いて本発明に係るX線CT装置の実施の形態について詳説する。
まず、本発明が適用されるコーンビームX線CT装置の概略的な構成について説明する。
Hereinafter, embodiments of the X-ray CT apparatus according to the present invention will be described in detail with reference to the accompanying drawings.
First, a schematic configuration of a cone beam X-ray CT apparatus to which the present invention is applied will be described.
図1は、Cアーム方式と称されるコーンビームX線CT装置1の概略構成図であり、被検体2に対してX線を照射し、被検体2のX線透過像を撮影してX線画像データを得る撮影部10と、撮影部10の各構成要素を制御したり、X線画像データに基づいて被検体2の3次元的X線CT像を再構成したりする制御演算部20を備える。また、画像を表示する表示装置80と、表示装置80に表示された画像の位置を移動するための指示を入力するための、マウス、キーボード、あるいはトラックボール等からなる情報入力装置70とを備える。 FIG. 1 is a schematic configuration diagram of a cone-beam X-ray CT apparatus 1 called a C-arm method. X-rays are irradiated onto a subject 2 and an X-ray transmission image of the subject 2 is taken to obtain an X-ray. An imaging unit 10 that obtains line image data, and a control calculation unit 20 that controls each component of the imaging unit 10 and reconstructs a three-dimensional X-ray CT image of the subject 2 based on the X-ray image data. Is provided. Further, the display device 80 includes an image display device 80 and an information input device 70 including a mouse, a keyboard, a trackball, or the like for inputting an instruction for moving the position of the image displayed on the display device 80. .
図1aは、座位方式と称されるコーンビームX線CT装置1aの概略構成図であり、撮影部10aと、撮影部10aの各構成要素を制御したり、3次元的X線CT像を再構成したりする制御演算部20aとを備える。 FIG. 1a is a schematic configuration diagram of a cone beam X-ray CT apparatus 1a called a sitting system. The imaging unit 10a and each component of the imaging unit 10a are controlled, and a three-dimensional X-ray CT image is reproduced. Or a control arithmetic unit 20a configured.
図1bは、いわゆる被検体回転方式と称されるコーンビームX線CT装置1bを示している。X線源11と2次元X線検出器12が、被検体2の周りを回転する代わりに、X線源11、2次元X線検出器12が床に固定され、回転台19に載置された被検体2が自転するようになっている。図1aの場合と比較して上記回転動作を除いてほぼ類似の構成となっている。 FIG. 1 b shows a cone beam X-ray CT apparatus 1 b called a so-called subject rotation method. Instead of the X-ray source 11 and the two-dimensional X-ray detector 12 rotating around the subject 2, the X-ray source 11 and the two-dimensional X-ray detector 12 are fixed to the floor and placed on the turntable 19. The subject 2 rotates. Compared to the case of FIG. 1a, the configuration is substantially similar except for the above-described rotation operation.
図1cは、移動型X線装置に搭載されたCアーム方式コーンビームX線CT装置1cを示している。コーンビームX線CT装置1cに、車輪9が搭載され、検査室、手術室を移動できるようになっている。図1の場合と比較してCアームが上記移動型X線装置に搭載されている点を除いてほぼ類似の構成となっている。 FIG. 1c shows a C-arm type cone beam X-ray CT apparatus 1c mounted on a mobile X-ray apparatus. A wheel 9 is mounted on the cone beam X-ray CT apparatus 1c so that it can move between an examination room and an operating room. Compared to the case of FIG. 1, the configuration is substantially similar except that the C-arm is mounted on the mobile X-ray apparatus.
以下、図1に示す各構成要素を主に説明し、必要に応じ、図1a、図1b及び図1cに示す構成要素を説明する。 Hereinafter, each component shown in FIG. 1 will be mainly described, and the components shown in FIGS. 1a, 1b, and 1c will be described as necessary.
(撮影部10)
撮影部10は、寝台17と、該寝台17に横臥された被検体2にX線を照射するX線源11と、該X線源11に対向して設置され被検体2を透過したX線を検出することによりX線画像データを出力する2次元X線検出器12と、X線源11及び2次元X線検出器12を機械的に接続するC型アーム13と、該C型アーム13を保持するC型アーム保持体14と、該C型アーム保持体14を天井に取り付ける天井支持体15と、該天井支持体15を図示の状態で前後左右の2次元方向に移動可能に支持する天井レール16と、被検体2に造影剤を注入するインジェクタ18とを備える。
(Shooting unit 10)
The imaging unit 10 includes a bed 17, an X-ray source 11 that irradiates the subject 2 lying on the bed 17 with X-rays, and an X-ray that is installed facing the X-ray source 11 and passes through the subject 2. A two-dimensional X-ray detector 12 that outputs X-ray image data by detecting the X-ray, a C-type arm 13 that mechanically connects the X-ray source 11 and the two-dimensional X-ray detector 12, and the C-type arm 13 A C-type arm holding body 14 for holding the ceiling, a ceiling support 15 for attaching the C-type arm holding body 14 to the ceiling, and supporting the ceiling support 15 movably in the front, rear, left and right directions in the illustrated state. A ceiling rail 16 and an injector 18 for injecting a contrast medium into the subject 2 are provided.
図1aに示す、座位方式のコーンビームX線CT装置1aの撮影部10aは、椅子7と、椅子7に載せられた被検体2にX線を照射するX線源11と、該X線源11に対向して設置され被検体2を透過したX線を検出することによりX線画像データを出力する2次元X線検出器12と、X線源11及び2次元X線検出器12を機械的に接続する旋回アーム5と、該旋回アーム5を保持する支柱6と、支柱6及び椅子7を保持する固定架台8とを備える。 An imaging unit 10a of a sitting-type cone beam X-ray CT apparatus 1a shown in FIG. 1a includes a chair 7, an X-ray source 11 that irradiates the subject 2 placed on the chair 7, and the X-ray source. A two-dimensional X-ray detector 12 that outputs X-ray image data by detecting X-rays installed opposite to the subject 2 and transmitted through the subject 2, an X-ray source 11, and a two-dimensional X-ray detector 12. The revolving arm 5 is connected, the support 6 holding the revolving arm 5, and the fixed base 8 holding the support 6 and the chair 7.
X線源11は、X線を発生するX線管11tと、X線管11tからのX線照射の方向を円錐または四角錐状に制御するコリメータ11cとを備える。
2次元X線検出器12は、たとえばTFT素子を用いるフラットパネルディテクター「FPD」(図1、図1a、図1b)、または、X線透過像を可視光像に変換するX線イメージインテンシファイア12i(図1c)とX線イメージインテンシファイア12iの像を結像する光学レンズ(図示しない)、及び光学レンズにより結像されたX線イメージインテンシファイア12iの可視光像を撮影するCCDテレビカメラ12c(図1c)の組み合わせ等からなる2次元X線検出器などが用いられる。そして、2次元X線検出器の形状は円形、方形いかなる形状であってもよい。
The X-ray source 11 includes an X-ray tube 11t that generates X-rays, and a collimator 11c that controls the direction of X-ray irradiation from the X-ray tube 11t in a cone or quadrangular pyramid shape.
The two-dimensional X-ray detector 12 is, for example, a flat panel detector “FPD” (FIG. 1, FIG. 1a, FIG. 1b) using a TFT element, or an X-ray image intensifier that converts an X-ray transmission image into a visible light image. 12i (FIG. 1c) and an optical lens (not shown) that forms an image of the X-ray image intensifier 12i, and a CCD television that captures a visible light image of the X-ray image intensifier 12i formed by the optical lens A two-dimensional X-ray detector composed of a combination of cameras 12c (FIG. 1c) or the like is used. The shape of the two-dimensional X-ray detector may be circular or square.
上記旋回アーム5及びC型アーム13は、被検体2の撮影に際して、所定の投影角度毎に回転中心軸4を中心として回転移動する。これにより、上記X線源11及び2次元X線検出器12は、ほぼ同一の平面上にある円軌道で回転移動しながら、X線撮影を行う。この回転移動については、画像再構成演算に使用される幾何学パラメータが存在する。幾何学パラメータとして、旋回アーム5またはC型アーム13が回転移動することにより、X線源11と2次元X線検出器12とが描く円軌道を含む面である回転軌道面(ミッドプレーン)3と、回転中心軸4が挙げられる。 When the subject 2 is imaged, the swivel arm 5 and the C-shaped arm 13 rotate around the rotation center axis 4 at every predetermined projection angle. Thereby, the X-ray source 11 and the two-dimensional X-ray detector 12 perform X-ray imaging while rotating and moving on a circular orbit on substantially the same plane. For this rotational movement, there are geometric parameters used for the image reconstruction operation. As a geometric parameter, a rotational trajectory plane (midplane) 3 which is a plane including a circular trajectory drawn by the X-ray source 11 and the two-dimensional X-ray detector 12 by the rotational movement of the swivel arm 5 or the C-shaped arm 13. And the rotation center axis 4.
(制御演算部20)
制御演算部20は、撮影部10を制御する撮影部制御手段100と、撮影部10が出力したX線画像データを収集して格納する画像収集手段110と、収集されたX線画像データに基づいて3次元的X線CT像を再構成する再構成手段200と、本発明の特徴である高吸収体補正投影データを生成する高吸収体補正手段300と、再構成手段200で生成した3次元的X線CT像を表示する画像表示手段120とを備える。
(Control operation unit 20)
The control calculation unit 20 is based on the imaging unit control unit 100 that controls the imaging unit 10, the image collection unit 110 that collects and stores the X-ray image data output by the imaging unit 10, and the collected X-ray image data. The reconstruction unit 200 for reconstructing a three-dimensional X-ray CT image, the high-absorber correction unit 300 for generating the high-absorber correction projection data, which is a feature of the present invention, and the three-dimensional generated by the reconstruction unit 200 And an image display means 120 for displaying a standard X-ray CT image.
(撮影部制御手段100)
撮影部制御手段100は、C型アーム13の、回転中心軸4の回りの回転移動(以下「プロペラ回転」という)を制御する撮影系回転制御手段101と、天井支持体15の天井レール16上での位置を制御してC型アーム13の被検体2に対する位置を2次元的に制御する撮影系位置制御手段102と、X線管11tに流す管電流のON、OFF等を制御するX線照射制御手段103と、インジェクタ18が被検体2に注入する造影剤の注入量及び注入タイミングを制御するインジェクタ制御手段104と、寝台17の位置を制御して被検体2の位置を調整するための寝台制御手段106と、2次元X線検出器12によるX線透過像の撮影を制御する検出系制御手段107とを備える。
(Shooting unit control means 100)
The imaging unit control unit 100 includes an imaging system rotation control unit 101 that controls the rotational movement of the C-arm 13 around the rotation center axis 4 (hereinafter referred to as “propeller rotation”), and the ceiling support 16 on the ceiling rail 16. Imaging system position control means 102 that controls the position of the C-arm 13 relative to the subject 2 in a two-dimensional manner, and X-rays that control ON / OFF of the tube current flowing through the X-ray tube 11t. The irradiation control means 103, the injector control means 104 for controlling the injection amount and injection timing of the contrast medium injected into the subject 2 by the injector 18, and the position of the bed 17 for adjusting the position of the subject 2 A couch control means 106 and a detection system control means 107 for controlling radiographic image capturing by the two-dimensional X-ray detector 12 are provided.
図1aに示す、座位方式のコーンビームX線CT装置1aの撮影部制御手段100aは、旋回アーム6の、回転中心軸4の回りの回転移動を制御する撮影系回転制御手段101と、X線照射制御手段103と、椅子7の位置を制御して被検体2の位置を調整するための椅子位置制御手段105と、検出系制御手段107とを備える。 The imaging unit control means 100a of the sitting-type cone beam X-ray CT apparatus 1a shown in FIG. 1a includes an imaging system rotation control means 101 for controlling the rotational movement of the swivel arm 6 around the rotation center axis 4, and an X-ray. Irradiation control means 103, chair position control means 105 for adjusting the position of the subject 2 by controlling the position of the chair 7, and detection system control means 107 are provided.
(再構成手段200)
再構成手段200は、前処理手段210と、第1のビームハードニング補正手段230と、フィルタリング手段270と、逆投影手段280とを備える。
(Reconstruction means 200)
The reconstruction unit 200 includes a preprocessing unit 210, a first beam hardening correction unit 230, a filtering unit 270, and a back projection unit 280.
前処理手段210は、画像収集手段110が収集したX線画像データをX線吸収係数の分布像(以下「投影データ220」という)に変換する。本実施の形態では、まず、被検体2及び寝台17を撮影視野内に配置しない状態で予め撮影された空気のX線透過像の各画素データに対して自然対数変換演算を施す。次に被検体2を寝台17に載せた状態で撮影したX線透過像の各画素データに対して自然対数変換演算を施す。そして、上記2つの自然対数変換演算を施したX線透過画像の差分を取ることにより、被検体2及び寝台17の投影データ220を得る。 The pre-processing unit 210 converts the X-ray image data collected by the image collecting unit 110 into an X-ray absorption coefficient distribution image (hereinafter referred to as “projection data 220”). In the present embodiment, first, natural logarithmic conversion calculation is performed on each pixel data of an X-ray transmission image of air that has been imaged in advance without the subject 2 and the bed 17 being placed in the field of view. Next, a natural logarithmic conversion operation is performed on each pixel data of an X-ray transmission image taken with the subject 2 placed on the bed 17. Then, the projection data 220 of the subject 2 and the bed 17 is obtained by taking the difference between the two X-ray transmission images subjected to the natural logarithmic transformation calculation.
第1のビームハードニング補正手段230は、投影データ220にある非線形変換を施し、第1のビームハードニング補正投影データ240を生成する。ここで「第1の」ビームハードニング補正と記しているのは、後述の「物質別の」ビームハードニング補正と名称を区別するためであり、便宜的なものである。なお、X線のビームハードニングを補正するための非線形変換の方法は、撮影管電圧、X線源11のX線スペクトル特性、及び2次元X線検出器12の応答特性により異なり、個々に実験的に求めるべきものである。したがって、非線形変換の具体的な方法を明らかにすることは本明細書の目的ではない。 The first beam hardening correction unit 230 performs nonlinear conversion on the projection data 220 to generate first beam hardening correction projection data 240. The term “first” beam hardening correction is used here to distinguish the name from “substance-specific” beam hardening correction, which will be described later, for convenience. Note that the nonlinear conversion method for correcting X-ray beam hardening differs depending on the tube voltage, the X-ray spectral characteristics of the X-ray source 11, and the response characteristics of the two-dimensional X-ray detector 12. Should be required. Therefore, it is not the purpose of this specification to clarify the specific method of nonlinear transformation.
そして、オリジナルの投影データ220(Prj)と、第1のビームハードニング補正投影データ240(BHC1)は、後述する高吸収体補正手段300の中で、物質含有率を算出するために使用するため、ワークエリアに保存しておく。 Then, the original projection data 220 (Prj) and the first beam hardening correction projection data 240 (BHC 1 ) are used for calculating the substance content rate in the high-absorber correction unit 300 described later. Therefore, save it in the work area.
また、この第1のビームハードニング補正手段230は、省略することも可能である。この場合、後述する第1の再構成CT像291の一様性が低下することにより、再投影係数入力手段310において設定する、再投影係数の精度が低下する場合も想定されるが、省略は可能である。第1のビームハードニング補正手段230を省略する場合は、上述の投影データ220(Prj)と、投影データ240(BHC1)は同一の画像データとなり、後述する図5のステップS331、S332、S333における物質の成分ごとの含有率算出工程と、ステップS334、S335におけるメタル成分投影データ生成処理及び骨成分投影データ生成処理を省略することができる。すなわち、ステップS323で生成されるメタル再投影データ323(FpjM)が、ステップS334におけるメタル成分投影データ(M)となり、ステップS330で生成される骨成分再投影データ(FpjB)がステップS335における骨成分投影データ(B)に相当する。また、軟部組織成分の投影データは、ステップS336における軟部組織成分投影データ生成処理において、メタル再投影データ323(FpjM)と骨成分再投影データ(FpjB)とを投影データ220(Prj)から減算することにより、軟部組織成分投影データ(N)が生成される。 Further, the first beam hardening correction means 230 can be omitted. In this case, it is assumed that the accuracy of the reprojection coefficient set in the reprojection coefficient input unit 310 is reduced due to a decrease in the uniformity of the first reconstructed CT image 291 to be described later. Is possible. When the first beam hardening correction unit 230 is omitted, the projection data 220 (Prj) and the projection data 240 (BHC 1 ) are the same image data, and steps S331, S332, and S333 in FIG. It is possible to omit the content rate calculation step for each component of the substance and the metal component projection data generation processing and bone component projection data generation processing in steps S334 and S335. That is, the metal reprojection data 323 (Fpj M ) generated in step S323 becomes the metal component projection data (M) in step S334, and the bone component reprojection data (Fpj B ) generated in step S330 is in step S335. This corresponds to the bone component projection data (B). In addition, the soft tissue component projection data is obtained from the metal reprojection data 323 (Fpj M ) and the bone component reprojection data (Fpj B ) from the projection data 220 (Prj) in the soft tissue component projection data generation processing in step S336. By subtracting, soft tissue component projection data (N) is generated.
フィルタリング手段270は、X線CT画像再構成におけるフィルタリング処理を行う。
逆投影手段280は、フィルタリング処理後の投影データに対し、例えば非特許文献2に記載のフェルドカンプの方法に基づいて逆投影演算を行い、3次元的CT像(以下「再構成CT像290」という)を生成する。再構成CT像290の出力は、通常はアキシャル断面像を積層したものとして行われる。
The filtering unit 270 performs filtering processing in X-ray CT image reconstruction.
The back projection means 280 performs a back projection operation on the projection data after filtering processing based on, for example, the Feldkamp method described in Non-Patent Document 2, and performs a three-dimensional CT image (hereinafter, “reconstructed CT image 290”). Generated). The output of the reconstructed CT image 290 is normally performed as a stack of axial cross-sectional images.
(高吸収体補正手段300)
高吸収体補正手段300は、本発明の特徴である高吸収体補正投影データを生成する手段であり、再投影係数入力手段310と、再投影手段320と、物質含有率算出手段330と、物質別ビームハードニング補正手段340とを備える。
(High Absorber Correction Unit 300)
The high-absorber correction unit 300 is a unit that generates high-absorber correction projection data, which is a feature of the present invention, and includes a reprojection coefficient input unit 310, a reprojection unit 320, a substance content rate calculation unit 330, And another beam hardening correction means 340.
再投影係数入力手段310より、再投影係数を設定する。再投影係数の入力には、後述するような画面を使用する場合と、あらかじめ高吸収体補正前の「第1の」再構成CT値と再投影値の変換テーブルを生成する再投影データの個数分用意しておき、読み込ませるという形で実現してもよい。 A reprojection coefficient is set from the reprojection coefficient input means 310. When inputting a reprojection coefficient, a screen as described later is used, and the number of reprojection data for generating a conversion table of “first” reconstructed CT values and reprojection values before correction of the high absorber. It may be realized in the form of preparing and reading in minutes.
再投影手段320は、再投影係数入力手段310が設定した再投影係数に基づいて再投影演算を行う。 The reprojection unit 320 performs a reprojection calculation based on the reprojection coefficient set by the reprojection coefficient input unit 310.
物質含有率算出手段330は、再投影データと、上記再構成手段200で保存したオリジナルの投影データ220(Prj)、第1のビームハードニング補正投影データ240(BHC1)に基づいて、投影データの各点の物質含有率を算出し、各物質成分の投影データを生成する。 Based on the reprojection data, the original projection data 220 (Prj) stored in the reconstruction unit 200, and the first beam hardening corrected projection data 240 (BHC 1 ), the substance content rate calculation unit 330 generates projection data. The substance content rate of each point is calculated, and projection data of each substance component is generated.
物質別ビームハードニング補正手段340は、物質含有率算出手段330が生成した各物質成分の投影データに、各物質固有のビームハードニング補正を施す。そして、ビームハードニング補正が施された各物質成分の投影データを加算し、高吸収体補正投影データ360(BHC2)を生成する。 The material-specific beam hardening correcting unit 340 performs beam hardening correction specific to each material on the projection data of each material component generated by the material content rate calculating unit 330. Then, the projection data of each material component that has undergone beam hardening correction is added to generate high-absorber correction projection data 360 (BHC 2 ).
生成された高吸収体補正投影データ360(BHC2)は、フィルタリング手段270に送られ、逆投影手段280により、高吸収体成分アーチファクトが補正された再構成CT像290を出力する。 The generated superabsorber corrected projection data 360 (BHC 2 ) is sent to the filtering unit 270, and the back projection unit 280 outputs a reconstructed CT image 290 in which the superabsorber component artifact is corrected.
上記のコーンビームX線CT装置1の仕様例は次のとおりである。X線管11tと回転中心軸4との距離は800mm、回転中心軸4と2次元X線検出器FPDのX線入射面との距離は400mm、FPDのX線入射面の大きさは400mm×300mmの長方形状であって、画像サイズは2048×1536、そして画素ピッチは0.2mmである。FPDにX線が入射すると、まずX線入射面でCsI等の発光体により光に変換され、光信号はフォトダイオードにより電荷に変換される。蓄積した電荷は一定のフレームレートごとにTFT素子によりデジタル信号に変換され、読み出される。回転撮影モードでは、毎秒30フレーム、画像サイズ1024×768で2次元X線画像データを読み出す。撮影系回転制御手段101は、2次元X線検出器12を、被検体2の左手の方向(−100度)から天井方向(0度)を通過し、被検体2の右手方向(+100度)まで移動させる。これにより、200度の投影角度にわたる被検体2の2次元X線画像データを撮影する。C型アーム13の回転速度の代表例は1秒当たり40度で、スキャン時間は例えば5秒である。 A specification example of the cone beam X-ray CT apparatus 1 is as follows. The distance between the X-ray tube 11t and the rotation center axis 4 is 800 mm, the distance between the rotation center axis 4 and the X-ray incidence surface of the two-dimensional X-ray detector FPD is 400 mm, and the size of the FPD X-ray incidence surface is 400 mm × It has a rectangular shape of 300 mm, the image size is 2048 × 1536, and the pixel pitch is 0.2 mm. When X-rays enter the FPD, the light is first converted into light by a light emitter such as CsI on the X-ray incident surface, and the optical signal is converted into electric charge by a photodiode. The accumulated charge is converted into a digital signal by a TFT element at a constant frame rate and read out. In the rotational imaging mode, two-dimensional X-ray image data is read at 30 frames per second and an image size of 1024 × 768. The imaging system rotation control means 101 passes through the two-dimensional X-ray detector 12 from the direction of the left hand of the subject 2 (−100 degrees) to the ceiling direction (0 degrees) and to the right hand direction of the subject 2 (+100 degrees). To move. Thereby, two-dimensional X-ray image data of the subject 2 over a projection angle of 200 degrees is captured. A typical example of the rotational speed of the C-arm 13 is 40 degrees per second, and the scan time is, for example, 5 seconds.
また図1aに示す、座位方式コーンビームX線CT装置1aの回転の仕様例は次のとおりである。撮影系回転制御手段101は、X線源11と2次元X線検出器12を被検体2の周りに回転させ、検出系制御手段107により、被検体2の360°の投影データを撮影する。旋回アーム5の回転速度の代表例は1秒当たり37.5度で、スキャン時間は9.6秒である。 An example of the rotation specification of the sitting-type cone beam X-ray CT apparatus 1a shown in FIG. 1a is as follows. The imaging system rotation control means 101 rotates the X-ray source 11 and the two-dimensional X-ray detector 12 around the subject 2, and the detection system control means 107 images 360 ° projection data of the subject 2. A typical example of the rotation speed of the swivel arm 5 is 37.5 degrees per second, and the scan time is 9.6 seconds.
次に、コーンビームX線CT装置1による撮影における動作の概要について説明する。 Next, an outline of operations in imaging by the cone beam X-ray CT apparatus 1 will be described.
コーンビームX線CT装置1では、先ず撮影系回転制御手段101はC型アーム13のプロペラ回転を開始する。回転加速期間を経たのち、X線照射制御手段103はX線管11tからX線を照射し、検出系制御手段107は2次元X線検出器12による撮像を開始する。X線管11tから照射されたX線は、被検体2を透過した後、2次元X線検出器12に取り込まれる。2次元X線検出器12の信号は、A/D変換を経た後、デジタル信号からなる2次元のX線画像データとして画像収集手段110に記録される。2次元X線検出器FPDの標準走査モードは毎秒30フレームである。撮影による投影角度間隔は1.33度で、5秒間に150枚のX線透過像を取得する。200度の回転撮影が完了すると、X線照射制御手段103はX線管11tのX線照射を終了し、撮影系回転制御手段101は回転減速期間を経たのち回転を停止する。 In the cone beam X-ray CT apparatus 1, the imaging system rotation control means 101 first starts propeller rotation of the C-arm 13. After the rotation acceleration period, the X-ray irradiation control means 103 emits X-rays from the X-ray tube 11t, and the detection system control means 107 starts imaging by the two-dimensional X-ray detector 12. X-rays irradiated from the X-ray tube 11t pass through the subject 2 and are taken into the two-dimensional X-ray detector 12. The signal from the two-dimensional X-ray detector 12 undergoes A / D conversion and is then recorded in the image collecting unit 110 as two-dimensional X-ray image data composed of a digital signal. The standard scanning mode of the two-dimensional X-ray detector FPD is 30 frames per second. The projection angle interval by imaging is 1.33 degrees, and 150 X-ray transmission images are acquired in 5 seconds. When the rotation imaging at 200 degrees is completed, the X-ray irradiation control unit 103 ends the X-ray irradiation of the X-ray tube 11t, and the imaging system rotation control unit 101 stops the rotation after a rotation deceleration period.
座位方式のコーンビームX線CT装置1aによる撮影の動作は次のとおりである。先ず撮影系回転制御手段101は旋回アーム5の回転を開始する。22.5度の回転加速期間を経たのち、X線照射制御手段103は、X線管11tからX線を照射し、検出系制御手段107は撮像を開始する。撮影による投影角度間隔は1.25度で、9.6秒間に288枚の2次元X線画像データを取得する。360度の回転撮影が完了すると、X線照射制御手段103はX線管11tからのX線照射を終了し、撮影系回転制御手段101は22.5度の回転減速期間を経たのち回転を停止する。 The imaging operation by the sitting-type cone beam X-ray CT apparatus 1a is as follows. First, the imaging system rotation control means 101 starts to rotate the revolving arm 5. After passing through the rotation acceleration period of 22.5 degrees, the X-ray irradiation control means 103 emits X-rays from the X-ray tube 11t, and the detection system control means 107 starts imaging. The projection angle interval by photographing is 1.25 degrees, and 288 pieces of two-dimensional X-ray image data are acquired in 9.6 seconds. When the 360-degree rotation imaging is completed, the X-ray irradiation control means 103 ends the X-ray irradiation from the X-ray tube 11t, and the imaging system rotation control means 101 stops the rotation after passing through a rotation decelerating period of 22.5 degrees. To do.
また、たとえば図1cに示すように、2次元X線検出器として、X線イメージインテンシファイア12iと、光学レンズ、CCDテレビカメラ12cの組み合わせを用いる場合の仕様例は、X線イメージインテンシファイア12iの直径は300mm、CCDテレビカメラ12cの標準走査モードは毎秒60フレーム、走査線数512本、あるいは毎秒30フレーム、走査線数1024本で、光学レンズにより結像されたX線イメージインテンシファイア12iの可視光像を撮影する。CCDテレビカメラ12cが撮影した2次元X線画像は、ビデオ信号に変換した後にA/D変換され、512×512あるいは1024×1024のデジタル画像として画像収集手段110に記録される。 For example, as shown in FIG. 1c, a specification example in the case of using a combination of an X-ray image intensifier 12i, an optical lens, and a CCD TV camera 12c as a two-dimensional X-ray detector is an X-ray image intensifier. The diameter of 12i is 300 mm, and the standard scanning mode of the CCD TV camera 12c is 60 frames per second and 512 scanning lines, or 30 frames per second and 1024 scanning lines, and an X-ray image intensifier imaged by an optical lens. A 12i visible light image is taken. A two-dimensional X-ray image taken by the CCD TV camera 12c is converted into a video signal, A / D converted, and recorded as a 512 × 512 or 1024 × 1024 digital image in the image collecting means 110.
再構成手段200は、以上のような撮影に並行し、あるいは撮影終了後に画像収集手段110から2次元のX線画像データを読み出し、このX線画像データに基づいて画像再構成演算を行い、被検体2の再構成CT像を生成する。画像表示手段120は、再構成CT像を、CRT装置や液晶ディスプレイ装置等からなる表示装置80に表示する。なお画像表示手段120は、画像収集手段110に記録された2次元X線画像データの表示、及び再投影係数入力手段310において再構成CT像を階調表示して再投影係数を入力するためにも使用される。 The reconstruction unit 200 reads the two-dimensional X-ray image data from the image collection unit 110 in parallel with the above-described imaging or after the imaging is completed, performs an image reconstruction calculation based on the X-ray image data, A reconstructed CT image of the specimen 2 is generated. The image display unit 120 displays the reconstructed CT image on a display device 80 including a CRT device or a liquid crystal display device. The image display means 120 displays the two-dimensional X-ray image data recorded in the image collection means 110, and displays the reconstructed CT image in gradation in the reprojection coefficient input means 310 to input the reprojection coefficient. Also used.
図2及び図2cに、複数の異なるX線吸収係数の高吸収体成分を含む被検体を撮影した場合のX線画像データの例を示す。図2は歯科撮影のX線画像の例を示し、23は歯列を示しエナメル質等の高吸収体成分を含んでいる。24は治療歯に埋め込まれている金属であり、21は回転中心軸4の2次元X線検出器12への投影、22はこれらの高吸収体成分を含む断面を示す。また、図2cに、整形外科の処置の一つである、人工股関節置換術と呼ばれる処置を施したX線画像の例を示す。なお、図2cは、主に図1cに記載したCアーム方式コーンビームX線CT装置1cにより得られる画像を示す。25は寛骨臼と呼ばれる骨盤の一部を示し、26は大転子と呼ばれる大腿骨の突起部である。27は人工股関節で、例えば多孔性表面の金属を材料に形成され、骨頭ボール部と金属ステム部から構成される。図2cは、磨耗した大腿骨頭を人工股関節27の骨頭ボールで置き換え、金属ステムを大腿骨に差し込んで固定する、人工股関節置換術の処置を施した画像を示す。22cは、高吸収体成分である、寛骨臼25、大転子26、人工股関節27を含む断面を示す。 2 and 2c show examples of X-ray image data when a subject including a plurality of high-absorber components having different X-ray absorption coefficients is imaged. FIG. 2 shows an example of an X-ray image of dental imaging. Reference numeral 23 denotes a dentition, which contains a high absorbent component such as enamel. Reference numeral 24 denotes a metal embedded in the treatment tooth, 21 denotes a projection of the rotation center axis 4 on the two-dimensional X-ray detector 12, and 22 denotes a cross section including these superabsorbent components. FIG. 2c shows an example of an X-ray image subjected to a procedure called hip replacement, which is one of orthopedic procedures. FIG. 2c mainly shows an image obtained by the C-arm type cone beam X-ray CT apparatus 1c described in FIG. 1c. Reference numeral 25 denotes a part of the pelvis called an acetabulum, and 26 denotes a femoral protrusion called a greater trochanter. Reference numeral 27 denotes an artificial hip joint, which is made of, for example, a metal having a porous surface, and is composed of a bone head ball portion and a metal stem portion. FIG. 2c shows an image of a hip replacement procedure in which the worn femoral head is replaced with a head ball of the hip 27 and a metal stem is inserted into the femur and secured. 22c shows a cross section including the acetabulum 25, the greater trochanter 26, and the artificial hip joint 27, which are superabsorbent components.
図3及び図3cは、図2及び図2cにおける、高吸収体成分を含む断面22及び22cの、再構成CT像30及び30cを示す。29は、金属と金属を結ぶ方向に発生したアーチファクトを示す。以下、図3の、歯列23と治療歯に埋め込まれている金属24がある場合の例を中心に、再構成CT像の高吸収体補正の説明を行うが、本発明は、整形外科撮影(例えば、23を腰椎に、24を金属製人工骨に置き換えた(例:図3cの人工股関節置換術)の場合にも、23を造影剤、24をステントに置き換えた循環器IVR撮影時(図示しない)の再構成X線CT像の高吸収体補正を行う場合にも、そのまま適用できる。 3 and 3c show reconstructed CT images 30 and 30c of the sections 22 and 22c containing the superabsorbent component in FIGS. 2 and 2c. Reference numeral 29 denotes an artifact generated in the direction connecting the metals. Hereinafter, the high-absorber correction of the reconstructed CT image will be described focusing on an example in which there is the dentition 23 and the metal 24 embedded in the treatment tooth in FIG. (For example, in the case where 23 is replaced with a lumbar spine and 24 is replaced with a metal artificial bone (eg, hip replacement shown in FIG. 3c), 23 is contrast medium, and 24 is replaced with a stent during circulatory IVR imaging ( The present invention can also be applied to the case of performing high-absorber correction of a reconstructed X-ray CT image (not shown).
次に、図4−1、図4−2、図4−1c、図4−2cを用い、再投影係数入力手段310が行う再投影係数入力処理を説明する。図4−1、図4−2は一般的な入力画面例を示し、図4−1c、図4−2cは、特に図1cに記載されたCアーム方式コーンビームX線CT装置において表示される画面表示例を示す。 Next, reprojection coefficient input processing performed by the reprojection coefficient input unit 310 will be described with reference to FIGS. 4-1, 4-2, 4-1c, and 4-2c. FIGS. 4-1 and 4-2 show typical input screen examples, and FIGS. 4-1c and 4-2c are particularly displayed in the C-arm type cone beam X-ray CT apparatus described in FIG. 1c. An example of screen display is shown.
図4−1は、再投影係数1を設定する画面を示す図であり、再投影係数1選択タブ41を選択すると現れる。43はユーザが編集可能な再投影係数名表示欄であり、図では再投影係数1に対して「メタル」以上のX線吸収係数を再投影するための係数であることを示す「メタル以上再投影係数」と命名されている。50は再投影下限値設定ポイントであり、この値以上のCT値を持つ再構成領域が再投影演算に使用される。再投影下限値設定ポイント50をマウス等の情報入力装置70を使用して左右にドラッグし、再投影下限CT値を変更する事ができる。51は再投影係数設定ポイントであり、上下左右にドラッグする事で再投影係数を変化させる事ができる。 FIG. 4A is a diagram illustrating a screen for setting the reprojection coefficient 1 and appears when the reprojection coefficient 1 selection tab 41 is selected. Reference numeral 43 denotes a reprojection coefficient name display field that can be edited by the user. In the figure, the reprojection coefficient 1 is a coefficient for reprojecting an X-ray absorption coefficient greater than or equal to “metal” with respect to the reprojection coefficient 1. It is named “projection coefficient”. Reference numeral 50 denotes a reprojection lower limit setting point, and a reconstruction area having a CT value equal to or greater than this value is used for the reprojection calculation. The reprojection lower limit CT value can be changed by dragging the reprojection lower limit setting point 50 to the left or right using an information input device 70 such as a mouse. Reference numeral 51 denotes a reprojection coefficient setting point, which can be changed by dragging up, down, left, and right.
ここで再投影係数とは、X線吸収係数がゼロ(HU、所謂ハウスフィールドの定義によるCT値が−1000)となる点を原点として所定のCT値(を、HUで定義されているCT値という数値をプラス1000して、正しい物理量に修正して)から計算される各ボクセル点の再投影値に掛算する係数を表しており、通常0〜1の範囲の値を設定する。52は再投影上限値設定ポイントであり、図では閾値Th3以上の再投影係数は1に固定されている。 Here, the reprojection coefficient is a predetermined CT value (the CT value defined by the HU) with the point where the X-ray absorption coefficient is zero (the CT value defined by the so-called house field is -1000). This is a coefficient that is multiplied by the reprojection value of each voxel point calculated by adding 1000 to the correct physical quantity and normally setting a value in the range of 0-1. Reference numeral 52 denotes a reprojection upper limit setting point. In the figure, a reprojection coefficient equal to or greater than the threshold Th3 is fixed to 1.
ポイント50−52を結ぶ線は、「メタル」領域を再投影するための再投影係数カーブを形成する。図では、折れ線で結ばれた再投影係数カーブを表示しているが、ポイント50−52をスプライン曲線で結んで再投影係数カーブを形成するようにしてもよい。53、54はそれぞれ、ポイント50−52の位置でのCT値を設定する再投影閾値設定欄、再投影係数設定欄であり、ポイント50−52をドラッグする代わりに53、54の値を直接編集する事により、再投影係数カーブを変更する事ができるようになっている。 The line connecting points 50-52 forms a reprojection coefficient curve for reprojecting the “metal” region. In the figure, a reprojection coefficient curve connected by a broken line is displayed, but a point 50-52 may be connected by a spline curve to form a reprojection coefficient curve. 53 and 54 are a reprojection threshold value setting field and a reprojection coefficient setting field for setting the CT value at the position of the point 50-52, respectively. Instead of dragging the point 50-52, the values of 53 and 54 are directly edited. By doing so, the reprojection coefficient curve can be changed.
55は再投影係数設定ポイント追加ボタンであり、クリックすると、現在選択されている再投影係数カーブの線分の間に、新たなポイントが現れるようになっている。逆に折れ線カーブを形成するポイントを減らすときは、ポイントの上を右クリックし削除を選択(図示しない)する。図の上限値(Th3)以上の再投影係数は1とするのが通常ではあるが、高吸収体成分による再投影値の寄与を軽減するためにTh3以上の再投影係数を0に近づけたり、Th3以上の再投影値と再投影係数の積を一定するために、Th3以上の再投影係数カーブ1/Xにしたり、そもそも図の縦軸を再投影係数(微分係数)でなく、CT値−再投影値(積分値)に変更したグラフにしてもよい。 55 is a reprojection coefficient setting point addition button, and when clicked, a new point appears between the line segments of the currently selected reprojection coefficient curve. Conversely, to reduce the number of points that form a polygonal curve, right-click on the point and select Delete (not shown). Normally, the reprojection coefficient equal to or greater than the upper limit (Th3) in the figure is set to 1, but in order to reduce the contribution of the reprojection value due to the high absorber component, the reprojection coefficient equal to or greater than Th3 is brought close to 0, In order to make the product of the reprojection value greater than Th3 and the reprojection coefficient constant, the reprojection coefficient curve 1 / X greater than Th3 is used. In the first place, the vertical axis is not the reprojection coefficient (differential coefficient) but the CT value− The graph may be changed to a reprojection value (integral value).
図4−2は、再投影係数2を設定する画面を示す図であり、再投影係数2選択タブ42を選択すると現れる。再投影係数名表示欄43には、「骨」以上のX線吸収係数を再投影するための係数であることを示す「骨以上再投影係数」と命名されている。57は現在の再投影係数2の設定カーブであり、56の白抜き線は再投影係数1の設定カーブを現わす。再投影係数2設定カーブ57を構成するポイントは、後述するステップS330により、「骨以上」再投影データからメタル再投影データを引算する事により骨成分再投影データを生成するため、再投影係数2で設定される「骨以上」再投影値が、再投影係数1で設定されるメタル再投影値を下回らないように制限されている。 FIG. 4B is a diagram illustrating a screen for setting the reprojection coefficient 2 and appears when the reprojection coefficient 2 selection tab 42 is selected. The reprojection coefficient name display field 43 is named “bone or more reprojection coefficient” indicating that it is a coefficient for reprojecting an X-ray absorption coefficient of “bone” or more. Reference numeral 57 denotes a current setting curve for the reprojection coefficient 2, and 56 white lines indicate the setting curve for the reprojection coefficient 1. The points constituting the reprojection coefficient 2 setting curve 57 are reprojection coefficients in order to generate bone component reprojection data by subtracting metal reprojection data from the “projection above bone” reprojection data in step S330 described later. The “projection above bone” reprojection value set at 2 is limited so as not to fall below the metal reprojection value set at the reprojection coefficient 1.
図4−1c、図4−2cは、整形外科の、人工股関節置換術の場合の、再投影係数入力画面を示す。寛骨臼25、大転子26、人工股関節27を含む断面の再構成CT像を、図4−1cではメタル閾値表示、図4−2cでは骨閾値表示している。この場合も、再投影係数を同様に設定し、再投影データを生成することができる。 FIGS. 4-1c and 4-2c show reprojection coefficient input screens in the case of orthopedic hip replacement. A reconstructed CT image of a cross section including the acetabulum 25, greater trochanter 26 and artificial hip joint 27 is displayed as a metal threshold value in FIG. 4-1c and as a bone threshold value in FIG. 4-2c. Also in this case, the reprojection data can be generated by setting the reprojection coefficient in the same manner.
なお、図4−1、図4−2において、再投影係数追加ボタン44をクリックすると、再投影係数選択タブが1つ増え、「第3」の再投影係数を入力する画面(例えば「繊維質」以上のX線吸収係数を「再投影係数設定」をする)を生成することができる。以下は、メタルと、骨以上の2つの再投影係数を入力して、高吸収体補正投影データを生成する場合について、詳しく説明する。 In FIGS. 4A and 4B, when the reprojection coefficient addition button 44 is clicked, the reprojection coefficient selection tab is incremented by one, and a screen for inputting the “third” reprojection coefficient (for example, “fiber structure” The above-mentioned X-ray absorption coefficient can be “re-projected coefficient setting”). In the following, the case where the high-absorber correction projection data is generated by inputting the metal and two reprojection coefficients equal to or higher than the bone will be described in detail.
次に、図5のフローチャートを用いて、本発明の特徴である高吸収体補正手段300が高吸収体補正投影データ360を生成する演算処理を具体的に説明する。 Next, with reference to the flowchart of FIG. 5, a specific description will be given of the arithmetic processing in which the high-absorber correction unit 300, which is a feature of the present invention, generates the high-absorber correction projection data 360.
上述した前処理手段210が施された投影データ(Prj)220が入力画像データとなり、ステップS230が開始される。 The projection data (Prj) 220 that has been subjected to the preprocessing means 210 described above becomes input image data, and step S230 is started.
(ステップS230)
第1のビームハードニング補正手段230により、投影データ220(Prj)にある非線形変換を施し、第1のビームハードニング補正投影データ240(BHC1)を生成する(S230)。ここで行う非線形変換により、図4−1、図4−2を用いて上述した再投影係数入力手段310において、第1の再構成CT像291のメタル閾値や骨閾値表示する場合に表示画像を一様にし、かつCT値をある程度正しい値に補正しておく事で、再投影係数の調整精度を向上する事ができる。
(Step S230)
The first beam hardening correction unit 230 performs nonlinear transformation on the projection data 220 (Prj) to generate first beam hardening correction projection data 240 (BHC 1 ) (S230). By the non-linear transformation performed here, the display image is displayed when the metal threshold value or the bone threshold value of the first reconstructed CT image 291 is displayed in the reprojection coefficient input unit 310 described above with reference to FIGS. It is possible to improve the reprojection coefficient adjustment accuracy by making the CT values uniform and correct the CT values to some extent.
上記「第1の再構成CT像291」とは、後述する高吸収体補正投影データ(BHC2)を再構成して得られる再構成CT像290と区別するために定義した再構成CT像であって、物質成分別のビームハードニング補正をすることなく、投影データ220(Prj)に対してビームハードニング補正をした第1のビームハードニング補正投影データ240(BHC1)を再構成して得られたCT像を意味する。同様に、「第1のビームハードニング補正」とは、後述する物質成分別のビームハードニング補正と区別するために定義したビームハードニング補正であって、物質成分を考慮することなく投影データ220(Prj)に対して施すビームハードニング補正を意味する。 The “first reconstructed CT image 291” is a reconstructed CT image defined to distinguish it from a reconstructed CT image 290 obtained by reconstructing high-absorber corrected projection data (BHC 2 ) described later. Then, the first beam hardening corrected projection data 240 (BHC 1 ) obtained by performing the beam hardening correction on the projection data 220 (Prj) without performing the beam hardening correction for each material component is reconstructed. It means the obtained CT image. Similarly, the “first beam hardening correction” is a beam hardening correction defined to be distinguished from a beam hardening correction for each material component, which will be described later, and includes the projection data 220 without considering the material component. This means beam hardening correction to be applied to (Prj).
非線形変換は、図では原点を通る3次関数を例に示してあるが、N次の多項式関数に限定される事なく、入力値に対して出力値が一意的に決まる関数であれば、有理関数(例:1/X)、代数関数(例:√X)、あるいは超越関数(例:SIN(X)、LOG(X))を含んでいても良いし、関数形に表さずに入力値、出力値のテーブルとその補間式からなる形で用意してもよい。また、ここで使用する非線形変換は、X線が平均的な物質を透過する場合のビームハードニングを較正するものが良いが、第1近似として後述する軟部組織に適用するビームハードニング補正関数で代用してもよい。 In the figure, the non-linear transformation is illustrated by taking a cubic function passing through the origin as an example, but is not limited to an Nth-order polynomial function, and can be rational as long as the output value is uniquely determined with respect to the input value. Functions (eg 1 / X), algebraic functions (eg √X), transcendental functions (eg SIN (X), LOG (X)) may be included, and input is not shown in the function form It may be prepared in the form of a table of values and output values and its interpolation formula. In addition, the nonlinear transformation used here is preferably one that calibrates beam hardening when X-rays pass through an average substance, but is a beam hardening correction function that is applied to soft tissue, which will be described later, as a first approximation. You may substitute.
(ステップS270)
第1のビームハードニング補正投影データ240に、フィルタリング手段270によるフィルタリング処理を施す(S270)。
(Step S270)
The first beam hardening corrected projection data 240 is subjected to filtering processing by the filtering means 270 (S270).
(ステップS280)
フィルタリング処理後の第1のビームハードニング補正投影データ240に、逆投影手段280に基づいた逆投影演算を行い、第1の再構成CT像291を生成する(S280)。
(Step S280)
A back projection operation based on the back projection unit 280 is performed on the first beam hardening corrected projection data 240 after the filtering process to generate a first reconstructed CT image 291 (S280).
(ステップS311)
再投影係数入力手段310により、上述した図4−1のような入力画面を使用して、またはあらかじめ登録してあるテーブルを読み込んで、メタルのX線吸収係数に相当する物質の再投影データを生成するための再投影係数1を入力する(S311)。
(Step S311)
The reprojection coefficient input means 310 uses the input screen as shown in FIG. 4-1 described above or reads a pre-registered table to obtain the reprojection data of the substance corresponding to the metal X-ray absorption coefficient. A reprojection coefficient 1 for generation is input (S311).
(ステップS312)
再投影係数入力手段310により、上述した図4−2のような入力画面を使用して、またはあらかじめ登録してあるテーブルを読み込んで、骨以上のX線吸収係数の物質に対応する再投影データを生成するための再投影係数2を入力する(S312)。
(Step S312)
Reprojection data corresponding to a substance having an X-ray absorption coefficient equal to or greater than bone is obtained by using the input screen as shown in FIG. Is input as a reprojection coefficient 2 (S312).
(ステップS321)
再投影手段320は、ステップS311で設定した再投影係数1を用いて第1の再構成CT像291の再投影演算を行い、メタル再投影データ323(FpjM)を生成する(S321)。
(Step S321)
The reprojection means 320 performs a reprojection operation on the first reconstructed CT image 291 using the reprojection coefficient 1 set in step S311 and generates metal reprojection data 323 (Fpj M ) (S321).
(ステップS322)
再投影手段320は、ステップS312で設定した再投影係数2を用いて第1の再構成CT像291の再投影演算を行い、骨以上のX線吸収係数の物質に対応する再投影データ324(FpjM+B)を生成する(S322)。
(Step S322)
The reprojection means 320 performs a reprojection calculation of the first reconstructed CT image 291 using the reprojection coefficient 2 set in step S312, and reprojection data 324 ( Fpj M + B ) is generated (S322).
(ステップS330)
再投影手段320は、ステップS322が生成した骨以上再投影データ324(FpjM+B)から、ステップS321で生成されたメタル再投影データ323(FpjM)を引算し、骨成分再投影データ(FpjB)を生成する(S330)。
(Step S330)
The reprojection means 320 subtracts the metal reprojection data 323 (Fpj M ) generated in step S321 from the bone or more reprojection data 324 (Fpj M + B ) generated in step S322, and the bone component reprojection data (Fpj). B ) is generated (S330).
(ステップS331)
投影データの各点において、物質含有率算出手段330が、メタル再投影データ323(FpjM)を第1のビームハードニング補正投影データ240(BHC1)で割算して、メタル含有率(RatioM)を算出(推定)する(S331)。
(Step S331)
At each point of the projection data, the substance content rate calculation means 330 divides the metal reprojection data 323 (Fpj M ) by the first beam hardening corrected projection data 240 (BHC 1 ), and the metal content rate (Ratio). M ) is calculated (estimated) (S331).
(ステップS332)
同様に、物質含有率算出手段330は、骨成分再投影データ(FpjB)を第1のビームハードニング補正投影データ240(BHC1)で割算して、骨含有率(RatioB)を算出する(S332)。
(Step S332)
Similarly, the substance content rate calculation unit 330 divides the bone component reprojection data (Fpj B ) by the first beam hardening corrected projection data 240 (BHC 1 ) to calculate the bone content rate (Ratio B ). (S332).
(ステップS333)
さらに、物質含有率算出手段330が、投影データの各点において、1.0(全物質成分の含有率に相当)から、メタル含有率(RatioM)と、骨含有率(RatioB)を引算し、軟部組織含有率(RatioN)を算出する(S333)。
(Step S333)
Further, the substance content rate calculation means 330 subtracts the metal content rate (Ratio M ) and the bone content rate (Ratio B ) from 1.0 (corresponding to the content rate of all substance components) at each point of the projection data. The soft tissue content (Ratio N ) is calculated (S333).
(ステップS334)
投影データの各点において、物質含有率算出手段330は、元の投影データ220(Prj)にメタル含有率(RatioM)を掛算し、メタル成分投影データ(M)を生成する(S334)。
(Step S334)
In each point of the projection data, substance content rate calculating means 330 multiplies the metal content (Ratio M) based on projection data 220 (Prj), to produce a metal component projection data (M) (S334).
(ステップS335)
同様に、物質含有率算出手段330は、投影データの各点において、投影データ220(Prj)に骨含有率(RatioB)を掛算し、骨成分投影データ(B)を生成する(S335)。
(Step S335)
Similarly, the substance content rate calculation unit 330 multiplies the projection data 220 (Prj) by the bone content rate (Ratio B ) at each point of the projection data to generate bone component projection data (B) (S335).
(ステップS336)
さらに、物質含有率算出手段330は、投影データの各点において、投影データ220(Prj)に軟部組織含有率(RatioN)を掛算し、軟部組織成分投影データ(N)を生成する(S336)。
(Step S336)
Further, the substance content rate calculating means 330 multiplies the projection data 220 (Prj) by the soft tissue content rate (Ratio N ) at each point of the projection data to generate soft tissue component projection data (N) (S336). .
(ステップS341)
物質別ビームハードニング補正手段340は、メタル成分投影データ(M)に非線形変換を施し、メタル成分のビームハードニング補正投影データ344(BHCM)を生成する(S341)。
(Step S341)
The material-specific beam hardening correction means 340 performs nonlinear transformation on the metal component projection data ( M ) to generate metal component beam hardening correction projection data 344 (BHC M ) (S341).
(ステップS342)
同様に、物質別ビームハードニング補正手段340は、骨成分投影データ(B)に非線形変換を施し、骨成分ビームハードニング補正投影データ345(BHCB)を生成する(S342)。
(Step S342)
Similarly, the material-specific beam hardening correction means 340 performs nonlinear transformation on the bone component projection data ( B ) to generate bone component beam hardening correction projection data 345 (BHC B ) (S342).
(ステップS343)
さらに、物質別ビームハードニング補正手段340は、軟部組織成分投影データ(N)に非線形変換を施し、軟部組織成分のビームハードニング補正投影データ346(BHCN)を生成する(S343)。
(Step S343)
Further, the material-specific beam hardening correction means 340 performs nonlinear transformation on the soft tissue component projection data ( N ) to generate soft tissue component beam hardening correction projection data 346 (BHC N ) (S343).
(ステップS350)
物質別ビームハードニング補正手段340は、ステップS341〜ステップS343にて生成した物質別ビームハードニング補正投影データを加算し、高吸収体補正投影データ360(BHC2)を生成する(S350)。
(Step S350)
The material-specific beam hardening correction unit 340 adds the material-specific beam hardening correction projection data generated in steps S341 to S343 to generate high-absorber correction projection data 360 (BHC 2 ) (S350).
以下、生成された高吸収体補正投影データ360(BHC2)をフィルタリング手段270に入力し、フィルタリング処理と逆投影処理を行い、高吸収体補正された再構成CT像290を得ることができる。この再構成CT像290は、画像表示手段120により表示装置80に表示される。 Thereafter, the generated high-absorber corrected projection data 360 (BHC 2 ) is input to the filtering unit 270, and filtering processing and backprojection processing are performed to obtain a reconstructed CT image 290 corrected for high-absorber. The reconstructed CT image 290 is displayed on the display device 80 by the image display means 120.
なお、上記のステップS321及びS322で実施される再投影演算は、共通の幾何学系で実行されるため、同一のサブルーチンの中で演算するようにコードディングする事により、演算量の増大を抑える事ができる。 Note that the reprojection calculation performed in steps S321 and S322 described above is performed in a common geometrical system, so that the amount of calculation is suppressed by coding so as to be calculated in the same subroutine. I can do things.
上記第一の実施形態では、メタルと、骨以上の2つの再投影係数を入力し、元の投影データ220(Prj)を、メタル成分投影データ(M)、骨成分投影データ(B)と、軟部組織成分投影データ(N)の3つに分けて、物質別に異なるビームハードニング補正を施す場合を例に挙げ、その具体的な処理を説明したが、再投影係数の入力が1つの場合であっても高吸収体成分領域と軟部組織とで異なるビームハードニング補正を施こすことである程度の改善があるし、またメタルと、骨以上に加え、第3番目にX線吸収係数が大きい組織(例えば、繊維質)に対応する再投影係数を追加し、元の投影データ220(Prj)を、メタル成分と、骨成分、繊維質成分、軟部組織成分の4成分に分けて、ビームハードニング補正を施せば、さらなる高吸収体成分アーチファクトの低減と、CT値定量性の向上が見込める。再投影係数の数が異なる場合についてもその演算方法は、上述の説明により容易に想像できるため、その詳細は省略する。 In the first embodiment, metal and two reprojection coefficients greater than or equal to bone are input, and the original projection data 220 (Prj) is converted into metal component projection data (M), bone component projection data (B), The example of the case where the softening component projection data (N) is divided into three parts and different beam hardening corrections are performed for each material is described as an example, and the specific processing has been described. However, when the input of the reprojection coefficient is one Even with this, there is a certain degree of improvement by applying different beam hardening corrections in the superabsorbent component region and soft tissue, and in addition to metal and bone, the tissue with the third largest X-ray absorption coefficient Re-projection coefficient corresponding to (for example, fiber) is added, and the original projection data 220 (Prj) is divided into a metal component, a bone component, a fiber component, and a soft tissue component, and beam hardening is performed. If you make corrections, That reduction of the superabsorbent component artifact, the improvement of CT values quantitative property can be expected. Even when the number of reprojection coefficients is different, the calculation method can be easily imagined from the above description, and the details are omitted.
また、物質別ビームハードニング補正手段340において、物質別ビームハードニング補正投影データの加算を行わず、ステップS341〜ステップS343にて生成した物質別のビームハードニング補正投影データを、個々にフィルタリング手段270に入力し、フィルタリング処理と逆投影処理を行い、物質別の再構成CT像を出力するようにしてもよい。 Further, the material-specific beam hardening correction projection data 340 does not add the material-specific beam hardening correction projection data, and individually filters the material-specific beam hardening correction projection data generated in steps S341 to S343. The image may be input to 270, filtering processing and back projection processing may be performed, and a reconstructed CT image for each substance may be output.
最後に、図6−1乃至図6−4に、本発明を実際に適用した再構成CT像を示す。図6−1は、高吸収体補正前のコーンビーム再構成CT像の例であり、画像上半分にある2つの丸い断面は金属であり、画像下半分の2つの長方形断面はセラミックスである。金属と金属、金属とセラミックスを結ぶ間に黒いアーチファクト領域があり、その接線に白いアーチファクトを引いている。 Finally, FIGS. 6-1 to 6-4 show reconstructed CT images to which the present invention is actually applied. FIG. 6A is an example of a cone beam reconstructed CT image before superabsorber correction, in which two round cross sections in the upper half of the image are metal, and two rectangular cross sections in the lower half of the image are ceramics. There is a black artifact area between metal and metal, and between metal and ceramics, and white artifacts are drawn at the tangent line.
図6−2は、従来の方法で、金属を再投影する1つの再投影閾値を設定して、再投影データを生成し補正演算を行った結果を示しており、アーチファクトの低減が不十分なものとなっている。 FIG. 6B shows a result of performing a recalculation data by setting one reprojection threshold value for reprojecting metal by a conventional method, and artifact reduction is insufficient. It has become a thing.
図6−3は、従来の方法で、セラミックス以上のX線吸収係数を有する吸収体を再投影する、低い再投影閾値を設定して、再投影データを生成し補正演算を行った場合を示す。2つの金属を結ぶ方向のアーチファクトが信号として認識され、逆に偽像として強調される結果となった。 FIG. 6-3 shows a case where a reprojection data is generated and a correction operation is performed by setting a low reprojection threshold value by re-projecting an absorber having an X-ray absorption coefficient higher than that of ceramics by a conventional method. . Artifacts in the direction connecting two metals were recognized as signals, and conversely, they were emphasized as false images.
図6−4に、本実施形態に係る高吸収体補正を適用した再構成CT像の例を示す。2つのX線吸収係数の異なる高吸収体成分から生じるアーチファクトが低減されCT値の一様性が向上していること事が分かる。 FIG. 6-4 shows an example of a reconstructed CT image to which the high absorber correction according to the present embodiment is applied. It can be seen that artifacts caused by two high-absorber components having different X-ray absorption coefficients are reduced and the uniformity of the CT value is improved.
以上、本発明の実施形態を用いて説明してきたが、これまでの各実施形態で説明した構成はあくまで一例であり、例えば、一度、高吸収体補正した再構成CT像を再び再投影し、高吸収体成分アーチファクト低減精度、CT値定量性精度を向上させるなど、本発明は、技術的思想を逸脱しない範囲内で適宜変更が可能である。 As described above, the embodiments of the present invention have been described. However, the configuration described in each of the embodiments so far is merely an example.For example, the reconstructed CT image corrected once with the high absorber is re-projected, The present invention can be modified as appropriate within a range that does not depart from the technical idea, such as improving the high-absorber component artifact reduction accuracy and CT value quantification accuracy.
本発明によれば、高吸収体成分アーチファクトの低減に要する演算時間を短くし、かつCT値定量性に優れた高度な高吸収体成分アーチファクト低減機能を有するコーンビームX線CT装置を提供することができ、頭部、腹部等の造影撮影、並びに歯顎、腰椎、四肢の整形外科撮影の診断性能を向上させることが期待できる。 According to the present invention, it is possible to provide a cone beam X-ray CT apparatus having an advanced high-absorber component artifact reduction function that shortens the calculation time required for reducing the high-absorber component artifact and is excellent in CT value quantification. It can be expected to improve the diagnostic performance of contrast imaging of the head, abdomen, etc. and orthopedic imaging of the teeth, jaws, lumbar vertebrae, and limbs.
1…コーンビームX線CT装置、1a…座位方式コーンビームX線CT装置、1b…被検体回転方式コーンビームX線CT装置、1c…移動型X線装置に搭載されたCアーム方式コーンビームX線CT装置、2…被検体、3…回転軌道面(ミッドプレーン)、4…回転中心軸、5…旋回アーム、6…支柱、7…椅子、8…固定架台、9…車輪、10…撮影部、10a…座位方式コーンビームX線CT1aの撮影部、10c…移動型X線装置に搭載されたCアーム方式コーンビームX線CT装置1cの撮影部、11…X線源、11t…X線管、11c…コリメータ、12…2次元X線検出器、12i…X線イメージインテンシファイア、12c…CCDテレビカメラ、13…C型アーム、14…C型アーム保持体、15…天井支持体、16…天井レール、17…寝台、18…インジェクタ、19…回転台、20…制御演算部、20a…座位方式コーンビームX線CT1aの制御演算部、21…回転中心軸4の2次元X線検出器12への投影、22…高吸収体成分を含む断面、22c…高吸収体成分を含む断面、23…歯列、24…金属、25…寛骨臼、26…大転子、27…人工股関節、29…アーチファクト、30…高吸収体成分を含む断面の再構成CT像、30c…高吸収体成分を含む断面の再構成CT像、31…再構成CT像30のメタル閾値表示、31c…再構成CT像30cのメタル閾値表示、32…再構成CT像30の骨閾値表示、32c…再構成CT像30cの骨閾値表示、40…再投影係数入力画面、41…再投影係数1選択タブ、42…再投影係数2選択タブ、43…再投影係数名表示欄、44…再投影係数追加ボタン、50…再投影下限値設定ポイント、51…再投影係数設定ポイント、52…再投影上限値設定ポイント、53…再投影閾値設定欄、54…再投影係数設定欄、55…再投影係数設定ポイント追加ボタン、56…再投影係数1設定カーブ、57…再投影係数2設定カーブ、70…情報入力装置、80…表示装置、100…撮影部制御手段、100a…座位方式コーンビームX線CT装置1aの撮影部制御手段、100b…被検体回転方式コーンビームX線CT装置1bの撮影部制御手段、100c…移動型X線装置に搭載されたCアーム方式コーンビームX線CT装置1cの撮影部制御手段、101…撮影系回転制御手段、102…撮影系位置制御手段、103…X線照射制御手段、104…インジェクタ制御手段、105…椅子位置制御手段、106…寝台制御手段、107…検出系制御手段、110…画像収集手段、120…画像表示手段、200…再構成手段、210…前処理手段、220…投影データ、230…第1のビームハードニング補正手段、240…第1のビームハードニング補正投影データ、270…フィルタリング手段、280…逆投影手段、290…再構成CT像、291…第1の再構成CT像、300…高吸収体補正手段、310…再投影係数入力手段、320…再投影手段、323…メタル再投影データ、324…骨以上再投影データ、330…物質含有率算出手段、340…物質別ビームハードニング補正手段、344…メタル成分ビームハードニング補正投影データ、345…骨成分ビームハードニング補正投影データ、344…軟部組織ビームハードニング補正投影データ、360…高吸収体補正投影データ DESCRIPTION OF SYMBOLS 1 ... Cone beam X-ray CT apparatus, 1a ... Sitting type cone beam X-ray CT apparatus, 1b ... Subject rotation type cone beam X-ray CT apparatus, 1c ... C-arm type cone beam X mounted in a mobile X-ray apparatus Line CT apparatus, 2 ... Subject, 3 ... Rotating orbital plane (midplane), 4 ... Rotation center axis, 5 ... Swivel arm, 6 ... Pole, 7 ... Chair, 8 ... Fixed mount, 9 ... Wheel, 10 ... Imaging , 10a... Sitting type cone beam X-ray CT1a imaging unit, 10c... C-arm type cone beam X-ray CT apparatus 1c imaging unit mounted on mobile X-ray apparatus, 11... X-ray source, 11t. Tube, 11c ... collimator, 12 ... two-dimensional X-ray detector, 12i ... X-ray image intensifier, 12c ... CCD TV camera, 13 ... C-type arm, 14 ... C-type arm holder, 15 ... ceiling support, 16 ... Heaven Rail, 17 ... bed, 18 ... injector, 19 ... turntable, 20 ... control calculation unit, 20a ... control calculation unit of sitting-type cone beam X-ray CT1a, 21 ... to the two-dimensional X-ray detector 12 of the rotation center axis 4 , 22 ... cross section including superabsorbent component, 22c ... cross section including superabsorbent component, 23 ... dentition, 24 ... metal, 25 ... acetabulum, 26 ... greater trochanter, 27 ... artificial hip joint, 29 ... artifact, 30 ... reconstructed CT image of a cross section including a high absorber component, 30c ... reconstructed CT image of a cross section including a high absorber component, 31 ... metal threshold value display of the reconstructed CT image 30, 31c ... reconstructed CT Metal threshold display of image 30c, 32 ... Bone threshold display of reconstructed CT image 30, 32c ... Bone threshold display of reconstructed CT image 30c, 40 ... Reprojection coefficient input screen, 41 ... Reprojection coefficient 1 selection tab, 42 ... Reprojection factor 2 selection tab, 4 ... Reprojection coefficient name display field, 44 ... Reprojection coefficient addition button, 50 ... Reprojection lower limit value setting point, 51 ... Reprojection coefficient setting point, 52 ... Reprojection upper limit value setting point, 53 ... Reprojection threshold value setting field, 54 ... Reprojection coefficient setting field, 55 ... Reprojection coefficient setting point addition button, 56 ... Reprojection coefficient 1 setting curve, 57 ... Reprojection coefficient 2 setting curve, 70 ... Information input device, 80 ... Display device, 100 ... Shooting Mounted on the control unit, 100a ... imaging unit control unit of the sitting-type cone beam X-ray CT apparatus 1a, 100b ... imaging unit control unit of the subject rotation type cone-beam X-ray CT apparatus 1b, 100c ... mounted on the mobile X-ray apparatus Imaging unit control means of the C-arm type cone-beam X-ray CT apparatus 1c, 101 ... Imaging system rotation control means, 102 ... Imaging system position control means, 103 ... X-ray irradiation control means, 10 4 ... injector control means, 105 ... chair position control means, 106 ... bed control means, 107 ... detection system control means, 110 ... image collection means, 120 ... image display means, 200 ... reconstruction means, 210 ... preprocessing means, 220 ... projection data, 230 ... first beam hardening correction means, 240 ... first beam hardening correction projection data, 270 ... filtering means, 280 ... back projection means, 290 ... reconstructed CT image, 291 ... first Reconstructed CT image, 300 ... high absorber correction means, 310 ... reprojection coefficient input means, 320 ... reprojection means, 323 ... metal reprojection data, 324 ... reprojection data over bone, 330 ... substance content calculation means 340 ... Beam hardening correction means for each material, 344 ... Metal component beam hardening correction projection data, 345 ... Bone component beam heart Training correction projection data, 344 ... soft tissue beam hardening correction projection data, 360 ... high-absorber correction projection data
Claims (2)
前記X線源に対向して配置され、高吸収体成分を含む被検体を透過した前記X線を検出して前記被検体の投影データを出力するX線検出器と、
前記X線源と前記X線検出器とを回転移動させる回転手段と、
前記被検体の投影データに基づいて再構成演算を行い、前記被検体のX線CT像を生成する再構成手段と、
前記再構成手段が生成した画像を表示する表示手段と、
前記X線CT像を前記高吸収体成分に対応させた再投影係数に基づいて再投影演算をすることにより前記高吸収体成分の再投影データを生成し、その高吸収体成分の再投影データと前記被検体の投影データとに基づいて、前記高吸収体成分の投影データと前記被検体に含まれる前記高吸収体成分を除く物質である軟部組織成分の投影データとを生成し、それら前記高吸収体成分の投影データと前記軟部組織成分の投影データとの各々に対して前記X線の線質硬化に起因するアーチファクトを低減するための補正を施す高吸収体補正手段と、を備え、
前記再構成手段は、前記補正後の前記高吸収体成分の投影データ及び前記軟部組織成分の投影データに基づいて再構成演算を行い、前記アーチファクトを補正したX線CT像を生成し、
前記表示手段は、前記アーチファクトを補正したX線CT像を表示し、
前記再投影係数は、CT値に応じて0から1の間で連続的に変化するものであり、
前記高吸収体成分の投影データは、前記被検体の投影データと前記投影データに含まれる前記高吸収体成分の割合を示す高吸収体成分含有率とを乗算することにより算出され、
前記高吸収体成分含有率は、前記高吸収体成分の再投影データを前記被検体の投影データで除算することにより算出され、
前記軟部組織成分の投影データは、前記被検体の投影データと前記被検体の投影データに含まれる前記軟部組織成分の割合を示す軟部組織成分含有率とを乗算することにより算出され、
前記軟部組織成分含有率は、全物質成分の含有率から前記高吸収体成分含有率を減算することにより算出される、
ことを特徴とするX線CT装置。 An X-ray source for generating X-rays;
An X-ray detector that is disposed opposite to the X-ray source, detects the X-ray transmitted through the subject including a high-absorber component, and outputs projection data of the subject;
Rotating means for rotating the X-ray source and the X-ray detector;
Reconstruction means for performing reconstruction calculation based on the projection data of the subject and generating an X-ray CT image of the subject;
Display means for displaying the image generated by the reconstruction means;
Reprojection data of the high absorber component is generated by performing a reprojection operation based on a reprojection coefficient corresponding to the high absorber component of the X-ray CT image, and the reprojection data of the high absorber component And the projection data of the subject to generate projection data of the high-absorber component and projection data of the soft tissue component that is a substance excluding the high-absorber component contained in the subject, High-absorber correction means for performing correction for reducing artifacts caused by X-ray radiation hardening on each of the projection data of the high-absorber component and the projection data of the soft tissue component,
The reconstruction means performs a reconstruction operation based on the corrected projection data of the superabsorber component and the projection data of the soft tissue component, and generates an X-ray CT image in which the artifact is corrected,
The display means displays an X-ray CT image in which the artifact is corrected,
The reprojection coefficient changes continuously between 0 and 1 according to the CT value,
The projection data of the high-absorber component is calculated by multiplying the projection data of the subject and the high-absorber component content indicating the ratio of the high-absorber component included in the projection data,
The high-absorber component content is calculated by dividing the re-projection data of the high-absorber component by the projection data of the subject,
The projection data of the soft tissue component is calculated by multiplying the projection data of the subject and the soft tissue component content indicating the ratio of the soft tissue component included in the projection data of the subject,
The soft tissue component content is calculated by subtracting the high absorbent component content from the content of all substance components,
An X-ray CT apparatus characterized by that.
前記表示手段は、入力された再投影係数に基づいて生成される再投影データの高吸収体成分を示すX線CT像を表示する、
ことを特徴とする請求項1記載のX線CT装置。 Re-projection coefficient input means for inputting an arbitrary number of the re-projection coefficients corresponding to the number of superabsorber components;
The display means displays an X-ray CT image indicating a high-absorber component of reprojection data generated based on the input reprojection coefficient.
The X-ray CT apparatus according to claim 1 .
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008235985A JP5243160B2 (en) | 2008-09-16 | 2008-09-16 | X-ray CT system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008235985A JP5243160B2 (en) | 2008-09-16 | 2008-09-16 | X-ray CT system |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2010068832A JP2010068832A (en) | 2010-04-02 |
JP5243160B2 true JP5243160B2 (en) | 2013-07-24 |
Family
ID=42201184
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008235985A Expired - Fee Related JP5243160B2 (en) | 2008-09-16 | 2008-09-16 | X-ray CT system |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5243160B2 (en) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101531440B1 (en) * | 2010-04-20 | 2015-06-24 | 이미징 사이언시즈 인터내셔널 엘엘씨 | Reduction and removal of artifacts from a three-dimensional dental x-ray data set using surface scan information |
US9684952B2 (en) | 2010-04-20 | 2017-06-20 | Dental Imaging Technologies Corporation | Alignment of mixed-modality data sets for reduction and removal of imaging artifacts |
IN2014CN02779A (en) * | 2011-10-24 | 2015-07-03 | Koninkl Philips Nv | |
JP6142172B2 (en) * | 2014-06-03 | 2017-06-07 | 朝日レントゲン工業株式会社 | Beam hardening correction apparatus, beam hardening correction method, and X-ray imaging apparatus |
JP7258474B2 (en) * | 2018-05-10 | 2023-04-17 | 株式会社日立製作所 | X-ray CT device and radiotherapy system |
CN110840479B (en) * | 2019-12-03 | 2023-05-26 | 上海联影医疗科技股份有限公司 | Bone component determination method, device, equipment and storage medium |
CN113100803B (en) * | 2021-04-20 | 2024-07-19 | 西门子数字医疗科技(上海)有限公司 | Method, apparatus, computer device and medium for displaying venous thrombosis |
DE112022003708T5 (en) | 2021-07-27 | 2024-05-16 | Rigaku Corporation | Correction device, system, procedure and program |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3223195B2 (en) * | 1992-01-23 | 2001-10-29 | 株式会社日立メディコ | X-ray CT system |
JP3414471B2 (en) * | 1994-01-24 | 2003-06-09 | 株式会社東芝 | X-ray CT system |
JP3556292B2 (en) * | 1994-10-28 | 2004-08-18 | 株式会社東芝 | Image processing device |
JP4854137B2 (en) * | 2001-06-21 | 2012-01-18 | 株式会社東芝 | Medical diagnostic imaging equipment |
JP4584550B2 (en) * | 2003-06-05 | 2010-11-24 | 株式会社日立メディコ | X-ray measuring device |
JP5107568B2 (en) * | 2006-12-12 | 2012-12-26 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X-ray CT system |
-
2008
- 2008-09-16 JP JP2008235985A patent/JP5243160B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2010068832A (en) | 2010-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5243160B2 (en) | X-ray CT system | |
JP5248031B2 (en) | X-ray CT system | |
Scarfe et al. | What is cone-beam CT and how does it work? | |
JP4854137B2 (en) | Medical diagnostic imaging equipment | |
JP4350738B2 (en) | X-ray tomography apparatus and artifact reduction method | |
JP4414420B2 (en) | X-ray tomography apparatus and artifact reduction method | |
JP2013099700A (en) | X-ray computed tomography apparatus, reconstruction processing apparatus, and image processing apparatus | |
Sheth et al. | Mobile C‐Arm with a CMOS detector: Technical assessment of fluoroscopy and Cone‐Beam CT imaging performance | |
WO2005011502A1 (en) | Tomograph | |
JP4989473B2 (en) | Three-dimensional reproduction method and apparatus using an inclined configuration | |
JP3763967B2 (en) | X-ray equipment | |
JP6566887B2 (en) | Radiation image processing apparatus, method and program | |
JP5487172B2 (en) | Medical image diagnostic apparatus and medical image processing apparatus | |
US20220071578A1 (en) | Improved method of acquiring a radiographic scan of a region-of-interest in a metal containing object | |
JP5808672B2 (en) | X-ray CT apparatus and X-ray CT imaging method | |
JP3742193B2 (en) | X-ray computed tomography system | |
JP5303154B2 (en) | X-ray CT system | |
JP4610304B2 (en) | X-ray CT system | |
JP4703221B2 (en) | X-ray CT system | |
JP5433283B2 (en) | X-ray computed tomography system | |
JP2010075553A (en) | X-ray computed tomography system | |
JP2000083946A (en) | Method and device for correction projection and radiation tomography apparatus | |
JP5031335B2 (en) | X-ray CT system | |
JP5854658B2 (en) | X-ray CT system | |
JP2019187814A (en) | Medical image processing apparatus, x-ray diagnostic apparatus and medical image processing program |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20110831 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20121120 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130115 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130212 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130304 |
|
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: 20130326 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20130404 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20160412 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
LAPS | Cancellation because of no payment of annual fees |