WO2012011414A1 - Dispositif de diagnostic par ultrasons, procédé et logiciel de mise en œuvre correspondants - Google Patents

Dispositif de diagnostic par ultrasons, procédé et logiciel de mise en œuvre correspondants Download PDF

Info

Publication number
WO2012011414A1
WO2012011414A1 PCT/JP2011/065910 JP2011065910W WO2012011414A1 WO 2012011414 A1 WO2012011414 A1 WO 2012011414A1 JP 2011065910 W JP2011065910 W JP 2011065910W WO 2012011414 A1 WO2012011414 A1 WO 2012011414A1
Authority
WO
WIPO (PCT)
Prior art keywords
specimen
tissue property
feature
tissue
feature amount
Prior art date
Application number
PCT/JP2011/065910
Other languages
English (en)
Japanese (ja)
Inventor
浩仲 宮木
弘孝 江田
裕雅 野口
忠明 神原
安広 和田
Original Assignee
オリンパスメディカルシステムズ株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by オリンパスメディカルシステムズ株式会社 filed Critical オリンパスメディカルシステムズ株式会社
Priority to JP2012524951A priority Critical patent/JPWO2012011414A1/ja
Publication of WO2012011414A1 publication Critical patent/WO2012011414A1/fr

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023—Details of receivers
    • G01S7/52025—Details of receivers for pulse systems
    • G01S7/52026—Extracting wanted echo signals
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023—Details of receivers
    • G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52042—Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/5205—Means for monitoring or calibrating

Definitions

  • the present invention relates to an ultrasonic diagnostic apparatus for diagnosing a tissue property of a specimen using ultrasonic waves, an operating method of the ultrasonic diagnostic apparatus, and an operating program for the ultrasonic diagnostic apparatus.
  • Ultrasonic elastography is a technique that utilizes the fact that the hardness of cancer and tumor tissue in a living body varies depending on the progress of the disease and the living body.
  • the amount of strain and elastic modulus of the biological tissue at the examination location are measured using ultrasound while the examination location is pressed from the outside, and the measurement result is displayed as a tomographic image.
  • the ultrasonic elastography described above has a problem in that the pressure applied to the lower part of a blood vessel or a lymph vessel is difficult to be transmitted. Therefore, when a tumor is formed in the vicinity of the blood vessel, the boundary of the tumor is unclear and it is difficult to distinguish the invasion of the tumor into the blood vessel. As described above, in ultrasonic elastography, there are cases in which it is not possible to accurately distinguish tissue properties.
  • the present invention has been made in view of the above, and enables an ultrasonic diagnostic apparatus and an ultrasonic diagnostic apparatus capable of accurately distinguishing tissue properties and improving the reliability of measurement results. It is an object to provide an operating method and an operating program for an ultrasonic diagnostic apparatus.
  • an ultrasonic diagnostic apparatus transmits ultrasonic waves to a specimen to be diagnosed and receives ultrasonic waves reflected by the specimen.
  • An ultrasonic diagnostic apparatus for diagnosing the tissue characteristics of the specimen based on the received ultrasound, the frequency analyzer calculating the frequency spectrum by analyzing the frequency of the received ultrasound, and the frequency analyzer calculating
  • a feature quantity extraction unit for extracting the feature quantity of the frequency spectrum by approximating the frequency spectrum, and the feature quantity of the frequency spectrum extracted based on the ultrasonic waves respectively reflected by the plurality of known specimens.
  • a storage unit that stores information in association with tissue properties of known samples, and a feature that the storage unit stores in association with tissue properties of the plurality of known samples. And characterized in that and a tissue characterization determining unit determines the tissue properties of a predetermined area of the specimen by using a feature quantity the feature amount extraction unit and extracted.
  • the ultrasonic diagnostic apparatus is characterized in that, in the above invention, the feature amount extraction unit approximates the frequency spectrum by a polynomial by regression analysis.
  • the feature amount extraction unit approximates the frequency spectrum with a linear expression, the slope of the linear expression, the intercept of the linear expression, and the frequency spectrum.
  • a plurality of feature quantities including a frequency included in the frequency band, an intensity determined by the slope and the intercept are extracted.
  • the storage unit stores an average of each feature amount in a group classified for each tissue property with respect to the plurality of known samples, and the tissue property
  • the determination unit sets a feature amount space including at least one of the plurality of feature amounts as a component, and sets the feature amount constituting the component of the feature amount space among the feature amounts of the frequency spectrum in the predetermined region of the specimen.
  • a specimen average point having an average as coordinates of the feature quantity space, and an average of feature quantities constituting components of the feature quantity space among the feature quantities in the group of the plurality of known specimens as coordinates of the feature quantity space
  • the tissue property of the specimen is determined based on the distance in the feature amount space with the known specimen average point.
  • the tissue property determination unit corresponds to a group in which a distance in the feature amount space between the specimen average point and the known specimen average point is minimum.
  • the tissue property is the tissue property of the specimen.
  • the tissue property determination unit has a larger value as a distance in the feature amount space between the specimen average point and the known specimen average point is smaller.
  • the probability that the sum of the values corresponding to all the distances is 1 is calculated with respect to all the distances, whereby the tissue properties of the specimen are determined probabilistically.
  • the ultrasonic diagnostic apparatus in the above invention, generates visual information corresponding to the feature amount of the specimen, the generated visual information, an image generated based on the received ultrasonic wave, and A determination result display image data generation unit for generating determination result display image data for displaying the determination result of the tissue property of the predetermined region of the specimen by using the result determined by the tissue property determination unit; To do.
  • the visual information is a variable constituting a color space.
  • the operation method of the ultrasonic diagnostic apparatus includes transmitting the ultrasonic wave to the specimen to be diagnosed and receiving the ultrasonic wave reflected by the specimen, whereby the specimen based on the received ultrasonic wave
  • the method of operating an ultrasonic diagnostic apparatus for diagnosing the tissue properties of a frequency analysis step of calculating a frequency spectrum by a frequency analysis unit by analyzing the frequency of the received ultrasonic wave, and the frequency calculated in the frequency analysis step A feature amount extraction step for extracting a feature amount of the frequency spectrum by approximating the spectrum by a feature amount extraction unit, and a feature amount of the frequency spectrum extracted based on ultrasonic waves respectively reflected by a plurality of known specimens
  • the operation program of the ultrasonic diagnostic apparatus includes transmitting the ultrasonic wave to the sample to be diagnosed and receiving the ultrasonic wave reflected by the sample, whereby the sample based on the received ultrasonic wave
  • a feature amount extraction step for extracting a feature amount of the frequency spectrum by a feature amount extraction unit; a feature amount of the frequency spectrum extracted based on ultrasonic waves respectively reflected by a plurality of known specimens;
  • the feature amount read from the storage unit that is stored in association with the property and the feature amount extraction step Characterized in that to perform tissue characterization determining step determines the tissue characterization determining unit tissue characterization of a predetermined area of said specimen by in using the extracted feature quantity.
  • FIG. 1 is a block diagram showing the configuration of the ultrasonic diagnostic apparatus according to Embodiment 1 of the present invention.
  • FIG. 2 is a flowchart showing an outline of processing of the ultrasonic diagnostic apparatus according to Embodiment 1 of the present invention.
  • FIG. 3 is a diagram showing a display example of a B-mode image on the display unit of the ultrasonic diagnostic apparatus according to Embodiment 1 of the present invention.
  • FIG. 4 is a flowchart showing an outline of processing performed by the frequency analysis unit of the ultrasonic diagnostic apparatus according to Embodiment 1 of the present invention.
  • FIG. 5 is a diagram schematically showing a data array of one sound ray.
  • FIG. 6 is a diagram illustrating an example (first example) of a frequency spectrum calculated by the frequency analysis unit of the ultrasonic diagnostic apparatus according to Embodiment 1 of the present invention.
  • FIG. 7 is a diagram illustrating an example (second example) of a frequency spectrum calculated by the frequency analysis unit of the ultrasonic diagnostic apparatus according to Embodiment 1 of the present invention.
  • FIG. 8 is a flowchart showing an outline of processing performed by the tissue property determination unit of the ultrasonic diagnostic apparatus according to Embodiment 1 of the present invention.
  • FIG. 9 is a diagram illustrating an example of a feature amount space set by the tissue property determination unit of the ultrasonic diagnostic apparatus according to Embodiment 1 of the present invention.
  • FIG. 10 is a diagram illustrating a display example of a determination result display image displayed by the display unit of the ultrasonic diagnostic apparatus according to Embodiment 1 of the present invention.
  • FIG. 11 is a diagram illustrating another display example of the determination result display image displayed by the display unit of the ultrasonic diagnostic apparatus according to Embodiment 1 of the present invention.
  • FIG. 12 is a diagram illustrating an overview of the tissue property determination process performed by the tissue property determination unit of the ultrasonic diagnostic apparatus according to Embodiment 4 of the present invention.
  • FIG. 1 is a block diagram showing the configuration of the ultrasonic diagnostic apparatus according to Embodiment 1 of the present invention.
  • An ultrasonic diagnostic apparatus 1 shown in FIG. 1 is an apparatus for diagnosing the tissue properties of a specimen to be diagnosed using ultrasonic waves.
  • the ultrasonic diagnostic apparatus 1 transmits and receives electrical signals between the ultrasonic probe 2 that outputs an ultrasonic pulse to the outside and receives an ultrasonic echo reflected from the outside, and the ultrasonic probe 2. Transmitting / receiving unit 3 for performing the calculation, calculation unit 4 for performing a predetermined calculation on the electrical echo signal converted from the ultrasonic echo, and generation of image data corresponding to the electrical echo signal converted from the ultrasonic echo An image processing unit 5 to be performed, an interface such as a keyboard, a mouse, a touch panel, and the like.
  • the image processing unit 5 is realized by using an input unit 6 that receives input of various information and a display panel made of liquid crystal or organic EL.
  • a display unit 7 that displays various types of information including the image generated
  • a storage unit 8 that stores various types of information including information related to the tissue properties of the known specimen
  • a control unit that controls the operation of the ultrasonic diagnostic apparatus 1. And, equipped with a.
  • the ultrasonic probe 2 converts an electrical pulse signal received from the transmission / reception unit 3 into an ultrasonic pulse (acoustic pulse signal), and converts an ultrasonic echo reflected by an external specimen into an electrical echo signal.
  • a signal conversion unit 21 for conversion is included.
  • the ultrasonic probe 2 may be one that mechanically scans an ultrasonic transducer, or one that electronically scans a plurality of ultrasonic transducers.
  • the transmitting / receiving unit 3 is electrically connected to the ultrasonic probe 2, transmits a pulse signal to the ultrasonic probe 2, and receives an echo signal from the ultrasonic probe 2. Specifically, the transmission / reception unit 3 generates a pulse signal based on a preset waveform and transmission timing, and transmits the generated pulse signal to the ultrasound probe 2. Further, the transmission / reception unit 3 performs processing such as amplification and filtering on the received echo signal, and then performs A / D conversion to generate and output a digital RF signal.
  • the transmission / reception unit 3 has a multi-channel circuit for beam synthesis corresponding to the plurality of ultrasonic transducers.
  • the calculation unit 4 is obtained by a frequency analysis unit 41 that performs a frequency analysis of an echo signal by performing a fast Fourier transform (FFT) on the digital RF signal output from the transmission / reception unit 3, and a frequency analysis performed by the frequency analysis unit 41.
  • the feature quantity extraction unit 42 that extracts the feature quantity of the frequency spectrum by approximating the obtained frequency spectrum (power spectrum), and the tissue quantity of the predetermined region of the specimen are determined using the feature quantity extracted by the feature quantity extraction unit 42 And a tissue property determination unit 43.
  • the frequency analysis unit 41 calculates a frequency spectrum by performing a fast Fourier transform on an FFT data group having a predetermined amount of data for each sound ray (line data).
  • the frequency spectrum shows different tendencies depending on the tissue properties of the specimen. This is because the frequency spectrum has a correlation with the size, density, acoustic impedance, and the like of the specimen as a scatterer that scatters ultrasonic waves.
  • “intensity” refers to any of parameters such as voltage, power, sound pressure, and acoustic energy.
  • the inclination a has a correlation with the size of the ultrasonic scatterer, and it is generally considered that the larger the scatterer, the lower the inclination.
  • the intercept b has a correlation with the size of the scatterer, the difference in acoustic impedance, the density (concentration) of the scatterer, and the like.
  • the intensity c is an indirect parameter derived from the slope a and the intercept b, and gives the spectral intensity at the center in the effective frequency band. Therefore, the intensity c is considered to have a certain degree of correlation with the brightness of the B-mode image in addition to the size of the scatterer, the difference in acoustic impedance, and the density of the scatterer.
  • the approximate polynomial calculated by the feature amount extraction unit 42 is not limited to a linear expression, and it is possible to use a quadratic or higher approximate polynomial.
  • the tissue property determination unit 43 calculates the average and standard deviation of the feature amounts of the frequency spectrum extracted by the feature amount extraction unit 42 for each feature amount.
  • the tissue property determination unit 43 determines the tissue property of a predetermined region of the sample by using the calculated average and standard deviation and the average and standard deviation of the feature amounts of the frequency spectrum of the known sample stored in the storage unit 8. .
  • the “predetermined region” referred to here is a region in the image (hereinafter referred to as “region of interest”) designated by the input unit 6 by the operator of the ultrasonic diagnostic apparatus 1 who viewed the image generated by the image processing unit 5. That is.
  • tissue property is, for example, any of cancer, endocrine tumor, mucinous tumor, normal tissue, vascular and the like.
  • tissue properties include chronic pancreatitis, autoimmune pancreatitis and the like.
  • the average and standard deviation of the feature values calculated by the tissue characterization determining unit 43 is a systematic change such as a change in cell level such as nuclear enlargement or anomaly, an increase in fibers in the stroma, or a replacement of the real tissue with fibers. This is a unique value depending on the tissue properties. Therefore, by using such an average and standard deviation of the feature amounts, it is possible to accurately determine the tissue properties of a predetermined region of the specimen.
  • the image processing unit 5 outputs the B-mode image data generating unit 51 that generates the B-mode image data to be displayed by converting the amplitude of the echo signal into the luminance, the B-mode image data generating unit 51, and the calculation unit 4, respectively.
  • a determination result display image data generation unit 52 that generates determination result display image data for displaying the determination result of the tissue property of the region of interest and information related to the determination result using the data.
  • the B-mode image data generation unit 51 performs signal processing using a known technique such as a bandpass filter, logarithmic conversion, gain processing, contrast processing, and the like on the digital signal, and also according to the image display range on the display unit 7.
  • B-mode image data is generated by thinning out data in accordance with the data step width determined in advance.
  • the determination result display image data generation unit 52 uses the B mode image data generated by the B mode image data generation unit 51, the feature amount calculated by the feature amount extraction unit 42, and the determination result determined by the tissue property determination unit 43.
  • determination result display image data including the determination result of the tissue property of the region of interest and the tissue property emphasized image that emphasizes the tissue property is generated.
  • the storage unit 8 includes a known sample information storage unit 81 that stores information on known samples, and a window function storage unit 82 that stores a window function used in the frequency analysis process performed by the frequency analysis unit 41.
  • the known specimen information storage unit 81 stores the characteristic amount of the frequency spectrum extracted by the frequency analysis for the known specimen in association with the tissue property of the known specimen.
  • the known specimen information storage unit 81 uses the average and standard deviation calculated for each group classified for each tissue property of the known specimen with respect to the feature quantity of the frequency spectrum related to the known specimen. It is memorized with all the data.
  • the information on the known specimen stored in the known specimen information storage unit 81 is desirably information with high reliability related to the tissue properties.
  • the window function storage unit 82 stores at least one of window functions such as Hamming, anning, and Blackman.
  • the storage unit 8 includes a ROM in which an operation program for the ultrasonic diagnostic apparatus according to the first embodiment, a program for starting a predetermined OS, and the like are stored in advance, and a RAM in which calculation parameters and data for each process are stored. It is realized using.
  • Components other than the ultrasound probe 2 of the ultrasound diagnostic apparatus 1 having the above functional configuration are realized using a computer having a CPU having a calculation and control function.
  • the CPU provided in the ultrasonic diagnostic apparatus 1 reads out various programs including information stored and stored in the storage unit 8 and the above-described operation program of the ultrasonic diagnostic apparatus from the storage unit 8, so that the ultrasonic according to the first embodiment is obtained. Arithmetic processing related to the operation method of the ultrasonic diagnostic apparatus is executed.
  • the operation program of the ultrasonic diagnostic apparatus may be recorded on a computer-readable recording medium such as a hard disk, a flash memory, a CD-ROM, a DVD-ROM, or a flexible disk and widely distributed. Is possible.
  • FIG. 2 is a flowchart showing an outline of processing of the ultrasonic diagnostic apparatus 1 having the above configuration.
  • the ultrasound diagnostic apparatus 1 first measures a new specimen by the ultrasound probe 2 (step S1). Thereafter, the B mode image data generation unit 51 generates B mode image data (step S2).
  • FIG. 3 is a diagram illustrating a display example of the B-mode image on the display unit 7.
  • a B-mode image 100 shown in the figure is a grayscale image in which values of red (R), green (G), and blue (B), which are variables when the RGB color system is adopted as a color space, are matched. .
  • step S4: Yes the frequency analysis unit 41 calculates a frequency spectrum by performing frequency analysis by FFT calculation (step S5). In this step S5, it is also possible to set the entire region of the image as the region of interest.
  • step S5 it is also possible to set the entire region of the image as the region of interest.
  • step S6: Yes when an instruction to end the process is input by the input unit 6 (step S6: Yes), the ultrasound diagnostic apparatus 1 ends the process. .
  • step S6: No when the region of interest is not designated (step S4: No), when the instruction to end the process is not input by the input unit 6 in step S6 (step S6: No), the ultrasonic diagnostic apparatus 1 Return to step S4.
  • the frequency analysis unit 41 sets the sound ray number L of the sound ray to be analyzed first as an initial value L 0 (step S11).
  • the initial value L 0 may be given, for example, to a sound ray that is first received by the transmission / reception unit 3, or for a sound ray corresponding to one of the left and right boundary positions of the region of interest set by the input unit 6. May be given.
  • FIG. 5 is a diagram schematically showing a data array of one sound ray.
  • a white or black rectangle means one piece of data.
  • the sound ray LD is discretized at a time interval corresponding to a sampling frequency (for example, 50 MHz) in A / D conversion performed by the transmission / reception unit 3.
  • FIG. 5 shows a case where the first data of the sound ray LD is set as the initial value Z 0 of the data position Z.
  • FIG. 5 is merely an example, and the position of the initial value Z 0 can be arbitrarily set.
  • the data position Z corresponding to the upper end position of the region of interest may be set as the initial value Z 0 .
  • the frequency analysis unit 41 acquires the FFT data group at the data position Z (step S13), and causes the window function stored in the window function storage unit 82 to act on the acquired FFT data group (step S14). In this way, by applying the window function to the FFT data group, it is possible to avoid the FFT data group from becoming discontinuous at the boundary and to prevent the occurrence of artifacts.
  • the frequency analysis unit 41 determines whether or not the FFT data group at the data position Z is a normal data group (step S15).
  • the FFT data group needs to have a power number of 2 data.
  • the number of data in the FFT data group is 2 n (n is a positive integer). That the FFT data group is normal means that the data position Z is the 2 n-1 th position from the front in the FFT data group.
  • the FFT data groups F 2 , F 3 and F K-1 are normal, while the FFT data groups F 1 and F K are abnormal.
  • step S15 If the result of determination in step S15 is that the FFT data group at the data position Z is normal (step S15: Yes), the frequency analysis unit 41 proceeds to step S17 described later.
  • step S15 If the result of determination in step S15 is that the FFT data group at the data position Z is not normal (step S15: No), the frequency analysis unit 41 generates a normal FFT data group by inserting zero data for the shortage (Ste S16).
  • the FFT function group determined to be not normal in step S15 is subjected to a window function before adding zero data. For this reason, discontinuity of data does not occur even if zero data is inserted into the FFT data group.
  • step S16 the frequency analysis unit 41 proceeds to step S17 described later.
  • step S17 the frequency analysis unit 41 obtains a frequency spectrum by performing an FFT operation using the FFT data group (step S17).
  • 6 and 7 are diagrams illustrating examples of frequency spectra calculated by the frequency analysis unit 41.
  • FIG. 6 and 7, the horizontal axis f is the frequency
  • the vertical axis I is the intensity.
  • the lower limit frequency f LOW and the upper limit frequency f HIGH of the frequency spectrum are the frequency band of the ultrasonic probe 2 and the pulse signal transmitted by the transmitting / receiving unit 3.
  • f LOW 3 MHz
  • f HIGH 10 MHz.
  • the straight line L 1 shown in FIG. 6 and the straight line L 2 shown in FIG. 7 will be described in a feature amount extraction process described later.
  • the curve and the straight line are composed of a set of discrete points. This also applies to the embodiments described later.
  • the frequency analysis unit 41 adds a predetermined data step width D to the data position Z to calculate the data position Z of the next FFT data group to be analyzed (step S18).
  • the data step width D here is preferably the same as the data step width used when the B-mode image data generation unit 51 generates the B-mode image data. However, when it is desired to reduce the amount of calculation in the frequency analysis unit 41 In this case, a value larger than the data step width used by the B-mode image data generation unit 51 may be set.
  • the frequency analysis unit 41 determines whether or not the data position Z is greater than the final data position Z max (step S19).
  • the final data position Z max may be the data length of the sound ray LD, or may be the data position corresponding to the lower end of the region of interest.
  • the frequency analysis unit 41 increments the sound ray number L by 1 (step S20).
  • the frequency analysis unit 41 returns to step S13.
  • [X] represents the maximum integer not exceeding X.
  • step S21: Yes When the sound ray number L after being incremented in step S20 is larger than the final sound ray number Lmax (step S21: Yes), the frequency analysis unit 41 returns to the main routine shown in FIG. On the other hand, when the sound ray number L after being incremented in step S20 is equal to or less than the final sound ray number Lmax (step S21: No), the frequency analysis unit 41 returns to step S12.
  • the frequency analysis unit 41 performs K FFT operations for each of (L max ⁇ L 0 +1) sound rays.
  • the final sound ray number L max may be given to the last sound ray received by the transmission / reception unit 3, for example, or may be given to the sound ray corresponding to either the left or right boundary of the region of interest.
  • P be the total number of FFT operations (L max ⁇ L 0 +1) ⁇ K performed by the frequency analysis unit 41 for all sound rays.
  • the feature amount extraction unit 42 extracts the feature amount by performing regression analysis on the P frequency spectra calculated by the frequency analysis unit 41 (step S7). Specifically, the feature amount extraction unit 42 calculates three feature amounts a, b, and c by calculating a linear expression that approximates the frequency spectrum of the frequency band f LOW ⁇ f ⁇ f HIGH by regression analysis. .
  • a straight line L 1 shown in FIG. 6 and a straight line L 2 shown in FIG. 7 are regression lines obtained in step S7.
  • the tissue property determination unit 43 determines the tissue property in the region of interest of the sample based on the feature amount extracted by the feature amount extraction unit 42 and the known sample information stored in the known sample information storage unit 81 (step S8).
  • the tissue characterization determining unit 43 calculates the average and standard deviation of each of the three feature values a, b, and c of the Q ( ⁇ P) sets of FFT data located inside the region of interest (step S31). .
  • the tissue property determining unit 43 sets a feature amount space used when determining the tissue property (step S32).
  • the tissue property determining unit 43 sets a feature amount space used when determining the tissue property (step S32).
  • the tissue property determining unit 43 sets a feature amount space used when determining the tissue property (step S32).
  • the feature amount space to be set is determined in advance, but the operator may select a desired feature amount space by the input unit 6.
  • FIG. 9 is a diagram illustrating an example of the feature amount space set by the tissue property determination unit 43.
  • the horizontal axis is the feature quantity b
  • the vertical axis is the feature quantity c.
  • the point Sp shown in FIG. 9 has, as coordinates of the feature amount space, the average of the feature amounts b and c of the frequency spectrum of the FFT data group included in the region of interest of the specimen calculated by the feature amount extracting unit 42 in step S31.
  • a point hereinafter, this point is referred to as “specimen average point”.
  • tissue properties of the known samples stored in the known sample information storage unit 81 are A, B, and C, respectively.
  • the three groups SA, SB, and SC exist in regions that do not intersect with other groups in the feature amount space.
  • the tissue characterization determining unit 43 uses the average of the specimen average point Sp and the feature quantities b and c of the frequency spectrum of the FFT data group included in each of the groups SA, SB, and SC as coordinates in the feature quantity space.
  • the distances ⁇ , ⁇ , ⁇ on the feature amount space between the points A 0 , B 0 , C 0 (hereinafter referred to as “known specimen average points”) are calculated (step S33).
  • known specimen average points are calculated (step S33).
  • the scales of the b-axis component and the c-axis component in the feature amount space are greatly different, it is desirable to appropriately perform weighting for making the contribution of each distance substantially equal.
  • the tissue property determination unit 43 determines the tissue property of the specimen average point Sp based on the distance calculated in step S33 (step S34).
  • the distance ⁇ is the minimum. Therefore, the tissue property determination unit 43 determines that the tissue property of the specimen is A.
  • the tissue property determination unit 43 may output an error signal.
  • the tissue property determination unit 43 may select all tissue properties corresponding to the minimum value as candidates, or any one according to a predetermined rule. One tissue property may be selected. In the latter case, for example, a method of setting a high priority for highly malignant tissue properties such as cancer can be mentioned.
  • the tissue property determination unit 43 may output an error signal.
  • the tissue property determination unit 43 outputs the distance calculation result in step S33 and the determination result in step S34 (step S35). Thereby, the tissue property determination process in step S8 ends.
  • the determination result display image data generation unit 52 includes the B mode image data generated by the B mode image data generation unit 51, the feature amount calculated by the feature amount extraction unit 42, and the tissue property determination unit 43.
  • the determination result display image data is generated by using the determination result determined by (step S9).
  • FIG. 10 is a diagram illustrating a display example of the determination result display image displayed on the display unit 7.
  • the determination result display image 200 shown in the figure is an image display that displays an information display unit 201 that displays various related information including a determination result of tissue properties, and a tissue property emphasized image that emphasizes tissue properties based on a B-mode image. Part 202.
  • the tissue property emphasizing image 300 displayed on the image display unit 202 is equivalent to the B mode image 100 shown in FIG. 3 in which the slices b are equal to R (red), G (green), and B (blue). It is the assigned grayscale image.
  • the display unit 7 displays the determination result display image 200 having the above configuration, the operator can more accurately grasp the tissue properties of the region of interest.
  • the tissue property enhancement image 300 shown in FIG. 10 is merely an example.
  • the tissue property enhanced image can be displayed as a color image.
  • the color space may be configured with complementary color system variables such as cyan, magenta, and yellow, and a feature amount may be assigned to each variable.
  • the tissue property enhanced image data may be generated by mixing B-mode image data and color image data at a predetermined ratio.
  • tissue property enhanced image data may be generated by replacing only the region of interest with color image data.
  • the frequency spectrum feature quantity in the predetermined region of the specimen is extracted by approximating the frequency spectrum obtained by analyzing the frequency of the received ultrasonic wave.
  • the strain amount and elastic modulus of the biological tissue It is possible to clearly distinguish the difference between tissues without using. Therefore, it is possible to accurately distinguish the tissue properties and improve the reliability of the measurement results.
  • FIG. 11 is a diagram illustrating another display example of the determination result display image on the display unit 7.
  • the determination result display image 400 shown in the figure includes an information display unit 401, a first image display unit 402 that displays a B-mode image, and a second image display unit 403 that displays a tissue property emphasized image.
  • the B-mode image 100 is displayed on the first image display unit 402
  • the tissue property emphasized image 300 is displayed on the second image display unit 403.
  • the difference between the two images can be recognized on one screen.
  • the image displayed on the first image display unit 402 and the image displayed on the second image display unit 403 may be interchanged.
  • the display between the determination result display image 200 and the determination result display image 400 may be switched by an input from the input unit 6.
  • the second embodiment of the present invention is different from the first embodiment in the tissue property determination process in the tissue property determination unit.
  • the configuration of the ultrasonic diagnostic apparatus according to the second embodiment is the same as the configuration of the ultrasonic diagnostic apparatus 1 described in the first embodiment. Therefore, in the following description, the same reference numerals are given to the components corresponding to the components of the ultrasonic diagnostic apparatus 1.
  • the tissue property determination unit 43 uses the feature quantities (a, b, c) of the Q sets of FFT data groups located inside the region of interest as the groups SA, SB, and SC constituting the tissue properties A, B, and C (see FIG. 9), a new population is formed, and then a standard deviation for each feature amount of data constituting each tissue property is obtained.
  • the tissue property determination unit 43 performs standard deviation of each feature amount of the groups SA, SB, SC in the original population consisting only of known samples, and groups SA, SB in the new population to which new samples are added, respectively.
  • SC and the standard deviation of each feature quantity (hereinafter, simply referred to as “standard deviation difference”), and the tissue characteristics corresponding to the group including the feature quantity having the smallest standard deviation difference are calculated as the tissue of the specimen. Judged as a property.
  • the tissue property determination unit 43 may calculate the difference of the standard deviation only with respect to the standard deviation of the feature amount selected in advance from a plurality of feature amounts.
  • the feature amount may be selected arbitrarily by the operator, or may be automatically performed by the ultrasonic diagnostic apparatus 1.
  • the tissue property determination unit 43 calculates a value obtained by appropriately weighting and adding the standard deviation differences of all feature amounts for each group, and the tissue property corresponding to the group having the minimum value is determined as the tissue property of the specimen. May be determined.
  • tissue characterization determining unit 43 three feature amounts a, b, respectively corresponding weights c w a, w b, as w c w a Calculate (difference of standard deviation of a) + w b (difference of standard deviation of b) + w c (difference of standard deviation of c), and determine the tissue properties of the specimen based on the calculated values It will be.
  • the values of the weights w a , w b , and w c may be set arbitrarily by the operator, or may be automatically set by the ultrasonic diagnostic apparatus 1.
  • the tissue property determination unit 43 calculates a square root of a value obtained by appropriately weighting and adding the square of the difference between the standard deviations of all feature amounts for each group, and the tissue property corresponding to the group having the minimum square root. May be determined as the tissue property of the specimen. In this case, for example, when the feature amounts are a, b, and c, the tissue property determination unit 43 assigns weights corresponding to the three feature amounts a, b, and c to w ′ a , w ′ b , and w ′ c, respectively.
  • tissue property is determined based on the calculated value.
  • the values of the weights w ′ a , w ′ b , and w ′ c may be arbitrarily set by the operator, or may be automatically set by the ultrasonic diagnostic apparatus 1. May be.
  • the tissue property determination unit 43 determines the tissue property based on a change in standard deviation of each feature amount between the original population and the population to which a new specimen is added. This was just an example.
  • the tissue property determination unit 43 may determine the tissue property based on an average change of each feature amount between the original population and the population to which a new specimen is added.
  • the third embodiment of the present invention is different from the first embodiment in the tissue property determination process in the tissue property determination unit.
  • the configuration of the ultrasonic diagnostic apparatus according to the third embodiment is the same as the configuration of the ultrasonic diagnostic apparatus 1 described in the first embodiment. Therefore, in the following description, components corresponding to those of the ultrasonic diagnostic apparatus 1 are denoted by the same reference numerals.
  • the tissue property determination unit 43 calculates the probability of belonging to each tissue property by using the distance between the average point of the sample in the feature amount space and the average point of the tissue property of the known sample. Specifically, in the case of the feature amount space (b, c) shown in FIG. 9, by using the distances ⁇ , ⁇ , ⁇ between the sample average point Sp and the known sample average points A 0 , B 0 , C 0 , The probability of belonging to each tissue property is calculated. The probability of belonging to each known specimen is set so that the smaller the distance, the larger the probability.
  • the probability belonging to the tissue property A is ⁇ / ⁇
  • the probability belonging to the tissue property B is ⁇ / ⁇
  • the probability belonging to the tissue property C is ⁇ / It can be defined as ⁇ .
  • the information display unit displays the probability of belonging to each tissue property.
  • the fourth embodiment of the present invention differs from the first embodiment in the tissue property determination process in the tissue property determination unit.
  • the configuration of the ultrasonic diagnostic apparatus according to the fourth embodiment is the same as the configuration of the ultrasonic diagnostic apparatus 1 described in the first embodiment. Therefore, in the following description, the same reference numerals are given to the components corresponding to the components of the ultrasonic diagnostic apparatus 1.
  • FIG. 12 is a diagram for explaining the outline of the tissue property determination process performed by the tissue property determination unit 43 in the fourth embodiment.
  • the horizontal axis is the feature quantity b
  • the vertical axis is the feature quantity c.
  • regions are grouped according to organizational properties.
  • the tissue property determination unit 43 determines the tissue property according to the position of the specimen average point.
  • FIG. 12 shows a case where the specimen average point Sp ′ belongs to the group SB ′ (region where the tissue property is B). In this case, the tissue property determination unit 43 determines that the tissue property of the region of interest of the specimen is B.

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

La présente invention concerne un dispositif de diagnostic par ultrasons, qui diagnostique des caractéristiques de tissus chez un sujet sur la base d'ondes ultrasonores reçues, et qui pour cela envoi des ondes ultrasonores à un sujet de diagnostic et reçoit aussi les ondes ultrasonores qui ont été renvoyées par le sujet. Ce dispositif de diagnostic par ultrasons procède par extraction de caractéristiques de fréquences d'un spectre de fréquences dans une zone prescrite chez un sujet. Cette extraction implique une approximation d'un spectre de fréquences résultant de l'analyse de la fréquence des ondes ultrasonores reçues. Le dispositif de diagnostic par ultrasons détermine ensuite les caractéristiques des tissus de la zone prescrite du sujet, utilisant à cet effet ces caractéristiques de fréquence et les caractéristiques de fréquence d'un spectre de fréquences extraites sur la base, d'une part des caractéristiques des tissus d'une pluralité de sujets connus, et d'autre part des ondes ultrasonores renvoyées par la pluralité de sujet connus.
PCT/JP2011/065910 2010-07-20 2011-07-12 Dispositif de diagnostic par ultrasons, procédé et logiciel de mise en œuvre correspondants WO2012011414A1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2012524951A JPWO2012011414A1 (ja) 2010-07-20 2011-07-12 超音波診断装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2010163290 2010-07-20
JP2010-163290 2010-07-20

Publications (1)

Publication Number Publication Date
WO2012011414A1 true WO2012011414A1 (fr) 2012-01-26

Family

ID=45496841

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2011/065910 WO2012011414A1 (fr) 2010-07-20 2011-07-12 Dispositif de diagnostic par ultrasons, procédé et logiciel de mise en œuvre correspondants

Country Status (2)

Country Link
JP (1) JPWO2012011414A1 (fr)
WO (1) WO2012011414A1 (fr)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013162810A (ja) * 2012-02-09 2013-08-22 Ge Medical Systems Global Technology Co Llc 超音波診断装置及びその制御プログラム
WO2015190180A1 (fr) * 2014-06-11 2015-12-17 オリンパス株式会社 Dispositif de diagnostic médical, procédé de fonctionnement d'un dispositif de diagnostic médical et programme de fonctionnement d'un dispositif de diagnostic médical
WO2015198712A1 (fr) * 2014-06-24 2015-12-30 オリンパス株式会社 Appareil d'observation ultrasonore, procédé de fonctionnement d'appareil d'observation ultrasonore, et programme de fonctionnement d'appareil d'observation ultrasonore
WO2017069068A1 (fr) * 2015-10-23 2017-04-27 オリンパス株式会社 Appareil d'observation ultrasonore, procédé de fonctionnement de l'appareil d'observation ultrasonore et programme de fonctionnement pour l'appareil d'observation ultrasonore
US9662093B2 (en) 2014-12-22 2017-05-30 Olympus Corporation Ultrasound observation apparatus, method for operating ultrasound observation apparatus, and computer-readable recording medium
WO2017098931A1 (fr) 2015-12-08 2017-06-15 オリンパス株式会社 Appareil de diagnostic à ultrasons, procédé de fonctionnement de l'appareil de diagnostic à ultrasons et programme de fonctionnement pour l'appareil de diagnostic à ultrasons
JP2022510333A (ja) * 2018-12-07 2022-01-26 ヴェラン メディカル テクノロジーズ，インコーポレイテッド 肺疾患の迅速診断用の気管支内カテーテルシステムおよび方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005536289A (ja) * 2002-08-26 2005-12-02 ザ クリーブランド クリニック ファウンデーション 血管組織を特徴付けするシステムおよび方法
JP2007524431A (ja) * 2003-02-03 2007-08-30 ザ クリーヴランド クリニック ファウンデーション 非侵襲性組織特性顕示システムと方法
JP2009530068A (ja) * 2006-03-22 2009-08-27 ヴォルケイノウ・コーポレーション 分類基準に従った自動プラーク特性決定に基づく自動病変解析
JP2010051553A (ja) * 2008-08-28 2010-03-11 Konica Minolta Medical & Graphic Inc 超音波診断装置および方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6200268B1 (en) * 1999-09-10 2001-03-13 The Cleveland Clinic Foundation Vascular plaque characterization
JP4575737B2 (ja) * 2004-09-29 2010-11-04 富士フイルム株式会社 超音波撮像装置
US8449465B2 (en) * 2005-10-14 2013-05-28 Cleveland Clinic Foundation System and method for characterizing vascular tissue
WO2007072720A1 (fr) * 2005-12-19 2007-06-28 Tohoku University Appareil d’imagerie diagnostique pour une utilisation medicale et procede d’identification de tissu biologique
JP2008200074A (ja) * 2007-02-16 2008-09-04 Hitachi Medical Corp 超音波診断装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005536289A (ja) * 2002-08-26 2005-12-02 ザ クリーブランド クリニック ファウンデーション 血管組織を特徴付けするシステムおよび方法
JP2007524431A (ja) * 2003-02-03 2007-08-30 ザ クリーヴランド クリニック ファウンデーション 非侵襲性組織特性顕示システムと方法
JP2009530068A (ja) * 2006-03-22 2009-08-27 ヴォルケイノウ・コーポレーション 分類基準に従った自動プラーク特性決定に基づく自動病変解析
JP2010051553A (ja) * 2008-08-28 2010-03-11 Konica Minolta Medical & Graphic Inc 超音波診断装置および方法

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013162810A (ja) * 2012-02-09 2013-08-22 Ge Medical Systems Global Technology Co Llc 超音波診断装置及びその制御プログラム
WO2015190180A1 (fr) * 2014-06-11 2015-12-17 オリンパス株式会社 Dispositif de diagnostic médical, procédé de fonctionnement d'un dispositif de diagnostic médical et programme de fonctionnement d'un dispositif de diagnostic médical
JP5897227B1 (ja) * 2014-06-11 2016-03-30 オリンパス株式会社 医用診断装置、医用診断装置の作動方法および医用診断装置の作動プログラム
CN105828726A (zh) * 2014-06-11 2016-08-03 奥林巴斯株式会社 医用诊断装置、医用诊断装置的工作方法以及医用诊断装置的工作程序
US9655593B2 (en) 2014-06-11 2017-05-23 Olympus Corporation Medical diagnostic apparatus, method for operating medical diagnostic apparatus, and computer-readable recording medium
WO2015198712A1 (fr) * 2014-06-24 2015-12-30 オリンパス株式会社 Appareil d'observation ultrasonore, procédé de fonctionnement d'appareil d'observation ultrasonore, et programme de fonctionnement d'appareil d'observation ultrasonore
US9517054B2 (en) 2014-06-24 2016-12-13 Olympus Corporation Ultrasound observation apparatus, method for operating ultrasound observation apparatus, and computer-readable recording medium
US9662093B2 (en) 2014-12-22 2017-05-30 Olympus Corporation Ultrasound observation apparatus, method for operating ultrasound observation apparatus, and computer-readable recording medium
WO2017069068A1 (fr) * 2015-10-23 2017-04-27 オリンパス株式会社 Appareil d'observation ultrasonore, procédé de fonctionnement de l'appareil d'observation ultrasonore et programme de fonctionnement pour l'appareil d'observation ultrasonore
JP6157790B1 (ja) * 2015-10-23 2017-07-05 オリンパス株式会社 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
US10617389B2 (en) 2015-10-23 2020-04-14 Olympus Corporation Ultrasound observation apparatus, method of operating ultrasound observation apparatus, and computer-readable recording medium
WO2017098931A1 (fr) 2015-12-08 2017-06-15 オリンパス株式会社 Appareil de diagnostic à ultrasons, procédé de fonctionnement de l'appareil de diagnostic à ultrasons et programme de fonctionnement pour l'appareil de diagnostic à ultrasons
JP2022510333A (ja) * 2018-12-07 2022-01-26 ヴェラン メディカル テクノロジーズ，インコーポレイテッド 肺疾患の迅速診断用の気管支内カテーテルシステムおよび方法
JP7408658B2 (ja) 2018-12-07 2024-01-05 ヴェラン メディカル テクノロジーズ，インコーポレイテッド 肺疾患の迅速診断用の気管支内カテーテルシステムおよび方法

Also Published As

Publication number Publication date
JPWO2012011414A1 (ja) 2013-09-09

Similar Documents

Publication Publication Date Title
JP5054254B2 (ja) 超音波診断装置、超音波診断装置の作動方法および超音波診断装置の作動プログラム
JP5433097B2 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP5307939B2 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP5568199B1 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP5079177B2 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP5054253B2 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
WO2012063930A1 (fr) Dispositif de diagnostic à ultrasons, procédé de fonctionnement pour dispositif de diagnostic à ultrasons et programme de fonctionnement pour dispositif de diagnostic à ultrasons
JP5430809B1 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP5066306B2 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
WO2012011414A1 (fr) Dispositif de diagnostic par ultrasons, procédé et logiciel de mise en œuvre correspondants

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 11809588

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2012524951

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 11809588

Country of ref document: EP

Kind code of ref document: A1