JPH0458941A - デジタル放射線画像の撮影体位判別装置 - Google Patents
デジタル放射線画像の撮影体位判別装置Info
- Publication number
- JPH0458941A JPH0458941A JP2168537A JP16853790A JPH0458941A JP H0458941 A JPH0458941 A JP H0458941A JP 2168537 A JP2168537 A JP 2168537A JP 16853790 A JP16853790 A JP 16853790A JP H0458941 A JPH0458941 A JP H0458941A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- level distribution
- image
- body position
- radiation image
- 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
- 230000005855 radiation Effects 0.000 title claims abstract description 39
- 238000009499 grossing Methods 0.000 claims abstract description 20
- 238000003384 imaging method Methods 0.000 claims description 32
- 238000001514 detection method Methods 0.000 claims description 17
- 238000013139 quantization Methods 0.000 claims description 6
- 238000000034 method Methods 0.000 description 22
- 238000004364 calculation method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 210000004072 lung Anatomy 0.000 description 5
- 210000000115 thoracic cavity Anatomy 0.000 description 5
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 4
- 210000000038 chest Anatomy 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000002902 bimodal effect Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 210000001370 mediastinum Anatomy 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Landscapes
- Image Analysis (AREA)
- Radiography Using Non-Light Waves (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〈産業上の利用分野・〉
本発明はデジタル放射線画像の撮影体位判別装置に関し
、詳しくは、被写体を透過した放射線量として撮影され
る放射線画像をデジタル信号化して各種の信号処理を行
うときの処理条件決定のために撮影体位を自動的に判別
する装置に関する。
、詳しくは、被写体を透過した放射線量として撮影され
る放射線画像をデジタル信号化して各種の信号処理を行
うときの処理条件決定のために撮影体位を自動的に判別
する装置に関する。
〈従来の技術〉
X線画像のような放射線画像は医療用として多(用いら
れており、この放射線画像を電気画像信号として得る方
法として、例えば以下のようなものがある。
れており、この放射線画像を電気画像信号として得る方
法として、例えば以下のようなものがある。
即ち、被写体としての人体を通過した放射線をある種の
蛍光体に吸収させ、その後、この蛍光体を光又は熱エネ
ルギーて励起することにより、この蛍光体が前記吸収に
より蓄積していた放射線エネルギーを蛍光として放射さ
せ、この蛍光を光電変換素子で検出して放射線画像情報
を電気的に得るものであり、かかる放射線画像信号をデ
ジタル化してから、階調処理や空間周波数処理を施して
CRT等に出力して可視化するようにしている(特開昭
63−189853号公報等参照)。
蛍光体に吸収させ、その後、この蛍光体を光又は熱エネ
ルギーて励起することにより、この蛍光体が前記吸収に
より蓄積していた放射線エネルギーを蛍光として放射さ
せ、この蛍光を光電変換素子で検出して放射線画像情報
を電気的に得るものであり、かかる放射線画像信号をデ
ジタル化してから、階調処理や空間周波数処理を施して
CRT等に出力して可視化するようにしている(特開昭
63−189853号公報等参照)。
ところで、上記のようにデジタル放射線画像信号を可視
像として出力する前の処理段階においては、上記のよう
に階調処理や空間周波数処理を施して読影に適した可視
像とする必要があるが、人体における同一部位の画像デ
ータであっても、撮影体位(被写体である人体の向き)
を変えて撮影する場合には、それぞれの再生画像におい
て該部位中の関心領域の濃度か変化してしまうことかあ
る。
像として出力する前の処理段階においては、上記のよう
に階調処理や空間周波数処理を施して読影に適した可視
像とする必要があるが、人体における同一部位の画像デ
ータであっても、撮影体位(被写体である人体の向き)
を変えて撮影する場合には、それぞれの再生画像におい
て該部位中の関心領域の濃度か変化してしまうことかあ
る。
例えば胸椎を診断するために、第5図に示すように人体
胸部を正面から撮影した場合には、関心領域である胸椎
は放射線の透過し難い縦隔部と重なるが、第6図に示す
ように側面から撮影した場合には、放射線の透過し易い
肺野と重なることになる。従って、これらのような画像
データを同一の画像処理条件で処理してから再生させる
と、正面画像においては胸椎部分か比較的低濃度となり
、側面画像においては比較的高濃度となってしまう。
胸部を正面から撮影した場合には、関心領域である胸椎
は放射線の透過し難い縦隔部と重なるが、第6図に示す
ように側面から撮影した場合には、放射線の透過し易い
肺野と重なることになる。従って、これらのような画像
データを同一の画像処理条件で処理してから再生させる
と、正面画像においては胸椎部分か比較的低濃度となり
、側面画像においては比較的高濃度となってしまう。
そこで、前記画像データかどのような体位で撮影された
かを自動的に判別することによって、それぞれの体位に
最適な画像処理を施す方法及び装置が、種々提案されて
いる。
かを自動的に判別することによって、それぞれの体位に
最適な画像処理を施す方法及び装置が、種々提案されて
いる。
この方法及び装置としては、画像データの濃度ヒス1〜
グラムを用いるものと、画像データ中の画像所定方向に
沿った信号レベル分布(以下、プロジェクションという
。)を用いるものとかあり、それぞれについて、累積操
作を行うもの(特開昭63−262128号公報,特開
昭61−262132号公報1特開昭63−26213
4号公報特開昭63−262139号公報等参照)、双
峰性ピークの分離度を求めるもの(特開昭63−262
129号公報,特開昭63−262135号公報等参照
)、関数近似を行うもの(特開昭63106642号公
報、特開昭63−262130号公報、特開昭63−2
62131号公報、特開昭63−262137号公報、
特開昭63−262138号公報等参照)かある。
グラムを用いるものと、画像データ中の画像所定方向に
沿った信号レベル分布(以下、プロジェクションという
。)を用いるものとかあり、それぞれについて、累積操
作を行うもの(特開昭63−262128号公報,特開
昭61−262132号公報1特開昭63−26213
4号公報特開昭63−262139号公報等参照)、双
峰性ピークの分離度を求めるもの(特開昭63−262
129号公報,特開昭63−262135号公報等参照
)、関数近似を行うもの(特開昭63106642号公
報、特開昭63−262130号公報、特開昭63−2
62131号公報、特開昭63−262137号公報、
特開昭63−262138号公報等参照)かある。
更に、濃度ヒストグラムから代表値や特性値を求めるも
の(特開昭63−106641号公報等参照)や、濃度
ヒストグラムの分散を求めるもの(特開昭63−262
133号公報等参照)、プロジェクションの左右方向の
中央部付近の信号値の総和を求めるもの(特開昭63−
262136号公報等参照)などが提案されている。
の(特開昭63−106641号公報等参照)や、濃度
ヒストグラムの分散を求めるもの(特開昭63−262
133号公報等参照)、プロジェクションの左右方向の
中央部付近の信号値の総和を求めるもの(特開昭63−
262136号公報等参照)などが提案されている。
〈発明か解決しようとする課題〉
しかしながら、上記の方法では、濃度ヒストダラムを用
いる場合でも、また、プロジェクションを用いる場合で
も、どのような体位で撮影されたかを自動的に判別する
ための計算量の多さか問題となっていた。例えば、多く
の画像データから正面画像データの標準パターンと側面
画像データの標準パターンとを得ておいて、この標準パ
ターンとその後入力されるデータとの類似度を調べるも
の(特開昭63−262130号公報、特開昭63−2
62137号公報等参照)では、標準パターンを得るた
めに重回帰分析等の複雑な手法を用いなけれはならず、
計算量か多くなってしまうものである。
いる場合でも、また、プロジェクションを用いる場合で
も、どのような体位で撮影されたかを自動的に判別する
ための計算量の多さか問題となっていた。例えば、多く
の画像データから正面画像データの標準パターンと側面
画像データの標準パターンとを得ておいて、この標準パ
ターンとその後入力されるデータとの類似度を調べるも
の(特開昭63−262130号公報、特開昭63−2
62137号公報等参照)では、標準パターンを得るた
めに重回帰分析等の複雑な手法を用いなけれはならず、
計算量か多くなってしまうものである。
また、濃度ヒストグラムを用いて撮影体位を判別する場
合には、濃度ヒストグラムの形状的な差異を捉えて撮影
体位を判別させることになるが、胸部正面画像データの
濃度ヒストグラムと、胸部側面画像データの濃度ヒスト
グラムとの間の形状的な差異か小さいため、判断の誤り
か多くなってしまう。
合には、濃度ヒストグラムの形状的な差異を捉えて撮影
体位を判別させることになるが、胸部正面画像データの
濃度ヒストグラムと、胸部側面画像データの濃度ヒスト
グラムとの間の形状的な差異か小さいため、判断の誤り
か多くなってしまう。
更に、プロジェクションを用いて撮影体位を判別させる
場合には、撮影範囲内における人体の横方向の位置ずれ
かあると、該画像データに基づくプロジェクションのパ
ターンか同じ撮影体位であっても変化してしまうため、
撮影体位の判別精度を確保することかできなくなってし
まうという問題かある。
場合には、撮影範囲内における人体の横方向の位置ずれ
かあると、該画像データに基づくプロジェクションのパ
ターンか同じ撮影体位であっても変化してしまうため、
撮影体位の判別精度を確保することかできなくなってし
まうという問題かある。
本発明は上記問題点に鑑みなされたものであり、デジタ
ル放射線画像における被写体の撮影体位を少ない計算量
で然も精度良く判別できる装置を提供することを目的と
する。
ル放射線画像における被写体の撮影体位を少ない計算量
で然も精度良く判別できる装置を提供することを目的と
する。
〈課題を解決するための手段〉
そのため本発明では、第1図に示すように、デジタル放
射線画像信号の画像所定方向に沿ったレベル分布を検出
するレベル分布検出手段と、このレベル分布検出手段で
検出されたレベル分布の信号レベルと複数種の信号レベ
ル閾値とをそれぞれに比較し、それぞれの閾値に対する
信号レベルの大小の変化パターンを前記画像所定方向に
沿って求める大小変化パターン検出手段と、この大小変
化パターン検出手段で検出されたそれぞれの閾値に対す
る大小変化パターンの種類別の数に基づいてデジタル放
射線画像信号における被写体の撮影体位を判別する撮影
体位判別手段と、を含んでデジタル放射線画像の撮影体
位判別装置を構成するようにした。
射線画像信号の画像所定方向に沿ったレベル分布を検出
するレベル分布検出手段と、このレベル分布検出手段で
検出されたレベル分布の信号レベルと複数種の信号レベ
ル閾値とをそれぞれに比較し、それぞれの閾値に対する
信号レベルの大小の変化パターンを前記画像所定方向に
沿って求める大小変化パターン検出手段と、この大小変
化パターン検出手段で検出されたそれぞれの閾値に対す
る大小変化パターンの種類別の数に基づいてデジタル放
射線画像信号における被写体の撮影体位を判別する撮影
体位判別手段と、を含んでデジタル放射線画像の撮影体
位判別装置を構成するようにした。
ここで、撮影体位判別手段が、予め記憶されている基準
の大小変化パターンと同じであるパターンの数と所定数
とを比較することによって被写体の撮影体位を判別する
よう構成することかできる。
の大小変化パターンと同じであるパターンの数と所定数
とを比較することによって被写体の撮影体位を判別する
よう構成することかできる。
また、撮影体位判別手段が、特定大小変化パターンの数
と所定数との比較による判別を複数組み合わせて被写体
の撮影体位を判別するよう構成することが好ましい。
と所定数との比較による判別を複数組み合わせて被写体
の撮影体位を判別するよう構成することが好ましい。
一方、レベル分布検出手段が、レベル分布を検出する画
像所定方向に略直交する画素列それぞれの信号レベルの
合計値又は平均値を用いて信号レベルの分布を検出する
ようにすると良い。
像所定方向に略直交する画素列それぞれの信号レベルの
合計値又は平均値を用いて信号レベルの分布を検出する
ようにすると良い。
また、レベル分布検出手段が、レベル分布を検出する前
記画像所定方向に沿った画素列それぞれの信号レベルを
そのまま用いて信号レベルの分布を検出するよう構成す
ることもてきる。
記画像所定方向に沿った画素列それぞれの信号レベルを
そのまま用いて信号レベルの分布を検出するよう構成す
ることもてきる。
更に、第1図点線示のように、デジタル放射線画像の画
素数と信号の量子化ステップ数との少なくとも一方を縮
減させる信号間引き手段を備えるようにして、レベル分
布検出手段かこの信号間引き手段で画素数と信号の量子
化ステップ数との少なくとも一方か縮減されたデジタル
放射線画像信号に基づいてレベル分布を検出するよう構
成すると良い。
素数と信号の量子化ステップ数との少なくとも一方を縮
減させる信号間引き手段を備えるようにして、レベル分
布検出手段かこの信号間引き手段で画素数と信号の量子
化ステップ数との少なくとも一方か縮減されたデジタル
放射線画像信号に基づいてレベル分布を検出するよう構
成すると良い。
また、第1図点線示のように、デジタル放射線画像の濃
度ヒストグラムの分布か広がるようにコントラストを強
調する補正をデジタル放射線画像信号に施すコン1へラ
スト強調補正手段を備えるようにして、レベル分布検出
手段かこのコントラスト強調補正手段で補正されたデジ
タル放射線画像信号に基ついてレベル分布を検出するよ
うにすると良い。
度ヒストグラムの分布か広がるようにコントラストを強
調する補正をデジタル放射線画像信号に施すコン1へラ
スト強調補正手段を備えるようにして、レベル分布検出
手段かこのコントラスト強調補正手段で補正されたデジ
タル放射線画像信号に基ついてレベル分布を検出するよ
うにすると良い。
また、第1図点線示のように、レベル分布検出手段で検
出されたレベル分布を平滑化処理するレベル分布平滑化
手段を備えるようにして、大小変化パターン検出手段か
このレベル分布平滑化手段で平滑化されたレベル分布に
基づいて大小の変化パターンを求めるようにすることが
好ましい。
出されたレベル分布を平滑化処理するレベル分布平滑化
手段を備えるようにして、大小変化パターン検出手段か
このレベル分布平滑化手段で平滑化されたレベル分布に
基づいて大小の変化パターンを求めるようにすることが
好ましい。
〈作用〉
かかる構成のデジタル放射線画像の撮影体位判別装置に
よると、所定画像方向における信号レベルの分布の様子
、より具体的には、前記画像所定方向における信号レベ
ルの大小変化の様子を、1つの閾値で切って判別するの
ではなく、複数種の閾値て切ることによって、レベル分
布における起伏の様子を2次元的に推定できるようにし
たものであり、かかるレベル分布の起伏の様子を示す大
小変化パターンの種類別の数に基づいて撮影体位か判別
される。
よると、所定画像方向における信号レベルの分布の様子
、より具体的には、前記画像所定方向における信号レベ
ルの大小変化の様子を、1つの閾値で切って判別するの
ではなく、複数種の閾値て切ることによって、レベル分
布における起伏の様子を2次元的に推定できるようにし
たものであり、かかるレベル分布の起伏の様子を示す大
小変化パターンの種類別の数に基づいて撮影体位か判別
される。
ここで、大小変化パターンの種類別の数に基ついて被写
体の撮影体位を判別するに当たっては、予め記憶させて
おいた基準の大小変化パターンと同じであるパターンの
数と所定数とを比較すれば、前記基準の大小変化パター
ンを多く示すことになるレベル分布のときの撮影体位を
予め知っておくことで、基準の大小変化パターンに対応
する撮影体位であるか否かを判別てき、このような特定
大小変化パターンの数と所定数との比較による判別を組
み合わせることて、複数種の撮影体位の判別を精度良く
行うことか可能となる。
体の撮影体位を判別するに当たっては、予め記憶させて
おいた基準の大小変化パターンと同じであるパターンの
数と所定数とを比較すれば、前記基準の大小変化パター
ンを多く示すことになるレベル分布のときの撮影体位を
予め知っておくことで、基準の大小変化パターンに対応
する撮影体位であるか否かを判別てき、このような特定
大小変化パターンの数と所定数との比較による判別を組
み合わせることて、複数種の撮影体位の判別を精度良く
行うことか可能となる。
デジタル画像信号の所定画像方向のレベル分布を検出す
るときには、前記画像所定方向に略直交する画素列それ
ぞれの信号レベルの合計値又は平均値を用いても良いし
、前記画像所定方向に沿った画素列それぞれの信号レベ
ルをそのまま用いても良い。
るときには、前記画像所定方向に略直交する画素列それ
ぞれの信号レベルの合計値又は平均値を用いても良いし
、前記画像所定方向に沿った画素列それぞれの信号レベ
ルをそのまま用いても良い。
また、デジタル放射線画像の画素数と信号の量子化ステ
ップ数との少なくとも一方を縮減させたり、デジタル放
射線画像の濃度ヒストグラムの分布か広がるようにコン
トラストを強調する補正を施したりしてから、レベル分
布を検出するようにすれば、少ないデータ量に基づいて
大小変化パターンを検出でき、また、撮影体位の特徴を
容易に捉えることができる。
ップ数との少なくとも一方を縮減させたり、デジタル放
射線画像の濃度ヒストグラムの分布か広がるようにコン
トラストを強調する補正を施したりしてから、レベル分
布を検出するようにすれば、少ないデータ量に基づいて
大小変化パターンを検出でき、また、撮影体位の特徴を
容易に捉えることができる。
また、検出されたレベル分布を平滑化処理すれば、微小
なレベル変化に撮影体位の判別が影響されることを防止
できる。
なレベル変化に撮影体位の判別が影響されることを防止
できる。
〈実施例〉
以下に本発明の詳細な説明する。
第2図は本発明にかかるデジタル放射線画像の撮影体位
判別装置の一実施例を示すものである。
判別装置の一実施例を示すものである。
ここで、画像入力部1は、蓄積性蛍光体上に記憶(蓄積
)された放射線画像の潜像をレーザ光で走査して読み取
るよう構成されたデジタルX線撮影装置や、X線フィル
ムをデジタル信号化するイメージスキャナ、或いは、既
にデジタル放射線画像データか複数記録されているファ
イリンクシステムなどである。
)された放射線画像の潜像をレーザ光で走査して読み取
るよう構成されたデジタルX線撮影装置や、X線フィル
ムをデジタル信号化するイメージスキャナ、或いは、既
にデジタル放射線画像データか複数記録されているファ
イリンクシステムなどである。
画像入力部1から入力されるオリジナルのデジタル放射
線画像データf (x、 y)は、−旦フレームメモ
リ2にストアされ、画像処理部3に渡される。前記画像
データf (x、 y) (x=L 2゜*+、
Nx、y=1.2. ・+、Ny)は、X。
線画像データf (x、 y)は、−旦フレームメモ
リ2にストアされ、画像処理部3に渡される。前記画像
データf (x、 y) (x=L 2゜*+、
Nx、y=1.2. ・+、Ny)は、X。
yによって番地指定される画素毎の信号レベルを示すも
のであり、本実施例では、X方向(第3図における左右
方向)に伸びる画素列を「行」、X方向(第3図におけ
る上下方向)に伸びる画素列を「列」 と呼ぶことにす
る。
のであり、本実施例では、X方向(第3図における左右
方向)に伸びる画素列を「行」、X方向(第3図におけ
る上下方向)に伸びる画素列を「列」 と呼ぶことにす
る。
画像処理部3は、本発明にかかる撮影体位判別装置に必
要な構成を備えると共に、撮影体位の判別結果に基づい
て読影に適した可視像に再生されるように階調変換や周
波数強調2輪郭抽出などの種々の画像処理・解析を行う
出力用画像処理部4を備えている。
要な構成を備えると共に、撮影体位の判別結果に基づい
て読影に適した可視像に再生されるように階調変換や周
波数強調2輪郭抽出などの種々の画像処理・解析を行う
出力用画像処理部4を備えている。
特徴パターン記憶部5は、後で詳細に説明するように撮
影体位判別に必要なデータである基準パターンや該パタ
ーンの数を判別するための閾値Th、更に、信号レベル
判別用の閾値値T等を記憶している。
影体位判別に必要なデータである基準パターンや該パタ
ーンの数を判別するための閾値Th、更に、信号レベル
判別用の閾値値T等を記憶している。
前記出力用画像処理部4で処理された画像データが入力
される画像出力部6では、画像処理後の画像データ(デ
ジタル放射線画像信号)を、再びファイリングシステム
内に記録したり、CRT上に再生したり、プリンタによ
ってフィルム上に記録したりする。ここで、被写体(本
実施例では人体)の撮影体位を判別した後の画像データ
を、ファイリングシステムにファイリングする際には、
患者名・撮影部位名・撮影年月日などの画像情報と共に
、撮影体位を画像データに対応させて記録しておけば、
後の検索に便利であると共に、繰り返し撮影体位の判別
を行う必要がなくなる。
される画像出力部6では、画像処理後の画像データ(デ
ジタル放射線画像信号)を、再びファイリングシステム
内に記録したり、CRT上に再生したり、プリンタによ
ってフィルム上に記録したりする。ここで、被写体(本
実施例では人体)の撮影体位を判別した後の画像データ
を、ファイリングシステムにファイリングする際には、
患者名・撮影部位名・撮影年月日などの画像情報と共に
、撮影体位を画像データに対応させて記録しておけば、
後の検索に便利であると共に、繰り返し撮影体位の判別
を行う必要がなくなる。
次に、本発明にかかるデジタル放射線画像の撮影体位判
別装置の実施例を、前記画像処理部3に含まれる各処理
部11〜16に沿って説明する。
別装置の実施例を、前記画像処理部3に含まれる各処理
部11〜16に沿って説明する。
フレームメモリ2に一旦スドアされるオリジナルのデジ
タル放射線画像データf (x、 y)は、ます、画
像処理部3内の信号間引き部(信号間引き手段)11に
入力される。
タル放射線画像データf (x、 y)は、ます、画
像処理部3内の信号間引き部(信号間引き手段)11に
入力される。
この信号間引き部11は、オリジナルのデジタル画像デ
ータf (x、 y)の画素数及び階調数(量子化ス
テップ数)を間引いて、撮影体位判別に用いる画像デー
タ量を縮減させて以下の計算量を減少させるようにする
。
ータf (x、 y)の画素数及び階調数(量子化ス
テップ数)を間引いて、撮影体位判別に用いる画像デー
タ量を縮減させて以下の計算量を減少させるようにする
。
具体的には、オリジナルのデジタル画像データf(x、
y)を複数の小領域に分割し、各小領域毎の信号値の平
均値や代表値を求めることで画素数及び階調数を縮減す
る。例えは、2048 X 2464画素。
y)を複数の小領域に分割し、各小領域毎の信号値の平
均値や代表値を求めることで画素数及び階調数を縮減す
る。例えは、2048 X 2464画素。
1024階調のオリジナルのデジタル画像データf(x
、 y)を、128X154画素、256階調にまで
信号を間引く処理を行っても、本発明にかかる撮影体位
判別を支障なく行わせることかできる。
、 y)を、128X154画素、256階調にまで
信号を間引く処理を行っても、本発明にかかる撮影体位
判別を支障なく行わせることかできる。
尚、上記信号間引き部11で画素数及び階調数の間引き
が行われた画像データをf“(x’、 y’) (X′
=1.2.・・・、Nx’、 y′= 1.2.・・・
、Ny’) (Nx’ < Nx、 Ny’ < N
y)とする。
が行われた画像データをf“(x’、 y’) (X′
=1.2.・・・、Nx’、 y′= 1.2.・・・
、Ny’) (Nx’ < Nx、 Ny’ < N
y)とする。
但し、上記の間引き処理は、計算量か少なくて済む本発
明にかかる撮影体位判別装置における計算量を更に少な
くするためのものであり、間引き処理を省略しても良い
。
明にかかる撮影体位判別装置における計算量を更に少な
くするためのものであり、間引き処理を省略しても良い
。
次のコントラスト強調補正部(コントラスト強調補正手
段)12では、間引き処理後の画像データf’(x“、
y′)のコントラスト強調補正を必要に応じて行ってコ
ントラスト強調データf” (x’。
段)12では、間引き処理後の画像データf’(x“、
y′)のコントラスト強調補正を必要に応じて行ってコ
ントラスト強調データf” (x’。
y’)を得る。これは、原画像データの信号値(濃淡値
)には、あまり濃淡の差がないことが多いため、プロジ
ェクション(画像データ中の画像所定方向に沿った信号
レベル分布)を作成した際に、分布曲線か描く山部と谷
部との高低差か小さくなることを防止するために行う処
理であり、コントラストの強調補正を施すことによって
、濃度差を強調して撮影体位判別に関わる画像データの
濃淡変化の特徴か明確に表れるようにする。
)には、あまり濃淡の差がないことが多いため、プロジ
ェクション(画像データ中の画像所定方向に沿った信号
レベル分布)を作成した際に、分布曲線か描く山部と谷
部との高低差か小さくなることを防止するために行う処
理であり、コントラストの強調補正を施すことによって
、濃度差を強調して撮影体位判別に関わる画像データの
濃淡変化の特徴か明確に表れるようにする。
コントラストの強調は、例えば以下のような公知の技術
を用いて行う(「コンピュータ画像処理入門」田村秀行
監修、総研出版(株)発行等参照)。
を用いて行う(「コンピュータ画像処理入門」田村秀行
監修、総研出版(株)発行等参照)。
強調前の信号のレベル範囲を[a□in+ amax
l、強調後の信号範囲を[b、、、、、b、、、、]と
した場合、強調前の信号値S o r :□を以下の■
式又は0式によってコントラスト強調信号値S pr。
l、強調後の信号範囲を[b、、、、、b、、、、]と
した場合、強調前の信号値S o r :□を以下の■
式又は0式によってコントラスト強調信号値S pr。
。に変換する。
(logs oriz loga n1lo) +
b +nin”’■ここで、■式を用いてコントラスト
強調を行った場合には濃度か線形に変換され、0式では
非線形に濃度変換されることになる。
b +nin”’■ここで、■式を用いてコントラスト
強調を行った場合には濃度か線形に変換され、0式では
非線形に濃度変換されることになる。
また、下式■のように平均濃度値m、標準偏差σの濃度
ヒストグラムを、平均値mn。7.標準偏差σ。。、。
ヒストグラムを、平均値mn。7.標準偏差σ。。、。
の濃度ヒストグラムに正規化するように、強調前の信号
値S o t + tをコントラスト強調信号値S p
r。。に変換するようにしても良い。
値S o t + tをコントラスト強調信号値S p
r。。に変換するようにしても良い。
σ
更に、全濃度域の画素数か均等になるように、濃度ヒス
1〜クラムを平坦化することによってコントラスト 次のプロジェクション作成部(レベル分布検出手段)1
3では、信号間引き処理及びコントラス1〜強調補正処
理を施された画像データf” (x’ 。
1〜クラムを平坦化することによってコントラスト 次のプロジェクション作成部(レベル分布検出手段)1
3では、信号間引き処理及びコントラス1〜強調補正処
理を施された画像データf” (x’ 。
y′)に基づいてプロジェクションp(x’)(画像X
方向に沿ったレベル分布)を作成する。
方向に沿ったレベル分布)を作成する。
前記プロジェクションp (x’)は、画像データf
(x’ 、 y’)における次式で表されるような
分布を用いることができる(但し、第3図に示すように
1≦n1≦n2≦Ny’ として、信号のサンプリン
グ範囲をn1≦y′≦n2とする。)。
(x’ 、 y’)における次式で表されるような
分布を用いることができる(但し、第3図に示すように
1≦n1≦n2≦Ny’ として、信号のサンプリン
グ範囲をn1≦y′≦n2とする。)。
P + ( x ’ ) −’): f ” (
x ’ 、y ’ )・・・■即ち、■式では画像デー
タf” (x’ 、y’)の画像上下方向(X方向)に
伸びる画素列毎の信号値の総和が求められ、0式では画
像データf″(x.y’)の画像上下方向に伸びる各画
素列毎の信号値の平均値が求められる。尚、場合によっ
ては、画像データf′’ (x’ 、 y’)の左右
方向(X方向)に伸びる各画素行毎の信号値の総和や平
均値を用いることもできる。
x ’ 、y ’ )・・・■即ち、■式では画像デー
タf” (x’ 、y’)の画像上下方向(X方向)に
伸びる画素列毎の信号値の総和が求められ、0式では画
像データf″(x.y’)の画像上下方向に伸びる各画
素列毎の信号値の平均値が求められる。尚、場合によっ
ては、画像データf′’ (x’ 、 y’)の左右
方向(X方向)に伸びる各画素行毎の信号値の総和や平
均値を用いることもできる。
本実施例では、上記のように、各画素列の信号値の総和
又は平均値に基づいてプロジェクションを作成するよう
にしたが、特定画素行における各画素の信号値をそのま
ま用いて画像のX方向(左右方向)に沿った信号レベル
変化を示すプロジェクションを作成するようにしても良
い。但し、撮影体位の特徴が明確に表れるように、前記
特定画素行を適宜選択する必要がある。
又は平均値に基づいてプロジェクションを作成するよう
にしたが、特定画素行における各画素の信号値をそのま
ま用いて画像のX方向(左右方向)に沿った信号レベル
変化を示すプロジェクションを作成するようにしても良
い。但し、撮影体位の特徴が明確に表れるように、前記
特定画素行を適宜選択する必要がある。
このプロジェクション部13て得られるプロジェクショ
ンp (x’)には、第4図に示すように、信号値の微
小な変動が含まれているのて、後述する閾値Tに対する
大小変化のパターンを検出するときに、この微小な変化
を撮影体位の特徴を示す大小変化として捉えないように
前記プロジェクションp (x’)を平滑化(スムージ
ング)する処理を、次の平滑化処理部(レベル分布平滑
化手段)14で行う。
ンp (x’)には、第4図に示すように、信号値の微
小な変動が含まれているのて、後述する閾値Tに対する
大小変化のパターンを検出するときに、この微小な変化
を撮影体位の特徴を示す大小変化として捉えないように
前記プロジェクションp (x’)を平滑化(スムージ
ング)する処理を、次の平滑化処理部(レベル分布平滑
化手段)14で行う。
前記平滑化処理は、例えばプロジェクションp(X′)
内の注目点を、その近傍内における平均値で置き換える
ことによって行える。第4図に示すように、プロジェク
ションp (x’)内の各点の横方向における位置はx
′、その信号値はp (x’)という関数であるから、
左右n個ずつの近傍処理を行う場合には、位置iにおけ
る平滑処理後の値p’ (i)は、 p’(i)−Σ p (x’)/ (2xn+1)とい
う式で表すことができ、かかる処理を複数回実行するこ
とて第4図に示すような信号値の微小な変動を平滑化す
る。尚、前記平滑化処理部14における平滑化処理には
、近傍内の中央値を用いることもできる。
内の注目点を、その近傍内における平均値で置き換える
ことによって行える。第4図に示すように、プロジェク
ションp (x’)内の各点の横方向における位置はx
′、その信号値はp (x’)という関数であるから、
左右n個ずつの近傍処理を行う場合には、位置iにおけ
る平滑処理後の値p’ (i)は、 p’(i)−Σ p (x’)/ (2xn+1)とい
う式で表すことができ、かかる処理を複数回実行するこ
とて第4図に示すような信号値の微小な変動を平滑化す
る。尚、前記平滑化処理部14における平滑化処理には
、近傍内の中央値を用いることもできる。
次に、変化パターン検出部(大小変化パターン検出手段
)15では、前記平滑化処理されたプロジェクションp
’(x’)の各位と複数種の信号閾値下とを順次比較す
ることによって、各閾値Tに対する信号レベルの大小変
化パターンを求め、次の撮影体位判別部(撮影体位判別
手段)16では、この変動パターンの中に、特定の撮影
体位に特徴的に表れる特徴パターンと同じパターンの数
かいくつあるかによって撮影体位を判別する。
)15では、前記平滑化処理されたプロジェクションp
’(x’)の各位と複数種の信号閾値下とを順次比較す
ることによって、各閾値Tに対する信号レベルの大小変
化パターンを求め、次の撮影体位判別部(撮影体位判別
手段)16では、この変動パターンの中に、特定の撮影
体位に特徴的に表れる特徴パターンと同じパターンの数
かいくつあるかによって撮影体位を判別する。
即ち、人体胸部放射線画像の正面と側面との場合、第5
図及び第6図に直線して示す画像左右方向の信号値の分
布において、−殻内に正面画像においては、左右方向の
中央部に放射線の透過し難い胸椎を表す信号レベルの谷
か存在し、その両脇に放射線を透過し易い肺野を示す信
号レベルの山が存在し、更に、左右端に放射線を透過し
難い胸郭を表す信号レベルの谷が存在する。一方、側面
画像では、左右端に放射線の素抜は部を示す信号レベル
の山が存在し、中央部付近では体制を示す谷となって、
体側部においては特徴的な波形は表れ難い。かかる特徴
は、直線りに直交する方向の各画素列の信号値の総和又
は平均値の分布を取っても略同様である。
図及び第6図に直線して示す画像左右方向の信号値の分
布において、−殻内に正面画像においては、左右方向の
中央部に放射線の透過し難い胸椎を表す信号レベルの谷
か存在し、その両脇に放射線を透過し易い肺野を示す信
号レベルの山が存在し、更に、左右端に放射線を透過し
難い胸郭を表す信号レベルの谷が存在する。一方、側面
画像では、左右端に放射線の素抜は部を示す信号レベル
の山が存在し、中央部付近では体制を示す谷となって、
体側部においては特徴的な波形は表れ難い。かかる特徴
は、直線りに直交する方向の各画素列の信号値の総和又
は平均値の分布を取っても略同様である。
ここで、信号レベル閾値Tを変化させながらプロジェク
ションの信号レベルと閾値Tとを順次比較していくと、
プロジェクションの各位か閾値Tよりも小さい部分が続
いた後、大きい部分に変わり、再び小さい部分か続くと
いった変動のパターンか閾値T数だけサンプリングされ
、これらの変動パターンは第5図及び第6図から明らか
なように、撮影体位か正面であるか側面であるかによっ
て大きく異なる。
ションの信号レベルと閾値Tとを順次比較していくと、
プロジェクションの各位か閾値Tよりも小さい部分が続
いた後、大きい部分に変わり、再び小さい部分か続くと
いった変動のパターンか閾値T数だけサンプリングされ
、これらの変動パターンは第5図及び第6図から明らか
なように、撮影体位か正面であるか側面であるかによっ
て大きく異なる。
このように、正面画像に多く表れるパターンと、側面画
像に多く表れるパターンとは、それぞれ特定されるから
、これらの撮影体位に対応する特徴的なパターンを予め
求めておけば、入力した画像データのプロジェクション
と閾値Tとの比較によって得られる閾値T数分のパター
ンの中に、前記特徴パターンと同じパターンのものがい
くつあるかによって撮影体位を判別できるものである。
像に多く表れるパターンとは、それぞれ特定されるから
、これらの撮影体位に対応する特徴的なパターンを予め
求めておけば、入力した画像データのプロジェクション
と閾値Tとの比較によって得られる閾値T数分のパター
ンの中に、前記特徴パターンと同じパターンのものがい
くつあるかによって撮影体位を判別できるものである。
かかる大小変化パターンの検出及び撮影体位の判別を次
に具体的に説明する。
に具体的に説明する。
変化パターン検出部(大小変化パターン検出手段)15
では、前記プロジェクションp’(x’)の各点を端か
ら順に閾値Tと比較して、例えばプロジェクション上の
値か閾値Tよりも大きければ1、小さけれは0を与えて
いく。ここで、1或いは0が続く長さ、即ち、■及びO
のラン長を無視して、1のランに関して符号1を、0の
ランに符号0を割当てれは、符号lとOの繰り返しパタ
ーンが、そのプロジェクションの波形の閾値Tに対する
変動の様子を表していることになる。更に、閾値Tを予
め設定されたステップで変化させて、信号値のレンジ全
体にわたって上記のような1.0符号の設定を行えは、
そのプロジェクションにおける波形の様子を2次元的に
捉えることかでき、そのプロジェクションの特徴となる
1、0の変化パターンが多く表れることになる。
では、前記プロジェクションp’(x’)の各点を端か
ら順に閾値Tと比較して、例えばプロジェクション上の
値か閾値Tよりも大きければ1、小さけれは0を与えて
いく。ここで、1或いは0が続く長さ、即ち、■及びO
のラン長を無視して、1のランに関して符号1を、0の
ランに符号0を割当てれは、符号lとOの繰り返しパタ
ーンが、そのプロジェクションの波形の閾値Tに対する
変動の様子を表していることになる。更に、閾値Tを予
め設定されたステップで変化させて、信号値のレンジ全
体にわたって上記のような1.0符号の設定を行えは、
そのプロジェクションにおける波形の様子を2次元的に
捉えることかでき、そのプロジェクションの特徴となる
1、0の変化パターンが多く表れることになる。
例えは、第7図に示すように、信号値のレンジt0から
t6まての間て、閾値TをΔTの間隔で変化させていっ
た場合、下記表1に示すように、tn≦’I’< t、
+1 (n=0.t、 ・・・、5)の各区間で、
符号0と符号lの並びの特徴パターンと、その特徴パタ
ーンの個数が得られる。
t6まての間て、閾値TをΔTの間隔で変化させていっ
た場合、下記表1に示すように、tn≦’I’< t、
+1 (n=0.t、 ・・・、5)の各区間で、
符号0と符号lの並びの特徴パターンと、その特徴パタ
ーンの個数が得られる。
表1
かかる処理を、第8図〜第12図に示すそれぞれの放射
線画像データについて行うと、正面画像である第8図で
は(01010)、側面画像である第9図では(101
) 、正面右ずれ画像である第10図では(10101
) 、側面右ずれ画像である第11図では(10)、子
供や痩せた人の正面画像である第12図では(1010
101)のパターンかそれぞれ多く検出されることにな
り、それぞれの画像ては、上記以外のパターンは僅かに
しか存在しない。従って、(01010)。
線画像データについて行うと、正面画像である第8図で
は(01010)、側面画像である第9図では(101
) 、正面右ずれ画像である第10図では(10101
) 、側面右ずれ画像である第11図では(10)、子
供や痩せた人の正面画像である第12図では(1010
101)のパターンかそれぞれ多く検出されることにな
り、それぞれの画像ては、上記以外のパターンは僅かに
しか存在しない。従って、(01010)。
(10101)、 (1010101)のパターンを正
面画像用の特徴パターン(正面の基準大小変化パターン
)とし、(101)、 (10)のパターンを側面画像
用の特徴バタン(側面の基準大小変化パターン)として
予め前記特徴パターン記憶部5に記憶させておけば、入
力される画像データのプロジェクションが、どのパター
ンを多くもっているかという情報と、特徴パターン記憶
部5に記憶されている撮影体位に対応する特徴パターン
の情報とを比較することによって、撮影体位判別部(撮
影体位判別手段)16において、人体を正面から撮影し
た画像であるが、又は、人体を側面から撮影した画像で
あるかの撮影体位の判別か行えるものである。
面画像用の特徴パターン(正面の基準大小変化パターン
)とし、(101)、 (10)のパターンを側面画像
用の特徴バタン(側面の基準大小変化パターン)として
予め前記特徴パターン記憶部5に記憶させておけば、入
力される画像データのプロジェクションが、どのパター
ンを多くもっているかという情報と、特徴パターン記憶
部5に記憶されている撮影体位に対応する特徴パターン
の情報とを比較することによって、撮影体位判別部(撮
影体位判別手段)16において、人体を正面から撮影し
た画像であるが、又は、人体を側面から撮影した画像で
あるかの撮影体位の判別か行えるものである。
上記のプロジェクションの波形判別に用いる信号レベル
の閾値Tは、特徴パターン毎に予め求めておいて、前記
特徴パターン記憶部5に特徴パターンと共に記憶させて
おくと良い。そして、判別の際には、第14図のフロー
チャートに示すように、入力画像データの中のそれぞれ
の特徴パターンの数を、それぞれの閾値(Thl−Th
5)と順次比較していき、その人力画像データかある特
徴パターンを所定数以上もっているかどうかによって判
別を行っても良いし、第15図のフローチャートに示す
判別のようにいくつかの特徴パターンを一纏めにして閾
値(Th6,7)と比較して撮影体位を判別するように
しても良い。
の閾値Tは、特徴パターン毎に予め求めておいて、前記
特徴パターン記憶部5に特徴パターンと共に記憶させて
おくと良い。そして、判別の際には、第14図のフロー
チャートに示すように、入力画像データの中のそれぞれ
の特徴パターンの数を、それぞれの閾値(Thl−Th
5)と順次比較していき、その人力画像データかある特
徴パターンを所定数以上もっているかどうかによって判
別を行っても良いし、第15図のフローチャートに示す
判別のようにいくつかの特徴パターンを一纏めにして閾
値(Th6,7)と比較して撮影体位を判別するように
しても良い。
第14図のフローチャートに示す判別では、標準正面(
第8図参照)、標準側面(第9図参照)正面右ずれ(第
10図参照)、側面右ずれ(第11図参照)、正面小(
第12図参照)それぞれの特徴パターンの数を個別に判
別していって、前記5種類の撮影体位を判別できるもの
である。
第8図参照)、標準側面(第9図参照)正面右ずれ(第
10図参照)、側面右ずれ(第11図参照)、正面小(
第12図参照)それぞれの特徴パターンの数を個別に判
別していって、前記5種類の撮影体位を判別できるもの
である。
即ち、第14図のフローチャートは、撮影体位判別を行
うための前処理も含めて示してあり、Slで原画像デー
タを入力し、S2で信号間引き処理を行い、S3では更
にコントラストを強調させる処理を行う。そして、これ
らの処理を施された画像データに基づいてS4でプロジ
ェクションの作成を行って、次の85でこのプロジェク
ションを平滑化する処理を行う。
うための前処理も含めて示してあり、Slで原画像デー
タを入力し、S2で信号間引き処理を行い、S3では更
にコントラストを強調させる処理を行う。そして、これ
らの処理を施された画像データに基づいてS4でプロジ
ェクションの作成を行って、次の85でこのプロジェク
ションを平滑化する処理を行う。
S6では平滑化されたプロジェクションの信号レベルと
閾値Tとを、閾値Tを変化させながら比較していって各
閾値T毎に信号レベルの大小変化パターンを検出する。
閾値Tとを、閾値Tを変化させながら比較していって各
閾値T毎に信号レベルの大小変化パターンを検出する。
検出された各閾値T毎の大小変化パターンは、同じパタ
ーンであるものの数か求められ、S7以降の撮影体位判
別に用いられる。
ーンであるものの数か求められ、S7以降の撮影体位判
別に用いられる。
S7ては、第8図に示すような標準正面体位であるとき
に特徴的に表れるパターンである(01010)と同じ
であったパターンの数と、パターン数の閾値、Thlと
を比較し、Th1以上に(01010)のパターンが表
れたときには、S8で撮影体位か第8図に示すような標
準正面であると判別する。
に特徴的に表れるパターンである(01010)と同じ
であったパターンの数と、パターン数の閾値、Thlと
を比較し、Th1以上に(01010)のパターンが表
れたときには、S8で撮影体位か第8図に示すような標
準正面であると判別する。
同様な判別か89〜S1.6で行われ、S9では標準側
面の特徴パターン(101)の数が判別され、S11で
は正面右ずれの特徴パターン(10101)の数か判別
され、S13ては側面右ずれの特徴パターン(10)の
数か判別され、更に、S15では正面小の特徴パターン
である( +、 01.0101 )の数か判別され、
それぞれ閾値Th2〜Th5よりも多いときには、対応
する特徴パターンの体位(側面標準、正面右ずれ、側面
右ずれ、正面小のいずれか)で撮影された画像であると
判別する( S 10.12.14.16)。
面の特徴パターン(101)の数が判別され、S11で
は正面右ずれの特徴パターン(10101)の数か判別
され、S13ては側面右ずれの特徴パターン(10)の
数か判別され、更に、S15では正面小の特徴パターン
である( +、 01.0101 )の数か判別され、
それぞれ閾値Th2〜Th5よりも多いときには、対応
する特徴パターンの体位(側面標準、正面右ずれ、側面
右ずれ、正面小のいずれか)で撮影された画像であると
判別する( S 10.12.14.16)。
そして、上記5種類の撮影体位に特徴的に表れるパター
ンと同じであるパターンの数かいずれも閾値Thl〜5
よりも小さいときには、正面標準。
ンと同じであるパターンの数かいずれも閾値Thl〜5
よりも小さいときには、正面標準。
側面標準、正面右ずれ、側面右ずれ、正面小のいずれで
もなく、S17で「その他」の撮影体位であると判別す
る。
もなく、S17で「その他」の撮影体位であると判別す
る。
一方、第15図のフローチャートに示す判別では、判別
を簡略化し、位置ずれの有無に関わらず撮影体位を正面
と側面との2種類に判別する。即ち、前処理の321−
326については第14図のフローチャートにおける8
1〜S6と同様であるが、撮影体位判別におけるS27
において、標準正面に対応するパターン数と、正面右ず
れに対応するパターンの数と、正面小に対応するパター
ンの数とか加算され、これらのいずれも正面であるとき
に表れるパターンの総和か所定の閾値Th6を越えると
きには、S28で標準正面、正面右ずれ、正面小の区別
を行うことなく、撮影体位が正面であると判別する。
を簡略化し、位置ずれの有無に関わらず撮影体位を正面
と側面との2種類に判別する。即ち、前処理の321−
326については第14図のフローチャートにおける8
1〜S6と同様であるが、撮影体位判別におけるS27
において、標準正面に対応するパターン数と、正面右ず
れに対応するパターンの数と、正面小に対応するパター
ンの数とか加算され、これらのいずれも正面であるとき
に表れるパターンの総和か所定の閾値Th6を越えると
きには、S28で標準正面、正面右ずれ、正面小の区別
を行うことなく、撮影体位が正面であると判別する。
そして、S27で正面パターンの数か閾値Th6以下で
あると判別されたときには、S29で標準側面に対応す
るパターン数と側面右ずれに対応するパターン数との総
和と閾値Th7とを比較し、閾値Th7を越えるときに
は、標準側面、側面右ずれの区別を行うことなくS30
で撮影体位が側面であると判別する。
あると判別されたときには、S29で標準側面に対応す
るパターン数と側面右ずれに対応するパターン数との総
和と閾値Th7とを比較し、閾値Th7を越えるときに
は、標準側面、側面右ずれの区別を行うことなくS30
で撮影体位が側面であると判別する。
S29で撮影体位か側面でもないと判別されると、S3
1て撮影体位が、正面及び側面のいずれでもない「その
他」であると判別する。
1て撮影体位が、正面及び側面のいずれでもない「その
他」であると判別する。
更に、例えは第13図に示すように、片肺か切除されて
いるような特殊なケースの画像データであっても、その
特徴パターン(010)を、特徴パターン記憶5におけ
る特徴パターンの集合に付加して、かかる片肺の人体に
対応する特徴パターンの数と該パターン数の閾値とを比
較することで、片肺の人体を正面から撮影した体位につ
いても判別が可能である。
いるような特殊なケースの画像データであっても、その
特徴パターン(010)を、特徴パターン記憶5におけ
る特徴パターンの集合に付加して、かかる片肺の人体に
対応する特徴パターンの数と該パターン数の閾値とを比
較することで、片肺の人体を正面から撮影した体位につ
いても判別が可能である。
上記のようにして撮影体位の判別か行われると、最後に
、撮影体位の判別結果を表す信号か撮影体位判別部16
から出力用画像処理部4に出力される。
、撮影体位の判別結果を表す信号か撮影体位判別部16
から出力用画像処理部4に出力される。
そして、出力用画像処理部4では、フレームメモリ2内
のオリジナル画像データf (x、y)を直接取り込み
、撮影体位に適合する変換方法で階調変換や周波数強調
なとの処理を行って画像データF (x、 y)を得
て、この画像データF(xy)を画像出力部6に出力す
る。
のオリジナル画像データf (x、y)を直接取り込み
、撮影体位に適合する変換方法で階調変換や周波数強調
なとの処理を行って画像データF (x、 y)を得
て、この画像データF(xy)を画像出力部6に出力す
る。
画像出力部6では、出力用画像処理部4から撮影体位に
応じた処理を施された画像データF(xy)をファイリ
ングシステム内に記録したり、CRT上に再生したり、
プリンタによってX線フィルム上に記録したりする。こ
こで、画像出力部6ては、画像データかその画像の撮影
体位に応じた処理を施されているから、再生される画像
を読影に適した可視像とすることかできる。
応じた処理を施された画像データF(xy)をファイリ
ングシステム内に記録したり、CRT上に再生したり、
プリンタによってX線フィルム上に記録したりする。こ
こで、画像出力部6ては、画像データかその画像の撮影
体位に応じた処理を施されているから、再生される画像
を読影に適した可視像とすることかできる。
〈発明の効果〉
以上説明したように、本発明にかかるデジタル放射線画
像の撮影体位判別装置によると、デジタル放射線画像デ
ータにおける被写体の大きさや位置ずれに影響されずに
、広範囲に渡る種々の画像データから被写体の撮影体位
を簡単かつ精度良く判別することかでき、かかる撮影体
位の判別結果に基づいてその撮影体位に最適な画像処理
を施して、再生画像における関心領域の濃度を一定に揃
えることが可能となり、放射線画像に基づく医療診断性
を高めることかできるという効果かある。
像の撮影体位判別装置によると、デジタル放射線画像デ
ータにおける被写体の大きさや位置ずれに影響されずに
、広範囲に渡る種々の画像データから被写体の撮影体位
を簡単かつ精度良く判別することかでき、かかる撮影体
位の判別結果に基づいてその撮影体位に最適な画像処理
を施して、再生画像における関心領域の濃度を一定に揃
えることが可能となり、放射線画像に基づく医療診断性
を高めることかできるという効果かある。
第1図は本発明の構成を示すフロック図、第2図は本発
明の一実施例を示すシステムブロック図、第3図は同上
実施例における画像データの表し方を示す線図、第4図
は同上実施例における平滑処理前のプロジェクションの
様子を示す線図、第5図及び第′6図はそれぞれ撮影体
位による信号レベル分布の特徴の違いを説明するための
線図、第7図は同上実施例におけるプロジェクションに
基つく大小変化パターンの検出の様子を示す線図、第8
図〜第13図はそれぞれ撮影体位による信号レベルの大
小変化パターンの違いを説明するための線図、第14図
及び115図はそれぞれ大小変化パターンに違いに基つ
く撮影体位判別の内容を示すフロ−チャートである。 l・・・画像入力部 2・・・フレームメモリ3・・
・画像処理部 4・・・出力用画像処理部5・・・特
徴パターン記憶部 6・・・画像出力部11・・・信
号間引き部 12・・・コントラスI・強調補正部
13・・・プロジェクション作成部 14・・・平
滑化処理部 15・・・変化パターン検出部 16
・・撮影体位判別部
明の一実施例を示すシステムブロック図、第3図は同上
実施例における画像データの表し方を示す線図、第4図
は同上実施例における平滑処理前のプロジェクションの
様子を示す線図、第5図及び第′6図はそれぞれ撮影体
位による信号レベル分布の特徴の違いを説明するための
線図、第7図は同上実施例におけるプロジェクションに
基つく大小変化パターンの検出の様子を示す線図、第8
図〜第13図はそれぞれ撮影体位による信号レベルの大
小変化パターンの違いを説明するための線図、第14図
及び115図はそれぞれ大小変化パターンに違いに基つ
く撮影体位判別の内容を示すフロ−チャートである。 l・・・画像入力部 2・・・フレームメモリ3・・
・画像処理部 4・・・出力用画像処理部5・・・特
徴パターン記憶部 6・・・画像出力部11・・・信
号間引き部 12・・・コントラスI・強調補正部
13・・・プロジェクション作成部 14・・・平
滑化処理部 15・・・変化パターン検出部 16
・・撮影体位判別部
Claims (8)
- (1)デジタル放射線画像信号の画像所定方向に沿った
レベル分布を検出するレベル分布検出手段と、該レベル
分布検出手段で検出されたレベル分布の信号レベルと複
数種の信号レベル閾値とをそれぞれに比較し、それぞれ
の閾値に対する信号レベルの大小の変化パターンを前記
画像所定方向に沿って求める大小変化パターン検出手段
と、 該大小変化パターン検出手段で検出されたそれぞれの閾
値に対する大小変化パターンの種類別の数に基づいてデ
ジタル放射線画像信号における被写体の撮影体位を判別
する撮影体位判別手段と、を含んで構成されたことを特
徴とするデジタル放射線画像の撮影体位判別装置。 - (2)前記撮影体位判別手段が、予め記憶されている基
準の大小変化パターンと同じであるパターンの数と所定
数とを比較することによって被写体の撮影体位を判別す
ることを特徴とする請求項1記載のデジタル放射線画像
の撮影体位判別装置。 - (3)前記撮影体位判別手段が、特定大小変化パターン
の数と所定数との比較による判別を複数組み合わせて被
写体の撮影体位を判別することを特徴とする請求項1又
は2のいずれかに記載のデジタル放射線画像の撮影体位
判別装置。 - (4)前記レベル分布検出手段が、レベル分布を検出す
る前記画像所定方向に略直交する画素列それぞれの信号
レベルの合計値又は平均値を用いて信号レベルの分布を
検出することを特徴とする請求項1、2又は3のいずれ
かに記載のデジタル放射線画像の撮影体位判別装置。 - (5)前記レベル分布検出手段が、レベル分布を検出す
る前記画像所定方向に沿った画素列それぞれの信号レベ
ルをそのまま用いて信号レベルの分布を検出することを
特徴とする請求項1、2又は3のいずれかに記載のデジ
タル放射線画像の撮影体位判別装置。 - (6)デジタル放射線画像の画素数と信号の量子化ステ
ップ数との少なくとも一方を縮減させる信号間引き手段
を備え、前記レベル分布検出手段が該信号間引き手段で
画素数と信号の量子化ステップ数との少なくとも一方が
縮減されたデジタル放射線画像信号に基づいてレベル分
布を検出することを特徴とする請求項1、2、3、4又
は5のいずれかに記載のデジタル放射線画像の撮影体位
判別装置。 - (7)デジタル放射線画像の濃度ヒストグラムの分布が
広がるようにコントラストを強調する補正をデジタル放
射線画像信号に施すコントラスト強調補正手段を備え、
前記レベル分布検出手段が前記コントラスト強調補正手
段で補正されたデジタル放射線画像信号に基づいてレベ
ル分布を検出することを特徴とする請求項1、2、3、
4、5又は6のいずれかに記載のデジタル放射線画像の
撮影体位判別装置。 - (8)前記レベル分布検出手段で検出されたレベル分布
を平滑化処理するレベル分布平滑化手段を備え、前記大
小変化パターン検出手段が前記レベル分布平滑化手段で
平滑化されたレベル分布に基づいて大小の変化パターン
を求めるようにしたことを特徴とする請求項1、2、3
、4、5、6又は7のいずれかに記載のデジタル放射線
画像の撮影体位判別装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2168537A JP2929030B2 (ja) | 1990-06-28 | 1990-06-28 | デジタル放射線画像の撮影体位判別装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2168537A JP2929030B2 (ja) | 1990-06-28 | 1990-06-28 | デジタル放射線画像の撮影体位判別装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0458941A true JPH0458941A (ja) | 1992-02-25 |
| JP2929030B2 JP2929030B2 (ja) | 1999-08-03 |
Family
ID=15869859
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2168537A Expired - Fee Related JP2929030B2 (ja) | 1990-06-28 | 1990-06-28 | デジタル放射線画像の撮影体位判別装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2929030B2 (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0486820A (ja) * | 1990-07-31 | 1992-03-19 | Fuji Photo Film Co Ltd | 放射線画像読取条件及び/又は画像処理条件決定方法 |
| JP2010082296A (ja) * | 2008-10-01 | 2010-04-15 | Infocom Corp | 胸部x線画像の変換方法、装置、システム、コンピュータプログラムおよび媒体 |
-
1990
- 1990-06-28 JP JP2168537A patent/JP2929030B2/ja not_active Expired - Fee Related
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0486820A (ja) * | 1990-07-31 | 1992-03-19 | Fuji Photo Film Co Ltd | 放射線画像読取条件及び/又は画像処理条件決定方法 |
| JP2010082296A (ja) * | 2008-10-01 | 2010-04-15 | Infocom Corp | 胸部x線画像の変換方法、装置、システム、コンピュータプログラムおよび媒体 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2929030B2 (ja) | 1999-08-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6775399B1 (en) | ROI segmentation image processing system | |
| US11672471B2 (en) | Osteoporosis diagnostic support apparatus | |
| JP2005296605A (ja) | 放射線写真画像を診断関連領域と診断非関連領域とにセグメント化する方法 | |
| JPH08335271A (ja) | 画像診断支援装置の画像処理方法 | |
| US6608915B2 (en) | Image processing method and apparatus | |
| US5283736A (en) | Radiographic image processing apparatus | |
| Zhang et al. | Automatic background recognition and removal (ABRR) in computed radiography images | |
| JP2001076141A (ja) | 画像認識方法および画像処理装置 | |
| US5933518A (en) | Method and device for image processing for automatically detecting objects in digitized images | |
| US20050069187A1 (en) | Image processing method, image processing apparatus and image processing program | |
| US6985614B2 (en) | Image processing method, apparatus and program | |
| US7418122B2 (en) | Image processing apparatus and method | |
| US20050161617A1 (en) | Image processing method, apparatus, and program | |
| JP2001148787A (ja) | 画像処理装置 | |
| JPH0458941A (ja) | デジタル放射線画像の撮影体位判別装置 | |
| JP4129598B2 (ja) | 画像処理装置及び画像処理方法 | |
| JP3731400B2 (ja) | 画像処理方法および画像処理装置 | |
| US6714674B1 (en) | Method for converting digital image pixel values | |
| JP3765320B2 (ja) | 放射線画像の画像処理条件決定装置及び画像処理装置 | |
| JP2929031B2 (ja) | デジタル放射線画像の撮影体位判別装置 | |
| JP2005167773A (ja) | 画像処理方法及び装置 | |
| JPH04341247A (ja) | 胸部放射線画像の撮影体位判別装置 | |
| JP2004086598A (ja) | 診断支援装置 | |
| JP4699166B2 (ja) | 画像処理装置、画像処理方法およびそのプログラム | |
| JP2001056854A (ja) | 医用ディジタル画像のダイナミック・レンジを決定するための方法および装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| 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: 20080521 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090521 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100521 Year of fee payment: 11 |
|
| LAPS | Cancellation because of no payment of annual fees |