JPH021080A - 放射線照射野輪郭候補点正誤判定方法 - Google Patents

放射線照射野輪郭候補点正誤判定方法

Info

Publication number
JPH021080A
JPH021080A JP63186084A JP18608488A JPH021080A JP H021080 A JPH021080 A JP H021080A JP 63186084 A JP63186084 A JP 63186084A JP 18608488 A JP18608488 A JP 18608488A JP H021080 A JPH021080 A JP H021080A
Authority
JP
Japan
Prior art keywords
point
contour
contour candidate
radiation
reading
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
JP63186084A
Other languages
English (en)
Other versions
JP2582625B2 (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 JP63186084A priority Critical patent/JP2582625B2/ja
Priority to US07/324,879 priority patent/US4992663A/en
Priority to DE8989104828T priority patent/DE68902255T2/de
Priority to EP89104828A priority patent/EP0335204B1/en
Publication of JPH021080A publication Critical patent/JPH021080A/ja
Application granted granted Critical
Publication of JP2582625B2 publication Critical patent/JP2582625B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

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

Description

【発明の詳細な説明】 (産業上の利用分野) 本発明は、放射線画像情報が記録された蓄積性蛍光体シ
ートなどの記録媒体から上記放射線画像情報を読み取る
際に、該記録媒体における放射線照射野を認識するため
に、この照射野の輪郭部分上に位置するとして検出した
輪郭候補点が正しいか否かを判定する方法に関するもの
である。
(従来の技術) 記録された放射線画像を読み取って画像信号を得、この
画像信号に適切な画像処理を施した後、画像を再生記録
することは種々の分野で行なわれている。たとえば、後
の画像処理に適合するように設計されたガンマ値の低い
X線フィルムを用いてX線画像を記録し、このX線画像
が記録されたフィルムからX線画像を読み取って電気信
号に変換し、この電気信号(画像信号)に画像処理を施
した後コピー写真等に可視像として再生することにより
、コントラスト シャープネス、粒状性等の画質性能の
良好な再生画像を得ることのできるシステムか開発され
ている(特公昭61−5193号公報参照)。
また本願出願人により、放射線(X線、α線。
β線、γ線、電子線、紫外線等)を照射するとこの放射
線エネルギーの一部が蓄積され、その後可視光等の励起
光を照射すると蓄積されたエネルギーに応じて輝尽発光
を示す蓄積性蛍光体(輝尽性蛍光体)を利用して、人体
等の被写体の放射線画像を一部シート状の蓄積性蛍光体
に撮影記録し、この蓄積性蛍光体シートをレーザー光等
の励起光で走査して輝尽発光光を生ぜしめ、得られた輝
尽発光光を光電的に読み取って画像信号を得、この画像
信号に基づき被写体の放射線画像を写真感光材料等の記
録材料、CRT等に可視像として出力させる放射線画像
記録再生システムがすでに提案されている(特開昭55
−12429号、同5B−11395号。
同55−163472号、同513−104645号、
同55−116340号等)。
このシステムは、従来の銀塩写真を用いる放射線写真シ
ステムと比較して極めて広い放射線露出域にわたって画
像を記録しうるという実用的な利点を有している。すな
わち、蓄積性蛍光体においては、放射線露光量に対して
蓄積後に励起によって輝尽発光する発光光の光量が極め
て広い範囲にわたって比例することが認められており、
従って種々の撮影条件により放射線露光量がかなり大幅
に変動しても、蓄積性蛍光体シートより放射される輝尽
発光光の光量を読取ゲインを適当な値に設定して充電変
換手段により読み取って電気信号に変換し、この電気信
号を用いて写真感光材料等の記録材料、CRT等の表示
装置に放射線画像を可視像として出力させることによっ
て、放射線露光量の変動に影響されない放射線画像を得
ることができる。
上記システムにおいて、蓄積性蛍光体シートに照射され
た放射線の線量等に応じて最適な読取条件で読み取って
画像信号を得る前に、予め低レベルの光ビームにより蓄
積性蛍光体シートを走査してこのシートに記録された放
射線画像の概略を読み取る先読みを行ない、この先読み
により得られた先読画像信号を分析し、その後上記シー
トに上記先読みの際の光ビームよりも高レベルの光ビー
ムを照射して走査し、この放射線画像に最適な読取条件
で読み取って画像信号を得る本読みを行なうように構成
されたシステムもある(特開昭58−67240号、同
58−67241号、同58−67242号等)。
ここで読取条件とは、読取りにおける輝尽発光光の光量
と読取装置の出力との関係に影響を与える各種の条件を
総称するものであり、例えば人出力の関係を定める読取
ゲイン、スケールファクタあるいは、読取りにおける励
起光のパワー等を意味するものである。
また、光ビームの高レベル/低レベルとは、それぞれ、
上記シートの単位面積当りに照射される光ビームの強度
の大/小、もしくは上記シートから発せられる輝尽発光
光の強度が上記光ビームの波長に依存する(波長感度分
布を有する)場合は、上記シートの単位面積当りに照射
される光ビームの強度を上記波長感度で重みづけした後
の重みづけ強度の大〆小をいい、光ビームのレベルを変
える方法としては、異なる波長の光ビームを用いる方法
、レーザ光源等から発せられる光ビームの強度そのもの
を変える方法、光ビームの光路上にNDフィルター等を
挿入、除去することにより光ビームの強度を変える方法
、光ビームのビーム径を変えて走査密度を変える方法、
走査速度を変える方法等、公知の種々の方法を用いるこ
とができる。
また、この先読みを行なうシステムか先読みを行なわな
いシステムかによらず、得られた画像信号(先読画像信
号を含む)を分析し、画像信号に画像処理を施す際の最
適な画像処理条件を決定するようにしたシステムもある
。この画像信号に基づいて最適な画像処理条件を決定す
る方法は、蓄積性蛍光体シートを用いるシステムに限ら
れず、たとえば従来のX線フィルム等の記録媒体に記録
された放射線画像から画像信号を得るシステムにも適用
されている。
上記画像信号(先読画像信号を含む)を分析して最適な
読取条件、画像処理条件を求める方法は種々提案されて
いるが、その方法のひとつとして、画像信号のヒストグ
ラムを作成する方法が知られている(たとえば、特開昭
60−1513055号)。画像信号のヒストグラムを
求めることにより、たとえば画像信号の最大値、最小値
や、頻度が最大となる点の画像信号の値等を知ることが
でき、これらの6値から蓄積性蛍光体シート、X線フィ
ルム等の記録媒体に記録された放射線画像の特徴を把握
することができる。そこでこのヒストグラムに基づいて
最適な読取条件1画像処理条件を求めることにより、観
察適正のすぐれた放射線画像を再生出力することが可能
となる。
一方、記録媒体に放射線画像を撮影記録するに際しては
、被写体の観察に必要の無い部分に放射線を照射しない
ようにするため、あるいは観察に不要な部分に放射線を
照射するとその部分から観察に必要な部分に散乱線が入
り画質性能が低下するため、放射線が被写体の必要な部
分および記録媒体の一部にのみ照射されるように放射線
の照射域を制限する照射野絞りを使用して撮影を行なう
ことも多い。
ところが、前述のようにして画像信号を分析して読取条
件1画像処理条件を求めるにあたって、分析に用いた画
像信号が、照射野絞りを用いて撮影した記録媒体から得
られた画像信号である場合、この照射野の存在を無視し
て画像信号を分析しても撮影記録された放射線画像が正
しく把握されず、誤った読取条件、画像処理条件が求め
られ観察適正の優れた放射線画像が再生記録されない場
合が生ずる。
これを解決するためには、読取条件1画像処理条件を求
める前に、照射野を認識し、照射野内の画像信号に基づ
いて読取条件1画像処理条件を求める必要がある。
本出願人は既に、放射線照射野を認識する方法をいくつ
か提案しており(例えば特開昭61−39039号)、
このような方法によって照射野を認識し、その認識領域
のみに対応する画像信号に基づいて読取条件、画像処理
条件を求めるようにすれば、上述の不具合は解消可能で
ある。
上述のような放射線照射野を認識する方法においては多
くの場合、まず照射野の輪郭上にあると考えられる点、
すなわち輪郭候補点をいくつか求めるようにしている。
そしてこのような輪郭候補点がいくつか求まったならば
、次にそれらの点に沿う直線あるいは曲線を求めれば、
これらの直線あるいは曲線の内側を放射線照射野と認識
することができる。
上記の輪郭候補点を検出する方法としては例えば、記録
媒体に記録された放射線画像を読み取って該記録媒体上
の各位置に対応するデジタル画像データを求め、次にこ
の画像データを記録媒体上の1本のラインに沿って微分
処理し、この処理によって得られた微分値の絶対値が所
定のしきい値を超える点を、あるいはそのような点がい
くつか存在したならばそれらのうち最も記録媒体の端部
に近い点を輪郭候補点として検出する、という方法が既
に提案されている(例えば特開昭62−15538号参
照)。また上記輪郭候補点はその他、パターンマツチン
グによる手法や、直線あてはめを行なってその傾きから
輪郭部を判別する方法等によって検出することもできる
(発明が解決しようとする課題) ところで、上述のような種々の方法によって輪郭候補点
を検出する際、画像内に骨の辺縁部等照射野輪郭部と同
様に濃度が急激に変化する部分が存在したり、あるいは
照射野外に散乱放射線のエネルギーが蓄積されていたり
すると、それらの部分の点が照射野輪郭候補点として誤
検出されることがある。このような誤検出を完全に無く
すことは困難であるが、ある輪郭候補点が誤検出である
旨が正確に分かれば、その輪郭候補点をキャンセルした
り、あるいは他の方法で検出した輪郭候補点に訂正する
ことにより、照射野を誤認識することは避けられる。
そこで本発明は、前述のようにして検出した輪郭候補点
が正しいか、あるいは誤りであるかを正確に判定するこ
とができる方法を提供することを目的とするものである
(課題を解決するための手段及び作用)本発明による放
射線照射野輪郭候補点正誤判定方法は、前述のようにし
て検出した輪郭候補点が所定の判定基準を満たしている
か否かを調べ、満たしていない場合はその輪郭候補点を
照射野輪郭上には存在しない誤検出の点であると判定す
ることを特徴とするものである。
また本発明のもう一つの放射線照射野輪郭候補点正誤判
定方法は、上述のような判定基準を調べるとともにさら
に、各輪郭候補点から、それ以外の輪郭候補点に基づい
て照射野輪郭上にあると想定される点までの平均画像濃
度を求め、輪郭候補点が上記判定基準を満たさず、かつ
上記平均画像濃度が所定のしきい値を下回る場合は、こ
の輪郭候補点が照射野輪郭上に位置しない誤検出の点で
あると判定することを特徴とするものである。
上記の判定基準としてより具体的には、例えば正誤判定
する輪郭候補点から画像内の所定点までの距離と、それ
以外の輪郭候補点から上記所定点までの距離との関係や
、正誤判定する輪郭候補点からそれ以外の輪郭候補点ま
での距離、さらには正誤判定する輪郭候補点からそれ以
外の輪郭候補点2つに向けてそれぞれ延ばした2本の直
線がなす角度等を用いることができる。
(実 施 例) 以下、図面に示す実施例に基づいて本発明の詳細な説明
する。
第1図は本発明の方法によって放射線照射野輪郭候補点
の正誤を判定するようにした放射線画像情報記録再生シ
ステムを示すものである。この放射線画像情報記録再生
システムは基本的に、放射線画像撮影部20、先読み用
読取部30、本読み用読取部40、および画像再生部5
0から構成されている。
放射線画像撮影部20においては、例えばX線管球等の
放射線源100から被写体(被検者)101に向けて、
放射線102が照射される。この被写体101を透過し
た放射線102が照射される位置には、先に述べたよう
に放射線エネルギーを蓄積する蓄積性蛍光体シート10
3が配置され、この蓄積性蛍光体シー) 103に被写
体lotの透過放射線画像情報が蓄積記録される。なお
放射線源100と被写体101との間には、放射線10
2の照射野を絞る絞り104が配されている。
このようにして被写体101の放射線画像情報が記録さ
れた蓄積性蛍光体シート103は、移送ローラ等のシー
ト移送手段110により、先読み用読取部30に送られ
る。先読み用読取部30において先読み用レーザ光源2
01から発せられたレーザ光202は、このレーザ光2
02の励起によって蓄積性蛍光体シート103から発せ
られる輝尽発光光の波長領域をカットするフィルター2
03を通過した後、ガルバノメータミラー等の光偏向器
204により直線的に偏向され、平面反射鏡205を介
して蓄積性蛍光体シート103上に入射する。ここでレ
ーザ光源201は、励起光としてのレーザ光202の波
長域が、蓄積性蛍光体シート103が発する輝尽発光光
の波長域と重複しないように選択されている。他方、蛍
光体シート103は移送ローラ等のシート移送手段21
0により矢印206の方向に移送されて副走査がなされ
、その結果、蛍光体シート103の全面にわたってレー
ザ光202が照射される。ここで、レーザ光源201の
発光強度、レーザ光202のビーム径、レーザ光202
の走査速度、蓄積性蛍光体シート103の移送速度は、
先読みの励起光(レーザ光202)のエネルギーが、後
述する本読み用読取部40で行なわれる本読みのそれよ
りも小さくなるように選択されている。
上述のようにレーザ光202が照射されると、蓄積性蛍
光体シート103は、それに蓄積記録されている放射線
エネルギーに対応した光量の輝尽発光光を発し、この発
光光は先読み用光ガイド207に入射する。輝尽発光光
はこの光ガイド207内を導かれ、射出面から射出して
フォトマルチプライヤ−等の光検出器208によって受
光される。該光検出器208の受光面には、輝尽発光光
の波長域の光のみを透過し、励起光の波長域の光を力・
ソトするフィルターが貼着されており、輝尽発光光のみ
を検出し得るようになっている。検出された輝尽発光光
は蓄積記録情報を担持する電気信号に変換され、増幅器
209により増幅される。増幅器209から出力された
信号はA/D変換器211によりデジタル化され、先読
み画像信号Spとして本読み用読取部40の本読み制御
回路314に入力される。この本読み制御回路314は
、先読み画像信号spが示す蓄積記録情報に基づいて、
読取ゲイン設定値a1収録スケールファクター設定値b
1再生画像処理条件設定値Cを決定する。また上記先読
み画像信号spは、後に詳述する照射野認識回路220
にも人力される。
以上のようにして先読みを完了した蓄積性蛍光体シート
103は本読み用読取部40へ移送される。
本読み用読取部40において本読み用レーザ光源301
から発せられたレーザ光302は、このレーザ光302
の励起によって蓄積性蛍光体シート103から発せられ
る輝尽発光光の波長領域をカットするフィルター303
を通過した後、ビームエクスパンダ−304によりビー
ム径の大きさが厳密に調整され、ガルバノメータミラー
等の光偏向器305によって直線的に偏向され、平面反
射鏡306を介して蓄積性蛍光体シート103上に入射
する。光偏向器305と平面反射鏡306との間にはf
θレンズ307が配され、蓄積性蛍光体シート103上
を走査するレザ光302のビーム径が均一となるように
されている。他方、蓄積性蛍光体シート103は移送ロ
ーラなどのシート移送手段320により矢印308の方
向に移送されて副走査がなされ、その結果、蓄積性蛍光
体シート103の全面にわたってレーザ光が照射される
。このようにレーザ光302が照射されると、蓄積性蛍
光体シート103はそれに蓄積記録されている放射線エ
ネルギーに対応した光量の輝尽発光光を発し、この発光
光は本読み用光ガイド309に入射する。本読み用光ガ
イド309の中を全反射を繰返しつつ導かれた輝尽発光
光はその射出面から射出され、フォトマルチプライヤ−
等の光検出器310によって受光される。光検出器31
0の受光面には、輝尽発光光の波長域のみを選択的に透
過するフィルターが貼着され、光検出器310が輝尽発
光光のみを検出するようになっている。
蓄積性蛍光体シート103に記録されている放射線画像
を示す輝尽発光光を光電的に検出した光検出器310の
出力は、前記制御回路314が決定した読取ゲイン設定
値aに基づいて読取ゲインが設定された増幅器311に
より、適正レベルの電気信号に増幅される。増幅された
電気信号はA/D変換器312に入力され、収録スケー
ルファクター設定値すに基づいて、信号変動幅に適した
収録スケールファクターでデジタル信号に変換されて信
号処理回路313に入力される。上記デジタル信号は、
この信号処理回路313において、観察読影適性の優れ
た放射線画像が得られるように再生画像処理条件設定値
Cに基づいて信号処理(画像処理)され、出力される。
信号処理回路313から出力された読取画像信号(本読
み画像信号)Soは、画像再生部50の光変調器401
に入力される。この画像再生部50においては、記録用
レーザ光源402からのレーザ光403が光変調器40
1により、上記信号処理回路313から人力される本読
み画像信号Soに基づいて変調され、走査ミラー404
によって偏向されて写真フィルム等の感光材料405上
を走査する。そして感光材料405は上記走査の方向と
直交する方向(矢印406方向)に走査と同期して移送
され、感光材料405上に、上記本読み画像信号Soに
基づく放射線画像が出力される。放射線画像を再生する
方法としては、このような方法の他、前述したCRTに
よる表示等、種々の方法を採用することができる。
次に、前記第2図に示されるように蓄積性蛍光体シート
103において放射線照射野Bが絞られている場合にも
、前記読取ゲイン設定値a1収録スケールファクター設
定値b1画像処理条件設定値Cが適正に決定される仕組
みについて、第5図を参照して説明する。この第5図に
示されるように前記制御回路314は、信号抽出部35
0、ヒストグラム解析部3511読出部352および記
憶部353からなる。先読み画像信号Spは上記信号抽
出部350に入力され、該信号抽出部350において、
後述するようにして指定される領域のみについての先読
み画像信号Sp’が抽出される。この信号抽出部350
から出力される先読み画像信号Sp’ はヒストグラム
解析部351に入力される。ヒストグラム解析部351
は先読み画像信号Sp’ のヒストグラムを作成し、例
えばその最大値、最小値、最大頻度値等を求め、それら
の値を示す情報Srを読出部352に送る。記憶部35
3にはこれら最大値、最小値等に対応する最適の読取ゲ
イン設定値a1収録スケールファクター設定値すおよび
画像処理条件設定値Cが記憶されており、読出部352
は上記情報Srに対応する設定値a、bScを記憶部3
53から読み出して、前述のようにそれぞれ増幅器31
1 、A/D変換器312および信号処理回路313に
送る。
次に信号抽出部350における信号抽出について説明す
る。照射野認識回路220は微分処理部221、しきい
値設定部222、輪郭候補点信号検出部223、正誤判
定部224および演算部225からなる。先読み画像信
号Spはこの照射野認識回路220において、微分処理
部221と輪郭候補点信号検出部223とに入力される
。微分処理部221はデジタル化されているこの先読み
画像信号Spを、まず第3図に示す方向D1に沿って微
分処理し、以下同様に方向D2、D3・・・・・・Dn
に沿って微分処理する。
この微分の方法は、1次元の1次微分でも高次の微分で
もよいし、また2次元の1次微分や高次の微分でもよい
。また、離散的に標本化された画像の場合、微分すると
は近傍に存在する画像データ同志の差分を求めることと
等価であり、本例ではこの差分を求める。上記複数の方
向D1〜Dnは、蓄積性蛍光体シート103の中心0か
らシート端部に向かう放射状の方向であり、本例では、
各方向D1〜Dnが互いに等角度間隔に設定されている
またこのような放射状の方向D1〜Dnは、例えば蓄積
性蛍光体シート103のサイズが半切サイズ(256x
L92 mm)の場合、64方向程度設定される。
このような微分処理を行なうことにより、上記の差分が
求められる。この差分を示す情報5Illは、輪郭候補
点信号検出部223に送られる。輪郭候補点信号検出部
223は上記差分を示す情報5Illと、しきい値設定
部222が出力するしきい値Thを示す情報sthとか
ら、シート103上の放射線照射野Bの輪郭部分にある
と考えられる輪郭候補点を求める。すなわち、照射野B
内についての画像信号のレベルは、照射野B外の領域に
ついての画像信号のレベルに比べて全体的に明らかに高
い値をとるので、ある方向D1に沿った先読み画像信号
Spの値は、第4図(a)に示すような分布をとる。
したがって上記差分の値は第4図(b)に示すように、
照射野エツジ部分において特異的に大きく変化する。そ
こで輪郭候補点信号検出部223は、この差分が負とな
りその絶対値が前記所定のしきい値Thを超える点を検
出して、輪郭候補点を求める。
輪郭候補点信号検出部223は、上述のようにして求め
た輪郭候補点についての先読み画像信号Spを抽出し、
その抽出された各先読み画像信号Spに対応する画素位
置を求め、その画素位置を示す情報Seを正誤判定部2
24に送る。上述のようにして抽出された先読み画像信
号Spは、大部分が蓄積性蛍光体シート103上の放射
線照射野B(第2図参照)のエツジ部分を担う画像信号
、つまり輪郭候補点信号となる。本例において、上記画
素位置は第2図に示すように、蓄積性蛍光体シト103
上のx−y直交座標系で表わされる。
正誤判定部224は、上記情報Seが示す画素位置の輪
郭候補点がそれぞれ、本当に照射野輪郭部に位置するも
のであるか否かを、以下に述べるようにして判定する。
正誤判定部224は、第13図に示すように、正誤を判
定する輪郭候補点をEmlその両端の輪郭候補点をそれ
ぞれErTl+I、Em−+としたとき、蓄積性蛍光体
シート103の中心Oから輪郭候補点Em、ETTl+
1、およびE m−1までの各距離免、Q、lおよび!
2.2について、 であるか否かを調べる(なおnは、前述した微分処理の
方向の数である)。通常そうであるように、本例におい
ても照射野Bの形状は、星形のように内方に入り込む部
分は持たない多角形であり、またシート中心Oを必ずこ
の照射野B内に位置させるということを前提として照射
野絞りがなされるようになっている。上記(1)式が満
足されないということは、第13図に示すように輪郭候
補点Emがその両隣の輪郭候補点EIIl+1およびE
 m−1よりもシート中心側に入り込んでいることを示
すが、照射野Bか上記の通りの多角形である限り、照射
野輪郭上にこのような入り込んだ点Emは存在しないは
ずである。したがって正誤判定部224は、上記(1)
式か満たされない場合は、輪郭候補点Emは照射野内の
点が誤検出されたものであるとみなして、前記画素位置
情報Seの中からこの輪郭候補点En1についての情報
を削除する。一方上記(1)式が満たされれば、輪郭候
補点Emについての画素位置情報はそのまま生かされる
。正誤判定部224は、検出されたすべての輪郭候補点
について上記の判定処理、および必要であれば上記削除
処理を行ない、この処理済みの画素位置情報Se’ を
演算部225に送る。
なおこの例では、輪郭候補点Enlが誤検出点であると
判定したとき、この点をキャンセルしているが、そのよ
うな点を、他の輪郭候補点から導かれるより適正な点に
訂正するようにしてもよい。
上記の例にあっては、そのような適正な点として例えば
、シート中心Oから誤検出の輪郭候補点Emと同方向に
あって、シート中心0からだけ離れた点等を用いること
ができる。
以上述べたようにして正しい輪郭候補点を幾つか求めた
後、これらの点に沿った線を求めれば、その線が照射野
の輪郭となることは、先に述べた通りである。この輪郭
候補点に沿った線は、例えばそれらの点を平滑化処理し
た後残った点を連結する方法、局所的に最小二乗法を適
用して複数の直線を求め、それらを連結する方法、スプ
ライン曲線等を当てはめる方法等によって求めることが
できるが、本例において演算部225は、特にHoug
h変換を利用して輪郭候補点に沿った複数の直線を求め
るように構成されている。以下、この直線を求める処理
について詳しく説明する。
演算部225は上記情報Se’が示す画素位置(輪郭候
補点)の座標を(Xo、Yo)としたとき、これらのx
oSyoを定数として p −X Q  CO8θ+yOsinθで表わされる
曲線を、すべての輪郭候補点座標(Xo、Yo)につい
て求める。この曲線は第6図に示すようなものとなり、
輪郭候補点座標(X。、Yo)の数だけ存在する。
次いで演算部225は、上記複数の曲線のうちの所定数
0以上の曲線が互いに交わる交点(ρ0゜θ0)を求め
る。なお、輪郭候補点座標(xo。
yo)の誤差等のため、多数の曲線が厳密に一点で交わ
ることは少ないので、実際には例えば2本の曲線の交点
が互いに微小所定値以下の間隔で存在するとき、それら
の交点群の中心を上記交点(ρ0.θ0)とする。次に
演算部225は、交点(ρ0.θ0)から前記x−y直
交座標系において次式 %式% で規定される直線を求める。この直線は、複数の輪郭候
補点座標(Xo、Vo)に沿って延びる直線となる。な
お放射線照射野B内において急激に濃度が変化する骨辺
縁部等も、上記輪郭候補点として検出されることがある
。したがって第2図にも示すように、このような輪郭候
補点と照射野輪郭部の輪郭候補点とを結ぶ直線りが求め
られる可能性があるが、前述の所定数Qを十分に大きく
(例えば20本以上等)設定しておけば、上記のような
直線りは求められない。つまり多数の輪郭候補点に沿う
、照射野輪郭を示す直線のみが求められる。
上述した直線は、輪郭候補点が第2図図示のように分布
している場合、第7図図示のようなものとなる。演算部
225は次に、こうして求めた複数の直線L1、L2、
L3・・・Lnによって囲まれる領域を求め、この領域
を放射線照射野Bとして認識する。この領域は、詳しく
は例えば以下に述べるようにして認識される。演算部2
25は蓄積性蛍光体シート103の隅部と中心Gとを結
ぶ線分M1、M2 、M3・・・M+a(蓄積性蛍光体
シート103が矩形の場合は4本)を記憶しており、こ
れら各線分M1〜MIIlと上記各直線t、、−L、と
の交点の有無を調べる。この交点が存在した場合、演算
部225は上記直線によって2分される平面のうち、シ
ート隅部を含む側の平面を切り捨てる。この操作をすべ
ての直線L1〜Ln、線分M1〜Mmに関して行なうこ
とにより、直線LH−Lnによって囲まれる領域が残さ
れる。この残された領域は、すなわち放射線照射野Bで
ある。
演算部225は以上のようにして認識した放射線照射野
βを示す情報Stを、制御回路314の信号抽出部35
0に送る。信号抽出部350は先読み画像信号Spから
、この情報Stが示す領域についての信号のみを抽出し
てヒストグラム解析部351に送る。したがって該ヒス
トグラム解析部351におけるヒストグラム解析は、蓄
積性蛍光体シートl。
a上の実際に放射線が照射された領域のみに関して行な
われることになるので、前述の設定値a1bおよびCは
、実際の蓄積記録情報に対して最適のものとなる。
以上説明した実施例においては、微分処理の方向の起点
となる点を照射野内のシート中心0としているか、この
起点はシート中心点に限らず、シート上に存在する点な
らばどのような点が利用されてもよい。例えば放射線照
射野が極めて小さく絞られる場合は、シート中心点が照
射野外に位置することもあるので、その場合は蓄積性蛍
光体シート内の濃度最大点、濃度重心点、さらには画像
濃度を2値化した際の高濃度側領域の重心等、照射野内
に存在することになる点を利用することができる。
また上記実施例では、微分処理の方向DI −Dnをシ
ート中心Oのまわりに等角度間隔で設定しているが、こ
れらの方向は特に等角度間隔に設定されなくても構わな
い。すなわち例えば第8図に示すように蓄積性蛍光体シ
ート103の辺部に等距離間隔の点を複数設定し、照射
野B内の点Pからそれらの点に向かう各方向D1〜Dn
を微分処理の方向とするようにしてもよい。
また第9図に示すように、照射野B内の点Pから輪郭候
補点Eまでの距離gがさほど変化しない所では微分処理
の方向りを比較的粗く設定しく図中h1の範囲)、上記
距離gがかなり変化するようになったら微分処理の方向
りを比較的精細に設定する(図中h2の範囲)ようにし
てもよい。
なお、正誤の判定に供する輪郭候補点から画像内の所定
点までの距離と、それ以外の輪郭候補点から上記所定点
までの距離との関係に基づいて正誤を判定するには、前
記(1)式以外の判定基準を用いることもできる。例え
ば照射野Bの形状が第14図に示すように矩形であるこ
とが前提として分がっているような場合は、照射野各辺
の最端部の輪郭候補点以外の輪郭候補点Emについて、
とならなければ、その輪郭候補点Emを誤検出点である
と判定することができる。すなわち、上述のような位置
の輪郭候補点Emは正しいものであれば、その両隣の輪
郭候補点Em+lとEm−1とを結ぶ線分よりシート端
部側に飛び出してしまうことは無いからである。
また照射野Bの形状が第15図に示すように円形である
場合は、正誤判定に供する輪郭候補点Emから照射野中
心Oまでの距離をzlその他のすべての輪郭候補点から
照射野中心Oまでの距離の平均値を免、所定のしきい値
をαとして、免−免1くα なる関係が成立しなければ、つまり距離λが他の輪郭候
補点から照射野中心Oまでの距離と大きくかけ離れてい
れば、その輪郭候補点Emは誤検出点であると判定する
ことができる。
また、照射野形状が上記と同様の場合は、判定に供する
輪郭候補点Emの両隣の輪郭候補点ETIl+1および
E m−1から照射野中心Oまでの距離をそれぞれ9,
1.9J2、所定のしきい値をβとして、免−(9Jt
 +9.z ) /21 <βなる関係が成立しなけれ
ば、つまり距離えが距離9.1と9.zの平均値と大き
くかけ離れていれば、その輪郭候補点Enlは誤検出点
であると判定することができる。
また本発明においては、求められた輪郭候補点間の距離
を判定基準に用いることもできる。例えば第16図に示
すように照射野Bの形状が円形である場合は、判定に供
する輪郭候補点Emとその両隣の輪郭候補点Eill+
lsEm→までの各距離をそれぞれQ13.9Ja 、
所定のしきい値をγとして、(51!、3+免ム)/2
くγ なる関係が成立しなければ、つまりこれらの距離Q、a
 、Q、、の平均値が特異的に大きくなっていれば、そ
の輪郭候補点Emは誤検出点であると判定することがで
きる。
さらに本発明においては、判定に供する輪郭候補点とそ
れ以外の輪郭候補点2つとをそれぞれ結ぶ2本の線分が
なす角度を判定基準として用いることもできる。例えば
第17図に示すように照射野Bが円形であって、判定に
供する輪郭候補点Emとその両隣の輪郭候補点E m+
I SE m−lを結ぶ2本の線分Jt 、Jzを考え
、これらの線分Jt 、Jzがなす角度をθ、所定のし
きい値をφとして、θ〉φ なる関係が成立しなければ、つまりこの角度θが特異的
に小さくなっていれば、その輪郭候補点En1は誤検出
点であると判定することができる。また上記のような角
度θをすべての輪郭候補点について求め、その平均値を
θ、所定のしきい値をδとしたとき、   1θ−θ1
くδ なる関係が成立しなければ、すなわちこの角度θが平均
値と大きくかけ離れていれば、その輪郭候補点Emは誤
検出点であると判定することができる。
次に第18図を参照して、輪郭候補点正誤判定の精度を
より高めるようにした本発明の別の方法について説明す
る。この第18図は、先に説明した第5図図示の照射野
認識回路220と同じように使用される照射野認識回路
250を示している。なおこの第18図において、第5
図中の要素と同等のものには同番号を付してあり、それ
らについては特に必要の無い限り説明を省略する。この
回路250において、輪郭候補点信号検出部223が出
力する輪郭候補点画素位置を示す情報Smは、正誤判定
部251に送られる。またこの正誤判定部251には、
先読み画像信号Spも入力される。本例においては、前
記第14図に示したように照射野Bの形状は、必ず矩形
とされるようになっており、正誤判定部251は先に述
べた(2式を判定基準として、該(2式が満たされなけ
れば、輪郭候補点Emがその両隣の輪郭候補点Em+I
とE m−1とを結ぶ線分よりも外側に飛び出している
とみなす。そして正誤判定部251は、このように輪郭
候補点Emが飛び出していると見なした場合、その両隣
の輪郭候補点En’?)lおよびE m−1に基づいて
照射野輪郭上に位置すると導き出される点Enl″を求
める。本例でこの点Emは、輪郭候補点Emと微分処理
の起点(シート中心)0を結ぶ線分上にあって、該起点
Oからだけ離れた点とされる(第19図参照)。
次いで正誤判定部251は、先読み画像信号Spに基づ
いて、上記点Em° と最初に検出した輪郭候補点Em
との間の平均画像濃度りを求め、この平均向1象濃度り
が所定のしきい値を下回る場合は、輪郭候補点Emが誤
検出点であると判定する。すなわち、輪郭候補点Emか
前述のように飛び出した点であれば、両点EmとEm゛
との間は照射野外領域であって、この部分はいかなる場
合も概して低濃度であるから、上述のように判定できる
のである。こうして輪郭候補点Emが照射野外まで飛び
出した点であると判定した場合、正誤判定部251は、
前述したようにその点Emをキャンセルしたり、あるい
は前記輪郭候補点En++およびE m−1に基づいて
導き出された点Ent゛に置き換える。
上述の正誤判定を受けた後の輪郭候補点の画素位置を示
す情報Se’ は演算部225に送られ、前述の場合と
同様にして、この情報Se’ に基づいて照射野輪郭が
求められる。
なお、検出した輪郭候補点が照射野外にあるものである
か否かを、以上述べたような平均画像濃度に基づいて判
定するステップは、特に前記(2)式を判定基準として
正誤判定する場合に限らず、その他の判定基■で輪郭候
補点の正誤判定を行なう場合においても付加的に実施さ
れうるちのである。
また以上説明したような「先読み」は、通常「本読み」
におけるよりも粗い画素単位で行なわれる。前述の微分
処理は、このような比較的粗い読取り操作によって得ら
れた画像データそのものに対して行なってもよいし、こ
れらの画像データを補間してより精細な画像データを得
てからそれらの画像データに対して行なってもよい。さ
らには、複数画素の画像信号を平均した画像データに対
して上記微分処理を行なうようにしても構わない。
また、第2図に示すような形状の照射野Bを認識する場
合、照射野輪郭部の輪郭候補点は通常1回の微分処理で
1つだけ求められるが、例えば第10図、第11図に示
すような形状の照射野Bを認識する場合は、照射野輪郭
部の輪郭候補点が複数検出されることもある。このよう
な場合でも、前述の差分の値がしきい値を超える点はす
べて輪郭候補点として検出するようにしておけば、照射
野輪郭部の輪郭候補点をすべて検出可能で、複雑な形状
の照射野Bも正しく認識できることになる。また例えば
照射野Bの形状が必ず矩形であると決まっているような
場合は、第12図図示のように、微分処理の起点となる
点Pを移動させつつ微分処理を行なうようにしてもよい
。このような場合は蓄積性蛍光体シート103の左端側
の照射野輪郭部と、右端側の照射野輪郭部とで差分の符
号(正負)が反対になるが、いずれの場合も、差分の絶
対値としきい値とを比較すればよい。本発明による輪郭
候補点正誤判定方法は、このようにして輪郭候補点を求
める場合においても当然実施可能である。
尚、第1図に示される装置は、本読み用読取部と先読み
用読取部とを個別に有しているが、例えば特開昭58−
67242号に示されるように本読み用読取系と先読み
用読取系とを兼用し、先読みが終了したならばシート移
送手段により蓄積性蛍光体シートを読取系に戻して本読
みを行ない、先読み時には励起光エネルギー調整手段に
より、励起光エネルギーか本読み時のそれよりも小さく
なるように調整してもよく、本発明方法はそのような装
置によって放射線画像情報読取りを行なう場合において
も適用可能である。
さらに、以上述べた実施例においては、先読み画像信号
から放射線照射野を認識するようにしているが、本読み
画像信号、または先読みを行なわずに本読みに相当する
読取りを直接行なって得た画像信号を利用して本発明方
法により照射野輪郭候補点を認識することも勿論可能で
ある。このような場合は、認識した照射野情報を、例え
ば画像処理条件設定値Cを適正に設定するために利用す
ることができる。
さらに、以上述べた実施例においては、蓄積性蛍光体シ
ートを放射線画像情報の記録媒体として利用しているが
、本発明方法は、従来から知られているX線撮影用銀塩
写真フィルムから放射線画像を読み取って画像信号を得
る場合においても、同様に実施されうるちのである。
(発明の効果) 以上詳細に説明した通り本発明の放射線照射野輪郭候補
点正誤判定方法においては、微分処理などによって検出
された輪郭候補点が所定の判定基準を満足しているか否
かを調べ、さらには、正誤判定に供する輪郭候補点から
、それ以外の輪郭候補点に基づいて照射野輪郭上にある
と想定される点までの平均画像濃度に基づいて、該輪郭
候補点が照射野外領域にあるものであるか否かを調べる
ようにしているので、輪郭候補点の正誤を正確に判定可
能となっている。したがって本方法によれば、誤った輪
郭候補点を照射野認識のために供することを防止し、放
射線照射野を正確に認識して、被写体に関する情報を正
しく把握し、本読みの読取条件や画像処理条件を最適に
設定することができる。よって本発明方法によれば、常
に観察読影適性の優れた放射線画像を再生することが可
能となる。
【図面の簡単な説明】
第1図は、本発明方法により照射野輪郭候補点の正誤を
判定して放射線画像情報読取りを行なう装置の概略構成
図、 第2図は本発明に係る蓄積性蛍光体シートへの放射線画
像情報記録状態を示す説明図、第3図は本発明に係る微
分処理の方向を説明する説明図、 第4図は本発明に係る画像信号の分布状態と画像信号差
分値の分布状態を示すグラフ、第5図は第1図の装置の
一部を詳しく示すブロック図、 第6図は輪郭候補点に沿った直線を検出する方法を説明
するためのグラフ、 第7図は輪郭候補点に沿った直線で囲まれる領域を抽出
する方法を説明するための説明図、第8.9.10.1
1および12図は、本発明に係る微分処理の方向設定の
別の例を示す説明図、第13.14.15.16および
17図は、本発明方法における判定基準を説明する説明
図、 第18図は本発明方法を実施する装置の別の例を示すブ
ロック図、 第19図は本発明方法における平均画像濃度による輪郭
候補点正誤判定を説明する説明図である。 20・・・放射線画像撮影部 40・・・本読み用読取部 lOl・・・被写体 103・・・蓄積性蛍光体シート 201・・・先読み用レーザ光源 202・・・先読み用レーザ光 204・・・先読み用光偏向器 208・・・先読み用光検出器 210・・・先読み用シート移送手段 220.250・・・照射野認識回路 222・・・しきい値設定部 223・・・輪郭候補点信号検出部 224.251・・・正誤判定部 301・・・本読み用レーザ光源 302・・・本読み用レーザ光 305・・・本読み用光偏向器 310・・・本読み用光検出器  311・・・増幅器
312・・・A/D変換器    313・・・信号処
理回路314・・・制御回路 320・・・本読み用シ
ート移送手段B・・・放射線照射野    a・・・読
取ゲイン設定値221・・・微分処理部 225・・・演算部 30・・・先読み用読取部 100・・・放射線源 102・・・放射線 104・・・絞り b・・・収録スケールファクター設定値C・・・再生画
像処理条件設定値 D1〜Dn・・・微分処理の方向 0、P・・・微分処理の起点となる点 So・・・本読み画像信号  Sp・・・先読み画像信
号Se・・・輪郭候補点を示す情報

Claims (2)

    【特許請求の範囲】
  1. (1)照射野絞りをかけて放射線が照射されて放射線画
    像情報が記録された記録媒体から前記画像情報を読み取
    り、この画像情報に基づいて、照射野の輪郭上にあると
    考えられる輪郭候補点を検出する方法において、 検出された輪郭候補点が所定の判定基準を満たしている
    か否かを調べ、満たしていない場合は、該輪郭候補点が
    前記輪郭上に位置しない誤検出の点であると判定するこ
    とを特徴とする放射線照射野輪郭候補点正誤判定方法。
  2. (2)照射野絞りをかけて放射線が照射されて放射線画
    像情報が記録された記録媒体から前記画像情報を読み取
    り、この画像情報に基づいて、照射野の輪郭上にあると
    考えられる輪郭候補点を検出する方法において、 検出された輪郭候補点が所定の判定基準を満たしている
    か否かを調べるとともに、 この輪郭候補点から、それ以外の輪郭候補点に基づいて
    照射野輪郭上にあると想定される点までの平均画像濃度
    を求め、 前記輪郭候補点が前記判定基準を満たさず、かつ前記平
    均画像濃度が所定のしきい値を下回る場合は、この輪郭
    候補点が前記輪郭上に位置しない誤検出の点であると判
    定することを特徴とする放射線照射野輪郭候補点正誤判
    定方法。
JP63186084A 1988-03-19 1988-07-26 放射線照射野輪郭候補点正誤判定方法 Expired - Fee Related JP2582625B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP63186084A JP2582625B2 (ja) 1988-03-19 1988-07-26 放射線照射野輪郭候補点正誤判定方法
US07/324,879 US4992663A (en) 1988-03-19 1989-03-17 Method of judging the correctness or incorrectness of a prospective contour point of an irradiation field
DE8989104828T DE68902255T2 (de) 1988-03-19 1989-03-17 Verfahren zur bestimmung der kontur eines bestrahlungsfelds.
EP89104828A EP0335204B1 (en) 1988-03-19 1989-03-17 Method for determining the contour of an irradiation field

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP6673788 1988-03-19
JP63-66737 1988-03-19
JP63186084A JP2582625B2 (ja) 1988-03-19 1988-07-26 放射線照射野輪郭候補点正誤判定方法

Publications (2)

Publication Number Publication Date
JPH021080A true JPH021080A (ja) 1990-01-05
JP2582625B2 JP2582625B2 (ja) 1997-02-19

Family

ID=26407934

Family Applications (1)

Application Number Title Priority Date Filing Date
JP63186084A Expired - Fee Related JP2582625B2 (ja) 1988-03-19 1988-07-26 放射線照射野輪郭候補点正誤判定方法

Country Status (1)

Country Link
JP (1) JP2582625B2 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6853740B1 (en) 1998-04-07 2005-02-08 Canon Kabushiki Kaisha Image processing method apparatus and storage medium for recognition of irradiation area
WO2010041731A1 (ja) * 2008-10-10 2010-04-15 日本電気株式会社 照合装置、照合方法、及びプログラム
CN115546107A (zh) * 2022-09-06 2022-12-30 西安大医集团股份有限公司 射野识别方法、计算机设备和计算机可读存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS58225488A (ja) * 1982-06-25 1983-12-27 Hitachi Ltd 特徴点自動抽出方式
JPS60179881A (ja) * 1984-02-15 1985-09-13 Komatsu Ltd 輪郭線の認識方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS58225488A (ja) * 1982-06-25 1983-12-27 Hitachi Ltd 特徴点自動抽出方式
JPS60179881A (ja) * 1984-02-15 1985-09-13 Komatsu Ltd 輪郭線の認識方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6853740B1 (en) 1998-04-07 2005-02-08 Canon Kabushiki Kaisha Image processing method apparatus and storage medium for recognition of irradiation area
US7046836B2 (en) 1998-04-07 2006-05-16 Canon Kabushiki Kaisha Image processing method apparatus and storage medium for recognition of irradiation area
WO2010041731A1 (ja) * 2008-10-10 2010-04-15 日本電気株式会社 照合装置、照合方法、及びプログラム
US8712190B2 (en) 2008-10-10 2014-04-29 Nec Corporation Collating device, collating method, and program
JP5556663B2 (ja) * 2008-10-10 2014-07-23 日本電気株式会社 照合装置、照合方法、及びプログラム
CN115546107A (zh) * 2022-09-06 2022-12-30 西安大医集团股份有限公司 射野识别方法、计算机设备和计算机可读存储介质

Also Published As

Publication number Publication date
JP2582625B2 (ja) 1997-02-19

Similar Documents

Publication Publication Date Title
US4967079A (en) Method of recognizing irradiation field
US4962539A (en) Method of recognizing layout and subdivision patterns of radiation images
JPH0584505B2 (ja)
JPS63257879A (ja) 放射線画像の分割パタ−ン認識方法
US4992663A (en) Method of judging the correctness or incorrectness of a prospective contour point of an irradiation field
US4995093A (en) Radiation field recognition method
JPH021080A (ja) 放射線照射野輪郭候補点正誤判定方法
JP2525648B2 (ja) 放射線照射野輪郭候補点正誤判定方法
US5533142A (en) Method for adjusting read-out and processing conditions for Magen images
JPH0279177A (ja) 放射線画像における乳房輪郭検出方法
JPH0470645A (ja) 放射線画像の分割パターン認識方法および装置
JPH021070A (ja) 照射野絞り有無判定方法
JPH01238654A (ja) 放射線画像の分割パターン認識方法
JPH02272682A (ja) 被写体像内画像点決定方法
JP2761691B2 (ja) 放射線照射野認識方法
JPH041745A (ja) 画像記録位置認識装置
JPH0470644A (ja) 放射線画像の照射野検出方法および装置
JPS63100437A (ja) 照射野認識方法
US4864133A (en) Method of adjusting radiation image read-out conditions
JP2764492B2 (ja) 放射線画像の撮影方向認識方法
JPH0469647A (ja) 放射線画像の分割パターン認識方法
JPH05192320A (ja) 放射線画像の照射野検出方法
JPH0271249A (ja) 放射線照射野判定方法
JPH0266682A (ja) 照射野輪郭候補点正誤判定方法
JPH01238638A (ja) 放射線照射野輪郭候補点検出方法

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees