JP2008073076A - 乳房x線画像処理装置 - Google Patents
乳房x線画像処理装置 Download PDFInfo
- Publication number
- JP2008073076A JP2008073076A JP2006252326A JP2006252326A JP2008073076A JP 2008073076 A JP2008073076 A JP 2008073076A JP 2006252326 A JP2006252326 A JP 2006252326A JP 2006252326 A JP2006252326 A JP 2006252326A JP 2008073076 A JP2008073076 A JP 2008073076A
- Authority
- JP
- Japan
- Prior art keywords
- line
- breast
- image data
- image
- pectoral muscle
- 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.)
- Pending
Links
- 210000000481 breast Anatomy 0.000 claims abstract description 136
- 238000012545 processing Methods 0.000 claims abstract description 71
- 210000002976 pectoralis muscle Anatomy 0.000 claims abstract description 48
- 238000000605 extraction Methods 0.000 claims abstract description 11
- 238000013519 translation Methods 0.000 claims abstract description 3
- 239000013598 vector Substances 0.000 claims description 15
- 210000000062 pectoralis major Anatomy 0.000 claims description 14
- 238000005314 correlation function Methods 0.000 claims 1
- 238000001514 detection method Methods 0.000 description 31
- 230000002159 abnormal effect Effects 0.000 description 23
- 238000000034 method Methods 0.000 description 19
- 210000002445 nipple Anatomy 0.000 description 14
- 210000001519 tissue Anatomy 0.000 description 13
- 238000005311 autocorrelation function Methods 0.000 description 12
- 238000003384 imaging method Methods 0.000 description 12
- 238000003745 diagnosis Methods 0.000 description 8
- 210000005075 mammary gland Anatomy 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 206010006187 Breast cancer Diseases 0.000 description 5
- 208000026310 Breast neoplasm Diseases 0.000 description 5
- 206010028980 Neoplasm Diseases 0.000 description 5
- 238000006073 displacement reaction Methods 0.000 description 5
- 210000000577 adipose tissue Anatomy 0.000 description 4
- 238000012217 deletion Methods 0.000 description 3
- 230000037430 deletion Effects 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 230000003902 lesion Effects 0.000 description 3
- 238000003825 pressing Methods 0.000 description 3
- 208000004434 Calcinosis Diseases 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000010191 image analysis Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 238000004195 computer-aided diagnosis Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000009607 mammography Methods 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 238000002559 palpation Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
Images
Abstract
【課題】左右の乳房のX線画像について、より正確に位置合わせ処理を行う乳房X線画像処理装置を提供する。
【解決手段】乳房X線画像処理装置は、乳房の内外斜位方向のX線画像の画像データを取り込む画像データ取込部と、画像データからスキンラインを検出し乳房領域を抽出する乳房領域抽出手段21、大胸筋ラインを検出し大胸筋領域を抽出する大胸筋領域抽出手段22、大胸筋ラインを直線近似し大胸筋直線を設定する大胸筋直線設定手段23、大胸筋直線とスキンライン上の画素点との最短距離を検出し、該最短距離が最長の画素点を基準点とする基準点設定手段25、左右のX線画像の基準点が中心線に対して線対称となるように画像データの座標を平行移動させる平行移動手段26、及び左右のX線画像の大胸筋直線が中心線に対して線対称となるように画像データの座標を基準点を中心に回転させる回転手段27を備える位置合わせ手段20とを具備する。
【選択図】図2
【解決手段】乳房X線画像処理装置は、乳房の内外斜位方向のX線画像の画像データを取り込む画像データ取込部と、画像データからスキンラインを検出し乳房領域を抽出する乳房領域抽出手段21、大胸筋ラインを検出し大胸筋領域を抽出する大胸筋領域抽出手段22、大胸筋ラインを直線近似し大胸筋直線を設定する大胸筋直線設定手段23、大胸筋直線とスキンライン上の画素点との最短距離を検出し、該最短距離が最長の画素点を基準点とする基準点設定手段25、左右のX線画像の基準点が中心線に対して線対称となるように画像データの座標を平行移動させる平行移動手段26、及び左右のX線画像の大胸筋直線が中心線に対して線対称となるように画像データの座標を基準点を中心に回転させる回転手段27を備える位置合わせ手段20とを具備する。
【選択図】図2
Description
本発明は、乳房X線画像処理装置に関するものであり、特に、内外斜位方向から撮影された乳房X線画像の処理を行う乳房X線画像処理装置に関するものである。
乳房X線写真(マンモグラム)による乳がんの診断は、触診では発見が困難な初期段階のがん病変を発見することができ、乳がんの早期発見のために有効な診断方法である。近年では、わが国における乳がんの羅患率が増加傾向にあることから、集団検診へのマンモグラムの導入が推進されている。
ここで、マンモグラムによる乳がんの診断は、医師がX線写真を観察し、微小石灰化クラスタや腫瘤による異常陰影を検出(読影)することにより行われる。そのため、集団検診へのマンモグラムの導入の拡大に伴い、撮影されるX線写真の枚数が増加し、読影する医師の負担が著しく増加することが懸念される。そこで、マンモグラムをコンピュータで自動的に解析し、診断に役立つ情報を医師に提示して、医師による診断を支援するコンピュータ支援診断(Computer−aided diagnosis,以下「CAD」という)システムの開発が進められている。例えば、既に確定的な診断がなされたX線画像のデータベースを検索し、診断の対象である画像と類似する画像を抽出して、医師に提示するCADシステムが提案されている。
ところで、乳房X線写真の読影に際しては、左右の乳房を対比観察する手法がとられることが多い。これは、正常な左右の乳房における組織はほぼ等しいため、一方の乳房で異常の疑いのある陰影が発見された場合、他方の乳房の同じ位置を観察して同様の陰影の有無を確認することが、異常陰影であるか正常な組織であるかの判断に役立つからである。
そこで、従来のCADシステムにおいても、左右一対の乳房のX線画像を対比して解析を行うために、左右のX線画像について乳房の位置を合わせる処理が試みられている。例えば、スキンラインを検出してその形状により位置合わせを行うもの(非特許文献1,非特許文献2参照)や、乳頭を検出し、これを基準として位置合わせを行うもの(非特許文献3,非特許文献4参照)が報告されている。また、従来のCADシステムでは、正常な組織が異常陰影として検出される「偽陽性候補」がまだまだ多いため、いかに偽陽性候補の削除を行うかが、CADシステムの実用化のための大きな課題となっている。そのため、これらの非特許文献においても、左右の乳房のX線画像に対してエッジ強調フィルタ処理を行い、対応する点の画素値の差を利用することによって偽陽性候補を削除する技術について、検討が行われている。
笠井 聡,藤田廣志,原 武史,畑中裕司,遠藤登喜子、「腫瘤陰影自動検出アルゴリズムにおける左右乳房画像の比較による偽陽性候補の削除」、Medical Imaging Technology、1998年、第16巻、第6号、p.655−666
二村 仁,畑中裕司,原 武史,藤田廣志,遠藤登喜子,岩瀬拓士、「マンモグラムCADシステムにおける腫瘤陰影の検出法」、1999年度医用画像工学研究会講演論文集、2000年、p.134−137
平井大輔,根本充貴,清水昭伸,小畑秀文,縄野 繁、「左右乳房X線像の特徴量比較による腫瘤影判別能の高度化」、信学技報、2003年、p.1−5
陸 苗,中川俊明,原 武史,藤田廣志,遠藤登喜子,岩瀬拓士、「ディジタルマンモグラムにおける乳頭領域の自動検出」、医用画像情報学会雑誌、2002年、第19巻、第3号、p.141−144
しかしながら、診断対象の画像と類似する画像をデータベースから抽出する方法は、画像解析の精度が、充実したデータベースが構築されているか否かに大きく依存するという問題があった。また、データベースが膨大であれば、より類似した画像を抽出できるものの、高速に処理できないという問題があった。
また、左右の乳房のX線画像を対比する方法では、データベースを必要とせず高速な処理が可能であるが、前提となる左右の乳房のX線画像の位置合わせが正確に行えないという問題があった。すなわち、従来の位置合わせ処理では、X線撮影時に左右両側から乳房を挟み込む位置が左右の乳房によって異なることによって、左右のX線画像間に生じる撮影された乳房位置の回転方向のずれは、補正することができないという問題があった。また、スキンラインの形状による位置合わせは(非特許文献1,非特許文献2参照)、一般的な日本人のように、スキンライン形状に特徴が乏しい小さな乳房の処理には適さないという問題があった。更に、スキンラインの先端を乳頭と定義する処理は(非特許文献3参照)、X線撮影時の乳房の押圧の仕方により、スキンラインの先端に乳頭が存在しないこともあるという問題があった。加えて、スキンラインを一次微分し、その変曲点を乳頭と定義する処理は(非特許文献4)、乳頭がスキンラインから突出して撮影されない場合には適用できないという問題があった。
また、正常な組織が異常陰影として誤って検出される可能性の高いものとして、乳腺組織により形成される索状の陰影が挙げられる。しかしながら、乳腺組織の微細な構造までは、左右の乳房でも完全に一致するものではないため、左右のX線画像の画素値の差に基づいて偽陽性候補を削除するという従来のシステムでは、偽陽性候補を充分に削除することができないという問題があった。
そこで、本発明は、上記の実情に鑑み、左右の乳房のX線画像について、より正確に位置合わせ処理を行う乳房X線画像処理装置の提供を、第一の課題とするものである。また、偽陽性候補を充分に削除することができる乳房X線画像処理装置の提供を、第二の課題とするものである。
上記の課題を解決するため、本発明にかかる乳房X線画像処理装置は、「乳房の内外斜位方向のX線画像の画像データを取り込む画像データ取込部と、取り込まれた前記画像データからスキンラインを検出し乳房領域を抽出する乳房領域抽出手段、取り込まれた前記画像データから大胸筋ラインを検出し大胸筋領域を抽出する大胸筋領域抽出手段、前記大胸筋ラインを直線に近似し大胸筋直線を設定する大胸筋直線設定手段、前記大胸筋直線と前記スキンライン上の全ての画素点との最短距離を検出し、検出された前記最短距離が最長の前記スキンライン上の前記画素点を基準点として設定する基準点設定手段、左右の乳房の前記X線画像についてぞれぞれ設定された前記基準点が、中心線に対して線対称となるように、左右の乳房の少なくとも何れか一方の画像データの座標を平行移動させる平行移動手段、及び、左右の乳房の前記X線画像についてぞれぞれ設定された前記大胸筋直線が、前記中心線に対して線対称となるように、左右の乳房の少なくとも何れか一方の画像データの座標を前記基準点を中心に回転させる回転手段を備え、左右の乳房の前記X線画像の位置合わせ処理を行う位置合わせ手段とを」具備している。
乳房のX線撮影には、乳房を左右両側から挟み込むように内外斜位方向(MLO;Medio−Lateral Oblique)で撮影する場合と、乳房を上下から挟み込むように頭尾方向(CC;Cranio−caudal)で撮影する場合があるが、本発明ではMLO方向で撮影された画像を対象とする。
「画像データ」は、フィルム/増感紙システムでX線撮影されたアナログ画像を、レーザスキャナ等のデジタル化装置によって電子データに変換したもの、或いは、CR(Computer−Radiology)やFPD(Flat Panel Detector)などのデジタル系撮像装置を用いた撮影により、直接電子データとして得られたものを、使用することができる。この画像データを乳房X線画像処理装置外から取り込む「画像データ取込部」としては、CD−Rやフレキシブルディスク等の記憶媒体に格納された画像データを読取るドライブ装置を例示することができる。或いは、レーザスキャナ等のデジタル化装置やデジタル系撮像装置から乳房X線画像処理装置まで、ハードウェアインターフェースを介して画像データを伝送することとし、ハードウェアインターフェースと接続されるコネクタ部を画像データ取込部とすることもできる。なお、デジタル化装置は、乳房X線画像処理装置の一構成であっても良い。
「位置あわせ処理」とは、左右一対の乳房のそれぞれをX線撮影する場合、撮影された画像全体における乳房の位置は通常同一とはならないため、対比による左右の乳房の画像解析を可能とするために、X線撮影後の画像データについて左右の画像における乳房の位置を対応させる処理である。
「乳房領域」は、被写体である乳房が撮影された領域であり、領域内には乳腺組織、脂肪組織が含まれる。そして「スキンライン」は、この乳房領域と、X線が直接透過して黒く表れる乳房外領域との境界となるラインである。また、「大胸筋領域」は乳腺組織や脂肪組織よりX線の吸収率が高い大胸筋組織の領域であり、右乳房の場合は画像の上端から右斜め下方向に向って、左乳房の場合は画像の上端から左斜め下方向に向って、白く表れる。そして、「大胸筋ライン」は、大胸筋領域と乳房領域との境界となるラインである。
「乳房領域抽出手段」及び「大胸筋領域抽出手段」は、何れも周囲との画素値の差異を利用してそれぞれの領域を抽出することができる。すなわち、X線画像上を走査するように画素値を測定し、互いに隣接する画素において画素値が大きく変化する箇所を認識することにより、スキンライン及び大胸筋ラインを検出でき、これと同時に、それぞれ乳房領域及び大胸筋領域を抽出することができる。なお、画素値の差異をより明確にするために、適宜の閾値を設定して二値化処理を行うこともできる。
「大胸筋直線決定手段」は、わずかに湾曲している大胸筋ラインを直線に近似する手段であり、曲線を直線に変換する公知のアルゴリズムを使用することができる。
「中心線」は、左右の乳房のX線画像を線対称に位置合せをするための中心となる線であり、データ処理上で仮想される線であっても良い。或いは、X線画像をモニタ等の表示装置に表示しながら位置合わせを行う場合、中心線が実際に画面に表示されても良い。
「位置合わせ手段」は、中央処理装置(CPU)や主記憶装置を具備して構成され、主記憶装置には、乳房領域抽出手段、大胸筋領域抽出手段、大胸筋直線設定手段、基準点設定手段、平行移動手段、及び、回転手段を備えたプログラムが記憶されている。そして、CPUはこれらのプログラムを読み出して、画像データ取込部から取り込まれた画像データに対し、左右の乳房について位置合わせ処理を行う。
従って、本発明によれば、X線撮影時に乳房を左右両側から挟み込んで押圧する位置が、左右の乳房について同一ではなく、左右のX線画像における乳房の位置に回転方向のずれが生じても、大胸筋直線の角度に基づいて、その回転方向のずれを補正することができる。すなわち、X線撮影時の押圧操作には影響を受けない大胸筋を位置合せに利用することにより、より正確に左右の乳房のX線画像の位置合わせを行うことができる。
また、スキンラインはX線撮影時の押圧操作によって影響を受けやすく、乳頭がスキンラインの先端の位置になかったり、乳頭がスキンラインから突出しなかったりすることがあり得るが、大胸筋直線との関係によってスキンライン上に基準点を定めることにより、左右の乳房の基準点をより正確に一致させることができる。加えて、大胸筋とスキンラインとの関係により定めた基準点と、大胸筋直線との双方を位置合わせの基準とすることにより、スキンラインに特徴が表れにくいためにその形状によっては位置合わせが困難である小さい乳房についても、より正確に位置合わせを行うことができる。
また、本発明にかかる乳房X線画像処理装置は、上記構成に加え、「前記位置合わせ手段によって位置合わせ処理された左右の乳房のX線画像のそれぞれについて、対応する座標により特定される左対象領域及び右対象領域の画像の特徴を表現する特徴ベクトルを、高次局所自己相関関数により求める特徴抽出手段と、求められた前記左対象領域及び前記右対象領域のそれぞれの前記特徴ベクトル間のユークリッド距離を算出し、前記左対象領域及び前記右対象領域の画像の類似度を判断する類似度判断手段とを」具備するものとすることができる。
「類似度」は、対比されるX線画像が、どの程度似ているかを示す尺度であり、後述のユークリッド距離Dによって定義される。上記のように、正常な組織は左右の乳房でほぼ等しいため、検出された陰影が異常陰影であるか正常組織であるかを、左右の乳房のX線画像の類似度によって判断することができる。
従って、次数や変位方向(a1,a2,・・・・aN)の取り方により高次自己相関関数は無数に存在する。ここで、画像の認識においては、離れた点間の相関よりも局所的な領域内での点間の相関が重要であると考えられる。そこで、変位方向を参照点rを中心とした局所領域に限定したものが、「高次局所自己相関関数」である。そして、局所領域内の全画素に対し、それぞれの局所パターンに対応する画素値の積をとって総和し、これを特徴量とする。すなわち、次数と変位方向の取り方により存在する局所パターンがn種類であるとすると、n個の要素をもつ特徴ベクトルによって、画像の特徴が表現される。
ここで、類似度を比較しようとする二つの対象領域A,Bについて、高次局所自己相関関数により求めた特徴ベクトルを、それぞれFAi,FBi(i=1,2,・・・,n)とすると、ユークリッド距離Dは次の数2で表される。そして、対象領域A及びBが同一であれば、FAi=FBiであるからユークリッド距離D=0であり、Dが短い(小さい)ほど、二つの対象領域A,Bの「類似度」は高いと判断される。
従って、本発明によれば、予め類否の境界となるD値を設定しておくことにより、位置合わせ処理された左右の乳房のX線画像において、特定の対象領域の画像が類似しているか否かを判断することができる。これにより、正常な左右の乳房における組織はほぼ等しいことを利用し、例えば、別途の検出手段により、一方の乳房、例えば左乳房のX線画像において異常陰影とみられる陰影が検出されたとき、その陰影が検出された左対象領域と、右の乳房のX線画像において左対象領域と対応する座標によって特定される右対象領域について、それぞれ特徴ベクトルを求め、ユークリッド距離を算出することにより、類似度が低ければ異常陰影であり、類似度が高ければ偽陽性候補であると判断し、偽陽性候補を削除することができる。
そして、局所パターンは画像の局所的な形状による特徴をよく反映し、これらが統計的に処理されて特徴ベクトルが求められるため、乳腺組織の微細な構造までは左右の乳房で完全に一致していなくても、左右の画像の類否を判断することが可能となる。これにより、従来、左右のX線画像の画素値の差のみによっては判断することが困難であった、乳腺組織の索状の陰影による偽陽性候補を、効果的に削除することができる。また、特徴ベクトルひいてはユークリッド距離は、対比されるX線画像のそれぞれの画素値から、直接的に算出することができるため、極めて簡易に類似度を判断することができる。
以上のように、本発明の第一の効果として、左右の乳房のX線画像について、より正確に位置合わせ処理を行う乳房X線画像処理装置を提供することができる。更に、第二の効果として、偽陽性候補を充分に削除することができる乳房X線画像処理装置を提供することができる。
以下、本発明の最良の一実施形態である乳房X線画像処理装置について、図1乃至図5に基づいて説明する。ここで、図1は本実施形態の乳房X線画像処理装置の構成を示すブロック図であり、図2は図1の乳房X線画像処理装置の画像処理部の機能的構成を示すブロック図であり、図3は図1の乳房X線画像処理装置による画像処理の流れを示すフローチャートであり、図4及び図5は図1の乳房X線画像処理装置による画像処理を説明する説明図である。
本実施形態の乳房X線画像処理装置1は、図1に示すように、内外斜位方向でX線撮影された乳房のX線画像の画像データを取り込む画像データ取込部2と、取り込まれた画像データを処理する画像処理部10と、画像処理部10に対する種々の指令の入力を行うキーボード、マウス、ポンティングデバイス等の操作入力部3と、画像データに基づいてX線画像を液晶ディスプレイ等による画面に表示する画像表示部5、及び画像処理部10によって処理された画像データに基づいてX線画像をフィルムや紙媒体に印刷する画像印刷部6によって構成される画像出力部4とを具備している。
更に詳細に説明すると、画像データ取込部2はCD−R等の記憶媒体に格納された画像データを読取るドライブ装置である。本実施形態の画像データは、フィルム/増感紙システムでX線撮影されたアナログ画像を、レーザスキャナでサンプリング間隔0.05mm、濃度分解能12bitでデジタル化したものであり、画像データ取込部2によって取り込んだ後0.4mm相当に縮小してから用いている。なお、サンプリング間隔は、微小石灰化像の検出のためには0.05〜0.1mm、腫瘤陰影の検出のためには0.1〜0.4mmとするのが一般的である。また、濃度分解能は10〜12bit(1024〜4096階調)とするのが一般的であるが、アナログ画像の淡い陰影を忠実にデジタル画像として再現するためには、充分な階調数によりデジタル化されることが望ましい。
画像処理部10は、ハード構成として、主記憶装置12と、主記憶装置12に記憶されたプログラムに従って演算や処理を行う中央処理装置(CPU)11と、ハードディスク等の補助記憶装置15とを具備している。ここで、補助記憶装置15には、画像データ取込部2から取り込まれた画像データや、画像処理部10によって処理された画像データを記憶することができる。
また、画像処理部10は、機能的構成として、図2に示すように、左右の乳房のX線画像の位置合わせ処理を行う位置合わせ手段20と、異常陰影を検出する異常陰影検出手段30とを具備している。ここで、位置合わせ手段20は、画像データ取込部2が取り込んだ画像データからスキンラインを検出して乳房領域を抽出する乳房領域抽出手段21と、大胸筋ラインを検出して大胸筋領域を抽出する大胸筋領域抽出手段22と、大胸筋ラインを直線に近似し大胸筋直線を設定する大胸筋直線設定手段23と、大胸筋直線とスキンライン上の全ての画素点との最短距離を検出し、検出された最短距離が最長のスキンライン上の画素点を基準点として設定する基準点設定手段25と、左右の乳房のX線画像についてぞれぞれ設定された基準点の座標が、中心線に対して線対称となるように、左右の乳房の少なくとも何れか一方の画像データの座標を平行移動させる平行移動手段26と、左右の乳房のX線画像についてぞれぞれ設定された大胸筋直線が中心線に対して線対称となるように、左右の乳房の少なくとも何れか一方の画像データの座標を基準点を中心に回転させる回転手段27とを備えている。
更に、異常陰影検出手段30は、位置合わせ手段20によって位置合わせされた左右の乳房のX線画像について、病変部をできるだけ漏れなく抽出する一次検出手段31と、一次検出手段31により一方の乳房の画像で異常とみられる陰影が検出されたとき、検出された陰影を含む対象領域を決定する対象領域決定手段32と、一方のX線画像で決定された対象領域、及び、他方のX線画像で対応する座標により特定される対象領域のそれぞれについて、高次局所自己相関関数により特徴ベクトルを求める特徴抽出手段33と、左右の画像のそれぞれの対象領域(左対象領域及び右対象領域)について求められた特徴ベクトル間のユークリッド距離を算出し、左対象領域と右対象領域との類似度を判断する類似度判断手段34と、判断された類似度に基づいて、異常陰影と判断された場合はその判断を一次検出結果に付加し、偽陽性候補であると判断された場合はこれを削除する偽陽性候補削除手段35とを備えている。
次に、本実施形態の乳房X線画像処理装置1における処理の流れについて、主に図3乃至図5に基づいて説明する。ここで、図3のフローチャートにおいてステップS3からステップS9が位置合わせ手段20に相当し、ステップS10からステップS19が本発明の異常陰影検出手段30に相当する。
まず、左右の乳房についてMLO方向からX線撮影された画像データが画像取込部2から取り込まれ(ステップS1)、補助記憶装置15に記憶される(ステップS2)。この画像データによるX線画像は、右の乳房について図4(a)に概略図を示したように、X線が直接透過し黒く表れた乳房外領域R3と、乳腺組織R1a及び脂肪組織R1bを含む乳房領域R1と、乳腺組織R1aや脂肪組織R1bよりX線の吸収率が高く低濃度に描出される大胸筋領域R2とに大別される。
そこで、まず乳房外領域R3内の点を抽出してから乳房のある方向に走査し、画素値が大きく変化する画素を検出し、そのような点の集合としてスキンラインN1を検出し、乳房領域R1を抽出する(ステップS3,図4(b)参照)。更に、同様に画素値の変化を利用して、大胸筋ラインN2を検出し、大胸筋領域R2を抽出する(ステップS4,図4(c)参照)。なお、本実施形態では、画素値の差異をより明確にするために、適宜の閾値によって二値化処理を行うことによって、大胸筋ラインN2を検出している。
大胸筋ラインN2はわずかに湾曲しているため、直線に近似して大胸筋直線N3を設定する(ステップS5、図4(d)参照)。更に、スキンラインN1上の全ての画素点について、大胸筋直線N3との最短距離を算出する(ステップS6、図4(e)参照)。この最短距離は、大胸筋直線N3上の点と、その点を通り大胸筋直線N3に直交する線がスキンラインN2と交差する点との距離から求めることができる。或いは、スキンラインN2上の点から大胸筋直線N3におろした垂線の長さから求めることができる。そして、求められた最短距離の長さが最大LmaxであるスキンラインN2上の画素点を、位置合せの基準点Pとして設定する(ステップS7、図4(f)参照)。通常、この基準点Pは乳頭の位置と一致するが、スキンラインN2の先端に乳頭が位置しない場合、スキンラインN2から乳頭が突出していない場合、手術により乳頭を切除した患者の場合であっても、位置合せの基準点Pが定められる。
なお、ステップS3からステップS7の処理は、左右の乳房の画像データのそれぞれについて行う。また、これらのステップは、画像データに基づいたX線画像を、画像表示部5によって画面に表示しながら行うこともできる。
以降は、ぞれぞれ基準点Pが設定された左右の乳房の画像データを一対として、位置合わせ処理を行う。以下では、左右のX線画像を画面表示しつつ位置合わせ処理を行う場合を例示するが、画面に表示することなく画像データのみを処理するものであっても良い。まず、左右の乳房のX線画像を単純に並列表示すると、図5(a)に示すように、左画像LI及び右画像RIでは、それぞれの基準点Pの座標が、中心線Mに対して上下及び左右方向で線対称とはなっていないのが通常である。そこで、それぞれの基準点Pが中心線Mに対して線対称となるように、左画像LI及び右画像RIの少なくとも何れか一方の画像データの座標を上下、左右方向に平行移動させる(ステップS8、図5(b)参照)。
更に、左画像LI及び右画像RIのぞれぞれについて設定された大胸筋直線N3が、中心線Mに対して線対称となるように、左画像LI及び右画像RIの少なくとも何れか一方の画像データの座標を、基準点Pを中心に回転させる(ステップS9、図5(c)参照)。以上の処理により、左右の乳房の画像データにおいて、乳房の位置合わせが完了する。なお、図5では、一方の画像(左画像LI)のみを平行移動させ、回転させる場合を図示したが、双方の画像を変位させても良い。
このようにして位置合わせが行われた左右のX線画像は、左右の乳房の対比観察による読影がし易くなっているため、例えば、左右に並列表示された状態で画像印刷部6によりフィルム印刷し、医師の診断に供することができる。本実施形態の乳房X線画像処理装置1は、医師による診断を更に高度に支援するため、後述のように、異常陰影の検出処理が可能な構成となっている。
異常陰影の検出処理には、画像に表れる病変部をできるだけ漏れなく抽出することを目的とする一次検出処理と、正常組織が異常陰影として誤って検出された偽陽性候補を削除する偽陽性候補削除処理が含まれる。そこで、まず、一次検出が行われる(ステップS10)。ここで、一次検出手段31としては、公知のアルゴリズムを利用することが可能であり、例えば、複数の閾値を用いて高濃度領域を逐次探索するもの、中心点に向かう勾配の集中度によって検出するもの、フィルタリングによって背景となる正常組織による濃度勾配の影響を低減し画像を強調して検出するもの等がある。
一次検出により、一方の乳房、例えば左の乳房の画像データにおいて、異常とみられる陰影(以下、「一次検出陰影」という)が検出された場合は(ステップS11においてYES)、その一次検出陰影が含まれる左対象領域を特定する(ステップS12)。本実施形態では、一次検出手段32として、一次検出陰影を含み面積が最小となる矩形の領域を、対象領域として決定するアルゴリズムを用いている。そして、他方の右の乳房の画像データにおいて、左対象領域と対応する座標として右対象領域を特定し、左対象領域及び右対象領域のそれぞれについて、高次局所自己相関関数により特徴ベクトルを求める(ステップS13)。
ここで、本実施形態では、高次局所自己相関関数として、高次自己相関関数の次数を2(N=0,1,2)に限定し、変位方向を参照点r周りの局所的な3×3画素領域に限定したものを考える。この場合、平行移動により等価となる特徴を除き、二値画像の場合は25種類、濃淡画像の場合は35種類の局所パターンが得られる。従って、各局所パターンについて、各画素の濃淡に対応する画素値との組合わせの数だけ異なる局所パターンが存在することとなるが、本実施形態では、画像データの階調数を4096階調(12bit)から64階調(6bit)に変換する前処理を行うことにより、画素値の範囲を設定している。
そして、左対象領域及び右対象領域のそれぞれの領域内の全画素について、それぞれの局所パターンに対応する画素値の積をとり総和して得られた左特徴ベクトルFAi、、及び右特徴ベクトルFBi(i=1,2,・・・n)から、上記の数2を用いてユークリッド距離Dを算出する(ステップS14)。ここで、予め、類似と非類似との判断の境界となるユークリッド距離Dbを定めておけば、算出されたDをDbと対比し、類否の判断をすることができる。例えば、D<Dbであれば(ステップS15においてYES)、左対象領域と右対象領域の画像は類似しており、一方の乳房で発見された一次検出陰影は、正常な組織が誤って検出された偽陽性候補と判断される。この偽陽性候補は、偽陽性候補削除手段35によって、一次検出結果からそのデータが削除される(ステップS16)。
一方、D≧Dbである場合は(ステップS15においてNO)、左対象領域と右対象領域の画像は類似していないことから、一方の乳房で発見された一次検出陰影は、他方の乳房にはみられない陰影であり、異常陰影と判断される。そして、この判断による情報は、偽陽性候補削除手段35によって、一次検出の結果に付加される(ステップS17)。なお、ここで言う異常陰影であるとの判断は、医師に情報を提示するための、乳房X線画像処理装置1としての判断であり、最終的な診断は医師によって行われる。
その後、類似度の判断結果による情報が付加されていない一次検出陰影の有無が確認され(ステップS18)、類似度について未判断である一次検出陰影がある場合は(ステップS18においてNO)、ステップS12に戻って上記と同様の類似度の判断が行われる。一方、全ての一次検出陰影について判断がなされた場合は(ステップS18においてYES)、位置合わせ処理、一次検出処理及び類似度の判断結果に関する情報も含めて、左右の乳房の画像データが補助記憶装置15に記憶される(ステップS19)。また、必要に応じて、上記の処理が行われた画像データを、画像印刷部6によりフィルムや紙媒体に印刷しても良い。そして、終了指令の有無が確認され、終了指令がなければ(ステップS20においてNO)、ステップS1に戻り新たな画像データの取り込みが行われ、終了指令がある場合は(ステップS20においてYES)、一連の画像処理を終了する。なお、一次検出によって異常とみられる陰影が検出されなかった場合は(ステップS11においてNO)、その情報が画像データに付加されて記憶される(ステップS19)。
上記のように、本実施形態の乳房X線画像処理装置1によれば、X線撮影時に乳房を左右両側から挟み込んで押圧する位置が左右の乳房で異なり、左右のX線画像における乳房の位置に回転方向のずれが生じても、押圧操作には影響を受けない大胸筋直線N3の角度に基づいて、その回転方向のずれを補正することにより、より正確に左右のX線画像の位置合わせ処理を行うことができる。
また、大胸筋直線N3との関係によってスキンラインN1上に基準点Pを定めることにより、乳頭がスキンラインの先端にない場合や、乳頭がスキンラインから突出しない場合であっても、左右の乳房のX線画像における基準点Pをより正確に一致させることができる。加えて、大胸筋直線N3とスキンラインN2との関係により定めた基準点Pと、大胸筋直線N3との双方を位置合わせの基準とすることにより、スキンラインの形状によっては位置合わせが困難である小さい乳房についても、より正確に位置合わせを行うことができる。
更に、高次局所自己相関関数による局所パターンは、画像の局所的な形状による特徴をよく反映し、これらが統計的に処理されて特徴ベクトルが求められるため、乳腺組織の微細な構造までは左右の乳房で完全に一致していなくても、左右の画像の類否を判断することが可能となる。これにより、一次検出手段によって異常陰影の可能性があるとして検出された一次検出陰影から、乳腺組織の索状の陰影を含む偽陽性候補を効果的に削除することができる。そのため、X線画像の画素値から直接的に算出できる極めて簡易な処理でありながら、偽陽性候補を削減し、異常陰影である可能性が高い陰影に関する情報を医師に提示して、医師による乳がんの診断を高度に支援することができる。
以上、本発明について好適な実施形態を挙げて説明したが、本発明は上記の実施形態に限定されるものではなく、以下に示すように、本発明の要旨を逸脱しない範囲において、種々の改良及び設計の変更が可能である。
例えば、上記の実施形態では、アナログ画像をデジタル画像に変換して記憶媒体に格納された画像データを、ドライブ装置から取り込む場合を例示したが、CRやFPDなどのデジタル系撮像装置から、直接画像データを取り込んでも良い。或いは、レーザスキャナ等のデジタル化装置から、ハードウェアインターフェースを介して伝送された画像データを取り込むものであっても良い。また、本発明の乳房X線画像処理装置自体が、デジタル化装置を具備するものとすることができる。
更に、上記の実施形態では、一次検出された陰影を含む対象領域を自動的に決定する手段を備える場合を例示したが、これに限定されない。例えば、一次検出の結果を画面表示し、作業者がマウスやポインティングデバイス等を使用して画面上で対象領域を指定し、或いは、座標値を入力することによっても、対象領域を特定することができる。
1 乳房X線画像処理装置
2 画像データ取込部
20 位置合わせ手段
21 乳房領域抽出手段
22 大胸筋領域抽出手段
23 大胸筋直線設定手段
25 基準点設定手段
26 平行移動手段
27 回転手段
33 特徴抽出手段
34 類似度判断手段
2 画像データ取込部
20 位置合わせ手段
21 乳房領域抽出手段
22 大胸筋領域抽出手段
23 大胸筋直線設定手段
25 基準点設定手段
26 平行移動手段
27 回転手段
33 特徴抽出手段
34 類似度判断手段
Claims (2)
- 乳房の内外斜位方向のX線画像の画像データを取り込む画像データ取込部と、
取り込まれた前記画像データからスキンラインを検出し乳房領域を抽出する乳房領域抽出手段、
取り込まれた前記画像データから大胸筋ラインを検出し大胸筋領域を抽出する大胸筋領域抽出手段、
前記大胸筋ラインを直線に近似し大胸筋直線を設定する大胸筋直線設定手段、
前記大胸筋直線と前記スキンライン上の全ての画素点との最短距離を検出し、検出された前記最短距離が最長の前記スキンライン上の前記画素点を基準点として設定する基準点設定手段、
左右の乳房の前記X線画像についてぞれぞれ設定された前記基準点が、中心線に対して線対称となるように、左右の乳房の少なくとも何れか一方の画像データの座標を平行移動させる平行移動手段、及び、
左右の乳房の前記X線画像についてぞれぞれ設定された前記大胸筋直線が、前記中心線に対して線対称となるように、左右の乳房の少なくとも何れか一方の画像データの座標を前記基準点を中心に回転させる回転手段を備え、
左右の乳房の前記X線画像の位置合わせ処理を行う位置合わせ手段と
を具備することを特徴とする乳房X線画像処理装置。 - 前記位置合わせ手段によって位置合わせ処理された左右の乳房のX線画像のそれぞれについて、対応する座標により特定される左対象領域及び右対象領域の画像の特徴を表現する特徴ベクトルを、高次局所自己相関関数により求める特徴抽出手段と、
求められた前記左対象領域及び前記右対象領域のそれぞれの前記特徴ベクトル間のユークリッド距離を算出し、前記左対象領域及び前記右対象領域の画像の類似度を判断する類似度判断手段と
を更に具備することを特徴とする請求項1に記載の乳房X線画像処理装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006252326A JP2008073076A (ja) | 2006-09-19 | 2006-09-19 | 乳房x線画像処理装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006252326A JP2008073076A (ja) | 2006-09-19 | 2006-09-19 | 乳房x線画像処理装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2008073076A true JP2008073076A (ja) | 2008-04-03 |
Family
ID=39345774
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2006252326A Pending JP2008073076A (ja) | 2006-09-19 | 2006-09-19 | 乳房x線画像処理装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2008073076A (ja) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013248140A (ja) * | 2012-05-31 | 2013-12-12 | Ge Medical Systems Global Technology Co Llc | 医用装置、表示装置、およびプログラム |
JP2018034014A (ja) * | 2011-10-05 | 2018-03-08 | ニューヴェイジヴ,インコーポレイテッド | 画像表示の生成方法 |
JP2018038609A (ja) * | 2016-09-07 | 2018-03-15 | キヤノン株式会社 | 画像処理装置、その制御方法、およびプログラム |
US11941179B2 (en) | 2010-10-06 | 2024-03-26 | Nuvasive, Inc. | Imaging system and method for use in surgical and interventional medical procedures |
-
2006
- 2006-09-19 JP JP2006252326A patent/JP2008073076A/ja active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11941179B2 (en) | 2010-10-06 | 2024-03-26 | Nuvasive, Inc. | Imaging system and method for use in surgical and interventional medical procedures |
JP2018034014A (ja) * | 2011-10-05 | 2018-03-08 | ニューヴェイジヴ,インコーポレイテッド | 画像表示の生成方法 |
JP2013248140A (ja) * | 2012-05-31 | 2013-12-12 | Ge Medical Systems Global Technology Co Llc | 医用装置、表示装置、およびプログラム |
JP2018038609A (ja) * | 2016-09-07 | 2018-03-15 | キヤノン株式会社 | 画像処理装置、その制御方法、およびプログラム |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107480677B (zh) | 一种识别三维ct图像中感兴趣区域的方法及装置 | |
US7916912B2 (en) | Efficient border extraction of image feature | |
US8233692B2 (en) | Method of suppressing obscuring features in an image | |
JP4891637B2 (ja) | 画像解析装置 | |
JP6967602B2 (ja) | 検査支援装置、内視鏡装置、内視鏡装置の作動方法、及び検査支援プログラム | |
JP4640845B2 (ja) | 画像処理装置およびそのプログラム | |
JPWO2007000940A1 (ja) | 異常陰影候補検出方法、異常陰影候補検出装置 | |
JP2002162705A (ja) | 位相コントラスト放射線画像処理装置 | |
JP3758894B2 (ja) | マンモグラム画像診断支援装置 | |
JP2000316837A (ja) | 画像診断支援装置 | |
JP6492553B2 (ja) | 画像処理装置及びプログラム | |
JP3852125B2 (ja) | 2次元多階調画像におけるパタン検出方法及び装置 | |
JP4453321B2 (ja) | 医用画像処理装置及びプログラム | |
JP2008073076A (ja) | 乳房x線画像処理装置 | |
JP2006325640A (ja) | 異常陰影候補の表示方法及び医用画像処理システム | |
JP2005253685A (ja) | 画像診断支援装置及び画像診断支援プログラム | |
JP5541989B2 (ja) | 放射線画像処理装置および方法並びにプログラム | |
JP2006340835A (ja) | 異常陰影候補の表示方法及び医用画像処理システム | |
JP2006230910A (ja) | 画像処理装置及び画像処理方法 | |
JP4631260B2 (ja) | 画像診断支援装置、画像診断支援方法及びプログラム | |
JP2006334140A (ja) | 異常陰影候補の表示方法及び医用画像処理システム | |
JP2005080758A (ja) | 画像処理装置 | |
JP6642048B2 (ja) | 医療画像表示システム、医療画像表示プログラム及び医療画像表示方法 | |
Noviana et al. | Axial segmentation of lungs CT scan images using canny method and morphological operation | |
JP4483250B2 (ja) | 画像診断支援装置、画像診断支援方法及びプログラム |