Next Article in Journal
Analysis of the Photogrammetric Use of 360-Degree Cameras in Complex Heritage-Related Scenes: Case of the Necropolis of Qubbet el-Hawa (Aswan Egypt)
Previous Article in Journal
Synthetic Iris Images: A Comparative Analysis between Cartesian and Polar Representation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improved Immune Moth–Flame Optimization Based on Gene Correction for Automatic Reverse Parking

1
College of Engineering, Inner Mongolia Minzu University, Tongliao 028000, China
2
Department of Automation, Shanghai Jiao Tong University, Shanghai 200240, China
3
School of Mechanical Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China
4
School of Automation and Electrical Engineering, Dalian Jiaotong University, Dalian 116026, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(7), 2270; https://doi.org/10.3390/s24072270
Submission received: 18 January 2024 / Revised: 24 March 2024 / Accepted: 28 March 2024 / Published: 2 April 2024
(This article belongs to the Section Intelligent Sensors)

Abstract

:
During the process of reverse parking, it is difficult to achieve the ideal reference trajectory while avoiding collision. In this study, with the aim of establishing reference trajectory optimization for automatic reverse parking that smooths and shortens the trajectory length and ensures the berthing inclination angle is small enough, an improved immune moth–flame optimization method based on gene correction is proposed. Specifically, based on the standard automatic parking plane system, a reasonable high-quality reference trajectory optimization model for automatic parking is constructed by combining the cubic spline-fitting method and a boundary-crossing solution based on gene correction integrated into moth–flame optimization. To enhance the model’s global optimization performance, nonlinear decline strategies, including crossover and variation probability and weight coefficient, and a high-quality solution-set maintenance mechanism based on fusion distance are also designed. Taking garage No.160 of the Dalian Shell Museum located in Dalian, Xinghai Square, as the experimental site, experiments on automatic parking reference trajectory optimization and tracking control were carried out. The results show that the proposed optimization algorithm provides higher accuracy for reference trajectory optimization and can achieve better tracking control of the reference trajectory.

1. Introduction

Automatic parking behavior whereby vehicles meet the reverse parking conditions from a starting area within an allowable range is called automatic reverse parking. With the rapid development of transportation, the demand for comprehensive performance indicators for intelligent vehicles is increasing, and the performance quality of intelligent ARP functions is also increasingly scrutinized [1]. The reverse automatic parking system is indispensable for intelligent vehicles to ensure optimal reverse parking performance indicators, such as safety, comfort, speed, and accuracy. Its reference trajectory optimizer should exhibit reverse parking trajectory optimization performance to provide collision-free, accurate automatic parking trajectories, and the path length should be smooth and short and the berthing inclination angle small [2,3]. However, given the optimization abilities of existing optimization algorithms, the ideal conditions are difficult to achieve [4].
Regarding the problems encountered when applying automatic parking systems to intelligent vehicles, scholars have carried out a series of studies. As early as 1996, a novel neural network architecture for trajectory generation and the control of automated car parking was proposed [5]. In [6], the authors undertook the design, testing, and validation of a deep neural network (DNN)-based control scheme for autonomous ground vehicles (AGVs) during the parking process. An automatic vacant parking space management system using multicamera vehicle detection was presented in [7]. A graph-based algorithm for the optimal control of switched systems for application to car parking was developed in [8]. A distributed control method for coordinating multiple vehicles within the framework of an automated valet parking system was introduced in [9]. However, there are relatively few effective improvement strategies for the reference trajectory optimization of intelligent ARP.
The parking trajectory optimizer has a strong influence on the performance quality of automatic parking systems. In recent years, intelligent optimization algorithms have received significant attention from researchers in the field of automatic parking, with many scholars around the world having conducted studies on this topic. A computational trajectory optimization framework for multiobjective automatic parking motion planning was proposed in [10]. An improved A* and dynamic window approach was proposed in [11] to enhance the accuracy and speed of automatic parking in tracked vehicles. Improved ant colony optimization (ACO) for path planning in an AGV-based intelligent parking system was introduced in [12]. Improved immune moth–flame optimization (IIMFO) based on the immune mechanism, Gaussian mutation mechanism, and opposition-based learning strategy was proposed for automatic parking path optimization in intelligent vehicles in [13]. However, the high-performance optimization frameworks using improved immune moth–flame optimization based on gene correction (IIMFO-GC) for ARP trajectory optimization lack guidelines for their design, so it may be difficult to obtain optimal ARP reference trajectories.
With the rapid development of the public transportation system, the quality requirements for daily travel are also increasing. However, for intelligent ARP, the reference trajectory optimization performance of traditional parking trajectory optimizers is often unsatisfactory. Based on a study [13] about reference trajectory optimization for automatic reverse parking in intelligent vehicles, the major contributions of this work are as follows:
(I) A reasonable high-quality reference trajectory optimization model for ARP is constructed by combining the cubic spline-fitting method and a boundary-crossing solution based on gene correction integrated into moth–flame optimization.
(II) An improved immune moth–flame optimization method is designed. Specifically, nonlinear decline strategies for weight coefficient and crossover and mutation probabilities, and a high-quality solution-set maintenance mechanism based on fusion distance, are designed to enhance the global optimization performance of immune moth–flame optimization.
(III) A gene correction method is designed to establish a highly effective boundary-crossing solution for ARP; on this basis, an IIMFO-GC is proposed and used in the trajectory optimizer of the ARP system to further improve reference trajectory optimization performance.
This paper is organized as follows: Section 2 introduces the reference trajectory optimization model for ARP; Section 3 introduces the gene correction method for reverse parking; Section 4 introduces the IIMFO; Section 5 presents the ARP experimental results and simulation analysis; and Section 6 concludes the article.

2. Reference Trajectory Optimization Model for ARP

2.1. Basic Principles of ARP

ARP generally adopts the method of automatic reversing. A schematic diagram of the ARP principle is shown in Figure 1.
In Figure 1, it is assumed that the vehicle coverage areas are purposefully collected in sequence n s times in the ARP process, which is denoted as the collected vehicle coverage area set ( A 1 , , A i s , , A n s ) . In this figure, A 1 represents the parking starting area, A n s is the parking stop area, and P 1 , P i s , and P n s are the center points of the collected vehicle coverage areas. L d represents the bottom line of parking garage, and L n e and L f e are the near and far side lines of the parking garage, respectively. P n t , P f t , P n b and P f b are the near and far corner points of the garage top and the near and far corner points of garage bottom, respectively. T R and T C are the reference parking trajectory and the actual tracking control trajectory, respectively. Due to the uncertain interference in the parking process, it is impossible to achieve a tracking control error of zero, that is, T C cannot completely coincide with T R . T S is the tracking control trajectory monitoring value obtained by the sensor and image system. Due to the limitation in sensor accuracy and measurement error, it is impossible to achieve an error value of zero for tracking control trajectory monitoring, that is, T S cannot completely coincide with T C .

2.2. Constraint Judgment of ARP Feasibility

ARP usually requires that there are no obstacles in the parking environment, that the parking garage is known, and that the bottom and side lines of the parking garage, as well as the bottom and top corners of the parking garage, are clearly identifiable. In addition, the constraint on the starting feasibility of parking in the garage should be satisfied: the center point P 1 of the starting area should be within the feasible starting area Ω C , 1 . The angle between the line connecting the center point P 1 of the starting area with its bottom midpoint P b 1 and the bottom line L d of the garage is called the starting inclination angle Δ B . If its absolute value is less than the threshold Δ B max , the constraint expression of the starting feasibility of reverse parking is as follows :
P 1 Ω C , 1 B = T P 1 , P b 1 , L d < Δ B max

2.3. Conversion Principle of ARP Plane Coordinate System

The standard ARP plane coordinate system takes the far corner point of the bottom of the garage as the origin and the bottom line of the garage as the positive horizontal axis (x-axis). Adopting a standard ARP plane coordinate system is not only beneficial for simplifying reverse parking calculations, but also enables drivers to evaluate real-time conditions and make more rational parking decisions. To convert the real scene’s ARP plane coordinate system obtained by the parking image system into a standard ARP plane system, we must carry out three steps: first, determine the largest range of the real rectangular parking scene Ω ; then, rotate the coordinate system according to the real scene’s rotation angle α ; finally, translate the origin according to the coordinate x 0 , y 0 of the far corner point of the garage bottom in the real scene’s coordinate system.
Assume there is a point s in the parking area whose coordinate in the standard ARP plane system is x s , y s and whose angle with the horizontal axis is β . Then, β = arcsin y s x s 2 + y s 2 .
If the coordinate x s , y s is known, the conversion formula for the coordinate x r , y r of the real scene’s reverse parking plane system can be obtained as follows:
x r = x s 2 + y s 2 · cos β α + x 0 y r = x s 2 + y s 2 · sin β α + y 0
Similarly, if the coordinate x r , y r is known, it is easy to establish the coordinate x s , y s .
In the standard parking coordinate system, the set of decision variables is set as x P 2 , y P 2 , , x P n s , y P n s T , and x P i s and y P i s are the abscissa and ordinates of the central point P i s of the i s collected vehicle coverage areas in the standard reverse parking coordinate system. On the premise that the starting inclination angle Δ B is small enough ( B Δ B max ), to obtain the reverse parking reference trajectory with the shortest length, it is known from the parking experience that the abscissa and ordinates about the center points of the collected vehicle coverage areas should decrease one by one. The expression is as follows:
x P n s , , x P i s , , x P 1
y P n s , , y P i s , , y P 1

2.4. Optimization Model for ARP Reference Trajectory

The number of center points of the collected vehicle coverage areas is set to n s , and the ARP reference trajectory can be obtained by fitting n s center points P 1 , P 2 , , P n s . Among them, the center point P 1 of the starting area is fixed, and the center points P 2 , , P n s of the other n 1 collected vehicle coverage areas can be adjusted within a certain region. The periodic time of the ARP process is set to Δ t , and the optimization model of the ARP reference trajectory is expressed as follows:
min L T R = 0 T max v ( t ) d t i t = 1 n t Δ L i t Δ S = T P n s , P b m , L f e S . T . T R = t f P 1 , P 2 , , P n s Δ L i t Δ L max Δ i t Δ max P i s Ω C , i s P g ( i g ) Ω c , i t P g ( i g ) Ω g Δ S = Δ S max Δ y n s = y P n s y H Δ y H max
where Ω g is the coverage area of the edges on both sides of the garage, P g is the collection of detected edge collision avoidance point positions, n g is the detected point position, P g ( i g ) is the i g -th detected edge collision avoidance point position, i 1 , 2 , , n g (then, P g ( i g ) Ω g ), and Δ L i t , Δ i t , and Ω c , i t are the parking interval distance, attitude angle difference, and vehicle coverage area for the i t -th period, respectively. In this study, the vehicle coverage area Ω c , i t is equalized as a rectangle with the length and width of the vehicle, Δ L max and Δ max are the threshold values for parking interval distance and attitude angle difference, T max is the parking time, n t is the number of periods, n t = T max Δ t , i t 1 , 2 , , n t , is the rounding up operator, since the vehicle cannot collide with the bottom line of the garage and the edges on both sides, P g ( i ) Ω c , i t , P i s is the center point of the i s -th collected vehicle coverage area, i s 2 , , n s , Ω C , i s represents its adjustable region (then, P i s Ω C , i s ), and t f is the reverse parking reference trajectory obtained through fitting, with the cubic spline-fitting method adopted in this study; Δ S is the parking inclination angle, Δ S = T P n s , P b m , L f e and is the angle between the line T P n s , P b m connecting the center point P n s of the parking area with the midpoint P b m at the bottom of the garage and the far side line L f e of the garage, and Δ S max is the threshold for its absolute value; Δ y n s is parking position error and is the error between the ordinate y P n s of the center point P n s of the parking area and the expected ordinate y H , Δ y H max is the threshold, and T R is the reference trajectory.
The cubic spline-fitting method is a spline curve-fitting method using the cubic spline interpolation function, and it has been widely used in various industrial applications. Assuming that the closed interval Δ = a , b and n + 1 interpolation points ( x 0 , x 1 , , x n ) within the interval Δ are known, in terms of the cubic spline interpolation function s ( x ) , there are two kinds of features.
Firstly, regarding the arbitrary subinterval x i , x i + 1 , i [ 0 , 1 , · · · · · · , n 1 ] , if s ( x ) is a cubic polynomial, the cubic polynomial formula is as follows:
s ( x ) = s 0 x x [ x 0 , x 1 ] s 1 x x [ x 1 , x 2 ] · · · · · · s n 1 x x [ x n 1 , x n ]
s i x = a i + 1 · x 3 + b i + 1 · x 2 + c i + 1 · x + d i + 1
where s i x is the segmented cubic polynomial function for the cubic spline interpolation function s ( x ) , and ( a 0 , a 1 , , a n ) , ( b 0 , b 1 , , b n ) , ( c 0 , c 1 , , c n ) , and ( d 0 , d 1 , , d n ) are polynomial coefficients for the segmented cubic polynomial function s 0 x , , s n x .
Secondly, s ( x ) is a second-order continuous and differentiable within interval Δ .

3. Gene Correction Method for ARP

3.1. Basic Concepts and Application Principles of Gene Correction

Gene correction is a method of repairing problematic genes in a diseased individual so that they can be correctly expressed [14]. The appropriate use of gene correction is beneficial for enhancing the preservation effect of optimization, but if used excessively, it will lead to too many adjustable parameters, which will increase the amount of calculation required [15]. Considering computational efficiency, this paper applies genetic modifications to three situations generated in the optimization of the reference trajectory for parking: collision avoidance on the far side line of the parking garage, berthing inclination, and parking dislocation. The other situations are replaced with new boundary-crossing solutions for reverse parking.

3.2. Gene Correction for Collision on the Far Side Line of the Parking Garage

In the standard parking coordinate system, the far side line of the garage is equal to the vertical axis (y-axis). To solve the problem of collision on the far side line of parking garage, it is necessary to rotate the coordinate system counterclockwise based on the center points of the collected vehicle coverage areas within a certain range. Whether or not it is necessary to rotate the center point of the i s -th collected vehicle coverage area counterclockwise, the empirical determination formula is as follows:
y P i s y t P g e t x l b , I T , y l b , I T Δ Y f c max
where t P g e t x l b , I T , y l b , I T is the function for the bottom left corner x l b , I T , y l b , I T , and Δ Y f c max is the maximum ordinate offset.
The empirical estimation formula for correcting the abscissa movement Δ x i s , f c for collision avoidance on the far side line of parking garage is as follows:
Δ x i s , f c = x t P g e t x l b , I T , y l b , I T x P i s · c f c · x l b , I T
where c f c is the empirical coefficient of collision avoidance on the far side line of the parking garage.
A schematic diagram of gene correction for collision avoidance on the far side line of the parking garage is shown in Figure 2.
In Figure 2, P i s 2 , = x f c , y f c , where P i s 2 is the center point of the vehicle about the crossing boundary point for collision on the far side line of garage during the parking process.

3.3. Berthing Inclination Gene Correction

If the condition stating that the absolute value of the berthing inclination angle Δ S must be less than the threshold value Δ S max is not met, it is called berthing inclination. If berthing inclination is to be corrected, the abscissa of some of the center points of the collected vehicle coverage areas should be moved slightly to the left. The empirical judgment formula for this left shift is as follows:
y P i s y P n s Δ Y S I max
where Δ Y S I max is the maximum ordinate offset for correcting the berthing inclination.
The empirical estimation formula for obtaining the left-shift value Δ x i s , S I is as follows:
Δ x i s , S I = y P i s y P n s · tan ( Δ S r · Δ S max )
A schematic diagram of gene correction for berthing inclination is shown in Figure 3.
In Figure 3, the abscissas of center points P n s 1 and P n s 2 of the collected vehicle coverage areas within the interval y n s , y n s + Δ Y S I max are shifted to the left by Δ x n s 1 , S I and Δ x n s 2 , S I , achieving the correction of berthing inclination.

3.4. Berthing Dislocation Gene Correction

If the condition stating that the absolute value of the berthing stop position error Δ y n s must be less than the threshold value is not met, it is called berthing dislocation. In order to correct berthing dislocation, the ordinate of the center point of the parking area should be appropriately moved. The empirical estimation formula for the ordinate movement Δ y S P of gene correction for berthing dislocation is as follows:
Δ y S P = Δ y n s r · Δ y H max
A schematic diagram of gene correction for berthing dislocation is shown in Figure 4.
In Figure 4, since the center point P n s of the parking area is not in the allowable parking range y H Δ y H max , y H + Δ y H max , it is above the parking area and is moved down to Δ y S P to achieve the goal of correcting berthing dislocation.
Based on the gene correction method proposed in this study, only three empirical parameters need to be added: the maximum ordinate offset Δ Y f c max and the empirical coefficient c f c of abscissa movement for correcting collision on the far side line of the parking garage, and the maximum ordinate offset Δ Y S I max for correcting berthing inclination.

4. Improved Immune Moth–Flame Optimization

4.1. Moth–Flame Optimization

Moth–flame optimization (MFO), proposed in 2015, imitates the behavior of moths searching for a light source through lateral positioning [16]. The optimization process of MFO can be abstractly represented by a triplet [17]. The mapping expressions are as follows:
M F O = I , P , T
I : ϕ M , O M , F , O F
P : M M
T : M i = ω · D i j · e τ t · cos ( 2 π t ) + F j D i j = | M i F j |
where I represents initializing the moth population and flame set and calculating the fitness value of each moth and flame, M is the moth population, O M is its fitness value vector matrix, F is the flame set, O F is its fitness value vector matrix, M = M 11 M 12 M 1 d M 21 M 22 M 2 d M n 1 M n 2 M n d , O M = O M 1 O M 2 O M n , F = F 11 F 12 F 1 d F 21 F 22 F 2 d F n f 1 F n f 2 F n f d , O F = O F 1 O F 2 O F n f , n is the population size, d is the solution dimension, and n f is the number of flames (at the beginning, the number of flames is equal to the population size), n f = n ; P represents updating the running trajectory of the moth population (after updating the moth population, the fitness value of each individual moth should be recalculated; if the fitness value is better than the current flame, the flame set will be updated), T represents updating position of moths according to the logarithmic spiral law, D i j · e τ t · cos ( 2 π t ) is the logarithmic spiral term, M i is the position of the i-th moth, F j is the position of the j-th flame, D i j is the Euclidean distance between the i-th moth and the j-th flame, τ is the logarithmic spiral shape adjustment coefficient, t is a random number in the interval 1 , 1 , and ω is the weight coefficient.
The number of flames will be reduced gradually during the iteration [18]. The formula for updating the number of flames is as follows:
n f = r o u n d ( n t n 1 t max )
where t is the current iteration number, t max is the maximum number of iterations, and r o u n d ( x ) represents x rounded to the nearest integer.

4.2. Immune Mechanism

In this study, the concentration selection mechanism in the artificial immune algorithm for adjusting the distribution of moths in space is introduced into the MFO to improve its global optimization performance [19]. In the concentration selection mechanism, the formulas for calculating antibody concentration and concentration probability are as follows:
D x i = 1 j = 1 m + N f x i f x j
P d x i = 1 D x i i = 1 m + N 1 D x i = j = 1 m + N f x i f x j i = 1 m + N j = 1 m + N f x i f x j
where f x i , D x i , and P d x i are the fitness function value, antibody concentration, and concentration probability of the i-th moth, respectively, and i = 1 , 2 , , m + N .

4.3. Nonlinear Decline Strategy of Weight Coefficient

If the weight coefficient ω is large, the swarm intelligence optimization algorithm offers strong global exploration ability and is suitable for the early stage of the optimization process; contrarily, if the weight coefficient is small, the algorithm’s local development ability is strong and it is suitable for the later stage [20]. Several nonlinear decline strategies for weight coefficients ω with strong plasticity are conducive improving the global optimization ability of MFO [21]. The two weight coefficient calculation formulas ω are as follows:
ω = ω min + ω d · cos ( 1 t r ) · π 2 β
where ω d is the decreasing weight coefficient; ω d = ω m ax ω min , ω min , and ω max are the minimum and maximum weight coefficient values, respectively; t r is the iteration progress; t r = t · t max 1 ; and β is the optimization factor.
ω = ω max ω d l N int max α
where α is the decline rate and exhibits a significant difference throughout the whole iteration process.
According to the above formulas, the nonlinear declining trend for the weight coefficient ω can be optimized and adjusted. The two calculation formulas above have excellent effects on improving global optimization ability; in different scenarios the improvement effects are discrepant, and to choose the best option, extensive tests should be carried out.

4.4. Nonlinear Decline Strategy for Crossover and Mutation Probability

In the early stage of iterative optimization, population diversity should be emphasized, the intensity of gene change should be enhanced, and global exploration should be strengthened to improve individual adaptability. Following iteration optimization, preserving the existing optimization results should be emphasized, the intensity of gene changes should be reduced, and local development should be strengthened to avoid the destruction of existing advantageous genes. Therefore, mutation and crossover probability should decrease gradually with an increasing number of iterations. This paper presents a nonlinear decline strategy for directional variation and crossover probabilities p c and p m with strong smoothness based on exponential formal decline. The formulas for calculating mutation and crossover probabilities p c and p m are as follows:
p c = p c , min + p c , max · e t r β c · ln p c , min p c , max
p m = p m , min + p m , max · e t r β m · ln p m , min p m , max
where p c , min , p c , max , p m , min , and p m , max are the minimum and maximum crossover probability and the minimum and maximum mutation probability, respectively, and β m and β c are the the optimization factors for the nonlinear decline of mutation and crossover probabilities p m and p c , respectively.
According to the above formulas, the exponential form of the nonlinear decline strategy is adopted, and the nonlinear decline trend can be optimized and adjusted by selecting the most suitable optimization factors, β m and β c . This strategy has strong smoothness, and its nonlinear decline is particularly smooth. This strategy can adjust the decline rates of mutation and cross probability in real time during iteration, which is conducive to improving the balance between global exploration and local exploitation, and thus, improving the global optimization ability of the algorithm.

4.5. High-Quality Solution-Set Maintenance Mechanism Based on Fusion Distance

In the evolution process of each iteration, a high-quality solution set should be established to preserve the existing optimization achievements. To prevent adverse effects on the algorithm’s computational efficiency, the high-quality solution-set size H S s should be less than the allowable maximum H S A S , that is, H S S H S A S . So, the design of an elite set maintenance mechanism is very significant for MFO. In terms of the definition of the Mahalanobis distance M d , the M d results between two samples under identical states in different reference spaces will differ because of the difference in their correlation of variables. Since the Euclidean distance E d considers the independence of characteristic variables, and M d considers the independence of characteristic variables, this study combines M d and E d to effectively measure distance in the case of fuzzy correlation, and a new calculation method called “Fusion distance F d ” that considers the correlation and independence between feature variables is given [22]. The “Fusion distance F d ” equation is as follows:
F d = μ · M d Y b , X a + 1 μ · E d Y b , X a
μ = 1 C M b a
C M b a = ρ y 1 , x 1 ρ y 1 , x 2 ρ y 1 , x a ρ y 2 , x 1 ρ y 2 , x 2 ρ y 2 , x a       ρ y b , x 1 ρ y b , x 2 ρ y b , x a
where Y b and Y a denote the reference vectors with b and a elements to be tested; M d Y b , X a denotes the Mahalanobis distance between vectors Y b and X a ; E d Y b , X a denotes the Ed between vectors Y b and X a ; C M b a denotes the correlation coefficient matrix between vectors Y b and X a ; and μ represents the parameters used to ensure that the sum is meaningful and the distance adjustment is on the same scale.
Due to the introduction of the correlation coefficient calculation in μ , μ reflects the adaptability of distance and the additivity on the distance scale [23].

5. ARP Experiment

5.1. Description of Experimental ARP Scene

In this study, we chose automatic parking garage No. 160 of the Dalian Shell Museum in Dalian, Xinghai Square, for the parking experiment. The near side and far side lines of the experimental garage are equal in length at 5 m long, the bottom line of the garage is 2.5 m long, and the width of the marker line is 0.1 m. The experimental vehicle model was a Toyota Corolla 1.2T S-CVT GL Pioneer Version (2019 model), with a length of 4635 mm, a width of 1780 mm, and a height of 1455 mm. The horizontal and longitudinal distances from the bottom right corner of the starting area of the experimental vehicle to the near corner point at the top of the garage were 2 m and 2.1 m, respectively.

5.2. Design and Configuration of ARP Experiment

The experimental vehicle used in this study is equipped with an ARP system, which has a real parking data acquisition sensor, a parking feasibility judgment device, a reference trajectory optimizer, a parking controller, a parking stop device, an emergency brake device, a parking image system, and a parking monitoring upper computer. The onboard sensors are mainly used for parking feasibility judgment and reference trajectory tracking control. In the parking process, the reference trajectory optimizer transmits the optimized reference trajectory to the parking controller and the upper parking monitoring computer, which not only provides the necessary reference trajectory for tracking control, but also enables safe parking, as the driver has a more comprehensive view of the real-time parking conditions. The ARP system exhibits powerful performance, especially regarding safety; once an accident risk is perceived, the emergency brake will be activated immediately. However, research on ARP is still in its early stages; security is critically important for this technology, and there is no room for error. The ARP experiment in this study was closely supervised by safety assurance personnel; the driver had a C2 driver’s license and acted as both a safety officer and an enhanced track control corrector, the ground safety inspector was able to contact the driver, and we employed an aerial drone photographer with good photography skills. A diagram of the design of the ARP experiment is shown in Figure 5.
As can be seen from Figure 5, in the ARP process, the ARP driver did not interfere with the normal parking tracking control, and would only activate the experimental vehicle’s stop or emergency brake function if necessary; some of the experimental results were directly obtained from the ARP process environment using an unmanned aerial camera, and some were indirectly obtained using ARP data acquisition sensors and an upper ARP monitoring computer.
The main parameters of the automatic parking experiment were as follows: The parking data acquisition, parking feasibility judgment, reference trajectory optimization, reference trajectory tracking control, emergency parking stop, and parking stop limit times were 3.5 s, 1.5 s, 12 s, 25 s, 0.4 s, and 1.5 s, respectively. One unmanned aerial camera (Mavic 2 DJI’s, ’Yu’, from Dajiang, Shenzhen, China) was used, and the upper computer for parking monitoring was a ’MacBook Pro 2016 Core i5 @ 2.9 GHZ’, which used Bluetooth for vehicle communication and was fixed to the mobile phone holder of the vehicle. Two MPC555LFMZP40 chips were used for the trajectory optimizer (used for parking data acquisition, parking feasibility judgment, and reference trajectory optimization) and tracking controller (used for reference trajectory tracking control, emergency parking stop, and parking stop), and two DSP28335 chips were used for the auxiliary control unit (ACU, used for the air-conditioner, overload protection device, communication, etc.) and actuator power unit (APU, used for the breaker, automobile steering gear, breaking/dynamical system, gear control system, etc.).

5.3. Experimental Results and Analysis of ARP

In order to verify the effectiveness of the optimization algorithm proposed in this study based on a real ARP scene in clear and windless weather, the IIMFO-GC proposed in this study, the IIMFO proposed in [13], the levy flight-based moth–flame optimization method proposed in [24], the MFO proposed in [25], the IPSO proposed in [4], and the PSO proposed in [26] were adopted to implement ARP reference trajectory optimization. Fuzzy proportional–integral–derivative (Fuzzy PID) control was adopted to implement parking tracking control for the reference trajectory. In order to achieve a fair comparison between the six algorithms, the same parameters were used for the optimization model as for the ARP reference trajectory, and are shown in Table 1. The parameters of IIMFO-GC for the ARP reference trajectory are shown in Table 2.
The experimental results for actual reverse parking given in this study include not only the optimization and tracking control results obtained using the parking monitoring host computer, but also the unmanned aerial photography images of several fixed points of the collected vehicle coverage areas during reverse parking. In order to achieve a fair comparison between the six algorithms, the same three fixed points of the collected vehicle coverage areas were selected: parking point P 4 followed by y ( P 4 ) [ 5.95 , 6.05 ] m, parking point P 7 followed by y ( P 7 ) [ 4.45 , 4.55 ] m, and parking stop point P 10 followed by y ( P 10 ) [ 2.30 , 2.50 ] m. The results obtained using the parking monitoring host computer for optimization and tracking control at fixed points of the collected vehicle coverage areas are shown in Figure 6 and Figure 7; the unmanned aerial images of the actual scene at fixed points of the collected vehicle coverage areas are shown in Figure 8; optimization and tracking control curves for the reference trajectories of automatic reverse parking are shown in Figure 9 and Figure 10; and tracking control trajectory error curves of the whole ARP process are shown in Figure 11. The position, path length, and berthing inclination angle results for automatic reverse parking optimization and tracking control are shown in Table 3 and Table 4. The tracking control errors in the position, path length, and berthing inclination angle results are shown in Table 5. The computational costs of each algorithm are shown in Table 6.
As can be seen from Table 3 and Figure 9, compared with the other five optimization algorithms, the IIMFO-GC proposed in this study exhibits a better reference trajectory; in particular, as shown in Figure 9, the reference trajectory path of IIMFO-GC is shorter. As can be seen from Table 4 and Figure 10, compared with the other five optimization algorithms, based on the same tracking control algorithm (Fuzzy PID), the IIMFO-GC proposed in this study obtains better tracking control results; in particular, as shown in Figure 10, the actual parking path of IIMFO-GC is shorter. As can be seen from Table 5 and Figure 11, compared with the other five optimization algorithms, based on the same tracking control algorithm (Fuzzy PID), the IIMFO-GC proposed in this study obtains tracking control results with fewer errors; in particular, as shown in Figure 11, the error in the tracking control of IIMFO-GC is smaller, and its tracking control trajectory errors for whole ARP process are fewer. Figure 7 shows that compared with the other five optimization algorithms, the IIMFO-GC proposed in this study exhibits a better reference trajectory; moreover, the optimized reference trajectory is closer to the near side line of parking garage and has a smaller expected berth inclination angle. As can be seen from Table 5, compared with the other five optimization algorithms, the IIMFO-GC proposed in this study has a better berthing inclination angle for both optimization and tracking control. Table 6 and Figure 6 and Figure 8 show that compared with the other five optimization algorithms, based on the same tracking control algorithm (Fuzzy PID), the IIMFO-GC proposed in this study obtains a better tracking control path; moreover, the tracking control path is closer to the near side line of the parking garage, and the actual berth inclination angle is smaller. As shown in Table 7, compared with the other five optimization algorithms, the IIMFO-GC proposed in this study is completed in a shorter time.

6. Conclusions

This paper proposes an IIMFO-GC for intelligent ARP function. Specifically, a novel reference trajectory optimization model for ARP is constructed, and a conversion principle for the ARP plane coordinate system is introduced; considering both the reference trajectory length and berthing inclination angle, a gene correction method for ARP is proposed for improving reference trajectory optimization performance. The proposed IIMFO-GC for the ARP trajectory optimizer effectively improves optimization performance.
The main contributions of this paper are as follows:
(I) Improvements to the ARP trajectory optimizer: The traditional reference trajectory optimization model is not suitable, and it is diffucult to establish effective boundary-crossing solutions for reverse parking. Therefore, in this study, based on the standard automatic parking plane system, a reasonable high-quality reference trajectory optimization model for ARP is constructed by combining the cubic spline-fitting method and a highly effective gene correction method.
(II) Improvements to the MFO algorithm: this paper proposes using IIMFO-GC to obtain a more appropriate reference trajectory for ARP; the following improvements are proposed: (II) integrate the immune mechanism with a high-quality solution-set maintenance mechanism based on fusion distance to strengthen global optimization quality; (I) introduce nonlinear decline strategies for weight coefficients and crossover and mutation probabilities to further strengthen global optimization ability.
These experimental results indicate that using IIMFO-GC for intelligent ARP function provides a better reference trajectory optimization effect than traditional optimization algorithms. Further, the comparison of the experimental results based on intelligent ARP function indicates that IIMFO-GC has a better control effect than anticipated.
However, there are several shortcomings in this research:
(I) Sensor accuracy has certain effects on ARP control.
(II) The effectiveness of the strategies for integrating the immune mechanism with moth–flame optimization is finite.
(III) The ARP scene used for the experimental comparison represents an ideal scenario, so the effectiveness of our experimental comparison is limited and one-sided.
We propose several future research directions for the use of MFO for ARP:
(I) Conduct further research on precise ARP trajectory data (such as position, length, and inclination angle) and measuring algorithms.
(II) Investigate existing improved MFOs, such as EMFO or Lévy-fight MFO, to further improve MFO optimization performance [18,24].
(III) Enhance the MFO algorithm by integrating other mechanisms or methods (such as the dual-population genetic mechanism) to further improve optimization quality [27].
(IV) Use several (more than two) typical ARP scenes for experimental comparison to further improve verification precision [28].

Author Contributions

The work presented here was performed in collaboration among all authors. G.L.: wrote this paper, completed the ARP experiment, and validated, investigated, edited, etc. X.X.: conceived the idea, and validated, investigated, data censored, writing—reviewed. L.W.: validated, data censored. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (62341313), the General Project of Natural Science Foundation of Inner Mongolia (2023MS06013), the Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region (NJYT-17-B34), Basic scientific research business expenses of colleges and universities directly under Inner Mongolia Autonomous Region - Plan to Improve the Scientific Research and Innovation Ability of Young Teachers (GXKY22125), Inner Mongolia Minzu University doctoral research initiation Fund Project (BS416), the Liaoning Provincial Department of Transportation Scientific Research Project (202344), the Liaoning Provincial Department of Education Scientific Research Project (JYTMS20230038).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

This research was supported by the Key Laboratory of Intelligent Manufacturing Technology in Inner Mongolia Minzu University. Thanks are due to Yang Guo, Zhang Wu and Guibo Qi for English polish work.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ARPautomatic reverse parking
IIMFO-GCimproved immune moth–flame optimization based on gene correction
MFOmoth–flame optimization

References

  1. Liu, H.; Luo, S.; Lu, J. Method for Adaptive Robust Four-Wheel Localization and Application in Automatic Parking Systems. IEEE Sens. J. 2019, 19, 10644–10653. [Google Scholar] [CrossRef]
  2. Liu, K.; Dao, M.; Inoue, T. An Exponentially epsilon-Convergent Control Algorithm for Chained Systems and Its Application to Automatic Parking Systems. IEEE Trans. Control Syst. Technol. 2006, 14, 1113–1126. [Google Scholar] [CrossRef]
  3. Suhr, J.K.; Jung, H.G. Automatic Parking Space Detection and Tracking for Underground and Indoor Environments. IEEE Trans. Ind. Electron. 2016, 63, 5687–5698. [Google Scholar] [CrossRef]
  4. Chen, Y.; Qian, Y. Improved Particle Swarm Optimization Algorithm for Automatic Entering Parking Space Based on Spline Theory. In Proceedings of the 2nd International Conference on Computing and Data Science (CONF-CDS 2021), Stanford, CA, USA, 28–30 January 2021; Volume 28, pp. 136–143. [Google Scholar]
  5. Gorinevsky, D.; Kapitanovsky, A. Neural network architecture for trajectory generation and control of automated car parking. IEEE Trans. Control Syst. Technol. 1996, 4, 50–56. [Google Scholar] [CrossRef]
  6. Chai, R.; Tsourdos, A.; Savvaris, A.; Chai, H.; Xia, Y.; Chen, C.L.P. Design and Implementation of Deep Neural Network-Based Control for Automatic Parking Maneuver Process. IEEE Trans. Neural Netw. Learn. Syst. 2022, 33, 1400–1413. [Google Scholar] [CrossRef] [PubMed]
  7. Nieto, R.M.; Garcia-Martin, A.; Hauptmann, A.G.; Martß′nez, J.M. Automatic Vacant Parking Places Management System Using Multicamera Vehicle Detection. IEEE Trans. Intell. Transp. Syst. 2019, 20, 1069–1080. [Google Scholar] [CrossRef]
  8. Laurini, M.; Consolini, L.; Locatelli, M. A graph-based algorithm for optimal control of switched systems: An application to car parking. IEEE Trans. Autom. Control 2021, 66, 6049–6055. [Google Scholar] [CrossRef]
  9. Maximilian, K.; Anil, K.M.; Adam, M.; Hasan, E.; Sandra, H. A Multi-Vehicle Control Framework With Application to Automated Valet Parking. IEEE Trans. Intell. Transp. Syst. 2021, 22, 5697–5707. [Google Scholar]
  10. Chai, R.; Tsourdos, A.; Savvaris, A.; Chai, S.; Xia, Y.; Chen, C.L.P. Multi-objective optimal parking maneuver planning of autonomous wheeled vehicles. IEEE Trans. Ind. Electron. 2020, 67, 10809–10821. [Google Scholar] [CrossRef]
  11. Yang, Y.; Xu, X.; Hong, J. Automatic Parking Path Planning of Tracked Vehicle Based on Improved A* and DWA Algorithms. IEEE Trans. Transp. Electrif. 2023, 9, 283–292. [Google Scholar] [CrossRef]
  12. Wang, X.; Shi, H.; Zhang, C. Path Planning for Intelligent Parking System Based on Improved Ant Colony Optimization. IEEE Access 2020, 8, 65267–65273. [Google Scholar] [CrossRef]
  13. Chen, Y.; Wang, L.; Liu, G.; Xia, B. Automatic Parking Path Optimization based on Immune Moth Flame Algorithm for Intelligent Vehicles. Symmetry 2022, 14, 1923. [Google Scholar] [CrossRef]
  14. Zhang, B.; Chen, S.; Gao, S.; Gao, Z.; Wang, D.; Zhang, X. Combination Balance Correction of Grinding Disk Based on Improved Quantum Genetic Algorithm. IEEE Trans. Instrum. Meas. 2023, 72, 1000112. [Google Scholar] [CrossRef]
  15. Liu, K.; Wang, X.; Wang, L.; Liu, G. Train Operation Strategy Optimization Based on Improved Genetic Algorithm. Int. J. Innov. Comput. Inf. Control 2019, 15, 1947–1965. [Google Scholar]
  16. Mirjalili, S. Moth-flame optimization algorithm: A novel nature-inspired heuristic paradigm. Knowl.-Based Syst. 2015, 89, 228–249. [Google Scholar] [CrossRef]
  17. Aziz, M.A.E.; Ewees, A.A.; Hassanien, A.E. Whale Optimization Algorithm and Moth-Flame Optimization for multilevel thresholding image segmentation. Expert Syst. Appl. 2017, 83, 242–256. [Google Scholar] [CrossRef]
  18. Xu, L.; Li, Y.; Li, K.; Beng, G.H.; Jiang, Z.; Wang, C.; Liu, N. Enhanced Moth-flame Optimization Based on Cultural Learning and Gaussian Mutation. J. Bionic Eng. 2018, 15, 751–763. [Google Scholar] [CrossRef]
  19. Qian, S.; Ye, Y.; Jiang, B.; Wang, J. Constrained Multiobjective Optimization Algorithm Based on Immune System Model. IEEE Trans. Cybern. 2016, 46, 2056–2069. [Google Scholar] [CrossRef] [PubMed]
  20. Parsopoulos, K.E.; Vrahatis, M.N. Recent approaches to global optimization problems through Particle Swarm Optimization. Nat. Comput. 2002, 1, 235–306. [Google Scholar] [CrossRef]
  21. Wang, L.; Wang, X.; Liu, G.; Li, Y. Improved Auto Disturbance Rejection Control Based on Moth Flame Optimization for Permanent Magnet Synchronous Motor. IEEJ Trans. Electr. Electron. Eng. 2021, 16, 1124–1135. [Google Scholar] [CrossRef]
  22. Liu, G.; Wang, X. Fault diagnosis of diesel engine based on fusion distance calculation. In Proceedings of the Advanced Information Management, Communicates, Electronic & Automation Control Conference, Xi’an, China, 3–5 October 2016; pp. 1621–1627. [Google Scholar]
  23. Liu, G.; Zhou, X.; Xu, X.; Wang, L.; Zhang, W. Fault Diagnosis of Diesel Engine Information Fusion based on Adaptive Dynamic Weighted Hybrid Distance-Taguchi method (ADWHD-T). Appl. Intell. 2022, 52, 10307–10329. [Google Scholar] [CrossRef]
  24. Li, Z.; Zhou, Y.; Zhang, S.; Song, J. Lévy-fight Moth-fame algorithm for function optimization and engineering design problems. Math. Probl. Eng. 2016, 2016, 1423930. [Google Scholar]
  25. Li, C.; Li, S.; Liu, Y. A least squares support vector machine model optimized by moth-flame optimization algorithm for annual power load forecasting. Appl. Intell. 2016, 45, 1166–1178. [Google Scholar] [CrossRef]
  26. Kennedy, J.; Eberhart, R.C. Particle swarm optimization. In Proceedings of the ICNN’95–International Conference on Neural Networks, Perth, WA, Australia, 27 November–1 December 1995; Volume 4, pp. 1942–1948. [Google Scholar]
  27. Che, S.; Chen, Y.; Wang, L.; Xu, C. Electric Vehicle Ordered Charging Planning Based on Improved Dual-Population Genetic Moth-Flame Optimization. Algorithms 2024, 17, 110. [Google Scholar] [CrossRef]
  28. Zhang, P.; Xiong, L.; Yu, Z.; Fang, P.; Yan, S.; Yao, J.; Zhou, Y. Reinforcement Learning-Based End-to-End Parking for Automatic Parking System. Sensors 2019, 19, 3996. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic diagram of the ARP principle.
Figure 1. Schematic diagram of the ARP principle.
Sensors 24 02270 g001
Figure 2. Schematic diagram of gene correction for collision avoidance on the far side line of the parking garage.
Figure 2. Schematic diagram of gene correction for collision avoidance on the far side line of the parking garage.
Sensors 24 02270 g002
Figure 3. Schematic diagram of gene correction for berthing inclination.
Figure 3. Schematic diagram of gene correction for berthing inclination.
Sensors 24 02270 g003
Figure 4. Schematic diagram of gene correction for berthing dislocation.
Figure 4. Schematic diagram of gene correction for berthing dislocation.
Sensors 24 02270 g004
Figure 5. Diagram showing the design of the ARP experiment.
Figure 5. Diagram showing the design of the ARP experiment.
Sensors 24 02270 g005
Figure 6. Results obtained using parking monitoring host computer for optimization of fixed points of collected vehicle coverage areas.
Figure 6. Results obtained using parking monitoring host computer for optimization of fixed points of collected vehicle coverage areas.
Sensors 24 02270 g006
Figure 7. Results obtained using the parking monitoring host computer for tracking control of fixed points of collected vehicle coverage areas.
Figure 7. Results obtained using the parking monitoring host computer for tracking control of fixed points of collected vehicle coverage areas.
Sensors 24 02270 g007
Figure 8. Unmanned aerial images of actual scene for fixed points of collected vehicle coverage areas.
Figure 8. Unmanned aerial images of actual scene for fixed points of collected vehicle coverage areas.
Sensors 24 02270 g008
Figure 9. Optimization curves for the reference trajectories of ARP.
Figure 9. Optimization curves for the reference trajectories of ARP.
Sensors 24 02270 g009
Figure 10. Tracking control curves for the reference trajectories of ARP.
Figure 10. Tracking control curves for the reference trajectories of ARP.
Sensors 24 02270 g010
Figure 11. Tracking control trajectory error curves of whole ARP process.
Figure 11. Tracking control trajectory error curves of whole ARP process.
Sensors 24 02270 g011
Table 1. Parameters of optimization model for ARP reference trajectory.
Table 1. Parameters of optimization model for ARP reference trajectory.
NameSymbolValue
Starting inclination angle threshold Δ B max π 75
Vehicle coverage area collection times n s 10
Parking inclination angle threshold Δ S max π 20
Parking position error threshold Δ y n s 0.15 m
Parking process period Δ t 1.5 × 10 3 s
Attitude angle difference threshold Δ max 5 × 10 4 π
Parking interval threshold Δ L max 10 3
Parking position expected ordinate y H 2.35 m
Table 2. Parameter configuration of IIMFO-GC for ARP reference trajectory.
Table 2. Parameter configuration of IIMFO-GC for ARP reference trajectory.
NameSymbolValue
Moth population sizen30
Iteration number t max 80
Logarithmic spiral shape adjustment coefficient τ π 20
Minimum weight coefficient ω min 0.4
Maximum weight coefficient ω max 0.9
Weight coefficient optimization factor β 1.85
Minimum crossover and mutation probabilities p c , min , p m , min 0.4,0.1
Maximum crossover and mutation probabilities p c , max , p m , max 0.8,0.3
Crossover probability optimization factor β m 1.66
Mutation probability optimization factor β c 1.49
Collision correction maximum ordinate offset Δ Y f c max 0.92 m
Collision correction empirical coefficient c f c 0.58
Berthing inclination maximum ordinate offset y H 2.35 m
Table 3. Position and path length results of reference trajectory optimization for ARP.
Table 3. Position and path length results of reference trajectory optimization for ARP.
AlgorithmsCoordinates of p 10 , p 7 , p 4 (m)Path Length (m)
IIMFO-GC(1.52,2.29), (1.55,4.43), (1.84,5.97)9.480
IIMFO(1.37,2.29), (1.44,4.43), (1.78,5.96)9.659
LMFO(1.02,2.30), (1.20,4.44), (1.69,5.91)9.694
MFO(1.29,2.31), (1.37,4.48), (1.80,5.87)9.619
IPSO(1.21,2.43), (1.29,4.46), (1.64,5.90)9.732
PSO(1.14,2.27), (1.32,4.40), (1.61,5.91)9.829
Table 4. Position and path length results of reference trajectory tracking control for ARP.
Table 4. Position and path length results of reference trajectory tracking control for ARP.
AlgorithmsCoordinates of p 10 , p 7 , p 4 (m)Path Length (m)
IIMFO-GC(1.51,2.29), (1.50,4.42), (1.90,5.96)9.493
IIMFO(1.39,2.30), (1.38,4.43), (1.74,5.95)9.658
LMFO(1.16,2.29), (1.30,4.45), (1.80,5.89)9.741
MFO(1.28,2.31), (1.52,4.39), (1.84,5.88)9.635
IPSO(1.31,2.27), (1.34,4.42), (1.70,5.88)9.728
PSO(1.27,2.25), (1.40,4.50), (1.66,5.91)9.785
Table 5. Berthing inclination angle results of reference trajectory optimization and tracking control for ARP.
Table 5. Berthing inclination angle results of reference trajectory optimization and tracking control for ARP.
AlgorithmsOptimizationTracking Control
IIMFO-GC0.01690.0177
IIMFO0.01900.0204
LMFO0.01810.0193
MFO0.11130.0530
PSO0.10930.1085
PSO0.10980.1078
Table 6. Tracking control errors in position, path length, and berthing inclination angle results for ARP.
Table 6. Tracking control errors in position, path length, and berthing inclination angle results for ARP.
AlgorithmsCoordinates of p 10 , p 7 , p 4 (m)Path Length (m)Berthing Inclination Angle
IIMFO-GC(0.01,0), (0.05,0.01), (0.06,0.01)0.0130.0008
IIMFO(0.02,0.01), (0.06,0), (0.04,0.01)0.0010.0014
LMFO(0.14,0.01), (0.10,0.01), (0.11,0.02)0.0470.0012
MFO(0.01,0), (0.15,0.09), (0.04,0.01)0.0160.0583
IPSO(0.10,0.16), (0.05,0.04), (0.06,0.02)0.0040.0008
PSO(0.13,0.02), (0.08,0.10), (0.04,0)0.0440.0020
Table 7. Computational cost of each algorithm performed for ARP.
Table 7. Computational cost of each algorithm performed for ARP.
AlgorithmsComputational Cost (s)
IIMFO-GC9.335
IIMFO9.754
LMFO10.108
MFO10.227
PSO9.935
PSO10.040
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, G.; Xu, X.; Wang, L. Improved Immune Moth–Flame Optimization Based on Gene Correction for Automatic Reverse Parking. Sensors 2024, 24, 2270. https://doi.org/10.3390/s24072270

AMA Style

Liu G, Xu X, Wang L. Improved Immune Moth–Flame Optimization Based on Gene Correction for Automatic Reverse Parking. Sensors. 2024; 24(7):2270. https://doi.org/10.3390/s24072270

Chicago/Turabian Style

Liu, Gang, Xinli Xu, and Longda Wang. 2024. "Improved Immune Moth–Flame Optimization Based on Gene Correction for Automatic Reverse Parking" Sensors 24, no. 7: 2270. https://doi.org/10.3390/s24072270

APA Style

Liu, G., Xu, X., & Wang, L. (2024). Improved Immune Moth–Flame Optimization Based on Gene Correction for Automatic Reverse Parking. Sensors, 24(7), 2270. https://doi.org/10.3390/s24072270

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop