Geometrical Assessment of Sunlit and Shaded Area of Urban Trees Based on Aligned Orthographic Views

: To quantify the ecosystem services of trees in urban environments, it is necessary to assess received direct solar radiation of each tree. While the Sky View Factor (SVF) is suitable for assessing the total incoming short- and longwave radiation ﬂuxes, its information is limited to speciﬁc points in space. For a spatial analysis, it is necessary to sample the area for SVF. A new geometrical method, Area View Factor (AVF), for the calculation of sunlit areas is proposed. AVF is the ratio of the unhidden, projected surface of an object to the whole projected surface of an object in a complex environment. Hereby, a virtual, orthographic camera is oriented in accordance to the sun’s position in the 3D model domain. The method is implemented in the microscale model SkyHelios, utilizing efﬁcient rendering techniques to assess AVF of all urban trees in parallel. The method was applied to Rieselfeld in Freiburg , Germany. The assessed sunlit area is compared to the SVF at the top of each tree and solar altitude angle, revealing a strong relationship between sunlit areas to solar altitude angles. This study shows that AVF is an efﬁcient methodology to assess received direct radiation of urban trees. Based on AVF, it is possible to identify urban areas with shaded and sunlit trees, but it can also be applied to other objects in complex environments. Therefore, AVF is applicable for urban architecture or energetic research questions.


Introduction
The urban morphology describing the composition, layout, and form of a city has great impact on the radiation absorption of urban trees [1][2][3]. Trees in streets and parks can be shaded by urban obstacles (e.g., by buildings and other trees), so that the absorption of direct and diffuse shortwave radiation is limited [4][5][6].
Radiation fluxes (especially photosynthetically active radiation (PAR; 0.34-0.71 µm)), which are modified by the processes of absorption, emission, reflection, and transmission in the complex urban environment, have great impact on the ecosystem services of urban trees and therefore on the micro-and bioclimate in cities [4,[7][8][9].
Quantifying the impact of ecosystem services on the urban microclimate as well as the assessment of ecological and financial benefits [18,19] is required for improved and optimized urban planning, architecture, human thermal comfort and health.
In order to quantify the impact of trees on the urban microclimate, it is necessary to quantify radiation fluxes modified by the urban morphology. Interception of solar irradiance by urban trees depends on the position of the sun (solar azimuth and zenith angle) in relation to the observed model domain. The absorption may be limited due to attenuation by the urban morphology.
The urban morphology can be characterized by various geometrical indicators for layout and density [20,21]. These provide an explicit or parametric specification of the morphology. It can be explicitly specified by geometrical properties of buildings and trees (position, shape, and height) describing the volume and location of green and gray infrastructure in urban districts. The orientation of buildings and vegetation in relation to the sunpath, depict the exposure of building facades and leafs surfaces towards the sun. The Height/Width ratio (H/W-ratio) is the aspect ratio of streets determined by the height of surrounding buildings in relation to the width of streets indicate the openness of urban streets. Sky View Factors (SVF) in open spaces [22,23] is a dimensionless parameter describe the exposure of a point in space to the open sky vault.
To specify the urban morphology in a parametric way, the mean distance between buildings is an indicator for the horizontal, spatial density of an urban district. Frontal and planar area densities (e.g., complexity and site coverage) are indicators for the vertical and horizontal density of urban forms, respectively.
To determine direct and diffuse shortwave irradiance on a single point in space the local, spherical Sky View Factor (SVF) [22,24,25] is the most commonly used parameter. Spherical SVF is the fraction of the visible sky of the upper hemisphere and is dependent on the Point of View (POV). SVF can be assessed from rendered fish-eye imagery [22,24,26,27] or real photographs as well as from approaches based on Monte Carlo computational ray tracing [28][29][30].
These techniques have been integrated into numerous numerical Urban Climate Models (UCMs) on the micro-and local scale to assess the thermal effects of urban structures [31]. Discrete Urban Canopy models (UCs) represent the urban morphology of typical street canyons in two-dimensional tiles. The physical processes can be modeled in one (e.g., Town Energy Balance model (TEB) [32], Arizona State University Single-Layer Urban Canopy Model (ASLUM) [33]), or more layers (see, e.g., in [30,33,34]) within the canopy to depict the surface energy balance (SEB) in the complex urban environment. However, individual building structures (e.g., from CityGML) cannot be represented. Besides Urban Canopy models, there are models for numerical fluid simulations (CFDs) in urban areas (e.g., ENVI-met [35,36] and PALM-4U [37][38][39][40]). Turbulence-resolving models are highly accurate but at the same time computationally costly and require a high-performance computing infrastructure (HPC), which is usually not available to municipalities and are at the same time difficult to operate. The subsequently presented method was integrated into SkyHelios, a three-dimension microscale model, which is capable of continuous Sky View Factor estimation in urban areas. SkyHelios applies vectorized geospatial data to resolve the urban morphology as individual buildings and trees of the respective study area in a vectorized three-dimensional scenery [41,42].
In this article, POV is changed to the source of shortwave radiation (e.g., sun) in order to calculate sun-lit areas of urban trees on an object-oriented basis, which are exposed to the point source of the shortwave radiation. The sunlit area is hereby expressed as the ratio between sunlit projected area in relation to fully projected area of the tree canopy. The ratio between these projected areas is called Area View Factor (AVF). This information is available on an object-oriented basis and represents the exposure of the whole object towards the sun, in contrast to SVF which considers the complete upper hemisphere, including incoming longwave reflectance radiation of the surrounding surfaces for a single point in space. By calculating AVF on an object-based basis, it is then possible to perform calculations for the entire urban tree canopy without having to resolve individual leaves. Therefore, AVF can be utilized in the (Two-) big-leaf model for transpiration and photosynthesis calculations [8,9,43].
The calculation of AVF is implemented by aligned orthographic views, assuming the rays of the sun to be in parallel. These synthetically rendered views are analyzed using image processing techniques to map pixel values to individual trees, allowing for the calculation of AVF for each tree in the whole study area at once.
With the aid of this new method, first the sunlit area of each tree canopy in a study area can be calculated in order to subsequently describe the tree physiological processes driven by the absorption of direct solar radiation. Furthermore, it is also possible to calculate the sunlit area of any object in an urban environment in relation to seasonal and daily conditions. This could be, e.g., photovoltaic systems, but also building facades or horizontal surfaces in a street canyon. Based on the calculation of the Area View Factor for each tree canopy, the following research questions will be addressed: (1) How does the sun-lit area of a tree canopy change over the course of the day? (2) What is the influence of urban morphology on tree canopy shading? (3) How are AVF and SVF related to each other?

Lambert's Cosine Law
A way to describe the irradiance on a horizontal surface is given by Lambert's cosine law. The law states that irradiance to a horizontal surface varies in relation to the orientation of the surface towards a point source of shortwave radiation, usually the sun [9]: This relationship is illustrated in Figure 1, showing the irradiance (I 0 , W m −2 ) for given solar radiation input and solar altitude angle (β,°). Θ is of same magnitude as β for geometrical reason. In the complex urban environment it is not possible to directly relate the irradiance to a given surface, because the surface might be obscured by other objects.
Center of the scene β θ I 0 Figure 1. Sun rays are incident in parallel for given solar altitude angle (β,°), depending on daytime and day of year. The flux density of these parallel beams (orange) to an orthogonal surface is I 0 (W m −2 ). The flux density on a surface I 0 is in relation to by the Lambert's cosine law Equation (1) (modified after the work in [9], p. 15).

SkyHelios
SkyHelios [24,25,41] is an urban microscale model with the ability to estimate crucial parameters of the urban environment running on desktop machines. It is a steady-state, time-independent model capable of estimating several thermal indices [42]: Physiologically Equivalent Temperature (PET [44,45]), Universal Thermal Climate Index (UTCI [46]) and Perceived Temperature (PT [47]). Furthermore, SkyHelios can compute different environ-mental factors: sunshine duration, shading, wind speed, and direction for areas of interest (AOIs) as well as points of interest (POIs). Beside that, SkyHelios can estimate the aerodynamic roughness, based on different approaches [48]. The development of SkyHelios is focused on some key aspects, to keep the software efficient, fast, and free of charge: • Reduction of time complexity by heavy utilization of cheap graphics processing units (GPU) and parallelism on 64 bit machines. Strong compliance with open standards of the Open Geospatial Consortium (OGC) by supporting well standardized data formats [49].
By creating virtual, three-dimensional (3D) scenes from 3D urban geoinformation (e.g., CityGML LOD 0, 1 and 3D-Shapefiles), SkyHelios can be classified as a 3D city model, which bears great potential in general [50]. All 3D entities in the 3D scene are based on vectorial geodata, allowing for the specification of spatial resolution for all results on demand. The 3D scene is managed by the Object-Oriented Graphics Rendering Engine [51] (MOGRE), which is licensed under the MIT license. For more information see [52,53]. SkyHelios is written in C# programming language, targeting the .NET Framework 4.6.1. and is therefore only available for computers running Microsoft Windows.

Geometrical Assessment of Sunlit and Shaded Area
The idea of Lambert's cosine law (Equation (1)) was adapted for the use in 3D scenes of the microscale model SkyHelios. The sunlit area of an urban tree is calculated based on a rendered image from the point of view of the sun (POV). The rendered image is hereby taken by an orthographic camera [52,53].
An orthographic camera is capable to render projected images with the qualities of an equal area projection, so that rendered areas are independent of their distances to the camera. Additionally, all parallel light beams are focused at the camera lens at infinite distance ( Figure 2). The view of the orthographic camera is recorded on a plane, called render target, in 3D space. An orthographic projection is contrary to the classical perspective projection, there all light beams are focused at the camera lens on the focal point at a specific distance (compare to focal length). The perspective projection does not preserve the equal area property, compared to the orthographic projection. Apart from three-dimensional modeling in computer graphics, orthographic imagery is well known in remote sensing [54] as well as for orthographic projections in classical geography [55]. The orthographic camera as well as the corresponding render target are setup with the required properties as specified in Table 1. The position of the camera, and width and height of the render target are determined by the position of the sun, dimensions of the scene, as well as orientation of the scene towards the sun. The height of the camera as well as the properties for the render target can be geometrically assessed by the properties of the scene, especially by its bounding box. The bounding box of a scene is the enclosing cuboid, defined by the minimal and maximal vector of that scene. The required properties can be geometrically calculated according to Figure 3.
Center of the scene α β Figure 3. Calculation of orthographic camera properties: height of the render target (green), camera height (blue), based on the scene bounding box (red dashed) and sun altitude angle (β,°). The scene is drawn according a 40°height angle of the sun. Table 1. Comparison between the geometrical and technical properties of a perspective camera projection and an orthographic camera projection. The render target is positioned along the line given by the position of the sun and the center of the scene. The size of the render target is calculated, so that the complete scene is covered by a single rendered picture. Afterwards, the spatial resolution of the render target is, based on its size, multiplied by a scaling factor in order to increase the accuracy of the rendered picture. The size of the render target is given in real-world dimensions.

Object-Based Analysis of the Rendered Orthographic View
In order to get object-based information from the rendered view, the scene has to be modified to distinguish between different objects in the scene (buildings, tree crowns, tree trunks, etc.) for further analysis. All elements (buildings, tree trunks, and the ground surface) in the scene, except crowns of urban trees, are rendered in solid black. These elements with color RGB (0,0,0) in the RGB color space are ignored in further analysis. Tree crowns are colored, based on unique identifiers (id), which get also encoded in the RGB color space (Figure 4). The encoding is implemented in the following form: The steps of scene modification are illustrated in Figure 4 for the study area (Section 2.6) at 05:00 UTC. The orthographic view of the scene based on POV of the sun shows all buildings and trees in the scene in orthographic projection (Figure 4a). Parts of the trees are obscured by buildings and other trees and therefore not visible in this rendered image. It is not possible to distinguish between single trees and map the unobscured area to each corresponding tree. When the whole scene beside tree crowns is rendered in black, buildings and unobscured areas of trees can be separated (Figure 4b). In order to map the unobscured, sunlit areas to single trees, each tree is encoded in accordance to its unique identifier ( Figure 4c). Each tree has an unique color in the orthographic rendered image and can be distinguished in further analysis. The rendered view is analyzed on a pixel basis. Pixels with specific pixel values (e.g., color values) are counted and mapped to the original unique identifiers of the urban trees. The size of the pixels are hereby given in real world dimensions. The number of pixels with a specific color corresponds to the sunlit, unobstructed area of a single tree crown. The unobstructed area of the tree crown is the part of the tree which receives direct shortwave radiation.

Orthogonal Area View Factor
The sunlit unobstructed area of each single urban tree crown is related to the complete, projected area of the tree crown. The projected area is hereby calculated as the area of a circle. This method is named orthogonal Area View Factor (AVF). The AVF is the fraction of projected visible surface of an object in a complex environment (A visible, projected , m 2 ) to the projected surface of an object in unobstructed, free space (A complete, projected , m 2 ). The AVF is dependent on solar azimuth angle (ψ,°), solar altitude angle (β,°), and surrounding urban morphology: In contrast to SVF, AVF is object based and knowledge about the projected surface of the object (e.g., projected crown area) is required.

Dataset: Rieselfeld in Freiburg, Germany
The implemented module is tested for a study area in Freiburg, Germany. The urban district Rieselfeld is located at the western fringe of the city (located at 47.9991 N, 7.7921 E; see Figure 5). A detailed three-dimensional city model with level of detail one and two (LOD1 and LOD2), as well as an urban tree cadastre was provided: • The city model is provided in the 3D CityGML Format specified after [56]. Therefore, the data for the buildings cover footprint and heights of all buildings in the study area.

Results
Urban trees in Rieselfeld get shaded by buildings and by other trees. The sunlit areas are calculated based on the assumption for clear sky condition, without consideration of diffuse radiation. The only source of direct shortwave radiation in this study is the sun. Reflections of shortwave radiation are neglected. The shading is dependent on the sunpath, and therefore on the time of the day but also on the day of the year (DOY). At first, the mean AVF between a typical spring and summer day are compared (Section 3.1) and the mean AVF is linked to the solar altitude angle (Section 3.2). Subsequently, the spatial, seasonal and daily dependencies of AVF are analyzed throughout the study area (Section 3.3). Finally, the spatial results for AVF are compared with SVF (Section 3.4) at the top of each tree crown.

Spatially Averaged Analysis of the Mean Area View Factor
An averaged overview of AVF of all urban trees in the scene is given by Figure 6 on a hourly basis for the 14th of March (spring) and June (summer). The timeseries for AVF follow a standard bell-shaped curve with its maximum during midday at 11:00 UTC in accordance to maximum solar altitude angle. Throughout the day (from 04:00 to 19:00 UTC) urban trees are more sunlit in summer than in spring. The lighting in summer starts earlier and lasts longer compared to spring. The gradients during dawn and dusk are greater in spring than in summer. The area of sun-exposed surfaces in summer is of greater magnitude compared to spring ( Table 2).  The standard deviation for each hour varies during the day, with highest values during dawn and dusk in spring.
A detailed insight for the temporal and spatial variation of AVF throughout the district is given in Section 3.3.

Relationship between the Mean Area View Factor and the Solar Altitude Angle
The solar altitude angle for the rising and setting phase of the sunpath was related to the averaged results of the sunlit area to assert the models results under clear sky conditions (Figure 7). The relationship shows strong correlation according to the Pearson's correlation coefficient. There was positive correlation observed between the solar altitude angle and the sun-lit area in spring (14 March) and summer (14 June): correlation coefficient of r spring = 0.964 with significance level of p spring = 4.56 × 10 −9 in spring and correlation coefficient of r summer = 0.959, with significance level of p summer = 1.89 × 10 −9 in summer. The maximum of the ratio of sunlit to max. sunlit area in summer is of greater magnitude compared to spring. The lighting of the scene is independent of the air temperature and relative humidity, such that there is no difference in the two scenarios for lighting.
The sunpath for each date of the study area was calculated with SkyHelios and is given by Figure 8 in a fish-eye image of the upper hemisphere at the center of the study area. The SVF is obviously the same, but the sunpath between 14th of March and June are differing, due to different solar altitude and azimuth angles.

Seasonally, Daily and Spatially Analysis of Area View Factors
The spatial variation of the sunlit area in Rieselfeld is given by Figures 9 and 10 for spring and summer, respectively. The figures show the ratio of sunlit area to max. sun-lit area (AVF) of urban trees for times between 04:00 and 19:00 UTC for the 14th of March and the 14th of June. Solar altitude and solar azimuth are given for each time step.
The maximum in June is of greater magnitude compared to March. The lighting of the areas occurs earlier and lasts longer in June, compared to March. It can be observed, that the lighting of the scene starts in the east of the district with the outermost trees. During the day, when the solar altitude angle reaches its maximum value, most trees of the scene are fully sunlit in summer but not with the same magnitude compared to spring. Different effects are observable in the scene: • Trees in open spaces are more sunlit compared to others, which are in streets with a great H/W-ratio.

•
The location of the trees within a street canyon has a great effect on lighting. Trees at the sun exposed side of the canyon are more sunlit, compared to trees on the shaded side. • Street canyons with high density of planted trees lead to more shaded trees.
By comparison of Figures 9 and 10 it is also observable, that the variation of the sun-lit area during dawn and dusk is greater in spring compared to summer. Figure 9. Ratio between sunlit area and max. sunlit area of urban trees for the 14th of March under clear-sky conditions. It is observable that the lighting of the study area starts according to the solar azimuth angle at the east side and ends at the west side of the district. It is also conspicuous, that the lighting of the trees is dependent on the orientation of the streets and buildings for small solar altitude angles, but less dependent for greater solar altitude angles. Figure 10. Ratio between sunlit area and max. sunlit area of urban trees for the 14th of June under clear-sky conditions. It is observable that the effect of the solar altitude angle is decreased compared to March because the solar altitude angle is of greater magnitude. It is also conspicuous, that the lighting of the trees starts earlier and lasts longer, compared to March.

Comparison of the Mean Area View Factor to the Sky View Factor
Daily mean Area View Factor can be utilized to analyze differences within a street in terms of lighting. While the Sky View Factor in the cross section of a typical east-west (E-W) oriented street is mirror-symmetrical, AVF is lower for trees on the side facing away from the sun than on the side facing the sun (Figure 11). This is observable by trees on the south side of a street which are shown in purple (AVF < SVF) and trees on the north side of the same street in green (SVF ≤ AVF). Trees in a north-south (N-S)-oriented street generally have smaller mean AVF, as the highest irradiation into such a street occurs only in a narrow time interval, namely, exactly when the orientation of the street corresponds to the azimuth angle of the sun. In open areas with small H/W-ratio in relation to the building in the immediate vicinity, the deviation between SVF and AVF is very small. If the horizon is not obscured, then the projected area of a tree crown is completely sun-lit, so that AVF and SVF are both equal to 1. In direct comparison, AVF and SVF are linked. It can be observed that SVF is larger at the top of the crown of tall trees than at the top of small trees; the same holds for AVF ( Figure 12). The taller the trees, the smaller the influence of the urban structures in the neighborhood around the tree.

Discussion
The results for sunlit areas based on the calculation of AVFs are clearly dependent on the day of the year and time of the day by design (Section 3.1). The fluctuations of the standard deviation throughout the day suggest a dependency of AVF on solar elevation. The relationship of AVF to the solar altitude angle (Section 3.2) is not linear and differs between dawn (rising sun) and dusk (setting sun). This can be explained by the fact that the orthographic view is constructed and placed in accordance to the solar altitude angle and aligned with the solar azimuth angle. The spatial distribution of AVF as described in Section 3.3 show typical pattern for different street orientations, H/W-ratios, SVF factors, and areas with dense vegetation. Five different patterns can be distinguished and are further discussed for streets in a city of mid-latitude:

1.
Urban trees are shading each other in dense arrangements. The effectiveness of urban trees to provide adequate cooling can be restricted due to self-shading by the trees [57]. Self-shading can be caused by intersecting canopies of individual trees or adverse arrangement of tall to small trees.

2.
Urban trees on the north-side of an east-west oriented street are more sun-lit compared to trees on the south-side of the same street. This result comply with the results by [58], that microclimatic modifications of the urban atmosphere by trees are dependent on the location of the trees within streets ("side of street") and the time of the day. This dependency can be explained by the linkage of solar exposure and urban morphology. It is in good agreement with the studies based on modeled results by [59,60], which show different thermal sensation within street canyons due to varying mean radiant temperature in E-W-and N-S-oriented street canyons. It can be assumed that this only applies to cities of mid-latitude, and not for cities on the southern hemisphere, there orientations are swapped.

3.
Urban trees in streets with great H/W-ratio are more sun-lit compared to trees in streets with low H/W-ratio.

4.
Urban trees in streets, which are aligned orthogonal to the azimuth angle of the sun are more shaded, compared to urban trees in streets, which are oriented in parallel to the azimuth angle. It is obvious, that an observed street is fully sunlit if it is aligned to the solar azimuth angle. But this only holds for single points in time. Therefore, street trees in an E-W-oriented street are more sun-lit during the complete course of the day than those in a N-S-oriented street.

5.
Urban trees in open spaces absorb more radiation compared to trees in narrow street canyons. Obscuration of urban trees by the urban morphology has its greatest magnitude for low solar altitude angle (i.e., in the morning and in the evening as well as in spring and autumn) and becomes less dependent for high solar altitude angles (during midday as well as in summer) due to the dependency of solar irradiance on horizon heightening and solar altitude angle. These findings are comparable to the findings of an earlier work showing that the solar availability in the urban context is depending on solar altitude angles [61]. The comparison between solar altitude angle and mean sunlit area of urban trees shows strong correlation (Figure 7). Trees which are fully exposed to solar radiation provide stronger cooling effect based on transpiration, compared to trees in the shade of buildings [62]. This leads to the optimal cooling potential of urban trees, under the assumption of unlimited water availability. Therefore, trees in open spaces with high solar exposure should be prioritized for mitigation of heat stress [57].
This, in general, is in good agreement, that the spatial configuration of green structures has a huge impact to the availability of sunlight and to its cooling effect [63][64][65].
SVF is particularly suitable for the quantification of emitted long-wave radiation of surrounding surfaces, because the upper hemisphere encompasses all surfaces in the neighborhood towards POV. Short-wave radiation is also accessible by SVF, but heterogeneous shading of surfaces is difficult to assess by SVF since the information of SVF is limited to single points in space. On the other hand, AVF is suitable for the short-wave radiation towards a complete surface, including heterogeneous shading, because the complete surface is evaluated as an integral entity. The relationship between SVF and mean AVF was analyzed in detail (Section 3.4) and shows the dependency between SVF and AVF.
This new method combines approaches from remote sensing and computer graphics and applies them to three dimensional modeling. As the sunlit area is computed directly on the graphics card, modern graphics engines are mighty and cheap tools for correct modeling of the distribution of solar irradiance in complex urban environments. Although, modern graphics engines are mainly designed for the calculation of dynamic lighting of complex scenes with complex shading in computer graphics, they are also useful for scientific purposes. In order to perform the calculations on the graphics card, a threedimensional city model [50] is needed to represent the three-dimensional scene using the graphics engine.
However, the proposed method also has some limitations due to negligence of diffuse soft shading, reflections of radiation by building facades and inbound diffuse shortwave radiation. These shortcomings could be improved by a more sophisticated sky model based on an anisotropic sky vault instead of the proposed solar irradiance with sun rays in parallel. Furthermore, ray tracing could be utilized to account for multiple reflections [5] with the disadvantage of costly calculations. However, just the last generations of graphics cards also allows for the real-time calculation of ray tracing in complex three-dimensional sceneries based on hardware acceleration [66]. As with any numerical spatial analysis, the boundary areas of the study area also pose a challenge here; this is because the spatial influence of structures outside the study area cannot be accounted for in an numeric analysis using either Sky View Factor or Area View Factor [67,68].
The spatial temporal distribution of solar radiation intensity is important for urban planers and landscape architects in order to select optimal species (in terms of resilience to drought stress, photosynthetic productivity, and cooling potential) to improve human thermal comfort in urban areas [1]. Most approaches mapping solar irradiance in the complex urban environment are motivated from an energetic and building perspective: photovoltaic power plants, passive heating of buildings and the availability of sunlight [5,21,69,70]. Little attention is payed to the aspects of the urban green infrastructure. This new method might be also useful for assessment of the photovoltaic potential on roofs and building facades in the urban area which is affected by the surrounding urban morphology.
Due to their schematic structure (e.g., 2D/3D Tile Approach), classical UCMs are suitable for scientific questions related to individual, typified street canyons, but not for questions related to specific urban neighborhoods [30,33,34]. In contrast to these classical UCMs, the newly developed method can be used to calculate the short-wave irradiance in a "real" urban morphology based on spatial geodata provided by local authorities. Thus, SkyHelios is not only suitable for scientific questions, but also for urban planning and design. To apply UCMs in urban planning, it is necessary to focus on the requirements of the users in the local authorities [71]. This means that the models should use standardized data (e.g., CityGML [56]) that have already been collected frequently in the municipalities in accordance to the FAIR principles [72] and choose the appropriate model which meets the requirements of the specific question and the available resources. In general, the advantages of different model types should be combined to get the benefits from all types that are necessary to assess the thermal impacts of urban structures.

Conclusions
The estimation of the fraction receiving direct radiation of urban trees, as calculated with the methodology proposed in this framework, seems like a reasonable approach in order to split total leaf area into shaded and sun-lit leaf area. The spatial, geometrical assessment of sun-lit and shaded areas of urban entities based on aligned orthographic views with respect to the position of the sun, is an additional method for modeling the urban radiation regime. AVF is the complement to the local Sky View Factor, as the POV for AVF is aligned according to the position of the sun. In contrast to SVF, AVF is suitable to assess the variability of heterogeneous sunlit surfaces under seasonal and daily conditions. In contrast, SVF is independent of time and its information is limited to single points in space.
Our proposed method is the first stage in our model pipeline to assess the transpirational cooling potential of trees in the urban environment. Based on this first stage, further work has to be done by determining species specific photosynthetic productivity and transpiration rates, quantifying the transpirative cooling effect and finally estimating the linkage of the urban trees to the urban atmosphere.

Acknowledgments:
The authors want to thank the surveyor's office of the city of Freiburg.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: