Urban Sprawl Simulation Mapping of Urmia (Iran) by Comparison of Cellular Automata–Markov Chain and Artiﬁcial Neural Network (ANN) Modeling Approach

: Considering urbanization can lead to irreversible land transformations, it is crucial to provide city managers, environmental resources managers, and even people with accurate predicting land use/land cover (LULC) to accomplish sustainable development goals. Although many methods have been used to predict land use/land cover (LULC), few studies have compared them. Therefore, by analyzing the results of various prediction models and, consequently, recognizing the most accurate and reliable ones, we can assist city managers, environmental resources managers, and researchers.. In this regard, this research compares Cellular Automata–Markov Chain and Artiﬁcial Neural Network (ANN) as frequently used models to overcome this gap and help those concerned about sustainable development to predict urban sprawl with the most reliable accuracy. In the ﬁrst step, Landsat satellite images acquired in 2000, 2010, and 2020 were classiﬁed with Maximum Likelihood Classiﬁcation (MLC), and LULC maps were prepared for each year. In the second step, to investigate the LULC prediction, validation of the CA–Markov and ANN methods was performed. In this way, the LULC simulation map of 2020 was prepared based on the LULC map of 2000 and 2010; next, the predicted LULC map of 2020 and the actual LULC map for 2020 were compared using correctness, completeness, and quality indices. Finally, the LULC map for 2030 was generated using both algorithms, and the corresponding change map was extracted, showing a reduction in soil and vegetation areas (respectively, 39% and 12%) and an expansion (58%) in built-up regions. Moreover, the validation test of the methods showed that the two algorithms were closer to each other; however, ANN had the highest completeness (96.21%) and quality (93.8%), while CA–Markov had the most correctness (96.47%). This study showed that the CA–Markov algorithm is more accurate in predicting the future of larger areas with higher allocations (urban and vegetation cover) while the ANN algorithm is more accurate in predicting the future of small areas with fewer allocations (soil and rock).


Introduction
The concept of sustainability refers to the long-term preservation of ecosystems. In other words, sustainable development emphasizes meeting the needs of today's people without endangering the future in meeting their needs [1]. Sustainability includes five dimensions: public health, environmental quality, economic vitality, countermeasures for areas [12]; hence, in these countries, urbanization and associated demographic changes pose unprecedented challenges in terms of hunger, food insecurity, and malnutrition [14]. Urbanization is a change process from a rural society to a civilized urban life with dense population and migration [42], affecting central sustainable policies such as water supply, green areas provision, and transportation [6]. Moreover, from a social perspective, it also affects urban sustainability [5] through marginalization, crime [43], and poverty [44].
Remote sensing is a rigorous and practical surface monitoring tool specifically used for creating maps of LULC [45][46][47], observing land surfaces, and extracting data [48]. In addition, compared to measurements taken only from a specific location, it provides large-scale, high-resolution, continuous information for LULCs [48,49]. Therefore, LULC classification based on remote sensing has an essential role in evaluating the results of management interventions and how the changes will occur in the future [50].
Changing natural environments, agricultural lands, etc., to urban areas is one of the most demanding environmental challenges in every country [51], which underlines the importance of spatial predictions to forecast and manage regional and global changes and the ability to improve environmental sustainability [52]. While some patterns are based on the prediction rate of change, others rely on spatial patterns that focus on required data and validation strategies [53]. In recent years, some spatial models have combined remote sensing (RS) and geographic information systems (GISs) to simulate and predict future scenarios of LULC [54], for instance, Markov chain [55,56], cellular automata-Markov (CA-Markov) [57], logistic regression [58,59], cellular automata model [60,61], SLEUTH model [62,63], and artificial neural network model [64,65]. Many studies in the literature have also separately applied the Fractal [66,67], CA-Markov [68][69][70][71] and ANN [72][73][74], logistic regression [75,76], and agent-based [77,78] models, which are popular in prediction studies for LULC. Among these, Markov chain analysis (MCA), cellular automata (CA), cellular automata-Markov model (CA-Markov), artificial neural network (ANN), binary logistic regression, and the fractal model can be considered the most common [79]. However, despite these methods being used to predict land use/land cover (LULC), few studies have compared them to recognize the most reliable results. Moreover, since each research can only evaluate a small area of science, this research intends to compare Cellular Automata-Markov Chain and Artificial Neural Network (ANN) as frequently used models to mitigate this problem and provide various stakeholders with more reliable and accurate models. More specifically, this paper seeks to answer two basic questions: (1) Which model has more validity for predicting LU/LC? (2) Which LU/LC does each model predict more reliably?
For this purpose, the LULC map in 2030 was projected with the help of changes in LULC in 2000, 2010, and 2020. However, before that, the 2020 map was predicted (with both models) and compared with the 2020 map prepared by using the maximum likelihood method. The accuracy of the produced model was investigated with different validation methods, including correctness, completeness, and quality indexes.

Study Area
The city of Urmia (37 • 33 N, 45 • 04 E) is the center of the West Azerbaijan province in the northwest of Iran, is located at a distance of 18 km from Lake Urmia [80], and is spread over an area of 7548 ha. The climate of this region is cold and semiarid, with an average of 360 mm annual precipitation and 11 • C annual temperature [81], making it one of the coldest cities in Iran [82]. Urmia is one of Iran's most important historical and growing metropolises and has grown considerably in recent years [83]. According to the 2016 census, it had a population of about 750,000 [84,85], which according to the 2018 census, has increased to 800,000 [86]. Rapid urbanization has not only led to the reduction of plants and biodiversity in the city, but has also been responsible for environmental and climate changes [87,88]. Figure 1 shows the study area location. metropolises and has grown considerably in recent years [83]. According to the 2016 census, it had a population of about 750,000 [84,85], which according to the 2018 census, has increased to 800,000 [86]. Rapid urbanization has not only led to the reduction of plants and biodiversity in the city, but has also been responsible for environmental and climate changes [87,88]. Figure 1 shows the study area location.

Data and Methods
Evaluating the rate of changes from one phase to another in a specific time and place of spatial data is significant for predicting future change scenarios [89]. Using noncommercial satellite images is an inexpensive and fast method of predicting LULC, and is a rigorous tool for land planners. The Landsat satellite series provides the longest record of satellite observations. Accordingly, Landsat is a precious source for monitoring global changes and observing the earth with medium-spatial resolution in decision-making procedures.
This study retrieved Landsat images (7 and 8) covering Urmia on 03/06/2000, 30/05/2010, and 02/06/2020. The images were downloaded from the USGS website. Initially, these images were pr-processed for geometric and radiometric corrections to make them suitable for information extraction. Then, in the categorization stage, , high-resolution images (such as Google Earth and World Imagery), normalized difference built-up index (NDBI), and normalized difference vegetation index (NDVI) were used (sample points size are in Table 1) for choosing sample points from the old maps. In order to extract the LULC map, this study applied MLC due to its highest quality results compared to others [74]. Finally, predictions for the year 2030 were carried out with the methods detailed below, and their accuracy was investigated.

Data and Methods
Evaluating the rate of changes from one phase to another in a specific time and place of spatial data is significant for predicting future change scenarios [89]. Using noncommercial satellite images is an inexpensive and fast method of predicting LULC, and is a rigorous tool for land planners. The Landsat satellite series provides the longest record of satellite observations. Accordingly, Landsat is a precious source for monitoring global changes and observing the earth with medium-spatial resolution in decision-making procedures.
This study retrieved Landsat images (7 and 8) covering Urmia on 03/06/2000, 30/05/2010, and 02/06/2020. The images were downloaded from the USGS website. Initially, these images were pr-processed for geometric and radiometric corrections to make them suitable for information extraction. Then, in the categorization stage, high-resolution images (such as Google Earth and World Imagery), normalized difference built-up index (NDBI), and normalized difference vegetation index (NDVI) were used (sample points size are in Table 1) for choosing sample points from the old maps. In order to extract the LULC map, this study applied MLC due to its highest quality results compared to others [74]. Finally, predictions for the year 2030 were carried out with the methods detailed below, and their accuracy was investigated.  [26,90,91]. CA model behaviors are influenced by the uncertainty stemming from the interaction between the model elements, structures, and quality of data sources that are considered model input. This model is often focused on the local interaction of the local cells, distinct spatial and temporal properties, and the rigorous computational ability of space that is proper for dynamic simulation and displays with self-made features. Thus, the CA model can be described as follows (Equation (1)): S is a collection of limited and discrete cellular models, N is a cellular field, t and t + 1 stand for different times, and f is the transformation of cellular states law in local space.
Markov chain is insufficient to actively simulate or predict LULC because it does not consider the spatial distribution in each land category or the growth direction [42,92]. CA-Markov is a Cellular Automata, Multiobjective Land Allocation predicting method that adds elements of spatial proximity and the knowledge of spatial distribution to the Markov chain [93][94][95][96][97][98].

Artificial Neural Network
To detect the probability of LULC transformation, ANN uses several output neurons to simulate LULC changes. In the first step, the neural network inputs are defined for the simulation. The simulation is cellular-based (pixel-based), and each cell has a set of natural features (spatial variables) as input to the neural network, defined as follows: where x i is the i-th property and T is the transition. Each correlation between the spatial variables is evaluated by the mutual comparison of the raster, choosing the first raster from one variable and the second from the other. Next, the LULC region and the changes related to each group are measured between the initial and final times. In the next step, the probability of transformation through ANN is simulated. The neural network structure is created with three layers: the input layer, the latent layer, and the output layer. In the latent layer, the received signal by the j-th cell and net j (k,t) is calculated from the input layer for the k-th cell at t time and is defined as follows: w i,k is the weight between the input and latent layers, and x i (k,t) is the i-th saleable attribute related to i-th neuron in the latent layer according to the k-th cell at t time. The probability of relocation, with the performance of the output of a neural network being considered, is calculated as follows: P(k,t,l) is the probability of changing from the current state of l LULC for k cell in t time, and w j,l is the weight between the latent and output layers [99].

Model Assessment
In this study, indices, including completeness, correctness, and quality, were used to quantitatively assess the forecast findings of these methods. The details of the indices used in the study are described below.
This index of completeness describes what percentages of the features shown in the source data are considered in the result. In this index, the feature units related to other features and distinguished wrongly have no impact on the value of this index. Therefore, this index is defined as follows: The correctness index is used for the correctness of classification. It refers to the percentage of the features that are detected in the results that are the same as the reference features. In this index, feature units that exist in the source data but that were not distinguished in the result do not influence the value of this index. This index is defined as follows: The quality that pertains to the evaluation of the findings of both correctness and completeness and is thus defined as follows: Quality = (TP/TP + FP + FN) * 100 (7) True positive (TP) is the number of units of features that exist in the source data and in the findings or the number of features that have been successfully detected correctly as a feature. False positive (FP) is the number of features that do not exist in the source data but have been identified in the results as a feature. False negative (FN) is the number of negative features that exist in source data but have not been identified in the results [100][101][102]. Figure 2 shows the research flowchart.
The correctness index is used for the correctness of classification. It refers to the percentage of the features that are detected in the results that are the same as the reference features. In this index, feature units that exist in the source data but that were not distinguished in the result do not influence the value of this index. This index is defined as follows: Correctness = (TP / FP TP) * 100 (6) The quality that pertains to the evaluation of the findings of both correctness and completeness and is thus defined as follows: True positive (TP) is the number of units of features that exist in the source data and in the findings or the number of features that have been successfully detected correctly as a feature. False positive (FP) is the number of features that do not exist in the source data but have been identified in the results as a feature. False negative (FN) is the number of negative features that exist in source data but have not been identified in the results [100][101][102]. Figure 2 shows the research flowchart.

Results
This section presents the findings of the methods used. Firstly, the findings of the classification with MLC are presented. Figure 3 shows the LULC maps obtained as a result of classification for 2000, 2010, and 2020.

Results
This section presents the findings of the methods used. Firstly, the findings of the classification with MLC are presented. Figure 3 shows the LULC maps obtained as a result of classification for 2000, 2010, and 2020. Before using the achieved maps as the input of the LULC prediction methods, we needed verify each class's classification. Table 2 depicts the result of the verification for the classification of each class in the three studied times. According to Table 2, the classification results in each of the three times appear suitable for all classes. The Kappa coeffi- Before using the achieved maps as the input of the LULC prediction methods, we needed verify each class's classification. Table 2 depicts the result of the verification for the classification of each class in the three studied times. According to Table 2, the classification results in each of the three times appear suitable for all classes. The Kappa coefficient in all three times is higher than 90%, and in the 2020 map, it is approximately 100%. Although the accuracy for soil was less than 90% in 2000 and 2010, it is acceptable due to being above 83%. For accuracy assessment, the prepared maps were compared with the samples seen from Google Earth.   According to Figure 4, the rocky areas and vegetation changed little over the 20 years, but soil and built-up areas changes extensively. The soil class shrunk from 5700 hectares to 4650 in 2010 and to 3450 hectares in 2020. According to Figure 4, rocky areas increased in 2010 (the reduction of vegetation in rocky areas that had more vegetation in 2000; this could be due to logging in mountainous areas or decreased rainfall), but it decreased again in 2020. Figure 5 shows the comparison of the LULC for the 2020 map (a) obtained as a result of the classification with LULC prediction maps for 2020 as produced by the CA-Markov (b) and the ANN (c) method. According to Figure 4, the rocky areas and vegetation changed little over the 20 years, but soil and built-up areas changes extensively. The soil class shrunk from 5700 hectares to 4650 in 2010 and to 3450 hectares in 2020. According to Figure 4, rocky areas increased in 2010 (the reduction of vegetation in rocky areas that had more vegetation in 2000; this could be due to logging in mountainous areas or decreased rainfall), but it decreased again in 2020. Figure 5 shows the comparison of the LULC for the 2020 map (a) obtained as a result of the classification with LULC prediction maps for 2020 as produced by the CA-Markov (b) and the ANN (c) method. to 4650 in 2010 and to 3450 hectares in 2020. According to Figure 4, rocky areas increased in 2010 (the reduction of vegetation in rocky areas that had more vegetation in 2000; this could be due to logging in mountainous areas or decreased rainfall), but it decreased again in 2020. Figure 5 shows the comparison of the LULC for the 2020 map (a) obtained as a result of the classification with LULC prediction maps for 2020 as produced by the CA-Markov (b) and the ANN (c) method.  In the LULC classes shown in Figure 5, we can see that the appearance of the CA-Markov and ANN are pretty similar to the current situation. To better understand and compare the methods, the results of completeness, accuracy, and quality indicators are shown in Table 3.  Table 3 presents the completeness, correctness, and quality of the validation methods in order to evaluate and compare the outcomes of the LULC prediction with MLC for 2020. The blue numbers shown in the table indicate the most promising results. Accordingly, although the accuracy of the CA-Markov method was close to that of ANN in the LULC prediction models, the best results were obtained with the ANN technique. The ANN technique map had the lowest vicinity of rocky lands and, despite having the lowest completeness, had the highest correctness and forecast quality.
In the case of built-up areas, the CA-Markov method had the highest completeness and quality, while the ANN method had the highest correctness. For soil, the highest completeness and quality belonged to the ANN method, and the highest correctness belonged to the CA-Markov method. In the case of vegetation, the CA-Markov method had the best indices. In terms of average (all LU/LCs), the highest completeness and quality were achieved with ANN method, and the highest correctness with the CA-Markov method.
In conclusion, the results of the indices and subsequently, the typical mean of the two algorithms, ANN and CA-Markov, are very close to each other, but the ANN technique had the perfect mean in the two indices of completeness and quality, and the CA-Markov algorithm had the best correctness. Therefore, the urban and vegetation LU/LCs in the CA-Markov algorithm and the soil and rock LCs in the ANN algorithm were better simulated, as can be seen in Table 3. The most significant results of this research are that the CA-Markov algorithm is better able to predict those phenomena with wider and more continuous surfaces, while the ANN algorithm performs better in simulating phenomena that are smaller areas on the map and that include a lower percentage of the area. Figure 6 shows the LULC prediction map for 2030 of the CA-Markov and ANN (first row) and the map of changes in 2020-2030 (bottom row) as obtained by comparing the map predicted by the algorithms for 2030 and the LULC map for 2020 using MLC.  Figure 6 show that CA-Markov could reveal more changes compared to ANN. In both methods, the east side of the study region changes more than does the west side. To depict the degree of the LULC change, Figure 6 shows the area of each LULC in 2030. Figure 7 shows that, considering that the built-up LU area in this study will reach more than 9000 hectares in 2030, which is less than the 6000 hectares in 2000, the city of Urmia will extend by 50% over the next 30 years. Table 3 shows the area and the map of LULCs of urban (built-up) and vegetation as generated through CA-Markov, as well as the map and area of soil and rocky lands as generated through ANN. Figure 7 also shows that LU/LCs are similar in the area and more reliable in evaluation parameters (smaller zones in ANN and larger zones in CA-Markov).  Figure 6 show that CA-Markov could reveal more changes compared to ANN. In both methods, the east side of the study region changes more than does the west side. To depict the degree of the LULC change, Figure 6 shows the area of each LULC in 2030. Figure 7 shows that, considering that the built-up LU area in this study will reach more than 9000 hectares in 2030, which is less than the 6000 hectares in 2000, the city of Urmia will extend by 50% over the next 30 years. Table 3 shows the area and the map of LULCs of urban (built-up) and vegetation as generated through CA-Markov, as well as the map and area of soil and rocky lands as generated through ANN. Based on the two methods, the built-up environment appears to be growing, resulting in the decline of the vegetation area. The most significant difference in area between LU/LCs derived from the 2030 forecast algorithms is related to soil LC. It can be understood from the validation results of Table 3 that the precision, quality, and completeness statistics of this LU are lower than those of other LU/LCs, representing a significant statistical discrepancy between the algorithms.

Discussion
Today, the world is facing the challenge of environmental instability, which is the product of human activity. For example, one of the environmental concerns in the suburbs of Urmia is the drying up of a significant part of its lake [103][104][105]. Urban sprawl causes irreversible changes in the earth's surface because even with the destruction of built-up areas, achieving a sustainable environment (natural environment) or one identical to what had been before construction is impossible. Therefore, managing urban sprawl is the most practical way to cope with this issue and can be accomplished by having reliable maps of the city expansion forecast. In this regard, various prediction methods such as CA-Markov and ANN have been employed in many regions all over the world, including the city of Thimphu, in Bhutan [106], the Majang Forest Biosphere Reserves of Gambella in southwestern Ethiopia [107], Orkhon Province in Mongolia [108], and the northeast of Iran [109]. However, unlike in these studies, in our study, which was conducted in the city of Urmia, two popular LULC prediction methods widely used in the literature were implemented and compared in a single study. Thus, this study comparing LULC prediction methods in Urmia significantly contributes to the literature.
This investigation indicated that the Urmia urban area has expanded too greatly [110]. It indicates that the built-up areas will occupy many natural and agricultural resources and that this land change situation threatens natural life [111]. It also confirmed the results of Zare-Naghadehi et al.'s study [74]. Table 2 shows that the evaluation of the accuracy of LULC maps is very suitable (except for soil in the 2000 and 2010 LULC maps). Moreover, comparing the methods revealed that both methods have a high capability to predict LULC; however, according to the evaluation criteria, ANN performed slightly better than did CA-Markov. Based on the results of this research, the CA-Markov algorithm is more accurate in predicting urban and vegetation areas, while the ANN algorithm is more accurate in predicting soil and rock cover. Based on the two methods, the built-up environment appears to be growing, resulting in the decline of the vegetation area. The most significant difference in area between LU/LCs derived from the 2030 forecast algorithms is related to soil LC. It can be understood from the validation results of Table 3 that the precision, quality, and completeness statistics of this LU are lower than those of other LU/LCs, representing a significant statistical discrepancy between the algorithms.

Discussion
Today, the world is facing the challenge of environmental instability, which is the product of human activity. For example, one of the environmental concerns in the suburbs of Urmia is the drying up of a significant part of its lake [103][104][105]. Urban sprawl causes irreversible changes in the earth's surface because even with the destruction of built-up areas, achieving a sustainable environment (natural environment) or one identical to what had been before construction is impossible. Therefore, managing urban sprawl is the most practical way to cope with this issue and can be accomplished by having reliable maps of the city expansion forecast. In this regard, various prediction methods such as CA-Markov and ANN have been employed in many regions all over the world, including the city of Thimphu, in Bhutan [106], the Majang Forest Biosphere Reserves of Gambella in southwestern Ethiopia [107], Orkhon Province in Mongolia [108], and the northeast of Iran [109]. However, unlike in these studies, in our study, which was conducted in the city of Urmia, two popular LULC prediction methods widely used in the literature were implemented and compared in a single study. Thus, this study comparing LULC prediction methods in Urmia significantly contributes to the literature.
This investigation indicated that the Urmia urban area has expanded too greatly [110]. It indicates that the built-up areas will occupy many natural and agricultural resources and that this land change situation threatens natural life [111]. It also confirmed the results of Zare-Naghadehi et al.'s study [74]. Table 2 shows that the evaluation of the accuracy of LULC maps is very suitable (except for soil in the 2000 and 2010 LULC maps). Moreover, comparing the methods revealed that both methods have a high capability to predict LULC; however, according to the evaluation criteria, ANN performed slightly better than did CA-Markov. Based on the results of this research, the CA-Markov algorithm is more accurate in predicting urban and vegetation areas, while the ANN algorithm is more accurate in predicting soil and rock cover.
The first thing that should be considered for a reliable prediction is the accuracy of the input data. Table 2 shows that the average accuracy of the built-up class classification is 97.47%, which is high. Furthermore, by examining other classes, it is clear that the accuracy of classification with MLC is quite suitable as an input for predicting changes (only the accuracy of the soil class in the 2000 and 2010 maps is less than 90%).
Examining the change and transformation of LULC ( Figure 4) shows that over 20 years (2000-2020), the area of the city grew by nearly 60%, which means the loss of other classes of LULCs. During this time, bare soil decreased to about 61% of its area (a loss of 2222 hectares during the study). Additionally, more than 800 hectares of vegetation cover were lost, meaning that there is low environmental stability in the area. In short, Table 3 presents the main results of this study: three indicators were used to assess the effectiveness of the CA-Markov and ANN methods and revealed that both algorithms have high capabilities for LULC prediction. Although, according to the results in Table 3, the more accurate prediction was achieved by ANN, the CA-Markov algorithm's performance regarding urban sprawl, which is the focus of this study, performed better than did ANN in two indicators (completeness and quality) for predicting build-up.

Conclusions
Urban sprawl can prevent sustainable development goals and cause considerable impairment, especially in the suburbs, which are more vulnerable to changes because of human exploitation. Moreover, predicting the urban development pattern and the region's changes can be crucial for urban investors, residents, and landowners, highlighting the importance of LULCC analysis. Therefore, the timely and efficient decisions are essential but impossible to make without accurate information. More specifically, a proper LULC map will pave the road for local planners and decision-makers to monitor environmental hazards such as soil erosion, flooding, landslides, and degradation of pastures. To this end, Landsat images with a long imaging record can be an excellent archive to research the shift in an area's LULC. Furthermore, various algorithms exist to forecast LU and are slightly or generally different from each other in their structures. Knowing which algorithm functions better can help researchers use the most efficient one in their future research.
Due to its high development rate in recent years, this study was conducted in Urmia to evaluate and compare two common algorithms (CA-Markov and ANN) and recognize which has the better performance. Consequently, researchers and authorities can use the most efficient one in future research and projects. For this purpose, satellite images were used to assess four land-use (constructed, rocky areas, vegetation cover, and soil) in Urmia in 2000, 2010, and 2020, and their LULC images were obtained. Then, the results of the CA-Markov and ANN algorithms were compared, indicating that all parameters in the two algorithms were higher than 90% (96.47 correctness in the CA-Markov algorithm and 96.21 completeness and 93.8 quality in the ANN algorithm) and that these algorithms have excellent prediction ability. Preliminary findings from the LULC map indicate that the city's growth was positive, while on the other hand, the areas covered by soil and vegetation were limited, and the rocky areas almost remained unchanged. The CA-Markov algorithm had the highest completeness for rock cover, and the ANN algorithm, which shows rocky terrain with less accuracy, had low completeness but high correctness and quality. On the other hand, in soil land cover, the ANN method had the highest completeness and quality, while the CA-Markov algorithm had the highest correctness. However, in vegetation, all three statistics pointed to the superiority of the CA-Markov method. Regarding the average criteria used in all LULCs, the CA-Markov algorithm had the highest completeness, and the ANN algorithm had the highest correctness and quality.
In summary, both CA-Markov and the ANN algorithms perform well in a way that one cannot decide which algorithm excels over the other, and thus more testing is required. However, the limitation of the study is that there is insufficient research comparing different prediction models in other regions to achieve a more comprehensive view of their efficiency and accuracy. Therefore, although the CA-Markov method is more effective in predicting LULCs in extensive areas and the ANN method is more effective for smaller areas, we suggest that other researchers evaluate other algorithms in other regions to better understand these prediction models. We also recommend that researchers compare the results from other methods, such as the Genetic Algorithm Optimized Neural Network Model [112], CA-Based SLEUTH [113], and CycleGANs-based CNN [114], using several evaluation indices.