US11020093B2 - Method and apparatus for improving the measurement of flow velocity of blood - Google Patents

Method and apparatus for improving the measurement of flow velocity of blood Download PDF

Info

Publication number
US11020093B2
US11020093B2 US16/082,305 US201716082305A US11020093B2 US 11020093 B2 US11020093 B2 US 11020093B2 US 201716082305 A US201716082305 A US 201716082305A US 11020093 B2 US11020093 B2 US 11020093B2
Authority
US
United States
Prior art keywords
ultrasound transducer
subject
spatial time
correction
ultrasound
Prior art date
Legal status (The legal status 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 status listed.)
Active, expires
Application number
US16/082,305
Other languages
English (en)
Other versions
US20190090852A1 (en
Inventor
Bart Kroon
Rick BEZEMER
Denny Mathew
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
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 Koninklijke Philips NV filed Critical Koninklijke Philips NV
Assigned to KONINKLIJKE PHILIPS N.V. reassignment KONINKLIJKE PHILIPS N.V. ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: MATHEW, DENNY, BEZEMER, Rick, KROON, BART
Publication of US20190090852A1 publication Critical patent/US20190090852A1/en
Application granted granted Critical
Publication of US11020093B2 publication Critical patent/US11020093B2/en
Active legal-status Critical Current
Adjusted expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/488Diagnostic techniques involving Doppler signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06Measuring blood flow
    • A61B8/065Measuring blood flow to determine blood output from the heart
    • 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/4245Details of probe positioning or probe attachment to the patient involving determining the position of the probe, e.g. with respect to an external reference frame or to the patient
    • 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
    • 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/20ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems

Definitions

  • the invention relates to a method and apparatus for improving the measurement of flow velocity of blood, and in particular relates to improving measurements of the flow velocity of blood using an ultrasound transducer.
  • Cardiac output is a measure of the amount of blood that is pumped out of the heart in one minute. Cardiac output is determined as stroke volume (i.e. the volume of blood pumped per beat of the heart) multiplied by the number of heart beats per minute (i.e. the heart rate). This parameter (particularly beat-to-beat cardiac output) is difficult to measure non-invasively. However, it is possible to measure cardiac output from the flow velocity of the blood (which can be obtained using ultrasound Doppler) and the cross-sectional area of the aorta.
  • GB 2257254 describes a method and apparatus for measuring the cross-sectional area of blood vessels and blood flow velocity (and hence cardiac output) using Doppler ultrasound.
  • An ultrasound probe is directed downwardly via the suprasternal notch and Doppler signals from the moving blood in the aorta are detected.
  • Signal processing means determine the Doppler power at successive ranges and the resulting power curve is correlated against stored curves representing a range of known aorta cross-sectional areas to find the closest fit.
  • the operator has to manually aim the probe at the center of the aorta, and the cross-sectional area of the aorta is estimated from the weight or length of a patient, or estimated according to other techniques, such as those described in GB 2257254.
  • This approach of calculating cardiac output is unreliable due to difficulty in locating the aorta, and interpersonal variations, especially at young and old ages.
  • multiple measurements of the cardiac output from the same patient can be inconsistent since the measurement is sensitive to the positioning of the probe.
  • the accuracy of the conventional solutions very much depend on the skill of the operator to identify the ascending aorta and determine the correct position of the transducer.
  • a method of improving a measurement of the flow velocity of blood in a blood vessel of a subject comprising: using an ultrasound transducer to emit an ultrasound beam to measure flow velocity in a part of a body of a subject; forming a spatial time velocity profile for the part of the body from the measured flow velocity; and analyzing the spatial time velocity profile to determine a correction to the angle of the ultrasound beam with respect to the subject, the correction being based on a difference between the position of a peak in the spatial time velocity profile and the center of the spatial time velocity profile. Therefore, the method determines how the angle of the ultrasound beam should be corrected in order to center a flow velocity peak and thus improve a measurement of flow velocity.
  • the peak in the spatial time velocity profile comprises the maximum flow velocity in the spatial time velocity profile.
  • the method further comprises the step of providing feedback to an operator of the ultrasound transducer based on the determined correction to the angle.
  • the feedback to the operator can comprise an indication of a rotation of the ultrasound transducer with respect to the subject that would position the peak in the spatial time velocity profile in, or nearer to, the center of the spatial time velocity profile.
  • the indication can be a visual, aural and/or tactile indication, or a combination thereof.
  • the method further comprises the steps of determining a first control signal for the ultrasound transducer based on the determined correction to the angle; and adjusting the direction of the ultrasound beam emitted by the ultrasound transducer according to the determined first control signal such that the position of the peak in the spatial time velocity profile is in, or nearer to, the center of the spatial time velocity profile.
  • the method further comprises the step of determining a second control signal for an orientation actuator for the ultrasound transducer based on the determined correction to the angle, wherein the orientation actuator is capable of adjusting the orientation of the ultrasound transducer with respect to the subject.
  • the method can further comprise the step of controlling the orientation actuator using the second control signal to adjust the orientation of the ultrasound transducer with respect to the subject such that the position of the peak in the spatial time velocity profile would be in, or nearer to, the center of the spatial time velocity profile.
  • the method further comprises the step of determining a correction to a value of a flow velocity in the spatial time velocity profile based on the determined correction to the angle.
  • the step of using an ultrasound transducer comprises using the ultrasound transducer to measure flow velocity at a plurality of depths in the part of the body; the step of forming a spatial time velocity profile comprises forming a respective spatial time velocity profile for each of the plurality of depths, each spatial time velocity profile being formed from the measured flow velocity at the respective depth; and the step of analyzing the spatial time velocity profile to determine a correction to the angle of the ultrasound beam comprises analyzing one or more of the plurality of spatial time velocity profiles to determine the correction to the angle of the ultrasound beam.
  • the method further comprises the step of analyzing the spatial time velocity profiles to determine a correction to the position of the ultrasound transducer with respect to the subject,
  • the step of analyzing the spatial time velocity profiles to determine a correction to the position of the ultrasound transducer with respect to the subject comprises determining the correction based on the fitting of a geometric model of the arterial path through peaks in the spatial time velocity profiles.
  • the method further comprises the step of providing feedback to an operator of the ultrasound transducer based on the determined correction to the position.
  • the feedback to the operator can comprise an indication of a movement of the ultrasound transducer with respect to the subject that would improve the alignment of the ultrasound transducer with the artery.
  • the indication can be a visual, aural and/or tactile indication, or a combination thereof.
  • the method further comprises the step of determining a third control signal for a position actuator for the ultrasound transducer based on the determined correction to the position, wherein the position actuator is capable of adjusting the position of the ultrasound transducer with respect to the subject.
  • the method can further comprise the step of controlling the position actuator using the third control signal to adjust the position of the ultrasound transducer with respect to the subject such that the alignment of the ultrasound transducer with the artery is improved.
  • a computer program product comprising a computer readable medium having computer readable code embodied therein, the computer readable code being configured such that, on execution by a suitable computer, control unit or processor, the computer, control unit or processor is caused to control an ultrasound transducer connected to the computer, control unit or processor to emit an ultrasound beam to measure flow velocity in a part of a body of a subject; form a spatial time velocity profile for the part of the body from the measured flow velocity; and analyze the spatial time velocity profile to determine a correction to the angle of the ultrasound beam with respect to the subject, the correction being based on a difference between the position of a peak in the spatial time velocity profile and the center of the spatial time velocity profile.
  • the computer readable code according to the second aspect is configured such that a suitable computer, control unit or processor is caused to perform any of the methods described above.
  • an apparatus for measuring flow velocity of blood in a blood vessel of a subject comprising a control unit that is configured to obtain measurements of flow velocity in a part of a body of a subject, wherein the measurements are made using an ultrasound transducer that emits an ultrasound beam; form a spatial time velocity profile for the part of the body from the measured flow velocity; and analyze the spatial time velocity profile to determine a correction to the angle of the ultrasound beam with respect to the subject, the correction being based on a difference between the position of a peak in the spatial time velocity profile and the center of the spatial time velocity profile. Therefore, the apparatus determines how the angle of the ultrasound beam should be corrected in order to center a flow velocity peak and thus improve a measurement of flow velocity.
  • the apparatus further comprises an ultrasound transducer.
  • the peak in the spatial time velocity profile comprises the maximum flow velocity in the spatial time velocity profile.
  • control unit is further configured to determine feedback for an operator of the ultrasound transducer based on the determined correction to the angle.
  • the feedback to the operator can comprise an indication of a rotation of the ultrasound transducer with respect to the subject that would position the peak in the spatial time velocity profile in, or nearer to, the center of the spatial time velocity profile.
  • the indication can be a visual, aural and/or tactile indication, or a combination thereof.
  • the apparatus can further comprise a user interface component for providing the determined feedback to the operator. This embodiment has an advantage that the operator can be guided as to an improved or correct positioning of the ultrasound transducer with respect to the subject.
  • control unit is further configured to determine a first control signal for the ultrasound transducer based on the determined correction to the angle; and output the first control signal for adjusting the direction of the ultrasound beam emitted by the ultrasound transducer such that the position of the peak in the spatial time velocity profile is in, or nearer to, the center of the spatial time velocity profile.
  • control unit is further configured to determine a second control signal for an orientation actuator for the ultrasound transducer based on the determined correction to the angle, wherein the orientation actuator is capable of adjusting the orientation of the ultrasound transducer with respect to the subject.
  • control unit is further configured to control the orientation actuator using the second control signal to adjust the orientation of the ultrasound transducer with respect to the subject such that the position of the peak in the spatial time velocity profile would be in, or nearer to, the center of the spatial time velocity profile.
  • the apparatus further comprises an orientation actuator that is configured to adjust the orientation of the ultrasound transducer with respect to the subject. This embodiment has an advantage that the angle of the ultrasound beam can be automatically adjusted to improve the flow velocity measurement.
  • control unit is further configured to determine a correction to a value of a flow velocity in the spatial time velocity profile based on the determined correction to the angle.
  • the measurements of the flow velocity are indicative of flow velocity at a plurality of depths in the part of the body; and wherein the control unit is configured to form a respective spatial time velocity profile for each of the plurality of depths, each spatial time velocity profile being formed from the measured flow velocity at the respective depth; and the control unit is configured to analyze the spatial time velocity profile to determine a correction to the angle of the ultrasound beam by analyzing one or more of the plurality of spatial time velocity profiles to determine a correction to the angle of the ultrasound beam.
  • control unit is further configured to analyze the spatial time velocity profiles to determine a correction to the position of the ultrasound transducer with respect to the subject,
  • control unit is configured to analyze the spatial time velocity profiles to determine a correction to the position of the ultrasound transducer with respect to the subject by determining the correction based on the fitting of a geometric model of the arterial path through peaks in the spatial time velocity profiles.
  • control unit is further configured to determine feedback for an operator of the ultrasound transducer based on the determined correction to the position.
  • the feedback to the operator can comprise an indication of a movement of the ultrasound transducer with respect to the subject that would improve the alignment of the ultrasound transducer with the artery.
  • the indication can be a visual, aural and/or tactile indication.
  • the apparatus can further comprise a user interface component for providing the determined feedback to the operator. This embodiment has an advantage that the operator can be guided as to an improved or correct positioning of the ultrasound transducer with respect to the subject.
  • control unit is further configured to determine a third control signal for a position actuator for the ultrasound transducer based on the determined correction to the position, wherein the position actuator is capable of adjusting the position of the ultrasound transducer with respect to the subject.
  • control unit is further configured to control the position actuator using the third control signal to adjust the position of the ultrasound transducer with respect to the subject such that the alignment of the ultrasound transducer with the artery is improved.
  • the apparatus further comprises a position actuator that is capable of adjusting the position of the ultrasound transducer with respect to the subject. This embodiment has an advantage that the position of the ultrasound transducer can be automatically adjusted to improve the flow velocity measurement.
  • FIG. 1 is an illustration of the type of measurement performed using the invention
  • FIG. 2 is a diagram illustrating an exemplary spatial time velocity profile for an artery
  • FIG. 3 is a block diagram of an apparatus according to an embodiment of the invention.
  • FIG. 4 illustrates a scenario in which the invention can be used
  • FIG. 5 is a flow chart illustrating a method of improving measurements of flow velocity of blood according to an embodiment
  • FIG. 6 illustrates another scenario in which the invention can be used.
  • the invention is for use in measuring the flow velocity of blood in a blood vessel, for example an artery or vein, particularly the ascending aorta or pulmonary artery, of a subject.
  • aorta is the main artery in the human body and starts at the left cardiac ventricle and extends down the chest and abdomen of the body.
  • the ascending aorta is the part of the aorta extending out of the top of the left ventricle before the aorta curves downwards.
  • the pulmonary artery is the artery that connects the right ventricle of the heart to the lungs.
  • an ultrasound transducer can be used on, in or near the suprasternal notch of a subject in order to measure the flow velocity of blood in the ascending aorta. In other embodiments, an ultrasound transducer can be used in a parasternal position to measure the flow velocity of the blood in the pulmonary artery.
  • FIG. 1 illustrates a measurement scenario in which the invention can be used.
  • FIG. 1 shows an aorta 2 (ascending or pulmonary) which is oriented in a generally vertical direction with blood flowing upwards, and an ultrasound transducer 4 that is positioned above the aorta 2 .
  • aorta 2 ascending or pulmonary
  • ultrasound transducer 4 that is positioned above the aorta 2 .
  • transverse, lateral and axial directions relative to the transducer 4 are shown in FIG.
  • the axial direction also referred to as the z-axis
  • the transverse and lateral directions also referred to the x- and y-axes comprising orthogonal axes in the plane of the transducer 4 .
  • the transducer 4 measures flow velocity in the axial direction across at least one cross-section using ultrasound Doppler techniques to form a spatial time velocity profile.
  • FIG. 1 shows two exemplary cross-sections at a distance of 100 mm and 150 mm from the transducer 4 although these should not be considered to be limiting. At a depth of around 150 mm into the body the aorta 4 is generally aligned with the transducer 4 .
  • FIG. 2 An exemplary arterial blood spatial time velocity profile (i.e. a spatial distribution of velocities) is shown in FIG. 2 , and it can be seen that, as expected for a flow of viscous fluid through a tube, the velocity increases towards the center of the artery.
  • a spatial time velocity profile comprises a velocity measurement or estimate at each of a plurality of positions on a two-dimensional plane/cross-section.
  • the invention aims to improve the accuracy of an ultrasound measurement of flow velocity by analyzing a spatial time velocity profile obtained by the ultrasound transducer and determining a correction to the angle of the ultrasound transducer with respect to the subject (which is also referred to herein as the orientation of the transducer).
  • spatial time velocity profiles for a plurality of depths can be analyzed in order to determine a correction to the angle and position of the ultrasound transducer.
  • the determined correction can then be indicated to an operator in order for the operator to make an adjustment to the positioning of the ultrasound transducer, and/or an actuator can be controlled according to the determined correction to automatically adjust the positioning of the transducer.
  • FIG. 3 An exemplary apparatus for measuring the flow velocity of blood according to an embodiment of the invention is shown in FIG. 3 .
  • the apparatus 10 comprises an ultrasound transducer 12 that is for emitting ultrasound and measuring ultrasound reflections.
  • the ultrasound transducer 12 can comprise a plurality of transducer elements that generate ultrasound and measure the reflected ultrasound.
  • the ultrasound transducer 12 can be a continuous-wave (CW) Doppler ultrasound transducer 12 or a pulsed-wave (PW) Doppler ultrasound transducer 12 .
  • CW continuous-wave
  • PW pulsed-wave
  • the ultrasound transducer 12 can comprise a plurality of transducer elements arranged as a linear array of transducer elements, and the ultrasound transducer 12 can be provided with an acoustic lens or other acoustic beam shaping element in order to create an annulus-shaped or toroidal-shaped ultrasound beam.
  • This ultrasound beam can be used to obtain spatial-time-velocity information from the subject (e.g. time-velocity information in the lateral and transverse directions).
  • the ultrasound transducer 12 can comprise a plurality of transducer elements, and the transducer elements can be arranged in a two-dimensional (2D) array.
  • the transducer elements can be driven individually or in combination in order to generate and steer an ultrasound beam and thereby obtain spatial-time-velocity information from the subject.
  • This type of ultrasound transducer 12 can be used to obtain spatial-time-velocity information in a three-dimensional (3D) volume of the subject, and thus spatial-time-velocity information at a plurality of depths in the subject (e.g. time-velocity information in the lateral and transverse directions at a plurality of depths).
  • the ultrasound transducer 12 is connected to a control unit 14 that controls the operation of the transducer 12 and that receives the output of the transducer 12 .
  • the ultrasound transducer 12 can output measurements of the received ultrasound reflections to the control unit 14 , and the control unit 14 can analyze the Doppler shift in the ultrasound reflections to determine measurements of flow velocity.
  • the control unit 14 can comprise one or more processors, processing units, multi-core processors or modules that are configured or programmed to control the apparatus 10 and components thereof to determine a measurement of flow velocity of blood in an artery of the subject as described below.
  • the apparatus 10 can also comprise a memory unit 16 that can be used for storing program code that can be executed by the control unit 14 to perform the method described herein.
  • the memory unit 16 can also be used to store signals and measurements made or obtained by the apparatus 10 and/or by the ultrasound transducer 12 .
  • the apparatus 10 can further comprise a user interface 18 for providing feedback to the operator of the apparatus 10 and/or one or more actuators 20 for adjusting the angle (orientation) and/or position of the ultrasound transducer 12 with respect to the subject based on a control signal from the control unit 14 .
  • the user interface 18 is to provide feedback to the operator to instruct, enable or encourage the operator to adjust the angle and/or position of the ultrasound transducer 12 in order to improve the accuracy and/or reliability of the measurement of flow velocity.
  • the user interface 18 can therefore comprise any component that is suitable for providing this feedback/information, and can be, for example, any one or more of a display screen or other visual indicator, a speaker, one or more lights, and a component for providing tactile feedback (e.g. a vibration function).
  • the user interface 18 can be used to provide the subject or operator of the apparatus 10 (for example a healthcare provider) with information resulting from the method according to the invention.
  • the user interface 18 can provide an indication of the flow velocity of the subject.
  • the one or more actuators 20 can comprise any suitable means for enabling the angle and/or position of the ultrasound transducer 12 to be adjusted with respect to the subject.
  • the one or more actuators 20 can be any type of actuator 20 that is capable of moving or controlling a mechanism or system to change the orientation and/or position of the ultrasound transducer 12 .
  • the actuator can be based on electric, magnetic, pneumatic or hydraulic principles, or any combination thereof.
  • the most suitable types of actuator are electric or magnetic actuators.
  • An electric actuator can be a motor or piezoelectric component.
  • a magnetic actuator can be actuated by applying magnetic energy. These types of actuators are typically compact, lightweight, economical, and with high power density.
  • Magnetic actuators can use shape memory materials (SMMs), such as shape memory alloys (SMAs), or magnetic shape-memory alloys (MSMAs).
  • SMMs shape memory materials
  • SMAs shape memory alloys
  • MSMAs magnetic shape-memory alloys
  • EAPs electroactive polymers
  • adjustment of the ultrasound transducer 12 can be effected by adjusting the direction of an ultrasound beam from the ultrasound transducer 12 , for example by adjusting the phase of individual transducer array elements in the ultrasound transducer 12 .
  • FIG. 3 only shows the components required to illustrate this aspect of the invention, and in a practical implementation the apparatus 10 will comprise additional components to those shown.
  • the apparatus 10 may comprise a battery or other power supply for powering the apparatus 10 or means for connecting the apparatus 10 to a mains power supply, and/or a communication module for enabling the flow velocity measurements to be communicated to a remote computer.
  • FIG. 4 illustrates a scenario in which the invention can be used
  • FIG. 5 is a flow chart illustrating a method according to an embodiment.
  • FIG. 4( a ) shows an aorta 2 and an ultrasound transducer 12 that is positioned generally above the aorta 2 (for example in the suprasternal notch of the subject).
  • Vectors 28 indicate the direction of blood flow in the aorta 2
  • vectors 29 indicate the component of the velocity towards the ultrasound transducer 12 that can be measured by ultrasound transducer 12 .
  • the ultrasound transducer 12 may not have the ability to obtain spatial-time-velocity information at a plurality of depths, and thus the spatial-time-velocity profile is obtained by integrating the axial component of the velocity of the blood flow at all depths for each (2D) spatial position.
  • the ultrasound beam 30 from the ultrasound transducer 12 is shown in FIG.
  • the ultrasound transducer 12 is not correctly aligned with the aorta 2 (as shown by the center axis 32 of the ultrasound beam being offset by an angle ⁇ from the axis 34 of the aorta 2 ).
  • FIG. 4( b ) shows a spatial-time-velocity profile obtained using the ultrasound transducer 12 where the center axis 32 of the ultrasound beam is offset by an angle ⁇ from the axis 34 of the aorta 2 , and it can be seen that the center of the part of the spatial-time-velocity profile relating to the flow of blood in the aorta 2 is offset from the center axis 32 of the ultrasound beam 30 by a distance d.
  • This offset leads to off-center angles interacting with the high-velocity center of the aorta 2 , thereby ‘smearing’ the obtained spatial-time-velocity profile.
  • the method according to this embodiment of the invention therefore aims to detect this offset to enable an operator and/or the apparatus 10 to take corrective action.
  • an ultrasound transducer 12 that emits an ultrasound beam is used to measure flow velocity in a part of a body of a subject.
  • the ultrasound transducer 12 can be used on, in or near the suprasternal notch of a subject in order to measure the flow velocity of blood in the ascending aorta.
  • an ultrasound transducer can be used in a parasternal position to measure the flow velocity of the blood in the pulmonary artery.
  • flow velocity is measured in a plurality of spatial directions.
  • the ultrasound transducer 12 can measure the flow velocities at a particular depth, or the ultrasound transducer 12 can measure the flow velocities as the integral of the axial component of the velocity of the blood flow at all depths for each spatial position.
  • a spatial time velocity profile (for example as shown in FIGS. 2 and 4 ( b )) is formed for the part of the body from the measured flow velocities. That is, a spatial time velocity profile indicating the flow velocity for each of a number of different positions in the lateral and transverse directions is formed from the flow velocity measurements made in a plurality of directions from the ultrasound transducer 12 . Step 103 can be performed by control unit 14 .
  • the spatial time velocity profile is analyzed to determine a correction to the angle of the ultrasound beam emitted by ultrasound transducer 12 with respect to the subject.
  • the spatial-time-velocity profile is analyzed to identify a peak in the profile, and the difference between the position of the peak in the spatial time velocity profile and the center of the spatial time velocity profile (corresponding to the intersection of the axis 32 in FIG. 4( a ) with the profile) is determined.
  • the correction to the angle of the ultrasound beam is based on the determined difference.
  • Step 105 can be performed by control unit 14 .
  • the peak in the spatial time velocity profile identified in step 105 can be the maximum flow velocity in the spatial time velocity profile.
  • the method further comprises the step of providing feedback to an operator of the ultrasound transducer based on the determined correction to the angle.
  • the feedback can be provided in any suitable form using any suitable feedback means in user interface 18 .
  • the feedback can simply comprise an indication that the orientation of the ultrasound transducer 12 is not correct, in which case the operator can adjust the orientation of the ultrasound transducer 12 until the feedback indicates that the orientation is correct (i.e. when the peak in the spatial time velocity profile is within a threshold distance of or is at the center of the spatial time velocity profile).
  • the indication can be a visual, aural and/or tactile indication, or a combination thereof.
  • the feedback can be an indication of the rotation of the ultrasound transducer 12 with respect to the subject and/or an indication of the direction of the rotation of the ultrasound transducer 12 that is required in order to position the peak in the time velocity profile in, or nearer to, the center of the spatial time velocity profile.
  • the indication can be a visual, aural and/or tactile indication, or a combination thereof.
  • the method can further comprise determining a control signal (also referred to herein as a first control signal) for the ultrasound transducer 12 based on the determined correction to the angle that can be used to steer the ultrasound beam 30 produced by the ultrasound transducer 12 .
  • the ultrasound transducer 12 can then be controlled using the first control signal to adjust the direction in which the ultrasound beam 30 is emitted with respect to the subject such that the position of the peak in the time velocity profile is in, or nearer to, the center of the spatial time velocity profile.
  • the first control signal can be used to adjust the phase of individual transducer elements in the ultrasound transducer 12 in order to steer the ultrasound beam 30 in the required direction.
  • the method can further comprise determining a control signal (also referred to herein as a second control signal) for an orientation actuator 20 for the ultrasound transducer 12 based on the determined correction to the angle.
  • the orientation actuator 20 can then be controlled using the second control signal to adjust the orientation of the ultrasound transducer 12 with respect to the subject such that the position of the peak in the time velocity profile is in, or nearer to, the center of the spatial time velocity profile.
  • steps 101 - 105 can be repeated. If in step 105 it is determined that no correction is required to the angle of the ultrasound transducer 12 , then the spatial time velocity profile determined in step 103 can be analyzed to determine the flow velocity of the blood in the aorta 2 , and/or any other cardiac parameter that can be derived from the flow velocity of the blood, for example using an aortic flow model.
  • the method can further comprise the step of determining a correction to a value of a flow velocity in the spatial time velocity profile based on the determined correction to the angle.
  • This step can be performed in addition or alternatively to the providing feedback and determining control signal steps set out above.
  • information from the spatial time velocity profile can be used to estimate the residual angle error of the aorta 2 with respect to the ultrasound transducer 12 , and the blood flow velocity estimate from the profile can be corrected using this residual angle error (for example by multiplying by the secant of the angle).
  • this step is performed as an alternative to the providing feedback and determining control signal steps, this step has the advantage that an automatic ‘correction’ to the angle (or more directly the flow velocity itself) can be made without requiring the apparatus 10 to include an actuator 20 .
  • the ultrasound transducer 12 may be able to measure flow velocity at a plurality of depths so that spatial-time-velocity profiles can be determined for specific depths.
  • the step of analyzing can comprise analyzing one or more of the plurality of spatial time velocity profiles to determine a correction to the angle of the ultrasound beam. Where different corrections are determined from each spatial time velocity profile, the correction for one of them, or an average of the corrections, can be selected as the output of step 105 .
  • FIG. 6 illustrates another scenario in which the invention can be used.
  • FIG. 6 shows an aorta 2 and an ultrasound transducer 12 that is positioned in position A generally above the aorta 2 (for example in the suprasternal notch of the subject).
  • Vectors 28 indicate the direction of blood flow in the aorta 2
  • vectors 29 indicate the component of the velocity towards the ultrasound transducer 12 that can be measured by ultrasound transducer 12 .
  • the ultrasound transducer 12 can measure the spatial-time-velocity information individually at a plurality of depths (e.g. it can be a 3D ultrasound system), and thus separate spatial-time-velocity profiles are obtained. Three depths are shown in FIG.
  • the resulting spatial-time-velocity profiles are shown in the figure (labelled 36 - 1 , 36 - 2 and 36 - 3 ).
  • the ultrasound beam 30 from the ultrasound transducer 12 is shown, and in this illustration, the ultrasound transducer 12 is not correctly positioned with respect to the aorta 2 , as shown by the ultrasound transducer 12 being offset in the lateral and/or transverse directions by a distance P from an ideal measurement position B. It can be seen that in each of the profiles 36 - 1 , 36 - 2 and 36 - 3 , the peak in the flow velocity is offset from the center of the profile (with the peak being increasingly offset as the measurement depth decreases).
  • the method in FIG. 5 can be used to identify a correction to the angle of the ultrasound beam from the obtained spatial time velocity profiles.
  • the spatial time velocity profiles can be analyzed to determine a correction to the position of the ultrasound transducer 12 (in the lateral and/or transverse directions) with respect to the subject.
  • the correction can be based on the fitting of a geometric model of the arterial path (i.e. the line through the center of the artery) to the peaks (e.g. maximum) of the spatial time profiles, and determining the correction from the deviation of the model and the peaks in the spatial time profiles.
  • the complexity of the geometric model can depend on the number of profiles, with just two profiles, the model is a line, with three profiles, the model is a parabola, etc.
  • feedback can be provided to an operator of the ultrasound transducer based on the determined correction to the position.
  • the feedback can be provided in any suitable form using any suitable feedback means in user interface 18 .
  • the feedback can simply comprise an indication that the orientation and/or position of the ultrasound transducer 12 is not correct, in which case the operator can adjust the orientation and/or position of the ultrasound transducer 12 until the feedback indicates that the orientation and/or orientation is correct.
  • the indication can be a visual, aural and/or tactile indication, or a combination thereof.
  • the feedback can be an indication of the translation of the ultrasound transducer 12 with respect to the subject in the lateral and/or transverse directions and/or an indication of the direction of the translation of the ultrasound transducer 12 in the lateral and/or transverse direction that is required to improve the alignment of the probe with the artery (i.e. by reducing the distance between the axis 32 of the ultrasound transducer 12 and the center line of the artery).
  • the indication can be a visual, aural and/or tactile indication, or a combination thereof.
  • the method can further comprise determining a control signal (also referred to herein as a third control signal) for a position actuator 20 for the ultrasound transducer 12 based on the determined correction to the position.
  • the position actuator 20 can then be controlled using the control signal to adjust the position of the ultrasound transducer 12 with respect to the subject in the lateral and/or transverse directions and thereby improve the alignment of the ultrasound transducer 12 with the artery.
  • the second and third control signal can be the same signal.
  • steps 101 - 105 can be repeated. If in step 105 it is determined that no correction is required to the position of the ultrasound transducer 12 , then one or more of the spatial time velocity profiles determined in step 103 can be analyzed to determine the flow velocity of the blood in the aorta 2 , and/or any other cardiac parameter that can be derived from the flow velocity of the blood, for example using an aortic flow model.
  • a computer program may be stored/distributed on a suitable medium, such as an optical storage medium or a solid-state medium supplied together with or as part of other hardware, but may also be distributed in other forms, such as via the Internet or other wired or wireless telecommunication systems. Any reference signs in the claims should not be construed as limiting the scope.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Radiology & Medical Imaging (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Hematology (AREA)
  • Cardiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Epidemiology (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Primary Health Care (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
US16/082,305 2016-03-23 2017-03-23 Method and apparatus for improving the measurement of flow velocity of blood Active 2038-01-27 US11020093B2 (en)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
EP16162037.2 2016-03-23
EP16162037 2016-03-23
EP16162037 2016-03-23
PCT/EP2017/056951 WO2017162802A2 (fr) 2016-03-23 2017-03-23 Procédé et appareil pour améliorer la mesure de la vitesse d'écoulement du sang

Publications (2)

Publication Number Publication Date
US20190090852A1 US20190090852A1 (en) 2019-03-28
US11020093B2 true US11020093B2 (en) 2021-06-01

Family

ID=55589764

Family Applications (1)

Application Number Title Priority Date Filing Date
US16/082,305 Active 2038-01-27 US11020093B2 (en) 2016-03-23 2017-03-23 Method and apparatus for improving the measurement of flow velocity of blood

Country Status (5)

Country Link
US (1) US11020093B2 (fr)
EP (1) EP3432801B1 (fr)
JP (1) JP6721707B2 (fr)
CN (1) CN109640827B (fr)
WO (1) WO2017162802A2 (fr)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3669787A1 (fr) * 2018-12-19 2020-06-24 Koninklijke Philips N.V. Unité de transducteur à ultrasons dotée d'une fonction de guidage de friction
CN110988853B (zh) * 2019-12-23 2023-08-01 中煤科工集团重庆研究院有限公司 基于有效峰值的超声波到达时间计算方法
CN113558659B (zh) * 2021-07-30 2023-07-04 重庆安酷科技有限公司 一种高精度超声波肺功能检测仪及其检测方法

Citations (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4509526A (en) 1983-02-08 1985-04-09 Lawrence Medical Systems, Inc. Method and system for non-invasive ultrasound Doppler cardiac output measurement
US4796634A (en) 1985-08-09 1989-01-10 Lawrence Medical Systems, Inc. Methods and apparatus for monitoring cardiac output
US5052395A (en) 1987-11-16 1991-10-01 Waters Instruments, Inc. Non-invasive ultrasonic pulse doppler cardiac output monitor
GB2257254A (en) 1991-06-14 1993-01-06 British Tech Group Monitoring fluid-carrying conduits
US5409010A (en) 1992-05-19 1995-04-25 Board Of Regents Of The University Of Washington Vector doppler medical devices for blood velocity studies
US5575289A (en) 1992-07-17 1996-11-19 British Technology Group Limited Flowmeters
US5928153A (en) 1998-11-24 1999-07-27 Hui-Hua Chiang Ultrasound flow velocity and Doppler angle measurement method and system
WO1999066835A1 (fr) 1998-06-24 1999-12-29 Uscom Pty Limited Moniteur ultrasonore de sorties cardiaques
US6261233B1 (en) 1996-01-05 2001-07-17 Sunlight Medical Ltd. Method and device for a blood velocity determination
WO2004004557A1 (fr) 2002-07-03 2004-01-15 Uscom Pty Ltd Procedes et appareil de mesure pression volume integres
WO2004004559A1 (fr) 2002-07-03 2004-01-15 Uscom Pty Ltd Transducteur d'ondes entretenues a profil bas dispose sur une poitrine
WO2004012618A1 (fr) 2002-08-06 2004-02-12 Uscom Pty Ltd Systeme d'analyse de flux sanguin
WO2004064631A1 (fr) 2003-01-22 2004-08-05 Uscom Limited Procede et systeme pour determiner des caracteristiques sanguines
WO2004064642A1 (fr) 2003-01-22 2004-08-05 Uscom Limited Transducteur ultrasonique
WO2004084735A1 (fr) 2003-03-25 2004-10-07 Uscom Limited Determination d'information de parametres a partir d'un signal de flux sanguin
US20050159664A1 (en) 2004-01-20 2005-07-21 Uscom Limited Method and apparatus for determining blood flow characteristics in small infants
WO2006024088A1 (fr) 2004-08-30 2006-03-09 Uscom, Limited Procede et appareil pour definir des intervalles cardiaques de temps
WO2006096915A1 (fr) 2005-03-15 2006-09-21 Uscom Limited Systeme et procede automatiques de suivi de flux
WO2007134394A1 (fr) 2006-05-22 2007-11-29 Uscom Limited Système et méthode de mesure cardiaque
US20090292208A1 (en) 2008-03-03 2009-11-26 Jeffrey Jr R Brooke Automated detection of asymptomatic carotid stenosis
US8473239B2 (en) 2009-04-14 2013-06-25 Maui Imaging, Inc. Multiple aperture ultrasound array alignment fixture
US20130345566A1 (en) * 2012-05-09 2013-12-26 The Regents Of The University Of Michigan Linear magnetic drive transducer for ultrasound imaging
WO2014107769A1 (fr) 2013-01-14 2014-07-17 Uscom Limited Système et procédé de surveillance combinée de la circulation sanguine et de la pression artérielle
US20140276072A1 (en) 2013-03-15 2014-09-18 B-K Medical Aps Doppler ultrasound imaging

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1287740C (zh) * 1995-08-08 2006-12-06 通用电器横河医疗系统株式会社 测定组织或血液的运动速度的方法和超声波诊断装置
JP4311904B2 (ja) * 2000-05-08 2009-08-12 ブレインスゲート リミテッド 血液脳関門及び脳血流の性質を改善するための翼口蓋神経節の刺激方法及び装置
JP4217023B2 (ja) * 2002-02-25 2009-01-28 一郎 佐久間 血管内皮計測装置
US7591787B2 (en) * 2005-09-15 2009-09-22 Piero Tortoli Method for removing Doppler angle ambiguity
JP5231768B2 (ja) * 2006-08-28 2013-07-10 株式会社東芝 超音波診断装置および超音波診断装置のデータ処理プログラム
JP5276465B2 (ja) * 2009-02-16 2013-08-28 日立アロカメディカル株式会社 超音波診断装置及び医療システム
CN101884526B (zh) * 2009-05-13 2012-12-26 深圳先进技术研究院 基于超声血流信息的动脉血压测量装置
ITGE20090070A1 (it) * 2009-08-31 2011-03-01 Esaote Spa Metodo e dispositivo per il rilevamento e la visualizzazione di informazioni emodinamiche in particolare del flusso ematico nelle vene, mediante ultrasoni
KR20120090170A (ko) * 2011-02-07 2012-08-17 삼성전자주식회사 초음파 측정장치 및 그 제어방법
CN103415246B (zh) * 2011-03-01 2015-11-25 皇家飞利浦有限公司 用于针对血流测量定位多普勒超声换能器的方法和设备以及用于血流测量的系统
CN103635144B (zh) * 2011-06-30 2015-12-02 皇家飞利浦有限公司 用于自动超声多普勒角及流速估计的方法与装置
WO2014003071A1 (fr) * 2012-06-27 2014-01-03 株式会社東芝 Dispositif d'échographie et procédé de correction de données d'image
JP6366589B2 (ja) * 2012-09-27 2018-08-01 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 超音波による狭窄評価のための自動化されたバイプレーンpwワークフロー
KR20150031091A (ko) * 2013-09-13 2015-03-23 삼성메디슨 주식회사 가이드 라인을 이용한 초음파 정보 제공 방법 및 장치

Patent Citations (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4509526A (en) 1983-02-08 1985-04-09 Lawrence Medical Systems, Inc. Method and system for non-invasive ultrasound Doppler cardiac output measurement
US4796634A (en) 1985-08-09 1989-01-10 Lawrence Medical Systems, Inc. Methods and apparatus for monitoring cardiac output
US5052395A (en) 1987-11-16 1991-10-01 Waters Instruments, Inc. Non-invasive ultrasonic pulse doppler cardiac output monitor
GB2257254A (en) 1991-06-14 1993-01-06 British Tech Group Monitoring fluid-carrying conduits
US5406948A (en) 1991-06-14 1995-04-18 British Technology Group Ltd. Flowmeters
US5409010A (en) 1992-05-19 1995-04-25 Board Of Regents Of The University Of Washington Vector doppler medical devices for blood velocity studies
US5575289A (en) 1992-07-17 1996-11-19 British Technology Group Limited Flowmeters
US6261233B1 (en) 1996-01-05 2001-07-17 Sunlight Medical Ltd. Method and device for a blood velocity determination
WO1999066835A1 (fr) 1998-06-24 1999-12-29 Uscom Pty Limited Moniteur ultrasonore de sorties cardiaques
US6565513B1 (en) 1998-06-24 2003-05-20 Uscom Pty Ltd. Ultrasonic cardiac output monitor
US5928153A (en) 1998-11-24 1999-07-27 Hui-Hua Chiang Ultrasound flow velocity and Doppler angle measurement method and system
WO2004004559A1 (fr) 2002-07-03 2004-01-15 Uscom Pty Ltd Transducteur d'ondes entretenues a profil bas dispose sur une poitrine
WO2004004557A1 (fr) 2002-07-03 2004-01-15 Uscom Pty Ltd Procedes et appareil de mesure pression volume integres
WO2004012618A1 (fr) 2002-08-06 2004-02-12 Uscom Pty Ltd Systeme d'analyse de flux sanguin
WO2004064631A1 (fr) 2003-01-22 2004-08-05 Uscom Limited Procede et systeme pour determiner des caracteristiques sanguines
WO2004064642A1 (fr) 2003-01-22 2004-08-05 Uscom Limited Transducteur ultrasonique
WO2004084735A1 (fr) 2003-03-25 2004-10-07 Uscom Limited Determination d'information de parametres a partir d'un signal de flux sanguin
US20050159664A1 (en) 2004-01-20 2005-07-21 Uscom Limited Method and apparatus for determining blood flow characteristics in small infants
WO2006024088A1 (fr) 2004-08-30 2006-03-09 Uscom, Limited Procede et appareil pour definir des intervalles cardiaques de temps
WO2006096915A1 (fr) 2005-03-15 2006-09-21 Uscom Limited Systeme et procede automatiques de suivi de flux
WO2007134394A1 (fr) 2006-05-22 2007-11-29 Uscom Limited Système et méthode de mesure cardiaque
US20090292208A1 (en) 2008-03-03 2009-11-26 Jeffrey Jr R Brooke Automated detection of asymptomatic carotid stenosis
US8473239B2 (en) 2009-04-14 2013-06-25 Maui Imaging, Inc. Multiple aperture ultrasound array alignment fixture
US20130345566A1 (en) * 2012-05-09 2013-12-26 The Regents Of The University Of Michigan Linear magnetic drive transducer for ultrasound imaging
WO2014107769A1 (fr) 2013-01-14 2014-07-17 Uscom Limited Système et procédé de surveillance combinée de la circulation sanguine et de la pression artérielle
US20140276072A1 (en) 2013-03-15 2014-09-18 B-K Medical Aps Doppler ultrasound imaging

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Bercoff, et al., "Ultrafast Compound Doppler Imaging: Providing Full Blood Flow Characterization"; IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 58, No. 1, Jan. 2011.
Girish Shirali, "Early Experience with a Miniaturized Three-dimensional Matrix Transducer in Children"; US Cardiovascular Disease, 2006.
Haugen, et al., "Blood Flow Velocity Profiles in the Aortic Annulus: A 3-Dimensional Freehand Color Flow Doppler Imaging Study"; Journal of the American Society of Echocardiography, vol. 15 No. 4, Apr. 2002.
J∅rgen Arendt Jensen "Linear description of ultrasound imaging systems"; Technical University of Denmark, Jun. 10, 1999.
Snider, et al., "Suprasternal Notch Echocardiography: A Two-dimensional Technique for Evaluating Congenital Heart Disease"; vol. 63, No. 1, Jan. 1981.

Also Published As

Publication number Publication date
JP6721707B2 (ja) 2020-07-15
WO2017162802A2 (fr) 2017-09-28
WO2017162802A3 (fr) 2019-01-03
CN109640827A (zh) 2019-04-16
JP2019509827A (ja) 2019-04-11
US20190090852A1 (en) 2019-03-28
CN109640827B (zh) 2022-04-29
EP3432801A2 (fr) 2019-01-30
EP3432801B1 (fr) 2020-05-06

Similar Documents

Publication Publication Date Title
CN109788940B (zh) 跟踪介入设备的特征
US11020093B2 (en) Method and apparatus for improving the measurement of flow velocity of blood
US6544181B1 (en) Method and apparatus for measuring volume flow and area for a dynamic orifice
JP7252130B2 (ja) 血管内流量及び圧力の測定
US11185305B2 (en) Intertial device tracking system and method of operation thereof
JPH11319107A (ja) 動物体内のカテ―テルの位置を測定するための方法、動物体内のカテ―テル先端部の位置を測定するための方法及びカテ―テル装置
JP2021118891A (ja) 介入デバイス認識
EP3471624B1 (fr) Système pour déterminer les paramètres hémodynamiques d'un patient
Prinz et al. Can echocardiographic particle image velocimetry correctly detect motion patterns as they occur in blood inside heart chambers? A validation study using moving phantoms
WO2013136573A1 (fr) Dispositif de détection de débit d'écoulement de liquide
WO2020030746A1 (fr) Positionnement de dispositif d'intervention à l'aide de signaux ultrasonores
EP3307171B1 (fr) Appareil d'imagerie à ultrasons
CN109561880B (zh) 用于确定心输出量的系统
US11241172B2 (en) Ultrasonic diagnostic apparatus and method of controlling the same
CN112566557A (zh) 相对于超声图像平面的介入设备定位
JP2016055093A (ja) 血圧計測装置及び血圧計測方法
JP4350994B2 (ja) 超音波診断装置
EP3833266B1 (fr) Positionnement d'un dispositif d'intervention utilisant des signaux ultrasonores
US11986339B2 (en) Tracking an interventional device within an ultrasound image plane
KR20230109628A (ko) 경정맥압 파형 측정 장치 및 방법
US20140121523A1 (en) Measuring apparatus and method thereof

Legal Events

Date Code Title Description
AS Assignment

Owner name: KONINKLIJKE PHILIPS N.V., NETHERLANDS

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:KROON, BART;BEZEMER, RICK;MATHEW, DENNY;SIGNING DATES FROM 20171013 TO 20171031;REEL/FRAME:046786/0698

FEPP Fee payment procedure

Free format text: ENTITY STATUS SET TO UNDISCOUNTED (ORIGINAL EVENT CODE: BIG.); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY

STPP Information on status: patent application and granting procedure in general

Free format text: APPLICATION UNDERGOING PREEXAM PROCESSING

STPP Information on status: patent application and granting procedure in general

Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION

STPP Information on status: patent application and granting procedure in general

Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER

STPP Information on status: patent application and granting procedure in general

Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS

STPP Information on status: patent application and granting procedure in general

Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT VERIFIED

STCF Information on status: patent grant

Free format text: PATENTED CASE