PET device, method and program

Abstract

To quantitatively and objectively evaluate influence on image quality due to a defective channel.SOLUTION: A PET device comprises an acquisition unit and a generation unit. The acquisition unit obtains information on a defective channel of a PET detector at a second timing after a first timing corresponding to a first sensitivity map, which is a sensitivity map of the PET detector at the first timing stored in a storage unit. The generation unit generates a second sensitivity map, which is a sensitivity map of the PET detector at the second timing, based on the information on the defective channel.SELECTED DRAWING: Figure 1

Description

[Technical Field]

[0001] The embodiments disclosed in the present specification and drawings relate to a PET device, a method, and a program. [Background technology]

[0002] A PET detector in a PET (Positron Emission Tomography) device includes a scintillator array and a photodetector array. The scintillator array is composed of a plurality of scintillators. The photodetector array is composed of a plurality of photodetectors (photodetecting elements). Here, for example, the photodetector array may be composed of a plurality of PMTs (Photomultiplier Tubes) as the plurality of photodetectors. PET devices equipped with such photodetector arrays employ, for example, a multiplexing system. In this case, the total number of channels (total number of signal output channels) in the PET detector is about several hundred channels because the number of PMTs is smaller than the number of scintillators that make up the scintillator array.

[0003] Because the fluorescence from many scintillators is incident on a small number of PMTs and converted into signals, a loss of signal output from one PMT has a non-negligible effect on a wide area of ​​the image. Therefore, if a phantom is actually imaged in a health check to diagnose (inspect) a PET detector malfunction, the presence or absence of a malfunction can be qualitatively determined by the user visually inspecting the image. Therefore, if a PET detector malfunctions, there is no option other than to repair or replace the PET detector, and in this sense, no special judgment is required.

[0004] On the other hand, photodetector arrays may be composed of multiple silicon photomultipliers (SiPMs) as multiple photodetectors. In this case, one scintillator is optically coupled to one SiPM (i.e., a scintillator and a SiPM are optically coupled one-to-one), or a small number of scintillators are optically bonded to SiPMs. As a result, the number of independent channels (signal output channels) in the entire PET detector can reach tens of thousands to hundreds of thousands. Unlike PET detectors that use multiple photomultiplier tubes (PMTs) as a photodetector array, if a defective channel is included, the image quality change and its impact are not so noticeable that a user can immediately notice it visually. Instead, there are variations depending on the number of defective channels and their relative positions within the PET detector. Here, a defective channel refers, for example, to a channel corresponding to a malfunctioning photodetector. Image uniformity changes do not necessarily occur smoothly over a wide range and can occur on a pixel-by-pixel or component-by-component basis. For this reason, visual confirmation of statistically noisy images obtained by imaging a phantom lacks quantitativeness and objectivity. In other words, when one scintillator is optically coupled to one SiPM, or when a small number of scintillators are optically bonded to SiPMs, the impact of individual channel defects on image quality (or the image) is relatively minor, and the criteria for determining how to respond are not clear.

[0005] Therefore, for PET detectors with many independent channels, it is necessary to quantitatively and objectively evaluate the impact of defective channels on image quality. It is important for users to make correct decisions based on this evaluation from the perspective of ensuring safety and performance in clinical settings, and from the perspective of appropriate cost allocation in manufacturing and service settings. [Prior art documents] [Patent documents]

[0006] [Patent Document 1] US Patent Application Publication No. 2018 / 0059267 [Patent Document 2] US Patent Application Publication No. 2009 / 0224164 [Patent Document 3] U.S. Patent No. 9,606,245 Summary of the Invention [Problem to be solved by the invention]

[0007] One of the problems to be solved by the embodiments disclosed in this specification and the drawings is to quantitatively and objectively evaluate the impact of a bad channel on image quality. However, the problem to be solved by the embodiments disclosed in this specification and the drawings is not limited to the above problem. Problems corresponding to the effects of each configuration shown in the embodiments described below can also be positioned as other problems. [Means for solving the problem]

[0008] According to an embodiment, the PET device includes an acquisition unit and a generation unit. The acquisition unit acquires information about a defective channel of the PET detector at a second timing that is later than a first timing corresponding to a first sensitivity map that is a sensitivity map of the PET detector at the first timing stored in a storage unit. The generation unit generates a second sensitivity map that is a sensitivity map of the PET detector at the second timing based on the information about the defective channel. [Brief explanation of the drawings]

[0009] [Figure 1] FIG. 1 is a diagram showing an example of the configuration of a PET apparatus according to the first embodiment. [Figure 2] FIG. 2 is a diagram illustrating an example of a sensitivity map according to the first embodiment. [Figure 3] FIG. 3 is a diagram for explaining an example of a method for generating a sensitivity map according to the first embodiment. [Figure 4] FIG. 4 is a flowchart showing an example of the flow of processing executed by the PET apparatus according to the first embodiment. [Figure 5]FIG. 5 is a diagram showing an example of the sensitivity map generated in step S103 according to the first embodiment. [Figure 6] FIG. 6 is a diagram showing an example of the sensitivity map generated in step S103 according to the first embodiment. [Figure 7] FIG. 7 is a flowchart showing an example of the flow of the evaluation process executed by the PET apparatus according to the second embodiment. [Figure 8] FIG. 8 is a flowchart showing an example of the flow of the evaluation process executed by the PET apparatus according to the third embodiment. [Figure 9] FIG. 9 is a flowchart showing an example of the flow of the specific processing executed in step S301 by the PET apparatus according to the third embodiment. [Figure 10] FIG. 10 is a flowchart showing an example of the flow of the process of determining the number of allowable defective channels executed by the PET apparatus according to the fourth embodiment. DETAILED DESCRIPTION OF THE INVENTION

[0010] Hereinafter, embodiments of a PET device, a method, and a program will be described in detail with reference to the drawings. Note that the PET device, a method, and a program according to the present application are not limited to the embodiments shown below. Furthermore, the content described in one embodiment is, in principle, similarly applied to other embodiments and modifications. Furthermore, the content described in one modification is, in principle, similarly applied to the embodiment and other modifications.

[0011] (First embodiment) Fig. 1 is a diagram showing an example of the configuration of a PET device 100 according to the first embodiment. As shown in Fig. 1, the PET device 100 according to the first embodiment includes a gantry device 10 and a console device 20. The PET device 100 is, for example, an example of a nuclear medicine diagnosis device.

[0012] The gantry 10 includes a PET detector 101 , a counting information generating circuit 102 , a tabletop 103 , a bed 104 , and a bed driver 105 .

[0013] The PET detector 101 includes a plurality of detector modules. The plurality of detector modules detect radiation by detecting scintillation light (fluorescence) that is re-emitted when a substance that has been excited by pair annihilation gamma rays emitted from positrons in the subject P interacts with a light emitter (scintillator) and transitions back to the ground state. The plurality of detector modules detects radiation energy information of the pair annihilation gamma rays emitted from the positrons in the subject P. The plurality of detector modules are arranged to surround the subject P in a ring shape.

[0014] The detector module includes, for example, a scintillator array, a photodetector array, and a light guide.

[0015] The scintillator array includes a plurality of scintillators arranged two-dimensionally. The scintillators convert incident pair annihilation gamma rays emitted from positrons in the subject P into scintillation photons (optical photons) and output the scintillation light. The scintillators are formed of scintillator crystals suitable for TOF measurement and energy measurement, such as LaBr3 (Lanthanum Bromide), LYSO (Lutetium Yttrium Oxyorthosilicate), LSO (Lutetium Oxyorthosilicate), and LGSO (Lutetium Gadolinium Oxyorthosilicate).

[0016] The light guide is made of a highly optically transparent plastic material or the like, and transmits the scintillation light output from the scintillator to the photodetector. Specifically, the light guide transmits the scintillation light to a SiPM (Silicon Photomultiplier), which will be described later.

[0017] The photodetector array includes a plurality of photodetectors (photodetecting elements) arranged two-dimensionally. For example, SiPMs are used as the photodetectors. In this embodiment, one scintillator is optically coupled to one SiPM. That is, the scintillator and the SiPM are optically coupled one-to-one. In this case, a combination of one scintillator and one SiPM optically coupled to this scintillator corresponds to one channel. Note that a small number of scintillators may be optically bonded to the SiPM. In this case, a combination of a small number of scintillators and one SiPM optically coupled to these small number of scintillators corresponds to one channel. In any case, the total number of channels (number of signal output channels) of the PET detector 101 ranges from tens of thousands to hundreds of thousands. The SiPM detects scintillation light emitted by the scintillator, generates an electrical signal corresponding to the scintillation light, and outputs the generated electrical signal as an output signal.

[0018] The PET detector 101 includes a counting information generating circuit 102 for each detector module.

[0019] The counting information generating circuit 102 converts the output signal from the PET detector 101 into digital data, thereby generating the digital data as counting information. This counting information includes the detection position, energy value, and detection time of the pair annihilation gamma ray. For example, the counting information generating circuit 102 identifies multiple photodetectors that converted scintillation light into electrical signals at the same time. Then, the counting information generating circuit 102 identifies a scintillator number (P) indicating the position of the scintillator onto which the pair annihilation gamma ray is incident (the position of the scintillator). The counting information generating circuit 102 can use various methods to identify the position of the scintillator onto which the pair annihilation gamma ray is incident. For example, since one scintillator corresponds to one SiPM, the counting information generating circuit 102 may identify the position of the scintillator corresponding to the SiPM from which the output is obtained as the position of the scintillator onto which the pair annihilation gamma ray is incident.

[0020] Furthermore, the counting information generating circuit 102 identifies the energy value (E) of the pair annihilation gamma ray incident on the PET detector 101 by integrating the intensity of the electrical signal output from each photodetector. The counting information generating circuit 102 also identifies the detection time (T) at which the PET detector 101 detects scintillation light due to the pair annihilation gamma ray. The detection time (T) may be an absolute time or the elapsed time from the start of imaging. In this manner, the counting information generating circuit 102 generates counting information including the scintillator number (P), the energy value (E), and the detection time (T). The counting information generating circuit 102 then stores the generated counting information in the memory 130 of the console device 20.

[0021] The count information generating circuit 102 is realized by, for example, a processor.

[0022] The top board 103 is a bed on which the subject P is placed, and is placed on the bed 104. The bed driving unit 105 moves the top board 103 under the control of the bed control function 106d of the processing circuitry 106. For example, the bed driving unit 105 moves the subject P into the imaging opening of the gantry device 10 by moving the top board 103.

[0023] The console device 20 accepts user operations on the PET device 100, controls the acquisition of PET image data, and reconstructs (generates) the PET image data using counting information collected by the gantry device 10. As shown in Fig. 1, the console device 20 includes a processing circuit 106, an input interface 110, a display 120, and a memory 130. The processing circuit 106, the input interface 110, the display 120, and the memory 130 are connected via a bus.

[0024] The processing circuitry 106 includes a coincidence information generation function 106a, a reconstruction processing function 106b, a system control function 106c, a bed control function 106d, an acquisition function 106e, a generation function 106f, a calculation function 106g, a determination function 106h, and a display control function 106i. Each of the coincidence information generation function 106a, the reconstruction processing function 106b, the system control function 106c, the bed control function 106d, the acquisition function 106e, the generation function 106f, the calculation function 106g, the determination function 106h, and the display control function 106i is stored in the memory 130 in the form of a computer-executable program. The processing circuitry 106 is a processor that reads each program from the memory 130 and executes the read program to realize the corresponding function. In other words, the processing circuitry 106, after reading each program, has the functions shown in the processing circuitry 106 of FIG. 1. 1 has been described as realizing the coincidence counting information generating function 106a, the reconstruction processing function 106b, the system control function 106c, the bed control function 106d, the acquisition function 106e, the generating function 106f, the calculation function 106g, the determination function 106h, and the display control function 106i in a single processing circuit 106. However, the processing circuit 106 may be configured with a plurality of independent processors, and each processor may execute a respective program to realize the functions.

[0025] The reconstruction processing function 106b is, for example, an example of a reconstruction processing unit. The acquisition function 106e is, for example, an example of an acquisition unit. The generation function 106f is, for example, an example of a generation unit. The calculation function 106g is, for example, an example of a calculation unit. The determination function 106h is, for example, an example of a determination unit. The display control function 106i is, for example, an example of a display control unit.

[0026] The term "processor" used in the above description refers to circuits such as a CPU (Central Processing Unit), a GPU (Graphical Processing Unit), an Application Specific Integrated Circuit (ASIC), a programmable logic device (e.g., a Simple Programmable Logic Device (SPLD), a Complex Programmable Logic Device (CPLD), and a Field Programmable Gate Array (FPGA)). The processor implements functions by reading and executing programs stored in memory 130.

[0027] The coincidence counting information generating function 106a acquires the counting information stored in the memory 130 and generates coincidence counting information based on the acquired counting information. The coincidence counting information generating function 106a then arranges the generated coincidence counting information in approximately chronological order based on the detection time (T) and stores it in the memory 130. As a result, a chronological list of the coincidence counting information is stored in the memory 130.

[0028] The reconstruction processing function 106b reconstructs PET image data. For example, the reconstruction processing function 106b acquires a time series list of coincidence counting information stored in the memory 130 and reconstructs PET image data using the acquired time series list of coincidence counting information. The reconstruction processing function 106b then stores the reconstructed PET image data in the memory 130.

[0029] The system control function 106c controls the entire PET device 100 by controlling the gantry device 10 and the console device 20. For example, the system control function 106c controls imaging in the PET device 100.

[0030] The bed control function 106d controls the bed driving unit 105 to control the movement of the tabletop 103. The acquisition function 106e, the generation function 106f, the calculation function 106g, and the determination function 106h will be described later.

[0031] The display control function 106i displays various images and various information on the display 120. For example, the display control function 106i displays a PET image based on PET image data on the display 120. The display control function 106i also displays a GUI (Graphical User Interface) on the display 120 for receiving various instructions and settings from a user (operator) of the PET device 100.

[0032] The input interface 110 accepts various instructions and settings input from a user and outputs the accepted instructions and settings to the processing circuitry 106. For example, the input interface 110 converts the input operation accepted by the user into an electrical signal and transmits the electrical signal to the processing circuitry 106. For example, the input interface 110 may be realized by a trackball, a switch button, a mouse, a keyboard, a touchpad that performs input operations by touching the operation surface, a touchscreen that integrates a display screen and a touchpad, a non-contact input circuit using an optical sensor, an audio input circuit, or the like. Note that in this specification, the input interface 110 is not limited to an interface that includes physical operation components such as a mouse and a keyboard. For example, an electrical signal processing circuit that receives an electrical signal corresponding to an input operation from an external input device provided separately from the PET device 100 and transmits the electrical signal to the processing circuitry 106 is also included as an example of the input interface 110. The input interface 110 is, for example, an example of a receiving unit.

[0033] The display 120 is connected to the processing circuitry 106 and displays various types of information and images. For example, the display 120 converts information and data transmitted from the processing circuitry 106 into an electrical signal for display and outputs the signal. To give a specific example, under the control of the display control function 106i, the display 120 displays a PET image based on PET image data and a GUI for receiving various instructions and settings from a user. For example, the display 120 is realized by a liquid crystal monitor, a CRT (Cathode Ray Tube) monitor, a touch panel, or the like. The display 120 is an example of a display unit.

[0034] The memory 130 stores various data used in the PET device 100. The memory 130 is an example of a storage unit. The memory 130 is realized by, for example, a semiconductor memory element such as a random access memory (RAM) or a flash memory, a hard disk, an optical disk, etc. The memory 130 stores counting information, which is information in which a scintillator number (P), an energy value (E), and a detection time (T) are associated with each other, a time-series list of coincidence counting information in which sets of counting information are associated with coincidence numbers, which are serial numbers of the coincidence counting information, reconstructed PET image data, etc.

[0035] The memory 130 also stores a sensitivity map 130a. FIG. 2 is a diagram illustrating an example of the sensitivity map 130a according to the first embodiment. The sensitivity map 130a indicates the sensitivity of the PET detector 101 for each pixel in a region to be reconstructed (reconstruction region). The pixel size is, for example, 4 mm×4 mm, but is not limited to this. The sensitivity map 130a is obtained when the PET detector 101 is not malfunctioning and all channels of the PET detector 101 are normal channels rather than faulty channels. That is, the sensitivity map 130a is a sensitivity map of the PET detector 101 at a timing when all channels of the PET detector 101 are normal channels. The sensitivity map 130a is, for example, an example of a first sensitivity map. The timing when all channels of the PET detector 101 are normal channels is also an example of a first timing.

[0036] Here, an example of a method for generating the sensitivity map 130a will be described. The sensitivity map 130a may be generated by the PET device 100, or may be generated by a device other than the PET device 100. In the following description, a case will be described in which the generation function 106f of the PET device 100 generates the sensitivity map 130a.

[0037] 3 is a diagram for explaining an example of a method for generating a sensitivity map 130a according to the first embodiment. As shown in FIG. 3, the generation function 106f generates a model 30 of the PET detector 101 in which models of multiple detector modules are arranged in a ring shape.

[0038] In the example of FIG. 3 , the generation function 106f generates a model 30 of the PET detector 101 in which models of 48 detector modules are arranged in a ring shape. One detector module model includes at least one combination 31 of one scintillator model and one SiPM model. For example, the number and arrangement of combinations 31 included in the model 30 are preferably the same as the number and arrangement of combinations of one scintillator and one SiPM included in the PET detector 101. For example, the number of combinations 31 included in the model 30 is preferably between tens of thousands and hundreds of thousands, which is the same as the number of combinations of one scintillator and one SiPM included in the PET detector 101.

[0039] However, for convenience of explanation, a case will be described here where one detector module model has only one combination 31. That is, in the model 30 of the PET detector 101 shown in FIG. 3, 48 combinations 31 are arranged in a ring shape. In this case, the number of channels in the model 30 of the PET detector 101 is "48". Also, all channels are normal channels. Each of the 48 channels is identified by an identifier "1" to "48".

[0040] Then, the generation function 106f draws LORs (Lines of Response) 32 connecting two of the 48 channels for each combination of two channels, for all combinations. In the example of Fig. 3, the generation function 106f draws 1128 ((48 x 47) / 2) LORs 32. The pattern of LORs 32 drawn in this way is also called an LOR pattern.

[0041] As shown in Figure 3, the density of LOR32 is not uniform, and there are areas where the density of LOR32 is relatively sparse and areas where it is relatively dense. Pair annihilation gamma rays emitted from dense areas are detected more frequently than pair annihilation gamma rays emitted from sparse areas. In other words, pair annihilation is more likely to occur in dense areas than in sparse areas.

[0042] Then, the generation function 106f generates the sensitivity map 130a by converting the LOR pattern in the reconstruction region 33 of the model 30 into image data expressed in terms of LOR density for each pixel. That is, the generation function 106f calculates the LOR density for each pixel based on the LOR pattern in the reconstruction region 33 of the model 30, and generates image data in which each pixel indicates the LOR density as the sensitivity map 130a. That is, the sensitivity map 130a is data in which the LOR density is registered as sensitivity for each pixel.

[0043] The sensitivity map 130a generated in this manner is stored in advance in the memory 130, and is used to normalize non-uniformity in sensitivity due to the arrangement geometry of channels (combinations of one scintillator model and one SiPM model) in the reconstruction processing of PET image data by the reconstruction processing function 106b. For example, the reconstruction processing function 106b uses the sensitivity map 130a to reconstruct PET image data according to the following equation (1):

[0044]

number

[0045] For example, equation (1) is a function of the pixel value λ of pixel j. k , the pixel value S of pixel j in the sensitivity map S j Reciprocal of (1 / S j ) to obtain the pixel value λ of pixel j of the PET image data finally obtained as a result of reconstruction. k+1 Equation (1) is a recurrence formula for finding the probability of detection C im is also called the system matrix and represents the probability that a photon generated from the mth pixel is counted by the ith LOR. iis the number of events counted at the i-th LOR obtained in the actual collection, i.e., the actual measured value. In the ML-EM method, some radioactivity distribution is assumed as the initial value. In other words, λ(k=0) is assumed. In the ML-EM method, for the assumed distribution, the number of events that should be counted at the i-th LOR is calculated as C im and the contribution of pixel m. In the ML-EM method, the ratio of the calculated number of events to the number of actually measured events (the number of actually measured events indicated by the numerator in formula (1)) is taken for each LOR. If this ratio is 1 for all LORs, the assumed pixel values ​​will reproduce the actually measured events. In this case, the assumed pixel values ​​will reproduce the actual radioactivity distribution of the subject. The operation of formula (1) is as follows: j If there is no channel defect, Σ i C ij Since it is defined as such, the factors on the right side after Σ are S j and the k-th pixel value and the k+1-th pixel value become equal. In other words, the calculation has converged.

[0046] In this way, the reconstruction processing function 106b reconstructs PET image data using the inverse of the sensitivity for each pixel indicated by the sensitivity map 130a.

[0047] An example of the configuration of the PET device 100 according to this embodiment has been described above. For example, the PET device 100 is installed in a medical institution such as a hospital or a clinic, and is used for various image diagnoses using PET image data generated by the PET device 100, with a patient who is admitted to or visits the medical institution as a subject P. In this embodiment, the PET device 100 executes various processes described below so that the influence of a poor channel on image quality (or an image) can be quantitatively and objectively evaluated.

[0048] Next, an example of evaluation processing executed by the PET device 100 will be described. Fig. 4 is a flowchart showing an example of the flow of the evaluation processing executed by the PET device 100 according to the first embodiment. The PET device 100 executes the evaluation processing when, for example, a user inputs an instruction (execution instruction) to execute the evaluation processing shown in Fig. 4 to the processing circuitry 106 via the input interface 110. Here, the user inputs the execution instruction to the processing circuitry 106 via the input interface 110 to determine whether the influence of a defective channel on image quality satisfies a predetermined quality standard (criteria), for example, after installing the PET device 100 in an imaging room of a medical institution, after regular maintenance of the PET device 100, after repair of the PET device 100, or after emergency response in the event of a malfunction of the PET device 100. In addition to the above-mentioned timings, the user inputs an execution instruction to the processing circuit 106 via the input interface 110 to determine whether the effect of a defective channel on image quality satisfies a predetermined quality standard, for example, during an acceptance test of the PET detector 101 alone, a pre-shipment test of the PET detector 101 alone, a pre-start-up inspection of the PET device 100, or a decision on a repair policy for the PET device 100.

[0049] (Step S101) As shown in FIG. 4, in step S101, the acquisition function 106e acquires failure information (failure information), which is information about a failed channel. Note that, at the timing when the evaluation process is executed, at least one failed channel may occur among the channels of the PET detector 101, or all channels may be normal. The following description will be given taking as an example a case where at least one failed channel occurs among the channels of the PET detector 101 at the timing when the evaluation process is executed. The failure information about the failed channel may include, for example, the location and number of the failed channel. For example, the acquisition function 106e performs a health check to check whether the PET detector 101 is operating normally, and acquires the failure information about the failed channel. The acquisition function 106e then stores the acquired failure information in the memory 130. Thus, in step S101, the acquisition function 106e acquires the failure information about the failed channel of the PET detector 101 at a timing when a failed channel occurs after a timing when all channels of the PET detector 101 corresponding to the sensitivity map 130a stored in the memory 130 are normal. The timing when a defective channel occurs is an example of the second timing.

[0050] (Step S102) Next, in step S102, the acquisition function 106e acquires from the memory 130 the sensitivity map 130a, which is the sensitivity map of the PET detector 101 when all channels of the PET detector 101 are normal channels.

[0051] (Step S103) Next, in step S103, the generation function 106f generates a sensitivity map of the PET detector 101 at the timing when a faulty channel occurs, based on the fault information. That is, in step S103, a sensitivity map reflecting the faulty channel is generated in consideration of the faulty channel that actually occurred. The sensitivity map generated in step S103 is, for example, an example of a second sensitivity map. A specific example of the processing in step S103 will be described below. For example, the generation function 106f first acquires fault information from the memory 130.

[0052] Then, the generation function 106f generates a sensitivity map that takes into account the faulty channels indicated by the fault information, using a method similar to the method for generating the sensitivity map 130a described with reference to Fig. 3. For example, the generation function 106f identifies a plurality of normal channels from among the 48 channels of the model 30 of the PET detector 101, excluding the faulty channels identified by the fault information.

[0053] The generation function 106f then draws an LOR 32 connecting two of the identified normal channels for each combination of two channels, for all combinations of two channels. The generation function 106f then generates a sensitivity map by converting the LOR pattern in the reconstruction region 33 of the model 30 into image data represented by the LOR density for each pixel. That is, the generation function 106f calculates the LOR density for each pixel as the sensitivity of the PET detector 101 based on the LOR pattern in the reconstruction region 33 of the model 30, and generates image data in which each pixel indicates the sensitivity as a sensitivity map. Using the method described above, in step S103, the generation function 106f generates a sensitivity map in which the sensitivity is registered for each pixel, taking into account the defective channels.

[0054] The arrangement and number of pixels constituting the sensitivity map generated in step S103 are the same as the arrangement and number of pixels constituting the sensitivity map 130a. Therefore, each of the pixels constituting the sensitivity map generated in step S103 corresponds to each of the pixels constituting the sensitivity map 130a.

[0055] Here, the sensitivity map 130a used when reconstructing PET image data is a sensitivity map generated in advance and stored in the memory 130 when all channels are normal channels. Therefore, a deviation of the sensitivity indicated by the sensitivity map generated in step S103 from the sensitivity indicated by the sensitivity map 130a appears in the sensitivity map generated in step S103 as a deviation of the sensitivity due to a defective channel.

[0056] It will be explained that even if the number of defective channels is the same, if the arrangement of the defective channels is different, the sensitivity maps generated in step S103 may be different. Figures 5 and 6 are diagrams showing examples of the sensitivity maps generated in step S103 according to the first embodiment.

[0057] The sensitivity map 35 shown in Fig. 5 is a sensitivity map generated in step S103 when the channel indicated by the identifier "1" and the channel indicated by the identifier "5" are faulty channels. On the other hand, the sensitivity map 35 shown in Fig. 6 is a sensitivity map generated in step S103 when the channel indicated by the identifier "1" and the channel indicated by the identifier "10" are faulty channels. As can be seen from a comparison of the sensitivity map 35 shown in Fig. 5 and the sensitivity map 35 shown in Fig. 6, when the number of faulty channels is the same but the arrangement of the faulty channels is different, the two sensitivity maps 35 are different. One reason why the two sensitivity maps 35 are different in this way is that the relative positions of the faulty channels are different between the sensitivity map 35 shown in Fig. 5 and the sensitivity map 35 shown in Fig. 6.

[0058] (Step S104) Next, in step S104, the calculation function 106g calculates an index value relating to the influence of a poor channel on image quality based on the sensitivity map 130a and the sensitivity map generated in step S103. A specific example of the processing in step S104 will be described below. For example, the calculation function 106g generates image data by dividing the sensitivity map generated in step S103 by the sensitivity map 130a. The arrangement and number of pixels constituting the image data generated by the calculation function 106g are the same as the arrangement and number of pixels constituting the sensitivity map 130a. Furthermore, the pixel value of each of the pixels constituting the image data generated by the calculation function 106g is a value obtained by dividing the sensitivity indicated by each of the pixels constituting the sensitivity map generated in step S103 by the sensitivity indicated by each of the pixels constituting the sensitivity map 130a. That is, the calculation function 106g calculates an index value for each pixel based on the ratio between the sensitivity map 130a and the sensitivity map generated in step S103.

[0059] In step S104, the calculation function 106g generates such image data and calculates the pixel values ​​of the pixels constituting the image data as an index value relating to the degree of influence of the defective channel on the image quality. That is, the calculation function 106g evaluates the sensitivity deviation due to the defective channel for each pixel and calculates the evaluation result for each pixel as an index value.

[0060] Here, the index value is a value in the range of 0 to 1. The closer the index value is to 1, the smaller the impact of a bad channel on image quality. On the other hand, the closer the index value is to 0, the greater the impact of a bad channel on image quality.

[0061] The calculation function 106g may calculate an index value for each pixel based on the difference between the sensitivity map 130a and the sensitivity map generated in step S103. For example, the calculation function 106g may generate image data by subtracting the sensitivity map generated in step S103 from the sensitivity map 130a. In this case, the pixel value of each of the multiple pixels constituting the image data generated by the calculation function 106g is a value obtained by subtracting the sensitivity indicated by each of the multiple pixels constituting the sensitivity map generated in step S103 from the sensitivity indicated by each of the multiple pixels constituting the sensitivity map 130a. In this case, the index value is a value equal to or greater than 0. The closer the index value is to 0, the smaller the influence of the poor channel on image quality. On the other hand, the larger the index value, the greater the influence of the poor channel on image quality.

[0062] (Step S105) Next, in step S105, the determination function 106h determines whether the sensitivity of the PET detector 101 satisfies a predetermined standard based on the index value. A specific example of the processing in step S105 will be described below. For example, a case will be described in which the required specifications of the PET device 100 are that it is possible to quantitatively detect accumulations of a pixel size of 4 mm x 4 mm and to correctly perform SUV (Standardized Uptake Value) evaluation.

[0063] In this case, for example, a necessary condition for the PET detector 101 is set such that all index values ​​calculated for each pixel in step S104 are greater than a threshold value α. While the threshold value α is, for example, 0.9, the threshold value α is not limited to this. Therefore, in this case, in step S105, the determination function 106h determines whether the sensitivity of the PET detector 101 satisfies a predetermined standard by determining whether all index values ​​calculated for each pixel are greater than the threshold value α. For example, if the determination function 106h determines that all index values ​​calculated for each pixel are greater than the threshold value α, it determines that the sensitivity of the PET detector 101 satisfies the predetermined standard. If the determination function 106h determines that at least one index value among all index values ​​calculated for each pixel is equal to or less than the threshold value α, it determines that the sensitivity of the PET detector 101 does not satisfy the predetermined standard.

[0064] If the determination function 106h determines that the sensitivity of the PET detector 101 satisfies the predetermined standard (step S105: Yes), the process proceeds to step S106. On the other hand, if the determination function 106h determines that the sensitivity of the PET detector 101 does not satisfy the predetermined standard (step S105: No), the process proceeds to step S107.

[0065] (Step S106) In step S106, the display control function 106i displays a message indicating the determination result by the determination function 106h, i.e., a message indicating that the sensitivity of the PET detector 101 meets a predetermined standard, on the display 120. For example, in step S106, the display control function 106i displays the character string "The sensitivity of the PET detector 101 meets the standard" as such a message on the display 120. Then, the display control function 106i ends the processing shown in FIG. 4.

[0066] (Step S107) In step S107, the display control function 106i displays a message indicating the determination result by the determination function 106h, i.e., a message indicating that the sensitivity of the PET detector 101 does not satisfy a predetermined standard, on the display 120. For example, in step S107, the display control function 106i displays the character string "The sensitivity of the PET detector 101 does not satisfy the standard" as such a message on the display 120. Then, the display control function 106i ends the processing shown in FIG.

[0067] As described above, in steps S106 and S107, the display control function 106i causes the display 120 to display the determination result by the determination function 106h.

[0068] 4 is limited to geometric calculations using the results of a health check that checks for defective channels in the PET detector 101. Furthermore, the evaluation process is limited to calculations at the PET detector 101 level, and is not dependent on the magnitude of the collected count, the strength of statistical noise due to the magnitude of the collected count, the experience of the user who interprets the images, etc. Therefore, the PET device 100 according to the first embodiment can quantitatively and objectively evaluate the influence of defective channels on image quality and diagnostic performance in the PET detector 101 having multiple channels.

[0069] Furthermore, the PET device 100 allows the user to confirm a message indicating that the sensitivity of the PET detector 101 meets a predetermined standard or a message indicating that the sensitivity of the PET detector 101 does not meet the predetermined standard. Then, the PET device 100 allows the user who has confirmed the message to correctly determine how to deal with the PET device 100. As a result, the user can take measures such as ensuring the safety and performance of the PET device 100 in clinical practice. Furthermore, the user can allocate appropriate costs in manufacturing and servicing.

[0070] The above has described the PET device 100 according to the first embodiment. As described above, the PET device 100 according to the first embodiment can quantitatively and objectively evaluate the influence of a poor channel on image quality (or an image).

[0071] 4, the reconstruction processing function 106b may reconstruct PET image data using, instead of the sensitivity map 130a, the sensitivity map generated in step S103 which more accurately reflects the channel state of the PET detector 101 than the sensitivity map 130a. For example, the reconstruction processing function 106b may reconstruct PET image data in accordance with the following equation (2) based on the sensitivity map generated in step S103.

[0072]

number

[0073] In addition, in equation (2), S' j is the pixel value of pixel j in the sensitivity map S' generated in step S103, which reflects the channel states of the PET detector 101 more accurately than the sensitivity map 130a. In this way, the reconstruction processing function 106b reconstructs the PET image data using the inverse of the sensitivity for each pixel indicated by the sensitivity map generated in step S103. As a result, the reconstruction processing function 106b reconstructs the PET image data using a sensitivity map which reflects the channel states of the PET detector 101 more accurately than the sensitivity map 130a, and therefore the PET image data can be reconstructed with even greater accuracy.

[0074] If it is determined in step S105 that the sensitivity of the PET detector 101 does not satisfy a predetermined standard (step S105: No), the reconstruction processing function 106b may reconstruct PET image data using the sensitivity map generated in step S103. When such a PET device 100 is operated in a state in which it has a defective channel, image uniformity is restored by renormalization that reflects the defective channel, making it easier to interpret the PET image.

[0075] (Modification of the first embodiment) Here, after the PET device 100 executes the evaluation process (first evaluation process) shown in Fig. 4, the user may replace the defective channel with a normal channel or change the position of the defective channel. In this case, the PET device 100 may execute the evaluation process (second evaluation process) again to generate, in step S103 of the second evaluation process, a sensitivity map having improved sensitivity compared to the sensitivity map generated in step S103 of the first evaluation process. Therefore, such a modification will be described as a modification of the first embodiment. Note that in the description of the modification of the first embodiment, differences from the first embodiment will be mainly described, and a description of the same configuration as the first embodiment may be omitted.

[0076] In the modification of the first embodiment, when the PET device 100 executes the first evaluation process and a message indicating that the sensitivity of the PET detector 101 does not satisfy the predetermined standard is displayed in step S107, the user checks the message and understands that the sensitivity of the PET detector 101 does not satisfy the predetermined standard.

[0077] In this case, the user can take various measures. For example, as a first measure, the user replaces the combination of scintillator and SiPM corresponding to the faulty channel with a combination of normal scintillator and SiPM that is not faulty, in order to make all channels normal. That is, the user replaces the faulty channel with a normal channel. Alternatively, as a second measure, the user changes the position of the combination of scintillator and SiPM corresponding to the faulty channel to another position, and changes the position of the normal combination of scintillator and SiPM that was placed at the other position to the position of the faulty channel before the change, so that some channels remain faulty channels but the sensitivity of the PET detector 101 satisfies a predetermined standard. That is, the user exchanges the position of the faulty channel with the position of the normal channel.

[0078] Then, the PET device 100 executes the second evaluation process. For example, when an execution instruction to execute the second evaluation process is input to the processing circuit 106 from the user via the input interface 110, the PET device 100 executes the second evaluation process. Then, in step S106 of the second evaluation process, if a message indicating that the sensitivity of the PET detector 101 meets a predetermined standard is displayed, the user confirms the message. Then, the user understands that the sensitivity of the PET detector 101 meets the predetermined standard. Then, after executing the second evaluation process, the reconstruction processing function 106b reconstructs PET image data according to the above equation (2) using the sensitivity map generated in step S103 of the second evaluation process. In this case, in equation (2), S' j is the pixel value of pixel j in the sensitivity map S' generated in step S103 of the second evaluation processing. Here, the sensitivity map generated in step S103 of the second evaluation processing is a sensitivity map with improved sensitivity compared to the sensitivity map generated in step S103 of the first evaluation processing. Therefore, in the modification of the first embodiment, PET image data with good image quality can be reconstructed.

[0079] Moreover, the sensitivity map generated in step S103 of the second evaluation process is the sensitivity map of the PET detector 101 after the defective channel has been changed to a normal channel that is not defective. Alternatively, the sensitivity map generated in step S103 of the second evaluation process is the sensitivity map of the PET detector 101 after the position of the defective channel has been changed to a position different from the position of the defective channel at the timing of executing the first evaluation process. The timing of executing the first evaluation process is, for example, an example of the first timing.

[0080] As described above, in the modified example of the first embodiment, in step S101 of the first evaluation process, the acquisition function 106e acquires failure information related to a faulty channel of the PET detector 101 at the timing when the first evaluation process is executed.

[0081] Then, in step S103 of the first evaluation process, the generation function 106f generates a sensitivity map of the PET detector 101 at the timing when a faulty channel occurs, based on the failure information acquired in step S101 of the first evaluation process. Here, the sensitivity map generated in step S103 of the first evaluation process is, for example, an example of a first sensitivity map.

[0082] Then, in step S103 of the second evaluation process, the generation function 106f generates a sensitivity map having improved sensitivity compared to the sensitivity map generated in step S103 of the first evaluation process, based on the failure information acquired in step S101 of the second evaluation process. Here, the sensitivity map generated in step S103 of the second evaluation process is, for example, an example of the second sensitivity map.

[0083] The PET device 100 according to the modification of the first embodiment has been described above. As described above, the PET device 100 according to the modification of the first embodiment can reconstruct PET image data with good image quality. Furthermore, the PET device 100 according to the modification of the first embodiment performs the same processing as the processing performed by the PET device 100 according to the first embodiment, and therefore achieves the same effects as the PET device 100 according to the first embodiment.

[0084] (Second embodiment) In the first embodiment, a case where the PET device 100 displays the determination result has been described. However, the PET device 100 may display, in addition to the determination result, an image based on image data generated when calculating the index value in step S104. Therefore, such an embodiment will be described as a second embodiment. Note that in the description of the second embodiment, differences from the first embodiment will be mainly described, and a description of the same configuration as the first embodiment may be omitted.

[0085] An example of evaluation processing executed by the PET device 100 according to the second embodiment will be described. FIG. 7 is a flowchart showing an example of the flow of evaluation processing executed by the PET device 100 according to the second embodiment. In the second embodiment, the PET device 100 executes the evaluation processing when an execution instruction for executing the evaluation processing is input to the processing circuitry 106 by the user via the input interface 110 at the same timing as in the first embodiment, for example. Note that the processing in steps S101 to S104 and step S105 according to the second embodiment shown in FIG. 7 is the same as the processing in steps S101 to S104 and step S105 according to the first embodiment shown in FIG. 4. Therefore, a description of the processing in steps S101 to S104 and step S105 according to the second embodiment will be omitted.

[0086] However, in step S104 according to the second embodiment, the calculation function 106g stores image data generated when calculating the index value in the memory 130. This image data is obtained by dividing the sensitivity map generated in step S103 by the sensitivity map 130a, or by subtracting the sensitivity map generated in step S103 from the sensitivity map 130a.

[0087] (Step S201) Then, as shown in FIG. 7, in step S201 between step S104 and step S105, the acquisition function 106e acquires the image data stored in the memory 130 in the previous step S104.

[0088] If the determination function 106h determines that the sensitivity of the PET detector 101 satisfies the predetermined standard (step S105: Yes), the display control function 106i proceeds to step S202. On the other hand, if the determination function 106h determines that the sensitivity of the PET detector 101 does not satisfy the predetermined standard (step S105: No), the display control function 106i proceeds to step S203.

[0089] (Step S202) In step S202, the display control function 106i displays a message indicating that the sensitivity of the PET detector 101 satisfies a predetermined standard, and an image based on the image data acquired in step S201, on the display 120. Then, the display control function 106i ends the evaluation process shown in FIG.

[0090] (Step S203) In step S203, the display control function 106i displays a message indicating that the sensitivity of the PET detector 101 does not satisfy a predetermined standard, and an image based on the image data acquired in step S201, on the display 120. Then, the display control function 106i ends the evaluation process shown in FIG.

[0091] The PET device 100 according to the second embodiment has been described above. The images displayed in steps S202 and S203 show a two-dimensional distribution of index values. Therefore, the PET device 100 according to the second embodiment allows the user to visually understand the state of the PET detector 101 of the PET device 100. Therefore, the PET device 100 according to the second embodiment allows the user to understand the detailed state of the PET detector 101. Furthermore, the PET device 100 according to the second embodiment allows the user to fully understand the state of the PET detector 101, thereby allowing the user to operate the PET device 100 in a state where the user is satisfied.

[0092] In addition, during pre-work inspections, macro performance such as the energy resolution and time resolution of the PET detector 101 is usually checked, but the PET device 100 of the second embodiment allows the user to understand not only such macro performance but also micro performance such as the two-dimensional distribution of index values.

[0093] Furthermore, the PET device 100 according to the second embodiment performs the same processing as the processing executed by the PET device 100 according to the first embodiment, and therefore provides the same effects as the PET device 100 according to the first embodiment.

[0094] (Third embodiment) In the first embodiment, a case has been described in which, when it is determined that the sensitivity of the PET detector 101 does not satisfy the predetermined standard (step S105: No), the PET device 100 displays a message indicating that the sensitivity of the PET detector 101 does not satisfy the predetermined standard. In the second embodiment, a case has been described in which, when it is determined that the sensitivity of the PET detector 101 does not satisfy the predetermined standard (step S105: No), the PET device 100 displays a message indicating that the sensitivity of the PET detector 101 does not satisfy the predetermined standard and an image based on the image data acquired in step S201. However, when it is determined that the sensitivity of the PET detector 101 does not satisfy the predetermined standard (step S105: No), the PET device 100 may search for an arrangement of defective channels that would result in the sensitivity of the PET device 100 satisfying the predetermined standard, and display information indicating the arrangement of the defective channels on the display 120. Therefore, such an embodiment will be described as a third embodiment.

[0095] In the following, a PET device 100 according to the third embodiment will be described, which executes the processes in steps S301 and S302 described below instead of the process in step S107 in the first embodiment. However, the PET device 100 according to the third embodiment may execute the processes in steps S301 and S302 instead of the process in step S203 in the second embodiment.

[0096] In addition, in the description of the third embodiment, differences from the first embodiment will be mainly described, and a description of the same configuration as the first embodiment may be omitted.

[0097] An example of evaluation processing executed by the PET device 100 according to the third embodiment will be described. FIG. 8 is a flowchart showing an example of the flow of evaluation processing executed by the PET device 100 according to the third embodiment. In the third embodiment, the PET device 100 executes the evaluation processing when an execution instruction for executing the evaluation processing is input to the processing circuit 106 by the user via the input interface 110 at the same timing as in the first embodiment, for example. Note that the processing in steps S101 to S106 according to the third embodiment shown in FIG. 8 is the same as the processing in steps S101 to S106 according to the first embodiment shown in FIG. 4. Therefore, a description of the processing in steps S101 to S106 according to the third embodiment will be omitted.

[0098] As shown in FIG. 8, in the third embodiment, if it is determined that the sensitivity of the PET detector 101 does not satisfy a predetermined standard (step S105: No), the processing circuitry 106 proceeds to step S301.

[0099] (Step S301) In step S301, the processing circuitry 106 executes a specification process for specifying the position (arrangement) of a defective channel when the sensitivity of the PET detector 101 satisfies a predetermined standard. Fig. 9 is a flowchart showing an example of the flow of the specification process executed in step S301 by the PET device 100 according to the third embodiment.

[0100] (Step S303) As shown in FIG. 9, in step S303, the generating function 106f sets the variable K to "1".

[0101] (Step S304) Next, in step S304, the generation function 106f changes the position (placement) of the faulty channel indicated by the failure information in the model 30 of the PET detector 101. For example, the generation function 106f randomly changes the position of the faulty channel. Specifically, for example, the generation function 106f randomly generates random numbers from information such as the current time. Here, each of the generated random numbers corresponds one-to-one to each of the multiple positions at which the channel can be placed. The generation function 106f then identifies the position corresponding to the generated random number. The generation function 106f then sets the faulty channel at the identified position. This changes the position of the faulty channel. The generation function 106f also sets a normal channel to a channel at which no faulty channel is set, among the multiple positions at which channels can be placed. In this way, the generation function 106f changes the position of the faulty channel in the model of the PET detector 101 without changing the number of faulty channels.

[0102] (Step S305) Next, in step S305, the generation function 106f draws LORs 32 connecting two channels from among multiple normal channels excluding the defective channel for all combinations of two channels. Then, the generation function 106f generates a sensitivity map by converting the LOR pattern in the reconstruction region 33 of the model 30 into image data represented by the LOR density for each pixel. That is, the generation function 106f calculates the LOR density for each pixel as the sensitivity of the PET detector 101 based on the LOR pattern in the reconstruction region 33 of the model 30, and generates image data indicating the sensitivity of each pixel as a sensitivity map.

[0103] In this way, when it is determined that the sensitivity of the PET detector 101 does not satisfy the predetermined standard (step S105: No), the generation function 106f changes the position of the defective channel in the PET detector 101 in step S303, and generates a sensitivity map of the PET detector 101 based on the changed position of the defective channel in step S305. The sensitivity map generated in step S305 is, for example, an example of a third sensitivity map.

[0104] (Step S306) Next, in step S306, the calculation function 106g calculates an index value related to the influence of the poor channel on image quality based on the sensitivity map 130a and the sensitivity map generated in step S305. For example, the calculation function 106g generates image data by dividing the sensitivity map generated in step S305 by the sensitivity map 130a, in a manner similar to the method for generating image data in the previous step S104. That is, the calculation function 106g calculates an index value for each pixel based on the ratio between the sensitivity map 130a and the sensitivity map generated in step S305.

[0105] Note that the calculation function 106g may calculate an index value for each pixel based on the difference between the sensitivity map 130a and the sensitivity map generated in step S305. For example, the calculation function 106g may generate image data by subtracting the sensitivity map generated in step S305 from the sensitivity map 130a in a manner similar to the method of generating image data in the previous step S104 by subtracting the sensitivity map generated in step S103 from the sensitivity map 130a.

[0106] Thus, in step S306, the calculation function 106g calculates an index value relating to the degree of influence of the defective channel whose position has been changed on the image quality, based on the sensitivity map 130a and the sensitivity map generated in step S305.

[0107] (Step S307) Next, in step S307, the determination function 106h determines, based on the index value, whether or not the sensitivity of the PET detector 101 satisfies a predetermined standard. For example, the determination function 106h determines, based on the index value, whether or not the sensitivity of the PET detector 101 satisfies a predetermined standard, in a manner similar to the method for determining, in the previous step S105, whether or not the sensitivity of the PET detector 101 satisfies a predetermined standard.

[0108] In this way, the determination function 106h determines whether the sensitivity of the PET detector 101 satisfies a predetermined standard based on the index value calculated based on the sensitivity map 130a and the sensitivity map generated in step S305.

[0109] If it is determined that the sensitivity of the PET detector 101 does not satisfy the predetermined standard (step S307: No), the generation function 106f proceeds to step S308. On the other hand, if it is determined that the sensitivity of the PET detector 101 satisfies the predetermined standard (step S307: Yes), the generation function 106f proceeds to step S310.

[0110] (Step S308) In step S308, the generating function 106f increments the value of the variable K by one.

[0111] (Step S309) Next, in step S309, the generation function 106f determines whether the value of the variable K is equal to or greater than a threshold value Th1. The threshold value Th1 is a positive integer value that determines the upper limit of the number of times that the processes in steps S303 to S308 are repeated. In this embodiment, the processes in steps S303 to S308 are repeated a maximum of (threshold value Th1-1) times.

[0112] If the value of the variable K is equal to or greater than the threshold Th1 (step S309: Yes), the generating function 106f returns. On the other hand, if the value of the variable K is less than the threshold Th1 (step S309: No), the generating function 106f returns to the previous step S304 and executes each process from step S304 onwards.

[0113] (Step S310) In step S310, the generation function 106f stores information indicating the position (arrangement) of the bad channel when the sensitivity of the PET detector 101 satisfies a predetermined standard in the memory 130. For example, the generation function 106f generates information indicating the position of the bad channel changed in the most recent step S304, and stores the generated information in the memory 130. Then, the generation function 106f returns.

[0114] Here, if the determination in step S309 is affirmative and the process is returned, memory 130 does not store information indicating the position of a defective channel when the sensitivity of PET detector 101 satisfies the predetermined standard. On the other hand, if the determination in step S307 is affirmative and the process is returned, memory 130 stores information indicating the position of a defective channel when the sensitivity of PET detector 101 satisfies the predetermined standard.

[0115] (Step S302) Returning to the description of Fig. 8, if information indicating the positions of defective channels when the sensitivity of the PET detector 101 satisfies the predetermined standard is stored in the memory 130, in step S302, as shown in Fig. 8, the display control function 106i causes the display 120 to display a message indicating that the sensitivity of the PET detector 101 does not satisfy the predetermined standard, and information indicating the positions of defective channels when the sensitivity of the PET detector 101 satisfies the predetermined standard. Then, the display control function 106i ends the evaluation process shown in Fig. 8.

[0116] In this way, the display control function 106i causes the display 120 to display information indicating the position of the defective channel in the PET detector 101 after the change when the judgment function 106h judges that the sensitivity of the PET detector 101 satisfies the predetermined standard.

[0117] On the other hand, if the memory 130 does not store information indicating the arrangement of defective channels when the sensitivity of the PET detector 101 satisfies the predetermined standard, in step S302, the display control function 106i displays a message indicating that the sensitivity of the PET detector 101 does not satisfy the predetermined standard on the display 120. Then, the display control function 106i ends the evaluation process shown in FIG.

[0118] The PET device 100 according to the third embodiment has been described above. As a measure to deal with a defective channel, the PET device 100 according to the third embodiment allows the user to understand the arrangement of the defective channel that reduces the impact on image quality. As a result, clinical downtime due to repair work by the user can be reduced.

[0119] Furthermore, the PET device 100 according to the third embodiment performs processing similar to the processing performed by the PET device 100 according to the first or second embodiment, and therefore achieves the same effects as the PET device 100 according to the first or second embodiment.

[0120] 8, the reconstruction processing function 106b may reconstruct PET image data using the sensitivity map generated in step S305 when the determination function 106h determines that the sensitivity of the PET detector 101 satisfies the predetermined standard, instead of the sensitivity map 130a. For example, the reconstruction processing function 106b may reconstruct PET image data according to the above formula (2) based on the sensitivity map generated in step S305 when the determination function 106h determines that the sensitivity of the PET detector 101 satisfies the predetermined standard. In this case, in formula (2), S' j is the pixel value of pixel j in the sensitivity map S' generated in step S305 when it is determined that the sensitivity of the PET detector 101 satisfies the predetermined criterion.

[0121] (Fourth embodiment) Here, in the above-described embodiment and modified examples, the PET device 100 may further determine the number of allowable defective channels in the PET detector 101. Such an allowable number of defective channels is used, for example, in an acceptance test or the like performed when receiving the PET detector 101. Therefore, such an embodiment will be described as a fourth embodiment.

[0122] The PET device 100 according to the fourth embodiment executes various processes described in any of the above-mentioned embodiments or modifications. The PET device 100 according to the fourth embodiment further executes a process for determining the number of allowable defective channels. An example of the process for determining the number of allowable defective channels executed by the PET device 100 according to the fourth embodiment will be described. FIG. 10 is a flowchart showing an example of the flow of the process for determining the number of allowable defective channels executed by the PET device 100 according to the fourth embodiment. In the fourth embodiment, the PET device 100 executes the process for determining the number of allowable defective channels when, for example, a user inputs an instruction to execute the process for determining the number of allowable defective channels to the processing circuitry 106 via the input interface 110.

[0123] (Step S401) 10, in step S401, the acquisition function 106e sets "1" to each of variables n and m. The variable n indicates the number of bad channels in a model of one detector module. For example, the model 30 of the PET detector 101 includes models of 48 detector modules. Therefore, the model 30 of the PET detector 101 has (48×n) bad channels.

[0124] (Step S402) Next, in step S402, the acquisition function 106e acquires from the memory 130 the sensitivity map 130a of the PET detector 101 when all channels of the PET detector 101 are normal channels.

[0125] (Step S403) Next, in step S403, the generation function 106f randomly determines the positions (arrangements) of each of the (48×n) failed channels in the model 30 of the PET detector 101. Specifically, for example, the generation function 106f randomly generates random numbers based on information about the current time and an identifier for identifying the detector module. Here, each of the generated random numbers corresponds one-to-one to each of the multiple positions where the channels can be arranged. The generation function 106f then identifies the positions of each of the (48×n) failed channels, which correspond to the generated random numbers. The generation function 106f then sets the failed channels at the positions of the identified (48×n) failed channels. This determines the positions of the failed channels. The generation function 106f also sets normal channels at positions where no failed channels are set, among the multiple positions where channels can be arranged. In this way, the generation function 106f determines the positions of the failed channels in the model of the PET detector 101 without changing the number of failed channels.

[0126] (Step S404) Next, in step S404, the generation function 106f generates a sensitivity map based on the position of the defective channel. An example of the processing in step S404 will be described below. First, the generation function 106f draws LORs 32 connecting two channels selected from multiple normal channels other than the defective channel for each combination of two channels. Then, the generation function 106f generates a sensitivity map by converting the LOR pattern in the reconstruction region 33 of the model 30 into image data represented by the LOR density for each pixel. That is, the generation function 106f calculates the LOR density for each pixel as the sensitivity of the PET detector 101 based on the LOR pattern in the reconstruction region 33 of the model 30, and generates image data indicating the sensitivity of each pixel as a sensitivity map.

[0127] (Step S405) Next, in step S405, the calculation function 106g calculates an index value related to the influence of the poor channel on image quality based on the sensitivity map 130a and the sensitivity map generated in step S404. For example, the calculation function 106g generates image data by dividing the sensitivity map generated in step S404 by the sensitivity map 130a, in a manner similar to the method for generating image data in the previous step S104. That is, the calculation function 106g calculates an index value for each pixel based on the ratio between the sensitivity map 130a and the sensitivity map generated in step S404.

[0128] Note that the calculation function 106g may calculate an index value for each pixel based on the difference between the sensitivity map 130a and the sensitivity map generated in step S404. For example, the calculation function 106g may generate image data by subtracting the sensitivity map generated in step S404 from the sensitivity map 130a in the same manner as the method of generating image data in the previous step S104 by subtracting the sensitivity map generated in step S103 from the sensitivity map 130a.

[0129] (Step S406) Next, in step S406, the determination function 106h determines, based on the index value, whether or not the sensitivity of the PET detector 101 satisfies a predetermined standard. For example, the determination function 106h determines, based on the index value, whether or not the sensitivity of the PET detector 101 satisfies a predetermined standard, in the same manner as in the previous step S105, whether or not the sensitivity of the PET detector 101 satisfies a predetermined standard.

[0130] If it is determined that the sensitivity of the PET detector 101 satisfies the predetermined standard (step S406: Yes), the generation function 106f proceeds to step S408.

[0131] (Step S408) In step S408, the generating function 106f determines whether the value of the variable m is equal to or greater than the threshold Th2. If the value of the variable m is equal to or greater than the threshold Th2 (step S408: Yes), the generating function 106f proceeds to step S409. On the other hand, if the value of the variable m is less than the threshold Th2 (step S408: No), the generating function 106f proceeds to step S410.

[0132] (Step S409) In step S409, the generating function 106f increments the value of the variable n by 1 and sets the value of the variable m to 1. Then, the generating function 106f returns to step S403 and executes the processes from step S403 onwards.

[0133] (Step S410) In step S410, the generating function 106f increments the value of the variable m by 1. Then, the generating function 106f returns to step S403 and executes the processes from step S403 onwards.

[0134] Therefore, the processing in steps S403 to S406 and S408 to S410 is repeated a maximum of (Th2-1) times while the value of the variable n remains constant.

[0135] Thus, in steps S403 to S406 and S408 to S410, the processing circuit 106 according to the fourth embodiment changes the positions of the defective channels in the PET detector 101 or in the multiple detector modules included in the PET detector 101, and generates a sensitivity map of the PET detector 101 based on the changed positions of the defective channels multiple times, thereby repeatedly performing the process of generating multiple sensitivity maps corresponding to multiple defective channel arrangement patterns in which the positions of the defective channels differ from one another while increasing the number of defective channels, until the sensitivity of the PET detector 101 no longer satisfies a predetermined standard in at least one of the multiple defective channel arrangement patterns.

[0136] On the other hand, if it is determined that the sensitivity of the PET detector 101 does not satisfy the predetermined standard (step S406: No), the generation function 106f proceeds to step S407.

[0137] (Step S407) In step S407, the determination function 106h determines the number of allowable defective channels in one detector module to be (n-1), and determines the number of allowable defective channels in the entire PET detector 101 to be (48×(n-1)). Then, the determination function 106h ends the process of determining the number of allowable defective channels shown in FIG.

[0138] Thus, in step S407, the processing circuit 106 according to the fourth embodiment determines the number of faulty channels allowed in the PET detector 101 or each of the multiple detector modules based on the number of faulty channels when the sensitivity of the PET detector 101 does not meet a predetermined standard in at least one arrangement pattern.

[0139] The number of allowable defective channels (48×(n-1)) determined by the process for determining the number of allowable defective channels is used in an acceptance test or the like when outsourcing the PET detector 101 and receiving the PET detector 101 from the outsourcer. For example, if the number of defective channels in one detector module notified by the outsourcer is n' and n' is equal to or less than (n-1), the user accepts the PET detector 101 manufactured by the outsourcer. On the other hand, if n' is greater than (n-1), the user does not accept the PET detector 101 manufactured by the outsourcer. Then, the user again outsources a PET detector 101 in which the number of defective channels in one detector module is equal to or less than (n-1).

[0140] Also, for example, if the number of defective channels of the PET detector 101 notified by the outsourcer is N and N is equal to or less than (48×(n−1)), the user accepts the PET detector 101 manufactured by the outsourcer. On the other hand, if N is greater than (48×(n−1)), the user does not accept the PET detector 101 manufactured by the outsourcer. Then, the user re-outsources a PET detector 101 with the number of defective channels equal to or less than (48×(n−1)).

[0141] The above has described the PET device 100 according to the fourth embodiment. The PET device 100 according to the fourth embodiment can determine the number of allowable defective channels in the PET detector 101 alone or in the detector module alone.

[0142] Furthermore, the PET device 100 according to the fourth embodiment performs the same processing as the processing performed by the PET device 100 according to the first embodiment, the second embodiment, the third embodiment, or the modified example of the first embodiment. Therefore, the PET device 100 according to the fourth embodiment achieves the same effects as the PET device 100 according to the first embodiment, the second embodiment, the third embodiment, or the modified example of the first embodiment.

[0143] In steps S106, S107, S202, S203, and S302, the display control function 106i may further display the sensitivity map generated in step S103 on the display 120. The sensitivity indicated by the sensitivity map generated in step S103 is closer to the actual state of the PET detector 101 than the sensitivity indicated by the sensitivity map 130a. Therefore, by displaying the sensitivity map generated in step S103 on the display 120, the user can grasp the sensitivity closer to the actual state of the PET detector 101. Furthermore, because the sensitivity map 130a is stored in advance in the memory 130, the user can grasp the distribution of the sensitivity indicated by the sensitivity map 130a at any desired timing. Therefore, by comparing the sensitivity map 130a with the sensitivity map generated in step S103, the user can grasp the locations where the sensitivity differs, and ultimately the range affected by the bad channel. Therefore, by displaying the sensitivity map generated in step S103 on the display 120, the influence of a poor channel on image quality can be evaluated quantitatively and objectively.

[0144] Furthermore, in steps S106, S107, S202, S203, and S302, the display control function 106i may further cause the sensitivity map 130a acquired in step S102 and the sensitivity map generated in step S103 to be displayed side by side so as to be comparable on the display 120. By displaying the sensitivity map 130a and the sensitivity map generated in step S103 so as to be comparable, the user can grasp the locations where the sensitivities differ, and ultimately the user can grasp the range affected by a poor channel.

[0145] Furthermore, in the above-described embodiment or modified example, the PET device 100 performs various processes using a two-dimensional sensitivity map, but the PET device 100 may perform similar processes using a three-dimensional sensitivity map.

[0146] The program executed by the processor is provided in advance in a read-only memory (ROM) or a storage circuit. The program may be provided in a format installable or executable by these devices, recorded on a non-transitory computer-readable storage medium such as a compact disk (CD)-ROM, a flexible disk (FD), a recordable CD-R, or a digital versatile disk (DVD). The program may also be stored on a computer connected to a network such as the Internet and provided or distributed by downloading it via the network. For example, the program may be composed of modules including the above-described processing functions. In actual hardware, a CPU reads and executes the program from a storage medium such as a ROM, whereby each module is loaded into a main memory device and generated on the main memory device.

[0147] According to at least one of the embodiments described above, the influence of a poor channel on image quality can be quantitatively and objectively evaluated.

[0148] Although several embodiments have been described, these embodiments are presented as examples and are not intended to limit the scope of the invention. These embodiments can be implemented in various other forms, and various omissions, substitutions, and modifications can be made without departing from the spirit of the invention. These embodiments and their modifications are included within the scope and spirit of the invention, as well as within the scope of the invention and its equivalents as defined in the claims. [Explanation of symbols]

[0149] 100 PET equipment 106e Acquisition Function 106f generation function

Claims

1. an acquisition unit that acquires information about a defective channel of the PET detector at a second timing that is later than the first timing corresponding to a first sensitivity map that is represented by a density of a line of response (LOR) for each pixel and that is a sensitivity map of the PET detector at the first timing stored in a storage unit; a generating unit that generates a second sensitivity map, which is a sensitivity map of the PET detector at the second timing, based on information about the failed channel, and is expressed by a density of LORs for each pixel; a calculation unit that calculates an index value relating to the degree of influence of the poor channel on image quality based on the first sensitivity map and the second sensitivity map; a determination unit that determines whether or not the sensitivity of the PET detector satisfies a predetermined standard based on the index value; a display control unit that causes a determination result by the determination unit to be displayed on a display unit; Equipped with when the determination unit determines that the sensitivity of the PET detector does not satisfy a predetermined standard, the generation unit changes a position of the defective channel in the PET detector, and generates a third sensitivity map that is a sensitivity map of the PET detector based on the position of the defective channel after the change; the calculation unit calculates an index value relating to the degree of influence of the defective channel whose position has been changed on image quality based on the first sensitivity map and the third sensitivity map; the determination unit determines whether or not the sensitivity of the PET detector satisfies a predetermined standard based on the index value calculated based on the first sensitivity map and the third sensitivity map; The display control unit causes the display unit to display information indicating the position of the defective channel in the PET detector after the change when the judgment unit judges that the sensitivity of the PET detector satisfies a predetermined standard.

2. The PET device according to claim 1 , further comprising a display control unit that causes the second sensitivity map to be displayed on a display unit.

3. The PET apparatus according to claim 1 , wherein the calculator calculates the index value based on a ratio or a difference between the first sensitivity map and the second sensitivity map.

4. 2. The PET device according to claim 1, further comprising: a display control unit that causes a display unit to display an image based on image data obtained by dividing the second sensitivity map by the first sensitivity map, or an image based on image data obtained by subtracting the second sensitivity map from the first sensitivity map.

5. 4. The PET device according to claim 1, further comprising a reconstruction processing unit that reconstructs PET image data based on the second sensitivity map.

6. The PET device according to claim 1 , further comprising a reconstruction processing unit that reconstructs PET image data based on the third sensitivity map when the determination unit determines that the sensitivity of the PET detector satisfies a predetermined standard.

7. The PET device according to claim 5 , wherein the reconstruction processing unit reconstructs the PET image data based on an inverse of the second sensitivity map.

8. an acquisition unit that acquires information about a failed channel of the PET detector at a first timing; a generator that generates a first sensitivity map, which is a sensitivity map of the PET detector at the first timing, based on information about the defective channel, and generates a second sensitivity map having improved sensitivity compared to the first sensitivity map, based on the information about a change in the position of the defective channel; A PET device comprising:

9. The PET device according to claim 8 , wherein the generator generates, as the second sensitivity map, a sensitivity map of the PET detector after the defective channel is replaced with a normal channel that is not defective.

10. 9. The PET device according to claim 8, wherein the generator generates, as the second sensitivity map, a sensitivity map of the PET detector after the position of the defective channel has been changed to a position different from the position of the defective channel at the first timing.

11. Acquire information about a defective channel of the PET detector at a second timing that is later than the first timing corresponding to a first sensitivity map represented by the density of LORs for each pixel, the first sensitivity map being a sensitivity map of the PET detector at the first timing stored in a storage unit; generating a second sensitivity map, which is a sensitivity map of the PET detector at the second timing based on information about the failed channel, and is represented by a density of LORs for each pixel; calculating an index value relating to the degree of influence of the defective channel on image quality based on the first sensitivity map and the second sensitivity map; determining whether the sensitivity of the PET detector satisfies a predetermined standard based on the index value; The judgment result is displayed on the display unit. if it is determined that the sensitivity of the PET detector does not satisfy a predetermined standard, changing the position of the defective channel in the PET detector, and generating a third sensitivity map that is a sensitivity map of the PET detector based on the position of the defective channel after the change; calculating an index value relating to the degree of influence of the defective channel whose position has been changed on image quality based on the first sensitivity map and the third sensitivity map; determining whether or not the sensitivity of the PET detector satisfies a predetermined standard based on the index value calculated based on the first sensitivity map and the third sensitivity map; a display unit displaying information indicating the position of the defective channel in the PET detector after the change when it is determined that the sensitivity of the PET detector satisfies a predetermined standard.

12. Obtaining information about a failed channel of the PET detector at a first timing; generating a first sensitivity map, which is a sensitivity map of the PET detector at the first timing, based on information about the failed channel; A method for generating a second sensitivity map having improved sensitivity compared to the first sensitivity map based on the information that the positions of the defective channels have been changed.

13. The computer a process of changing the positions of defective channels in a PET detector model or models of a plurality of detector modules included in the PET detector model and generating sensitivity maps of the PET detector model based on the changed positions of the defective channels a plurality of times, thereby generating a plurality of sensitivity maps corresponding to a plurality of defective channel arrangement patterns in which the positions of the defective channels are different from each other, while increasing the number of the defective channels, until the sensitivity of the PET detector model no longer satisfies a predetermined standard in at least one of the plurality of defective channel arrangement patterns; A method for determining the number of defective channels allowed in the PET detector or each of a plurality of detector modules of the PET detector based on the number of defective channels when the sensitivity of the PET detector model in the at least one arrangement pattern does not satisfy a predetermined criterion.

14. On the computer, Acquire information about a defective channel of the PET detector at a second timing that is later than the first timing corresponding to a first sensitivity map represented by the density of LORs for each pixel, the first sensitivity map being a sensitivity map of the PET detector at the first timing stored in a storage unit; generating a second sensitivity map, which is a sensitivity map of the PET detector at the second timing based on information about the failed channel, and is represented by a density of LORs for each pixel; calculating an index value relating to the degree of influence of the defective channel on image quality based on the first sensitivity map and the second sensitivity map; determining whether the sensitivity of the PET detector satisfies a predetermined standard based on the index value; The judgment result is displayed on the display unit. if it is determined that the sensitivity of the PET detector does not satisfy a predetermined standard, changing the position of the defective channel in the PET detector, and generating a third sensitivity map that is a sensitivity map of the PET detector based on the position of the defective channel after the change; calculating an index value relating to the degree of influence of the defective channel whose position has been changed on image quality based on the first sensitivity map and the third sensitivity map; determining whether or not the sensitivity of the PET detector satisfies a predetermined standard based on the index value calculated based on the first sensitivity map and the third sensitivity map; A program for executing a process of displaying, on the display unit, information indicating the position of the defective channel in the PET detector after the change when it is determined that the sensitivity of the PET detector satisfies a predetermined standard.

15. On the computer, Obtaining information about a failed channel of the PET detector at a first timing; generating a first sensitivity map, which is a sensitivity map of the PET detector at the first timing, based on information about the failed channel; a program for executing a process of generating a second sensitivity map having improved sensitivity compared to the first sensitivity map, based on the information on the change in the position of the defective channel;

16. On the computer, a process of changing the positions of defective channels in a PET detector model or models of a plurality of detector modules included in the PET detector model and generating sensitivity maps of the PET detector model based on the changed positions of the defective channels a plurality of times, thereby generating a plurality of sensitivity maps corresponding to a plurality of defective channel arrangement patterns in which the positions of the defective channels are different from each other, while increasing the number of the defective channels, until the sensitivity of the PET detector model no longer satisfies a predetermined standard in at least one of the plurality of defective channel arrangement patterns; A program for executing a process of determining the number of defective channels allowed in the PET detector or each of the multiple detector modules of the PET detector based on the number of defective channels when the sensitivity of the PET detector model in the at least one arrangement pattern does not satisfy a predetermined standard.

Images