JP2004061256A - Computerized tomography - Google Patents

Computerized tomography Download PDF

Info

Publication number
JP2004061256A
JP2004061256A JP2002219099A JP2002219099A JP2004061256A JP 2004061256 A JP2004061256 A JP 2004061256A JP 2002219099 A JP2002219099 A JP 2002219099A JP 2002219099 A JP2002219099 A JP 2002219099A JP 2004061256 A JP2004061256 A JP 2004061256A
Authority
JP
Japan
Prior art keywords
center axis
rotation center
projection position
ray
data
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.)
Granted
Application number
JP2002219099A
Other languages
Japanese (ja)
Other versions
JP3846577B2 (en
Inventor
Taketo Kishi
岸 武人
Shuhei Onishi
大西 修平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shimadzu Corp
Original Assignee
Shimadzu Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shimadzu Corp filed Critical Shimadzu Corp
Priority to JP2002219099A priority Critical patent/JP3846577B2/en
Publication of JP2004061256A publication Critical patent/JP2004061256A/en
Application granted granted Critical
Publication of JP3846577B2 publication Critical patent/JP3846577B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Abstract

<P>PROBLEM TO BE SOLVED: To provide a computerized tomography for accurately estimating a projection position onto an X-ray detector at a rotary center axis and greatly contributing to the improvement of a tomographic image by placing an arbitrary object on a turntable , rotating the object, and sampling X-ray transmission data without using any exclusive phantoms. <P>SOLUTION: The computerized tomography has a rotary center axis projection position estimation means 53 that changes a rotary center axis projection position used for converting X-ray transmission data rows detected by the X-ray detector 2 in a plurality of times to parallel projection data rows, and estimates the rotary center axis projection position based on each data row being converted by using each rotary center axis projection position. The tomogram is reconfigured by the data rows being converted by using the estimated rotary center axis projection position, thus obtaining the stable, high-quality sectional image without being affected by manufacturing errors in the phantoms. <P>COPYRIGHT: (C)2004,JPO

Description

【0001】
【発明の属する技術分野】
本発明はコンピュータ断層撮影装置に関し、更に詳しくは、被写体をX線源とX線検出器の間で回転させる方式のコンピュータ断層撮影装置に関する。
【0002】
【従来の技術】
例えば産業用の被破壊検査などに用いられるX線CT装置、つまりX線を用いたコンピュータ断層撮影装置においては、一般に、ファンビーム状のX線を発生するX線源と、そのX線源に対向して、X線の広がり方向に複数の素子が配列されたX線検出器を配置し、これらの間に、被写体を搭載してX線の光軸方向とX線検出器の素子の配列方向の双方に直交する回転軸の回りに回転させるターンテーブルを配置した構成が採用される。
【0003】
そして、被写体を回転させつつ、所定の微小角度ごとに被写体を透過したX線強度の検出データを360°にわたってX線検出器の各素子出力から採取し、その各回転角度ごとの各素子による透過X線検出データをもとに、回転軸に直交する平面に沿った被写体の断層像を再構成する。
【0004】
ここで、X線検出器には、各素子ごとの感度むらや歪みが存在し、また、X線検出器に入射する透過X線は、点状のX線源から出たファンビーム状のX線が被写体を透過したものであるため、X線検出器の各素子出力をそのまま用いては断層像を再構成することはできない。
【0005】
そこで、従来、各素子出力を感度補正および歪み補正を行うとともに、ファン−パラ変換と称される変換によって、X線検出器の各素子による透過X線データ列を、平行投影データ列に変換した後のデータ列を用いて、断層像の再構成演算を行う。このデータの変換には、被写体の回転中心軸のX線検出器上への投影位置(以下、回転中心軸投影位置と称する)が正確に判っていなければ、正確なデータ変換を行うことができず、ひいては高品質の断層像を得ることができない。
【0006】
回転中心軸投影位置を推定する手法として、従来、例えばアクリルパイプ内にタングステンなどのX線吸収率の高い材質からなるワイヤを通す等の構造を有する専用のファントムを用意し、そのファントムをターンテーブル上に搭載して回転を与え、そのファントムを透過したX線をX線検出器で検出して得た360°分のデータ列より、ワイヤが最も左に寄った点と右に寄った点を検出し、それら2点の中心点を回転中心軸投影位置として推定していた。
【0007】
【発明が解決しようとする課題】
しかしながら、上記した従来の回転中心投影位置の推定手法によると、専用のファントムの製造誤差が、そのまま回転中心軸投影位置の推定結果に影響を及ぼし、その出来具合に応じて断層像の品質が左右されるという問題があった。
【0008】
本発明はこのような実情に鑑みてなされたもので、専用のファントム等を用いることなく、任意の物体をターンテーブル上に搭載して回転を与えてX線透過データを採取することにより、回転中心軸投影位置を正確に推定することができ、コンピュータが断層画像の高品質化に大きく寄与することのできるコンピュータ断層撮影装置の提供を目的としている。
【0009】
【課題を解決するための手段】
上記の目的を達成するため、本発明のコンピュータ断層撮影装置は、ファン状に広がりを持つX線を出力するX線源と、そのX線源に対向配置され、当該X線源からのX線光軸に直交し、かつ、X線の広がり方向に複数の素子が配列されたX線検出器と、そのX線源とX線検出器の間に配置され、被写体を搭載して上記X線光軸および素子の配列方向の双方に直交する回転中心軸の回りを回転するターンテーブルと、上記ターンテーブルを駆動しつつ被写体にX線を照射して所定角度ごとに検出した透過X線データ列を、回転中心軸のX線検出器上への投影位置を用いて平行投影データ列に変換するデータ変換手段と、その変換後のデータを用いて、上記回転中心軸に直交する面で被写体をスライスした断層像生成する再構成演算手段を備えたコンピュータ断層撮影装置において、上記データ変換手段で用いる回転中心軸投影位置を複数にわたって変化させ、その各回転中心軸投影位置を用いて変換した各データ列に基づいて回転中心軸投影位置を推定する回転中心軸投影位置推定手段を備え、上記再構成演算手段は、その推定手段により推定された回転中心軸投影位置を用いて変換した平行投影データから断層像を構成することによって特徴づけられる(請求項1)。
【0010】
ここで、本発明における回転中心軸投影位置推定手段による最適の回転中心軸投影位置の決定手法の具体的で好適な手法としては、請求項2に係る発明のように、複数に変化させた回転中心軸投影位置のそれぞれを用いて変換した複数の平行投影データ列による各サイノグラムについて、180°対向するデータどうしの対称性を評価し、対称性が最も高くなる平行投影データ列を得た回転中心軸投影位置を最適と推定する手法を好適に採用することができる。
【0011】
本発明は、被写体にファンビーム状のX線を照射してX線検出器により検出した被写体のX線透過データ列を、X線検出器の感度補正および歪み補正を行うとともに平行投影データ列に変換する際、回転中心軸投影位置を種々に変化させ、その各回転中心軸投影位置を用いて変換した平行投影データ列から、真の回転中心軸投影位置またはそれに最も近い回転中心軸投影位置を推定することにより、専用のファントムを用いることなく、任意の被写体を用いた透過X線データから、正確に回転中心軸投影位置を求めるものである。
【0012】
具体的には、ファンビーム状のX線が被写体を透過したX線の検出データ列から、平行投影データ列に変換する際に用いる回転中心軸投影位置として、真の回転中心軸投影位置を用いた場合、変換後のデータ列においては、例えばそのサイノグラム中で180°対向するデータどうしが対称となるはずである。そこで、任意の被写体を360°回転させて得た透過X線のデータ列を、平行投影データ列に変換する際に用いる回転中心軸投影位置を種々に変化させ、それぞれの回転中心軸投影位置を用いて変換したデータ列について、それぞれのサイノグラムにおいて180°対向するデータどうしの対称性を評価することによって、最も対称性の高いデータ列が得られた回転中心軸投影位置が、真の回転中心軸投影位置もしくはそれに最も近い回転中心軸投影位置であると推定することができる。そして、このように推定された回転中心軸投影位置を用いて変換したデータ列を用いて断層像を再構成することにより、得られる断層像は常に高品質なものとなる。
【0013】
【発明の実施の形態】
以下、図面を参照しつつ本発明の実施の形態について説明する。
図1は本発明の実施の形態の構成図で、要部の機械的構成を表す模式図と、電気的な機能構成を表すブロック図とを併記して示す図である。
【0014】
X線源1で発生したX線はコリメータ11を介してファン状のビームBとなり、対向配置されているX線検出器2側に向かう方向(x方向)に照射される。X線検出器2は、X線の広がり方向(y方向)に複数の素子が円弧状に配列された構造を有している。
【0015】
X線源1とX線検出器2の間には、X線の光軸方向であるx方向とX線検出器2の素子の配列方向であるy方向の双方に直交する方向(z方向)に沿った回転中心軸Raの回りを回転するターンテーブル3が配置されており、被写体Wはこのターンテーブル3上に搭載されて回転が与えられる。なお、ターンテーブル3は、x,y,z方向に移動可能なx−y−zテーブル31上に載せられている。
【0016】
X線検出器2の各素子の出力は、ターンテーブル3により被写体Wが微小角度回転するごとに、データサンプリング回路4によりデジタル化されたうえで、演算装置5に取り込まれる。すなわち、演算装置5には、X線検出器2の各素子ごとに、被写体Wの360°分のX線透過データが取り込まれる。
【0017】
演算装置5は、データサンプリング回路4からのX線透過データのガンマ補正、画像歪み補正、更には各素子の感度補正や、対数変換などを行う前処理部51と、その前処理後のX線透過データ列を、ファン−パラ変換して平行投影データ列に変換するデータ変換部52と、その変換に必要な、回転中心軸RaのX線検出器2上への投影位置、つまり回転中心軸投影位置を推定する回転中心軸投影位置推定部53と、データ変換部52による変換後のデータ列を用いて断層像を生成する断層像再構成演算部54を備えている。そして、この断層像再構成演算部54により生成された被写体Wの断層像は表示器6に表示される。
【0018】
回転中心軸投影位置推定部53では、データサンプリング回路4により取り込んだX線検出器2の各素子からの360°分のX線透過データを用いて、以下の手法により回転中心軸投影位置を推定する。
【0019】
すなわち、回転中心軸投影位置を複数にわたって変化させ、その各回転中心軸投影位置を用いて、取り込んだデータ列を平行投影データ列に変換する。そして、変換後の各データ列によるサイノグラムについて、180°対向するデータどうしを比較し、その対称性を評価する。そして、最も対称性の高いデータ列を得た回転中心軸投影位置を、真の、あるいは最も真の位置に近い回転中心軸投影位置として推定する。
【0020】
図2にサイノグラムを例示する。この例では、横軸を変換後のチャンネルナンバー(y方向位置)、縦軸を回転角度φとして、360°分のデータを並べている。このようなサイノグラム中における180°対向するデータどうしの対称性の評価の具体的方法については、幾つかの方法を挙げることができ、それぞれの方法について説明する。
【0021】
一つ目は、各サイノグラム中で、回転角度φ°のデータ列と(φ+180)°のデータ列において、それぞれに対応する(図2に矢印で例示するように、y方向位置中心を挟んで互いに対称の位置に存在する)チャンネルのデータの値(X線透過量、換言すれば明るさ)どうしの差分の分散値もしくは標準偏差値を求める。そして、この分散値もしくは標準偏差値が最も小さくなるサイノグラムを形成しているデータ列を算出するのに用いた回転中心軸投影位置を真の投影位置、あるいは最も真に近い投影位置と推定する。
【0022】
二つ目は、サイノグラム中で回転角度φ°のデータ列と(φ+180)°のデータ列において、上記と同様、それぞれに対応するチャンネルのデータの値の加算値の最大値と最小値の差を計算し、その値が最も大きくなるサイノグラムを形成しているデータ列を算出するのに用いた回転中心軸投影位置を真の投影位置、あるいは最も真に近い投影位置と推定する。
【0023】
三つ目は、各サイノグラム中におけるφ°のデータ列と(φ+180)°のデータ列のピーク検出処理を行い、ピーク位置およびそのピーク位置におけるデータ値(X線透過量)を求める。そして、φ°のデータ列におけるピーク位置並びにそのデータ値と、(φ+180)°のデータ列におけるピーク位置並びにそのデータ値を比較し、互いのピーク位置とそのデータ値が最も一致するサイノグラムを形成しているデータ列を算出するのに用いた回転中心軸投影位置を真の投影位置、あるいは最も真に近い投影位置と推定する。
【0024】
なお、以上の各方法において、φは任意の単一の角度であってもよいし、φの値を0〜180°まで変化させて、それぞれの評価結果の加算もしくは平均処理を行って最終的な評価を行ってもよい。
【0025】
再構成演算部54では、以上のようにして推定された回転中心軸投影位置を用いて変換された平行投影データ列を用いて被写体Wの断層像再構成する。
従って、以上の実施の形態によって得られる断層像は、常に正確な回転中心軸を用いて変換された平行投影データ列を用いて生成されるので、常に高い品質の像となる。
【0026】
【発明の効果】
以上のように、本発明によれば、任意の被写体のX線透過像データから回転中心軸投影位置を正確に推定するので、従来のように専用のファントムを用いる必要がなく、しかもそのファントムの製造誤差に起因する回転中心軸投影位置の推定誤差を生じる恐れがなくなり、常に安定して正確に回転中心軸投影位置を推定することが可能となって、得られる断層像は常に安定して高品質なものとなる。
【図面の簡単な説明】
【図1】本発明の実施の形態の構成図で、要部の機械的構成を表す模式図と電気的な機能構成を表すブロック図とを併記して示す図である。
【図2】本発明の実施の形態における回転中心軸投影位置推定部53で行う回転中心軸の投影位置の推定処理に用いるサイノグラムの例を示す図である。
【符号の説明】
1 X線源
2 X線検出器
3 ターンテーブル
4 データサンプリング回路
5 演算装置
51 前処理部
52 データ変換部
53 回転中心軸投影位置推定部
54 断層像再構成演算部
6 表示器
Ra 回転中心軸
W 被写体
[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to a computer tomography apparatus, and more particularly, to a computer tomography apparatus that rotates an object between an X-ray source and an X-ray detector.
[0002]
[Prior art]
For example, in an X-ray CT apparatus used for industrial destructive inspection, that is, a computer tomography apparatus using X-rays, generally, an X-ray source that generates X-rays in the form of a fan beam, and an X-ray source An X-ray detector in which a plurality of elements are arranged in the direction in which the X-ray spreads is arranged opposite to the object, and an object is mounted between the X-ray detectors and the direction of the optical axis of the X-ray and the arrangement of the elements of the X-ray detector. A configuration in which a turntable that rotates around a rotation axis orthogonal to both directions is employed.
[0003]
Then, while rotating the subject, the detection data of the X-ray intensity transmitted through the subject at every predetermined minute angle is collected from the output of each element of the X-ray detector over 360 °, and the transmission by each element at each rotation angle is obtained. Based on the X-ray detection data, a tomographic image of the subject along a plane orthogonal to the rotation axis is reconstructed.
[0004]
Here, the X-ray detector has sensitivity unevenness and distortion for each element, and transmitted X-rays incident on the X-ray detector are fan beam-like X-rays emitted from a point-like X-ray source. Since the line is transmitted through the subject, it is not possible to reconstruct a tomographic image using the element outputs of the X-ray detector as they are.
[0005]
Therefore, conventionally, sensitivity correction and distortion correction have been performed on the output of each element, and a transmission X-ray data string by each element of the X-ray detector has been converted into a parallel projection data string by a conversion called fan-para conversion. Reconstruction calculation of a tomographic image is performed using the subsequent data sequence. In this data conversion, accurate data conversion can be performed unless the projection position of the rotation center axis of the subject on the X-ray detector (hereinafter referred to as the rotation center axis projection position) is accurately known. Therefore, a high quality tomographic image cannot be obtained.
[0006]
Conventionally, as a method of estimating the rotation center axis projection position, a dedicated phantom having a structure such as passing a wire made of a material having a high X-ray absorptivity such as tungsten in an acrylic pipe is prepared, and the phantom is turned on a turntable. From the 360 ° data train obtained by detecting the X-rays transmitted through the phantom with the X-ray detector mounted on the top and rotating it, the points where the wires were shifted to the left and right were determined. The center point of these two points is detected and estimated as the rotation center axis projection position.
[0007]
[Problems to be solved by the invention]
However, according to the above-described conventional method of estimating the rotation center projection position, the manufacturing error of the dedicated phantom directly affects the estimation result of the rotation center axis projection position, and the quality of the tomographic image depends on the degree of performance. There was a problem that was.
[0008]
The present invention has been made in view of such circumstances, and without using a dedicated phantom or the like, an arbitrary object is mounted on a turntable and rotated to collect X-ray transmission data, thereby obtaining a rotation. It is an object of the present invention to provide a computer tomography apparatus capable of accurately estimating a central axis projection position and enabling a computer to greatly contribute to high quality of a tomographic image.
[0009]
[Means for Solving the Problems]
In order to achieve the above object, a computer tomography apparatus according to the present invention includes an X-ray source that outputs an X-ray having a fan-like spread, an X-ray source that is disposed to face the X-ray source, and an X-ray from the X-ray source. An X-ray detector having a plurality of elements arranged in a direction orthogonal to the optical axis and in the direction of spread of the X-rays; and an X-ray detector arranged between the X-ray source and the X-ray detector and mounted with a subject. A turntable rotating around a rotation center axis orthogonal to both the optical axis and the arrangement direction of the elements, and a transmitted X-ray data string detected at predetermined angles by irradiating the subject with X-rays while driving the turntable Data conversion means for converting the rotation center axis into a parallel projection data sequence using the projection position of the rotation center axis on the X-ray detector, and using the data after the conversion, the object is projected on a plane orthogonal to the rotation center axis. Equipped with reconstruction operation means for generating sliced tomographic images In the computer tomography apparatus, the rotation center axis projection position used by the data conversion means is changed over a plurality of positions, and the rotation center axis projection position is estimated based on each data sequence converted using the rotation center axis projection position. A center axis projection position estimating unit, wherein the reconstruction calculating unit is characterized by forming a tomographic image from parallel projection data converted using the rotation center axis projection position estimated by the estimating unit. 1).
[0010]
Here, as a specific and preferred method of determining the optimum rotation center axis projection position by the rotation center axis projection position estimating means in the present invention, a plurality of rotations are changed as in the invention according to claim 2. For each sinogram of a plurality of parallel projection data sequences converted using each of the center axis projection positions, the symmetry of the 180 ° opposing data was evaluated, and the parallel projection data sequence with the highest symmetry was obtained. A method of estimating the axial projection position as optimal can be suitably adopted.
[0011]
According to the present invention, an X-ray transmission data sequence of a subject detected by an X-ray detector by irradiating the subject with fan-beam-shaped X-rays is converted into a parallel projection data sequence while performing sensitivity correction and distortion correction of the X-ray detector. At the time of conversion, the rotation center axis projection position is variously changed, and the true rotation center axis projection position or the rotation center axis projection position closest to the true rotation center axis projection position is calculated from the parallel projection data sequence converted using each rotation center axis projection position. By estimating, without using a dedicated phantom, the rotation center axis projection position is accurately obtained from transmitted X-ray data using an arbitrary subject.
[0012]
Specifically, a true rotation center axis projection position is used as a rotation center axis projection position used when converting a detection data sequence of X-rays in which a fan beam-shaped X-ray has transmitted through a subject into a parallel projection data sequence. In such a case, in the data string after the conversion, for example, data that is opposed by 180 ° in the sinogram should be symmetric. Therefore, the rotation center axis projection position used when converting a transmission X-ray data string obtained by rotating an arbitrary subject by 360 ° into a parallel projection data string is variously changed, and each rotation center axis projection position is changed. By evaluating the symmetry of the 180 ° -opposed data in each sinogram of the data sequence converted using the above, the rotation center axis projected position at which the data sequence with the highest symmetry was obtained becomes the true rotation center axis. It can be estimated that it is the projection position or the rotation center axis projection position closest thereto. Then, by reconstructing the tomographic image using the data sequence converted using the rotation center axis projected position estimated in this way, the obtained tomographic image always has high quality.
[0013]
BEST MODE FOR CARRYING OUT THE INVENTION
Hereinafter, embodiments of the present invention will be described with reference to the drawings.
FIG. 1 is a configuration diagram of an embodiment of the present invention, in which a schematic diagram showing a mechanical configuration of a main part and a block diagram showing an electrical functional configuration are shown together.
[0014]
The X-rays generated by the X-ray source 1 are converted into a fan-shaped beam B via the collimator 11 and emitted in the direction (x direction) toward the X-ray detector 2 which is disposed opposite. The X-ray detector 2 has a structure in which a plurality of elements are arranged in an arc shape in the X-ray spreading direction (y-direction).
[0015]
Between the X-ray source 1 and the X-ray detector 2, a direction (z-direction) orthogonal to both the x-direction, which is the optical axis direction of the X-ray, and the y-direction, which is the arrangement direction of the elements of the X-ray detector 2. A turntable 3 that rotates around a rotation center axis Ra along the line is disposed, and the subject W is mounted on the turntable 3 and is given rotation. The turntable 3 is mounted on an xyz table 31 that can move in the x, y, and z directions.
[0016]
The output of each element of the X-ray detector 2 is digitized by the data sampling circuit 4 every time the subject W is rotated by a small angle by the turntable 3 and is taken into the arithmetic unit 5. That is, the arithmetic unit 5 takes in 360-degree X-ray transmission data of the subject W for each element of the X-ray detector 2.
[0017]
The arithmetic unit 5 includes a preprocessing unit 51 that performs gamma correction, image distortion correction, sensitivity correction of each element, logarithmic conversion, and the like of the X-ray transmission data from the data sampling circuit 4, and an X-ray after the preprocessing. A data conversion unit 52 for converting the transmission data sequence into a parallel projection data sequence by performing a fan-para conversion, and a projection position of the rotation center axis Ra on the X-ray detector 2 necessary for the conversion, that is, a rotation center axis A rotation center axis projection position estimating unit 53 for estimating a projection position and a tomographic image reconstruction calculating unit 54 for generating a tomographic image using the data sequence converted by the data converting unit 52 are provided. Then, the tomographic image of the subject W generated by the tomographic image reconstruction calculation unit 54 is displayed on the display 6.
[0018]
The rotation center axis projection position estimating unit 53 estimates the rotation center axis projection position by the following method using 360 ° X-ray transmission data from each element of the X-ray detector 2 captured by the data sampling circuit 4. I do.
[0019]
That is, a plurality of rotation center axis projection positions are changed, and the captured data sequence is converted into a parallel projection data sequence using the respective rotation center axis projection positions. Then, with respect to the sinogram based on each data string after the conversion, the data which is opposed by 180 ° are compared with each other to evaluate the symmetry. Then, the rotation center axis projection position at which the data string with the highest symmetry is obtained is estimated as the true or closest rotation center axis projection position.
[0020]
FIG. 2 illustrates a sinogram. In this example, 360-degree data are arranged with the converted channel number (position in the y direction) on the horizontal axis and the rotation angle φ on the vertical axis. As a specific method of evaluating the symmetry of data that is 180 ° facing each other in such a sinogram, there are several methods, and each method will be described.
[0021]
First, in each sinogram, the data sequence of the rotation angle φ ° and the data sequence of (φ + 180) ° correspond to each other (as illustrated by an arrow in FIG. A variance value or a standard deviation value of a difference between data values (X-ray transmission amounts, in other words, brightness) of the channel data (located at symmetric positions) is obtained. Then, the rotation center axis projection position used for calculating the data sequence forming the sinogram with the smallest variance value or standard deviation value is estimated as the true projection position or the closest projection position.
[0022]
The second is that, in the sinogram, the difference between the maximum value and the minimum value of the sum of the data values of the corresponding channels in the data sequence of the rotation angle φ ° and the data sequence of (φ + 180) °, as described above. Then, the rotation center axis projection position used to calculate the data sequence forming the sinogram having the largest value is estimated as the true projection position or the closest projection position.
[0023]
Third, peak detection processing is performed on the data sequence of φ ° and the data sequence of (φ + 180) ° in each sinogram, and the peak position and the data value (X-ray transmission amount) at the peak position are obtained. Then, the peak position and its data value in the data sequence of φ ° are compared with the peak position and its data value in the data sequence of (φ + 180) ° to form a sinogram in which the mutual peak position and the data value are the best. The projection position of the rotation center axis used to calculate the data sequence is estimated as the true projection position or the closest projection position.
[0024]
In each of the above methods, φ may be an arbitrary single angle, or the value of φ may be changed from 0 to 180 °, and each evaluation result may be added or averaged to perform final processing. May be evaluated.
[0025]
The reconstruction calculation unit 54 reconstructs a tomographic image of the subject W using the parallel projection data sequence converted using the rotation center axis projection position estimated as described above.
Therefore, the tomographic image obtained by the above-described embodiment is always generated using the parallel projection data sequence converted using the accurate rotation center axis, and thus always has a high quality image.
[0026]
【The invention's effect】
As described above, according to the present invention, since the rotation center axis projection position is accurately estimated from the X-ray transmission image data of an arbitrary subject, it is not necessary to use a dedicated phantom as in the related art, There is no danger of an error in estimating the rotation center axis projection position due to manufacturing errors, and the rotation center axis projection position can always be stably and accurately estimated. It will be of high quality.
[Brief description of the drawings]
FIG. 1 is a configuration diagram of an embodiment of the present invention, and is a diagram illustrating both a schematic diagram illustrating a mechanical configuration of a main part and a block diagram illustrating an electrical functional configuration.
FIG. 2 is a diagram illustrating an example of a sinogram used for estimating a projection position of a rotation center axis performed by a rotation center axis projection position estimation unit 53 according to the embodiment of the present invention.
[Explanation of symbols]
REFERENCE SIGNS LIST 1 X-ray source 2 X-ray detector 3 Turntable 4 Data sampling circuit 5 Computing device 51 Preprocessing unit 52 Data conversion unit 53 Rotation center axis projection position estimating unit 54 Tomographic image reconstruction calculation unit 6 Display Ra Rotation center axis W subject

Claims (2)

ファン状に広がりを持つX線を出力するX線源と、そのX線源に対向配置され、当該X線源からのX線光軸に直交し、かつ、X線の広がり方向に複数の素子が配列されたX線検出器と、そのX線源とX線検出器の間に配置され、被写体を搭載して上記X線光軸および素子の配列方向の双方に直交する回転中心軸の回りを回転するターンテーブルと、上記ターンテーブルを駆動しつつ被写体にX線を照射して所定角度ごとに検出した透過X線データ列を、回転中心軸のX線検出器上への投影位置を用いて平行投影データ列に変換するデータ変換手段と、その変換後のデータを用いて、上記回転中心軸に直交する面で被写体をスライスした断層像生成する再構成演算手段を備えたコンピュータ断層撮影装置において、
上記データ変換手段で用いる回転中心軸投影位置を複数にわたって変化させ、その各回転中心軸投影位置を用いて変換した各データ列に基づいて回転中心軸投影位置を推定する回転中心軸投影位置推定手段を備え、上記再構成演算手段は、その推定手段により推定された回転中心軸投影位置を用いて変換した平行投影データから断層像を構成することを特徴とするコンピュータ断層撮影装置。
An X-ray source that outputs an X-ray having a fan-like spread, and a plurality of elements that are arranged to face the X-ray source, are orthogonal to the X-ray optical axis from the X-ray source, and extend in the X-ray spread direction. Are arranged between the X-ray detector and the X-ray source and the X-ray detector. The X-ray detector is mounted around a rotation center axis orthogonal to both the X-ray optical axis and the arrangement direction of the elements. Using a turntable that rotates, and a transmission X-ray data sequence detected at predetermined angles by irradiating the subject with X-rays while driving the turntable, using the projection position of the rotation center axis on the X-ray detector. Tomography data conversion means for converting the data into a parallel projection data sequence, and a reconstruction operation means for generating a tomographic image obtained by slicing a subject on a plane perpendicular to the rotation center axis using the converted data. At
Rotation center axis projection position estimating means for changing a plurality of rotation center axis projection positions used by the data conversion means and estimating a rotation center axis projection position based on each data string converted using each rotation center axis projection position. Wherein the reconstruction calculating means forms a tomographic image from the parallel projection data converted using the rotation center axis projection position estimated by the estimating means.
上記回転中心軸投影位置推定手段は、複数に変化させた回転中心軸投影位置のそれぞれを用いて変換した複数の平行投影データ列による各サイノグラムについて、180°対向するデータどうしの対称性を評価し、対称性が最も高くなる平行投影データ列を得た回転中心軸投影位置を最適と推定することを特徴とする請求項1に記載のコンピュータ断層像撮影装置。The rotation center axis projection position estimating means evaluates the symmetry of 180 ° facing data with respect to each sinogram obtained by converting a plurality of parallel projection data strings converted using each of the plurality of rotation center axis projection positions. 2. The computer tomographic imaging apparatus according to claim 1, wherein the rotation center axis projection position at which the parallel projection data string with the highest symmetry is obtained is estimated to be optimal.
JP2002219099A 2002-07-29 2002-07-29 Computed tomography equipment Expired - Lifetime JP3846577B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2002219099A JP3846577B2 (en) 2002-07-29 2002-07-29 Computed tomography equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2002219099A JP3846577B2 (en) 2002-07-29 2002-07-29 Computed tomography equipment

Publications (2)

Publication Number Publication Date
JP2004061256A true JP2004061256A (en) 2004-02-26
JP3846577B2 JP3846577B2 (en) 2006-11-15

Family

ID=31940085

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2002219099A Expired - Lifetime JP3846577B2 (en) 2002-07-29 2002-07-29 Computed tomography equipment

Country Status (1)

Country Link
JP (1) JP3846577B2 (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006003200A (en) * 2004-06-17 2006-01-05 Toshiba It & Control Systems Corp Computerized tomograph
JP2007064876A (en) * 2005-09-01 2007-03-15 Toshiba Corp Nuclear medicine inspection device
JP2007078469A (en) * 2005-09-13 2007-03-29 Toshiba It & Control Systems Corp Computerized-tomographic scanning device
US20070274456A1 (en) * 2006-05-23 2007-11-29 Bio-Imaging Research, Inc. Method and apparatus to facilitate determination of a parameter that corresponds to a scanning geometry characteristic
JP2009074997A (en) * 2007-09-21 2009-04-09 Toshiba It & Control Systems Corp Computerized tomographic apparatus
US8777485B2 (en) 2010-09-24 2014-07-15 Varian Medical Systems, Inc. Method and apparatus pertaining to computed tomography scanning using a calibration phantom

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009045092A (en) 2007-08-13 2009-03-05 Canon Inc Ct equipment and control method thereof

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006003200A (en) * 2004-06-17 2006-01-05 Toshiba It & Control Systems Corp Computerized tomograph
JP4504743B2 (en) * 2004-06-17 2010-07-14 東芝Itコントロールシステム株式会社 Computed tomography equipment
JP2007064876A (en) * 2005-09-01 2007-03-15 Toshiba Corp Nuclear medicine inspection device
JP2007078469A (en) * 2005-09-13 2007-03-29 Toshiba It & Control Systems Corp Computerized-tomographic scanning device
JP4733484B2 (en) * 2005-09-13 2011-07-27 東芝Itコントロールシステム株式会社 Computed tomography equipment
US20070274456A1 (en) * 2006-05-23 2007-11-29 Bio-Imaging Research, Inc. Method and apparatus to facilitate determination of a parameter that corresponds to a scanning geometry characteristic
US9052265B2 (en) * 2006-05-23 2015-06-09 Varian Medical Systems, Inc. Method and apparatus to facilitate determination of a parameter that corresponds to a scanning geometry characteristic
JP2009074997A (en) * 2007-09-21 2009-04-09 Toshiba It & Control Systems Corp Computerized tomographic apparatus
US8777485B2 (en) 2010-09-24 2014-07-15 Varian Medical Systems, Inc. Method and apparatus pertaining to computed tomography scanning using a calibration phantom

Also Published As

Publication number Publication date
JP3846577B2 (en) 2006-11-15

Similar Documents

Publication Publication Date Title
JP5142664B2 (en) X-ray computed tomography system
JP2004121446A (en) X-ray ct scanner
JP2009507544A (en) Direct measurement and scatter correction for CT
JP2011512187A (en) X-ray detector for phase contrast imaging
US6876718B2 (en) Scatter correction methods and apparatus
JP2011504395A (en) X-ray phase contrast imaging detection setup
JP2005087592A (en) X-rays measuring instrument
JP2006516433A (en) Voxel driven spiral reconstruction method for cone beam CT
JP2007252898A (en) Image display and x-ray ct scanner
JP3846576B2 (en) Computed tomography equipment
JP2004113785A (en) Image formation method and ct unit for implementing the same in computerized tomography
JP3846577B2 (en) Computed tomography equipment
US6411671B2 (en) Method for reducing line artifacts in a CT image and device for implementing the method
Yang et al. A review of geometric calibration for different 3-D X-ray imaging systems
WO2015054946A1 (en) Method and device for determining center of rotation of ct
JP2005233760A (en) Tomosynthesis device
JP2011147763A (en) Method for reconstructing image and x-ray computer tomography apparatus
JP2010204060A (en) X-ray inspection device, and inspection method of the same
JP2001170042A (en) Background monitoring method for ct data about presence and position of defective detector
JP2009150667A (en) Inspection device, alignment method of inspection device and method of manufacturing phantom for calibration
US7221730B2 (en) Multi-row detector x-ray CT apparatus and method for creating tomogram
JP2009276142A (en) Radiographic inspection system and imaging method of radiographic inspection
JP2022145494A (en) Image processing device, correction method, and program
JP2013116213A (en) X-ray ct apparatus
JP4299749B2 (en) CT image reconstruction method, CT apparatus, and program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20041125

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20060411

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060612

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: 20060802

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20060815

R150 Certificate of patent or registration of utility model

Ref document number: 3846577

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090901

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100901

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110901

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110901

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120901

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120901

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130901

Year of fee payment: 7

EXPY Cancellation because of completion of term