JPH01503490A - 2次元の放射線減衰測定から3次元の光画像形成を得る方法および装置 - Google Patents
2次元の放射線減衰測定から3次元の光画像形成を得る方法および装置Info
- Publication number
- JPH01503490A JPH01503490A JP63504482A JP50448288A JPH01503490A JP H01503490 A JPH01503490 A JP H01503490A JP 63504482 A JP63504482 A JP 63504482A JP 50448288 A JP50448288 A JP 50448288A JP H01503490 A JPH01503490 A JP H01503490A
- Authority
- JP
- Japan
- Prior art keywords
- point
- dimensional
- plane
- attenuation
- radiation
- 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
- 230000005855 radiation Effects 0.000 title claims description 100
- 238000005259 measurement Methods 0.000 title claims description 56
- 238000000034 method Methods 0.000 title claims description 47
- 238000012634 optical imaging Methods 0.000 title description 3
- 229910052704 radon Inorganic materials 0.000 claims description 46
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims description 46
- 238000004364 calculation method Methods 0.000 claims description 31
- 230000009471 action Effects 0.000 claims description 18
- 230000003287 optical effect Effects 0.000 claims description 16
- 230000009466 transformation Effects 0.000 claims description 14
- 239000013598 vector Substances 0.000 claims description 13
- 230000033001 locomotion Effects 0.000 claims description 12
- 150000003257 radon Chemical class 0.000 claims description 12
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 230000007246 mechanism Effects 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 230000008707 rearrangement Effects 0.000 claims description 2
- 238000003384 imaging method Methods 0.000 claims 4
- 206010011416 Croup infectious Diseases 0.000 claims 1
- 201000010549 croup Diseases 0.000 claims 1
- 230000004069 differentiation Effects 0.000 claims 1
- 238000000691 measurement method Methods 0.000 claims 1
- 238000001514 detection method Methods 0.000 description 13
- 238000005070 sampling Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 238000013016 damping Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- 238000007689 inspection Methods 0.000 description 5
- 101150105133 RRAD gene Proteins 0.000 description 4
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 239000004020 conductor Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000002457 bidirectional effect Effects 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 102000016917 Complement C1 Human genes 0.000 description 1
- 108010028774 Complement C1 Proteins 0.000 description 1
- 102100024061 Integrator complex subunit 1 Human genes 0.000 description 1
- 101710092857 Integrator complex subunit 1 Proteins 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 238000012952 Resampling Methods 0.000 description 1
- 241001655798 Taku Species 0.000 description 1
- 230000002745 absorbent Effects 0.000 description 1
- 239000002250 absorbent Substances 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000036621 balding Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- ZNNLBTZKUZBEKO-UHFFFAOYSA-N glyburide Chemical compound COC1=CC=C(Cl)C=C1C(=O)NCCC1=CC=C(S(=O)(=O)NC(=O)NC2CCCCC2)C=C1 ZNNLBTZKUZBEKO-UHFFFAOYSA-N 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 239000012780 transparent material Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
- A61B6/4435—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure
- A61B6/4447—Tiltable gantries
-
- 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/027—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Image Processing (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるため要約のデータは記録されません。
Description
【発明の詳細な説明】
発明の名称
2次元の放射線減衰測定り1ら3次元の光画像形成會得る方法および装置
技術分野
本発明は放射線によって物体の3次元の光画像形部を得る方法、ならびにこの方
法の適用Vc1i2!用される装置IC関するものである。
3次元の再構成は物体を通る放り線の一連の2次元減衰測定の処理に;りて行な
われ、2i!H41i!の入射は1IIl記測定によって変更される。
従来技術
2次元の放射l減衰測定からの2次元の元画像形成は公印でわる。その装置は憔
康診断用X線のごとき適宜な放射線源り)らなり、検査されるべき物体は照射源
と紙片または感光性フィルムとの間にtかれ、該紙片まtは感光性フィルムの点
が物体の出口での放射線の強さに応じて発生される像を有し、該僚に関して観察
され几;ントラストは物体の吸収性領域の位置を示す。
しかしながら、このようにして得られ几情報扛幾つかの用途には不適切でらv7
Dhりしたがって物体の3次元の再構成の方法が!!果された。
磁気共囁トモグラフィ法は高価な設備お工ひ検査されるべき物体が置かれる磁界
の極めて広い均一性の使用上要求する。さらに、3次元再構成を具体化するのに
必要とされる測定時間は極めて長い。したがって、これらの欠点はこれらの方法
の利点全制限する。
同様に、2次元再得成または物体の部分の重ね合せにより3次元再構成全確立す
ることが提案されている。X線源のコリメーショ/は物体の1部分を横切りで1
つ久いでセンサからなるラインの像を発生するファン形状放射機の獲得を可能に
する。X線源は異なる角Ifにより同一部分t−照射するように物体のまわりに
回転する。連続測定は記憶されかつコンピュータは部分の整合の各点における減
衰についての局部的寄与の決定を可能にする。X線源およびセンサは次いでずら
され2)−り他の部分上に叉持されこの部分のまわりにX−源およびセンサは前
記の軌道に対して平行な軌道にし九がって動く。
したがって、検量時間は採用される軌道の数に依存する。実際には、残念ながら
、部分内のサンプリングと同じ位コンバク)Kされた軸方向のサンプリングを軒
容することはできない。ま九、物体は局部の不確実性を増大する、2りの部分の
検査間で移動する危険がわる。
さらに他の欠点はX@源のエネルギ効率管減じかつ装置がこれを冷却させるよ1
51IC検査の間中時々停止させられること全要求するかも仰れないコリメーシ
ョンiC関遅づけられる。
ま几、物体のまわりに回転する円錐ビームを使用する方法が提案されておりかつ
それICより物体の多くの2方向画像を供給する幾つかの照射を実施することが
できる。
もしもこれらの画儂が十分な数でおるならば、コンピュータは物体の3次元画*
t−再構成するためにこれらの画偉を分析しかつ結合することができる。これら
の方法は物体の整合の各点における放射減衰の変換ラドンとして矧られているも
のt使用する。1点での作用の変換ラドンは作用が伝えられる範囲の少なくとも
1つの点を通過する各平面上のこの作用の局部値のすべての合計に等しい。実際
に、これは変換ラドンを説明する丸めに限定数の平面および作用を説明するため
IC限定数の点による別個のトポロジー全開らかに満足させる。
物体の背部に置かれた平らな検出器の12インのセンサ、上での放射線減衰の測
定(sLラドン窒間の平面円に収容される光線ビームIC応じて減衰の合計はこ
の平面に対して減衰の変換ラドンの[全付与する。変換ラドンの数値的反転は作
用の限定範囲のすべての点において減衰管付与する。この方法全体が複雑でおり
かつ提案された幾つかの方法は得られている誤った結果またはそれにも拘らず不
正確であるすべての結果音生ずるということが認められることが必要である。
利用し得る文献の中で、ここでは、変換ラドンを使用する方法がその複雑さの九
めに拒絶される、シュリ/ドクエイノVcよる「ツインコー7ビーム投射による
反復3次元再構成」と題する論文にュークリャ・サイエンスに関するXEBB会
報、N5−25巻、1978年10月、第4号第1135〜1143頁)、およ
び変換ラドンを実行する方法全使用する、ミネルボlCよる「コーンビーム投射
データによる回旋状再構匠」と題する論文にュークリヤーサイエンスに関′する
XBEE会報、N8−25巻、1979年4月、第2号第2682〜2684頁
)を参照することができる。
本発明によって結果として示されるように、1選ラドン自体の使用はそれlrも
拘らず数値計算における近似の使用を要求し、そしてさらに従来技術の論文は3
次元画像の良好な品質の再僕成會許容する具体的な装置を提供しない。
発明の開示
本発明はこれらの欠点に打ち勝つことができる。まず本発明は、そのすべてが物
体の前方の単−円錐放射線源および物体の後方の2次元ビーム検出器からなり、
放射線源およびネットワークが物体に関連して種々の入射に沿って可動でらる3
次元元画像形#:装宜に関する。また本発明は物体の点についての放射線減衰の
変換ラドンの第1誘導体、その計算および反転を実行する方法に関する。
許容し得る結果を得る友めに要求されるサンプリングが紀s!される。とくに補
間法を説明する詳細なフローチャートが提案される。
本発明の他の目的は物体のまわりicなされるような測定のための放射線源と装
置の軌道を付与することにaりこれらの軌道は方法と矛盾しない。
不発明は、まず、物体の照#gVcよる3次元光画像形成妄會1C関し、該装置
は物体が置かれる円錐形状徂閲を照射する放射線源%物体を横切った放射機の減
衰を測定する2方向装重を含む検出器、種々の作用下での物体の、−遍の照射の
実施を可能にする機構、ならびに測定チェーン2よび放射線減衰lCiぬして物
体を示す種々の1合点において局部的fl:可与全厚出するように照射の間中2
次装置の清報を分析しかつ処理するコンピュータからfr、t)、前記放射線源
はただ1つでるりかつコンピュータハ計算を行ないかつ作用力;定義される範囲
の少なくとも1つの点ヲー1!A遺する名平面VC関するこの作用のすべての局
部的値として定義される作用の変換ラドンの誘導体の反転を実行するのIC適す
るユニットからなり、そして変換ラドンの誘導体は運動が球状座像系によって定
義されるA常のベクトルの方向において前記平面に対して垂[iC発生するなら
ば前記平面の各々VC関する変化率の合計として定義される。
考え得る実施例によれば、優薄は原点上に一ア心が合さり、かつ放射線源が円錐
11L道に沿りて照射を実施する円形レールからなり、検出器はこの同一軌道上
で移動しかつ原点に関連して対向位置を占める。
より入念軽冥庭例に↓れば、@構は2本の平行な円形レール、原点に対して対向
位ff1iVC配置された放#g緑源および検出器がそれぞれ摺動する2りの部
分からなる。これらの部分は照射が実施されるとき円形軌道全通過する。
さらに、これらは放射線源および検出器が原点から一定距離にあるような方法に
おいて円弧の形で後方に曲げられる。
本発明はま几、焦点を含む円錐放射線源お工びセンサネットワークから形放され
る2次元検出器から構成される装置の使用によt)′41I体を通る放射線の減
衰の2次元測定からの物の3次元光画像形成方法IC関し、物体を示す整合の第
1点の各々において、減衰が物体の変換ラドンに関連する第2整合の点の放射線
の減衰のに換ラドンを示す多a′t−計算することにより放射線から計算され、
作用の変換ラドンFs、該作用が定義される範囲の少なくとも1つの点を通過す
る各平面に関するこの作用のすべての局部的値であるとして定義され、そして変
換ラドンの誘導体は運動が球状座像系によって定義される通常のベクトルの方向
において前記平面に対して垂直に発生するならば前記平面の各々に関する変化皐
の合計として定義され、量は各第2点に関して円錐放射線の焦点を通過する平面
との検出器の1差ICよって得られる少なくども1つの線に沿う放射稼減衰の変
化、焦点を通過する平面に対してほぼ直交する原点および第2点を通過する11
次いでこれらの加算の直線の組合せを加算することIcよって計算され、第1点
の各々における放射線減設は元の距離に関連するこれらの量の導出によって得ら
れ、かつ最後に導出され友量のN@の結合により、間挿が第2点から第1点に通
る九めKさらに行なわれる。
本号法は、放射線源および検出器か物体のまわりの完全な回転から少なくとも2
つの周期を含むほぼ正弦曲−の形で原点から一定の距*において2りの軌道全通
過する場合に野都合に使用されることができ、正弦曲線の振幅は原点と物体の点
との間の距離に等しいかまたはそれより大きく、軌道の点と原点との間の距離は
さらに物体を通過する平面が軌道に、遭遇するように十分になっている。2つの
周期1c関して、この距離がF「を掛は九この振@に等しいかま′fcはそれよ
り大きいならば証明される。
本発明による画像再り成力法は特殊な計算方法を含んでいる。そのデカルトg標
が一般に均一に分布される物体の点の第1の3次元型合の第1点Vl:、関する
パラメータの計算シー関して、この計算はまず、特徴点の第2および第3の3次
元型合を構成する第2および第3点全足義することからなり、篤2点の球状座標
は均一に分布さり、かつ第2点C2と<IC@線に集中する子午線平面に槁し、
円柱座標が均一に分布する第6点は第1点金倉みかり補1に対し、で直交する子
午線平面および平行な平面の両方に興し;次VC第2点に関する情報七得り1り
導出された情報を計算し;次いで第3点に関する中間のff報を導き出すように
同一子午@平面に属する72点のダルーグに関するこの情報全結合し;そして最
後IC第1点に、橘するパラメータを導き出す工うに同一平行平面に属する第3
点のダルーブ!fc関する中間情報を結合重Sことからなる。
本発明は欠して制限しない以下の添付口面を読むことによりより容易に理解され
ることができる。
図面の簡単な説明
第1図は本発明による装置の主要部、ならびに本方伝を説明するのに使用される
表記を示す概略図;第2図は本発明による考え得る第2の装mを示す概略図;
第3図は本発明による考え得る袈tを示す概略図;第4図は本発明1cよる方法
の1工程を示す幾何学的構造を示す概略図;
第5因は本方法の応用に必要とさノする間挿工8を示す幾伺学的構造を示す概略
図;
第6図は第3図の装置1cよる分析を行なう九めの好都合な二重軌道を示す概略
図;
第7因は本発明による装ffi全操縦する装置全示″f′鋲略図;
第8図および第9図は本発明による他の2つの@激を示す概略図;
第10図は使用される方法の70−tヤードである。
発f3At−夾施する九めの好適な形態本発明による@tはその焦点が8形状で
89かつ分析されるべき物11を減衰により横切りそして減衰時の局部を与(ロ
ーカル拳コントリビュージョン〕を再構成することが望まれる発散円錐ビーム七
発出する放射線源1Dからなる(81図参照)。放射線の性質は、その他方の要
素によって行なわれる一方の要素と異なる減衰ICよる前記要素の存在ま九は集
中?:給断することが望まれる物体110要= 11+ 特徴づけることができ
るその条件でX線のごときものにしても良い。
個々の放射@R1は入力点Meiと出力点Msiとの間で物体11t−横切り、
とくに想像される中間点Mを横切りかつ露光されるT:x射線用のペースとして
役立つ。放射線は最後に検出器121C達しかり点itC&かれたスクリーン1
4上の2次元ネットワークのセンサのうちの特定なセンサ1BVC影響を及ぼす
。それはORy:対して直交する原点0を通過したつ第1図にPa@%で示され
る検出平面上でその軌道Atvm標(p、q)によってマークさね、る。この笑
体のない平面はスクリーン14の遠隔およびその最終的な曲線を考慮する必要な
しに再構成のよ17簡囃な説8Aを許容するように導入される。
もL7も放射線の減衰時点Mにおける局部寄与かで(ロ)で表記されるならば、
センサ13にかくして放射壱強度Ms。
I 1 = X O*xp −f(Mi) dMiei
Tx 抑; r L、減衰は物体11の外部で無視し得ると仮定されかう10に
物体11の不存在において測定された放射@Rユの既匂でかつ一定O強疲である
。
IOはまた、2次元ネットワークのセンサ13上の直接フシックスま几は放kg
e源10(D8力において物体11と接触して測定領域の外部に局部化され几7
ラツクスを評価する監視方法によって供給されることができる。
対数変換後、測定はかくしてB:?ら導出されかり入(第10図の工程101〕
を通過する放射線R1に沿りて物体を通る放射線の減衰について行なわれる。
si
r(Mt)aMi = X(8,A)
ei
その減衰が公匂ておる物体に関する測定の間中最初の段階で計算され几すンプリ
ング係数番;センサによって供給されるデータの補正を可能にする。
2次元画像15はスクリーン14上の物体11から得られる。通常のX線撮影法
はこのような画@rc工って形成される。3次元再構成を得る友めに、物体11
を種々の入射により照射しかつ次いで得られた表示をそれらの間で結合するよう
に物体11のまわりの原点0上に中心が合された円Cs上での物体11のまわり
の放射線源10の焦点8の回転運り′t−プリントすることのみ必要である。
検出器12がビームの運動に追随することは明らかに必要である。ここで、また
ビームは円Ce上で動くと仮定される。ま九、他の方法において、例えば原点0
上に中心が合されるが直径が異なる他の円上で動くこともでき単に画像15の大
きさの増加が異なる。円Ce1C沿う軌道は放射線源10および検出器12に−
その直径が採用され友フッキングモードに応じて円C・の直径に等しいかま友は
それと異なっても良い円形レール20に7ツキングすることにより実現されるこ
とができる。
しかしながら、1つの大きな問題が現われる。プなわち、迅速に記憶されるよう
な情報量が全く無視できないので、再構成の問題の数値的解明は複雑なマトリク
ス計算を使用することに頼る必要がある。しかしながら、本発明蛙この呈の方法
に頼る必要がない。
まず、放射柵の減衰t(M)からの変換ラドンは点M(Rf(ロ)〕を通過する
それらの平WJIC関して定義され5点M(Rt(M)) u各平面上の点にお
ける放射線の減衰にりいて局部寄与を説明する作用のすべての合計vi′定義す
る。これらの平面はnが当該平面に対して垂直な単一ペクト、ルでありかりηD
nが原点OK対する平面の距離と関連づけられる代表的測定であるベアリング(
OM・n、n)によりて各々!!#徴づけられる。これらの平面kW認する丸め
に、独特の点Cが使用され、平面上の原点0の1交投影、すなわち、
OC= (OMon) !!
でらる。
次に放射線の減衰で(ロ)の変換ラドンの第1の誘導体は点M(R’r(M))
’を通過するそれらの平面に関して定mされ。
点M(i’r(M))はベクトル二の方向に関連して減衰についての局部寄与を
説明する作用の誘導体の平面のすべての点VC関する合計に各平面P(OMon
、n)に関しての定義により等しい。
これらすべての平面中で、独特の点CMと関連づけられるラドン平面p(aMi
n、n)を省略PMを介して識別しかつ表わすことができる。nは円C・の平面
において測定されるその長さIVcよ!72>り円coの平面に関連して測定さ
れた全緯閤θによりnが円C・の平面に属するときθ=π/2お工び二がこの平
面に対して垂Iであるときθ=十〇または+πであるように物体11に関連づけ
られるマーキングVCおいて定躾される。
変換ラドンif(M)のま九は点Mでのその第1の誘導体x’t(幻の値を決定
するために、とくに、減衰f(財)についての局部寄与ま几は点MI/c関連づ
けられるラドン平面PMのすべての点に関するその誘導体を説明する作用(ファ
ンクション)の合計を計算することが必要である。これはもしも放射線源10自
体の焦点6がラドン平面PMK属するならば、放射線の1部分またはその近似が
その場合にこのラドン平面PMまたはその近似に残らずりきりしたがって平面ま
たはその近似の外部の物体11の点より減衰を受けないかも仰れないので可能で
ある。減衰についての局部寄与を説明する作用のま穴にラドン平面PMに属して
いる物体11のすべての点に関するその誘導体の合計を評価するためにはこの部
分を検討するだけである(大気は事実上吸収性でない。関連の物体が吸収性環境
内に収容されるならば、外部環境Fi妨害要因とみなされる)。
しかしながら、変換ラドンif(M)の計算はたくして単に近似的になされるこ
とができる。出願人線変換ラドンR’f(M)C+第18導体が他方で正確に計
算されることができ、それが3次元画像の改良され九品質を結果として生しるこ
とを示した。
その焦点8が円C@に沿って物体10のまわりに回転する放射線源10の場合に
おいて、第4図はそのラドン平面PMが2つの点、TGおよび、TDにおいて円
C@j切断する点Mを示す。し7t:J:って放射線源10の焦点8は減衰の合
計ま几はラドン平面PM上のその誘導体の評価を可能にするためにこれらの点の
1つに置り為れなければならない。
物体自体の変換2トンR’fcD第1!!導体は減衰についての局部寄与を説明
する作用の通常のベクトルnlC沿う誘導体の物体の少なくとも1つの点Mを通
過するラドン平面PMK関するすべての合計として定義される。ラドン平面IM
はそれが収容する物体11のすべての点かつとくに独特な点CMに共通でbる。
ここでこの点に関連づけられるは変換ラドンit(M)または平(iiiPM上
のその誘導体u’r(M)の値である。
物体の独特な量は物体の少なくとも1つの点Mを通過するすべてのラドン平面P
Mと関連づけられるすべての独特な点CMK関連している。3次元画像の再構成
扛減衰がこれらの千面丁ぺてに関して利用できるとき行なわれることができる。
しかしながら、測定は焦点81Cよりて通過される軌道に遭遇する平面PM#C
関するR’fの値へのアクセスを単に許容する。測定の独特なf#、Fi軌道の
少なくとも1点を通過するラドン平面に関連づけられるすべての!!!!特な点
0MIC関連している。物体の独特の量は測定の独特な量に出来るだけ含まれね
ばならない。
原点0に中心が合されかつ半径Robおよび円形軌道C@を有する球状物体の場
合において、物体の独特な量は0に中心が合されかつ半径Robk有する。同−
球である。測定の独fFな量は第4図に示されかつ焦点B七通過する平面に収容
される円の回転によって得られるトーラス′!′0でありそして円C・の軸線は
原点0および直径SOのレベルにある軸!!に玉受する。
かくして、測定の独411Fな量扛物体の独4!Fな量全体のカバーを可能にす
ることが観察されることがで色る。丁なわち、円C・ではない物体に遭遇する平
面の独特なシャドー領域が残る。この領域に関するR’fの計算は間挿によって
のみ行なわれることができる。
物体の位tVC関係なく、物体を通過するが円C’eK遭遇しない平面に対して
関連づけられるシャドー領域がある。このシャドー領域を満たすために、平らな
円形軌道を放棄しかつ物体の少なくとも1点をA過する平面が軌道に遭遇するよ
うに軌道を選択することが必要である。
この条件は円形軌道Coが円形レール20に沿う放射組10および検出器120
付随運動によって発生されるならば具体的に満足させられることができそして第
2図に関連づけられるならば、レール20は原点Qを通過する軸線に沿りて角度
ξを枢動するための機傳21でありかくして円CalおよびC112によって実
現される2りの位置を取ることができる。ブロッキングはこれら2つの位置のた
めに2!!互な手段によって設けられることができる。これは角度ξがπ/2に
より近い程、その独特な倉がシャドー領域を持たない原点OiC中心が合され友
物体球の最大手径Rot)がエフ増大するとき容易に実現される。π/2の値で
あるξに関して、すなわち2つの垂直軌道Vr−関して、Robの最大量はRe
が円ご・1またはCe2の半径を示すならばHe/ 、/Tである。
第3図はこのような軌道を具体化する九めの他の可能性を示す。ここで、装置は
2つの円形平行レール30および50′からなり、これらは各々2つの直径的に
対向し九支持体32および34と32′訃よび34′によって通過される同一半
径を有する。支持体32および62′は放射線源10がそれに沿ってハンドル5
6によりて摺動する部分31の末端部を構成し、同様に、支持体34および34
′は検出器12がそれに沿ってハンドル354よって摺動する部分33の末端部
を一匡する。部分31および35は原点OK中心が合された円弧で6る。
し友がって、2つの同中心球Vc@する軌道1c関して回転を置き換える几めに
放射線源10および検出器12の2つの回転管結合することができる。結果とし
て、それらを同一であると仮定する。
1つの考え得る例は第6図に示されかつ式・=xnt・cos2φt−有する軌
道′I′Ct−放射@源1放射熱源8が通過する場合である。ここでψはレール
30および30′に沿う回転角度を示し、e位置点がI点0全通過しかクレール
30および30′に対して平行な平面に属するとき・=ot−有するレール30
おLび30’の軸線に対して平行なりリヤランスを示しそしてInsは1つの振
幅である。
その場合に、ラドンRf(ロ)またはその誘導体R’f(M)c1変換状態が1
1f≧Int1および工nt≧Robでゎるならば物体11のすべての独4ii
Fな点eMK関して得られることができることを示すことができる。ここでRl
fは原点Oに中心が合されかり物体11の丁べてを収容する最小球の半径である
。
放射線源10の移動の間中、検出器12の取看点は同時に原点0に4遅して対称
の軌道rc′の点上に動く。
第2図の装置に対する第5図の装置の利点は二重回転および停止時間および同様
に一定の精密な速度に関して少なくとも2倍長い検査時間を付与するときに測定
の1連続回転icおける笑施を許容することでおる。さらに、円形二重軌道は2
つの軌道の両方に遭遇する平面の高い割合があるので測定において顕著な冗長性
を導入する。
第3図の装置に対する変更例は焦点8が円Ce上で動くとき角変ξを変化させる
ようにライフ23によってコンピュータ50(第7図)Kよりて制ai1される
モータ22を有する第2図の点21t−設けることからなる。軌道は3仄元空間
において得られることができる。放射線源Xおよび光度増幅器のごとき2次元検
出器のそれぞれの質量を考慮して、第3図の装置が好ましいと思われる。
第3図はまた、部分33上の検出器12の運動を制御する装置i示す。該層ff
1iはライン44を介してコノピユータ50(第7図)K甜Vcされる電動機4
oからなり、該電動機の軸は部分33のシック餉の内部に噛合するギヤ421C
よって終端される。運動が大きなnll1ft−要求するとき、部分33の他側
は光学センサ41Vcよってマークされる目盛り43からなる。その場合にセッ
サ41は2イン51を介してコンピュータ50へ信号を送りかり電1ilI7機
40の停止はライン441Cよって制御される。・同じt勤FM訃よび光学セン
サ扛また部分51上の放射線源10の運動を1WIJ飯する。他方で、支持坏3
2 、32’および54.54’に27≧くしてまたラックおよび目盛りり為ら
なるレール30 、50’に沿りてmかす几めにライン144および1451C
よってコンピュータ5C1ICそれぞれ接続された電動fiM14[IDよび光
学センサ141を有する同様な装置がbる。
M1図および第2図の装置はすべてまたここでは示されないけれどもこれらの5
E fltIcよって操縦されることができる。
放射線源およびスクリーンの運動は一般には独立しかり同期さぜられることがで
きる。それらはまた単−モータおよび例えばバーまたはより一般的VCは竪固な
機械的構造によりて行なわれる憬械的運結#′cよりて得られることができる。
物体11を検登するために、また他の装置も可能でわまず、第8図に示されるよ
−うIC5物体11のまわりにコンピュータ50#Cよって制御されるモータ6
oの作用により回転する堅固な機械的構造65に工つて接続される放射線源10
および検出器12のジヨイント運動を伴なう装fkk考えることができる。この
構造はまtモータ66の出力軸である1厘コラム67からなりり為り該コラムか
ら2つの対向する半径方向アーム68および69が突出する。第17−ム68は
放射線源10がそれに垂下される第1ポール70#Cよって、かつ第2アーム6
9は検出器12が垂下される第2ボール71によって終端する。前述の実施例に
おけると同様にこの実施例において物体11は放射Sを透過しかりここでは符号
72によって示される支持体上に置かれる。
第9図に示されるように、さらに、放射線源10および検出器12が不動構造7
5に固定されかつフックされそしてこの図ではそれぞれ物体11の頂部および底
部において示される装置に本発明を適用することができる。
物体11はそのスリーブ80が放射線源10によって発出される放射線に対して
ほぼ垂直な軸線のまわりIc構造75円に設けられt軸受77円で回転する7オ
ーク76によって回転させられる。それぞれ符号81および82として参照され
る2つの軸受はフォーク76の2つの分岐部78と79内に設けられかつ平らな
フレーム83の両@に配置された2つの枢動ビンを受容し、フレーム83にかく
してフォーク76の回転細砂に対して垂直な軸線に対して垂illな軸線のまわ
t)iC2りの分岐部78および79との間で回転する。
物体11は例えば圧縮によって、伸縮自在なロッド85の端部において2枚の圧
縮プレート84間に竪固に固定さ九、ロッドの他端はフレーム83と一体になっ
ている。
伸縮自在のロッド85は共通のfs線ヲ有しかつ互いに広がt)かつロッドのま
わりでまたは内部でばね86によって押されかつフレーム83と圧縮グレートと
の間でフレーム83に関連して圧mされる。かくして、如何なる回転にも従がい
り1りこれまで記憶された他の装置によると同じ検f’に行なうことができる。
圧縮プレートおよび多分ロッドは放射線源10pらの5射紗を透過するものと仮
定される。
フォーク76およびフレーム83の回転速!17はモータ(図示せず)Vcよっ
てコンピュータ50によって制御される。
この装置は人体の検斎に適さないが、小さな物体の非破壊検葺に関心がある。し
クーしながら、フレーム83および分岐部78および79’を横切る放射@を回
避するように回転を制御するかまたはこれらの部分を透明な材料において適!設
けることが必要でわる。
変換ラドンの誘導体R’fKよりて画像の再構成方法を説明する丸めに、新友な
表記(第1図ンが定義されることが必要である。放射線源10の焦点8の位置に
関連づけられる検出平面Pa@tが参照され、この平面は軸@08に対して垂直
でかつ原点0を通過する。この平面はマーク(u、v、 oslos)が真亘ぐ
でらるようにデカルトマーク(u、v)t−備えている。ベクトルUは円C@の
平面に対して平行に選択される。この検量平面pastは検出スクリーン14上
のセンサ13の座llpおよびqk定義するのに役立り。D(OMsn、n)
まtはDMとして示される合計直線は検出平面との点Mと関連づけられるラドン
平面P(OMII!!、!りまたはPMの交点でら9、そして0′はDMの点で
おり、厘繍DM上の点oo*xzする投影でるる。厘@DMは、
(n * v 1t O’ 810 ’ 8 )がM接マークを形成するように
単一ベクトルv1によって方向付けられる。これは値Vとvlとの間の角度と呼
ばれる。
筐た、2つの重みを付けられた減衰作用を定義する必要かわる。すなわち、
Y(8,A) = X(S、A) e Re/8AZ(13,A) = X(8
,A) ’ Re /BAここで、Reは焦点Bが放射線源10から動く軌道の
牛後、すなわち距離os”t’6る。これらの作用は@該うドン平面から独立し
ている。
各[@DMIC関してそれらに関連づけられるのは、合計!線についての作用Y
および2の全体をそれぞれ示す作用BY(8,M)および5z(s、M) テロ
る。
数値計算の開始点である。これは
0S二
ここで、βはベクトルnと08との間の角lit示す。
第1式(1)は変換ラドンRfK対する測定Kg連する近似計算である。第2式
(2)は変換ラドンR’fの第1誘導坏に対する測定に関連する正確な関係であ
り、それは放射線源10と物体11との闇の距離VcPA係なく確認される;り
正確な結果および残部の獲得を可能にする。かくしてともに最大に持ち米たされ
ることができる一方式(1)の近似は放射g源10が不利な形状の装置の使用を
付与する物体11から間隔が!かれるときより正iでわる。
式(2)は以下が主張されるならば数学的に等しい。すなわち、
ts+ Y’(S、A) = cos a x/p ・ (8,A) +5in
a Y/q 拳 (8yA)’2つの式により
ここでA(p)およびA (Q)は横座標pおよび縦座標qの合計[@DMの点
入を示す。
好ましくは、式(6)はπを法として0に近いaVc関して選択されそして式(
4)はπを法としてπ/2に近いaに関して選択される。
第10図の70−チャートによって示されるように、放射線源10の付与された
位置に関して、tX(B、A〕およびY(8、A)¥1検検出面Flatの点A
Vcついて計算される(工程101および102)。検出平面Pa@t (Q座
alpおよびqVc関連してY(B 、 A)の誘導体は2つの差動フィルタに
よる作用Y(B 、 A)f)画旋を使用することによりf過作業(工程103
)によV連成される。放射線源1oはかくしてその次の位置(工程104〕に限
り勤かされかつサイクルはすべての獲得が行なわれるまでやV厘丁。
次の工程は物体の独特な量の独特な点U′に関して計算を行なうことからなる。
これらの計算はそれらに関連づけられるラドン平面PM、ならびに検出平面Pd
@%とのラドン平面PMの交差ラドン[線DMt−導く。計算は以下のように開
始する(工程105)。
式(4)の場合において、
式(5)の場合において、
すなわち、ラドン直線DMおよび久いで直線結合(工程106ンに属する合計に
より以下が得られるようにする。
式(番)の場合において、
式(5)の場合において、
その後これらの量はそれらを
で、次いで1 / sin aまたは1/cesαで乗算することにより標準化
される(工程107〕。その結果はR’ r(o)に等しい。
もしも座標系(ppq)が検出器12上のセンサの分析に:#)一定のサンプリ
ングに適合させられないならば、これらの式は結果として通亘な座標系VC置き
換えられる。
t fC座標系(p、q) IC減衰を間挿することにより逆に再計算すること
ができ、かくして測定の再補正ま友は再サンプリングと呼ばれることを行なう。
ラドン平面が幾つ〃1の点、例えば第4図の場合において2つ(JG>よびJD
)の点において軌道に遭遇する場合において、これらの測定に関するノイズに関
連づけられる統計的なエラーを減少するようにこれらの点の各々lc6連づけら
れるR’fの値の丁べてまkは少なくとも一部に関する平均を実施するのが有用
である。一般に、平均はそれにも拘らず2つの優先的な点に関して選択される。
独特な点CMに関連づけられる2トン平面が測定の独特な量が物体の独%な倉内
にシャドー領域を浅丁ならば発生する点において軌道に遭遇しない場合において
、それにも拘らず間挿方法七介して値R’fが割り当てられることが望ましい。
これを行なうために、はぼゼロの間挿を設定することができる。すなわち、測定
の独特な量の外部で、独特な点CMに関連して、Oで心出しされかり測定O独特
な量[よりCMを通過する球の交点に対応する表面である。仄いてこの表面上で
1点C’M (図示せず〕はCMから最小の距離において選択される。測定の独
特な量の足義IC工り、独特な点C’Mとして許容する平面P’M (図示せず
)は1またはより以上の点番ζ2いて軌道に遭遇し、はとんどの時間がそれらに
対して正接として残る。前記方法はC’Mと関連づけられる値R’fの定義上可
能にする。はぼゼロの間挿はこの同一値を点CMに割り当てることからなる。
り為<シて、値R’fに物体の独特な量のすべての点に割り当てられる。既述の
ごとく、もちろん、シャドー領域に関する間挿を回避しり為りしたがって第3図
に示されるような装置上使用することが好ましい。しり島しながら、装置のより
機械的な簡単の丸めに、簡単な円形軌道C・により満足させられることができか
つ関連のシャドー領域に関する@[ki定するために間挿を許容する。
物体の独特なilC関する変換ラドンの第1誘導体R’fを情@ptからf(M
)t−導出することが単VC残る。この変換作業は7N接でありかつ高容量計算
#′c工って実行される。
一定の単一のベクトル二に関して、代表的測定(半径)pは原点の細繊oおよび
方向を示すベクトルnK11lする点σ′について言及される。
変換ラドンR/ tの第1誘導体の理論的変換式は以下のように書かれる。
ψおよびθは当該点(ここではtr/ )の経lおよび全緯度を示しそして9
:CO(r’= 9 nによって定義される。
しかしながら、この明細書に記載されてない表出する別の問題を調べる必要がる
る。とくに、放射線源10が限定数の決められた入射により単に露光を行なうこ
とができ、検出スクリーン14がま一!′c限定数として検出平面上にマークさ
れるセンサ13のネットワークによって測定を行ない、そして物体11が別個I
Cされるかま几は整合されねばならないことは明らり為である。
必然的に、間挿〔インタポレーション〕の問題が生起する。勧められかつまた本
発明1c属する解決は以下に説明される。検葺は正しい再構成が得られることが
できるように表出する別の問題を解決することに関する条件についてなされる。
原点0上に心出しされかり半径Rot)を有する球面に完全に含まれるとみなさ
れる物体11は偶数の平行六面体の表示また(霊表示メンシュッグを別個にさせ
る結果を生じる点MICより記εされることができ、その座標は、!(1) =
((21−1−Nz)/Nz)sRooここで1≦1≦Hz。
y(j )÷((2j−1−Ny)/Ny)拳Rot)ここで1≦j≦’7tz
(k) = ((2に−1−Nz)/Nz)eRobここで1≦に≦Nz1−[
認する。
イメージ系3Dのために必要とされている遮断周波はがあるならば、Hz 、N
yおよびN1は好ましく扛4・Rot)ニジ大きい〃為ま′fcは等しいとして
取られる。
検出スクリーン14はtfc好都合には、そのうえすでに導入されているデカル
ト座標において目rk?付けされる。検出平面pastよでそれらの座標によっ
て定義6れるセンサ13の配置は、
p(a) = ((2a・1−NP)/(NP−2)) ・Rcbここで1≦a
≦Np。
q(b) = [:(2b−1−Nq)/(Nq−2))・RQりここで1≦b
≦Nqを111gする。NpおよびNqは4・ 争ROCIに等しいPまたはそ
れエリ大きい。
変換ラドンの第1誇導体のサンプリングの点Uは原点Oに心出しされる球状整合
によって定aliされる。それらの球面座標は。
半径: P(!l) = ((2m−1−Nn)/(Nn−2))拳Rradこ
こで1≦n≦Nn。
全緯#L二〇(L) = ((21−1)/2N1)π力ここで1<l≦Nl。
経度: F(m) = ((1!1−1)/Na+)−2gここで1≦!II<
−N!!1
″に確認する。ここでRradは物体の独特な蓋を含む0上で心出しされた球の
半径である。半径RookVする物体の球面IC胸してRrad = Ro b
である。
Nnシて関して、偶数は、これが限定数の平面を特徴づけるので%食合において
好ましくは原点0’z組み込まないように取られる。他方において、以下が選択
される。
丁なわち、
2yrmV@Rob≦Nn≦4z*V*RobK関して、N m :2・Nnお
:びNユニ N rs / 2が選択される。
これは許容し得るレベルへの変換ラドンの第1誘導体の変換に関連づけられる人
為構造の減少を可能にする。
減衰測定がそれによりなされる放射@源10の位置鉱軌遭Ce4たはTcと経変
(、)の子午線面に対して直交する子午線面の交点である。
したがって、実際の問題は、放射1源10によって想像されることができるNm
位fxpら、点Uの完全なネットワーク&C関しての変換ラドンa’t (σ)
の第1誘導俸かり次いで点Mの完全なネットワークに関して減衰f(ロ)につい
ての局部寄与を計算することり1らなる。これは−遍の間挿工程によって達成さ
れる。考え得る1つの方法はこの明細書の延長において詳細に示される。
独特な一定量の一部を形成する点Uに関して、第5図は点σの変換ラドンn’t
(σ)の誘導体が焦点8の少なくとも2つの一般に異なる位置81および82
に関して得られることができ、その直径り;原点0および位&81および8υU
りて制限される球が焦点8に関して点σで合流する。しかしながら、実際には、
位置81およびS2は焦点Sの測定位置に対応するが、各々これらの位置の2つ
の間、すなわち811および812と821および82’2との間に配意される
。
独特な測定量の外部の点σに関して、放射線源の位置81およびs2はmFcf
#照されたゼロ順序の間挿を定義するのに使用される特特な測定量の点からマー
クされる。
実際に、変換ラドンの第1誘導体は丁べての点σ′の一部を形成しかつ点σの平
行線に対応する円とのそれぞれ0811.0812,0821および0B22の
直径を育する琢のσの交差の最も近い点でおる点σ11.σ12σ21およびσ
22に関して計Xされる(同一半径(n)および同−全緯度θ(1)をiするす
べての点〕。もしもこの交差が空ならば、交差円の点にこれらの球とOで心出し
され〃為つ点σQ干行−に対応する円の攻も近い口?通過する球との間で取られ
る。点σのそれらの間隔はそれぞれdll、!12,421および422である
。
変換ラドンは点σに+’らまたはその誘4俸から間挿式(7)%式%(22
を介して計算される。この作業はフローチャートの工程108に対応する。
方法を処理するコノピユータがサンプリングの点σに関しての点811,812
,821,822.σ11゜σ12,021お=びσ22、ならびにメモリIC
西己憚される間隔at1.a12.d2tおよびa22に@以って匂っているこ
とは明らかでおる。点σ11.σ12゜σ21およびσ22に関連づけられる合
計ilI線の位1fはまに前身って知られρ)りその点に沿う減衰値x(s 、
AEこのM線lCよって横切られるセンサ13上で測定さ九る値の間挿に工つ
て得られる。合計厘−〇尚該照に関する合irt t hlブ゛るために、コノ
ピユータ50に〃為〈シて各センサ13上ic関連づけられる重与付σ係はを所
Mする。
提案されtサンプリングVCよれば、点σζ数NIXN!IIX N n :=
N n 3で刀為つ211NΩの放射線源10の位置がある。これらの位置の
各々に関して、変惧ラドンまたはその誘導体はNn2の合計y4. 綴D M
ic関して計算される。
そのる合に、平均は2つづつ給金される値VC関して計算される。
式(6)の数fil計xk行なうためIC1全体を切り隨丁ことができplり好
都合でおる。
この計算の5A施を説明するえめに、再配置され九投射の平面の表記框原点ot
−通過して導入されかつその点は一定の経度ψに6る。かくして、これに子午一
平面を含む。
この平面のすべての点8および経度ψを有するすべてのベクトル!llIc点a
Sが関連づけられ、幅縁上の点Bの直交投射は原点Oを有過しかつ方向性ベクト
ルnt有する。Bが再起#Lされ几投射の平面を説明するとき、点CBは経度v
lC関連づけられる子午砂平面に対応する独特な点の平面を説明する。
もしもBが再配置された投射の平面上の点MC1亘又投射であるならば、点CB
はyおよび!!によって定義される平面DMと関連づけられる独%な点CMと同
一でるることが指摘されることができる。
かくして、物体の少なくとも1点Mの直交投射でらるすべでの点BK関して、以
下の童、丁なわち、を関連O子午1平面に−する値R’fη為ら再配置の投射の
各平面に関して計算することにより開始することができる。
第2に、減衰についての局部計算は物体の各点MiC関して計算される。
この方法は物体0A51Wなiが測定の独特な量に完全に含まれるとしても物体
11の厳密なη樗成の獲得を可能にする。
さらに、正′Rな式の使用は物体11’に投射線源10のより近くに愕ち来た丁
ことt可能にする。したがって、装量の空間的要求を減じかつ匡長要因(焦点8
の距離−検出器12/焦点8の距離−物体11)を増大しかりしたがって装置の
窒間的解像度金改菩することができる。
さらIC5これは照射が物体の固足照射角を拡大する範囲に放射線のより良好な
使用を設けることに寄与する。放射線源が角度単位ごとのフォトンの出力におい
て制限されるので、測定の周期の間中物体を通過するのに要求される7オトンの
全体の数に関して、これは検査時間が減じられかつ一時的に解像髪を増大するこ
とができる。ま7′C1同一検査時間を保持することにより、再lll底され九
物体VC関して乱計的な精度が改菩させられる。
再構成ダイアダラムは、例えば放射線#10からのまたは検出器12からの焦点
の距離が異なる、実際に同様lC変化し得ることを是認することにより広範な種
類の軌道に一般化さ九ることかできる。前式において、底長に関連づけられる重
み付は係数にその場合に測定に依存する。
説明されるべきことは、放射線源10および物体11は軌道が第6図によって選
択されf)hり式・: I n t((2)noψ)に一致するならばともに持
ち米几されることができるということである。この場合にnはより高い。
上述し次式によって定義される変換ラドンの誘導体の反転はサンプリング、すな
わち、
r(c) = ((2e−1−Nn)/(Nn−2))Rradここで、1≦C
≦Nn
z(a)= ((2(li−Nz)/1Jz)*Robここで、1≦d≦Nz
(コ)= C(rn−1)/′Nn・2〕ここで、1≦m (h Nm
によって点σの子午線と関連づけられる再配置の平面上にプロットされた再記f
It整合として知られる整合の点8によって具体的に行なわれる。r (c)は
円C・の平面に対して平行な@@ ic沿う子午線平面の点Bの座標を示しかり
Z (a)は円Cmの回転m線に沿う再配置点80座像を示す。各子午一平面上
で、点8は矩形ネットワーク上vc、均一に分布される。
@VC定義されえような量Q(ψ、a)の計算は、したがって、間挿が子午−〇
各方向に関して誘導体aR’f(σ)/apから得られるように誘導体θR’f
(CB)/c?p f可能にしされる(工;!it 1 j )。
誘導体θR’f/clpの計算は好都合1cは、例えば関連のフィルタによる回
旋のようなデジタル処理技術によって行なわれた。データはま几、まず要因si
I+θ[jりで重み付けされる〔重み付けおよびf過工程109および110〕
。
公匂の技術Vc工れば、θR’ j/&pの計算に関連づけられる一過作業およ
び余緯度θに関する合計ice迷づけられる再投射作業は本発明の文脈から逸脱
することなく1りの濾過作業が研いて行なわれる再投射によってtき換えられる
ことができる。f過作業はま友多分角度θに関する合計後物体11上でなされる
ことができり1つ実2I!!iされた。
実施されるべき最後の間挿は子午線平面に対して垂亘な物体の各平面の内部でな
される。点Mの減衰f(ロ)についての局部寄与の計算が考慮された各配置投射
面に関する点MO直交投射を導入することが判った。これらの投射はそれらを加
算する前に責Q(ψe B)k間挿することが必要である点B間ic同き、かく
して再投射作莱を寅現する(工程112)。
これらの最終作業1ctsu、て%装置は1几、異なる整合の点の位置を前極っ
て気付き%かり表またはマトリクスの形で所有しかつ考慮されるべき係数を所有
する。かくして従来のプログラミング作業は重要で8!7刀1つ幾つかの異なる
サンプリング七茗することが望まれるならば繰り返されねばならないが、検量の
間中なされるべき計算は相変らず妥当でらりたり三としてl−結合およびf過か
らなる。
コンピュータ5υは河薄成方法を制御する(第7図〕。
該コンピュータは同期ユニット51、メモリユニット52、計算ユニット53お
よび周辺ユニット54たらなる。最初に、同期ユニット51は測定の獲得を保証
し、光学セ/す41または141の指示VC応じて、X動機40まfcは140
t−作動しかつこれらを放射線源10および検出器12が予め定めた測定位@に
、おるとき停止する。
次いでタイン571Cよって、その処理され九情報がライン58によってメモリ
ユニット52に供給されるセ/す16の各々にライン56T/Cよりて戻続され
る通常の測定回路55での取得の開始および終了全指示する。
取得が放射線源10の一定の位置について終了されるとき、同期ユニット51は
電動!!&40または140會再始動し、一方センサ13はゼロにリセットされ
る。次いで次の測定が行なわれかつ記憶されることができる。
さらに、放射線源10および検出器12を連続的に動かすことができる。測足ヲ
実禿するのに伴われる合計時間はモータがもはや停止されないので、その場合に
明らかに減じられるが、各照射の間中め回転により画像の焦点の不足を受容する
必要がるる。
記憶されるべきことに、同期ユニット51上の照射の位置を示す目盛り43は必
要でなく、それらはまた除去されることかで色そして同期ユニット51自体が該
同期ユニットに組み込まれ友クロックVCよって照射のサイクルを決定するとい
うことである。
すべての測定の取得後、物体の実際の検査が完了さ?する。すべての情報はライ
ン59全介して計算ユニット53に通り、変換ラドンまたはその第1誘導体の計
算ζそれらの反転と同様になされる。
減衰f(M)Kりいての局部寄与の値はとくにグラフ出力および表示スクリー7
を含んでいる周辺ユニット54に向りてライン60により最後に送出される。
ここに記載され次微分および間挿は単に計算手段としてみなされるべきでなくか
つ本発明の前後関係〃〉ら逸脱すること で他の演算VCよって置き換えられて
も良い。
同様に、提案さ−またサンプリングは単なる例で89pつ比較的広い限界内で変
化することができる。
しり)しながら、留意されるべきことは、変換ラドンの誘導体、円形軌導>よび
円錐ビームの表現は、分析的ρ為り数値的稍変の理由のため、平らな部分が定義
されることができる点MO矩形ネットワークによって、擬似球形状によると同様
に、物体を示すことが当り前でおる点りの球状ビームiCxる作用を付与すると
いうことでおる。
幾つかの間挿演頁はしたがって一方のネットワークから他方のネットワークに通
すために必須である。他方において、物体110点りの整合の音饗に数Nmにエ
フで定義される測定の数に厳しく依存しない。すなわち、このレベルにおいてし
7tがって一定の柔軟性が必要である。
最後に、図示によってここに説明された運動と異なる放射線源、物体ンよび検出
器の運動?伴なう装置が本発明の文脈内にま)ることは明らかでおる。同様に、
測定は同時に作用する幾つかの放射線源によってより迅速に行なわれることがで
き、図はそれらの間に分布される。
したがって、不発FIAは放射線の減衰の2次元測定から3次元再構成?得る個
めて有用な方法を提供する。1つのもつともらしい用途は明らかに医用イメージ
でありかつ第2は非破壊制御であるが、装置の寸法および利用し得る放射@に矛
盾のないすべての物体11は検査されることができる。使用されるアルゴリズム
はとくに反復する方法と比較される処理Yr7JrI速する解法に直接計算する
。
最後に、本方法は測定される情報の正確な局部化を保証する。イメージが従来技
術の他の方法を使用することにより作られる場合でないイメージのひずみは現わ
れない。
測定装置および使用される慣例によれば、種々の標準化要因ま几は尺度変化がな
されることができる(70−チャートの工程113〕。したしま九、例えば、同
等の水の長さとして減衰測定を表わし、かつ次いで再構成後医用スキャナXの慣
例IC応じてハウンズフィールドヲ使用することにより計算の結果を香き換える
ことができる。
国際調査報告
mAm5w52.p(7/ji ε8100254国際調査報告
Claims (20)
- 1.物体(11)を通る放射線の2次元減衰測定からの3次元画像形成装置にお いて、前記物体が量かれる焦点(S)からの円錐形状空間を照射する放射線源( 10)、前記物体を横切つた放射線の減衰を測定する2方向装置(18)を含む 検出器(12)、種々の入射により一連の測定(11)の実施を可能にする機構 (20,30,30′)、ならびに測定チェーン(55)および放射線減衰に関 して前記物体(11)を示す種々の整合点(M)においてこれから局部的な寄与 f(M)を導出するように照射の間中2次元装置(13)の情報を分析しかつ処 理するコンピュータ(50)からなり、前記放射線源(10)はただ1つであり かつ前記コンピュータ(50)は放射線の減衰の変換ラドンの誘導体の計算およ び反転を行なうのに適するユニツト(58)を含み、作用の変換ラドンは作用が 定義される範囲の少なくとも1つの点を通過する各平面に関するこの作用のすべ ての局部的値として定義され、そして変換ラドンの誘導体はもしも運動が球面座 標系によつて定義される通常のベクトルの方向において前記平面に対して垂直に 発生するならば前記各平面上の変化率の合計として定義されることを特徴とする 2次元の減衰測定からの3次元面像形成装置。
- 2.前記放射線源(10)および検出器(12)は種々の入射に関して固定の原 点(O)から一定の距離に残ることを特徴とする請求の範囲第1項に記載の2次 元の減衰測定からの3次元面像形成装置。
- 3.前記機構が原点(O)上に中心が合された円形レール(20)からなり、該 レールは前記放射線源(10)と前記検出器(12)が一定の距離で原点Oと常 に整列されているように測定の間中前記放射線源(10)と前記検出器(12) を案内し、前記放射線源が円形軌道Ceを描くことを特徴とする請求の範囲第2 項に記載の2次元の減衰測定からの3次元面像形成装置。
- 4.前記円形レール(20)が前記放射線源(10)を一定角度にずらされた2 つの円形軌道(Ce1,Ce2)に応じて動かすことができる枢軸(21)を備 えていることを特徴とする請求の範囲第3項に記載の2次元の減衰測定からの3 次元画像形庇装置。
- 5.前記円形レール(20)はモータ(22)を時間制御の可変回転に従わせる ことができる前記コンピュータ(50)によつて制御される前記モータによつて 制御される枢軸(21)を備えていることを特徴とする請求の範囲第3項に記載 の2次元の減衰測定からの3次元画像形成装置。
- 6.前記機構は2本の平行円形レール(30,30′)、前記放射線源(10) および前記検出器(12)がそれぞれその上で枢動する2つの部分(31,33 )からなり、これらの部分(31、33)は原点(O)に関連して対向位置に配 置されかつ測定が行なわれるとき前記円形レール(30,30′)を横切り、前 記部分(31,33)はさらにそれらの点が原点(O)から一定の距離にあるよ うに円弧として後方に屈曲され、前記放射線源および前記検出器は永久的に整列 されかつ原点(O)から一定距離に残ることを特徴とする請求の範囲第2項の2 次元の減衰測定からの3次元面像形成装置。
- 7.前記放射線源(10)および前記検出器(12)はマーキング装置(41お よび43)により前記コンピユータ(50)により操縦される電動機(40)に よつて前記レールまたは部分(31、33)上で動くことを特徴とする請求の範 囲第3項ないし第6項のいずれか1項に記載の2次元の減衰測定からの3次元面 像形成装置。
- 8.前記部分(31,33)は前記マーキング装置(41,43)により前記コ ンピュータ(50)によつて操縦される電動機によつて前記円形レール(30, 30′)上で動くことを特徴とする請求の範囲第6項に記載の2次元の減衰測定 からの3次元面像形成装置。
- 9.前記放射線源(10)および前記検出器(12)は堅固な機械的構造(65 ,75)によつて接続されることを特徴とする請求の範囲第2項に記載の2次元 の減衰測定からの3次元面像形成装置。
- 10.前記堅固な機械的構造(65)は原点(O)を通過する軸線(67)のま わりに枢動することを特徴とする請求の範囲第9項に記載の2次元の減衰測定か らの3次元画像形成装置。
- 11.前記堅固な機械的構造(75)は不動でありかつ物体(11)がそれに固 定されかつ原点(O)を通過する少なくとも1つの軸線に応じて回転を行なわせ ることができる他の機械的構造(76〜86)を含むことを特徴とする請求の範 囲第9項に記載の2次元の減衰測定からの3次元面像形成方法。
- 12.焦点(S)を含んでいる円錐状放射線源(10)から形成される装置およ びセンサ(19)のネツトワークから形成される2次元検出器(12)を使用す ることにより物体(11)を通る放射線の2次元の減衰測定から前記物体の3次 元光画像形成方法において、前記物体(11)を示す整合の複数の第1点(M) の各々において、放射線の減衰f(M)が前記物体の第2整合の複数の点(U) の放射線減衰の変換ラドンの誘導体を示す量(R′f(U))を計算することに より計算され、作用の変換ラドンが該作用が定義される範囲の少なくとも1点を 通過する各面に関するこの作用のすべての局部的な値として定義され、そして変 換ラドンの誘導体がもしも運動が球面摩擦系によつて定義される通常のベクトル の方向において前記平面に対して垂直に発生するならば前記平面の各々に関する 変化率の合計として定義され、量(R′f(U))は前記円錐状放射線の焦点( S)を通過する平面との前記検出器(12)の交差により得られた少なくとも1 つのラインに沿つてかつ前記第2点(U)に近接して前記放射線減衰の変化の各 2点(U)についての合計を行なうことにより計算され、原点(O)および第2 点(U)を通過する面線は焦点(S)を通過する平面に対してほぼ直交しており 、次いでこれらの合計の直線結合および第1点(M)の各々における放射線減衰 は原点(O)に対する距離に関連するこれらの量の微分により、かつ最後に導出 された量の直線結合(式6)によつて得られ、さらに、間挿が前記第2点(U) から前記第1点(M)へ通過するように行なわれることを特徴とする3次元光画 像形成方法。
- 13.前記量(Rf(U))は各々前記検出器(12)および対応する第2点( U)に近接する間挿点(U11,U12,U21,U22)ならびに円錐状放射 線の焦点(S)の種々の位置を通過する平面の交差によつて得られるラインに沿 う減衰変化の合計により決定され、原点(O)および各間挿点(U11,U12 ,U21,U22)によつて定義される直線は前記間挿点を通過する平面に対し て直交することを特徴とする請求の範囲第12項に記載の3次元光画像形成方法 。
- 14.前記測定は円錐状放射線の焦点(S)が規則的な角度(2π/Nm)によ り離れておりかつ原点(O)を通過する軸線において集束する子午線平面上に発 生するとき行なわれることを特徴とする請求の範囲第12項または第13項に配 置の3次元光画像形成装置。
- 15.前記放射線源(10)および前記検出器(12)は前記物体(11)のま わりの完全な回転にわたつて少なくとも2つの周期を含んでいるほぼ正弦曲線の 形で原点(O)から一定の距離において2つの軌道(Tc,Tc′)を横切るこ とを特徴とする請求の範囲第14項に記載の3次元光面像形成方法。
- 16.前記放射線源(10)および前記検出器(12)は物体(11)のまわり の完全な回転にわたつて少なくとも2つの周期を含んでいるほぼ正弦曲線の形で 原点(O)から一定の距離において2つの軌道(Tc,Tc′)を横切りそして その振幅(Int)は原点(O)と物体(11)のすべての点との間の距離に等 しいかまたはそれより大きく、さらに、前記軌道(Tc)と原点(O)との間の 距離は√3を掛けたこの振幅に等しいかまたはそれより大きいことを特徴とする 請求の範囲第14項に記載の3次元光画像形成方法。
- 17.それに関して量R′f(U)が決定される第2点(U)は等しく分布され た球面座標((n),θ(1),(m))を有しかつ円錐状放射線の焦点(S) が照射のときに発生する子午線面に対して垂直な平面に属することを特徴とする 請求の範囲第14項ないし第16項のいずれか1項に記載の3次元光面像形成方 法。
- 18.前記導出された量は重み付けされ、間挿されかつ各子午線面の内部で該子 午線面と合流される平面の再配置投射の平面にわたつて均一に分布された矩形座 標((r(a),z(b))を有する第3整合に属しかつ物体を示す整合点(M )と子午線面に対して直交する同一平面に属する第3点(B)において量Q(ψ ,B)を得るように加算され、そしてまた、第3整合のいずれかの点(B)に関 して、該点(B)を特徴づける量(Q,(ψ,B))の計算が再配置投影の同一 平面に属する第3整合の他の点の計算と結合されることを特徴とする請求の範囲 第17項に記載の3次元光画像形成装置。
- 19.代表的な整合の点(M)での減衰f(M)は再配置投射による平面上の点 (M)の直交投射の前記第3点(B)に関連づけられる量(Q,ψ,B)の直線 結合によつて得られることを特徴とする請求の範囲第18項に記載の3次元光面 像形成装置。
- 20.そのデカルト座標が一般に均一に分布される物体(11)の第1の3次元 整合の第1点(M)上のパラメータ(f(M))の計算方法において、まず、独 特な点の第2および第3の3次元整を構成する第2点(U)および第3点(B) を定定し前記第2点(U)の球面座標が均一に分布されかつ前記第2点(U)が 軸線において集束する子午線面にとくに属し、その円柱摩擦が均一に分布する前 記第3点(B)が前記第1点(M)を含みかつ軸線に対して直交する子午線面お よび平行面双方に属し;次いで前記第2点(U)に関する情報(R′f(U)) を獲得し;次いでこれからデジタル処理により導出情報を引き出しかつ次いで前 記第3点(B)に関する中間情報(Q(ψ,B))をこれから引き出すように同 一子午線面に属する前記第2点(U)のクループに関する前記導出情報を結し; そして最後に、これから前記第1点(M)に関するパラメータ(f(M))を引 き出すために同一平行平面に属する前記第3点(B)のグループに関する中間情 報(Q(ψ,B))を結合することからなることを特徴とするパラメータ計算方 法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR8707134A FR2615619B1 (fr) | 1987-05-21 | 1987-05-21 | Procede et dispositif d'imagerie tridimentionnelle a partir de mesures bidimensionnelles de l'attenuation d'un rayonnement |
FR87/07134 | 1987-05-21 |
Publications (1)
Publication Number | Publication Date |
---|---|
JPH01503490A true JPH01503490A (ja) | 1989-11-22 |
Family
ID=9351297
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP63504482A Pending JPH01503490A (ja) | 1987-05-21 | 1988-05-20 | 2次元の放射線減衰測定から3次元の光画像形成を得る方法および装置 |
Country Status (8)
Country | Link |
---|---|
US (1) | US5124914A (ja) |
EP (1) | EP0292402B1 (ja) |
JP (1) | JPH01503490A (ja) |
CA (1) | CA1305566C (ja) |
DE (1) | DE3885650T2 (ja) |
ES (1) | ES2047570T3 (ja) |
FR (1) | FR2615619B1 (ja) |
WO (1) | WO1988009148A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011516870A (ja) * | 2008-04-07 | 2011-05-26 | カール ザイス インダストリエル メステクニーク ゲゼルシャフト ミット ベシュレンクテル ハフツング | 機械的ワークを断層撮影法によって測定するための方法 |
Families Citing this family (57)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0454682A (ja) * | 1990-06-22 | 1992-02-21 | Toshiba Corp | 立体画像処理方法及びその装置 |
US5075554A (en) * | 1990-09-27 | 1991-12-24 | Siemens Gammasonics, Inc. | Scintillation camera gantry supporting a plurality of detector heads between two parallel plates |
FR2668829B1 (fr) * | 1990-11-05 | 1993-10-22 | Commissariat Energie Atomique | Dispositif et procede de controle non destructif a acquisition simultanee de donnees radiographiques et de donnees tomographiques. |
FR2670038B1 (fr) * | 1990-11-29 | 1993-12-24 | Commissariat A Energie Atomique | Procede et dispositif de reconstruction d'images tridimentionnelles d'un objet en utilisant deux trajectoires circulaires d'axe commun. |
FR2670039B1 (fr) * | 1990-11-29 | 1993-12-24 | Commissariat A Energie Atomique | Procede et dispositif de reconstruction d'images tridimentionnelles d'un objet en utilisant deux trajectoires circulaires d'acquisition. |
US5404293A (en) * | 1991-06-11 | 1995-04-04 | The University Of Utah | Cone beam reconstruction using helical data collection paths |
US5463721A (en) * | 1991-07-29 | 1995-10-31 | General Electric Company | Method for constructing a three dimensional scanning trajectory capable of acquiring a complete set of radon data for exact image reconstruction of an object irradiated by a cone beam source |
US5309520A (en) * | 1991-12-16 | 1994-05-03 | General Electric Co. | Signal processing system using recursive radon video transform processor |
US5828723A (en) * | 1995-02-15 | 1998-10-27 | Tomografia De Hormigon Armado S.A. | Process for determining the internal three-dimensional structure of a body opaque to visible light by means of radiations from a single source, specially suitable for reinforced concrete parts |
US5647360A (en) * | 1995-06-30 | 1997-07-15 | Siemens Corporate Research, Inc. | Digital subtraction angiography for 3D diagnostic imaging |
US5690106A (en) * | 1995-06-30 | 1997-11-25 | Siemens Corporate Research, Inc. | Flexible image registration for rotational angiography |
US6256369B1 (en) * | 1999-03-31 | 2001-07-03 | Analogic Corporation | Computerized tomography scanner with longitudinal flying focal spot |
JP2000315833A (ja) * | 1999-04-30 | 2000-11-14 | Fuji Photo Film Co Ltd | 単一縦モード固体レーザー |
US6822658B1 (en) * | 1999-10-25 | 2004-11-23 | Intel Corporation | Rendering a silhouette edge |
US6650324B1 (en) * | 1999-10-29 | 2003-11-18 | Intel Corporation | Defining surface normals in a 3D surface mesh |
US6798411B1 (en) * | 1999-10-29 | 2004-09-28 | Intel Corporation | Image processing |
US6762759B1 (en) | 1999-12-06 | 2004-07-13 | Intel Corporation | Rendering a two-dimensional image |
US7098925B1 (en) | 2000-03-10 | 2006-08-29 | Intel Corporation | Shading of images using texture |
US7180523B1 (en) | 2000-03-31 | 2007-02-20 | Intel Corporation | Trimming surfaces |
US7061501B1 (en) | 2000-11-07 | 2006-06-13 | Intel Corporation | Rendering a pencil-sketch image |
DE10063442A1 (de) * | 2000-12-20 | 2002-07-04 | Philips Corp Intellectual Pty | Verfahren und Röntgeneinrichtung zur Ermittlung eines Satzes von Projektionsabbildungen eines Untersuchungsobjektes |
US6603471B2 (en) | 2001-01-31 | 2003-08-05 | Intel Corporation | Reducing detail in animated three-dimensional models |
US6593927B2 (en) | 2001-02-01 | 2003-07-15 | Intel Corporation | Squash and stretch three-dimensional rendering |
US7190374B2 (en) | 2001-02-28 | 2007-03-13 | Intel Corporation | Shading polygons from a three-dimensional model |
US7116330B2 (en) * | 2001-02-28 | 2006-10-03 | Intel Corporation | Approximating motion using a three-dimensional model |
US6924798B2 (en) * | 2001-05-22 | 2005-08-02 | Intel Corporation | Real-time multi-resolution shadows |
US20020175911A1 (en) * | 2001-05-22 | 2002-11-28 | Light John J. | Selecting a target object in three-dimensional space |
US6772175B2 (en) | 2001-05-31 | 2004-08-03 | Intel Corporation | Database that stores data for a three-dimensional mesh |
US6980206B2 (en) * | 2001-06-07 | 2005-12-27 | Intel Corporation | Rendering a three-dimensional model using a dither pattern |
US6798415B2 (en) | 2001-06-21 | 2004-09-28 | Intel Corporation | Rendering collisions of three-dimensional models |
US6924804B2 (en) * | 2001-09-25 | 2005-08-02 | Intel Corporation | Reducing the resolution of bones in a three-dimensional model |
US6865246B2 (en) | 2001-09-26 | 2005-03-08 | Massachusetts Institute Of Technology | True 3D cone-beam imaging method and apparatus |
JP2005504571A (ja) | 2001-09-26 | 2005-02-17 | マサチューセッツ インスティチュート オブ テクノロジー | 多機能コーンビーム結像装置とその方法 |
US6906724B2 (en) * | 2001-10-17 | 2005-06-14 | Lntel Corporation | Generating a shadow for a three-dimensional model |
US20030080937A1 (en) * | 2001-10-30 | 2003-05-01 | Light John J. | Displaying a virtual three-dimensional (3D) scene |
US7548241B2 (en) * | 2002-01-04 | 2009-06-16 | Intel Corporation | Determining a node path through a node graph |
US7301547B2 (en) * | 2002-03-22 | 2007-11-27 | Intel Corporation | Augmented reality system |
US7146297B2 (en) * | 2002-03-27 | 2006-12-05 | Intel Corporation | Detecting collisions of three-dimensional models |
US6975318B2 (en) * | 2002-06-25 | 2005-12-13 | Intel Corporation | Polygon binning process for tile-based rendering |
US6982715B2 (en) | 2002-07-26 | 2006-01-03 | Intel Corporation | Mesh compression process |
US7215810B2 (en) * | 2003-07-23 | 2007-05-08 | Orametrix, Inc. | Method for creating single 3D surface model from a point cloud |
US7170966B2 (en) * | 2003-08-05 | 2007-01-30 | Gioietta Kuo-Petravic | Practical implementation of a CT cone beam algorithm for 3-D image reconstruction as applied to nondestructive inspection of baggage, live laboratory animal and any solid materials |
JP5129480B2 (ja) | 2003-09-25 | 2013-01-30 | パイエオン インコーポレイテッド | 管状臓器の3次元再構成を行うシステム及び血管撮像装置の作動方法 |
EP1706702A2 (en) * | 2003-12-21 | 2006-10-04 | KREMEN, Stanley H. | System and apparatus for recording, transmitting, and projecting digital three-dimensional images |
JP2008541982A (ja) * | 2005-06-07 | 2008-11-27 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | コーンビームct用高速再構成アルゴリズム |
US7953260B2 (en) * | 2006-06-09 | 2011-05-31 | Craniosim Solutions, Inc. | Predicting movement of soft tissue of the face in response to movement of underlying bone |
GB2475722B (en) * | 2009-11-30 | 2011-11-02 | Mirada Medical | Measurement system for medical images |
US8437522B2 (en) | 2011-02-17 | 2013-05-07 | Kabushiki Kaisha Toshiba | Motion index for medical imaging data based upon Grangeat's formula |
US20180217051A1 (en) * | 2015-06-02 | 2018-08-02 | Centre National De La Recherche Scientifique - Cnrs | Acoustic-optical imaging methods and systems |
US10336465B2 (en) * | 2016-01-08 | 2019-07-02 | The Regents Of The University Of Michigan | Ice crystals and volcanic ash detection system |
US10845615B2 (en) | 2016-09-29 | 2020-11-24 | Marel Iceland Ehf | Method of generating a three dimensional surface profile of a food object |
US10621865B2 (en) | 2018-03-29 | 2020-04-14 | The Regents Of The University Of Michigan | Road condition monitoring system |
US10508952B1 (en) | 2018-10-31 | 2019-12-17 | The Regents Of The University Of Michigan | Optimum spectral bands for active vision systems |
US11127106B2 (en) | 2019-06-28 | 2021-09-21 | Intel Corporation | Runtime flip stability characterization |
US11409341B2 (en) | 2019-10-01 | 2022-08-09 | Intel Corporation | Repeating graphics render pattern detection |
US11809161B2 (en) * | 2020-07-13 | 2023-11-07 | Lawrence Livermore National Security, Llc | Computed axial lithography optimization system |
CN114354659B (zh) * | 2021-12-21 | 2024-01-09 | 昆山善思光电科技有限公司 | 一种x光射线光源在x光机检测仪中的衰减检测方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE1033844B (de) * | 1951-09-04 | 1958-07-10 | Ialicenciaia Talalmanyokat Ert | Roentgenapparat fuer Bewegungsbestrahlung |
FR1195082A (fr) * | 1958-05-09 | 1959-11-13 | Installation pour examen radiologique | |
DE2726635C3 (de) * | 1977-06-14 | 1981-05-27 | Philips Patentverwaltung Gmbh, 2000 Hamburg | Tomographische Vorrichtung zur Ermittlung der Absorption einer Strahlung in einem dreidimensionalen Untersuchungsbereich |
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 |
FR2519445B1 (fr) * | 1981-12-31 | 1985-12-27 | Thomson Csf | Procede d'imagerie tomodensitometrique en mode dynamique, et tomodensitometre mettant en oeuvre un tel procede |
US4578753A (en) * | 1983-08-29 | 1986-03-25 | Elscint Ltd. | Systems and methods for minimizing noncoplanarity artifacts |
-
1987
- 1987-05-21 FR FR8707134A patent/FR2615619B1/fr not_active Expired - Lifetime
-
1988
- 1988-05-20 JP JP63504482A patent/JPH01503490A/ja active Pending
- 1988-05-20 ES ES88401234T patent/ES2047570T3/es not_active Expired - Lifetime
- 1988-05-20 US US07/310,726 patent/US5124914A/en not_active Expired - Lifetime
- 1988-05-20 EP EP88401234A patent/EP0292402B1/fr not_active Expired - Lifetime
- 1988-05-20 WO PCT/FR1988/000254 patent/WO1988009148A1/fr unknown
- 1988-05-20 CA CA000567376A patent/CA1305566C/fr not_active Expired - Lifetime
- 1988-05-20 DE DE3885650T patent/DE3885650T2/de not_active Expired - Lifetime
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011516870A (ja) * | 2008-04-07 | 2011-05-26 | カール ザイス インダストリエル メステクニーク ゲゼルシャフト ミット ベシュレンクテル ハフツング | 機械的ワークを断層撮影法によって測定するための方法 |
Also Published As
Publication number | Publication date |
---|---|
DE3885650T2 (de) | 1994-05-19 |
DE3885650D1 (de) | 1993-12-23 |
US5124914A (en) | 1992-06-23 |
EP0292402A1 (fr) | 1988-11-23 |
FR2615619B1 (fr) | 1991-07-19 |
CA1305566C (fr) | 1992-07-21 |
ES2047570T3 (es) | 1994-03-01 |
WO1988009148A1 (fr) | 1988-12-01 |
EP0292402B1 (fr) | 1993-11-18 |
FR2615619A1 (fr) | 1988-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JPH01503490A (ja) | 2次元の放射線減衰測定から3次元の光画像形成を得る方法および装置 | |
US7186023B2 (en) | Slice image and/or dimensional image creating method | |
KR101252010B1 (ko) | 라돈데이터로부터 (n+1)차원 영상 함수를 재구성하는방법과 장치 | |
CN105806858B (zh) | Ct检测方法和ct设备 | |
US5598453A (en) | Method for X-ray fluoroscopy or radiography, and X-ray apparatus | |
JP4415762B2 (ja) | 断層撮影装置 | |
US20070116177A1 (en) | Imaging system | |
CN100457040C (zh) | 同步辐射x射线相位衬度ct成像装置及实验方法 | |
US6038282A (en) | X-ray imaging system | |
JP7198457B2 (ja) | インテリアct位相イメージングx線顕微鏡装置 | |
JPH07120411B2 (ja) | 円錐ビーム投射データから物体の三次元計算機断層撮影映像を生成する方法および装置 | |
JPH0910202A (ja) | 診断撮像方法と装置およびヘリカル部分コーンビームデータからの画像再構成 | |
US6944263B2 (en) | Apparatus and methods for multiple view angle stereoscopic radiography | |
JP2007216052A (ja) | 胸部撮像のための断層合成装置 | |
CN101542240A (zh) | C型臂计算机断层成像*** | |
JP2007512034A (ja) | 発散ビームスキャナのための画像再構成方法 | |
WO2003027954B1 (en) | Versatile cone-beam imaging apparatus and method | |
JPS61181443A (ja) | 回転式扇形ビ−ムct装置の回転中心のずれを補正する方法と装置 | |
JPH05508093A (ja) | 特にトモグラフィーのための、体の検査を扶助する方法および装置 | |
JPH0943354A (ja) | 物体の三次元画像の再構成方法 | |
US5705819A (en) | Emission CT apparatus | |
CN103505237A (zh) | 经由计算机断层摄影的定量二维荧光透视 | |
IL105732A (en) | Process for reconstructing three-dimensional images of bone by measurements with the help of conical radiation and a two-dimensional detector system | |
McMahon et al. | X-ray tomographic imaging of the complex refractive index | |
JPH04307036A (ja) | 3次元像再構成方法およびその装置 |