Source: https://patents.google.com/patent/JP4745133B2/en
Timestamp: 2020-02-19 06:49:35
Document Index: 71819452

Matched Legal Cases: ['art 7', 'art 9', 'art 4', 'art 41', 'art 5', 'art 51', 'art 52', 'art 6', 'art 61', 'art 62', 'art 7']

JP4745133B2 - Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing program - Google Patents
JP4745133B2
JP4745133B2 JP2006149414A JP2006149414A JP4745133B2 JP 4745133 B2 JP4745133 B2 JP 4745133B2 JP 2006149414 A JP2006149414 A JP 2006149414A JP 2006149414 A JP2006149414 A JP 2006149414A JP 4745133 B2 JP4745133 B2 JP 4745133B2
sectional positions
JP2006149414A
JP2007319190A (en
2006-05-30 Application filed by 東芝メディカルシステムズ株式会社, 株式会社東芝 filed Critical 東芝メディカルシステムズ株式会社
2006-05-30 Priority to JP2006149414A priority Critical patent/JP4745133B2/en
2007-12-13 Publication of JP2007319190A publication Critical patent/JP2007319190A/en
2011-08-10 Publication of JP4745133B2 publication Critical patent/JP4745133B2/en
The present invention relates to an ultrasonic diagnostic apparatus that transmits ultrasonic waves into a subject using an ultrasonic probe and acquires a medical image of the subject based on the reflected waves, and a medical image acquired by the ultrasonic diagnostic apparatus. The present invention relates to a medical image processing apparatus and a medical image processing program. In particular, the present invention relates to an ultrasonic diagnostic apparatus, a medical image processing apparatus, and a medical image processing program that are used for evaluating the motion function of a living tissue.
Patent Document 1 discloses an ultrasonic diagnostic apparatus that is used for evaluation of motor function of living tissue. This ultrasonic diagnostic apparatus extracts feature points in a two-dimensional image of a living tissue and designates measurement points based on the feature points. Then, a template having a size including a predetermined number or more of feature points around each designated measurement point is set, and cross correlation processing between the two images is performed using this template, and physical displacement such as displacement and velocity of the tissue is determined. Measurement accuracy is improved by measuring parameters.
JP 2004-313291 A
By the way, the living tissue generally moves three-dimensionally. For example, the heart performs movements such as expansion and contraction and distortion in a direction orthogonal to the movement in the wall thickness direction (wall thickness change, wall thickness distortion, etc.). Since conventional motor function evaluation measures the two-dimensional displacement and velocity of a living tissue, there is a limit to grasping such a three-dimensional movement.
On the other hand, when a three-dimensional motion is to be evaluated using a three-dimensional image of a biological tissue, that is, a plurality of three-dimensional measurement points are specified and a three-dimensional displacement of each measurement point is analyzed. Thus, when trying to evaluate the three-dimensional motion of the living tissue, there arises a problem that the amount of data to be processed increases and the processing time becomes long.
In addition, when an attempt is made to evaluate a motor function using a three-dimensional image, there arises a problem that it becomes difficult to specify an area to be measured for movement on the image. For example, when it is desired to measure a change in the wall thickness of the heart, it is necessary to find a region corresponding to the intima and a region corresponding to the adventitia in a three-dimensional image of the heart, and specify these regions. However, in order to perform such work accurately, a skillful technique is required, and considerable work time is required.
The present invention has been made to solve the above problems, and an ultrasonic diagnostic apparatus, a medical image processing apparatus, and a medical apparatus that can measure a three-dimensional motion of a living tissue in a short time. An object is to provide an image processing program.
Another object of the present invention is to provide an ultrasonic diagnostic apparatus, a medical image processing apparatus, and a medical image processing program that make it possible to easily specify a region for measuring the movement of a living tissue.
In order to achieve the above object, the invention according to claim 1 is directed to an ultrasonic probe, a transmission / reception means for transmitting / receiving ultrasonic waves to the ultrasonic probe, and a reception signal obtained as a result of transmission / reception of the ultrasonic waves. An image generating means for generating image data of a series of tomographic images along a time series for each of the two or more cross-sectional positions of the biological tissue, and the generated image data for each of the two or more cross-sectional positions. Display means for displaying one tomographic image of the series of tomographic images, and specifying means for designating a measurement image region in the displayed one tomographic image for each of the two or more cross-sectional positions. And, for each of the two or more cross-sectional positions, based on the image data of the series of tomographic images, the movement of the living tissue in the designated measurement image region A computing means for computing local motion information indicating a state, and calculating motion information indicating a motion state of the living tissue based on the calculated local motion information for each of the two or more cross-sectional positions. The ultrasonic diagnostic apparatus is characterized in that the display means displays the movement information calculated by the calculation means.
The invention according to claim 15 is a medical image processing apparatus that processes image data of a medical image of a biological tissue obtained by an ultrasonic diagnostic apparatus, and each of the two or more cross-sectional positions of the biological tissue. Storage means for storing image data of a series of tomographic images along a time series, and one tomographic image of the series of tomographic images based on the stored image data for each of the two or more cross-sectional positions Display means for displaying, for each of the two or more cross-sectional positions, designation means for designating a measurement image region in the displayed one tomographic image, and for each of the two or more cross-sectional positions, the series of Based on the image data of the tomographic image, the local motion information indicating the motion state of the living tissue in the designated measurement image region is calculated, and the two or more Calculation means for calculating movement information indicating a movement state of the living tissue based on the calculated local movement information for each of the surface positions, and the display means is calculated by the calculation means It is characterized by displaying exercise information.
The invention according to claim 16 is a computer having storage means for storing image data of a series of tomographic images along a time series at each of two or more cross-sectional positions of a living tissue, and a display means. For each of the two or more cross-sectional positions, based on the stored image data, the display unit functions so as to display one tomographic image of the series of tomographic images. Corresponding to the measurement region designated in the displayed one tomographic image for each, the state of motion of the living tissue in the designated measurement image region based on the image data of the series of tomographic images Based on the calculated local motion information for each of the two or more cross-sectional positions. To function so as to calculate the motion information indicating the state of motion, the so computed the motion information function to display on the display means, a medical image processing program characterized by.
In the invention according to claim 1, claim 15 or claim 16 , one tomographic image is displayed for each of two or more cross-sectional positions of a living tissue, and a measurement image region is designated in the displayed tomographic image. In response to this, the local motion information indicating the motion state of the biological tissue in the measurement image region is calculated, and based on the calculated local motion information in each of the two or more cross-sectional positions, The exercise information is calculated and the exercise information is displayed.
According to the present invention, it is possible to measure the three-dimensional motion of a living tissue by obtaining motion information based on local motion information at each of two or more cross-sectional positions of the living tissue. In addition, since the motion information is obtained considering only the local motion information at two or more cross-sectional positions, it is not necessary to calculate the displacement between the cross-sectional positions. Can be obtained at.
In addition, since the measurement image area is designated for the displayed tomographic image, it is possible to easily designate the measurement image area for acquiring motion information.
Exemplary embodiments of an ultrasonic diagnostic apparatus, a medical image processing apparatus, and a medical image processing program according to the present invention will be described in detail with reference to the drawings.
An embodiment of an ultrasonic diagnostic apparatus according to the present invention will be described. FIG. 1 shows an example of the overall configuration of an ultrasonic diagnostic apparatus according to the present invention. An ultrasonic diagnostic apparatus 1 shown in FIG. 1 is an apparatus used to acquire an image representing a form of a living tissue such as a heart or an image representing a blood flow state, and is a two-dimensional ultrasonic probe 2 and a transmission / reception unit 3. , A signal processing unit 4, an image processing unit 5, an arithmetic processing unit 6, a storage unit 7, a user interface 8, and a control unit 9. Hereinafter, a specific example of each part constituting the ultrasonic diagnostic apparatus 1 will be described.
[Storage unit, user interface, control unit]
First, the memory | storage part 7, the user interface 8, and the control part 9 are demonstrated. The storage unit 7 is configured by a storage device such as a hard disk drive. The storage unit 7 stores in advance a medical image processing program 71 for causing the ultrasonic diagnostic apparatus 1 to execute operations characteristic of the present invention. The storage unit 7 stores various data such as image data of an ultrasonic image and incidental information (DICOM (Digital Imaging and Communications in Medicine) incidental information) attached to the image data.
The user interface 8 is provided with a display unit 81 and an operation unit 82. The display unit 81 corresponds to an example of the “display unit” of the present invention, and is configured by an arbitrary display device such as a liquid crystal display or a CRT (Cathode Ray Tube) display. The display unit 81 displays an image such as an ultrasound image acquired by the ultrasound diagnostic apparatus 1 and information such as DICOM supplementary information of the image.
The operation unit 82 includes an arbitrary operation device or input device such as a mouse, a trackball, a joystick, a control panel, or a keyboard.
In particular, the operation unit 82 functions as a “designating unit” for designating a measurement image region on an ultrasonic image (tomographic image) displayed on the display unit 81. This measurement image region is a region on the tomographic image (actually, image data corresponding to this region) that serves as a reference for measuring the motion state of the living tissue. Although details will be described later, the ultrasonic diagnostic apparatus 1 scans a living tissue three-dimensionally with ultrasound, and repeats this three-dimensional scan to thereby repeat a series of volumes in time series. Data is generated and the image data of the tomographic image of the living tissue is generated based on one of the series of volume data. The operation unit 82 is used to designate a measurement image region for this tomographic image. The ultrasonic diagnostic apparatus 1 acts to measure the motion state of the living tissue by analyzing how the designated measurement image region changes along the time series.
The control unit 9 includes a microprocessor such as a CPU, and controls each unit of the ultrasonic diagnostic apparatus 1 based on the medical image processing program 71. In particular, the control unit 9 performs a process for displaying an image or a screen on the display unit 81. Further, a process for causing the ultrasonic diagnostic apparatus 1 to execute an operation according to the operation signal from the operation unit 82 is performed.
The two-dimensional ultrasonic probe 2 (sometimes simply referred to as the ultrasonic probe 2) includes a plurality of ultrasonic transducers arranged two-dimensionally (for example, in a matrix (lattice)), as in the past. (Not shown). The plurality of ultrasonic transducers are individually driven by a transmission / reception unit 3 described later.
2 and 3 show an aspect of ultrasonic scanning by the two-dimensional ultrasonic probe 2. As shown in FIG. 2A, the ultrasonic probe 2 scans ultrasonic waves (beams) output from the array surface of the ultrasonic transducers in the main scanning direction X, thereby causing a two-dimensional radial (fan-shaped) shape. A scan plane P is formed. Further, the ultrasonic probe 2 scans the ultrasonic waves in the sub-scanning direction Y orthogonal to the main scanning direction X, so that a plurality of fan shapes arranged in the sub-scanning direction Y as shown in FIG. , Pn of the two-dimensional scan planes P1, P2,. Thereby, a three-dimensional scan region R as shown in FIG. 3 is formed.
The sub-scanning direction Y corresponds to the “predetermined direction” of the present invention, and the main scanning direction X corresponds to the “direction orthogonal to the predetermined direction” of the present invention. The two-dimensional scan planes P1 to Pn are formed at “plurality (n) positions along a predetermined direction” of the present invention.
The transmission / reception unit 3 receives an echo signal (reception signal) output from the transmission unit that generates an ultrasonic wave by supplying an electrical signal to the ultrasonic probe 2 and the ultrasonic probe 2 that has received the reflected wave of the ultrasonic wave. (Not shown).
The reception unit in the transmission / reception unit 3 includes a preamplifier circuit, an A / D conversion circuit, and a reception delay / addition circuit (not shown). The preamplifier circuit amplifies the echo signal output from each ultrasonic transducer of the ultrasonic probe 2 for each reception channel. The A / D conversion circuit performs A (analog) / D (digital) conversion on the amplified echo signal. The reception delay / adder circuit gives a delay time necessary for determining the reception directivity to the echo signal after A / D conversion, and adds the delay time. By this addition processing, the reflection component from the direction corresponding to the reception directivity is emphasized. The added signal may be referred to as “RF data (or raw data)”. The transmission / reception unit 3 inputs the acquired RF data to the signal processing unit 4.
The signal processing unit 4 performs signal processing for imaging the amplitude information of the echo signal based on the RF data input from the transmission / reception unit 3. Data generated by the signal processing unit 4 is sent to the control unit 9 and displayed on the display unit 81 of the user interface 8 or input to the image processing unit 5. The signal processing unit 4 mainly includes a B-mode processing unit 41, a Doppler processing unit 42, and a CMF processing unit 43.
The B (Brightness) mode processing unit 41 generates B-mode ultrasonic raster data based on the RF data. More specifically, the B-mode processing unit 41 performs band pass filter processing on the RF data, detects an envelope of the output signal, and compresses the detected data by logarithmic transformation. Apply. As a result, for each of the two-dimensional scan planes P1 to Pn, image data of a tomographic image in which the signal intensity is expressed by brightness is generated. The B-mode processing unit 41 corresponds to an example of the “first tomographic image generation unit” of the present invention.
The Doppler processing unit 42 generates blood flow information in the living tissue by, for example, a pulse Doppler method (PW Doppler method) or a continuous wave Doppler method (CW Doppler method).
In the pulse Doppler method, a pulse wave is used to detect an ultrasonic frequency displacement (Doppler displacement frequency component) caused by a Doppler effect due to blood flow at a specific depth (distance from the ultrasonic probe 2). Can do. As described above, since the pulse Doppler method has a good distance resolution, the pulse Doppler method is suitably used for tissue measurement of a specific site or blood flow depth measurement. When this pulse Doppler method is applied, the Doppler processing unit 42 extracts a Doppler displacement frequency component from the RF data input from the transmission / reception unit 3 by phase detection of a signal in a blood flow observation region having a predetermined size. Further, FFT (Fast Fourier Transform) processing is performed to generate data indicating the Doppler frequency distribution representing the blood flow velocity in the blood flow observation region.
Further, in the continuous wave Doppler method, unlike the pulse Doppler method, by using a continuous wave, the ultrasonic wave transmission / reception direction (radial direction in the fan-shaped two-dimensional scan plane P shown in FIG. Although a signal on which the Doppler displacement frequency component is superimposed, that is, a signal reflecting all the blood flow conditions on the ultrasonic path, is obtained, there is an advantage that the measurement speed is excellent. When this continuous wave Doppler method is applied, the Doppler processing unit 42 detects the Doppler displacement frequency component by performing phase detection on the RF data input from the transmitting / receiving unit 3 on the sample line for blood flow observation. Extraction is performed, and further FFT processing is performed to generate data indicating a Doppler frequency distribution representing the blood flow velocity on the sample line.
A CFM (Color Flow Mapping) processing unit 43 operates when performing a color flow mapping method in which blood flow information of a living tissue is superimposed on a monochrome B-mode image in color and displayed in real time. The displayed blood flow information includes blood flow velocity, dispersion, power, and the like. This blood flow information is obtained as binarized information. More specifically, the CFM processing unit 43 includes a phase detection circuit, an MTI (Moving Target Indication) filter, an autocorrelator, a flow velocity / dispersion calculator, and the like. The CFM processing unit 43 separates the morphological signal reflecting the morphology of the biological tissue and the blood flow signal reflecting the blood flow by high-pass filter processing (MTI filter processing), and performs blood flow velocity by autocorrelation processing. Blood flow information such as dispersion and power is obtained for a plurality of positions. In addition, non-linear processing for reducing the shape signal may be performed.
The image processing unit 5 performs various image processing based on the data generated by the signal processing unit 4. For example, the image processing unit 5 has a DSC (Digital Scan Converter), and converts data synchronized with ultrasonic scanning generated by the signal processing unit 4 into display data (TV scanning method data). That is, scan conversion processing is performed.
The volume data generation unit 51 performs interpolation processing on the image data of each of the two-dimensional scan planes P1 to Pn generated by the B mode processing unit 41 of the signal processing unit 4 to generate volume data (voxel data). The volume data generation unit 51 corresponds to an example of the “volume data generation unit” of the present invention, and includes, for example, a DSC or a microprocessor.
When displaying a pseudo three-dimensional image based on volume data, the image processing unit 5 performs volume rendering processing, MIP (Maximum Intensity Projection) processing, and the like on the volume data.
An MPR (MultiPlanar Reconstruction) processing unit 52 executes cross-section conversion processing based on the volume data generated by the volume data generation unit 51, thereby generating image data of a tomographic image at an arbitrary cross-section. The MPR processing unit 52 corresponds to an example of the “second tomographic image generation unit” of the present invention, and includes, for example, a DSC or a microprocessor. The B-mode processing unit 41, the volume data generation unit 51, and the MPR processing unit 52 of the signal processing unit 4 function as an example of the “image generation unit” of the present invention.
[Operation processing unit]
Based on the image data of the tomographic image generated by the MPR processing unit 52 of the image processing unit 5, the arithmetic processing unit 6 performs local motion information indicating the state of local motion of the biological tissue, or more global motion. It calculates motion information indicating the state and functions as an example of the “calculation means” of the present invention.
For example, when the living tissue is the heart, the local motion information obtained by the arithmetic processing unit 6 includes a change in the thickness of the heart wall, the speed of the change, a strain of the heart wall motion, a strain, and the like. There are a rate, a rotation angle of the intima and epicardium of the heart wall, a speed (rotation speed) of the rotation angle, a relative rotation angle of the intima and adventitia (details will be described later).
The exercise information includes, for example, the torsional motion of the heart wall, the speed of the torsional motion, the expansion / contraction (shortening), the speed of the expansion / contraction, the strain of the heart wall motion, the strain rate, the relative rotational gradient, etc. (It will be described later.)
The arithmetic processing unit 6 includes a microprocessor such as a CPU. The calculation processing unit 6 is provided with a displacement calculation unit 61 and a motion information calculation unit 62.
(Displacement calculation unit)
The displacement calculation unit 61 calculates the displacement along the time series of the measurement image region (the living tissue) by tracking the measurement image region specified in the tomographic image by the operation unit 82 along the time series. The displacement of the measurement image region corresponds to an example of “local motion information” of the present invention.
The displacement speed can be obtained by dividing the displacement between two-dimensional or three-dimensional images along the time series by the time interval (frame interval) between the images. Conversely, the displacement between the images can be obtained by multiplying the displacement speed of the measurement image region by the time interval between the images. That is, if the time interval between images is known, the displacement and the speed can be considered synonymous. In this sense, in the present invention, displacement and speed may be identified with each other.
The operation of the displacement calculator 61 will be described in more detail. As described above, the ultrasound diagnostic apparatus 1 generates a series of volume data along a time series of biological tissue, and a tomography based on one of the volume data (volume data at a certain time (time phase)). Generate image data of an image. A measurement image region is designated for this tomographic image. The displacement calculation unit 61 tracks the amount of displacement of the measurement image area in the volume data at other times (time phases) along a time series. Such tracking processing of the displacement of the measurement image region in the volume data along the time series may be referred to as “three-dimensional tracking”.
In addition, the displacement calculation unit 61 tracks in time series how much the measurement image region specified in the tomographic image is displaced in the tomographic image at another time (time phase) at the same cross-sectional position as the tomographic image. You can also Such tracking processing is sometimes referred to as “two-dimensional tracking”. In this two-dimensional tracking, tomographic image data along the time series at the cross-sectional position is generated based on volume data along the time series, and displacement is tracked in the image data of the tomographic images along the time series. Can be executed. Further, two-dimensional tracking may be performed by tracking the displacement in the cross-sectional position in the volume data along the time series.
Such tracking processing can be performed by a method similar to the conventional method. For example, in the two-dimensional tracking, similar to the method described in Patent Document 1 described above, feature points are extracted from measurement image regions designated for each of a plurality of tomographic images along a time series, and based on the feature points. Specify the measurement point. Then, a template having a size including a predetermined number or more of feature points is set around the measurement point, and cross correlation processing (pattern matching processing) is performed on two tomographic images (measurement image regions) using the template. Calculate the displacement of each measurement point.
In three-dimensional tracking, a three-dimensional template is similarly set based on volume data, and pattern matching processing is performed on two volume data having different times (time phases) using the three-dimensional template. To calculate the displacement of each measurement point.
Note that feature points and measurement points may be specified only on a measurement image region (for example, an intima position image m1 shown in FIG. 5 described later), or within a region having the measurement image region as a boundary ( For example, it may be configured to designate an area other than the measurement image area, such as an image area corresponding to a cross section of the heart wall surrounded by the intima position image m1 and the epicardial position image M1, or a vicinity area of the measurement image area. Good. In any case, it is possible to apply any tracking method capable of calculating the displacement along the time series of the designated measurement image region.
The displacement of the measurement point obtained by two-dimensional tracking or three-dimensional tracking can be used as it is as the displacement of the measurement image region. Further, based on the displacement of these measurement points, the displacement of the boundary of the measurement image region (for example, the intima position images m1 to m3 and the epicardial position images M1 to M3 shown in FIG. 6 described later) is calculated, and this boundary Can be used as the displacement of the measurement image region.
As described above, the displacement calculation unit 61 designates a plurality of measurement points in the measurement image region designated for one tomographic image, and obtains the positions of these measurement points for the tomographic image of each frame. Then, the displacement along the time series of each measurement point is calculated based on the position of the measurement point of the tomographic image of each frame, and the displacement of the designated measurement image area is calculated based on the displacement of this measurement point. It acts to do.
(Exercise information calculator)
The motion information calculation unit 62 performs a process of calculating motion information indicating the (more global) motion state of the living tissue based on the displacement of the measurement image region calculated by the displacement calculation unit 61. A specific example of the exercise information calculation process by the exercise information calculation unit 62 will be described later.
[Mode of operation]
An example of the operation mode of the ultrasonic diagnostic apparatus 1 according to this embodiment having the above-described configuration will be described with reference to FIGS. Here, a case where the motion state of the heart is evaluated will be described. Hereinafter, acquisition of an ultrasonic image and designation of a measurement image region will be described with reference to FIGS. 4 to 6, and then tracking processing of the measurement image region and calculation processing of motion information will be described with reference to FIGS. 7 to 14. I will explain.
[Acquisition of ultrasonic image and specification of measurement image area]
First, an ultrasound image of the heart to be evaluated for motion status is acquired. For this purpose, first, the transmitting / receiving unit 3 controls the control unit 9 by applying an ultrasonic probe to the body surface in the vicinity of the subject's heart (generally, in the vicinity of the heart apex) and performing a predetermined operation. Based on the above, the ultrasonic probe 2 is controlled to perform a three-dimensional ultrasonic scan of the heart (the ultrasonic scan shown in FIGS. 2 and 3) (S01).
The ultrasonic diagnostic apparatus 1 repeatedly executes this three-dimensional ultrasonic scan. At this time, it is desirable that the three-dimensional ultrasonic scan is repeated for a time equal to or longer than one cardiac cycle (one cycle of the heart beat (for example, from the R wave to the next R wave in the electrocardiogram)).
Next, the transmission / reception unit 3 converts the echo signals sequentially transmitted from the ultrasonic probe 2 into RF data, and sequentially inputs them to the signal processing unit 4. The B-mode processing unit 41 of the signal processing unit 4 sequentially generates image data of tomographic images on the respective two-dimensional scan planes P1 to Pn based on the RF data (S02). The generated image data is input to the image processing unit 5. When the volume data is obtained first, MPR processing or the like may be performed on the volume data to generate tomographic image data.
Subsequently, the volume data generation unit 51 of the image processing unit 5 sequentially performs interpolation processing on the image data of the tomographic images on the two-dimensional scan planes P1 to Pn, thereby corresponding to one three-dimensional ultrasonic scan. Volume data in the three-dimensional scan region R is sequentially generated (S03). Thereby, a plurality of volume data in time series are obtained.
Next, the MPR processing unit 52 generates image data of a tomographic image at a predetermined cross-sectional position of the heart based on one of the generated plurality of volume data. In this embodiment, image data of an apex four-chamber image (also referred to as apex four-chamber tomography) and image data of an apex two-chamber image (also referred to as apex two-chamber tomography) are generated (S04). Here, the apex four-chamber image and the apex two-chamber image are tomographic images at cross-sectional positions along the long axis direction of the heart, and the cross-sectional positions thereof are orthogonal to each other.
The control unit 9 causes the display unit 81 to display a tomographic image (apical four-chamber image and apical two-chamber image) based on the image data generated in step S04 (S05). FIG. 5 shows an example of the display mode at this time. The display unit 81 (display screen thereof) in the figure is provided with cross-sectional position designation image display units 81A and 81B on which tomographic images based on the image data generated in step S04 are displayed. The apex four-chamber image is displayed on the cross-sectional position designation image display unit 81A of the figure, and the apex two-chamber image is displayed on the cross-section position designation image display unit 81B.
In the case of ultrasonic diagnosis based on electrocardiogram synchronization, an electrocardiogram is displayed on the electrocardiogram display unit 81F of the display unit 81. The electrocardiogram display unit 81F displays a time cursor T indicating the time (time phase, time) when the tomographic images displayed on the cross-section position designation image display units 81A and 81B are acquired. In FIG. 5, the time cursor T is arranged at the time phase of the R wave of the electrocardiogram. Here, the time cursor T is configured to move (drag and drop, etc.) in the time direction (lateral direction) of the electrocardiogram, and a tomographic image at the time (time phase) to which the time cursor T is moved is generated from the volume data. Then, it may be configured to display on the cross-section position designation image display portions 81A and 81B.
Cross-sectional position specifying cursors C1, C2, and C3 are provided in the horizontal position of the cross-sectional position specifying image display portion 81B. The user, for example, operates the mouse of the operation unit 82 while observing the apex four-chamber image or the apex two-chamber image, and drags and drops the cross-section position designation cursors C1 to C3 in the vertical direction (long axis direction of the heart). Thus, the cross-sectional position is designated (S06).
For example, based on the coordinates on the display screen of the apex four-chamber image (and / or apex two-chamber image) and the coordinates on the display screen of the cross-section position designation cursors C1 to C3, the control unit 9 images the apex four-chamber image. The coordinates of the designated cross-sectional position in the data are determined and sent to the image processing unit 5.
Based on the coordinate information sent from the control unit 9 and the plurality of volume data generated in step S03, the MPR processing unit 52 is a tomographic image having a cross section in the short axis direction of the heart at each specified cross section position. Are respectively generated (S07). At this time, the MPR processing unit 52 generates image data of a tomographic image at each designated cross-sectional position for each of a plurality of volume data.
The control unit 9 displays on the tomographic image display units 81C, 81D, and 81E of the display unit 81 the tomographic image based on the image data generated from the same volume data as in step S04 among the image data generated in step S07. (S08). In step S07, the image data related to the display process in step S08 may be preferentially generated, and the other image data generation processes may be executed in the background of the process after step S08.
In the example shown in FIG. 5, the cross-sectional position designation cursors C1 to C3 are designated at the heart apex level, the papillary muscle level, and the heart base level, respectively. In this case, the tomographic image display unit 81C displays a tomographic image (apex apex short axis image) G1 at the apex level specified by the cross-section position specifying cursor C1. Further, the tomographic image display portion 81D displays a tomographic image (papillary muscle short axis image) G2 at the papillary muscle level designated by the cross-section position designation cursor C2. Further, the tomographic image display portion 81E displays a tomographic image (cardiac base short axis image) G3 at the heart base level designated by the cross-section position designation cursor C3.
The user operates the operation unit 82 to designate a measurement image region on the tomographic images displayed on the tomographic image display units 81C, 81D, 81E (S09). This operation is performed, for example, by dragging the mouse and inputting a boundary indicating the measurement image region on the tomographic image.
FIG. 6 shows an example of a measurement image region designation mode for the apex short-axis image G1, the papillary muscle short-axis image G2, and the base basal short-axis image G3 in FIG. In addition, the inner dotted line shown in each short-axis image G1 to G3 indicates the intima of the heart wall, and the outer dotted line indicates the adventitia. The user operates the operation unit 82 while observing the displayed short-axis images G1 to G3, and inputs a line indicating the measurement image region so as to trace the inner membrane and the outer membrane, respectively.
Thereby, as shown in FIG. 6, the apex short-axis image G1 includes an intima position image m1 obtained by tracing the intima of the heart wall, and an adventitial position image M1 obtained by tracing the adventitia. Is entered. In addition, the intima position image m2 obtained by tracing the intima of the heart wall and the adventitia position image M2 obtained by tracing the epicardium are input to the papillary muscle short-axis image G2. In addition, an intima position image m3 obtained by tracing the intima of the heart wall and an adventitia position image M3 obtained by tracing the epicardium are input to the short base image G3.
This completes the acquisition of the ultrasound image and the designation of the measurement image region.
[Measurement image area tracking processing and motion information calculation processing]
Next, the tracking process (displacement calculation process) of the measurement image region and the calculation process of motion information will be described with reference to FIGS. Hereinafter, it demonstrates for every exercise | movement information used as acquisition object. Each process described below is executed following step S09 in the flowchart of FIG. A plurality of processes described below may be performed sequentially or in parallel. Further, the same processes executed in different processes do not need to be executed individually.
(Exercise information: Twist exercise)
First, processing for acquiring the state of the torsional motion of the heart wall as motion information will be described with reference to FIGS. When a measurement image region is designated on the tomographic images displayed on the tomographic image display units 81C, 81D, and 81E (S09), the displacement calculation unit 61 is generated for each of a plurality of volume data along the time series. By performing two-dimensional tracking of the intima position image m1 based on the image data (S07) of the tomographic image, the inner axis with the direction orthogonal to the cross section of the apex short axis image G1 (the major axis direction of the heart) as an axis. The rotation angle (local motion information) of the film position image m1 is calculated. Similarly, for the papillary muscle short-axis image G2 and the cardiac base short-axis image G3, the rotation angles (local motion information) of the intima position images m2 and m3 about the long-axis direction of the heart are calculated (S11). Instead of the intima position images m1, m2, and m3, the rotation angles of the intima position images M1, M2, and M3 may be calculated.
At this time, for example, for each time phase, the displacement calculation unit 61 rotates the intima position images m1, m2, and m3 as rotation angles with respect to the time phase (reference time phase) to which the intima position image m1 or the like is input in step S09. Calculate the angle. Further, the rotation angles of the intima position images m1, m2, and m3 in adjacent frames (that is, continuous frames) along the time series may be sequentially calculated.
The motion information calculation unit 62 calculates the difference (relative rotation angle) between the rotation angle of the intima position image m1 and the rotation angle of the intima position image m2 (S12). Similarly, the difference (relative rotation angle) between the rotation angle of the intima position image m2 and the rotation angle of the intima position image m3 is calculated (S13). These relative rotation angles correspond to an example of “difference information” of the invention.
The processing of steps S12 and S13 will be specifically described with reference to FIG. In the tomographic image display portions 81C, 81D, 81E, for example, the counterclockwise direction is defined as a positive rotation direction (+ θ direction). The rotation angle of the intima position image m1 is θ1, the rotation angle of the intima position image m2 is θ2, and the rotation angle of the intima position image m3 is θ3.
At this time, the relative rotation angle Δθ12 calculated in step S12 is calculated by Δθ12 = θ1-θ2 (or θ2-θ1). Further, the relative rotation angle Δθ23 calculated in step S13 is calculated by Δθ23 = θ2-θ3 (or θ3-θ2).
The relative rotation angle Δθ12 obtained in step S12 is information reflecting the state (magnitude) of the torsional motion of the heart wall between the cross-sectional position of the apex short-axis image G1 and the cross-sectional position of the papillary muscle short-axis image G2. is there. That is, when the relative rotation angle Δθ12 = 0 (θ1 = θ2), the heart wall is rotated by the same angle in the same direction at any position between these cross-sectional positions, and there is no twist in the rotation direction. Can be considered a thing.
On the other hand, when | Δθ12 | ≠ 0, there is a difference in rotation angle between these cross-sectional positions, and the heart wall is twisted in the rotation angle direction. This heart wall twist increases as the absolute value of the relative rotation angle Δθ12 increases. For example, when the sign of θ1 and the sign of θ2 are different, that is, when the rotation direction of the intima position image m1 is opposite to the rotation direction of the intima position image m2, the absolute value of the relative rotation angle Δθ12 is compared. Become bigger.
Similarly, the relative rotation angle Δθ23 obtained in step S13 is information that reflects the magnitude of the torsional motion of the heart wall between the cross-sectional position of the papillary muscle short-axis image G2 and the cross-sectional position of the base proximal short-axis image G3. is there.
The control unit 9 causes the display unit 81 to display the relative rotation angles Δθ12 and Δθ23 calculated in steps S12 and S13 as exercise information indicating the magnitude of the torsional motion of the heart wall (S14). By referring to the displayed relative rotation angles Δθ12 and Δθ23, the user can grasp the magnitude of the torsional motion of the heart wall. Here, relative rotation angles can be calculated for the intima and epicardium of the heart wall, and the magnitude of the twisting motion can be evaluated based on the two relative rotation angles (for example, the average value of the two relative rotation angles). Etc.).
The speed of the torsional motion of the heart wall between the intima position images m1 and m2 can be obtained by differentiating the relative rotation angle Δθ12 with respect to time. Similarly, the speed of the torsional motion of the heart wall between the intima position images m2 and m3 can be obtained by differentiating the relative rotation angle Δθ23 with respect to time. It is possible to configure the display unit 81 to display these speeds. Here, “differentiation” includes processing for dividing the relative rotation angle by a time interval between frames for which the relative rotation angle is obtained, in addition to a normal differentiation operation.
(Motion information: Relative rotational gradient)
Processing for acquiring the relative rotation gradient of the heart wall as motion information will be described with reference to FIGS. This relative rotational gradient is motion information indicating the degree of torsional motion of the heart wall.
First, the displacement calculation unit 61 performs the rotation angle θ1 of the intima position image m1 of the apex short axis image G1 and the rotation of the intima position image m2 of the papillary muscle short axis image G2 in the same manner as in step S11 of FIG. The angle θ2 and the rotation angle θ3 of the intima position image m3 of the core base short axis image G3 are respectively calculated (S21).
Next, the motion information calculation unit 62 calculates the relative rotation angle Δθ12 between the rotation angle θ1 of the intima position image m1 and the rotation angle θ2 of the intima position image m2 in the same manner as steps S12 and S13 (S22). ), The rotation angle θ2 of the intima position image m2, the rotation angle θ3 of the intima position image m3, and the relative rotation angle Δθ23 are calculated (S23).
Subsequently, the motion information calculation unit 62 calculates a distance d12 between the apex short axis image G1 and the papillary muscle short axis image G2 (S24), and the papillary muscle short axis image G2 and the base proximal short axis image G3 are calculated. The distance d23 is calculated (S25). The distances d12 and d23 are calculated based on, for example, the coordinates of the cross-sectional positions of the apex short axis image G1, the papillary muscle short axis image G2, and the base proximal short axis image G3 obtained by the control unit 9 after step S06. can do.
Further, as shown in FIG. 10, the motion information calculation unit 62 divides the relative rotation angle Δθ12 obtained in step S22 by the distance d12 obtained in step S24 to obtain an intima position image m1 and an intima position image m2. Relative rotational gradient δθ12 = Δθ12 / d12 is calculated (S26). Similarly, the motion information calculation unit 62 divides the relative rotation angle Δθ23 obtained in step S23 by the distance d23 obtained in step S25, and performs relative rotation between the intima position image m2 and the intima position image m3. The gradient δθ23 = Δθ23 / d23 is calculated (S27).
The control unit 9 causes the display unit 81 to display the relative rotation gradients δθ12 and δθ23 calculated in steps S26 and S27 as exercise information indicating the degree of torsional motion of the heart wall (S28).
The relative rotational gradient δθ12 indicates the magnitude of twist per unit distance between the intima at the apex level and the intima at the papillary muscle level. The relative rotational gradient δθ23 indicates the magnitude of twist per unit distance between the intima of the intima at the papillary muscle level and the intima at the heart base level. That is, the relative rotational gradients δθ12 and δθ23 are motion information reflecting the degree of twist of the heart wall (intima). The user can grasp the degree of the torsional motion of the heart wall by referring to the displayed relative rotational gradients δθ12 and δθ23. It is also possible to calculate a relative rotational gradient for each of the intima and outer membranes of the heart wall and evaluate the degree of twisting motion based on the two relative rotational gradients (for example, taking an average value of the two relative rotational gradients). Such.).
(Exercise information: Stretching in the long axis direction)
Processing for acquiring expansion / contraction (shortening) in the major axis direction of the heart wall as motion information will be described with reference to FIGS. 11 and 12.
First, the displacement calculation unit 61, based on a plurality of time-series volume data (S03), an intima position image m1 of the apex short axis image G1, an intima position image m2 of the papillary muscle short axis image G2, By performing three-dimensional tracking of the intima position image m3 of the basal part short axis image G3, the measurement image region is designated for each of the intima position image m1, the intima position image m2, and the intima position image m3. The three-dimensional displacements (Δx1, Δy1, Δz1), (Δx2, Δy2, Δz2), (Δx3, Δy3, Δz3) of the tomographic image are calculated (S31). These displacements correspond to an example of “local motion information”. Instead of the intima position images m1, m2, and m3, three-dimensional displacements of the epicardial position images M1, M2, and M3 may be calculated.
At this time, for example, for each time phase, the displacement calculation unit 61 sets 3 of intima position images m1, m2, and m3 as three-dimensional displacements with respect to the reference time phase to which the intima position image m1 or the like is input in step S09. Calculate the dimensional displacement. Alternatively, the three-dimensional displacement of the intima position images m1, m2, and m3 in successive frames may be calculated sequentially.
Here, Δx and Δy represent displacements in the X direction (any direction is defined as + X direction) and Y direction shown in FIGS. 2 and 3, respectively. The plane including the X direction and the Y direction is parallel to the cross sections of the apex short axis image G1, the papillary muscle short axis image G2, and the basal base short axis image G3. Δz represents a displacement in a Z direction orthogonal to the X direction and the Y direction (for example, the direction of the apex as viewed from the papillary muscle level is the −Z direction and the direction of the base is the + Z direction). This Z direction is parallel to the long axis direction of the heart.
The displacement calculation unit 61 further moves in the Z direction (major axis direction) from each of these three-dimensional displacements (Δx1, Δy1, Δz1), (Δx2, Δy2, Δz2), (Δx3, Δy3, Δz3). The displacements Δz1, Δz2, and Δz3 are extracted (S32).
Here, it is assumed that the cross section of the apex short axis image G1 and the like is parallel to the XY plane. However, even if the cross section of the apex short axis image G1 and the like is not parallel to the xy plane, By projecting the displacement (vector) in the Z direction, the displacement in the Z direction can be easily calculated.
Next, as shown in FIG. 12, the motion information calculating unit 62 calculates the difference between the displacement Δz1 in the Z direction of the intima position image m1 and the displacement Δz2 in the Z direction of the intima position image m2. Thus, the expansion and contraction Δz12 = Δz1−Δz2 (or Δz2−Δz1) of the heart wall between the apex level and the papillary muscle level is obtained (S33). Similarly, the motion information calculation unit 62 calculates the difference between the displacement Δz2 in the Z direction of the intima position image m2 and the displacement Δz3 in the Z direction of the intima position image m3, thereby calculating the papillary muscle level and the heart. Stretching of the heart wall between the base level and Δz23 = Δz2−Δz3 (or Δz3−Δz2) is obtained (S34). These expansions Δz12 and Δz13 correspond to an example of “difference information”.
The control unit 9 causes the display unit 81 to display the heart wall stretching Δz12 and Δz23 calculated in steps S33 and S34 as exercise information indicating the size of the heart wall stretching (S35). The user can grasp the magnitude of the expansion and contraction of the heart wall by referring to the displayed expansion and contraction Δz12 and Δz23 of the heart wall.
It is also possible to calculate the expansion and contraction of the heart wall for each of the intima and the epicardium of the heart wall and evaluate the expansion / contraction magnitude based on the two expansion / contraction values (for example, the average value of the two expansion / contraction values is Etc.)
Further, by differentiating the expansion / contraction Δz12 with respect to time, the speed of the expansion / contraction motion of the heart wall between the intima position images m1 and m2 can be obtained. Similarly, by differentiating the expansion / contraction Δz23 with time, the speed of the expansion / contraction motion of the heart wall between the intima position images m2 and m3 can be obtained. It is possible to configure the display unit 81 to display these speeds. Here, “differentiation” has the same meaning as described above.
(Exercise information: Long axis strain)
Processing for acquiring the strain in the long axis direction of the heart wall as motion information will be described with reference to FIGS. This strain is information indicating the degree of distortion of the heart wall, and indicates the distortion state of the heart wall.
First, similarly to the case of obtaining expansion and contraction in the long axis direction, the displacement calculation unit 61 designates a measurement image region for each of the intima position image m1, the intima position image m2, and the intima position image m3. The tomographic image three-dimensional displacement (Δx1, Δy1, Δz1), (Δx2, Δy2, Δz2), (Δx3, Δy3, Δz3) is calculated (S41), and the Z direction (long) is calculated from these three-dimensional displacements. The displacements Δz1, Δz2, Δz3 in the axial direction are extracted (S42).
Next, the motion information calculation unit 62 calculates the expansion / contraction Δz12 = Δz1−Δz2 of the heart wall between the apex level and the papillary muscle level in the same manner as described above (S43), and the papillary muscle level and the heart The expansion / contraction Δz23 = Δz2−Δz3 of the heart wall between the base level is calculated (S44).
Further, the motion information calculation unit 62 performs the apex short axis image G1, the papillary muscle short axis image G2, and the base basal short axis image G3 in which the measurement image regions are designated in the same manner as in steps S24 and S25 of FIG. A distance d12 between the apex short axis image G1 and the papillary muscle short axis image G2 is calculated (S45), and a distance d23 between the papillary muscle short axis image G2 and the base basal short axis image G3 is calculated (S46). ).
The motion information calculation unit 62 further divides the expansion / contraction Δz12 calculated in step S43 by the distance d12 calculated in step S45, thereby causing a strain δz12 = Δz12 / in the major axis direction between the apex level and the papillary muscle level. d12 is calculated (S47). In addition, the motion information calculation unit 62 divides the expansion / contraction Δz23 calculated in step S44 by the distance d23 calculated in step S46, so that the strain δz23 in the long axis direction between the papillary muscle level and the basal level is equal to Δz23 / d23 is calculated (S48).
The control unit 9 causes the display unit 81 to display the heart wall strains δz12 and δz23 calculated in steps S47 and S48 as motion information indicating the size of the heart wall strain (S49). The user can grasp the magnitude of the distortion of the heart wall by referring to the displayed strains δz12 and δz23 of the heart wall.
It is also possible to calculate a strain for each of the intima and outer membranes of the heart wall and evaluate the magnitude of the distortion based on the two strain values (for example, taking the average value of the two strain values). ).
(Exercise information: Long axis strain rate)
Processing for acquiring the strain rate in the long axis direction of the heart wall as motion information will be described. This strain rate is information indicating the temporal change rate of the strain (strain) of the heart wall, and indicates the strain state of the heart wall.
When obtaining the strain rate, the same processing as in steps S41 to S48 in the flowchart of FIG. 13 is performed, and the strain δz12 in the major axis direction between the apex level and the papillary muscle level, the papillary muscle level, and the base level are determined. The strain δz23 in the major axis direction between the two is calculated.
Here, the strain δz12 and the strain δz23 are calculated for the apex short-axis image G1, the papillary muscle short-axis image G2, and the base-base short-axis image G3 in two time phases t1 and t2 (t1 ≠ t2). . The motion information calculation unit 62 calculates the strain rate in the long axis direction between the apex level and the papillary muscle level by dividing the strain δz12 by the time interval Δt = | t1-t2 |. Further, by dividing the strain δz23 by the time interval Δt, the strain rate in the major axis direction between the papillary muscle level and the base portion level is calculated. Note that the strain rate may be calculated from the strain by executing a normal differential operation.
The control unit 9 causes the display unit 81 to display the calculated heart wall strain rates δz12 / Δt and δz23 / Δt as exercise information indicating the temporal rate of change of the heart wall strain. The user can grasp the time change rate of the distortion of the heart wall by referring to the displayed strain rate of the heart wall.
It is also possible to calculate a strain rate for each of the intima and outer membranes of the heart wall and evaluate the time change rate of the strain based on these two strain rate values (for example, the average value of the two strain rate values). Etc.).
According to the ultrasonic diagnostic apparatus 1 that operates as described above, the following operations and effects are achieved.
First, the ultrasonic diagnostic apparatus 1 generates image data of moving images for two or more cross-sectional positions (for example, three cross-sectional positions of apex level, papillary muscle level, and cardiac base level) of a living tissue (heart). To do. The image data of the moving image is a series of tomographic image data (frames) in time series. In the case of displaying an image in real time, the moving image may not be displayed in consideration of the operation speed of the CPU or the like, but in that case, for example, it can be obtained by performing MPR processing on the volume data. It is possible to display only the tomographic image to be displayed.
Next, the ultrasound diagnostic apparatus 1 displays one tomographic image of a series of tomographic images for each of two or more cross-sectional positions. The user operates the operation unit 82 to designate a measurement image region for each displayed tomographic image. When the measurement image region is designated, the ultrasonic diagnostic apparatus 1 calculates a displacement along the time series of the measurement image region designated in the tomographic image at each cross-sectional position. Based on the displacement of the measurement image area calculated for each of the two or more cross-sectional positions, motion information indicating the motion state of the biological tissue is calculated.
Thus, unlike the conventional method in which the movement of the living tissue is measured based on the displacement (local movement information) of the measurement image region at one cross-sectional position, the ultrasonic diagnostic apparatus 1 according to this embodiment has two or more. Since the measurement image area is designated for each of the cross-sectional positions and the movement information is obtained by obtaining the displacement (local movement information) of the two or more measurement image areas, the three-dimensional movement of the living tissue is obtained. Measurement can be performed.
In particular, based on the displacement of the measurement image region at each cross-sectional position, a relative displacement (relative displacement (difference information); the above-described relative rotation angle, expansion and contraction, etc.) between different cross-sectional positions is obtained, and based on this relative displacement. Thus, it is possible to measure the three-dimensional movement of the living tissue. In addition, by performing measurement based on such relative displacement, it is not necessary to perform data analysis on the portion between the cross-sectional positions (that is, it is not necessary to perform three-dimensional data analysis). There is also the merit that it is possible to execute a typical measurement in a short time.
In addition, since the measurement image region specified in the two-dimensional tomographic image is three-dimensionally tracked based on the volume data, the three-dimensional displacement of the measurement image region can be obtained, so that highly accurate measurement can be performed. It can be carried out.
Note that the processing efficiency can be improved by switching between two-dimensional tracking and three-dimensional tracking according to the motion information to be acquired. It is also possible to configure an ultrasonic diagnostic apparatus that can execute only two-dimensional tracking or an ultrasonic diagnostic apparatus that can execute only three-dimensional tracking according to the motion information to be acquired.
Further, according to the ultrasonic diagnostic apparatus 1, unlike the conventional configuration in which a measurement image region is specified in a pseudo three-dimensional image to be displayed, a tomographic image (MPR image) based on volume data is displayed, Since the measurement image area is designated for the tomographic image, the measurement image area can be easily designated.
According to this ultrasonic diagnostic apparatus 1, in addition to the motion information described above, for example, a change (velocity) in the thickness of the heart wall, a strain or strain rate in the thickness direction (short axis direction) of the heart wall, a long axis Rotation angle (rotation speed) of the intima and epicardium of the heart wall with the direction as the axis, strain and strain rate of the heart wall in the rotation direction with the major axis as the axis, and heart in the rotation direction with the major axis as the axis It is also possible to acquire various (local) motion information such as the relative rotation angle between the inner membrane and outer membrane of the wall. These pieces of motion information can be obtained by two-dimensional tracking in one tomographic image as in the conventional case. Even when three-dimensional tracking is used, for example, motion information can be acquired from an image obtained by projecting volume data in a predetermined line-of-sight direction (an image obtained by rendering processing). Hereinafter, an example of a method for calculating these pieces of exercise information will be briefly described.
The change in the thickness (wall thickness) of the heart wall can be obtained by calculating the thickness of the heart wall and calculating the difference between the two tomographic images having different time phases. Here, the wall thickness of the heart wall is obtained at any position of the intima (or adventitia) by obtaining a straight line perpendicular to the tangent line that touches the intima (or adventitia) at that position, and this straight line is the adventitia (or The position (intersection point) that intersects the inner membrane) can be obtained, and the distance between the arbitrary position and the intersection point can be calculated. The wall thickness change rate can be calculated by dividing the change in wall thickness by the time between two tomographic images, or subjecting the change in wall thickness to normal differentiation (differential processing with time as a variable). This can be easily obtained.
The strain in the thickness direction of the heart wall obtains the change in the wall thickness as described above, and the change in the wall thickness is determined based on the wall thickness in one of the two tomographic images (a tomographic image at a certain time phase). It can be obtained by dividing by the value of. The strain rate can be obtained by dividing the strain value by the time interval between the time phases of the two tomographic images (or by differentiating with time).
The rotation angle of the intima (outer membrane) of the heart wall with the major axis as the axis is obtained for the position of the rotational direction in the minor axis direction of the intima (outer membrane) for two tomographic images having different time phases. The position of the intima (outer membrane) in the other tomographic image with respect to the position of the intima (outer membrane) in one tomographic image can be calculated. Also, the rotational speed can be easily obtained in the manner described above.
The strain of the heart wall in the rotation direction with the major axis as the axis calculates the distance in the rotation direction between the two positions of the intima (outer membrane) for two tomographic images having different time phases. Calculate the difference between Then, it can be obtained by dividing the difference value by the distance calculated for one of the two tomographic images. Further, the strain rate can be obtained by dividing the strain value by the time interval between the time phases of the two tomographic images.
The relative rotation angle between the intima and the epicardium of the heart wall in the rotation direction about the major axis direction is obtained as the intima rotation angle and the epicardial rotation angle in two tomographic images having different time phases (described above). ) And calculating the difference between them.
In this embodiment, the measurement image region parallel to the short axis direction of the heart is specified, but the measurement image region can also be specified to a cross section parallel to the long axis direction. In that case, the torsional motion and the relative rotational gradient in the major axis direction can be obtained by two-dimensional tracking in the tomographic image of the cross section. Further, the expansion / contraction, the strain, and the strain rate in the minor axis direction can be obtained by executing three-dimensional tracking using volume data.
Further, it is possible to designate a measurement image region parallel to an arbitrary cross section of the heart. In this case, the motion information that can be acquired only from the displacement in the direction parallel to the measurement image region can be obtained by two-dimensional tracking in the tomographic image in the direction parallel to the cross section. Further, motion information that requires displacement in a direction orthogonal to the cross section can be obtained by executing three-dimensional tracking using volume data.
The same motion information acquisition process can also be performed when a three-dimensional ultrasonic scan is performed in an arbitrary scan mode for generating volume data. For example, even when a two-dimensional scan surface is rotated and scanned, the same processing as in this embodiment can be executed. That is, as long as volume data of a living tissue can be generated, an ultrasonic scan mode can be arbitrarily selected.
The ultrasonic diagnostic apparatus 1 described in detail above is only a specific example for suitably implementing the ultrasonic diagnostic apparatus according to the present invention. Hereinafter, various modifications related to the ultrasonic diagnostic apparatus according to the present invention will be described.
This modification is intended to facilitate the designation operation of the measurement image region for the image based on the volume data obtained by the three-dimensional ultrasonic scan for the living tissue. The ultrasonic diagnostic apparatus according to this modification has the same configuration as the ultrasonic diagnostic apparatus 1 according to the above embodiment. Hereinafter, processing according to this modification will be described.
Processing until the MPR processing unit 52 of the image processing unit 5 generates image data of a series of tomographic images along a time series is the same as that in the above embodiment. In the above embodiment, the image data of the tomographic image is generated for each of two or more cross-sectional positions. However, in this modification, the image data of the tomographic image may be generated for only one cross-sectional position. .
The control unit 9 causes the display unit 81 to display one tomographic image of the series of tomographic images. The user operates the operation unit 82 to designate a measurement image region for the displayed tomographic image. This makes it possible to easily specify the measurement image region as compared with the conventional configuration in which the pseudo three-dimensional image is displayed and the measurement image region is specified.
According to this modification, for example, when a measurement image region parallel to the short-axis cross section of the heart is designated, for example, the change (velocity) of the thickness of the heart wall in the motion information described above, the thickness of the heart wall Strain and strain rate in the direction, rotation angle (rotation speed) of the intima and epicardium of the heart wall with the major axis as the axis, strain and strain rate of the heart wall in the rotation direction with the major axis as the axis, and major axis It is possible to calculate the relative rotation angle between the intima and the epicardium of the heart wall in the rotation direction about the direction.
In the ultrasonic diagnostic apparatus 1 according to the above-described embodiment, the user operates the operation unit 82 to designate a measurement image region. In this modification, a measurement image region designation operation is automatically performed.
An example of an ultrasonic diagnostic apparatus according to this modification is shown in FIG. The ultrasonic diagnostic apparatus 100 shown in the figure has substantially the same configuration as the ultrasonic diagnostic apparatus 1 of the above embodiment, but is different in that the image processing unit 5 includes a measurement image region specifying unit 53. Further, in order to execute processing characteristic of this modification, a medical image processing program 72 different from that of the above embodiment is provided.
The measurement image region designation unit 53 analyzes the image data of the tomographic image (for example, apex short axis image, papillary muscle short axis image, base base short axis image) generated by the MPR processing unit 52, and measures this tomographic image. Specifies the image area.
More specifically, the measurement image region designating unit 53 applies a general boundary extraction process, for example, and analyzes the pixel value of each pixel of the image data to thereby obtain a boundary portion (cardiac wall) of the living tissue. Extract the inner membrane, outer membrane, etc.). Then, the measurement image region specifying unit 53 specifies the extracted boundary portion as the measurement image region in the tomographic image. The measurement image area designating unit 53 corresponds to an example of the “designating unit” of the present invention, and includes, for example, a microprocessor that operates based on the medical image processing program 72.
The process performed by the measurement image area specifying unit 53 corresponds to the procedure of step S09 in the flowchart of FIG. Note that when the measurement image region specifying unit 53 performs the measurement image region specifying process, the tomographic image may or may not be displayed on the display unit 81.
In addition, it is desirable to display the tomographic image in which the measurement image region is input on the display unit 81 so that the user can confirm the measurement image region automatically designated. At this time, it is desirable to configure so that the measurement image region can be appropriately changed using the operation unit 82.
According to this modification, since the measurement image area is automatically specified, it is easy to specify the measurement image area (in fact, when the automatically specified measurement image area is used as it is, the measurement image area is used as it is. (There is no need to specify the area.)
In the ultrasonic diagnostic apparatus 1 of the above embodiment, a three-dimensional ultrasonic scan is electronically performed by using a two-dimensional ultrasonic transducer in which ultrasonic transducers are two-dimensionally arranged. This modification relates to an ultrasonic diagnostic apparatus equipped with a one-dimensional ultrasonic probe in which ultrasonic transducers are arranged in a one-dimensional direction.
When a one-dimensional ultrasonic probe is used, only an ultrasonic scan in the one-dimensional direction (main scanning direction X in FIGS. 2 and 3) can be electronically executed. Scanning is performed manually or mechanically.
Even when such a one-dimensional ultrasonic probe is used, volume data based on a three-dimensional ultrasonic scan is generated, image data of a tomographic image is generated based on the volume data, and the tomographic image is displayed. The measurement image area can be specified.
The ultrasonic diagnostic apparatus 1 of the above embodiment calculates the displacement along the time series of the measurement image region specified in the tomographic image at two or more cross-sectional positions such as the apex level and the papillary muscle level, and the two or more The motion information of the living tissue is calculated based on the displacement of the measurement image region at the cross-sectional position. That is, in the above-described embodiment, motion information is calculated based only on the displacement of the designated measurement image region.
In this modified example, a measurement image area other than the designated measurement image area is automatically designated separately, and motion information is calculated in consideration of the displacement of the automatically designated measurement image area. Hereinafter, an ultrasonic diagnostic apparatus according to this modification will be described. Similar to the ultrasonic diagnostic apparatus 100 shown in FIG. 15, the ultrasonic diagnostic apparatus of this modification includes a measurement image region specifying unit 53.
The operation of the ultrasonic diagnostic apparatus according to this modification will be described. As in the above-described embodiment, the user can use, for example, an intima position image m1, a measurement image region of the apex short axis image G1, the papillary muscle short axis image G2, and the base basal short axis image G3 displayed on the display unit 81. m2 and m3 and outer membrane position images M1, M2 and M3 are designated as measurement image regions, respectively (see FIGS. 5 and 6). Note that the measurement image area may be designated by the automatic designation process of [Modification 2].
The measurement image area designation unit 53 designates measurement image areas other than the designated measurement image areas m1 to m3 and M1 to M3. As an example, the measurement image region designation unit 53 designates the intima position image m4 and the epicardium position image M4 at positions outside the apex short axis image G1, as shown in FIG. The intima position image m5 and the epicardial position image M5 are designated at a position between the image G1 and the papillary muscle short axis image G2, and the intima position is set between the papillary muscle short axis image G2 and the basal base short axis image G3. An image m6 and an outer membrane position image M6 are designated.
This measurement image area automatic designation processing will be described more specifically. First, the measurement image region specifying unit 53 is based on the cross-sectional positions (coordinates) of the apex short-axis image G1, the papillary muscle short-axis image G2, and the basal-portion short-axis image G3 in which the measurement image region is specified by the user or the like. A cross-sectional position (coordinates) for designating a new measurement image area is determined. This processing is performed, for example, by obtaining coordinates (first Z coordinates) at a position away from the apex short axis image G1 by a predetermined distance in the direction opposite to the papillary muscle short axis image G2, and then calculating the Z coordinates of the apex short axis image G1. And the Z coordinate of the papillary muscle short axis image G2 and the Z coordinate of the basal part short axis image G3. The center coordinate (third Z coordinate) is obtained.
Next, the MPR processing unit 52 generates image data of tomographic images in directions parallel to the apex short axis image G1 and the like at the first to third Z coordinates based on the volume data.
The measurement image region designating unit 53 extracts the boundary portion of the heart wall by analyzing the pixel value of the image data of the tomographic image at the first Z coordinate, sets the inner boundary portion as the intima position image m4, and The boundary portion of is an outer membrane position image M4. Similarly, the boundary portion of the heart wall is extracted by analyzing the pixel value of the image data of the tomographic image at the second Z coordinate, the inner boundary portion is set as the intima position image m5, and the outer boundary portion is set as the outer boundary portion. The boundary position of the heart wall is extracted by analyzing the pixel value of the image data of the tomographic image at the third Z coordinate as the membrane position image M5, the inner boundary portion is set as the intima position image m6, and the outer boundary The portion is referred to as an outer membrane position image M6. Then, the intima position images m4 to m6 and the epicardial position images M4 to M6 are designated as new measurement image areas, respectively.
The displacement calculation unit 61 of this modification includes all of the measurement image areas m1 to m3 and M1 to M3 designated by the user and the measurement image areas m4 to m6 and M4 to M6 designated by the measurement image area designation unit 53. For each, the displacement along the time series is calculated.
The exercise information calculation unit 62 calculates exercise information based on all the displacements of the measurement image regions m1 to m6 and M1 to M6. For example, when evaluating the torsional motion of the heart wall, the relative rotation angle between the intima position image m4 and the intima position image m1 is calculated, and the relative rotation angle between the intima position image m1 and the intima position image m5 is calculated. Then, the relative rotation angle between the intima position image m5 and the intima position image m2 is calculated, the relative rotation angle between the intima position image m2 and the intima position image m6 is calculated, and the intima position image m6 and the intima A relative rotation angle with the position image m3 is calculated. The relative rotation angle can be similarly calculated for the outer membrane position images M1 to M6.
The control unit 9 causes the display unit 81 to display motion information based on the calculated relative rotation angle. Thus, according to this modification, it is possible to acquire exercise information with higher accuracy than in the above embodiment.
In the above embodiment, the user designates the cross-sectional position for designating the measurement image area (see the cross-section position designation cursors C1 to C3 in FIG. 5), but the cross-sectional position is automatically designated. It is also possible to do.
For example, when there is image data of ultrasonic images acquired for the same living tissue in the past, such as preoperative and postoperative observations and follow-up observations, by storing the designated cross-sectional position, the past The information of the cross-sectional position specified in (1) can be read and automatically specified as the current cross-sectional position.
In addition, typical cross-sectional positions such as apex level, papillary muscle level, and cardiac base level are set in advance, and for example, a typical cross-sectional position is determined based on image data or volume data of a B-mode image. Can do.
In addition, based on the image data of the B-mode image, etc., the size of the living tissue (for example, the length in the long axis direction of the heart) is analyzed, and each cross-sectional position when the size is divided into a plurality of parts is measured. It is also possible to configure so that the image area is designated as a cross-sectional position.
In the above embodiment, the displacement of the living tissue is obtained and the displacement is differentiated (divided by the time interval) to obtain the velocity. First, the velocity is obtained, and the velocity is integrated to obtain the displacement. You may make it ask.
In the above embodiment, the motion information indicating the (more global) motion state of the living tissue is obtained and displayed, but the local motion information in each measurement image region is displayed. It may be configured. In this case, when the motion information is not calculated and displayed, it is possible to apply an ultrasonic diagnostic apparatus having only a configuration necessary for calculating and displaying the local motion information (this will be described later). The same applies to the medical image processing apparatus and the medical image processing program.
A medical image processing apparatus according to the present invention will be described. This medical image processing apparatus includes, for example, a computer connected to an ultrasonic diagnostic apparatus and a computer connected to a database such as PACS (Picture Archiving and Communication System) that stores image data of ultrasonic images. The FIG. 17 shows an example of a medical image processing apparatus according to the present invention. In the figure, the same components as those in the first embodiment are denoted by the same reference numerals.
A medical image processing apparatus 1000 illustrated in FIG. 17 includes an image processing unit 5, an arithmetic processing unit 6, a storage unit 7, a user interface 8, and a control unit 9 similar to those of the ultrasonic diagnostic apparatus 1 illustrated in FIG. 1. The storage unit 7 functions as an example of the “storage unit” of the present invention. The medical image processing apparatus 1000 is connected to the ultrasonic diagnostic apparatus 2000 and the medical image database 3000 via a network N such as a LAN (Local Area Network). The control unit 9 of this embodiment includes a network adapter that performs data communication through the network N.
The operation of the medical image processing apparatus 1000 will be described. The medical image processing apparatus 1000 receives ultrasonic image data from the ultrasonic diagnostic apparatus 2000 and the medical image database 3000.
When the input image data is image data of a B-mode image, the volume data generation unit 51 generates volume data based on this image data. Similar to the first embodiment, the MPR processing unit 52 generates image data of a series of tomographic images along a time series based on the volume data. The series of tomographic image data is stored in the storage unit 7.
When the input image data is volume data, the MPR processing unit 52 generates a series of tomographic image data along the time series based on the volume data. The series of tomographic image data is stored in the storage unit 7.
The medical image processing apparatus 1000 executes processing similar to that in the first embodiment based on a series of tomographic image data (and volume data) stored in the storage unit 7 (see FIGS. 4 to 14). ). Thereby, the measurement of the three-dimensional movement of the living tissue can be performed in a short time. In addition, it is possible to easily specify a region for measuring the movement of the living tissue.
As a modification of this embodiment, a configuration in which the measurement image region specifying unit 53 is provided in the image processing unit 5 can be applied (see FIG. 15, processing based on the medical image processing program 72). Thereby, the same processing as that of the second modification and the fourth modification of the first embodiment can be executed.
A medical image processing program according to the present invention will be described. The medical image processing programs 71 and 72 described in the first and second embodiments are examples of the medical image processing program according to the present invention. The medical image processing program causes a computer to execute the processing described in the first embodiment and its modifications. This medical image processing program may be stored in advance in a storage device such as a hard disk drive built in the computer, or stored in advance in a server on a network such as a LAN, and the computer reads it out. May be executed.
This medical image processing program can be stored in an arbitrary storage medium so that it can be read by a computer. Examples of the storage medium include optical disks, magneto-optical disks (CD-ROM / DVD-RAM / DVD-ROM / MO, etc.), magnetic storage media (hard disk / floppy (registered trademark) disk / ZIP, etc.), semiconductor memory, and the like. There is.
1 is a schematic block diagram showing an example of the overall configuration of a preferred embodiment of an ultrasonic diagnostic apparatus according to the present invention. It is a schematic explanatory drawing for demonstrating an example of the scanning aspect of the ultrasonic wave in suitable embodiment of the ultrasonic diagnosing device concerning this invention. It is a schematic explanatory drawing for demonstrating an example of the scanning aspect of the ultrasonic wave in suitable embodiment of the ultrasonic diagnosing device concerning this invention. It is a flowchart showing an example of the operation | movement aspect of suitable embodiment of the ultrasound diagnosing device based on this invention. It is the schematic showing an example of the display mode of the display screen by suitable embodiment of the ultrasound diagnosing device concerning this invention. It is the schematic showing an example of the display mode of the display screen by suitable embodiment of the ultrasound diagnosing device concerning this invention. It is a flowchart showing an example of the process by suitable embodiment of the ultrasound diagnosing device concerning this invention. It is a schematic explanatory drawing for demonstrating the process by suitable embodiment of the ultrasound diagnosing device concerning this invention. It is a flowchart showing an example of the process by suitable embodiment of the ultrasound diagnosing device concerning this invention. It is a schematic explanatory drawing for demonstrating the process by suitable embodiment of the ultrasound diagnosing device concerning this invention. It is a flowchart showing an example of the process by suitable embodiment of the ultrasound diagnosing device concerning this invention. It is a schematic explanatory drawing for demonstrating the process by suitable embodiment of the ultrasound diagnosing device concerning this invention. It is a flowchart showing an example of the process by suitable embodiment of the ultrasound diagnosing device concerning this invention. It is a schematic explanatory drawing for demonstrating the process by suitable embodiment of the ultrasound diagnosing device concerning this invention. It is a schematic block diagram showing an example of the whole structure of the modification of suitable embodiment of the ultrasound diagnosing device concerning this invention. It is a schematic explanatory drawing for demonstrating the process by the modification of suitable embodiment of the ultrasound diagnosing device concerning this invention. 1 is a schematic block diagram showing an example of the overall configuration of a preferred embodiment of an ultrasonic diagnostic apparatus according to the present invention.
DESCRIPTION OF SYMBOLS 1 Ultrasonic diagnostic apparatus 2 Two-dimensional ultrasonic probe 3 Transmission / reception part 4 Signal processing part 41 B mode processing part 5 Image processing part 51 Volume data generation part 52 MPR processing part 6 Operation processing part 61 Displacement calculation part 62 Motion information calculation part 7 Storage unit 71, 72 Medical image processing program 8 User interface 81 Display unit 82 Operation unit 9 Control unit P 2D scan plane R 3D scan area X Main scan direction Y Subscan direction m1 to m6 Intima position images M1 to M6 Outside Membrane Position Image 1000 Medical Image Processing Device
Transmitting and receiving means for transmitting and receiving ultrasonic waves to the ultrasonic probe;
Based on the reception signal obtained as a result of the transmission and reception of the ultrasonic wave, an image generation means for generating image data of a series of tomographic images along a time series for each of two or more cross-sectional positions of the living tissue;
Display means for displaying one tomographic image of the series of tomographic images based on the generated image data for each of the two or more cross-sectional positions;
For each of the two or more cross-sectional positions, designation means for designating a measurement image region in the displayed one tomographic image;
For each of the two or more cross-sectional positions, based on the image data of the series of tomographic images, local motion information indicating the state of motion of the living tissue in the designated measurement image region is calculated, and the two or more cross-sectional positions are calculated. Based on the calculated local motion information for each of the cross-sectional positions, calculating means for calculating motion information indicating the state of motion of the living tissue;
The display means displays the exercise information calculated by the calculation means.
The ultrasonic probe scans the transmission / reception direction of ultrasonic waves along a direction orthogonal to the predetermined direction at each of a plurality of positions along the predetermined direction by the transmission / reception means, Repeating transmission and reception of ultrasonic waves along the orthogonal direction,
An image of a first tomographic image in a cross section including the transmission / reception direction of the living tissue and the orthogonal direction based on reception signals sequentially obtained when the ultrasonic transmission / reception direction is scanned in the orthogonal direction. First tomographic image generation means for generating data;
Volume data generation means for generating volume data of the living tissue based on the generated image data of the first tomographic image for each of the plurality of positions;
Second tomographic image generation means for generating image data of a second tomographic image at each of the two or more cross-sectional positions based on the generated volume data;
For each of the two or more cross-sectional positions, a plurality of image data of the second tomographic image obtained by repeating the transmission / reception of the ultrasonic wave, image data of a series of tomographic images along the time series, and To
The calculation means calculates the difference information indicating the difference of the local movement information at any two of the two or more cross-sectional positions to obtain the movement information.
The cross section of the series of tomographic images at one of the two cross sectional positions and the cross section of the series of tomographic images at the other cross sectional position are parallel to each other,
For each of the two cross-sectional positions, calculating the rotation angle of the biological tissue in the measurement image region with the direction orthogonal to the cross-section of the series of tomographic images as the local motion information,
The difference between the calculated rotation angles for each of the two cross-sectional positions is calculated as the difference information, and the calculated difference between the rotation angles is used as the motion information indicating the torsional motion of the living tissue.
The ultrasonic diagnostic apparatus according to claim 3.
The difference between the calculated rotation angles for each of the two cross-sectional positions is calculated as the difference information, and the calculated rotation angle difference is differentiated with respect to time, whereby the speed of the torsional motion of the living tissue is determined. Calculating the motion information shown,
The difference between the calculated rotation angles for each of the two cross-sectional positions is calculated as the difference information, and the calculated rotation angle difference is divided by the distance between the two cross-sectional positions. Obtaining the motion information indicating the relative rotational gradient of the tissue;
For each of the two or more cross-sectional positions, based on the volume data, calculating the local motion information indicating a three-dimensional motion state of the biological tissue in the measurement image region specified by the specifying means,
The motion information is calculated based on local motion information indicating the calculated three-dimensional motion state for each of the two or more cross-sectional positions.
The calculation means calculates difference information indicating a difference in local movement information indicating a state of the three-dimensional movement at any two of the two or more cross-section positions to obtain the movement information.
The ultrasonic diagnostic apparatus according to claim 7.
For each of the two cross-sectional positions, based on the calculated local motion information indicating the state of the three-dimensional motion, in a direction orthogonal to the cross section of the series of tomographic images of the living tissue in the measurement image region. And calculating the difference in the calculated displacement in the orthogonal direction as the difference information for each of the two cross-sectional positions, and calculating the calculated difference in the orthogonality of the living tissue. The motion information indicating the stretching motion in the direction to perform,
The ultrasonic diagnostic apparatus according to claim 8.
For each of the two cross-sectional positions, based on the calculated local motion information indicating the state of the three-dimensional motion, in a direction orthogonal to the cross section of the series of tomographic images of the living tissue in the measurement image region. And calculating the difference of the calculated displacement in the orthogonal direction as the difference information for each of the two cross-sectional positions, and differentiating the calculated difference of the time by time. Calculating the motion information indicating the speed of expansion and contraction motion in the orthogonal direction of the biological tissue,
For each of the two cross-sectional positions, based on the calculated local motion information indicating the state of the three-dimensional motion, in a direction orthogonal to the cross section of the series of tomographic images of the living tissue in the measurement image region. And calculating the difference between the calculated displacements in the orthogonal direction as the difference information for each of the two cross-sectional positions, and calculating the difference between the calculated displacements of the two cross-sectional positions. Dividing by the original distance between, the quotient is the motion information indicating the strain in the orthogonal direction of the biological tissue,
For each of the two cross-sectional positions, based on the calculated local motion information indicating the state of the three-dimensional motion, in a direction orthogonal to the cross section of the series of tomographic images of the living tissue in the measurement image region. And calculating the difference between the calculated displacements in the orthogonal direction as the difference information for each of the two cross-sectional positions, and calculating the difference between the calculated displacements of the two cross-sectional positions. Dividing by the original distance between and calculating the motion information indicating the strain rate in the orthogonal direction of the biological tissue by differentiating the quotient by time,
The designation unit designates a new measurement image region at a cross-sectional position different from the two or more cross-sectional positions based on the designated measurement image region with respect to the one tomographic image displayed on the display unit. Including automatic designation means to
The computing means computes local motion information indicating a motion state of the living tissue in the designated new measurement image region, and the computation is performed for each of the two or more cross-sectional positions and the different cross-sectional positions. Calculating the motion information based on the local motion information;
The measurement image area specified by the specifying means is an image area corresponding to the heart wall of the heart in the one tomographic image displayed on the display means.
The ultrasonic diagnostic apparatus according to any one of claims 1 to 13, wherein the apparatus is an ultrasonic diagnostic apparatus .
A medical image processing apparatus for processing image data of a medical image of a biological tissue obtained by an ultrasonic diagnostic apparatus,
Storage means for storing image data of a series of tomographic images along a time series at each of two or more cross-sectional positions of a biological tissue;
Display means for displaying one tomographic image of the series of tomographic images based on the stored image data for each of the two or more cross-sectional positions;
A medical image processing apparatus .
A computer having storage means for storing image data of a series of tomographic images along a time series at each of two or more cross-sectional positions of a biological tissue, and a display means,
For each of the two or more cross-sectional positions, based on the stored image data, function to display one tomographic image of the series of tomographic images on the display means,
For each of the two or more cross-sectional positions, in response to the measurement area being designated in the displayed one tomographic image, based on the image data of the series of tomographic images, in the designated measurement image area Function to calculate local motion information indicating the state of motion of the biological tissue,
Based on the calculated local motion information for each of the two or more cross-sectional positions, function to calculate motion information indicating the motion state of the biological tissue,
Causing the display means to display the calculated exercise information;
A medical image processing program characterized by that .
JP2006149414A 2006-05-30 2006-05-30 Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing program Active JP4745133B2 (en)
JP2006149414A JP4745133B2 (en) 2006-05-30 2006-05-30 Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing program
PCT/JP2007/000587 WO2007138751A1 (en) 2006-05-30 2007-05-30 Ultrasonograph, medical image processing device, and medical image processing program
US12/302,506 US8343052B2 (en) 2006-05-30 2007-05-30 Ultrasonograph, medical image processing device, and medical image processing program
CN 200780024890 CN101484074B (en) 2006-05-30 2007-05-30 Ultrasonograph, and medical image processing apparatus and program
EP20070737243 EP2022404B1 (en) 2006-05-30 2007-05-30 Ultrasonograph, medical image processing device, and medical image processing program
JP2007319190A JP2007319190A (en) 2007-12-13
JP4745133B2 true JP4745133B2 (en) 2011-08-10
ID=38778278
JP2006149414A Active JP4745133B2 (en) 2006-05-30 2006-05-30 Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing program
US (1) US8343052B2 (en)
EP (1) EP2022404B1 (en)
JP (1) JP4745133B2 (en)
CN (1) CN101484074B (en)
WO (1) WO2007138751A1 (en)
JP5624258B2 (en) 2007-04-26 2014-11-12 株式会社東芝 Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image processing program
JP5414157B2 (en) 2007-06-06 2014-02-12 株式会社東芝 Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image processing program
JP5468909B2 (en) * 2008-01-31 2014-04-09 株式会社日立メディコ Medical diagnostic imaging equipment
JP5259267B2 (en) 2008-06-19 2013-08-07 株式会社東芝 Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image processing program
US8538103B2 (en) 2009-02-10 2013-09-17 Hitachi Medical Corporation Medical image processing device, medical image processing method, medical image diagnostic apparatus, operation method of medical image diagnostic apparatus, and medical image display method
WO2010092919A1 (en) * 2009-02-13 2010-08-19 株式会社 日立メディコ Medical image display method, medical image diagnosis device, and medical image display device
JP5658296B2 (en) * 2013-03-11 2015-01-21 株式会社東芝 Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image processing program
KR101590876B1 (en) 2014-02-21 2016-02-02 삼성전자주식회사 Method and apparatus for smoothing motion vector
CN106308848A (en) * 2015-07-10 2017-01-11 通用电气公司 Method and device for measuring ultrasonic image
KR20190085342A (en) * 2018-01-10 2019-07-18 삼성메디슨 주식회사 Method for controlling ultrasound imaging apparatus and ultrasound imaging aparatus thereof
JP3219280B2 (en) 1991-02-27 2001-10-15 株式会社日立メディコ Spect tomography apparatus
JP3382831B2 (en) * 1997-11-11 2003-03-04 ジーイー横河メディカルシステム株式会社 Method of manufacturing ultrasonic transducer array, ultrasonic transducer array, ultrasonic probe, and ultrasonic imaging apparatus
JP4116122B2 (en) * 1997-11-28 2008-07-09 株式会社東芝 Ultrasonic diagnostic apparatus and ultrasonic image processing apparatus
JP2003319939A (en) * 2002-04-26 2003-11-11 Ge Medical Systems Global Technology Co Llc Ultrasonic imaging device
JP2004313291A (en) * 2003-04-14 2004-11-11 Toshiba Corp Ultrasonograph, and medical image analysis instrument and method
JP2006061581A (en) 2004-08-30 2006-03-09 Toshiba Corp Medical motion analysis apparatus, medical motion analysis method and medical motion analysis program
EP1744181B8 (en) 2005-07-15 2014-04-09 Kabushiki Kaisha Toshiba Ultrasonic diagnostic apparatus and ultrasonic image processing method
2006-05-30 JP JP2006149414A patent/JP4745133B2/en active Active
2007-05-30 US US12/302,506 patent/US8343052B2/en active Active
2007-05-30 CN CN 200780024890 patent/CN101484074B/en active IP Right Grant
2007-05-30 WO PCT/JP2007/000587 patent/WO2007138751A1/en active Application Filing
2007-05-30 EP EP20070737243 patent/EP2022404B1/en active Active
US20090198133A1 (en) 2009-08-06
WO2007138751A1 (en) 2007-12-06
EP2022404A1 (en) 2009-02-11
US8343052B2 (en) 2013-01-01
CN101484074A (en) 2009-07-15
JP2007319190A (en) 2007-12-13
CN101484074B (en) 2010-12-22
EP2022404A4 (en) 2010-08-18
EP2022404B1 (en) 2011-09-07
US20050049502A1 (en) 2005-03-03 Method and apparatus for correcting a volumetric scan of an object moving at an uneven period
2009-02-14 RD02 Notification of acceptance of power of attorney
2011-02-08 A871 Explanation of circumstances concerning accelerated examination
2011-02-25 A975 Report on accelerated examination