JPH0312785A

JPH0312785A - 画像処理方法

Info

Publication number: JPH0312785A
Application number: JP1147537A
Authority: JP
Inventors: Kazuo Shimura; 一男志村; Wataru Ito; 渡伊藤
Original assignee: Fuji Photo Film Co Ltd
Current assignee: Fujifilm Holdings Corp
Priority date: 1989-06-09
Filing date: 1989-06-09
Publication date: 1991-01-21

Abstract

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

Description

【発明の詳細な説明】（産業上の利用分野）本発明は、グリッドを使用して撮影を行なうことにより
、上記グリッドに対応する縞模様が被写体像に重畳され
た放射線画像の画像処理方法に関するものである。

（従来の技術）記録された放射線画像を読み取って画像データを得、こ
の画像データに適切な画像処理を施した後、画像を再生
記録することは種々の分野で行なわれている。たとえば
ζ後の画像処理に適合するように設計されたガンマ値の
低いＸ線フィルムを用いてＸ線画像を記録し、このＸ線
画像が記録されたフィルムからＸ線画像を読み取って電
気信号に変換し、この電気信号（画像データ）に画像処
理を施した後コピー写真等に可視像として再生すること
により、コントラスト、シャープネス、粒状性等の画質
性能の良好な再生画像を得ることのできるシステムが開
発されている（特公昭８１−５１９３号公報参照）。

また本願出願人により、放射線（Ｘ線、α線。

β線、γ線、電子線、紫外線等）を照射するとこの放射
線エネルギーの一部が蓄積され、その後可視光等の励起
光を照射すると蓄積されたエネルギーに応じて輝尽発光
を示す蓄積性蛍光体（輝尽性蛍光体）を利用して、人体
等の被写体の放射線画像を一部シート状の蓄積性蛍光体
に撮影記録し、この蓄積性蛍光体シートをレーザー光等
の励起光で走査して輝尽発光光を生せしめ、得られた輝
尽発光光を光電的に読み取って画像データを得、この画
像データに基づき被写体の放射線画像を写真感光材料等
の記録材料、ＣＲＴ等に可視像として出力させる放射線
画像記録再生システムがすでに提案されている（特開昭
５５−１２４２９号、同５６−１１３９５号、同５５−
１６３４７２号、同５Ｂ−１０４８４５号、同５５−１
１６３４０号等）。

このシステムは、従来の銀塩写真を用いる放射線写真シ
ステムと比較して極めて広い放射線露出域にわたって画
像を記録しうるという実用的な利点を有している。すな
わち、蓄積性蛍光体においては、放射線露光量に対して
蓄積後に励起によって輝尽発光する発光光の光量が極め
て広い範囲にわたって比例することが認められており、
従って種々の撮影条件により放射線露光量がかなり大幅
に変動しても、蓄積性蛍光体シートより放射される輝尽
発光光の光量を読取ゲインを適当な値に設定して光電変
換手段により読み取って電気信号に変換し、この電気信
号を用いて写真感光材料等の記録材料、ＣＲＴ等の表示
装置に放射線画像を可視像として出力させることによっ
て、放射線露光量の変動に影響されない放射線画像を得
ることができる。

（発明が解決しようとする課題）上記Ｘ線フィルムや蓄積性蛍光体シート等の記録シート
に被写体の放射線画像を撮影記録する際に、被写体で散
乱された放射線が記録シートに照射されないように、４
．０本／ｍ程度の細かなピッチで放射線の透過しないた
とえば鉛と透過しやすいたとえばアルミニウムや木材等
が交互に配置されたグリッドを被写体と記録シートとの
間に配置して撮影を行なうことがある。グリッドを用い
て撮影を行なうと被写体により散乱された放射線が記録
シートに照射されにくいため、被写体の放射線画像のコ
ントラストを向上させることができるが、一方被写体像
とともにグリッドに対応した細かな縞模様が記録される
。

一方、放射線画像が記録された記録シートから画像デー
タを得る放射線画像読取装置は、通常、放射線画像情報
を担持した光を光電的に読み取り、画像情報として必要
な最高の空間周波数（以下、この必要な最高の空間周波
数をｆｓｓで表わす。）に対応したサンプリング間隔Δ
ｘ−１／（２・ｆ　ｓｓ）でサンプリングしてディジタ
ル信号としての画像データを得るように構成される。こ
のようにして得られた画像データには、被写体の放射線
画像を表わす有用な情報だけでなく、たとえ上記縞模様
の空間周波数が必要な最高の空間周波数ｆｓｓより高周
波であっても上記縞模様に起因するノイズも含まれるこ
とがある。

第５Ａ図の実線で示すグラフ３３は、記録シートに記録
された放射線画像（縞模様を含む）の、縞模様と直交す
る方向の空間周波数特性を表わすものである。

ここでは４．０本／＃のグリッドを用いて撮影を行ない
、縞模様の空間周波数が４ｃｙｃｌｅ　／ｍｍにあるも
のとし、被写体の放射線画像の再現に必要な最高の空間
周波数ｆｓｓは２．５ｅｙｃｌｅ／ｍであるとして説明
する。

第５Ａ図の破線で示すグラフ３７は、被写体の放射線画
像の再現に必要なｆ　ｓｓ　−２，５（ｃｙｃｌｅ　／
ｍ）以下の空間周波数帯の情報を得るために、ｆｓｓ−
２，５（ｃｙｃｌｅ　／ｍ）に対応したサンプリング間
隔Δｘ　−１／　（２拳ｆ　ｓｓ）　−０，２（＃）　
、すなわち１閣につき５回サンプリングを行なった場合
のノイズの混入の状態を示している。この場合実線３３
で示すグラフがｆ　５ｓ−２，５（ｃｙｃｌｅ　／Ｉｎ
ＩＲ）で折り返した位置にいわゆるエリアジングと呼ば
れるノイズが混入する。したがってこの場合縞模様の空
間周波数４ｃｙｃｌｅ／ａｍのエリアジングは１ｃｙｅ
ｌｅ／閣に生じる。

第５Ａ図に一点鎖線で示すグラフ３８は、上記サンプリ
ング間隔Δｘ＝１／　（２・ｆｓｓ）　−０，２（＃）
でサンプリングして得られた画像データが担持する放射
線画像の空間周波数特性を示している。

１ｃｙｃｌｅ／ｉｓの位置に上記縞模様に対応したノイ
ズが混入し、可視画像を再生記録した際、該可視画像が
記録シートに記録された放射線画像と等倍率である場合
、該可視画像にｌ　ｃｙｃｌｅ　／ｒｔｎの縞模様とし
て現われる結果となる。放射線画像を観察するに際し、
たとえグリッドに起因する縞模様がそれほど目立たない
空間周波数帯にある場合であっても、サンプリングして
画像データを得ると目立ちやすい空間周波数帯に縞模様
が現われ、この画像データに基づいて再生した可視画像
が非常に見にくいものとなってしまうことがあった。

これまでこの縞模様を取り除く方法として、種々の方法
が試みられているが汎用性に乏しかった。

たとえば特開昭６２−１８１５７８号に記載された方法
等がある。この方法は、画像データの圧縮処理を行なう
に際し、原画像のデータ数を減らす前処理を行なうブロ
ックを各ブロックライン毎に位相をずらすようにしたも
のであるが、この方法は画像圧縮処理の際にしか用いる
ことができず、しかもこの方法は一方向のみにしか効果
がないため、この方法を用いる際は使用するグリッドの
方向（記録シート上に形成される縞模様の方向）を制限
する必要があった。本発明は、上記事情に鑑みグリッド
を使用して撮影記録された放射線画像から得られた画像
データから、該グリッドに起因するノイズを低減または
除去することのできる汎用性のある画像処理方法を提供
することを目的とするものである。

（課題を解決するための手段）本発明の第一の画像処理方法は、グリッドを使用した撮影により得られた該グリッドに対
応する縞模様を含む画像を読み取って得た画像データに
フーリエ変換を施し、このフーリエ変換により得られた
周波数データから前記縞模様に対応する周波数データを
低減または除去し、この低減または除去の行なわれた周
波数データに逆フーリエ変換を施すことを特徴とするも
のである。

また、本発明の第二の画像処理方法は、グリッドを使用
した撮影により得られた該グリッドに対応する縞模様を
含む画像を読み取って得た画像データに、前記縞模様の
空間周波数成分を低減または除去す実空間フィルタリン
グ処理を施すことを特徴とするものである。

（作　、用）本発明の第一の画像処理方法は、画像データにフーリエ
変換を施し、グリッドに起因する縞模様に対応する周波
数データを低減または除去した後、逆フーリエ変換を施
すようにしたため、逆フーリエ変換後の画像データはグ
リッドに起因する縞模様のノイズが低減または除去され
たものとなり、したがってこの画像データに基づいて、
上記縞模様の低減または除去された観察のしやすい可視
画像を得ることができる。

また、本発明の第二の画像処理方法は、画像データに、
グリッドに起因する縞模様の空間周波数成分を低減また
は除去する実空間フィルタリング処理を施すようにした
ため、上記第一の画像処理方法と同様に、上記縞模様の
低減または除去された観察のしやすい可視画像を得るこ
とができる。

（実　施　例）以下、本発明の実施例について、図面を参照して説明す
る。

第１図は、放射線画像撮影装置の一例の概略を示した図
である。ここでは記録シートとして前述した蓄積性蛍光
体シートを用いる例について説明する。

放射線源１から放射された放射線２は、被写体３を経由
して、さらにグリッド４を経由して蓄積性蛍光体シート
１１を照射する。グリッド４は４本／麿のピッチで鉛４
ａとアルミニウム４ｂとが交互に配置されている。放射
線２は、鉛４ａには遮ぎられ、アルミニウム４ｂは透過
してシート１１を照射する。

このためシート１１には被写体３の放射線画像とともに
４本／ｍの縞模様が蓄積記録される。被写体３内で散乱
された放射線２ａはグリッド４に斜めに入射するため該
グリッドに遮られ、または該グリッド４により反射され
、シート１１には照射されず、したがってシート１１に
は散乱放射線の照射の少ない鮮鋭な放射線画像が蓄積記
録される。

第２図は、グリッドを使用して撮影を行なうことにより
蓄積性蛍光体シート１１に蓄積記録された、被写体像（
図の斜線部）にグリッドに起因する縞模様（図の縦縞）
が重畳された放射線画像を示した図である。このように
シート１１には被写体像５と縞模様６とが重畳された放
射線画像が記録される。

第３図は、放射線画像読取装置の一例の斜視図である。

所定位置にセットされた放射線画像が記録された蓄積性
蛍光体シート１１は、図示しない駆動手段により駆動さ
れるエンドレスベルト等のシート搬送手段１５により、
矢印Ｙ方向に搬送（副走査）される。一方、レーザー光
源１６から発せられた光ビーム１７はモータ２４により
駆動され矢印方向に高速回転する回転多面鏡１８によっ
て反射偏向され、ｆθレンズ等の集束レンズ１９を通過
した後、ミラー２０により光路を変えて前記シート１１
に入射し副走査の方向（矢印Ｙ方向）と略垂直な矢印Ｘ
方向に主走査する。光ビーム１７が照射されたシート１
１の箇所からは、蓄積記録されている放射線画像情報に
応じた光量の輝尽発光光２１が発散され、この輝尽発光
光２１は光ガイド２２によって導かれ、フォトマルチプ
ライヤ（光電子増倍管）２３によって光電的に検出され
る。上記光ガイド２２はアクリル板等の導光性材料を成
形して作られたものであり、直線状をなす入射端面２２
ａが蓄積性蛍光体シート１１上の主走査線に沿って延び
るように配され、円環状に形成された射出端面２２ｂに
フォトマルチプライヤ２３の受光面が結合されている。

入射端面２２ａから光ガイド２２内に入射した輝尽発光
光２１は、該光ガイド２２の内部を全反射を繰り返して
進み、射出端面２２ｂから射出してフォトマルチプライ
ヤ２３に受光され、放射線画像を表わす輝尽発光光２１
がフォトマルチプライヤ２３によって電気信号に変換さ
れる。

フォトマルチプライヤ２３から出力されたアナログ出力
信号Ｓには、良好な放＃、［ｌ可視画像を再生出力する
ために必要な空間周波数帯のうちの最高の第１の空間周
波数ｆ　ｓｓ　−２，５ｃｙｃｌｅ／ｍより高い空間周
波数帯域の情報、特に第２図に示す縞模様６の情報も含
まれている。この縞模様の情報は可視画像を観察する際
にその可視画像を見にくくする原因のひとつとなるもの
であり、取り除く必要のある情報である。

アナログ出力信号Ｓはログアンプ２Ｂで対数的に増幅さ
れた後、Ａ／Ｄ変換器２８で所定のサンプリング間隔で
サンプリングされてディジタル化され、ディジタルの画
像データＳＤが得られる。この画像データＳＤは一旦記
憶部２９に記憶された後、画像処理装置５０に送られる
。

以下、画像処理装置５０で実施される、グリッド４に起
因する縞模様状のノイズの低減または除去の方法につい
て説明する。

第４図は、画像処理装置５０で実施される本発明の画像
処理方法の一例を説明するための図である。

実線のグラフ３３は、縞模様６（第２図参照）と直交す
る方向の放射線画像の空間周波数特性を表わした図であ
る。この放射線画像の空間周波数は高周波側は５ｃｙｃ
ｌｅ　／ｍｒｓまで延びており、４Ｃｙｅｌｅ／ｒｒａ
の位置に縞模様に対応するビーク３３ａがある。破線の
グラフ３４は第３図の放射線画像読取装置において、こ
の放射線画像を０．１４ｍのサンプリング間隔で読み取
ったときのエリアジングを模式的に示したもの、−点鎖
線で示すグラフ３５は、上記０．１４ａ＋ｍのサンプリ
ング間隔で読み取って得られた画像データＳＤが担持す
る放射線画像の空間周波数特性を示したものである。縞
模様のピーク３３ａのエリアジングが約３ｃｙｃｌｅ　
／ｓの位置に発生している。

ところでこの画像データＳＤに基づいて可視画像を再生
するに必要な最高の空間周波数は２．５ｃｙｃＩｅ／ｍ
である。そこで、画像データＳＤにフーリエ変換を施し
、図に二点鎖線で示すように、２，５ｃｙｃ　Ｉ　ｅ／
　ｒｒｍ以下の空間周波数成分を残し、それ以上の空間
周波数成分が除去される。その後逆フーリエ変換を施す
ことにより、必要な空間周波数成分を残すとともに、グ
リッドに起因した縞模様をなくすことができる。

第５Ａ図、第５Ｂ図は、本発明の画像処理方法の他の例
を説明するための図である。

実線で示すグラフ３３は第４図のグラフ３３と同一のも
のであり、放射線画像の縞模様を横切る方向の空間周波
数成分を表わしている。破線のグラフ３７はエリアジン
グを模式的に示したもの、−点鎖線のグラフ３８は読み
取られた画像データＳＤが担持する画像の空間周波数成
分を示したものである。

この例は、０．２　ｒｔａのサンプリング間隔で放射線
画像を読み取った場合を示しており、したがって１ｃｙ
ｃｌｅ　／ｙｍの位置にピーク３３ａのエリアジングが
発生している。

良好な可視画像を再生するためには、上記実施例（第４
図）の場合と同様に２．５ｃｙｃｌｅ／ｍ＋＋までの情
報が必要であるため、上記実施例と同様な方法を採用す
ることはできない。そこでここでは以下の方法が実施さ
れる。

第５Ｂ図は、第５Ａ図のグラフ３８（画像データＳ、が
担持する放射線画像の空間周波数特性を表わすグラフ）
を横に拡大して示した図である。

まず画像データにフーリエ変換を施して、図に示すよう
な周波数データが生成される。その後、各空間周波数ｆ
ｏに対し、その前後にΔｆｏの幅を設定し、ｆｏ−Δｆ
ｏ％ｆｏ＋Δｆｏの間の空間周波数成分の強度Ｐ（縦軸
）平均値Ｐを求め、ｆｏに対応する強度Ｐとの差Ｐ−Ｐ
が、所定のしきい値Ｔｈ１−　（１＋α）φ丁＋Ｔｈ２
−（１−α）・Ｐ（α：１より小さな正数）と比較され
、Ｐ＞ＴｈｌまたはＰ　＜　Ｔ　ｈ　２のときは、そこにグリッドに対応するピークがあるもの
としてその空間周波数ｆ、に対する強度Ｐがｐに置き換
えられる。図のピーク３８ａは上記のようにして図に破
線で示すグラフ３９に示すように取り除かれる。

このようにしてピーク３８ａを取り除いたあと逆フーリ
エ変換を施すことにより縞模様の低減された画像データ
を得ることができる。尚この実施例ではｆ、−Δｆｏ〜
ｆｏ＋Δｆ０の間の強度Ｐの平均値丁を用いているが、
これは平均値である必要はなくたとえばメジアン値等で
あってもよい。

また、しきい値も上記Ｔｈ１　、’ｒｈ２のように平均
値に基づいて定めるものではなくてもよい。

また他の実施例として、画像データをフーリエ変換した
後、各空間周波数ｆ０を中心としたｆ。

−Δｆｏ−ｆｏ＋Δｆｏの間の強度Ｐ（周波数データ）
のヒストグラムを求め、空間周波数ｆｏに対応する強度
Ｐがヒストグラムの端にある場合にその強度Ｐを平均値
Ｐやメジアン値等に置き換える方法を採用することもで
きる。

また本発明は、フーリエ変換を用いる方法に限られるも
のではない。

画像データから縞模様の情報を取り除くために実空間に
おいてフィルタリング処理を施すこともできる。

第６図は、蓄積性蛍光体シート上の各サンプリング点を
示した図である。図の・印が各サンプリング点を表わし
ており、サンプリング間隔Δは０．１　ｒｔｒｍである
。このサンプリング間隔の調整は、Ｘ方向（主走査方向
）についてはＡ／Ｄ変換器２８（第１図参照）のサンプ
リング時間間隔を調整すること等より行なわれ、Ｘ方向
（副走査方向）については、シート搬送手段１５（第１
図参照）によるシート１１の搬送速度を調整すること等
により行なわれる。このようにして得られた画像データ
は、サンプリング間隔ΔがΔ−０，１ｍｒｒｔであるた
め最高５ｃｙｃｌｅ　／ｒｒｍまでの情報を忠実に担持
しており、４ｃｙｃｌｅ／ｍにある縞模様の情報もその
まま担持している。この４ｃｙｃｌｅ／ｍの縞模様は可
視画像の再生に必要とする最高の空間周波数２．５ｃｙ
ｃｌｅ／間より高周波であるが、たとえば得られた画像
データの保存に必要なメモリ容量の削減を図るために画
像データにデータ圧縮処理を施す場合等において、第４
図、第５Ａ図を用いて説明したように低周波領域に４ｃ
ｙｃｌｅ／ｍの縞模様のエリアジングが生ずることが考
えられ、したがって上記縞模様が可視画像の再生に必要
とする最高の空間周波数２．５ｃｙｃｌｅ／ｍより高周
波であっても取り除いておく必要がある。

そこでここでは以下のようなフィルタリング処理が施さ
れる。

第４図に示すａ＋１　（１，ｊ　−１，２，・・・・・
・）は各サンプリング点のフィルタリング処理前の画像
データの値を表わし、ｂ＋、、（１，■−１，２．・・
・・・・）はフィルタリング処理後の画像データの値を
表わしている。ここでフィルタリング処理前の画像デー
タａ＋＋　（１，ｊ−１，２，・・・・・・）からフィ
ルタリング処理後の画像データｂ　Ｉｓ　（１，ｓ　−
１、２，・・・・・・）が、ｂ、− ６００Ｘ　（１２１ａ＋１”１９８”　ａ＋＋１、ｊ”１２１
−ａ＋＋２．　ｒ＋１９８−ａ＋、　Ｈ＋＋　”３２４
Φａ＋＋＋、　１＋＋　”１９８−ａ＋＋２．　ｒ＋１
２１　・ａ＋、　ｌ＋２　”１９８”　　ａ＋＋＋、　
１４２　”１２１　”　ａ　ｌ＋２．１４２）　　　　
　　　　　　　　　　　　　　　　　・・・・・・（１
）の演算式に従って求められる。

（１）式の演算を各サンプリング点について施すことに
より、フィルタリング処理が施される。

第７図は、上記（１）式のフィルタリング特性を表わし
た図である。

グラフ３１で示される特性のフィルタ（上記（１）式）
は４　ｅｙｃｌｅ　／ａｍの伝達特性が零であり、した
がって（１）式に従ってフィルタリング処理を施すこと
により求められた画像データは、縞模様の情報を担持し
ないものとなる。したがってこのフィルタリング処理の
施された画像データにたとえばデータ圧縮処理等を施し
てもグリッドに起因した縞模様が発生することはない。

尚、このフィルタリング処理に用いる演算式は（１）式
に限られるものではなく、縞模様の空間周波数、可視画
像の再生に必要な空間周波数等によって種々の演算式の
中がら適宜選択され得るものである。

また、フィルタリングの方法として、特異的なピークが
あった場合に、このピークのデータ値を周囲の所定範囲
内のデータのメジアン値で置き換えるメジアンフィルタ
を用いてもよい。上記特異的なピークとして縞模様に照
準を合わせこの縞模様を取り除くようにメジアンフィル
タを適切に施すことにより、グリッドに起因した縞模様
を除去または低減することができる。

尚、これまで述べた各実施例においてグリッドの方向（
画像上の縞模様の方向）が不明の場合には、上記各フィ
ルタリングを二次元的に施してもよいが、グリッドの方
向（縞模様の方向）が既知の場合には、縞模様と直交す
る方向に一次元的なフィルタリングを施すとこのフィル
タリング処理に要する時間が短くて済み効率的である。

また、縞模様の方向が未知の場合、主走査方向（第３図
のＸ方向）に並んだ１ラインないし数ライン分の画像デ
ータと副走査方向（第３図のＹ方向）に並んだ１ライン
ないし数ライン分の画像データにそれぞれフーリエ変換
を施し、ピークの存在を探索することにより、このピー
クの存在により縞模様の方向を見い出すことができ、前
述したように一次元的なフィルタリングにより縞模様を
取り除くことができる。さらに、縞模様の空間周波数が
不明の場合、縞模様と交叉する方向に並んだ１ラインな
いし数ライン分の画像データをフーリエ変換しピークを
形成している空間周波数を探索することにより、このピ
ークの空間周波数により縞模様の空間周波数を知ること
ができ、これにより適切なフィルタを選択することがで
きる。

第３図に示す画像処理装置５０において、これまで説明
した方法のひとつを用いて、グリッドに起因するノイズ
の情報を画像データから除去または低減した後、該画像
データにその他の画像処理が施され、−旦記憶され、そ
の後図示しない画像再生装置に送られて該画像データに
基づく可視画像が再生される。

以上本発明の種々の実施例について説明したが、本発明
は上記実施例に限られるものではなく、周波数空間、実
空間における公知の各種タイルタをグリッドに対応した
縞模様の除去または低減に用いることができる。

（発明の効果）以上詳細に説明したように、本発明の画像処理方法は、
画像データに実空間または周波数空間におけるフィルタ
リングを施すことにより、該画像データからグリッドに
起因した縞模様の情報を除去または低減するようにした
ため、このフィルタリング処理の後、データ圧縮処理、
画像拡大縮小処理等の様々な処理を行なう場合において
もグリッドに起因した縞模様（アーティファクト）が生
ずることはない。また、上記フィルタリングを二次元的
に施す場合、又は縞模様の方向やピッチを自動的に認識
して縞模様と直交する方向に一次元的に施す場合は、撮
影の際にグリッドのピッチやその方向を気にする必要も
ない。

【図面の簡単な説明】

第１図は、放射線画像撮影装置の一例の概略を示した図
、第２図は、グリッドを使用して撮影を行なうことにより
、蓄積性蛍光体シートに蓄積記録された、被写体像にグ
リッドに起因する縞模様が重畳された画像を示した図、第３図は、放射線画像読取装置の一実施例の斜視図、第４図は、蓄積性蛍光体シートに記録された放射線画像
（縞模様を含む）の、縞模様と直交する方向の空間周波
数特性と、該画像を横方向に０．１４ｍ間隔でサンプリ
ングして得た画像データが担持する空間周波数特性を示
した図、第５Ａ図は、蓄積性蛍光体シートに記録された放射線画
像（縞模様を含む）の、縞模様と直交する方向の空間周
波数特性と、該画像を該方向に０．２＃間隔でサンプリ
ングして得た画像データが担持する空間周波数特性を示
した図、第５Ｂ図は、第５Ａ図の画像データが担持する空間周波
数特性を横方向に拡大して示した図、第６図は、蓄積性
蛍光体シート上の各サンプリング点を示した図、第７図は、フィルタ特性の一例を示した図である。１・・・放射線源　　　　３・・・被写体５・・・被写
体像１１・・・蓄積性蛍光体シート４・・・グリッド６・・・縞模様１７・・・光ビーム２３・・・フォトマルチプライヤ２６・・・ログアンプ　　　２８・・・Ａ／Ｄ変換器２
９・・・記憶部　　　　　５０・・・画像処理装置第１
図第２図こ龜

Claims

【特許請求の範囲】

（１）グリッドを使用した撮影により得られた該グリッ
ドに対応する縞模様を含む画像を読み取って得た画像デ
ータにフーリエ変換を施し、このフーリエ変換により得
られた周波数データから前記縞模様に対応する周波数デ
ータを低減または除去し、この低減または除去の行なわ
れた周波数データに逆フーリエ変換を施すことを特徴と
する画像処理方法。
（２）グリッドを使用した撮影により得られた該グリッ
ドに対応する縞模様を含む画像を読み取って得た画像デ
ータに、前記縞模様の空間周波数成分を低減または除去
す実空間フィルタリング処理を施すことを特徴とする画
像処理方法。