JPH04233377A - ディジタル影像装置の変調伝達関数決定方法 - Google Patents
ディジタル影像装置の変調伝達関数決定方法Info
- Publication number
- JPH04233377A JPH04233377A JP3131414A JP13141491A JPH04233377A JP H04233377 A JPH04233377 A JP H04233377A JP 3131414 A JP3131414 A JP 3131414A JP 13141491 A JP13141491 A JP 13141491A JP H04233377 A JPH04233377 A JP H04233377A
- Authority
- JP
- Japan
- Prior art keywords
- detected signal
- transfer function
- determined
- input screen
- value
- 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
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/26—Measuring, controlling or protecting
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N17/00—Diagnosis, testing or measuring for television systems or their details
- H04N17/002—Diagnosis, testing or measuring for television systems or their details for television cameras
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N23/00—Cameras or camera modules comprising electronic image sensors; Control thereof
- H04N23/80—Camera processing pipelines; Components thereof
- H04N23/81—Camera processing pipelines; Components thereof for suppressing or minimising disturbance in the image signal generation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N5/00—Details of television systems
- H04N5/30—Transforming light or analogous information into electric information
- H04N5/32—Transforming X-rays
Landscapes
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Toxicology (AREA)
- Biomedical Technology (AREA)
- Image Analysis (AREA)
- Closed-Circuit Television Systems (AREA)
- Radiography Using Non-Light Waves (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Transforming Light Signals Into Electric Signals (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【０００１】
【産業上の利用分野】本発明は、試験対象により空間変
調された放射強度分布が、入力スクリーンの画素に対応
するメモリー位置における検出信号値を記憶するためデ
ィジタルデータ処理装置に接続された検出器の入力スク
リーンに表示されるディジタル影像装置の変調伝達関数
を決定する方法に係る。
調された放射強度分布が、入力スクリーンの画素に対応
するメモリー位置における検出信号値を記憶するためデ
ィジタルデータ処理装置に接続された検出器の入力スク
リーンに表示されるディジタル影像装置の変調伝達関数
を決定する方法に係る。
【０００２】本発明はまた、上記の方法を実行する装置
と、かかる装置を含むＸ線影像装置と、かかるＸ線装置
で使用される試験対象とに係る。
と、かかる装置を含むＸ線影像装置と、かかるＸ線装置
で使用される試験対象とに係る。
【０００３】
【従来の技術】上述の種類の方法は、１９８４年３・４
月，Ｍｅｄ．Ｐｈｙｓ．１１（２），Ｒ．Ａ．ソーンズ
、Ｇ．Ｔ．バーンズの「ディジタルＸ線装置のＭＴＦを
測定する方法」に開示されている。
月，Ｍｅｄ．Ｐｈｙｓ．１１（２），Ｒ．Ａ．ソーンズ
、Ｇ．Ｔ．バーンズの「ディジタルＸ線装置のＭＴＦを
測定する方法」に開示されている。
【０００４】この論文は、ディジタル影像装置、例えば
放射線写真又は蛍光透視装置において、変調伝達関数、
以下略してＭＴＦがどのように決定されるかを記載して
いる。影像装置のＭＴＦはその影像の質の客観的な判断
基準である。幾何学的歪みがなければ、影像装置は検出
器入力面に、そのコントラストが減じられその位相がも
との強度分布に対してシフトされた正弦波強度分布とし
て、正弦波強度分布を表示する。表示された強度分布に
おける比率を検出器の入力面の強度分布のコントラスト
で除することにより、検出器のＭＴＦが複数の空間周波
数について測定可能である。空間周波数零についてＭＴ
Ｆは１であり、周波数の増大と反対に０まで減少する。
放射線写真又は蛍光透視装置において、変調伝達関数、
以下略してＭＴＦがどのように決定されるかを記載して
いる。影像装置のＭＴＦはその影像の質の客観的な判断
基準である。幾何学的歪みがなければ、影像装置は検出
器入力面に、そのコントラストが減じられその位相がも
との強度分布に対してシフトされた正弦波強度分布とし
て、正弦波強度分布を表示する。表示された強度分布に
おける比率を検出器の入力面の強度分布のコントラスト
で除することにより、検出器のＭＴＦが複数の空間周波
数について測定可能である。空間周波数零についてＭＴ
Ｆは１であり、周波数の増大と反対に０まで減少する。
【０００５】
【発明が解決しようとする課題】影像装置のＭＴＦをよ
り効率的に測定する方法は、ＭＴＦが、影像装置のライ
ンスプレッドファンクションの１次元フーリエ変換の法
としてあらわすことができるという事実に基づいている
。ラインスプレッドファンクションは、検出器により表
示された検出器入力面のラインのイメージを表す。その
検出器入力面が個々の検出小分割面のマトリックスに小
分割されるディジタル検出器、又は検出信号値がディジ
タルメモリーに記憶されるアナログ検出器を有する影像
装置を用いると、ＭＴＦがラインスプレッドファンクシ
ョンのフーリエ変換から決定される際に標本化周波数が
低すぎるためにエイリアシングが発生するという問題が
生じる。エイリアシングは、スペクトル要素が実際に属
する周波数よりも低い周波数で表示イメージのスペクト
ル要素が発生することである。その理由は、周波数領域
では標本化された信号のスペクトルは実際のスペクトル
の周期的変形であるということである。試験対象を、第
１及び最後のスロット間の間隔の半分以下の相互スロッ
ト間隔を有する複数の等距離平行スロットとして構成し
、間隔を空間標本周波数の逆元の（ｎ＋１／２）倍に等
しくし、ｎは自然数とすると、エイリアシング問題は防
止されて、ＭＴＦは表示さるべきイメージに比例する検
出信号をフーリエ変換することにより決定しうる。例え
ば、テレビジョン撮像管に光学的に結合されるＸ線影像
増倍管が検出器として用いられる時に生ずる問題は、テ
レビジョン撮像管の入力スクリーンを走査することがラ
インスプレッドファンクションに非対称をもたらすとい
うことである。
り効率的に測定する方法は、ＭＴＦが、影像装置のライ
ンスプレッドファンクションの１次元フーリエ変換の法
としてあらわすことができるという事実に基づいている
。ラインスプレッドファンクションは、検出器により表
示された検出器入力面のラインのイメージを表す。その
検出器入力面が個々の検出小分割面のマトリックスに小
分割されるディジタル検出器、又は検出信号値がディジ
タルメモリーに記憶されるアナログ検出器を有する影像
装置を用いると、ＭＴＦがラインスプレッドファンクシ
ョンのフーリエ変換から決定される際に標本化周波数が
低すぎるためにエイリアシングが発生するという問題が
生じる。エイリアシングは、スペクトル要素が実際に属
する周波数よりも低い周波数で表示イメージのスペクト
ル要素が発生することである。その理由は、周波数領域
では標本化された信号のスペクトルは実際のスペクトル
の周期的変形であるということである。試験対象を、第
１及び最後のスロット間の間隔の半分以下の相互スロッ
ト間隔を有する複数の等距離平行スロットとして構成し
、間隔を空間標本周波数の逆元の（ｎ＋１／２）倍に等
しくし、ｎは自然数とすると、エイリアシング問題は防
止されて、ＭＴＦは表示さるべきイメージに比例する検
出信号をフーリエ変換することにより決定しうる。例え
ば、テレビジョン撮像管に光学的に結合されるＸ線影像
増倍管が検出器として用いられる時に生ずる問題は、テ
レビジョン撮像管の入力スクリーンを走査することがラ
インスプレッドファンクションに非対称をもたらすとい
うことである。
【０００６】
【課題を解決するための手段】本発明はとりわけ、その
目的として、ラインスプレッドファンクションの非対称
の程度についての情報を含むディジタル影像装置の伝達
関数が決定できる方法を提供する。この目的のため、本
発明による方法は、放射強度分布はブロック状に空間変
調され、まずデータ処理装置において変調方向を横断し
て延びる画像線にそって配置された画像線の検出信号値
の誘導値が決定され、その後に、正及び負の誘導値から
夫々、局所的周波数連続最大のシーケンスが計算され、
その包絡線から変調伝達関数が決定されることを特徴と
する。
目的として、ラインスプレッドファンクションの非対称
の程度についての情報を含むディジタル影像装置の伝達
関数が決定できる方法を提供する。この目的のため、本
発明による方法は、放射強度分布はブロック状に空間変
調され、まずデータ処理装置において変調方向を横断し
て延びる画像線にそって配置された画像線の検出信号値
の誘導値が決定され、その後に、正及び負の誘導値から
夫々、局所的周波数連続最大のシーケンスが計算され、
その包絡線から変調伝達関数が決定されることを特徴と
する。
【０００７】ブロック状の変調放射強度分布を検出器に
用いることにより、検出器のエッジスプレッドファンク
ションが決定される。エッジスプレッドファンクション
は急な強度遷移のイメージの形状を示す。エッジスプレ
ッドファンクションを微分すると上昇及び下降エッジに
つきそれぞれ正及び負の誘導値が得られる。正及び負の
誘導値をフーリエ変換することにより、ラインスプレッ
ドファンクションの非対称の測定値である２つの異なる
ＭＴＦを決定しうる。複数の平行な等距離のスロットに
よりラインスプレッドファンクションが既知の方法で変
調伝達関数に変換されると、非対称ラインスプレッドフ
ァンクションを特徴づけるために位相スペクトルについ
ての情報がさらに必要となる。この場合において変調伝
達関数のみがラインスプレッドファンクションの測定値
として採用されるときは、対称ラインスプレッドファン
クションが対称信号ｆ（ｘ）として存在すると思われ、
これにつきｆ（ｘ）＝ｆ（−ｘ）が正しいフーリエ変換
を有しその法により完全に特徴づけられる。ブロック状
の放射強度分布を用いることの別の利点は、線状放射強
度分布に比べて検出器により放射が入射するということ
である。その結果、検出信号は高信号対雑音比を有し、
ＭＴＦはより少数の測定から計算可能である。
用いることにより、検出器のエッジスプレッドファンク
ションが決定される。エッジスプレッドファンクション
は急な強度遷移のイメージの形状を示す。エッジスプレ
ッドファンクションを微分すると上昇及び下降エッジに
つきそれぞれ正及び負の誘導値が得られる。正及び負の
誘導値をフーリエ変換することにより、ラインスプレッ
ドファンクションの非対称の測定値である２つの異なる
ＭＴＦを決定しうる。複数の平行な等距離のスロットに
よりラインスプレッドファンクションが既知の方法で変
調伝達関数に変換されると、非対称ラインスプレッドフ
ァンクションを特徴づけるために位相スペクトルについ
ての情報がさらに必要となる。この場合において変調伝
達関数のみがラインスプレッドファンクションの測定値
として採用されるときは、対称ラインスプレッドファン
クションが対称信号ｆ（ｘ）として存在すると思われ、
これにつきｆ（ｘ）＝ｆ（−ｘ）が正しいフーリエ変換
を有しその法により完全に特徴づけられる。ブロック状
の放射強度分布を用いることの別の利点は、線状放射強
度分布に比べて検出器により放射が入射するということ
である。その結果、検出信号は高信号対雑音比を有し、
ＭＴＦはより少数の測定から計算可能である。
【０００８】本発明による方法の実施例は、入力スクリ
ーンの均一照射において、検出信号値は、入力スクリー
ンの画素に対応する記憶位置に記憶され、検出信号値の
誘導値がまず画像線に沿った画素につき決定され、その
後に雑音スペクトルが計算され、これは包絡線から減じ
られることを特徴とする。
ーンの均一照射において、検出信号値は、入力スクリー
ンの画素に対応する記憶位置に記憶され、検出信号値の
誘導値がまず画像線に沿った画素につき決定され、その
後に雑音スペクトルが計算され、これは包絡線から減じ
られることを特徴とする。
【０００９】放射が検出器により検出され、放射は入力
スクリーンにおいて、例えばＸ線検出器により光に変換
され、又は放射はテレビジョン撮像管により電荷パター
ンに変換されると、検出器により形成された画像は低放
射強度で雑音要素を含む。この、通例量子化歪みと呼ば
れるものは、放射線量の検出が主に平らなスペクトル（
ホワイトノイズ）での確率過程であるという事実により
生ずる。線状変調放射強度が検出器に印加されると、ホ
ワイトノイズの存在下で、検出信号値のフーリエ変換の
後にＭＴＦが得られ、これはノイズのない状態で得られ
るＭＴＦに比してホワイトノイズの標準偏差に等しい係
数だけ異なる。しかし、フーリエ変換に先立ってブロッ
ク状に変調された放射強度の場合には、検出信号値は微
分され、これは雑音スペクトルを変え、これはもはや一
定ではなく０Ｈｚで零から増加する。正確なＭＴＦは、
雑音スペクトルの別な計算及びこの雑音スペクトルをフ
ーリエ変換検出信号から減算することにより得られる。
スクリーンにおいて、例えばＸ線検出器により光に変換
され、又は放射はテレビジョン撮像管により電荷パター
ンに変換されると、検出器により形成された画像は低放
射強度で雑音要素を含む。この、通例量子化歪みと呼ば
れるものは、放射線量の検出が主に平らなスペクトル（
ホワイトノイズ）での確率過程であるという事実により
生ずる。線状変調放射強度が検出器に印加されると、ホ
ワイトノイズの存在下で、検出信号値のフーリエ変換の
後にＭＴＦが得られ、これはノイズのない状態で得られ
るＭＴＦに比してホワイトノイズの標準偏差に等しい係
数だけ異なる。しかし、フーリエ変換に先立ってブロッ
ク状に変調された放射強度の場合には、検出信号値は微
分され、これは雑音スペクトルを変え、これはもはや一
定ではなく０Ｈｚで零から増加する。正確なＭＴＦは、
雑音スペクトルの別な計算及びこの雑音スペクトルをフ
ーリエ変換検出信号から減算することにより得られる。
【００１０】本発明による方法に用いるのに適当な装置
は、検出器に結合可能なデータ処理装置で、これは検出
信号値を記憶するための記憶部分と、検出信号値の誘導
値を計算し誘導値から変調伝達関数を計算する演算部分
とを含む。
は、検出器に結合可能なデータ処理装置で、これは検出
信号値を記憶するための記憶部分と、検出信号値の誘導
値を計算し誘導値から変調伝達関数を計算する演算部分
とを含む。
【００１１】マイクロプロセッサや多目的コンピュータ
が検出信号の標本化や検出信号の例えば（５１２）２
バイトのＲＡＭでの記憶を制御する。例えばＲＯＭに記
憶されたアルゴリズムは画像線に沿って配置された検出
信号値を微分し、ファーストフーリエ変換アルゴリズム
に従って微分値からフーリエ変換値を計算する。雑音ス
ペクトルを減じた後、ＭＴＦ値はＤ／Ａ変換器を介して
ディスプレイ装置に印加される。銅製のプレートにより
形成される試験対象と協同して、かかる装置は、検出器
がＸ線影像増倍管及びそれに協同するテレビジョン撮像
装置により構成されるＸ線影像装置において使用可能で
ある。Ｘ線源により放出されたＸ線放射は銅製プレート
のスロットを通ってＸ線影像増倍管にぶつかり、そこで
ブロック状のＸ線影像は発光イメージに変換され、さら
にテレビジョン撮像装置により検出信号に変換される。 この状況下では、Ｘ線影像増倍管の入力スクリーン状に
表示された際の銅製プレートのスロットの幅は、画素に
比して大きい。実質的に０．５ｍｍの厚さと、５ｍｍの
スロット幅と、５ｍｍの相対スロット間隔を有する銅製
の試験対称を用いると、良い試験結果が得られ、試験対
称はＸ線影像増倍管の入力スクリーンに対して設置され
る。
が検出信号の標本化や検出信号の例えば（５１２）２
バイトのＲＡＭでの記憶を制御する。例えばＲＯＭに記
憶されたアルゴリズムは画像線に沿って配置された検出
信号値を微分し、ファーストフーリエ変換アルゴリズム
に従って微分値からフーリエ変換値を計算する。雑音ス
ペクトルを減じた後、ＭＴＦ値はＤ／Ａ変換器を介して
ディスプレイ装置に印加される。銅製のプレートにより
形成される試験対象と協同して、かかる装置は、検出器
がＸ線影像増倍管及びそれに協同するテレビジョン撮像
装置により構成されるＸ線影像装置において使用可能で
ある。Ｘ線源により放出されたＸ線放射は銅製プレート
のスロットを通ってＸ線影像増倍管にぶつかり、そこで
ブロック状のＸ線影像は発光イメージに変換され、さら
にテレビジョン撮像装置により検出信号に変換される。 この状況下では、Ｘ線影像増倍管の入力スクリーン状に
表示された際の銅製プレートのスロットの幅は、画素に
比して大きい。実質的に０．５ｍｍの厚さと、５ｍｍの
スロット幅と、５ｍｍの相対スロット間隔を有する銅製
の試験対称を用いると、良い試験結果が得られ、試験対
称はＸ線影像増倍管の入力スクリーンに対して設置され
る。
【００１２】
【実施例】以下図面を参照しながら本発明による方法及
び装置の実施例を幾つか詳細に示す。
び装置の実施例を幾つか詳細に示す。
【００１３】図１は、試験対象５を通るＸ線ビーム３を
放出するＸ線源１を示す。試験対象５はＸ線ビームの強
度を空間的にブロック状に変調する。放射強度分布は、
Ｘ線影像増倍管９の入力スクリーン７で発光イメージに
変換され、このイメージはテレビジョン撮像装置１によ
り検出され、検出信号に変換される。検出信号はデータ
処理装置１３に印加され、入力スクリーン７の試験対象
５の影像の画素に対応する記憶要素に記憶される。
放出するＸ線源１を示す。試験対象５はＸ線ビームの強
度を空間的にブロック状に変調する。放射強度分布は、
Ｘ線影像増倍管９の入力スクリーン７で発光イメージに
変換され、このイメージはテレビジョン撮像装置１によ
り検出され、検出信号に変換される。検出信号はデータ
処理装置１３に印加され、入力スクリーン７の試験対象
５の影像の画素に対応する記憶要素に記憶される。
【００１４】図２は約０．５ｍｍの厚さの銅製ディスク
の形の試験対象５を示す。ディスクは、５ｍｍのスロッ
ト幅ｌと、５ｍｍの相対スロット間隔ｄを有する４つの
スロットを含み、これらは次の４つのスロット群に直交
する。Ｘ線は、スロットにより空間的及びブロック状に
変調され、ＭＦＴは２つの相互に直交する方向において
決定される。
の形の試験対象５を示す。ディスクは、５ｍｍのスロッ
ト幅ｌと、５ｍｍの相対スロット間隔ｄを有する４つの
スロットを含み、これらは次の４つのスロット群に直交
する。Ｘ線は、スロットにより空間的及びブロック状に
変調され、ＭＦＴは２つの相互に直交する方向において
決定される。
【００１５】図３は、対象５のイメージにおける画像線
Ｍに沿って配置された画素に対応する記憶位置において
データ処理装置に記憶された検出信号値を示す。画素の
数は水平上にプロットされ、検出信号値に割り当てられ
０〜２５５の値を有するグレーレベルは垂直線上にプロ
ットされる。図は検出器の急な強度遷移への応答を示し
、エッジスプレッドファンクションを表す。
Ｍに沿って配置された画素に対応する記憶位置において
データ処理装置に記憶された検出信号値を示す。画素の
数は水平上にプロットされ、検出信号値に割り当てられ
０〜２５５の値を有するグレーレベルは垂直線上にプロ
ットされる。図は検出器の急な強度遷移への応答を示し
、エッジスプレッドファンクションを表す。
【００１６】図４は検出信号値から決定される誘導値を
示す。ラインスプレッドファンクションＬＳＦ（ｘｋ
）は、ここでｘｋ を画像線における画素の数とすると
、以下のようにエッジスプレッドファンクションＥＳＦ
（ｘｋ ）の誘導値により与えられる。
示す。ラインスプレッドファンクションＬＳＦ（ｘｋ
）は、ここでｘｋ を画像線における画素の数とすると
、以下のようにエッジスプレッドファンクションＥＳＦ
（ｘｋ ）の誘導値により与えられる。
【００１７】
【数１】
上述のような図３における検出信号値の差は、４つの等
距離の正及び４つの等距離の負のピークに帰結する。正
のピークと負のそれとを別々に考えることにより、ライ
ンスプレッドファンクションの非対象を表す２つのＭＴ
Ｆが計算可能である。
距離の正及び４つの等距離の負のピークに帰結する。正
のピークと負のそれとを別々に考えることにより、ライ
ンスプレッドファンクションの非対象を表す２つのＭＴ
Ｆが計算可能である。
【００１８】図５の（ａ）は、ＭＴＦに等しいラインス
プレッドファンクションのフーリエ変換の概略を表す。 ラインスプレッドファンクションの有限シーケンス（こ
こでは４つ）はパルス列によるラインスプレッドファン
クションのコンボリューションにより表され、この後に
有限の幅の方形関数により乗算される。周波数領域では
、有限のシーケンスのフーリエ変換はＭＦＴにパルス列
を乗ずることにより表され、この後にｓｉｎｃ関数［ｓ
ｉｎｃ（ｘ）＝ｓｉｎ（ｘ）／ｘ］でのコンボリューシ
ョンが続く。図４から、ラインスプレッドファンクショ
ンが離散点でデータ処理装置１３に記憶されることが明
らかである。離散点において抽出された信号の正確な再
構成のため、ナイキスト基準によると、２つの画素間の
間隔Ｔｎは、ｆｍ をＭＴＦにおいて生ずる最高周波数
とすると、以下を満たすことを要する：Ｔｎ≦［１／（
２ｆｍ ）］。しかし、試験対象のイメージにおける２
つの近接画素間の間隔がΔｘに相当し、Δｘ＞Ｔｎとす
ると、ラインスプレッドファンクションから正確なＭＴ
Ｆを決定することは一定の条件下ではまだ可能である。 これは図５の（ａ）−（ｄ）に示されている。図４に示
されるようなラインスプレッドファンクションはフーリ
エ変換によりＭＴＦに変換され、変換の後、図５の（ｄ
）に示されるような形状をとる。ラインスプレッドファ
ンクションはラインスプレッドファンクションでコンボ
リュートされたパルス関数のシーケンスとして示される
。周波数領域では、これは図５の（ａ）及び図５の（ｂ
）の乗算として示される。ラインスプレッドファンクシ
ョンを相互間隔Δｘとして離散点で抽出することは、周
波数領域においては、相対間隔ｆｓ でのパルス関数の
シーケンスでのコンボリューションにより表され、ここ
でｆｓ ＞２ｆｍ を満たす。これら処理作用の結果は
、図５の（ｄ）に示される。元のＭＦＴは離散ＭＴＦ値
の包絡線から再構成可能である。明確を期すために、図
５の（ｄ）は、図５の（ｂ）におけるパルス列が無限に
長いと仮定したパルス関数を示す。試験対象におけるス
ロットの数が有限なのでこれはあり得ず、図５の（ｄ）
におけるパルス関数は所与の幅を有し（それらは実際ｓ
ｉｎｃ関数である）、ＭＴＦは局所的最大の包絡線から
得られる。
プレッドファンクションのフーリエ変換の概略を表す。 ラインスプレッドファンクションの有限シーケンス（こ
こでは４つ）はパルス列によるラインスプレッドファン
クションのコンボリューションにより表され、この後に
有限の幅の方形関数により乗算される。周波数領域では
、有限のシーケンスのフーリエ変換はＭＦＴにパルス列
を乗ずることにより表され、この後にｓｉｎｃ関数［ｓ
ｉｎｃ（ｘ）＝ｓｉｎ（ｘ）／ｘ］でのコンボリューシ
ョンが続く。図４から、ラインスプレッドファンクショ
ンが離散点でデータ処理装置１３に記憶されることが明
らかである。離散点において抽出された信号の正確な再
構成のため、ナイキスト基準によると、２つの画素間の
間隔Ｔｎは、ｆｍ をＭＴＦにおいて生ずる最高周波数
とすると、以下を満たすことを要する：Ｔｎ≦［１／（
２ｆｍ ）］。しかし、試験対象のイメージにおける２
つの近接画素間の間隔がΔｘに相当し、Δｘ＞Ｔｎとす
ると、ラインスプレッドファンクションから正確なＭＴ
Ｆを決定することは一定の条件下ではまだ可能である。 これは図５の（ａ）−（ｄ）に示されている。図４に示
されるようなラインスプレッドファンクションはフーリ
エ変換によりＭＴＦに変換され、変換の後、図５の（ｄ
）に示されるような形状をとる。ラインスプレッドファ
ンクションはラインスプレッドファンクションでコンボ
リュートされたパルス関数のシーケンスとして示される
。周波数領域では、これは図５の（ａ）及び図５の（ｂ
）の乗算として示される。ラインスプレッドファンクシ
ョンを相互間隔Δｘとして離散点で抽出することは、周
波数領域においては、相対間隔ｆｓ でのパルス関数の
シーケンスでのコンボリューションにより表され、ここ
でｆｓ ＞２ｆｍ を満たす。これら処理作用の結果は
、図５の（ｄ）に示される。元のＭＦＴは離散ＭＴＦ値
の包絡線から再構成可能である。明確を期すために、図
５の（ｄ）は、図５の（ｂ）におけるパルス列が無限に
長いと仮定したパルス関数を示す。試験対象におけるス
ロットの数が有限なのでこれはあり得ず、図５の（ｄ）
におけるパルス関数は所与の幅を有し（それらは実際ｓ
ｉｎｃ関数である）、ＭＴＦは局所的最大の包絡線から
得られる。
【００１９】図６において、雑音スペクトルの減算の後
、局所的最大の包絡線から、ＭＴＦが決定される。図４
の誘導値はラインスプレッドファンクションＬＳＦ（ｘ
）のパルス列ｃｏｍｂ（ｘ／ｄ）でのコンボリューショ
ンとして示される。試験対象の有限の幅は、方形関数ｒ
ｅｃｔ（ｘ／ｃ）による乗算により与えられ、相互間隔
Δｘでの離散点における抽出はパルス列ｃｏｍｂ（ｋΔ
ｘ）による乗算により与えられる。データ処理装置１３
では、画像線での画素に対応する記憶位置における微分
の後、検出信号値Ｉ（ｘ）は記憶され、それにつき以下
が当てはまる： Ｉ（ｘ）＝［ＬＳＦ（ｘ）＊Ｃｏｍｂ（ｘ／ｄ）］・ｒ
ｅｃｔ（ｘ／ｃ） ・ｃｏ
ｍｂ（ｋΔｘ）微分された検出信号値のフーリエ変換の
後、以下で与えられる信号が得られる： ＤＦＴ［Ｉ（ｘ）］＝［ＭＴＦ（ｆ）・ｃｏｍｂ（ｆｄ
）］＊ｓｉｎｃ（ｆｃ）
＊ｃｏｍｂ（ｆｓ ）ここで
、演算ＤＦＴは、離散フーリエ変換を意味する。 上記の公式から、ＭＴＦは複数の連続ｓｉｎｃ関数（ｃ
ｏｍｂ（ｆｄ）＊ｓｉｎｃ（ｆｃ）；これは標本化周波
数ｆｓ で繰り返す）の包絡線であることがわかる。
、局所的最大の包絡線から、ＭＴＦが決定される。図４
の誘導値はラインスプレッドファンクションＬＳＦ（ｘ
）のパルス列ｃｏｍｂ（ｘ／ｄ）でのコンボリューショ
ンとして示される。試験対象の有限の幅は、方形関数ｒ
ｅｃｔ（ｘ／ｃ）による乗算により与えられ、相互間隔
Δｘでの離散点における抽出はパルス列ｃｏｍｂ（ｋΔ
ｘ）による乗算により与えられる。データ処理装置１３
では、画像線での画素に対応する記憶位置における微分
の後、検出信号値Ｉ（ｘ）は記憶され、それにつき以下
が当てはまる： Ｉ（ｘ）＝［ＬＳＦ（ｘ）＊Ｃｏｍｂ（ｘ／ｄ）］・ｒ
ｅｃｔ（ｘ／ｃ） ・ｃｏ
ｍｂ（ｋΔｘ）微分された検出信号値のフーリエ変換の
後、以下で与えられる信号が得られる： ＤＦＴ［Ｉ（ｘ）］＝［ＭＴＦ（ｆ）・ｃｏｍｂ（ｆｄ
）］＊ｓｉｎｃ（ｆｃ）
＊ｃｏｍｂ（ｆｓ ）ここで
、演算ＤＦＴは、離散フーリエ変換を意味する。 上記の公式から、ＭＴＦは複数の連続ｓｉｎｃ関数（ｃ
ｏｍｂ（ｆｄ）＊ｓｉｎｃ（ｆｃ）；これは標本化周波
数ｆｓ で繰り返す）の包絡線であることがわかる。
【００２０】図７は、正及び負の誘導値から決定された
２つのＭＦＴを示す。ラインスプレッドファンクション
の非対称性ゆえに、ＭＦＴ間に相当の違いがあることが
わかる。
２つのＭＦＴを示す。ラインスプレッドファンクション
の非対称性ゆえに、ＭＦＴ間に相当の違いがあることが
わかる。
【００２１】図８は、データ処理装置１３及びディスプ
レイ装置３４からなる変調伝達関数の決定用の装置を示
す。ゲート回路２０は検出器への接続用の端子１９を有
し、アナログ−ディジタル変換器２２に接続される。ア
ナログ−ディジタル変換器はデータバス２４に接続され
、このバスＲＯＭ２６と、ＲＡＭ２８と、ＣＰＵ３０と
に接続される。データバスはまた、テレビジョンモニタ
のようなディスプレイ装置３４に接続されたディジタル
−アナログ変換器３２に接続される。制御線３６を介し
て、影像装置のＭＦＴの決定の間、ＣＰＵはゲート回路
２０を活性化し、端子１９に接続された検出器の検出信
号値はＡ／Ｄ変換器２２に到着する。データバスを介し
て、例えば２８ の段階でディジタル化された検出信号
値はメモリー２８のデータ入力に印加され、このメモリ
ーは、例えば、夫々８ビットを有する５１２２ のマト
リックスの記憶位置を有する。メモリー２８における各
記憶位置は、検出器により検出されたイメージにおける
画素に対応する。メモリー２８における検出信号値の記
憶は、アドレスバス３８を介してメモリー２８の記憶位
置にアドレスを送るＣＰＵ３０により制御される。メモ
リー２６はＣＰＵ３０に対する指示を含んでおり、これ
はまず、メモリー２８に記憶された信号値の微分からな
り、この後に、微分値のフーリエ変換と、フーリエ変換
された微分検出信号値の局所的最大の包絡線からのＭＴ
Ｆの決定が続く。行又は列に沿って配置された記憶位置
の内容はその後にデータバス２４を介してＤ／Ａ変換器
３２に印加され、例えばテレビジョンモニタのようなデ
ィスプレイ装置３４に表示されたり、プリンタで印刷さ
れたりする。
レイ装置３４からなる変調伝達関数の決定用の装置を示
す。ゲート回路２０は検出器への接続用の端子１９を有
し、アナログ−ディジタル変換器２２に接続される。ア
ナログ−ディジタル変換器はデータバス２４に接続され
、このバスＲＯＭ２６と、ＲＡＭ２８と、ＣＰＵ３０と
に接続される。データバスはまた、テレビジョンモニタ
のようなディスプレイ装置３４に接続されたディジタル
−アナログ変換器３２に接続される。制御線３６を介し
て、影像装置のＭＦＴの決定の間、ＣＰＵはゲート回路
２０を活性化し、端子１９に接続された検出器の検出信
号値はＡ／Ｄ変換器２２に到着する。データバスを介し
て、例えば２８ の段階でディジタル化された検出信号
値はメモリー２８のデータ入力に印加され、このメモリ
ーは、例えば、夫々８ビットを有する５１２２ のマト
リックスの記憶位置を有する。メモリー２８における各
記憶位置は、検出器により検出されたイメージにおける
画素に対応する。メモリー２８における検出信号値の記
憶は、アドレスバス３８を介してメモリー２８の記憶位
置にアドレスを送るＣＰＵ３０により制御される。メモ
リー２６はＣＰＵ３０に対する指示を含んでおり、これ
はまず、メモリー２８に記憶された信号値の微分からな
り、この後に、微分値のフーリエ変換と、フーリエ変換
された微分検出信号値の局所的最大の包絡線からのＭＴ
Ｆの決定が続く。行又は列に沿って配置された記憶位置
の内容はその後にデータバス２４を介してＤ／Ａ変換器
３２に印加され、例えばテレビジョンモニタのようなデ
ィスプレイ装置３４に表示されたり、プリンタで印刷さ
れたりする。
【００２２】
【発明の効果】ブロック状の変調放射強度分布を検出器
に用いることにより、検出器のエッジスプレッドファン
クションが決定される。エッジスプレッドファンクショ
ンは急な強度遷移のイメージの形状を示す。エッジスプ
レッドファンクションを微分すると上昇及び下降エッジ
につきそれぞれ正及び負の誘導値が得られる。正及び負
の誘導値をフーリエ変換することにより、ラインスプレ
ッドファンクションの非対称の測定値である２つの異な
るＭＴＦを決定しうる。
に用いることにより、検出器のエッジスプレッドファン
クションが決定される。エッジスプレッドファンクショ
ンは急な強度遷移のイメージの形状を示す。エッジスプ
レッドファンクションを微分すると上昇及び下降エッジ
につきそれぞれ正及び負の誘導値が得られる。正及び負
の誘導値をフーリエ変換することにより、ラインスプレ
ッドファンクションの非対称の測定値である２つの異な
るＭＴＦを決定しうる。
【図１】影像装置のＭＴＦを決定するための装置が設け
られたＸ線影像装置を示す。
られたＸ線影像装置を示す。
【図２】本発明による試験対象を示す。
【図３】図２の試験対象の影像における検出信号値を示
す。
す。
【図４】図３で示された検出信号の誘導値を示す。
【図５】（ａ）から（ｄ）は離散画素において抽出され
たラインスプレッドファンクションのフーリエ変換を示
す。
たラインスプレッドファンクションのフーリエ変換を示
す。
【図６】図５の（ｄ）によるスペクトルでの局所的最大
の包絡線からのＭＴＦの決定を示す。
の包絡線からのＭＴＦの決定を示す。
【図７】誘導値から決定された正及び負のエッジ遷移の
ＭＴＦを示す。
ＭＴＦを示す。
【図８】影像装置のＭＴＦを決定するための装置を示す
。
。
１
Ｘ線源３
Ｘ線ビーム５
試験対象７
入力スクリーン１３
デー
タ処理装置１９
端子２０
ゲート回路２２
アナログ−ディジタル変換器 ２４
データバス２６
ＲＯＭ２８
ＲＡ
Ｍ３０
ＣＰＵ３２
ディジタル−アナログ変
換器
Ｘ線源３
Ｘ線ビーム５
試験対象７
入力スクリーン１３
デー
タ処理装置１９
端子２０
ゲート回路２２
アナログ−ディジタル変換器 ２４
データバス２６
ＲＯＭ２８
ＲＡ
Ｍ３０
ＣＰＵ３２
ディジタル−アナログ変
換器
Claims (6)
- 【請求項１】試験対象により空間変調された放射強度分
布が、入力スクリーンの画素に対応するメモリー位置に
おける検出信号値を記憶するためディジタルデータ処理
装置に接続された検出器の入力スクリーンに表示される
ディジタル影像装置の変調伝達関数を決定する方法であ
って、放射強度分布はブロック状に空間変調され、まず
データ処理装置において変調方向を横断して延びる画像
線にそって配置された画像線の検出信号値の誘導値が決
定され、その後に、正及び負の誘導値から夫々、局所的
周波数連続最大のシーケンスが計算され、その包絡線か
ら変調伝達関数が決定されることを特徴とする変調伝達
関数決定方法。 - 【請求項２】入力スクリーンの均一照射において、検出
信号値は、入力スクリーンの画素に対応する記憶位置に
記憶され、検出信号値の誘導値がまず画像線に沿った画
素につき決定され、その後に雑音スペクトルが計算され
、これは包絡線から減じられることを特徴とする請求項
１記載の方法。 - 【請求項３】検出信号値を記憶するための記憶部分と、
検出信号値の誘導値を計算し誘導値から変調伝達関数を
計算する演算部分とを含む、検出器に結合可能なデータ
処理装置からなる請求項１又は２記載の方法が応用可能
な装置。 - 【請求項４】請求項３記載の装置に接続されたＸ線源と
Ｘ線検出器を含むＸ線影像装置。 - 【請求項５】スロットのついた放射吸収板を有すること
を特徴とする請求項４記載のＸ線影像装置に用いられる
のに適当な試験対象。 - 【請求項６】放射吸収板は銅製であり実質的に０．５ｍ
ｍの厚さを有することを特徴とする請求項５記載の試験
対象。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| NL9001265 | 1990-06-05 | ||
| NL9001265A NL9001265A (nl) | 1990-06-05 | 1990-06-05 | Werkwijze voor bepaling van een modulatie-overdrachtfunctie van een digitaal beeldvormend systeem. |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH04233377A true JPH04233377A (ja) | 1992-08-21 |
Family
ID=19857188
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP3131414A Pending JPH04233377A (ja) | 1990-06-05 | 1991-06-03 | ディジタル影像装置の変調伝達関数決定方法 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US5191621A (ja) |
| EP (1) | EP0460749A1 (ja) |
| JP (1) | JPH04233377A (ja) |
| NL (1) | NL9001265A (ja) |
Families Citing this family (35)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2633793A1 (fr) * | 1988-07-01 | 1990-01-05 | Gen Electric Cgr | Procede de correction de la distorsion d'images radiologiques |
| US5661773A (en) * | 1992-03-19 | 1997-08-26 | Wisconsin Alumni Research Foundation | Interface for radiation therapy machine |
| US6529639B1 (en) * | 1992-04-21 | 2003-03-04 | Raytheon Company | Modulation transfer function characterization employing phased slit reticle |
| DE4238268C2 (de) * | 1992-11-12 | 1998-07-02 | Siemens Ag | Verfahren und Vorrichtung zur Abnahme- und Konstanzprüfung filmloser Dental-Röntgengeräte |
| JPH06304159A (ja) * | 1993-04-27 | 1994-11-01 | Fuji Photo Film Co Ltd | 画像の撮影体位認識方法 |
| JPH08511452A (ja) * | 1993-06-09 | 1996-12-03 | ウイスコンシン アラムナイ リサーチ フオンデーシヨン | 抑制された回転自由度を有する放射治療システムおよび動的コリメータ |
| US5420441A (en) * | 1993-11-23 | 1995-05-30 | Eastman Kodak Company | Automated technique for calibrating a storage phosphor reader |
| US5515447A (en) * | 1994-06-07 | 1996-05-07 | United Parcel Service Of America, Inc. | Method and apparatus for locating an acquisition target in two-dimensional images by detecting symmetry in two different directions |
| US5581631A (en) * | 1994-09-20 | 1996-12-03 | Neopath, Inc. | Cytological system image collection integrity checking apparatus |
| EP0710014A1 (en) * | 1994-10-27 | 1996-05-01 | Agfa-Gevaert N.V. | Method for verification of the frequency response of a digital radiographic read out system |
| JP2001503922A (ja) * | 1995-03-03 | 2001-03-21 | ノースロップ グラマン コーポレーション | Ｍｔｆ試験用の適応ウィンドによる固定パターン補償 |
| US5867211A (en) * | 1995-10-27 | 1999-02-02 | Northrop Grumman Corporation | Apparatus and method for removing offset and shading from a video signal |
| US5859930A (en) * | 1995-12-06 | 1999-01-12 | Fpr Corporation | Fast pattern recognizer utilizing dispersive delay line |
| US5760403A (en) * | 1996-04-18 | 1998-06-02 | Loral Fairchild Corp. | High modulation transfer function CCD X-ray image sensor apparatus and method |
| US6233349B1 (en) * | 1997-06-20 | 2001-05-15 | General Electric Company | Apparata and methods of analyzing the focal spots of X-ray tubes |
| US6231231B1 (en) | 1999-06-24 | 2001-05-15 | General Electric Company | Modular interchangeable phantoms for multiple x-ray systems |
| DE10153237A1 (de) * | 2001-10-31 | 2003-05-15 | Lfk Gmbh | Verfahren und Vorrichtung zur automatisierten Bestimmung der Modulations-Transfer-Funktion (MTF) von Focal-Plane-Array (FPA)- Kameras |
| US7334613B2 (en) | 2002-10-15 | 2008-02-26 | Black & Decker Inc. | Router base securing mechanism |
| US7073993B2 (en) * | 2002-10-15 | 2006-07-11 | Porter-Cable Corporation | Switch assembly |
| US7334614B2 (en) * | 2002-10-15 | 2008-02-26 | Black & Decker Inc. | Depth adjustment mechanism |
| US7451791B2 (en) | 2002-10-15 | 2008-11-18 | Black & Decker Inc. | Handle assembly |
| US20060191597A1 (en) * | 2002-10-15 | 2006-08-31 | Black & Decker Inc. | Handle assembly |
| US7316528B2 (en) | 2002-10-15 | 2008-01-08 | Black & Decker Inc. | Ergonomic router assembly |
| US6986369B1 (en) | 2002-11-12 | 2006-01-17 | Porter-Cable Corporation | Router height adjustment apparatus |
| US7089979B2 (en) * | 2003-05-01 | 2006-08-15 | Black & Decker Inc. | Ergonomic router |
| US20060102249A1 (en) * | 2003-05-01 | 2006-05-18 | Cooper Randy G | Router with drive shaft lock mechanism |
| US7275900B1 (en) | 2003-07-25 | 2007-10-02 | Black & Decker Inc. | Router elevating mechanism |
| CN1306454C (zh) * | 2003-11-06 | 2007-03-21 | Ge医疗系统环球技术有限公司 | 调制传递函数测量方法和系统 |
| USD546654S1 (en) | 2004-01-29 | 2007-07-17 | Black & Decker Inc. | Router with plunge base |
| USD531871S1 (en) | 2005-06-16 | 2006-11-14 | Black & Decker Inc. | Router |
| US9477330B2 (en) * | 2013-11-05 | 2016-10-25 | Microsoft Technology Licensing, Llc | Stylus tilt tracking with a digitizer |
| US10067029B2 (en) | 2016-02-12 | 2018-09-04 | Google Llc | Systems and methods for estimating modulation transfer function in an optical system |
| CN106841115A (zh) * | 2017-02-07 | 2017-06-13 | 南京理工大学 | 基于空间频率特征的海面溢油检测方法 |
| CN109862351B (zh) * | 2019-03-05 | 2020-07-17 | 横店集团东磁有限公司 | 一种摄像头模组解析度检验方法 |
| US12165341B2 (en) * | 2020-11-24 | 2024-12-10 | Applied Materials, Inc. | Optical resolution measurement method for optical devices |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4193089A (en) * | 1976-03-05 | 1980-03-11 | Picker Corporation | Television radiation imaging system and method |
| US4259579A (en) * | 1979-11-01 | 1981-03-31 | The United States Of America As Represented By The Secretary Of The Air Force | Waveguide line spread function analyzing apparatus |
| IL70214A (en) * | 1983-11-13 | 1987-10-20 | Elscint Ltd | Image contrast enhancement arrangement |
| US4947323A (en) * | 1986-05-22 | 1990-08-07 | University Of Tennessee Research Corporation | Method and apparatus for measuring small spatial dimensions of an object |
| DD259692A1 (de) * | 1986-12-31 | 1988-08-31 | Mfnv | Verfahren zur lasergestuetzten ermittlung der modulationsuebertragungsfunktion aus roentgenstrahlungsverteilungen |
| US5056130A (en) * | 1989-10-06 | 1991-10-08 | The United States Of America As Represented By The Administrator, National Aeronautics And Space Administration | Computerized tomography calibrator |
-
1990
- 1990-06-05 NL NL9001265A patent/NL9001265A/nl not_active Application Discontinuation
-
1991
- 1991-06-03 US US07/709,665 patent/US5191621A/en not_active Expired - Fee Related
- 1991-06-03 EP EP91201335A patent/EP0460749A1/en not_active Withdrawn
- 1991-06-03 JP JP3131414A patent/JPH04233377A/ja active Pending
Also Published As
| Publication number | Publication date |
|---|---|
| NL9001265A (nl) | 1992-01-02 |
| US5191621A (en) | 1993-03-02 |
| EP0460749A1 (en) | 1991-12-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH04233377A (ja) | ディジタル影像装置の変調伝達関数決定方法 | |
| Cunningham et al. | A method for modulation transfer function determination from edge profiles with correction for finite‐element differentiation | |
| Fujita et al. | Investigation of basic imaging properties in digital radiography. 6. MTFs of II–TV digital imaging systems | |
| JP3843169B2 (ja) | 記録部材における放射線画像の欠陥を補正する方法 | |
| Williams et al. | Noise power spectra of images from digital mammography detectors | |
| US4574311A (en) | Random array sensing devices | |
| EP1265194A2 (en) | Method, system and program for processing and storing radiation images | |
| US4653909A (en) | Sinusoidal determination of limiting optical resolution | |
| JPH067694B2 (ja) | ビデオカメラの伝達関数決定法 | |
| US11915402B2 (en) | Method and system to calculate the point spread function of a digital image detector system based on a MTF modulated quantum-noise measurement | |
| US4972451A (en) | Method of determining a modulation transfer function of a digital imaging system | |
| Ichikawa et al. | MTF measurement method for medical displays by using a bar‐pattern image | |
| JPH06124423A (ja) | トラック幅誤差検査装置 | |
| US6058220A (en) | Device for scanning an X-ray image | |
| Davies | Digital mammography—the comparative evaluation of film digitizers | |
| EP0818020A4 (en) | STABLE COMPENSATION FOR FTM MEASUREMENT | |
| Bunch et al. | Interpreting the noise power spectra of digitally printed images | |
| Goldman | Fluoroscopic performance tests using a portable computer/frame grabber: Wiener spectra measurements | |
| Ishimitsu et al. | Radiological laser film scanner sampling artifact | |
| Ichikawa et al. | Novel MTF measurement method for medical image viewers using a bar pattern image | |
| Ishimitsul et al. | Radiological laser film scanner sampling artifact | |
| JPH05240807A (ja) | 電子回折画像のピーク抽出方法 | |
| Wallingford et al. | Real-time X-ray image processing; Techniques for sensitivity improvement using low-cost equipment | |
| EP1168827A2 (en) | Apparatus and method for removing offset and shading from a video signal | |
| Ukishima et al. | An evaluation of sharpness in different image displays used for medical imaging |