Embodiments of the invention relate generally to imaging systems, and more specifically, to a system and method for performing fault-tolerant reconstruction of an image of an object.
Positron Emission Tomography (“PET”) imaging is a non-invasive nuclear imaging technique that involves the creation of tomographic images from positron emitting radionuclides injected into a subject of interest. In many PET imaging systems, a radionuclide-labeled agent is administered to the subject while positioned within a detector ring. As the radionuclides decay, positively charged anti-electrons known as “positrons” are emitted therefrom. As these positrons travel through the tissues of the subject, they lose kinetic energy and ultimately combine with an electron, resulting in mutual annihilation. The positron annihilation results in a pair of oppositely directed annihilation photons, e.g., gamma rays, being emitted at approximately 511 keV, which are subsequently detected by PET detectors/sensors, i.e., scintillators such as crystals, within the detector ring. When struck by an annihilation photon, each PET detector converts the energy of the photon into a burst of light photons of lower energy; this burst of light is then detected by a photovoltaic component such as a photodiode.
The signals from the photovoltaics are processed as incidences of annihilation photons typically referred to as “PET events.” When two PET events occur at oppositely positioned PET detectors at approximately the same time, a “coincidence” is registered. Data sorting units process the coincidences to determine true coincidence events and to sort out data representing dead times and single gamma ray detections. The coincidence events are binned and integrated to form frames of PET data, which may be reconstructed into images depicting the distribution of the radionuclide-labeled agent and/or metabolites thereof in the subject.
Many PET detectors, however, are subject to various modes of failure that typically generate artifacts and/or quantitative inaccuracies in reconstructed images, which in turn may result in false and/or inconclusive diagnosis. For example, PET detector failures may cause one or more of the following errors: insufficient counts for a PET detector; excessive counts for a PET detector; and/or misplaced counts for a PET detector, i.e., where a gamma ray strikes one PET detector, but is recorded against a different PET detector.
While some PET imaging systems seek to correct the effects of PET detector failures, such PET imaging systems typically are able to only detect and correct for intermittent PET detector failures that occur when a PET detector “flickers”, i.e., toggles between providing good data and no data. Moreover, many such PET imaging systems require the PET data to be acquired and/or stored in what is known as “List Mode,” which often entails storing coincidences in chronological order within an array for the purpose of calculating a fractional weight of valid PET data with respect to a period of time. This approach, however, is limited to correcting the very specific type of PET detector failure, mentioned above, in which a PET detector is either working perfectly or providing no data.
As such, there exist types of PET detector failures that can corrupt PET data while typically going undetected by prior art PET imaging systems. For example, many PET imaging systems usually fail to detect and/or correct for corrupted detector crystal maps, large energy gain errors, false detector triggering causing significant detector dead time, and/or time stamp corruption.
What is needed, therefore, is an improved system and method for performing fault-tolerant reconstruction of an image of an object.
In an embodiment, a method for performing fault-tolerant reconstruction of an image of an object is provided. The method includes acquiring a dataset corresponding to the object via a plurality of sensors, and detecting anomalous data within the dataset based at least in part on a statistical difference between the anomalous data and reference data within the dataset. The anomalous data is acquired by at least a first sensor of the plurality, and the reference data is acquired by at least a second sensor of the plurality that neighbors the first sensor.
In another embodiment, a system for performing fault-tolerant reconstruction of an image of an object is provided. The system includes a plurality of sensors and a controller. The plurality of sensors is operative to acquire a dataset corresponding to the object. The controller is in electronic communication with the plurality of sensors and operative to detect anomalous data within the dataset based at least in part on a statistical difference between the anomalous data and reference data within the dataset. The anomalous data is acquired by at least a first sensor of the plurality, and the reference data is acquired by at least a second sensor of the plurality that neighbors the first sensor.
In yet another embodiment, a non-transitory computer readable medium storing instructions is provided. The stored instructions are configured to adapt a controller to detect anomalous data within a dataset based at least in part on a statistical difference between the anomalous data and reference data within the dataset. The dataset corresponds to an object and is acquired by a plurality of sensors. The anomalous data is acquired by at least a first sensor of the plurality, and the reference data is acquired by at least a second sensor of the plurality that neighbors the first sensor.
The present invention will be better understood from reading the following description of non-limiting embodiments, with reference to the attached drawings, wherein below:
Reference will be made below in detail to exemplary embodiments of the invention, examples of which are illustrated in the accompanying drawings. Wherever possible, the same reference characters used throughout the drawings refer to the same or like parts, without duplicative description.
As used herein, the terms “substantially,” “generally,” and “about” indicate conditions within reasonably achievable manufacturing and assembly tolerances, relative to ideal desired conditions suitable for achieving the functional purpose of a component or assembly. As used herein, “electrically coupled”, “electrically connected”, and “electrical communication” mean that the referenced elements are directly or indirectly connected such that an electrical current may flow from one to the other. The connection may include a direct conductive connection, i.e., without an intervening capacitive, inductive or active element, an inductive connection, a capacitive connection, and/or any other suitable electrical connection. Intervening components, such as a fiber optical connection, may be present. The term “real-time,” as used herein, means a level of processing responsiveness that a user senses as sufficiently immediate or that enables the processor to keep up with an external process.
Further, while the embodiments disclosed herein are described with respect to PET imaging systems, it is to be understood that embodiments of the present invention may be applicable to other Nuclear Medicine (“NM”) imaging systems, e.g., Single Photon Emission Computed Tomography (“SPECT”); MRI systems; Computed Tomography (“CT”) imaging systems; ultrasound imaging systems; and imaging systems which utilize more than one detector/sensor or view to reconstruct an image of an object, i.e., a photosensitive chip within a digital camera. Further still, as will be appreciated, embodiments of the present invention related imaging systems may be used to analyze tissue generally and are not limited to human tissue.
Referring now to
As will be appreciated, in embodiments, the sensors 24 may be disposed within the detector assembly 14 as a ring (best seen in
As shown in
The detector assembly 14 may further include a PET scan sub-controller 60 and one or more acquisition circuits 62.
The HMI interface 26 may include an input device 64, e.g., a keyboard, mouse, etc., a display screen 66, a processor 68, and/or a memory device 70, and may electronically communicate with the PET scanner controller 60, data acquisition sub-controller 32, and/or the image reconstruction processor 34 via communication link 72.
During operation of the system 10, in accordance with an embodiment, the patient 12 is injected with a radionuclide-labeled agent and positioned within the central opening 30. As positrons are emitted from the radionuclide-labeled agent, they eventually annihilate with electrons to produce/generate gamma rays/photons having energies of about 511 keV. Specifically, each positron/electron annihilation generates two gamma rays having generally opposite velocities, and under ideal conditions, i.e., in the absence of scattering and/or other attenuation effects, are detected nearly simultaneously by two oppositely disposed sensors 24, i.e., a single annihilation event may result in two (2) PET events occurring at approximately the same time to form a coincidence. The location of the annihilation event that generated a coincidence can be derived from lines of response (“LORs”) corresponding to the two gamma rays that generated the coincidence.
Accordingly, as PET events occur/are detected by the sensors 24, the acquisition circuits 62 receive signals, e.g., analog voltages, generated by the sensors 24 in response to detecting the gamma rays. As will be appreciated, the signals convey three-dimensional (“3D”) positional information and total energy for each gamma ray of each PET event, and the acquisition circuits 62 may produce an event detection pulse that indicates the time/moment the PET events were detected. The signals are then transmitted through the communication links 22 to the data acquisition sub-controller 32, which performs various data enhancing techniques on the signals prior to transmitting the signals to the image reconstruction processor 34.
For example, in embodiments, the signals received by the data acquisition sub-controller 32 from the acquisition circuits 62 may be sent to the event locator module 40, which in turn processes the 3D positional and total energy information of the gamma rays to identify valid PET events. As will be appreciated, in embodiments, the event locator module 40 may determine the exact sensor 24 that detected a particular PET event, and then transmit information, e.g., location, energy, detecting sensor 24 identifier, and/or other information regarding the valid PET event, within an event data packet that is sent to the coincidence detector module 42 and/or image reconstruction module 44.
The coincidence detector module 42 receives event data packets from the event locator module 40 and determines if any two valid PET events form a coincidence, i.e., there is a probability, which in embodiments may be a high probability, that the two valid PET events were generated by gamma rays from the same positron/electron annihilation. Coincident PET event pairs are then recorded as a dataset 74 that corresponds to the object 12, which is then communicated to the image reconstruction module 44. In embodiments, the dataset 74 may be stored in one or more of the memory devices 38, 50, and/or 70.
After having been processed by the image reconstruction module 44, the dataset 74 is electronically communicated to the sorter/histogrammer module 48, which generates a histogram and/or other data structure capable of readily conveying statistical data of the dataset 74. In embodiments, the histogram may be stored as a data array within memory device 50, and accessed by the image processor 54, which reconstructs one or more images of the object 12 based on the dataset 74. In embodiments, the dataset and/or histogram 74 may be/include an array having a plurality of indices each corresponding to a different combination/pair of sensors 24, wherein the value of the array at a particular index represents data corresponding to the two sensors 24, e.g., PET events, coincidence count, etc. As will be appreciated, in embodiments, the coincidence count at an index of the array may be the total number of coincidences for the corresponding pair of sensor 24 during an image frame, which as used herein, refers to the length of time that PET data was collected to generate a single image of the one or more images. As will be understood, such images may be individual still images and/or form part of a series of moving images. Accordingly, in embodiments, the dataset 74 may be at least one of a coincidence count detector map, a singles map, a coincidence-qualified singles count detector map, and a block busy map. The one or more images may then be electronically communicated to the display screen 66 on the HMI 26 for viewing by a technologist and/or physician.
Turning now to
Referring now to
As will be understood, embodiments of the system 10 take advantage of the spatially slowly changing nature of the dataset 74, i.e., the tendency of neighboring sensors 24 to have similar PET event counts, by looking for outlying/anomalous data 82 within the dataset 74 to determine the likelihood that a particular sensor 24 has failed.
As will be explained in greater detail below, the dataset 74 may be transformed into one or more two-dimensional (“2D”) maps which may be organized so as to retain the neighborhoods 100 of the sensors 24. For example, in embodiments, the count statistics of sensors 24 disposed around the ring 92 at a given axial location may be plotted on a single row, while sensors 24 located along the axial direction may be plotted on a single column. As will be appreciated, the values across the locations within such 2D maps tend to be smoothly varying, with sudden changes in value indicating likely abnormalities in the corresponding sensors 24.
For example, illustrated in
Next, outliers may be searched for by selecting elements within the maps that are more than a certain number of standard deviations from the mean. The outliers may then be removed so that the standard deviation of the remaining elements may be recomputed. As will be appreciated, the searching/removing of the outliers, and the computation of the standard deviation of the remaining elements may be repeated in an iterative manner, referred to herein as “significant z-processing,” until new outliers are no longer found.
After the completion of significant z-processing, elements that fall outside of a predetermined threshold may be found for both the neighbor difference and the z map 114. As will be appreciated, this may generate two masks: M1 and M3, where elements are=1 if the corresponding unit 98 is an outlier. Each neighbor difference may then be checked to see if it is greater than some practical threshold, e.g., mask M2, =1 if difference is significant. Thus, the final expression for outliers may be given by (M1&&M2)∥M3. Accordingly, as shown above, anomalous data 82 may initially be considered merely “suspect data” if its statistical difference 84 with respect to the reference data 86 exceeds an initial threshold, and then subsequently classified as “faulty data” if the statistical difference 84 exceed a subsequent threshold beyond the initial threshold. Further, while the above method of detecting outliers was described on a per/unit 98 basis, is should be understood that other groupings of sensors, e.g., modules, blocks, etc., may be used.
Returning back to
For example, in embodiments, reconstruction 130 of the image may be accomplished via Ordered Subset Expectation Maximization (“OSEM”) with an image update equation of:
where {circumflex over (f)}j(n+1) the next estimate of voxel j based on the current voxel estimate {circumflex over (f)}j(n), S represents the projection subsets, p represents the measured data, Hij are the elements of the system model H and represent the probability that an emission from voxel j is detected in projection, or LOR, i. As will be appreciated, the system model accounts for factors of detector geometry, normalization, point spread function, and/or attenuation. As such, expectation maximization may be performed by forwarding projection {circumflex over (f)}j(n) into the projection (sinogram) domain by multiplying it with elements Hij. The ratio of the measured data pi may then be taken with the projection of the current voxel estimate. Next, the resulting ratio for the projections i in subset S may be backprojected into the image domain to obtain an update factor for the current voxel estimate. The update factor may then be multiplied with the current estimate and divided by a weighting term based on the system model. The above process may be repeated until convergence and/or for a preset/specified number of iterations.
As will be understood, embodiments of the above resilient reconstruction may rely on the weighting of measured data based on the description of the system model. Thus, the weighting term in such embodiments indicates confidence in the data being backprojected. Therefore, a low value in the system matrix Hij gives a smaller weight of the corresponding projection i in the measured data pi on image voxel j, while a larger value will give pi a larger impact on voxel j. Thus, as will also be understood, by setting Hij equal to zero for the projections pi for detected anomalous data 82 acquired by failed sensors, e.g., sensor 88, the corresponding LORs are removed from the image update equation, thus minimizing the impact of the failed sensors on the reconstructed image. In other words, as shown in
Accordingly, in embodiments, the dataset 74 may be a singles histogram that is inspected to identify outlying PET events counts. As will be appreciated, such inspection may be performed on the dataset 74 in each PET imaging bad position, which in turn ensures that any intermittent sensors 24 are accounted for, even if they are not flagged during a daily quality assurance (“DQA”) process/systems check. As will be understood, DQA system checks evaluate sensor 24 integrity so as to distinguish between general problems in sensor 24 calibration and inconsistent data from a few sensors 24 within a well-calibrated system. As such, in embodiments, the method of detecting 80 anomalous data 82 described above may be performed on DQA acquired data. Thus, if the system 10 fails conventional DQA criteria, but passes DQA when one or more sensors 24, e.g., five (5), are excluded, further scanning via the system 10 may be permitted.
Further, in embodiments, the method 76 may also include generating 136 a notification, e.g., e-mail, sms message, computer popup, audio/visual alarm, etc., that includes an indication/indicator that the first sensor 88 has acquired anomalous data 82. In other words, when the system 10 determines that the statistical difference between the anomalous data 82 and the reference data 84 indicates a high likelihood that one or more sensors 88 have failed, a message may be generated by the controller 16 and sent to appropriate repair personnel, and/or technologist operating the system 10.
Finally, it is also to be understood that the system 10 may include the necessary electronics, software, memory, storage, databases, firmware, logic/state machines, microprocessors, communication links, displays or other visual or audio user interfaces, printing devices, and any other input/output interfaces to perform the functions described herein and/or to achieve the results described herein. For example, as previously mentioned, the system may include at least one processor and system memory/data storage structures, which may include random access memory (RAM) and read-only memory (ROM). The at least one processor of the system 10 may include one or more conventional microprocessors and one or more supplementary co-processors such as math co-processors or the like. The data storage structures discussed herein may include an appropriate combination of magnetic, optical and/or semiconductor memory, and may include, for example, RAM, ROM, flash drive, an optical disc such as a compact disc and/or a hard disk or drive.
Additionally, a software application that adapts the controller to perform the methods disclosed herein may be read into a main memory of the at least one processor from a computer-readable medium. Thus, embodiments of the present invention may perform the methods disclosed herein in real-time. The term “computer-readable medium,” as used herein, refers to any medium that provides or participates in providing instructions to the at least one processor of the system 10 (or any other processor of a device described herein) for execution. Such a medium may take many forms, including but not limited to, non-volatile media and volatile media. Non-volatile media include, for example, optical, magnetic, or opto-magnetic disks, such as memory. Volatile media include dynamic random access memory (DRAM), which typically constitutes the main memory. Common forms of computer-readable media include, for example, a floppy disk, a flexible disk, hard disk, solid state drive (SSD), magnetic tape, any other magnetic medium, a CD-ROM, DVD, any other optical medium, a RAM, a PROM, an EPROM or EEPROM (electronically erasable programmable read-only memory), a FLASH-EEPROM, any other memory chip or cartridge, or any other medium from which a computer can read.
While in embodiments, the execution of sequences of instructions in the software application causes at least one processor to perform the methods/processes described herein, hard-wired circuitry may be used in place of, or in combination with, software instructions for implementation of the methods/processes of the present invention. Therefore, embodiments of the present invention are not limited to any specific combination of hardware and/or software.
It is further to be understood that the above description is intended to be illustrative, and not restrictive. For example, the above-described embodiments (and/or aspects thereof) may be used in combination with each other. Additionally, many modifications may be made to adapt a particular situation or material to the teachings of the invention without departing from its scope.
For example, in an embodiment, a method for performing fault-tolerant reconstruction of an image of an object is provided. The method includes acquiring a dataset corresponding to the object via a plurality of sensors, and detecting anomalous data within the dataset based at least in part on a statistical difference between the anomalous data and reference data within the dataset. The anomalous data is acquired by at least a first sensor of the plurality, and the reference data is acquired by at least a second sensor of the plurality that neighbors the first sensor. In certain embodiments, the method further includes reconstructing the image based at least in part on the dataset. In such embodiments, the anomalous data contributes less to the reconstructed image than the reference data. In certain embodiments, the anomalous data does not contribute to the reconstructed image. In certain embodiments, reconstructing the image based at least in part on the dataset includes applying a scattering model to the dataset. In such embodiments, a scaling process of the scattering model that determines an amplitude of a scatter contribution of the scattering model ignores lines of response corresponding to the first sensor. In certain embodiments, the statistical difference exceeds an initial threshold. In certain embodiments, the statistical difference exceeds a subsequent threshold beyond the initial threshold. In certain embodiments, the dataset corresponds to positron emissions from the object. In certain embodiments, the dataset includes at least one of a coincidence count detector map, a singles map, a coincidence-qualified singles count detector map, and a block busy map. In certain embodiments, the coincidence count detector map has a plurality of indices each corresponding to a coincidence count for a different combination of sensors of the plurality. In certain embodiments, the coincidence count is a total coincidence count for an imaging scan during which the dataset is acquired. In certain embodiments, the method further includes generating a notification that includes an indication that the first sensor has acquired the anomalous data.
Yet other embodiments provide for a system for performing fault-tolerant reconstruction of an image of an object. The system includes a plurality of sensors and a controller. The plurality of sensors is operative to acquire a dataset corresponding to the object. The controller is in electronic communication with the plurality of sensors and operative to detect anomalous data within the dataset based at least in part on a statistical difference between the anomalous data and reference data within the dataset. The anomalous data is acquired by at least a first sensor of the plurality, and the reference data is acquired by at least a second sensor of the plurality that neighbors the first sensor. In certain embodiments, the controller is further operative to reconstruct the image based at least in part on the dataset. In such embodiments, the anomalous data contributes less to the reconstructed image than the reference data. In certain embodiments, the anomalous data does not contribute to the reconstructed image. In certain embodiments, the controller is further operative to reconstruct the image by applying a scattering model to the dataset. In such embodiments, a scaling process of the scattering model that determines an amplitude of a scatter contribution of the scattering model ignores lines of response corresponding to the first sensor. In certain embodiments, the statistical difference exceeds an initial threshold. In certain embodiments, the statistical difference exceeds a subsequent threshold beyond the initial threshold. In certain embodiments, the dataset includes at least one of a coincidence count detector map, a singles map, a coincidence-qualified singles count detector map, and a block busy map. In certain embodiments, the coincidence count detector map has a plurality of indices each corresponding to a coincidence count for a different combination of sensors of the plurality. In certain embodiments, the coincidence count is a total coincidence count for an imaging scan during which the dataset is acquired. In certain embodiments, the controller is further operative to generate a notification that includes an indication that the first sensor has acquired the anomalous data.
Yet still other embodiments provide for a non-transitory computer readable medium storing instructions. The stored instructions are configured to adapt a controller to detect anomalous data within a dataset based at least in part on a statistical difference between the anomalous data and reference data within the dataset. The dataset corresponds to an object and is acquired by a plurality of sensors. The anomalous data is acquired by at least a first sensor of the plurality, and the reference data is acquired by at least a second sensor of the plurality that neighbors the first sensor. In certain embodiments, the stored instructions are further configured to adapt the controller to reconstruct the image based at least in part on the dataset. In such embodiments, the anomalous data contributes less to the reconstructed image than the reference data. In certain embodiments, the anomalous data does not contribute to the reconstructed image. In certain embodiments, the controller reconstructs the image by applying a scattering model to the dataset. In such embodiments, a scaling process of the scattering model that determines an amplitude of a scatter contribution of the scattering model ignores lines of response corresponding to the first sensor. In certain embodiments, the statistical difference exceeds an initial threshold. In certain embodiments, the statistical difference exceeds a subsequent threshold beyond the initial threshold. In certain embodiments, the dataset includes at least one of a coincidence count detector map, a singles map, a coincidence-qualified singles count detector map, and a block busy map. In certain embodiments, the coincidence count detector map has a plurality of indices each corresponding to a coincidence count for a different combination of sensors of the plurality. In certain embodiments, the coincidence count is a total coincidence count for an imaging scan during which the dataset is acquired. In certain embodiments, the stored instructions are further configured to adapt the controller to generate a notification that includes an indication that the first sensor has acquired the anomalous data.
Accordingly, by detecting failed sensors by comparing PET data acquired by a sensor to PET data acquired by neighboring sensors, some embodiments of the present invention may provide for the detection and correction of a wider variety of PET sensor failures than currently provided by prior art PET systems. In other words, and as described above, some embodiments of the present invention take advantage of the spatially slowly changing nature of acquired PET data, i.e., the tendency of neighboring sensors to have similar PET event counts, by looking for outlying/anomalous data within an acquired PET dataset to determine the likelihood that a particular sensor has failed.
Thus, unlike traditional PET systems, some embodiments of the present invention provide for the ability to detect failed sensors without requiring the PET data to be acquired and/or formatted in accordance with List Mode. As such, by detecting anomalous data within each acquired dataset, some embodiments may provide for the detection of intermittent sensor failures that typically go undetected by traditional DQA system checks. Therefore, some embodiments of the present invention provide for the ability to detect and to correct for a wider variety of sensor failures than many traditional PET systems.
Further, by generating a notification message for maintenance personnel that identifies one or more failed PET sensors, some embodiments of the present invention provide for improved diagnostic and repair times, and/or reduced maintenance costs. For example, some embodiments may prevent scenarios in which a patient is injected with a radionuclide-labeled agent only to be prevented from completing a PET scan due to a failed sensor.
Additionally, while the dimensions and types of materials described herein are intended to define the parameters of the invention, they are by no means limiting and are exemplary embodiments. Many other embodiments will be apparent to those of skill in the art upon reviewing the above description. The scope of the invention should, therefore, be determined with reference to the appended claims, along with the full scope of equivalents to which such claims are entitled. In the appended claims, the terms “including” and “in which” are used as the plain-English equivalents of the respective terms “comprising” and “wherein.” Moreover, in the following claims, terms such as “first,” “second,” “third,” “upper,” “lower,” “bottom,” “top,” etc. are used merely as labels, and are not intended to impose numerical or positional requirements on their objects. Further, the limitations of the following claims are not written in means-plus-function format are not intended to be interpreted as such, unless and until such claim limitations expressly use the phrase “means for” followed by a statement of function void of further structure.
This written description uses examples to disclose several embodiments of the invention, including the best mode, and also to enable one of ordinary skill in the art to practice the embodiments of invention, including making and using any devices or systems and performing any incorporated methods. The patentable scope of the invention is defined by the claims, and may include other examples that occur to one of ordinary skill in the art. Such other examples are intended to be within the scope of the claims if they have structural elements that do not differ from the literal language of the claims, or if they include equivalent structural elements with insubstantial differences from the literal languages of the claims.
As used herein, an element or step recited in the singular and proceeded with the word “a” or “an” should be understood as not excluding plural of said elements or steps, unless such exclusion is explicitly stated. Furthermore, references to “one embodiment” of the present invention are not intended to be interpreted as excluding the existence of additional embodiments that also incorporate the recited features. Moreover, unless explicitly stated to the contrary, embodiments “comprising,” “including,” or “having” an element or a plurality of elements having a particular property may include additional such elements not having that property.
Since certain changes may be made in the above-described invention without departing from the spirit and scope of the invention herein involved, it is intended that all of the subject matter of the above description shown in the accompanying drawings shall be interpreted merely as examples illustrating the inventive concept herein and shall not be construed as limiting the invention.
This application claims priority to and the benefit of, U.S. Provisional Pat. App. Ser. No. 62/488,345, filed Apr. 21, 2017, the disclosure of which is incorporated herein by reference in its entirety.
Number | Date | Country | |
---|---|---|---|
62488345 | Apr 2017 | US |