JPH02268749A - Image processing apparatus - Google Patents
Image processing apparatusInfo
- Publication number
- JPH02268749A JPH02268749A JP9140089A JP9140089A JPH02268749A JP H02268749 A JPH02268749 A JP H02268749A JP 9140089 A JP9140089 A JP 9140089A JP 9140089 A JP9140089 A JP 9140089A JP H02268749 A JPH02268749 A JP H02268749A
- Authority
- JP
- Japan
- Prior art keywords
- sound ray
- data value
- spatial filter
- smoothing
- unidimensional
- 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.)
- Pending
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Description
【発明の詳細な説明】
[産業上の利用分野]
この発明は、画像の平滑化等の演算を行なう画像処理装
置に関し、特に、超音波画像診断装置に有用である。DETAILED DESCRIPTION OF THE INVENTION [Industrial Field of Application] The present invention relates to an image processing device that performs calculations such as image smoothing, and is particularly useful for ultrasonic diagnostic imaging devices.
[従来の技術]
従来の画像処理装置を組み込んだ超音波血流イメージン
グ装置の一例を第5図に示す。[Prior Art] FIG. 5 shows an example of an ultrasonic blood flow imaging device incorporating a conventional image processing device.
この超音波血流イメージング装置51では、プローブ2
および送受信部3で音線(第6図に示すSl、S2.S
3)ごとに超音波エコー信号を得、B/Mモード処理部
4およびDSC5を経て2次元の画像とし、カラー表示
部6で87Mモードの画像表示が行なわれる。また、直
交検波部7で血流ドプラ信号を得、ドプラスペクトラム
処理部8を経て、ドプラスペクトラム表示が行なわれる
。In this ultrasonic blood flow imaging device 51, the probe 2
and sound rays (Sl, S2.S shown in FIG. 6) in the transmitting and receiving unit 3.
3), an ultrasonic echo signal is obtained, passed through the B/M mode processing section 4 and the DSC 5, and converted into a two-dimensional image, which is then displayed on the color display section 6 in 87M mode. Further, a blood flow Doppler signal is obtained by the orthogonal detection section 7, and a Doppler spectrum is displayed through the Doppler spectrum processing section 8.
また、MTIフィルタ9で固定信号成分を除去し、演算
部10で各観測点の流速値1分散値、パワーを算出し、
平滑部61で空間フィルタにより平滑化し、比較部12
で所定の閾値Th以上の領域のみを取り出し、DSC5
を経て、カラー表示部6で血流速断層像等を表示する。In addition, the MTI filter 9 removes the fixed signal component, and the calculation unit 10 calculates the flow velocity value 1 variance value and power at each observation point,
Smoothing is performed by a spatial filter in a smoothing unit 61, and the comparison unit 12
extracts only the area above a predetermined threshold Th, and performs DSC5
After that, a blood flow velocity tomographic image and the like are displayed on the color display section 6.
ここで平滑部61で行なわれる演算を例示すると、第6
図に示す原画像の3×3の2次元領域についての原デー
タ値P1.〜Pjmがあるとき、第7図に示す平滑化画
像の平滑化データ値R1!は、で算出される。上記[α
□]のマトリクスが3×3の2次元領域についての空間
フィルタである。Here, to illustrate the calculations performed by the smoothing section 61, the sixth
Original data value P1 for the 3x3 two-dimensional area of the original image shown in the figure. ~Pjm, the smoothed data value R1 of the smoothed image shown in FIG. 7! is calculated by Above [α
□] matrix is a spatial filter for a 3×3 two-dimensional area.
上記の平滑部61では、3×3の2次元領域についての
空間フィルタを用いており、1つの平滑化データ値を得
るためには、9回の乗算と8回の加算を必要とする。The above smoothing unit 61 uses a spatial filter for a 3×3 two-dimensional area, and requires nine multiplications and eight additions to obtain one smoothed data value.
一般に、(2m+1)X (2m+1)の2次元領域に
ついての空間フィルタを用いて平滑化や尖鋭化等の画像
処理を行なうと、1つの平滑化データ値を得るために、
(2m+1)”回の乗算と4m(m+1)回の加算とが
必要となる。Generally, when performing image processing such as smoothing and sharpening using a spatial filter for a two-dimensional area of (2m+1) x (2m+1), in order to obtain one smoothed data value,
(2m+1)" multiplications and 4m(m+1) additions are required.
すなわち、従来の空間フィルタを用いた画像処理装置で
は、空間フィルタを大きくしたり画像の面積を広くする
と演算量が極めて多くなるため、ハードウェア、ソフト
ウェアの負担が大きくなり、処理時間も長くかかる問題
点がある。In other words, with conventional image processing devices that use spatial filters, increasing the size of the spatial filter or widening the image area requires an extremely large amount of calculations, which increases the burden on hardware and software and increases processing time. There is a point.
そこで、この発明の目的は、演算量を低減できるように
した画像処理装置を提供することにある。SUMMARY OF THE INVENTION Therefore, an object of the present invention is to provide an image processing device that can reduce the amount of calculations.
[課題を解決するための手段]
この発明の画像処理装置は、2次元領域についての空間
フィルタ[αlk]に基づいて画像の平滑化等を行なう
画像処理装置において、縦または横方向の1次元領域に
ついての空間フィルタ[βb]を用い原画像に対して演
算を行なう第1演算手段と、前記1次元領域と交差する
方向の1次元領域についての空間フィルタ[γ、コを用
い前記第1演算手段の演算結果の画像に対して演算を行
なう第2演算手段とを具備したことを構成上の特徴とす
るものである。[Means for Solving the Problems] An image processing device of the present invention is an image processing device that performs image smoothing based on a spatial filter [αlk] for a two-dimensional region. a first calculation means for performing calculations on the original image using a spatial filter [βb] for , and a first calculation means for performing calculations on the original image using a spatial filter [γ, The configuration is characterized in that it includes a second calculation means for performing calculations on the image of the calculation result.
[作用]
第1演算手段でたとえば縦方向に空間フィルタ[β−で
演算し、次いで第2演算手段で横方向に空間フィルタ[
γc]で演算すると、結局は2次元的に空間フィルタ[
αlk]で演算することと等価となる。つまり、空間フ
ィルタ[βb]と空間フィルタ[γc]とを(2m+1
)の1次元領域についての空間フィルタとすれば、(2
m−1−1)X(2m+1)の2次元領域についての空
間フィルタを用いて画像処理を行なう場合と同じ結果を
得られる。[Operation] The first calculation means calculates the spatial filter [β-] in the vertical direction, and then the second calculation means calculates the spatial filter [β-] in the horizontal direction.
γc], the result is a two-dimensional spatial filter [
αlk]. In other words, the spatial filter [βb] and the spatial filter [γc] are (2m+1
), then (2
The same result as when image processing is performed using a spatial filter for a two-dimensional area of m-1-1) x (2m+1) can be obtained.
ところで、(2m+1)の1次元領域についての空間フ
ィルタを用いて画像処理を行なうと、1つの平滑化デー
タ値を得るのに(2m+1)回の乗算42m回の加算が
必要になる。従って、これを2回行なうときは、2 (
2m+1)回の乗算と4m回の加算とが必要になる。By the way, when image processing is performed using a spatial filter for a one-dimensional area of (2m+1), it takes (2m+1) multiplications and 42m additions to obtain one smoothed data value. Therefore, when doing this twice, 2 (
2m+1) multiplications and 4m additions are required.
これは従来の演算量(先述)よりも少ない演算量である
。This is a smaller amount of calculation than the conventional calculation amount (described earlier).
従って、ハードウェア、ソフトウェアの負担が軽減され
、処理時間も短縮できることとなる。Therefore, the burden on hardware and software is reduced, and processing time can also be shortened.
[実施例]
以下、図に示す実施例によりこの発明を更に詳細に説明
する。なお、これによりこの発明が限定されるものでは
ない。[Example] Hereinafter, the present invention will be explained in more detail with reference to Examples shown in the drawings. Note that this invention is not limited to this.
第1図はこの発明の一実施例の画像処理装置である平滑
部11を組み込んだ超音波血流イメージング装置1を示
すものである。FIG. 1 shows an ultrasonic blood flow imaging apparatus 1 incorporating a smoothing section 11, which is an image processing apparatus according to an embodiment of the present invention.
この超音波血流イメージング装置1において、プローブ
2.送受信部3.B/Mモード処理部4゜DSC5,カ
ラー表示部6.直交検波部7.ドプラスペクトラム処理
部8.MTIフィルタ9.演算部10および比較部12
は、第5図で示した装置51と同様の構成要素であり、
同じ参照番号を付している。In this ultrasonic blood flow imaging apparatus 1, a probe 2. Transmitting/receiving section 3. B/M mode processing section 4゜DSC5, color display section 6. Quadrature detection section 7. Doppler spectrum processing section 8. MTI filter9. Arithmetic unit 10 and comparison unit 12
are the same components as the device 51 shown in FIG.
They bear the same reference numbers.
平滑部11は、音線上平滑部11sと音線開平滑部li
dとからなっている。The smooth part 11 includes a sound ray smooth part 11s and a sound ray open smooth part li.
It consists of d.
音線上平滑部11sは、3×1の一次元領域についての
空間フィルタ[βb]を用いて、一つの音線上の観測点
に対して演算を行ない、中間平滑化データ値QIkを算
出する。The acoustic ray smoothing unit 11s calculates an intermediate smoothed data value QIk by using a spatial filter [βb] for a 3×1 one-dimensional area to perform calculations on observation points on one acoustic ray.
この演算を例示すると、第2図に示す各音線S1、S2
.S3上の原データ値pH〜P81+ P II〜P
*** P 18〜psaがあるとき、第3図に示
す中間平滑化データ値Q宮l。To illustrate this calculation, each sound ray S1, S2 shown in FIG.
.. Original data value on S3 pH~P81+ P II~P
*** When there is P 18~psa, the intermediate smoothed data value Q as shown in FIG.
Q1!+ Q10・ は、 で算出される。Q1! + Q10・ teeth, It is calculated by
音線開平滑部lidは、IX3の一次元領域についての
空間フィルタ[γc]を用いて、各音線の同−深さの中
間平滑化データ値に対して演算を行ない、平滑化データ
値R,kを算出する。The sound ray open smoothing unit lid performs calculations on the intermediate smoothed data value of the same depth of each sound ray using the spatial filter [γc] for the one-dimensional region of IX3, and calculates the smoothed data value R. , k.
この演算を例示すると、第3図に示す各音線S1、S2
.S3の同−深さの中間平滑化データ値Q□+ Q**
+ Q□があるとき、第4図に示す平滑化データ値R■
は、
で算出される。To illustrate this calculation, each sound ray S1, S2 shown in FIG.
.. Intermediate smoothed data value of S3 at the same depth Q□+ Q**
+ Q□, the smoothed data value R■ shown in Figure 4
is calculated as .
βb・γ、;α、1.βb・γ1=α28.・・・とす
れば、従来の2次元領域についての空間フィルタ[α+
hlを用いたのと同じ結果が得られる。すなわち、同じ
平滑化画像を得ることができる。βb・γ, ;α, 1. βb・γ1=α28. ..., then the conventional spatial filter [α+
The same result is obtained using hl. That is, the same smoothed image can be obtained.
R□を得るために必要な演算は、見かけ上は12回の乗
算と8回の加算になるが、R0以外の平滑化データ値の
算出にも中間平滑化データ値Q、1゜Q□、Q3.が使
われることを考慮すれば、実質的には6回の乗算と4回
の加算になる。すなわち、従来の9回の乗算と8回の加
算(先述)に比較して演算量が低減されている。The calculations required to obtain R□ appear to be 12 multiplications and 8 additions, but calculations of smoothed data values other than R0 also require intermediate smoothed data values Q, 1°Q□, Q3. Considering that , is used, it is actually 6 multiplications and 4 additions. That is, the amount of calculation is reduced compared to the conventional nine multiplications and eight additions (described above).
なお、原データ値P、には各音線ごとに順に得られるの
で、上記のように第1に音線上の原データ値から中間平
滑化データ値を算出し、第2に音線間の中間平滑化デー
タ値から平滑化データ値を算出する構成とするのが処理
のタイミングの観点から好ましい。しかし、逆に第1に
音線間の原データ値から中間平滑化データ値を算出し、
第2に音線上の中間平滑化データ値から平滑化データ値
を算出する構成としても演算量を低減させる効果は同じ
である。Note that the original data value P is obtained sequentially for each sound ray, so as described above, firstly, the intermediate smoothed data value is calculated from the original data value on the sound ray, and secondly, the intermediate smoothed data value between the sound rays is calculated from the original data value on the sound ray. From the viewpoint of processing timing, it is preferable to have a configuration in which the smoothed data value is calculated from the smoothed data value. However, conversely, first, intermediate smoothed data values are calculated from the original data values between sound rays,
Second, the effect of reducing the amount of calculation can be obtained even if the smoothed data value is calculated from the intermediate smoothed data value on the sound ray.
[発明の効果]
この発明の画像処理装置によれば、従来より少ない演算
量で従来と同じ空間フィルタを用いたときと同じ画像処
理を行なうことができる。[Effects of the Invention] According to the image processing device of the present invention, it is possible to perform the same image processing as when using the same spatial filter as in the past with a smaller amount of calculation than in the past.
従って、ハードウェアやソフトウェアの負担が軽減され
、処理時間も短縮されるので高速処理が可能となる。Therefore, the burden on hardware and software is reduced, and processing time is also shortened, making high-speed processing possible.
第1図はこの発明の一実施例の画像処理装置を組み込ん
だ超音波血流イメージング装置の要部ブロック図、第2
図は音線上の原データ値の例示概念図、第3図は音線上
の中間平滑化データ値の例示概念図、第4図は音線上の
平滑化データ値の例示概念図、第5図は従来の画像処理
装置を組み込んだ超音波血流イメージング装置の一例の
要部ブロック図、第6図は音線上の原データ値の例示概
念図、第7図は音線上の平滑化データ値の例示概念図で
ある。
(符号の説明)
1・・・・・・超音波血流イメージング装置2・・・・
・・プローブ
10・・・・・・演算部
11・・・・・・平滑部
11s・・・音線上平滑部
lid・・・音線開平滑部
pH〜P口・・・原データ値
Q+I−Qs*・・・中間平滑化データ値Rs I=
Rs s・・・平滑化データ値。
出願人 横河メディカルシステム株式会社第
図
l
D!
第
図
l
lFIG. 1 is a block diagram of main parts of an ultrasonic blood flow imaging device incorporating an image processing device according to an embodiment of the present invention, and FIG.
The figure is a conceptual diagram of an example of original data values on a sound ray, FIG. 3 is a conceptual diagram of an example of intermediate smoothed data values on a sound ray, FIG. 4 is a conceptual diagram of an example of smoothed data values on a sound ray, and FIG. A block diagram of essential parts of an example of an ultrasonic blood flow imaging device incorporating a conventional image processing device, FIG. 6 is a conceptual diagram illustrating original data values on a sound ray, and FIG. 7 is an example of smoothed data values on a sound ray. It is a conceptual diagram. (Explanation of symbols) 1... Ultrasonic blood flow imaging device 2...
... Probe 10 ... Calculation section 11 ... Smoothing section 11s ... Sound ray smoothing section lid ... Sound ray open smoothing section pH~P mouth ... Original data value Q+I- Qs*...Intermediate smoothed data value Rs I=
Rs s...Smoothed data value. Applicant Yokogawa Medical System Co., Ltd. Figure l D! Figure l l
Claims (1)
に基づいて画像の平滑化等を行なう画像処理装置におい
て、 縦または横方向の1次元領域についての空間フィルタ[
β_b]を用い原画像に対して演算を行なう第1演算手
段と、前記1次元領域と交差する方向の1次元領域につ
いての空間フィルタ[γ_c]を用い前記第1演算手段
の演算結果の画像に対して演算を行なう第2演算手段と
を具備したことを特徴とする画像処理装置。[Claims] Spatial filter [α_j_k] for one- and two-dimensional regions
In an image processing device that performs image smoothing etc. based on the spatial filter [
β_b] is used to perform calculations on the original image, and a spatial filter [γ_c] for a one-dimensional area in a direction intersecting the one-dimensional area is used to calculate the image of the calculation result of the first calculation means. An image processing device comprising: a second calculation means for performing calculations on the image processing device.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP9140089A JPH02268749A (en) | 1989-04-11 | 1989-04-11 | Image processing apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP9140089A JPH02268749A (en) | 1989-04-11 | 1989-04-11 | Image processing apparatus |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH02268749A true JPH02268749A (en) | 1990-11-02 |
Family
ID=14025331
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP9140089A Pending JPH02268749A (en) | 1989-04-11 | 1989-04-11 | Image processing apparatus |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH02268749A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004033732A (en) * | 2002-03-19 | 2004-02-05 | Zonare Medical Systems Inc | System and method for post-processing of ultrasonic color doppler imaging |
| JP2006255322A (en) * | 2005-03-18 | 2006-09-28 | Matsushita Electric Ind Co Ltd | Ultrasonic Doppler blood flow meter |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS54136135A (en) * | 1978-04-13 | 1979-10-23 | Iryo Gijutsu Kenkyu Kaihatsu Zaidan | Picture information processor |
| JPS6238981A (en) * | 1985-08-15 | 1987-02-19 | Canon Inc | Space filter |
-
1989
- 1989-04-11 JP JP9140089A patent/JPH02268749A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS54136135A (en) * | 1978-04-13 | 1979-10-23 | Iryo Gijutsu Kenkyu Kaihatsu Zaidan | Picture information processor |
| JPS6238981A (en) * | 1985-08-15 | 1987-02-19 | Canon Inc | Space filter |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004033732A (en) * | 2002-03-19 | 2004-02-05 | Zonare Medical Systems Inc | System and method for post-processing of ultrasonic color doppler imaging |
| JP2006255322A (en) * | 2005-03-18 | 2006-09-28 | Matsushita Electric Ind Co Ltd | Ultrasonic Doppler blood flow meter |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2934402B2 (en) | Three-dimensional ultrasonic image creating method and image processing apparatus | |
| JP5048175B2 (en) | An ultrasound system that executes computer instructions in parallel | |
| Nelson et al. | Visualization of 3D ultrasound data | |
| McCann et al. | Multidimensional ultrasonic imaging for cardiology | |
| JP7041147B2 (en) | Systems and methods that characterize hepatic perfusion of contrast medium flow | |
| US5779641A (en) | Method and apparatus for three-dimensional ultrasound imaging by projecting filtered pixel data | |
| US6482161B1 (en) | Medical diagnostic ultrasound system and method for vessel structure analysis | |
| TW447215B (en) | Ultrasound scan conversion with spatial dithering | |
| JP4204095B2 (en) | 3D imaging system and method for subject volume | |
| KR870000638B1 (en) | Ultrasonic Fluid Observation Device | |
| US20040081340A1 (en) | Image processing apparatus and ultrasound diagnosis apparatus | |
| CN103156637A (en) | Ultrasound size image data processing method and equipment | |
| Lovstakken et al. | Blood flow imaging-a new real-time, flow imaging technique | |
| JPH0428354A (en) | Ultrasonic diagnosing device | |
| Christofides et al. | Automatic respiratory gating for contrast ultrasound evaluation of liver lesions | |
| JPH02268749A (en) | Image processing apparatus | |
| CN114680940A (en) | Ultrasonic image-based vascular plaque presenting method and ultrasonic imaging system | |
| US20230329670A1 (en) | Ultrasonic measurement of vessel stenosis | |
| JP2001252270A (en) | Ultrasonic imaging method and ultrasound diagnostic apparatus | |
| US7118532B2 (en) | Data dependent color wall filters | |
| JPH04348745A (en) | Ultrasonic image processor | |
| JP3039808B2 (en) | Ultrasound blood flow drawing device | |
| JPH07328001A (en) | Ultrasonic bloodstream image generating method and ultrasonic diagnostic device | |
| Gee et al. | Sequential 3d diagnostic ultrasound using the stradx system | |
| CN109598718B (en) | Equipment replacement demand analysis agency |