JPH05168621A - X-ray ct system - Google Patents

X-ray ct system

Info

Publication number
JPH05168621A
JPH05168621A JP3355856A JP35585691A JPH05168621A JP H05168621 A JPH05168621 A JP H05168621A JP 3355856 A JP3355856 A JP 3355856A JP 35585691 A JP35585691 A JP 35585691A JP H05168621 A JPH05168621 A JP H05168621A
Authority
JP
Japan
Prior art keywords
noise
data
measuring
measurement data
measuring data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP3355856A
Other languages
Japanese (ja)
Inventor
Shinichi Uda
晋一 右田
Tetsuo Nakazawa
哲夫 中澤
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 JP3355856A priority Critical patent/JPH05168621A/en
Publication of JPH05168621A publication Critical patent/JPH05168621A/en
Pending legal-status Critical Current

Links

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)

Abstract

PURPOSE:To constitute a measuring instrument for handling a minute signal such as an X-ray CT system, etc., so that it is prevented that a minute disturbance noise is mixed into measuring data, and this noise appears as a beat noise on the measuring data by synchronizing with a measuring sampling period, and a moire pattern artifact generated on an image due to this noise is reduced. CONSTITUTION:A weighted mean processing between measuring data of every ch and the view direction being adjacent to this ch or ch direction data is executed, and by replacing this value with the measuring data and forming it in the same shape as the original measuring data, all the same image processing as a conventional one is executed. Accordingly, a mixed noise appears as an oblique line-like noise (beat noise) on the measuring data by synchronizing a measuring sampling period. When a suitable weighted mean processing of the measuring data and the adjacent data is executed, this beat noise is scattered suitably and subjected to smoothing, and a moire pattern on an image comes not to be generated.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】本発明は、計測データへのビート
パターンの混入除去をはかるX線CT装置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an X-ray CT apparatus for removing a beat pattern from measurement data.

【0002】[0002]

【従来の技術】X線管装置やX線スキャナの機械的な振
動ノイズやX線発生装置のスィッチングノイズ、交流電
圧源からの誘導ノイズが、X線CT装置の計測データD
ijに混入することがある。ここでiとはX線検出器のチ
ャンネル番号,jとは被検体の円周方向のビュー方向番
号である。振動ノイズや誘導ノイズの混入例を図3に示
す。図3(イ)は横軸に計測時間の推移、縦軸は信号振
幅を示す。時間軸上の黒マル印は、サンプル点を示す。
図3(イ)では、微小高周波成分の振幅部分N1が振動
ノイズやN1誘導ノイズより成る外部ノイズであり、正
規の計測信号にこの外乱ノイズ4が乗っている様子が示
されている。この図3(イ)によれば、外乱ノイズN1
は時間推移に対して一様に乗っていることがわかる。
2. Description of the Related Art Mechanical vibration noise of an X-ray tube device or X-ray scanner, switching noise of an X-ray generator, and induction noise from an AC voltage source are measured data D of an X-ray CT device.
May be mixed in ij . Here, i is the channel number of the X-ray detector, and j is the view direction number in the circumferential direction of the subject. An example of mixing vibration noise and induction noise is shown in FIG. In FIG. 3A, the horizontal axis represents the transition of the measurement time, and the vertical axis represents the signal amplitude. Black circle marks on the time axis indicate sample points.
FIG. 3A shows that the amplitude portion N 1 of the minute high-frequency component is external noise including vibration noise and N 1 induction noise, and the disturbance noise 4 is superimposed on the regular measurement signal. According to FIG. 3A, the disturbance noise N 1
It can be seen that is uniformly applied to the time transition.

【0003】一方、図3(b)は、各サンプル点に外乱
ノイズ周期が一致した例である。図では全サンプル点一
致としているが、2サンプル点や3サンプル点毎の一致
もありうる。外乱ノイズN2は時間推移の平均でみれば
大きな振幅ではないが、外乱ノイズのピークがサンプル
点に一致すると、サンプル点が計測点である故に計測信
号に対する外乱ノイズN2の割合は大きくなる。従っ
て、平均的にみては微小な外乱ノイズN2であっても、
ノイズピーク値がサンプル点に一致すると、無視できな
いものとなる。
On the other hand, FIG. 3 (b) is an example in which the disturbance noise period matches each sample point. In the figure, all sample points match, but there may be matching at every 2 sample points or every 3 sample points. The disturbance noise N 2 does not have a large amplitude on the average of the time transition, but when the peak of the disturbance noise coincides with the sample point, the ratio of the disturbance noise N 2 to the measurement signal increases because the sample point is the measurement point. Therefore, even if the disturbance noise N 2 is small on average,
When the noise peak value matches the sample point, it cannot be ignored.

【0004】図4(イ)には、図3(イ)の如き計測で
のビュー方向番号jとチャンネル番号iとで定義された
計測データDijの様子を示す(Dijはバッファメモリに
格納されている)。図4(イ)によれば、ランダムノイ
ズは黒マル点の如くランダムに配置される。図4(ロ)
は、図3(ロ)の如き計測での計測データDijの様子を
示す図である。図4(ロ)によれば、サンプル周期と外
乱ノイズ周期が特定の周期で一致した例であり、これに
よって計測データDij中には、斜線上のノイズ(ビート
パターンノイズ)が現われる。
FIG. 4A shows a state of the measurement data D ij defined by the view direction number j and the channel number i in the measurement as shown in FIG. 3A, where D ij is stored in the buffer memory. Has been). According to FIG. 4A, the random noises are randomly arranged like black circle points. Figure 4 (b)
FIG. 4 is a diagram showing a state of measurement data D ij in the measurement shown in FIG. FIG. 4B shows an example in which the sampling period and the disturbance noise period coincide with each other at a specific period, which causes noise (beat pattern noise) on the diagonal line in the measurement data D ij .

【0005】図4(ハ)は、図4(イ)のランダムノイ
ズを持つ計測システムで計測が行われこの時のデータを
再構成した場合のCT像である。このCT像では、ラン
ダムノイズが大きいとノイズが多くなるといった影響は
あるが、図形上特に明確に現われない。図4(ニ)は、
図4(ロ)のビートパターンを持つ計測システムで計測
されたデータを再構成した場合のCT像である。図4
(ニ)によれば、渦巻き状のアーチファクト(モアレパ
ターンアーチアァクトと呼ばれる)が発生する。
FIG. 4C is a CT image when the measurement is performed by the measurement system having random noise of FIG. 4A and the data at this time is reconstructed. In this CT image, if the random noise is large, the noise is increased, but it does not appear clearly on the figure. Figure 4 (d) shows
It is a CT image when the data measured by the measurement system having the beat pattern of FIG. 4B is reconstructed. Figure 4
According to (d), a spiral artifact (called a moire pattern artifact) occurs.

【0006】従来、ビートパターンノイズ除去のために
は、各ユニットの共振周波数・計測周期・X線インバー
タ周期が同期しない様に各周波数の設定をずらしたり、
検出器の振動防止用装置の装着や検出回路とノイズ源と
の距離を離したり、シールド処理などの機械的処理の他
に、各ユニット間I/Fフィルター等の電気処理を併用
していた。
Conventionally, in order to remove beat pattern noise, the setting of each frequency is shifted so that the resonance frequency of each unit, the measurement cycle, and the X-ray inverter cycle are not synchronized.
In addition to mounting a device for preventing vibration of the detector, increasing the distance between the detection circuit and the noise source, and mechanical processing such as shield processing, electrical processing such as I / F filter between units has been used together.

【0007】[0007]

【発明が解決しようとする課題】上記従来の対策である
装置の小型化とX線装置の高出力に伴い通常な方法だけ
では完全にノイズを除去するのが困難な状況になってい
る。本発明は従来技術を行っても除去できないビートノ
イズが混入した計測データから画像再構成して発生する
モアレパターンアーチファクトをなくし良好なCT画像
を得るX線CT装置を提供することを目的としている。
Due to the downsizing of the apparatus and the high output of the X-ray apparatus, which are the above-mentioned conventional measures, it is difficult to completely remove the noise only by the usual method. SUMMARY OF THE INVENTION It is an object of the present invention to provide an X-ray CT apparatus that eliminates moire pattern artifacts generated by image reconstruction from measurement data mixed with beat noise that cannot be removed by the conventional technique and obtains a good CT image.

【0008】本発明は、外部から計測データDijに混入
する外乱ノイズ周期に同期しないように、ビュー方向番
号(又はチャンネル方向番号)に沿っての計測データ間
で加重平均を行い、この加重平均値をビュー方向番号
(又はチャンネル方向番号)に沿っての計測データとし
て置き換えることとした。
The present invention performs weighted averaging between measurement data along the view direction number (or channel direction number) so as not to be synchronized with the disturbance noise period mixed into the measurement data D ij from the outside, and the weighted average is calculated. It was decided to replace the value as the measurement data along the view direction number (or channel direction number).

【0009】[0009]

【作用】本発明によれば、計測データ間で、ビュー方向
番号(又はチャンネル方向番号)に沿って加重平均値を
とり、これを計測データに代って置き換えるようにす
る。かくして、外乱ノイズによるビートパターンノイズ
はスムージングされることとなり実質上の除去をはか
る。
According to the present invention, a weighted average value is taken along the view direction number (or channel direction number) between measurement data, and this is replaced in place of the measurement data. Thus, the beat pattern noise due to the disturbance noise is smoothed and is effectively removed.

【0010】[0010]

【実施例】図1は、計測データDijのビュー方向加重平
均を行う実施例図を示す。データ群1としては、ビュー
番号j、チャンネル番号iとする計測データDijであ
り、ビュー区間は、0°〜360°(又はハーフスキャ
ンでは0°〜180°)であり、これをビューピッチで
区分けし、jで番号付けている。チャンネル番号iと
は、多チャンネルX線検出器の検出器番号を示す。図で
は、チャンネル番号iはi=1,2,……,n,ビュー
番号jはj=1,2,……,Nとしている。
EXAMPLE FIG. 1 shows an example diagram in which a view direction weighted average of measurement data D ij is performed. The data group 1 is the measurement data D ij with the view number j and the channel number i, and the view section is 0 ° to 360 ° (or 0 ° to 180 ° in half scan), which is the view pitch. It is classified and numbered by j. The channel number i indicates the detector number of the multi-channel X-ray detector. In the figure, the channel number i is i = 1, 2, ..., N, and the view number j is j = 1, 2 ,.

【0011】データ群2は、データ群1をビュー方向加
重平均して得た加重平均データ群Sijである。この加重
平均処理及び計測データへの置換は以下に従う。 (1)、i=1……,n−1の場合。
The data group 2 is a weighted average data group S ij obtained by weighting the data group 1 in the view direction. The weighted average processing and the replacement with the measured data are as follows. (1), when i = 1 ..., N-1.

【数1】Sij=a×Dij+b×Di+1,j ここで、a+b=1.0 aとbとは加重平均の重みである。 (2)、i=nの場合。## EQU1 ## S ij = a × D ij + b × D i + 1, j where a + b = 1.0 a and b are weighted average weights. (2) When i = n.

【数2】Sij=a×Dnj+b×D1j ここで、a+b=1.0 重みaとbとは、以下の考え方で選ぶ。即ち、従来の加
重平均を行わない例はa=1、b=0の例であり、aを
小さくしてbを大きくすることで外乱ノイズが低減され
ることがわかる。但し、あまりbを大きくすると、画像
のぼけが大きくなるとの欠点がある。こうしたノイズを
少なくする目的とぼけをどこまで許すかとの2点からa
とbとは決定される。以上の実施例では、計測データ数
は、データ群1と2とでは同じである。
## EQU2 ## S ij = a × D nj + b × D1j where a + b = 1.0 The weights a and b are selected according to the following concept. That is, the conventional example in which the weighted average is not performed is an example in which a = 1 and b = 0, and it is understood that the disturbance noise is reduced by decreasing a and increasing b. However, if b is set too large, there is a drawback that the blur of the image becomes large. From the two points of the purpose of reducing such noise and how much blur is allowed
And b are determined. In the above embodiment, the number of measurement data is the same in the data groups 1 and 2.

【0012】また、チャンネル方向の加重平均処理を行
う場合は、j=2……,N−1の範囲で
When the weighted average processing in the channel direction is performed, in the range of j = 2 ...

【数3】Sij=a×Dij+b×Di,j+1 の処理を行えばよい。j=1、i=nのデータは処理し
ないか、処理ずみチャンネルデータとの差が大きくなら
ない程度の補正を行えばよい。また、本例では隣合うデ
ータ間の処理の例を示したが、数chまたは数ビューに
わたるデータ間の処理も可能であるし、ch方向にわた
るマトリックス状の加重平均処理値に置き換えることも
できる。
## EQU3 ## The process of S ij = a × D ij + b × D i, j + 1 may be performed. The data of j = 1 and i = n may not be processed, or may be corrected to such an extent that the difference from the processed channel data does not become large. In addition, although an example of processing between adjacent data is shown in this example, processing between data over several channels or several views is also possible, and the weighted average processing value in matrix form over the ch direction can be replaced.

【0013】また、本例はオリジナル計測データと置き
換えデータが同じch数、同じビュー数の形を示した
が、計測ビュー数が多い計測の場合は図2に示すように
数ビュー毎の平均値に置き換えて計測データを圧縮した
後に(数1)……,(数3)の処理を追加しても効果が
ある。また計測データはオフセット補正と感度補正済み
後の場合も、補正前でも良い。ここでaとbの重み付け
は任意に変えることができ、bの値が大きいほど画像の
ぼけが大きくなると先に述べたがこのぼけは周辺chデ
ータ部ほど顕著になるため図5(イ)、(ロ)に示すよ
うに周辺chほどbの重みを小さくして画像のぼけを少
なくすることも可能である。本例では双曲線の形状を示
したがこの形は周辺から中心部に徐々に変化していれば
よく、中心部Hの重みも任意に変えられる事はいうまで
もない。上記の補正を行った後に通常の画像再構成を行
って画像を得ることになるがX線CT装置による再構成
法等の詳細については省略する。
Also, in this example, the original measurement data and the replacement data have the same number of channels and the same number of views, but in the case of measurement with a large number of measurement views, as shown in FIG. It is also effective to add the processing of (Equation 1) ..., (Equation 3) after compressing the measurement data by replacing with. Further, the measurement data may be after offset correction and sensitivity correction, or before correction. Here, the weighting of a and b can be arbitrarily changed, and the blur of the image becomes larger as the value of b becomes larger. However, since this blur becomes more remarkable in the peripheral ch data portion, FIG. As shown in (b), it is also possible to reduce the blur of the image by reducing the weight of b for the peripheral ch. In this example, the shape of the hyperbola is shown, but it is needless to say that this shape may be gradually changed from the periphery to the central portion, and the weight of the central portion H can be arbitrarily changed. After performing the above correction, an ordinary image reconstruction is performed to obtain an image, but details of the reconstruction method by the X-ray CT apparatus and the like will be omitted.

【0014】以上の各種の加重平均処理は、コンピュー
タによって実施する。この場合、コンピュータは、計測
データDijを格納したバッファに対して、アクセスを行
い、aとbとを与えて加重平均処理を行い、この結果得
たSijをi,jでの計測データとし、これを別のバッフ
ァに次々に格納してゆく。
The above various weighted average processing is executed by a computer. In this case, the computer accesses the buffer storing the measurement data D ij , gives a and b, performs the weighted averaging process, and sets the resulting S ij as the measurement data at i, j. , Store it in another buffer one after another.

【0015】尚、加重平均処理を行った場合、計測デー
タのピーク状ノイズも小さくなるため、ストリークアー
チファクト状のノイズを併せて低減できる利点を持つ。
When the weighted averaging process is performed, the peak noise of the measurement data is also reduced, which has the advantage that streak artifact noise can be reduced.

【0016】本実施例は、X線CT装置の例としたが、
二次元的な定義可能で且つ時系列でデータを計測するよ
うな計測系の計測データに対してはすべて適用可能であ
る。
In this embodiment, the X-ray CT apparatus is used as an example.
It can be defined two-dimensionally and is applicable to all measurement data of a measurement system that measures data in time series.

【0017】[0017]

【発明の効果】本発明によれば、計測データに混入する
計測サンプリングに同期したビートノイズを押さえ、画
像上に発生するモアレパターン状アーチファクトをなく
すことができる。更に加重平均処理により計測データの
ピーク状ノイズも小さくなるため、ストリークアーチフ
ァクト状のノイズも低減し良好な画像を得ることができ
る。
According to the present invention, it is possible to suppress beat noise that is mixed with measurement data and that is synchronized with measurement sampling, and eliminate moire pattern-like artifacts generated on an image. Further, since the peak noise of the measurement data is reduced by the weighted averaging process, streak artifact noise is also reduced and a good image can be obtained.

【図面の簡単な説明】[Brief description of drawings]

【図1】本発明のビュー方向加重平均処理による計測デ
ータ置き換えの実施例図である。
FIG. 1 is a diagram illustrating an example of replacement of measurement data by a view direction weighted average processing according to the present invention.

【図2】本発明のビュー方向加重平均処理による計測デ
ータ置き換えの他の実施例図である。
FIG. 2 is a diagram illustrating another embodiment of replacement of measurement data by the view direction weighted average processing according to the present invention.

【図3】外乱ノイズ混入を示す図である。FIG. 3 is a diagram showing mixing of disturbance noise.

【図4】外乱ノイズ混入によるCT画像への影響の説明
図である。
FIG. 4 is an explanatory diagram of an influence on a CT image due to mixing of disturbance noise.

【図5】加重平均時重み係数の関係を示す図である。FIG. 5 is a diagram showing a relationship between weighted average weighting factors.

【符号の説明】[Explanation of symbols]

1 計測データ群Dij 2 加重平均処理後の置換計測データ群Sij 1 Measurement data group D ij 2 Replacement measurement data group S ij after weighted average processing

Claims (1)

【特許請求の範囲】[Claims] 【請求項1】 X線検出器チャンネル番号iとビュー方
向番号jとで定義された計測データDijを用いて再構成
するX線CT装置において、外部から計測データDij
混入する外乱ノイズ周期に同期しないように、ビュー方
向番号(又はチャンネル方向番号)に沿っての計測デー
タ間で加重平均を行い、この加重平均値をビュー方向番
号(又はチャンネル方向番号)に沿っての計測データと
して置き換えこの置き換え後の計測データから再構成す
るX線CT装置。
1. In an X-ray CT apparatus for reconstructing using measurement data D ij defined by an X-ray detector channel number i and a view direction number j, a disturbance noise period mixed into the measurement data D ij from the outside. In order not to synchronize with, the weighted average is performed between the measurement data along the view direction number (or channel direction number), and this weighted average value is replaced as the measurement data along the view direction number (or channel direction number). An X-ray CT device that reconstructs from the measurement data after this replacement.
JP3355856A 1991-12-20 1991-12-20 X-ray ct system Pending JPH05168621A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP3355856A JPH05168621A (en) 1991-12-20 1991-12-20 X-ray ct system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP3355856A JPH05168621A (en) 1991-12-20 1991-12-20 X-ray ct system

Publications (1)

Publication Number Publication Date
JPH05168621A true JPH05168621A (en) 1993-07-02

Family

ID=18446082

Family Applications (1)

Application Number Title Priority Date Filing Date
JP3355856A Pending JPH05168621A (en) 1991-12-20 1991-12-20 X-ray ct system

Country Status (1)

Country Link
JP (1) JPH05168621A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007275125A (en) * 2006-04-03 2007-10-25 Ge Medical Systems Global Technology Co Llc X-ray ct equipment and its x-ray ct image reconstruction method, and x-ray ct image photographic method

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007275125A (en) * 2006-04-03 2007-10-25 Ge Medical Systems Global Technology Co Llc X-ray ct equipment and its x-ray ct image reconstruction method, and x-ray ct image photographic method

Similar Documents

Publication Publication Date Title
US6115487A (en) Correction algorithm for bone-induced spectral artifacts in computed tomograph imaging
JP3342414B2 (en) Video quality degradation measuring method and apparatus
EP0407213B1 (en) Image processing apparatus
JPH0731606A (en) Magnetic resonance tomographic camera
JP2955723B2 (en) Imaging device
US4136388A (en) Data acquisition system for computed tomography
JPH08336522A (en) Ct scanner and production of ct scanner picture
King et al. Use of a nonstationary temporal Wiener filter in nuclear medicine
US4752879A (en) Method and apparatus for medical imaging
JPH05168621A (en) X-ray ct system
US4654796A (en) System for computing correction factors for sensitivity differences in detector channels of an X-ray tomography apparatus
JP2001276055A (en) Method for reducing linear artifact in ct image and device for performing this method
JPH04233817A (en) Method for low-pass filtering and device for executing the method
US5761267A (en) Methods and apparatus for simplified filtering of scan data in an imaging system
US5383032A (en) Image processing system
CN115063503B (en) Magnetic resonance image convolution artifact removing method, device, equipment and medium
JP3583836B2 (en) X-ray CT apparatus with view number conversion means
JP3884929B2 (en) Radiation image acquisition apparatus and design method
CN1320500C (en) Method of treating radiation image sequence obtained from dynamic image converter
JP3657689B2 (en) X-ray CT system
SU1192164A1 (en) Device or processing optical images
JP3779398B2 (en) Projection data correction method and X-ray CT apparatus
JPH0435466A (en) Picture signal processor
JP2996086B2 (en) Ghost removal device
JP2881903B2 (en) Image reading device