A method and system for charged particle CT imaging

Abstract

The application relates to a kind of charged particle CT imaging detection method and system. Including: obtaining the beam parameter of charged particle beam, and calculating the relative position of the charged particle beam out of the beam parameter and the adjustment position of each detector, the detector includes the upstream detector before the charged particle beam passes through the object to be measured and the digital detector after passing through the object to be measured;The position of the upstream detector and the digital detector is controlled to the corresponding adjustment position;Receive the charged particle beam, sequentially through the upstream detector, the object to be measured and the digital detector, and read the detection information of the upstream detector and the digital detector in the form of digital signal;According to the detection information, generate Sinogram image, and carry out preliminary image processing;According to the Sinogram image of preliminary image processing, generate tomographic image based on image reconstruction algorithm.The application can improve data processing efficiency and accurately measure particle trajectory.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of radiotherapy, in particular to a charged particle CT imaging detection method and system. BACKGROUND

[0002] At present, charged particle radiotherapy such as proton and heavy ion is an advanced tumor treatment method, which can achieve precise dose delivery to the tumor target area while protecting normal tissues as much as possible by using Bragg peak effect. In the process of making a proton or heavy ion treatment plan, the relative stopping power (RSP) of heavy ions in the human body is usually obtained by converting the HU value of the collected patient CT image. However, the CT image is obtained based on X-rays, and the interaction between heavy ions and matter is greatly different from the interaction between photons and matter. Therefore, the conversion process inevitably introduces uncertainty, which directly leads to inaccurate heavy ion range calculation and reduces the accuracy of heavy ion treatment, thereby affecting the treatment effect.

[0003] The present application relates to the field of radiotherapy, in particular to a charged particle CT imaging detection method and system. SUMMARY

[0004] To solve the above problems, the present application provides a charged particle CT imaging detection method and system, which uses a digital measurement scheme, unifies data processing methods, improves data processing efficiency, accurately measures the trajectory of charged particles, and accurately obtains the RSP of the measured object.

[0005] To achieve the above-mentioned purpose, the present application adopts the following technical solutions:

[0006] In a first aspect, the present application provides a charged particle CT imaging detection method, comprising:

[0007] Obtaining the beam parameters of the charged particle beam, and calculating the relative position of the charged particle beam outflow and the adjustment position of each detector according to the beam parameters, wherein the detector includes an upstream detector before the charged particle beam passes through the measured object and a digital detector after the charged particle beam passes through the measured object;

[0008] Controlling the positions of the upstream detector and the digital detector to the corresponding adjustment positions;

[0009] receive the charged particle beam, sequentially pass through the upstream detector, the object to be measured and the digital detector, and read the detection information of the upstream detector and the digital detector in the form of a digital signal;

[0010] generate a Sinogram image according to the detection information, and perform preliminary image processing;

[0011] generate a tomographic image based on an image reconstruction algorithm according to the Sinogram image after preliminary image processing.

[0012] In an implementation scheme of the present application, the beam current parameters of the charged particle beam are obtained by one or a combination of the following:

[0013] The relative position coordinates or absolute position coordinates of the outgoing beam current of the charged particle beam are obtained;

[0014] The energy of the outgoing beam current of the charged particle beam is obtained;

[0015] The full width at half maximum of the outgoing beam current of the charged particle beam is obtained;

[0016] The beam spot position of the outgoing beam current of the charged particle beam is obtained;

[0017] The deflection angle of the outgoing beam current of the charged particle beam is obtained;

[0018] The particle information of the outgoing beam current of the charged particle beam is obtained;

[0019] The frequency of the outgoing beam current of the charged particle beam is obtained;

[0020] The momentum of the outgoing beam current of the charged particle beam is obtained; or

[0021] The cross-sectional size of the outgoing beam current of the charged particle beam is obtained.

[0022] In an implementation scheme of the present application, the method further comprises: calculating the required water equivalent compensation path according to the thickness of the object to be measured, and setting a compensation object corresponding to the equivalent water thickness, so that the charged particle beam is completely captured by the detector after passing through the object to be measured.

[0023] In an implementation scheme of the present application, the Sinogram image is generated according to the detection information, comprising:

[0024] Based on the detection information of the upstream detector and the digital detector, particles meeting a predetermined condition are screened out;

[0025] Based on the calculated relative position of the outgoing charged particle beam and the adjusted position of each detector, the corresponding point coordinates in the Sinogram are calculated.

[0026] Based on the ray exit coordinates of the charged particle beam and the case information, the ray range under each case is calculated;

[0027] Based on the ray exit coordinates and the ray range under each case, a Sinogram image is solved.

[0028] In an implementation scheme of the present application, the preliminary image processing includes Mapping, tracing, non-tracing direct imaging, filtering or smoothing.

[0029] In an implementation scheme of the present application, the generating a tomographic image based on an image reconstruction algorithm includes:

[0030] Generating a tomographic image based on a FBP algorithm, a SIRT algorithm, a SART algorithm or a machine learning algorithm.

[0031] In a second aspect, the present application provides a detection system for charged particle CT imaging, comprising:

[0032] A beam and detector control module is configured to obtain beam parameters of a charged particle beam, and calculate relative positions of the charged particle beam exit and adjustment positions of each detector based on the beam parameters, wherein the detector includes an upstream detector before the charged particle beam passes through a to-be-detected object and a digital detector after the charged particle beam passes through the to-be-detected object;

[0033] The beam and detector control module is further configured to control positions of the upstream detector and the digital detector to the corresponding adjustment positions;

[0034] A charged particle beam range detection and readout module is configured to receive detection information of digital signals of the upstream detector and the digital detector when the charged particle beam sequentially passes through the upstream detector, the to-be-detected object and the digital detector;

[0035] A data processing module is configured to generate a Sinogram image based on the detection information, and perform preliminary image processing;

[0036] An image display and post-processing module is configured to generate a tomographic image based on an image reconstruction algorithm based on the Sinogram image after preliminary image processing.

[0037] In an implementation scheme of the present application, the beam and detector control module includes:

[0038] A beam and detector position calculation unit is configured to calculate relative positions of the charged particle beam exit and adjustment positions of each detector based on the beam parameters; and

[0039] A detector position control unit is configured to control the positions of the upstream detector and the digital detector to the corresponding adjustment positions by the motion unit.

[0040] In an implementation of the present application, the digital detector comprises a downstream detector, an adaptive compensator and a detection module connected in sequence.

[0041] The adaptive compensator comprises a controller and a plurality of compensator units of the same structure connected to the controller, each of the compensator units comprising a motion unit and a plurality of compensator layers of fixed width.

[0042] The detection module comprises a plurality of detection units, each of the detection units comprising an absorption layer and a detection layer.

[0043] In an implementation of the present application, the motion unit is a motion unit driven in a pneumatic, hydraulic or electric manner.

[0044] The present application has the following advantages due to the above technical solutions: the detection method and system for charged particle CT imaging provided in the present application obtain the beam parameters of the charged particle beam, calculate the relative positions of the charged particle beam out of the beam and the adjustment positions of the respective detectors according to the beam parameters, further control the respective detectors to the corresponding adjustment positions, receive the charged particle beam, sequentially pass through the upstream detector, the object to be measured and the digital detector, read the detection information of the respective detectors in the form of digital signals, generate a Sinogram image according to the detection information, perform preliminary image processing, generate a tomographic image based on an image reconstruction algorithm according to the Sinogram image after preliminary image processing, and thus improve the data processing efficiency and more accurately measure the trajectory of the charged particle compared to the prior art. BRIEF DESCRIPTION OF DRAWINGS

[0045] Figure 1 is a structural schematic diagram of a detection system for charged particle CT imaging provided by an embodiment of the present application;

[0046] Figure 2 is a detailed structural schematic diagram of part of the detection system for charged particle CT imaging in the embodiment of the present application;

[0047] Figure 3 is a structural schematic diagram of an adaptive compensator in the embodiment of the present application;

[0048] Figure 4 is a flowchart of a detection method for charged particle CT imaging in the embodiment of the present application. DETAILED DESCRIPTION

[0049] In order to make the objects, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions of the embodiments of the present application will be described clearly and completely below with reference to the drawings of the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, rather than all the embodiments of the present application. Based on the described embodiments of the present application, all other embodiments obtained by a person of ordinary skill in the art belong to the scope of protection of the present application.

[0050] The prior art detection system has the technical problems of limited energy range of charged particle detection, low data processing efficiency and detection accuracy. The technical solutions of the present application correspondingly provide a detection method and system for charged particle CT imaging. The method comprises: obtaining beam parameters of a charged particle beam, and calculating relative positions of the charged particle beam out of the beam and adjustment positions of each detector according to the beam parameters, wherein the detector comprises an upstream detector before the charged particle beam passes through an object to be detected and a digital detector after the charged particle beam passes through the object to be detected; controlling the positions of the upstream detector and the digital detector to the corresponding adjustment positions; receiving the charged particle beam, sequentially passing through the upstream detector, the object to be detected and the digital detector, reading detection information of the upstream detector and the digital detector in the form of a digital signal; generating a sinogram image according to the detection information, and performing preliminary image processing; generating a tomographic image based on an image reconstruction algorithm according to the sinogram image of the preliminary image processing. The technical solutions of the present application can adjust the positions of the detectors according to the obtained beam parameters of the charged particle beam, so that the detectors can adaptively adjust the adjustment positions of the detection, and further read the detection information of the detectors in the form of a digital signal, perform image processing according to the detection information, and finally generate a tomographic image, thereby improving the detection range and data processing efficiency compared with the prior art, and making the RSP reconstruction accuracy of the object to be detected higher.

[0051] Reference is made to Figure 1 In an embodiment of the present application, a detection system for charged particle CT imaging (referred to as a detection system in the embodiments of the present application) is provided.

[0052] The detection system in the embodiments of the present application can be applied in a medical radiotherapy device, but is not limited thereto. The detection system can be matched with a miniaturized accelerator. According to different radiotherapy modes, the type of charged particles can be selected, for example, can include protons, heavy ion beams, etc. The miniaturized accelerator is used to accelerate charged particles and emit from a beam outlet device. The emitted charged particle beam can be applied to directly act on a human body or other to-be-detected object. The detection system is used to adjust the detector, obtain detection information before and after the particle beam passes through the to-be-detected object, process data according to the detection information, generate a tomographic image, and calculate the RSP of the to-be-detected object.

[0053] Specifically, Figure 1 The detection system 100 for charged particle CT imaging shown in the embodiments includes:

[0054] A beam and detector control module (a) is configured to obtain a beam parameter of a charged particle beam, and calculate a relative position of the charged particle beam out of the beam and an adjustment position of each detector according to the beam parameter, the detector including an upstream detector before the charged particle beam passes through a to-be-detected object and a digital detector after the charged particle beam passes through the to-be-detected object.

[0055] The beam and detector control module (a) is further configured to control the positions of the upstream detector and the digital detector to the corresponding adjustment positions.

[0056] A charged particle beam range detection and readout module (b) is configured to receive detection information of digital signals of the upstream detector and the digital detector when the charged particle beam sequentially passes through the upstream detector, the to-be-detected object and the digital detector.

[0057] A data processing module (c) is configured to generate a Sinogram image according to the detection information, and perform preliminary image processing.

[0058] An image display and post-processing module (d) is configured to generate a tomographic image based on an image reconstruction algorithm according to the Sinogram image after preliminary image processing.

[0059] The detection system for charged particle CT imaging provided in the embodiments of the present application acquires the beam parameters of the charged particle beam, calculates the relative position of the charged particle beam out of the beam and the adjustment position of each detector according to the beam parameters, further controls each detector to the corresponding adjustment position, receives the charged particle beam in turn through the upstream detector, the object to be measured and the digital detector, reads the detection information of each detector in the form of digital signal, generates a Sinogram image according to the detection information, performs preliminary image processing, and generates a tomographic image based on an image reconstruction algorithm according to the Sinogram image processed preliminarily, thereby improving the data processing efficiency compared with the prior art and more accurately measuring the trajectory of the charged particle.

[0060] The above Figures 1 to 3 The structure of the detection system is further described in a more detailed embodiment of the present application.

[0061] In the embodiments of the present application, the beam and detector control module (a) of the detection system 100 comprises

[0062] The beam and detector position calculation unit (a1) is configured to calculate the relative position of the charged particle beam out of the beam and the adjustment position of each detector according to the beam parameters; and

[0063] The detector position control unit (a2) is configured to control the positions of the upstream detector and the digital detector to the corresponding adjustment positions by the motion unit.

[0064] Specifically, the beam parameters corresponding to the beam and detector control system (a) are provided before the charged particle beam is emitted. The beam parameters include but are not limited to the following information, and can also be a combination of the following information: 1) relative position coordinates or absolute position coordinates of the emitted beam; 2) energy of the emitted beam; 3) full width at half maximum of the emitted beam; 4) beam spot position of the emitted beam; 5) deflection angle information of the emitted beam; 6) particle type of the emitted beam; 7) frequency of the emitted beam; 8) momentum of the emitted beam; 9) cross-sectional size of the emitted beam, etc.

[0065] After the beam and detector position calculation unit (a1) calculates the relative position of the beam out of the beam and the adjustment position of the detector by the above information or a combination thereof, the information is transmitted to the detector position control unit (a2), which controls the motion unit inside to adjust the detection unit to the corresponding adjustment position. The motion unit is a motion unit driven by air, liquid or electricity.

[0066] As Figure 2In the embodiment of the present application, after adjusting the position of the detector, the charged particle beam emitted by the beam outlet device is received, and the capture, detection and reading of the detection information are performed by the charged particle beam range detection and reading module (b).

[0067] Specifically, the charged particle beam passes through the upstream detector (b1), the object to be detected (b2) and the digital detector (b3) in turn.

[0068] The object to be detected (b2) can be, but is not limited to, a human body, an animal or other models.

[0069] The digital detector (b3) includes a downstream detector (b31), an adaptive compensator (b32) and a detection module (b33).

[0070] The main function of the adaptive compensator (b32) is to set a compensation object with a standard equivalent water thickness according to the thickness of the object to be detected, such as a patient's body, so that the charged particles of different energies after passing through the patient's body can be completely captured by the detector. The adaptive compensator (b32) is composed of a controller and a plurality of compensation units with the same structure. Each compensation unit is composed of a motion unit and a series of compensation layers with fixed widths. The controller is responsible for communication between the software and the compensation unit, and drives the motion unit in the compensation unit to move after receiving the instruction from the software. The motion unit can use one or a combination of pneumatic, hydraulic or electric drive. The position of the driving mechanism can be accurately read by a potentiometer or an optical system. The driving mode of the motion unit includes, but is not limited to, the following combinations: 1) cylinder and push rod; 2) hydraulic cylinder and push rod; 3) motor and screw, etc.

[0071] During the imaging process, the beam outlet device gives a charged particle pencil beam that meets the requirements of clinical treatment. In the case of known beam range or energy, the required water equivalent path of compensation is calculated according to the water equivalent thickness of the object to be detected (b2) used in clinical treatment, combined with the position information provided by the beam and detector position calculation unit (a1), the adaptive compensator (b32) is adjusted to the corresponding water equivalent thickness by the controller and the driving control motion unit.

[0072] The detection module (b33) is composed of a series of detection units composed of absorption layers and detection layers. The absorption layer can be made of metals such as aluminum, copper, tungsten and uranium, and the detection layer can be composed of silicon microstrips, pixelated digital chips, silicon photomultiplier tubes and other types of detectors through separate or splicing.

[0073] In the embodiment of the present application, the data output by the digital detector (b3) includes, but is not limited to, the following information or a combination of the following information: 1) residual range; 2) residual energy; 3) scattering angle; 4) angular momentum; 5) hit position; 6) beam relative position; 7) detector relative position, etc.

[0074] Further, the data processing module (c) receives the above data or combinations thereof, determines the position of the particle stopping, obtains the accurate range of the particle, generates the Sinogram image using the remaining range and position information, and performs preliminary image processing, such as including: 1) Mapping; 2) Tracing; 3) Non-tracing direct imaging; 4) Filtering; and 5) Smoothing, etc.

[0075] In a feasible embodiment, the specific processing process of the above-mentioned data processing module (c) can include:

[0076] (1) Based on the transmission data of the upstream detector (b1) and the digital detector (b3), the particles meeting the conditions are screened, including but not limited to the following methods, and can also be a combination of the following methods: 1) screening part exceeding a certain energy threshold; 2) screening part in a unit time window; 3) screening part of a certain charge; 4) screening part of a certain momentum; 5) screening part of a certain speed; 6) directly screening part of a certain particle type, etc.

[0077] (2) Based on the transmission data of the beam and detector position calculation unit (a1), the corresponding point coordinates in the Sinogram are calculated, including but not limited to the following methods, which can be a combination of the following methods, and the following coordinates can be relative coordinates, absolute coordinates or a combination thereof: 1) directly calculating the corresponding Sinogram coordinates based on the upstream detector coordinates; 2) calculating the corresponding Sinogram coordinates based on the coordinates of the beam out of the beam device; 3) directly calculating the corresponding Sinogram coordinates based on the digital detector coordinates; 4) calculating the corresponding Sinogram coordinates based on the coordinates of the object to be measured; 5) calculating the corresponding cases in the phantom based on the upstream detector coordinates and the digital detector coordinates, and then calculating the Sinogram coordinates; 6) calculating the corresponding cases in the phantom based on the upstream detector coordinates, the coordinates of the object to be measured, and the digital detector coordinates, and then calculating the Sinogram coordinates; 7) calculating the corresponding cases in the phantom based on the coordinates of the beam out of the beam device and the upstream detector coordinates, the coordinates of the object to be measured, and the digital detector coordinates (b3), and then calculating the Sinogram coordinates; 8) based on the combination of the above calculation schemes, taking the statistical solution or the optimization solution (including but not limited to global optimization, local optimization, Pareto optimization (Pareto Optimality), etc.) and calculating the Sinogram coordinates; 9) directly mapping the Sinogram coordinates based on other reference images, etc.

[0078] (3) Based on the ray exit coordinates and the case information, the ray range under each case is calculated, including but not limited to the following methods, which can be a combination of the following methods: 1) directly using the maximum layer depth with response under each case as the ray range; 2) using the maximum layer depth after the main track to determine the ray range; 3) using the distribution of the number of hits in each detector to determine the ray range; 4) using the momentum distribution in each detector to determine the ray range; 5) using the distribution of the number of hits in each detector to determine the ray range; 6) using the energy deposition in each detector to determine the ray range, etc.

[0079] Wherein, the maximum layer depth refers to the layer number of the detection unit arranged in sequence in the detection module (b33).

[0080] Further, the tracking algorithm includes but is not limited to the following algorithms, which can be a combination of the following algorithms, or a statistical solution or an optimized solution (including but not limited to global optimization, local optimization, Pareto optimization, etc.) of the results of the following algorithms: 1) based on the deflection angle calculated from the coordinates of the hits in each layer, the tracking is performed based on the minimum deflection angle; 2) based on the energy deposition of the hits in each layer, the tracking is performed based on the minimum energy change; 3) based on the momentum of the hits in each layer, the tracking is performed based on the minimum momentum change; 4) based on the scattering angle of the hits in each layer, the tracking is performed based on the minimum scattering angle, etc.

[0081] The above tracking schemes can be further optimized, including but not limited to the following methods, which can be a combination of the following optimization methods: 1) set a deflection angle threshold, and only track within a specific deflection angle range; 2) set an energy threshold, and only track within a specific energy change range; 3) set a momentum threshold, and only track within a specific momentum change range; 4) set a scattering angle threshold, and only track within a specific scattering angle range, etc.

[0082] (4) Based on the ray exit coordinates and the statistical value of the ray range of each coordinate under the case, the Sinogram image is solved. The calculation of the range of each coordinate includes but is not limited to the following algorithms, which can be a combination of the following algorithms, or a statistical solution or an optimized solution (including but not limited to global optimization, local optimization, Pareto optimization, etc.) of the results of the following algorithms, and the statistical values described below include but are not limited to median, maximum, minimum, average, quartile, weighted value, etc., which can be a combination of the above values: 1) based on the distribution of the range and the deflection angle of all cases, the range of the coordinate is solved by the statistical value of the peak points of the two distributions; 2) based on the distribution of the range and the energy of all cases, the range of the coordinate is solved by the statistical value of the peak points of the two distributions; 3) based on the distribution of the range and the momentum of all cases, the range of the coordinate is solved by the statistical value of the peak points of the two distributions; 4) based on the distribution of the range of all cases, the peak point is taken as the range of the coordinate, etc.

[0083] The peak value statistics solving coordinate range includes but is not limited to the following algorithms, can be a combination of the following algorithms, can be a statistical solution or an optimized solution (including but not limited to global optimization, local optimization, Pareto optimization, etc.) of the results of the following algorithms: 1) by its Histogram distribution, filter the coordinates exceeding the given threshold as its peak point, and select the maximum value of the peak point as its range; 2) by its Histogram distribution, filter the coordinates exceeding the given threshold as its peak point, and select the minimum value of the peak point as its range; 3) by its Histogram distribution, filter the coordinates exceeding the given threshold as its peak point, and select the statistical value (including but not limited to mean, median, weighted value, etc.) of the maximum of the two values as its range; 4) by its Histogram distribution, filter the coordinates exceeding the given threshold as its peak point, and select the statistical value (including but not limited to mean, median, weighted value, etc.) of the minimum of the two values as its range.

[0084] In the embodiment of the present application, the image display and post-processing module (d) in the detection system calculates the image domain image based on the image reconstruction algorithm according to the Sinogram calculated by the data processing module (c). The tomographic image of the charged particle after passing through the object to be detected generated by the image display and post-processing module (d) allows further image processing operations, including but not limited to the following operations, or a combination of the following operations: 1) automatic delineation of the target area; 2) super-resolution; 3) automatic diagnosis; and 4) statistical analysis of abnormal regions, etc.

[0085] The detection system of the embodiment of the present application has the following advantages:

[0086] 1) By increasing the adaptive adjustment compensator, the dynamic adjustment of the detection range can be realized without adjusting the detector structure, the use efficiency is improved, and the application range of the imaging system is expanded.

[0087] 2) By increasing the adaptive detection system with a motion-driven adjustment mechanism, the positions of the detector and the detection unit can be adaptively adjusted, the scanning range of the pencil beam can be synchronized, the use amount of the digitization chip can be reduced, the device can be miniaturized, and the data processing efficiency is improved while the equipment cost is significantly reduced.

[0088] 3) The charged particle track tracing algorithm can efficiently and accurately reconstruct the track, further optimizing the image quality.

[0089] In another aspect of the embodiment of the present application, a detection method for charged particle CT imaging is also provided.

[0090] For details, please refer to Figure 4The application provides a charged particle CT imaging detection method. Figures 1 to 3 The detection system in the embodiment.

[0091] The charged particle CT imaging detection method in the embodiment comprises:

[0092] S41, acquiring a beam parameter of a charged particle beam, and calculating a relative position of the charged particle beam out of the beam and an adjustment position of each detector according to the beam parameter, wherein the detector comprises an upstream detector before the charged particle beam passes through an object to be detected and a digital detector after the charged particle beam passes through the object to be detected;

[0093] S42, controlling the positions of the upstream detector and the digital detector to the corresponding adjustment positions;

[0094] S43, receiving the charged particle beam, sequentially passing through the upstream detector, the object to be detected and the digital detector, and reading detection information of the upstream detector and the digital detector in a digital signal mode;

[0095] S44, generating a sinogram image according to the detection information, and performing preliminary image processing;

[0096] S45, generating a tomographic image based on an image reconstruction algorithm according to the sinogram image after preliminary image processing.

[0097] The specific process of the charged particle CT imaging detection method can be referred to the related description in the foregoing system embodiment, and will not be repeated here.

[0098] The method in the foregoing embodiment can acquire a beam parameter of a charged particle beam, calculate a relative position of the charged particle beam out of the beam and an adjustment position of each detector according to the beam parameter, further control each detector to the corresponding adjustment position, then receive the charged particle beam, sequentially pass through the upstream detector, the object to be detected and the digital detector, read detection information of each detector in a digital signal mode, generate a sinogram image according to the detection information, perform preliminary image processing, and generate a tomographic image based on an image reconstruction algorithm according to the sinogram image after preliminary image processing, so that the data processing efficiency is improved, and the trajectory of the charged particle can be more accurately measured compared with the prior art.

[0099] In several embodiments provided in the present application, it should be understood that the disclosed system, device and method can be implemented in other manners. For example, the embodiments of the device described above are merely schematic, and the division of the units is merely logical function division. There can be another division manner for the actual implementation, for example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the displayed or discussed mutual couplings or direct couplings or communication connections between the units can be indirect couplings or communication connections through some interfaces, devices or units, and can be in electrical, mechanical or other forms.

[0100] The above merely provides preferred embodiments of the application, and is not used to limit the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A method of detecting for charged particle CT imaging, characterized by, The method comprises the following steps: acquiring beam parameters of the charged particle beam, and calculating a relative position of the charged particle beam and an adjusting position of each detector according to the beam parameters, the detectors including an upstream detector before the charged particle beam passes through a to-be-detected object and a digital detector after the charged particle beam passes through the to-be-detected object; controlling the positions of the upstream detector and the digital detector to the corresponding adjusting positions, wherein the method comprises the following steps: setting a compensation object corresponding to an equivalent water thickness according to a thickness of the to-be-detected object, so that the charged particle beam is completely captured by the detectors after passing through the to-be-detected object; receiving the charged particle beam, and sequentially passing the charged particle beam through the upstream detector, the to-be-detected object and the digital detector, and reading detection information of the upstream detector and the digital detector in the form of a digital signal; generating a sinogram image according to the detection information, and performing preliminary image processing; generating a tomographic image based on an image reconstruction algorithm according to the sinogram image after preliminary image processing; wherein the beam parameters include a combination of the following information: 1) relative position coordinates or absolute position coordinates of the outgoing beam; 2) energy of the outgoing beam; 3) full width at half maximum of the outgoing beam; 4) beam spot position of the outgoing beam; 5) deflection angle information of the outgoing beam; 6) particle type of the outgoing beam; 7) frequency of the outgoing beam; 8) momentum of the outgoing beam; and 9) cross-sectional size of the outgoing beam; in the case of known beam range or energy, calculating a required water equivalent path of compensation according to a water equivalent thickness of a to-be-detected object used in clinical application, combining position information provided by a beam and detector position calculation unit, and adjusting an adaptive compensator to a corresponding water equivalent thickness through a controller and a driving control motion unit; determining the adjusting position of each detector based on data output by the digital detector, wherein the data output by the digital detector includes a combination of the following information: 1) remaining range; 2) remaining energy; 3) scattering angle; 4) angular momentum; 5) hit position; 6) relative position of the beam; and 7) relative position of the detector.

2. The method of charged particle CT imaging according to claim 1, characterized in that The acquisition of the beam parameters of the charged particle beam comprises a combination of one or more of the following: acquiring relative position coordinates or absolute position coordinates of the outgoing beam of the charged particle beam; acquiring energy of the outgoing beam of the charged particle beam; acquiring full width at half maximum of the outgoing beam of the charged particle beam; acquiring beam spot position of the outgoing beam of the charged particle beam; acquiring deflection angle of the outgoing beam of the charged particle beam; acquiring particle information of the outgoing beam of the charged particle beam; acquiring frequency of the outgoing beam of the charged particle beam; acquiring momentum of the outgoing beam of the charged particle beam; or acquiring cross-sectional size of the outgoing beam of the charged particle beam.

3. The method of charged particle CT imaging according to claim 1, wherein, The generation of the sinogram image according to the detection information comprises the following steps: screening out particles meeting a preset condition based on the detection information of the upstream detector and the digital detector; calculating corresponding point coordinates in the sinogram based on the calculated relative position of the charged particle beam and the adjusting position of each detector. Based on the ray exit coordinates of the charged particle beam and the case information, the ray range under each case is calculated; Based on the ray exit coordinates and the ray range under each case, a Sinogram image is solved.

4. The method of charged particle CT imaging according to claim 3, characterized in that, The preliminary image processing includes Mapping, tracing, non-tracing direct imaging, filtering or smoothing.

5. The method of claim 1, wherein, The generation of the tomographic image based on the image reconstruction algorithm includes: The generation of the tomographic image based on the FBP algorithm, the SIRT algorithm, the SART algorithm or the machine learning algorithm.

6. A detection system for charged particle CT imaging, characterized by It includes: A beam and detector control module is configured to obtain beam parameters of a charged particle beam, and calculate relative positions of the charged particle beam out of the beam and adjustment positions of each detector based on the beam parameters, the detector including an upstream detector before the charged particle beam passes through a to-be-detected object and a digital detector after the charged particle beam passes through the to-be-detected object; The beam and detector control module is further configured to control positions of the upstream detector and the digital detector to the corresponding adjustment positions; A charged particle beam range detection and readout module is configured to receive detection information of digital signals of the upstream detector and the digital detector when the charged particle beam sequentially passes through the upstream detector, the to-be-detected object and the digital detector; A data processing module is configured to generate a Sinogram image based on the detection information and perform preliminary image processing; An image display and post-processing module is configured to generate a tomographic image based on an image reconstruction algorithm based on the Sinogram image after preliminary image processing; The digital detector includes a downstream detector, an adaptive compensator and a detection module connected in sequence; The adaptive compensator includes a controller and a plurality of compensator units of the same structure connected with the controller, each compensator unit including a motion unit and a plurality of compensator layers of fixed width; The detection module includes a plurality of detection units, each detection unit including an absorption layer and a detection layer; The detection system of the charged particle CT imaging realizes the method of claim 1.

7. The charged particle CT imaging detection system of claim 6, wherein, The beam and detector control module includes: A beam and detector position calculation unit is configured to calculate relative positions of the charged particle beam out of the beam and adjustment positions of each detector based on the beam parameters; and A detector position control unit is configured to control positions of the upstream detector and the digital detector to the corresponding adjustment positions by the motion unit.

8. The charged particle CT imaging detection system of claim 6, wherein, The motion unit is a motion unit driven in a pneumatic, hydraulic or electric manner.

Images