A method and system for calculating spatial efficiency based on 3D model and voxel partitioning

Abstract

This invention discloses a method and system for calculating spatial interaction efficiency based on a 3D model and voxel partitioning, relating to the field of 3D modeling and data processing technology. The method includes the following steps: voxelizing the 3D model of the target object into a material mesh; adaptively partitioning the source space into non-uniform voxels based on solid angle gradients; obtaining an anisotropic response correction operator through a lookup table; performing path tracing on the emitted rays of each source voxel and calculating the total path transmission weight using the correction operator; finally, numerically integrating based on the solid angle, transmission weight, and volume of each voxel to obtain the total interaction efficiency value, and fitting it into a continuous scale curve. This invention balances computational accuracy and speed through adaptive partitioning, and improves model realism through full path tracing and anisotropic correction, achieving efficient and accurate calculation of spatial interaction efficiency under complex conditions. It is applicable to scenarios such as passive efficiency calibration of radiation detectors.

Description

Technical Field

[0001] This invention relates to the field of 3D modeling and data processing technology, specifically to a method and system for calculating spatial efficiency based on 3D models and voxel partitioning. Background Technology

[0002] In the field of three-dimensional spatial efficiency calculation, such as radiation protection, optical rendering, and acoustic simulation, the core problem can be summarized as solving high-dimensional numerical integrals under complex geometric scenes. Traditional solution schemes are mainly divided into two categories: Monte Carlo methods and analytical integration methods. Monte Carlo methods approximate integration through random sampling. Although they can handle arbitrarily complex geometric and physical processes, they have inherent drawbacks such as slow convergence speed and extremely high computational cost, making it difficult to meet the needs of real-time or rapid analysis. The other category is analytical or semi-analytical methods based on deterministic mathematical models. Although they are computationally fast, they usually rely on oversimplification of geometric and physical processes, such as assuming homogeneity of the medium and ignoring anisotropy. When dealing with irregular three-dimensional shapes, uneven spatial distribution of materials, and sensitivity to the direction of the action target in real-world scenes, the modeling error is large and the accuracy is limited.

[0003] Therefore, there is an irreconcilable contradiction between the solution accuracy and computational efficiency in existing technologies, and there is an urgent need for a new and efficient numerical integration framework that can intelligently adapt to spatial complexity and finely integrate multiple physical effects. Summary of the Invention

[0004] The technical problem to be solved by the present invention is to address the shortcomings of the prior art by providing a method and system for calculating spatial interaction efficiency based on three-dimensional models and voxel division.

[0005] To achieve the above objectives, the technical solution adopted by the present invention is as follows: The spatial interaction efficiency calculation method based on 3D model and voxel division includes the following steps: Step S1: Obtain the three-dimensional surface model of the target and perform voxelization on the three-dimensional surface model to generate a voxel mesh with material properties. Step S2: For the source space, recursively divide it according to the preset adaptive subdivision criteria to generate a three-dimensional voxel set that is not uniformly distributed in space, wherein the adaptive subdivision criteria are based on the rate of change of the contribution of spatial position to efficiency. Step S3: Based on the spatial orientation of each voxel relative to the target and the current interaction energy, obtain the corresponding anisotropic response correction operator by querying the pre-built anisotropic response correction lookup table; Step S4: For each voxel in the three-dimensional voxel set, the voxel's full path transport weight is obtained by emitting a ray and traversing the voxel mesh generated in step S1, and by combining the anisotropic response correction operator. Step S5: Calculate the total spatial efficiency value based on the solid angle of each voxel to the target, the total path transport weight, and the volume of the voxel. Step S6: Based on the total spatial efficiency value, use an empirical formula to fit the data and generate a continuous efficiency-energy relationship curve.

[0006] Furthermore, in step S1, the voxelization process specifically includes the following steps: A uniform three-dimensional mesh is established within the bounding box space of the three-dimensional surface model; For each grid cell, calculate the signed distance function value between its center point and the three-dimensional surface model; Based on the sign of the signed distance function value, it is determined whether the mesh cell is located inside the model entity, and the mesh cell determined to be located inside the entity is assigned a corresponding material label to form a voxel mesh with material properties.

[0007] Furthermore, in step S2, the recursive division based on the preset adaptive subdivision criterion specifically includes the following steps: The axial bounding box of the aforementioned source space is used as the initial voxel unit; Calculate the geometric contribution gradient of the current voxel unit, which is evaluated by calculating the solid angle change of the target at the boundary vertex of the voxel unit; Based on a hybrid criterion of geometric contribution gradient, spatial distance, and voxel volume, a decision is made on whether to divide the current voxel node into eight equal sub-units. For each newly generated sub-unit, the steps of calculating the geometric contribution gradient and determining whether it is equally divided are recursively repeated until all units no longer meet the equal division condition, thus forming the non-uniformly distributed three-dimensional voxel set.

[0008] Furthermore, in step S3, the pre-constructed anisotropic response correction lookup table is generated through the following steps: At multiple discrete energy points, the actual response of the target to the action from different angles is measured or simulated. The actual response value is compared with the theoretical response value under the isotropic assumption to obtain a series of discrete anisotropic response correction operator data. Using angle and energy as independent variables, the anisotropic response correction operator data is interpolated to generate a continuous lookup table for querying.

[0009] Furthermore, step S4 specifically includes the following steps: Step S4.1: Fire a virtual ray from the current voxel center toward the geometric center of the target. Step S4.2: Using the three-dimensional digital differential analyzer algorithm, calculate the travel path of the virtual ray in the voxel mesh, and record the path length of the ray passing through voxels of different materials; Step S4.3: Calculate the total attenuation of the ray based on the energy-related attenuation coefficient of each material and the path length to obtain the full path transmission weight.

[0010] Furthermore, step S5 specifically includes the following steps: Step S5.1: Calculate the corresponding solid angle value for each voxel; Step S5.2: For each voxel, multiply the solid angle, the total path transport weight, and the volume of the voxel itself to obtain the individual contribution value of the voxel to the total efficiency. Step S5.3: Sum the individual contribution values ​​of all leaf node voxels to obtain the total spatial efficiency value of the action source space on the action target under the current action energy.

[0011] Furthermore, step S6 specifically includes the following steps: Step S6.1: For each energy point within the target energy range, calculate the corresponding total spatial efficiency value to obtain a set of discrete energy-efficiency data pairs. Step S6.2: Perform function fitting on the discrete energy-efficiency data pairs to generate a continuous efficiency-energy relationship curve; Step S6.3: Systematically store and manage the efficiency-energy relationship curves calculated for different sample shapes, volumes and different placement positions relative to the target, forming a standardized efficiency scale database.

[0012] A spatial interaction efficiency calculation system based on 3D models and voxel partitioning is provided to implement any of the spatial interaction efficiency calculation methods based on 3D models and voxel partitioning described in the preceding paragraph, including: The model processing module is used to convert the 3D surface model of the target into a voxel mesh with material properties. An adaptive meshing module is used to recursively subdivide the source space into a non-uniform set of voxels, with the subdivision criterion based on the geometric contribution gradient of the voxels to the target. The physical transmission weight calculation module is used to calculate the full-path transmission weight of each voxel in the voxel set. The anisotropic response module is used to look up a pre-set lookup table based on the voxel's angle and the current energy of the interaction to obtain the corresponding anisotropic response correction operator. The numerical integration module is used to calculate and sum the total spatial efficiency value based on the solid angle, total path transport weight, energy intrinsic efficiency and volume of each voxel; The central processing unit is communicatively connected to the model processing module, adaptive mesh generation module, physical transfer weight calculation module, anisotropic response module, and numerical integration module, and is used to coordinate and control the sequential execution and data transmission of each module.

[0013] Furthermore, the physical transmission weight calculation module specifically includes: A ray generation unit is used to define a ray path from the voxel center to the target center; Mesh traversal cells are used to calculate the traversal sequence of rays in a uniform voxel model and record the path lengths through voxels of different materials. The weighting synthesis unit is used to synthesize the full-path transmission weight based on the path length, the energy-related attenuation coefficient of the corresponding material, and the anisotropic response correction operator.

[0014] Furthermore, the system also includes a data storage and management module, which is connected to the model processing module, the anisotropic response module, the physical transmission weight calculation module, and the numerical integration module; the data storage and management module is used to centrally store and provide the following data: anisotropic response correction operator lookup table, attenuation coefficient lookup table for different materials and energies, energy intrinsic efficiency data, and generated efficiency-energy relationship curve data.

[0015] Compared with the prior art, the beneficial effects of the present invention are as follows: 1. This invention introduces adaptive voxel partitioning based on contribution gradients, which automatically refines the mesh in critical regions and sparses the mesh in secondary regions, achieving computational accuracy close to that of the Monte Carlo method with far less computational cost than uniform mesh.

[0016] 2. This invention employs full-path transmission weight calculation, accurately accumulates attenuation in different materials through ray tracing, and introduces anisotropic response correction operators to effectively correct the directional response differences of the detector itself, significantly improving the calculation accuracy under complex geometry and non-uniform media.

[0017] 3. The continuous efficiency calibration curve and the standard database that can be constructed are finally output by this invention, which are easy to integrate into measurement and analysis software and can be directly called by users, realizing the rapid transformation from theoretical algorithm to practical application.

[0018] 4. The method framework of this invention can be extended to other fields, such as optical and acoustic simulation, by replacing physical parameters. At the same time, as a passive method, it avoids the risks and costs of operating a radioactive source and simplifies the calibration process. Attached Figure Description

[0019] Other features, objects, and advantages of the invention will become more apparent from the following detailed description of non-limiting embodiments with reference to the accompanying drawings: Figure 1 This is a schematic diagram of the spatial efficiency calculation method provided in an embodiment of the present invention; Figure 2 This is a schematic diagram of adaptive voxel partitioning provided in an embodiment of the present invention; Figure 3 This is a schematic diagram of the spatial efficiency calculation system provided in an embodiment of the present invention. Detailed Implementation

[0020] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be described in detail below with reference to the accompanying drawings and specific embodiments.

[0021] like Figure 1 As shown, the method for calculating spatial interaction efficiency based on a 3D model and voxel partitioning includes the following steps: Step S1: Obtain the three-dimensional surface model of the target and perform voxelization on the three-dimensional surface model to generate a voxel mesh with material properties. Step S2: For the source space, recursively divide it according to the preset adaptive subdivision criteria to generate a three-dimensional voxel set that is not uniformly distributed in space, wherein the adaptive subdivision criteria are based on the rate of change of the contribution of spatial position to efficiency. Step S3: Based on the spatial orientation of each voxel relative to the target and the current interaction energy, obtain the corresponding anisotropic response correction operator by querying the pre-built anisotropic response correction lookup table; Step S4: For each voxel in the three-dimensional voxel set, the voxel's full path transport weight is obtained by emitting a ray and traversing the voxel mesh generated in step S1, and by combining the anisotropic response correction operator. Step S5: Calculate the total spatial efficiency value based on the solid angle of each voxel to the target, the total path transport weight, and the volume of the voxel. Step S6: Based on the total spatial efficiency value, use an empirical formula to fit the data and generate a continuous efficiency-energy relationship curve.

[0022] In step S1, the voxelization process specifically includes the following steps: A uniform three-dimensional mesh is established within the bounding box space of the three-dimensional surface model; For each grid cell, calculate the signed distance function value between its center point and the three-dimensional surface model; Based on the sign of the signed distance function value, it is determined whether the mesh cell is located inside the model entity, and the mesh cell determined to be located inside the entity is assigned a corresponding material label to form a voxel mesh with material properties.

[0023] Obtain an accurate three-dimensional surface mesh model of the target, which should distinguish different physical components, such as the core sensor, shielding layer, and outer shell; An initial uniform 3D mesh is established within the bounding box of the model. For each mesh cell, i.e., voxel, the signed distance function SDF value from its center point to the target surface is calculated. If the SDF value is ≤ 0, the voxel is located inside the target entity and is marked as valid; If the SDF value is >0, the voxel is located outside the target or outside the surface and is marked as invalid.

[0024] For each valid voxel, a corresponding material label is assigned based on the physical component to which its spatial location belongs. Each material type is associated with a physical property table, including a linear attenuation coefficient, which represents the material's ability to attenuate an action of energy E.

[0025] like Figure 2 As shown, in step S2, the recursive division based on the preset adaptive subdivision criterion specifically includes the following steps: The axial bounding box of the aforementioned source space is used as the initial voxel unit; Calculate the geometric contribution gradient of the current voxel unit, which is evaluated by calculating the solid angle change of the target at the boundary vertex of the voxel unit; Based on a hybrid criterion of geometric contribution gradient, spatial distance, and voxel volume, a decision is made on whether to divide the current voxel node into eight equal sub-units. For each newly generated sub-unit, the steps of calculating the geometric contribution gradient and determining whether it is equally divided are recursively repeated until all units no longer meet the equal division condition, thus forming the non-uniformly distributed three-dimensional voxel set.

[0026] Define the axial bounding box of the action source space as the initial root voxel; for the voxel to be evaluated, assess the rate of change of its spatial position's contribution to the efficiency of the action target, and define its geometric contribution gradient as the maximum value of the solid angle change at the voxel's boundary, specifically including: The solid angles of the eight vertices of the voxel with respect to the geometric center of the target are calculated. The solid angles are calculated by approximately projecting the target onto a unit sphere centered at the target's geometric center and then performing an area integral. The specific formula is as follows: in, Indicates the geometric center of the target region. This indicates that the target region is perpendicular to the vector. The projected area on the plane, Point Approximate solid angle of the target region; Calculate the ratio of the solid angle variation difference of the voxel along each axis to the side length, and take the maximum value as the gradient. The specific formula is as follows: in, , and These represent the differences in solid angles between the two vertices of a voxel along the x, y, and z axes, respectively. , and Let x, y, and z represent the side lengths of the voxel, respectively. This represents the geometric contribution gradient.

[0027] Set gradient threshold Near-field distance threshold and maximum voxel volume threshold For the current voxel, if any of the following conditions are met, it is divided into eight smaller sub-voxels, specifically including: 1. Gradient criterion: ,in, Indicates the average side length of a voxel; 2. Near-field criterion: ; 3. Volume criterion: ; For each newly generated sub-voxel, the geometric contribution gradient and adaptive subdivision criterion are recursively calculated until all voxels no longer meet the subdivision conditions, ultimately generating a three-dimensional voxel set with non-uniform depth. Voxels are dense in regions near the target or where contribution changes drastically, and sparse in regions away from the target.

[0028] In step S3, the pre-constructed anisotropic response correction lookup table is generated through the following steps: At multiple discrete energy points, the actual response of the target to the action from different angles is measured or simulated. The actual response value is compared with the theoretical response value under the isotropic assumption to obtain a series of discrete anisotropic response correction operator data. Using angle and energy as independent variables, the anisotropic response correction operator data is interpolated to generate a continuous lookup table for querying.

[0029] The anisotropic response correction operator data were pre-obtained through experimental calibration and data fitting. Specifically, in a laboratory environment, a standard point source emitting a single-energy action was used, moving around the target center within a fixed step size range of 0° to 180° while measuring the target's response. The actual response efficiency measured at different angles was compared with the theoretical response value of an ideal homogeneous target at the same energy to obtain a set of discrete anisotropic response correction operator data. Finally, a continuous lookup table for energy and angle was generated through function interpolation and embedded into the system for rapid lookup during calculation.

[0030] Step S4 specifically includes the following steps: Step S4.1: Fire a virtual ray from the current voxel center toward the geometric center of the target. Step S4.2: Using the three-dimensional digital differential analyzer algorithm, calculate the travel path of the virtual ray in the voxel mesh, and record the path length of the ray passing through voxels of different materials; Step S4.3: Based on the energy-related attenuation coefficients of each material and the path length, calculate the total attenuation of the ray to obtain the full path transmission weight. The specific formula is as follows: in, This represents the total path transport weight of voxel i. This represents the attenuation coefficient of the j-th medium on the path for the energy E. Let represent the path length along the line connecting the center of source voxel i to the geometric center of the target, passing through the j-th medium. This represents an anisotropic response correction operator, used to correct the target's response as a function of the incident direction due to its internal physical structure. The changing phenomena, among which, The angle between the direction of the action from the voxel incident to the center of the target and the direction of the target axis. Indicates the total number of media.

[0031] Step S5 specifically includes the following steps: Step S5.1: Calculate the corresponding solid angle value for each voxel; Step S5.2: For each voxel, multiply the solid angle, the total path transport weight, and the voxel's own volume to obtain the voxel's individual contribution to the overall efficiency. The specific formula is as follows: in, This represents the individual contribution value of voxel i. Let represent the solid angle of the i-th voxel relative to the detector's sensitive crystal. The intrinsic energy efficiency of a target, i.e., the probability that an action with energy E will be recorded under ideal conditions, is an inherent property of the detector. Represents the volume of voxel i; The intrinsic energy efficiency is obtained by calling the detector's factory calibration data, consulting a standard database, or pre-calibrating it through a standard point source experiment. Step S5.3: Sum the individual contribution values ​​of all leaf node voxels to obtain the total spatial efficiency value of the action source space on the action target at the current action energy. The specific formula is as follows: in, This represents the overall space efficiency value. Indicates the number of voxels.

[0032] Step S6 specifically includes the following steps: Step S6.1: For each energy point within the target energy range, calculate the corresponding total spatial efficiency value to obtain a set of discrete energy-efficiency data pairs. Step S6.2: Perform function fitting on the discrete energy-efficiency data pairs to generate a continuous efficiency-energy relationship curve. The specific formula for the function fitting is as follows: in, Represents the fitting coefficient. The order of the polynomial is represented by the least squares method, and the coefficients are determined to obtain the optimal fitting curve. Step S6.3: Systematically store and manage the efficiency-energy relationship curves calculated for different sample shapes, volumes and different placement positions relative to the target, forming a standardized efficiency scale database.

[0033] like Figure 3 As shown, a spatial interaction efficiency calculation system based on a 3D model and voxel partitioning is used to implement any of the spatial interaction efficiency calculation methods based on a 3D model and voxel partitioning, including: The model processing module is used to convert the 3D surface model of the target into a voxel mesh with material properties. An adaptive meshing module is used to recursively subdivide the source space into a non-uniform set of voxels, with the subdivision criterion based on the geometric contribution gradient of the voxels to the target. The physical transmission weight calculation module is used to calculate the full-path transmission weight of each voxel in the voxel set. The anisotropic response module is used to look up a pre-set lookup table based on the voxel's angle and the current energy of the interaction to obtain the corresponding anisotropic response correction operator. The numerical integration module is used to calculate and sum the total spatial efficiency value based on the solid angle, total path transport weight, energy intrinsic efficiency and volume of each voxel; The central processing unit is communicatively connected to the model processing module, adaptive mesh generation module, physical transfer weight calculation module, anisotropic response module, and numerical integration module, and is used to coordinate and control the sequential execution and data transmission of each module.

[0034] The physical transmission weight calculation module specifically includes: A ray generation unit is used to define a ray path from the voxel center to the target center; Mesh traversal cells are used to calculate the traversal sequence of rays in a uniform voxel model and record the path lengths through voxels of different materials. The weighting synthesis unit is used to synthesize the full-path transmission weight based on the path length, the energy-related attenuation coefficient of the corresponding material, and the anisotropic response correction operator.

[0035] The data storage and management module is connected to the model processing module, the anisotropic response module, the physical transmission weight calculation module, and the numerical integration module. The data storage and management module is used to centrally store and provide the following data: anisotropic response correction operator lookup table, attenuation coefficient lookup table for different materials and energies, energy intrinsic efficiency data, and generated efficiency-energy relationship curve data.

[0036] Any combination of one or more computer-readable media may be used. A computer-readable medium can be a computer-readable signal medium or a computer-readable storage medium. A computer-readable storage medium can be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples (a non-exhaustive list) of computer-readable storage media include: an electrical connection having one or more wires, a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage device, magnetic storage device, or any suitable combination thereof. In this document, a computer-readable storage medium can be any tangible medium that contains or stores a program that can be used by or in connection with an instruction execution system, apparatus, or device.

[0037] The examples described herein are merely preferred embodiments of the invention and are not intended to limit the concept and scope of the invention. Any modifications and improvements made by those skilled in the art to the technical solutions of the invention without departing from the design concept of the invention should fall within the protection scope of the invention.

Claims

1. A method for calculating spatial interaction efficiency based on 3D models and voxel partitioning, characterized in that, Includes the following steps: Step S1: Obtain the three-dimensional surface model of the target and perform voxelization on the three-dimensional surface model to generate a voxel mesh with material properties. Step S2: For the source space, recursively divide it according to the preset adaptive subdivision criteria to generate a three-dimensional voxel set that is not uniformly distributed in space, wherein the adaptive subdivision criteria are based on the rate of change of the contribution of spatial position to efficiency. Step S3: Based on the spatial orientation of each voxel relative to the target and the current interaction energy, obtain the corresponding anisotropic response correction operator by querying the pre-built anisotropic response correction lookup table; Step S4: For each voxel in the three-dimensional voxel set, the voxel's full path transport weight is obtained by emitting rays and traversing the voxel mesh generated in step S1, and by combining the anisotropic response correction operator. Step S5: Calculate the total spatial efficiency value based on the solid angle of each voxel to the target, the total path transport weight, and the volume of the voxel. Step S6: Based on the total spatial efficiency value, use an empirical formula to fit the data and generate a continuous efficiency-energy relationship curve.

2. The method according to claim 1, characterized in that, In step S1, the voxelization process specifically includes the following steps: A uniform three-dimensional mesh is established within the bounding box space of the three-dimensional surface model; For each grid cell, calculate the signed distance function value between its center point and the three-dimensional surface model; Based on the sign of the signed distance function value, it is determined whether the mesh cell is located inside the model entity, and the mesh cell determined to be located inside the entity is assigned a corresponding material label to form a voxel mesh with material properties.

3. The method according to claim 2, characterized in that, In step S2, the recursive division based on the preset adaptive subdivision criterion specifically includes the following steps: The axial bounding box of the aforementioned source space is used as the initial voxel unit; Calculate the geometric contribution gradient of the current voxel unit, which is evaluated by calculating the solid angle change of the target at the boundary vertex of the voxel unit; Based on a hybrid criterion of geometric contribution gradient, spatial distance, and voxel volume, a decision is made on whether to divide the current voxel node into eight equal sub-units. For each newly generated sub-unit, the steps of calculating the geometric contribution gradient and determining whether it is equally divided are recursively repeated until all units no longer meet the equal division condition, thus forming the non-uniformly distributed three-dimensional voxel set.

4. The method according to claim 3, characterized in that, In step S3, the pre-constructed anisotropic response correction lookup table is generated through the following steps: At multiple discrete energy points, the actual response of the target to the action from different angles is measured or simulated. The actual response value is compared with the theoretical response value under the isotropic assumption to obtain a series of discrete anisotropic response correction operator data. Using angle and energy as independent variables, the anisotropic response correction operator data is interpolated to generate a continuous lookup table for querying.

5. The method according to claim 4, characterized in that, Step S4 specifically includes the following steps: Step S4.1: Fire a virtual ray from the current voxel center toward the geometric center of the target. Step S4.2: Using the three-dimensional digital differential analyzer algorithm, calculate the travel path of the virtual ray in the voxel mesh, and record the path length of the ray passing through voxels of different materials; Step S4.3: Calculate the total attenuation of the ray based on the energy-related attenuation coefficient of each material and the path length to obtain the full path transmission weight.

6. The method according to claim 5, characterized in that, Step S5 specifically includes the following steps: Step S5.1: Calculate the corresponding solid angle value for each voxel; Step S5.2: For each voxel, multiply the solid angle, the total path transport weight, and the volume of the voxel itself to obtain the individual contribution value of the voxel to the total efficiency. Step S5.3: Sum the individual contribution values ​​of all leaf node voxels to obtain the total spatial efficiency value of the action source space on the action target under the current action energy.

7. The method according to claim 6, characterized in that, Step S6 specifically includes the following steps: Step S6.1: For each energy point within the target energy range, calculate the corresponding total spatial efficiency value to obtain a set of discrete energy-efficiency data pairs. Step S6.2: Perform function fitting on the discrete energy-efficiency data pairs to generate a continuous efficiency-energy relationship curve; Step S6.3: Systematically store and manage the efficiency-energy relationship curves calculated for different sample shapes, volumes and different placement positions relative to the target, forming a standardized efficiency scale database.

8. A spatial interaction efficiency calculation system based on a three-dimensional model and voxel partitioning, used to implement the spatial interaction efficiency calculation method based on a three-dimensional model and voxel partitioning as described in any one of claims 1-7, characterized in that, include: The model processing module is used to convert the 3D surface model of the target into a voxel mesh with material properties. An adaptive meshing module is used to recursively subdivide the source space into a non-uniform set of voxels, with the subdivision criterion based on the geometric contribution gradient of the voxels to the target. The physical transmission weight calculation module is used to calculate the full-path transmission weight of each voxel in the voxel set. The anisotropic response module is used to look up a pre-set lookup table based on the voxel's angle and the current energy of the interaction to obtain the corresponding anisotropic response correction operator. The numerical integration module is used to calculate and sum the total spatial efficiency value based on the solid angle, total path transport weight, energy intrinsic efficiency and volume of each voxel; The central processing unit is communicatively connected to the model processing module, adaptive mesh generation module, physical transfer weight calculation module, anisotropic response module, and numerical integration module, and is used to coordinate and control the sequential execution and data transmission of each module.

9. The system according to claim 8, characterized in that, The physical transmission weight calculation module specifically includes: A ray generation unit is used to define a ray path from the voxel center to the target center; Mesh traversal cells are used to calculate the traversal sequence of rays in a uniform voxel model and record the path lengths through voxels of different materials. The weighting synthesis unit is used to synthesize the full-path transmission weight based on the path length, the energy-related attenuation coefficient of the corresponding material, and the anisotropic response correction operator.

10. The system according to claim 9, characterized in that, The system also includes a data storage and management module, which is connected to the model processing module, the anisotropic response module, the physical transmission weight calculation module, and the numerical integration module. The data storage and management module is used to centrally store and provide the following data: anisotropic response correction operator lookup table, attenuation coefficient lookup table for different materials and energies, energy intrinsic efficiency data, and generated efficiency-energy relationship curve data.

Images