JPS5911150B2 - 透過性放射線に対する被検査体における平面状仮想輪切り部分の吸収の変化に関する表示を構成する方法および装置 - Google Patents

透過性放射線に対する被検査体における平面状仮想輪切り部分の吸収の変化に関する表示を構成する方法および装置

Info

Publication number
JPS5911150B2
JPS5911150B2 JP49047032A JP4703274A JPS5911150B2 JP S5911150 B2 JPS5911150 B2 JP S5911150B2 JP 49047032 A JP49047032 A JP 49047032A JP 4703274 A JP4703274 A JP 4703274A JP S5911150 B2 JPS5911150 B2 JP S5911150B2
Authority
JP
Japan
Prior art keywords
signal
absorption
storage
signals
beams
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired
Application number
JP49047032A
Other languages
English (en)
Other versions
JPS5028385A (ja
Inventor
ア−チボルド ゴ−ドン レメイ クリストフア
エドワ−ド ベスト ジヨン
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
SOON II EMU AI PLC
Original Assignee
SOON II EMU AI PLC
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 SOON II EMU AI PLC filed Critical SOON II EMU AI PLC
Publication of JPS5028385A publication Critical patent/JPS5028385A/ja
Publication of JPS5911150B2 publication Critical patent/JPS5911150B2/ja
Expired legal-status Critical Current

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/40—Arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/4064—Arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
    • A61B6/4071—Pencil beams
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/032—Transmission computed tomography [CT]
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T12/00—Tomographic reconstruction from projections
    • G06T12/20—Inverse problem, i.e. transformations from projection space into object space
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S378/00—X-ray or gamma ray systems or devices
    • Y10S378/901—Computer tomography program or processor

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Veterinary Medicine (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Public Health (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Pulmonology (AREA)
  • General Physics & Mathematics (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Description

【発明の詳細な説明】 本発明はX線またはγ線のような透過性放射線に対する
例えば人体等の被検査体における平面状の仮想輪切り部
分即ち平面状スライス部分の吸収の変化zご関する表示
を構成する放射線検査装置に関するものである。
本発明による検査装置は、陰極線管や他の映像形成装置
上の画像や、そのような画像の写真や、あるいはデジタ
ル計算機によつて発生されうる吸収係数の分布図(これ
について、後に輪部を描くこともできる)のような任意
適当な形態の放射線写真を作成するために用いられうる
ものである。
透過性放射線に対する被検査体における平面状の仮想輪
切り部分の吸収の変化に関する表示を構成するための装
置の一例が米国特許第3778614号に記載されてい
る。この特許のものにおいて−は、外部放射線源から放
射線がペンシル・ビニム状をなして被検査体の一部分を
透過せしめられる。そのビームには走査運動が与えられ
、従つてそのビームは多数の異なる位置を順次採るよう
になさべかつそのビームが被検査体を通過した後に、そ
のような位置のそれぞれにおけるビームの吸収の測定値
を与えるために検知器が用いられる。ビームがこれら種
々の位置を採るように、放射線源と検知器が1つの平面
内で往復運動をなさしめられかつその平面に直交する軸
線のまわりで軌道運動をなさしめられる。従つて、被検
査体中の一平面内に存在しており、この前記種々の位置
は、被検査体中の一平面内に存在しており、この平面内
の吸収係数の分布が検知器によつて与えられるビ一仏吸
収データを処理することによつて得られる。この場合の
処理は、最終的に表示される吸収の分布が一連の近似の
結果として得られるようにして行われる。その目的を実
現するため上述の特許によるX線検゛査装置11丸上述
の走査運動をしている間検出装置によつて発生された走
査運動信号を各角度位置に対応する一出力信号の組に組
合せる制御手段を有する。
当該出力信号が被検査体のスライス部分を走査する−異
なる方向のビームの組の吸収量を表わし1つの組内の出
力信号の数が50以上である。上記X線検査装置はさら
゛に、出力信号の組を対数信号の組に変換する対数変換
器と、スライス部分に重ね合わされた2次元的な観念上
のマトリクスを構成する基本領域に対応した記憶要素を
有する合計出力マトリクス記憶装置と、処理装置を含む
計算手段と、アトレース選択器とを有し、各対数出力信
号から他の対数出力信号に応じて修正゛された修正信号
を形成しかつ対応するビーム通路に属するマトリクスの
基本領域に対応する出力マトリクス記憶装置の記憶要素
に−上記修正信号を分配することにより出力マートリク
ス記憶装置に補正されたレイヤ゛グラムを形成する信号
を蓄積させる回路を有する。前記米国特許に記載されて
いる−X線検査装置はたとえば頭のような人体の一部分
の横断面表面を発生するのには十分であることが認めら
れた。
本発明の1つの目的は、ビーム吸収データを処理するた
めの新規にして優れた技術を利用し7て、X線やγ線の
ような透過性放射線に対する被検査体内における仮想輪
切り部分の吸収の変化に関する表示を構成するための改
良された検査装置を提供することである。「ザ・プロシ
ーデイングス・オブ・ザ・ナシヨナル・アーカデミ一・
オブ・サイエンスイズ」(米国)第68巻、第9号、第
2236ページ(The一PrOceedingsOf
theNatiOnalAcademy′U.S.A.
,VOl.68,E9,p・Age2236)に197
1年9月に発表されたジ一・エヌ・ラマチヤナドラン(
G−N−Ramachanaram)およびエイ・フィ
・ラクシミナラヤナン(A.V.Lakshnllna
rayan一An)による[放射線写真および電子顕微
鏡写真からの三次元再構成:フーリエ変換にかわるコン
ポリユーシヨンの適用゛」(Three−Divens
iOnaIrecOnstructiOnfrOmRa
diOgraphsandElectrOnmicrO
grapHs:ApplicatiOnOfcOnvO
lutiOnsinsteadOfFOurierTr
ansfOrms)という論文には、前記米国特許に記
載されている装置から得られたもののような吸収データ
から横断面表示を再構成するために原理゛的に適用可能
な数学的手法が記載されている。
すなわちこの論文には3次元の被検査体のスライス部の
面についての線積分を用いて、該スライス部内の小体積
部分の密度を算定するための数式が示されている。また
関連する理論の解析によれば算定は線積分値に係数を乗
算することによつて行われ、係数は小体積部分から各線
までの距離の関数となる。しかしこの論文には数学的理
論を実際に適用した放射線装置は示ざれていない。本発
明の目的は改良した検査装置を提供しようとするもので
、予定の補正係数を用いて画像全体に分散した基本領域
について信号値を順次更新しで行くことによつて当該画
像を形成するものである。本発明の検査装置は、前述の
米国特許第 3778614号の装置に比へ処理時間を短くし、記憶
容量を小さくすることができるものである。
以下図面を参照Lながら本発明の実施例について詳細に
説明しよう。図面を参照すると、患者の頭部を検査する
のに適した形態の装置が示されており、この装置は、そ
の主フレームの一部分を構成している固定ケーシング2
の内側で回転可能な回転部材1を有している。
この回転部材1は中心開孔3を有しており、その中に患
者の頭部を入れるようになされている。中心開孔3は、
封止用フランジ5に固着された可撓性の材料よりなる袋
状のカバー4によつて水密状に閉塞されている。フラン
ジ5は部材1の遠い方の面に対して密封状態であるが回
転可能な関係に保持されている。第1図には袋の断面が
示されている。患者の頭部は開孔3を通じてカバー4の
袋に挿入される。なお、図示されてはいないが、その頭
部を袋内に保持するために附加的な頭置きを設けてもよ
い。検査中患者の身体を支持するために適当な椅子また
はベツトが設けられている。頭部は、開孔3を通じて袋
4に挿入されると、側壁7を有する水溜め6に突入する
。この場合、袋4が頭部を水から分離している。水溜め
6(丸正面を部材1及びカバー4で、側部をプラスチツ
ク壁7で、背面を基壁(図示せず)でそれぞれ閉塞され
ている。壁7と基壁は部材1とともに回転する力ζ カ
バー4はそれのフランジ5とともに静止している。なお
、そのフランジは装置のフランジに固着されている。管
8は、溜めに水を送込み且つそれから送り出すためのポ
ンプに連結されている。患者が袋に挿入されて後に、水
が溜め6に送り込まれて、その袋と患者の頭との間から
空気を追い出す。モータ10によつて1駆動される歯車
9が設けられており、それが回転部材1を駆動してその
部材をしてそれの軸線(開孔3の軸線でもみろ)のまわ
りで軌道走査運動をなさしめる。
歯車9はケーシング部材2の内周面のまわりに形成され
た歯に係合する。回転部材1は透過性放射線源11(こ
の実施例ではX線発生管)を担持しており、開孔3の他
の側には放射線源11に対面してX線検知器12が設け
られている。シンチレーテイング・クリスタル及び二次
電子増倍管よりなる検知器12はコリメータ13を有し
ている。放射線源11は実効的に点放射線源であるよう
になされており、且つその放射線源はコリメータ14を
有している。これらのコリメータ13及び14は、検知
器12に到達する放射線を回転部材1の軸線Lこ直交す
る平面断面内に位置づ゛けられた単一の細いビーム21
にする。その平面は溜め6内にある。放射線源11は、
回転部材1に承支された歯付駆動軸16によつて駆動さ
れる歯付ベルト15に固着されている。ベルト15は軸
16と部材1に承支された第2の軸17との間に架張さ
れている。軸16は可逆モータ18によつて駆動される
ようになされており、そのモータの制御はモータ10の
制御にインターロツクされている−。放射線源11は質
量が大きいから、釣合錘(図示せず一)が設けられベル
トの他の走行部に固着され放射線源といつしよに往復運
動しうるようになつている。この装置の動作時には、放
射線源11とコリメータ14は、モータ18により、回
転部材1の軸線ノに垂直な前述した平面内で横方向の走
査運動をなさしめられる。
検知器12はそのコリメータ13とともにヨーク19に
よつて放射線源11に連結されており、それらの検知器
及びコリメータも同じ横方向の走査運動をなさしめられ
る。その横方向走査時に放射線源とヨークを支持するた
めにガイド20が設けられている。各横方向走査時に検
知器12から出力信号が得られるのであるが、それらの
信号は検査中の平面状断面内における互いに近接する平
行なビーム通路のサンプリングの組に沿つたビーム21
の透過又は吸収を表わす。モータ10及び18間の不ン
ターロツクは、いずれかの方向に一おける各横方向走査
に続いて、例えば1゜の軌道運動の増分がモータ10に
よつて回転部材1に与えられるようになされている。然
る後、モータ18の制御のもとで他の横方向走査が発生
するのであるが、この場合の走査は前の横方向走査とは
反対方向である。近接する平行なビーム通路の他の組に
沿つたビーム21の透過を表わす他の組の信号が得られ
、そのビーム通路の組は前の組に関して1゜だけ角度が
異なる。横方向走査の変位を監視し且つ出力信号のタイ
ミJャOを決定するために、ヨーク19に連結された計数
線(図示せず)と協働するフオトセル装置22が設けら
れている。全体として180゜の軌道運動が完了するま
で、交互に軌道及び横方向の走査運動が継続される。第
1図に示されているように、溜め6の横方向の大きさは
横方向走査のそれに等しく、その横方向の大きさの両端
が点線ビーム21′及び21′lこよつて示されている
。
それは開孔3の両側に突出していて、各横方向走査の開
始点においてビーム21が、或る時間の間、溜め内の水
中を既知の通路長だけ透過するようになされている。か
くして、溜めは、それtこ水が注入されると、被検査体
にビームが達する以前に各横方向走査の開始時において
ビーム21に既知の減衰を与えるように前記位置決め手
段に関して位置づけられた基準減衰器を与える。次に述
べるところから明らかなように、ビームが水溜めによつ
てさえぎられている間に検知手段から基準信号が得られ
、この基準信号は、ビームが被検査体によつてさえぎら
れた時に得られた出力信号を修正するために利用される
。カバー4を除いて溜めが部材1といつしよに回転する
ので、その溜め6の側部分によつて与えられている基準
減衰器中のビームの通路は、角度に関係なく、各横方向
走査について実質的に同一となる。部材1上は鉛のプロ
ツク23も取付けられており、その鉛プロツクは、放射
線源11及び検知器12によつて行われる横方向の走査
の一端に配置されている。鉛プロツク23はX線を実質
的に完全に吸収するものであり、且つビームがその鉛プ
ロツクによつてさえぎられた場合に検知器12から得ら
れた出力信号は第2の基準信号を与えるものであり、そ
の第2の基準信号は、ビーム21が被検査体によつてさ
えぎられた場合のみならず、それが基準減衰器として作
用する溜めの部分によつてさえぎられた場合にも、検知
器から得られた信号を修正するために用いられる。溜め
6は横方向走査運動全体に亘つてビ一仏21に減衰を与
えるが、ビームが被検査体によつてさえぎられやすい領
域では減衰が軽減されることに注意すべきである。水の
減衰又は吸収係数は、被検査体が袋内に位置づけられて
いる間(−勿論、ビームが鉛プロツク23によつてさえ
ぎ一られる時を除き)、ビーム21の全吸収が各横方向
走査全体に亘つてほぼ等しくなるようなものである。検
知器12の出力信号を処理するための回路においては、
各横方向走査の開始点においてビームが既知の減衰をう
けた場合に発生する基準信号の対数が他の出力信号の対
数から差引かれ、そしてその結果得られた出力信号は実
質的に水を透過する場合の減衰と被検査体内におけるビ
ームの減衰の差のみを表わす。基準検知器24はX線源
11に近接して取付けられていて、コリメータ25を通
じて放射線源から放射線を直接受取る。検知器22はX
線のエネルギ―を監視するために設けられている。第2
図に示されている装置の回路は第1図に関して説一明し
た機構の検知器12及び24からはじまつている。
検知器12の出力信号はゲート30に与えられるが、そ
のゲートは主制御回路31からのサンプリング・パルス
により、予め定められた時点で開放される。この主制御
回路はフオトセル装置からの入力信号を受取り、且つゲ
ート30に対してだけではなく、モータ10及び可逆モ
ータ18に対しても適当な制御信号を与える。その主制
御回路がゲート30に与えるサンプリング・パルスは前
述した計数線によつて決定された時間で発生され、平行
なビーム通路のサンプリングの組即ち群を通じてのビー
ム21の透過に対応する一連の出力信号を検知器12か
ら得るようになされている。その通路の組に対する角度
は回転部材1の角度位置によつて決定される。各サンプ
リングの期間の間、検知器12の出力は積分器32で積
分され、そして次にアナログ・デジタル変換器33にお
いてデジタルコードに変換される。各サンプリング・パ
ルス時に発生される信号は記憶装置34にデジタル信号
として記憶される。X線ビーム21は2つの横方向走査
で一回だけ鉛プロツク23によつてさえぎられ、従つて
、検知器12からの対応する出力信号が2回の横方向運
動の期間の間記憶される。記憶装置34内の特定の平行
な組の通路の信号は、−X線ビーム21が溜め6内の基
準通路を通過すべきことがわかつたとき得られた信号と
、被検査体を通過するときに得られた信号とを含んでい
る。ゲート35は、ビーム21が鉛プロツク23によつ
て実質的に−さえぎ一られる時点で検知器から得られた
出力信号を選択するために設けられている。他のゲート
36は他の運動時点において得られた信号を選択するた
めに設けられている。その選択は、鉛による事実上完全
な減衰を表わす基準信号が引算回路37においてサンプ
リング組のそれぞれの他の信号から差引かれるようにし
てサンプリング−・パルスによつて制御され、従つて、
その引算の後では、その結果得られた信号が、基準とし
ての鉛の吸収に関係づけられた被検査体内のビームの透
過又は吸収を表わす。このようにして、検知器1?内に
おける「アフタ一・グロー」の影響はほとんど除去され
る。その結果得られた信号が割算回路38に入れられる
。前述した基準検知器24は、ゲート30に与えられる
サンプリング・パルスと一致した主制御回路31からの
サンプリングパルスを受取る出力ゲート40を有してい
る。ゲート40を通過する信号は積分器41内で積分さ
れ、そして変換器42内でデジタル信号に変換される。
それらの積分器41及び変換器42は積分器32及び変
換器33にそれぞれ対応している。検知器24から得ら
れたデジタル化された信号は次に記憶装置43に送られ
、そして前記割算回路38に与えられる。その割算回路
では、検知器からの各信号が記憶装置43からの対応す
る信号によつて割算され、放射線源11のエネルギーの
変化を補償する。そのように補償された信号は、基準と
しての鉛上で基準化された検知器12からの各信号を対
数に変換し且つそれらをこの形態に保持する対数変換回
路45に送られる。これらの信号は、主制御回路31か
らのパルスによつて制御濾れる2つのゲート46及び4
7に与えられる。ゲート46は、被検査体が位置づけら
れうる領域を特定の横方向運動時にX線ビーム21が通
過する−時点に対応する時点で開放され、他方、ゲート
47は、水の中の基準通路をビーム21が通過しつつあ
り、従つて既知の減衰を受ける時点に対応する時点にお
いて開放される。従つて、“ゲート4−7からの信号は
基準信号と呼ばれうるものであり、ゲート46からの信
号は出力信号として区別される。更に、基準信号は、特
定の走査の各出力信号に一致するように反復的に読み出
され、且つ引算回路48においてそれらの出力信号から
差引かれ、従つてその出力信号が、溜め内の水によつて
生ぜしめられた既知の減衰に対するビーム21の減衰の
比を表わすことになる。検知器12の感度の迅速な変動
に基因する出力信号のスプリアスな変化がかくして実質
的に補償される。これらの修正後における出力信号は信
号処理器50に供給され、他のすべての組の出力信号と
ともに、検査中の体の部分における放射線の吸収の分布
についての映像再構成に関与する。第3図に示されてい
るようlζ前述した実施例においては、第2図に一般的
に示されたゲートおよび積分器回路30,32は、2つ
の積分器32a,32bおよび2つのゲート30a,3
0bよりなつている。
公知の構成のミラー積分器であるこれらの積分器は検知
器12に対する連続入力を受取りかつそれらの入力を積
分するように接続されている。しかしながら、積分器3
2aは、放射線源11と検知器12が距離2aだけ移動
する一定の間隔をもつて存続期間の短いパルス51aを
受取る。積分器32bは同様に存続期間の短いパルス5
1bを受取る力ζそれらのパルスは積分憫2aに与えら
れるパルスの間に発生するようになされている。これら
のパルスは両積分器における読取りおよびりセツト動作
を開始する。それらの動作は、隣接パルス間の時間間隔
に比較して無視できる程度に短い時間間隔を占める。前
記両積分器に供給されるパルスは、読取り動作時に、ゲ
ート30aおよび30bにもそれぞれ与えられて各積分
器をアナログ−・デジタル変換器33に結合する。この
動作のモード(丸 コーリメータ13及び14の寸法を
適当に選定することにより、変換器33によつてデジタ
ル型式に変換された各信号をして、変換器33の出力リ
ードの近傍におけるグラフによつて示されているように
、幅が2aであつて、その幅を横切つてほぼ正弦波状の
強度分布を有する事実上のビームの透過を表示せしめる
作用を有する。コリメータ13および14は横断面が矩
形状であるが、放射線源と放射線自体の性質によつて、
事実上のビーム強度分布はそれの中心最大値からなめら
かに落下する。さらに、前記のグラフから明らかなよう
に;変換器33からの出力信号はaというサンプリング
間隔だけ離間された平行なビームによる被検査横断面内
の吸収サンプリングを表わす。第2図のゲートおよび積
分器回路40,41は第3図に示されているものと同様
の構成を有している。この実施例においては、積分器3
2aおよび32bは、各横方向走査運動において160
またはそれ以上のデータ信号を与えるように制御され、
その場合、このような信号の組が180組存在している
。一組当りのビーム・データ信号の数およびそのような
組即ち群の数は勿論変更しうる力(常に鮮明度の高い映
像を与えるのに十分なだけ大きく選定されている。実際
には、各数は、種々の寸法の映像を再構成できるように
するため、50またはそれ以上であろう。第4−図は第
2図に示された信号処理器50の全体の構成を概略的に
示している。
第4図によれば、処理器50は、ビーム・データ記憶器
239を有しており、引算回路48からの各信号がこの
記憶器に書き込まれ、この場合、前述の信号はデジタル
符号の対数として表わされる。この記憶段においては、
各信号に対する記憶器内のアドレスがアドレス・セレク
タ234によつて選定される。記憶器239においてこ
のような対数記憶が完了すると、アドレス・セレクタ2
51によつてこの記憶器からデータが平行な組ご゛とに
引き出され、そしてそのデータがデータ処理器252に
与えられる。データ処理器252は補正されたレイヤー
グラム(COrrectedlayergrarrl)
と称し得るものを発生するもので、その性質および動作
については、後に詳しく説明する。各組についての処理
が行なわれると、処理されたデータは項ごとに(Ter
m−TO−Tenn)処理済みデータ記憶器即ち合計デ
ニタ記憶器253に記憶される。この記憶器253は;
それぞれ各1つの組から得られるデータだけを受取るた
めの異なる−部分を有−している。本発明の装置はさら
に、すべての処理が完了した場合に、データが、被検査
横断面の領域上の吸収係数の分布を直接表わす形態で保
持される出力マトリクス記憶器254を包含している。
その記憶器254のアドレスは;たとえばデカルトの網
(Cartefsian−MeshwOrk′X!)メ
ツシユに対応しており、ごの場合各メツシユは被検査横
断面の特定の基本領域を直接表わしており、かつすべて
のメツシユは連続的にともに延長して、検査横断面の少
なくとも関心のある領域のすべてを包含するようになさ
れている−。各メツシユのア下レス(ごは、関連するメ
ツシユの基本領域内にある被検査体の物質の吸収係数を
、装置によつて許容される程度の精度をもつて、表わす
信号が最終的に記憶される。す−ベてのメツシユにつき
記憶が完了すると、たとえば陰極線管またはプリンタあ
るいはそれらと組合せたまたはそれらに代替せしめた磁
気テープ型記憶装置を有する装置134によつて表示さ
れうる。映像再構成には高い精度が要求されるから、記
憶゛器253に記憶された処理済みデニタを出力マトリ
クス記憶器254に転送する際に、補間器255によつ
て補間が行なわれる。この補間は;後述するように、゛
アドレス・セレクタ256とビーム通路データ記憶器2
57との協働によつて実現される。信号処理器50につ
いてさらに説明する゛まえに、本発明による補正された
レイヤーグラム技術の理論的な点につき第5図を参照し
て説明しよう。
第5図において、点0は、被検査体の物質の吸収が要求
される探索ビーム21の平面内の任意の点である。いま
詳細に説明しようとしている態様でその点における吸収
を確定すると、その平面内の他の点における吸収も同様
にし7て検知されうるので、かくしてその平面内におけ
る吸収の映像を確立することができるのである。点0は
出力マトリクス記憶器の前記デカルトの網を基準として
いるものとする。第5図においては、その網(メツシユ
ワーク)の軸線0xはX軸の方向にあり、0yはy軸の
方向にある。簡単のために、透過がサンプリングされる
各横方向走査時の一連の位置におけるビーム21を一組
の平行なビームよりなるものと考える。既に述べたよう
に、これらのビームはすべて同様の横断面を有するもの
であり、ビームの間隔は任意の組において均一で−あり
かつすべての組につき同一である。さらに、それらの組
は互いに対して均一に規則的に離間されたシーケンスで
存在−している。このような探索ビームを用いて上記点
0のような点における吸収を確認す−るための不正確な
方法としては、その点0を通過するすべてのビームがう
ける吸収を合計する方法がある。
点0を通過するすべてのビームに共通な比較的小さい領
域の外側における物質による吸収が比較的小さくて無視
できるような特別の状況においては、この和は点Oのす
ぐ近傍における吸収に近似する。しかしながら、実際に
は、点0のまわりの小さい領域の外側における吸収がこ
の方法に誤差を導入することになる。この方法はレイヤ
ーグラム法(1aYerg一Ramtechnique
)として知られているものであり、一定の点に関する1
組のビームのレイヤーグラムはその点を通過するすべて
のビ―ムの吸収の和に比例するものと理解される。本発
明の基礎となる補正されたレイヤーグラム法は、点0に
関する単純なレイヤーグラムをその点0における吸収の
表示として取り出す場合に生ずる誤差が点0を通るビー
ムのその点以外のすべての点における吸収による寄与を
包含しているという原理に基づくものである。さらに、
全誤差が点0を通らないビームの吸収から得られるが、
このようにして得られた誤差は前記単純なレイヤーグラ
ムから差引かれ、これによりそのレイヤーグラムの補正
がなされる。f(θν)が軸0xに関して角度θッをも
つて点Oを通過するビームが受ける吸収であり、かつ異
なる角度をなすこのようなビームが全体でN+1あると
すると、探索される平面のビームの点0に関するレイヤ
ーグラムの値Gは次の式で定義される。
ω(X,y)が探索平面内の任意の点(X,y)(ただ
しXおよびyは観念的なデカルト座標系に関するその点
のデカルト座標である)のすぐ近傍における、探索に使
用された放射線の単位長当りの吸収であるとすると、ω
(X,y)はその平面内の吸収係数の分布を表わす。
さらに、次の式が成立する。ただしDsはビームに沿つ
た距離の要素であり、かつこの積分は添字νで識別され
ているような角度θッにおけるビームの軌導、に沿つて
なされる。
DAが点(X,y)に位置した探索平面の面積の要凛で
ありかつsが点0から点(X,y)までの距離を示すも
のであるとすると、次の式が得られる。ただしこの積分
は原点における特異点を除いて前記平面の全体に亘つて
いる。
明らかなごとく、もし所望されれは、その積分は、前記
平面の任意の数の別々の領域が積分の全領域をカバーし
ていれば、そ右らの別々の領域に亘つてとつたものとみ
なすことができる。特に、その積分は、点0を中心とす
る円によつて境界づけられた第1の領域をカバーすると
ともに、点0に対して同心状の多数の薄い環状体部分を
もカバーしておりかつ面積として前記平面の残部のそれ
を総計したものとみなされうる。点0から第k番目の環
状体部分についての積分は次の式で与えられる。ただ−
しSkは点0からのその環状体部分の平均距離である。
これは、その環状体部分力叶分に薄いと仮定している。
この積分は点0を通過するすべてのビームの、それらの
ビームが環状体部分kを通過することに基因する吸収の
和を表わしている。ω(X,y)が第1の領域の全体に
亘つて略同じ値ω。を有していると考えることができる
とすれは、この領域についての積分は次の値をとる。た
だしaはその領域の直径であり、すべてのピ゛ームは対
称的に点0を通過−かつすべてのビームがaに比較して
小″さい幅を有しているものとする。しかしながら、ビ
ームがaなる幅を有しており、放射線の強度がそれらの
ビームの横断面にわたつて感知できる程度に均一である
とすると、αは次のような異なる値をとる。この場合、
ビームは探索平面に直交する方向に単位厚さを有するも
のとみなされている。
いずれにしても、もしであわば、Gに対して次の式が得
られる。
― V 点0を通る探索ビ々ムの通過に関する対称性の仮定は、
点0が探索平面内の異なる場所に選定されるから、一般
的に完全に成立する−わけではない。
しかしながら、近似的には、後でさらに詳細に説明する
ように、その仮定は適当な補間法によつてすべての場所
に対して実効的に成立せしめられうると一般的にいえる
であろう。ω(X,y)がaという大きさの距離ではど
こでも感知できる程度には変化しない(このことはω(
X,y)が周波数帯域を制限された関数であることを意
味する)と一般的に仮定すると、幅a(この幅は薄さの
条件に従う)を有するすべての環状体部分をとりだすこ
とができ、したがつてと書くと、Gに対する式は次のよ
うになる。
このようにして得られた式において、αという項は吸収
関数の所定の値ω。の大きさ、すなわち点0における値
である。他の項、すなわちkに関する総和の項はωoの
値を確定するレイヤーグラム法の誤差を表わしている。
しかしながら、その誤挙は、探索の他のビームのすべて
、すなわち項αに寄与する中央の円形領域を通過しない
ビームが受ける吸収を考慮することによつて打消されう
る。補正されたレイヤーグラム法の詳細を展開するため
には、直径aの中央円形領域を通過しない残りのビーム
を異なるクラスに分割する必要がある。
帯域kはインターセプトされ得ないことがわかるであろ
う。したがつて、さらに、クラス1のビームが受ける全
吸収はであるが、クラス2のビームが受ける吸収は であり、 クラス3のビームに対する全吸収は であり、以下同様にして各クラスにつき与えられる。
であれば、 次の式 (たマし であり、かつL,で表わされる1組の係数はノである。
)の値は単に値αとなる。
係数Lttの組は明らかにこれらの式によつて決定され
る。
これらの方程式は、インターセプト係数(その典型的な
ものがFJkである)が決定されれば、Ll,L2,L
3・・・・・・について解くことができる。
極限士1におけるL,SfiO)ような項の総和を次に
得て、点0における吸収関数の値を得ることができる。
係数Lァは帯域重み付け係数と称し得るものであり、単
にL係数とも呼ばれ得るものである。なぜならば、たと
えば、Lァは帯域μによつて区別されるクラスμのビー
ムが受ける吸収に掛けられるものだからである。原理的
には、インターセプトすなわちF係数が正確にわかつて
いれば、L係数を正確に計算することができる。
F係数はビームの幾何学形状のみならず、そのビームの
横断面上における放射線の強度の分布に関してなされた
仮定にも依存する。前述のように、放射線の強度は、ほ
とんどの場合において「走査開口]効果があるから、ビ
ームの幅方向には一般的に均一ではない力( ビームの
厚さ方向には略均一であると考えることができる。実際
的な目的のためには、F係数、したがつてL係数は、走
査系統の物理的なパラメータに依存する異なる仮定に基
づいて種々の態様で十分に近似されうる。それらの態様
のうちの幾つかについて説明する。第5図においては、
帯域0が点0を中心にして描かれた円周囲の最も内側の
ものの中の領域として示されており その帯域は前述の
ように−aにS2等しい周囲半径を有している。
帯域0のすぐ外側にζ人平均半径aとaに等しい幅とを
有する第1の環状領域すなわち領域1が存在している。
帯域k以外の他の周囲領域は示されていないが、その帯
域は平均半径Kaを有しており、かつすべての環状領域
と同様に半径方向の幅aを有している。クラスOの1つ
のビームは、それが非常に薄いビームであると仮定する
と、第5図における0y軸によつて表わされているもの
と考えられうる。このビームE。FOで表わされている
長さで示されているように第k番目の帯域をインターセ
プトしている。同じ薄い横断面寸法を有するタラスjの
ビームも示されているが、これは0y軸からJaの距離
のところに引かれたその軸と平行な破線によつて表わさ
れている。このビームはEFで表わされた長さで示され
ているように第k番目の帯域をインターセプトしている
。これらの場合に6ζインターセプト係数Fjkはイン
ターセプトEFの長さをインターセプトE。FOの長さ
で割つたもので与えられる。この比は、クラスOのビー
ムに比較した場合の、帯域kを通過するときにクラスj
のビームが受ける吸収の増加を表わしている。この比が
有効であるためにはビームが無限に薄くなければならな
い必要はないが、F係数を近似するこの方法は、隣接し
た平行ビームの分離によつて決定される環状体の幅aに
比較して実効ビーム幅が小さければ、L係数に対して良
好な結果を与える。ビームが均一な横断面強度を有しか
つ環状幅aに等しい幅を有すると仮定しうるとするなら
ば、F係数は次の根拠に基づいて計算されうる。
第5図において、クラスOのビームは、帯域0の周囲円
に正接する平行な線P。QOおよびS。ROによつて境
界づけられて示されている。クラスjの同様なビームの
一例力ζそれぞれ座標距離(j+2)A.lおよび(j
−)aにおいて0x軸と交差している境界平行線PQお
よびSRによつて表わされている。
上記ビームのうちの最初のものが領域POQOROSO
上において第k番目の帯域をインターセプトし、2番目
のビームはインターセプト領域PQRSを有している。
このようなビームの場合にtl人インターセプト係数F
jkの値は後者の領域PQRSの面積と前者の領域P。
QOROSOの面積との比である。これは、前述したよ
うに、ビームの強度がそれらの横断面に亘つて略均一で
ありかつ吸収の程度が十分に小さいという仮定のもとに
いえることである。補正されたレイヤーグラム法による
吸収データの処理への他のアプローチにつき第6図を参
照して説明しよう。
これはここに記載されかつ図示されている装置において
採用されている方法である。第6図においては、第5図
におけると同様に、点0は被検査体の物質の吸収が求め
られているものと仮定されている探索ビームの平面内の
点である。ここでも、この点は、前述のデカルトの網を
基準として位置が特定されているものと考える。なお、
第6図における0x軸はその網(メツシユワータ)のX
軸の方向にあり、0y軸は網のy軸の方向にある。第6
図は、軸線が0y軸に平行あるいは略平行である−よう
な角度を有するビームの平行な組の中の5つのビームA
B,CD,CD′,EFおよびE/F′を示している。
まず、これらのビームはずべて同一の横断面を有しかつ
幅aを有していると仮定しよう。さらに、ビームABの
中心線は軸0yであると仮定する。他のビームの中心線
は破線で示されている。また、点0と同心状に描かれた
円13C0,C1およびC2はそれぞれ−2a,−2a
および−aなる半径を有しており、かつそれらの順序は
この通りであると仮定する。
第6図は1つの特定の平行な組のビームだけしか示して
いないが、他のすべての角度を持つビーム組も存在して
いることを考えなければならない。
かくして、もしN+1がこのような平行な組の総 .数
であるとすると、1つのビームのみが円C。の領域すな
わち帯域0にインターセプトしている如く図示されてい
るが、N+1のビームが同様にして(角度の差は別にし
て)この帯域にインターセプトしていることを想起しな
ければならない。数1N+1は1つの組当りのビームの
数のπ倍のオーダを有するように選定されうるものであ
り、映像再構成に関してはそのようにするのが有利であ
る。このように仮定されたインターセプシヨンは、帯域
0をインーターセプトするすべてのビームの中心線が点
0を通過するという意味において理想的である。探索平
面内において点0の選定を変えれば、隣接ビームの吸収
データ間の適当な補間によつて、この理想が一般的に近
似されうる。円C。
の外側で円C1の内側の領域をここでもZ帯域1と称し
、帯域0をインターセプトしないすべてのビームのうち
、2(N+1)のビームが帯域1をインターセプトする
。これらは、ここでもまたクラス1のビームである。帯
域0をインターセプトするビームに関してなされた理想
的な仮定に基づいて、クラス1のビームの中心線は、中
心が点0で半径がaの円に接するであろう。さらに、平
行な組が均一な角度配置関係にあるという仮定に基づい
て、ビームの組の角度配置は均一な角度系列のそれ即ち
であり、かつ接点はこの円のまわりに均一にかつ対応し
て分布しているであろう。
円C1の外側で円C2の内側に存任する領域をここでも
帯域2と称する。
クラス2のビームは、クラス0および1のビームを除き
、帯域2をインターセプトするすべてのビームよりなつ
ている。このクラスのビームもその数は2(N+1)で
あり、かつクラス1のビームと同様にして、半径2a1
中心0の円に関して理想的に分布されている。探索面内
にあるビームの分類も、第5図の場合と同様に、このよ
うに継続され、而して順序的に分類される一連のビーム
・クラスすなわち帯域0,1,2,3,4・・・・・・
に対応するクラス0,1,2,3,4、・・・・・・を
確立する。
すでに説明したように、探索平面内の吸収領域のみが帯
域0内に位置しているならC丸クラスOの種々のビーム
のすべてが受ける吸収の和がこの吸収領域の真の吸収値
となるであろう。
しかしながら、実際の場合には、それらのビームは帯域
0の外側の他の吸収領域を通過するから、クラスOのビ
ームの吸収の和は正確な目安とはならない。しかしなが
ら、対応するL係数により全帯域ビーム吸収を重み付け
しかつ加算するという前述した一般原理を適用すること
ができ、かつその場合には、次のようにしてL係数を近
似することができる。ビームの平面内における被検査体
の物質の吸収が、ビームの幅寸法aに比して非常に小さ
いと考えることのできる程度の小さい寸法の単一の領域
(この領域は帯域0内にはなくて他の任意の帯域内にあ
る)を除くあらゆる場所でゼロであると仮定する。
このように仮定すると、L係数が正しく選定されている
場合には、重み即ちウエイトをつけられた吸収の全総和
の加算は、その領域がどこに位置していても、ゼロとい
う結果を与えるべきである。このことから、L係数を解
として与える1組の方程式が得られる。前記小さい領域
の場所即ち位置が与えられたとき、全総冴田こ吸収を寄
与せしめる任意のクラスのビームの数はその領域をイン
ターセプトするそのクラスのビームの数であることがわ
かるであろう。
ビームはその横断面に亘つて略均一であり、かつ吸収の
程度が十分に小さいと仮定すると、重みを付けない吸収
の総和は領域によつてインターセプトされるビームの数
に単純に比例するであろう。たとえば、その領域が帯域
2内にあるならば、その領域がそのように位置決めされ
ている場合において全総和に吸収を寄与せしめるクラス
lのビームの数は、その領域に重なるすなわちその領域
をインターセプトするクラス1のビームの数であり、か
つクラス1のビームに対するL係数による重み付けの前
における吸収寄与はこのビームの数に単純に比例する。
これによつて、中心0に対して同心状でかつ前記小さい
吸収領域を通過し、クラス1のビーム内にある円弧長に
比例するような重み付けされない吸収寄与の式が得られ
る。これに基づいて、中心0からの小さい吸収領域の半
径方向の距離についてある仮定がなされれば、L係数を
得るための方程式を導出することができる。たとえば、
前記小さい領域がどの帯域内に位置していようとも、そ
れはその帯域に対する平均半径のところにあると考える
こととしてもよい。小さい領域の位置として帯域の両極
端の半径の間の或る他の半径が選択されたとすると、イ
ンターセプシヨンの円弧長が異なり、かつそれに対応し
て、平均半径に位置するという仮定に基づく場合と同様
ではあるが同一ではないL係数が得られることがわかる
であろう。
このような差異を説明するに当り、第6図に示されてい
る半径2aの円弧JKLlV7)象限について考えると
ともに、前記小さい領域がこの円弧上の帯域2内にある
と仮定しよう。
クラス2のビームは、円弧JKの長さ即ちクラス2のビ
ームによつてインターセプトされる円弧JKLMの長さ
に比例する重み付けされない全吸収与与をなすであろう
。同時に、クラス1のビームは、円弧KLの長さ即ちク
ラス1のビームCDによつてインターセプトされる円弧
JKLMO)長さに比例する重み付けされない全吸収寄
与をなすであろう。クラス1のビームは、円弧LMの長
さ即ちクラスOのビームの半幅によつてインターセプト
される円弧JKLMの長さに比例する対応した寄与をな
すであろう。そのインターセプトとして半分の長さを取
るのは、クラスOのビームの数が他のすべてのクラスの
数の半分にすぎないという事実を考慮したためである。
以上のような理由づけに基づいて、円弧長JKlKL,
LMがF係数Fl2およびF22の値ならびにすべての
F係数に対するような値を決定する。しかしながら、前
記小さい領域が平均半径0J以外のある中間帯域半径、
たとえば円弧J″K″L″ぱの半径0J′(ただしJ!
,K,L′およびぱはインターセプトJ,K,L,Mに
対応する新しいインターセプトを示している)に位値し
ているならば、新しい円弧長J/K′,KL′,Cイ(
まF係数Fl2およびF22を異なる値に決定し、たと
えL係数の値が同様であつても、結果は異なるであろう
。実際には、前dウ」\さい領域が平均帯域半径に位置
する場合には、一連の帯域係数L,,L2,L3,・・
・・・・は、境界と遷移領域を強調する傾向のある再構
成された映像を生ずる。
他方、もし前記小さい領域がより大ぎい半径に位置して
おれぱ、たとえば、もしその領域が帯域の平均半径と外
側半径との中間に位置しておれば、この強調は事実上除
去される。そのように平均半径と外側半径との中間に位
置している場合には、低次のL係数は次のよう・な値を
とる。一方、前記小さい領域が平均帯域半径にあるなら
ば、最初の3つのL係数は1.000、−0500、+
0.7060となる。
1つの環状部分内の小さい領域の半径方向の場所を、そ
の環状部分の内側境界からの前記小さい領域の半径方向
の偏位の、内側境界および外側境界間の半径の差に対す
る比として定義づけられたβで規定すると、最後に述べ
た系列はβの値が0.50に等しいそれとして説明する
ことができ、一方、最初に述べた系列はβを0.75に
等しく設定することによつて得られる。
関係する帯域の全範囲に亘つてβの値が一定である必要
はない。たとえば、帯域の順位が高くなれば、そのβの
値は漸次減少するもの、例えば帯域1に対しては0.7
5の値を有し、然る後、それが最高次の帯域の場合に有
する0.5という値まで漸近的に減少するものであつて
もよい。第5図および第6図についての以上の解析は、
検査の対象となつている平面の吸収係数の分布をあられ
す関数ω(X,y)が周波数帯域制限されているので任
意の1つの横方向走査時における薄いと仮定されたビー
ム21の吸収の変化が、距離aに関係付けられた特定の
空間周波数をこえる割合では発生しないという仮定のも
とでなされた。
実際の場合には、勿論、関数ω(X,y)はそのような
程度には必らずしも帯域制限されておらず、かつサンプ
リング間隔aの2倍に等しい波長に相当する周波数以上
の空間周波数が存在する場合には、処理時に誤差が発生
して、再構成された映像に実質的な誤差を生ぜしめるこ
とになりうる。しかしながら、第1図、第2図および第
3図に示された本発明の実施例においては、このような
誤差を減少させるようにしてサンプリングが実施される
。前述したように、サンプリングは、隣接ビーム間の間
隔の2倍に略等しい幅を有する事実上のビームによつて
なされ、隣接ビームの間には大きなオーバラツプが存在
する。さらに。実際上のビームはその幅方向に略正弦波
状の強度分布を有し !ている。これが所要のデータに
周波数帯域制限を導入することを示すことができる。正
弦波分布は、次のような形(ここでaはサンプリング間
隔であり、)で表わされる強度分布を有する中央の正の
領域に対する近似を与えるものであり、かつそれは、こ
の近似が該当する程度において、高い周波数を有する望
まれない成分を除去するために被探索データを帯域制限
する作用を有している。
前述したところで仮定されたようにaが帯域幅であり、
かつXがビーム中心線からの距離をあられすとすると、
の範囲においては、第1図、第2図および第3図による
ビーム強度の分布は近似的に次のごとき形式をとる。各
ビームの通路を横切る吸収物質の適当な分布の導入によ
り、ビームを上記の形または他の所望の形に精密に合致
させ得るが、実際通常の目的のためには、正弦波状であ
ると仮定して適当なL係数5を計算するに際し、この附
加一的な方法を利用しなくても十分な精度が得られるこ
とをが実際に認められた。
βの値を0.625に等しい値にとると、上述のように
ほぼ正弦波状の放射線強度分布を有するビームの場合に
は、再構成された映像に十分満足しうる程度に高い精度
、たとえば各映像点において1/1000のオーダの精
度が得られることが認められた。
βをこの値にすると、映像再構成に転換即ちエツジをと
もなう或る程度のオーバーシユートが存在するが、その
程度はほんの僅少である。系列の最初の200L系数が
第7図の表に示されている。この系列は、再構成を受け
ている映像の基本領域の位置に関係なく、その映像の境
界までの吸収を補償するのに要するだけの個数の係数ま
で継続される。第7図の系列においては、係数L4Oは 0.0001より少し大きいが、たとえそのように非常
に小さくても、この値の係数ならびにそれよりも小さい
値の系数は、それらが多数(いま考えている例の場合に
は約180)のビームの受ける吸収の和に掛けられるの
であるから、正確な映像再構成にとつて重要であること
に変りない。
般の画像を得るための比較的不正確な映像再構成の場合
でさえ、L係数の数は各組におけるビームの数の約半分
でなければならない。一般的に言つて、かつ系列の最初
のほうの項に影響を及ぼす小さい動揺を別にすれば、β
が0.625に等しい場合の係数(LOは別として)は
、常に負であつてかつ漸近的にゼロに近づき、最初の数
項の後にゆつくりと変化する単調関数に適合する。上述
した動揺の種類は振動成分であり、そればL4が数値的
にL5よりも小さいという事実によりβ−0.625に
つき立証されている。βが0.75に等しい場合にはそ
れに対応する動揺はみられないが、βが0.5に等しい
場合には一層顕著にさえなる。このような動揺はL5以
降のL係数に対しては特に影響を与えない。一般的な特
徴として、ビーム幅が隣接した平行ビームの間隔に比較
して小さいかあるいはその間隔に匹適しうるかに関係な
く、L係数(LOは別として)は負の単調関係に近似す
る。前述のようにL−ゅ−L。という原則について厳密
にいえば、この関数はL。に対して対称であり、かつ両
半分が恐らく前述した小さい動揺は別として、LOのま
わりで単調である。この値の変化は開口補正と同様の作
用を有することが明らかである。
かくして、0.5に近接したβの値は、被検査体の吸収
のステツプ関数に応答して映像再構成に著しいオーバー
シユートを与える傾向がある。ステツプ関数の近傍にお
ける数値精度は或る程度低下することにはなるが、細部
を描写するという見地Z)ら、即ち横断面内の或る点が
吸収の急激に変化する部分のすぐ近傍にある場合にその
点の近傍における吸収の大きさを特定の限界内で決定す
るという見地から、所望の程度のオーバーシユートが望
ましいことが多い。他方、一般原則として、βの値が0
.75を大きく超えると、アンダーシユートに相当する
条件が存在する傾向のあることが認められる。βの実際
的な範囲は0.5と0.75との限昇値間にあり、まれ
なケースではあるが、特別の効果が所程される場合には
、βはこれらの特定された両限岑値の外側にあつてもよ
い。本明細書で述べられている他の近似に基づいて発生
されるL係数を適当に調節することにより、開口補正も
シユートされ得る。以上レイヤーグラム補正について述
べた一般的な手法は、帯域ごとになされる一連の近似の
連続と見ることができる。
たとえば横断面の物質の吸収が帯域1の外側ではゼロで
あるとすると、帯域O内の吸収に対する単純なレイヤー
グラム評価(数値を定めること)における誤差は、係数
L1にクラス1のビームの全吸収を掛けたものにより補
償される。他方、有限の吸収が帯域2の外側境界まで延
長しているならば、L1補正は不干努でり、それは同一
のパターンについてのL2補正によつて補足されなけれ
ばならない。このようにして継続しかつ帯域を次から次
へと考慮に入れると、この方法は各帯域が計算に導入さ
れるにつれて真】の結果に常に一層近密に接近する方法
であると一般的に考えられ得ることが明らかである。第
n番目の帯域の外側で吸収が消える場合には、真の結果
を得るための係数Lnまで一連のL係数補正が継続され
なければならない。再構成された映像のう境岑の外側で
は吸収は勿論ゼロであると仮定されている。第7図に示
された系列を実際に決定する場合には、正弦波状分布の
各所定ビームは、そのビームの全幅に亘つて延長してい
る多数の隣接した、比こ較的細いビームと同じ効果を有
するものと仮定された。
各細いビームはそれの幅方向に均一な強度を有するもの
と考えられたが、そのビームの強度は前記所定のビーム
の強度分布曲線上の一連の点に対応すると仮定された。
概略的に述べた理論に5より、各場合においてβ−0.
625であると仮定して、L係数が、対応する細かい帯
域に附与された。これらの係数にも各細いビームの強度
が掛けられて一連の寄与が発生され、これらの寄与は加
算されて正弦波状の強度分布を有する上記所定の4ビー
ムのそれぞれに対するL係数が得られる。第1図および
第2図に関して前述したように、特定の系列のL係数を
用いることを決めかつ平行な組をなすビーム吸収データ
を対数の形で得ると、第3図のユニツト252によつて
実施されるべき処理が適当なプログラムを与えられた計
算機によつて行なわれる。上述の理論によれば、物質の
吸収係数が評価されるべき基本領域における所定の点に
関連して定義された所定のクラスのすべてのビームが受
ける吸収が総計され、かつその総計に各L係数が掛けら
れることに注意すべきである。すべてのクラス(前記所
定の点に関連して定義された)に対するこのような重み
をつけられた総計の和は所要の吸収係数に比例する。こ
れらについての手法は原理的には前述したラマチヤナド
ラン他の論文に記載された処理に類似している。しかし
ながら、本発明の理論的な基礎は前記論文に記載されて
いるものとは異なるものであり、本発明に従つて用いら
れるL係数は、前記論文において到達した関数とは実質
的に異なつた関数に従つている。また、それらの係数は
実質的に改善された結果を生ずる。さらに、本発明の特
徴によれば、各組のすべてのビーム吸収値が個々にL係
数系列を掛げられるのである。このようにして得られた
積は処理済データ記憶器253の異なる区画における異
なるアドレスに分配即ち仕分けされる(この場合、各平
行な組のビ・−ムに対して1つの区画が存在しており、
かつ各区画内の1つのアドレスが各平行な組内の各ビー
ムの中心線に対応している)。各アドレスは1つの積の
系列の全部を蓄積する。記憶器253に蓄積された総計
は出力マトリクス記憶器254内の異なるアドレスに分
配され、検査領域における吸収分布に関する所要の表示
を形成する。評価点を一度に1つずつ取り出す基本的に
二次元的な処理法に代えて、すべての評価点に関するす
べての組のすべてのビームに対する一連の線型処理およ
びこれに続く選択的加算を行なえば、動作速度の点で特
別な利点があり、かつまた線型処理を実施するための次
に述べるような特殊目的用の計算機が使用できる。この
特殊計算機は動作速度を上げる作用をもする。また、こ
の方法は、出力マトリクス記憶器254内のアドレスで
表示されるデカルトの網の各メツシユに対する吸収を評
価する場合の補間を容易にする。第8図を参照して、ま
ず線型手法のシーケンスについて考えてみよう。第8図
において、数字239はビーム吸収データが前述した平
行な組ないし群をなして記憶される記憶器を示しており
、記憶器は第4図に示すように対数の形のデータを受入
れる。
記憶器239から出力を取り出すためのアドレス・セレ
クタ251ならびに、線型処理法に従つて動作している
と理解されるべぎ処理器252からの処理されたデータ
を保持する処理済データ記憶器253も同様に示されて
いる。記憶器253の使用について考えるために、1つ
の平行な組の探索ビームに関する吸収データ値が、それ
ぞれそれらのビームがその組を横切つて発生する順序に
従つて、EO?E1ツE2ツE3?●◆●●●●●●●
●●●●●●のごとく考えられると仮定する。
また、次の項LO,Ll,L2,L3,・・・・・・・
・・・・・・・・の列がL係数列を構成する列であると
仮定するとともに、値E。が記憶器239から引き出さ
れていると考える。この値にはデータ処理器252内に
おいてL係数列の各項がそれぞれ別々に掛けられ、かつ
その結果得られた積EOLO,EOLl,EOL2,E
OL3,EOL4・・・・・・・・・EOLO)がアド
レス・セレクタ251の制御のもとで処理済データ記憶
器253の区画253′内のアドレスAOツA1?A2
7A3ツA4゜O゜゜OO′1″?Aωにそれぞれ記憶
される。
次に、値E1が記憶器239から取り出されそしてそれ
に同様にL係数列252の各項が掛けられて次のような
積が得られるものと考える。ElLOツElLl?El
L2?ElL3゜゜゛゛ElLω−1これらの積は記憶
器253の区画253′におけるアドレスAl7A2?
A3?A4ツ゜1゜゜03゜゜゜?Aωにすでに記憶さ
れている積に加算されそれらのアドレスに記憶される。
次に、記憶器239から値E2を取り出しかつ前と同じ
掛算の手法に従つて、次の積E2LOラE2Ll?E2
L2ラ゜゜゜゜゜゛E2Lω−2がアドレスA2,A3
,A4,・・・・・・・・・,AOにそれぞれ加算して
記憶される。
前記平行な組のすべての吸収値が使用されるまで同様に
してプロセスが続けられる。すべての吸収値が使用され
た時には、アドレスA。に記憶された大きさはSO−E
OLOであり、アドレスA1のものは S1=ElLO+EOLl であり、アドレスA2のものは S2−E2LO+ElLl+EOL2 であり、アドレスA3のものは S3=E3LO+E2Ll+ElL2+EOL3であり
、それ以降のアドレスのものもこのパターンに従つたも
のとなる。
この手法の次の段階においては、上述と同様にしてでは
あるが、いままでとは逆の順序で記憶器239から吸収
値が引き出される。
そのように逆順に取り出される場合には、それらの吸収
値は、記憶器253の区画2535におけるアドレスも
逆順序で使用されるという点を除けば前述したのと同一
の態様で処理される。これは第8図に示されている。最
後に、次の値 EOLO?ElLO7E2LO7E3LO7E4LOF
゜゛゜゜゜゜E0)LOがアドレス AOラAl7A2ツA37A4゜33″″′YAωに既
に蓄積されているものからそれぞれ差し引かれて、これ
らのアドレスに最終的な1組の値を生ずる。
この手法はすべての平行なビームの組に対応する吸収値
に対してなされ、平行なビームの各組に対する値の組は
記憶器253の異なる区画に蓄積される。
所要の映像再構成を与えるための蓄積された値(大きさ
)に関する好ましい最終処理について次に説明しよう。
値(大きさ)を蓄積するための上述したプロセスは、各
ビームデータ信号につき合計信号を得るプロセスと呼ぶ
ことができよう。このプロセスはレイヤーグラフを形成
するために用いられた場合に補償されたレイヤーグラム
を発生する。第8図に示された一連の動作は適当なプロ
グラムを与えられたデジタル計算機によつて行なわれ得
るが、本発明のこの実施例では、L係数ダータ処理器は
、特殊な計算機の形態をなす。
この計算機では、記憶器253に転送されるべき積を発
生する掛算をも実施する回路によつて発生される時系列
の項に内在的にL係数が発生される。これは、L係数が
一般にゆつくりと変化する単調関数に近似するものであ
り、従つてL係数のすべてまたはほとんどが有限個の指
数項の組として近似されうる故になされうる。第9図は
このような考え方に従つて動作する回路の一般的な形を
示すものである。この図において、210はアトkス・
セレクタ251によつて第8図のビーム・データ記憶器
239から引き出されたヂジタル型のビーム吸収データ
Ejを運ぶ通路である。
この通路21からそのデータが、破線で示されたプロツ
ク212?よび2122内に含まれているような複数の
同様の回路の入力に通路211を通つて分配される。こ
れらの回路からのデジタル型の出力が加算回路214を
通じて共通の出力通路213に供給される。通路213
上に現われる変換されたデータは、第8図に関連して概
説された手法に従つて第8図の記憶器253の任意の区
画例えば253′に転送されかつ記憶されるのに適した
ものである。これらの寄与の発生および性質は2121
および2122のような回路についての説明から一層明
らかとなるであろう。回路2121をこれらの回路のう
ち代表的なものと考えると、通路210に現われる吸収
データ信号Ejは通路211を通つて加算回路2151
に供給され、そしてデジタル遅延線2161の入力に供
給され、然る後、その信号は、すべての回路2121,
2122,・・・・・・を制御するために用いられる主
クロツク・パルス源のクロツク・パルス間の期間に等し
い遅延を有して遅延線2161から出て来る。
このようなりロツク・パルス源を用いることは公知であ
るから図示しない。最初の吸収データ信号はそれの通路
210上への出現が1つの平行な組に対する最初の即ち
基準のクロツク・パルスに同期されるように調時される
。これが前述した信号E。であるとすると、その平行な
組の他の信号El,E2等は一連のクロツク・パルスと
同期して順次通路210上に現われる。遅延線2161
から出て来た信号は掛算回路即ち係数回路2171に与
えられる。その掛算回路は入力信号に予め定められた1
よりも小さい数を掛け、そしてそのように掛算された信
号を加算回路2151の2つの入力のうちの他方の入力
に与え、然る後、その信号がデジタル型遅延線2161
の入力に再び与えられる。次にその信号はループ216
1,2171,2151を再循環し、平行な組を形成す
るデータ信号の数、従つて利用されるL係数の数より1
だけ少ない数に等しい回数だけ再循環をくりかえす。そ
の回路の出力は加算回路2151の出力から取り出され
、そして導体2181を通じて予め定められた定数を掛
ける掛算回路即ち係数回路2191に与えられる。掛算
回路2191の出力は通路2201を通じて加算回路2
14の入力に与えられ、回路2122等の出力と加算さ
れる。
一連のクロツク・パルスの発生にともなつて、信号が加
算回路2151の出力に順次現われ、一連の信号が通路
213を通じて順次取り出される。このような一連の動
作が、利用されるL係数の数によつて定められる時間継
続した時、回路2121,2122,・・・・・・の遅
延線をタリアすることによつてその動作が終了され、次
に、新しい系列についての動作が開始される。回路21
212122,・・・・・・・・・のグループの動作に
関して、それらのうちの1つ、たとえば212kの動作
について考えよう。前述のように、グループ出力導体2
13における一連の信号の発生を開始させるときにその
グループに与えられる吸収データ信号はE。という値を
有している。そこで、まず、信号E。のその後の経歴の
みについて考えると、信号E。が与えられてからr個の
クロツク・パルスの期間後における出力導体220kを
通る回路212kの出力は、なる値を有する。
ここで、であり、Akは回路212kにおける掛算回路
217kに設定された1よりも小さい掛算係数であり、
そしてBkは回路212kにおける掛算回路219kに
設定された掛算係数である。
回路212kを典型的なものとするn個の回路が存在し
ているとすると、時間t(r)においては、その回路に
信号E。を印加したことに基因する導体213における
出力(の成分)は次の値を有する。かくして、第9図の
回路は信号E。の印加に基づいて時間t(r)において
一連の信号E。(r)を発生せしめられる。Akおよび
Bkのような乗算係数を適当に選択することによつて、
一連のクロツク・パルスと同期した信号E。(r)の系
列が、少なくとも数個の初期項を除くすべての値に対し
て、高い精度をもつて次の系列EOLO,EOLl,E
OL2,・・・・・・,EOL,,・・・・・・に近似
せしめられうることがわかる。
ただしLは系統に対する一連のL系数である。系列E。
(r)の5うちの数個の初期項、すなわちたとえばE。
(1),EO(2),EO(3)は、初期系数Ll,L
2およびL3が十分正確に近似されていないから、所要
の値EOLl,EOL2,EOL3に十分正確には対応
しないかも知れないが、調節1によつてそのような対応
を得ることができるものであり、このような限定された
調節はデジタル計算機により簡単に実現できる。
212kのような回路が5つ用いられており、実際には
、n=5であれば、上述のように補正されうる最初の数
個の・項を除いて、すべての項が約1/1000という
最終的な映像再構成における精度に相当する精度をもつ
て発生されうることが認められた。
前述のように、第9図の回路を使用した場合には、平行
なビーム組の吸収データ信号が順次第9二図の回路の入
力に与えられる。
この場合、新しい値は各クロツク・パルスのタイミング
をもつて与えられる。その回路がそのように動作される
と、各クロツク・パルスにともなつて、前述したSOフ
S1ツS2FS32llOのような系列の項が順次出力
される。
1つの組のすべての吸収データ信号が入力されかつ遅延
線がクリアされた時、同じ組の吸収データ信号が逆の順
序で供給され、他の同様な系列を発生する。
最後に、EOLO2ElLO?E2LO?E3LO?E
4LO?゛0″11″に比例した信号が記憶器の鋳応す
るアドレスに供給されて引き算される。
これらの演算動作はすべて、第9図の回路が一部分をな
しているデジタル計算器によつて制御されかつ実施され
る。かくして、関係する記憶器の区画、たとえば253
′がマトリクス記憶器254への転送のための最終的に
満たされ、第9図の回路におけるデータの再循環によつ
てステツプ・バイ・ステツプ蓄積が実現される。Ll,
L2,L3という値の発生しうる誤差については、それ
らは、大きい場合には、次のように処理することが好ま
しい。
フアントムを検査して、検査結果を期待値と比較するこ
とにより得られる上記誤差がそれぞれΔL,,ΔL2,
△L3 であるとすると、それらはデジタル計算機に記憶される
。
また、或るビーム・データ値Ejが第8図の回路に与え
られたとする(添字jはデータが1つの平行な組の第j
番目のビームに対応するものであることを示す)。次に
、計算機により補正値EjΔL1 が、その組に対応する記憶器区画のアドレスAj−1お
よびAj+1の内容に加算されて、それらのアドレスに
記憶される。
さらに、補正値E・△L が同様にしてアドレスAj−3およびAj+3に供給さ
れる。
この手法の原理は、その組のすべての値およびすべての
組に対して適用される。第9図によるデータの組の処理
は、吸収データの平行な組を順次取り出しかつそれらを
第9図におけるような再循環型式の単一の処理器に与え
て、直列モードで実施されうるが、再循環型式の処理を
複数の平行なデータの組に1つずつ与えかつそれをプロ
グラムすることによつて、複数の平行なデータの組を並
行に同時的に処理することもできる。
これにより(必要に応じて、)データ処理を非常に迅速
することができる。本発明は第9図のデジタル型遅延線
を用いて実施されるのが好ましいが、そのような回路を
実際に用いなくても、それらをシミユレートする標準的
な計算機を用いることによつても実施されうる。
この場合の処理時間の節約ほそんなに大きくはないが、
それでも有用である。第10図は、第8図を参照して説
明した値、例えばアドレスAO,A,,・・・・・・に
最終的に記憶された値を記憶器253から出力マトリク
ス記憶器254(第4図)に転送するための回路を示す
プロツクダイヤグラムである。
この転送時に、第4図に関連して説明した補間が実施さ
れる。この補間は正弦法則に従う。線型補間を用いるこ
とも可能ではあるが、ビームの幅方向の放射線強度分布
がほぼ正弦波状をなしていることを考えて、正弦法則を
選ぶほうが好ましい。第10図において、プロツク25
2はここでもL係数処理回路を表わしている。
処理されるべき吸収データがN+1の平行な組に入ると
すると、処理回路252は、第9図に示された再循環型
回路をN+1個、各組に1個ずつ包含しており、従つて
N+1の組がすべて同時に処理されうる。それら種々の
組に対して重みをつけたL係数の積は記憶器253の区
画253,,2532,2533,・・・・・・253
。に記憶される。1つの区画は1つの再循環型L係数重
み付け回路に関連して作用する。
これらの記憶器に記憶されたデータは、正弦法則による
補間の後に、出力マトリクス記憶器254に転送されな
ければならない。この転送は、補間係数の早見表即ちル
ツク・アツプ・テーブル(100k−Uptable)
が記憶されている基準記憶器255に連結されたアドレ
ス・セレクタ256の制御のもとで行なわれる。これら
の係数は、探索ビームを横切る方向の放射線強度の分布
を表わす正弦波状の重み付け曲線上の点に合致している
。ユニツト134は第3図に関連して前に述べたような
表示、印字記憶のいずれか1つまたは2つ以上をなすユ
ニツトである。所要のデータを記憶器254に入れる目
的のために、セレクタ256が記録器254のアドレス
を順次に選択する。
この際、1つのアドレスを選択する毎にデータ源例えば
ビーム通路データ記憶器257を参照して、選択したア
ドレスを記憶器区画2531,2532,2533,2
534,・・・・・・のアドレスに関係づける。記憶器
254のアドレスは、前述のように最終的な映像再構成
の基礎をなすデカルトの網のメツシユ要素の中点の座標
を示している。そのデカルトの網の小部分が第11図に
示されており、そのメツシユにおいては、代表的なメツ
シユの中点0が座標X,,Xqによつて示されている。
記憶器253のアドレスは、第8図および第9図に関す
る説明から理解されるように、平行な組のビームの中心
線を示しており、1つの記・滝器区画は各1つの平行な
組のビームに関係している。たとえば、角度θ2をなす
平行な組の2つのビームB(Ka)およびB(k+1a
)が第11図に示されており、かつ記憶器253にはこ
れら2つのビームに対応するアドレスがある。一般的に
、任意の平行な組のビームの中心線は、アドレス・セレ
クタ256によつて選択された記憶器254における1
つのアドレスに対応する特定のメツシユ要素の中点をい
ずれも通らないが、そのような点の非常に近くを通るこ
とがありうる。たとえば、第11図においては、B(K
a)およびB(k+1a)の中心線の点0からの距離は
dおよび(a−d)である。ビーム通路データ記憶器2
57は、記憶器254の各アドレスに対応して、そのア
ドレスの中点の各側にあり、その中点に最も近い2つの
ビームの表示を記憶している。これら2つのビームは中
点をはさんでいるということができる。中点アドレスま
での各ビーム中心線の距離の表示も記憶されている。こ
の情報は装置の構成によつて決定され、したがつて計算
機に永久的に記憶されうるものである。アドレス・セレ
クタ256は、出力マトリクス記憶器254内の任意の
アドレスを選択すると同時に記憶器254内の選択され
たアドレスをはさんでいるあらゆる組の2つのビームの
表示をビーム通路データ記憶器257から取り出す。た
とえば、記憶器254内のアドレスX,,y,に関して
は、ビーム通路データ記憶器は、角度θ2をなす組のは
さむビームが距離dおよび(a−d)にあるビームB(
Ka)およびB(k+1a)であるという情報を含んで
いる。その場合、アドレス・セレクタ256が記憶器2
53の各区画における2つの適当なアドレスを選択しか
つそのような各一対のアドレスにおける値(大きさ)を
補間ユニツト255に送らせる。記憶器254内の選択
されたアドレスからのビーム中心線の距離についての情
報が同時に補間器255に与えられる。記憶器253の
各区画からそのように選択された2つの値(大きさ)は
、257から受け取られた距離表示に基いて補間器25
5で選択された適当な正弦波的重み付け係数によつてそ
れぞれ正弦波的に重みをつけられる。これらの重みづけ
は、あらゆる場合に、2つの選択された値(大きさ)の
それぞれによつてなされる寄与を、今問題にしているメ
ツシユ要素に割当てる作用を有する。記憶器253の各
区画からの選択されたデータは、補間器255による重
みづけの後に、記憶器254内の適当なアドレスX,,
yqに転送され、そしてそこで加算される。補間器25
5は、ビーム中心線に関して等しく分散されかつ正弦補
間関数に従うたとえば20の重み付け係数を含むルツク
・アツプ・テーブルを包含しており、その補間で用いら
れる重み付け係数はこのテーブルから得られるものであ
る。実際には第9図に示されたすべての構成要素はプロ
グラムされたデジタル計算機によつて実現される。第1
0図の構成によつて得られるものよりも精度の悪い映像
構成でも差支えないのであれば、補間器255は省略し
てもよい。これに示すために、記憶器253における各
区画のアドレスAOツAl7A2ツA37A4?゜″−
0゜ツAωに記憶された値(大きさ)について前述した
方法で考えてみる。
アドレス、AOは、今考えている平行な組の最初のビー
ム(リード23によつてインターセプトされるビームは
含まれていない)に・対応するものと考えられ、かつア
ドレスA。に最終的に記憶される値(大きさ)は、この
ビームの中心線によつてインターセプトされるデカルト
の網のメツシユに対応する網記憶器(MeshwOrk
stOre)254のすべてのアドレスに転送される。
アドレスA1はその平行な組の2番目のビームに対応す
るものと考えられ、かつそのアドレスA,に記憶された
値(大きさ)は2番目のビームの中心線によつてインタ
ーセプトされるメツシユに対応する記憶器254のすべ
てのアドレスに転ト送される。この転送はアドレスA。
,Al,A2,A3,A4・チ・・・・・・・AO)に
記憶されたすべての値(大きさ)が使用されるまで継続
する。同様な動作が、記憶器254における加算によつ
て適当な最終映像再構成を生ずるために、記憶器253
の他のす−ベての区画についても行なわれる。もし所望
されれば、記憶器254への転送は、検査される横断面
の限定された領域内だけの場所に対応する記憶器254
の限定された数のアドレスに限定されうる。本発明は多
くの組の平行なビーム通路と交差する探索ビームの受け
る対数的吸収に関係するデータを得るためのシステムの
物理的なパラメータに正しく関係づけられた方法で映像
再構成を実現するための装置に主として関するものであ
る。
本発.明の装置はビームの組のパラメータを正しく考慮
したものであり、かつ近似がなされる点においても、こ
の種のシステムにおいて妥当なものである。しかしなが
ら、ビーム吸収データを得るための第1図および第2図
に示された装置は例示のためのものにすぎず、本発明に
よる映像再構成のための装置は多数の組の平行なビーム
通路から吸収データを得るための任意のシステムに適用
できる。前述のように、図示されている映像再構成のた
めに用いられる装置は、第2図に示された構成要素のう
ちの幾つかと第8図に示された構成要素とを除けば適当
なプログラムを与えられたデジタル計算機例えば汎用の
計算機によつてほとんど実現される。しかしながら、そ
の計算機は多かれ少なかれ、特別な集積回路、印刷回路
、永久記憶器等によつて構成されうるものである。以上
の説明および特許請求の範囲において平面状の仮想輪切
り体における吸収の変化に関する表示を発生または構成
するという場合には、表示される量は明確に吸収である
必要はなく、その表示は吸収の或る関数(放射線に対す
る透過のような)であつてもよくあるいは吸収の差であ
つてもよいことが理解されるであろう。
【図面の簡単な説明】
第1図は本発明の1つの実施例によるX線により物体を
検査するための装置の走査機構を概略的に示す図、第2
図はその走査機構から得られた信号を処理してそれらの
信号を、多数の組の平行な通路を通過する探索ビームの
受ける対数吸収を表示する信号に変換するための回路の
概略図、第3図は第2図の一部分を示す詳細なプロツク
図、第4図は第2図の回路からの信号出力を利用して映
像を再構成する回路の概略図、第5図および第6図は第
3図の回路の動作態様を説明するための図、第7図は第
3図の回路で利用される係数の表、第8図は第3図の一
部の詳細なプロツク図、第9図は第7図の回路に用いら
れているデータ処理器の一例を示す図、第10図は第3
図の他の部分の詳細なプロツク図、第11図は第10図
の動作を説明するための図である。 3・・・・・・開孔、11・・・・・・放射線源、12
・・・・・・X線検知器、45・・・・・・対数変換器
、50・・・・・・信号処理器、252・・・・・・処
理装置、253・・・・・・処理済データ記憶器(合計
データ記憶器)、254・・・・・・出力マトリクス記
憶器、256・・・・・・アドレス・セレクタ、257
・・・・・・ビーム通路データ記憶器(データ源)。

Claims (1)

  1. 【特許請求の範囲】 1 被検査体の平面状スライス部分の吸収量を表わす信
    号を発生する放射検査装置であつて、以下の特徴を有す
    るもの。 〔A〕上記検査装置は、 放射線源と、 少なくとも1つのビーム通路に沿つて、 上記放射線源から投射された放射線を、検査の間被検査
    体を配置することにより当該被検査体の平面状スライス
    部分に放射線を受けさせるために設けられた開孔を、横
    切らせた後に検出する検知器と、上記スライス部分に異
    なる角度から放射線を照射するように走査運動を生じさ
    せる走査手段と、上記走査運動の間上記検知器によつて
    発生された信号を組合わせて各角度位置に対応する出力
    信号の組を形成する制御手段とを備え、 上記出力信号は上記スライス部分を走査する異なる方向
    のビームの吸収量を表わし、1つの組の出力信号の数が
    50以上であり、〔B〕上記検査装置はさらに、 出力信号の組を対数出力信号の組に変換する対数変換器
    45と、上記スライス部分に重ね合わされた2次元的な
    観念上のマトリクスを構成する基本領域に対応する記憶
    要素を有する出力マトリクス記憶器254と、処理装置
    252を含み、各対数出力信号から他の対数出力信号に
    応じて修正された修正信号を形成する計算手段と、修正
    信号を、修正信号に対応するビーム通路に属するマトリ
    クスの基本領域に対応する出力マトリクス記憶器254
    の記憶要素に分配するためのアドレス・セレクタ256
    とを備え、 これらにより、出力マトリクス記憶器254に補正され
    たレイヤグラムを形成する。 〔C〕(a)上記計算手段は、上記処理装置252に内
    蔵された係数信号回路を備え、この係数信号回路は、項
    数が1組の対数出力信号の数に対応する一連の区域重み
    信号(L_0、L_1…L_k)を与え、上記一連の区
    域重み信号の各項は、レイヤグラムのための寄与を与え
    、上記区域重み信号は、上記基本領域を順次取り囲む環
    状領域を、当該基本領域からの距離に対応して分類され
    たビームで切取ることによつて形成した円弧の長さを考
    慮して計算されたものであり、上記寄与はマトリクスの
    対応する基本領域からのビーム通路の距離の増大に伴つ
    て減少する傾向を有する。 (b)上記計算手段は、上記合計出力マトリクス記憶器
    254に加えて、上記処理装置252と上記アドレス・
    セレクタ256との間に接続された合計データ記憶器2
    53を備え、この合計データ記憶器253は複数組の上
    記対数出力信号を記憶するのに十分な数の記憶領域を有
    する。 (c)上記処理装置252は、各対数出力信号と区域重
    み信号とを乗算するように構成されている。 (d)上記処理装置252と上記合計データ記憶器25
    3とは、各組の対数出力信号に対して上記合計データ記
    憶器253の一組の記憶領域が選択され、上記(c)項
    における乗算によつて発生された積を、選択した組の記
    憶領域に分配し、これにより各記憶領域に特定の対数出
    力信号と区域重み信号列の第1項(L_0)とを乗算し
    たものと、 上記特定の対数出力信号に対応するビーム通路に1番近
    いビーム通路に対応する対数出力信号と区域重み信号列
    の第2項(L_1)とを乗算したものと、上記特定の対
    数出力信号に対応するビーム通路に2番目に近いビーム
    通路に対応する対数出力信号と区域重み信号列の第3項
    (L_2)とを乗算したもの等々 の和を形成することによつて、当該組の各対応出力信号
    に対応して修正信号が形成されるように結合されている
    。 (e)上記係数信号回路を含む上記処理装置252およ
    び上記合計データ記憶器253は、対数出力信号を個々
    に順次に選択し、上記(d)項の動作を各対数出力信号
    組について繰返し、これにより対応する修正信号の組を
    発生して、上記合計データ記憶器253の1組の記憶領
    域に記憶するように構成されている。 (f)上記計算手段は、マトリクス内に各ビーム通路の
    位置を規定するデータ源257を含み、上記アドレス・
    セレクタ256は、上記データ源257に接続されて、
    対数出力信号の組から形成された修正信号が、上記合計
    データ記憶器253の各記憶領域から、出力マトリクス
    記憶器254の修正信号に対応するビーム通路に属する
    すべての基本領域に対応する記憶要素に分配されるよう
    にする。
JP49047032A 1973-04-25 1974-04-25 透過性放射線に対する被検査体における平面状仮想輪切り部分の吸収の変化に関する表示を構成する方法および装置 Expired JPS5911150B2 (ja)

Applications Claiming Priority (10)

Application Number Priority Date Filing Date Title
GB1952873A GB1471531A (en) 1973-04-25 1973-04-25 Radiography
GB1952873 1973-04-25
GB3881773 1973-08-11
GB3881773 1973-08-11
GB3907073 1973-08-17
GB3907073 1973-08-17
GB3942073 1973-08-21
GB3942073 1973-08-21
GB4750773 1973-10-11
GB4750773 1973-10-11

Publications (2)

Publication Number Publication Date
JPS5028385A JPS5028385A (ja) 1975-03-22
JPS5911150B2 true JPS5911150B2 (ja) 1984-03-13

Family

ID=27516210

Family Applications (3)

Application Number Title Priority Date Filing Date
JP49047032A Expired JPS5911150B2 (ja) 1973-04-25 1974-04-25 透過性放射線に対する被検査体における平面状仮想輪切り部分の吸収の変化に関する表示を構成する方法および装置
JP14275778A Pending JPS54119892A (en) 1973-04-25 1978-11-18 Method of inspecting radiant ray
JP53142756A Expired JPS5926302B2 (ja) 1973-04-25 1978-11-18 放射線検査装置

Family Applications After (2)

Application Number Title Priority Date Filing Date
JP14275778A Pending JPS54119892A (en) 1973-04-25 1978-11-18 Method of inspecting radiant ray
JP53142756A Expired JPS5926302B2 (ja) 1973-04-25 1978-11-18 放射線検査装置

Country Status (11)

Country Link
US (1) US3924129A (ja)
JP (3) JPS5911150B2 (ja)
CA (1) CA994011A (ja)
DE (1) DE2420500B2 (ja)
FR (1) FR2227582B1 (ja)
GB (1) GB1471531A (ja)
HK (1) HK57677A (ja)
IT (1) IT1010063B (ja)
MY (1) MY7800030A (ja)
NL (1) NL175027C (ja)
SE (1) SE399480B (ja)

Families Citing this family (97)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4076985A (en) * 1973-08-18 1978-02-28 Emi Limited Computerized tomographic scanner with beam distribution control
US4103169A (en) * 1973-08-31 1978-07-25 Emi Limited Apparatus for examining bodies by means of penetrating radiation
US4206359A (en) * 1974-01-31 1980-06-03 E M I Limited Radiography
GB1493243A (en) * 1974-02-15 1977-11-30 Emi Ltd Radiographic apparatus
US4091287A (en) * 1974-05-08 1978-05-23 Emi Limited Scanning radiology with initial scan for adjusting system so that detector means operates within its preferred range
US4379329A (en) * 1974-07-20 1983-04-05 Emi Limited Radiology
US4398251A (en) * 1974-12-19 1983-08-09 Emi Limited Radiography
DE2503789C3 (de) * 1975-01-30 1980-10-09 Philips Patentverwaltung Gmbh, 2000 Hamburg Einrichtung zum Ermitteln der Absorption einer Strahlung in einer Ebene eines Körpers, mit einer Anordnung einer Strahlenquelle und einer Vielzahl von Strahlendetektoren, die relativ zu dem Körper kontinuierlich gedreht wird
DE2611706C2 (de) * 1975-02-21 1984-09-27 Emi Ltd., Hayes Computer-Tomograph
DE2511231A1 (de) * 1975-03-14 1976-09-23 Philips Patentverwaltung Anordnung zur ermittlung der absorption einer strahlung in einer ebene eines koerpers
DE2521796C3 (de) * 1975-05-16 1979-12-13 Philips Patentverwaltung Gmbh, 2000 Hamburg Anordnung zur Ermittlung der räumlichen Verteilung der Absorption oder der Emission von Strahlung in einer Ebene eines Körpers
GB1551736A (en) * 1975-06-11 1979-08-30 Emi Ltd Scanning radiographic apparatus
GB1547566A (en) * 1975-07-04 1979-06-20 Emi Ltd Information display arrangements
DE2537333C3 (de) * 1975-08-21 1982-10-28 Siemens AG, 1000 Berlin und 8000 München Verfahren zum Herstellen eines Körperschnittbildes aus der Absorptionsverteilung ionisierender Strahlen nach dem Faltungsverfahren und Einrichtung zur Durchführung des Verfahrens
JPS5242085A (en) * 1975-09-30 1977-04-01 Hitachi Medical Corp Inspecting unit of cross section
US4137775A (en) * 1975-10-02 1979-02-06 Emi Limited Ultrasonic apparatus
US4206653A (en) * 1975-10-02 1980-06-10 E M I Limited Ultrasonic apparatus
GB1558062A (en) * 1975-10-25 1979-12-19 Emi Ltd Radiology
GB1562196A (en) * 1975-11-05 1980-03-05 Emi Ltd Radiography
IT1069995B (it) * 1975-11-25 1985-03-25 Philips Med Syst Inc Procedimento ed apparecchio di tomografia computerizzata
US4042811A (en) * 1975-11-28 1977-08-16 Picker Corporation Tomography system having an ultrahigh-speed processing unit
US4149248A (en) * 1975-12-23 1979-04-10 Varian Associates, Inc. Apparatus and method for reconstructing data
US4149247A (en) * 1975-12-23 1979-04-10 Varian Associates, Inc. Tomographic apparatus and method for reconstructing planar slices from non-absorbed and non-scattered radiation
GB1571800A (en) * 1976-01-15 1980-07-16 Emi Ltd Radiography
GB1569413A (en) * 1976-02-05 1980-06-18 Emi Ltd Radography
GB1572599A (en) * 1976-02-10 1980-07-30 Emi Ltd Radiography
GB1571510A (en) * 1976-02-25 1980-07-16 Emi Ltd Radiography
GB1572600A (en) * 1976-03-04 1980-07-30 Emi Ltd Radiology
GB1578754A (en) * 1976-03-25 1980-11-12 Emi Ltd Patient positioning and supporting arrangement
GB1575411A (en) * 1976-04-15 1980-09-24 Emi Ltd Radiography
GB1577615A (en) * 1976-05-13 1980-10-29 Emi Ltd Radiography
FR2403059A1 (fr) * 1976-05-18 1979-04-13 Emi Ltd Appareil de tomographie pour l'examen du corps ayant des organes en mouvement
GB1536448A (en) * 1976-06-01 1978-12-20 Emi Ltd Radiography
GB1578442A (en) * 1976-07-01 1980-11-05 Emi Ltd Radiography
US4149079A (en) * 1976-07-14 1979-04-10 Elscint, Ltd. Method of and means for scanning a body to enable a cross-section thereof to be reconstructed
GB1584954A (en) * 1976-08-17 1981-02-18 Emi Ltd Radiography
GB1589469A (en) * 1976-09-23 1981-05-13 Emi Ltd Radiography
US4326252A (en) * 1976-11-29 1982-04-20 Hitachi Medical Corporation Method of reconstructing cross-section image
CA1111575A (en) * 1977-01-19 1981-10-27 General Electric Company Apparatus and method for reconstructing data
GB1594751A (en) * 1977-01-31 1981-08-05 Tokyo Shibaura Electric Co Method and apparatus for tomography by means of penetrating radiation
GB1595805A (en) * 1977-02-24 1981-08-19 Emi Ltd Radiography
JPS53126891A (en) * 1977-04-12 1978-11-06 Toshiba Corp Computer tomographic diagnosis apparatus
GB1598685A (en) * 1977-04-28 1981-09-23 Emi Ltd Radiography
GB1576286A (en) * 1977-07-12 1980-10-08 Emi Ltd Radiograhy
DE2834306C2 (de) * 1977-08-05 1987-04-02 EMI Ltd., Hayes, Middlesex Einrichtung zum Nachweis von Röntgen- oder Gamma-Strahlung
US4135247A (en) * 1977-08-15 1979-01-16 Siemens Aktiengesellschaft Tomography signal processing system
CA1104727A (en) * 1977-08-18 1981-07-07 Godfrey N. Hounsfield Rotation-only ct scanner with beam deflection
US4293912A (en) * 1977-09-30 1981-10-06 Technicare Corporation Tomographic apparatus
US4670892A (en) * 1977-11-15 1987-06-02 Philips Medical Systems, Inc. Method and apparatus for computed tomography of portions of a body plane
US4339799A (en) * 1977-11-15 1982-07-13 New York University Tomographic scanner
US4181939A (en) * 1977-12-30 1980-01-01 Union Carbide Corporation Scanner data multiplexer for interfacing a radiation detector array and a computer
US4437006A (en) 1978-02-21 1984-03-13 Scientific Measurement Systems, Inc. Method and apparatus for measuring radiation in computer-assisted tomography and radiographic applications
US4263916A (en) * 1978-03-27 1981-04-28 University Of Southern California Image averaging for angiography by registration and combination of serial images
US4217641A (en) * 1978-04-28 1980-08-12 U.S. Philips Corporation Correction for polychromatic X-ray distortion in CT images
GB1602072A (en) * 1978-05-12 1981-11-04 Emi Ltd Processing arrangements for constructing representations of parts of bodies
GB1602071A (en) * 1978-05-12 1981-11-04 Emi Ltd Processing arrangements for medical examination
GB1584949A (en) * 1978-05-25 1981-02-18 Emi Ltd Imaging systems
GB1584950A (en) * 1978-05-25 1981-02-18 Emi Ltd Imaging systems
GB1578910A (en) * 1978-05-25 1980-11-12 Emi Ltd Imaging systems
US4204123A (en) * 1978-06-20 1980-05-20 Union Carbide Corporation Radionuclide body function imager
GB2027262A (en) * 1978-08-04 1980-02-13 Emi Ltd Onisation chambers
JPS5532534A (en) * 1978-08-31 1980-03-07 Hitachi Medical Corp Radiation detector in ct
EP0010885B1 (en) * 1978-10-24 1983-02-23 EMI Limited Computerized tomographic apparatus
US4222104A (en) * 1978-11-02 1980-09-09 E M I Limited Radiography
US4340862A (en) * 1978-12-13 1982-07-20 Picker International Limited Imaging systems
US4339716A (en) 1979-05-23 1982-07-13 Picker International Limited Nuclear magnetic resonance systems
US4250387A (en) * 1979-10-12 1981-02-10 Emi Limited Medical radiographic apparatus and method
US4316091A (en) * 1979-11-15 1982-02-16 Emi Limited CT Scanner
US4288264A (en) * 1979-11-21 1981-09-08 Emi Limited Detector construction
NL7908545A (nl) * 1979-11-23 1981-06-16 Philips Nv Inrichting voor het bepalen van een stralingsabsorptie- verdeling in een vlak van een lichaam.
GB2076541B (en) 1980-05-21 1984-05-31 Emi Ltd Nmr imaging apparatus and methods
GB2076540B (en) 1980-05-21 1984-02-01 Emi Ltd Sampling arrangement for nmr imaging apparatus
FR2486269B1 (fr) * 1980-07-04 1986-03-28 Thomson Csf Systeme de traitement et de memorisation de donnees numeriques, notamment pour tomodensitometre, et tomodensitometre comportant un tel systeme
US4428059A (en) 1980-07-30 1984-01-24 Honeywell Inc. Real time fill circuit
CA1159957A (en) * 1980-07-30 1984-01-03 James M. Gessert Real time fill circuit
FR2496934A1 (fr) * 1980-12-18 1982-06-25 Thomson Csf Appareil de tomographie axiale transverse a calculateur analogique
JPS57134142A (en) * 1981-02-12 1982-08-19 Tokyo Shibaura Electric Co Computed tomography apparatus
US4504438A (en) * 1981-12-07 1985-03-12 Levy Richard H Method and apparatus for determining the density characteristics of underground earth formations
US4608635A (en) * 1982-08-03 1986-08-26 Thomas Jefferson University Method and apparatus for tomographic diagnosis
US4532591A (en) * 1982-08-03 1985-07-30 Thomas Jefferson University Method and apparatus for tomographic diagnosis
US4559639A (en) * 1982-11-22 1985-12-17 General Electric Company X-Ray detector with compensation for height-dependent sensitivity and method of using same
JPS59190643A (ja) * 1983-04-14 1984-10-29 Hitachi Ltd 核磁気共鳴を用いた検査装置
JPS6085357A (ja) * 1983-10-17 1985-05-14 Toshiba Corp 核磁気共鳴映像装置
US4763008A (en) * 1983-12-27 1988-08-09 General Electric Company Ionization detector with conductive signal and ground traces
FR2562371A1 (fr) * 1984-03-30 1985-10-04 Thomson Cgr Procede de reconstruction d'une image a haute resolution, par tomodensitometrie
NL8601004A (nl) * 1986-04-21 1987-11-16 Philips Nv Computertomografie-inrichting met gekombineerd recursief- en convolutiefilter.
US5319547A (en) * 1990-08-10 1994-06-07 Vivid Technologies, Inc. Device and method for inspection of baggage and other objects
FR2675927B1 (fr) * 1991-04-25 1993-07-30 Inst Nat Rech Inf Automat Procede et dispositif d'aide a l'inspection d'un corps, notamment pour la tomographie.
US5651047A (en) * 1993-01-25 1997-07-22 Cardiac Mariners, Incorporated Maneuverable and locateable catheters
US5682412A (en) * 1993-04-05 1997-10-28 Cardiac Mariners, Incorporated X-ray source
US5550378A (en) * 1993-04-05 1996-08-27 Cardiac Mariners, Incorporated X-ray detector
US5485528A (en) * 1993-10-15 1996-01-16 Star Technologies, Inc. Apparatus and method for back-projecting an image
RU2120122C1 (ru) * 1996-10-29 1998-10-10 Войсковая часть 75360 Способ томографического контроля
WO2003081220A2 (en) 2002-03-19 2003-10-02 Breakaway Imaging, Llc Computer tomograph with a detector following the movement of a pivotable x-ray source
JP5572386B2 (ja) * 2006-10-18 2014-08-13 コーニンクレッカ フィリップス エヌ ヴェ 関心領域を撮像する撮像システム、撮像方法及びコンピュータプログラム
US9417340B2 (en) * 2012-07-06 2016-08-16 Morpho Detection, Llc Compact geometry CT system
CN118275958A (zh) * 2022-12-29 2024-07-02 通用电气精准医疗有限责任公司 磁共振成像系统及方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1283915A (en) * 1968-08-23 1972-08-02 Emi Ltd A method of and apparatus for examination of a body by radiation such as x or gamma radiation

Also Published As

Publication number Publication date
JPS54119892A (en) 1979-09-18
DE2420500B2 (de) 1979-12-06
SE399480B (sv) 1978-02-13
AU6813274A (en) 1975-10-23
MY7800030A (en) 1978-12-31
NL175027C (nl) 1984-09-17
GB1471531A (en) 1977-04-27
JPS5028385A (ja) 1975-03-22
DE2420500A1 (de) 1974-11-07
NL7405610A (ja) 1974-10-29
FR2227582B1 (ja) 1979-04-13
HK57677A (en) 1977-11-25
IT1010063B (it) 1977-01-10
US3924129A (en) 1975-12-02
JPS54119891A (en) 1979-09-18
NL175027B (nl) 1984-04-16
FR2227582A1 (ja) 1974-11-22
CA994011A (en) 1976-07-27
JPS5926302B2 (ja) 1984-06-26

Similar Documents

Publication Publication Date Title
JPS5926302B2 (ja) 放射線検査装置
US4149249A (en) Apparatus and method for reconstructing data
US3778614A (en) Method and apparatus for measuring x- or {65 -radiation absorption or transmission at plural angles and analyzing the data
US4233662A (en) Radiography
US5375156A (en) Method and apparatus for 3-D computer tomography
US5909477A (en) Computed tomography scanning target detection using non-parallel slices
US5293312A (en) Method and apparatus for computing tomographic scans
EP1677253A1 (en) Method and device of reconstructing an (n+1)-dimensional image function from radon data
US6256365B1 (en) Apparatus and method for reconstruction of images in a computed tomography system using oblique slices
GB1577015A (en) Tomographic apparatus for reconstructing planar slices from non-absorbed and non-scattered radiation
JPH08509408A (ja) 円すい状ビームデータからの画像の再構成
JPS6340538B2 (ja)
JPH0919425A (ja) X線コンピュータ断層撮影装置
JPS61181443A (ja) 回転式扇形ビ−ムct装置の回転中心のずれを補正する方法と装置
US4052619A (en) Method and apparatus for measuring and analyzing radiation transmitted at plural angles
JP2007512034A (ja) 発散ビームスキャナのための画像再構成方法
IL116258A (en) Image reconstruction apparatus and method for helical scanning
US4639941A (en) Radiography
US4002911A (en) Data acquisition in tomography
US4583241A (en) X-ray tomographic apparatus
CN101785029B (zh) 根据能量相关投影数据对感兴趣区域进行成像的成像系统
JPH0998968A (ja) 物体の断層写真像を作成する方法及び装置
US4032761A (en) Tomography
US4114042A (en) Radiography
US4191892A (en) Spread beam overlap method