JPH0675576B2 - Ultrasonic image display device - Google Patents

Ultrasonic image display device

Info

Publication number
JPH0675576B2
JPH0675576B2 JP3048301A JP4830191A JPH0675576B2 JP H0675576 B2 JPH0675576 B2 JP H0675576B2 JP 3048301 A JP3048301 A JP 3048301A JP 4830191 A JP4830191 A JP 4830191A JP H0675576 B2 JPH0675576 B2 JP H0675576B2
Authority
JP
Japan
Prior art keywords
echo data
dimensional
median
ultrasonic
image display
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
Application number
JP3048301A
Other languages
Japanese (ja)
Other versions
JPH04285545A (en
Inventor
正安 伊東
剛 望月
憲吉 立川
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 Ltd
Original Assignee
Aloka 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 Aloka Co Ltd filed Critical Aloka Co Ltd
Priority to JP3048301A priority Critical patent/JPH0675576B2/en
Publication of JPH04285545A publication Critical patent/JPH04285545A/en
Publication of JPH0675576B2 publication Critical patent/JPH0675576B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)
  • Image Processing (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】本発明は、超音波画像表示装置、
特に空間的なフィルタを備えた超音波画像表示装置に関
する。
BACKGROUND OF THE INVENTION The present invention relates to an ultrasonic image display device,
In particular, it relates to an ultrasonic image display device provided with a spatial filter.

【0002】[0002]

【従来の技術】超音波を生体に送受波し、生体内の断層
画像などを表示する超音波診断装置が知られており、医
療の分野において活用されている。
2. Description of the Related Art An ultrasonic diagnostic apparatus for transmitting and receiving ultrasonic waves to and from a living body and displaying a tomographic image in the living body is known and used in the medical field.

【0003】近年では、生体内三次元領域に超音波を送
受波し、取り込まれたエコーデータを用いて、前記三次
元領域の三次元画像を表示する超音波画像表示装置が提
案されており、その実用化が要望されている。
In recent years, an ultrasonic image display device has been proposed which transmits and receives ultrasonic waves to and from a three-dimensional region in a living body and uses the echo data taken in to display a three-dimensional image of the three-dimensional region. There is a demand for its practical application.

【0004】この超音波画像表示装置においては、一般
に、表示される物体の表面を抽出し、その物体表面を仮
想視点からの距離に応じて白黒濃淡処理することによ
り、前記物体の空間的な表現を行っている。
In this ultrasonic image display apparatus, generally, the surface of an object to be displayed is extracted, and the surface of the object is subjected to black-and-white shading processing according to the distance from a virtual viewpoint, thereby spatially expressing the object. It is carried out.

【0005】ところで、従来の二次元断層画像を表示す
る超音波診断装置においては、その表示すべき画像を鮮
明なものとするためいわゆるフィルタ処理が行われてお
り、そのフィルタ処理としては、例えば平滑化法や中央
値抽出法などが用いられている。
By the way, in the conventional ultrasonic diagnostic apparatus for displaying a two-dimensional tomographic image, so-called filter processing is performed in order to make the image to be displayed clear. For example, the filter processing is smoothing. The conversion method and median extraction method are used.

【0006】[0006]

【発明が解決しようとする課題】しかしながら、上述し
た超音波画像表示装置においてフィルタ処理を行う場
合、いわゆるエッジを有効に保存しつつ効果的なノイズ
除去を行うことができないという問題があった。
However, when performing the filtering process in the above-mentioned ultrasonic image display apparatus, there is a problem that so-called edges cannot be effectively saved while effective noise removal cannot be performed.

【0007】すなわち、平滑化法を適用した場合には、
ある局所領域において全体的に平均化が行われるため、
ノイズ除去効果はあるがその反面、フィルタ処理対象で
ある注目エコーデータが保存されない結果となり、これ
により画像がぼやけたり不鮮明となるという問題があっ
た。
That is, when the smoothing method is applied,
Since the averaging is performed globally in a local region,
Although there is a noise removal effect, on the other hand, the result is that the target echo data to be filtered is not saved, which causes a problem that the image becomes blurred or unclear.

【0008】一方、エッジを保存できるものとしていわ
ゆるメディアンフィルタ(中央値抽出フィルタ)が知ら
れているが、このメディアフィルタは、そのほとんどが
二次元的な処理を行うものであり、超音波三次元画像表
示装置に適用することができないという問題があった。
On the other hand, a so-called median filter (median extraction filter) is known as a device capable of preserving edges. Most of the median filters perform two-dimensional processing, and ultrasonic three-dimensional processing is performed. There is a problem that it cannot be applied to an image display device.

【0009】本発明は上記従来の課題に鑑みなされたも
のであり、その目的は、超音波三次元画像処理という特
有の画像処理において、エッジを保存しつつ効果的にノ
イズの除去を行えるフィルタを備えた超音波画像表示装
置を提供することにある。
The present invention has been made in view of the above conventional problems, and an object thereof is to provide a filter capable of effectively removing noise while preserving edges in a specific image processing called ultrasonic three-dimensional image processing. An object is to provide a provided ultrasonic image display device.

【0010】[0010]

【課題を解決するための手段】上記目的を達成するため
に、本発明は、超音波を生体内三次元領域に送受波して
取り込まれたエコーデータを格納する三次元データメモ
リを含み、前記エコーデータを用いて前記生体内三次元
領域の超音波三次元画像を形成する超音波画像表示装置
において、フィルタ処理対象である注目エコーデータを
中心として三次元に広がった局所領域に含まれるエコー
データを前記三次元データメモリから抽出する局所領域
内データ抽出手段と、前記注目エコーデータを通過する
互いに交差設定されたn個(n≧3)の参照方向それぞ
れについて、その参照方向上のエコーデータの中から中
央値を判定する個別中央値判定手段と、前記判定された
n個の中央値の中から更に代表中央値を判定する代表中
央値判定手段と、を含み、前記判定された代表中央値を
前記注目エコーデータの更新値とすることを特徴とす
る。
In order to achieve the above object, the present invention includes a three-dimensional data memory for storing echo data that is transmitted and received by transmitting and receiving ultrasonic waves to and from a three-dimensional region in a living body. In an ultrasonic image display device for forming an ultrasonic three-dimensional image of the in-vivo three-dimensional region using echo data, echo data included in a three-dimensional local region centered on the target echo data to be filtered Local area data extraction means for extracting from the three-dimensional data memory, and n (n ≧ 3) reference directions that are set to intersect each other and that pass through the echo data of interest. Individual median value determining means for determining a median value from inside, representative median value determining means for further determining a representative median value from the determined n median values, Wherein, and wherein the determined representative central value be updated value of the target echo data.

【0011】[0011]

【作用】上記構成によれば、各参照方向上にあるエコー
データの中から個別中央値が判定され、更に各個別中央
値の中から代表中央値が判定されるので、注目エコーデ
ータ近傍の中で最も標準的、換言すれば最も確からしい
値を注目エコーデータの更新値とすることができる。こ
の結果、例えばスパイクノイズなどを有効に排除して、
かつ、いわゆるエッジをぼやかさずに有効に保存するこ
とが可能となる。特に、本発明によれば、従来行われて
いた二次元的なフィルタ処理では得られない、注目エコ
ーデータの前後のスライス面(断層面)からのエコーデ
ータをも考慮して三次元的に更新値を求められるという
利点がある。
According to the above structure, the individual median value is determined from the echo data in each reference direction, and the representative median value is determined from each individual median value. The most standard value, in other words, the most probable value can be used as the updated value of the echo data of interest. As a result, for example, spike noise is effectively eliminated,
Moreover, it is possible to effectively save so-called edges without blurring. In particular, according to the present invention, three-dimensional updating is performed in consideration of echo data from slice planes (tomographic planes) before and after the echo data of interest, which cannot be obtained by the conventional two-dimensional filter processing. There is an advantage that the value can be obtained.

【0012】[0012]

【実施例】以下、本発明の好適な実施例を図面に基づい
て説明する。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT A preferred embodiment of the present invention will be described below with reference to the drawings.

【0013】図1には、本発明に係る超音波画像表示装
置10(以下、装置という)の構成が示されている。な
お、この装置10は、図においては超音波診断装置12
に外部接続されているが、超音波診断装置12の内部に
組み込むこともできる。
FIG. 1 shows the configuration of an ultrasonic image display apparatus 10 (hereinafter referred to as apparatus) according to the present invention. It should be noted that this device 10 is shown as an ultrasonic diagnostic device 12 in the figure.
Although it is externally connected to the ultrasonic diagnostic apparatus 12, it can be incorporated into the ultrasonic diagnostic apparatus 12.

【0014】超音波診断装置12においては、超音波探
触子により生体に対して超音波が送受波される。これに
より得られたエコーデータは、装置10に設けられたI
/O部14を介して、更にバス16を介して三次元デー
タメモリ18に格納される。なお、この三次元データメ
モリ18への格納は、従来の超音波画像表示装置と同様
に行われている。
In the ultrasonic diagnostic apparatus 12, ultrasonic waves are transmitted / received to / from a living body by the ultrasonic probe. The echo data obtained by this is the I data provided in the device 10.
The data is stored in the three-dimensional data memory 18 via the / O unit 14 and the bus 16. The storage in the three-dimensional data memory 18 is performed in the same manner as in the conventional ultrasonic image display device.

【0015】バス16に接続された演算処理部20は、
本実施例において三次元メディアンフィルタ22と、三
次元画像形成部24と、を含む。なお、この演算処理部
20は、例えばCPUや高速演算プロセッサなどから構
成される。
The arithmetic processing unit 20 connected to the bus 16 is
In this embodiment, the three-dimensional median filter 22 and the three-dimensional image forming unit 24 are included. The arithmetic processing unit 20 is composed of, for example, a CPU and a high-speed arithmetic processor.

【0016】三次元メディアンフィルタ22は、三次元
データメモリ18に格納された各エコーデータに対し
て、空間的なフィルタ処理を行うものであり、その詳細
については後述する。
The three-dimensional median filter 22 performs a spatial filter process on each echo data stored in the three-dimensional data memory 18, the details of which will be described later.

【0017】三次元画像形成部24は、フィルタ処理さ
れた各エコーデータを用いて生体内三次元領域の白黒濃
淡画像を形成するものである。具体的には、表示すべき
物体の表面を抽出し、その表面を仮想視点からの距離に
応じて白黒濃淡処理することにより、空間的なすなわち
疑似立体的な物体の画像を形成している。ここで、その
表面抽出処理に当たっては、ある閾値を境とした2値化
処理が行われている。前記三次元画像形成部24にて形
成された画像データは、バス16を介して、フレームメ
モリ26に一旦格納された後、CRT28に送出され
て、画像表示が行われる。
The three-dimensional image forming section 24 forms a black-and-white gray-scale image of a three-dimensional region in the living body by using each echo data which has been filtered. Specifically, the surface of the object to be displayed is extracted, and the surface is subjected to grayscale processing according to the distance from the virtual viewpoint to form a spatial or pseudo-stereoscopic image of the object. Here, in the surface extraction processing, binarization processing is performed with a certain threshold as a boundary. The image data formed by the three-dimensional image forming unit 24 is temporarily stored in the frame memory 26 via the bus 16 and then sent to the CRT 28 for image display.

【0018】次に、三次元メディアンフィルタ22につ
いて詳述する。
Next, the three-dimensional median filter 22 will be described in detail.

【0019】まず、フィルタ処理の原理について図2〜
図4を用いて説明する。
First, the principle of the filtering process will be described with reference to FIGS.
This will be described with reference to FIG.

【0020】図2には、三次元データメモリ18に格納
されたエコーデータの一部集合体が示され、図において
は、注目エコーデータXを中心として、三次元的に広が
った局所領域Vが示されている。ここで、各エコーデー
タを“ボクセル”と定義する。従って、局所領域Vは、
注目すべきボクセルXを中心として3×3×3個のボク
セルで構成されている。なお、この局所領域Vに対し
て、直交座標系X−Y−Zが定義されている。
FIG. 2 shows a partial aggregate of the echo data stored in the three-dimensional data memory 18. In the figure, a local region V that spreads three-dimensionally around the echo data X of interest is shown. It is shown. Here, each echo data is defined as a "voxel". Therefore, the local region V is
It is composed of 3 × 3 × 3 voxels centering on the notable voxel X. An orthogonal coordinate system XYZ is defined for this local region V.

【0021】このような局所領域Vに対して、本発明に
おいては、注目すべきボクセルXを中心として、互いに
交差する複数の参照方向が設定されている。本実施例に
おいては、注目すべきボクセルXを互いに交叉する13
個の参照方向が設定されている。その参照方向の例につ
いて図3及び図4を用いて説明する。
In the present invention, a plurality of reference directions intersecting each other with respect to the voxel X to be noted are set for such a local region V. In this embodiment, the voxels X to be noted are crossed with each other 13
The reference direction of each is set. An example of the reference direction will be described with reference to FIGS. 3 and 4.

【0022】図3(A)及び(B)には、2つの参照方
向100a及び100bが例示されている。参照方向1
00aは、X−Y平面内対角方向に設定され、一方、参
照方向100bは、X−Z平面内対角方向に設定されて
いる。この図3(A)及び(B)に示す各参照方向が通
過するボクセルは図4(A)及び(B)に実線で示され
ている。
Two reference directions 100a and 100b are illustrated in FIGS. 3A and 3B. Reference direction 1
00a is set diagonally in the XY plane, while reference direction 100b is set diagonally in the XY plane. The voxels through which the respective reference directions shown in FIGS. 3A and 3B pass are indicated by solid lines in FIGS. 4A and 4B.

【0023】図4において、このような参照方向が通過
するボクセル(図においては3つのボクセル)の中で、
本実施例においては、まず中央値が判定される。すなわ
ち、参照方向が縦横斜めに13方向あるとすると、13
個の中央値が判定されることになる。
In FIG. 4, among the voxels (three voxels in the figure) through which such a reference direction passes,
In this embodiment, the median value is first determined. That is, if there are 13 reference directions in the vertical and horizontal directions, then 13
The median value of each will be determined.

【0024】そして、本実施例においては、この判定さ
れた13個の中央値の中から、更に、中央値が判定され
る。この判定された中央値が、注目すべきボクセルX、
すなわち注目エコーデータの更新値として決定されるの
である。
In the present embodiment, the median value is further determined from the 13 determined median values. This determined median is the voxel X of interest,
That is, it is determined as the updated value of the echo data of interest.

【0025】ここで、その中央値判定の計算について述
べる。各ボクセルの座標を(l,m,n)とし、各参照
方向における中央値判定処理をmedi[ ]として表
すと、13個の参照方向における各中央値(個別中央
値)M1〜M13は、次の(1)〜(13)式で表され
る。
Now, the calculation of the median value judgment will be described. When the coordinates of each voxel are (l, m, n) and the median value determination process in each reference direction is represented by medi [], the median values (individual median values) M1 to M13 in the 13 reference directions are as follows. It is represented by the equations (1) to (13).

【0026】 M1 =medi[d(l-1,m+1,n),d(l,m,n),d(l+1,m-1,n)] …(1) M2 =medi[d(l-1,m-1,n),d(l,m,n),d(l+1,m+1,n)] …(2) M3 =medi[d(l-1,m,n),d(l,m,n),d(l+1,m,n)] …(3) M4 =medi[d(l,m-1,n),d(l,m,n),d(l,m+1,n)] …(4) M5 =medi[d(l,m+1,n-1),d(l,m,n),d(l,m-1,n+1)] …(5) M6 =medi[d(l,m-1,n-1),d(l,m,n),d(l,m+1,n+1)] …(6) M7 =medi[d(l,m,n-1),d(l,m,n),d(l,m,n+1)] …(7) M8 =medi[d(l-1,m,n+1),d(l,m,n),d(l+1,m,n-1)] …(8) M9 =medi[d(l+1,m,n+1),d(l,m,n),d(l-1,m,n-1)] …(9) M10 =medi[d(l-1,m+1,n+1),d(l,m,n),d(l+1,m-1,n-1)] …(10) M11 =medi[d(l+1,m+1,n-1),d(l,m,n),d(l-1,m-1,n+1)] …(11) M12 =medi[d(l+1,m+1,n+1),d(l,m,n),d(l-1,m-1,n-1)] …(12) M13 =medi[d(l-1,m+1,n-1),d(l,m,n),d(l+1,m-1,n+1)] …(13) また、この13個の中央値から更に中央値(代表中央
値)を求める場合は、次の(14)式が用いられ、この
式によりフィルタ処理としての出力値M(l,m,n)
が得られる。
M1 = medi [d (l-1, m + 1, n), d (l, m, n), d (l + 1, m-1, n)] (1) M2 = medi [ d (l-1, m-1, n), d (l, m, n), d (l + 1, m + 1, n)] (2) M3 = medi [d (l-1, m , n), d (l, m, n), d (l + 1, m, n)] (3) M4 = medi [d (l, m-1, n), d (l, m, n) ), d (l, m + 1, n)] (4) M5 = medi [d (l, m + 1, n-1), d (l, m, n), d (l, m-1) , n + 1)] (5) M6 = medi [d (l, m-1, n-1), d (l, m, n), d (l, m + 1, n + 1)] (6) M7 = medi [d (l, m, n-1), d (l, m, n), d (l, m, n + 1)] (7) M8 = medi [d (l- 1, m, n + 1), d (l, m, n), d (l + 1, m, n-1)] (8) M9 = medi [d (l + 1, m, n + 1) ), d (l, m, n), d (l-1, m, n-1)] (9) M10 = medi [d (l-1, m + 1, n + 1), d (l , m, n), d (l + 1, m-1, n-1)] (10) M11 = medi [d (l + 1, m + 1, n-1), d (l, m, n), d (l-1, m-1, n + 1)] (11) M12 = medi [d (l + 1, m + 1, n + 1), d (l, m, n), d (l-1, m-1, n-1)] (12) M13 = medi [d (l-1, m + 1, n-1), d (l, m, n), d (l + 1, m-1, n + 1)] (13) Further, when the median value (representative median value) is calculated from the 13 median values, the following formula (14) is used. Output as a filter by M (l, m, n)
Is obtained.

【0027】 M(l,m,n)=medi[M1,M2,M3, ……………M13] …(14) 以上のように、本発明においては、局所領域に対して、
注目エコーデータを中心として互いに交差する複数の参
照方向を設定し、各参照方向について個別中央値を求め
た後、その求められた複数の個別中央値の中から代表中
央値を決定するものである。なお、本実施例においては
3×3×3の局所領域Vを設定したが、これには限られ
ず、例えば5×5×5としてもよい。この場合には、参
照方向の数を13方向以上に設定してもよく、また例え
ばその処理を行うプロセッサの能力に応じて、参照方向
の数を増減しても好適である。いずれにしても、注目エ
コーデータを中心として均等に参照方向が分散されるこ
とが好適と言える。ここで、局所領域Vは、本実施例に
おいて、立方体としたが、場合によっては球状であって
もかまわない。
M (l, m, n) = medi [M1, M2, M3, ... M13] (14) As described above, in the present invention, the local region is
A plurality of reference directions intersecting each other around the echo data of interest are set, individual median values are obtained for each reference direction, and then a representative median value is determined from the obtained individual median values. . In addition, although the local region V of 3 × 3 × 3 is set in the present embodiment, the present invention is not limited to this and may be, for example, 5 × 5 × 5. In this case, the number of reference directions may be set to 13 or more, and the number of reference directions may be increased or decreased depending on, for example, the capability of the processor performing the processing. In any case, it can be said that it is preferable that the reference directions are evenly distributed around the echo data of interest. Here, although the local region V is a cube in this embodiment, it may be spherical in some cases.

【0028】次に、図5を用いて、以上説明した原理を
用いてフィルタ処理を行う場合の具体例について説明す
る。このフィルタ処理は、図1に示した三次元メディア
ンフィルタ22によって行われるものである。
Next, with reference to FIG. 5, a specific example in the case of performing the filtering process using the principle described above will be described. This filtering process is performed by the three-dimensional median filter 22 shown in FIG.

【0029】図5において、ステップ101では、三次
元データメモリ18から、フィルタ処理すべき注目エコ
ーデータを中心とした局所領域V内の各エコーデータが
抽出される。
In FIG. 5, in step 101, each echo data in the local area V centering on the echo data of interest to be filtered is extracted from the three-dimensional data memory 18.

【0030】ステップ102では、本実施例において1
3方向に設定された各参照方向における中央値、すなわ
ち13個の個別中央値が判定される。勿論、参照方向が
それ以上に設定されている場合には、参照方向と同数の
個別中央値が判定されることになる。なお、本発明にお
いて特徴とする三次元メディアンフィルタ22は、空間
的なフィルタ処理を行うものであるため、参照方向は少
なくとも3方向(互いに直交した場合を含む)以上に設
定される。
In step 102, 1 is set in this embodiment.
The median value in each reference direction set in the three directions, that is, 13 individual median values are determined. Of course, when the reference direction is set to be more than that, the same number of individual median values as the reference direction are determined. Since the three-dimensional median filter 22 that is a feature of the present invention performs spatial filter processing, the reference directions are set to at least three directions (including the case where they are orthogonal to each other).

【0031】ステップ103では、ステップ102で求
められた全ての個別中央値の中から、更に中央値が求め
られる。すなわち、代表中央値が判定される。この代表
中央値は、前記注目エコーデータのフィルタ処理後の更
新値である。
In step 103, the median value is further calculated from all the individual median values calculated in step 102. That is, the representative median value is determined. The representative median value is an updated value of the echo data of interest after filtering.

【0032】ステップ104では、前記三次元データメ
モリ18に格納された全てのエコーデータについて、こ
のようなフィルタ処理が終了したか否かが判定されてい
る。そして、終了していないと判断されるとステップ1
05に移行し、ステップ101で設定された局所領域が
1つ隣の位置へ移動し、上記過程が繰り返されることに
なる。
In step 104, it is judged whether or not such filter processing has been completed for all the echo data stored in the three-dimensional data memory 18. Then, if it is determined that the processing is not completed, step 1
The flow shifts to 05, the local area set in step 101 moves to the next adjacent position, and the above process is repeated.

【0033】一方、ステップ104で、全てのエコーデ
ータについてフィルタ処理が終了したと判断されると、
このルーチンは終了する。
On the other hand, when it is judged in step 104 that the filtering process is completed for all the echo data,
This routine ends.

【0034】この処理終了後には、三次元画像形成部2
4でフィルタ処理され、ノイズなどが排除されたエコー
データによって、鮮明な三次元画像が形成され、この形
成された三次元画像がCRT28に表示されることにな
る。
After completion of this processing, the three-dimensional image forming unit 2
A clear three-dimensional image is formed by the echo data that has been filtered by No. 4 and noise or the like has been removed, and this formed three-dimensional image is displayed on the CRT 28.

【0035】従って、本発明に係るこのような三次元メ
ディアンフィルタ22によれば、従来においては平面的
な処理でしかなかったメディアンフィルタを空間的なフ
ィルタに拡張し、特に超音波三次元画像処理において有
効なフィルタを構成することが可能となる。つまり、超
音波三次元画像においては、平面的に存在するエコーデ
ータの他にその平面前後のエコーデータも考慮してフィ
ルタ処理を行う必要があるが、このような三次元メディ
アンフィルタによれば、二次元フィルタでは欠落してし
まうエコーデータもその前後のスライス面(断層面)の
エコーデータによって復元することが可能となり、同時
に、いわゆるスパイクノイズなどを有効に除去すること
が可能となる。
Therefore, according to such a three-dimensional median filter 22 of the present invention, the median filter which has been conventionally only a planar process is expanded to a spatial filter, and particularly, an ultrasonic three-dimensional image process is performed. It becomes possible to construct an effective filter in. That is, in the ultrasonic three-dimensional image, in addition to the echo data existing in a plane, it is necessary to consider the echo data before and after that plane, but according to such a three-dimensional median filter, The echo data that is missing in the two-dimensional filter can be restored by the echo data of the slice plane (tomographic plane) before and after it, and at the same time, so-called spike noise and the like can be effectively removed.

【0036】[0036]

【発明の効果】以上説明したように、本発明によれば、
超音波三次元画像処理において効果的な空間的フィルタ
を構成することができ、欠落したエコーデータの復元や
スパイクノイズなどの排除を効果的に行うことが可能と
なる。また、たとえ画像表示を二次元的にすなわち断層
面で表示する場合においても、このような本発明に係る
フィルタ処理を行うことにより、その断層面の前後のエ
コーデータによって、表示される断層面内のエコーデー
タの補間などを行えるという利点を有する。
As described above, according to the present invention,
A spatial filter effective in ultrasonic three-dimensional image processing can be configured, and it is possible to effectively restore missing echo data and eliminate spike noise. Further, even when the image display is displayed two-dimensionally, that is, in the tomographic plane, by performing such filter processing according to the present invention, the echo data before and after the tomographic plane causes the inside of the tomographic plane to be displayed. It has an advantage that interpolation of echo data can be performed.

【図面の簡単な説明】[Brief description of drawings]

【図1】本発明に係る超音波画像表示装置の構成を示す
ブロック図である。
FIG. 1 is a block diagram showing a configuration of an ultrasonic image display device according to the present invention.

【図2】局所領域Vを示す概念図である。FIG. 2 is a conceptual diagram showing a local area V.

【図3】参照方向を示す説明図である。FIG. 3 is an explanatory diagram showing a reference direction.

【図4】参照方向に沿うボクセルを示す説明図である。FIG. 4 is an explanatory diagram showing voxels along a reference direction.

【図5】本発明に係るフィルタ処理の工程を示すフロー
チャート図である。
FIG. 5 is a flow chart diagram showing steps of a filtering process according to the present invention.

【符号の説明】[Explanation of symbols]

10 超音波画像表示装置 12 超音波診断装置 18 三次元データメモリ 22 三次元メディアンフィルタ 10 ultrasonic image display device 12 ultrasonic diagnostic device 18 three-dimensional data memory 22 three-dimensional median filter

Claims (1)

【特許請求の範囲】[Claims] 【請求項1】超音波を生体内三次元領域に送受波して取
り込まれたエコーデータを格納する三次元データメモリ
を含み、前記エコーデータを用いて前記生体内三次元領
域の超音波三次元画像又は断層画像を形成する超音波画
像表示装置において、フィルタ処理対象である注目エコ
ーデータを中心として三次元に広がった局所領域に含ま
れるエコーデータを前記三次元データメモリから抽出す
る局所領域内データ抽出手段と、前記注目エコーデータ
を通過する互いに交差設定されたn個(n≧3)の参照
方向それぞれについて、その参照方向上のエコーデータ
の中から中央値を判定する個別中央値判定手段と、前記
判定されたn個の中央値の中から更に代表中央値を判定
する代表中央値判定手段と、を含み、前記判定された代
表中央値を前記注目エコーデータの更新値とすることを
特徴とする超音波画像表示装置。
1. A three-dimensional data memory for storing echo data obtained by transmitting and receiving ultrasonic waves to and from a three-dimensional region in a living body, and using the echo data, a three-dimensional ultrasonic wave in the three-dimensional region in a living body. In an ultrasonic image display device for forming an image or a tomographic image, local area data for extracting echo data contained in a local area spread three-dimensionally around the target echo data to be filtered from the three-dimensional data memory Extraction means and individual median value determination means for determining the median value from the echo data in the reference direction for each of the n (n ≧ 3) reference directions that are set to intersect with each other and that pass through the echo data of interest. A representative median value determining means for further determining a representative median value from among the determined n median values, and Ultrasonic image display apparatus characterized by the updated value of the echo data.
JP3048301A 1991-03-13 1991-03-13 Ultrasonic image display device Expired - Fee Related JPH0675576B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP3048301A JPH0675576B2 (en) 1991-03-13 1991-03-13 Ultrasonic image display device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP3048301A JPH0675576B2 (en) 1991-03-13 1991-03-13 Ultrasonic image display device

Publications (2)

Publication Number Publication Date
JPH04285545A JPH04285545A (en) 1992-10-09
JPH0675576B2 true JPH0675576B2 (en) 1994-09-28

Family

ID=12799614

Family Applications (1)

Application Number Title Priority Date Filing Date
JP3048301A Expired - Fee Related JPH0675576B2 (en) 1991-03-13 1991-03-13 Ultrasonic image display device

Country Status (1)

Country Link
JP (1) JPH0675576B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007171051A (en) * 2005-12-22 2007-07-05 Honda Electronic Co Ltd Sonic velocity measuring method and device, and ultrasonic image inspecting device

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09179962A (en) * 1995-12-22 1997-07-11 Matsushita Electric Ind Co Ltd Picture processor
US6102858A (en) * 1998-04-23 2000-08-15 General Electric Company Method and apparatus for three-dimensional ultrasound imaging using contrast agents and harmonic echoes
JP5179801B2 (en) * 2007-08-24 2013-04-10 株式会社東芝 Ultrasonic image display method and apparatus

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007171051A (en) * 2005-12-22 2007-07-05 Honda Electronic Co Ltd Sonic velocity measuring method and device, and ultrasonic image inspecting device

Also Published As

Publication number Publication date
JPH04285545A (en) 1992-10-09

Similar Documents

Publication Publication Date Title
US4831528A (en) Apparatus and method for improvement of 3D images derived from tomographic data
JP2845995B2 (en) Region extraction method
EP1757955B1 (en) Apparatus and method for processing an ultrasound image
EP1791086A1 (en) Multi-resolution adaptive filtering
US20040165766A1 (en) Method and apparatus for forming and displaying projection image from a plurality of sectional images
US20040081340A1 (en) Image processing apparatus and ultrasound diagnosis apparatus
JPH1128211A (en) Three-dimensional imaging system and method for ultrasonic scattering medium
JP2000279416A (en) Three-dimensional imaging method and system
JPH03119489A (en) Method and apparatus for generating three- dimensional image from plural two- dimensional images
JP2002045358A (en) Automatic detection of cystic substance, and method and apparatus for dimensional measurement
DE19611990A1 (en) Method and device for producing large, composite ultrasound images
JP2007149092A (en) Multiple resolution adaptive filtering
JPS63125241A (en) Image processor
JPH08173414A (en) Method and equipment to decrease artifact in tomography picture
US8047991B2 (en) Automatic identification of orientation in medical diagnostic ultrasound
EP0318293B1 (en) Apparatus and method for processing tomographic data
CN104224230B (en) Three-dimensional and four-dimensional ultrasonic imaging method and device based on GPU (Graphics Processing Unit) platform and system
JPH0675576B2 (en) Ultrasonic image display device
KR20060074908A (en) Four-dimensional labeling apparatus, n-dimensional labeling apparatus, four-dimensional spatial filter apparatus, and n-dimensional spatial filter apparatus
US20100054567A1 (en) Method and apparatus for interactive ct reconstruction
CN116958481A (en) Point cloud reconstruction method and device, electronic equipment and readable storage medium
JP2004267584A (en) Ultrasonic diagnostic equipment
Kwon et al. A fast 3D adaptive bilateral filter for ultrasound volume visualization
JP2768627B2 (en) Ultrasound image processing device
CN109662689A (en) A kind of health early warning system of the hospital based on electrocardiogram

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090928

Year of fee payment: 15

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100928

Year of fee payment: 16

LAPS Cancellation because of no payment of annual fees