WO2011007439A1 - 超音波血管検査装置 - Google Patents

超音波血管検査装置 Download PDF

Info

Publication number
WO2011007439A1
WO2011007439A1 PCT/JP2009/062909 JP2009062909W WO2011007439A1 WO 2011007439 A1 WO2011007439 A1 WO 2011007439A1 JP 2009062909 W JP2009062909 W JP 2009062909W WO 2011007439 A1 WO2011007439 A1 WO 2011007439A1
Authority
WO
WIPO (PCT)
Prior art keywords
blood vessel
image
ultrasonic
reflected wave
index value
Prior art date
Application number
PCT/JP2009/062909
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 US13/383,866 priority Critical patent/US8939908B2/en
Priority to JP2011522663A priority patent/JP5411272B2/ja
Priority to EP09847337.4A priority patent/EP2454995B1/en
Priority to CN200980160468.4A priority patent/CN102469988B/zh
Priority to PCT/JP2009/062909 priority patent/WO2011007439A1/ja
Priority to KR1020127003580A priority patent/KR101624846B1/ko
Publication of WO2011007439A1 publication Critical patent/WO2011007439A1/ja

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • A61B8/14Echo-tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0858Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/02007Evaluating blood vessel condition, e.g. elasticity, compliance
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/42Details of probe positioning or probe attachment to the patient
    • A61B8/4209Details of probe positioning or probe attachment to the patient by using holders, e.g. positioning frames
    • A61B8/4218Details of probe positioning or probe attachment to the patient by using holders, e.g. positioning frames characterised by articulated arms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461Displaying means of special interest
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/54Control of the diagnostic device
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8909Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • G01S15/8915Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
    • G01S15/8925Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array the array being a two-dimensional transducer configuration, i.e. matrix or orthogonal linear arrays
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8934Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a dynamic transducer configuration
    • G01S15/8936Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a dynamic transducer configuration using transducers mounted for mechanical movement in three dimensions
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52053Display arrangements
    • G01S7/52057Cathode ray tube displays
    • G01S7/52073Production of cursor lines, markers or indicia by electronic means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52053Display arrangements
    • G01S7/52057Cathode ray tube displays
    • G01S7/52074Composite displays, e.g. split-screen displays; Combination of multiple images or of images and alphanumeric tabular information
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52079Constructional features
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment

Definitions

  • the present invention relates to a technique for generating a blood vessel cross-sectional image using an ultrasonic probe.
  • Patent Document 1 uses an H-shaped ultrasonic probe including a first and second ultrasonic array probe that are parallel to each other and a third ultrasonic array probe that connects the central portions thereof.
  • the ultrasonic wave which measures the lumen diameter of the blood vessel, the thickness of the intima-media complex, and the like by positioning the third ultrasonic array probe on the center line of the blood vessel in a state parallel to the blood vessel A blood vessel inspection device has been proposed.
  • An apparatus of Patent Document 1 showing an example of such an ultrasonic blood vessel inspection apparatus includes a multi-axis drive device that controls the positioning state of the ultrasonic probe, and the first ultrasonic array probe is used as the first ultrasonic array probe.
  • the distance from the first ultrasonic array probe to the center of the blood vessel based on the first short-axis ultrasonic image and the second short-axis ultrasonic image by the second ultrasonic array probe.
  • the multi-axis driving device is operated to position the ultrasonic probe so that the distances from the second ultrasonic array probe to the center of the blood vessel are equal to each other.
  • the multi-axis drive device is operated to position the ultrasonic probe so that the blood vessel image is positioned at the center in the width direction of both the first and second short-axis ultrasonic images. .
  • the blood vessel immediately below the third ultrasonic array probe may be curved, and the ultrasonic probe is automatically positioned by the ultrasonic blood vessel inspection apparatus of Patent Document 1.
  • the long-axis ultrasonic image obtained by the third ultrasonic array probe does not always become a clear image capable of measuring the blood vessel diameter. Therefore, for example, after the ultrasonic probe is automatically positioned, the position of the ultrasonic probe is finely adjusted by an operator manually so that the long-axis ultrasonic image becomes a clear image capable of measuring the blood vessel diameter. I was adjusting. And the positioning (fine adjustment) operation
  • the present invention has been made against the background of the above circumstances, and the object of the present invention is to efficiently obtain a clear ultrasonic image of a blood vessel even by an operator operation with a low level of skill. It is to provide a new ultrasonic blood vessel inspection device.
  • the gist of the present invention is as follows: (a) A longitudinal section of a blood vessel in the epidermis of an organism based on the reflected wave signal of the ultrasound using an ultrasound probe disposed on the epidermis of the organism.
  • An ultrasonic angiography apparatus comprising a blood vessel longitudinal section image generating means for generating a plane image, wherein (b) the clarity of an image representing an intima complex of the blood vessel in the longitudinal section image of the blood vessel is shown.
  • the operator of the ultrasonic blood vessel inspection apparatus does not need to directly determine the clarity of the image from the longitudinal cross-sectional image (blood vessel longitudinal cross-sectional image) of the blood vessel, and the intima-complex of the blood vessel Therefore, it is possible to objectively determine the clarity of the image from the index value indicating the clarity of the image representing the image. Therefore, the position of the ultrasonic probe is corrected so that the index value is further improved, so that the skill of the operator Even if the degree is low, a clear blood vessel longitudinal cross-sectional image can be obtained efficiently.
  • the index value calculation means includes the intima composite in the front wall which is a blood vessel wall cross section close to the ultrasonic probe among a pair of blood vessel wall cross sections in the longitudinal cross section image of the blood vessel.
  • the front wall image clarity index value indicating the clarity of the image representing the body and the intima complex in the rear wall which is the blood vessel wall cross section far from the ultrasound probe of the pair of blood vessel wall cross sections
  • a rear wall image clarity index value indicating the clarity of the image is calculated. In this way, it is possible to position the ultrasonic probe so that each of the front wall and the rear wall in the blood vessel longitudinal section image becomes a clearer image.
  • Reflected wave recognition control for detecting a second peak that occurs within a range that does not exceed the value and exceeds the second peak determination threshold value that is determined in advance on each of the front wall and the rear wall.
  • Reflection to run (B) the index value calculating means includes the first peak of the plurality of reflected wave signals when the reflected wave recognition control means executes the reflected wave recognition control on the front wall.
  • the front wall image intelligibility index value is calculated based on the number of reflected wave signals in which all of the valley and the second peak are detected, and the reflected wave recognition control means performs the reflected wave recognition control on the rear wall.
  • the rear wall image clarity index value is calculated based on the number of reflected wave signals that have detected all of the first peak, the valley, and the second peak among the plurality of reflected wave signals. In this way, since the front wall image clarity index value and the rear wall image clarity index value are calculated based on a standard that does not vary, the front wall can be compared between a plurality of blood vessel examinations. The image clarity index value and the rear wall image clarity index value can be obtained.
  • the reflected wave recognition control means performs the reflected wave recognition control on the reflected wave signal received within a predetermined observation target range in the longitudinal direction of the blood vessel. In this way, compared to the case where the reflected wave recognition control is executed on all of the reflected wave signals received by the ultrasonic probe for generating the blood vessel longitudinal cross-sectional image, It is possible to reduce a control load when calculating the front wall image clarity index value and the rear wall image clarity index value.
  • the magnitude of the reflected wave signal is a value obtained by replacing the amplitude of the reflected wave signal or the amplitude of the reflected wave signal with the luminance for displaying the longitudinal cross-sectional image of the blood vessel. That's it. This facilitates detection of the first peak, the valley, and the second peak.
  • the ultrasonic probe includes a pair of first short-axis ultrasonic waves parallel to each other in which a plurality of ultrasonic oscillators are linearly arranged perpendicularly to the longitudinal direction of the blood vessel.
  • the blood vessel longitudinal section image generation The means generates a longitudinal sectional image of the blood vessel based on the reflected wave signal of the ultrasonic wave received by the long-axis ultrasonic array probe. In this way, the blood vessel longitudinal cross-sectional image can be generated using a practically used ultrasonic probe.
  • the ultrasonic angiography apparatus comprises: (a) a first short axis image display area for displaying an ultrasonic image by the first short axis ultrasonic array probe; and the second short axis use.
  • An image display device having a second short-axis image display area for displaying an ultrasonic image by an ultrasonic array probe; and a long-axis image display area for displaying a longitudinal cross-sectional image of the blood vessel; and (b) the ultrasonic wave.
  • a multi-axis drive device for controlling the positioning state of the probe, and (c) the distance from the first short-axis ultrasonic array probe to the center of the blood vessel and the second short-axis ultrasonic wave The distance from the array probe to the center of the blood vessel is equal to each other, and both the first short-axis image display area and the second short-axis image display area are in the center in the width direction.
  • the multi-axis drive device is constructed to position the blood vessel image. Short axis image position determining means for positioning the ultrasonic probe, and (d) the front wall image clarity index value and the rear wall image after completion of positioning of the ultrasonic probe by the short axis image position determining means.
  • Ultrasonic probe position correcting means for operating the multi-axis driving device to position the ultrasonic probe so that a calculated value based on the clarity index value falls within a predetermined target range. In this way, the operational burden on the operator can be reduced, and a clear blood vessel longitudinal cross-sectional image can be obtained even if the level of skill of the operator is even lower.
  • the resting diameter of the blood vessel is measured in advance before the blood vessel is released from the ischemia, the maximum diameter of the blood vessel is measured after the ischemia is released, and the diameter change after the ischemia is released relative to the resting diameter of the blood vessel.
  • a blood vessel diameter measuring means for calculating a maximum value of the ratio is provided;
  • the index value calculating means includes the front wall image clarity index value and the rear wall image clarity index at the time of measuring the rest diameter; Value and the front wall image clarity index value and the rear wall image clarity index value at the time of measuring the maximum diameter, after ischemia is released to the rest diameter of the blood vessel calculated by the blood vessel diameter measuring means An index value indicating the reliability of the maximum value of the diameter change ratio is calculated.
  • the index value calculating means continuously changes each of the front wall image clarity index value and the rear wall image clarity index value according to the magnitude of the index value.
  • the image is displayed on the image display device as a comparable image. In this way, the operator can intuitively determine the clarity of the image compared to the case where the index values are displayed numerically, and the images of the front wall and the rear wall can be clearly and quickly displayed. Can be made.
  • FIG. 2 It is a figure showing the whole structure of the ultrasonic blood vessel inspection apparatus which is one Example of this invention. It is a figure explaining the xyz-axis orthogonal coordinate axis for expressing the attitude
  • FIG. 3 is a time chart illustrating an example of a change in blood vessel lumen diameter after release of ischemia measured by ultrasound from the ultrasound probe of FIG. 2. It is the figure which showed an example of the blood vessel longitudinal cross-section image displayed on the long-axis image display area of the monitor screen display apparatus of FIG. 1, Comprising: It is a figure which shows the image with high intelligibility of the intima complex of the blood vessel wall. .
  • FIG. 7 is a view different from FIG. 6 showing an example of a blood vessel longitudinal section image displayed in the long-axis image display area of the monitor screen display device of FIG. 1, and the clarity of the intima complex of the blood vessel wall is low It is a figure which shows an image.
  • FIG. 7 is a view different from FIG. 6 showing an example of a blood vessel longitudinal section image displayed in the long-axis image display area of the monitor screen display device of FIG. 1, and the clarity of the intima complex of the blood vessel wall is low It is a figure which shows an image.
  • FIG. 2 is a functional block diagram of a first embodiment for explaining a main part of a control function provided in an electronic control device included in the ultrasonic blood vessel inspection device of FIG. 1. It is the figure which illustrated the some reflected wave signal from which the receiving position received with the long axis ultrasonic array probe of FIG. 2 differs mutually in the blood vessel longitudinal direction.
  • FIG. 10 is a diagram showing one reflected wave signal among a plurality of reflected wave signals shown in FIG. 9 in a two-dimensional coordinate system in which the signal amplitude is the vertical axis and the radial position of the blood vessel is the horizontal axis.
  • FIG. 10 is a diagram showing one reflected wave signal among a plurality of reflected wave signals shown in FIG. 9 in a two-dimensional coordinate system in which the signal amplitude is the vertical axis and the radial position of the blood vessel is the horizontal axis.
  • FIG. 11 is a diagram schematically showing a blood vessel longitudinal cross-sectional image displayed in a long-axis image display region of the monitor screen display device of FIG. 1 in order to explain the mutual relationship between FIG. 10 and the blood vessel longitudinal cross-sectional image.
  • It is the figure which showed the same reflected wave signal as FIG. 10, Comprising: It is a figure explaining the vascular tissue which the waveform of the reflected wave signal shows.
  • It is a figure which illustrates the blood vessel longitudinal cross-sectional image, blood vessel cross-sectional image, front wall image clarity index value, and back wall image clarity index value which were displayed on the monitor screen display apparatus of FIG. It is the figure which expanded the display part of the front wall image clarity index value of FIG. 13, and a rear wall image clarity index value.
  • FIG. 10 is a diagram schematically showing a blood vessel longitudinal cross-sectional image displayed in a long-axis image display region of the monitor screen display device of FIG. 1 in order to explain the mutual relationship between FIG. 10 and the blood vessel longitudinal cross-section
  • FIG. 8 illustrates a first control operation for explaining the main part of the control operation of the electronic control unit, that is, the control operation for generating a blood vessel longitudinal section image and calculating an index value indicating the clarity of the image representing the intima complex of the blood vessel.
  • It is a flowchart of an Example. 8 is a flowchart of the first embodiment for explaining the main part different from FIG. 15 of the control operation of the electronic control unit, that is, the control operation for calculating the index value indicating the measurement accuracy of the FMD evaluation.
  • FIG. 9 is a functional block diagram of a second embodiment for explaining a main part of a control function provided in an electronic control device included in the ultrasonic blood vessel inspection device of FIG. 1, corresponding to FIG. 8.
  • FIG. 8 illustrates a first control operation for explaining the main part of the control operation of the electronic control unit, that is, the control operation for generating a blood vessel longitudinal section image and calculating an index value indicating the clarity of the image representing the intima complex of the blood vessel. It is
  • FIG. 17 is a flowchart of the second embodiment for explaining the main part of the control operation of the electronic control unit, that is, the control operation for automatically positioning the ultrasonic probe at a position where FMD measurement is possible, corresponding to FIG. It is a figure to do.
  • 14 is a diagram showing the display of the front wall image clarity index value and the rear wall image clarity index value displayed on the monitor screen display device in FIG. 14, and shows another display pattern example different from FIG. 14.
  • FIG. 14 is a diagram showing the display of the front wall image clarity index value and the rear wall image clarity index value displayed on the monitor screen display device in FIG. 14, and shows another display pattern example different from FIG. 14.
  • FIG. 1 shows a non-circulation of a blood vessel 20 such as an artery located directly below the skin 18 of the upper arm 16 of the living body 14 from the upper surface 16 of the living body 14 using the probe unit 12 held by the sensor holder 10. It is a figure showing the whole structure of the ultrasonic blood vessel inspection apparatus 22 (henceforth "the blood vessel inspection apparatus 22") which performs an invasive ultrasonic diagnosis.
  • the probe unit 12 functions as a sensor for detecting biological information related to the blood vessel 20, that is, a blood vessel parameter.
  • the probe unit 12 includes a pair of first short-axis ultrasonic array probes 24a and second short-axis ultrasonic transducers 24a that are parallel to each other.
  • An H-shaped super array comprising an axial ultrasonic array probe 24b and a long-axis ultrasonic array probe 24c connecting the central portions in the longitudinal direction on a flat or flat probe surface 27.
  • the acoustic probe 24 and a multi-axis drive device (positioning device) 26 for positioning the ultrasonic probe 24 in the xyz direction and positioning the rotation angles around the x-axis and the z-axis are provided.
  • the first short-axis ultrasonic array probe 24a, the second short-axis ultrasonic array probe 24b, and the long-axis ultrasonic array probe 24c are, for example, a large number of piezoelectric ceramics. It is configured respectively in the longitudinal shape by the ultrasonic transducer (ultrasonic oscillator) a 1 ⁇ a n are linearly arranged.
  • FIG. 2 is for explaining the xyz-axis orthogonal coordinate axes used in the present embodiment, and is parallel to the longitudinal direction of the first short-axis ultrasonic array probe 24a and the first short-axis ultrasonic array.
  • the direction directly below the probe 24a and passing through or near the blood vessel 20 is defined as the x axis
  • the direction parallel to the longitudinal direction of the long axis ultrasonic array probe 24c and perpendicular to the x axis is defined as the y axis, 1
  • the direction passing through the intersection of the longitudinal direction of the short-axis ultrasonic array probe 24a and the longitudinal direction of the long-axis ultrasonic array probe 24c and perpendicular to the x-axis direction and the y-axis direction is defined as the z-axis.
  • the ultrasonic probe 24 is translated in the x-axis direction and rotated about the x-axis and the z-axis by the multi-axis drive device 26.
  • the blood vessel 20 that is, for example, the brachial artery has a three-layer structure including an intima L 1 , a media L 2 , and an adventitia L 3 . Since the reflection of ultrasonic waves occurs at different parts of the acoustic impedance, the boundary surface between the blood in the blood vessel lumen and the intima L 1 and the boundary surface between the media L 2 and the outer membrane L 3 are actually displayed in white. Organizations are displayed in black and white groups.
  • the blood vessel examination device 22 includes an electronic control device 28 including a so-called microcomputer having a CPU that processes an input signal in accordance with a program stored in the ROM in advance using the temporary storage function of the RAM, A monitor screen display device (image display device) 30, an ultrasonic drive control circuit 32, and a triaxial drive motor control circuit 34 are provided.
  • a drive signal is supplied from the ultrasonic drive control circuit 32 by the electronic control unit 28, and the first short axis ultrasonic array probe 24a and the second short axis ultrasonic array probe of the ultrasonic probe 24 of the probe unit 12 are supplied.
  • Beam-like ultrasonic waves are sequentially emitted from the transducer 24b and the long-axis ultrasonic array probe 24c by the well-known beam forming drive, and the first short-axis ultrasonic array probe 24a and the second short-axis ultrasonic array probe 24a.
  • the ultrasonic wave under the skin 18 is processed.
  • a sound wave image is generated and displayed on the monitor screen display device 30.
  • the monitor screen display device 30 includes a first short axis image display area G1 for displaying an ultrasonic image by the first short axis ultrasonic array probe 24a, and a first short axis image display area G1.
  • a second short-axis image display region G2 for displaying an ultrasonic image by the two-short-axis ultrasonic array probe 24b and an ultrasonic image (blood vessel longitudinal cross-sectional image) by the long-axis ultrasonic array probe 24c are displayed.
  • first short-axis image display area G1, the second short-axis image display area G2, and the long-axis image display area G3 are provided with a common vertical axis indicating the depth dimension from the skin 18. is there.
  • the monitor screen display device 30 displays the rate of change of the intima diameter, that is, the rate of expansion R of the lumen diameter in time series when evaluating the FMD (blood flow dependent vasodilator response).
  • the ultrasonic probe 24 is controlled by the electronic control unit 28 so as to be at a predetermined measurement position PT1 with respect to the blood vessel 20.
  • the multi-axis drive device 26 to which the drive signal is supplied is driven to be positioned.
  • the predetermined measurement position PT1 is a position where the first short-axis ultrasonic array probe 24a and the second short-axis ultrasonic array probe 24b are orthogonal to the blood vessel 20, and for the long axis. This is a position where the ultrasonic array probe 24 c is parallel to the blood vessel 20.
  • the distance from the first short-axis ultrasonic array probe 24a to the center of the blood vessel 20 is equal to the distance from the second short-axis ultrasonic array probe 24b to the center of the blood vessel 20, and the first Both the first short-axis image display region G1 and the second short-axis image display region G2 are measurement positions at which the image of the blood vessel 20 is positioned at the center in the width direction.
  • the sensor holder 10 is in a state in which the blood vessel 20 positioned immediately below the skin 18 of the upper arm 16 of the living body 14 is lightly brought into contact with the desired position in the three-dimensional space, that is, the predetermined measurement position PT1 so as not to deform.
  • the probe unit 12 is held in a desired posture.
  • a coupling agent such as jelly, olive oil or glycerin, or a water bag containing water in a resin bag is interposed.
  • the sensor holder 10 includes, for example, a magnet base 36 fixed to a desk, a pedestal, etc. by magnetic attraction, a unit fixture 38 to which the probe unit 12 is fixed, and one end to the magnet base 36 and the unit fixture 38. Are connected to and supported by the magnet base 36 and the unit fixture 38 via the connecting members 44 and 45, and the connecting members 44 and 45 having the tip 42 formed in a spherical shape.
  • a universal arm 40 is provided.
  • the universal arm 40 includes two links 46 and 47 that are rotatably connected to each other, and a distal end portion 42 of which one end of the links 46 and 47 is energized with a predetermined resistance. And the other ends of the links 46 and 47 are connected to each other so that they can rotate relative to each other.
  • a rotating joint portion 54 that is made relatively unrotatable by a fastening force obtained by tightening a male threaded fixing knob 52 that is screwed into a screw hole penetrating the connecting portion.
  • the multi-axis drive device 26 includes an x-axis rotation (yawing) mechanism fixed to the unit fixture 38 in order to position the rotation position around the x-axis of the ultrasonic probe 24 by an x-axis rotation actuator, An x-axis translation mechanism for positioning the translation position of the ultrasonic probe 24 in the x-axis direction by the rotation actuator, and a z-axis rotation for positioning the rotation position of the ultrasonic probe 24 around the z-axis by the z-axis actuator. And a moving mechanism.
  • the multi-axis drive device 26 controls the positioning state of the ultrasonic probe 24 with such a configuration.
  • the ultrasonic drive control circuit 32 is arranged in accordance with a command from the electronic control unit 28, for example, a large number of ultrasonic transducers arranged in a line constituting the first short axis ultrasonic array probe 24a.
  • the ultrasonic transducer a 1 et of the end, a certain number of the ultrasonic transducer groups for example 15 a 1 to a 15 a frequency of about 10MHz while applying a predetermined phase difference for each
  • the beam forming drive that is driven at the same time causes the convergent ultrasonic beam to be sequentially emitted toward the blood vessel 20 in the arrangement direction of the ultrasonic transducers, and the ultrasonic beam is scanned while shifting the ultrasonic transducers one by one.
  • a reflected wave for each radiation at the time of (scanning) is received and input to the electronic control unit 28.
  • the electronic control unit 28 synthesizes an image based on the reflected wave, generates a cross-sectional image (short axis image) or a vertical cross-sectional image (long axis image) of the blood vessel 20 under the skin 18, and displays the monitor screen.
  • the information is displayed on the device (image display device) 30. Further, from the image, the diameter of the blood vessel 20 or the diameter of the endothelium 70 (endothelium diameter (lumen diameter) d 1 ) is calculated.
  • FIG. 5 is a ischemia (avascularization) a time chart illustrating the change in the vessel lumen diameter d 1 after opening.
  • t1 point represents the time ischemic opened, it begins to expand vascular lumen diameter d 1 from time t2, vascular lumen diameter d 1 at the time t3 has reached its maximum value d MAX It is shown. Therefore, the expansion rate R of the blood vessel lumen diameter calculated by the electronic control unit 28 is maximized at time t3.
  • the ischemia for FMD evaluation is controlled by a cuff pressure control unit 56 (cuff pressure control means 56) provided in the electronic control device 28 using a pressure control valve 60 to control the original pressure from the air pump 58. Then, the pressure is supplied to the cuff 62 wound around the upper arm 16 and the pressure (cuff pressure) of the cuff 62 is increased to a predetermined ischemic cuff pressure exceeding the maximum blood pressure of the living body 14. At this time, the cuff pressure control unit 56 detects the cuff pressure by a signal from the pressure sensor 64 for detecting the cuff pressure.
  • a cuff pressure control unit 56 cuff pressure control means 56
  • the cuff pressure control unit 56 maintains the cuff pressure at the ischemic cuff pressure for a predetermined time before the release of the ischemia, that is, a predetermined time before the time t1, and at the time of releasing the ischemia (time t1). Immediately reduce the cuff pressure to atmospheric pressure.
  • the ultrasonic probe 24 is positioned at the predetermined measurement position PT1, so that it is basically displayed in the long-axis image display region G3 in FIG.
  • An image having a high clarity of the intima-media complex IMC (intima-media complex) that combines the intima L 1 and the intima L 2 of the blood vessel 20 is obtained like the longitudinal cross-sectional image of the blood vessel.
  • IMC intima-media complex
  • BL IMC of the intima-media complex IMC is clearly displayed with continuity in the longitudinal direction of the blood vessel 20 in any of the rear wall BR B which is FIG. 6 can be said to be an image with high clarity of the inner-media complex IMC. If blood vessel image diagnosis such as FMD evaluation is performed based on such a clear image, it can be said that the image diagnosis is sufficiently accurate.
  • the intima-media complex IMC as shown in FIG. 6 has high clarity due to the blood vessel 20 being curved or the like. In some cases, a blood vessel longitudinal section image cannot be obtained. For example, it is a case where the intima-media complex IMC as shown in FIG. In FIG. 7, unlike FIG. 6, the boundary between the blood vessel lumen and the blood vessel wall is discontinuous in the longitudinal direction of the blood vessel 20 in both the front wall BR F and the rear wall BR B. hypoechoic band BL IMC body IMC is not the most clearly can not be recognized from the image. Ambiguous inner in membrane complex IMC as in FIG.
  • the blood vessel inspection device 22 of this embodiment has a control function that supports the fine adjustment, and the main part of the control function will be described with reference to FIG.
  • FIG. 8 is a functional block diagram for explaining the main part of the control function provided in the blood vessel inspection device 22.
  • a blood vessel image evaluation unit 100 (see FIG. 1) provided in the electronic control unit 28 includes an image generation execution determination unit 80 as an image generation execution determination unit and a blood vessel as a blood vessel cross-sectional image generation unit.
  • a cross-sectional image generation unit 82, an index value calculation unit 84 as an index value calculation unit, and a blood vessel diameter measurement unit 88 as a blood vessel diameter measurement unit are provided.
  • the image generation execution determining means 80 determines whether or not to acquire an ultrasonic image of the blood vessel 20 using the ultrasonic probe 24. That is, the acquisition start and acquisition end of the ultrasound image of the blood vessel 20 are determined.
  • the image generation execution determination unit 80 determines that the ultrasonic wave is turned on when the switch is turned ON. It is determined that acquisition of an image is started, and when the switch is turned off, it is determined that acquisition of the ultrasonic image is completed.
  • the blood vessel cross-sectional image generating means 82 functions as a blood vessel longitudinal cross-sectional image generating means, and uses the ultrasonic probe 24 arranged on the epidermis of the living body 14 and the subepidermal surface of the living body 14 based on the reflected wave signal SG EC of the ultrasonic wave.
  • a longitudinal section image (blood vessel longitudinal section image) of the blood vessel 20 is sequentially generated. Specifically, sequentially generates the blood vessel longitudinal sectional image based on the ultrasonic reflected wave signal SG EC which are successively and repeatedly received by the ultrasonic detector array 24c for the long axis.
  • the longitudinal cross-sectional image of the blood vessel 20 is an image displayed in the long-axis image display region G3 as shown in FIG. 6 or FIG. 7, and the blood vessel cross-sectional image generation unit 82 uses a general B mode method or the like. Generate using.
  • the blood vessel cross-sectional image generating means 82 when one of the blood vessel longitudinal cross-sectional images is generated, the blood vessel cross-sectional image generating means 82 in detail, the ultrasonic reflected wave signal SG EC received by the long-axis ultrasonic array probe 24c.
  • the blood vessel longitudinal section image is generated based on the above. That is, the blood vessel cross-sectional image generation unit 82 scans the reflected wave signal SG EC ultrasound in the longitudinal direction at a predetermined reflection wave reception interval (line pitch) PC RV in ultrasonic array probe 24c for the long axis It received, to generate the blood vessel longitudinal sectional image based on the plurality of reflected wave signal SG EC which is the reception.
  • the blood vessel cross-sectional image generating means 82 also functions as a blood vessel cross-sectional image generating means, and the ultrasonic wave reflected wave signal SG EC sequentially and repeatedly received by the first short-axis ultrasonic array probe 24a.
  • the cross-sectional images (blood vessel cross-sectional images) of the blood vessels 20 to be displayed in the first short-axis image display region G1 are sequentially generated, and sequentially received by the second short-axis ultrasonic array probe 24b. sequentially generating a cross sectional image of the blood vessel 20 to be displayed on the second short axis image display area G2 based on the reflected wave signal SG EC of sound waves.
  • the blood vessel cross-sectional image generation unit 82 flickers in the image displayed on the monitor screen display device 30 until the image generation execution determination unit 80 determines that the acquisition of the ultrasonic image is started and then ends.
  • the generation of the blood vessel longitudinal cross-sectional image and the blood vessel cross-sectional image is repeatedly performed at a predetermined period T A that does not occur.
  • the index value calculation unit 84 calculates an index value indicating the clarity of an image representing the intima-media complex IMC of the blood vessel 20 in the blood vessel longitudinal cross-sectional image generated by the blood vessel cross-sectional image generation unit 82. For example, every time the blood vessel cross-sectional image generating unit 82 generates a blood vessel longitudinal cross-sectional image, the index value is sequentially calculated. Specifically, the index value is calculated for each of the front wall BR F and the rear wall BR B in the blood vessel longitudinal sectional image.
  • the index value calculating means 84 a front wall image clarity index value XCR F indicating the clarity of the image representing the intima-media complex IMC in the front wall BR F in the blood vessel longitudinal sectional image, the blood vessel longitudinal section A rear wall image clarity index value XCR B indicating the clarity of the image representing the intima-media complex IMC in the rear wall BR B in the image is calculated.
  • the index value calculation means 84 includes a reflected wave recognition control means 86 as a reflected wave recognition control section in order to calculate the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B. Yes.
  • the clarity of the image representing the intima-media complex IMC means that the intima-media complex IMC exists along the longitudinal direction of the blood vessel 20, so It can be said that the low-echo band BL IMC and the high-echo band image representing the membrane complex IMC can be recognized continuously in the longitudinal direction of the blood vessel 20.
  • the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B are collectively expressed, they represent an image representing the intima complex IMC. It can be explained that it is an index value indicating the degree of clarity.
  • the reflected wave recognition control unit 86 receives the plurality of reflected wave signals SG EC received by the blood vessel cross-sectional image generation unit 82 to generate the blood vessel vertical cross-sectional image, for example, the blood vessel cross-sectional image generation unit 82 generates the blood vessel vertical cross-sectional image. Get every time it is generated. For example, it is acquired from the blood vessel cross-sectional image generation means 82. Looking at the above plurality of the reflected wave signal SG EC for generating one of a blood vessel longitudinal sectional image, each of the plurality of reflected wave signal SG EC has a length of axis ultrasonic array probe 24c vessels longitudinal The signals are received at different reception positions. That is, the reflected wave recognition control means 86, as illustrated in FIG. 9, to obtain a plurality of reflected wave signal SG EC reception positions mutually different in the vessel longitudinal length axis ultrasonic array probe 24c .
  • Reflected wave recognition control means 86 obtained for each of said plurality of reflected wave signal SG EC, the relationship between the radial position PT R of the reflected wave signal SG EC signal amplitude AM SG and the blood vessel 20 as shown in FIG. 10 , The first peak PK1 in which the signal amplitude AM SG of the reflected wave signal SG EC exceeds a predetermined first peak determination threshold value LT1 PK , and radially outside the blood vessel 20 from the position where the first peak PK1 occurs.
  • a small valley BTM than valleys determination threshold LT BTM signal amplitude AM SG of the generated reflected wave signal SG EC is predetermined, the radially outer vessel 20 through the valley BTM from the generation position of the first peak PK1 Reflection that occurs within a range that does not exceed a predetermined peak interval threshold LT SP and that detects a second peak PK2 in which the signal amplitude AM SG of the reflected wave signal SG EC exceeds a predetermined second peak determination threshold LT2 PK wave recognition control is performed on each of the front wall BR F and the rear wall BR B This reflected wave recognition control will be specifically described with reference to FIG.
  • Figure 10 is a reflected wave of the signal amplitude AM SG of the reflected wave signal SG EC radial position PT R of the (signal strength AM SG) was used as a vertical axis vessel 20 to the coordinate system of the horizontal axis, one (one) It is a figure showing the part applicable to front wall BR F or rear wall BR B of signal SG EC .
  • a portion of the reflected wave signal SG EC represented in FIG. 10 corresponds to the front wall BR F
  • arrow AR01 of FIG. 10 while corresponding to the arrow AR02 in vascular longitudinal sectional image of FIG. 11
  • the first peak PK1, valleys BTM, and the waveform of the reflected wave signal SG EC in each of the second peak PK2 is pointed at an acute angle which is one example, for example, the first peak PK1, In any of the valley BTM and the second peak PK2, the waveform may be substantially flat like a trapezoid without being sharpened.
  • Reflected wave recognition control unit 86 first, recognizes a portion corresponding to the blood vessel lumen in the reflected wave signal SG EC. For example, it may be recognized from cross-sectional images of the blood vessel 20 (short axis image), an intermediate position of that portion of the reflected wave signal SG EC indicating each of the front wall BR F and the rear wall BR B in the vessel lumen You may recognize that there is. Further, in the present embodiment, as described later, the reflection wave recognition control with respect to the longitudinal direction of the observation target range AOB received in (see FIG.
  • the reflected wave signal SG EC of the blood vessel 20 defined in advance Since the portion corresponding to the lumen of the blood vessel is recognized, the reflection received in the observation target range AOB among the plurality of reflected wave signals SG EC to be a basis for generating the blood vessel longitudinal section image is recognized. Only the wave signal SG EC is sufficient.
  • the reflective wave recognition control means 86 before either in one wall BR F and the rear wall BR B, toward the vessel outwardly from the vessel cavity side as indicated by an arrow AR01 in FIG. 10
  • the signal amplitude AM SG of the reflected wave signal SG EC is recognized, and if there is a maximum point where the signal amplitude AM SG exceeds the first peak determination threshold LT1 PK , it is detected as the first peak PK1.
  • the valley determination threshold value LT BTM outside the radial direction of the blood vessel 20 from the position where the first peak PK1 is generated (the horizontal axis position in FIG. 10).
  • the signal amplitude AM SG satisfies the second peak determination threshold value LT2 PK within a range that does not exceed the peak interval threshold value LT SP radially outside the blood vessel 20 through the valley BTM from the position where the first peak PK1 occurs. If there is a local maximum exceeding it, it is detected as the second peak PK2.
  • Reflected wave recognition control unit 86 executes such a reflected wave recognition control at each of the front wall BR F and the rear wall BR B for each of the plurality of reflected wave signal SG EC.
  • the comparison target is not limited to the signal amplitude AM SG , and may be any size as long as the reflected wave signal SG EC associated with the blood vessel longitudinal cross-sectional image.
  • the signal amplitude AM SG is luminance-modulated. The brightness of the blood vessel when it is replaced with the brightness for displaying the blood vessel longitudinal section image (B-mode long axis image) may be used.
  • first peak determination threshold LT1 PK , the valley determination threshold LT BTM , and the second peak determination threshold LT2 PK are, for example, super-corresponding to the first peak PK1, the valley BTM, and the second peak PK2, respectively.
  • the sound wave images are experimentally determined so that they can be visually recognized as different patterns, and the valley determination threshold value LT BTM is smaller than the first peak determination threshold value LT1 PK and the second peak determination threshold value LT2 PK. Is set. In FIG.
  • the second peak PK2 is shown to be larger than the first peak PK1, but conversely the second peak PK2 may be smaller than the first peak PK1,
  • the first peak determination threshold LT1 PK and the second peak determination threshold LT2 PK are not particularly limited in their magnitude relationship, and the determination thresholds LT1 PK and LT2 PK may be different from each other or the same. It may be a value.
  • the low echo band BL of the intima-media complex IMC in the blood vessel longitudinal cross-sectional image is shown to be larger than the first peak PK1, but conversely the second peak PK2 may be smaller than the first peak PK1,
  • the first peak determination threshold LT1 PK and the second peak determination threshold LT2 PK are not particularly limited in their magnitude relationship, and the determination thresholds LT1 PK and LT2 PK may be different from each other or the same. It may be a value.
  • the peak interval threshold value LT SP is experimentally set based on the thickness of the intima-media complex IMC of the living body, for example.
  • the reflected wave recognition control unit 86 may execute the reflected wave recognition control for all of the plurality of reflected wave signal SG EC underlying the vessel longitudinal section images, calculation load in this embodiment for reduction, as shown in FIG. 11, for the received reflected wave signal SG EC in the longitudinal direction of the observation target range of the vessel 20 a predetermined AOB, executes the reflected wave recognition control.
  • This observation target range AOB is, for example, a part of the longitudinal width of the long-axis ultrasonic array probe 24c and may be fixed or changeable by an operator. It is experimentally set so that the clarity of the intima complex IMC in the image can be judged.
  • first peak PK1 and valleys BTM and second peak PK2 indicates an example of a reflected wave signal SG EC detected
  • first peak PK1 and valleys also exists a reflected wave signal SG EC either or both are not detected in the BTM and the second peak PK2.
  • the reflected wave recognition control means 86 the observation target range AOB of the plurality of reflected wave signal SG EC received in, the execution of the reflected wave recognition control the first peak PK1 and valleys BTM and second peak PK2 number QL SG of the reflected wave signal SG EC detecting a whole, i.e., the peak recognition line number QL SG, stores for each of the front wall BR F and the rear wall BR B.
  • observation target range AOB / reflected wave reception interval PC RV Suppose there is.
  • the front wall image clarity index value XCR F is calculated based on the number QL SG (peak recognition line number QL SG ).
  • the index value calculation means 84 is configured such that the reflected wave recognition control means 86 performs the reflected wave recognition control on the rear wall BR B , so that the first peak PK1, the valley BTM, and the peak among the plurality of reflected wave signals SG EC
  • the rear wall image intelligibility index value XCR B is calculated based on the number of reflected wave signals QL SG (peak recognition line number QL SG ) that has detected all of the second peak PK2.
  • the index value calculation means 84 calculates the front wall image clarity index value XCR F, that is, the front wall score XCR F by the following formula (1), and the rear wall image clarity by the following formula (2).
  • the index value XCR B that is, the rear wall score XCR B is calculated.
  • XCR F indicates the front wall score XCR F
  • XCR B indicates the rear wall score XCR B
  • QLF SG indicates the front wall.
  • QLB SG represents a peak recognition line number QL SG in the rear wall BR B
  • PC RV indicates the reflected wave reception interval (line pitch) PC RV
  • AOB indicates the observation target range (observation area width) AOB.
  • XCR F (QLF SG ⁇ PC RV ) / AOB ⁇ 100 (1)
  • XCR B (QLB SG x PC RV ) / AOB x 100 (2)
  • the index value calculation means 84 sequentially calculates the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B.
  • the image clarity index values XCR F and XCR B are relative values whose change ranges are “0 to 100”.
  • the index value calculation means 84 surrounds each of the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B by a two-dot chain line L02 in FIG.
  • the monitor screen display device 30 sequentially displays the blood vessel longitudinal section image and the blood vessel transverse section image. Specifically, as shown in FIG. 14 which is an enlarged view of a portion surrounded by a two-dot chain line L02 in FIG.
  • the index value calculation means 84 includes the front wall image clarity index value XCR F and the rear wall image clarity.
  • Each of the degree index values XCR B is displayed as a numerical value and is continuously displayed on the monitor screen display device 30 as an image or a figure that changes continuously according to the size of the index values XCR F and XCR B and can be compared with each other. Display.
  • each of the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B has a center angle and an area that increase as the size thereof increases, and exhibits line symmetry with respect to a common straight line. It is expressed by a pair of fan-shaped fan graphs (black portions in FIG. 14). When the maximum value is “100”, the fan-shaped graph has a maximum area and becomes a semicircle. When the index values XCR F and XCR B are both maximum, the semicircle forms a complete circle.
  • Vascular diameter measuring unit 88 non-invasively from the blood vessel longitudinal sectional image to measure the blood vessel lumen diameter d 1.
  • the blood vessel diameter measuring means 88 in FMD evaluation, in order to measure the diameter change rate of the blood vessel 20 after ischemia opening of the vessel 20 (expansion ratio of the vessel lumen diameter d 1 R), the vessel 20 previously measuring the vessel lumen diameter d a (rest diameter d a) at rest in prior ischemia open.
  • the blood vessel diameter measuring means 88 for FMD evaluation, the vessel lumen diameter d 1 sequentially measured ischemia open predetermined blood vessel diameter measuring the periodic inter TIME1 after further blood vessel lumen diameter was the measured
  • the diameter change rate R of the blood vessel 20 is sequentially calculated from d 1 and the rest diameter d a .
  • the blood vessel diameter measuring means 88 sequentially and continuously measures the blood vessel lumen diameter d 1 that changes thereafter from when the ischemia is released over the blood vessel diameter measuring period TIME1.
  • a measurement time of the blood vessel lumen diameter d 1 is expected to be substantially maximized one or more experimentally based on the time of ischemic opened, by measuring the blood vessel lumen diameter d 1 at the measurement point Also good.
  • the blood vessel diameter measuring time period TIME1 is intravascular stored in vessel diameter measuring means 88 is experimentally set on the basis of the time ischemic opened to measure the maximum value d MAX vascular lumen diameter d 1 after ischemia open lumen diameter is the measurement period d 1, starting from the vessel lumen diameter d 1 is and when ischemia opening includes a time point (t3 time) to reach its maximum value d MAX (t1 time) as shown in FIG. Therefore, the blood vessel diameter measuring means 88 measures the maximum diameter (maximum lumen diameter) d MAX of the blood vessel 20 after the ischemia is released.
  • the vascular diameter measuring unit 88 after a lapse of the blood vessel diameter measuring time period TIME1, the maximum diameter change ratio after ischemia opening for said resting diameter d a of the blood vessel 20 R MAX, i.e., relative to the its rest diameter d a
  • Index value calculating means 84 in addition to the above functions, and the front wall image clarity index XCR F and the rear wall image clarity index value during measurement of the resting diameter d a XCR B, the maximum size after ischemia open The front wall image clarity index value XCR F and the rear wall image clarity index value XCR B at the time of measuring d MAX are calculated.
  • the index value calculating means 84 calculates an average value of the front wall image clarity index XCR F and the rear wall image clarity index value XCR B during measurement of the resting diameter d a as resting diameter measuring scheduled score XCR1
  • the average value of the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B at the time of measuring the maximum diameter d MAX is calculated as the maximum diameter measurement score XCR2.
  • the FMD reliability index value XCR FMD is calculated from the rest diameter measurement score XCR1 and the maximum diameter measurement score XCR2 by the following equation (3).
  • the index value calculation means 84 displays the calculated FMD reliability index value XCR FMD on the monitor screen display device 30 together with the vascular dilation rate evaluation value RMAX .
  • the FMD reliability index value XCR FMD is an index value indicating the reliability of the vascular dilation rate evaluation value R MAX as described above, in other words, an index indicating the measurement accuracy of the vascular dilation rate evaluation value R MAX.
  • a value, FMD reliability index value XCR FMD is an indication that a high measurement accuracy of the vasodilating factor evaluation value R MAX larger the value.
  • the display control unit 90 (display control unit 90) provided in the electronic control unit 28 is configured to generate the blood vessel longitudinal cross-sectional image and the blood vessel cross-sectional image generated by the blood vessel cross-sectional image generating unit 82.
  • the monitor screen display device 30 is sequentially displayed on the monitor screen display device 30, and the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B calculated by the index value calculation means 84 are displayed on the monitor screen as numerical values and figures (fan charts). The information is sequentially displayed on the device 30.
  • the display control unit 90 determines the vascular dilation rate evaluation value R MAX and the FMD reliability.
  • the degree index value XCR FMD is displayed on the monitor screen display device 30.
  • FIG. 15 shows an index value indicating the clarity of an essential part of the control operation of the blood vessel inspection device 22 (electronic control device 28), that is, an image representing the intima complex IMC of the blood vessel 20 generated by the blood vessel longitudinal section image. It is a flowchart for demonstrating the control action
  • step it is determined whether or not acquisition of an ultrasound image of the blood vessel 20 is started. If the determination of SA1 is affirmative, that is, if acquisition of the ultrasonic image is started, the process proceeds to SA2. On the other hand, if the determination of SA1 is negative, the flowchart of FIG. 15 ends.
  • the reflected wave recognition control is performed for the front wall BR F for each of the plurality of reflected wave signals SG EC received in the observation target range AOB.
  • the front wall image clarity index value (front wall score) XCR F is the valley between the first peak PK1 of the front wall BR said plurality of reflected wave signal by the execution of the reflected wave recognition control in F SG EC It is calculated based on the number QL SG (QLF SG ) of the reflected wave signals SG EC that have detected all of the BTM and the second peak PK2.
  • the reflected wave recognition control the rear wall BR B for each of the plurality of reflected wave signal SG EC received by said observation target range AOB
  • the rear wall image clarity index value (rear wall score) XCR B is obtained by executing the reflected wave recognition control on the rear wall BR B , among the plurality of reflected wave signals SG EC and the first peak PK1 and the valley. It is calculated based on the number QL SG (QLB SG ) of reflected wave signals SG EC that have detected all of the BTM and the second peak PK2.
  • the calculated rear wall score XCR B is displayed on the monitor screen display device 30 in real time.
  • the front wall score XCR F and the rear wall score XCR B are each displayed as a numerical value that changes between 0 and 100, and the front wall score XCR F and the rear wall score XCR If both B and the maximum value are “100”, they are displayed as a pair of sector graphs that form a complete circle.
  • SA6 corresponding to the image generation execution determining means 80, it is determined whether or not to end the acquisition of the ultrasound image of the blood vessel 20.
  • the determination of SA6 is affirmed, that is, when the acquisition of the ultrasonic image is finished, the flowchart of FIG. 15 is finished.
  • the determination of SA6 is negative, that is, when the acquisition of the ultrasonic image is continued, the process returns to SA2. Accordingly, the steps from SA2 to SA5 are repeatedly executed from the start of acquisition of the ultrasonic image to the end thereof. For example, it is repeatedly executed in a very short cycle.
  • FIG. 16 is a flowchart for explaining the main part of the control operation of the blood vessel inspection device 22 (electronic control device 28) different from FIG. 15, that is, the control operation for calculating the index value indicating the measurement accuracy of the FMD evaluation. is there.
  • the control operation shown in FIG. 16 is executed alone or in parallel with other control operations.
  • the reflected wave recognition control the blood vessel longitudinal section for measuring the resting diameter d a of the front the ischemia opening in FMD evaluation for each of the plurality of reflected wave signal SG EC underlying the image is performed in each of the front wall BR F and the rear wall BR B.
  • the reflected wave recognition control is executed for each of the plurality of reflected wave signal SG EC received by said observation target range AOB.
  • the reflected wave recognition control of the blood vessel longitudinal sectional image for determining the maximum diameter d MAX after the ischemia opening in FMD evaluation for each of the plurality of reflected wave signal SG EC underlying, it is executed in each of the front wall BR F and the rear wall BR B. Similar to the above in detail SB1, rather than to all of the plurality of reflected wave signal SG EC, for each of the plurality of reflected wave signal SG EC received by said observation target range AOB, the reflected wave recognition Control is executed.
  • the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B at the time of measuring the maximum diameter d MAX are calculated, and the average value thereof is the maximum value. Calculated as the diameter measurement score XCR2.
  • the FMD reliability index value XCR FMD is calculated from the above-mentioned rest diameter measurement score XCR1 and maximum diameter measurement score XCR2 by the above equation (3).
  • the FMD reliability index value XCR FMD is, FMD evaluation result for example the vasodilation factor evaluation value (% FMD value) together with R MAX, is displayed on a monitor screen display apparatus 30.
  • the blood vessel cross-sectional image generating means (blood vessel longitudinal cross-sectional image generating means) 82 uses the ultrasonic probe 24 arranged on the epidermis of the living body 14 to reflect the reflected wave signal SG EC of the ultrasonic wave. Based on the above, the longitudinal cross-sectional images of the subepidermal blood vessel 20 of the living body 14 are sequentially generated. Then, the index value calculation unit 84 calculates an index value indicating the clarity of the image representing the intima complex IMC of the blood vessel 20 in the longitudinal cross-sectional image of the blood vessel 20 generated by the blood vessel cross-sectional image generation unit 82. .
  • the operator of the blood vessel inspection device 22 does not need to directly determine the clarity of the image from the longitudinal cross-sectional image (blood vessel longitudinal cross-sectional image) of the blood vessel 20, and represents the intima-media complex IMC of the blood vessel 20. Since the intelligibility of the image can be objectively determined from the index value indicating the intelligibility of the image, the position of the ultrasonic probe 24 can be easily fine-tuned so that the index value is further improved. It is possible to efficiently obtain a clear blood vessel longitudinal section image even if the degree of skill is low.
  • the index value calculating means 84 the blood vessel longitudinal sectional wall image clarity before indicating the clarity of the image representing the intima-media complex IMC in the front wall BR F in the image Since the index value XCR F and the rear wall image clarity index value XCR B indicating the clarity of the image representing the intima-media complex IMC in the rear wall BR B in the blood vessel longitudinal cross-sectional image are calculated, for example, an operator However, based on the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B , each of the front wall BR F and the rear wall BR B in the blood vessel longitudinal section image becomes a clearer image. Thus, the ultrasonic probe 24 can be positioned efficiently.
  • the reflected wave recognition control means 86 has a plurality of reflected wave signals SG received at different reception positions in the blood vessel longitudinal direction of the long-axis ultrasonic array probe 24c.
  • the signal amplitude AM SG of the reflected wave signal SG EC is predetermined
  • the first peak PK1 exceeding the first peak determination threshold LT1 PK , and the signal amplitude AM SG of the reflected wave signal SG EC generated in the radial direction outside the blood vessel 20 from the generation position of the first peak PK1 are predetermined.
  • a valley BTM smaller than the valley determination threshold LT BTM and the first peak PK1 are generated within a range that does not exceed a predetermined peak interval threshold LT SP radially outside the blood vessel 20 through the valley BTM from the generation position of the first peak PK1.
  • second peak determination threshold LT2 PK signal amplitude AM SG of the reflected wave signal SG EC is predetermined The reflected wave recognition control for detecting a second peak PK2 exceeding, run on each of the front wall BR F and the rear wall BR B.
  • the front wall image clarity index value XCR F is calculated based on the number QL SG of reflected wave signals SG EC that have detected all of the second peak PK2, and the reflected wave recognition control means 86 performs the reflection on the rear wall BR B.
  • the rear wall An image clarity index value XCR B is calculated. Accordingly, since the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B are calculated based on a standard that does not vary, the front wall image can be compared between a plurality of blood vessel examinations. The clarity index value XCR F and the rear wall image clarity index value XCR B can be obtained.
  • the reflected wave recognition control means 86 the reflection wave signal SG EC received in the longitudinal direction in the observation scope AOB of the blood vessel 20 to a predetermined, the reflection Since the wave recognition control is executed, the reflected wave recognition control is performed on all of the plurality of reflected wave signals SG EC received by the long-axis ultrasonic array probe 24c for generating the blood vessel longitudinal section image. Compared with the case where it is executed, it is possible to reduce the control load when calculating the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B.
  • the reflected wave signal SG EC signal amplitude AM SG is the first peak determination threshold LT1 PK, the valleys determination threshold LT BTM, and the second since it is compared to the peak decision threshold LT2 PK, the in the reflected wave signal SG EC first peak PK1, the valleys BTM, and it is easy detection of the second peak PK2.
  • the luminance height after luminance modulation of the signal amplitude AM SG is the first peak determination threshold LT1 PK , the valley determination threshold LT BTM , and the above. The same applies when compared with the second peak determination threshold LT2 PK .
  • the ultrasonic probe 24 includes a plurality of ultrasonic oscillators arranged in a straight line perpendicular to the longitudinal direction of the blood vessel 20.
  • a long-axis ultrasonic array probe 24c provided adjacent to the center of the probe 24b and in which a plurality of ultrasonic oscillators are linearly arranged in the longitudinal direction of the blood vessel 20 is provided on a single plane.
  • the blood vessel longitudinal sectional image generation unit 82 generates the blood vessel longitudinal sectional image based on the reflected wave signal SG EC of the received ultrasonic its long axis ultrasonic array probe 24c. Therefore, the blood vessel longitudinal section image can be generated using an ultrasonic probe that has been put into practical use.
  • the blood vessel diameter measuring means 88 before ischemia opening of the blood vessel 20 measured in advance rest diameter d a of the blood vessel 20, the maximum diameter d of the blood vessel 20 after ischemia opening of the vessel 20 measuring the MAX, and calculates the maximum value R MAX of diameter change rate after ischemia open for resting diameter d a of the vessel 20.
  • the index value calculating means 84 a front wall image clarity index XCR F and the rear wall image clarity index value XCR B during measurement of the resting diameter d a, before during measurement of the maximum diameter d MAX wall based on the image clarity index XCR F and the rear wall image clarity index value XCR B, the maximum diameter change ratio R MAX after ischemia open for resting diameter d a of the blood vessel 20 that has been calculated by the blood vessel diameter measuring means 88
  • An index value indicating reliability (FMD reliability index value) XCR FMD is calculated.
  • the index value calculation means 84 calculates the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B as numerical values.
  • the image is displayed on the monitor screen display device 30 as an image or a figure (fan graph) that continuously changes according to the index values XCR F and XCR B and can be compared with each other. Therefore, the operator can intuitively determine the clarity of the image as compared with the case where the index values XCR F and XCR B are displayed only by numerical values, and the images of the front wall BR F and the rear wall BR B can be determined. It becomes possible to clarify more quickly.
  • the operator of the blood vessel inspection device 22 makes fine adjustment of the position of the ultrasonic probe 24 arranged on the skin of the living body 14 by manual operation, and the front wall image clarity index value XCR F and The rear wall image clarity index value XCR B is calculated and displayed, and the multi-axis drive device 26 is controlled based on the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B. It is possible to automatically fine-tune the position of the ultrasonic probe 24 so that the clarity of the intima-media complex IMC in the blood vessel longitudinal cross-sectional image is greater than a certain limit by driving. In the second embodiment, control for automatically fine-tuning the position of the ultrasonic probe 24 will be described.
  • the control of the second embodiment may be executed in place of the display control of the first embodiment described above, but may be executed simultaneously with the display control of the first embodiment.
  • description of points common to the first embodiment will be omitted, and differences from the first embodiment will be mainly described.
  • FIG. 17 is a functional block diagram illustrating a main part of the control function provided in the blood vessel inspection device 22 (blood vessel image evaluation unit 100), and corresponds to FIG. 8 of the first embodiment.
  • the blood vessel image evaluation unit 100 of the present embodiment includes a blood vessel cross-sectional image generation means 82, an index value calculation means 84, and a blood vessel diameter measurement means 88, and further serves as a short-axis image position determination section.
  • a short-axis image position determining unit 210 and an ultrasonic probe position correcting unit 212 as an ultrasonic probe position correcting unit are provided.
  • the ultrasonic probe position correcting unit 212 includes a long-axis image score calculating unit 214 as a long-axis image score calculating unit.
  • Blood vessel cross-sectional image generation unit 82 similarly to the first embodiment, the first short axis image display area based on the ultrasonic reflected wave signal SG EC received in the first short axis ultrasonic array probe 24a G1 cross-sectional image or the first vessel cross-sectional image of the blood vessel 20 to be displayed sequentially generate the second short based on the ultrasonic reflected wave signal SG EC received at the second short axis ultrasonic array probe 24b
  • a cross-sectional image of the blood vessel 20 to be displayed in the axial image display region G2, that is, a second blood vessel cross-sectional image is sequentially generated, and based on the reflected wave signal SG EC of the ultrasonic wave received by the long-axis ultrasonic array probe 24c.
  • a longitudinal cross-sectional image of the blood vessel 20 to be displayed in the long-axis image display region G3, that is, a blood vessel longitudinal cross-sectional image is sequentially generated.
  • the short-axis image position determining unit 210 includes a center position CR BV (blood vessel cross-sectional center CR) of the blood vessel 20 in each of the first blood vessel cross-sectional image and the second blood vessel cross-sectional image generated by the blood vessel cross-sectional image generating unit 82. BV ). Then, the short axis image position determining means 210 determines the distance from the first short axis ultrasonic array probe 24a to the center of the blood vessel 20 (blood vessel cross-sectional center CR BV ) and the second short axis ultrasonic array probe.
  • CR BV blood vessel cross-sectional center CR
  • the distance from the child 24b to the center of the blood vessel 20 (blood vessel cross-sectional center CR BV ) is equal to each other, and both of the first short-axis image display region G1 and the second short-axis image display region G2
  • the shift (the shift amount and the shift direction) of the blood vessel center position CR BV with respect to the predetermined measurement position PT1, and the multi-axis drive device 26
  • the relationship between the amount of operation for positioning the ultrasonic probe 24 of each actuator at the predetermined measurement position PT1 is experimentally obtained in advance and set in the short-axis image position determining means 210.
  • the short-axis image position determining means 210 recognizes the blood vessel center position CR BV in each of the first blood vessel cross-sectional image and the second blood vessel cross-sectional image and then performs the predetermined measurement position of the blood vessel center position CR BV.
  • the calculated deviation of the blood vessel center position CR BV is within a predetermined allowable range that can be regarded as zero or substantially zero, and the deviation of the blood vessel center position CR BV is zero or within the predetermined allowable range. If it is determined that it is not, the calculated deviation of the blood vessel center position CR BV is calculated from the relationship between the deviation of the previously set blood vessel center position CR BV and the operation amount of each actuator included in the multi-axis drive device 26. Based on this, the operation amount of each actuator is determined and the multi-axis drive device 26 is operated.
  • Short-axis image position determining means 210 the recognition of such vessel center position CR BV, the calculation of the deviation relative to the predetermined measurement position PT1 vascular center position CR BV, and the operation of the multi-axis drive device 26, the blood vessel Until the deviation of the center position CR BV becomes zero or within the predetermined allowable range, for example, every time the blood vessel cross-sectional image generating unit 82 generates the first blood vessel cross-sectional image and the second blood vessel cross-sectional image, it is repeated.
  • the operation of the multi-axis drive device 26 is terminated, and the ultrasonic probe 24 based on the blood vessel cross-sectional image is terminated.
  • the fact that the positioning is completed is output to the ultrasonic probe position correcting means 212.
  • the long-axis image score calculation unit 214 calculates the inner / middle in the blood vessel longitudinal section image based on the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B calculated by the index value calculation unit 84.
  • calculating the calculated value BR FB i.e. overall intelligibility index BR FB indicates the clarity of the image representing the membrane composite IMC for the whole of the front wall BR F and the rear wall BR B.
  • the comprehensive clarity index value BR FB is an average value or a total value of the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B.
  • these index values XCR F , XCR B average that is, index average.
  • the long axis image score calculating unit 214 for example, every time the index value calculating unit 84 calculates the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B , in other words, the blood vessel cross section image generating means 82 every time the reflected wave signal SG EC of ultrasound to generate the blood vessel longitudinal sectional image is received (acquired), sequentially calculates the overall intelligibility index BR FB.
  • the ultrasonic probe position correcting unit 212 calculates a calculated value based on the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B after the positioning of the ultrasound probe 24 by the short-axis image position determination unit 210 is completed.
  • the multi-axis drive device 26 is operated to position the ultrasonic probe 24 so that the comprehensive clarity index value BR FB ) is equal to or greater than a predetermined determination value (index average determination value) LBR1 FB.
  • the position of the probe 24 is finely adjusted.
  • the index average determination value LBR1 FB is experimentally set, for example, to such a value that the intima-media complex IMC in the blood vessel longitudinal cross-sectional image is clear enough to measure the blood vessel diameter.
  • the ultrasound probe position correcting unit 212 determines whether or not the overall clarity index value BR FB calculated by the long-axis image score calculating unit 214 is equal to or greater than the index average determination value LBR1 FB. If the overall clarity index value BR FB is determined not to be greater than or equal to the index average judgment value LBR1 FB , the multi-axis drive device 26 is operated in a direction to increase the overall clarity index value BR FB. The position of the sonic probe 24 is shifted by a predetermined minute amount.
  • the ultrasonic probe position correcting means 212 repeatedly executes the shift of the position of the ultrasonic probe 24 by a predetermined minute amount until the total intelligibility index value BR FB becomes equal to or greater than the index average determination value LBR1 FB.
  • the operation of the multi-axis drive device 26 is terminated, and the fine adjustment of the position of the ultrasonic probe 24 is completed. Is output to the blood vessel diameter measuring means 88.
  • the ultrasonic probe position correcting means 212 when finely adjusting the position of the ultrasonic probe 24, for example, operates the multi-axis drive device 26 to set the position of the ultrasonic probe 24 by the predetermined minute amount.
  • the total clarity index value BR FB decreases with respect to the movement before the movement in the next cycle, the position of the ultrasonic probe 24 is returned before the movement and further moved to the opposite side by the predetermined minute amount.
  • a range equal to or greater than the index average determination value LBR1 FB corresponds to a predetermined target range of the present invention.
  • the blood vessel diameter measuring means 88 performs FMD measurement, in particular, when it has received from the ultrasonic probe position correcting means 212 that the fine adjustment of the position of the ultrasonic probe 24 has been completed.
  • measurement of the blood vessel lumen diameter d 1 (d a , d MAX ) is started.
  • FIG. 18 is a diagram for explaining the main part of the control operation of the blood vessel inspection device 22 (electronic control device 28) of the present embodiment, that is, the control operation for automatically positioning the ultrasonic probe 24 at a position where FMD measurement is possible.
  • FIG. 16 is a flowchart corresponding to FIG. 15 of the first embodiment. Note that SC6, SC7, and SC8 in FIG. 18 are the same as SA2, SA3, and SA4 in FIG. The control operation shown in FIG. 18 is executed alone or in parallel with other control operations.
  • B-mode short axis image data is acquired. That is, the reflected wave signal SG EC ultrasound at a predetermined reflection wave reception interval PC RV in the longitudinal direction in the first ultrasonic detector array 24a for minor axis is received is scanned (acquired), a second short axis the longitudinal direction at a predetermined reflection wave reception interval PC RV ultrasonic reflected wave signal SG EC is scanned is received (acquired) in ultrasonic detector array 24b.
  • the acquired B-mode short-axis image data (reflected wave signal SG EC ) is stored in a storage device of the electronic control unit 28 or the like.
  • the first said first vessel cross-sectional image based on the reflected wave signal SG EC of the received ultrasonic wave short axis ultrasonic array probe 24a is generated, the ultrasonic array probe for the second short axis the second blood vessel cross-sectional image is generated based on the reflected wave signal SG EC of the received ultrasonic wave 24b.
  • SC2 corresponding to the short axis image position determining means 210, the center position of the blood vessel 20 in the left and right short axis images of the monitor screen display device 30, that is, the first blood vessel cross section image and the second blood vessel cross section image, respectively.
  • CR BV is recognized.
  • SC3 corresponding to the short-axis image position determining unit 210, a shift (shift amount) of the blood vessel center position CR BV with respect to the predetermined measurement position PT1 in each of the first blood vessel cross-sectional image and the second blood vessel cross-sectional image. And the displacement direction) are calculated.
  • SC4 corresponding to the short axis image position determining means 210, it is determined whether or not the deviation of the blood vessel center position CR BV calculated in SC3 is zero or substantially zero. If the determination in SC4 is affirmative, that is, if the deviation of the blood vessel center position CR BV is zero or substantially zero, the process proceeds to SC6. On the other hand, if the determination at SC4 is negative, the operation goes to SC5.
  • each actuator that is, the probe holder motor included in the multi-axis drive device 26 is controlled so as to reduce the deviation of the blood vessel center position CR BV .
  • the multi-axis drive device 26 is operated in a direction to reduce the deviation based on the deviation of the blood vessel center position CR BV based on the relationship set experimentally in advance. After SC5, the process returns to SC1.
  • SC10 corresponding to the ultrasonic probe position correcting means 212, it is determined whether or not the comprehensive clarity index value BR FB calculated in SC9 is equal to or greater than the index average determination value (threshold value) LBR1 FB . If the determination in SC10 is affirmative, that is, if the comprehensive clarity index value BR FB is equal to or greater than the index average determination value LBR1 FB , the process proceeds to SC12. On the other hand, if the determination of SC10 is negative, the process proceeds to SC11.
  • each actuator included in the multi-axis drive device 26 is controlled in a direction to increase the total intelligibility index value BR FB .
  • each actuator probe holder motor included in the multi-axis drive device 26 is controlled in a direction to increase the total intelligibility index value BR FB .
  • the process returns to SC6.
  • FMD measurement specifically, measurement of the blood vessel lumen diameter d 1 (d a , d MAX ) is started.
  • the short-axis image position determining means 210 uses the distance from the first short-axis ultrasonic array probe 24a to the center of the blood vessel 20 (blood vessel cross-sectional center CR BV ) and the second short-axis image position determining means 210.
  • the first short-axis image display region G1 and the second short-axis image display region 30 of the monitor screen display device 30 are set so that the distance from the ultrasonic array probe 24b to the center of the blood vessel 20 (blood vessel cross-sectional center CR BV ) is equal to each other.
  • the multi-axis drive device 26 is operated to position the ultrasonic probe 24 so that the image of the blood vessel 20 is positioned at the center in the width direction thereof.
  • the ultrasonic probe position correcting means 212 calculates based on the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B after the positioning of the ultrasound probe 24 by the short axis image position determination means 210 is completed.
  • the multi-axis drive device 26 is operated to position the ultrasonic probe 24 so that the value (total clarity index value BR FB ) is equal to or greater than a predetermined index average determination value LBR1 FB . Therefore, the operation burden on the operator can be reduced.
  • the operator finely adjusts the position of the ultrasonic probe 24 so as to increase the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B , the operator's A clear blood vessel longitudinal cross-sectional image can be obtained even if the degree of skill is even lower. Further, it is possible to realize automation of FMD measurement including the clarification of the blood vessel longitudinal section image by fine adjustment of the position of the ultrasonic probe 24.
  • the front wall image clarity index value XCR F is calculated by the equation (1)
  • the rear wall image clarity index value XCR B is calculated by the equation (2).
  • the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B may be calculated by other methods.
  • each of the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B is displayed on the monitor screen display device 30 as the sector graph.
  • Other display methods such as a bar graph as illustrated in FIG.
  • each of the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B is numerically displayed and also displayed on the monitor screen display device 30 as a figure (fan graph). There is no problem even if only one of the graphic displays is used.
  • each of the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B is a monitor screen that displays the blood vessel longitudinal section image and the blood vessel transverse section image. Although displayed on the display device 30, it may be displayed on a display device independent of the monitor screen display device 30.
  • the FMD reliability index value XCR FMD has an average value of the front wall image clarity index XCR F and the rear wall image clarity index value XCR B during measurement of the resting diameter d a
  • the average value of the front wall image articulation index value XCR F and the rear wall image articulation index value XCR B at the time of measuring the maximum diameter d MAX is further averaged by the above equation (3).
  • it is not limited to the calculation.
  • it may be a total value of all the index values XCR F and XCR B.
  • the reflected wave recognition control unit 86 is executing the reflected wave recognition control to the reflected wave signal SG EC as shown in FIG. 10, prior to its execution, the reflected wave signal SG A sharpening filter well known for EC may be applied, and sharpening processing may be performed to make the peak and valley of the reflected wave signal SG EC more prominent and clear.
  • a sharpening filter well known for EC may be applied, and sharpening processing may be performed to make the peak and valley of the reflected wave signal SG EC more prominent and clear.
  • the unsharp masking method is a method of sharpening an image by subtracting a second-order differential image (Laplacian image) of the image from a blurred image.
  • the selective image sharpening method is a method of selectively sharpening an edge portion of an image while suppressing the influence of noise by extracting only a Laplacian image for the edge portion of the image.
  • the present invention observes an IMT (Intima-media thickness) test and the carotid artery with ultrasound.
  • the present invention can also be effectively applied to such a control device.
  • the front wall image clarity index value XCR F and the rear wall image clarity index value XCR B are calculated based on the execution result of the reflected wave recognition control. As long as it is calculated in association with the clarity of the inner-inner membrane complex IMC, it may be calculated by other methods.
  • the electronic control device 28 includes the blood vessel diameter measuring means 88, but an electronic control device 28 that does not include the blood vessel diameter measuring means 88 is also conceivable.
  • the upper arm 16 in FIG. 1 is, for example, the upper arm of a human body.

Abstract

 熟練の程度が低いオペレータの操作であっても明瞭な血管の超音波画像を得ることが可能な超音波血管検査装置を提供する。 血管断面画像生成手段82は、生体14の表皮上に配置された超音波プローブ24を用いてその表皮下の血管縦断面画像を逐次生成する。そして、指標値算出手段84は、その血管縦断面画像内の前壁における内中膜複合体を表す画像の明瞭度を示す前壁画像明瞭度指標値XCRFと、血管縦断面画像内の後壁における内中膜複合体を表す画像の明瞭度を示す後壁画像明瞭度指標値XCRBとを算出する。従って、オペレータは、その前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBに基づき、上記血管縦断面画像内の前壁および後壁がより明瞭な画像となるように超音波プローブ24を位置決めすることが可能となる。その結果、そのオペレータの熟練の程度が低くても明瞭な血管縦断面画像を得ることが可能となる。

Description

超音波血管検査装置
 本発明は、超音波プローブを用いて血管断面画像を生成する技術に関するものである。
 複数個の超音波発振子が直線的に配列された超音波アレイ探触子を用いて生体の皮膚下の血管(動脈等)の径などを測定することが行われている。たとえば、特許文献1では、互いに平行な第1および第2の超音波アレイ探触子とそれらの中央部を連結する第3の超音波アレイ探触子とから成るH型の超音波プローブを用い、その第3の超音波アレイ探触子を前記血管に平行な状態でその血管の中心線上に位置させることにより、その血管の内腔径や内中膜複合体の厚みなどを計測する超音波血管検査装置が提案されている。そのような超音波血管検査装置の一例を示す特許文献1の装置は、前記超音波プローブの位置決め状態を制御する多軸駆動装置を備えており、前記第1の超音波アレイ探触子による第1の短軸超音波画像と前記第2の超音波アレイ探触子による第2の短軸超音波画像とに基づき、前記第1の超音波アレイ探触子から前記血管の中心までの距離と前記第2の超音波アレイ探触子から前記血管の中心までの距離とが互いに等しくなるように、前記多軸駆動装置を作動させて前記超音波プローブを位置決めする。更に、前記第1および第2の短軸超音波画像の何れでもそれらの幅方向中央部に前記血管の画像を位置させるように、前記多軸駆動装置を作動させて前記超音波プローブを位置決めする。
特開2009-089911号公報
 しかし、実際には、前記第3の超音波アレイ探触子の直下の血管は湾曲している場合もあり、前記特許文献1の超音波血管検査装置により自動的に超音波プローブが位置決めされただけでは、必ずしも前記第3の超音波アレイ探触子による長軸超音波画像が血管径の計測可能な明瞭な画像にならない場合があった。そのため、例えば、自動的に超音波プローブが位置決めされた後に、オペレータが手動操作で、前記長軸超音波画像を血管径の計測可能な明瞭な画像になるように、超音波プローブの位置を微調整していた。そして、このようなオペレータの手動操作による超音波プローブの位置決め(微調整)作業は熟練を必要とするとともに作業能率を低いものとしているという課題があった。
 本発明は、以上の事情を背景として為されたものであり、その目的とするところは、熟練の程度が低いオペレータの操作であっても明瞭な血管の超音波画像を能率良く得ることが可能な超音波血管検査装置を提供することにある。
 前記目的を達成するための本発明の要旨とするところは、(a)生体の表皮上に配置された超音波プローブを用いて超音波の反射波信号に基づきその生体の表皮下の血管の縦断面画像を生成する血管縦断面画像生成手段を備えた超音波血管検査装置であって、(b)前記血管の縦断面画像内でその血管の内中膜複合体を表す画像の明瞭度を示す指標値を算出する指標値算出手段を含むことにある。
 このようにすれば、超音波血管検査装置のオペレータは、前記血管の縦断面画像(血管縦断面画像)からその画像の明瞭度を直接に判断する必要が無く、前記血管の内中膜複合体を表す画像の明瞭度を示す指標値から客観的にその画像の明瞭度を判断できるので、その指標値がより改善するように超音波プローブの位置が修正されることにより、そのオペレータの熟練の程度が低くても明瞭な血管縦断面画像を能率良く得ることが可能である。
 ここで、好適には、前記指標値算出手段は、前記血管の縦断面画像内の一対の血管壁断面のうち前記超音波プローブに近い側の血管壁断面である前壁における前記内中膜複合体を表す画像の明瞭度を示す前壁画像明瞭度指標値と、前記一対の血管壁断面のうち前記超音波プローブから遠い側の血管壁断面である後壁における前記内中膜複合体を表す画像の明瞭度を示す後壁画像明瞭度指標値とを算出する。このようにすれば、前記血管縦断面画像内の前記前壁および前記後壁の各々がより明瞭な画像となるように、前記超音波プローブを位置決めすることが可能である。
 また、好適には、(a)前記超音波プローブの前記血管の長手方向の相互に異なる複数の受信位置で受信された複数の前記反射波信号の各々に対し、その反射波信号の大きさと前記血管の径方向位置との関係において、その反射波信号の大きさが予め定められた第1ピーク判定閾値を超える第1ピークと、その第1ピークの発生位置よりも前記血管の径方向外側に発生し前記反射波信号の大きさが予め定められた谷判定閾値よりも小さい谷と、前記第1ピークの発生位置から前記谷を介して前記血管の径方向外側の予め定められたピーク間隔閾値を超えない範囲内で発生し前記反射波信号の大きさが予め定められた第2ピーク判定閾値を超える第2ピークとを検出する反射波認識制御を、前記前壁および前記後壁のそれぞれで実行する反射波認識制御手段が設けられており、(b)前記指標値算出手段は、前記反射波認識制御手段が前記前壁において前記反射波認識制御の実行により前記複数の反射波信号のうち前記第1ピークと前記谷と前記第2ピークとの全部を検出した反射波信号の数に基づいて前記前壁画像明瞭度指標値を算出し、前記反射波認識制御手段が前記後壁において前記反射波認識制御の実行により前記複数の反射波信号のうち前記第1ピークと前記谷と前記第2ピークとの全部を検出した反射波信号の数に基づいて前記後壁画像明瞭度指標値を算出する。このようにすれば、前記前壁画像明瞭度指標値と前記後壁画像明瞭度指標値とがばらつきの無い基準に基づき算出されるので、複数の血管検査の相互間で対比可能な前記前壁画像明瞭度指標値と前記後壁画像明瞭度指標値とを得ることができる。
 また、好適には、前記反射波認識制御手段は、予め定められた前記血管の長手方向の観測対象範囲内で受信された前記反射波信号に対して、前記反射波認識制御を実行する。このようにすれば、前記血管縦断面画像の生成のために前記超音波プローブで受信された複数の前記反射波信号の全てに対して上記反射波認識制御が実行される場合に比較して、前記前壁画像明瞭度指標値と前記後壁画像明瞭度指標値とを算出する際の制御負荷を軽減することが可能である。
 また、好適には、前記反射波信号の大きさとは、その反射波信号の振幅またはその反射波信号の振幅を前記血管の縦断面画像を表示するための輝度に置き換えたときのその輝度の高さである。このようにすれば、前記第1ピーク、前記谷、および前記第2ピークの検出が容易になる。
 また、好適には、(a)前記超音波プローブは、複数個の超音波発振子が直線的に前記血管の長手方向とは直交して配列された互いに平行な一対の第1短軸超音波アレイ探触子および第2短軸超音波アレイ探触子と、その第1短軸超音波アレイ探触子およびその第2短軸超音波アレイ探触子の一方又は両方の中央部に隣接して設けられ複数個の超音波発振子が直線的に前記血管の長手方向に配列された長軸用超音波アレイ探触子とを一平面に備えており、(b)前記血管縦断面画像生成手段は、前記長軸用超音波アレイ探触子で受信された前記超音波の反射波信号に基づき前記血管の縦断面画像を生成する。このようにすれば、実用化されている超音波プローブを用いて、前記血管縦断面画像を生成できる。
 また、好適には、前記超音波血管検査装置は(a)前記第1短軸用超音波アレイ探触子による超音波画像を表示する第1短軸画像表示領域と、前記第2短軸用超音波アレイ探触子による超音波画像を表示する第2短軸画像表示領域と、前記血管の縦断面画像を表示する長軸画像表示領域とを有する画像表示装置と、(b)前記超音波プローブの位置決め状態を制御する多軸駆動装置とが設けられており、(c)前記第1短軸用超音波アレイ探触子から前記血管の中心までの距離と前記第2短軸用超音波アレイ探触子から前記血管の中心までの距離とが互いに等しくなるように、且つ、前記第1短軸画像表示領域および前記第2短軸画像表示領域の何れでもそれらの幅方向中央部に前記血管の画像を位置させるように、前記多軸駆動装置を作動させて前記超音波プローブを位置決めする短軸画像位置確定手段と、(d)前記短軸画像位置確定手段による前記超音波プローブの位置決め完了後に、前記前壁画像明瞭度指標値および前記後壁画像明瞭度指標値に基づく算出値が予め定められた目標範囲内になるように、前記多軸駆動装置を作動させて前記超音波プローブを位置決めする超音波プローブ位置修正手段とを、含む。このようにすれば、オペレータの操作負担を軽減することができ、そのオペレータの熟練の程度が更に低くても明瞭な血管縦断面画像を得ることが可能である。
 また、好適には、(a)前記血管の阻血開放前にその血管の安静径を予め測定し、阻血開放後に前記血管の最大径を測定し、前記血管の安静径に対する阻血開放後の径変化割合の最大値を算出する血管径測定手段が設けられており、(b)前記指標値算出手段は、前記安静径の測定時における前記前壁画像明瞭度指標値および前記後壁画像明瞭度指標値と、前記最大径の測定時における前記前壁画像明瞭度指標値および前記後壁画像明瞭度指標値とに基づいて、前記血管径測定手段により算出された前記血管の安静径に対する阻血開放後の径変化割合の最大値の信頼度を示す指標値を算出する。このようにすれば、前記血管の安静径に対する阻血開放後の径変化割合の最大値を複数対比する場合に、その最大値毎に得られた上記信頼度を示す指標値を用いることで、測定精度の悪いデータ(最大値)を排除することが可能となり、例えば、FMD評価結果の信頼性をより引き上げることが可能である。
 また、好適には、前記指標値算出手段は、前記前壁画像明瞭度指標値および前記後壁画像明瞭度指標値の各々を、その指標値の大きさに応じて連続的に変化し相互に対比可能な画像として画像表示装置に表示させる。このようにすれば、オペレータは、それらの指標値が数値で表示されている場合と比較して、直感的に画像の明瞭度を判断でき前記前壁および前記後壁の画像をより迅速に明瞭なものとすることが可能となる。
本発明の一実施例である超音波血管検査装置の全体的な構成を表した図である。 図1の超音波血管検査装置に用いられる超音波プローブの血管に対する姿勢を表すためのxyz軸直交座標軸を説明する図である。 図2の超音波プローブから超音波が放射される計測対象である血管の多層膜構成を説明するための拡大図である。 血管の超音波画像が生成される際に所定の計測位置に位置決めされた図2の超音波プローブと血管との位置関係を示し、そのときのモニタ画面表示装置に表示される血管の超音波画像を示した図である。 図2の超音波プローブからの超音波で測定される阻血開放後の血管内腔径の変化を例示したタイムチャートである。 図1のモニタ画面表示装置の長軸画像表示領域に表示された血管縦断面画像の一例を示した図であって、血管壁の内中膜複合体の明瞭度が高い画像を示す図である。 図1のモニタ画面表示装置の長軸画像表示領域に表示された血管縦断面画像の一例を示した図6とは別の図であって、血管壁の内中膜複合体の明瞭度が低い画像を示す図である。 図1の超音波血管検査装置に含まれる電子制御装置に備えられた制御機能の要部を説明する第1実施例の機能ブロック線図である。 図2の長軸用超音波アレイ探触子で受信されたその受信位置が血管長手方向で相互に異なる複数の反射波信号を例示した図である。 図9に示された複数の反射波信号のうちの1本の反射波信号を、信号振幅を縦軸とし血管の径方向位置を横軸とする2次元座標系に示した図である。 図10と血管縦断面画像との相互関係を説明するために、図1のモニタ画面表示装置の長軸画像表示領域に表示された血管縦断面画像を模式的に示した図である。 図10と同一の反射波信号を示した図であって、その反射波信号の波形が示す血管組織を説明する図である。 図1のモニタ画面表示装置に表示された血管縦断面画像と血管横断面画像と前壁画像明瞭度指標値および後壁画像明瞭度指標値とを例示する図である。 図13の前壁画像明瞭度指標値および後壁画像明瞭度指標値の表示部分を拡大した図である。 図8電子制御装置の制御作動の要部、すなわち、血管縦断面画像を生成し血管の内中膜複合体を表す画像の明瞭度を示す指標値を算出する制御作動を説明するための第1実施例のフローチャートである。 図8電子制御装置の制御作動の図15とは別の要部、すなわち、FMD評価の計測精度を示す指標値を算出する制御作動を説明するための第1実施例のフローチャートである。 図1の超音波血管検査装置に含まれる電子制御装置に備えられた制御機能の要部を説明する第2実施例の機能ブロック線図であって、図8に相当する図である。 図17電子制御装置の制御作動の要部、すなわち、超音波プローブをFMD計測可能な位置に自動的に位置決めする制御作動を説明するための第2実施例のフローチャートであって、図15に相当する図である。 図14においてモニタ画面表示装置に表示されている前壁画像明瞭度指標値および後壁画像明瞭度指標値の表示を示した図であって、図14とは異なる別の表示パターン例を示した図である。
 以下、本発明の実施例を図面を参照しつつ詳細に説明する。
 図1は、センサ保持器10に保持されたプローブユニット12を用いて、生体14の上腕16の皮膚18(厳密には表皮)の上からその皮膚18直下に位置する動脈等の血管20の非侵襲的な超音波診断を行う超音波血管検査装置22(以下、「血管検査装置22」という)の全体的な構成を表した図である。
 プローブユニット12は、血管20に関連する生体情報すなわち血管パラメータを検出するためのセンサとして機能するものであって、互いに平行な一対の第1短軸用超音波アレイ探触子24aおよび第2短軸用超音波アレイ探触子24bとそれらの長手方向中央部を連結する長軸用超音波アレイ探触子24cとを一平面上すなわち平坦な探触面27に有して成るH型の超音波プローブ24と、その超音波プローブ24をxyz方向において位置決めし且つx軸及びz軸まわりの回転角度を位置決めするための多軸駆動装置(位置決め装置)26とを備えている。それら第1短軸用超音波アレイ探触子24a、第2短軸用超音波アレイ探触子24b、および長軸用超音波アレイ探触子24cは、たとえば圧電セラミックスから構成された多数個の超音波振動子(超音波発振子)a~aが直線的に配列されることにより長手状にそれぞれ構成されている。
 図2は、本実施例で用いられるxyz軸直交座標軸を説明するためのものであり、第1短軸用超音波アレイ探触子24aの長手方向と平行でその第1短軸用超音波アレイ探触子24aの直下に位置し血管20又はその付近を通る方向をx軸とし、長軸用超音波アレイ探触子24cの長手方向と平行でx軸と直交する方向をy軸とし、第1短軸用超音波アレイ探触子24aの長手方向と長軸用超音波アレイ探触子24cの長手方向との交点を通り且つ前記x軸方向およびy軸方向に直交する方向をz軸とする。超音波プローブ24は、多軸駆動装置26によりx軸方向に並進、および、x軸およびz軸まわりに回動させられるようになっている。
 図3に示すように、たとえば上腕動脈である血管20は、内膜L、中膜L、外膜Lから成る3層構造を備えている。超音波の反射は音響インピーダンスの異なる部分で発生することから、実際は血管内腔の血液と内膜Lの境界面、および中膜Lと外膜Lとの境界面が白く表示され、組織が白黒の班で表示される。
 図1に戻って、血管検査装置22は、RAMの一時記憶機能を利用しつつ予めROMに記憶されたプログラムに従って入力信号を処理するCPUを有する所謂マイクロコンピュータから構成された電子制御装置28と、モニタ画面表示装置(画像表示装置)30と、超音波駆動制御回路32と、3軸駆動モータ制御回路34とを備えている。上記電子制御装置28によって超音波駆動制御回路32から駆動信号が供給されてプローブユニット12の超音波プローブ24の第1短軸用超音波アレイ探触子24a、第2短軸用超音波アレイ探触子24b、および長軸用超音波アレイ探触子24cからよく知られたビームフォーミング駆動によりビーム状の超音波が順次放射され、その第1短軸用超音波アレイ探触子24aおよび第2短軸用超音波アレイ探触子24bおよび長軸用超音波アレイ探触子24cにより検知された超音波反射信号を受けてその超音波反射信号の処理が行われることによって、皮膚18下の超音波画像が発生させられモニタ画面表示装置30に表示される。
 ここで、モニタ画面表示装置30は、図4(a)に示すように、第1短軸用超音波アレイ探触子24aによる超音波画像を表示する第1短軸画像表示領域G1と、第2短軸用超音波アレイ探触子24bによる超音波画像を表示する第2短軸画像表示領域G2と、長軸用超音波アレイ探触子24cによる超音波画像(血管縦断面画像)を表示する長軸画像表示領域G3とを有している。さらには、それらの第1短軸画像表示領域G1と第2短軸画像表示領域G2と長軸画像表示領域G3とは、皮膚18からの深さ寸法を示す共通の縦軸を備えたものである。なお、図4(a)内の「ImA,ImB」はそれぞれ血管20の横断面である。
 また、モニタ画面表示装置30は、FMD(血流依存性血管拡張反応)の評価に際しては、その内膜の径の変化率すなわち内腔径の拡張率Rを時系列的に表示する。
 また、上記FMDの評価および血管20の超音波画像が生成されるに際して、超音波プローブ24は、血管20に対して所定の計測位置PT1となるよう電子制御装置28によって3軸駆動モータ制御回路34から駆動信号を供給された多軸駆動装置26が駆動することにより位置決めさせられる。上記所定の計測位置PT1とは、上記第1短軸用超音波アレイ探触子24aおよび第2短軸用超音波アレイ探触子24bが血管20に対して直交する位置、且つ、長軸用超音波アレイ探触子24cが血管20に対して平行となる位置である。図4(a)(b)を用いて説明すれば、上記所定の計測位置PT1とは、その図4において「a=b,c=d,e=f」となる位置であり、すなわち、上記第1短軸用超音波アレイ探触子24aから血管20の中心までの距離と上記第2短軸用超音波アレイ探触子24bから血管20の中心までの距離とが互いに等しく、且つ、第1短軸画像表示領域G1および第2短軸画像表示領域G2の何れでもそれらの幅方向中央部に血管20の画像が位置させられた計測位置である。
 センサ保持器10は、三次元空間内の所望の位置すなわち所定の計測位置PT1において生体14の上腕16の皮膚18の上からその皮膚18直下に位置する血管20を変形させない程度に軽く接触させる状態でプローブユニット12を所望の姿勢で保持する。上記プローブユニット12の超音波プローブ24の端面と皮膚18との間には、通常、超音波の減衰、境界面における反射や散乱を抑制して超音波画像を明瞭とするためのよく知られたゼリー、オリーブ油、グリセリン等のカップリング剤や、水を樹脂製袋内に閉じ込めた水袋が介在させられる。
 上記センサ保持器10は、たとえば磁気的吸着力により机、台座等に固定されるマグネット台36と、前記プローブユニット12が固定されるユニット固定具38と、マグネット台36およびユニット固定具38に一端が固定され且つ球状に形成された先端部42を備えた連結部材44、45と、それら連結部材44、45を介して、マグネット台36とユニット固定具38とを相対移動可能に連結し支持する自在アーム40とを備えている。上記自在アーム40は、相互に回動可能に連結された2つのリンク46、47と、そのリンク46、47の一端にて前記各先端部42に所定の抵抗が付勢されつつその先端部42に対して回曲可能に嵌め入れられた勘合穴48をそれぞれ有する回曲関節部50、51と、各リンク46、47の他端にてその他端を相互に相対回動可能に連結し且つその連結箇所を貫設するねじ穴に螺合されたおねじ付き固定ノブ52が締め付けられることで得られる締着力により相対回動不能にされる回動関節部54とを、有する。
 多軸駆動装置26は、x軸回動アクチュエータにより超音波プローブ24のx軸まわりの回動位置を位置決めするためにユニット固定具38に固定されるx軸回動(ヨーイング)機構と、x軸回動アクチュエータによって超音波プローブ24のx軸方向の並進位置を位置決めするためのx軸並進機構と、z軸アクチュエータにより超音波プローブ24のz軸まわりの回動位置を位置決めするためのz軸回動機構とから構成されている。多軸駆動装置26は、このような構成により、超音波プローブ24の位置決め状態を制御する。
 図1において、超音波駆動制御回路32は、電子制御装置28からの指令に従って、たとえば上記第1短軸用超音波アレイ探触子24aを構成する一列に配列された多数個の超音波振動子a乃至aのうち、その端の超音波振動子aら、一定数の超音波振動子群たとえば15個のa乃至a15毎に所定の位相差を付与しつつ10MHz程度の周波数で同時駆動するビームフォーミング駆動することにより超音波振動子の配列方向において収束性の超音波ビームを血管20に向かって順次放射させ、超音波振動子を1個ずつずらしながらその超音波ビームをスキャン(走査)させたときの放射毎の反射波を受信して電子制御装置28へ入力させる。
 電子制御装置28は、上記反射波に基づいて画像を合成し、皮膚18下における血管20の横断面画像(短軸画像)、あるいは縦断面画像(長軸画像)を生成させて、モニタ画面表示装置(画像表示装置)30にそれぞれ表示させる。また、その画像から、血管20の径或いは内皮70の直径である内皮径(内腔径)dなどが算出される。また、血管内皮機能を評価するために、虚血反応性充血後のFMD(血流依存性血管拡張反応)を表す血管内腔径の拡張率(変化率)R(%)[=100×(d-d)/d]が算出される。なお、上記式の「d」は安静時の血管内腔径(ベース径、安静径)を示している。
 図5は、阻血(駆血)開放後の血管内腔径dの変化を例示したタイムチャートである。図5では、t1時点が阻血開放時を表しており、t2時点から血管内腔径dが拡張し始め、t3時点で血管内腔径dがその最大値dMAXに達していることが示されている。従って、電子制御装置28が算出する血管内腔径の拡張率Rは、t3時点で最大になる。
 FMD評価のための前記阻血は、図1に示すように、電子制御装置28が備えるカフ圧制御部56(カフ圧制御手段56)が空気ポンプ58からの元圧を圧力制御弁60で制御して上腕16に巻回されたカフ62に供給し、そのカフ62の圧力(カフ圧)を生体14の最高血圧を超える所定の阻血カフ圧にまで昇圧することにより行われる。このとき、上記カフ圧制御部56は上記カフ圧を検出するための圧力センサ64からの信号でそのカフ圧を検出する。そして、図5においては、例えば、カフ圧制御部56は、阻血開放前の所定時間すなわちt1時点前の所定時間にわたって上記カフ圧を上記阻血カフ圧で維持し、阻血開放時(t1時点)に上記カフ圧を直ちに大気圧にまで減圧する。
 本実施例では、多軸駆動装置26が駆動されことにより、超音波プローブ24が前記所定の計測位置PT1に位置決めされるので、基本的には、図6の長軸画像表示領域G3に表示された血管縦断面画像のように、血管20の内膜Lと中膜Lとを併せた内中膜複合体IMC(intima-media complex)の明瞭度が高い画像が得られる。この図6について説明すると、上記血管縦断面画像内の一対の血管壁断面のうち超音波プローブ24に近い側の血管壁断面である前壁BRFと超音波プローブ24から遠い側の血管壁断面である後壁BRBとの何れでも内中膜複合体IMCの低エコー帯(図6の黒いライン)BLIMCが血管20の長手方向に連続性を有して明瞭に表示されているので、図6は内中膜複合体IMCの明瞭度が高い画像であると言える。このような明瞭な画像に基づきFMD評価のような血管画像診断がなされれば、それは十分に精度の高い画像診断であると言える。
 一方で、超音波プローブ24が前記所定の計測位置PT1に位置決めされたとしても、血管20が湾曲している等に起因して、図6のような内中膜複合体IMCの明瞭度が高い血管縦断面画像が得られない場合もある。例えば、図7のような内中膜複合体IMCの明瞭度が低い血管縦断面画像となる場合である。図7では、図6と異なり、前壁BRFおよび後壁BRBの何れでも血管内腔と血管壁との境界が血管20の長手方向に不連続であり、更に、それらの内中膜複合体IMCの低エコー帯BLIMCが殆ど画像から認識できず明瞭ではない。この図7のような不明瞭な内中膜複合体IMCの血管縦断面画像である場合には、例えばオペレータが超音波プローブ24の位置を微調整して前壁BRFおよび後壁BRBの内中膜複合体IMCをより明瞭にする操作を行う。本実施例の血管検査装置22はその微調整を支援する制御機能を有しており、その制御機能の要部を図8を用いて説明する。
 図8は、血管検査装置22に備えられた制御機能の要部を説明する機能ブロック線図である。図8に示すように、電子制御装置28に設けられた血管画像評価部100(図1参照)は、画像生成実行判断部としての画像生成実行判断手段80と、血管断面画像生成部としての血管断面画像生成手段82と、指標値算出部としての指標値算出手段84と、血管径測定部としての血管径測定手段88とを備えている。
 画像生成実行判断手段80は、血管20の超音波画像を超音波プローブ24を用いて取得するか否かをを判断する。すなわち、その血管20の超音波画像の取得開始と取得終了とを判断する。
 例えば、血管20の超音波画像を取得するために操作されるスイッチが血管検査装置22に設けられている場合には、画像生成実行判断手段80は、そのスイッチがONに切り替えられると上記超音波画像の取得を開始すると判断し、そのスイッチがOFFに切り替えられると上記超音波画像の取得を終了すると判断する。
 血管断面画像生成手段82は、血管縦断面画像生成手段として機能し、生体14の表皮上に配置された超音波プローブ24を用いて超音波の反射波信号SGECに基づき生体14の表皮下の血管20の縦断面画像(血管縦断面画像)を逐次生成する。具体的には、長軸用超音波アレイ探触子24cで順次且つ繰返し受信された超音波の反射波信号SGECに基づき上記血管縦断面画像を逐次生成する。例えば、上記血管20の縦断面画像とは、図6または図7に示されるような長軸画像表示領域G3に表示される画像であり、血管断面画像生成手段82が一般的なBモード法などを用いて生成する。
 ここで、1つの前記血管縦断面画像の生成の際、詳細には、血管断面画像生成手段82は、長軸用超音波アレイ探触子24cで受信された前記超音波の反射波信号SGECに基づき前記血管縦断面画像を生成する。すなわち、血管断面画像生成手段82は、長軸用超音波アレイ探触子24cにおいてその長手方向に所定の反射波受信間隔(ラインピッチ)PCRVで超音波の反射波信号SGECを走査して受信し、その受信された複数の反射波信号SGECに基づいて前記血管縦断面画像を生成する。
 更に、血管断面画像生成手段82は、血管横断面画像生成手段としても機能し、第1短軸用超音波アレイ探触子24aで順次且つ繰返し受信された前記超音波の反射波信号SGECに基づき第1短軸画像表示領域G1に表示させる血管20の横断面画像(血管横断面画像)を逐次生成し、第2短軸用超音波アレイ探触子24bで順次且つ繰返し受信された前記超音波の反射波信号SGECに基づき第2短軸画像表示領域G2に表示させる血管20の横断面画像を逐次生成する。血管断面画像生成手段82は、例えば、画像生成実行判断手段80が上記超音波画像の取得を開始すると判断してから終了すると判断するまで、モニタ画面表示装置30上に表示される画像にちらつきが生じない所定以上の周期TAで上記血管縦断面画像および上記血管横断面画像の生成を繰り返し行う。
 指標値算出手段84は、血管断面画像生成手段82が生成する前記血管縦断面画像内で血管20の内中膜複合体IMCを表す画像の明瞭度を示す指標値を算出する。例えば、血管断面画像生成手段82が血管縦断面画像を生成する毎に、その指標値を逐次算出する。具体的には、その指標値は上記血管縦断面画像内の前壁BRFおよび後壁BRBのそれぞれについて算出される。すなわち、指標値算出手段84は、上記血管縦断面画像内の前壁BRFにおける内中膜複合体IMCを表す画像の明瞭度を示す前壁画像明瞭度指標値XCRFと、上記血管縦断面画像内の後壁BRBにおける内中膜複合体IMCを表す画像の明瞭度を示す後壁画像明瞭度指標値XCRBとを算出する。指標値算出手段84は、その前壁画像明瞭度指標値XCRFと後壁画像明瞭度指標値XCRBとを算出するために、反射波認識制御部としての反射波認識制御手段86を備えている。なお、上記内中膜複合体IMCを表す画像の明瞭度とは、内中膜複合体IMCは血管20の長手方向に沿って存在しているのであるから、上記血管縦断面画像においてその内中膜複合体IMCを表す低エコー帯BLIMC及び高エコー帯の画像が血管20の長手方向において連続的に認識できる度合いであると言える。また、確認的に述べるが、上記前壁画像明瞭度指標値XCRFと上記後壁画像明瞭度指標値XCRBとを総括して表現すれば、それらは前記内中膜複合体IMCを表す画像の明瞭度を示す指標値であると説明できる。
 反射波認識制御手段86は、血管断面画像生成手段82が前記血管縦断面画像を生成するために受信された複数の反射波信号SGECを、例えば血管断面画像生成手段82が血管縦断面画像を生成する毎に取得する。例えば血管断面画像生成手段82から取得する。1つの血管縦断面画像を生成するための上記複数の反射波信号SGECについて見れば、その複数の反射波信号SGECの各々は、長軸用超音波アレイ探触子24cの血管長手方向の相互に異なる受信位置で受信されたものである。すなわち、反射波認識制御手段86は、図9に例示するように、長軸用超音波アレイ探触子24cでの受信位置が血管長手方向で相互に異なる複数の反射波信号SGECを取得する。
 反射波認識制御手段86は、取得した上記複数の反射波信号SGECの各々に対し、図10のような反射波信号SGECの信号振幅AMSGと血管20の径方向位置PTRとの関係において、その反射波信号SGECの信号振幅AMSGが予め定められた第1ピーク判定閾値LT1PKを超える第1ピークPK1と、その第1ピークPK1の発生位置よりも血管20の径方向外側に発生し反射波信号SGECの信号振幅AMSGが予め定められた谷判定閾値LTBTMよりも小さい谷BTMと、前記第1ピークPK1の発生位置から谷BTMを介して血管20の径方向外側の予め定められたピーク間隔閾値LTSPを超えない範囲内で発生し反射波信号SGECの信号振幅AMSGが予め定められた第2ピーク判定閾値LT2PKを超える第2ピークPK2とを検出する反射波認識制御を、前記前壁BRFおよび前記後壁BRBのそれぞれで実行する。この反射波認識制御を図10を用いて具体的に説明する。
 図10は、反射波信号SGECの信号振幅AMSG(信号強度AMSG)を縦軸とし血管20の径方向位置PTRを横軸とする座標系に、1つ(1本)の反射波信号SGECの前壁BRFまたは後壁BRBに該当する部分を表した図である。例えば、図10に表された反射波信号SGECが前壁BRFに該当する部分であるとすれば、図10の矢印AR01は図11の血管縦断面画像における矢印AR02に対応する一方で、図10に表された反射波信号SGECが後壁BRBに該当する部分であるとすれば、図10の矢印AR01は図11の矢印AR03に対応する。なお、図10では、第1ピークPK1、谷BTM、及び、第2ピークPK2のそれぞれで反射波信号SGECの波形が鋭角に尖っているがこれは一例であり、例えば、第1ピークPK1、谷BTM、及び、第2ピークPK2の何れかではその波形が鋭角に尖らずに台形のように略平坦になっている場合もある。
 反射波認識制御手段86は、先ず、反射波信号SGECにおいて血管内腔に該当する部分を認識する。例えば、血管20の横断面画像(短軸画像)から認識してもよいし、前壁BRFおよび後壁BRBのそれぞれを示す反射波信号SGECの該当部分の中間位置を血管内腔であると認識してもよい。また、本実施例では、後述するように、予め定められた血管20の長手方向の観測対象範囲AOB(図11参照)内で受信された反射波信号SGECに対して前記反射波認識制御が実行されるので、上記血管内腔に該当する部分が認識されるのは、前記血管縦断面画像を生成する基となる複数の反射波信号SGECのうち前記観測対象範囲AOBで受信された反射波信号SGECのみでよい。
 そして、前記反射波認識制御において、反射波認識制御手段86は、前壁BRFおよび後壁BRBの何れか一方で、図10の矢印AR01のように血管内腔側から血管外側に向けて反射波信号SGECの信号振幅AMSGを認識して行き、信号振幅AMSGが前記第1ピーク判定閾値LT1PKを超える極大点が存在すればそれを第1ピークPK1として検出する。次に、その第1ピークPK1の発生位置(図10の横軸位置)よりも血管20の径方向外側で信号振幅AMSGが前記谷判定閾値LTBTMよりも小さい極小点が存在すればそれを谷BTMとして検出する。次に、前記第1ピークPK1の発生位置から谷BTMを介して血管20の径方向外側の前記ピーク間隔閾値LTSPを超えない範囲内で信号振幅AMSGが前記第2ピーク判定閾値LT2PKを超える極大点が存在すればそれを第2ピークPK2として検出する。前壁BRFおよび後壁BRBの一方で第1ピークPK1と谷BTMと第2ピークPK2との検出が完了した後、その前壁BRFおよび後壁BRBの他方でも同一の反射波信号SGECに対し同様の検出を行う。反射波認識制御手段86は、このような前壁BRFおよび後壁BRBのそれぞれでの前記反射波認識制御を前記複数の反射波信号SGECの各々に対し実行する。
 ここで、前記反射波認識制御において、反射波信号SGECの信号振幅AMSGが前記第1ピーク判定閾値LT1PK、前記谷判定閾値LTBTM、及び前記第2ピーク判定閾値LT2PKと比較されるが、その比較対象は上記信号振幅AMSGに限られるわけではなく、前記血管縦断面画像と関連する反射波信号SGECの大きさであればよく、例えば、上記信号振幅AMSGが輝度変調されて血管縦断面画像(Bモード長軸画像)を表示するための輝度に置き換えられたときのその輝度の高さなどでもよい。また、前記第1ピーク判定閾値LT1PK、前記谷判定閾値LTBTM、及び前記第2ピーク判定閾値LT2PKは、例えば、第1ピークPK1と谷BTMと第2ピークPK2とのそれぞれに対応する超音波画像が相互に異なった模様として視認できるように実験的に定められており、前記谷判定閾値LTBTMは前記第1ピーク判定閾値LT1PK及び前記第2ピーク判定閾値LT2PKよりも小さい値に設定されている。また、図10では第1ピークPK1よりも第2ピークPK2の方が大きいように図示されているが逆に第1ピークPK1よりも第2ピークPK2の方が小さい場合もあり得るので、前記第1ピーク判定閾値LT1PK及び前記第2ピーク判定閾値LT2PKはそれら相互の大小関係に特に制限はなく、それらの判定閾値LT1PK,LT2PKが相互に異なる値とされていてもよいし同一の値とされていても差し支えない。また、図12に示すように、反射波信号SGECの谷BTMを含む第1ピークPK1と第2ピークPK2との間は前記血管縦断面画像では内中膜複合体IMCの前記低エコー帯BLIMCとして認識されるので、前記ピーク間隔閾値LTSPは、例えば、生体の内中膜複合体IMCの厚さに基づいて実験的に設定されている。また、反射波認識制御手段86は、前記血管縦断面画像の基となる複数の反射波信号SGECの全部に対して前記反射波認識制御を実行してもよいが、本実施例では演算負荷の軽減のため、図11に示すように、予め定められた血管20の長手方向の観測対象範囲AOB内で受信された反射波信号SGECに対して、前記反射波認識制御を実行する。この観測対象範囲AOBは、例えば、長軸用超音波アレイ探触子24cの長手方向幅の一部であり固定されていてもオペレータなどにより設定変更可能とされていてもよく、前記血管縦断面画像内の内中膜複合体IMCの明瞭度を判断できるように実験的に設定されている。
 図10は、前記反射波認識制御の実行により第1ピークPK1と谷BTMと第2ピークPK2との全てが検出される反射波信号SGECの例を示しているが、第1ピークPK1と谷BTMと第2ピークPK2との何れか又は何れもが検出されない反射波信号SGECも存在する。そこで、反射波認識制御手段86は、前記観測対象範囲AOBで受信された複数の反射波信号SGECのうち、上記反射波認識制御の実行により第1ピークPK1と谷BTMと第2ピークPK2との全部を検出した反射波信号SGECの数QLSG、すなわち、ピーク認識ライン数QLSGを、前壁BRFと後壁BRBとの各々について記憶する。例えば、図9に示すように、前記反射波認識制御の実行対象となる反射波信号SGECの数(ライン数、本数)が13本(=観測対象範囲AOB/反射波受信間隔PCRV)であるとする。この場合、反射波認識制御手段86は、その13本の反射波信号SGECの各々に対し前記反射波認識制御を前壁BRFおよび後壁BRBのそれぞれで実行する。そして、図9のNo.3及びNo.4の反射波信号SGECのように破線L01で示す箇所で第1ピークPK1を検出できなかった場合には、これら2本をピーク認識ライン数QLSGに含めずに、ピーク認識ライン数QLSGを11本(=13-2)と記憶する。
 反射波認識制御手段86が前記観測対象範囲AOBで受信された複数の反射波信号SGECに対する前記反射波認識制御の実行を完了した後、指標値算出手段84は、反射波認識制御手段86が前壁BRFにおいて前記反射波認識制御の実行により前記複数の反射波信号SGECのうち前記第1ピークPK1と前記谷BTMと前記第2ピークPK2との全部を検出した反射波信号SGECの数QLSG(ピーク認識ライン数QLSG)に基づいて前記前壁画像明瞭度指標値XCRFを算出する。そして、指標値算出手段84は、反射波認識制御手段86が後壁BRBにおいて前記反射波認識制御の実行により前記複数の反射波信号SGECのうち前記第1ピークPK1と前記谷BTMと前記第2ピークPK2との全部を検出した反射波信号の数QLSG(ピーク認識ライン数QLSG)に基づいて前記後壁画像明瞭度指標値XCRBを算出する。具体的には、指標値算出手段84は、下記式(1)により前記前壁画像明瞭度指標値XCRFすなわち前壁スコアXCRFを算出し、下記式(2)により前記後壁画像明瞭度指標値XCRBすなわち後壁スコアXCRBを算出する。なお、下記式(1)及び下記式(2)において、「XCRF」は上記前壁スコアXCRFを示し、「XCRB」は上記後壁スコアXCRBを示し、「QLFSG」は前壁BRFにおけるピーク認識ライン数QLSGを示し、「QLBSG」は後壁BRBにおけるピーク認識ライン数QLSGを示し、「PCRV」は前記反射波受信間隔(ラインピッチ)PCRVを示し、「AOB」は前記観測対象範囲(観測エリア幅)AOBを示す。
  XCRF=(QLFSG×PCRV)/AOB×100 ・・・(1)
  XCRB=(QLBSG×PCRV)/AOB×100 ・・・(2)
 このようにして、指標値算出手段84は、前記前壁画像明瞭度指標値XCRFおよび前記後壁画像明瞭度指標値XCRBを逐次算出する。上記式(1)、式(2)から判るように、何れの画像明瞭度指標値XCRF,XCRBも、その変化範囲が「0~100」である相対値である。そして、指標値算出手段84は、それらの算出をする毎に、その前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBの各々を、図13の二点鎖線L02で囲んで示すように、モニタ画面表示装置30に前記血管縦断面画像および前記血管横断面画像とともに逐次表示させる。具体的には、図13の二点鎖線L02で囲まれた部分の拡大図である図14に示すように、指標値算出手段84は、前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBの各々を、数値表示するとともに、その指標値XCRF,XCRBの大きさに応じて連続的に変化し相互に対比可能な画像乃至は図形として、モニタ画面表示装置30に表示させる。図14では、前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBの各々は、その大きさが大きくなるほど中心角及び面積が拡大し共通の直線を境とする線対称を成す一対の扇形の扇形グラフ(図14の黒色部分)で表現されており、最大値の「100」になるとその扇形グラフはその面積が最大となって半円形となる。そして、それらの指標値XCRF,XCRBが何れも最大値になると上記半円形は1つの完全な円形を構成する。
 血管径測定手段88は、前記血管縦断面画像から非侵襲的に血管内腔径dを測定する。具体的には、血管径測定手段88は、FMD評価で、血管20の阻血開放後の血管20の径変化割合(血管内腔径dの拡張率R)を測定するために、その血管20の阻血開放前に安静時の血管内腔径d(安静径d)を予め測定する。また、血管径測定手段88は、FMD評価のために、阻血開放後の予め定められた血管径測定期間TIME1内に血管内腔径dを逐次測定し、更に、その測定した血管内腔径dと前記安静径dとから血管20の径変化割合Rを逐次算出する。例えば、血管径測定手段88は、図5に示されるように阻血開放時からその後変化する血管内腔径dを上記血管径測定期間TIME1にわたって時間経過に従って逐次連続的に測定する。或いは、阻血開放時を基準として血管内腔径dが略最大になると予測される測定時点を1又は2以上実験的に求めておき、その測定時点で血管内腔径dを測定してもよい。前記血管径測定期間TIME1は、阻血開放後の血管内腔径dの最大値dMAXを測定するために阻血開放時を基準として実験的に設定され血管径測定手段88に記憶された血管内腔径dの測定期間であり、図5に示すように血管内腔径dがその最大値dMAXに到達する時点(t3時点)を含み且つ阻血開放時(t1時点)から開始する。従って、血管径測定手段88は、阻血開放後に血管20の最大径(最大内腔径)dMAXを測定することになる。
 更に、血管径測定手段88は、上記血管径測定期間TIME1の経過後に、血管20の前記安静径dに対する阻血開放後の最大径変化割合RMAX、すなわち、その安静径dを基準とした阻血開放後の血管20の径変化割合(拡張率R)の最大値RMAX(%)[=100×(dMAX-d)/d]を血管拡張率評価値(%FMD値)として算出する。そして、その算出した血管拡張率評価値RMAXをモニタ画面表示装置30に表示させる。
 指標値算出手段84は、前述の機能に加え、前記安静径dの測定時における前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBと、阻血開放後の前記最大径dMAXの測定時における前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBとを算出する。そして、その安静径dの測定時における前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBと、その最大径dMAXの測定時における前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBとに基づいて、血管径測定手段88により算出された血管20の安静径dに対する阻血開放後の最大径変化割合(血管拡張率評価値)RMAXの信頼度を示す指標値XCRFMDすなわちFMD信頼度指標値XCRFMDを算出する。例えば、指標値算出手段84は、前記安静径dの測定時における前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBの平均値を安静径測定時スコアXCR1として算出し、前記最大径dMAXの測定時における前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBの平均値を最大径測定時スコアXCR2として算出する。そして、その安静径測定時スコアXCR1と最大径測定時スコアXCR2とから下記式(3)により前記FMD信頼度指標値XCRFMDを算出する。指標値算出手段84は、その算出したFMD信頼度指標値XCRFMDを、モニタ画面表示装置30に前記血管拡張率評価値RMAXとともに表示させる。そのFMD信頼度指標値XCRFMDは、前述したように上記血管拡張率評価値RMAXの信頼度を示す指標値であり、換言すれば、その血管拡張率評価値RMAXの計測精度を示す指標値であって、FMD信頼度指標値XCRFMDはその値が大きいほど上記血管拡張率評価値RMAXの計測精度が高いことを示すものである。
  XCRFMD=(XCR1+XCR2)/2 ・・・(3)
 電子制御装置28に設けられた表示制御手段90(表示制御部90)は、図13および図14に示すように、血管断面画像生成手段82が生成した前記血管縦断面画像および前記血管横断面画像をモニタ画面表示装置30に逐次表示し、指標値算出手段84が算出した前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBを数値および図形(扇形グラフ)でモニタ画面表示装置30に逐次表示する。
 また、表示制御手段90は、FMD評価が実施された場合、すなわち、血管径測定手段88が血管拡張率評価値RMAXを算出した場合には、その血管拡張率評価値RMAXと前記FMD信頼度指標値XCRFMDとをモニタ画面表示装置30に表示する。
 図15は、血管検査装置22(電子制御装置28)の制御作動の要部、すなわち、前記血管縦断面画像を生成し血管20の内中膜複合体IMCを表す画像の明瞭度を示す指標値を算出する制御作動を説明するためのフローチャートである。この図15に示す制御作動は、単独で或いは他の制御作動と並列的に実行される。
 先ず、画像生成実行判断手段80に対応するステップ(以下、「ステップ」を省略する)SA1においては、血管20の超音波画像の取得を開始するか否かが判断される。このSA1の判断が肯定された場合、すなわち、前記超音波画像の取得を開始する場合には、SA2へ移る。一方、このSA1の判断が否定された場合には、図15のフローチャートは終了する。
 血管断面画像生成手段82に対応するSA2においては、長軸用超音波アレイ探触子24cにおいてその長手方向に所定の反射波受信間隔PCRVで超音波の反射波信号SGECが走査され受信(取得)される。そして、その受信された複数の反射波信号SGECすなわちBモード長軸画像データが電子制御装置28の記憶装置等に記憶され、その受信された複数の反射波信号SGECに基づいて前記血管縦断面画像が生成される。
 指標値算出手段84および反射波認識制御手段86に対応するSA3では、前記反射波認識制御が、前記観測対象範囲AOBで受信された複数の反射波信号SGECの各々に対し前記前壁BRFにおいて実行される。そして、前記前壁画像明瞭度指標値(前壁スコア)XCRFが、前壁BRFにおいて上記反射波認識制御の実行により上記複数の反射波信号SGECのうち前記第1ピークPK1と前記谷BTMと前記第2ピークPK2との全部を検出した反射波信号SGECの数QLSG(QLFSG)に基づいて算出される。
 指標値算出手段84および反射波認識制御手段86に対応するSA4では、前記反射波認識制御が、前記観測対象範囲AOBで受信された複数の反射波信号SGECの各々に対し前記後壁BRBにおいて実行される。そして、前記後壁画像明瞭度指標値(後壁スコア)XCRBが、後壁BRBにおいて上記反射波認識制御の実行により上記複数の反射波信号SGECのうち前記第1ピークPK1と前記谷BTMと前記第2ピークPK2との全部を検出した反射波信号SGECの数QLSG(QLBSG)に基づいて算出される。
 表示制御手段90に対応するSA5では、図13および図14に示すように、前記SA2にて生成された前記血管縦断面画像とともに、前記SA3で算出された前記前壁スコアXCRFと前記SA4で算出された前記後壁スコアXCRBとがモニタ画面表示装置30にリアルタイムで表示される。例えば、図14に示すように、その前壁スコアXCRFと後壁スコアXCRBとはそれぞれ、0~100の間で変化する数値で表示されるとともに、前壁スコアXCRFと後壁スコアXCRBとの両方が最大値の「100」になれば完全な円形となる一対の扇形グラフで表示される。
 画像生成実行判断手段80に対応するSA6においては、血管20の超音波画像の取得を終了するか否かが判断される。このSA6の判断が肯定された場合、すなわち、上記超音波画像の取得を終了する場合には、図15のフローチャートは終了する。一方、このSA6の判断が否定された場合、すなわち、上記超音波画像の取得を継続する場合には、SA2に戻る。従って、前記SA2からSA5までのステップが、上記超音波画像の取得開始時から終了時まで繰り返し実行される。例えば、極めて短いサイクルで繰り返し実行される。
 図16は、図15とは別の血管検査装置22(電子制御装置28)の制御作動の要部、すなわち、FMD評価の計測精度を示す指標値を算出する制御作動を説明するためのフローチャートである。この図16に示す制御作動は、単独で或いは他の制御作動と並列的に実行される。
 先ず、指標値算出手段84および反射波認識制御手段86に対応するSB1においては、前記反射波認識制御が、FMD評価での前記阻血開放前の安静径dを測定するための前記血管縦断面画像の基となる複数の反射波信号SGECの各々に対し、前壁BRFおよび後壁BRBのそれぞれにおいて実行される。詳細には、その複数の反射波信号SGECの全部に対してではなく、前記観測対象範囲AOBで受信された複数の反射波信号SGECの各々に対し、前記反射波認識制御が実行される。そして、その反射波認識制御の実行結果から上記安静径dの測定時における前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBが算出され、それらの平均値が前記安静径測定時スコアXCR1として算出される。
 指標値算出手段84および反射波認識制御手段86に対応するSB2においては、前記反射波認識制御が、FMD評価での前記阻血開放後の最大径dMAXを測定するための前記血管縦断面画像の基となる複数の反射波信号SGECの各々に対し、前壁BRFおよび後壁BRBのそれぞれにおいて実行される。詳細には前記SB1と同様に、その複数の反射波信号SGECの全部に対してではなく、前記観測対象範囲AOBで受信された複数の反射波信号SGECの各々に対し、前記反射波認識制御が実行される。そして、その反射波認識制御の実行結果から上記最大径dMAXの測定時における前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBが算出され、それらの平均値が前記最大径測定時スコアXCR2として算出される。
 指標値算出手段84に対応するSB3においては、FMD信頼度指標値XCRFMDが、上記安静径測定時スコアXCR1と最大径測定時スコアXCR2とから前記式(3)により算出される。
 表示制御手段90に対応するSB4においては、上記FMD信頼度指標値XCRFMDが、FMD評価結果例えば前記血管拡張率評価値(%FMD値)RMAXとともに、モニタ画面表示装置30に表示される。
 本実施例には次のような効果(A1)乃至(A8)がある。(A1)本実施例によれば、血管断面画像生成手段(血管縦断面画像生成手段)82は、生体14の表皮上に配置された超音波プローブ24を用いて超音波の反射波信号SGECに基づき生体14の表皮下の血管20の縦断面画像を逐次生成する。そして、指標値算出手段84は、血管断面画像生成手段82が生成する上記血管20の縦断面画像内でその血管20の内中膜複合体IMCを表す画像の明瞭度を示す指標値を算出する。従って、血管検査装置22のオペレータは、上記血管20の縦断面画像(血管縦断面画像)からその画像の明瞭度を直接に判断する必要が無く、上記血管20の内中膜複合体IMCを表す画像の明瞭度を示す指標値から客観的にその画像の明瞭度を判断できるので、その指標値がより改善するように超音波プローブ24の位置を容易に微調整することができ、そのオペレータの熟練の程度が低くても明瞭な血管縦断面画像を能率良く得ることが可能である。
 (A2)また、本実施例によれば、指標値算出手段84は、前記血管縦断面画像内の前壁BRFにおける内中膜複合体IMCを表す画像の明瞭度を示す前壁画像明瞭度指標値XCRFと、上記血管縦断面画像内の後壁BRBにおける内中膜複合体IMCを表す画像の明瞭度を示す後壁画像明瞭度指標値XCRBとを算出するので、例えば、オペレータが、その前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBに基づき、上記血管縦断面画像内の前壁BRFおよび後壁BRBの各々がより明瞭な画像となるように能率良く超音波プローブ24を位置決めすることが可能である。
 (A3)また、本実施例によれば、反射波認識制御手段86は、長軸用超音波アレイ探触子24cの血管長手方向の相互に異なる受信位置で受信された複数の反射波信号SGECの各々に対し、図10のような反射波信号SGECの信号振幅AMSGと血管20の径方向位置PTRとの関係において、その反射波信号SGECの信号振幅AMSGが予め定められた第1ピーク判定閾値LT1PKを超える第1ピークPK1と、その第1ピークPK1の発生位置よりも血管20の径方向外側に発生し反射波信号SGECの信号振幅AMSGが予め定められた谷判定閾値LTBTMよりも小さい谷BTMと、前記第1ピークPK1の発生位置から谷BTMを介して血管20の径方向外側の予め定められたピーク間隔閾値LTSPを超えない範囲内で発生し反射波信号SGECの信号振幅AMSGが予め定められた第2ピーク判定閾値LT2PKを超える第2ピークPK2とを検出する前記反射波認識制御を、前記前壁BRFおよび前記後壁BRBのそれぞれで実行する。そして、指標値算出手段84は、反射波認識制御手段86が前壁BRFにおいて前記反射波認識制御の実行により前記複数の反射波信号SGECのうち前記第1ピークPK1と前記谷BTMと前記第2ピークPK2との全部を検出した反射波信号SGECの数QLSGに基づいて前記前壁画像明瞭度指標値XCRFを算出し、反射波認識制御手段86が後壁BRBにおいて前記反射波認識制御の実行により前記複数の反射波信号SGECのうち前記第1ピークPK1と前記谷BTMと前記第2ピークPK2との全部を検出した反射波信号の数QLSGに基づいて前記後壁画像明瞭度指標値XCRBを算出する。従って、その前壁画像明瞭度指標値XCRFと後壁画像明瞭度指標値XCRBとがばらつきの無い基準に基づき算出されるので、複数の血管検査の相互間で対比可能な上記前壁画像明瞭度指標値XCRFと上記後壁画像明瞭度指標値XCRBとを得ることができる。
 (A4)また、本実施例によれば、反射波認識制御手段86は、予め定められた血管20の長手方向の観測対象範囲AOB内で受信された反射波信号SGECに対して、前記反射波認識制御を実行するので、前記血管縦断面画像の生成のために長軸用超音波アレイ探触子24cで受信された複数の反射波信号SGECの全てに対して上記反射波認識制御が実行される場合に比較して、前記前壁画像明瞭度指標値XCRFと前記後壁画像明瞭度指標値XCRBとを算出する際の制御負荷を軽減することが可能である。
 (A5)また、本実施例によれば、前記反射波認識制御において、反射波信号SGECの信号振幅AMSGが前記第1ピーク判定閾値LT1PK、前記谷判定閾値LTBTM、及び前記第2ピーク判定閾値LT2PKと比較されるので、その反射波信号SGECにおける前記第1ピークPK1、前記谷BTM、および前記第2ピークPK2の検出が容易である。なお、上記反射波信号SGECの信号振幅AMSGに替えてその信号振幅AMSGの輝度変調後の輝度の高さが、上記第1ピーク判定閾値LT1PK、上記谷判定閾値LTBTM、及び上記第2ピーク判定閾値LT2PKと比較されても同様である。
 (A6)また、本実施例によれば、図2に示すように、超音波プローブ24は、複数個の超音波発振子が直線的に血管20の長手方向とは直交して配列された互いに平行な一対の第1短軸超音波アレイ探触子24aおよび第2短軸超音波アレイ探触子24bと、その第1短軸超音波アレイ探触子24aおよびその第2短軸超音波アレイ探触子24bの中央部に隣接して設けられ複数個の超音波発振子が直線的に血管20の長手方向に配列された長軸用超音波アレイ探触子24cとを一平面に備えている。そして、血管縦断面画像生成手段82は、その長軸用超音波アレイ探触子24cで受信された超音波の反射波信号SGECに基づき前記血管縦断面画像を生成する。従って、実用化されている超音波プローブを用いて、上記血管縦断面画像を生成できる。
 (A7)また、本実施例によれば、血管径測定手段88は、血管20の阻血開放前に血管20の安静径dを予め測定し、血管20の阻血開放後に血管20の最大径dMAXを測定し、その血管20の安静径dに対する阻血開放後の径変化割合の最大値RMAXを算出する。そして、指標値算出手段84は、上記安静径dの測定時における前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBと、上記最大径dMAXの測定時における前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBとに基づいて、血管径測定手段88により算出された血管20の安静径dに対する阻血開放後の最大径変化割合RMAXの信頼度を示す指標値(FMD信頼度指標値)XCRFMDを算出する。従って、血管20の安静径dに対する阻血開放後の径変化割合の最大値RMAXを複数対比する場合に、その最大値RMAX毎に得られた上記FMD信頼度指標値XCRFMDを用いることで、測定精度の悪いデータ(最大値RMAX)を排除することが可能となり、例えば、FMD評価結果の信頼性をより引き上げることが可能である。また、臨床の現場において、FMD信頼度指標値XCRFMDが低い場合には再検査をするなどの客観的な判断をすることが可能となる。また、FMD信頼度指標値XCRFMDは、オペレータがFMD評価において血管内腔径dの計測を練習する場合に、その上達度を客観的に評価できる目安にもなる。
 (A8)また、本実施例によれば、図14に示すように、指標値算出手段84は、前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBの各々を、数値表示するとともに、その指標値XCRF,XCRBの大きさに応じて連続的に変化し相互に対比可能な画像乃至は図形(扇形グラフ)として、モニタ画面表示装置30に表示させる。従って、オペレータは、それらの指標値XCRF,XCRBが数値だけで表示されている場合と比較して、直感的に画像の明瞭度を判断でき前壁BRFおよび後壁BRBの画像をより迅速に明瞭なものとすることが可能となる。
 続いて、本発明の他の実施例を説明する。なお、以下の説明において実施例相互に共通する部分には同一の符号を付して説明を省略する。
 前述の第1実施例は、血管検査装置22のオペレータが生体14の表皮上に配置された超音波プローブ24の位置を手動操作で微調整するために、前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBが算出されてそれらが表示されるものであるが、前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBに基づき多軸駆動装置26を駆動して、前記血管縦断面画像内の内中膜複合体IMCの明瞭度が一定限度以上となるように、超音波プローブ24の位置を自動的に微調整するも可能である。本第2実施例では、その超音波プローブ24の位置を自動的に微調整する制御について説明する。この第2実施例の制御は、前述の第1実施例の表示制御に替えて実行されてもよいが、上記第1実施例の表示制御と同時に実行されてもよい。以下、第1実施例と共通する点の説明は省略し、主として第1実施例と異なる点について説明する。
 図17は、血管検査装置22(血管画像評価部100)に備えられた制御機能の要部を説明する機能ブロック線図であって、第1実施例の図8に相当する図である。本実施例の血管画像評価部100は、第1実施例と同様に血管断面画像生成手段82と指標値算出手段84と血管径測定手段88とを備え、更に、短軸画像位置確定部としての短軸画像位置確定手段210と、超音波プローブ位置修正部としての超音波プローブ位置修正手段212とを備えている。そして、超音波プローブ位置修正手段212は、長軸画像スコア算出部としての長軸画像スコア算出手段214を備えている。
 血管断面画像生成手段82は、第1実施例と同様に、第1短軸用超音波アレイ探触子24aで受信された超音波の反射波信号SGECに基づき第1短軸画像表示領域G1に表示させる血管20の横断面画像すなわち第1血管横断面画像を逐次生成し、第2短軸用超音波アレイ探触子24bで受信された超音波の反射波信号SGECに基づき第2短軸画像表示領域G2に表示させる血管20の横断面画像すなわち第2血管横断面画像を逐次生成し、長軸用超音波アレイ探触子24cで受信された超音波の反射波信号SGECに基づき長軸画像表示領域G3に表示させる血管20の縦断面画像すなわち血管縦断面画像を逐次生成する。
 短軸画像位置確定手段210は、血管断面画像生成手段82により生成された前記第1血管横断面画像および前記第2血管横断面画像のそれぞれにおいて血管20の中心位置CRBV(血管横断面中心CRBV)を認識する。そして、短軸画像位置確定手段210は、第1短軸用超音波アレイ探触子24aから血管20の中心(血管横断面中心CRBV)までの距離と第2短軸用超音波アレイ探触子24bから血管20の中心(血管横断面中心CRBV)までの距離とが互いに等しくなるように、且つ、第1短軸画像表示領域G1および第2短軸画像表示領域G2の何れでもそれらの幅方向中央部に血管20の画像を位置させるように、多軸駆動装置26を作動させて超音波プローブ24を位置決めする。図4を用いて説明すれば、その図4において「a=b,c=d,e=f」となるように、すなわち、超音波プローブ24が前記所定の計測位置PT1に配置されるように、多軸駆動装置26を作動させて超音波プローブ24を位置決めする。例えば、前記第1血管横断面画像および前記第2血管横断面画像のそれぞれでの上記所定の計測位置PT1に対する血管中心位置CRBVのずれ(ずれ量及びずれ方向)と、多軸駆動装置26が有する各アクチュエータの超音波プローブ24を上記所定の計測位置PT1に位置決めするための作動量との関係を予め実験的に求め短軸画像位置確定手段210に設定しておく。短軸画像位置確定手段210は、前記第1血管横断面画像および前記第2血管横断面画像のそれぞれにおいて、血管中心位置CRBVを認識した上でその血管中心位置CRBVの上記所定の計測位置PT1に対するずれ(ずれ量及びずれ方向)を算出する。そして、その算出した血管中心位置CRBVのずれが零または略零とみなせる所定の許容範囲内であるか否かを判断し、その血管中心位置CRBVのずれが零または上記所定の許容範囲内ではないと判断した場合には、前記予め設定された血管中心位置CRBVのずれと多軸駆動装置26が有する各アクチュエータの作動量との関係から、上記算出した血管中心位置CRBVのずれに基づき上記各アクチュエータの作動量を決定し多軸駆動装置26を作動させる。短軸画像位置確定手段210は、このような血管中心位置CRBVの認識、血管中心位置CRBVの上記所定の計測位置PT1に対するずれの算出、及び、多軸駆動装置26の作動を、前記血管中心位置CRBVのずれが零または前記所定の許容範囲内となるまで、例えば血管断面画像生成手段82が第1血管横断面画像及び第2血管横断面画像を生成する毎に繰り返し行う。そして、その血管中心位置CRBVのずれが零または前記所定の許容範囲内となったと判断した場合には、多軸駆動装置26の作動を終了し、血管横断面画像に基づく超音波プローブ24の位置決めが完了した旨を超音波プローブ位置修正手段212に出力する。
 長軸画像スコア算出手段214は、指標値算出手段84により算出された前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBに基づいて、血管縦断面画像内での内中膜複合体IMCを表す画像の明瞭度を前壁BRF及び後壁BRBの全体について示す算出値BRFBすなわち総合明瞭度指標値BRFBを算出する。例えば、その総合明瞭度指標値BRFBは、前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBの平均値或いは合計値であるが、本実施例ではそれら指標値XCRF,XCRBの平均値すなわち指標平均値とされる。また、長軸画像スコア算出手段214は、例えば、指標値算出手段84が前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBを算出する毎に、換言すれば、血管断面画像生成手段82が血管縦断面画像を生成するための超音波の反射波信号SGECが受信(取得)される毎に、上記総合明瞭度指標値BRFBを逐次算出する。
 超音波プローブ位置修正手段212は、短軸画像位置確定手段210による超音波プローブ24の位置決め完了後に、前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBに基づく算出値(総合明瞭度指標値BRFB)が予め定められた判定値(指標平均判定値)LBR1FB以上になるように、多軸駆動装置26を作動させて超音波プローブ24を位置決めする、すなわち、超音波プローブ24の位置を微調整する。上記指標平均判定値LBR1FBは、例えば、血管径を測定可能な程度に血管縦断面画像内の内中膜複合体IMCが明瞭となる値に実験的に設定されている。具体的に言えば、超音波プローブ位置修正手段212は、長軸画像スコア算出手段214により算出された総合明瞭度指標値BRFBが前記指標平均判定値LBR1FB以上であるか否かを判断し、その総合明瞭度指標値BRFBがその指標平均判定値LBR1FB以上ではないと判断した場合には、その総合明瞭度指標値BRFBを増加させる方向に多軸駆動装置26を作動させて超音波プローブ24の位置を所定の微小量だけずらす。そして、超音波プローブ位置修正手段212は、この超音波プローブ24の位置を所定の微小量だけずらすことを、総合明瞭度指標値BRFBが前記指標平均判定値LBR1FB以上となるまで繰り返し実行し、総合明瞭度指標値BRFBが前記指標平均判定値LBR1FB以上になったと判断した場合には、多軸駆動装置26の作動を終了し、超音波プローブ24の位置の微調整が完了した旨を血管径測定手段88に出力する。なお、超音波プローブ位置修正手段212は、超音波プローブ24の位置の微調整をする場合には、例えば、多軸駆動装置26を作動させて超音波プローブ24の位置を上記所定の微小量だけ移動させ、その次のサイクルで総合明瞭度指標値BRFBがその移動前に対して減少した場合は、超音波プローブ24の位置をその移動前に戻し更に反対側に上記所定の微小量だけ移動させる。また、前記指標平均判定値LBR1FB以上の範囲が、本発明の予め定められた目標範囲に対応する。
 血管径測定手段88は、第1実施例で説明した機能に加え、超音波プローブ位置修正手段212から超音波プローブ24の位置の微調整が完了した旨を受けた場合に、FMD計測、具体的には、血管内腔径d(d,dMAX)の計測を開始する。
 図18は、本実施例の血管検査装置22(電子制御装置28)の制御作動の要部、すなわち、超音波プローブ24をFMD計測可能な位置に自動的に位置決めする制御作動を説明するためのフローチャートであって、第1実施例の図15に相当する図である。なお、図18のSC6、SC7、SC8はそれぞれ、図15のSA2、SA3、SA4と同じであるので、その説明を省略する。この図18に示す制御作動は、単独で或いは他の制御作動と並列的に実行される。
 血管断面画像生成手段82に対応するSC1においては、Bモード短軸画像データが取得される。すなわち、第1短軸用超音波アレイ探触子24aにおいてその長手方向に所定の反射波受信間隔PCRVで超音波の反射波信号SGECが走査され受信(取得)され、第2短軸用超音波アレイ探触子24bにおいてその長手方向に所定の反射波受信間隔PCRVで超音波の反射波信号SGECが走査され受信(取得)される。その取得された上記Bモード短軸画像データ(反射波信号SGEC)は電子制御装置28の記憶装置等に記憶される。そして、第1短軸用超音波アレイ探触子24aで受信された超音波の反射波信号SGECに基づき前記第1血管横断面画像が生成され、第2短軸用超音波アレイ探触子24bで受信された超音波の反射波信号SGECに基づき前記第2血管横断面画像が生成される。
 短軸画像位置確定手段210に対応するSC2においては、モニタ画面表示装置30の左右の短軸画像すなわち前記第1血管横断面画像および前記第2血管横断面画像のそれぞれにおいて、血管20の中心位置CRBVが認識される。
 短軸画像位置確定手段210に対応するSC3においては、前記第1血管横断面画像および前記第2血管横断面画像のそれぞれにおいて、血管中心位置CRBVの前記所定の計測位置PT1に対するずれ(ずれ量及びずれ方向)が算出される。
 短軸画像位置確定手段210に対応するSC4においては、上記SC3にて算出された血管中心位置CRBVのずれが零または略零であるか否かが判断される。このSC4の判断が肯定された場合、すなわち、上記血管中心位置CRBVのずれが零または略零である場合には、SC6に移る。一方、このSC4の判断が否定された場合には、SC5に移る。
 短軸画像位置確定手段210に対応するSC5においては、多軸駆動装置26が有する各アクチュエータすなわちプローブホルダモーターが上記血管中心位置CRBVのずれを減少させるように制御される。例えば、予め実験的に設定された関係から、上記血管中心位置CRBVのずれに基づきそのずれを減少させる方向に多軸駆動装置26が作動させられる。SC5の次はSC1に戻る。
 SC8の次のSC9においては、SC7で算出された前壁画像明瞭度指標値(前壁スコア)XCRFおよびSC8で算出された後壁画像明瞭度指標値(後壁スコア)XCRBの平均値が算出され、その平均値が、前記総合明瞭度指標値BRFBすなわち前記血管縦断面画像(長軸画像)のスコアとされる。このSC9は、長軸画像スコア算出手段214に対応する。
 超音波プローブ位置修正手段212に対応するSC10においては、上記SC9で算出された総合明瞭度指標値BRFBが前記指標平均判定値(閾値)LBR1FB以上であるか否かが判断される。このSC10の判断が肯定された場合、すなわち、上記総合明瞭度指標値BRFBが上記指標平均判定値LBR1FB以上である場合には、SC12に移る。一方、このSC10の判断が否定された場合には、SC11に移る。
 超音波プローブ位置修正手段212に対応するSC11においては、多軸駆動装置26が有する各アクチュエータ(プローブホルダモーター)が上記総合明瞭度指標値BRFBを増加させる方向に制御される。例えば、超音波プローブ24の位置が多軸駆動装置26の作動により前記所定の微小量だけ移動させられ、その次のサイクルで総合明瞭度指標値BRFBがその移動前に対して減少した場合は、超音波プローブ24の位置がその移動前に戻され更に反対側に上記所定の微小量だけ移動させられる。SC11の次はSC6に戻る。
 血管径測定手段88に対応するSC12においては、FMD計測、具体的には、血管内腔径d(d,dMAX)の計測が開始される。
 本実施例には、前述の第1実施例の効果に加え、更に次のような効果がある。本実施例によれば、短軸画像位置確定手段210は、第1短軸用超音波アレイ探触子24aから血管20の中心(血管横断面中心CRBV)までの距離と第2短軸用超音波アレイ探触子24bから血管20の中心(血管横断面中心CRBV)までの距離とが互いに等しくなるように、且つ、モニタ画面表示装置30の第1短軸画像表示領域G1および第2短軸画像表示領域G2の何れでもそれらの幅方向中央部に血管20の画像を位置させるように、多軸駆動装置26を作動させて超音波プローブ24を位置決めする。そして、超音波プローブ位置修正手段212は、短軸画像位置確定手段210による超音波プローブ24の位置決め完了後に、前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBに基づく算出値(総合明瞭度指標値BRFB)が予め定められた指標平均判定値LBR1FB以上になるように、多軸駆動装置26を作動させて超音波プローブ24を位置決めする。従って、オペレータの操作負担を軽減することができる。また、オペレータが前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBを増加させるように手動操作で超音波プローブ24の位置を微調整する場合と比較して、そのオペレータの熟練の程度が更に低くても明瞭な血管縦断面画像を得ることが可能である。また、超音波プローブ24の位置の微調整による血管縦断面画像の明瞭化までをも含めたFMD計測の自動化を実現できる。
 以上、本発明の実施例を図面に基づいて詳細に説明したが、これはあくまでも一実施形態であり、本発明は当業者の知識に基づいて種々の変更、改良を加えた態様で実施することができる。
 例えば、前述の実施例において、前記前壁画像明瞭度指標値XCRFは前記式(1)により算出され、前記後壁画像明瞭度指標値XCRBは前記式(2)により算出されるが、それら前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBはその他の方法により算出されることも考え得る。
 また、前述の実施例の図14において、前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBの各々は前記扇形グラフでモニタ画面表示装置30に表示されるが、図19に例示するような棒グラフなどその他の表示方法が採用されても差し支えない。また、前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBの各々は数値表示されるとともに図形(扇形グラフ)としてもモニタ画面表示装置30に表示されるが、数値表示と図形表示の何れか一方のみの表示方法であっても差し支えない。
 また、前述の実施例の図13において、前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBの各々は、前記血管縦断面画像および前記血管横断面画像を表示するモニタ画面表示装置30に表示されているが、そのモニタ画面表示装置30とは別の独立した表示装置に表示されても差し支えない。
 また、前述の実施例において、前記FMD信頼度指標値XCRFMDは、前記安静径dの測定時における前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBの平均値と、前記最大径dMAXの測定時における前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBの平均値とを更に前記式(3)により平均して算出されるが、そのようして算出されることに限定されるものではない。例えば、それらの指標値XCRF,XCRBの全ての合計値であっても差し支えない。
 また、前述の実施例において、反射波認識制御手段86は、図10に示すような反射波信号SGECに対し前記反射波認識制御を実行するが、その実行に先立って、その反射波信号SGECによく知られた先鋭化(鮮鋭化)フィルタを適用し、その反射波信号SGECのピークと谷とをより際立たせ明瞭なものとする先鋭化処理を行ってもよい。その先鋭化処理には、例えば、よく知られたアンシャープマスキング法や選択的画像鮮鋭化法などが用いられる。そのアンシャープマスキング法とは、ボケた画像からその画像の2次微分画像(ラプラシアン画像)を引くことにより画像の鮮鋭化を行う方法である。また、選択的画像鮮鋭化法とは、画像のエッジ部分に対するラプラシアン画像のみを抽出することにより、ノイズの影響を抑えて画像のエッジ部分を選択的に鮮鋭化する方法である。
 また、前述の実施例においては、超音波プローブ24の位置決め後にFMD計測が実施される例が説明されているが、本発明は、IMT(Intima-media thickness)検査や頚動脈を超音波で観測するような制御装置においても有効に適用できる。
 また、前述の実施例において、前壁画像明瞭度指標値XCRFおよび後壁画像明瞭度指標値XCRBは、前記反射波認識制御の実行結果に基づいて算出されるが、前記血管縦断面画像内の内中膜複合体IMCの明瞭度と関連付けられて算出されるのであれば、他の方法により算出されても差し支えない。
 また、前述の実施例において、電子制御装置28は血管径測定手段88を備えているが、これを備えていない電子制御装置28も考え得る。
 また、前述の実施例において、図1の上腕16は、例えば人体の上腕である。
 また前述した複数の実施例はそれぞれ、例えば優先順位を設けるなどして、相互に組み合わせて実施することができる。
 その他、一々例示はしないが、本発明はその趣旨を逸脱しない範囲内において種々の変更が加えられて実施されるものである。
14:生体
20:血管
22:血管検査装置(超音波血管検査装置)
24:超音波プローブ
24a:第1短軸用超音波アレイ探触子
24b:第2短軸用超音波アレイ探触子
24c:長軸用超音波アレイ探触子
26:多軸駆動装置
30:モニタ画面表示装置(画像表示装置)
82:血管断面画像生成手段(血管縦断面画像生成手段)
84:指標値算出手段
86:反射波認識制御手段
88:血管径測定手段
210:短軸画像位置確定手段
212:超音波プローブ位置修正手段
IMC:内中膜複合体
BRF:前壁
BRB:後壁
SGEC:反射波信号

Claims (9)

  1.  生体の表皮上に配置された超音波プローブを用いて超音波の反射波信号に基づき該生体の表皮下の血管の縦断面画像を生成する血管縦断面画像生成手段を備えた超音波血管検査装置であって、
     前記血管の縦断面画像内で該血管の内中膜複合体を表す画像の明瞭度を示す指標値を算出する指標値算出手段
     を含むことを特徴とする超音波血管検査装置。
  2.  前記指標値算出手段は、前記血管の縦断面画像内の一対の血管壁断面のうち前記超音波プローブに近い側の血管壁断面である前壁における前記内中膜複合体を表す画像の明瞭度を示す前壁画像明瞭度指標値と、前記一対の血管壁断面のうち前記超音波プローブから遠い側の血管壁断面である後壁における前記内中膜複合体を表す画像の明瞭度を示す後壁画像明瞭度指標値とを算出する
     ことを特徴とする請求項1に記載の超音波血管検査装置。
  3.  前記超音波プローブの前記血管の長手方向の相互に異なる複数の受信位置で受信された複数の前記反射波信号の各々に対し、該反射波信号の大きさと前記血管の径方向位置との関係において、該反射波信号の大きさが予め定められた第1ピーク判定閾値を超える第1ピークと、該第1ピークの発生位置よりも前記血管の径方向外側に発生し前記反射波信号の大きさが予め定められた谷判定閾値よりも小さい谷と、前記第1ピークの発生位置から前記谷を介して前記血管の径方向外側の予め定められたピーク間隔閾値を超えない範囲内で発生し前記反射波信号の大きさが予め定められた第2ピーク判定閾値を超える第2ピークとを検出する反射波認識制御を、前記前壁および前記後壁のそれぞれで実行する反射波認識制御手段が設けられており、
     前記指標値算出手段は、前記反射波認識制御手段が前記前壁において前記反射波認識制御の実行により前記複数の反射波信号のうち前記第1ピークと前記谷と前記第2ピークとの全部を検出した反射波信号の数に基づいて前記前壁画像明瞭度指標値を算出し、前記反射波認識制御手段が前記後壁において前記反射波認識制御の実行により前記複数の反射波信号のうち前記第1ピークと前記谷と前記第2ピークとの全部を検出した反射波信号の数に基づいて前記後壁画像明瞭度指標値を算出する
     ことを特徴とする請求項2に記載の超音波血管検査装置。
  4.  前記反射波認識制御手段は、予め定められた前記血管の長手方向の観測対象範囲内で受信された前記反射波信号に対して、前記反射波認識制御を実行する
     ことを特徴とする請求項3に記載の超音波血管検査装置。
  5.  前記反射波信号の大きさとは、該反射波信号の振幅または該反射波信号の振幅を前記血管の縦断面画像を表示するための輝度に置き換えたときの該輝度の高さである
     ことを特徴とする請求項3又は4に記載の超音波血管検査装置。
  6.  前記超音波プローブは、複数個の超音波発振子が直線的に前記血管の長手方向とは直交して配列された互いに平行な一対の第1短軸超音波アレイ探触子および第2短軸超音波アレイ探触子と、該第1短軸超音波アレイ探触子および該第2短軸超音波アレイ探触子の一方又は両方の中央部に隣接して設けられ複数個の超音波発振子が直線的に前記血管の長手方向に配列された長軸用超音波アレイ探触子とを一平面に備えており、
     前記血管縦断面画像生成手段は、前記長軸用超音波アレイ探触子で受信された前記超音波の反射波信号に基づき前記血管の縦断面画像を生成する
     ことを特徴とする請求項2乃至5の何れか1項に記載の超音波血管検査装置。
  7.  前記第1短軸用超音波アレイ探触子による超音波画像を表示する第1短軸画像表示領域と、前記第2短軸用超音波アレイ探触子による超音波画像を表示する第2短軸画像表示領域と、前記血管の縦断面画像を表示する長軸画像表示領域とを有する画像表示装置と、
     前記超音波プローブの位置決め状態を制御する多軸駆動装置とが設けられており、
     前記第1短軸用超音波アレイ探触子から前記血管の中心までの距離と前記第2短軸用超音波アレイ探触子から前記血管の中心までの距離とが互いに等しくなるように、且つ、前記第1短軸画像表示領域および前記第2短軸画像表示領域の何れでもそれらの幅方向中央部に前記血管の画像を位置させるように、前記多軸駆動装置を作動させて前記超音波プローブを位置決めする短軸画像位置確定手段と、
     前記短軸画像位置確定手段による前記超音波プローブの位置決め完了後に、前記前壁画像明瞭度指標値および前記後壁画像明瞭度指標値に基づく算出値が予め定められた目標範囲内になるように、前記多軸駆動装置を作動させて前記超音波プローブを位置決めする超音波プローブ位置修正手段と
     を、含むことを特徴とする請求項6に記載の超音波血管検査装置。
  8.  前記血管の阻血開放前に該血管の安静径を予め測定し、阻血開放後に前記血管の最大径を測定し、前記血管の安静径に対する阻血開放後の径変化割合の最大値を算出する血管径測定手段が設けられており、
     前記指標値算出手段は、前記安静径の測定時における前記前壁画像明瞭度指標値および前記後壁画像明瞭度指標値と、前記最大径の測定時における前記前壁画像明瞭度指標値および前記後壁画像明瞭度指標値とに基づいて、前記血管径測定手段により算出された前記血管の安静径に対する阻血開放後の径変化割合の最大値の信頼度を示す指標値を算出する
     ことを特徴とする請求項2乃至7の何れか1項に記載の超音波血管検査装置。
  9.  前記指標値算出手段は、前記前壁画像明瞭度指標値および前記後壁画像明瞭度指標値の各々を、その指標値の大きさに応じて連続的に変化し相互に対比可能な画像として画像表示装置に表示させる
     ことを特徴とする請求項2乃至8の何れか1項に記載の超音波血管検査装置。
PCT/JP2009/062909 2009-07-16 2009-07-16 超音波血管検査装置 WO2011007439A1 (ja)

Priority Applications (6)

Application Number Priority Date Filing Date Title
US13/383,866 US8939908B2 (en) 2009-07-16 2009-07-16 Ultrasonic blood vessel inspecting apparatus
JP2011522663A JP5411272B2 (ja) 2009-07-16 2009-07-16 超音波血管検査装置
EP09847337.4A EP2454995B1 (en) 2009-07-16 2009-07-16 Ultrasonic blood vessel examination apparatus
CN200980160468.4A CN102469988B (zh) 2009-07-16 2009-07-16 超声波血管检查装置
PCT/JP2009/062909 WO2011007439A1 (ja) 2009-07-16 2009-07-16 超音波血管検査装置
KR1020127003580A KR101624846B1 (ko) 2009-07-16 2009-07-16 초음파 혈관 검사장치

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2009/062909 WO2011007439A1 (ja) 2009-07-16 2009-07-16 超音波血管検査装置

Publications (1)

Publication Number Publication Date
WO2011007439A1 true WO2011007439A1 (ja) 2011-01-20

Family

ID=43449055

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2009/062909 WO2011007439A1 (ja) 2009-07-16 2009-07-16 超音波血管検査装置

Country Status (6)

Country Link
US (1) US8939908B2 (ja)
EP (1) EP2454995B1 (ja)
JP (1) JP5411272B2 (ja)
KR (1) KR101624846B1 (ja)
CN (1) CN102469988B (ja)
WO (1) WO2011007439A1 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012249844A (ja) * 2011-06-03 2012-12-20 Fujifilm Corp 超音波診断装置
JP2013017721A (ja) * 2011-07-13 2013-01-31 Seiko Epson Corp 超音波測定装置及び超音波測定装置の制御方法
JP2017205482A (ja) * 2016-05-17 2017-11-24 東芝メディカルシステムズ株式会社 支援装置及び支援方法

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8439839B2 (en) * 2007-06-04 2013-05-14 Panasonic Corporation Ultrasonic diagnosis device and ultrasonic probe for use in ultrasonic diagnosis device
ES2651898T3 (es) 2007-11-26 2018-01-30 C.R. Bard Inc. Sistema integrado para la colocación intravascular de un catéter
US8781555B2 (en) 2007-11-26 2014-07-15 C. R. Bard, Inc. System for placement of a catheter including a signal-generating stylet
US9521961B2 (en) 2007-11-26 2016-12-20 C. R. Bard, Inc. Systems and methods for guiding a medical instrument
US9532724B2 (en) 2009-06-12 2017-01-03 Bard Access Systems, Inc. Apparatus and method for catheter navigation using endovascular energy mapping
EP2912999B1 (en) 2010-05-28 2022-06-29 C. R. Bard, Inc. Apparatus for use with needle insertion guidance system
KR20150000627A (ko) * 2013-06-25 2015-01-05 삼성전자주식회사 초음파 영상 장치 및 그 제어 방법
JP2015139629A (ja) * 2014-01-30 2015-08-03 セイコーエプソン株式会社 超音波測定装置および超音波測定方法
US20190015076A1 (en) * 2015-12-21 2019-01-17 Koninklijke Philips N.V. Ultrasound imaging apparatus and ultrasound imaging method for inspecting a volume of a subject
KR102425141B1 (ko) * 2016-03-09 2022-07-26 에코너스 인코퍼레이티드 인공지능 네트워크를 이용한 초음파 이미지 인식 시스템 및 방법
JP6371334B2 (ja) * 2016-05-27 2018-08-08 株式会社ユネクス 超音波断面画像測定装置
US11020563B2 (en) 2016-07-14 2021-06-01 C. R. Bard, Inc. Automated catheter-to-vessel size comparison tool and related methods
US10928494B1 (en) * 2017-07-10 2021-02-23 Interson Corporation Ultrasound plural display apparatus with light indicator for positional access
WO2020081373A1 (en) 2018-10-16 2020-04-23 Bard Access Systems, Inc. Safety-equipped connection systems and methods thereof for establishing electrical connections
CN109771245B (zh) * 2019-02-18 2021-07-13 中山大学附属第一医院 利用体外反搏技术评价个体血管内皮功能的方法和系统
CA3152545A1 (en) 2019-09-20 2021-03-25 Bard Access Systems, Inc. Automatic vessel detection tools and methods
WO2022020351A1 (en) 2020-07-21 2022-01-27 Bard Access Systems, Inc. System, method and apparatus for magnetic tracking of ultrasound probe and generation of 3d visualization thereof
WO2022051657A1 (en) 2020-09-03 2022-03-10 Bard Access Systems, Inc. Portable ultrasound systems and methods
KR102250219B1 (ko) * 2020-09-15 2021-05-11 주식회사 아이도트 초음파 진단 시스템
US11925505B2 (en) 2020-09-25 2024-03-12 Bard Access Systems, Inc. Minimum catheter length tool
JP2022074392A (ja) * 2020-11-04 2022-05-18 コニカミノルタ株式会社 超音波診断装置、超音波診断装置の制御方法、及び、超音波診断装置の制御プログラム

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005020821A1 (ja) * 2003-09-01 2005-03-10 Matsushita Electric Industrial Co., Ltd. 生体信号モニタ装置
JP2007222291A (ja) * 2006-02-22 2007-09-06 Yunekusu:Kk 動脈血管判定方法および装置
JP2009089911A (ja) 2007-10-09 2009-04-30 Yunekusu:Kk 血管超音波画像測定方法
JP2009153573A (ja) * 2007-12-25 2009-07-16 Panasonic Corp 超音波診断装置および断層画像処理方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2889568B1 (ja) * 1998-05-18 1999-05-10 正男 伊藤 血管膜厚測定装置及び動脈硬化診断装置
US6817982B2 (en) * 2002-04-19 2004-11-16 Sonosite, Inc. Method, apparatus, and product for accurately determining the intima-media thickness of a blood vessel
US6835177B2 (en) * 2002-11-06 2004-12-28 Sonosite, Inc. Ultrasonic blood vessel measurement apparatus and method
US7727153B2 (en) 2003-04-07 2010-06-01 Sonosite, Inc. Ultrasonic blood vessel measurement apparatus and method
US20060184034A1 (en) * 2005-01-27 2006-08-17 Ronen Haim Ultrasonic probe with an integrated display, tracking and pointing devices
US7955265B2 (en) * 2005-08-15 2011-06-07 General Electric Company Method and apparatus for measuring anatomic structures
WO2007108359A1 (ja) * 2006-03-20 2007-09-27 Matsushita Electric Industrial Co., Ltd. 超音波診断装置
JP4799276B2 (ja) 2006-05-30 2011-10-26 パナソニック株式会社 超音波診断装置
JP4918369B2 (ja) * 2007-01-16 2012-04-18 パナソニック株式会社 超音波診断装置
JP5346440B2 (ja) * 2007-02-15 2013-11-20 富士フイルム株式会社 超音波診断装置及びデータ計測プログラム
JPWO2009013871A1 (ja) * 2007-07-20 2010-09-30 パナソニック株式会社 超音波診断装置
JP5412242B2 (ja) * 2009-11-05 2014-02-12 伸治 久米 超音波断層画像の処理装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005020821A1 (ja) * 2003-09-01 2005-03-10 Matsushita Electric Industrial Co., Ltd. 生体信号モニタ装置
JP2007222291A (ja) * 2006-02-22 2007-09-06 Yunekusu:Kk 動脈血管判定方法および装置
JP2009089911A (ja) 2007-10-09 2009-04-30 Yunekusu:Kk 血管超音波画像測定方法
JP2009153573A (ja) * 2007-12-25 2009-07-16 Panasonic Corp 超音波診断装置および断層画像処理方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012249844A (ja) * 2011-06-03 2012-12-20 Fujifilm Corp 超音波診断装置
JP2013017721A (ja) * 2011-07-13 2013-01-31 Seiko Epson Corp 超音波測定装置及び超音波測定装置の制御方法
JP2017205482A (ja) * 2016-05-17 2017-11-24 東芝メディカルシステムズ株式会社 支援装置及び支援方法
JP7019301B2 (ja) 2016-05-17 2022-02-15 キヤノンメディカルシステムズ株式会社 支援装置

Also Published As

Publication number Publication date
KR101624846B1 (ko) 2016-05-27
EP2454995B1 (en) 2019-01-23
JP5411272B2 (ja) 2014-02-12
EP2454995A1 (en) 2012-05-23
US8939908B2 (en) 2015-01-27
EP2454995A4 (en) 2017-03-15
CN102469988A (zh) 2012-05-23
US20120116227A1 (en) 2012-05-10
KR20120041756A (ko) 2012-05-02
JPWO2011007439A1 (ja) 2012-12-20
CN102469988B (zh) 2014-12-03

Similar Documents

Publication Publication Date Title
JP5411272B2 (ja) 超音波血管検査装置
JP5219228B2 (ja) 血管機能検査装置
EP2198784B1 (en) Blood-vessel ultrasonographing method
JP5224221B2 (ja) 血管機能検査装置
JP5474986B2 (ja) 血管機能検査装置
JP5474966B2 (ja) 生体血管状態測定装置
JP2007195662A (ja) 生体内管状体の内腔径測定装置
JP5998197B2 (ja) 生体の血管径連続測定装置
WO2017204176A1 (ja) 超音波断面画像測定装置
JP5527953B2 (ja) 血管拡張反応画像表示装置
JP6671065B2 (ja) 上腕動脈用血管内皮機能測定装置
JP5896759B2 (ja) 生体の動脈内皮機能測定装置
JP6243719B2 (ja) 生体血管状態測定装置
JP6192490B2 (ja) 生体血管状態測定装置
Zhang et al. Vibration mode imaging
JP6013377B2 (ja) 生体血管パラメータ測定装置
JP2008212548A (ja) 超音波診断装置
JP5723960B2 (ja) 血管機能検査装置
JP2017131374A (ja) 血管内皮機能検査装置

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 200980160468.4

Country of ref document: CN

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

Ref document number: 09847337

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 13383866

Country of ref document: US

Ref document number: 2011522663

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 2009847337

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 875/CHENP/2012

Country of ref document: IN

ENP Entry into the national phase

Ref document number: 20127003580

Country of ref document: KR

Kind code of ref document: A