Nuclear energy spectrum stabilizing method and electronic equipment

Abstract

The invention provides a spectrum stabilization method of a nuclear energy spectrum and electronic equipment, and the method comprises the steps: determining a channel address feature value corresponding to each channel address in an energy spectrum diagram; determining a target reference source characteristic peak according to the plurality of channel address characteristic values; determining a reference position with the maximum slope in a falling edge region on the right side of the target reference source characteristic peak; according to a predetermined proportionality coefficient and the reference position, determining a central position of a target reference source characteristic peak; wherein the predetermined proportionality coefficient is the ratio of the reference distance to the full width at half maximum; the reference distance is the distance between the center position of the reference source characteristic peak and the reference position; and based on the central position of the characteristic peak of the target reference source, correcting the energy spectrum diagram to obtain a spectrum-stabilized energy spectrum. According to the technical scheme of the embodiment of the invention, the central position of the characteristic peak of the target reference source can be accurately determined, the influence of statistical fluctuation and the interference of characteristic peaks of other radionuclides are avoided, and the stability of a spectrum stabilization result is ensured.

Description

Technical Field

[0001] The present application relates to the technical field of nuclear radiation monitoring, and in particular to a method for stabilizing a nuclear energy spectrum and electronic equipment. Background Art

[0002] In the related technology, in the process of nuclear energy spectrum stabilization, it is first necessary to find the characteristic peak of potassium-40 (K-40) in the natural background, and then adjust the peak position of the characteristic peak to the channel expected by the user, and use the pre-calibrated energy curve to perform energy mapping on other channels. An important reason why the traditional spectrum stabilization method uses the characteristic peak of K-40 as the reference peak is that the characteristic peak energy of K-40 is higher than the energy of most artificial radioactive nuclides. However, due to the limitation of the resolution ability of the detector, the characteristic peak energy of some nuclides is close to the characteristic peak energy of K-40, and the broadened characteristic peak will interfere with the left half of the K-40 characteristic peak. At this time, it is difficult to accurately locate the peak position of the K-40 characteristic peak using the traditional spectrum stabilization method. Summary of the invention

[0003] The embodiment of the present application provides a nuclear energy spectrum stabilization method and electronic device to solve the problems existing in the related technology. The technical solution is as follows:

[0004] In a first aspect, an embodiment of the present application provides a method for stabilizing a nuclear energy spectrum, comprising: determining a channel address characteristic value corresponding to each channel address in an energy spectrum diagram; determining a target reference source characteristic peak based on multiple channel address characteristic values; determining a reference position with the largest slope in the right falling edge area of ​​the target reference source characteristic peak; determining a center position of the target reference source characteristic peak based on a predetermined proportionality coefficient and a reference position; wherein the predetermined proportionality coefficient is a ratio of a reference distance to a half-width at half-maximum; the reference distance is a distance between the center position of the reference source characteristic peak and the reference position; based on the center position of the target reference source characteristic peak, the energy spectrum diagram is corrected to obtain a stabilized energy spectrum.

[0005] In one embodiment, determining a reference position with the largest slope in a falling edge region on the right side of a characteristic peak of a target reference source includes: calculating original count derivative values ​​corresponding to each channel address in the falling edge region on the right side of a characteristic peak of a target reference source; smoothing the original count derivative values ​​to obtain mean count derivative values ​​corresponding to each channel address; and determining a position corresponding to the maximum value in the mean count derivative values ​​as a reference position with the largest slope.

[0006] In one embodiment, in the falling edge area on the right side of the characteristic peak of the target reference source, the original count derivative value corresponding to each channel address is calculated, including: in the falling edge area on the right side of the characteristic peak of the target reference source, a first count value located at a predetermined distance to the left of the target channel address and a second count value located at a predetermined distance to the right of the target channel address are calculated; wherein the predetermined distance is determined based on the half-height width corresponding to the channel address; and the original count derivative value corresponding to the target channel address is determined based on the first count value and the second count value.

[0007] In one embodiment, the raw count derivative value is smoothed to obtain the mean count derivative value corresponding to each track address, including: calculating the total number of tracks within a preset track address range of the target track address and the sum of the raw count derivative values; wherein the preset track address range is a range between a predetermined distance to the left of the target track address and a predetermined distance to the right of the target track address; and determining the mean count derivative value corresponding to the target track address based on the total number of tracks within the preset track address range and the sum of the raw count derivative values.

[0008] In one embodiment, determining a target reference source characteristic peak based on multiple channel address characteristic values ​​includes: determining a local maximum value among the multiple channel address characteristic values; determining a characteristic peak corresponding to a local maximum value greater than a preset threshold as a pre-selected reference source characteristic peak; and determining a target reference source characteristic peak among the pre-selected reference source characteristic peaks.

[0009] In one embodiment, determining a target reference source characteristic peak from among preselected reference source characteristic peaks includes: determining the count rate of the preselected reference source characteristic peak; determining a preselected reference source characteristic peak whose count rate is within a preset count rate range as a candidate reference source characteristic peak; and determining a candidate reference source characteristic peak that is closest to a reference source characteristic peak after a previous spectrum stabilization adjustment as a target reference source characteristic peak.

[0010] In one embodiment, determining the channel address characteristic value corresponding to each channel address in the energy spectrum diagram includes: determining the count value corresponding to each channel address in the energy spectrum diagram; performing convolution calculation with the count value corresponding to each channel address using a symmetric window function to obtain the channel address characteristic value corresponding to each channel address.

[0011] In a second aspect, an embodiment of the present application provides an electronic device, including a memory, a processor, and a computer program stored in the memory, and the processor implements any method of the embodiment of the present application when executing the computer program.

[0012] In a third aspect, an embodiment of the present application provides a computer-readable storage medium, in which a computer program is stored. When the computer program is executed by a processor, any method of the embodiment of the present application is implemented.

[0013] In a fourth aspect, an embodiment of the present application provides a computer program product, including a computer program, which implements any method of the embodiments of the present application when executed by a processor.

[0014] The advantages or beneficial effects of the above technical solution include at least: by determining the reference position with the largest slope in the falling edge area on the right side of the characteristic peak of the target reference source, and accurately determining the center position of the characteristic peak of the target reference source according to the predetermined proportional coefficient and the reference position, the influence of statistical fluctuations and interference from the characteristic peaks of other radioactive nuclides are avoided, thereby ensuring the stability of the stabilization spectrum results.

[0015] The above summary is for illustrative purposes only and is not intended to be limiting in any way. In addition to the illustrative aspects, embodiments and features described above, further aspects, embodiments and features of the present application will be readily apparent by reference to the accompanying drawings and the following detailed description. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] In the accompanying drawings, unless otherwise specified, the same reference numerals throughout the multiple drawings represent the same or similar parts or elements. These drawings are not necessarily drawn to scale. It should be understood that these drawings only depict some embodiments disclosed in the present application and should not be regarded as limiting the scope of the present application.

[0017] Figure 1 An example diagram showing an application of a spectrum stabilization method in the related art;

[0018] Figure 2 A schematic diagram showing the characteristic peaks of lanthanum (La) and K-40 in the related art;

[0019] Figure 3 A schematic diagram showing a flow chart of a method for stabilizing a nuclear energy spectrum according to an embodiment of the present application;

[0020] Figure 4 A block diagram of an electronic device according to an embodiment of the present application is shown. DETAILED DESCRIPTION

[0021] In the following, only some exemplary embodiments are briefly described. As those skilled in the art will appreciate, the described embodiments may be modified in various ways without departing from the spirit or scope of the present application. Therefore, the drawings and descriptions are considered to be exemplary and non-restrictive in nature.

[0022] Figure 1 The following is an example diagram showing the application of spectrum stabilization methods in related technologies. Figure 1 (a) is the energy spectrum of cobalt-60 (Co-60); Figure 1 (b) is the peak search result diagram. Figure 1In the process of nuclear energy spectrum stabilization, we first need to find the characteristic peak of K-40 in the natural background (such as Figure 1 Then adjust the peak position of the characteristic peak to the channel address expected by the user, and use the pre-calibrated energy curve to perform energy mapping on other channels. However, due to the limitation of the detector resolution, the characteristic peak energy of some nuclides is close to the characteristic peak energy of K-40, and the broadened characteristic peak will interfere with the characteristic peak of K-40, making it difficult to accurately locate the peak position of K-40 characteristic peak.

[0023] Taking the energy spectrum of lanthanum bromide (Labr3) detector as an example, the Labr3 detector contains the weak radioactive nuclide lanthanum (La). Figure 2 Schematic diagram showing the characteristic peaks of La and K-40 in the related art. Figure 2 As shown in the figure, the blue arrow points to the characteristic peak of La on the left, and the green arrow points to the characteristic peak of K-40. The energy of the characteristic peak of La is similar to that of the characteristic peak of K-40. Therefore, the characteristic peak of La will have a significant interference in the calculation of the center channel address of the characteristic peak of K-40, resulting in an error in the calculated center channel address of the characteristic peak of K-40.

[0024] The activity components of the radioactive nuclide La in different Labr3 detectors vary greatly, and the activity of K-40 also varies in different measurement environments. Even in the same environment, due to the influence of statistical fluctuations, the superposition of the two characteristic peaks (i.e., the characteristic peak of La mentioned above and the characteristic peak of K-40) also has many different forms. Therefore, no matter what calculation algorithm is used (such as the maximum value method, the first-order derivative method or the multi-order derivative method, the center of gravity method, the symmetrical zero area method, the Gaussian fitting method, etc.), the peak position of the K-40 characteristic peak cannot be accurately calculated.

[0025] For other scintillator detectors (such as sodium iodide (NaI) detectors), although they do not carry their own radioactive nuclides, the resolution of the detectors is poor, the half-width of their characteristic peaks is large, and they will also be interfered by the characteristic peaks of other external radioactive nuclides (such as the 1332Kev characteristic peak of Co-60, the 1408Kev characteristic peak of europium-152 (Eu-152), etc.).

[0026] Figure 3 FIG. 1 is a flow chart showing a method for stabilizing a nuclear energy spectrum according to an embodiment of the present application. Figure 3 As shown, the stabilization method of the nuclear energy spectrum includes:

[0027] Step S301: determining the channel characteristic value corresponding to each channel in the energy spectrum;

[0028] Step S302: Determine a target reference source characteristic peak according to a plurality of channel address characteristic values.

[0029] For example, the address characteristic value S i It can be understood as the local characteristic intensity of the energy spectrum data at channel address i, which reflects the signal change characteristics of channel address i and its adjacent areas, and can be used to identify characteristic peaks. Each channel address has a corresponding channel address characteristic value, and the target reference source characteristic peak can be determined based on multiple channel address characteristic values ​​corresponding to multiple channels.

[0030] Exemplarily, the reference source may be the natural radionuclide K-40, and the characteristic peak of the reference source is the characteristic peak of K-40.

[0031] Step S303: Determine the reference position with the largest slope in the falling edge region on the right side of the characteristic peak of the target reference source.

[0032] The count value of the signal corresponding to the top of the characteristic peak is the highest, so its statistical fluctuation is relatively significant, resulting in the peak position being greatly affected by noise. In the right falling edge area, the signal strength is lower than the peak, the relative impact of noise is also smaller, and the position with the largest slope is determined by the geometric change of the signal, rather than the simple count value, so the reference position is less affected by statistical fluctuations. Therefore, compared with the peak position or the left half position, the position with the largest slope in the right falling edge area of ​​the characteristic peak is more stable.

[0033] Among them, statistical fluctuations refer to signal fluctuations caused by the random emission of radiation source particles and the uncertainty of the detector (such as counting noise).

[0034] Step S304: determining the center position of the characteristic peak of the target reference source according to the predetermined scale factor and the reference position;

[0035] Step S305: Based on the center position of the characteristic peak of the target reference source, the energy spectrum is corrected to obtain a stable spectrum.

[0036] Among them, the predetermined proportional coefficient ratio is the reference distance d k40-edge The ratio of the half-width FWHM is Reference distance d k40-edge is the center position d of the characteristic peak of the reference source k40 With reference position d edge Therefore, the reference distance d k40-edge is the product of the predetermined ratio and the half-height width FWHM. Wherein, the half-height width FWHM is the product of the resolution of the detector and the energy corresponding to the current channel address. k40-edge After that, the reference position d can be calculated edge Distance d from reference k40-edge The difference between the center position d of the target reference source characteristic peak is obtained k40 For example, the predetermined proportional coefficient ratio may be a calibration parameter of the detector.

[0037] According to the nuclear energy spectrum stabilization method of the embodiment of the present application, by determining the reference position with the largest slope in the falling edge area on the right side of the characteristic peak of the target reference source, and accurately determining the center position of the characteristic peak of the target reference source according to a predetermined proportional coefficient and the reference position, the influence of statistical fluctuations and interference from the characteristic peaks of other radioactive nuclides are avoided, thereby ensuring the stability of the stabilization result.

[0038] In one embodiment, in step S303, determining the reference position with the largest slope in the falling edge area on the right side of the characteristic peak of the target reference source may include: calculating the original count derivative value corresponding to each channel address in the falling edge area on the right side of the characteristic peak of the target reference source; smoothing the original count derivative value to obtain the mean count derivative value corresponding to each channel address; and determining the position corresponding to the maximum value in the mean count derivative value as the reference position with the largest slope.

[0039] It should be noted that the above-mentioned "count derivative value" can be understood as the derivative value of the count value, which is used to indicate the speed of the change of the count value. The "count derivative value" in this application is an absolute value.

[0040] In this embodiment, by smoothing the original count derivative value to obtain the mean count derivative value corresponding to each channel address, the error caused by statistical fluctuations can be reduced, thereby reducing noise and local fluctuations, making the obtained reference position more stable and reliable.

[0041] In one embodiment, calculating the original count derivative value corresponding to each track address in the right falling edge region of the target reference source characteristic peak may include: calculating a first count value located at a predetermined distance to the left of the target track address and a second count value located at a predetermined distance to the right of the target track address in the right falling edge region of the target reference source characteristic peak; wherein the predetermined distance is determined based on the half-height width corresponding to the track address; and determining the original count derivative value corresponding to the target track address based on the first count value and the second count value.

[0042] For example, the half-width FWHM may be an integer multiple of the predetermined distance. Take the example for illustration. Then the predetermined distance to the left of the target address i is The first count value at a predetermined distance to the left of the target track i is The predetermined distance to the right of the target address i is The second count value at a predetermined distance to the right of the target track i is The original count derivative value corresponding to the target address i

[0043] In one example, in the falling edge region on the right side of the target reference source characteristic peak, the partial channel address i and the count value y iThe relationship is shown in Table 1:

[0044] Table 1

[0045] Address <![CDATA[Count value y i > 12 600 13 550 14 400 15 250 16 100

[0046] For example, when the half-width FWHM is 3, the predetermined distance is 1. When calculating the raw count derivative value corresponding to track address 13, the target track address i is 13. The predetermined distance to the left of the target track address 13 is 13-1=12, that is, track address 12. The first count value y at the predetermined distance to the left of the target track address 13 is 12 is 600. The predetermined distance to the right of the target track address 13 is 13+1=14, that is, track address 14. The second count value y at the predetermined distance to the right of the target track address 13 is 14 The original count derivative value corresponding to the target address 13 is 400. Similarly, the raw count derivative value corresponding to each channel address in the falling edge area on the right side of the characteristic peak of the target reference source can be calculated.

[0047] In this embodiment, the raw count derivative value corresponding to each channel address can be effectively determined, so that the raw count derivative value corresponding to each channel address can be smoothed to obtain the mean count derivative value corresponding to each channel address, and then the reference position with the largest slope can be determined.

[0048] In one embodiment, smoothing the raw count derivative values ​​to obtain the mean count derivative values ​​corresponding to each track address may include: calculating the total number of tracks within a preset track address range of the target track address and the sum of the raw count derivative values; wherein the preset track address range is a range between a predetermined distance to the left of the target track address and a predetermined distance to the right of the target track address; and determining the mean count derivative value corresponding to the target track address based on the total number of tracks within the preset track address range and the sum of the raw count derivative values.

[0049] For example, the predetermined distance is Take the example for illustration. Then the predetermined distance to the left of the target address i is The predetermined distance to the right of the target address i is The default address range is (Including the endpoint value). The total number of addresses within the preset address range is The sum of the raw count derivative values ​​within the preset address range is The mean count derivative value corresponding to the target address i

[0050] In one example, in the falling edge region on the right side of the target reference source characteristic peak, the partial channel address i and the raw count derivative value The relationship is shown in Table 2:

[0051] Table 2

[0052]

[0053] For example, when the half-width FWHM is 6, the predetermined distance is 2. When calculating the mean count derivative value corresponding to track address 14, the target track address i is 14. The predetermined distance to the left of the target track address 14 is 14-2=12, that is, track address 12. The predetermined distance to the right of the target track address 14 is 14+2=16, that is, track address 16. The preset track address range is track address 12 to track address 16 (including the endpoint values). The total number of tracks within the preset track address range is 5.

[0054] The mean count derivative value f corresponding to the target address 14 avg-14 =100 / 5=20. Similarly, the mean count derivative value corresponding to each channel address in the falling edge area on the right side of the target reference source characteristic peak can be calculated. The position corresponding to the maximum value in the mean count derivative value is determined as the reference position with the maximum slope.

[0055] In this embodiment, the original count derivative values ​​within the preset track address range of the target track address can be average filtered, thereby effectively reducing the error caused by statistical fluctuations, reducing noise and local fluctuations, and making the obtained reference position more stable and reliable.

[0056] In one embodiment, in step S302, determining a target reference source characteristic peak based on multiple channel address characteristic values ​​may include: determining a local maximum value among the multiple channel address characteristic values; determining a characteristic peak corresponding to a local maximum value greater than a preset threshold as a pre-selected reference source characteristic peak; and determining a target reference source characteristic peak among the pre-selected reference source characteristic peaks.

[0057] The local maximum can be understood as the maximum value in a local range, which is greater than the address characteristic value corresponding to the adjacent address. i >S i-1 And S i >S i+1 In the case of S i is a local maximum.

[0058] In one example, the partial address i and the address feature value S i The relationship is shown in Table 3:

[0059] Table 3

[0060] Address <![CDATA[Road address eigenvalue S i > 1 5 2 10 3 8 4 12 5 7

[0061] For example, in combination with Table 3, the channel address characteristic value S2 corresponding to channel address 2 is 10, S2>S1 and S2>S3, so the channel address characteristic value S2 is a local maximum. The channel address characteristic value S4 corresponding to channel address 4 is 12, S4>S3 and S4>S5, so the channel address characteristic value S4 is also a local maximum. Similarly, all local maxima in the energy spectrum can be determined. Then, the local maximum values ​​greater than the preset threshold are screened out from all local maxima, and are determined as pre-selected reference source characteristic peaks. For example, when the preset threshold is 9, the channel address characteristic value S2 and the channel address characteristic value S4 are both determined as pre-selected reference source characteristic peaks; when the preset threshold is 11, only the channel address characteristic value S4 is determined as the pre-selected reference source characteristic peak; when the preset threshold is 13, neither the channel address characteristic value S2 nor the channel address characteristic value S4 is determined as the pre-selected reference source characteristic peak. Finally, the target reference source characteristic peak is determined from the pre-selected reference source characteristic peaks.

[0062] In this embodiment, by determining the characteristic peak corresponding to the local maximum value greater than a preset threshold as the pre-selected reference source characteristic peak, the insignificant local maximum value can be effectively filtered out, and the significant characteristic peaks can be retained as the pre-selected reference source characteristic peaks, and finally the target reference source characteristic peak can be determined from the pre-selected reference source characteristic peaks.

[0063] In one embodiment, determining a target reference source characteristic peak from among preselected reference source characteristic peaks includes: determining the count rate of the preselected reference source characteristic peak; determining a preselected reference source characteristic peak whose count rate is within a preset count rate range as a candidate reference source characteristic peak; and determining a candidate reference source characteristic peak that is closest to a reference source characteristic peak after a previous spectrum stabilization adjustment as a target reference source characteristic peak.

[0064] Exemplarily, the count rate doserate of the characteristic peak of the preselected reference source can be calculated by the following formula: Among them, y i is the count value of the peak region channel i; FWHM is the half maximum width; is the sum of all count values ​​in the peak area; t is the data acquisition time. The preset count rate range can be 0.3Dose~3Dose (including the endpoint values). Among them, Dose is the standard count rate of the detector in the passive spectrum stabilization state, which can be measured when the detector is powered on and calibrated. After determining the characteristic peak of the candidate reference source, the distance between each candidate reference source characteristic peak and the reference source characteristic peak recorded after the previous spectrum stabilization adjustment can be calculated, and the nearest candidate reference source characteristic peak is determined as the target reference source characteristic peak.

[0065] In this embodiment, by determining the characteristic peak of the preselected reference source whose counting rate is within the preset counting rate range as the characteristic peak of the candidate reference source, the characteristic peaks of the preselected reference source with too low intensity and too high intensity can be excluded, and the influence of noise or background fluctuations and interference caused by other types of radionuclides can be avoided. By determining the characteristic peak of the candidate reference source closest to the characteristic peak of the reference source after the previous spectrum stabilization adjustment as the characteristic peak of the target reference source, the characteristic peak with the smallest deviation from the characteristic peak of the historical reference source can be selected from the characteristic peaks of the candidate reference sources as the characteristic peak of the target reference source, thereby improving the positioning accuracy of the characteristic peak of the target reference source.

[0066] In one embodiment, in step S301, determining the channel address characteristic value corresponding to each channel address in the energy spectrum diagram may include: determining the count value corresponding to each channel address in the energy spectrum diagram; performing convolution calculation with the count value corresponding to each channel address using a symmetric window function to obtain the channel address characteristic value corresponding to each channel address.

[0067] Exemplarily, the symmetric window function C j It is a symmetrical window function with a total area of ​​zero, that is, a symmetrical window function C j satisfy: In this way, the contribution of the background signal can be eliminated, thereby enhancing the ability to identify the characteristic peak. The prototype of the window function can be a Gaussian function, such as a Gaussian function with a height of 1 and a half-height width equal to the half-height width of the detector. Symmetric window function C j satisfy: Among them, σ=FWHM / 2.355, is the standard deviation of the Gaussian function; FWHM is the half-height width corresponding to channel address i. The convolution calculation formula can be satisfied: in, is the weighted sum of the symmetric window function and the count value; To normalize the weighted sum and reduce the impact of noise. j The formula and convolution formula can calculate the address characteristic value corresponding to each address.

[0068] In this embodiment, the channel address characteristic value corresponding to each channel address can be obtained by using the symmetric zero area method, so that the target reference source characteristic peak can be determined according to the multiple channel address characteristic values.

[0069] As an implementation of the above methods, an embodiment of the present application also provides a nuclear energy spectrum stabilization device, which may include: a characteristic value determination module, used to determine the channel address characteristic value corresponding to each channel address in the energy spectrum diagram; a characteristic peak determination module, used to determine the target reference source characteristic peak according to multiple channel address characteristic values; a reference position determination module, used to determine the reference position with the largest slope in the right falling edge area of ​​the target reference source characteristic peak; a center position determination module, used to determine the center position of the target reference source characteristic peak according to a predetermined proportionality coefficient and a reference position; wherein the predetermined proportionality coefficient is the ratio of the reference distance to the half-height width; the reference distance is the distance between the center position of the reference source characteristic peak and the reference position; a correction module, used to correct the energy spectrum diagram based on the center position of the target reference source characteristic peak to obtain a stabilized spectrum.

[0070] In one embodiment, the reference position determination module includes: a calculation submodule, which is used to calculate the original count derivative value corresponding to each channel address in the falling edge area on the right side of the characteristic peak of the target reference source; a processing submodule, which is used to smooth the original count derivative value to obtain the mean count derivative value corresponding to each channel address; and a position determination submodule, which is used to determine the position corresponding to the maximum value in the mean count derivative value as the reference position with the largest slope.

[0071] In one embodiment, the calculation submodule is also used to: calculate a first count value located at a predetermined distance to the left of the target track address and a second count value located at a predetermined distance to the right of the target track address in the falling edge area on the right side of the characteristic peak of the target reference source; wherein the predetermined distance is determined based on the half-height width corresponding to the track address; and determine the original count derivative value corresponding to the target track address based on the first count value and the second count value.

[0072] In one embodiment, the processing submodule is further used to: calculate the total number of addresses within a preset address range of the target address and the sum of the original count derivative values; wherein the preset address range is a range between a predetermined distance to the left of the target address and a predetermined distance to the right of the target address; and determine the mean count derivative value corresponding to the target address based on the total number of addresses within the preset address range and the sum of the original count derivative values.

[0073] In one embodiment, the characteristic peak determination module includes: a maximum value determination submodule, which is used to determine the local maximum value among multiple channel address characteristic values; a pre-selected characteristic peak determination submodule, which is used to determine the characteristic peak corresponding to the local maximum value greater than a preset threshold as a pre-selected reference source characteristic peak; and a target characteristic peak determination submodule, which is used to determine the target reference source characteristic peak among the pre-selected reference source characteristic peaks.

[0074] In one embodiment, the target characteristic peak determination submodule is also used to: determine the count rate of the preselected reference source characteristic peak; determine the preselected reference source characteristic peak whose count rate is within a preset count rate range as the candidate reference source characteristic peak; determine the candidate reference source characteristic peak that is closest to the reference source characteristic peak after the previous spectrum stabilization adjustment as the target reference source characteristic peak.

[0075] In one embodiment, the characteristic value determination module is also used to: determine the count value corresponding to each channel address in the energy spectrum; use a symmetric window function and the count value corresponding to each channel address to perform convolution calculation to obtain the channel address characteristic value corresponding to each channel address.

[0076] The functions of each module in the device of the embodiment of the present invention can refer to the corresponding description in the above method, which will not be repeated here.

[0077] According to an embodiment of the present application, the present application also provides an electronic device, a computer-readable storage medium and a computer program product.

[0078] Figure 4 FIG. 2 shows a structural block diagram of an electronic device according to an embodiment of the present invention. Figure 4 As shown, the electronic device includes: a memory 410 and a processor 420. The memory 410 stores a computer program that can be run on the processor 420. When the processor 420 executes the computer program, the radiation source positioning method in the above embodiment is implemented. The number of the memory 410 and the processor 420 can be one or more.

[0079] The electronic device also includes:

[0080] The communication interface 430 is used to communicate with external devices and perform data exchange transmission.

[0081] If the memory 410, the processor 420 and the communication interface 430 are implemented independently, the memory 410, the processor 420 and the communication interface 430 can be connected to each other through a bus and communicate with each other. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus. The bus can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 4 Only one thick line is used in the diagram, but this does not mean that there is only one bus or only one type of bus.

[0082] Optionally, in a specific implementation, if the memory 410, the processor 420 and the communication interface 430 are integrated on a chip, the memory 410, the processor 420 and the communication interface 430 can communicate with each other through an internal interface.

[0083] An embodiment of the present invention provides a computer-readable storage medium storing a computer program, which implements the method provided in the embodiment of the present application when the program is executed by a processor.

[0084] An embodiment of the present application also provides a chip, which includes a processor for calling and executing instructions stored in the memory from the memory, so that a communication device equipped with the chip executes the method provided by the embodiment of the present application.

[0085] An embodiment of the present application also provides a chip, including: an input interface, an output interface, a processor and a memory, wherein the input interface, the output interface, the processor and the memory are connected via an internal connection path, and the processor is used to execute the code in the memory. When the code is executed, the processor is used to execute the method provided in the embodiment of the application.

[0086] It should be understood that the processor may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or any conventional processor, etc. It is worth noting that the processor may be a processor supporting the advanced reduced instruction set machine (ARM) architecture.

[0087] Further, optionally, the above-mentioned memory may include a read-only memory and a random access memory, and may also include a non-volatile random access memory. The memory may be a volatile memory or a non-volatile memory, or may include both volatile and non-volatile memories. Among them, the non-volatile memory may include a read-only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable programmable read-only memory (EEPROM), or a flash memory. The volatile memory may include a random access memory (RAM), which is used as an external cache. By way of example but not limitation, many forms of RAM are available. For example, static random access memory (static RAM, SRAM), dynamic random access memory (dynamic random access memory, DRAM), synchronous dynamic random access memory (synchronous DRAM, SDRAM), double data rate synchronous dynamic random access memory (double data rate SDRAM, DDR SDRAM), enhanced synchronous dynamic random access memory (enhanced SDRAM, ESDRAM), synchronous link dynamic random access memory (synchlink DRAM, SLDRAM) and direct memory bus random access memory (direct RAM bus, DR RAM).

[0088] In the above embodiments, it can be implemented in whole or in part by software, hardware, firmware or any combination thereof. When implemented using software, it can be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on a computer, the process or function according to the present application is generated in whole or in part. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable devices. The computer instructions can be stored in a computer-readable storage medium, or transmitted from one computer-readable storage medium to another computer-readable storage medium.

[0089] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine different embodiments or examples described in this specification and the features of different embodiments or examples, unless they are contradictory.

[0090] In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Therefore, a feature defined as "first" or "second" may explicitly or implicitly include at least one of the features. In the description of this application, the meaning of "plurality" is two or more, unless otherwise clearly and specifically defined.

[0091] Any process or method description in the flow chart or otherwise described herein can be understood to represent a module, fragment or portion of a code including one or more executable instructions for implementing the steps of a specific logical function or process. And the scope of the preferred embodiment of the present application includes other implementations, in which the functions may not be performed in the order shown or discussed, including in a substantially simultaneous manner or in a reverse order according to the functions involved.

[0092] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing logical functions, which can be embodied in any computer-readable medium for use by an instruction execution system, apparatus or device (such as a computer-based system, a system including a processor or other system that can fetch instructions from an instruction execution system, apparatus or device and execute instructions), or used in combination with these instruction execution systems, apparatuses or devices.

[0093] It should be understood that the various parts of the present application can be implemented with hardware, software, firmware or a combination thereof. In the above embodiments, multiple steps or methods can be implemented with software or firmware stored in a memory and executed by a suitable instruction execution system. All or part of the steps of the above embodiment method can be completed by instructing the relevant hardware through a program, which can be stored in a computer-readable storage medium, and when the program is executed, it includes one of the steps of the method embodiment or a combination thereof.

[0094] In addition, each functional unit in each embodiment of the present application can be integrated into a processing module, or each unit can exist physically separately, or two or more units can be integrated into one module. The above-mentioned integrated module can be implemented in the form of hardware or in the form of a software functional module. If the above-mentioned integrated module is implemented in the form of a software functional module and sold or used as an independent product, it can also be stored in a computer-readable storage medium. The storage medium can be a read-only memory, a disk or an optical disk, etc.

[0095] The above is only a specific implementation of the present application, but the protection scope of the present application is not limited thereto. Any technician familiar with the technical field can easily think of various changes or substitutions within the technical scope disclosed in the present application, which should be included in the protection scope of the present application. Therefore, the protection scope of the present application should be based on the protection scope of the claims.

Claims

1. A method for stabilizing a nuclear energy spectrum, characterized in that: include: Determine the channel address characteristic value corresponding to each channel address in the energy spectrum; Determine a target reference source characteristic peak according to a plurality of channel address characteristic values; Determine a reference position with the largest slope in the falling edge region on the right side of the characteristic peak of the target reference source; Determine the center position of the characteristic peak of the target reference source according to a predetermined proportionality coefficient and the reference position; wherein the predetermined proportionality coefficient is the ratio of the reference distance to the half-height width; and the reference distance is the distance between the center position of the characteristic peak of the reference source and the reference position; Based on the center position of the characteristic peak of the target reference source, the energy spectrum is corrected to obtain a stable spectrum.

2. The method according to claim 1, characterized in that Determining a reference position with the largest slope in a falling edge region on the right side of the characteristic peak of the target reference source comprises: In the falling edge region on the right side of the characteristic peak of the target reference source, calculating the raw count derivative value corresponding to each of the channel addresses; Smoothing the original count derivative values ​​to obtain mean count derivative values ​​corresponding to the channel addresses; The position corresponding to the maximum value in the mean count derivative value is determined as the reference position with the maximum slope.

3. The method according to claim 2, characterized in that In the falling edge region on the right side of the characteristic peak of the target reference source, the raw count derivative value corresponding to each of the channel addresses is calculated, including: In the falling edge region on the right side of the characteristic peak of the target reference source, a first count value located at a predetermined distance to the left of the target track address and a second count value located at a predetermined distance to the right of the target track address are calculated; wherein the predetermined distance is determined based on the half-height width corresponding to the track address; The original count derivative value corresponding to the target address is determined according to the first count value and the second count value.

4. The method according to claim 2, characterized in that: The raw count derivative values ​​are smoothed to obtain the mean count derivative values ​​corresponding to the channel addresses, including: Calculate the total number of tracks within a preset track range of the target track and the sum of the raw count derivative values; wherein the preset track range is a range between a predetermined distance to the left of the target track and a predetermined distance to the right of the target track; The mean count derivative value corresponding to the target address is determined according to the total number of addresses within the preset address range and the sum of the original count derivative value.

5. The method according to claim 1, characterized in that: Determining a target reference source characteristic peak according to a plurality of channel address characteristic values ​​comprises: Determining a local maximum among a plurality of said track address characteristic values; Determine the characteristic peak corresponding to the local maximum value greater than a preset threshold as the pre-selected reference source characteristic peak; A target reference source characteristic peak is determined from the preselected reference source characteristic peaks.

6. The method according to claim 5, characterized in that Determining a target reference source characteristic peak from the preselected reference source characteristic peaks comprises: determining a count rate of a characteristic peak of the preselected reference source; Determine the preselected reference source characteristic peak whose counting rate is within the preset counting rate range as the reference source characteristic peak to be selected; The characteristic peak of the candidate reference source closest to the characteristic peak of the reference source after the previous spectrum stabilization adjustment is determined as the characteristic peak of the target reference source.

7. The method according to any one of claims 1 to 6, characterized in that Determine the channel characteristic value corresponding to each channel in the energy spectrum, including: Determine the count value corresponding to each channel address in the energy spectrum; The symmetric window function is used to perform convolution calculation with the count value corresponding to each channel address to obtain the channel address characteristic value corresponding to each channel address.

8. An electronic device, comprising: A processor, and a memory storing a program, the program comprising instructions which, when executed by the processor, cause the processor to perform the method according to any one of claims 1 to 7.

9. A computer-readable storage medium, wherein a computer program is stored in the computer-readable storage medium, and when the computer program is executed by a processor, the method according to any one of claims 1 to 7 is implemented.

10. A computer program product, comprising a computer program, which, when executed by a processor of a computer, is used to cause the computer to perform the method of any one of claims 1 to 7.

Images