Modeling of Pharmaceutical Propagation

• US 20080097197A1
• Filed: 11/16/2005
• Published: 04/24/2008
• Est. Priority Date: 11/16/2004
• Status: Active Grant

• Current Assignee: Bayer Healthcare LLC (Bayer AG)
• Original Assignee: Medrad Incorporated (Bayer AG)
• Inventors: Kalafut, John, Uber III, Arthur
• Granted Patent: US 8,295,914 B2
• Time in Patent Office: Days
• US Class Current: 600/431
• CPC Class Codes: A61B 5/0044 for the heart
• CPC Class Codes: A61B 5/02028 Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
• CPC Class Codes: A61B 5/0205 Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
• CPC Class Codes: A61B 5/026 Measuring blood flow A61B3/1233, A61B3/1241 take precedence
• CPC Class Codes: A61B 5/029 Measuring or recording blood output from the heart, e.g. minute volume A61B8/065 takes precedence
• CPC Class Codes: A61B 5/0295 using plethysmography, i.e. measuring the variations in the volume of a body part as modified by the circulation of blood therethrough, e.g. impedance plethysmography
• CPC Class Codes: A61B 5/055 involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
• CPC Class Codes: A61B 5/7257 using Fourier transforms
• CPC Class Codes: A61B 6/032 Transmission computed tomography [CT]
• CPC Class Codes: A61B 6/481 involving the use of contrast agents
• CPC Class Codes: A61B 6/504 for diagnosis of blood vessels, e.g. by angiography
• CPC Class Codes: A61B 6/507 for determination of haemodynamic parameters, e.g. perfusion CT
• CPC Class Codes: A61B 6/583 using calibration phantoms
• CPC Class Codes: A61B 8/13 Tomography A61B8/10, A61B8/12 take precedence; tomography for radiation diagnosis A61B6/02
• CPC Class Codes: A61B 8/481 involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream analysis of echo signal exploiting properties of a contrast enhancer in devices for short-range imaging with acoustic waves G01S7/52039
• CPC Class Codes: A61K 49/0002 General or multifunctional contrast agents, e.g. chelated agents
• CPC Class Codes: A61M 5/007 for contrast media
• CPC Class Codes: G01R 33/5601 involving use of a contrast agent for contrast manipulation, e.g. a paramagnetic, super-paramagnetic, ferromagnetic or hyperpolarised contrast agent
• CPC Class Codes: G16H 30/20 for handling medical images, e.g. DICOM, HL7 or PACS
• CPC Class Codes: G16H 30/40 for processing medical images, e.g. editing
• CPC Class Codes: G16H 50/50 : for simulation or modelling of medical disorders

Abstract

A method of modeling propagation of a pharmaceutical fluid in a patient, includes: collecting data corresponding to a time response curve resulting from injection of the fluid; and determining at least one mathematical model describing the data. The mathematical model can, for example, be a model which is not determined by a continuous or a discrete-time Fourier transform of the data. A method of controlling injection of a pharmaceutical fluid into a patient using an injector in a medical procedure, includes: collecting data corresponding to a patient response curve resulting from injection of the fluid; determining at least one mathematical model describing the data; and controlling the injector during the medical procedure to control injection of the fluid into the patient to create patient response at least in part on the basis of the mathematical model. A method of controlling injection of a contrast medium into a patient using an injector in a medical imaging procedure using an imaging scanner, includes: determining at least one mathematical model to predict a time enhancement response resulting from injection of the contrast medium; determining an injection protocol to approximate a predetermined time enhancement response in the patient by determining a constrained input solution to the mathematical model; and using the injection protocol to control the injector during the medical imaging procedure to control injection of the contrast medium into the patient to create an image of a region of interest.

Claims

1. A method of controlling an injector system for delivering a contrast enhancing fluid to a patient, comprising:

determining at least one patient transfer function for the patient based upon data specific to the patient, the at least one patient transfer function providing a time enhancement output for a given input;

determining a desired time enhancement output;

using the at least one patient transfer function to determine an injection procedure input; and

controlling the injector system at least in part on the basis of the determined injection procedure input.

2. The method of claim 1 wherein the injection procedure input is determined considering at least one operational limitation or constraint of the injector system.

3. The method of claim 2 wherein the at least one patient transfer function is determined using a system identification model comprising parameters related to physiological parameters of the patient.

4. The method of claim 3 wherein the system identification model is discretizable.

5. The method of claim 3 further comprising the steps of:

developing an initial patient transfer function using estimates of at least one physiological parameter of the patient;

performing an injection; and

revising the patient transfer function based upon at least one time enhancement output of the injection.

6. The method of claim 5 wherein at least one patient physiological parameter is measured from the at least one time enhancement output.

7. The method of claim 5 wherein the injection is a test injection performed prior to a diagnostic imaging procedure.

8. The method of claim 7 wherein time enhancement outputs resulting from the test injection are measured for at least two different regions of interest.

9. The method of claim 8 wherein differences between the time enhancement outputs provide a measure of at least one patient physiological parameter.

10. The method of claim 9 wherein the at least one patient physiological parameter is a parameter of the cardiopulmonary system

11. The method of claim 10 wherein the at least one patient physiological parameter is cardiac output, blood volume in a region, a rate transfer term or a transit delay.

12. The method of claim 8 wherein a first time enhancement output is measured in the ascending aorta or the descending aorta and a second time enhancement output is measured in the pulmonary artery trunk.

13. The method of claim 9 wherein a first time enhancement output is measured in the ascending aorta or the descending aorta and a second time enhancement output is measured in the pulmonary artery trunk.

14. The method of claim 13 wherein the at least one patient physiological parameter is cardiac output, blood volume in a region, a rate transfer term or a transit delay.

15. The method of claim 2 wherein the at least one patient transfer function is determined by the steps of:

collecting data corresponding to a time response curve resulting from injection of the fluid; and

determining at least one mathematical model describing the data.

16. The method of claim 15 wherein the mathematical model is not determined by a continuous or a discrete-time Fourier transform of the data.

17. The model of claim 15 wherein the model is a parametric model.

18. The method of claim 17 wherein the model is a moving average or an autoregressive moving average.

19. The method of claim 17 wherein the mathematical model assumes linearity and time invariance.

20. The method of claim 15 wherein the model is a non-parametric model determined by a spectral estimation technique.

21. The method of claim 20 wherein the spectral estimation technique is Welch'"'"'s method, Bartlett'"'"'s method, a multiple signal classification (MUSIC) method, or the Periodogram method.

22. The method of claim 15 wherein data is collected during at least one test injection prior to an imaging injection.

23. The method of claim 15 wherein the at least one patient transfer function is updated with data collected during the imaging injection.

24. The method of claim 3 wherein the at least one patient transfer function is updated with data collected during the imaging injection.

25. The method of claim 2 wherein the at least one patient transfer function is determined at least in part on the basis of at least one injection.

26. The method of claim 25 wherein the at least one injection is a test injection performed prior to a diagnostic imaging procedure.

27. The method of claim 26 wherein the test injection comprises injection of contrast medium followed by injection of a non-contrast fluid.

28. The method of claim 27 wherein the non-contrast fluid is injected at substantially the same volumetric flow rate as a flow rate of contrast medium preceding the injection of non-contrast fluid.

29. The method of claim 28 wherein the non-contrast fluid is saline.

30. The method of claim 26 wherein more than one test injection is performed.

31. The method of claim 30 where one test injection comprises injection of contrast medium only and another test injection comprises injection of contrast medium followed by injection of a non-contrast fluid.

32. The method of claim 2 wherein the injection procedure input is determined using an analytical solution or using a numerical, constrained optimization technique.

33. The method of claim 32 wherein the numerical, constrained optimization technique is a weighted least-squared numerical optimization.

34. The method of claim 32 wherein the injection procedure input is optimized to minimize the mass of a contrast enhancing agent in the contrast enhancing fluid delivered to the patient.

35. The method of claim 34 wherein the contrast enhancing agent is iodine, xenon or gadolinium.

36. The method of claim 1 wherein the contrast enhancing fluid is a CT contrast enhancing fluid, a MRI contrast enhancing fluid, an ultrasound enhancing imaging fluid or a radioactive contrast enhancing fluid.

37. The method of claim 1 wherein at least two patient transfer functions are determined and the injection procedure input is determined on the basis of one of the patient transfer functions.

38. The method of claim 37 wherein a first patient transfer function is determined using a system identification model comprising parameters related to physiological parameters of the patient and a second patient transfer function is determined using a mathematical model determined by collecting data corresponding to a time enhancement curve resulting from injection, the mathematical model describing the data.

39. The method of claim 37 wherein a determination is made as to which patient transfer function provides the best correlation between a given input and a resulting output.

40. A method of determining at least one patient physiological parameter from an imaging procedure, comprising the following steps:

measuring time enhancement outputs for at least two different regions of interest; and

determining at least one difference between the time enhancement outputs to provide a measure of the at least one patient physiological parameter.

41. The method of claim 40 wherein the at least one patient physiological parameter is a parameter of the cardiopulmonary system

42. The method of claim 41 wherein the at least one patient physiological parameter is cardiac output, blood volume in a region, a rate transfer term or a transit delay.

43. The method of claim 40 wherein a first time enhancement output is measured in the ascending aorta or the descending aorta and a second time enhancement output is measured in the pulmonary artery trunk.

44. The method of claim 43 wherein the at least one patient physiological parameter is cardiac output, blood volume in a region, a rate transfer term or a transit delay.

45. An injector system for the delivery of a fluid to a patients comprising:

an injector; and

a controller in communicative connection with the injector, the controller comprising;

at least one patient transfer function determined for the patient based upon data specific to the patient, the at least one patient transfer function providing a time enhancement output for a given input; and

a processor to determine an injection procedure input for a desired time enhancement output using the at least one patient transfer function.

46. The injector system of claim 45 wherein the injection procedure input is determined considering at least one operational limitation or constraint of the injector.

47. The injector system of claim 45 wherein the injection procedure input is determined using an analytical solution or a numerical, constrained optimization technique.

48. The injector system of claim 47 wherein the numerical, constrained optimization technique is a weighted least-squared numerical optimization.

49. The injector system of claim 47 wherein the injection procedure input is optimized to minimize the mass of a contrast enhancing agent in the fluid delivered to the patient.

50. The injector system of claim 49 wherein the contrast enhancing agent is iodine, xenon or gadolinium.

51. The injector system of claim 45 wherein the fluid is a CT contrast enhancing fluid, a MRI contrast enhancing fluid, an ultrasound enhancing imaging fluid or a radioactive contrast enhancing fluid.

52-86. -86. (canceled)

87. A method of controlling an injector in a medical imaging procedure using an imaging scanner, comprising:

determining at least one mathematical model to predict a time enhancement response resulting from injection of a contrast medium;

determining an injection protocol to approximate a predetermined time enhancement response in the patient by determining a constrained input solution to the mathematical model; and

using the injection protocol to control the injector during the medical imaging procedure to control injection of the contrast medium into the patient to create an image of a region of interest.

88. The method of claim 87 further comprising the step of changing the injection protocol as a result of feedback regarding the time enhancement response during the imaging procedure.

89. The method of claim 87 further comprising the step of changing the injection protocol as a result of data regarding at least one patient physiological parameter during the imaging procedure.

90. The method of claim 87 wherein the step of determining an injection protocol to approximate the predetermined time enhancement response is accomplished using a numerical solver or a numerical optimizer.

91. The method of claim 87 wherein the constrained input solution to the mathematical model is constrained by at least one operational limitation of the injector.

92. The method of claim 87 wherein the constrained input solution to the mathematical model is constrained by at least one operational limitation related to patient safety or comfort.

93. The method of claim 87 wherein the injection of the contrast medium is commenced at one time and an image scan of the region of interest is commenced at a second time determined at least in part on the basis of the mathematical model.

94. The method of claim 93 wherein the second time is determined on the basis of a prediction of a time of attainment of a predetermined enhancement level as determined by the mathematical model.

95. The method of claim 87 wherein the at least one mathematical model is a patient transfer function for the patient based upon data specific to the patient, the at least one patient transfer function providing a time enhancement output for a given input.

96. The method of claim 95 wherein the patient transfer function is determined using a system identification model comprising parameters related to physiological parameters of the patient or is determined using a mathematic identification model determined by collecting data corresponding to a time enhancement curve resulting from injection of the patient, wherein the mathematical identification model describes the data.

97-98. -98. (canceled)

99. A method of controlling injection of a pharmaceutical fluid into a patient using an injector having a controller in communicative connection with a computer memory in a medical procedure, comprising:

collecting data corresponding to a patient response curve resulting from injection of the fluid;

choosing at least one mathematical model from a plurality of mathematical models stored in the computer memory to describe the data;

adapting the model to the collected data; and

controlling the injector via the controller during the medical procedure to control injection of the fluid into the patient to create patient response at least in part on the basis of the mathematical model.

100. A system for creating an image of a region of interest of a patient comprising:

an imaging device for measuring a property of the patient over a region of interest;

an injector for injecting a pharmaceutical into the patient;

at least one standard region that is also measured by the imaging device; and

a computation algorithm that corrects measurements of the property over the patient'"'"'s region of interest based upon the measurements in the standard region.

101. The system of claim 100 wherein the at least one standard region is outside of the patient.

102. The system of claim 100 where wherein the at least one standard region is a region of the patient.

103. A method of creating an image of a region of interest in a patient comprising:

measuring a property of the patient over the region of interest using an imaging device;

measuring at least one standard region with the imaging device; and

correcting measurements of the property over the patient'"'"'s region of interest based upon the measurements in the standard region.

104. The method of claim 103 wherein the at least one standard region is outside of the patient.

105. The method of claim 103 wherein the at least one standard region is a region of the patient.

106. An imaging system, comprising:

an imager to create an image of a region of interest of a patient;

an injector adapted to inject a contrast medium; and

a controller in operative communication with the injector to control the injector, the controller comprising.

at least one patient transfer function determined for the patient based upon data specific to the patient, the at least one patient transfer function providing a time enhancement output for a given input; and

a processor to determine an injection procedure input for a desired time enhancement output using the at least one patient transfer function.