A method and system for designing a lens with cumulative prism deviation

Abstract

The present application relates to the technical field of optical lens design, and discloses a lens design method and system for interlayer angle cumulative offset, comprising: constructing polar coordinate layered topology parameters; generating pre-optimization angle coefficients and each layer starting phase sequence through geometric shading conditions deduction and safety factor adjustment; constructing spiral point array coordinate set based on target line density method and polar coordinate rotation transformation; calling radial ray collision detection algorithm for shading evaluation, combining numerical gradient iterative optimization and rational number avoidance strategy to extract optimization angle coefficients of maximum corridor shading rate; generating lens processing coordinate file through smooth transition constraint verification and boundary clipping. The present application can effectively shield the radial visual corridor and improve the night vision quality of myopia prevention lenses.

Description

Technical Field

[0001] This invention relates to the field of optical lens design technology, and more specifically, to a lens design method and system for interlayer angle cumulative offset. Background Technology

[0002] Myopia control lenses achieve peripheral defocus control by incorporating microstructure arrays on the lens surface. These lenses are widely used in the field of myopia control for teenagers. The microstructure arrays typically employ a multi-layered circular layout, with each layer containing a different number of microstructure units. The arrangement of these layers directly affects the optical performance and visual quality of the lens.

[0003] Existing microstructure array layout schemes mainly include concentric ring layout, golden angle arrangement, and random layout. In the concentric ring layout, each layer of units is strictly aligned in the radial direction, forming a radial visual corridor that runs through the entire array; the golden angle arrangement uses a fixed interlayer offset angle of about 2.4 radians, which cannot be adapted to specific layers and unit sizes; although the random layout can break the radial alignment, the uncertainty of the unit position makes it difficult to guarantee optical uniformity.

[0004] The existing technologies have the following technical problems: the radial corridors of the concentric circle layout produce stray light radial inhomogeneity and visual interference in the nighttime point light source environment; the golden angle scheme may still have periodic alignment forming local visual corridors under a certain number of layers; the random layout causes random fluctuations in the intensity of defocus signals in different directions, affecting the uniformity and stability of myopia prevention and control effects. Summary of the Invention

[0005] This invention provides a lens design method and system for interlayer angular cumulative offset, which solves the technical problems of nighttime visual interference caused by radial visual corridors and periodic alignment and optical inhomogeneity in existing layout schemes under specific configurations.

[0006] This invention provides a lens design method for interlayer angle cumulative offset, comprising the following steps:

[0007] S1, obtain the diameter of the optical region of the lens and the nominal diameter of the microstructure unit, and construct the polar coordinate hierarchical topology parameters through geometric mapping and discretization;

[0008] S2, for polar coordinate layered topology parameters and microstructure unit diameter, generates pre-optimized angle coefficients and initial phase sequences for each layer through geometric occlusion condition derivation and safety factor adjustment strategy;

[0009] S3 integrates pre-optimized angle coefficients, initial phase sequences of each layer, and topology parameters. Based on the target line density method and polar coordinate rotation transformation mechanism, it constructs a set of spiral lattice coordinates and a sequence of points in each layer.

[0010] S4, the radial ray collision detection algorithm is called to evaluate the occlusion of the spiral lattice. Combined with numerical gradient iterative optimization and rational number avoidance strategy, the optimized angle coefficient and optimized lattice coordinate set that maximize the corridor occlusion rate are extracted.

[0011] S5 optimizes the angle coefficients, the set of optimized point coordinates, and the sequence of points in each layer. After smooth transition constraint verification and boundary clipping, the final lens processing coordinate file is generated.

[0012] In a preferred embodiment, S1 includes:

[0013] The effective area of ​​the microstructure array is determined according to the lens design requirements. The minimum radius is set as the radius of the central visible area of ​​the lens, and the maximum radius is determined by the diameter of the optical zone of the lens.

[0014] Determine the interlayer radial step size, which determines the radial distance between two adjacent annular layers. Its value matches the radial dimension of the nominal diameter of the microstructure unit. Calculate the total number of layers based on the array radial range and the interlayer radial step size. Calculate the number of reference circumference points based on the innermost layer circumference and the nominal diameter of the microstructure unit.

[0015] In a preferred embodiment, S2 includes:

[0016] The angle of the microstructure unit in the innermost layer is calculated. The angle is determined by the ratio of the arc length of the microstructure unit to the radius of the ring it is in. The theoretical minimum offset is derived. In order to make the microstructure units of the two adjacent layers have a shading effect on the radial projection, the interlayer angle offset causes the centers of the upper and lower layers of units to be misaligned in the radial direction. When the offset reaches half of the angle of the unit, the two layers of units begin to intersect effectively.

[0017] A safety factor is introduced, which is set according to the range of array layers and the uniformity of cell size. The initial angle accumulation factor is calculated by multiplying the safety factor by the theoretical minimum offset.

[0018] In a preferred embodiment, S2 further includes:

[0019] The safety factor is set according to the number of array layers. A first safety factor threshold is set for large-scale arrays within a preset first layer range, a second safety factor threshold is set for medium-scale arrays within a preset second layer range, and a third safety factor threshold is set for small-scale arrays within a preset third layer range. The first safety factor threshold is greater than the second safety factor threshold, and the second safety factor threshold is greater than the third safety factor threshold.

[0020] Based on a simplified theoretical model, the occlusion performance is estimated and the initial angle accumulation coefficient is adjusted. If the estimated number of occlusions corresponding to the current initial angle accumulation coefficient is less than the target number, the pre-optimized angle coefficient is obtained by adjusting it proportionally. The initial phase sequence of each layer is generated, and the initial phase of the i-th layer is equal to the product of the layer index and the pre-optimized angle coefficient.

[0021] In a preferred embodiment, S3 includes:

[0022] The linear density of layer 0 is calculated as a baseline. Linear density is defined as the number of points distributed per unit arc length. The baseline linear density is equal to the number of baseline points divided by the perimeter of layer 0. The target number of points for each layer is calculated. For layer i, in order to make the linear density of the corresponding layer close to the baseline linear density, the required number of points for the corresponding layer is obtained by rounding up the product of the baseline linear density and the perimeter of the corresponding layer.

[0023] In a preferred embodiment, S3 further includes:

[0024] For the j-th point on the i-th layer, calculate its angular coordinates on the corresponding layer's annulus. The angular coordinates consist of two parts: first, the base angle of the corresponding point in the circumference division to ensure uniform distribution of points within the same layer; second, the inter-layer angular offset introduced by the initial phase of the corresponding layer. Convert the polar coordinates to Cartesian coordinates. For a point located on the i-th layer, its Cartesian coordinates are obtained by trigonometric function transformation of the corresponding layer's radius and angular coordinates.

[0025] Traverse all layers and all points to generate a complete set of point coordinates. The set of point coordinates contains a total of coordinates for each layer, and each point contains its Cartesian coordinates and the index of the layer to which it belongs.

[0026] In a preferred embodiment, S4 includes:

[0027] A set of radial rays for occlusion detection is generated within a preset angle range, starting from the geometric center of the lens. The number of rays determines the angular resolution for occlusion rate evaluation. For each ray, a layer-by-layer collision detection is performed from the inner layer to the outer layer to determine the number of times the corresponding ray is blocked by microstructure units when passing through the multi-layer array. For the k-th ray and the i-th layer, the Euclidean distance between the radial position of the ray in the corresponding layer and the center of each unit in the corresponding layer is calculated. If the corresponding distance is less than the effective occlusion radius of the microstructure unit, it is determined that the ray is blocked by the corresponding unit in the corresponding layer.

[0028] The number of completely blocked rays is counted, and a threshold for the number of blocking times is set to define the minimum number of blocking times required to determine that a ray is completely blocked. The number of rays that meet the requirement of not less than the threshold number of blocking times is counted. The corridor blocking rate is calculated, which is defined as the proportion of completely blocked rays to the total number of rays.

[0029] In a preferred embodiment, S4 further includes:

[0030] Set the target occlusion rate. If the current occlusion rate is greater than or equal to the target occlusion rate, the current angle coefficient meets the requirements. Otherwise, start the gradient optimization iteration. Initialize the iteration parameters. Calculate the occlusion rate corresponding to the current angle coefficient in each iteration. Then calculate the gradient using the numerical difference method. Apply a preset perturbation amount to the current angle coefficient to calculate the perturbed occlusion rate. The gradient is equal to the difference between the two divided by the perturbation amount.

[0031] The adjustment direction of the angle coefficient is determined based on the difference between the target occlusion rate and the current occlusion rate, as well as the gradient. The angle coefficient is then updated. The current step size is calculated by dividing the initial step size by the decay factor to achieve adaptive decay. Boundary constraints are applied to the updated angle coefficient to determine the convergence condition. If the current occlusion rate has reached the target, or the change in the angle coefficient between two adjacent iterations is less than the preset convergence threshold, or the preset maximum number of iterations has been reached, the iteration is terminated.

[0032] In a preferred embodiment, S4 further includes:

[0033] After gradient optimization converges, a rational number avoidance test is performed. When the ratio of the cumulative inter-layer angle coefficient to twice pi is equal to a simple rational number, the cell positions in the spiral lattice are periodic. Every few layers in the angle space, a complete alignment occurs, forming equally spaced radial alignment bands, and the gaps between them constitute a visual corridor. The candidate angle coefficient is divided by twice pi to obtain a dimensionless ratio. The dimensionless ratio is compared one by one with all simple rational numbers whose denominators do not exceed the preset upper limit of the denominator.

[0034] If the absolute value of the difference between the ratio of a candidate angle coefficient to twice pi and a certain simple rational number is less than or equal to a preset judgment threshold, it is determined to fall into the neighborhood of rational numbers. A preset offset is applied to the candidate angle coefficient, and the occlusion rate is recalculated after the offset. If it still falls into the neighborhood, the offset continues until it passes the test. The criterion for passing the test is that the difference with all examined simple rational numbers is greater than the judgment threshold. The final optimized angle coefficient, the corresponding occlusion rate, and the set of optimized point coordinates corresponding to the optimized angle coefficient are recorded.

[0035] In a preferred embodiment, setting the target occlusion rate includes:

[0036] The target occlusion rate is set in different levels according to the application scenario of the lens. The first target occlusion rate threshold is set for applications with high requirements for nighttime vision, and the second target occlusion rate threshold is set for applications that are mainly used during the day. The first target occlusion rate threshold is greater than the second target occlusion rate threshold.

[0037] In a preferred embodiment, the step of extracting the optimization angle coefficients by combining numerical gradient iterative optimization further includes:

[0038] Boundary constraints are applied to the updated angle coefficients. The lower bound of the angle coefficients is determined based on the minimum effective interleaving requirement, and the upper bound of the angle coefficients is determined based on the optical area loss limit to avoid excessive overlap. The adaptive attenuation is calculated by dividing the current step size by the initial step size and the attenuation factor related to the number of iterations.

[0039] In a preferred embodiment, the smooth transition constraint verification and boundary clipping process includes:

[0040] For each coordinate point in the optimized dot matrix coordinate set, perform boundary mapping of the actual lens contour, check whether each coordinate point is within the boundary of the actual lens contour, and perform clipping and removal processing on coordinate points that exceed the actual contour boundary to adapt to the non-circular frame contour.

[0041] In a preferred embodiment, generating the final lens processing coordinate file further includes:

[0042] Optical parameters, including additional diopter, surface radius of curvature, and aspheric coefficient, are added to each microstructure unit in the coordinate file to generate a complete machining file that can be directly used in the manufacturing process.

[0043] In a preferred embodiment, a lens design system for interlayer angle cumulative offset is used to perform the steps in the lens design method for interlayer angle cumulative offset described above, including:

[0044] The topology parameter construction module is used to obtain the diameter of the optical region of the lens and the nominal diameter of the microstructure unit. After geometric mapping and discretization, the polar coordinate hierarchical topology parameters are constructed.

[0045] The angle coefficient pre-optimization module is used to generate pre-optimized angle coefficients and initial phase sequences for each layer based on the polar coordinate layered topology parameters and the nominal diameter of the microstructure unit, through geometric occlusion condition derivation and safety factor adjustment strategy.

[0046] The spiral lattice construction module is used to integrate pre-optimized angle coefficients, initial phase sequences of each layer, and polar coordinate layered topology parameters. Based on the target line density method and polar coordinate rotation transformation mechanism, it constructs a spiral lattice coordinate set and a sequence of points in each layer.

[0047] The occlusion assessment and optimization module is used to call the radial ray collision detection algorithm to assess the occlusion of the spiral lattice coordinate set. It combines numerical gradient iterative optimization and rational number avoidance strategy to extract the optimized angle coefficient and optimized lattice coordinate set that maximizes the corridor occlusion rate.

[0048] The coordinate file generation module is used to generate the final lens processing coordinate file by taking the optimized angle coefficients, optimized point matrix coordinate set and the number of points in each layer, verifying the smooth transition constraint and performing boundary clipping.

[0049] The beneficial effects of this invention are as follows:

[0050] By deriving geometric occlusion conditions and adjusting safety factors, pre-optimized angle coefficients are generated. The optimal inter-layer angle accumulation coefficient is determined by combining numerical gradient iterative optimization and rational number avoidance strategies. This causes systematic misalignment of each layer unit in the radial direction, which geometrically blocks the formation of radial visual corridors and avoids periodic alignment caused by the ratio of the angle coefficient to twice the pi falling into the neighborhood of simple rational numbers, thereby improving the visual quality of the lens in nighttime point light source environments.

[0051] The distribution of points in each layer is determined by the target line density method, and a spiral lattice is constructed by polar coordinate rotation transformation. This spiral quasi-symmetric layout maintains rotational symmetry mathematically, ensuring the uniformity of defocus signals in all directions. At the same time, the corridor occlusion rate is quantified by the radial ray collision detection algorithm, providing a reliable evaluation mechanism and optimization basis for the optical design of myopia control lenses. Attached Figure Description

[0052] Figure 1 This is a flowchart of a lens design method for interlayer angle cumulative offset according to the present invention;

[0053] Figure 2 This is a flowchart of a lens design method for interlayer angle cumulative offset according to the present invention. Detailed Implementation

[0054] The subject matter described herein will now be discussed with reference to exemplary embodiments. It should be understood that these embodiments are discussed only to enable those skilled in the art to better understand and implement the subject matter described herein, and changes may be made to the function and arrangement of the elements discussed without departing from the scope of this specification. Various processes or components may be omitted, substituted, or added as needed in the examples. Furthermore, some features described in the examples may be combined in other examples.

[0055] At least one embodiment of the present invention discloses a lens design method for interlayer angle cumulative offset, such as Figures 1 to 2 As shown, it includes the following steps:

[0056] S1, obtain the diameter of the optical region of the lens and the nominal diameter of the microstructure unit, and construct the polar coordinate hierarchical topology parameters through geometric mapping and discretization;

[0057] Diameter of the optical zone of the receiving lens and nominal diameter of microstructure unit Using geometric mapping and discretization methods, polar coordinate layered topology parameters (total number of layers N, number of reference points) are obtained. Radial step size (and the radius sequence of each layer).

[0058] Specifically, the effective area of ​​the microstructure array is determined according to the lens design requirements; the minimum radius of the array... This value is typically set to the radius of the central visual field of the lens. This visual field is used to ensure clear central vision. The default value is set to a range of one-tenth to one-fifth of the optical zone diameter. This embodiment selects... mm; maximum radius of the array The value is determined by the diameter of the optical zone of the lens, with the default value set to half the diameter of the optical zone. This embodiment applies to standard lenses. mm take mm;

[0059] Determine the interlayer radial step size This step size determines the radial distance between two adjacent ring layers, and its value should match the radial dimension of the microstructure unit; in this embodiment, To ensure sufficient overlap between adjacent layers to maintain optical continuity, while avoiding excessive density to prevent fabrication difficulties, the default value for the interlayer radial step size is set to 0.85 times the element diameter. In this embodiment, the step size is... mm;

[0060] Calculate the total number of layers N; based on the array radial range and inter-layer step size, the total number of layers is... Substituting the values ​​into the equations is worthwhile. Layers; the radius sequence of each layer is as follows ,in to The corresponding radii are 1.5mm, 2.35mm, 3.2mm, 4.05mm, 4.9mm, 5.75mm, 6.6mm, and 7.45mm, respectively.

[0061] Calculate the number of points on the reference circle This parameter determines the number of microstructural units distributed on the 0th layer annulus; the circumference of the 0th layer is determined by... The number of reference points is calculated by dividing the perimeter by the unit diameter and rounding down to ensure adequate spacing between units. In this embodiment, the perimeter of layer 0 is approximately 9.42 mm, and the number of reference points... This baseline number serves as a reference value for calculating the number of points in subsequent layers.

[0062] The above calculations completed the initialization of the polar coordinate hierarchical topology, yielding the total number of layers. Number of benchmark points Radial step size mm, the radius sequence of each layer; these parameters constitute the basic topological framework for the subsequent generation of the spiral lattice.

[0063] S2, for polar coordinate layered topology parameters and microstructure unit diameter, generates pre-optimized angle coefficients and initial phase sequences for each layer through geometric occlusion condition derivation and safety factor adjustment strategy;

[0064] Receive topology parameters (total number of layers N, number of reference points) output by S1 (radius sequence of each layer) and nominal diameter of microstructural unit By employing the derivation of geometric occlusion conditions and the adjustment of safety factors, the pre-optimized angle coefficients are obtained. and the initial phase sequence of each layer;

[0065] Calculate the angular range of the microstructural unit in the innermost layer; the angular range occupied by a single microstructural unit on the annulus is determined by the ratio of the arc length of the unit to the radius of the annulus; for the innermost layer radius... mm, unit diameter mm configuration, angle subtended The radius is approximately 38.2 degrees.

[0066] Derivation of the theoretical minimum offset: To ensure that adjacent microstructural units create an occlusion effect in the radial projection, the interlayer angular offset must be large enough to cause radial misalignment between the centers of upper and lower unit layers. According to geometric relationships, when the offset reaches half the unit angle, the two layers begin to effectively interlock; therefore, the theoretical minimum offset... radian;

[0067] However, setting α solely based on the theoretical minimum may result in insufficient shading, especially in the outer ring where the relative shading efficiency between cells can vary. To ensure effective shading across all layers, a safety factor β is introduced, which is set based on the array's layer count and the uniformity of cell size. For large-scale arrays with 50 to 80 layers, the default safety factor can be set to 1.5; for medium-scale arrays with 10 to 30 layers, the default safety factor can be set to 1.3; and for small-scale arrays with fewer than 10 layers, the default safety factor can be set to 1.2. The array in this embodiment has 8 layers, making it a small-scale demonstrative array, and the default safety factor is set to 1.2 to obtain a moderate shading safety margin.

[0068] Calculate the initial angle accumulation factor by multiplying the safety factor by the theoretical minimum offset. Radius; to avoid blind optimization in subsequent iterations, occlusion performance is predicted and initial values ​​are adjusted based on a simplified theoretical model; the simplified model used is: when a radial ray passes through a layer, if the angle difference between the center of a cell and the ray on that layer is less than half of the angle subtended by the cells of that layer, it is considered occluded. For phase offsets of... The spiral layout, when This ensures that any radial ray is blocked at least N / 2 times when passing through N layers; based on the prediction results, if the current If the estimated number of occlusions is less than the target number, the pre-optimized angle coefficient is obtained by adjusting it proportionally. This implementation method has been estimated and adjusted. radian;

[0069] Generate the initial phase sequence for each layer, for the first layer... Layers and their initial phase In this embodiment, the initial phases of each layer are 0, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, and 3.5 radians, respectively. This phase sequence represents the cumulative effect of the spiral rotation between layers.

[0070] Through the above processing, the initial value calculation and pre-optimization of the inter-layer angle accumulation coefficient were completed, and the pre-optimized angle coefficient was output. The radian and the corresponding initial phase sequence of each layer lay the parameter foundation for the subsequent precise lattice construction;

[0071] S3 integrates pre-optimized angle coefficients, initial phase sequences of each layer, and topology parameters. Based on the target line density method and polar coordinate rotation transformation mechanism, it constructs a set of spiral lattice coordinates and a sequence of points in each layer.

[0072] Receive the pre-optimized angle coefficients output by S2 And the initial phase sequence of each layer, and the topology parameters (number of reference points) output by S1. The distribution of points in each layer is determined by the target line density method, and the coordinate set of the spiral lattice and the number sequence of points in each layer are obtained by combining polar coordinate rotation transformation and coordinate system transformation methods.

[0073] Determine the target line density and calculate the number of points in each layer; calculate the line density of layer 0 as the baseline; line density is defined as the number of points distributed per unit arc length, the perimeter of layer 0 is the product of the radius of layer 0 and 2π, and the number of points is the baseline number of points. Therefore, baseline density equal Divide by the perimeter of layer 0; in this embodiment, the perimeter of layer 0 is approximately 9.42 mm, the number of reference points is 9, and the reference line density is calculated as follows: The linear density per millimeter serves as the target density standard that each layer should maintain.

[0074] Calculate the number of target points in each layer; for the i-th layer, its perimeter is... To make the linear density of this layer close to The number of points required for this layer is according to The product of the product with the perimeter of the layer is rounded up; in this embodiment, the calculated sequence of points for each layer is 9, 14, 19, 24, 29, 34, 39, 44;

[0075] Generate the polar and Cartesian coordinates of the spiral lattice; for the... The first layer For each point, calculate its angular coordinates on the annular layer; the angular coordinates of this point... It consists of two parts: one is the basic angle of the point in the equal division of the circle. To ensure uniform distribution of points within the same layer, and secondly, the initial phase of that layer. Introducing interlayer angle offset; adding the two parts yields... All points in layer 1 have changed relative to layer 0. The overall rotational offset of the arc reflects the characteristics of spiral accumulation;

[0076] Convert polar coordinates to Cartesian coordinates; for a layer located at angle , The point has Cartesian coordinates as and ; Traverse all layers and all points and perform the above transformation to generate a complete set of point coordinates; This set is a list containing the total number of coordinates of points in each layer, and each point contains its Cartesian coordinates and the index of the layer to which it belongs; The total number of points generated in this embodiment is 212;

[0077] Technical characteristics of the spiral lattice analyzed and compared with existing schemes; observing the coordinate distribution of the generated lattice reveals a spiral arrangement; taking the lens center as the origin and observing outward along any fixed angle, the points in each layer are not on the same radial line, but shift angularly with increasing radius, forming a spiral arm-like array of points; this spiral arrangement is formed by the linear accumulation of the initial phases of each layer, mathematically maintaining rotational symmetry, meaning that after rotating the entire lattice around the center by any angle, its topological structure is the same as the original lattice; this spiral quasi-symmetric lattice is fundamentally different from several existing typical layout schemes; compared with the traditional concentric ring layout, the traditional layout... In traditional arrays, each layer of units is strictly aligned radially, forming a radial visual corridor running through the entire array. However, this invention uses cumulative interlayer angular offset to create a systematic misalignment of each layer of units in the radial direction, geometrically blocking the formation of this corridor. Compared to the known golden angle arrangement schemes in botany and optics, which use a fixed interlayer offset angle of approximately 2.4 radians (approximately 137.5 degrees), this fixed value is indiscriminate across all lens specifications and cannot be adapted to specific layer counts and unit sizes. In certain layer configurations, localized radial alignment may still occur. For example, when the lens design has 34 layers, the golden angle of 2.4 radians... The ratio is approximately 0.382, which is close to a simple rational number. That is, 0.375. The deviation between the two is much smaller than the engineering safety threshold. At this time, the units of the 8th, 16th, 24th and 32nd layers will show quasi-periodic alignment in the radial direction, forming 8 equally spaced local visual corridors. The angular accumulation coefficient α of this invention is dynamically determined for specific lens parameters through the gradient optimization process of S4, and combined with rational number avoidance test to ensure that the coefficient is consistent with the actual lens parameters. The ratio does not fall into the neighborhood of simple rational numbers, thus ensuring the effectiveness of corridor shading under any number of layers. Compared with a completely random layout, although the random layout can break radial alignment, the uncertainty of its unit position makes it difficult to guarantee optical uniformity, and the intensity of defocus signal in different directions fluctuates randomly. The spiral quasi-symmetric layout of the present invention maintains rotational symmetry mathematically. This characteristic ensures the uniformity of defocus signal in all directions, meeting the strict requirements of myopia control lenses for optical uniformity in all directions.

[0078] Through the above processing, the construction of the spiral quasi-symmetric lattice model was completed, and the spiral lattice coordinate set (containing 212 points in this embodiment) and the point number sequence of each layer (9, 14, 19, 24, 29, 34, 39, 44) were output, providing a geometric basis for subsequent occlusion evaluation and optimization.

[0079] S4, the radial ray collision detection algorithm is called to evaluate the occlusion of the spiral lattice. Combined with numerical gradient iterative optimization and rational number avoidance strategy, the optimized angle coefficient and optimized lattice coordinate set that maximize the corridor occlusion rate are extracted.

[0080] Receive the spiral lattice coordinate set and the sequence of points in each layer output by S3, as well as the pre-optimized angle coefficients output by S2. By employing a radial ray collision detection algorithm and an iterative optimization method based on numerical gradients, the optimal angle coefficient that maximizes the corridor occupancy rate is obtained. and the corresponding optimized set of lattice coordinates;

[0081] Generate a radial detection ray and configure its geometric parameters; starting from the geometric center of the lens, within the angular range of 0 to... A set of radial rays is generated internally for occlusion detection; the number of rays, K, determines the angular resolution of the occlusion rate assessment. Based on the angular scale of the microstructural units, to ensure detection accuracy, K should be no less than three times the product of the number of reference points and the total number of layers; the default value for the number of rays is set to... The corresponding angular resolution is 1 degree; the first The angle of a ray is the product of its index and the angular interval, where the ray index is... From 0 to Each ray is represented as a semi-straight line extending infinitely outward from the origin along the corresponding angle;

[0082] Perform ray collision detection and count the number of occlusions; for each ray, perform layer-by-layer collision detection from the inner layer to the outer layer to determine the number of times the ray is occluded by microstructure units while traversing the multilayer array; for the ... ray and the first Layer, obtain the center coordinates of all microstructural units in that layer, and filter out the layers with the index as... All points; rays In the The radial position corresponding to the layer is This position is represented in Cartesian coordinates as ; Calculate the relationship between this point and the first The first layer Euclidean distance between the centers of each unit If this distance is less than the effective shielding radius of the microstructure unit Then determine the ray at the th . The layer is blocked by this unit; the effective shading radius takes into account the physical size and optical range of the unit, and the default value is set to 1.2 times the nominal radius of the unit. In this embodiment, it is taken as... mm; traversing the first For all elements in the layer, determine if at least one element satisfies the occlusion condition; if so, determine that the ray is in the [th]th layer. When a layer is obscured, increment the obscuration count for that layer by one, up to the total obscuration count counter. If it does not exist, then the ray is determined to be in the th position. The layer is not occluded; repeat the above detection for all N layers to obtain the th layer. Total number of times a ray is blocked as it passes through the entire array ;

[0083] Count the number of rays that are completely blocked; set a threshold for the number of blocking attempts. Define the minimum number of blocking operations required to completely block a ray. Considering the interlayer gaps and finite element dimensions, the standard for complete blocking can be appropriately relaxed. The default threshold for the number of blocking operations is set to be 0.75 times the total number of layers, rounded up. In this implementation, it is 6. Traverse all K rays and count the number of rays that satisfy the requirement of not less than the threshold number of blocking operations, denoted as . ;

[0084] Calculate the corridor shading rate; shading rate Defined as the proportion of rays that are completely blocked out of the total number of rays, i.e. This metric quantifies the array's ability to occlude radial visual corridors; the closer η is to 1, the more fully the corridor is occluded.

[0085] Iterative gradient optimization of angle coefficients and rational number avoidance; setting target occlusion rate This value is determined based on the visual quality requirements of the application scenario. For applications with high requirements for nighttime visual scenes, the default value can be set to 0.92, and for applications that are mainly daytime, the default value can be set to 0.88. In this embodiment, the default value for target occlusion rate is set to... If the current occlusion rate is greater than or equal to the target value, then the current α value meets the requirements, and we directly proceed to the rational number avoidance test; otherwise, we start the gradient optimization iteration.

[0086] The gradient optimization iterative process is as follows; initialize the iteration parameters: current iteration number. Current α value Equal to S2 pre-optimization The value is 0.5 radians in this implementation; the initial step size is determined based on the effective range of α and the target accuracy, with a default value of 0.05 radians; the maximum number of iterations is used to prevent infinite loops, with a default value of 20; in each iteration, the current... Corresponding occlusion rate Then, the gradient is calculated using the numerical difference method: in Apply a small perturbation to the base layer, with the default perturbation value set to 0.01 radians, and calculate the occlusion rate and gradient after the perturbation. It equals the difference between the two divided by the disturbance; based on the difference between the target occlusion rate and the current occlusion rate. and gradient Determine the adjustment direction of α: when Increase α when it is positive, decrease α when it is negative; update the value of α using the formula , where the sign function This indicates that the sign of the extracted value is either positive or negative 1, and the current step size is... Divide the initial step size by the decay factor The calculation achieves adaptive attenuation; the updated α value is subject to boundary constraints, with a lower bound default value of 0.2 radians based on the minimum effective interleaving requirement and an upper bound default value of 0.8 radians based on the limitation of optical area loss to avoid excessive overlap; the convergence condition is determined as follows: if the current occlusion rate has reached the target, or the change in α value between two adjacent iterations is less than the convergence threshold (the threshold default value is set to 0.001 radians), or the maximum number of iterations has been reached, then the iteration is terminated; otherwise, the next round continues; this implementation converges after several iterations, obtaining candidate values. radian;

[0087] After gradient optimization converges, a rational number avoidance test is performed; the physical basis of this test is: when the inter-layer angle cumulative coefficient α and... The ratio is exactly equal to some simple rational number. At that time, the positions of the elements in the spiral lattice have strict periodicity, after... After the first layer The angular positions of all units on the upper layer completely coincide with those of the 0th layer, and in the angular space, every... The layers then achieve a complete alignment, forming Equally spaced radial alignment bands, the gaps between which form new visual corridors; even if this ratio is only close to rather than exactly equal to a simple rational number, as long as the deviation is small enough, the above-mentioned periodic alignment effect will still approximately appear within a finite number of layers; this phenomenon has particular harmfulness in the long-term wearing scenario of myopia control lenses for teenagers: myopia control lenses need to continuously play a defocusing role in the wearer's all-day natural visual environment. If there is a periodic corridor direction, the microlens defocus signal in that direction will be weakened by the corridor gap, resulting in a systematic deficiency of peripheral retinal defocus in a specific direction. Long-term accumulation may affect the uniformity and stability of myopia control effect. Therefore, rational number avoidance is a key technical measure to ensure that the myopia control function is uniform and effective in all directions; the test method is: divide the candidate value by The dimensionless ratio is obtained, and this ratio is compared one by one with all simple rational numbers whose denominators do not exceed a specified upper limit. The default upper limit for the denominator is set to 20, and the default threshold is set to 0.005. The determination of this threshold takes into account the parameter range in actual lens design where the number of layers usually does not exceed 80 and the denominator of rational numbers does not exceed 20. Within this range, the cumulative phase error corresponding to this threshold is still much smaller than the critical condition for periodic alignment, and can be used as an engineering safety threshold. If the candidate value is compared with... If the absolute value of the ratio of a candidate value to a simple rational number is less than or equal to a threshold, the candidate value is considered to fall within the neighborhood of rational numbers. An offset is then applied to the candidate value, with the default offset range being 0.01 to 0.03 radians. The offset direction is preferably positive; if the occlusion rate drops below the target value or α exceeds the upper bound after a positive offset, a negative offset is used. The occlusion rate is recalculated after the offset. If the candidate value still falls within the neighborhood, the offset continues until it passes the test. The criterion for passing the test is that the difference between the candidate value and all examined simple rational numbers is clearly greater than the threshold. In this embodiment, the candidate value... radians divided by The difference between the result and a certain simple rational number falls within the judgment threshold range, and after rational number avoidance adjustment, the result is obtained. The radian value passed the test; the shading rate was recalculated to be 0.911, which meets the target requirement.

[0088] Output the optimization results; record the final optimized angle coefficients. Corresponding occlusion rate And the The corresponding optimized set of point coordinates; the final implementation of this method radian, ;

[0089] The spatial distribution uniformity of the shading performance was analyzed. In one embodiment of the present invention, in order to further understand the distribution characteristics of the shading effect in different radial directions, the optimized dot matrix was subjected to regional shading rate analysis. The generated rays were uniformly divided into several sectors according to the angle, and the shading rate of each sector was calculated. In this embodiment, the standard deviation of the shading rate of each sector was about 0.023, indicating that the shading performance of the spiral layout was relatively uniform in all directions, and there were no weak shading areas in specific directions, which verified the advantage of quasi-symmetry.

[0090] Through the above processing, a quantitative assessment of the visual corridor's occlusion and gradient optimization of the angle coefficients were completed, outputting optimized angle coefficients that meet the occlusion rate requirements. The set of radians and corresponding optimized lattice coordinates, which has achieved effective geometric blocking of the radial visual corridor;

[0091] S5 optimizes the angle coefficients, optimizes the set of point coordinates and the sequence of points in each layer, and generates the final lens processing coordinate file after smooth transition constraint verification and boundary clipping.

[0092] Optimized angle coefficients received from S4 output The corresponding optimized set of point coordinates, and the point number sequence of each layer output by S3, are used to verify the smooth transition constraint between adjacent layers and the boundary clipping method to obtain the final lens processing coordinate file;

[0093] Verify the smooth transition constraint of the number of points in adjacent layers; check whether the sequence of point numbers in each layer output by S3 satisfies the smooth transition constraint of adjacent layers; for the layer( ), calculate the rate of change of its number of points relative to the previous layer. Set the maximum allowable rate of change, with a default value of 0.6 to constrain the jump in the number of points between adjacent layers; traverse each layer, if the rate of change of a certain layer exceeds the maximum allowable value, then the number of points in that layer is growing too fast, and the number of points in that layer needs to be adjusted to the result of rounding up the product of the number of points in the previous layer and the upper limit of the rate of change. After adjustment, re-execute the lattice generation and optimization process of S3.2 and S4; if the adjustment of the number of points results in insufficient line density in a certain layer, it can be compensated by finely adjusting the unit size of that layer in the tangential direction; in this embodiment, the number of points in each layer is 9, 14, 19, 24, 29, 34, 39, 44. After checking, the rate of change of all adjacent layers meets the smooth transition constraint and no adjustment is required;

[0094] Generate the final processing coordinate file and attach process information; generate the lens processing coordinate file; the file is organized in a standard data format (such as CSV or GDSII), with each line recording the information of a microstructure unit, including the x-coordinate and y-coordinate of the unit center, the layer index, the nominal diameter of the unit, and other parameters; directly convert the optimized lattice coordinates output by S4 into the processing coordinate format. The coordinate file generated in this embodiment contains complete information of 212 units;

[0095] Perform boundary mapping of the actual lens contour; since the actual lens may have a non-circular frame contour, it is necessary to map the circular array area into the actual contour; for each point in the coordinate file, check whether it is within the boundary of the actual lens contour, if it exceeds the boundary, remove the point from the coordinate file to complete the boundary clipping.

[0096] Additional processing information; depending on the type of microstructure unit (microlens or defocus ring), add the optical parameters of the unit to the coordinate file, such as additional diopter, surface radius of curvature, aspheric coefficient, etc.; in this embodiment, the default value of the additional diopter of all microlens units is set to positive 3.0D;

[0097] Through the above S5 processing, the smooth transition constraint verification and the engineering encapsulation of the coordinate file are completed, and the final lens processing coordinate file is output. This file contains the complete geometric definition and optical parameters of all microstructure units after interlayer angle cumulative offset optimization and boundary clipping, which can be directly used in downstream manufacturing processes such as CNC milling, injection mold making or photolithography mask design.

[0098] In one embodiment of the present invention, a detailed process of the lens design method and system for interlayer angle cumulative offset described above is provided:

[0099] S1 receives the optical region diameter and nominal diameter of the microstructure unit of the receiving lens, and outputs polar coordinate layered topology parameters (total number of layers N, number of reference points). (Radial step size, radius sequence of each layer), these parameters are the geometric basis of the entire design;

[0100] S2 receives the topology parameters and cell diameter output by S1, and derives the pre-optimized angle coefficients through geometric occlusion conditions. And the initial phase sequence of each layer, which directly determines the intensity of the spiral offset;

[0101] S3 receives the output of S2. The initial phase sequence of each layer and the topology parameters output by S1 are used to determine the point distribution of each layer according to the target line density in order to achieve edge density compensation. Then, polar coordinate transformation is performed to generate a set of point coordinates and a sequence of points of each layer. This point matrix is ​​the geometric entity for subsequent evaluation.

[0102] S4 receives the set of matrix coordinates and the sequence of points for each layer output by S3, as well as the output of S2. The corridor occlusion performance is quantified using ray collision detection, and the α value is adjusted to [value missing] through gradient optimization and rational number avoidance. Output the optimized angle coefficients and the corresponding set of optimized point coordinates; the optimized... The data is then fed back to S3 to reconstruct the lattice, forming a computational closed loop to ensure the effectiveness of the optimization results under the final point configuration.

[0103] S5 receives the output of S4. The system optimizes the set of lattice coordinates and the sequence of points for each layer output by S3, performs smooth transition constraint verification and engineering processing, and outputs the final lens processing coordinate file. The inputs and outputs of each step are strictly connected to ensure the logical integrity and feasibility of the technical solution.

[0104] In one embodiment of the present invention, an application example is provided, which focuses on the application of myopia prevention and control lenses for teenagers. The design was verified for daily wear scenarios of school children aged 8 to 14. A typical configuration with a lens optical zone diameter of 16 mm, a microstructure unit diameter of 1.0 mm, and an array coverage radius ranging from 1.5 mm to 8 mm was selected, and the design was completed using the method of the present invention.

[0105] Examples of the obtained initial design parameters are shown in Table 1:

[0106] Table 1, Examples of Initial Design Parameters

[0107]

[0108] Table 2 summarizes the topology parameters obtained after S1 initialization and the number of points in each layer after density calculation using S3.

[0109] Table 2, Example of Summary of Points for Each Layer

[0110]

[0111] The number of points in each layer in the table is S3, based on the target line density. The calculated point distribution reflects the effect of the edge density compensation mechanism in increasing the number of points in the outer ring with the circumference.

[0112] After completing the design using the method of this invention, an optimized inter-layer angle cumulative coefficient is obtained. The arc corresponding to the corridor occlusion rate reaches 0.911. The design parameters were applied to the lens manufacturing and the performance was verified by the optical testing platform. Under simulated nighttime point light source illumination conditions, compared with the traditional concentric circle layout lens, the lens designed in this invention improved the stray light radial non-uniformity index of retinal imaging and improved the contrast sensitivity, verifying the effectiveness of the interlayer angle cumulative offset technology for radial visual corridor occlusion.

[0113] To further verify the applicability of the method of the present invention on a real-world product scale, a second embodiment is provided; the lens optical area diameter is selected as 16 mm, the nominal diameter of the microstructure unit is 0.9 mm, and the interlayer radial step size is [not specified]. mm, array coverage radius mm to Large-scale configuration of mm; execute S1 initialization to obtain the total number of layers. Layers, number of reference points ; Perform S2 to calculate the angle of the 0th layer. radians, theoretical minimum offset The array has 55 layers, which is considered a large-scale array. The default safety factor is set to [value missing]. Initial angle cumulative coefficient The radian, after being estimated and adjusted, yields the pre-optimized angle coefficient. Radius; Execute S3 to calculate the number of points in each layer based on the target line density and generate a lattice; Execute S4 for gradient optimization, and after several iterations, converge to obtain candidate values. After rational number avoidance testing and adjustment, the final value is obtained. Radiance, corresponding to occlusion rate The target requirements are met; S5 smoothing constraint verification and coordinate file generation are performed, with a total of approximately 3180 elements; this second embodiment can still stably converge to the effective angle accumulation coefficient under a large-scale array configuration, verifying the scalability and engineering practical value of the method of the present invention.

[0114] The embodiments of the present invention have been described above. However, the embodiments are not limited to the specific implementation methods described above. The specific implementation methods described above are merely illustrative and not restrictive. Those skilled in the art can make more equivalent embodiments under the guidance of the present embodiments, and all of them are within the protection scope of the present embodiments.

Claims

1. A lens design method for interlayer angle cumulative offset, characterized in that, Includes the following steps: S1, obtain the diameter of the optical region of the lens and the nominal diameter of the microstructure unit, and construct the polar coordinate hierarchical topology parameters through geometric mapping and discretization; S2, for polar coordinate layered topology parameters and microstructure unit diameter, generates pre-optimized angle coefficients and initial phase sequences for each layer through geometric occlusion condition derivation and safety factor adjustment strategy; calculates the angle subtended by the microstructure unit in the innermost layer, which is determined by the ratio of the arc length of the microstructure unit to the radius of the ring it is in; derives the theoretical minimum offset, in order to make the microstructure units of adjacent layers produce an occlusion effect on the radial projection, the interlayer angle offset causes the centers of the upper and lower layer units to be misaligned in the radial direction, and when the offset reaches half of the unit angle subtended, the two layers of units begin to effectively intersect; A safety factor is introduced, which is set according to the range of array layers and the uniformity of cell size. The initial angle accumulation factor is calculated by multiplying the safety factor by the theoretical minimum offset. S3 integrates pre-optimized angle coefficients, initial phase sequences of each layer, and topology parameters. Based on the target line density method and polar coordinate rotation transformation mechanism, it constructs a set of spiral lattice coordinates and a sequence of points in each layer. S4, the radial ray collision detection algorithm is called to evaluate the occlusion of the spiral lattice. Combined with numerical gradient iterative optimization and rational number avoidance strategy, the optimized angle coefficient and optimized lattice coordinate set that maximize the corridor occlusion rate are extracted. S5 optimizes the angle coefficients, the set of optimized point coordinates, and the sequence of points in each layer. After smooth transition constraint verification and boundary clipping, the final lens processing coordinate file is generated.

2. The lens design method for interlayer angle cumulative offset according to claim 1, characterized in that, S1 includes: The effective area of ​​the microstructure array is determined according to the lens design requirements. The minimum radius is set as the radius of the central visible area of ​​the lens, and the maximum radius is determined by the diameter of the optical zone of the lens. Determine the interlayer radial step size, which determines the radial distance between two adjacent annular layers. Its value matches the radial dimension of the nominal diameter of the microstructure unit. Calculate the total number of layers based on the array radial range and the interlayer radial step size. Calculate the number of reference circumference points based on the innermost layer circumference and the nominal diameter of the microstructure unit.

3. The lens design method for interlayer angle cumulative offset according to claim 1, characterized in that, S2 further includes: The safety factor is set according to the number of array layers. A first safety factor threshold is set for large-scale arrays within a preset first layer range, a second safety factor threshold is set for medium-scale arrays within a preset second layer range, and a third safety factor threshold is set for small-scale arrays within a preset third layer range. The first safety factor threshold is greater than the second safety factor threshold, and the second safety factor threshold is greater than the third safety factor threshold. Based on a simplified theoretical model, the occlusion performance is estimated and the initial angle accumulation coefficient is adjusted. If the estimated number of occlusions corresponding to the current initial angle accumulation coefficient is less than the target number, the pre-optimized angle coefficient is obtained by adjusting it proportionally. The initial phase sequence of each layer is generated, and the initial phase of the i-th layer is equal to the product of the layer index and the pre-optimized angle coefficient.

4. The lens design method for interlayer angle cumulative offset according to claim 1, characterized in that, S3 includes: The linear density of layer 0 is calculated as a baseline. Linear density is defined as the number of points distributed per unit arc length. The baseline linear density is equal to the number of baseline points divided by the perimeter of layer 0. The target number of points for each layer is calculated. For layer i, in order to make the linear density of the corresponding layer close to the baseline linear density, the required number of points for the corresponding layer is obtained by rounding up the product of the baseline linear density and the perimeter of the corresponding layer.

5. The lens design method for interlayer angle cumulative offset according to claim 4, characterized in that, S3 further includes: For the j-th point on the i-th layer, calculate its angular coordinates on the corresponding layer's annulus. The angular coordinates consist of two parts: first, the base angle of the corresponding point in the circumference division to ensure uniform distribution of points within the same layer; second, the inter-layer angular offset introduced by the initial phase of the corresponding layer. Convert the polar coordinates to Cartesian coordinates. For a point located on the i-th layer, its Cartesian coordinates are obtained by trigonometric function transformation of the corresponding layer's radius and angular coordinates. Traverse all layers and all points to generate a complete set of point coordinates. The set of point coordinates contains a total of coordinates for each layer, and each point contains its Cartesian coordinates and the index of the layer to which it belongs.

6. The lens design method for interlayer angle cumulative offset according to claim 1, characterized in that, S4 includes: A set of radial rays for occlusion detection is generated within a preset angle range, starting from the geometric center of the lens. The number of rays determines the angular resolution for occlusion rate evaluation. For each ray, a layer-by-layer collision detection is performed from the inner layer to the outer layer to determine the number of times the corresponding ray is blocked by microstructure units when passing through the multi-layer array. For the k-th ray and the i-th layer, the Euclidean distance between the radial position of the ray in the corresponding layer and the center of each unit in the corresponding layer is calculated. If the corresponding distance is less than the effective occlusion radius of the microstructure unit, it is determined that the ray is blocked by the corresponding unit in the corresponding layer. The number of completely blocked rays is counted, and a threshold for the number of blocking times is set to define the minimum number of blocking times required to determine that a ray is completely blocked. The number of rays that meet the requirement of not less than the threshold number of blocking times is counted. The corridor blocking rate is calculated, which is defined as the proportion of completely blocked rays to the total number of rays.

7. The lens design method for interlayer angle cumulative offset according to claim 6, characterized in that, S4 further includes: Set the target occlusion rate. If the current occlusion rate is greater than or equal to the target occlusion rate, the current angle coefficient meets the requirements. Otherwise, start the gradient optimization iteration. Initialize the iteration parameters. Calculate the occlusion rate corresponding to the current angle coefficient in each iteration. Then calculate the gradient using the numerical difference method. Apply a preset perturbation amount to the current angle coefficient to calculate the perturbed occlusion rate. The gradient is equal to the difference between the two divided by the perturbation amount. The adjustment direction of the angle coefficient is determined based on the difference between the target occlusion rate and the current occlusion rate, as well as the gradient. The angle coefficient is then updated. The current step size is calculated by dividing the initial step size by the decay factor to achieve adaptive decay. Boundary constraints are applied to the updated angle coefficient to determine the convergence condition. If the current occlusion rate has reached the target, or the change in the angle coefficient between two adjacent iterations is less than the preset convergence threshold, or the preset maximum number of iterations has been reached, the iteration is terminated.

8. The lens design method for interlayer angle cumulative offset according to claim 7, characterized in that, S4 further includes: After gradient optimization converges, a rational number avoidance test is performed. When the ratio of the cumulative inter-layer angle coefficient to twice pi is equal to a simple rational number, the cell positions in the spiral lattice are periodic. Every few layers in the angle space, a complete alignment occurs, forming equally spaced radial alignment bands, and the gaps between them constitute a visual corridor. The candidate angle coefficient is divided by twice pi to obtain a dimensionless ratio. The dimensionless ratio is compared one by one with all simple rational numbers whose denominators do not exceed the preset upper limit of the denominator. If the absolute value of the difference between the ratio of a candidate angle coefficient to twice pi and a certain simple rational number is less than or equal to a preset judgment threshold, it is determined to fall into the neighborhood of rational numbers. A preset offset is applied to the candidate angle coefficient, and the occlusion rate is recalculated after the offset. If it still falls into the neighborhood, the offset continues until it passes the test. The criterion for passing the test is that the difference with all examined simple rational numbers is greater than the judgment threshold. The final optimized angle coefficient, the corresponding occlusion rate, and the set of optimized point coordinates corresponding to the optimized angle coefficient are recorded.

9. The lens design method for interlayer angle cumulative offset according to claim 7, characterized in that, The set target occlusion rate includes: The target occlusion rate is set in different levels according to the application scenario of the lens. The first target occlusion rate threshold is set for applications with high requirements for nighttime vision, and the second target occlusion rate threshold is set for applications that are mainly used during the day. The first target occlusion rate threshold is greater than the second target occlusion rate threshold.

10. The lens design method for interlayer angle cumulative offset according to claim 7, characterized in that, The extraction of optimization angle coefficients by combining numerical gradient iterative optimization also includes: Boundary constraints are applied to the updated angle coefficients. The lower bound of the angle coefficients is determined based on the minimum effective interleaving requirement, and the upper bound of the angle coefficients is determined based on the optical area loss limit to avoid excessive overlap. The adaptive attenuation is calculated by dividing the current step size by the initial step size and the attenuation factor related to the number of iterations.

11. The lens design method for interlayer angle cumulative offset according to claim 1, characterized in that, The smooth transition constraint verification and boundary clipping process includes: For each coordinate point in the optimized dot matrix coordinate set, perform boundary mapping of the actual lens contour, check whether each coordinate point is within the boundary of the actual lens contour, and perform clipping and removal processing on coordinate points that exceed the actual contour boundary to adapt to the non-circular frame contour.

12. The lens design method for interlayer angle cumulative offset according to claim 1, characterized in that, The process of generating the final lens processing coordinate file also includes: Optical parameters, including additional diopter, surface radius of curvature, and aspheric coefficient, are added to each microstructure unit in the coordinate file to generate a complete machining file that can be directly used in the manufacturing process.

13. A lens design system for interlayer angle cumulative offset, used to perform the steps in the lens design method for interlayer angle cumulative offset as described in any one of claims 1-12, characterized in that, include: The topology parameter construction module is used to obtain the diameter of the optical region of the lens and the nominal diameter of the microstructure unit. After geometric mapping and discretization, the polar coordinate hierarchical topology parameters are constructed. The angle coefficient pre-optimization module is used to generate pre-optimized angle coefficients and initial phase sequences for each layer based on the polar coordinate layered topology parameters and the nominal diameter of the microstructure unit, through geometric occlusion condition derivation and safety factor adjustment strategy. The spiral lattice construction module is used to integrate pre-optimized angle coefficients, initial phase sequences of each layer, and polar coordinate layered topology parameters. Based on the target line density method and polar coordinate rotation transformation mechanism, it constructs a spiral lattice coordinate set and a sequence of points in each layer. The occlusion assessment and optimization module is used to call the radial ray collision detection algorithm to assess the occlusion of the spiral lattice coordinate set. It combines numerical gradient iterative optimization and rational number avoidance strategy to extract the optimized angle coefficient and optimized lattice coordinate set that maximizes the corridor occlusion rate. The coordinate file generation module is used to generate the final lens processing coordinate file by taking the optimized angle coefficients, optimized point matrix coordinate set and the number of points in each layer, verifying the smooth transition constraint and performing boundary clipping.

Images