JP2791231B2 - Ｘ線画像診断装置の糸巻歪補正方法 - Google Patents
Ｘ線画像診断装置の糸巻歪補正方法Info
- Publication number
- JP2791231B2 JP2791231B2 JP3084956A JP8495691A JP2791231B2 JP 2791231 B2 JP2791231 B2 JP 2791231B2 JP 3084956 A JP3084956 A JP 3084956A JP 8495691 A JP8495691 A JP 8495691A JP 2791231 B2 JP2791231 B2 JP 2791231B2
- Authority
- JP
- Japan
- Prior art keywords
- center
- ray image
- marker
- distortion
- ray
- 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 - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 39
- 238000002059 diagnostic imaging Methods 0.000 title claims description 8
- 240000000136 Scabiosa atropurpurea Species 0.000 title claims 3
- 241000226585 Antennaria plantaginifolia Species 0.000 claims description 74
- 239000003550 marker Substances 0.000 claims description 35
- 239000000463 material Substances 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 2
- 238000003745 diagnosis Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 1
- 235000013372 meat Nutrition 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
- H05G—X-RAY TECHNIQUE
- H05G1/00—X-ray apparatus involving X-ray tubes; Circuits therefor
- H05G1/08—Electrical details
- H05G1/60—Circuit arrangements for obtaining a series of X-ray photographs or for X-ray cinematography
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Nuclear Medicine (AREA)
- Image-Pickup Tubes, Image-Amplification Tubes, And Storage Tubes (AREA)
Description
【０００１】
【産業上の利用分野】この発明は、Ｘ線画像を用いて被
写体の診断を行う装置の糸巻歪補正方法に関し、特に高
精度の画像診断を実現するためのＸ線画像診断装置の糸
巻歪補正方法に関するものである。
写体の診断を行う装置の糸巻歪補正方法に関し、特に高
精度の画像診断を実現するためのＸ線画像診断装置の糸
巻歪補正方法に関するものである。
【０００２】
【従来の技術】一般に、Ｘ線医療診断においては、Ｘ線
管とＸ線像増倍管との間に被写体を配置し、被写体を透
過したＸ線をディジタル化した画像が診断に用いられて
いる。しかし、Ｘ線像増倍管から得られる画像は、電子
レンズ系に特有の幾何学的非一様性（正方形が糸巻形に
歪曲するため糸巻歪と呼ばれる）を含んでいる。これは
Ｘ線画像の高精度化に対する大きな障害となるため、従
来より、糸巻歪を補正するための対策が種々提案されて
いる。
管とＸ線像増倍管との間に被写体を配置し、被写体を透
過したＸ線をディジタル化した画像が診断に用いられて
いる。しかし、Ｘ線像増倍管から得られる画像は、電子
レンズ系に特有の幾何学的非一様性（正方形が糸巻形に
歪曲するため糸巻歪と呼ばれる）を含んでいる。これは
Ｘ線画像の高精度化に対する大きな障害となるため、従
来より、糸巻歪を補正するための対策が種々提案されて
いる。
【０００３】図５は例えば特開平2−10636号公報に記載
された一般的なＸ線画像診断装置を示す構成図である。
図において、１は例えば実線及び破線位置のように矢印
方向に移動可能に配置されてＸ線Ｒを放射するＸ線管、
２はＸ線管１と同軸上に対向配置されてＸ線Ｒを受ける
Ｘ線像増倍管、３はＸ線管１とＸ線像増倍管２との間に
介在される被写体、４はＸ線像増倍管２に接続されてＸ
線Ｒの画像をディジタル化するイメージプロセッサ、５
はイメージプロセッサ４でディジタル化された画像を表
示する表示装置である。
された一般的なＸ線画像診断装置を示す構成図である。
図において、１は例えば実線及び破線位置のように矢印
方向に移動可能に配置されてＸ線Ｒを放射するＸ線管、
２はＸ線管１と同軸上に対向配置されてＸ線Ｒを受ける
Ｘ線像増倍管、３はＸ線管１とＸ線像増倍管２との間に
介在される被写体、４はＸ線像増倍管２に接続されてＸ
線Ｒの画像をディジタル化するイメージプロセッサ、５
はイメージプロセッサ４でディジタル化された画像を表
示する表示装置である。
【０００４】図６は図５内のＸ線像増倍管２を詳細に示
す断面図であり、10はＸ線管１の線源に相当するＸ線焦
点、20はＸ線像増倍管２を構成する以下の各要素21〜26
を封入する真空容器である。21は例えば板状の被写体３
を透過したＸ線Ｒを受けて光電子Ｅを発生する光電陰
極、22は光電子Ｅをクロスオーバ点Ｐに収束するための
複数のグリッド電極、23及び24はクロスオーバ点Ｐを通
過した光電子Ｅに対する電子レンズ系となる前段陽極及
び後段陽極、25は前段陽極23及び後段陽極24の間に介在
された補正用の中間部電極、26は光電子Ｅを受けて発光
する蛍光膜、26ａは蛍光膜26の出力面、ＦＩＤはＸ線焦
点10と光電陰極21との間のＸ線軌道距離である。
す断面図であり、10はＸ線管１の線源に相当するＸ線焦
点、20はＸ線像増倍管２を構成する以下の各要素21〜26
を封入する真空容器である。21は例えば板状の被写体３
を透過したＸ線Ｒを受けて光電子Ｅを発生する光電陰
極、22は光電子Ｅをクロスオーバ点Ｐに収束するための
複数のグリッド電極、23及び24はクロスオーバ点Ｐを通
過した光電子Ｅに対する電子レンズ系となる前段陽極及
び後段陽極、25は前段陽極23及び後段陽極24の間に介在
された補正用の中間部電極、26は光電子Ｅを受けて発光
する蛍光膜、26ａは蛍光膜26の出力面、ＦＩＤはＸ線焦
点10と光電陰極21との間のＸ線軌道距離である。
【０００５】図７はＸ線焦点10と光電陰極21との距離Ｆ
ＩＤを変化させたときの糸巻歪の大きさを示す特性図で
あり、横軸は半径を100%としたときの糸巻歪の歪中心か
ら画像点までの距離、縦軸は糸巻歪の大きさに相当する
幾何学歪の積分値である。ここでは、Ｘ線焦点10と光電
陰極21との距離ＦＩＤを1000mmから700mmまで100mmずつ
近づけたときの特性曲線を示している。図７より、糸巻
歪は、歪中心(電子レンズ軸と光電陰極21との交点)から
の距離が大きくなるにつれて増大し、Ｘ線Ｒの軌道距離
ＦＩＤが小さくなるにつれて増大することが分かる。
ＩＤを変化させたときの糸巻歪の大きさを示す特性図で
あり、横軸は半径を100%としたときの糸巻歪の歪中心か
ら画像点までの距離、縦軸は糸巻歪の大きさに相当する
幾何学歪の積分値である。ここでは、Ｘ線焦点10と光電
陰極21との距離ＦＩＤを1000mmから700mmまで100mmずつ
近づけたときの特性曲線を示している。図７より、糸巻
歪は、歪中心(電子レンズ軸と光電陰極21との交点)から
の距離が大きくなるにつれて増大し、Ｘ線Ｒの軌道距離
ＦＩＤが小さくなるにつれて増大することが分かる。
【０００６】次に、図５〜図７を参照しながら、従来の
Ｘ線画像診断装置の糸巻歪補正方法について説明する。
Ｘ線管１のＸ線焦点10から放射されたＸ線Ｒは、被写体
３を透過してＸ線像増倍管２の光電陰極21に入射する。
これにより、光電陰極21から発生した光電子Ｅは、電子
レンズ系により収束され、クロスオーバ点Ｐを越えた
後、蛍光膜26に照射されて結像される。蛍光膜26の出力
面26ａから発生したＸ線画像は、イメージプロセッサ４
によりディジタル処理されて表示装置５に表示される。
Ｘ線画像診断装置の糸巻歪補正方法について説明する。
Ｘ線管１のＸ線焦点10から放射されたＸ線Ｒは、被写体
３を透過してＸ線像増倍管２の光電陰極21に入射する。
これにより、光電陰極21から発生した光電子Ｅは、電子
レンズ系により収束され、クロスオーバ点Ｐを越えた
後、蛍光膜26に照射されて結像される。蛍光膜26の出力
面26ａから発生したＸ線画像は、イメージプロセッサ４
によりディジタル処理されて表示装置５に表示される。
【０００７】このとき、Ｘ線軌道距離ＦＩＤの変化によ
り、図７のように糸巻歪の積分値が変化するためＸ線画
像(特に外側の画像)の画質が変化し、診断が困難にな
る。そこで、Ｘ線軌道距離ＦＩＤに応じて、中間部電極
25に印加する電位を制御することにより、糸巻歪の積分
値が一定となるようする。しかし、糸巻歪そのものが無
くなるのではなく、又、歪中心がＸ線軸の中心と同軸上
にない場合は、糸巻歪が非対象に残ることになる。
り、図７のように糸巻歪の積分値が変化するためＸ線画
像(特に外側の画像)の画質が変化し、診断が困難にな
る。そこで、Ｘ線軌道距離ＦＩＤに応じて、中間部電極
25に印加する電位を制御することにより、糸巻歪の積分
値が一定となるようする。しかし、糸巻歪そのものが無
くなるのではなく、又、歪中心がＸ線軸の中心と同軸上
にない場合は、糸巻歪が非対象に残ることになる。
【０００８】
【発明が解決しようとする課題】従来のＸ線画像診断装
置の糸巻歪補正方法は以上のように、単に中間部電極25
の印加電位を制御するのみなので、糸巻歪が残存したま
まのＸ線画像を用いて診断しなければならない。又、歪
中心がＸ線中心と同一軸上にない場合は、糸巻歪が非対
象に残存するため、糸巻歪の積分値を一定にすることが
できないという問題点があった。
置の糸巻歪補正方法は以上のように、単に中間部電極25
の印加電位を制御するのみなので、糸巻歪が残存したま
まのＸ線画像を用いて診断しなければならない。又、歪
中心がＸ線中心と同一軸上にない場合は、糸巻歪が非対
象に残存するため、糸巻歪の積分値を一定にすることが
できないという問題点があった。
【０００９】この発明は上記のような問題点を解決する
ためになされたもので、糸巻歪の大きさが歪中心に対し
て対称となる性質に注目し、歪中心を求めて糸巻歪の補
正を容易にしたＸ線画像診断装置の糸巻歪補正方法を得
ることを目的とする。
ためになされたもので、糸巻歪の大きさが歪中心に対し
て対称となる性質に注目し、歪中心を求めて糸巻歪の補
正を容易にしたＸ線画像診断装置の糸巻歪補正方法を得
ることを目的とする。
【００１０】又、この発明は、歪中心に対して対称に補
正係数を算出し、糸巻歪を補正して高精度のＸ線画像を
実現することができるＸ線画像診断装置の糸巻歪補正方
法を得ることを目的とする。
正係数を算出し、糸巻歪を補正して高精度のＸ線画像を
実現することができるＸ線画像診断装置の糸巻歪補正方
法を得ることを目的とする。
【００１１】
【課題を解決するための手段】この発明の請求項１に係
るＸ線画像診断装置の糸巻歪補正方法は、Ｘ線を透過し
にくい材料により被写体に形成されたマーカ中心及びマ
ーカ中心から等距離にある３つ以上のサンプル点を用
い、Ｘ線画像におけるマーカ中心から各サンプル点まで
の距離に基づいてＸ線画像の糸巻歪の歪中心を推定する
ものである。
るＸ線画像診断装置の糸巻歪補正方法は、Ｘ線を透過し
にくい材料により被写体に形成されたマーカ中心及びマ
ーカ中心から等距離にある３つ以上のサンプル点を用
い、Ｘ線画像におけるマーカ中心から各サンプル点まで
の距離に基づいてＸ線画像の糸巻歪の歪中心を推定する
ものである。
【００１２】又、この発明の請求項２に係るＸ線画像診
断装置の糸巻歪補正方法は、Ｘ線を透過しにくい材料に
より被写体に形成されたマーカ中心及びこのマーカ中心
から直交座標軸に沿ってそれぞれ異なる距離の複数のサ
ンプル点を用い、マーカ中心から各サンプル点までの実
際の距離及びＸ線画像上での距離に基づいてＸ線画像の
糸巻歪の歪中心を推定するものである。
断装置の糸巻歪補正方法は、Ｘ線を透過しにくい材料に
より被写体に形成されたマーカ中心及びこのマーカ中心
から直交座標軸に沿ってそれぞれ異なる距離の複数のサ
ンプル点を用い、マーカ中心から各サンプル点までの実
際の距離及びＸ線画像上での距離に基づいてＸ線画像の
糸巻歪の歪中心を推定するものである。
【００１３】又、この発明の請求項５に係るＸ線画像診
断装置の糸巻歪補正方法は、Ｘ線画像の歪中心を求める
ステップと、Ｘ線を透過しにくい材料により被写体に形
成されたマーカ中心及びこのマーカ中心からそれぞれ異
なる距離の複数のサンプル点群を用い、マーカ中心から
各サンプル点までの実際の距離とＸ線画像上での距離と
を比較してサンプル点における歪係数を求めるステップ
と、歪係数に基づいてサンプル点における糸巻歪補正係
数を求めるステップと、Ｘ線画像の歪中心及び糸巻歪補
正係数を用いてＸ線画像を補正するステップとを含むも
のである。
断装置の糸巻歪補正方法は、Ｘ線画像の歪中心を求める
ステップと、Ｘ線を透過しにくい材料により被写体に形
成されたマーカ中心及びこのマーカ中心からそれぞれ異
なる距離の複数のサンプル点群を用い、マーカ中心から
各サンプル点までの実際の距離とＸ線画像上での距離と
を比較してサンプル点における歪係数を求めるステップ
と、歪係数に基づいてサンプル点における糸巻歪補正係
数を求めるステップと、Ｘ線画像の歪中心及び糸巻歪補
正係数を用いてＸ線画像を補正するステップとを含むも
のである。
【００１４】
【作用】この発明においては、既知形状のマーカ中心及
び複数のサンプル点からなる被写体のＸ線画像が元の形
状と一致するように、糸巻歪の歪中心を求める。
び複数のサンプル点からなる被写体のＸ線画像が元の形
状と一致するように、糸巻歪の歪中心を求める。
【００１５】又、この発明においては、既知形状のマー
カ中心及び複数のサンプル点からなる被写体のＸ線画像
の相対位置が元の位置関係と一致するように、各位置毎
に対応させて糸巻歪補正係数を設定し、Ｘ線画像を補正
して糸巻歪を除去する。
カ中心及び複数のサンプル点からなる被写体のＸ線画像
の相対位置が元の位置関係と一致するように、各位置毎
に対応させて糸巻歪補正係数を設定し、Ｘ線画像を補正
して糸巻歪を除去する。
【００１６】
【実施例】以下、この発明の一実施例を図について説明
する。図１は歪中心を求めるためのこの発明の一実施例
を示すフローチャートであり、図２は図１内のステップ
S3を説明するための表示画面を示す平面図である。尚、
図１の歪中心推定ルーチンは、例えばイメージプロセッ
サ４又は他のプロセッサに予めプログラムとして組み込
まれているものとする。又、この発明に使用される装置
は、例えば図５及び図６に示した通りである。
する。図１は歪中心を求めるためのこの発明の一実施例
を示すフローチャートであり、図２は図１内のステップ
S3を説明するための表示画面を示す平面図である。尚、
図１の歪中心推定ルーチンは、例えばイメージプロセッ
サ４又は他のプロセッサに予めプログラムとして組み込
まれているものとする。又、この発明に使用される装置
は、例えば図５及び図６に示した通りである。
【００１７】図２において、5aは表示装置５の表示画面
となる円形のブラウン管、Ｍはブラウン管5a上に表示さ
れた円形マーカ、Ｍoは円形マーカＭのマーカ中心、Ｍ1
〜Ｍ4は円形マーカＭと直交座標軸ＸＹ(マーカ中心Ｍｏ
を原点とする)との交点により決定する４つのサンプル
点である。
となる円形のブラウン管、Ｍはブラウン管5a上に表示さ
れた円形マーカ、Ｍoは円形マーカＭのマーカ中心、Ｍ1
〜Ｍ4は円形マーカＭと直交座標軸ＸＹ(マーカ中心Ｍｏ
を原点とする)との交点により決定する４つのサンプル
点である。
【００１８】次に、図１、図２及び図５〜図７を参照し
ながら、円形マーカＭを用いて歪中心を求めるようにし
たこの発明の一実施例について説明する。まず、Ｘ線Ｒ
を透過しにくい物質(例えば、鉛)で円形マーカＭ及びマ
ーカ中心Ｍｏが形成された板状の被写体３を準備し、Ｘ
線管１及びＸ線像増倍管２の中心軸とマーカ中心Ｍｏと
が一致するように、被写体３をＸ線像増倍管２の前面に
配置する(ステップS1)。
ながら、円形マーカＭを用いて歪中心を求めるようにし
たこの発明の一実施例について説明する。まず、Ｘ線Ｒ
を透過しにくい物質(例えば、鉛)で円形マーカＭ及びマ
ーカ中心Ｍｏが形成された板状の被写体３を準備し、Ｘ
線管１及びＸ線像増倍管２の中心軸とマーカ中心Ｍｏと
が一致するように、被写体３をＸ線像増倍管２の前面に
配置する(ステップS1)。
【００１９】次に、Ｘ線管１からＸ線Ｒを放射して被写
体３に照射し、Ｘ線像増倍管２上に円形マーカＭのＸ線
画像を撮影し(ステップS2)、このＸ線画像をイメージプ
ロセッサ４によりディジタル化して図２のように表示装
置５に表示する。そして、表示装置５のブラウン管5aに
表示された円形マーカＭの線上に、例えば４つのサンプ
ル点Ｍ1〜Ｍ4を指定する(ステップS3)。このとき、マー
カ中心Ｍoも自動的に表示されて指定される。
体３に照射し、Ｘ線像増倍管２上に円形マーカＭのＸ線
画像を撮影し(ステップS2)、このＸ線画像をイメージプ
ロセッサ４によりディジタル化して図２のように表示装
置５に表示する。そして、表示装置５のブラウン管5aに
表示された円形マーカＭの線上に、例えば４つのサンプ
ル点Ｍ1〜Ｍ4を指定する(ステップS3)。このとき、マー
カ中心Ｍoも自動的に表示されて指定される。
【００２０】次に、マーカ中心Ｍｏと各サンプル点Ｍ1
〜Ｍ4との表示画面上での距離を算出する(ステップS
4)。このとき、糸巻歪の大きさは、図７のように、歪中
心に対して距離の関数でほぼ決定するので、歪中心から
離れているサンプル点の距離は長くなる。従って、同一
軸上の２つのサンプル点の各中点、即ち、Ｍ1及びＭ3の
中点と、Ｍ2及びＭ4の中点とを求め、更に、これら中点
間の中点を求めて、これを歪中心として推定する(ステ
ップS5)。
〜Ｍ4との表示画面上での距離を算出する(ステップS
4)。このとき、糸巻歪の大きさは、図７のように、歪中
心に対して距離の関数でほぼ決定するので、歪中心から
離れているサンプル点の距離は長くなる。従って、同一
軸上の２つのサンプル点の各中点、即ち、Ｍ1及びＭ3の
中点と、Ｍ2及びＭ4の中点とを求め、更に、これら中点
間の中点を求めて、これを歪中心として推定する(ステ
ップS5)。
【００２１】次に、マーカ中心Ｍｏから各サンプル点Ｍ
1〜Ｍ4までのＸ線画像上での距離と既知の円形マーカＭ
の実半径とに基づいて、ステップS5で推定された歪中心
がＸ線像増倍管２の光電陰極21上のどの位置に相当する
かを算出する。そして、光電陰極21上に配置された実際
のマーカ中心Ｍｏの位置と、推定された歪中心の光電陰
極21上での位置とを比較し、各中心位置の差異が所定の
許容誤差範囲内であるか否かを判定する(ステップS6)。
1〜Ｍ4までのＸ線画像上での距離と既知の円形マーカＭ
の実半径とに基づいて、ステップS5で推定された歪中心
がＸ線像増倍管２の光電陰極21上のどの位置に相当する
かを算出する。そして、光電陰極21上に配置された実際
のマーカ中心Ｍｏの位置と、推定された歪中心の光電陰
極21上での位置とを比較し、各中心位置の差異が所定の
許容誤差範囲内であるか否かを判定する(ステップS6)。
【００２２】もし、実際のマーカ中心Ｍｏの位置が、光
電陰極21上での推定された歪中心位置から許容誤差範囲
以上離れていれば、ステップS5で推定した歪中心位置の
信頼性が低いと見なされる。従って、円形マーカＭを有
する被写体３を光電陰極21上で移動させて、マーカ中心
が推定された歪中心と一致するように再度位置決めし
(ステップS7)、ステップS2〜S6を繰り返す。
電陰極21上での推定された歪中心位置から許容誤差範囲
以上離れていれば、ステップS5で推定した歪中心位置の
信頼性が低いと見なされる。従って、円形マーカＭを有
する被写体３を光電陰極21上で移動させて、マーカ中心
が推定された歪中心と一致するように再度位置決めし
(ステップS7)、ステップS2〜S6を繰り返す。
【００２３】そして、ステップS6において、実際のマー
カ中心位置と推定された歪中心位置との差異が許容誤差
範囲内と判定された時点で、推定された歪中心を歪中心
と決定し、この歪中心位置を光電陰極21上に印を付け、
歪中心推定ルーチン(図１)を終了する。
カ中心位置と推定された歪中心位置との差異が許容誤差
範囲内と判定された時点で、推定された歪中心を歪中心
と決定し、この歪中心位置を光電陰極21上に印を付け、
歪中心推定ルーチン(図１)を終了する。
【００２４】その後、人体等の被写体３を診断するとき
は、光電陰極21上の印を目標にしてＸ線管１及びＸ線像
増倍管２の軸合わせを行うことができる。又、ブラウン
管5a上の歪中心位置をイメージプロセッサ４に記憶させ
ておくことができる。こうして、歪中心が確定すること
により、例えば中間部電極25(図６参照)を用いた糸巻歪
の補正が対称に行われるので、歪補正の信頼性が向上す
る。
は、光電陰極21上の印を目標にしてＸ線管１及びＸ線像
増倍管２の軸合わせを行うことができる。又、ブラウン
管5a上の歪中心位置をイメージプロセッサ４に記憶させ
ておくことができる。こうして、歪中心が確定すること
により、例えば中間部電極25(図６参照)を用いた糸巻歪
の補正が対称に行われるので、歪補正の信頼性が向上す
る。
【００２５】尚、上記実施例では、ステップS5におい
て、２つのサンプル点間の中点を歪中心として推定した
が、同一軸上の２つのサンプル点のマーカ中心Ｍｏから
の距離に逆比例した点をＸ軸及びＹ軸上でそれぞれ求
め、各々の座標をもつ点を歪中心として推定してもよ
い。例えば、マーカ中心ＭｏからＸ軸上のサンプル点Ｍ
1及びＭ3までの距離の比が２：１であったとすると、サ
ンプル点Ｍ1及びＭ2から距離の比が１：２となる点をＸ
軸上の歪中心と推定する。
て、２つのサンプル点間の中点を歪中心として推定した
が、同一軸上の２つのサンプル点のマーカ中心Ｍｏから
の距離に逆比例した点をＸ軸及びＹ軸上でそれぞれ求
め、各々の座標をもつ点を歪中心として推定してもよ
い。例えば、マーカ中心ＭｏからＸ軸上のサンプル点Ｍ
1及びＭ3までの距離の比が２：１であったとすると、サ
ンプル点Ｍ1及びＭ2から距離の比が１：２となる点をＸ
軸上の歪中心と推定する。
【００２６】又、４つのサンプル点Ｍ1〜Ｍ4の距離のバ
ラツキが最小となる位置を最小２乗法で推定してもよ
く、４点Ｍ1〜Ｍ4の距離のバラツキが最小となる位置を
シンプレックス法で推定してもよい。
ラツキが最小となる位置を最小２乗法で推定してもよ
く、４点Ｍ1〜Ｍ4の距離のバラツキが最小となる位置を
シンプレックス法で推定してもよい。
【００２７】又、サンプル点を４点としたが、３点以上
設定すればよい。又、Ｘ線画像上のサンプル点の距離を
歪中心の推定に用いたが、ブラウン管5aでのピクセル値
を用いてもよい。又、ブラウン管5a上での円形マーカＭ
の目視により、真円となる位置を試行錯誤で求めてもよ
い。
設定すればよい。又、Ｘ線画像上のサンプル点の距離を
歪中心の推定に用いたが、ブラウン管5aでのピクセル値
を用いてもよい。又、ブラウン管5a上での円形マーカＭ
の目視により、真円となる位置を試行錯誤で求めてもよ
い。
【００２８】又、円形マーカの代わりに例えばグリッド
マーカを用いることもでき、この場合、ブラウン管5a上
で対称なグリッド図形が再現される位置を歪中心と推定
することができる。更に、グリッドマーカを用いた場
合、以下のように、画像位置に対応した糸巻歪の補正係
数を演算により求めることができる。
マーカを用いることもでき、この場合、ブラウン管5a上
で対称なグリッド図形が再現される位置を歪中心と推定
することができる。更に、グリッドマーカを用いた場
合、以下のように、画像位置に対応した糸巻歪の補正係
数を演算により求めることができる。
【００２９】図３はグリッドマーカが形成された被写体
３を用いて糸巻歪を補正するための演算処理ルーチンを
示すフローチャート、図４は表示装置５のブラウン管5a
上に表示されたグリッドマーカ及びその格子点を示す平
面図である。尚、図３の歪補正ルーチンは、例えばイメ
ージプロセッサ４又は他のプロセッサに予めプログラム
として組み込まれているものとする。
３を用いて糸巻歪を補正するための演算処理ルーチンを
示すフローチャート、図４は表示装置５のブラウン管5a
上に表示されたグリッドマーカ及びその格子点を示す平
面図である。尚、図３の歪補正ルーチンは、例えばイメ
ージプロセッサ４又は他のプロセッサに予めプログラム
として組み込まれているものとする。
【００３０】又、図４において、Ｍｏはグリッドマーカ
の中心となる格子点即ちマーカ中心であり、マーカ中心
Ｍｏを原点とする各直交グリッド線は直交座標軸Ｘ及び
Ｙに対応している。従って、この場合、直交座標系の回
転方向の表示状態も確定される。又、各グリッドマーカ
線の交点からなる格子点Ｍ11〜Ｍ48は、実際には等間隔
のマトリクス状に配列されているものとする。
の中心となる格子点即ちマーカ中心であり、マーカ中心
Ｍｏを原点とする各直交グリッド線は直交座標軸Ｘ及び
Ｙに対応している。従って、この場合、直交座標系の回
転方向の表示状態も確定される。又、各グリッドマーカ
線の交点からなる格子点Ｍ11〜Ｍ48は、実際には等間隔
のマトリクス状に配列されているものとする。
【００３１】次に、図３〜図７を参照しながら、糸巻歪
を補正したＸ線画像が取得可能にしたこの発明の他の実
施例について説明する。この場合、歪中心は前述のルー
チンにより予め決定され且つ光電陰極21上に印が付され
ており、Ｘ線中心軸が歪中心と一致するように位置合わ
せが行われているものとする。
を補正したＸ線画像が取得可能にしたこの発明の他の実
施例について説明する。この場合、歪中心は前述のルー
チンにより予め決定され且つ光電陰極21上に印が付され
ており、Ｘ線中心軸が歪中心と一致するように位置合わ
せが行われているものとする。
【００３２】まず、Ｘ線Ｒを透過しにくい物質でグリッ
ドマーカが形成された板状の被写体３を準備し、マーカ
中心ＭｏがＸ線管１及びＸ線像増倍管２の軸上に一致す
るように、且つ、マーカ中心Ｍｏが歪中心位置とほぼ一
致するように被写体３をＸ線像増倍管２上に設置する
(ステップＳ11)。このとき、マーカ中心Ｍｏと歪中心と
が必ずしも一致しなくてもよい。
ドマーカが形成された板状の被写体３を準備し、マーカ
中心ＭｏがＸ線管１及びＸ線像増倍管２の軸上に一致す
るように、且つ、マーカ中心Ｍｏが歪中心位置とほぼ一
致するように被写体３をＸ線像増倍管２上に設置する
(ステップＳ11)。このとき、マーカ中心Ｍｏと歪中心と
が必ずしも一致しなくてもよい。
【００３３】次に、Ｘ線管１からＸ線Ｒを放射して被写
体３に照射し、Ｘ線像増倍管２上にグリッドマーカのＸ
線画像を撮影し(ステップＳ12)、このＸ線画像をイメー
ジプロセッサ４によりディジタル化して図４のように表
示装置５に表示する。このとき、マーカ中心Ｍｏと歪中
心とが既に位置合わせされているので、糸巻歪は図示し
たようにマーカ中心Ｍｏから対称に生じる。
体３に照射し、Ｘ線像増倍管２上にグリッドマーカのＸ
線画像を撮影し(ステップＳ12)、このＸ線画像をイメー
ジプロセッサ４によりディジタル化して図４のように表
示装置５に表示する。このとき、マーカ中心Ｍｏと歪中
心とが既に位置合わせされているので、糸巻歪は図示し
たようにマーカ中心Ｍｏから対称に生じる。
【００３４】次に、表示装置５のブラウン管5aに表示さ
れたグリッドマーカ上で、マーカ中心Ｍoを指定すると
共に、マーカ中心Ｍoから等距離の格子点(例えば、Ｍ11〜
Ｍ14、Ｍ21〜Ｍ24、Ｍ31〜Ｍ34及びＭ41〜Ｍ48)を、それ
ぞれ等距離のサンプル点群として指定する(ステップS1
3)。そして、マーカ中心Ｍoと各サンプル点Ｍ11〜Ｍ48
との表示画面上での距離を、例えばピクセル数により算
出する(ステップS14)。
れたグリッドマーカ上で、マーカ中心Ｍoを指定すると
共に、マーカ中心Ｍoから等距離の格子点(例えば、Ｍ11〜
Ｍ14、Ｍ21〜Ｍ24、Ｍ31〜Ｍ34及びＭ41〜Ｍ48)を、それ
ぞれ等距離のサンプル点群として指定する(ステップS1
3)。そして、マーカ中心Ｍoと各サンプル点Ｍ11〜Ｍ48
との表示画面上での距離を、例えばピクセル数により算
出する(ステップS14)。
【００３５】このとき、マーカ中心Ｍｏから各サンプル
点Ｍ11〜Ｍ48までの実際の距離は既知である。又、糸巻
歪の大きさは、図７のように、歪中心に対して距離の関
数でほぼ決定するので、マーカ中心Ｍｏから等距離のサ
ンプル点群の各位置は、ブラウン管上においても等距離
として算出される。従って、マーカ中心Ｍｏから各サン
プル点Ｍ11〜Ｍ48までのブラウン管5a上の距離と実際の
距離とに基づいて、各サンプル点の位置での糸巻歪の大
きさを表わす歪係数を算出することができる(ステップS
15)。
点Ｍ11〜Ｍ48までの実際の距離は既知である。又、糸巻
歪の大きさは、図７のように、歪中心に対して距離の関
数でほぼ決定するので、マーカ中心Ｍｏから等距離のサ
ンプル点群の各位置は、ブラウン管上においても等距離
として算出される。従って、マーカ中心Ｍｏから各サン
プル点Ｍ11〜Ｍ48までのブラウン管5a上の距離と実際の
距離とに基づいて、各サンプル点の位置での糸巻歪の大
きさを表わす歪係数を算出することができる(ステップS
15)。
【００３６】次に、ステップS15で算出された歪係数の
逆数に基づいて、マーカ中心Ｍo及び各サンプル点Ｍ11
〜Ｍ48の表示位置を実際の位置と一致させるための糸巻
歪補正係数を算出し(ステップS16)、この糸巻歪補正係
数に基づいて、グリッドマーカのＸ線画像を補正する(ス
テップS17)
逆数に基づいて、マーカ中心Ｍo及び各サンプル点Ｍ11
〜Ｍ48の表示位置を実際の位置と一致させるための糸巻
歪補正係数を算出し(ステップS16)、この糸巻歪補正係
数に基づいて、グリッドマーカのＸ線画像を補正する(ス
テップS17)
【００３７】最後に、補正後のＸ線画像が良好であるか
否かを表示装置５上に表示して判定し(ステップS18)、
もし良好でなければ、ステップS16に戻って糸巻歪補正
係数の算出を再度行い、Ｘ線画像が元のグリッドマーカ
と同一となるまで繰り返す。又、補正後のＸ線画像が良
好であることが確認されれば、その時点の糸巻歪補正係
数をデータとして格納し(ステップＳ19)、図３の糸巻歪
補正ルーチンを終了する。この結果、高精度のＸ線画像
が得られ、被写体３に対して信頼性の高い診断を行うこ
とができる。
否かを表示装置５上に表示して判定し(ステップS18)、
もし良好でなければ、ステップS16に戻って糸巻歪補正
係数の算出を再度行い、Ｘ線画像が元のグリッドマーカ
と同一となるまで繰り返す。又、補正後のＸ線画像が良
好であることが確認されれば、その時点の糸巻歪補正係
数をデータとして格納し(ステップＳ19)、図３の糸巻歪
補正ルーチンを終了する。この結果、高精度のＸ線画像
が得られ、被写体３に対して信頼性の高い診断を行うこ
とができる。
【００３８】尚、各サンプル点間のピクセル位置での糸
巻歪補正係数の算出法としては、マーカ中心Ｍｏから各
サンプル点までの実際の距離を横軸ｘとし、ブラウン管
5a上での距離を縦軸ｙとしたときに、各算出距離をプロ
ットした点を結ぶ曲線の関数の次数ｎから求めることが
できる。
巻歪補正係数の算出法としては、マーカ中心Ｍｏから各
サンプル点までの実際の距離を横軸ｘとし、ブラウン管
5a上での距離を縦軸ｙとしたときに、各算出距離をプロ
ットした点を結ぶ曲線の関数の次数ｎから求めることが
できる。
【００３９】即ち、補正関数をｘ＝ｆ(ｙ)とすれば、各
点を通る関数はｙ＝ｆ-1(ｘ)で表わされ、マーカ中心Ｍ
ｏが歪中心と一致している場合、関数はｙ＝aｘnとな
り、補正関数はｘ＝(ｘ／ａ)-nで表わされる。この場
合、糸巻歪は原点を通る高次式により決定され、次数ｎ
の値は、各サンプル点を通る曲線において、各点との誤
差が最小となるように最小２乗法で求めることができ
る。
点を通る関数はｙ＝ｆ-1(ｘ)で表わされ、マーカ中心Ｍ
ｏが歪中心と一致している場合、関数はｙ＝aｘnとな
り、補正関数はｘ＝(ｘ／ａ)-nで表わされる。この場
合、糸巻歪は原点を通る高次式により決定され、次数ｎ
の値は、各サンプル点を通る曲線において、各点との誤
差が最小となるように最小２乗法で求めることができ
る。
【００４０】又、ステップS11においてマーカ中心Ｍoが
歪中心と一致していない場合、関数がｎ次の多項式ｙ＝
anｘn＋an-1ｘn-1＋…＋a1ｘ＋a0で表わされるが、各係
数ａ0〜ａnは、同様に最小２乗法を用いて求めることが
できる。但し、次数ｎは、マーカ中心Ｍｏからの距離が
異なる各サンプル点群の数（図４の場合、Ｍ11〜Ｍ14、
Ｍ21〜Ｍ24、Ｍ31〜Ｍ34及びＭ41〜Ｍ48の４組）以下に
設定されることは言うまでもない。
歪中心と一致していない場合、関数がｎ次の多項式ｙ＝
anｘn＋an-1ｘn-1＋…＋a1ｘ＋a0で表わされるが、各係
数ａ0〜ａnは、同様に最小２乗法を用いて求めることが
できる。但し、次数ｎは、マーカ中心Ｍｏからの距離が
異なる各サンプル点群の数（図４の場合、Ｍ11〜Ｍ14、
Ｍ21〜Ｍ24、Ｍ31〜Ｍ34及びＭ41〜Ｍ48の４組）以下に
設定されることは言うまでもない。
【００４１】尚、最小２乗法のみならず、シンプレック
ス法によりｎ次の多項式を推定してもよい。又、マーカ
中心Ｍｏから各サンプル点Ｍ11〜Ｍ48までの距離を、ピ
クセル数の代わりに表示画面上の距離を用いてもよい。
又、補正係数を求めるためにグリッドマーカを用いた
が、マーカ中心Ｍｏからの距離が既知のサンプル点であ
れば何でもよく、例えば複数の円形マーカと直交座標軸
との交点を用いてもよい。
ス法によりｎ次の多項式を推定してもよい。又、マーカ
中心Ｍｏから各サンプル点Ｍ11〜Ｍ48までの距離を、ピ
クセル数の代わりに表示画面上の距離を用いてもよい。
又、補正係数を求めるためにグリッドマーカを用いた
が、マーカ中心Ｍｏからの距離が既知のサンプル点であ
れば何でもよく、例えば複数の円形マーカと直交座標軸
との交点を用いてもよい。
【００４２】更に、歪中心を予め求めておき、ステップ
S11において歪中心にマーカ中心Ｍoを位置合わせする場
合について説明したが、歪中心が予め推定されていない
場合でも、グリッドマーカを用いて歪中心を推定するこ
とができる。即ち、図４のように表示画面上でマーカ中
心Ｍｏを原点とする直交座標系ＸＹを考え、マーカ中心
ＭｏからＸ軸上の各サンプル点Ｍ12、Ｍ14、Ｍ32及びＭ
34までの距離、並びに、マーカ中心ＭｏからＹ軸上の各
サンプル点Ｍ11、Ｍ13、Ｍ31及びＭ33までの距離を求め
る。
S11において歪中心にマーカ中心Ｍoを位置合わせする場
合について説明したが、歪中心が予め推定されていない
場合でも、グリッドマーカを用いて歪中心を推定するこ
とができる。即ち、図４のように表示画面上でマーカ中
心Ｍｏを原点とする直交座標系ＸＹを考え、マーカ中心
ＭｏからＸ軸上の各サンプル点Ｍ12、Ｍ14、Ｍ32及びＭ
34までの距離、並びに、マーカ中心ＭｏからＹ軸上の各
サンプル点Ｍ11、Ｍ13、Ｍ31及びＭ33までの距離を求め
る。
【００４３】そして、Ｘ軸上の各サンプル点に対して、
実際の距離ｐと画面上の距離ｑとの関係を、ｎ次の多項
式ｑ＝ｇ(ｐ)＝ｂnｐn＋ｂn-1ｐn-1＋…＋ｂ1ｐ＋ｂ0で
近似する。この多項式の各係数ｂ0〜ｂnは、例えば最小
２乗法又はシンプレックス法で決定される。又、Ｙ軸上
のサンプル点に対しても同様の多項式が決定される。
実際の距離ｐと画面上の距離ｑとの関係を、ｎ次の多項
式ｑ＝ｇ(ｐ)＝ｂnｐn＋ｂn-1ｐn-1＋…＋ｂ1ｐ＋ｂ0で
近似する。この多項式の各係数ｂ0〜ｂnは、例えば最小
２乗法又はシンプレックス法で決定される。又、Ｙ軸上
のサンプル点に対しても同様の多項式が決定される。
【００４４】次に、マーカ中心Ｍoを歪中心と一致させ
るための座標値ｐ0を考え、ｎ次の関数ｑ＝ｇ(ｐ)＝ｂn
(ｐ−ｐ0)nを満たすｐ0を最小２乗法又はシンプレック
ス法で決定する。これにより、Ｘ軸及びＹ軸に対して決
定されるｐ0の値をそれぞれＸ0及びＹ0とすれば、歪中
心の座標(Ｘ0,Ｙ0)を推定することができる。以下、前
述と同様に、推定された歪中心に基づいて全ての格子点
に対する補正係数を求めることができる。
るための座標値ｐ0を考え、ｎ次の関数ｑ＝ｇ(ｐ)＝ｂn
(ｐ−ｐ0)nを満たすｐ0を最小２乗法又はシンプレック
ス法で決定する。これにより、Ｘ軸及びＹ軸に対して決
定されるｐ0の値をそれぞれＸ0及びＹ0とすれば、歪中
心の座標(Ｘ0,Ｙ0)を推定することができる。以下、前
述と同様に、推定された歪中心に基づいて全ての格子点
に対する補正係数を求めることができる。
【００４５】
【発明の効果】以上のようにこの発明によれば、Ｘ線を
透過しにくい材料により被写体に形成されたマーカ中心
とマーカ中心から等距離にある３つ以上のサンプル点又
は直交座標軸に沿ってそれぞれ異なる距離の複数のサン
プル点とを用い、Ｘ線画像におけるマーカ中心から各サ
ンプル点までの距離又は実際の距離及びＸ線画像上での
距離に基づいてＸ線画像の糸巻歪の歪中心を推定し、既
知形状のマーカ中心及び複数のサンプル点からなる被写
体のＸ線画像が元の形状と一致するように歪中心を求め
るようにしたので、糸巻歪の補正を容易にしたＸ線画像
診断装置の糸巻歪補正方法が得られる効果がある。
透過しにくい材料により被写体に形成されたマーカ中心
とマーカ中心から等距離にある３つ以上のサンプル点又
は直交座標軸に沿ってそれぞれ異なる距離の複数のサン
プル点とを用い、Ｘ線画像におけるマーカ中心から各サ
ンプル点までの距離又は実際の距離及びＸ線画像上での
距離に基づいてＸ線画像の糸巻歪の歪中心を推定し、既
知形状のマーカ中心及び複数のサンプル点からなる被写
体のＸ線画像が元の形状と一致するように歪中心を求め
るようにしたので、糸巻歪の補正を容易にしたＸ線画像
診断装置の糸巻歪補正方法が得られる効果がある。
【００４６】又、この発明によれば、Ｘ線画像の歪中心
を求めるステップと、Ｘ線を透過しにくい材料により被
写体に形成されたマーカ中心及びこのマーカ中心からそ
れぞれ異なる距離の複数のサンプル点群を用い、マーカ
中心から各サンプル点までの実際の距離とＸ線画像上で
の距離とを比較してサンプル点における歪係数を求める
ステップと、歪係数に基づいてサンプル点における糸巻
歪補正係数を求めるステップと、Ｘ線画像の歪中心及び
糸巻歪補正係数を用いてＸ線画像を補正するステップと
を含み、既知形状のマーカ中心及び複数のサンプル点か
らなる被写体のＸ線画像の相対位置が元の位置関係と一
致するように各位置毎に対応させて糸巻歪補正係数を設
定するようにしたので、簡単な補正係数を用いて糸巻歪
を有効に補正し、高精度のＸ線画像を取得できるＸ線画
像診断装置の糸巻歪補正方法が得られる効果がある。
を求めるステップと、Ｘ線を透過しにくい材料により被
写体に形成されたマーカ中心及びこのマーカ中心からそ
れぞれ異なる距離の複数のサンプル点群を用い、マーカ
中心から各サンプル点までの実際の距離とＸ線画像上で
の距離とを比較してサンプル点における歪係数を求める
ステップと、歪係数に基づいてサンプル点における糸巻
歪補正係数を求めるステップと、Ｘ線画像の歪中心及び
糸巻歪補正係数を用いてＸ線画像を補正するステップと
を含み、既知形状のマーカ中心及び複数のサンプル点か
らなる被写体のＸ線画像の相対位置が元の位置関係と一
致するように各位置毎に対応させて糸巻歪補正係数を設
定するようにしたので、簡単な補正係数を用いて糸巻歪
を有効に補正し、高精度のＸ線画像を取得できるＸ線画
像診断装置の糸巻歪補正方法が得られる効果がある。
【図１】歪中心を求めるためのこの発明の一実施例を示
すフローチャートである。
すフローチャートである。
【図２】図１内のステップS3を説明するためのブラウン
管画面を示す平面図である。
管画面を示す平面図である。
【図３】糸巻歪を補正するためのこの発明の他の実施例
を示すフローチャートである。
を示すフローチャートである。
【図４】図２内のステップS13を説明するためのブラウ
ン管画面を示す平面図である。
ン管画面を示す平面図である。
【図５】一般的なＸ線画像診断装置を示す構成図であ
る。
る。
【図６】図５内のＸ線像増倍管を詳細に示す断面図であ
る。
る。
【図７】図６のＸ線像増倍管の糸巻歪を説明するための
特性図である。
特性図である。
１ Ｘ線管 ２ Ｘ線像増倍管 ３ 被写体 ５ 表示装置 Ｒ Ｘ線 Ｍｏ マーカ中心 Ｍ1〜Ｍ4 サンプル点 Ｍ11〜Ｍ48 サンプル点群 Ｘ、Ｙ 直交座標軸 S4、S14 マーカ中心からサンプル点までの距離を算
出するステップ S5 糸巻歪の歪中心を推定するステップ S6 マーカ中心と推定歪中心との差異が誤差範囲内か
を判定するステップ S11 マーカ中心を歪中心と一致させるステップ S15 サンプル点における歪係数を求めるステップ S16 歪係数に基づいてサンプル点での糸巻歪補正係
数を求めるステップ S17 糸巻歪補正係数を用いてＸ線画像を補正するス
テップ
出するステップ S5 糸巻歪の歪中心を推定するステップ S6 マーカ中心と推定歪中心との差異が誤差範囲内か
を判定するステップ S11 マーカ中心を歪中心と一致させるステップ S15 サンプル点における歪係数を求めるステップ S16 歪係数に基づいてサンプル点での糸巻歪補正係
数を求めるステップ S17 糸巻歪補正係数を用いてＸ線画像を補正するス
テップ
Claims (9)
- 【請求項１】 被写体を透過したＸ線によるＸ線像増倍
管からのＸ線画像を用いたＸ線画像診断装置の糸巻歪補
正方法において、 前記Ｘ線を透過しにくい材料により前記被写体に形成さ
れたマーカ中心及びこのマーカ中心から等距離にある３
つ以上のサンプル点を用い、 前記Ｘ線画像における前記マーカ中心から前記各サンプ
ル点までの距離に基づいて前記Ｘ線画像の糸巻歪の歪中
心を推定することを特徴とするＸ線画像診断装置の糸巻
歪補正方法。 - 【請求項２】 被写体を透過したＸ線によるＸ線像増倍
管からのＸ線画像を用いたＸ線画像診断装置の糸巻歪補
正方法において、 前記Ｘ線を透過しにくい材料により前記被写体に形成さ
れたマーカ中心及びこのマーカ中心から直交座標軸に沿
ってそれぞれ異なる距離の複数のサンプル点を用い、 前記マーカ中心から前記各サンプル点までの実際の距離
及び前記Ｘ線画像上での距離に基づいて前記Ｘ線画像の
糸巻歪の歪中心を推定することを特徴とするＸ線画像診
断装置の糸巻歪補正方法。 - 【請求項３】 前記マーカ中心の位置と前記推定された
歪中心の位置との差異が許容誤差範囲内であることを判
定したときに前記推定された歪中心を実際の歪中心とし
て決定することを特徴とする請求項１又は請求項２のＸ
線画像診断装置の糸巻歪補正方法。 - 【請求項４】 前記歪中心の前記Ｘ線像増倍管上での位
置に印を付けるステップと、前記Ｘ線の中心軸を前記印
に一致させるステップとを含むことを特徴とする請求項
１乃至請求項３のいずれかのＸ線画像診断装置の糸巻歪
補正方法。 - 【請求項５】 被写体を透過したＸ線によるＸ線像増倍
管からのＸ線画像を用いたＸ線画像診断装置の糸巻歪補
正方法において、 前記Ｘ線画像の歪中心を求めるステップと、 Ｘ線を透過しにくい材料により前記被写体に形成された
マーカ中心及びこのマーカ中心からそれぞれ異なる距離
の複数のサンプル点群を用い、 前記マーカ中心から前記各サンプル点までの実際の距離
と前記Ｘ線画像上での距離とを比較して前記サンプル点
における歪係数を求めるステップと、 前記歪係数に基づいて前記サンプル点における糸巻歪補
正係数を求めるステップと、前記Ｘ線画像の歪中心及び 前記糸巻歪補正係数を用いて
前記Ｘ線画像を補正するステップと、 を含むことを特徴とするＸ線画像診断装置の糸巻歪補正
方法。 - 【請求項６】 前記歪係数を求めるステップの前に、前
記マーカ中心を前記歪中心と一致させるステップを含む
ことを特徴とする請求項５のＸ線画像診断装置の糸巻歪
補正方法。 - 【請求項７】 前記マーカ中心から前記サンプル点まで
の間の実際の距離と前記Ｘ線画像上の距離との関係をｎ
次の多項式で近似することにより、前記サンプル点の各
々の間のＸ線画像位置での糸巻歪補正係数を求めること
を特徴とする請求項５又は請求項６のＸ線画像診断装置
の糸巻歪補正方法。 - 【請求項８】 マトリクス状の複数の格子点を含むグリ
ッドマーカを前記被写体に形成し、前記格子点を前記サ
ンプル点としたことを特徴とする請求項１乃至請求項７
のいずれかのＸ線画像診断装置の糸巻歪補正方法。 - 【請求項９】 前記マーカ中心と前記サンプル点との間
の距離を、前記Ｘ線画像上のピクセル数により求めるこ
とを特徴とする請求項１乃至請求項８のいずれかのＸ線
画像診断装置の糸巻歪補正方法。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP3084956A JP2791231B2 (ja) | 1991-04-17 | 1991-04-17 | Ｘ線画像診断装置の糸巻歪補正方法 |
| US07/869,506 US5263074A (en) | 1991-04-17 | 1992-04-16 | Spool distortion correction method for an x-ray radiograph diagnosis device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP3084956A JP2791231B2 (ja) | 1991-04-17 | 1991-04-17 | Ｘ線画像診断装置の糸巻歪補正方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH04319236A JPH04319236A (ja) | 1992-11-10 |
| JP2791231B2 true JP2791231B2 (ja) | 1998-08-27 |
Family
ID=13845087
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP3084956A Expired - Fee Related JP2791231B2 (ja) | 1991-04-17 | 1991-04-17 | Ｘ線画像診断装置の糸巻歪補正方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US5263074A (ja) |
| JP (1) | JP2791231B2 (ja) |
Families Citing this family (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5600701A (en) * | 1993-05-11 | 1997-02-04 | Hitachi Medical Corporation | X-ray imaging system and method therefor |
| DE4316847A1 (de) * | 1993-05-19 | 1994-11-24 | Philips Patentverwaltung | Verfahren zum Entzerren von Röntgenaufnahmen und Anordnung zur Durchführung des Verfahrens |
| US5526442A (en) * | 1993-10-04 | 1996-06-11 | Hitachi Medical Corporation | X-ray radiography method and system |
| US5730129A (en) * | 1995-04-03 | 1998-03-24 | General Electric Company | Imaging of interventional devices in a non-stationary subject |
| US5642395A (en) * | 1995-08-07 | 1997-06-24 | Oec Medical Systems, Inc. | Imaging chain with miniaturized C-arm assembly for mobile X-ray imaging system |
| US6282261B1 (en) * | 1996-02-21 | 2001-08-28 | Lunar Corporation | Multi-mode x-ray image intensifier system |
| EP0882275A1 (en) * | 1996-02-21 | 1998-12-09 | Lunar Corporation | X-ray imaging system |
| US5757880A (en) * | 1997-01-08 | 1998-05-26 | Colomb; Denis | Apparatus, article of manufacture, and method for creation of an uncompressed image of compressed matter |
| NL1005496C2 (nl) * | 1997-03-11 | 1998-09-14 | Optische Ind Oede Oude Delftoe | Werkwijze voor het corrigeren van beeldvertekening in een röntgenbeeld. |
| US6260999B1 (en) * | 1999-07-26 | 2001-07-17 | Siemens Medical Systems, Inc. | Isocenter localization using electronic portal imaging |
| FR2823968B1 (fr) * | 2001-04-27 | 2005-01-14 | Ge Med Sys Global Tech Co Llc | Procede d'etalonnage d'un systeme d'imagerie, support de memoire et dispositif associe |
| JP2004234379A (ja) * | 2003-01-30 | 2004-08-19 | Sony Corp | 画像処理方法、画像処理装置及び画像処理方法を適用した撮像装置、表示装置 |
| US8484001B2 (en) * | 2003-08-26 | 2013-07-09 | Voyant Health Ltd. | Pre-operative medical planning system and method for use thereof |
| US7564999B2 (en) * | 2005-07-25 | 2009-07-21 | Carestream Health, Inc. | Method for identifying markers in radiographic images |
| JP4854314B2 (ja) * | 2006-01-27 | 2012-01-18 | キヤノン株式会社 | 情報処理装置及びその制御方法、プログラム |
| JP5065779B2 (ja) * | 2007-06-27 | 2012-11-07 | 株式会社日立製作所 | Ｘ線荷物検査装置の評価方法 |
| JP5239585B2 (ja) * | 2008-07-28 | 2013-07-17 | 株式会社島津製作所 | Ｘ線撮像装置 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4048507A (en) * | 1976-02-11 | 1977-09-13 | Gaston Alexis Neal De | X-ray beam perpendicular finder |
| JPS63229419A (ja) * | 1987-03-18 | 1988-09-26 | Komatsu Ltd | レンズ歪補正装置 |
| JPH0210636A (ja) * | 1988-06-28 | 1990-01-16 | Toshiba Corp | Ｘ線画像診断装置 |
-
1991
- 1991-04-17 JP JP3084956A patent/JP2791231B2/ja not_active Expired - Fee Related
-
1992
- 1992-04-16 US US07/869,506 patent/US5263074A/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH04319236A (ja) | 1992-11-10 |
| US5263074A (en) | 1993-11-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2791231B2 (ja) | Ｘ線画像診断装置の糸巻歪補正方法 | |
| Chakraborty | Image intensifier distortion correction | |
| EP0363907B1 (en) | System for correcting extension rate of images | |
| EP1672672A2 (en) | Charged particle beam apparatus, method of displaying sample image, and method of measuring image shift sensitivity | |
| JPS5825042A (ja) | 陰極線管 | |
| US4173720A (en) | Method and apparatus for image construction | |
| JP3041202B2 (ja) | デジタルコンバージェンス補正方法およびその装置 | |
| GB1597170A (en) | Apparatus for making laminar radiograms | |
| JPH0982257A (ja) | 荷電粒子光学鏡筒における非点収差の補正及び焦点合わせ方法 | |
| Reimann et al. | Automated distortion correction of x-ray image intensifier images | |
| Cosby et al. | Computer-aided radiation therapy simulation: image intensifier spatial distortion correction for large field of view digital fluoroscopy | |
| JP3728328B2 (ja) | カラー表示装置 | |
| EP0974120A1 (en) | Method for correcting image distortions in an x-ray image | |
| JPH02203681A (ja) | Ｘ線ｔｖ装置 | |
| Wang et al. | Accurate and efficient image intensifier distortion correction algorithm for volume tomographic angiography | |
| JPH01241735A (ja) | コンバーゼンス測定装置 | |
| JP3344329B2 (ja) | ランディング調整方法及びその装置 | |
| JPH10215467A (ja) | 位置測定方法及び位置測定システム | |
| KR20020021954A (ko) | 평판 디지털 엑스-레이 검출기의 특성 오류 정정 방법 | |
| JP3463654B2 (ja) | 陰極線管 | |
| JP4127901B2 (ja) | 放射画像における対称線の決定方法 | |
| JPH0528316U (ja) | Ｘ線画像診断装置 | |
| JPH01313895A (ja) | 放射線照射装置 | |
| KR0123915B1 (ko) | 페라이트판 부착위치 결정방법 | |
| JPH01267783A (ja) | 放射線画像処理装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080612 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080612 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090612 Year of fee payment: 11 |
|
| LAPS | Cancellation because of no payment of annual fees |