Method, device, equipment and medium for determining energy calibration of radioactive detector

Abstract

The present disclosure provides a method, apparatus, device and medium for determining the energy scale of a radioactive detector. The method includes: obtaining background energy spectra of multiple known nuclides, wherein the background energy spectra are obtained by detecting the multiple known nuclides with the radioactive detector; performing differential processing on the background energy spectra to obtain a channel address-difference map, wherein the channel address-difference map includes multiple effective peak information determined from multiple full-energy peak information; determining a target coefficient based on the channel address-difference map, wherein the target coefficient represents a scale coefficient in an initial energy channel address function; and determining a target energy channel address function based on the channel address-difference map and the initial energy channel address function, wherein the target energy channel address function represents the energy scale of the radioactive detector.

Description

Technical Field

[0001] The present disclosure relates to the field of environmental monitoring, and in particular to a method, apparatus, device, medium and program product for determining the energy calibration of a radioactive detector. Background Art

[0002] Radiation environmental monitoring is an important component of environmental monitoring and the foundation of radiation environmental management. The prevention and control of environmental radioactive contamination is also a hot topic of public concern. To ensure the accuracy of detection results, the detector system must be regularly energy calibrated.

[0003] Energy calibration typically involves detecting and measuring gamma-emitting nuclides in environmental samples using common detectors. This provides information about the sample's gamma energy spectrum, specifically the relationship between photon counts and channel addresses. This information is then used to calibrate the detector based on the quadratic nonlinear relationship between energy and channel addresses. However, due to changes in ambient conditions and the system's own influence, regular energy calibration of the detector system is necessary to ensure the accuracy of the detector's detection results. Summary of the Invention

[0004] In view of the above problems, the present disclosure provides a method, apparatus, device, medium and program product for determining the energy calibration of a radioactive detector, which improve the energy calibration correction speed of the radioactive detector.

[0005] According to a first aspect of the present disclosure, a method for determining an energy calibration of a radioactive detector is provided, comprising:

[0006] Acquiring background energy spectra of a plurality of known nuclides, wherein the background energy spectra are obtained by detecting the plurality of known nuclides by the radioactivity detector;

[0007] Performing differential processing on the background energy spectrum to obtain a channel address-difference map, wherein the channel address-difference map includes a plurality of effective peak information determined from a plurality of full energy peak information;

[0008] determining a target coefficient according to the track address-difference map, wherein the target coefficient represents a scale coefficient in the initial energy track address function;

[0009] A target energy track address function is determined based on the track address-difference map and the initial energy track address function, wherein the target energy track address function represents the energy scale of the radioactive detector.

[0010] According to an embodiment of the present disclosure, determining a target coefficient according to the track address-difference map includes:

[0011] Determining a plurality of target intervals of a target nuclide from the track address-difference map, wherein the target intervals include a plurality of valid peak information of the target nuclide;

[0012] The target coefficient is determined based on the multiple target intervals.

[0013] According to an embodiment of the present disclosure, determining the target coefficient according to the plurality of target intervals includes:

[0014] For each target interval, an initial weight is obtained according to a plurality of full energy peak weighting coefficients and a plurality of differential coefficients of full energy peak information;

[0015] Based on the preset screening rules, a target coefficient is determined according to the multiple initial weights.

[0016] According to an embodiment of the present disclosure, obtaining an initial weight according to the plurality of full energy peak weighting coefficients and the differential coefficients of the plurality of full energy peak information includes:

[0017] For each effective peak information, determining a product value according to the full energy peak weighting coefficient corresponding to the effective peak information and the differential coefficient at the full energy peak energy channel address in the effective peak information;

[0018] A plurality of the product values ​​are summed to obtain the initial weight of the target interval.

[0019] According to an embodiment of the present disclosure, the differential processing includes second-order differential processing, and the differential coefficient includes a second-order differential coefficient.

[0020] According to an embodiment of the present disclosure, determining the target coefficient based on the preset screening rules and the multiple initial weights includes:

[0021] Drawing a weight-coefficient relationship diagram based on the multiple initial weights, wherein the horizontal axis of the weight-coefficient relationship diagram is the coefficient scale and the vertical axis is the weight scale;

[0022] In the weight-coefficient relationship diagram, the horizontal coordinate corresponding to the maximum weight value is determined as the target coefficient.

[0023] According to an embodiment of the present disclosure, determining a target energy track address function according to the track address-difference map and the initial energy track address function includes:

[0024] Constructing a plurality of coordinate points in a channel address-energy coordinate system according to the channel address corresponding to each of the total energy peak information in each of the target intervals and the known total energy peak energy of the target nuclide;

[0025] Based on the target coefficient, the initial energy address function is solved according to the plurality of coordinate points to obtain the target energy address function.

[0026] According to an embodiment of the present disclosure, solving the initial energy address function according to the plurality of coordinate points to obtain the target energy address function includes:

[0027] Performing fitting processing on the plurality of coordinate points in the track address-energy coordinate system to obtain a target curve;

[0028] Solving the initial energy address function using the target curve to obtain a constant term coefficient and a quadratic term coefficient, wherein the target coefficient represents the linear term coefficient;

[0029] Based on the initial energy address function, the target energy address function is obtained according to the constant term coefficient, the target coefficient and the quadratic term coefficient.

[0030] According to an embodiment of the present disclosure, the initial weight is as shown in the first formula, and the initial energy address function is as shown in the second formula;

[0031]

[0032] E=c0+c1×ch+c2×ch 2

[0033] Among them, W j is the initial weight of the jth target interval; w i，j is the weighting coefficient of the i-th full energy peak in the j-th target interval; sec i，j is the differential coefficient at the energy address of the ith full-energy peak in the jth target interval, E is the energy of the γ-ray full-energy peak, ch is the address of the full-energy peak of the ray, c0, c1, and c2 are the constant term coefficient, target coefficient, and quadratic term coefficient to be solved respectively.

[0034] According to a second aspect of the present disclosure, there is provided an energy calibration determination device for a radioactive detector, comprising:

[0035] an acquisition module, configured to acquire background energy spectra of a plurality of known nuclides, wherein the background energy spectra are obtained by detecting the plurality of known nuclides by the radioactive detector;

[0036] a differential module, configured to perform differential processing on the background energy spectrum to obtain a channel address-difference map, wherein the channel address-difference map includes a plurality of effective peak information determined from a plurality of full energy peak information;

[0037] a first determining module, configured to determine a target coefficient based on the track address-difference map, wherein the target coefficient represents a scale coefficient in an initial energy track address function;

[0038] The second determination module is configured to determine a target energy track address function according to the track address-difference map and the initial energy track address function, wherein the target energy track address function represents the energy scale of the radioactive detector.

[0039] The third aspect of the present disclosure provides an electronic device comprising: one or more processors; and a memory for storing one or more programs, wherein when the one or more programs are executed by the one or more processors, the one or more processors execute the above-mentioned method for determining the energy scale of the radioactive detector.

[0040] The fourth aspect of the present disclosure further provides a computer-readable storage medium having executable instructions stored thereon, which, when executed by a processor, causes the processor to execute the above-mentioned method for determining the energy scale of a radioactive detector.

[0041] The fifth aspect of the present disclosure further provides a computer program product, comprising a computer program, which implements the above-mentioned energy calibration determination method of the radioactive detector when executed by a processor.

[0042] According to the embodiments of the present disclosure, the background energy spectrum obtained by detecting a plurality of the known nuclides through the radioactive detector to be calibrated is differentially processed to obtain a channel address-difference diagram. Based on the channel address-difference diagram, a target coefficient in the initial energy channel address function can be determined. By combining the channel address-difference diagram and the initial energy channel address function, the target energy channel address function is determined. The target energy channel address function can then be used to perform energy calibration correction processing on the radioactive detector to be calibrated, thereby avoiding the problem of inaccurate detection results of the calibrated radioactive detector caused by the manual use of a radioactive source for energy calibration in related technologies, thereby improving the correction efficiency and detection accuracy of the radioactive detector. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] The above contents and other objects, features and advantages of the present disclosure will become more apparent through the following description of the embodiments of the present disclosure with reference to the accompanying drawings, in which:

[0044] Figure 1 A flowchart of a method for determining the energy calibration of a radioactive detector according to an embodiment of the present disclosure is schematically shown;

[0045] Figure 2 A schematic diagram of a background energy spectrum according to an embodiment of the present disclosure is schematically shown;

[0046] Figure 3Schematically shows a track address-difference map according to an embodiment of the present disclosure;

[0047] Figure 4 Schematically shows a weight-coefficient relationship diagram according to an embodiment of the present disclosure;

[0048] Figure 5 Schematically shows a schematic diagram of the channel address-energy coordinate system according to an embodiment of the present disclosure;

[0049] Figure 6 A block diagram schematically illustrates a structure of an energy calibration determination device for a radioactive detector according to an embodiment of the present disclosure;

[0050] Figure 7 A block diagram of an electronic device suitable for implementing a method for determining the energy calibration of a radioactive detector according to an embodiment of the present disclosure is schematically shown. DETAILED DESCRIPTION

[0051] Hereinafter, embodiments of the present disclosure will be described with reference to the accompanying drawings. However, it should be understood that these descriptions are merely exemplary and are not intended to limit the scope of the present disclosure. In the detailed description below, for ease of explanation, many specific details are set forth to provide a comprehensive understanding of the embodiments of the present disclosure. However, it is apparent that one or more embodiments may also be implemented without these specific details. In addition, in the following description, descriptions of well-known structures and technologies are omitted to avoid unnecessary confusion of the concepts of the present disclosure.

[0052] The terms used herein are only for describing specific embodiments and are not intended to limit the present disclosure. The terms "comprise," "include," etc. used herein indicate the presence of the features, steps, operations, and / or components, but do not exclude the presence or addition of one or more other features, steps, operations, or components.

[0053] All terms used herein (including technical and scientific terms) have the meanings commonly understood by those skilled in the art unless otherwise defined. It should be noted that the terms used herein should be interpreted as having a meaning consistent with the context of this specification and should not be interpreted in an idealized or overly rigid manner.

[0054] When expressions such as "at least one of A, B and C, etc." are used, they should generally be interpreted in accordance with the meaning of the expression commonly understood by those skilled in the art (for example, "a system having at least one of A, B and C" should include but is not limited to a system having A alone, B alone, C alone, A and B, A and C, B and C, and / or A, B, C, etc.).

[0055] Embodiments of the present disclosure provide a method, apparatus, device and medium for determining the energy scale of a radioactive detector, the method comprising: acquiring background energy spectra of a plurality of known nuclides, wherein the background energy spectra are obtained by detecting the plurality of known nuclides by the radioactive detector; performing differential processing on the background energy spectra to obtain a channel address-difference map, wherein the channel address-difference map includes a plurality of effective peak information determined from a plurality of full-energy peak information; determining a target coefficient based on the channel address-difference map, wherein the target coefficient represents a scale coefficient in an initial energy channel address function; determining a target energy channel address function based on the channel address-difference map and the initial energy channel address function, wherein the target energy channel address function represents the energy scale of the radioactive detector.

[0056] Figure 1 A flow chart of a method for determining the energy calibration of a radioactive detector according to an embodiment of the present disclosure is schematically shown. Figure 2 A schematic diagram of a background energy spectrum according to an embodiment of the present disclosure is schematically shown. Figure 3 The figure schematically shows a track address-difference map according to an embodiment of the present disclosure.

[0057] like Figure 1 As shown, the method for determining the energy calibration of a radioactive detector includes operations S110 to S140.

[0058] In operation S110, background energy spectra of a plurality of known nuclides are acquired, wherein the background energy spectra are obtained by detecting the plurality of known nuclides using a radioactive detector;

[0059] In operation S120, a differential process is performed on the background energy spectrum to obtain a channel address-difference map, wherein the channel address-difference map includes a plurality of effective peak information determined from the plurality of full energy peak information;

[0060] In operation S130, a target coefficient is determined according to the track address-difference map, wherein the target coefficient represents a scale coefficient in the initial energy track address function.

[0061] In operation S140, a target energy track address function is determined based on the track address-difference map and the initial energy track address function, wherein the target energy track address function represents the energy scale of the radioactive detector.

[0062] According to embodiments of the present disclosure, energy calibration refers to the experimental work of determining the correspondence between the energy of incident particles in a radioactive detector and the number of channels in a multichannel analyzer. Radioactive detectors are typically installed outdoors in the environment to be detected. After a period of operation, the energy calibration of the radioactive detector must be recalibrated to ensure detection accuracy.

[0063] According to the embodiments of the present disclosure, a channel address refers to an energy range, and its energy index is the channel address. Full-energy peak information refers to the spectral peak formed by the detector output pulse when all the energy of the incident radiation is lost within the detector's sensitive volume. Effective peak information refers to the full-energy peak information of a known nuclide determined from multiple full-energy peak information of multiple known nuclides. Known nuclides may include thorium (Th-232), radium (Ra-226), potassium (K-40), etc.

[0064] According to the embodiment of the present disclosure, the radioactive detector is used to detect multiple known nuclides, and the following can be obtained: Figure 2 As shown in FIG, the background energy spectrum of multiple known nuclides is subjected to differential processing, for example, using the Mariscotti method for second-order differential processing, to obtain the following: Figure 3 The channel address-difference diagram shown clearly shows the effective peak information of multiple known nuclides. Based on the channel address-difference diagram, the target coefficient can be determined. The target coefficient is the second-order coefficient of the initial energy channel address function. Then, based on the target coefficient, the channel address-difference diagram and the initial energy channel address function, the target energy channel address function is determined. Based on the target energy channel address function, the energy scale of the radioactive detector can be calibrated.

[0065] According to the embodiments of the present disclosure, the background energy spectrum obtained by detecting multiple known nuclides through the radioactive detector to be calibrated is differentially processed to obtain a channel address-difference diagram. Based on the channel address-difference diagram, a target coefficient in the initial energy channel address function can be determined. By combining the channel address-difference diagram and the initial energy channel address function, the target energy channel address function is determined, so that the target energy channel address function can be used to perform energy calibration correction processing on the radioactive detector to be calibrated, avoiding the problem of inaccurate detection results of the calibrated radioactive detector caused by the manual use of a radioactive source for energy calibration in related technologies, thereby improving the correction efficiency and detection accuracy of the radioactive detector.

[0066] According to an embodiment of the present disclosure, determining a target coefficient based on a track address-difference map includes the following operations:

[0067] Determining a plurality of target intervals of the target nuclide from the track address-difference map, wherein the target intervals include a plurality of valid peak information of the target nuclide;

[0068] A target coefficient is determined based on multiple target intervals.

[0069] According to an embodiment of the present disclosure, the target interval is determined based on the effective peak information of a known nuclide before calibration. For example, the target intervals of thorium are known to be [c1, c2], [c3, c4] and [c5, c6], where c iIt is the horizontal axis scale of the channel address-difference diagram. The specific value is related to the specific full energy peak information of the target nuclide.

[0070] According to an embodiment of the present disclosure, after a plurality of target intervals are determined, target coefficients are determined based on the track address-difference maps corresponding to the plurality of target intervals.

[0071] According to an embodiment of the present disclosure, determining a target coefficient based on multiple target intervals includes the following operations:

[0072] For each target interval, an initial weight is obtained based on multiple full energy peak weighting coefficients and multiple full energy peak information difference coefficients;

[0073] Based on the preset screening rules and multiple initial weights, the target coefficient is determined.

[0074] According to an embodiment of the present disclosure, for each target interval, an initial weight of the target interval is obtained based on multiple effective peak information within the target interval, such as the full-energy peak weighting coefficient corresponding to each effective peak information and the differential coefficient at the effective peak information during differential processing. Finally, the target coefficient is determined from the multiple initial weights using preset screening rules.

[0075] It should be noted that the full energy peak weighting coefficient is the weighting coefficient of the full energy peak energy, and the weighting coefficient is specifically set by the staff according to actual needs.

[0076] According to an embodiment of the present disclosure, obtaining an initial weight according to a plurality of full energy peak weighting coefficients and a plurality of differential coefficients of full energy peak information includes the following operations:

[0077] For each effective peak information, a product value is determined according to the full energy peak weighting coefficient corresponding to the effective peak information and the differential coefficient at the full energy peak energy channel address in the effective peak information;

[0078] The multiple product values ​​are summed up to obtain the initial weight of the target interval.

[0079] According to an embodiment of the present disclosure, when calculating the initial weight of each target interval, for each valid peak information within the target interval, the full-energy peak weighting coefficient and the differential coefficient corresponding to the valid peak information are multiplied to obtain a product value, and the multiple product values ​​corresponding to multiple valid peak information are summed, and the sum result is used as the initial weight of the target interval.

[0080] According to an embodiment of the present disclosure, the differential processing includes a second-order differential processing, and the differential coefficient includes a second-order differential coefficient.

[0081] It should be noted that the second-order difference processing method used in the present disclosure is the Mariscotti method, but the present disclosure is not limited to using only this method.

[0082] Figure 4 The figure schematically shows a weight-coefficient relationship diagram according to an embodiment of the present disclosure.

[0083] According to an embodiment of the present disclosure, based on a preset screening rule and according to a plurality of initial weights, determining a target coefficient includes the following operations:

[0084] Based on the multiple initial weights, a weight-coefficient relationship diagram is drawn, wherein the horizontal axis of the weight-coefficient relationship diagram is the coefficient scale and the vertical axis is the weight scale;

[0085] In the weight-coefficient relationship diagram, the horizontal coordinate corresponding to the maximum weight is determined as the target coefficient.

[0086] According to an embodiment of the present disclosure, after determining the initial weight of each target interval, a weight-coefficient relationship diagram is drawn based on multiple initial weights, such as Figure 4 As shown, each initial weight is plotted in the weight-coefficient relationship diagram, and the points are connected. Finally, the maximum weight value is determined according to the plotted diagram, and the horizontal coordinate corresponding to the maximum weight value is determined as the target coefficient, for example Figure 4 The horizontal coordinate X=0.185 corresponding to Y=26.0295 is determined as the target coefficient, and then the target energy address function is determined according to the address-difference diagram and the initial energy address function.

[0087] Figure 5 The figure schematically shows a schematic diagram of the track address-energy coordinate system according to an embodiment of the present disclosure.

[0088] According to an embodiment of the present disclosure, determining a target energy-address function according to a track-difference map and an initial energy-address function includes the following operations:

[0089] According to the channel address corresponding to each full-energy peak information in each target interval and the known full-energy peak energy of the target nuclide, multiple coordinate points are constructed in the channel address-energy coordinate system;

[0090] Based on the target coefficient, the initial energy address function is solved according to multiple coordinate points to obtain the target energy address function.

[0091] According to the embodiments of the present disclosure, since the known total energy peak energy of the target nuclide is a value determined before the correction begins, combined with the channel address corresponding to each total energy peak information in the target interval, the energy coordinates of each total energy peak can be marked in the channel address-energy coordinate system, thereby obtaining multiple coordinate points about the total energy peak information.

[0092] According to the embodiment of the present disclosure, since the target coefficient has been determined, the other coefficients in the initial energy address function can be obtained according to the plotted Figure 5 The multiple coordinate points of the address-energy coordinate system shown in FIG are solved to obtain the final target energy address function, as shown in FIG. Figure 5 The formula shown in .

[0093] According to an embodiment of the present disclosure, solving the initial energy address function according to multiple coordinate points to obtain the target energy address function includes the following operations:

[0094] Fitting multiple coordinate points in the track address-energy coordinate system to obtain the target curve;

[0095] The target curve is used to solve the initial energy address function to obtain the constant term coefficient and the quadratic term coefficient, where the target coefficient represents the linear term coefficient;

[0096] Based on the initial energy address function, the target energy address function is obtained according to the constant term coefficient, the target coefficient and the quadratic term coefficient.

[0097] According to an embodiment of the present disclosure, a plurality of coordinate points are fitted in the channel address-energy coordinate system to obtain a target curve, wherein the purpose of fitting is to make the final target energy channel address function as accurate as possible to ensure the calibration accuracy of the final calibrated radioactive detector when it is used.

[0098] According to an embodiment of the present disclosure, after obtaining the target curve, the target curve can be used to solve the constant term coefficient and the quadratic term coefficient in the initial energy address function, and finally the final target energy address function is constructed based on the determined target coefficient, constant term coefficient and quadratic term coefficient.

[0099] According to an embodiment of the present disclosure, the initial weight is shown in formula (1), and the initial energy address function is shown in formula (2);

[0100]

[0101] E=c0+c1×ch+c2×ch 2 (2)

[0102] Among them, W j is the initial weight of the jth target interval; w i，j is the weighting coefficient of the i-th full energy peak in the j-th target interval; sec i，j is the differential coefficient at the energy address of the ith full-energy peak in the jth target interval, E is the energy of the γ-ray full-energy peak, ch is the address of the full-energy peak of the ray, c0, c1, and c2 are the constant term coefficient, target coefficient, and quadratic term coefficient to be solved respectively.

[0103] In some embodiments, under normal circumstances, the quadratic term coefficient c2 in the energy-channel address scaling function approaches 0, and the value of the constant term coefficient c0 is also near 0, and hardly changes with changes in external conditions. Therefore, in the process of solving the target coefficient of the present disclosure, the energy scaling function can be approximately written as formula (3).

[0104] E=c1×ch (3)

[0105] According to an embodiment of the present disclosure, after determining the target coefficient, the constant term coefficient and the quadratic term coefficient in formula (2) are solved using the channel address-energy coordinate system, and finally the target energy channel address function can be obtained.

[0106] Figure 6 The structural block diagram of the energy calibration determination device of the radioactive detector according to the embodiment of the present disclosure is schematically shown.

[0107] like Figure 6 As shown, the energy calibration determination device 600 of the radioactive detector includes an acquisition module 610 , a difference module 620 , a first determination module 630 and a second determination module 640 .

[0108] The acquisition module 610 is used to acquire background energy spectra of a plurality of known nuclides, wherein the background energy spectra are obtained by detecting the plurality of known nuclides using a radioactive detector.

[0109] The difference module 620 is used to perform difference processing on the background energy spectrum to obtain a channel address-difference map, wherein the channel address-difference map includes multiple effective peak information determined from multiple full-energy peak information.

[0110] The first determination module 630 is configured to determine a target coefficient according to the track address-difference map, wherein the target coefficient represents a scale coefficient in the initial energy track address function.

[0111] The second determination module 640 is configured to determine a target energy track address function according to the track address-difference map and the initial energy track address function, wherein the target energy track address function represents the energy scale of the radioactive detector.

[0112] According to the embodiments of the present disclosure, the background energy spectrum obtained by detecting a plurality of the known nuclides through the radioactive detector to be calibrated is differentially processed to obtain a channel address-difference diagram. Based on the channel address-difference diagram, a target coefficient in the initial energy channel address function can be determined. By combining the channel address-difference diagram and the initial energy channel address function, the target energy channel address function is determined. The target energy channel address function can then be used to perform energy calibration correction processing on the radioactive detector to be calibrated, thereby avoiding the problem of inaccurate detection results of the calibrated radioactive detector caused by the manual use of a radioactive source for energy calibration in related technologies, thereby improving the correction efficiency and detection accuracy of the radioactive detector.

[0113] According to an embodiment of the present disclosure, the first determining module 630 includes a first determining submodule and a second determining submodule.

[0114] The first determination submodule is used to determine multiple target intervals of the target nuclide from the channel address-difference map, wherein the target interval includes multiple valid peak information of the target nuclide.

[0115] The second determination submodule is used to determine the target coefficient according to multiple target intervals.

[0116] According to an embodiment of the present disclosure, the second determining submodule includes a first obtaining unit and a determining unit.

[0117] The first obtaining unit is used to obtain an initial weight for each target interval according to multiple full energy peak weighting coefficients and multiple full energy peak information difference coefficients.

[0118] The determination unit is used to determine the target coefficient based on a preset screening rule and a plurality of initial weights.

[0119] According to an embodiment of the present disclosure, the first obtaining unit includes a first determining subunit and an obtaining subunit.

[0120] The first determining subunit is used to determine a product value for each effective peak information according to the full energy peak weighting coefficient corresponding to the effective peak information and the differential coefficient at the full energy peak energy channel address in the effective peak information.

[0121] A subunit is obtained, which is used to sum a plurality of product values ​​to obtain an initial weight of a target interval.

[0122] According to an embodiment of the present disclosure, the differential processing includes a second-order differential processing, and the differential coefficient includes a second-order differential coefficient.

[0123] According to an embodiment of the present disclosure, the determining unit includes a drawing subunit and a second determining subunit.

[0124] The drawing subunit is used to draw a weight-coefficient relationship diagram based on multiple initial weights, wherein the horizontal axis in the weight-coefficient relationship diagram is the coefficient scale and the vertical axis is the weight scale.

[0125] The second determining subunit is configured to determine the horizontal coordinate corresponding to the maximum weight value in the weight-coefficient relationship diagram as the target coefficient.

[0126] According to an embodiment of the present disclosure, the second determining module 640 includes a constructing submodule and a solving submodule.

[0127] The construction submodule is used to construct multiple coordinate points in the channel address-energy coordinate system according to the channel address corresponding to each full energy peak information in each target interval and the known full energy peak energy of the target nuclide.

[0128] The solving submodule is used to solve the initial energy address function based on the target coefficient and multiple coordinate points to obtain the target energy address function.

[0129] According to an embodiment of the present disclosure, the solution submodule includes a fitting unit, a solution unit and a second obtaining unit.

[0130] The fitting unit is used to perform fitting processing on multiple coordinate points in the channel address-energy coordinate system to obtain a target curve.

[0131] The solving unit is used to solve the initial energy address function using the target curve to obtain the constant term coefficient and the quadratic term coefficient, wherein the target coefficient represents the linear term coefficient.

[0132] The second obtaining unit is used to obtain the target energy address function based on the initial energy address function according to the constant term coefficient, the target coefficient and the quadratic term coefficient.

[0133] According to an embodiment of the present disclosure, any multiple modules among the acquisition module 610, the difference module 620, the first determination module 630, and the second determination module 640 can be combined into one module for implementation, or any one of the modules can be split into multiple modules. Alternatively, at least part of the functions of one or more of these modules can be combined with at least part of the functions of other modules and implemented in one module. According to an embodiment of the present disclosure, at least one of the acquisition module 610, the difference module 620, the first determination module 630, and the second determination module 640 can be at least partially implemented as a hardware circuit, such as a field programmable gate array (FPGA), a programmable logic array (PLA), a system on a chip, a system on a substrate, a system on a package, an application specific integrated circuit (ASIC), or can be implemented by hardware or firmware such as any other reasonable way of integrating or packaging the circuit, or implemented in any one of the three implementation methods of software, hardware, and firmware, or in an appropriate combination of any of them. Alternatively, at least one of the acquisition module 610 , the difference module 620 , the first determination module 630 , and the second determination module 640 may be at least partially implemented as a computer program module, which may perform corresponding functions when executed.

[0134] Figure 7 A block diagram of an electronic device suitable for implementing a method for determining the energy calibration of a radioactive detector according to an embodiment of the present disclosure is schematically shown.

[0135] like Figure 7 As shown, the electronic device 700 according to an embodiment of the present disclosure includes a processor 701, which can perform various appropriate actions and processes according to a program stored in a read-only memory (ROM) 702 or a program loaded from a storage part 708 into a random access memory (RAM) 703. The processor 701 may, for example, include a general-purpose microprocessor (e.g., a CPU), an instruction set processor and / or a related chipset and / or a special-purpose microprocessor (e.g., an application-specific integrated circuit (ASIC)), etc. The processor 701 may also include an onboard memory for caching purposes. The processor 701 may include a single processing unit or multiple processing units for performing different actions of the method flow according to an embodiment of the present disclosure.

[0136] Various programs and data required for the operation of the electronic device 700 are stored in the RAM 703. The processor 701, ROM 702, and RAM 703 are connected to each other via a bus 704. The processor 701 executes the various operations of the method flow according to the embodiment of the present disclosure by executing the programs in the ROM 702 and / or RAM 703. It should be noted that the programs may also be stored in one or more memories other than the ROM 702 and RAM 703. The processor 701 may also execute the various operations of the method flow according to the embodiment of the present disclosure by executing the programs stored in the one or more memories.

[0137] According to an embodiment of the present disclosure, the electronic device 700 may further include an input / output (I / O) interface 705, which is also connected to the bus 704. The electronic device 700 may further include one or more of the following components connected to the I / O interface 705: an input portion 706 including a keyboard, a mouse, etc.; an output portion 707 including devices such as a cathode ray tube (CRT), a liquid crystal display (LCD), and a speaker; a storage portion 708 including a hard disk; and a communication portion 709 including a network interface card such as a LAN card or a modem. The communication portion 709 performs communication processing via a network such as the Internet. A drive 710 is also connected to the I / O interface 705 as needed. A removable medium 711, such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, etc., is installed in the drive 710 as needed, so that a computer program read therefrom can be installed into the storage portion 708 as needed.

[0138] The present disclosure also provides a computer-readable storage medium, which may be included in the device / apparatus / system described in the above embodiments, or may exist independently and not be incorporated into the device / apparatus / system. The computer-readable storage medium carries one or more programs, and when executed, implements the method according to the embodiments of the present disclosure.

[0139] According to an embodiment of the present disclosure, a computer-readable storage medium may be a non-volatile computer-readable storage medium, for example, it may include but is not limited to: a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination thereof. In the present disclosure, a computer-readable storage medium may be any tangible medium containing or storing a program that can be used by or in combination with an instruction execution system, apparatus, or device. For example, according to an embodiment of the present disclosure, a computer-readable storage medium may include the ROM 702 and / or RAM 703 described above and / or one or more memories other than ROM 702 and RAM 703.

[0140] The embodiments of the present disclosure also include a computer program product, which includes a computer program containing program code for executing the method shown in the flowchart. When the computer program product is executed in a computer system, the program code is used to cause the computer system to implement the item recommendation method provided by the embodiments of the present disclosure.

[0141] The computer program executes the above functions defined in the system / device of the embodiment of the present disclosure when the processor 701 executes the computer program. According to the embodiment of the present disclosure, the system, device, module, unit, etc. described above can be implemented by a computer program module.

[0142] In one embodiment, the computer program may be stored on a tangible storage medium such as an optical storage device or a magnetic storage device. In another embodiment, the computer program may be transmitted and distributed in the form of a signal on a network medium, downloaded and installed via the communication portion 709, and / or installed from a removable medium 711. The program code contained in the computer program may be transmitted using any appropriate network medium, including but not limited to wireless, wired, or any suitable combination thereof.

[0143] In such an embodiment, the computer program can be downloaded and installed from a network via the communication section 709, and / or installed from a removable medium 711. When the computer program is executed by the processor 701, the above-described functions defined in the system of the embodiment of the present disclosure are performed. According to the embodiment of the present disclosure, the systems, devices, means, modules, units, etc. described above can be implemented by computer program modules.

[0144] According to an embodiment of the present disclosure, the program code for executing the computer program provided by the embodiment of the present disclosure can be written in any combination of one or more programming languages. Specifically, these computer programs can be implemented using high-level procedural and / or object-oriented programming languages, and / or assembly / machine languages. Programming languages ​​include, but are not limited to, languages ​​such as Java, C++, python, "C" or similar programming languages. The program code can be executed entirely on the user computing device, partially on the user device, partially on a remote computing device, or entirely on a remote computing device or server. In cases involving a remote computing device, the remote computing device can be connected to the user computing device through any type of network, including a local area network (LAN) or a wide area network (WAN), or can be connected to an external computing device (for example, using an Internet service provider to connect via the Internet).

[0145] The flowcharts and block diagrams in the accompanying drawings illustrate the possible implementation architecture, functions and operations of the systems, methods and computer program products according to various embodiments of the present disclosure. In this regard, each box in the flowchart or block diagram can represent a module, program segment, or a part of code, and the above-mentioned module, program segment, or a part of code contains one or more executable instructions for implementing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the box can also occur in an order different from that marked in the accompanying drawings. For example, two boxes represented in succession can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram or flowchart, and the combination of boxes in the block diagram or flowchart, can be implemented with a dedicated hardware-based system that performs the specified function or operation, or can be implemented with a combination of dedicated hardware and computer instructions.

[0146] Those skilled in the art will appreciate that the features described in the various embodiments and / or claims of this disclosure may be combined and / or coupled in various ways, even if such combinations and / or couplings are not explicitly described in this disclosure. In particular, the features described in the various embodiments and / or claims of this disclosure may be combined and / or coupled in various ways without departing from the spirit and teachings of this disclosure. All such combinations and / or couplings are intended to fall within the scope of this disclosure.

[0147] The embodiments of the present disclosure are described above. However, these embodiments are for illustrative purposes only and are not intended to limit the scope of the present disclosure. Although each embodiment has been described separately above, this does not mean that the measures in each embodiment cannot be used in combination to advantage. The scope of the present disclosure is defined by the appended claims and their equivalents. Without departing from the scope of the present disclosure, those skilled in the art may make various substitutions and modifications, which should all fall within the scope of the present disclosure.

Claims

1. A method for determining the energy calibration of a radioactive detector, comprising: Acquiring background energy spectra of a plurality of known nuclides, wherein the background energy spectra are obtained by detecting the plurality of known nuclides by the radioactivity detector; Performing differential processing on the background energy spectrum to obtain a channel address-difference map, wherein the channel address-difference map includes a plurality of effective peak information determined from a plurality of full energy peak information; determining a target coefficient according to the track address-difference map, wherein the target coefficient represents a scale coefficient in the initial energy track address function; Determine a target energy track address function according to the track address-difference map and the initial energy track address function, wherein the target energy track address function represents the energy scale of the radioactive detector. The determining of the target coefficient according to the track address-difference map includes: determining a plurality of target intervals of the target nuclide from the track address-difference map, wherein the target intervals include a plurality of effective peak information of the target nuclide; determining the target coefficient according to the plurality of target intervals, The method of determining the target energy channel address function based on the channel address-difference map and the initial energy channel address function includes: constructing a plurality of coordinate points in a channel address-energy coordinate system based on the channel address corresponding to each of the total energy peak information in each of the target intervals and the known total energy peak energy of the target nuclide; and solving the initial energy channel address function based on the plurality of coordinate points based on the target coefficient to obtain the target energy channel address function.

2. The method according to claim 1, wherein The determining the target coefficient according to the plurality of target intervals includes: For each target interval, an initial weight is obtained according to a plurality of full energy peak weighting coefficients and a plurality of differential coefficients of full energy peak information; Based on the preset screening rules, a target coefficient is determined according to the multiple initial weights.

3. The method according to claim 2, wherein: The step of obtaining an initial weight according to the plurality of full energy peak weighting coefficients and the differential coefficients of the plurality of full energy peak information comprises: For each effective peak information, determining a product value according to the full energy peak weighting coefficient corresponding to the effective peak information and the differential coefficient at the full energy peak energy channel address in the effective peak information; A plurality of the product values ​​are summed to obtain the initial weight of the target interval.

4. The method according to claim 3, wherein: The differential processing includes second-order differential processing, and the differential coefficients include second-order differential coefficients.

5. The method according to claim 2, wherein: The step of determining a target coefficient based on a plurality of the initial weights based on a preset screening rule includes: Drawing a weight-coefficient relationship diagram based on the multiple initial weights, wherein the horizontal axis of the weight-coefficient relationship diagram is the coefficient scale and the vertical axis is the weight scale; In the weight-coefficient relationship diagram, the horizontal coordinate corresponding to the maximum weight value is determined as the target coefficient.

6. The method according to claim 1, wherein Solving the initial energy address function according to the plurality of coordinate points to obtain the target energy address function includes: Performing fitting processing on the plurality of coordinate points in the track address-energy coordinate system to obtain a target curve; Solving the initial energy address function using the target curve to obtain a constant term coefficient and a quadratic term coefficient, wherein the target coefficient represents the linear term coefficient; Based on the initial energy address function, the target energy address function is obtained according to the constant term coefficient, the target coefficient and the quadratic term coefficient.

7. The method according to claim 2, wherein: The initial weight is shown in formula (1), and the initial energy address function is shown in formula (2); （1） （2） in, is the initial weight of the j-th target interval; is the weighting coefficient of the i-th full energy peak in the j-th target interval; is the differential coefficient at the energy address of the ith full-energy peak in the jth target interval, E is the energy of the γ-ray full-energy peak, ch is the energy address of the γ-ray full-energy peak, are the constant term coefficient, target coefficient and quadratic term coefficient to be solved respectively.

8. A device for determining the energy calibration of a radioactive detector, comprising: an acquisition module, configured to acquire background energy spectra of a plurality of known nuclides, wherein the background energy spectra are obtained by detecting the plurality of known nuclides by the radioactive detector; a differential module, configured to perform differential processing on the background energy spectrum to obtain a channel address-difference map, wherein the channel address-difference map includes a plurality of effective peak information determined from a plurality of full energy peak information; a first determining module, configured to determine a target coefficient based on the track address-difference map, wherein the target coefficient represents a scale coefficient in an initial energy track address function; a second determining module, configured to determine a target energy track address function according to the track address-difference map and the initial energy track address function, wherein the target energy track address function represents the energy scale of the radioactive detector; The first determining module includes: a first determining submodule, configured to determine a plurality of target intervals of a target nuclide from the track address-difference map, wherein the target intervals include a plurality of valid peak information of the target nuclide; A second determining submodule, configured to determine the target coefficient according to the plurality of target intervals; The second determining module includes: A construction submodule, configured to construct a plurality of coordinate points in a channel address-energy coordinate system according to the channel address corresponding to each of the total energy peak information in each of the target intervals and the known total energy peak energy of the target nuclide; The solving submodule is used to solve the initial energy address function based on the target coefficient and the plurality of coordinate points to obtain the target energy address function.

9. An electronic device comprising: one or more processors; a storage device for storing one or more programs, When the one or more programs are executed by the one or more processors, the one or more processors are enabled to execute the method according to any one of claims 1 to 7.

10. A computer-readable storage medium having executable instructions stored thereon, which, when executed by a processor, causes the processor to perform the method according to any one of claims 1 to 7.

11. A computer program product, comprising a computer program, wherein when the computer program is executed by a processor, the method according to any one of claims 1 to 7 is implemented.

Images