Combined X-ray and nuclear imaging

Abstract

The invention relates to a combined imaging detector (110) for the detection of x-ray and gamma quanta. The combined imaging detector (110) is adapted for simultaneous detection of gamma and x-ray quanta. The combined imaging detector (110) includes an x-ray anti-scatter grid (111), a layer of x-ray scintillator elements (112), a first photodetector array (113), a layer of gamma scintillator elements (114), and a second photodetector array (115) that are arranged in a stacked configuration along a radiation-receiving direction (116). The x-ray anti-scatter grid (111) comprises a plurality of septa (117.sub.A, B, C) that define a plurality of apertures (118) which are configured to collimate both x-ray quanta and gamma quanta received from the radiation receiving direction (116) such that received gamma quanta are collimated only by the x-ray anti-scatter grid (111). The use of the x-ray anti-scatter grid as a collimator for received gamma quanta results in a significantly lighter combined imaging detector.

Claims

1. An imaging detector for detecting gamma quanta and x-ray quanta, comprising: an x-ray anti-scatter grid; a layer of x-ray scintillator elements configured to generate x-ray scintillation light signals in response to the x-ray quanta; a first photodetector array configured to convert the x-ray scintillation light signals into a first group of electrical signals; a layer of gamma scintillator elements configured to generate gamma scintillation light signals in response to the gamma quanta; a second photodetector array configured to convert the gamma scintillation light signals into a second group of electrical signals; wherein the x-ray anti-scatter grid, the layer of x-ray scintillator elements, the first photodetector array, the layer of gamma scintillator elements, and the second photodetector array are arranged in a predetermined stacked configuration along a radiation receiving direction; wherein the x-ray anti-scatter grid comprises a plurality of septa that define a plurality of apertures which are configured to collimate the x-ray quanta and the gamma quanta received from the radiation receiving direction, such that the gamma quanta are collimated by the x-ray anti-scatter grid only; and a reconstruction unit configured to generate a reconstructed nuclear scintigraphy image from a measured nuclear scintigraphy image, which is based on the second group of electrical signals, by subtracting a reference pixel value from a pixel value in the measured nuclear scintigraphy image.

2. The imaging detector of claim 1, wherein the plurality of septa comprises a first set of septa that extend along a first direction, and a second set of septa that extend along a second direction; wherein the first set of septa and the second set of septa are arranged in a stacked configuration along the radiation receiving direction, such that the first direction and the second direction are mutually transverse.

3. The imaging detector of claim 1, wherein each septum in the plurality of septa has a depth along the radiation receiving direction that is less than or equal to approximately 5 millimeters, and/or each septum in the plurality of septa has a thickness in a direction perpendicular to the radiation receiving direction that is less than or equal to approximately 100 microns.

4. The imaging detector of claim 1, wherein the reconstruction unit is further configured to generate the reconstructed nuclear image based on a point spread function model that represents a distribution of the gamma scintillation light signals for each aperture of the x-ray anti-scatter grid.

5. The imaging detector of claim 4, wherein the point spread function model includes, in a direction perpendicular to the radiation receiving direction, a central lobe portion and a non-zero portion beyond the central lobe portion; and wherein the reconstructed nuclear image is generated by subtracting the non-zero portion from the point spread function.

6. The imaging detector of claim 4, wherein the reconstructed nuclear image is a SPECT image corresponding to a radiotracer distribution in a volume of interest, and wherein the reconstruction unit is further configured to generate the reconstructed nuclear image using an iterative reconstruction algorithm that includes: estimating the radiotracer distribution in the volume of interest; projecting the estimated radiotracer distribution onto the combined imaging detector to provide a projected estimated scintillation light distribution; receiving, based on the second group of electrical signals, distribution of the measured scintillation light signals corresponding to the radiotracer distribution; comparing the distribution of the measured scintillation light signals with the projected estimated scintillation light distribution; and updating the estimated radiotracer distribution in the volume of interest based on the comparing; wherein the point spread function model is applied when projecting the estimated radiotracer distribution.

7. The imaging detector of claim 1, wherein the reconstructed nuclear scintigraphy image corresponds to an object having a boundary; and wherein the reference pixel value is at a point beyond the boundary in the measured nuclear scintigraphy image.

8. The imaging detector of claim 1, further comprising a field of view for detected gamma quanta; wherein the reference pixel value is in the measured nuclear scintigraphy image at a point within approximately 10 pixels of a pixel at the edge of the field of view.

9. The imaging detector of claim 1, wherein the subtracting is performed in the frequency domain.

10. A C-arm comprising the imaging detector of claim 1.

11. A method for detecting gamma quanta and x-ray quanta, comprising: providing an x-ray anti-scatter grid; generating x-ray scintillation light signals by a layer of x-ray scintillator elements in response to the x-ray quanta; converting the x-ray scintillation light signals by a first photodetector array into a first group of electrical signals; generating gamma scintillation light signals by a layer of gamma scintillator elements in response to the gamma quanta; converting the gamma scintillation light signals by a second photodetector array into a second group of electrical signals; wherein the x-ray anti-scatter grid, the layer of x-ray scintillator elements, the first photodetector array, the layer of gamma scintillator elements, and the second photodetector array are arranged in a predetermined stacked configuration along a radiation receiving direction; wherein the x-ray anti-scatter grid comprises a plurality of septa that define a plurality of apertures which are configured to collimate the x-ray quanta and the gamma quanta received from the radiation receiving direction, such that received gamma quanta are collimated by the x-ray anti-scatter grid only; and generating a reconstructed nuclear scintigraphy image from a measured nuclear scintigraphy image, which is generated based on the second group of electrical signals, by subtracting a reference pixel value from a pixel value in the measured nuclear scintigraphy image.

12. A non-transitory computer-readable medium having one or more executable instructions stored thereon, which, when executed by a processor, cause the processor to perform a method for detecting gamma quanta and x-ray quanta, the method comprising: providing an x-ray anti-scatter grid; generating x-ray scintillation light signals by a layer of x-ray scintillator elements in response to the x-ray quanta; converting the x-ray scintillation light signals by a first photodetector array into a first group of electrical signals; generating gamma scintillation light signals by a layer of gamma scintillator elements in response to the gamma quanta; converting the gamma scintillation light signals by a second photodetector array into a second group of electrical signals; wherein the x-ray anti-scatter grid, the layer of x-ray scintillator elements, the first photodetector array, the layer of gamma scintillator elements, and the second photodetector array are arranged in a predetermined stacked configuration along a radiation receiving direction; wherein the x-ray anti-scatter grid comprises a plurality of septa that define a plurality of apertures which are configured to collimate the x-ray quanta and the gamma quanta received from the radiation receiving direction, such that received gamma quanta are collimated by the x-ray anti-scatter grid only; and generating a reconstructed nuclear scintigraphy image from a measured nuclear scintigraphy image, which is generated based on the second group of electrical signals, by subtracting a reference pixel value from a pixel value in the measured nuclear scintigraphy image.

Description

BRIEF DESCRIPTION OF THE FIGURES

(1) FIG. 1 illustrates a first embodiment of a combined imaging detector 110 that is adapted for simultaneous detection of gamma and x-ray quanta.

(2) FIG. 2 illustrates a portion of a 2D x-ray anti-scatter grid 211 that is formed from two crossed 1D anti-scatter grids 217.sub.x and 217.sub.y.

(3) FIG. 3 illustrates a combined imaging detector 310 in combination with a reconstruction unit 321 that is configured to generate a reconstructed nuclear image.

(4) FIG. 4 illustrates three simulated point spread functions in arbitrary units for a 1D x-ray anti-scatter grid for three different source-ASG ranges, z.

(5) FIG. 5 illustrates simulated point spread functions in arbitrary units for a 2D x-ray anti-scatter grid formed from two crossed 1D anti-scatter grids (A) for a source-ASG range, z=100 mm (B) for a source-ASG range, z=150 mm.

(6) FIG. 6 illustrates various method steps that may be used in combination with the combined imaging detector during an iterative reconstruction process to generate a SPECT image.

DETAILED DESCRIPTION OF THE INVENTION

(7) As described above, the present invention provides a system for providing both a nuclear image and an x-ray image of a region of interest. Thereto a combined imaging detector is provided.

(8) FIG. 1 illustrates a first embodiment of a combined imaging detector 110 that is adapted for simultaneous detection of gamma and x-ray quanta. The combined imaging detector includes an x-ray anti-scatter grid 111, a layer of x-ray scintillator elements 112 for generation of x-ray scintillation light signals in response to detected x-ray quanta, and a first photodetector array 113 configured to receive the x-ray scintillation light signals and to convert said x-ray scintillation light signals into a first group of electrical signals. The first group of electrical signals may be output by means of x-ray data bus 119. The combined imaging detector further includes a layer of gamma scintillator elements 114 for generation of gamma scintillation light signals in response to detected gamma quanta, and a second photodetector array 115 configured to receive the gamma scintillation light signals and to convert said gamma scintillation light signals into a second group of electrical signals. The second group of electrical signals may be output by means of gamma data bus 120. Moreover, the x-ray anti-scatter grid 111, the layer of x-ray scintillator elements 112, the first photodetector array 113, the layer of gamma scintillator elements 114, and the second photodetector array 115 are arranged in a stacked configuration along a radiation-receiving direction 116. Combined detector 110 is thus configured to receive gamma and x-ray radiation from radiation-receiving direction 116. Radiation-receiving direction 116 is perpendicular to a plane comprising any of items 111, 112, 113, 114, 115. The term perpendicular as used throughout this specification is intended to mean within approximately 10 degrees, or within 5 degrees of exactly perpendicular. X-ray anti-scatter grid 111 comprises a plurality of septa 117.sub.A, B, C that define a plurality of apertures 118 which are configured to collimate both x-ray quanta and gamma quanta received from radiation receiving direction 116 such that received gamma quanta are collimated only by the x-ray anti-scatter grid 111. As described above, the combined detector benefits from a significant weight reduction owing to the use of x-ray anti-scatter grid 111 as an anti-scatter grid for the x-ray quanta and as a collimator for the gamma quanta. Furthermore, the combined detector provides simultaneous detection of x-ray and gamma quanta, and the resulting x-ray and nuclear images are self-registered.

(9) As described above, an x-ray anti-scatter grid differs from a gamma camera collimator primarily in that the septa of an x-ray anti-scatter grid have a much shorter depth than the septa of a gamma camera collimator along radiation receiving direction 116. In the examples given above the depth of the septa along radiation receiving direction 116 for an x-ray anti-scatter grid is typically approximately 1/10.sup.th of that for a gamma camera collimator. This reduction alone gives rise to a corresponding weight reduction in the combined detector. Thus, with reference to x-ray anti-scatter grid 111 in FIG. 1, preferably each septum 117.sub.A,B,C in the plurality of septa has a depth z.sub.d1 along the radiation receiving direction 116 that is less than or equal to 5, or 4, or 3, or 2 millimeters.

(10) Another key difference between an x-ray anti-scatter grid and a gamma camera collimator is that the thickness of each septum in a direction perpendicular to the radiation receiving direction is smaller for an x-ray anti-scatter grid. Typically the septum thickness for an x-ray anti-scatter grid is approximately 40 microns. By contrast a gamma camera collimator has a septum thickness of between 150 microns and 1700 microns for the detection of 140 keV gamma quanta. Thus, with reference to x-ray anti-scatter grid 111 in FIG. 1, preferably each septum 117.sub.A,B,C in the plurality of septa has a thickness in a direction perpendicular to the radiation receiving direction 116 that is less than or equal to 100 microns, or less than or equal to 75 microns, or less than or equal to 50 microns, or less than or equal to 40 microns. These constraints on the septa thickness may be used as an alternative, or in addition to the above depth constraint to define x-ray anti-scatter grid 111.

(11) Another key difference between an x-ray anti-scatter grid and a gamma camera collimator is that an x-ray anti-scatter grid conventionally includes only a 1D array of septa. Such slits that are provided by the septa adequately reduce scatter. By contrast, a collimator for gamma imaging is typically formed from a 2D array of septa. As is described in more detail below, either a 1D x-ray anti-scatter grid that comprises a first set of septa that extend only along a first direction (x), or a 2D x-ray anti-scatter grid may be used as x-ray anti-scatter grid 111.

(12) Further notable differences between an x-ray anti-scatter grid and a gamma camera collimator include i) that the septum pitch of an x-ray anti-scatter grid is typically smaller than the septum pitch in a gamma camera collimator and ii) that the apertures defined by the septa of an x-ray anti-scatter grid are conventionally in the form of slits whereas the apertures defined by the septa of a gamma collimator are typically circular or hexagonal and thus have rotational symmetry around radiation receiving direction 116 with an order exceeding 2, or are circularly symmetric. Exemplary x-ray anti-scatter grids include the N44 r15 grid; i.e. 44 line pair per cm, or 227 micron pitch, with an interspace ratio of 15, and the N24 r12; i.e. 24 line pair per cm, or 417 micron pitch, with an interspace ration of 12. Thus the pitch of an x-ray anti-scatter grid may be less than or equal to 1 mm; and by contrast the pitch of a gamma camera collimator is conventionally greater than 1 mm.

(13) The use of x-ray anti-scatter grid 111 as a collimator for received gamma quanta essentially provides slightly poorer collimation of received gamma photons as compared to a gamma camera collimator, and thus results in a slight degradation in nuclear image quality. In essence the use of the x-ray anti-scatter grid as a collimator for the gamma camera results in a slightly elevated background signal in the nuclear image due to increased septal penetration from obliquely-incident gamma quanta. However the inventors have found that this slight degradation in nuclear image quality is acceptable for many nuclear imaging applications.

(14) The performance of a gamma camera collimator is conventionally described by a point spread function, i.e. a PSF. A PSF describes a detector's response to gamma quanta originating from a point source at a predetermined source-collimator range. This is typically represented by a graphical plot of the detected gamma scintillation light as a function of distance from an axis passing through the source and the centre of the collimator, i.e. off-axis distance. The detected gamma scintillation light is recorded in a plane that is perpendicular to the axis. In such a plot, nuclear image resolution is determined in-part by the width of the PSF.

(15) Such a plot can be generated for the response of an x-ray anti-scatter grid, i.e. an ASG, to gamma photons. FIG. 4 illustrates three simulated point spread functions in arbitrary units for a 1D x-ray anti-scatter grid for three different source-ASG ranges, z. The simulation parameters were as follows: 44 line pair per cm anti-scatter grid with an interspace ratio of 15 (code: N44 r15); lead ASG strip depth along the radiation receiving direction 2.87 mm; lead ASG strip width perpendicular to the radiation receiving direction 36 m; i.e. a typical thickness for low line pair grids; septa pitch 216 microns, modified from the typical value of 227 microns to simplify simulation. An additional plot labelled Str. Abs. in FIG. 4 represents the simulated point spread function for infinite septal absorption for a source-ASG range z=100 mm.

(16) As illustrated in FIG. 4; the three PSF curves reveal a background signal in the order of 20-30% outside the Gaussian-shaped signal. The low-frequency ripples seen on all curves are due to single, twice, triple, etc. penetration of septa. The simulated spatial resolution for this x-ray anti-scatter grid when used as a SPECT collimator for 140 keV radiation (Tc-99m) is 7.34 mm, 9.06 mm, 13.38 mm FWHM for a corresponding source-ASG distance of 80 mm, 100 mm, 150 mm. Maximum penetration of the thin septa was 51%. As noted above, this performance which is achieved with only a 1D x-ray anti-scatter grid, is considered adequate for many nuclear imaging applications. Thus, whilst FIG. 1 illustrates x-ray anti-scatter grid 111 as a 2D array of septa, x-ray anti scatter grid 111 may alternatively be provided by a 1D x-ray anti-scatter grid (211) that comprises a first set of septa that extend only along a first direction (x).

(17) The simulated point spread functions illustrated in FIG. 4 were for a 1D x-ray anti-scatter grid. Similar PSFs are also plotted for a 2D x-ray anti-scatter grid formed from two crossed 1D anti-scatter grids such as anti-scatter grid 211 illustrated in FIG. 2, in FIG. 5. FIG. 2 illustrates a portion of a 2D x-ray anti-scatter grid 211 that is formed from two crossed 1D anti-scatter grids 217.sub.x and 217.sub.y. In FIG. 2, a first set of septa 217.sub.x extend along a first direction, x, and a second set of septa 217.sub.y extend along a second direction, y. Moreover the first set of septa 217.sub.x and the second set of septa 217.sub.y are arranged in a stacked configuration along radiation receiving direction 216 such that the first direction, x, and the second direction, y, are mutually transverse. Each septum in FIG. 2 has a depth z.sub.d2 along radiation receiving direction 216. FIG. 5 illustrates simulated point spread functions in arbitrary units for a 2D x-ray anti-scatter grid formed from two crossed 1D anti-scatter grids (A) for a source-ASG range, z=100 mm (B) for a source-ASG range, z=150 mm. One quadrant of the complete PSF for the gamma scintillation light distribution in a plane perpendicular to the radiation receiving direction is illustrated in FIG. 5. The x and y axes of FIG. 5 represent in millimeters the off-axis distance in orthogonal planes x and y, in the range 0 to 10 mm. The off-axis distance, x, in FIG. 5 labelled x-DIST./mm thus corresponds to the off-axis distance in FIG. 4 labelled DIST./mm. Except for the use of crossed ASG, the simulation parameters were the same as those in FIG. 4. As compared to FIG. 4, in FIG. 5 it can be seen that the magnitude of the background signal is reduced from the 20-30% in FIG. 4 to 5.5% and 6.2% in FIG. 5A and FIG. 5B respectively due to use of the crossed ASG. Thus, the crossed ASG offers improved spatial resolution.

(18) FIG. 5A and FIG. 5B, each representing one quadrant of the whole planar PSF, indicate that without any additional processing, the projection images from point sources results in both a near-constant background signal overlaid with the typical 4-star pattern for rays parallel to the grid axes. The additional low-frequency ripple pattern in the PSF will largely average-out for typical SPECT image resolution during a subsequent nuclear image reconstruction process.

(19) Thus, as with the simulated 1D x-ray anti-scatter grid simulated in FIG. 4, FIG. 5 illustrates that adequate performance for many nuclear imaging applications can also be obtained from a gamma camera in which received gamma quanta are collimated only by a 2D x-ray anti-scatter grid that is formed from two crossed 1D anti-scatter grids. With reference to anti-scatter grid 211 in FIG. 2, preferably each septum 217.sub.x, 217.sub.y in the plurality of septa has a depth z.sub.d2 along radiation receiving direction 216 that is less than or equal to 5, or 4, or 3, or 2 millimeters. Alternatively, or additionally, each septum 217.sub.x, y in the plurality of septa may have a thickness in a direction perpendicular to the radiation receiving direction 216 that is less than or equal to 100 microns, or less than or equal to 75 microns, or less than or equal to 50 microns, or less than or equal to 40 microns. Moreover, the septa of the two crossed 1D x-ray anti-scatter grids in FIG. 2 can be combined in a single layer to provide a 2D array of septa in a single layer. The performance of such a single layer is expected to be much better than the performance of the 1D and the 2D x-ray anti-scatter grids in FIG. 4 and FIG. 5 in terms of penetration through the septa.

(20) Returning to FIG. 1; additional processing of the second electrical signals from the combined imaging detector, specifically those output at gamma data bus 120, may also be used based on the principles outlined above to further improve the quality of a subsequently-reconstructed nuclear image. Such processing is conventionally performed during image reconstruction by an image reconstruction unit.

(21) FIG. 3 illustrates a combined imaging detector 310 in combination with a reconstruction unit 321 that is configured to generate a reconstructed nuclear image. Combined imaging detector 310 corresponds to item 110 in FIG. 1 and is likewise configured to receive gamma and x-ray radiation from radiation-receiving direction 316. The nuclear image is generated based on the second group of electrical signals that are received from the second photodetector array by means of gamma bus 320. The image may subsequently be displayed by means of display device 322 that is optionally in communication with reconstruction unit 321.

(22) The form in which the PSF is accounted-for the image reconstruction process may depend to some extent on the type of nuclear image that is acquired. The image may for example be a scintigraphy image; i.e. an image acquired from a gamma detector in a fixed position, or a tomographic, i.e. a SPECT image, indicative of a cross section of an imaged object. In the case of a scintigraphy image, i.e. an image generated during a planar image acquisition mode, the non-zero portion of the PSF, i.e. a background signal, can be subtracted from a measured image. By contrast, in the case of a SPECT image it is preferable to take account of the PSF shape during the image reconstruction process.

(23) Thereto, in one configuration, reconstruction unit 321 is configured to generate a reconstructed nuclear image based on a point spread function model that represents the gamma scintillation light signal distribution for each aperture 118, 218 of x-ray anti-scatter grid 111, 211. The gamma scintillation light signal distribution is modelled for a direction or a plane that is perpendicular to radiation receiving direction 316. By taking account of the PSF during the nuclear image reconstruction process, the resulting image quality, specifically the spatial resolution can be improved. It is noted that, conventionally, when a PSF is used in a nuclear image reconstruction process it is the PSF of a SPECT collimator that is used and not the PSF of an x-ray anti-scatter grid.

(24) Optionally the point spread function model used in the reconstruction process may include, in a direction perpendicular to radiation receiving direction 316, a central lobe portion and a non-zero portion beyond the central lobe. The central lobe portion typically has a peak value in the centre of the aperture. Examples of such a lobe portion, and a non-zero portion, or background portion are illustrated in FIG. 4 between distance=0 and approximately 7 mm, and between distance=approximately 7 and 20 mm. Although FIG. 4 only plots the PSF for the positive x-axis, the central lobe portion is inherent when it is considered that the gamma scintillation light signal distribution in the negative x-axis is a mirrored version of that plotted for the positive x-axis. Similar PSFs can also be plotted for the orthogonal direction to the distance axis in FIG. 4. Moreover, in the reconstruction process, reconstruction unit 321 may be further configured to generate the reconstructed nuclear image by subtracting the non-zero portion from the point spread function. By accounting for the non-zero portion, i.e. the background portion of the PSF in this way, an improved nuclear image can be generated. As shown in FIG. 4 and FIG. 5, typically this non-zero portion can be approximately 5% to 30% of the peak value of the PSF. Conventionally when a PSF is used in an image reconstruction process it is modelled as only a central lobe similar to the straight line in FIG. 4; i.e. there is no non-zero portion in the PSF model beyond the central lobe at all.

(25) Moreover, in SPECT acquisition mode, the PSF may be applied in each loop of an iterative reconstruction process. Each time the loop is executed the noise in the resulting image is reduced. The PSF is preferably applied on the forward projected estimate of the spatial distribution of the radiotracer onto the detector.

(26) FIG. 6 illustrates various method steps that may be used in combination with the combined imaging detector during an iterative reconstruction process to generate a SPECT image. Thus, with reference to FIG. 6, when the reconstructed nuclear image is a SPECT image corresponding to an actual radiotracer distribution in a volume of interest, the reconstruction unit may be configured to generate a reconstructed nuclear image using an iterative reconstruction algorithm that includes the steps of:

(27) providing (M1) an estimate of the radiotracer distribution in the volume of interest;

(28) projecting (M2) the estimated radiotracer distribution onto the combined imaging detector to provide a projected estimated scintillation light distribution;

(29) receiving (M3), based on the second group of electrical signals, signals indicative of a measured scintillation light distribution corresponding to the actual radiotracer distribution;

(30) comparing (M4) the measured scintillation light distribution with the projected estimated scintillation light distribution; and

(31) updating (M5) the estimate of the radiotracer distribution in the volume of interest based on the comparing step;

(32) wherein the point spread function model is applied during the step of projecting (M2) the estimated radiotracer distribution.

(33) With reference to FIG. 1, the estimated radiotracer distribution is typically projected onto a plane that includes second photodetector array 115.

(34) Thus, whilst the performance of the x-ray ASG provides adequate collimation of the gamma quanta for some nuclear imaging applications in the absence of the above PSF model, further improved performance, for example that suited to diagnostic SPECT imaging, can be obtained by further taking account of the performance of the ASG through the PSF model.

(35) Furthermore, a computer program product comprising instructions for controlling the combined imaging detector described above, is disclosed. The instructions, when executed on a processor cause the processor to generate a reconstructed SPECT nuclear image corresponding to an actual radiotracer distribution in a volume of interest using an iterative reconstruction algorithm based on the second group of electrical signals and based on a point spread function model that represents the gamma scintillation light signals for each septum in the plurality of septa, comprise the steps of:

(36) providing (M1) an estimate of the radiotracer distribution in the volume of interest;

(37) projecting (M2) the estimated radiotracer distribution onto the combined imaging detector to provide a projected estimated scintillation light distribution;

(38) receiving (M3), based on the second group of electrical signals, signals indicative of a measured scintillation light distribution corresponding to the actual radiotracer distribution;

(39) comparing (M4) the measured scintillation light distribution with the projected estimated scintillation light distribution; and

(40) updating (M5) the estimate of the radiotracer distribution in the volume of interest based on the comparing step;

(41) wherein the point spread function model is applied during the step of projecting (M2) the estimated radiotracer distribution.

(42) The computer program product may be provided by dedicated hardware as well as hardware capable of executing software in association with appropriate software. When provided by a processor, the functions can be provided by a single dedicated processor, by a single shared processor, or by a plurality of individual processors, some of which can be shared. Moreover, explicit use of the term processor or controller should not be construed to refer exclusively to hardware capable of executing software, and can implicitly include, without limitation, digital signal processor DSP hardware, read only memory ROM for storing software, random access memory RAM, non-volatile storage, etc. Furthermore, embodiments of the present invention can take the form of a computer program product accessible from a computer-usable or computer-readable storage medium providing program code for use by or in connection with a computer or any instruction execution system. For the purposes of this description, a computer-usable or computer readable storage medium can be any apparatus that may include, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device. The medium can be an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, or apparatus or device, or a propagation medium. Examples of a computer-readable medium include a semiconductor or solid state memory, magnetic tape, a removable computer diskette, a random access memory RAM, a read-only memory ROM, a rigid magnetic disk and an optical disk. Current examples of optical disks include compact disk read only memory CD-ROM, compact disk read/write CD-R/W, Blu-Ray and DVD.

(43) In an alternative configuration, reconstruction unit 321 in FIG. 3 may be configured to generate a reconstructed nuclear scintigraphy image from a measured nuclear scintigraphy image. In this configuration the measured nuclear scintigraphy image comprises a plurality of pixel values that are generated based on the second group of electrical signals. In this configuration reconstruction unit 321 is further configured to generate the reconstructed nuclear scintigraphy image by subtracting a background pixel value, i.e. a reference pixel value, from the pixel values in the measured nuclear scintigraphy image. By subtracting the background pixel value, an improved nuclear image can be obtained. The subtraction may be performed either in the spatial or frequency domain. Preferably the reconstructed nuclear scintigraphy image corresponds to an object having a boundary, and the background pixel value is the value of a pixel in the measured nuclear scintigraphy image at a point beyond the boundary. Alternatively the background pixel value, i.e. the reference pixel value, that is subtracted may be an estimated value corresponding to a modelled or expected value, or the value of a pixel at or near the edge of the field of view of the gamma camera, for example within 20 or 10 or 5 pixels of the edge of the field of view. Again, when the reconstructed nuclear scintigraphy image corresponds to an object having a boundary, the background pixel value that is subtracted is, furthermore, preferably the value of a pixel in the measured nuclear scintigraphy image at a point beyond the boundary; i.e. the subtracted pixel position lies between the boundary and the edge of the field of view. Clearly the non-negativity signal constraint should be preserved. Other techniques such as 2D de-convolution filtering using an averaged, i.e. locationand depth-independentPSF may also be applied to improve compensation for directional star pattern artefacts outlined above. The image reconstruction process for a nuclear scintigraphy image can therefore be a more straightforward process than that for a SPECT image.

(44) Various example implementations are now enumerated below:

Example 1

(45) Combined imaging detector (110) for simultaneous detection of gamma and x-ray quanta; the combined imaging detector comprising:

(46) an x-ray anti-scatter grid (111);

(47) a layer of x-ray scintillator elements (112) for generation of x-ray scintillation light signals in response to detected x-ray quanta;

(48) a first photodetector array (113) configured to receive the x-ray scintillation light signals and to convert said x-ray scintillation light signals into a first group of electrical signals;

(49) a layer of gamma scintillator elements (114) for generation of gamma scintillation light signals in response to detected gamma quanta;

(50) a second photodetector array (115) configured to receive the gamma scintillation light signals and to convert said gamma scintillation light signals into a second group of electrical signals;

(51) wherein the x-ray anti-scatter grid (111), the layer of x-ray scintillator elements (112), the first photodetector array (113), the layer of gamma scintillator elements (114), and the second photodetector array (115) are arranged in a stacked configuration along a radiation-receiving direction (116);

(52) wherein the x-ray anti-scatter grid (111) comprises a plurality of septa (117A, B, C) that define a plurality of apertures (118) which are configured to collimate both x-ray quanta and gamma quanta received from the radiation receiving direction (116) such that received gamma quanta are collimated only by the x-ray anti-scatter grid (111).

Example 2

(53) The combined imaging detector of Example 1 wherein the x-ray anti-scatter grid (211) comprises a first set of septa (217x) that extend along a first direction (x), and a second set of septa (217y) that extend along a second direction (y);

(54) wherein the first set of septa (217x) and the second set of septa (217y) are arranged in a stacked configuration along the radiation receiving direction (116, 216) such that the first direction (x) and the second direction (y) are mutually transverse.

Example 3

(55) The combined imaging detector of any one of Examples 1-2 wherein i) each septum (117A,B,C, 217x, 217y) in the plurality of septa has a depth (zd1, zd2) along the radiation receiving direction (116, 216) and the depth (zd1, zd2) is less than or equal to 5 millimeters and/or ii) each septum (117A,B,C, 217x, 217y) in the plurality of septa has a thickness (t2) in a direction perpendicular to the radiation receiving direction (116, 216) and the thickness (t2) is less than or equal to 100 microns.

Example 4

(56) The combined imaging detector of any one of Examples 1-3 further comprising a reconstruction unit (321) configured to generate a reconstructed nuclear image based on the second group of electrical signals;

(57) wherein the reconstruction unit (321) is configured to generate the reconstructed nuclear image based on a point spread function model that represents the gamma scintillation light signal distribution for each aperture (118, 218) of the x-ray anti-scatter grid (111, 211).

Example 5

(58) The combined imaging detector of Example 4 wherein the point spread function model includes, in a direction perpendicular to the radiation receiving direction, a central lobe portion and a non-zero portion beyond the central lobe; and wherein the reconstruction unit is further configured to generate the reconstructed nuclear image by subtracting the non-zero portion from the point spread function.

Example 6

(59) The combined imaging detector of Example 4 or Example 5 wherein the reconstructed nuclear image is a SPECT image corresponding to an actual radiotracer distribution in a volume of interest and wherein the reconstruction unit is configured to generate the reconstructed nuclear image using an iterative reconstruction algorithm that includes the steps of:

(60) providing (M1) an estimate of the radiotracer distribution in the volume of interest;

(61) projecting (M2) the estimated radiotracer distribution onto the combined imaging detector to provide a projected estimated scintillation light distribution;

(62) receiving (M3), based on the second group of electrical signals, signals indicative of a measured scintillation light distribution corresponding to the actual radiotracer distribution;

(63) comparing (M4) the measured scintillation light distribution with the projected estimated scintillation light distribution; and

(64) updating (M5) the estimate of the radiotracer distribution in the volume of interest based on the comparing step;

(65) wherein the point spread function model is applied during the step of projecting the estimated radiotracer distribution.

Example 7

(66) The combined imaging detector of any one of Examples 1-3 further comprising a reconstruction unit (321) configured to generate a reconstructed nuclear scintigraphy image from a measured nuclear scintigraphy image comprising a plurality of pixel values, the measured nuclear scintigraphy image being generated based on the second group of electrical signals;

(67) wherein the reconstruction unit (321) is further configured to generate the reconstructed nuclear scintigraphy image by subtracting a reference pixel value from the pixel values in the measured nuclear scintigraphy image.

Example 8

(68) The combined imaging detector of Example 7, wherein the reconstructed nuclear scintigraphy image corresponds to an object having a boundary; and wherein the reference pixel value is the value of a pixel in the measured nuclear scintigraphy image at a point beyond the boundary.

Example 9

(69) The combined imaging detector of Example 7 or Example 8 wherein the combined imaging detector has a field of view for detected gamma quanta; and wherein the reference pixel value is the value of a pixel in the measured nuclear scintigraphy image at a point within 10 pixels of a pixel at the edge of the field of view.

Example 10

(70) The combined imaging detector of any one of Examples 7-9 wherein the subtracting is performed in the frequency domain.

Example 11

(71) C-arm comprising the combined imaging detector of any one of Examples 1-10.

Example 12

(72) Computer program product comprising instructions for controlling the combined imaging detector of any one of Examples 1-3; which instructions, when executed on a processor cause the processor to either:

(73) i) generate a reconstructed SPECT nuclear image corresponding to an actual radiotracer distribution in a volume of interest using an iterative reconstruction algorithm based on the second group of electrical signals and based on a point spread function model that represents the gamma scintillation light signals for each septum in the plurality of septa, wherein the instructions comprise the steps of:

(74) providing (M1) an estimate of the radiotracer distribution in the volume of interest;

(75) projecting (M2) the estimated radiotracer distribution onto the combined imaging detector to provide a projected estimated scintillation light distribution; receiving (M3), based on the second group of electrical signals, signals indicative of a measured scintillation light distribution corresponding to the actual radiotracer distribution;

(76) comparing (M4) the measured scintillation light distribution with the projected estimated scintillation light distribution; and

(77) updating (M5) the estimate of the radiotracer distribution in the volume of interest based on the comparing step;

(78) wherein the point spread function model is applied during the step of projecting the estimated radiotracer distribution;

(79) or to

(80) ii) generate a reconstructed nuclear scintigraphy image from a measured nuclear scintigraphy image comprising a plurality of pixel values, the measured nuclear scintigraphy image being generated based on the second group of electrical signals, wherein the instructions comprise the step of:

(81) subtracting a reference pixel value from the pixel values in the measured nuclear scintigraphy image.

(82) In summary, a combined imaging detector has been described. The combined imaging detector 110 is adapted for simultaneous detection of gamma and x-ray quanta. The combined imaging detector 110 includes an x-ray anti-scatter grid 111, a layer of x-ray scintillator elements 112, a first photodetector array 113, a layer of gamma scintillator elements 114, and a second photodetector array 115 that are arranged in a stacked configuration along a radiation-receiving direction 116. The x-ray anti-scatter grid 111 comprises a plurality of septa 117A, B, c that define a plurality of apertures 118 which are configured to collimate both x-ray quanta and gamma quanta received from the radiation receiving direction 116 such that received gamma quanta are collimated only by the x-ray anti-scatter grid 111.