特開 2024-76773 乳腺密度の算出方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公開特許公報(A)

(11)【公開番号】P2024076773

(43)【公開日】2024-06-06

(54)【発明の名称】乳腺密度の算出方法

(51)【国際特許分類】

   A61B 6/00 20240101AFI20240530BHJP

【ＦＩ】

A61B6/00 330Z

A61B6/00 350Z

【審査請求】未請求

【請求項の数】4

【出願形態】ＯＬ

(21)【出願番号】P 2022188516

(22)【出願日】2022-11-25

【新規性喪失の例外の表示】特許法第３０条第２項適用申請有り 集会名 ：第３１回日本乳癌検診学会学術総会 開催日 ：２０２１年１１月２６日（２０２１年１１月２６日・２７日） 〔刊行物等〕 配信年月日：２０２１年１２月１０日～２４日（オンデマンド配信） 配信アドレス：ｈｔｔｐｓ：／／ｃｏｎｖｅｎｔｉｏｎ．ｊｔｂｃｏｍ．ｃｏ．ｊｐ／ｊａｂｃｓ２０２１／

(71)【出願人】

【識別番号】000125347

【氏名又は名称】学校法人近畿大学

(74)【代理人】

【識別番号】100118924

【弁理士】

【氏名又は名称】廣幸 正樹

(72)【発明者】

【氏名】浅井 義行

(72)【発明者】

【氏名】山室 美佳

【テーマコード（参考）】

4C093

【Ｆターム（参考）】

4C093AA01

4C093AA07

4C093CA18

4C093DA06

4C093FD03

4C093FD09

4C093FD11

4C093FF16

4C093FF28

4C093FG13

(57)【要約】

【課題】マンモグラフィの画像から乳腺密度を算出する際に基準値となる脂肪領域の画素値を求めるのは容易でなかった。
【解決手段】検査乳房厚（Ｄｍｔ）と検査Ｘ線量（ＤＲｍｔ）と検査乳房Ｘ線画像を得る工程と、
予め複数の乳房Ｘ線サンプル撮影より得られる乳房厚（ｍｔ）とＸ線量（Ｒｍｔ）と脂肪画素値（Ｙｍｔ）の情報から、前記脂肪画素値（Ｙｍｔ）を前記Ｘ線量（Ｒｍｔ）で除して得られた単位Ｘ線量あたりの単位脂肪画素値（Ｙｍｔ／Ｒｍｔ＝Ｚｍｔ）と前記乳房厚（ｍｔ）を平準化することで求めた規格化乳房厚脂肪画素値（ＺＥｍｔ）のうち前記検査乳房厚（Ｄｍｔ）に対応する前記規格化乳房厚脂肪画素値（ＺＥｍｔ）を前記検査Ｘ線量（ＤＲｍｔ）に乗じて検査乳房厚脂肪画素値（ＤＹ）を求める工程と
前記検査乳房厚脂肪画素値（ＤＹ）を用いて乳腺密度を求める工程を有する乳腺密度の算出方法は、妥当な脂肪画素値を容易に求めることができる。
【選択図】図１

【特許請求の範囲】

【請求項1】

乳房Ｘ線撮影装置で得られるマンモグラフィ画像から乳腺密度を算出する際に用いられる検査乳房厚脂肪画素値（ＤＹ）を前記マンモグラフィ画像の撮影時に得られる検査乳房厚（Ｄｍｔ）に対応し、前記撮影時の検査Ｘ線量（ＤＲｍｔ）に乗ずることで前記検査脂肪画素値（ＤＹ）を算出できる規格化乳房厚脂肪画素値（ＺＥｍｔ）の算出方法であって、
乳房厚（ｍｔ）とＸ線量（Ｒｍｔ）と脂肪画素値（Ｙｍｔ）の情報を有する複数の乳房Ｘ線サンプル画像データに対して前記脂肪画素値（Ｙｍｔ）を前記Ｘ線量（Ｒｍｔ）で除し、単位Ｘ線量あたりの脂肪画素値（Ｙｍｔ／Ｒｍｔ＝Ｚｍｔ）を算出する工程と、
前記単位Ｘ線量あたりの脂肪画素値（Ｚｍｔ）と前記乳房厚（ｍｔ）を平準化し、前記乳房厚（ｍｔ）毎の前記規格化乳房厚脂肪画素値（ＺＥｍｔ）を算出する工程を有する規格化乳房厚脂肪画素値の算出方法。

【請求項2】

乳房Ｘ線撮影時の検査乳房厚（Ｄｍｔ）と検査Ｘ線量（ＤＲｍｔ）と検査乳房Ｘ線画像を得る工程と、
予め複数の乳房Ｘ線サンプル撮影より得られる乳房厚（ｍｔ）とＸ線量（Ｒｍｔ）と脂肪画素値（Ｙｍｔ）の情報から、前記脂肪画素値（Ｙｍｔ）を前記Ｘ線量（Ｒｍｔ）で除して得られた単位Ｘ線量あたりの脂肪画素値（Ｙｍｔ／Ｒｍｔ＝Ｚｍｔ）と前記乳房厚（ｍｔ）を平準化することで求めた規格化乳房厚脂肪画素値（ＺＥｍｔ）のうち前記検査乳房厚（Ｄｍｔ）に対応する前記規格化乳房厚脂肪画素値（ＺＥｍｔ）を前記検査Ｘ線量（ＤＲｍｔ）に乗じて検査乳房厚脂肪画素値（ＤＹ）を求める工程と
前記検査乳房厚脂肪画素値（ＤＹ）を用いて乳腺密度を求める工程を有する乳腺密度の算出方法。

【請求項3】

前記乳腺密度を求める工程は、
前記検査乳房Ｘ線画像において前記検査乳房厚脂肪画素値（ＤＹ）より胸筋側の領域で乳腺密度を測定する工程である請求項２に記載された乳腺密度の算出方法。

【請求項4】

少なくとも圧迫板で挟んだ乳房をＸ線撮影し、前記乳房のＸ線画像を表示する表示装置を持った乳房Ｘ線撮影装置であって、
予め複数の乳房Ｘ線サンプル撮影より得られる乳房厚（ｍｔ）とＸ線量（Ｒｍｔ）と脂肪画素値（Ｙｍｔ）の情報から、前記脂肪画素値（Ｙｍｔ）を前記Ｘ線量（Ｒｍｔ）で除して得られた単位Ｘ線量あたりの脂肪画素値（Ｙｍｔ／Ｒｍｔ＝Ｚｍｔ）と前記乳房厚（ｍｔ）を平準化することで求めた規格化乳房厚脂肪画素値（ＺＥｍｔ）と前記乳房厚（ｍｔ）を対応付ける参照テーブルと、
撮影された乳房のデータのうち前記検査乳房厚（Ｄｍｔ）に対応する前記規格化乳房厚脂肪画素値（ＺＥｍｔ）を前記検査Ｘ線量（ＤＲｍｔ）に乗じて算出した検査乳房厚脂肪画素値（ＤＹ）を算出し、
前記撮影された乳房のＸ線画像上に前記検査乳房厚脂肪画素値（ＤＹ）を有する画素のうち、前記検査乳房厚脂肪画素値（ＤＹ）より画素値の大きな画素と隣接する画素をマークする制御部を有する乳房Ｘ線撮影装置。

【発明の詳細な説明】

【技術分野】

【0001】

本発明はマンモグラフィで得られる乳房のＸ線画像から乳腺密度を算出する方法に関する。

【背景技術】

【0002】

乳癌は乳腺密度と関連するといわれており、乳腺密度を観測できる乳房Ｘ線装置（マンモグラフィ）が利用されている。マンモグラフィ画像から３次元乳腺密度（乳腺密度と呼ぶ）を計測するためには画像に含まれる脂肪領域と乳腺領域の画素値が要求される。このうち、基準となるのは脂肪領域の画素値（以後「脂肪画素値」と呼ぶ）であり、その精度が乳腺密度計測の精度を左右する重要な要素となる。

【0003】

特許文献１には、Ｘ線画像を所定の表示部において表示する際の表示条件として、当該Ｘ線画像を構成する画素の画素値のうち前記所定の表示部にて表示される画素値の範囲と、当該画素値の範囲における中央値とを、被検体の乳房において推定された乳腺密度に応じて算出する表示条件算出手段を有する画像表示装置が開示されている。

【先行技術文献】

【特許文献】

【0004】

【特許文献1】特開２００９－２６８７２６号公報

【発明の概要】

【発明が解決しようとする課題】

【0005】

特許文献１は撮影したＸ線画像において乳房部分を常にコントラスト良く表示する手法であるが、却ってＸ線画像毎の脂肪画素値がわかりにくくなり、乳腺密度を数値的に求め比較しようとすると容易でないという課題を有する。

【0006】

現在、市販の乳腺密度自動計測ソフトは基準値の算出法及び算出値を開示しておらず、独自の値を設定して乳腺密度値だけを表示している。そのため、同一のマンモグラフィ画像の乳腺密度を異なる自動解析ソフトで計測すると３～４倍も異なる結果が得られており、臨床の場においてどの値が正しいのか判断できず、乳腺密度と関連する乳癌罹患の危険性や病変見落としの危険性を患者へ詳しく説明することができていない。

【課題を解決するための手段】

【0007】

本発明は上記の課題に鑑みて想到されたものであり、異なる条件で撮影されたＸ線画像であっても、精度よく脂肪画素値を得ることができ、安定した乳腺密度を算出するための係数となる規格化乳房厚脂肪画素値を与える方法を提供する。

【0008】

より具体的に本発明に係る規格化乳房厚脂肪画素値の算出方法は、
乳房Ｘ線撮影装置で得られるマンモグラフィ画像から乳腺密度を算出する際に用いられる検査乳房厚脂肪画素値（ＤＹ）を前記マンモグラフィ画像の撮影時に得られる検査乳房厚（Ｄｍｔ）に対応し、前記撮影時の検査Ｘ線量（ＤＲｍｔ）に乗ずることで前記検査乳房厚脂肪画素値（ＤＹ）を算出できる規格化乳房厚脂肪画素値（ＺＥｍｔ）の算出方法であって、
乳房厚（ｍｔ）とＸ線量（Ｒｍｔ）と脂肪画素値（Ｙｍｔ）の情報を有する複数の乳房Ｘ線サンプル画像データに対して前記脂肪画素値（Ｙｍｔ）を前記Ｘ線量（Ｒｍｔ）で除し、単位Ｘ線量あたりの脂肪画素値（Ｙｍｔ／Ｒｍｔ＝Ｚｍｔ）を算出する工程と、
前記単位Ｘ線量あたりの単位脂肪画素値（Ｚｍｔ）と前記乳房厚（ｍｔ）を平準化し、前記乳房厚（ｍｔ）毎の前記規格化乳房厚脂肪画素値（ＺＥｍｔ）を算出する工程を有することを特徴とする。

【0009】

また、この規格化乳房厚脂肪画素値を用いた乳腺密度算出方法も提供する。

【0010】

より具体的に本発明に係る乳腺密度算出方法は、
乳房Ｘ線撮影時の検査乳房厚（Ｄｍｔ）と検査Ｘ線量（ＤＲｍｔ）と検査乳房Ｘ線画像を得る工程と、
予め複数の乳房Ｘ線サンプル撮影より得られる乳房厚（ｍｔ）とＸ線量（Ｒｍｔ）と脂肪画素値（Ｙｍｔ）の情報から、前記脂肪画素値（Ｙｍｔ）を前記Ｘ線量（Ｒｍｔ）で除して得られた単位Ｘ線量あたりの脂肪画素値（Ｙｍｔ／Ｒｍｔ＝Ｚｍｔ）と前記乳房厚（ｍｔ）を平準化することで求めた規格化乳房厚脂肪画素値（ＺＥｍｔ）のうち前記検査乳房厚（Ｄｍｔ）に対応する前記規格化乳房厚脂肪画素値（ＺＥｍｔ）を前記検査Ｘ線量（ＤＲｍｔ）に乗じて検査乳房厚脂肪画素値（ＤＹ）を求める工程と
前記検査乳房厚脂肪画素値（ＤＹ）を用いて乳腺密度を求める工程を有することを特徴とする。

【発明の効果】

【0011】

本発明に係る乳腺密度算出方法は、乳房厚毎に求められた撮影時のＸ線量当たりの脂肪画素値である規格化乳房厚脂肪画素値を用いるので、実際の撮影時の乳房厚（これを検査乳房厚と呼ぶ。）に対応する規格化乳房厚脂肪画素値をＸ線撮影時のＸ線量（これを検査Ｘ線量と呼ぶ。）に乗じることで、撮影されたＸ線画像における妥当な脂肪画素値（これを検査脂肪画素値ＤＹｍｔ若しくは単にＤＹと呼ぶ。）を得ることができ、安定して乳腺密度を算出することができる。

【図面の簡単な説明】

【0012】

【図1】マンモグラフィの撮影の概略を示す図である。

【図2】マンモグラフィデータを示す図である。

【図3】乳房厚ｍｔと単位脂肪画素値Ｚｍｔを近似した式を規格化乳房厚脂肪画素式１００を示す図である。

【図4】乳房厚ｍｔと単位脂肪画素値Ｚｍｔの実測値例を示すグラフである。

【図5】参照テーブルを例示する図である。

【図6】検査脂肪画素値の算出方法を示す図である。

【図7】乳腺密度を測定する処理を示すフロー図である。

【図8】圧迫不足で乳腺密度の算出から削除する部分を示す図である。

【発明を実施するための形態】

【0013】

以下に本発明に係る乳腺密度の算出方法について図面を用いて説明を行う。なお、以下の説明は、本発明の一実施形態および一実施例を例示するものであり、本発明が以下の説明に限定されるものではない。以下の説明は本発明の趣旨を逸脱しない範囲で改変することができる。

【0014】

本発明に係る乳腺密度の算出方法では、いかなる条件でのＸ線撮影画像であっても、同一厚みの脂肪領域の単位Ｘ線量当たりの画素値は同じであるという発想に基づく。

【0015】

図１にマンモグラフィの撮影の概略図を示す。図１（ａ）は、典型的に撮影される乳房のＸ線画像１０を示す。図１（ｂ）は撮影時の様子を示す。図１（ｂ）を参照して、乳房１２は、被験者の左右から押し当てられる圧迫板３０で平たく伸ばされ、左右方向の何れか一方からＸ線が照射され、他方向の受光部３２で画像化される。

【0016】

この時にＸ線画像１０の他、撮影時の乳房厚（以後「検査乳房厚Ｄｍｔ」と呼ぶ。）、撮影時のＸ線量（以後「検査Ｘ線量ＤＲｍｔ」と呼ぶ。単位はｍＡｓであり、ｍＡｓ値とも呼ばれる。）も同時に計測される。すなわち、マンモグラフィの撮影では、少なくともＸ線画像１０と検査乳房厚Ｄｍｔと検査Ｘ線量ＤＲｍｔは同一被験者からのデータとして一括に扱われる。これをマンモグラフィデータ２０と呼ぶ。またＸ線画像１０はデジタルデータとして取得され、各画素での値（画素値）はデジタルデータとして扱うことができる。

【0017】

乳房１２は、脂肪領域１４と、乳腺領域１６とニップル１８で構成される。乳腺密度は、脂肪領域１４の、面積当たりの重量と乳腺領域１６の面積当たりの質量の比で表される。これには、乳腺領域（乳腺組織）１６の厚さおよび脂肪領域１４の厚さが必要となる。Ｘ線画像１０からこれらの厚みを算出する手法として（１）式のＥｎｇｅｌａｎｄの式がよく知られている。

【0018】

【数1】

Ｔ_ｇｌａ：乳腺組織の厚さ
Ｐ_ｇｌａ：乳腺領域の画素値
Ｐ_ａｄｉ：脂肪領域の画素値
μ（ｅｆｆ）：脂肪および乳腺組織の実効線減弱係数
Ｎ：撮影管電圧での全光子数
ｎ：光子エネルギー毎の光子数
Ｅ：光子エネルギー

【0019】

（２）式は、（１）式の分母を算出するための式で、脂肪領域１４および乳腺領域１６がそれぞれどれだけのＸ線を減弱させるかを表す。この部分は乳房Ｘ線撮影装置毎で決まる値となる。

【0020】

（１）式において、Ｐ_ａｄｉが脂肪領域１４の画素値に当たる。この脂肪領域１４の画素値を検査乳房厚「Ｄｍｔ」の時の脂肪画素値という意味で以後「検査脂肪画素値ＤＹｍｔ」と呼ぶ。また、「撮影で得られたマンモグラフィ画像の脂肪画素値ＤＹ」もしくは「検査乳房厚脂肪画素値ＤＹ」と呼んでもよい。

【0021】

検査脂肪画素値ＤＹｍｔは従来「閾値法」などの公知の方法で、Ｘ線画像１０毎にエキスパートがＸ線画像１０の画素値を特定表示させながら、乳腺領域１６の形状分布などを考慮し逐次マニュアルで決められていた。本発明は、検査乳房厚Ｄｍｔ毎に求めた規格化乳房厚脂肪画素値ＺＥｍｔを検査Ｘ線量ＤＲｍｔに乗ずるだけで検査脂肪画素値ＤＹｍｔすなわち撮影で得られたマンモグラフィ画像の脂肪画素値ＤＹを算出する方法を提供する。

【0022】

規格化乳房厚脂肪画素値ＺＥｍｔは以下のように求める。図２を参照して、まず、多くのサンプル用マンモグラフィデータ２４を収集する。これらは乳房Ｘ線サンプル撮影によって得られるが、通常のマンモグラフィの撮影であってよい。この際（ｉ）脂肪の含有量が顕著に少ないマンモグラフィ画像、（ｉｉ）ポジショニングが適切に行われていないマンモグラフィ画像、（ｉｉｉ）乳房の形状が標準的でないケース、においては脂肪画素値Ｙｍｔを求める際の誤差が大きくなる場合があるので、削除するのが望ましい。

【0023】

サンプル用マンモグラフィデータ２４には、Ｘ線画像ＸＧ、乳房厚ｍｔ、Ｘ線量Ｒｍｔ、脂肪画素値Ｙｍｔが含まれる。すなわち、エキスパートによる脂肪画素値Ｙｍｔが付与されている。サンプル用マンモグラフィデータ２４はｎ個あるとして、それぞれのパラメータの右下に添え字で示す。例えば、１番目のサンプル用マンモグラフィデータ２４には、Ｘ線画像ＸＧ_１、乳房厚ｍｔ_１、Ｘ線量Ｒｍｔ_１、脂肪画素値Ｙｍｔ_１が含まれる。

【0024】

また、これらのサンプル用マンモグラフィデータ２４の集合をマンモグラフィデータ母集団２６と呼ぶ。マンモグラフィデータ母集団２６は同一機種の乳房Ｘ線撮影装置のデータが好ましい。機種が異なると、Ｘ線量の設定や受光部３２の構成や感度が異なる場合があるからである。

【0025】

図２（ｂ）には、乳房厚ｍｔ、Ｘ線量Ｒｍｔ、脂肪画素値Ｙｍｔをまとめたテーブル形式で示す。さらに脂肪画素値ＹｍｔをＸ線量Ｒｍｔで除したものを示す。これは単位Ｘ線量あたりの脂肪画素値であり、単位脂肪画素値Ｚｍｔと呼ぶ。乳房厚ｍｔは図１（ｂ）で示す様に、撮影時には乳房を一定の圧力で押さえつけており、乳腺領域１６のない脂肪領域１４の厚みとほぼ考えられる。したがって、いかなる条件でのＸ線撮影画像であっても、同一厚みの脂肪領域の単位Ｘ線量当たりの画素値は同じであるならば、単位脂肪画素値Ｚｍｔは、乳房厚ｍｔと比例関係にあるはずである。

【0026】

実際は、使用するＸ線のエネルギーに幅があり、また、脂肪画素値Ｙｍｔの決定にも誤差が入る。そこで、乳房厚ｍｔと単位脂肪画素値Ｚｍｔを最小二乗法により平準化する。このように乳房厚ｍｔと単位脂肪画素値Ｚｍｔを近似した式を規格化乳房厚脂肪画素式１００と呼ぶ。

【0027】

図３にその様子を示す。図３において、横軸は乳房厚ｍｔであり、縦軸は単位脂肪画素値Ｚｍｔである。黒丸は、乳房厚ｍｔ毎の単位脂肪画素値Ｚｍｔの平均を表す。単位脂肪画素値Ｚｍｔは、乳房厚ｍｔに対して明らかに不連続となる点ＤＣＰを有する。これは、乳房厚ｍｔの変化に伴い撮影装置の管電流のレンジが切り替わることで生じると考えられる。このような場合は、不連続点ＤＣＰを除いた領域に分け、領域ごとに平準化を行う。図３では、領域ＡＭ１とＡＭ２とＡＭ３である。なお、不連続点ＤＣＰはそれ毎に扱ってよいし、それ以降の領域の先頭データとして扱ってもよい。

【0028】

図４には、ある乳房Ｘ線撮影装置において、サンプル用マンモグラフィデータ２４が３００枚のマンモグラフィデータ母集団２６を用いて乳房厚ｍｔと単位脂肪画素値Ｚｍｔをプロットしたものである。各黒丸はその乳房厚ｍｔにおける平均値である。この装置では、乳房厚ｍｔが１０ｍｍから８０ｍｍの間に７つの不連続点ＤＣＰがあり、７つの領域に分けて平準化を行う。

【0029】

再度図３を参照し、乳房厚ｍｔと単位脂肪画素値Ｚｍｔを近似する近似式の型式は特に限定されないが、指数関数や多項式が好適に利用できる。図３では、乳房厚ｍｔの２次多項式で近似した様子を示す。言い換えれば図３では、規格化乳房厚脂肪画素式１００を乳房厚ｍｔの２次多項式として表した。図３では、乳房厚ｍｔの３つの領域ごとに規格化乳房厚脂肪画素式１００が求められている。

【0030】

この多項式は、乳房厚ｍｔを与えると、単位脂肪画素値Ｚｍｔの推定値を算出する式である。この推定値を規格化乳房厚脂肪画素値ＺＥｍｔと呼ぶ。この規格化乳房厚脂肪画素値ＺＥｍｔは、その乳房Ｘ線撮影装置において、「同一厚みの脂肪領域の単位Ｘ線量当たりの画素値は同じ」という思想から求められたもので、ある検査乳房厚Ｄｍｔに対して検査Ｘ線量ＤＲｍｔを照射したとすれば、ＤＲｍｔ×ＺＥｍｔが検査脂肪画素値ＤＹｍｔであり、撮影で得られたマンモグラフィ画像の脂肪画素値ＤＹとなる。

【0031】

実際の乳房Ｘ線撮影装置では、規格化乳房厚脂肪画素式１００を内蔵しておき、その都度規格化乳房厚脂肪画素値ＺＥｍｔを算出してもよいが、参照テーブルにしておくとアクセスも早い。特に不連続点ＤＣＰが存在する場合は参照テーブルにして保持するのが望ましい。図５には、この参照テーブルを例示する。図５（ａ）は図２（ｂ）の情報も含めた参照テーブルＴ１である。参照テーブルとしては、図５（ｂ）のように乳房厚ｍｔと規格化乳房厚脂肪画素値ＺＥｍｔを有した参照テーブルＴ２であれば足りる。

【0032】

図６に本発明による検査脂肪画素値ＤＹｍｔの算出方法を示す。新たなマンモグラフィデータ２０が得られたとする。ここにはＸ線画像１０と検査乳房厚Ｄｍｔと検査Ｘ線量ＤＲｍｔが含まれる。まず検査乳房厚Ｄｍｔと同じ乳房厚ｍｔ_ｉを参照テーブルＴ２から探す（Ｐ１）。そして、乳房厚ｍｔ_ｉに対応する規格化乳房厚脂肪画素値ＺＥｍｔ_ｉを参照テーブルＴ２から選出する。そして、撮影で得たＸ線画像１０の検査脂肪画素値ＤＹｍｔを検査Ｘ線量ＤＲｍｔと規格化乳房厚脂肪画素値ＺＥｍｔ_ｉの積として求める（Ｐ２）。これが撮影で得られたマンモグラフィ画像の脂肪画素値ＤＹである。

【0033】

また、規格化乳房厚脂肪画素値ＺＥｍｔ_ｉを規格化乳房厚脂肪画素式１００から求めてもよい。Ｐ１１には、規格化乳房厚脂肪画素式１００に検査乳房厚Ｄｍｔを代入して規格化乳房厚脂肪画素値ＺＥｍｔ＊を算出し、Ｐ１２でこれに検査Ｘ線量ＤＲｍｔを乗ずることで検査脂肪画素値ＤＹｍｔを求めるプロセスを示した。なお、不連続点ＤＣＰが規格化乳房厚脂肪画素式１００に組み込まれていない場合は、その時の検査乳房厚Ｄｍｔについては、別途参照テーブルから引いてくる。

【0034】

以上にようにして検査脂肪画素値ＤＹｍｔ（検査乳房厚脂肪画素値ＤＹ）を求め、これを（１）式のＰ_ａｄｉとすることで、乳腺密度を求めることができる。

【0035】

図７には、（１）式のＥｎｇｅｌａｎｄの式を用いた乳腺密度の具体的な算出フローの一例を示す。乳腺密度ＭａＭｄを算出する際には、対象となるマンモグラフィデータ２０が得られている事、規格化乳房厚脂肪画素値ＺＥｍｔが求められている事、撮影した乳房Ｘ線撮影装置の設定から（２）式について求められている事が必要である。また、規格化乳房厚脂肪画素値ＺＥｍｔは、図５（ｂ）に示した参照テーブルＴ２の状態で算出されているものとする。

【0036】

図７を参照して、処理が開始されると（ステップＳ１００）、脂肪領域の画素値Ｐａｄｉに撮影されたマンモグラフィデータ２０の検査乳房厚Ｄｍｔに対応する乳房厚ｍｔの規格化乳房厚脂肪画素値ＺＥｍｔが入力される（ステップＳ１０２）。

【0037】

次に乳腺密度を算出したい領域の画素値（乳腺領域１６の画素値）を入力する（ステップＳ１０４）。実際は乳腺領域１６と考えられる一定の領域を指定し、その領域の画素値の平均値を入力するようにしてもよい。

【0038】

次に（１）式のＥｎｇｅｌａｎｄの式を用いて乳腺組織の厚さＴ_ｇｌａ（ｍｍ）を算出する（ステップＳ１０６）。そして、脂肪組織の厚さＴ_ａｄｉ（ｍｍ）を乳房厚ｍｔから乳腺組織の厚さＴ_ｇｌａを引くことで求める（ステップＳ１０８）。

【0039】

乳腺密度ＭａＭｄは、乳房に対する乳腺組織の質量比として表される。乳腺組織と脂肪組織の比重を、それらをρ_ａｄｉ（ｇ／ｍｍ^３）、ρ_ｇｌａ（ｇ／ｍｍ^３）とすれば、Ｔ_ｇｌａ＊ρ_ｇｌａとＴ_ａｄｉ＊ρ_ａｄｉ、はそれぞれ乳腺組織の面積質量（ｇ／ｍｍ^２）と脂肪組織の面積質量（ｇ／ｍｍ^２）となる。従って、乳腺密度はこれらの比として算出される（ステップＳ１１０）。その後処理は停止する（ステップＳ１２０）が、他の処理に移ることも可能である。

【0040】

このフローは乳房Ｘ線撮影装置の制御装置で行われてもよい。またこのフローを実行できるプログラムとして提供されることもできる。乳房Ｘ線撮影装置としては、図５の参照テーブルＴ２と上記のフローを実行する制御装置を有していれば本発明を実施する乳房Ｘ線撮影装置と言える。

【0041】

図８には、再度Ｘ線画像１０（図８（ａ））と圧迫板３０で圧迫した乳房１２（図８（ｂ））を示した。図８（ｂ）を参照して、領域Ｓ２は圧迫板３０で両側から圧迫している領域で、しかも脂肪だけの領域である。この領域Ｓ２の画素値は検査脂肪画素値ＤＹｍｔとして求めることができる。一方領域Ｓ１は圧迫板３０では押しきれていない部分である。すなわち、脂肪画素値Ｙｍｔを正しく適用できない範囲であるので、領域Ｓ１は乳腺密度の算出としては、省くべき領域である。

【0042】

領域Ｓ１と領域Ｓ２の境界ＳＰは、検査脂肪画素値ＤＹｍｔを境界値として識別することができる。つまり、検査脂肪画素値ＤＹｍｔより画素値が高い部分が領域Ｓ１となる。このように、本発明によって与えられる規格化乳房厚脂肪画素値ＺＥｍｔを用いて検査脂肪画素値ＤＹｍｔを算出すると、乳腺密度を算出するのに好適な領域Ｓ２を得ることができる。

【0043】

図８（ａ）のＸ線画像１０でこの境界ＳＰを例示した。つまり、検査脂肪画素値ＤＹｍｔで作られるニップル側の境界ＳＰよりニップル側の領域Ｓ１は乳腺密度を算出する際には除外する。言い換えると、境界ＳＰより胸筋側の領域が乳腺密度を測定するのに好適な領域である。

【0044】

すなわち、少なくとも圧迫板で挟んだ乳房をＸ線撮影し、そのＸ線画像を表示する表示装置を持った乳房Ｘ線撮影装置において、表示された乳房のＸ線画像に本発明によって得られる規格化乳房厚脂肪画素値ＺＥｔｍを用いた検査脂肪画素値ＤＹｍｔと、検査脂肪画素値ＤＹｍｔより画素値の大きな領域との境界線を表示できれば、乳腺密度の算出の際に誤った部分を指定することを排除できる。

【0045】

具体的に境界ＳＰは、検査乳房厚脂肪画素値（ＤＹ）を有する画素のうち、検査乳房厚脂肪画素値（ＤＹ）より画素値の大きな画素と隣接する画素をマーキングすることで、得ることができる。

【産業上の利用可能性】

【0046】

本発明は乳腺密度の算出ソフトだけでなく乳腺Ｘ線撮影装置にも好適に利用できる。

【符号の説明】

【0047】

１０ Ｘ線画像
１２ 乳房
１４ 脂肪領域
１６ 乳腺領域
１８ ニップル
２０ マンモグラフィデータ
３０ 圧迫板
３２ 受光部
１００ 規格化乳房厚脂肪画素式
Ｄｍｔ 検査乳房厚
ＤＲｍｔ 検査Ｘ線量
ＤＹｍｔ 検査脂肪画素値
ＺＥｍｔ 規格化乳房厚脂肪画素値
２４ サンプル用マンモグラフィデータ
ｍｔ 乳房厚
Ｒｍｔ Ｘ線量
Ｙｍｔ 脂肪画素値
２６ マンモグラフィデータ母集団
Ｚｍｔ 単位脂肪画素値
ＤＣＰ 不連続点

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】