JPH10275213A - 放射線画像の照射野認識方法および装置、並びに黒化処理方法および装置 - Google Patents
放射線画像の照射野認識方法および装置、並びに黒化処理方法および装置Info
- Publication number
- JPH10275213A JPH10275213A JP9079696A JP7969697A JPH10275213A JP H10275213 A JPH10275213 A JP H10275213A JP 9079696 A JP9079696 A JP 9079696A JP 7969697 A JP7969697 A JP 7969697A JP H10275213 A JPH10275213 A JP H10275213A
- Authority
- JP
- Japan
- Prior art keywords
- irradiation field
- candidate
- candidate line
- lines
- reference candidate
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/25—Determination of region of interest [ROI] or a volume of interest [VOI]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S128/00—Surgery
- Y10S128/92—Computer assisted medical diagnostics
- Y10S128/922—Computer assisted medical diagnostics including image analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
(57)【要約】
【課題】 放射線画像の照射野認識装置において、照射
野認識の精度をより向上させる。 【解決手段】 エッジ候補点検出手段10により、照射野
の認識対象である画像中の所定の点について設定した放
射状の直線方向に沿ってエッジ候補点を検出し、これら
のエッジ候補点について基準候補線検出手段20によりハ
フ（ Hough）変換を適用して所定数の基準候補線を求
め、これらの基準候補線で囲まれる領域を照射野決定手
段30により照射野領域として定める。
野認識の精度をより向上させる。 【解決手段】 エッジ候補点検出手段10により、照射野
の認識対象である画像中の所定の点について設定した放
射状の直線方向に沿ってエッジ候補点を検出し、これら
のエッジ候補点について基準候補線検出手段20によりハ
フ（ Hough）変換を適用して所定数の基準候補線を求
め、これらの基準候補線で囲まれる領域を照射野決定手
段30により照射野領域として定める。
Description
【０００１】
【発明の属する技術分野】本発明は放射線画像の照射野
認識方法および認識装置並びに黒化処理方法および黒化
処理装置に関するものである。
認識方法および認識装置並びに黒化処理方法および黒化
処理装置に関するものである。
【０００２】
【従来の技術】記録された放射線画像を読み取って画像
データを得、この画像データに適切な画像処理を施した
後、処理済みの画像データに基づいて読影に適した可視
像を再生することは種々の分野で行なわれている。例え
ば、本出願人による蓄積性蛍光体シートを利用した放射
線画像記録再生システムもその一つである。
データを得、この画像データに適切な画像処理を施した
後、処理済みの画像データに基づいて読影に適した可視
像を再生することは種々の分野で行なわれている。例え
ば、本出願人による蓄積性蛍光体シートを利用した放射
線画像記録再生システムもその一つである。
【０００３】ところでＸ線フイルムや蓄積性蛍光体シー
トに放射線画像を撮影記録するに際しては、放射線の照
射による生体への影響を極力小さくする観点、および例
えば観察に不要な部分からの散乱光による画質性能の低
下等を防止する観点から、放射線が被写体の必要な部分
にのみ照射されるように照射域を制限する鉛などで作ら
れた照射野絞りを使用することが多い。
トに放射線画像を撮影記録するに際しては、放射線の照
射による生体への影響を極力小さくする観点、および例
えば観察に不要な部分からの散乱光による画質性能の低
下等を防止する観点から、放射線が被写体の必要な部分
にのみ照射されるように照射域を制限する鉛などで作ら
れた照射野絞りを使用することが多い。
【０００４】照射野絞りを用いて撮影を行なった場合、
蓄積性蛍光体シート等の記録媒体には、照射野絞りの開
口輪郭の内部領域（照射野領域）に被写体等の画像が記
録され、開口輪郭の外側領域（照射野外領域）には放射
線が到達せず未露光状態となる。したがってこの開口輪
郭に対応する画像の照射野輪郭はエッジ線となる。
蓄積性蛍光体シート等の記録媒体には、照射野絞りの開
口輪郭の内部領域（照射野領域）に被写体等の画像が記
録され、開口輪郭の外側領域（照射野外領域）には放射
線が到達せず未露光状態となる。したがってこの開口輪
郭に対応する画像の照射野輪郭はエッジ線となる。
【０００５】そして、このように照射野領域内にのみ画
像が記録された記録媒体から画像データを読み取って画
像処理を行なう場合、照射野領域内の画像データについ
てのみ上記階調処理等の処理を施せばよく、処理の回数
が大幅に低減され処理負荷の低減、処理スピードの向上
を図ることができる。
像が記録された記録媒体から画像データを読み取って画
像処理を行なう場合、照射野領域内の画像データについ
てのみ上記階調処理等の処理を施せばよく、処理の回数
が大幅に低減され処理負荷の低減、処理スピードの向上
を図ることができる。
【０００６】一方、照射野外領域は未露光状態であるた
め、医用Ｘ線フイルムのようなネガ画像においては最低
濃度領域となるが、例えば医用Ｘ線フイルムをシャーカ
スに掛けて蛍光灯の光による透過画像を観察するとき等
には、このような最低濃度領域は非常に明るい領域とな
るため、照射野領域のうち特に照射野外領域に近い部分
については当該明るい照射野外領域に眩惑されて観察に
支障のある場合がある。同様にそのような画像を電気信
号として読み取って得られた画像データを用いてＣＲＴ
等の画像表示装置にその画像を表示する場合には、当該
照射野外領域は最高輝度となるため必要な照射野内の画
像の読取りに支障を生じる。
め、医用Ｘ線フイルムのようなネガ画像においては最低
濃度領域となるが、例えば医用Ｘ線フイルムをシャーカ
スに掛けて蛍光灯の光による透過画像を観察するとき等
には、このような最低濃度領域は非常に明るい領域とな
るため、照射野領域のうち特に照射野外領域に近い部分
については当該明るい照射野外領域に眩惑されて観察に
支障のある場合がある。同様にそのような画像を電気信
号として読み取って得られた画像データを用いてＣＲＴ
等の画像表示装置にその画像を表示する場合には、当該
照射野外領域は最高輝度となるため必要な照射野内の画
像の読取りに支障を生じる。
【０００７】そこで放射線画像記録再生システムにおい
ては、このような照射野外領域についての各画像データ
を一律に最高濃度（若しくは最低輝度）に相当する値に
強制的に置換する処理が行なわれる。そしてこの処理は
一般に黒化処理と称されているが、この黒化処理を行な
う場合の空間的基準となる照射野輪郭を精度よく認識す
ることが非常に重要である。
ては、このような照射野外領域についての各画像データ
を一律に最高濃度（若しくは最低輝度）に相当する値に
強制的に置換する処理が行なわれる。そしてこの処理は
一般に黒化処理と称されているが、この黒化処理を行な
う場合の空間的基準となる照射野輪郭を精度よく認識す
ることが非常に重要である。
【０００８】すなわち、この照射野輪郭を実際のものよ
り小さく（すなわち実際よりも内側に）認識してしまう
と、本来照射野内のものとして観察したい画像部分が照
射野外のものとして黒化処理されて読み取れなくなり、
一方、実際のものより大きく（すなわち実際よりも外側
に）認識してしまうと、最低濃度（最高輝度）の部分が
残ってしまい、黒化処理を行う意義が失われてしまうの
である。
り小さく（すなわち実際よりも内側に）認識してしまう
と、本来照射野内のものとして観察したい画像部分が照
射野外のものとして黒化処理されて読み取れなくなり、
一方、実際のものより大きく（すなわち実際よりも外側
に）認識してしまうと、最低濃度（最高輝度）の部分が
残ってしまい、黒化処理を行う意義が失われてしまうの
である。
【０００９】このため、精度よく照射野（輪郭）の認識
処理を行う方法が種々研究されている。
処理を行う方法が種々研究されている。
【００１０】例えば、上記照射野輪郭が画像の濃度変化
が急峻に変化するエッジ線になることを利用して、画像
データの変化が急峻な部分を探索することによって、照
射野輪郭を求める方法がよく知られているが、このエッ
ジ線を求める具体的な方法として、画像の所定の点（た
とえば画像の中心点等）から画像端部に向かう放射状の
複数の直線を設定し、これらの各直線の方向に沿った画
像データに基づいて各方向ごとにデータの差分が大きい
エッジ候補点をそれぞれ検出し、これらのエッジ候補点
に基づいてそのエッジ線を得る方法等の技術が提案され
ている（特開昭63-259538 号等）。
が急峻に変化するエッジ線になることを利用して、画像
データの変化が急峻な部分を探索することによって、照
射野輪郭を求める方法がよく知られているが、このエッ
ジ線を求める具体的な方法として、画像の所定の点（た
とえば画像の中心点等）から画像端部に向かう放射状の
複数の直線を設定し、これらの各直線の方向に沿った画
像データに基づいて各方向ごとにデータの差分が大きい
エッジ候補点をそれぞれ検出し、これらのエッジ候補点
に基づいてそのエッジ線を得る方法等の技術が提案され
ている（特開昭63-259538 号等）。
【００１１】
【発明が解決しようとする課題】しかしながら、上述の
放射状の方向に沿って得られたエッジ候補点に基づいて
エッジ線、すなわち照射野輪郭を認識する方法では、必
ずしも実際の照射野輪郭に合致したものが得られるとは
限られず、この場合、上述したように必要な画像部分が
黒化処理され、または低濃度部分が残存する、という問
題が生じる。
放射状の方向に沿って得られたエッジ候補点に基づいて
エッジ線、すなわち照射野輪郭を認識する方法では、必
ずしも実際の照射野輪郭に合致したものが得られるとは
限られず、この場合、上述したように必要な画像部分が
黒化処理され、または低濃度部分が残存する、という問
題が生じる。
【００１２】本発明は上記事情に鑑みなされたものであ
って、従来の照射野認識処理よりも認識精度を向上させ
た放射線画像の照射野認識方法および装置を提供するこ
とを目的とするものである。
って、従来の照射野認識処理よりも認識精度を向上させ
た放射線画像の照射野認識方法および装置を提供するこ
とを目的とするものである。
【００１３】また本発明の他の目的は、従来の黒化処理
よりも処理の精度を向上させた放射線画像の黒化処理方
法および装置を提供することである。
よりも処理の精度を向上させた放射線画像の黒化処理方
法および装置を提供することである。
【００１４】
【課題を解決するための手段】本発明の放射線画像の照
射野認識方法および照射野認識装置は、まず画像中の所
定の点について設定した放射状の直線方向に沿ってエッ
ジ候補点を検出し、これらのエッジ候補点についてハフ
（ Hough）変換を適用して所定数の基準候補線を求め、
これらの基準候補線で囲まれる領域を照射野領域として
定めることを特徴とするものである。
射野認識方法および照射野認識装置は、まず画像中の所
定の点について設定した放射状の直線方向に沿ってエッ
ジ候補点を検出し、これらのエッジ候補点についてハフ
（ Hough）変換を適用して所定数の基準候補線を求め、
これらの基準候補線で囲まれる領域を照射野領域として
定めることを特徴とするものである。
【００１５】すなわち本発明の照射野認識方法は、照射
野絞りを使用して撮影された照射野を有する放射線画像
の該照射野を、該放射線画像を表す画像データに基づい
て認識処理する照射野認識方法において、前記照射野内
に含まれる所定の点から前記放射線画像の端部にそれぞ
れ向う複数の放射状の直線を設定し、該設定された各放
射状の直線ごとに、該各直線上における前記照射野の輪
郭と考えられる各エッジ候補点を、該各直線に沿った前
記画像データに基づいて求め、前記放射線画像面内に直
交座標系における前記各エッジ候補点の座標を（ｘi，
ｙi ）としたとき、これらのｘi ，ｙi を定数として次
式（１）で表わされる曲線を各エッジ候補点ごとに求
め、 ρ＝ｘi cos θ＋ｙi sin θ （１） 該求められた曲線同士の各交点（ρj ，θj ）を求める
とともに、該各交点（ρj ，θj ）ごとの交わる前記曲
線の数（ボーディング値）をカウントし、このカウント
値の大きい順に所定の数の前記交点（ρj ，θj ）を抽
出し、該抽出された各交点（ρj ，θj ）に対応する、
次式（２）で規定される前記直交座標系（ｘｙ座標系）
上の各基準候補線を求め、 ρj ＝ｘ cosθj ＋ｙ sinθj （２） 前記式（２）で規定される複数の基準候補線によって囲
まれる領域を前記照射野として定めることを特徴とする
ものである。
野絞りを使用して撮影された照射野を有する放射線画像
の該照射野を、該放射線画像を表す画像データに基づい
て認識処理する照射野認識方法において、前記照射野内
に含まれる所定の点から前記放射線画像の端部にそれぞ
れ向う複数の放射状の直線を設定し、該設定された各放
射状の直線ごとに、該各直線上における前記照射野の輪
郭と考えられる各エッジ候補点を、該各直線に沿った前
記画像データに基づいて求め、前記放射線画像面内に直
交座標系における前記各エッジ候補点の座標を（ｘi，
ｙi ）としたとき、これらのｘi ，ｙi を定数として次
式（１）で表わされる曲線を各エッジ候補点ごとに求
め、 ρ＝ｘi cos θ＋ｙi sin θ （１） 該求められた曲線同士の各交点（ρj ，θj ）を求める
とともに、該各交点（ρj ，θj ）ごとの交わる前記曲
線の数（ボーディング値）をカウントし、このカウント
値の大きい順に所定の数の前記交点（ρj ，θj ）を抽
出し、該抽出された各交点（ρj ，θj ）に対応する、
次式（２）で規定される前記直交座標系（ｘｙ座標系）
上の各基準候補線を求め、 ρj ＝ｘ cosθj ＋ｙ sinθj （２） 前記式（２）で規定される複数の基準候補線によって囲
まれる領域を前記照射野として定めることを特徴とする
ものである。
【００１６】ここで上記放射状の直線の起点となる所定
の点としては、例えば画像の物理的な中心の点が選択さ
れる。画像の中心部に照射野外領域を配置することは一
般的でないためである。ただし本発明においては、常に
中心の点とする必要はなく、これ以外の他の点を用いる
ようにしてもよい。
の点としては、例えば画像の物理的な中心の点が選択さ
れる。画像の中心部に照射野外領域を配置することは一
般的でないためである。ただし本発明においては、常に
中心の点とする必要はなく、これ以外の他の点を用いる
ようにしてもよい。
【００１７】また放射状の各直線方向に沿ったエッジ候
補点は、その方向線上において互いに隣接する２つの画
素の画像データの差分を順次求め、その差分が最大とな
るときの両画素について、その画像データ値の小さい方
の画素をエッジ候補点とすればよい。
補点は、その方向線上において互いに隣接する２つの画
素の画像データの差分を順次求め、その差分が最大とな
るときの両画素について、その画像データ値の小さい方
の画素をエッジ候補点とすればよい。
【００１８】なお上記放射線画像としては、読み取った
画像そのもの（原画像）の他、濃度・階調等が所定の範
囲となるように一定の規格化処理を施した画像、原画像
から画素数を間引いて得られる縮小画像、および原画像
を得るためのいわゆる本読みに先立って本読みよりも画
素数を減らして粗く読み取って得られた先読み画像等種
々の画像を適用することができる。
画像そのもの（原画像）の他、濃度・階調等が所定の範
囲となるように一定の規格化処理を施した画像、原画像
から画素数を間引いて得られる縮小画像、および原画像
を得るためのいわゆる本読みに先立って本読みよりも画
素数を減らして粗く読み取って得られた先読み画像等種
々の画像を適用することができる。
【００１９】また、式（２）で規定される複数の基準候
補線によって囲まれる領域を照射野として定めるに当た
って、このような閉領域が２以上形成されることも考え
られるが、この場合、１つの閉領域を除いて他の閉領域
の面積は非常に小さいものとなる。したがってこのよう
な場合は各閉領域の面積を求め、各閉領域の面積を比較
することによって１つの照射野を定めることができる。
補線によって囲まれる領域を照射野として定めるに当た
って、このような閉領域が２以上形成されることも考え
られるが、この場合、１つの閉領域を除いて他の閉領域
の面積は非常に小さいものとなる。したがってこのよう
な場合は各閉領域の面積を求め、各閉領域の面積を比較
することによって１つの照射野を定めることができる。
【００２０】基準候補線を抽出する本数を決定付ける上
記所定の数は、少なくとも照射野輪郭を構成する直線の
辺の数以上であることが必要である。
記所定の数は、少なくとも照射野輪郭を構成する直線の
辺の数以上であることが必要である。
【００２１】以上のように式（２）により規定された各
基準候補線に関し、さらに所定の評価関数を用いた評価
を行ない、この評価結果に基づいて前記基準候補線のう
ちから又は前記基準候補線に代えて、所定の候補線を求
め、上記基準候補線に代えて該候補線によって囲まれる
領域を前記照射野として定めるようにしてもよい。
基準候補線に関し、さらに所定の評価関数を用いた評価
を行ない、この評価結果に基づいて前記基準候補線のう
ちから又は前記基準候補線に代えて、所定の候補線を求
め、上記基準候補線に代えて該候補線によって囲まれる
領域を前記照射野として定めるようにしてもよい。
【００２２】そのような所定の評価関数を用いた評価と
しては、第１に、式（２）で規定された各基準候補線上
に、照射野輪郭上の点として予め設定された所定の代表
点が含まれるか否かを該各基準候補線ごとに判定するこ
とを適用することができ、この場合、所定の候補線は基
準候補線のうち所定の代表点を含む基準候補線である。
しては、第１に、式（２）で規定された各基準候補線上
に、照射野輪郭上の点として予め設定された所定の代表
点が含まれるか否かを該各基準候補線ごとに判定するこ
とを適用することができ、この場合、所定の候補線は基
準候補線のうち所定の代表点を含む基準候補線である。
【００２３】照射野輪郭上の点として予め設定された所
定の代表点とは、実際の画像における照射野輪郭上に存
在する蓋然性の高い点を意味し、例えば前述した画像デ
ータの差分が極めて大きい値を示す等の点が該当する。
定の代表点とは、実際の画像における照射野輪郭上に存
在する蓋然性の高い点を意味し、例えば前述した画像デ
ータの差分が極めて大きい値を示す等の点が該当する。
【００２４】所定の評価関数を用いた評価としては、第
２に、式（２）で規定された各基準候補線を、該基準候
補線を中心とする±ｍ画素、±α度の範囲内でシフトお
よび／または回転させて、｛（２ｍ＋１）（２α＋１）
−１｝個の変換候補線を設定し、前記基準候補線および
前記各変換候補線ごとに、各候補線に関する微分値の、
各候補線の延在方向についての平均値を求めることであ
り、前記所定の候補線は、基準候補線に代えて、この微
分値の平均値が大きい順に選択された所定数の基準候補
線および変換候補線である。
２に、式（２）で規定された各基準候補線を、該基準候
補線を中心とする±ｍ画素、±α度の範囲内でシフトお
よび／または回転させて、｛（２ｍ＋１）（２α＋１）
−１｝個の変換候補線を設定し、前記基準候補線および
前記各変換候補線ごとに、各候補線に関する微分値の、
各候補線の延在方向についての平均値を求めることであ
り、前記所定の候補線は、基準候補線に代えて、この微
分値の平均値が大きい順に選択された所定数の基準候補
線および変換候補線である。
【００２５】なお上記ｍの値およびθの値は、例えばｍ
＝１（画素），θ＝１／（ｘ2 ＋ｙ2 ）1/2 （度）によ
り設定するのが好ましい。ここで（ｘ，ｙ）は、縮小画
像のサイズを画素数（ｘ画素×Ｙ画素）で表したもので
あり、例えば縮小画像のサイズが 100画素×80画素の場
合は、ｘ＝ 100（画素），Ｙ＝80（画素）より、ｍ＝１
（画素），θ＝１／（ 1002 ＋802 ）1/2 ＝ 7.8×10-3
（度）である。
＝１（画素），θ＝１／（ｘ2 ＋ｙ2 ）1/2 （度）によ
り設定するのが好ましい。ここで（ｘ，ｙ）は、縮小画
像のサイズを画素数（ｘ画素×Ｙ画素）で表したもので
あり、例えば縮小画像のサイズが 100画素×80画素の場
合は、ｘ＝ 100（画素），Ｙ＝80（画素）より、ｍ＝１
（画素），θ＝１／（ 1002 ＋802 ）1/2 ＝ 7.8×10-3
（度）である。
【００２６】また、縮小率＝１／18〜１倍（拡大率＝１
〜18倍）程度の範囲内であれば、ｍの値およびθの値と
して、縮小率の逆数の値（すなわち拡大率と同一の値）
をそのまま適用することができる。例えば縮小率が１／
９倍の場合は、その逆数である９をそのまま適用して、
ｍ＝９（画素），θ＝９（度）とすることもできる。
〜18倍）程度の範囲内であれば、ｍの値およびθの値と
して、縮小率の逆数の値（すなわち拡大率と同一の値）
をそのまま適用することができる。例えば縮小率が１／
９倍の場合は、その逆数である９をそのまま適用して、
ｍ＝９（画素），θ＝９（度）とすることもできる。
【００２７】上記評価関数において、各候補線に関する
微分値とは、各候補線を挟んで相対向する２つの画素間
の画像データの差分を意味する。この差分値が大きい場
合は、この２つの画素間で濃度（または輝度）の差が大
きいことを意味し、画像のエッジが存在することを意味
する。そして図９（１）に示すように、このような候補
線を挟んで相対向する２つの画素の組は、候補線の延び
る方向にそれぞれ存在するため、上記微分値を、候補線
の延びる方向にそれぞれ存在する画素の組について加算
し、その和を加算した組数で除することにより上記平均
値を求める。この平均値が大きい程、当該候補線が放射
線画像の現実のエッジである可能性が高く、最も濃度差
（または輝度差）のあるエッジが照射野輪郭を構成する
線である。
微分値とは、各候補線を挟んで相対向する２つの画素間
の画像データの差分を意味する。この差分値が大きい場
合は、この２つの画素間で濃度（または輝度）の差が大
きいことを意味し、画像のエッジが存在することを意味
する。そして図９（１）に示すように、このような候補
線を挟んで相対向する２つの画素の組は、候補線の延び
る方向にそれぞれ存在するため、上記微分値を、候補線
の延びる方向にそれぞれ存在する画素の組について加算
し、その和を加算した組数で除することにより上記平均
値を求める。この平均値が大きい程、当該候補線が放射
線画像の現実のエッジである可能性が高く、最も濃度差
（または輝度差）のあるエッジが照射野輪郭を構成する
線である。
【００２８】所定の評価関数を用いた評価としては、第
３に、式（２）で規定された各基準候補線を、該基準候
補線を中心とする±ｍ画素、±α度の範囲内でシフトお
よび／または回転させて、｛（２ｍ＋１）（２α＋１）
−１｝個の変換候補線を設定し、前記基準候補線および
前記各変換候補線ごとに、各候補線に関する濃度勾配ベ
クトルの向きの、各候補線の延在方向についてのエント
ロピーを求めることであり、前記所定の候補線は、前記
基準候補線に代えて、前記エントロピーが小さい順に選
択された所定数の基準候補線および変換候補線である。
３に、式（２）で規定された各基準候補線を、該基準候
補線を中心とする±ｍ画素、±α度の範囲内でシフトお
よび／または回転させて、｛（２ｍ＋１）（２α＋１）
−１｝個の変換候補線を設定し、前記基準候補線および
前記各変換候補線ごとに、各候補線に関する濃度勾配ベ
クトルの向きの、各候補線の延在方向についてのエント
ロピーを求めることであり、前記所定の候補線は、前記
基準候補線に代えて、前記エントロピーが小さい順に選
択された所定数の基準候補線および変換候補線である。
【００２９】上記評価関数において、各候補線に関する
濃度勾配ベクトルとは、各候補線の一方の側の各画素か
ら濃度勾配（画像データの差分値）が最も大きい方向を
向くベクトルを意味する。そして、この濃度勾配ベクト
ルの向いた方向を示す指標値（例えば候補線に直交する
方向となす角度βの正弦値（ sinβ））を求め、図９
（２）に示すように、候補線の延びる方向にそれぞれ存
在する画素からの各濃度勾配ベクトルの向きについて、
上記指標値ごとのヒストグラムを作成し、このヒストグ
ラムの確率密度Piに基づいてエントロピー（−Σ{(Pi)l
og Pi}）を求める。このエントロピーが小さい程、濃度
勾配ベクトルの向きが一定の方向に揃っていることを意
味し、したがってそのような候補線が放射線画像の現実
の照射野輪郭を構成する線である可能性が高いものであ
る。
濃度勾配ベクトルとは、各候補線の一方の側の各画素か
ら濃度勾配（画像データの差分値）が最も大きい方向を
向くベクトルを意味する。そして、この濃度勾配ベクト
ルの向いた方向を示す指標値（例えば候補線に直交する
方向となす角度βの正弦値（ sinβ））を求め、図９
（２）に示すように、候補線の延びる方向にそれぞれ存
在する画素からの各濃度勾配ベクトルの向きについて、
上記指標値ごとのヒストグラムを作成し、このヒストグ
ラムの確率密度Piに基づいてエントロピー（−Σ{(Pi)l
og Pi}）を求める。このエントロピーが小さい程、濃度
勾配ベクトルの向きが一定の方向に揃っていることを意
味し、したがってそのような候補線が放射線画像の現実
の照射野輪郭を構成する線である可能性が高いものであ
る。
【００３０】なお、このエントロピーを評価関数として
用いる場合には、濃度勾配ベクトルの向きが一定の方向
に揃うか否かに基づく評価であるため、照射野輪郭を構
成する線が直線である場合にのみ適用でき、円弧や曲線
には適用できない。
用いる場合には、濃度勾配ベクトルの向きが一定の方向
に揃うか否かに基づく評価であるため、照射野輪郭を構
成する線が直線である場合にのみ適用でき、円弧や曲線
には適用できない。
【００３１】また上記濃度勾配ベクトルとは必ずしも画
像データが濃度値を示す場合にのみ適用しうる用語では
なく、濃度値や輝度値を含む濃淡や明暗を表す画像デー
タに基づく勾配ベクトルを広く含むものである。
像データが濃度値を示す場合にのみ適用しうる用語では
なく、濃度値や輝度値を含む濃淡や明暗を表す画像デー
タに基づく勾配ベクトルを広く含むものである。
【００３２】本発明の放射線画像の照射野認識装置は、
上記本発明の放射線画像の照射野認識方法を実施するた
めの装置であって、照射野絞りを使用して撮影された照
射野を有する放射線画像の該照射野を、該放射線画像を
表す画像データに基づいて認識処理する照射野認識装置
において、前記照射野内に含まれる所定の点から前記放
射線画像の端部にそれぞれ向う複数の放射状の直線を設
定し、該設定された各放射状の直線ごとに、該各直線上
における前記照射野の輪郭と考えられる各エッジ候補点
を、該各直線に沿った前記画像データに基づいて求める
エッジ候補点検出手段と、前記エッジ候補点検出手段に
より検出された前記各エッジ候補点の座標を（ｘi ，ｙ
i ）としたとき、これらのｘi ，ｙi を定数として次式
（１）で表わされる曲線を各エッジ候補点ごとに求め、 ρ＝ｘi cos θ＋ｙi sin θ （１） 該求められた曲線同士の各交点（ρj ，θj ）を求める
とともに、該各交点（ρj ，θj ）ごとの交わる前記曲
線の数をカウントし、このカウント値の大きい順に所定
の数の前記交点（ρj ，θj ）を抽出し、該抽出された
各交点（ρj ，θj ）に対応する、次式（２）で規定さ
れる前記直交座標系上の各基準候補線を求める基準候補
線検出手段と、 ρj ＝ｘ cosθj ＋ｙ sinθj （２） 前記検出された複数の基準候補線によって囲まれる領域
を照射野として定める照射野決定手段とを備えたことを
特徴とするものである。
上記本発明の放射線画像の照射野認識方法を実施するた
めの装置であって、照射野絞りを使用して撮影された照
射野を有する放射線画像の該照射野を、該放射線画像を
表す画像データに基づいて認識処理する照射野認識装置
において、前記照射野内に含まれる所定の点から前記放
射線画像の端部にそれぞれ向う複数の放射状の直線を設
定し、該設定された各放射状の直線ごとに、該各直線上
における前記照射野の輪郭と考えられる各エッジ候補点
を、該各直線に沿った前記画像データに基づいて求める
エッジ候補点検出手段と、前記エッジ候補点検出手段に
より検出された前記各エッジ候補点の座標を（ｘi ，ｙ
i ）としたとき、これらのｘi ，ｙi を定数として次式
（１）で表わされる曲線を各エッジ候補点ごとに求め、 ρ＝ｘi cos θ＋ｙi sin θ （１） 該求められた曲線同士の各交点（ρj ，θj ）を求める
とともに、該各交点（ρj ，θj ）ごとの交わる前記曲
線の数をカウントし、このカウント値の大きい順に所定
の数の前記交点（ρj ，θj ）を抽出し、該抽出された
各交点（ρj ，θj ）に対応する、次式（２）で規定さ
れる前記直交座標系上の各基準候補線を求める基準候補
線検出手段と、 ρj ＝ｘ cosθj ＋ｙ sinθj （２） 前記検出された複数の基準候補線によって囲まれる領域
を照射野として定める照射野決定手段とを備えたことを
特徴とするものである。
【００３３】ここで前記基準候補線検出手段により検出
された各基準候補線に関し、所定の評価関数を用いた評
価を行なう評価手段と、該評価手段による評価結果に基
づいて前記基準候補線のうちから又は前記基準候補線に
代えて、所定の候補線を求める候補線検出手段とをさら
に備え、前記照射野決定手段を、前記基準候補線に代え
て前記候補線によって囲まれる領域を前記照射野として
定めるものとした構成を採用することもできる。
された各基準候補線に関し、所定の評価関数を用いた評
価を行なう評価手段と、該評価手段による評価結果に基
づいて前記基準候補線のうちから又は前記基準候補線に
代えて、所定の候補線を求める候補線検出手段とをさら
に備え、前記照射野決定手段を、前記基準候補線に代え
て前記候補線によって囲まれる領域を前記照射野として
定めるものとした構成を採用することもできる。
【００３４】ここで、評価手段とは、基準候補線検出手
段により検出された各基準候補線に関し、照射野輪郭と
しての適格性若しくは適合性を評価するものをいい、候
補線検出手段はこの適格性等が高い候補線、すなわち照
射野輪郭を構成する可能性が基準候補線よりも少なくと
も同じか又は高い候補線を求めるものを意味する。
段により検出された各基準候補線に関し、照射野輪郭と
しての適格性若しくは適合性を評価するものをいい、候
補線検出手段はこの適格性等が高い候補線、すなわち照
射野輪郭を構成する可能性が基準候補線よりも少なくと
も同じか又は高い候補線を求めるものを意味する。
【００３５】そのような所定の評価関数を用いた評価を
行なう評価手段としては、第１に、前記式（２）で規定
された各基準候補線上に、前記照射野輪郭上の点として
予め設定された所定の代表点または評価手段自体により
設定された代表点が含まれるか否かを該各基準候補線ご
とに評価（判定）するものを適用することができ、この
場合、前記候補線検出手段は、前記基準候補線のうちか
ら前記代表点が含まれる基準候補線を求めるものであ
る。
行なう評価手段としては、第１に、前記式（２）で規定
された各基準候補線上に、前記照射野輪郭上の点として
予め設定された所定の代表点または評価手段自体により
設定された代表点が含まれるか否かを該各基準候補線ご
とに評価（判定）するものを適用することができ、この
場合、前記候補線検出手段は、前記基準候補線のうちか
ら前記代表点が含まれる基準候補線を求めるものであ
る。
【００３６】照射野輪郭上の点として予め設定された所
定の代表点とは、実際の画像における照射野輪郭上に存
在する蓋然性の高い点を意味し、例えば前述した画像デ
ータの差分が極めて大きい値を示す等の点が該当する。
定の代表点とは、実際の画像における照射野輪郭上に存
在する蓋然性の高い点を意味し、例えば前述した画像デ
ータの差分が極めて大きい値を示す等の点が該当する。
【００３７】所定の評価関数を用いた評価を行なう評価
手段としては、第２に、式（２）で規定された各基準候
補線を、該基準候補線を中心とする±ｍ画素、±α度の
範囲内でシフトおよび／または回転させて、｛（２ｍ＋
１）（２α＋１）−１｝個の変換候補線を設定する変換
候補線設定手段と、該変換候補線設定手段により設定さ
れた前記基準候補線および前記各変換候補線ごとに、各
候補線に関する微分値の、各候補線の延在方向について
の平均値を求める評価値算出手段とをさらに備え、前記
候補線検出手段が、前記基準候補線に代えて、前記微分
値の平均値が大きい順に所定数の候補線を求めるもので
あることを特徴とするものである。
手段としては、第２に、式（２）で規定された各基準候
補線を、該基準候補線を中心とする±ｍ画素、±α度の
範囲内でシフトおよび／または回転させて、｛（２ｍ＋
１）（２α＋１）−１｝個の変換候補線を設定する変換
候補線設定手段と、該変換候補線設定手段により設定さ
れた前記基準候補線および前記各変換候補線ごとに、各
候補線に関する微分値の、各候補線の延在方向について
の平均値を求める評価値算出手段とをさらに備え、前記
候補線検出手段が、前記基準候補線に代えて、前記微分
値の平均値が大きい順に所定数の候補線を求めるもので
あることを特徴とするものである。
【００３８】所定の評価関数を用いた評価を行なう評価
手段としては、第３に、式（２）で規定された各基準候
補線を、該基準候補線を中心とする±ｍ画素、±α度の
範囲内でシフトおよび／または回転させて、｛（２ｍ＋
１）（２α＋１）−１｝個の変換候補線を設定する変換
候補線設定手段と、該変換候補線設定手段により設定さ
れた前記基準候補線および前記各変換候補線ごとに、各
候補線に関する濃度勾配ベクトルの向きの、各候補線の
延在方向についてのエントロピーを求める評価値算出手
段とをさらに備え、前記候補線検出手段が、前記基準候
補線に代えて、前記エントロピーが小さい順に所定数の
候補線を求めるものであることを特徴とするものであ
る。
手段としては、第３に、式（２）で規定された各基準候
補線を、該基準候補線を中心とする±ｍ画素、±α度の
範囲内でシフトおよび／または回転させて、｛（２ｍ＋
１）（２α＋１）−１｝個の変換候補線を設定する変換
候補線設定手段と、該変換候補線設定手段により設定さ
れた前記基準候補線および前記各変換候補線ごとに、各
候補線に関する濃度勾配ベクトルの向きの、各候補線の
延在方向についてのエントロピーを求める評価値算出手
段とをさらに備え、前記候補線検出手段が、前記基準候
補線に代えて、前記エントロピーが小さい順に所定数の
候補線を求めるものであることを特徴とするものであ
る。
【００３９】本発明の放射線画像の黒化処理方法および
黒化処理装置は、上記本発明の照射野認識方法または認
識装置により得られた照射野の外側部分について黒化処
理を行うことを特徴とするものである。
黒化処理装置は、上記本発明の照射野認識方法または認
識装置により得られた照射野の外側部分について黒化処
理を行うことを特徴とするものである。
【００４０】すなわち本発明の放射線画像の黒化処理方
法は、照射野絞りを使用して撮影された照射野を有する
放射線画像の該照射野の外側部分について、略最高濃度
または略最低輝度となるように、前記放射線画像を表す
画像データを処理する放射線画像の黒化処理方法におい
て、前記照射野として、本発明の放射線画像の照射野認
識方法により得られた照射野を適用することを特徴とす
るものである。
法は、照射野絞りを使用して撮影された照射野を有する
放射線画像の該照射野の外側部分について、略最高濃度
または略最低輝度となるように、前記放射線画像を表す
画像データを処理する放射線画像の黒化処理方法におい
て、前記照射野として、本発明の放射線画像の照射野認
識方法により得られた照射野を適用することを特徴とす
るものである。
【００４１】また本発明の放射線画像の黒化処理装置
は、照射野絞りを使用して撮影された照射野を有する放
射線画像の該照射野の外側部分について、略最高濃度ま
たは略最低輝度となるように、前記放射線画像を表す画
像データを処理する放射線画像の黒化処理装置におい
て、前記照射野として、本発明の放射線画像の照射野認
識装置により得られた照射野を適用することを特徴とす
るものである。
は、照射野絞りを使用して撮影された照射野を有する放
射線画像の該照射野の外側部分について、略最高濃度ま
たは略最低輝度となるように、前記放射線画像を表す画
像データを処理する放射線画像の黒化処理装置におい
て、前記照射野として、本発明の放射線画像の照射野認
識装置により得られた照射野を適用することを特徴とす
るものである。
【００４２】
【発明の効果】本発明の放射線画像の照射野認識方法お
よび照射野認識装置は、まず画像中の所定の点について
設定した放射状の直線方向に沿ってエッジ候補点を検出
し、これらのエッジ候補点についてハフ（ Hough）変換
を利用してエッジ線となるべき所定数の基準候補線を求
め、これらの基準候補線で囲まれる領域を照射野領域と
して定める。
よび照射野認識装置は、まず画像中の所定の点について
設定した放射状の直線方向に沿ってエッジ候補点を検出
し、これらのエッジ候補点についてハフ（ Hough）変換
を利用してエッジ線となるべき所定数の基準候補線を求
め、これらの基準候補線で囲まれる領域を照射野領域と
して定める。
【００４３】このように本発明の放射線画像の照射野認
識方法・装置は、第１段階で画像面内の放射状の直線方
向に沿って多数のエッジ候補点を検出し、この得られた
多数のエッジ候補点に基づいて第２段階でハフ変換を利
用してこれら多数のエッジ候補点を結ぶ直線成分である
基準候補線を抽出するため、第１段階と同じ方法で検出
された多数のエッジ候補点に基づいて、予め設定された
形状の範囲内で照射野輪郭を検出していた従来の照射野
認識方法・装置に比べて、抽出し得る照射野輪郭の形状
の自由度を拡大することができ、実際の照射野輪郭によ
り一層適合した、すなわち精度を一層向上させて照射野
を認識することができる。
識方法・装置は、第１段階で画像面内の放射状の直線方
向に沿って多数のエッジ候補点を検出し、この得られた
多数のエッジ候補点に基づいて第２段階でハフ変換を利
用してこれら多数のエッジ候補点を結ぶ直線成分である
基準候補線を抽出するため、第１段階と同じ方法で検出
された多数のエッジ候補点に基づいて、予め設定された
形状の範囲内で照射野輪郭を検出していた従来の照射野
認識方法・装置に比べて、抽出し得る照射野輪郭の形状
の自由度を拡大することができ、実際の照射野輪郭によ
り一層適合した、すなわち精度を一層向上させて照射野
を認識することができる。
【００４４】また本発明の放射線画像の黒化処理方法・
装置は、上述した本発明の放射線画像の照射野認識方法
または認識装置により、精度を向上させて検出された照
射野輪郭を利用して、その照射野外領域について黒化処
理を行うことにより、従来の黒化処理方法・装置に比し
て、照射野外領域がより精度よく黒化処理された放射線
画像を得ることができる。
装置は、上述した本発明の放射線画像の照射野認識方法
または認識装置により、精度を向上させて検出された照
射野輪郭を利用して、その照射野外領域について黒化処
理を行うことにより、従来の黒化処理方法・装置に比し
て、照射野外領域がより精度よく黒化処理された放射線
画像を得ることができる。
【００４５】
【発明の実施の形態】以下、本発明の放射線画像の照射
野認識装置の具体的な実施の形態について図面を用いて
説明する。
野認識装置の具体的な実施の形態について図面を用いて
説明する。
【００４６】図１は本発明の照射野認識装置の一実施形
態の構成を示す図、図２は照射野絞りを用いた放射線画
像の撮影装置を示す図（（１）参照）および照射野絞り
の開口輪郭に対応する照射野輪郭が形成された蓄積性蛍
光体シートを示す図（（２）参照）、図３は図１に示し
た照射野認識装置による照射野認識の作用を説明するた
めの図である。
態の構成を示す図、図２は照射野絞りを用いた放射線画
像の撮影装置を示す図（（１）参照）および照射野絞り
の開口輪郭に対応する照射野輪郭が形成された蓄積性蛍
光体シートを示す図（（２）参照）、図３は図１に示し
た照射野認識装置による照射野認識の作用を説明するた
めの図である。
【００４７】図１に示した照射野認識装置により読み取
られる画像の取得方法についてまず説明する。すなわ
ち、図２（１）に示すように、Ｘ線源と被写体との間
に、矩形等の開口輪郭を有し、この開口の外側の部分
は、Ｘ線が被写体および蓄積性蛍光体シートにＸ線が到
達するのを防止する鉛板である照射野絞りが配され、こ
の状態でＸ線源から被写体にＸ線が照射される。この結
果、被写体を透過したＸ線が蓄積性蛍光体シートに照射
されるが、その際、図２（１）に示すように、照射野絞
りの開口輪郭より外側に対応する領域（照射野外領域）
Ｐout にはＸ線が全く照射されず、一方、照射野絞りの
開口輪郭より内側に対応する領域（照射野領域）Ｐinに
は被写体の透過Ｘ線情報が記録される。そして、照射野
絞りの開口輪郭に対応する部分は、開口輪郭と略同一形
状で、濃度が急峻に変化する複数のエッジ線からなる照
射野輪郭ＰＳが形成される。
られる画像の取得方法についてまず説明する。すなわ
ち、図２（１）に示すように、Ｘ線源と被写体との間
に、矩形等の開口輪郭を有し、この開口の外側の部分
は、Ｘ線が被写体および蓄積性蛍光体シートにＸ線が到
達するのを防止する鉛板である照射野絞りが配され、こ
の状態でＸ線源から被写体にＸ線が照射される。この結
果、被写体を透過したＸ線が蓄積性蛍光体シートに照射
されるが、その際、図２（１）に示すように、照射野絞
りの開口輪郭より外側に対応する領域（照射野外領域）
Ｐout にはＸ線が全く照射されず、一方、照射野絞りの
開口輪郭より内側に対応する領域（照射野領域）Ｐinに
は被写体の透過Ｘ線情報が記録される。そして、照射野
絞りの開口輪郭に対応する部分は、開口輪郭と略同一形
状で、濃度が急峻に変化する複数のエッジ線からなる照
射野輪郭ＰＳが形成される。
【００４８】このようにして蓄積性蛍光体シートに記録
された放射線画像（以下、原画像という）Ｐは図示しな
い公知の放射線画像情報読取装置によって読み取られ、
この原画像Ｐを表すデジタル画像データＳが出力され
る。図示の照射野認識装置100には、この画像データＳ
が入力される。
された放射線画像（以下、原画像という）Ｐは図示しな
い公知の放射線画像情報読取装置によって読み取られ、
この原画像Ｐを表すデジタル画像データＳが出力され
る。図示の照射野認識装置100には、この画像データＳ
が入力される。
【００４９】この照射野認識装置100 は、入力された画
像データＳに基づいて、上記照射野輪郭ＰＳ上の点と考
えられる多数のエッジ候補点Ｅを検出するエッジ候補点
検出手段10と、多数のエッジ候補点Ｅに基づいてハフ変
換を利用してそれらのエッジ候補点Ｅが構成するエッジ
線（基準候補線）Ｌを求める基準候補線検出手段20と、
得られた複数に基準候補線に基づいて照射野を決定する
照射野決定手段30とを備えた構成である。
像データＳに基づいて、上記照射野輪郭ＰＳ上の点と考
えられる多数のエッジ候補点Ｅを検出するエッジ候補点
検出手段10と、多数のエッジ候補点Ｅに基づいてハフ変
換を利用してそれらのエッジ候補点Ｅが構成するエッジ
線（基準候補線）Ｌを求める基準候補線検出手段20と、
得られた複数に基準候補線に基づいて照射野を決定する
照射野決定手段30とを備えた構成である。
【００５０】これらの構成の詳細を以下に説明する。
【００５１】エッジ候補点検出手段10は，入力された画
像データＳについて、図３（１）に示すように、そのデ
ータＳが表す画像Ｐの中心点Ｋから放射線画像の端部に
それぞれ向う等角度間隔（α＝３度） 120本の放射状の
直線Ｌ０を設定する。この放射状の直線の設定本数およ
び間隔、放射状の直線群の中心の設定位置等は適宜変更
することができる。
像データＳについて、図３（１）に示すように、そのデ
ータＳが表す画像Ｐの中心点Ｋから放射線画像の端部に
それぞれ向う等角度間隔（α＝３度） 120本の放射状の
直線Ｌ０を設定する。この放射状の直線の設定本数およ
び間隔、放射状の直線群の中心の設定位置等は適宜変更
することができる。
【００５２】さらにエッジ候補点検出手段10は、これら
の 120本の放射状の直線Ｌ０にそれぞれ沿った方向にお
いて互いに隣接する画素間の画像データを比較し、その
差分が最も大きい画素組をその線上で探索する。この探
索された画素組の間には、濃度差が大きいエッジ部に相
当するエッジ候補点が存在することを意味する。このよ
うにして各直線方向ごとにそれぞれ探索されたエッジ候
補点Ｅが検出され、合計 120個のエッジ候補点Ｅが得ら
れる。
の 120本の放射状の直線Ｌ０にそれぞれ沿った方向にお
いて互いに隣接する画素間の画像データを比較し、その
差分が最も大きい画素組をその線上で探索する。この探
索された画素組の間には、濃度差が大きいエッジ部に相
当するエッジ候補点が存在することを意味する。このよ
うにして各直線方向ごとにそれぞれ探索されたエッジ候
補点Ｅが検出され、合計 120個のエッジ候補点Ｅが得ら
れる。
【００５３】基準候補線検出手段20は、エッジ候補点検
出手段10により検出された各エッジ候補点Ｅの、図３
（１）に示すｘｙ座標系における座標を（ｘi ，ｙi ）
（ i＝１，２，…，120 ）としたとき、これらのｘi ，
ｙi を定数として次式（１）で表わされるサインカーブ
ＬＳを各エッジ候補点ごとに求める。
出手段10により検出された各エッジ候補点Ｅの、図３
（１）に示すｘｙ座標系における座標を（ｘi ，ｙi ）
（ i＝１，２，…，120 ）としたとき、これらのｘi ，
ｙi を定数として次式（１）で表わされるサインカーブ
ＬＳを各エッジ候補点ごとに求める。
【００５４】 ρ＝ｘi cos θ＋ｙi sin θ （１） この式（１）は、図３（２）に示すように、中心座標を
（ｘi ，ｙi ）に固定してこの座標（ｘi ，ｙi ）を通
過する直線の式を表し、ρは座標（ｘi ，ｙi）を通過
する各直線Ｌに向けてｘｙ座標系の原点Ｏから下ろした
垂線の長さ、θはこの垂線とｘ軸とのなす角度を表して
いる。なお、このθ 度間隔とする。
（ｘi ，ｙi ）に固定してこの座標（ｘi ，ｙi ）を通
過する直線の式を表し、ρは座標（ｘi ，ｙi）を通過
する各直線Ｌに向けてｘｙ座標系の原点Ｏから下ろした
垂線の長さ、θはこの垂線とｘ軸とのなす角度を表して
いる。なお、このθ 度間隔とする。
【００５５】ここで座標（ｘi ，ｙi ）を通過する各直
線Ｌは式（１）においてρ、θを変化させたものであ
り、θρ空間（ハフ空間）では、座標（ｘi ，ｙi ）を
通過するこれら直線Ｌは図３（３）に示す１本のサイン
カーブとして表現される。
線Ｌは式（１）においてρ、θを変化させたものであ
り、θρ空間（ハフ空間）では、座標（ｘi ，ｙi ）を
通過するこれら直線Ｌは図３（３）に示す１本のサイン
カーブとして表現される。
【００５６】各エッジ候補点についてそれぞれ同様の操
作を行なうと、ハフ空間に 120本のサインカーブＬＳi
（ i＝１，２，…，120 ）が表される。
作を行なうと、ハフ空間に 120本のサインカーブＬＳi
（ i＝１，２，…，120 ）が表される。
【００５７】ここで基準候補線検出手段20は、 120本の
サインカーブＬＳi 同士の各交点位置（ρj ，θj ）を
求めるとともに、該各交点位置（ρj ，θj ）ごとの交
わる曲線の数をカウントする。なお交点位置（ρj ，θ
j ）において２つのサインカーブが交差した場合には、
両サインカーブとも座標（ρj ，θj ）を通過している
ことを意味し、これはすなわち互いに異なる２つのエッ
ジ候補点Ｅ１，Ｅ２が、（ρj ，θj ）に対応するｘｙ
座標系における座標（ｘ，ｙ）を通過する下記式（２）
で表される同一直線上に存在することを意味する。
サインカーブＬＳi 同士の各交点位置（ρj ，θj ）を
求めるとともに、該各交点位置（ρj ，θj ）ごとの交
わる曲線の数をカウントする。なお交点位置（ρj ，θ
j ）において２つのサインカーブが交差した場合には、
両サインカーブとも座標（ρj ，θj ）を通過している
ことを意味し、これはすなわち互いに異なる２つのエッ
ジ候補点Ｅ１，Ｅ２が、（ρj ，θj ）に対応するｘｙ
座標系における座標（ｘ，ｙ）を通過する下記式（２）
で表される同一直線上に存在することを意味する。
【００５８】 ρj ＝ｘcos θj ＋ｙsin θj （２） したがって交差するサインカーブの数（カウント値）が
多い交点こそ、多くのエッジ候補点Ｅをその線上に含む
照射野輪郭ＰＳであることを意味する。
多い交点こそ、多くのエッジ候補点Ｅをその線上に含む
照射野輪郭ＰＳであることを意味する。
【００５９】基準候補線検出手段20は、全てのエッジ候
補点Ｅについて上記カウントを行ない、予め設定された
個数に到達するまでカウント値が上位のものから順に交
点位置（ρj ，θj ）を抽出し、それらが表すｘｙ座標
系における直線を求める。この求められた所定数の直線
が基準候補線ＬＬである。
補点Ｅについて上記カウントを行ない、予め設定された
個数に到達するまでカウント値が上位のものから順に交
点位置（ρj ，θj ）を抽出し、それらが表すｘｙ座標
系における直線を求める。この求められた所定数の直線
が基準候補線ＬＬである。
【００６０】照射野決定手段30は、得られた複数に基準
候補線ＬＬに基づいて照射野輪郭ＰＳを決定し、この照
射野輪郭ＰＳに相当する画像データを外部の画像処理装
置等に出力する。
候補線ＬＬに基づいて照射野輪郭ＰＳを決定し、この照
射野輪郭ＰＳに相当する画像データを外部の画像処理装
置等に出力する。
【００６１】次に本実施形態の放射線画像の照射野認識
装置100 の作用について説明する。
装置100 の作用について説明する。
【００６２】図示の照射野認識装置100 に、照射野輪郭
ＰＳを有する画像Ｐを表す画像データＳが入力される
と、まずエッジ候補点検出手段10が、入力された画像デ
ータＳについての中心点Ｋを設定し、この中心点Ｋから
等角度間隔の 120本の放射状の直線Ｌ０を設定する。
ＰＳを有する画像Ｐを表す画像データＳが入力される
と、まずエッジ候補点検出手段10が、入力された画像デ
ータＳについての中心点Ｋを設定し、この中心点Ｋから
等角度間隔の 120本の放射状の直線Ｌ０を設定する。
【００６３】次いでエッジ候補点検出手段10は、放射状
の各直線ごとにそれらの直線にそれぞれ沿った方向にお
いて互いに隣接する画素間の画像データを比較し、その
差分が最も大きいエッジ候補点Ｅを探索する。放射線の
全ての方向においてエッジ候補点Ｅが検出されると、こ
れらの各エッジ候補点Ｅの位置を表す座標が基準候補線
検出手段20に入力される。このときの座標系は放射線画
像面内におけるｘｙ座標系によるものである。
の各直線ごとにそれらの直線にそれぞれ沿った方向にお
いて互いに隣接する画素間の画像データを比較し、その
差分が最も大きいエッジ候補点Ｅを探索する。放射線の
全ての方向においてエッジ候補点Ｅが検出されると、こ
れらの各エッジ候補点Ｅの位置を表す座標が基準候補線
検出手段20に入力される。このときの座標系は放射線画
像面内におけるｘｙ座標系によるものである。
【００６４】基準候補線検出手段20は、入力された 120
点のエッジ候補点に基づいて、基準候補線検出手段20に
予め設定された本数の基準候補線ＬＬを求める。図４
（１）はこの所定の本数が５本の場合を示すもので、得
られた基準候補線ＬＬ（ LL1〜LL5 ）は、その一部が照
射野輪郭ＰＳの辺を構成し得る直線となる。
点のエッジ候補点に基づいて、基準候補線検出手段20に
予め設定された本数の基準候補線ＬＬを求める。図４
（１）はこの所定の本数が５本の場合を示すもので、得
られた基準候補線ＬＬ（ LL1〜LL5 ）は、その一部が照
射野輪郭ＰＳの辺を構成し得る直線となる。
【００６５】そしてこれらの求められた複数の基準候補
線ＬＬは照射野決定手段30に入力され、照射野決定手段
30はこれら入力された複数の基準候補線ＬＬに囲まれる
領域を照射野領域Ｐinと決定し、この照射野領域Ｐinの
輪郭を照射野輪郭ＰＳと定め、照射野輪郭ＰＳを構成す
る画素を示す位置のデータ（図４（２））を外部の画像
処理装置等に出力する。
線ＬＬは照射野決定手段30に入力され、照射野決定手段
30はこれら入力された複数の基準候補線ＬＬに囲まれる
領域を照射野領域Ｐinと決定し、この照射野領域Ｐinの
輪郭を照射野輪郭ＰＳと定め、照射野輪郭ＰＳを構成す
る画素を示す位置のデータ（図４（２））を外部の画像
処理装置等に出力する。
【００６６】外部の画像処理装置等はこのように出力さ
れた照射野輪郭ＰＳを構成する画素の位置データに基づ
いて、他の経路から入力された放射線画像データＳのう
ち、照射野輪郭ＰＳよりも外側領域である照射野外領域
Ｐout に対応する画像データを一律に最高濃度の画像デ
ータに置換する黒化処理等を施す。そして照射野外領域
が黒化処理された画像は観察読影に供される。
れた照射野輪郭ＰＳを構成する画素の位置データに基づ
いて、他の経路から入力された放射線画像データＳのう
ち、照射野輪郭ＰＳよりも外側領域である照射野外領域
Ｐout に対応する画像データを一律に最高濃度の画像デ
ータに置換する黒化処理等を施す。そして照射野外領域
が黒化処理された画像は観察読影に供される。
【００６７】このように本実施形態の放射線画像の照射
野認識装置によれば、画像中の所定の点について設定し
た放射状の直線方向に沿って検出したエッジ候補点につ
いて、さらにハフ（ Hough）変換を適用して所定数の基
準候補線を求め、これらの基準候補線で囲まれる領域を
照射野領域として定めることによって、従来の照射野認
識方法または装置よりも照射野輪郭の認識精度を向上さ
せることができる。
野認識装置によれば、画像中の所定の点について設定し
た放射状の直線方向に沿って検出したエッジ候補点につ
いて、さらにハフ（ Hough）変換を適用して所定数の基
準候補線を求め、これらの基準候補線で囲まれる領域を
照射野領域として定めることによって、従来の照射野認
識方法または装置よりも照射野輪郭の認識精度を向上さ
せることができる。
【００６８】図５は本発明の放射線画像の照射野認識装
置の第２の実施形態を示す図である。
置の第２の実施形態を示す図である。
【００６９】図示の照射野認識装置は、図１に示した実
施形態の照射野認識装置100 において、基準候補線検出
手段20によって、図１に示した実施形態よりも多くの基
準候補線を求めておき、このように求められたより多く
の基準候補線に関し、さらに所定の評価関数を用いた評
価を行なう評価手段40と、この評価手段40による評価結
果に基づいて基準候補線のうちから又は基準候補線に代
えて、所定の照射野輪郭である精度がより一層高い候補
線ＬＬ′を求める候補線検出手段50とを備え、照射野決
定手段として、基準候補線ＬＬによって囲まれる領域を
照射野として定めるものに代えて候補線検出手段50によ
り求められた候補線ＬＬ′によって囲まれる領域を照射
野として定めるものとした構成を採用したものである。
施形態の照射野認識装置100 において、基準候補線検出
手段20によって、図１に示した実施形態よりも多くの基
準候補線を求めておき、このように求められたより多く
の基準候補線に関し、さらに所定の評価関数を用いた評
価を行なう評価手段40と、この評価手段40による評価結
果に基づいて基準候補線のうちから又は基準候補線に代
えて、所定の照射野輪郭である精度がより一層高い候補
線ＬＬ′を求める候補線検出手段50とを備え、照射野決
定手段として、基準候補線ＬＬによって囲まれる領域を
照射野として定めるものに代えて候補線検出手段50によ
り求められた候補線ＬＬ′によって囲まれる領域を照射
野として定めるものとした構成を採用したものである。
【００７０】このような構成を採用した本実施形態の照
射野認識装置によれば、図１に示した実施形態の照射野
認識装置よりもさらに照射野輪郭としての適格性が高い
候補線を求めることができ、したがって照射野輪郭の認
識精度を向上させることができる。
射野認識装置によれば、図１に示した実施形態の照射野
認識装置よりもさらに照射野輪郭としての適格性が高い
候補線を求めることができ、したがって照射野輪郭の認
識精度を向上させることができる。
【００７１】ここで評価手段40としては、例えば、放射
線画像を図６に示すように画像中心点等の所定の点（所
定の点がエッジ候補点検出手段10により設定された場合
は、当該設定された点であってもよい）周りに角度間隔
45度ごとの８つの領域Ｐ１，Ｐ２，…，Ｐ８に分割し、
この分割された各領域の中に含まれる、エッジ候補点検
出手段10により設定された各15個の各エッジ候補点（図
中・印で示す）のうち差分が最も大きい値を示すエッジ
候補点（図中×印で示す）を代表点として設定し、各基
準候補線上にこれらの代表点が含まれるか否かを各基準
候補線ごとに判定するものなどである。
線画像を図６に示すように画像中心点等の所定の点（所
定の点がエッジ候補点検出手段10により設定された場合
は、当該設定された点であってもよい）周りに角度間隔
45度ごとの８つの領域Ｐ１，Ｐ２，…，Ｐ８に分割し、
この分割された各領域の中に含まれる、エッジ候補点検
出手段10により設定された各15個の各エッジ候補点（図
中・印で示す）のうち差分が最も大きい値を示すエッジ
候補点（図中×印で示す）を代表点として設定し、各基
準候補線上にこれらの代表点が含まれるか否かを各基準
候補線ごとに判定するものなどである。
【００７２】そして、この場合、候補線検出手段50とし
ては、基準候補線のうちから代表点が含まれると判定さ
れた基準候補線を候補線として設定するものである。
ては、基準候補線のうちから代表点が含まれると判定さ
れた基準候補線を候補線として設定するものである。
【００７３】なお代表点は評価手段40により設定される
場合の他、予め他の構成手段等により設定されるもので
あってもよい。
場合の他、予め他の構成手段等により設定されるもので
あってもよい。
【００７４】図７は、評価手段および候補線検出手段を
他の態様とした実施形態を示す図である。
他の態様とした実施形態を示す図である。
【００７５】すなわち本実施形態の照射野認識装置にお
ける評価手段40′は、基準候補線検出手段20によって求
められた各基準候補線を、図８に示すように基準候補線
を中心とする±ｍ画素、±α度の範囲内でシフトおよび
／または回転させて、｛（２ｍ＋１）（２α＋１）−
１｝個の変換候補線を設定する候補線作成手段41と、元
の基準候補線およびこれら各基準候補線ごとにそれぞれ
設定された各変換候補線ごとに、 図９（１）に示すように、各候補線（基準候補線およ
び変換候補線）に関する微分値の、各候補線の延在方向
についての平均値を求める関数、または、 図９（２）に示すように、各候補線（基準候補線およ
び変換候補線）に関する濃度勾配ベクトルの向きの、各
候補線の延在方向についてのエントロピーを求める関
数、の演算を行なう評価値算出手段42とを備え、候補線
検出手段50′は、の場合は微分値の平均値が大きい順
に所定数の候補線を求めるものであり、の場合はエン
トロピーの値が小さい順に所定数の候補線を求めるもの
である。
ける評価手段40′は、基準候補線検出手段20によって求
められた各基準候補線を、図８に示すように基準候補線
を中心とする±ｍ画素、±α度の範囲内でシフトおよび
／または回転させて、｛（２ｍ＋１）（２α＋１）−
１｝個の変換候補線を設定する候補線作成手段41と、元
の基準候補線およびこれら各基準候補線ごとにそれぞれ
設定された各変換候補線ごとに、 図９（１）に示すように、各候補線（基準候補線およ
び変換候補線）に関する微分値の、各候補線の延在方向
についての平均値を求める関数、または、 図９（２）に示すように、各候補線（基準候補線およ
び変換候補線）に関する濃度勾配ベクトルの向きの、各
候補線の延在方向についてのエントロピーを求める関
数、の演算を行なう評価値算出手段42とを備え、候補線
検出手段50′は、の場合は微分値の平均値が大きい順
に所定数の候補線を求めるものであり、の場合はエン
トロピーの値が小さい順に所定数の候補線を求めるもの
である。
【００７６】上記各候補線に関する微分値とは、各候補
線を挟んで相対向する２つの画素間の画像データの差分
を意味する。この差分値が大きい場合は、この２つの画
素間で濃度の差が大きいことを意味し、画像のエッジが
存在することを意味する。そして図９（１）に示すよう
に、このような候補線を挟んで相対向する２つの画素の
組（図中破線で示す）は、候補線の延びる方向にそれぞ
れ存在するため、上記微分値を、候補線の延びる方向に
それぞれ存在する画素の組について加算し、その和を加
算した組数で除することにより上記平均値を求める。こ
の平均値が大きい程、当該候補線が放射線画像の現実の
エッジである可能性が高い。
線を挟んで相対向する２つの画素間の画像データの差分
を意味する。この差分値が大きい場合は、この２つの画
素間で濃度の差が大きいことを意味し、画像のエッジが
存在することを意味する。そして図９（１）に示すよう
に、このような候補線を挟んで相対向する２つの画素の
組（図中破線で示す）は、候補線の延びる方向にそれぞ
れ存在するため、上記微分値を、候補線の延びる方向に
それぞれ存在する画素の組について加算し、その和を加
算した組数で除することにより上記平均値を求める。こ
の平均値が大きい程、当該候補線が放射線画像の現実の
エッジである可能性が高い。
【００７７】したがって候補線検出手段50′が評価値算
出手段42から出力された微分値の平均値が大きい順に求
めた所定数の候補線は、図１に示した実施形態の照射野
認識装置よりもさらに照射野輪郭としての適格性が高い
候補線となり、照射野輪郭の認識精度を向上させること
ができる。
出手段42から出力された微分値の平均値が大きい順に求
めた所定数の候補線は、図１に示した実施形態の照射野
認識装置よりもさらに照射野輪郭としての適格性が高い
候補線となり、照射野輪郭の認識精度を向上させること
ができる。
【００７８】一方、上記各候補線に関する濃度勾配ベク
トルとは、各候補線の一方の側の各画素から濃度勾配
（画像データの差分値）が最も大きい方向を向くベクト
ルを意味する。そして、この濃度勾配ベクトルの向いた
方向を示す指標値（例えば候補線に直交する方向となす
角度βの正弦値（ sinβ））を求め、図９（２）に示す
ように、候補線の延びる方向にそれぞれ存在する画素か
らの各濃度勾配ベクトルの向きについて、上記指標値ご
とのヒストグラムを作成し、このヒストグラムの確率密
度Piに基づいてエントロピー（−Σ{(Pi)log Pi}）を求
める。このエントロピーが小さい程、濃度勾配ベクトル
の向きが一定の方向に揃っていることを意味し、そのよ
うな候補線が放射線画像の現実のエッジである可能性が
高い。
トルとは、各候補線の一方の側の各画素から濃度勾配
（画像データの差分値）が最も大きい方向を向くベクト
ルを意味する。そして、この濃度勾配ベクトルの向いた
方向を示す指標値（例えば候補線に直交する方向となす
角度βの正弦値（ sinβ））を求め、図９（２）に示す
ように、候補線の延びる方向にそれぞれ存在する画素か
らの各濃度勾配ベクトルの向きについて、上記指標値ご
とのヒストグラムを作成し、このヒストグラムの確率密
度Piに基づいてエントロピー（−Σ{(Pi)log Pi}）を求
める。このエントロピーが小さい程、濃度勾配ベクトル
の向きが一定の方向に揃っていることを意味し、そのよ
うな候補線が放射線画像の現実のエッジである可能性が
高い。
【００７９】したがって候補線検出手段50′が評価値算
出手段42から出力されたエントロピーが小さい順に求め
た所定数の候補線は、図１に示した実施形態の照射野認
識装置よりもさらに照射野輪郭としての適格性が高い候
補線となり、照射野輪郭の認識精度を向上させることが
できる。
出手段42から出力されたエントロピーが小さい順に求め
た所定数の候補線は、図１に示した実施形態の照射野認
識装置よりもさらに照射野輪郭としての適格性が高い候
補線となり、照射野輪郭の認識精度を向上させることが
できる。
【００８０】なお、本発明の放射線画像の照射野認識方
法および装置は、必ずしもこの実施形態に限定されるも
のではなく、評価手段も上記の態様に限るものではな
い。
法および装置は、必ずしもこの実施形態に限定されるも
のではなく、評価手段も上記の態様に限るものではな
い。
【００８１】さらに、入力される画像データが原画像デ
ータＳではなく、他の何らかの目的のために、既に縮小
された画像のデータ（例えば、蓄積性蛍光体シートに蓄
積記録された放射線画像情報のダイナミックレンジ等の
概要を把握するために、いわゆる本読みに先だって行な
われる、本読みよりも粗い画素密度で画像情報を読み取
る先読みによって読み取られた画像データ等）Ｓ′であ
ってもよい。
ータＳではなく、他の何らかの目的のために、既に縮小
された画像のデータ（例えば、蓄積性蛍光体シートに蓄
積記録された放射線画像情報のダイナミックレンジ等の
概要を把握するために、いわゆる本読みに先だって行な
われる、本読みよりも粗い画素密度で画像情報を読み取
る先読みによって読み取られた画像データ等）Ｓ′であ
ってもよい。
【００８２】図10は本発明の放射線画像の黒化処理装置
の一実施形態の構成を示す図である。
の一実施形態の構成を示す図である。
【００８３】図示の黒化処理装置200 は、図７に示した
実施形態の照射野認識装置により出力された照射野輪郭
ＰＳを構成する画素を示す位置のデータ（図４（２））
に基づいて、他の経路から入力された放射線画像データ
Ｓのうち、照射野輪郭ＰＳよりも外側領域である照射野
外領域Ｐout に対応する画像データを一律に最高濃度の
画像データに置換する黒化処理等を施すものである。
実施形態の照射野認識装置により出力された照射野輪郭
ＰＳを構成する画素を示す位置のデータ（図４（２））
に基づいて、他の経路から入力された放射線画像データ
Ｓのうち、照射野輪郭ＰＳよりも外側領域である照射野
外領域Ｐout に対応する画像データを一律に最高濃度の
画像データに置換する黒化処理等を施すものである。
【００８４】この実施形態の黒化処理装置200 によれ
ば、従来の照射野認識方法または装置よりも精度の高い
照射野輪郭に基づいて黒化処理がなされるため、従来の
黒化処理装置よりも照射野外領域がより精度よく黒化処
理された放射線画像を得ることができる。
ば、従来の照射野認識方法または装置よりも精度の高い
照射野輪郭に基づいて黒化処理がなされるため、従来の
黒化処理装置よりも照射野外領域がより精度よく黒化処
理された放射線画像を得ることができる。
【００８５】なお本発明の放射線画像の黒化処理装置
は、図10に示す実施形態の構成に限るものではなく、図
１または図５に示す本発明の照射野認識装置の各実施形
態から、または本発明の照射野認識装置のその他の態様
のものから出力された照射野輪郭に基づいて黒化処理を
なすものであれば、如何なる形態のものであってもよ
い。
は、図10に示す実施形態の構成に限るものではなく、図
１または図５に示す本発明の照射野認識装置の各実施形
態から、または本発明の照射野認識装置のその他の態様
のものから出力された照射野輪郭に基づいて黒化処理を
なすものであれば、如何なる形態のものであってもよ
い。
【図１】本発明の放射線画像の照射野認識装置の一実施
形態の構成を示す図
形態の構成を示す図
【図２】図１に示した照射野認識装置によって照射野領
域の検出処理の対象となる、照射野絞りを用いて撮影さ
れた蓄積性蛍光体シートを示す図
域の検出処理の対象となる、照射野絞りを用いて撮影さ
れた蓄積性蛍光体シートを示す図
【図３】（１）エッジ候補点検出手段10によるエッジ候
補点の検出作用を示す図、（２），（３）基準候補線検
出手段20による基準候補線の検出作用を示す図
補点の検出作用を示す図、（２），（３）基準候補線検
出手段20による基準候補線の検出作用を示す図
【図４】（１）基準候補線検出手段20により検出された
基準候補線ＬＬを示す図、（２）照射野決定手段30によ
り決定された照射野輪郭ＰＳを示す図
基準候補線ＬＬを示す図、（２）照射野決定手段30によ
り決定された照射野輪郭ＰＳを示す図
【図５】本発明の放射線画像の照射野認識装置の第２の
実施形態の構成を示す図
実施形態の構成を示す図
【図６】図５に示した実施形態において用いられる代表
点を示す図
点を示す図
【図７】本発明の放射線画像の照射野認識装置の第３の
実施形態の構成を示す図
実施形態の構成を示す図
【図８】複数の候補線の設定作用を説明するための図
【図９】評価値の設定作用を説明するための図
【図１０】本発明の放射線画像の黒化処理装置の一実施
形態の構成を示す図
形態の構成を示す図
10 エッジ候補点検出手段 20 基準候補線検出手段 30 照射野決定手段 40 評価手段 41 変換候補線設定手段 42 評価値算出手段 50 候補線設定手段 100 放射線画像の照射野認識装置 200 放射線画像の黒化処理装置
Claims (12)
- 【請求項１】 照射野絞りを使用して撮影された照射野
を有する放射線画像の該照射野を、該放射線画像を表す
画像データに基づいて認識処理する照射野認識方法にお
いて、 前記照射野内に含まれる所定の点から前記放射線画像の
端部にそれぞれ向う複数の放射状の直線を設定し、 該設定された各放射状の直線ごとに、該各直線上におけ
る前記照射野の輪郭と考えられる各エッジ候補点を、該
各直線に沿った前記画像データに基づいて求め、 前記放射線画像面内に直交座標系における前記各エッジ
候補点の座標を（ｘi，ｙi ）としたとき、これらのｘi
，ｙi を定数として次式（１）で表わされる曲線を各
エッジ候補点ごとに求め、 ρ＝ｘi cos θ＋ｙi sin θ （１） 該求められた曲線同士の各交点（ρj ，θj ）を求める
とともに、該各交点（ρj ，θj ）ごとの交わる前記曲
線の数をカウントし、 このカウント値の大きい順に所定の数の前記交点（ρj
，θj ）を抽出し、 該抽出された各交点（ρj ，θj ）に対応する、次式
（２）で規定される前記直交座標系上の各基準候補線を
求め、 ρj ＝ｘ cosθj ＋ｙ sinθj （２） 前記式（２）で規定される複数の基準候補線によって囲
まれる領域を前記照射野として定めることを特徴とする
放射線画像の照射野認識方法。 - 【請求項２】 前記式（２）で規定された各基準候補線
に関し、所定の評価関数を用いた評価を行ない、 この評価結果に基づいて前記基準候補線のうちから又は
前記基準候補線に代えて、所定の候補線を求め、 前記基準候補線に代えて該候補線によって囲まれる領域
を前記照射野として定めることを特徴とする請求項１記
載の放射線画像の照射野認識方法。 - 【請求項３】 前記所定の評価関数を用いた評価が、 前記式（２）で規定された各基準候補線上に、前記照射
野輪郭上の点として予め設定された所定の代表点が含ま
れるか否かを該各基準候補線ごとに判定することであ
り、 前記所定の候補線は、前記基準候補線のうち、前記所定
の代表点を含む基準候補線であることを特徴とする請求
項２記載の放射線画像の照射野認識方法。 - 【請求項４】 前記所定の評価関数を用いた評価が、 前記式（２）で規定された各基準候補線を、該基準候補
線を中心とする±ｍ画素、±α度の範囲内でシフトおよ
び／または回転させて、｛（２ｍ＋１）（２α＋１）−
１｝個の変換候補線を設定し、前記基準候補線および前
記各変換候補線ごとに、各候補線に関する微分値の、各
候補線の延在方向についての平均値を求めることであ
り、 前記所定の候補線は、前記基準候補線に代えて、前記微
分値の平均値が大きい順に選択された所定数の基準候補
線および変換候補線であることを特徴とする請求項２記
載の放射線画像の照射野認識方法。 - 【請求項５】 前記所定の評価関数を用いた評価が、 前記式（２）で規定された各基準候補線を、該基準候補
線を中心とする±ｍ画素、±α度の範囲内でシフトおよ
び／または回転させて、｛（２ｍ＋１）（２α＋１）−
１｝個の変換候補線を設定し、前記基準候補線および前
記各変換候補線ごとに、各候補線に関する濃度勾配ベク
トルの向きの、各候補線の延在方向についてのエントロ
ピーを求めることであり、 前記所定の候補線は、前記基準候補線に代えて、前記エ
ントロピーが小さい順に選択された所定数の基準候補線
および変換候補線であることを特徴とする請求項２記載
の放射線画像の照射野認識方法。 - 【請求項６】 照射野絞りを使用して撮影された照射野
を有する放射線画像の該照射野を、該放射線画像を表す
画像データに基づいて認識処理する照射野認識装置にお
いて、 前記照射野内に含まれる所定の点から前記放射線画像の
端部にそれぞれ向う複数の放射状の直線を設定し、該設
定された各放射状の直線ごとに、該各直線上における前
記照射野の輪郭と考えられる各エッジ候補点を、該各直
線に沿った前記画像データに基づいて求めるエッジ候補
点検出手段と、 前記エッジ候補点検出手段により検出された前記各エッ
ジ候補点の座標を（ｘi ，ｙi ）としたとき、これらの
ｘi ，ｙi を定数として次式（１）で表わされる曲線を
各エッジ候補点ごとに求め、 ρ＝ｘi cos θ＋ｙi sin θ （１） 該求められた曲線同士の各交点（ρj ，θj ）を求める
とともに、該各交点（ρj ，θj ）ごとの交わる前記曲
線の数をカウントし、 このカウント値の大きい順に所定の数の前記交点（ρj
，θj ）を抽出し、 該抽出された各交点（ρj ，θj ）に対応する、次式
（２）で規定される前記直交座標系上の各基準候補線を
求める基準候補線検出手段と、 ρj ＝ｘ cosθj ＋ｙ sinθj （２） 前記検出された複数の基準候補線によって囲まれる領域
を照射野として定める照射野決定手段とを備えたことを
特徴とする放射線画像の照射野認識装置。 - 【請求項７】 前記基準候補線検出手段により検出され
た各基準候補線に関し、所定の評価関数を用いた評価を
行なう評価手段と、該評価手段による評価結果に基づい
て前記基準候補線のうちから又は前記基準候補線に代え
て、所定の候補線を求める候補線検出手段とをさらに備
え、 前記照射野決定手段は、前記基準候補線に代えて前記候
補線によって囲まれる領域を前記照射野として定めるも
のとしたことを特徴とする請求項６記載の放射線画像の
照射野認識装置。 - 【請求項８】 前記評価手段が、前記式（２）で規定さ
れた各基準候補線上に、前記照射野輪郭上の点として予
め設定された所定の代表点が含まれるか否かを該各基準
候補線ごとに評価するものであり、 前記候補線検出手段が、前記基準候補線のうちから前記
代表点が含まれる基準候補線を求めるものであることを
特徴とする請求項７記載の放射線画像の照射野認識装
置。 - 【請求項９】 前記評価手段が、前記式（２）で規定さ
れた各基準候補線を、該基準候補線を中心とする±ｍ画
素、±α度の範囲内でシフトおよび／または回転させ
て、｛（２ｍ＋１）（２α＋１）−１｝個の変換候補線
を設定する変換候補線設定手段と、該変換候補線設定手
段により設定された前記基準候補線および前記各変換候
補線ごとに、各候補線に関する微分値の、各候補線の延
在方向についての平均値を求める評価値算出手段とをさ
らに備え、 前記候補線検出手段が、前記基準候補線に代えて、前記
微分値の平均値が大きい順に所定数の候補線を求めるも
のであることを特徴とする請求項７記載の放射線画像の
照射野認識装置。 - 【請求項１０】 前記評価手段が、前記式（２）で規定
された各基準候補線を、該基準候補線を中心とする±ｍ
画素、±α度の範囲内でシフトおよび／または回転させ
て、｛（２ｍ＋１）（２α＋１）−１｝個の変換候補線
を設定する変換候補線設定手段と、該変換候補線設定手
段により設定された前記基準候補線および前記各変換候
補線ごとに、各候補線に関する濃度勾配ベクトルの向き
の、各候補線の延在方向についてのエントロピーを求め
る評価値算出手段とをさらに備え、 前記候補線検出手段が、前記基準候補線に代えて、前記
エントロピーが小さい順に所定数の候補線を求めるもの
であることを特徴とする請求項７記載の放射線画像の照
射野認識装置。 - 【請求項１１】 照射野絞りを使用して撮影された照射
野を有する放射線画像の該照射野の外側部分について、
略最高濃度または略最低輝度となるように、前記放射線
画像を表す画像データを処理する放射線画像の黒化処理
方法において、 前記照射野として、請求項１から５のうちいずれか１項
に記載の放射線画像の照射野認識方法により得られた照
射野を適用することを特徴とする放射線画像の黒化処理
方法。 - 【請求項１２】 照射野絞りを使用して撮影された照射
野を有する放射線画像の該照射野の外側部分について、
略最高濃度または略最低輝度となるように、前記放射線
画像を表す画像データを処理する放射線画像の黒化処理
装置において、 前記照射野として、放射線画像のうち、請求項６から１
０のうちいずれか１項に記載の放射線画像の照射野認識
装置により得られた照射野を適用することを特徴とする
放射線画像の黒化処理装置。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP07969697A JP3923131B2 (ja) | 1997-03-31 | 1997-03-31 | 放射線画像の照射野認識方法および装置、並びに黒化処理方法および装置 |
| US09/050,353 US6356651B2 (en) | 1997-03-31 | 1998-03-31 | Method and apparatus for recognizing irradiation fields on radiation images |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP07969697A JP3923131B2 (ja) | 1997-03-31 | 1997-03-31 | 放射線画像の照射野認識方法および装置、並びに黒化処理方法および装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH10275213A true JPH10275213A (ja) | 1998-10-13 |
| JP3923131B2 JP3923131B2 (ja) | 2007-05-30 |
Family
ID=13697382
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP07969697A Expired - Fee Related JP3923131B2 (ja) | 1997-03-31 | 1997-03-31 | 放射線画像の照射野認識方法および装置、並びに黒化処理方法および装置 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US6356651B2 (ja) |
| JP (1) | JP3923131B2 (ja) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003265462A (ja) * | 2002-03-19 | 2003-09-24 | Hitachi Ltd | 関心領域抽出方法及び画像処理サーバ |
| JP2006181362A (ja) * | 2004-12-24 | 2006-07-13 | General Electric Co <Ge> | ディジタル画像放射線撮像においてコリメーション・エッジを検出するシステム、方法及び装置 |
| JP2006293522A (ja) * | 2005-04-07 | 2006-10-26 | Fuji Photo Film Co Ltd | 直線検出装置、直線検出方法およびそのプログラム |
| US7796584B2 (en) | 2003-09-30 | 2010-09-14 | Nec Corporation | Method for connection between communication networks of different types and gateway apparatus |
| US7991212B2 (en) | 2008-01-09 | 2011-08-02 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, and computer program |
| EP2372395A2 (en) | 2010-03-31 | 2011-10-05 | Fujifilm Corporation | Radiographic imaging system |
| US8238633B2 (en) | 2010-09-29 | 2012-08-07 | Fujifilm Corporation | Image processing apparatus, radiographic image system, and image processing method and program |
| JP2015123157A (ja) * | 2013-12-26 | 2015-07-06 | コニカミノルタ株式会社 | 画像処理装置及び照射野認識方法 |
| JP2016517100A (ja) * | 2014-02-22 | 2016-06-09 | 小米科技有限責任公司Ｘｉａｏｍｉ Ｉｎｃ． | 直線検出方法、装置、プログラム、及び記録媒体 |
| JP2018110945A (ja) * | 2018-04-19 | 2018-07-19 | コニカミノルタ株式会社 | 画像処理装置及び照射野認識方法 |
| CN110490887A (zh) * | 2019-06-26 | 2019-11-22 | 河南埃尔森智能科技有限公司 | 一种基于3d视觉的对矩形包裹边缘快速识别定位方法 |
| JP2022008868A (ja) * | 2018-08-14 | 2022-01-14 | キヤノン株式会社 | 画像処理装置、画像処理方法、及びプログラム |
| JP2022020078A (ja) * | 2020-06-25 | 2022-02-01 | キヤノン株式会社 | 画像処理装置、画像処理方法、及び、プログラム |
| US12062181B2 (en) | 2018-08-14 | 2024-08-13 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, and storage medium for extracting an irradiation field of a radiograph |
Families Citing this family (15)
| 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 |
| US7558440B2 (en) * | 1998-09-30 | 2009-07-07 | Canon Kabushiki Kaisha | Image processing apparatus and method, photographing system, controlling method for the system, and computer-readable memory |
| JP4219571B2 (ja) * | 2001-04-26 | 2009-02-04 | 富士フイルム株式会社 | 放射線画像の黒化処理方法および装置 |
| US7397953B2 (en) * | 2001-07-24 | 2008-07-08 | Hewlett-Packard Development Company, L.P. | Image block classification based on entropy of differences |
| US7813576B2 (en) * | 2002-11-27 | 2010-10-12 | Fujifilm Corporation | Image processing apparatus |
| JP4708712B2 (ja) * | 2004-02-04 | 2011-06-22 | キヤノン株式会社 | 画像処理装置及びその制御方法、プログラム |
| DE102004008979B4 (de) * | 2004-02-24 | 2006-12-28 | Siemens Ag | Verfahren zur Filterung tomographischer 3D-Darstellungen nach erfolgter Rekonstruktion von Volumendaten |
| JP4280729B2 (ja) * | 2005-05-31 | 2009-06-17 | キヤノン株式会社 | 照射野領域抽出方法及び放射線撮影装置 |
| US7801344B2 (en) * | 2006-12-01 | 2010-09-21 | Carestream Health, Inc. | Edge boundary definition for radiographic detector |
| US8340387B2 (en) * | 2007-04-13 | 2012-12-25 | Three Palm Software | Fast preprocessing algorithms for digital mammography CAD and workstation |
| US20090154782A1 (en) * | 2007-12-17 | 2009-06-18 | Three Palm Software | Dual-magnify-glass visualization for soft-copy mammography viewing |
| CN101540040B (zh) * | 2008-03-21 | 2012-12-12 | 深圳迈瑞生物医疗电子股份有限公司 | 自动检测限束器边界的方法与装置 |
| WO2009124130A1 (en) * | 2008-04-01 | 2009-10-08 | Infimed, Inc. | Multi-pass edge detection of collimator blades in digital radiography images |
| JP5914404B2 (ja) | 2012-04-12 | 2016-05-11 | 富士フイルム株式会社 | Ｘ線露出制御装置、ｘ線画像検出装置及びｘ線画像撮影システム |
| CN111915569B (zh) * | 2020-07-09 | 2022-04-22 | 西安交通大学 | 自由曲面类零件数字射线图像区域筛选方法、设备及介质 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5512429A (en) * | 1978-07-12 | 1980-01-29 | Fuji Photo Film Co Ltd | Radioactive image reader |
| JP2613049B2 (ja) * | 1987-03-31 | 1997-05-21 | 富士写真フイルム株式会社 | 放射線照射野認識方法 |
| JP2596744B2 (ja) * | 1987-04-16 | 1997-04-02 | 富士写真フイルム株式会社 | 放射線照射野認識方法 |
| DE69429952T2 (de) * | 1993-07-22 | 2002-10-24 | Koninklijke Philips Electronics N.V., Eindhoven | Bildverarbeitungsverfahren und Vorrichtung zur Durchführung dieses Verfahrens |
| DE69518829T2 (de) * | 1995-05-11 | 2001-03-15 | Agfa-Gevaert N.V., Mortsel | Bestrahlungsfeldererkennungsverfahren |
| US5901240A (en) * | 1996-12-12 | 1999-05-04 | Eastman Kodak Company | Method for detecting the collimation field in a digital radiography |
-
1997
- 1997-03-31 JP JP07969697A patent/JP3923131B2/ja not_active Expired - Fee Related
-
1998
- 1998-03-31 US US09/050,353 patent/US6356651B2/en not_active Expired - Lifetime
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003265462A (ja) * | 2002-03-19 | 2003-09-24 | Hitachi Ltd | 関心領域抽出方法及び画像処理サーバ |
| US7796584B2 (en) | 2003-09-30 | 2010-09-14 | Nec Corporation | Method for connection between communication networks of different types and gateway apparatus |
| JP2006181362A (ja) * | 2004-12-24 | 2006-07-13 | General Electric Co <Ge> | ディジタル画像放射線撮像においてコリメーション・エッジを検出するシステム、方法及び装置 |
| JP2006293522A (ja) * | 2005-04-07 | 2006-10-26 | Fuji Photo Film Co Ltd | 直線検出装置、直線検出方法およびそのプログラム |
| US7991212B2 (en) | 2008-01-09 | 2011-08-02 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, and computer program |
| EP2372395A2 (en) | 2010-03-31 | 2011-10-05 | Fujifilm Corporation | Radiographic imaging system |
| US8238633B2 (en) | 2010-09-29 | 2012-08-07 | Fujifilm Corporation | Image processing apparatus, radiographic image system, and image processing method and program |
| JP2015123157A (ja) * | 2013-12-26 | 2015-07-06 | コニカミノルタ株式会社 | 画像処理装置及び照射野認識方法 |
| JP2016517100A (ja) * | 2014-02-22 | 2016-06-09 | 小米科技有限責任公司Ｘｉａｏｍｉ Ｉｎｃ． | 直線検出方法、装置、プログラム、及び記録媒体 |
| US9449245B2 (en) | 2014-02-22 | 2016-09-20 | Xiaomi Inc. | Method and device for detecting straight line |
| JP2018110945A (ja) * | 2018-04-19 | 2018-07-19 | コニカミノルタ株式会社 | 画像処理装置及び照射野認識方法 |
| JP2022008868A (ja) * | 2018-08-14 | 2022-01-14 | キヤノン株式会社 | 画像処理装置、画像処理方法、及びプログラム |
| US12062181B2 (en) | 2018-08-14 | 2024-08-13 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, and storage medium for extracting an irradiation field of a radiograph |
| CN110490887A (zh) * | 2019-06-26 | 2019-11-22 | 河南埃尔森智能科技有限公司 | 一种基于3d视觉的对矩形包裹边缘快速识别定位方法 |
| CN110490887B (zh) * | 2019-06-26 | 2023-02-28 | 河南埃尔森智能科技有限公司 | 一种基于3d视觉的对矩形包裹边缘快速识别定位方法 |
| JP2022020078A (ja) * | 2020-06-25 | 2022-02-01 | キヤノン株式会社 | 画像処理装置、画像処理方法、及び、プログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| US6356651B2 (en) | 2002-03-12 |
| JP3923131B2 (ja) | 2007-05-30 |
| US20010046312A1 (en) | 2001-11-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3923131B2 (ja) | 放射線画像の照射野認識方法および装置、並びに黒化処理方法および装置 | |
| EP0610605B1 (en) | Method of recognising an irradiation field | |
| EP2079054B1 (en) | Detection of blobs in images | |
| US20030161506A1 (en) | Face detection computer program product for redeye correction | |
| JP2574181B2 (ja) | 異常陰影検出装置 | |
| JP3617751B2 (ja) | 放射線画像の照射野認識方法および装置 | |
| JP2961437B2 (ja) | 画像出力方法 | |
| JP3347865B2 (ja) | 放射線画像の位置合せ方法 | |
| CN117935246A (zh) | 基于深度学习技术的仪表示数识别方法及装置 | |
| JP2582667B2 (ja) | 線状パターン幅演算装置 | |
| JP2006293522A (ja) | 直線検出装置、直線検出方法およびそのプログラム | |
| JP2571132B2 (ja) | 異常陰影検出装置 | |
| JP3165538B2 (ja) | 放射線画像の位置合せ方法および装置 | |
| JP3661797B2 (ja) | 放射線画像の位置合せ用対応点決定方法 | |
| JP2582666B2 (ja) | 異常陰影検出装置 | |
| JPH08215183A (ja) | 石灰化検出方法および装置 | |
| JP2867036B2 (ja) | 異常陰影検出装置 | |
| JP2582665B2 (ja) | 異常陰影検出装置 | |
| JPH04341247A (ja) | 胸部放射線画像の撮影体位判別装置 | |
| JP2987633B2 (ja) | 異常陰影検出装置 | |
| JP2006288467A (ja) | 照射野判定装置、照射野判定方法およびそのプログラム | |
| JP2006254934A (ja) | 照射野判定装置、照射野判定方法およびそのプログラム | |
| JP2582660B2 (ja) | 異常陰影検出装置 | |
| JPH0327485A (ja) | 異常陰影判定装置 | |
| Matthieu et al. | Image augmentation and automated measurement of endotracheal-tube-to-carina distance on chest radiographs in intensive care unit using a deep learning model with external validation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20060714 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060725 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20061024 |
|
| A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20061128 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20061225 |
|
| A911 | Transfer to examiner for re-examination before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20070115 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20070220 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070221 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100302 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110302 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110302 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120302 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120302 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130302 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130302 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140302 Year of fee payment: 7 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |