JPH1033542A - Reconstituting method for optical computed tomography - Google Patents
Reconstituting method for optical computed tomographyInfo
- Publication number
- JPH1033542A JPH1033542A JP8196440A JP19644096A JPH1033542A JP H1033542 A JPH1033542 A JP H1033542A JP 8196440 A JP8196440 A JP 8196440A JP 19644096 A JP19644096 A JP 19644096A JP H1033542 A JPH1033542 A JP H1033542A
- Authority
- JP
- Japan
- Prior art keywords
- light
- subject
- weight
- light source
- detector
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000003287 optical effect Effects 0.000 title claims abstract description 20
- 238000000034 method Methods 0.000 title claims description 21
- 238000002591 computed tomography Methods 0.000 title description 2
- 238000010521 absorption reaction Methods 0.000 claims description 3
- 230000007423 decrease Effects 0.000 abstract description 2
- 230000015654 memory Effects 0.000 abstract description 2
- 238000003325 tomography Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 238000007796 conventional method Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000031700 light absorption Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
Description
【0001】[0001]
【発明の属する技術分野】本発明は、光を用いたコンピ
ュータ断層撮影法(光CT;光ComputedTom
ography)における、被検体を透過した光を受光
して得た受光信号に基づいて、その被検体の、光の光路
を含む断層面内の光の吸収分布を表わす断層画像を得る
再構成方法に関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a computed tomography method using light (optical CT; optical Computed Tom).
a reconstruction method for obtaining a tomographic image representing a light absorption distribution in a tomographic plane including an optical path of light of the subject based on a light reception signal obtained by receiving light transmitted through the subject in the method of FIG. .
【0002】[0002]
【従来の技術】従来、X線を用いたX線CTが、医療等
の分野で実現されており、近年では、光を使った断層撮
影法(光CT)を生体に適用するための研究が盛んに行
なわれている。生体は高散乱体であるため生体に照射さ
れた光は多重散乱を繰り返しながら吸収されて生体の外
部に出てくる。このようにして生体の外部に出射した光
のうち、生体内で散乱を受けずに前方方向に出てくる光
のみを検出することにより生体の吸収プロファイルを得
る方法として、時間幅の極めて狹いパルス光を被検体に
照射しその被検体からの出射光のうちの最も初期の出射
光のみを検出対象とするいわゆる時間ゲイト法や、被検
体の出射光と所定の参照光とを重ね合わせて被検体内を
直進した成分のみの干渉信号を得るヘテロダイン検出法
等が知られている。しかしながら、これらの方法を用い
て生体内部を直進した光のみを検出対象としたのでは極
めて微弱な光を検出する必要があり、実用的な大きさの
被検体の場合、光の減衰が大き過ぎ得て測定は不能とな
る。実用的な大きさの被検体の測定を実施するために
は、入射軸から少し離れた領域の前方散乱情報、すなわ
ち、完全な直進光のみではなく、被検体内で散乱を受け
て前方に出射する光も含めて検出対象とする必要があ
る。2. Description of the Related Art Conventionally, X-ray CT using X-rays has been realized in the field of medical treatment and the like. In recent years, research has been conducted on applying tomography using light (optical CT) to living bodies. It is being actively conducted. Since the living body is a high scatterer, the light applied to the living body is absorbed while repeating multiple scattering, and comes out of the living body. As a method of obtaining an absorption profile of a living body by detecting only light that is emitted in the forward direction without being scattered in the living body among the light emitted outside the living body in this way, the time width is extremely narrow. A so-called time-gate method in which only the earliest outgoing light of the outgoing light from the subject is irradiated with the pulsed light and the outgoing light of the subject and a predetermined reference light are superimposed. A heterodyne detection method or the like that obtains an interference signal of only a component that travels straight in a subject is known. However, if only light that has traveled straight inside the living body is to be detected using these methods, it is necessary to detect extremely weak light.In the case of a practically sized subject, the attenuation of light is too large. This makes measurement impossible. In order to measure a practically sized subject, forward scattering information in a region slightly away from the incident axis, that is, not only perfectly straight light, but also scattered within the subject and emitted forward It is necessary to include the light to be detected.
【0003】しかしながら、この場合、検出対象とする
出射光被検体内で散乱した光が含まれることから、その
出射光が担持する被検体情報の質が劣化することとな
り、光CT法により再構成して断層像を得ても分解能の
劣化した画像しか得られないおそれがある。However, in this case, since the emitted light to be detected includes light scattered in the object, the quality of the object information carried by the emitted light is degraded. There is a possibility that even if a tomographic image is obtained in such a manner, only an image with reduced resolution can be obtained.
【0004】[0004]
【発明が解決しようとする課題】本発明は、上記事情に
鑑み、光CTにおいて、分解能の優れた断層像を得るた
めの再構成方法を提供することを目的とする。SUMMARY OF THE INVENTION In view of the above circumstances, an object of the present invention is to provide a reconstruction method for obtaining a tomographic image having excellent resolution in optical CT.
【0005】[0005]
【課題を解決するための手段】上記目的を達成する本発
明の光CTの再構成方法は、被検体を透過した光を受光
して得た受光信号に基づいて被検体の断層面内の光の吸
収分布を表わす断層像を得る光CTの再構成方法におい
て、断層像の各ピクセルに、被検体を透過する光の光路
に沿って下流側ほど大きなウェイトを与えて、再構成の
ための逆投影演算を実行することをことを特徴とする。According to the present invention, there is provided a method of reconstructing an optical CT, which achieves the above-mentioned object, based on a light-receiving signal obtained by receiving light transmitted through a subject. In the reconstruction method of the optical CT for obtaining the tomographic image representing the absorption distribution of the tomographic image, a greater weight is given to each pixel of the tomographic image on the downstream side along the optical path of the light transmitted through the subject, and the inverse for the reconstruction is performed. Executing a projection operation.
【0006】ここで、上記本発明の光CTの再構成方法
において、断層像にあらわれる被検体の輪郭を抽出し、
その輪郭の内部領域についてウェイトを与えて、再構成
のための逆投影演算を実行することが好ましい。被検体
から前方(検出器側)に出射した前方散乱光は、検出器
に近い側では散乱されずに直進した可能性が高く、被検
体への光の入射側ほど散乱された可能性が高まる。すな
わち、その前方散乱光に含まれる被検体情報のうち、被
検体の検出器に近い側い部分の情報は信頼性が高く、検
出器から離れた部分(入射側に近い部分)の情報ほど情
報の信頼性が低下する。Here, in the optical CT reconstruction method of the present invention, the contour of the subject appearing in the tomographic image is extracted,
It is preferable to execute a back projection operation for reconstruction by giving a weight to an inner region of the contour. Forward scattered light emitted forward (detector side) from the subject is more likely to have traveled straight without being scattered on the side closer to the detector, and is more likely to be scattered toward the incident side of the light on the subject. . In other words, of the object information included in the forward scattered light, the information of the portion of the object closer to the detector has high reliability, and the information of the portion farther from the detector (the portion closer to the incident side) has more information. The reliability of the device decreases.
【0007】そこで、本発明では、逆投影(バックプロ
ジェクション)を行なうにあたり、断層像の各ピクセル
に、被検体を透過する光の光路に沿って下流側ほど大き
なウェイトを与える。これにより、再構成により得られ
た断層像は、一律のウェイトを持ったこれまでの再構成
により得られる断層像と比べ分解能の優れた断層像を得
ることができる。Therefore, in the present invention, when performing back projection (back projection), a greater weight is given to each pixel of the tomographic image toward the downstream side along the optical path of light transmitted through the subject. Thereby, the tomographic image obtained by the reconstruction can obtain a tomographic image having a uniform weight and higher resolution than the tomographic image obtained by the previous reconstruction.
【0008】ウェイトを与えるにあたっては、散乱は被
検体内部で生じるため、被検体の輪郭を求めて被検体内
部領域のみウェイトを与えると、一層正確なウェイトを
与えることができ、一層高画質の断層像を得ることがで
きる。尚、本発明は、上述した時間ゲイト法や光ヘテロ
ダイン法と組み合わせて用いることもできる。時間ゲイ
ト法や光ヘテロダイン法は直進光を抽出することができ
る点で優れているが、時間ゲイト法によるゲイト時間を
長くしたり、すなわち被検体から出射した光のうちある
程度の長時間の出射光を検出対象としたり、光ヘテロゲ
イン法において参照光軸に幅を持たせるなどにより、検
出光量をあげることによる情報の劣化を回復し、高画質
の断層像を得ることができる。When the weight is given, since scattering occurs inside the subject, if the contour of the subject is determined and the weight is given only to the region inside the subject, a more accurate weight can be given, and a higher-quality tomographic image can be given. An image can be obtained. The present invention can be used in combination with the above-described time gate method or optical heterodyne method. The time-gate method and the optical heterodyne method are excellent in that they can extract straight-ahead light.However, the time-gate method can be used to lengthen the gate time, that is, the light emitted from the subject for a certain long time. , Or by giving a width to the reference optical axis in the optical hetero-gain method, it is possible to recover information deterioration due to an increase in the detected light amount and obtain a high-quality tomographic image.
【0009】[0009]
【発明の実施の形態】以下、本発明の実施形態について
説明する。図1は、本発明の一実施形態の説明図であ
る。被検体10を光源21と検出器22との間に配置
し、光源21から出射した光ビーム21aを被検体10
に入射させて、被検体10からの出射光21bの、断層
面内(図1の紙面内)の一次元的な強度分布(プロファ
イル)40を検出する。1つのプロファイルを得た後回
転中心Oを中心として、図示の矢印の方向に回転させ
る。この操作を繰り返して全方向のプロファイルを得
る。Embodiments of the present invention will be described below. FIG. 1 is an explanatory diagram of one embodiment of the present invention. The subject 10 is arranged between the light source 21 and the detector 22, and the light beam 21 a emitted from the light source 21 is
And a one-dimensional intensity distribution (profile) 40 of the emitted light 21b from the subject 10 in the tomographic plane (in the paper plane of FIG. 1) is detected. After obtaining one profile, it is rotated about the rotation center O in the direction of the arrow shown in the figure. This operation is repeated to obtain profiles in all directions.
【0010】逆投影にあたっては、光源21と検出器2
2との間を直線で結んだときの、その直線の、光源21
と検出器22との間を、例えば10等分し、光源21側
ほど小さな、検出器22側ほど大きなウェイトを、再構
成領域30内の、メモリ上の各ピクセルに与えて、検出
器22から光源21に向かう図1に破線で示す矢印の方
向の逆投影(バックプロジェクション)を全データにつ
いて実行する。これを、各角度毎ウェイトを変更しなが
ら順次実行し、断層像を再構成する。In back projection, the light source 21 and the detector 2
When the light source 21 is connected to the light source 21 by a straight line
And the detector 22 are divided into, for example, 10 equal parts, and a smaller weight on the light source 21 side and a larger weight on the detector 22 side are given to each pixel on the memory in the reconstruction area 30, and Back projection in the direction of the arrow indicated by the broken line in FIG. 1 toward the light source 21 is executed for all data. This is sequentially performed while changing the weight for each angle to reconstruct a tomographic image.
【0011】こうすることにより、ウェイトを与えずに
逆投影(バックプロジェクション)を行なった場合と比
べ、高分解能の断層像を得ることができる。図2は、本
発明のもう1つの実施形態の説明図である。ここでは、
例えばウェイトを与えずに、あるいは図1を参照して説
明した、被検体10の断面形状を無視して光源と検出器
との間を等分してウェイトを与える方法により、ウェイ
トを与えて一旦再構成を行ない、断層像を得る。その
後、その断層像から被検体10の輪郭10aを抽出し、
その輪郭10a上の出射点pから光源側に向かう距離x
の関数としてウェイトを定義する。By doing so, a high-resolution tomographic image can be obtained as compared with the case where back projection is performed without giving a weight. FIG. 2 is an explanatory diagram of another embodiment of the present invention. here,
For example, the method of giving weights by giving equal weight between the light source and the detector ignoring the cross-sectional shape of the subject 10 without giving weights or with reference to FIG. Reconstruction is performed to obtain a tomographic image. Thereafter, the contour 10a of the subject 10 is extracted from the tomographic image,
The distance x from the emission point p on the contour 10a to the light source side
Define the weight as a function of
【0012】すなわち、ウェイトwを、 w=f(x) ……(1) として、再度再構築を行なう。こうすることにより、一
層正確なウェイトを与えることができ、一層高画質の断
層像を得ることができる。ただし(1)式のようにウェ
イトを与えると、輪郭10aの内部のピクセルの各角度
のウェイトの総和は、各ピクセルごとに異なることにな
る。そこでこの場合、ピクセルごとのウェイトの総和が
等しくなるように再構成後にピクセルごとにウェイトの
逆数が掛け合わせる。こうすることによりウェイトの総
和の違いが補正される。That is, the weight w is set to w = f (x) (1), and reconstruction is performed again. By doing so, a more accurate weight can be given, and a higher quality tomographic image can be obtained. However, if weights are given as in equation (1), the sum of the weights at each angle of the pixels inside the contour 10a will be different for each pixel. Therefore, in this case, the reciprocal of the weight is multiplied for each pixel after reconstruction so that the sum of the weights for each pixel becomes equal. By doing so, the difference in the sum of the weights is corrected.
【0013】尚、プロファイルに適当なコンボリュウシ
ョン関数を掛けて補正した後、逆投影を行なうこと等は
従来方式と同様である。It is to be noted that back projection after correcting the profile by multiplying the profile by an appropriate convolution function is the same as the conventional method.
【0014】[0014]
【発明の効果】以上説明したように、本発明によれば、
従来と比べ高画質の断層像を再構成することができる。As described above, according to the present invention,
It is possible to reconstruct a high-quality tomographic image as compared with the related art.
【図1】本発明の一実施形態の説明図である。FIG. 1 is an explanatory diagram of one embodiment of the present invention.
【図2】本発明のもう1つの実施形態の説明図である。FIG. 2 is an explanatory diagram of another embodiment of the present invention.
10 被検体 21 光源 21a 光ビーム 21b 出射光 22 検出器 30 再構成領域 40 プロファイル Reference Signs List 10 subject 21 light source 21a light beam 21b emission light 22 detector 30 reconstruction area 40 profile
Claims (2)
信号に基づいて該被検体の断層面内の光の吸収分布を表
わす断層像を得る光CTの再構成方法において、 前記断層像の各ピクセルに、被検体を透過する光の光路
に沿って下流側ほど大きなウェイトを与えて、再構成の
ための逆投影演算を実行することをことを特徴とする光
CTの再構成方法。1. A method of reconstructing an optical CT for obtaining a tomographic image representing an absorption distribution of light in a tomographic plane of a subject based on a light reception signal obtained by receiving light transmitted through the subject. A method of reconstructing an optical CT, wherein each pixel of the image is given a larger weight toward the downstream side along the optical path of light transmitted through the subject, and a back projection operation for reconstruction is performed. .
抽出し、該輪郭の内部領域についてウェイトを与えて、
再構成のための逆投影演算を実行することを特徴とする
請求項1記載の光CTの再構成方法。2. Extracting a contour of a subject appearing in the tomographic image, giving a weight to an inner region of the contour,
2. The method for reconstructing optical CT according to claim 1, wherein a back projection operation for reconstruction is performed.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP8196440A JP2826297B2 (en) | 1996-07-25 | 1996-07-25 | Optical CT reconstruction method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP8196440A JP2826297B2 (en) | 1996-07-25 | 1996-07-25 | Optical CT reconstruction method |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH1033542A true JPH1033542A (en) | 1998-02-10 |
JP2826297B2 JP2826297B2 (en) | 1998-11-18 |
Family
ID=16357866
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP8196440A Expired - Lifetime JP2826297B2 (en) | 1996-07-25 | 1996-07-25 | Optical CT reconstruction method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2826297B2 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7016523B1 (en) | 1999-04-21 | 2006-03-21 | Hiroyuki Ogawa | Method for observing object by projection, method for detecting microorganisms and projection detecting system |
-
1996
- 1996-07-25 JP JP8196440A patent/JP2826297B2/en not_active Expired - Lifetime
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7016523B1 (en) | 1999-04-21 | 2006-03-21 | Hiroyuki Ogawa | Method for observing object by projection, method for detecting microorganisms and projection detecting system |
Also Published As
Publication number | Publication date |
---|---|
JP2826297B2 (en) | 1998-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100457039C (en) | X-ray scatter correction | |
US7424089B2 (en) | System and method for reconstructing image by using straight-line trajectory scan | |
US7072436B2 (en) | Volumetric computed tomography (VCT) | |
US8483471B2 (en) | Method and system for scatter correction in X-ray imaging | |
US10568588B2 (en) | Tiled detector arrangement for differential phase contrast CT | |
US10878602B2 (en) | Image reconstruction device, X-ray CT device, and image reconstruction method | |
US20040202280A1 (en) | Signal profiling for medical imaging systems | |
US20040202279A1 (en) | Method and apparatus for blocking radiographic scatter | |
US20080118024A1 (en) | Method for reconstructing a local high resolution x-ray ct image and apparatus for reconstructing a local high resolution x-ray ct image | |
JP2004509689A (en) | X-ray computed tomography device | |
JP5052281B2 (en) | Method for estimating scattered ray intensity distribution in X-ray CT and X-ray CT apparatus | |
JP5214942B2 (en) | Method for scaling scattered ray intensity distribution in multi-tube X-ray CT and multi-tube X-ray CT apparatus | |
KR20070051758A (en) | X-ray ct apparatus and x-ray ct fluoroscopic apparatus | |
JP2007529738A (en) | Beam hardening correction and attenuation correction of coherent scattering CT | |
EP2002286A2 (en) | Effective dual-energy x-ray attenuation measurement | |
Zhu et al. | Cone beam micro-CT system for small animal imaging and performance evaluation | |
US7564938B2 (en) | Scatter rejection for composite medical imaging systems | |
US8229060B2 (en) | Medical X-ray examination apparatus and method for k-edge imaging | |
Yang et al. | Combined system of fluorescence diffuse optical tomography and microcomputed tomography for small animal imaging | |
JPH0257B2 (en) | ||
US11134907B2 (en) | Signal processing method and imaging system for scatter correction in computed tomography | |
JP4330881B2 (en) | Method for reducing artifacts in target images | |
JP2004357969A (en) | X-ray measuring instrument | |
JP2826297B2 (en) | Optical CT reconstruction method | |
JP2008154669A (en) | X-ray ct apparatus |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 19980818 |