Nothing Special   »   [go: up one dir, main page]

JPH10108073A - 骨部画像処理方法および装置 - Google Patents

骨部画像処理方法および装置

Info

Publication number
JPH10108073A
JPH10108073A JP8256943A JP25694396A JPH10108073A JP H10108073 A JPH10108073 A JP H10108073A JP 8256943 A JP8256943 A JP 8256943A JP 25694396 A JP25694396 A JP 25694396A JP H10108073 A JPH10108073 A JP H10108073A
Authority
JP
Japan
Prior art keywords
processing
image
bone
image data
trabecular
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.)
Withdrawn
Application number
JP8256943A
Other languages
English (en)
Inventor
Kazuo Shimura
一男 志村
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujifilm Holdings Corp
Original Assignee
Fuji Photo Film Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fuji Photo Film Co Ltd filed Critical Fuji Photo Film Co Ltd
Priority to JP8256943A priority Critical patent/JPH10108073A/ja
Priority to US08/937,602 priority patent/US5910972A/en
Publication of JPH10108073A publication Critical patent/JPH10108073A/ja
Withdrawn legal-status Critical Current

Links

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)
  • Transforming Light Signals Into Electric Signals (AREA)
  • Image Analysis (AREA)

Abstract

(57)【要約】 【課題】 骨部画像処理装置において、骨梁の経時変化
をより観察しやすい状態で示す。 【解決手段】 軟部組織Ａと骨部組織Ｂとを含む被検体
を被写体とする、撮影時期が相異なる複数の画像をそれ
ぞれ表す複数の画像データを記憶する記憶手段10、記憶
手段10から比較しようとする対象の２つの画像について
の画像データ（第１Ｘ線画像データＳ１、第２Ｘ線画像
データＳ２）を読み出す読出手段20、読み出された２つ
の各画像データＳ１，Ｓ２に対してそれぞれ骨部組織Ｂ
における骨梁Ｃの強調処理を施す骨梁強調手段30、およ
び強調処理された画像データＳ１′，Ｓ２′を、両画像
間で位置対応がとれるように位置合わせ処理を施し、両
画像データＳ１′，Ｓ２′間で、相対応する位置の画素
についてのデータ毎に減算処理を行なって、上記骨梁の
経時変化を示す差分データΔＳ12を求める画像処理手段
40を備える。

Description

【発明の詳細な説明】
【０００１】
【発明の属する技術分野】本発明は骨部画像処理方法お
よび装置に関し、詳細には、同一被写体の時間的に相異
なる画像データに基づいて、この被写体の骨梁の変化状
態を観察するための骨梁画像処理方法および装置に関す
るものである。
【０００２】
【従来の技術】放射線画像やＣＴ画像等の医用画像によ
る診断、研究においては、腫瘍等の病変部の病状進行状
態や治癒状態を調べる等のために、同一被写体を撮影時
点を変えて複数回撮影し、そうして得られた複数の画像
から被写体の経時変化を判断することが多い。
【０００３】従来、このようにして被写体の経時変化を
判断するために医師等の観察者は、Ｘ線写真フィルムに
記録された撮影時点の相異なる複数の画像を比較観察す
るようにしていた。なお、例えば特開平１−107739号公
報等に示されているように、医用画像をＣＲＴ表示装置
等の画像表示手段に表示することもあるが、その場合も
事情は基本的に同じである。つまりこの場合、画像表示
手段には、撮影時点の相異なる複数の画像が同時表示あ
るいは切換え表示され、医師等の観察者がそれら複数の
表示画像を比較観察して、被写体の経時変化を判断する
ようにしていた。
【０００４】しかし、上記のように複数の画像を比較す
る場合、医師等の観察者が不慣れであったりすると、経
時変化がどの程度なのか良く分からないこともある。
【０００５】そこで本願出願人は、撮影時点の相異なる
複数の画像のうち比較対象である２つの撮影時期におい
て撮影された各画像を、これら画像を示す画像データ上
で被写体の位置合わせを行ない、対応する画素の画像デ
ータ間で減算処理することにより、両画像データの差分
を求め、この差分データを可視像として表示することに
より、両撮影時点間における経時的変化を画像表示手段
に表示し、医師等の観察者が被写体の経時変化を容易か
つ明確に観察できる技術を提案している（特願平７−61
040 号、同７−61041 号）。
【０００６】
【発明が解決しようとする課題】しかしながら上記本願
出願人による上記技術によっても、軟部組織の経時変化
やアーティファクトを完全に取り除けるものではないた
め、特に骨粗鬆症の診断に有効な骨部組織の骨梁の経時
的変化を、軟部組織からの影響等を極力排除して、より
観察しやすくする技術が望まれている。
【０００７】本発明は上記事情に鑑みなされたものであ
って、被写体、特に骨梁の経時変化をより観察しやすい
状態で示すことができる骨部画像処理方法および装置を
提供することを目的とするものである。
【０００８】
【課題を解決するための手段】本発明の骨部画像処理方
法および装置は、同一被写体についての撮影時点が相異
なる２つの放射線画像を表す画像データをそれぞれ骨梁
強調することにより、両撮影時点間における軟部組織の
厚さの変化等による骨部組織（骨梁を含む）の画像に与
えるビームハードニングの差および両撮影時点における
被写体の撮影姿勢の差異により生じる軟部組織の位置ず
れに起因する軟部組織のアーティファクトを相対的に低
下せしめ、そのうえでエネルギーサブトラクション処理
を施すことにより、軟部組織による影響が極めて少ない
骨梁の経時的変化分の画像データを得ることを特徴とす
るものである。
【０００９】すなわち、本発明の骨部画像処理方法は、
撮影時点が相異なる複数の、軟部組織と骨部組織とを含
む被検体を被写体とする画像を表す画像データ、が記憶
せしめられた画像記憶手段から、比較対象の２つの画像
についての画像データをそれぞれ読み出し、該読み出さ
れた各画像データに対してそれぞれ前記骨部組織におけ
る骨梁の強調処理を施し、これらの強調処理された画像
データを、該２つの画像の間で位置対応が取れるように
位置合わせ処理し、この位置合わせ処理後の、前記強調
処理された２つの画像データ間で、相対応する位置の画
素についてのデータ毎に減算処理を行なって、前記骨梁
の経時変化を示す差分データを求めることを特徴とする
ものである。
【００１０】ここで上記比較対象の２つの画像について
の画像データが外部から入力される場合は、記憶手段か
らこれらの画像データを読み出す操作に代えて、入力さ
れた２つの画像データを認識するという操作としてもよ
い。
【００１１】また上記位置合わせ処理としては、対象の
画像を撮影する際に所定の位置合わせ用マーカーを写し
込んでいる場合はそのマーカーを一致させる処理や、文
献 Medical Imaging Technology （メディカル・イメー
ジング・テクノロジー）Vol.11 No.3 July 1993 pp.373
〜374 に示されている２次元の非線形画像変形による位
置合わせ処理等種々の公知の方法を適用することができ
る。
【００１２】そしてこのようにして得られた上記差分デ
ータは記録手段や表示手段等に出力すればよい。
【００１３】なお、前記各撮影時点ごとの画像データと
して、該各撮影時点においてエネルギーサブトラクショ
ン処理のために撮影されたエネルギー分布が互いに異な
る２つの画像をそれぞれ表す２つの画像データを適用す
ることもできる。この場合、前記骨梁の強調処理を施す
のに先だって、前記各撮影時点ごとの前記２つの画像デ
ータに基づいて前記骨部組織を抽出するエネルギーサブ
トラクション処理を施して骨部画像データを得、各撮影
時点ごとに得られた骨部画像データに対して前記骨梁の
強調処理、前記位置合わせ処理および前記減算処理を施
せばよい。ただし、画像記憶手段に記憶せしめられてい
る画像データが、すでにエネルギーサブトラクション処
理を施されて得られた骨部画像データである場合は、再
度のエネルギーサブトラクション処理を行なう必要はな
い。
【００１４】また前記骨梁の強調処理後、前記位置合わ
せ処理前に、位置合わせ処理の対象となる各画像データ
に対して平滑化処理を施し、この平滑化処理がなされた
画像データに対して位置合わせ処理および減算処理を施
すのが望ましい。位置合わせ処理前に平滑化処理を行な
うことにより、両撮影時点における被写体のわずかな撮
影姿勢の差異に起因する骨梁の位置ずれのアーティファ
クトを低減することができるからである。
【００１５】なお、骨梁の強調処理を施す対象の画像デ
ータが高濃度高信号レベルの画像データである場合に
は、上記骨梁の強調処理として、下記式（１）で示され
るモーフォロジー演算に基づくスケルトン処理を適用す
るのが好ましく、
【００１６】
【数１】
【００１７】一方、処理の対象の画像データが高輝度高
信号レベルの画像データである場合には、上記骨梁の強
調処理として、下記式（２）で示されるモーフォロジー
演算に基づくスケルトン処理を適用するのが好ましい。
【００１８】
【数２】
【００１９】ここでモーフォロジー演算（以下、モーフ
ォロジー処理ともいう）およびスケルトン処理について
説明する。
【００２０】モーフォロジー処理とは、原画像のうち異
常な陰影等の特定の画像部分だけを選択的に抽出する、
モーフォロジー（Morphology；モフォロジーまたはモル
フォロジーとも称する）のアルゴリズムに基づく処理で
あり、特に乳癌における特徴的形態である微小石灰化像
を検出するのに有効な手法として研究されているが、対
象画像としてはこのようなマンモグラムにおける微小石
灰化像に限るものではない。
【００２１】そしてこのモーフォロジー処理は、抽出し
ようとする画像部分の大きさ、形状に対応した構造要素
Ｂを用いた処理を行なうものであり、複雑なバックグラ
ウンド情報に影響されにくい、抽出した画像が歪まな
い、などの特徴がある。
【００２２】すなわち、この手法は一般の微分処理に比
べて、石灰化像のサイズ・形状・濃度分布などの幾何学
的情報をよりよく保って検出することができる。
【００２３】以下、このモーフォロジー処理の概要を、
マンモグラムにおける微小石灰化像の検出に適用した例
について説明する。
【００２４】（モーフォロジーの基本演算）モーフォロ
ジー処理は一般的にはＮ次元空間における集合論として
展開されるが、直感的な理解のために２次元の濃淡画像
を対象として説明する。
【００２５】濃淡画像を座標（ｘ，ｙ）の点が濃度値ｆ
（ｘ，ｙ）に相当する高さをもつ空間とみなす。ここ
で、濃度値ｆ（ｘ，ｙ）は、濃度が低い（ＣＲＴに表示
した場合には輝度が高い）程大きな画像信号値となる高
輝度高信号レベルの信号とする。
【００２６】まず、簡単のため上記２次元の濃淡画像の
断面に相当する１次元の関数ｆ（ｘ）を考える。モーフ
ォロジー演算に用いる構造要素ｇは次式（３）に示すよ
うに、原点について対称な対称関数
【００２７】
【数３】
【００２８】であり、定義域内で値が０で、その定義域
Ｇが下記式（４）であるとする。
【００２９】
【数４】
【００３０】このとき、モーフォロジー演算の基本形は
式（５）〜（８）に示すように、非常に簡単な演算とな
る。
【００３１】
【数５】
【００３２】すなわち、ダイレーション（dilation）処
理は、注目画素を中心とした、±ｍ（構造要素Ｂに応じ
て決定される値であって、図６中のマスクサイズに相
当）の幅の範囲内の最大値を探索する処理であり（同図
（Ａ）参照）、一方、エロージョン（erosion ）処理
は、注目画素を中心とした、±ｍの幅の範囲内の最小値
を探索する処理である（同図（Ｂ）参照）。また、オー
プニング（opening ）処理はエロージョン処理後にダイ
レーション処理を行なう処理、すなわち最小値の探索の
後に最大値を探索する処理であり、クロージング（clos
ing ）処理は、ダイレーション処理後にエロージョン処
理を行なう処理、すなわち最大値の探索の後に最小値を
探索する処理に相当する。
【００３３】つまりオープニング処理は、低輝度側から
濃度曲線ｆ（ｘ）を滑らかにし、マスクサイズ２ｍより
空間的に狭い範囲で変動する凸状の濃度変動部分（周囲
部分よりも輝度が高い部分）を取り除くことに相当する
（同図（Ｃ）参照）。
【００３４】一方、クロージング処理は、高輝度側から
濃度曲線ｆ（ｘ）を滑らかにし、マスクサイズ２ｍより
空間的に狭い範囲で変動する凹状の濃度変動部分（周囲
部分よりも輝度が低い部分）を取り除くことに相当する
（同図（Ｄ）参照）。
【００３５】なお、構造要素ｇが原点に対して対称では
ない場合の、式（５）に示すダイレーション演算をミン
コフスキー（Minkowski ）和、式（６）に示すエロージ
ョン演算をミンコフスキー差という。
【００３６】ここで、濃度の高いもの程大きな値となる
高濃度高信号レベルの信号の場合においては、濃度値ｆ
（ｘ）の画像信号値が高輝度高信号レベルの場合に対し
て大小関係が逆転するため、高濃度高信号レベルの信号
に対するダイレーション処理と高輝度高信号レベルに対
するエロージョン処理（同図（Ｂ））とは一致し、高濃
度高信号レベルの信号に対するエロージョン処理と高輝
度高信号レベルに対するダイレーション処理（同図
（Ａ））とは一致し、高濃度高信号レベルの信号に対す
るオープニング処理と高輝度高信号レベルに対するクロ
ージング処理（同図（Ｄ））とは一致し、高濃度高信号
レベルの信号に対するクロージング処理と高輝度高信号
レベルに対するオープニング処理（同図（Ｃ））とは一
致する。
【００３７】なお、本項では高輝度高信号レベルの画像
信号（輝度値）の場合について説明する。
【００３８】（石灰化陰影検出への応用）石灰化陰影の
検出には、原画像から平滑化した画像を引き去る差分法
が考えられる。単純な平滑化法では石灰化陰影と細長い
形状の非石灰化陰影（乳腺、血管および乳腺支持組織
等）との識別が困難であるため、東京農工大の小畑ら
は、多重構造要素を用いたオープニング処理に基づく下
記式（９）で表されるモーフォロジーフィルターを提案
している（「多重構造要素を用いたモルフォロジーフィ
ルタによる微小石灰化像の抽出」電子情報通信学会論文
誌 D-II Vol.J75-D-II No.7 P1170 〜1176 1992年７
月、「モルフォロジーの基礎とそのマンモグラム処理へ
の応用」MEDICAL IMAGING TECHNOLOGY Vol.12 No.1 Jan
uary 1994 等）。
【００３９】
【数６】
【００４０】ここでＢi （ｉ＝１，２，…，ｎ）は、直
線状の大きさがｍ画素でｎ個（例えば図７に示すもので
は、９画素４方向でありｍ＝９，ｎ＝４）の構造要素
（これらを全体として以下、ｍ画素ｎ方向の多重構造要
素という）である。構造要素Ｂi を検出対象である石灰
化陰影よりも大きく設定すれば、上記オープニング演算
による処理で、構造要素Ｂi よりも細かな信号変化部分
（空間的に狭い範囲で信号が変動する画像部分）であっ
て周囲よりも輝度値の大きい凸状の部分である石灰化陰
影は取り除かれる。一方、細長い形状の乳腺の陰影等の
非石灰化陰影はその長さが構造要素Ｂi よりも長く、そ
の傾きが４つの構造要素Ｂi のいずれかに一致すればオ
ープニング処理（式（９）の第２項の演算）をしてもそ
のまま残る。したがってオープニング処理によって得ら
れた平滑化画像（石灰化陰影のみが取り除かれた画像）
を原画像ｆから引き去ることで、小さな石灰化陰影のみ
が含まれる画像が得られる。これが式（９）の考え方で
ある。
【００４１】なお、前述したように、高濃度高信号レベ
ルの信号の場合においては、石灰化陰影は周囲の画像部
分よりも濃度値が低くなり、石灰化陰影は周囲部分に対
して濃度値の小さい凹状の信号変化部分となるため、オ
ープニング処理に代えてクロージング処理を適用し、式
（９）に代えて式（10）を適用する。
【００４２】
【数７】
【００４３】なお、モーフォロジー演算の一例である式
（10）のクロージング処理を具体的に説明する。
【００４４】すなわち、高濃度高信号レベルの画像信号
である濃度値Ｓorg についてのモーフォロジー演算によ
れば、例えば図８（１）の実線に示すような濃度値Ｓor
g の分布を有する画像データに対して、同図（２）に示
すような直線状の３画素の構造要素Ｂで、最大値処理
（ダイレーション処理）を行うことにより、ある注目画
素の濃度値Ｓi は、その注目画素を中心として互いに隣
接する３画素（構造要素Ｂにより決定される）の中の最
大値Ｓi+1 を採用したＳi ′に変換される。この演算を
全画素について行うことにより、濃度値Ｓorg ′の分布
を有する同図（１）の破線で示す最大値信号に変換され
る。
【００４５】次に、この最大値処理により得られた最大
値信号に対してさらに構造要素Ｂによる最小値処理（エ
ロージョン処理）を考えると、同図（１）の破線で示さ
れた注目画素の最大値信号Ｓi ′は、その注目画素を中
心として互いに隣接する３画素の中の最小値Ｓi-1 ′を
採用したＳi ″（＝Ｓi ）に変換される。この演算を全
画素について行うことにより、最大値処理後の最小値信
号Ｓorg ″の分布は同図（１）の一点鎖線で示すものと
される。この一点鎖線で示された画像信号は、もとの実
線のオリジナルの画像データに対して、構造要素Ｂより
も空間的に狭い範囲で信号が変動する画像部分が消え、
構造要素Ｂよりも空間的に広い範囲で変動する信号値の
変化部分である画像部分や変動のない画像部分はもとの
形状のまま残っていることを示している。すなわち、以
上の処理（クロージング処理）は、画像濃度の分布を高
濃度側から平滑化する処理として作用する。
【００４６】このようにクロージング処理で得られた値
（Ｓorg に対して最大値処理を行なった後にさらに最小
値処理を行なった値）を原画像信号Ｓorg から差し引く
ことにより得られた値Ｓmor は、上記クロージング処理
で消された空間的に狭い範囲で変動する信号値の変化部
分である画像部分を表す。
【００４７】ここで、本来、画像信号は２次元の要素で
ある位置（ｘ，ｙ）と、３次元目の要素である信号値ｆ
（ｘ，ｙ）を有するが、上記説明においては、理解の容
易化のために、この２次元上に展開された画像の所定の
断面に現れた、１次元状の画像信号分布曲線について説
明した。
【００４８】したがって実際には、以上の説明を２次元
画像に適用する必要があり、多重構造要素を用いるのも
２次元画像に対応させるためである。
【００４９】次にモーフォロジー演算に基づくスケルト
ン処理について説明する。
【００５０】スケルトン処理とは一般に図形の骨格（sk
eleton）を抽出する処理であり、骨格とは図形に内接す
る円盤の中心の集合としてとらえることができる。すな
わち例えば図９（Ａ）〜（Ｅ）に示す各図形（中太の実
線で表記）の骨格はそれぞれ太い実線で示すものとな
る。
【００５１】以下、このスケルトン処理を、上述したモ
ーフォロジー演算により行なう場合について説明する。
この場合、スケルトン処理は下記式（１）または（２）
で表すことができる。
【００５２】
【数１】
【００５３】
【数２】
【００５４】ここで式（１）と式（２）とは前述した、
画像を高濃度高信号レベルの画像信号として表すか、ま
たは高輝度高信号レベルの画像信号として表すかの差異
によるものであり、高濃度高信号レベルの画像信号で表
された画像から低濃度（高輝度）の画像部分の骨格を抽
出する場合は式（１）を適用し、一方、高輝度高信号レ
ベルの画像信号で表された画像から低輝度（高濃度）の
画像部分の骨格を抽出する場合は式（２）を適用するも
のであり、作用自体の実質的な差はない。
【００５５】例えばネガフイルム（高濃度高信号レベ
ル）上においては骨部は他の画像部分に比べてその濃度
は低いものとなり、骨梁の存在する部分は濃度が低く、
存在しない部分は濃度が高くなる。したがって、周囲よ
りも濃度の低い部分となる骨梁を対象としてスケルトン
処理を行なうことに該当するため式（１）を適用すれば
よい。
【００５６】ここで式（１）における構造要素Ｂを半径
ｒの円とし、図10に示す図形に対してスケルトン処理を
施した状態を示す。図10に示す図形は、その輪郭よりも
外側領域が濃度の高い部分であり、内側が濃度の低い部
分である。
【００５７】この図形に対してまず図示上段では構造要
素Ｂによるエロージョン処理が施される。λ＝０（構造
要素Ｂによる０回のエロージョン処理）では図形は何ら
変化がない。
【００５８】λ＝１（構造要素Ｂによる１回のエロージ
ョン処理）では図形が構造要素Ｂの半径ｒ分だけ内側に
埋め込められる。
【００５９】λ＝２（構造要素Ｂによる２回のエロージ
ョン処理）では図形の円から突出した部分が完全に消失
する。
【００６０】同様の操作を繰り返すことにより、λ＝Ｎ
−１（構造要素ＢによるＮ−１回のエロージョン処理）
で図形は半径ｒ以下の円のみとなる。
【００６１】一方、図示中段は、この構造要素Ｂによる
各回（λ＝０，１，２，…，Ｎ−１，Ｎ）のエロージョ
ン処理を施した画像に対して、さらにそれぞれ構造要素
Ｂによるオープニング処理を施した図形である。
【００６２】この図示上段の図形から図示中段の図形
を、処理の回数を対応させて差し引いた図形が図示下段
である。
【００６３】λ＝１において、元の図形の円から突出し
た部分の骨格要素が抽出され、λ＝Ｎ−１において、元
の図形の円の骨格要素が抽出されていることが分かる。
【００６４】このように元の図形に対してエロージョン
処理を施し、さらにオープニング処理を施し、処理回数
を対応させて差し引き、この結果の和集合を求めたもの
が式（１）の意味するところである。
【００６５】式（２）では濃度の高低が式（１）とは逆
転した図形から骨格要素を抽出するのに有効であり、元
の図形に対してダイレーション処理を施し、さらにクロ
ージング処理を施し、処理回数を対応させて差し引き、
この結果の和集合を求めることにより骨格要素を抽出す
ることを意味する。
【００６６】このように骨梁の強調処理として、式
（１）または（２）で示されるモーフォロジー演算に基
づくスケルトン処理を適用すれば、複雑なバックグラン
ド情報からの影響を抑制して骨梁（若しくはその骨梁の
骨格）だけを効果的に強調することができる。
【００６７】なお上記式（１）、（２）中の和集合とし
て、ｎ1 の値が比較的大きいもののみの和集合、すなわ
ち例えばｎ1 ＝０，１の場合を除いたλ＝２，３，４，
５等の和集合を採用し、この場合に得られる骨格要素を
表示することにより、骨梁の変化をより見易くすること
ができ、好ましい。これは、λが０や１のときは特に、
非常に空間周波数の高いノイズ成分も抽出されるため、
和集合からこれらのノイズ成分を除外することによっ
て、より読影性能の高い画像を得ることができるからで
ある。以下の発明においても同様である。
【００６８】本発明の骨部画像処理装置は、本発明の骨
部画像処理方法を実施するための装置であって、撮影時
点が相異なる複数の、軟部組織と骨部組織とを含む被検
体を被写体とする画像を表す画像データを記憶する記憶
手段と、この記憶手段から比較対象の２つの画像につい
ての画像データを読み出す読出手段と、これらの読み出
された各画像データに対してそれぞれ前記骨部組織にお
ける骨梁の強調処理を施す骨梁強調手段と、これらの強
調処理された画像データを、該２つの画像の間で位置対
応が取れるように位置合わせ処理を施し、この位置合わ
せ処理後の、前記強調処理された２つの画像データ間
で、相対応する位置の画素についてのデータ毎に減算処
理を行なって、前記骨梁の経時変化を示す差分データを
求める画像処理手段とを備えたことを特徴とするもので
ある。
【００６９】ここで、前記各撮影時点ごとの画像データ
に、該各撮影時点においてエネルギーサブトラクション
処理のために撮影されたエネルギー分布が互いに異なる
２つの画像をそれぞれ表す２つの画像データを適用する
こともできる。この場合、前記各撮影時点ごとの前記２
つの画像データに基づいて前記骨部組織を抽出するエネ
ルギーサブトラクション処理を施すエネルギーサブトラ
クション手段をさらに備え、前記骨梁強調手段が、該エ
ネルギーサブトラクション処理により得られた骨部画像
データに対して、前記骨梁の強調処理を施すものとすれ
ばよいなお画像記憶手段に記憶せしめられている画像デ
ータが、すでにエネルギーサブトラクション処理を施さ
れて得られた骨部画像データである場合は、再度のエネ
ルギーサブトラクション処理を行なう必要はないため、
上記エネルギーサブトラクション手段をさらに備える必
要はない。
【００７０】また撮影姿勢の差異に起因する骨梁の位置
ずれのアーティファクトを低減するため、上記位置合わ
せ処理の対象となる各画像データに対して、平滑化処理
を施す平滑化処理手段をさらに備えた構成とするのが好
ましい。
【００７１】なお、骨梁強調手段による骨梁の強調処理
を施す対象の画像データが高濃度高信号レベルの画像デ
ータである場合には、上記骨梁の強調処理として、式
（１）で示されるモーフォロジー演算に基づくスケルト
ン処理を適用するのが好ましく、一方、処理の対象の画
像データが高輝度高信号レベルの画像データである場合
には、上記骨梁の強調処理として、式（２）で示される
モーフォロジー演算に基づくスケルトン処理を適用する
のが好ましい。
【００７２】
【発明の効果】本発明の骨部画像処理方法および装置
は、同一被写体についての撮影時点が相異なる２つの放
射線画像を表す画像データをそれぞれ骨梁強調すること
により、両撮影時点間における軟部組織の厚さの変化等
による骨部組織（骨梁を含む）の画像に与えるビームハ
ードニングの差および両撮影時点における被写体の撮影
姿勢の差異により生じる軟部組織の位置ずれに起因する
軟部組織のアーティファクトを相対的に低下せしめ、そ
のうえでエネルギーサブトラクション処理を施すことに
より、軟部組織による影響が極めて少ない骨梁の経時変
化分の画像データを得ることができる。したがってこの
得られた画像データに基づいて、フイルム等の記録媒体
あるいはＣＲＴ等の表示手段に可視像として出力すれ
ば、上記軟部組織による影響が極めて少ない骨梁の経時
変化分を観察することができる。
【００７３】さらに上記放射線画像を表す画像データと
して、高エネルギー画像と低エネルギー画像とをエネル
ギーサブトラクション処理により骨部組織を強調または
抽出した骨部画像データを用いることもでき、上記放射
線画像を表す画像データをそのまま使用する場合より
も、軟部組織による影響を抑制することができる。
【００７４】また、骨梁の強調処理後、前記位置合わせ
処理前に、位置合わせ処理の対象となる各画像データに
対して平滑化処理を施し、この平滑化処理がなされた画
像データに対して位置合わせ処理および減算処理を施す
のが望ましい。位置合わせ処理前に平滑化処理を行なう
ことにより、両撮影時点における被写体のわずかな撮影
姿勢の差異に起因する骨梁の位置ずれのアーティファク
トを低減することができる。
【００７５】なお、骨梁の強調処理を施す対象の画像デ
ータが高濃度高信号レベルの画像データである場合にお
いて、上記骨梁の強調処理として、式（１）または式
（２）で示されるモーフォロジー演算に基づくスケルト
ン処理を適用すれば、複雑なバックグランド情報からの
影響を抑制して骨梁（若しくはその骨梁の骨格）だけを
効果的に強調することができる。
【００７６】
【発明の実施の形態】以下、本発明の骨部画像処理装置
の具体的な実施の形態について図面を用いて説明する。
【００７７】図１は本発明の骨部画像処理装置の第１の
実施形態の構成を示すブロック図、図２は図１に示す骨
部画像処理装置の作用の流れを示す図、である。
【００７８】図示の骨部画像処理装置100 は、軟部組織
Ａと骨部組織Ｂとを含む被検体を被写体とする、撮影時
期が相異なる複数の画像をそれぞれ表す複数の画像デー
タを記憶する記憶手段10と、この記憶手段10から経時的
な骨梁の変化を比較しようとする対象の２つの画像（第
１Ｘ線画像Ｐ１と第２Ｘ線画像Ｐ２）についての画像デ
ータ（第１Ｘ線画像データＳ１、第２Ｘ線画像データＳ
２）を読み出す読出手段20と、これらの読み出された２
つの各画像データＳ１，Ｓ２に対してそれぞれ骨部組織
Ｂにおける骨梁Ｃの強調処理を施す骨梁強調手段30と、
これらの強調処理された画像データＳ１′，Ｓ２′を、
これら２つの画像の間で位置対応がとれるように位置合
わせ処理を施し、この位置合わせ処理後の上記強調処理
された２つの画像データＳ１′，Ｓ２′間で、相対応す
る位置の画素についてのデータ毎に減算処理を行なっ
て、上記骨梁の経時変化を示す差分データを求める画像
処理手段40とを備えた構成である。
【００７９】ここで骨梁強調手段30としては、一例とし
て下記式（１）で示されるモーフォロジー演算に基づく
スケルトン処理を施すモーフォロジー演算に基づくスケ
ルトン処理手段を適用するものとする。
【００８０】
【数１】
【００８１】図３は本実施形態の骨部画像処理装置の作
用によりその経時変化が抽出される様子を可視化した模
式図である。
【００８２】以下、本実施形態の骨部画像処理装置100
の作用について説明する。
【００８３】図１に示した骨部画像処理装置100 の記憶
手段10には、同一被写体について撮影時期を変えて撮影
した複数のＸ線画像Ｐｉをそれぞれ表す画像データＳｉ
が記憶されている。
【００８４】ここで読出手段20が、例えば医師または放
射線技師等の観察読影者から入力された、比較しようと
する撮影時期の相異なる２つの画像データＳ１（第１Ｘ
線画像データ）とＳ２（第２Ｘ線画像データ）とを読み
出す。
【００８５】この読み出された第１Ｘ線画像データＳ１
と第２Ｘ線画像データＳ２は、骨梁強調手段30に入力さ
れる。骨梁強調手段30は入力された第１Ｘ線画像データ
Ｓ１と第２Ｘ線画像データＳ２に対して、前記式（１）
で示されるモーフォロジー演算に基づくスケルトン処理
を施し、各画像データＳ１，Ｓ２から骨部組織Ｂの骨梁
Ｃを強調処理する。なお、このモーフォロジー演算に用
いられる構造要素は、骨梁の骨格を抽出するのに適した
大きさのものを予め設定しておく。
【００８６】骨梁Ｃは骨の内部に縦横に延びる海綿質で
あり、放射線画像においては縦横に延びる海綿質が幾重
にも重なって写し出されるが、骨梁強調手段30による前
述したモーフォロジー演算に基づくスケルトン処理によ
って、図３に示すように海綿質の骨格だけが幾重にも重
なった画像として抽出される。なおこの強調処理された
第１Ｘ線画像データＳ１′が表す第１Ｘ線画像Ｐ１′お
よび第２Ｘ線画像データＳ２′が表す第２Ｘ線画像Ｐ
２′はいずれも、軟部組織Ａが殆ど消去された画像とな
る。
【００８７】次に、この骨梁Ｃが強調処理された第１Ｘ
線画像データＳ１′と第２Ｘ線画像データＳ２′は画像
処理手段40に入力される。
【００８８】画像処理手段40には、この骨梁強調手段30
を経由した第１Ｘ線画像データＳ１′および第２Ｘ線画
像データＳ２′、並びに骨梁強調手段30を経由しない第
１Ｘ線画像データＳ１および第２Ｘ線画像データＳ２が
入力される。
【００８９】ここで、記憶手段10に記憶された各画像デ
ータＳｉがそれぞれ表す画像Ｐｉには、画像処理手段40
により２つの画像を位置合わせするための指標となるマ
ーカーＤがそれぞれ写し込まれているが、骨梁強調手段
30によって強調処理された第１Ｘ線画像データＳ１′、
第２Ｘ線画像データＳ２′においては、このマーカーＤ
が消去されてしまう。そこで画像処理手段40は、骨梁強
調手段30を経由しないで入力された第１Ｘ線画像データ
Ｓ１および第２Ｘ線画像データＳ２におけるマーカーＤ
の位置に基づいて位置合わせのための情報（位置合わせ
情報という）を算出し、この位置合わせ情報に基づいて
第１Ｘ線画像データＳ１′と第２Ｘ線画像データＳ２′
とを位置合わせする。
【００９０】さらにこの位置合わせがなされた第１Ｘ線
画像データＳ１′と第２Ｘ線画像データＳ２′とについ
て、画像処理手段40が相対応する位置の画素についての
データ毎に減算処理を行なって、上記骨梁の経時変化を
示す差分データΔＳ12を求める。
【００９１】求められた差分データΔＳ12は本実施形態
の骨部画像処理装置100 から画像表示装置200 等の画像
出力装置に出力されて、差分データΔＳ12が表す可視
像、すなわち骨梁の経時変化が表示される。
【００９２】この表示された可視像は、原画像である第
１Ｘ線画像データＳ１、第２Ｘ線画像データＳ２をそれ
ぞれ骨梁強調することにより、両撮影時点間における軟
部組織の厚さの変化等による骨梁の画像に与えるビーム
ハードニングの差および両撮影時点における被写体の撮
影姿勢の差異により生じる軟部組織の位置ずれに起因す
る軟部組織のアーティファクトを相対的に低下させるこ
とができ、このアーティファクトが相対的に低下した画
像に基づいてサブトラクション処理を施すことにより、
得られた画像は軟部組織による影響が極めて少ないもの
となる。
【００９３】そしてこの得られた画像は骨梁の経時変化
を示す画像データΔＳ12であるから、この画像データΔ
Ｓ12に基づいて、フイルム等の記録媒体あるいはＣＲＴ
等の表示手段に可視像として出力することにより、上記
軟部組織による影響が極めて少ない骨梁の経時変化が表
示されるため、診断が容易な観察しやすい画像を得るこ
とができる。
【００９４】本実施形態の骨部画像処理装置は、骨梁強
調手段としてモーフォロジー演算に基づくスケルトン処
理手段を用いたが、本発明の骨部画像処理装置はこの態
様に限るものではなく、骨梁を強調し得る他の種々の手
段を用いることができる。
【００９５】また、画像処理手段における位置合わせ処
理は、本実施形態のようにマーカーを用いて行なう方法
に限らず、前述した２次元の非線形画像変形による位置
合わせ処理（文献 Medical Imaging Technology Vol.11
No.3 July 1993 pp.373〜374 ）等種々の公知の方法を
適用することができる。
【００９６】図４は本発明の骨部画像処理装置の第２の
実施形態の構成を示すブロック図である。図示の骨部画
像処理装置は、図１に示した実施形態の装置に対して、
骨梁強調手段30と画像処理手段40との間に平滑化手段50
を介在せしめた点と、本装置が処理対象とする画像デー
タとして、原画像データ（第１Ｘ線画像データＳ１およ
び第２Ｘ線画像データＳ２）に代えて、エネルギーサブ
トラクション（図示では「エネサブ」と略記）処理（特
開昭59-83486号等参照）された画像データ（第１エネル
ギーサブトラクション画像データおよび第２エネルギー
サブトラクション画像データ）を用いる点が相違する以
外は図１に示した骨部画像処理装置と同一である。
【００９７】この第２の実施形態の骨部画像処理装置で
は、エネルギーサブトラクション画像データを処理対象
とすることにより、予め軟部組織をある程度消去してお
くことができ、軟部組織によるアーティファクト等の影
響をさらに低減した骨梁Ｃの経時的変化を示す画像デー
タを得ることができる。
【００９８】また、骨梁強調手段30により強調処理して
得られた各撮影時点における骨梁Ｃを表す画像データに
対して、平滑化手段50による平滑化処理を施すことによ
り、各撮影時点における被写体の撮影姿勢の差異により
生じる骨部組織の位置ずれに起因する骨梁のアーティフ
ァクトをも低減させることができ、このようなアーティ
ファクトが低減された骨梁強調画像間で画像処理装置40
による減算処理を行なうことにより、第１の実施形態の
骨部画像処理装置よりもさらに観察読影性能に優れた骨
梁の経時変化を示す画像データを得ることができる。
【００９９】なお、上記平滑化手段による平滑化処理と
しては、非鮮鋭マスク処理（マトリクスによる空間フィ
ルタリング）等の種々の公知の方法を適用することがで
き、また第１の実施形態においても同様に適用すること
ができる。
【０１００】第２の実施形態（図４）においては、記憶
手段10に記憶されている画像データが、既にエネルギー
サブトラクション処理がなされた結果であるエネサブ画
像データである場合の構成を示したが、エネルギーサブ
トラクション処理前の高圧画像および低圧画像が記憶手
段10に記憶されているとする構成の形態（第３の実施形
態）を採ることもできる（図５参照）。
【０１０１】この第３の実施形態の場合、骨梁強調処
理、平滑化処理、位置合わせ処理、減算処理はエネルギ
ーサブトラクション画像データを対象とするが、記憶さ
れている画像データはエネルギーサブトラクション処理
前の画像データであるため、第２の実施形態の装置に加
えてさらにエネルギーサブトラクション処理手段（図示
では「エネサブ手段」と略記）60を備えており、エネル
ギーサブトラクション処理手段60が、各撮影時点ごとの
高圧画像および低圧画像に基づいて各撮影時点ごとのエ
ネルギーサブトラクション画像データを求め、求められ
たエネルギーサブトラクション画像データに対して上記
骨梁強調処理、平滑化処理、位置合わせ処理、減算処理
がなされる。その他の構成、作用、効果は第２の実施形
態の装置と同様である。
【図面の簡単な説明】
【図１】本発明の骨部画像処理装置の第１の実施形態の
構成を示すブロック図
【図２】図１に示す骨部画像処理装置の作用の流れを示
す図
【図３】図１に示した実施形態の骨部画像処理装置の作
用によりその経時変化が抽出される様子を可視化した模
式図
【図４】本発明の骨部画像処理装置の第２の実施形態の
構成を示すブロック図
【図５】本発明の骨部画像処理装置の第３の実施形態の
構成を示すブロック図
【図６】モーフォロジー処理の基本的な作用を示す図、
（Ａ）ダイレーション（dilation）処理、（Ｂ）エロー
ジョン（erosion ）処理、（Ｃ）オープニング（openin
g ）処理、（Ｄ）クロージング（closing ）処理
【図７】モーフォロジー処理に用いる構造要素Ｂi の一
例を示す図
【図８】クロージング処理を具体的に説明する図
【図９】種々の図形とその骨格（skeleton：スケルト
ン）を示す図
【図１０】スケルトン処理を具体的に説明する図
【符号の説明】
10 記憶手段 20 読出手段 30 骨梁強調手段 40 画像処理手段 50 平滑化手段 60 エネルギーサブトラクション処理手段 100 骨部画像処理手段 200 表示装置

Claims (14)

    【特許請求の範囲】
  1. 【請求項１】 撮影時点が相異なる複数の、軟部組織と
    骨部組織とを含む被検体を被写体とする画像を表す画像
    データが記憶せしめられた画像記憶手段から、比較対象
    の２つの画像についての画像データをそれぞれ読み出
    し、 該読み出された各画像データに対してそれぞれ前記骨部
    組織における骨梁の強調処理を施し、 これらの強調処理された画像データを、該２つの画像の
    間で位置対応が取れるように位置合わせ処理し、 この位置合わせ処理後の、前記強調処理された２つの画
    像データ間で、相対応する位置の画素についてのデータ
    毎に減算処理を行なって、前記骨梁の経時変化を示す差
    分データを求めることを特徴とする骨部画像処理方法。
  2. 【請求項２】 前記各撮影時点ごとの画像データが、該
    各撮影時点においてエネルギーサブトラクション処理の
    ために撮影されたエネルギー分布が互いに異なる２つの
    画像をそれぞれ表す２つの画像データからなり、 前記骨梁の強調処理を施すのに先だって、前記各撮影時
    点ごとの前記２つの画像データに基づいて前記骨部組織
    を抽出するエネルギーサブトラクション処理を施して骨
    部画像データを得、 該各撮影時点ごとに得られた骨部画像データに対して前
    記骨梁の強調処理、前記位置合わせ処理および前記減算
    処理を施すことを特徴とする請求項１記載の骨部画像処
    理方法。
  3. 【請求項３】 前記骨梁の強調処理後、前記位置合わせ
    処理前に、該位置合わせ処理の対象となる各画像データ
    に対して平滑化処理を施し、 該平滑化処理がなされた画像データに対して前記位置合
    わせ処理および前記減算処理を施すことを特徴とする請
    求項１または２記載の骨部画像処理方法。
  4. 【請求項４】 前記骨梁の強調処理が、下記式（１）で
    示されるモーフォロジー演算に基づくスケルトン処理で
    あり、該モーフォロジー演算を適用する画像が高濃度高
    信号レベルの画像信号で表されるものであることを特徴
    とする請求項１から３のうちいずれか１項に記載の骨部
    画像処理方法。 【数１】
  5. 【請求項５】 前記式（１）におけるｎ1 の値が２以上
    であることを特徴とする請求項４記載の骨部画像処理方
    法。
  6. 【請求項６】 前記骨梁の強調処理が、下記式（２）で
    示されるモーフォロジー演算に基づくスケルトン処理で
    あり、該モーフォロジー演算を適用する画像が高輝度高
    信号レベルの画像信号で表されるものであることを特徴
    とする請求項１から３のうちいずれか１項に記載の骨部
    画像処理方法。 【数２】
  7. 【請求項７】 前記式（２）におけるｎ1 の値が２以上
    であることを特徴とする請求項６記載の骨部画像処理方
    法。
  8. 【請求項８】 撮影時点が相異なる複数の、軟部組織と
    骨部組織とを含む被検体を被写体とする画像を表す画像
    データを記憶する記憶手段と、 この記憶手段から比較対象の２つの画像についての画像
    データを読み出す読出手段と、 これらの読み出された各画像データに対してそれぞれ前
    記骨部組織における骨梁の強調処理を施す骨梁強調手段
    と、 これらの強調処理された画像データを、該２つの画像の
    間で位置対応が取れるように位置合わせ処理を施し、こ
    の位置合わせ処理後の、前記強調処理された２つの画像
    データ間で、相対応する位置の画素についてのデータ毎
    に減算処理を行なって、前記骨梁の経時変化を示す差分
    データを求める画像処理手段とを備えたことを特徴とす
    る骨部画像処理装置。
  9. 【請求項９】 前記各撮影時点ごとの画像データが、該
    各撮影時点においてエネルギーサブトラクション処理の
    ために撮影されたエネルギー分布が互いに異なる２つの
    画像をそれぞれ表す２つの画像データからなり、 前記各撮影時点ごとの前記２つの画像データに基づいて
    前記骨部組織を抽出するエネルギーサブトラクション処
    理を施すエネルギーサブトラクション手段をさらに備
    え、 前記骨梁強調手段が、該エネルギーサブトラクション処
    理により得られた骨部画像データに対して、前記骨梁の
    強調処理を施すものであることを特徴とする請求項８記
    載の骨部画像処理装置。
  10. 【請求項１０】 前記位置合わせ処理の対象となる各画
    像データに対して平滑化処理を施す平滑化処理手段をさ
    らに備えたことを特徴とする請求項８または９記載の骨
    部画像処理装置。
  11. 【請求項１１】 前記骨梁強調手段による骨梁の強調処
    理が、下記式（１）で示されるモーフォロジー演算に基
    づくスケルトン処理であり、該モーフォロジー演算を適
    用する画像が高濃度高信号レベルの画像信号で表される
    ものであることを特徴とする請求項８から１０のうちい
    ずれか１項に記載の骨部画像処理装置。 【数１】
  12. 【請求項１２】 前記式（１）におけるｎ1 の値が２以
    上であることを特徴とする請求項１１記載の骨部画像処
    理装置。
  13. 【請求項１３】 前記骨梁強調手段による骨梁の強調処
    理が、下記式（２）で示されるモーフォロジー演算に基
    づくスケルトン処理であり、該モーフォロジー演算を適
    用する画像が高輝度高信号レベルの画像信号で表される
    ものであることを特徴とする請求項８から１０のうちい
    ずれか１項に記載の骨部画像処理装置。 【数２】
  14. 【請求項１４】 前記式（２）におけるｎ1 の値が２以
    上であることを特徴とする請求項１３記載の骨部画像処
    理装置。
JP8256943A 1996-09-25 1996-09-27 骨部画像処理方法および装置 Withdrawn JPH10108073A (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP8256943A JPH10108073A (ja) 1996-09-27 1996-09-27 骨部画像処理方法および装置
US08/937,602 US5910972A (en) 1996-09-25 1997-09-25 Bone image processing method and apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP8256943A JPH10108073A (ja) 1996-09-27 1996-09-27 骨部画像処理方法および装置

Publications (1)

Publication Number Publication Date
JPH10108073A true JPH10108073A (ja) 1998-04-24

Family

ID=17299520

Family Applications (1)

Application Number Title Priority Date Filing Date
JP8256943A Withdrawn JPH10108073A (ja) 1996-09-25 1996-09-27 骨部画像処理方法および装置

Country Status (1)

Country Link
JP (1) JPH10108073A (ja)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002197440A (ja) * 2000-12-27 2002-07-12 Fuji Photo Film Co Ltd 放射線画像データ処理方法および装置
JP2003515376A (ja) * 1999-12-03 2003-05-07 コミツサリア タ レネルジー アトミーク 二重エネルギーコーンビームｘ線放射を使用した骨濃度測定システムの使用方法
JP2006346459A (ja) * 2005-06-14 2006-12-28 General Electric Co <Ge> 超音波空間複合イメージングのリアルタイム運動補正のための方法及び装置
JP2008200103A (ja) * 2007-02-16 2008-09-04 Konica Minolta Medical & Graphic Inc 骨疾患評価システム
JP2008200114A (ja) * 2007-02-16 2008-09-04 Konica Minolta Medical & Graphic Inc 骨梁評価システム
WO2011114733A1 (ja) * 2010-03-17 2011-09-22 富士フイルム株式会社 医用画像変換装置および方法並びにプログラム
JP2011193997A (ja) * 2010-03-19 2011-10-06 Fujifilm Corp 医用画像変換装置および方法並びにプログラム
JP2011193887A (ja) * 2010-03-17 2011-10-06 Fujifilm Corp 医用画像変換装置および方法並びにプログラム
WO2013005533A1 (ja) * 2011-07-06 2013-01-10 富士フイルム株式会社 内視鏡システム、内視鏡システムのプロセッサ装置、及び画像表示方法
JP2014518125A (ja) * 2011-07-06 2014-07-28 コーニンクレッカ フィリップス エヌ ヴェ フォローアップ画像の取得計画及び／又は後処理
JP2015080555A (ja) * 2013-10-22 2015-04-27 コニカミノルタ株式会社 医用画像処理装置及びプログラム
JP2017119056A (ja) * 2015-12-29 2017-07-06 ラトックシステムエンジニアリング株式会社 骨代謝解析支援プログラム及び方法

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003515376A (ja) * 1999-12-03 2003-05-07 コミツサリア タ レネルジー アトミーク 二重エネルギーコーンビームｘ線放射を使用した骨濃度測定システムの使用方法
JP2002197440A (ja) * 2000-12-27 2002-07-12 Fuji Photo Film Co Ltd 放射線画像データ処理方法および装置
JP2006346459A (ja) * 2005-06-14 2006-12-28 General Electric Co <Ge> 超音波空間複合イメージングのリアルタイム運動補正のための方法及び装置
JP2008200103A (ja) * 2007-02-16 2008-09-04 Konica Minolta Medical & Graphic Inc 骨疾患評価システム
JP2008200114A (ja) * 2007-02-16 2008-09-04 Konica Minolta Medical & Graphic Inc 骨梁評価システム
US9295442B2 (en) 2010-03-17 2016-03-29 FUJIFILM Corporation ko Medical image conversion apparatus, method and program
WO2011114733A1 (ja) * 2010-03-17 2011-09-22 富士フイルム株式会社 医用画像変換装置および方法並びにプログラム
JP2011193887A (ja) * 2010-03-17 2011-10-06 Fujifilm Corp 医用画像変換装置および方法並びにプログラム
JP2011193997A (ja) * 2010-03-19 2011-10-06 Fujifilm Corp 医用画像変換装置および方法並びにプログラム
JP2013013656A (ja) * 2011-07-06 2013-01-24 Fujifilm Corp 内視鏡システム、内視鏡システムのプロセッサ装置、及び画像表示方法
JP2014518125A (ja) * 2011-07-06 2014-07-28 コーニンクレッカ フィリップス エヌ ヴェ フォローアップ画像の取得計画及び／又は後処理
WO2013005533A1 (ja) * 2011-07-06 2013-01-10 富士フイルム株式会社 内視鏡システム、内視鏡システムのプロセッサ装置、及び画像表示方法
US9675311B2 (en) 2011-07-06 2017-06-13 Koninklijke Philips N.V. Follow up image acquisition planning and/or post processing
JP2015080555A (ja) * 2013-10-22 2015-04-27 コニカミノルタ株式会社 医用画像処理装置及びプログラム
JP2017119056A (ja) * 2015-12-29 2017-07-06 ラトックシステムエンジニアリング株式会社 骨代謝解析支援プログラム及び方法

Similar Documents

Publication Publication Date Title
JP4130661B2 (ja) 時間的に連続する胸部画像間の経時変化を検出する装置
JP5753791B2 (ja) ノイズ除去された所定の解像度の医療画像を提供する方法、所定の解像度のノイズ除去された所定の解像度の医療画像を提供するシステム
JPH1031745A (ja) 自動画像解析方法
US5910972A (en) Bone image processing method and apparatus
JPH11112877A (ja) 骨計測方法および装置
JPH08294479A (ja) 計算機支援画像診断装置
JP3590216B2 (ja) 異常陰影候補の検出方法および装置
JP3669455B2 (ja) 画像処理方法および画像処理装置
JP2002099896A (ja) 異常陰影候補検出方法および装置
CN110610498A (zh) 乳腺钼靶图像处理方法、系统、存储介质及设备
JPH10108073A (ja) 骨部画像処理方法および装置
JP3931792B2 (ja) 時系列処理画像の表示装置及び表示方法
JP3731760B2 (ja) 画像処理方法
JP3590183B2 (ja) 乳房画像表示装置
JPH09106448A (ja) 異常陰影候補の検出方法
JPH11155849A (ja) 骨画像情報取得方法
JPH08287230A (ja) 計算機支援画像診断装置
JP3824749B2 (ja) 骨部画像処理方法および装置
CN116433695A (zh) 一种乳腺钼靶图像的乳腺区域提取方法及系统
JP3807541B2 (ja) 乳房画像表示方法および装置
JP2001325583A (ja) 画像処理方法および画像処理装置
JP2022052210A (ja) 情報処理装置、情報処理方法及びプログラム
JPH10108859A (ja) 異常陰影候補の検出方法および装置
JP3696324B2 (ja) 画像処理方法
JP3690844B2 (ja) 画像処理方法および装置

Legal Events

Date Code Title Description
A300 Application deemed to be withdrawn because no request for examination was validly filed

Free format text: JAPANESE INTERMEDIATE CODE: A300

Effective date: 20031202