Energy spectrum acquisition method and system and storage medium

Abstract

The invention discloses an energy spectrum acquisition method and system and a storage medium. The acquisition method comprises the following steps: acquiring a first depth dose curve corresponding to first energy and a second depth dose curve corresponding to second energy; based on the first depth dose curve and the second depth dose curve, respectively acquiring a first group of dose values and a second group of dose values along the depth direction; determining a relationship between the first depth dose curve and the second depth dose curve based on the first set of dose values and the second set of dose values; acquiring a first energy spectrum of the first depth dose curve based on the first depth dose curve; and determining a second energy spectrum of the second depth dose curve based on the first energy spectrum of the first depth dose curve and a relationship between the first depth dose curve and the second depth dose curve. The energy spectrum acquisition system (200) comprises a curve acquisition module (210), a dose value acquisition module (220), a relationship determination module (230), a first energy spectrum acquisition module (240) and a second energy spectrum determination module (250).

Description

Energy spectrum acquisition method, system and storage medium Technical Field

[0001] This specification relates to the field of radiotherapy, and in particular to an energy spectrum acquisition method, system, and storage medium. Background Art

[0002] Radiation therapy is an important treatment for malignant tumors. In many radiotherapy applications, such as predicting the dose distribution within a patient, the energy spectrum of the radiation beam is crucial for accurate calculation. For example, the Monte Carlo method is used clinically to predict dose distribution. The Monte Carlo algorithm is a stochastic simulation method that accurately determines the dose distribution of particles within a phantom by stochastically simulating the interactions between particles in the radiation beam and atoms within the phantom. However, in the Monte Carlo algorithm, the energy spectrum is the foundation for accurate dose calculations.

[0003] Therefore, an energy spectrum acquisition scheme is needed to improve the calculation efficiency of the energy spectrum.

[0004] Summary of the Invention

[0005] One or more embodiments of the present specification provide an energy spectrum acquisition method, which includes: acquiring a first depth dose curve corresponding to a first energy and a second depth dose curve corresponding to a second energy; acquiring a first group of dose values ​​and a second group of dose values ​​based on the first depth dose curve and the second depth dose curve, respectively; determining a relationship between the first depth dose curve and the second depth dose curve based on the first group of dose values ​​and the second group of dose values; acquiring a first energy spectrum corresponding to the first depth dose curve; and determining a second energy spectrum of the second depth dose curve based on the first energy spectrum and the relationship between the first depth dose curve and the second depth dose curve.

[0006] One or more embodiments of the present specification provide an energy spectrum acquisition system, which includes: a curve acquisition module for acquiring a first depth dose curve corresponding to a first energy and a second depth dose curve corresponding to a second energy; a dose value acquisition module for acquiring a first group of dose values ​​and a second group of dose values, respectively, based on the first depth dose curve and the second depth dose curve; a relationship determination module for determining a relationship between the first depth dose curve and the second depth dose curve based on the first group of dose values ​​and the second group of dose values; a first energy spectrum acquisition module for acquiring a first energy spectrum corresponding to the first depth dose curve; and a second energy spectrum determination module for determining a second energy spectrum of the second depth dose curve based on the first energy spectrum and the relationship between the first depth dose curve and the second depth dose curve.

[0007] One or more embodiments of this specification provide a computer-readable storage medium, wherein the storage medium stores computer instructions. When a computer reads the computer instructions in the storage medium, the computer executes the energy spectrum acquisition method.

[0008] One or more embodiments of this specification provide a device for acquiring an energy spectrum, characterized in that the device includes at least one processor and at least one memory; the at least one memory is used to store computer instructions; and the at least one processor is used to execute at least part of the computer instructions to implement a method for acquiring an energy spectrum. BRIEF DESCRIPTION OF THE DRAWINGS

[0009] This specification will be further described in the form of exemplary embodiments, which will be described in detail with reference to the accompanying drawings. These embodiments are not limiting, and in these embodiments, like numbers represent like structures, wherein:

[0010] FIG1 is a schematic diagram of an application scenario of an energy spectrum acquisition system according to some embodiments of this specification;

[0011] FIG2 is an exemplary module diagram of an energy spectrum acquisition system according to some embodiments of this specification;

[0012] FIG3 is an exemplary flow chart of a method for acquiring an energy spectrum according to some embodiments of this specification;

[0013] FIG4 a is a schematic diagram of obtaining a first depth dose curve and / or a second depth dose curve according to some embodiments of this specification;

[0014] FIG4 b is a schematic diagram of a first depth dose curve and a second depth dose curve according to some embodiments of this specification;

[0015] FIG5 is a schematic diagram of a first depth dose curve and a second depth dose curve based on depth discretization according to some embodiments of this specification;

[0016] FIG6 is an exemplary flowchart of a method for acquiring a first energy spectrum of a first depth dose curve according to some embodiments of the present specification. DETAILED DESCRIPTION

[0017] To more clearly illustrate the technical solutions of the embodiments of this specification, the following briefly describes the drawings required for describing the embodiments. Obviously, the drawings described below are merely examples or embodiments of this specification. Those skilled in the art can apply this specification to other similar scenarios based on these drawings without inventive effort. Unless otherwise apparent from the context or otherwise noted, the same reference numerals in the figures represent the same structure or operation.

[0018] It should be understood that the terms "system," "device," "unit," and / or "module" used herein are a method for distinguishing different components, elements, parts, portions, or assemblies at different levels. However, if other terms can achieve the same purpose, the terms may be replaced by other expressions.

[0019] As used in this specification and claims, unless the context clearly indicates otherwise, the words "a," "an," "an," and / or "the" do not refer to the singular but also include the plural. Generally speaking, the terms "comprises" and "include" only indicate the inclusion of the steps and elements specifically identified, and these steps and elements do not constitute an exclusive list. A method or apparatus may also include other steps or elements.

[0020] Flowcharts are used throughout this specification to illustrate the operations performed by systems according to embodiments of this specification. It should be understood that the preceding and following operations do not necessarily need to be performed in exact order. Instead, the steps may be processed in reverse order or simultaneously. Furthermore, other operations may be added to these processes.

[0021] FIG1 is a schematic diagram of an application scenario of an energy spectrum acquisition system according to some embodiments of this specification.

[0022] The energy spectrum acquisition system can acquire the energy spectrum of a multi-energy radiation beam by implementing the methods and / or processes disclosed in this specification. As shown in FIG1 , an application scenario 100 may include: a device 110 , a processing device 120 , a terminal device 130 , a network 140 , and / or a storage device 150 .

[0023] The components of the energy spectrum acquisition system can be connected in one or more various ways. By way of example only, as shown in FIG1 , device 110 can be connected to processing device 120 via network 140. For another example, device 110 can be directly connected to processing device 120 (as indicated by the dashed double-headed arrow connecting device 110 and processing device 120). As a further example, storage device 150 can be connected to processing device 120 directly or via network 140. As a further example, terminal device 130 can be connected to processing device 120 directly (as indicated by the dashed double-headed arrow connecting terminal device 130 and processing device 120) and / or via network 140.

[0024] The device 110 can utilize one or more radiation beams to treat and / or image a target object. In some embodiments, the device 110 may include a radiation source. In some embodiments, the radiation source can emit radiation beams toward a phantom. In some embodiments, the radiation beams may include, but are not limited to, radioisotope (e.g., beta, gamma) beams, X-ray beams, photon beams, electron beams, heavy particle (e.g., proton, neutron) beams, and the like. In some embodiments, the device 110 may be a radiotherapy device that utilizes one or more radiation beams to treat a target object. Examples include an X-ray therapy machine, a gamma ray afterloader, an electron accelerator, a neutron accelerator, and a proton therapy machine. In some embodiments, the device 110 may be an imaging device that utilizes one or more radiation beams to image a target object. Examples include an X-ray imaging device, a computed tomography (CT) imaging device, a cone beam CT (CBCT) scanning device, and the like. In some embodiments, the energy spectrum acquisition system 200 may utilize measurement equipment such as an ionization chamber measurement device or a semiconductor measurement device to measure the absorbed dose of the radiation beam at different depths within the phantom.

[0025] The processing device 120 can process data and / or information obtained from the device 110, the terminal device 130 and / or the storage device 150. For example, the processing device 120 can obtain the absorbed dose of the radiation beam at different depths in the phantom from the device 110, and obtain a depth dose curve, such as a percentage depth dose (PDD) curve. For another example, the processing device 120 can discretize a first depth dose curve and a second depth dose curve, and determine the relationship between the first depth dose curve and the second depth dose curve. For another example, the processing device 120 can determine a second energy spectrum of the second depth dose curve. In some embodiments, the processing device 120 may include a central processing unit (CPU), a digital signal processor (DSP), a system on a chip (SoC), a microcontroller unit (MCU), etc. and / or any combination thereof. In some embodiments, the processing device 120 may include a computer, a user console, a single server or a server group, etc. The server group may be centralized or distributed. In some embodiments, the processing device 120 may be local or remote. For example, processing device 120 can access information and / or data stored in device 110, terminal device 130, and / or storage device 150 via network 140. For another example, processing device 120 can directly connect to device 110, terminal device 130, and / or storage device 150 to access stored information and / or data. In some embodiments, processing device 120 can be implemented on a cloud platform. By way of example only, a cloud platform can include a private cloud, a public cloud, a hybrid cloud, a community cloud, a distributed cloud, an inter-cloud, a multi-cloud, or any combination thereof. In some embodiments, processing device 120, or a portion of processing device 120, can be integrated into device 110.

[0026] The terminal device 130 can display the depth dose curve and / or the energy spectrum of the radiation beam to the user. The terminal device 130 can include a mobile device 131, a tablet computer 132, a notebook computer 133, etc., or any combination thereof. In some embodiments, the terminal device 130 can be part of the processing device 120.

[0027] Network 140 may comprise any suitable network that facilitates information and / or data exchange within the energy spectrum acquisition system. In some embodiments, one or more components of the energy spectrum acquisition system (e.g., device 110, processing device 120, storage device 150, terminal device 130) may communicate information and / or data with one or more other components of the energy spectrum acquisition system via network 140. For example, processing device 120 may obtain a monoenergetic depth dose curve from storage device 150 via the network. For another example, terminal device 130 may obtain a second energy spectrum of a second depth dose curve from processing device 120 via network 140. Network 140 may be and / or comprise a public network, a private network, a wide area network (WAN), a wired network, a wireless network, a cellular network, a frame relay network, a virtual private network, a satellite network, a telephone network, a router, a hub, a switch, a server computer, and / or any combination thereof. In some embodiments, network 140 may comprise one or more network access points. For example, the network 140 may include wired and / or wireless network access points such as base stations and / or Internet exchange points, through which one or more components of the energy spectrum acquisition system may connect to the network 140 to exchange data and / or information.

[0028] The storage device 150 can store data, instructions, and / or any other information. In some embodiments, the storage device 150 can store data obtained from the device 110, the terminal device 130, and / or the processing device 120. In some embodiments, the storage device 150 can include a mass storage device, a removable storage device, a volatile read-write memory, a read-only memory (ROM), or the like, or any combination thereof. In some embodiments, the storage device 150 can be executed on a cloud platform. In some embodiments, the storage device 150 can be connected to a network 140 to communicate with one or more other components of the energy spectrum acquisition system (e.g., the device 110, the processing device 120, the terminal device 130). One or more components of the energy spectrum acquisition system can access the data or instructions stored in the storage device 150 via the network 140. In some embodiments, the storage device 150 can be directly connected to or communicate with one or more other components of the energy spectrum acquisition system (e.g., the device 110, the processing device 120, the storage device 150, the terminal device 130). In some embodiments, the storage device 150 can be part of the processing device 120.

[0029] FIG2 is an exemplary module diagram of an energy spectrum acquisition system according to some embodiments of this specification. As shown in FIG2 , the energy spectrum acquisition system 200 may include a curve acquisition module 210 , a dose value acquisition module 220 , a relationship determination module 230 , a first energy spectrum acquisition module 240 , and a second energy spectrum determination module 250 .

[0030] Curve acquisition module 210 can be used to acquire a first depth dose curve corresponding to a first energy and a second depth dose curve corresponding to a second energy. In some embodiments, the maximum measurement depths corresponding to the first and second depth dose curves can be determined based on the thickness and / or material of the phantom. In some embodiments, the first energy can be no less than the second energy. A detailed description of curve acquisition module 210 can be found in the description of step 310 and is not repeated here.

[0031] The dose value acquisition module 220 can be configured to acquire a first set of dose values ​​and a second set of dose values ​​along the depth direction based on the first depth dose curve and the second depth dose curve, respectively. In some embodiments, the dose value acquisition module 220 can be configured to perform one or more of the following operations: determining a discrete point interval; determining discrete points based on the maximum measured depth and the discrete point interval; and discretizing the first depth dose curve and the second depth dose curve based on the discrete points to acquire a first set of dose values ​​and a second set of dose values. In some embodiments, the first set of dose values ​​and the second set of dose values ​​can be represented by a first vector and a second vector, respectively. A detailed description of the dose value acquisition module 220 can be found in the description of step 320 and is not further elaborated here.

[0032] The relationship determination module 230 can be configured to determine a relationship between the first depth dose curve and the second depth dose curve based on the first set of dose values ​​and the second set of dose values. In some embodiments, the relationship determination module 230 can determine the relationship between the first depth dose curve and the second depth dose curve based on a ratio of at least one dose value in the first set of dose values ​​and a corresponding at least one dose value in the second set of dose values. In some embodiments, the relationship between the first depth dose curve and the second depth dose curve can be represented by a diagonal matrix. In some embodiments, the non-zero elements in each row of the diagonal matrix can be determined based on a ratio of the elements in the corresponding row in the second vector and the first vector. A detailed description of the relationship determination module 230 can be found in the description of step 330 and is not further described here.

[0033] The first energy spectrum acquisition module 240 can be used to acquire a first energy spectrum corresponding to a first depth dose curve. In some embodiments, the first energy spectrum acquisition module 240 is further used to perform one or more of the following operations: determining multiple monoenergetic depth dose curves based on the first energy; acquiring a first energy spectrum of the first depth dose curve based on the first depth dose curve and the multiple monoenergetic depth dose curves. In some embodiments, the first energy spectrum acquisition module 240 is further used to perform one or more of the following operations: determining the monoenergetic energy intervals corresponding to the multiple monoenergetic depth dose curves; determining the multiple monoenergetic depth dose curves based on the first energy and the monoenergetic energy intervals. In some embodiments, the first energy spectrum acquisition module 240 is further used to perform one or more of the following operations: simulating multiple monoenergetic depth dose curves based on a Monte Carlo algorithm to acquire a phase space file; analyzing the first depth dose curve based on the phase space file to acquire a first energy spectrum of the first depth dose curve. A detailed description of the first energy spectrum acquisition module 240 can be found in the relevant description of step 340 and will not be repeated here.

[0034] The second energy spectrum determination module 250 can be used to determine the second energy spectrum of the second depth dose curve based on the first energy spectrum and the relationship between the first depth dose curve and the second depth dose curve. A detailed description of the second energy spectrum determination module 250 can be found in the relevant description of step 350 and will not be repeated here.

[0035] FIG3 is an exemplary flow chart of a method for acquiring an energy spectrum according to some embodiments of the present specification. As shown in FIG3 , process 300 includes the following steps:

[0036] Step 310 : Acquire a first depth dose curve corresponding to the first energy and a second depth dose curve corresponding to the second energy. Specifically, step 310 may be performed by the curve acquisition module 210 .

[0037] The monoenergetic ray may be a ray composed of particles of the same energy, for example, a monoenergetic ray X1 with an energy of 0.1 MeV composed of photons with an energy of 0.1 MeV.

[0038] Radiation beams used for radiotherapy and imaging can be polyenergetic rays composed of a continuous series of particles with different energies. In some embodiments, the energy of the radiation beam can be obtained by weighted summing the energies of the different particles in the radiation beam. It will be appreciated that in some embodiments, the energy of the radiation beam can be obtained by superimposing monoenergetic rays of different energies based on their corresponding weights.

[0039] The energy spectrum of a radiation beam can describe the energies and corresponding weights of the particles contained in the radiation beam. In some embodiments, the energy spectrum can be represented by a distribution curve of the particle weights in the radiation beam as a function of particle energy. Specifically, the abscissa and ordinate of the curve can represent the energy and corresponding weights of the particles in the radiation beam, respectively. For example, the energy spectrum corresponding to a radiation beam with an energy of 5 MeV can include the distribution of particle weights in the radiation beam over energies between 0 and 5 MeV.

[0040] The first radiation beam may be a radiation beam used to measure the energy spectrum of any other radiation beam.

[0041] The first energy may be the energy of the first radiation beam. For example, the first radiation beam may be an X-ray beam I1, and the first energy may be the energy of the X-ray beam I1: 10 MeV.

[0042] The second radiation beam can be any radiation beam with any energy spectrum to be measured. In some embodiments, the second radiation beam and the first radiation beam can contain the same particle type. For example, both the second radiation beam and the first radiation beam are X-ray beams. In some embodiments, the first radiation beam and the second radiation beam can originate from the same or different devices. For example, the first radiation beam and the second radiation beam can both originate from the first device 111. In another example, the first radiation beam originates from the first device 111, and the second radiation beam originates from the second device 112.

[0043] The second energy may be the energy corresponding to the second radiation beam. Continuing with the above example, the second radiation beam may be the X-ray beam I2, and the second energy may be the energy of the X-ray beam I2: 5 MeV.

[0044] In some embodiments, the first energy may be not lower than the second energy. For example, the first energy of 10 MeV may be not lower than the second energy of 5 MeV.

[0045] As can be seen from the foregoing, the first radiation beam can be used to measure the energy spectrum of the second radiation beam. Therefore, in some embodiments of the present specification, the energy range of the particles contained in the first radiation beam can be greater than the energy range of the particles contained in the second radiation beam, so that the energy of the particles contained in the first radiation beam can be used to describe the energy of the particles in the second radiation beam. Continuing with the above example, the X-ray beam I1 corresponding to the first energy of 10MeV can contain photons in the range of 0 to 10MeV, and the X-ray beam I2 corresponding to the second energy of 5MeV can contain photons in the range of 0 to 5MeV. Then, the photons in the range of 0 to 5MeV within the range of 0 to 10MeV contained in the X-ray beam I1 can be used to describe the photons contained in the X-ray beam I2.

[0046] The absorbed dose can be the radiation energy absorbed per unit mass of the phantom.

[0047] The phantom can be a model used to replace human tissue. In some embodiments, the phantom material can have the same or similar absorption and scattering properties as the human body. Examples include water, polystyrene, plexiglass, paraffin wax, and / or polyethylene. In some embodiments, the phantom can have a thickness similar to that of the human body. For example, the phantom can be a 40 cm thick water phantom.

[0048] The central axis of the radiation field within the phantom can be the central symmetry axis of the radiation beam. In some embodiments, the central axis of the radiation field within the phantom can be parallel to the thickness direction of the phantom. Figure 4a is a schematic diagram of obtaining a first depth dose curve and / or a second depth dose curve according to some embodiments of this specification. As shown in Figure 4a, the central axis OO' of the radiation field is parallel to the thickness direction of the phantom, where the thickness direction of the phantom is the depth direction.

[0049] In some embodiments, the curve acquisition module 210 may measure the absorbed dose at each depth on the central axis of the radiation field in the phantom under the same measurement conditions, respectively.

[0050] In some embodiments, the measurement conditions may include irradiation field, source-to-skin distance and measurement method, etc. The irradiation field may be the range of the radiation beam irradiating the phantom after passing through the collimator. In some embodiments, the curve acquisition module 210 may use the same collimator to obtain the same irradiation field. The source-to-skin distance (SSD) may be the distance from the radiation source to the center of the irradiation field on the surface of the phantom. In some embodiments, the curve acquisition module 210 may respectively measure the absorbed dose at each depth on the central axis of the field in the corresponding phantom at the same source-to-skin distance for the first radiation beam and the second radiation beam. In some embodiments, the measurement method may include but is not limited to ionization chamber measurement method, semiconductor measurement method, calorimetry, chemical dose method and thermoluminescence dosimeter, etc.

[0051] For example, the curve acquisition module 210 can measure the absorbed doses at different depths of the phantom under the measurement conditions of "source-skin distance equal to 30 cm, using collimator Z1, and based on ionization chamber measurement method" respectively. For example, the curve acquisition module 210 can measure the absorbed doses at different depths of the phantom under the measurement conditions of "source-skin distance equal to 30 cm, using collimator Z1, and based on semiconductor measurement method" under the measurement conditions of "source-skin distance equal to 30 cm, using collimator Z1, and based on semiconductor measurement method" rb …d ra ...corresponding absorbed dose D0 a 、D1 a …D rb a …D ra a...; Under the same measurement conditions, measure the depths d0, d1...d rb …d ra ...corresponding absorbed dose D0 b 、D1 b …D rb b …D ra b ….

[0052] The depth dose may represent the absorbed dose at different depths within the phantom. In some embodiments, the depth dose may include a percentage depth dose (PDD). A percentage depth dose curve may be a percentage of the absorbed dose at a depth on the central axis of the field within the phantom relative to the absorbed dose at a reference depth.

[0053] The reference point can be the position in the phantom along the central axis of the radiation field where the depth dose is determined to be 100%. When the energy of the radiation beam is less than 400 keV, the reference point can be located on the surface of the phantom, as shown in FIG4a , and the depth of the reference point located on the surface of the phantom is d0 = 0 mm; when the energy of the radiation beam is greater than 400 keV, the reference point can be the point of maximum absorbed dose on the central axis of the radiation field, as shown in FIG4a , and the depth of the reference point as the maximum absorbed dose point is d r .

[0054] In some embodiments, the curve acquisition module 210 may use the maximum absorbed dose position on the central axis of the radiation field in the phantom under irradiation by the first radiation beam and the second radiation beam as the first reference point and the second reference point, respectively.

[0055] Continuing with the above example, under the irradiation of the first radiation beam, the maximum absorbed dose D on the central axis of the phantom field can be ra a The corresponding position is used as the first reference point; under the irradiation of the second radiation beam, the maximum absorbed dose D on the central axis of the phantom can be rb b The corresponding position is used as the second reference point.

[0056] As another example, under the irradiation of the first radiation beam less than 400 keV, the maximum absorbed dose D0 on the central axis of the phantom field can be a The corresponding position is used as the first reference point; under the irradiation of the second radiation beam greater than 400keV, the maximum absorbed dose D on the central axis of the phantom can be rb b The corresponding position is used as the second reference point.

[0057] As can be seen from the foregoing, the selection criteria for the first reference point corresponding to the first radiation beam and the second reference point corresponding to the second radiation beam may be the same or different.

[0058] Furthermore, in some embodiments, the curve acquisition module 210 can calculate the percentage depth dose at each depth based on the absorbed dose at each depth and the absorbed dose corresponding to the reference point. For example, the depths d0, d1, ..., d rb …d ra ...The corresponding percentage depth doses are (D0 a / D ra a )×100%,(D1 a / D ra a )×100%...(D rb a / D ra a )×100%...(D ra a / D ra a )×100%...; the depths d0, d1...d on the central axis of the field of the second radiation beam rb …d ra ...The corresponding percentage depth doses are (D0 b / D rb b )×100%,(D1 b / D rb b )×100%...(D rb b / D rb b )×100%...(D ra b / D rb b )×100%...

[0059] A depth dose curve can be a curve describing how the depth dose within a phantom irradiated by a radiation beam varies with depth. In some embodiments, the depth dose curve can include a percentage depth dose (PDD) curve. Specifically, the abscissa of the PDD curve can be the depth from the phantom surface along the central axis of the radiation field within the phantom, and the ordinate can be the percentage depth dose.

[0060] The first depth dose curve can be a curve describing the variation of the depth dose of the phantom irradiated by the first radiation beam with depth. FIG4 b is a schematic diagram of the first depth dose curve and the second depth dose curve according to some embodiments of the present specification. As shown in FIG4 b , the first depth dose curve corresponding to the X-ray beam I1 (the first radiation beam) is curve a, where d ra It is the depth corresponding to the maximum dose of the phantom irradiated by the X-ray beam I1.

[0061] The second depth dose curve can be a curve describing the variation of the depth dose of the phantom irradiated by the second radiation beam with depth. As shown in FIG4b , the second depth dose curve corresponding to the X-ray beam I2 (the second radiation beam) is curve b, where d rb It is the depth corresponding to the maximum dose of the phantom irradiated by the X-ray beam I2.

[0062] The maximum measurement depth may be the maximum depth within the measurement range of the absorbed dose, that is, the maximum value corresponding to the abscissa of the depth dose curve.

[0063] It is understood that the radiation beam propagates along the thickness of the phantom and may be scattered at the bottom of the phantom (for example, when the radiation beam is incident on the phantom wall and is scattered), resulting in inaccurate measurement of the absorbed dose at the bottom of the phantom. Optionally, the maximum measurement depth can be less than the phantom thickness, for example, less than 5 cm of the phantom thickness.

[0064] In some embodiments, when the maximum measurement depths corresponding to the first depth dose curve and the second depth dose curve are different, the smaller value can be used as the maximum measurement depth. For example, if the maximum measurement depth corresponding to the first depth dose curve is 30 cm and the maximum measurement depth corresponding to the second depth dose curve is 35 cm, 30 cm is used as the maximum measurement depth for the first depth dose curve and the second depth dose curve.

[0065] Some embodiments of the present specification determine the maximum measurement depth corresponding to the first depth dose curve and the second depth dose curve based on the thickness and / or material of the phantom, which can ensure that the dose absorbed by the phantom within the measurement range is not deficient due to scattering, thereby reducing the impact of scattering on the measurement of the first depth dose curve and the second depth dose curve.

[0066] Step 320 , based on the first depth dose curve and the second depth dose curve, respectively obtain a first set of dose values ​​and a second set of dose values. Specifically, step 320 may be performed by the dose value acquisition module 220 .

[0067] In some embodiments, the dose value acquisition module 220 may discretize the first depth dose curve and the second depth dose curve to respectively acquire a first set of dose values ​​and a second set of dose values.

[0068] The discrete points may be points corresponding to discrete depths on the depth dose curve (the first depth dose curve and / or the second depth dose curve).

[0069] In some embodiments, the abscissa of the discrete point can be the discrete depth, and the ordinate can be the depth dose corresponding to the discrete depth (such as the percentage depth dose). Figure 5 is a schematic diagram of a first depth dose curve and a second depth dose curve based on depth discretization according to some embodiments of the present specification. As shown in Figure 5, the abscissa of the discrete point a1 can be the discrete depth of 1 cm, and the ordinate can be 87% of the percentage depth dose corresponding to 1 cm on the first depth dose curve a; the abscissa of the discrete point b1 can be the discrete depth of 1 cm, and the ordinate can be 97% of the percentage depth dose corresponding to 1 cm on the second depth dose curve b.

[0070] Discretization can be to represent infinite and continuous points on the depth dose curve (the first depth dose curve and / or the second depth dose curve) with finite discrete points. For example, infinite and continuous points on curve a can be represented by m discrete points on curve a.

[0071] The discrete point interval may be the interval between two adjacent discrete points on the horizontal axis, that is, the interval between two adjacent discrete points in depth.

[0072] In some embodiments, the dose value acquisition module 220 may determine the discrete point interval based on discrete accuracy and / or discrete speed.

[0073] Discretization accuracy can be a parameter used to assess the degree of discretization of the depth dose curve. It will be appreciated that a higher discretization accuracy means a greater number of discrete points representing the same depth range, and the closer the obtained first and / or second vectors are to the first and / or second depth dose curves. A detailed description of the first and second vectors is provided below and is omitted here.

[0074] In some embodiments, the discrete precision can represent the ratio of the discrete point interval to the depth measurement range. For example, a discrete precision of 0.03 indicates that the ratio of the discrete point interval to the depth measurement range is 0.03. In some embodiments, the dose value acquisition module 220 can determine the discrete point interval based on the product of the discrete precision and the depth measurement range. For example, if the discrete precision is 0.03 and the depth measurement range is 30 cm, the discrete point interval can be 0.9 cm.

[0075] The discrete speed can be a parameter for evaluating the discretization speed of the depth dose curve. In some embodiments, the discrete speed can represent the time required to complete a discretization of the depth dose curve. For example, the discrete speed can include 1s, indicating that it takes 1s to complete the depth-based discretization of the depth dose curve. In some embodiments, the dose value acquisition module 220 can determine the initial number of discrete points based on the discrete speed, thereby determining the discrete point interval. Specifically, the dose value acquisition module 220 can pre-acquire a first relationship curve between the discrete speed and the initial number of discrete points, and then determine the initial number of discrete points based on the discrete speed and the first relationship curve, and then determine the discrete point interval based on the depth measurement range and the initial number of discrete points. For example, if the discrete speed is 1s and the corresponding initial number of discrete points is 31, the discrete point interval can be determined to be 1cm by dividing the depth measurement range 30 by the difference 30 obtained by subtracting 1 from the initial number of discrete points.

[0076] In some embodiments, the dose value acquisition module 220 may average the discrete intervals determined based on the discrete precision and the discrete intervals determined based on the discrete speed to determine a final discrete interval. For example, averaging 0.9 cm and 1 cm may determine a final discrete interval of 0.95 cm.

[0077] In some embodiments, the discrete precision and / or discrete speed may be determined based on user input.

[0078] In some embodiments, the dose value acquisition module 220 may also determine a method for obtaining discrete intervals based on user-entered parameters to thereby obtain discrete intervals. For example, if the user only enters the discrete precision of "0.03," the dose value acquisition module 220 may determine the final discrete interval to be a discrete interval of 1 cm obtained based on the discrete precision.

[0079] In some embodiments, the discrete precision and / or discrete speed may be default values ​​set in advance.

[0080] In some embodiments, the dose value acquisition module 220 may obtain the discrete point interval input by the user from the terminal device 130 and determine whether the discrete point interval input by the user satisfies the discrete precision and / or discrete speed. For example, if the discrete point interval input by the user does not satisfy the discrete precision and / or discrete speed, the dose value acquisition module 220 may prompt the user to re-enter the discrete point interval or directly display a recommended discrete point interval to the user based on the discrete precision and / or discrete speed on the terminal device 130.

[0081] As can be seen from the foregoing, if the maximum measurement depth of the first depth dose curve and the second depth dose curve are the same, then the depth measurement range of the first depth dose curve and the second depth dose curve are the same. Therefore, in some embodiments, based on the same discrete accuracy and discrete speed, the discrete point intervals of the first depth dose curve and the second depth dose curve can be the same.

[0082] Furthermore, in some embodiments, the dose value acquisition module 220 may determine discrete points based on the maximum measurement depth and the discrete point interval.

[0083] Specifically, the dose value acquisition module 220 can use 0 cm as the horizontal coordinate of the first discrete point, determine multiple horizontal coordinates in sequence based on the discrete point intervals until the horizontal coordinate exceeds the maximum measurement depth, and use multiple horizontal coordinates that do not exceed the maximum measurement depth as the horizontal coordinates of the discrete points, thereby determining the discrete points of the depth dose curve (the first depth dose curve and / or the second depth dose curve).

[0084] As can be seen from the foregoing, the maximum measurement depth and discrete point spacing of the first depth dose curve and the second depth dose curve are the same. Therefore, in some embodiments, the number of discrete points and the horizontal coordinate of each discrete point in the first depth dose curve and the second depth dose curve can be the same. In some embodiments, the discretization methods of the first depth dose curve and the second depth dose curve can be the same or different, as long as the number of discrete points in the first depth dose curve and the second depth dose curve are the same.

[0085] As shown in FIG5 , 0 cm is used as the abscissa of the first discrete point. Based on the discrete point interval of 1 cm, the abscissas can be determined in sequence: 0, 1, 2, 3…30, 31 cm. Among them, 31 cm is greater than the maximum measurement depth of 30 cm. The abscissas 0, 1, 2, 3…30 cm are used as the abscissas of the first depth dose curve a and the second depth dose curve b, thereby determining m=31 discrete points of curve a: a0(0, 80), a1(1, 87),…a 30 (30, 80), and m=31 discrete points of curve b: b0(0, 100), b1(1, 97), ...b 30 (30, 45).

[0086] Furthermore, in some embodiments, the dose value acquisition module 220 may discretize the first depth dose curve and the second depth dose curve based on discrete points, respectively, to acquire a first set of dose values ​​and a second set of dose values.

[0087] The first set of dose values ​​can be a set of dose values ​​representing a first depth dose curve using dose depths corresponding to discrete points. In some embodiments, the dose values ​​of the first set of dose values ​​can be the ordinates of the discrete points on the first depth dose curve. The number of dose values ​​in the first set of dose values ​​can be the same as the number of discrete points on the first depth dose curve. For example, if the number of discrete points on the first depth dose curve is m, then the number of dose values ​​in the first set of dose values ​​is m. For example, as shown in FIG5 , the first set of dose values ​​can include the ordinates [80, 87, …, 80] corresponding to the 31 discrete points on curve a.

[0088] In some embodiments, the first set of dose values ​​can be represented by a first vector. The dimension of the first vector can be the same as the number of discrete points on the first depth dose curve. For example, if the number of discrete points on the first depth dose curve is m, the dimension of the first vector is m×1. For example, continuing with FIG. 5 as an example, the elements of the first vector can be 31 discrete points on curve a: a0(0,80), a1(1,87), ...a 30 The vertical coordinate of (30, 80) is the first vector A = [80, 87, ..., 80] T , the first vector can be a 31-dimensional vector.

[0089] The second set of dose values ​​can be a set of dose values ​​representing a second depth dose curve using dose depths corresponding to discrete points. In some embodiments, the dose values ​​of the second set of dose values ​​can be the ordinates of the discrete points on the second depth dose curve. The number of dose values ​​in the second set of dose values ​​can be the same as the number of discrete points on the second depth dose curve. For example, if the number of discrete points on the second depth dose curve is m, then the number of dose values ​​in the second set of dose values ​​is m. For example, as shown in Figure 5, the second dose values ​​can include the ordinates [100, 97, ..., 45] corresponding to the 31 discrete points on curve b.

[0090] In some embodiments, the second set of dose values ​​can be represented by a second vector. The second vector can be a vector that represents the second depth dose curve using the percentage dose depth corresponding to the discrete points. Similar to the first vector, the dimension of the second vector can be the same as the number of discrete points on the second depth dose curve. Exemplarily, if the number of discrete points on the second depth dose curve is m, the dimension of the second vector is m×1. For example, continuing with Figure 5 as an example, the elements of the second vector can be 31 discrete points on curve b: b0(0, 100), b1(1, 97), ...b 30 The ordinate of (30, 45), then the second vector B = [100, 97, ..., 45] T , the second vector can be a 31-dimensional vector.

[0091] As can be seen from the above, if the number of discrete points on the first depth dose curve and the second depth dose curve is the same, then the dimensions of the first vector and the second vector are the same. For example, if the first depth dose curve and the second depth dose curve both contain 31 discrete points, then the dimensions of the first vector A and the second vector B are both 31.

[0092] It can be understood that the smaller the discrete accuracy, the slower the discrete speed, the smaller the discrete point interval, and the greater the number of discrete points, the larger the dimension of the corresponding first vector and / or second vector, and the lower the overall computational efficiency of the system.

[0093] Some embodiments of the present specification determine the discrete point interval based on discrete precision and / or discrete speed, and determine the discrete points based on the discrete point interval. On the one hand, this can ensure the accuracy with which the first group of dose values ​​and the second group of dose values ​​represent the first depth dose curve and the second depth dose curve, respectively. On the other hand, it can limit the dimensionality of the first group of dose values ​​and the second group of dose values ​​in terms of the number of rows based on the discrete speed, thereby improving the computational efficiency of the entire system.

[0094] In some embodiments, the dose value acquisition module 220 may utilize a discrete model to acquire a first set of dose values ​​and a second set of dose values ​​based on the first depth dose curve and the second depth dose curve, respectively.

[0095] Specifically, the dose value acquisition module 220 may input the first depth dose curve and the second depth dose curve into the discrete model, and the discrete model may output the first set of dose values ​​and the second set of dose values, respectively. In some embodiments, the input of the discrete model may further include discrete accuracy and / or discrete speed.

[0096] In some embodiments, the discrete model may include but is not limited to a support vector machine model, a logistic regression model, a naive Bayes classification model, a Gaussian distribution Bayes classification model, a decision tree model, a random forest model, a KNN classification model, and a neural network model.

[0097] In some embodiments, a discrete model can be trained based on a large number of labeled training samples. Specifically, the labeled training samples are input into the discrete model, and the parameters of the discrete model are updated through training. In some embodiments, the training samples can include sample depth-dose curves. In some embodiments, the labels can be manually annotated sample arrays.

[0098] In some embodiments, the first depth dose curve and the second depth dose curve may be discretized in the same manner. For example, the first depth dose curve and the second depth dose curve may both be discretized based on discrete point intervals. In some embodiments, the first depth dose curve and the second depth dose curve may be discretized in different manners. For example, the first depth dose curve may be discretized based on discrete point intervals, and the second depth dose curve may be discretized based on a discrete model.

[0099] Step 330 : Determine the relationship between the first depth dose curve and the second depth dose curve based on the first set of dose values ​​and the second set of dose values. Specifically, step 330 may be performed by the relationship determination module 230 .

[0100] The relationship between the first depth dose curve and the second depth dose curve can be a conversion relationship between the first depth dose curve and the second depth dose curve. In some embodiments, the relationship determination module 230 can determine the relationship between the first depth dose curve and the second depth dose curve based on a ratio of at least one dose value in the first set of dose values ​​and a corresponding at least one dose value in the second set of dose values.

[0101] In some embodiments, the relationship between the first depth dose curve and the second depth dose curve may be represented by relationship data.

[0102] In some embodiments, each numerical value in the relationship data can be determined based on the ratio between the dose value corresponding to each discrete point in the first set of dose values ​​and the dose value corresponding to each discrete point in the second set of dose values. Continuing with the aforementioned example, the first numerical value in the relationship data can be the ratio (1.25) of the first dose value 100 corresponding to the first discrete point in the second set of dose values ​​[100, 97, ..., 45] to the first dose value 80 corresponding to the first set of dose values ​​[80, 87, ..., 80]. The second numerical value can be the ratio (1.11) of the second dose value 97 corresponding to the second discrete point in the second set of dose values ​​to the second dose value 87 corresponding to the first set of dose values. The mth numerical value can be the ratio (0.56) of the mth dose value 45 corresponding to the mth discrete point in the second set of dose values ​​to the mth dose value 80 corresponding to the first set of dose values. That is, the relationship data can be [1.25, 1.11, ..., 0.56].

[0103] In some embodiments, the relationship data (i.e., the relationship between the first depth dose curve and the second depth dose curve) can be represented by a diagonal matrix. Continuing with the above example, the relationship between curve a and curve b can be represented by a diagonal matrix C. In some embodiments, the relationship between the first vector, the second vector, and the diagonal matrix can be represented by formula (1).

[0104]

[0105] Among them, A is the first vector, B is the second vector, diag(C 1,1 , C 2,2 ,...,C m,m ) is a diagonal matrix.

[0106] In some embodiments, the number of rows and columns of the diagonal matrix can be the same as the dimensions of the first vector and / or the second vector. For example, if the dimensions of the first vector and the second vector are m, the dimensions of the diagonal matrix are m×m. For example, if the first vector and the second vector are both 31-dimensional vectors, the dimensions of the diagonal matrix are 31×31.

[0107] In some embodiments, the non-zero elements in each row of the diagonal matrix can be determined based on the ratio of the elements in the corresponding row in the second vector and the first vector. For example, the first vector A = [80, 87, ..., 80] T , the second vector B = [100,97,…,45] T , the non-zero elements C in the first row of the diagonal matrix 1,1 The ratio of the first row element 100 in B to the first row element 80 in the first vector A is determined to be 1.25, and the non-zero element C in the second row is 2,2 The non-zero element C in the 31st row can be determined as 1.11 based on the ratio of the 2nd row element 97 in B and the 2nd row element 87 in the first vector A. 31,31 The ratio of the 31st row element 45 in B and the first row element 80 in the first vector A can be determined to be 0.56, thereby determining the diagonal matrix C=diag(1.25, 1.11, . . . , 0.56).

[0108] Some embodiments of the present specification respectively obtain the corresponding depth doses (such as percentage depth doses) on the first depth dose curve and the second depth dose curve based on the same measurement depth interval (i.e., discrete point interval), and represent the relationship between the first vector and the second vector based on a diagonal matrix, using each element in the diagonal matrix to represent the ratio of the element in the first vector and the element in the second vector corresponding to the same horizontal coordinate, thereby simplifying the relationship between the first depth dose curve and the second depth dose curve based on a simple vector operation relationship between the first vector and the second vector, thereby improving the calculation efficiency of the subsequent second energy spectrum.

[0109] Step 340 : Based on the first depth dose curve, a first energy spectrum of the first depth dose curve is acquired. Specifically, step 340 may be executed by the first energy spectrum acquisition module 240 .

[0110] As mentioned above, the energy spectrum can be a description of the energy and weight of the particles contained in the radiation beam.

[0111] The first depth dose spectrum may be a description of the energies and first weights of particles contained in the first radiation beam.

[0112] In some implementations, the first energy spectrum acquisition module 240 may determine multiple monoenergetic depth dose curves based on the first energy; and acquire a first energy spectrum of the first depth dose curve based on the first depth dose curve and the multiple monoenergetic depth dose curves. For example, the first energy spectrum of the first depth dose curve may include first weights corresponding to 100 monoenergetic energies: 0.5% for a monoenergetic energy of 0.1 MeV, 1% for a monoenergetic energy of 0.2 MeV, ..., 0.2% for a monoenergetic energy of 10 MeV. A detailed description of acquiring the first energy spectrum of the first depth dose curve can be found in FIG. 6 and its related description and is not further elaborated here.

[0113] Step 350 : Determine a second energy spectrum of the second depth dose curve based on the first energy spectrum of the first depth dose curve and the relationship between the first depth dose curve and the second depth dose curve. Specifically, step 350 may be performed by the second energy spectrum determination module 250 .

[0114] The second energy spectrum may be a description of the monoenergetic energy contained in the second radiation beam and its corresponding second weight.

[0115] In some embodiments, the relationship between the second weight corresponding to the monoenergetic energy and the corresponding monoenergetic depth dose curve and the second depth dose curve corresponding to the second radiation beam can be expressed by formula (2):

[0116]

[0117] Among them, PDD b It can be the second depth dose curve, n can represent n single energy depth dose curves, PDD i It can be the i-th monoenergetic depth dose curve, It can be the second weight of the depth dose (e.g., percentage depth dose) corresponding to the i-th mono-energy depth dose curve in the depth dose (e.g., percentage depth dose) corresponding to the second depth dose curve. A detailed description of the mono-energy depth dose curve can be found in the relevant description of step 610 and is not repeated here.

[0118] As can be seen from formula (2), based on the second weight corresponding to the monoenergetic energy, the depth dose (such as the percentage depth dose) corresponding to the same depth on multiple monoenergetic depth dose curves is weighted and summed to obtain the depth dose (such as the percentage depth dose) corresponding to the depth on the second depth dose curve. Among them, the second weight corresponding to each monoenergetic energy can represent the contribution of the corresponding monoenergetic ray in the second radiation beam to the depth dose (such as the percentage depth dose) of the phantom.

[0119] As can be seen from the foregoing, the first energy spectrum of the first depth dose curve can be determined based on multiple first weights of the depth doses (such as percentage depth doses) corresponding to the multiple mono-energetic depth dose curves in the depth dose (such as percentage depth dose) corresponding to the first depth dose curve, and the relationship between the first depth dose curve and the second depth dose curve can be represented by the relationship between the first vector and the second vector. Therefore, in some embodiments, the second energy spectrum determination module 250 can use the multiple mono-energetic depth dose curves to represent the first vector and the second vector, respectively, so as to obtain multiple second weights of the depth doses (such as percentage depth doses) corresponding to the multiple mono-energetic depth dose curves in the depth dose (such as percentage depth dose) corresponding to the second depth dose curve based on the first energy spectrum of the first depth dose curve.

[0120] In some embodiments, the second energy spectrum determination module 250 may determine discrete points based on the same discrete point spacing as that on the first depth dose curve and / or the second depth dose curve, and discretize the multiple monoenergetic depth dose curves based on the discrete points. Specifically, the second energy spectrum determination module 250 may represent each of the n monoenergetic depth dose curves using multiple discrete points. A detailed description of the discrete points can be found in the description of step 320 and is not further elaborated here.

[0121] In some embodiments, the number of discrete points and the horizontal coordinates of each discrete point of each monoenergetic depth dose curve can be the same as those of the first depth dose curve and / or the second depth dose curve. For example, the horizontal coordinates of m=31 discrete points of the first depth dose curve and / or the second depth dose curve are 0 cm, 1 cm, ... 30 cm, respectively. Then, the second energy spectrum determination module 250 can obtain 31 discrete points corresponding to n (e.g., 100) monoenergetic depth dose curves based on the 31 horizontal coordinates. For example, the 31 discrete points on the monoenergetic depth dose curve E1 corresponding to 0.1 MeV are: e 1,0 (0, 40), e 1,1 (1, 42),…e 1,30 (30, 45), 31 discrete points on the monoenergetic depth dose curve E2 corresponding to 0.2 MeV: e 2,0 (0, 50), e 2,1 (1,54),…e 2,30 (30, 55), ....10MeV corresponding to the single energy depth dose curve E 100 31 discrete points on: e 100,0 (0,90),e 100,1 (1,92),…e 100,30 (30, 96).

[0122] The monoenergetic vector can be a vector representing a monoenergetic depth dose curve using the percentage dose depth corresponding to a discrete point. In some embodiments, the elements of the monoenergetic vector can be the ordinates of the discrete points on the monoenergetic depth dose curve, and the dimension of the monoenergetic vector can be the same as the number of discrete points on the monoenergetic depth dose curve. For example, if the number of discrete points on the monoenergetic depth dose curve is m, then the dimension of the monoenergetic vector is m×1.

[0123] In some embodiments, the second energy spectrum determination module 250 can express the monoenergetic vector corresponding to the monoenergetic depth dose curve according to formula (3):

[0124]

[0125] Among them, D i represents the monoenergetic vector corresponding to the i-th monoenergetic depth dose curve, They respectively represent the ordinates of the first, second, ...mth discrete points on the i-th monoenergetic depth dose curve (i.e., the depth dose corresponding to the discrete points).

[0126] For example, the first monoenergetic depth dose curve E1 corresponds to a monoenergetic vector D1 = [40 42 … 45] T The monoenergetic vector D2 corresponding to the second monoenergetic depth dose curve E2 is [50 54 … 55] T ...The 100th single energy depth dose curve E 100 The corresponding single energy vector D 100 =[90 92 … 96] T .

[0127] In some embodiments, the second energy spectrum determination module 250 may represent the first vector using a monoenergetic vector based on the relationship between the monoenergetic depth dose curve and the first depth dose curve corresponding to the first radiation beam. A detailed description of the relationship between the monoenergetic depth dose curve and the first depth dose curve corresponding to the first radiation beam can be found in the description of formula (8) in step 620 and is not further described here.

[0128] In some embodiments, the second energy spectrum determination module 250 may represent the first vector using a single energy vector according to formula (4).

[0129]

[0130] Where A is the first vector corresponding to the first depth dose curve, D1, D2…D n They are the monoenergetic vectors corresponding to the first monoenergetic depth dose curve, the second monoenergetic depth dose curve...the nth monoenergetic depth dose curve, They may be the first weights of the depth doses (such as percentage depth doses) corresponding to the first monoenergetic depth dose curve, the second monoenergetic depth dose curve...the nth monoenergetic depth dose curve in the depth dose (such as percentage depth dose) corresponding to the first depth dose curve.

[0131] In some embodiments, the second energy spectrum determining module 250 may represent the second vector using a monoenergetic vector based on the relationship between the monoenergetic depth dose curve and the second depth dose curve corresponding to the second radiation beam. A detailed description of the relationship between the monoenergetic depth dose curve and the second depth dose curve corresponding to the second radiation beam can be found in the description of formula (2), which is not repeated here.

[0132] In some embodiments, the second energy spectrum determination module 250 may represent the second vector using a single energy vector according to formula (5).

[0133]

[0134] Where B is the second vector corresponding to the second depth dose curve, D1, D2…D n They are the monoenergetic vectors corresponding to the first monoenergetic depth dose curve, the second monoenergetic depth dose curve...the nth monoenergetic depth dose curve, They may be the second weights of the depth doses (such as percentage depth doses) corresponding to the first mono-energy depth dose curve, the second mono-energy depth dose curve...the nth mono-energy depth dose curve in the depth dose (such as percentage depth dose) corresponding to the second depth dose curve.

[0135] Furthermore, in some embodiments, the second energy spectrum determination module 250 may represent the relationship between the first depth dose curve and the second depth dose curve using a monoenergetic vector of the monoenergetic depth dose curve. Specifically, in some embodiments, the second energy spectrum determination module 250 may substitute Formula (4) for representing the first vector using a monoenergetic vector and Formula (5) for representing the second vector using a monoenergetic vector into Formula (1) for representing the relationship between the first vector and the second vector, thereby representing the relationship between the first vector and the second vector using a monoenergetic vector according to Formula (6).

[0136]

[0137] Among them, D1, D2…D n They are the monoenergetic vectors corresponding to the first monoenergetic depth dose curve, the second monoenergetic depth dose curve...the nth monoenergetic depth dose curve; The first weights of the depth doses (such as percentage depth doses) corresponding to the first mono-energy depth dose curve, the second mono-energy depth dose curve, ... the nth mono-energy depth dose curve in the depth dose (such as percentage depth dose) corresponding to the first depth dose curve may be respectively; These may be second weights of the depth doses (e.g., percentage depth doses) corresponding to the first monoenergetic depth dose curve, the second monoenergetic depth dose curve, ..., the nth monoenergetic depth dose curve, in the depth dose (e.g., percentage depth dose) corresponding to the second depth dose curve; diag(1.25, 1.11, ..., 0.56) is a diagonal matrix C determined based on the first depth dose curve and the second depth dose curve. A detailed description of the diagonal matrix C can be found in the description of step 330 and is not repeated here.

[0138] Furthermore, in some embodiments, the second energy spectrum determination module 250 may obtain a second weight based on formula (6), thereby obtaining a second energy spectrum of the second depth dose curve. In some embodiments, the second weight may be expressed according to formula (7).

[0139]

[0140] For example, continuing the above example, the first energy spectrum of the first depth dose curve corresponding to the first energy 10MeV includes: 0.5% of the monoenergetic energy 0.1MeV, 1% of the monoenergetic energy 0.2MeV, ..., 0.2% of the monoenergetic energy 10MeV, then C=diag(C 1,1 , C 2,2 ,...,C m,m )=diag(1.25, 1.11,..., 0.56); Then the second weight corresponding to the second energy 5MeV can be obtained based on formula (7): For example, the second weight Thus, a second energy spectrum of the second depth dose curve is obtained: 0.1% monoenergetic energy 0.1 MeV, 1% monoenergetic energy 0.25 MeV, ..., 0% monoenergetic energy 10 MeV.

[0141] Similarly, the energy spectrum acquisition system 200 can acquire second energy spectra of second depth dose curves corresponding to other second energies based on the first energy spectrum of the first depth dose curve. For example, the energy spectrum acquisition system 200 can acquire second energy spectra of second depth dose curves corresponding to any second energy in the range of 0 to 10 MeV based on the first energy spectrum of the first depth dose curve corresponding to 10 MeV.

[0142] Some embodiments of this specification are based on the relationship between the quickly acquired first depth dose curve and the second depth dose curve. According to the first energy spectrum of the first depth dose curve, the second energy spectrum of the second depth dose curve can be quickly acquired, thereby improving the acquisition efficiency of the second energy spectrum.

[0143] In some embodiments, the second energy spectrum can be used to effectively assess / evaluate the impact of radiation dose on the patient.

[0144] In some embodiments, the second energy spectrum can be used for dose calculation, thereby improving the computational efficiency of predicting radiation dose distribution in a patient.

[0145] In some embodiments, the second energy spectrum can also be used for image correction. For example, the second energy spectrum can be used to accurately calculate the energy of particles before and after scattering, thereby improving the accuracy of image scattering correction.

[0146] FIG6 is an exemplary flow chart of a method for obtaining a first energy spectrum of a first depth dose curve according to some embodiments of this specification. Specifically, FIG6 can be executed by the first energy spectrum acquisition module 240. As shown in FIG6, process 600 includes the following steps:

[0147] Step 610: Determine a plurality of monoenergetic depth dose curves based on the first energy.

[0148] The mono-energetic depth dose curve may be a curve describing how the depth dose within a phantom irradiated with mono-energetic radiation varies with depth.

[0149] In some embodiments, after the first radiation beam irradiates the phantom, particles of different energies in the first radiation beam interact with atoms in the phantom, causing atoms at different depths in the phantom to absorb different energies from particles of different energies. That is, particles of different energies in the first radiation beam contribute different energies to the depth dose (or absorbed dose) at different depths in the phantom. For example, after the first radiation beam I1 irradiates the phantom, the absorbed dose D1 at a depth d1=1 cm in the phantom is a =8Gy, the corresponding percentage depth dose is 87%, wherein the absorbed dose D1 of the particles with energy of 0.1MeV in the first radiation beam to the depth d1=1cm of the phantom is a The contribution of the particle with energy of 0.2MeV is 0.05Gy, and the absorbed dose D1 at the depth of the phantom d1 = 1cm is a The contribution is 0.04 Gy... Therefore, in some embodiments, at the same depth of the phantom, the percentage depth dose corresponding to the first radiation beam can be expressed by the percentage depth dose corresponding to multiple monoenergetic rays, and the first depth dose curve of the first radiation beam can be expressed by the monoenergetic depth dose curves corresponding to multiple monoenergetic energies.

[0150] The monoenergetic energy may be the energy corresponding to the monoenergetic ray. For example, the monoenergetic energy of a monoenergetic ray composed of particles with an energy of 0.1 MeV is 0.1 MeV.

[0151] The monoenergetic energy interval may be an interval between a plurality of monoenergetic energies corresponding to a plurality of monoenergetic depth dose curves used to represent the first depth dose curve.

[0152] In some embodiments, the first energy spectrum acquisition module 240 may determine the monoenergetic energy intervals corresponding to the plurality of monoenergetic depth dose curves based on the energy spectrum accuracy and / or the energy spectrum calculation speed.

[0153] The energy spectrum accuracy may be the degree of dispersion of multiple monoenergetic energies corresponding to multiple monoenergetic depth doses over the energy range of the first radiation beam.

[0154] In some embodiments, the energy spectrum accuracy may be determined based on user input.

[0155] In some embodiments, energy spectrum accuracy can represent the ratio of the monoenergetic energy interval to the first energy corresponding to the first radiation beam. For example, an energy spectrum accuracy of 0.01 indicates that the ratio of the monoenergetic energy interval to the first energy corresponding to the first radiation beam is 0.01. It will be appreciated that a higher energy spectrum accuracy reduces the dispersion of multiple monoenergetic energies across the energy range of the first radiation beam, resulting in a greater number of monoenergetic depth dose curves used to represent the first depth dose curve, and a closer approximation of the obtained first depth dose curve to the true value.

[0156] In some embodiments, the first energy spectrum acquisition module 240 may determine the monoenergetic energy interval based on the product of the energy spectrum accuracy and the first energy of the first radiation beam. For example, if the energy spectrum accuracy is 0.01 and the first energy is 10 MeV, the monoenergetic energy interval may be 0.1 MeV.

[0157] The energy spectrum calculation speed may be a parameter for evaluating the first energy spectrum calculation speed.

[0158] In some embodiments, the energy spectrum calculation speed may represent the time required to complete the calculation of the first energy spectrum. For example, the energy spectrum calculation speed may be 2s, indicating that it takes 2 seconds to complete the calculation of the first energy spectrum. For more information on calculating the first energy spectrum, see step 620 and will not be repeated here.

[0159] In some embodiments, the energy spectrum calculation speed may be determined based on user input.

[0160] In some embodiments, the first energy spectrum acquisition module 240 can determine the initial monoenergetic energy quantity based on the energy spectrum calculation speed, thereby determining the monoenergetic energy interval. Specifically, the first energy spectrum acquisition module 240 can pre-acquire a second relationship curve between the energy spectrum calculation speed and the initial monoenergetic energy quantity, then determine the initial monoenergetic energy quantity based on the energy spectrum calculation speed and the second relationship curve input by the user, and then determine the monoenergetic energy interval based on the first energy and the initial monoenergetic energy quantity. For example, if the energy spectrum calculation speed is 2s and the corresponding initial monoenergetic energy is 50, then the monoenergetic energy interval can be determined to be 0.2MeV by dividing the first energy 10MeV by the initial monoenergetic energy quantity.

[0161] In some embodiments, the first energy spectrum acquisition module 240 may average the monoenergetic energy interval determined based on the energy spectrum accuracy and the monoenergetic energy interval determined based on the energy spectrum calculation speed to determine a final monoenergetic energy interval. For example, averaging 0.1 MeV and 0.2 MeV may determine a final monoenergetic energy interval of 0.15 MeV.

[0162] In some embodiments, the first energy spectrum acquisition module 240 may also determine a method for acquiring the monoenergetic energy interval based on user-entered parameters, thereby acquiring the monoenergetic energy interval. For example, if the user only enters the energy spectrum accuracy of "0.01," the first energy spectrum acquisition module 240 may determine the final monoenergetic energy interval to be 0.1 MeV, the monoenergetic energy interval acquired based on the energy spectrum accuracy.

[0163] In some embodiments, the energy spectrum accuracy and / or the energy spectrum calculation speed may be default values ​​set in advance.

[0164] In some embodiments, the first energy spectrum acquisition module 240 may obtain the monoenergetic energy interval input by the user from the terminal device 130 and determine whether the monoenergetic energy interval input by the user satisfies the energy spectrum accuracy and / or energy spectrum calculation speed. For example, if the monoenergetic energy interval input by the user does not meet the energy spectrum accuracy and / or energy spectrum calculation speed, the first energy spectrum acquisition module 240 may prompt the user to re-enter the interval or directly display a suggested monoenergetic energy interval based on the energy spectrum accuracy and / or energy spectrum calculation speed on the terminal device 130.

[0165] It can be understood that the lower the energy spectrum accuracy, the slower the energy spectrum calculation speed, the smaller the monoenergetic energy interval, the greater the number of monoenergetic energies, the greater the number of monoenergetic depth dose curves used to represent the first depth dose curve, and the correspondingly greater number of monoenergetic depth dose curves used to represent the second depth dose curve, the lower the overall calculation efficiency of the system and the higher the calculation accuracy. A detailed description of how to represent the second depth dose curve using a monoenergetic depth dose curve can be found in the description of step 350 and will not be repeated here.

[0166] In some embodiments, the first energy spectrum acquisition module 240 may determine a plurality of monoenergetic depth dose curves based on the first energy and the monoenergetic energy interval.

[0167] Specifically, in some embodiments, the first energy spectrum acquisition module 240 may use an energy at a distance from the monoenergetic energy of 0 MeV as the first monoenergetic energy, and then sequentially determine multiple energies based on the monoenergetic energy intervals until the energy exceeds the first energy, thereby acquiring multiple energies that do not exceed the first energy as the monoenergetic energy. For example, using 0.1 MeV as the first monoenergetic energy, and then sequentially determining multiple energies based on monoenergetic energy intervals of 0.1 MeV: 0.1 MeV, 0.2 MeV, ..., 10 MeV, 10.1 MeV, where 10.1 MeV is greater than the first energy of 10 MeV, multiple energies that do not exceed 10 MeV: 0.1 MeV, 0.2 MeV, ..., 10 MeV may be acquired as the monoenergetic energy.

[0168] Furthermore, in some embodiments, the first energy spectrum acquisition module 240 may determine multiple corresponding monoenergetic depth dose curves based on the multiple monoenergetic energies. In some embodiments, the first energy spectrum acquisition module 240 may simulate, using a Monte Carlo algorithm, the monoenergetic depth dose curves corresponding to the multiple monoenergetic energies under the same measurement conditions as those used to acquire the first depth dose curve. A detailed description of the same measurement conditions can be found in the description of step 310 and is not further elaborated here.

[0169] For example, the first energy spectrum acquisition module 240 can acquire corresponding n=100 monoenergetic depth dose curves: E1, E2, ..., E based on n=100 monoenergetic energies: 0.1 MeV, 0.2 MeV, ..., 10 MeV. 100 .

[0170] Some embodiments of this specification determine the monoenergetic energy interval based on energy spectrum accuracy and / or energy spectrum calculation speed, and determine the monoenergetic depth dose curve based on the monoenergetic energy interval. This can ensure the accuracy of the monoenergetic depth dose curve used to represent the first depth dose curve, and can also limit the number of columns of the first and second vectors subsequently represented by the monoenergetic energy based on the energy spectrum calculation speed, thereby improving the computational efficiency of the entire system. A detailed description of representing the first and second vectors by monoenergetic energy can be found in the description of step 350 and will not be repeated here.

[0171] Step 620 : Acquire a first energy spectrum of the first depth dose curve based on the first depth dose curve and the multiple monoenergetic depth dose curves.

[0172] The first energy spectrum may be a description of the monoenergetic energy contained in the first radiation beam and its corresponding first weight.

[0173] In some embodiments, the relationship between the first weight corresponding to the monoenergetic energy and the corresponding monoenergetic depth dose curve and the first depth dose curve corresponding to the first radiation beam can be expressed by formula (8):

[0174]

[0175] Among them, PDD a It can be the first depth dose curve, n can represent n single energy depth dose curves, PDD i It can be the i-th monoenergetic depth dose curve, It may be the first weight of the depth dose (such as the percentage depth dose) corresponding to the i-th monoenergetic depth dose curve in the depth dose (such as the percentage depth dose) corresponding to the first depth dose curve.

[0176] As can be seen from formula (8), based on the first weight corresponding to the monoenergetic energy, the depth doses (such as percentage depth doses) corresponding to the same depth on multiple monoenergetic depth dose curves are weighted and summed to obtain the depth dose (such as percentage depth dose) corresponding to the depth on the first depth dose curve. The first weight corresponding to each monoenergetic energy can represent the contribution of the corresponding monoenergetic ray in the first radiation beam to the depth dose (such as percentage depth dose) of the phantom.

[0177] In some embodiments, the first energy spectrum may be a Monte Carlo spectrum. Accordingly, in some embodiments, the first energy spectrum acquisition module 240 may acquire the first energy spectrum of the first depth dose curve by using, but not limited to, an ant colony optimization algorithm (ACO), a simulated annealing algorithm (SAA), and a Monte Carlo method (MC method).

[0178] For example, the Monte Carlo method is used as an example. The Monte Carlo algorithm is a stochastic simulation method that can obtain the dose distribution of particles in the phantom by stochastically simulating the interaction between particles in the radiation beam and atoms in the phantom.

[0179] In some embodiments, the first energy spectrum acquisition module 240 can simulate multiple monoenergetic depth dose curves based on a Monte Carlo algorithm to obtain a phase space file. The phase space file can be simulated information about the coupled transport of particles of different energies in the radiation beam. In some embodiments, the phase space file can include, but is not limited to, information such as particle energy at a specific plane, particle type, particle position coordinates, particle incidence direction, and particle statistical weight.

[0180] Furthermore, in some embodiments, the first energy spectrum acquisition module 240 may analyze the first depth dose curve based on the phase space file of the device 110 to acquire a first energy spectrum of the first depth dose curve. For example, the first energy spectrum of the first depth dose curve may include first weights corresponding to 100 monoenergetic energies: 0.5% of monoenergetic energy 0.1 MeV, 1% of monoenergetic energy 0.2 MeV, ..., 0.2% of monoenergetic energy 10 MeV.

[0181] Some embodiments of this specification determine the first energy spectrum of the first depth dose curve based on a Monte Carlo algorithm, which can improve the accuracy of the first energy spectrum, thereby improving the accuracy of the second energy spectrum determined based on the first energy spectrum.

[0182] In some embodiments, the first energy spectrum may be a gamma spectrum. Accordingly, in some embodiments, the first energy spectrum acquisition module 240 may simulate multiple monoenergetic depth dose curves based on a convolution algorithm to acquire the first energy spectrum.

[0183] In some embodiments, the first energy spectrum may also be an electron energy spectrum. Accordingly, in some embodiments, the first energy spectrum acquisition module 240 may simulate multiple monoenergetic depth dose curves based on the AXB algorithm to acquire the first energy spectrum.

[0184] The beneficial effects that may be brought about by the embodiments of this specification include but are not limited to: (1) discretizing the first depth dose curve and the second depth dose curve, using the first dose value group and the second dose value group to represent the first depth dose curve and the second depth dose curve respectively, using each numerical value in the relationship data to represent the ratio of the dose value in the corresponding first dose value group and the dose value in the second dose value group, thereby using a simple operation relationship (such as a vector operation relationship) between the first dose value group and the second dose value group to represent the relationship between the first depth dose curve and the second depth dose curve, thereby improving the operation efficiency; (2) based on the relationship between the first depth dose curve and the second depth dose curve that is quickly acquired, the energy spectrum of any other energy level can be quickly acquired based on the known first energy spectrum (such as a Monte Carlo energy spectrum with high accuracy), without relying on the energy spectrum optimization algorithm, thereby improving the adaptability of acquiring energy spectra of other energy levels and saving the time of acquiring energy spectra. It should be noted that different embodiments may have different beneficial effects. In different embodiments, the beneficial effects that may be produced may be any one or a combination of the above, or any other beneficial effects that may be obtained.

[0185] While the basic concepts have been described above, it will be apparent to those skilled in the art that the detailed disclosure is merely illustrative and does not limit this specification. Although not explicitly stated herein, various modifications, improvements, and revisions to this specification may be made by those skilled in the art. Such modifications, improvements, and revisions are suggested in this specification and remain within the spirit and scope of the exemplary embodiments of this specification.

[0186] This specification also uses specific terms to describe the embodiments of this specification. For example, "one embodiment," "an embodiment," and / or "some embodiments" refer to a feature, structure, or characteristic associated with at least one embodiment of this specification. Therefore, it should be emphasized and noted that references to "one embodiment," "an embodiment," or "an alternative embodiment" two or more times in different locations in this specification do not necessarily refer to the same embodiment. Furthermore, certain features, structures, or characteristics of one or more embodiments of this specification may be appropriately combined.

[0187] In addition, unless expressly stated in the claims, the order of the processing elements and sequences, the use of alphanumeric characters, or the use of other names described in this specification are not intended to limit the order of the processes and methods of this specification. Although the above disclosure discusses some of the invention embodiments currently considered useful through various examples, it should be understood that such details are for illustrative purposes only, and the appended claims are not limited to the disclosed embodiments. On the contrary, the claims are intended to cover all modifications and equivalent combinations that are consistent with the spirit and scope of the embodiments of this specification. For example, although the system components described above can be implemented by hardware devices, they can also be implemented only by software solutions, such as installing the described system on an existing server or mobile device.

[0188] Similarly, it should be noted that, in order to simplify the presentation of this specification and thus facilitate understanding of one or more embodiments of the invention, the foregoing descriptions of the embodiments of this specification sometimes combine multiple features into a single embodiment, figure, or description thereof. However, this disclosure method does not imply that the subject matter of this specification requires more features than those recited in the claims. In fact, an embodiment may have fewer features than all of the features of a single disclosed embodiment.

[0189] In some embodiments, numbers are used to describe the quantity of components and attributes. It should be understood that such numbers used in the description of the embodiments are modified by the modifiers "about", "approximately" or "substantially" in some examples. Unless otherwise stated, "about", "approximately" or "substantially" indicate that the numbers are allowed to vary by ±20%. Accordingly, in some embodiments, the numerical parameters used in the description and claims are approximate values, which may change according to the required characteristics of individual embodiments. In some embodiments, the numerical parameters should take into account the specified significant digits and adopt the general method of retaining digits. Although the numerical domains and parameters used to confirm the breadth of their range in some embodiments of this specification are approximate values, in specific embodiments, the settings of such numerical values ​​are as accurate as possible within the feasible range.

[0190] Each patent, patent application, patent application publication, and other materials, such as articles, books, specifications, publications, and documents, cited in this specification is hereby incorporated by reference in its entirety. This excludes any application history that is inconsistent with or conflicts with the content of this specification, as well as any document (currently or subsequently appended to this specification) that limits the broadest scope of the claims of this specification. It should be noted that if the descriptions, definitions, and / or terminology used in the accompanying materials are inconsistent with or conflict with the content of this specification, the descriptions, definitions, and / or terminology used in this specification will control.

[0191] Finally, it should be understood that the embodiments described in this specification are intended only to illustrate the principles of the embodiments of this specification. Other variations may also fall within the scope of this specification. Therefore, by way of example and not limitation, alternative configurations of the embodiments of this specification may be considered consistent with the teachings of this specification. Accordingly, the embodiments of this specification are not limited to the embodiments explicitly described and illustrated in this specification.

Claims

1. A method for obtaining an energy spectrum, characterized in that: The method comprises: Acquire a first depth dose curve corresponding to the first energy and a second depth dose curve corresponding to the second energy; Based on the first depth dose curve and the second depth dose curve, respectively acquiring a first set of dose values ​​and a second set of dose values ​​along a depth direction; determining a relationship between the first depth dose curve and the second depth dose curve based on the first set of dose values ​​and the second set of dose values; Acquire a first energy spectrum corresponding to the first depth dose curve; Based on the first energy spectrum and the relationship between the first depth dose curve and the second depth dose curve, a second energy spectrum corresponding to the second depth dose curve is determined.

2. The method according to claim 1, wherein A relationship between the first depth dose curve and the second depth dose curve is determined based on a ratio of at least one dose value in the first set of dose values ​​to a corresponding at least one dose value in the second set of dose values.

3. The method according to claim 1, wherein The maximum measurement depths corresponding to the first depth dose curve and the second depth dose curve are determined based on the thickness and / or material of the phantom.

4. The method according to claim 3, wherein The acquiring, based on the first depth dose curve and the second depth dose curve, a first set of dose values ​​and a second set of dose values, respectively, comprises: Determine the discrete point interval; determining discrete points based on the maximum measurement depth and the discrete point interval; Based on the discrete points, the first depth dose curve and the second depth dose curve are discretized respectively to obtain the first group of dose values ​​and the second group of dose values ​​corresponding to the discrete points.

5. The method according to claim 1, wherein The first energy is not lower than the second energy.

6. The method according to claim 1, wherein The acquiring, based on the first depth dose curve, a first energy spectrum of the first depth dose curve includes: determining a plurality of monoenergetic depth dose curves based on the first energy; A first energy spectrum of the first depth dose curve is acquired based on the first depth dose curve and the plurality of monoenergetic depth dose curves.

7. The method according to claim 6, wherein The determining of a plurality of monoenergetic depth dose curves based on the first energy includes: determining the monoenergetic energy intervals corresponding to the plurality of monoenergetic depth dose curves; The plurality of monoenergetic depth dose curves are determined based on the first energy and the monoenergetic energy interval.

8. The method according to claim 6, wherein The acquiring, based on the first depth dose curve and the plurality of monoenergetic depth dose curves, a first energy spectrum of the first depth dose curve includes: Simulating the multiple monoenergetic depth dose curves based on a Monte Carlo algorithm to obtain a phase space file; A first energy spectrum of the first depth dose curve is acquired based on the phase space file and the first depth dose curve.

9. The method according to claim 2, wherein: The first set of dose values ​​and the second set of dose values ​​are represented by a first vector and a second vector, respectively. The relationship between the first depth dose curve and the second depth dose curve is represented by a diagonal matrix, and the non-zero elements of each row in the diagonal matrix are determined based on the ratio of the elements in the corresponding rows in the second vector and the first vector.

10. An energy spectrum acquisition system, characterized in that: The system comprises: A curve acquisition module, configured to acquire a first depth dose curve corresponding to a first energy and a second depth dose curve corresponding to a second energy; a dose value acquisition module, configured to acquire a first group of dose values ​​and a second group of dose values ​​along a depth direction based on the first depth dose curve and the second depth dose curve, respectively; a relationship determining module, configured to determine a relationship between the first depth dose curve and the second depth dose curve based on the first set of dose values ​​and the second set of dose values; A first energy spectrum acquisition module, configured to acquire a first energy spectrum corresponding to the first depth dose curve; The second energy spectrum determining module is configured to determine a second energy spectrum of the second depth dose curve based on the first energy spectrum and a relationship between the first depth dose curve and the second depth dose curve.

11. The system according to claim 10, wherein: A relationship between the first depth dose curve and the second depth dose curve is determined based on a ratio of at least one dose value in the first set of dose values ​​to a corresponding at least one dose value in the second set of dose values.

12. The system according to claim 10, wherein: The maximum measurement depths corresponding to the first depth dose curve and the second depth dose curve are determined based on the thickness and / or material of the phantom.

13. The system according to claim 12, wherein: The dosage value acquisition module is further used for: Determine the discrete point interval; determining discrete points based on the maximum measurement depth and the discrete point interval; Based on the discrete points, the first depth dose curve and the second depth dose curve are discretized respectively to obtain the first group of dose values ​​and the second group of dose values ​​corresponding to the discrete points.

14. The system according to claim 10, wherein: The first energy is not lower than the second energy.

15. The system according to claim 10, wherein: The first energy spectrum acquisition module is further configured to: determining a plurality of monoenergetic depth dose curves based on the first energy; A first energy spectrum of the first depth dose curve is acquired based on the first depth dose curve and the plurality of monoenergetic depth dose curves.

16. The system according to claim 15, wherein: The first energy spectrum acquisition module is further configured to: determining the monoenergetic energy intervals corresponding to the plurality of monoenergetic depth dose curves; The plurality of monoenergetic depth dose curves are determined based on the first energy and the monoenergetic energy interval.

17. The system according to claim 15, wherein: The first energy spectrum acquisition module is further configured to: Simulating the multiple monoenergetic depth dose curves based on a Monte Carlo algorithm to obtain a phase space file; The first depth dose curve is analyzed based on the phase space file to obtain a first energy spectrum of the first depth dose curve.

18. The system of claim 11, wherein: The first set of dose values ​​and the second set of dose values ​​are represented by a first vector and a second vector, respectively. The relationship between the first depth dose curve and the second depth dose curve is represented by a diagonal matrix, and the non-zero elements of each row in the diagonal matrix are determined based on the ratio of the elements in the corresponding rows in the second vector and the first vector.

19. A computer-readable storage medium storing computer instructions, wherein when a computer reads the computer instructions in the storage medium, the computer executes the energy spectrum acquisition method according to any one of claims 1 to 9.

20. An energy spectrum acquisition device, characterized in that: The device includes at least one processor and at least one memory; the at least one memory is used to store computer instructions; and the at least one processor is used to execute at least part of the computer instructions to implement the method according to claims 1 to 9.

Images