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
Application number
JP1306569A
Other languages
English (en)
Other versions
JP2993687B2 (ja
Inventor
Yoshihiro Goto
良洋 後藤
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.)
Hitachi Healthcare Manufacturing Ltd
Original Assignee
Hitachi Medical Corp
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 Hitachi Medical Corp filed Critical Hitachi Medical Corp
Priority to JP1306569A priority Critical patent/JP2993687B2/ja
Publication of JPH03167677A publication Critical patent/JPH03167677A/ja
Application granted granted Critical
Publication of JP2993687B2 publication Critical patent/JP2993687B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

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図(イ)のA−B切断線上での脳部画像52を抽
出したい。そのために、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のみの抽出は不可能と
なる。
本発明の目的は、不必要な領域は除外して、目的とする
臓器領域のみを抽出可能とする臓器領域抽出方法を提供
するものである。
[R題を解決するための手段] 本発明は、画素の抽出の際に、その画素だけで?<r#
Jtmの画素の情報を考慮するようにした。
[作用コ 本発明によれば、画素の抽出の際に、その画素だけでな
く、周囲の画素の情報を考慮する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は抽出対象外として扱う
ステップ3では、PX+yに隣接する4点の画素を含む
演算を行い、その演算結果から抽出条件を満足するか否
かをチェックし、満足していれば、?!Px+yは臓器
画素として抽出する.満足していなければ、ステップ5
に移り,抽出対象外とする。
以上の第1図の処理は、全画素に対して行う。
勿論、抽出領域をしぼっておけば、その領域内だけの全
画素について第1図の処理をすればよい。
第3図は、第1図のステップ2,3を具体化した実施例
図である。ステップ2は、ステップ2aと2bとより成
り、対象画素PX+YがT1<Px+Y<T,の範囲に
あるか否かをチェックする。この条件を満足すれば、ス
テップ3aに移る。
ステップ3は、ステップ38〜3dより成る。これらの
ステップは、対象画素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と
してよい。
第4図,第5図,第6図は、核磁気共鳴装置でスピンエ
コー法を利用して得た画像への、適用例を示す。
第4図は、第1次エコー画像P1〜第4次エコー画像P
4の例を示す。この第工次エコー〜第4次エコーとは、
被検体を磁場励起した場合に、繰返し計測されるエコー
であり、第1次エコーの振幅が最大で、第4次エコーの
振幅が最小である。
第5次以降のエコーも存在するが、今回の例では採用し
ていない。
座標(x,y)での第t次エコー〜第4次エコーまでの
濃度特性を第5図に示す。この濃度特性は、臓器によっ
て変る。従って、臓器によって基準となる特性(曲4!
)を予め求めておき、この基準特性と同一又は類似の特
性となる位置の画素は、この臓器を構或すると判定でき
る。
そこで、本実施例では、第6図に示す如き処理を行う。
先ず、ステップ2Cで、計測して得た画素位!(x.y
)での第1次〜第4次エコーPIx,y〜P4Xl1/
の曲線が事前設定した基準曲線と一致するか否かチェッ
クする。一致には同一又は若干の差をも含む。
一致していれば、ステップ3eに移り、隣接する4点の
画素との差の絶対値が閾値よりも小さいか否かチェック
する。この4点の各点については、それぞれ第1次〜第
4次エコーまで画素値があるため、その全エコーについ
て差分の絶対値を求める。しかし、第1次エコーのみを
代表させ、この差分を求めるようにしてもよい。
閾値より小さければ、ステップ4に移り、(X,y)点
を臓器として抽出する。
第7図,第8図は本発明の他の実施例を示す。
この実施例は、差分ではなく標準偏差値を利用した。
第7図(イ)は原画とする。この原画の任意の点?n(
その濃度をTa,とする)に対して、そのa。
点を中心とする正方形(又は矩形)(斜線で示した)を
考え、正方形の中に含まれる全位置について濃度(画素
値)を統計処理してその正方形内でのa0点の濃度の標
準偏差値σ。を求める。この標準偏差値σ。を第7図に
示す処理画像上のa。点の濃度とする。同様に、a1点
に関しても、点線で示す正方形内での81点の濃度の標
準偏差値σ、を求め、これを第7図(口)のaエ点の濃
度とする。
以上の処理を原画の全位置について行って、第7図(口
)に示す標準偏差値を画素濃度とする処理画像を得る。
次に、A−B切断線上で臓器抽出を行う例を第7図(ハ
),(二)に示す。第7図(ハ)は第7図(イ)の原画
のA−B切断線上での濃度特性図、第7図(二)は第7
図(口)の処理画像のA−B切断線上での濃度特性図を
示す。
そこで、第7図(ハ)に関しては,抽出臓器用の2つの
閾値T■,T,、第7図(二)に関しては抽出臓器用の
標準偏差値に対する2つの閾値jx+ t,を指定し、
閾値処理を行う。閾値処理とは、2つの閾値内に入って
いる濃度を取出すことである。
この結果、第7図(ホ)に示すように、第7図(ハ)か
らはP+t p2+  ・・+ Psが得られ、第7図
(二)からはQt+ 92+  93が得られる。
ここで、p,〜Pill qt〜q,の両者の論理積を
とれば、P4,PSを抽出できる。このP4s Psが
求めるべきg器領域となる(第lO図と同様な臓器抽出
例としたため)。
以上の処理手順は第8図に示されている。尚、第7図で
標準偏差の代りに濃度勾配を計算し、同様な手順を用い
てもよい。
第9図は、抽出装置の実施例を示す。抽出装置は、CP
UIO、主メモリ11.トラックボール12、高速演算
回路13、演算用メモリ15、表示メモリI4、CR’
rl6、記録媒体17,共通バス18より或る。
CPUl.Oは,主メモリ1lを利用して全体システム
の管理を行う。トランクボール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・・・高
速演算回路。 特許

Claims (1)

  1. 【特許請求の範囲】 1、原画素に対する抽出条件と、前記原画素に隣接する
    画素を含む演算結果に対する抽出条件とを同時に満す場
    合に前記原画素を臓器領域中の画素として抽出する臓器
    領域抽出方法。 2、原画素に対して第1の閾値処理をし、前記原画素に
    隣接する画素を含めて演算した結果に対する第2の閾値
    を処理し、該2つの閾値処理の両者の条件を満足する原
    画素を、臓器領域中の画素として抽出する臓器領域抽出
    方法。 3、請求項1又は2記載の原画素に隣接する画素を含む
    演算は、原画素とその周りの画素との差の絶対値を求め
    る演算とする臓器領域抽出方法。 4、請求項1又は2記載の原画素に隣接する画素を含む
    演算は、各画素付近の画素を用いての標準偏差を求める
    演算とする臓器領域抽出方法。 5、請求項1又は2記載の原画素に隣接する画素を含む
    演算は、各画素の濃度勾配を求める演算とする臓器領域
    抽出方法。 6、請求項1乃至5のいずれか1つの臓器領域抽出方法
    において、原画素に対する抽出条件として、核磁気共鳴
    法での高次のスピンエコー像まで用いた曲線比較法を含
    む臓器領域抽出方法。
JP1306569A 1989-11-28 1989-11-28 臓器領域抽出方法 Expired - Lifetime JP2993687B2 (ja)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006180987A (ja) * 2004-12-27 2006-07-13 Hitachi Medical Corp 医用画像表示装置

Cited By (1)

* Cited by examiner, † Cited by third party
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