CN113032732B - Dose rate fitting method and system based on relative error segmentation - Google Patents
Dose rate fitting method and system based on relative error segmentation Download PDFInfo
- Publication number
- CN113032732B CN113032732B CN202110314794.9A CN202110314794A CN113032732B CN 113032732 B CN113032732 B CN 113032732B CN 202110314794 A CN202110314794 A CN 202110314794A CN 113032732 B CN113032732 B CN 113032732B
- Authority
- CN
- China
- Prior art keywords
- fitting
- dose rate
- relative error
- value
- values
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/02—Dosimeters
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Evolutionary Biology (AREA)
- Algebra (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Molecular Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Measurement Of Radiation (AREA)
Abstract
The invention discloses a dose rate fitting method and a dose rate fitting system based on relative error segmentation, wherein the method comprises the following steps of: acquiring detector count values under different dose rate values, and performing data fitting according to the detector count values and the dose rate values to obtain a first fitting curve; acquiring a first fitting interval according to the first fitting curve, and performing secondary fitting on data in the first fitting interval; wherein, the relative error value of the fitting point in the first fitting interval is greater than the relative error threshold value; and repeating the steps until the relative error of each fitting point on the obtained fitting curve is less than or equal to the relative error threshold. The invention aims to provide a dose rate fitting method and system based on relative error segmentation, which can automatically judge the reasonability of each interval of a fitting curve by taking relative errors as judgment bases, and re-fit unreasonable intervals, wherein the relative error between the dose rate value after the segmentation fitting and a true value in the whole dose rate interval is very small.
Description
Technical Field
The invention relates to the technical field of dose rate measurement, in particular to a dose rate fitting method and system based on relative error segmentation.
Background
For dose rate measurement, detectors such as a GM counting tube, a semiconductor and an ionization chamber are mostly adopted, the detectors mostly adopt a pulse counting mode, and in the conversion process of a counting value and dose rate, a least square method is mostly adopted for curve fitting. Because the least square method carries out data fitting with minimum variance, the fitting result is in a high dose rate area, and the relative error is very small; in the low dose rate region, the relative error between the partial regions can be very large.
In order to make the count-to-dose converted values closer to the true values at low dose rates, it is desirable to reduce the relative error of the fit curve in the fractional region, especially the low dose rate region. To reduce the relative error, it is common to employ:
1. performing data fitting in a pre-segmentation mode, segmenting a fitting area, performing metering rate fitting after segmentation, wherein the relative error of a fitting result is related to the segmentation size and experience of the fitting area;
2. the least square method fitting based on relative errors is very complex for function fitting of more than 2 orders, and causes large relative errors in high dose rate areas.
Disclosure of Invention
The invention aims to provide a dose rate fitting method and system based on relative error segmentation, which can automatically judge the reasonability of each interval of a fitting curve by taking relative errors as judgment bases, and re-fit unreasonable intervals, wherein the relative error between the dose rate value after the segmentation fitting and a true value in the whole dose rate interval is very small.
The invention is realized by the following technical scheme:
a dose rate fitting method based on relative error segmentation comprises the following steps:
s1: acquiring detector count values under different dose rate values, and performing data fitting according to the detector count values and the dose rate values to obtain a first fitting curve;
s2: acquiring a first fitting interval according to the first fitting curve, and performing secondary fitting on data in the first fitting interval;
wherein the relative error value of the fitting points in the first fitting interval is greater than a relative error threshold value;
s3: and repeating the step S2 until the relative error of each fitting point on the obtained fitting curve is less than or equal to the relative error threshold.
Preferably, the data fitting is performed using a least squares method.
Preferably, said S1 comprises the following sub-steps:
s11: collecting detector count values under different dose rate values; wherein a plurality of detector count values are collected at any one dose rate value;
s12: acquiring the average counting value of detectors under different dose rate values;
s13: and performing data fitting on the detector counting average value and the dose rate value by using a least square method to obtain the first fitting curve.
Preferably, said S2 comprises the following sub-steps:
s21: obtaining relative error values of all fitting points on the first fitting curve;
s22: comparing the relative error value of each fitting point with the relative error threshold value, and acquiring a first fitting interval according to the comparison result;
s23: and performing data fitting on the detector counting average value and the dose rate value of the first fitting interval by adopting a least square method to obtain a second fitting curve.
A dose rate fitting system based on relative error segmentation comprises a fitting module, a correction module and a judgment module;
the fitting module is used for obtaining detector counting values under different dose rate values and performing data fitting according to the detector counting values and the dose rate values to obtain a first fitting curve;
the correction module is used for acquiring a first fitting interval according to the first fitting curve and performing quadratic fitting on data in the first fitting interval;
wherein the relative error value of the fitting points in the first fitting interval is greater than a relative error threshold value;
the judging module is used for judging whether the relative error of each fitting point on the fitting curve transmitted by the correcting module is larger than a relative error threshold value, if so, the fitting curve is transmitted to the correcting module, and the correcting module corrects the fitting curve again.
Preferably, the fitting module comprises the following processes:
collecting detector count values under different dose rate values; wherein a plurality of detector count values are collected at any one dose rate value;
acquiring the average counting value of detectors under different dose rate values;
and performing data fitting on the detector counting average value and the dose rate value by adopting a least square method to obtain the first fitting curve.
Preferably, the correction module comprises the following processing procedures:
obtaining relative error values of all fitting points on the first fitting curve;
comparing the relative error value of each fitting point with the relative error threshold value, and acquiring a first fitting interval according to the comparison result;
and performing data fitting on the average detector count value and the dose rate value in the first fitting interval by using a least square method to obtain a second fitting curve.
Compared with the prior art, the invention has the following advantages and beneficial effects:
1. the relative error is used as the judgment basis for the fitting rationality of the dose rate conversion data, so that the fitting data result is more appropriate to the actual application requirement;
2. the rationality judgment basis of the dose rate fitting data is a relative error, a relative error threshold value can be preset, dose rate data fitting is carried out according to the set relative error judgment threshold value, and the relative error of a dose rate fitting result is controllable;
3. compared with the traditional pre-segmentation method, the method eliminates the influence of human factors, can realize automatic segmentation of the dose rate value through programming, and has more rationality than the pre-segmentation method;
4. carrying out dose rate numerical fitting from a high dose rate region to a low dose rate region in sequence, and fully utilizing the premise that a least square method has small relative error of a relatively high dose rate point in the fitting region, so that the relative error of a high dose rate section in the fitting dose rate value is very close to a true value;
5. by presetting a relative error judgment threshold value, the final dosage rate relative to the inherent error is mainly determined by the statistical fluctuation of the count after smoothing.
Drawings
The accompanying drawings, which are included to provide a further understanding of the embodiments of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the principles of the invention. In the drawings:
FIG. 1 is a schematic flow diagram of the present invention;
FIG. 2 is a first fit graph of the present invention;
FIG. 3 is a plot of the residuals of a first fitted curve according to the present invention;
FIG. 4 is a second fit graph of the present invention;
FIG. 5 is a plot of the residuals of a second fitted curve according to the present invention;
FIG. 6 is a third fit graph of the present invention;
FIG. 7 is a plot of the residuals of a third fitted curve according to the present invention;
FIG. 8 is a prior art fitting graph
FIG. 9 is a residual plot of a prior art fitted curve.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to examples and the accompanying drawings, and the exemplary embodiments and descriptions thereof are only used for explaining the present invention and are not used as limiting the present invention.
Example 1
A dose rate fitting method based on relative error segmentation, as shown in fig. 1, comprising the steps of:
s1: acquiring detector count values under different dose rate values, and performing data fitting according to the detector count values and the dose rate values to obtain a first fitting curve;
specifically, collecting detector count values at different dose rate values; in order to enable the acquired data to be more accurate, the counting value of the detector is acquired for multiple times under any dosage rate value, the average value is obtained, and the average value is used as the technical value of the detector under the dosage rate value, so that the influence caused by inaccurate acquisition of single data is avoided; at the same time, in order to make the fitted curve closer to the true value, the number of samples should be large enough, preferably larger than 1000 sample values.
S2: acquiring a first fitting interval according to the first fitting curve, and performing quadratic fitting on data in the first fitting interval;
wherein, the relative error value of the fitting point in the first fitting interval is greater than the relative error threshold value; the relative error threshold is set in advance.
Since the fitting curve is fit according to a large amount of data, most data can only be located on or relatively to the fitting curve, and therefore, the error of the data located outside or relatively to the fitting curve is large, so that the fitting correction needs to be performed again. Therefore, in the present application, by obtaining the relative error value of each fitting point on the first fitting curve; and comparing the relative error value of each fitting point with a relative error threshold value, selecting the fitting point with the relative error larger than the relative error threshold value so as to form a first fitting interval, and then performing data refitting on the detector counting average value and the dose rate value of the first fitting interval by adopting a least square method so as to obtain a second fitting curve.
Preferably, in the process of quadratic fitting, dose rate numerical fitting can be carried out from a high dose rate region to a low dose rate region in sequence, and in the quadratic fitting process, the premise that relative errors of relatively high dose rate points in the fitting region are small by a least square method is fully utilized, so that relative errors of high dose rate sections in the fitted dose rate value are very close to a true value.
It should be noted that there may be more than one first fitting interval on one fitting interval, and the first fitting interval may be reasonably set according to the actual distribution of the fitting points.
S3: and S2, repeating the step until the relative error of each fitting point on the obtained fitting curve is less than or equal to the relative error threshold.
In the prior art, when the counting and the dose rate curve are fitted, the least square method is mostly adopted for curve fitting. Since the least square method performs data fitting with the minimum variance (residual variance), the relative error is very large in a partial interval in a low dose rate region, as shown in fig. 8 and 9; if the least square method fitting based on the relative error is adopted, the weight of the least square method fitting becomes the relative error, so the weight of the least square method fitting in the high dosage rate interval is the same as that of the least square method fitting in the low dosage rate interval, and the relative error of the high dosage rate interval is larger.
Based on this, in order to solve the problem that the relative error of the sub-regions is large after data fitting, in this embodiment, a method for fitting the dose rate by using a least square method is constructed, and the relative error is used as a judgment basis, so that reasonability of each region of a fitting curve can be automatically judged, and refitting is performed on an unreasonable region again, and the relative error between the dose rate value after the piecewise fitting and a true value in the whole dose rate region is very small. The invention is implemented with the relative error between the fitted dose rate value and the true value. On one hand, the fitted dose rate value is closer to the true value than the dose rate value fitted by the traditional method, and the relative error is obviously improved; on the other hand, the fitting relative error judgment threshold value can be adjusted to obtain a controllable relative error.
Example 2
The embodiment provides a relative error segmentation-based dose rate fitting system which comprises a fitting module, a correcting module and a judging module;
the fitting module is used for obtaining detector count values under different dose rate values and performing data fitting according to the detector count values and the dose rate values to obtain a first fitting curve;
the correction module is used for acquiring a first fitting interval according to the first fitting curve and performing quadratic fitting on data in the first fitting interval;
wherein, the relative error value of the fitting point in the first fitting interval is greater than the relative error threshold value;
and the judging module is used for judging whether the relative error of each fitting point on the fitting curve transmitted by the correcting module is greater than the relative error threshold value or not, and if so, transmitting the fitting curve to the correcting module.
Further, the fitting module includes the following processes:
collecting detector count values under different dose rate values; wherein, a plurality of detector count values are collected under any one dose rate value;
acquiring the average counting value of detectors under different dose rate values;
and performing data fitting on the detector counting average value and the dose rate value by adopting a least square method to obtain the first fitting curve.
Further, the correction module comprises the following processing procedures:
obtaining relative error values of all fitting points on the first fitting curve;
comparing the relative error value of each fitting point with a relative error threshold value, and acquiring a first fitting interval according to the comparison result;
and performing data fitting on the average value of the counting number of the detector and the dose rate value in the first fitting interval by adopting a least square method to obtain a second fitting curve.
Example 3
The technical scheme of the application is further explained by taking a GM tube ZP1202 as an example as follows:
the detector count values for different dose rate values are shown in table 1:
TABLE 1 Detector count values at different dose rate values
According to the data, a first fitting curve fitted by a least square method is as follows:
f(x)=p1*x^4+p2*x^3+p3*x^2+p4*x+p5;
wherein p1=6.978e-12; p2= -2.816e-08; p3= -3.421e-05; p4=0.9428; p5= -3.529;
the first fitted graph is shown in fig. 2, and a fitted residual graph can be obtained according to the first fitted graph or the first fitted graph, which is shown in fig. 3; from FIG. 3, it can be determined that the first fit interval is [0.38, 326.789]; fitting again based on the data in this interval to obtain a second fitted curve:
f(x)=p1*x^4+p2*x^3+p3*x^2+p4*x+p5;
wherein p1=3.951e-08; p2= -1.408e-05; p3=0.001631; p4=0.5095; p5= -0.01991;
the second fitted graph is shown in fig. 4, and a fitted residual error map of the second fitted graph can be obtained according to the second fitted graph or the second fitted graph, which is shown in fig. 5; from FIG. 5, it can be determined that the quadratic fit interval is [0.38, 87.47]; fitting again based on the data in this interval to obtain a third fitted curve:
f(x)=p1*x^4+p2*x^3+p3*x^2+p4*x+p5;
wherein, p1=2.899e-06, p2= -0.000431, p3=0.01705, p4=0.439, p5= -0.06927;
as can be seen from fig. 7, the residual error after the third fitting is negligible, and the relative error is obtained from the residual error, so the relative error is also very small and negligible.
Finally, the transformation consists of three segments of data fitting, respectively: [0.38, 87.47],[87.47, 326.789],[326.789, 4604.68].
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are merely exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (5)
1. A dose rate fitting method based on relative error segmentation is characterized by comprising the following steps:
s1: acquiring detector count values under different dose rate values, and performing data fitting according to the detector count values and the dose rate values to obtain a first fitting curve;
wherein, the step S1 includes the following substeps:
s11: collecting detector count values under different dose rate values; wherein, a plurality of detector count values are collected under any one dose rate value;
s12: acquiring the average counting value of detectors under different dose rate values;
s13: performing data fitting on the detector count average value and the dose rate value by using a least square method to obtain a first fitting curve;
s2: acquiring a first fitting interval according to the first fitting curve, and performing secondary fitting on data in the first fitting interval; the relative error value of the fitting points in the first fitting interval is greater than a relative error threshold value;
wherein, the step S2 includes the following substeps:
s21: obtaining relative error values of all fitting points on the first fitting curve;
s22: comparing the relative error value of each fitting point with the relative error threshold value, and acquiring the first fitting interval according to the comparison result;
s23: performing data fitting on the average detector count value and the dose rate value in the first fitting interval by using a least square method to obtain a second fitting curve;
s3: and S2, repeating the step until the relative error of each fitting point on the obtained fitting curve is less than or equal to the relative error threshold.
2. A relative error segmentation based dose rate fitting method as claimed in claim 1 wherein the data fitting is performed using a least squares method.
3. A dose rate fitting system based on relative error segmentation is characterized by comprising a fitting module, a correcting module and a judging module;
the fitting module is used for obtaining detector counting values under different dose rate values and performing data fitting according to the detector counting values and the dose rate values to obtain a first fitting curve;
the correction module is used for acquiring a first fitting interval according to the first fitting curve and performing quadratic fitting on data in the first fitting interval;
wherein the relative error value of the fitting points in the first fitting interval is greater than a relative error threshold value;
the judging module is used for judging whether the relative error of each fitting point on the fitting curve transmitted by the correcting module is larger than a relative error threshold value, if so, the fitting curve is transmitted to the correcting module, and the correcting module corrects the fitting curve again.
4. A relative error segment based dose rate fitting system as defined in claim 3, wherein the fitting module comprises the processes of:
collecting detector count values under different dose rate values; wherein, a plurality of detector count values are collected under any one dose rate value;
acquiring the average counting value of detectors under different dose rate values;
and performing data fitting on the detector counting average value and the dose rate value by adopting a least square method to obtain the first fitting curve.
5. The relative error segment based dose rate fitting system of claim 4 wherein the correction module comprises the processes of:
obtaining relative error values of all fitting points on the first fitting curve;
comparing the relative error value of each fitting point with the relative error threshold value, and acquiring a first fitting interval according to a comparison result;
and performing data fitting on the average detector count value and the dose rate value in the first fitting interval by using a least square method to obtain a second fitting curve.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110314794.9A CN113032732B (en) | 2021-03-24 | 2021-03-24 | Dose rate fitting method and system based on relative error segmentation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110314794.9A CN113032732B (en) | 2021-03-24 | 2021-03-24 | Dose rate fitting method and system based on relative error segmentation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113032732A CN113032732A (en) | 2021-06-25 |
CN113032732B true CN113032732B (en) | 2023-01-24 |
Family
ID=76473954
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110314794.9A Active CN113032732B (en) | 2021-03-24 | 2021-03-24 | Dose rate fitting method and system based on relative error segmentation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113032732B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113686280B (en) * | 2021-09-08 | 2023-12-19 | 中国航发沈阳黎明航空发动机有限责任公司 | Method for evaluating sectional fitting of air inlet and outlet edges of blade profile |
CN116593548A (en) * | 2023-05-18 | 2023-08-15 | 深圳市法莱茵科技有限公司 | Diesel mechanical sulfide detection method and system |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2044461B1 (en) * | 2006-07-10 | 2020-05-06 | University Health Network | Apparatus and methods for real-time verification of radiation therapy |
CN101797423B (en) * | 2010-04-12 | 2011-06-15 | 江苏海明医疗器械有限公司 | Self-correcting multi-cure-fitting digital dosage monitoring method and system |
GB201318995D0 (en) * | 2013-10-28 | 2013-12-11 | Kromek Ltd | Method of spectral data detection and manipulation |
CN104166153B (en) * | 2014-07-25 | 2016-08-24 | 北京永新医疗设备有限公司 | The measuring method of the two dimension angular distribution of radioactive substance radiation dose rate and equipment |
CN104636634B (en) * | 2015-01-04 | 2017-08-25 | 中国科学院光电技术研究所 | A kind of method screened signal data point in cavity attenuation and vibration technique and correct fitting result |
CN104716903A (en) * | 2015-03-17 | 2015-06-17 | 上海交通大学 | Automatic segmentation polynomial fitting method applicable to photovoltaic array output characteristic curve |
CN106814384B (en) * | 2015-11-27 | 2019-03-15 | 华北电力大学 | Nuclear power plant's point source radiation source strength backstepping method and point source radiation source strength backstepping system |
CN106052670B (en) * | 2016-05-27 | 2018-12-21 | 北京航空航天大学 | A kind of interference formula digital closed-loop optic fiber gyroscope Active anti-radiation design method |
CN109350098B (en) * | 2018-08-27 | 2021-02-26 | 苏州瑞派宁科技有限公司 | Method for determining fitting mode of signal, method for reconstructing signal and device |
CN110376937B (en) * | 2019-06-26 | 2021-06-25 | 中国船舶重工集团公司第七一九研究所 | Time sequence control logic and signal processing method suitable for quasi-real-time dose rate measuring device based on pulse light release technology |
CN110988967B (en) * | 2019-12-17 | 2021-10-26 | 北京中科核安科技有限公司 | Method for expanding dosage rate range of environment X and gamma radiation detector |
-
2021
- 2021-03-24 CN CN202110314794.9A patent/CN113032732B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN113032732A (en) | 2021-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113032732B (en) | Dose rate fitting method and system based on relative error segmentation | |
Crow | Confidence intervals for a proportion | |
CN105676817B (en) | The statistical process control method of different size sample average-standard deviation control figure | |
CN110134919B (en) | Method for cleaning abnormal data of wind turbine generator | |
CN102841233B (en) | Secondary current compensation method for current mutual inductors | |
CN105117725B (en) | A kind of car plate self-correction recognition methods based on car plate confidence level | |
CN101674398B (en) | Real-time correcting method for dark current of CCD digital image data | |
CN107911185B (en) | Method for calculating highest available frequency of short wave link in ionosphere storm period | |
CN106018223A (en) | Method for improving measurement precision of air quality detection equipment by segmentation optimal calibration | |
CN106529582B (en) | For the prior probability appraisal procedure for introducing expert opinion in Bayesian network | |
CN113610331A (en) | Cost accounting method and device based on information waterfall and storage medium | |
CN107238442A (en) | Without the detector thermometric method for regulating temperature of temperature stabilizer | |
US20060243916A1 (en) | Spectrometer with charge-carrier-trapping correction | |
CN106303309B (en) | Image correction system and method | |
CN110987875B (en) | Water turbidity detection method for spatial fusion based on support degree and confidence degree | |
CN110763808B (en) | Factory calibration-free method for gas detector | |
CN115510373A (en) | Power supply calibration method based on least square quadratic polynomial fitting | |
CN114779731A (en) | Intelligent manufacturing-oriented production data dynamic monitoring and analyzing system and method | |
CN112285635A (en) | Terminal correction method and system of external current transformer | |
CN112485748B (en) | Single-phase ammeter phase judging method | |
CN113673759A (en) | Real-time marshalling method and terminal for hydrological data | |
CN109975845B (en) | Time service positioning system based on Beidou system and time service positioning method thereof | |
CN109459925B (en) | Method for hand-drawing identification and parameter setting of PID control data | |
Verschueren et al. | Order extraction and background subtraction on CCD CASPEC echelle spectra | |
CN113139405A (en) | Cell classification result correction method based on multimodal fitting analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |