WO2011099102A1 - 超音波診断装置および内中膜の厚さを計測する方法 - Google Patents

超音波診断装置および内中膜の厚さを計測する方法 Download PDF

Info

Publication number
WO2011099102A1
WO2011099102A1 PCT/JP2010/006658 JP2010006658W WO2011099102A1 WO 2011099102 A1 WO2011099102 A1 WO 2011099102A1 JP 2010006658 W JP2010006658 W JP 2010006658W WO 2011099102 A1 WO2011099102 A1 WO 2011099102A1
Authority
WO
WIPO (PCT)
Prior art keywords
boundary
template
intensity distribution
similarity
echo
Prior art date
Application number
PCT/JP2010/006658
Other languages
English (en)
French (fr)
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 JP2011553664A priority Critical patent/JP5532057B2/ja
Priority to CN201080060241.5A priority patent/CN102695457B/zh
Priority to US13/520,214 priority patent/US8740796B2/en
Publication of WO2011099102A1 publication Critical patent/WO2011099102A1/ja

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
    • A61B8/0858—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
    • 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
    • A61B8/0891—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of blood vessels
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/12—Edge-based segmentation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10132—Ultrasound image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30048—Heart; Cardiac
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30101—Blood vessel; Artery; Vein; Vascular

Definitions

  • the present invention relates to an ultrasonic diagnostic apparatus for calculating the thickness of an intima complex of a blood vessel such as a carotid artery.
  • Diagnosis of arteriosclerosis is performed using an ultrasonic diagnostic apparatus that transmits and receives ultrasonic waves to and from a subject and acquires biological information based on echo signals from the subject.
  • the intima-media complex thickness of the carotid artery is known as an important index for knowing the initial atherosclerosis.
  • “Intima complex” or “intima” is a collective term for the intima and media of the vascular wall that exist between the vascular lumen and the outer membrane.
  • the inner-media complex thickness or inner-media thickness is also called “IMT” (Intima-Media Thickness), and is described as “IMT” in this specification.
  • FIG. 1 shows the intima 3 present between the blood vessel lumen 1 and the outer membrane 2.
  • a lumen boundary 4 that is a boundary between the blood vessel lumen 1 and the intima and an outer membrane boundary 5 that is a boundary between the media and the outer membrane 2 are detected, and the thickness thereof is measured.
  • the technique of Patent Document 1 is known.
  • FIG. 2 shows an echo intensity distribution used when the outer membrane boundary and the lumen boundary are determined in the technique of Patent Document 1.
  • Patent Document 1 in the echo intensity distribution in the depth direction, a position where the intensity change is a predetermined value or more is set as a boundary candidate, a larger one is defined as an outer membrane boundary, and a smaller one is defined as a lumen boundary. It is disclosed.
  • Patent Document 1 since processing is performed uniformly according to whether the intensity change is large or small, there is a case where a correct IMT cannot be measured.
  • FIG. 3 shows an example in which the intensity change at the lumen boundary is larger than the intensity change at the epicardial boundary. This is an example in which the position where the intensity change is large is not always the outer membrane boundary.
  • the present invention has been made in view of the above problems, and its purpose is to accurately determine the lumen boundary and the outer membrane boundary without depending on the magnitude relationship of the echo intensity change between the lumen side and the outer membrane side. It is to detect and realize accurate IMT measurement.
  • the ultrasonic diagnostic apparatus transmits ultrasonic waves to a subject and acquires biological information based on the received echo signal from the subject.
  • the ultrasonic diagnostic apparatus includes an intensity distribution generation unit that generates an echo intensity distribution in a depth direction related to a blood vessel wall from the echo signal, and a template that represents a reference pattern of the echo intensity distribution prepared in advance.
  • a template generation unit that generates a template used for boundary detection, and a pattern similarity calculation unit that calculates the similarity between the template and the intensity distribution, the template being changed in the depth direction
  • a pattern similarity calculation unit that calculates a similarity between the template and the intensity distribution for each change by changing a pattern corresponding to the inner-medium film thickness of the template, and the similarity calculated for each change
  • a boundary determining unit that determines a lumen boundary and an adventitia boundary based on the degree.
  • Another ultrasonic diagnostic apparatus is based on an intensity distribution generation unit that generates an echo intensity distribution in the depth direction related to a blood vessel wall from the echo signal, and a template that represents a reference pattern of the echo intensity distribution prepared in advance.
  • a pattern similarity calculation unit that calculates a similarity between the template and the intensity distribution for each change by changing a pattern corresponding to the inner and middle film thickness of the template, and a plurality of adjacent acoustic lines Boundary continuity for calculating the boundary continuity based on the difference in the intensity value of the echo signal of each acoustic line
  • a calculation unit a boundary evaluation unit that generates an evaluation value obtained by integrating the pattern similarity and the boundary continuity for the combination, and a boundary that determines a lumen boundary and an adventitia boundary based on the integrated evaluation value
  • a determination unit that calculates a similarity between the template and the intensity distribution for each change by changing a pattern corresponding to the inner and middle film thickness of the template, and a plurality of adjacent acoustic lines Boundary continuity for calculating the boundary continuity based on the difference in the intensity value of the echo signal of each acoustic line
  • a calculation unit a boundary evaluation unit that generates an evaluation value obtained by integrating the pattern similarity and the
  • the template generated by the template generation unit may include at least three regions, and the at least three regions may include at least one variable-length region that can change the range of the region.
  • the at least three regions may be arranged in order, and the variable length region may be a central region of the at least three regions.
  • the one end region of the at least three regions may correspond to a low echo tissue region, and the other end region of the at least three regions may correspond to a high echo tissue region.
  • the low-echo tissue region may correspond to a blood flow region, and the high-echo tissue region may correspond to an epicardial wall and body tissue.
  • variable length region may correspond to an intima region of a blood vessel.
  • the pattern similarity calculation unit may calculate a similarity obtained by normalizing the length of the variable length region.
  • the pattern similarity calculation unit may normalize the length of the variable length region by using the average strength of the variable length region as the strength of the variable length region.
  • the template generated by the template generation unit includes two types of sub-templates, and each of the two types of sub-templates includes at least two regions, and the distance between the two types of templates can be changed.
  • the pattern similarity calculation unit may calculate the similarity using each of the two types of sub-templates, and may set a value obtained by combining the two as the similarity.
  • the two regions constituting the two types of templates may correspond to low and high echo tissue regions, respectively.
  • the low echo and high echo tissue regions of the two types of templates correspond to the blood flow region and the intima region, respectively, or to the lumen boundary and the predetermined tissue region, and the predetermined tissue region is May be the adventitia border, the adventitia wall and the body tissue.
  • the template generation unit may change the coefficient value of the template to be generated according to the contrast of the intensity distribution generated by the intensity distribution generation unit.
  • the similarity evaluated by the pattern similarity calculation unit may be defined by an intensity difference between boundaries specified by a template generated by the template generation unit.
  • the similarity evaluated by the pattern similarity calculation unit may be defined by a normalized correlation between the template generated by the template generation unit and the intensity distribution generated by the intensity distribution generation unit.
  • the ultrasonic diagnostic apparatus further includes a rear wall detection unit that detects a region of high echo from the intensity distribution generated by the intensity distribution generation unit, and the boundary determination unit narrows down to the detected rear wall range, May be determined.
  • the high echo area detected by the rear wall detection unit may correspond to a blood vessel wall.
  • the boundary evaluation unit may calculate the evaluation value using dynamic programming.
  • the boundary evaluation unit calculates the evaluation value using an intensity distribution with reduced resolution, and the boundary determination unit provisionally determines a candidate based on the intensity distribution with reduced resolution, and then the original resolution In the intensity distribution, a lumen boundary and an epicardial boundary may be determined from a region around the tentatively determined candidate.
  • the boundary evaluation unit may calculate the evaluation value in the acoustic line direction while narrowing down to a candidate having a good evaluation value or a candidate having a predetermined evaluation value by using dynamic programming.
  • the boundary continuity calculation unit may calculate a boundary continuity defined by a difference in boundary position in the acoustic line direction.
  • the pattern similarity calculation unit, the boundary continuity calculation unit, the boundary evaluation unit, and the boundary determination unit evaluate the possible boundary candidates for two boundaries, a lumen boundary and an epicardial boundary.
  • the boundary determination unit may simultaneously determine the lumen boundary and the epicardial boundary.
  • the boundary continuity calculation unit may calculate a boundary continuity defined by a difference in thickness determined by a position of a lumen boundary and an outer membrane boundary in the acoustic ray direction.
  • the pattern similarity calculation unit calculates the similarity when the other boundary position is not determined , One boundary is fixed, a similarity is calculated for a template having a plurality of thicknesses corresponding to the other boundary position, and the maximum value or the average value of the calculated similarity is It is good also as a similarity of a boundary.
  • the boundary determining unit detects the lumen boundary and then narrows down below the lumen boundary and detects the epicardial boundary, and after detecting the epicardial boundary, narrows down above the epicardial boundary.
  • the evaluation value used for each detection may be compared with the case where the lumen boundary is detected, and a detection boundary having a good evaluation value may be adopted as the final boundary.
  • An IMT calculation unit that calculates the thickness of the intima from the lumen boundary and the epicardial boundary may be further provided.
  • the method according to the present invention is a method for transmitting an ultrasonic wave to a subject and measuring the thickness of the intima of a blood vessel wall based on the received echo signal from the subject, wherein the blood vessel wall is obtained from the echo signal.
  • Generating an echo intensity distribution in the depth direction with respect to the template generating a template used for boundary detection based on a template representing a reference pattern of the echo intensity distribution prepared in advance, the template and the intensity distribution
  • the template is changed in the depth direction, and the pattern corresponding to the inner and middle film thickness of the template is changed, and the template and the intensity distribution are changed for each change.
  • a step of calculating a luminal boundary and an adventitia boundary based on the similarity calculated for each change From the determined said lumen boundary and the adventitia boundary was includes calculating the thickness of the intima-media.
  • the method according to the present invention is a method for transmitting an ultrasonic wave to a subject and measuring the thickness of the intima of a blood vessel wall based on the received echo signal from the subject, wherein the blood vessel wall is obtained from the echo signal.
  • Generating an echo intensity distribution in the depth direction with respect to the template generating a template used for boundary detection based on a template representing a reference pattern of the echo intensity distribution prepared in advance, the template and the intensity distribution For each candidate combination of boundaries existing at different positions in the depth direction, the pattern corresponding to the inner and middle film thickness of the template is changed, and the template is changed for each change. And calculating the similarity between the intensity distribution and the combination determined for a plurality of adjacent acoustic lines.
  • the present invention it is possible to correctly detect the IMT even when the intensity change of the lumen boundary is larger than the intensity change of the epicardial boundary by detecting two boundaries as a set instead of detecting the lumen boundary and the epicardial boundary independently. Can be measured.
  • FIG. 1 is a configuration block diagram of an ultrasonic diagnostic apparatus 10 according to a first embodiment of the present invention. It is a figure which shows the ultrasonic vibration element group 6 incorporated in the ultrasonic probe 11 for transmitting an ultrasonic pulse.
  • FIG. 3 is a block diagram showing a configuration of an IMT measurement unit 16.
  • FIG. It is a flowchart which shows the procedure of IMT calculation processing. It is a figure which shows a coefficient pattern (inner membrane pattern). It is a figure which shows the intensity distribution 90 of the target acoustic line (line), and the template group 91 which changed depth. It is a figure which shows the intensity distribution 90 of the target acoustic line (line), and the template group 92 which changed thickness.
  • FIG. 4 is a diagram illustrating an example of route objects 130a to 130e determined in advance with respect to a target sample 130.
  • (A)-(c) is a figure which shows the example of the cost calculation method using the dynamic programming in 2nd embodiment of this invention.
  • FIG. 1 It is a figure which shows the example of the pair i line of the lumen
  • FIG. 4 is a configuration block diagram of the ultrasonic diagnostic apparatus 10 according to the first embodiment of the present invention.
  • the ultrasonic diagnostic apparatus 10 includes a probe 11, a transmission unit 12, a reception unit 13, a B image generation unit 14, an ROI setting unit 15, an IMT measurement unit 16, a synthesis unit 17, and a display unit 18. I have.
  • the transmission unit 12 generates an ultrasonic pulse and outputs it to the probe 11.
  • the probe 11 transmits the ultrasonic pulse output from the transmission unit 12 into the living body, and receives the echo signal that is reflected back from the living body.
  • the probe 11 outputs the received echo signal to the receiving unit 13.
  • FIG. 5 shows an ultrasonic vibration element group 6 built in the ultrasonic probe 11 for transmitting ultrasonic pulses.
  • the ultrasonic vibration element group 6 (hereinafter, referred to as “vibrator 6”), for example, the ultrasonic vibration elements are arranged along one direction to constitute a so-called 1D array vibrator.
  • the vibrator 6 is made of, for example, a piezoelectric body, and transmits an ultrasonic wave by driving the piezoelectric body, and converts the ultrasonic wave into an electric signal when the piezoelectric body receives the ultrasonic wave.
  • the vibrator 6 can scan each of the predetermined ranges by sequentially driving the respective ultrasonic vibration elements to transmit and receive ultrasonic waves.
  • the transducer 6 can also receive signals reflected at the focal position by superimposing the phases of the ultrasonic waves from the plurality of ultrasonic transducer elements at a predetermined position (focal position). An example of the latter is shown in FIG.
  • FIGS. 6 (a) and 6 (b) schematically show a focused wave of ultrasonic waves when a focal point is formed using a plurality of ultrasonic vibration elements arranged along the x direction.
  • the focused ultrasonic wave has a predetermined width as shown, and has a focal point at a predetermined depth in the z-axis direction.
  • the central axis of the focused wave is also called “acoustic line” as shown in FIGS. 6 (c) and 6 (d).
  • the receiving unit 13 bundles the echo signals output from the probe 11 and outputs them to the B image generation unit 14 and the IMT measurement unit 16.
  • the B image generation unit 14 generates a cross-sectional image from the echo signal, and outputs the cross-sectional image data to the synthesis unit 17.
  • the ROI setting unit 15 receives an instruction for specifying the position of the ROI from the user, and sets a region of interest (ROI) in the position of the artery wall based on the instruction.
  • ROI region of interest
  • the IMT measuring unit 16 measures IMT. More specifically, the IMT measurement unit 16 detects the lumen boundary / outer membrane boundary from the echo signal output from the reception unit 13 and the ROI position information output from the ROI setting unit 15. Further, the IMT measurement unit 16 calculates an IMT value based on the detected positions of the lumen boundary and the outer membrane boundary. The IMT measurement unit 16 outputs information (boundary position information) for specifying the boundary position to the combining unit 17 and the IMT value to the display unit 18.
  • the synthesizing unit 17 synthesizes the IMT value and the B image based on the B image data output from the B image generating unit 14 and the boundary position information output from the IMT measuring unit 16. More specifically, the synthesizing unit 17 superimposes the lumen boundary and the epicardial boundary on the B image, and outputs the obtained synthesized image data.
  • the display unit 18 displays the IMT value output from the IMT measurement unit 16 and the synthesis result output from the synthesis unit 17 on the monitor.
  • the display unit 18 displays the combined image data output from the combining unit 17 on a display device such as a liquid crystal display.
  • FIG. 7 is a block diagram showing a configuration of the IMT measurement unit 16.
  • the IMT measurement unit 16 includes an intensity distribution generation unit 21, a template generation unit 22, a pattern similarity calculation unit 23, a boundary determination unit 24, and an IMT calculation unit 25.
  • the intensity distribution generation unit 21 generates an echo intensity distribution based on the echo signal output from the reception unit 13 and the ROI position information output from the ROI position setting unit 15. Then, the intensity distribution generation unit 21 outputs the generated intensity distribution to the pattern similarity calculation unit 23.
  • This echo intensity distribution is an echo intensity distribution related to a direction (depth direction) from the epidermis of the subject to the inside of the living body.
  • the template generation unit 22 generates a template used for boundary detection based on a template representing a reference pattern of echo intensity distribution prepared and set in advance.
  • the template generation unit 22 outputs a template used for boundary detection to the pattern similarity calculation unit 23.
  • the pattern similarity calculation unit 23 calculates the similarity between the intensity distribution of the received signal output from the intensity distribution generation unit 21 and the template output from the template generation unit 22. This similarity is hereinafter referred to as “pattern similarity”.
  • pattern similarity is hereinafter referred to as “pattern similarity”.
  • the pattern similarity calculation unit 23 outputs the calculated similarity to the boundary determination unit 24.
  • the boundary determination unit 24 determines a lumen boundary and an adventitia boundary based on the pattern similarity output from the pattern similarity calculation unit 23.
  • the boundary determination unit 24 outputs information on the determined boundary position to the synthesis unit 18 and the IMT calculation unit 25.
  • the IMT calculation unit 25 calculates an IMT value based on the determined lumen boundary and outer membrane boundary, and outputs the calculated IMT value to the display unit 17.
  • FIG. 8 is a flowchart showing a procedure of IMT calculation processing.
  • the intensity distribution generation unit 21 generates an echo intensity distribution by performing envelope detection and logarithmic compression on the received echo signal. This is done for each point in the ROI.
  • the intensity distribution generation unit 21 applies, for example, a 3 ⁇ 3 average value filter to the intensity distribution and smoothes it for the purpose of noise reduction.
  • step S13 the pattern similarity calculation unit 23 evaluates the similarity between the intensity distribution in the ROI and the template indicating the intima pattern.
  • the template is generated by the template generation unit 22.
  • the coefficient pattern (inner-media pattern) shown in FIG. 9 is used as a template.
  • a range l 1 corresponds to the blood vessel lumen
  • a range l 2 and a range l 3 correspond to the inner media
  • a range l 4 corresponds to the outer membrane.
  • the template shown in FIG. 9 is a differential type filter, and as a reference coefficient pattern, a range of l 1 and l 3 that takes a negative value in the depth direction in advance and a range of l 2 and l 4 that takes a positive value are shown. Alternately provided.
  • the template is provided with at least three regions, that is, a region (or pattern) corresponding to the blood vessel lumen, a region corresponding to the inner media, and a region corresponding to the outer membrane.
  • the pattern similarity calculation unit 23 evaluates the similarity while moving in the depth direction as shown in FIG. 10 using this template as a detection window, and searches for a suitable depth.
  • the pattern similarity calculation unit 23 also evaluates the similarity to the inner and middle film thicknesses. That is, as shown in FIG. 11, the template generation unit 22 creates a plurality of templates with various changes in inner and middle film thicknesses, calculates the similarity between these templates and the intensity distribution, and searches for a suitable thickness. .
  • FIG. 10 shows an intensity distribution 90 of a target acoustic line (line) and a template group 91 with a changed depth.
  • Each template of the template group 91 is a pattern in one row in the depth direction parallel to the acoustic line.
  • templates of depths p1, p2, p3, p4,..., P10 are shown.
  • only the depth differs in each template shown in figure, The pattern is the same.
  • the pattern similarity calculation unit 23 first calculates the similarity between the intensity distribution 90 of the target acoustic line (line) and the template at the depth p1 and the inner-medium thickness q.
  • the similarity is given by the inner product value of the template T and the intensity distribution X, the intensity distribution at the depth p is X p , the template inner and middle film thickness is q, the template T q corresponding to the inner and middle film thickness q , the template If the size is L, the similarity Y is calculated by the following equation.
  • the pattern of the inner-medium film thickness q included in the template coincides with the pattern of the inner-medium film thickness in the intensity distribution 90 of the target line at the most applicable depth p4.
  • this is an example for convenience of understanding.
  • the pattern of the inner-median film thickness q of the template does not always match the pattern of the inner-median film thickness in the intensity distribution 90 of the target line.
  • the pattern similarity calculation unit 23 calculates the similarity while changing the pattern of the inner and middle film thickness q at the depth p1.
  • FIG. 11 shows an intensity distribution 90 of target acoustic lines (lines) and a template group 92 with a changed thickness.
  • the pattern similarity calculation unit 23 obtains the similarity Y corresponding to each thickness while changing the inner-medium film thickness q.
  • the pattern similarity calculation unit 23 changes the depth to p2 and obtains the similarity Y while changing the inner / middle film thickness q.
  • step S14 of FIG. 8 the boundary determining unit 24 searches for the depth p_max and the inner-media thickness q_max that maximize the similarity Y (p, q), and in step S15, the lumen boundary li and the outer membrane boundary ma. Respectively. That is, when the blood vessel lumen length of the template is l1, it is given by the following equation.
  • the coefficient values a 1 ′, a 2 ′, a 3 ′, a 4 ′ correspond to the corresponding range lengths l 1 , l 2 , l 3 so that the similarity is not biased to a specific inner and middle film thickness.
  • L 4 is recommended. That is, the average intensity is calculated and weighted in each region corresponding to the range lengths l 1 , l 2 , l 3 , and l 4 .
  • the similarity at this time is calculated by the following formula.
  • a 1 ′, a 2 ′, a 3 ′, and a 4 ′ are weighting factors.
  • l 2 when changing the inner and middle film thicknesses, l 2 may be fixed and l 3 may be variable, or both l 2 and l 3 may be variable based on a predetermined ratio.
  • Step S15 is a process for switching to the next line. Thereafter, the above-described processing is performed for each line.
  • step S16 the boundary position is smoothed in the line direction.
  • the similarity is the inner product value between the intensity distribution and the template in the form of a differential filter, but may be a correlation.
  • the coefficient value of the template need not include a negative value.
  • the similarity may be a combined value by preparing sub-templates that match the lumen boundary and the epicardial boundary, calculating the similarity for each, and adding weights.
  • each sub-template patterns corresponding to boundaries and other organizations are defined, but these patterns are variable. That is, the position of the pattern corresponding to the boundary is variable.
  • the boundary region and the region corresponding to the other tissue correspond to the high-echo and low-echo tissue regions, respectively.
  • the coefficient value of the template may be changed based on the rendering degree in the ROI. For example, when the variance of the intensity value is calculated and the variance is large, the contrast of the template is increased, and conversely, the variance is small. May reduce the contrast of the template.
  • the lumen boundary / outer membrane boundary is detected directly, but first, the rear wall is detected using a differential filter or the like, and the inner wall is narrowed to the inside of the rear wall.
  • the outer membrane boundary may be detected. This works effectively when multiple echoes are present in the blood vessel.
  • the template that most closely resembles the intensity distribution of the target line is specified, and the lumen boundary / outer membrane boundary is detected.
  • the template defines a pattern of intensity distribution that includes a lumen boundary and an epicardial boundary. This is not to detect the lumen boundary and epicardial boundary independently, but to detect the lumen boundary and epicardial boundary while considering the intensity change of the lumen boundary and the intensity change of the epicardial boundary. I mean. Therefore, no false detection occurs even when the intensity of the lumen boundary is greater than the epicardial boundary.
  • the difference of this embodiment from the first embodiment is the processing in the IMT measurement unit in this embodiment. Specifically, before smoothing in the line direction, the IMT measurement unit in the present embodiment selects a boundary in consideration of continuity in the line direction in addition to the pattern similarity. When this process is performed, a stable boundary can be detected even when the intensity corresponding to the lumen boundary is not clear.
  • FIG. 4 is used as the configuration of the ultrasonic diagnostic apparatus of the present embodiment, and descriptions other than the IMT measurement unit are omitted.
  • FIG. 13 is a block diagram illustrating a configuration of the IMT measurement unit 30 according to the present embodiment.
  • the IMT measurement unit 30 includes an intensity distribution generation unit 31, a template generation unit 32, a pattern similarity calculation unit 33, a boundary continuity calculation unit 34, a boundary evaluation unit 35, a boundary determination unit 36, and an IMT calculation unit. 37.
  • the intensity distribution generation unit 31 generates an echo intensity distribution based on the echo signal output from the reception unit 13 and the ROI position information output from the ROI position setting unit 15. Then, the intensity distribution generation unit 31 outputs the generated intensity distribution to the pattern similarity calculation unit 23.
  • the template generation unit 32 generates a template used for boundary detection based on a preset template (representing a reference pattern).
  • the template generation unit 32 outputs a template used for boundary detection to the pattern similarity calculation unit 33.
  • the pattern similarity calculation unit 33 receives the received signal intensity distribution output from the intensity distribution generation unit 31 and the template output from the template generation unit 22. Further, the pattern similarity calculation unit 33 receives an instruction regarding the position of the lumen boundary and the position of the epicardial boundary from the boundary evaluation unit 35. The pattern similarity calculation unit 33, based on the position of the lumen boundary and the position of the epicardial boundary specified by the boundary evaluation unit 35, the intensity distribution of the received signal at the position and the template output from the template generation unit 22 The similarity is calculated. The pattern similarity calculation unit 33 outputs the calculated similarity to the boundary evaluation unit 35.
  • the boundary continuity calculation unit 34 receives an instruction regarding the position of the lumen boundary and the position of the epicardial boundary from the boundary evaluation unit 35.
  • the boundary continuity calculation unit 34 evaluates the continuity of the boundary position in the line direction and the thickness using the target line and the adjacent line for the designated lumen boundary position and epicardial boundary position.
  • the boundary continuity calculation unit 34 outputs the evaluation result to the boundary evaluation unit 35 as “continuity”.
  • the boundary evaluation unit 35 calculates lumen boundaries for all combinations of predetermined boundaries. Evaluate the uniqueness and the epithelial boundary. The evaluation value is calculated based on the dynamic programming method described later. The boundary evaluation unit 35 outputs the calculated evaluation value to the boundary determination unit 36.
  • the boundary determination unit 36 determines a lumen boundary and an outer membrane boundary based on the evaluation value output from the boundary evaluation unit 34.
  • the boundary determination unit 36 outputs the determined boundary position to the synthesis unit 18 and the IMT calculation unit 37.
  • the IMT calculation unit 37 calculates an IMT value based on the determined lumen boundary and outer membrane boundary, and outputs the calculated IMT value to the display unit 17. This process is the same as the process of the IMT calculation unit 25 shown in FIG.
  • boundary continuity calculation unit 34 the boundary evaluation unit 35, and the boundary determination unit 36.
  • FIG. 14 is a flowchart showing a procedure of boundary detection processing based on dynamic programming.
  • step S ⁇ b> 21 the boundary evaluation unit 35 initializes the path cost C path for all lumen boundary / epicardial boundary pairs i line .
  • the initial value is given by the similarity C pattern to the inner-media pattern. If the weight is W pattern , the path cost C path is given by the following equation.
  • the similarity C pattern is the same value as in the first embodiment. Note, however, that the values are normalized within the line. However, it does not have to be normalized all over the line. For example, at least a value in a line corresponding to the inner media may be normalized.
  • a pair of lumen boundary / outer membrane boundary corresponds to a combination of lumen boundary / outer membrane boundary.
  • the sample point of ROI is N points, the number of pairs is N squared.
  • the lumen boundary is shallower than the epicardial boundary, or the number of pairs actually used for calculation is less than the N-square pair due to constraints such as the maximum measured thickness determined by the system.
  • step S ⁇ b> 22 the pattern similarity calculation unit 33 calculates the similarity between the inner line pattern and the lumen line / outer membrane boundary pair i line-1 of the previous line.
  • the boundary continuity calculation unit 34 calculates the continuity of the inner median thickness and the boundary position for each lumen boundary / outer membrane boundary pair i line-1 in the previous line.
  • the boundary evaluation unit 35 receives the calculation results and evaluates the continuity of the inner and middle film thicknesses in the line direction and the boundary position.
  • the similarity C pattern with the inner-inner membrane pattern is given by an intensity distribution, a template, and an inner product value, but may be a correlation.
  • the continuity C Thickness of the intima-media thickness is given by the following formula, where the lumen boundary positions of the target line and the previous line are LI line and LI line-1 and the epicardial boundary is MA line and MA line-1 , respectively. It is done.
  • the boundary evaluation unit 35 normalizes the evaluation values C Thickness and C distance in the line, and then integrates them into the cost C total of the route of the selection candidate as in the following equation.
  • the weights W thickness and W distance are weight coefficients for the respective evaluation values.
  • the boundary evaluation unit 35 performs the above calculation for a predetermined lumen target / outer membrane boundary pair i line-1 .
  • FIG. 15 shows route targets 130a to 130e that are determined in advance with respect to the target sample 130 (for example, a lumen boundary selection candidate).
  • the boundary evaluation unit 35 calculates the above-described cost for the route determined by the combination of the target sample 130 and each of the route targets 130a to 130e.
  • FIGS. 16A to 16C show examples of cost calculation methods using dynamic programming in this embodiment.
  • the numerical values described in the route objects 130a to 130e are the accumulated values of route costs.
  • a numerical value written in the vicinity of a line segment connecting the target sample 130 and each of the path targets 130a to 130e is a cost of distance.
  • the numerical value written along the depth direction on the right hand of the target sample 130 is the cost of the pattern.
  • the boundary evaluation unit 35 obtains the minimum value (in FIG. 16B, the sum of the route costs of the route objects 130a to 130e plus the distance cost. “66”) is adopted.
  • the boundary evaluation unit 35 records a route that minimizes the cost.
  • the boundary evaluation unit 35 determines the value of the target sample 130 by adding the cost of the pattern to the above-mentioned minimum value.
  • step S23 a pair of the minimum value of the cost C total and the lumen boundary / outer membrane boundary at that time is searched and selected as a route. At this time, the selected route is recorded and the route cost is updated.
  • the update is given by the following formula.
  • steps S22 and S23 are performed for all lumen boundary / outer membrane boundary pairs i line of the target line.
  • the boundary evaluation unit 35 sends the obtained result to the boundary determination unit 36 as an evaluation value.
  • the boundary determining unit 36 searches for a lumen i / outer membrane boundary pair i line that minimizes the path cost C path .
  • the pair with the minimum path cost C path is the boundary of the rightmost line of the ROI.
  • the boundary of the left line is set by referring to the recorded route.
  • FIG. 17 shows an example of a pair i line of a lumen boundary / outer membrane boundary where the path cost C path is minimized.
  • smoothing in the line direction may be performed as in the first embodiment. Thereby, a smooth boundary line can be obtained.
  • the route search based on dynamic programming is performed on the intensity distribution with reduced resolution (particularly in the depth direction) to temporarily determine candidates, then return to the original resolution,
  • the search may be performed by narrowing down to the peripheral (neighboring) region, and the lumen boundary and the epicardial boundary may be finally determined. Thereby, the processing amount can be reduced while maintaining the detection performance.
  • the route search may be limited to a higher pair with a low route cost or a pair with a route cost of a predetermined value or less. Thereby, the processing amount can be reduced while maintaining the detection performance.
  • the route search is performed by pairing the lumen boundary and the outer membrane boundary, but the route search may be performed separately.
  • FIG. 18 is a flowchart showing a procedure of a process for performing a route search for each of the lumen boundary and the outer membrane boundary.
  • the lumen boundary is first detected, and then the epicardial boundary is detected by narrowing down to a deeper position. Next, the epicardial boundary and the lumen boundary at a position shallower than that are detected in the reverse order. Then, the one with the smaller total route cost is selected.
  • step S31 the boundary evaluation unit 35, the lumen boundary tentatively detected based on dynamic programming, the route cost A 1 is recorded on the basis of the detection result.
  • step S32 the boundary evaluation unit 35 temporarily detects the epicardial boundary based on the dynamic programming for the part below the temporarily detected lumen boundary, and calculates the path cost A 2 based on the detection result. Record.
  • step S34 the boundary evaluation unit 35 temporarily detects the outer membrane boundary based on the dynamic programming, and records the path cost B 1 based on the detection result.
  • step S35 the boundary evaluation unit 35 temporarily detects the lumen boundary based on the dynamic programming for the portion above the temporarily detected outer membrane boundary, and calculates the path cost B 2 based on the detection result. Record.
  • the boundary determination unit 36 determines the epicardial boundary and the lumen boundary using the temporary detection result of the cost A and the cost B with the smaller cost.
  • one boundary is fixed, and the provisional evaluation is performed with a plurality of templates in which the inner and middle film thicknesses are changed by changing the other temporarily detected boundary, and the maximum value of the evaluation value, Alternatively, the average value may be used as the similarity. Thereby, the processing amount can be reduced while maintaining the detection performance.
  • the rear wall may be further detected, and the detection result may be used for determining the boundary.
  • the rear wall is a region existing deeper than the outer membrane boundary 5.
  • FIG. 19 is a block diagram showing the configuration of the IMT measurement unit 16a having the rear wall detection unit 26. Here, it demonstrates as a modification of the IMT measurement part 16 which concerns on 1st embodiment described in FIG.
  • the rear wall detection unit 26 is positioned after the intensity distribution generation unit 21, detects the rear wall using the intensity distribution output from the intensity distribution generation unit 21, and sends the detected rear wall range to the boundary determination unit 24. Output.
  • the rear wall is a blood vessel wall, and exists further in the depth direction than the lumen boundary 4, the intima 3, and the adventitia boundary 5.
  • the rear wall is obtained as a high echo area. Therefore, the rear wall detection unit 26 detects a high echo area from the intensity distribution generated by the intensity distribution generation unit 21. As a result, the boundary determining unit 26 can determine the outer membrane boundary 5 by narrowing it down to the detected rear wall range.
  • the IMT measuring apparatus automatically measures IMT with high accuracy, can reduce the examination time, and is useful, for example, in a screening examination for arteriosclerosis.

Abstract

　内中膜の厚さ計測における内腔境界と外膜境界の検出において、内腔側と外膜側とのエコー強度変化の大小関係に依存しない境界検出方法を提供する。　超音波診断装置は、超音波を被検体に送信し、受信した被検体からのエコー信号に基づいて生体情報を取得する。超音波診断装置は、エコー信号から血管壁に関する深さ方向のエコー強度分布を生成する強度分布生成部と、予め用意された、エコー強度分布の基準パタンを表すテンプレートに基づいて、境界検出に用いるテンプレートを生成するテンプレート生成部と、テンプレートと強度分布との類似度を計算するパタン類似度計算部であって、テンプレートを深さ方向に変化させ、かつ、テンプレートの内中膜厚に対応するパタンを変化させて、変化ごとにテンプレートと強度分布との類似度を計算するパタン類似度計算部と、変化ごとに計算された類似度に基づいて、内腔境界および外膜境界を決定する境界決定部とを備えている。

Description

超音波診断装置および内中膜の厚さを計測する方法
　本発明は、たとえば頸動脈などの血管の内中膜複合体の厚さを計算する超音波診断装置に関する。
　超音波を被検体に送受信し、前記被検体からのエコー信号を基に生体情報を取得する超音波診断装置を用いて、動脈硬化の診断が行われている。動脈硬化の診断において、頸動脈の内中膜複合体厚は初期の粥状硬化を知る重要な指標として知られている。「内中膜複合体」または「内中膜」は血管内腔と外膜との間に存在する、血管壁の内膜および中膜をまとめた呼称である。そして内中膜複合体厚または内中膜厚は「ＩＭＴ」（Ｉｎｔｉｍａ－Ｍｅｄｉａ　Ｔｈｉｃｋｎｅｓｓ）とも呼ばれており、本願明細書では「ＩＭＴ」と記述する。
　図１は、血管内腔１と外膜２との間に存在する内中膜３を示している。検査では、この血管内腔１と内膜との境界である内腔境界４と、中膜と外膜２との境界である外膜境界５とを検出し、その厚さを計測する。このＩＭＴを自動的に計測する方法としては、例えば、特許文献１の技術が知られている。
　図２は、特許文献１の技術において、外膜境界および内腔境界を定める際に利用されるエコー強度分布を示している。特許文献１では、図２に示す通り、深さ方向のエコー強度分布において、強度変化が所定値以上の位置を境界候補とし、大きい方を外膜境界、小さい方を内腔境界と定める方法が開示されている。
特開２００８－１６８０１６号公報
　ところが、特許文献１では、強度変化が大きいか小さいかに応じて画一的に処理が行われているため、正しいＩＭＴを計測できない場合がある。
　たとえば図３は、内腔境界の強度変化の方が、外膜境界の強度変化よりも大きくなる例を示している。強度変化が大きい位置が常に外膜境界とは限らない例である。
　この結果、特許文献１に開示されている方法では、内腔境界を外膜境界として検出することになり、正しいＩＭＴを計測できない。
　本発明は上記課題に鑑みてなされたものであり、その目的は、内腔側と外膜側とのエコー強度変化の大小関係に依存することなく、正確に内腔境界と外膜境界とを検出し、正確なＩＭＴ計測を実現することにある。
　本発明による超音波診断装置は、超音波を被検体に送信し、受信した前記被検体からのエコー信号に基づいて生体情報を取得する。本発明にかかる、ある前記超音波診断装置は、前記エコー信号から血管壁に関する深さ方向のエコー強度分布を生成する強度分布生成部と、予め用意された、エコー強度分布の基準パタンを表すテンプレートに基づいて、境界検出に用いるテンプレートを生成するテンプレート生成部と、前記テンプレートと前記強度分布との類似度を計算するパタン類似度計算部であって、前記テンプレートを前記深さ方向に変化させ、かつ、前記テンプレートの内中膜厚に対応するパタンを変化させて、変化ごとに前記テンプレートと前記強度分布との類似度を計算するパタン類似度計算部と、前記変化ごとに計算された前記類似度に基づいて、内腔境界および外膜境界を決定する境界決定部とを備えている。
　本発明による他の超音波診断装置は、前記エコー信号から血管壁に関する深さ方向のエコー強度分布を生成する強度分布生成部と、予め用意された、エコー強度分布の基準パタンを表すテンプレートに基づいて、境界検出に用いるテンプレートを生成するテンプレート生成部と、前記テンプレートと前記強度分布との類似度を計算するパタン類似度計算部であって、前記深さ方向の異なる位置に存在する境界の各候補の組み合わせについて、前記テンプレートの内中膜厚に対応するパタンを変化させて、変化ごとに前記テンプレートと前記強度分布との類似度を計算するパタン類似度計算部と、隣接する複数の音響線に関して定まる前記組み合わせについて、各音響線のエコー信号の強度値の差分に基づいて境界連続度を計算する境界連続度計算部と、前記組み合わせについて、前記パタン類似度と前記境界連続度とを統合した評価値を生成する境界評価部と、前記統合した評価値に基づき、内腔境界と外膜境界を決定する境界決定部とを備えている。
　前記テンプレート生成部が生成するテンプレートは、少なくとも３つの領域を含み、前記少なくとも３つの領域には、前記領域の範囲を変えることができる可変長領域が少なくとも１つ含まれていてもよい。
　前記少なくとも３つの領域は順に配列されており、前記可変長領域は、前記少なくとも３つの領域のうちの中央の領域であってもよい。
　前記少なくとも３つの領域の一方の端部の領域は、低エコーの組織領域に対応し、前記少なくとも３つの領域の他方の端部の領域は、高エコーの組織領域に対応してもよい。
　前記低エコーの組織領域は血流領域に対応し、前記高エコーの組織領域は血管外膜壁および体組織に対応してもよい。
　前記可変長領域は、血管の内中膜領域に対応してもよい。
　前記パタン類似度計算部は、前記可変長領域の長さを正規化した類似度を計算してもよい。
　前記パタン類似度計算部は、前記可変長領域の平均強度を前記可変長領域の強度として用いることにより、前記可変長領域の長さを正規化してもよい。
　前記テンプレート生成部が生成するテンプレートは、２種類のサブテンプレートから構成されており、前記２種類のサブテンプレートは、それぞれ、少なくとも２つの領域から構成され、２種類のテンプレートの距離を変えることができ、前記パタン類似度計算部は、前記２種類のサブテンプレートを用いて、それぞれ類似度を計算し、両者を合わせた統合した値を類似度としてもよい。
　前記２種類のテンプレートを構成する２つの領域は、低エコーおよび高エコーの組織領域にそれぞれ対応してもよい。
　前記２種類のテンプレートの低エコーおよび高エコーの組織領域は、それぞれ、血流領域および内中膜領域に対応し、または、内腔境界および所定の組織領域に対応し、前記所定の組織領域は、外膜境界、動脈外膜壁および体組織であってもよい。
　前記テンプレート生成部は、前記強度分布生成部が生成した強度分布のコントラストに応じて、生成するテンプレートの係数値を変えてもよい。
　前記パタン類似度計算部が評価する類似度は、前記テンプレート生成部が生成したテンプレートで指定される境界間の強度差で定義されてもよい。
　前記パタン類似度計算部が評価する類似度は、前記テンプレート生成部が生成したテンプレートと前記強度分布生成部が生成した強度分布との正規化相関で定義されてもよい。
　前記超音波診断装置は、前記強度分布生成部が生成した強度分布から高エコーの領域を検出する後壁検出部を更に備え、前記境界決定部は、検出した後壁範囲内に絞って、境界を決定してもよい。
　前記後壁検出部が検出する高エコー領域は、血管壁に対応してもよい。
　前記境界評価部は、動的計画法を用いて前記評価値を計算してもよい。
　前記境界評価部は、解像度を下げた強度分布を用いて前記評価値を計算し、前記境界決定部は、前記解像度を下げた強度分布に基づいて候補を仮決定し、その後、元の解像度の強度分布において、前記仮決定した候補の周辺の領域から、内腔境界と外膜境界とを決定してもよい。
　前記境界評価部は、動的計画法を利用して、前記評価値の良い上位候補、または、所定の評価値を有する候補に絞りながら音響線方向に前記評価値計算してもよい。
　前記境界連続度計算部は、音響線方向の境界位置の差で定義される境界連続度を計算してもよい。
　前記パタン類似度計算部、前記境界連続度計算部、前記境界評価部、および、前記境界決定部は、前記取り得る境界候補を、内腔境界および外膜境界の二つの境界を対象に評価し、前記境界決定部は、前記内腔境界および外膜境界を同時に決定してもよい。
　前記境界連続度計算部は、音響線方向の内腔境界および外膜境界の位置で決まる厚さの差で定義される境界連続度を計算してもよい。
　前記境界決定部が、内腔境界と外膜境界を決められた順序に基づいて決定する際に、前記パタン類似度計算部は、他方の境界位置が決まっていない段階で類似度を計算する場合、一方の境界を固定し、他方の境界位置に対応した複数の厚さを有するテンプレートに対して、類似度を計算し、前記計算した類似度の最大値、又は、平均値を、前記一方の境界の類似度としてもよい。
　前記境界決定部は、内腔境界を検出した後、前記内腔境界の下方に絞って、外膜境界を検出した場合と、外膜境界を検出した後、前記外膜境界の上方に絞って、内腔境界を検出した場合とで、それぞれの検出に用いた評価値を比較し、評価値の良い検出境界を最終的な境界として採用してもよい。
　前記内腔境界と外膜境界から内中膜の厚さを計算するＩＭＴ計算部をさらに備えていてもよい。
　本発明による方法は、超音波を被検体に送信し、受信した前記被検体からのエコー信号に基づいて血管壁の内中膜の厚さを計測する方法であって、前記エコー信号から血管壁に関する深さ方向のエコー強度分布を生成するステップと、予め用意された、エコー強度分布の基準パタンを表すテンプレートに基づいて、境界検出に用いるテンプレートを生成するステップと、前記テンプレートと前記強度分布との類似度を計算するステップであって、前記テンプレートを前記深さ方向に変化させ、かつ、前記テンプレートの内中膜厚に対応するパタンを変化させて、変化ごとに前記テンプレートと前記強度分布との類似度を計算するステップと、前記変化ごとに計算された前記類似度に基づいて、内腔境界および外膜境界を決定するステップと、決定された前記内腔境界および前記外膜境界から、内中膜の厚さを計算するステップとを包含する。
　本発明による方法は、超音波を被検体に送信し、受信した前記被検体からのエコー信号に基づいて血管壁の内中膜の厚さを計測する方法であって、前記エコー信号から血管壁に関する深さ方向のエコー強度分布を生成するステップと、予め用意された、エコー強度分布の基準パタンを表すテンプレートに基づいて、境界検出に用いるテンプレートを生成するステップと、前記テンプレートと前記強度分布との類似度を計算するステップであって、前記深さ方向の異なる位置に存在する境界の各候補の組み合わせについて、前記テンプレートの内中膜厚に対応するパタンを変化させて、変化ごとに前記テンプレートと前記強度分布との類似度を計算するステップと、隣接する複数の音響線に関して定まる前記組み合わせについて、各音響線のエコー信号の強度値の差分に基づいて境界連続度を計算するステップと、前記組み合わせについて、前記パタン類似度と前記境界連続度とを統合した評価値を生成するステップと、前記統合した評価値に基づき、内腔境界と外膜境界を決定するステップと、決定された前記内腔境界および前記外膜境界から、内中膜の厚さを計算するステップとを包含する。
　本発明により、内腔境界・外膜境界を独立に検出するのではなく２つの境界をセットとして検出することで、内腔境界の強度変化が外膜境界の強度変化より大きい場合でも、正しくＩＭＴを計測できる。
血管内腔１と外膜２との間に存在する内中膜３を示す図である。 特許文献１の技術において、外膜境界および内腔境界を定める際に利用されるエコー強度分布を示す図である。 内腔境界の強度変化の方が、外膜境界の強度変化よりも大きくなる例を示す図である。 本発明の第一の実施形態に係る超音波診断装置１０の構成ブロック図である。 超音波パルスを送信するための、超音波プローブ１１に内蔵された超音波振動素子群６を示す図である。 （ａ）および（ｂ）は、ｘ方向に沿って配列された複数の超音波振動子を用いて焦点を形成したときの超音波集束波の模式図であり、（ｃ）および（ｄ）は超音波集束波の簡略図である。 ＩＭＴ計測部１６の構成を示すブロック図である。 ＩＭＴ計算処理の手順を示すフローチャートである。 係数パタン（内中膜パタン）を示す図である。 対象となる音響線（ライン）の強度分布９０と、深さを変えたテンプレート群９１とを示す図である。 対象となる音響線（ライン）の強度分布９０と、厚さを変えたテンプレート群９２とを示す図である。 内腔境界・外膜境界にそれぞれ合ったサブテンプレートを用いて類似度を計算し、重み付け加算等によって類似度を求める例を示す図である。 本発明の第二の実施形態によるＩＭＴ計測部３０の構成を示すブロック図である。 動的計画法に基づく境界検出処理の手順を示すフローチャートである。 対象サンプル１３０に関して予め決められた経路対象１３０ａ～１３０ｅの例を示す図である。 （ａ）～（ｃ）は、本発明の第二の実施形態における動的計画法を用いたコスト計算方法の例を示す図である。 経路コストＣpathが最小となる内腔境界・外膜境界のペアｉlineの例を示す図である。 内腔境界および外膜境界それぞれに関して経路検索を行う処理の手順を示すフローチャートである。 後壁検出部２６を有するＩＭＴ計測部１６ａの構成を示すブロック図である。
　以下、添付の図面を参照しながら、本発明による超音波診断装置の実施形態を説明する。
　（第一の実施形態）
　図４は、本発明の第一の実施形態に係る超音波診断装置１０の構成ブロック図である。
　超音波診断装置１０は、プローブ１１と、送信部１２と、受信部１３と、Ｂ画像生成部１４と、ＲＯＩ設定部１５と、ＩＭＴ計測部１６と、合成部１７と、表示部１８とを備えている。
　送信部１２は超音波パルスを生成し、プローブ１１へ出力する。
　プローブ１１は送信部１２より出力される超音波パルスを生体内へ送信し、生体内から反射して帰ってきたエコー信号を受信する。プローブ１１は受信したエコー信号を受信部１３に出力する。
　図５は、超音波パルスを送信するための、超音波プローブ１１に内蔵された超音波振動素子群６を示す。超音波振動素子群６（以下「振動子６」と記述する。）は、たとえば各超音波振動素子が一方向に沿って並べられており、いわゆる１Ｄアレイ振動子を構成している。
　振動子６は例えば圧電体によって構成され、圧電体を駆動することにより、超音波を送信し、また、超音波を圧電体が受けることによって超音波を電気信号に変換する。振動子６は、各超音波振動素子を順次駆動させて超音波の送受信を行い、所定範囲を走査することができる。また振動子６は、複数の超音波振動素子からの各超音波の位相を所定の位置（焦点位置）において重ね合わせ、焦点位置において反射した信号を受信することもできる。後者の例を図６に示す。
　図６（ａ）および（ｂ）は、ｘ方向に沿って配列された複数の超音波振動素子を用いて焦点を形成したときの超音波の集束波を模式的に示している。超音波の集束波は図示されるような所定の幅を持っており、ｚ軸方向の所定の深さにおいて焦点を有する。ただし、その集束波の中心軸は、図６（ｃ）および（ｄ）において示されるように「音響線」とも呼ばれる。
　再び図４を参照する。受信部１３は、プローブ１１から出力されるエコー信号を束ねて、Ｂ画像生成部１４とＩＭＴ計測部１６へ出力する。
　Ｂ画像生成部１４は、エコー信号から断面画像を生成し、その断面画像のデータを合成部１７に出力する。
　ＲＯＩ設定部１５は、ユーザからＲＯＩの位置を特定する指示を受け取り、その指示に基づいて動脈壁の位置に関心領域（Ｒｅｇｉｏｎ　ｏｆ　Ｉｎｔｅｒｅｓｔ：ＲＯＩ）を設定する。
　ＩＭＴ計測部１６は、ＩＭＴを計測する。より具体的には、ＩＭＴ計測部１６は、受信部１３より出力されるエコー信号とＲＯＩ設定部１５より出力されるＲＯＩ位置の情報とから、内腔境界・外膜境界を検出する。更に、ＩＭＴ計測部１６は検出した内腔境界および外膜境界の位置に基づいてＩＭＴ値を計算する。ＩＭＴ計測部１６は、境界の位置を特定する情報（境界位置情報）を合成部１７へ、ＩＭＴ値を表示部１８へ、それぞれ出力する。
　合成部１７は、Ｂ画像生成部１４より出力されるＢ画像のデータと、ＩＭＴ計測部１６より出力される境界位置情報とに基づいて、ＩＭＴ値とＢ画像とを合成する。より具体的には、合成部１７は、内腔境界および外膜境界をＢ画像に重畳し、得られた合成画像データを出力する。
　表示部１８は、ＩＭＴ計測部１６より出力されるＩＭＴ値と合成部１７より出力される合成結果とをモニタに表示する。
　表示部１８は、合成部１７から出力された合成画像データを、たとえば液晶ディスプレイなどの表示機器に表示する。
　次に、ＩＭＴ計測部１６を詳細に説明する。
　図７は、ＩＭＴ計測部１６の構成を示すブロック図である。
　ＩＭＴ計測部１６は、強度分布生成部２１と、テンプレート生成部２２と、パタン類似度計算部２３と、境界決定部２４と、ＩＭＴ計算部２５とを有している。
　強度分布生成部２１は、受信部１３から出力されたエコー信号、および、ＲＯＩ位置設定部１５から出力されたＲＯＩ位置の情報とに基づいてエコー強度分布を生成する。そして、強度分布生成部２１は、生成した強度分布をパタン類似度計算部２３へ出力する。このエコー強度分布は、被検体の表皮から生体内部への方向（深さ方向）に関するエコー強度分布である。
　テンプレート生成部２２は、あらかじめ用意され設定された、エコー強度分布の基準パタンを表すテンプレートを基に、境界検出に用いるテンプレートを生成する。テンプレート生成部２２は、境界検出に用いるテンプレートをパタン類似度計算部２３に出力する。
　パタン類似度計算部２３は、強度分布生成部２１から出力された、受信信号の強度分布と、テンプレート生成部２２から出力されたテンプレートとの類似度を計算する。この類似度を、以下「パタン類似度」と記述する。パタン類似度計算部２３は、計算した類似度を境界決定部２４に出力する。
　境界決定部２４は、パタン類似度計算部２３から出力されたパタン類似度に基づいて、内腔境界および外膜境界を決定する。境界決定部２４は、決定した境界位置の情報を合成部１８とＩＭＴ計算部２５とに出力する。
　ＩＭＴ計算部２５は、決定された内腔境界および外膜境界に基づいてＩＭＴ値を計算し、計算したＩＭＴ値を、表示部１７へ出力する。
　次に、　本実施形態によるＩＭＴ計算処理の手順を説明する。
　図８は、ＩＭＴ計算処理の手順を示すフローチャートである。
　最初に、ステップＳ１１において、強度分布生成部２１は、受信したエコー信号を包絡線検波、対数圧縮し、エコー強度分布を生成する。これは、ＲＯＩ内の各点について行う。
　更に、ステップＳ１２において、強度分布生成部２１は、ノイズ低減を目的に、強度分布に対して、例えば、３ｘ３の平均値フィルタを適用し、平滑化する。
　続いて処理はパタン類似度計算部２３に移る。なお、以下の処理は、音響線毎（ライン単位）に行われる。
　まず、ステップＳ１３において、パタン類似度計算部２３は、ＲＯＩ中の強度分布と内中膜パタンを示すテンプレートとの類似度を評価する。テンプレートはテンプレート生成部２２によって生成される。
　テンプレートとして、図９に示す係数パタン（内中膜パタン）を用いる。図９において、範囲ｌ1は血管内腔に対応し、範囲ｌ2と範囲ｌ3は内中膜、範囲ｌ4は外膜に対応する。図９に示すテンプレートは、微分型のフィルタであり、基準の係数パタンとして、あらかじめ深さ方向に負値をとるｌ1、ｌ3の範囲と、正値をとるｌ2、ｌ4の範囲を交互に設ける。図示されるように、テンプレートには、少なくとも３つの領域、すなわち、血管内腔に対応する領域（またはパタン）、内中膜に対応する領域、および、外膜に対応する領域が設けられる。
　パタン類似度計算部２３は、このテンプレートを検出窓として、図１０に示す通り、深さ方向に移動させながら、類似度を評価し、当てはまりの良い深さを探索する。
　加えて、パタン類似度計算部２３は、内中膜厚に対する類似度も評価する。すなわち、図１１に示す通り、テンプレート生成部２２で内中膜厚を様々に変えたテンプレートを複数作成し、これらのテンプレートと強度分布との類似度を計算し、当てはまりの良い厚さを探索する。
　以下、より詳しく説明する。
　図１０は、対象となる音響線（ライン）の強度分布９０と、深さを変えたテンプレート群９１とを示している。テンプレート群９１の各テンプレートは、音響線に平行な深さ方向の１列のパターンである。図面には、深さｐ１、ｐ２、ｐ３、ｐ４、・・・、ｐ１０のテンプレートが示されている。なお、図示される各テンプレートは深さのみが異なっており、そのパターンは同じである。
　パタン類似度計算部２３は、まず、対象となる音響線（ライン）の強度分布９０と、深さｐ１で、かつ内中膜厚ｑのときのテンプレートとの類似度を計算する。
　類似度は、テンプレートＴと強度分布Ｘとの内積値で与えられ、深さｐの強度分布をＸp、テンプレートの内中膜厚をｑ、内中膜厚ｑに対応したテンプレートＴq、テンプレートサイズをＬとすると、類似度Ｙは、以下の式によって計算される。
Figure JPOXMLDOC01-appb-M000001
　図１０に示す例では、最も当てはまりのよい深さｐ４において、テンプレートに含まれる内中膜厚ｑのパタンと、対象ラインの強度分布９０における内中膜厚のパタンとが一致していた。しかしながら、これは理解の便宜のための一例である。実際には、最も当てはまりのよい深さにおいて、テンプレートの内中膜厚ｑのパタンと、対象ラインの強度分布９０における内中膜厚のパタンとが一致するとは限らない。
　そこで、本実施形態では、次に、パタン類似度計算部２３は、深さｐ１において、内中膜厚ｑのパターンを変えながら、類似度を計算する。図１１は、対象となる音響線（ライン）の強度分布９０と、厚さを変えたテンプレート群９２とを示している。パタン類似度計算部２３は、内中膜厚ｑを変化させながら、それぞれの厚さに対応する類似度Ｙを求める。
　続いてパタン類似度計算部２３は、深さをｐ２に変更し、内中膜厚ｑを変化させながら類似度Ｙを求める。
　このような処理を深さｐ１０まで継続する。
　図８のステップＳ１４において、境界決定部２４は、類似度Ｙ（ｐ、ｑ）が最大となる深さｐ＿ｍａｘと内中膜厚ｑ＿ｍａｘを探し、ステップＳ１５において、内腔境界ｌｉと外膜境界ｍａをそれぞれ決定する。すなわち、テンプレートの血管内腔長をｌ１とすると、以下の式で与えられる。
Figure JPOXMLDOC01-appb-M000002
　以上が、ライン単位の処理である。
　上述の処理において、特定の内中膜厚に類似度が偏らないよう、係数値ａ1’、ａ2’、ａ3’、ａ4’は、対応する範囲長ｌ1、ｌ2、ｌ3、ｌ4の逆数値に設定すると良い。すなわち、範囲長ｌ1、ｌ2、ｌ3、ｌ4に対応する各領域で、平均強度を計算し、重み付けする形にする。このときの類似度は下記の式によって計算される。
Figure JPOXMLDOC01-appb-M000003
 ここで、ａ1’、ａ2’、ａ3’、ａ4’は、重み係数である。
　更に、内中膜厚を変える場合に、ｌ2を固定、ｌ3を可変、としても良いし、ｌ2、ｌ3を所定の比率に基づき、共に可変にしても良い。
　ステップＳ１５は、次のラインへの切り替え処理である。以下、各ラインに関して上述の処理を行う。
　ライン単位の処理終了後、ステップＳ１６において、境界位置をライン方向に平滑化する。
　以上が、第一の実施形態の動作の流れである。
　なお、上記説明では、類似度を、強度分布と微分型フィルタの形を取るテンプレートとの内積値としたが、相関でも良い。この場合、テンプレートの係数値に負値を含む必要はない。
　また、類似度は、図１２に示す通り、内腔境界・外膜境界にそれぞれ合ったサブテンプレートを用意して、それぞれで類似度を計算し、重み付け加算等により、組み合わせた値としても良い。各サブテンプレートにおいては、境界とそれ以外の組織等に対応するパタンが規定されているが、これらのパタンは可変である。すなわち境界に対応するパタンの位置は可変である。上述の２種類のテンプレートにおける、境界領域とそれ以外の組織等に対応する領域は、それぞれ高エコーおよび低エコーの組織領域にそれぞれ対応している。
　また、テンプレートの係数値は、ＲＯＩ中の描出度に基づいて変更しても良く、例えば、強度値の分散を計算し、分散が大きい場合は、テンプレートのコントラストを強め、逆に分散が小さい場合は、テンプレートのコントラストを弱める、としても良い。
　最後に、上記説明では、直接、内腔境界・外膜境界を検出するとしたが、まず、微分型のフィルタ等を用いて、後壁を検出し、後壁内部に絞って、内腔境界・外膜境界を検出しても良い。これは、血管中に多重エコーが出ている場合に有効に機能する。
　以上、第一の実施形態では、対象ラインの強度分布と最もよく類似するテンプレートを特定して、内腔境界・外膜境界を検出する。テンプレートは、内腔境界および外膜境界を含む強度分布のパターンを規定している。これは、内腔境界・外膜境界を独立して検出するのではなく、内腔境界の強度変化および外膜境界の強度変化を考慮しながら、内腔境界および外膜境界を検出することを意味している。よって、内腔境界の強度が外膜境界より大きい場合においても、誤検出が生じない。
　（第二の実施形態）
　本実施形態が第一の実施形態と異なる点は、本実施形態ではＩＭＴ計測部における処理である。具体的には、本実施形態におけるＩＭＴ計測部は、ライン方向に平滑化する前に、パタン類似度に加え、ライン方向の連続性を考慮して境界を選択する。この処理を行うと、特に、内腔境界に対応する強度が明確でない場合においても、安定した境界を検出可能である。
　本実施形態による超音波診断装置の構成および動作は、ＩＭＴ計測部以外は図４に示す超音波診断装置１０と同じである。したがって、図４を本実施形態の超音波診断装置の構成として援用するとともに、ＩＭＴ計測部以外の説明は省略する。
　図１３は、本実施形態によるＩＭＴ計測部３０の構成を示すブロック図である。
　ＩＭＴ計測部３０は、強度分布生成部３１と、テンプレート生成部３２と、パタン類似度計算部３３と、境界連続度計算部３４と、境界評価部３５と、境界決定部３６と、ＩＭＴ計算部３７とを有している。
　強度分布生成部３１は、受信部１３から出力されたエコー信号、および、ＲＯＩ位置設定部１５から出力されたＲＯＩ位置の情報とに基づいてエコー強度分布を生成する。そして、強度分布生成部３１は、生成した強度分布をパタン類似度計算部２３へ出力する。
　テンプレート生成部３２は、あらかじめ設定したテンプレート（基準パタンを表す）を基に、境界検出に用いるテンプレートを生成する。テンプレート生成部３２は、境界検出に用いるテンプレートをパタン類似度計算部３３に出力する。
　パタン類似度計算部３３は、強度分布生成部３１から出力された、受信信号の強度分布と、テンプレート生成部２２から出力されたテンプレートを受け取る。さらにパタン類似度計算部３３は、境界評価部３５から内腔境界の位置および外膜境界の位置に関する指示を受け取る。パタン類似度計算部３３は、境界評価部３５から指定された内腔境界の位置および外膜境界の位置に基づいて、その位置における受信信号の強度分布と、テンプレート生成部２２から出力されたテンプレートとの類似度を計算する。パタン類似度計算部３３は、計算した類似度を、境界評価部３５へ出力する。
　境界連続度計算部３４は、境界評価部３５から内腔境界の位置および外膜境界の位置に関する指示を受け取る。境界連続度計算部３４は、指定された内腔境界の位置および外膜境界の位置について、対象ラインと隣接ラインとを利用してライン方向の境界位置と厚さの連続性を評価する。境界連続度計算部３４は、評価結果を「連続度」として境界評価部３５へ出力する。
　境界評価部３５は、パタン類似度計算部３３より出力されるパタン類似度と境界連続度計算部３４より出力される境界の連続度に基づき、予め定めた全ての境界の組み合わせについて、内腔境界らしさ、および、外膜境界らしさを評価する。評価値の計算は、後述する動的計画法に基づいて行う。境界評価部３５は、計算した評価値を境界決定部３６へ出力する。
　境界決定部３６は、境界評価部３４より出力される評価値に基づき、内腔境界および外膜境界を決定する。境界決定部３６は、決定した境界位置を合成部１８およびＩＭＴ計算部３７に出力する。
　ＩＭＴ計算部３７は、決定された内腔境界および外膜境界に基づいてＩＭＴ値を計算し、計算したＩＭＴ値を、表示部１７へ出力する。この処理は、図７に示すＩＭＴ計算部２５の処理と同じである。
　次に、本実施形態による境界検出処理の手順を説明する。下記の処理は主として、境界連続度計算部３４、境界評価部３５および境界決定部３６によって行われる。
　図１４は、動的計画法に基づく境界検出処理の手順を示すフローチャートである。
　本実施形態においては、パタン類似性とライン方向の境界の連続性で定義されるコストが最小となるライン方向の経路（＝境界）を探索する形で境界検出処理が実施される。
　ステップＳ２１において、境界評価部３５は、経路コストＣpathを全ての内腔境界・外膜境界のペアｉlineについて初期化する。初期値は、内中膜パタンとの類似度Ｃpatternで与える。重みをＷpatternとすると、経路コストＣpathは下記の式で与えられる。
Figure JPOXMLDOC01-appb-M000004
　ここで、類似度Ｃpatternは、第一の実施形態と同様の値である。しかしながら、ライン内で正規化された値であることに留意されたい。ただし、ライン内全てにおいて正規化されている必要はない。たとえば少なくとも内中膜に対応するライン内の値だけ正規化してもよい。
　「内腔境界・外膜境界のペア」とは、内腔境界・外膜境界の組み合わせに相当する。ＲＯＩのサンプル点をＮ点とすると、ペアの数はＮの二乗存在する。但し、内腔境界は外膜境界より浅い位置にある、または、システムが定める最大計測厚、等の拘束条件により、実際に計算に使うペア数はＮの二乗ペアより少ない。
　次に、２ライン目へ進む。
　そして、ステップＳ２２において、パタン類似度計算部３３は、前ラインの各内腔境界・外膜境界のペアｉline-1に対する内中膜パタンとの類似度を計算する。また境界連続度計算部３４は、前ラインの各内腔境界・外膜境界のペアｉline-1について、内中膜厚と境界位置の連続度をそれぞれ計算する。境界評価部３５は、それらの計算結果を受け取り、ライン方向の内中膜厚と境界位置の連続性を評価する。
　内中膜パタンとの類似度Ｃpatternは、第一の実施形態と同様に、強度分布とテンプレートと内積値で与えられるが、相関でも良い。
　内中膜厚の連続度ＣThicknessは、対象ラインと前ラインの内腔境界位置をそれぞれＬＩline、ＬＩline-1、外膜境界をＭＡline、ＭＡline-1とすると、下記の式で与えられる。
Figure JPOXMLDOC01-appb-M000005
　また境界位置の連続度Ｃdistanceは、下記の式で与えられる。
Figure JPOXMLDOC01-appb-M000006
　境界評価部３５は、評価値ＣThicknessとＣdistanceをライン内でそれぞれ正規化した後、それらを下記の式のように選択候補の経路のコストＣtotalに統合する。
Figure JPOXMLDOC01-appb-M000007
 ここで、重みＷthicknessおよびＷdistanceは、各評価値に対する重み係数である。なお、数７には評価値ＣThicknessとＣdistanceとが含まれているが、境界位置の連続度Ｃdistanceが含まれていればよく、内中膜厚の連続度ＣThicknessはなくてもよい。
　境界評価部３５は、以上の計算を、予め決められた経路対象の内腔境界・外膜境界のペアｉline-1について実施する。
　たとえば図１５では、対象サンプル１３０（たとえば内腔境界の選択候補）に関して予め決められた経路対象１３０ａ～１３０ｅが示されている。境界評価部３５は、対象サンプル１３０と各経路対象１３０ａ～１３０ｅのそれぞれの組み合わせで定まる経路に関し、上述のコストを計算する。
　図１６（ａ）～（ｃ）は、本実施形態における動的計画法を用いたコスト計算方法の例を示す。図１６（ａ）において、経路対象１３０ａ～１３０ｅの内部に記載されている数値は、経路コストの累積値である。対象サンプル１３０と、各経路対象１３０ａ～１３０ｅとを結ぶ線分の近傍に記載されている数値は、距離のコストである。そして、対象サンプル１３０の右手に深さ方向に沿って記載されている数値は、パタンのコストである。
　図１６（ｂ）に示されるように、境界評価部３５は、各経路対象１３０ａ～１３０ｅの経路コストの累積値に距離のコストを加えたもののうちの、最小の値（図１６（ｂ）では「６６」）を与える経路を採用する。境界評価部３５は、コストが最小となる経路を記録しておく。
　さらに境界評価部３５は、上述の最小の値にパタンのコストを加算して、対象サンプル１３０の値として決定する。
　再び図１４を参照する。ステップＳ２３において、コストＣtotalの最小値とそのときの内腔境界・外膜境界のペアを探索、経路として選択する。このとき、選択した経路を記録しておくと共に、経路コストを更新する。更新は、下記の式であたえられる。
Figure JPOXMLDOC01-appb-M000008
　これら、ステップＳ２２とＳ２３の処理を対象ラインの全ての内腔境界・外膜境界のペアｉlineについて実施する。境界評価部３５は、得られた結果を評価値として境界決定部３６に送る。
　全ラインについて計算が完了した後、ステップＳ２６において境界決定部３６は、経路コストＣpathが最小となる内腔境界・外膜境界のペアｉlineを探す。ここで、経路コストＣpathが最小となるペアがＲＯＩの右端ラインの境界となる。以降、左側のラインの境界については、記録しておいた、経路を参照していくことで、設定される。
　図１７は、経路コストＣpathが最小となる内腔境界・外膜境界のペアｉlineの例を示す。
　以上が、処理の流れである。
　尚、この後、第一の実施形態と同様に、ライン方向の平滑化を実施しても良い。これにより、滑らかな境界線を得ることができる。
　また、動的計画法に基づく経路探索は、解像度を落とした（特に、深さ方向）強度分布に対して実施して候補を仮決定した後、元の解像度に戻して、仮決定した候補の周辺（近傍）の領域に絞って探索を行い、内腔境界と外膜境界とを最終的に決定してもよい。これにより、検出性能を維持しつつ、処理量を削減できる。
　また、動的計画法に基づき、経路探索していく際に、経路コストが小さい上位ペア、または、経路コストが所定値以下のペア、に絞って経路探索しても良い。これにより、検出性能を維持しつつ、処理量を削減できる。
　最後に、上記説明では、内腔境界・外膜境界をペアにして経路探索しているが、別々に経路探索しても良い。
　たとえば図１８は、内腔境界および外膜境界それぞれに関して経路検索を行う処理の手順を示すフローチャートである。この処理では、まず先に内腔境界を検出し、その後、それより深い位置に絞って外膜境界を検出する。次に、その逆の順序で外膜境界およびそれより浅い位置の内腔境界を検出する。そして、トータルの経路コストが小さい方を選択する。
　ステップＳ３１において、境界評価部３５は、動的計画法に基づいて内腔境界を仮検出し、検出結果に基づいて経路コストＡ1を記録する。次のステップＳ３２において、境界評価部３５は、仮検出した内腔境界よりも下方を対象に、動的計画法に基づいて外膜境界を仮検出し、検出結果に基づいて経路コストＡ2を記録する。ステップＳ３３において、境界評価部３５は、内腔境界および外膜境界の順に仮検出したときの経路コストＡ（＝Ａ1＋Ａ2）を計算し記録する。
　次に、ステップＳ３４では、境界評価部３５は、動的計画法に基づいて外膜境界を仮検出し、検出結果に基づいて経路コストＢ1を記録する。次のステップＳ３５において、境界評価部３５は、仮検出した外膜境界よりも上方を対象に、動的計画法に基づいて内腔境界を仮検出し、検出結果に基づいて経路コストＢ2を記録する。ステップＳ３３において、境界評価部３５は、外膜境界および内腔境界の順に仮検出したときの経路コストＢ（＝Ｂ1＋Ｂ2）を計算し記録する。
　そして、境界決定部３６は、コストＡおよびコストＢのうち、コストの小さい方の仮検出結果を用いて、外膜境界および内腔境界を決定する。
　但し、一方の境界を検出する場合、他方の境界は決まっていない。そのため、図１１に示すように、一方の境界を固定し、仮検出した他方の境界を変化させることで内中膜厚を変えた複数のテンプレートで当てはまりを評価し、その評価値の最大値、又は、平均値を類似度としても良い。これにより、検出性能を維持しつつ、処理量を削減できる。
　以上、第二の実施形態では、ライン方向の境界連続性を考慮するため、特に、内腔境界に対応する強度が弱い場合においても、安定した境界を検出できる。
　上述の第一および第二の実施形態においては、さらに後壁を検出し、その検出結果を境界の決定に利用してもよい。後壁は、図１に示されるように、外膜境界５よりも深い位置に存在する領域である。
　図１９は、後壁検出部２６を有するＩＭＴ計測部１６ａの構成を示すブロック図である。ここでは、図７に記載された第一の実施形態に係るＩＭＴ計測部１６の変形例として説明する。
　後壁検出部２６は、強度分布生成部２１の後段に位置し、強度分布生成部２１より出力される強度分布を用いて、後壁を検出し、検出した後壁範囲を境界決定部２４へ出力する。
　より具体的に説明する。図１に示されるように、後壁は血管壁であり、深さ方向にみて、内腔境界４、内中膜３、外膜境界５よりもさらに奥に存在する。そして、後壁は高エコーの領域として得られる。そこで後壁検出部２６は、強度分布生成部２１が生成した強度分布から高エコーの領域を検出する。その結果、境界決定部２６は、検出した後壁範囲内に絞って、外膜境界５を決定することができる。
　本発明に係るＩＭＴ計測装置は、精度の良いＩＭＴを自動的に計測するものであり、検査時間を短縮でき、例えば、動脈硬化のスクリーニング検査で有用である。
　１１　プローブ
　１２　送信部
　１３　受信部
　１４　Ｂ画像生成部
　１５　ＲＯＩ設定部
　１６　ＩＭＴ計測部
　１７　合成部
　１８　表示部
　２１　強度分布生成部
　２２　テンプレート生成部
　２３　パタン類似度計算部
　２４　境界決定部
　２５　ＩＭＴ計算部
　２６　後壁検出部
　２７　コントラスト計算部
　２８　検出窓
　３１　強度分布生成部
　３２　テンプレート生成部
　３３　パタン類似度計算部
　３４　境界連続度計算部
　３５　境界評価部
　３６　境界決定部
　３７　ＩＭＴ計算部

Claims (28)

  1. 　超音波を被検体に送信し、受信した前記被検体からのエコー信号に基づいて生体情報を取得する超音波診断装置であって、
    　前記エコー信号から血管壁に関する深さ方向のエコー強度分布を生成する強度分布生成部と、
    　予め用意された、エコー強度分布の基準パタンを表すテンプレートに基づいて、境界検出に用いるテンプレートを生成するテンプレート生成部と、
    　前記テンプレートと前記強度分布との類似度を計算するパタン類似度計算部であって、前記テンプレートを前記深さ方向に変化させ、かつ、前記テンプレートの内中膜厚に対応するパタンを変化させて、変化ごとに前記テンプレートと前記強度分布との類似度を計算するパタン類似度計算部と、
    　前記変化ごとに計算された前記類似度に基づいて、内腔境界および外膜境界を決定する境界決定部と
    　を備えた、超音波診断装置。
  2. 　超音波を被検体に送信し、受信した前記被検体からのエコー信号に基づいて生体情報を取得する超音波診断装置であって、
    　前記エコー信号から血管壁に関する深さ方向のエコー強度分布を生成する強度分布生成部と、
    　予め用意された、エコー強度分布の基準パタンを表すテンプレートに基づいて、境界検出に用いるテンプレートを生成するテンプレート生成部と、
    　前記テンプレートと前記強度分布との類似度を計算するパタン類似度計算部であって、前記深さ方向の異なる位置に存在する境界の各候補の組み合わせについて、前記テンプレートの内中膜厚に対応するパタンを変化させて、変化ごとに前記テンプレートと前記強度分布との類似度を計算するパタン類似度計算部と、
    　隣接する複数の音響線に関して定まる前記組み合わせについて、各音響線のエコー信号の強度値の差分に基づいて境界連続度を計算する境界連続度計算部と、
    　前記組み合わせについて、前記パタン類似度と前記境界連続度とを統合した評価値を生成する境界評価部と、
    　前記統合した評価値に基づき、内腔境界と外膜境界を決定する境界決定部と
    　を備えた、超音波診断装置。
  3. 　前記テンプレート生成部が生成するテンプレートは、少なくとも３つの領域を含み、前記少なくとも３つの領域には、前記領域の範囲を変えることができる可変長領域が少なくとも１つ含まれる、請求項１または請求項２に記載の超音波診断装置。
  4. 　前記少なくとも３つの領域は順に配列されており、
    　前記可変長領域は、前記少なくとも３つの領域のうちの中央の領域である、請求項３に記載の超音波診断装置。
  5.  前記少なくとも３つの領域の一方の端部の領域は、低エコーの組織領域に対応し、前記少なくとも３つの領域の他方の端部の領域は、高エコーの組織領域に対応する、請求項３に記載の超音波診断装置。
  6. 　前記低エコーの組織領域は血流領域に対応し、前記高エコーの組織領域は血管外膜壁および体組織に対応する、請求項５に記載の超音波診断装置。
  7. 　前記可変長領域は、血管の内中膜領域に対応する、請求項４に記載の超音波診断装置。
  8. 　前記パタン類似度計算部は、前記可変長領域の長さを正規化した類似度を計算する、請求項３に記載の超音波診断装置。
  9. 　前記パタン類似度計算部は、前記可変長領域の平均強度を前記可変長領域の強度として用いることにより、前記可変長領域の長さを正規化する、請求項８に記載の超音波診断装置。
  10. 　前記テンプレート生成部が生成するテンプレートは、２種類のサブテンプレートから構成されており、前記２種類のサブテンプレートは、それぞれ、少なくとも２つの領域から構成され、２種類のテンプレートの距離を変えることができ、
    　前記パタン類似度計算部は、前記２種類のサブテンプレートを用いて、それぞれ類似度を計算し、両者を合わせた統合した値を類似度とする、請求項１または請求項２に記載の超音波診断装置。
  11. 　前記２種類のテンプレートを構成する２つの領域は、低エコーおよび高エコーの組織領域にそれぞれ対応する、請求項１０に記載の超音波診断装置。
  12. 　前記２種類のテンプレートの低エコーおよび高エコーの組織領域は、それぞれ、血流領域および内中膜領域に対応し、または、内腔境界および所定の組織領域に対応し、前記所定の組織領域は、外膜境界、動脈外膜壁および体組織である、請求項１０に記載の超音波診断装置。
  13. 　前記テンプレート生成部は、前記強度分布生成部が生成した強度分布のコントラストに応じて、生成するテンプレートの係数値を変える、請求項１または請求項２記載の超音波診断装置。
  14. 　前記パタン類似度計算部が評価する類似度は、前記テンプレート生成部が生成したテンプレートで指定される境界間の強度差で定義される、請求項１または請求項２に記載の超音波診断装置。
  15. 　前記パタン類似度計算部が評価する類似度は、前記テンプレート生成部が生成したテンプレートと前記強度分布生成部が生成した強度分布との正規化相関で定義される、請求項１または請求項２に記載の超音波診断装置。
  16. 　前記強度分布生成部が生成した強度分布から高エコーの領域を検出する後壁検出部を更に備え、前記境界決定部は、検出した後壁範囲内に絞って、境界を決定する、請求項１または請求項２に記載の超音波診断装置。
  17. 　前記後壁検出部が検出する高エコー領域は、血管壁に対応する、請求項１６に記載の超音波診断装置。
  18. 　前記境界評価部は、動的計画法を用いて前記評価値を計算する、請求項２に記載の超音波診断装置。
  19. 　前記境界評価部は、解像度を下げた強度分布を用いて前記評価値を計算し、
    　前記境界決定部は、前記解像度を下げた強度分布に基づいて候補を仮決定し、その後、元の解像度の強度分布において、前記仮決定した候補の周辺の領域から、内腔境界と外膜境界とを決定する、請求項１８に記載の超音波診断装置。
  20. 　前記境界評価部は、動的計画法を利用して、前記評価値の良い上位候補、または、所定の評価値を有する候補に絞りながら音響線方向に前記評価値計算する、請求項１８に記載の超音波診断装置。
  21. 　前記境界連続度計算部は、音響線方向の境界位置の差で定義される境界連続度を計算する、請求項１８に記載の超音波診断装置。
  22. 　前記パタン類似度計算部、前記境界連続度計算部、前記境界評価部、および、前記境界決定部は、前記取り得る境界候補を、内腔境界および外膜境界の二つの境界を対象に評価し、
    　前記境界決定部は、前記内腔境界および外膜境界を同時に決定する、請求項１８に記載の超音波診断装置。
  23. 　前記境界連続度計算部は、音響線方向の内腔境界および外膜境界の位置で決まる厚さの差で定義される境界連続度を計算する、請求項１８に記載の超音波診断装置。
  24. 　前記境界決定部が、内腔境界と外膜境界を決められた順序に基づいて決定する際に、
    　前記パタン類似度計算部は、他方の境界位置が決まっていない段階で類似度を計算する場合、一方の境界を固定し、他方の境界位置に対応した複数の厚さを有するテンプレートに対して、類似度を計算し、前記計算した類似度の最大値、又は、平均値を、前記一方の境界の類似度とする、請求項１８に記載の超音波診断装置。
  25. 　前記境界決定部は、内腔境界を検出した後、前記内腔境界の下方に絞って、外膜境界を検出した場合と、外膜境界を検出した後、前記外膜境界の上方に絞って、内腔境界を検出した場合とで、それぞれの検出に用いた評価値を比較し、評価値の良い検出境界を最終的な境界として採用する、請求項１８に記載の超音波診断装置。
  26. 　前記内腔境界と外膜境界から内中膜の厚さを計算するＩＭＴ計算部をさらに備えた、請求項１から２５のいずれかに記載の超音波診断装置。
  27. 　超音波を被検体に送信し、受信した前記被検体からのエコー信号に基づいて血管壁の内中膜の厚さを計測する方法であって、
    　前記エコー信号から血管壁に関する深さ方向のエコー強度分布を生成するステップと、
    　予め用意された、エコー強度分布の基準パタンを表すテンプレートに基づいて、境界検出に用いるテンプレートを生成するステップと、
    　前記テンプレートと前記強度分布との類似度を計算するステップであって、前記テンプレートを前記深さ方向に変化させ、かつ、前記テンプレートの内中膜厚に対応するパタンを変化させて、変化ごとに前記テンプレートと前記強度分布との類似度を計算するステップと、
    　前記変化ごとに計算された前記類似度に基づいて、内腔境界および外膜境界を決定するステップと、
    　決定された前記内腔境界および前記外膜境界から、内中膜の厚さを計算するステップと
    　を包含する、内中膜の厚さを計測する方法。
  28. 　超音波を被検体に送信し、受信した前記被検体からのエコー信号に基づいて血管壁の内中膜の厚さを計測する方法であって、
    　前記エコー信号から血管壁に関する深さ方向のエコー強度分布を生成するステップと、
    　予め用意された、エコー強度分布の基準パタンを表すテンプレートに基づいて、境界検出に用いるテンプレートを生成するステップと、
    　前記テンプレートと前記強度分布との類似度を計算するステップであって、前記深さ方向の異なる位置に存在する境界の各候補の組み合わせについて、前記テンプレートの内中膜厚に対応するパタンを変化させて、変化ごとに前記テンプレートと前記強度分布との類似度を計算するステップと、
    　隣接する複数の音響線に関して定まる前記組み合わせについて、各音響線のエコー信号の強度値の差分に基づいて境界連続度を計算するステップと、
    　前記組み合わせについて、前記パタン類似度と前記境界連続度とを統合した評価値を生成するステップと、
    　前記統合した評価値に基づき、内腔境界と外膜境界を決定するステップと、
    　決定された前記内腔境界および前記外膜境界から、内中膜の厚さを計算するステップと
    　を包含する、内中膜の厚さを計測する方法。
PCT/JP2010/006658 2010-02-10 2010-11-12 超音波診断装置および内中膜の厚さを計測する方法 WO2011099102A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2011553664A JP5532057B2 (ja) 2010-02-10 2010-11-12 超音波診断装置および内中膜の厚さを計測する方法
CN201080060241.5A CN102695457B (zh) 2010-02-10 2010-11-12 超声波诊断装置及测量内中膜的厚度的方法
US13/520,214 US8740796B2 (en) 2010-02-10 2010-11-12 Ultrasonic diagnostic device, and method for measuring intima-media thickness

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2010-027245 2010-02-10
JP2010027245 2010-02-10

Publications (1)

Publication Number Publication Date
WO2011099102A1 true WO2011099102A1 (ja) 2011-08-18

Family

ID=44367416

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2010/006658 WO2011099102A1 (ja) 2010-02-10 2010-11-12 超音波診断装置および内中膜の厚さを計測する方法

Country Status (4)

Country Link
US (1) US8740796B2 (ja)
JP (1) JP5532057B2 (ja)
CN (1) CN102695457B (ja)
WO (1) WO2011099102A1 (ja)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013161228A1 (ja) * 2012-04-23 2013-10-31 パナソニック株式会社 超音波診断装置および超音波診断装置の制御方法
WO2014034148A1 (ja) * 2012-09-03 2014-03-06 パナソニック株式会社 超音波診断装置、超音波診断装置の制御方法および超音波診断装置の制御器
EP2848201A1 (en) 2013-09-13 2015-03-18 Fujifilm Corporation Ultrasound diagnostic apparatus of producing ultrasound image
JP2015073584A (ja) * 2013-10-07 2015-04-20 コニカミノルタ株式会社 超音波診断装置、超音波診断装置の制御器及び超音波診断装置の制御方法
JP2015104477A (ja) * 2013-11-29 2015-06-08 セイコーエプソン株式会社 超音波測定装置および超音波測定方法
JP2016013472A (ja) * 2015-09-17 2016-01-28 富士フイルム株式会社 超音波画像生成方法
WO2019187649A1 (ja) 2018-03-28 2019-10-03 富士フイルム株式会社 超音波診断装置および超音波診断装置の制御方法
WO2019187647A1 (ja) 2018-03-28 2019-10-03 富士フイルム株式会社 超音波診断装置および超音波診断装置の制御方法
JP2021041029A (ja) * 2019-09-12 2021-03-18 テルモ株式会社 診断支援装置、診断支援システム、及び診断支援方法

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9592027B2 (en) * 2013-03-14 2017-03-14 Volcano Corporation System and method of adventitial tissue characterization
US20150216500A1 (en) * 2014-01-31 2015-08-06 Seiko Epson Corporation Ultrasonic measurement apparatus and ultrasonic measurement method
CN107106125B (zh) * 2014-12-22 2024-02-13 皇家飞利浦有限公司 用于测量动脉参数的系统和方法
CN109688939A (zh) 2016-09-12 2019-04-26 富士胶片株式会社 超声波诊断系统及超声波诊断系统的控制方法
EP3749215A4 (en) * 2018-02-07 2021-12-01 Atherosys, Inc. DEVICE AND METHOD FOR CONTROLLING THE ULTRASONIC RECORDING OF THE PERIPHERAL ARTERIES IN THE TRANSVERSAL PLANE
JP7199972B2 (ja) 2019-01-11 2023-01-06 キヤノン株式会社 情報処理装置、情報処理方法、プログラム
CN110680398A (zh) * 2019-10-11 2020-01-14 深圳开立生物医疗科技股份有限公司 一种超声回波成像方法、系统及超声设备

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11327A (ja) * 1997-06-11 1999-01-06 Toshiba Iyou Syst Eng Kk 超音波診断装置
JPH11318896A (ja) * 1998-05-18 1999-11-24 Masao Ito 血管膜厚測定装置及び動脈硬化診断装置
WO2005002446A1 (ja) * 2003-07-03 2005-01-13 Matsushita Electric Industrial Co., Ltd. 超音波診断装置
JP2005534462A (ja) * 2002-08-02 2005-11-17 ダイアグノスティック・アルトラサウンド・コーポレーション 液体が充填した構造物および液体が充填していない構造物の非侵襲性測定用三次元超音波ベース機器
JP2007512862A (ja) * 2003-11-13 2007-05-24 サントル・オスピタリエ・ドゥ・リュニヴェルシテ・ドゥ・モントリオール 自動多次元血管内超音波画像セグメンテーションの方法
JP2008006188A (ja) * 2006-06-30 2008-01-17 Fujifilm Corp 医用画像表示処理装置、及び、医用画像表示処理プログラム
WO2008023618A1 (fr) * 2006-08-21 2008-02-28 Tohoku University échographe
JP2008161220A (ja) * 2006-12-26 2008-07-17 Hitachi Medical Corp 医用画像診断装置
JP2008168016A (ja) * 2007-01-15 2008-07-24 Fujifilm Corp 超音波診断装置、ｉｍｔ計測方法及びｉｍｔ計測プログラム

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7611466B2 (en) 2002-06-07 2009-11-03 Verathon Inc. Ultrasound system and method for measuring bladder wall thickness and mass
JP2004305377A (ja) * 2003-04-04 2004-11-04 Shimadzu Corp 超音波診断装置
US7727153B2 (en) * 2003-04-07 2010-06-01 Sonosite, Inc. Ultrasonic blood vessel measurement apparatus and method
WO2008001928A1 (fr) 2006-06-30 2008-01-03 Fujifilm Corporation Appareil d'affichage d'image médicale et programme d'affichage d'image médicale
EP2072013A4 (en) * 2006-10-10 2014-12-03 Hitachi Medical Corp MEDICAL PICTURE DIAGNOSTIC DEVICE, MEDICAL IMAGING METHOD AND MEDICAL PICTURE MEASUREMENT PROGRAM
JP5166154B2 (ja) * 2008-07-18 2013-03-21 富士フイルム株式会社 超音波診断装置

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11327A (ja) * 1997-06-11 1999-01-06 Toshiba Iyou Syst Eng Kk 超音波診断装置
JPH11318896A (ja) * 1998-05-18 1999-11-24 Masao Ito 血管膜厚測定装置及び動脈硬化診断装置
JP2005534462A (ja) * 2002-08-02 2005-11-17 ダイアグノスティック・アルトラサウンド・コーポレーション 液体が充填した構造物および液体が充填していない構造物の非侵襲性測定用三次元超音波ベース機器
WO2005002446A1 (ja) * 2003-07-03 2005-01-13 Matsushita Electric Industrial Co., Ltd. 超音波診断装置
JP2007512862A (ja) * 2003-11-13 2007-05-24 サントル・オスピタリエ・ドゥ・リュニヴェルシテ・ドゥ・モントリオール 自動多次元血管内超音波画像セグメンテーションの方法
JP2008006188A (ja) * 2006-06-30 2008-01-17 Fujifilm Corp 医用画像表示処理装置、及び、医用画像表示処理プログラム
WO2008023618A1 (fr) * 2006-08-21 2008-02-28 Tohoku University échographe
JP2008161220A (ja) * 2006-12-26 2008-07-17 Hitachi Medical Corp 医用画像診断装置
JP2008168016A (ja) * 2007-01-15 2008-07-24 Fujifilm Corp 超音波診断装置、ｉｍｔ計測方法及びｉｍｔ計測プログラム

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013161228A1 (ja) * 2012-04-23 2013-10-31 パナソニック株式会社 超音波診断装置および超音波診断装置の制御方法
JPWO2014034148A1 (ja) * 2012-09-03 2016-08-08 コニカミノルタ株式会社 超音波診断装置、超音波診断装置の制御方法および超音波診断装置の制御器
WO2014034148A1 (ja) * 2012-09-03 2014-03-06 パナソニック株式会社 超音波診断装置、超音波診断装置の制御方法および超音波診断装置の制御器
EP2848201A1 (en) 2013-09-13 2015-03-18 Fujifilm Corporation Ultrasound diagnostic apparatus of producing ultrasound image
JP2015054179A (ja) * 2013-09-13 2015-03-23 富士フイルム株式会社 超音波診断装置および超音波画像生成方法
US10085714B2 (en) 2013-09-13 2018-10-02 Fujifilm Corporation Ultrasound diagnostic apparatus and method of producing ultrasound image
US9538988B2 (en) 2013-09-13 2017-01-10 Fujifilm Corporation Ultrasound diagnostic apparatus and method of producing ultrasound image
JP2015073584A (ja) * 2013-10-07 2015-04-20 コニカミノルタ株式会社 超音波診断装置、超音波診断装置の制御器及び超音波診断装置の制御方法
JP2015104477A (ja) * 2013-11-29 2015-06-08 セイコーエプソン株式会社 超音波測定装置および超音波測定方法
JP2016013472A (ja) * 2015-09-17 2016-01-28 富士フイルム株式会社 超音波画像生成方法
WO2019187649A1 (ja) 2018-03-28 2019-10-03 富士フイルム株式会社 超音波診断装置および超音波診断装置の制御方法
WO2019187647A1 (ja) 2018-03-28 2019-10-03 富士フイルム株式会社 超音波診断装置および超音波診断装置の制御方法
EP4238504A1 (en) 2018-03-28 2023-09-06 FUJI-FILM Corporation Ultrasound diagnostic apparatus and control method of ultrasound diagnostic apparatus
JP2021041029A (ja) * 2019-09-12 2021-03-18 テルモ株式会社 診断支援装置、診断支援システム、及び診断支援方法
WO2021048834A1 (ja) * 2019-09-12 2021-03-18 テルモ株式会社 診断支援装置、診断支援システム、及び診断支援方法
JP7300352B2 (ja) 2019-09-12 2023-06-29 テルモ株式会社 診断支援装置、診断支援システム、及び診断支援方法
US11957509B2 (en) 2019-09-12 2024-04-16 Terumo Kabushiki Kaisha Diagnosis assistance device, diagnosis assistance system, and diagnosis assistance method

Also Published As

Publication number Publication date
US20120310086A1 (en) 2012-12-06
JPWO2011099102A1 (ja) 2013-06-13
CN102695457B (zh) 2014-12-31
CN102695457A (zh) 2012-09-26
US8740796B2 (en) 2014-06-03
JP5532057B2 (ja) 2014-06-25

Similar Documents

Publication Publication Date Title
JP5532057B2 (ja) 超音波診断装置および内中膜の厚さを計測する方法
CN101505664B (zh) 超声波诊断装置
JPWO2006082966A1 (ja) 超音波診断装置
JP5844325B2 (ja) 超音波診断装置
CN101431942A (zh) 超声波诊断装置
US9402600B2 (en) 3-dimensional elastic image generation method and ultrasonic diagnostic apparatus
US20180360417A1 (en) Imaging system and method of determining a translation speed of a catheter
JP5166154B2 (ja) 超音波診断装置
US20220084496A1 (en) Method and system for ultrasonic characterization of a medium
US11761928B2 (en) Method and system for ultrasonic characterization of a medium
WO2019187649A1 (ja) 超音波診断装置および超音波診断装置の制御方法
JP4870449B2 (ja) 超音波診断装置及び超音波画像処理方法
CN115151193A (zh) 用于胎儿心脏评估的方法和系统
AU2021209270A1 (en) Method and system for ultrasonic characterization of a medium
JP4598652B2 (ja) 超音波診断装置
EP3840660A1 (en) Systems and method for performing pulse wave velocity measurements
WO2021161683A1 (ja) 超音波診断装置、超音波診断装置の制御方法および超音波診断装置用プロセッサ
CN111388012B (zh) 用于检测组织硬度的方法、设备及系统
JP4921816B2 (ja) 超音波診断装置及びその制御プログラム
EP3777697B1 (en) Ultrasonic diagnostic device and control method for ultrasonic diagnostic device
JP6041954B2 (ja) 超音波画像生成方法
JP4918344B2 (ja) 画像生成方法および超音波診断装置
JP5806189B2 (ja) 超音波診断装置および超音波画像生成方法
WO2013161228A1 (ja) 超音波診断装置および超音波診断装置の制御方法
CN116077092A (zh) 一种声速与声衰减的双工测量方法及测量设备

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

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2011553664

Country of ref document: JP

WWE Wipo information: entry into national phase

Ref document number: 13520214

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 10845707

Country of ref document: EP

Kind code of ref document: A1