Marine crab radiation dose estimation method and device based on voxel model

Abstract

The invention discloses a marine crab radiation dose estimation method based on a voxel model, which obtains detailed anatomical data of marine crabs by using a CT scanning technology, constructs an accurate voxel model, and can effectively distinguish organs such as crab shells, crab meat, hearts, gills and crab spawns. The process avoids the defect that a traditional simplified model neglects organ differences, so that the accuracy of radiation dose estimation is improved. According to the method, the radiation absorption characteristic of each organ can be accurately simulated, and a more real radiation dose calculation result is provided. The three-dimensional anatomical data of the living body can be accurately obtained based on the CT and image processing technology, so that the technical scheme has high adaptability, is suitable for marine crabs, can be flexibly adjusted according to anatomical features of different species, and can be widely applied to the fields of marine crabs, marine crabs and marine crabs. The method is widely applied to radiation dose estimation of other aquatic organisms or land species.

Description

Technical Field

[0001] The present invention relates to the field of radiation estimation of non-human species, and in particular to a method, device and computer-readable storage medium for estimating radiation dose of marine crabs based on a voxel model. Background Art

[0002] As society pays more and more attention to the protection of non-human species, the possible impact of ionizing radiation on non-human species has gradually become the focus of global attention. International organizations such as the International Atomic Energy Agency (IAEA), the United Nations Commission on the Effects of Radiation (UNSCEAR), and the International Commission on Radiological Protection (ICRP) have successively published relevant publications to discuss and study the radiation effects of ionizing radiation on non-human species and dose estimation issues. In particular, in my country, the "Radioactive Pollution Prevention and Control Law" (issued in 2003) and the "Nuclear Safety Law" (implemented in 2018) have made clear provisions for the protection of non-human species, emphasizing the importance of evaluating the impact of radiation on these species.

[0003] Internationally, there are currently a variety of assessment methods and models for radiation doses and radiation effects of non-human species, such as Canada's AECL approach, the UK's R&D128, the Netherlands' LAKECOB, the EU's ERICA, the US's RESRAD-BIOTA, France's EDEN2, CATEAUR, etc. These models classify organisms in different habitats (such as freshwater, marine, and terrestrial), and establish radiation dose assessment models based on biological morphology for each type of organism. Taking the ERICA program as an example, it divides organisms into freshwater, marine, and terrestrial organisms according to their habitats. Each type is further divided according to the family and growth stage of the plant, with a total of 40 types of biological models. The RESRAD-BIOTA program divides organisms into four major categories: aquatic animals, coastal animals, terrestrial animals, and terrestrial plants. The biological models are divided by body size, and all use ellipsoids as representative models.

[0004] However, these existing models have significant shortcomings:

[0005] Simplification of biological models: Existing models such as ERICA and RESRAD-BIOTA use an ellipsoid as the reference biological morphology, assuming that the density, elemental composition and radionuclide distribution in the organism are uniform throughout the body. Although this simplified model can estimate radiation dose to a certain extent, it fails to accurately reflect the actual morphology and complex internal structure of the organism, especially the radiation absorption characteristics of different tissues and organs.

[0006] Lack of consideration of radiation differences between organs: Traditional ellipsoid models ignore the density differences, elemental composition and radionuclide distribution of sensitive organs in organisms. For example, there are significant differences in the radiation absorption characteristics of organs such as crab shells, crab meat, heart, gills and crab roe of marine crabs, but these differences are not reflected in existing models. Therefore, there are uncertainties and errors in the dose values ​​calculated using these simplified models, making it difficult to accurately assess the actual impact of radiation on non-human species.

[0007] In recent years, with the rapid development of tomography technology (such as computed tomography CT) and computer image processing technology, voxel models have emerged and become an important technical means to solve the defects of traditional models. Voxel models are based on CT scanning, magnetic resonance imaging (MRI) or frozen sections, and use the acquired image sets to reconstruct three-dimensional biological models. Compared with traditional simplified geometric models, voxel models can more accurately display the various tissues and organs in the body, and take into account factors such as their specific morphology, size, and density.

[0008] Although voxel models have significant advantages in radiation dose estimation, existing voxel modeling technologies still have problems such as difficulty in obtaining biological anatomical data, complexity in the processing process, and insufficient accuracy in organ segmentation. For more complex organisms such as marine crabs, the construction of their voxel models often faces the following technical difficulties:

[0009] Incomplete data acquisition: The cross-sectional image data of marine crabs may be noisy or have insufficient resolution, making it difficult to accurately identify and separate individual tissues during organ segmentation.

[0010] Tissue and organ identification and segmentation are difficult: Traditional image segmentation technology cannot accurately distinguish complex organs, especially between organs with large differences in density and elemental composition (such as between crab shell and crab meat, heart and gills). Summary of the invention

[0011] The purpose of the present invention is to solve the problem that simplified models in the prior art ignore the differences in organs in the organism, and to provide a more accurate and reliable method and device for estimating radiation dose of marine crabs based on a voxel model. At the same time, a computer-readable storage medium capable of executing the method is provided.

[0012] In order to solve the above technical problems, the technical method adopted by the present invention is as follows: The present invention discloses a method for estimating the radiation dose of marine crabs based on a voxel model, comprising the following steps:

[0013] Step S1. Acquisition of tomographic sequence images of organisms: performing tomographic scanning on marine crabs to obtain tomographic sequence images of organisms;

[0014] Step S2. Processing of tomographic sequence images: using image processing software to process the tomographic images, segmenting the organs / tissues of the organism, and assigning different RGB values ​​to each tissue and organ for labeling;

[0015] Step S3. Voxel matrix acquisition: converting the processed image into a Monte Carlo repeating structure card input matrix according to the RGB value of each color to obtain a voxel matrix for Monte Carlo simulation, specifically including the following steps:

[0016] S31. Image size definition: Define the length and width of each CT image to ensure consistent image size;

[0017] S32. Pixel reading and classification: setting the category of the current pixel to identify which tissue or organ the pixel belongs to;

[0018] Initialize the current pixel category counter currentNum to 0, which is used to record the number of pixels of each category;

[0019] S33. Read each pixel in turn and obtain the RGB value: read each pixel of the image and obtain the RGB value of the pixel; the RGB value determines which organ or tissue category the pixel belongs to. According to the RGB value, determine which tissue / organ category the pixel belongs to, and set currentClass to the number of the category;

[0020] S33. Determine whether the next pixel is a pixel of the current category and update the counter: Determine whether the next pixel belongs to the current category currentClass; if it is a pixel of the same category, add 1 to the counter currentNum and continue counting; if the pixel does not belong to the current category currentClass, store the current category and its number in the result list res, update currentClass to the new pixel category, and then reset the counter currentNum to 1;

[0021] S34. Output the results and generate a voxel matrix: All classification results including the number of pixels in each category and their types are continuously output. The number of pixels in each category currentNum and the tissue category currentClass to which it belongs will be used as input data of the voxel matrix. The voxel matrix will be used for subsequent Monte Carlo simulation. Ultimately, the distribution, size and morphology of all tissues and organs are accurately recorded in the voxel matrix.

[0022] Step S4. Calculate the deposition energy of organisms using MCNP simulation: input the voxel matrix into the MCNP software, and calculate the deposition energy of radionuclides on organisms through Monte Carlo simulation; according to the deposition energy, use the internal and external radiation dose calculation formula to calculate the internal and external radiation dose coefficients of marine crabs, and further calculate the dose rate.

[0023] Furthermore, the tissues or organs include crab shell, crab meat, heart, gills and crab roe.

[0024] Further, the organ or tissue that concentrates a certain radionuclide is set as the source organ / tissue; then, the organ or tissue that is irradiated is set as the target organ / tissue; finally, the energy in the source organ / tissue k is E k The calculation formula of the internal radiation dose coefficient DC of the target organ / tissue j is:

[0025] In the formula, The energy in the source organ / tissue k is E k The radiation corresponds to the internal dose coefficient of the target organ / tissue j, in units of (μGy·h -1 ) / (Bq·kg -1 );Y k The energy of each decay of the radionuclide is E k The radiation yield; M k is the mass of the source organ / tissue in kg; F k,j (E) is the energy of the radionuclide decay once, E k The fraction of radiation energy absorbed by the organism;

[0026] M j is the mass of the target organ / tissue in kg;

[0027] The external radiation dose coefficient is calculated by treating all environmental media in which the organism is located as radiation sources. The calculation formula is:

[0028]

[0029] In the formula, The mass is M k The energy in the external environment medium is E k The radiation corresponds to the external radiation dose coefficient of the target organ / tissue j, in units of (μGy·h -1 ) / (Bq·kg -1 );

[0030] Y k The energy of each decay of the radionuclide is E k The radiation yield; F k,j (E) is the energy of the radionuclide decay once, Ek The fraction of radiation absorbed by the organism; M k is the mass of the external environment medium k, in kg; M j is the mass of target organ / tissue j in kg.

[0031] Furthermore, the activity concentration of each part of the organism is weighted according to the mass to obtain the average activity concentration of the organism, and this concentration and the geometric dimensions of the organism are input into the BiotaDC program to obtain the dose coefficient of the radionuclide for internal and external irradiation of the organism.

[0032] Furthermore, the activity concentration of each part of the organism is weighted according to the mass to obtain the average activity concentration of the organism, including the following steps: S41. Calculating the activity concentration of each organ;

[0033] S42. Weight the activity concentration of each organ by mass;

[0034] S43. Calculate the average activity concentration of organisms:

[0035] Where, activity concentration i is the activity concentration of organ / tissue i, in Bq / kg; mass i is the mass of organ / tissue i, in kg; ∑i mass is the total mass of all organs / tissues;

[0036] S44. Input the average activity concentration and geometric dimensions of the organism into the BiotaDC program and run it to obtain the dose coefficient of the radionuclide on the internal and external irradiation of the organism.

[0037] Furthermore, in MCNP, LAT cards, FILL cards, U cards, cell cards, surface cards, and data cards are used to write Monte Carlo voxel model input files.

[0038] Furthermore, in step S1, a tomographic scan is performed on the marine crab to obtain a tomographic sequence image of the organism, specifically including using a CT scanning device, setting the X-ray tube voltage to 120.0 kV, the tube current to 210.0 mA, the image scanning interval to 0.625 mm, the resolution to 28.346 pixels / cm, and the voxel size to 0.3527 mm×0.3527 mm×0.625 mm.

[0039] Furthermore, in step S2, the following steps are also included:

[0040] S21. To ensure that the subsequent processing steps can more accurately distinguish different tissues and organs, the acquired marine crab tomography images are first preprocessed in the image processing software, wherein the preprocessing includes denoising, contrast enhancement and image quality improvement;

[0041] S22. Use image segmentation technology to identify and extract different organs / tissues of marine crabs, and segment the crab shell, crab meat, heart, gills and crab roe;

[0042] S23. After the image segmentation is completed, a unique RGB value is assigned to the crab shell, crab meat, heart, gills and crab roe; the RGB value represents the color identification of the area;

[0043] S24. The segmented and labeled image data are saved as image files or structured data, which are used as input data for subsequent voxel matrix generation;

[0044] S25. Generate a three-dimensional matrix containing each tissue organ by reading and classifying the RGB value of each pixel in the image.

[0045] The present invention also discloses a device for estimating radiation dose of marine crabs based on a voxel model, comprising:

[0046] A processing unit, used to execute the above-mentioned method for estimating radiation dose of marine crabs based on voxel model;

[0047] A memory, used to store programs or instructions, wherein the programs or instructions are loaded and executed by the processing unit to implement the above steps;

[0048] An input module, for receiving tomographic image data of marine crabs and generating a simplified anatomical model;

[0049] A simulation module, for performing Monte Carlo simulation based on the simplified anatomical model and calculating radiation dose;

[0050] Output module for displaying radiation dose estimation results for marine crabs.

[0051] The present invention also discloses a computer-readable storage medium, on which a program or instruction is stored. When the program or instruction is executed by a processor, the method for estimating radiation dose of marine crabs based on a voxel model as described above is implemented.

[0052] Beneficial effects:

[0053] 1. Compared with the existing actual environmental assessment, conservative estimates are more accepted by people, but in scientific research, it is also very important to obtain accurate results. The model established this time is conducive to avoiding excessive protection in nuclear industry environmental assessment due to conservative estimates. Compared with the BiotaDC program, the crab voxel model provided by this patent has the advantage of distinguishing tissues and organs, and it is convenient to identify and segment tissues / organs according to CT tomographic sequence images, and the data is easy to obtain, which has good promotion and application value.

[0054] 2. The present invention uses CT scanning technology to obtain detailed anatomical data of marine crabs and construct an accurate voxel model, which can effectively distinguish organs such as crab shells, crab meat, heart, gills and crab roe. This process avoids the defect of traditional simplified models (such as ellipsoid models) that ignore organ differences, thereby improving the accuracy of radiation dose estimation. This method can accurately simulate the radiation absorption characteristics of each organ and provide more realistic radiation dose calculation results.

[0055] 3. Since the present invention is based on CT tomography and image processing technology, it can accurately obtain the three-dimensional anatomical data of the organism, so its technical solution has strong adaptability. It is not only suitable for marine crabs, but also can be flexibly adjusted according to the anatomical characteristics of different species, and is widely used in radiation dose estimation of other aquatic organisms or terrestrial species. This technical solution can be applied to multiple fields, such as ecological environment radiation monitoring, nuclear power plant radiation protection, etc., and has good promotion and application value. BRIEF DESCRIPTION OF THE DRAWINGS

[0056] Figure 1 This is a cross-sectional view of the crab voxel model in Example 1 of the present invention;

[0057] Figure 2 This is a comparison of the crab operation result data in Example 1 of the present invention. DETAILED DESCRIPTION

[0058] The present invention is further described in detail below in conjunction with the accompanying drawings and specific implementation methods.

[0059] Example 1. A method for estimating radiation dose of marine crabs based on a voxel model.

[0060] Step 1: Acquisition of biological tomographic sequence images First, perform CT scanning on marine crabs to obtain their tomographic sequence images. The parameters of CT scanning are as follows:

[0061] Tube voltage: 120.0 kV; tube current: 210.0 mA; image scanning interval: 0.625 mm; resolution: 28.346 pixels / cm; voxel size: 0.3527 mm × 0.3527 mm × 0.625 mm;

[0062] The scanned tomographic images contain detailed structural information of each organ and tissue in the crab's body. In order to ensure the accuracy of subsequent processing steps, each slice image obtained by scanning has been scanned at high resolution so that the tissues and organs inside the organism can be extracted more accurately.

[0063] Step 2: Tomographic sequence image processing

[0064] Image preprocessing: Preprocess the acquired CT images, including:

[0065] De-noising: Remove the noise caused by scanning equipment or other factors to make the image clearer.

[0066] Contrast enhancement: Enhance image contrast to make tissue boundaries more obvious, which helps with subsequent segmentation operations.

[0067] Image segmentation and labeling: Use image processing software to segment the image and extract different organs / tissues of marine crabs. According to the characteristics of each organ in the image (such as density, shape, texture, etc.), threshold segmentation or region growing method is used for segmentation. The specific segmented tissues and organs include: crab shell, crab meat, heart, gills, and crab roe.

[0068] Each segmented tissue and organ is assigned a unique RGB value, the crab shell is red (255, 0, 0), the crab meat is green (0, 255, 0), the heart is blue (0, 0, 255), the gills are yellow (255, 255, 0), and the crab roe is cyan (0, 255, 255) for subsequent pixel classification and identification.

[0069] Step 3: Voxel matrix acquisition

[0070] Image processing and pixel classification: Each pixel in each CT image is read one by one. Each pixel contains an RGB value, and based on these RGB values, it is determined which organ or tissue the pixel belongs to. The specific steps are as follows:

[0071] Crab shell: currentClass=1; Crab meat: currentClass=2; Heart: currentClass=3; Gills: currentClass=4; Crab roe: currentClass=5;

[0072] For each pixel in the image, first get its RGB value. If the RGB value matches the current class currentClass, then the counter currentNum is incremented by 1, indicating that the number of pixels in this class has increased. If the RGB value does not match the current class, then save the current class and its number of pixels, switch to the next class, and reset currentNum to 1 to start counting the number of pixels in the new class.

[0073] Generate voxel matrix: record the number of pixels of each category currentNum and its category identifier currentClass in the result list res. Finally, these data are converted into voxel matrix, each voxel represents a tissue or organ in a specific area of ​​the marine crab body, and records its spatial distribution, size, shape and other information. This information provides important basic data for subsequent Monte Carlo simulation.

[0074] Step 4: Generate MCNP input file by MCNP simulation to calculate biomass deposition energy: Write the input file required for MCNP simulation based on the generated voxel matrix data. The input file contains the spatial information, morphological information and assigned radiation source information of each organ / tissue in the marine crab, such as radionuclides.

[0075] The input files include: LAT card: material geometry definition, FILL card: definition of the distribution of structural units, U card: definition of radiation source, grid element card: definition of space division in the organism, surface card: definition of model surface, data card: setting calculation parameters.

[0076] MCNP simulation run:

[0077] Input the prepared MCNP input file into the MCNP software to perform Monte Carlo simulation. The deposition energy of radionuclides in the marine crab body is calculated through simulation, and the energy deposition value of each organ / tissue is obtained.

[0078] Activity concentration weighting and calculation:

[0079] The organ or tissue that concentrates a certain radionuclide is set as the source organ / tissue; then, the organ or tissue that is irradiated is set as the target organ / tissue; finally, the energy in the source organ / tissue k is E k The calculation formula of the internal radiation dose coefficient DC of the target organ / tissue j is:

[0080] In the formula, The energy in the source organ / tissue k is E k The radiation corresponds to the internal dose coefficient of the target organ / tissue j, in units of (μGy·h -1 ) / (Bq·kg -1 );Y k The energy of each decay of the radionuclide is E k The radiation yield; M k is the mass of the source organ / tissue in kg; F k,j (E) is the energy of the radionuclide decay once, E k The fraction of radiation energy absorbed by the organism;

[0081] M j is the mass of the target organ / tissue in kg;

[0082] The external radiation dose coefficient is calculated by treating all environmental media in which the organism is located as radiation sources. The calculation formula is:

[0083] In the formula, The mass is M kThe energy in the external environment medium is E k The radiation corresponds to the external radiation dose coefficient of the target organ / tissue j, in units of (μGy·h -1 ) / (Bq·kg -1 );Y k The energy of each decay of the radionuclide is E k The radiation yield; F k,j (E) is the energy of the radionuclide decay once, E k The fraction of radiation absorbed by the organism; M k is the mass of the external environment medium k, in kg; M j is the mass of target organ / tissue j in kg.

[0084] The activity concentration of each part of the organism is weighted according to mass to obtain the average activity concentration of the organism, including the following steps:

[0085] S41. Calculate the activity concentration of each organ;

[0086] S42. Weight the activity concentration of each organ by mass;

[0087] S43. Calculate the average activity concentration of organisms:

[0088] Where, activity concentration i is the activity concentration of organ / tissue i, in Bq / kg; mass i is the mass of organ / tissue i, in kg; ∑i mass is the total mass of all organs / tissues;

[0089] S44. Input the average activity concentration and geometric dimensions of the organism into the BiotaDC program and run it to obtain the dose coefficient of the radionuclide on the internal and external irradiation of the organism.

[0090] The average activity concentration and geometric dimensions of the organism are input into the BiotaDC program to calculate the internal and external radiation dose coefficients for each organ / tissue. These calculation results provide data support for subsequent radiation dose assessment.

[0091] Combination Figure 1 and 2 As shown, this embodiment constructs a voxel model of marine crabs and uses it for radiation dose estimation, and achieves the following key results:

[0092] Consistency between model and entity: The accuracy of the model was verified by comparing it with the mass and morphology of actual marine crabs. Comparison of radiation dose estimation with the BiotaDC program: The calculation results of the established voxel model were comparable to those of the BiotaDC program, indicating that the model can effectively reflect the actual situation.

[0093] Detailed analysis

[0094] 1. Comparison of the morphology and quality of voxel models with actual marine crabs

[0095] Combination Figure 1 and 2 In this embodiment, the mass of the voxel model of the marine crab is 106.49g, while the mass of the actual marine crab is 104.6g. The difference between the two is only 1.89g, indicating that the mass of the organism established by the voxel model is very close to the actual marine crab. This result shows that the model can accurately reproduce the size and volume of the marine crab, which is crucial for radiation dose estimation.

[0096] In addition, by comparing with the 3D images of the actual crab, the voxel model's shape is basically consistent with the actual shape, indicating that the model can accurately capture the crab's external structure. This result further verifies the effectiveness of the voxel model, especially in the application of radiation dose estimation.

[0097] 2. Mass and volume of organs / tissues

[0098] The volume, density and mass of each organ / tissue were extracted and calculated from the tomographic sequence images of marine crabs obtained from CT scans. The following are the relevant data of different organs / tissues:

[0099]

[0100] According to the table data, it can be seen that the crab shell and crab meat are the most significant parts of the marine crab body in terms of volume and mass. By accurately calculating the mass and volume of each organ, the model of the present invention can not only effectively estimate the dose in the organism, but also accurately provide the contribution of each organ.

[0101] 3. Comparison of model calculation results with BiotaDC program

[0102] By comparing the calculation results of the voxel model of the present invention with those of the BiotaDC program, we draw the following conclusions:

[0103] like Figure 2 As shown in the figure, the calculation results of the voxel model are comparable to those of the BiotaDC program, with little difference between the two. This shows that the results of this model have high accuracy and can reflect the actual radiation dose of marine crabs.

[0104] As the accuracy of the model increases, the calculation results become closer to the actual situation. In particular, when the details of each organ / tissue (such as tissue density, shape, etc.) are accurately divided in the model, the accuracy of the results is significantly improved.

[0105] It is worth noting that the BiotaDC program uses a relatively simplified model and fails to accurately distinguish between different organs and tissues of marine crabs. Therefore, its calculation results are relatively conservative and often overestimate the radiation dose. The present invention can provide a more accurate dose estimate by distinguishing between various organs and tissues.

[0106] The marine crab voxel model provided by this patent has the following advantages compared to the BiotaDC program:

[0107] Differentiating tissue organs: Through precise image segmentation, the present invention can distinguish different organs of marine crabs (such as crab shell, crab meat, heart, gills and crab roe), which is crucial for radiation dose assessment.

[0108] Based on CT tomographic sequence images: Using CT scanning technology and combining it with modern image processing algorithms, the accuracy of the model has been greatly improved. This method can not only effectively improve the calculation accuracy, but also make the data acquisition process simpler and more efficient.

[0109] Easy to promote and apply: Since CT scanning technology has been widely used in medical imaging and biological research, the technical solution used in this model has good promotion value and can be applied to the radiation dose estimation of other marine organisms or ecosystems.

[0110] Through the voxel model establishment and calculation method of the present invention, a high-precision radiation dose estimation of marine crabs was successfully achieved. The model is not only highly consistent with the actual organism in terms of mass and morphology, but can also effectively reflect the differences in radiation absorption of different organs, providing reliable data support for radiation protection and environmental assessment. Compared with the traditional simplified model, the voxel model of the present invention can provide more accurate results and has good scientific research value and application prospects.

[0111] Example 2: Marine crab radiation dose estimation device based on voxel model

[0112] This embodiment provides a device for estimating radiation dose of marine crabs based on a voxel model, which can generate a high-precision voxel model by inputting CT tomographic image data of marine crabs, and obtain radiation dose estimation results of marine crabs through Monte Carlo simulation calculation. The device includes a processing unit, a memory, an input module, a simulation module, and an output module.

[0113] Processing unit: The processing unit is used in the above-mentioned method for estimating radiation dose of marine crabs based on voxel model. Specifically, the processing unit performs image processing, organ segmentation, voxel matrix generation, and Monte Carlo simulation calculation according to the acquired CT image data, so as to obtain the radiation dose estimation result.

[0114] Memory: Memory is used to store programs or instructions that are loaded and executed by the processing unit to implement steps including image processing, organ segmentation, voxel matrix acquisition, and Monte Carlo simulation calculations. Memory can also save temporary data and calculation results.

[0115] Input module: The input module is used to receive the tomographic image data of marine crabs and generate a simplified anatomical model. The input module can also receive other environmental information related to marine crabs, such as the radioactivity concentration of water bodies, and pass this information to the processing unit for subsequent calculations.

[0116] Simulation module: The simulation module is used to perform Monte Carlo simulation based on the generated simplified anatomical model and calculate the radiation dose of each organ / tissue in the marine crab. The simulation module uses MCNP or similar Monte Carlo simulation software for dose estimation.

[0117] Output Module: The output module is used to display or store the radiation dose estimation results of marine crabs. The output module can display the results through a graphical user interface (GUI) or output them in the form of data files for further analysis.

[0118] The specific operation process is as follows:

[0119] Step 1. Data input: The user inputs the CT scan data of marine crabs into the input module, and the system automatically reads the CT image data.

[0120] Step 2. Image processing and organ segmentation: The processing unit preprocesses the image data, including denoising, contrast enhancement, etc. Then, the image segmentation algorithm is used to segment different organs in the image (such as crab shell, crab meat, heart, gills, and crab roe), and the RGB value of each organ is marked.

[0121] Step 3. Voxel matrix generation: The voxel matrix generation module converts the RGB value of each pixel into the corresponding organ / tissue category based on the processed image data. Then a three-dimensional voxel matrix is ​​generated, in which each voxel corresponds to a tissue organ and records its position, volume, density and other information.

[0122] Step 4. Monte Carlo simulation calculation: The simulation module uses Monte Carlo simulation software such as MCNP to calculate according to the voxel matrix and the input radionuclide information to obtain the radiation deposition energy of different organs / tissues of marine crabs.

[0123] Step 5. Calculate the dose: The processing unit calculates the radiation dose in the marine crab based on the deposition energy and the internal and external irradiation dose calculation formula, and the result is displayed in the output module.

[0124] Step 6. Result output and display: The output module displays the calculated radiation dose results in digital or graphical form for user analysis and decision-making.

[0125] Embodiment 3: Computer readable storage medium

[0126] This embodiment provides a computer-readable storage medium, in which a program or instruction is stored, and when the program or instruction is executed by a computer processor, a method for estimating the radiation dose of marine crabs based on a voxel model is implemented. The computer-readable storage medium can be in the form of a hard disk, a solid-state hard disk, an optical disk, a USB flash memory, etc., for storing program codes for executing the method of the present invention.

[0127] Includes the following parts:

[0128] Computer readable storage medium: The storage medium stores a program or instruction, which includes steps such as image processing, organ segmentation, voxel matrix generation, Monte Carlo simulation calculation, etc. The function of the program is to execute these steps in the above-mentioned embodiment 1.

[0129] Executing the program: When the program is loaded and executed by the processor, it controls the computer's hardware resources to perform tasks related to radiation dose estimation. The execution process includes receiving CT image data, performing image segmentation and processing, generating voxel matrices, running Monte Carlo simulations, and calculating radiation doses.

[0130] Usage scenario: The computer-readable storage medium is suitable for different types of computer systems, including personal computers, servers, embedded devices, etc. Whether it is offline analysis or online evaluation, the program code in the storage medium can be run in a variety of hardware environments to perform the task of estimating the radiation dose of marine crabs.

[0131] Storage and transmission: Users can use the storage medium for transmission and storage between different devices, facilitating data exchange, backup and sharing. The storage medium can be used not only for local computer execution, but also for distributed computing on cloud computing platforms or computing clusters.

[0132] The above shows and describes the basic principles, main features and advantages of the present invention. It should be understood by those skilled in the art that the present invention is not limited to the above embodiments, and the above embodiments and descriptions are only for explaining the principles of the present invention. Without departing from the spirit and scope of the present invention, the present invention may have various changes and improvements, which fall within the scope of the present invention to be protected. The scope of protection of the present invention is defined by the attached claims and their equivalents.

Claims

1. A method for estimating radiation dose of marine crabs based on voxel model, characterized in that: The following steps are involved: Step S1. Acquisition of tomographic sequence images of organisms: performing tomographic scanning on marine crabs to obtain tomographic sequence images of organisms; Step S2. Processing of tomographic sequence images: using image processing software to process the tomographic images, segmenting the organs / tissues of the organism, and assigning different RGB values ​​to each tissue and organ for labeling; Step S3. Voxel matrix acquisition: converting the processed image into a Monte Carlo repeating structure card input matrix according to the RGB value of each color to obtain a voxel matrix for Monte Carlo simulation, specifically including the following steps: S31. Image size definition: Define the length and width of each CT image to ensure consistent image size; S32. Pixel reading and classification: Set the category of the current pixel to identify which tissue or organ the pixel belongs to; initialize the current pixel category counter currentNum to 0 to record the number of pixels of each category; S33. Read each pixel in turn and obtain the RGB value: read each pixel of the image and obtain the RGB value of the pixel; the RGB value determines which organ or tissue category the pixel belongs to. According to the RGB value, determine which tissue / organ category the pixel belongs to, and set currentClass to the number of the category; S33. Determine whether the next pixel is a pixel of the current category and update the counter: Determine whether the next pixel belongs to the current category currentClass; If it is a pixel of the same category, the counter currentNum is increased by 1 and the counting continues; if the pixel does not belong to the current category currentClass, the current category and its number are stored in the result list res, and the currentClass is updated to the new pixel category, and then the counter currentNum is reset to 1; S34. Output the results and generate a voxel matrix: All classification results including the number of pixels in each category and their types are continuously output. The number of pixels in each category currentNum and the tissue category currentClass to which it belongs will be used as input data of the voxel matrix. The voxel matrix will be used for subsequent Monte Carlo simulation. Ultimately, the distribution, size and morphology of all tissues and organs are accurately recorded in the voxel matrix. Step S4. Calculate the deposition energy of organisms using MCNP simulation: input the voxel matrix into the MCNP software, and calculate the deposition energy of radionuclides on organisms through Monte Carlo simulation; according to the deposition energy, use the internal and external radiation dose calculation formula to calculate the internal and external radiation dose coefficients of marine crabs, and further calculate the dose rate.

2. The method for estimating radiation dose of marine crabs based on voxel model according to claim 1, characterized in that: The tissues or organs include crab shell, crab meat, heart, gills and crab roe.

3. The method for estimating radiation dose of marine crabs based on voxel model according to claim 1, characterized in that: The organ or tissue that concentrates a certain radionuclide is set as the source organ / tissue; then, the organ or tissue that is irradiated is set as the target organ / tissue; finally, the energy in the source organ / tissue k is E k The calculation formula of the internal radiation dose coefficient DC of the target organ / tissue j is: In the formula, The energy in the source organ / tissue k is E k The radiation corresponds to the internal dose coefficient of the target organ / tissue j, in units of (μGy·h -1 ) / (Bq·kg -1 );Y k The energy of each decay of the radionuclide is E k The radiation yield; M k is the mass of the source organ / tissue in kg; F k,j (E) is the energy of the radionuclide decay once, E k The fraction of radiation energy absorbed by the organism; M j is the mass of the target organ / tissue in kg; The external radiation dose coefficient is calculated by treating all environmental media in which the organism is located as radiation sources. The calculation formula is: In the formula, The mass is M k The energy in the external environment medium is E k The radiation corresponds to the external radiation dose coefficient of the target organ / tissue j, in units of (μGy·h -1 ) / (Bq·kg -1 ); Y k The energy of each decay of the radionuclide is E k The radiation yield; F k,j (E) is the energy of the radionuclide decay once, E k The fraction of radiation absorbed by the organism; M k is the mass of the external environment medium k, in kg; M j is the mass of target organ / tissue j in kg.

4. The method for estimating radiation dose of marine crabs based on voxel model according to claim 1, characterized in that: The activity concentration of each part of the organism is weighted according to mass to obtain the average activity concentration of the organism. This concentration and the geometric dimensions of the organism are input into the BiotaDC program to obtain the dose coefficient of radionuclides for internal and external irradiation of the organism.

5. The method for estimating radiation dose of marine crabs based on voxel model according to claim 4, characterized in that: The activity concentration of each part of the organism is weighted according to the mass to obtain the average activity concentration of the organism, including the following steps: S41. Calculating the activity concentration of each organ; S42. Weight the activity concentration of each organ by mass; S43. Calculate the average activity concentration of organisms: Where, activity concentration i is the activity concentration of organ / tissue i, in Bq / kg; Mass i is the mass of organ / tissue i in kg; ∑i mass is the total mass of all organs / tissues; S44. Input the average activity concentration and geometric dimensions of the organism into the BiotaDC program and run it to obtain the dose coefficient of the radionuclide on the internal and external irradiation of the organism.

6. The method for estimating radiation dose of marine crabs based on voxel model according to claim 1, characterized in that: In MCNP, LAT cards, FILL cards, U cards, cell cards, surface cards, and data cards are used to write Monte Carlo voxel model input files.

7. The method for estimating radiation dose of marine crabs based on voxel model according to claim 1, characterized in that: In the step S1, a tomographic scan is performed on the marine crab to obtain a tomographic sequence image of the organism, specifically including using a CT scanning device, setting the X-ray tube voltage to 120.0 kV, the tube current to 210.0 mA, the image scanning interval to 0.625 mm, the resolution to 28.346 pixels / cm, and the voxel size to 0.3527 mm×0.3527 mm×0.625 mm.

8. The method for estimating radiation dose of marine crabs based on voxel model according to claim 2, characterized in that: In the step S2, the following steps are also included: S21. To ensure that the subsequent processing steps can more accurately distinguish different tissues and organs, the acquired marine crab tomography images are first preprocessed in the image processing software, wherein the preprocessing includes denoising, contrast enhancement and image quality improvement; S22. Use image segmentation technology to identify and extract different organs / tissues of marine crabs, and segment the crab shell, crab meat, heart, gills and crab roe; S23. After the image segmentation is completed, a unique RGB value is assigned to the crab shell, crab meat, heart, gills and crab roe; the RGB value represents the color identification of the area; S24. The segmented and labeled image data are saved as image files or structured data, which are used as input data for subsequent voxel matrix generation; S25. Generate a three-dimensional matrix containing each tissue organ by reading and classifying the RGB value of each pixel in the image.

9. A device for estimating radiation dose of marine crabs based on voxel model, characterized in that: include: A processing unit, configured to execute the method for estimating radiation dose of marine crabs based on voxel model as claimed in any one of claims 1 to 8; A memory, used to store programs or instructions, wherein the programs or instructions are loaded and executed by the processing unit to implement the above steps; An input module, for receiving tomographic image data of marine crabs and generating a simplified anatomical model; A simulation module, for performing Monte Carlo simulation based on the simplified anatomical model and calculating radiation dose; Output module for displaying radiation dose estimation results for marine crabs.

10. A computer-readable storage medium, characterized in that: The storage medium stores a program or instruction, and when the program or instruction is executed by the processor, the method for estimating radiation dose of marine crabs based on a voxel model as described in any one of claims 1 to 8 is implemented.

Images