EP1505908A2 - Nonuniform rotational distortion (nurd) reduction - Google Patents

Nonuniform rotational distortion (nurd) reduction

Info

Publication number
EP1505908A2
EP1505908A2 EP03737552A EP03737552A EP1505908A2 EP 1505908 A2 EP1505908 A2 EP 1505908A2 EP 03737552 A EP03737552 A EP 03737552A EP 03737552 A EP03737552 A EP 03737552A EP 1505908 A2 EP1505908 A2 EP 1505908A2
Authority
EP
European Patent Office
Prior art keywords
image
instruction
image vector
computing
vector
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.)
Granted
Application number
EP03737552A
Other languages
German (de)
French (fr)
Other versions
EP1505908B1 (en
Inventor
Shashidhar Sathyanarayana
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.)
Boston Scientific Ltd Barbados
Boston Scientific Corp
Original Assignee
Boston Scientific Ltd Barbados
Boston Scientific Corp
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 Boston Scientific Ltd Barbados, Boston Scientific Corp filed Critical Boston Scientific Ltd Barbados
Publication of EP1505908A2 publication Critical patent/EP1505908A2/en
Application granted granted Critical
Publication of EP1505908B1 publication Critical patent/EP1505908B1/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/10—Image enhancement or restoration using non-spatial domain filtering
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/12—Diagnosis using ultrasonic, sonic or infrasonic waves in body cavities or body tracts, e.g. by using catheters
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4444—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
    • A61B8/4461—Features of the scanning mechanism, e.g. for moving the transducer within the housing of the probe
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8934—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a dynamic transducer configuration
    • G01S15/8938—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a dynamic transducer configuration using transducers mounted for mechanical movement in two dimensions
    • G01S15/894—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a dynamic transducer configuration using transducers mounted for mechanical movement in two dimensions by rotation about a single axis
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8977—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/80—Geometric correction
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/40—Analysis of texture
    • G06T7/41—Analysis of texture based on statistical description of texture
    • G06T7/42—Analysis of texture based on statistical description of texture using transform domain methods
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20048—Transform domain processing
    • G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30101—Blood vessel; Artery; Vein; Vascular

Definitions

  • the invention relates generally to medical imaging, and more particularly to reducing Nonuniform Rotational Distortion (NURD) in medical images.
  • NURD Nonuniform Rotational Distortion
  • imaging techniques such as ultrasound imaging are commonly used in medical procedures to obtain images of the inside of a patient's body.
  • IVUS intravascular ultrasound
  • images revealing the internal anatomy of blood vessels are obtained by inserting a catheter with an ultrasound transducer mounted on or near its tip into the blood vessel.
  • the ultrasound transducer is positioned in a region of the blood vessel to be imaged, where it emits pulses of ultrasound energy into the blood vessel and surrounding tissue.
  • a portion of the ultrasound energy is reflected off of the blood vessel wall and surrounding tissue back to the transducer.
  • the reflected ultrasound energy (echo) impinging on the transducer produces an electrical signal, which is used to form an image of the blood vessel.
  • FIG. 1 is a representation of an axial view of a rotating transducer 10 mounted on the tip of a prior art catheter 20.
  • the transducer 10 is coupled to a drive motor (not shown) via a drive cable 30 and rotates within a sheath 35 of the catheter 20.
  • the blood vessel 40 being imaged typically includes a blood region 45 and wall structures (blood-wall interface) 50 and the surrounding tissue.
  • a cross-sectional image of the blood vessel is obtained by having the transducer 10 emit a plurality of ultrasound pulses, e.g., 256, at different angles as it is rotated over one revolution.
  • FIG. 1 illustrates one exemplary ultrasound pulse 60 being emitted from the transducer 10.
  • the echo pulse 65 for each emitted pulse 60 received by the transducer is used to compose one radial line or "image vector" in the image of the blood vessel.
  • the transducer 10 is rotated at a uniform angular velocity so that the image vectors are taken at evenly spaced angles within the blood vessel 40.
  • An image processor (not shown) assembles the image vectors acquired during one revolution of the transducer 10 into a cross-sectional image of the blood vessel 40.
  • the image processor assembles the image vectors based on the assumption that the image vectors were taken at evenly spaced angles within the blood vessel 40, which occurs when the transducer 10 is rotated at a uniform angular velocity.
  • a uniform angular velocity for the transducer 10 This is because the transducer 10 is mechanically coupled to a drive motor (not shown), which may be located one to two meters from the transducer, via the drive cable 30.
  • the drive cable 30 must follow all the bends along the path of the blood vessel to reach the region of the blood vessel 40 being imaged. As a result, the drive cable 30 typically binds and/or whips around as it is rotated in the blood vessel 40.
  • FIG. 1 is a representation of a rotating transducer of a prior art catheter inside a blood vessel.
  • FIG. 2 is a flowchart illustration of an example embodiment of a new image processing method for reducing NURD in IVUS images acquired using a rotating transducer.
  • Described below is a new image processing method that reduces NURD in IVUS images acquired using a rotating transducer.
  • an IVUS image of a blood vessel the blood inside the blood vessel and the tissue surrounding the blood vessel have texture, which appear as speckles in the IVUS image.
  • the blood typically has a fine image texture and the surrounding tissue has a course image texture.
  • the image texture of the blood and the surrounding tissue should be fairly consistent throughout the image.
  • the transducer rotates at a nonuniform angular velocity, the image texture in the blood and the surrounding tissue becomes nonuniform.
  • the new imaging processing method corrects for NURD in an image, as explained further below.
  • an image processor receives an input image comprising a plurality of image vectors, e.g., 256 vectors.
  • the image vectors are mapped onto angles in the image based on the assumption that the image vectors were taken at uniformly spaced angles.
  • Each of the image vectors further comprises a plurality of pixels.
  • the value of each pixel corresponds to the amplitude of a received echo pulse that is reflected back to the transducer from a certain angle and radial distance with respect to the transducer.
  • the values of the pixels may be scaled according to a gray scale and/or a color scale.
  • a spectral measure of texture around each pixel is computed in the azimuthal direction. This may be accomplished by performing a one-dimensional Fourier transform on a set of pixels within a weighted window centered at the pixel.
  • the Fourier transform may be performed using standard signal processing techniques known to those of ordinary skill in the art.
  • the Fourier transform for each pixel produces a frequency spectrum that contains local textural information for the pixel.
  • the weight of the window used in the Fourier transform may be computed using the following equation:
  • Weight where w is the width of the window, ⁇ determines the drop off rate of the weight from the center of the window, and n is incremented from 1 to w.
  • the width w may be 16 pixels and ⁇ may be 4.
  • step 230 the mean frequency of the Fourier transform for each pixel is computed.
  • the mean frequency for each pixel provides a textural measure for the pixel with higher values indicating textural compression and lower values indicating textural blurring.
  • step 240 for each image vector, the average value of the mean frequency for the pixels in the image vector is computed.
  • the average frequency value for each image vector correlates with the relative angular velocity for the transducer at the image vector.
  • a high average frequency value indicates a relatively high angular velocity for the transducer at the image vector and a low average frequency value indicates a relatively low angular velocity for the transducer at the image vector.
  • the average frequency values for the image vectors is noted to be fairly constant.
  • step 250 the integral of the average frequency values for all the image vectors is computed with the integral normalized to a value of 2 ⁇ radians, which is the angle of one revolution of the transducer.
  • step 260 an estimate of the actual angle for each image vector is computed using the running value of the normalized integral at the image vector. This estimated angle for each image vector takes into account the fact that image vectors are not taken at uniformly spaced angles.
  • step 270 each image vector is remapped to its respective estimated angle to produce a NURD corrected image. In other words, NURD is reduced or eliminated by deriving an estimated angle for each image vector and using that estimated angle instead of the inaccurately assumed uniformly spaced angle.
  • the value of the width w and ⁇ used to compute weight of the window in step 220 may be optimized through normal experimentation.
  • a phantom e.g., made of rubber, having a known cross-sectional profile may be imaged using a rotating transducer.
  • the NURD algorithm may then be applied to the image of the phantom while adjusting the values of w and ⁇ until the NURD corrected image exhibits the least amount of NURD.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Theoretical Computer Science (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Medical Informatics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Probability & Statistics with Applications (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Image Processing (AREA)
  • Amplifiers (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Prostheses (AREA)
  • Oscillators With Electromechanical Resonators (AREA)

Abstract

A new image processing method reduces Nonuniform Rotational Distortion (NURD) in a medical image acquired using a rotating transducer. The image comprises a plurality of image vectors having texture. In a preferred embodiment, the image processing technique computes an average frequency of the texture for each image vector and estimates an angle for each image vector based on the average frequency for the respective image vector. The image processing technique then corrects for NURD by remapping each image vector to the estimated angle for the respective image vector.

Description

NONUNIFORM ROTATIONAL DISTORTION (NURD) REDUCTION
FIELD OF THE INVENTION
The invention relates generally to medical imaging, and more particularly to reducing Nonuniform Rotational Distortion (NURD) in medical images.
BACKGROUND
For purposes of diagnosis and treatment planning, imaging techniques such as ultrasound imaging are commonly used in medical procedures to obtain images of the inside of a patient's body. In intravascular ultrasound (IVUS) imaging, images revealing the internal anatomy of blood vessels are obtained by inserting a catheter with an ultrasound transducer mounted on or near its tip into the blood vessel. The ultrasound transducer is positioned in a region of the blood vessel to be imaged, where it emits pulses of ultrasound energy into the blood vessel and surrounding tissue. A portion of the ultrasound energy is reflected off of the blood vessel wall and surrounding tissue back to the transducer. The reflected ultrasound energy (echo) impinging on the transducer produces an electrical signal, which is used to form an image of the blood vessel. To obtain a cross-sectional image or "slice" of the blood vessel, the transducer must interrogate the vessel in all directions. This can be accomplished by mechanically rotating the transducer during imaging. FIG. 1 is a representation of an axial view of a rotating transducer 10 mounted on the tip of a prior art catheter 20. The transducer 10 is coupled to a drive motor (not shown) via a drive cable 30 and rotates within a sheath 35 of the catheter 20. The blood vessel 40 being imaged typically includes a blood region 45 and wall structures (blood-wall interface) 50 and the surrounding tissue.
A cross-sectional image of the blood vessel is obtained by having the transducer 10 emit a plurality of ultrasound pulses, e.g., 256, at different angles as it is rotated over one revolution. FIG. 1 illustrates one exemplary ultrasound pulse 60 being emitted from the transducer 10. The echo pulse 65 for each emitted pulse 60 received by the transducer is used to compose one radial line or "image vector" in the image of the blood vessel. Ideally, the transducer 10 is rotated at a uniform angular velocity so that the image vectors are taken at evenly spaced angles within the blood vessel 40. An image processor (not shown) assembles the image vectors acquired during one revolution of the transducer 10 into a cross-sectional image of the blood vessel 40. The image processor assembles the image vectors based on the assumption that the image vectors were taken at evenly spaced angles within the blood vessel 40, which occurs when the transducer 10 is rotated at a uniform angular velocity. Unfortunately, it is difficult to achieve and maintain a uniform angular velocity for the transducer 10. This is because the transducer 10 is mechanically coupled to a drive motor (not shown), which may be located one to two meters from the transducer, via the drive cable 30. The drive cable 30 must follow all the bends along the path of the blood vessel to reach the region of the blood vessel 40 being imaged. As a result, the drive cable 30 typically binds and/or whips around as it is rotated in the blood vessel 40. This causes the transducer 10 to rotate at a nonuniform angular velocity even though the motor rotates at a uniform angular velocity. This is a problem because the angles assumed by the image processor in assembling the image vectors into the cross-sectional image of the blood vessel 40 are different from the actual angles at which the image vectors were taken. This causes the cross-sectional image of the blood vessel to be distorted in the azimuthal direction. The resulting distortion is referred as Nonuniform Rotational Distortion (NURD).
Therefore, there is need for an image processing technique that reduces NURD in IVUS images acquired using a rotating transducer BRIEF DESCRIPTION OF THE DRAWINGS
The components in the figures are not necessarily to scale, emphasis instead being placed upon illustrating the concepts being discussed. All illustrations are intended to convey concepts, where relative sizes, shapes and other detailed attributes may be illustrated schematically rather than literally or precisely. FIG. 1 is a representation of a rotating transducer of a prior art catheter inside a blood vessel.
FIG. 2 is a flowchart illustration of an example embodiment of a new image processing method for reducing NURD in IVUS images acquired using a rotating transducer. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
Described below is a new image processing method that reduces NURD in IVUS images acquired using a rotating transducer. In an IVUS image of a blood vessel, the blood inside the blood vessel and the tissue surrounding the blood vessel have texture, which appear as speckles in the IVUS image. The blood typically has a fine image texture and the surrounding tissue has a course image texture. For an IVUS image taken with a transducer rotating at a uniform angular velocity, the image texture of the blood and the surrounding tissue should be fairly consistent throughout the image. However, when the transducer rotates at a nonuniform angular velocity, the image texture in the blood and the surrounding tissue becomes nonuniform. In regions of the image where the angular velocity of the transducer speeds up, the image texture becomes compressed in the azimuthal direction. In regions of the image where the angular velocity of the transducer slows down, the image texture becomes expanded, e.g., smeared out, in the azimuthal direction. Therefore, the degree of texture compression/expansion in the image yields information about the relative angular velocity of the transducer during imaging. Using this principle, the new imaging processing method corrects for NURD in an image, as explained further below.
Turning now to FIG. 2, an example embodiment of a new image processing method for reducing NURD will be described. In step 210, an image processor receives an input image comprising a plurality of image vectors, e.g., 256 vectors. The image vectors are mapped onto angles in the image based on the assumption that the image vectors were taken at uniformly spaced angles. Each of the image vectors further comprises a plurality of pixels. The value of each pixel corresponds to the amplitude of a received echo pulse that is reflected back to the transducer from a certain angle and radial distance with respect to the transducer. The values of the pixels may be scaled according to a gray scale and/or a color scale.
In step 220, a spectral measure of texture around each pixel is computed in the azimuthal direction. This may be accomplished by performing a one-dimensional Fourier transform on a set of pixels within a weighted window centered at the pixel. The Fourier transform may be performed using standard signal processing techniques known to those of ordinary skill in the art. The Fourier transform for each pixel produces a frequency spectrum that contains local textural information for the pixel. The weight of the window used in the Fourier transform may be computed using the following equation:
Weight = where w is the width of the window, χ determines the drop off rate of the weight from the center of the window, and n is incremented from 1 to w. As an example, the width w may be 16 pixels and χ may be 4.
In step 230, the mean frequency of the Fourier transform for each pixel is computed. The mean frequency for each pixel provides a textural measure for the pixel with higher values indicating textural compression and lower values indicating textural blurring.
In step 240, for each image vector, the average value of the mean frequency for the pixels in the image vector is computed. The average frequency value for each image vector correlates with the relative angular velocity for the transducer at the image vector. A high average frequency value indicates a relatively high angular velocity for the transducer at the image vector and a low average frequency value indicates a relatively low angular velocity for the transducer at the image vector. For a transducer rotating at a constant angular velocity, the average frequency values for the image vectors is noted to be fairly constant.
In step 250, the integral of the average frequency values for all the image vectors is computed with the integral normalized to a value of 2π radians, which is the angle of one revolution of the transducer. In step 260, an estimate of the actual angle for each image vector is computed using the running value of the normalized integral at the image vector. This estimated angle for each image vector takes into account the fact that image vectors are not taken at uniformly spaced angles. In step 270, each image vector is remapped to its respective estimated angle to produce a NURD corrected image. In other words, NURD is reduced or eliminated by deriving an estimated angle for each image vector and using that estimated angle instead of the inaccurately assumed uniformly spaced angle.
The value of the width w and χ used to compute weight of the window in step 220 may be optimized through normal experimentation. For example, a phantom, e.g., made of rubber, having a known cross-sectional profile may be imaged using a rotating transducer. The NURD algorithm may then be applied to the image of the phantom while adjusting the values of w and χ until the NURD corrected image exhibits the least amount of NURD.
In the foregoing specification, the invention has been described with reference to a specific embodiment thereof. It will, however, be evident that various modifications and changes may be made thereto without departing from the broader spirit and scope of the invention. For example, the reader is to understand that the specific ordering and combination of process actions shown in the process flow diagrams described herein is merely illustrative, and the invention can be performed using different or additional process actions, or a different combination or ordering of process actions. As another example, features known to those of skill in the art can be added to the embodiment. Other processing steps known to those of ordinary skill in the art may similarly be incorporated as desired. Additionally and obviously, features may be added or subtracted as desired. Accordingly, the invention is not to be restricted except in light of the attached claims and their equivalents.

Claims

CLAIMS WHAT IS CLAIMED IS:
1. A method for reducing Nonuniform Rotational Distortion (NURD) in an image, said image comprising a plurality of image vectors, each image vector having texture and each image vector being mapped to an angle in the image, the method comprising: computing an average frequency of the texture for each image vector; estimating an angle for each image vector based on the average frequency for the respective image vector; and remapping each image vector to the estimated angle for the respective image vector.
2. The method of claim 1 , wherein each image vector comprises a plurality of pixels, each pixel representing the amplitude of an echo pulse reflected from a certain image depth.
3. The method of claim 2, wherein the step of computing the average frequency for each image vector further comprises: computing a mean frequency of the texture for each pixel in each image vector; and computing an average of the mean frequency for the pixels in each image vector.
4. The method of claim 3, wherein the step of computing the mean frequency for each pixel further comprises: performing a Fourier transform around each pixel; and computing a mean of the Fourier transform for each pixel.
5. The method of claim 1 , wherein the step of estimating the angle for each image vector further comprises: computing an integral of the average frequency for all of the image vectors normalizing the integral to a predetermined value; and estimating the angle for the each image vector based on the value of the normalized integral at the respective image vector.
6. The method of claim 5, wherein the predetermined value is 2π radians.
7. The method of claim 3, wherein the step of estimating the angle for each image vector further comprises: computing an integral of the average frequency for all of the image vectors; normalizing the integral to a predetermined value; and estimating the angle for the each image vector based on the value of the normalized integral at the respective image vector.
8. The method of claim 7, wherein the predetermined value is 2π radians.
9. A computer program product that includes a medium useable by a processor, the medium comprising a sequence of instructions which, when executed by the processor, causes the processor to execute a method for reducing Nonuniform Rotational Distortion (NURD) in an image, the computer program product comprising: an instruction for receiving an input image, the input image comprising a plurality of image vectors, each image vector having texture and each image vector being mapped to an angle in the image; an instruction for computing an average frequency of the texture for each image vector in the input image; an instruction for estimating an angle for each image vector based on the average frequency for the respective image vector; and an instruction for producing an output image by remapping each image vector to the estimated angle for the respective image vector.
10. The computer program product of claim 9, wherein each image vector comprises a plurality of pixels, each pixel representing the amplitude of an echo pulse reflected from a certain image depth.
11. The computer program product of claim 10, wherein the instruction for computing the average frequency for each image vector further comprises: an instruction for computing a mean frequency of the texture for each pixel in each image vector; and an instruction for computing an average of the mean frequency for the pixels in each image vector.
12. The computer program product of claim 11 , wherein the instruction computing the mean frequency for each pixel further comprises: an instruction for performing a Fourier transform on the image around each pixel; and an instruction for computing a mean of the Fourier transform for each pixel.
13. The computer program product of claim 9, wherein the instruction for estimating the angle for each image vector further comprises: an instruction for computing an integral of the average frequency for all of the image vectors; an instruction for normalizing the integral to a predetermined value; and an instruction for estimating the angle for the each image vector based on the value of the normalized integral at the respective image vector.
14. The computer program product of claim 13, wherein the predetermined value is 2π radians.
15. A medical imaging system comprising:
(a) a processor;
(b) an interface to receive data for the processor to use to create a medical image; and
(c) a medium useable by the processor, the medium comprising a sequence of instructions which, when executed by the processor, causes the processor to create a medical image with reduced Nonuniform Rotational Distortion (NURD), the medium including
(i) an instruction for receiving an input image, the input image comprising a plurality of image vectors, each image vector having texture and each image vector being mapped to an angle in the image;
(ii) an instruction for computing an average frequency of the texture for each image vector in the input image;
(iii) an instruction for estimating an angle for each image vector based on the average frequency for the respective image vector; and
(iv) an instruction for producing an output image by remapping each image vector to the estimated angle for the respective image vector.
16. The medical imaging system of claim 15 further comprising a display to display the output image.
17. The medical imaging system of claim 15 further comprising a printer to print the output image.
18. The medical imaging system of claim 15 further comprising: a catheter; and an ultrasound transducer mounted on the catheter, the ultrasound transducer to be rotated by a motor, the ultrasound transducer to emit ultrasound waves and to receive reflected ultrasound waves, the ultrasound transducer to send reflected ultrasound waves to the interface.
19. The medical imaging system of claim 15, wherein each image vector comprises a plurality of pixels, each pixel representing the amplitude of an echo pulse reflected from a certain image depth.
20. The medical imaging system of claim 19, wherein the instruction for computing the average frequency for each image vector further comprises: an instruction for computing a mean frequency of the texture for each pixel in each image vector; and an instruction for computing an average of the mean frequency for the pixels in each image vector.
21. The medical imaging system of claim 20, wherein the instruction for computing the mean frequency for each pixel further comprises: an instruction for performing a Fourier transform on the image around each pixel; and an instruction for computing a mean of the Fourier transform for each pixel.
22. The medical imaging system of claim 15, wherein the instruction for estimating the angle for each image vector further comprises: an instruction for computing an integral of the average frequency for all of the image vectors; an instruction for normalizing the integral to a predetermined value; and an instruction for estimating the angle for the each image vector based on the value of the normalized integral at the respective image vector.
23. The medical imaging system of claim 22, wherein the predetermined value is 2π radians.
EP03737552A 2002-02-05 2003-01-23 Nonuniform rotational distortion (nurd) reduction Expired - Lifetime EP1505908B1 (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US72355 2002-02-05
US10/072,355 US7024025B2 (en) 2002-02-05 2002-02-05 Nonuniform Rotational Distortion (NURD) reduction
PCT/US2003/002182 WO2003067526A2 (en) 2002-02-05 2003-01-23 Nonuniform rotational distortion (nurd) reduction

Publications (2)

Publication Number Publication Date
EP1505908A2 true EP1505908A2 (en) 2005-02-16
EP1505908B1 EP1505908B1 (en) 2008-10-15

Family

ID=27659457

Family Applications (1)

Application Number Title Priority Date Filing Date
EP03737552A Expired - Lifetime EP1505908B1 (en) 2002-02-05 2003-01-23 Nonuniform rotational distortion (nurd) reduction

Country Status (9)

Country Link
US (1) US7024025B2 (en)
EP (1) EP1505908B1 (en)
JP (1) JP4242292B2 (en)
AT (1) ATE410956T1 (en)
AU (1) AU2003210648A1 (en)
CA (1) CA2471733C (en)
DE (1) DE60324137D1 (en)
ES (1) ES2315507T3 (en)
WO (1) WO2003067526A2 (en)

Families Citing this family (65)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4489770B2 (en) * 2003-07-21 2010-06-23 パイエオン インコーポレイテッド Method and system for identifying optimal images in a series of images depicting moving organs
US8047996B2 (en) 2005-10-31 2011-11-01 Volcano Corporation System and method for reducing angular geometric distortion in an imaging device
US9867530B2 (en) 2006-08-14 2018-01-16 Volcano Corporation Telescopic side port catheter device with imaging system and method for accessing side branch occlusions
US20080123911A1 (en) 2006-09-26 2008-05-29 Duc Lam Systems and Methods for Restoring a Medical Image Affected by Nonuniform Rotational Distortion
US8460195B2 (en) 2007-01-19 2013-06-11 Sunnybrook Health Sciences Centre Scanning mechanisms for imaging probe
CN101686827B (en) 2007-01-19 2014-08-13 桑尼布鲁克健康科学中心 Imaging probe with combined ultrasound and optical imaging device
WO2009009799A1 (en) 2007-07-12 2009-01-15 Volcano Corporation Catheter for in vivo imaging
US10219780B2 (en) 2007-07-12 2019-03-05 Volcano Corporation OCT-IVUS catheter for concurrent luminal imaging
US9596993B2 (en) 2007-07-12 2017-03-21 Volcano Corporation Automatic calibration systems and methods of use
US8197413B2 (en) * 2008-06-06 2012-06-12 Boston Scientific Scimed, Inc. Transducers, devices and systems containing the transducers, and methods of manufacture
US8351667B2 (en) * 2008-08-15 2013-01-08 Sti Medical Systems, Llc Methods of contrast enhancement for images having blood vessel structures
US9039626B2 (en) 2009-03-31 2015-05-26 Sunnybrook Health Sciences Centre Medical device with means to improve transmission of torque along a rotational drive shaft
GB2470942B (en) * 2009-06-11 2014-07-16 Snell Ltd Detection of non-uniform spatial scaling of an image
EP2632338B1 (en) * 2010-10-28 2015-07-01 Boston Scientific Scimed, Inc. Systems and methods for reducing non-uniform rotation distortion in ultrasound images
US11141063B2 (en) 2010-12-23 2021-10-12 Philips Image Guided Therapy Corporation Integrated system architectures and methods of use
US11040140B2 (en) 2010-12-31 2021-06-22 Philips Image Guided Therapy Corporation Deep vein thrombosis therapeutic methods
US9360630B2 (en) 2011-08-31 2016-06-07 Volcano Corporation Optical-electrical rotary joint and methods of use
US9292918B2 (en) 2012-10-05 2016-03-22 Volcano Corporation Methods and systems for transforming luminal images
US9324141B2 (en) 2012-10-05 2016-04-26 Volcano Corporation Removal of A-scan streaking artifact
US11272845B2 (en) 2012-10-05 2022-03-15 Philips Image Guided Therapy Corporation System and method for instant and automatic border detection
US9858668B2 (en) 2012-10-05 2018-01-02 Volcano Corporation Guidewire artifact removal in images
US9286673B2 (en) 2012-10-05 2016-03-15 Volcano Corporation Systems for correcting distortions in a medical image and methods of use thereof
US20140100454A1 (en) 2012-10-05 2014-04-10 Volcano Corporation Methods and systems for establishing parameters for three-dimensional imaging
US10070827B2 (en) 2012-10-05 2018-09-11 Volcano Corporation Automatic image playback
US9367965B2 (en) 2012-10-05 2016-06-14 Volcano Corporation Systems and methods for generating images of tissue
US10568586B2 (en) 2012-10-05 2020-02-25 Volcano Corporation Systems for indicating parameters in an imaging data set and methods of use
US9307926B2 (en) 2012-10-05 2016-04-12 Volcano Corporation Automatic stent detection
EP2904671B1 (en) 2012-10-05 2022-05-04 David Welford Systems and methods for amplifying light
US9840734B2 (en) 2012-10-22 2017-12-12 Raindance Technologies, Inc. Methods for analyzing DNA
EP2931132B1 (en) 2012-12-13 2023-07-05 Philips Image Guided Therapy Corporation System for targeted cannulation
US9730613B2 (en) 2012-12-20 2017-08-15 Volcano Corporation Locating intravascular images
JP6785554B2 (en) 2012-12-20 2020-11-18 ボルケーノ コーポレイション Smooth transition catheter
US10942022B2 (en) 2012-12-20 2021-03-09 Philips Image Guided Therapy Corporation Manual calibration of imaging system
US9709379B2 (en) 2012-12-20 2017-07-18 Volcano Corporation Optical coherence tomography system that is reconfigurable between different imaging modes
US10939826B2 (en) 2012-12-20 2021-03-09 Philips Image Guided Therapy Corporation Aspirating and removing biological material
US11406498B2 (en) 2012-12-20 2022-08-09 Philips Image Guided Therapy Corporation Implant delivery system and implants
EP2936626A4 (en) 2012-12-21 2016-08-17 David Welford SYSTEMS AND METHODS FOR REDUCING LIGHT WAVE LENGTH TRANSMISSION
EP2934280B1 (en) 2012-12-21 2022-10-19 Mai, Jerome Ultrasound imaging with variable line density
CA2896004A1 (en) 2012-12-21 2014-06-26 Nathaniel J. Kemp Power-efficient optical buffering using optical switch
US10058284B2 (en) 2012-12-21 2018-08-28 Volcano Corporation Simultaneous imaging, monitoring, and therapy
WO2014099672A1 (en) 2012-12-21 2014-06-26 Andrew Hancock System and method for multipath processing of image signals
US9612105B2 (en) 2012-12-21 2017-04-04 Volcano Corporation Polarization sensitive optical coherence tomography system
EP2934653B1 (en) 2012-12-21 2018-09-19 Douglas Meyer Rotational ultrasound imaging catheter with extended catheter body telescope
US10413317B2 (en) 2012-12-21 2019-09-17 Volcano Corporation System and method for catheter steering and operation
US10332228B2 (en) 2012-12-21 2019-06-25 Volcano Corporation System and method for graphical processing of medical data
US9486143B2 (en) 2012-12-21 2016-11-08 Volcano Corporation Intravascular forward imaging device
US10226597B2 (en) 2013-03-07 2019-03-12 Volcano Corporation Guidewire with centering mechanism
EP2965263B1 (en) 2013-03-07 2022-07-20 Bernhard Sturm Multimodal segmentation in intravascular images
US20140276923A1 (en) 2013-03-12 2014-09-18 Volcano Corporation Vibrating catheter and methods of use
JP2016521138A (en) 2013-03-12 2016-07-21 コリンズ，ドナ System and method for diagnosing coronary microvascular disease
WO2014159819A1 (en) 2013-03-13 2014-10-02 Jinhyoung Park System and methods for producing an image from a rotational intravascular ultrasound device
US9301687B2 (en) 2013-03-13 2016-04-05 Volcano Corporation System and method for OCT depth calibration
US11026591B2 (en) 2013-03-13 2021-06-08 Philips Image Guided Therapy Corporation Intravascular pressure sensor calibration
US10292677B2 (en) 2013-03-14 2019-05-21 Volcano Corporation Endoluminal filter having enhanced echogenic properties
US12343198B2 (en) 2013-03-14 2025-07-01 Philips Image Guided Therapy Corporation Delivery catheter having imaging capabilities
US10219887B2 (en) 2013-03-14 2019-03-05 Volcano Corporation Filters with echogenic characteristics
CN105208947B (en) 2013-03-14 2018-10-12 火山公司 Filter with echoing characteristic
JP6981967B2 (en) 2015-08-31 2021-12-17 ジェンテュイティ・リミテッド・ライアビリティ・カンパニーGentuity, LLC Imaging system including imaging probe and delivery device
US10952702B2 (en) 2016-06-21 2021-03-23 Canon U.S.A., Inc. Non-uniform rotational distortion detection catheter system
US12262872B2 (en) 2018-09-17 2025-04-01 Gentuity, Llc Imaging system with optical pathway
EP3962346A4 (en) 2019-04-30 2023-04-19 Gentuity LLC IMAGING PROBE WITH FLUID PRESSURE ELEMENT
WO2020237024A1 (en) 2019-05-21 2020-11-26 Gentuity, Llc Systems and methods for oct-guided treatment of a patient
CN111671468B (en) * 2020-05-28 2023-04-28 苏州博动戎影医疗科技有限公司 Asynchronous focusing dynamic compensation method for intravascular ultrasound single-array element transducer
US20230000321A1 (en) * 2021-03-01 2023-01-05 Gentuity, Llc Optical imaging system
CN117152026B (en) * 2023-10-30 2024-02-02 天津恒宇医疗科技有限公司 Intravascular ultrasound image processing method, device and equipment

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5361768A (en) * 1992-06-30 1994-11-08 Cardiovascular Imaging Systems, Inc. Automated longitudinal position translator for ultrasonic imaging probes, and methods of using same
US5383460A (en) * 1992-10-05 1995-01-24 Cardiovascular Imaging Systems, Inc. Method and apparatus for ultrasound imaging and atherectomy
US6424725B1 (en) * 1996-05-16 2002-07-23 Digimarc Corporation Determining transformations of media signals with embedded code signals
US5485845A (en) * 1995-05-04 1996-01-23 Hewlett Packard Company Rotary encoder for intravascular ultrasound catheter
US6095976A (en) * 1997-06-19 2000-08-01 Medinol Ltd. Method for enhancing an image derived from reflected ultrasound signals produced by an ultrasound transmitter and detector inserted in a bodily lumen
US5921934A (en) * 1997-11-25 1999-07-13 Scimed Life Systems, Inc. Methods and apparatus for non-uniform rotation distortion detection in an intravascular ultrasound imaging system
US5989191A (en) * 1998-06-19 1999-11-23 Hewlettt-Packard Company Using doppler techniques to measure non-uniform rotation of an ultrasound transducer
US6450964B1 (en) * 2000-09-05 2002-09-17 Advanced Cardiovascular Systems, Inc. Imaging apparatus and method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See references of WO03067526A3 *

Also Published As

Publication number Publication date
US7024025B2 (en) 2006-04-04
ATE410956T1 (en) 2008-10-15
WO2003067526A3 (en) 2004-11-25
US20030147551A1 (en) 2003-08-07
AU2003210648A8 (en) 2003-09-02
EP1505908B1 (en) 2008-10-15
CA2471733A1 (en) 2003-08-14
JP4242292B2 (en) 2009-03-25
DE60324137D1 (en) 2008-11-27
JP2006505293A (en) 2006-02-16
CA2471733C (en) 2011-06-07
ES2315507T3 (en) 2009-04-01
WO2003067526A2 (en) 2003-08-14
AU2003210648A1 (en) 2003-09-02

Similar Documents

Publication Publication Date Title
CA2471733C (en) Nonuniform rotational distortion (nurd) reduction
US9585636B2 (en) Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing method
EP1690230B1 (en) Automatic multi-dimensional intravascular ultrasound image segmentation method
US7856130B2 (en) Object recognition system for medical imaging
US5654509A (en) Control system that distinguishes between imaging and nonimaging environments in an ultrasound system
US20210022714A1 (en) Non-uniform rotation distortion (nurd) reduction in ultrasound imaging devices, systems, and methods
CN101467892B (en) Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image processing method
US20020167533A1 (en) Compound image display system and method
EP2613170A1 (en) Clutter suppression in ultrasonic imaging systems
US7852334B2 (en) Ultrasonic imaging apparatus, an image-processing apparatus, and an ultrasonic image-processing method
CN103391749B (en) Analysis of Mitral Regurgitation by Ultrasound Imaging
EP2082369B1 (en) Systems and methods for restoring a medical image affected by nonuniform rotational distortion
CN117094917B (en) Cardiovascular 3D printing data processing method
JP5069022B2 (en) Method and system for accurate time delay estimation for use in ultrasound imaging
CN103379864B (en) For the wall filter that ultrasonic mitral reflux is analyzed
CN112826535B (en) Method, device and equipment for automatically positioning blood vessel in ultrasonic imaging
EP4666103A1 (en) Adaptively weighted spatial compounding for ultrasound image contrast enhancement
US9918699B2 (en) Seam elimination and motion compensation in imaging data
CN114066767A (en) Processing method for enhancing contrast of blood vessel ultrasonic image and image processing system
CN113491535A (en) Ultrasonic diagnostic apparatus, signal processing apparatus, and signal processing program
CN114010228A (en) Self-adaptive image acquisition method
KR20110039506A (en) Ultrasound system and method for correcting volume data
CN121616489A (en) Transcranial ultrasonic image self-adaptive optimization method and system based on noise morphology
CN115844457A (en) Intracardiac ultrasonic imaging system and method based on phase calibration
CN118986413A (en) Method for estimating included angle between ultrasonic beam and blood flow movement direction, calibration system and method

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20040707

AK Designated contracting states

Kind code of ref document: A2

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LI LU MC NL PT SE SI SK TR

AX Request for extension of the european patent

Extension state: AL LT LV MK RO

17Q First examination report despatched

Effective date: 20070628

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LI LU MC NL PT SE SI SK TR

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

Ref country code: GB

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

REF Corresponds to:

Ref document number: 60324137

Country of ref document: DE

Date of ref document: 20081127

Kind code of ref document: P

REG Reference to a national code

Ref country code: ES

Ref legal event code: FG2A

Ref document number: 2315507

Country of ref document: ES

Kind code of ref document: T3

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20090115

Ref country code: AT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20081015

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20081015

Ref country code: FI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20081015

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20090316

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: DK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20081015

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20081015

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20090115

Ref country code: MC

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20090131

Ref country code: CZ

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20081015

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

26N No opposition filed

Effective date: 20090716

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20081015

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LI

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20090131

Ref country code: CH

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20090131

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20090123

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20090116

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 20101215

Year of fee payment: 9

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LU

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20090123

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: HU

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20090416

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: TR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20081015

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CY

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20081015

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: FR

Payment date: 20120111

Year of fee payment: 10

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: BE

Payment date: 20120222

Year of fee payment: 10

Ref country code: IT

Payment date: 20120124

Year of fee payment: 10

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: NL

Payment date: 20120118

Year of fee payment: 10

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: ES

Payment date: 20120116

Year of fee payment: 10

BERE Be: lapsed

Owner name: BOSTON SCIENTIFIC LIMITED

Effective date: 20130131

REG Reference to a national code

Ref country code: NL

Ref legal event code: V1

Effective date: 20130801

GBPC Gb: european patent ceased through non-payment of renewal fee

Effective date: 20130123

REG Reference to a national code

Ref country code: FR

Ref legal event code: ST

Effective date: 20130930

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: NL

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20130801

Ref country code: BE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20130131

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GB

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20130123

Ref country code: FR

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20130131

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IT

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20130123

REG Reference to a national code

Ref country code: ES

Ref legal event code: FD2A

Effective date: 20140321

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: ES

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20130124

REG Reference to a national code

Ref country code: DE

Ref legal event code: R082

Ref document number: 60324137

Country of ref document: DE

Representative=s name: PFENNING MEINIG & PARTNER GBR, DE

REG Reference to a national code

Ref country code: DE

Ref legal event code: R082

Ref document number: 60324137

Country of ref document: DE

Representative=s name: PFENNING MEINIG & PARTNER GBR, DE

Effective date: 20150202

Ref country code: DE

Ref legal event code: R081

Ref document number: 60324137

Country of ref document: DE

Owner name: BOSTON SCIENTIFIC LIMITED, BM

Free format text: FORMER OWNER: BOSTON SCIENTIFIC LIMITED, ST. MICHAEL, BB

Effective date: 20150202

Ref country code: DE

Ref legal event code: R082

Ref document number: 60324137

Country of ref document: DE

Representative=s name: PFENNING, MEINIG & PARTNER MBB PATENTANWAELTE, DE

Effective date: 20150202

REG Reference to a national code

Ref country code: DE

Ref legal event code: R081

Ref document number: 60324137

Country of ref document: DE

Owner name: BOSTON SCIENTIFIC MEDICAL DEVICE LIMITED, IE

Free format text: FORMER OWNER: BOSTON SCIENTIFIC LIMITED, HAMILTON, BM

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: DE

Payment date: 20211130

Year of fee payment: 20

REG Reference to a national code

Ref country code: DE

Ref legal event code: R071

Ref document number: 60324137

Country of ref document: DE