JPH03167677A - 臓器領域抽出方法 - Google Patents
臓器領域抽出方法Info
- Publication number
- JPH03167677A JPH03167677A JP1306569A JP30656989A JPH03167677A JP H03167677 A JPH03167677 A JP H03167677A JP 1306569 A JP1306569 A JP 1306569A JP 30656989 A JP30656989 A JP 30656989A JP H03167677 A JPH03167677 A JP H03167677A
- Authority
- JP
- Japan
- Prior art keywords
- pixel
- calculation
- original pixel
- pixels
- organ region
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000000605 extraction Methods 0.000 claims abstract description 30
- 210000000056 organ Anatomy 0.000 claims description 25
- 239000000284 extract Substances 0.000 claims description 4
- 238000005481 NMR spectroscopy Methods 0.000 claims description 3
- 210000001835 viscera Anatomy 0.000 abstract description 2
- 230000002093 peripheral effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 10
- 210000004556 brain Anatomy 0.000 description 6
- 238000002592 echocardiography Methods 0.000 description 4
- 210000000988 bone and bone Anatomy 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 239000013256 coordination polymer Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000006187 pill Substances 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
[産業上の利用分野コ
本発明は、3次元表示方法、特に3次元表示の前段階と
しての臓器領域抽出方法に関する。
しての臓器領域抽出方法に関する。
[従来の技術コ
X@CTで脳の3次元像を構成するためには、エッジ抽
出方法を確立する必要がある。従来のエッジ等の臓器抽
出方法にあっては,原画素の画素濃度に対して閾値を設
定し、閾値との大小を比較して閾値以内ならばその画素
を抽出するやり方をとる。
出方法を確立する必要がある。従来のエッジ等の臓器抽
出方法にあっては,原画素の画素濃度に対して閾値を設
定し、閾値との大小を比較して閾値以内ならばその画素
を抽出するやり方をとる。
[発明が解決しようとする課題]
この従来例では、正確な領域抽出ができないことがある
。その例を第10図に示す。第10図(イ)は脳に対す
るXgCT像を示す。画像50が皮膚、画像51が骨、
画像52. 53が脳本体を示す。
。その例を第10図に示す。第10図(イ)は脳に対す
るXgCT像を示す。画像50が皮膚、画像51が骨、
画像52. 53が脳本体を示す。
第10図(イ)のA−B切断線上での脳部画像52を抽
出したい。そのために、A−B切断線上での画像濃度分
布を示すと第10図(口)になる。濃度の最も高い部分
Q.,Q!は骨部であり、最も低い部分Q,, Q.は
空気部分である。
出したい。そのために、A−B切断線上での画像濃度分
布を示すと第10図(口)になる。濃度の最も高い部分
Q.,Q!は骨部であり、最も低い部分Q,, Q.は
空気部分である。
第10図(口)の特性上で、扇部画像52. 53を検
出すべくこの脳部画像抽出用の2つの閾値Tエ,T2を
設定する。そして、T1<T<T,なる範囲の濃度Tを
抽出すると,第10図(ハ)となる。この結果、脳部画
像52. 53用の抽出P,, P.の他に、Pエ,P
.,P,,P,,P7も抽出されてしまう。p ,,P
4以外は不要であり、p,,p4のみの抽出は不可能と
なる。
出すべくこの脳部画像抽出用の2つの閾値Tエ,T2を
設定する。そして、T1<T<T,なる範囲の濃度Tを
抽出すると,第10図(ハ)となる。この結果、脳部画
像52. 53用の抽出P,, P.の他に、Pエ,P
.,P,,P,,P7も抽出されてしまう。p ,,P
4以外は不要であり、p,,p4のみの抽出は不可能と
なる。
本発明の目的は、不必要な領域は除外して、目的とする
臓器領域のみを抽出可能とする臓器領域抽出方法を提供
するものである。
臓器領域のみを抽出可能とする臓器領域抽出方法を提供
するものである。
[R題を解決するための手段]
本発明は、画素の抽出の際に、その画素だけで?<r#
Jtmの画素の情報を考慮するようにした。
Jtmの画素の情報を考慮するようにした。
[作用コ
本発明によれば、画素の抽出の際に、その画素だけでな
く、周囲の画素の情報を考慮する6これにより、目的と
する臓器領域のみを限定抽出する。
く、周囲の画素の情報を考慮する6これにより、目的と
する臓器領域のみを限定抽出する。
[実施例]
第1図は、本実施例の処理フローであり、第2図はその
ための画素の位置関係を示す図である。
ための画素の位置関係を示す図である。
画像1中の画素PX+Fが抽出臓器に含まれるか否かを
チェックしたいものとする。そのために,その周囲に隣
接する4つの画素px,y−■,Pうry++ +P
X”itl/t PX+11Fを選ぶ。この4つの画素
に対して、第1図では、先ず抽出対象画素Px,yに対
する抽出条件を満すか否かをチェックする(ステップ2
)6条件を満足すればステップ3に移り、満足しなけれ
ばステップ5に移り.px,yは抽出対象外として扱う
。
チェックしたいものとする。そのために,その周囲に隣
接する4つの画素px,y−■,Pうry++ +P
X”itl/t PX+11Fを選ぶ。この4つの画素
に対して、第1図では、先ず抽出対象画素Px,yに対
する抽出条件を満すか否かをチェックする(ステップ2
)6条件を満足すればステップ3に移り、満足しなけれ
ばステップ5に移り.px,yは抽出対象外として扱う
。
ステップ3では、PX+yに隣接する4点の画素を含む
演算を行い、その演算結果から抽出条件を満足するか否
かをチェックし、満足していれば、?!Px+yは臓器
画素として抽出する.満足していなければ、ステップ5
に移り,抽出対象外とする。
演算を行い、その演算結果から抽出条件を満足するか否
かをチェックし、満足していれば、?!Px+yは臓器
画素として抽出する.満足していなければ、ステップ5
に移り,抽出対象外とする。
以上の第1図の処理は、全画素に対して行う。
勿論、抽出領域をしぼっておけば、その領域内だけの全
画素について第1図の処理をすればよい。
画素について第1図の処理をすればよい。
第3図は、第1図のステップ2,3を具体化した実施例
図である。ステップ2は、ステップ2aと2bとより成
り、対象画素PX+YがT1<Px+Y<T,の範囲に
あるか否かをチェックする。この条件を満足すれば、ス
テップ3aに移る。
図である。ステップ2は、ステップ2aと2bとより成
り、対象画素PX+YがT1<Px+Y<T,の範囲に
あるか否かをチェックする。この条件を満足すれば、ス
テップ3aに移る。
ステップ3は、ステップ38〜3dより成る。これらの
ステップは、対象画素Px,yと隣接する画素との差分
の絶対値が所定の閾値に入るか否かのチェックを行い、
すべての隣接画素の差分の絶対値がその閾値内ならばス
テップ4に移り. Px,yを画素抽出する。
ステップは、対象画素Px,yと隣接する画素との差分
の絶対値が所定の閾値に入るか否かのチェックを行い、
すべての隣接画素の差分の絶対値がその閾値内ならばス
テップ4に移り. Px,yを画素抽出する。
具体的には以下の関係式を使った。
tuo< l Px+y Px+y−x l < t
i■t,。< l P x+y P x+Y+1l
< t 2x?3n< l Px+y Px+t+y
l < t,、t4。< l Px+y Px−1
ty t < t4、ここで、tlo−t41が閾値で
あり、通常は、txo=t霊0=t,。”” t4 a
+ j I. z ” j 2■=t,1=t41と
してよい。
i■t,。< l P x+y P x+Y+1l
< t 2x?3n< l Px+y Px+t+y
l < t,、t4。< l Px+y Px−1
ty t < t4、ここで、tlo−t41が閾値で
あり、通常は、txo=t霊0=t,。”” t4 a
+ j I. z ” j 2■=t,1=t41と
してよい。
第4図,第5図,第6図は、核磁気共鳴装置でスピンエ
コー法を利用して得た画像への、適用例を示す。
コー法を利用して得た画像への、適用例を示す。
第4図は、第1次エコー画像P1〜第4次エコー画像P
4の例を示す。この第工次エコー〜第4次エコーとは、
被検体を磁場励起した場合に、繰返し計測されるエコー
であり、第1次エコーの振幅が最大で、第4次エコーの
振幅が最小である。
4の例を示す。この第工次エコー〜第4次エコーとは、
被検体を磁場励起した場合に、繰返し計測されるエコー
であり、第1次エコーの振幅が最大で、第4次エコーの
振幅が最小である。
第5次以降のエコーも存在するが、今回の例では採用し
ていない。
ていない。
座標(x,y)での第t次エコー〜第4次エコーまでの
濃度特性を第5図に示す。この濃度特性は、臓器によっ
て変る。従って、臓器によって基準となる特性(曲4!
)を予め求めておき、この基準特性と同一又は類似の特
性となる位置の画素は、この臓器を構或すると判定でき
る。
濃度特性を第5図に示す。この濃度特性は、臓器によっ
て変る。従って、臓器によって基準となる特性(曲4!
)を予め求めておき、この基準特性と同一又は類似の特
性となる位置の画素は、この臓器を構或すると判定でき
る。
そこで、本実施例では、第6図に示す如き処理を行う。
先ず、ステップ2Cで、計測して得た画素位!(x.y
)での第1次〜第4次エコーPIx,y〜P4Xl1/
の曲線が事前設定した基準曲線と一致するか否かチェッ
クする。一致には同一又は若干の差をも含む。
)での第1次〜第4次エコーPIx,y〜P4Xl1/
の曲線が事前設定した基準曲線と一致するか否かチェッ
クする。一致には同一又は若干の差をも含む。
一致していれば、ステップ3eに移り、隣接する4点の
画素との差の絶対値が閾値よりも小さいか否かチェック
する。この4点の各点については、それぞれ第1次〜第
4次エコーまで画素値があるため、その全エコーについ
て差分の絶対値を求める。しかし、第1次エコーのみを
代表させ、この差分を求めるようにしてもよい。
画素との差の絶対値が閾値よりも小さいか否かチェック
する。この4点の各点については、それぞれ第1次〜第
4次エコーまで画素値があるため、その全エコーについ
て差分の絶対値を求める。しかし、第1次エコーのみを
代表させ、この差分を求めるようにしてもよい。
閾値より小さければ、ステップ4に移り、(X,y)点
を臓器として抽出する。
を臓器として抽出する。
第7図,第8図は本発明の他の実施例を示す。
この実施例は、差分ではなく標準偏差値を利用した。
第7図(イ)は原画とする。この原画の任意の点?n(
その濃度をTa,とする)に対して、そのa。
その濃度をTa,とする)に対して、そのa。
点を中心とする正方形(又は矩形)(斜線で示した)を
考え、正方形の中に含まれる全位置について濃度(画素
値)を統計処理してその正方形内でのa0点の濃度の標
準偏差値σ。を求める。この標準偏差値σ。を第7図に
示す処理画像上のa。点の濃度とする。同様に、a1点
に関しても、点線で示す正方形内での81点の濃度の標
準偏差値σ、を求め、これを第7図(口)のaエ点の濃
度とする。
考え、正方形の中に含まれる全位置について濃度(画素
値)を統計処理してその正方形内でのa0点の濃度の標
準偏差値σ。を求める。この標準偏差値σ。を第7図に
示す処理画像上のa。点の濃度とする。同様に、a1点
に関しても、点線で示す正方形内での81点の濃度の標
準偏差値σ、を求め、これを第7図(口)のaエ点の濃
度とする。
以上の処理を原画の全位置について行って、第7図(口
)に示す標準偏差値を画素濃度とする処理画像を得る。
)に示す標準偏差値を画素濃度とする処理画像を得る。
次に、A−B切断線上で臓器抽出を行う例を第7図(ハ
),(二)に示す。第7図(ハ)は第7図(イ)の原画
のA−B切断線上での濃度特性図、第7図(二)は第7
図(口)の処理画像のA−B切断線上での濃度特性図を
示す。
),(二)に示す。第7図(ハ)は第7図(イ)の原画
のA−B切断線上での濃度特性図、第7図(二)は第7
図(口)の処理画像のA−B切断線上での濃度特性図を
示す。
そこで、第7図(ハ)に関しては,抽出臓器用の2つの
閾値T■,T,、第7図(二)に関しては抽出臓器用の
標準偏差値に対する2つの閾値jx+ t,を指定し、
閾値処理を行う。閾値処理とは、2つの閾値内に入って
いる濃度を取出すことである。
閾値T■,T,、第7図(二)に関しては抽出臓器用の
標準偏差値に対する2つの閾値jx+ t,を指定し、
閾値処理を行う。閾値処理とは、2つの閾値内に入って
いる濃度を取出すことである。
この結果、第7図(ホ)に示すように、第7図(ハ)か
らはP+t p2+ ・・+ Psが得られ、第7図
(二)からはQt+ 92+ 93が得られる。
らはP+t p2+ ・・+ Psが得られ、第7図
(二)からはQt+ 92+ 93が得られる。
ここで、p,〜Pill qt〜q,の両者の論理積を
とれば、P4,PSを抽出できる。このP4s Psが
求めるべきg器領域となる(第lO図と同様な臓器抽出
例としたため)。
とれば、P4,PSを抽出できる。このP4s Psが
求めるべきg器領域となる(第lO図と同様な臓器抽出
例としたため)。
以上の処理手順は第8図に示されている。尚、第7図で
標準偏差の代りに濃度勾配を計算し、同様な手順を用い
てもよい。
標準偏差の代りに濃度勾配を計算し、同様な手順を用い
てもよい。
第9図は、抽出装置の実施例を示す。抽出装置は、CP
UIO、主メモリ11.トラックボール12、高速演算
回路13、演算用メモリ15、表示メモリI4、CR’
rl6、記録媒体17,共通バス18より或る。
UIO、主メモリ11.トラックボール12、高速演算
回路13、演算用メモリ15、表示メモリI4、CR’
rl6、記録媒体17,共通バス18より或る。
CPUl.Oは,主メモリ1lを利用して全体システム
の管理を行う。トランクボールl2は各種の操作指示を
行う。高速演算回路l3は、3次元画像処理用の専用処
理をメモリl5を利用して行う。表示メモリ14は、3
次元画像やX線CT像や本実施例で扱った各種画像を格
納し,CRT16がこれを表示する。
の管理を行う。トランクボールl2は各種の操作指示を
行う。高速演算回路l3は、3次元画像処理用の専用処
理をメモリl5を利用して行う。表示メモリ14は、3
次元画像やX線CT像や本実施例で扱った各種画像を格
納し,CRT16がこれを表示する。
記録媒体l7は各種CT画像や3次元画像を格納する。
以上の実施例で、臓器抽出処理は、高速演算回路13に
行わせる。
行わせる。
[発明の効果]
本発明によれば、抽出すべきか否かを判定する際に、そ
の画素だけでなく周囲の画素の情報も考慮に入れるため
、目的とする臓器領域のみの限定抽出が可能となる。
の画素だけでなく周囲の画素の情報も考慮に入れるため
、目的とする臓器領域のみの限定抽出が可能となる。
第■図は本発明の処理手順を示す図,第2図は抽出画素
とその隣接画素とを示す図、第3図は画素の具体的抽出
手順例図、第4図はNMRでのスピンエコーによる画像
例図、第5図はその濃度特性図、第6図は第5図の例で
の抽出手順を示す図、第7図は標準偏差による抽出例を
示す図、第8図は第7図の抽出手順を示す図、第9図は
抽出装置の実施例図、第10図は従来例図である。 lO・・・cpu.tt・・・主メモリ、13・・・高
速演算回路。 特許
とその隣接画素とを示す図、第3図は画素の具体的抽出
手順例図、第4図はNMRでのスピンエコーによる画像
例図、第5図はその濃度特性図、第6図は第5図の例で
の抽出手順を示す図、第7図は標準偏差による抽出例を
示す図、第8図は第7図の抽出手順を示す図、第9図は
抽出装置の実施例図、第10図は従来例図である。 lO・・・cpu.tt・・・主メモリ、13・・・高
速演算回路。 特許
Claims (1)
- 【特許請求の範囲】 1、原画素に対する抽出条件と、前記原画素に隣接する
画素を含む演算結果に対する抽出条件とを同時に満す場
合に前記原画素を臓器領域中の画素として抽出する臓器
領域抽出方法。 2、原画素に対して第1の閾値処理をし、前記原画素に
隣接する画素を含めて演算した結果に対する第2の閾値
を処理し、該2つの閾値処理の両者の条件を満足する原
画素を、臓器領域中の画素として抽出する臓器領域抽出
方法。 3、請求項1又は2記載の原画素に隣接する画素を含む
演算は、原画素とその周りの画素との差の絶対値を求め
る演算とする臓器領域抽出方法。 4、請求項1又は2記載の原画素に隣接する画素を含む
演算は、各画素付近の画素を用いての標準偏差を求める
演算とする臓器領域抽出方法。 5、請求項1又は2記載の原画素に隣接する画素を含む
演算は、各画素の濃度勾配を求める演算とする臓器領域
抽出方法。 6、請求項1乃至5のいずれか1つの臓器領域抽出方法
において、原画素に対する抽出条件として、核磁気共鳴
法での高次のスピンエコー像まで用いた曲線比較法を含
む臓器領域抽出方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1306569A JP2993687B2 (ja) | 1989-11-28 | 1989-11-28 | 臓器領域抽出方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1306569A JP2993687B2 (ja) | 1989-11-28 | 1989-11-28 | 臓器領域抽出方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH03167677A true JPH03167677A (ja) | 1991-07-19 |
| JP2993687B2 JP2993687B2 (ja) | 1999-12-20 |
Family
ID=17958638
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP1306569A Expired - Lifetime JP2993687B2 (ja) | 1989-11-28 | 1989-11-28 | 臓器領域抽出方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2993687B2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5748801A (en) * | 1993-04-22 | 1998-05-05 | Hitachi Medical Corporation | Method of setting threshold values for extracting image data |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006180987A (ja) * | 2004-12-27 | 2006-07-13 | Hitachi Medical Corp | 医用画像表示装置 |
-
1989
- 1989-11-28 JP JP1306569A patent/JP2993687B2/ja not_active Expired - Lifetime
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5748801A (en) * | 1993-04-22 | 1998-05-05 | Hitachi Medical Corporation | Method of setting threshold values for extracting image data |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2993687B2 (ja) | 1999-12-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113223015B (zh) | 血管壁图像分割方法、装置、计算机设备和存储介质 | |
| US20070280518A1 (en) | Method and Apparatus for Identifying Pathology in Brain Images | |
| US20100290689A1 (en) | Discriminating infarcts from artifacts in mri scan data | |
| Fan et al. | An automatic estimation of arterial input function based on multi-stream 3D CNN | |
| WO2023124830A1 (zh) | 一种基于中心线提取的血管壁影像自动曲面重建方法 | |
| US8229189B2 (en) | Visual enhancement of interval changes using temporal subtraction, convolving, and non-rigid transformation field mapping | |
| JP2000308628A (ja) | 曲率を用いた表面のレジスタリング方法 | |
| Chen et al. | Automatic brain extraction from 3D fetal MR image with deep learning-based multi-step framework | |
| US20090252418A1 (en) | Detection of edges in an image | |
| CN112785541A (zh) | 一种ct和核磁共振医学图像融合方法 | |
| TWI679652B (zh) | 個人化大腦影像評估之方法、非暫時性電腦可讀媒體及設備 | |
| Kim et al. | Automatic localization of anatomical landmarks in cardiac MR perfusion using random forests | |
| US8224046B2 (en) | Visual enhancement of interval changes using rigid and non-rigid transformations | |
| US20080193005A1 (en) | Visual Enhancement of Interval Changes Using Temporal Subtraction | |
| CN120112936A (zh) | 用于拼接图像的方法及其设备 | |
| Liu et al. | Multiple sclerosis medical image analysis and information management | |
| JPH03167677A (ja) | 臓器領域抽出方法 | |
| JP2024161515A (ja) | 放射線治療における解剖学的変化の自動定性的記述 | |
| Han et al. | Automatic registration of brain magnetic resonance images based on Talairach reference system | |
| CN118918131A (zh) | 用于医学影像的肿瘤边缘的确定方法和装置 | |
| CN116862930A (zh) | 适用于多模态的脑血管分割方法、装置、设备和存储介质 | |
| CN115619803A (zh) | 一种图像分割方法、系统及电子设备 | |
| Mo et al. | Focusing intracranial aneurysm lesion segmentation by graph Mask2Former with local refinement in DSA images | |
| CN118823751B (zh) | 一种智能硬件的药品识别方法及装置 | |
| CN111080572A (zh) | 白质高信号定位方法、装置、设备和存储介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091022 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091022 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101022 Year of fee payment: 11 |
|
| EXPY | Cancellation because of completion of term | ||
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101022 Year of fee payment: 11 |