JP4644292B2 - X-ray CT apparatus and image display method thereof - Google Patents

X-ray CT apparatus and image display method thereof Download PDF

Info

Publication number
JP4644292B2
JP4644292B2 JP2009099614A JP2009099614A JP4644292B2 JP 4644292 B2 JP4644292 B2 JP 4644292B2 JP 2009099614 A JP2009099614 A JP 2009099614A JP 2009099614 A JP2009099614 A JP 2009099614A JP 4644292 B2 JP4644292 B2 JP 4644292B2
Authority
JP
Japan
Prior art keywords
subject
ray
tube current
dose
value
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 - Lifetime
Application number
JP2009099614A
Other languages
Japanese (ja)
Other versions
JP2009154021A (en
Inventor
孝行 長岡
宮崎  靖
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.)
Hitachi Healthcare Manufacturing Ltd
Original Assignee
Hitachi Medical 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 Hitachi Medical Corp filed Critical Hitachi Medical Corp
Priority to JP2009099614A priority Critical patent/JP4644292B2/en
Publication of JP2009154021A publication Critical patent/JP2009154021A/en
Application granted granted Critical
Publication of JP4644292B2 publication Critical patent/JP4644292B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Description

本発明は、撮影中のX線管電流(以下、管電流と略称する)を制御して被検体の被曝線量を抑制するX線CT装置に係り、特に被検体の被曝線量と画質とを考慮して、撮影中の管電流の変化曲線を再設定できるX線CT装置に関する。   The present invention relates to an X-ray CT apparatus that controls an X-ray tube current (hereinafter abbreviated as tube current) during imaging to suppress the exposure dose of the subject, and in particular, considers the exposure dose and image quality of the subject. The present invention relates to an X-ray CT apparatus that can reset a change curve of tube current during imaging.

従来のX線CT装置は、同一断層面では、同一のCTスキャン(以下、スキャンと略称する)条件(X線管電圧(以下、管電圧と略称する)や管電流等)で、撮影を行うようにしている。また、近年被検体を螺旋状にスキャンして撮影するヘリカルスキャンが広く用いられるようになってきているが、体軸方向のスキャン条件もスキャン中一定である。   Conventional X-ray CT apparatuses perform imaging under the same CT scan (hereinafter abbreviated as scan) conditions (X-ray tube voltage (hereinafter abbreviated as tube voltage), tube current, etc.) on the same tomographic plane. I am doing so. In recent years, a helical scan that scans and images a subject in a spiral manner has been widely used, but the scan condition in the body axis direction is also constant during the scan.

従って、例えば被検体の断面が、CTスキャナ(以下、スキャナと略称する)の回転軸に対して同心円ではなく楕円である場合、X線源の回転角度位置によって、被検体におけるX線の透過長が大きく変化するため、同一断層面内で、透過するX線量の過不足が発生する問題点を有していた。   Therefore, for example, when the cross section of the subject is an ellipse rather than a concentric circle with respect to the rotation axis of the CT scanner (hereinafter abbreviated as a scanner), the X-ray transmission length of the subject depends on the rotational angle position of the X-ray source. However, there was a problem that excessive and insufficient X-ray doses were transmitted within the same fault plane.

また、肺等の胸部の低密度の臓器と、肝臓等の腹部の高密度の臓器とでは、X線の吸収係数が大きく異なるため、胸部から上腹部へ連続的にスキャンを行う場合、肺に適するようなX線量を設定すると、肝臓では不足し、肝臓に適するようなX線量を設定すると、肺では過剰となる事態が生じていた。   In addition, the X-ray absorption coefficient differs greatly between low-density organs in the chest such as the lung and high-density organs in the abdomen such as the liver, so when scanning continuously from the chest to the upper abdomen, When the appropriate X-ray dose was set, there was a shortage in the liver, and when an X-ray dose appropriate for the liver was set, there was an excess in the lungs.

透過するX線量が不足する場合、X線検出器(以下、検出器と略称する)によって検出されるX線光子の量の減少によりS/N(SN比)が悪化し、結果として画像再構成によって得られる断層画像全体のS/Nが悪化する。逆に、透過するX線量が多すぎる場合には、被検体に対して、無効なX線被曝がなされていることになる。   When the transmitted X-ray dose is insufficient, the S / N (SN ratio) deteriorates due to a decrease in the amount of X-ray photons detected by an X-ray detector (hereinafter abbreviated as a detector), resulting in image reconstruction. The S / N of the entire tomographic image obtained by On the other hand, when the transmitted X-ray dose is too large, an invalid X-ray exposure is applied to the subject.

これらの問題点を解決する方法として、特許文献1で開示された管電圧を制御する方法や、特許文献2、特許文献3で開示された管電流を制御する方法が提案されている。   As a method for solving these problems, a method for controlling the tube voltage disclosed in Patent Document 1 and a method for controlling the tube current disclosed in Patent Document 2 and Patent Document 3 have been proposed.

特開昭53-110495号公報JP-A-53-110495 特開平9-108209号公報JP-A-9-108209 特開平10-309271号公報Japanese Patent Laid-Open No. 10-309271 特開2001-276040号公報Japanese Patent Laid-Open No. 2001-276040

しかしながら、特許文献1の管電圧を制御する方法は、スキャン中に管電圧を変化させるためにX線のスペクトルが変化し、CT値が決定できないという問題点がある。このため、現在では、管電流を制御する方法が主流となっている。   However, the method of controlling the tube voltage in Patent Document 1 has a problem that the CT value cannot be determined because the X-ray spectrum changes because the tube voltage is changed during scanning. For this reason, at present, the method of controlling the tube current has become the mainstream.

管電流を被検体に応じて最適に制御しようとする方法として、特許文献3のように、スキャナ回転の半周期前の透過X線量データを用いて管電流を制御する方法と、特許文献2のように、異なる2方向から撮影したスキャノグラムをもとにして、予め管電流を制御するパターンを被検体の位置に応じて決定しておく方法がある。   As a method of optimally controlling the tube current according to the subject, as in Patent Document 3, a method of controlling the tube current using transmitted X-ray dose data before a half cycle of scanner rotation, and Patent Document 2 Thus, there is a method in which a pattern for controlling the tube current is determined in advance according to the position of the subject based on scanograms taken from two different directions.

しかし、特許文献3のスキャナ回転の半周期前の透過X線量データを用いる方法は、特にヘリカルスキャンにて、スキャンピッチを大きくする場合に、透過X線量データのずれが大きくなる問題点がある。また、横隔膜前後のように、被検体のX線吸収特性が大きく変わる領域では、対応することができない。   However, the method using the transmitted X-ray dose data half a cycle before the rotation of the scanner in Patent Document 3 has a problem that a deviation of the transmitted X-ray dose data becomes large particularly when the scan pitch is increased in the helical scan. In addition, it is impossible to cope with a region where the X-ray absorption characteristics of the subject greatly change, such as before and after the diaphragm.

特許文献2の異なる2方向からのスキャノグラムを取得する方法は、スキャノグラム撮影を2回行うことにより、被検体への無用なX線被曝を増加させることになり、管電流制御による被曝線量の低減の目的に相反するものである。   The method of acquiring scanograms from two different directions in Patent Document 2 increases unnecessary X-ray exposure to the subject by performing scanogram imaging twice, and reduces the exposure dose by tube current control. It is contrary to the purpose.

また、本発明の発明者達も、管電流を被検体に応じて最適に制御する一方法を特許文献4にて提案している。特許文献4の発明は、被検体への無用なX線被曝を抑えて、被検体の低被曝化を実現したX線CT装置に関するもので、被検体に関するスキャナの回転角度とX線透過長との関係を示すモデルをメモリに格納しておき、被検体に対してのスキャン計測時に、このモデルから定まるスキャナの回転角度毎の設定管電流によるX線を、同じ被検体に照射してスキャン計測を行い、断層画像を再構成するものである。   The inventors of the present invention have also proposed in Patent Document 4 one method for optimally controlling the tube current according to the subject. The invention of Patent Document 4 relates to an X-ray CT apparatus that realizes low exposure of a subject by suppressing unnecessary X-ray exposure to the subject. A model showing the relationship between the two is stored in memory, and scan measurement is performed by irradiating the same subject with X-rays based on the set tube current for each rotation angle of the scanner determined from this model during scan measurement of the subject. And tomographic images are reconstructed.

特許文献4の発明では、被検体のX線透過長モデルの生成と、このX線透過長モデルのX線透過長に基づいて管電流を設定することに重点が置かれ、設定された管電流によって被検体に照射されるX線量や、被検体内の臓器へのX線被曝などについては考慮されていなかった。   In the invention of Patent Document 4, emphasis is placed on the generation of an X-ray transmission length model of a subject and setting the tube current based on the X-ray transmission length of the X-ray transmission length model. The X-ray dose applied to the subject by X-ray and the X-ray exposure to the organ in the subject were not considered.

上記の問題点を考慮し、本発明では、被検体のX線透過長モデルから自動的に設定した管電流の変化パターンによる被検体の内外における被曝線量を算定し、この被検体のX線被曝を考慮に入れて操作者が管電流の変化パターンを再設定することができるX線CT装置を提供することを目的とする。   In consideration of the above problems, in the present invention, the exposure dose inside and outside the subject is calculated from the tube current change pattern automatically set from the X-ray transmission length model of the subject, and the X-ray exposure of this subject is calculated. It is an object of the present invention to provide an X-ray CT apparatus in which the operator can reset the change pattern of the tube current in consideration of

上記目的を達成するため、本発明のX線CT装置は、X線源と対向配置され被検体の透過X線量を検出するX線検出器と、被検体のスキャノグラム像を表示する表示装置と、表示されたスキャノグラム像上に被検体の断層像の撮影位置を設定する撮影位置設定手段と、X線源を回転させながらそのX線源からX線を照射し、X線検出器により検出された被検体の透過X線から断層像を再構成する断層像再構成手段と、再構成された断層像を表示する表示装置とを備えたX線CT装置において、被検体の撮影位置に応じたX線をX線源に発生させるための管電流パターンを、管電流値と撮影位置とを対比させるためにスキャノグラム像と並置して、又は重ねて表示装置に表示させることを特徴とする。
To achieve the above object, the X-ray CT apparatus of the present invention, the X-ray detector for detecting transmitted X-ray amount of the subject arranged to face the X-ray source, a display device for displaying a scanogram image of the subject The X-ray source emits X-rays while rotating the X-ray source, and the X-ray detector detects the X-ray source while setting the X- ray source tomographic position on the displayed scanogram. and tomographic image reconstruction means for reconstructing a tomographic image from the transmitted X-ray amount of the subject was, in X-ray CT apparatus that includes a display device for displaying the reconstructed tomographic image, depending on the imaging position of the object In order to compare the tube current value and the imaging position, the tube current pattern for generating the X-rays generated in the X-ray source is displayed in parallel with the scanogram image or displayed on the display device.

また、本発明のX線CT装置では、前記管電流パターン計算手段は、スキャノグラム像に基づき被検体のX線透過長を計算する透過長計算手段を有し、X線透過長とスキャン条件に基づき管電流パターンを計算してもよい。In the X-ray CT apparatus of the present invention, the tube current pattern calculation means has transmission length calculation means for calculating the X-ray transmission length of the subject based on the scanogram image, and based on the X-ray transmission length and the scanning conditions. A tube current pattern may be calculated.

また、本発明のX線CT装置では、管電流パターンに応じた線量を計算する線量計算手段をさらに備え、計算された線量をスキャノグラム像と並置して、又は重ねて表示装置に表示させてもよい
Further, the X-ray CT apparatus of the present invention further includes a dose calculation means for calculating a dose according to the tube current pattern, and the calculated dose may be displayed on the display device in parallel with or superimposed on the scanogram image. Good .

本発明のX線CT装置の画像表示方法は、被検体のスキャノグラム像を表示する表示ステップと、表示されたスキャノグラム像上被検体の断層像の撮影位置設定された結果を受付ける撮影位置設定ステップと、X線源回転しながらX線を放射し、X線検出器により検出された被検体の透過X線から断層像を再構成する断層像再構成ステップと、再構成された断層像を表示する表示ステップとを備えたX線CT装置の画像表示方法において、被検体の撮影位置に応じたX線をX線源に発生させるための管電流パターンを、管電流値と撮影位置とを対比させるためにスキャノグラム像と並置して、又は重ねて表示するステップを備えることを特徴とする。
The image display method of the X-ray CT apparatus of the present invention, photographing position set to accept the results and display step of displaying a scanogram image of the subject, the imaging position of the tomographic image of the subject on the displayed scanogram image is set a step, the X-ray source emits a rotating article al X-ray, and the tomographic image reconstruction step of reconstructing a tomographic image from the transmitted X-ray amount of the subject detected by the X-ray detector, reconstituted In the image display method of the X-ray CT apparatus provided with a display step for displaying a tomographic image, a tube current pattern for generating an X-ray according to the imaging position of the subject in the X-ray source is expressed as a tube current value. In order to contrast with a photographing position, the method includes a step of displaying juxtaposed or superposed on a scanogram image.

また、本発明のX線CT装置の画像表示方法では、前記管電流パターン計算ステップは、スキャノグラム像に基づき被検体のX線透過長を計算する透過長計算ステップを有し、X線透過長とスキャン条件に基づき管電流パターンを計算してもよい。Further, in the image display method of the X-ray CT apparatus of the present invention, the tube current pattern calculation step includes a transmission length calculation step of calculating the X-ray transmission length of the subject based on the scanogram image, The tube current pattern may be calculated based on the scan condition.

また、本発明のX線CT装置の画像表示方法では、管電流パターンに応じた線量を計算する線量計算ステップをさらに備え、計算された線量をスキャノグラム像と並置して、又は重ねて表示装置に表示させてもよい。The X-ray CT apparatus image display method of the present invention further includes a dose calculation step of calculating a dose according to the tube current pattern, and the calculated dose is juxtaposed with the scanogram image or superimposed on the display device. It may be displayed.

本発明のX線CT装置の画像表示方法は、被検体の断層像を得るためのスキャノグラム像を得て、そのスキャノグラム像を表示装置に表示するステップと、表示されたスキャノグラム像と対応づけてX線源への供給電流の変化のパターンを表示するステップと、そのスキャノグラム像を参照しながら供給電流の変化パターンを編集するステップと、編集された供給電流の変化パターンに基づいて被検体の断層像を得て、その断層像を表示装置に表示するステップとを備える。The X-ray CT apparatus image display method of the present invention obtains a scanogram image for obtaining a tomographic image of a subject, displays the scanogram image on the display device, and associates the displayed scanogram image with the displayed scanogram image. A step of displaying a change pattern of the supply current to the radiation source, a step of editing the change pattern of the supply current while referring to the scanogram, and a tomographic image of the subject based on the edited change pattern of the supply current And displaying the tomographic image on a display device.

本発明のX線CT装置は、被検体の断層像を得るためのスキャノグラム像を取得するスキャノグラム像取得手段と、そのスキャノグラム像を表示する表示装置であって、表示されたスキャノグラム像と対応づけてX線源への供給電流の変化のパターンを表示する表示装置と、そのスキャノグラム像を参照しながら供給電流の変化パターンを編集する編集手段と、編集された供給電流の変化パターンに基づいて被検体の断層像を得て、その断層像を表示装置に表示する手段とを備える。An X-ray CT apparatus of the present invention is a scanogram image acquisition means for acquiring a scanogram image for obtaining a tomographic image of a subject, and a display device for displaying the scanogram image, in association with the displayed scanogram image. A display device that displays a change pattern of the supply current to the X-ray source, an editing means for editing the change pattern of the supply current while referring to the scanogram, and a subject based on the edited supply current change pattern Means for obtaining the tomographic image and displaying the tomographic image on a display device.

以上説明した如く、本発明のX線CT装置では、装置のスキャン条件を設定する操作手段や被検体のスキャノグラム画像データから被検体の3次元的X線透過長モデルを生成するスキャノグラム解析手段や、スキャン条件と被検体の3次元的X線透過長モデルとから被検体の撮影部位に応じた管電流の変化パターンを自動的に設定する管電流設定手段や、管電流の変化パターンに基づき被検体に照射される線量を計算し、表示する線量計算手段などを備えているので、スキャン条件として管電流の最大値、最小値を入力することにより、スキャン中の管電流の変化パターンを自動的に設定することができ、かつ被検体へのX線被曝も評価することができる。更に、被検体へのX線被曝が過剰になる恐れがある場合には、管電流の変化パターンを再設定することも可能である。   As described above, in the X-ray CT apparatus of the present invention, operation means for setting the scanning conditions of the apparatus, scanogram analysis means for generating a three-dimensional X-ray transmission length model of the subject from the scanogram image data of the subject, Tube current setting means for automatically setting the tube current change pattern according to the imaging region of the subject from the scan conditions and the three-dimensional X-ray transmission length model of the subject, and the subject based on the tube current change pattern Since the dose calculation means for calculating and displaying the dose irradiated to the tube is provided, the tube current change pattern during the scan is automatically entered by inputting the maximum and minimum tube current values as the scan conditions. It can be set and X-ray exposure to the subject can also be evaluated. Furthermore, when there is a possibility that the X-ray exposure to the subject becomes excessive, it is possible to reset the change pattern of the tube current.

また、本発明のX線CT装置では、スキャン中の管電流の変化パターンと予め生成した被検体の3次元CT値モデルデータに基づいて、被検体の体内の線量分布を計算し、表示する線量分布計算手段を備えているので、CT撮影による被検体の体内での線量分布を予め知ることができ、注目する臓器のX線被曝の程度を考慮して、スキャン実行の可否を判断することができる。   Further, in the X-ray CT apparatus of the present invention, the dose distribution in the body of the subject is calculated and displayed based on the change pattern of the tube current during the scan and the three-dimensional CT value model data of the subject generated in advance. Since it has a distribution calculation means, it is possible to know in advance the dose distribution in the body of the subject by CT imaging, and to determine whether or not to perform scanning in consideration of the degree of X-ray exposure of the organ of interest it can.

また、本発明のX線CT装置では、上記の被検体の3次元CT値モデルデータを人体ファントムなどをCT撮影して取得した標準人体CT値モデルデータと被検体のスキャノグラム画像データに基づいて生成する被検体CT値モデル生成手段を備えているので、以前にCT撮影を実施したことのない被検体についても、予備撮影としての1回のスキャノグラム画像データの取得のみによって、被検体CT値モデルデータの生成が可能である。   Further, in the X-ray CT apparatus of the present invention, the above three-dimensional CT value model data of the subject is generated based on the standard human CT value model data obtained by CT imaging of a human phantom or the like and the scanogram image data of the subject. CT image model data is generated only by acquiring one scanogram image data as a preliminary image, even for a sample that has not been subjected to CT imaging before. Can be generated.

また、本発明のX線CT装置では、表示手段の同一画面に被検体のスキャノグラム画像と管電流の変化パターンを並置又は重ねて表示しているので、操作者などは被検体の撮影部位を見ながら管電流の変化パターンの編集を行うことが可能となり、撮影部位に適した管電流の設定を容易に行うことができる。   In the X-ray CT apparatus of the present invention, the scanogram image of the subject and the change pattern of the tube current are displayed side by side or superimposed on the same screen of the display means. However, the change pattern of the tube current can be edited, and the tube current suitable for the imaging region can be easily set.

本発明に係るX線CT装置の全体構成を示すブロック図。1 is a block diagram showing the overall configuration of an X-ray CT apparatus according to the present invention. 本発明に係るX線CT装置の要部構成要素。The principal part component of the X-ray CT apparatus which concerns on this invention. 本発明に係るX線CT装置を使用した第1のスキャン操作例の一連の動作のフローチャート。6 is a flowchart of a series of operations in a first scan operation example using the X-ray CT apparatus according to the present invention. スキャノグラム画像とスライス位置及び投影データ例との対応を示す図。The figure which shows a response | compatibility with a scanogram image, a slice position, and the example of projection data. 3次元的X線透過長モデルの1スライス位置でのモデルを示す図。The figure which shows the model in the 1 slice position of a three-dimensional X-ray transmission length model. スライス位置ZにおけるX線透過長モデル。X-ray transmission length model at slice position Z. 管電流の変化パターンの表示例。A display example of a change pattern of tube current. 管電流の変化パターンの例。An example of a change pattern of tube current. 本発明に係るX線CT装置を使用した第2のスキャン操作例の一連の動作のフローチャート。9 is a flowchart of a series of operations in a second scan operation example using the X-ray CT apparatus according to the present invention. 被検体内の線量分布の計算結果の表示例。A display example of a calculation result of a dose distribution in a subject. 被検体CT値モデルの作成手順を説明するための図。The figure for demonstrating the preparation procedure of a subject CT value model. 照射X線量分布計算の手順を説明するための図。The figure for demonstrating the procedure of irradiation X dose distribution calculation. 被検体のスキャノグラム画像上に管電流の変化パターンを重畳して表示したもの。A superposed display of the tube current change pattern on the scanogram image of the subject. 管電流と被検体の厚さとの関係を示した図。The figure which showed the relationship between tube current and the thickness of a subject.

以下、添付図面を用いて、本発明の実施例について説明する。図1は、本発明に係るX線CT装置の全体構成を示すブロック図である。図1に示すように、このX線CT装置は、主としてX線源12と検出器13等を搭載し、被検体15に対してその周囲を連続回転可能なスキャナ11を内蔵するガントリ10と、装置全体を総括するホストコンピュータ20と、X線源12に高電圧を供給する高電圧発生装置22と、画像データの前処理や画像再構成処理、あるいは各種の解析処理を行う画像処理装置24と、画像を表示する表示装置25と、被検体15を載せるテーブル装置18と、操作者がスキャン条件等を入力する操作手段21等から成る。なお、スキャナ11と被検体15とは、相対的に回転することができればよいので、被検体15が静止してスキャナ11が回転してもよいし、スキャナ11が静止して被検体15の方が回転するとしてもよい。   Embodiments of the present invention will be described below with reference to the accompanying drawings. FIG. 1 is a block diagram showing the overall configuration of an X-ray CT apparatus according to the present invention. As shown in FIG. 1, the X-ray CT apparatus mainly includes an X-ray source 12 and a detector 13, and a gantry 10 including a scanner 11 that can continuously rotate around the subject 15. A host computer 20 for generalizing the entire apparatus, a high voltage generator 22 for supplying a high voltage to the X-ray source 12, an image processing apparatus 24 for performing preprocessing of image data, image reconstruction processing, and various analysis processes The display device 25 includes an image display device 25, a table device 18 on which the subject 15 is placed, an operation means 21 for an operator to input scan conditions and the like. Since the scanner 11 and the subject 15 need only be able to rotate relatively, the subject 15 may be stationary and the scanner 11 may be rotated, or the scanner 11 may be stationary and the subject 15 may be rotated. May rotate.

図2は、本発明に係るX線CT装置の要部構成要素を示す。図2を用いて、先ず、スキャナ11の詳細について説明する。図2において、スキャナ11には、X線源12と検出器13とが180度対向した位置関係で配置されている。このX線源12から放射されたX線ビーム14は、コリメータ19によってビームの幅及び厚さが制限されたファン状のX線ビーム14となって、被検体15に照射される。X線源12は、高電圧発生装置22を介してホストコンピュータ20によって制御される。スキャナ11全体は、スキャナ角度検出手段17によって回転角度を検出し、検出した回転角度に基づいてホストコンピュータ20が、スキャナ駆動手段16を制御し、スキャナ11を駆動する。検出器13は、被検体15を透過したX線14を検出し、検出データは被検体15によるX線の減衰量を示す投影データとして取り込まれる。投影データは、画像処理装置24において、ホストコンピュータ20の持つスキャナ角度等のデータと照合され、画像再構成等の処理をされた後に、表示装置25にて断層画像として表示される。   FIG. 2 shows the main components of the X-ray CT apparatus according to the present invention. First, details of the scanner 11 will be described with reference to FIG. In FIG. 2, the X-ray source 12 and the detector 13 are arranged in the scanner 11 in a positional relationship facing each other by 180 degrees. The X-ray beam 14 emitted from the X-ray source 12 is irradiated onto the subject 15 as a fan-shaped X-ray beam 14 whose beam width and thickness are limited by a collimator 19. The X-ray source 12 is controlled by the host computer 20 via the high voltage generator 22. The entire scanner 11 detects the rotation angle by the scanner angle detection means 17, and the host computer 20 controls the scanner drive means 16 based on the detected rotation angle to drive the scanner 11. The detector 13 detects the X-rays 14 that have passed through the subject 15, and the detection data is captured as projection data indicating the amount of X-ray attenuation by the subject 15. The projection data is collated with data such as the scanner angle of the host computer 20 in the image processing device 24, and after processing such as image reconstruction, it is displayed on the display device 25 as a tomographic image.

次に、本発明に係るスキャナ11以外の要部構成要素について説明する。図2において、装置全体を総括するホストコンピュータ20には、スキャナ11はスキャナ駆動手段16を介して、X線源12は高電圧発生装置22を介して、検出器13は画像処理装置24を介して、それぞれ間接的に接続されており、操作手段21と、管電流設定手段22と、画像処理装置24と、スキャノグラム解析手段26は直接接続されている。ホストコンピュータ20とスキャナ11、X線源12、検出器13との接続により、ホストコンピュータ20はX線源12による被検体15へのX線照射と、検出器13による投影データ(検出データ)の取り込みを制御する。画像処理装置24はホストコンピュータ20の指令に応じて、取り込まれた投影データに基づき、断層画像を順次再構成する。   Next, components other than the scanner 11 according to the present invention will be described. In FIG. 2, a host computer 20 that generalizes the entire apparatus includes a scanner 11 via a scanner driving means 16, an X-ray source 12 via a high voltage generator 22, and a detector 13 via an image processor 24. The operation means 21, tube current setting means 22, image processing device 24, and scanogram analysis means 26 are directly connected to each other. By connecting the host computer 20 to the scanner 11, the X-ray source 12, and the detector 13, the host computer 20 allows the X-ray irradiation to the subject 15 by the X-ray source 12 and the projection data (detection data) of the detector 13. Control uptake. The image processing device 24 sequentially reconstructs tomographic images based on the captured projection data in response to a command from the host computer 20.

本発明に係るX線CT装置では、被検体の断層画像を取得する本スキャンの前に、スキャン条件を設定するために種々の準備操作を行う。この準備操作としては、被検体の位置決めのためのスキャノグラム画像の撮影、管電流設定のためのスキャノグラム画像データの解析、スキャン条件としての管電流の変化パターンの決定などが、ホストコンピュータ20の介在のもとで行われる。   In the X-ray CT apparatus according to the present invention, various preparation operations are performed to set scan conditions before a main scan for acquiring a tomographic image of a subject. This preparatory operation includes scanning of a scanogram image for positioning the subject, analysis of scanogram image data for tube current setting, determination of a change pattern of tube current as a scan condition, etc. Done in the original.

これらの準備操作に関与する主な構成要素としては、図2において、ホストコンピュータ20と、操作手段21と、スキャノグラム解析手段26と、管電流設定手段22と、X線源12と、検出器13などである。この準備操作において、先ず、操作手段21は主として管電流の設定値(最大値、最小値)などのスキャン条件をシステムに入力する。X線源12と検出器13はスキャナ11を回転させずに、スキャノグラム画像の撮影を行い、画像データをホストコンピュータ20に保存する。スキャノグラム解析手段26はスキャノグラム画像データを解析し、被検体のX線透過長を体軸方向のスライス位置毎及びスキャナの回転角度毎に算出可能な、3次元形状データとしてモデル化し、このモデル(以下、被検体の3次元的X線透過長モデルという)のデータをホストコンピュータ20に保存する。管電流設定手段22は操作手段21から入力された管電流設定値と被検体の3次元的X線透過長モデルのデータを基にして、スキャン中に被検体の撮影部位のX線透過長の変化に応じて経時的に変化する一連の管電流値すなわち、管電流の変化パターンを自動的に決定する。このように決定された管電流の変化パターンは、ホストコンピュータ20に保存され、本スキャン時に被検体の撮影部位に応じて順次呼び出されて、X線源12の管電流を変化させる。   Main components involved in these preparatory operations are, in FIG. 2, a host computer 20, an operation means 21, a scanogram analysis means 26, a tube current setting means 22, an X-ray source 12, and a detector 13. Etc. In this preparation operation, first, the operation means 21 mainly inputs a scan condition such as a set value (maximum value, minimum value) of the tube current to the system. The X-ray source 12 and the detector 13 take a scanogram image without rotating the scanner 11 and store the image data in the host computer 20. The scanogram analysis means 26 analyzes the scanogram image data, and models the X-ray transmission length of the subject as three-dimensional shape data that can be calculated for each slice position in the body axis direction and for each rotation angle of the scanner. (Referred to as a three-dimensional X-ray transmission length model of the subject) is stored in the host computer 20. Based on the tube current setting value input from the operation means 21 and the data of the three-dimensional X-ray transmission length model of the subject, the tube current setting means 22 determines the X-ray transmission length of the imaging region of the subject during the scan. A series of tube current values that change with time according to the change, that is, a change pattern of the tube current is automatically determined. The tube current change pattern determined in this way is stored in the host computer 20 and is sequentially called according to the imaging region of the subject during the main scan to change the tube current of the X-ray source 12.

本発明では、更に上記で決定された管電流の変化パターンを基に、被検体15に照射されるX線量を本スキャン前に予め計算する。X線装置では、通常照射したX線量に対応する量としてmAs値が使用されているが、この計算でもmAs値を採用している。このmAs値は、管電流(mA)と照射時間(s)の積であり、管電圧が一定である場合(X線CT装置では管電圧が一定で使用される場合が多い)には、X線源12より照射されるX線量の総和に比例するために、X線量の基準として用いられる。ここで計算された被検体15に照射されるX線量については被検体15への予測被曝線量として操作者によって評価されることになる。   In the present invention, based on the change pattern of the tube current determined as described above, the X-ray dose irradiated to the subject 15 is calculated in advance before the main scan. In X-ray equipment, the mAs value is used as the amount corresponding to the normal X-ray dose, but this calculation also uses the mAs value. This mAs value is the product of the tube current (mA) and the irradiation time (s) .If the tube voltage is constant (the tube voltage is often used at a constant value in the X-ray CT system), X Since it is proportional to the total X-ray dose irradiated from the radiation source 12, it is used as a reference for the X-ray dose. The X-ray dose irradiated to the subject 15 calculated here is evaluated by the operator as the predicted exposure dose to the subject 15.

図3に、本発明に係るX線CT装置を使用した第1のスキャン操作の一連の動作のフローチャートを示す。このスキャン操作では、ステップ103、106、108の3次元データ生成、ステップ110の管電流パターン生成、ステップ111のmAs計算、ステップ114のmAs表示に特徴がある。以下、図2を参照しながら、図3の第1のスキャン操作の詳細について説明する。   FIG. 3 shows a flowchart of a series of operations of the first scanning operation using the X-ray CT apparatus according to the present invention. This scanning operation is characterized by three-dimensional data generation in steps 103, 106, and 108, tube current pattern generation in step 110, mAs calculation in step 111, and mAs display in step 114. Hereinafter, the details of the first scanning operation of FIG. 3 will be described with reference to FIG.

図3において、先ず、ステップ101のスキャノグラム撮影の工程では、被検体15のスキャノグラム画像を撮影する。被検体15のスキャノグラム画像を撮影する構成と断層画像を撮影する構成とは基本的には同じである。本ステップでは、スキャノグラム画像データは、スキャナ11を回転させずに、被検体15に正面方向からX線14を照射して、検出器13によって検出データを取り込むことによって得られる。このとき得られるスキャノグラム画像は正面方向のものである。このスキャノグラム画像データは、検出器13からホストコンピュータ20に送られる。このスキャノグラム画像データは、本スキャン時の被検体15の位置決めのために利用される他、本発明では特に管電流制御のための管電流の変化パターンの決定のために利用される。   In FIG. 3, first, in the scanogram imaging process of step 101, a scanogram image of the subject 15 is captured. The configuration for capturing a scanogram image of the subject 15 and the configuration for capturing a tomographic image are basically the same. In this step, scanogram image data is obtained by irradiating the subject 15 with X-rays 14 from the front direction without rotating the scanner 11 and capturing the detection data by the detector 13. The scanogram image obtained at this time is in the front direction. This scanogram image data is sent from the detector 13 to the host computer 20. The scanogram image data is used not only for positioning the subject 15 during the main scan, but also for determining the change pattern of the tube current for tube current control in the present invention.

次に、ステップ102のスキャノグラムデータ解析の工程及びステップ103の第1の3次元データ生成の工程では、スキャノグラム画像データがホストコンピュータ20に接続されたスキャノグラム解析手段26によって解析され、被検体15の3次元的X線透過長モデルが生成される。この3次元的X線透過長モデルは、被検体15をCT撮影する場合の、被検体15の位置とX線透過長との関係を示すモデルである。被検体15の3次元的X線透過長モデルの作成方法については特許文献4にも開示されている。   Next, in the scanogram data analysis step of step 102 and the first three-dimensional data generation step of step 103, the scanogram image data is analyzed by the scanogram analysis means 26 connected to the host computer 20, and the subject 15 A three-dimensional X-ray transmission length model is generated. This three-dimensional X-ray transmission length model is a model showing the relationship between the position of the subject 15 and the X-ray transmission length when the subject 15 is subjected to CT imaging. A method for creating a three-dimensional X-ray transmission length model of the subject 15 is also disclosed in Patent Document 4.

以下、被検体15の3次元的X線透過長モデルの作成方法の一例について説明する。図4はスキャノグラム画像とスライス位置及び投影データ例との対応を示す図、図5は3次元的X線透過長モデルの1スライス位置でのモデルを示す図である。図4(a)は、ステップ101で撮影された被検体のスキャノグラム画像29を示す。このスキャノグラム画像29は胸部から腹部の中間位置までの領域を撮影領域としている。このようなスキャノグラム画像の撮影領域の中からスライス位置が選ばれる。図示の場合n個のスライス位置が選ばれている。図中、P1、・・・、Pi、・・・、Pj、・・・、Pnがスライス位置である。 Hereinafter, an example of a method for creating a three-dimensional X-ray transmission length model of the subject 15 will be described. FIG. 4 is a diagram showing a correspondence between a scanogram image, a slice position, and an example of projection data, and FIG. 5 is a diagram showing a model at one slice position of a three-dimensional X-ray transmission length model. FIG. 4 (a) shows a scanogram image 29 of the subject imaged in step 101. In this scanogram image 29, an area from the chest to the middle position of the abdomen is an imaging area. A slice position is selected from the imaging region of such a scanogram image. In the case shown in the figure, n slice positions are selected. In the figure, P 1, ···, P i , ···, P j, ···, is P n is the slice position.

図4(b),(c),(d)は3次元的X線透過長モデルのモデル決定の説明図である。任意の2つのスライス位置Pi、PjでのCT断層像が図4(b)の如くなっているものと仮定すると、その縦方向(図示の上下方向)のX線減衰量の投影データは図4(c)の如くなるはずである。人体の体幹部の横断面は通常楕円形に近いものであるので、任意のスライス位置Pi、PjのCT断層像を図4(b)で仮定するのは、大きな誤差はないものと判断される。そこで、図4(c)の投影データについては、X線透過長データに換算し、その後横軸に沿って積分して面積を求める。このとき、X線減衰量の投影データのX線透過長への換算にあたっては、簡単のため人体が水と等価であるとみなしてデータの変換を行う。X線減衰量をc、X線透過長をb、水の線減弱係数をμwとしたとき、両者の関係はb=log c/μwで表される。また、横軸については、X線減衰量データの存在する領域全体の幅が人体の幅寸法に一致するように変換する。図4(d)は、図4(c)の投影データから変換したスライス位置Pi、Pjにおける被検体15のX線透過長データの分布図である。図4(d)においてスライス位置Pi、Pjにおける最大X線透過長はbi、bj、面積はSi、Sjとなる。図4(d)のX線透過長データについて、最大X線透過長bi、bj及び面積Si、Sjに注目すると、bi、Siはスライス位置Piでの断層画像のX線透過状況を反映し、bj、Sjはスライス位置Pjでの断層画像のX線透過状況を反映した値と見なすことができる。 4 (b), 4 (c), and 4 (d) are explanatory diagrams of model determination of a three-dimensional X-ray transmission length model. Assuming that the CT tomograms at any two slice positions P i and P j are as shown in Fig. 4 (b), the projection data of the X-ray attenuation in the vertical direction (vertical direction in the figure) is It should look like Figure 4 (c). Since the cross section of the trunk of the human body is usually close to an ellipse, assuming that CT tomograms at arbitrary slice positions P i and P j in Fig. 4 (b) are considered to have no significant error. Is done. Therefore, the projection data in FIG. 4 (c) is converted into X-ray transmission length data, and then integrated along the horizontal axis to obtain the area. At this time, in converting the X-ray attenuation amount into the X-ray transmission length of the projection data, the data is converted assuming that the human body is equivalent to water for simplicity. X-ray attenuation to c, and X-ray transmission length b, when the linear attenuation coefficient of water was mu w, the relationship between them is expressed by b = log c / μ w. Further, the horizontal axis is converted so that the width of the entire region where the X-ray attenuation data exists matches the width dimension of the human body. FIG. 4D is a distribution diagram of X-ray transmission length data of the subject 15 at the slice positions P i and P j converted from the projection data of FIG. 4C. In FIG. 4D, the maximum X-ray transmission lengths at the slice positions P i and P j are b i and b j , and the areas are S i and S j . In the X-ray transmission length data of FIG. 4D, when attention is paid to the maximum X-ray transmission lengths b i and b j and the areas S i and S j , b i and S i are X of the tomographic image at the slice position P i. Reflecting the line transmission state, b j and S j can be regarded as values reflecting the X-ray transmission state of the tomographic image at the slice position P j .

そこで、被検体15の3次元的X線透過長モデルとして、各スライス位置でのスライス断面を図5に示すような楕円形30でモデル化することにした。このモデル化では、スライス位置Pi、Pjにおける楕円形モデル30i、30jの面積をSi、Sj、短軸をbi、bjとしている。この結果、楕円形モデル30i、30jの長軸をai、ajとすると、楕円形モデル30i、30jの面積が[数1]で表されることから、 Therefore, as a three-dimensional X-ray transmission length model of the subject 15, the slice cross section at each slice position is modeled as an ellipse 30 as shown in FIG. In this modeling, the areas of the elliptic models 30 i and 30 j at the slice positions P i and P j are S i and S j , and the short axes are b i and b j . As a result, if the major axes of the elliptic models 30 i and 30 j are a i and a j , the area of the elliptic models 30 i and 30 j is expressed by [Equation 1]

Figure 0004644292
長軸ai、ajは[数2]によって求められる。
Figure 0004644292
The major axes a i and a j are obtained by [Equation 2].

Figure 0004644292

上記によって、各スライス位置での断層画像に対応するX線透過長モデルとしての楕円形モデル30が求められたので、これらの楕円形モデル30を体軸方向に配列することによって3次元的X線透過長モデル30を作成することができる。体軸方向のスライス位置のピッチが粗いときには、例えば隣り合う楕円形モデル間で最小自乗法により、その途中の1つ又は2つ以上の楕円形モデルを補間で求める。以上の如き手順により被検体15の3次元的X線透過長モデル30のデータとして、3次元座標(X、Y、Z)系での被検体15のX線透過長データT=T(X、Y、Z)が生成される。
Figure 0004644292

As described above, the elliptical model 30 as the X-ray transmission length model corresponding to the tomographic image at each slice position is obtained. By arranging these elliptical models 30 in the body axis direction, three-dimensional X-rays are obtained. A transmission length model 30 can be created. When the pitch of slice positions in the body axis direction is coarse, for example, one or two or more elliptical models in the middle are obtained by interpolation by the least square method between adjacent elliptical models. As described above, the data of the three-dimensional X-ray transmission length model 30 of the subject 15 as the data of the subject 15 in the three-dimensional coordinate (X, Y, Z) system X-ray transmission length data T = T (X, Y, Z) is generated.

次に、ステップ104からステップ110の工程では、上記の3次元的X線透過長モデル30を使用してX線源12に付与する管電流の変化パターンを設定することになるが、その前に3次元的X線透過長モデル30を使用した管電流の求め方について説明する。図4、図5で求めた3次元的X線透過長モデル30は、被検体の各スライス位置での断層画像のX線透過長を反映したものである。3次元的X線透過長モデル30のデータは一旦ホストコンピュータ20のレジスタを含むメモリに格納されているので、撮影範囲及びテーブルピッチなどのスキャン条件が決定するとその範囲のモデルのデータがメモリから取り出され、第2、第3の3次元データの生成及び管電流の変化パターンの決定に使用される。   Next, in the process from Step 104 to Step 110, the change pattern of the tube current to be applied to the X-ray source 12 is set using the above three-dimensional X-ray transmission length model 30, but before that, A method of obtaining the tube current using the three-dimensional X-ray transmission length model 30 will be described. The three-dimensional X-ray transmission length model 30 obtained in FIGS. 4 and 5 reflects the X-ray transmission length of the tomographic image at each slice position of the subject. Since the data of the three-dimensional X-ray transmission length model 30 is once stored in the memory including the register of the host computer 20, when scanning conditions such as the imaging range and table pitch are determined, the model data in that range is retrieved from the memory. The second and third three-dimensional data are generated and the change pattern of the tube current is determined.

管電流は、管電流設定手段23によって各スライス位置でスキャナ回転角度毎に3次元的X線透過長モデル30から得られるX線透過長に基づいて決定される。図6には、スライス位置(体軸方向の位置)ZにおけるX線透過長モデル30を示す。また、あるスキャナ回転角度における管電流は、通常そのスキャナ回転角度における3次元的X線透過長モデルのX線透過長のうちの最大値に対応付けて決定される。この最大値を示すX線透過長は図6の楕円形モデル30の中心0を通過するパスで得られるので、管電流を設定するにあたっては、スキャナ回転角度毎にこの楕円形モデル30の中心0を通過するパスのX線透過長のみを考慮すればよい。従って、図6においてスライス位置をZ、スキャナ回転角度をθ(θの始点は楕円形モデル30の短軸方向とする)としたとき、その位置における最大X線透過長Tは、Zとθの関数としてT=T(Z、θ)と表すことができる。   The tube current is determined by the tube current setting means 23 based on the X-ray transmission length obtained from the three-dimensional X-ray transmission length model 30 for each scanner rotation angle at each slice position. FIG. 6 shows the X-ray transmission length model 30 at the slice position (position in the body axis direction) Z. Further, the tube current at a certain scanner rotation angle is usually determined in association with the maximum value of the X-ray transmission lengths of the three-dimensional X-ray transmission length model at the scanner rotation angle. Since the X-ray transmission length indicating the maximum value is obtained by a path passing through the center 0 of the elliptic model 30 in FIG. 6, when setting the tube current, the center 0 of the elliptic model 30 is set for each scanner rotation angle. Only the X-ray transmission length of the path passing through the image needs to be considered. Accordingly, in FIG. 6, when the slice position is Z and the scanner rotation angle is θ (the start point of θ is the short axis direction of the elliptical model 30), the maximum X-ray transmission length T at that position is Z and θ. It can be expressed as T = T (Z, θ) as a function.

この最大X線透過長T(Z、θ)は、楕円形モデル30の中心位置0を通るパスの長さであるので、長軸をa、短軸をb、スキャナ回転角度をθとした場合、[数3]の如く表すことができる。   This maximum X-ray transmission length T (Z, θ) is the length of the path passing through the center position 0 of the elliptical model 30. Therefore, when the major axis is a, the minor axis is b, and the scanner rotation angle is θ , [Equation 3].

Figure 0004644292
ここで、a、bは[数1]、[数2]のai、aj、及びbi、bjと対応する。
Figure 0004644292
Here, a and b correspond to a i and a j and b i and b j in [Equation 1] and [Equation 2].

次に、管電流の設定方法の一例について説明する。先ず、被検体をスキャンする全範囲におけるパスの最大値(全スライス位置P1〜Pnの中でのパスの最大値)をTmax、最小値(同じくパスの最小値)をTminとする。これらの値は3次元的X線透過長モデル30を作るとき既知である。管電流を最大値Imax(mA)と最小値Imin(mA)の範囲で変化させる場合、本例では管電流の最大値、最小値とパスの最大値、最小値をそれぞれ対応させて、管電流とパスとの間に直線関係を持たせるものである。管電流IとパスTとの関係は[数4]の如く表される。 Next, an example of a tube current setting method will be described. First, let T max be the maximum value of the path (maximum value of the path among all slice positions P1 to Pn) in the entire range in which the subject is scanned, and T min be the minimum value (similarly, the minimum value of the path). These values are known when the three-dimensional X-ray transmission length model 30 is created. When changing the tube current in the range of the maximum value I max (mA) and the minimum value I min (mA), in this example, the maximum value of the tube current, the minimum value and the maximum value of the path, respectively, correspond to each other, A linear relationship is established between the tube current and the path. The relationship between the tube current I and the path T is expressed as [Equation 4].

Figure 0004644292
ここで、パスTはT(Z、θ)に対応するので、管電流Iは(Z、θ)の1次関数となり、スライス位置Z及びスキャナ回転角度θごとに求められる。
Figure 0004644292
Here, since the path T corresponds to T (Z, θ), the tube current I is a linear function of (Z, θ), and is obtained for each slice position Z and scanner rotation angle θ.

次に、図3のフローチャートに戻って説明する。ステップ104及びステップ105の工程では、操作者がスキャノグラム画像を参照して操作手段21からスキャン条件としてのテーブルピッチ及びスキャン開始位置を入力する。これらのデータにより、被検体のCT撮影範囲とスライス位置とスキャン回転角度が決定される。このときの座標系としては、上記の如く(Z、θ)座標系がよく、スキャン条件のデータも、(Z、θ)座標系のデータで入力するのがよい。   Next, returning to the flowchart of FIG. In steps 104 and 105, the operator inputs a table pitch and a scan start position as scan conditions from the operation means 21 with reference to the scanogram image. With these data, the CT imaging range, slice position, and scan rotation angle of the subject are determined. The coordinate system at this time is preferably the (Z, θ) coordinate system as described above, and the scan condition data is preferably input as data in the (Z, θ) coordinate system.

次に、ステップ106の工程では、第2の3次元的X線透過長モデルのデータを生成する。この工程で生成するデータは、各スライス位置Z、スキャン回転角度θ毎の最大X線透過長であり、第1の3次元的X線透過長モデルのデータから[数3]によって求めることができるので、ホストコンピュータ20のメモリから第1の3次元的X線透過長モデルのデータを呼び出して演算する。この演算結果は、T=T(Z、θ)で表される。   Next, in the process of step 106, data of a second three-dimensional X-ray transmission length model is generated. The data generated in this step is the maximum X-ray transmission length for each slice position Z and scan rotation angle θ, and can be obtained from [Expression 3] from the data of the first three-dimensional X-ray transmission length model. Therefore, the data of the first three-dimensional X-ray transmission length model is called from the memory of the host computer 20 and calculated. The calculation result is expressed as T = T (Z, θ).

次に、ステップ107の工程では、スキャン条件としてのスキャン時間を操作手段21より入力する。スキャン開始位置と、テーブルピッチと、スキャン時間が決定すると、スキャン中のX線源12の位置(Z、θ)はスキャン開始後の経過時間tの関数として表すことができるので、各スキャン位置での被検体15の第2の3次元的X線透過長モデル、すなわち最大X線透過長Tも時間tの関数T=T(t)として表すことができる。このため、ステップ108の第3の3次元的X線透過長モデル生成の工程では、最大X線透過長Tの関数を、T=T(Z、θ)からT=T(t)に変換する。   Next, in step 107, a scan time as a scan condition is input from the operation means 21. When the scan start position, table pitch, and scan time are determined, the position (Z, θ) of the X-ray source 12 during the scan can be expressed as a function of the elapsed time t after the start of the scan. The second three-dimensional X-ray transmission length model of the subject 15, that is, the maximum X-ray transmission length T can also be expressed as a function T = T (t) of time t. Therefore, in the third three-dimensional X-ray transmission length model generation step 108, the function of the maximum X-ray transmission length T is converted from T = T (Z, θ) to T = T (t). .

次に、ステップ109の管電流設定値入力の工程では、操作者が操作手段21より、管電流の設定値、例えばスキャン中の管電流の最大値Imaxと最小値Iminを入力する。ステップ110の管電流パターン生成の工程では、管電流設定手段23がホストコンピュータ20から3次元的X線透過長モデルのデータT(t)を呼び出し、上記の管電流設定値に基づき、被検体15の撮影部位に応じた管電流の変化パターンを自動的に決定する。このとき、X線透過長T(t)に対応付けてスキャン中の管電流の値を設定することになるが、X線透過長T(t)が最小の時には最小の管電流を設定し、X線透過長T(t)が最大の時には最大の管電流を設定するように管電流の変化パターンを決定する。また、X線透過長T(t)と管電流の値との関係として[数4]に示した1次関数の他、種々のものがある(但し、[数4]のTはT(t)に対応する)。 Next, in the step of inputting the tube current set value in step 109, the operator inputs the set value of the tube current, for example, the maximum value I max and the minimum value I min of the tube current during scanning, from the operation means 21. In the tube current pattern generation process of step 110, the tube current setting means 23 calls the data T (t) of the three-dimensional X-ray transmission length model from the host computer 20, and based on the tube current setting value, the subject 15 The tube current change pattern corresponding to the imaging region is automatically determined. At this time, the value of the tube current during scanning is set in association with the X-ray transmission length T (t), but when the X-ray transmission length T (t) is minimum, the minimum tube current is set, When the X-ray transmission length T (t) is the maximum, the change pattern of the tube current is determined so as to set the maximum tube current. In addition to the linear function shown in [Equation 4], there are various relationships between the X-ray transmission length T (t) and the tube current value (however, T in [Equation 4] is T (t )).

上記の如くして被検体15の3次元的X線透過長モデルに合わせて、管電流が時間tの関数として決定される。従って、管電流の変化パターンは、I=I(t)と表すことができる。このように決定された管電流の変化パターンI=I(t)はホストコンピュータ20に保存され、本スキャン時に被検体15の撮影部位に応じて順次呼び出されて、高電圧発生装置22を介してスキャン中の管電流を制御する。   As described above, the tube current is determined as a function of time t in accordance with the three-dimensional X-ray transmission length model of the subject 15. Therefore, the change pattern of the tube current can be expressed as I = I (t). The tube current change pattern I = I (t) determined in this way is stored in the host computer 20 and is sequentially called in accordance with the imaging region of the subject 15 during the main scan, via the high voltage generator 22. Control the tube current during the scan.

図7に、管電流の変化パターンの表示例を示す(本例も含めて以下の管電流の変化パターンの表示例では、簡単のためスキャナ回転角度θの変化に伴う大略周期的な管電流の変化については省略し、スライス位置による変化のみ示している)。これは、表示装置25の画面上にスキャノグラム画像29と対比して表示したものである。管電流の変化パターン31では、縦軸に管電流値(mA)、横軸にスキャン開始後の経過時間tをとっている。表示例の場合、管電流はスキャンの初期(腹部)には中程度の値で、中期(腹部と胸部との間)には最小値となり、終期(胸部)には最大値となっている。本表示例の如く、管電流の変化パターン31とスキャノグラム画像29を同一画面上に並置することにより、管電流値と撮影部位との対比を一目で行うことができるので、管電流値の妥当性の判断に有効である。   FIG. 7 shows a display example of the tube current change pattern (in the following display examples of the tube current change pattern including this example, for the sake of simplicity, an approximately periodic tube current change due to the change in the scanner rotation angle θ is shown. The change is omitted and only the change due to the slice position is shown). This is displayed in contrast to the scanogram image 29 on the screen of the display device 25. In the tube current change pattern 31, the vertical axis represents the tube current value (mA), and the horizontal axis represents the elapsed time t after the start of scanning. In the case of the display example, the tube current has a medium value at the beginning of the scan (abdomen), a minimum value at the middle period (between the abdomen and the chest), and a maximum value at the end (chest). As shown in this display example, the tube current change pattern 31 and the scanogram image 29 are juxtaposed on the same screen so that the tube current value and the imaging region can be compared at a glance. It is effective for judgment.

次に、ステップ111のmAs計算の工程では、ステップ110で決定された管電流の変化パターンに基づき、スキャン中に被検体15に照射されるX線量を計算する。ここで、被検体15に照射されるX線量の基準としては上述の如く管電流(mA)と照射時間(s)との積であるmAsが用いられる。従って、本ステップでは、管電流の変化パターンI=I(t)を時間で積分して、被検体15に照射されるX線量mAsを求める。この積分は[数5]によって行われる。   Next, in the step of calculating mAs in step 111, the X-ray dose irradiated to the subject 15 during the scan is calculated based on the change pattern of the tube current determined in step 110. Here, as a reference for the X-ray dose irradiated to the subject 15, mAs, which is the product of the tube current (mA) and the irradiation time (s), is used as described above. Therefore, in this step, the tube current change pattern I = I (t) is integrated over time to obtain the X-ray dose mAs irradiated to the subject 15. This integration is performed by [Equation 5].

Figure 0004644292
ここで求めたmAs値は、あくまでも被検体15に照射されるX線量に相当する量であるので、実験等によってX線量とmAs値との正確な対応をとり、両者の間の換算ができるようにしておく必要がある。
Figure 0004644292
Since the mAs value obtained here is an amount corresponding to the X-ray dose irradiated to the subject 15 only, the X-ray dose and the mAs value can be accurately matched by experiments, etc. It is necessary to keep it.

図7に示した管電流の変化パターン31を例に上げて、被検体15への照射X線量mAsを計算することを考える。この場合、管電流の変化パターン31は、I=I(t)で、時間の関数であるから、これを積分してmAsを計算することは、管電流の変化パターン31の図の面積を計算することになり、管電流の変化パターン31の図の面積SがmAsに相当する。   Taking the tube current change pattern 31 shown in FIG. 7 as an example, let us consider calculating the irradiation X-ray dose mAs to the subject 15. In this case, the tube current change pattern 31 is I = I (t) and is a function of time, so calculating the mAs by integrating this calculates the area of the tube current change pattern 31 figure. Therefore, the area S in the diagram of the tube current change pattern 31 corresponds to mAs.

次に、ステップ112のmAs計算値表示の工程では、ステップ111にて求めたmAsの計算値を表示装置25の画面に表示する。ステップ111では、管電流の変化パターンが生成された段階で、被検体15のスキャン範囲の全領域にわたって被検体に照射されるX線量に対応するmAsが計算されるので、このステップ112でこのmAs値がスキャンを開始してもよいか否かを操作者が判断するための資料として操作者に提示される。   Next, in the step of displaying the calculated mAs value in step 112, the calculated value of mAs obtained in step 111 is displayed on the screen of the display device 25. In step 111, when a change pattern of the tube current is generated, mAs corresponding to the X-ray dose irradiated to the subject over the entire region of the scan range of the subject 15 is calculated. The value is presented to the operator as data for the operator to determine whether or not to start scanning.

ステップ113のmAs判断の工程では、操作者が全体としてのmAsの妥当性について判断する。すなわち、操作者はこれから行われるCT撮影による利益と被検体15へのX線被曝による不利益とを比較考量して、全体としてのmAsが大きすぎないか否かを判断し、全体としてのmAsが大きすぎると判断した場合には管電流の設定値を下げることになる。この場合には、ステップ109に戻り、管電流設定値を再度入力し、管電流の変化パターンを再設定する。   In the step of determining mAs in step 113, the operator determines the validity of the mAs as a whole. That is, the operator weighs the benefits of CT imaging to be performed in the future and the disadvantages of X-ray exposure to the subject 15 to determine whether or not the overall mAs is too large. If it is determined that is too large, the set value of the tube current is lowered. In this case, the process returns to step 109, the tube current set value is input again, and the change pattern of the tube current is reset.

図8には、管電流の変化パターンの例を示す。図8(a)に示した管電流の変化パターン31aは、通常の管電流一定(I0)の場合で、管電流の最小値(Imin)と最大値(Imax)が同じI0である場合の例である。図8(a)の管電流の変化パターン31aのグラフでは、mAsは面積Saである。次に、図8(b)示した管電流の変化パターン31bでは、管電流は初期には図8(a)と同じI0で、中期にI0より低い最小値Iminとなり、終期にI0より高い最大値Imaxとなっている。図8(b)の管電流の変化パターン31aのグラフでのmAsは面積Sbであるが、この面積Sbは図8(a)の面積Saより小さくなっており、被検体15の被曝線量は低減している。 FIG. 8 shows an example of a change pattern of the tube current. The tube current change pattern 31a shown in Fig. 8 (a) is the case where the normal tube current is constant (I0), and the minimum value (I min ) and the maximum value (I max ) of the tube current are the same I0. It is an example. In the graph of the tube current change pattern 31a in FIG. 8A, mAs is the area Sa. Next, in FIG. 8 (b) shows the tube current change pattern 31b, the tube current in the same I0 and 8 early (a), the mid-term to the minimum value I min next lower I0, higher than I0 at the end The maximum value is I max . The mAs in the graph of the tube current change pattern 31a in FIG. 8 (b) is the area Sb, but this area Sb is smaller than the area Sa in FIG. 8 (a), and the exposure dose of the subject 15 is reduced. is doing.

図8(c)に示した管電流の変化パターン31cでは、図8(b)のグラフと比べて、終期の最大値Imaxを小さくし、I0とほぼ同じとしたものである。この場合、スキャン全体として管電流が低目に抑えられている。図8(c)の管電流の変化パターン31cのグラフでのmAsは面積Scであるが、この面積Scは図8(b)の面積Sbより更に小さくなっており、被検体15の被曝線量は更に低減している。 In the tube current change pattern 31c shown in FIG. 8 (c), the maximum value I max at the end is made smaller than that in the graph of FIG. 8 (b) and is almost the same as I0. In this case, the tube current is suppressed to a low level for the entire scan. The mAs in the graph of the tube current change pattern 31c in FIG.8 (c) is the area Sc, but this area Sc is smaller than the area Sb in FIG.8 (b), and the exposure dose of the subject 15 is Further reduction.

図8(d)に示した管電流の変化パターン31dでは、図8(b)に比べて中期の最小値Iminを更に小さくし、IminがI0に対し大幅に低下するようにしたものである。図8(d)の管電流の変化パターン31dでのmAsは面積Sdであるが、図8(c)の面積Scより更に小さくなっている。 In the tube current change pattern 31d shown in FIG. 8 (d), the minimum value I min in the middle period is further reduced as compared with FIG. 8 (b), and I min is greatly reduced with respect to I0. is there. The mAs in the tube current change pattern 31d in FIG. 8 (d) is the area Sd, but is smaller than the area Sc in FIG. 8 (c).

図8(c)の場合の如く管電流設定値の最大値を下げた場合には、被検体15の厚さが厚い部分、すなわちX線透過長が大きい部分での管電流を減らすことになる。従って、図8(c)の管電流の変化パターン31cは、特に肺のように低密度の部位での画質を重視し、腹部のように高密度の領域での被曝線量を低減したい場合に適する。   When the maximum value of the tube current set value is lowered as in FIG. 8 (c), the tube current is reduced at the portion where the subject 15 is thick, that is, the portion where the X-ray transmission length is large. . Therefore, the tube current change pattern 31c in FIG. 8 (c) is particularly suitable when the image quality in a low-density region such as the lung is emphasized and the exposure dose in a high-density region such as the abdomen is desired to be reduced. .

図8(d)の場合の如く管電流設定値の最小値を下げた場合には、被検体15の厚さが薄い部分、すなわちX線透過長が小さい部分での管電流を減らすことになる。従って、図8(d)の管電流の変化パターン31dは、特に骨周辺や実質部のように高密度の部位での画質を重視し、低密度の領域での被曝線量を低減したい場合に適する。   When the minimum value of the tube current set value is lowered as in the case of FIG. 8 (d), the tube current is reduced in the portion where the subject 15 is thin, that is, the portion where the X-ray transmission length is small. . Therefore, the tube current change pattern 31d in FIG. 8 (d) is particularly suitable when the image quality in a high-density region such as the bone periphery or the substantial part is emphasized and the exposure dose in a low-density region is desired to be reduced. .

図8(c),図8(d)では、管電流設定値の最大値Imax又は最小値Iminのいずれか一方のみを小さくしているが、両方を小さくして、被検体15の全領域について平均的に被曝線量を低減することも可能である。上記の如く新たに設定した管電流の変化パターンについては、再度mAsを計算し、問題が無い場合には、そのまま採用する。 In FIG. 8 (c) and FIG. 8 (d), only one of the maximum value I max and the minimum value I min of the tube current set value is reduced, but both are reduced so that the whole of the subject 15 is reduced. It is also possible to reduce the exposure dose on average for the region. For the tube current change pattern newly set as described above, mAs is calculated again, and if there is no problem, it is adopted as it is.

次に、ステップ114のスキャンの工程では、操作者は上記で決定した管電流の変化パターンを含めたスキャン条件にてスキャンを実行する。   Next, in the scanning process of step 114, the operator executes scanning under the scanning conditions including the tube current change pattern determined above.

次に、ステップ115のmAs積算値表示の工程では、上記の[数5]に従って、スキャン中のmAsの積算値を逐次計算し、表示装置25の画面にリアルタイムで表示し、操作者に提示する。mAs積算値表示方法としては、全体としてのmAs値に対する比率としての相対値を表示する方法と、mAs積算値の絶対値を表示する方法のいずれかを選択することができる。もちろん、その両方を同時に表示することも可能である。   Next, in the step of displaying the mAs integrated value in step 115, the integrated value of mAs during scanning is sequentially calculated according to the above [Formula 5], displayed in real time on the screen of the display device 25, and presented to the operator. . As the mAs integrated value display method, either a method of displaying a relative value as a ratio to the mAs value as a whole or a method of displaying an absolute value of the mAs integrated value can be selected. Of course, both can be displayed simultaneously.

次に、本発明に係るX線CT装置の第2のスキャン操作例について説明する。図9は、第2のスキャン操作例の一連の動作のフローチャートを示したものである。このスキャン操作例では、図3の第1のスキャン操作例のフローチャートに対し、ステップ111のmAs計算及びステップ112のmAs表示の工程の後に、被検体15内の線量分布の計算及び表示の工程を追加し、操作者が被検体15内の線量分布をも見てスキャン実行の判断を下せるようにしたものである。このため、本操作例の説明では、ステップ201〜205の被検体15内の線量分布の計算及び表示の工程を重点に説明する。   Next, a second scanning operation example of the X-ray CT apparatus according to the present invention will be described. FIG. 9 shows a flowchart of a series of operations of the second scan operation example. In this scan operation example, with respect to the flowchart of the first scan operation example in FIG. 3, the steps of calculating and displaying the dose distribution in the subject 15 are performed after the mAs calculation in step 111 and the mAs display step in step 112. In addition, the operator can also determine the scan execution by looking at the dose distribution in the subject 15. For this reason, in the description of this operation example, the process of calculating and displaying the dose distribution in the subject 15 in steps 201 to 205 will be mainly described.

以下、図9の第2のスキャン操作例の中のステップ201からステップ205について説明するが、ステップの内容説明に入る前に、被検体15内の線量分布の計算結果の表示例を図10に示す。図10(a)は被検体の横断面35での線量分布の表示例、図10(b)は被検体の体軸方向の側面36での線量分布の表示例である。両図とも被検体15内の分布線量の等しい等線量線38a〜38c、39a〜39cが示されており、体表に近いほど高線量になっている。本操作例では、被検体15の線量分布が操作者に提示されるため、操作者は被検体のX線被曝についてより詳細な評価を行うことができるということが大きな特徴となる。   Hereinafter, step 201 to step 205 in the second scan operation example of FIG. 9 will be described. Before entering the description of the contents of the step, a display example of the calculation result of the dose distribution in the subject 15 is shown in FIG. Show. FIG. 10 (a) is a display example of the dose distribution on the cross section 35 of the subject, and FIG. 10 (b) is a display example of the dose distribution on the side surface 36 in the body axis direction of the subject. In both figures, isodose lines 38a to 38c and 39a to 39c having the same distributed dose in the subject 15 are shown, and the closer to the body surface, the higher the dose. In the present operation example, since the dose distribution of the subject 15 is presented to the operator, the operator can perform a more detailed evaluation on the X-ray exposure of the subject.

本操作例では、ステップ203で被検体内の線量分布の計算を行う前に、準備として、ステップ201で被検体のCT値モデルの生成、ステップ202で被検体のμモデルの生成を行い、その後で被検体のμモデルデータと被検体に照射されたX線量データを基にして被検体内の線量分布を計算する。   In this operation example, before calculating the dose distribution in the subject in step 203, as a preparation, the CT value model of the subject is generated in step 201, the μ model of the subject is generated in step 202, and then The dose distribution in the subject is calculated based on the μ model data of the subject and the X-ray dose data irradiated on the subject.

先ず、ステップ201の被検体CT値モデル生成の工程では、予め標準的な人体のCT値分布モデル(以下、標準人体CT値モデルという)データを取得して、ホストコンピュータ20の記憶手段に保存しておき、この標準人体CT値モデルデータを、ステップ101で取得した被検体15のスキャノグラム画像のデータに基づいて補正を行うことによって、被検体15のCT値モデル(以下、被検体CT値モデルという)を作成する。上記の標準人体CT値モデルデータとしては、例えば標準的な人体ファントムなどをCT撮影した断層画像から得られる3次元的CT値分布データが用いられる。断層画像はCT値の分布を表し、実効エネルギー(通常60keV)のX線に対する線減弱係数の分布を表しているので、この断層画像を3次元に再構成をした3次元的CT値分布データは線減弱係数の3次元的な空間分布のデータであり、被検体に照射されたX線の減弱量の計算に利用することができる。   First, in the subject CT value model generation step of step 201, standard human CT value distribution model (hereinafter referred to as standard human CT value model) data is acquired in advance and stored in the storage means of the host computer 20. The standard human body CT value model data is corrected based on the scanogram image data of the subject 15 acquired in step 101, whereby the CT value model of the subject 15 (hereinafter referred to as the subject CT value model). ). As the standard human CT value model data, for example, three-dimensional CT value distribution data obtained from a tomographic image obtained by CT imaging of a standard human phantom or the like is used. The tomographic image represents the distribution of CT values, and the distribution of the linear attenuation coefficient for effective energy (usually 60 keV) X-rays, so the three-dimensional CT value distribution data obtained by reconstructing this tomographic image in three dimensions is This is data of the three-dimensional spatial distribution of the line attenuation coefficient, and can be used to calculate the amount of attenuation of X-rays irradiated to the subject.

次に、被検体CT値モデルの生成方法の一例について、図11を用いて説明する。この例では、ステップ101で被検体15を撮影して取得して実測したスキャノグラム画像データと、上記の標準人体CT値モデルデータとから、被検体15の3次元的CT値分布を表す被検体CT値モデルデータを生成する。この被検体CT値モデルの生成にあたっては、被検体15のスキャノグラム画像データと標準的人体のスキャノグラム画像データが媒体として利用される。   Next, an example of a method for generating a subject CT value model will be described with reference to FIG. In this example, the subject CT representing the three-dimensional CT value distribution of the subject 15 from the scanogram image data obtained by photographing the subject 15 in step 101 and actually measured, and the standard human CT value model data described above. Generate value model data. In generating the subject CT value model, the scanogram image data of the subject 15 and the scanogram image data of the standard human body are used as a medium.

図11は、ステップ201の被検体CT値モデルの作成手順を説明するための図である。図11において、図11(a)は標準人体CT値モデルデータ41の例を、図11(b)は図11(a)の標準人体CT値モデルデータ41から計算で求めた標準人体のスキャノグラム画像データ42の例を、図11(c)は被検体の実測のスキャノグラム画像データ43の例を、図11(d)は計算で求めた被検体CT値モデルデータ44の例を示す。図11(a)の標準人体CT値モデルデータ41は、人体ファントムなどの標準的な人体の体幹部のCT値分布モデルで、肩から腹部までのスライス位置ごとの断面のCT値分布モデルを示している。   FIG. 11 is a diagram for explaining a procedure for creating a subject CT value model in step 201. In FIG. 11, FIG. 11 (a) is an example of the standard human CT value model data 41, and FIG. 11 (b) is a scanogram image of the standard human body calculated from the standard human CT value model data 41 of FIG. 11 (a). FIG. 11 (c) shows an example of the data 42, FIG. 11 (d) shows an example of the actually measured scanogram image data 43, and FIG. 11 (d) shows an example of the subject CT value model data 44 obtained by calculation. The standard human CT value model data 41 in Fig. 11 (a) is a CT value distribution model of the trunk of a standard human body such as a human phantom, and shows the CT value distribution model of the cross section at each slice position from the shoulder to the abdomen. ing.

スキャノグラム画像は、その3次元的CT値分布モデルから計算によって生成することができるので、図11(b)の標準人体スキャノグラム画像データ42は図11(a)の標準人体CT値モデルデータ41について正面方向から投影したデータを求めることによって得られる。図11(c)の被検体の実測スキャノグラム画像データ43は、被検体15の体幹部について標準人体スキャノグラム画像データ42と同じ領域を正面方向から撮影したスキャノグラム画像データである。この画像については、以下被検体スキャノグラム画像と呼ぶことにする。   Since the scanogram image can be generated by calculation from the three-dimensional CT value distribution model, the standard human body scanogram data 42 in FIG. 11 (b) is a front view of the standard human CT value model data 41 in FIG. 11 (a). It is obtained by obtaining data projected from the direction. The measured scanogram image data 43 of the subject in FIG. 11 (c) is scanogram image data obtained by photographing the same area as the standard human scanogram image data 42 from the front direction of the trunk of the subject 15. This image is hereinafter referred to as a subject scanogram image.

図11においては、体幹部の標準人体スキャノグラム画像データ42と被検体スキャノグラム画像データ43とを並置して対比できるように示してあるが、両者は普通寸法及びCT値とも異なるものである。このため、標準人体スキャノグラム画像データ42と被検体スキャノグラム画像データ43とを対比しながら、両者の差異に基づいて、一致する部分はそのままとし、異なる部分については変形させて、被検体15に合うように標準人体CT値モデルデータ41を補正して、被検体CT値モデルデータ44を生成する。   In FIG. 11, the standard human body scanogram image data 42 of the trunk and the subject scanogram image data 43 are shown so as to be juxtaposed, but they are different from each other in normal dimensions and CT values. For this reason, while comparing the standard human body scanogram image data 42 and the subject scanogram image data 43, based on the difference between the two, the matching portions are left as they are, and the different portions are deformed to fit the subject 15. Then, the standard human body CT value model data 41 is corrected to generate the subject CT value model data 44.

図11の体幹部の例では、先ず体軸方向に関して、標準人体スキャノグラム画像データ42と被検体スキャノグラム画像データ43の肩から横隔膜までの長さAと、横隔膜から腸管までの長さBとに分けて、それぞれの差異を基に、標準人体CT値モデルデータ41を補間、伸長したり、あるいは間引き、短縮したりすることで、標準人体CT値モデルデータ41の体軸方向のCT値分布を被検体15の実状に近似させる。左右方向に関しても、同様に、体軸を基準にして左と右に分けて、それぞれの差異を基に、左右の広がりを補正し、被検体15の実状に近似させる。前後方向に関しては、被検体スキャノグラム画像データ43から推定される前後方向のX線透過長を基に、標準人体CT値モデルデータ41の前後方向のデータを線形に補間する。このようにして、2つのスキャノグラム画像データ42、43を基に、標準人体CT値モデルデータ41を実際の被検体15に合わせ込むことにより、被検体CT値モデルデータ44を生成する。   In the example of the trunk in FIG. 11, first, the body axis direction is divided into the length A from the shoulder to the diaphragm and the length B from the diaphragm to the intestinal tract of the standard human scanogram image data 42 and the subject scanogram image data 43. Based on these differences, the standard human CT value model data 41 is interpolated, expanded, thinned out, or shortened, so that the CT value distribution in the body axis direction of the standard human CT value model data 41 is covered. The actual condition of the specimen 15 is approximated. Similarly, in the left-right direction, the left and right sides are divided with respect to the body axis, and the left and right spreads are corrected based on the respective differences to approximate the actual state of the subject 15. With respect to the front-rear direction, based on the X-ray transmission length in the front-rear direction estimated from the subject scanogram image data 43, the data in the front-rear direction of the standard human CT value model data 41 is linearly interpolated. In this manner, the subject CT value model data 44 is generated by matching the standard human CT value model data 41 with the actual subject 15 based on the two scanogram image data 42 and 43.

次に、ステップ202の被検体μモデル生成の工程では、ステップ201で生成した被検体CT値モデルデータ44のCT値を線減弱係数μに変換して、被検体15の3次元的μ値分布モデルを生成する。CT値から線減弱係数μへの変換は下記の如く行われる。   Next, in the process of generating the subject μ model in step 202, the CT value of the subject CT value model data 44 generated in step 201 is converted into a linear attenuation coefficient μ, and the three-dimensional μ value distribution of the subject 15 is converted. Generate a model. Conversion from the CT value to the linear attenuation coefficient μ is performed as follows.

CT値は実効エネルギー(通常60keVを使用)のX線に対する線減弱係数によって決定され、水=0、空気=1000、平均的な骨=1000として定義されている。今、被検体CT値モデルの位置xにおけるCT値をCTxとすると、その位置xにおける実効エネルギー(60keV)での線減弱係数μxは[数6]、[数7]で表される。 The CT value is determined by the linear attenuation coefficient for X-rays of effective energy (usually using 60 keV), and is defined as water = 0, air = 1000, average bone = 1000. Now, assuming that the CT value at the position x of the subject CT value model is CT x , the linear attenuation coefficient μ x at the effective energy (60 keV) at the position x is expressed by [Equation 6] and [Equation 7].

Figure 0004644292
Figure 0004644292

Figure 0004644292
ここで、μwは水の線源弱係数(=0.206cm-1)、μairは空気の線減弱係数(=0.00025cm-1)、μboneは骨の線減弱係数(=0.567cm-1、ただし密度1.8g/cm3の場合)である。
Figure 0004644292
Here, mu w is the source weak coefficient of water (= 0.206cm-1), μ air is linear attenuation coefficient of air (= 0.00025cm-1), μ bone is linear attenuation coefficient of bone (= 0.567cm-1 However, the density is 1.8 g / cm3).

次に、ステップ203の線量分布計算の工程では、図12に示す如く、ステップ111で求めた被検体に照射されたX線量のデータ(図12(a))とステップ202で求めた被検体μモデルデータ45(図12(b))を用いて、被検体15内の線量分布(図12(c))を計算する。このステップ203の計算では、被検体に照射されるX線のエネルーギースペクトルを考慮して、被検体15内でのX線の減衰を計算し、被検体15内の線量の空間的分布を計算する。被検体15に任意の方向からX線を照射した場合のX線の減弱量は、被検体の3次元的線減弱係数(μ)のモデルである被検体μモデルデータを用いることにより解析的に計算可能であり、このような計算手法は既に他の分野、例えば放射線治療計画装置などの分野でも行われている(参考文献1、稲邑清也、放射線治療計画システム、P.90〜92、P.113〜115、篠原出版、平成4月20日発行)。   Next, in the dose distribution calculation process of step 203, as shown in FIG. 12, data of the X-ray dose irradiated to the subject obtained in step 111 (FIG. 12 (a)) and the subject μ obtained in step 202 are obtained. Using the model data 45 (FIG. 12 (b)), the dose distribution in the subject 15 (FIG. 12 (c)) is calculated. In the calculation of step 203, the attenuation of X-rays in the subject 15 is calculated in consideration of the energy spectrum of the X-rays irradiated to the subject, and the spatial distribution of the dose in the subject 15 is calculated. . The amount of X-ray attenuation when subject 15 is irradiated with X-rays from any direction can be analyzed analytically by using subject μ model data that is a model of the subject's three-dimensional linear attenuation coefficient (μ). Such calculation methods have already been carried out in other fields, for example, in the field of radiation therapy planning devices (Reference 1, Kiyoya Inagi, Radiation therapy planning system, P. 90-92, P. 113-115, Shinohara Publishing, April 20, 1992).

X線の減衰の計算にあたっては、先ず、被検体15内の注目する位置xからX線源12に向かってX線が透過する実効距離δを計算する。実効距離は、X線が透過する媒質によって1/eに減衰する距離を1と定義される。X線のエネルギースペクトルを考慮しない場合、注目位置xからX線源12までに、組成iの媒質が実距離でdi(cm)あり、かつ組成iの線減弱係数がμi(cm-1)であるとき、実効距離δは[数8]で表される。 In calculating the attenuation of X-rays, first, an effective distance δ through which X-rays pass from the target position x in the subject 15 toward the X-ray source 12 is calculated. The effective distance is defined as 1 that is attenuated to 1 / e by the medium through which X-rays pass. When the X-ray energy spectrum is not taken into consideration, the medium of composition i has a real distance d i (cm) from the target position x to the X-ray source 12, and the linear attenuation coefficient of composition i is μ i (cm-1 ), The effective distance δ is expressed by [Equation 8].

Figure 0004644292
しかし、X線のエネルギースペクトルを考慮すると実効距離はX線のエネルギーに応じて異なる値となる。組成iのX線のエネルギーjに対する線減弱係数をμij(cm-1)とすると、X線のエネルギーjに対する実効距離δjは[数9]で表される。
Figure 0004644292
However, considering the X-ray energy spectrum, the effective distance varies depending on the X-ray energy. When the linear attenuation coefficient with respect to the energy j of the X-ray of the composition i is μ ij (cm−1), the effective distance δ j with respect to the energy j of the X-ray is expressed by [Equation 9].

Figure 0004644292
ここで、μijは組成iのエネルギーjのX線に対する線減弱系数(cm-1)、diは組成i中のX線の透過距離(cm)である。μijについては、X線のエネルギーjに応じて被検体15のμ値モデルから求める必要がある。
Figure 0004644292
Here, μ ij is the linear attenuation coefficient (cm −1) for X-rays of energy j of composition i , and d i is the transmission distance (cm) of X-rays in composition i . μ ij needs to be obtained from the μ value model of the subject 15 according to the energy j of the X-ray.

次に、被検体15内の注目位置xでの線量を計算する。X線源12から注目位置xまでの距離をrx(m)、距離1mにおける線量をI0(C/kg:Cはクーロン)とする。I0については、例えば実験的に求める。X線のエネルギースペクトル、つまりエネルギーjの成分比をSjとすると、注目位置xでの線量Ix(C/kg)は[数10]で表される。 Next, the dose at the target position x in the subject 15 is calculated. The distance from the X-ray source 12 to the target position x is r x (m), and the dose at a distance of 1 m is I0 (C / kg: C is Coulomb). For example, I0 is obtained experimentally. When the energy spectrum of X-rays, that is, the component ratio of energy j is S j , the dose I x (C / kg) at the target position x is expressed by [Equation 10].

Figure 0004644292
ここで、I0は線量(X線源12から単位距離における空中の線量、単位C/kg)、rxはX線源12から注目位置xまでの距離、単位m)、SjはX線のエネルギースペクトルである。
Figure 0004644292
Where I0 is the dose (dose in the air at a unit distance from the X-ray source 12, unit C / kg), r x is the distance from the X-ray source 12 to the target position x, unit m), and S j is the X-ray It is an energy spectrum.

[数10]による計算の結果、X線源12がある位置Q(Z、θ)にある場合の被検体15内の任意に位置xにおける線量が求められる。X線源12を1回転した場合の任意の位置xにおける線量は、X線源12の位置Q(Z、θ)を被検体15の周囲で回転させ、上記の線量を1回転分(θ=0〜2π)積算することによって得られる。上記の如き手順で、被検体15内の設定された各位置での計算を行うことにより、スライス位置Zでの被検体15内の線量分布を求めることができる。また、他のスライス位置における線量分布も同様に計算によって求めることができるので、被検体15の体軸方向について撮影部位の全領域の計算を進めることにより、被検体15内の3次元的な線量分布が得られる。   As a result of the calculation by [Equation 10], a dose at an arbitrary position x in the subject 15 when the X-ray source 12 is at a certain position Q (Z, θ) is obtained. The dose at an arbitrary position x when the X-ray source 12 is rotated once is obtained by rotating the position Q (Z, θ) of the X-ray source 12 around the subject 15 and the above dose for one rotation (θ = 0 to 2π) obtained by integrating. The dose distribution in the subject 15 at the slice position Z can be obtained by performing calculation at each set position in the subject 15 by the procedure as described above. In addition, dose distributions at other slice positions can be obtained by calculation in the same manner, so that the three-dimensional dose in the subject 15 can be calculated by proceeding with the calculation of the entire region of the imaging region in the body axis direction of the subject 15. Distribution is obtained.

1スライスの断層画像を撮影する場合(2次元の場合)には上記の計算によって被検体15内の線量の分布が精度良く求められるが、1回のスキャン中に複数スライスのCT撮影を行う場合(3次元の場合)には散乱X線の撮影を考慮しないと計算精度が低下する恐れがある。X線CT装置では、X線のエネルギーが100keV以下のオーダーであるため、散乱線としてはコンプトン(Compton)散乱のみ考慮すればよい(参考文献1参照)。このコンプトン散乱を考慮することにより、計算精度をより高くすることができる。   When taking a 1-slice tomographic image (in the case of 2D), the distribution of the dose in the subject 15 can be accurately obtained by the above calculation, but when performing CT imaging of multiple slices during one scan. In the case of (three-dimensional), there is a risk that the calculation accuracy may be lowered unless taking the scattered X-ray image. In the X-ray CT apparatus, since the energy of the X-ray is on the order of 100 keV or less, only Compton scattering needs to be considered as the scattered radiation (see Reference 1). Considering this Compton scattering, the calculation accuracy can be further increased.

上記の如き手順で、被検体15内の各スライス位置での各設定点における線量を計算することにより、計算された被検体15の3次元的線量分布が得られる。この被検体15の線量分布の計算結果はホストコンピュータ20に一時的に保存され、操作者などにとって見やすい図象、例えば図10に示したような図象で表示される。   By calculating the dose at each set point at each slice position in the subject 15 by the procedure as described above, the calculated three-dimensional dose distribution of the subject 15 is obtained. The calculation result of the dose distribution of the subject 15 is temporarily stored in the host computer 20, and is displayed in a diagram that is easy to see for an operator, for example, a diagram as shown in FIG.

次に、ステップ204の線量分布表示の工程ではステップ203の計算結果が表示装置25に表示される。本実施例での表示例としては、図10に示したような被検体15の断層面35内の線量分布(図10(a))又は被検体15の側面36の線量分布(図10(b))が上げられる。これらの図では、被検体15の臓器と線量分布を示す等線量線を重ねて表示しているので、一目で各臓器への被曝線量を認識することができるので、被検体15へのX線被曝を評価する上では有効である。   Next, in the dose distribution display process of step 204, the calculation result of step 203 is displayed on the display device 25. As a display example in the present embodiment, the dose distribution in the tomographic plane 35 of the subject 15 as shown in FIG. 10 (FIG. 10 (a)) or the dose distribution on the side surface 36 of the subject 15 (FIG. 10 (b) )) Is raised. In these figures, since the isodose lines indicating the organs and dose distribution of the subject 15 are displayed in an overlapping manner, the doses to each organ can be recognized at a glance, so the X-rays to the subject 15 It is effective in evaluating exposure.

次に、ステップ205の被曝線量判断の工程では、ステップ204で表示された被検体15内の線量分布の計算結果を操作者が見て、被検体15内の臓器へのX線被曝が過剰になる恐れがないかどうかを判断し、Yesと判断した場合にはステップ114のスキャン実行の工程に進み、スキャンを開始することになり、Noと判断した場合にはステップ109の管電流設定値入力の工程に戻り、管電流設定値の再入力、管電流パターンの再検討を行うことになる。   Next, in the step of determining the exposure dose in step 205, the operator looks at the calculation result of the dose distribution in the subject 15 displayed in step 204, and the X-ray exposure to the organ in the subject 15 is excessive. If yes, proceed to the scan execution process in step 114 and start scanning. If no, enter the tube current setting value in step 109. Returning to the process, the tube current set value is re-input and the tube current pattern is re-examined.

また、ステップ201での被検体15のCT値モデルデータの生成方法としては、上記の標準人体CT値モデルデータを使う方法以外に、過去に撮影した同じ被検体15のCT撮影データを使用する方法も実施可能である。この場合には、実際に同一の被検体15のCT値分布データを使用するために、標準人体CT値モデルデータ41の形状を補正する手順が必要なくなるという利点がある。しかし初回のCT撮影には適さないため、過去にCT撮影を行った被検体に関して、2回目以降のCT撮影を行う場合が対象となる。   Further, as a method of generating CT value model data of the subject 15 in step 201, in addition to the method of using the standard human CT value model data described above, a method of using CT imaging data of the same subject 15 taken in the past Can also be implemented. In this case, since the CT value distribution data of the same subject 15 is actually used, there is an advantage that a procedure for correcting the shape of the standard human CT value model data 41 is not necessary. However, since it is not suitable for the first CT imaging, the case where the CT imaging for the second time or later is performed on the subject that has been CT imaging in the past is targeted.

上記したように、CTスキャン前に被検体15の体内の線量分布をシミュレーション計算して、計算結果を、例えば図10に示す如く表示することにより、操作者は事前に撮影手技に応じた被検体の体内の線量分布を近似的に知ることが可能となる。   As described above, the dose distribution in the body of the subject 15 is calculated by simulation before the CT scan, and the calculation result is displayed as shown in FIG. 10, for example. It becomes possible to know the dose distribution in the body approximately.

この結果、例えば、単純に被検体の全ての組織に関して一律に被曝線量を減らすのではなく、骨髄や肺等の放射線感受性の高い組織に関しては特に被曝線量を低減し、逆に脂肪や筋肉等の放射線感受性の比較的低い組織には、画質が満足できる程度に被曝線量のレベルを維持するというような詳細な設定が可能となる。   As a result, for example, instead of simply reducing the dose uniformly for all tissues of the subject, the dose is reduced especially for tissues with high radiation sensitivity such as bone marrow and lungs. For tissues with relatively low radiation sensitivity, it is possible to make detailed settings such as maintaining the exposure dose level to the extent that image quality is satisfactory.

次に、図13を用いて、被検体をCT撮影する管電流の変化パターンの編集例を説明する。図13は、被検体のスキャノグラム画像上に管電流の変化パターンを重畳して表示したものである。図12において、スキャノグラム画像データ29aは体幹部のもの、管電流の変化パターンは、編集前の初期の管電流の変化パターン46aと編集後の修正された管電流の変化パターン46bである。   Next, an editing example of a change pattern of tube current for CT imaging of a subject will be described with reference to FIG. FIG. 13 shows a tube current change pattern superimposed on a scanogram image of a subject. In FIG. 12, the scanogram image data 29a is for the trunk, and the tube current change patterns are an initial tube current change pattern 46a before editing and a modified tube current change pattern 46b after editing.

この管電流の変化パターンの編集工程では、表示装置25の画面においてスキャノグラム画像29a上に表示された初期に設定された管電流の変化パターン46aに対し、スキャノグラム画像データ29aを参照しながら、また場合によっては被検体15の内部の照射線量分布を参照して、操作手段21によって修正を加えて、新しい管電流の変化パターン46bを編集する。この編集操作によって任意の部位の管電流の変化パターンを再設定する。   In the editing process of the tube current change pattern, with respect to the initially set tube current change pattern 46a displayed on the scanogram image 29a on the screen of the display device 25, while referring to the scanogram image data 29a, In some cases, the irradiation dose distribution inside the subject 15 is referred to, the operation means 21 is used for correction, and a new tube current change pattern 46b is edited. By this editing operation, the change pattern of the tube current at an arbitrary part is reset.

この編集操作において、例えば、自動的な管電流の変化パターンの設定では、横隔膜付近のように密度が大きく変化する領域では、管電流を平均的な値に設定するが、被曝線量が増えても画質を向上させる必要がある領域などでは、管電流を部分的に高く設定する。管電流の変化パターンは、上記の如くスキャン条件が設定されていれば、時間tのみの関数になるので、任意時刻の管電流の値を変化させることができる。図13の例では、初期の管電流の変化パターン46aに対し、肺の領域の管電流を少し低下させ、横隔膜の領域の管電流を少し増加させることで、修正後の管電流の変化パターン46bに編集している。   In this editing operation, for example, in the automatic tube current change pattern setting, the tube current is set to an average value in a region where the density changes greatly as in the vicinity of the diaphragm, but even if the exposure dose increases. In areas where image quality needs to be improved, the tube current is set to be partially high. Since the change pattern of the tube current is a function of only the time t if the scan condition is set as described above, the value of the tube current at an arbitrary time can be changed. In the example of FIG. 13, the tube current change pattern 46b after correction is obtained by slightly reducing the tube current in the lung region and slightly increasing the tube current in the diaphragm region with respect to the initial tube current change pattern 46a. I am editing it.

図14は、管電流と被検体の厚さ(X線透過長に相当)との関係を示した図である。上記のステップ110の工程の説明では、管電流の最大値、最小値と被検体の厚さの最大値、最小値とを一致させて、両者間で線形の関係を持つものとして説明したが、両者の関係については操作者の設定により非線形の関係を持つようにすることが可能である。図14に示した管電流Iと被検体の厚さTとの関係は、[数11]で表される。   FIG. 14 is a diagram showing the relationship between the tube current and the thickness of the subject (corresponding to the X-ray transmission length). In the description of the step 110 above, the maximum value and the minimum value of the tube current are matched with the maximum value and the minimum value of the thickness of the object, and it is described as having a linear relationship between the two. Regarding the relationship between the two, it is possible to have a non-linear relationship by setting of the operator. The relationship between the tube current I and the thickness T of the subject shown in FIG. 14 is expressed by [Equation 11].

Figure 0004644292
ここで、Imax、Iminは管電流の最大値と最小値、Tmax、Tminは被検体の厚さの最大値と最小値、γは定数である。γについては以下ガンマと呼ぶことにする。
Figure 0004644292
Here, I max and I min are the maximum and minimum values of the tube current, T max and T min are the maximum and minimum values of the thickness of the subject, and γ is a constant. Hereinafter, γ will be referred to as gamma.

図14において、グラフ50はガンマ=1の場合で、管電流と被検体の厚さとの関係は線形であり、グラフ51はガンマ<1の場合、グラフ52はガンマ>1の場合で、共に管電流と被検体の厚さとの関係は非線形である。図14の場合、ガンマの値を決めることによって管電流と被検体の厚さとの関係が一義的に決まるので、[数11]のような関係式を装置に組み込んでおくことにより、操作者は操作手段21からガンマの値を入力することによって、図14の如く、管電流Iと被検体の厚さTとの関係を変化させることができる。また、実際の操作にあたっては、例えば操作者が特別な設定をしない初期設定では管電流の最大値、最小値と被検体の厚さの最大値、最小値を一致させて、線形の関係を持たせることにし、操作者の設定によりガンマを入力することで、非線形の関係を持たせることができる。図14において、ガンマ=1を基準にした場合、ガンマ>1のときは被検体の被曝低減を重視する場合とみられ、ガンマ<1のときは画質を重視する場合とみられる。   In FIG. 14, the graph 50 shows a case where gamma = 1, and the relationship between the tube current and the thickness of the subject is linear, the graph 51 shows a case where gamma <1, and the graph 52 shows a case where gamma> 1, both in the tube. The relationship between the current and the thickness of the subject is non-linear. In the case of FIG. 14, the relationship between the tube current and the thickness of the subject is uniquely determined by determining the value of gamma. Therefore, by incorporating a relational expression such as [Equation 11] into the apparatus, the operator can By inputting the value of gamma from the operating means 21, the relationship between the tube current I and the thickness T of the subject can be changed as shown in FIG. In actual operation, for example, in an initial setting in which no special setting is made by the operator, the maximum value and the minimum value of the tube current are matched with the maximum value and the minimum value of the thickness of the subject, and a linear relationship is established. By inputting gamma according to the operator's setting, a non-linear relationship can be obtained. In FIG. 14, when gamma = 1 is used as a reference, when gamma> 1, it is considered that the reduction in exposure of the subject is emphasized, and when gamma <1, it is considered that the image quality is emphasized.

10 ガントリ、11 CTスキャナ(スキャナ)、12 X線源、13 X線検出器(検出器)、14 X線ビーム(X線)、15 被検体、16 スキャナ駆動手段、17 スキャナ角度検出手段、18 テーブル、19 コリメータ、20 ホストコンピュータ、21 操作手段、22 高電圧発生装置、23 管電流設定手段、24 画像処理装置、25 表示装置、26 スキャノグラム解析手段、29、29a スキャノグラム画像、30 X線透過長モデル(楕円形モデル)、31、31a、31b、31c、31d、47a、47b 管電流の変化パターン、35 被検体横断面、36 被検体側面、38a、38b、38c、39a、39b、39c、39d、46a、46b、46c 等線量線、41 標準人体CT値モデルデータ、42 計算スキャノグラム画像データ、43 実測スキャノグラム画像データ、44 被検体CT値モデルデータ、45 被検体μモデルデータ、51、52、53 グラフ   10 Gantry, 11 CT scanner (scanner), 12 X-ray source, 13 X-ray detector (detector), 14 X-ray beam (X-ray), 15 Subject, 16 Scanner drive means, 17 Scanner angle detection means, 18 Table, 19 Collimator, 20 Host computer, 21 Operating means, 22 High voltage generator, 23 Tube current setting means, 24 Image processing device, 25 Display device, 26 Scanogram analysis means, 29, 29a Scanogram image, 30 X-ray transmission length Model (elliptical model), 31, 31a, 31b, 31c, 31d, 47a, 47b Change pattern of tube current, 35 Cross section of subject, 36 Side surface of subject, 38a, 38b, 38c, 39a, 39b, 39c, 39d 46a, 46b, 46c isodose line, 41 standard human CT value model data, 42 calculated scanogram image data, 43 measured scanogram image data, 44 subject CT value model data, 45 subject μ model data, 51, 52, 53 Graph

Claims (1)

被検体にX線を照射するX線源と、
X線源と対向配置され被検体の透過X線量を検出するX線検出器と、
被検体のスキャノグラム像を表示する表示装置と、
表示されたスキャノグラム像上に被検体の断層像の撮影位置を設定する撮影位置設定手段と、
X線源を回転させながらそのX線源からX線を照射し、X線検出器により検出された被検体の透過X線量から断層像を再構成する断層像再構成手段と、
再構成された断層像を表示する表示装置とを備え、
被検体の撮影位置に応じたX線をX線源に発生させるための管電流パターンを、管電流値と撮影位置とを対比させるためにスキャノグラム像と重畳または並置して表示装置に表示させ、撮影位置を参照しながら前記管電流パターン上の特定の撮影位置における当該管電流値を編集できるようにし、編集された当該管電流値に基づいて管電流パターンを表示できるようにしたX線CT装置において、
管電流パターンに応じた線量を、前記管電流パターンを時間で積分して計算する線量計算手段をさらに備え、計算された線量をスキャノグラム像と並置して、又は重ねて表示装置に表示できるようにしたことを特徴とするX線CT装置。
An X-ray source for irradiating the subject with X-rays;
An X-ray detector disposed opposite to the X-ray source to detect the transmitted X-ray dose of the subject;
A display device for displaying a scanogram image of the subject;
An imaging position setting means for setting an imaging position of a tomographic image of the subject on the displayed scanogram image;
A tomographic image reconstruction means for irradiating an X-ray from the X-ray source while rotating the X-ray source and reconstructing a tomographic image from the transmitted X-ray dose of the subject detected by the X-ray detector;
E Bei a display device for displaying the reconstructed tomographic image,
A tube current pattern for generating X-rays corresponding to the imaging position of the subject on the X-ray source is displayed on the display device so as to be superimposed or juxtaposed with the scanogram image in order to compare the tube current value and the imaging position, X-ray CT apparatus capable of editing the tube current value at a specific imaging position on the tube current pattern while referring to the imaging position, and displaying the tube current pattern based on the edited tube current value In
Dose calculating means for calculating the dose according to the tube current pattern by integrating the tube current pattern with time so that the calculated dose can be displayed on the display device in juxtaposition with or overlapping the scanogram image. An X-ray CT apparatus characterized by that.
JP2009099614A 2009-04-16 2009-04-16 X-ray CT apparatus and image display method thereof Expired - Lifetime JP4644292B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2009099614A JP4644292B2 (en) 2009-04-16 2009-04-16 X-ray CT apparatus and image display method thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009099614A JP4644292B2 (en) 2009-04-16 2009-04-16 X-ray CT apparatus and image display method thereof

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2001067150A Division JP4532005B2 (en) 2001-03-09 2001-03-09 X-ray CT apparatus and image display method thereof

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2010174976A Division JP5027909B2 (en) 2010-08-04 2010-08-04 X-ray CT system

Publications (2)

Publication Number Publication Date
JP2009154021A JP2009154021A (en) 2009-07-16
JP4644292B2 true JP4644292B2 (en) 2011-03-02

Family

ID=40958630

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009099614A Expired - Lifetime JP4644292B2 (en) 2009-04-16 2009-04-16 X-ray CT apparatus and image display method thereof

Country Status (1)

Country Link
JP (1) JP4644292B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5479069B2 (en) * 2009-12-17 2014-04-23 株式会社東芝 X-ray CT apparatus and control method of X-ray CT apparatus
JP7055709B2 (en) * 2018-06-27 2022-04-18 キヤノンメディカルシステムズ株式会社 X-ray CT device and imaging planning device

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6152860A (en) * 1984-08-24 1986-03-15 株式会社 日立メデイコ X-ray ct apparatus
JPH05220135A (en) * 1992-02-14 1993-08-31 Toshiba Corp Ct device
JPH07124152A (en) * 1993-11-01 1995-05-16 Toshiba Corp X-ray ct scanner
JPH07204196A (en) * 1994-01-24 1995-08-08 Hitachi Medical Corp X-ray ct system
JPH08168486A (en) * 1994-08-03 1996-07-02 General Electric Co <Ge> Method of modulating dose of x-ray beam in x-ray ct system
JPH08206107A (en) * 1994-09-06 1996-08-13 General Electric Co <Ge> X-ray volume reducing method
JPH09199292A (en) * 1995-11-20 1997-07-31 General Electric Co <Ge> System for modulating x-ray tube current supplied to x-ray source
JP2002263097A (en) * 2001-03-09 2002-09-17 Hitachi Medical Corp Radiographic tomograph

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6152860A (en) * 1984-08-24 1986-03-15 株式会社 日立メデイコ X-ray ct apparatus
JPH05220135A (en) * 1992-02-14 1993-08-31 Toshiba Corp Ct device
JPH07124152A (en) * 1993-11-01 1995-05-16 Toshiba Corp X-ray ct scanner
JPH07204196A (en) * 1994-01-24 1995-08-08 Hitachi Medical Corp X-ray ct system
JPH08168486A (en) * 1994-08-03 1996-07-02 General Electric Co <Ge> Method of modulating dose of x-ray beam in x-ray ct system
JPH08206107A (en) * 1994-09-06 1996-08-13 General Electric Co <Ge> X-ray volume reducing method
JPH09199292A (en) * 1995-11-20 1997-07-31 General Electric Co <Ge> System for modulating x-ray tube current supplied to x-ray source
JP2002263097A (en) * 2001-03-09 2002-09-17 Hitachi Medical Corp Radiographic tomograph

Also Published As

Publication number Publication date
JP2009154021A (en) 2009-07-16

Similar Documents

Publication Publication Date Title
JP4532005B2 (en) X-ray CT apparatus and image display method thereof
JP4937927B2 (en) X-ray CT apparatus and imaging condition determination method in X-ray CT apparatus
US8483363B2 (en) Movable wedge for improved image quality in 3D X-ray imaging
US7818045B2 (en) Method for reduced-artifact radiological 3D imaging, medical imaging device and method for creating a therapy plan
JP6066596B2 (en) Method and system for scatter correction in X-ray imaging
CN107252318A (en) Acquisition protocols valuator device
JP2007181623A (en) X-ray ct apparatus
JP4468352B2 (en) Reconstruction of local patient dose in computed tomography
JP2006312047A (en) Method of calculating x-ray attenuation in orthogonal cross direction of object, computer tomographic photography device, and correction module therefor
JP4429694B2 (en) X-ray CT system
JP4554185B2 (en) X-ray CT system
JP2008220452A (en) Method and system for radiography having radiation profile setting based on organ
JP5027909B2 (en) X-ray CT system
WO2006090321A1 (en) Determination of the coverage of a ct scan
US10346957B2 (en) Method for image improvement of image data from a dental image generation system
JP4644292B2 (en) X-ray CT apparatus and image display method thereof
KR20210146384A (en) A medical image processing apparatus, a storage medium, a medical apparatus, and a treatment system
CN107809954B (en) Display of depth position of computer tomography slice image relative to object to be imaged
US11937970B2 (en) System and method for calibrating a camera feature detection system of an x-ray system
JP5384293B2 (en) X-ray CT system
JP7258474B2 (en) X-ray CT device and radiotherapy system
JP4056968B2 (en) X-ray CT apparatus and image processing method

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090507

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20090507

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20090507

A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20090609

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090615

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090812

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20091028

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20091225

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20100308

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100428

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20100607

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100804

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

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20101203

R150 Certificate of patent or registration of utility model

Ref document number: 4644292

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

Year of fee payment: 3

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

EXPY Cancellation because of completion of term