JPH0228782A - 被写体像内画像点決定方法 - Google Patents

被写体像内画像点決定方法

Info

Publication number
JPH0228782A
JPH0228782A JP1029179A JP2917989A JPH0228782A JP H0228782 A JPH0228782 A JP H0228782A JP 1029179 A JP1029179 A JP 1029179A JP 2917989 A JP2917989 A JP 2917989A JP H0228782 A JPH0228782 A JP H0228782A
Authority
JP
Japan
Prior art keywords
image
image signal
value
recording sheet
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.)
Granted
Application number
JP1029179A
Other languages
English (en)
Other versions
JP2631737B2 (ja
Inventor
Hideya Takeo
英哉 武尾
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujifilm Holdings Corp
Original Assignee
Fuji Photo Film Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fuji Photo Film Co Ltd filed Critical Fuji Photo Film Co Ltd
Priority to JP1029179A priority Critical patent/JP2631737B2/ja
Priority to EP89107156A priority patent/EP0340553B1/en
Priority to DE68921786T priority patent/DE68921786T2/de
Priority to US07/340,744 priority patent/US5179597A/en
Publication of JPH0228782A publication Critical patent/JPH0228782A/ja
Application granted granted Critical
Publication of JP2631737B2 publication Critical patent/JP2631737B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/66Analysis of geometric attributes of image moments or centre of gravity
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10116X-ray image

Abstract

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

Description

【発明の詳細な説明】 (産業上の利用分野) 本発明は、被写体像を含む放射線画像が記録された記録
シート上の各画素がら得られた画像信号に基づいて、上
記被写体像内の画像点を求める方法に関するものである
(従来の技術) 記録された放射線画像を読み取って画像信号を得、この
画像信号に適切な画像処理を施した後、画像を再生記録
することは種々の分野で行なわれている。たとえば、後
の画像処理に適合するように設計されたガンマ値の低い
X線フィルムヲ用いてX線画像を記録腰このX線画像が
記録されたフィルムからX線画像を読み取って電気信号
に変換腰この電気信号(画像信号)に画像処理を施した
後コピー写真等に可視像として再生することにより、コ
ントラスト、シャープネス、粒状性等の画質性能の良好
な再生画像を得ることのできるシステムが開発されてい
る(特公昭61−5193号公報参照)。
また本願出願人により、放射線(X線、α線。
β線、γ線、電子線、紫外線等)を照射するとこの放射
線エネルギーの一部か蓄積され、その後可視光等の励起
光を照射すると蓄積されたエネルギーに応じて輝尽発光
を示す蓄積性蛍光体(輝尽性蛍光体)を利用して、人体
等の被写体の放射線画像を一部シート状の蓄積性蛍光体
に撮影記録し、この蓄積性蛍光体シートをレーザー光等
の励起光で走査して輝尽発光光を生ぜしめ、得られた輝
尽発光光を光電的に読み取って画像信号を得、この画像
信号に基づき被写体の放射線画像を写真感光材料等の記
録材料、CRT等に可視像として出力させる放射線画像
記録再生システムかすでに提案されている(特開昭55
−12429号、同5B−11395号。
同55−183472号、同5B−LO4[i45号、
同55−116340号等)。
このシステムは、従来の銀塩写真を用いる放射線写真シ
ステムと比較して極めて広い放射線露出域にわたって画
像を記録しうるという実用的な利点を有している。すな
わち、蓄積性蛍光体においては、放射線露光量に対して
蓄積後に励起によって輝尽発光する発光光の光量が極め
て広い範囲にわたって比例することが認められており、
従って種々の撮影条件により放射線露光量がかなり大幅
に変動しても、蓄積性蛍光体シートより放射される輝尽
発光光の光量を読取ゲインを適当な値に設定して光電変
換手段により読み取って電気信号に変換し、この電気信
号を用いて写真感光材料等の記録材料、CRT等の表示
装置に放射線画像を可視像として出力させることによっ
て、放射線露光量の変動に影響されない放射線画像を得
ることができる。
上記システムにおいて、蓄積性蛍光体シートに照射され
た放射線の線量等に応じて最適な読取条件で読み取って
画像信号を得る前に、予め低レベルの光ビームにより蓄
積性蛍光体シートを走査してこのシートに記録された放
射線画像の概略を読み取る先読みを行ない、この先読み
により得られた先読画像信号を分析し、その後上記シー
トに上記先読みの際の光ビームよりも高レベルの光ビム
を照射して走査し、この放射線画像に最適な読取条件で
読み取って画像信号を得る本読みを行なうように構成さ
れたシステムもある(特開昭58−67240号、同5
g−67241号、同58−67242号等)。
ここで読取条件とは、読取りにおける輝尽発光光の光量
と読取装置の出力との関係に影響を与える各種の条件を
総称するものであり、例えば入出力の関係を定める読取
ゲイン、スケールファクタあるいは、読取りにおける励
起光のパワー等を意味するものである。
また、光ビームの高レベル/低レベルとは、それぞれ、
上記シートの単位面積当りに照射される光ビームの強度
の大/小、もしくは上記シートから発せられる輝尽発光
光の強度が上記光ビームの波長に依存する(波長感度分
布を有する)場合は、上記シートの単位面積当りに照射
される光ビームの強度を上記波長感度で重みづけした後
の重みづけ強度の大/小をいい、光ビームのレベルを変
える方法としては、異なる波長の光ビームを用いる方法
、レーザ光源等から発せられる光ビームの強度そのもの
を変える方法、光ビームの光路上にNDフィルター等を
挿入、除去することにより光ビムの強度を変える方法、
光ビームのビーム径を変えて走査密度を変える方法、走
査速度を変える方法等、公知の種々の方法を用いること
ができる。
また、この先読みを行なうシステムか先読みを行なわな
いシステムかによらず、得られた画像信号(先読画像信
号を含む)を分析し、画像信号に画像処理を施す際の最
適な画像処理条件を決定するようにしたシステムもある
。この画像信号に基づいて最適な画像処理条件を決定す
る方法は、蓄積性蛍光体シートを用いるシステムに限ら
れず、たとえば従来のX線フィルム等の記録シートに記
録された放射線画像から画像信号を得るシステムにも適
用されている。
上記画像信号(先読画像信号を含む)を分析して最適な
読取条件、画像処理条件を求める方法は種々提案されて
いるが、その方法のひとつとして、画像信号のヒストグ
ラムを作成する方法か知られている(たとえば、特願昭
59−12658号)。画像信号のヒストグラムを求め
ることにより、たとえば画像信号の最大値、最小値や、
頻度が最大となる点の画像信号の値等を知ることができ
、これらの各個から蓄積性蛍光体シート、X線フィルム
等の記録シートに記録された放射線画像の特徴を把握す
ることができる。そこでこのヒストグラムに基づいて最
適な読取条件1画像処理条件を求めることにより、観察
適正のすぐれた放射線画像を再生出力することが可能と
なる。
一方、記録シートに放射線画像を撮影記録するに際して
は、被写体の観察に必要の無い部分に放射線を照射しな
いようにするため、あるいは観察に不要な部分に放射線
を照射するとその部分から観察に必要な部分に散乱線が
入り画質性能が低下するため、放射線が被写体の必要な
部分および記録シートの一部にのみ照射されるように放
射線の照射域を制卜艮する照射野絞りを使用して撮影を
行なうことも多い。
ところが、前述のようにして画像信号を分析して読取条
件1画像処理条件を求めるにあたって、分析に用いた画
像信号が、照射野絞りを用いて撮影した記録シートから
得られた画像信号である場合、この照射野の存在を無視
して画像信号を分析しても撮影記録された放射線画像か
正しく把握されず、誤った読取条件、画像処理条件が求
められ観察適正の優れた放射線画像が再生記録されない
場合が生ずる。
これを解決するためには、読取条件1画像処理条件を求
める前に、照射野を認識し、照射野内の画像信号に基づ
いて読取条件1画像処理条件を求める必要がある。
照射野を認識する方法のうち、放射線照射野が不規則な
形状をしていても正確に照射野を認識することのできる
汎用性のある方法としては、例えば、照射野内に含まれ
る所定の点からシート端部に向かう放射状の複数の方向
に沿った各画素に対応する画像信号に基づいて、照射野
の輪郭上にあると考えられる複数の輪郭点を求め、これ
らの輪郭点に沿った線で囲まれる領域を照射野と認識す
る方法が、本出願人により既に提案されている(特願昭
82−93633号)。
この方法において、上記所定の点は照射野内の点である
必要があり、さらには照射野内のうちの被写体像内の点
であることが望ましい。ただし、照射野絞りは被写体の
うちの必要な部分のみを撮影するために使用するもので
あるため、照射野絞りを使用した撮影では、照射野内(
特に照射野の略中心)の画像点はほとんどの場合被写体
像内の画像点となる。
(発明が解決しようとする課題) 上記のようにしてまず放射線照射野を求め、その後求め
られた放射線照射野内に対応する画像信号を分析するこ
とにより、適切な読取条件1画像処理条件が求められる
しかし、上記特願昭62−93633号で提案した方法
で照射野を認識するためには、その前提として照射野内
(望ましくは被写体像内)の画像点を見つけ出すことが
必要となる。
記録シートの中心点を上記照射野内(被写体像内)の画
像点とする方法も考えられる。例えば第5A図に示すよ
うに、記録シート1の中央に照射野2かある場合は、記
録シート1の中心点Cは照射野2内の画像点となり、前
述したように、はとんどの場合被写体像3内の画像点と
なる。
しかし、第5B図に示すように照射野が記録シート1の
中心Cから外れている場合は、この方法では、照射野2
内の画像点を求めることはできない。また、照射野絞り
を使用しないで撮影を行なった場合は、第5C図に示す
ように記録シート1の中心点を含む全点か照射野2内の
画像点ではあるが、この場合照射野絞りが使用されてお
らず、被写体像3が照射野2の略中心(この場合は記録
シート1の中心Cと同一)にない確率は、第5A図、第
5B図のように照射野絞りを使用して撮影を行なった場
合と比べて高く、したがって記録シートの中心点Cが被
写体像内の画像点ではない場合も多い。
また、被写体像内の画像点を求めることは、上記のよう
に照射野を認識することのほか、放射線画像の一部を例
えば拡大して可視画像として再生出力する場合に、放射
線画像のどの部分を中心とした可視画像を出力するかを
決定する際にも有効である。
本発明は、上記事情に鑑み、照射野絞りを使用しもしく
は使用しないで撮影することにより得られた被写体像を
含む放射線画像の、被写体像内の画像点を求める方法を
提供することを目的とするものである。
(課題を解決するための手段) 本発明の被写体像白画像点決定方法のひとつは、被写体
像を含む放射線画像が記録された、蓄積性蛍光体シート
、写真フィルム等の記録シート上の各画素から得られた
画像信号に基づいて、各画素に対応する画像信号値また
はこの逆数でそれぞれ対応する各画素を重み付けして記
録シートの重心を求め、この重心を、被写体像内の画像
点とすることを特徴とするものである。
また、本発明の被写体像白画像点決定方法の他のひとつ
は、上記と同様にして画像信号を得た後、この画像信号
に基づいて、各画素に対応する画像信号値またはこの画
像信号値の逆数をそれぞれ対応する各画素に対応づけた
ときに、記録シート上の互いに異なる2つの方向のそれ
ぞれについて、画像信号値またはこの逆数を上記各方向
に累積してプロットした累積分布を求め、これらの累積
分布のそれぞれについて最大の累積値の略半分の値に対
応する各方向についての座標点を求め、これらの座標点
により定まる記録シート上の位置を、被写体像内の画像
点とすることを特徴とするものである。
また、上記2つの方法において、画像信号値またはこの
逆数のいずれを用いるかについては、以下の方法を用い
ることができる。すなわち、画像信号を得た後、この画
像信号に基づいて、記録シートの周辺部に対応する画像
信号値を代表する第1の代表値と、記録シートの全部ま
たは略中央部に対応する画像信号値を代表する第2の代
表値とを求め、これら第1の代表値と第2の代表値との
大小を比較し、この比較結果に応じて画像信号値または
この逆数の一方を選択することができる。
ここで、上記「記録シート上の各画素から得られた画像
信号」には、たとえば蓄積性蛍光体シートから発せられ
た輝尽発光光を光電的に読み取って得た画像信号や写真
フィルムを透過しまたは写真フィルムから反射された光
を光電的に読み取って得た画像信号等が含まれるが、か
ならずしもこれらに限られるものではない。
また、上記「画像信号値」は、上記画像信号の各画素毎
の値をいう。
また、「記録シ、−トの周辺部に対応する画像信号値を
代表する第1の代表値」は、かならずしも記録シートの
周辺部に対応する画像信号全てに基づく必要はなく、周
辺部に対応する画像信号値を代表しているとみなし得る
範囲で変動があってもよい。すなわち周辺部に対応する
画像信号値を代表しているとみなし得る範囲で、周辺部
の一部に対応する画像信号に基づいて第1の代表値を求
めてもよい。
また、上記「記録シートの全部または略中央部に対応す
る画像信号値を代表する第2の代表値」は、上記周辺部
と同様に、かならずしも記録シート全体または略中央部
に対応する画像信号全てに基づく必要はなく、記録シー
トの全部または略中央部に対応する画像信号値を代表し
ているとみなし得る範囲で変動かあってもよい。すなわ
ち、記録シートの略全領域または略中央付近のまわりに
所定の広がりをもつ領域に対応する画像信号に基づいて
第2の代表値を求めてもよい。
また、上記第1の代表値、第2の代表値としては、対応
する画像信号のたとえば平均値、中央値。
(最大値十最小値)/21画像信号の値の小さいものか
らその画像信号の数を加算してプロットした累積ヒスト
グラム(後述する第2C図のグラフB参照)から求めた
所定の累積値に対応する画像信号の値等が用いられるが
、第1の代表値と第2の代表値は同一の演算方法により
求めたもの(たとえば両者とも平均値)でなくともよく
、それぞれ適切な代表値が選択される。
(作  用) 本発明の被写体像内画像点決定方法のひとつは、各画素
に対応する画像信号値またはこの画像信号値の逆数でそ
れぞれ対応する各画素を重み付けして記録シートの重心
を求めるものであるため、この重心が実用上十分な高確
率をもって被写体像内の画像点として求められる。
また、本発明の他の被写体像内画像点決定方法は、各画
素に対応する画像信号値またはこの画像信号値の逆数を
それぞれ対応する各画素に対応づけたときに、記録シー
ト上の互いに異なる2つの方向のそれぞれについて上記
画像信号値または上記逆数を上記各方向に累積してプロ
ットした累積分布を求め、これらの累積分布のそれぞれ
について最大の累積値の略半分の値に対応する各方向に
ついての座標点を求めるものであるため、これらの座標
点により定まる記録シート上の位置が、上記重心と同様
に、実用上十分な高確率をもって被写体像内の画像点と
して求められる。
ここで、上記2つの方法において、上記画像信号値また
は上記逆数のいずれを用いるかについては、撮影の際に
照射野絞りが使用され記録シートのうちの狭い範囲にの
み放射線画像が記録されているか、あるいは照射野絞り
を使用せずにまたは使用しても記録シートの広い範囲に
わたって放射線画像が記録されているか、また、記録シ
ートに記録された放射線画像がいわゆるネガ像かポジ像
かによって異なる。このため、これらの情報を直接用い
て、撮影または読取りの都度上記画像信号値または上記
逆数の一方を選択することも考えられるが、上記ネガ像
又はポジ像の一方のみを取り扱うシステム、または記録
シートとして蓄積性蛍光体シートを用いるシステム等、
ネガ、ポジの区別のないシステムなどにおいては、記録
シートの周辺部に対応する画像信号値を代表する第1の
代表値と、記録シートの全部または略中央部に対応する
画像信号値を代表する第2の代表値とを求め、これら第
1の代表値と第2の代表値との大小を比較することによ
り、照射野絞りが使用され記録シートのうちの狭い範囲
にのみ放射線画像が記録されているか、あるいは照射野
絞りを使用せずにまたは使用しても記録シートの広い範
囲にわたって放射線画像が記録されているかを判定する
ことができ、この比較結果に応じて上記画像信号値また
は上記逆数の一方を選択して、この選択された上記画像
信号値または上記逆数を用いて上記いずれかの方法によ
り被写体像内の画像点を定めることができる。
(実 施 例) 以下、本発明の実施例について、図面を参照して説明す
る。
第1図は、本発明の被写体像内画像点決定方法の一例を
使用した、放射線画像読取再生装置の一実施例の斜視図
である。この放射線画像読取再生装置は前述した蓄積性
蛍光体シートを用いる装置である。
放射線画像が記録された蓄積性蛍光体シート11は、読
取手段100の所定位置にセットされる。この所定位置
にセットされた蓄積性蛍光体シート11は、図示しない
駆動手段により駆動されるエンドレスベルト等のシート
搬送手段15により、矢印Y方向に搬送(副走査)され
る。一方、レーザー光源16から発せられた光ビーム1
7はモータ24により駆動され矢印方向に高速回転する
回転多面鏡18によって反射偏向され、fθレンズ等の
集束レンズ19を通過した後、ミラー20により光路を
変えて前記シート11に入射し副走査の方向(矢印Y方
向)と略垂直な矢印X方向に主走査する。光ビーム17
が照射されたシート11の箇所からは、蓄積記録されて
いる放射線画像情報に応じた光量の輝尽発光光21が発
散され、この輝尽発光光21は光ガイド22によって導
かれ、フォトマルチプライヤ(光電子増倍管)23によ
って光電的に検出される。上記光ガイド22はアクリル
板等の導光性材料を成形して作られたものであり、直線
状をなす入射端面22aが蓄積性蛍光体シート11上の
主走査線に沿って延びるように配され、円環状に形成さ
れた出射端面22bにフォトマルチプライヤ23の受光
面が結合されている。入射端面22aから光ガイド22
内に入射した輝尽発光光21は、該光ガイド22の内部
を全反射を繰り返して進み、出射端面22bから出射し
てフォトマルチプライヤ23に受光され、放射線画像を
表わす輝尽発光光21の光量がフォトマルチプライヤ2
3によって電気信号に変換される。
フォトマルチプライヤ23から出力されたアナログ出力
信号Sはログアンプ26で対数的に増幅され、A/D変
換器27でディジタル化され、画像信号S。が得られる
。得られた画像信号SQは一旦記憶手段28に記憶され
た後、画像処理手段29によって読み出される。
画像処理手段29では、被写体像内の画像点が求められ
、この画像点から蓄積性蛍光体シート1■の端部に向か
う放射状の複数の方向に沿った各画素に対応する画像信
号に基づいて、照射野の輪郭上にあると考えられる複数
の輪郭点が求められ、これらの輪郭点に沿った線で囲ま
れる領域が照射野として認識される。照射野が認識され
ると、この照射野に対応する画像信号S0に適切な画像
処理が施される。
画像処理の施された画像信号SOは再生手段30に送ら
れ、再生手段30ではこの画像信号SQに基づく放射線
画像が再生記録される。
ここで、画像処理手段29で行なわれる、被写体像内の
画像点を画像信号SQに基づいて決定する方法について
さらに説明する。
画像処理手段29では、画像信号値So  (簡単のた
め画像信号値を表わす場合も画像信号を表わす場合と同
様に80を用いる。)またはこの画像信号値の逆数17
S0でそれぞれ対応する各画素を重み付けして、蓄積性
蛍光体シート11の重心が求められ、この重心が被写体
像内の画像点とされる。
この演算を行なうために、まず上記画像信号値SQを用
いて演算を行なうか、それとも上記逆数1/S0を用い
て演算を行なうかを決定する必要がある。この実施例に
示すシステムでは、前述したネガ像、ポジ像の区別はな
く、蓄積性蛍光体シート11に蓄積された放射線エネル
ギーに略比例する光量の輝尽発光光が得られるものであ
るため、上記画像信号値S0を用いるかまたは上記逆数
1/Soを用いるかは、撮影に際して照射野絞りが使用
され、蓄積性蛍光体シート11のうちの狭い範囲にのみ
放射線画像が記録されているか、あるいは照射野絞りを
使用せずにまたは使用しても蓄積性蛍光体シート11の
広い範囲にわたって放射線画像が記録されているかによ
り定められる。この照射野に関する情報は、たとえば撮
影の都度、図示しないキーボード等の入力手段により画
像処理手段29に入力してもよいが、本実施例では、以
下に述べる方法によりこの装置自身で上記画像信号値S
0を用いるか上記逆数1/S0を用いるかの判定が行な
われる。
第2A図は、照射野絞りを使用して撮影され、蓄積性蛍
光体シート11上の狭い領域にのみ放射線が照射された
様子を表わす図、第2B図は、照射野絞りを使用しない
で撮影され、蓄積性蛍光体シート全面にわたって放射線
が照射された様子を表した図である。
第2A図に示す照射野12外の領域は、放射線はほとん
ど照射されていないため、この領域の各画素から得られ
た画像信号SOは、非常に小さな値を有している。第2
A図、第2B図に示す照射野12のうち、被写体像13
以外の領域は、放射線が被写体を経由(透過又は反射)
せずに、直接蓄積性蛍光体シート11に照射された領域
であり、この領域の各画素から得られた画像信号SQは
非常に大きな値を有している。被写体像13に対応する
各画素から得られた画像信号S0は、被写体を経由した
後の放射線が蓄積性蛍光体シートに照射された領域に対
応するものであるため、中間的な値を有している。
第2A図、第2B図に示す周辺部14の各画素に対応す
る画像信号に基づいて、この領域14を代表する第1の
代表値か求められる。この第1の代表値としては、領域
14に対応する画像信号の平均値か採用される。また、
蓄積性蛍光体シート11の全体に対応する画像信号に基
づいて、このシート11の放射線画像の全体を代表する
第2の代表値が求められる。
第2C図を参照してこの第2の代表値の求め方について
説明する。
第2C図は、シート11の全面に対応する画像信号値S
0のヒストグラム(グラフA)と累積ヒストグラム(グ
ラフB)の−例を示したグラフである。横軸は画像信号
値S0を示しており、縦軸は、グラフAについては画像
信号SOの各個の発生頻度、グラフBについては、グラ
フAの頻度を順次累積した累積値(%)を示している。
このように、第2の代表値を求めるために第2C図に示
すような累積ヒストグラムが求められ、この累積ヒスト
グラムから累積値(%)30%に対応する画像信号値S
30が求められ、このS30が第2の代表値とされる。
このようにして第1の代表値と第2の代表値が求められ
ると、この2つの代表値の大小が判断される。第1の代
表値は、前述したように領域14に対応する画像信号S
0の平均値である。第2A図に示す放射線画像の場合は
、領域14は照射野絞りにより放射線が遮ぎられた領域
であるため、この領域14に対応する画像信号の平均値
は小さな値を示す。また、この領域14は第2B図に示
す放射線画像の場合は、被写体を経由(透過又は反射)
せずに直接放射線が照射された領域であるため、この領
域14に対応する画像信号の平均値は大きな値を示す。
また、第2の代表値は前述のようにして、第2A図に対
応する第1の代表値と第2B図に対応する第1の代表値
との略中夫の値が選ばれている。
したがって、上記第1の代表値と第2の代表値とを比較
し、第1の代表値く第2の代表値の場合は、撮影に際し
て照射野絞りが使用され、蓄積性蛍光体シート1のうち
の狭い範囲にのみ放射線画像が記録されている場合に対
応し、この場合は後述するように画像信号値SQを用い
て被写体像内の画像点を求めるべき旨の決定が行なわれ
る。また、第1の代表値〉第2の代表値の場合は、撮影
に際して照射野絞りが使用されていないかもしくは使用
されていても蓄積性蛍光体シート11の広い範囲にわた
って放射線画像が記録されている場合に対応し、この場
合は後述するように画像信号値S0の逆数1/SQを用
いて被写体像内の画像点を求めるべき旨の決定がなされ
る。
第3A図、第3B図は、第2A図、第2B図にそれぞれ
対応した画像信号値S0、逆数1/S0のグラフの一例
を放射線画像とともに示す図である。
第3A図は、第2A図に対応し、蓄積性蛍光体シート1
1のうち狭い領域にのみ照射野12が形成された放射線
画像を示している。この場合に直線に沿った各画素に対
応する画像信号値S0をプロットすると、グラフCに示
すように、照射野12内の各画素に対応する画像信号が
大きな値を有する。
したがって、この場合、画像信号値S。で対応する各画
素を重み付けして重心を求めることにより、この重心は
照射野12内の画像点として求められ、かつ、前述した
ように照射野絞りは被写体の必要な部分のみを撮影する
目的で使用するものであるため、この重心は実用上十分
な精度をもって被写体像13内の画像点とに求められる
第3B図は、第2B図に対応し、蓄積性蛍光体シート1
1の全面が照射野12に対応した放射線画像を表わして
いる。この場合に、直線に沿った各画素に対応する画像
信号値SQの逆数1/s0をプロットすると、グラフD
に示すように、被写体像内の各画素に対応する上記逆数
1/SQが大きな値を有する。したがって、この場合、
上記逆数1/ S oで対応する各画素を重み付けして
重心を求めることにより、この重心が被写体像内の画像
点として求められる。
ここで、上記画像信号値S。またはその逆数1/SOで
対応する各画素を重み付けして重心を求めることは、以
下に示す演算を行なうことをいう。
X方向に並ぶ画素の座標を] (正の整数)、X方向に
並ぶ画素の座標をj (正の整数)で表わし、各画素(
i、j)に対応する上記画像信号値SQまたはその逆数
1/S0をpQ、j)で表わす。また、Px(j)=Σ
P(i、j)      ・・・・・・(1)P y 
(i)−ΣP(i、j)      ・・・・・・(2
)で表わし、X方向、X方向の重心をXc+V。とする
と、 Σi  −Py(i) Σj −PX(j) により、重心の座標(xe、ye)が求められる。
また、第1図に示す画像処理手段29で行なわれる演算
は、前述した重心を求める演算に代えて、上記画像信号
値S0またはその逆数1/S0でそれぞれ対応する前記
各画素に対応づけたときに、前記記録シート上の互いに
異なる2つの方向(ここではX方向とX方向)のそれぞ
れについて、画像信号値Soまたはその逆数1/S0を
X方向。
X方向についてそれぞれ累積してプロットした累積分布
を求め、これらの累積分布のそれぞれについて最大の累
積値の略半分の値に対応するX方向。
X方向についての座標点Xc  、yc’を求めるもの
であってもよい。この場合、これらの座標点xc 、y
。′により定まる蓄積性蛍光体シート11上の位置(x
c  +ye’)が、被写体像内の画像点として求めら
れる。
第4図は、この演算を説明するための、x 方向。
X方向の画像信号の値の分布およびそれらの分布の累積
分布の一例を表わした図である。
前述した方法により選択された画像信号値sQまたはそ
の逆数1 / S oを各画素(i、j)に対応させて
P (i、j)で表わす。
このとき、グラフEは、 Py(i)=Σp(i、j) を、X方向(i方向)にプロットしたグラフ、グラフG
は、 px(j)=ΣP(i、j) をX方向(j方向)にプロットしたグラフを表わしてお
り、グラフF、グラフHは、それぞれグラフE、グラフ
Gに示すPy(i)、Px(j)をそれぞれX方向、X
方向に累積した値(%)、を表わしている。すなわち、
グラフF、グラフGはそれそ”れ を示している。
ここで、X方向、X方向の座標点X。
′がそれぞれ C として求められる。これらの座標点xc 、ye′は、
第4図におけるグラフF、グラフHが50%に達する点
Xc 、yc′である。
このようにして座標点(xe  rye’)を求めるこ
とにより、この座標点(xc’+ye’)が、上記重心
を求めた場合と同様、被写体像内の画像点となる。
このようにして、被写体像内の画像点が求められた後、
前述したように、この画像点から蓄積性蛍光体シート1
1の端部に向かう放射状の複数の方向に沿った各画素に
対応する画像信号に基づいて、照射野の輪郭上にあると
考えられる複数の輪郭点が求められ、これらの輪郭点に
沿った線で囲まれる領域が照射野として認識され、この
照射野に対応する画像信号SOに適切な画像処理が施さ
れる。
上記実施例においては、求められた被写体像白画像点を
、照射野の認識に用いる場合について説明したが、被写
体像白画像点を求めることは、照射野の認識にのみ用い
られるものではない。たとえば、第1図に示す再生手段
30で可視画像を得る場合に、上記のようにして画像処
理手段29において被写体像白画像点を求めることによ
り、被写体像が蓄積性蛍光体シート11の端に記録され
た場合であっても、この被写体像を中央に配置した可視
画像を得ることができる。
また、上記実施例は、先読みを行なわない放射線画像読
取装置について説明したが、先読みを行なって先読画像
信号を求め、この先読画像信号に基づいて照射野を求め
、この照射野内に対応する先読画像信号に基づいて本読
みの際の読取条件を求めるシステムにも本発明の被写体
像白画像点決定方法を用いることができることはいうま
でもない。
また、本発明は、蓄積性蛍光体シートを用いる装置のほ
か、従来のX線フィルムを用いる装置等にも用いること
ができる。
第6図は、X線フィルムに記録されたX線画像を読み取
るX線画像読取装置の一実施例の斜視図である。
所定位置にセットされた、被写体像を含むX線画像が記
録されたX線フィルム40がフィルム搬送手段41によ
り図に示す矢印Y′方向に搬送される。
また、−次元的に長く延びた光源42から発せられた読
取光43は、シリンドリカルレンズ44により収束され
、X線フィルム上を矢印Y′方向と略直角なX′力方向
直線状に照射する。読取光43が照射されたX線フィル
ム40の下方には、X線フィルム40を透過しX線フィ
ルム40に記録されたX線画像により強度変調された読
取光43を受光する位置に、上記X線画像のX′力方向
各画素間隔に対応した多数の固体光電変換素子が直線状
に配置されたMOSセンサ45が設けられている。この
MOSセンサ45は、X線フィルムが読取光48により
照射されながら矢印Y′方向に搬送される間、X線フィ
ルム40を透過した読取光をX線画像のY′方向の各画
素間隔に対応した所定の時間間隔で受光する。
第7図は、上記MOSセンサ45の等価回路を示した回
路図である。
多数の固体光電変換素子46に読取光43が当たって発
生するフォトキャリアによる信号は、固体光電変換素子
46内のキャパシタCi  (i =1,2.・・・・
・・n)に蓄積される。蓄積されたフォトキャリアの信
号は、シフトレジスタ47によって制御されるスイッチ
部48の順次開閉により順次読み出され、これにより時
系列化された画像信号が得られる。この画像信号は、そ
の後増幅器49で増幅されてその出力端子50から出力
される。
出力されたアナログの画像信号はサンプリングされてデ
ィジタルの画像信号に変換され、その後、この画像信号
に基づいて、前述した実施例と同様にして、被写体像内
の画像点が決定され、この画像点を用いて照射野が認識
される。尚、本実施例において、MOSセンサ45の代
わりにCCD、CP D (Charge T’rim
ing Device)等を用いることができることは
いうまでもない。またX線フィルムの読取りにおいても
、前述した蓄積性蛍光体シートの読取りと同様に光ビー
ムで2次元的に走査して読取りを行なってもよいことも
もちろんである。
薫た上記実施例ではX線フィルム40を透過した光を受
光しているが、X線フィルム40から反射した光を受光
するように構成することができることももちろんである
このように、本発明の被写体像内画像点決定方法は、被
写体の放射線画像が記録された記録シートから得られた
放射線画像を表わす光を読み取って画像信号を得、この
画像信号に基づいて放射線画像を再生出力する放射線画
像読取再生装置一般に適用することができる。
(発明の効果) 以上詳細に説明したように、本発明の被写体像内画像点
決定方法のひとつは、各画素に対応する画像信号値また
はこの逆数でそれぞれ対応する各画素を重み付けして記
録シートの重心を求めることにより、この重心を被写体
像内の画像点として定めることができる。また、本発明
の他の被写体像内画像点決定方法は、各画素に対応する
画像信号値またはこの画像信号値の逆数をそれぞれ対応
する各画素に対応づけられたときに、記録シート上の互
いに異なる2つの方向のそれぞれについて上記画像信号
値または上記逆数を上記各方向に累積してプロットした
累積分布を求め、これらの累積分布のそれぞれについて
最大の累積値の略半分の値に対応する各方向についての
座標点を求めるものであるため、これらの座標点により
定まる記録シート上の位置が、上記重心と同様に、被写
体像内の画像点として求められる。
また、上記各方法を、ネガ像又はポジ像の一方のみを取
り扱うシステム、または記録シートとして蓄積性蛍光体
シートを用いるシステム等、ネガ。
ポジの区別のないシステムなどに適用する場合において
は、記録シートの周辺部に対応する画像信号値を代表す
る第1の代表値と、記録シートの全部または略中央部に
対応する画像信号値を代表する第2の代表値とを求め、
これら第1の代表値と第2の代表値との大小を比較する
ことにより、照射野絞りが使用され記録シートのうちの
狭い範囲にのみ放射線画像が記録されているか、あるい
は照射野絞りを使用せずにまたは使用しても記録シート
の広い範囲にわたって放射線画像が記録されているかを
判定することができ、この比較結果に応じて上記画像信
号値または上記逆数の一方を選択して、この選択された
上記画像信号値または逆数を用いて上記いずれかの方法
により被写体像内の画像点を定めることができる。
【図面の簡単な説明】
第1図は、本発明の被写体像内画像点決定方法の一例を
使用した、放射線画像読取再生装置の一実施例の斜視図
、 第2A図は、蓄積性蛍光体シート上の狭い領域にのみ放
射線か照射された様子を表わした図、第2B図は、蓄積
性蛍光体シート全面にわたって放射線が照射された様子
を表わした図、第2C図は、蓄積性蛍光体シートの全面
に対応する画像信号S0のヒストグラム(グラフA)と
累積ヒストグラム(グラフB)の−例を示したグラフ、 第3A図、第3B図は、画像信号値S Q +逆数1 
/ S oのグラフの一例を、放射線画像とともに示す
図、 第4図は、X方向、X方向の画像信号の値の分布および
それらの分布の累積分布の一例を表わした図、 第5A図は、記録シートの中央に照射野のある放射線画
像を表わした図、 第5B図は、照射野が記録シートの中心から外れた放射
線画像を表わした図、 第5C図は、照射野絞りを使用せずに撮影したことによ
り記録シート全面が照射野である放射線画像を表わした
図、 第6図は、X線フィルムに記録されたX線画像を読み取
るX線画像読取装置の一実施例の斜視図、第7図は、M
OSセンサの等価回路を示した回路図である。 1・・・記録シート2.12・・・照射野3.13・・
・被写体像  11・・・蓄積性蛍光体シート14・・
・周辺部     12・・・輝尽発光光23・・・フ
ォトマルチプライヤ 26・・・ログアンプ   27・・・A/D変換器2
8・・・記憶手段    29・・・画像処理手段30
・・・再生手段    40・・・X線フィルム45・
・・MOSセンサ 第 3A 図 第 B 図 区 Uノ 沫

Claims (2)

    【特許請求の範囲】
  1. (1)被写体像を含む放射線画像が記録された記録シー
    ト上の各画素から得られた画像信号に基づいて、前記各
    画素に対応する画像信号値または画像信号値の逆数でそ
    れぞれ対応する前記各画素を重み付けして前記記録シー
    トの重心を求め、この重心を、前記被写体像内の画像点
    とすることを特徴とする被写体像内画像点決定方法。
  2. (2)被写体像を含む放射線画像が記録された記録シー
    ト上の各画素から得られた画像信号に基づいて、前記各
    画素に対応する画像信号値または画像信号値の逆数をそ
    れぞれ対応する前記各画素に対応づけたときに、前記記
    録シート上の互いに異なる2つの方向のそれぞれについ
    て、画像信号値または画像信号値の逆数を各前記方向に
    累積してプロットした累積分布を求め、これらの累積分
    布のそれぞれについて最大の累積値の略半分の値に対応
    する各前記方向についての座標点を求め、これらの座標
    点により定まる前記記録シート上の位置を、前記被写体
    像内の画像点とすることを特徴とする被写体像内画像点
    決定方法。(3)前記画像信号を得た後、この画像信号
    に基づいて、前記記録シートの周辺部に対応する前記画
    像信号値を代表する第1の代表値と、前記記録シートの
    全部または略中央部に対応する前記画像信号値を代表す
    る第2の代表値とを求め、これら第1の代表値と第2の
    代表値との大小を比較し、この比較結果に応じて画像信
    号値または画像信号値の逆数の一方を選択し、このよう
    にして選択された前記画像信号値または画像信号値の逆
    数を用いて前記被写体像内の画像点を定めることを特徴
    とする請求項1または2記載の被写体像内画像点決定方
    法。
JP1029179A 1988-04-20 1989-02-08 被写体像内画像点決定方法 Expired - Lifetime JP2631737B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP1029179A JP2631737B2 (ja) 1988-04-20 1989-02-08 被写体像内画像点決定方法
EP89107156A EP0340553B1 (en) 1988-04-20 1989-04-20 Method for determining an image point in an object image
DE68921786T DE68921786T2 (de) 1988-04-20 1989-04-20 Verfahren zur Ermittlung eines Bildpunktes eines Objekts.
US07/340,744 US5179597A (en) 1988-04-20 1989-04-20 Method for determining an image point in an object image

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP9789788 1988-04-20
JP63-97897 1988-04-20
JP1029179A JP2631737B2 (ja) 1988-04-20 1989-02-08 被写体像内画像点決定方法

Publications (2)

Publication Number Publication Date
JPH0228782A true JPH0228782A (ja) 1990-01-30
JP2631737B2 JP2631737B2 (ja) 1997-07-16

Family

ID=26367336

Family Applications (1)

Application Number Title Priority Date Filing Date
JP1029179A Expired - Lifetime JP2631737B2 (ja) 1988-04-20 1989-02-08 被写体像内画像点決定方法

Country Status (4)

Country Link
US (1) US5179597A (ja)
EP (1) EP0340553B1 (ja)
JP (1) JP2631737B2 (ja)
DE (1) DE68921786T2 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05276440A (ja) * 1991-04-30 1993-10-22 Fuji Photo Film Co Ltd 被写体像内画像点決定方法

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5506917A (en) * 1990-07-13 1996-04-09 Nippon Telegraph And Telephone Corporation Thresholding method for segmenting gray scale image, method for determining background concentration distribution, and image displacement detection method
JP2694580B2 (ja) * 1991-03-07 1997-12-24 富士写真フイルム株式会社 被写体像内画像点決定方法
DE4129656C2 (de) * 1991-09-06 1994-02-10 Siemens Ag Wiedergabevorrichtung für Videosignale auf einem Monitor
US5272760A (en) * 1992-05-29 1993-12-21 Cimpiter Corporation Radiographic image evaluation apparatus and method
US5384862A (en) * 1992-05-29 1995-01-24 Cimpiter Corporation Radiographic image evaluation apparatus and method
EP0599345B1 (en) * 1992-11-27 2002-06-05 Fuji Photo Film Co., Ltd. Method for adjusting positions of radiation images
US5566245A (en) * 1993-03-09 1996-10-15 United Parcel Service Of America, Inc. The performance of a printer or an imaging system using transform-based quality measures
US5933524A (en) * 1993-09-27 1999-08-03 Siemens Aktiengesellschaft Method for segmentation of digital color images
US5588071A (en) * 1994-10-26 1996-12-24 Minnesota Mining And Manufacturing Company Identifying an area of interest using histogram data arranged in predetermined sequence

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61180376A (ja) * 1985-02-05 1986-08-13 Mitsubishi Heavy Ind Ltd 重心点演算装置
JPS62254277A (ja) * 1986-04-28 1987-11-06 Toshiba Corp 画像計測装置
JPS63182503A (ja) * 1987-01-26 1988-07-27 Matsushita Electric Ind Co Ltd 立体物測定装置

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1980001002A1 (en) * 1978-10-30 1980-05-15 Fujitsu Ltd Pattern inspection system
US4516263A (en) * 1982-04-30 1985-05-07 The Charles Stark Draper Laboratory, Inc. Spatially integral, video signal processor
US4701807A (en) * 1983-09-22 1987-10-20 Canon Kabushiki Kaisha Method and apparatus for processing an image
US4593406A (en) * 1984-01-16 1986-06-03 The United States Of America As Represented By The United States Department Of Energy Automatic image acquisition processor and method
DE3588169T2 (de) * 1984-07-31 1998-02-19 Fuji Photo Film Co Ltd Verfahren zur Detektion einer Bildbelichtungsfläche in einem Bildauslesevorgang
JPS6276358A (ja) * 1985-09-27 1987-04-08 Fuji Photo Film Co Ltd 放射線画像情報の読取条件決定方法
JPH0671300B2 (ja) * 1987-03-20 1994-09-07 富士写真フイルム株式会社 放射線画像情報の読取処理条件決定装置
JPH0721621B2 (ja) * 1987-04-16 1995-03-08 富士写真フイルム株式会社 放射線画像の分割パタ−ン認識方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61180376A (ja) * 1985-02-05 1986-08-13 Mitsubishi Heavy Ind Ltd 重心点演算装置
JPS62254277A (ja) * 1986-04-28 1987-11-06 Toshiba Corp 画像計測装置
JPS63182503A (ja) * 1987-01-26 1988-07-27 Matsushita Electric Ind Co Ltd 立体物測定装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05276440A (ja) * 1991-04-30 1993-10-22 Fuji Photo Film Co Ltd 被写体像内画像点決定方法

Also Published As

Publication number Publication date
DE68921786T2 (de) 1995-07-13
EP0340553B1 (en) 1995-03-22
EP0340553A1 (en) 1989-11-08
US5179597A (en) 1993-01-12
DE68921786D1 (de) 1995-04-27
JP2631737B2 (ja) 1997-07-16

Similar Documents

Publication Publication Date Title
US5067163A (en) Method for determining a desired image signal range from an image having a single background
JPH083843B2 (ja) 放射線画像読取装置
JP2631737B2 (ja) 被写体像内画像点決定方法
US4994662A (en) Radiation image read-out apparatus and method for operating the same
JPH0296883A (ja) 放射線照射野認識方法
JP2582640B2 (ja) 放射線画像の分割パターン認識方法
US4963739A (en) Radiation image read-out apparatus
US5091970A (en) Method for judging the presence or absence of limited irradiation field
JP2631742B2 (ja) 被写体像内画像点決定方法
JP2707355B2 (ja) 画像記録位置認識装置
JPH07104944B2 (ja) 放射線画像分割パターン認識正誤判定方法
JP2704288B2 (ja) 被写体像内画像点決定方法
JPH02111346A (ja) 乳房撮影画像における被写体部認識方法
JPH0214376A (ja) 放射線画像読取装置
JP2631741B2 (ja) 放射線画像の分割パターン認識方法
JP2524208B2 (ja) 未露光認識方法
JP3165530B2 (ja) 放射線画像読取条件決定方法および装置
JPH0271247A (ja) 照射野認識方法
JPH0271249A (ja) 放射線照射野判定方法
JPH0453540A (ja) 胸壁入り乳房放射線画像における読取条件及び/又は画像処理条件決定方法
JPH0267690A (ja) 放射線照射野判定方法
JPH02110787A (ja) 放射線画像読取条件及び/又は画像処理条件決定方法
JPH05181953A (ja) 被写体像内画像点決定方法
JPH0752284B2 (ja) 未露光認識方法
JPH02271343A (ja) 放射線画像の分割パターン認識方法

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080425

Year of fee payment: 11

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090425

Year of fee payment: 12

EXPY Cancellation because of completion of term