Abstract:
A method and apparatus for evaluating the condition of the heart muscle utilizing a single roughly collimated radiation detector which is non-critically positioned over the heart muscle. A small amount of radioactive material is injected into the bloodstream and detected during the first transit through the heart. The background radiation level is empirically derived from, and may be utilized to correct, the aforementioned first transit data. Numerous parameters affecting the heart muscle may be determined by analysis of the resultant data.

Description:
BACKGROUND OF THE INVENTION 
     This invention relates to a method and apparatus for evaluating the condition of the heart muscle in living beings, and more specifically to methods which utilize radioactive materials. 
     DESCRIPTION OF THE PRIOR ART 
     Attempts to evaluate the condition of the heart muscle, utilizing radioactive materials, are known in the prior art. One problem associated with such attempts involves the &#34;background&#34; radiation generated by radioactive material which is deposited into the body tissues and organs surrounding the heart during such testing. The inability to accurately categorize a portion of the radioactivity detected near the heart as attributable to this &#34;background&#34; level has long prevented accurate evaluation of the heart muscle utilizing these techniques. Several solutions to this problem have been attempted in the art. One example of an attempted solution may be seen in U.S. Pat. No. 3,221,731, issued to M. Annis et al. The Annis et al system utilizes two radiation detectors placed over the patient&#39;s body. One detector is utilized to detect radioactive activity in the heart, while the other detector is positioned away from the heart, to detect background radiation. The amount of background radiation is then subtracted from the radiation detected at the heart to achieve a more accurate reading. 
     A second example of a method utilized to distinguish background radiation may be seen in the &#34;nuclear stethoscope&#34; manufactured by BIOS, of Valhalla, New York. The &#34;nuclear stethoscope&#34; attempts to avoid background radiation problems as much as possible by utilizing a tightly focused or collimated radiation detector. While this tight focus does tend to decrease the background radiation, there are several problems inherent in this approach. The &#34;nuclear stethoscope&#34; requires a second reading in a location away from the heart to obtain a background level. Additionally, the sharply focused detector requires a critical positioning on the heart, a task made doubly difficult by individual variations in heart location and attitude. For a discussion of these problems see &#34;Potential Pitfalls of the Nuclear Stethoscope,&#34; Michael J. Zema M.D. et al, J. B. Lippincott Company, 1980. 
     SUMMARY OF THE INVENTION 
     It is, therefore, one object of the invention to provide an apparatus for more accurate evaluation of the condition of the heart muscle. 
     It is another object of the invention to provide an apparatus for evaluation of the heart muscle which does not require critical positioning. 
     It is yet another object of the invention to provide an apparatus for evaluation of the heart muscle which requires a single detector positioned on the patient&#39;s body for a short period of time. 
     It is another object of this invention to provide an apparatus for evaluation of the heart muscle which requires injection of a small amount of radioactive material. 
     The foregoing objects are achieved as is now described. A roughly collimated radiation detector is non-critically positioned over the heart muscle. A small amount of radioactive material is injected into the bloodstream of the patient and the first transit of the material through the patient&#39;s heart is detected and analyzed. The background radiation level is empirically derived from the first transit data, and may be utilized to correct the data. The resultant data may then be analyzed to determine many parameters utilized to analyze heart functions. 
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS 
     The novel features believed characteristic of the invention are set forth in the appended claims. The invention itself, however, as well as a preferred mode of use, further objects and advantages thereof, will best be understood by reference to the following detailed description of an illustrative embodiment when read in conjunction with the accompanying drawings, wherein: 
     FIG. 1 is a partially diagrammatic and partially cross sectional view of the apparatus of the invention; 
     FIG. 2 is a block diagram of the major components of the scaler of FIG. 1; 
     FIG. 3a is an average cardiac curve obtained by the apparatus of FIG. 1; 
     FIG. 3b is an average cardiac curve for a single heartbeat; 
     FIG. 4 is a block diagram of the major components of the processor of FIG. 1. 
     FIG. 5 is a flow chart diagram of the method and apparatus of the invention. 
    
    
     DETAILED DESCRIPTION 
     With reference now to FIG. 1, a patient&#39;s couch 10 is indicated, above which is located a single scanning detector. The scanning detector of a preferred embodiment utilizes a roughly focused lead collimator 12, capable of collimating incoming radiation from an area roughly the size of the human heart, at the average depth of the heart. By utilizing a roughly collimated scanning detector, the apparatus of this invention becomes relatively insensitive to variations in heart position and attitude. This non-critical positioning capability greatly enhances the utility of the apparatus of this invention by greatly diminishing the amount of training required to successfully operate the system. Above collimator 12 is a scintillation crystal 14, preferably of the sodium iodide, thallium activated type (NaI(Tl)), capable of converting the gamma radiation emitted by radioactive material into flashes of light. In the disclosed embodiment, an intravenous injection of approximately one mCi of Tc-99m is utilized to provide the radioactive agent. Associated with scintillation crystal 14, at the end of a light tight light pipe 16, is photomultiplier 18. Each light flash generated by scintillation crystal 14 is converted into an electron pulse and multiplied by photomultiplier 18. High gain rates possible in such photomultipliers yield a current pulse well above background. 
     These current pulses are applied to scaler 22 for scaling and amplification, and then coupled to processor 24. Processor 24 determines the background radiation level in a manner to be described herein in detail, and generates a desired parameter of heart muscle function. Desired parameters may be electronically visually displayed, recorded, or printed via output device 26. 
     Referring now to FIG. 2, there is depicted a block diagram of scaler 22. Scaler 22 includes two major sections; a high voltage supply 58 and an amplifier/discriminator 44. Amplifier/discriminator 44 includes a preamp 36 which integrates and differentiates the electron current generated by the scanning detector, generating a voltage pulse, the height of which is proportional to the energy of the gamma radiation interacting in scintillation crystal 14. This voltage pulse is further amplified by amplifier 38 and is applied to discriminator 40. Discriminator 40 acts as a pulse height analyzer and and accepts those pulses which vary a selected amount as valid inputs. In the embodiment disclosed, the median pulse height is representative of unscattered Tc-99m; however those skilled in the art will realize that this may vary for alternatively selected radioactive agents. The output of discriminator 40 is applied to an interface circuit 42 to create a TTL compatible pulse, in a manner well known in the art, for processor 24. In the embodiment disclosed, the minimum pulse width was constrained to 0.5 microseconds. The counting system utilized should be capable of attaining 200,000 counts per second without saturation. While various interface circuits may be utilized to compensate for alternative processor circuitry, the remaining elements of scaler 22 may be provided by commercially available systems. One example of such a system is the Model MS-2, Mini-Scaler, manufactured by the Eberline Instrument Corporation of Santa Fe, New Mexico. 
     With reference now to FIG. 3a, there is depicted an average cardiac curve obtained utilizing the apparatus of this invention. The cardiac curve of FIG. 3a plots the relative count rate versus time. The first peak of radioactivity, occurring at time t 1  represents the bolus of radioactive material arriving at the right ventricle of the subject&#39;s heart. The radioactive level diminishes as the blood containing the injected material leaves the heart for the lungs, and then increases again as the material reenters the heart, into the left ventricle at time t 2 . Those skilled in this art will appreciate that the difference between time t 2  and t 1  approximates the so-called &#34;pulmonary transit&#34; time. Individual variations present on the cardiac curve represent variations due to individual heartbeats. The method of determining the background radiation disclosed involves a calculation of the left ventricle ejection fraction (LVEF); however, in addition to the aforementioned pulmonary transit time, other cardiac functions may be determined utilizing the apparatus of this invention. 
     The left ventricle ejection fraction (LVEF) is a particularly useful parameter for determining the condition of the heart muscle. Referring now to FIG. 3b, we see a cardiac curve typical of a single heartbeat. The high points of the curve (when radiation is at the highest level) represent those periods when the heart is relaxed, and full of blood. This is referred to as the diastolic volume. The low point represents the contraction of the heart, and the radioactivity present in the heart at that point is proportional to the systolic volume, or the volume of blood present in the ventricle at the end of contraction. The LVEF is the ratio of the volume of blood pumped by the pumping chamber to the maximum volume of blood held by the pumping chamber. Thus, it should be clear that the LVEF will be a ratio, expressed as a decimal between zero and one. Referring again to FIG. 3b, we see that after correcting for background radiation due to the surrounding tissue and blood, (assumed to be a constant value B for this example) the LVEF may be expressed in terms of the labels of FIG. 3b as equation (1). ##EQU1## 
     It should now be apparent that background correction is necessary for an accurate evaluation of the heart muscle. 
     Referring again to FIG. 3a, the method whereby the first transit data may be corrected to account for background radiation will be explained in detail. As explained earlier, at time t 2 , the activity is seen to peak in the left ventricle. An assumption is made that at time t 2  the activity in the right ventricle is negligible. Therefore, at any point between time t a  and time t b , the detected radiation may be said to be attributable to the activity in the left ventricle of the heart at that point and the background radiation present in surrounding tissues and organs. Examining the cardiac curve of FIG. 3a, four representative relative peaks, P 1  through P 4 , which occur during the time period t a  -t b , are selected. These relative peaks represent those times when the heart is relaxed and full of blood. Similarly, four relative valleys, V 1  through V 4 , are selected. Also seen in FIG. 3a is a curve B(t) which is assumed to approximate the values of background radiation during the selected time period. B(t) may be assumed to be expressible as a general second order function, as seen in Equation (2). 
     
         B(t)=B.sub.2 t.sup.2 +B.sub.1 t+B.sub.0                    (2) 
    
     Next, by any convenient method, the peak values, and valley values are each fitted to a similar second order function to create a single curve. These upper and lower curves are also assumed to be expressible as general second order functions as seen in Equations (3) and (4). 
     
         U(t)=U.sub.2 t.sup.2 +U.sub.1 t+U.sub.0                    (3) 
    
     
         L(t)=L.sub.2 t.sup.2 +L.sub.1 t+L.sub.0                    (4) 
    
     In the embodiment disclosed, processor 24 implements a simple least squares approximation of the desired functions, in the manner well known in the art. Having approximated the second order curves which best fit the peak values and the valley values, and having assumed a second order curve which fits the background function, the left ventricle ejection fraction may be expressed, as was discussed with respect to FIG. 3b, as seen in Equation (5). ##EQU2## Further, the left ventricle ejection fraction may be assumed to be constant during the period in question, and therefore, the coefficients of the second order function B(t) may be expressed in terms of Equation (6), as simply derived from Equation (5). ##EQU3## It will be appreciated by anyone skilled in the mathematics art that for every assumed value of LVEF in Equation (6), a corresponding value may be found for B n . Therefore, a second relationship is necessary to correctly derive the background relationship. An advantageous relationship to consider at this point is a relationship of background activity to the activity taking place in the chamber of the heart. Recalling that background activity is generated as a result of radioactive material being deposited into the tissues and organs surrounding the heart, it should be apparent that the background activity at any given time period should be minimally correlated with, or independent of, the chamber activity during that time period. The instantaneous chamber activity may be easily expressed as the difference between the peak curve and the valley curve, of Equation (3) minus Equation (4). This relationship is expressed as the envelope function E(t), as seen in Equation (7). 
     
         E(t)=U(t)-L(t)                                             (7) 
    
     Assuming that no correlation exists between the envelope function E(t) and the background function B(t), then the relationship seen in Equation (8) is true. 
     
         E{B(t)·E(t)}=E{B(t)}·E{E(t)}             (8) 
    
     Where E {X} represents the expectation value of X over some unspecified range. 
     By forming a ratio R of the two sides of Equation (8) and integrating from time t a  to time t b , it can be seen that where B(t) and E(t) are minimally correlated, as assumed, the ratio R will approach unity. Utilizing these assumptions, it is a simple task for processor 24 to assume a value for LVEF and to calculate the corresponding values for B(t) and R, over the time period in question. A simple iterative approach will allow each possible value of LVEF and the corresponding value of B(t) to be derived. Having determined each possible corresponding set of LVEF values and B(t) values, ratio R is examined for each set to select the value of LVEF and B(t) for which ratio R most nearly approaches unity (indicating the minimal correlation relationship initially assumed). 
     Referring now to FIG. 4, a block diagram of processor 24 is depicted. The count data generated by scaler 22 is stored in a Random-Access-Memory 60. A Central Processing Unit 64 operates upon the data in Random-Access-Memory 60 in a manner specified by an appropriate algorithm stored in Read-Only-Memory 62. Central Processing Unit 64 is preferably utilized to determine the second order functions necessary to describe the peak and valley curves, and to calculate and evaluate the equations required to determine the LVEF and background functions. A desired parameter may be electronically displayed, recorded or printed via output device 26. The prototype device, constructed for evaluation, utilized a standard micro or home-type computer to implement processor 24; however it is contemplated that a special purpose microprocessor programmed with a similar algorithm will possess those characteristics of processor 24 which are necessary for commercial implementation. 
     Referring now to FIG. 5, there is depicted a flow chart of the operation of the disclosed system. The flow chart begins with the intravenous injection of the radioactive agent into the patient. The detector/scaler is then utilized to detect and quantize the levels of radiation present in the patient&#39;s heart during the first transit of the radioactive material. The next two functions, acquiring data and locating the appropriate peaks and valleys, are controlled, in the embodiment disclosed, with the general purpose micro-computer. One example of a program for such data acquisition and plotting may be seen in Table I. The fitting of the peak and valley values to appropriate second order functions and the determination of which value for LVEF yields a minimal correlation between background activity and chamber activity may be accomplished by a program similar to that contained in Table II. It should be noted that in the embodiment disclosed, LVEF has been constrained to those values which fall between 0.10 and 0.90. It has been determined that the majority of patients will fall within these constraints; however, it would be a simple matter for those skilled in the art to narrow those constraints, or remove them entirely, with small adjustments to the program disclosed in Table II. 
     ALTERNATIVE EMBODIMENTS 
     Although the invention has been described with reference to a specific embodiment, this description is not meant to be construed in a limiting sense. Numerous technical advantages taught in the disclosed system will, upon reference to this description by persons skilled in the art, lend themselves to alternative embodiments of the invention. As an example, while the disclosed embodiment is described with reference to the calculation of the left ventricle ejection fraction, it should be apparent that the right ventricle ejection fraction may be similarly derived. Additionally, the pulmonary transit time parameter described herein may be derived utilizing the system disclosed, as well as a cardiac output index or other similar heart function parameters. It is therefore contemplated that the appended claims will cover any such modifications or embodiments that fall within the true scope of the invention. ##SPC1## ##SPC2##