WO2009150565A2 - Reverse data reconstruction for optimal time sampling of counts in physiological list-mode nuclear imaging - Google Patents

Reverse data reconstruction for optimal time sampling of counts in physiological list-mode nuclear imaging Download PDF

Info

Publication number
WO2009150565A2
WO2009150565A2 PCT/IB2009/052290 IB2009052290W WO2009150565A2 WO 2009150565 A2 WO2009150565 A2 WO 2009150565A2 IB 2009052290 W IB2009052290 W IB 2009052290W WO 2009150565 A2 WO2009150565 A2 WO 2009150565A2
Authority
WO
WIPO (PCT)
Prior art keywords
image representation
set forth
events
merit
image
Prior art date
Application number
PCT/IB2009/052290
Other languages
French (fr)
Other versions
WO2009150565A3 (en
Inventor
Daniel Gagnon
Sameer Tipnis
Raymond Muzic
Original Assignee
Koninklijke Philips Electronics N.V.
University Hospitals Medical Group, Inc.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips Electronics N.V., University Hospitals Medical Group, Inc. filed Critical Koninklijke Philips Electronics N.V.
Priority to EP09762092.6A priority Critical patent/EP2288291B1/en
Priority to US12/995,685 priority patent/US8437836B2/en
Priority to CN200980121926.3A priority patent/CN102098964B/en
Priority to JP2011513080A priority patent/JP2011523076A/en
Priority to RU2011100843/14A priority patent/RU2517586C2/en
Publication of WO2009150565A2 publication Critical patent/WO2009150565A2/en
Publication of WO2009150565A3 publication Critical patent/WO2009150565A3/en

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/50Clinical applications
    • A61B6/507Clinical applications involving determination of haemodynamic parameters, e.g. perfusion CT
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computerised tomographs
    • A61B6/037Emission tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/486Diagnostic techniques involving generating temporal series of image data
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/006Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/412Dynamic

Definitions

  • the present application relates to the diagnostic imaging arts. It finds particular application in a nuclear medicine scanner that utilizes list-mode data acquisition, and will be described with particular reference thereto. It is to be understood, however, that it also finds application in other tracer perfusion studies, and is not necessarily limited to the aforementioned application.
  • the quality of the reconstructed images is dependent on the choice of the reconstruction interval. Choosing a longer reconstruction time interval, such as 90 - 360 seconds, can maximize the signal strength, that is, it can maximize the counts received. This comes at the expense of contrast. For example, in imaging the heart, the above selected time range might result in low contrast since blood in the ventricles, not yet absorbed by the tissues of the heart, may still be active. Because the early data has the highest count rate, it tends to dominate the lower count rate later data. On the other hand, if a window of 300 - 360 seconds is selected, the image will most likely have good contrast, as most of the tracer has been absorbed by the tissue of interest. There may be significant noise, however, as most of the tracer will have decayed by that time, leading to a low event count rate.
  • the present application provides a new and improved event processing method that is able to leverage list-mode data to efficiently optimize the image quality.
  • a diagnostic imaging apparatus In accordance with one aspect, a diagnostic imaging apparatus is provided.
  • a detector array including individual detectors senses photons emitted by radioactivity decay within the patient.
  • a triggering processor assigns a time stamp to received potential events.
  • An event verification processor applies verification criteria to received potential events.
  • a reconstruction processor reconstructs valid events into an image representation of the imaging region.
  • a figure of merit analyzer analyzes a reconstructed image to determine a figure of merit.
  • a method of diagnostic imaging is provided.
  • a set of data points indicative of nuclear decay events is collected and the data points are sorted according to a time that the data points were detected.
  • a reference point in time is selected.
  • An image representation is reconstructed from data points occurring before the reference point.
  • a figure of merit associated with the image representation is determined.
  • a time interval is selected and applied backwards from the reference point creating a new reference point.
  • a method of diagnostic imaging is provided. Photons emitted in radioactive decay are detected.
  • a time stamp is assigned to received potential events. Verification criteria are applied to received events.
  • Valid events are reconstructed into an image representation of the imaging region.
  • a reconstructed image is analyzed to determine a figure of merit.
  • One advantage lies in the ability to select an optimal combination of signal strength and image contrast in a nuclear tracer perfusion study.
  • Another advantage lies in the ability to iteratively add data to or remove data from a reconstructed image until a maximum image quality is achieved. Another advantage is that it is easily retrofitted into present scanners provided that the scanners have compatible data generation and data format (e.g. list mode). Still further advantages of the present invention will be appreciated to those of ordinary skill in the art upon reading and understand the following detailed description.
  • the invention may take form in various components and arrangements of components, and in various steps and arrangements of steps.
  • the drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention.
  • FIGURE 1 is a diagrammatic illustration of a nuclear imaging device in accordance with the present application
  • FIGURE 2 is a graph showing count rate and tissue perfusion as functions of time
  • FIGURE 3 is a flow diagram of an iterative determination of the most valuable data to be used in reconstruction, in accordance with the present application.
  • a diagnostic imaging device 10 includes a housing 12 and a subject support 14. Enclosed within the housing 12 is a detector array.
  • the detector array includes a plurality of individual detector elements 16. While one particular embodiment is described with reference to a positron emission tomography (PET) scanner, it is to be understood that the present application is also useful in other medical applications, such as single photon emission computed tomography (SPECT) as well as x-ray astrophysics, gamma ray telescopes, radiography, security, and industrial applications. Generally, the present application finds use in imaging x-rays, gamma rays, or charged particles with high energy and spatial resolution.
  • the array is arranged so that detector elements 16 are disposed adjacent to an imaging region 18.
  • the detector array can be a ring of detectors 16, multiple rings, one or more discrete flat or arced panels, or the like.
  • PET positron emission tomography
  • pairs of gamma rays are produced by a positron annihilation event in the imaging region and travel in approximately opposite directions. Such an event may be produced from the nuclear decay of 82 Rb.
  • These gamma rays are detected as pairs, with a slight time difference (on the order of nanoseconds or fractions thereof) between detections if one gamma ray travels farther to reach a detector than the other. Accordingly, in PET scanners, the detector arrays typically encircle the imaging region.
  • the radiopharmaceutical contains a radioactive element, such as 82 Rb, coupled to a tag molecule.
  • the tag molecule is associated with the region to be imaged, and tends to gather there through body processes. For example, rapidly multiplying cancer cells tend to expend abnormally high amounts of energy duplicating themselves.
  • the radiopharmaceutical can be linked to a molecule, such as glucose, or an analog thereof, that a cell typically metabolizes to create energy, which gathers in such regions and appear as "hot spots" in the image.
  • a tag is also useful in cardiac perfusion imaging, since the heart expends relatively large amounts of energy.
  • a triggering processor 20 monitors each detector 16 for an energy spike, e.g., integrated area under the pulse, characteristic of the energy of the gamma rays generated by the radiopharmaceutical.
  • the triggering processor 20 checks a clock 22 and stamps each detected gamma ray with a time of leading edge receipt stamp.
  • the time stamp, energy estimate and position estimation is first used by an event verification processor 24 to determine if the event data is valid, e.g., if the pair of events are coincident, have the proper energy, and the like. Accepted pairs define lines of response (LORs). Because gamma rays travel at the speed of light, if detected gamma rays arrive more than several nanoseconds apart, they probably were not generated by the same annihilation event and are usually discarded. Timing is especially important in time of flight PET (TOF-PET), as the minute difference in substantially simultaneous coincident events is used to further localize the annihilation event along the LOR. As the temporal resolution of events becomes more precise, so too does the accuracy with which an event can be localized along its LOR. LORs are stored in an event storage buffer 26. In one embodiment, the
  • LORs are stored in a list-mode format. That is, the events are stored in temporal order with time indicators periodically inserted. Alternatively, the events can be individually time stamped.
  • a reconstruction processor 28 reconstructs all or a portion of the LORs into an image representation of the subject using filtered backprojection or other appropriate reconstruction algorithms.
  • An analyzer 29 analyzes a reconstructed image to determine a figure of merit or other indicator of image quality. The analyzer 29 indexes a temporal window selector 31 to adjust the temporal window that defines a portion of the LORs that are reconstructed until the image quality is optimized or reaches a preselected level. This process is described in more detail hereinbelow.
  • the reconstruction can then be displayed for a user on a display device 30, printed, saved for later use, and the like.
  • event data is collected in a "list-mode" format. Recording the relevant properties of each detected event in a list has become a common practice in emission tomography applications and has become known as list-mode data acquisition and storage.
  • the list-mode reconstruction approach differs in several ways from binned or histogrammed-mode methods.
  • List-mode data acquisitions provide extremely high temporal resolution with full spatial resolution and allows frame durations to be determined after acquisition. Acquiring the data in list-mode format, the interaction location can be stored to a high degree of accuracy with greater efficiency than achievable with frame mode acquisition.
  • Gantry angles do not have to be binned into predefined frames, but can be recorded as the actual angle, thereby removing the impact of angular blurring with continuous acquisition.
  • the actual energy of the interaction can be recorded instead of attributing the event's energy to one of a limited number of pre-defined windows.
  • the data is stored in a list where it can be arranged and sorted by several different parameters, e.g. time of receipt, instead of a bin.
  • List-mode can also store gating signals without temporally framing of the data before this information is completely available. The result is a significant increase in the fidelity of recording the projection data with list-mode acquisition, without a tremendous increase in storage space.
  • Another advantage is the ability to identify events by the time of their occurrence, and beneficially to the present application, being able to select a window of time in which the most optimal combination of contrast and signal strength is observed.
  • the most optimal quality of a reconstructed image is achieved by a judicious choice of the reconstruction interval or window 33 within the time span of the data acquisition.
  • the count rate 32 is highest at the beginning, decreasing toward the end.
  • the contrast 34 is lowest at the beginning of the scan.
  • the LORs collected in the selectable temporal window 33 are reconstructed.
  • the start 36 of the window and/or the end 38 are selectively adjusted to optimize the balance between contrast and count rate.
  • multiple windows including overlapping windows, can be defined and an image reconstructed for each.
  • the subject is injected with the radiopharmaceutical 40.
  • the data is collected 42 and time stamped during the scan.
  • the time of the data collection can vary depending on the radioactive material used to produce the annihilation events. In the case of 82 Rb, the scan typically lasts about six minutes (360 seconds). In this time period, physically only one-sixteenth of the original amount of 82 Rb will remain, and it is expected that the count rate will have become too low to be diagnostically useful.
  • the window as defined by the initial or start point and the terminal or end point are selected 44.
  • the terminal or end point 38 is selected to be a few seconds from the end of the scan. As mentioned previously, this data is most likely to provide the highest contrast image reconstruction, as the tracer has been given ample time to profuse to the tissue of interest. Using data from the very end of the scan is a universally applicable approach, and is useful especially when little is known of the physiology of the subject (human or animal) with respect to profusion.
  • the window can be selected to be on or about that time.
  • the ratio of the myocardium activity compared to the blood stream e.g., cavity of the ventricle
  • a baseline image is reconstructed 46.
  • the baseline image is reconstructed using data from the start point until the end of the scan.
  • This baseline image includes high-contrast data counts, but if the start point is set late, the volume of the counts is low, leading to a low signal-to-noise ratio.
  • the baseline image is evaluated to determine a figure of merit 48.
  • the figure of merit is the contrast-to-noise ratio.
  • Other figures of merit are certainly contemplated, such as raw event count per unit time, and others.
  • the start point is adjusted 50, e.g., moved earlier to increase the number of counts in an effort to improve the specified figure of merit.
  • LORs in the increased region of the window are added to the LORs that were used to reconstruct the baseline image, and an updated image is reconstructed 52.
  • the figure of merit is recalculated for the new, updated image 54. If the figure of merit improves or remains stable within some selected statistical confines of the test, shown as decision block 56 in FIGURE 3, the start event may be further adjusted (e.g., moving a small time period backwards, using a binary search, or a selected optimization algorithm) and the process repeated.
  • the image that achieves the best figure of merit is used for further analysis 58, display, interpretation, or retained for future use.
  • the start and end times also be displayed to the user, as uptake times can have significant diagnostic value.
  • the end point for the window may optionally be adjusted. Data at the end of the scan may become less valuable as fewer counts occur.
  • the end point of the window may be stepped backwards in time, and the updated image 52 is reconstructed. The new figure of merit is again re-calculated 54 until it ceases to improve. Optimization of the end point can occur either before or after optimization of the start point.
  • the window 33 can be displayed to a user so the user can use it for aid in analyses, quantification, or diagnoses.
  • the optimized window 33 can also be used by the system as a parameter in computer assisted diagnosis, as the uptake times and efficiencies can be indicative of body processes and function. In one embodiment, the above-described process is automated, performed by the analyzer 29.
  • the beginning and/or end time can be set by a user with a user input 60.
  • the user could drag the start and end time indicators in a display analogous to FIGURE 2.
  • the user may be prompted to decide whether the figure of merit has improved in a particular image relative to its predecessor.
  • the analyzer 29 comes to a decision on which image is the best, it may present it and several prior and subsequent images to the user for review.
  • the process is still largely automated, but with user analysis to supplement.
  • the levels of user input and feedback desired would be selectable, and able to be toggled on or off at will.
  • the data can be compiled in a non-list-mode format. This embodiment present an additional problem of sorting events by time after they have already been binned, requiring extra processing power and time.
  • Some radiopharmaceuticals are absorbed by different tissue types within the body at different rates.
  • the optimum image for different parts of the body might have different reference points.
  • a plurality of optimized images may be displayed.
  • a series of cine images could be displayed that show the uptake and washout over time. The time window corresponding to each of the cine images can be optimized in the same way.

Abstract

A diagnostic imaging device includes detector elements (16) for detecting ƴ-rays indicative of nuclear decay events. Pairs of concurrently detected γ-rays define lines of response (LORs) which are collected, time stamped, and compiled in list-mode. In tissue perfusion studies, it is beneficial to use the data that concurrently maximizes contrast and signal-to-noise ratio in the reconstructed images. Using the list-mode data, events in an adjustable temporal window (33) are reconstructed and the reconstructed images are analyzed to determine a figure of merit based on contrast and signal-to-noise properties of the image. By iteratively adjusting the temporal window, extending its start point (36) backwards in time, and repeating the reconstructing, analyzing, and adjusting steps, an image with an optimal figure of merit is obtained.

Description

REVERSE DATA RECONSTRUCTION FOR OPTIMAL TIME SAMPLING OF COUNTS IN PHYSIOLOGICAL LIST-MODE NUCLEAR IMAGING
DESCRIPTION
The present application relates to the diagnostic imaging arts. It finds particular application in a nuclear medicine scanner that utilizes list-mode data acquisition, and will be described with particular reference thereto. It is to be understood, however, that it also finds application in other tracer perfusion studies, and is not necessarily limited to the aforementioned application.
The capture of physiological processes with nuclear dynamic imaging is a complex process requiring optimization of the sampling for the particular imaging system and an understanding of the phenomena involved in the tracer uptake. For example, cardiac perfusion studies using a rubidium-82 (82Rb) injection and PET imaging combines the challenge of dynamically imaging a high flux, short half life isotope (about 90 seconds) and an uptake process that can extract Rb from the blood stream in as short as 30 seconds or as long as 200 seconds, depending on the subject, in a typical 360 second scan. Initially, there is a large amount of tracer in the blood, but little taken up by the area of interest. This generates a large amount of data, but with poor contrast. As time progresses in the scan, more of the tracer reaches and is taken up by the area of interest and is cleared from the blood, but concurrently, signal strength wanes as time goes on due to the natural decay of the tracer. This results in good contrast, but a low amount of data. There is difficulty in finding the optimal time for imaging that harmonizes signal strength and contrast.
The quality of the reconstructed images is dependent on the choice of the reconstruction interval. Choosing a longer reconstruction time interval, such as 90 - 360 seconds, can maximize the signal strength, that is, it can maximize the counts received. This comes at the expense of contrast. For example, in imaging the heart, the above selected time range might result in low contrast since blood in the ventricles, not yet absorbed by the tissues of the heart, may still be active. Because the early data has the highest count rate, it tends to dominate the lower count rate later data. On the other hand, if a window of 300 - 360 seconds is selected, the image will most likely have good contrast, as most of the tracer has been absorbed by the tissue of interest. There may be significant noise, however, as most of the tracer will have decayed by that time, leading to a low event count rate.
The present application provides a new and improved event processing method that is able to leverage list-mode data to efficiently optimize the image quality.
In accordance with one aspect, a diagnostic imaging apparatus is provided. A detector array including individual detectors senses photons emitted by radioactivity decay within the patient. A triggering processor assigns a time stamp to received potential events. An event verification processor applies verification criteria to received potential events. A reconstruction processor reconstructs valid events into an image representation of the imaging region. A figure of merit analyzer analyzes a reconstructed image to determine a figure of merit.
In accordance with another aspect, a method of diagnostic imaging is provided. A set of data points indicative of nuclear decay events is collected and the data points are sorted according to a time that the data points were detected. A reference point in time is selected. An image representation is reconstructed from data points occurring before the reference point. A figure of merit associated with the image representation is determined. A time interval is selected and applied backwards from the reference point creating a new reference point. In accordance with another aspect, a method of diagnostic imaging is provided. Photons emitted in radioactive decay are detected. A time stamp is assigned to received potential events. Verification criteria are applied to received events. Valid events are reconstructed into an image representation of the imaging region. A reconstructed image is analyzed to determine a figure of merit.
One advantage lies in the ability to select an optimal combination of signal strength and image contrast in a nuclear tracer perfusion study.
Another advantage lies in the ability to iteratively add data to or remove data from a reconstructed image until a maximum image quality is achieved. Another advantage is that it is easily retrofitted into present scanners provided that the scanners have compatible data generation and data format (e.g. list mode). Still further advantages of the present invention will be appreciated to those of ordinary skill in the art upon reading and understand the following detailed description.
The invention may take form in various components and arrangements of components, and in various steps and arrangements of steps. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention.
FIGURE 1 is a diagrammatic illustration of a nuclear imaging device in accordance with the present application; FIGURE 2 is a graph showing count rate and tissue perfusion as functions of time;
FIGURE 3 is a flow diagram of an iterative determination of the most valuable data to be used in reconstruction, in accordance with the present application.
With reference to FIGURE 1, a diagnostic imaging device 10 includes a housing 12 and a subject support 14. Enclosed within the housing 12 is a detector array. The detector array includes a plurality of individual detector elements 16. While one particular embodiment is described with reference to a positron emission tomography (PET) scanner, it is to be understood that the present application is also useful in other medical applications, such as single photon emission computed tomography (SPECT) as well as x-ray astrophysics, gamma ray telescopes, radiography, security, and industrial applications. Generally, the present application finds use in imaging x-rays, gamma rays, or charged particles with high energy and spatial resolution. The array is arranged so that detector elements 16 are disposed adjacent to an imaging region 18. The detector array can be a ring of detectors 16, multiple rings, one or more discrete flat or arced panels, or the like. In positron emission tomography (PET), pairs of gamma rays are produced by a positron annihilation event in the imaging region and travel in approximately opposite directions. Such an event may be produced from the nuclear decay of 82Rb. These gamma rays are detected as pairs, with a slight time difference (on the order of nanoseconds or fractions thereof) between detections if one gamma ray travels farther to reach a detector than the other. Accordingly, in PET scanners, the detector arrays typically encircle the imaging region.
Before the PET scan commences, a subject is injected with a radiopharmaceutical. In one common exam, the radiopharmaceutical contains a radioactive element, such as 82Rb, coupled to a tag molecule. The tag molecule is associated with the region to be imaged, and tends to gather there through body processes. For example, rapidly multiplying cancer cells tend to expend abnormally high amounts of energy duplicating themselves. The radiopharmaceutical can be linked to a molecule, such as glucose, or an analog thereof, that a cell typically metabolizes to create energy, which gathers in such regions and appear as "hot spots" in the image. Such a tag is also useful in cardiac perfusion imaging, since the heart expends relatively large amounts of energy. Other techniques monitor tagged molecules flowing in the circulatory system. In such a technique, it is beneficial to tag a molecule that is not quickly absorbed by tissues of the body. When a gamma ray strikes the detector array, a time signal is generated. A triggering processor 20 monitors each detector 16 for an energy spike, e.g., integrated area under the pulse, characteristic of the energy of the gamma rays generated by the radiopharmaceutical. The triggering processor 20 checks a clock 22 and stamps each detected gamma ray with a time of leading edge receipt stamp. The time stamp, energy estimate and position estimation is first used by an event verification processor 24 to determine if the event data is valid, e.g., if the pair of events are coincident, have the proper energy, and the like. Accepted pairs define lines of response (LORs). Because gamma rays travel at the speed of light, if detected gamma rays arrive more than several nanoseconds apart, they probably were not generated by the same annihilation event and are usually discarded. Timing is especially important in time of flight PET (TOF-PET), as the minute difference in substantially simultaneous coincident events is used to further localize the annihilation event along the LOR. As the temporal resolution of events becomes more precise, so too does the accuracy with which an event can be localized along its LOR. LORs are stored in an event storage buffer 26. In one embodiment, the
LORs are stored in a list-mode format. That is, the events are stored in temporal order with time indicators periodically inserted. Alternatively, the events can be individually time stamped. A reconstruction processor 28 reconstructs all or a portion of the LORs into an image representation of the subject using filtered backprojection or other appropriate reconstruction algorithms. An analyzer 29 analyzes a reconstructed image to determine a figure of merit or other indicator of image quality. The analyzer 29 indexes a temporal window selector 31 to adjust the temporal window that defines a portion of the LORs that are reconstructed until the image quality is optimized or reaches a preselected level. This process is described in more detail hereinbelow. The reconstruction can then be displayed for a user on a display device 30, printed, saved for later use, and the like.
In one embodiment, event data is collected in a "list-mode" format. Recording the relevant properties of each detected event in a list has become a common practice in emission tomography applications and has become known as list-mode data acquisition and storage. The list-mode reconstruction approach differs in several ways from binned or histogrammed-mode methods. List-mode data acquisitions provide extremely high temporal resolution with full spatial resolution and allows frame durations to be determined after acquisition. Acquiring the data in list-mode format, the interaction location can be stored to a high degree of accuracy with greater efficiency than achievable with frame mode acquisition. Gantry angles do not have to be binned into predefined frames, but can be recorded as the actual angle, thereby removing the impact of angular blurring with continuous acquisition. The actual energy of the interaction can be recorded instead of attributing the event's energy to one of a limited number of pre-defined windows. When increasing the dimensionality in this way, the data is stored in a list where it can be arranged and sorted by several different parameters, e.g. time of receipt, instead of a bin. List-mode can also store gating signals without temporally framing of the data before this information is completely available. The result is a significant increase in the fidelity of recording the projection data with list-mode acquisition, without a tremendous increase in storage space. Another advantage is the ability to identify events by the time of their occurrence, and beneficially to the present application, being able to select a window of time in which the most optimal combination of contrast and signal strength is observed.
With reference to FIGURE 2, the most optimal quality of a reconstructed image is achieved by a judicious choice of the reconstruction interval or window 33 within the time span of the data acquisition. During the data acquisition, e.g., 360 seconds, the count rate 32 is highest at the beginning, decreasing toward the end. Conversely, the contrast 34 is lowest at the beginning of the scan. The LORs collected in the selectable temporal window 33 are reconstructed. The start 36 of the window and/or the end 38 are selectively adjusted to optimize the balance between contrast and count rate. Of course, multiple windows, including overlapping windows, can be defined and an image reconstructed for each.
With reference now to FIGURE 3, a flowchart describing image optimization is provided. First, the subject is injected with the radiopharmaceutical 40. The data is collected 42 and time stamped during the scan. The time of the data collection can vary depending on the radioactive material used to produce the annihilation events. In the case of 82Rb, the scan typically lasts about six minutes (360 seconds). In this time period, physically only one-sixteenth of the original amount of 82Rb will remain, and it is expected that the count rate will have become too low to be diagnostically useful.
After the data has been collected, the window as defined by the initial or start point and the terminal or end point are selected 44. In one embodiment, the terminal or end point 38 is selected to be a few seconds from the end of the scan. As mentioned previously, this data is most likely to provide the highest contrast image reconstruction, as the tracer has been given ample time to profuse to the tissue of interest. Using data from the very end of the scan is a universally applicable approach, and is useful especially when little is known of the physiology of the subject (human or animal) with respect to profusion.
If more is known, more significant points can be selected. For instance, with 82Rb in cardiac imaging, 200 seconds should be enough time for healthy human subjects. More generally, if there is a time at which it is known that the physiology of the subject should be in a steady state with respect to the variable to be measured, the window can be selected to be on or about that time. Further, there may be a parameter that can differentiate an intermediate state from the final steady state. Again in an example, in 82Rb cardiac imaging, the ratio of the myocardium activity compared to the blood stream (e.g., cavity of the ventricle) should be an indicator of the process.
Next, a baseline image is reconstructed 46. In one embodiment, the baseline image is reconstructed using data from the start point until the end of the scan. This baseline image includes high-contrast data counts, but if the start point is set late, the volume of the counts is low, leading to a low signal-to-noise ratio. The baseline image is evaluated to determine a figure of merit 48. In one embodiment, the figure of merit is the contrast-to-noise ratio. Other figures of merit are certainly contemplated, such as raw event count per unit time, and others.
Given the specified figure of merit for the baseline image, the start point is adjusted 50, e.g., moved earlier to increase the number of counts in an effort to improve the specified figure of merit. LORs in the increased region of the window are added to the LORs that were used to reconstruct the baseline image, and an updated image is reconstructed 52. Once the updated image has been reconstructed, the figure of merit is recalculated for the new, updated image 54. If the figure of merit improves or remains stable within some selected statistical confines of the test, shown as decision block 56 in FIGURE 3, the start event may be further adjusted (e.g., moving a small time period backwards, using a binary search, or a selected optimization algorithm) and the process repeated. The image that achieves the best figure of merit is used for further analysis 58, display, interpretation, or retained for future use. The start and end times also be displayed to the user, as uptake times can have significant diagnostic value.
In the same fashion, the end point for the window may optionally be adjusted. Data at the end of the scan may become less valuable as fewer counts occur. The end point of the window may be stepped backwards in time, and the updated image 52 is reconstructed. The new figure of merit is again re-calculated 54 until it ceases to improve. Optimization of the end point can occur either before or after optimization of the start point. During the described window selection process, the window 33 can be displayed to a user so the user can use it for aid in analyses, quantification, or diagnoses. The optimized window 33 can also be used by the system as a parameter in computer assisted diagnosis, as the uptake times and efficiencies can be indicative of body processes and function. In one embodiment, the above-described process is automated, performed by the analyzer 29. It is contemplated that varying levels of user input can be accepted to aid in the decision. For example, the beginning and/or end time can be set by a user with a user input 60. For example, the user could drag the start and end time indicators in a display analogous to FIGURE 2. The user may be prompted to decide whether the figure of merit has improved in a particular image relative to its predecessor. When the analyzer 29 comes to a decision on which image is the best, it may present it and several prior and subsequent images to the user for review. In this embodiment, the process is still largely automated, but with user analysis to supplement. The levels of user input and feedback desired would be selectable, and able to be toggled on or off at will.
In an alternate embodiment, the data can be compiled in a non-list-mode format. This embodiment present an additional problem of sorting events by time after they have already been binned, requiring extra processing power and time.
Some radiopharmaceuticals are absorbed by different tissue types within the body at different rates. Thus, the optimum image for different parts of the body might have different reference points. Accordingly, a plurality of optimized images may be displayed. As another alternative, a series of cine images could be displayed that show the uptake and washout over time. The time window corresponding to each of the cine images can be optimized in the same way.
The invention has been described with reference to the preferred embodiments. Modifications and alterations may occur to others upon reading and understanding the preceding detailed description. It is intended that the invention be construed as including all such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.

Claims

CLAIMSHaving thus described the preferred embodiments, the invention is now claimed to be:
1. A diagnostic imaging apparatus comprising: a detector array for receiving events from an imaging region (18); a triggering processor (20) for assigning a time stamp to the received events; a reconstruction processor (28) that reconstructs an adjustable portion of the events into an image representation; an analyzer (29) that analyzes the reconstructed image representation to determine a parameter indicative of image quality and adjusts the adjustable portion of the verified events for reconstruction into a subsequent image representation in accordance with the determined parameter.
2. The diagnostic imaging apparatus as set forth in claim 1, wherein the determined parameter is a figure of merit based on a selected parameter of the image representation.
3. The diagnostic imaging apparatus as set forth in either of claims 1 or 2, further including: an event storage buffer (26) that stores verified events in a list that is sortable by time of receipt of the event; and a temporal window selector (31) which is controlled by the analyzer (29) to adjust a temporal window (33) that defines the events which are reconstructed.
4. The diagnostic imaging apparatus as set forth in claim 3, further including: a display (30) for displaying at least the temporal window (33) to a user.
5. A method of diagnostic imaging comprising: collecting a set of data points indicative of received nuclear decay events; time stamping the received events; reconstructing an image representation from an adjustable portion of the collected data points; analyzing the image representation to determine a value of a parameter indicative of image representation quality; and adjusting the adjustable portion of the collected data points which are used in reconstructing a subsequent image representation.
6. The method as set forth in claim 5, wherein the parameter determining step includes: determining a figure of merit for the image representation.
7. The method as set forth in 6, further including: comparing the figure of merit for the subsequent image representation to the figure of merit for the previous image representation, and if the figure of merit for the subsequent image representation is greater than or equal to the figure of merit for the previous image representation, repeating further adjusting the portion of the collected data points that are used to reconstruct a next image representation.
8. The method as set forth in claim 6, wherein the figure of merit is a contrast-to-noise ratio of the image representation.
9. The method as set forth in any one of claims 5-8, further including: sorting the data points by time of detection, the portion of the data points being in an adjustable temporal window (33).
10. The method as set forth in claim 9, further including: analyzing the subsequent image representations to determine a value of the parameter indicative of the second image representation quality; in response to the values of the subsequent and previous image representation parameters indicating that the subsequent image is better than the previous image, adjusting the temporal window (33) and generating another image representation.
11. The method as set forth in claim 10, wherein the temporal window has a start point and an end point and further including: in response to the subsequent image being better than the previous image representation, moving the start point toward a beginning of the scan.
12. The method as set forth in claim 10, further including: iteratively repeating the steps of adjusting, reconstructing, and analyzing until an image representation with value of the parameter is one of optimized or exceeds a selected threshold.
13. The method as set forth in any one of claims 9-12, further including: displaying the window (33) to a user.
14. The method as set forth in any one of claims 9-13, further including: analyzing the window (33) to aid in one of a diagnosis, quantification, and a measure of a bodily process
15. A nuclear imaging system including a processor programmed to control the nuclear imaging system to perform the method of any one of claims 5-14.
16. A computer readable medium carrying software to control a processor to perform the method of any one of claims 5-14.
PCT/IB2009/052290 2008-06-13 2009-05-29 Reverse data reconstruction for optimal time sampling of counts in physiological list-mode nuclear imaging WO2009150565A2 (en)

Priority Applications (5)

Application Number Priority Date Filing Date Title
EP09762092.6A EP2288291B1 (en) 2008-06-13 2009-05-29 Reverse data reconstruction for optimal time sampling of counts in physiological list-mode nuclear imaging
US12/995,685 US8437836B2 (en) 2008-06-13 2009-05-29 Reverse data reconstruction for optimal time sampling of counts in physiological list-mode nuclear imaging
CN200980121926.3A CN102098964B (en) 2008-06-13 2009-05-29 Reverse data reconstruction for optimal time sampling of counts in physiological list-mode nuclear imaging
JP2011513080A JP2011523076A (en) 2008-06-13 2009-05-29 Inverse data reconstruction for optimal time sampling of counts in physiological list mode nuclear medicine imaging
RU2011100843/14A RU2517586C2 (en) 2008-06-13 2009-05-29 Reverse reconstruction of data for optimal time generation of count pulses in radiological physiological visualisation in list mode

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US6128908P 2008-06-13 2008-06-13
US61/061,289 2008-06-13

Publications (2)

Publication Number Publication Date
WO2009150565A2 true WO2009150565A2 (en) 2009-12-17
WO2009150565A3 WO2009150565A3 (en) 2010-06-24

Family

ID=41417195

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/IB2009/052290 WO2009150565A2 (en) 2008-06-13 2009-05-29 Reverse data reconstruction for optimal time sampling of counts in physiological list-mode nuclear imaging

Country Status (6)

Country Link
US (1) US8437836B2 (en)
EP (1) EP2288291B1 (en)
JP (2) JP2011523076A (en)
CN (1) CN102098964B (en)
RU (1) RU2517586C2 (en)
WO (1) WO2009150565A2 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2464212A (en) * 2008-10-09 2010-04-14 Siemens Medical Solutions Analyzing Medical Image Data
CN102028494A (en) * 2011-01-24 2011-04-27 海纳医信(北京)软件科技有限责任公司 Method and system for processing cerebral perfusion image sequence
CN104598356A (en) * 2014-12-23 2015-05-06 沈阳东软医疗系统有限公司 Event ordering method and device
US10133952B2 (en) 2014-12-24 2018-11-20 Shenyang Neusoft Medical Systems Co., Ltd. Computer tomography image reconstruction
EP3494885A1 (en) * 2017-12-07 2019-06-12 Koninklijke Philips N.V. Apparatus for presentation of dark field x-ray image information

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9332955B2 (en) 2011-10-04 2016-05-10 Koninklijke Philips N.V. Adaptive dual-pass targeted reconstruction and acquisition
JP6209520B2 (en) 2011-10-24 2017-10-04 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. Perfusion imaging
JP6495923B2 (en) * 2013-09-17 2019-04-03 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. Descattering method via energy calibration
US20150198684A1 (en) * 2014-01-13 2015-07-16 Tamer BASHA System and Method For Flexible Automated Magnetic Resonance Imaging Reconstruction
US20150269753A1 (en) * 2014-03-18 2015-09-24 Siemens Medical Solutions Usa, Inc. Multi-Output With Acquisition Replay System and Method
CN105761293B (en) * 2016-03-17 2019-01-08 上海联影医疗科技有限公司 Medical imaging procedure and system
EP3817663A1 (en) * 2018-08-07 2021-05-12 Siemens Medical Solutions USA, Inc. Adaptive compton camera for medical imaging
US11037338B2 (en) * 2018-08-22 2021-06-15 Nvidia Corporation Reconstructing image data

Family Cites Families (35)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IL66127A (en) * 1982-06-24 1987-11-30 Israel State Method and apparatus for measuring the index of refraction of fluids
US5841140A (en) * 1997-01-08 1998-11-24 Smv America, Inc. Gamma camera for pet and spect studies
US7209779B2 (en) * 2001-07-17 2007-04-24 Accuimage Diagnostics Corp. Methods and software for retrospectively gating a set of images
US6664543B2 (en) * 2002-01-28 2003-12-16 Cti Pet Systems, Inc. Continuous sampling and digital integration for PET scintillation
US7132664B1 (en) * 2002-11-09 2006-11-07 Crosetto Dario B Method and apparatus for improving PET detectors
US7983735B2 (en) * 2003-04-15 2011-07-19 General Electric Company Simulation of nuclear medical imaging
RU2233119C1 (en) * 2003-06-24 2004-07-27 Центральный научно-исследовательский рентгено-радиологический институт Method for predicting cardiac sarcoidosis
CN1846231A (en) * 2003-07-21 2006-10-11 派昂公司 Method and system for identifying optimal image within a series of images that depict a moving organ
EP1665125A4 (en) * 2003-09-02 2007-10-03 Ludwig Inst Cancer Res Data driven motion correction for nuclear imaging
US7045802B2 (en) * 2003-11-26 2006-05-16 General Electric Company Method and apparatus for coincidence imaging digital triggering
US9040016B2 (en) * 2004-01-13 2015-05-26 Biosensors International Group, Ltd. Diagnostic kit and methods for radioimaging myocardial perfusion
US7968851B2 (en) * 2004-01-13 2011-06-28 Spectrum Dynamics Llc Dynamic spect camera
WO2007010537A2 (en) * 2005-07-19 2007-01-25 Spectrum Dynamics Llc Reconstruction stabilizer and active vision
WO2007074467A2 (en) * 2005-12-28 2007-07-05 Spectrum Dynamics Llc Gating with anatomically varying durations
US8586932B2 (en) * 2004-11-09 2013-11-19 Spectrum Dynamics Llc System and method for radioactive emission measurement
WO2006051531A2 (en) * 2004-11-09 2006-05-18 Spectrum Dynamics Llc Radioimaging
EP1778957A4 (en) * 2004-06-01 2015-12-23 Biosensors Int Group Ltd Radioactive-emission-measurement optimization to specific body structures
CN100583155C (en) * 2004-06-18 2010-01-20 皇家飞利浦电子股份有限公司 Artifact reduction
JP4353040B2 (en) * 2004-09-22 2009-10-28 株式会社島津製作所 Positron CT system
US8423125B2 (en) * 2004-11-09 2013-04-16 Spectrum Dynamics Llc Radioimaging
US20060239398A1 (en) * 2005-03-07 2006-10-26 Fused Multimodality Imaging, Ltd. Breast diagnostic apparatus for fused SPECT, PET, x-ray CT, and optical surface imaging of breast cancer
JP2008539414A (en) * 2005-04-27 2008-11-13 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Temporal sampling of ECG gates in the modeling of cardiac dynamics
BRPI0611707A2 (en) 2005-06-29 2012-04-24 Compumedics Ltd conductive bridge sensor assembly
EP1902328B1 (en) * 2005-07-08 2015-11-11 Wisconsin Alumni Research Foundation Constrained backprojection reconstruction method for undersampled mri
US8644910B2 (en) * 2005-07-19 2014-02-04 Biosensors International Group, Ltd. Imaging protocols
EP1952180B1 (en) * 2005-11-09 2017-01-04 Biosensors International Group, Ltd. Dynamic spect camera
RU2312602C2 (en) * 2005-11-23 2007-12-20 Государственное образовательное учреждение высшего профессионального образования "Московский государственный институт электронной техники" (технический университет) Method for producing tomographic images
JP5254810B2 (en) * 2006-02-28 2013-08-07 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Local motion compensation based on list mode data
GB0612234D0 (en) 2006-06-20 2006-08-02 Mallinckrodt Inc Method of making and using rubidium-81-containing compositions
US7569828B2 (en) * 2006-09-21 2009-08-04 Siemens Medical Solutions Usa, Inc. Application-driven optimization of acquisition and reconstruction of SPECT/PET projection data
US8047715B2 (en) * 2006-11-03 2011-11-01 Koninklijke Philips Electronics N.V. Multiple rotation C-arm
DE102007060689B4 (en) * 2007-12-17 2017-07-13 Siemens Healthcare Gmbh Method for recording angiographic data sets and magnetic resonance system therefor
EP2286274B1 (en) * 2008-05-28 2019-07-24 Koninklijke Philips N.V. Geometrical transformations preserving list-mode format
US7732774B2 (en) * 2008-09-19 2010-06-08 Jefferson Science Associates, Llc High resolution PET breast imager with improved detection efficiency
US7858944B2 (en) * 2009-01-23 2010-12-28 Jefferson Science Associates, Llc Dedicated mobile high resolution prostate PET imager with an insertable transrectal probe

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
None

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2464212A (en) * 2008-10-09 2010-04-14 Siemens Medical Solutions Analyzing Medical Image Data
GB2464212B (en) * 2008-10-09 2013-03-06 Siemens Medical Solutions Methods and apparatus for analyzing medical imaging data
US8693741B2 (en) 2008-10-09 2014-04-08 Siemens Medical Solutions Usa, Inc. Methods and apparatus for analyzing medical imaging data
CN102028494A (en) * 2011-01-24 2011-04-27 海纳医信(北京)软件科技有限责任公司 Method and system for processing cerebral perfusion image sequence
CN104598356A (en) * 2014-12-23 2015-05-06 沈阳东软医疗系统有限公司 Event ordering method and device
US10133952B2 (en) 2014-12-24 2018-11-20 Shenyang Neusoft Medical Systems Co., Ltd. Computer tomography image reconstruction
EP3494885A1 (en) * 2017-12-07 2019-06-12 Koninklijke Philips N.V. Apparatus for presentation of dark field x-ray image information
WO2019110313A1 (en) 2017-12-07 2019-06-13 Koninklijke Philips N.V. Apparatus for presentation of dark field x-ray image information
US11350897B2 (en) 2017-12-07 2022-06-07 Koninklijke Philips N.V. Apparatus for presentation of dark field X-ray image information

Also Published As

Publication number Publication date
CN102098964B (en) 2015-04-08
US20110105887A1 (en) 2011-05-05
WO2009150565A3 (en) 2010-06-24
RU2517586C2 (en) 2014-05-27
JP2014149308A (en) 2014-08-21
JP5864658B2 (en) 2016-02-17
RU2011100843A (en) 2012-07-20
US8437836B2 (en) 2013-05-07
JP2011523076A (en) 2011-08-04
EP2288291A2 (en) 2011-03-02
EP2288291B1 (en) 2018-11-14
CN102098964A (en) 2011-06-15

Similar Documents

Publication Publication Date Title
US8437836B2 (en) Reverse data reconstruction for optimal time sampling of counts in physiological list-mode nuclear imaging
DePuey Advances in SPECT camera software and hardware: currently available and new on the horizon
EP1952180B1 (en) Dynamic spect camera
US7507968B2 (en) Systems and methods for correcting a positron emission tomography emission image
US7026623B2 (en) Efficient single photon emission imaging
EP1971257B1 (en) Gating with anatomically varying durations
US8507869B2 (en) Methods and systems for adaptive tomographic imaging
EP1940290B1 (en) Patient scan time optimization for pet/spect imaging
US20110228999A1 (en) System and method of image reconstruction using a temporal subset of sparsified image data
US7907698B2 (en) Gated CT with irregular sampling for slow CT acquisition
JP2004530915A (en) Dual isotope testing in nuclear medicine imaging
JP2011200656A (en) System and method for tomographic data acquisition and image reconstruction
US8058625B2 (en) Limiting viewing angles in nuclear imaging
JP2021512335A (en) Correction method for improving quantification accuracy in list mode reconstruction
DePuey III Advances in cardiac processing software
JP4895080B2 (en) Myocardial blood flow quantification method, blood flow quantification program, and blood flow quantification system
GB2497834A (en) Analyzing pet medical imaging data
Nichols et al. Prospects for advancing nuclear cardiology by means of new detector designs
Gullberg et al. Cardiac single-photon emission-computed tomography using combined cone-beam/fan-beam collimation
Lagopati Nuclear Medicine Imaging Essentials
Germano et al. Physics and technical aspects of gated myocardial perfusion SPECT
Khatami The Impact of Simulated Image Acquisition Time Reduction on Image Quality and Parameters of Myocardial Perfusion Imaging Using a Dedicated Cardiac Camera
DiBella et al. Compartmental modeling of/sup 18/FDG in the heart using dynamic SPECT
Fatemi Respiratory Motion Detection and Compensation for PET Imaging of Patients with Lung Cancer
Liu et al. Quantitative cardiac SPECT/CT

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 200980121926.3

Country of ref document: CN

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 09762092

Country of ref document: EP

Kind code of ref document: A2

WWE Wipo information: entry into national phase

Ref document number: 2009762092

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 12995685

Country of ref document: US

WWE Wipo information: entry into national phase

Ref document number: 2011513080

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 172/CHENP/2011

Country of ref document: IN

WWE Wipo information: entry into national phase

Ref document number: 2011100843

Country of ref document: RU