JPS5980234A - Ｃｔ像を補償する装置 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 この発明は計算機式断層写真法、更に具体的に云えば、
物体が走査器の視野を越えることによってＣＴ像に生ず
るアーティファクト（人為効果）を少なくする方法並び
に装置に据する。

ＣＴ走査器が正確な再生像を発生する為には、検査する
身体が走査器の視野以内に局限されること、並ひに正確
な基準チャンネル・データを利用し得ることが必要であ
る。典型的には、減衰データを測定するチャンネルに隣
接している基準チャンネルは、入射放射ビーム強度の監
視装置として作用する。視野は、身体を露出し、身体の
周シの複数個の角度から放射を検出する為に、源か、或
いは源と検出器がその周りに回転する区域であり、典型
的には円形である。身体が視野を越えて伸びると、若干
の測定される投影は、（イ）基準検出器が遮蔽され、こ
うして若干の図に対する基準チャンネル情報がくずれる
こと、並びに（ロ）身体の内、視野を越えて伸びる部分
に関する減衰情報を収集することが出来なくて、截頭し
た図になる為に、誤シがある。前者の場合、最終的な画
像は不正な基準レベルのシフトによシ、低い周波数の陰
影を持つ１．後者の場合、人為効果が典型的には低い空
間周波数の陰影又は力、ピングとなって現われ、時には
縞となって現われる。

成る１ｌｉｉの走査では、身体が視野の中にとソする様
（（保証しながら、所望のスライス全行ることが簡単に
シよ出来ない。例えば、大きな人の胴体の上側部分を走
査する時、肩が視野の中にとソまる様に保証することが
出来ないことがある。

以上に述べたことから、この発明の目的は、截頭した一
組の投影から作られた画像の人為効果を最小限に抑える
ことである。

この目的を達成する為、この発明の目的は、不良の基準
情報を持つ投影ケ同定し、正規化する前に、この基準情
報全改善することである。

仁の発明の別の目的は、截頭した投影の絹を同定すると
共に、截頭しなかった組から取出した最善の計画情報を
用いていう云う組を完全にすることである。最後に、こ
の発明の目的は、ＣＴ像を発生するのに必要な時間を許
容し離い程延ばさない様な能率の良い形で前述の目的を
達成することである。

この他の目的並びに利点は、以下図面について詳しく説
明する所から明らかになろう。

この発明の成る好ましい実施例を説明するが、この発明
をこう云う実施例Ｖこ制限するつもりはないことをこと
わっておきたい。むしろ以下の説明の意図する所は、特
許請求の範囲の記載によって定められたこの発明の範囲
内に含まれる全ての変更並びに均等物を包括することで
ある。

第１図にＣ’Ｔ走査器の主な素子が概略的に示されてい
る。、例として、「回転専用」又は「第３世代」走査器
が示でれているが、この発明の考えがこの様な形状に制
限されないことに注意されたい。走査器が透過性放射の
源１０を持っており、これは回転陽極形ｘｍ管の形にす
る場合が非常に多い、Ｘ線管１０によって発生された放
射が１１でコリメートされ、薄い扇形の放射ビーム１２
を発生し、これが患者用開口１３を介してＸ線検出器の
配列１４に向けて投射さｆする。患者１５の様な検査し
ようとする身体を扇形Ｘ線ビーム１２の通路内で患者用
開口１３の中に６＞置ぎめし、身体を通過するビームが
、出会った物体の密度に応じて減衰する様にする。その
結果、各々の検出セル１４ａ、１４ｂ等が、そのセルが
受取った放射の強度に依存する電気信号を発生する。従
って、こうして発生される信号は、Ｘ線ビームが通過し
た身体の部分によるＸ線ビームの減衰の目安である。。

動作について説明すると、源及び検出器の配列を患者用
開口の周りに回転する時、患者に対する複数個の角度位
置で、各々のセルからのＸ線の読みを求める。特定の角
度に於ける夫々−絹の読みは、投影又は図と呼ぶ場合が
多い。利用し得る１つの走査器に対する典型的な図は、
５１２個の個別の検出器の読みで構成される。各々の図
に対する読みをディジタル化し、再生計算機１６に送る
。計算機１６ば、利用し得る多数のアルゴリズムの内の
１つを用いて、扇形ビームが横切った断面の画像を発生
することが出来る。この画像はＣＲＴ１７に表示するこ
とが出来るし、或いはその代りに診断医が更に検討する
為のフィルムを作るのに使うことが出来る。

患者用開口１５の中にあって、源及び検出器の回転軸線
を中心とする円２０ば、第１図に示した形状に対する視
野を表わす１．正確な再生を行なおうとする場合、何よ
りも信頼性のある基準チャンネル情報を得る為に、そし
て２番目に截頭しない投影を発生する為に、患者１５は
視野２０内に局限されていなければならない。

図示の実施例では、検出器配列の末端の検出セル１ｄｘ
、１４ｙけ患者の減衰を測定する為に使わず、その代り
にＸ線源の減衰しない強度に特有な基準の読、みを発生
する為に使われる。従って、患者用開口１３を通過する
が、視野２０の外側を通過するＸ線２１．２２によって
表わされる放射が、何の干渉も受けずに、セル１４に、
１４７に入射することが重要である。こう云う読みが、
全体を２３で示した扇形の残りの部分からの読みを正規
化する為に再生計算機１６で使われる。。

第２図に問題を図式的に示しである。視野２０が、視野
の中にある部分１５ａと視野を越える部分１５ｂ、ｊ５
ｃ’！ｉ７持つ被検体１５と共に示されている。１対の
投影を２５．２６に示１〜である。投影２５は基準Ｘ線
２１．２２が視野２０を妨げなしに通過して、関連した
検出セルによって受取られると云う点で、「正常」であ
ることが判る。同様に、身体全体が検出される放射を通
過させる点で、身体１５の全体が走査される。。

然し、図２６は２通りに截頭される投影を示している。

先づ、両方の基準ビーム２１．２２が身体を通過するが
、基準チャンネル１４Ｘ、１４ｙの読みは、この図に対
するＸ線源の減衰しない強度の真の目安ではない。第２
に、身体全体の減衰が図２６では測定されない。これは
部分１５ｂ。

１５Ｑが視野の中になく、それらがＸ線に露出したとし
ても、こういうＸ線は検出されないからである１、従っ
て、図２６に対しては、走査器は部分１５ａだけで構成
逼れた別の身体を測定していた場合の様に作用する、。

上Ｖこ述べた２つの状態は、ＯＴ像を普通の様に単に処
理しｌξ場合、ＧＴ像に重大な影響を持つことがある。

先づ、若干の投影で基準チャンネルがぼけることにより
、この様にぼけた各々の投影は一様な直流シフトを含む
こと（でなる。この直流シフトは、ぼけた投影を用いて
得られた再生に低い周波数の陰影を発生する傾向がある
。この発明では、単に基準チャンネルがぼけたことによ
って生じた直流シフトは、典型的な臨床では７０Ｈ，Ｕ
・程度になることが判ったが、こればかなりの量である
。

この発明を実施する時、基準チャンネルがぼけた截頭図
に対する対数投影の直流シフトは、基準チャンネル情報
がぼけた各々の図を同定し、截頭しなかった投影から取
出した一層よい情報をそれＶＣ，置き換えることにより
、実質的に減少する。

このやり方は、第３図について説明するのが一番判り易
い。第３図は第１図の装置を用いて実施される一連の工
程を示している。。

この手順を工程３０で開始した後、最初に装置は工程５
１で、源及び検出器が患者の周り全成る範囲のＯＶこわ
たって回転する時、各々の角度θに対する投影又は図を
発生するサイクルを進む、１各々の投影に対して約５１
１個の読みを発生すると共に、源及び検出器が身体の周
りを回転する時、約５７６個の投影をとることの出来る
走査器を利用することが出来る１、この各六の投影には
第６図で工、と示した基準セル強度情報が関連している
。

工程３１で全ての投影を発生したら、装置は工程６２に
進む。この工程５２は、全ての基準チャンネル情報の内
の最大強度を探して、各々の投影を調べる。この最大強
度情報が得られたら、装置は最大強度の成る端数である
閾値強度丁Ｔ　を決定する、１次に各々の基準強度を閾
値と比較して、何れかの図の基準チャンネルがぼけた慣
れがあるかどうかの判定を下すことが出来る。３４に示
す様に、最初の図を選択し、３５の所で、その図に対す
る基準チャンネル情報を求める。工程３６で、基準チャ
ンネル情報を閾値強度と比較する。基準チャンネルが閾
値より低いことが判れば、工程６７によシ、後で走査す
る為に、この投影にフラグをつける１、これに反して、
基準チャンネル情報が閾値よシ亮ければ、この図けほや
けていない基準チャンネルによって撮影されたと結論す
る。その場合、工程３８で強度情報を貯蔵し、工程３９
で指数を増数して、後で使う為の典型的な基準情報を発
生する。工程４Ｄを使って、次の図？選択し、この手順
を繰返す１．各々の場合、閾値強度に対する関係に応じ
て、図は工程３７でフラグをつけるか、或いは工程３８
で和に加算される。全ての図を処理した後、工程４１で
典型的な強度を計算する。

図示の実施例では、よい図に対する平均強度であ勺、こ
れは加算した強度の値をこうして加算した数値の数で除
すことによって得られる。次にプログラムは工程４２で
、工程３７でフラグを付けた各々の図又は投影を検索し
、そういう各々の図に対し、工程４３で、基準強度情報
を工程４１で計算した平均情報に置き換える。

全ての図に対してこういうことを行なうと、最良の基準
チャンネル情報が利用可能になり、工程４４を実施して
、普通の様に各々の図を正規化することが出来る。

基準チャンネル情報をこの様に処理すると、約７０　Ｈ
，Ｕ、　（Ｄ　＃４差？１＝ＳＯＨ，Ｕ以上補正するこ
とｒこなシ、従って２石重な利点になることが判った。

残りの誤差の大部分は、前に述べた２番目の形式の収差
、即ち、身体の一部分が一組の投影に寄与しない様な図
が存在することによるものである。３問題のこういう面
は、第４図について説明するのが一番判り易い。

第４図は成る図又は投影に対して４１１定はれた対数減
衰値のグラフである。縦軸に振幅をとり、扇形の中での
検出器の位置を横＠１に沿ってφとして示す。中心の検
出器ハ配列の中心のφ＝０１Ｃ示す。仁のグラフは、検
出器自体がφｌがらφ２寸での扇形の角度だけを受入れ
るという点で、截頭した図を示している。Ｐｍ（φ）は
検出器の読みによって得られた値を表わす。−組の投影
がφ１及びφ２　を越えて夫々φ貫′　及びφ２′まで
拡がっていることが水爆れている。関数ｇｔ（φ）及び
ｇｚ（φ）は截頭した一組の投影を示す。第２図の身体
の部分１５ｂ、１５ｃが視野より突出していることを前
に述べたが、−組の投影のｇＩ（φ）及びｇ２（φテの
部分が失なわれていることが理解されよう。従って、問
題の投影を収集する際に得られた全ての情報がＰｍ（φ
）で表わきれており、と７″Ｌを普通の様にた＼み込め
ば、ｇｔ（φ）及びｇ２ｃφ）情報が無視されているこ
とによシ、この結果前られる画像に人為効果が生ずる、
。

第４図に示す様に投影が截頭になっている時、即ち、Ｐ
ｍ（φ）情報しか利用することが出来ない時、物体の一
意的な再生は不可能で々）ることか知られている３、然
し、妥当な数の投影（例えば少なくとも半分）が截頭さ
れていなければ、「完備関数−」の低い空間周波数の特
徴の妥当な評価をし、更に処理する為にこの最善の評価
の完備投影を使うことが可能である。

この発明を実施する時、截頭したものも、截頭しなかっ
たものも、投影の複数個のモーメント、実施例でば０次
及び１次モーメントを使って、截頭の程度を決定すると
共に、各々の図に対し、モーメントの他Ｖこ、問題の区
内にある情報だけを使って、截頭情報全補正関数に割振
る。

第６Ａ図で、装置が７０で動作を開始し、最初に工程７
１で、全ての角度θを含む一組の投影を前に述べた様に
発生することが判る。工程７２を使って、はやけた基準
チャンネルの影＠を取去る為に、第６図の手順全体を行
なう。次に工程７３を実施して、この−組内にある各々
の投影に対し、０次及び１次モーメント。。及び。ｌを
決定する。次に第５図について、−組の投影のモーメン
トの性格を簡単に説明する１、モーメントを計算する式
は次の通りである、。

こ＼でＱけモーメントであシ、ｋけその次数を示す指数
でろし、残シの変数は前に説明した通シである。従って
、０次°モーメントは、投影Ｐ（φ）の振幅の積分と云
う性質を持ち、投影内にある合計「質量」に関するもの
とみなすことが出来る１、この質量は、物体全体が視野
の中にとソまっていれば、あまシ変化しないから、０次
モーメントは全ての図に対してかなり一定のま＼にと！
まる筈である１、ビームが平行な場合、０次モーメント
は実際に定数であるが、第１図に示す様にビームが扇形
である場合、０次モーメント５５は正弦状の関数になる
ことがあり、−組全体にわたって１サイクル以上にはな
らず、振幅は比較的小さい。１次モーメントを５６に示
してあり、これは０次モーメントによって表わされる質
量の図心の場所に関係するものと考えることが出来る。

０次モーメント七同じく、これば角度と共に変化するが
、−組の投影の３６０°あたり、完全な２周期しか持た
ない。

この発明では、上に述べたモーメント情報を調べて、予
想きれる挙動からの変化を検出し、この様な変化を生じ
た図を截頭したものと同定し、完備した図で考えられる
更によいモーメントを評価し、その情報を利用して截頭
した投影を完成する。第５図はモーメントの予想しない
様な挙動が起る２つの区域を示しである。これらの区域
は５ン及び５８の括弧で示しである１２区域５ンでは、
０次モーメントＱｏ　　の振幅が予想しない形で増力口
し、最初は挙動のよかった１次モーメント。１　が予想
しない起伏を持つことが判る。この為、区域５７内に入
る各々の図は截頭したものと云うことが出来、図を完備
したものにする為に、実際のモーメントと比較すべき予
想モーメントを決定すべきである。この為、０次モーメ
ントでは、多少とも真直ぐな線６０ｔ１完備しプこデー
タの組のモーメントがそうなると思われる値のよシよい
評価として決定し、真直ぐな線分６１を一組の投影が完
全であれば、その区域に対して考えられる更によい１次
モーメントになると評価する。同様な判断基準によシ、
５８の括弧に囲まれた図は截頭したものと同定逼れ、真
直ぐな線６２が０次モーメントになり得るものと評価さ
れ、真直ぐな線６５が、図が完備していれば、この区域
の１次モーメントとして考えられるものど評価される。

真直ぐな線を使うのが便利であるが、希望によっては他
の曲線を用いてもよい。

第６図に戻って、工程７４は区域５７及び５８の同定を
行ない、この様な括弧の区域内におる各々の図にフラグ
を付ける。次に工程７５が、−組の投影が前に述べた様
に完備していた場合、モーメントがどの様になったかを
決定する。この様な決定をすると、装置はフラグを付し
た各々の図に戻シ、前に述べたモーメント情報をその図
の情報と共に用いて、延長関数（第４図のｇ（φ））を
決定することにより、この図を完成する１、これは、利
用し得る情報に低次多項式をはめることによって行なわ
れる。ｇｌ（φ）を表わすのに選んだ多項式の次数が、
表わ心なければならない制約の数を決定する。次にこの
数が表示で計算しなければならないモーメントの数を決
定する１、こ＼で説明した特定の実施例は１次多項式（
即ち直線）を用いており、この為にはＱｏ及びＱ＋だけ
を考慮すればよい。

特定の低次多項式の係数を決定する為、−組体が成る点
で終るという肉体的な条件に類似する。

第４図は延長した投影の組がゼロになる点φｌ′及びφ
２′ヲ示している。φ１′及びφ２′の実際の位置もや
はシ決定しなければならない、。

０次及び１次モーメントを使う関数には別の２つの制約
が加えられる。上に述べたモーメントに対して最善の評
価値を決定したら、截頭した投影の組に対する実際のモ
ーメントが判っていれば、予想モーメントと同形になる
様に測定したモーメントを変えることが条件である。

有用であることが判った最後の制約は、截頭した投影の
組と延長した投影の組の間の境界にわたって、投影並び
にその微係数の連続性があることである。これはφに対
する成る次数のｇ□　の微係数がφｌ及びφ２で、φに
対するＰｌ　の微係数と等しいと云ってもよい。この制
約に要求てれる微係数の次数は、多項式の次数によって
決定される。。

１次多項式、即ち直線を選んだ場合、完備した投影が境
界で連続的であることを保証しさえすればよい。単に截
頭した組の最後の点の振幅Ｐ１又けＰｌ　　を、夫々延
長した組ｇ１又はｇ２の最初の点に割当てることによっ
て、連続性を達成することが出粱る。

第６Ａ図及び第６Ｂ図で、工程７６．７７が第１の截頭
図（第４図のＰｍ（φ）に対応する）を得る様に作用す
る。次に工程７８が前に説明した様に行なわれ、延長し
た投影の組の最初の点の値をこの組の各々の端に於ける
截頭しなかった投影の組の最後の点の値に設定する、３
この動作が行なわれる点が第４図の７９．．８Ｇに示し
である。

次に装置は工程８２に進み、この工程で０次モーメント
の情報を使うことによって、截頭された質量の量を決定
する。更に具体的に云うと、第５図で述べたことを念頭
において、装置が括弧を付した部分５７内の図に対して
作用していることを仮定すれば、装置は評価した０次モ
ーメント６０と截頭データからの対応する測定モーメン
トとの間の差を決定する。

次に工程８！Ｉを実施して、工程７５からの１次モーメ
ントの情報を用いて、延長した投影の組に於ける質量の
合計量の図心を決定する。このことの物理的な音吐は第
４図から理解はれよう。点８０．８１が判っているから
、工程８２が曲線ｇＩ（φ）及びｇｚ（φ）の下にある
合計面積を決定する。。

こういう関数を更に限定する為、工程８３が投影の組全
体の質１代の中心が何処にあるべきかを決定し、この情
報から、曲線ｇ＋（φ）及びｇｚｃφ）の下にある面積
を割振ることが出来る。これは工程８５で、延長した投
影の組に低次多項式をはめ合せることによって行なわれ
、こうして図示の実施例では、ｇＩ及びｇｚを直線であ
るとして、その勾配及び交点を決定することにより、ｇ
Ｉ（φ）及びｇｚ（φ）を限定する。次に工程８６が次
の截頭した投影を選択し、この手順全体を繰返す。最終
結果として、再生アルゴリズムの少なくとも低周波数の
連続性の条件を一番よく充たすと思われる評価データで
、各々の截頭した投影が完成される。

第６Ａ図及び第６Ｂ図について説明した操作工程は数学
的に説明することも出来る。次に、この発明の方法を実
施する様に、再生用計算機をプログラムする助けとして
、この数学的な説明をする１、最初に、延長関数ｇ１（
φ）又はｇｚ（φ）、又は１を１又は２として一般的１
ｃｇｌ（φ）に対して定めるべき直線は、次の様に定義
することが出来る。。

ｇｌ（φ）＝ａ工φ十ｂ　１（２） 関数を限定する為に、係数ａ１及びｂｉミラ定しなけれ
ばならない。投影の連続性からｇｌ（φｉ）＝Ｐｍ（φ
ｉ）三Ｐ　１　　　　　　　　（３）これＶよ数学的に
は、截頭曲線が第４図の点８０゜８１で延長曲線と滑ら
かに接続することを述べている１、この条件から、振幅
がゼロの時のφ１　の値が次の式で足められることか判
る。

φ１’　−ｂ１／ａエ　　　　　　　　　　　（４）こ
れ全念頭におけば、式（３）は次の様に書き直すと〜ゝ
゛ とを出来る。。

ｇ１ｃφ）＝ａｉ（φ−φ、）　＋　Ｐ、　　　　　　
　　　　（５）更に、延長関数（φ１′−φ□）の水平
１里標に関係する関数δ１　を定義することが出来る１
、これは式（４）に述べた交点及び勾配に等し７いとお
くことが出来る。

従って δ、＝φ、′−φ１＝ｐｌ／　ａｌ　　　　　　　（６
）Ｒｋ、：　Ｑ、に−ｆφ２Ｐｍ（φ）φｋｄφφ１ この式は全ての次数のモーメントに一般的に適用される
が、指数ｋを０及び１に等しいとおくことにより、０次
及び１次モーメント？用いる。量Ｒｋ　は截頭した組の
モーメントとその胡が完全であった場合の予想モーメン
トとの間の差である。

量Ｑｋ　がこの後者を表わし、φ１及びφ２の限界の間
のＰｍ　の積分が截頭データの組からの測定モ〜メント
１表わす。この差は、延長関数に対するｇＩ　（φ）及
びｐ２（φ）の積分について述べた様に、延長関数の間
で分割しなければならない量である。。

式（旧は、延長関数の係数を決定する際にモーメント及
び図の情報全史に利用し易い形になる様に操作すること
が出来る。成る操作の後、０次及び１次モーメントの両
方分使うと、式（７）は次の様になる。

ξｉ　＝：　ｉ／ａｉ　　　　　　　　　　　　　　　
　　（］０ξ１　を消去すると ｖｌξ２”　＋　ｖ２ξｔ＋　Ｖ３＝　ＯＱｌ）こ＼で ２次式の解を求める普通の式全使い、式αηをξ２　に
ついて解くと ξ２＝（ｖ２　％呂−４Ｖ、Ｖ、　）／　２Ｖ＋　　　
　θ罎最後に、ξ２　を使ってξ２　、即ち延長関数の
勾配の逆数を解くと、 最初に実際の截頭モーメントと、問題の投影を截頭しな
い場合の予ｒｔ’ｙモーメントと決定すれｌｄ’、単に
モーメント情報及び問題の図のデータを使うことにより
、任意の図をアクセスして完成することが出来ることは
明らかである１、最初に、Ｐｌ及びＰ２　　で表わした
点８ｏ及び８１（第４図）の１辰幅、０次及び１次モー
メントＲ６，Ｒ１及び切断点φｌ及びφ２ＶＣ於けるチ
ャンネルの角度を使って、式θ→に示した係数全引算す
ることが必要である。３次にこういう係数２式０；ヤで
使って、ξ２　分割算し、その後代（１→でξ２　を使
ってξｌ　を言１算する。交点は式（４）　Ｋよって与
えられ、その図の中にある情報だけを使って、関数８１
（φ）及びｇ２（φ）に対する係数が決定きれ、截頭デ
ータの組が完成する。

投影の片側だけが截頭でれる特別な場合が存在し得るこ
とに注意されたい。この場合、勾配を見出すのに、式（
３）と１次モーメントの式（７）だけを使えばよい１．
投影の左側の端が脱落している場合、勾配の逆数は次の
式によって決定される、。

ξ！＝　２　Ｒｏ　／　Ｐ＋２０！ｉ これに対して右側の端では、 ξ２　＝　　２　Ｒｏ　／　Ｐ２２　　　　　　　　　
　　　αＱ従って、２つの延長関数の係数が完全に決定
謬れ、元の截頭した投影と組合さって、再生アルゴリズ
ムの少なくとも低い周波数の連続性の条件ケ充たす様な
、延長した投影が限定される。

この後、この延長関数を使い、それを標準型のた＼み込
み戸波器でたメみ込んで、画像を再生する為の逆向き投
影に使われる原波データの組を発生することが出来る。

然し、多くの今日ＯＣＴ走査器では、た５み込みはフー
リエ領域で行なわれ、最初に速いフーリエ変換を使った
後、乗算及び逆変換をしてから逆向き投影を行なう。

特定のＯＴ走査器では、循環を避ける為に、約５１１個
の要素を持つ非截頭データの組に１０２４点のＦＦ’Ｔ
を使った。従って、上に述べた様に投影のｌ；＋１を延
長することにより、データの絹が拡大して、１０２４点
のＦＦＴで十分でないと想定するのけ考えられないこと
でもない。従って、この延長データの組に標準型のた＼
み込み方法を使うには、少なくとも２０４８点のＦＦＴ
が必要であり、変換を行なうのに必要な時間は大体２倍
になる。

この発明の特定の１面として、データの組の非截頭部分
に対して普通のＦＦＴによるた＼み込みを使い、延長部
分に対して近似的な閉じた形の解を用い、その後２つの
た＼み込み全組合せて、逆向き投影の為の変更投影の組
を発生することによシ、この難点ケ避ける。即ち、第６
Ｂ図で、工程９０が、截頭したものであってもなくても
、各々の測定した投影の組をた＼み込む。次に、工程９
１が各々の截頭した組？選択して、延長のた＼み込みの
為、閉じた形の解を求める１１次に工程９２が、各々の
截頭した投影の糾に対し、閉じた形のた＼み込みを工程
９０のＦＦＴのた＼み込みと組合せ、身体の一部分が視
野を越えたことを考慮した更に正確な一組のＰ波データ
を発生する。次評価しようとする延長部分の挟み込みは
次の様に表わすことが出来る。

ｑｌ（φ）−ｆｌ（φ）　ｃｏｓφ〕餐Ｋ（φ）、　　
φ１〈φ〈φ２（１りこ＼でＫ（φ）は核であり、この
式は視野の中、即ちφｌとφ２の間だけで評価しでえす
ればよい。特定の核が選ばれた時、今の例で（は、ＬＡ
ＫＳ核と知られるものを選んだ時、弐〇っは閉じた形で
評価することが出来る。最初に、これは次の様に近似す
ることが出来る。

ｑｌ（φ）−ｃｏｅφｉ［ｇｌ（φ）簀Ｋ（φ）〕０槌
次に φ−φ２　　　　　　　　（至） ｇ・（φ）−２・（ｉ−−７７） に対するフーリエ変換は 従って、ＬＡＫＳ核に対して、た＼み込みは次の様に表
わすことが出来る。

そのＭ（は次の様になる。。

これは視野の限界内に対するものである１、ｑｌ　　に
対して同じ様な成分計算することが出来、視野の外側に
ある区域のた＼み込みの、視野の内部にある要素に対す
る影響を評価することが出来る。

扇形のビームの場合Ｋ　”ＬＡＫＳ”核を使うことは厳
密には有効ではない１．この発明では、誤差が生じても
、この延長過程で回りするのけ低い９４１１周波数だけ
であるから、無視し得ることが判った。

延長投影の最初の２つのモーメントラ他の図と合う様に
作る上述の手順を用いた実験で、截頭図が原因で起る陰
影の人為効果が目立って減少することが判った。更に、
投影の組の延長のた＼み込みの評価に閉じた形の解孕使
うことにより、処理時間を許容し難い程長くせずに、こ
の結果が得られる。。

【図面の簡単な説明】

第１図けＣＴ走査器の主な素子を示す略図、第２図は物
体が視野を越えることを示す略図、第３図はこの発明に
従って基準チャンネル情報を改善する様子を示すフロー
チャート、第４図は截頭投影の組？示す図、第１図は截
頭投影の組に対する０次及び１次モーメントを示す図、
第６図はこの発明を実施して截頭投影の組を完備する場
合を示スフ０−テヤートである。 主な符号の説明 １４　検出器配列 １５：患者 ２０：視野 ２１．２２：基準チャンネルのＸ線

Claims (1)

1. 【特許請求の範囲】 １）身体を該身体の眉漫の複数個の角度の放射に対して
   露出し、ｅｆＪ記複数個の角度で身体を通過した放射を
   検出して、−組の投影を作り、若干の投影は身体が視野
   を越えることによって影響を受け、各々の投影には基準
   チャンネル情報が関連しておシ、予定のレベルより低い
   基準チャンネル情報を検出して、基準チャンネル情報が
   身体の影響を受けた投影を同定し、影響を受けなかった
   投影から典型的な基準チャンネル・レベルを決定し、前
   記典型的な基準ナヤンネル・レベルを前記同定場れた投
   影の基準チャンネルに割当て、前記基準レベルを割当て
   て各々の投影を正規化し、−組の正却、化てれた投影か
   ら画像を発生する工程から成るＣＴ像を補償する方法。 ２、特許請求の範囲１）に記載した方法に於て、若干の
   投影は截頭してあシ、更に、前記−組の投影の０次及び
   １次モーメントを決定し、該０次モーメントから截頭し
   た投影から脱落している質量の敵を決定し、１次モーメ
   ントから各々の截頭した投影に関係する延長した投影の
   図心を決定し、各々の截頭した投影に対する脱落した質
   量及び図心の情報を用いて前記投影を延長し、各々の投
   影をた＼み込んで一組のた＼み込み投影を発生し、該−
   組のた＼み込み投影から画像を発生する工程から成る方
   法。 ５）身体を該身体の周りの複数個の角度の放射に対して
   露出し、前記複数個の角度で身体を通過した放射を検出
   して一組の投影を作り、若干の投影は身体が視野を越え
   ることによって影響を受け、各々の投影には基準検出器
   情報が関連しており、前記−組の投影の内の最大基準チ
   ャンネル・レベルをみつけ、前記最大レベルよシ予定量
   だけ小感い基準チャンネル−レベルを持つ各々の投影を
   同定して、基準チャンネル情報が身体の影響を受けた投
   影を同定し、同定しなかった投影に対して基準チャンネ
   ル・レベルの典型的な値を決定し、同定式ＩＬ　ｆｃ投
   影の基準チャンネル・レベルを前記典型的な値に置き換
   え、各々の投影を割当てられた基準レベルに対して正規
   化し、−、＋：ｎの変更し／こ投影から画像を発生する
   工程から成るＣＴ像？袖償する方法、。 ４）特許請求の範囲３）　ＶＣ記載した方法Ｖこ於て、
   若干の投影は截頭してあり、史に、前記−組の投影の０
   次及び１次モーメントを決定し、該０次モーメントから
   截頭しブＣ投影から脱落している質量の量を決定し、１
   次モーメントから各々の截頭し／ヒ投影ｒ（関係する延
   長した投影の図心を決定し、各々の截頭した投影に対す
   る脱落した臀部及び図心の情報を用いて前記投影を延長
   し、各々の投影をた＼み込んで一組のた＼み込み投影を
   発生し、該−１ｌｉＦＩのた＼み込み投影から画像ケ発
   生する工程から成る方法。 ５）特許請求の範囲４）に記載した方法に於て、た＼み
   込む工程が、速いフーリエ変換を用いて各々の截頭しな
   い投影並びに各々の截頭した投影の截頭しなかった部分
   のた＼み込みをし、閉じた形のＭを用いて截頭した投影
   の延長部分をた＼み込み、各々の截頭した投影の組ＶＣ
   対し、閉じた形の解をた＼み込みと組合せる工程を含ん
   でいる方法。 ６）身体を該身体の周シの複数個の角度の放射に露出し
   、前記複数個の角１ｙで身体を通過した放射を検出して
   一組の投影を作シ、若干の投影は身体が視野を越えてい
   ることによって截頭され、前記−組の投影の複数個のモ
   ーメントラ決定し、該モーメントから各々の截頭した投
   影から脱落した質量の量並びに各々の截頭した投影の脱
   落した質量の分布を決定し、低次多項式にモーメント及
   び境界条件をはめて各々の截頭した投影を延長し、各々
   の投影を核とた＼み込んで一組のた＼み込み投影を発生
   し、該−組のた＼み込み投影から画像を発生する工程か
   ら成るＣＴ像を補償する方法。 ７）特許請求の範囲６）に記載した方法に於て、複数個
   のモーメントを決定する工程が、０次及び１次モーメン
   トを決定することを含み、脱落した質量の量を決定する
   工程が０次モーメントを使い、脱落している質量の分布
   を決定する工程が１次モーメントヲ使う方法。 ８）特許請求の範囲６）に記載した方法に於て、た＼み
   込む工程が、速いフーリエ変換を用いて各々の截頭しな
   い投影並びに各々の截頭した投影の截頭しなかった部分
   をた＼み込み、閉じた形の解を使って截頭した投影の延
   長部分ケた＼み込み、各々の截頭した投影に対し閉じた
   形の解をた＼み込みと絹合せる工程を含む方法。 ９）４・￥許請求の範；ノ１ゴロ）、７）又け８）のい
   ず゛れか一項に記載した方法に於て、低次多項式ケはめ
   る工程が、投影の截頭しなかった部分と滑らかに合流す
   る点から延長部分を開始し、ｏ？′に及び１次モーメン
   トを用いて前記、開始する点と延長□した投影がゼロに
   なる点との間で低次多項式をはめることを含む方法。 １０）特許請求の範囲９）に記載した方法に於て、低次
   多項式が１次多項式である方法。