Ultrasonic imaging apparatus indicating dispersion of characteristic parameter of medium.

BACKGROUND OF THE INVENTION

(1) Field of the Invention

0001The present invention relates to an ultrasonic imaging apparatus which functions to display an image of a two-dimensional distribution of a tissue characterization parameter of a tissue, for example, an ultrasonic attenuation coefficient, in a realtime mode.

0002In the field of the ultrasonic diagnosis, the two-dimensional distribution of the tissue characterization parameters, for example, an ultrasonic attenuation coefficient, an ultrasonic reflection coefficient, an elasticity parameter, are used to diagnose a condition of the tissue. For example, the value of the ultrasonic attenuation coefficient of tissue of a liver suffering from hepatitis is small compared with a normal tissue, and the value of the ultrasonic attenuation coefficient of a fatty liver is large compared with a normal tissue. Since high accuracy is required in ultrasonic diagnosis, the tissue characterization parameters must be determined precisely.

(2) Description of the Related Art

0003Conventionally, the ultrasonic imaging apparatus which functions to display an image of a two-dimensional distribution of a tissue characterization parameter of a tissue in a realtime mode is disclosed in the Japanese Unexamined Patent Publication No.62-109553, and the U.S.Patent No.4,836,210 of Shiba et al, dated June 6, 1989.

0004In ultrasonic imaging apparatuses, a pulsed beam of an ultrasound is generated to inject the beam into a tissue in a certain direction, and ultrasound which is reflected at various depths of the tissue in the same direction are detected. The direction of the pulsed ultrasound is scanned so that a plane is scanned by the pulsed ultrasound during one plane scan cycle. Thus, a two-dimensional distribution of the tissue characterization parameter on the plane (cross section of the tissue) is obtained. The above injection and detection of the ultrasound is carried out through an ultrasonic transducer array. The operator (doctor) applies a probe containing the transducer array on a skin of a human body to diagnose the condition of the tissue under the skin.

0005However, the obtained values of the above tissue characterization parameters fluctuate due to subtle movements of the probe or the tissue which is subject to the diagnosis. Further, for example, the values of the ultrasonic attenuation coefficient, which are obtained as above, greatly vary when a structure which is different from the tissue which is to be diagnosed, e.g., a vascular tract, exists in the scanned cross section. Generally, it is considered that the obtained values of the tissue characterization parameters are significant only when the obtained values are stable. Thus, generally, the values of the tissue characterization parameters must be obtained on a cross section of the tissue where the obtained values are stable. Nevertheless, conventionally, it is difficult for the operator to determine whether or not an obtained value of a tissue characterization parameter is stable. Therefore, in conventional ultrasonic imaging apparatuses, the credibility of the obtained values of the tissue characterization parameters is low.

SUMMARY OF THE INVENTION

0006An object of the present invention is to provide an ultrasonic imaging apparatus whereby an operator can easily determine stability of a value of a tissue characterization parameter which is obtained from the apparatus.

0007Another object of the present invention is to provide an ultrasonic imaging apparatus whereby an operator can easily determine a cross section which gives the most stable value of a tissue characterization parameter which is to be obtained, among a plurality of cross sections from which values of the tissue characterization parameter are most recently obtained.

0008According to the aspect of the present invention, there is provided an ultrasonic imaging apparatus which comprises: an ultrasound generating unit for generating a pulsed ultrasound to externally output in a certain direction for each output of the pulsed ultrasound; an ultrasound detecting unit for detecting reflected ultrasound which is generated by reflection of the pulsed ultrasound at various points outside of the ultrasonic imaging apparatus in the above direction, and reaches the ultrasound detecting unit after the above output of pulsed ultrasound; a scanning unit for scanning the direction of the pulsed ultrasound which is generated by the ultrasound generating unit and the reflected ultrasound which reaches the ultrasound detecting unit, so that a plane is scanned by the pulsed ultrasound which is generated by the ultrasound generating unit during one plane scan cycle; a two-dimensional distribution obtaining unit for obtaining a two-dimensional distribution of a predetermined characteristic parameter of a medium of the ultrasound at various points in the plane, wherein distances from the ultrasound detecting unit to the various points in the plane at which the reflected ultrasound are generated, are determined by measuring times which elapsed after the pulsed ultrasound is output from the ultrasound generating unit until the corresponding reflected ultrasound are detected in the ultrasound detecting unit; a region of interest designating unit for designating a region of interest in the plane; and a statistic information obtaining unit for obtaining statistic information on a dispersing characteristic of the predetermined characteristic parameter of the medium at the various points in the region of interest of the plane.

0009The ultrasonic imaging apparatus further comprises a statistic information indicating unit for indicating the statistic information.

0010The ultrasonic imaging apparatus further comprises a statistic information storing unit for storing the statistic information for a predetermined number of planes.

0011The ultrasonic imaging apparatus further comprises a stored statistic information indicating unit for indicating the statistic information which is stored in the statistic information storing unit.

0012The ultrasonic imaging apparatus further comprises a two-dimensional distribution storing unit for storing the two-dimensional distributions of the predetermined characteristic parameter of the medium for a predetermined number of planes.

0013The ultrasonic imaging apparatus further comprises a stored distribution indicating unit for indicating the two-dimensional distributions which are stored in the two-dimensional distribution storing unit.

0014The statistic information storing unit stores the statistic information for a predetermined number of planes for which the two-dimensional distribution of the predetermined characteristic parameter of the medium are most recently obtained by the two-dimensional distribution obtaining unit.

0015The ultrasonic imaging apparatus further comprises a stored statistic information indicating unit for indicating the statistic information which is stored in the statistic information storing unit.

0016The ultrasonic imaging apparatus further comprises a two-dimensional distribution storing unit for storing the two-dimensional distributions of the predetermined characteristic parameter of the medium for a predetermined number of planes which are most recently obtained by the two-dimensional distribution obtaining unit.

0017The ultrasonic imaging apparatus further comprises a stored distribution indicating unit for indicating one of the two-dimensional distributions which are stored in the two-dimensional distribution storing unit.

0018The ultrasonic imaging apparatus further comprises a storage freezing unit for freezing an operation renewing the statistic information storing unit and the two-dimensional distribution storing unit.

0019The ultrasonic imaging apparatus further comprises a frozen plane designating unit for designating one of the planes for which the statistic information and the two-dimensional distributions are stored in the statistic information storing unit and the two-dimensional distribution storing unit, respectively, after the statistical information storing means and the two-dimensional distribution storing means are frozen.

0020The ultrasonic imaging apparatus further comprises a designated information indicating unit for indicating the statistic information and the two-dimensional distribution for the plane which are designated by the frozen plane designating unit.

0021The ultrasonic imaging apparatus further comprises: a plane designating unit for designating one of the planes for which the statistic information and the two-dimensional distributions are stored in the statistic information storing unit and the two-dimensional distribution storing unit, and a history information indicating unit for indicating a recent history of the statistic information from a plane for which the two-dimensional distribution is most recently obtained, to the plane which is designated by the plane designating unit.

0022When the statistic information obtaining unit comprises a standard deviation obtaining unit for obtaining a standard deviation of the predetermined characteristic parameter of the medium at the various points in the region of interest of the plane, the ultrasonic imaging apparatus may further comprise: a minimum standard deviation plane determining unit for determining a plane for which the minimum standard deviation is obtained by the standard deviation obtaining unit, among planes for which the statistic information and the two-dimensional distributions are stored in the statistic information storing unit and the two-dimensional distribution storing unit; and a minimum standard deviation plane indicating unit for indicating the plane for which the minimum standard deviation is obtained by the standard deviation obtaining unit.

0023The ultrasonic imaging apparatus further comprises an optimum two-dimensional distribution indicating unit for indicating a two-dimensional distribution for the plane for which the minimum standard deviation is obtained by the standard deviation obtaining unit.

0024The statistic information includes an average of the predetermined characteristic parameter of the medium at the various points in the region of interest of the plane.

0025The statistic information includes a standard deviation of the predetermined characteristic parameter of the medium at the various points in the region of interest of the plane.

0026The statistic information includes a skewness of a frequency distribution of values of the tissue characterization parameter which was obtained in the parameter calculation circuit within the region of interest of the predetermined characteristic parameter of the medium at the various points in the region of interest of the plane.

0027The statistic information includes a kurtosis of a frequency distribution of the values of the tissue characterization parameter of the medium at the various points in the region of interest of the plane.

0028The ultrasonic imaging apparatus further comprises: a region of interest redesignating unit for designating a region of interest other than the region of interest which is designated by the region of interest designating unit, in the plane; and a statistic information reobtaining unit for obtaining statistic information on a dispersing characteristic of the predetermined characteristic parameter of the medium at the various points in the region of interest which is designated by the region of interest redesignating unit.

BRIEF DESCRIPTION OF THE DRAWINGS

0029In the drawings: <ul id="ul0001" list-style="none"><li>Figure 1 is a block diagram showing a construction of the embodiment of the present invention;</li><li>Figure 2 is a block diagram showing a construction of an example of the parameter calculation circuit 100 in Fig. 1;</li><li>Figure 3 is a block diagram showing a construction of the characteristic portion of the embodiment of the present invention;</li><li>Figure 4 is a block diagram showing a construction of the frequency distribution obtaining circuit 33 in Fig. 3;</li><li>Figure 5 is a histogram showing a frequency distribution of values of a tissue characterization parameter in a region of interest;</li><li>Figure 6 is a diagram showing an example of the display of the two-dimensional distribution of detected values of a tissue characterization parameter and an indication of the designated region of interest;</li><li>Figure 7 is a diagram showing an example indication of a recent history of the average of detected values of a tissue characterization parameter in a region of interest;</li><li>Figure 8 is a diagram showing an example indication of a recent history of the standard deviation of detected values of a tissue characterization parameter in a region of interest;</li><li>Figure 9 is a diagram showing an example indication of a recent history of the skewness of a frequency distribution of detected values of the tissue characterization parameter in a region of interest;</li><li>Figure 10 is a diagram showing an example indication of a recent history of the kurtosis of a frequency distribution of detected values of the tissue characterization parameter in a region of interest;</li><li>Figure 11 is a diagram showing an example indication of a recent history of the frequency distribution of detected values of a tissue characterization parameter in a region of interest;</li><li>Figure 12 is a diagram showing an example indication in an operation designating a frame by an operator using a cursor;</li><li>Figure 13 is a diagram showing an example indication of a frame which is designated as a frame which gives the minimum standard deviation among a plurality of frames which has been most recently observed; and</li><li>Figure 14 is a diagram showing an example indication of a frequency distribution of values of a tissue characterization parameter in a current frame, where a frequency distribution of values of the tissue characterization parameter in the frame which gives the minimum standard deviation is superimposed on the frequency distribution in the current frame.</li></ul>

DESCRIPTION OF THE PREFERRED EMBODIMENTS

Fig. 1

0030Figure 1 is a block diagram showing a construction of the embodiment of the present invention. In Fig. 1, reference numeral 1 denotes an ultrasonic transducer array, 2₁, 2₂, 2₃, ... 2_n and 6₁, 6₂, 6₃, ... 6_n each denote an amplifier, 3₁, 3₂, 3₃, ... 3_n and 7₁, 7₂, 7₃, ... 7_n each denote a delay circuit, 4 denotes a pulse generator, 5 denotes a timing generator, 8 denotes a summing amplifier, 9 denotes a LOG amplifier, 9' denotes a time gain control circuit, 10 denotes a detector, 11 denotes a low pass filter, 12 denotes an analog to digital converter, 100 denotes a parameter calculation circuit, 200 denotes a characteristic portion of the embodiment of the present invention, and 300 denotes a display unit.

0031The timing generator 5 generates a timing signal at a constant rate, and the pulse generator 4 generates a pulse signal to generate an ultrasound pulse. Each delay circuit 3_i and 7_i (i=1, 2, ... n) is connected in series to a corresponding one of the amplifiers 2_i and 6_i (i=1, 2, ... n). The pulse signal is input in parallel to the respective series connections of the delay circuits 3_i (i=1, 2, ... n) and the amplifiers 2_i (i=1, 2, ... n). The ultrasonic transducer array 1 comprises a plurality of transducer elements. Each transducer element comprises a piezo-electric element which converts an electric signal to a mechanical deformation to generate an ultrasound, and generates an electric signal when the element receives an ultrasound. The transducer array 1 is manually applied to a skin over a region of tissue of the human body in which a doctor is interested, to diagnose the region. The above outputs of the series connections parallel connection are respectively applied to the plurality of transducer elements of the ultrasonic transducer array 1. Delay times in the above delay circuits 3_i (i=1, 2, ... n) are controlled so that the ultrasound which is generated at the above transducer elements are synthesized to form an ultrasound beam directed in a certain direction for each cycle of the output of the ultrasound, and the direction is moved to scan a cross section of the tissue for each scanning cycle.

0032The electric signals which are output from the plurality of transducer elements are respectively input in parallel into the series connections of the delay circuits 7_i (i=1, 2, ... n) and the amplifiers 6_i (i=1, 2, ... n). The outputs of the series connections of the delay circuits 7_i (i=1, 2, ... n) and the amplifiers 6_i (i=1, 2, ... n) are summed in the summing amplifier 8 to apply the summed signal to the LOG amplifier 9. Delay times in the above delay circuits 7_i (i=1, 2, ... n) are controlled so that the summed signal corresponds to ultrasound which is received in the transducer array 1 from the above direction of the above ultrasound beam which is generated in the transducer array 1, and the direction is moved to scan a cross section of the tissue for each scanning cycle, synchronized with the above scanning operation by the delay circuits 3_i (i=1, 2, ... n). The LOG amplifier 9 amplifies the summed signal. The gain of the LOG amplifier 9 is controlled by the time gain control circuit 9' so that the gain of the LOG amplifier 9 is increased with a time which is elapsed from each time of the above generation of the ultrasound in the transducer array 1 in accordance with a logarithm line. The output of the LOG amplifier 9 is applied in parallel to the detector 10 and the parameter calculation circuit 100.

0033The envelope of the above output of the LOG amplifier 9 is detected in the detector 10, and the detected signal is input into the low pass filter 11. The output of the low pass filter 11 corresponds to image data of ultrasonograph. The output of the low pass filter 11 is converted to a digital form by the analog to digital converter 12. The digital data is supplied to the characteristic portion 200.

Parameter Calculation Circuit (Fig. 2)

0034Figure 2 is a block diagram showing a construction of an example of the parameter calculation circuit 100 in Fig. 1. In Fig. 2, reference numeral 54 denotes a gradient calculation circuit, 55 denotes a center frequency detector, 56 denotes a divider, and 57 denotes a multiplier.

0035The gradient calculation circuit 54 obtains a gradient of an envelope of the above output of the LOG amplifier 9 in a digital form, and the center frequency detector 55 detects a center frequency of the output of the LOG amplifier 9 in a digital form. The divider 56 divides the output of the gradient calculation circuit 54 by the output of the center frequency detector 55, and the multiplier 57 multiplies the output of the divider 56 by a constant -1/2c to obtain the ultrasonic attenuation coefficient. The obtained ultrasonic attenuation coefficient is supplied to the above characteristic portion 200. The detail of the construction for obtaining the ultrasonic attenuation coefficient is disclosed in the Japanese Unexamined Patent Publication No. 62-109553. The parameter calculation circuit 100 may be another construction for obtaining an ultrasonic reflection coefficient, or another elasticity parameter. These constructions for obtaining the ultrasonic reflection coefficient, and the other elasticity parameters are conventionally known.

Characteristic Portion (Fig. 3)

0036Figure 3 is a block diagram showing a construction of the characteristic portion of the embodiment of the present invention. In Fig. 3, reference numeral 31 denotes a frame memory for the two-dimensional distribution of a tissue characterization parameter, 32 denotes a statistic calculation circuit, 33 denotes a frequency distribution obtaining circuit, 34 denotes a scroll memory, 35 denotes a memory, 36 denotes a digital scan converter, 37 denotes a region of interest designation circuit, 38 denotes a frame designation circuit, 39 denotes a minimum dispersion frame obtaining circuit, 40 denotes a freeze control circuit, 41 denotes a frame memory for image data of ultrasonograph, and 42 denotes a CRT.

0037The above digital image data of ultrasonograph is stored in the frame memory 41, and the above digital tissue characterization parameter is stored in the frame memory 31. These frame memories respectively hold the data of a plurality of frames which have been most recently observed. The statistic calculation circuit 32 obtains statistic information on a dispersing characteristic of the characteristic parameter of the tissue at the various points in a region of interest of each cross section. The region of interest is designated by the region of interest designation circuit 37 as explained later.

Statistic Information

0038The statistic information is, for example, an average of the values of the tissue characterization parameter which was obtained in the parameter calculation circuit 100 within the region of interest in each cross section, a standard deviation of the values of the tissue characterization parameter which was obtained in the parameter calculation circuit 100 within the region of interest in each cross section, a skewness of a frequency distribution of the values of the tissue characterization parameter which was obtained in the parameter calculation circuit 100 within the region of interest in each cross section, or the kurtosis of a frequency distribution of the values of the tissue characterization parameter which was obtained in the parameter calculation circuit 100 within the region of interest in each cross section. When the values of the tissue characterization parameter in the region of interest are denoted by x_i (i=1,2, ... N) where N denotes the number of the obtained values in the region of interest.

0039The average m is obtained by the following equation (1).<maths id="math0001" num=""><img file="EP0459813A2_D0001.tif" /></maths>

0040The variance σ² is obtained by the following equation (2).<maths id="math0002" num=""><img file="EP0459813A2_D0002.tif" /></maths>

0041The standard deviation σ is obtained by the following equation (3).<maths id="math0003" num="(3)"><math display="block"><mrow><mtext>σ = √σ²</mtext></mrow></math><img file="EP0459813A2_D0003.tif" /></maths>

0042The third order center moment µ₃ is obtained as the following equation (4).<maths id="math0004" num=""><img file="EP0459813A2_D0004.tif" /></maths>

0043The fourth order center moment µ₄ is obtained by the following equation (5).<maths id="math0005" num=""><img file="EP0459813A2_D0005.tif" /></maths>

0044The skewness S of the frequency distribution of the values of the tissue characterization parameter in the region of interest indicates asymmetry of the frequency distribution, and is obtained by the following equation (6).<maths id="math0006" num="(6)"><math display="block"><mrow><mtext>S=µ₃/σ³</mtext></mrow></math><img file="EP0459813A2_D0006.tif" /></maths>

0045The kurtosis K of the frequency distribution of the values of the tissue characterization parameter in the region of interest indicates flatness of the frequency distribution, and is obtained by the following equation (6).<maths id="math0007" num="(7)"><math display="block"><mrow><mtext>K=µ₄/σ⁴ -3</mtext></mrow></math><img file="EP0459813A2_D0007.tif" /></maths>

0046The results of the statistic calculations as above are input into the scroll memory 34 for each frame (cross section). The scroll memory 34 stores the statistic information for a predetermined number of frames which are most recently observed.

Frequency Distribution Obtaining Circuit (Fig. 4)

0047The frequency distribution obtaining circuit 33 obtains a frequency distribution of the values of the tissue characterization parameter in the region of interest in each cross section. Figure 4 is a block diagram showing a construction of the frequency distribution obtaining circuit 33 in Fig. 3. In Fig. 4, reference numeral 50 denotes a decoder, and 51₁, 51₂, 51₃, ... 51_n each denote an accumulator.

0048The values of the tissue characterization parameter are sorted into a plurality of ranks (ranges of the values) to obtain a histogram of the frequency distribution. The accumulators 51₁, 51₂, 51₃, ... 51_n are provided for the respective ranks. The decoder 50 receives respective values of the tissue characterization parameter which are stored in the frame memory 31. The decoder 50 has a plurality of outputs corresponding to the plurality of ranks to realize the function of the sorting. One of the outputs of the decoder 50 corresponding to one of the ranks to which the received value is to be sorted, becomes active when the enable signal is active. The enable signal is generated by the region of interest designation circuit 37, and is made active when values of the tissue characterization parameter in the region of interest are output from the frame memory 31. A reset signal is applied to the plurality of accumulators 51₁, 51₂, 51₃, ... 51_n when the accumulators should be reset. The accumulation is carried out for example, for each scanning cycle (each cross section), or for a predetermined number of successive scanning cycles. Figure 5 is an example histogram showing a frequency distribution of values of the tissue characterization parameter in the region of interest. The above obtained frequency distribution is input into the memory 35 for each frame (cross section). The memory 35 stores the data of the frequency distributions for a predetermined number of frames which are most recently observed.

Fig. 6

0049Figure 6 is a diagram showing an example of the display of the two-dimensional distribution of detected values of a tissue characterization parameter and an indication of the designated region of interest. Viewing the display, the operator can designate a desired region of interest, for example, by a pointing device. The region of interest designation circuit 37 controls the operation reading the values of the tissue characterization parameter from the frame memory 31 and the operation writing the values into the statistic calculation circuit 32 and the frequency distribution obtaining circuit 33 by supplying corresponding addresses and reading and writing signals. The technique of designating a region of interest is well-known.

Indication of History of Statistics (Figs.3, 7, 8, 9, and 10)

0050In an indication mode, a recent history of one or more of the aforementioned statistic information in a predetermined number of successive frames ending with a current frame is indicated in a real time mode as shown in Figs. 7 to 11. The above history data is supplied from the scroll memory 34 and the armory 35. In the above indication, each statistic information in the history may be calculated for a predetermined number of successive frames instead of one frame. Figure 7 is a diagram showing an example indication of a recent history of an average of a frequency distribution of detected values of a tissue characterization parameter in a region of interest, Figure 8 is a diagram showing an example indication of a recent history of a standard deviation of a frequency distribution of detected values of a tissue characterization parameter in a region of interest, Figure 9 is a diagram showing an example indication of a recent history of the skewness of a frequency distribution of detected values of the tissue characterization parameter in a region of interest, Figure 10 is a diagram showing an example indication of a recent history of the kurtosis of a frequency distribution of detected values of the tissue characterization parameter in a region of interest, and Figure 11 is a diagram showing an example indication of a recent history of the frequency distribution of detected values of a tissue characterization parameter in a region of interest. In Fig. 11, the lines for the frequency distributions of the different frames may be indicated by the different types of lines, respectively, for example, by using different colors. The operator can view the above indication while manipulating the probe. Therefore, for example, the operator can recognize whether or not a most recently variation of a position of the probe gives a better (more stable and credible) result, and can avoid to take unstable data based on the above history indications.

Freeze Control and Frame Designation (Figs. 3 and 12)

0051When the operator finds in the above history indication that credible data has been obtained, the operator can freeze an operation renewing the contents of the scroll memory 34 and the memory 35, and thus, can freeze the history indication. The freezing operation can be carried out through the freeze control circuit 40. When the freeze control circuit 40 receives a command to effect the freeze operation, the freeze control circuit 40 supplies a freeze control signal to the frame memories 31 and 41, the scroll memory 34, and the memory 35 to freeze these elements, as shown in Fig. 3.

0052Figure 12 is a diagram showing an example indication in an operation designating a frame by an operator using a cursor. As shown in Fig. 12, the operator can designate one of the frames for which the history of the statistic information is indicated, based on the indication, by moving a cursor to the frame in the history indication using a pointing device (not shown). This operation is realized through the frame designation circuit 38. Responding to the designation by the operator, the frame designation circuit 38 can control the frame memories 31 and 41, the scroll memory 34, and the memory 35, so that these memories output contents thereof regarding the designated frame, to the digital scan converter 36. Further, the frame designation circuit 38 can send the information on the designated frame to the digital scan converter 36 so that the above cursor indication is realized.

0053Further, the aforementioned region of interest designation circuit 37 can be operated again to reobtain statistic information on the dispersing characteristic of the predetermined characteristic parameter in a region of interest which is different from the region of interest which was first designated. Redesignating the new region of interest through the region of interest designation circuit 37 after freezing the above memories in Fig. 3, the above statistical information can be obtained by using the frozen data of the tissue characterization parameter, and operating the statistic calculation circuit 32 and/or the frequency distribution obtaining circuit 33.

Calculation of Most Stable Frame (Figs. 3, 13, and 14)

0054The most stable frame among a plurality of frames which have been most recently observed, can be automatically obtained by using the recent history data of the statistic information.

0055The minimum dispersion frame obtaining circuit 39 obtains a frame of a minimum standard deviation. Although not shown, the minimum dispersion frame obtaining circuit 39 may comprise registers for holding a frame number of a temporarily minimum standard deviation, and the temporarily minimum standard deviation, and a comparator for comparing a standard deviation of a current frame with the held standard deviation, and the minimum dispersion frame obtaining circuit 39 replace the held frame number and the standard deviation with the current frame number and the standard deviation of the current frame when the current standard deviation is not larger than the held standard deviation. When a command to indicate the frame of the minimum standard deviation which is determined by the minimum dispersion frame obtaining circuit 39 is input, the minimum dispersion frame obtaining circuit 39 sends a marker indication data to the digital scan converter 36. The marker indication data indicates how many frames before the current frame the frame of the minimum standard deviation is. When the aforementioned control circuit of the digital scan converter 36 receives the marker indication data, a marker indication is shown superimposed on the above history indication of the statistic information. Figure 13 is a diagram showing an example indication of a frame which is designated as a frame which gives the minimum standard deviation among a plurality of frames which has been most recently observed.

0056Further, a histogram indicating a frequency distribution of the values of the tissue characterization parameter in a current frame can be indicated and superimposed on another histogram indicating a frequency distribution of the values of the tissue characterization parameter in the above frame of the minimum standard deviation, as shown in Fig. 14. This indication is realized by controlling the memory 35 through the frame designation circuit 38 by another output of the minimum dispersion frame obtaining circuit 39. The output of the minimum dispersion frame obtaining circuit 39 includes information of the frame of the minimum standard deviation. When the frame designation circuit 38 receives the above output, the frame designation circuit 38 receives the above output, the frame designation circuit 38 can control the frame memories 31 and 41, the scroll memory 34, and the memory 35, so that these memories output contents thereof regarding the frame of the minimum standard deviation, to the digital scan converter 36.

0057Otherwise, using the function of the minimum dispersion frame obtaining circuit 39, it is possible to set the ultrasonic imaging apparatus of Fig. 1 so that the apparatus automatically indicates the image of the ultrasonograph and the two-dimensional distribution image of the tissue characterization parameter of the frame of the minimum standard deviation by controlling the frame memories 31 and 41 through the frame designation circuit 38.

0058Although not shown, all the above operations in the embodiment can be realized by either a hardware logic circuit or software on a microcomputer.