JP2011056160A - Ultrasonic diagnostic apparatus - Google Patents

Ultrasonic diagnostic apparatus Download PDF

Info

Publication number
JP2011056160A
JP2011056160A JP2009211577A JP2009211577A JP2011056160A JP 2011056160 A JP2011056160 A JP 2011056160A JP 2009211577 A JP2009211577 A JP 2009211577A JP 2009211577 A JP2009211577 A JP 2009211577A JP 2011056160 A JP2011056160 A JP 2011056160A
Authority
JP
Japan
Prior art keywords
image
diagnostic apparatus
ultrasonic diagnostic
images
virtual
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
JP2009211577A
Other languages
Japanese (ja)
Other versions
JP5475376B2 (en
Inventor
Eiji Kasahara
英司 笠原
Masaru Murashita
賢 村下
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 JP2009211577A priority Critical patent/JP5475376B2/en
Priority to US12/755,045 priority patent/US8852106B2/en
Priority to EP10003745.6A priority patent/EP2241256B1/en
Priority to CN201010148456.4A priority patent/CN101897601B/en
Publication of JP2011056160A publication Critical patent/JP2011056160A/en
Application granted granted Critical
Publication of JP5475376B2 publication Critical patent/JP5475376B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

<P>PROBLEM TO BE SOLVED: To improve reliability of a display image of an object tissue in unstable periodic motion. <P>SOLUTION: A representative part memory 15 stores echo data collected from a representative part of an object tissue. A virtual period setting unit 22 sets a virtual period as a temporary period of the motion of the object tissue on the basis of the echo data of a plurality of time phase stored in the representative part memory 15. A mutual difference value is used to calculate the virtual period. A pre-memory 14 stores a plurality of tomographic image data in time-series order. A reference image searching section 24 extracts a plurality of reference images from a plurality of tomographic images using the set virtual period. A reconstruction processing unit 20 divides the plurality of tomographic image data into a plurality of image groups by using each of the plurality of reference images as a division unit. A three-dimensional image is reconstructed by extracting, from the plurality of image groups, a plurality of sets of tomographic image data that correspond periodically to each other. <P>COPYRIGHT: (C)2011,JPO&amp;INPIT

Description

本発明は、周期的に運動する対象組織の表示画像を形成する超音波診断装置に関する。   The present invention relates to an ultrasonic diagnostic apparatus that forms a display image of a target tissue that moves periodically.

心臓などの運動を伴う組織の三次元超音波画像を形成する超音波診断装置が知られている。例えば、三次元空間内において超音波ビームをスキャン(走査)して三次元空間内からエコーデータを収集し、収集したエコーデータに基づいて三次元超音波画像を形成してリアルタイム表示する技術が知られている。但し、リアルタイム表示の場合には、スキャンレートとビーム密度とビーム範囲が互いにトレードオフの関係になるという原理的な制約がある。   2. Description of the Related Art An ultrasonic diagnostic apparatus that forms a three-dimensional ultrasonic image of a tissue that accompanies motion such as the heart is known. For example, there is a technology that scans an ultrasonic beam in a three-dimensional space, collects echo data from the three-dimensional space, forms a three-dimensional ultrasonic image based on the collected echo data, and displays it in real time. It has been. However, in the case of real-time display, there is a principle restriction that the scan rate, the beam density, and the beam range are in a trade-off relationship with each other.

三次元超音波画像のリアルタイム表示における原理的な制約を回避するための技術も提案されている。例えば、特許文献1には、心電信号などに同期させて三次元空間内において走査面を少しずつ移動させながら、走査面の各位置において複数の時相に亘って複数の断層画像データを収集し、収集された複数の断層画像データを並べ替えて再構築して三次元画像データを形成する技術(再構築処理)が記載されている。この技術は、直接的に心電信号を得ることが困難な胎児などに適用することが難しい。   Techniques for avoiding the fundamental limitations in real-time display of 3D ultrasound images have also been proposed. For example, Patent Document 1 collects a plurality of tomographic image data over a plurality of time phases at each position on the scanning plane while moving the scanning plane little by little in the three-dimensional space in synchronization with an electrocardiogram signal or the like. A technique (reconstruction processing) is described in which a plurality of collected tomographic image data is rearranged and reconstructed to form three-dimensional image data. This technique is difficult to apply to a fetus or the like for which it is difficult to directly obtain an electrocardiogram signal.

また、特許文献2には、心電信号に換えて、ある時間間隔ごとにスキャンして再構築する技術が記載されている。しかし、この技術では、データ収集中における心臓などの周期が一定と仮定しており、そのため、例えば心臓の周期が一定ではない場合に、再構築後の画像における心臓の形態が実際のものから歪められて信頼性が低くなる可能性がある。   Patent Document 2 describes a technique for scanning and reconstructing at certain time intervals instead of an electrocardiographic signal. However, this technique assumes that the period of the heart, etc. during data collection is constant, so that, for example, if the period of the heart is not constant, the shape of the heart in the reconstructed image is distorted from the actual one. Reliability may be reduced.

特許第3537594号公報Japanese Patent No. 3537594 特開2005−74225号公報JP 2005-74225 A

上述した背景技術に鑑み、本願の発明者は、再構築処理により超音波画像を形成する技術について研究開発を重ねてきた。特に、胎児の心臓などのように運動の周期が不安定な対象組織についての超音波画像を形成する技術に注目した。   In view of the background art described above, the inventor of the present application has repeatedly researched and developed a technique for forming an ultrasonic image by reconstruction processing. In particular, we focused on techniques for forming ultrasound images of target tissues with unstable movement cycles, such as the fetal heart.

本発明は、その研究開発の過程において成されたものであり、その目的は、運動の周期が不安定な対象組織についての表示画像の信頼性を高めることにある。   The present invention has been made in the course of research and development, and an object thereof is to increase the reliability of a display image of a target tissue having an unstable motion cycle.

上記目的にかなう好適な超音波診断装置は、周期的に運動する対象組織を含む三次元空間内で超音波を送受するプローブと、プローブを制御することにより、前記運動の複数の周期に亘って対象組織の代表部位から受信信号を収集し、さらに、前記運動の複数の周期に亘って走査面を移動させつつ前記三次元空間内で複数の走査面を形成する送受信制御部と、前記代表部位から収集された受信信号に基づいて対象組織に関する運動の仮想周期を設定する仮想周期設定部と、前記複数の走査面に対応した複数の画像の中から前記仮想周期に対応した間隔で複数の基準画像を探索する基準画像探索部と、前記複数の基準画像の各々を分割の単位とすることにより前記複数の画像を複数の画像群に分割し、複数の画像群の各々から互いに周期的に対応した複数の画像を抽出する画像再構成部と、互いに周期的に対応した複数の画像に基づいて対象組織の表示画像を形成する表示画像形成部と、を有することを特徴とする。   An ultrasonic diagnostic apparatus suitable for the above-described object includes a probe that transmits and receives ultrasonic waves in a three-dimensional space including a target tissue that periodically moves, and a plurality of cycles of the movement by controlling the probe. A transmission / reception control unit that collects reception signals from a representative part of a target tissue and further forms a plurality of scanning planes in the three-dimensional space while moving the scanning plane over a plurality of cycles of the movement, and the representative part A virtual period setting unit for setting a virtual period of motion related to the target tissue based on the received signal collected from the plurality of references at intervals corresponding to the virtual period from among a plurality of images corresponding to the plurality of scanning planes A reference image search unit for searching for an image, and dividing each of the plurality of reference images into a plurality of image groups by using each of the plurality of reference images as a unit of division, and periodically matching each other from each of the plurality of image groups. An image reconstruction unit for extracting a plurality of images, characterized in that and a display image forming unit for forming a display image of the target tissue based on a plurality of images periodically correspond to one another.

望ましい具体例として、前記仮想周期設定部は、前記代表部位から収集された複数時相の受信信号に基づいて、時相間における受信信号の変化を反映させた特徴量を算出し、当該特徴量に基づいて前記仮想周期を決定する、ことを特徴とする。   As a desirable specific example, the virtual cycle setting unit calculates a feature amount that reflects a change in the received signal between time phases based on the received signals of a plurality of time phases collected from the representative part, and calculates the feature amount. The virtual period is determined on the basis of this.

望ましい具体例として、前記仮想周期設定部は、前記特徴量として各時相ごとに相互差分値を算出し、複数時相に亘って得られた相互差分値の時間的変化に基づいて前記仮想周期を決定する、ことを特徴とする。   As a desirable specific example, the virtual cycle setting unit calculates a mutual difference value for each time phase as the feature amount, and the virtual cycle based on temporal changes of the mutual difference values obtained over a plurality of time phases. It is characterized by determining.

望ましい具体例として、前記仮想周期設定部は、複数時相に亘って得られた相互差分値の極大値に対応した時相の間隔に基づいて前記仮想周期を決定する、ことを特徴とする。   As a desirable specific example, the virtual cycle setting unit is characterized in that the virtual cycle is determined based on a time phase interval corresponding to a maximum value of a mutual difference value obtained over a plurality of time phases.

望ましい具体例として、前記仮想周期設定部は、前記代表部位に設定された超音波ビームから1次元的に収集された複数時相の受信信号を利用する、ことを特徴とする。   As a desirable specific example, the virtual cycle setting unit uses reception signals of a plurality of time phases collected one-dimensionally from an ultrasonic beam set in the representative part.

望ましい具体例として、前記仮想周期設定部は、前記代表部位に設定された走査面から2次元的に収集された複数時相の受信信号を利用する、ことを特徴とする。   As a desirable specific example, the virtual cycle setting unit uses reception signals of a plurality of time phases collected two-dimensionally from the scanning plane set in the representative part.

望ましい具体例として、前記送受信制御部は、前記代表部位から受信信号を収集した後に、前記三次元空間内で複数の走査面を形成する、ことを特徴とする。   As a preferred specific example, the transmission / reception control unit forms a plurality of scanning planes in the three-dimensional space after collecting reception signals from the representative part.

望ましい具体例として、前記送受信制御部は、前記三次元空間内で複数の走査面を形成した後に、前記代表部位から受信信号を収集する、ことを特徴とする。   As a desirable specific example, the transmission / reception control unit collects reception signals from the representative portion after forming a plurality of scanning planes in the three-dimensional space.

望ましい具体例として、前記基準画像探索部は、前記複数の画像の各々に対応した相互差分値に基づいて、最大の相互差分値に対応した一つの画像を代表基準画像とし、その代表基準画像を起点として、極大の相互差分値に対応した複数の画像の中から、前記仮想周期だけ離れた位置に最も近い画像を次々に探索し、前記複数の基準画像とする、ことを特徴とする。   As a preferred specific example, the reference image search unit sets one image corresponding to the maximum mutual difference value as a representative reference image based on the mutual difference value corresponding to each of the plurality of images, and the representative reference image is used as the representative reference image. As a starting point, an image closest to a position separated by the virtual period is searched one after another from a plurality of images corresponding to the maximum mutual difference value, and is set as the plurality of reference images.

本発明により、運動の周期が不安定な対象組織についての表示画像の信頼性を高めることが可能になる。   According to the present invention, it is possible to improve the reliability of a display image for a target tissue whose movement cycle is unstable.

本発明に係る超音波診断装置の好適な実施形態を示すブロック図である。1 is a block diagram showing a preferred embodiment of an ultrasonic diagnostic apparatus according to the present invention. 本実施形態における三次元的な走査を説明するための図である。It is a figure for demonstrating the three-dimensional scanning in this embodiment. 断面差分値の変化を示す図である。It is a figure which shows the change of a cross-sectional difference value. 代表部位に関する相互差分値の変化を示す図である。It is a figure which shows the change of the mutual difference value regarding a representative site | part. 基準画像の探索を説明するための図である。It is a figure for demonstrating the search of a reference | standard image. 再構築処理部による処理を説明するための図である。It is a figure for demonstrating the process by the reconstruction process part. 再構築処理部による別の好適な処理を説明するための図である。It is a figure for demonstrating another suitable process by the reconstruction process part.

以下に本発明の好適な実施形態を説明する。   Hereinafter, preferred embodiments of the present invention will be described.

図1は、本発明の実施において好適な超音波診断装置の全体構成を示す図である。プローブ10は、対象組織を含む三次元空間内において超音波を送受する。プローブ10は、各々が超音波を送受する複数の振動素子を備えており、複数の振動素子がビームフォーマ12によって送信制御されて送信ビームが形成される。また、複数の振動素子が対象組織から反射された超音波を受波し、これにより得られた信号がビームフォーマ12へ出力され、ビームフォーマ12が受信ビームを形成する。   FIG. 1 is a diagram showing an overall configuration of an ultrasonic diagnostic apparatus suitable for implementing the present invention. The probe 10 transmits and receives ultrasonic waves in a three-dimensional space including the target tissue. The probe 10 includes a plurality of vibration elements each of which transmits / receives an ultrasonic wave, and transmission of the plurality of vibration elements is controlled by the beam former 12 to form a transmission beam. Further, the plurality of vibration elements receive the ultrasonic waves reflected from the target tissue, and signals obtained thereby are output to the beam former 12, and the beam former 12 forms a reception beam.

本実施形態のプローブ10は、超音波ビーム(送信ビームと受信ビーム)を三次元空間内において走査して立体的にエコーデータを収集する3Dプローブである。例えば、一次元的に配列された複数の振動素子(1Dアレイ振動子)によって電子的に形成される走査面を機械的に動かすことにより超音波ビームが三次元的に走査される。また、二次元的に配列された複数の振動素子(2Dアレイ振動子)を電子的に制御して超音波ビームを三次元的に走査してもよい。   The probe 10 of the present embodiment is a 3D probe that collects echo data three-dimensionally by scanning an ultrasonic beam (a transmission beam and a reception beam) in a three-dimensional space. For example, the ultrasonic beam is scanned three-dimensionally by mechanically moving a scanning surface formed electronically by a plurality of vibration elements (1D array transducers) arranged one-dimensionally. Alternatively, the ultrasonic beam may be scanned three-dimensionally by electronically controlling a plurality of vibration elements (2D array transducers) arranged two-dimensionally.

ビームフォーマ12は、プローブ10が備える複数の振動素子の各々に対応した送信信号を供給することにより超音波の送信ビームを形成する。また、ビームフォーマ12は、プローブ10が備える複数の振動素子の各々から得られる受信信号に対して整相加算処理などを施すことにより超音波の受信ビームを形成し、受信ビームに沿って得られるエコーデータを出力する。本実施形態において、対象組織は、周期的に運動する組織であり、例えば胎児の心臓などである。   The beam former 12 forms an ultrasonic transmission beam by supplying a transmission signal corresponding to each of a plurality of vibration elements included in the probe 10. In addition, the beam former 12 forms an ultrasonic reception beam by performing a phasing addition process or the like on a reception signal obtained from each of a plurality of vibration elements included in the probe 10, and is obtained along the reception beam. Output echo data. In the present embodiment, the target tissue is a periodically moving tissue, such as a fetal heart.

図2は、本実施形態における三次元的な走査を説明するための図である。図2において対象組織を含む三次元空間はXYZ直交座標系で表現されている。本実施形態では、XY平面に対してほぼ平行となるように走査面Sが形成され、その走査面SをZ軸方向にゆっくりと移動させつつ、Z軸方向に沿って複数の走査面Sが形成される。走査面Sは、胎児の心臓などの周期的な運動に関する複数の周期に亘って、例えば約8秒で約20心拍を含む期間に亘って、Z軸方向にゆっくりと移動する。   FIG. 2 is a diagram for explaining three-dimensional scanning in the present embodiment. In FIG. 2, the three-dimensional space including the target tissue is expressed in an XYZ orthogonal coordinate system. In the present embodiment, the scanning surface S is formed so as to be substantially parallel to the XY plane, and a plurality of scanning surfaces S are formed along the Z-axis direction while slowly moving the scanning surface S in the Z-axis direction. It is formed. The scanning plane S moves slowly in the Z-axis direction over a plurality of periods related to periodic movements such as the fetal heart, for example, over a period including about 20 heartbeats in about 8 seconds.

図1に戻り、胎児の心拍の複数の周期に亘ってZ軸方向に沿って複数の走査面が形成されると、各走査面ごとに断層画像データが収集され、複数の走査面に対応した複数の断層画像データが次々に前メモリ14に記憶される。   Returning to FIG. 1, when a plurality of scanning planes are formed along the Z-axis direction over a plurality of periods of the fetal heartbeat, tomographic image data is collected for each scanning plane and corresponds to the plurality of scanning planes. A plurality of tomographic image data is sequentially stored in the previous memory 14.

エラー判定部16は、前メモリ14に記憶された複数の断層画像データから得られる画像間の差分量に基づいて、複数の断層画像データが良好か否かを判定する。例えば胎児や母体やプローブの動きにより画像内で胎児の心臓が大きく動いてしまい、良好な画像が得られない可能性がある。そこで、エラー判定部16は、診断のための良好な画像が得られるか否かを判定する。その判定にあたって、エラー判定部16は、次式で定義される断面差分値を利用する。   The error determination unit 16 determines whether or not the plurality of tomographic image data is good based on the difference amount between the images obtained from the plurality of tomographic image data stored in the previous memory 14. For example, the fetal heart may move greatly in the image due to the movement of the fetus, mother or probe, and a good image may not be obtained. Therefore, the error determination unit 16 determines whether a good image for diagnosis is obtained. In the determination, the error determination unit 16 uses a cross-sectional difference value defined by the following equation.

Figure 2011056160
Figure 2011056160

数1式におけるx,y,zは、図2のXYZ直交座標系における各軸上の座標値でありpは断層画像データ内の各座標に対応した画素値である。数1式により、Z軸方向に隣接する2つの断層画像データ間の差分値が算出される。   X, y, and z in Equation 1 are coordinate values on each axis in the XYZ orthogonal coordinate system of FIG. 2, and p is a pixel value corresponding to each coordinate in the tomographic image data. The difference value between two pieces of tomographic image data adjacent in the Z-axis direction is calculated by Equation (1).

図3は、断面差分値の変化を示す図であり、図3の横軸は、各断層画像データの位置を示している。つまり、図3の横軸は、各走査面の位置と各走査面が得られた時間に対応しており、図2のZ軸(時間の経過に伴う位置の変化方向)に対応している。   FIG. 3 is a diagram showing changes in cross-sectional difference values, and the horizontal axis of FIG. 3 shows the position of each tomographic image data. That is, the horizontal axis in FIG. 3 corresponds to the position of each scanning plane and the time at which each scanning plane is obtained, and corresponds to the Z axis in FIG. 2 (the direction of change in position over time). .

胎児の心臓が大きく移動してしまうことが無ければ、隣接する断層画像データは互いに似たものとなり、数1式により得られる差分値は比較的小さくなる。一方、例えば胎児自身の動き、母体の呼吸動作、プローブの位置の大きなずれなどがあると、断層画像内において胎児の心臓が大きく動いてしまい、隣接する断層画像データ間の差分値が比較的大きくなる。そこで、エラー判定部16は、断面差分値が所定の閾値を超えた場合に、画像内において心臓が大きくずれてしまったと判断する。   If the fetal heart does not move significantly, the adjacent tomographic image data will be similar to each other, and the difference value obtained from Equation 1 will be relatively small. On the other hand, for example, if there is a movement of the fetus itself, the mother's breathing movement, a large displacement of the probe position, the fetal heart moves greatly in the tomographic image, and the difference value between adjacent tomographic image data is relatively large. Become. Therefore, the error determination unit 16 determines that the heart has greatly shifted in the image when the cross-sectional difference value exceeds a predetermined threshold value.

図1に戻り、エラー判定部16により心臓が大きくずれてしまったと判断されると、制御部40は、例えば、ビームフォーマ12などを制御して、断層画像データの収集を中止させる。なお、制御部40は、図1内の各部を集中的に制御しており、例えば、エラー判定部16によりエラーであると判断された場合に、エラーである旨を示す表示や警告などを表示部30に表示させてもよい。エラー判定部16によりエラーの判定が成されなければ、前メモリ14に記憶された複数の断層画像データに基づいて、後述する処理が実行される。   Returning to FIG. 1, when the error determination unit 16 determines that the heart has greatly shifted, the control unit 40 controls, for example, the beamformer 12 to stop collecting tomographic image data. Note that the control unit 40 controls each unit in FIG. 1 in a concentrated manner. For example, when the error determination unit 16 determines that an error has occurred, a display or a warning indicating an error is displayed. You may display on the part 30. If the error determination unit 16 does not determine an error, a process to be described later is executed based on a plurality of tomographic image data stored in the previous memory 14.

また、本実施形態においては、上述した三次元的なエコーデータ(受信信号)の収集に加えて、対象組織の運動に関する複数の周期に亘って、対象組織の代表部位からエコーデータが収集される。例えば、対象組織が胎児の心臓の場合には、胎児の心臓の中心付近において固定的に走査面が形成され、その走査面内から、複数の心拍を含む期間(例えば、1〜2秒間)に亘って、複数時相のエコーデータ(受信信号)が収集される。収集されたエコーデータは、代表部位メモリ15に記憶される。例えば、前メモリ14に記憶される三次元的なエコーデータが収集されてから、代表部位メモリ15に記憶される代表部位からのエコーデータが収集される。もちろん、代表部位メモリ15に記憶されるエコーデータが収集されてから、前メモリ14に記憶されるエコーデータが収集されてもよい。   In the present embodiment, in addition to the collection of the above-described three-dimensional echo data (reception signal), echo data is collected from the representative portion of the target tissue over a plurality of periods related to the motion of the target tissue. . For example, when the target tissue is a fetal heart, a scan plane is formed in the vicinity of the center of the fetal heart, and within a period (for example, 1 to 2 seconds) including a plurality of heartbeats from within the scan plane. Over a period of time, echo data (reception signals) of a plurality of time phases are collected. The collected echo data is stored in the representative part memory 15. For example, after the three-dimensional echo data stored in the previous memory 14 is collected, the echo data from the representative part stored in the representative part memory 15 is collected. Of course, the echo data stored in the previous memory 14 may be collected after the echo data stored in the representative part memory 15 is collected.

仮想周期設定部22は、代表部位メモリ15に記憶された複数時相のエコーデータに基づいて、胎児の心臓に関する仮の周期となる仮想周期を算出する。仮想周期の算出にあたって、仮想周期設定部22は、次式で定義される相互差分値を利用する。   The virtual cycle setting unit 22 calculates a virtual cycle serving as a temporary cycle related to the fetal heart based on the echo data of a plurality of time phases stored in the representative part memory 15. In calculating the virtual period, the virtual period setting unit 22 uses a mutual difference value defined by the following equation.

Figure 2011056160
Figure 2011056160

数2式おけるx,yは、図2のXY直交座標系における各軸上の座標値であり、pは、代表部位に固定的に設定された走査面内(断層画像内)の各座標に対応した画素値である。また、仮想周期設定部22による仮想周期の設定にあたって、数2式におけるzは時刻である。つまり、数2式において、隣接する2つの時刻に対応した2つの断層画像データ間の2つの画素値の差分に対して、一方の画素値が乗算されている。これにより、心臓が収縮する場合に比べて心臓が拡張する場合に相互差分値が比較的大きな値となり、単純な差分値では識別が難しい拡張と収縮を相互差分値により識別することが可能になる。   X and y in Equation 2 are coordinate values on each axis in the XY orthogonal coordinate system of FIG. 2, and p is a coordinate in the scanning plane (in the tomographic image) fixedly set to the representative part. The corresponding pixel value. In setting the virtual cycle by the virtual cycle setting unit 22, z in Equation 2 is time. That is, in Equation 2, one pixel value is multiplied by the difference between two pixel values between two tomographic image data corresponding to two adjacent times. As a result, the mutual difference value is relatively large when the heart expands compared to when the heart contracts, and it is possible to identify expansion and contraction that are difficult to identify with a simple difference value by the mutual difference value. .

例えば、ある断層画像データz内において、画素p(x,y,z)が心臓内壁の近傍の心筋であると仮定し、その画素値をp(x,y,z)=100とする。心臓が拡張して心腔が大きくなると、断層画像データzに続いて得られる断層画像データz+1内において画素p(x,y,z+1)が心腔の画素となる。心筋に比べて心腔の画素値は小さいためその画素値をp(x,y,z+1)=10とする。この例において、数2式の右辺の絶対値を算出すると100×(100−10)=9000となる。心臓が拡張する場合には、心臓内壁の周辺において、心筋から心腔に変化する画素が多く発生するため、数2式の相互差分値の値が比較的大きくなる。   For example, in a certain tomographic image data z, it is assumed that the pixel p (x, y, z) is a myocardium in the vicinity of the inner wall of the heart, and the pixel value is set to p (x, y, z) = 100. When the heart expands and the heart chamber becomes larger, the pixel p (x, y, z + 1) becomes the heart chamber pixel in the tomographic image data z + 1 obtained following the tomographic image data z. Since the pixel value of the heart chamber is smaller than that of the myocardium, the pixel value is set to p (x, y, z + 1) = 10. In this example, the absolute value of the right side of Equation 2 is calculated to be 100 × (100−10) = 9000. When the heart expands, there are many pixels that change from the myocardium to the heart chamber around the inner wall of the heart, so that the mutual difference value of Equation 2 becomes relatively large.

一方、心臓が収縮する場合には、上記の例とは反対の現象が発生する。つまり、心臓が収縮して心腔が小さくなるため、心腔に対応した画素p(x,y,z)=10から、心筋に対応した画素p(x,y,z+1)=100に変化する。この例において、数2式の右辺の絶対値を算出すると|10×(10−100)|=900となり、拡張の場合における値9000よりも小さくなる。そのため、拡張と収縮を相互差分値により識別することが可能になる。   On the other hand, when the heart contracts, a phenomenon opposite to the above example occurs. That is, since the heart contracts and the heart chamber becomes smaller, the pixel p (x, y, z) = 10 corresponding to the heart chamber changes from the pixel p (x, y, z + 1) = 100 corresponding to the heart muscle. . In this example, when the absolute value of the right side of Equation 2 is calculated, | 10 × (10−100) | = 900, which is smaller than the value 9000 in the case of expansion. Therefore, expansion and contraction can be identified by the mutual difference value.

図4は、代表部位に関する相互差分値の変化を示す図である。図4の横軸は、各断層画像データが得られる時刻(時間)を示している。数2式を利用して各時刻(z)において相互差分値が算出されると、心臓が拡張する場合に相互差分値が比較的大きな値となる。そこで、仮想周期設定部22は、相互差分値のピーク値(極大値)を検出し、隣接するピーク値の間隔を心臓の周期(心拍の周期)と判断する。   FIG. 4 is a diagram illustrating a change in the mutual difference value regarding the representative site. The horizontal axis in FIG. 4 indicates the time (time) at which each tomographic image data is obtained. When the mutual difference value is calculated at each time (z) using Equation 2, the mutual difference value becomes a relatively large value when the heart expands. Therefore, the virtual cycle setting unit 22 detects the peak value (maximum value) of the mutual difference value, and determines the interval between adjacent peak values as the heart cycle (heartbeat cycle).

但し、例えば胎児の心臓は、心拍の周期が変動する場合があり、心拍の周期が変動するとピーク値の間隔も変動する。そこで、仮想周期設定部22は、例えば、ピーク値の間隔のうちの代表的な(例えば2番目に大きな)間隔を仮想周期に設定する。なお、ピーク値の間隔の平均値を仮想周期としてもよいし、ピーク値の間隔のヒストグラムから得られる最多頻度の値や重心値などを仮想周期としてもよい。   However, for example, in the fetal heart, the heartbeat period may fluctuate, and when the heartbeat period fluctuates, the interval between peak values also fluctuates. Therefore, the virtual cycle setting unit 22 sets, for example, a representative (for example, second largest) interval among the peak value intervals as the virtual cycle. Note that an average value of peak value intervals may be a virtual cycle, or a most frequently used value or a centroid value obtained from a histogram of peak value intervals may be a virtual cycle.

また、仮想周期を設定するにあたって、Mモード計測が利用されてもよい。例えば、対象組織が胎児の心臓の場合には、胎児の心臓の中心付近において固定的に超音波ビームが設定され、その超音波ビームから、複数の心拍を含む期間に亘って複数時相のエコーデータ(ビームデータ)が収集される。この場合においても、収集されたエコーデータは、代表部位メモリ15に記憶される。ビームデータの場合には、数2式おけるxは、超音波ビームの位置に応じた固定値となる。そして、複数の時刻(z)に亘って、数2式に基づいて相互差分値を算出することにより、ビームデータの場合においても、図4と同様な周期的な波形が形成され、断層画像データの場合と同様に、ピーク値の間隔のうちの代表的な間隔やピーク値の間隔の平均値などが仮想周期とされる。   In setting the virtual period, M-mode measurement may be used. For example, when the target tissue is a fetal heart, an ultrasonic beam is fixedly set near the center of the fetal heart, and echoes of a plurality of time phases are generated from the ultrasonic beam over a period including a plurality of heartbeats. Data (beam data) is collected. Also in this case, the collected echo data is stored in the representative part memory 15. In the case of beam data, x in Equation 2 is a fixed value corresponding to the position of the ultrasonic beam. Then, by calculating the mutual difference value based on Equation 2 over a plurality of times (z), a periodic waveform similar to that in FIG. 4 is formed even in the case of beam data, and tomographic image data As in the case of, a typical interval among the peak value intervals, an average value of the peak value intervals, and the like are used as the virtual period.

仮想周期が設定されると、図1の基準画像探索部24は、複数の断層画像データの中から、仮想周期を利用して、複数の基準画像を探索する。この探索においても、数2式で定義される相互差分値が利用される。   When the virtual period is set, the reference image search unit 24 in FIG. 1 searches for a plurality of reference images from the plurality of tomographic image data using the virtual period. Also in this search, the mutual difference value defined by Equation 2 is used.

基準画像探索部24は、前メモリ14に記憶された複数の断層画像データに対して、数2式を適用して相互差分値を算出する。基準画像探索部24による相互差分値の算出において、数2式おけるx,y,zは、図2のXYZ直交座標系における各軸上の座標値であり、pは断層画像データ内の各座標に対応した画素値である。この場合、数2式において、Z軸方向に隣接する2つの断層画像データ間の2つの画素値の差分に対して、一方の画素値が乗算されている。これにより、心臓が収縮する場合に比べて心臓が拡張する場合に相互差分値が比較的大きな値となり、単純な差分値では識別が難しい拡張と収縮を相互差分値により識別することが可能になる。   The reference image search unit 24 calculates the mutual difference value by applying Equation 2 to the plurality of tomographic image data stored in the previous memory 14. In the calculation of the mutual difference value by the reference image search unit 24, x, y, z in Equation 2 are coordinate values on each axis in the XYZ orthogonal coordinate system of FIG. 2, and p is each coordinate in the tomographic image data. Is a pixel value corresponding to. In this case, in Equation 2, one pixel value is multiplied by a difference between two pixel values between two tomographic image data adjacent in the Z-axis direction. As a result, the mutual difference value is relatively large when the heart expands compared to when the heart contracts, and it is possible to identify expansion and contraction that are difficult to identify with a simple difference value by the mutual difference value. .

図5は、基準画像の探索を説明するための図である。図5(A)〜(C)の各々には、基準画像探索部24において算出された相互差分値の変化が図示されている。つまり、図5において、横軸は、各断層画像データの位置(各走査面の位置と時間)を示しており、図2のZ軸(時間の経過に伴う位置の変化方向)に対応している。数2式を利用してZ軸上の各位置(z)において相互差分値が算出されると、心臓が拡張する場合に相互差分値が比較的大きな値となる。   FIG. 5 is a diagram for explaining the search for the reference image. Each of FIGS. 5A to 5C illustrates a change in the mutual difference value calculated by the reference image search unit 24. That is, in FIG. 5, the horizontal axis indicates the position of each tomographic image data (position and time of each scanning plane), and corresponds to the Z axis of FIG. 2 (change direction of position with the passage of time). Yes. When the mutual difference value is calculated at each position (z) on the Z-axis using Equation 2, the mutual difference value becomes a relatively large value when the heart expands.

基準画像の探索にあたって、基準画像探索部24は、まず、複数の断層画像の中から代表となる基準画像(代表基準画像)を探索する。基準画像探索部24は、図5(A)に示すように、例えば相互差分値が最大となる位置に対応した断層画像データを代表基準画像(代表基準断面)とする。そして、基準画像探索部24は、代表基準画像を起点として、極大の相互差分値に対応した複数の断層画像の中から、仮想周期設定部22(図1)により設定された仮想周期だけ離れた位置に最も近い断層画像を次々に探索する。   In searching for a reference image, the reference image search unit 24 first searches a representative reference image (representative reference image) from a plurality of tomographic images. As shown in FIG. 5A, the reference image search unit 24 sets, for example, tomographic image data corresponding to a position where the mutual difference value is maximum as a representative reference image (representative reference cross section). Then, the reference image search unit 24 is separated from the plurality of tomographic images corresponding to the maximum mutual difference values by the virtual cycle set by the virtual cycle setting unit 22 (FIG. 1), with the representative reference image as a starting point. The tomographic image closest to the position is searched one after another.

まず、図5(A)に示すように、代表基準画像からZ軸方向の正方向と負方向に仮想周期(VHR)だけ離れた位置に最も近い断層画像が探索されて基準画像とされる。次に、基準画像探索部24は、図5(B)に示すように、探索された基準画像から仮想周期(VHR)だけ離れた位置に最も近い断層画像を探索して新たな基準画像とする。図5(B)において、破線の矢印が複数の基準画像(基準断面)の位置を示している。   First, as shown in FIG. 5A, a tomographic image closest to a position separated from the representative reference image by a virtual period (VHR) in the positive and negative directions in the Z-axis direction is searched and used as a reference image. Next, as shown in FIG. 5B, the reference image search unit 24 searches for a tomographic image closest to a position separated from the searched reference image by a virtual period (VHR) to be a new reference image. . In FIG. 5B, broken arrows indicate the positions of a plurality of reference images (reference cross sections).

基準画像探索部24は、代表基準画像を起点として次々に複数の基準画像を探索する。こうして、極大の相互差分値に対応した複数の断層画像の中から、図5(C)に示すように複数の基準画像が探索される。図5(C)において、破線の矢印が複数の基準画像(基準断面)の位置を示している。   The reference image search unit 24 searches a plurality of reference images one after another using the representative reference image as a starting point. In this way, a plurality of reference images are searched from a plurality of tomographic images corresponding to the maximum mutual difference values as shown in FIG. In FIG. 5C, dashed arrows indicate the positions of a plurality of reference images (reference cross sections).

図1に戻り、複数の基準画像が探索されると、再構築処理部20は、複数の基準画像の各々を分割の単位とすることにより、複数の断層画像をいくつかの画像群に分割する。そして、再構築処理部20は、複数の画像群の各々から互いに周期的に対応した複数の断層画像を抽出することにより再構築処理を実現する。再構築処理部20は、前メモリ14に記憶された複数の断層画像データを再構築して後メモリ26に記憶する。   Returning to FIG. 1, when a plurality of reference images are searched, the reconstruction processing unit 20 divides the plurality of tomographic images into several image groups by using each of the plurality of reference images as a unit of division. . And the reconstruction process part 20 implement | achieves a reconstruction process by extracting the some tomographic image corresponding to each other periodically from each of several image groups. The reconstruction processing unit 20 reconstructs a plurality of tomographic image data stored in the previous memory 14 and stores it in the rear memory 26.

図6は、再構築処理部20による処理を説明するための図であり、図6には、前メモリ14に記憶されるデータと後メモリ26に記憶されるデータの対応関係が示されている。図6において、「断層画像Zn(n=1,2,3,・・・,60)」は、Z軸(図2参照)上における座標Znの位置の断層画像データを意味している。   FIG. 6 is a diagram for explaining processing by the reconstruction processing unit 20. FIG. 6 shows a correspondence relationship between data stored in the previous memory 14 and data stored in the rear memory 26. . In FIG. 6, “tomographic image Zn (n = 1, 2, 3,..., 60)” means tomographic image data at the position of the coordinate Zn on the Z axis (see FIG. 2).

前メモリ14には、Z軸方向に沿って次々に形成される複数の走査面に対応した複数の断層画像データが形成された順に記憶されている。つまり、前メモリ14には、いくつかの断層画像に続いて、断層画像Z1,断層画像Z2,・・・,断層画像Z60,・・・の順に複数の断層画像データが記憶されている。   The previous memory 14 stores a plurality of tomographic image data corresponding to a plurality of scan planes formed one after another along the Z-axis direction. That is, the previous memory 14 stores a plurality of tomographic image data in the order of tomographic images Z1, tomographic images Z2,..., Tomographic images Z60,.

再構築処理部20は、複数の基準画像の各々を分割の単位とすることにより、前メモリ14に記憶された複数の断層画像データを複数の画像群に分割する。そして、複数の画像群の各々から、互いに周期的に対応した複数の断層画像データが抽出される。   The reconstruction processing unit 20 divides the plurality of tomographic image data stored in the previous memory 14 into a plurality of image groups by using each of the plurality of reference images as a unit of division. Then, a plurality of tomographic image data periodically corresponding to each other is extracted from each of the plurality of image groups.

図6において、断層画像Z1,断層画像Z15,・・・,断層画像Z51が基準画像探索部24により探索された複数の基準画像である。再構築処理部20は、互いに周期的に対応した複数の断層画像データとして、まず、基準画像である断層画像Z1,断層画像Z15,・・・,断層画像Z51を抽出する。そして、抽出された断層画像Z1,断層画像Z15,・・・,断層画像Z51が一つのデータブロックとなって後メモリ26内に記憶される。   In FIG. 6, the tomographic image Z1, the tomographic image Z15,..., The tomographic image Z51 are a plurality of reference images searched by the reference image searching unit 24. The reconstruction processing unit 20 first extracts a tomographic image Z1, a tomographic image Z15,..., A tomographic image Z51, which are reference images, as a plurality of tomographic image data periodically corresponding to each other. The extracted tomographic image Z1, tomographic image Z15,..., And tomographic image Z51 are stored in the rear memory 26 as one data block.

次に、再構築処理部20は、互いに周期的に対応した複数の断層画像データとして、複数の基準画像の各々に対してZ軸方向の正方向に隣接する複数の断層画像を抽出する。つまり、断層画像Z2,断層画像Z16,・・・,断層画像Z52が抽出され、これらが一つのデータブロックとなって後メモリ26内に記憶される。   Next, the reconstruction processing unit 20 extracts a plurality of tomographic images adjacent in the positive direction in the Z-axis direction to each of the plurality of reference images as a plurality of tomographic image data periodically corresponding to each other. That is, the tomographic image Z2, the tomographic image Z16,..., And the tomographic image Z52 are extracted and stored as one data block in the rear memory 26.

さらに、再構築処理部20は、断層画像Z2,断層画像Z16,・・・,断層画像Z52の各々に対してZ軸方向の正方向に隣接する複数の断層画像を抽出する。こうして、複数の基準画像の各々を起点として、互いに周期的に対応した複数の断層画像のデータブロックが次々に抽出されて、後メモリ26内に記憶される。   Further, the reconstruction processing unit 20 extracts a plurality of tomographic images adjacent to each of the tomographic image Z2, the tomographic image Z16,. In this way, data blocks of a plurality of tomographic images periodically corresponding to each other from each of the plurality of reference images are sequentially extracted and stored in the rear memory 26.

なお、上述した再構築処理によれば、前メモリ14に記憶された複数の断層画像の中で再構築処理に利用されない断層画像がある。例えば、前メモリ14の中の断層画像Z10と断層画像Z15の間の断層画像(Z11〜Z14)などが再構築処理に利用されない断層画像である。   Note that, according to the reconstruction process described above, there are tomographic images that are not used for the reconstruction process among a plurality of tomographic images stored in the previous memory 14. For example, the tomographic images (Z11 to Z14) between the tomographic image Z10 and the tomographic image Z15 in the previous memory 14 are tomographic images that are not used for the reconstruction process.

また、上述した再構築処理では、再構築処理後の後メモリ26内に、複数のデータブロックが形成される。例えば、断層画像Z1,断層画像Z15,・・・,断層画像Z51が一つのデータブロックとなり、断層画像Z2,断層画像Z16,・・・,断層画像Z52が次の一つのデータブロックとなる。再構築処理によって後メモリ26内に形成されるデータブロックの個数は、基準画像に基づいて分割された複数の画像群のうちの、断層画像数が最も少ない画像群の断層画像数に対応している。例えば、図5(C)に示すように複数の基準画像が探索された場合においては、互いに隣り合う2つの基準画像の間隔が最も短い区間に対応した画像群の断層画像数と、図6に示す後メモリ26内のデータブロックの個数が一致する。   In the reconstruction process described above, a plurality of data blocks are formed in the post-memory 26 after the reconstruction process. For example, the tomographic image Z1, the tomographic image Z15,..., And the tomographic image Z51 become one data block, and the tomographic image Z2, the tomographic image Z16,. The number of data blocks formed in the rear memory 26 by the reconstruction process corresponds to the number of tomographic images of the image group having the smallest number of tomographic images among the plurality of image groups divided based on the reference image. Yes. For example, when a plurality of reference images are searched as shown in FIG. 5C, the number of tomographic images in the image group corresponding to the section where the interval between two reference images adjacent to each other is the shortest is shown in FIG. As shown, the number of data blocks in the memory 26 matches.

そこで、例えば、図5(C)に示すように複数の基準画像が探索された後に、互いに隣り合う2つの基準画像の間隔が最も短い区間に対応した画像群の断層画像数eを確認し、それをデータブロック数eとし、再構築処理の結果としてデータブロック数がeに達した時に再構築処理を終了させることにより、再構築処理の無駄を低減(望ましくは完全に除去)することが可能になる。   Therefore, for example, as shown in FIG. 5C, after a plurality of reference images are searched, the number e of tomographic images in the image group corresponding to the interval where the interval between two adjacent reference images is the shortest is confirmed. It is possible to reduce (preferably completely remove) the rebuilding process by setting the number of data blocks as e and terminating the rebuilding process when the number of data blocks reaches e as a result of the rebuilding process. become.

また、図6に示す例においては、基準画像に対応したデータブロックを複数のデータブロックの先頭としているが、基準画像に対応したデータブロックを中心として、複数のデータブロックを形成してもよい。   In the example shown in FIG. 6, the data block corresponding to the reference image is the head of the plurality of data blocks. However, a plurality of data blocks may be formed around the data block corresponding to the reference image.

図7は、再構築処理部20による別の好適な処理を説明するための図であり、図6と同様に、図7には、前メモリ14に記憶されるデータと後メモリ26に記憶されるデータの対応関係が示されている。   FIG. 7 is a diagram for explaining another preferable process by the reconstruction processing unit 20. Like FIG. 6, FIG. 7 shows data stored in the front memory 14 and stored in the rear memory 26. Correspondences between data are shown.

図7に示す例においても、再構築処理部20は、複数の基準画像の各々を分割の単位とすることにより、前メモリ14に記憶された複数の断層画像データを複数の画像群に分割する。その際、予め確認されるデータブロック数eを利用し、各基準画像の位置を中心として各基準画像からZ軸の負方向へe/2枚の断層画像とZ軸の正方向へe/2枚の断層画像とを一つの画像群とする。   Also in the example illustrated in FIG. 7, the reconstruction processing unit 20 divides the plurality of tomographic image data stored in the previous memory 14 into a plurality of image groups by using each of the plurality of reference images as a unit of division. . At that time, using the number of data blocks e confirmed in advance, e / 2 tomographic images in the negative direction of the Z-axis and e / 2 in the positive direction of the Z-axis from each reference image around the position of each reference image One tomographic image is defined as one image group.

例えば、図7に示す前メモリ14内において、基準画像である断層画像Z15を中心として、Z軸の負方向にあるe/2枚の断層画像Z14,断層画像Z13,・・・と、Z軸の正方向にあるe/2枚の断層画像Z16,断層画像Z17,・・・と、により、断層画像Z15を中心としたe枚の断層画像からなる画像群が形成される。同様に、基準画像である断層画像Z32を中心としたe枚の断層画像からなる画像群が形成され、基準画像である断層画像Z51を中心としたe枚の断層画像からなる画像群が形成される。   For example, in the previous memory 14 shown in FIG. 7, e / 2 tomographic images Z14, tomographic images Z13,... In the negative direction of the Z axis centered on the tomographic image Z15 that is the reference image, and the Z axis The e / 2 tomographic images Z16, tomographic images Z17,... In the positive direction form an image group including e tomographic images centered on the tomographic image Z15. Similarly, an image group composed of e tomographic images centered on the tomographic image Z32 that is the reference image is formed, and an image group composed of e sheet tomographic images centered on the tomographic image Z51 that is the reference image. The

そして、図7に示す例においても、複数の画像群の各々から、互いに周期的に対応した複数の断層画像データが抽出される。つまり、再構築処理部20は、各画像群の各々から、まず、Z軸上において最も小さい位置にある断層画像を抽出して一つのデータブロックとし、次に、Z軸上において2番目に小さい位置にある断層画像を抽出して一つのデータブロックとする。この抽出処理をZ軸上において最も小さい位置から順に最も大きい位置まで続けることにより、総数e個のデータブロックが後メモリ26内に形成される。   In the example shown in FIG. 7 as well, a plurality of tomographic image data periodically corresponding to each other is extracted from each of the plurality of image groups. That is, the reconstruction processing unit 20 first extracts the tomographic image at the smallest position on the Z axis from each of the image groups to form one data block, and then the second smallest on the Z axis. A tomographic image at a position is extracted to form one data block. By continuing this extraction process from the smallest position to the largest position on the Z axis, a total number of e data blocks are formed in the rear memory 26.

図7には、総数e個のうちの3個のデータブロックが図示されている。つまり、各々が基準画像の一つ前にある断層画像Z14,断層画像31,・・・,断層画像50で構成されるデータブロックと、基準画像である断層画像Z15,断層画像32,・・・,断層画像51で構成されるデータブロックと、各々が基準画像の一つ後にある断層画像Z16,断層画像33,・・・,断層画像52で構成されるデータブロックが図示されている。   FIG. 7 shows three data blocks out of the total number e. That is, a data block composed of a tomographic image Z14, a tomographic image 31,..., A tomographic image 50 each preceding the reference image, and a tomographic image Z15, a tomographic image 32,. , A data block composed of a tomographic image 51, and a data block composed of a tomographic image Z16, a tomographic image 33,.

このように、図7の例においては、複数の基準画像の各々を中心として、互いに周期的に対応した複数の断層画像のデータブロックが次々に抽出されて、後メモリ26内に記憶される。   As described above, in the example of FIG. 7, data blocks of a plurality of tomographic images periodically corresponding to each other around each of the plurality of reference images are sequentially extracted and stored in the rear memory 26.

なお、図7に示す例では、各基準画像の位置を中心として各基準画像からZ軸の負方向へe/2枚の断層画像とZ軸の正方向へe/2枚の断層画像とを一つの画像群としているが、各基準画像の位置を中心からずらして一つの画像群を形成してもよい。つまり、各基準画像に対してその前後に異なる枚数の断層画像を加えて一つの画像群を形成してもよい。この場合、後メモリ26内において、基準画像に対応したデータブロックが全データブロックの中心からずれることになる。   In the example shown in FIG. 7, e / 2 tomographic images in the negative direction of the Z axis and e / 2 tomographic images in the positive direction of the Z axis from each reference image with the position of each reference image as the center. Although one image group is provided, one image group may be formed by shifting the position of each reference image from the center. That is, one image group may be formed by adding different numbers of tomographic images before and after each reference image. In this case, the data block corresponding to the reference image is shifted from the center of all the data blocks in the rear memory 26.

図1に戻り、三次元画像形成部28は、後メモリ26に記憶された再構築後の複数の断層画像データに基づいて、胎児の心臓を立体的に映し出す三次元画像データを形成する。三次元画像形成部28は、後メモリ26に記憶された一つのデータブロックに基づいて各時相の三次元画像データを形成する。例えば、図6に示す後メモリ26に記憶された断層画像Z1,断層画像Z15,・・・,断層画像Z51に基づいて時相T1の三次元画像データが形成され、断層画像Z2,断層画像Z16,・・・,断層画像Z52に基づいて時相T2の三次元画像データが形成される。または、図7に示す後メモリ26に形成された総数e個のデータブロックの各々に基づいて各時相の三次元画像データが形成されることにより、時相T1〜Teまでの三次元画像データが形成される。図7を利用して説明した再構築処理においては、心拍などの周期的な運動の時相的な基準となる基準画像に対応した三次元画像データが、時相T1〜Teまでの三次元画像データの中心に配置される。つまり、互いの時相関係が最も一致していると考えられる基準画像による三次元画像データを時相T1〜Teの中心に配置することが可能になる。   Returning to FIG. 1, the three-dimensional image forming unit 28 forms three-dimensional image data that three-dimensionally displays the fetal heart based on a plurality of reconstructed tomographic image data stored in the post-memory 26. The three-dimensional image forming unit 28 forms three-dimensional image data for each time phase based on one data block stored in the post-memory 26. For example, three-dimensional image data of time phase T1 is formed based on the tomographic image Z1, tomographic image Z15,..., Tomographic image Z51 stored in the post-memory 26 shown in FIG. ,..., Three-dimensional image data of time phase T2 is formed based on the tomographic image Z52. Alternatively, the three-dimensional image data of the time phases T1 to Te is formed by forming the three-dimensional image data of each time phase based on each of the total number of e data blocks formed in the post-memory 26 shown in FIG. Is formed. In the reconstruction process described with reference to FIG. 7, three-dimensional image data corresponding to a reference image serving as a temporal reference of a periodic motion such as a heartbeat is a three-dimensional image from time phases T1 to Te. Placed in the center of the data. That is, it becomes possible to arrange the three-dimensional image data based on the reference image that is considered to have the most coincident temporal relationship with each other at the center of the time phases T1 to Te.

三次元画像形成部28は、例えば、ボリュームレンダリング法や積算法や投影法などの各種の手法を適用して、各時相ごとに複数の時相に亘って三次元画像データを形成する。こうして、複数の時相に亘って形成された三次元画像データに対応した画像が表示部30に表示され、擬似的にリアルタイムの三次元動画像が表示される。例えば、時相T1から最終時相Teまでの三次元画像データに対応した画像が繰り返し表示されてループ再生が実行されてもよい。   The three-dimensional image forming unit 28 applies various methods such as a volume rendering method, an integration method, and a projection method to form three-dimensional image data over a plurality of time phases for each time phase. Thus, an image corresponding to the three-dimensional image data formed over a plurality of time phases is displayed on the display unit 30, and a pseudo real-time three-dimensional moving image is displayed. For example, an image corresponding to three-dimensional image data from the time phase T1 to the final time phase Te may be repeatedly displayed and loop reproduction may be executed.

上述した実施形態によれば、心拍の周期が不安定な胎児の心臓などを診断対象とする場合においても、周期の変動に応じて適切な基準画像が探索されてデータブロックが再構成されるため、周期の変動に伴う画像の乱れなどが軽減されて(望ましくは完全に除去されて)信頼性の高い表示画像を得ることが可能になる。特に、仮想周期を決定する際に、対象組織の代表部位から、つまり固定的な位置から収集されるエコーデータを利用しているため、移動的に収集されるエコーデータから仮想周期を設定する場合に比べて、仮想周期の設定の精度が向上する。   According to the above-described embodiment, even when a fetal heart or the like whose heartbeat cycle is unstable is to be diagnosed, an appropriate reference image is searched according to the cycle variation and the data block is reconstructed. The disturbance of the image due to the fluctuation of the period is reduced (desirably completely removed), and a display image with high reliability can be obtained. In particular, when determining the virtual period, since the echo data collected from the representative part of the target tissue, that is, from a fixed position is used, the virtual period is set from the echo data collected in a mobile manner Compared to, the accuracy of setting the virtual period is improved.

以上、本発明の好適な実施形態を説明したが、上述した実施形態は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。本発明は、その本質を逸脱しない範囲で各種の変形形態を包含する。例えば、上述した実施形態においては、Z軸方向に沿って走査面をゆっくりと移動させつつ複数の走査面を形成して三次元空間内からエコーデータを収集する低速スキャンを実行している。この低速スキャンの前に、往復約2秒程度で三次元空間内を試験的に高速スキャンさせて、プローブ位置指定の支援を行うようにしてもよい。   As mentioned above, although preferred embodiment of this invention was described, embodiment mentioned above is only a mere illustration in all the points, and does not limit the scope of the present invention. The present invention includes various modifications without departing from the essence thereof. For example, in the above-described embodiment, a low-speed scan is performed in which a plurality of scan planes are formed while slowly moving the scan plane along the Z-axis direction, and echo data is collected from within a three-dimensional space. Prior to this low-speed scanning, the probe position designation may be supported by performing a high-speed scan on a trial basis in about two seconds for a round trip.

10 プローブ、12 ビームフォーマ、16 エラー判定部、20 再構築処理部、22 仮想周期設定部、24 基準画像探索部、28 三次元画像形成部。   10 probe, 12 beam former, 16 error determination unit, 20 reconstruction processing unit, 22 virtual cycle setting unit, 24 reference image search unit, 28 three-dimensional image formation unit.

Claims (9)

周期的に運動する対象組織を含む三次元空間内で超音波を送受するプローブと、
プローブを制御することにより、前記運動の複数の周期に亘って対象組織の代表部位から受信信号を収集し、さらに、前記運動の複数の周期に亘って走査面を移動させつつ前記三次元空間内で複数の走査面を形成する送受信制御部と、
前記代表部位から収集された受信信号に基づいて対象組織に関する運動の仮想周期を設定する仮想周期設定部と、
前記複数の走査面に対応した複数の画像の中から前記仮想周期に対応した間隔で複数の基準画像を探索する基準画像探索部と、
前記複数の基準画像の各々を分割の単位とすることにより前記複数の画像を複数の画像群に分割し、複数の画像群の各々から互いに周期的に対応した複数の画像を抽出する画像再構成部と、
互いに周期的に対応した複数の画像に基づいて対象組織の表示画像を形成する表示画像形成部と、
を有する、
ことを特徴とする超音波診断装置。
A probe for transmitting and receiving ultrasonic waves in a three-dimensional space including a target tissue that moves periodically;
By controlling the probe, received signals are collected from representative regions of the target tissue over a plurality of periods of the movement, and further, the scanning plane is moved over the plurality of periods of the movement while moving within the three-dimensional space. A transmission / reception control unit for forming a plurality of scanning planes;
A virtual period setting unit that sets a virtual period of motion related to the target tissue based on the received signal collected from the representative site;
A reference image search unit that searches a plurality of reference images at intervals corresponding to the virtual period from a plurality of images corresponding to the plurality of scanning planes;
Image reconstruction that divides the plurality of images into a plurality of image groups by extracting each of the plurality of reference images as a unit of division, and extracts a plurality of images periodically corresponding to each of the plurality of image groups And
A display image forming unit that forms a display image of the target tissue based on a plurality of images periodically corresponding to each other;
Having
An ultrasonic diagnostic apparatus.
請求項1に記載の超音波診断装置において、
前記仮想周期設定部は、前記代表部位から収集された複数時相の受信信号に基づいて、時相間における受信信号の変化を反映させた特徴量を算出し、当該特徴量に基づいて前記仮想周期を決定する、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 1,
The virtual cycle setting unit calculates a feature value reflecting a change in the received signal between time phases based on the received signals of a plurality of time phases collected from the representative part, and based on the feature amount, the virtual cycle setting unit To decide,
An ultrasonic diagnostic apparatus.
請求項2に記載の超音波診断装置において、
前記仮想周期設定部は、前記特徴量として各時相ごとに相互差分値を算出し、複数時相に亘って得られた相互差分値の時間的変化に基づいて前記仮想周期を決定する、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 2,
The virtual cycle setting unit calculates a mutual difference value for each time phase as the feature amount, and determines the virtual cycle based on a temporal change of the mutual difference value obtained over a plurality of time phases.
An ultrasonic diagnostic apparatus.
請求項3に記載の超音波診断装置において、
前記仮想周期設定部は、複数時相に亘って得られた相互差分値の極大値に対応した時相の間隔に基づいて前記仮想周期を決定する、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 3.
The virtual cycle setting unit determines the virtual cycle based on a time phase interval corresponding to a maximum value of a mutual difference value obtained over a plurality of time phases.
An ultrasonic diagnostic apparatus.
請求項4に記載の超音波診断装置において、
前記仮想周期設定部は、前記代表部位に設定された超音波ビームから1次元的に収集された複数時相の受信信号を利用する、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 4,
The virtual cycle setting unit uses reception signals of a plurality of time phases collected one-dimensionally from an ultrasonic beam set in the representative part,
An ultrasonic diagnostic apparatus.
請求項4に記載の超音波診断装置において、
前記仮想周期設定部は、前記代表部位に設定された走査面から2次元的に収集された複数時相の受信信号を利用する、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 4,
The virtual cycle setting unit uses reception signals of a plurality of time phases that are two-dimensionally collected from the scanning plane set in the representative part.
An ultrasonic diagnostic apparatus.
請求項1から6のいずれか1項に記載の超音波診断装置において、
前記送受信制御部は、前記代表部位から受信信号を収集した後に、前記三次元空間内で複数の走査面を形成する、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to any one of claims 1 to 6,
The transmission / reception control unit forms a plurality of scanning planes in the three-dimensional space after collecting reception signals from the representative part.
An ultrasonic diagnostic apparatus.
請求項1から6のいずれか1項に記載の超音波診断装置において、
前記送受信制御部は、前記三次元空間内で複数の走査面を形成した後に、前記代表部位から受信信号を収集する、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to any one of claims 1 to 6,
The transmission / reception control unit collects reception signals from the representative part after forming a plurality of scanning planes in the three-dimensional space,
An ultrasonic diagnostic apparatus.
請求項1から8のいずれか1項に記載の超音波診断装置において、
前記基準画像探索部は、前記複数の画像の各々に対応した相互差分値に基づいて、最大の相互差分値に対応した一つの画像を代表基準画像とし、その代表基準画像を起点として、極大の相互差分値に対応した複数の画像の中から、前記仮想周期だけ離れた位置に最も近い画像を次々に探索し、前記複数の基準画像とする、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to any one of claims 1 to 8,
The reference image search unit, based on a mutual difference value corresponding to each of the plurality of images, uses one image corresponding to the maximum mutual difference value as a representative reference image, and starts from the representative reference image as a starting point. From the plurality of images corresponding to the mutual difference value, the image closest to the position separated by the virtual period is searched one after another, and the plurality of reference images are used.
An ultrasonic diagnostic apparatus.
JP2009211577A 2009-04-13 2009-09-14 Ultrasonic diagnostic equipment Active JP5475376B2 (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2009211577A JP5475376B2 (en) 2009-09-14 2009-09-14 Ultrasonic diagnostic equipment
US12/755,045 US8852106B2 (en) 2009-04-13 2010-04-06 Ultrasound diagnostic apparatus
EP10003745.6A EP2241256B1 (en) 2009-04-13 2010-04-07 Ultrasound diagnostic apparatus
CN201010148456.4A CN101897601B (en) 2009-04-13 2010-04-13 Ultrasound diagnostic apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009211577A JP5475376B2 (en) 2009-09-14 2009-09-14 Ultrasonic diagnostic equipment

Publications (2)

Publication Number Publication Date
JP2011056160A true JP2011056160A (en) 2011-03-24
JP5475376B2 JP5475376B2 (en) 2014-04-16

Family

ID=43944499

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009211577A Active JP5475376B2 (en) 2009-04-13 2009-09-14 Ultrasonic diagnostic equipment

Country Status (1)

Country Link
JP (1) JP5475376B2 (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005074225A (en) * 2003-08-28 2005-03-24 Ge Medical Systems Global Technology Co Llc Method and device for obtaining volumetric scan to periodically moving object
JP2007054635A (en) * 2005-08-24 2007-03-08 Medison Co Ltd Apparatus and method for processing three-dimensional ultrasonic image
JP2008142362A (en) * 2006-12-12 2008-06-26 Aloka Co Ltd Ultrasonic diagnostic equipment
JP2009119250A (en) * 2007-10-23 2009-06-04 Aloka Co Ltd Ultrasonic diagnostic apparatus

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005074225A (en) * 2003-08-28 2005-03-24 Ge Medical Systems Global Technology Co Llc Method and device for obtaining volumetric scan to periodically moving object
JP2007054635A (en) * 2005-08-24 2007-03-08 Medison Co Ltd Apparatus and method for processing three-dimensional ultrasonic image
JP2008142362A (en) * 2006-12-12 2008-06-26 Aloka Co Ltd Ultrasonic diagnostic equipment
JP2009119250A (en) * 2007-10-23 2009-06-04 Aloka Co Ltd Ultrasonic diagnostic apparatus

Also Published As

Publication number Publication date
JP5475376B2 (en) 2014-04-16

Similar Documents

Publication Publication Date Title
JP5576036B2 (en) Ultrasonic diagnostic equipment
US8565504B2 (en) Ultrasonic image processing apparatus and ultrasonic image processing method
JP5753798B2 (en) Ultrasonic diagnostic apparatus and its operating method
US7758507B2 (en) Blood flow imaging
JP5172310B2 (en) Ultrasonic diagnostic apparatus and control method thereof
JPH1147132A (en) Three-dimensional image forming method and ultrasonic image processing device
JP2003175041A (en) Ultrasound diagnostic apparatus and image processing method
JP2007020908A (en) Ultrasonic diagnostic equipment and control program of ultrasonic diagnostic equipment
JP5770175B2 (en) Ultrasonic diagnostic apparatus and ultrasonic transmission / reception method
JP2005152648A (en) Method and system for motion-adaptive type spatial synthesis
JP5558931B2 (en) Ultrasonic diagnostic equipment
JP2009011468A (en) Ultrasound diagnosis apparatus
JP5525748B2 (en) Ultrasonic diagnostic equipment
CN106461766B (en) Synchronized phased array data acquisition from multiple acoustic windows
JP5369005B2 (en) Ultrasonic diagnostic equipment
JP5550931B2 (en) Ultrasonic diagnostic equipment
JP5801995B2 (en) Ultrasonic diagnostic equipment
JP5331313B2 (en) Ultrasonic diagnostic equipment
JP2009112374A (en) Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image processing program
JP5475376B2 (en) Ultrasonic diagnostic equipment
JP5396285B2 (en) Ultrasonic diagnostic equipment
JP5461934B2 (en) Ultrasonic diagnostic equipment
JP6193449B1 (en) Ultrasonic diagnostic equipment
JP5461931B2 (en) Ultrasonic diagnostic equipment
JP6063154B2 (en) Ultrasonic diagnostic equipment

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20120807

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130827

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20130830

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20131023

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: 20140128

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140206

R150 Certificate of patent or registration of utility model

Ref document number: 5475376

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250