JPH119589A - X-ray ct tomograph and image recomposing method - Google Patents
X-ray ct tomograph and image recomposing methodInfo
- Publication number
- JPH119589A JPH119589A JP12051198A JP12051198A JPH119589A JP H119589 A JPH119589 A JP H119589A JP 12051198 A JP12051198 A JP 12051198A JP 12051198 A JP12051198 A JP 12051198A JP H119589 A JPH119589 A JP H119589A
- Authority
- JP
- Japan
- Prior art keywords
- image
- divided
- reconstructed
- reconstructed images
- ray
- 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
Links
- 238000000034 method Methods 0.000 title claims description 24
- 230000015654 memory Effects 0.000 claims abstract description 46
- 238000005259 measurement Methods 0.000 claims description 17
- 239000000470 constituent Substances 0.000 claims description 5
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 11
- 238000012937 correction Methods 0.000 abstract description 4
- 238000007792 addition Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000005055 memory storage Effects 0.000 description 2
- 101000741289 Homo sapiens Calreticulin-3 Proteins 0.000 description 1
- 101000969621 Homo sapiens Monocarboxylate transporter 12 Proteins 0.000 description 1
- 102100021444 Monocarboxylate transporter 12 Human genes 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000010562 histological examination Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Medical Treatment And Welfare Office Work (AREA)
Abstract
Description
【0001】[0001]
【発明の属する技術分野】本発明は、ほぼ同一断面を連
続的に計測する際に、アーチファクトの軽減された断層
像を高速で得るのに好適なX線CT装置及び画像再構成
方法に関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an X-ray CT apparatus and an image reconstruction method suitable for obtaining a tomographic image with reduced artifacts at a high speed when continuously measuring substantially the same cross section.
【0002】[0002]
【従来の技術】X線CT装置は広く利用され、様々な利
用方法がユーザーによってなされている。最近では、病
巣の組織検査や治療を経皮的に実施する場合にCTを穿
刺のガイドとして用いることが行われている。CTガイ
ド下で実施することによって手術時間も少なくなり、精
度も上がるとして有効視されている。2. Description of the Related Art X-ray CT apparatuses are widely used, and various uses are made by users. Recently, CT has been used as a puncture guide when performing histological examination or treatment of a lesion percutaneously. It is regarded as effective because the operation time is reduced and the accuracy is improved by performing the procedure under the CT guide.
【0003】このようなCT画像の連続的な観察には、
高速な画像再構成が必要となる。現在、いくつかの提案
がなされている。特公平1−23136号では第1の画
像を再構成するのに用いた投影データと第2の画像を得
るのに必要な投影データとの差分を求め、該差分データ
を再構成した差分画像を第1の画像に加えることで時間
差を有する第2の画像を得るものである。[0003] For continuous observation of such CT images,
High-speed image reconstruction is required. Currently, several proposals have been made. In Japanese Patent Publication No. 23136/1990, the difference between the projection data used to reconstruct the first image and the projection data necessary to obtain the second image is determined, and the difference image obtained by reconstructing the difference data is calculated. A second image having a time difference is obtained by adding the second image to the first image.
【0004】特開平8−24252号では1枚目の画像
を再構成した後、1枚の画像を再構成するのに満たない
投影データ(1/Nビュー:Nは分割数)で再構成した
新たな画像とから新たな断層像を得るとしている。特開
平8−24253号も同様な内容を開示する。In Japanese Patent Application Laid-Open No. Hei 8-24252, after reconstructing a first image, the image is reconstructed with projection data (1 / N view: N is the number of divisions) that is less than one image. A new tomographic image is to be obtained from the new image. Japanese Patent Application Laid-Open No. Hei 8-24253 discloses a similar content.
【0005】これらは共に部分的な再構成をすることで
高速に画像を更新することを可能にしたものである。通
常のスキャンは360°/スキャンで、360°分のデ
ータから1枚の断層像を再構成する。それに対し、特開
平8−24252号等は例えば30°毎の投影データを
再構成した画像(分割画像)を必要枚数(30°の場合
は12枚)加算することで完成した断層像を得る手法で
ある。[0005] Both of them make it possible to update an image at high speed by performing partial reconstruction. A normal scan is 360 ° / scan, and one tomographic image is reconstructed from data for 360 °. On the other hand, Japanese Patent Application Laid-Open No. H8-24252 discloses a method of obtaining a completed tomographic image by adding a required number of images (divided images) (for example, 12 images in the case of 30 °) obtained by reconstructing projection data at every 30 ° It is.
【0006】即ち、図5(a)は、2つのビュー(図中
には0°と30°との2つのビューを図示している)で
の平行投影データの様子を示す図である。ここで平行投
影データ化は、ファンビームX線を数式的に平行ビーム
化することであり、図5(b)には、0°〜30°、3
0°〜60°、…で示す30°単位に仕分けた投影デー
タP(P1、P2…)のメモリ格納例を示す。That is, FIG. 5A is a diagram showing a state of parallel projection data in two views (two views of 0 ° and 30 ° are shown in the figure). Here, the parallel projection data conversion is to mathematically convert the fan beam X-ray into a parallel beam, and FIG.
0 ° to 60 °, a memory storage of a projection data P (P 1, P 2 ... ) which sorted into 30 ° unit represented by ....
【0007】図5(c)は、分割再構成画像g(g1、
g2…)、正規の再構成画像r(r1、r2…)のメモリ
格納例を示す。ここで、分割再構成画像gとは、30°
幅の平行投影データから得た各再構成画像のことであ
り、分割再構成画像g1とは、0°〜30°の平行投影
データから得た再構成画像、分割再構成画像g2とは3
0°〜60°の平行投影データから得た再構成画像であ
り、他の分割再構成画像も同様である。FIG. 5C shows a divided reconstructed image g (g 1 ,
g 2 ...) and an example of memory storage of the normal reconstructed image r (r 1 , r 2 ...). Here, the divided reconstructed image g is 30 °
And that of the reconstructed images obtained from parallel projection data width divided the reconstructed image g 1, 0 ° ~30 ° parallel projection reconstruction image obtained from the data, a dividing reconstructed image g 2 is 3
This is a reconstructed image obtained from parallel projection data of 0 ° to 60 °, and the same applies to other divided reconstructed images.
【0008】正規の再構成画像rとは、複数の分割再構
成画像gを合成(通常、加算)して得た360°分の再
構成画像を指す。通常のCT画像のことである。30°
幅の平行投影データから得た分割再構成画像gは、デー
タ不足のため画像の体をなしていない。そこで、30°
幅で得た12個の分割再構成画像gを加算して正規の再
構成画像を得ることとした。更に、30°幅の分割再構
成画像の加算は逐次的に行うものとする。逐次的に加算
する利点は、リアルタイムで次々に正規の再構成画像を
得ることができる点である。尚、加算することで再構成
できることは周知の事項である。[0008] The normal reconstructed image r refers to a 360 ° reconstructed image obtained by combining (usually adding) a plurality of divided reconstructed images g. It is a normal CT image. 30 °
The divided reconstructed image g obtained from the parallel projection data having the width does not form an image due to insufficient data. So 30 °
The normal reconstructed image is obtained by adding the twelve divided reconstructed images g obtained in the width. Further, the addition of the divided reconstructed images having a width of 30 ° is performed sequentially. The advantage of sequential addition is that regular reconstructed images can be obtained one after another in real time. It is a well-known matter that reconfiguration can be performed by adding.
【0009】即ち、図5(c)で、第1の正規の再構成
画像r1はg1〜g12を加算して得る。第2の正規の再構
成画像r2は、g2〜g13を加算して得る。第2の正規の
再構成画像r2のための画像g2〜g12はすでに第1の正
規の再構成画像r1を得るときに求めておいたものであ
り、新たに必要なのは画像g13のみである。この画像g
13は次の360°全てのスキャンの完了を待つことな
く、次のスキャンの0°〜360°のうちの0〜30°
で得たものであり、データの取得時間は1スキャン分に
比して1/12の時間であり、且つこの30°分の分割
再構成画像を加算するだけで正規の画像r2を得ること
ができた。That is, in FIG. 5C, the first normal reconstructed image r 1 is obtained by adding g 1 to g 12 . Reconstructed image r 2 of the second normal is obtained by adding the g 2 to g 13. The second image g 2 to g 12 for normalization of the reconstructed image r 2 are those already previously asked when obtaining a reconstructed image r 1 of the first normal, the newly required image g 13 Only. This image g
13 is 0-30 ° of 0 ° -360 ° of the next scan without waiting for completion of all next 360 ° scans.
The data acquisition time is 1/12 of the time for one scan, and a regular image r 2 can be obtained only by adding the 30 ° divided reconstructed images. Was completed.
【0010】このように、2枚目の画像r2は、さら
に、13枚目の分割画像g13が再構成されればg2〜g
12の加算で得た加算分割画像に13枚目を加算すること
で求めることができる。従って、1枚の分割画像g13を
再構成する時間Δt1と画像加算する時間Δt2で新たな
断層像を得ることができる。再構成時間は逆投影ビュー
数に比例するので、1枚の分割画像は通常再構成に比べ
約1/12の短い時間で再構成でき、高速な再構成が可
能となる。Δt1+Δt2がスキャンタイムの1/12で
あれば計測と平行してリアルタイムに断層像を連続的に
得ることができる。3枚目以降のr3、r4でも同様であ
る。As described above, the second image r 2 is further converted into g 2 to g 2 if the thirteenth divided image g 13 is reconstructed.
It can be obtained by adding the thirteenth image to the added divided image obtained by the addition of twelve . Therefore, a new tomographic image can be obtained in the time Δt 1 for reconstructing one divided image g 13 and the time Δt 2 for image addition. Since the reconstruction time is proportional to the number of backprojected views, one divided image can be reconstructed in a shorter time, about 1/12 compared to normal reconstruction, and high-speed reconstruction is possible. If Δt 1 + Δt 2 is 1/12 of the scan time, a tomographic image can be continuously obtained in real time in parallel with the measurement. The same applies to r 3 and r 4 for the third and subsequent sheets.
【0011】[0011]
【発明が解決しようとする課題】これらいずれのアルゴ
リズムも高速化を基本的な考えとしている。リアルタイ
ムCTでは穿刺針や螺旋スキャン時の計測断面のモニタ
リングなど動いているもの撮影対象とする場合が多い。
しかし、これらの公知例はこのような動いているものを
撮影したときに現れるCT装置特有のアーチファクトの
抑制は考慮されていない。All of these algorithms have a basic idea of speeding up. In real-time CT, a moving object such as monitoring of a puncture needle or a measurement cross section at the time of a spiral scan is often taken as an imaging target.
However, these known examples do not consider suppression of an artifact peculiar to a CT apparatus which appears when capturing an image of a moving object.
【0012】また、CT装置のアーチファクトは計測開
始と終了時のデータの不連続性が最も顕著に現れること
は同業者の間では最もよく理解されている点で、180
°分のデータを用いた再構成時には投影方向が互いに逆
の関係にある0°と180°のデータの不連続性が問題
となる。このような不連続性の補正アルゴリズムは従来
から用いられているが演算時間の点で実施が難しく、リ
アルタイム再構成画像では省略されているのが現状であ
った。It is well understood by those skilled in the art that the artifact of the CT apparatus has the most remarkable discontinuity of data at the start and end of the measurement.
At the time of reconstruction using data of °, the discontinuity of the data at 0 ° and 180 ° where the projection directions are opposite to each other becomes a problem. Such a discontinuity correction algorithm has been conventionally used, but it is difficult to implement in terms of calculation time, and at present, it has been omitted in real-time reconstructed images.
【0013】本発明は、連続的にまたは継続的に撮影を
実施した場合において、高速に連続的な画像を再構成す
ることを可能とすると共にアーチファクトの低減を可能
にするX線CT装置と画像再構成方法を提供することを
目的とする。The present invention relates to an X-ray CT apparatus and an X-ray CT apparatus which enable a continuous image to be reconstructed at a high speed and reduce artifacts when imaging is performed continuously or continuously. It is intended to provide a reconstruction method.
【0014】[0014]
【課題を解決するための手段】上記目的を達成するため
に本発明によるX線CT装置は、ファンビームを平行ビ
ームに変換して得られた平行ビーム投影データ群を所定
の投影角度で分割し、各分割区分それぞれに対応した複
数個の分割構成画像を記憶するメモリと、該複数個の分
割構成画像のそれぞれに所定の重み係数を乗算する複数
の乗算器と、該複数の乗算器の出力を合成して正規の再
構成画像を得る合成器とを含んで構成されている。In order to achieve the above object, an X-ray CT apparatus according to the present invention divides a parallel beam projection data group obtained by converting a fan beam into a parallel beam at a predetermined projection angle. A memory for storing a plurality of divided constituent images corresponding to each of the divided sections, a plurality of multipliers for multiplying each of the plurality of divided constituent images by a predetermined weighting coefficient, and outputs of the plurality of multipliers And a synthesizer for obtaining a normal reconstructed image by synthesizing.
【0015】本発明に於いては、連続した分割再構成画
像を一部の分割再構成画像が重なるように逐次的に複数
個選び、選んだ分割再構成画像に所定の重み係数を乗算
するようにしている。その際に、選ばれた複数個の分割
再構成画像については、ビュー角度が互いに対向関係に
ある分割再構成画像の重みの和が、その他の分割再構成
画像の重みと等しくなるように重み係数が設定される。In the present invention, a plurality of continuous divided reconstructed images are sequentially selected so that some of the divided reconstructed images overlap, and the selected divided reconstructed image is multiplied by a predetermined weighting coefficient. I have to. At this time, for the selected plurality of divided reconstructed images, the weighting factor is set such that the sum of the weights of the divided reconstructed images whose view angles are opposite to each other is equal to the weight of the other divided reconstructed images. Is set.
【0016】また、180゜分の分割再構成画像を用
い、リカーシブフィルタを用いて再構成画像を得る場
合、ビュー角度が互いに対向関係にある分割再構成画像
の重みの和が1になるように画像間に施すリカーシブフ
ィルタの重みを決定している。When a reconstructed image is obtained by using a recursive filter by using a 180 ° divided reconstructed image, the sum of the weights of the divided reconstructed images whose view angles are opposed to each other is set to one. The weight of the recursive filter applied between images is determined.
【0017】更に本発明は、所定範囲の投影角の再構成
画像を所定角度づつ順次ずらして連続的な再構成画像を
得るX線CT装置に於いて、ファンビーム計測データを
平行ビームに変換して得られた平行ビーム投影データを
所定の投影角度で分割し、各分割区分にそれぞれが対応
する複数個の分割再構成画像を得るステップと、前記複
数個の分割再構成画像のそれぞれに所定の重み係数を乗
算するステップと、重み係数を乗算された前記複数個の
分割再構成画像を合成して1枚の再構成画像を得るステ
ップと、を含む画像再構成方法を開示する。Further, according to the present invention, in an X-ray CT apparatus in which a reconstructed image having a projection angle in a predetermined range is sequentially shifted by a predetermined angle to obtain a continuous reconstructed image, fan beam measurement data is converted into a parallel beam. Dividing the obtained parallel beam projection data at a predetermined projection angle to obtain a plurality of divided reconstructed images corresponding to the respective divided sections; and An image reconstruction method including a step of multiplying a weight coefficient and a step of combining the plurality of divided reconstructed images multiplied by the weight coefficient to obtain one reconstructed image is disclosed.
【0018】更に本発明は、所定範囲の投影角の再構成
画像を所定角度づつ順次ずらして連続的な再構成画像を
得るX線CT装置に於いて、ファンビーム計測データを
平行ビームに変換して得られた180゜分の平行ビーム
投影データを所定の投影角度で分割し、各分割区分にそ
れぞれが対応した複数個の分割構成画像を得るステップ
と、新たに構成された複数の分割再構成画像に第1の重
みを付与する第1の重み乗算ステップと、前記直前の画
像フレームに第2の重みを付与する第2の重み乗算ステ
ップと、前記第1の重みを付与された直前の画像フレー
ムと前記第2の重みを付与された新たに構成された複数
の分割再構成画像を合成して1つの正規の再構成画像を
構成するステップと、を含む画像再構成方法に於いて、
前記第1、第2の重み乗算ステップは、ビュー角度が互
いに対向関係にある分割再構成画像の重みの和が1にな
るように前記第1、第2の重みを設定するステップと、
を含む画像再構成方法を開示する。更に本発明は、所定
範囲の投影角の再構成画像を所定角度づつずらして連続
的な再構成画像を得るX線CT装置に於いて、ファンビ
ーム計測データを平行ビームに変換して得られた平行ビ
ーム投影データのうち、所定の投影角に相当するm個の
分割再構成画像を得るステップと、前記m個の分割再構
成画像のそれぞれに所定の重みを付与するステップと、
前記重み付けされたm個の分割再構成画像を加算して1
枚の画像を再構成するステップとを含む画像再構成方法
を開示する。Further, according to the present invention, in an X-ray CT apparatus for sequentially reconstructing a reconstructed image having a projection angle within a predetermined range by a predetermined angle to obtain a continuous reconstructed image, the fan beam measurement data is converted into a parallel beam. Dividing the obtained 180 ° parallel beam projection data at a predetermined projection angle to obtain a plurality of divided constituent images respectively corresponding to each divided section, and a plurality of newly formed divided reconstructions A first weight multiplying step for assigning a first weight to an image, a second weight multiplying step for assigning a second weight to the immediately preceding image frame, and an immediately preceding image to which the first weight is assigned Composing a frame and a plurality of newly constructed divided reconstructed images given the second weight to form one regular reconstructed image.
The first and second weight multiplying steps include: setting the first and second weights so that the sum of the weights of the divided reconstructed images whose view angles are opposite to each other becomes 1;
An image reconstruction method is disclosed. Further, the present invention provides an X-ray CT apparatus which obtains a continuous reconstructed image by shifting a reconstructed image having a projection angle in a predetermined range by a predetermined angle, and is obtained by converting fan beam measurement data into a parallel beam. Obtaining, from the parallel beam projection data, m divided reconstructed images corresponding to a predetermined projection angle; and assigning a predetermined weight to each of the m divided reconstructed images;
The weighted m divided reconstructed images are added to obtain 1
Reconstructing a single image.
【0019】[0019]
【発明の実施の形態】図4に本実施の形態のファンビー
ムX線CT装置を示した。全体を統括するホストコンピ
ュータ1と、X線発生系、X線検出系などを搭載し、ス
リップリングによって1秒/1回転の連続スキャンが可
能なスキャナ2、前処理、画像再構成処理や各種解析処
理を担当する画像処理装置3、X線に高電圧を供給する
高電圧発生装置4、患者テーブル5、表示装置6からな
る。ホストコンピュータ1は設定された撮影のプランニ
ングに応じ高電圧発生装置4やテーブル5を制御する。
撮影準備が完了し被検体の計測が開始されると計測され
た投影データは逐次、画像処理装置3に転送され、再構
成が行われ、この再構成画像が表示装置6に表示され
る。FIG. 4 shows a fan beam X-ray CT apparatus according to the present embodiment. A host computer 1 that controls the whole, a scanner 2 equipped with an X-ray generation system, an X-ray detection system, etc., capable of continuous scanning of 1 second / 1 rotation by a slip ring, preprocessing, image reconstruction processing and various analysis It comprises an image processing device 3 in charge of processing, a high voltage generator 4 for supplying a high voltage to X-rays, a patient table 5, and a display device 6. The host computer 1 controls the high-voltage generator 4 and the table 5 according to the set imaging plan.
When the preparation for imaging is completed and measurement of the subject is started, the measured projection data is sequentially transferred to the image processing device 3 and reconstructed, and the reconstructed image is displayed on the display device 6.
【0020】図1は、画像処理装置3及び表示装置6の
内部構成例図を示す。画像処理装置3は、再構成演算器
11、重み付け画像加算器12より成り、表示装置6
は、入力インターフェース26、LUT(ルックアップ
テーブル)回路27、表示用画像メモリ28、CRT2
9より成る。再構成演算器11は、投影データメモリ2
0、前処理演算器21、ファンビーム一平行ビーム変換
器22、フィルタ補正演算器23、逆投影演算器24よ
り成る。重み付け画像加算器12は、7個の画像メモリ
10(#1〜#7)、7個の重み係数乗算器13(W1
〜W7が重み係数)、画像合成器25より成る。FIG. 1 shows an example of the internal configuration of the image processing device 3 and the display device 6. The image processing device 3 includes a reconstruction arithmetic unit 11 and a weighted image adder 12, and the display device 6
Are input interface 26, LUT (lookup table) circuit 27, display image memory 28, CRT2
9 The reconstruction arithmetic unit 11 includes a projection data memory 2
0, a pre-processing calculator 21, a fan beam-to-parallel beam converter 22, a filter correction calculator 23, and a back projection calculator 24. The weighted image adder 12 includes seven image memories 10 (# 1 to # 7) and seven weight coefficient multipliers 13 (W 1
WW 7 are weighting factors), and an image synthesizer 25.
【0021】図1の画像処理装置3は、1画像当たり例
えば1秒未満で再構成できるようにしてある。これは、
分割再構成画像を30°幅で次々に逐次的に得ること、
更に、1つの再構成画像は、前回のものに新たに得た最
新の30°幅の分割再構成画像を加算することで達成で
きることによる。例えば、30゜幅での角度更新では1
2枚/秒の再構成画像を得ることができる。The image processing apparatus 3 shown in FIG. 1 is capable of reconstructing an image in less than one second, for example. this is,
Sequentially obtaining divided reconstructed images one after another with a width of 30 °;
Furthermore, one reconstructed image can be achieved by adding the latest 30 ° -width divided reconstructed image newly obtained to the previous one. For example, if the angle is updated at a width of 30 °, 1
Two reconstructed images per second can be obtained.
【0022】図1の再構成演算器11は、演算器21に
よる前処理、変換器22による平行ビーム化、演算器2
3によるフィルタ補正処理、演算器24による逆投影演
算を行うことで再構成像を得る。この再構成演算は、3
60°分の一括再構成ではなく、部分的な角度幅(例え
ば30°幅)での平行ビームデータから得た分割再構成
画像を得る演算である。重み付け画像加算器12は、次
々に演算器11で得られる分割再構成画像(g1、g2、
…)を、順次、画像メモリ10の各メモリ#1〜#7に
割り当てる。例えば、分割再構成画像g1は#1へ、g2
は#2へ、…、g7は#7へ割り当てて格納する。それ
以降の分割再構成画像g8、g9、…は、g1に代わって
g8が#1、g2に代わってg9が#2、…の如く割り当
てて格納を行う。乗算器13は、各メモリ#1〜#7の
画像に対して、重み係数W1〜W7を対応付けて乗算して
各分割画像の重み付けを行う。合成器25では、総加算
を行って、1つの再構成画像を得る。The reconstruction operation unit 11 shown in FIG. 1 includes a pre-processing unit 21, a parallel beam conversion by a converter 22, and an operation unit 2.
The reconstructed image is obtained by performing the filter correction processing by 3 and the back projection calculation by the arithmetic unit 24. This reconstruction operation is 3
This is an operation for obtaining a divided reconstructed image obtained from parallel beam data with a partial angular width (for example, a 30 ° width), instead of a batch reconstruction for 60 °. The weighted image adder 12 sequentially generates the divided reconstructed images (g 1 , g 2 ,
..) Are sequentially assigned to the memories # 1 to # 7 of the image memory 10. For example, dividing the reconstructed image g 1 is the # 1, g 2
Is to # 2, ..., g 7 stores assigned to # 7. After that, the divided reconstructed images g 8 , g 9 ,... Are stored by assigning g 8 as # 1 instead of g 1 and g 9 as # 2 instead of g 2 . The multiplier 13 for each memory # 1 to # 7 of the image, performs weighting of each divided image by multiplying in association with weighting coefficients W 1 to W-7. The synthesizer 25 performs total addition to obtain one reconstructed image.
【0023】ここで、メモリ10の数を7個としたが、
基本的には分割再構成を得るための幅が30°分であれ
ば最小限6個でよい。最小限6個とする理由は、ハーフ
スキャン(180°)によって1つの再構成画像を得る
ことができるためである。(180°=30°分×
6)。ハーフスキャンとすることで、360°分に比し
て、メモリ容量が少なくてよく、また合成器25での加
算数も少なくてよいために再構成の高速化を達成できる
こと、ハーフスキャンが半分の回転角度(180°分)
であることから360°分に比して半分の時間での計測
で再構成が可能であること、の利点を持つ。この場合、
1回転1秒(1秒スキャンとも呼ばれる)のスキャナに
あっては、1秒間で12枚の、重み付け加算した再構成
像、の出力が可能となる(即ち、フレームレートが12
枚/秒となる。フレームレートとは単位時間に得ること
ができる再構成画像の枚数のこと)。尚、180°分と
は、平行ビームで180°分とのことであり、ファンビ
ームでの180°分では再構成するにデータが不足する
故、ファンビーム上では180°+αとする。また、平
行ビームは等間隔と不等間隔とのいずれでもよい。Here, the number of the memories 10 is set to seven.
Basically, if the width for obtaining the divided reconstruction is 30 °, the minimum is six. The reason why the number is set to a minimum of six is that one reconstructed image can be obtained by half scan (180 °). (180 ° = 30 ° minute ×
6). By using half scan, the memory capacity can be reduced and the number of additions in the combiner 25 can be reduced as compared with 360 °, so that high-speed reconstruction can be achieved. Rotation angle (180 °)
Therefore, there is an advantage that reconstruction can be performed by measurement in half the time as compared with 360 °. in this case,
In the case of a scanner with one rotation per second (also referred to as a one-second scan), it is possible to output twelve weighted and added reconstructed images per second (that is, a frame rate of 12).
Sheets / sec. The frame rate is the number of reconstructed images that can be obtained in a unit time). Note that 180 ° means 180 ° for a parallel beam, and 180 ° for a fan beam is insufficient for reconstructing data, so it is 180 ° + α on the fan beam. Further, the parallel beams may be either at regular intervals or irregular intervals.
【0024】本発明において基本的なメモリ10(#1
〜#7)の使い方、及び重みW1〜W 7の設定例について
以下で述べる。 (1)、W1=W2=…=W6=1とし、W7=0とした
例。この例は図2(a)の事例である。図2(a)中、
平行ビーム投影データPは、ファン−平行ビーム変換器
の出力であり、各P1、P2、…から分割再構成画像
g1、g2、…を得る。gn〜gn+5の重み係数Wは、正規
の再構成画像rnを得るために用いられる分割再構成画
像gn〜gn+5に対する重み係数である。 (a)、第1再構成画像r1を得るための手順は以下の
通りである。メモリ10の#1に0°〜30°の分割再
構成画像g1を格納する。次に、30°〜60°の分割
再構成画像g2をメモリ10の#2に格納する。同様
に、画像g3を#3、画像g4を#4、画像g5を#5、
画像g6を#6に格納する。#7には何も格納しない。
#1〜#6への格納が終わった後で、#1〜#6の画像
を乗算器13を介して取り出して合成器25で加算し
て、0〜180°のハーフスキャンでの正規の再構成画
像r1を得る。 (b)、次に180°〜210°の分割再構成画像7を
#1に格納する。この時、(a)で格納した画像g1は
事前にクリアしておくことで、画像g1に代わって画像
g7が格納される。次に、再び#1〜#6の画像g7、g
2〜g6を読み出し乗算器13を得て合成器25で加算し
て次の正規の再構成画像r2を得る。 (c)、次に210°〜240°の分割再構成画像g8
を、#2に格納する。当然に、#2の先に格納した画像
2はクリアにしておく。同様に、合成器25で画像
g7、g8、g3〜g6の加算を行い、次の正規の再構成画
像r3を得る。 (d)、以下、同様にして、次々に1つの再構成画像r
4、r5、…を得る。以上の手順で、第2の再構成画像r
2は、r1を得てから180°〜210°の30°分の平
行データを得た直後であって画像g7を得るまでの時間
と加算時間とで定まる、短い時間の間に算出できる。g
3、g4、…も同様である。In the present invention, the basic memory 10 (# 1
To # 7) and weight W1~ W 7Setting example
It is described below. (1), W1= WTwo= ... = W6= 1 and W7= 0
Example. This example is the case of FIG. In FIG. 2A,
The parallel beam projection data P is a fan-parallel beam converter.
Output of each P1, PTwoReconstructed image from…
g1, GTwo, ... gn~ Gn + 5Weight coefficient W is
Reconstructed image rnSegmentation image used to obtain
Statue gn~ Gn + 5Is a weighting factor for. (A), first reconstructed image r1The steps to get
It is on the street. In the memory 10 # 1, a 0 ° to 30 ° split
Composition image g1Is stored. Next, a division of 30 ° to 60 °
Reconstructed image gTwoIs stored in # 2 of the memory 10. As well
And image gThreeTo # 3, image gFourTo # 4, image gFiveTo # 5,
Image g6Is stored in # 6. Nothing is stored in # 7.
After storing in # 1 to # 6, images of # 1 to # 6
Is taken out through the multiplier 13 and added by the combiner 25.
And a regular reconstructed image in a half scan of 0 to 180 °
Image r1Get. (B) Then, the divided reconstructed image 7 of 180 ° to 210 ° is
Stored in # 1. At this time, the image g stored in (a)1Is
By clearing in advance, the image g1Picture on behalf of
g7Is stored. Next, the images g of # 1 to # 6 again7, G
Two~ G6Is read out, the multiplier 13 is obtained, and added by the synthesizer 25.
The next regular reconstructed image rTwoGet. (C) Next, a divided reconstructed image g of 210 ° to 240 °8
Is stored in # 2. Naturally, the image stored before # 2
Keep 2 clear. Similarly, the image is synthesized by the synthesizer 25.
g7, G8, GThree~ G6Of the next regular reconstructed image
Image rThreeGet. (D) Hereinafter, similarly, one reconstructed image r one after another
Four, RFive, ... By the above procedure, the second reconstructed image r
TwoIs r1From 180 ° to 210 ° for 30 ° flat
Immediately after obtaining the row data, image g7Time to get
It can be calculated in a short time determined by the time and the addition time. g
Three, GFour, ... are the same.
【0025】(2)、ビュー角度が対向するデータの平
均化をはかった例。この例を図3(a)、図9(a)に
示す。図9(a)は図3(a)の分割再構成画像に対す
る重みの変化を示す。 (a)、この例は、30°分の分割再構成画像を得るご
とに重みW1〜W7を変更する例である。変更の仕方は、
互いに対向する関係にある画像に対して、重みの和が
“1”になるように変更する。この重み付けは、対向関
係にあるデータ相互を加算して加重平均をとり、端部ア
ーチファクトの減少をはかるために行う。 (b)、先ず0°〜180°の第1の再構成画像r
1(これは図2(a)(b)のr1とは内容的に異なる
が、同一記号を使う。以下同じ)を得るには、W1=W7
=1/2、W2〜W6は1とする。この場合は、g1、g7
が対向関係にある。この設定のもとで、0°〜30°の
分割再構成画像g1を#1に格納する。以下、30°〜
60°、60°〜90°、…、180°〜210°の各
分割再構成画像g2、g3…、g7を#2、#3、…、#
7に格納する。そして、これらの画像を読み出して乗算
器13で重み付け乗算を行い、合成器25で加算し、正
規の再構成画像r1を得る。この画像r1で、画像g1と
画像g7とはビュー角度が互いに対向関係にあるデータ
であり、この対向データを加算し、1/2化したことで
加重平均値となる。かくして端部の不連続性が軽減さ
れ、端部不連続によるアーチファクトを軽減できた。 (c)、次に第2の再構成画像r2を得るには、W1=W
2=1/2、W3〜W7は1とする。この場合は、g2とg
8が対向関係にある。この設定のもとで210°〜24
0°の分割再構成画像g8を#1に格納する。この際、
#1に格納してあった画像g1はクリアする。そして、
#1〜#7の画像g8、g2、g3、…g7について、乗算
器13を介して合成器25で加算する。この際、画像g
8とg2とが対向関係になり、加算平均されて、端部の不
連続性を軽減した画像r2を得る。 (d)、画像r3を得るには、W2=W3=1/2、W1、
W4〜W7は1にする。この場合は、g3とg9が対向関係
にある。画像r4を得るにはW3=W4=1/2、W1、W
2、W5〜W7は1とする。以下同様に、逐次的に再構成
画像を得る毎にその切替の短い時間帯の中で重みをその
都度変更する。(2) An example in which data having opposite view angles is averaged. This example is shown in FIGS. 3A and 9A. FIG. 9A shows a change in weight for the divided reconstructed image of FIG. 3A. (A), This example is an example in which the weights W 1 to W 7 are changed every time a divided reconstructed image for 30 ° is obtained. How to change
Changes are made so that the sum of the weights becomes "1" for the images that are in a relationship facing each other. This weighting is performed in order to reduce edge artifacts by taking a weighted average by adding data in the opposite relationship to each other. (B) First, a first reconstructed image r of 0 ° to 180 °
In order to obtain 1 (this is different in content from r 1 in FIGS. 2A and 2B, the same symbol is used, the same applies hereinafter), W 1 = W 7
= 1 /, and W 2 to W 6 are 1. In this case, g 1 , g 7
Are opposed to each other. Under this setting, 0 is stored in ° to 30 ° divided reconstructed image g 1 to # 1. Below, 30 ° ~
The divided reconstructed images g 2 , g 3, ..., G 7 of 60 °, 60 ° to 90 °,..., 180 ° to 210 ° are represented by # 2, # 3,.
7 is stored. Then, a weighting multiplier in the multiplier 13 reads out these images are added in combiner 25 to obtain a reconstructed image r 1 regular. In this image r 1, the image g 1 and the image g 7 is data in the opposite relationship to each other view angle, adds the counter data, a weighted average value by ized 1/2. Thus, the discontinuity at the end was reduced, and the artifact due to the discontinuity at the end could be reduced. (C) Then, to obtain the second reconstructed image r 2 , W 1 = W
2 = 1/2 and W 3 to W 7 are set to 1. In this case, g 2 and g
8 are facing each other. Under this setting 210 ° ~ 24
The 0 ° divided reconstructed images g 8 stored in # 1. On this occasion,
# Image g 1, which had been stored in the 1 to clear. And
# Image g 8 of 1~ # 7, g 2, g 3, ... the g 7, adding at the combiner 25 via the multiplier 13. At this time, the image g
And 8 and g 2 become opposite relationship, are averaged to obtain an image r 2 with reduced discontinuities in the end. (D) To obtain an image r 3 , W 2 = W 3 = 1 /, W 1 ,
W 4 to W 7 are set to 1. In this case, g 3 and g 9 are in opposed relationship. To obtain an image r 4 W 3 = W 4 = 1/2, W 1, W
2, W 5 ~W 7 is set to 1. Similarly, every time a reconstructed image is successively obtained, the weight is changed each time in a short time period of the switching.
【0026】(3)、図2(b)は2分割画像を得る毎
に1つの再構成画像を得る例。W1〜W6は1、W7=0
にする。 (a)、第1再構成画像r1を得るために、0°〜30
°の画像g1を#1、30°〜60°の画像g2を#2、
…、150°〜180°の画像g6を#6へ格納する。
そして、読み出して加算することで、画像r1を得る。
これまでの手順は、図2(a)の例と同じである。 (b)、第2再構成画像r2を得るために、180°〜
210°の画像g7を#1へ格納する。図2(a)では
この後で再構成のための加算を行ったが、本例では、こ
の時点で加算を行わずに、210°〜240°の画像g
8を#2へ格納する。この新しい画像g7、g8を得た時
点で、画像g7、g8、g3〜g6を加算して再構成画像r
2を得る。 (c)、以下、同様に、連続する2つの分割画像を得る
ごとに再構成のための加算を行い、再構成画像r3、
r4、…を得る。(3) FIG. 2B shows an example in which one reconstructed image is obtained every time a two-divided image is obtained. W 1 to W 6 are 1, W 7 = 0
To (A) In order to obtain the first reconstructed image r 1 , 0 ° to 30 °
An image g 1 of ° # 1,30 ° to 60 ° image g 2 and # 2,
..., and stores the image g 6 of 0.99 ° to 180 ° to # 6.
Then, by reading and adding, an image r 1 is obtained.
The procedure so far is the same as the example of FIG. (B) In order to obtain the second reconstructed image r 2 , 180 °
An image g 7 of 210 ° and stored in the # 1. In FIG. 2A, the addition for the reconstruction is performed after this, but in this example, the addition is not performed at this time, and the image g of 210 ° to 240 ° is not added.
8 is stored in # 2. When the new images g 7 and g 8 are obtained, the images g 7 , g 8 and g 3 to g 6 are added and the reconstructed image r
Get two . (C) Hereinafter, similarly, each time two continuous divided images are obtained, addition for reconstruction is performed, and a reconstructed image r 3 ,
r 4 , ...
【0027】(4)、図3(b)、図9(b)は端部の
不連続性の軽減を一層はかった例を示す。図9(b)は
図3(b)の分割再構成画像に対する重みの変化を示
す。この例では、メモリ10は8個(#1〜#8)、及
び重み係数W1〜W8の乗算器13も8個用意する。 (a)、先ず第1再構成画像r1を得る際には、W1=
1/4、W7=3/4、W2=1/2、W8=1/2、W3
〜W6は1に設定する。この場合はg1とg7、g2とg8
がそれぞれ対向関係にある。そして0°〜30°の画像
g1を#1、30°〜60°の画像g2を#2、…、18
0°〜210°の画像g7を#7へ、210°〜240
°の画像g8を#8へ、格納する。そして、乗算器13
を介して読み出し加算することで画像r1を得る。 (b)、第2再構成画像r2を得る際には、W1=W3=
1/2、W2=1/4、W8=3/4、W4〜W7は1に設
定する。この場合は、g2とg8、g3とg9がそれぞれ対
向関係にある。そして、240°〜270°の画像g9
を#1へ格納し、加算し、画像r2を得る。以下同様
に、r3、r4、…を得る。互いに対向関係にある分割画
像データに対する重みの例としては、和が1となる条件
を満たしながら、様々な配分が考えられるがリアルタイ
ム性が重視される場合は最新のデータに対しての重みを
より大きい値にする方が望ましい。(4), FIG. 3 (b) and FIG. 9 (b) show examples in which the discontinuity at the end is further reduced. FIG. 9B shows a change in weight for the divided reconstructed image of FIG. 3B. In this example, the memory 10 is 8 (# 1 to # 8), and a multiplier 13 of the weight coefficients W 1 to W-8 is also eight prepared. (A) First, when the first reconstructed image r1 is obtained, W 1 =
4, W 7 = 3, W 2 = 1 /, W 8 = 1 /, W 3
WW 6 is set to 1. In this case, g 1 and g 7 , g 2 and g 8
Are opposed to each other. The 0 ° to 30 images g 1 of ° # 1,30 ° to 60 ° image g 2 and # 2, ..., 18
0 ° to 210 ° image g 7 to # 7, 210 ° to 240
° pictures g 8 to # 8, are stored. And the multiplier 13
The image r 1 is obtained by reading out and adding via the. (B) When obtaining the second reconstructed image r2, W 1 = W 3 =
、, W 2 1 /, W 8 = 3, and W 4 to W 7 are set to 1. In this case, g 2 and g 8, g 3 and g 9 are in opposed relationship respectively. Then, the image g 9 of 240 ° to 270 °
Was stored in the # 1 adds, to obtain an image r 2. Similarly, r 3 , r 4 ,... Are obtained. As an example of the weight for the divided image data that is opposed to each other, various distributions can be considered while satisfying the condition that the sum is 1, but when real-time performance is important, the weight for the latest data is increased. It is desirable to set a large value.
【0028】(5)、30°区分ではなく60°区分の
例。図3(c)に示す例である。これは、60°区分と
することが特徴である。ビュー角度が対向関係にあるデ
ータの加算を行うために、対向位置の平行投影データ
(例えば、0°〜60°に対して180°〜240°)
をメモリ10に格納するようにしてある。但し、メモリ
10は、#1〜#4の4個でよい。(5) An example of a 60 ° section instead of a 30 ° section. This is an example shown in FIG. This is characterized by a 60 ° section. In order to add data in which the view angles are opposed, parallel projection data at opposed positions (for example, 180 ° to 240 ° for 0 ° to 60 °)
Is stored in the memory 10. However, four memories # 1 to # 4 may be used.
【0029】以上の例の中で、(2)、(4)、(5)
は各フレームに重み係数を設定する必要があった。こう
したやり方をとらずに、同様の目的を達成するやり方も
ある。リカーシブフィルタ化する例である。図6がその
構成図である。図1と異なる点は、メモリ10を6個と
し(#1〜#6)、7個目のメモリ#7を直前フレーム
メモリ33としたこと、乗算器30、31、合成器32
を新たに設けたこと、メモリ10の出力側の乗算器13
をなくしたこと、またメモリ10等を表示装置6側に設
けたこと等である。乗算器30は、単純加算平均のため
の手段であり、この単純加算平均結果と、直前の再構成
画像riとの加算を合成器32で行うことにしている。
合成器32の出力は、表示用画像メモリ28に送られる
と共に、乗算器31データ重み付け(W2)をした上で
メモリ33へ、次の再構成画像ri+1の演算のために格
納する。重みW1とW2についてはW1+W2=1の関係が
ある。W2の値を大きくすると前の画像との平均化は高
まるが残像が多くなるため、画像を見ながら、W1とW2
の値を調整する。In the above examples, (2), (4), (5)
Needed to set a weighting factor for each frame. There are ways to achieve similar goals without doing so. This is an example of recursive filtering. FIG. 6 is a configuration diagram thereof. 1 in that the number of memories 10 is six (# 1 to # 6), the seventh memory # 7 is a previous frame memory 33, multipliers 30 and 31, and a combiner 32.
And the multiplier 13 on the output side of the memory 10
And that the memory 10 and the like are provided on the display device 6 side. The multiplier 30 is a means for simple averaging, this a simple averaging results, and to perform addition of the reconstructed image r i immediately before at the synthesizer 32.
The output of the synthesizer 32 is sent to the display image memory 28 and, after the multiplier 31 data weighting (W 2 ), is stored in the memory 33 for calculation of the next reconstructed image ri + 1. . The weights W 1 and W 2 have a relation of W 1 + W 2 = 1. The average of the previous image A larger value of W 2 is increased becomes much residual image, while looking at the image, W 1 and W 2
Adjust the value of.
【0030】この構成によれば、単純平均でリカーシブ
な結合であるため、重みはW1とW2との2つだけですむ
との利点を持つ。図7(a)〜(c)で図6の構成の動
作を説明する。図7(a)〜(c)で上段はメモリ33
の内容であり、直前の画像フレームを示す。下段は乗算
器30の出力であり、新たに構成された複数の分割再構
成画像から成る新たな画像フレームを表す。先ず、重み
W1=1、W2=1/2に初期設定する。この初期設定時
にはメモリ33はクリアされている。そして計測シーケ
ンスに入る。分割画像メモリ10には、#1に0°〜3
0°、#2に30°〜60°、…、#6に150°〜1
80°の各分割再構成画像g1〜g6を次々に計測のシー
ケンスに送って算出した順序に従って格納する。これら
の画像g1〜g6を読み出し合成器25で加算し、乗算器
30を通す。W1=1である故に、画像g1〜g6の加算
画像がそのまま合成器32に入力する。この段階ではメ
モリ33の内容はクリアされているため、合成器32か
らは、上記乗算器30の出力がそのまま出力となる。こ
の出力は正規の再構成画像r1であり、メモリ28側に
送られると共に、乗算器31を介してメモリ33に格納
される。W2=1/2である故に、再構成画像r1を2分
の1化し、画像メモリ33に格納される。1枚目表示の
時も、高速化等の理由により、W1=W2=1/2と固定
重みにしても良い。According to this configuration, there is an advantage that only two weights, W 1 and W 2 , are required because the recursive connection is a simple average. The operation of the configuration of FIG. 6 will be described with reference to FIGS. 7A to 7C, the upper part is the memory 33.
And indicates the immediately preceding image frame. The lower stage is the output of the multiplier 30 and represents a new image frame composed of a plurality of newly formed divided reconstructed images. First, weights W 1 = 1 and W 2 = 1 / are initially set. At the time of this initial setting, the memory 33 is cleared. Then, the measurement sequence starts. In the divided image memory 10, 0 ° to 3 are stored in # 1.
0 °, 30 ° to 60 ° for # 2, ..., 150 ° to 1 for # 6
The 80 ° divided reconstructed images g 1 to g 6 are sequentially sent to the measurement sequence and stored in the calculated order. These images g 1 to g 6 are read out, added by the synthesizer 25, and passed through the multiplier 30. Since W 1 = 1, the added image of the images g 1 to g 6 is input to the synthesizer 32 as it is. At this stage, since the contents of the memory 33 have been cleared, the output of the multiplier 30 is directly output from the synthesizer 32. This output is a regular reconstructed image r 1, which is sent to the memory 28 and stored in the memory 33 via the multiplier 31. Since W 2 = 1 /, the reconstructed image r 1 is halved and stored in the image memory 33. At the time of displaying the first image, a fixed weight of W 1 = W 2 = 1 / may be used for reasons such as speeding up.
【0031】次に、第2の再構成画像r2を得る処理に
入るが、この際には、重み係数W1=W2=1/2に設定
する。先ず、180°〜210°の平行ビームデータ及
びこれに基づく分割再構成画像g7を、再構成演算器1
1で算出し、メモリ10の#1に格納する。この際、先
に格納しておいたメモリ10の#1の画像g1はリセッ
トしておくことで、新画像g7のみが#1に格納され
る。そして、#1〜#6の画像g7、g2〜g6を読み出
し合成器25で加算した上で乗算器30で1/2化す
る。合成器32では、この乗算器30の出力とメモリ3
3(#7)の出力(前回の再構成画像r1の2分の1の
値)とを加算し(図7(a)参照)し、第2の再構成画
像r2としてメモリ28に送ると共に、乗算器31で1
/2化した上でメモリ33(#7)に、前回の再構成画
像r1に代わって次の再構成画像r3算出用に格納する。
次の再構成画像r3を含む以降の再構成画像の算出のた
めには、W1=W2=1/2としておいてよく、重みの変
更はない。Next, the process for obtaining the second reconstructed image r 2 is started. At this time, the weight coefficient W 1 = W 2 == is set. First, the parallel beam data and the divided reconstructed image g 7 thereunder of 180 ° to 210 °, reconstruction operation unit 1
1 and stored in # 1 of the memory 10. At this time, the image g 1 of # 1 of the memory 10 which has been stored earlier than that you reset, only the new image g 7 is stored in # 1. Then, 1/2 of a multiplier 30 in terms of the sum in the synthesizer 25 reads out the image g 7, g 2 to g 6 of # 1 to # 6. In the synthesizer 32, the output of the multiplier 30 and the memory 3
3 (# 7) output (a value of 1 half of the previous reconstructed image r 1) and the sum was (see FIG. 7 (a)), and sends the second reconstructed image r 2 as the memory 28 At the same time, 1
/ 2 phased on the memory 33 in (# 7), on behalf of the previous reconstructed image r 1 stores for the next reconstructed image r 3 calculated.
For calculation of a reconstructed image including the next reconstructed image r 3 , W 1 = W 2 = 1 / may be set, and the weight is not changed.
【0032】このような動作によれば、分割再構成画像
g2〜g6は1枚目と2枚目の画像に共通に含まれている
ため、図7(a)に示すように画像g1とg7のみが1/
2の重み付けがなされたことになる。第2の再構成画像
r2を得る時には図7(b)に示すように、g1とg7に
ついては、図7(a)の場合と同様に1/4の重みとな
り、第3の再構成画像r3を得る時には図7(c)に示
すようにg1とg7については 1/8の重みとなり、徐
々に画像としては減衰してゆく。画像g2とg8、g3と
g9でも同様に逐次的に減衰してゆく。更に図7(b)
では対向関係にあるg2とg8が1/2の重みとなり、図
7(c)では同様に対向関係にあるg3とg9が1/2の
重みとなる。従って、180°再構成画像にリカーシブ
フィルタを加えたことで、対向データ間の平均化による
アーチファクト、ノイズ低減を実現できる。According to such an operation, the divided reconstructed images g 2 to g 6 are commonly included in the first and second images, and therefore, as shown in FIG. Only 1 and g 7 are 1 /
This means that weighting of 2 has been performed. When obtaining the second reconstructed image r 2, as shown in FIG. 7B, g 1 and g 7 have a weight of 4 as in the case of FIG. when obtaining a structure image r 3 becomes a weight of 1/8 for g 1 and g 7 as shown in FIG. 7 (c), slide into decay as slowly image. Image g 2 and g 8, g 3 and g 9 even slide into attenuated similarly sequentially. Further, FIG.
In FIG. 7, g 2 and g 8 in the opposite relationship have a weight of 、, and similarly, in FIG. 7C, g 3 and g 9 in the opposite relationship have a weight of 1 /. Therefore, by adding the recursive filter to the 180-degree reconstructed image, it is possible to reduce artifacts and noise by averaging the facing data.
【0033】尚、図6では、メモリ10の#1〜#6の
出力側に重み係数を乗算する乗算器を設けていないが、
#1〜#6のすべての出力側に図1と同じように設けて
もよい。図6のメモリ33の出力側でも同じである。ま
た、図1と図6、並びに上述の変形側にあっても、各重
みの設定は、ノイズやアーチファクトの低減効果によっ
て任意に設定可能である。In FIG. 6, a multiplier for multiplying a weight coefficient is not provided on the output side of # 1 to # 6 of the memory 10.
All the output sides of # 1 to # 6 may be provided in the same manner as in FIG. The same applies to the output side of the memory 33 in FIG. 1 and 6, and also on the deformed side described above, the setting of each weight can be arbitrarily set by the effect of reducing noise and artifacts.
【0034】図8(a)、(b)は、重み設定部の構成
図を示す。重み設定部とは、各乗算器13、31、32
等の重み係数を決定するための回路である。図8(a)
は、1個の重みテーブル50より成る例であり、分割区
分番号をアドレスとし、各番号に対応する重み係数
W1、W2、…をデータとして格納するメモリである。分
割区分番号とは、0°〜30°を番号1、30°〜60
°を番号2、…の如く設定されたものであり、これはス
キャナの管理信号から直接又はエンコーダ等を介して得
ることができる。FIGS. 8A and 8B show the configuration of the weight setting unit. The weight setting unit includes the multipliers 13, 31, 32
It is a circuit for determining a weight coefficient such as. FIG. 8 (a)
Is an example of a single weight table 50, which is a memory that stores the division coefficient numbers as addresses and the weight coefficients W 1 , W 2 ,... Corresponding to the numbers as data. The division section numbers are 0 ° to 30 ° as numbers 1 and 30 ° to 60 °.
Is set as number 2,..., Which can be obtained directly from the control signal of the scanner or via an encoder or the like.
【0035】図8(b)は、重みテーブル51、52、
…53の如く、再構成モードに応じた重みテーブルを持
つ例である。例えば、図2(a)の動作を行うためのメ
モリがテーブル51、図2(b)の動作を行うためのメ
モリがテーブル52、図3(a)の動作を行うためのメ
モリがテーブル52、…の如く設けてある。従って、再
構成モード信号(これは入力設定されるものであり、図
2(a)や図2(b)や図3(a)や(b)や(c)の
いずれかの選択指定信号)でテーブル51〜53の1つ
を選択し、且つ分割区分番号信号(30°毎の区分であ
れば、0°〜30°、30°〜60°、…を番号1、
2、…の如く設定され、60°の区分であれば0°〜6
0°、60°〜120°を番号1、2、…の如く設定さ
れたもの)でアドレスを選択し、対応する重み係数
W1、W2、…を読み出す。勿論、重み係数も図6の例で
は、W1、W2の2つのみとか、図1では、W1〜W6の6
個とか種々である。FIG. 8B shows weight tables 51 and 52,
53 is an example having a weight table according to the reconstruction mode. For example, the memory for performing the operation in FIG. 2A is the table 51, the memory for performing the operation in FIG. 2B is the table 52, the memory for performing the operation in FIG. ... is provided. Therefore, a reconstruction mode signal (this is an input setting signal, and any one of the selection designation signals shown in FIGS. 2A, 2B, 3A, 3B, and 3C) , One of the tables 51 to 53 is selected, and the division section number signal (in the case of division by 30 °, 0 to 30 °, 30 ° to 60 °,.
2, etc., and 0 ° to 6 if the section is 60 °
0 °, 60 ° to 120 ° are set as numbers 1, 2,...), And the corresponding weight coefficients W 1 , W 2 ,. Of course, in the example also the weighting factor 6, Toka only two W 1, W 2, in Figure 1, the W 1 to W-6 6
Individual and various.
【0036】[0036]
【発明の効果】本発明によれば、X線CT装置で連続的
にまたは継続的に撮影を実施した場合において、高速に
連続的な画像を再構成することを可能とし、CTガイド
下の生検などの超低線量撮影時においてもアーチファク
トやノイズの低減を実現した。更に、フレームレート
(枚/秒)も、自在に変更でき、表示目的や診断目的等
に応じたフレームレートの自在な選択設定が可能になっ
た。According to the present invention, it is possible to reconstruct a continuous image at a high speed when the image is continuously or continuously taken by the X-ray CT apparatus. It also reduces artifacts and noise during ultra-low-dose radiography such as medical examinations. Further, the frame rate (sheets / second) can be freely changed, and the frame rate can be freely selected and set in accordance with the purpose of display or diagnosis.
【図1】本発明のX線CT装置の主要部の構成例図であ
る。FIG. 1 is a structural example of a main part of an X-ray CT apparatus according to the present invention.
【図2】本発明のX線CT装置での重み係数の設定例図
である。FIG. 2 is a diagram illustrating a setting example of a weight coefficient in the X-ray CT apparatus according to the present invention.
【図3】本発明のX線CT装置での重み係数の他の設定
例図である。FIG. 3 is a diagram illustrating another example of setting a weight coefficient in the X-ray CT apparatus according to the present invention.
【図4】本発明のX線CT装置の全体構成図である。FIG. 4 is an overall configuration diagram of the X-ray CT apparatus of the present invention.
【図5】従来例の説明図である。FIG. 5 is an explanatory diagram of a conventional example.
【図6】本発明のX線CT装置の主要部の他の構成例図
である。FIG. 6 is another configuration example of a main part of the X-ray CT apparatus of the present invention.
【図7】本発明のX線CT装置の他の重み係数の例を示
す図である。FIG. 7 is a diagram showing an example of another weighting factor of the X-ray CT apparatus of the present invention.
【図8】本発明の重みテーブル例図である。FIG. 8 is an example of a weight table according to the present invention.
【図9】図3の分割再構成画像に対する重み変化を示す
図である。FIG. 9 is a diagram showing a weight change for the divided reconstructed image of FIG. 3;
1 ホストコンピュータ 2 スキャナ 3 画像処理装置 4 高電圧発生装置 5 患者テーブル 6 表示装置 10 画像メモリ 11 再構成演算器 12 重み付け画像処理装置 13 乗算器 DESCRIPTION OF SYMBOLS 1 Host computer 2 Scanner 3 Image processing device 4 High voltage generator 5 Patient table 6 Display device 10 Image memory 11 Reconstruction arithmetic unit 12 Weighted image processing device 13 Multiplier
Claims (19)
度づつ順次ずらして連続的な再構成画像を得るX線CT
装置であって、該X線CT装置は、 ファンビーム計測データを平行ビームに変換して得られ
た平行ビーム投影データ群を所定の投影角度で分割し、
各分割区分にそれぞれ対応した複数個の分割再構成画像
を記憶する記憶装置と、 該複数個の分割再構成画像のそれぞれに所定の重み係数
を乗算する複数の乗算器と、 該複数の乗算器の出力を合成して正規の再構成画像を得
る合成器と、 を含むX線CT装置。1. An X-ray CT for obtaining a continuous reconstructed image by sequentially shifting a reconstructed image having a projection angle in a predetermined range by a predetermined angle.
An X-ray CT apparatus that divides a group of parallel beam projection data obtained by converting fan beam measurement data into parallel beams at a predetermined projection angle,
A storage device for storing a plurality of divided reconstructed images respectively corresponding to the respective divided sections; a plurality of multipliers for multiplying each of the plurality of divided reconstructed images by a predetermined weight coefficient; and the plurality of multipliers An X-ray CT apparatus, comprising: a synthesizer that synthesizes an output of (1) to obtain a normal reconstructed image.
記憶装置は、連続した分割再構成画像を一部の分割再構
成画像が重なるように逐次的に複数個選んで記憶するメ
モリを含むX線CT装置。2. The X-ray CT apparatus according to claim 1, wherein the storage device sequentially selects and stores a plurality of continuous divided reconstructed images such that some divided reconstructed images overlap. X-ray CT apparatus including:
重み係数を乗算する複数の乗算器は、ビュー角度が互い
に対向関係にある分割再構成画像の重みの和が、その他
の分割再構成画像の重みと等しくなるように設定された
重み係数を含むX線CT装置。3. The X-ray CT apparatus according to claim 2, wherein the plurality of multipliers for multiplying the weighting coefficient are configured such that the sum of the weights of the divided reconstructed images whose view angles are opposite to each other is the other divided image. An X-ray CT apparatus including a weight coefficient set to be equal to a weight of a reconstructed image.
記憶装置は、180゜分の投影データに対応した複数個
の分割構成画像をそれぞれ記憶するメモリを含むX線C
T装置。4. The X-ray CT apparatus according to claim 1, wherein the storage device includes a memory for storing a plurality of divided constituent images corresponding to 180 ° projection data.
T device.
以上連続的に回転し、この回転で逐次的に得たファンビ
ーム計測投影データを平行ビーム化するファン−平行ビ
ーム変換器と、 前記ファン−平行ビーム変換器の出力を区分し、各区分
毎に分割再構成画像を構成する演算器と、 連続した前記分割再構成画像を一部の分割再構成画像が
重なるように逐次的に複数個選んで記憶する記憶装置
と、 前記記憶装置の出力側に接続され、前記記憶された複数
個の分割再構成画像に重み係数を乗算する複数の乗算器
と、 前記乗算器の出力を合成して正規の再構成画像を構成す
る合成器と、 を含むX線CT装置。5. A fan-parallel beam converter for continuously rotating the fan beam X-ray source at least two times or more and converting the fan beam measurement projection data sequentially obtained by the rotation into a parallel beam; An arithmetic unit that divides the output of the parallel beam converter and forms a divided reconstructed image for each section; and sequentially selects a plurality of continuous divided reconstructed images so that some of the divided reconstructed images overlap. A plurality of multipliers connected to the output side of the storage device for multiplying the stored plurality of divided reconstructed images by a weighting factor; and An X-ray CT apparatus, comprising: a synthesizer configured to form a reconstructed image.
演算器は、180゜分の投影データに対してN分割した
区分毎に前記分割再構成画像を構成する手段を含み、前
記複数の乗算器はビュー角度が互いに対向関係にある分
割再構成画像の重みの和が、その他の分割再構成画像の
重みと等しくなるように設定された重み係数を含み、前
記合成器はN+α個の前記分割再構成画像を用いて前記
正規の再構成画像を構成する手段を含むX線CT装置。6. The X-ray CT apparatus according to claim 5, wherein the computing unit includes a unit configured to form the divided reconstructed image for each of N divided sections of the projection data of 180 °. The plurality of multipliers include a weight coefficient set so that the sum of the weights of the divided reconstructed images whose view angles are opposite to each other is equal to the weight of the other divided reconstructed images. An X-ray CT apparatus including a unit configured to form the normal reconstructed image using the divided reconstructed image.
変換して得られた180゜分の平行ビーム投影データを
所定の投影角度で分割し、各分割区分にそれぞれ対応し
た複数個の分割再構成画像を記憶する記憶装置と、 直前の画像フレームを記憶するフレームメモリと、 新たに構成された複数の分割再構成画像に重みを付与す
る第1の乗算器と、 前記直前の画像フレームに重みを付与する第2の乗算器
と、 重みを付与された前記直前の画像フレームと前記新たに
構成された複数の分割再構成画像を合成して1つの正規
の再構成画像を構成する合成器とを含むX線CT装置で
あって、前記第1、第2の乗算器は、ビュー角度が互い
に対向関係にある分割再構成画像の重みの和が1になる
ように重みを付与する手段を含むX線CT装置。7. A 180 ° parallel beam projection data obtained by converting fan beam measurement data into a parallel beam is divided at a predetermined projection angle, and a plurality of divided reconstructed images respectively corresponding to each divided section. A frame memory for storing the immediately preceding image frame; a first multiplier for weighting a plurality of newly-constructed divided and reconstructed images; and a weighting unit for weighting the immediately preceding image frame. A second multiplier, and a combiner that combines the weighted image frame immediately before and the plurality of newly formed divided reconstructed images to form one normal reconstructed image. An X-ray CT apparatus, wherein the first and second multipliers include means for assigning weights so that the sum of the weights of the divided reconstructed images whose view angles are opposed to each other becomes one. CT device.
度づつ順次ずらして連続的な再構成画像を得るX線CT
装置であって、該X線CT装置は、 ファンビーム計測データを平行ビームに変換して得られ
た平行ビーム投影データのうち、所定の投影角に相当す
るm個の分割再構成画像を得る手段と、 前記m個の分割再構成画像に重みを付与する重み付け手
段と、 前記重み付けされたm個の分割再構成画像を加算して1
枚の画像を再構成する合成手段と、 前記再構成された画像を表示する手段とをX線CT装置
であって、前記合成手段は前記投影角をα゜づつずらし
て得られた前記分割再構成画像から連続的な再構成画像
を得る手段を含むX線CT装置。8. An X-ray CT for obtaining a continuous reconstructed image by sequentially shifting a reconstructed image having a projection angle in a predetermined range by a predetermined angle.
An X-ray CT apparatus for obtaining m divided reconstructed images corresponding to a predetermined projection angle out of parallel beam projection data obtained by converting fan beam measurement data into parallel beams. Weighting means for assigning weights to the m divided reconstructed images; and adding the weighted m divided reconstructed images to 1
An X-ray CT apparatus comprising: synthesizing means for reconstructing a single image; and means for displaying the reconstructed image, wherein the synthesizing means is configured to shift the projection angle by α ゜. An X-ray CT apparatus including means for obtaining a continuous reconstructed image from a constituent image.
重み付け手段は、所定の分割再構成画像に対しては重み
係数を変えて重みを付与する手段を含むX線CT装置。9. An X-ray CT apparatus according to claim 8, wherein said weighting means includes means for changing a weighting factor and assigning weight to a predetermined divided reconstructed image.
記分割再構成画像を得る手段は、前記所定の投影角とし
て180゜を越える角度を設定する手段を含むX線CT
装置。10. The X-ray CT apparatus according to claim 8, wherein the means for obtaining the divided reconstructed image includes a means for setting an angle exceeding 180 ° as the predetermined projection angle.
apparatus.
角度づつ順次ずらして連続的な再構成画像を得るX線C
T装置に於いて、 ファンビーム計測データを平行ビームに変換して得られ
た平行ビーム投影データを所定の投影角度で分割し、各
分割区分にそれぞれが対応する複数個の分割再構成画像
を得るステップと、 前記複数個の分割再構成画像のそれぞれに所定の重み係
数を乗算するステップと、 重み係数を乗算された前記複数個の分割再構成画像を合
成して1枚の再構成画像を得るステップと、 を含む画像再構成方法。11. An X-ray C that obtains a continuous reconstructed image by sequentially shifting a reconstructed image of a predetermined range of projection angles by a predetermined angle.
In the T device, the parallel beam projection data obtained by converting the fan beam measurement data into parallel beams is divided at a predetermined projection angle, and a plurality of divided reconstructed images corresponding to each divided section are obtained. Multiplying each of the plurality of divided reconstructed images by a predetermined weighting factor; and combining the plurality of divided reconstructed images multiplied by the weighting factor to obtain one reconstructed image. An image reconstruction method, comprising:
て、前記分割再構成画像を得るステップは、連続した分
割再構成画像を一部の分割再構成画像が重なるように逐
次的に複数個選んで重み係数を乗算するステップに渡す
ステップを含む画像再構成方法。12. The image reconstructing method according to claim 11, wherein the step of obtaining the divided reconstructed images is performed such that a plurality of consecutive divided reconstructed images are successively arranged so that some of the divided reconstructed images overlap. An image reconstruction method comprising the step of selecting and passing to a step of multiplying by a weighting factor.
て、前記重み係数を乗算するステップは、ビュー角度が
互いに対向関係にある分割再構成画像の重みの和が、そ
の他の分割再構成画像の重みと等しくなるように設定す
るステップを含む画像再構成方法。13. The image reconstructing method according to claim 12, wherein the step of multiplying the weighting coefficient is performed such that a sum of weights of the divided reconstructed images whose view angles are opposite to each other is equal to that of another divided reconstructed image. Image reconstructing method including the step of setting the weight to be equal to
て、前記複数個の分割再構成画像を得るステップは、1
80゜分の投影データに対応する複数の分割再構成画像
を得るステップを含む画像再構成方法。14. The image reconstructing method according to claim 11, wherein the step of obtaining the plurality of divided reconstructed images comprises:
An image reconstruction method including a step of obtaining a plurality of divided reconstructed images corresponding to 80 ° of projection data.
て、前記複数個の分割再構成画像を合成して1枚の再構
成画像を得るステップは、180゜+α゜分の投影デー
タに対応する複数個の分割再構成画像を合成して1枚の
再構成画像を得るステップを含む画像再構成方法。15. The image reconstructing method according to claim 14, wherein the step of combining the plurality of divided reconstructed images to obtain one reconstructed image corresponds to 180 ° + α ° of projection data. An image reconstruction method including a step of combining a plurality of divided reconstructed images to obtain one reconstructed image.
角度づつ順次ずらして連続的な再構成画像を得るX線C
T装置に於いて、 ファンビーム計測データを平行ビームに変換して得られ
た180゜分の平行ビーム投影データを所定の投影角度
で分割し、各分割区分にそれぞれが対応した複数個の分
割構成画像を得るステップと、 新たに構成された複数の分割再構成画像に第1の重みを
付与する第1の重み乗算ステップと、 前記直前の画像フレームに第2の重みを付与する第2の
重み乗算ステップと、 前記第1の重みを付与された直前の画像フレームと前記
第2の重みを付与された新たに構成された複数の分割再
構成画像を合成して1つの正規の再構成画像を構成する
ステップと、 を含む画像再構成方法に於いて、前記第1、第2の重み
乗算ステップは、ビュー角度が互いに対向関係にある分
割再構成画像の重みの和が1になるように前記第1、第
2の重みを設定するステップと、 を含む画像再構成方法。16. An X-ray C that obtains a continuous reconstructed image by sequentially shifting a reconstructed image having a projection angle in a predetermined range by a predetermined angle.
In the T device, the parallel beam projection data of 180 ° obtained by converting the fan beam measurement data into a parallel beam is divided at a predetermined projection angle, and a plurality of divided configurations corresponding to each divided section. Obtaining an image, a first weight multiplying step of giving a first weight to the plurality of newly formed divided reconstructed images, and a second weight of giving a second weight to the immediately preceding image frame Multiplying the image frame immediately before the first weighted image and the plurality of newly constructed divided reconstructed images to which the second weight is applied to form one normal reconstructed image And a weighting multiplying step, wherein the first and second weight multiplying steps are performed such that the sum of weights of the divided reconstructed images whose view angles are opposed to each other becomes one. 1st, 2nd Image reconstructing method comprising the steps of: setting a weight, a.
角度づつずらして連続的な再構成画像を得るX線CT装
置に於いて、 ファンビーム計測データを平行ビームに変換して得られ
た平行ビーム投影データのうち、所定の投影角に相当す
るm個の分割再構成画像を得るステップと、 前記m個の分割再構成画像のそれぞれに所定の重みを付
与するステップと、 前記重み付けされたm個の分割再構成画像を加算して1
枚の画像を再構成するステップとを含む画像再構成方
法。17. An X-ray CT apparatus which obtains a continuous reconstructed image by shifting a reconstructed image having a projection angle in a predetermined range by a predetermined angle and converting the fan beam measurement data into a parallel beam. Obtaining, from the parallel beam projection data, m divided reconstructed images corresponding to a predetermined projection angle; and applying a predetermined weight to each of the m divided reconstructed images; Addition of m divided reconstructed images gives 1
Reconstructing a number of images.
て、前記分割再構成画像を得るステップは、前記投影角
をα゜づつずらす毎に前記m個の分割再構成画像を得る
ステップを含む画像再構成方法。18. The image reconstruction method according to claim 17, wherein the step of obtaining the divided reconstructed images includes the step of obtaining the m divided reconstructed images each time the projection angle is shifted by α ゜. Image reconstruction method.
て、前記分割再構成画像を得るステップは、180゜を
越える投影角に対応する平行ビーム投影データから前記
m個の分割再構成画像を得るステップを含む画像再構成
方法。19. The image reconstruction method according to claim 18, wherein the step of obtaining the divided reconstructed images comprises the step of reconstructing the m divided reconstructed images from parallel beam projection data corresponding to a projection angle exceeding 180 °. An image reconstruction method comprising the step of obtaining.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP12051198A JPH119589A (en) | 1997-04-30 | 1998-04-30 | X-ray ct tomograph and image recomposing method |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP11230297 | 1997-04-30 | ||
JP9-112302 | 1997-04-30 | ||
JP12051198A JPH119589A (en) | 1997-04-30 | 1998-04-30 | X-ray ct tomograph and image recomposing method |
Publications (1)
Publication Number | Publication Date |
---|---|
JPH119589A true JPH119589A (en) | 1999-01-19 |
Family
ID=26451503
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP12051198A Pending JPH119589A (en) | 1997-04-30 | 1998-04-30 | X-ray ct tomograph and image recomposing method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPH119589A (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0969414A2 (en) * | 1998-07-01 | 2000-01-05 | General Electric Company | Computerized tomographic multi-frame image reconstruction method and apparatus for helical scanning |
JP2000271113A (en) * | 1999-01-20 | 2000-10-03 | Toshiba Corp | Computerized tomograph |
JP2003528683A (en) * | 2000-03-30 | 2003-09-30 | ザ、ボード、オブ、トラスティーズ、オブ、ザ、ユニバシティー、オブ、イリノイ | Fast hierarchical backprojection method for 3D Radon transform |
JP2011036671A (en) * | 2009-08-12 | 2011-02-24 | Toshiba Corp | Method for generating computed tomogram and computed tomography apparatus |
JP2011067659A (en) * | 1999-01-20 | 2011-04-07 | Toshiba Corp | Computer tomographic apparatus |
WO2011122613A1 (en) * | 2010-03-30 | 2011-10-06 | 株式会社 日立メディコ | Reconstruction computing device, reconstruction computing method, and x-ray ct device |
CN107260195A (en) * | 2016-02-29 | 2017-10-20 | 通用电气公司 | The system and method that artifact for Computed tomography is removed |
JP2018000618A (en) * | 2016-07-04 | 2018-01-11 | 株式会社日立製作所 | Image reconstruction device, x-ray ct device, and image reconstruction method |
JP2019516083A (en) * | 2016-03-25 | 2019-06-13 | エルスリー・セキュリティー・アンド・ディテクション・システムズ・インコーポレイテッドL−3 Communications Security and Detection Systems,Inc. | System and method for reconstructing a projection image from a computed tomography volume |
-
1998
- 1998-04-30 JP JP12051198A patent/JPH119589A/en active Pending
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0969414A3 (en) * | 1998-07-01 | 2003-01-29 | General Electric Company | Computerized tomographic multi-frame image reconstruction method and apparatus for helical scanning |
EP0969414A2 (en) * | 1998-07-01 | 2000-01-05 | General Electric Company | Computerized tomographic multi-frame image reconstruction method and apparatus for helical scanning |
JP2014208300A (en) * | 1999-01-20 | 2014-11-06 | 株式会社東芝 | Computed tomography apparatus |
JP2000271113A (en) * | 1999-01-20 | 2000-10-03 | Toshiba Corp | Computerized tomograph |
JP2011067659A (en) * | 1999-01-20 | 2011-04-07 | Toshiba Corp | Computer tomographic apparatus |
JP2016168508A (en) * | 1999-01-20 | 2016-09-23 | 東芝メディカルシステムズ株式会社 | Computed tomography apparatus |
JP2003528683A (en) * | 2000-03-30 | 2003-09-30 | ザ、ボード、オブ、トラスティーズ、オブ、ザ、ユニバシティー、オブ、イリノイ | Fast hierarchical backprojection method for 3D Radon transform |
JP2011036671A (en) * | 2009-08-12 | 2011-02-24 | Toshiba Corp | Method for generating computed tomogram and computed tomography apparatus |
US8687871B2 (en) | 2009-08-12 | 2014-04-01 | Kabushiki Kaisha Toshiba | Image domain based noise reduction for low dose computed tomography fluoroscopy |
JPWO2011122613A1 (en) * | 2010-03-30 | 2013-07-08 | 株式会社日立メディコ | Reconstruction calculation device, reconstruction calculation method, and X-ray CT apparatus |
US9050003B2 (en) | 2010-03-30 | 2015-06-09 | Hitachi Medical Corporation | Reconstruction computing device, reconstruction computing method, and X-ray CT apparatus |
JP5960048B2 (en) * | 2010-03-30 | 2016-08-02 | 株式会社日立製作所 | Reconstruction calculation device, reconstruction calculation method, and X-ray CT apparatus |
WO2011122613A1 (en) * | 2010-03-30 | 2011-10-06 | 株式会社 日立メディコ | Reconstruction computing device, reconstruction computing method, and x-ray ct device |
CN107260195A (en) * | 2016-02-29 | 2017-10-20 | 通用电气公司 | The system and method that artifact for Computed tomography is removed |
CN107260195B (en) * | 2016-02-29 | 2022-10-14 | 通用电气公司 | System and method for artifact removal for computed tomography imaging |
JP2019516083A (en) * | 2016-03-25 | 2019-06-13 | エルスリー・セキュリティー・アンド・ディテクション・システムズ・インコーポレイテッドL−3 Communications Security and Detection Systems,Inc. | System and method for reconstructing a projection image from a computed tomography volume |
JP2018000618A (en) * | 2016-07-04 | 2018-01-11 | 株式会社日立製作所 | Image reconstruction device, x-ray ct device, and image reconstruction method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5907593A (en) | Image reconstruction in a CT fluoroscopy system | |
US6628742B2 (en) | Cardiac helical half scan reconstructions for multiple detector row CT | |
US5559847A (en) | Systems, methods and apparatus for reconstructing images in a CT system implementing a helical scan | |
US20020136350A1 (en) | Phase-driven multisector reconstruction for multislice helical CT imaging | |
US6324247B1 (en) | Partial scan weighting for multislice CT imaging with arbitrary pitch | |
US5606585A (en) | Methods and apparatus for multislice helical image reconstruction in a computer tomography system | |
EP2030170B1 (en) | Dynamic computed tomography imaging | |
JP2000271112A (en) | Method for forming image of object and imaging system | |
US5513236A (en) | Image reconstruction for a CT system implementing a dual fan beam helical scan | |
JPH08308824A (en) | System to create tomographic image of article | |
US6332013B1 (en) | Methods and apparatus for tilted helical reconstruction multislice CT | |
EP1372115B1 (en) | Methods and apparatus for reconstructing an image of an object | |
JPH10243942A (en) | System for forming basic tomographic image of object and succeeding tomographic image and method for reconstituting image of object | |
JPH119589A (en) | X-ray ct tomograph and image recomposing method | |
JP4509255B2 (en) | Perspective image creation method and apparatus | |
JP3917684B2 (en) | Method and apparatus for creating a tomographic image of an object | |
US6061422A (en) | X-ray ct apparatus and an image reconstructing method for obtaining a normal reconstructed image by adding weights to a plurality of partial reconstructed images | |
JP2000023966A (en) | Spiral reconstituting algorithm | |
US6055291A (en) | Systems, methods and apparatus for reconstructing images | |
US7221730B2 (en) | Multi-row detector x-ray CT apparatus and method for creating tomogram | |
JP2001212133A (en) | Synthetic tomographic image generating method | |
JPH07222741A (en) | Computer aideo tomography | |
JP2006068229A (en) | Image reconstruction method of x-ray ct image and x-ray ct apparatus | |
JP2000271113A (en) | Computerized tomograph | |
JPH0779964A (en) | Computed tomograph |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20040427 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070207 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070529 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20071009 |