JP3647970B2 - Region extraction device - Google Patents
Region extraction device Download PDFInfo
- Publication number
- JP3647970B2 JP3647970B2 JP11648596A JP11648596A JP3647970B2 JP 3647970 B2 JP3647970 B2 JP 3647970B2 JP 11648596 A JP11648596 A JP 11648596A JP 11648596 A JP11648596 A JP 11648596A JP 3647970 B2 JP3647970 B2 JP 3647970B2
- Authority
- JP
- Japan
- Prior art keywords
- region
- threshold
- contour
- time
- threshold value
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、超音波診断装置やＭＲＩ等により得られた、被検体（主に人体）の臓器の断層像を表わす画像データに基づいて臓器の領域もしくは臓器内の組織の領域を抽出する領域抽出装置に関する。
【０００２】
【従来の技術】
従来より、被検体、特に人体内に超音波を送信し被検体内の組織で反射して戻ってきた超音波を受信して受信信号を得、この受信信号に基づいて被検体内の断層像を表示することにより内臓等の疾患の診断に供する超音波診断装置が用いられており、近年では、疾患の診断をより容易ならしめるために、心臓の断層像に現れた左心室の領域を抽出して表示し、あるいはその左心室の面積や、その左心室の大きさから推定される左心室の体積、あるいはそれらの時間的変化を表示することが要求されている。
【０００３】
図１１は、超音波診断装置による心臓の断層像の表示例を表わした図である。この図１１に示すように、心筋部と血流部とでは画像データのデータ値の平均値が大きく異なり、したがって、基本的な処理として所定の閾値を用いてその画像データを二値化処理することにより、心筋部と血流部とを分離し、血流部の一部分である左心室の領域を求めることができる。
【０００４】
図１１に示すような１つの断層像を表わす画像データに基づいてその断層像上にあらわれた左心室の領域をを求めるに際し、その領域をより高精度に求めるために、単なる二値化処理に加え種々の工夫が提案されている（例えば、特願平６−２９７０７１号、特願平７−２８５０５３号参照）。
【０００５】
【発明が解決しようとする課題】
図１２、図１３は、閾値を用いた二値化処理を含む処理により求められた左心室の輪郭を表わした図である。
画像データと閾値との比較により左心室の領域を抽出する場合、閾値の値が適当でないと図１２のように右心室の領域まで左心室の領域とみなされ本来よりも大きな領域として認識されてしまう場合がある。また、図１３のように本来より小さな領域として認識されてしまう場合もある。
【０００６】
本発明は、上記事情に鑑み、適切な閾値を用いて領域を高精度に抽出することのできる領域抽出装置を提供することを目的とする。
【０００７】
【課題を解決するための手段】
上記目的を達成する本発明の領域抽出装置は、所定のフレーム時間間隔毎に順次得られた、所望の臓器もしくは組織を含む被検体の断層像を担持する画像データに基づいて、上記断層像上にあらわれた臓器もしくは組織の領域を抽出し、その領域の面積もしくはその領域から推定される臓器もしくは組織の体積の時間変化を求める領域抽出装置において、
上記画像データに所定の第１の閾値を用いた閾値処理を含む処理を施すことにより、その画像データが担持する断層像上にあらわれた所望の臓器もしくは組織の輪郭を抽出する抽出手段と、
抽出手段により抽出された輪郭に囲まれた領域の面積もしくはその領域から推定される臓器もしくは組織の体積の時間的変化を表わす時系列データを求める演算手段と、
演算手段で求められた時系列データを平滑化することにより平滑化データを求める平滑化手段と、
上記時系列データと上記平滑化データとの差分が第２の閾値を越えるか否かを判定する判定手段とを備え、
上記抽出手段が、上記判定手段により上記差分が第２の閾値を越えた旨判定された場合に、第１の閾値を変更して上記輪郭を抽出し直すものであることを特徴とする。
【０００８】
ここで、上記本発明の領域抽出装置において、上記判定手段が、前記差分と複数段階の第２の閾値のうちのいずれの第２の閾値を越えるか、もしくはいずれの第２の閾値も越えないかを判定するものであって、上記抽出手段が、複数段階の第２の閾値のうち上記差分が越えた第２の閾値に応じて前記第１の閾値を用いて上記輪郭を抽出し直すものであってもよい。
【０００９】
本発明の領域抽出装置は、時系列に順次得られる断層像のうち近隣の時刻に得られた断層像上にあらわれられる所望の臓器もしくは組織の領域の大きさはほぼ近似しているという点に着目したものであり、上記差分と第２の閾値とを比較してその差分が第２の閾値を越えたとき、すなわち近隣の断層像上の領域の面積ないし体積とは極端に異なる面積ないし体積が求められた場合に、その領域の輪郭を求めるための第１の閾値を変更して輪郭を再度求め直すものであるため、最終的には、過大もしくは過小の領域として認識されることが防止され、適切な第１の閾値により適切な広さの領域が抽出される。
【００１０】
【発明の実施の形態】
以下、本発明の実施形態について説明する。
図１は、本発明の領域抽出装置の一実施形態の構成を示すブロック図である。例えば超音波診断装置やＭＲＩ等の診断装置１０で、左心室の断層面全域を内部に含む心臓の断層像を担持する画像データが、所望のフレーム時間間隔毎に得られ、得られた画像データが、本発明の一実施形態としての領域抽出装置２０に入力されて、画像メモリ２１に一旦格納される。
【００１１】
画像メモリ２１に一旦格納された画像データは画像メモリ２１から読み出されて抽出手段２２に入力される。この抽出手段２２では、送られてきた画像データが担持する心臓の断層像上にあらわれた左心室の輪郭が抽出される。
この輪郭の抽出は、基本的には、送られてきた画像データを所定の第１の閾値で二値化することにより行なわれるが、二値化処理のみでなく、例えば前掲の提案等を採用し、できるだけ正確に左心室を抽出することが望ましい。
【００１２】
抽出手段１により抽出された左心室の輪郭を表わすデータは演算手段２３に入力される。この演算手段２３は、左心室の輪郭に囲まれた領域の面積、もしくは、例えばその領域を所定の中心軸のまわりに回転させることにより定められる回転体の体積を表わすデータが求められる。この面積ないし体積を表わすデータは、診断装置１０で所望のフレーム時間間隔毎に順次得られる画像データそれぞれに対応して求められ、全体として面積ないし体積の時間的変化を表わす時系列データが求められる。
【００１３】
演算手段２３で求められた時系列データは平滑化手段２４に入力される。平滑化手段２４では、入力された時系列データが、例えば移動平均法により平滑化され、これにより平滑化データが求められる。
演算手段２３で求められた時系列データおよび平滑化手段２４で求められた平滑化データは、判定手段２５に入力される。この判定手段２５では、時系列データと平滑化データとの、各断層像毎の差分が所定の第２の閾値を越えるか否かが判定される。この判定結果は抽出手段２２にフィードバックされ、抽出手段２２では、判定手段２５により上記差分が第２の閾値を越えた旨判定された場合に、上記第１の閾値が変更されてその変更された第１の閾値を用いて左心室の輪郭の再抽出が行なわれる。
【００１４】
ここで、判定手段２５において、上記差分と複数段階の第２の閾値それぞれとを比較し、抽出手段において、複数段階の第２の閾値のうち上記差分が越えた第２の閾値に応じた第１の閾値を用いて左心室の輪郭を抽出し直してもよい。
画像メモリ２１に一旦格納された画像データ、抽出手段２２で抽出された左心室の輪郭を表わすデータ、演算手段２３で求められた時系列データ、平滑化手段２４で求められた平滑化データ、および判定手段２５による判定結果は、必要に応じて表示装置３０に送られ、その表示装置３０の図示しない表示画面上に表示される。
【００１５】
尚、抽出手段２２、演算手段２３、平滑化手段２４および判定手段２５は、それぞれの機能を満足するハードウェアで構成されていてもよく、あるいは、この領域抽出装置２０がＣＰＵを搭載しており、それらの手段２２〜２５のうちの１つもしくは複数、もしくは全部が、そのＣＰＵで実行されるソフトウェアで構成されていてもよい。
【００１６】
以下、具体的な各実施形態について説明する。
図２，図３は、本発明の具体的な実施形態の中の第１実施形態における、それぞれ、抽出手段２２による左心室輪郭を最初に抽出した時点、抽出し直した時点の、断層像の時系列変化（ａ）と輪郭内領域の面積の時系列変化（ｂ）を示す図である。
【００１７】
左心室輪郭部を最初に抽出した時点の状態を示す図２において、図２（ａ）に示すように、抽出手段２２で左心室輪郭を抽出するための第１の閾値として、時系列に配列されたいずれの断層像についても、閾値ａが採用されている。
図１に示す抽出手段２３においてこれらの断層像それぞれから求められた輪郭内領域の面積は、例えば図２（ｂ）の瞬時データ（本発明にいう時系列データ）のグラフに示すようなものとなる。
【００１８】
平滑化手段２４では、この瞬時データを以下の移動平均法を用いて平滑化することにより、図２（ｂ）に示す平滑化データを得る。
y（t ）を、時刻ｔにおける瞬時データ、すなわち各断層像毎に求められた輪郭内領域の面積を表わすデータであるとし、その平滑化データをy'(t）としたとき、移動平均法である。
【００１９】
【数１】
但し、ｋi （ｉ＝０，１，……，２ｎ）は各定数である。
により、平滑化データが求められる。
判定手段２５では、第２の閾値として、図２（ｂ）に示すように、閾値差±αが採用され、瞬時データと平滑化データとの差分が±αの範囲を越えるか否かが判定される。
【００２１】
抽出手段２２で第１の閾値として採用した閾値ａは適切な閾値ではなく、右心室まで輪郭内に取り込んだ瞬時データが求められ、その瞬時データは±αの範囲を越えている。そこで、判定手段２５によ、±αの範囲を越えたという判定結果が抽出手段２２に伝達され、抽出手段２２ではその判定結果を受けて、図３（ａ）に示すように、第１の閾値を閾値ｂに変更し、全ての断層像について左心室輪郭の再抽出が行なわれる。その結果、図３（ｂ）に示すように、全ての断層像について、瞬時データと平滑化データとの差分が±αの領域内に入るデータが得られている。
【００２２】
図３，図４は、本発明の具体的な実施形態中の第２実施形態における、それぞれ図１，図２と同様の図である。
この第２実施形態では、図５（ａ）に示すように、瞬時データと平滑化データとの差分が±αの領域から外れたときの断層像のみ、抽出手段２２における輪郭抽出のための第１の閾値を、閾値ｃに変更している。このように、過大な、もしくは過小な輪郭が求められた断層像のみ、閾値を変更して輪郭を抽出し直してもよい。
【００２３】
図６，図７は、本発明の具体的な実施形態中の第３実施形態における、それぞれ図２，図３と同様な図である。
この第３実施形態では、判定手段２５において、第２の閾値として、図６（ｂ）に示すように２段階の閾値α，βが用いられており、ある１つの断層像についての瞬時データは閾値αのみ越え、別の１つの断層像についての瞬時データは閾値βを越えている。
【００２４】
そこで、抽出手段２２では、閾値αのみ越えた断層像については、輪郭抽出のための閾値を閾値ａから少しだけ値の異なる閾値ｂに変更し、閾値βを越えた断層像については閾値αから値が大きくなる閾値ｃに変更して、それらの断層像について輪郭の再抽出を行なっている。このように、判定手段２５において複数段階の閾値を採用し、どの閾値を越えたか応じて輪郭抽出の際の閾値を変更してもよい。
【００２５】
図８〜図１０は、本発明の具体的な実施形態中の第４実施形態における、それぞれ、第１回目の輪郭抽出時点、第２回目の輪郭抽出時点、第３回目の輪郭抽出時点における断層像（輪郭）の時系列変化（ａ）と輪郭内領域の面積の時系列変化（ｂ）を示す図である。
第１回の輪郭抽出時点では、判定手段２５では、図８（ｂ）に示すように、値の大きな閾値βが用いられている。このとき、１つの断層像に関し、瞬時データと平滑化データとの差分が閾値βを越えている。
【００２６】
そこで、次に抽出手段２２において、図９（ａ）に示すように、その閾値βを越えた断層像について、輪郭抽出の閾値を、閾値βに対応した閾値ｂに変更して輪郭抽出が行なわれ、判定手段２５では、今度は、図９（ｂ）に示すように閾値βよりも小さな値の閾値αを用いて判定が行なわれる。今度は別の断層像が閾値αを越えたので、今度は、抽出手段２２では、図１０（ａ）に示すように、閾値αを越えた断層像について、輪郭抽出の閾値を、閾値αに対応した閾値ｃに変更して輪郭抽出が行なわれる。このように、判定手段２５では、一回には１つの閾値のみを用いて判定を行ない、その判定を、閾値を変更しながら繰り返してもよい。この場合に、無限に繰り返されるのを防ぐために、繰り返し回数の上限を定めておくことが好ましい。
【００２７】
以上では、抽出された輪郭に囲まれた領域の面積が過大であった場合を例に挙げて説明したが、過小の場合についても同様に本発明を適用することができる。また、上記では、輪郭に囲まれた領域の面積を求める場合について説明したが、面積に代えて例えば輪郭に囲まれた領域を断面にもつ回転体等を仮定しその回転体の体積を求めてもよい。
【００２８】
さらに、上記では超音波診断装置で得られた心臓の断層像を想定した説明を行なったが、本発明は、超音波診断装置以外の、例えばＭＲＩ等で得られた断層像にも適用することができ、心臓以外の他の臓器等にも適用することができる。
【００２９】
【発明の効果】
以上説明したように、本発明によれば、輪郭を抽出しようとする所望の臓器ないし組織が過大に、あるいは過小に抽出されることが防止され、高精度の抽出が実現する。
【図面の簡単な説明】
【図１】本発明の領域抽出装置の一実施形態の構成を示すブロック図である。
【図２】本発明の具体的な実施形態中の第１実施形態における、左心室輪郭を最初に抽出した時点の、断層像の時系列変化（ａ）と輪郭内領域の面積の時系列変化（ｂ）を示す図である。
【図３】本発明の具体的な実施形態中の第１実施形態における、左心室輪郭を抽出し直した時点の、断層像の時系列変化（ａ）と輪郭内領域の面積の時系列変化（ｂ）を示す図である。
【図４】本発明の具体的な実施形態中の第２実施形態における、左心室輪郭を最初に抽出した時点の、断層像の時系列変化（ａ）と輪郭内領域の面積の時系列変化（ｂ）を示す図である。
【図５】本発明の具体的な実施形態中の第２実施形態における、左心室輪郭を抽出し直した時点の、断層像の時系列変化（ａ）と輪郭内領域の面積の時系列変化（ｂ）を示す図である。
【図６】本発明の具体的な実施形態中の第３実施形態における、左心室輪郭を最初に抽出した時点の、断層像の時系列変化（ａ）と輪郭内領域の面積の時系列変化（ｂ）を示す図である。
【図７】本発明の具体的な実施形態中の第３実施形態における、左心室輪郭を抽出し直した時点の、断層像の時系列変化（ａ）と輪郭内領域の面積の時系列変化（ｂ）を示す図である。
【図８】本発明の具体的な実施形態中の第４実施形態における、第１回目の輪郭抽出時点における断層像（輪郭）の時系列変化（ａ）と輪郭内領域の面積の時系列変化（ｂ）を示す図である。
【図９】本発明の具体的な実施形態中の第４実施形態における、第２回目の輪郭抽出時点における断層像（輪郭）の時系列変化（ａ）と輪郭内領域の面積の時系列変化（ｂ）を示す図である。
【図１０】本発明の具体的な実施形態中の第４実施形態における、第３回目の輪郭抽出時点における断層像（輪郭）の時系列変化（ａ）と輪郭内領域の面積の時系列変化（ｂ）を示す図である。
【図１１】超音波診断装置による心臓の断層像の表示例を表わした図である。
【図１２】二値化処理を含む処理により求められた左心室の輪郭を表わした図である。
【図１３】二値化処理を含む処理により求められた左心室の輪郭を表わした図である。
【符号の説明】
１０診断装置
２０領域抽出装置
２１画像メモリ
２２抽出手段
２３演算手段
２４ 平滑化手段
２５ 判定手段
３０ 表示装置[0001]
BACKGROUND OF THE INVENTION
The present invention relates to region extraction for extracting an organ region or a region of a tissue in an organ based on image data representing a tomographic image of an organ of a subject (mainly a human body) obtained by an ultrasonic diagnostic apparatus or MRI. Relates to the device.
[0002]
[Prior art]
Conventionally, an ultrasonic wave transmitted to an object, particularly a human body, and reflected by a tissue in the object and received back to obtain a received signal, and a tomogram in the object based on the received signal In recent years, the left ventricular region that appears in the tomographic image of the heart is extracted to make it easier to diagnose the disease. Display of the left ventricle, the volume of the left ventricle estimated from the size of the left ventricle, or the temporal change thereof is required.
[0003]
FIG. 11 is a diagram illustrating a display example of a tomographic image of the heart by the ultrasonic diagnostic apparatus. As shown in FIG. 11, the average value of the image data differs greatly between the myocardial part and the blood flow part. Therefore, as a basic process, the image data is binarized using a predetermined threshold value. Thus, the myocardial part and the blood flow part can be separated, and the region of the left ventricle that is a part of the blood flow part can be obtained.
[0004]
When obtaining the region of the left ventricle appearing on the tomographic image based on the image data representing one tomographic image as shown in FIG. 11, in order to obtain the region with higher accuracy, a simple binarization process is performed. In addition, various devices have been proposed (see, for example, Japanese Patent Application No. 6-297071 and Japanese Patent Application No. 7-285053).
[0005]
[Problems to be solved by the invention]
12 and 13 are diagrams showing the contour of the left ventricle obtained by processing including binarization processing using a threshold value.
When the left ventricle region is extracted by comparing the image data with the threshold value, if the threshold value is not appropriate, the left ventricular region is regarded as the left ventricular region as shown in FIG. May end up. Further, as shown in FIG. 13, it may be recognized as an area smaller than the original.
[0006]
In view of the above circumstances, an object of the present invention is to provide a region extraction device that can extract a region with high accuracy using an appropriate threshold value.
[0007]
[Means for Solving the Problems]
The region extracting apparatus of the present invention that achieves the above object is based on image data carrying a tomographic image of a subject including a desired organ or tissue sequentially obtained at predetermined frame time intervals. In a region extraction device that extracts an area of an organ or tissue that appears, and obtains a temporal change in the area of the region or the volume of the organ or tissue estimated from the region,
Extraction means for extracting a contour of a desired organ or tissue appearing on a tomographic image carried by the image data by performing processing including threshold processing using a predetermined first threshold on the image data;
Calculating means for obtaining time-series data representing a time change of the area of the region surrounded by the contour extracted by the extracting means or the volume of the organ or tissue estimated from the region;
Smoothing means for obtaining smoothed data by smoothing the time-series data obtained by the computing means;
Determining means for determining whether or not a difference between the time series data and the smoothed data exceeds a second threshold;
The extraction means is characterized in that, when the determination means determines that the difference exceeds a second threshold, the first threshold is changed and the contour is extracted again.
[0008]
Here, in the region extraction apparatus according to the present invention, the determination means does not exceed any second threshold value among the difference and a plurality of second threshold values, or does not exceed any second threshold value. Wherein the extraction means re-extracts the contour using the first threshold value according to a second threshold value that exceeds the difference among the second threshold values in a plurality of stages. It may be.
[0009]
The region extracting apparatus of the present invention is that the size of a desired organ or tissue region appearing on a tomographic image obtained at a nearby time among tomographic images sequentially obtained in time series is approximately approximate. When the difference exceeds the second threshold when the difference is compared with the second threshold, that is, the area or volume that is extremely different from the area or volume of the region on the neighboring tomographic image. Is obtained, the first threshold value for obtaining the contour of the region is changed and the contour is obtained again, so that it is finally prevented from being recognized as an over or under region. Then, an area having an appropriate width is extracted based on the appropriate first threshold value.
[0010]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, embodiments of the present invention will be described.
FIG. 1 is a block diagram showing the configuration of an embodiment of the region extraction apparatus of the present invention. For example, in the
[0011]
The image data once stored in the
This contour extraction is basically performed by binarizing the sent image data with a predetermined first threshold. In addition to the binarization processing, for example, the above proposal is adopted. It is desirable to extract the left ventricle as accurately as possible.
[0012]
Data representing the contour of the left ventricle extracted by the extraction means 1 is input to the calculation means 23. The calculating means 23 obtains data representing the area of the region surrounded by the outline of the left ventricle, or the volume of the rotating body determined by, for example, rotating the region around a predetermined central axis. The data representing the area or volume is obtained corresponding to each image data sequentially obtained at every desired frame time interval by the
[0013]
The time series data obtained by the calculation means 23 is input to the smoothing means 24. In the smoothing means 24, the input time-series data is smoothed by, for example, a moving average method, thereby obtaining smoothed data.
The time series data obtained by the computing means 23 and the smoothed data obtained by the smoothing means 24 are input to the determining means 25. In this determination means 25, it is determined whether or not the difference between the time series data and the smoothed data for each tomographic image exceeds a predetermined second threshold value. The determination result is fed back to the
[0014]
Here, the determination means 25 compares the difference with each of the second threshold values in a plurality of stages, and the extraction means determines the second threshold value corresponding to the second threshold value exceeding the difference among the second threshold values in the plurality of stages. A left ventricular contour may be re-extracted using a threshold value of one.
Image data once stored in the
[0015]
The extraction means 22, the calculation means 23, the smoothing means 24, and the determination means 25 may be configured by hardware satisfying the respective functions, or the area extraction device 20 is equipped with a CPU. One, a plurality, or all of the
[0016]
Hereinafter, specific embodiments will be described.
2 and 3 show the tomographic images at the time when the left ventricular contour is first extracted and re-extracted by the extraction means 22 in the first embodiment among the specific embodiments of the present invention. It is a figure which shows the time series change (a) and the time series change (b) of the area of the area | region in an outline.
[0017]
In FIG. 2 showing the state at the time when the left ventricular contour is first extracted, as shown in FIG. 2 (a), the first means for extracting the left ventricular contour by the extracting
The area of the contour region obtained from each of these tomographic images in the extracting means 23 shown in FIG. 1 is as shown in the graph of instantaneous data (time series data referred to in the present invention) in FIG. Become.
[0018]
The smoothing means 24 smoothes the instantaneous data using the following moving average method to obtain the smoothed data shown in FIG.
When y (t) is instantaneous data at time t, that is, data representing the area of the region in the contour obtained for each tomographic image, and the smoothed data is y ′ (t), the moving average method It is.
[0019]
[Expression 1]
However, k i (i = 0, 1,..., 2n) is a constant.
Thus, smoothed data is obtained.
In the determination means 25, as shown in FIG. 2B, a threshold difference ± α is adopted as the second threshold, and it is determined whether or not the difference between the instantaneous data and the smoothed data exceeds the range of ± α. Is done.
[0021]
The threshold value a adopted as the first threshold value by the extracting
[0022]
3 and 4 are views similar to FIGS. 1 and 2, respectively, in the second embodiment among the specific embodiments of the present invention.
In the second embodiment, as shown in FIG. 5A, only the tomographic image when the difference between the instantaneous data and the smoothed data deviates from the range of ± α is used for the contour extraction in the extracting
[0023]
6 and 7 are views similar to FIGS. 2 and 3, respectively, in the third embodiment among the specific embodiments of the present invention.
In the third embodiment, the determination means 25 uses two-stage threshold values α and β as the second threshold value as shown in FIG. 6B, and the instantaneous data for a certain tomographic image is as follows. Only the threshold value α is exceeded, and the instantaneous data for another tomographic image exceeds the threshold value β.
[0024]
Therefore, the extraction means 22 changes the threshold for contour extraction from the threshold value a to a slightly different threshold value b for the tomographic image exceeding only the threshold value α, and from the threshold value α for the tomographic image exceeding the threshold value β. The threshold value c is increased so that the contour is re-extracted for those tomographic images. As described above, the determination unit 25 may employ a plurality of threshold levels, and may change the threshold for contour extraction depending on which threshold is exceeded.
[0025]
8 to 10 show faults at the first contour extraction time point, the second contour extraction time point, and the third contour extraction time point, respectively, in the fourth embodiment of the specific embodiments of the present invention. It is a figure which shows the time series change (a) of an image (contour), and the time series change (b) of the area of the area | region in an outline.
At the first contour extraction time point, the determination means 25 uses a threshold value β having a large value as shown in FIG. At this time, for one tomographic image, the difference between the instantaneous data and the smoothed data exceeds the threshold value β.
[0026]
Therefore, next, the extracting
[0027]
In the above description, the case where the area of the region surrounded by the extracted contour is excessive has been described as an example. However, the present invention can be similarly applied to an excessively small area. In the above description, the area of the region surrounded by the contour is obtained. However, instead of the area, for example, a rotating body having a cross section of the region surrounded by the contour is assumed, and the volume of the rotating body is obtained. Also good.
[0028]
Furthermore, although the above description has been made on the assumption of a tomographic image of the heart obtained by the ultrasonic diagnostic apparatus, the present invention is also applicable to a tomographic image obtained by MRI or the like other than the ultrasonic diagnostic apparatus. It can be applied to other organs other than the heart.
[0029]
【The invention's effect】
As described above, according to the present invention, it is possible to prevent a desired organ or tissue from which a contour is to be extracted from being excessively or excessively extracted, thereby realizing highly accurate extraction.
[Brief description of the drawings]
FIG. 1 is a block diagram showing a configuration of an embodiment of a region extraction device of the present invention.
FIG. 2 is a time-series change (a) of a tomographic image and a time-series change of the area of the region in the contour at the time when the left ventricle contour is first extracted in the first embodiment of the specific embodiments of the present invention. It is a figure which shows (b).
FIG. 3 is a time-series change (a) of a tomographic image and a time-series change of the area of the region in the contour at the time when the left ventricle contour is re-extracted in the first embodiment of the specific embodiments of the present invention. It is a figure which shows (b).
FIG. 4 is a time-series change (a) of a tomographic image and a time-series change of the area of the region in the contour at the time when the left ventricle contour is first extracted in the second embodiment of the specific embodiments of the present invention. It is a figure which shows (b).
FIG. 5 is a time-series change (a) of a tomographic image and a time-series change of the area of the region in the contour at the time when the left ventricle contour is re-extracted in the second embodiment of the specific embodiments of the present invention. It is a figure which shows (b).
FIG. 6 is a time-series change (a) of a tomographic image and a time-series change of the area of the region in the contour at the time when the left ventricular contour is first extracted in the third embodiment of the specific embodiments of the present invention. It is a figure which shows (b).
FIG. 7 is a time-series change (a) of a tomographic image and a time-series change of the area of the area in the contour at the time when the left ventricle contour is re-extracted in the third embodiment of the specific embodiments of the present invention. It is a figure which shows (b).
FIG. 8 is a time-series change (a) of a tomographic image (contour) and a time-series change of the area of the region in the contour at the time of the first contour extraction in the fourth embodiment among the specific embodiments of the present invention; It is a figure which shows (b).
FIG. 9 is a time-series change (a) of a tomographic image (contour) and a time-series change of the area of the region in the contour at the second contour extraction time in the fourth embodiment of the specific embodiments of the present invention; It is a figure which shows (b).
FIG. 10 is a time-series change (a) of a tomographic image (contour) and a time-series change of the area of the region in the contour at the time of the third contour extraction in the fourth embodiment among the specific embodiments of the present invention; It is a figure which shows (b).
FIG. 11 is a diagram showing a display example of a tomographic image of the heart by the ultrasonic diagnostic apparatus.
FIG. 12 is a diagram showing an outline of a left ventricle obtained by processing including binarization processing.
FIG. 13 is a diagram showing an outline of a left ventricle obtained by processing including binarization processing.
[Explanation of symbols]
10 diagnosis device 20
Claims (2)
前記画像データに所定の第１の閾値を用いた閾値処理を含む処理を施すことにより、該画像データが担持する断層像上にあらわれた前記臓器もしくは組織の輪郭を抽出する抽出手段と、
前記抽出手段により抽出された前記輪郭に囲まれた領域の面積もしくは該領域から推定される前記臓器もしくは組織の体積の時間的変化を表わす時系列データを求める演算手段と、
前記演算手段で求められた時系列データを平滑化することにより平滑化データを求める平滑化手段と、
前記時系列データと前記平滑化データとの差分が第２の閾値を越えるか否かを判定する判定手段とを備え、
前記抽出手段が、前記判定手段により前記差分が前記第２の閾値を越えた旨判定された場合に、前記第１の閾値を変更して前記輪郭を抽出し直すものであることを特徴とする領域抽出装置。Based on image data carrying a tomographic image of a subject including a desired organ or tissue sequentially obtained at predetermined frame time intervals, the region of the organ or tissue appearing on the tomographic image is extracted, In the region extraction device for obtaining the time change of the area of the region or the volume of the organ or tissue estimated from the region,
Extraction means for extracting a contour of the organ or tissue appearing on a tomographic image carried by the image data by performing processing including threshold processing using a predetermined first threshold on the image data;
Computing means for obtaining time series data representing an area of the region surrounded by the contour extracted by the extracting means or a temporal change in the volume of the organ or tissue estimated from the region;
Smoothing means for obtaining smoothed data by smoothing the time-series data obtained by the computing means;
Determining means for determining whether a difference between the time series data and the smoothed data exceeds a second threshold;
The extraction unit is configured to change the first threshold and re-extract the contour when the determination unit determines that the difference exceeds the second threshold. Region extraction device.
前記抽出手段が、複数段階の第２の閾値のうち前記差分が越えた第２の閾値に応じた前記第１の閾値を用いて前記輪郭を抽出し直すものであることを特徴とする請求項１記載の領域抽出装置。The determination means determines whether the difference exceeds any second threshold among the second thresholds of a plurality of stages, or does not exceed any second threshold;
The said extraction means re-extracts the said outline using the said 1st threshold value according to the 2nd threshold value which the said difference exceeded among the 2nd threshold values of several steps, The said outline is characterized by the above-mentioned. 1. The area extracting apparatus according to 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP11648596A JP3647970B2 (en) | 1996-05-10 | 1996-05-10 | Region extraction device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP11648596A JP3647970B2 (en) | 1996-05-10 | 1996-05-10 | Region extraction device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH09299366A JPH09299366A (en) | 1997-11-25 |
| JP3647970B2 true JP3647970B2 (en) | 2005-05-18 |
Family
ID=14688290
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP11648596A Expired - Fee Related JP3647970B2 (en) | 1996-05-10 | 1996-05-10 | Region extraction device |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3647970B2 (en) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3534009B2 (en) | 1999-09-24 | 2004-06-07 | 日本電気株式会社 | Outline extraction method and apparatus |
| JP3986866B2 (en) * | 2002-03-29 | 2007-10-03 | 松下電器産業株式会社 | Image processing apparatus and ultrasonic diagnostic apparatus |
| JP2005073817A (en) * | 2003-08-29 | 2005-03-24 | Hitachi Medical Corp | Medical image displaying apparatus |
| CN101151638B (en) * | 2005-04-01 | 2012-05-16 | 皇家飞利浦电子股份有限公司 | A method, and a system for segmenting a structure in a dataset |
| JP4327139B2 (en) * | 2005-09-05 | 2009-09-09 | ザイオソフト株式会社 | Image processing method and image processing program |
| JP4912807B2 (en) * | 2006-09-22 | 2012-04-11 | 株式会社東芝 | Ultrasound diagnostic imaging equipment |
| JP5714711B2 (en) * | 2011-09-27 | 2015-05-07 | 日本メジフィジックス株式会社 | Myocardial contour determination technology |
| CN111080584B (en) * | 2019-12-03 | 2023-10-31 | 上海联影智能医疗科技有限公司 | Quality control method for medical image, computer device and readable storage medium |
-
1996
- 1996-05-10 JP JP11648596A patent/JP3647970B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH09299366A (en) | 1997-11-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6217520B1 (en) | Diagnostic medical ultrasound system and method for object of interest extraction | |
| CN109035284B (en) | Heart CT image segmentation method, device, equipment and medium based on deep learning | |
| JP3486615B2 (en) | Region extraction method for medical images | |
| CN111667467B (en) | Clustering algorithm-based lower limb vascular calcification index multi-parameter accumulation calculation method | |
| EP2365471B1 (en) | Diagnosis assisting apparatus, coronary artery analyzing method and recording medium having a coronary artery analzying program stored therein | |
| CN103026379B (en) | The method calculating image noise level | |
| CN110838114B (en) | Pulmonary nodule detection method, device and computer storage medium | |
| CN108805923B (en) | Physiological parameter acquisition method and device | |
| JPH10229979A (en) | Method for determining epicardial boundary of myocardium and system therefor | |
| JP3194741B2 (en) | Automatic digital image quantification system | |
| CN109166104A (en) | A kind of lesion detection method, device and equipment | |
| US11751823B2 (en) | Image processing apparatus, image processing method, and program | |
| CN111862038B (en) | Plaque detection method, plaque detection device, plaque detection equipment and plaque detection medium | |
| US7689018B2 (en) | Anomaly detection in volume data structure information | |
| US20060211940A1 (en) | Blood vessel structure segmentation system and method | |
| CN112950737A (en) | Fundus fluorescence radiography image generation method based on deep learning | |
| JP3647970B2 (en) | Region extraction device | |
| JP3195136B2 (en) | Image processing device | |
| Loizou et al. | Segmentation of atherosclerotic carotid plaque in ultrasound video | |
| CN117152026B (en) | Intravascular ultrasound image processing method, device and equipment | |
| JP2001070247A (en) | Method of detecting blood vessel crossing part in image of eyeground | |
| US20230222668A1 (en) | Image processing apparatus, image processing method, and recording medium | |
| CN116091444A (en) | Side branch evaluation method and device, storage medium and terminal | |
| Benrabha et al. | Automatic ROI detection and classification of the achilles tendon ultrasound images | |
| CN110827255A (en) | Plaque stability prediction method and system based on coronary artery CT image |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050202 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20050208 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050210 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110218 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140218 Year of fee payment: 9 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |