METHOD AND APPARATUS FOR EXAMINING FLUID FLOW

Abstract: A system (B) introduces an enhancement agent into a fluid. The agent concentration in the fluid is monitored and the characteristics of the absorption curve representing the projected increase of concentration of the agent in the fluid are determined (34). A look-up table array (60) includes a plurality of look-up tables, each look-up table corresponding to one of a plurality of preselected absorption curves. A look-up table selection circuit (36) selects the look-up table which most closely corresponds to the projected absorption curve. An imaging means (A) generates a plurality of image representations at preselected intervals after commencement of enhancement agent introduction, which images are stored in an image memory array (144). Enhancement agent concentration values from corresponding pixels of each image representation are utilized to address the look-up table corresponding most closely to the fluid's absorption curve in order to retrieve precalculated partition coefficient, fluid flow rate, and confidence values. The retrieved values may be interpolated (80, 82) to compensate for the fluid's absorption curve falling among a plurality of curves or the measured enhancement agent concentration values falling between the preselected addressable values of the look-up tables. A display (96) displays images indicative of the partition coefficient, the fluid flow rate, and the confidence value.

Description: This invention relates to methods and apparatus for examining fluid flow. More particularly the invention relates to methods and apparatus for producing images indicative of fluid flow in a subject. The invention finds particular application in conjunction with computed tomography (CT) blood flow mapping of the brain and will be described with particular reference thereto. However, it should be understood that the invention also finds application in conjunction with other imaging modalities, such as digital x-ray, magnetic resonance, radiation and positron emission, ultrasound, and the like. Furthermore, the invention is applicable to examining blood flow in regions of a human or animal body other than the brain, and to examining fluid flow in subjects other than the human or animal body, e.g. inanimate objects. In the human brain, blood reaches the tissue in two modes, directly through the arteries and indirectly through other tissue. In normal, healthy brain tissue, blood reaches gray matter at 40 to 140 ml per 100 ml per minute. Gray matter tissue which receives less than 30 ml per 100 ml per minute is not adequately fed for proper functioning and may suffer irreperable damage. In white matter, cerebral blood flows are typically about one third of those for gray matter, with flows under about 10 ml per 100 ml per minute being considered inadequate. The early detection of brain regions with subnormal blood flows enables corrective action to be taken before the blood tissue is irreversibly damaged. One of the most common causes of insufficient feeding of the tissue is a blockage in the arterial blood flow. In the past, iodine was utilized as an enhancement agent injected into the blood to facilitate the location of arterial blockages, see for example EP-A-0088935. However, brain tissue membrane blocked the iodine enhancement agent from permeating the tissue. Because the iodine was unable to pass from the blood flow into the tissue, iodine was only able to enhance images of blood in arteries, capillaries, and veins. Iodine was unable to enhance representations of the actual profusion of blood into the tissues. Unlike iodine, xenon passes from the blood into the brain tissue. Thus, utilizing the xenon as an enhancement agent facilitates the imaging and measurement of blood profusion into the tissue. As the concentration of xenon gas in the patient's blood rises, the concentration of xenon gas in the brain tissue increases, asymptopically approaching an equilibrium concentration. The rate of the exponential xenon gas concentration increase in the tissue is indicative of the blood flow rate. The equilibrium concentration which is asymtopically approached is indicative of a partition coefficient, λ. The partition coefficient, which is different for different kinds of tissue, is defined as the ratio of the quantity of xenon in each unit volume or voxel of tissue to the quantity of xenon per like volume in blood. For gray matter, the partition coefficient is typically about 0.95 and for white matter is typically about 1.35. Partition coefficients which differ significantly from these values are indicative of sick or dying tissue. The xenon concentration in the tissue of the ith unit volume or voxel at a time t is described by a formula known as the Kety equation: where C is the tissue xenon concentration, C a is the blood xenon concentration, K is the tissue clearance or build-up rate, and f is the flow rate. The partition coefficient, λ, is related to the flow rate and the clearance or build-up rate by the equation: f=λK (2), where λ is the tissue-blood partition coefficient. The blood xenon concentration is readily monitorable, see for example EP-A-0180876. The tissue xenon concentration for a tissue in a given voxel can be calculated from the CT number or other data value of a pixel of a CT image corresponding to the given voxel. By taking several CT images at different times, with the blood xenon concentration known for times preceeding each image, one can theoretically solve the Kety equation to determine the partition coefficient and blood flow for the tissue compartment corresponding to each pixel. Typically, three to six images were taken. More particularly, the values or CT numbers from the corresponding pixels of each of the three to six images were iteratively fit to the "best" flow f and partition coefficient λ which, with the known C a (w), allowed comparative C(t) to be calculated using any of various conventional curve fitting techniques. Perhaps the most common approximation implemented was the "minimum chi-square" curve fitting criterion which required extended and time consuming calculations. The chi-square curve fitting technique determined a best fit flow, a best fit partition coefficient, and a fit or confidence value indicative of the closeness of the best fit. The curve fitting technique was repeated for each pixel of the images. It is to be appreciated that chi-square and other curve fitting techniques for fitting three to six data points with a curve, then deriving the slope, the end point which the curve is asymtopically approaching, and the degree of conformity to the curve or best fit was a time consuming operation. When this operation was repeated over 65,000 times to fit the CT numbers of corresponding pixels of a conventional 256 x 256 image to corresponding curves, the computational time became excessive, even on a high speed computer. To reduce the computation to an acceptable time, the image resolution was commonly reduced from 256 x 256 pixels to as little as 32 x 32. However, calculating the flow, partition coefficient, and fit for each of the over 1000 pixels of a 32 x 32 image still required up to ten minutes. US-A-4610258 discloses a method of determining tissue blood flow using xenon concentration measurements and the Kety equation wherein calculations are carried out in a simplified manner which can speed up production of an image considerably e.g. by a factor of three compared with prior art methods. The present invention contemplates a new and improved technique for more rapidly and more accurately determining the flow, the partition coefficient, and the fit or confidence value from CT or other image data. In accordance with one aspect of the present invention there is provided an apparatus for producing images representing partition coefficients and fluid flow rates in a region of interest of a subject, the apparatus comprising: an enhancement agent means for introducing an enhancement agent into the fluid; a means for providing an indication of a concentration of the enhancement agent in the fluid; an imaging means for generating electronic image representations of at least the region of interest, each image representation including a plurality of pixel values, each pixel value being indicative of enhancement agent concentration in a corresponding subregion of the region of interest; an image representation storing means for storing a plurality of image representations of the selected region of interest, the stored image representations being generated at preselected sampling intervals; and means for using said pixel values and said enhancement agent concentrations in the fluid to obtain corresponding fluid flow and partition coefficient values, characterised in that said apparatus further includes: an absorption curve projecting means for computing parameters representative of an absorption curve indicative of enhancement agent absorption by the fluid projected from the provided indications of the enhancement agent concentration in the fluid over time; a look-up table means including a plurality of look-up tables, each preprogrammed with fluid flow and partition coefficient values which have a previously calculated relationship to preselected enhancement agent concentrations at the preselected sampling intervals for respective preselected fluid enhancement agent absorption curves; and a look-up table access means for selectively accessing the look-up table means with at least the computed parameters to select at least the look-up table which corresponds most closely to the projected absorption curve and the image pixel value for selected subregions to retrieve the corresponding preprogrammed partition coefficient and fluid flow values relating to the selected subregions from the selected look-up table. The apparatus normally further includes means for displaying images representing the retrieved partition coefficient and/or fluid flow values. In an apparatus according to the invention the imaging means is suitably an x-ray examination apparatus, e.g. a CT scanner. In one particular apparatus according to the invention the the look-up table selecting means includes a multidimensional look-up table which is addressed by the computed parameters to retrieve a designation of at least the most nearly corresponding preselected absorption curve. Preferably an apparatus according to the invention further includes: a reference image storage means for storing a reference image representation of the region of interest generated without the enhancement agent in the fluid; and a subtraction means for subtracting the reference image representation from the image representation generated at each of the selected sampling intervals to create difference image representations each having difference image pixel values that are indicative of enhancement agent concentration in the corresponding subregion of the region of interest at the corresponding selected sampling interval, the look-up tables being addressed by the corresponding difference pixel values. An apparatus according to the invention may also further include an interpolating means for interpolating retrieved partition coefficient and fluid flow values to compensate for at least one of the differences in the preselected enhancement agent concentrations, sampling intervals, absorption curves, and the image pixel values, the projected absorption curve, and the sampling intervals at which the image representations were actually generated. According to a second aspect of the invention there is provided a method of determining at least a selected one of partition coefficient and flow value of each of a plurality of voxels of a region of interest of a subject, the method comprising: at a commencement time, starting introducing an enhancement agent into the subject; measuring concentration of the enhancement agent in a fluid of the subject; generating a plurality of image representations each at a preselected sampling interval subsequent to the commencement time, each image representation being defined by a plurality of pixels corresponding to voxels of the region of interest and including a pixel value for each pixel, which pixel value is indicative of an enhancement agent concentration in the corresponding voxel; and using said pixel values and enhancement agent concentrations to obtain corresponding values of the selected one of partition coefficient and flow value; characterised by computing parameters representative of an absorption curve indicative of enhancement agent absorption by the fluid over time projected from the measured enhancement agent concentrations; and accessing a look-up table array including a plurality of look-up tables, each look-up table corresponding to a precalculated absorption curve and being addressable by pixel values generated at preselected sampling intervals, with at least the computed parameters to select the look-up table which most closely corresponds to the projected absorption curve and pixel values of image representations generated at the preselected sampling intervals to retrieve the selected one of precalculated partition coefficient and flow values relating to the voxels corresponding to said pixel values from the selected look-up table. The method suitably further includes prior to the commencement time: calculating corresponding partition coefficient flow rate values, absorption curves, sampling intervals, and pixel values in accordance with a Kety equation; storing the precalculated partition coefficient and flow values in the look-up table array; and defining at least the precalculated absorption curves and pixel values as addresses to index the look-up table array. In one particular method according to the invention the method further includes conducting the image representation generating step at the preselected sampling intervals. The method suitably further includes interpolating the retrieved partition coefficient and flow values from the plurality of look-up tables which most closely correspond to the projected absorption curve. In a further particular method according to the invention the look-up table array is addressable by preselected pixel values; and the method further includes interpolating partition coefficient and fluid flow values which are retrieved by addressing the look-up table array with preselected pixel values that are close to the pixel values from the generated image representations. A method according to the invention will normally further include the step of displaying man-readable images defined by pixels that correspond to the voxels of the region, a first image representing the partition coefficients and a second image representing flow values. A method according to the invention may also suitably include imaging the subject prior to the commencement time to create a reference image and subtracting the reference image from each subsequent image such that the stored pixel values are indicative of a difference image attributable substantially solely to the presence of the enhancement agent in each voxel. An apparatus and method according to the invention finds particular application wherein the fluid is blood and/or the enhancing agent is xenon. A primary advantage of the present invention is that it generates flow, partition coefficient, and confidence images quickly. Another advantage of the present invention is that it generates flow, partition coefficient, and confidence images with full, detailed resolution. In particular, flow, partition coefficient and confidence images for 256 × 256 pixel images can be done in as little as seven seconds as compared to ten minutes for 32 × 32 pixel images of the prior art curve fitting techniques. Another advantage of the present invention is that it enables interpolated and data improved images of 256 × 256 pixel images to be generated in less than four minutes as compared to 646 minutes to calculate 256 × 256 pixel images using a chi-square curve fitting technique. One method and apparatus in accordance with the invention will now be described, by way of example, with reference to the accompanying drawings in which:- Figure 1 is a diagrammatic illustration of the apparatus; Figure 2 illustrates a technique for determining a patient's blood absorption curve; Figure 3 illustrates a three dimensional data organisation for look-up tables of the appatatus of Figure 1; and, Figure 4 illustrates a data interpolation technique for interpolating data from the look-up tables of the apparatus of Figure 1. Referring to Figure 1, the apparatus includes an imaging means A, such as an axial tomographic scanner, arranged to selectively examine a subject and reconstruct image representations depicting properties of each voxel or subregion within an examined region of interest. An enhancement agent system B introduces selected amounts of an enhancement agent into the subject and derives an indication of the concentration of the enhancement agent in at least selected portions of the subject. The enhancement agent is selected such that its concentration or presence in the subject alters the image representations. In the preferred embodiment, the enhancement agent is xenon gas whose concentration is reflected in the reconstructed image representations. The resultant image representations are essentially the sum of two images -- an image depicting the xenon gas in the voxels of the regions of interest and an image of the tissue in the region of interest. A processing means C processes the image data from the scanner A collected with varying concentrations of the enhancement agent and processes the enhancement agent concentrations from the enhancement agent system B to derive images representing the partition coefficient or permeability, blood or other fluid flow rates, and a confidence value or degree of fit. The imaging means A may be a CT scanner, a digital x-ray scanner, a magnetic resonance imager, or other imaging device which generates data which can be reconstructed into a representation of an image of a region of interest. In the CT scanner embodiment, the scanner includes a gantry 10 which houses a rotating x-ray source and radiation detectors. Images from a CT scanner are indicative of x-ray absorption/transmission properties of tissue in each voxel of a planar region of interest or slice. In a magnetic resonance scanner, the gantry 10 houses appropriate electromagnetic coils and antennae. A patient supporting couch 12 selectively indexes a patient through the gantry such that one or more planar slices are selectively imageable. The enhancement agent supply system B includes a source 20 of the enhancement agent, such as a tank of xenon gas. A flow or pressure regulator 22 controls the supply rate of the xenon gas to be breathing mask 24 or other means for introducing the enhancement agent into the patient. An enhancement agent concentration means 26 determines concentrations of the enhancement agent within the patient. It is well known that there is a linear relationship between the xenon concentration in exhaled breath at the end of an exhalation cycle and the concentration of xenon gas in a patient's blood. In the preferred embodiment, the enhancement agent concentration means 26 measures the concentration of xenon gas at the end of a respiratory cycle and provides an output signal indicative thereof. Before the xenon gas or other enhancement agent is introduced to the patient, a reference scan is conducted to generate an image representation depicting the imaged region without the enhancement agent. The reference image is stored in a reference image memory 30. Upon completion of the reference image at a time t₀, the enhancement agent system B starts supplying xenon gas to the patient. For example, xenon may be substituted for 30% of the gases breathed by the patient. At the end of each respiratory cycle, generally every few seconds, the breath analyzer means 26 determines the concentration of xenon gas in the patient's blood. The increase in xenon concentration in the blood with time normally increases exponentially as illustrated in curve 32 of FIGURE 2. Generally, the data acquisition period lasts for several minutes, e.g. five to seven minutes. During this time, some 50 or more xenon concentration readings are commonly made by the breath analyzer means 26. A curve fitting means 34 utilizes known curve fitting techniques to project an absorption curve which is indicative of the absorption of xenon in the patient's blood. That is, the curve fitting means generates equations which describe the parameters of a curve 32 that most closely fits the discrete data points from the breath analyzer 26. Various conventional curve fitting techniques may be implemented. In a preferred embodiment, the curve fitting technique implements a double exponential fit. That is, from the breath analyzer data points, two exponential curves are determined, which two curves sum to produce a curve which follows the actual blood flow concentration. The double exponential fit technique requires the determination of three variables: a scaler multiplier f and two exponential values, a, b. That is: f(1-e -at ) + (1-f)(1-e -bt ) (3). Optionally, other curve fitting techniques may be implemented. However, more precise curve fitting techniques, such as those that describe the curve in terms of three or more parameters are preferred. The values of the three or other number of parameters which describe the absorption curve are used to address a blood curve selection means 36 such as a three dimensional look-up table. The three dimensional blood curve selection look-up table 36 is addressable only be preselected values of the three parameters. Normally, each of the three calculated parameter values falls between two preselected address values. In one embodiment, the preselected address value that is closest to the calculated value is used as the address. In another embodiment, each of the eight combinations of adjoining preselected address values are used to retrieve eight blood curves, which may all be the same. As the xenon concentration in the patient's blood is building, a plurality of images of the planar region are generated at preselected sampling intervals or times. As each image is reconstructed, it is stored in a temporary memory 40. A subtraction means 42 subtracts the xenon image from the reference image stored in the reference image memory means 30 to produce a difference or xenon concentration image representative only of xenon concentration. A difference image memory 44 stores a plurality of difference images, each corresponding to a different sampling interval. A first xenon enhanced scan is conducted a first preselected sampling interval t1 after the patient commences breathing xenon gas. This first produces a first image from which the reference image is subtracted to produce a first difference image that is stored in a first difference image memory 46. At a second preselected sampling interval or time after the patient starts breathing the xenon gas, t2, a second scan generates a second image from which the reference image is subtracted to produce a second difference image that is stored in a second image memory 48. In the preferred embodiment, a third difference image memory 50 stores a third difference image which is generated from a scan at a third sampling interval t3. Fourth and additional image memories are contemplated as may be appropriate to the degree of precision required. In the preferred embodiment, the first sampling interval is 1 1/2 minutes after t o , the second sampling interval t2 is 3 minutes after t o , and the third sampling interval t3 is 5 minutes after t o . Optionally, a fourth image may be taken at a fourth sampling interval of 7 minutes. Each difference image is defined by an array of pixels which each correspond to a voxel of the region of interest. Corresponding pixels of each difference image stored are those that depict the same voxel but at a different sampling interval or time. The concentration of xenon in the corresponding pixels can be expected to increase logarithmically with time, analogous to curve 32. The CT number or other value stored for the corresponding pixel of each difference memories 46, 48, 50 is proportional to the concentration of xenon gas in the corresponding voxel at the sampling interval t1, t2, t3, etc. If a voxel is totally within an unblocked artery, the concentration would be expected to follow the blood-xenon absorption curve 32 exactly. If the pixel is in tissue, the increase in xenon concentration might be expected to somewhat less. The exact concentration is determined by the Kety equation (1) discussed above. A look-up table means 60 is preprogrammed with a blood flow value and partition coefficient stored in each memory element. Each memory element is addressed by a corresponding combination of preselected addressing adsorption curves, pixel values, and sampling times. By addressing the look-up table means with the preselected addressing absorption curve, pixel values, and sampling times that are closest to the absorption curve projected by the curve projecting means 34, the corresponding pixel values of difference memory means 44, and the actual sampling interval, corresponding blood flow, and partition coefficient values are retrieved. Various techniques may be utilized to calculate the correspondence of the flow values and partition coefficients stored in the look-up table with the addressing values. In the preferred method for loading the look-up table, partition coefficients, blood flow rates, and absorption curves are selected which span the range normally encountered in human patients. The number of partition coefficients, blood flow rates, and absorption curves selected determines the precision with which the final answer is reached but increases the size of the look up table and the complexity of the calculations to fill it. For a given blood flow rate, partition coefficient and absorption curve, xenon concentration curves are calculated from the Kety equation. The calculated xenon concentration curves give the exact concentration at each of the preselected sampling intervals or times t1, t2, and t3, i.e. the concentration or pixel values which address each stored blood flow value and partition coefficient. Because all possible concentrations which might be measured at t1, t2, and t3, will not fall exactly on a theoretically calculated xenon concentration curve, many of the look-up table entries will remain unfilled. To fill these remaining look-up table entries, the CT number values theoretically predicted at each sampling time are varied randomly a small amount. These stochastic variations generate additional look-up table entries with varying likelihoods. The most likely flow and lambda values are then used to fill this matrix value and the degree of likelihood or "confidence" value is also stored. It is to be appreciated that the reliability of flow and partition coefficient values generated by large deviations in observed CT absorption are less reliable than flow and lambda data which corresponds precisely to one of the theoretical curves. To advise the viewer of the degree of fit to the theoretical curves, the confidence value is stored in the corresponding memory element along with the stored partition coefficient and flow rate. Looking now to a preferred organization of the look-up table means 60, a plurality of look-up tables are provided, each look-up table corresponding to one of a plurality of preselected blood curves, i.e. C a (w) in Equation (1), and xenon enhanced measurement times t1, t2, t3. The look-up table selection means 36 selects at least the look-up table which was precalculated in accordance with the absorption curve which is most like the patient's blood curve C a (w) and the selected sampling times t1, t2, t3. Each of the look-up tables is an n dimensional look-up table, where n is the number of difference images generated and stored in the image memory means 44 during a single examination. For simplicity of illustration, the organization of each individual look-up table is described in conjunction with three images per study, i.e. a three dimensional look-up table. As illustrated in Figure 3, C(t1), C(t2), and C(t3) are the coordinates of a three dimensional look-up table, where C(tn) is the xenon concentration at sampling time tn. The magnitude of the preselected addressing xenon concentrations, CT numbers, or pixel values along each axis function as the look-up table. It is to be appreciated that additional sampling intervals, e.g. t4, may also be accommodated. Further, additional look-up tables may be provided for different sampling intervals or numbers of sampling intervals. During a patient scan, a read control means 62 causes the pixel value or CT number corresponding to the same pixel in each image to be read from the difference images taken at sampling intervals t1, t2, and t3. As illustrated in FIGURE 3, the pixel value from the pixel of the t1 image which normally falls between two precalculated xenon concentration values defines a planar region of possible solutions 64. The thickness of the planar region is determined by the distance between adjacent precalculated addressing xenon concentrations or pixel values. The data value from the t2 image similarly defines a planar region 66 and the data value taken from the t3 image determines a planar region of possible solutions 68. With reference to FIGURES 3 and 4, the intersection of these three planar regions defines a volume or cube 70 within which point 72 defined by the three actual pixel values falls. The precalculated blood flow, partition coefficient and confidence values are stored in memory elements corresponding to each corner of the volume. The partition coefficient, flow value, and confidence values at the closest corner are retrieved and read out of the look-up tables 60. For greater accuracy, the partition coefficient, flow rate, and fit may be interpolated from values stored in the eight memory elements at the corners of the volume 70 by interpolating means 80. In a preferred embodiment, the distance between the point 72 and each corner is calculated. The partition coefficient, flow rate, confidence value from the memory element at each corner is weighted in inverse proportion to the calculated distance. The weighted retrieved values are summed. As indicated above, the theoretical blood-xenon absorption curve projected by absorption curve calculating means 34 in many instances will not describe one of the preselected addressing absorption curves precisely. Rather, the value may tend to fall between two or more of the precalculated absorption curves. Because the absorption curve is described by three variables in the preferred embodiment, variables may address a point within a volume defined by eight corners as discussed above in conjunction with FIGURE 4. In the preferred embodiment, the look-up tables corresponding to each of the eight corners absorption curves are also addressed by each pixel value to generate eight additional sets of partition coefficients, flow rates, and fit values. Each of the eight sets is interpolated by the first interpolating means 80. A second interpolating means 82 performs a weighted averaging of eight partition coefficients, flow rates, and confidence values from the first interpolating means. The weighting is again in inverse proportion to the distance between 72 and the corners. It is to be appreciated that the processing may be simplified when several cornes of absorption curve look-up table 36 define the same absorption curve. The sampling times t1, t2, and t3 may also be varied. To vary the sampling times, additional sets of look-up tables are provided for each additional sampling time combination. A third interpolating means may be provided to adjust for differences between the actual sampling time and the look-up tables actually provided. The interpolated partition coefficient values corresponding each pixel are stored in corresponding pixels of a partition coefficient image memory means 90. Analogously, the blood flow rates retrieved for each pixel are stored in a flow image memory means 92 and the retrieved fit or confidence values are stored in a fit or confidence image memory means 94. A display means 96 displays partition coefficient, flow rate, and fit or confidence images. The look-up tables 60, as well as the absorption curve select look-up table 36, are all described as three dimensional look-up tables. Optionally, other dimensions may be provided. Single dimensional look-up tables could be utilized but may produce suspect partition coefficients, flow values, and confidence value images. Dimensions higher than those can produce greater accuracy and higher confidence values. The invention has been described with reference to the preferred embodiment. Obviously, modifications and alterations will occur to others upon reading and understanding the preceding specification. It is intended that the invention include all such alterations and modifications insofar as they come within the scope of the appended claims or the equivalents thereof.

Claims: An apparatus for producing images representing partition coefficients and fluid flow rates in a region of interest of a subject, the apparatus comprising: an enhancement agent means (B, 20, 22, 24) for introducing an enhancement agent into the fluid; a means (26) for providing an indication of a concentration of the enhancement agent in the fluid; an imaging means (A, 30, 40, 42) for generating electronic image representations of at least the region of interest, each image representation including a plurality of pixel values, each pixel value being indicative of enhancement agent concentration in a corresponding subregion of the region of interest; an image representation storing means (44) for storing a plurality of image representations of the selected region of interest, the stored image representations being generated at preselected sampling intervals (t₁, t₂ etc); and means for using said pixel values and said enhancement agent concentrations in the fluid to obtain corresponding fluid flow and partition coefficient values characterised in that said apparatus further includes: an absorption curve projecting means (34) for computing parameters representative of an absorption curve indicative of enhancement agent absorption by the fluid projected from the provided indications of the enhancement agent concentration in the fluid over time; a look-up table means (60) including a plurality of look-up tables, each preprogrammed with fluid flow and partition coefficient values which have a previously calculated relationship to preselected enhancement agent concentrations at the preselected sampling intervals for respective preselected fluid enhancement agent absorption curves; and a look-up table access means (36, 62) for selectively accessing the look-up table means (60) with at least the computed parameters to select at least the look-up table which corresponds most closely to the projected absorption curve and the image pixel value for selected subregions to retrieve the corresponding preprogrammed partition coefficient and fluid flow values relating to the selected subregions from the selected look-up table. An apparatus according to Claim 1 further including a partition coefficient image memory means (90) for storing a partition coefficient image representation in which the retrieved partition coefficient values from the look-up table means (60) are stored as pixel values for pixels corresponding to the same subregion as the stored image representation pixels. An apparatus according to Claim 1 or Claim 2 further including a fluid flow image memory means (92) for storing a fluid flow image representation in which the retrieved fluid flow values from the look-up table means (60) are stored as pixel values for fluid flow image pixels corresponding to the same subregions as the stored image representation pixels. An apparatus according to Claim 2 or Claim 3 further including a display means (96) for displaying man-readable images, the display means (96) being operatively connected with the partition coefficient memory means (90) and/or the flow memory means (92) to display the partition coefficient and/or fluid flow image representations. An apparatus according to any one of the preceding claims wherein the imaging means (A, 30, 40, 42) includes an x-ray examination apparatus (A). An apparatus according to Claim 5 wherein the x-ray examination apparatus comprises a CT scanner. An apparatus according to any one of the preceding Claims wherein the look-up table selecting means (36) includes a multidimensional look-up table which is addressed by the computed parameters to retrieve a designation or at least the most nearly corresponding preselected absorption curve. An apparatus according to any one of the preceding Claims further including: a reference image storage means (30) for storing a reference image representation of the region of interest generated without the enhancement agent in the fluid; and a subtraction means (42) for subtracting the reference image representation from the image representation generated at each of the selected sampling intervals to create difference image representations each having difference image pixel values that are indicative of enhancement agent concentration in the corresponding subregion of the region of interest at the corresponding selected sampling interval, the look-up tables being addressed by the corresponding difference pixels values. An apparatus according to any one of the preceding claims further including an interpolating means (80, 82) for interpolating retrieved partition coefficient and fluid flow values to compensate for at least one of the differences in the preselected enhancement agent concentrations, sampling intervals, absorption curves, and the image pixel values, the projected absorption curve, and the sampling intervals at which the image representations were actually generated. An apparatus according to any one of the preceding claims wherein said fluid is blood. An apparatus according to any one of the preceding claims wherein said enhancement agent is xenon. A method of determining at least a selected one of partition coefficient and flow value of each of a plurality of voxels of a region of interest of a subject, the method comprising: at a commencement time, starting introducing an enhancement agent into the subject; measuring a concentration of the enhancement agent in a fluid of the subject; generating a plurality of image representations each at a preselected sampling interval subsequent to the commencement time, each image representation being defined by a plurality of pixels corresponding to voxels of the region of interest and including a pixel value for each pixel, which pixel value is indicative of an enhancement agent concentration in the corresponding voxel; and using said pixel values and enhancement agent concentrations to obtain corresponding values of the selected one of partition coefficient and flow value; characterised by computing parameters representative of an absorption curve indicative of enhancement agent absorption by the fluid over time projected from the measured enhancement agent concentrations; and accessing a look-up table array including a plurality of look-up tables, each look-up table corresponding to a precalculated absorption curve and being addressable by pixel values generated at preselected sampling intervals, with at least the computed parameters to select the look-up table which most closely corresponds to the projected absorption curve and pixel values of image representations generated at the preselected sampling intervals to retrieve the selected one of precalculated partition coefficient and flow value relating to the voxels corresponding to said pixel values from the selected look-up table. A method according to Claim 12 further including prior to the commencement time: calculating corresponding partition coefficient flow rate values, absorption curves, sampling intervals, and pixel values in accordance with a Kety equation; storing the precalculated partition coefficient and flow values in the look-up table array; and defining at least the precalculated absorption curves and pixel values as addresses to index the look-up table array. A method according to Claim 13 further including conducting the image representation generating step at the preselected sampling intervals. A method according to any one of Claims 12 to 14 further including interpolating the retrieved partition coefficient and flow values from the plurality of look-up tables which most closely correspond to the projected absorption curve. A method according to any one of Claims 12 to 15 wherein the look-up table array is addressable by preselected pixel values and further including interpolating partition coefficient and fluid flow values which are retrieved by addressing the look-up table array with preselected pixel values that are close to the pixel values from the generated image representations. A method according to any one of Claims 12 to 16 further including the step of displaying man-readable images defined by pixels that correspond to the voxels of the region, a first image representing the partition coefficients and a second image representing flow values. A method according to any one of Claims 12 to 17 including imaging the subject prior to the commencement time to create a reference image and subtracting the reference image from each subsequent image such that the stored pixel values are indicative of a difference image attributable substantially solely to the presence of the enhancement agent in each voxel. A method according to any one of Claims 12 to 18 wherein said fluid is blood. A method according to any one of Claims 12 to 19 wherein the enhancement agent is xenon. A method according to any one of Claims 12 to 20 wherein the image representation generating step includes radiographically examining the subject.

Inventors: LINDSTROM, WALTER W.
Applicants: PICKER INTERNATIONAL, INC.
Application year: 1987
Main classification: G06F15
IPC codes: A61B6, G06T5, A61B5