(57)【要約】【課題】処理速度を速くできるだけでなく、境界外の領域を示す電気信号が境界内の領域を示す電気信号の処理に影響を及ぼす範囲を削除しない場合にその範囲を縮小することもできる。【解決手段】画像処理すべき物体（１３０）を、Ｘ線源（１００）とＸ線画像装置（１４０）との間に配置する。このＸ線源（１００）を作動させてＸ線（１２０）を前記物体（１３０）の中及び周囲に向ける。前記Ｘ線画像装置（１４０）に衝突する前記Ｘ線（１２０）の強度を示す電気信号を用いて、前記物体（１３０）の境界を捜し出し、この境界を用いてＸ線画像を生成する。代わりに、又は追加して、電気信号のセット内で最高度数の大きさを用いてＸ線画像を生成する。【選択図】図１
2007313360.tif 000002
【特許請求の範囲】
【請求項１】
Ｘ線源とＸ線画像装置との間に物体を配置するステップと、 前記Ｘ線源を作動させてＸ線を前記Ｘ線画像装置まで前記物体の中及び周囲に向け、これにより、前記Ｘ線画像装置に衝突するＸ線の強度が前記Ｘ線画像装置の表面上で変わるステップと、 前記Ｘ線画像装置に衝突するＸ線の強度を示す電気信号を生成するステップと、 前記電気信号を用いて前記物体の境界を捜し出すステップと、 前記物体の境界を用いてＸ線画像を生成するステップと、を含む、物体のＸ線画像を生成する方法。
【請求項２】
前記物体が胸部である、請求項１記載の方法。
【請求項３】
前記物体の境界を捜し出すステップが、前記胸部の皮膚の境界を捜し出すことを含む、請求項２記載の方法。
【請求項４】
前記電気信号がディジタル値の配列を含む、請求項１記載の方法。
【請求項５】
前記物体の境界を捜し出すステップが、前記配列の線に沿って各々の複数の値の勾配を計算することを含む、請求項４記載の方法。
【請求項６】
前記物体の境界を捜し出すステップが、前記配列の線に沿ってディジタル値の最大勾配を判定することを含む、請求項４記載の方法。
【請求項７】
前記物体の境界を捜し出すステップが、前記配列の各々の複数の線に沿ってディジタル値の最大勾配を判定することを含む、請求項６記載の方法。
【請求項８】
前記物体の境界を捜し出すステップが、前記配列の各線に沿ってディジタル値の最大勾配を判定することを含む、請求項７記載の方法。
【請求項９】
前記物体の境界を捜し出すステップが、前記物体の境界に対応する配列のディジタル値のセットを選択することを含む、請求項４記載の方法。
【請求項１０】
前記Ｘ線画像を生成するステップが、 前記セットのディジタル値を平均化するステップと、 前記平均化ディジタル値を用いて境界値を計算するステップと、を含む、請求項９記載の方法。
【請求項１１】
前記Ｘ線画像を生成するステップが、 前記境界値以上の配列のディジタル値を特定するステップと、 前記特定されたディジタル値を境界値に設定するステップと、を含む、請求項１０記載の方法。
【請求項１２】
前記Ｘ線画像を生成するステップが、 前記境界値未満の配列のディジタル値を特定するステップと、 前記特定されたディジタル値の平均及び標準偏差を計算するステップと、 前記特定されたディジタル値にＳ字状変換を行うステップであって、このＳ字状変換関数が、前記計算された平均から得られる中央値を有する範囲と前記計算された標準偏差から得られる傾斜とを有するステップと、を含む、請求項１０記載の方法。
【請求項１３】
前記Ｓ字状変換の範囲の中央値が、前記計算された平均の０．９５倍である、請求項１２記載の方法。
【請求項１４】
前記Ｓ字状変換の傾斜が、前記計算された標準偏差の１．２倍である、請求項１２記載の方法。
【請求項１５】
前記Ｘ線画像を生成するステップが、 前記Ｓ字状変換関数の出力に整定数を掛けて中間出力セットを形成するステップと、 前記中間出力セットの各要素から前記中間出力セットの最小値を引くステップと、 出力装置での出力に適した出力値の範囲に前記中間出力セットを測定して出力セットを得るステップと、を含む、請求項１２記載の方法。
【請求項１６】
前記Ｘ線画像を生成するステップが、ディジタル値の範囲を圧縮することを含む、請求項４記載の方法。
【請求項１７】
前記ディジタル値の範囲を圧縮するステップが、対数関数を前記ディジタル値に適用することを含む、請求項１６記載の方法。
【請求項１８】
前記Ｘ線画像を生成するステップが、ディジタル値をフィルタリングすることを含む、請求項４記載の方法。
【請求項１９】
前記ディジタル値をフィルタリングするステップが、低域空間周波数を前記ディジタル値からフィルタリングすることを含む、請求項１８記載の方法。
【請求項２０】
前記ディジタル値をフィルタリングするステップが、緩いマスキングを前記ディジタル値に施すことを含む、請求項１８記載の方法。
【請求項２１】
前記緩いマスキングを前記ディジタル値に施すステップが、第１のマスクを前記ディジタル値から生成することを含む、請求項２０記載の方法。
【請求項２２】
前記第１のマスクが、カーナルサイズ約１００×１００を用いて前記ディジタル値から生成される、請求項２１記載の方法。
【請求項２３】
前記緩いマスキングを前記ディジタル値に施すステップが、第２のマスクを前記第１のマスクから生成することを含む、請求項２１記載の方法。
【請求項２４】
前記緩いマスキングを前記ディジタル値に施すステップが、前記第２のマスクに定数を掛けるステップと、この積を前記ディジタル値の配列から引いてマスクされたデータセットを生成するステップとを含む、請求項２１記載の方法。
【請求項２５】
前記定数が約０．５５である、請求項２４記載の方法。
【請求項２６】
前記定数が０．５５である、請求項２５記載の方法。
【請求項２７】
前記第２のマスクが、カーナルサイズ約７５×７５を用いて前記第１のマスクから生成される、請求項２３記載の方法。
【請求項２８】
Ｘ線源とＸ線画像装置との間に物体を配置する手段と、 前記Ｘ線源を作動させてＸ線を前記Ｘ線画像装置まで前記物体の中及び周囲に向け、これにより、前記Ｘ線画像装置に衝突するＸ線の強度が前記Ｘ線画像装置の表面上で変わる手段と、 前記Ｘ線画像装置に衝突するＸ線の強度を示す電気信号を生成する手段と、 前記電気信号を用いて前記物体の境界を捜し出す手段と、 前記物体の境界を用いてＸ線画像を生成する手段と、を含む、物体のＸ線画像を生成する装置。
【請求項２９】
Ｘ線源とＸ線画像装置との間に物体を配置するステップと、 前記Ｘ線源を作動させてＸ線を前記Ｘ線画像装置まで前記物体の中に向け、これにより、前記Ｘ線画像装置に衝突するＸ線の強度が前記Ｘ線画像装置の表面上で変わるステップと、 電気信号のセットを生成するステップであって、各電気信号が、前記Ｘ線画像装置の領域に対応付けられ、このような領域に衝突するＸ線の強度を示す大きさを有するステップと、 前記電気信号のセット内で最高度数の大きさを特定するステップと、 前記最高度数の大きさを用いてＸ線画像を生成するステップと、を含む、物体のＸ線画像を生成する方法。
【請求項３０】
前記最高度数の大きさを特定するステップが、前記セット内の電気信号
の大きさのヒストグラムを生成することを含む、請求項２９記載の方法。
【請求項３１】
前記Ｘ線画像を生成するステップが、前記最高度数の大きさからしきい値の大きさの値を判定することを含む、請求項２９記載の方法。
【請求項３２】
前記Ｘ線画像を生成するステップが、前記セット内の電気信号の大きさから前記しきい値の大きさの値を引くことを含む、請求項３１記載の方法。
【請求項３３】
前記しきい値の大きさの値を判定するステップが、前記最高度数の大きさに対応する度数、ＭａｘＦｒｅｑｕｅｎｃｙを判定することを含む、請求項３１記載の方法。
【請求項３４】
前記しきい値の大きさの値を判定するステップが、ＭａｘＦｒｅｑｕｅｎｃｙを第１の定数で割ってＬｏｗＦｒｅｑｕｅｎｃｙを判定することを含む、請求項３３記載の方法。
【請求項３５】
前記第１の定数が約１０００である、請求項３４記載の方法。
【請求項３６】
前記しきい値の大きさの値を判定するステップが、ＬｏｗＦｒｅｑｕｅｎｃｙに対応する大きさから第２の定数を引くことを含む、請求項３４記載の方法。
【請求項３７】
前記第２の定数が約３００である、請求項３６記載の方法。
【請求項３８】
Ｘ線源とＸ線画像装置との間に物体を配置する手段と、 前記Ｘ線源を作動させてＸ線を前記Ｘ線画像装置まで前記物体の中に向け、これにより、前記Ｘ線画像装置に衝突するＸ線の強度が前記Ｘ線画像装置の表面上で変わる手段と、 電気信号のセットを生成する手段であって、各電気信号が、前記Ｘ線画像装置の領域に対応付けられ、このような領域に衝突するＸ線の強度を示す大きさを有する手段と、 前記電気信号のセット内で最高度数の大きさを特定する手段と、 前記最高度数の大きさを用いてＸ線画像を生成する手段と、を含む、物体のＸ線画像を生成する装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、胸部のＸ線画像（「乳房Ｘ線像」）などのディジタルＸ線画像の向上に関する。
【背景技術】
【０００２】
乳房Ｘ線像は、人体の胸部内の異常を早期発見し易くするために胸部の内部構造の様相を明示する。Ｘ線乳房撮影においては、各種組織、例えば、健康な腺及び／又は脂肪組織、及び／又は腫瘍及び／又は微小石灰化などの異常組織から成る胸部の中にＸ線を向ける。
【０００３】
異なる組織のタイプでは、密度が異なる。高密度の組織は、低密度の組織よりもＸ線を減衰させる。密度が相対的に近い組織間の区別を行える性能を有すると、胸部内の異常又は疑わしい病巣が初期状態でまだ小さい場合でも、そのような病巣の検出を容易にすることができる。
【０００４】
その後、減衰されたＸ線は捕捉されて乳房Ｘ線像を生成する。あるシステムでは、１枚の感光フィルムが胸部の後ろに置かれる。露光後、前記フィルムを現像して診断のために見る。他のシステムでは、電気画像装置を用いてＸ線を電気信号又は画像データに変換し、これを例えば診断用レーザプリンタで印刷することができる。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
本発明の一面においては、画像処理すべき物体を、Ｘ線源とＸ線画像装置との間に配置する。このＸ線源を作動させてＸ線を前記物体の中及び周囲に向ける。前記Ｘ線画像装置に衝突する前記Ｘ線の強度を示す電気信号を用いて、前記物体の境界を捜し出し、この境界を用いてＸ線画像を生成する。
【０００６】
多くの利点の中でも、本発明のこの側面により、例えば、境界外の領域を示す電気信号とは別々に、境界内の領域を示す電気信号を全部又は一部処理することができる。これにより、処理速度を速くできるだけでなく、境界外の領域を示す電気信号が境界内の領域を示す電気信号の処理に影響を及ぼす範囲を削除しない場合にその範囲を縮小することもできる。
【課題を解決するための手段】
【０００７】
本発明の上記側面の実施形態は、次のような特徴を含むこともできる。前記電気信号は、ディジタル値の配列を含むことができる。物体の境界（例えば、胸部の皮膚境界）は、前記配列の各線に沿ってディジタル値の最大勾配を計算することにより捜し出すことができる。境界値は、物体の境界に対応する配列のディジタル値のセットを平均化することにより求められる。前記境界値以上の配列の前記ディジタル値は境界値に等しく設定され、境界値未満のディジタル値の平均及び標準偏差が計算される。次に、前記計算された平均により判定された中央値（median）を有する範囲（例えば、計算された平均の０．９５倍に等しい中央値を有する範囲）と、前記計算された標準偏差により判定された傾斜（例えば、計算された標準偏差の１．２倍に等しい傾斜）とを有するＳ字状変換は、前記ディジタル値の後者のセットで行われる。前記Ｓ字状変換関数の出力に整定数を掛けて中間出力セット（intermediary output set）を形成し、この中間出力セットの最小値を前記セットの各要素から引く。出力セットは、出力装置での出力に適した出力値の範囲に前記中間出力セットを測定（scaling）することにより得られる。
【０００８】
また、配列のディジタル値の範囲は、圧縮（例えば、対数関数を使用）及び／又はフィルタリングすることもできる。フィルタリングを用いて、例えば緩いマスキングを施すことにより、低域空間周波数を前記ディジタル値からフィルタリングすることができる。第１のマスク（例えば、カーナルサイズ（kernal size）１００×１００を使用）は、前記ディジタル値から生成し、第２のマスク（例えば、カーナルサイズ７５×７５を使用）は、前記第１のマスクから生成することができる。次に、前記第２のマスクに定数（例えば、０．５５）を掛け前記ディジタル値から引いて、マスクされたデータのセットを得ることができる。
【０００９】
本発明のもう一つの側面においては、画像処理すべき物体を、Ｘ線源とＸ線画像装置との間に配置する。このＸ線源を作動させてＸ線を前記物体の中に向け、前記Ｘ線画像装置の領域に対応付けられ、かつこのような領域に衝突する前記Ｘ線の強度を示す大きさを有する電気信号を生成する。前記電気信号のセット内の最高度数の大きさを用いて、Ｘ線画像を生成する。
【００１０】
多くの利点の中でも、本発明のこの側面により、例えば、電気信号の範囲を使用可能な解像度まで広げて、画像細部の品質向上及び／又は出力装置での出力に適した画像データの作成が可能となる。
【００１１】
本発明の上記側面の実施形態は、次のような特徴を含むこともできる。前記最高度数の大きさは、前記セットの電気信号の大きさのヒストグラムを生成することにより特定できる。しきい値の大きさの値は、前記最高度数の大きさに対応する度数、ＭａｘＦｒｅｑｕｅｎｃｙを判定し、このＭａｘＦｒｅｑｕｅｎｃｙを約１０００で割ってＬｏｗＦｒｅｑｕｅｎｃｙを判定し、このＬｏｗＦｒｅｑｕｅｎｃｙに対応する大きさから約３００を引くことにより計算できる。次いで、前記しきい値の大きさの値を、前記セットの電気信号の大きさから引いてもよい。
【００１２】
他の利益と特徴は、以下の好ましい実施形態と請求の範囲とから明白になるであろう。
【発明の効果】
【００１３】
本発明は、胸部のＸ線画像（「乳房Ｘ線像」）などのディジタルＸ線画像の向上を達成することができる。
【発明を実施するための最良の形態】
【００１４】
ディジタル乳房撮影の画像取得段階において、Ｘ線管１１０を含むＸ線源１００は、圧迫仕切板１５０ａ−ｂにより圧迫された人体の胸部１３０の中及び周囲にＸ線１２０を向ける。前記Ｘ線源１００は、前記胸部１３０の中及び周囲に下方へＸ線１２０を向けるよう胸部１３０の上方に位置決めすることができる。１つ以上の検出器を含むディジタル画像装置１４０は、胸部１３０の下に横設してある。前記画像装置の出力を用いて、胸部のＸ線画像を生成する。
【００１５】
図２に示すように、前記ディジタル画像装置１４０は、シンチレータ１６０、光ファイバーテーパ１７０の層、及び二次元配列１８０に配置された電荷結合デバイス(ＣＣＤ)の層を含む。各ＣＣＤは光電性素子のグリッドを含み、その各々が、得られた乳房Ｘ線画像２５０（図３）の画素２６５に対応する。（説明し易いように、胸部画像は図３に示していない。）典型的なＣＣＤ上の各素子のサイズは、約４０×４０ミクロンである。
【００１６】
画像処理中、Ｘ線管１１０を作動させて胸部１３０の中及び周囲にＸ線１２０を向ける。Ｘ線の一部は、胸部により散乱される。
【００１７】
胸部１３０の中及び周囲を通過した後、Ｘ線１２０は画像装置１４０のシンチレータ１６０に当たり、そのシンチレータが可視光を発する。この可視光は、光ファイバーテーパ１７０の中を通過してＣＣＤ配列１８０に達する。ＣＣＤ配列の各光電性素子は、前記素子の領域に当たった前記光の強度に対応する電圧を出力する。しかし、すべてのＣＣＤ素子が、同一のＸ線強度に等しく感応するとは限らない。このことは「白色フィールド補正」によりある程度説明されているが、素子の一部は一般的に、他の素子よりも低い露光で飽和する。
【００１８】
画像装置１４０は、ＣＣＤ１８０からの出力電圧データをアナログ／ディジタル変換器１９０に伝送し、この変換器は前記データをディジタル形式に変換する。その後、これらのディジタル画像データはオプションとして、プリプロセッサ２００の中を通過させ、例えば、ずれ除去、継ぎ合わせ(例えば、個々のＣＣＤからのデータを単一の二次元配列画像データに結合する)、欠陥マッピング、及び白色フィールディングを行うことができる。得られた画像２５０の全解像度は、例えば４８００×６４００画素でよい。その後、前記データは、後述のルーチンを実現するソフトウェアを実行するコンピュータワークステーション２１０に送られる。このコンピュータ２１０は通常、キーボード２２０又はその他入力装置に装着されている。コンピュータが画像データの処理を完了すると、この処理データは、モニタ２３０上に表示したり、プリンタ２４０で印刷したり、又はハードディスクやその他固定媒体に格納したりすることができる。
【００１９】
コンピュータ２１０がプリプロセッサ２００から受信する画像データは通常、約１４ビットの範囲を有する。各画像データは、乳房Ｘ線画像２５０の単一画素２６５に対応する。コンピュータ２１０上のアルゴリズムは、画像の左上隅に原点があって其処から右にｘの値が増加し下にｙの値が増加するデカルト座標系をデータセットに割り当てる。しかし、本発明はある特定の画素座標系に左右されることはない。更に、ここに記載したアルゴリズムは、胸部が正面の方を指して水平に向いていると仮定しているが、アルゴリズムを修正して胸部を異なる向きに動かすこともできる。
【００２０】
まず、アルゴリズムは、「胸部皮膚境界」を捜し出して、其処から「胸部境界値」を計算する。この胸部皮膚境界は、胸部の真の幾何学的境界に必ずしも一致しないが、その幾何学的境界の輪郭の代わりとなる。前記胸部境界値は、胸部皮膚境界上に存在する画像データ点の平均強度である。
【００２１】
図４は、胸部皮膚境界を捜し出して胸部境界値を計算するコンピュータ実現ルーチンのフロー図である。
【００２２】
このルーチンは、画像データの行に沿って最高強度の勾配を有する画像データの各行の画素を捜し出し、このような全画素の軌跡として胸部皮膚境界を特定する。まず、前記ルーチンは、行インデックス局所変数（local row index variable）ｙをゼロに初期設定する(ステップ３１０)。また、行ｙ上で生じる最大勾配を示す局所変数ｍａｘ＿ｇｒａｄをゼロに初期設定し、ｍａｘ＿ｇｒａｄに対応する画素のｘ位置を示す局所変数ｍａｘ＿ｇｒａｄ＿ｘもゼロに初期設定する。
【００２３】
行ｙに沿った最大勾配は、ステップ３４０〜３８０で見つけられる。具体的には、ステップ３４０において、座標ｘ，ｙでの画素の勾配は、ソーベル（sobel）演算を用いて｜Ｇx｜＋｜Ｇy｜を近似することにより計算される(ステップ３４０)。ｘ，ｙでの勾配がｍａｘ＿ｇｒａｄ（前述の最大勾配）に格納された値よりも大きい場合、ｍａｘ＿ｇｒａｄをそのｘ，ｙでの勾配に設定し、ｍａｘ＿ｇｒａｄ＿ｘをｘに設定する(ステップ３５０、３６０)。このような処理が、行ｙのあらゆる画素に対して繰り返される(ステップ３７０、３８０)。
【００２４】
最大勾配（ｍａｘ＿ｇｒａｄ）及びそれに対応する行ｙのｘ座標（ｍａｘ＿ｇｒａｄ＿ｘ）を見つけた後、ｍａｘ＿ｇｒａｄの値をｍａｘ＿ｇｒａｄ＿ａｒｒａｙに格納し、その画素の座標（ｍａｘ＿ｇｒａｄ＿ｘ，ｙ）をＢｒｅａｓｔ＿Ｂｏｕｎｄａｒｙ配列に追加し、この配列は胸部皮膚境界に対応する画素のセットを格納する(ステップ３９０)。変数ｙを１だけ増加させてから(ステップ４００)、前記処理をすべての行に対して繰り返す(ステップ４１０)。
【００２５】
上記方法で画像の各行を処理すると、胸部境界値がｍａｘ＿ｇｒａｄ＿ａｒｒａｙに格納された値の平均として算出される(ステップ４２０)。
【００２６】
胸部皮膚境界及び胸部境界値を算出した後、アルゴリズムは画像データの範囲を圧縮する。この圧縮を行うルーチンは、図５に示してある。まず、胸部マップを作成して、画像データの各データを胸部データ又は非胸部データに分類する(ステップ５００)。具体的には、画像データの大きさが胸部境界値よりも大きい場合、その画像データを非胸部データに分類する。残りの画像データを胸部データに分類する。胸部マップを用いて、各非胸部データの大きさを胸部境界値に設定する(ステップ５００)。これにより、胸部の外側に均質な場所が効果的に作成される。
【００２７】
次に、図６に示すように、各々の大きさにおける胸部データ点の数を示す胸部データのヒストグラムを計算する(ステップ５０５)。（説明し易いように、図６では少数のデータ点のみが示してある。）
【００２８】
その後、しきい値度数未満のデータを廃棄して、残りのデータをステップ５１０〜５３０に移行させる。具体的には、ステップ５１０において、ステップ５０５で計算されたヒストグラムの最高点を捜し出して、胸部データ点の最大数により示される大きさ（即ち、最高度数の大きさ）(１)と、その大きさを示す胸部データ点の数（即ち、最高度数）(２)との両方を特定する。後者の値を、変数ＭａｘＦｒｅｑｕｅｎｃｙに割り当てる(図６)。次に、変数ＬｏｗＦｒｅｑｕｅｎｃｙをＭａｘＦｒｅｑｕｅｎｃｙ／１０００として計算する(ステップ５１５)。再度ヒストグラムを用いて、ＬｏｗＦｒｅｑｕｅｎｃｙに対応する大きさを決定する。この大きさを変数ＬｏｗＦｒｅｑｕｅｎｃｙＭａｇに割り当てる(ステップ５２０)。次いで、ＴｈｒｅｓｈｏｌｄＭａｇを、３００とＬｏｗＦｒｅｑｕｅｎｃｙＭａｇのうち小さい方をＬｏｗＦｒｅｑｕｅｎｃｙＭａｇから引いて計算する(ステップ５２５)。ＴｈｒｅｓｈｏｌｄＭａｇ未満の大きさのデータ点を廃棄する。その後、残りの各データ点からＴｈｒｅｓｈｏｌｄＭａｇを引くことによりデータ範囲をシフトさせる（即ち、前記データ点はＴｈｒｅｓｈｏｌｄＭａｇより大なる大きさを有する）(ステップ５３０)。そして、浮動小数点出力付き対数関数を、前記シフトされたデータに実行する(ステップ５４０)。
【００２９】
次に、緩いマスキングを用いて、画像データから低域空間周波数をフィルタリングして画像細部の品質向上を図る。一般に、緩いマスキングルーチンは、ソース画像のマスク画像を作成してから、そのマスク画像の関数をソース画像から引いて結果画像を作成する。座標ｘ，ｙにおけるマスク画像の値は、座標ｘ，ｙを中心とした領域（「カーナル(kernal)）内のソース画像データの平均値に等しい。
【００３０】
図７は、画像データに緩いマスキングを施すコンピュータ実現ルーチンのフロー図である。まず、カーナルサイズ１００×１００画素を用いて、マスクＭａｓｋ１を画像データから作成する(ステップ６００)。後述の図８は、マスク値を計算するコンピュータ実現ルーチンのフロー図である。次に、カーナルサイズ７５×７５画素を用いて、マスクＭａｓｋ２をＭａｓｋ１から作成する(ステップ６１０)。次いで、Ｍａｓｋ２に０．５５を掛けたものを画像データから引くことにより結果画像Ｒｅｓｕｌｔを得る(ステップ６２０)。Ｒｅｓｕｌｔは、以後の全ルーチンで画像データとして使用される(ステップ６４０)。
【００３１】
図８に示すように、マスク値を計算するルーチンは、所望のマスク値のｘ座標とｙ座標、及び幅Ｗと高さＨにより表されるカーナルサイズを入力として採る。
【００３２】
前記ルーチンは、変数Ｓｕｍ、ｉ及びｊをゼロに初期設定して開始する(ステップ８００、８１０、８２０)。変数Ｓｕｍを用いて、前記カーナル内のデータ点の和を把握しておく。変数ｉ及びｊを用いて、カーナル内のデータ点のｘ座標及びｙ座標をそれぞれ把握しておく。
【００３３】
まず、前記ルーチンは、カーナルの開始点（即ち、カーナルの左上側）の座標を計算する(ステップ８３０)。この点が画像境界の外側にある場合、その代わりにルーチンは、画像境界内の点の鏡像を使用する(ステップ８４０〜８９０)。
【００３４】
次に、点（実際の点か又はその実際の点が画像境界の外側にある場合はその鏡像の何れか）の強度をステップ９００で求める。この強度を、強度の実行中の和に追加する(ステップ９１０)。前記処理を、カーナル内の各点に対して繰り返す(ステップ９３０、９４０、９５０)。最後に、強度の平均値を指定された座標に対するマスク値として戻す(ステップ９６０)。
【００３５】
図７の緩いマスキングルーチンにより低域空間周波数を画像からフィルタリングした後、結果画像データ（Ｒｅｓｕｌｔ）を、０から２５５におよぶ２５６グレーレベルの範囲を有するデータに変換して、普通のレーザプリンタ及びその他の出力装置での出力に適するようにする。図９は、この処理のフロー図である。
【００３６】
まず、ステップ５００で特定した非胸部データを引いた結果画像データ（Ｒｅｓｕｌｔ）の平均及び標準偏差を計算する(ステップ１０００)。
【００３７】
次に、計算された平均及び標準偏差をパラメータとして有するＳ字状変換関数を、画像データに実行する(ステップ１０１０)。具体的には、前記Ｓ字状変換関数の範囲は前記平均の０．９５倍に等しい中央値であり、前記関数の傾斜は前記標準偏差の１．２倍に等しく設定されている。また、得られたデータを、得られた各値に１６３８３を掛けることにより［０，１６３８３］の整数範囲に正規化する(ステップ１０２０)。次いで、正規化されたデータのヒストグラムを計算して(ステップ１０３０)、各データ値から最小データ値を引くことによりデータをゼロの方へシフトさせる(ステップ１０４０)。最後に、シフトされたデータを[０，２５５]の整数範囲に正規化する(ステップ１０５０)。
【００３８】
その後、得られた画像データを、レーザプリンタやその他の適当な出力装置に送信したり、ハードディスクドライブやその他の記憶媒体のファイルに格納したりすることができる。
【００３９】
他の実施形態では、画素配列から成る単一のフラットパネル型検出器を用いた画像装置の使用を含む。緩いマスキングの代わりに、又はこのマスキングに追加して、フーリエ変換を用いて、低域空間周波数を画像データからフィルタリングすることもできる。
【００４０】
本発明は、ソフトウェアにより実現する内容で記載されているけれども、ソフトウェア、ハードウェア、ファームウェア、又はこれら３つの組合せで本発明を実現することもできる。
【００４１】
本発明は、上述された実施形態に限定されない。もっと正確に言えば、本発明の範囲は請求項により規定される。
【図面の簡単な説明】
【００４２】
【図１】乳房撮影装置の略図である。
【図２】ディジタル乳房Ｘ線像用の画像データを収集及び処理する装置の略図である。
【図３】二次元配列の画素から成る乳房Ｘ線像空間の略図である。
【図４】画像データの胸部皮膚線を捜し出すコンピュータ実現方法のフロー図である。
【図５】画像データの範囲を圧縮するコンピュータ実現方法のフロー図である。
【図６】図５のルーチンにより生成されたヒストグラムの略図である。
【図７】画像データに緩いマスキングを施すコンピュータ実現方法のフロー図である。
【図８】図７の緩いマスキングルーチンにおけるマスク値を計算するコンピュータ実現方法のフロー図である。
【図９】画像データのグレースケール変換を行うコンピュータ実現方法のフロー図である。
【符号の説明】
【００４３】
１００    Ｘ線源
１２０    Ｘ線
１３０    物体
１４０    Ｘ線画像装置
【図１】
2007313360.tif 000003
【図２】
2007313360.tif 000004
【図３】
2007313360.tif 000005
【図４】
2007313360.tif 000006
【図５】
2007313360.tif 000007
【図６】
2007313360.tif 000008
【図７】
2007313360.tif 000009
【図８】
2007313360.tif 000010
【図９】
2007313360.tif 000011
【手続補正書】
【提出日】平成19年9月7日(2007.9.7)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
Ｘ線源とＸ線画像装置との間に物体を配置する手段と、 
Ｘ線が前記物体を通過して前記Ｘ線画像装置まで向かうように、前記Ｘ線源を作動させる手段であって、これにより、前記Ｘ線画像装置に衝突するＸ線の強度が前記Ｘ線画像装置の表面上で変わる手段と、
電気信号のセットを生成する手段であって、各電気信号が、前記Ｘ線画像装置の領域に対応付けられ、このような領域に衝突するＸ線の強度を示す大きさを有する手段と、
前記電気信号のセット内で最高頻度を有する信号の大きさを特定する手段と、 
前記最高頻度を有する信号の大きさに少なくとも一部基づいて定められる範囲に画素値を有するＸ線画像を生成する手段と、を含む、物体のＸ線画像を生成する装置。
【手続補正２】
【補正対象書類名】図面
【補正対象項目名】図１
【補正方法】変更
【補正の内容】
【図１】
2007313360.tif 000012