BNCT-based system and method for analyzing boron 10 in human body

Abstract

The invention relates to the technical field of radioactive isotope analysis, in particular to a BNCT-based system and method for analyzing boron 10 in a human body, and the system comprises a boron 10 distribution detection module, a voxel reconstruction imaging module, a concentration region division module and a concentration correlation analysis module. According to the method, decay frequency, position and time information of boron 10 is captured, comprehensiveness and spatial accuracy of data acquisition are improved, spatial positioning processing of frequency accumulation is refined, frequency values in voxels are associated to spatial coordinates, and it is ensured that concentration value mapping accurately reflects the actual distribution condition of the positions of the voxels; by combining threshold comparison frequency characteristics, marking high and low frequency regions, accurately depicting regional frequency characteristics, supporting refined calculation of concentration difference and frequency characteristics, analyzing regional concentration and frequency difference data through item-by-item operation, deepening differentiated understanding of boron 10 in different regions, remarkably improving boron 10 distribution detection precision and information interpretation depth, and improving boron 10 distribution detection accuracy and information interpretation accuracy. And data support is provided for application optimization.

Description

Technical Field

[0001] The present invention relates to the technical field of radioisotope analysis, and in particular to a system and method for analyzing boron 10 in the human body based on BNCT. Background Art

[0002] The field of radioisotope analysis technology involves the use of radioisotopes for quantitative and qualitative analysis of elements or compounds, and is mainly used to detect the distribution and concentration of trace elements in materials or organisms. This type of technology uses the decay characteristics and radiation characteristics of radioactive markers to generate data through specific imaging or detection methods under non-invasive conditions to achieve accurate measurement of the target substance. Radioisotope analysis is widely used in material composition analysis, environmental monitoring, and biological research due to its high sensitivity and high accuracy. It is particularly suitable for tracking and distribution evaluation of low-content elements or targets that are difficult to measure directly. This technology can also be combined with high-resolution imaging methods to visualize the analysis data, thereby providing detailed multi-dimensional data support for the distribution of chemical components in complex samples.

[0003] Among them, the human body boron 10 analysis system is a special technical system used to evaluate the distribution characteristics of boron 10 isotopes in the human body. The system uses non-invasive imaging equipment and radioactively labeled boron compounds to measure the concentration and distribution data of boron 10 in human tissues. Through multi-dimensional image processing technology, the system can reconstruct the distribution of boron 10 in different tissue areas and form a detailed three-dimensional distribution image. The main purpose of this system is to provide real-time and accurate data support for the distribution of boron 10, optimize the application process of boron compounds in organisms, and enable scientists to better evaluate the dynamic change characteristics of boron 10.

[0004] Although the existing technology is highly sensitive and non-invasive, it has limitations in terms of spatial accuracy and feature capture of boron 10 isotope distribution. The decay frequency information is separated from the position information, resulting in insufficient correlation between the decay frequency and spatial position between voxels, causing partial data distortion and weakened explanatory power of the distribution. The current method lacks detailed processing in the concentration threshold comparison, and cannot accurately mark high- and low-frequency areas, affecting the feature recognition of low-content areas, making the description of the differentiated distribution of different concentration areas unclear, and easily causing confusion of information in the concentration feature map. The existing system relies on the presentation of the distribution of the overall image, and it is difficult to decompose the concentration difference and frequency characteristics item by item, resulting in the easy neglect of subtle differences in regional analysis, affecting the accuracy and interpretability of concentration feature analysis, and thus poor application effect in dealing with complex distribution environments. Summary of the invention

[0005] The purpose of the present invention is to solve the shortcomings of the prior art and to propose a human body boron 10 analysis system and method based on BNCT.

[0006] In order to achieve the above object, the present invention adopts the following technical scheme: The human body boron 10 analysis system based on BNCT comprises: The boron 10 distribution detection module detects the boron 10 isotope in the human body based on BNCT, obtains the voxel decay frequency set by capturing the decay frequency, position parameters and decay event time information of boron 10 in the voxel, and integrates the frequency data and position information of the voxel according to the voxel decay frequency set to obtain the initial distribution data of boron 10; The voxel reconstruction imaging module accumulates the decay frequencies in each voxel item by item based on the initial distribution data of boron 10, associates and locates the voxel frequency data through spatial coordinates, processes the spatial position and frequency value of the voxel, generates a boron 10 spatial distribution model, maps the concentration distribution value of the voxel position according to the boron 10 spatial distribution model, and obtains a voxel boron 10 concentration structure diagram; The concentration region division module compares the concentration threshold with the voxel based on the voxel boron 10 concentration structure map, extracts the frequency characteristics of the region, marks the high-frequency and low-frequency regions, and generates a characteristic concentration map of the region by combining the voxel frequency characteristics and concentration information; The concentration correlation analysis module calculates the difference value and frequency characteristics based on the characteristic concentration map of the area according to the concentration value difference of boron 10 in the area, calculates the concentration mean and frequency difference data in the area item by item, compares the concentration and distribution characteristics of the area in turn, and generates the boron concentration characteristic analysis results.

[0007] As a further solution of the present invention, the step of acquiring the voxel decay frequency set is specifically: Based on BNCT, the boron-10 isotope in the human body is detected, and the frequency parameters, position parameters and time parameters of the decay events of the boron-10 isotope are recorded. The frequency information of the decay events is sorted, and the voxels with potential analysis value are screened to establish the initial voxel frequency distribution set. The initial voxel frequency distribution set is calibrated to remove abnormal data that deviates from the distribution, and weighted correction is performed through position, time and frequency weight parameters, using the formula: , Generate voxel decay frequency calibration results, where represents the voxel decay frequency calibration result, represents the initial voxel decay frequency, and are the spatial and temporal parameters of the voxel, Weight parameter for controlling decay frequency correction; Combined with the distribution of the voxel decay frequency calibration results, the matching degree of the decay frequency distribution parameters is analyzed, and according to the peak position of the voxel decay, the decay frequency data that meets the standard distribution is determined to obtain the voxel decay frequency set.

[0008] As a further solution of the present invention, the steps for obtaining the initial distribution data of boron 10 are specifically as follows: Based on the voxel decay frequency set, calculating the decay frequency and position parameter combination of each voxel, calling the frequency data of each voxel, and generating a frequency and position correlation matrix; The frequency and position association matrix is ​​combined with the voxel position parameter, the attenuation frequency is selected and multiplied with the position parameter, and the frequency weight set is called, using the formula: , Get the frequency and position weight results of the voxels, where represents the frequency and position weight results of the voxels, represents the frequency weight parameter, represents the voxel frequency, represents the voxel position parameter; Based on the frequency and position weight results of the voxels, the high-frequency weight position data of the voxels are called, and the distribution of boron 10 is identified through density weighting to obtain the initial distribution data of boron 10.

[0009] As a further solution of the present invention, the steps for obtaining the spatial distribution model of boron 10 are specifically as follows: Based on the initial distribution data of boron 10, the decay frequency of each voxel is accumulated item by item, and the average frequency of the voxel frequency array is calculated according to the spatial position of the voxel to obtain the initial voxel frequency distribution result; Based on the initial voxel frequency distribution result, the voxel frequency array and the average frequency are used to adopt the formula: , The frequency deviation of the voxel is calculated, where represents the frequency deviation of the voxel, is the total number of voxels, is the decay frequency for each voxel, is the average frequency, is the frequency weight parameter, is the frequency adjustment factor; Combined with the frequency deviation of the voxel, the spatial coordinate set and the frequency deviation value are analyzed item by item, the frequency deviation value of each voxel is calibrated, the frequency deviation is associated with the voxel spatial position, the spatial coordinates and the frequency deviation are operated in correspondence to generate a boron 10 spatial distribution model.

[0010] As a further solution of the present invention, the steps for obtaining the voxel boron 10 concentration structure map are specifically as follows: Mapping the initial concentration value of the voxel position according to the boron 10 spatial distribution model, calculating and mapping the concentration distribution parameter of each voxel position to the concentration distribution value, performing mapping calculation on the voxel position, and generating an initial voxel position concentration map; Based on the initial voxel position concentration map, the concentration difference is refined and analyzed, and the concentration difference of adjacent voxels is calculated by difference weights, and the concentration difference formula is applied: , Calculate the difference between adjacent voxels and generate a concentration difference parameter matrix, where represents the difference between adjacent voxels, represents the concentration value at the voxel location, is the assigned concentration difference weight, is the mapping value of the voxel position, is the total number of pixels, is the average concentration distribution value; The concentration difference parameter matrix is ​​used as a judgment standard to eliminate and screen the voxel positions that do not meet the local concentration standard, and the voxel distribution is gradually screened and optimized to obtain a voxel boron 10 concentration structure map.

[0011] As a further solution of the present invention, the steps of obtaining the characteristic concentration map of the region are specifically as follows: Based on the voxel boron 10 concentration structure diagram, the concentration of each voxel is compared with a preset concentration threshold, the high-frequency and low-frequency voxel distribution areas are extracted, the voxel frequency characteristics of each area are calculated, and the area is marked to obtain the high-frequency and low-frequency area extraction results; Based on the high-frequency and low-frequency region extraction results, the concentration characteristic values ​​of the voxels in the high-frequency and low-frequency regions are calculated using the formula: , The concentration feature extraction results are obtained, where represents the concentration characteristic value, represents the voxel frequency in the high-frequency region, represents the voxel concentration threshold, is the frequency in the low frequency region, represents the difference in voxel distribution, is the adjustment factor; In combination with the concentration feature extraction result, the frequency feature and concentration feature data of the region are used for comparison and classification to generate a characteristic concentration map of the region.

[0012] As a further solution of the present invention, the steps for obtaining the boron concentration characteristic analysis results are specifically as follows: Based on the characteristic concentration map of the region, the concentration values ​​of boron 10 in each region are compared item by item, the concentration data of the difference value is extracted, the concentration difference in the region is recorded, the mean of the difference item is calculated, and the preliminary concentration difference data is obtained by comparing the frequency characteristics of each region; Based on the preliminary concentration difference data, the concentration mean and frequency difference of each area are calculated area by area, using the formula: , The concentration mean and frequency difference characteristic data are calculated, where: represents the concentration mean and frequency difference eigenvalues, is the concentration value of a voxel in the region, represents the average concentration value in the area. is the frequency value of the current area, is the adjustment factor for the mean concentration; Combined with the concentration mean and frequency difference characteristic data, the frequency difference distribution of voxels in the region is compared, the concentration mean and frequency difference characteristic value of the voxels in each region are calculated, the concentration and frequency characteristic data of the differentiated regions are compared item by item, and the boron concentration characteristic analysis results are obtained.

[0013] The BNCT-based human body boron 10 analysis method is performed based on the above-mentioned BNCT-based human body boron 10 analysis system, and includes the following steps: S1: Based on BNCT, the boron 10 isotope in the human body is detected, the decay frequency and position parameters of boron 10 in the voxel are captured, the time information of the decay event is recorded, the decay frequency data of the voxel is matched with the position parameter, the frequency data of the voxel is integrated with the spatial coordinates, each data is classified and sorted, and the initial distribution data of boron 10 is generated; S2: Based on the initial distribution data of boron 10, the decay frequencies of the voxels are accumulated item by item, the frequency information and the position are associated through the spatial coordinates, the frequency values ​​are aggregated and compared with the positions according to the positions, the frequency data of each voxel is bound to the spatial coordinates, and the spatial distribution model of boron 10 is generated; S3: Based on the B10 spatial distribution model, the B10 concentration parameters of the voxel positions are mapped, the spatial coordinates of the voxels are matched with the concentration values, and the frequency values ​​and coordinates of each position are associated item by item by summarizing the concentration distribution data of the voxels to generate a voxel B10 concentration structure diagram; S4: Based on the voxel boron 10 concentration structure diagram, the voxel concentration value is compared with the preset concentration threshold value item by item, the voxel frequency characteristics are divided, high-frequency and low-frequency voxels are extracted according to the frequency characteristics, the spatial distribution of the characteristic data is sorted out, and a regional characteristic concentration map is generated; S5: Based on the regional characteristic concentration map, the concentration value difference and frequency difference values ​​within the region are calculated, and the concentration mean and frequency changes of each region are compared item by item according to the regional division, and the concentration and distribution differences within the region are processed and summarized to generate the boron concentration characteristic analysis results.

[0014] Compared with the prior art, the advantages and positive effects of the present invention are: In the present invention, by capturing the decay frequency, position and time information of boron 10, the comprehensiveness and spatial accuracy of data acquisition are effectively improved, the spatial positioning processing of decay frequency accumulation is refined in the process, the frequency value in the voxel is associated with the spatial coordinate, and the structural accuracy of boron 10 concentration distribution is significantly improved, so that the subsequent concentration value mapping can directly reflect the actual concentration distribution status of each voxel position. In the concentration division process, the frequency characteristics are compared by concentration threshold, and the high and low frequency areas are clearly marked to ensure the accurate characterization of the frequency characteristics of different regions, optimize the concentration distribution analysis of boron 10 in each region, and generate a characteristic concentration map by fusing voxel frequency characteristics and concentration information to achieve the structured expression of regional characteristic concentration, laying the foundation for the subsequent concentration difference and frequency feature calculation, and refining the regional concentration and frequency difference data in a term-by-term operation mode to ensure the accurate analysis of concentration distribution characteristics, deepen the differentiated understanding of boron 10 distribution characteristics in different regions, and significantly improve the accuracy of boron 10 distribution detection in the human body and the depth of information interpretation through multi-dimensional difference analysis and regional marking technology, so as to provide data support for the application optimization of boron 10. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] Figure 1 is a system flow chart of the present invention; Figure 2 is a flow chart of a voxel decay frequency set in the present invention; Figure 3 This is a flow chart of the initial distribution data of boron 10 in the present invention; Figure 4 This is a flow chart of the spatial distribution model of boron 10 in the present invention; Figure 5 It is a flow chart of the voxel boron 10 concentration structure diagram in the present invention; Figure 6 A flow chart of a characteristic concentration map of a region in the present invention; Figure 7 Flow chart of the boron concentration characteristic analysis results in the present invention. DETAILED DESCRIPTION

[0016] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0017] In the description of the present invention, it should be understood that the terms "length", "width", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inside", "outside" and the like indicate positions or positional relationships based on the positions or positional relationships shown in the drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as limiting the present invention. In addition, in the description of the present invention, "multiple" means two or more, unless otherwise clearly and specifically defined.

[0018] See also Figure 1 , the present invention provides a technical solution: a human body boron 10 analysis system based on BNCT includes; The Boron 10 distribution detection module detects the Boron 10 isotope in the human body based on BNCT. By capturing the decay frequency, position parameters and decay event time information of Boron 10 in the voxel, the voxel decay frequency set is obtained. According to the voxel decay frequency set, the frequency data and position information of the voxel are integrated to obtain the initial distribution data of Boron 10. The voxel reconstruction imaging module accumulates the decay frequencies in each voxel item by item based on the initial distribution data of boron 10, associates and locates the voxel frequency data through spatial coordinates, processes the spatial position and frequency value of the voxel, generates a boron 10 spatial distribution model, maps the concentration distribution value of the voxel position according to the boron 10 spatial distribution model, and obtains the voxel boron 10 concentration structure diagram; The concentration region division module is based on the voxel boron 10 concentration structure map, compares the concentration threshold with the voxel, extracts the frequency characteristics of the region, marks the high-frequency and low-frequency regions, and combines the voxel frequency characteristics and concentration information to generate the characteristic concentration map of the region; The concentration correlation analysis module is based on the characteristic concentration map of the region. It calculates the difference value and frequency characteristics according to the concentration value difference of boron 10 in the region, calculates the concentration mean and frequency difference data in the region item by item, compares the concentration and distribution characteristics of the region in turn, and generates the boron concentration characteristic analysis results.

[0019] The voxel decay frequency set includes decay frequency, position parameters, and decay event time information. The boron 10 initial distribution data includes the voxel decay frequency set and position information. The boron 10 spatial distribution model includes voxel spatial coordinates, frequency data, and spatial distribution position. The voxel boron 10 concentration structure diagram includes the concentration distribution value of the voxel position, spatial coordinate mapping, and concentration distribution structure. The characteristic concentration diagram of the region includes high-frequency areas, low-frequency areas, frequency characteristics, and concentration information. The boron concentration characteristic analysis results include concentration mean, frequency difference data, and distribution characteristics.

[0020] See also Figure 2, the specific steps for obtaining the voxel decay frequency set are: Based on BNCT, the boron-10 isotope in the human body is detected, and the frequency parameters, position parameters and time parameters of the decay events of the boron-10 isotope are recorded. The frequency information of the decay events is sorted, and the voxels with potential analysis value are screened to establish the initial voxel frequency distribution set. Detect and record the frequency parameters, position parameters and time parameters of each event in the decay event of boron 10 isotope, establish a frequency set of decay events in the voxel according to the acquired decay event information, sort according to the frequency data of the decay events, classify the decay frequency values ​​according to the fluctuation range and mark them as voxels to be analyzed, classify the frequency data in each category, and identify the distribution position of high-frequency decay events by analyzing the change trend of the decay frequency value and the frequency peak position in the voxel, further eliminate the extreme frequency values ​​that are not in a reasonable range, and establish a data parameter set for the effective voxel area to provide preliminary screening data for subsequent steps. It is necessary to comprehensively consider the decay event frequency, position change and time information to ensure that the preliminary screening of voxels is more in line with the actual distribution, and establish an initial parameter set of voxel decay frequency for voxel frequency calibration in subsequent steps.

[0021] The initial voxel frequency distribution set is calibrated to remove abnormal data that deviates from the distribution, and weighted correction is performed through position, time and frequency weight parameters using the formula: , Generate voxel decay frequency calibration results, where represents the voxel decay frequency calibration result, represents the initial voxel decay frequency, and are the spatial and temporal parameters of the voxel, respectively. Weight parameter for controlling decay frequency correction; The benefit of the formula is that by introducing the collaborative correction method of spatial position Pa, time parameter Ta and weight parameter Wa, the accuracy of voxel decay frequency data is improved, the spatial and temporal resolution of the data is enhanced, and the frequency calibration is made more consistent with the actual decay characteristics; Set the initial voxel decay frequency data Fa to 5.4 Hz, the position parameter Pa to 3.2 mm, the time parameter Ta to 0.8 s, and the weight parameter Wa to 2.5. Substitute them into the formula to obtain the calibration frequency data The specific calculation process is as follows: , Expand the calculation to get: , The results show that the calibrated voxel decay frequency data is 3.69 Hz, which is used to represent the corrected frequency value and lays a data foundation for subsequent frequency data distribution analysis.

[0022] Combined with the distribution of the voxel decay frequency calibration results, the matching degree of the decay frequency distribution parameters is analyzed, and according to the peak position of the voxel decay, the decay frequency data that meets the standard distribution is determined to obtain the voxel decay frequency set; Combined with the calibration voxel decay frequency data generated in the previous step, distribution analysis is performed according to the variation range of the spatial position and time parameters of the calibrated data. The calibrated frequency data is classified based on the numerical range of Ff and its rationality is further judged. The reference range of the distribution is set, and the voxel decay frequency data set is compared with other distributions in turn. The values ​​that deviate from the standard distribution are identified and eliminated. The decay frequency peaks that meet the characteristics of the standard distribution are extracted, and the corresponding peak voxel positions are used to construct the final decay frequency distribution results. In this process, the peak positions in the distribution results need to be compared for rationality to ensure that the control frequency is typical and meets the detection requirements, and finally the voxel decay frequency data that meets the expected standards are obtained.

[0023] See also Figure 3 The specific steps for obtaining the initial distribution data of boron 10 are as follows: Based on the voxel decay frequency set, the decay frequency and position parameter combination of each voxel is calculated, the frequency data of each voxel is called, and a frequency and position correlation matrix is ​​generated; Based on the voxel decay frequency set, the voxel decay frequency is first quantified. The frequency and spatial position data of each voxel are combined, and the distribution characteristics of each voxel frequency are extracted. The set operation of decay characteristics and spatial positions is calculated, and each frequency result is combined with its corresponding voxel position parameter. The distribution of frequency data is distinguished by frequency amplitude and the position is marked to form the frequency amplitude information corresponding to the voxel position. The frequency information is integrated with the position coordinates to obtain the frequency and position information integration matrix.

[0024] Combine the frequency and position association matrix with the voxel position parameters, select the attenuation frequency and position parameters for product operation, call the frequency weight set, and use the formula: , Get the frequency and position weight results of the voxels, where represents the frequency and position weight results of the voxels, represents the frequency weight parameter, represents the voxel frequency, represents the voxel position parameter; The benefit of the formula is that by introducing absolute value and square root operations, the calculation of frequency-position weights is smoother, while the stability and flexibility of the results are improved; In the formula, the frequency weight parameter is the weighted value of the voxel frequency extracted from the frequency-position association matrix, and the position parameter is the position coordinate of each voxel, the frequency It is obtained by sampling the frequency of each voxel; The calculation is performed using the following specific values: =2.5, =1.8, =3.2; Substitute into the formula: , , , ; The result shows that the frequency-position weight is 5.37, which will be used in subsequent voxel high-frequency integration processing to obtain the initial distribution data of boron 10.

[0025] Based on the frequency and position weight results of the voxels, the high-frequency weight position data of the voxels are called, and the distribution of boron 10 is identified by density weighting to obtain the initial distribution data of boron 10; Based on the frequency-position weight result, the frequency and position parameters of the voxel are called, and the density-weighted calculation method is selected. The frequency value and position data of the voxel are combined one by one. The frequency proportion of the voxel is determined through successive iterations, and the voxel set that best conforms to the density distribution is obtained. The density-weighted synthesis is performed on the frequency parameter and position parameter of each voxel. After the frequency proportion parameters in the comprehensive density distribution are combined, the voxel cluster with high frequency distribution is further screened out, and the data of the cluster is summarized to obtain the initial distribution data of boron 10.

[0026] See also Figure 4 The specific steps for obtaining the spatial distribution model of boron 10 are as follows: Based on the initial distribution data of boron 10, the decay frequency of each voxel is accumulated item by item, and the average frequency of the voxel frequency array is calculated according to the spatial position of the voxel to obtain the initial voxel frequency distribution result; Based on the initial distribution data of boron 10, the decay frequency of each voxel is accumulated item by item to generate a decay frequency array containing all voxel frequencies. The elements in the array represent the frequency of each voxel. The data in the decay frequency array are numbered according to the spatial position of the voxel to achieve accurate correspondence between the frequency value and the spatial position of the specific voxel. The generated frequency array contains complete spatial frequency information, which is sorted according to a specific spatial position to form an initial voxel frequency distribution. The overall average value of the initial frequency distribution is calculated based on the decay frequency data, which represents the statistical mean of all voxel frequencies. Through the association of this spatial position and frequency data, the decay frequency of each voxel is linked to the initial distribution of boron 10 in the area, and then the initial voxel frequency distribution result is obtained, providing basic data support for the subsequent calculation of frequency deviation.

[0027] Based on the initial voxel frequency distribution results, the voxel frequency array and the average frequency are used, using the formula: , The frequency deviation of the voxel is calculated, where represents the frequency deviation of the voxel, is the total number of voxels, is the decay frequency for each voxel, is the average frequency, is the frequency weight parameter, is the frequency adjustment factor; The formula is beneficial in that it introduces adjustment coefficients and weight parameters based on the original voxel frequency deviation calculation, thereby improving the dynamic adjustment ability of the frequency deviation result, making the frequency deviation calculation adaptable between different voxels; For ease of calculation, the number of voxels is =5, the decay frequency of each voxel They are 2, 3, 4, 5, and 6 respectively, and the average frequency Set the frequency weight to 4.0 =0.8 and adjustment coefficient =0.2, the specific calculation process of the formula is as follows: Calculate the adjustment result for each frequency: , ; Compute the sum of squares and take the square root: , Expand and calculate: , This result shows that the voxel frequency deviation in the spatial distribution of boron 10 is 2.33, which represents the decay frequency deviation level in this area and is used for deviation data calling in subsequent distribution model construction.

[0028] Combined with the frequency deviation of voxels, the frequency deviation value of each voxel is calibrated by analyzing the correspondence between the spatial coordinate set and the frequency deviation value, the frequency deviation is associated with the spatial position of the voxel, and the spatial coordinates are operated in correspondence with the frequency deviation to generate a spatial distribution model of boron 10; Combined with the initial distribution data of boron 10, the frequency deviation Δ of each voxel is used, and the voxel position in the spatial coordinates is used to correspond to the associated frequency deviation value Δ to form a bidirectional mapping relationship between the voxel coordinates and the frequency deviation, so as to establish the overall boron 10 spatial distribution model. The model generation process includes quantifying the frequency deviation value Δ of each voxel and mapping it to the corresponding spatial coordinate position, and at the same time establishing an overall distribution structure for the frequency deviation value Δ of all voxels. Through this model, the distribution information of boron 10 at different spatial positions can be accurately presented, so that the decay frequency of each voxel can be fully reflected in the entire spatial distribution model.

[0029] See also Figure 5 The specific steps for obtaining the voxel boron 10 concentration structure map are as follows: Mapping the initial concentration value of the voxel position according to the boron 10 spatial distribution model, calculating and mapping the concentration distribution parameter of each voxel position to the concentration distribution value, performing mapping calculation on the voxel position, and generating an initial voxel position concentration map; Based on the boron 10 spatial distribution model, the voxel position concentration data is mapped, and the concentration value is calculated for each voxel position through the concentration distribution parameters. In the calculation process, the voxel position coordinate value is used as input, and the corresponding concentration distribution function is called. The concentration value of each voxel is mapped according to the voxel position and distribution model. The concentration distribution value of each voxel is corrected based on the parameter structure of the model to form an initial concentration distribution map. The concentration value is then checked for each voxel position. The concentration values ​​that are obviously deviated from the distribution are eliminated through specific concentration matching rules to ensure the numerical integrity and accuracy of the initial voxel position concentration map, and generate the initial voxel position concentration map.

[0030] Based on the initial voxel position concentration map, the concentration difference is refined and analyzed, and the concentration difference of adjacent voxels is calculated by difference weights, and the concentration difference formula is applied: , Calculate the difference between adjacent voxels and generate a concentration difference parameter matrix, where represents the difference between adjacent voxels, represents the concentration value at the voxel location, is the assigned concentration difference weight, is the mapping value of the voxel position, is the total number of pixels, is the average concentration distribution value; The formula is useful because it converts the concentration value at the voxel location into With distribution weight Combined with the average difference offset, the concentration difference correction is performed to improve the adaptability and accuracy of the concentration difference calculation; For each parameter in the formula, Represents the concentration value at each voxel position, obtained by model distribution mapping calculation, The correction weight representing the concentration distribution difference needs to be calculated and determined based on the concentration difference fluctuation range between each voxel and its surrounding voxels. is the concentration map value at the voxel position, which is directly generated by the spatial distribution model. is the total number of pixels, obtained through data statistics, is the average concentration distribution value, which is obtained by averaging the concentration values ​​of all voxels in the initial mapping. Let’s take a numerical example: set up is 0.8, is 1.2, is 0.6, is 0.5, is 1000, the calculation is as follows: ; This result shows that Pd is the relative concentration difference value of each voxel, which can be used in the subsequent voxel concentration difference structure optimization.

[0031] The concentration difference parameter matrix is ​​used as the judgment standard to eliminate and screen the voxel positions that do not meet the local concentration standard, and the voxel distribution is gradually screened and optimized to obtain the voxel boron 10 concentration structure map; Using the concentration difference parameter matrix as the judgment criterion, the concentration difference of voxel positions is gradually screened, and conditional judgment is made on the concentration difference between each voxel position and the adjacent voxel. By setting the concentration difference threshold, the voxel positions that do not meet the requirements are gradually eliminated, and the voxel positions with difference values ​​lower than the threshold are retained. If the concentration difference exceeds the set upper limit, the voxel concentration is adjusted. By correcting and screening one by one, a set of voxel positions that meet the standards is formed, and a voxel boron 10 concentration structure diagram is generated.

[0032] See also Figure 6 , the steps for obtaining the characteristic concentration map of the region are as follows: Based on the voxel boron 10 concentration structure diagram, the concentration of each voxel is compared with the preset concentration threshold, the high-frequency and low-frequency voxel distribution areas are extracted, the voxel frequency characteristics of each area are calculated, and the area is marked to obtain the high-frequency and low-frequency area extraction results; Based on the voxel boron 10 concentration structure diagram, the concentration of each voxel is detected in turn and compared with the preset concentration threshold one by one. The voxel distribution information above or below the threshold is detected voxel by voxel. During the comparison process, the voxel concentration is higher or lower than the threshold, and it is marked as a high-frequency or low-frequency voxel area, forming a preliminary high-frequency and low-frequency voxel area classification. The frequency characteristics of each voxel are counted separately, and the voxel frequency characteristics of each area are recorded and calculated to determine the high-frequency or low-frequency attributes of each voxel area, and obtain the high-frequency and low-frequency area extraction results.

[0033] Based on the extraction results of high-frequency and low-frequency regions, the concentration characteristic values ​​of voxels in the high-frequency and low-frequency regions are calculated using the formula: , The concentration feature extraction results are obtained, where represents the concentration characteristic value, represents the voxel frequency in the high-frequency region, represents the voxel concentration threshold, is the frequency in the low frequency region, represents the difference in voxel distribution, is the adjustment factor; The benefit of the formula is that it uses frequency characteristics and adjustment coefficients to effectively extract concentration characteristic values, achieving the accuracy of concentration differences between high-frequency and low-frequency areas; The voxel frequency of the high-frequency area can be obtained by detecting the ratio of the number of voxels above the threshold to the total number of voxels in the area; is the voxel concentration threshold, which is directly introduced into the calculation according to the preset threshold and is set to 40; The frequency of the low-frequency region was calculated by the ratio of the number of voxels below the threshold to the total number of voxels and was set to 0.2; is the voxel distribution difference, which is calculated based on the standard deviation of the voxel distribution in the region and is set to 5; The adjustment coefficient is used to adjust the concentration difference and is dynamically set to 2 according to the fluctuation range of frequency characteristics and distribution differences; Specific operation process: Substitute into numerical calculations: ; Calculate the absolute difference: ; The final result is: ; The result shows that the concentration characteristic value is 33.13, which is the key basis for subsequent concentration feature extraction.

[0034] Combined with the concentration feature extraction results, the frequency characteristics and concentration feature data of the region are used for comparison and classification to generate a characteristic concentration map of the region; Combined with the concentration feature extraction results, the frequency features and concentration feature data are compared and classified, the high-frequency and low-frequency concentration feature values ​​are input into the comparison system, and the comparison operation is performed region by region according to the distribution of characteristic concentration in each region. The high- and low-frequency concentration division of the region is further verified through the concentration distribution of each voxel in the region, so as to confirm the characteristic concentration information of each region, mark and classify the high-frequency and low-frequency voxel concentrations in different regions, and generate the characteristic concentration map of the region.

[0035] See also Figure 7 , the specific steps for obtaining the boron concentration characteristic analysis results are: Based on the characteristic concentration map of the region, the concentration values ​​of boron 10 in each region were compared item by item, the concentration data of the difference values ​​were extracted, the concentration differences in the regions were recorded, the mean of the difference items was calculated, and the preliminary concentration difference data were obtained by comparing the frequency characteristics of each region; Based on the characteristic concentration map of the region, the boron 10 concentration value in each region is detected one by one. The voxels with larger difference values ​​are extracted from each region through the difference detection of concentration values, and a concentration difference data set for each region is constructed to facilitate separate analysis of high difference values ​​in subsequent comparisons. In the specific operation, the boron 10 concentration value of each voxel is extracted in turn, and the difference between it and the mean concentration in the region is calculated to ensure the accuracy of the difference value calculation. The concentration difference data are recorded and summarized one by one. By comparing the concentration difference values, the regions with obvious frequency characteristics are further screened out, and the high-frequency or low-frequency voxels are marked to form preliminary concentration difference data. By summarizing the concentration difference characteristic data of each region, the overall difference in the concentration characteristics of the voxels in each region is analyzed to establish the concentration distribution structure of each region, lay the foundation for further frequency characteristic comparison, obtain preliminary concentration difference data, and lay the data foundation for the next step of frequency characteristic analysis.

[0036] Based on the preliminary concentration difference data, the concentration mean and frequency difference of each area were calculated area by area using the formula: , The concentration mean and frequency difference characteristic data are calculated, where: represents the concentration mean and frequency difference eigenvalues, is the concentration value of a voxel in the region, represents the average concentration value in the area. is the frequency value of the current area, is the adjustment factor for the mean concentration; The benefit of the formula is that it introduces concentration difference and frequency characteristic adjustment coefficients during calculation, so that concentration difference and frequency characteristics are comprehensively balanced and the key characteristics of frequency and concentration differences within the region are effectively extracted; It represents the concentration value of a single voxel detected in the area. The actual measured value is 45, which is obtained by voxel-by-voxel concentration detection. is the mean value of the voxel concentration in the region, which is calculated based on the average value of the concentration values ​​of all voxels in the region, and the value is set to 42; is the voxel frequency value in the region, which is 1.2 according to the ratio of the number of voxels in each region to the total number of voxels; is the concentration mean adjustment coefficient, which is obtained by calculating the standard deviation of the concentration difference of voxels in each region and is set to 2; The specific operation process is as follows: Calculate the concentration difference term: ; Calculate the frequency difference term: ; Combine the two results: ; The results show that the concentration frequency characteristic value is 7.77, which reflects the differences in concentration differences and frequency distribution characteristics within the region, and is used as the input characteristic value for subsequent boron concentration characteristic analysis.

[0037] Combine the concentration mean and frequency difference characteristic data, compare the frequency difference distribution of voxels in the region, calculate the concentration mean and frequency difference characteristic value of the voxels in each region, compare the concentration and frequency characteristic data of the differentiated regions item by item, and obtain the boron concentration characteristic analysis results; Combined with the concentration-frequency characteristic results, the differences in regional concentration and frequency distribution characteristics are further analyzed. First, the frequency difference characteristic data is input into the analysis system to screen out regional frequency difference values ​​that meet the analysis conditions, and the concentration distribution of each region is compared and the concentration frequency characteristics of each voxel are marked to ensure the comprehensiveness of the results. In the specific processing, the concentration distribution of the high-frequency area and the low-frequency area are extracted and summarized respectively by calling the concentration frequency characteristic data, and the concentration characteristics of the high-frequency and low-frequency areas are marked according to the frequency difference to generate the boron concentration characteristic analysis results.

[0038] The BNCT-based human body boron 10 analysis method is performed based on the above-mentioned BNCT-based human body boron 10 analysis system, and includes the following steps: S1: Based on BNCT, the boron 10 isotope in the human body is detected, the decay frequency and position parameters of boron 10 in the voxel are captured, the time information of the decay event is recorded, the decay frequency data of the voxel is matched with the position parameter, the frequency data of the voxel is integrated with the spatial coordinates, each data is classified and sorted, and the initial distribution data of boron 10 is generated; S2: Based on the initial distribution data of boron 10, the decay frequencies of the voxels are accumulated item by item, the frequency information and the position are associated through the spatial coordinates, the frequency values ​​are aggregated and compared with the positions by position, the frequency data of each voxel is bound to the spatial coordinates, and the spatial distribution model of boron 10 is generated; S3: Based on the B10 spatial distribution model, the B10 concentration parameters of the voxel position are mapped, the spatial coordinates of the voxel are matched with the concentration value, and the frequency value and coordinate of each position are associated item by item by summarizing the concentration distribution data of the voxel to generate the B10 concentration structure diagram of the voxel; S4: Based on the voxel boron 10 concentration structure map, the voxel concentration value is compared with the preset concentration threshold item by item, the voxel frequency characteristics are divided, high-frequency and low-frequency voxels are extracted according to the frequency characteristics, the spatial distribution of the feature data is sorted out, and a regional feature concentration map is generated; S5: Based on the regional characteristic concentration map, calculate the concentration value difference and frequency difference value within the region, compare the concentration mean and frequency changes of each region item by item according to the regional division, process and summarize the concentration and distribution differences within the region, and generate the boron concentration characteristic analysis results.

[0039] The above are only preferred embodiments of the present invention and are not intended to limit the present invention in other forms. Any technician familiar with the profession may use the technical contents disclosed above to change or modify them into equivalent embodiments with equivalent changes and apply them to other fields. However, any simple modification, equivalent change and modification made to the above embodiments based on the technical essence of the present invention without departing from the technical solution of the present invention still falls within the protection scope of the technical solution of the present invention.

Claims

1. The human body boron 10 analysis system based on BNCT is characterized by: The system comprises: The boron 10 distribution detection module detects the boron 10 isotope in the human body based on BNCT, obtains the voxel decay frequency set by capturing the decay frequency, position parameters and decay event time information of boron 10 in the voxel, and integrates the frequency data and position information of the voxel according to the voxel decay frequency set to obtain the initial distribution data of boron 10; The voxel reconstruction imaging module accumulates the decay frequencies in each voxel item by item based on the initial distribution data of boron 10, associates and locates the voxel frequency data through spatial coordinates, processes the spatial position and frequency value of the voxel, generates a boron 10 spatial distribution model, maps the concentration distribution value of the voxel position according to the boron 10 spatial distribution model, and obtains a voxel boron 10 concentration structure diagram; The concentration region division module compares the concentration threshold with the voxel based on the voxel boron 10 concentration structure map, extracts the frequency characteristics of the region, marks the high-frequency and low-frequency regions, and generates a characteristic concentration map of the region by combining the voxel frequency characteristics and concentration information; The concentration correlation analysis module calculates the difference value and frequency characteristics based on the characteristic concentration map of the area according to the concentration value difference of boron 10 in the area, calculates the concentration mean and frequency difference data in the area item by item, compares the concentration and distribution characteristics of the area in turn, and generates the boron concentration characteristic analysis results.

2. The human body boron 10 analysis system based on BNCT according to claim 1, characterized in that: The steps for obtaining the voxel decay frequency set are specifically as follows: Based on BNCT, the boron-10 isotope in the human body is detected, and the frequency parameters, position parameters and time parameters of the decay events of the boron-10 isotope are recorded. The frequency information of the decay events is sorted, and the voxels with potential analysis value are screened to establish the initial voxel frequency distribution set. The initial voxel frequency distribution set is calibrated to remove abnormal data that deviates from the distribution, and weighted correction is performed through position, time and frequency weight parameters, using the formula: ， Generate voxel decay frequency calibration results, where represents the voxel decay frequency calibration result, represents the initial voxel decay frequency, and are the spatial and temporal parameters of the voxel, respectively. Weight parameter for controlling decay frequency correction; Combined with the distribution of the voxel decay frequency calibration results, the matching degree of the decay frequency distribution parameters is analyzed, and according to the peak position of the voxel decay, the decay frequency data that meets the standard distribution is determined to obtain the voxel decay frequency set.

3. The human body boron 10 analysis system based on BNCT according to claim 2, characterized in that: The steps for obtaining the initial distribution data of boron 10 are specifically as follows: Based on the voxel decay frequency set, calculating the decay frequency and position parameter combination of each voxel, calling the frequency data of each voxel, and generating a frequency and position correlation matrix; The frequency and position association matrix is ​​combined with the voxel position parameter, the attenuation frequency is selected and multiplied with the position parameter, and the frequency weight set is called, using the formula: ， Get the frequency and position weight results of the voxels, where represents the frequency and position weight results of the voxels, represents the frequency weight parameter, represents the voxel frequency, represents the voxel position parameter; Based on the frequency and position weight results of the voxels, the high-frequency weight position data of the voxels are called, and the distribution of boron 10 is identified through density weighting to obtain the initial distribution data of boron 10.

4. The BNCT-based human body boron 10 analysis system according to claim 3, characterized in that: The steps for obtaining the boron 10 spatial distribution model are specifically as follows: Based on the initial distribution data of boron 10, the decay frequency of each voxel is accumulated item by item, and the average frequency of the voxel frequency array is calculated according to the spatial position of the voxel to obtain the initial voxel frequency distribution result; Based on the initial voxel frequency distribution result, the voxel frequency array and the average frequency are used to adopt the formula: ， The frequency deviation of the voxel is calculated, where represents the frequency deviation of the voxel, is the total number of voxels, is the decay frequency for each voxel, is the average frequency, is the frequency weight parameter, is the frequency adjustment factor; Combined with the frequency deviation of the voxel, the spatial coordinate set and the frequency deviation value are analyzed item by item, the frequency deviation value of each voxel is calibrated, the frequency deviation is associated with the voxel spatial position, the spatial coordinates and the frequency deviation are operated in correspondence to generate a boron 10 spatial distribution model.

5. The BNCT-based human body boron 10 analysis system according to claim 4, characterized in that: The steps for obtaining the voxel boron 10 concentration structure diagram are specifically as follows: Mapping the initial concentration value of the voxel position according to the boron 10 spatial distribution model, calculating and mapping the concentration distribution parameter of each voxel position to the concentration distribution value, performing mapping calculation on the voxel position, and generating an initial voxel position concentration map; Based on the initial voxel position concentration map, the concentration difference is refined and analyzed, and the concentration difference of adjacent voxels is calculated by difference weights, and the concentration difference formula is applied: ， Calculate the difference between adjacent voxels and generate a concentration difference parameter matrix, where represents the difference between adjacent voxels, represents the concentration value at the voxel location, is the assigned concentration difference weight, is the mapping value of the voxel position, is the total number of pixels, is the average concentration distribution value; The concentration difference parameter matrix is ​​used as a judgment standard to eliminate and screen the voxel positions that do not meet the local concentration standard, and the voxel distribution is gradually screened and optimized to obtain a voxel boron 10 concentration structure map.

6. The BNCT-based human body boron 10 analysis system according to claim 5, characterized in that: The steps for obtaining the characteristic concentration map of the region are specifically as follows: Based on the voxel boron 10 concentration structure diagram, the concentration of each voxel is compared with a preset concentration threshold, the high-frequency and low-frequency voxel distribution areas are extracted, the voxel frequency characteristics of each area are calculated, and the area is marked to obtain the high-frequency and low-frequency area extraction results; Based on the extraction results of the high-frequency and low-frequency regions, the concentration characteristic values ​​of the voxels in the high-frequency and low-frequency regions are calculated using the formula: ， The concentration feature extraction results are obtained, where represents the characteristic value of concentration, represents the voxel frequency in the high-frequency region, represents the voxel concentration threshold, is the frequency in the low frequency region, represents the difference in voxel distribution, is the adjustment factor; In combination with the concentration feature extraction result, the frequency feature and concentration feature data of the region are used for comparison and classification to generate a characteristic concentration map of the region.

7. The BNCT-based human body boron 10 analysis system according to claim 6, characterized in that: The steps for obtaining the boron concentration characteristic analysis results are specifically as follows: Based on the characteristic concentration map of the region, the concentration values ​​of boron 10 in each region are compared item by item, the concentration data of the difference value is extracted, the concentration difference in the region is recorded, the mean of the difference item is calculated, and the preliminary concentration difference data is obtained by comparing the frequency characteristics of each region; Based on the preliminary concentration difference data, the concentration mean and frequency difference of each area are calculated area by area, using the formula: ， The concentration mean and frequency difference characteristic data are calculated, where: represents the concentration mean and frequency difference eigenvalues, is the concentration value of a voxel in the region, represents the average concentration value in the area. is the frequency value of the current area, is the adjustment factor for the mean concentration; Combined with the concentration mean and frequency difference characteristic data, the frequency difference distribution of voxels in the region is compared, the concentration mean and frequency difference characteristic value of the voxels in each region are calculated, the concentration and frequency characteristic data of the differentiated regions are compared item by item, and the boron concentration characteristic analysis results are obtained.

8. A method for analyzing boron 10 in the human body based on BNCT, characterized in that: The method according to any one of claims 1 to 7 is performed by the BNCT-based human body boron 10 analysis system, comprising the following steps: Based on BNCT, the boron 10 isotope in the human body is detected, the decay frequency and position parameters of boron 10 in the voxel are captured, the time information of the decay event is recorded, the decay frequency data of the voxel is matched with the position parameter, the frequency data of the voxel is integrated with the spatial coordinates, each data is classified and sorted, and the initial distribution data of boron 10 is generated; Based on the initial distribution data of boron 10, the decay frequencies of voxels are accumulated item by item, the frequency information and the position are associated through the spatial coordinates, the frequency values ​​are aggregated and compared with the positions according to the positions, the frequency data of each voxel is bound to the spatial coordinates, and the spatial distribution model of boron 10 is generated; Based on the B10 spatial distribution model, the B10 concentration parameters of the voxel positions are mapped, the spatial coordinates of the voxels are matched with the concentration values, and the frequency values ​​and coordinates of each position are associated item by item by summarizing the concentration distribution data of the voxels to generate a Voxel B10 concentration structure diagram; Based on the voxel boron 10 concentration structure diagram, the voxel concentration value is compared with the preset concentration threshold value item by item, the voxel frequency characteristics are divided, high-frequency and low-frequency voxels are extracted according to the frequency characteristics, the spatial distribution of the characteristic data is sorted out, and a regional characteristic concentration map is generated; Based on the regional characteristic concentration map, the concentration value difference and frequency difference values ​​within the region are calculated, the concentration mean and frequency changes of each region are compared item by item according to the regional division, the concentration and distribution differences within the region are processed and summarized, and the boron concentration characteristic analysis results are generated.

Images