Exploring Spatial Distribution of Pollen Allergenic Risk Zones in Urban China

Junfeng Xu 1,2, Zhanqing Cai 1,2, Tiantian Wang 1,2, Guang Liu 1,2, Peng Tang 1,2 and Xinyue Ye 3,* 1 Institute of Remote Sensing and Earth Sciences, Hangzhou Normal University, Hangzhou 311121, China; xjf11@zju.edu.cn (J.X.); fortyjames@163.com (Z.C.); wangtt1130@163.com (T.W.); yidoua@163.com (G.L.); iviptp@163.com (P.T.) 2 Zhejiang Provincial Key Laboratory of Urban Wetland and Regional Change Studies, Hangzhou 311121, China 3 Department of Geography, Kent State University, Kent, OH 44240, USA * Correspondence: xye5@kent.edu; Tel.: +1-419-494-7825; Fax: +1-330-672-4034


Introduction
With the growing trend of urbanization, it is estimated that the majority of the world's population will live in urban areas.The growth of urban areas and the impact of urban ecosystems on public health and urban sustainability have been leading issues of both academic and public interest.Pollen is the male reproductive cell of spermatophytes.Pollen grains are not easily oxidized and are small in size, occurring in large quantities [1,2].The main reason of inducing the seasonal allergic rhinitis and asthma is the airborne allergic pollen, which is largely caused by the main inhalation allergens in outdoor air [3][4][5].The impact of tree species and associated urban landscape change have raised growing concerns about the public health risk issues involved in daily activities of individuals.Some people have allergic reactions after inhaling pollen, commonly referred to as a pollen allergy.Pollen allergies show significant seasonal and regional features.In the USA, 2%-10% residents suffer from pollen allergies, and most cases are associated with Ambrosia artemisiifolia (ragweed) [6,7].In Mexico, 23%-31% of allergic reactions are due to pollen allergy [8][9][10].In Japan, about 6 million people are affected by pollen, mainly from Cryptomeria japonica (L.f.) D. Don [11].The allergenic pollen in northern China is mainly from Ambrosia artemisiifolia L. which belongs to the wild Compositae weeds in Artemisia [12,13].Typical manifestations of pollen allergy include sneezing, runny nose, watery eyes, and even severe itchiness of the nose, eyes, and external auditory canal.In severe cases, pollen allergy can induce bronchitis, bronchial asthma, pulmonary heart disease, etc., thereby affecting the daily life, work, and travel of the affected population [14].Without timely treatment, pollen allergy can develop into chronic asthma, conjunctivitis, pneumonia, and other respiratory diseases, and in severe cases can cause heart failure, kidney damage, or other fatal diseases, or even shock and sudden cerebral death [14].Exact data on the overall morbidity of pollen allergies in China have yet to be reported, but [15] estimated 0.9% morbidity (up to 5% in endemic regions) from pollen allergies in China, with an increasing trend.
Traditional methods to detect airborne pollen needs to be conducted through the point of distribution of pollen collection device to get pollen concentration.So a large amount of pollen collection devices are needed to get a more complete regional surface pollen distribution information, which requires much time and energy [16][17][18].Remote sensing technology has incomparable ground monitoring capabilities, with high time-effectiveness, large areal coverage, and traceability.By the complementary of the two methods, we can not only get a large area of the distribution of the pollen, but also can know some points of the pollen concentration information.
In this study, by utilizing high-resolution satellite remote sensing data, an object-oriented classification method was used to recognize tree species.A field investigation was conducted to verify the classification results, and then information on tree species distribution was extracted for different seasons.Spatial analysis was conducted to extract risk zones for high-level pollen concentration and to analyze the results.

Overview of Study Area
Though urban areas only account for a small proportion of the earth, the population flow is highly dynamic, bringing in a large number of visitors.Such intensive and rapid human movement exacerbates impacts of pollen allergenic risk zones in Urban China, which has been experiencing a restless and dramatic growth both vertically and horizontally.The study area is located in the West Lake Scenic Area of Hangzhou, Zhejiang Province, China (120 • 8 E, 30 • 15 N).This scenic area is of global importance for cultural heritage-especially the plant landscape, receiving many tourists each year from all over the world.This metropolitan area has been the preferred destination for millions of domestic migrants and overseas investors.While studies have indicated that the pollen of these plants can be caused by different degrees of pollen allergy [19].According to the West Lake Hangzhou Scenic Area Management Committee statistics, the total number of visitors to all pay-sight spots of the West Lake in 2014 was about 15 million person-time.Considering that pay-sight spots attract one-quarter of all visitors to the West Lake Scenic Area, so the total number of visitors to the West Lake Scenic Area in 2014 was estimated at more than 60 million, of which a large overall number would be expected to be pollen sensitive [14].According to survey results, more than 40% of the total number of tourists visited sites such as Breeze-ruffled Lotus at Quyuan Garden, Su Causeway, Broken Bridge with Thawing Snow, New Lakeside, etc., therefore those representative areas were selected for this study.Zoning and classification of a study area can improve classification accuracy [20][21][22]; since locations have differing vegetation distributions, the study area was divided into three subareas: Breeze-ruffled Lotus at Quyuan Garden and Su Causeway; Solitary Hill and Bai Causeway; and New Lakeside , as shown in Figure 1.

Study Materials
WorldView-2 data for 31 December 2009 were used, including panchromatic data with 0.5 m resolution and multispectral data with 2 m resolution.The survey day was sunny and cloudless, and the remote sensing images were clear with good imaging quality.Field-measured GPS data for the study area were used in combination with auxiliary data from Tencent "Street View Map" images.The data were analyzed using ArcGIS 10.1 (Esri, Redlines, CA, USA) and ENVI 5.0 (Exelis VIS, Boulder, CO, USA) software.

Study Methods
The dramatic urban transition in coastal China has led to accelerated urban expansion and a booming tourism industry, which has driven many people to seek various opportunities and move around in urban areas.China's unprecedented and multifaceted urban growth has generated many fascinating issues for scholarly research.Understanding the interaction between human dynamics and urban environment in China would allow appropriate strategies and policies to be formulated to facilitate the sustainability and environmental justice studies in China.Identifying pollen allergenic risk zones is of interest in urban landscape planning and management.Remote sensing technology was used to identify tree species in the study area.The occurrence and seasonal flowering phases of these species were used to map zones that present pollen allergenic risk.First, high-resolution satellite remote sensing images were pre-processed to extract vegetation cover for the study area [23,24].Tree types are judged by the information of vegetation spectra and texture-because the different tree species have different canopy shape features, different leaf size and density, and the canopy images reflected in the remote sensing images have different texture features.The gray-level co-occurrence matrix (GLCM) method is usually used to extract texture information [25].Combined with the field survey data, the classification and distribution range of the species were determined.From this, tree species were identified.Subsequently, the tree species category, distribution range and tree species were verified by means of field survey.Next, the spatial

Study Materials
WorldView-2 data for 31 December 2009 were used, including panchromatic data with 0.5 m resolution and multispectral data with 2 m resolution.The survey day was sunny and cloudless, and the remote sensing images were clear with good imaging quality.Field-measured GPS data for the study area were used in combination with auxiliary data from Tencent "Street View Map" images.The data were analyzed using ArcGIS 10.1 (Esri, Redlines, CA, USA) and ENVI 5.0 (Exelis VIS, Boulder, CO, USA) software.

Study Methods
The dramatic urban transition in coastal China has led to accelerated urban expansion and a booming tourism industry, which has driven many people to seek various opportunities and move around in urban areas.China's unprecedented and multifaceted urban growth has generated many fascinating issues for scholarly research.Understanding the interaction between human dynamics and urban environment in China would allow appropriate strategies and policies to be formulated to facilitate the sustainability and environmental justice studies in China.Identifying pollen allergenic risk zones is of interest in urban landscape planning and management.Remote sensing technology was used to identify tree species in the study area.The occurrence and seasonal flowering phases of these species were used to map zones that present pollen allergenic risk.First, high-resolution satellite remote sensing images were pre-processed to extract vegetation cover for the study area [23,24].Tree types are judged by the information of vegetation spectra and texture-because the different tree species have different canopy shape features, different leaf size and density, and the canopy images reflected in the remote sensing images have different texture features.The gray-level co-occurrence matrix (GLCM) method is usually used to extract texture information [25].Combined with the field survey data, the classification and distribution range of the species were determined.From this, tree species were identified.Subsequently, the tree species category, distribution range and tree species were verified by means of field survey.Next, the spatial ranges of pollen dispersal were determined from the flowering phases of the various tree species and their patterns of pollen dispersal.Given that pollen spread uses the whole canopy as the dispersal source, the GIS buffer function was used for spatial analysis of the possible pollen dispersal area.Pollen allergenic risk zones of different tree species were extracted separately.Finally, the spatial and seasonal distributions of pollen allergenic risk zones within the study area were determined by integrating the distribution information of pollen allergenic risk zones for different tree species.The detailed process is shown in Figure 2.
ranges of pollen dispersal were determined from the flowering phases of the various tree species and their patterns of pollen dispersal.Given that pollen spread uses the whole canopy as the dispersal source, the GIS buffer function was used for spatial analysis of the possible pollen dispersal area.Pollen allergenic risk zones of different tree species were extracted separately.Finally, the spatial and seasonal distributions of pollen allergenic risk zones within the study area were determined by integrating the distribution information of pollen allergenic risk zones for different tree species.The detailed process is shown in Figure 2.

Processing of Remotely Sensed Images
Recognition of tree species from remotely sensed images requires highly-specific processing.The image processing method utilized WorldView-2's advantages, high resolution and abundant spectral information, to maximize the information extracted from the remote sensing images [26,27].Within a study area, it is usually easy to distinguish evergreen species and deciduous species from images obtained in fall or winter seasons [28,29].Our images were usually unprocessed or only underwent simple parameter correction processing, and were specifically pre-processed for different research requirements.The images mainly underwent atmospheric correction to faithfully reflect the spectral information of ground objects.Fusion processing was also required for images to determine the spectral features of different ground objects, and to reduce human error during training sample selection and improve classification accuracy.Although high-resolution remote sensing images were used, the phenomenon of "different spectra for the same object, same spectrum for different objects" [30,31] was still observed.This phenomenon was even more pronounced in this study, because it was necessary to differentiate between tree species of the same genus.Hence, pre-classification partition and stratification processing of the study area were still necessary to improve classification accuracy.

Processing of Remotely Sensed Images
Recognition of tree species from remotely sensed images requires highly-specific processing.The image processing method utilized WorldView-2's advantages, high resolution and abundant spectral information, to maximize the information extracted from the remote sensing images [26,27].Within a study area, it is usually easy to distinguish evergreen species and deciduous species from images obtained in fall or winter seasons [28,29].Our images were usually unprocessed or only underwent simple parameter correction processing, and were specifically pre-processed for different research requirements.The images mainly underwent atmospheric correction to faithfully reflect the spectral information of ground objects.Fusion processing was also required for images to determine the spectral features of different ground objects, and to reduce human error during training sample selection and improve classification accuracy.Although high-resolution remote sensing images were used, the phenomenon of "different spectra for the same object, same spectrum for different objects" [30,31] was still observed.This phenomenon was even more pronounced in this study, because it was necessary to differentiate between tree species of the same genus.Hence, pre-classification partition and stratification processing of the study area were still necessary to improve classification accuracy.

Atmospheric Correction
Atmospheric correction was performed to eliminate the influence of light factors on ground object reflection; the influence of atmospheric molecules and aerosol scattering; and the influences of atmospheric substances such as water vapor, carbon dioxide, methane, ozone, etc. on ground object reflection.In this study, the ENVI Quick Atmospheric Correction module was employed to eliminate the influence of atmosphere and light on ground object reflection in images, restoring the actual spectral information of ground objects, obtaining the real reflectivity of ground objects and thereby allowing for quantitative analyses such as inversion of ground object parameters [32].Comparison between ground object spectral information, before and after quick atmospheric correction, shows that the corrected spectrum curve better fits the actual situation (Figure 3).

Atmospheric Correction
Atmospheric correction was performed to eliminate the influence of light factors on ground object reflection; the influence of atmospheric molecules and aerosol scattering; and the influences of atmospheric substances such as water vapor, carbon dioxide, methane, ozone, etc. on ground object reflection.In this study, the ENVI Quick Atmospheric Correction module was employed to eliminate the influence of atmosphere and light on ground object reflection in images, restoring the actual spectral information of ground objects, obtaining the real reflectivity of ground objects and thereby allowing for quantitative analyses such as inversion of ground object parameters [32].
Comparison between ground object spectral information, before and after quick atmospheric correction, shows that the corrected spectrum curve better fits the actual situation (Figure 3).

Image Band Fusion
For image fusion, panchromatic WorldView-2 satellite images have high resolution and can reflect text information of ground objects, which is however presented in grayscale, and thereby provides poor color representation.Its multispectral spatial resolution is relatively lower than that of the panchromatic band, with less prominent detailed information.However, it includes abundant spectral information and can reflect color conditions of different types of ground objects.In this study, data fusion methods in ENVI software were used to combine the advantages of high spatial resolution of panchromatic data and the broad color range of multi-spectra, to generate new data with both high spatial resolution and prominent color features that facilitated analysis and interpretation of images and to some extent improved the application value of images.The findings indicate that the Gram-Schmidt spectral sharpening fusion method improves image quality.Figure 4 compares the images before and after fusion.

Image Band Fusion
For image fusion, panchromatic WorldView-2 satellite images have high resolution and can reflect text information of ground objects, which is however presented in grayscale, and thereby provides poor color representation.Its multispectral spatial resolution is relatively lower than that of the panchromatic band, with less prominent detailed information.However, it includes abundant spectral information and can reflect color conditions of different types of ground objects.In this study, data fusion methods in ENVI software were used to combine the advantages of high spatial resolution of panchromatic data and the broad color range of multi-spectra, to generate new data with both high spatial resolution and prominent color features that facilitated analysis and interpretation of images and to some extent improved the application value of images.The findings indicate that the Gram-Schmidt spectral sharpening fusion method improves image quality.Figure 4 compares the images before and after fusion.Atmospheric correction was performed to eliminate the influence of light factors on ground object reflection; the influence of atmospheric molecules and aerosol scattering; and the influences of atmospheric substances such as water vapor, carbon dioxide, methane, ozone, etc. on ground object reflection.In this study, the ENVI Quick Atmospheric Correction module was employed to eliminate the influence of atmosphere and light on ground object reflection in images, restoring the actual spectral information of ground objects, obtaining the real reflectivity of ground objects and thereby allowing for quantitative analyses such as inversion of ground object parameters [32].
Comparison between ground object spectral information, before and after quick atmospheric correction, shows that the corrected spectrum curve better fits the actual situation (Figure 3).

Image Band Fusion
For image fusion, panchromatic WorldView-2 satellite images have high resolution and can reflect text information of ground objects, which is however presented in grayscale, and thereby provides poor color representation.Its multispectral spatial resolution is relatively lower than that of the panchromatic band, with less prominent detailed information.However, it includes abundant spectral information and can reflect color conditions of different types of ground objects.In this study, data fusion methods in ENVI software were used to combine the advantages of high spatial resolution of panchromatic data and the broad color range of multi-spectra, to generate new data with both high spatial resolution and prominent color features that facilitated analysis and interpretation of images and to some extent improved the application value of images.The findings indicate that the Gram-Schmidt spectral sharpening fusion method improves image quality.Figure 4 compares the images before and after fusion.

Partition and Stratification
Regarding stratification of the study area, before classification and extraction of remote sensing images, it was considered that since only the vegetation's canopy information was extracted, masking non-vegetation information from classification using vegetation index could effectively reduce misclassification and improve classification accuracy.Grass lawns were also extracted from vegetation based on different texture information, leaving only trees.

Recognition of Tree Species
Object-oriented classification is appropriate for differentiating tree species.One of its advantages is to classify based on the texture, shape, and proximal information of different image objects of different types that share similar spectral features [33].Classification of remote sensing imagery also relies on the accuracy of interpretation marks and the selection of a training sample.In this study, ENVI software was employed and object-oriented classification was used to perform ground object classification of images.Additionally, the classification results, class vector data, and classification statistical data were exported to facilitate data analysis and post-classification processing.
Visual interpretation markers were mainly determined from field GPS survey and verification, and through auxiliary recognition of Tencent "Street View Map" software (Tencent, Shenzhen, China).The interpretation markers are shown in Table 1.

Tree Species Interpretation Marker Tree Species Interpretation Marker
Platanus acerifolia Willd.Through field survey, we identified dominant tree species in the Breeze-ruffled Lotus at Quyuan Garden and Su Causeway subarea: camphor, metasequoia, salix, osmanthus, sapindus, platanus acerifolia, and quercus acutissima.Therefore, only those species were classified and extracted.The classification results are shown in Figure 5 and classification accuracy is shown in Table 2. User accuracy means the classification of the various land cover categories, corresponding to the ground true reference, the real number of pixels of the category of the percentage.Users do not know the true value graph, and therefore can only be based on the results given by the classifier to guess how much is reliable.So this calculation is called user accuracy.Producer means true ground reference data of the land cover classification of the test point, there are parts of the test point is the wrong classification, and the percentage of the number of pixels is correctly classified.Producers know the true value of the map, they use the evaluation of the classifier as the true value of the number of points in the graph of the number of points.

Recognition of Tree Species
4.1.1.Breeze-Ruffled Lotus at Quyuan Garden and Su Causeway Subarea Through field survey, we identified dominant tree species in the Breeze-ruffled Lotus at Quyuan Garden and Su Causeway subarea: camphor, metasequoia, salix, osmanthus, sapindus, platanus acerifolia, and quercus acutissima.Therefore, only those species were classified and extracted.The classification results are shown in Figure 5 and classification accuracy is shown in Table 2. User accuracy means the classification of the various land cover categories, corresponding to the ground true reference, the real number of pixels of the category of the percentage.Users do not know the true value graph, and therefore can only be based on the results given by the classifier to guess how much is reliable.So this calculation is called user accuracy.Producer means true ground reference data of the land cover classification of the test point, there are parts of the test point is the wrong classification, and the percentage of the number of pixels correctly Producers know the true value of the map, they use the evaluation of the classifier as the true value of the number of points in the graph of the number of points.It can be inferred from the classification results and estimation that the tree species in the subarea are mostly metasequoia, camphor, and platanus acerifolia, as the distribution of these species is relatively concentrated.Sapindus and salix have the highest producer's accuracy and the lowest user's accuracy, which may be due to their low proportions in tree species distribution.Other tree species have very high recognition accuracies, resulting in an overall classification accuracy of 86%.

Solitary Hill and Bai Causeway Subarea
Solitary Hill is a small peninsula on the north bank of the West Lake.In terms of plant allocation, the road along Bai Causeway is densely planted with salix, while different peach trees are planted extensively along the causeway bank.Tall platanus acerifolias are planted along the roads on Solitary Hill, with the dominant tree species being camphor, plum, quercus acutissima, liquidambar, sapindus, metasequoia, etc.Only these dominant tree species were classified and extracted.The classification results are shown in Figure 6 and classification accuracy is shown in Table 3.It can be inferred from the classification results and estimation that the tree species in the subarea are mostly metasequoia, camphor, and platanus acerifolia., as the distribution of these species is relatively concentrated.Sapindus and salix have the highest producer's accuracy and the lowest user's accuracy, which may be due to their low proportions in tree species distribution.Other tree species have very high recognition accuracies, resulting in an overall classification accuracy of 86%.

Solitary Hill and Bai Causeway Subarea
Solitary Hill is a small peninsula on the north bank of the West Lake.In terms of plant allocation, the road along Bai Causeway is densely planted with salix, while different peach trees are planted extensively along the causeway bank.Tall platanus acerifolias are planted along the roads on Solitary Hill, with the dominant tree species being camphor, plum, quercus acutissima, liquidambar, sapindus, metasequoia, etc.Only these dominant tree species were classified and extracted.The classification results are shown in Figure 6 and classification accuracy is shown in Table 3.According to the classification results, the Solitary Hill and Bai Causeway subarea is dominated by camphor, platanus acerifolia, and quercus acutissima.Salix is mainly distributed on Bai Causeway and along the lake shores.platanus acerifolia is mainly distributed along both sides of the main road.Metasequoia, platanus acerifolia, and quercus acutissima have the highest producer's accuracy, while plum, metasequoia, and sapindus have the highest user's accuracy.The overall accuracy is 90.36%.Liquidambar has the lowest classification accuracy.

New Lakeside Subarea
New Lakeside subarea is located on the east side of the West Lake, stretching to Beishan Street in the north and Wenying Pavilion in the south.The dominant tree species in this subarea are: camphor, platanus acerifolia, sapindus, salix, cedrus, magnolia grandiflora, and quercus acutissima, as shown in Figure 7.According to the classification results, the Solitary Hill and Bai Causeway subarea is dominated by camphor, platanus acerifolia, and quercus acutissima.Salix is mainly distributed on Bai Causeway and along the lake shores.platanus acerifolia is mainly distributed along both sides of the main road.Metasequoia, platanus acerifolia, and quercus acutissima have the highest producer's accuracy, while plum, metasequoia, and sapindus have the highest user's accuracy.The overall accuracy is 90.36%.Liquidambar has the lowest classification accuracy.

New Lakeside Subarea
New Lakeside subarea is located on the east side of the West Lake, stretching to Beishan Street in the north and Wenying Pavilion in the south.The dominant tree species in this subarea are: camphor, platanus acerifolia, sapindus, salix, cedrus, magnolia grandiflora, and quercus acutissima, as shown in Figure 7.According to the classification results, the New Lakeside subarea has a higher distribution of camphor and platanus acerifolia; the distributions of other tree species are comparatively similar.According to the classification results, the New Lakeside subarea has a higher distribution of camphor and platanus acerifolia; the distributions of other tree species are comparatively similar.Cedrus and sapindus have the highest producer's accuracy and the lowest user's accuracy.Quercus acutissima has the highest user's accuracy.The overall classification accuracy is 86.05%, as shown in Table 4.In this study, the testing sample set comprised one-on-one GPS points (one point corresponded to one tree) obtained through field survey, and supplementary points obtained from Tencent Street View mapping.ENVI software was used to estimate the accuracy of the classification results via a confusion matrix.The overall accuracy and Kappa coefficient can be calculated from the confusion matrix; classification accuracy is more than 85%, which is acceptable.
There are two peak pollen periods in eastern China, in the spring (March to May) and fall (August to October) seasons [34].In the present survey, the tree species flowering in the spring are camphor, quercus acutissima, salix, platanus acerifolia, liquidambar, plum, hovenia, celtis sinensis, and white oak; whereas osmanthus fragrans and cedrus deodara flower in the fall.The flowering phases of metasequoia, sapindus, and magnolia grandiflora do not occur during peak pollen periods, and therefore the distributions of their pollen allergenic zones were not extracted.The combined distributions of flowering tree species in spring and fall are shown in Figures 8 and 9, respectively.In this study, the testing sample set comprised one-on-one GPS points (one point corresponded to one tree) obtained through field survey, and supplementary points obtained from Tencent Street View mapping.ENVI software was used to estimate the accuracy of the classification results via a confusion matrix.The overall accuracy and Kappa coefficient can be calculated from the confusion matrix; classification accuracy is more than 85%, which is acceptable.
There are two peak pollen periods in eastern China, in the spring (March to May) and fall (August to October) seasons [34].In the present survey, the tree species flowering in the spring are camphor, quercus acutissima, salix, platanus acerifolia, liquidambar, plum, hovenia, celtis sinensis, and white oak; whereas osmanthus fragrans and cedrus deodara flower in the fall.The flowering phases of metasequoia, sapindus, and magnolia grandiflora do not occur during peak pollen periods, and therefore the distributions of their pollen allergenic zones were not extracted.The combined distributions of flowering tree species in spring and fall are shown in Figures 8 and 9, respectively.

Distribution of Pollen Allergenic Zones
Yang et al. pointed out that because pollen is dispersed via air currents, pollen concentrations differ according to the distance from the source plants.Excluding the influence of wind, the concentration of drifted pollen increases with distance, and reaches a maximum value in the area 10-15 m from the plant [35].Hence, in this study, 15 m was selected as the pollen allergenic radius of the risk zones.
The GIS software was used to vectorize the spatial distribution of tree species during their flowering phases, and the buffer function was used for spatial analysis, to establish a buffer zone of 15 m according to the vectorized point or line features of the tree species distribution.These buffer zones were considered as pollen allergenic zones.The spatial distributions of pollen allergenic zones in spring and fall are shown in Figures 10 and 11, respectively.It can be seen that in high pollen allergy incidence seasons, the spring flowering period, the entire West Lake Scenic Area forms a single, large allergenic zone that is larger during the spring than during fall.For this reason, visitors with a pollen allergy should avoid visiting the West Lake during these two periods.In the fall, tourists who suffer from a pollen allergy might be better advised to visit the Solitary Hill and Bai Causeway subareas.

Distribution of Pollen Allergenic Zones
Yang et al. pointed out that because pollen is dispersed via air currents, pollen concentrations differ according to the distance from the source plants.Excluding the influence of wind, the concentration of drifted pollen increases with distance, and reaches a maximum value in the area 10-15 m from the plant [35].Hence, in this study, 15 m was selected as the pollen allergenic radius of the risk zones.
The GIS software was used to vectorize the spatial distribution of tree species during their flowering phases, and the buffer function was used for spatial analysis, to establish a buffer zone of 15 m according to the vectorized point or line features of the tree species distribution.These buffer zones were considered as pollen allergenic zones.The spatial distributions of pollen allergenic zones in spring and fall are shown in Figures 10 and 11, respectively.It can be seen that in high pollen allergy incidence seasons, the spring flowering period, the entire West Lake Scenic Area forms a single, large allergenic zone that is larger during the spring than during fall.For this reason, visitors with a pollen allergy should avoid visiting the West Lake during these two periods.In the fall, tourists who suffer from a pollen allergy might be better advised to visit the Solitary Hill and Bai Causeway subareas.

Distribution of Pollen Allergenic Zones
Yang et al. pointed out that because pollen is dispersed via air currents, pollen concentrations differ according to the distance from the source plants.Excluding the influence of wind, the concentration of drifted pollen increases with distance, and reaches a maximum value in the area 10-15 m from the plant [35].Hence, in this study, 15 m was selected as the pollen allergenic radius of the risk zones.
The GIS software was used to vectorize the spatial distribution of tree species during their flowering phases, and the buffer function was used for spatial analysis, to establish a buffer zone of 15 m according to the vectorized point or line features of the tree species distribution.These buffer zones were considered as pollen allergenic zones.The spatial distributions of pollen allergenic zones in spring and fall are shown in Figures 10 and 11, respectively.It can be seen that in high pollen allergy incidence seasons, the spring flowering period, the entire West Lake Scenic Area forms a single, large allergenic zone that is larger during the spring than during fall.For this reason, visitors with a pollen allergy should avoid visiting the West Lake during these two periods.In the fall, tourists who suffer from a pollen allergy might be better advised to visit the Solitary Hill and Bai Causeway subareas.

Discussion and Conclusions
Most previous studies emphasized the vital benefits from ecosystem functions that underlie the foundation of human society in urban areas.However, the impact of tree species on public health should also be paid attention in order to facilitate urban sustainability [36][37][38].As the world's largest developing country with about three decades of double-digit economic growth, China's systemic shift from state socialism to a market economy has been the increasing human-environmental conflicts.Given the important role of metropolitan areas in China's economic and social transition, there is a growing need to address emerging environmental challenges and develop policies for sustainable human settlement development.The spatial relationship needs to be explored between health risk and the tree environment.However, it is still a great challenge to incorporate a data-driven approach into urban landscape planning, although the concepts of pollen allergenic risk have provided new perspectives for a better understanding of the roots of human-environment conflict.In this study, an object-oriented classification system was used to identify tree species in the upper canopy, from high-resolution remote sensing images; fourteen species, including camphor and metasequoia, were identified.In addition, the distribution of tree species that flower during the peak pollen periods of spring and fall were separately extracted.Based on pollen dispersal patterns, a GIS buffer function was used for spatial analysis to determine the distribution ranges of pollen allergenic zones and to map pollen allergenic risk zones during the spring and fall.
The findings are summarized as follows: (1) Using remote sensing technology, the overall accuracy of tree species recognition for the study area exceeds 86%.(2) Most plants in the study area flower during the spring, thus the pollen allergenic risk zones are also extensive during spring and sparse in fall.(3) As a complementary approach to traditional methods, the use of remote sensing technology is both time-and energy-efficient, and can provide highly effective, broad-range monitoring of the study area.
The accelerated growth of large-scale urban environmental data has attracted considerable attention from scholars, planners, and policy makers.This study proposes a tree species recognition process based on high-resolution remote sensing images.The results were verified by field GPS data and auxiliary Tencent Street View mapping data; spatial analysis of tree species image classification

Discussion and Conclusions
Most previous studies emphasized the vital benefits from ecosystem functions that underlie the foundation of human society in urban areas.However, the impact of tree species on public health should also be paid attention in order to facilitate urban sustainability [36][37][38].As the world's largest developing country with about three decades of double-digit economic growth, China's systemic shift from state socialism to a market economy has been the increasing human-environmental conflicts.Given the important role of metropolitan areas in China's economic and social transition, there is a growing need to address emerging environmental challenges and develop policies for sustainable human settlement development.The spatial relationship needs to be explored between health risk and the tree environment.However, it is still a great challenge to incorporate a data-driven approach into urban landscape planning, although the concepts of pollen allergenic risk have provided new perspectives for a better understanding of the roots of human-environment conflict.In this study, an object-oriented classification system was used to identify tree species in the upper canopy, from high-resolution remote sensing images; fourteen species, including camphor and metasequoia, were identified.In addition, the distribution of tree species that flower during the peak pollen periods of spring and fall were separately extracted.Based on pollen dispersal patterns, a GIS buffer function was used for spatial analysis to determine the distribution ranges of pollen allergenic zones and to map pollen allergenic risk zones during the spring and fall.
The findings are summarized as follows: (1) Using remote sensing technology, the overall accuracy of tree species recognition for the study area exceeds 86%.(2) Most plants in the study area flower during the spring, thus the pollen allergenic risk zones are also extensive during spring and sparse in fall.(3) As a complementary approach to traditional methods, the use of remote sensing technology is both time-and energy-efficient, and can provide highly effective, broad-range monitoring of the study area.
The accelerated growth of large-scale urban environmental data has attracted considerable attention from scholars, planners, and policy makers.This study proposes a tree species recognition process based on high-resolution remote sensing images.The results were verified by field GPS data and auxiliary Tencent Street View mapping data; spatial analysis of tree species image classification results was performed using GIS software.Methods for determining the pollen allergenic risk zones of different tree species in different seasons reflected specific conditions within the study area.This method is fast, simple, and practicable, and can indicate the macro-scale distributions of pollen allergenic risk zones.
Nevertheless, the present study has some limitations.The determination of pollen risk zones should consider multiple factors, such as weather, temperature, humidity, wind speed, topography, light intensity, etc. [39].This study presumed flat topography and windless weather conditions.Hence, future studies should aim to establish an analytical model that integrates landform, wind speed, temperature, and other factors.Grade specification of pollen allergenic risk zones should also be considered in further studies.In planning for urban sustainability, we need to move beyond examining the specific tree species and consider the entirety of the human-dominated urban environmental system, in order to seek a more systematic solution.

Figure 1 .
Figure 1.Overview of West Lake Scenic Area, Hangzhou, Zhejiang Province, China, showing study subareas.Figure legend: Breeze-ruffled Lotus at Quyuan Garden and Su Causeway; Solitary Hill and Bai Causeway; New Lakeside.

Figure 1 .
Figure 1.Overview of West Lake Scenic Area, Hangzhou, Zhejiang Province, China, showing study subareas.Figure legend: Breeze-ruffled Lotus at Quyuan Garden and Su Causeway; Solitary Hill and Bai Causeway; New Lakeside.

Figure 2 .
Figure 2. Flow diagram of the study process.

Figure 2 .
Figure 2. Flow diagram of the study process.

Figure 3 .
Figure 3. Spectral differences in ground objects before and after atmospheric correction.

Figure 4 .
Figure 4. Comparison of images before and after fusion.

Figure 3 .
Figure 3. Spectral differences in ground objects before and after atmospheric correction.

Figure 3 .
Figure 3. Spectral differences in ground objects before and after atmospheric correction.

Figure 4 .
Figure 4. Comparison of images before and after fusion.

Figure 4 .
Figure 4. Comparison of images before and after fusion.

Figure 10 .
Figure 10.Spatial distribution of pollen allergenic risk zone during spring (Figure legend: pollen allergenic zone in spring).

Figure 10 .
Figure 10.Spatial distribution of pollen allergenic risk zone during spring (Figure legend: pollen allergenic zone in spring).

Figure 10 .
Figure 10.Spatial distribution of pollen allergenic risk zone during spring (Figure legend: pollen allergenic zone in spring).

Figure 11 .
Figure 11.Spatial distribution of pollen allergenic risk zone during fall (Figure legend: pollen allergenic zone in fall).

Figure 11 .
Figure 11.Spatial distribution of pollen allergenic risk zone during fall (Figure legend: pollen allergenic zone in fall).

Table 2 .
Accuracy of tree species recognition: Breeze-ruffled Lotus at Quyuan Garden and Su Causeway subarea.

Table 2 .
Accuracy of tree species recognition: Breeze-ruffled Lotus at Quyuan Garden and Su Causeway subarea.

Table 3 .
Accuracy of tree species recognition: Solitary Hill and Bai Causeway subarea.

Table 3 .
Accuracy of tree species recognition: Solitary Hill and Bai Causeway subarea.

Table 4 .
Accuracy of tree species recognition: New Lakeside subarea.

Table 4 .Table 4 .
Cedrus and sapindus have the highest producer's accuracy and the lowest user's accuracy.Quercus acutissima has the highest user's accuracy.The overall classification accuracy is 86.05%, as shown in Accuracy of tree species recognition: New Lakeside subarea.