JPH07100060B2 - 肺野部輪郭線を求める画像計測方法 - Google Patents
肺野部輪郭線を求める画像計測方法Info
- Publication number
- JPH07100060B2 JPH07100060B2 JP62077353A JP7735387A JPH07100060B2 JP H07100060 B2 JPH07100060 B2 JP H07100060B2 JP 62077353 A JP62077353 A JP 62077353A JP 7735387 A JP7735387 A JP 7735387A JP H07100060 B2 JPH07100060 B2 JP H07100060B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- contour
- point
- column
- density
- 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.)
- Expired - Fee Related
Links
- 210000004072 lung Anatomy 0.000 title claims description 25
- 238000000691 measurement method Methods 0.000 title 1
- 238000000034 method Methods 0.000 claims description 17
- 238000005259 measurement Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 5
- 230000005856 abnormality Effects 0.000 description 12
- 230000002159 abnormal effect Effects 0.000 description 9
- 238000011976 chest X-ray Methods 0.000 description 9
- 238000012545 processing Methods 0.000 description 6
- 238000000605 extraction Methods 0.000 description 5
- 238000010191 image analysis Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000002834 transmittance Methods 0.000 description 3
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 206010007572 Cardiac hypertrophy Diseases 0.000 description 1
- 208000006029 Cardiomegaly Diseases 0.000 description 1
- 206010020880 Hypertrophy Diseases 0.000 description 1
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004040 coloring Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000000968 intestinal effect Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Description
【発明の詳細な説明】 〔産業上の利用分野〕 本発明は、胸部放射線画像データから肺野部の輪郭線を
求める計測方法に関する。
求める計測方法に関する。
従来から、胸部X線画像の画像処理により肺野の輪郭線
を抽出し、得られた輪郭線の形状から被写体に異常があ
るか否かを判定する方法が提案されている(例えば特開
昭62−26047号公報)。このような処理を行う場合の輪
郭線の抽出は、画像の濃度データを適当な閾値を用いて
2値化し、『1』と『0』の境界を追跡して輪郭線とし
ている。或いは別に、現画像の濃度データに対しラプラ
シアンやエッジ抽出のオペレータを用いたフィルタリン
グ処理を施してエッジ画像を作成し、得られたエッジ画
像に対してノイズや不要な部分のエッジを除去した後
に、細線化処理を行うことにより輪郭線画像を求めてい
る ところが、前者の輪郭線抽出方法は、胸部X線画像では
1つの部分であっても位置によって濃度レベルが異なる
ために、肺野とその他の部分とを分けるような完全な閾
値を求めることが困難乃至不可能である。また、後者の
方法では、原画像の種類、濃度レベル、画素の大きさ、
輪郭線を求める対象となる部位、その他の細かい条件に
より、オペレータの値等のフィルタリングの条件が大き
く異なり、最適の条件を求めるにはかなりの専門的知識
が必要とされる。また、人体のX線画像のように各種の
部位が複雑に重なっている場合には、単一のフィルタリ
ングのみで全ての輪郭線を抽出することは難しく、数種
類のフィルタリングを組み合わせる等の複雑な処理が必
要となってくる。
を抽出し、得られた輪郭線の形状から被写体に異常があ
るか否かを判定する方法が提案されている(例えば特開
昭62−26047号公報)。このような処理を行う場合の輪
郭線の抽出は、画像の濃度データを適当な閾値を用いて
2値化し、『1』と『0』の境界を追跡して輪郭線とし
ている。或いは別に、現画像の濃度データに対しラプラ
シアンやエッジ抽出のオペレータを用いたフィルタリン
グ処理を施してエッジ画像を作成し、得られたエッジ画
像に対してノイズや不要な部分のエッジを除去した後
に、細線化処理を行うことにより輪郭線画像を求めてい
る ところが、前者の輪郭線抽出方法は、胸部X線画像では
1つの部分であっても位置によって濃度レベルが異なる
ために、肺野とその他の部分とを分けるような完全な閾
値を求めることが困難乃至不可能である。また、後者の
方法では、原画像の種類、濃度レベル、画素の大きさ、
輪郭線を求める対象となる部位、その他の細かい条件に
より、オペレータの値等のフィルタリングの条件が大き
く異なり、最適の条件を求めるにはかなりの専門的知識
が必要とされる。また、人体のX線画像のように各種の
部位が複雑に重なっている場合には、単一のフィルタリ
ングのみで全ての輪郭線を抽出することは難しく、数種
類のフィルタリングを組み合わせる等の複雑な処理が必
要となってくる。
本発明の目的は、人体の胸部放射線画像のような複雑な
画像の中の特定の部位の輪郭線を容易にしかも正確に抽
出する方法を提供することである。
画像の中の特定の部位の輪郭線を容易にしかも正確に抽
出する方法を提供することである。
このために本発明は、行列で表現された胸部放射線画像
データを対象として、1つの行或いは列についてのみ注
目し、その1次元の全部の濃度データの統計量を計算し
ていくつかの濃度閾値を求め、該求めた濃度閾値いの範
囲内の値を持つ濃度データの中で、前後の濃度データの
値の関係が予め定めた特定のパターンとなる点をその行
或いは列における輪郭の点とし、必要な範囲の行或いは
列について該輪郭の点を求めてそれらの点を結んだ線を
肺野の輪郭線とするようにした。
データを対象として、1つの行或いは列についてのみ注
目し、その1次元の全部の濃度データの統計量を計算し
ていくつかの濃度閾値を求め、該求めた濃度閾値いの範
囲内の値を持つ濃度データの中で、前後の濃度データの
値の関係が予め定めた特定のパターンとなる点をその行
或いは列における輪郭の点とし、必要な範囲の行或いは
列について該輪郭の点を求めてそれらの点を結んだ線を
肺野の輪郭線とするようにした。
以下、本発明の実施例について説明する。第1図は本発
明を実施するための装置の構成ブロックを示す図であ
る。1は被検体のX線写真画像を読み取るための画像入
力装置、2は読み取られたX線画像の濃度データを記憶
する画像記憶装置、3は画像記憶装置に記憶された濃度
データを読み出して輪郭線の抽出処理を行う画像解析装
置、4は解析された結果から被検体の異常を判定する判
定装置、5は判定された結果を使用者に報知する表示装
置である。
明を実施するための装置の構成ブロックを示す図であ
る。1は被検体のX線写真画像を読み取るための画像入
力装置、2は読み取られたX線画像の濃度データを記憶
する画像記憶装置、3は画像記憶装置に記憶された濃度
データを読み出して輪郭線の抽出処理を行う画像解析装
置、4は解析された結果から被検体の異常を判定する判
定装置、5は判定された結果を使用者に報知する表示装
置である。
まず、画像入力装置1は、例えばポリゴンミラー等の光
偏向器を用いてレーザ光等を使用した光ビームのスポッ
トで胸部X線画像面(フィムル)をスキャンし、その時
の透過光の強度をフォトセンサ等の光検出器で光電変換
し、得られた電気信号を順次A/D変換してデジタル画像
の濃度データとして出力する装置である。ここでは光ビ
ームを機械的にスキャンしているがCCD等の固体撮像素
子や撮像管を使用して電子的にスキャンする方法を採用
することもできる。また、CCDラインセンサ等を用いれ
ば、複数の胸部X線画像フィルムを順次機械的に搬送し
ながら連続的に読み取ることも可能である。
偏向器を用いてレーザ光等を使用した光ビームのスポッ
トで胸部X線画像面(フィムル)をスキャンし、その時
の透過光の強度をフォトセンサ等の光検出器で光電変換
し、得られた電気信号を順次A/D変換してデジタル画像
の濃度データとして出力する装置である。ここでは光ビ
ームを機械的にスキャンしているがCCD等の固体撮像素
子や撮像管を使用して電子的にスキャンする方法を採用
することもできる。また、CCDラインセンサ等を用いれ
ば、複数の胸部X線画像フィルムを順次機械的に搬送し
ながら連続的に読み取ることも可能である。
画像記憶装置2は、処理の高速性から通常はRAM等の半
導体メモリが使用されるが、磁気ディスク、フロッピー
ディスク、或いは光ディスク等のディスクメモリ、又は
磁気テープ、磁気バブルメモリ等を用いることもでき
る。
導体メモリが使用されるが、磁気ディスク、フロッピー
ディスク、或いは光ディスク等のディスクメモリ、又は
磁気テープ、磁気バブルメモリ等を用いることもでき
る。
画像解析装置3は本発明の要部を構成する部分である。
この装置3は画像記憶装置2内に記憶された画像の濃度
データを順次読み出して、肺野に相当する部分と周辺と
の輪郭線を示すデータを求める処理を行う。この装置3
においては、読み出された画像濃度データの少なくとも
1行分或いは1列分を記憶できるだけのラインメモリ、
解析を行うプロセッサ、抽出された輪郭線データを記憶
する輪郭線データメモリ、解析のアルゴリズム等を記憶
し、また解析時の作業領域となる主記憶装置、補助記憶
装置等からなる。ここでは、画像濃度データや輪郭線デ
ータを記憶するメモリを持つ構成にせずに、単一の主記
憶装置で全てを処理することもできる。
この装置3は画像記憶装置2内に記憶された画像の濃度
データを順次読み出して、肺野に相当する部分と周辺と
の輪郭線を示すデータを求める処理を行う。この装置3
においては、読み出された画像濃度データの少なくとも
1行分或いは1列分を記憶できるだけのラインメモリ、
解析を行うプロセッサ、抽出された輪郭線データを記憶
する輪郭線データメモリ、解析のアルゴリズム等を記憶
し、また解析時の作業領域となる主記憶装置、補助記憶
装置等からなる。ここでは、画像濃度データや輪郭線デ
ータを記憶するメモリを持つ構成にせずに、単一の主記
憶装置で全てを処理することもできる。
判定装置4は、画像解析装置3で得られた輪郭線データ
の各種の特徴量を計算して、検診の項目毎に標準値の範
囲から出たか否かを判定する。
の各種の特徴量を計算して、検診の項目毎に標準値の範
囲から出たか否かを判定する。
表示装置5は判定装置4の出力によって、解析の行った
画像データの被検体に異常があることを医師に知らせ
る。
画像データの被検体に異常があることを医師に知らせ
る。
次に作用を説明する。
(X線写真画像の入力) X線画像は画像入力装置1によって約200μmのサンプ
リングピッチで読み取られ、濃度範囲0〜3を256階調
に量子化される。このときX線画像は正位して入力さ
れ、第2図に示すように画像の左上が(1,1)、右下が
(m,n)としたm行n列の行列形式により、画像記憶装
置2に格納される。
リングピッチで読み取られ、濃度範囲0〜3を256階調
に量子化される。このときX線画像は正位して入力さ
れ、第2図に示すように画像の左上が(1,1)、右下が
(m,n)としたm行n列の行列形式により、画像記憶装
置2に格納される。
なお、サンプリングピッチは200μmに固定されるもの
ではなく、場合に応じて25μmから1000μmの範囲で変
化させることができる。また読取濃度範囲や量子化レベ
ル数も変化させることができる。
ではなく、場合に応じて25μmから1000μmの範囲で変
化させることができる。また読取濃度範囲や量子化レベ
ル数も変化させることができる。
(肺野輪郭線の抽出) 肺野の輪郭線の抽出は、第2図に示すように、左右各肺
野のそれぞれについて外側輪郭線AL、AR、内側輪郭線B
L、BR、横隔膜輪郭線CL、CRについて行う。外側及び内
側輪郭線AL、AR、BL、BRについては、画像データの各行
における輪郭線の存在する列の値の配列として求めら
れ、横隔膜輪郭線CL、CRについては、各列における輪郭
線の存在する行の配列として求められる。
野のそれぞれについて外側輪郭線AL、AR、内側輪郭線B
L、BR、横隔膜輪郭線CL、CRについて行う。外側及び内
側輪郭線AL、AR、BL、BRについては、画像データの各行
における輪郭線の存在する列の値の配列として求めら
れ、横隔膜輪郭線CL、CRについては、各列における輪郭
線の存在する行の配列として求められる。
肺野の外側輪郭線AL、ARと内側輪郭線BL、BRの抽出は、
画像記憶装置2より横方向の1行分の画像濃度データを
画像解析装置3内に読み込み、読み込んだ全データの統
計量を計算して、それらの値からいくつかの濃度閾値を
求める。そして、それらの閾値の範囲内の値を持つ画像
濃度データ(点)の中で、ある点について前後の点との
値の関係が後記する特定のパターン(条件)を満たす場
合に、そのある点の位置をその行における肺野の輪郭の
列の位置として記憶する。特定のパターンとしては、 極小となる点…例えば第3図(b)のa 傾きが最小となる点…例えば第3図(b)のc 傾きが最大となる点…例えば第3図(b)のb 等があり、これらのパターンを満足する点をいくつか探
し出す。第3図は(a)で示す行m1〜m4の列方向のライ
ンの濃度変化を(b)に示す。Dは濃度(透過率)であ
る。更に、 解析の対象となっている行の画像上の位置 求められた点の位置 求められた点の濃度値 求められた点における値の傾き 用いた閾値の範囲 等を総合的に判断し、いずれの輪郭の点であるか、或い
はノイズであるかを決定する。例えば、右肺野外側の輪
郭ARを構成する点であれば、『画像の向かって左側1/3
の範囲に存在し、全濃度データの平均値以下の値(濃度
を透過率で示した場合、高濃度部は小さい透過率を示
す。)をもち、前後の点との値の関係は極小或いは傾き
が0となる点』といった条件を満足する点を探索する。
画像記憶装置2より横方向の1行分の画像濃度データを
画像解析装置3内に読み込み、読み込んだ全データの統
計量を計算して、それらの値からいくつかの濃度閾値を
求める。そして、それらの閾値の範囲内の値を持つ画像
濃度データ(点)の中で、ある点について前後の点との
値の関係が後記する特定のパターン(条件)を満たす場
合に、そのある点の位置をその行における肺野の輪郭の
列の位置として記憶する。特定のパターンとしては、 極小となる点…例えば第3図(b)のa 傾きが最小となる点…例えば第3図(b)のc 傾きが最大となる点…例えば第3図(b)のb 等があり、これらのパターンを満足する点をいくつか探
し出す。第3図は(a)で示す行m1〜m4の列方向のライ
ンの濃度変化を(b)に示す。Dは濃度(透過率)であ
る。更に、 解析の対象となっている行の画像上の位置 求められた点の位置 求められた点の濃度値 求められた点における値の傾き 用いた閾値の範囲 等を総合的に判断し、いずれの輪郭の点であるか、或い
はノイズであるかを決定する。例えば、右肺野外側の輪
郭ARを構成する点であれば、『画像の向かって左側1/3
の範囲に存在し、全濃度データの平均値以下の値(濃度
を透過率で示した場合、高濃度部は小さい透過率を示
す。)をもち、前後の点との値の関係は極小或いは傾き
が0となる点』といった条件を満足する点を探索する。
このような探索のアルゴリズムは、従来からの手続き型
のプログラム言語、例えばFORTRANやPascal等によって
記述できることもできるが、最近のルール指向言語、例
えばOPS5、OPS83等によって更に容易に記述することが
できる。
のプログラム言語、例えばFORTRANやPascal等によって
記述できることもできるが、最近のルール指向言語、例
えばOPS5、OPS83等によって更に容易に記述することが
できる。
このようにして画像データの各行についてそれぞれの輪
郭を構成する点の位置を求め、必要な部分について結び
合わせることにより、それぞれの輪郭線が得られる。
郭を構成する点の位置を求め、必要な部分について結び
合わせることにより、それぞれの輪郭線が得られる。
横隔膜の輪郭線については、画像読取装置1より縦方向
mの1列分の画像データを画像解析装置3のラインメモ
リに読み込み、上述した方法と同様にして列における輪
郭の存在する行の位置として輪郭線を得る。第4図は
(a)で示す列n1〜n5の行方向のラインの濃度変化を
(b)に示す。探索の条件としては、主にある点から前
後に数点離れた2つの点の間の値の差が最大となるその
ある点を用いる。なお、第4図では傾きが最大となる点
を用いている(第4図(b)のd参照)。
mの1列分の画像データを画像解析装置3のラインメモ
リに読み込み、上述した方法と同様にして列における輪
郭の存在する行の位置として輪郭線を得る。第4図は
(a)で示す列n1〜n5の行方向のラインの濃度変化を
(b)に示す。探索の条件としては、主にある点から前
後に数点離れた2つの点の間の値の差が最大となるその
ある点を用いる。なお、第4図では傾きが最大となる点
を用いている(第4図(b)のd参照)。
このようにして、左右各肺野について外側AL、AR、内側
BL、BR、横隔膜CL、CRの輪郭線を求め、画像解析装置3
内の輪郭線データメモリに格納する。輪郭線データは、
外側、内側のものについては、行の位置に対応した番地
のメモリ上に列の位置の値として記憶される。横隔膜輪
郭線データは列の位置に対応した番地のメモリ上に行の
位置の値として記憶される。なお、得られた輪郭線のそ
れぞれの点の行と列の値の組を、順次記憶させる方法を
用いることもできる。
BL、BR、横隔膜CL、CRの輪郭線を求め、画像解析装置3
内の輪郭線データメモリに格納する。輪郭線データは、
外側、内側のものについては、行の位置に対応した番地
のメモリ上に列の位置の値として記憶される。横隔膜輪
郭線データは列の位置に対応した番地のメモリ上に行の
位置の値として記憶される。なお、得られた輪郭線のそ
れぞれの点の行と列の値の組を、順次記憶させる方法を
用いることもできる。
以上のように、本実施例では、行列で表現された胸部X
線画像データを対象として、1つの行或いは列について
のみ注目し、その1次元の濃度データ列の中で前後のデ
ータの値の関係が予め定めた特定のパターンとなる点を
その行或いは列における輪郭の点とし、必要な範囲の行
或いは列について該輪郭の点を求めてそれらの点を結ん
だ線を肺野の輪郭線としている。
線画像データを対象として、1つの行或いは列について
のみ注目し、その1次元の濃度データ列の中で前後のデ
ータの値の関係が予め定めた特定のパターンとなる点を
その行或いは列における輪郭の点とし、必要な範囲の行
或いは列について該輪郭の点を求めてそれらの点を結ん
だ線を肺野の輪郭線としている。
そして、上記特定のパターンとなる点として、所定の濃
度閾値の範囲内にあり且つある点から所定の距離だけ離
れた2点間の値の差が最大となるそのある点、或いは所
定の濃度閾値の範囲内にある極小の点その他を使用する
ことができる。
度閾値の範囲内にあり且つある点から所定の距離だけ離
れた2点間の値の差が最大となるそのある点、或いは所
定の濃度閾値の範囲内にある極小の点その他を使用する
ことができる。
このような特定のパターンを複数用意して、計測を行っ
ている行或いは列の位置によって、又は目的とする輪郭
の部位によって、上記複数内から1つのパターンを選択
すること、又は特定のパターンを複数用意して行或いは
列における輪郭の点の候補を複数求め、該候補点の位置
よって、或いは計測を行っている行或いは列の位置によ
って、更には目的とする輪郭の部位によって、輪郭の点
を1つ選択することができる。
ている行或いは列の位置によって、又は目的とする輪郭
の部位によって、上記複数内から1つのパターンを選択
すること、又は特定のパターンを複数用意して行或いは
列における輪郭の点の候補を複数求め、該候補点の位置
よって、或いは計測を行っている行或いは列の位置によ
って、更には目的とする輪郭の部位によって、輪郭の点
を1つ選択することができる。
(輪郭線の形状による異常の判定) 肺野の輪郭線の形状を計測することにより、下記の異常
を検出することができる。
を検出することができる。
心臓の肥大…肺野の内側輪郭線BL、BRの心臓部分の外
への張り出し具合から判定する。
への張り出し具合から判定する。
大動脈の突出…左肺野内側輪郭線BLの上部の外への張
り出し具合から判定する。
り出し具合から判定する。
横隔膜の異常…横隔膜輪郭線CL、CRの滑らかさから判
定する。
定する。
また、肺野の輪郭線の情報を利用することにより、肺
野の異常陰影、肺尖部の異常、脊椎の湾曲等の異常
の検出の探索範囲を限定することができる。
野の異常陰影、肺尖部の異常、脊椎の湾曲等の異常
の検出の探索範囲を限定することができる。
(判定結果の表示) 以上のようにして求めた判定結果は、この支援装置を用
いている医師が判断できる形で表示される。ここでは、
操作を指示するためのCRT端末に、異常があった場合に
その胸部X線写真番号、被検体のデータ及び異常のあっ
た項目と計測した値とを表示する。ここで、表示方法と
しては、印字式の端末を用いることができる。表示項目
は変更することができる。
いている医師が判断できる形で表示される。ここでは、
操作を指示するためのCRT端末に、異常があった場合に
その胸部X線写真番号、被検体のデータ及び異常のあっ
た項目と計測した値とを表示する。ここで、表示方法と
しては、印字式の端末を用いることができる。表示項目
は変更することができる。
また、計測の対象となった画像データを画像表示端末に
表示し、異常と判定された部分に色をつけたり輝度を増
す等の表示を行うことも可能であり、更に視覚的な表示
だけではく、異常の際に警告音を発生させることもでき
る。また、胸部X線写真フィルムを自動的に搬送する機
構を持つ装置では、異常の見られたフィルムのみを正常
のフィルムとは別の搬送口から排出することにより、異
常を知らせることもできる。
表示し、異常と判定された部分に色をつけたり輝度を増
す等の表示を行うことも可能であり、更に視覚的な表示
だけではく、異常の際に警告音を発生させることもでき
る。また、胸部X線写真フィルムを自動的に搬送する機
構を持つ装置では、異常の見られたフィルムのみを正常
のフィルムとは別の搬送口から排出することにより、異
常を知らせることもできる。
なお、以上の説明では、X線写真のみを読み取って判定
を行う場合を述べたが、X線TVや輝尽性螢光体を用いた
X線撮影装置のように、撮影時に直接電気信号に変換さ
れる装置を用いる場合には、この信号を直接A/D変換し
て画像記憶装置2格納するとこができる。またX線写真
に限らずCT等の断層写真やMRI等の他のモダリティの画
像を用いることも可能である。更に、サブストラクショ
ン等の画像処理を行ない、骨の陰影を除去した画像等を
用いることも効果的である。
を行う場合を述べたが、X線TVや輝尽性螢光体を用いた
X線撮影装置のように、撮影時に直接電気信号に変換さ
れる装置を用いる場合には、この信号を直接A/D変換し
て画像記憶装置2格納するとこができる。またX線写真
に限らずCT等の断層写真やMRI等の他のモダリティの画
像を用いることも可能である。更に、サブストラクショ
ン等の画像処理を行ない、骨の陰影を除去した画像等を
用いることも効果的である。
(実施例の検証) 以上のように構成した装置を用いて、医師から提供され
た100例の胸部X線写真について肺野輪郭線の抽出の結
果を調べたところ、7例の写真を除いては、ほぼ良好な
輪郭であるとの評価を医師から得ることができた。7例
については、画像の濃度レベルの異常、腸のガスの像等
の影響、疾患による画像のボケ等の原因よっていた。ま
た、得られた輪郭線情報を用いて心臓肥大等の計測を行
ったところ、95%の異常例について、計測した値が標準
値の範囲を越え、異常と判定できる結果が得られた。
た100例の胸部X線写真について肺野輪郭線の抽出の結
果を調べたところ、7例の写真を除いては、ほぼ良好な
輪郭であるとの評価を医師から得ることができた。7例
については、画像の濃度レベルの異常、腸のガスの像等
の影響、疾患による画像のボケ等の原因よっていた。ま
た、得られた輪郭線情報を用いて心臓肥大等の計測を行
ったところ、95%の異常例について、計測した値が標準
値の範囲を越え、異常と判定できる結果が得られた。
以上のように本発明を用いることにより、胸部X線写真
の肺野像のように複雑な部分が重なり合った画像データ
から、肺野部分等の輪郭線を容易且つ正確に抽出するこ
とができる。そして、従来から行われていた胸部X写真
について画像計測を行ない被検体に異常があるか否か
を、医師が読影を行う前に、表示することができる。
の肺野像のように複雑な部分が重なり合った画像データ
から、肺野部分等の輪郭線を容易且つ正確に抽出するこ
とができる。そして、従来から行われていた胸部X写真
について画像計測を行ない被検体に異常があるか否か
を、医師が読影を行う前に、表示することができる。
従って、異常の疑いのある写真のみについて細かく観察
し、疑いのない写真については確認のみを行うことで済
ますことができる。よって、従来のX線写真を1枚1枚
丁寧に観察することによって行っていたX線写真検診の
能率を向上させることができる。
し、疑いのない写真については確認のみを行うことで済
ますことができる。よって、従来のX線写真を1枚1枚
丁寧に観察することによって行っていたX線写真検診の
能率を向上させることができる。
また、被検体や撮影・現像の条件により写真の仕上がり
濃度にバラツキがあっても、画像処理により一定の条件
での計測が可能となるので、肉眼では読影が困難な写真
であっても異常の有無を判定でき、疾病の見落しを減少
させることができるようになる。
濃度にバラツキがあっても、画像処理により一定の条件
での計測が可能となるので、肉眼では読影が困難な写真
であっても異常の有無を判定でき、疾病の見落しを減少
させることができるようになる。
第1図は本発明の方法を実施するための装置のブロック
図、第2図は画像データの行列による表現の説明図、第
3図は特定の行の列方向の濃度データを調べて肺野の内
外輪郭線を構成する点の抽出を行う説明図、第4図は特
定の列の行方向の濃度データを調べて肺野の横隔膜輪郭
線を構成する点の抽出を行う説明図である。
図、第2図は画像データの行列による表現の説明図、第
3図は特定の行の列方向の濃度データを調べて肺野の内
外輪郭線を構成する点の抽出を行う説明図、第4図は特
定の列の行方向の濃度データを調べて肺野の横隔膜輪郭
線を構成する点の抽出を行う説明図である。
Claims (5)
- 【請求項1】行列で表現された胸部放射線画像データを
対象として、1つの行或いは列についてのみ注目し、そ
の1次元の全部の濃度データの統計量を計算していくつ
かの濃度閾値を求め、該求めた濃度閾値の範囲内の値を
持つ濃度データの中で、前後の濃度データの値との関係
が予め定めた特定のパターンとなる点をその行或いは列
における輪郭の点とし、必要な範囲の行或いは列につい
て上記輪郭の点を求めてそれらの点を結んだ線を肺野の
輪郭線とすることを特徴とする画像計測方法 - 【請求項2】上記特定のパターンとなる点を、所定の距
離だけ離れた2点間の濃度値の差が最大となる点とした
ことを特徴とする特許請求の範囲第1項記載の画像計測
方法。 - 【請求項3】上記特定のパターンとなる点を、極小の点
としたことを特徴とする特許請求の範囲第1項記載の画
像計測方法。 - 【請求項4】上記特定のパターンを複数用意して、計測
を行なっている行或いは列の位置によって、又は目的と
する輪郭の部位によって、上記複数内から1つのパター
ンを選択するようにしたことを特徴とする特許請求の範
囲第1項乃至第3項記載の画像計測方法。 - 【請求項5】上記特定のパターンを複数用意して計測を
行なっている行或いは列における輪郭の点の候補を複数
求め、該候補点の位置によって、或いは計測を行なって
いる行或いは列の位置によって、更には目的とする輪郭
の部位によって、輪郭の点を1つ選択するようにしたこ
とを特徴とする特許請求の範囲第1項乃至第3項記載の
画像計測方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP62077353A JPH07100060B2 (ja) | 1987-03-30 | 1987-03-30 | 肺野部輪郭線を求める画像計測方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP62077353A JPH07100060B2 (ja) | 1987-03-30 | 1987-03-30 | 肺野部輪郭線を求める画像計測方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPS63240832A JPS63240832A (ja) | 1988-10-06 |
JPH07100060B2 true JPH07100060B2 (ja) | 1995-11-01 |
Family
ID=13631547
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP62077353A Expired - Fee Related JPH07100060B2 (ja) | 1987-03-30 | 1987-03-30 | 肺野部輪郭線を求める画像計測方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPH07100060B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005176896A (ja) * | 2003-12-16 | 2005-07-07 | Canon Inc | X線画像処理装置、x線画像処理方法、プログラム及びコンピュータ可読記憶媒体 |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3631095B2 (ja) | 2000-04-17 | 2005-03-23 | キヤノン株式会社 | 照射野領域抽出装置、放射線撮影装置、放射線画像用システム、照射野領域抽出方法、及びコンピュータ可読記憶媒体 |
JP4660998B2 (ja) * | 2001-07-31 | 2011-03-30 | トヨタ自動車株式会社 | 接合検査装置 |
JP4639338B2 (ja) * | 2007-03-20 | 2011-02-23 | 国立大学法人 大分大学 | 胸部x線画像からのリブケイジ境界検出方法 |
CA3031175C (en) | 2016-07-19 | 2021-11-09 | Radwisp Pte.Ltd. | Diagnosis support program |
WO2019135412A1 (ja) * | 2018-01-05 | 2019-07-11 | ラドウィスプ プライベート リミテッド | 診断支援プログラム |
CN116452591B (zh) * | 2023-06-16 | 2023-08-15 | 天津市肿瘤医院(天津医科大学肿瘤医院) | 一种基于超声引导的胆道造影图像位置标记方法 |
-
1987
- 1987-03-30 JP JP62077353A patent/JPH07100060B2/ja not_active Expired - Fee Related
Non-Patent Citations (1)
Title |
---|
「医用画像処理」(1982−2−25)朝倉書店P.115−118 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005176896A (ja) * | 2003-12-16 | 2005-07-07 | Canon Inc | X線画像処理装置、x線画像処理方法、プログラム及びコンピュータ可読記憶媒体 |
Also Published As
Publication number | Publication date |
---|---|
JPS63240832A (ja) | 1988-10-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3137561B2 (ja) | 画像品質を判定する方法及び装置、画像捕捉装置の性能を監視する方法及び装置 | |
JP4604451B2 (ja) | 医用画像処理装置及び悪性度の判定方法 | |
JPH11510419A (ja) | 定量的の歯科カリエス検査システムおよびその方法 | |
US10540765B2 (en) | Image processing device, image processing method, and computer program product thereon | |
JP4915737B2 (ja) | 画像解析システム、及び画像解析プログラム | |
JP4010375B2 (ja) | 放射線画像処理方法 | |
US5224036A (en) | Pattern recognition apparatus | |
JP2014221163A (ja) | 画像診断支援装置及び支援方法 | |
US3869602A (en) | Apparatus for automatic computation of cardiothoracic ratio | |
JPH06304159A (ja) | 画像の撮影体位認識方法 | |
JP4935895B2 (ja) | エッジ評価方法とエッジ検出方法と画像補正方法と画像処理システム | |
JPH07100060B2 (ja) | 肺野部輪郭線を求める画像計測方法 | |
US7248729B2 (en) | Image processing apparatus and image processing method | |
JP2001076141A (ja) | 画像認識方法および画像処理装置 | |
JPH09251535A (ja) | 2次元多階調画像におけるパタン検出方法及び装置 | |
JPH09106448A (ja) | 異常陰影候補の検出方法 | |
JPH057579A (ja) | 放射線画像の照射野領域抽出装置 | |
JP3901188B2 (ja) | 放射線画像処理方法 | |
JP3731400B2 (ja) | 画像処理方法および画像処理装置 | |
JPH1119077A (ja) | 放射線画像における腫瘤影の検出方法及び装置 | |
CN114581402A (zh) | 胶囊内窥镜质检方法、装置及存储介质 | |
JP2598636B2 (ja) | 画像処理条件決定方法 | |
JPH07100061B2 (ja) | 脊椎湾曲を求める画像計測方法 | |
JPH08215183A (ja) | 石灰化検出方法および装置 | |
JP2582667B2 (ja) | 線状パターン幅演算装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
LAPS | Cancellation because of no payment of annual fees |