Pet apparatus, method, and recording medium

A Positron Emission Tomography (PET) apparatus according to an embodiment includes processing circuitry. The processing circuitry is configured to obtain information about a defective channel of a PET detector at a second point in time later than a first point in time corresponding to a first sensitivity map that is a sensitivity map of the PET detector corresponding to the first point in time and being stored in a storage unit. The processing circuitry is configured to generate a second sensitivity map that is a sensitivity map of the PET detector corresponding to the second point in time, on the basis of the information about the defective channel.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is based upon and claims the benefit of priority from Japanese Patent Application No. 2021-034629, filed on Mar. 4, 2021; the entire contents of which are incorporated herein by reference.

FIELD

Embodiments described herein relate generally to a Positron Emission Tomography (PET) apparatus, a method, and a recording medium.

BACKGROUND

A Positron Emission Tomography (PET) detector employed in a PET apparatus includes a scintillator array and an optical detector array. The scintillator array is structured with a plurality of scintillators. The optical detector array is structured with a plurality of optical detectors (optical detecting elements). In this regard, for example, the optical detector array may be structured, in some situations, with a plurality of Photomultiplier Tubes (PMTs) serving as the plurality of optical detectors. A PET apparatus including such an optical detector array may adopt a multiplex method, for example. In this situation, in the entire PET detector, the number of channels (a total number of signal output channels) is approximately hundreds, because the quantity of the PMTs is smaller than the quantity of a large number of scintillators structuring the scintillator array.

Because fluorescent light from the large number of scintillators becomes incident to the smaller number of PMTs and is converted into signals, if a signal output from one of the PMTs has a loss, a non-negligible impact is made on a large part of an image. Accordingly, during a health check to diagnose (to inspect) a failure of PET detectors, taking an actual image of a phantom makes it qualitatively clearer whether a failure is present or absent than when a user visually inspects images. For this reason, when a PET detector fails, there is no other option besides repairing or replacing the PET detector. In this sense, there is no need to make special judgment.

In another example, the optical detector array may be structured, in some situations, with a plurality of Silicon Photomultipliers (SiPMs) serving as the plurality of optical detectors. In this example, each of the scintillators may be optically coupled with a different one of the SiPMs (i.e., the scintillators and the SiPMs are in one-to-one optical coupling) or each group of a small number of scintillators may be optically coupled with a different one of the SiPMs. Accordingly, the number of independent channels (the number of signal output channels) in the entire PET detector is as many as tens of thousands to hundreds of thousands. If the channels include one or more defective channels, changes in or impacts on image quality of images are not so significant as to be immediately noticed in visual inspection of users unlike in the example of the PET detector in which the plurality of PMTs are used in the optical detector array and may vary depending on the quantity of the defective channels and the relative positions of the defective channels within the PET detector. In this situation, the term “defective channels” denotes channels corresponding to failed optical detectors, for example. Changes in uniformity of an image do not necessarily occur evenly in a large area and may occur in units of pixels or in a component. For this reason, the visual inspection method to inspect an image that is obtained by imaging a phantom and contains statistical noise lacks in quantitativeness and objectiveness. In other words, when each of the scintillators is optically coupled with a different one of the SiPMs or when each group of a small number of scintillator is optically coupled with a different one of the SiPMs, because the impact on the image quality (or the image) imposed by individual channel defects is relatively minor, criteria used for determining countermeasures are not definitive.

Accordingly, with regard to PET detectors having a large number of independent channels, there is a demand for quantitatively and objectively evaluating impacts on image quality imposed by defective channels. It is important that users correctly determine countermeasures on the basis of such an evaluation, from the viewpoint of guaranteeing safety and capabilities in the clinical field and from the viewpoint of properly allocating costs in the fields of manufacturing and services.

DETAILED DESCRIPTION

One of the problems to be solved by the embodiments disclosed in the present specification and drawings is to quantitatively and objectively evaluate impacts on image quality imposed by defective channels. However, possible problems to be solved by the embodiments disclosed in the present specification and drawings are not limited to the problem stated above. It is also possible to consider any of the problems corresponding to advantageous effects achieved by the configurations described in the embodiments presented below, as other problems.

A PET apparatus according to an embodiment includes processing circuitry. The processing circuitry is configured to obtain information about a defective channel of a PET detector at a second point in time later than a first point in time corresponding to a first sensitivity map that is a sensitivity map of the PET detector corresponding to the first point in time and being stored in a storage unit. The processing circuitry is configured to generate a second sensitivity map that is a sensitivity map of the PET detector corresponding to the second point in time, on the basis of the information about the defective channel.

Exemplary embodiments of a PET apparatus, a method, and a computer program (hereinafter, “program”) will be explained in detail below, with reference to the accompanying drawings. The PET apparatus, the method, and the program of the present disclosure are not limited to the embodiments described below. Further, the description of each of the embodiments is, in principle, similarly applicable to any other embodiment or modification example. Also, the description of each of the modification examples is, in principle, similarly applicable to any of the embodiments and the other modification examples.

First Embodiment

FIG.1is a diagram illustrating an exemplary configuration of a PET apparatus100according to a first embodiment. As illustrated inFIG.1, the PET apparatus100according to the first embodiment includes a gantry device10and a console device20. For instance, the PET apparatus100is an example of a nuclear medicine diagnosis apparatus.

The gantry device10includes a PET detector101, a count information generating circuit102, a tabletop103, a table104, and a table driving unit105.

The PET detector101includes a plurality of detector modules. The plurality of detector modules are configured to detect radiation by detecting scintillation photons (fluorescent light) that are re-released when a substance in an excited state transitions back into a ground state as a result of an interaction between annihilation gamma rays released from positrons inside an examined subject (hereinafter, “patient”) P and a light emitting bodies (the scintillators). The plurality of detector modules are configured to detect radiation energy information of the annihilation gamma rays released from the positrons inside the patient P. The plurality of detector modules are arranged so as to surround the patient P in a ring formation.

Each of the detector modules includes, for example, a scintillator array, an optical detector array, and a light guide.

The scintillator array includes a plurality of scintillators arranged two-dimensionally. The scintillators are configured to convert annihilation gamma rays that have become incident thereto after being released from the positrons inside the patient P, into scintillation photons (or optical photons) and to output the scintillation photons. The scintillators are formed with scintillator crystals that are suitable for a Time Of Flight (TOF) measuring process or an energy measuring process, such as those of Lanthanum Bromide (LaBr3), Lutetium Yttrium Oxyorthosilicate (LYSO), Lutetium Oxyorthosilicate (LSO), or Lutetium Gadolinium Oxyorthosilicate (LGSO), for example.

The light guide is formed by using a plastic material or the like having excellent light transmittance and is configured to transfer the scintillation photons output from the scintillators to the optical detectors. More specifically, the light guide is configured to transfer the scintillation photons to Silicon Photomultipliers (SiPMs) explained later.

The optical detector array includes a plurality of optical detectors (optical detecting elements) arranged two-dimensionally. For example, the SiPMs may be used as the optical detectors. In the present embodiment, each of the scintillators is optically coupled with a different one of the SiPM. In other words, the scintillators and the SiPMs are in one-to-one optical coupling. In this situation, a set made up of one scintillator and one SiPM optically coupled with the one scintillator corresponds to one channel. Alternatively, a small number of scintillators may be optically coupled with a SiPM. In that situation, a set made up of the small number of scintillators and one SiPM optically coupled with the small number of scintillators corresponds to one channel. In either situation, the number of channels (the number of signal output channels) in the entirety of the PET detector101is as many as tens of thousands to hundreds of thousands. The SiPMs are configured to detect the scintillation photons emitted by the scintillators, to generate an electrical signal corresponding to the scintillation photons, and to output the generated electrical signal as an output signal.

The PET detector101is provided with the count information generating circuit102for each of the detector modules.

The count information generating circuit102is configured to generate digital data as count information, by converting the output signal from the PET detector101into the digital data. The count information includes detection positions of the annihilation gamma rays, energy values, and detection times. For example, the count information generating circuit102identifies a plurality of optical detectors that converted scintillation photons into electrical signals at mutually the same time. Further, the count information generating circuit102identifies scintillator numbers (P) indicating the positions of the scintillators (scintillator positions) to which the annihilation gamma rays became incident. The count information generating circuit102may use any of various methods, as the method for identifying the scintillator positions in which the annihilation gamma rays became incident. For example, because each one of the scintillators corresponds to a different one of the SiPMs, the count information generating circuit102may identify the position of a scintillator corresponding to a SiPM from which an output was obtained, as the scintillator position in which the annihilation gamma rays became incident.

Further, the count information generating circuit102is configured to identify energy values (E) of the annihilation gamma rays that became incident to the PET detector101, through an integral calculation on intensities of the electrical signals output from the optical detectors. Further, the count information generating circuit102is configured to identify detection times (T) at which the scintillation photons from the annihilation gamma rays were detected by the PET detector101. The detection times (T) may be absolute times or elapsed time periods since the start of the imaging process. As explained herein, the count information generating circuit102is configured to generate the count information including the scintillator numbers (P), the energy values (E), and the detection times (T). Further, the count information generating circuit102is configured to store the generated count information into a memory130of the console device20.

For example, the count information generating circuit102is realized by using a processor.

The tabletop103is a bed on which the patient P is placed and is arranged over the table104. The table driving unit105is configured to move the tabletop103under control of a table controlling function106dof processing circuitry106. For example, the table driving unit105is configured to move the patient P to the inside of an imaging opening of the gantry device10, by moving the tabletop103.

Upon receipt of an operation performed by a user on the PET apparatus100, the console device20is configured to control imaging of PET image data and to reconstruct (generate) the PET image data by using the count information acquired by the gantry device10. As illustrated inFIG.1, the console device20includes the processing circuitry106, an input interface110, a display device120, and the memory130. In this situation, the processing circuitry106, the input interface110, the display device120, and the memory130are connected together via a bus.

The processing circuitry106includes a coincidence information generating function106a, a reconstruction processing function106b, a system controlling function106c, the table controlling function106d, an obtaining function106e, a generating function106f, a calculating function106g, a judging function106h, and a display controlling function106i. The functions such as the coincidence information generating function106a, the reconstruction processing function106b, the system controlling function106c, the table controlling function106d, the obtaining function106e, the generating function106f, the calculating function106g, the judging function106h, and the display controlling function106iare stored in the memory130in the form of computer-executable programs. The processing circuitry106is a processor configured to realize the functions corresponding to the programs, by reading the programs from the memory130and executing the read programs. In other words, the processing circuitry106that has read the programs has the functions illustrated within the processing circuitry106inFIG.1. With reference toFIG.1, the example was explained in which the single piece of processing circuitry (i.e., the processing circuitry106) realizes the coincidence information generating function106a, the reconstruction processing function106b, the system controlling function106c, the table controlling function106d, the obtaining function106e, the generating function106f, the calculating function106g, the judging function106h, and the display controlling function106i. However, another arrangement is also acceptable in which the processing circuitry106is structured with a plurality of independent processors, so that the functions are realized as a result of the processors executing the programs.

In this situation, for instance, the reconstruction processing function106bis an example of a reconstruction processing unit. For instance, the obtaining function106eis an example of an obtaining unit. For instance, the generating function106fis an example of a generating unit. For instance, the calculating function106gis an example of a calculating unit. For instance, the judging function106his an example of a judging unit. For instance, the display controlling function106iis an example of a display controlling unit.

The term “processor” used in the above explanation denotes, for example, a Central Processing Unit (CPU), a Graphical Processing Unit (GPU), or a circuit such as an Application Specific Integrated Circuit (ASIC) or a programmable logic device (e.g., a Simple Programmable Logic Device [SPLD], a Complex Programmable Logic Device [CPLD], or a Field Programmable Gate Array [FPGA]). The one or more processors are configured to realize the functions by reading the programs saved in the memory130and executing the read programs.

The coincidence information generating function106ais configured to obtain the count information stored in the memory130and to generate pieces of coincidence information on the basis of the obtained count information. Further, the coincidence information generating function106ais configured to arrange the generated pieces of coincidence information in a substantially chronological order on the basis of the detection times (T) and to store the arranged information into the memory130. As a result, the memory130has stored therein a time-series list of the coincidence information.

The reconstruction processing function106bis configured to reconstruct the PET image data. For example, the reconstruction processing function106bis configured to obtain the time-series list of the coincidence information stored in the memory130and to reconstruct the PET image data by using the obtained time-series list of the coincidence information. Further, the reconstruction processing function106bis configured to store the reconstructed PET image data into the memory130.

The system controlling function106cis configured to control the entirety of the PET apparatus100, by controlling the gantry device10and the console device20. For example, the system controlling function106cis configured to control imaging processes performed by the PET apparatus100.

The table controlling function106dis configured to control the moving of the tabletop103, by controlling the table driving unit105. The obtaining function106e, the generating function106f, the calculating function106g, and the judging function106hwill be explained later.

The display controlling function106iis configured to cause the display device120to display various types of images and various types of information. For example, the display controlling function106iis configured to cause the display device120to display a PET image based on the PET image data. Also, the display controlling function106iis configured to cause the display device120to display a Graphical User interface (GUI) used for receiving various types of instructions and various types of settings from the user (an operator) of the PET apparatus100.

The input interface110is configured to receive the inputs of the various types of instructions and the various types of settings from the user and to output the received various types of instructions and various types of settings to the processing circuitry106. For example, the input interface110is configured to convert the input operations received from the user into electrical signals and to transmit the electrical signals to the processing circuitry106. For example, the input interface110is realized by using a trackball, a switch button, a mouse, a keyboard, a touchpad on which input operations can be performed by touching an operation surface thereof, a touch screen in which a display screen and a touchpad are integrally formed, a contactless input circuit using an optical sensor, an audio input circuit, and/or the like. In the present disclosure, the input interface110does not necessarily have to include physical operation component parts such as a mouse, a keyboard, and/or the like. For instance, possible examples of the input interface110include electrical signal processing circuitry configured to receive an electrical signal corresponding to an input operation from an external input device provided separately from the PET apparatus100and to transmit the electrical signal to the processing circuitry106. For instance, the input interface110is an example of a receiving unit.

The display device120is connected to the processing circuitry106and is configured to display the various types of information and the various types of information. For example, the display device120is configured to convert information and data transmitted thereto from the processing circuitry106into display-purpose electrical signals and to output the electrical signals. In a specific example, under the control of the display controlling function106i, the display device120is configured to display the PET image based on the PET image data and to display the GUT used for receiving the various types of instructions and the various types of settings from the user. For example, the display device120is realized by using a liquid crystal monitor, a Cathode Ray Tube (CRT) monitor, a touch panel, or the like. The display device120is an example of a display unit.

The memory130is configured to store therein various types of data used in the PET apparatus100. The memory130is an example of the storage unit. For example, the memory130is realized by using a semiconductor memory element such as a Random Access Memory (RAM) or a flash memory, or a hard disk, an optical disk, or the like. The memory130is configured to store therein the count information which is the information in which the scintillator numbers (P), the energy values (E), and the detection times (T) are kept in correspondence with one another, the time-series list of the coincidence information in which coincidence numbers serving as serial numbers of the pieces of coincidence information are kept in correspondence with sets of count information, the reconstructed PET image data, and the like.

Further, the memory130is configured to store therein a sensitivity map130a.FIG.2is a drawing illustrating an example of the sensitivity map130aaccording to the first embodiment. The sensitivity map130aindicates sensitivity of the PET detector101with respect to each of the pixels in a region (a reconstruction region) subject to a reconstructing process. The pixel size may be 4 mm×4 mm, for example, but the pixel size is not limited to this example. The sensitivity map130ais obtained when the PET detector101has no failure, so that all the channels of the PET detector101are normal channels without being defective channels. In other words, the sensitivity map130ais a sensitivity map of the PET detector101corresponding to a point in time when all the channels of the PET detector101are normal channels. For example, the sensitivity map130ais an example of the first sensitivity map. Further, the point in time when all the channels of the PET detector101are normal channels is an example of the first point in time.

Next, an example of a method for generating the sensitivity map130awill be explained. The sensitivity map130amay be generated by the PET apparatus100or may be generated by an apparatus other than PET apparatus100. In the following explanations, an example will be explained in which the generating function106fof the PET apparatus100is configured to generate the sensitivity map130a.

FIG.3is a drawing for explaining an example of the method for generating the sensitivity map130aaccording to the first embodiment. As illustrated inFIG.3, the generating function106fis configured to generate a model30of the PET detector101in which models of the plurality of detector modules are arranged in a ring formation.

In the example inFIG.3, the generating function106fis configured to generate the model30of the PET detector101in which the models of 48 detector modules are arranged in the ring formation. The model of each of the detector modules includes at least one set31made up of a model of one scintillator and a model of one SiPM. For example, it is preferable when the quantity and the positional arrangements of the sets31included in the model30are the same as the quantity and the positional arrangements of sets that are each made up of one scintillator and one SiPM and are included in the PET detector101. For example, it is preferable when the quantity of the sets31included in the model30is a number between tens of thousands and hundreds of thousands that is equal to the quantity of the sets that are each made of one scintillator and one SiPM and are included in the PET detector101.

However, for the sake of convenience in the explanation, the following will explain an example in which the model of each of the detector modules includes only one set31. In other words, in the model30of the PET detector101illustrated inFIG.3,48sets31are arranged in the ring formation. In this situation, the number of channels in the model30of the PET detector101is “48”. In the present example, all the channels are normal channels. Further, the48channels are identified with identifiers “1” to “48”, respectively.

Further, the generating function106fis configured to draw a Line Of Response (LOR)32connecting two channels among the48channels, with respect to every combination each made up of two channels. In the example inFIG.3, the generating function106fis configured to draw 1,128 (=48×47/2) LORs32. A pattern of the LORs32drawn in this manner may be referred to as an LOR pattern.

As illustrated inFIG.3, the density of the LORs32is not uniform. There are regions where the density of the LORs32is relatively low and where the density is relatively high. Annihilation gamma rays emitted from a high-density region are detected more than annihilation gamma rays emitted from a low-density region. In other words, a high-density region has a higher possibility of the occurrence of annihilation than a low-density region.

Further, the generating function106fis configured to generate the sensitivity map130a, by converting the LOR pattern in a reconstruction region33of the model30, into image data expressing the density of the LORs for each of the pixels. In other words, on the basis of the LOR pattern in the reconstruction region33of the model30, the generating function106fcalculates the density of the LORs for each of the pixels and generates the image data in which each of the pixels indicates the density of the LORs, as the sensitivity map130a. In other words, the sensitivity map130ais the data in which, with respect to each of the pixels, the density of the LORs is registered as sensitivity.

The sensitivity map130agenerated in this manner is stored in the memory130in advance and is used, during the PET image data reconstructing process performed by the reconstruction processing function106b, for normalizing non-uniformity of sensitivity caused by the geometric arrangement of the channels (the sets each made up of the model of one scintillator and the model of one SiPM). For example, by using the sensitivity map130a, the reconstruction processing function106bis configured to reconstruct the PET image data according to Expression (1) presented below:

For example, Expression (1) is a recurrence relation to calculate the pixel value λk+1of a pixel j in the PET image data to be eventually obtained as a result of the reconstruction, by using the pixel value λkof the pixel j and the reciprocal (1/Sj) of the pixel value Sjof the pixel in the sensitivity map S. Expression (1) is an expression indicating an algorithm implementing a Maximum Likelihood—Expectation Maximization (ML-EM) scheme in a statistical successive approximation method. In Expression (1), the detection probability Cimmay be called a system matrix and denotes a probability that a photon occurring from an m-th pixel is counted by an i-th LOR. The element yidenotes an actual counted value, i.e., the number of events counted by the i-th LOR during the actual acquisition. According to the ML-EM scheme, a certain radiation distribution is postulated as an initial value. In other words, λ (k=0) is assumed. According to the ML-EM scheme, with respect to the postulated distribution, the number of events expected to be counted by the i-th LOR is calculated by calculating the sum of Cimand contribution of the pixel m. According to the ML-EM scheme, a ratio between the calculated number of events and the number of actual counted events (the number of actual counted events indicated by the numerator in Expression (1)) is calculated with respect to each of the LORs. If the ratio is 1 with respect to each of all the LORs, it means that the postulated pixel values reproduce the actual counted events. In that situation, the postulated pixel values reproduce the actual radiation distribution of the imaged subject. As for the operation of Expression (1), the pixel value Sjof the sensitivity map S is defined as ΣiCijif there is no channel loss. As a result, the factor on the right starting with Σ is equal to Sj, while a k-th pixel value is equal to a (k+1)-th pixel value, which means that the calculation is converged.

As explained above, the reconstruction processing function106bis configured to reconstruct the PET image data, by using the reciprocal of the sensitivity of each of the pixels indicated in the sensitivity map130a.

The exemplary configuration of the PET apparatus100according to the present embodiment has thus been explained. For example, the PET apparatus100is installed in a medical facility such as a hospital or a clinic and is used in various types of image diagnosis processes using the PET image data generated by the PET apparatus100, for the patient P who is a patient hospitalized or visiting the medical facility. In the present embodiment, the PET apparatus100is configured to perform various types of processes described below, so as to be able to quantitatively and objectively evaluate impacts on image quality (or images) imposed by defective channels.

Next, an example of an evaluating process performed by the PET apparatus100will be explained.FIG.4is a flowchart illustrating an example of a flow in the evaluating process performed by the PET apparatus100according to the first embodiment. For example, the PET apparatus100is configured to perform the evaluating process, when the user inputs, to the processing circuitry106, an instruction (an execution instruction) to perform the evaluating process presented inFIG.4, via the input interface110. In this situation, for example, the user inputs the execution instruction to the processing circuitry106via the input interface110, so as to understand whether or not the impacts on image quality imposed by defective channels satisfy a predetermined quality criteria, after the PET apparatus100is installed in an imaging room of the medical facility, after regular maintenance of the PET apparatus100, after the PET apparatus100is repaired, after an emergency countermeasure is taken when the PET apparatus100had a failure, or the like. In other examples, besides the abovementioned occasions, the user may input the execution instruction to the processing circuitry106via the input interface110, so as to understand whether or not the impacts on image quality imposed by defective channels satisfy the predetermined quality criteria, when performing an acceptance test on the PET detector101by itself, when testing the PET detector101by itself before shipment, when inspecting the PET apparatus100before a work day starts, when determining a repair plan for the PET apparatus100, or the like.

As illustrated inFIG.4, at step S101, the obtaining function106eobtains failure information (defect information) serving as information about defective channels. At the time of performing the evaluating process, there may be one or more defective channels among the channels of the PET detector101in some situations, and all the channels may be normal channels in other situations. In the following sections, an example will be explained in which there are one or more defective channels among the channels of the PET detector101at the time of performing the evaluating process. The failure information about the defective channels include, for example, the positions in which the defective channels are arranged and the quantity thereof. For example, the obtaining function106eperforms a health check to inspect whether or not the PET detector101is working properly and obtains the failure information about the defective channels. After that, the obtaining function106estores the obtained failure information into the memory130. In this manner, at step S101, the obtaining function106eobtains the failure information about the defective channels of the PET detector101at a point in time when the defective channels are present, which is later than the point in time which corresponds to the sensitivity map130astored in the memory130and at which all the channels of the PET detector101were normal channels. The point in time when the defective channels are present is an example of the second point in time.

Subsequently, at step S102, the obtaining function106eobtains, from the memory130, the sensitivity map130athat is a sensitivity map of the PET detector101corresponding to the time when all the channels of the PET detector101were normal.

After that, at step S103, on the basis of the failure information, the generating function106fgenerates a sensitivity map of the PET detector101corresponding to the point in time when the defective channels are present. In other words, at step S103, the sensitivity map reflecting the defective channels is generated, while the defective channels that are actually present are taken into consideration. For instance, the sensitivity map generated at step S103is an example of the second sensitivity map. The following sections will explain a specific example of the process at step S103. To begin with, for example, the generating function106fobtains the failure information from the memory130.

After that, by using the same method as the method for generating the sensitivity map130aexplained with reference toFIG.3, the generating function106fgenerates the sensitivity map taking the defective channels indicated in the failure information into consideration. For example, the generating function106fidentifies the plurality of normal channels excluding the defective channels identified in the failure information, from among the48channels in the model30of the PET detector101.

After that, the generating function106fdraws a LOR32connecting two channels among the identified plurality of normal channels, with respect to every combination each made up of two channels. Further, the generating function106fgenerates the sensitivity map, by converting the LOR pattern in the reconstruction region33of the model30, into image data expressing the density of the LORs for each of the pixels. In other words, on the basis of the LOR pattern in the reconstruction region33of the model30, the generating function106fcalculates the density of the LORs for each of the pixels, as sensitivity of the PET detector101, and generates the image data in which each of the pixels indicates the sensitivity, as the sensitivity map. According to the method described herein, the generating function106fgenerates, at step S103, the sensitivity map registering the sensitivity with respect to each of the pixels, while taking the defective channels into consideration.

In this situation, the positional arrangements and the quantity of the plurality of pixels structuring the sensitivity map generated at step S103are the same as the positional arrangements and the quantity of the plurality of pixels structuring the sensitivity map130a. Accordingly, the plurality of pixels structuring the sensitivity map generated at step S103correspond, respectively, to the plurality of pixels structuring the sensitivity map130a.

In this situation, the sensitivity map130aused at the time of reconstructing the PET image data is a sensitivity map that was generated in advance and saved in the memory130and that corresponds to the time when all the channels were normal channels. For this reason, differences in sensitivity indicated in the sensitivity map generated at step S103in comparison to the sensitivity indicated in the sensitivity map130aare exhibited in the sensitivity map generated at step S103, as the differences in sensitivity caused by the defective channels.

The following sections will explain the situation where the sensitivity map generated at step S103may vary when the positional arrangements of the defective channels are different, even when the quantity of the defective channels is the same.FIGS.5and6are drawings illustrating an example of the sensitivity map generated at step S103according to the first embodiment.

A sensitivity map35illustrated inFIG.5is a sensitivity map generated at step S103when the channel identified with an identifier “1” and the channel identified with an identifier “5” are defective channels. In contrast, another sensitivity map35illustrated inFIG.6is a sensitivity map generated at step S103when the channel identified with the identifier “1” and the channel identified with an identifier “10” are defective channels. As understood when the sensitivity map35inFIG.5is compared with the sensitivity map35inFIG.6, it is observed that the two sensitivity maps35are different from each other, when the positional arrangements of the defective channels are different even through the quantity of the defective channels is the same. The reason why the two sensitivity maps35are different from each other in this manner is, for example, that the relative positions of the defective channels are different between the sensitivity map35inFIG.5and the sensitivity map35inFIG.6.

Subsequently, at step S104, the calculating function106gcalculates index values related to a degree of impacts on the image quality imposed by the defective channels, on the basis of the sensitivity map130aand the sensitivity map generated at step S103. The following sections will explain a specific example of the process at step S104. For example, the calculating function106ggenerates image data by dividing the sensitivity map generated at step S103by the sensitivity map130a. The positional arrangements and the quantity of the plurality of pixels structuring the image data generated by the calculating function106gare the same as the positional arrangements and the quantity of the plurality of pixels structuring the sensitivity map130a. Also, each of the pixel values of the plurality of pixels structuring the image data generated by the calculating function106gis a value obtained by dividing the sensitivity indicated by a corresponding one of the plurality of pixels structuring the sensitivity map generated at step S103, by the sensitivity indicated by a corresponding one of the plurality of pixels structuring the sensitivity map130a. In other words, the calculating function106gcalculates the index values each corresponding to a different one of the pixels, on the basis of ratios between the sensitivity map130aand the sensitivity map generated at step S103.

At step S104, by generating the image data described above, the calculating function106gcalculates the pixel values of the plurality of pixels structing the image data as the index values related to the degree of impacts on the image quality imposed by the defective channels. In other words, the calculating function106gevaluates, with respect to each of the pixels, the difference in sensitivity caused by the defective channels and calculates the evaluation result of each of the pixels as the index value.

In this situation, each of the index values falls in the range from 0 to 1 inclusive. As an index value approaches 1, the degree of impacts on the image quality imposed by the defective channels becomes smaller. Conversely, as an index value approaches 0, the degree of impacts on the image quality imposed by the defective channels becomes larger.

Alternatively, the calculating function106gmay calculate an index value for each of the pixels, on the basis of the differences between the sensitivity map130aand the sensitivity map generated at step S103. In an example, the calculating function106gmay generate image data by subtracting the sensitivity map generated at step S103from the sensitivity map130a. In that situation, the pixel value of each of the plurality of pixels structuring the image data generated by the calculating function106gis a value obtained by subtracting the sensitivity indicated by a corresponding one of the plurality of pixels structuring the sensitivity map generated at step S103, from the sensitivity indicated by a corresponding one of the plurality of pixels structuring the sensitivity map130a. In this situation, each of the index values is a value equal to or larger than 0. As an index value approaches 0, the degree of impacts on the image quality imposed by the defective channels becomes smaller. Conversely, as an index value increases, the degree of impacts on the image quality imposed by the defective channels becomes larger.

Subsequently, at step S105, the judging function106hjudges whether or not the sensitivity of the PET detector101satisfies a predetermined criterion, on the basis of the index values. The following sections will explain a specific example of the process at step S105. For example, an example will be explained in which the PET apparatus100is required to have specifications where it is possible to quantitatively detect an aggregation of pixels with a pixel size of 4 mm×4 mm, while it is also required that a Standardized Uptake Value (SUV) evaluation can properly be performed.

In this situation, for example, a necessary condition of the PET detector101is set as that all the index values calculated in correspondence with the pixels at step S104are each equal to or larger than a threshold value α. The threshold value α may be “0.9”, for example, but the threshold value α is not limited to this example. In this situation, at step S105, the judging function106hjudges whether or not the sensitivity of the PET detector101satisfies the predetermined criterion, by judging whether or not all the index values calculated in correspondence with the pixels are each larger than the threshold value α. For example, when determining that all the index values calculated in correspondence with the pixels are each larger than the threshold value α, the Judging function106hdetermines that the sensitivity of the PET detector101satisfies the predetermined criterion. On the contrary, when determining that at least one among all the index values calculated in correspondence with the pixels is equal to or smaller than the threshold value α, the judging function106hdetermines that the sensitivity of the PET detector101does not satisfy the predetermined criterion.

When having determined that the sensitivity of the PET detector101satisfies the predetermined criterion (step S105: Yes), the judging function106hproceeds to step S106. On the contrary, when having determined that the sensitivity of the PET detector101does not satisfy the predetermined criterion (step S105: No), the judging function106hproceeds to step S107.

As explained above, at step S106and step S107, the display controlling function106icauses the display device120to display the judgment result obtained by the judging function106h.

The evaluating process presented inFIG.4is closed only with the geometric calculation using the result of the health check to inspect the defective channels of the PET detector101. Further, the evaluating process is closed only with the calculation at the level of the PET detector101and is not dependent on whether an acquisition count is large or small, the magnitude of statistical noise varied by the acquisition count being large/small, or experience of the user who interprets the image. Accordingly, the PET apparatus100according to the first embodiment is able to quantitatively and objectively evaluate the impacts on the image quality and diagnosing capabilities imposed by the defective channels in the PET detector101having a large number of channels.

Further, the PET apparatus100is able to have the user acknowledge the message indicating that the sensitivity of the PET detector101satisfies the predetermined criterion or the message indicating that the sensitivity of the PET detector101does not satisfy the predetermined criterion. Furthermore, the PET apparatus100is able to make it possible for the user who acknowledged the message to correctly determine a countermeasure to be taken for the PET apparatus100. As a result, the user is able to take the countermeasure to ensure that safety and capabilities of the PET apparatus100are guaranteed in the clinical field. Also, the user is able to properly allocate costs in the field of manufacture or services.

The PET apparatus100according to the first embodiment has thus been explained. As explained above, the PET apparatus100according to the first embodiment is able to quantitatively and objectively evaluate the impacts on the image quality (or the image) imposed by the defective channels.

Further, after the PET apparatus100has performed the evaluating process presented inFIG.4, the reconstruction processing function106bmay reconstruct PET image data by using, instead of the sensitivity map130a, the sensitivity map generated at step S103that reflects the state of the channels of the PET detector101more accurately than the sensitivity map130adoes. For example, on the basis of the sensitivity map generated at step S103, the reconstruction processing function106bmay reconstruct the PET image data according to Expression (2) presented below.

In Expression (2), the element S′jdenotes the pixel value of the pixel in a sensitivity map S′ generated at step S103reflecting the state of the channels of the PET detector101more accurately than the sensitivity map130adoes. In this manner, the reconstruction processing function106bis configured to reconstruct the PET image data by using the reciprocal of the sensitivity of each of the pixels indicated in the sensitivity map generated at step S103. As a result, the reconstruction processing function106breconstructs the PET image data by using the sensitivity map reflecting the state of the channels of the PET detector101more accurately than the sensitivity map130adoes. It is therefore possible to reconstruct the PET image data with a higher level of precision.

Alternatively, when it is determined at step S105that the sensitivity of the PET detector101does not satisfy the predetermined criterion (step S105: No), the reconstruction processing function106bmay reconstruct PET image data by using the sensitivity map generated at step S103. When the technique is applied to the state in which the PET apparatus100has defective channels, interpreting the PET image is made easier because uniformity of the image is restored by re-normalization that reflects the defective channels.

A Modification Example of the First Embodiment

After the PET apparatus100has performed the evaluating process (an evaluating process for the first time) presented inFIG.4, there are some situations in which the user replaces the defective channels with normal channels or the user changes the positions of the defective channels. In those situations, by performing the evaluating process again (a second-time evaluating process), the PET apparatus100may generate, at step S103in the second-time evaluating process, a sensitivity map having enhanced sensitivity compared to the sensitivity map generated at step S103in the first-time evaluating process. Thus, this modification example will be explained as a modification example of the first embodiment. In the description of the modification example of the first embodiment, differences from the first embodiment will primarily be explained. Explanations of some of the elements that are the same as those in the first embodiment may be omitted.

In the modification example of the first embodiment, when the message is displayed at step S107indicating that the sensitivity of the PET detector101does not satisfy the predetermined criterion, as a result of the PET apparatus100performing the evaluating process for the first time, the user acknowledges this message. Thus, the user understands that the sensitivity of the PET detector101does not satisfy the predetermined criterion.

In this situation, the user may take any of various countermeasures. For example, in a first example of the countermeasures, to make all the channels normal channels, the user may replace the sets that are each made up of a scintillator and a SiPM and correspond to the defective channels, with normal sets each made up of a scintillator and a SiPM having no failure. In other words, the user replaces the defective channels with the normal channels. In another example, as a second example of the countermeasures, while some channels remain as defective channels, the user may change the positions of the sets that are each made up of a scintillator and a SiPM and correspond to the defective channels to different positions, so that the sensitivity of the PET detector101satisfies the predetermined criterion, and also change the positions of the normal sets that are each made up of a scintillator and a SiPM having no failure and used to be arranged in the abovementioned different positions, to the pre-change positions of the defective channels. In other words, the user interchanges the positions between the defective channels and the normal channels.

After that, the PET apparatus100performs the evaluating process for the second time. For example, when the user inputs, to the processing circuitry106, an execution instruction to perform the second-time evaluating process via the input interface110, the PET apparatus100performs the evaluating process for the second time. Further, at step S106in the second-time evaluating process, when a message is displayed indicating that the sensitivity of the PET detector101satisfies the predetermined criterion, the user acknowledges this message. Thus, the user understands that the sensitivity of the PET detector101satisfies the predetermined criterion. Subsequently, after the second-time evaluating process is performed, the reconstruction processing function106breconstructs PET image data according to Expression (2) presented above, by using the sensitivity map generated at step S103in the second-time evaluating process. In this situation, in Expression (2), the element S′jdenotes the pixel value of the pixel j in the sensitivity map S′ generated at step S103in the second-time evaluating process. In this situation, the sensitivity map generated at step S103in the second-time evaluating process is a sensitivity map having the enhanced sensitivity compared to the sensitivity map generated at step S103in the first-time evaluating process. Accordingly, in the modification example of the first embodiment, it is possible to reconstruct the PET image data having excellent image quality.

Further, the sensitivity map generated at step S103in the second-time evaluating process is a sensitivity map of the PET detector101in which the defective channels have been changed with the normal channels that have no defects. Alternatively, the sensitivity map generated at step S103in the second-time evaluating process is a sensitivity map of the PET detector101in which the positions of the defective channels have been changed to the positions different from the positions of the defective channels observed at the point in time when the first-time evaluating process was performed. For instance, the point in time when the first-time evaluating process was performed is an example of the first point in time.

As explained above, in the modification example of the first embodiment, at step S101in the first-time evaluating process, the obtaining function106eobtains the failure information about the defective channels of the PET detector101corresponding to the point in time when the first-time evaluating process is performed.

After that, at step S103in the first-time evaluating process, the generating function106fgenerates the sensitivity map of the PET detector101corresponding to the point in time when the defective channels are present, on the basis of the failure information obtained at step S101in the first-time evaluating process. In this situation, for instance, the sensitivity map generated at step S103in the first-time evaluating process is an example of the first sensitivity map.

Further, at step S103in the second-time evaluating process, on the basis of the failure information obtained at step S101in the second-time evaluating process, the generating function106fgenerates the sensitivity map having the enhanced sensitivity compared to the sensitivity map generated at step S103in the first-time evaluating process. In this situation, for instance, the sensitivity map generated at step S103in the second-time evaluating process is an example of the second sensitivity map.

The PET apparatus100according to the modification example of the first embodiment has thus been explained. As explained above, the PET apparatus100according to the modification example of the first embodiment is able to reconstruct the PET image data having excellent image quality. Further, because the PET apparatus100according to the modification example of the first embodiment performs the processes similar to those performed by the PET apparatus100according to the first embodiment, it is possible to achieve advantageous effects similar to those achieved by the PET apparatus100according to the first embodiment.

Second Embodiment

In the first embodiment, the example was explained in which the PET apparatus100is configured to display the judgment result. However, the PET apparatus100may display, in addition to the judgment result, an image based on image data generated at the time of calculating the index values at step S104. Accordingly, this embodiment will be explained as a second embodiment. In the description of the second embodiment, differences from the first embodiment will primarily be explained. Explanations of some of the elements that are the same as those in the first embodiment may be omitted.

An example of an evaluating process performed by the PET apparatus100according to the second embodiment will be explained.FIG.7is a flowchart illustrating an example of a flow in the evaluating process performed by the PET apparatus100according to the second embodiment. In the second embodiment, for example, when the user inputs, to the processing circuitry106, an execution instruction to perform the evaluating process via the input interface110with the same timing as that in the first embodiment, the PET apparatus100performs the evaluating process. The processes at steps S101through S104and step S105inFIG.7according to the second embodiment are the same processes as those at steps S101through S104and step S105inFIG.4according to the first embodiment. Thus, explanations of the processes at steps S101through3104and step S105according to the second embodiment will be omitted.

It should be noted, however, that the calculating function106gstores, at step S104in the second embodiment, the image data generated at the time of calculating the index values, into the memory130. The image data is obtained either by dividing the sensitivity map generated at step S103by the sensitivity map130aor by subtracting the sensitivity map generated at step S103from the sensitivity map130a.

Further, as illustrated inFIG.7, at step S201performed between step S104and step S105, the obtaining function106eobtains the image data previously stored in the memory130at step S104.

When the judging function106hdetermines that the sensitivity of the PET detector101satisfies the predetermined criterion (step S105: Yes), the display controlling function106iproceeds to step S202. On the contrary, when the judging function106hdetermines that the sensitivity of the PET detector101does not satisfy the predetermined criterion (step S105: No), the display controlling function106iproceeds to step S203.

At step S202, the display controlling function106icauses the display device120to display a message indicating that the sensitivity of the PET detector101satisfies the predetermined criterion and an image based on the image data obtained at step S201. After that, the display controlling function106iends the evaluating process presented inFIG.7.

At step S203, the display controlling function106icauses the display device120to display a message indicating that the sensitivity of the PET detector101does not satisfy the predetermined criterion and an image based on the image data obtained at step S201. After that, the display controlling function106iends the evaluating process presented inFIG.7.

The PET apparatus100according to the second embodiment has thus been explained. The images displayed at step S202and S203each indicate a two-dimensional distribution of the index values. Accordingly, the PET apparatus100according to the second embodiment is able to make it possible for the user to visually understand the state of the PET detector101of the PET apparatus100. Consequently, the PET apparatus100according to the second embodiment is able to make it possible for the user to understand the state of the PET detector101in detail. Further, the PET apparatus100according to the second embodiment is able to make it possible for the user to sufficiently understand the state of the PET detector101. As a result, it is possible have the user use the PET apparatus100with the user's consent.

Further, during an inspection performed before the start of a work day or the like, macro capabilities of the PET detector101such as an energy resolution and a time resolution are usually checked. However, in addition to the macro capabilities, the PET apparatus100according to the second embodiment is able to make it possible for the user to understand micro capabilities such as the two-dimensional distribution of the index values.

Furthermore, because the PET apparatus100according to the second embodiment performs the processes similar to those performed by the PET apparatus100according to the first embodiment, it is possible to achieve advantageous effects similar to those achieved by the PET apparatus100according to the first embodiment.

Third Embodiment

In the first embodiment, the example was explained in which, when it is determined that the sensitivity of the PET detector101does not satisfy the predetermined criterion (step S105: No), the PET apparatus100is configured to display the message indicating that the sensitivity of the PET detector101does not satisfy the predetermined criterion. Further, in the second embodiment, the example was explained in which, when it is determined that the sensitivity of the PET detector101does not satisfy the predetermined criterion (step S105: No), the PET apparatus100is configured to display the message indicating that the sensitivity of the PET detector101does not satisfy the predetermined criterion and the image based on the image data obtained at step S201. Further, another arrangement is also acceptable in which, when it is determined that the sensitivity of the PET detector101does not satisfy the predetermined criterion (step S105: No), the PET apparatus100is configured to search for positional arrangements of the defective channels corresponding to a time when sensitivity of the PET apparatus100satisfies the predetermined criterion and to further cause the display device120to display information indicating the positional arrangements of the defective channels found in the search. Thus, this embodiment will be explained as a third embodiment.

As the PET apparatus100according to the third embodiment, the following will describe a PET apparatus100configured to perform processes at steps S301and S302described below, in place of the process at step S107in the first embodiment. However, the PET apparatus100according to the third embodiment may alternatively be configured to perform the processes at steps S301and S302described below, in place of the process at step S203in the second embodiment.

Further, in the description of the third embodiment, differences from the first embodiment will primarily be explained. Explanations of some of the constituent elements that are the same as those in the first embodiment may be omitted.

An example of the evaluating process performed by the PET apparatus100according to the third embodiment will be explained.FIG.8is a flowchart illustrating an example of a flow in the evaluating process performed by the PET apparatus100according to the third embodiment. In the third embodiment, for example, when the user inputs, to the processing circuitry106, an execution instruction to perform the evaluating process via the input interface110with the same timing as that in the first embodiment, the PET apparatus100performs the evaluating process. The processes at steps S101through S106inFIG.8according to the third embodiment are the same processes as those at steps S101through S106inFIG.4according to the first embodiment. Thus, explanations of the processes at steps S101through S106according to the third embodiment will be omitted.

As illustrated inFIG.8, in the third embodiment, when it is determined that the sensitivity of the PET detector101does not satisfy the predetermined criterion (step S105: No), the processing circuitry106proceeds to step S301.

At step S301, the processing circuitry106performs an identifying process to identify the positions (the positional arrangements) of the defective channels corresponding to the time when the sensitivity of the PET detector101satisfies the predetermined criterion.FIG.9is a flowchart illustrating an example of a flow in the identifying process performed at step S301by the PET apparatus100according to the third embodiment.

As illustrated inFIG.9, at step S303, the generating function106fsets a variable K to “1”.

Subsequently, at step S304, the generating function106fchanges the positions (the positional arrangements) of the defective channels indicated in the failure information, within the model30of the PET detector101. For example, the generating function106frandomly changes the positions of the defective channels. More specifically, for example, the generating function106frandomly generates random numbers from information about the present time or the like. In this situation, the generated plurality of random numbers are in one-to-one correspondence with a plurality of positions in which the channels can possibly be arranged. Accordingly, the generating function106fidentifies the positions corresponding to the generated random numbers. Further, the generating function106fsets the defective channels in the identified positions. Thus, the positions of the defective channels have been changed. After that, among the plurality of positions in which the channels can possibly be arranged, the generating function106fsets normal channels in the channels where no defective channels are set. In this manner, in the model of the PET detector101, the generating function106fhas changed the positions of the defective channels without changing the quantity of the defective channels.

Subsequently, at step S305, the generating function106fdraws a LOR32connecting two channels among the plurality of normal channels other than the defective channels, with respect to every combination each made up of two channels. Further, the generating function106fgenerates a sensitivity map, by converting the LOR pattern in the reconstruction region33of the model30, into image data expressing the density of the LORs for each of the pixels. In other words, on the basis of the LOR pattern in the reconstruction region33of the model30, the generating function106fcalculates the density of the LORs for each of the pixels, as sensitivity of the PET detector101, and generates the image data in which each of the pixels indicates the sensitivity, as the sensitivity map.

In this manner, when it is determined that the sensitivity of the PET detector101does not satisfy the predetermined criterion (step S105: No), the generating function106fchanges the positions of the defective channels in the PET detector101at step S303and generates the sensitivity map of the PET detector101on the basis of the post-change positions of the defective channels at step S305. For instance, the sensitivity map generated at step S305is an example of the third sensitivity map.

After that, at step S306, the calculating function106gcalculates index values related to a degree of impacts on image quality imposed by the defective channels, on the basis of the sensitivity map130aand the sensitivity map generated at step S305. For example, the calculating function106ggenerates image data by dividing the sensitivity map generated at step S305by the sensitivity map130a, while using the same method as the method previously used at step S104for generating the image data. In other words, the calculating function106gcalculates the index value for each of the pixels, on the basis of the ratios between the sensitivity map130aand the sensitivity map generated at step S305.

Alternatively, the calculating function106gmay calculate an index value for each of the pixels, on the basis of differences between the sensitivity map130aand the sensitivity map generated at step S305. For example, the calculating function106gmay generate image data by subtracting the sensitivity map generated at step S305from the sensitivity map130a, while using the same method as the method previously used at step S104for generating the image data by subtracting the sensitivity map generated at step S103from the sensitivity map130a.

As explained above, at step S306, the calculating function106gcalculates the index values related to the degree of impacts on the image quality imposed by the defective channels of which the positions were changed, on the basis of the sensitivity map130aand the sensitivity map generated at step S305.

Subsequently, at step S307, the judging function106hjudges whether or not the sensitivity of the PET detector101satisfies the predetermined criterion, on the basis of the index values. For example, the judging function106hjudges whether or not the sensitivity of the PET detector101satisfies the predetermined criterion on the basis of the index values, while using the same method as the method previously used at step S105for judging whether or not the sensitivity of the PET detector101satisfies the predetermined criterion.

As explained above, the judging function106hjudges whether or not the sensitivity of the PET detector101satisfies the predetermined criterion, on the basis of the index values calculated on the basis of the sensitivity map130aand the sensitivity map generated at step S305.

When it is determined that the sensitivity of the PET detector101does not satisfy the predetermined criterion (step S307: No), the generating function106fproceeds to step S308. On the contrary, when it is determined that the sensitivity of the PET detector101satisfies the predetermined criterion (step S307: Yes), the generating function106fproceeds to step S310.

At step S308, the generating function106fincrements the value of the variable K by 1.

Subsequently, at step S309, the generating function106fjudges whether or not the value of the variable K is equal to or larger than a threshold value Th1. The threshold value Th1is a positive integer value used for determining an upper limit value for the number of times the processes at steps S303through S308are repeatedly performed. In the present embodiment, the processes at steps3303through S308may repeatedly be performed up to (“the threshold value Th1”−1) times at most.

When the value of the variable K is equal to or larger than the threshold value Th1(step S309: Yes), the generating function106fmakes a return. On the contrary, when the value of the variable K is smaller than the threshold value Th1(step S309: No), the generating function106freturns to step S304described earlier and performs the processes at step S304and thereafter.

At step S310, the generating function106fstores, in the memory130, information indicating the positions (the positional arrangements) of the defective channels corresponding to the time when the sensitivity of the PET detector101satisfies the predetermined criterion. For example, the generating function106fgenerates information indicating the positions of the defective channels that were changed in the most recent process at step S304and stores the generated information into the memory130. After that, the generating function106fmakes a return.

In this situation, when the judgment at step S309is in the affirmative and the process makes a return, the memory130does not store therein the information indicating the positions of the defective channels corresponding to the time when the sensitivity of the PET detector101satisfies the predetermined criterion. On the contrary, when the judgment at step S307is in the affirmative and the process makes a return, the memory130stores therein the information indicating the positions of the defective channels corresponding to the time when the sensitivity of the PET detector101satisfies the predetermined criterion.

Returning to the description ofFIG.8, when the memory130has stored therein the information indicating the positions of the defective channels corresponding to the time when the sensitivity of the PET detector101satisfies the predetermined criterion, the display controlling function106icauses the display device120to display, at step S302as illustrated inFIG.8, a message indicating that the sensitivity of the PET detector101does not satisfy the predetermined criterion and the information indicating the positions of the defective channels corresponding to the time when the sensitivity of the PET detector101satisfies the predetermined criterion. After that, the display controlling function106iends the evaluating process presented inFIG.8.

As explained above, the display controlling function106icauses the display device120to display the information indicating the post-change positions of the defective channels in the PET detector101corresponding to the time when the judging function106hdetermines that the sensitivity of the PET detector101satisfies the predetermined criterion.

On the contrary, when the information indicating the positional arrangements of the defective channels corresponding to the time when the sensitivity of the PET detector101satisfies the predetermined criterion is not stored in the memory130, the display controlling function106icauses the display device120to display, at step S302, a message indicating that the sensitivity of the PET detector101does not satisfy the predetermined criterion. After that, the display controlling function106iends the evaluating process presented inFIG.8.

The PET apparatus100according to the third embodiment has thus been explained. As a countermeasure against the defective channels, the PET apparatus100according to the third embodiment makes it possible for the user to understand the positional arrangements of the defective channels that alleviate the impacts on the image quality. As a result, it is possible to suppress clinical downtime that may be caused by repair work performed by the user.

Further, because the PET apparatus100according to the third embodiment performs the processes similar to those performed by the PET apparatus100according to the first embodiment or the second embodiment, it is possible to achieve advantageous effects similar to those achieved by the PET apparatus100according to the first embodiment or the second embodiment.

Alternatively, after the PET apparatus100has performed the evaluating process presented inFIG.8, the reconstruction processing function106bmay reconstruct PET image data by using, instead of the sensitivity map130a, the sensitivity map generated at step S305corresponding to the time when the judging function106hdetermines that the sensitivity of the PET detector101satisfies the predetermined criterion. For example, on the basis of the sensitivity map generated at step S305corresponding to the time when the sensitivity of the PET detector101is determined to satisfy the predetermined criterion, the reconstruction processing function106bmay reconstruct the PET image data according to Expression (2) presented above. In that situation, in Expression (2), the element S′jdenotes the pixel value of the pixel j in the sensitivity map S′ generated at step S305corresponding to the time when the sensitivity of the PET detector101is determined to satisfy the predetermined criterion.

Fourth Embodiment

Alternatively, in the embodiments and the modification examples described above, the PET apparatus100may be configured to further determine the quantity of defective channels tolerated in the PET detector101. The tolerable quantity of defective channels may be used, for example, during an acceptance test performed at the time of accepting the PET detector101, or the like. Thus, this embodiment will be explained as a fourth embodiment.

The PET apparatus100according to the fourth embodiment is configured to perform the various processes described in any of the embodiments and the modification examples presented above. Further, the PET apparatus100according to the fourth embodiment is configured to further perform a tolerable defective channel quantity determining process. An example of the tolerable defective channel quantity determining process performed by the PET apparatus100according to the fourth embodiment will be explained.FIG.10is a flowchart illustrating an example of a flow in the tolerable defective channel quantity determining process performed by the PET apparatus100according to the fourth embodiment. In the fourth embodiment, for example, when the user inputs, to the processing circuitry106, an instruction to perform the tolerable defective channel quantity determining process via the input interface110, the PET apparatus100performs the tolerable defective channel quantity determining process.

As illustrated inFIG.10, at step S401, the obtaining function106esets a variable n and a variable m each to “1”. The variable n denotes the quantity of defective channels in a model of one detector module. For example, the model30of the PET detector101includes models of the48detector modules. Accordingly, the model30of the PET detector101includes defective channels of which the quantity is equal to (48×n).

Subsequently, at step S402, the obtaining function106eobtains, from the memory130, the sensitivity map130aof the PET detector101corresponding to the time when all the channels of the PET detector101are normal channels.

After that, at step S403, the generating function106frandomly determines the position (a positional arrangement) of each of the “48×n” defective channels in the model30of the PET detector101. More specifically, for example, the generating function106frandomly generates random numbers from information about the present time and identifiers identifying the detector modules or the like. In this situation, the generated plurality of random numbers are in one-to-one correspondence with the plurality of positions in which the channels can possibly be arranged. Further, the generating function106fidentifies the position of each of the (48×n) defective channels corresponding to a pertinent one of the generated random numbers. After that, the generating function106fsets a defective channel in a corresponding one of the identified (48×n) positions for the defective channels. As a result, the positions of the defective channels have been determined. Further, among the plurality of positions in which the channels can possibly be arranged, the generating function106fsets normal channels in certain positions that are not set with the defective channels. In this manner, the generating function106fdetermines, in the model of the PET detector101, the positions of the defective channels without changing the quantity of the defective channels.

Subsequently, at step S404, the generating function106fgenerates a sensitivity map on the basis of the position of the defective channels. The following will describe an example of the process performed at step S404. To begin with, the generating function106fdraws a LOR32connecting two channels among the plurality of normal channels other than the defective channels, with respect to every combination each made up of two channels. Further, the generating function106fgenerates the sensitivity map, by converting the LOR pattern in the reconstruction region33of the model30, into image data expressing the density of the LORs for each of the pixels. In other words, on the basis of the LOR pattern in the reconstruction region33of the model30, the generating function106fcalculates the density of the LORs for each of the pixels, as sensitivity of the PET detector101, and generates the image data in which each of the pixels indicates the sensitivity, as the sensitivity map.

Subsequently, at step S405, the calculating function106gcalculates index values related to a degree of impacts on the image quality imposed by the defective channels, on the basis of the sensitivity map130aand the sensitivity map generated at step S404. For example, the calculating function106ggenerates the image data by dividing the sensitivity map generated at step S404by the sensitivity map130a, while using the same method as the method previously used at step S104for generating the image data. In other words, the calculating function106gcalculates the index value for each of the pixels, on the basis of the ratios between the sensitivity map130aand the sensitivity map generated at step S404.

Alternatively, the calculating function106gmay calculate an index value for each of the pixels, on the basis of differences between the sensitivity map130aand the sensitivity map generated at step S404. In an example, the calculating function106gmay generate image data by subtracting the sensitivity map generated at step S404from the sensitivity map130a, while using the same method as the method previously used at step S104for generating the image data by subtracting the sensitivity map generated at step S103from the sensitivity map130a.

Subsequently, at step S406, the judging function106hjudges whether or not the sensitivity of the PET detector101satisfies the predetermined criterion, on the basis of the index values. For example, the judging function106hJudges whether or not the sensitivity of the PET detector101satisfies the predetermined criterion on the basis of the index values, while using the same method as the method previously used at step S105for judging whether or not the sensitivity of the PET detector101satisfies the predetermined criterion.

When it is determined that the sensitivity of the PET detector101satisfies the predetermined criterion (step S406: Yes), the generating function106fproceeds to step S408.

At step S408, the generating function106fjudges whether or not the value of the variable m is equal to or larger than a threshold value Th2. When the value of the variable m is equal to or larger than the threshold value Th2(step S408: Yes), the generating function106fproceeds to step S409. On the contrary, when the value of the variable m is smaller than the threshold value Th2(step S408: No), the generating function106fproceeds to step S410.

At step S409, the generating function106fincrements the value of the variable n by 1 and sets the value of the variable m to “1”. After that, the generating function106freturns to step S403and performs the processes at step S403and thereafter.

At step S410, the generating function106fincrements the value of the variable m by 1. After that, the generating function106freturns to step S403and performs the processes at step S403and thereafter.

Accordingly, while the value of the variable n is constant, the processes at steps3403through3406and steps S408through S410are repeatedly performed up to (Th2−1) times at most.

In this manner, at steps S403through S406and steps S408through S410, the processing circuitry106according to the fourth embodiment generates a plurality of sensitivity maps corresponding to a plurality of defective channel arrangement patterns having mutually-different defective channel positions, by performing, multiple times, a process of changing the positions of the defective channels either in the PET detector101or in the plurality of detector modules included in the PET detector101and generating the sensitivity map of the PET detector101on the basis of the post-change positions of the defective channels, and further repeatedly performs the abovementioned sensitivity map generating process while increasing the quantity of the defective channels until the sensitivity of the PET detector101no longer satisfies the predetermined criterion in at least one of the plurality of defective channel arrangement patterns.

On the contrary, when it is determined that the sensitivity of the PET detector101does not satisfy the predetermined criterion (step S406: No), the generating function106fproceeds to step S407.

At step S407, the judging function106hdetermines the tolerable quantity of defective channels in each of the detector modules as (n−1) and determines the tolerable quantity of defective channels in the entirety of the PET detector101as (48×(n−1)). After that, the judging function106hends the tolerable defective channel quantity determining process presented inFIG.10.

As explained above, at step S407, the processing circuitry106according to the fourth embodiment determines the quantity of defective channels tolerated either in the PET detector101or in each of the plurality of detector modules, on the basis of the quantity of the defective channels corresponding to the time when the sensitivity of the PET detector101no longer satisfies the predetermined criterion in at least one of the arrangement patterns.

The tolerable quantity of defective channels (48×(n−1)) determined in the tolerable defective channel quantity determined process is used, when the PET detector101is ordered from an external source so that an acceptance test is performed at the time of accepting the PET detector101from the external source or the like. For example, when the quantity of defective channels in each of the detector modules of which the external source notified the user is n′ while n′ is equal to or smaller than (n−1), the user will accept the PET detector101manufactured by the external source. On the contrary, when n′ is larger than (n−1), the user will not accept the PET detector101manufactured by the external source. After that, the user will re-order a PET detector101in which the quantity of defective channels in each of the detector modules is (n−1) or smaller.

In another example, when the quantity of defective channels in the PET detector101of which the external source notified the user is N while N is equal to or smaller than (48×(n−1)), the user will accept the PET detector101manufactured by the external source. On the contrary, when N is larger than (48×(n−1)), the user will not accept the PET detector101manufactured by the external source. After that, the user will re-order a PET detector101in which the quantity of defective channel is (48×(n−1)) or smaller.

The PET apparatus100according to the fourth embodiment has thus been explained. The PET apparatus100according to the fourth embodiment is able to determine the tolerable quantity of defective channels either in the PET detector101by itself or each of the detector modules by itself.

Further, the PET apparatus100according to the fourth embodiment performs the processes similar to the processes performed by the PET apparatus100according to any of the first embodiment, the second embodiment, the third embodiment, and the modification example of the first embodiment. Consequently, the PET apparatus100according to the fourth embodiment achieves advantageous effects similar to those achieved by the PET apparatus100according to any of the first embodiment, the second embodiment, the third embodiment, and the modification example of the first embodiment.

At steps S106, S107, S202, S203, and S302, the display controlling function106imay be configured to further cause the display device120to display the sensitivity map generated at step S103. The sensitivity indicated in the sensitivity map generated at step S103is closer to the actual state of the PET detector101than the sensitivity indicated in the sensitivity map130ais. For this reason, by causing the display device120to display the sensitivity map generated at step S103, it is possible to enable the user to understand the sensitivity closer to the actual state of the PET detector101. Further, because the sensitivity map130ais stored in the memory130in advance, the user is able to understand the distribution of sensitivity indicated in the sensitivity map130awith desired timing. Consequently, the user is able to understand the sections having different sensitivity levels by comparing the sensitivity map130awith the sensitivity map generated at step S103. Thus, the user is able to understand a range in which the impacts of the defective channels are exerted. Consequently, causing the display device120to display the sensitivity map generated at step S103makes it possible to quantitatively and objectively evaluate the impacts on the image quality imposed by the defective channels.

Further, at step S106, S107, S202, S203, or S302, the display controlling function106imay be configured to further cause the display device120to display the sensitivity map130aobtained at step S102and the sensitivity map generated at step S103so as to be arranged side by side to make comparison possible. Causing the sensitivity map130aand the sensitivity map generated at step S103to be displayed in a comparable manner makes it possible for the user to understand the sections having mutually-different sensitivity levels. Consequently, it is possible to enable the user to understand a range in which the impacts of the defective channels are exerted.

Further, in the embodiments and the modification examples described above, the example was explained in which the PET apparatus100performs the various types of processes by using the two-dimensional sensitivity map. However, the PET apparatus100may perform the same processes by using a three-dimensional sensitivity map.

The program executed by the one or more processors are provided as being incorporated, in advance, in a Read Only Memory (ROM), a storage circuit, or the like. Further, the program may be provided as being recorded in a non-transitory computer-readable storage medium such as a Compact Disk Read-Only Memory (CD-ROM), a Flexible Disk (FD), a Compact Disk Recordable (CD-R), or a Digital Versatile Disk (DVD), in a file in a format that is installable or executable by those devices. Further, the program may be stored in a computer connected to a network such as the Internet so as to be provided or distributed as being downloaded via the network. For example, the program is structured with modules including the processing functions described above. In actual hardware, as a result of a CPU reading and executing the program from a storage medium such as a ROM, the modules are loaded into a main storage device and generated in the main storage device.

According to at least one aspect of the embodiments described above, it is possible to quantitatively and objectively evaluate the impacts on the image quality imposed by the defective channels.