Cutting height inversion method

Abstract

The invention relates to the field of atmospheric remote sensing detection, in particular to a tangent height inversion method, which comprises the following steps of: acquiring knee point radiation in satellite remote sensing limb scattered radiation data; generating a cut-high probability solution; defining a fitness function according to knee point radiation; performing multi-round iteration, and for each round of iteration, selecting from the cut-high possible solutions according to the fitness function to obtain a first intermediate possible solution; performing digital position exchange on the first intermediate possible solution and a randomly formed cross pair to obtain a second intermediate possible solution; randomly carrying out numerical modification on the second intermediate possible solution to obtain a third intermediate possible solution; and when a preset iteration completion condition is met, selecting a numerical value with the highest fitness value from the third intermediate possible solutions of the last round of iteration as a target cut height. According to the method, the calculation efficiency is improved, meanwhile, the requirement for calculation resources is reduced, and real-time or near-real-time tangent height inversion becomes possible.

Description

Technical Field

[0001] The invention relates to the field of atmospheric remote sensing detection, and in particular to a cutting height inversion method. Background Art

[0002] Cut-off height refers to a specific height in the atmosphere, which is a critical parameter in many scientific studies and applications. For example, weather forecasts, climate models, and environmental science research all require accurate knowledge of the cut-off height. Satellite remote sensing limb scatter radiation detection technology is an effective atmospheric detection mode to determine the cut-off height of the atmosphere, because the limb scatter radiation will have significant characteristic changes at a specific cut-off height. Therefore, inverting the cut-off height based on satellite remote sensing limb scatter radiation is of great significance to improving the accuracy of atmospheric science research.

[0003] Existing methods mainly calculate the cut height through model simulation and data matching. These methods usually require a lot of computing resources because they need to simulate a large number of possible situations and compare them with actual observation data. However, these methods have low computational efficiency and cannot meet real-time or near-real-time requirements. Summary of the invention

[0004] In order to solve the above problems, the present invention provides a cutting height inversion method.

[0005] The cutting height inversion method includes:

[0006] Obtain knee radiation from satellite remote sensing limb-scatter radiation data;

[0007] Generate M possible solutions for cutting height;

[0008] According to the knee point radiation, the fitness function f is defined;

[0009] Perform multiple rounds of iterations, and for each round of iteration:

[0010] According to the fitness function f, select from M possible solutions of cutting height to obtain M first intermediate

[0011] Possible solution;

[0012] For the M first intermediate possible solutions, randomly form crossover pairs and exchange the positions of the numbers to obtain M

[0013] The second intermediate possible solution;

[0014] Randomly modify the values ​​of the M second intermediate possible solutions to obtain M third intermediate possible solutions;

[0015] The M third intermediate possible solutions are used as the M cutting height possible solutions for the next round of iteration;

[0016] When the preset iteration completion condition is reached, the value with the highest fitness value is selected from the M third intermediate possible solutions of the last round of iteration as the target cutting height, and the fitness value refers to the value calculated according to the fitness function f.

[0017] Furthermore, the step of generating M possible solutions of cutting height includes:

[0018] Randomly generate M binary codes of the first length as the binary codes corresponding to the M possible solutions of the height cut;

[0019] Among them, the i-th cutting height may solve h i for:

[0020]

[0021] Among them, h max Indicates the maximum possible solution value of the preset cutting height, h min Indicates the minimum possible solution of the preset cutting height, L is the preset first length, x ik Represents the k-th binary code value of the possible solution of the i-th cutting height.

[0022] Furthermore, according to the knee radiation, the steps of defining the fitness function f include:

[0023] Define the fitness function f as:

[0024]

[0025] Among them, I obs (λ,h knee ) represents the limb scattered radiation in the satellite remote sensing limb scattered radiation data, with a wavelength of λ and a tangent height corresponding to the knee point, h i represents the possible solution of the i-th cutting height, I mod (λ,h i ) represents the wavelength λ and the height h i When the simulated radiation is calculated according to the radiation transfer model, h knee represents the tangent height corresponding to the knee point, and ε represents a preset positive real number.

[0026] Further, according to the fitness function f, the step of selecting from the M possible solutions of cutting height to obtain M first intermediate possible solutions includes:

[0027] Take the possible solution with the highest fitness value as a first intermediate possible solution;

[0028] Then, according to the fitness value corresponding to each possible solution of cutting height, one is selected from the M possible solutions of cutting height as a first intermediate possible solution, and a total of M-1 times are selected.

[0029] Further, according to the fitness value corresponding to each possible solution of cutting height, the step of selecting one of the M possible solutions of cutting height as a first intermediate possible solution includes:

[0030] Each time a selection is made, a random number rand is generated, and its value range is (0,1];

[0031] Select the possible solution of cutting height according to the following formula:

[0032]

[0033] Among them, h i represents the possible solution of the ith cutting height, Represents the cumulative probability of the possible solution of the i-th cutting height.

[0034] Furthermore, the cumulative probability of the possible solution of the i-th cutting height is The calculation method is:

[0035]

[0036] Among them, j represents the serial number of the possible solutions of the i-th and previous cutting heights, p j Represents the probability of the jth possible solution being selected;

[0037] Among them, the probability that the i-th possible solution is selected is p i for:

[0038]

[0039] Among them, f i Represents the fitness value of the i-th possible solution of cutting height calculated according to the fitness function.

[0040] Furthermore, the step of randomly forming crossover pairs and exchanging digital positions of the M first intermediate possible solutions to obtain M second intermediate possible solutions includes:

[0041] From the M first intermediate possible solutions, two first intermediate possible solutions are selected each time as a crossover pair, for a total of M / 2 times, and each first intermediate possible solution is selected at most once;

[0042] For each crossover pair:

[0043] Generate crossover random numbers;

[0044] Determine whether the crossover random number is less than or equal to the preset crossover threshold:

[0045] When the result is yes, the binary values ​​of the two first intermediate possible solutions in the crossover pair are exchanged in digital position;

[0046] When the answer is no, each first intermediate possible solution in the crossover pair is used as a second intermediate possible solution.

[0047] Furthermore, the step of exchanging the digital positions of the binary values ​​of the two first intermediate possible solutions in the cross pair includes:

[0048] Generate binary random numbers q, q for the crossover pair <L;

[0049] In the binary codes corresponding to the two first intermediate possible solutions in the cross pair, the qth bit and subsequent binary strings are interchanged, and each interchanged first intermediate possible solution is used as a second intermediate possible solution.

[0050] Furthermore, the step of randomly modifying the values ​​of the M second intermediate possible solutions includes:

[0051] For each bit of the binary code corresponding to each second intermediate possible solution, generate a mutated random number;

[0052] Determine whether the mutation random number is less than or equal to the preset mutation threshold:

[0053] When the result is yes, the value of the corresponding bit is inverted;

[0054] When the result is negative, the value of the corresponding bit remains unchanged.

[0055] Furthermore, the iteration completion condition refers to reaching a preset number of iterations.

[0056] The present invention provides a method for inversion of cutting height, by generating a batch of possible solutions of random cutting height, and assigning a fitness value to each possible solution according to the difference between the predicted radiation and the actually observed radiation. Then the cutting height with the highest fitness value is retained, and other possible cutting heights are selected to generate a new set of possible solutions. These new solutions generate more possible solutions through crossover and mutation operations. This process will continue until the set number of iterations is met, and then the cutting height with the highest fitness value is selected as the final result. This method not only improves the calculation efficiency, but also reduces the demand for computing resources, making real-time or near real-time cutting height inversion possible. BRIEF DESCRIPTION OF THE DRAWINGS

[0057] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0058] Figure 1It is the orbital cut height of OMPS limb observation 49270-49273 on May 1, 2021, inverted according to the method of the present invention disclosed in an embodiment of the present invention. DETAILED DESCRIPTION

[0059] In order to enable those skilled in the art to better understand the scheme of the present invention, the technical scheme in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0060] The terms "first", "second", etc. in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish different objects, rather than to describe a specific order. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions. For example, a process, method, device, or product end including a series of steps or units is not limited to the listed steps or units, but may optionally include steps or units that are not listed, or may optionally include other steps or units inherent to these processes, methods, products, or ends.

[0061] Reference to "embodiments" herein means that a particular feature, structure, or characteristic described in conjunction with the embodiments may be included in at least one embodiment of the present invention. The appearance of the phrase in various places in the specification does not necessarily refer to the same embodiment, nor is it an independent or alternative embodiment that is mutually exclusive with other embodiments. It is explicitly and implicitly understood by those skilled in the art that the embodiments described herein may be combined with other embodiments.

[0062] The present invention provides a method for inverting the cutting height, and generates a batch of random possible cutting height solutions for satellite remote sensing edge scattered radiation. Each possible cutting height solution is assigned a fitness value, which reflects the difference between the edge scattered radiation predicted by the cutting height and the actually observed edge scattered radiation. Subsequently, the cutting height with the highest fitness value is retained, and other cutting heights are selected according to the fitness value of each cutting height to form a new set of possible solutions. These new possible solutions will generate more possible cutting height solutions through mutual combination and random adjustment. This process will continue until the set number of iterations is met, and then the possible solution with the highest fitness value is selected as the cutting height corresponding to the edge scattered radiation.

[0063] The present invention is explained below in conjunction with a specific embodiment.

[0064] 1. Data Acquisition

[0065] Satellite remote sensing limb scattered radiation data is obtained by earth observation satellites. These satellites are equipped with high-resolution spectrometers that can observe the radiation intensity of different bands in the atmosphere. Limb scattered radiation refers to the part of the sunlight that is scattered in the atmosphere and radiates from the edge of the atmosphere (i.e. the earth's skyline) toward the satellite. This part of the light contains a lot of atmospheric information, such as gas concentration, aerosol, temperature, etc.

[0066] In the satellite remote sensing edge scattered radiation data, when the UV-B band (about 280-320nm) line passes through the stratosphere, due to the absorption of ozone, the intensity of its edge scattered radiation will show a peak at a certain cut height. This peak is called the knee point, and the cut height corresponding to the knee point is the target cut height that needs to be inverted in the present invention.

[0067] 2. Possible deinitialization

[0068] For the target cutting height to be solved, M possible cutting height solutions are generated, each of which represents a possible true value of the target cutting height. In this embodiment, M is 20.

[0069] A binary code with a first length is randomly generated as the binary code corresponding to the possible solution of the height cut. In this embodiment, the first length is 50.

[0070] The i-th cutting height may solve h i for:

[0071]

[0072] Among them, h max Indicates the maximum possible solution value of the preset cutting height, h min Indicates the minimum possible solution of the preset cutting height, L is the preset first length, x ik Represents the k-th binary code value of the possible solution of the i-th cutting height.

[0073] That is, the purpose of the present invention is to solve and obtain the target cutting height to satisfy the following formula:

[0074] F obj =max(I mod (λ,h i ));

[0075] Where λ represents the wavelength, I mod (λ,h i ) indicates that the wavelength λ is cut off at h iWhen the edge scattered radiation is calculated according to the radiation transfer model, the value of the wavelength λ must satisfy the corresponding knee point characteristic at high altitudes, and there are fewer or even no aerosols in high altitudes. In this embodiment, λ is 305nm. The radiation transfer models include MODTRAN (MODerateresolution atmospheric TRANsmission, medium-resolution atmospheric transfer model), SCIATRAN (radiative TRANsfer model for SCIAMACHY, SCIAMACHY radiation transfer model), etc. How to calculate the edge scattered radiation corresponding to the cut height according to the radiation transfer model is common knowledge in the field and will not be repeated in the present invention.

[0076] 3. Fitness function definition

[0077] The fitness function is a measure of the quality of possible solutions to the cutting height, and guides the iteration process. The present invention uses the absolute difference between the knee point radiation actually observed by the satellite and the simulated radiation to design the fitness function.

[0078] Knee point radiation is the strongest edge scattered radiation in satellite remote sensing edge scattered radiation data. Simulating radiation is to use the radiation transfer model. According to the i-th cut height, h i , the calculated limb scattered radiation.

[0079] The fitness function f is:

[0080]

[0081] Among them, I obs (λ,h knee ) represents the limb scattered radiation in the satellite remote sensing limb scattered radiation data, with a wavelength of λ and a tangent height corresponding to the knee point, I mod (λ,h i ) represents the wavelength λ and the height h i When the simulated radiation is calculated according to the radiation transfer model, h knee represents the cut height corresponding to the knee point, and ε represents a preset very small positive real number to prevent the denominator of the fitness function from being 0. In this embodiment, ε is 0.00001.

[0082] 4. Iterative update cutting height possible solution

[0083] 4.1 Select the first possible intermediate solution

[0084] In each round of iteration, the cutting height possible solution with the highest fitness value is used as a first intermediate possible solution, and then according to the fitness value corresponding to each cutting height possible solution, one is selected from the M cutting height possible solutions as a first intermediate possible solution, and a total of M-1 times are selected. In this step, a total of M first intermediate possible solutions are obtained. Since the selected cutting height possible solution will not be removed from the M cutting height possible solutions after each selection, the M first intermediate possible solutions may have the same value.

[0085] Each time a choice is made, the probability p of each possible solution being chosen is i for:

[0086]

[0087] The cumulative probability of each possible solution for:

[0088]

[0089] Among them, f i represents the fitness value of the i-th possible solution for cutting height calculated according to the fitness function, and j represents the sequence number of the i-th and previous possible solutions for cutting height.

[0090] Each time the first possible intermediate solution is selected from the M possible solutions for cutting height, a random number rand is generated, whose value range is (0,1], and the possible solution for cutting height is selected according to the following formula:

[0091]

[0092] For example, the four possible solutions are the first possible solution, the second possible solution, the third possible solution, and the fourth possible solution. The corresponding fitness values ​​are 2, 4, 6, and 8, respectively, and the corresponding probabilities are p i They are 0.1, 0.2, 0.3, and 0.4 respectively, and the corresponding cumulative probabilities They are 0.1, 0.3, 0.6, and 1 respectively. From the four possible solutions for cutting height, the fourth possible solution for cutting height with the highest fitness value is selected first, and then three random numbers are generated, which are 0.57, 0.78, and 0.04 respectively. The random number 0.3<0.57≤0.6, so the random number 0.57 corresponds to the third possible solution for cutting height. Similarly, the random number 0.78 corresponds to the fourth possible solution for cutting height, and the random number 0.04 corresponds to the first possible solution for cutting height. That is, the four first intermediate possible solutions are: the fourth possible solution for cutting height, the third possible solution for cutting height, the fourth possible solution for cutting height, and the first possible solution for cutting height.

[0093] In each iteration, select the one with the highest fitness value from the current set of possible solutions for the cutting height. Then, based on the fitness values of each possible solution for the cutting height, make multiple random selections. In this way, possible solutions for the cutting height with high fitness values have a higher probability of being selected multiple times. Although possible solutions for the cutting height with low fitness values have a lower probability of being selected, they still have a chance to be selected, thus maintaining the diversity of the solution space.

[0094] 4.2 Generate the second intermediate possible solutions

[0095] From the M first intermediate possible solutions, each time select two first intermediate possible solutions as a crossover pair, and select a total of M / 2 times. Each first intermediate possible solution can be selected at most once.

[0096] Generate a crossover random number for each crossover pair, and its value range is [0, 1]. When the crossover random number is less than or equal to the preset crossover threshold, generate a binary bit random number q for this crossover pair, and swap the binary strings after the q-th bit in the binary encodings corresponding to the two first intermediate possible solutions in this crossover pair. Take each swapped first intermediate possible solution as a second intermediate possible solution, where L / 2 ≤ q < L. When the crossover random number is greater than the crossover threshold, take each first intermediate possible solution in this crossover pair as a second intermediate possible solution. In this embodiment, the crossover threshold is set to 0.8.

[0097] For example: The binary encodings corresponding to the two first intermediate possible solutions in a crossover pair are 11110000 and 00001111 respectively, L is 8, q is 5, and the crossover random number is less than the preset crossover threshold, then the two second intermediate possible solutions are 11111111 and 00000000 respectively.

[0098] Through the crossover operation, information exchange can be carried out between possible solutions, exploring more new possible solutions, and improving the diversity of the search and the global search ability.

[0099] 4.3 Generate the third intermediate possible solutions

[0100] For each bit of the binary encoding corresponding to each second intermediate possible solution, generate a mutation random number, and its value range is [0, 1]. When the mutation random number is less than or equal to the preset mutation threshold, invert the value of the corresponding bit. After processing each bit, obtain the third intermediate possible solution. In this embodiment, the mutation threshold is set to 0.05.

[0101] For example: the binary code corresponding to the second intermediate possible solution is 10101010, and the mutant random numbers generated for the eight binary bits are: 0.12, 0.02, 0.87, 0.01, 0.55, 0.09, 0.31, 0.04. Only the second bit (random number 0.02) and the eighth bit (random number 0.04) are less than or equal to the mutant random number 0.05, so only the second and eighth bits are inverted. The binary code corresponding to the third intermediate possible solution is: 11101011.

[0102] The mutation operation can introduce randomness into possible solutions, prevent the search process from falling into the local optimal solution, and increase the global search capability. At the same time, by introducing new features, it is also possible to find a better solution.

[0103] After the mutation operation, the M third intermediate possible solutions are used as M possible solutions for cutting height, and the next round of iteration is performed until the preset number of iterations is reached.

[0104] Finally, the value with the highest fitness value among the third intermediate possible solutions of the last iteration is taken as the target cutting height.

[0105] 5. Effect comparison

[0106] The present invention uses the limb observation mode of the Ozone Mapping and Profiler (OMPS) as the experimental object. Figure 1 This is the orbital cut height of OMPS limb observation 49270-49273 on May 1, 2021, inverted by the method of the present invention. Figure 1 It can be seen that at the beginning of each orbit and in the latitude zone south of 30° in the southern hemisphere, the method of the present invention can still be inverted to obtain the inverted height. The beginning of the orbit generally refers to SZA>90°, SZA (Solar Zenith Angle).

[0107] Table 1 shows the average value of the OMPS limb observation 49270-49273 orbit cut height inverted by the method of the present invention and the comparison result with the cut height in the OMPS Level 1G data (L1G, Level 1G) provided by NASA (National Aeronautics and Space Administration). The inverted cut height in Table 1 represents the inverted cut height of the OMPS limb observation 49270-49273 orbit obtained by a cut height inversion method provided by the present invention, the NASA cut height in Table 1 represents the cut height of the 49270-49273 orbit in the Level 1G data, and ACC represents the accuracy of the cut height correction, which is calculated as follows:

[0108] Acc=|TH ret -TH NASA |;

[0109] Among them, TH ret is the cutting height obtained by the method provided by the present invention, TH NASA The cut height of Level 1G (L1G) data provided by NASA.

[0110] Table 1 Orbital inversion cut height, NASA cut height and accuracy of the method of the present invention

[0111] track Inversion cutting height (km) NASA Cut High Acc(km) 49270 47.56 47.36 0.20 49271 47.57 47.41 0.16 49272 47.54 47.39 0.15 49273 47.55 47.38 0.17 average 47.56 47.39 0.17

[0112] It can be seen from Table 1 that the average accuracy of all tracks is 0.17km, which is less than 250m, meeting the requirements of height correction.

[0113] In summary, the present invention discloses a shear height inversion method, which is suitable for various latitudes. It adopts the inverse of the absolute value of the logarithmic difference between the wavelength maximum observed radiation and the simulated radiation calculated by the radiation transfer model to invert the shear height as the fitness function to evaluate the advantages and disadvantages of possible shear height solutions. It has high inversion accuracy and only uses the wavelength of the UV-B band to reduce the interference of aerosols. The limb radiation of different wavelengths has different knee heights, and only uses the limb radiation characteristics of a single wavelength (305nm in the embodiment), which improves the shear height inversion accuracy and speed. Therefore, the present invention effectively overcomes the various shortcomings in the prior art and has a high industrial utilization value.

[0114] The above embodiments are merely illustrative of the principles and effects of the present invention, and are not intended to limit the present invention. Anyone familiar with the art may modify or alter the above embodiments without departing from the spirit and scope of the present invention. Therefore, all equivalent modifications or alterations made by a person of ordinary skill in the art without departing from the spirit and technical concept disclosed by the present invention shall still be covered by the claims of the present invention.

Claims

1. A cutting height inversion method, characterized in that: The method comprises: Obtain knee radiation from satellite remote sensing limb-scatter radiation data; Generate M possible solutions for cutting height; According to the knee point radiation, the fitness function f is defined; Perform multiple rounds of iterations, and for each round of iteration: According to the fitness function f, select from M possible solutions of cutting height to obtain M first intermediate possible solutions; For the M first intermediate possible solutions, randomly form crossover pairs and exchange the positions of the numbers to obtain M second intermediate possible solutions; Randomly modify the values ​​of the M second intermediate possible solutions to obtain M third intermediate possible solutions; The M third intermediate possible solutions are used as the M cutting height possible solutions for the next round of iteration; When the preset iteration completion condition is reached, the value with the highest fitness value is selected from the M third intermediate possible solutions of the last round of iteration as the target cutting height, and the fitness value refers to the value calculated according to the fitness function f.

2. A cutting height inversion method according to claim 1, characterized in that: The steps of generating M possible solutions for cutting height include: Randomly generate M binary codes of the first length as the binary codes corresponding to the M possible solutions of the height cut; Among them, the i-th cutting height may solve h i for: Among them, h max Indicates the maximum possible solution value of the preset cutting height, h min Indicates the minimum possible solution of the preset cutting height, L is the preset first length, x ik Represents the k-th binary code value of the possible solution of the i-th cutting height.

3. A cutting height inversion method according to claim 1, characterized in that: According to the knee radiation, the steps to define the fitness function f include: Define the fitness function f as: Among them, I obs (λ,h knee ) represents the limb scattered radiation in the satellite remote sensing limb scattered radiation data, with a wavelength of λ and a tangent height corresponding to the knee point, h i represents the possible solution of the i-th cutting height, I mod (λ,h i ) represents the wavelength λ and the height h i When the simulated radiation is calculated according to the radiation transfer model, h knee represents the tangent height corresponding to the knee point, and ε represents a preset positive real number.

4. A cutting height inversion method according to claim 2, characterized in that: The steps of selecting from M possible solutions for cutting height according to the fitness function f to obtain M first intermediate possible solutions include: Take the possible solution with the highest fitness value as a first intermediate possible solution; Then, according to the fitness value corresponding to each possible solution of cutting height, one is selected from the M possible solutions of cutting height as a first intermediate possible solution, and a total of M-1 times are selected.

5. A cutting height inversion method according to claim 4, characterized in that: The step of selecting one of the M possible solutions for cutting height as a first intermediate possible solution according to the fitness value corresponding to each possible solution for cutting height includes: Each time a selection is made, a random number rand is generated, and its value range is (0,1]; Select the possible solution of cutting height according to the following formula: Among them, h i represents the possible solution of the ith cutting height, Represents the cumulative probability of the possible solution of the i-th cutting height.

6. A cutting height inversion method according to claim 5, characterized in that: The cumulative probability of the possible solution of the i-th cutting height The calculation method is: Among them, j represents the serial number of the possible solutions of the i-th and previous cutting heights, p j Represents the probability of the jth possible solution being selected; Among them, the probability that the i-th possible solution is selected is p i for: Among them, f i Represents the fitness value of the i-th possible solution of cutting height calculated according to the fitness function.

7. A cutting height inversion method according to claim 2, characterized in that: The steps of randomly forming crossover pairs and exchanging digital positions of M first intermediate possible solutions to obtain M second intermediate possible solutions include: From the M first intermediate possible solutions, two first intermediate possible solutions are selected each time as a crossover pair, for a total of M / 2 times, and each first intermediate possible solution is selected at most once; For each crossover pair: Generate crossover random numbers; Determine whether the crossover random number is less than or equal to the preset crossover threshold: When the result is yes, the binary values ​​of the two first intermediate possible solutions in the crossover pair are exchanged in digital position; When the answer is no, each first intermediate possible solution in the crossover pair is used as a second intermediate possible solution.

8. A cutting height inversion method according to claim 7, characterized in that: The step of exchanging the digital positions of the binary values ​​of the two first intermediate possible solutions in the cross pair comprises: Generate binary random numbers q, q for the crossover pair <L; In the binary codes corresponding to the two first intermediate possible solutions in the cross pair, the qth bit and subsequent binary strings are interchanged, and each interchanged first intermediate possible solution is used as a second intermediate possible solution.

9. A cutting height inversion method according to claim 2, characterized in that: The steps of randomly modifying the values ​​of the M second intermediate possible solutions include: For each bit of the binary code corresponding to each second intermediate possible solution, generate a mutated random number; Determine whether the mutation random number is less than or equal to the preset mutation threshold: When the result is yes, the value of the corresponding bit is inverted; When the result is negative, the value of the corresponding bit remains unchanged.

10. A cutting height inversion method according to claim 1, characterized in that: The iteration completion condition refers to reaching a preset number of iterations.

Images