JPH04178591A - Ｒｉ分布イメージの位置合わせ方法 - Google Patents
Ｒｉ分布イメージの位置合わせ方法Info
- Publication number
- JPH04178591A JPH04178591A JP2306335A JP30633590A JPH04178591A JP H04178591 A JPH04178591 A JP H04178591A JP 2306335 A JP2306335 A JP 2306335A JP 30633590 A JP30633590 A JP 30633590A JP H04178591 A JPH04178591 A JP H04178591A
- Authority
- JP
- Japan
- Prior art keywords
- images
- addition
- image
- profile
- profiles
- 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
Landscapes
- Nuclear Medicine (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔産業上の利用分野〕
本発明は、ディジタル式のシンチレーションカメラ（デ
ィジタルカメラ）などから取り込んだ、同一被検体の同
一部位より時間を異にして得た任意の２つのディジタル
ＲＩ（ラジオアイソトープ）分布イメージの位置合わせ
方法に関するものである。
ィジタルカメラ）などから取り込んだ、同一被検体の同
一部位より時間を異にして得た任意の２つのディジタル
ＲＩ（ラジオアイソトープ）分布イメージの位置合わせ
方法に関するものである。
被検体にＲＩを注入し、特定の臓器への集積をシンチレ
ーションカメラによりイメージとして検出する、特に心
筋イメージングや骨シンチグラムにおいては、ＲＩの注
入直後、２〜３時間及び２４時間後といように、時間が
経過したＲＩ分布イメージ（以下、イメージと略称する
。）を収集する必要がある。
ーションカメラによりイメージとして検出する、特に心
筋イメージングや骨シンチグラムにおいては、ＲＩの注
入直後、２〜３時間及び２４時間後といように、時間が
経過したＲＩ分布イメージ（以下、イメージと略称する
。）を収集する必要がある。
また核医学検査では、手術前２手術後の経過、発病後の
経過などを診るため、時間が経過したイメージを収集す
る必要がある。
経過などを診るため、時間が経過したイメージを収集す
る必要がある。
このように、同一被検体の同一部位より時間を異にして
得た複数のイメージを比較観察する場合、従来は、それ
らのイメージ（写真又はフィルム）を単に並べて診断し
ていた。
得た複数のイメージを比較観察する場合、従来は、それ
らのイメージ（写真又はフィルム）を単に並べて診断し
ていた。
しかし近年は、上記イメージがディジタルデータで得ら
れる、ディジタル式のシンチレーションカメラ（ディジ
タルカメラ）が出現し、これによる定量的な診断方法が
増えるに従ってイメージ（ディジタルイメージをいう。
れる、ディジタル式のシンチレーションカメラ（ディジ
タルカメラ）が出現し、これによる定量的な診断方法が
増えるに従ってイメージ（ディジタルイメージをいう。
以下間し。）の位置合わせが重要になってきた。
従来、この種のイメージの位置合わせ方法として、１つ
はイメージ収集時にカメラと被検体の位置を前回収集時
と同一位置に合わせて収集する方法、もう１つは収集後
にイメージの位置合わせをする方法があった。
はイメージ収集時にカメラと被検体の位置を前回収集時
と同一位置に合わせて収集する方法、もう１つは収集後
にイメージの位置合わせをする方法があった。
しかし前者は、手術前１手術後の検査など、イメージ収
集間隔が比較的長期に亘る場合には実用的でなく、用い
られていない。
集間隔が比較的長期に亘る場合には実用的でなく、用い
られていない。
後者には、現在用いられている方法として、（１）２枚
のイメージを交互にモニタに表示し、画面を見ながら、
一方のイメージを上、下、左。
のイメージを交互にモニタに表示し、画面を見ながら、
一方のイメージを上、下、左。
右に移動し、位置合わせする方法。
（２）両者のイメージ上に、２点を設定し、その２点を
合わせるようにイメージを移動する方法、があるが、共
に手動での操作が必要になる。
合わせるようにイメージを移動する方法、があるが、共
に手動での操作が必要になる。
上記従来技術は、手動による操作が必要で、操作が面倒
である上に、定量化を行う上での再現性にも問題があっ
た。
である上に、定量化を行う上での再現性にも問題があっ
た。
本発明の目的は、同一被検体の同一部位より時間を異に
して得た任意の２つのイメージを簡単かつ再現性よく位
置合わせすることのできるＲＩ分布イメージの位置合わ
せ方法を提供することにある。
して得た任意の２つのイメージを簡単かつ再現性よく位
置合わせすることのできるＲＩ分布イメージの位置合わ
せ方法を提供することにある。
上記目的は、位置合わせする２つのイメージ（Ａ、Ｂ）
について、各々Ｘ方向及びＸ方向にデータを加算して両
イメージのＸ、Ｘ方向の加算プロフィール（ＡＸ、ＡＹ
、ＢＸ、ＢＹ）を得、その両イメージの加算プロフィー
ルをＸ、Ｘ方向毎に重ね合わせて各々Ｘ、Ｙ方向に相対
移動させつつ両プロフィールの差分（ＡＸ−ＢＸ、ＡＹ
−ＢＹ）の合計を求め、各々その差分の合計値が最小と
なる位置をもって位置合わせ完了とすることにより達成
される。
について、各々Ｘ方向及びＸ方向にデータを加算して両
イメージのＸ、Ｘ方向の加算プロフィール（ＡＸ、ＡＹ
、ＢＸ、ＢＹ）を得、その両イメージの加算プロフィー
ルをＸ、Ｘ方向毎に重ね合わせて各々Ｘ、Ｙ方向に相対
移動させつつ両プロフィールの差分（ＡＸ−ＢＸ、ＡＹ
−ＢＹ）の合計を求め、各々その差分の合計値が最小と
なる位置をもって位置合わせ完了とすることにより達成
される。
２つのイメージ（Ａ、Ｂ）について、各々Ｘ方向及びＸ
方向にデータを加算することにより、両イメージのＸ、
Ｘ方向の輪郭の投影である加算プロフィール（ＡＸ、Ａ
Ｙ、ＢＸ、ＢＹ）が得られる。これを、Ｘ、Ｘ方向毎に
重ね合わせて各々Ｘ。
方向にデータを加算することにより、両イメージのＸ、
Ｘ方向の輪郭の投影である加算プロフィール（ＡＸ、Ａ
Ｙ、ＢＸ、ＢＹ）が得られる。これを、Ｘ、Ｘ方向毎に
重ね合わせて各々Ｘ。
Ｘ方向に相対移動させつつ両プロフィールの差分（ＡＸ
−ＢＸ、ＡＹ−ＢＹ）の合計を求める。
−ＢＸ、ＡＹ−ＢＹ）の合計を求める。
ここで、この差分の合計値が大きいほど両イメージの重
なり程度が小さく、逆に差分の合計値が小さいほど両イ
メージの重なり程度が大きいことになる。したがって、
両プロフィールの差分の合計値が最小となる位置が求め
られたら、その位置をもって位置合わせ完了とすればよ
い。具体的には、その位置にいずれか一方のイメージを
移動させればよい。
なり程度が小さく、逆に差分の合計値が小さいほど両イ
メージの重なり程度が大きいことになる。したがって、
両プロフィールの差分の合計値が最小となる位置が求め
られたら、その位置をもって位置合わせ完了とすればよ
い。具体的には、その位置にいずれか一方のイメージを
移動させればよい。
以上の操作は、通常の画像処理、演算手段にて容易に可
能であり、したがって、本発明方法によれば、任意の２
つ（複数）のイメージを自動的に。
能であり、したがって、本発明方法によれば、任意の２
つ（複数）のイメージを自動的に。
すなわち簡単かつ再現性よく、位置合わせできることに
なる。
なる。
以下、図面を参照して本発明の詳細な説明する。第１図
は本発明によるＲＩ分布イメージの位置合わせ方法の一
例を説明するためのフローチャート、第２図はイメージ
ＡとそのＸ、Ｘ方向の加算プロフィールＡＸ、ＡＹを示
す図、第３図はイメージＢとそのＸ、Ｘ方向の加算プロ
フィールＢＸ、ＢＹを示す図、第４図はイメージＡ、Ｈ
の加算プロフィールＡＸ、ＢＸを重ね合わせてＸ。
は本発明によるＲＩ分布イメージの位置合わせ方法の一
例を説明するためのフローチャート、第２図はイメージ
ＡとそのＸ、Ｘ方向の加算プロフィールＡＸ、ＡＹを示
す図、第３図はイメージＢとそのＸ、Ｘ方向の加算プロ
フィールＢＸ、ＢＹを示す図、第４図はイメージＡ、Ｈ
の加算プロフィールＡＸ、ＢＸを重ね合わせてＸ。
Ｘ方向に相対移動させている状態を示す図、第５図は上
記両プロフィールの差分（ＡＸ−ＢＸ）の合計値が最小
Ｘとなった状態を示す図、第６図はイメージＡ、Ｂのプ
ロフィールの差分（ＡＸ−ＢＸ、ＡＹ−ＢＹ）の合計値
が最小ＸＴ３’となる位置にイメージＢを移動させる様
子を示す図である。
記両プロフィールの差分（ＡＸ−ＢＸ）の合計値が最小
Ｘとなった状態を示す図、第６図はイメージＡ、Ｂのプ
ロフィールの差分（ＡＸ−ＢＸ、ＡＹ−ＢＹ）の合計値
が最小ＸＴ３’となる位置にイメージＢを移動させる様
子を示す図である。
第７図は本発明方法が適用されたディジタルカメラ（デ
ィジタル式のシンチレーションカメラ）の−例を示すブ
ロック図である。このディジタルカメラは、被検体の体
内に投与したＲＩから放出されるγ線を検出して電気信
号に変換することにより上記被検体のイメージ信号を収
集するシンチレーションカメラ本体１と、このカメラ本
体１からのイメージ信号をディジタル化し、ディジタル
データとして出力するＡ／Ｄ変換器２とを備える。
ィジタル式のシンチレーションカメラ）の−例を示すブ
ロック図である。このディジタルカメラは、被検体の体
内に投与したＲＩから放出されるγ線を検出して電気信
号に変換することにより上記被検体のイメージ信号を収
集するシンチレーションカメラ本体１と、このカメラ本
体１からのイメージ信号をディジタル化し、ディジタル
データとして出力するＡ／Ｄ変換器２とを備える。
また、このＡ／Ｄ変換器２からのイメージデータを格納
しておく記憶装置３と、上記イメージデータを記憶する
イメージメモリ４と、本発明方法による加算動作を行う
加算回路５と、イメージメモリ４に記憶されたイメージ
データを表示するＣＲＴなどのデイスプレィ６と、これ
ら各部を制御するＣＰｔＪ７とを備える。
しておく記憶装置３と、上記イメージデータを記憶する
イメージメモリ４と、本発明方法による加算動作を行う
加算回路５と、イメージメモリ４に記憶されたイメージ
データを表示するＣＲＴなどのデイスプレィ６と、これ
ら各部を制御するＣＰｔＪ７とを備える。
いま、同一被検体の同一部位を時間を異にして撮影して
得た任意の２つのイメージ（データ）ＡとＢが、記憶装
置３に格納されているものとする。
得た任意の２つのイメージ（データ）ＡとＢが、記憶装
置３に格納されているものとする。
本実施例では、まずイメージＡを記憶装置３から読み出
しくステップ１０１参照）、加算回路５でそのＸ方向の
画素につきデータを加算して、Ｘ方向の加算プロフィー
ルＡＸを作る（ステップ１０２参照）。同様にして、Ｙ
方向の加算プロフィールＡＹを作る（ステップ１０３参
照）（以上。
しくステップ１０１参照）、加算回路５でそのＸ方向の
画素につきデータを加算して、Ｘ方向の加算プロフィー
ルＡＸを作る（ステップ１０２参照）。同様にして、Ｙ
方向の加算プロフィールＡＹを作る（ステップ１０３参
照）（以上。
第２図参照）。
次に、イメージＢを記憶装置３から読み出しくステップ
１０４参照）、上記と同様にして加算プロフィールＢＸ
、ＢＹを作る（ステップ１０５゜１０６参照）（以上、
第３図参照）。
１０４参照）、上記と同様にして加算プロフィールＢＸ
、ＢＹを作る（ステップ１０５゜１０６参照）（以上、
第３図参照）。
次に、加算プロフィールＡＸを基準とし、加算プロフィ
ールＢＸを１画素ずつ正方向及び負方向（図中、左右方
向）へ移動させ、それらの差分の合計値が最小となると
きの移動量Ｘを求める（ステップ１０７．第４図及び第
５図参照）。同様に。
ールＢＸを１画素ずつ正方向及び負方向（図中、左右方
向）へ移動させ、それらの差分の合計値が最小となると
きの移動量Ｘを求める（ステップ１０７．第４図及び第
５図参照）。同様に。
加算プロフィールＡＹを基準とし、加算プロフィールＢ
Ｙを移動させ、それらの差分の合計値が最小となるとき
の移動量ｙを求める（ステップ１０８参照）。
Ｙを移動させ、それらの差分の合計値が最小となるとき
の移動量ｙを求める（ステップ１０８参照）。
最後に、ステップ１０７，１０８で算出した移動量Ｘ、
ｙ画素分、イメージＢを移動することにより、両イメー
ジＡとＢが位置合わせされる（ステップ１０９．第６図
参照）。
ｙ画素分、イメージＢを移動することにより、両イメー
ジＡとＢが位置合わせされる（ステップ１０９．第６図
参照）。
ここで、ステップ１０７〜１０９は、全てＣＰＵ７によ
って処理される。
って処理される。
上述実施例では、イメージの加算をステップ１０２．１
０３，１０５及び１０６において加算回路５で行うこと
としたが、加算回路５を用いずＣＰＵ７に対するプログ
ラムにより行うようにしてもよい。
０３，１０５及び１０６において加算回路５で行うこと
としたが、加算回路５を用いずＣＰＵ７に対するプログ
ラムにより行うようにしてもよい。
なお、２つのイメージＡ、Ｂが完全に重ね合わされると
零になるのは、両イメージが全く同じ加算プＯ’フィー
ＪＬ／ （ＡＸ＝ＢＸ、ＡＹ＝ＢＹ）をもっている場合
である。しかし、時間を異にして得た２つのイメージは
、通常、全く同じプロフィールをもっていることはなく
、したがってここでは、両加算プロフィールの差分の合
計値が最小となった位置を２つのイメージが重ね合わさ
れた位置としている。時間を異にして得た２つのイメー
ジについても、両加算プロフィールの差分が零になった
ときに２つのイメージが重ね合わされたとすると、いず
れか一方のイメージのプロフィールを補正しなければな
らいが１本発明方法によれば、このような補正を必要と
しない。
零になるのは、両イメージが全く同じ加算プＯ’フィー
ＪＬ／ （ＡＸ＝ＢＸ、ＡＹ＝ＢＹ）をもっている場合
である。しかし、時間を異にして得た２つのイメージは
、通常、全く同じプロフィールをもっていることはなく
、したがってここでは、両加算プロフィールの差分の合
計値が最小となった位置を２つのイメージが重ね合わさ
れた位置としている。時間を異にして得た２つのイメー
ジについても、両加算プロフィールの差分が零になった
ときに２つのイメージが重ね合わされたとすると、いず
れか一方のイメージのプロフィールを補正しなければな
らいが１本発明方法によれば、このような補正を必要と
しない。
以上説明したように上述実施例によれば、従来、オペレ
ータが手動で行っていたイメージの位置合わせが自動で
行え、操作が容易になるばかりか、オペレヘータによる
ばらつきがなくなり、定量診断の上でも有効である。
ータが手動で行っていたイメージの位置合わせが自動で
行え、操作が容易になるばかりか、オペレヘータによる
ばらつきがなくなり、定量診断の上でも有効である。
本発明によれば、同一被検体の同一部位より時間を異に
して得た任意の２つのイメージの位置合わせが自動化で
きて手動による操作が省略できる。
して得た任意の２つのイメージの位置合わせが自動化で
きて手動による操作が省略できる。
したがって、イメージの位置合わせ操作が極めて簡単に
なり、また、オペレータによる解析結果のばらつきがな
くなって再現性が向上するという効果がある。
なり、また、オペレータによる解析結果のばらつきがな
くなって再現性が向上するという効果がある。
第１図は本発明方法の一例を説明するためのフローチャ
ート、第２図〜第６図は各々第１図に示す本発明方法の
説明図、第７図は本発明方法が適用されたディジタルカ
メラの一例を示すブロック図である。 Ａ、Ｂ・・・位置合わせされるイメージ（ディジタルの
ＲＩ分布イメージ）、ＡＸ、ＡＹ・・・イメージＡの加
算プロフィール、ＢＸ、ＢＹ・・・イメージＢの加算プ
ロフィール、Ｘ？ ｙ・・・位置合わせのためのイメー
ジ移動量、１・・・シンチレーションカメラ本体、２・
・・Ａ／Ｄ変換器、３・・・記憶装置、４・・・イメー
ジメモリ、５・・・加算回路、６・・・デイスプレィ、
第１区
ート、第２図〜第６図は各々第１図に示す本発明方法の
説明図、第７図は本発明方法が適用されたディジタルカ
メラの一例を示すブロック図である。 Ａ、Ｂ・・・位置合わせされるイメージ（ディジタルの
ＲＩ分布イメージ）、ＡＸ、ＡＹ・・・イメージＡの加
算プロフィール、ＢＸ、ＢＹ・・・イメージＢの加算プ
ロフィール、Ｘ？ ｙ・・・位置合わせのためのイメー
ジ移動量、１・・・シンチレーションカメラ本体、２・
・・Ａ／Ｄ変換器、３・・・記憶装置、４・・・イメー
ジメモリ、５・・・加算回路、６・・・デイスプレィ、
第１区
Claims (1)
- １、同一被検体の同一部位より時間を異にして得た任意
の２つのディジタルＲＩ分布イメージ（Ａ、Ｂ）の位置
合わせ方法において、位置合わせする２つの前記イメー
ジ（Ａ、Ｂ）について、各々Ｘ方向及びＹ方向にデータ
を加算して両イメージのＸ、Ｙ方向の加算プロフィール
（ＡＸ、ＡＹ、ＢＸ、ＢＹ）を得、その両イメージの加
算プロフィールをＸ、Ｙ方向毎に重ね合わせて各々Ｘ、
Ｙ方向に相対移動させつつ両プロフィールの差分（ＡＸ
−ＢＸ、ＡＹ−ＢＹ）の合計を求め、各々その差分の合
計値が最小となる位置をもつて位置合わせ完了とするＲ
Ｉ分布イメージの位置合わせ方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2306335A JPH04178591A (ja) | 1990-11-14 | 1990-11-14 | Ｒｉ分布イメージの位置合わせ方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2306335A JPH04178591A (ja) | 1990-11-14 | 1990-11-14 | Ｒｉ分布イメージの位置合わせ方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH04178591A true JPH04178591A (ja) | 1992-06-25 |
Family
ID=17955870
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2306335A Pending JPH04178591A (ja) | 1990-11-14 | 1990-11-14 | Ｒｉ分布イメージの位置合わせ方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH04178591A (ja) |
-
1990
- 1990-11-14 JP JP2306335A patent/JPH04178591A/ja active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109472829B (zh) | 一种物体定位方法、装置、设备和存储介质 | |
| JPH08131429A (ja) | 管状体像再生方法およびその装置 | |
| US6868172B2 (en) | Method for registering images in a radiography application | |
| CA2106146A1 (en) | Video-cad comparator system | |
| JPS5917332A (ja) | 医用画像重ね合わせ方式 | |
| JPH04178591A (ja) | Ｒｉ分布イメージの位置合わせ方法 | |
| JPS59183461A (ja) | 複合画像設定方法 | |
| CN105719235A (zh) | 基于圆周扫描视频图像拼接与分屏显示方法 | |
| US11948290B2 (en) | Individual channel characterization of collimator | |
| JPH07325908A (ja) | 投影画像処理方法及び投影画像処理装置 | |
| JPH0436684A (ja) | Ｓｐｅｃｔ装置 | |
| JPH0777749A (ja) | 放射線画像情報読取装置 | |
| JP2970674B2 (ja) | 断層像撮影装置 | |
| JP2838522B2 (ja) | ディジタル・サブトラクション・アンギオグラフィ装置 | |
| JPH04261279A (ja) | ディジタル・サブトラクション・アンギオグラフィ装置 | |
| JPS6230983A (ja) | シンチレ−シヨンカメラのイメ−ジ表示方法 | |
| US5391880A (en) | Gamma camera system | |
| JPH1027253A (ja) | 画像処理方法 | |
| JPH10201745A (ja) | Ｘ線撮影装置 | |
| JPS61146078A (ja) | 赤外線映像装置 | |
| Chandler | Height estimation using video security imagery | |
| WO1986007168A1 (fr) | Methode et appareil de photographie medicale | |
| JPH06323821A (ja) | 骨格形状データ作成装置 | |
| JPH1023334A (ja) | Ｘ線映像装置 | |
| CN116051606A (zh) | 影像跟踪采集方法、成像系统及计算机设备 |