36 results about "Radial coordinate" patented technology

A method for modeling acoustics in a delimited area, i.e., a computer-implemented method of modeling acoustics of an object as a numerical simulation model to determine field variables within a delimited control volume, the method comprising: (a) providing a coordinate system mapping positions in the control volume (b) providing a geometric specification of a sound source in the control volume based on the coordinate system. The objective of the invention is to improve the efficiency of numerical modeling of a transient (vibration) acoustic wave propagation phenomenon in a boundary-free domain. Additional steps are proposed: (c) constructing a convex shape around a sound source using a fast convex hull algorithm (d) stretching an infinite element starting from the convex shape (e) constructing a field shape function (g) using auxiliary coordinates by mapping radial coordinates of a coordinate system to auxiliary coordinates extending from the convex shape to the infinite element to transform the coordinate system into a field variable coordinate system (f1) Determining field variables by solving a numerical simulation model provides a computer system for simulating a system by applying a computer-implemented method.

A computer-implemented method for modeling acoustics of an object as a numerical simulation model in an unbounded domain to determine a field variable within a control volume of the unbounded domain includes providing a coordinate system mapping positions in the control volume. A geometrical specification of a sound source is provided in the control volume based on the coordinate system. To improve efficiency of numerically modelling transient acoustics wave propagation phenomena in unbounded domains, a convex shape is constructed around the sound source using Quickhull algorithm, infinite elements are extruded starting from the convex shape, and the coordinate system is transformed into a field variable coordinate system by mapping a radial coordinate of the coordinate system to an auxiliary coordinate extending from the convex shape to the infinite elements. Field shape functions are constructed using the auxiliary coordinate, and the field variable is determined by solving the numerical simulation model.

A computer-implemented method for modeling acoustics of an object (OBJ) as a model (MDL) in an unbounded domain to determine a field variable (FDV) within a control volume (CVL) of the unbounded domain, comprising: (a) Providing a coordinate system (COS) mapping the positions (POS) in said control volume (CVL), (b) Providing a geometrical model of a sound source (SDS) in said control volume (CVL) based on said coordinate system. To improve the efficiency of numerically modelling transient (vibro-) acoustics wave propagation phenomena in unbounded domains the additional steps are proposed: (c) Constructing a convex shape (CXS) around said sound source (SDS) using Quickhull algorithm (QHA), (d) Extruding infinite elements (IFE) starting from the convex shape (CXS), (e) Transforming the coordinate system (COS) into a field variable coordinate system (FVC) by mapping a radial coordinate of the coordinate system (COS) to an auxiliary coordinate (v(s, t, v)) extending from the convex shape (CXS) to the infinite elements (IFE), (f) Constructing of field shape functions (FSF) using the auxiliary coordinate (v(s,t,v)), (g) Determining the field variable (FDV) by solving the model (MDL). Furthermore, a computer-system (CPS) for simulating a system by applying the computer-implemented method is provided.

The invention belongs to the field of water pump runner blade design, and particularly relates to a water pump runner with an S swept shape, a design method and an optimization method thereof. The method comprises the steps that 1, a rotating wheel three-dimensional parameterized model is established, axial coordinates z, radial coordinates r and circumferential wrap angles theta of control points P1-P5 serve as design variables, hydraulic efficiency eta and lift H serve as optimization targets, and iterative optimization is conducted through a multi-target genetic algorithm; 2, constraint conditions in the optimization process include that the curvature radius change rate between adjacent control points is smaller than or equal to 15%; the maximum deflection of the S-shaped curve of the outlet edge is less than or equal to 0.1 D0; the rounding radius ratio RP2 / RP4 belongs to [1.5, 2.0]; and 3, fitting the relationship between the design variable and the target function through the response surface model, and finally outputting a Pareto optimal solution set. The problem that an existing water pump is single in design parameter is solved, blade surface pressure gradient distribution is improved, the vortex intensity of a wake flow area is reduced, and the performance and efficiency of the water pump are improved.

The invention relates to the technical field of screw head size measurement, in particular to a high-precision size measurement method and system for a screw head, and the method comprises the steps: obtaining a hardness factor of the screw head in each measurement process and a radial coordinate set at each rotation angle; obtaining the groove profile matching degree of each measurement process according to the difference between the fitting function of the radial coordinate set under each rotation angle in each measurement process and the ideal groove profile function of the screw head groove profile; according to the fusion result of the groove profile matching degree and the hardness factor of all the current measurement processes and the difference between the predicted value of the groove profile matching degree of the next measurement process and the average level of the groove profile matching degree of all the current measurement processes, the current groove profile risk coefficient of the measured screw is obtained, and whether a closed-loop control mechanism is triggered or not is judged; and the size of the screw head groove is obtained. According to the method, the assembly failure risk of the screw is analyzed, the angular resolution during laser scanning is adjusted, and the precision of screw head size measurement is improved.

The invention relates to a dressing tool for dressing a grinding worm for generating a workpiece with pre-formed teeth, which dressing tool is designed for rotation about a dressing rotation axis (A) and has a dressing side surface which forms a circular ring in a projection plane perpendicular to the dressing rotation axis of the dressing tool, wherein the circular ring defines a radial direction having radial coordinates (r) and a circumferential direction having angular coordinates (phi). The trimming side surface is provided with a trimming side modification which is produced in a targeted manner, and the trimming side modification can be described as a two-dimensional Fourier series related to a radial coordinate (r) and an angular coordinate (phi) in a manner projected onto a circular ring. The Fourier series has at least one non-zero Fourier coefficient for a non-zero circumferential frequency component with respect to the angular coordinates.

The application relates to a generation method of a three-dimensional body grid based on radial growth of a conical surface grid, and belongs to the field of electromagnetic engineering.The generation method comprises the following steps: reading geometric parameters and discretization data of the conical surface grid, screening triangular surface units needing to be processed according to target material numbers, and generating a unit list to be processed; a local coordinate system C is established with a conical axis as a reference, vertex coordinates are converted to the local coordinate system C, and transformed coordinates are obtained; radial coordinate mapping is performed on the triangular surface units in the unit list to be processed, each mapping point is inversely transformed to a world coordinate system, and grown vertex coordinates are obtained; a three-prism body grid is generated by connecting original triangular surface units and the grown vertex coordinates, the three-prism body grid is decomposed into tetrahedron units, and a final body grid is output.The application realizes geometric self-adaptive deformation, ensures unit quality, controls the conical surface curvature error within 0.1%, and supports efficient generation of ten-million-level grids.

A turbocharger turbine wheel can include a hub including an axis of rotation, a back disk, and a nose, wherein the axis of rotation defines an axial coordinate (z) in a cylindrical coordinate system including a radial coordinate (r) and a azimuthal coordinate ( ) along an intended direction of rotation about the axis of rotation; and blades extending outwardly from the hub, wherein each of the blades includes a hub profile, a shroud edge, a leading edge, a trailing edge, a pressure side, and a suction side, wherein the hub profile includes a global maximum radius and a global minimum radius, and wherein, between the global maximum radius and the global minimum radius, in an axial direction from the back disk to the nose, the hub profile includes a local minimum radius at a first axial coordinate location and a local maximum radius at a second axial coordinate location.

The application discloses a fan blade grid dynamic matching method for two-way fluid-structure coupling analysis. The blade discrete deformation data output from a structure solver is used to extract the spanwise coordinate and displacement component, the neutral axis deflection sample is obtained by adopting the robust statistics of the box, the main deflection function is constructed by adaptively fitting Chebyshev polynomial, and the monotone regression constraint is applied. Then, the second-order deformation strength function is identified based on the curvature and the radial coordinate of the cross section, the axial projection shortening term is constructed by combining the large deflection geometric relationship, and the parametric analytical displacement field is formed. The above function is written into a parameter file in a small number of polynomial coefficients, the displacement is evaluated in real time according to the initial coordinates of the nodes in the dynamic grid of the fluid solver, and the coordinate-driven mapping without grid consistency is realized. The application effectively solves the displacement mapping difficulty caused by the inconsistent interface grid, reduces the calculation cost, avoids the dynamic grid distortion, and is suitable for the two-way fluid-structure coupling analysis of the fan blade.

The invention discloses a PDC drill bit tooth arrangement method, and relates to the technical field of drilling tools. According to the expected drill bit mechanical drilling speed and the drill bit rotating speed; the drilling footage of a drill bit in each rotation of a circle, the total volume of rock cut off in each rotation of the drill bit, the average cutting volume and the average cutting area of each cutting tooth in each rotation of the drill bit, the maximum cutting area of a single cutting tooth, and the minimum rotation radius of each cutting tooth for realizing equal cutting volume are calculated in sequence; the average cutting area is larger than the average cutting volume of each cutting tooth and the average cutting area of each cutting tooth with the minimum rotation radius for achieving the equal cutting volume, and then the cutting section of each cutting tooth is equivalent to a quadrangle, so that the relation between the radial coordinates of the left and right boundaries of the section of each cutting tooth and the average cutting volume of each cutting tooth is obtained; the radial coordinate of the section centroid of the cutting tooth is calculated; and finally, the positioning radius and height of the cutting teeth are calculated by combining the crown profile curve equation of the drill bit and the radius of the cutting teeth, so that the service life of the arranged drill bit is long.

The present application relates to the technical field of optical detection, and particularly relates to a kind of holographic detection optical element processing error global inversion method, comprising: constructing the function relationship of etching depth error about local grating period reciprocal and radial coordinate, as initial etching depth error physical parameterization model;The function relationship of duty cycle error about local grating period reciprocal and radial coordinate is constructed, as initial duty cycle error physical parameterization model;Obtain multiple groups of measured parameters;Multiple groups of measured parameters are substituted into initial etching depth error physical parameterization model and initial duty cycle error physical parameterization model, and etching depth error physical parameterization model and duty cycle error physical parameterization model are obtained;Full aperture design data is brought into etching depth error physical parameterization model and duty cycle error physical parameterization model to obtain etching depth error global distribution and duty cycle error global distribution.The present application is at least beneficial to realize low-cost, high-precision error reconstruction.

The invention discloses a method for determining the maximum radial bending moment of a variable thickness circular plate under the action of a non-uniformly distributed load: a clamping device with an inner radius of a is used to fix and clamp the periphery of a variable thickness circular plate with a Young's elastic modulus of E, a Poisson's ratio of v, a thickness at a center point of h0, a thickness variation coefficient of α, and a thickness h that changes from the center along a radial coordinate r to the surroundings in a natural logarithmic function #imgabs0#, thereby forming a variable thickness circular plate structure with a radius of a and a fixed clamped periphery, and a load of q0 at a center point, a load variation coefficient of β, and a load that changes from the center along a radial coordinate r to the surroundings in a quadratic function q0(1+βr 2 / a 2 ) changes in the non-uniformly distributed load q, causing the variable thickness circular plate to produce axisymmetric deformation. Based on the static equilibrium analysis of the axisymmetric deformation problem of the variable thickness circular plate, the maximum radial bending moment M after the axisymmetric deformation of the variable thickness circular plate can be determined by using the measured value of the load q0 at the center point. rm .

The invention relates to the technical field of myopic lenses, and discloses a first-order differential offset-free concentric ring alternating type outer ring defocusing method and a myopic lens. The method comprises the following steps: establishing a radial coordinate system by taking the center of a lens as an original point, and dividing the surface of the lens into a central visual area and an outer ring defocus area; the M concentric rings of the outer ring defocus area are divided into odd rings and even rings to form a ring belt structure with alternating wide and thin rings; and at the boundary position rm, establishing constraint conditions of the mth ring and the (m + 1) th ring in the ring belt structure, and solving a curved surface function of the (m + 1) th ring to form a first-order differential offset-free annular out-of-focus curved surface. The problem of optical aberration caused by first-order differential offset at the junction of the ring belt in the prior art is solved, the first-order differential offset-free design is adopted, and the optical performance and wearing comfort of the lens are improved.

The present application relates to the technical field of data processing. Disclosed are a data encryption method and apparatus, a data decryption method and apparatus, and a storage medium. The encryption method comprises: converting into a polar coordinate a coordinate of each piece of vector data in a vector data set, wherein each piece of vector data comprises at least one vertex, and the polar coordinate comprises a polar angular coordinate and a polar radial coordinate which correspond to the at least one vertex; with regard to any polar coordinate, executing the following operations: sequentially rotating the polar angular coordinate by at least two encryption angles according to a first direction, so as to obtain an encrypted polar angular coordinate, and adding preset watermark information to a valid bit of the polar radial coordinate, so as to obtain a watermark-containing polar radial coordinate, wherein the at least two encryption angles are angle values respectively corresponding to at least two keys in a random rotation sequence, and the angle values in the random rotation sequence are generated on the basis of the keys; and obtaining encrypted vector data on the basis of the encrypted polar angular coordinate and the watermark-containing polar radial coordinate. Therefore, convenient and efficient data encryption and watermark embedding are achieved.

The invention discloses a shield tunnel monitoring method, and the method comprises the steps: obtaining the three-dimensional point cloud data of a measurement point in a shield tunnel, carrying out the coordinate transformation of the three-dimensional point cloud data, obtaining the target point cloud data, and enabling the target point cloud data to comprise the target point cloud data of the measurement point, the target point cloud data of the measurement point comprises the longitudinal coordinate, the azimuth angle coordinate and the radial coordinate of the measurement point, the structure of the shield tunnel is monitored based on the target point cloud data of the measurement point, and the monitoring result of the shield tunnel is obtained, so that the structure of the shield tunnel can be represented by the point cloud data of the longitudinal coordinate, the azimuth angle coordinate and the radial coordinate; the structure of the shield tunnel can be analyzed and monitored in the longitudinal direction, the radial direction and the rotation angle, and complex and indirect three-dimensional rectangular coordinate change analysis is simplified into direct comparison of longitudinal coordinates, radial coordinates and angle coordinates; and the intuition, accuracy and efficiency of longitudinal, circumferential, inter-ring slab staggering and torsional deformation identification are improved.

The invention discloses a magnetic resonance system fault diagnosis method based on multi-parameter radial coordinate visualization, and belongs to the technical field of magnetic resonance. The method comprises the following steps: firstly, performing contribution rate sorting on monitoring parameters under different fault types according to priori knowledge, and rearranging variables; arranging the sorted multiple parameters on the circumference at equal intervals to form a radial coordinate initial projection layout, and iteratively optimizing the weight and the anchor point angle of each dimension by adopting an intelligent optimization algorithm to maximize the Euclidean distance between the fault data and the normal data in a two-dimensional projection space; and carrying out radial coordinate projection on the current test data based on the optimized projection parameters, and judging whether the system has a fault and the specific fault type by calculating Euclidean distances between the current test data and various fault data and normal data projection points. The method can effectively fuse multi-dimensional parameter association information, improves early fault sensitivity and fault traceability accuracy, and has the advantages of high diagnosis efficiency, strong adaptability and small required sample size.

The present application relates to the technical field of screw head size measurement, and specifically to a high-precision screw head size measurement method and system, the method comprising: obtaining the hardness factor of the screw head in each measurement process and the radial coordinate set at each rotation angle; obtaining the groove matching degree of each measurement process based on the difference between the fitting function of the radial coordinate set at each rotation angle of each measurement process and the ideal groove function of the screw head groove; obtaining the current groove risk coefficient of the measured screw based on the difference between the predicted value of the groove matching degree of the next measurement process and the average level of the groove matching degree of all current measurement processes, judging whether to trigger a closed-loop control mechanism, and then obtaining the size of the screw head groove. The present application improves the accuracy of screw head size measurement by analyzing the risk of screw assembly failure and adjusting the angular resolution during laser scanning.

The invention discloses a superconducting motor air gap winding inductance calculation method. The method comprises the steps that a plane polar coordinate system is established; determining a first radial coordinate and a first angular coordinate of the actual position of the first coil on the plane polar coordinate system; determining a first equivalent radial coordinate and a first equivalent angular coordinate of a first equivalent coil of the stator core corresponding to the first coil on the plane polar coordinate system; calculating the first self-excitation inductance of the first coil; calculating a second self-excitation inductance of the first equivalent coil based on the polar distance of the winding, the radius of the inner ring, the first radial coordinate, the first angular coordinate, and a second preset radial coordinate and a second preset angular coordinate of a second preset point on the planar polar coordinate system, and determining a separate excitation inductance of the stator core corresponding to the first coil based on the second self-excitation inductance; and determining the total self-inductance of the first coil based on the first self-excitation inductance and the separate-excitation inductance. Through the method, the calculation difficulty and the test cost of the winding self-inductance are reduced, and the calculation efficiency and the test efficiency are improved.

The invention relates to the field of crystal bar machining, in particular to a crystal bar rounding machining method, device and system and a storage medium. The outer circumferences of the crystal bars are scanned, three-dimensional coordinate data # imgabs0 # of N crystal bars are obtained, an initial contour is formed, N is the number of sampling points, xi is circumferential coordinates, yi is axial coordinates, and zi is radial coordinates; establishing a first function model of the initial contour and the target contour about the removal amount; and the removal amount is determined based on the obtained initial contour and the target contour, and the crystal bar is ground. By adjusting the parameters of the grinding wheel, the technical effect of improving the crystal bar rolling machining efficiency is achieved.

The invention relates to the field of design of circulating circles of hydraulic torque converters, in particular to an egg-shaped circulating circle design method and structure based on arc connecting lines. Comprising the following steps: calculating an overflowing area according to an outer diameter, an inner diameter and an overflowing coefficient of a circulating circle; calculating radial coordinates of an upper dead center and a lower dead center of the middle streamline according to the outer diameter, the inner diameter and the open area; the upper end arc radius, the lower end arc radius and the connection arc radius of the middle streamline are obtained, and whether the upper end arc radius, the lower end arc radius and the connection arc radius meet constraint conditions or not is verified; according to the design method, the streamline fluency of the circulating circle of the hydraulic torque converter can be improved, and then the transmission efficiency of the torque converter is improved.

The invention relates to the technical field of optical detection, in particular to a computer-generated hologram detection optical element processing error global inversion method, which comprises the following steps of: constructing a function relationship of an etching depth error with respect to a local grating period reciprocal and a radial coordinate as an initial etching depth error physical parameterization model; constructing a function relationship of a duty ratio error with respect to a local grating period reciprocal and a radial coordinate, and taking the function relationship as an initial duty ratio error physical parameterization model; obtaining a plurality of groups of actual measurement parameters; substituting multiple groups of actually measured parameters into the initial etching depth error physical parameterized model and the initial duty ratio error physical parameterized model to obtain an etching depth error physical parameterized model and a duty ratio error physical parameterized model; and substituting the full-aperture design data into the etching depth error physical parameterization model and the duty ratio error physical parameterization model to obtain etching depth error global distribution and duty ratio error global distribution. According to the invention, low-cost and high-precision error reconstruction can be realized at least.

The present application relates to the field of petrochemical industry and environmental protection technology, and particularly relates to a multi-core pipe parallel type main and auxiliary cavity cyclone oil-gas-water separation system and a separation method. The method comprises obtaining the time of flight of reflected sound waves and the Doppler power spectral density diagram of the cyclone separation main cavity; determining the axial coordinate sequence and the radial coordinate sequence based on the time of flight; determining whether the cyclone intensity of the cyclone separation main cavity is qualified based on the height difference of the cyclone liquid column; in response to the unqualified cyclone intensity, determining the adjustment direction of the gas phase outlet back pressure based on the cyclone gas core ratio; based on the qualified cyclone intensity, determining whether the cyclone load of the cyclone separation auxiliary cavity is up to standard based on the frequency spectrum width value, and in response to the unqualified cyclone load, prolonging the liquid phase separation time of the cyclone separation main cavity; based on the qualified cyclone load, determining whether there is a downward trend of the oil-gas-water separation efficiency, and enabling sound early warning when the oil-gas-water separation efficiency has a downward trend. The present application solves the problem of lag in sensing the state fluctuation of the oil-gas-water mixed liquid flow field.

The application provides a motor cover plate processing detection system and method based on visual recognition. The system collects the scattered light image and the reflected light image of the motor cover plate mounting hole area. The gray scale distribution characteristics of multiple concentric ring areas with the mounting hole center as the center are extracted in the scattered light image, and the ring gray scale distribution curve of the motor cover plate mounting hole area is determined by all the gray scale distribution characteristics. The radial coordinate constraint is performed on the ring gray scale distribution curve based on the geometric center coordinates of the motor cover plate mounting hole in the reflected light image, and the confidence ring gray scale distribution curve is generated. According to the confidence ring gray scale distribution curve and the qualified product gray scale distribution curve, the local variance abnormal section of the motor cover plate mounting hole existing hidden cracks is identified. The risk level of the motor cover plate mounting hole area existing hidden cracks is determined based on the local variance abnormal section. The technical scheme provided by the application can accurately detect the hidden cracks of the motor cover plate mounting hole area in the motor cover plate processing process.

The invention relates to the field of ophthalmology and is intended for testing the differential light sensitivity of the retina in different areas of the visual field using SAP perimetry. The present device comprises a unit for generating stimuli and establishing a testing procedure, a unit for displaying stimuli, and a feedback unit. The unit for generating stimuli and establishing a testing procedure generates diagnostic stimuli and establishes a topographical model of diagnostic stimuli which is based on a radial coordinate system and is a grid of lines that is formed by radial rays originating from the centre of the retina and by concentric circles, the centre of which is also situated at the centre of the retina. The stimuli are situated at the nodes of intersection of the radial rays with the concentric circles. The result of the invention is an increase in the reliability of testing of the differential light sensitivity of the retina of the eye using SAP perimetry, and an increase in the accuracy of detection of potential central visual field defects caused by depression of the light sensitivity of the retina.

An optical system (100) comprises focusing optics (10) configured for generating a focused light beam (B) at a focusing region (Fr). The focusing optics (10) comprise a set of spherical surfaces (11s, 12s) configured for introducing respective spherical aberrations into the focused light beam (B). The spherical aberration is adjusted to maximize the length (Lr) of the focal region (Fr). The system comprises at least one rotatable element (22) configured for redirecting the focused beam (B) in a variable beam direction. The gradient index 10 plate (31) is configured to receive the light beam from the rotatable element (22) and linearize a correlation of a radial coordinate (R) of the focused light beam in a target plane (P) with a rotation angle (delta) of the rotatable element (22). Advantageously, the gradient index plate is flat and thus does not introduce additional aberrations that may interfere with the focusing function.

A dressing tool for dressing a grinding worm for generating a pre-toothed workpiece is configured to rotate about a dressing rotation axis A and has a dressing flank surface that forms an annulus in a projection plane perpendicular to the dressing rotation axis of the dressing tool, the annulus defining a radial direction by a radial coordinate r and a circumferential direction by an angular coordinate φ. The dressing flank surface is provided with an intentionally generated dress flank modification, and the dress flank modification projected onto the annulus can be described as a two-dimensional Fourier series depending on the radial coordinate r and the angular coordinate φ, the Fourier series having at least one non-zero Fourier coefficient for a circumferential frequency component that is non-zero with respect to the angular coordinate.

The application discloses a proportional boundary perfect matched layer method and system for numerical simulation of infinite domain wave motion, and belongs to the field of dynamic analysis, which comprises the following steps: irregular geometric infinite domain in a Cartesian coordinate system is mapped to a regular geometric domain in a proportional boundary local coordinate system based on proportional boundary coordinate transformation; the radial coordinate of the regular geometric domain is analytically extended to a complex space by using a complex coordinate stretching function, and a proportional boundary perfect matched layer is obtained by truncation processing according to a preset thickness; based on the Galerkin weighted residual principle, corresponding equivalent integral expressions are constructed for the balance equation and the physical equation of the wave motion problem in the proportional boundary perfect matched layer; mixed displacement-stress finite element technology is adopted to perform finite element discrete processing on the displacement, stress and auxiliary variables of stress integral in the constructed equivalent integral expressions, so that mixed displacement-stress finite element equations of the proportional boundary perfect matched layer are obtained; the method has the advantages of simple boundary implementation and good adaptability of the boundary to the inner domain.

The application relates to the technical field of myopic lenses, and discloses a concentric ring alternating type outer ring defocusing method with first-order differential without discontinuity and a myopic lens. The method comprises the following steps: taking a lens center as an origin to establish a radial coordinate system, and dividing a lens surface into a central vision area and an outer ring defocusing area; M concentric ring circles of the outer ring defocusing area are divided into odd ring circles and even ring circles to form a ring belt structure with alternating wide and narrow rings; at a boundary position r m , constraint conditions of an mth ring circle and an m+1th ring circle in the ring belt structure are established, and a surface function of the m+1th ring circle is solved to form a ring defocusing surface with first-order differential without discontinuity. The application solves the optical aberration problem caused by first-order differential discontinuity at the junction of the ring belt in the prior art, adopts the design of first-order differential without discontinuity, and improves the optical performance and wearing comfort of the lens.

The application relates to a positron trapping system and method. The system comprises a trapping light source and a magnetic field deflection device; the trapping light source is used for forming a trapping electromagnetic field, so that the radial coordinate and radial momentum of a charged particle beam entering the trapping electromagnetic field are reduced and limited in a trapping area; the trapping light source is a left-handed circularly polarized Laguerre-Gaussian laser or a vector polarized light beam; there is a phase difference of pi / 2 between a radial electric field and a longitudinal electric field in the trapping electromagnetic field; the longitudinal electric field comprises a plurality of alternating longitudinal acceleration electric fields and longitudinal deceleration electric fields; the magnetic field deflection device is arranged in front of the trapping area and is used for separating positrons and electrons in the trapped charged particle beam to obtain a positron beam.

The application discloses a kind of neural radiation field volume super-resolution reconstruction methods for brain microhemorrhage detection.First, the low-resolution SWI body data of patient is acquired, sample is constructed by multi-scale coordinate sampling, including: local sparse coordinate sampling is carried out on slice to obtain the first-level coordinate set of hollow patch, then three-dimensional radial coordinate sampling is carried out to each coordinate to obtain the second-level coordinate set;Subsequently, a neural radiation field model is constructed, the sample is input into the model and combined with the spherical body rendering technique to generate the final rendering value;Then, according to the preset super-resolution multiple a interval selection slice constitutes training set and to be reconstructed set, after training the model, the to be reconstructed set is predicted, and finally the high-resolution SWI body data with the number of Z-axis slices increased is obtained.The application does not need paired training data, and the local structure reconstruction fidelity is significantly improved by the innovative sampling strategy, and the unique magnetic susceptibility characteristics of SWI can be accurately reconstructed, which effectively overcomes the shortcomings of the prior art.