Landscape Ecological Analysis of Green Network in Urban Area Using Circuit Theory and Least ‐ Cost Path

: Quantitative securing of green space in already developed cities has many practical limi ‐ tations due to socio ‐ economic limitations. Currently, South Korea is planning a green network to secure and inject effective green space, but it is difficult to reflect it in the actual space plan due to the abstract plan. This study utilizes circuit theory and least ‐ cost path methods for presenting a green network that is objectified and applicable to spatial planning. First, an analysis of the Least ‐ cost Path revealed 69 least ‐ cost paths between 43 core green areas of the study site. Most least ‐ cost paths have been identified as passing through small green areas and streams in the city. Using the circuit theory, it was also possible to distinguish areas other than least ‐ cost paths from areas with high potential for development, areas where target species are concentrated within corridors. In particular, areas with relatively high green network improvement effects were derived within and around corridors. This study is most significant in establishing and evaluating existing urban green networks, overcoming the limitations discussed at the linear level and expanding to the area level. To increase the utilization of this study in the future, field surveys and monitoring studies on target species need to be supplemented.


Introduction
Over the past few decades, Korea has expanded its gray infrastructure, such as roads railways and buildings, due to economic growth centered on cities. On the contrary, the city's external expansion has continuously reduced its Green Infrastructure [1][2][3]. In particular, the large green space and connection corridor of the exurban, which supplies biodiversity to inner urban areas, are threatened with reduction and fragmentation [4,5]. To secure ecological quality by reflecting this trend of the times, discussions on quantitative securing of green space in the city are actively being held [6][7][8]. However, since the quantitative increase in green areas in overcrowded urban spaces leads to realistic limitations, research on green networks is being actively conducted to distinguish between relatively valuable green spaces and areas that require preferential input [9][10][11][12].
From a landscape ecological perspective, the concept of green networks began with the idea of how spatial structural features such as the arrangement, shape and size of green land would interconnect to species [13]. Among several features, inner green adjacency and connectivity can be used to quantify the relationship between complex ecological features, including habitat fragmentation [14], invasive species [15] and infectious [16,17]. From island-biogeography theory [18] to the current P-C-M model [19], it has been demonstrated that the closer the linkage between the green areas, the higher the similarity of species [20,21]. However, consideration of the difference in connectivity caused by wildlife moving across the substrate surface was considered more important than Euclidean distance, which has been theorized and demonstrated by several studies [22,23]. Therefore, the connectivity, mobility, etc., of most green networks have not utilized the proximity defined by Euclidean distance, but have utilized the proximity measured using advanced approaches such as least-cost distance.
Least-cost distance has been utilized to measure locality in the field of traffic geography, where various costs arise as the cargo travels and the optimal path needs to be identified [24]. Although the least-cost distance analysis was developed in the field of traffic geography, it was included in the initial development scope of the GIS program to address the fundamental problem with proximity. Since then, landscape ecologists have been able to use GISʹs Least-cost path analysis to demonstrate that the least-cost distance has advantages over the Euclidean distance in measuring ecological proximity [25,26]. Due to these advantages, many studies using the least-cost path are being conducted [27].
For example, Balbi et al. chose moths and birds as the target species to test the predictive efficiency of the last-cost path used to develop urban green networks [28]. Balbi et al. demonstrates the hypothesis that the mobility of the target species on last-cost paths is further increased [28]. In the context of urban areas, there is much discussion about target species when setting up the green network, but birds are considered suitable for target species as relatively mobile groups that can pass through the city [29,30]. In particular, as urbanization progresses and habitats are fragmented, bird species are reported to be homogenized [31][32][33].
The objectives of this work are: (a) distinguish important greenery in the construction of green networks; (b) conduct a quantitative analysis of current green networks using the least-cost distinction; (c) identify the relative importance of each green, path and corridor based on circuit theory; (d) additional greenery. By achieving these objectives, this study aims to present an analysis of the level of area that can be reflected in spatial planning, overcoming the existing green network limitations being prepared at the linear level. Through this, we believe that the results of this study will greatly contribute to increasing the possibility of green network reality.

Research Area
The study site is Suseong district (35°49′ N. 128°40′ E.), located east of Daegu Metropolitan City, South Korea ( Figure 1). The mountainous areas surrounding the southern part of the research areas are connected to Bisul Mountain to the west. The southern mountainous area consists of mountains with an average elevation of 600m, including Daeduk Mountain and Beopi Mountain. The northern region consists of low slopes and hills with an average elevation of 150 m. An alluvial plain formed by the flooding of Sin Stream is formed in the western region. In the eastern part of the country, a flat surface shaped like a fan is developed in the valley of Daeduk Mountain.
The western part of Daeguʹs downtown area and the eastern part of the country, which consists of residential areas, form an urbanized area and the hilly area between them is set as a development-restricted zone. However, over the past 30 years, urbanization has continued to expand around large-scale housing development projects and housing districts have been under development in areas adjacent to the development restricted areas. Over the past decade, the urbanized area has expanded by about 34.1% and, on the contrary, green areas continue to decline and, recently, pressure on the lifting of development-restricted zones has been increasing. In consideration of this serious situation, local governments have implemented several greenhouse policies. In 2017, a park-green plan was proposed with the aim of 2030. Daegu Metropolitan City designates the development-restricted zone of this research site as a core area for preserving the green network. In particular, the plan aims to restore the disconnected ecological axis and link the fragmented parks and forests within the city. This plan, when implemented, is expected to have a significant effect on realizing a sustainable city. Therefore, analyzing the green network to identify green areas that need connection and conservation is essential for preserving biodiversity and sustainable urban development.

Research Framework
In this study, we selected (1) viable target species in urban areas for analysis of green network and distinguished the core green areas where target species can inhabit. (2) A least-cost path analysis was conducted considering the characteristics of the substrate surface and the cost weighted distance between green areas. (3) By analyzing the centrality, pinch-points and restoration opportunities using circuit theory, the relative importance of each green area, path and corridor was identified and areas requiring additional green space were identified. The overall research process is presented in Figure 2 and the following sub-chapter details each step.

Selection of Target Species and Extraction of Core Green Areas
In the case of this study, due to the spatial specificity of urban build-up area and natural area appearing together, we selected the target species that can enhance citizensʹ aesthetic satisfaction along with the importance of preserving species. A total of 20 mammal species and 18 amphibian species were found to be emerging [34,35], but it was reported that there was little possibility of habitat and migration in urban areas, so it was limited to the target species for the establishment of green networks in urban areas [1,36]. However, birds can live in urban green areas, especially by creating a pleasant urban environment with their own appearance and sound, providing citizens with easy access to nature around their lives. In addition, birds are highly mobile and have active metabolism, which can immediately respond to environmental changes. Previous studies have shown 95 species of birds that have been confirmed to appear in the research area [37]. Among them, 12 species (Columba livia, Corvus corone, Dendrocopos leucotos, Eurystomus orientalis, Falco tinnunculus, Parus major, Passer montanus, Phasianus colchicus, Phoenicurus auroreus, Pica pica, Picus canus, Streptopelia orientalis) were selected as target species. This is because it is viable in cities and can be observed continuously by citizens throughout the year [38,39]. The target species mainly live in mountains near city, parks and they show the ecological characteristics of catching small seeds or insects in nearby farmland and open spaces [38][39][40][41]. These indicators are a good benchmark for species that effectively use the urban habitat and are motile at different spatial scales.
The core green areas where target species can be found were selected by reflecting the land use characteristics, minimum area standards and morphological characteristics. First, the selection of core green areas based on land use characteristics was based on the fact that small-sized grasshoppers lived evenly throughout Korea, but mainly in forests, rice paddies, fields and low hills. The minimum area standard for core green areas was set to 1 ha in this study by combining the prior study [42][43][44] on the habitat area of birds. In addition, in terms of morphology, the form of the core greenery was measured using the CIRCLE [45][46][47][48] among the landscape indices (Equation (1)).
Where ai is the area of patch i. ai s is the area of smallest circumscribing circle around patch i. CIRCLE's range is 0 to 1. CIRCLE is 0 for circular patches. In addition, the value approaches 1 for elongated, linear patches. FRAGSTATS 4.2 was used to measure CIRCLE values and, as a result, green areas with values below 0.8 were selected as core green areas.

Least-Cost Path Analysis
For the analysis of the least-cost path, quantified data are needed to estimate the resistance of the target species depending on the characteristics of various substrate surfaces [28,29,[49][50][51]. The cost to move the target species to habitat or an area rich in objects should be computerized. Here, the cost is affected by various landscape structures of a particular substrate plane due to natural and artificial factors [52]. The surface of the substrate between the core greenery may have factors that allow wild animals to move easily, or that hinder movement or change direction and the cost and resistance should be measured by considering the nature and arrangement of these factors.
The data used to prepare total resistance surfaces in this study are 1. biotope class, 2. NDVI (normalized difference vegetation index). First, the resistance values for each biotope class were set differently. In addition, the migration of the target species depends greatly on the vegetation density [29,[53][54][55]. Vegetation density can vary even with the same biotope class. There may also be vegetative and non-vegetable areas within the biotope of urban areas. Therefore, NDVI was utilized to determine the distribution and density of each vegetation. NDVI uses data from 27 July 2017, which show the least interference of clouds among Landsat8 satellite images [56], taken over the past five years, to set weights by classifying them into three classes (Table 1). Furthermore, since the least-cost path analysis is an analysis algorithm based on Raster data, it was necessary to convert resistance values in the form of vector data into Raster-based files [57]. Thus, in this work, we utilize the feature to raster function of ArcGIS 10.2 to convert the 30 m × 30 m cell size to a TIFF file to extract the resistance surface. The green network is identified by utilizing the analysis core greenery and resistance surface derived above. The analysis was conducted using Linkage Pathway to check the neighborhood between core green areas, analyze the least-cost path and analyze the connection path. First, the proximity between core green areas distinguished the substrate surface according to the nearest core green area on the distance and conducted Euclidean allocation, which only reflects the simple straight distance. If moving from one region to another must pass through another, it was determined that the two regions were not adjacent to each other. Leveraging these adjacency data, we exclude non-adjacent core green area pairs from the network configuration.
Next, the Linkage Pathway tool was used to analyze the cost weighted distance, the least-cost path and the connection path by reflecting the proximity data and resistance surface between the core green areas. The cost weighted distance represents the minimum value of the cumulative cost when the target species moves from the core greenery to the target point. Cost weighted distance maps written in this way provide estimates of relative access between each point and the closest core green areas and can be particularly useful in identifying areas that promote and inhibit wildlife migration [58]. The least-cost path analysis not only represents the most efficient path passing by the substrate plane between the base and the adjacent core greenery, but also analyzes cost weighted distance of the least-cost path by calculating the sum of all cost weighted distances on the path. The result is a relative comparison of each least-cost path, which enables the identification of areas with good connections and areas with poor connections. At this stage, the leastcost path from core green area A to B was automatically excluded if it passes through core green area C. The connection path analysis was previously analyzed by reflecting the cost weighted distance between the connected core green areas during the least-cost path analysis. The cumulative cost weighted distance value of the least-cost path and the cost weighted distance from a particular point to each core green area was derived to normalize the connecting path as shown below (Equation (2)).
where NLCCAB (Normalized Least-Cost Corridor) is the normalized connection between core green areas A and B, cost weighted distance (CWDA) is the cost weighted distance from core green areas A, CWDB is the cost weighted distance from core green areas B and LCDAB is the least-cost distance. Furthermore, the range of connecting path values starts with a value of zero, which means the optimal path. This process quantitatively identifies the reduced effectiveness of connectivity compared to the least-cost path and the derived route map can be useful for comparing the scope and effectiveness of the corridor area in establishing a green network. On the other hand, the existence of such a modeled connection path is not functionally connected but aims to present directions for further analysis and planning in the future.

Analysis of Green Network Using Circuit Theory
This study conducted a central, pinch-point and restore opportunity analysis based on McRaeʹs circuit theory [59] to identify areas of high importance, areas of high potential mobility, areas of limited development, areas of greatest restoration effect and core green areas of significantly affect population growth. The definitions and ecological interpretations of the basic concepts of resistance, conductivity, current and voltage used in the circuit theory are shown in Table 2 [59].

Electrical Term (Symbol, Unit) Ecological Interpretation
Resistance (R, ohm) the opposition that a resistor offers to the flow of electrical current.
Opposition of a habitat type to movement of organisms, similar to ecological concepts of landscape resistance or friction. Grid cells allowing less movement are assigned higher resistance.
Current (I, amp) flow of charge through a node or resistor in a circuit.
Current through nodes or resistors can be used to predict expected net movement probabilities for random walkers moving through corresponding core areas or paths. Voltage (V, volt) the potential difference in electrical charge between two nodes in an electrical circuit.
Voltages can be used to predict the probability that random walkers leaving any point on a green network will reach a given destination (representing, e.g., successful dispersal) before another.
Centrality Mapper [60] was used for centrality analysis, in which each core green area was treated as a node and each link was assigned a resistance corresponding to the cost weighted distance of the least-cost path. When a network of nodes and links is complete, the Centrality Mapper uses Circuitscape [61] to calculate current flow centrality. Current flow centrality iteratively works on all core green area pairs to inject 1 Amp of current into one core region and set the other core region to ground. The current values flowing across all nodes and links were then summed to calculate the centrality score for each. Centrality was classified as Class 4 by the Quantile classification method, of which the top 10% of the upper classifications with the highest current flow-mediated centricity were highlighted and plotted.
A pinch-point is a point where mobility is concentrated within a corridor when the target species moves from habitat to substrate surface to other habitats. Generally, these points appear in the shape of a funnel and due to the highly concentrated flow of species, the impact of disturbance on overall connectivity is very high [59], i.e., difficult to replace the path, so even a small area loss can eventually pose a huge risk to the entire green network. To identify these pinch-points, Pinch-point Mapper [62] was utilized in this work. Pinch-point Mapper has the advantage of using Circuitscape to schematize current flow over core green area pairs connected directly to the least-cost path. An electric current of 1 Amp was injected from one core green area and the other was set to ground. The current then flows through the gangway area, setting the point where the flow of charges within the gangway area was concentrated as a pinch-point.
The Barrier Mapper [63] was used for restoration opportunity analysis. Barrier Mapper utilizes the results of the previous analysis of the green network using cost weighted distances. Restoration opportunities are derived by Equations (3) and (4) based on the assumption that the detection area has been restored to a type with a minimum resistance value.
LCD' is the least-cost weighted distance through the center of the coverage area after removal and restoration of the barrier and CWDXMIN represents the least-cost weighted distance from the core green area X within the coverage area. D means the size of a particular area or diameter (50m diameter in this study) and R' means the adjusted resistance value (for this study, 1). The difference between the new least-cost weighted distance and the existing least-cost distance is defined as the opportunity to restore that pixel.

Core Green Area Extraction
A total of 43 greenery were found as a result of core greenery extraction ( Figure 3). The total area of 43 core green areas is about 36.18 km 2 , which is about 47.28% of the total area of the research site and 93.5% of the total green areas that are emerging within the research site. The core green area, which covers the largest area, is the No.9 green area, which covers an area of about 2156 ha. Most of the core green areas with large areas are set up as development-restricted zones and green areas 6, 7, 8, 30 and 31 adjacent to the urbanized areas showed relatively narrow areas. This is believed to have caused the phenomenon of the fragmentation of green areas around the city center due to the disorderly expansion of the city. In addition, even large-scale green areas are cut due to inter-fat connection roads, which are expected to become more severe in the future. The results are like Formanʹs claim that Koreaʹs development-restricted zones have recently been fragmented and cut into several large green spaces [1]. On the contrary, in the case of green areas No. 9 and 32, roads are passing by, but the area of green areas is maintained as it is constructed in the form of tunnels. As a result, the construction of roads was easing the fragmentation of green areas.

Resistance Value and Resistance Surface Calculation
The following Figure 4 shows the total resistance values of the urban ecological status map and NDVI combined. The highest resistance value was set at 300 and was found to be mainly concentrated in overcrowded residential and commercial areas. Areas with relatively low resistance values were centered on development-restricted zones. On the other hand, green-rich educational facilities and parks were significantly lower resistance values than those in the surrounding areas. However, the isolation level of these areas is expected to be relatively high because some neighborhood parks in the city center, such as No. 37 and No. 41, which are set as core green areas, are showing a wide range of resistance values.
The lowest resistance value in the core green area was 2.26, which was analyzed as No. 12. The area is considered to have a high density of vegetation due to the relatively large area of Quercus acutissima and Robinia pseudoacacia in the southeast, although the construction of roads in the northwest is undergoing slope recording. Conversely, the number 14 core green area had the highest average resistance value of 77.65 because it was used as a transportation island with a low density of vegetation on the highway. Overall resistance value map (Red has a high resistance value and blue has a low resistance value. The minimum resistance is 1 and the maximum resistance is 300.).

Result of Proximity Allocation Analysis between Core Green Areas
A total of 117 pairs of neighbors were found to be adjacent, of which 109 were adjacent according to a simple straight line using Euclidean allocation and 107 were found to be adjacent using cost allocation reflecting the resistance surface (Table A1 in  Appendix A).
For example, in the case of 11 cases, the simple direct election management was applied and found to be adjacent to a total of four regions, including 9, 13, 16 and 23. On the other hand, when applying cost weighted distances, it was analyzed that no. 16 and no. 23 of them were adjacent and also, areas 28 were found to be adjacent to a total of three regions ( Figure 5). Using these results, it is believed that the overall connectivity of urban green networks can be efficiently enhanced and the effectiveness of the plan can be enhanced by preferentially securing the connectivity of pairs adjacent to but not to the cost.

Least-Cost Distance, Least-Cost Path and Corridor Analysis Results
The results of cost distances and least-cost paths using the Linkage pathway are shown in Figure 6. The cost weighted distance is increasing as the market is getting bigger, especially in low-rise residential areas in the west. This area was mostly rice paddies and fields in the 1960s, but it seems that the value of the substrate surface has decreased in terms of ecology due to the reckless development. On the other hand, we analyzed the least-cost paths for 117 core green area pairs previously derived from the analysis of proximity. In addition, a total of 69 least-cost paths were derived except for 48 routes passing other core green areas (Table A1). The least-cost path with the highest cost weighted distance is Route 112, which connects Area 39 and 41, with a total cost weighted distance of approximately 149.5 km and the actual length of the route is approximately 2.2 km. This is attributed to the absence of overcrowded residential areas and a wide range of road-enclosed planting spaces, with an average resistance of 68. In addition, most of the long-established routes with a cost weighted distance of more than 30 kilometers were penetrating the marketplace of the research site, while the leastcost path passing inside the development-restricted zone was low.
Most of the least-cost paths were also identified as utilizing natural streams and small green areas. For example (Figure 7), in the case of No. 57 least-cost path linking No. 15 to No. 27 core green, it was analyzed that it passed through the forests and rivers behind the baseball field. In addition, the 33rd least-cost path connecting No. 9 to No. 16 was passing through a small forest near the high school and the least-cost path connecting No. 1 to No. 32 was via Childrenʹs Park. These results are judged to indicate the value of small green and linear green areas that were not previously selected in the target green area as connecting the core green area to the core green area. The least-cost corridor analysis results for verifying the extent and effectiveness of the gangway area are shown in Figure 8 below. A wide corridor with a lower cost than the least-cost corridor was formed around the wetlands and development-restricted zones around the Kumho River and a least-cost corridor was also found around the arable land near the No. 28 core green area. On the other hand, the corridor area across the street was quite narrow. Considering that the least-cost path and corridors are formed around the extrusion of the core green area, it is believed that extrusions of green areas are effective in attracting target species through the city to green areas. This shows significant similarities to the fence dispersion and fence effects of patches claimed by Forman [64]. In addition to the least-cost path, corridor areas have been identified and these areas are particularly likely to be used as alternative routes. For example, the least-cost path between core green areas No. 16 and No. 23 was set to pass through high schools and elementary schools, but the least-cost corridor were derived using apartment landscaping green areas, small rivers and roadside buffer green areas (Figure 9). Areas where this alternative corridor can be established are considered extremely valuable in terms of preserving connectivity.

Centrality
The centrality analysis was conducted to distinguish the relatively high-importance areas and routes for increasing and preserving the population of the target species and strengthening the connectivity of the green network in the target green area and least-cost path. The results are shown in Figure 10 below. The highest centrality of the target green area was analyzed at 32. The results are believed to be since the No. 32 green areas are located in the center of the research site and 15 green areas, which account for about 35% of the total green area and the least-cost path. In addition, the area is in the middle of the development-restricted zone and is believed to have the highest centrality as the least-cost path connected to the surrounding green area shows a relatively low-cost distance compared to other routes. Next, the No.9 target green area, which is analyzed for its high centrality, is in the south of the research site with a large area and is connected to 11 target green areas and the least-cost path, which is shown to have high values. However, compared to Area 32, it is believed that the connected routes passed through the marketplace, resulting in several high-cost routes being connected, resulting in relatively low values. The development of areas with high centralities, such as areas 9 and 32, is feared to isolate the target species to sub-object groups, so it needs to be set as a conservation area. In addition, large areas such as 17, 27, 29, 35 and 39 with relatively higher centrality values than other areas are highly useful as core areas for establishing green networks for research sites in conjunction with Area 32. However, even areas with low or moderate centrality may perform important functions to maintain a sub-object group of target species. For example, Beomeo Park No. 34 had a moderate centrality analysis across the study site but considering that it was the only area connecting Yashigol Park to Hwarang Park No. 41, this area could be considered a critical area for maintaining the target species of the surrounding green area. These centrality analysis results can be the basis for establishing the hierarchy of axes in the future planning of green networks.

Pinch-Point
A pinch-point analysis was required to determine where the movement of target species within the gangway region is concentrated. The results of the pinch-point analysis are shown in Figure 11. Overall, the frequency of pinch-points within the marketplace was high and the result is that the high resistance of the marketplace narrowed the width of the aisle area and the possibility of creating alternative routes was also reduced. First, green areas 11, 13 and 14 were analyzed to be vulnerable to isolation as the width of the corridor connecting to the adjacent green areas was narrow and the movement of species within the corridor was concentrated at certain points. Looking at the flow leading to 9-11-13-28, we can see pinch-point in the small-scale forest and Wooksu stream behind the middle ground and the pinch-pint was mitigated (Figure 12a) by widening the gangway width outside the market. In addition, if you look at the flow from area 29 to the Kumho River, which is located north of the target site, you can see that the surrounding substrate surface is dominated by agricultural land, but the pinch-point occurs along Maeho Stream (Figure 3-13b). In addition, to these reservoirs, the flow was found to be concentrated in the water purification plants near No. 5, 6, 10 and 12. It tells us that it is necessary to consider not only green but also water space together when planning a green network. If you look at the connections between 32 and 34, they are connected by corridor ①32-34, ②32-22-21-34, ③32-30-34 and the leastcost path is formed in corridor ① and ③ (Figures 3-13c). However, the pinch-point is high in the case of the two corridor. Therefore, it is believed that increasing the amount of green near the corridor of ② will alleviate the pinch-point of corridor of ① and ③.  Figure 11).

Restoration
Restoration opportunity analysis was conducted to identify areas that reduce the connectivity of green networks. The analysis results are shown in Figure 13 below. The area with the highest opportunity for restoration is the areas between No. 16 and No. 23. If green areas were to be increased, the least-cost distance could be reduced by up to 255 m. It was elucidated that if the street environment is improved as well as residential and commercial areas along Dal-guBeol road. This is expected to further strengthen the connectivity of green networks between the roadside park and Maeho Park, which lead to Daegu Agricultural Meister High School, DalguBeol road Street Green Area, Residential Area, Commercial Area and Shiji Elementary School. In addition, if you look at the areas between No. 34 and No. 37, you can see that improvement of public facilities such as Suseonggu Office, Suseong Police Station and Suseong District Beomeo Library is needed first. In the case of public facilities, there is an advantage that the feasibility of policies is relatively high compared to private lands such as houses and shops. In addition, the public notice in the land of the Suseong District Beomeo Library is expected to be highly needed not only in terms of ecology but also in terms of providing a pleasant resting place for citizens through improvement. In addition, the strengthening of the green network in this area is expected to significantly reduce the isolation of green areas No. 37 and 41, which have a single corridor.

Discussion
Existing prior studies mainly used graph theory [65][66][67] and gravity model [9,68,69] to approach the linear level in evaluating the connectivity and network of urban green areas. However, the analysis results of the linear level were limited in the application of spatial planning and an approach attempt to the area level with clear boundaries is required to realize the green network plan [70]. To address this issue, this study aims to effectively use urban habitats, select motile species at various spatial scales and analyze urban green networks using circuit theory and least-cost distance methodology.
Studies have been conducted in many ways to find restoration areas, especially in cities, to develop ecological integrity [71,72]. Kong et al. points out that setting up the restoration area under the supervision of the planners is ineffective. It claims that it is necessary to establish a restoration area by quantified evaluation methods [9]. In this study, restoration opportunities were presented divided into five grades. In the research area, it was analyzed that increasing the amount of green areas on the roadside and the built-up area and public facilities between No. 16 and No. 23 would be the most effective. These results can be used to make budget decisions to establish green network. There is also a similarity to [73] in that the establishment of the restoration area should consider the proximity to the minimum cost path and the potential for development as alternative paths. On the other hand, a number of studies have been conducted to establish conservation areas of high biodiversity value for sustainable cities [74][75][76]. In order to establish these conservation areas, the method and results of the centrality analysis in this study can be referred to. According to the centrality analysis, green area No. 32 with the highest value could be selected as a conservation area first. In the case of No. 32, it can be interpreted that the target species in the outer part of the city serve as gatekeepers to enter the research area. We believe that if the field survey of the target species is conducted in the future, it will be possible to hypothesize that No. 32 has a higher variety of target species than other green areas. Research also showed that most of the least-cost paths and corridors were concentrated in development-restricted zones. This is in line with Korean studies focusing on the ecological value of the restricted development zone [77,78]. It is believed that it will provide a basis for securing the ecological value of the city in Korea, where the pressure to lift the restricted development zones has been increasing recently.
We believe that applying the quantitative methodology used in this study to a 10year green network plan would be effective in determining the achievement rate of the plan. In addition, due to the characteristics of the analysis method, detailed variables such as habitat selection and resistance surface that are used as input data may vary when target species are applied differently, such as mammals and amphibians [79,80]. Subsequent analysis has the advantage of being able to automate. In particular, it is judged that this research can be utilized not only for improving connectivity but also for planning for reducing connectivity. For example, wild boars, rabbits and mice carrying germs could be used for space planning to prevent the invasion of harmful tides or infectious diseases [81,82].
In future studies, the results of the study can be verified by conducting a field survey to identify the species distribution of the target species. In addition, research can be conducted to explore the connectivity with the green network plan of wide-area units by performing analysis on various spatial scales.

Conclusions
In this study, we analyzed the green network using the least-cost distance and circuit theory for Suseong district, Daegu, South Korea. A total of 12 species of birds that can effectively use urban habitats and apply to various spaces are set as target species. A total of 43 core green areas were selected and the resistance surface considered the biotope type and NDVI. A total of 69 least-cost paths were derived. The least-cost paths were mainly through streams and small green areas. The corridors through the built-up areas was generally narrow, with high bottlenecks. The area with the highest centrality in the research area was found to be No. 32 green, located in the center. The area with the highest chance of restoration was the built-up area between No. 16 and No. 23. Previously, green networks were analyzed for urban areas, but our research was analyzed at the area level to make it easier to apply to spatial planning than previous studies. If the target species changes, the results can be different. Future researchers should consider conducting field surveys and comparing them with the results of this study.