(57)【特許請求の範囲】
【請求項１】
Ｘ線発生器から被検体の乳房に向けて曝射されたＸ線を検出器で検出することにより取得して画像としてＸＹ座標上に形成され、前記Ｘ線が前記乳房を透過せずに入射するＸ線の直接線領域の画素がＸ軸の一端側に配され、被検体の胸壁側の画素がＸ軸の他端側に配されたデジタルマンモグラフィ画像の医用画像処理装置であって、
前記デジタルマンモグラフィ画像を基に、Ｘ軸に沿って配された画素中の最大の画素値をＹ軸上の位置毎に求め、該求められた複数の前記最大の画素値を基に閾値を求める閾値算出部と、
前記デジタルマンモグラフィ画像に基づく画像に対し、Ｘ軸上の位置をＸ軸の一端側から他端側の方向へずらしながら検索し、前記閾値以上の画素値を有する画素の第１位置を検出する第１検索部と、
前記第１位置の周囲を検索することにより、前記閾値以上の画素値を有する画素の第２位置を検出し、それ以後は、周囲の検索により検出した画素の位置の周囲を検索することを繰り返すことにより、前記乳房の輪郭を求める第２検索部と、
を有する
ことを特徴とする医用画像処理装置。
【請求項２】
前記デジタルマンモグラフィ画像中で、前記直接線領域のみの画素値を予め定められた値以下の画素値に変換することにより、前記画像としての変換画像を生成する画像変換処理部を更に有し、
前記閾値算出部は、前記変換画像を基に閾値を求め、
前記第１検索部及び前記第２検索部は、前記変換画像に対して前記検出を行うものであることを特徴とする請求項１に記載の医用画像処理装置。
【請求項３】
前記デジタルマンモグラフィ画像をウェーブレット変換し、該ウエーブレット変換後の画像を高周波成分のみで逆変換することにより、前記画像としての変換画像を生成する画像変換処理部を更に有し、
前記閾値算出部は、前記変換画像を基に閾値を求め、
前記第１検索部及び前記第２検索部は、前記変換画像に対して前記検出を行うものであることを特徴とする請求項１に記載の医用画像処理装置。
【請求項４】
前記閾値算出部は、前記複数の前記最大の画素値の中央値を基に前記閾値を求めることを特徴とする請求項１から請求項３のいずれかに記載の医用画像処理装置。
【請求項５】
前記第１位置の周囲又は前記直前に検出した前記第２位置の周囲を前記ＸＹ座標上で検索するとき、前記第２検索部が検索する最初の方向は、前記第１位置又は前記直前に検出した前記第２位置から前記一端側の方向であることを特徴とする請求項１から請求項４のいずれかに記載の医用画像処理装置。
【請求項６】
前記第１位置の周囲又は前記直前に検出した前記第２位置の周囲を前記ＸＹ座標上で検索するとき、前記第２検索部が検索する方向は、前記被検体の右乳房に係る前記変換画像と、前記被検体の左乳房に係る前記変換画像とで互いに反対にしたことを特徴とする請求項５に記載の医用画像処理装置。
【請求項７】
前記検出された前記第１位置及び前記第２位置の各位置に隣り合う各々の位置の座標を基に、平滑化された位置を求める平滑化処理部を有することを特徴とする請求項１から請求項６のいずれかに記載の医用画像処理装置。
【請求項８】
前記第１位置の周囲は、前記第１位置にある画素に隣接する画素であり、
前記第２位置の周囲は、前記第２位置にある画素に隣接する画素であることを特徴とする請求項１から請求項７のいずれかに記載の医用画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、医用画像処理装置に関し、特に、Ｘ線発生器から被検体の乳房に向けて曝射されたＸ線を検出器で検出することにより取得されたデジタルマンモグラフィ画像の医用画像処理装置に関する。
【背景技術】
【０００２】
Ｘ線発生器から被検体の乳房に向けて曝射されたＸ線を検出器で検出することによりデジタルマンモグラフィ画像を取得し、取得した画像をＸＹ座標上に形成し、ＸＹ座標上に形成された画像のヒストグラム解析することにより、Ｘ線発生器から曝射され、乳房を透過した後に検出器に検出されたＸ線の領域（以下、「Ｘ線透過領域」という。）と、乳房を透過せずに検出器に検出されたＸ線の領域（以下、「直接線領域」という。）との境界に位置する画素の輝度値である境界値を求め（図１参照）、デジタルマンモグラフィ画像のＸＹ座標上で、Ｙ軸上の位置毎にＸ軸の一端側（胸壁側）から他端側（直接線領域側）の方向へ境界値以上の輝度値を最初に有する画素の位置を検索し（図２参照）、検索された位置に対して平滑化処理をするスキンライン検出法がある。ここで、スキンラインとは、乳房の皮膚の領域であって、Ｘ線透過領域が直接線領域と接する領域である。
【０００３】
デジタルマンモグラフィ画像のポジ表示をすると、直接線領域の画素の輝度値が高いため、直接線領域が真っ白に表示され、読影者の目が疲れる等の悪影響を及ぼす（図３（ａ）参照）。
【０００４】
この対策として、直接線領域を黒く表示するため、直線線領域のみの画素の輝度値を低くした黒化処理の技術がある（図３（ｂ）参照）。
【０００５】
また、Ｘ線撮影によって得られたデジタルマンモグラフィ画像において、ソーベルフィルタを用いて複数の画素間の輝度の変化に対応する輪郭データを取得し、輪郭データを２値化し、２値化された輪郭データを凸包処理を施すことにより最外輪郭を抽出し、抽出された最外輪郭の外部領域における画素の輝度を低減した画像データを生成する黒化処理の技術がある（例えば、特許文献１）。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００４−２８３２８１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、前記図１及び図２を参照して示す従来の技術では、Ｙ軸上の位置毎にＸ軸の一端側から他端側の方向へ境界値以上の輝度値を最初に有する画素の位置を検索するため、検索対象となるＸ軸上において、前記最初に有する画素の位置が必ずしもスキンラインの位置として最適でなく、スキンラインの位置として最適な位置が前記Ｘ軸上の他の位置に存在する場合があり、その結果、スキンラインの位置精度を低下させるおそれがある。
【０００８】
また、上記特許文献に記載された従来の技術では、ソーベルフィルタを用いることにより取得された輪郭データが、ボケやノイズの影響を受けやすいため、輪郭データが最外輪郭として最適でない場合があり、その結果、最外輪郭の位置精度を低下させるおそれがあるという問題点があった。
【０００９】
この発明は、上記の問題を解決するものであり、スキンラインの位置精度を向上させることが可能な医用画像処理装置を提供することを目的とする。
【課題を解決するための手段】
【００１０】
上記課題を解決するため、この発明は、直前に検出した直前位置の周囲中から次の位置を検出していくことに着目した。
具体的に、Ｘ線発生器から被検体の乳房に向けて曝射されたＸ線を検出器で検出することにより取得して画像としてＸＹ座標上に形成され、前記Ｘ線が前記乳房を透過せずに入射するＸ線の直接線領域の画素がＸ軸の一端側に配され、被検体の胸壁側の画素がＸ軸の他端側に配されたデジタルマンモグラフィ画像の医用画像処理装置であって、前記デジタルマンモグラフィ画像を基に、Ｘ軸に沿って配された画素中の最大の画素値をＹ軸上の位置毎に求め、該求められた複数の前記最大の画素値を基に閾値を求める閾値算出部と、前記デジタルマンモグラフィ画像に基づく画像に対し、Ｘ軸上の位置をＸ軸の一端側から他端側の方向へずらしながら検索し、前記閾値以上の画素値を有する画素の第１位置を検出する第１検索部と、前記第１位置の周囲を検索することにより、前記閾値以上の画素値を有する画素の第２位置を検出し、それ以後は、周囲の検索により検出した画素の位置の周囲を検索することを繰り返すことにより、前記乳房の輪郭を求める第２検索部と、を有するものである。
【発明の効果】
【００１１】
この発明によると、スキンラインの位置精度を向上させることが可能となる。
【００１２】
また、この発明の第１の形態によると、第１位置の周囲を検索するようにしたので、より適切な第２位置を検出することができ、例えば、Ｙ軸上の位置毎にＸ軸方向からスキンラインの位置を検索する場合に比較して、スキンラインの位置精度を向上させることが可能となる。
【図面の簡単な説明】
【００１３】
【図１】デジタルマンモグラフィ画像のヒストグラムを示す図である。
【図２】デジタルマンモグラフィ画像、及び、デジタルマンモグラフィ画像の一部を拡大した図であって、検出した直接線領域の輪郭の位置を斜めのハッチングで示した図である。
【図３】（ａ）黒化処理前を示す図、（ｂ）黒化処理後を示す図である。
【図４】本発明の実施形態に係るマンモグラフィ装置の機能ブロック図である。
【図５】（ａ）原画像を示す図、（ｂ）ウェーブレット変換後の画像を示す図、（ｃ）変換画像を示す図である。
【図６】全Ｘｐｒｏｆｉｌｅの最大値の中央値の算出方法を示す図である。
【図７】（ａ）右乳房の画像において、第１位置を検出する処理を示す図、（ｂ）右乳房の画像において、第２位置を検出する処理を示す図である。
【図８】（ａ）左乳房の画像において、第１位置を検出する処理を示す図、（ｂ）左乳房の画像において、第２位置を検出する処理を示す図である。
【図９】（ａ）〜（ｃ）右乳房の画像において、第２位置を順次検出する処理を示す図である。
【図１０】（ａ）〜（ｃ）左乳房の画像において、第２位置を順次検出する処理を示す図である。
【図１１】平滑化処理の一例を示す図である。
【図１２】スキンラインの位置検出を示すフローチャートである。
【図１３】第１位置の検索がＹ軸上の位置毎に行われることを示すフローチャートである。
【図１４】第１位置の検索を示すフローチャートである。
【図１５】第２位置の検索が第１位置の周囲又は直前に検出された第２位置の周囲で行われることを示すフローチャートである。
【図１６】第２位置の検索を示すフローチャートである。
【発明を実施するための形態】
【００１４】
（構成）
以下、本発明の実施形態に係る医用画像処理装置の構成について図４〜図１２を参照して説明する。図４は、本実施形態に係るマンモグラフィ装置１の機能ブロック図である。
【００１５】
マンモグラフィ装置１は、Ｘ線撮影台３とコンピュータ装置５とを備えている。
【００１６】
Ｘ線撮影台３は、基台部１０と、Ｃアーム１１を有する。Ｃアーム１１は、水平軸回りに回動可能なように基台部１０に支持されている。Ｃアーム１１を回転させることにより、頭尾方向、内外方向、内外斜方向等の撮影を行うことができる。
【００１７】
Ｃアーム１１には、Ｘ線発生装置１５、Ｘ線検出器１６、及び圧迫ユニット１７が取り付けられている。Ｘ線発生装置１５及びＸ線検出器１６は、Ｃアーム１１の両端部に配置されている。圧迫ユニット１７は、Ｘ線発生装置１５とＸ線検出器１６との中間に配置されている。
【００１８】
Ｘ線発生装置１５は、Ｘ線管１８と高電圧発生器１９とを有する。Ｘ線管１８は、高電圧発生器１９から管電圧の印加、及びフィラメント電流の供給を受けて圧迫ユニットに向けて所定のＸ線継続時間Ｘ線を発生する。印加する管電圧とＸ線継続時間とは、撮影制御部３７からの制御信号を受けて、撮影に適した値に調整されている。
【００１９】
Ｘ線管１８は、陰極フィラメントと陽極とを備えている。陽極は、Ｍｏ（モリブデン）を材質としたＭｏ陽極、Ｒｈ（ロジウム）を材質としたＲｈ陽極、ＭｏとＲｈとを混合してなるＭｏ・Ｒｈ陽極等である。これら陽極は、撮影制御部３７からの制御信号を受けて、随時切り替え可能である。フィラメント電流の供給を受けた陰極フィラメントは加熱せれ、熱電子を発生する。発生された熱電子は、陰極フィラメントと陽極との間に印加された管電圧によって、陽極に衝突されている。このように熱電子が陽極へ衝突することによりＸ線が発生されている。陽極に衝突する熱電子によって、管電流が流れる。管電流は、  フィラメント電流により調整されている。撮影時におけるＸ線線量の調節は、撮影制御部３７からの制御信号を受けて、管電流とＸ線継続時間との積である管電流時間積を調節することにより行なわれる。
【００２０】
Ｘ線管１８には、発生されたＸ線の線質を変更するための線質フィルタが取り付けられる。線質フィルタは、Ｍｏを材質としたＭｏフィルタや、Ｒｈを材質としたＲｈフィルタ、Ａｌ（アルミニウム）を材質としたＡｌフィルタ、或いはこれら材質を組み合わせてなるフィルタ等である。これら線質フィルタは、撮影制御部３７からの制御信号を受けて、随時切り替え可能である。
【００２１】
圧迫ユニット１７は、Ｘ線検出器１６の検出面に沿って接近及び離間可能なようにＣアーム１１によって支持されている。圧迫ユニット１７は、撮影制御部３７からの制御信号を受けて、動作することにより被検体の乳房を圧迫し、乳房厚を所定の状態にする。
【００２２】
Ｘ線検出器１６は、乳房を透過したＸ線を検出する検出素子を２次元に配列したフラット・パネル・ディテクタ等のデジタル検出器である。Ｘ線検出器１６は、その検出面中心とＸ線管１８の焦点とを結ぶＺ軸に沿って、Ｘ線管１８に接近及び離間可能なようにＣアーム１１によって支持されている。ここで、Ｙ軸をＸ軸及びＺ軸に直交する軸に規定する。つまり、ＸＹＺ座標系は、Ｘ軸を回転中心軸とした回転座標系である。Ｚ軸は乳房の厚さ方向を規定する軸であり、ＸＹ平面は乳房の厚さ方向に垂直な広がり方向を規定する軸である。
【００２３】
また、圧迫ユニット１７が上側圧迫板と下側圧迫板とを備え、Ｘ線検出器１６の検出面と、上側圧迫板と下側圧迫板とで圧迫された乳房との距離を離間させて拡大撮影を行う場合は、Ｚ軸方向のＸ線焦点−乳房間距離とＸ線焦点−Ｘ線検出器間距離とを調整することで、拡大率を撮影に適した状態にする。
【００２４】
Ｘ線検出器１６には、データ収集部（ＤＡＳ；Data Acquisition System）２１が接続されている。データ収集部２１は、Ｘ線検出器１６の各チャンネル（検出素子）の電流信号を収集する。データ収集部２１は、収集した電流信号をデジタル信号に変換し、前処理を行う。前処理されたデジタル信号は、投影データと呼ばれる。投影データは、データ収集部２１によってコンピュータ装置５に供給されている。
【００２５】
コンピュータ装置５は、画像発生部２３、操作部３５、撮影制御部３７、システム制御部３９、及び、医用画像処理装置４０を備えている。
【００２６】
画像発生部２３は、データ収集部２１からの投影データに基づいて画像のデータを発生する。通常、マンモグラフィ撮影によって得られた画像の生体領域には、乳房領域だけではなく大胸筋領域等の乳房領域外の領域を含む。
【００２７】
操作部３５は、撮影条件（管電圧、管電流時間積、陽極の材質、線質フィルタの材質、乳房厚、Ｘ線焦点−Ｘ線検出器間距離、拡大率等）を撮影制御部３７に設定するための操作パネルである。また、操作部３５は、撮影対象である乳房の左右いずれかを示すコードを撮影制御部３７に設定する。また、操作部３５は、Ｃアーム１１を動作させるためのインターフェースを備えており、その操作に応じてＣアーム１１はＺ軸回りに回動され任意の位置に設定されている。設定されたＣアーム１１の位置に応じて、撮影方向が決定されている。
【００２８】
撮影制御部３７は、操作部３５を介して設定された撮影条件（管電圧、管電流時間積、陽極の材質、線質フィルタの材質、乳房厚、Ｘ線焦点−Ｘ線検出器間距離、拡大率等）に基づいてＸ線撮影台３の各構成要素を制御することによって、Ｘ線撮影台３に設定に応じたＸ線撮影を行わせる。
【００２９】
次に、医用画像処理装置４０について、図４〜図１１を参照して説明する。図５は（ａ）原画像を示す図、（ｂ）ウェーブレット変換後の画像を示す図、（ｃ）変換画像を示す図、図６は、全Ｘｐｒｏｆｉｌｅ（各Ｙ軸についてのＸ軸方向の画素値）の最大値の中央値の算出方法を示す図、図７は（ａ）右乳房の画像において、第１位置を検出する処理を示す図、（ｂ）右乳房の画像において、第２位置を検出する処理を示す図、図８は（ａ）左乳房の画像において、第１位置を検出する処理を示す図、（ｂ）左乳房の画像において、第２位置を検出する処理を示す図、図９は、（ａ）〜（ｃ）右乳房の画像において、第２位置を順次検出する処理を示す図、図１０は、（ａ）〜（ｃ）左乳房の画像において、第２位置を順次検出する処理を示す図、図１１は平滑化処理の一例を示す図である。
【００３０】
医用画像処理装置４０は、画像記憶部２５、画像変換処理部２７、閾値算出部２９、第
１検索部３１、第２検索部３２、及び、平滑化処理部３３を備えている。
【００３１】
画像記憶部２５は、発生された画像のデータに、撮影条件と、撮影時の撮影方向を示すコードと、撮影した乳房の左右を示すコードとを関連付けて記憶する。
【００３２】
次に、画像変換処理部２７について図５を参照して説明する。画像変換処理部２７は、デジタルマンモグラフィ画像（原画像）をウェーブレット変換し、ウエーブレット変換後の画像を高周波成分のみで逆変換することにより、変換画像（高周波画像）を生成する。画像変換処理部２７は、変換画像（高周波画像）をＸＹ座標上に形成し、各画素をｎ×ｎのマトリクス状に配置する。ここで、Ｘ軸上のｎ番目の画素の位置をＸｎ、Ｙ軸上のｎ番目の画素の位置をＹｎとし、位置（Ｘｎ，Ｙｎ）に配置された画素の画素値をＧｎｎとする。デジタルマンモグラフィ画像（原画像）を図５（ａ）に示し、ウエーブレット変換後の画像（ＬＬ、ＬＨ、ＨＬ、ＨＨ）を図５（ｂ）に示し、ウエーブレット変換後の画像を高周波成分（ＬＨ、ＨＬ、ＨＨ）のみで逆変換した変換画像を図５（ｃ）に示した。ここで、Ｈは高周波を示し、Ｌは低周波を示す。
【００３３】
次に、閾値算出部２９について図６を参照して説明する。図６は、閾値の算出過程を説明するための模式図であって、ＸＹ座標上に形成された変換画像（高周波画像）の一部を示している。ＸＹ座標上の位置（Ｘ１，Ｙ１）、（Ｘ１，Ｙ２）、（Ｘ１，Ｙ３）、（Ｘ２，Ｙ１）、（Ｘ２，Ｙ２）、（Ｘ２，Ｙ３）、（Ｘ３，Ｙ１）、（Ｘ３，Ｙ２）、及び、（Ｘ３，Ｙ３）に配置された各画素の画素値２、４、３、５、３、６、８、２、及び、３を図６に示す。
【００３４】
閾値算出部２９は、変換画像（高周波画像）について、Ｘ軸に沿って配された画素中の最大の画素値をＹ軸上の位置毎に求め、求められた複数（ｎ個）の最大の画素値内の中央値を求め、中央値に予め定められたパラメータを乗算することにより閾値を求める。なお、ここで、中央値とは、画素値を大きさの順に１列に並べたときの真ん中の値をいう。
【００３５】
図６の例では、ＸＹ座標上の位置（Ｘ２，Ｙ１）、（Ｘ３，Ｙ２）、及び、（Ｘ１，Ｙ３）の画素値４、６、及び、８が最大値となる。また、ＸＹ座標上の位置（Ｘ３，Ｙ２）の画素値６が、上記３つの最大値（４、６、８）の内の中央値となる。閾値算出部２９は、中央値６に予め定められた係数３を乗算することにより、閾値１８を求める。係数３は、実験的に求められた数値である。３つの最大値（４、６、８）をＸｐｒｏｆｉｌｅの最大値として図６に示す。なお、閾値算出部２９は、変換画像（高周波画像）に代えてデジタルマンモグラフィ画像（原画像）を基に、全Ｘｐｒｏｆｉｌｅの最大値の内の中央値を求め、求めた中央値に同じく実験的に求められた係数を乗算することにより、閾値を求めるようにしても良い。
【００３６】
次に、第１検索部３１について図７（ａ）及び図８（ａ）を参照して説明する。図７（ａ）は、検索過程を説明するための模式図であり、右端側を胸壁側、左端側を直線線領域側とした右乳房の変換画像（高周波画像）の一部を示す。また、図８（ａ）は、同様に、左端側を胸壁側、右端側を直接線領域側とした左乳房の変換画像の一部を示す。なお、説明を簡明にするために、図７（ａ）を反転することにより図８（ａ）を示した。また、図７（ａ）及び図８（ａ）において、位置Ｘ１側を胸壁側とし、位置Ｘｎ側を直接線領域側とした。
【００３７】
以下、右乳房の画像及び左乳房の画像に対し、第１位置Ｓ１の検出をそれぞれ行う第１検索部３１の動作の説明において、直接線領域側を一端側とし、胸壁側を他端側とし、位置Ｙ１側を一端側とし、位置Ｙｎ側を他端側とした。なお、閾値算出部２９により求められた閾値を３０とした。また、Ｘ線透過領域と直接線領域との境界に位置する画素の輝度値である輝度値を実験的に求め３とした。
【００３８】
ＸＹ座標上の位置（Ｘ１，Ｙ１）、（Ｘ１，Ｙ２）、（Ｘ１，Ｙ３）、（Ｘ１，Ｙ４）、（Ｘ１，Ｙ５）、…、（Ｘ４，Ｙ１）、（Ｘ４，Ｙ２）、（Ｘ４，Ｙ３）、（Ｘ４，Ｙ４）、及び、（Ｘ４，Ｙ５）に配置された画素の画素値１、２、１、１、１、…、２０、２５、３２、３３、及び、２５を図７（ａ）に示す。また、図７（ａ）を反転することにより示した図８（ａ）において、ＸＹ座標上の位置に、同じ画素値を有する画素が配されている。
【００３９】
第１検索部３１は、ＸＹ座標上に形成された変換画像を、Ｙ軸上の位置をＹ軸の一端側から他端側の方向へ１つずつずらし、ずらしたＹ軸上の位置に対してＸ軸上を検索する場合に、Ｘ軸の一端側（直接線領域側）から他端側（胸壁側）の方向へ検索し、閾値３０以上の画素値を最初に有する画素の第１位置Ｓ１を検出する。
【００４０】
ここで、ずらしたＹ軸上の位置に対してＸ軸上を検索する場合とは、Ｘ軸上を検索するための所定の条件を満たした場合をいう。所定の条件とは、例えば、Ｙ軸上の位置に配された画素の画素値が境界値３以上である場合をいう。境界値３は、変換画像においてＸ線透過領域と直接線領域との境界に位置する画素の輝度値である。なお、所定の条件を、Ｙ軸上の位置に配された画素の画素値が０以上である場合としても良い。この場合、Ｙ軸上の位置に配された画素の画素値は０から２５５までのいずれかの数値を有するから、Ｘ軸上を検索するための所定の条件を必ず満たすこととなる。
【００４１】
さらに、Ｘ軸の一端側（直接線領域側）から他端側（胸壁側）の方向へ検索するとは、Ｘ軸の一端側から他端側の方向へ検索するときに、最終的に、第１位置を検索することをいい、例えば、始めに、Ｘ軸の他端側から一端側の方向に検索し、検索した位置に配置された画素の画素値が境界値３未満の場合に、Ｘ軸の一端側から他端側の方向へ検索することをいう。さらに、Ｘ軸上を検索するとは、Ｘ軸上の全位置又は一部位置の検索を含むことをいう。さらに、画素値を最初に有する画素とは、検索対象となったＸ軸上の位置に配された画素の内から、所定条件（閾値以上の画素値）を最初に有する画素をいう。
【００４２】
次に、具体的に、右乳房の変換画像において、第１位置Ｓ１を検出する第１検索部３１の動作を図７（ａ）を参照して説明する。第１検索部３１の動作の手順を（i）〜（ｘiii）のまでの番号を付して図７（ａ）に示した。ここで、上記のＸ軸上を検索するための所定の条件とは、Ｙ軸上の位置に配置された画素の画素値が前記境界値３以上である場合とした。また、閾値を３０とした。
【００４３】
（i）検索対象の位置の画素値が境界値３未満なら下（↓）へ移動する。（ii）検索対象の位置の画素値が境界値３以上なら左（←）へ移動する。（iii）移動した検索対象の位置の画素値が境界値３以上なら左（←）へ移動する。（iv）移動した検索対象の位置の画素値が境界値３未満なら右（→）へ移動する。
【００４４】
（v）移動した検索対象の位置の画素値が閾値３０未満なら右（→）へ移動する。（vi
）右の移動により検索対象の位置が右端に戻ったら下（↓）へ移動する。（vii）移動した検索対象の位置の画素値が境界値３以上なら左（←）へ移動する。（viii）移動した検索対象の位置の画素値が境界値３以上なら左（←）へ移動する。（ix）さらに、移動した検索対象の位置の画素値が境界値３以上なら左（←）へ移動する。（x）さらに、移動した検索対象の位置の画素値が境界値３以上なら左（←）へ移動する。（xi）移動した検索対象の位置の画素値が閾値３０未満なら右（→）へ移動する。（xii）さらに、移動した検索対象の位置の画素値が閾値３０未満なら右（→）へ移動する。（xiii）移動した検索対象の位置の画素値が閾値３０以上なら第１位置Ｓ１を検出する。
【００４５】
なお、第１位置Ｓ１の検索においては、Ｙ軸上の位置毎に、検索対象の位置を左端側（一端側）から右端側（他端側）に移動させるようにしても良い。本実施形態においては、第１位置Ｓ１の検索において、第１検索部３１は、Ｙ軸上の位置毎に、検索対象の位置を先ずＸ軸の右端側である他端側（胸壁側）から左端側である一端側（直接線領域側）に移動させ、検索対象の位置の画素値が３未満であると判断した後、検索対象の位置を左端側（一端側）から右端側（他端側）に移動させたので、Ｙ軸上の位置毎に、検索対象の位置を単に左端側（一端側）から右端側（他端側）に移動させた場合に比較して、迅速に第１位置Ｓ１を検出することが可能となる。
【００４６】
次に、図８（ａ）で示す左乳房の変換画像において、第１位置Ｓ１を検出する第１検索部３１の動作を説明する。なお、左乳房の変換画像において第１検索部３１の動作の手順は、図７（ａ）を用いて説明した第１位置Ｓ１を検出する場合の第１検索部３１の動作の手順と同じであり、その説明を省略する。左乳房の変換画像において、第１位置Ｓ１を検出する第１検索部３１の動作の手順を（i）〜（xiii）のまでの番号を付して図８（ａ）に示した。
【００４７】
次に、第２検索部３２について図７（ｂ）、図８（ｂ）、図９、及び、図１０を参照して説明する。図７（ｂ）は、右端側を胸壁側、左端側を直線線領域側とした右乳房の変換画像（高周波画像）の一部を示し、検出された第１位置Ｓ１（図７（ａ）参照）の周囲を検出する第２検索部３２の動作を示す。また、図８（ｂ）は、左端側を胸壁側、右端側を直接線領域側とした左乳房の変換画像の一部を示し、検出された第１位置Ｓ１（図８（ａ）参照）の周囲を検出する第２検索部３２の動作を示す。なお、図７（ａ）及び図８（ａ）と同様に、説明を簡明にするために、図７（ｂ）を反転することにより図８（ｂ）を示した。また、図７（ｂ）及び図８（ｂ）において、位置Ｘ１側を胸壁側とし、位置Ｘｎ側を直接線領域側とした。図７（ｂ）及び図８（ｂ）の場合においても、Ｘ軸上を検索するための所定の条件とは、Ｙ軸上の位置に配置された画素の画素値が前記境界値３以上である場合とした。また、閾値を３０とした。
【００４８】
第２検索部３２は、第１位置Ｓ１の周囲をＸＹ座標上で検索することにより、閾値３０以上の画素値を最初に有する画素の第２位置Ｓ２を検出し、それ以後は、直前に検出した第２位置Ｓ２の周囲をＸＹ座標上で検索することにより、閾値３０以上の画素値を最初に有する画素の第２位置Ｓ２をＹ軸の他端側（下端）まで検出していく。
【００４９】
ここで、第１位置Ｓ１の周囲とは、第１位置Ｓ１に対して１画素（１コマ）分だけ離れた周囲をいう。また、第２位置Ｓ２の周囲とは、第２位置Ｓ２に対して１画素（１コマ）分だけ離れた周囲をいう。第２位置Ｓ２の検索において、第１位置Ｓ１の周囲及び直前に検出された第２位置Ｓ２の周囲を検索する理由は、スキンラインの各位置が連続して存在し、第１位置Ｓ１の周囲及び直前に検出された第２位置Ｓ２の周囲に、連続する次の第２位置Ｓ２が存在することを実験的に分かっていることに拠る。
【００５０】
さらに、第１位置Ｓ１の周囲をＸＹ座標上で検索するとは、検索する方向を被検体の右乳房に係る前記変換画像と、被検体の左乳房に係る前記変換画像とで互いに反対に検索することをいい、例えば、図７（ｂ）に示す胸壁側を右端側とした右乳房の変換画像において、反時計方向に検索し、図８（ｂ）に示す胸壁側を左端側とした左乳房の変換画像において、時計方向に検索することをいう。さらに、第１位置Ｓ１の検出する場合と同様に、閾値３０以上の画素値となる。なお、所定値を第１位置Ｓ１の検出する場合と異ならせても良い。さらに、画素値を最初に有する画素とは、検索対象となった周囲の位置に配された画素の内から、所定条件（閾値以上の画素値）を最初に有する画素をいう。さらに、直前に検出した前記第２位置とは、既に検出された１又は複数の第２位置の内で、最新に検出された第２位置をいう。
【００５１】
次に、図７（ｂ）で示す右乳房の変換画像において、第１位置Ｓ１の周囲を検索することにより、第２位置Ｓ２を検出する第２検索部３２の動作を説明する。
【００５２】
（i）検索対象の位置を第１位置の左（←）へ移動する。（ii）検索対象の位置の画素値が閾値３０未満なら下（↓）へ移動する。（iii）検索対象の位置の画素値が閾値３０以上なら第２位置を検出する。
【００５３】
右端側を胸壁側とした右乳房の変換画像において、第２検索部３２の検索する最初の方向が第１位置Ｓ１の左方向（胸壁側から直接線領域側の方向）である理由は、第１位置Ｓ１に連続し、第１位置Ｓ１の次にスキンラインの位置となるべき第２位置Ｓ２が第１位置Ｓ１の左側に存在する可能性が大きいことが実験的に分かっていることに拠る。
【００５４】
次に、図８（ｂ）で示す左乳房の変換画像において、第１位置Ｓ１の周囲を検索することにより、第２位置Ｓ２を検出する第２検索部３２の動作を説明する。なお、左乳房の変換画像において第２検索部３２の動作の手順は、図８（ａ）を用いて説明した第２位置Ｓ２を検出する場合の第２検索部３２の動作の手順と同じであり、その説明を省略する。左乳房の変換画像において、第２位置Ｓ２を検出する第２検索部３２の動作の手順を（i）〜（iii）のまでの番号を付して図８（ｂ）に示した。
【００５５】
なお、左端側を胸壁側とした左乳房の変換画像において、第２検索部３２の検索する最初の方向が第１位置Ｓ１の右方向（胸壁側から直接線領域側の方向）である理由は、第１位置Ｓ１に連続し、第１位置Ｓ１の次にスキンラインの位置となるべき第２位置Ｓ２が第１位置Ｓ１の右側に存在する可能性が大きいことが実験的に分かっていることに拠る。
【００５６】
以上は、第１位置Ｓ１の周囲を検索することにより、第２位置Ｓ２を検出する第２検索部３２の動作を説明した。
【００５７】
次に、直前に検出された第２位置Ｓ２の周囲を検索することにより、新たな第２位置Ｓ２を検出する第２検索部３２の動作について図９及び図１０を参照して説明する。
【００５８】
図９は、右端側を胸壁側、左端側を直線線領域側とした右乳房の変換画像（高周波画像）の一部を示す。また、図１０は、左端側を胸壁側、右端側を直接線領域側とした左乳房の変換画像の一部を示す。なお、説明を簡明にするために、図９を反転することにより図１０を示した。また、図９及び図１０において、位置Ｘ１側を胸壁側とし、位置Ｘｎ側を直接線領域側とした。以下、右乳房の画像及び左乳房の画像に対し、第１位置Ｓ１の検出をそれぞれ行う第１検索部３１の動作の説明において、直接線領域側を一端側とし、胸壁側を他端側とし、位置Ｙ１側を一端側とし、位置Ｙｎ側を他端側とした。
【００５９】
先ず、図９で示す右乳房の変換画像において、直前に検出された第２位置Ｓ２の周囲を検索することにより、新たな第２位置Ｓ２を検出する第２検索部３２の動作について説明する。
【００６０】
第１位置をＸＹ座標上の位置（Ｘ２，Ｙ２）として図９（ａ）に示す。先ず、検索対象の位置を第１位置Ｓ１の左（←）へ移動する。移動した検索対象の位置（Ｘ３，Ｙ２）の画素値が閾値３０未満なら下（↓）へ移動する。移動した検索対象の位置（Ｘ３，Ｙ３）の画素値が閾値３０未満なら右（→）へ移動する。移動した検索対象の位置（Ｘ２，Ｙ３）の画素値が閾値３０以上なら第２位置を検出する。位置（Ｘ２，Ｙ３）が直前に検出された第２位置Ｓ２となる。
【００６１】
次に、直前に検出された第２位置Ｓ２の周囲をＸＹ座標上で検索する第２検索部３２の動作を図９（ｂ）に示す。
【００６２】
先ず、検索対象の位置を第２位置Ｓ２の左（←）へ移動する。移動した検索対象の位置（Ｘ３，Ｙ３）の画素値が閾値３０未満なら下（↓）へ移動する。移動した検索対象の位置（Ｘ３，Ｙ４）の画素値が閾値３０以上なら第２位置を検出する。位置（Ｘ３，Ｙ４）が直前に検出された第２位置Ｓ２となる。また、胸壁側を右端側とした右乳房の変換画像において、第２検索部３２による検索方向は、反時計回りとなる。直前に検出された第２位置Ｓ２を図９（ｃ）に位置（Ｘ３，Ｙ４）で示し、反時計回りの矢印で示した。
【００６３】
次に、図１０で示す左乳房の変換画像において、直前に検出された第２位置Ｓ２の周囲を検索することにより、新たな第２位置Ｓ２を検出する第２検索部３２の動作について説明する。なお、左乳房の変換画像において第２検索部３２の動作の手順は、図９を用いて説明した第２位置Ｓ２を検出する場合の第２検索部３２の動作の手順と同じであり、その説明を省略する。左乳房の変換画像において、第２位置Ｓ２を検出する第２検索部３２の動作の順を図１０に（ａ）〜（ｃ）を付して示した。
【００６４】
次に、平滑化処理部３３について図１１を参照して説明する。
平滑化処理部３３は、第１検索部３１及び第２検索部３２により検出された第１位置Ｓ１及び第２位置Ｓ２に対し平滑化をする。検出された第１位置Ｓ１及び第２位置Ｓ２を（Ｘ１００,Ｙ１００）、（Ｘ１０１,Ｙ１０１）、（Ｘ１０２,Ｙ１０２）、…として図１１に示す。例えば、平滑化処理部３３は、検出された第１位置又は第２位置（Ｘ１００,Ｙ１００）について、第１位置又は第２位置の両隣の位置（Ｘ９９,Ｙ９９）、（Ｘ１０１,Ｙ１０１）を基に、（（Ｘ９９＋Ｘ１００＋Ｘ１０１）／３，（Ｙ９９＋Ｙ１００＋Ｙ１０１）／３）の式により求められた平滑化した位置をスキンラインの位置とする。なお、上記平滑化処理において、重み付けを、求めようとする第１位置Ｓ１及び第２位置Ｓ２と、その第１位置Ｓ１又は第２位置Ｓ２の両隣の位置とで異ならせても良い。
【００６５】
システム制御部３９は、マンモグラフィ装置１の中枢として、各構成要素を制御する。
【００６６】
（動作）
次に、一連の動作について図１２から図１６を参照して説明する。図１２は、スキンラインの位置検出を示すフローチャート、図１３は、第１位置の検索がＹ軸上の位置毎に行われることを示すフローチャート、図１４は、第１位置の検索を示すフローチャート、図１５は、第２位置の検索が第１位置の周囲又は直前に検出された第２位置の周囲で行われることを示すフローチャート、図１６は、第２位置の検索を示すフローチャートである。
【００６７】
データ収集部２１は、デジタルマンモグラフィ画像（原画像）を取得する（ステップＳ１０１）。原画像には、被検体を透過したＸ線に基づく生体領域と、被検体を透過していないＸ線に基づく直接線領域とが含まれる。医用画像処理装置４０は、原画像のヒストグラムの曲線を微分することにより、曲線の極小点を算出する。極小点を算出することにより、被検体の生体領域と直接線領域との境となる境界値が特定されている。直接線領域の境界値を図１に示す。
【００６８】
次に、画像変換処理部２７は、デジタルマンモグラフィ画像（原画像）を前述したウェーブレット変換し、ウエーブレット変換後の画像を高周波成分のみで逆変換することにより、変換画像を生成する（ステップＳ１０２）。
【００６９】
次に、閾値算出部２９は、変換画像の全ＸＰｒｏｆｉｌｅの最大値の内の中央値に、予め定められた係数を乗算することにより閾値を算出する（ステップＳ１０３）。
【００７０】
次に、第１検索部３１は、変換画像を、Ｙ軸上の位置をＹ軸の一端側から他端側の方向へ１つずつずらし、ずらしたＹ軸上の位置に対してＸ軸上を検索する場合に、Ｘ軸の一端側から他端側の方向へ検索し、閾値以上の画素値を最初に有する画素の第１位置Ｓ１を検出する（ステップＳ１０４）。
【００７１】
さらに、第１検索部３１により第１位置Ｓ１の検出を図１３を参照して説明する。第１位置Ｓ１の検出は、Ｙ軸上の各位置をＹ軸の一端側の位置Ｙ１から他端側の位置Ｙｎの方向に１つずつずらして行われる。Ｙ軸上の位置をＹｊとし、Ｙ１からＹｎを１からｎとして図１３に示す。
【００７２】
さらに、第１検索部３１がＹ軸上の位置毎に行う第１位置Ｓ１の検出を図１４を参照して説明する。なお、Ｘ軸上の位置をＸｉとし、Ｘ１からＸｎを１からｎとして図１４に示し、また、画素値（輝度値）をＧとして図１４に示す。また、境界値を３とし、閾値を３０とした。
【００７３】
先ず、第１検索部３１は、検索対象の位置Ｘ１の画素値Ｇが３以上であるか否かを判断する（ステップＳ２０１）。第１検索部３１は、位置Ｘ１の画素値Ｇが３未満であると判断した場合（ステップＳ２０１；Ｎｏ）、Ｙ軸上の位置に対する第１位置の検索を終了し、Ｙ軸上で１つずらした位置に対する第１位置を検索する。
【００７４】
第１検索部３１は、位置Ｘ１の画素値Ｇが３以上であると判断した場合（ステップＳ２０１；Ｙｅｓ）、第１検索部３１は、検索対象の位置をＸ軸の他端側（胸壁側）から一端側（直接線領域側）へ１つ移動させ（ｉ＝ｉ＋１）（ステップＳ２０２）、さらに、移動させた検索対象の位置の画素値Ｇが３以上であるか否かを判断する（ステップＳ２０３）。
【００７５】
第１検索部３１は、移動させた検索対象の位置の画素値Ｇが３以上であると判断した場合（ステップＳ２０３；Ｙｅｓ）、検索対象の位置をＸ軸の他端側（胸壁側）から一端側（直接線領域側）へさらに１つ移動させる（ｉ＝ｉ＋１）（ステップＳ２０２）。
【００７６】
第１検索部３１は、位置Ｘｉの画素値Ｇが３未満であると判断した場合（ステップＳ２０３；Ｎｏ）、検索対象の位置をＸ軸の一端側（直接線領域側）から他端側（胸壁側）へ１つ移動させ（ｉ＝ｉ−１）（ステップＳ２０４）、検索対象の位置の画素値Ｇが３０以上であるか否かを判断する（ステップＳ２０５）。
【００７７】
第１検索部３１は、検索対象の位置の画素値Ｇが３０未満であると判断した場合（ステップＳ２０５；Ｎｏ）、検索対象の位置がＸ１であるか否かを判断する（ステップＳ２０７）。検索対象の位置がＸ１でないと判断した場合（ステップＳ２０７；Ｎｏ）、検索対象の位置をＸ軸の一端側（直接線領域側）から他端側（胸壁側）へ１つ移動させる（ｉ＝ｉ−１）ステップＳ２０４に戻る。第１検索部３１は、検索対象の位置がＸ１であると判断した場合（ステップＳ２０７；Ｙｅｓ）、Ｙ軸上の位置に対する第１位置の検索を終了し、Ｙ軸上で１つずらした位置に対する第１位置を検索する。
【００７８】
第１検索部３１は、検索対象の位置の画素値Ｇが３０以上であると判断した場合（ステップＳ２０５；Ｙｅｓ）、閾値３０以上である画素値Ｇを最初に有する画素としてその位置（Ｘｉ，Ｙｉ）を第１位置Ｓ１として、第１検索部３１は、医用画像処理装置４０の内部メモリに記憶させる（ステップＳ２０６）。次に、Ｙ軸上の位置に対する第１位置の検索を終了し、Ｙ軸上で１つずらした位置に対する第１位置を検索する。
【００７９】
以上のような第１位置の検索において、第１検索部３１は、検索対象の位置を先ずＸ軸の他端側（胸壁側）から一端側（直接線領域側）に移動させ、検索対象の位置の画素値が３未満であると判断した後、検索対象の位置を一端側から他端側に移動させたので、検索対象の位置を単に一端側から他端側に移動させた場合に比較して、迅速に第１位置を検出することが可能となる。
【００８０】
次に、第２検索部３２は、第１位置Ｓ１の周囲をＸＹ座標上で検索することにより、閾値３０以上の画素値を最初に有する画素の第２位置Ｓ２を検出し、それ以後は、直前に検出した第２位置Ｓ２の周囲をＸＹ座標上で検索することにより、閾値３０以上の画素値を最初に有する画素の第２位置Ｓ２をＹ軸の他端側（Ｘｎの位置）まで検出していく（ステップＳ１０５）。
【００８１】
さらに、第２検索部３２による第２位置Ｓ２の検索を図１５を参照して説明する。なお、第２検索部３２は、医用画像処理装置４０の内部メモリに記憶された第１位置及び第２位置が無くなるまで、第２位置Ｓ２の検索を続ける。なお、医用画像処理装置４０の内部メモリに記憶された第１位置及び第２位置をＳ１及びＳ２で図１５に示す。
【００８２】
さらに、第２検索部３２が行う第２位置Ｓ２の検出を図１６を参照して説明する。なお、Ｘ軸上の位置をＸｉとし、Ｙ軸上の位置をＹｊとし、また、画素値（輝度値）をＧとして図１６に示す。また、閾値を３０とした。
【００８３】
第２検索部３２は、第１位置Ｓ１又は直前に検出された第２位置Ｓ２を医用画像処理装置４０の内部メモリから読み出す（ステップＳ３０１）。次に、第２検索部３２は、読み出した位置（Ｘｉ，Ｙｊ）に対し、Ｘ軸上の位置（Ｘｉ）を他端側（胸壁側）から一端側（直接線領域側）に移動させ（ｉ＝ｉ＋１）（ステップＳ３０２）、第２検索部３２は、移動させた検索対象の位置（Ｘｉ＋１，Ｙｊ）の画素値Ｇが閾値３０以上であるか否かを判断する（ステップＳ３０３）。
【００８４】
第２検索部３２は、移動させた検索対象の位置（Ｘｉ＋１，Ｙｊ）の画素値Ｇが閾値３０以上であると判断した場合（ステップＳ３０３；Ｙｅｓ）、その位置（Ｘｉ＋１，Ｙｊ）を第２位置として、医用画像処理装置４０の内部メモリに記憶させ（ステップＳ３０４）、その後、終了する。
【００８５】
第２検索部３２は、移動させた検索対象の位置（Ｘｉ＋１，Ｙｊ）の画素値Ｇが閾値３０未満であると判断した場合（ステップＳ３０３；Ｎｏ）、その位置（Ｘｉ＋１，Ｙｊ）をＹ軸の一端側から他端側に１つ移動させ（ｊ＝ｊ＋１）（ステップＳ３０５）、移動させた位置（Ｘｉ＋１，Ｙｊ＋１）の画素値Ｇが閾値３０以上であるか否かを判断する（ステップＳ３０６）。第２検索部３２は、移動させた検索対象の位置（Ｘｉ＋１，Ｙｊ＋１）の画素値Ｇが閾値３０以上であると判断した場合（ステップＳ３０６；Ｙｅｓ）、その位置（Ｘｉ＋１，Ｙｊ＋１）を第２位置として、医用画像処理装置４０の内部メモリに記憶させる（ステップＳ３０４）。
【００８６】
第２検索部３２は、移動させた検索対象の位置（Ｘｉ＋１，Ｙｊ＋１）の画素値Ｇが閾値３０未満であると判断した場合（ステップＳ３０６；Ｎｏ）、その位置（Ｘｉ＋１，Ｙｊ＋１）をＸ軸の一端側から他端側に１つ移動させ（ｉ＝ｉ−１）（ステップＳ３０７）、その位置（Ｘｉ，Ｙｊ＋１）の画素値Ｇが閾値３０以上であるか否かを判断する（ステップＳ３０８）。第２検索部３２は、移動させた検索対象の位置（Ｘｉ，Ｙｊ＋１）の画素値Ｇが閾値３０以上であると判断した場合（ステップＳ３０８；Ｙｅｓ）、その位置（Ｘｉ，Ｙｊ＋１）を第２位置として、医用画像処理装置４０の内部メモリに記憶させる（ステップＳ３０４）。
【００８７】
第２検索部３２は、移動させた検索対象の位置（Ｘｉ，Ｙｊ＋１）の画素値Ｇが閾値３０未満であると判断した場合（ステップＳ３０８；Ｎｏ）、その位置（Ｘｉ，Ｙｊ＋１）をＸ軸の一端側から他端側に１つ移動させ（ｉ＝ｉ−１）（ステップＳ３０９）、その位置（Ｘｉ−１，Ｙｊ＋１）がＸ軸の他端側及び他端側より一端側である場合（ｉ≧１）（ステップＳ３１０；Ｙｅｓ）、その位置（Ｘｉ−１，Ｙｊ＋１）の画素値Ｇが閾値３０以上であるか否かを判断する（ステップＳ３１１）。その位置（Ｘｉ−１，Ｙｊ＋１）がＸ軸の他端側よりマイナスである場合（ｉ＜１）（ステップＳ３１０；Ｎｏ）、終了する。
【００８８】
第２検索部３２は、移動させた検索対象の位置の画素値Ｇが閾値３０以上であると判断した場合（ステップＳ３１１；Ｙｅｓ）、その位置（Ｘｉ−１，Ｙｊ＋１）を第２位置として、医用画像処理装置４０の内部メモリに記憶させる（ステップＳ３０４）。
【００８９】
第２検索部３２は、移動させた検索対象の位置（Ｘｉ−１，Ｙｊ＋１）の画素値Ｇが閾値３０未満であると判断した場合（ステップＳ３１１；Ｎｏ）、その位置（Ｘｉ−１，Ｙｊ＋１）をＹ軸の他端側から一端側に１つ移動させ（ｊ＝ｊ−１）（ステップＳ３１２）、その位置（Ｘｉ−１，Ｙｊ）の画素値Ｇが閾値３０以上であるか否かを判断する（ステップＳ３１３）。第２検索部３２は、移動させた検索対象の位置（Ｘｉ−１，Ｙｊ）の画素値Ｇが閾値３０以上であると判断した場合（ステップＳ３１３；Ｙｅｓ）、その位置（Ｘｉ−１，Ｙｊ）を第２位置として、医用画像処理装置４０の内部メモリに記憶させる（ステップＳ３０４）。第２検索部３２は、移動させた検索対象の位置（Ｘｉ−１，Ｙｊ）の画素値Ｇが閾値３０未満であると判断した場合（ステップＳ３１３；Ｎｏ）、終了する。
【００９０】
以上のように、第２検索部３２は、第１位置Ｓ１又は直前に検出された第２位置Ｓ２に対し、検索対象となる位置を、次の順番で移動させる。すなわち、先ず、Ｘ軸上の位置（Ｘｉ）を他端側（胸壁側）から一端側（直接線領域側）に１つ移動させる（ｉ＝ｉ＋１）（ステップＳ３０２）。その後、検索対象となる位置をＹ軸の一端側（上側）から他端側（下側）に１つ移動させる（ｊ＝ｊ＋１）（ステップＳ３０５）。その後、検索対象となる位置を一端側から他端側に１つ移動させる。さらに、検索対象となる位置を一端側から他端側に１つ移動させる。その後、検索対象となる位置をＹ軸の他端側（下側）から一端側（上側）に１つ移動させる。
【００９１】
したがって、検索対象となる位置をＹ軸の一端側（上側）から他端側（下側）に移動する場合において、胸壁側（他端側）を右端側とした右乳房の変換画像の検索においては、第２検索部３２は検索対象となる位置を反時計方向に移動させる。また、胸壁側（他端側）を左端側とした左乳房の変換画像の検索においては、第２検索部３２は検索対象となる位置を時計方向に移動させる。つまり、検索対象となる位置をＹ軸の一端側から他端側に移動する場合において、第２検索部３２による検索する回転方向は、右乳房の変換画像と左乳房の変換画像とでは互いに反対の方向となる。
【００９２】
以上のようにして、第２検索部３２は、第２位置Ｓ２をＹ軸の他端側の位置（Ｙｎ）まで、検索する。
【００９３】
次に、平滑化処理部３３は、前述したように、第１検索部３１及び第２検索部３２により検出された第１位置Ｓ１及び第２位置Ｓ２に対し平滑化をする（ステップＳ１０６）。平滑化により、スキンラインの位置精度を向上することが可能となる。
【００９４】
なお、上記実施形態では、デジタルマンモグラフィ画像（原画像）をウェーブレット変換し、ウエーブレット変換後の画像を高周波成分のみで逆変換することにより、変換画像を生成する画像変換処理部２７を示したが、これに限らず、デジタルマンモグラフィ画像（原画像）中で、直接線領域のみの画素値を予め定められた値以下の画素値に変換することにより、変換画像を生成する画像変換処理部２７であっても良い。
【符号の説明】
【００９５】
１  マンモグラフィ装置    ３  Ｘ線撮影台    ５  コンピュータ装置
１０  基台部    １１  Ｃアーム    １５  Ｘ線発生装置    １６  Ｘ線検出器
１７  圧迫ユニット    １８  Ｘ線管    １９  高電圧発生器
２１  データ収集部    ２３  画像発生部    ２５  画像記憶部
２７  画像変換処理部    ２９  閾値算出部    ３１  第１検索部
３２  第２検索部    ３３  平滑化処理部    ３５  操作部    ３７  撮影制御部
３９  システム制御部    ４０  医用画像処理装置
【図１】
0005457812.tif 000002
【図４】
0005457812.tif 000003
【図６】
0005457812.tif 000004
【図１１】
0005457812.tif 000005
【図１２】
0005457812.tif 000006
【図１３】
0005457812.tif 000007
【図１４】
0005457812.tif 000008
【図１５】
0005457812.tif 000009
【図２】
0005457812.tif 000010
【図３】
0005457812.tif 000011
【図５】
0005457812.tif 000012
【図７】
0005457812.tif 000013
【図８】
0005457812.tif 000014
【図９】
0005457812.tif 000015
【図１０】
0005457812.tif 000016
【図１６】
0005457812.tif 000017