<figure><img src="image_1.png" /></figure>

HAL
open science

# Detection of Old Diffusive Windthrow Using Low Cost Resources. the Case of Xynthia Storm in the Vosges Mountains, 28 February 2010

Ionel Haidu, Paula Furtuna, Sébastien Lebaut

► To cite this version:

Ionel Haidu, Paula Furtuna, Sébastien Lebaut. Detection of Old Diffusive Windthrow Using Low Cost Resources. the Case of Xynthia Storm in the Vosges Mountains, 28 February 2010. 2020. [hal-02488288](https://hal.univ-lorraine.fr/hal-02488288)

HAL Id: hal-02488288

https://hal.univ-lorraine.fr/hal-02488288

Preprint submitted on 22 Feb 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
1 Article

2 Detection of Old Diffusive Windthrow Using Low

3 Cost Resources. the Case of Xynthia Storm in the

4 Vosges Mountains, 28 February 2010

5 Ionel Haidu¹*, Paula Furtuna²* and Sébastien Lebaut¹,

6 ¹ Laboratoire LOTERR-EA7304, Université de Lorraine, 57045 Metz Cedex 01, France;

7 ionel.haidu@univ-lorraine.fr, sebastien.lebaut@univ-lorraine.fr

8 ² Faculty of Geography, Babes Bolyai University, 400334 Cluj-Napoca, Romania;

9 paula_roxana.furtuna@yahoo.com

10 * Correspondence: ionel.haidu@univ-lorraine.fr Tel.: +33-605-984-480; paula_roxana.furtuna@yahoo.com

11 Tel.: +40-754-493-296

12 **Abstract:** Unlike the contiguous wind throws, the diffuse wind throws occurred as a result of wind 13 gusts of lower speed (100-140 km/h) than in the first case (>140 km/h) are much more difficult to 14 detect due to their much lower areas and due to their very large number, of several hundreds in the 15 wooded mountain massifs. The objective of this research is to present a rapid procedure for the 16 detection of the diffuse wind throws based on low cost, Landsat type images, knowing that certain 17 sensors cannot be accessed without significant investments. Our application is based on the study 18 of effects caused by the Xynthia storm in the Vosges Mountains in the North-East of France, on 28 19 February 2010. Thus, based on two sets of Landsat satellite images, we used the “dark object” 20 approach and the Disturbance Index, as well as a classification of the images before and after the 21 storm, resulting in a change map. Following the detection process, 257 scattered polygons were 22 detected, totalling 229 ha. For validation purposes, high-resolution images and orthophotoplans 23 taken before and after storm were used. The error matrix was calculated, achieving an overall 24 accuracy of 86%, which confirms the quality of our analysis and supports this procedure for 25 detecting diffuse wind throw based on low cost resources.

26 **Keywords:** wind throw; Xynthia storm; Landsat imagery; limited data.

27

28 **1. Introduction**

29 Storms are known to be the most destructive natural disasters for forests in terms of frequency 30 of occurrence, damages caused and their spatial extent [1,2,3]. Such a meteorological phenomenon 31 occurring in France was the Xynthia storm, on 28 February 2010. Although it did not have the 32 exceptional nature of the catastrophic Lothar storm in 1999 [4,5] or Klaus storm in 2009 [6] by its 33 maximum speed (160 km/h), it caused major damages, i.e. thousands of cubic meters of forest 34 destroyed in Lorraine and Alsace. The Pôle interrégional de la Santé des Forêts du Nord-Est Institute 35 [7,8] mentions for the Lorraine region some 15,000 m³ of destabilized spruce at Louppy-sur-Loison 36 (department of Meuse), some pine forests in the Wardnt Plateau (department of Moselle), and, to a 37 lesser extent, some aged beech forests. In the Alsace region, the same source mentions some limited 38 destabilization in forest border areas or in areas where there is no longer any interest in logging. 39 Besides the official reports [9,10,11,12] the press of the time, i.e. The Vosges Martin newspaper, 40 LOR'Actu, France Culture, but also other local and regional media sources have presented more 41 dramatically the effects of this storm on the forests of Vosges, suggesting that many more points 42 would have been affected, as the wind throw had a diffuse character [13,14]. Even if several years
43

44

45

have passed since that event, through remote sensing images that stored the condition of vegetation,
in this paper we will try to identify the affected points that seem to be masked by the extent of the
forest crop.

46

47

48

49

50

51

The detection of the changes caused by storms, be they contiguous or diffused, is very important
in order to establish spatial and temporal trends in the rehabilitation and management of forest
resources. Papaik and Canham [15] show that storms can be seen as an integral part of the forest
ecosystems functioning, but that the detection of changes must be conducted in order to know and
manage ecosystems and forest resources. The use of remote sensing in forestry applications is today
an extremely useful tool to acquire knowledge and conduct the management of forest ecosystems.
The information provided by this tool on the interaction between the atmosphere and the forest for
the purpose, for example, of locating the damages caused by storms in the context of climate change
[16,17,18,19] or for carbon sequestration potential [20,21] is a must because it can help extrapolate this
information to larger spatial scales.

52

53

54

55

56

Huang et al. [22], proposed a method for the detection of forest areas affected by storms, the use
of vegetation indices and the Tasselled Cap Transformation. Based on this transformation, Healey et
al. [23] developed the Disturbance Index, used for windthrow detection. The idea behind the index
is that disturbance sites exhibit higher brightness, and lower greenness and wetness values compared
to the undisturbed forests [24]. This procedure was then optimized and adapted to the purpose of
continual monitoring of forest disturbances through time series images [25,26,27,28].

57

58

59

60

61

62

63

Another direction for the development of remote sensing to identify wind throws refers to semi-
automatic or automatic procedures. Stach et al. [29] proposed mapping by optical and radar remote
sensing. Huang et al. [30] developed dark object concept and support vector machines to automate
forest cover change analysis, while Kennedy et al. [31] used LandTrendr segmentation algorithms in
order to develop maps of forest disturbance and regrowth.

64

65

66

67

68

Regarding the methodology applied to assess the consequences of Martin and Lothar hurricanes
in 1999 on the forests of France, Stach et al. [32] concluded that the best results were those achieved
by radar remote sensing.

69

70

71

Due to the meteorological risk events causing significant damages throughout recent decades,
France, but also other countries in Western Europe experienced large windthrow areas. Therefore,
the researchers from the French National Forest Inventory [6,33] have developed a new semi-
automatic methodology, the “7 point method” to map the felling damages caused by the storms and
to monitor the land use adapted to varied landform conditions. The principle of the method is based
on the “pixel per pixel” comparison of two satellite images recorded before and after the storm. The
“pixel per pixel” and “spectral band with spectral band” difference is performed to highlight the
radiometric evolution between the images before and after the storm. In order to detect the
wind throws, the “7-point method” combines the change detection techniques with an automatic
image segmentation phase, namely the delimitation of homogeneous areas in terms of radiometry
and texture. The detected forest changes caused by the storm are highlighted, the damages identified
in each plot are classified and, in the end, the damage map for the analysed area is obtained by
combining all the previous results. This methodology has been used to map the wind throws and
damages caused by the Lothar and Martin storms in 1999 in Forêt des Landes, South-Western France.
The windthrow mapping in the Vosges department in the year 2002 was carried out by combining
remote sensing for the plain and hill areas with aerial photogrammetry for the mountain [33]. The
forest damages caused by the Klaus storm occurring in January 2009 were also evaluated using the
7-point method [6]. The working principle of these methods is the achievement of real-time
information by transforming the remote sensing images into vegetation state indices needed to
identify the disturbances or abnormalities occurring in the forest areas [34].

72

73

74

75

76

77

78

79

80

81

82

83

84

85

86

87

88

89

90

91

92

93

94

The wider and wider use of very high-resolution (VHR) remotely sensed images has stimulated
the development of object-based methods and data mining techniques, which were much more
effective for change detection [35]. Bauman et al. [24] have developed a methodology for detecting
and analysing wind throws for various meteorological phenomena. Moreover, they have managed to
separate the disturbances due to the storms as compared to those caused by massive forest cuttings
95

96

97

98

99

100

101

102

103

104

105

106

107

108

109

110

111

112

113

114

115

116

117

118

119

120

121

122

123

124

125

126

127

128

129

130

131

132

133

134

135

136

137

138

139

140

141

142

143

144

in the forest areas, as well as to the other land cover classes. This study was implemented in a temperate zone in Russia and in an area in the USA boreal forests. The methodology applied in this study was based on the detection of wind throws using the Disturbance Index that was designed considering the Tasselled Cap transformation and the "dark object" approach.

The purpose of our application is to detect the forest areas in the Vosges Mountains affected by the Xynthia storm in February 2010, as it is well known that there are diffuse wind throws, without taking into account more sophisticated but also more expensive imaging resources, such as those mentioned above, but, on the contrary, using Landsat's least expensive satellite imagery.

## 2. The study area and the storm event

The study area is represented by the Vosges Mountains (7,360 km², of which 4,400 km² are forest area), located in the North-Eastern part of France (Figure 1), in low altitude mountains, with 14 peaks exceeding 1,000 m and the maximum altitude, of 1,424 m, in the Grand Ballon Peak. The mountains expand longitudinally over approximately 200 km until the German border, along the west bank of the Rhine River. Their maximum width is 60 km in the South, 20 km in the North and only 4 km in their upper third part. The Southern part of the massif consists in crystalline schists, while the Northern part is composed largely of sandstones. The lithology of the area has been heavily influenced by glacial erosion, with many glacial lakes being found here.

Due to the geographic location, the Vosges Mountains experience two climate influences: oceanic and semi-continental. West from the Vosges Mountains, the relief has lower altitudes (Parisian Basin and Lorrain Plateau). Therefore, the air masses propelled by the oceanic disturbances are forced to escalate the Western flanks of the Vosges Mountains generating severe weather conditions and strong wind gusts. However, the storms in the Vosges are not as severe or devastating as those on the Atlantic coast, although they are reactivated by the escalation of the mountains. This is because they weaken in intensity due to the distance of about 1,000 km that these storms have to pass through until the North-East of France, in the exit area from the French territory. The dominant winds are in the Western sector, and in case of storms such as Lothar on 26 December 1999 and Xynthia on 28 February 2010, blow in gusts exceeding 170 km/h on the mountains peaks. This was also the case of some tornadoes, such as the one on 11 July 1984, when some extremely strong gusts were recorded, of 200 km/h.

The massif is predominantly covered by forest (60% of the surface of the mountains today), especially conifers (pine and spruce), which is up to the altitude of 800-1000 m, while towards the foothills of the mountains there is a transition to deciduous forests, as well as to open meadows that were created by man for agricultural practices. The forest cover of the Vosges Mountains ranks second in France in terms of extent, after the Landes forest. It covers the East of Lorraine, the Alsace mountain region and the North-East of the Franche Comté region. The Vosges forest is the first region of France in terms of pine, fir and spruce production, with a crop of 1.5 million m³ of logs, namely 1/3 of the country's total production. Without the human intervention, with the most intense deforestations in the 18th and 19th centuries, the forest would have covered the entire mountain area, except the highest peaks and areas exposed to the weather. Nowadays, the authorities apply the sustainable management principle, so that despite the significant harvest, the forest of Vosges is progressing, the annual harvest being much lower than the growth rate [6].

It is worth mentioning that the Vosges Mountains host 2 regional natural parks (Parc Naturel Régional des Vosges du Nord, Parc Naturel Régional des Ballons des Vosges), Réserve de biosphère transfrontalière des Vosges du Nord-Pfälzerwald (France-Germany), as well as some national natural reserves related to the presence of the forest.

But sometimes weather conditions are so severe that trees are broken or uprooted. The most frequent and extensive wind throws in the Vosges Mountains were caused by wind gusts either due to tornadoes formed during some hot summers or due to winter storms.

The tornadoes monitoring and research group Keraunos, Meteo France, as well as the historical contributions of the CLI.M.A.57-67-68 association and the works of Garnier [37] and [38], enable us
145 to list the most devastating tornadoes mentioned for the Vosges: 13.08.1770 at Gérardmer, 29.05.1807
146 at Lusse, 11.06.1880 at Eloyes, 10.07.1968 a tornado in the Eastern part of Lorraine, 30.06.2012 at
147 Gérardmer and 13.05.2015 at Gerbépal. It is worth mentioning that most often such tornadoes or mini-
148 tornadoes are produced in the Western part of the Vosges Mountains. In remote sensing studies
149 careful consideration should be given when selecting images for the consequences of certain episodes
150 with wind gusts in order to avoid confusion and to correctly identify the cause of wind throws, be
151 they summer tornadoes or winter storms.

152 A very violent convective phenomenon similar to a mini-tornado called "derecho" occurred on
153 16 September 2015 in the Centre and East of France, which affected also the Western part of the
154 Vosges Mountains and caused significant wind throws. Even in the winters with remarkable
155 snowfalls (1969/1970, 1994/1995, 2005/2006) some trees were damaged or blown over by "heavy
156 snow", sometimes in combination with a stronger wind. Rarely, in the case of sudden snow removal,
157 stone and mud torrents can be formed that damage or break trees in their path.

<figure><img src="image_3.png" /><figcaption>158
159 **Figure 1.** Location of Vosges Mountains and a sample area highlighted later.</figcaption></figure>

160 Due to location, climatic conditions and topography, the Vosges Mountains area has often been
161 affected by major winter storms that have led to wind throws on extensive forest areas, causing the
162 greatest damages if we also relate to other damaging causes for forests. Meteo France distinguishes
163 12 types of storms, [39] with the most frequent being the North and West oceanic storms. Hence, the
164 Vosges Mountains area, especially the Northern part, is on the storm's travel path. Obviously, the
165 French coastline is the most exposed and feels the maximum strength of any storm originating in the
166 Atlantic. Then, the storm's strength gradually weakens towards the interior of the country. This
167 explains why wind throws in the North-Eastern France and the ones from the Vosges Mountains have a local character or are often diffuse.

168
The Xynthia winter storm occurred from a low pressure system formed in the subtropical area
169
south of the Azores on 25 February 2010. The southern flow of cold air masses from the altitude
170
caused the formation of a depression over the Centre and Eastern part of the North Atlantic Ocean.
171
The high temperature difference between the extremely hot air flow above the Africa and the cold air
172
in the Eastern Atlantic Ocean has accentuated the decrease in atmospheric pressure that led to the
173
formation of the powerful Xynthia cyclone. The cyclone experienced a rapid intensification affecting
174
8 countries, and causing 65 deaths [40] in its movement from SW to NE (Figure 2). The depression
175
core reached the Gascogne Bay on the evening of 27 February bringing violent winds from the
176
Southern sector that blew in the Pyrenees Mountains at 242 km/h. On the morning of 28 February,
177
Xynthia furiously hit the French coastline, the pressure at the centre of the cyclone dropped to 968
178
hPa and the wind gusts reached 120-160 km/h. Following its path, Xynthia crossed the South-West
179
and Central France, approaching the Vosges Mountains from the South-West direction on 28
180
February 2010. Wind gusts of 172 km/h were recorded in the Markstein Resort (68) located at 1184 m
181
altitude, but also the other mountain or resorts in the neighbourhood of the Vosges Mountains
182
recorded notable wind gusts: Septsarges (55) 147 km/h, Goin (57) 136 km/h, Metz (57) 136 km/h,
183
Belmont (67) 134 km/h, and Vagney (88) 133 km/h. According to the classification cited above,
184
Xynthia was of the SW type, with an extremely southern origin (close to the Tropic of Cancer) for a
185
typical winter storm. It had an atypical trajectory and its effect was felt in North-Eastern France
186
exceeding its borders, until the Baltic Sea, which showed that winter storms could threaten the entire
187
country. From a statistical point of view, Bedach and Hofherr [40] estimate that storms as severe as
188
Xynthia have an average occurrence frequency in France once every two years.
189
The existing literature [4,37,38,39,41,42] allowed us to locate Xynthia in a series of devastating
190
storms that have hit the Vosges Mountains (Table 1).
191

<table><caption>**Table 1.** The main recorded winter storms that caused wind throws in the Vosges Mountains</caption><thead><tr><th>Date</th><th>Affected area</th><th>Notable characteristics</th></tr></thead><tbody><tr><td>1698</td><td>N-E of France</td><td></td></tr><tr><td>1716</td><td>N-E of France</td><td></td></tr><tr><td>1739, 14-18.01</td><td>N-E of France</td><td>Major forest damage in the Vosges Mountains</td></tr><tr><td>1836, 23.03</td><td>All France</td><td></td></tr><tr><td>1866, 10-11.01</td><td>1866, 10-11.01</td><td></td></tr><tr><td>1879, 20.02</td><td>All France</td><td>130 km/h gusts, Meurthe-et-Moselle very affected</td></tr><tr><td>1882, 26-27.03</td><td>All France</td><td></td></tr><tr><td>1886, 27.02</td><td>N-E France</td><td></td></tr><tr><td>1892, 29-30.03</td><td>Alsace</td><td>0,5 M m³ felled trees</td></tr><tr><td>1896, 04.12</td><td>N-E France</td><td></td></tr><tr><td>1900, 13.02</td><td>Lorraine</td><td></td></tr><tr><td>1902, 31.01-01.02</td><td>Vosges</td><td>2,2 M m³ felled trees</td></tr><tr><td>1920, 10-14.01</td><td>Moselle, Alsace</td><td>1,0 M m³ felled trees</td></tr><tr><td>1924, 08.10</td><td>N-E France</td><td></td></tr><tr><td>1935, 22-23.02</td><td>N-E France</td><td></td></tr><tr><td>1940, 14.03</td><td>Lorraine</td><td>180 km/h wind gusts in Nancy</td></tr><tr><td>1952, 12-13.12</td><td>Lorraine, Alsace</td><td>112 km/h wind gusts in Metz-Frescaty, major wind throws</td></tr><tr><td>1958, 02</td><td>Vosges</td><td>1.2 M m³ felled trees</td></tr><tr><td>1962, 16.02</td><td>Lorraine, Alsace</td><td>151 km/h wind gusts at Metz</td></tr><tr><td>1967, 02</td><td>N-E of France</td><td>130 km/h wind gusts at Bale-Moulhouse; 3.6 M m³ felled trees</td></tr><tr><td>1990, 25.01</td><td>N-E of France</td><td>DARIA storm, 137 km/h wind gusts at Nancy-Ochey</td></tr><tr><td>1990, 03.02</td><td>Almost all France</td><td>HERTA storm</td></tr></tbody></table>
<table><tr><td>1990, 26-28.02</td><td>Almost all France</td><td>VIVIANE storm</td></tr><tr><td>1990, 03</td><td>Lorraine, Alsace</td><td>WIEBKE storm, 144 km/h wind gusts at Metz</td></tr><tr><td>1994, 28.01</td><td>Landes, Alsace</td><td>1.2 M³ felled trees</td></tr><tr><td>1999, 26-28.12</td><td>All France</td><td>LOTHAR storm, 155 km/h wind gusts at Nancy, 140 M<br>m³ felled trees from which Lorraine 29.5 and Alsace 6.5</td></tr><tr><td>2006, 03.10</td><td>Alsace</td><td>RENATE storm, 149 km/h wind gusts at Markstein</td></tr><tr><td>2006,08.12</td><td>Alsace</td><td>VERA storm, 136 km/h wind gusts at Markstein</td></tr><tr><td>2009, 9-10.02</td><td>Alsace</td><td>QUINTEN storm, 166 km/h wind gusts at Markstein</td></tr><tr><td>2010, 28.02</td><td>Almost all France</td><td>XYNTHIA storm</td></tr><tr><td>2011, 15-18.12</td><td>N of France</td><td>JOACHIM storm, 192 km/h wind gusts at Markstein,<br>limited wind throws, dispersed in the coniferous<br>massifs</td></tr><tr><td>2013, 27-28.10</td><td>N-E of France</td><td>CHRISTIAN storm, 162 km/h wind gusts at Markstein</td></tr><tr><td>2014, 13.02</td><td>N-E of France</td><td>TINI II storm, 151 km/h wind gusts at Markstein,<br>17,000 m³ diffuse wind throws</td></tr></table>

193

194

195

196

197

198

199

200

201

202

As Bock et al. [43] stated, the institutes in France involved in the management of the forest resources - ONF (Office National des Forêts), IDF (Institut pour le Développement Forestier) and the INRA (Institut National pour la Recherche Agronomique), have studied the dependence between wind throws and wind speed. In regular beech forests, isolated wind throws occur at a wind speed ranging between 100 and 120 km/h, multiple diffuse wind throws occur at a 120-140 km/h wind speed, and large scale and frequent wind throws occur at speeds over 140 km/h. The same authors also found that at the height of the trees below 23-24 m the wind throws were rarer, but at wind speeds over 120 km/h the frequency of wind throws increased proportionally to the height of trees.

<figure><img src="image_4.png" /><figcaption>Figure. 2. The trajectory of the Xythia storm, (source [39])</figcaption></figure>

203

## 3. The alternative of using low cost resources

204

205

206

207

208

209

210

211

212

Nowadays, in the context of the impetuous development of remote sensing techniques, there is a natural preference for high and very-high-resolution (VHR) remote sensing images. Unfortunately, these are expensive and more difficult to access. We also have to admit that sometimes it is not possible to capture real-time images or create time series of images. Maybe this would be possible, but the cost of raw data purchases is too high. Therefore, we would like to assess if less expensive and easier-to-acquire images, such as Landsat, despite their lower resolution, would still be sufficient and capable of detecting diffuse wind throws. Starting with 2014 Sentinel Radar images are also free, but obviously these cannot be considered because the event we studied occurred in 2010. One should mention that in the case of the free images, very expensive software is often needed to process them.
But if the purpose of a research is to rebuild the old windthrow trees, obviously, only the sources produced by the older satellites or having an archive (Landsat from 1980) can be accessed.

The availability of satellite data varies and choosing the best sensor is very important in order to select the data which better suit our purpose and to get more accurate results. When detecting wind throws using optical satellite imagery it is important to know which sensor is appropriate. Schwarz et al. [44] summarized the criteria to be considered in choosing the appropriate sensor, namely: *data availability, spectral resolution, spatial resolution and temporal resolution* (to increase the probability of achieving useful cloud-free data). Besides availability and resolution, the cost of data represents also an important factor to be considered. Satellite imagery products vary in their spatial and spectral resolution, geographic and temporal coverage, cloud cover, security regulations, and price [45]. The most frequently used and inexpensive data sources for detecting wind throws are the remote sensing images offered by the Landsat sensor [48]. These are basic sources due to the free of charge images and to the spatial, spectral and temporal details they provide for the mapping of the area damaged by storms.

The free of charge availability of Landsat archive and the large geographic coverage provide the opportunities to conduct studies over a long period of time and over a large area [47-49]. Beside the medium resolution provided, Landsat images at a spatial resolution of 30 m were successfully used to detect wind throws in various areas [23,24,27,50-53]. However, the increasing availability of high-resolution satellite imagery offers the potential to acquire data with fine spatial resolution (e.g. WorldView, Ikonos, Quickbird), but the cost of these data is a hindering factor in using them. Jackson et al. [54] showed that high-resolution imagery, like LIDAR, offered a great potential to detect wind throws, but they required more basic operating instructions.

In this study we chose to use Landsat data to detect the wind throws along the Vosges Mountains. The first reason in choosing Landsat sensor was the wide swath and the low cost. The medium resolution Landsat data is not as sharp as high-resolution images, but the amount of information it contains along with the large area coverage make this data suitable for large- or small-scale studies. We studied the diffuse wind throws and the exact location of these disturbances was unknown. Therefore, the availability and the coverage represented another factor in choosing the sensor. The Landsat TM and ETM+ program covers the entire world within 16 days and the images are saved in an archive, so the possibility to find cloud-free images is increasing due to this availability of free data [44]. Each type of satellite images has advantages and disadvantages in detecting wind throws. Malthus et al. [55] have developed a comparative study to emphasize the advantages, limitations and even prices, where possible, of current remote sensing image sources used in *commercial forestry*. Similar comparisons have been developed by Suarez et al. [56], which compared optical techniques with RADAR and LIDAR sensors. Melesse et al. [57] presented the characteristics of 20 different satellite sensors, while [3] summarised the specification of sensors used for *forest disturbance detection*.

The table below (Table 2) summarizes the characteristics of the most important satellite data used in wind throw detection and forest cover monitoring. The most important differences between these data are the spatial resolutions, e.g. from 30 m for Landsat-7 multispectral data to 1 m and even under 1 m for the Ikonos, WorldView, and WuickBird data. But when detecting and assessing wind throws, the choice of the appropriate sensor will depend on the availability of data, cost, scale of the fragmented landscape and area of interest, because the high resolution imagery is not available globally, like Landsat data. If the exact small scale (< 0.5 ha) delineation of the wind throw is needed, then the high-resolution imagery is suitable. But if the goal is to detect large scale wind throws quickly, then the medium resolution data is more useful [44].

<table><caption>Table 2. Characteristics of a number of remote sensing sensors</caption><thead><tr><td>Sensor</td><td>Years</td><td>Number of Bands</td><td>Multispectral resolution (m)</td><td>Panchrom. resolution (m)</td><td>Repeat cycle</td><td>Scene size (km)</td><td>Cost per sq. km</td></tr></thead></table>
<table><thead><tr><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th></tr></thead><tbody><tr><td>Landsat 4-5</td><td>1982</td><td>4</td><td>80</td><td>—</td><td>18days</td><td>185 x 185</td><td>Free</td><td>8 of 34</td></tr><tr><td>MSS</td><td>1999</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Landsat 4-5 TM</td><td>1982</td><td>7</td><td>30</td><td>—</td><td>16 days</td><td>185 x 172</td><td>Free</td><td></td></tr><tr><td>Landsat 7 ETM+</td><td>1999</td><td>8</td><td>120 - thermal</td><td>15</td><td>16 days</td><td>185 x 172</td><td>Free</td><td></td></tr><tr><td>Landsat 8</td><td>2013</td><td>11</td><td>30</td><td>15</td><td>16 days</td><td>185 x 180</td><td>Free</td><td></td></tr><tr><td>ASTER</td><td>1999</td><td>14</td><td>15 - band 1-3<br>30 - band 4-9<br>90-band10-14</td><td>—</td><td>1-2 days</td><td>60</td><td>Free<br>(on-demand)</td><td></td></tr><tr><td>MODIS Terra</td><td>1999</td><td>36</td><td>250 - 2 bands</td><td>—</td><td>1-16 days</td><td>2330</td><td>Free</td><td></td></tr><tr><td>MODIS Aqua</td><td>2002</td><td></td><td>500 - 5 bands<br>1000-29bands</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>SPOT 1-3</td><td>1986<br>1997</td><td>3</td><td>20</td><td>10</td><td>26 days</td><td>60 x 60</td><td>$1,200<br>per<br>scene</td><td></td></tr><tr><td>SPOT 4</td><td>1998<br>2013</td><td>4</td><td>20</td><td>10</td><td>26 days</td><td>60 x 60</td><td>$1,200<br>per<br>scene</td><td></td></tr><tr><td>SPOT 5</td><td>2002</td><td>4</td><td>10</td><td>2.5/5</td><td>26 days</td><td>60 x 60</td><td>$2,700<br>per<br>scene</td><td></td></tr><tr><td>SPOT 6-7</td><td>2012</td><td>4</td><td>6</td><td>1.5</td><td>26 days</td><td>60 x 60</td><td>$5.15</td><td></td></tr><tr><td>IKONOS</td><td>1999</td><td>4</td><td>4</td><td>1</td><td>14 days<br>Less than<br>3 days</td><td>11.3</td><td>$30-53<br>per<br>scene</td><td></td></tr><tr><td>QuickBird</td><td>2001</td><td>4</td><td>2.4</td><td>0.61</td><td>Between<br>1 and 2.8<br>days<br>daily</td><td>16.8</td><td>$18</td><td></td></tr><tr><td>FORMOSAT 2</td><td>2004</td><td>5</td><td>8</td><td>2</td><td>24x24</td><td></td><td></td><td></td></tr><tr><td>Pléiades 1A-1B</td><td>2011</td><td>4</td><td>2</td><td>0.5</td><td>26 days<br>daily</td><td>20</td><td>$13-40</td><td></td></tr><tr><td>WorldView-1</td><td>2007</td><td>1</td><td>—</td><td>0.46</td><td>3.7 days (for 0.52 m)<br>1.1 days (for 1 m)</td><td>17.6 x 14</td><td>$14-48</td><td></td></tr><tr><td>WorldView-2</td><td>2009</td><td>8</td><td>1.84</td><td>0.46</td><td>Up to 1.1 days</td><td>16.4</td><td>$14-58</td><td></td></tr><tr><td>WorldView-3</td><td>2014</td><td>28</td><td>1.24</td><td>0.31</td><td>between<br>1 to 4.5<br>days</td><td>13.1</td><td>$19</td><td></td></tr><tr><td>WorldView-4</td><td>2016</td><td>4</td><td>1.24</td><td>0.31</td><td>less than<br>3 days</td><td>—</td><td>—</td><td></td></tr><tr><td>RapidEye</td><td>2008</td><td>5</td><td>5</td><td></td><td>between<br>1 to 5.5<br>days</td><td>25x25</td><td>$14-55</td><td></td></tr><tr><td>Sentinel 1A<br>and 1B</td><td>2014,<br>2016</td><td>1 SAR<br>band</td><td>variable</td><td></td><td>5 to 12<br>days</td><td>—</td><td>free</td><td></td></tr><tr><td>Sentinel 2A<br>and 2B</td><td>2015,<br>2017</td><td>13</td><td>variable</td><td></td><td>5 to 10<br>days</td><td>—</td><td>free</td><td></td></tr></tbody></table>

260

* Source [60,61]

261 ** Prices are only estimates based on online sources, and may vary and the price can increase due to

262 the demanding rules set by the providers. E.g. - usually the minimum order area for archive

263 imagery is 25 sq. /km with a 2 km (digital globe) or 3 km (GeoEye) minimum order width.
The free access policy for Landsat data has made this data more and more frequently used in remote sensing studies [60,61]. However, with the accessibility of these data, there is a challenge to find images without clouds or with very little cloud coverage in the analysed area, because clouds can significantly influence the reflection from different spectral bands [62]. This is an issue we have also faced and we will present next how we have solved this shortcoming

## 4. Materials and Methods

In this study we aimed to highlight the applicability and efficiency of Landsat data in detecting diffuse wind throws caused by the Xhyntia storm that crossed the Vosges Mountains on 28 February 2010. Wind throw detection depends on the affected area. If the affected area is larger, the detection capability increases, together with the results accuracy of the Landsat satellite images based analysis. Figure 3 illustrates the methodological workflow.

<figure><img src="image_5.png" /><figcaption>Figure 3. Methodological steps</figcaption></figure>

275

276

We started with the pre-processing of satellite imagery, namely converting them from digital numbers into reflectance. The next step consists in masking the water bodies and dark soil to avoid these pixels being confused with pixels corresponding to the vegetation areas, by using the NDVI index values. Then the images were classified into forest and non-forest areas. The classification was made for the two images, before and after the storm, respectively. Then the difference between the classified images was made, resulting into the identification of areas with lost forest vegetation.

In addition, to improve the result obtained in the previous step, the Disturbance Index (DI) was used. The calculation of DI is made based on the Tasselled Cap bands. After overlapping the two resulting images we obtained the forest changes map. Then, the validation of results was carried out. The evaluation of results depends on the data and information used. The accuracy can be directly influenced by the quality and availability of the ground data (the resolution of orthophotos, high-resolution satellite imagery). The first step in results validation is to calculate the accuracy assessment and the Kappa statistic. Another validation was made by overlapping the polygons resulting from the change map with the orthophotos from 2009 and with the RapidEye satellite image from 2010.

### 4.1. Pre-processing, land-water and cloud mask

Landsat satellite images were used for wind throw detection. The images were purchased for the previous year and the year after the meteorological phenomenon (28 February 2010) that has led to the fall of trees in the Vosges Mountains. Thereby, a multi-temporal set of Landsat 5 TM images (with 7 spectral bands - 30 m spatial resolution) was used. The images were acquired on different days 17.06.2009, before the Xynthia storm, and on 4.06.2010, after the phenomenon. The Landsat images were downloaded freely from the [63] website. Table 3 shows that between the two dates of image selection there were no meteorological events with extreme wind gusts, no storms or tornadoes. Therefore, the data reflected by the two images cannot be artefacts. We selected no clouds or very
few clouds images that were recorded during the growing season and these were purchased for the same vegetation season (June), to reduce the impact of vegetation in different seasons.

<table><caption>Table 3. The characteristics of the images that were used for the Vosges Mountains area</caption><thead><tr><td>Date of the Xynthia storm</td><td>Path/Row</td><td>Lat/Long</td><td>Previous to Xynthia</td><td>After Xynthia</td></tr></thead><tbody><tr><td>27-28.02.2010</td><td>196/026-027</td><td>48.9/6.9</td><td>17.06.2009</td><td>4.06.2010</td></tr></tbody></table>

In order to use the satellite images to detect changes occurring at the ground surface, several steps are required, namely: selecting the images that best fit the intended purpose, orthorectification, georeferencing, various corrections, calculating some indices to enhance the data information [64]. The Landsat satellite images provided by USGS have already been radiometrically and geometrically corrected, but for the proposed operations we needed the values of the reflective satellite images. This operation was done using the formulas given by Chander & Markham [65]. The first step was to pre-process the two sets of images, before and after the storm, by converting the digital numbers into spectral radiance using the equation:

$$L_{\lambda} = (\text{gain} * \text{DN}) + \text{bias}_{\lambda} \qquad (1)$$

Where:

$L_{\lambda}$ - the spectral radiance at the sensor level and referring to the wavelength $\lambda$ for the specific band;

DN - Digital Number;

gain and bias - are specific calibration parameters determined before launch.

For the relatively clear Landsat images, a reduction in variability of images can be made by converting spectral radiance, calculated using equation (1), into surface reflectance or albedo. To obtain surface reflectance, the fraction of solar energy reflected on a certain wavelength, the following equation was used:

$$R_{\lambda} = (\pi * L_{\lambda} * d^2) / (ESUN_{\lambda} * \sin\theta_{s}) \qquad (2)$$

Where:

$R_{\lambda}$ - Spectral radiance at the sensor's aperture;

$L_{\lambda}$ - Spectral radiance at the sensor's aperture calculated through equation 1;

$d$ - Earth-Sun distance [astronomical units];

$ESUN_{\lambda}$ - Mean exoatmospheric solar irradiance [66];

$\sin\theta_{s}$ - solar zenith angle in degrees

The next step consists in masking the clouds, as we mentioned previously. This step is essential because during the detection process, the pixels corresponding to the clouds or their shadows can be confused with darker vegetation pixels, such as the coniferous forests. Clouds are difficult to detect on Landsat images, due to the spectral properties of these images and the complexity of the clouds [67]. It may be more difficult to distinguish the shadow of clouds from "false shades" [64] given by the dark soils, the topographic shadows, water, etc. For example, Huang et al. [26] reported difficulties in separating shadows from the spectral response of waters (cited by [62]). On the other hand, the lack of adequate images may introduce errors in the process of masking the clouds and lead to the inability to differentiate between the cloud shadows and the topographic shadows where there are areas with shadow potential [61,68]. If the areas with changes are relatively small, this error time significantly decreases the detection accuracy of the changes [67]. The clouds and associated shadows were separated through a classification of associated pixels and excluded from the subsequent analyses of the changes occurred in the studied area. After applying this mask, the areas where clouds and shadows were located received the value "No data", not taking into account subsequent processing.
In addition to the masking of clouds and their associated shadows, it is necessary to apply another mask for some non-forest dark surfaces such as water, dark soil, and dark impervious surfaces because these can be as dark as or even darker than forest pixels [26]. To identify and remove these pixels we used the method proposed by Huang et al. [26], which is based on the spectral properties of water bodies and dark soil. Pixels are classified as water or dark soil if they have reduced reflection in the SWIR band and if they have low NDVI index values calculated based on the RED and NIR reflectance bands:

$$NDVI = (NIR-RED)/(NIR+RED) \quad (3)$$

After calculating the Normalized Vegetation Difference Index, a threshold was set to delineate the land-cover classes from the masked areas covered by water and soil, respectively. In this study the threshold 0.25 was used. To exemplify the result we performed a zoom on a 26.6 km x 20.5 km sample area, i.e. on a 54,613 ha area, which has the Remiremont city in the centre (Figure 4), a sample that was also used for the other examples. This step made it possible to identify and eliminate pixels that may be problematic throughout the next steps.

<figure><img src="image_6.png" /><figcaption>Figure 4. NDVI values (a) and masking of non-vegetated dark objects (b)</figcaption></figure>

## 4.2. Training and image classification

In order to obtain the forest areas, the image before the Xynthia storm must be classified into “forest” and “non-forest”. This is carried out following two steps: automatic extraction of forest pixels based on the histogram, and by conducting a supervised classification of the two images before and after the occurrence of the phenomenon.

### 4.2.1. Automatic extraction of forest pixels

The forest pixels have been identified using the “dark object” approach [30] which involves the delineation of the “forest peak”, based on a band local histogram (Figure 5). In the absence of non-vegetated dark objects (which were masked previously), the forest pixels will be located to the left of the histogram, while the first peak on the left will be considered the “peak forest”.

Huang et al. [30] recommends taking into account two major factors: i) delineation should be done for a local image window, because this area will contain fewer classes of land cover and the delineation of forest classes from other classes will be easier to achieve.

ii) selection should be made of the appropriate band to create the histogram which sets the basis for class delineation.
<figure><img src="image_7.png" /><figcaption>a) satellite image in 5-4-3 combination</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>b) delineation of "forest peak" using RED band histogram</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>c) forest pixels (white) after applying the 0.25 threshold</figcaption></figure>

<caption>Figure 5. The procedure for the delineation of forest pixels based on the histogram for the local image window, (adapted after Huang et al. [30])</caption>

In order to determine whether pixels have been classified properly into forest and non-forest, we have set samples of black pixels representing vegetation based on which we calculated the Integrated Forest Index (IFI) [24, 30].

$$IFI = \sqrt{\frac{1}{NB} \sum_{i=1}^{NB} \left( \frac{b_{pi} - \bar{b}_i}{SD_i} \right)^2} \quad (4)$$

where:

$\bar{b}_i$ and $SD_i$ are the mean and standard deviation of forest training pixels within the image for band i;

$I$;

$b_{pi}$ is the band i spectral value for pixel $p$;

NB is the number of spectral bands.

393

Low IFIs values indicate a high probability of wooded areas, while high values indicate high probability of other land cover categories. Certain dark pixels corresponding to crops may be spectrally similar to pixels for forest areas, but will have IFIs values lower than the forest areas. To eliminate them, it is necessary to set a threshold value, which separates the forest areas from the other areas. For the present study, the threshold was set to 3, so pixels smaller than or equal to 3 are considered as belonging to the wooded areas.

### 4.2.2. Supervised classification

The general objective of the image classification procedure is to classify all the pixels from an image in classes of ground cover or themes. The purpose of the classification in this study is to complement the result of extracting forest pixels by means of the histogram. First of all, land cover classes (water bodies), forest areas (conifers, mixed and deciduous trees), urban areas, pastures and agricultural areas were established and the pixels were selected for the classification algorithm to acknowledge the classes. After running the classification algorithm (Figure 6a), the resulting images have been reclassified into forest areas and non-forest areas. Finally, following the classification of both images from 2009 and 2010, the separation of the forest areas from the non-forest areas was achieved (Figure 6b). This will be used in the end to obtain a change map for the two analysed temporal frames.

411
<figure><img src="image_10.png" /><figcaption>Figure 6. Supervised classification (a) and forest/non-forest classification (b).</figcaption></figure>

412 4.3. Forest Disturbance Detection

413 The detection of forest disturbed areas was carried out using the Disturbance Index (DI) 414 developed by Healey et al. [23] and applied by [24, 69-75]. DI was developed to highlight the spectral 415 signatures of non-vegetated areas associated with disturbance zones and to separate them from all 416 other spectral signatures of the forest areas. The procedure is based on a combination of bands 417 resulting from the Tasselled Cap Transformation, which reduces the six bands of Landsat TM image 418 (excluding the thermal band) into three bands that render brightness, greenness, and wetness. 419

420 This index highlights the difference in spectral signatures and thus disturbed forest areas have 421 higher brightness values and lower values for greenness and wetness bands compared to 422 undisturbed forest areas. The index is calculated based on the formula: 423

424 425 426 427 428 429 430 431 432 433 DI = Br - (Gr + Wr) (5)

424 Where: Br, Gr, Wr are the Tasselled Cap bands, standardized (resized) around the scene's mean forest 425 value. 426

427 The result of the calculation is a single band, where positive values generally indicate areas with 428 disturbance. In order to separate the disturbed areas from the non-disturbed ones, it is necessary to 429 identify a limit threshold, depending on the environmental conditions in the studied area. In our case, 430 the value of DI = 2.5 best identifies the areas affected by wind throws. 431

432 The last step involved the overlapping of the initial image classified as “forest areas” and “non- 433 forest areas” with the map derived from the DI calculation, and so we obtained a change map for the 434 forest area in the Vosges Mountains in 2009-2010 (Figure 7).
<figure><img src="image_11.png" /><figcaption>Figure 7. Areas affected by diffuse windthrow.</figcaption></figure>

As seen in Figure 7, diffuse wind throws occur in the form of scattered polygons, or rather points, on the territory of the Vosges Mountains, some points being so small that they cannot be visually perceived at the scale of cartographic representation in the figure above. We identified 257 polygons (list available in Table S1: Polygons obtained after running the detection algorithm), affected by wind throws, with a total area of 229 ha. We attach additional material to indicate diffuse windthrow (Appendix A).

## 5. Results and discussions

In the algorithm we have developed, we included neither the object-based image analysis nor the segmentation procedures. We will explain next why we particularly made this decision. Segmentation techniques are recommended to detect possible changes, especially for high-resolution satellite imagery. Comparing the results obtained by the segmentation of the Landsat images with high-resolution images, it was revealed that Landsat results provide a low accuracy [76]. The classes are not well defined. The analysis is more difficult and requires a high user experience [77]. In addition to the resolution characteristics, there are other factors such as geomorphological structure, topography, and dead vegetation [78], that affect the applicability of segmentation on Landsat

434

435

436

437

438

439

440

441

442

443

444

445

446

447

448

449

450
images. The algorithm we propose is similar to the French methodology known as the "7-point method" [33]. However, this still includes different stages, proving to be effective in detecting diffuse wind throws based on Landsat images.

The difference between this method and the "7 points" method is given by the steps followed and the approaches used. The first difference is given by the satellite images used. This study is based on medium-resolution Landsat satellite images, while the "7 points" method uses high resolution images. In the "7 points" method, "pixel per pixel" is used in order to detect the changes. The same "pixel per pixel" approach was also used in this study, but it was applied to forest/non-forest images. The next stage in the "7 points" approach consisted in the segmentation of the disturbed areas. As Darwish et al. [76] stated, the segmentation techniques on the Landsat images emphasized the fact that the results provided a low accuracy. Therefore, this method is not recommended for medium resolution images. Instead, in this study we used the Disturbance Index to detect the affected areas. This index is calculated based on the bands obtained from the Tasselled Cap Transformation and it highlights the difference in spectral signatures. Thereby, the disturbed forest areas have higher brightness values and lower greenness and wetness values compared to undisturbed forest areas.

The last step in the "7 points method" is to integrate the results obtained by segmentation with the soil results. Thus, by classifying the disturbances produced, a map of damages was obtained. Due to lack of ground data, we have chosen to do a validation based on the error matrix. Thus, we have calculated the accuracy for the supervised classification into forest and non-forest categories, as well as for the map of changes in forest areas. This map was obtained by overlaying the two resulting maps (from the difference between the forest/non-forest image and the image resulting from the DI calculation). For validation, we overlapped the resulting polygons from change map on the 2009 orthophotos and the 2010 RapidEye satellite image (Figure 9).

## 5.1. Accuracy Assessment

To evaluate the quality of the analysis, we assessed the overall accuracy and the Kappa coefficient. Following the classification of the two images, the accuracy assessment was performed in two steps: a) the accuracy of supervised classification and of forest/non-forest classification and b) the accuracy of the forest change map.

To evaluate the accuracy of the supervised classification, we used 100 random points for each class we used. These points were manually labelled using Quickbird Google Earth and RapidEye high-resolution imagery [79]. These high-resolution images were used only for the confirmation and validation of the classification. Also, for results verification and validation purposes, orthophotoplans from 2009 were used, with a resolution of 2.5 and 5 m, available on [80].

The polygons resulting from the detection process were overlapped with the orthophotos to see whether there is vegetation in the detected areas. The result was summarized in an error matrix, while calculating the overall accuracy and the Kappa coefficient. Therefore, the accuracy assessment calculated for the forest/non-forest classification is given by the overall accuracy = 94.14% and by the Kappa coefficient = 0.92. The evaluation of the forest changes occurred between 2009 and 2010 was carried out based on a sample of 300 points. The previously defined samples considered as real elements in the field were used to calculate the error matrix, the overall accuracy, the kappa statistic and the user's and producer's accuracy for each class, as shown in Table 4.

For the change map we used three classes, namely *constant forest* (CF) that did not change in the two images, *disturbed areas* (D) and *other non-forest classes* (NF). The constant forest area is given by those areas that were represented by forest both in the *before Xynthia* storm image and in the *after Xynthia* image, while the areas with disturbances are those areas that were classified as forest in the *before Xynthia* image and that appeared as non-forest in the *after Xynthia* image. The overall accuracy obtained for the change map was high, i.e. 86%, and the value of the Kappa coefficient was 0.84. Also, the class represented by the constant forest had the user's and producer's accuracy higher than the disturbance class.

<caption>**Table 4.** The accuracy assessment calculated for the change map</caption>
501

(CF=constant forest class; D=disturbance class; CNF=constant non-forest class)

<table><thead><tr><th rowspan="2">Overall accuracy (%)</th><th rowspan="2">Kappa statistic</th><th colspan="3">User's accuracy (%)</th><th colspan="3">Producer's accuracy (%)</th></tr><tr><th>CF</th><th>D</th><th>CNF</th><th>CF</th><th>D</th><th>CNF</th></tr></thead><tbody><tr><td>86</td><td>0.84</td><td>89.10</td><td>78</td><td>88</td><td>81</td><td>80</td><td>88.47</td></tr></tbody></table>

502

## 5.2. Results validation

503

504

505

506

507

508

509

510

511

512

513

Once the accuracy of the thematic data has been set, the validation of the results was the next step. Detecting forest surface changes using satellite imagery is a difficult task and requires validation using ground data. The validation was conducted based on disturbance data occurred in 2010 in the mountains area (data available on [81], referring to all types of disturbances, not only wind throws), on the cadastral forest plots for the two years (data available on [82], as well as on the orthophotos of 2009 and the high-resolution satellite imagery of 2010. For the final validation we used the change map obtained previously and we combined it with data on disturbed forest areas, according to the French National Forest Inventory (IFN). The resulting polygons during the detection process were overlapped with the disturbed areas. Hence, we obtained a first validation of the results. Figure 8 highlights one of the polygons detected from the image processing, overlapped with the Landsat satellite image of 2009 (Figure 8b), and the Landsat image in 2010, respectively (Figure 8c).

<figure><img src="image_12.png" /></figure>

514

515

516

517

518

<caption>**Figure 8.** Forest change map (a) overlapped with the polygons resulting from the detection process, as well as the points where wind throws occurred in 2010 (according to INF): one polygon overlapped with the 2009 satellite image (b), and the 2010, respectively (c), and a zoom on the same polygon (d).</caption>

519

520

521

522

Also, in order to obtain a picture of the “ground-truth”, this polygon was overlapped with the orthophoto of 2009, and with the 2010 RapidEye high-resolution image (Figure 9), thus achieving a new validation.
<figure><img src="image_13.png" /><figcaption>**Figure 9.** Overlapping a polygon extracted through the detection process with the 2009 orthophoto (a) and with the satellite image in 2010 (b).</figcaption></figure>

Another validation was carried out using high-resolution satellite imagery available in Google Earth. For an easier interpretation of the effect of the Xynthia storm, we zoomed out to a sample territory of 26.6 km x 20.5 km, i.e. 54,613 ha, which is in the centre of the Remiremont city (Figure 10).

Hence, Figure 10(b) highlights another polygon obtained as a result of the disturbed areas detection process. This polygon is superimposed on the change map (in raster format). In the wind throws detection phase, raster data were used which were subsequently converted into vector data in order to overlap the final results. By converting these images, a comparison between the results detected on raster format (blue colour) and vector representation (purple polygon) was achieved. The geometric effect of this polygon obtained by switching from raster format to vector format suggests a possible underestimation of the actual area damaged by the wind. Figure 10 (d) shows the image *before* the Xhyntia storm, while Figure 10 (c) illustrates the forest state *after* the Xhyntia storm.

523

524

525

526

527

528

529

530

531

532

533

534

535

536

537
<figure><img src="image_14.png" /><figcaption>**Figure 10.** The forest change map overlaid with the polygons resulting from the detection process (a); zooming one polygon (b); the same polygon overlapped with the 2009 Goggle Earth satellite image (c) and from 2010, respectively (d)</figcaption></figure>

## 5. Conclusions

The tree wind throws caused by storms are much easier to detect by the authorities if they are contiguous, compared to the diffuse or isolated ones. The diffuse wind throws occur in the temperate zone at wind gusts of 100-140 km/h and their areas range between a few tens of square meters and a few hectares, but they are very numerous, of several hundreds in the wooded mountain massifs. For this reason, they are more difficult to detect, map and, obviously, it is much more difficult to characterize the extent of the damage. However, the detection of changes occurred as a result of the diffuse wind throws, even if they are several hundred, is also important in order to set the spatial and temporal rehabilitation and forest management trends. In this regard, this paper presents a rapid procedure for the detection of the diffuse wind throws based on low cost images, knowing that certain sensors cannot be accessed without significant investments. Our application is based on the study of effects caused by the Xynthia storm in the Vosges Mountains, on 28 February 2010.

The results of this research highlight the usefulness of the Landsat satellite imagery, more easily to access than other types of images for the analysis and detection of wind throws. The detection of areas damaged by the Xynthia storm by applying the pixels delineation method based on the histogram enabled us to automatically delimit the forest pixels from the non-forest ones, thus reducing the working time. The Disturbance Index enabled the more accurate delineation of the affected areas.

538

539

540

541

542

543

544

545

546

547

548

549

550

551

552

553

554

555

556

557

558

559
The classification of the area into forest and non-forest increases the precision level for the detection of affected areas. The Tasseled Cap Transformation and the calculation of the Disturbance Index enable us to highlight the affected areas. The accuracy assessment of the results obtained for the change map was high (86 % overall accuracy and the value of the Kappa coefficient was 0.84). Also, the constant forest class in the two images, before and after the phenomenon, had the user's and producer's accuracy higher than the disturbance class, i.e. 89.10% and 81% respectively, compared to 78% and 80%.

The spatial analysis enabled us to highlight the most affected areas, resulting in 257 affected polygons, totalling 229 ha. The largest polygon is 13 ha, located east from the Grand-Soldat hamlet, Abreschviller village, department of Moselle. Moreover, following the analysis conducted, it was concluded that most of the affected polygons had a predominantly coniferous vegetation composition (75%), with slopes facing N, NW and W. For validation purposes, high-resolution images and orthophotoplans taken before and after the storm were used. The error matrix was calculated, achieving an overall accuracy of 86%, which reveals that our methodology adapted for low cost image resources, can be used for the rapid detection of diffuse wind throws in the temperate zone.

### Supplementary Materials:

**Table S1.** Polygons obtained after running the detection algorithm.

**Acknowledgments:** This research was supported by Centre de recherche en Géographie, LOTTER-EA7304, Université de Lorraine, Metz. The second author mentions that his contribution is a product of the activity undertaken in 2015 at the Laboratory LOTERR-EA7304 during the Erasmus scholarship mobility.

**Author Contributions:** I. Haidu and P. Furtuna conceived and designed the experiments; P. Furtuna performed the experiments; I. Haidu and S. Lébaut analyzed the data; P. Furtuna, I. Haidu and S. Lébaut wrote the paper.

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

## **584 Appendix A**

The largest polygon has an area of 13 ha (450-550 m altitude), East of the Grand-Soldat hamlet, Abreschviller village, department of Moselle. In the neighbouring Southern village of Saint-Quirin, the second disturbed polygon is the size of 9.2 ha (450-330 m altitude). Also in Moselle there is a 7.7 ha polygon in the Northern Vosges Regional Natural Park, North of Baerenthal and a 7.9 hectare polygon South-West of Mouterhouse. In the department of Vosges there is a 6.1 hectare polygon South of Arrentès-de-Corcieux. In the department of Haut-Rhin in Alsace there is a 7.5 ha disturbed area North-West of Lake Grand Ballon, on the territory of the Lauteenbach-Zell village.

Following the conducted analysis, it was found that most of the affected polygons had predominantly coniferous vegetation (75%), at an average altitude of 950 m, with N-, NW- and W-facing slopes. Of the 257 polygons, only 7 have more than 5 hectares of disturbed area, but they are 51.12 % of the total wind throw area. 56 polygons range from 1 to 5 ha and the remaining 194 polygons are under 1 ha. The smallest affected polygon detected by the methodology presented above and based on the limited data images approach is 0.09 ha. This is much better than the detection threshold based on 1 ha of satellite imagery that was used to assess the damages caused by the 1999 storms [83].

The large number of small polygons detected, 194 polygons under 1 ha, come as a result of the particularities of interaction between high wind speeds, over 130 km/h and the relief features of the Vosges Mountains. Generally, the affected polygons appear on the slopes facing the direction of the Xynthia storm. At the escalation of the slopes, the speed of the wind gusts accelerated resulting in an effect similar to bearing, resulting in a greater pressure exerted on the forest cover area.

## **604 References**

1. Rauch, E. Wind and storms. In: Pfister, C., Schellnhuber, H., Rahmstorf, S., Graßl, H. (Eds.), Weather catastrophies and climate change. Is there still hope for us? MunichRe, Munich, 2005, pp. 132–143.
607 2. Bernardes, S.; Madden, M. (2016) Vegetation disturbance and recovery following a rare windthrow event in the Great Smoky Mountains National Park, Int. Arch. *Photogramm. Remote Sens. Spatial Inf. Sci.*, XLI-B8, **2016**, pp.571-575. DOI:10.5194/isprs-archives-XLI-B8-571-2016

608 3. Einzmann K.; Immitzer M.; Böck S.; Bauer O.; Schmitt, A.; Atzberger, A. Windthrow Detection in European Forests with Very High-Resolution Optical Data. Forests, **2017**, 8, 21, DOI10.3390/f8010021.

610 4. Vernier, F. L'ouragan Lothar, 26 décembre 1999 effets et perspectives sur la forêt Lorraine. Bulletin de l'Académie Lorraine des Sciences, **2001**, 40 (1-2), 40-50.

612 5. Birot, Y. ; Landmann, G.; Bonhême, I. *La forêt face aux tempêtes*. Editeur: Quae, France **2009**, ISBN 978-2-7592-0330-7

616 6. Inventaire Forestier National - IFN Tempête Klaus du 24 janvier 2009: 234 000 hectares de forêt affectés à plus de 40 % 42,5 millions de mètres cubes de dégât. L'IF, **2009**, nr. 21, 1er trimestre. Available online http://inventaire-forestier.ign.fr/spip/IMG/pdf/IF21_internet.pdf.

619 7. Pôle interrégional de la Santé des Forêts du Nord-Est. *Bilan de la santé des forets 2010 en Alsace*. DRAAF/SRAL, Département de la santé des forêts, Metz, **2011a**.

621 8. Pôle interrégional de la Santé des Forêts du Nord-Est. *Bilan de la santé des forets 2010 en Lorraine*. DRAAF/SRAL, Département de la santé des forêts, Metz, **2011b**.

623 9. Anziani, A. Rapport d'information sur les conséquences de la tempête Xynthia. Session ordinaire du Sénat, **2010**, No. 554/10 juin, 100p.

625 10. Hernu, H.; Kbaier, R.; Rochard, J.; Casteigts, M.; Furst, X.; Jullien, B.; Balay, L.P.; Guillet, M. Rapport sur l'évaluation des dommages causés par la tempête xvnthia des 27 et 28 février 2010 à prendre en compte a titre du fonds de s olidarité de l'union européenne Ministère de l'écologie, de l'énergie, du développement durable et de la mer, **2010**.

629 11. FFSA/GEMA, Fédération française des sociétés d'assurances / Groupement des entreprises mutuelles d'assurance. La tempête Xynthia du 28 février 2010 – Bilan chiffré au 31 décembre 2010. **2011**, 21 pp.

631 12. Calvet, M.M.; Manable, C. Rapport d'information sur *Xynthia 5 ans après: pour une véritable culture du risque dans les territoires*. Session ordinaire du Sénat, **2015**, No. 536/18 juin.

633 13. Office national des forêts - *Tempête Xynthia: des dégâts diffus en forêt, 2010*. Available online: http://www.foretpriveefrancaise.com/actualite/voir/515/page:24. Accessed: 12 July, 2016.

635 14. Forêt Privée Française (FPF). Available online http://www.foretpriveefrancaise.com/tempete-xynthia-des-degats-localise

637 15. Papaik, M.J; Canham, C.D. Species resistance and community response to wind disturbance regimes in northern temperate forests. *J Ecol*, **2006**, 94:1011-1026.

639 16. Bonan, G.B. Forests and climate change: forcings, feedbacks, and the climate benefits of forests. *Science*, **2008**, 320(5882) pp.1444-9. DOI: 10.1126/science.1155121.

641 17. Klaus, M.; Holsten, A.; Hostert, P.; Kropp, J. P. (2011) Integrated methodology to assess windthrow impacts on forest stands under climate change. *Forest Ecology and Management*, **2011**, 261, pp.1799-1810.

643 18. Merry, K.; Bettinger, P.; Hepinstall-Cymerman, J. A prototype model for estimating the location of forest damage from tropical cyclones. *Geographia Technica*, **2011**, 6, 2, pp.65-80.

645 19. Thom, D.; Rammer, W.; Dirnböck, T.; Müller, J.; Kobler, J.; Katzensteiner, K.; Helm, N.; Seidl, R. The impacts of climate change and disturbance on spatio-temporal trajectories of biodiversity in a temperate forest landscape. *J Appl Ecol*, **2017**, 54, pp. 28-38.

648 20. Seidl, R.; Schelhaas, M.J.; Rammer, W.; Verkerk, P.J. Increasing forest disturbances in Europe and their impact on carbon storage. *Nat Clim Chang*, **2014**, 4 (9), p. 806-810

650 21. Merry, K.; Bettinger, P.; Siry, J.; Bowker, J.M. (2015) Estimating urban forest carbon sequestration potential in the southern United States using current remote sensing imagery sources. *Geographia Technica*, **2015**, 10, 2, pp.78-89.

653 22. Huang, C.; Wylie, B.; Yang, L.; Homer, C.; Zylstra, G. Derivation of a tasselled cap transformation based on Landsat 7 at-satellite reflectance. *Int. J. Remote Sensing*, **2002**, 23, 8, pp.1741-1748. DOI: 10.1080/01431160110106113

656 23. Healey, S.P.; Cohen, W.B.; Zhiqiang, Y.; Krankina, O.N. Comparison of Tasseled Cap-based Landsat data structures for use in forest disturbance detection. *Remote Sens. Environ.*, **2005**, 97, pp.301-310. DOI: 10.1016/j.rse.2005.05.009
659 24. Baumann, M.; Ozdogan, M.; Wolter, P.T.; Krylov, A.; Vladimirova, N.; Radeloff, V.C. (2014), Landsat remote sensing of forest windfall disturbance. *Remote Sens. Environ.*, **2014**, *143*, pp.171–179. DOI: 10.1016/j.rse.2013.12.020.

660
661
662
663
664
25. Jin, S.; Sader, S.A. Comparison of time-series tasseled cap wetness and the normalized difference moisture index in detecting forest disturbances. *Remote Sens. Environ.* **2005**, *94*(3) pp.364-372. DOI: 10.1016/j.rse.2004.10.012

665
666
667
668
669
26. Huang, C.; Goward, S.N.; Masek, J.G.; Thomas, N.; Zhu, Z.; Vogelmann, J.E. An automated approach for reconstructing recent forest disturbance history using dense Landsat time series stacks. *Remote Sens. Environ.*, **2010**, *114*, pp.183–198. DOI: 10.1016/j.rse.2007.07.023

670
671
672
673
674
27. Zhu, Z.; Woodcock, C.E.; Olofsson, P. Continuous monitoring of forest disturbance using all available Landsat imagery. *Remote Sens. Environ.*, **2012**, *122*, pp.75–91.

675
676
677
678
679
680
28. Lambert, J.; Drenou, C.; Denux, J.P.; Balent, G.; Cheret, V. (2013) Monitoring forest decline through remote sensing time series analysis. *GIScience & Remote Sensing*, 2013, 50 (4), pp. 437-457.

681
682
683
684
685
29. Stach, N.; Deshayes, M.; Le Toan, T. Mapping storm damage to forests using optical and radar remote sensing – The case of the December 1999 storms in France, in: Olsson H. (Ed.), Proceedings of ForestSat, 2005, Borås May 31-Jun

690
691
692
693
694
695
30. Huang, C.Q.; Song, K.; Kim, S.; Townshend, J.R.G., Davis, P., Masek, J.G., Goward, S.N. Use of a dark object concept and support vector machines to automate forest cover change analysis. *Remote Sens. Environ.*, **2008**, *112*, pp.970–985. DOI: 10.1016/j.rse.2007.07.023

698
699
700
701
702
703
31. Kennedy, R.E.; Yang, Z.; Cohen, E.; Braaten, J.; Nelson, P., Spatial and temporal patterns of forest disturbance and regrowth within the area of the Northwest Forest Plan. *Remote Sens. Environ.* **2012**, *122*, 117–133

705
706
707
708
709
710
32. Stach, N.; Deshayes, M.; Le Toan, T. Évaluation des dégâts de tempête par télédétection satellitaire: Aspects méthodologiques et opérationnels. Rapport final, Convention DERF Nr. 2001-12-132, 2002, 148.

711
712
713
714
715
716
33. Inventaire Forestier National - IFN (2008), Cartographie des coupes rases des peuplements de pin maritime en Aquitaine par détection satellitaire. Comité interprofessionnel des bois d'Aquitaine (CIBA) et Inventaire Forestier National (IFN), février 2008.

717
718
719
720
721
722
34. Chéret, V. ; Denux, J.P. ; Gacherieu, C. ; Ortisset, J.P. Utilisation de séries temporelles d'images satellites pour cartographier le dépérissement des boisements résineux du Sud Massif Central. *Rendez-vous Techniques de l'ONF, 2010*, (31), pp.55-62. http://prodinra.inra.fr/record/179412

723
724
725
726
727
728
35. Hussain, M.; Chen, D.; Cheng, A.; Wei, H.; Stanley, D. Change detection from remotely sensed images: From pixel-based to object-based approaches, *IPRS Journal of Photogrammetry and Remote Sensing*, **2013**, *80*, 91-106.

729
730
731
732
733
734
36. Sélection Vosges. Available online. http://www.selection-vosges.com/fr/la-foret-vosgienne.html.

735
736
737
738
739
740
741
Garnier, E. Les forêts vosgiennes à l'épreuve des tempêtes sous l'Ancien Régime (XVIIe et XVIIIe siècles). *Revue Gògraphique de l'Est*, **2000**, *40* (3), pp. 1-12.

745
746
747
748
749
750
751
Doll, D. Statistiques historiques des grands chablis éoliens en Europe occidentale depuis le milieu du XIXe siècle: analyse critique. In: Drouineau S., Laroussinie O., Birot Y., Terrasson D., Formery T., Roman-Amat B. *Expertise collective sur les tempêtes, la sensibilité des forêts et sur leur reconstitution*. Dossier de l'environnement de l'INRA n°20, INRA-ME&S, Paris, 2000, pp. 38-41

755
756
757
758
759
760
39. Meteo France. Available online: http://tempetes.meteo.fr/Classification-des-tempetes.html

765
766
767
768
769
770
771
772
773
774
775
776
777
778
40. Bedacht, E.; Hofherr, T. *Winter Storm Xynthia in southwest Europe and Germany*, MUNICH RE Topics Geo, 2010, pp.18-21.

780
781
782
783
784
785
41. Doll, D. Les cataclysmes éoliens dans les forêts d'Europe, aperçu historique. Forêt Entreprise, **1991**, 77(5).

789
790
791
792
793
794
795
796
797
798
799
800
801
42. Dreveton, C. L'évolution du nombre de tempêtes en France sur la période 1950-1999. *La Météorologie*, **2002**, *37*, pp. 46-56.

805
806
807
808
809
810
43. Bock, J. ; Vinkler, I. ; Duplat, P. ; Renaud, J.P. Stabilité au vent des hêtraies: les enseignements de la tempête de 1999. Dossier Facteurs de stabilité des peuplements. Forêt Entreprise, **2004**, 156, pp. 34-38.

815
816
817
818
819
820
44. Schwarz, M.; Steinmeier, C.; Holecz, F.; Stebler, O.; Wagner, H. Detection of Windthrow in Mountainous Regions with Different Remote Sensing Data and Classification Methods. *Scandinavian Journal of Forest Research*, **2003**, 18:6, pp.525 – 536.

825
826
827
828
829
830
45. Boyle, S.A.; Kennedy, C.M.; Torres, J.; Colman, K., Perez-Estigarribia, P.E.; Sancha, N.U. High-Resolution Satellite Imagery Is an Important yet Underutilized Resource in Conservation Biology. *PLoS ONE*, **2014**, *9*(1). DOI:10.1371/journal.pone.0086908.
713 46. Goward, S.; Arvidson, T.; Williams, D.; Faundeen, J.; Irons, J.; Franks, F. Historical Record of Landsat Global Coverage: Mission Operations, NSLRSDA, and International Cooperator Stations. *Photogrammetric Engineering & Remote Sensing*, **2006**, 72, 10, 1155-1169.

714
715
716 47. Vogelmann, J. E., Gallant, A. L., Shi, H. and Zhu, Z. (2016) Perspectives on monitoring gradual change across the continuity of Landsat sensors using time-series data. *Remote Sens. Environ.*, 185, pp. 258–270

718 48. Turner, W.; Rondinini, C.; Pettorelli, N.; Mora, B.; Leidner, A.K.; Szantoi, Z.; Buchanan, G.; Dech, S.; Dwyer, J.; Herold, M.; Koh, L.P.; Leimgruber, P.; Taubenboeck, H.; Wegmann, M.; Wikelski, M.; Woodcock, C. Free and open-access satellite data are key to biodiversity conservation. *Biological Conservation*, **2015**, 182, pp. 173–176.

722 49. Cohen, W.; and Goward, S. Landsat's Role in Ecological Applications of Remote Sensing. *BioScience* **2004**, 54(6), pp535-545. DOI: 10.1641/0006-3568(2004)054[0535:LRIEAO]2.0.CO;2.

724 50. Mukai, Y.; Hasegawa, I. Extraction of damaged areas of windfall trees by typhoons using Landsat TM data. *International Journal of Remote Sensing*, **2000**, 21, pp.647-654;

726 51. Ramsey, E.W.; Hodgson, M.E.; Sapkota, S.K.; Nelson, G.A. Forest impact estimated with NOAA AVHRR and LandsatTM data related to an empirical hurricane wind-field distribution. *Remote Sensing Environment*, **2001**, 77, 279 – 292;

729 52. Jonikavicius, D.; Mozgeris, G.; Rapid assessment of wind storm-caused forest damage using satellite images and stand-wise forest inventory data. *iForest*, **2013**, 6, pp.50-55;

731 53. Hermosilla, T.; Wulder, M.A.; White, J.C.; Coops, N.C.; Hobart, G.W. Regional detection, characterization, and attribution of annual forest change from 1984 to 2012 using Landsat-derived time-series metrics. *Remote Sensing Environment*, **2015**, 170, pp. 121–132;

734 54. Jackson, R.G.; Foody, G.M.; Quine, C.P. Characterising windthrown gaps from fine spatial resolution remotely sensed data. *Forest Ecology and Management*, **2000**, 135, pp.253 - 260.

736 55. Malthus, T.J.; Suárez, J.C.; Woodhouse, I.H.; Shaw, D.T. Review of remote sensing in commercial forestry. Forestry Commission. Forest research internal report, 2002.

738 56. Suárez, J.C.; Ontiveros, C.; Smith, S.; Snape, S. Use of airborne LiDAR and aerial photography in the estimation of individual tree heights in forestry. *Computers and Geosciences*, **2005**, 31(2), pp.253-262.

740 57. Melesse, A.M.; Weng, Q.; Thenkabail, P.S.; Senay, G.B. Remote Sensing Sensors and Applications in Environmental Resources Mapping and Modelling. *Sensors*, **2007**, 7, pp. 3209-3241.

742 58. European Space Agency. Available online: https://earth.esa.int

743 59. Land info. Available online: http://www.landinfo.com/

744 60. Gancz, V.; Apostol, B.; Petrila, M.; Lorent, A. Detectarea cu ajutorul imaginilor satelitare a doborâturilor de vânt și evaluarea efectelor acestora. *Rev. pădur*. **2010**, Nov. 01. Available online: http://www.revistapadurilor.ro/(16713) accessed on 2016-05-19.

747 61. Zhu, Z.; Woodcock, C.E. Object-based cloud and cloud shadow detection in Landsat imagery. *Remote Sensing of the Environment*. **2012**, 118, pp.83-94.

749 62. Goodwin, N.R.; Collett, L.J.; Denham, R.J.; Flood, N.; Tindall, D. Cloud and cloud shadow screening across Queensland, Australia: An automated method for Landsat TM/ETM+ time series. *Remote Sensing of Environment*, **2013**, 134, pp. 50–65. DOI:10.1016/j.rse.2013.02.019

752 63. Glovis. Available online: http://glovis.usgs.gov/

753 64. Duché, Y.; Savazzi, R.; Reymond, B., Carnus, E., Crespy, L. Cartographie de la sensibilité de la végétation aux incendies de forêts en région méditerranéenne: apports des images satellites et limites d'usage. *RDV techniques*, **2011**, nr.32 - printemps 2011 - ONF.

756 65. Chander, G.; Markham, B.L. Revised Landsat-5 TM Radiometric Calibration Procedures and Post calibration Dynamic Ranges. *IEEE Transactions on Geosciences and Remote Sensing*, **2003**, 41, 11, pp.2674-2677.

758 66. Chander, G.; Markham, B.L.; Helder, D.L. Summary of current radiometric calibration coefficients for Landsat MSS, TM, ETM+, and EO-1 ALI sensors, *Remote Sens. Environ.*, **2009**, 113, 893–903.

760 67. Zhu, Z.; Woodcock, C.E. Continuous change detection and classification of land cover using all available Landsat data. *Remote Sens. Environ.* **2014**, 144, pp.152-171.

762 68. Martinuzzi, S.; Gould, W.A.; Ramos-Gonzalez, O.M. Creating Cloud-Free Landsat ETM+ Data Sets in Tropical Landscapes: Cloud and Cloud-Shadow Removal. General Technical Report, IITF-GTR-32. USDA, Forest Service, International Institute of Tropical Forestry. 2007.

765 69. Kennedy, R.E.; Cohen, W.B.; Moisen, G.G.; Goward, S.N.; Wulder, M.; Powell, S.L.; Masek, J.G.; Huang, C.; Healey, S P. A sample design for Landsat-based estimation of national trends in forest disturbance and
767 regrowth. NASA joint workshop on biodiversity, terrestrial ecology, and related applied sciences. August 21–25, 2006, College Park, MD.

769 70. Kuemmerle, T.; Hostert, P.; Radeloff, V.C.; Perzanowski, K.,; Kruhlov, I. Post- socialist forest disturbance in the Carpathian border region of Poland, Slovakia, and Ukraine. *Ecol. Appl.* **2007**, 17, pp.1279–1295.

771 71. Masek, J.G.; Huang, C.; Wolfe, R.; Cohen, W.; Hall, F.; Kutler, J.; Nelson, P. North American forest disturbance mapped from decadal Landsat record. *Remote Sens. Environ.* **2008**,112, pp.2914–2926. DOI: 10.1016/j.rse.2008.02.010

774 72. Kuemmerle, T.; Chaskovskyy, O.; Knorn, J.; Radeloff, V.C.; Kruhlov, I.; Keeton, W.S.; Hostert, P. Forest cover change and illegal logging in the Ukrainian Carpathians in the transition period from 1988 to 2007. *Remote Sens. Environ.* **2009**,113, pp.1194–1207.

777 73. Neigh, Christopher S R, Bolton, Douglas K., Diabate, Mouhamad, Williams, Jennifer J. Carvalhais, Nuno. An automated approach to map the history of forest disturbance from insect mortality and harvest with landsat time-series data. Remote Sens. 2014, 6, 2782-2808; doi:10.3390/rs6042782

780 74. Frantz, D.; Röder, A.; Udelhoven, T.; Schmidt, M. Forest disturbance mapping using dense synthetic Landsat/MODIS time-series and permutation-based disturbance index detection. *Remote Sensing*, 2016, 8(4). DOI: 10.3390/rs8040277.

783 75. Furtuna, P. Temporal and spatial variation of forest coverage in apuseni natural park, 2000-2014 period. *Geographia Technica*, 2017, 12(1). DOI: 10.21163/GT_2017.121.05

785 76. Darwish, A.; Leukert, K.; Reinhardt, W. Image Segmentation for The Purpose of Object-Based Classification, Geoscience and Remote Sensing Symposium. IGARSS '03. Proceedings. *IEEE International*. **2003**, 3 pp.2039 – 2041. DOI:10.1109/IGARSS.2003.1294332.

788 77. Pekkarinen, A., Image segmentation in multi-source forest inventory. Finnish Forest Research Institute. Research papers 926, 2004, ISBN 952-10-2003-2.

790 78. Yoshino, K.; Kawaguchi, S.; Kanda, F., Tsai, F. Very High Resolution Plant Community Mapping at High Moor, Kushiro Wetland. *Photographic engineering and remote sensing*, **2014**, 80 (9), pp.895-905.

792 79. Geoportail. Available online: http://www.geoportail.gouv.fr/

793 80. Institut national de l'information géographique et forestière. Available online: http://professionnels.ign.fr

794 81. Inventaire forestiere. Available online: http://inventaire-forestier.ign.fr/spip/

795 82. Carmen - L'application cartographique au service des données environnementales. Available online http://carmen.carmencarto.fr/

797 83. Inventaire Forestier National - IFN (2003), Les tempêtes de décembre 1999 bilan national et Enseignements, L'IF n° 2, décembre 2003.

799

800 Supplementary Materials:

801 Table S1. Polygons obtained after running the detection algorithm.
<table><thead><tr><th></th><th>OBJECTID</th><th>Id</th><th>Shap_Lg</th><th>Area_m²</th><th>Ara_ha</th><th>X (decimal degrees)</th><th>Y (decimal degrees)</th><th>Longitude</th><th>Latitude</th></tr></thead><tbody><tr><td>1</td><td>381</td><td>381</td><td>327.66</td><td>4921.4</td><td>0.4921</td><td>7.347278</td><td>49.147080</td><td>7° 20' 50.202"</td><td>49° 8' 49.486"</td></tr><tr><td>2</td><td>2461</td><td>2461</td><td>269.30</td><td>3290.6</td><td>0.3291</td><td>7.323385</td><td>49.070599</td><td>7° 19' 24.184"</td><td>49° 4' 14.156"</td></tr><tr><td>3</td><td>7037</td><td>7037</td><td>1454.88</td><td>77313.3</td><td>7.7313</td><td>7.509725</td><td>48.991177</td><td>7° 30' 35.010"</td><td>48° 59' 28.237"</td></tr><tr><td>4</td><td>7462</td><td>7462</td><td>169.76</td><td>1228.8</td><td>0.1229</td><td>7.515414</td><td>48.989874</td><td>7° 30' 55.490"</td><td>48° 59' 23.547"</td></tr><tr><td>5</td><td>7482</td><td>7482</td><td>120.00</td><td>900.0</td><td>0.0900</td><td>7.514831</td><td>48.989579</td><td>7° 30' 53.391"</td><td>48° 59' 22.484"</td></tr><tr><td>6</td><td>7601</td><td>7601</td><td>256.60</td><td>2963.2</td><td>0.2963</td><td>7.511610</td><td>48.988011</td><td>7° 30' 41.797"</td><td>48° 59' 16.838"</td></tr><tr><td>7</td><td>8227</td><td>8227</td><td>219.57</td><td>2034.1</td><td>0.2034</td><td>7.560925</td><td>48.979285</td><td>7° 33' 39.328"</td><td>48° 58' 45.428"</td></tr><tr><td>8</td><td>8925</td><td>8925</td><td>169.77</td><td>1229.0</td><td>0.1229</td><td>7.567033</td><td>48.970118</td><td>7° 34' 1.319"</td><td>48° 58' 12.423"</td></tr><tr><td>9</td><td>9353</td><td>9353</td><td>213.01</td><td>2132.2</td><td>0.2132</td><td>7.410459</td><td>48.961315</td><td>7° 24' 37.651"</td><td>48° 57' 40.735"</td></tr><tr><td>10</td><td>9363</td><td>9363</td><td>1748.00</td><td>78665.9</td><td>7.8666</td><td>7.407721</td><td>48.963144</td><td>7° 24' 27.794"</td><td>48° 57' 47.317"</td></tr><tr><td>11</td><td>10361</td><td>10361</td><td>224.28</td><td>2418.1</td><td>0.2418</td><td>7.188378</td><td>48.938636</td><td>7° 11' 18.162"</td><td>48° 56' 19.089"</td></tr><tr><td>12</td><td>11771</td><td>11771</td><td>489.00</td><td>10211.5</td><td>1.0211</td><td>7.404978</td><td>48.913743</td><td>7° 24' 17.922"</td><td>48° 54' 49.473"</td></tr><tr><td>13</td><td>11823</td><td>11823</td><td>432.87</td><td>8105.0</td><td>0.8105</td><td>7.400416</td><td>48.912271</td><td>7° 24' 1.499"</td><td>48° 54' 44.177"</td></tr><tr><td>14</td><td>11863</td><td>11863</td><td>224.28</td><td>2418.1</td><td>0.2418</td><td>7.396563</td><td>48.910951</td><td>7° 23' 47.627"</td><td>48° 54' 39.425"</td></tr><tr><td>15</td><td>11884</td><td>11884</td><td>224.28</td><td>2418.1</td><td>0.2418</td><td>7.394142</td><td>48.909838</td><td>7° 23' 38.911"</td><td>48° 54' 35.418"</td></tr><tr><td>16</td><td>11936</td><td>11936</td><td>213.01</td><td>2132.2</td><td>0.2132</td><td>7.388374</td><td>48.908921</td><td>7° 23' 18.148"</td><td>48° 54' 32.115"</td></tr><tr><td>17</td><td>11937</td><td>11937</td><td>362.36</td><td>7525.5</td><td>0.7526</td><td>7.392273</td><td>48.909277</td><td>7° 23' 32.181"</td><td>48° 54' 33.398"</td></tr><tr><td>18</td><td>11975</td><td>11975</td><td>120.00</td><td>900.0</td><td>0.0900</td><td>7.395239</td><td>48.908091</td><td>7° 23' 42.861"</td><td>48° 54' 29.127"</td></tr><tr><td>19</td><td>11992</td><td>11992</td><td>172.66</td><td>1046.0</td><td>0.1046</td><td>7.394726</td><td>48.907814</td><td>7° 23' 41.014"</td><td>48° 54' 28.130"</td></tr><tr><td>20</td><td>12014</td><td>12014</td><td>322.64</td><td>4991.2</td><td>0.4991</td><td>7.388727</td><td>48.907413</td><td>7° 23' 19.416"</td><td>48° 54' 26.686"</td></tr><tr><td>21</td><td>12063</td><td>12063</td><td>268.59</td><td>3297.5</td><td>0.3298</td><td>7.387090</td><td>48.906768</td><td>7° 23' 13.525"</td><td>48° 54' 24.364"</td></tr><tr><td>22</td><td>12084</td><td>12084</td><td>170.22</td><td>1214.9</td><td>0.1215</td><td>7.385695</td><td>48.906354</td><td>7° 23' 8.502"</td><td>48° 54' 22.875"</td></tr><tr><td>23</td><td>12205</td><td>12205</td><td>180.00</td><td>1800.0</td><td>0.1800</td><td>7.382662</td><td>48.904542</td><td>7° 22' 57.582"</td><td>48° 54' 16.350"</td></tr><tr><td>24</td><td>12453</td><td>12453</td><td>120.00</td><td>900.0</td><td>0.0900</td><td>7.216358</td><td>48.895476</td><td>7° 12' 58.889"</td><td>48° 53' 43.715"</td></tr><tr><td>25</td><td>13489</td><td>13489</td><td>169.77</td><td>1229.0</td><td>0.1229</td><td>7.284117</td><td>48.852362</td><td>7° 17' 2.823"</td><td>48° 51' 8.504"</td></tr><tr><td>26</td><td>13508</td><td>13508</td><td>322.66</td><td>4274.5</td><td>0.4275</td><td>7.282140</td><td>48.851920</td><td>7° 16' 55.705"</td><td>48° 51' 6.912"</td></tr></tbody></table>
<table><tr><td>27</td><td>13525</td><td>13525</td><td>376.43</td><td>5786.9</td><td>0.5787</td><td>7.285366</td><td>48.850807</td><td>7° 17' 7.317"</td><td>48° 51' 2.905"</td></tr><tr><td>28</td><td>13540</td><td>13540</td><td>267.12</td><td>2878.8</td><td>0.2879</td><td>7.281793</td><td>48.849646</td><td>7° 16' 54.453"</td><td>48° 50' 58.727"</td></tr><tr><td>29</td><td>13556</td><td>13556</td><td>120.00</td><td>900.0</td><td>0.0900</td><td>7.279782</td><td>48.848400</td><td>7° 16' 47.215"</td><td>48° 50' 54.239"</td></tr><tr><td>30</td><td>13559</td><td>13559</td><td>120.00</td><td>900.0</td><td>0.0900</td><td>7.280200</td><td>48.848136</td><td>7° 16' 48.719"</td><td>48° 50' 53.289"</td></tr><tr><td>31</td><td>13568</td><td>13568</td><td>120.00</td><td>900.0</td><td>0.0900</td><td>7.280618</td><td>48.847872</td><td>7° 16' 50.224"</td><td>48° 50' 52.340"</td></tr><tr><td>32</td><td>14882</td><td>14882</td><td>120.00</td><td>900.0</td><td>0.0900</td><td>7.275157</td><td>48.768712</td><td>7° 16' 30.565"</td><td>48° 46' 7.363"</td></tr><tr><td>33</td><td>15198</td><td>15198</td><td>265.37</td><td>3187.2</td><td>0.3187</td><td>7.349408</td><td>48.748613</td><td>7° 20' 57.868"</td><td>48° 44' 55.007"</td></tr><tr><td>34</td><td>15277</td><td>15277</td><td>209.69</td><td>1972.1</td><td>0.1972</td><td>7.281863</td><td>48.743617</td><td>7° 16' 54.708"</td><td>48° 44' 37.022"</td></tr><tr><td>35</td><td>15289</td><td>15289</td><td>360.00</td><td>7200.0</td><td>0.7200</td><td>7.282939</td><td>48.740911</td><td>7° 16' 58.580"</td><td>48° 44' 27.279"</td></tr><tr><td>36</td><td>15298</td><td>15298</td><td>334.52</td><td>6612.2</td><td>0.6612</td><td>7.282392</td><td>48.742647</td><td>7° 16' 56.610"</td><td>48° 44' 33.527"</td></tr><tr><td>37</td><td>16608</td><td>16608</td><td>833.26</td><td>32473.6</td><td>3.2474</td><td>7.140469</td><td>48.687928</td><td>7° 8' 25.688"</td><td>48° 41' 16.542"</td></tr><tr><td>38</td><td>16672</td><td>16672</td><td>718.37</td><td>18109.8</td><td>1.8110</td><td>7.142110</td><td>48.685686</td><td>7° 8' 31.594"</td><td>48° 41' 8.470"</td></tr><tr><td>39</td><td>16674</td><td>16674</td><td>299.45</td><td>4802.4</td><td>0.4802</td><td>7.143827</td><td>48.685477</td><td>7° 8' 37.779"</td><td>48° 41' 7.716"</td></tr><tr><td>40</td><td>16938</td><td>16938</td><td>120.00</td><td>900.0</td><td>0.0900</td><td>7.132053</td><td>48.674985</td><td>7° 7' 55.392"</td><td>48° 40' 29.944"</td></tr><tr><td>41</td><td>16939</td><td>16939</td><td>264.55</td><td>4093.0</td><td>0.4093</td><td>7.132919</td><td>48.675171</td><td>7° 7' 58.509"</td><td>48° 40' 30.615"</td></tr><tr><td>42</td><td>16963</td><td>16963</td><td>643.55</td><td>10507.4</td><td>1.0507</td><td>7.130703</td><td>48.674397</td><td>7° 7' 50.530"</td><td>48° 40' 27.830"</td></tr><tr><td>43</td><td>16981</td><td>16981</td><td>169.77</td><td>1228.9</td><td>0.1229</td><td>7.132906</td><td>48.673229</td><td>7° 7' 58.463"</td><td>48° 40' 23.623"</td></tr><tr><td>44</td><td>17028</td><td>17028</td><td>661.69</td><td>15097.1</td><td>1.5097</td><td>7.128864</td><td>48.671886</td><td>7° 7' 43.911"</td><td>48° 40' 18.788"</td></tr><tr><td>45</td><td>17106</td><td>17106</td><td>789.71</td><td>12482.8</td><td>1.2483</td><td>7.132993</td><td>48.670633</td><td>7° 7' 58.775"</td><td>48° 40' 14.278"</td></tr><tr><td>46</td><td>17144</td><td>17144</td><td>285.23</td><td>3842.1</td><td>0.3842</td><td>7.134045</td><td>48.668716</td><td>7° 8' 2.562"</td><td>48° 40' 7.377"</td></tr><tr><td>47</td><td>19017</td><td>19017</td><td>256.11</td><td>3961.9</td><td>0.3962</td><td>7.069609</td><td>48.621487</td><td>7° 4' 10.592"</td><td>48° 37' 17.354"</td></tr><tr><td>48</td><td>19411</td><td>19411</td><td>1777.12</td><td>129798.2</td><td>12.9798</td><td>7.157301</td><td>48.618670</td><td>7° 9' 26.284"</td><td>48° 37' 7.213"</td></tr><tr><td>49</td><td>19458</td><td>19458</td><td>212.17</td><td>2112.0</td><td>0.2112</td><td>7.183347</td><td>48.615975</td><td>7° 11' 0.048"</td><td>48° 36' 57.511"</td></tr><tr><td>50</td><td>19802</td><td>19802</td><td>120.00</td><td>900.0</td><td>0.0900</td><td>6.997850</td><td>48.606063</td><td>6° 59' 52.259"</td><td>48° 36' 21.825"</td></tr><tr><td>51</td><td>20815</td><td>20815</td><td>163.00</td><td>1255.7</td><td>0.1256</td><td>7.082725</td><td>48.589829</td><td>7° 4' 57.811"</td><td>48° 35' 23.383"</td></tr><tr><td>52</td><td>20860</td><td>20860</td><td>412.02</td><td>6571.0</td><td>0.6571</td><td>7.081275</td><td>48.589519</td><td>7° 4' 52.591"</td><td>48° 35' 22.268"</td></tr><tr><td>53</td><td>20890</td><td>20890</td><td>300.00</td><td>5400.0</td><td>0.5400</td><td>7.082554</td><td>48.588611</td><td>7° 4' 57.193"</td><td>48° 35' 19.000"</td></tr></table>
<table><tr><td>54</td><td>20948</td><td>20948</td><td>163.00</td><td>1255.7</td><td>0.1256</td><td>7.083998</td><td>48.587691</td><td>7° 5' 2.393"</td><td>48° 35' 15.686"</td></tr><tr><td>55</td><td>21005</td><td>21005</td><td>120.00</td><td>900.0</td><td>0.0900</td><td>7.085278</td><td>48.586497</td><td>7° 5' 7.001"</td><td>48° 35' 11.390"</td></tr><tr><td>56</td><td>21043</td><td>21043</td><td>120.00</td><td>900.0</td><td>0.0900</td><td>7.147099</td><td>48.586969</td><td>7° 8' 49.555"</td><td>48° 35' 13.089"</td></tr><tr><td>57</td><td>21076</td><td>21076</td><td>317.56</td><td>5114.4</td><td>0.5114</td><td>7.086101</td><td>48.585656</td><td>7° 5' 9.964"</td><td>48° 35' 8.363"</td></tr><tr><td>58</td><td>21155</td><td>21155</td><td>224.98</td><td>2781.7</td><td>0.2782</td><td>7.147809</td><td>48.585483</td><td>7° 8' 52.111"</td><td>48° 35' 7.738"</td></tr><tr><td>59</td><td>21213</td><td>21213</td><td>120.00</td><td>900.0</td><td>0.0900</td><td>7.146374</td><td>48.584528</td><td>7° 8' 46.947"</td><td>48° 35' 4.302"</td></tr><tr><td>60</td><td>21226</td><td>21226</td><td>200.25</td><td>2515.0</td><td>0.2515</td><td>7.127448</td><td>48.584163</td><td>7° 7' 38.813"</td><td>48° 35' 2.988"</td></tr><tr><td>61</td><td>21241</td><td>21241</td><td>256.11</td><td>3961.9</td><td>0.3962</td><td>7.149951</td><td>48.584210</td><td>7° 8' 59.823"</td><td>48° 35' 3.156"</td></tr><tr><td>62</td><td>21267</td><td>21267</td><td>6439.64</td><td>328767.0</td><td>9.1852</td><td>7.098719</td><td>48.587042</td><td>7° 5' 55.387"</td><td>48° 35' 13.352"</td></tr><tr><td>63</td><td>21274</td><td>21274</td><td>300.00</td><td>5400.0</td><td>0.5400</td><td>7.127102</td><td>48.583137</td><td>7° 7' 37.568"</td><td>48° 34' 59.292"</td></tr><tr><td>64</td><td>21277</td><td>21277</td><td>223.00</td><td>3055.7</td><td>0.3056</td><td>7.164146</td><td>48.583599</td><td>7° 9' 50.927"</td><td>48° 35' 0.955"</td></tr><tr><td>65</td><td>21452</td><td>21452</td><td>550.10</td><td>12540.6</td><td>1.2541</td><td>7.168475</td><td>48.581431</td><td>7° 10' 6.508"</td><td>48° 34' 53.153"</td></tr><tr><td>66</td><td>22356</td><td>22356</td><td>120.00</td><td>900.0</td><td>0.0900</td><td>7.121400</td><td>48.566849</td><td>7° 7' 17.042"</td><td>48° 34' 0.658"</td></tr><tr><td>67</td><td>22443</td><td>22443</td><td>577.47</td><td>15575.5</td><td>1.5575</td><td>7.119109</td><td>48.566457</td><td>7° 7' 8.794"</td><td>48° 33' 59.245"</td></tr><tr><td>68</td><td>22551</td><td>22551</td><td>120.00</td><td>900.0</td><td>0.0900</td><td>7.334484</td><td>48.567161</td><td>7° 20' 4.143"</td><td>48° 34' 1.779"</td></tr><tr><td>69</td><td>22958</td><td>22958</td><td>430.87</td><td>9278.3</td><td>0.9278</td><td>7.035810</td><td>48.556293</td><td>7° 2' 8.916"</td><td>48° 33' 22.654"</td></tr><tr><td>70</td><td>23024</td><td>23024</td><td>480.00</td><td>10800.0</td><td>1.0800</td><td>7.035847</td><td>48.555030</td><td>7° 2' 9.050"</td><td>48° 33' 18.107"</td></tr><tr><td>71</td><td>23038</td><td>23038</td><td>169.76</td><td>1228.8</td><td>0.1229</td><td>7.031971</td><td>48.554678</td><td>7° 1' 55.096"</td><td>48° 33' 16.842"</td></tr><tr><td>72</td><td>23094</td><td>23094</td><td>954.76</td><td>33320.0</td><td>3.3320</td><td>6.978291</td><td>48.552374</td><td>6° 58' 41.846"</td><td>48° 33' 8.547"</td></tr><tr><td>73</td><td>23116</td><td>23116</td><td>823.73</td><td>26656.1</td><td>2.6656</td><td>7.029786</td><td>48.553955</td><td>7° 1' 47.230"</td><td>48° 33' 14.238"</td></tr><tr><td>74</td><td>23143</td><td>23143</td><td>349.59</td><td>5576.1</td><td>0.5576</td><td>6.975188</td><td>48.553122</td><td>6° 58' 30.676"</td><td>48° 33' 11.238"</td></tr><tr><td>75</td><td>23176</td><td>23176</td><td>163.00</td><td>1255.7</td><td>0.1256</td><td>7.078820</td><td>48.553326</td><td>7° 4' 43.752"</td><td>48° 33' 11.974"</td></tr><tr><td>76</td><td>23187</td><td>23187</td><td>224.28</td><td>2418.2</td><td>0.2418</td><td>7.077782</td><td>48.553047</td><td>7° 4' 40.016"</td><td>48° 33' 10.970"</td></tr><tr><td>77</td><td>23220</td><td>23220</td><td>300.00</td><td>5400.0</td><td>0.5400</td><td>7.075992</td><td>48.552469</td><td>7° 4' 33.572"</td><td>48° 33' 8.889"</td></tr><tr><td>78</td><td>23288</td><td>23288</td><td>354.87</td><td>6187.5</td><td>0.6187</td><td>7.079457</td><td>48.551531</td><td>7° 4' 46.044"</td><td>48° 33' 5.512"</td></tr><tr><td>79</td><td>23302</td><td>23302</td><td>163.00</td><td>1255.7</td><td>0.1256</td><td>7.082152</td><td>48.551222</td><td>7° 4' 55.747"</td><td>48° 33' 4.400"</td></tr><tr><td>80</td><td>23380</td><td>23380</td><td>168.01</td><td>1393.6</td><td>0.1394</td><td>7.095778</td><td>48.549970</td><td>7° 5' 44.801"</td><td>48° 32' 59.891"</td></tr></table>
<table><tr><td>81</td><td>23389</td><td>23389</td><td>438.23</td><td>8861.8</td><td>0.8862</td><td>7.079972</td><td>48.549914</td><td>7° 4' 47.899"</td><td>48° 32' 59.690"</td></tr><tr><td>82</td><td>23458</td><td>23458</td><td>277.89</td><td>3387.5</td><td>0.3387</td><td>7.081279</td><td>48.548692</td><td>7° 4' 52.604"</td><td>48° 32' 55.291"</td></tr><tr><td>83</td><td>23526</td><td>23526</td><td>1684.01</td><td>53471.8</td><td>5.3472</td><td>7.095052</td><td>48.550939</td><td>7° 5' 42.189"</td><td>48° 33' 3.380"</td></tr><tr><td>84</td><td>23607</td><td>23607</td><td>621.52</td><td>14512.3</td><td>1.4512</td><td>7.095847</td><td>48.546374</td><td>7° 5' 45.050"</td><td>48° 32' 46.947"</td></tr><tr><td>85</td><td>23827</td><td>23827</td><td>2436.83</td><td>166363.4</td><td>3.8650</td><td>7.002467</td><td>48.541006</td><td>7° 0' 8.881"</td><td>48° 32' 27.620"</td></tr><tr><td>86</td><td>23843</td><td>23843</td><td>522.95</td><td>13622.9</td><td>1.3623</td><td>6.994330</td><td>48.540124</td><td>6° 59' 39.587"</td><td>48° 32' 24.445"</td></tr><tr><td>87</td><td>24072</td><td>24072</td><td>653.42</td><td>16176.8</td><td>1.6177</td><td>6.987321</td><td>48.536292</td><td>6° 59' 14.355"</td><td>48° 32' 10.652"</td></tr><tr><td>88</td><td>24075</td><td>24075</td><td>227.61</td><td>1857.7</td><td>0.1858</td><td>7.015204</td><td>48.536142</td><td>7° 0' 54.734"</td><td>48° 32' 10.112"</td></tr><tr><td>89</td><td>24091</td><td>24091</td><td>224.28</td><td>2418.2</td><td>0.2418</td><td>6.945620</td><td>48.534553</td><td>6° 56' 44.233"</td><td>48° 32' 4.389"</td></tr><tr><td>90</td><td>24116</td><td>24116</td><td>168.12</td><td>1349.7</td><td>0.1350</td><td>7.018118</td><td>48.535307</td><td>7° 1' 5.227"</td><td>48° 32' 7.105"</td></tr><tr><td>91</td><td>24196</td><td>24196</td><td>266.95</td><td>3199.2</td><td>0.3199</td><td>7.018890</td><td>48.533683</td><td>7° 1' 8.003"</td><td>48° 32' 1.259"</td></tr><tr><td>92</td><td>24226</td><td>24226</td><td>328.62</td><td>5129.5</td><td>0.5130</td><td>6.951903</td><td>48.531878</td><td>6° 57' 6.850"</td><td>48° 31' 54.762"</td></tr><tr><td>93</td><td>24286</td><td>24286</td><td>978.26</td><td>42216.8</td><td>4.2217</td><td>6.960387</td><td>48.531902</td><td>6° 57' 37.393"</td><td>48° 31' 54.848"</td></tr><tr><td>94</td><td>24322</td><td>24322</td><td>170.22</td><td>1214.9</td><td>0.1215</td><td>6.949084</td><td>48.529375</td><td>6° 56' 56.704"</td><td>48° 31' 45.751"</td></tr><tr><td>95</td><td>24740</td><td>24740</td><td>257.66</td><td>3102.6</td><td>0.3103</td><td>7.081205</td><td>48.522858</td><td>7° 4' 52.337"</td><td>48° 31' 22.288"</td></tr><tr><td>96</td><td>24906</td><td>24906</td><td>276.77</td><td>3300.4</td><td>0.3300</td><td>7.079757</td><td>48.520636</td><td>7° 4' 47.125"</td><td>48° 31' 14.289"</td></tr><tr><td>97</td><td>24963</td><td>24963</td><td>120.00</td><td>900.0</td><td>0.0900</td><td>7.078473</td><td>48.519177</td><td>7° 4' 42.502"</td><td>48° 31' 9.036"</td></tr><tr><td>98</td><td>27756</td><td>27756</td><td>390.60</td><td>6298.4</td><td>0.6298</td><td>6.908963</td><td>48.476141</td><td>6° 54' 32.266"</td><td>48° 28' 34.108"</td></tr><tr><td>99</td><td>27787</td><td>27787</td><td>411.21</td><td>9407.9</td><td>0.9408</td><td>6.906792</td><td>48.475905</td><td>6° 54' 24.452"</td><td>48° 28' 33.257"</td></tr><tr><td>100</td><td>28074</td><td>28074</td><td>1883.34</td><td>85002.9</td><td>4.1185</td><td>6.906583</td><td>48.473189</td><td>6° 54' 23.699"</td><td>48° 28' 23.480"</td></tr><tr><td>101</td><td>28208</td><td>28208</td><td>286.43</td><td>4871.0</td><td>0.4871</td><td>6.904374</td><td>48.470593</td><td>6° 54' 15.746"</td><td>48° 28' 14.134"</td></tr><tr><td>102</td><td>28330</td><td>28330</td><td>521.79</td><td>10777.3</td><td>1.0777</td><td>6.907999</td><td>48.469531</td><td>6° 54' 28.797"</td><td>48° 28' 10.312"</td></tr><tr><td>103</td><td>28533</td><td>28533</td><td>328.28</td><td>4408.7</td><td>0.4409</td><td>6.944120</td><td>48.468323</td><td>6° 56' 38.832"</td><td>48° 28' 5.964"</td></tr><tr><td>104</td><td>28662</td><td>28662</td><td>306.17</td><td>4614.9</td><td>0.4615</td><td>6.945653</td><td>48.467190</td><td>6° 56' 44.352"</td><td>48° 28' 1.885"</td></tr><tr><td>105</td><td>28837</td><td>28837</td><td>270.85</td><td>3900.6</td><td>0.3901</td><td>6.976134</td><td>48.465983</td><td>6° 58' 34.082"</td><td>48° 27' 57.538"</td></tr><tr><td>106</td><td>28967</td><td>28967</td><td>302.67</td><td>5076.5</td><td>0.5077</td><td>6.973708</td><td>48.464862</td><td>6° 58' 25.349"</td><td>48° 27' 53.503"</td></tr><tr><td>107</td><td>28969</td><td>28969</td><td>300.00</td><td>5400.0</td><td>0.5400</td><td>6.980086</td><td>48.464589</td><td>6° 58' 48.309"</td><td>48° 27' 52.520"</td></tr></table>
<table>
  <tr>
    <td>108</td>
    <td>29170</td>
    <td>29170</td>
    <td>485.16</td>
    <td>14181.3</td>
    <td>1.4181</td>
    <td>6.978799</td>
    <td>48.462471</td>
    <td>6° 58' 43.676"</td>
    <td>48° 27' 44.894"</td>
  </tr>
  <tr>
    <td>109</td>
    <td>29222</td>
    <td>29222</td>
    <td>684.90</td>
    <td>22449.5</td>
    <td>2.2450</td>
    <td>6.972086</td>
    <td>48.462778</td>
    <td>6° 58' 19.511"</td>
    <td>48° 27' 46.001"</td>
  </tr>
  <tr>
    <td>110</td>
    <td>29387</td>
    <td>29387</td>
    <td>324.64</td>
    <td>5958.0</td>
    <td>0.5958</td>
    <td>6.977524</td>
    <td>48.460374</td>
    <td>6° 58' 39.088"</td>
    <td>48° 27' 37.347"</td>
  </tr>
  <tr>
    <td>111</td>
    <td>29522</td>
    <td>29522</td>
    <td>172.66</td>
    <td>1046.0</td>
    <td>0.1046</td>
    <td>6.975760</td>
    <td>48.458542</td>
    <td>6° 58' 32.736"</td>
    <td>48° 27' 30.751"</td>
  </tr>
  <tr>
    <td>112</td>
    <td>31127</td>
    <td>31127</td>
    <td>169.76</td>
    <td>1228.8</td>
    <td>0.1229</td>
    <td>6.950751</td>
    <td>48.441517</td>
    <td>6° 57' 2.705"</td>
    <td>48° 26' 29.461"</td>
  </tr>
  <tr>
    <td>113</td>
    <td>32858</td>
    <td>32858</td>
    <td>360.00</td>
    <td>7200.0</td>
    <td>0.7200</td>
    <td>6.910702</td>
    <td>48.424617</td>
    <td>6° 54' 38.526"</td>
    <td>48° 25' 28.621"</td>
  </tr>
  <tr>
    <td>114</td>
    <td>33039</td>
    <td>33039</td>
    <td>501.15</td>
    <td>10310.9</td>
    <td>1.0311</td>
    <td>6.912671</td>
    <td>48.423554</td>
    <td>6° 54' 45.617"</td>
    <td>48° 25' 24.794"</td>
  </tr>
  <tr>
    <td>115</td>
    <td>33117</td>
    <td>33117</td>
    <td>406.52</td>
    <td>5990.2</td>
    <td>0.5990</td>
    <td>6.915411</td>
    <td>48.422283</td>
    <td>6° 54' 55.480"</td>
    <td>48° 25' 20.219"</td>
  </tr>
  <tr>
    <td>116</td>
    <td>33207</td>
    <td>33207</td>
    <td>163.00</td>
    <td>1255.7</td>
    <td>0.1256</td>
    <td>6.925611</td>
    <td>48.421648</td>
    <td>6° 55' 32.198"</td>
    <td>48° 25' 17.932"</td>
  </tr>
  <tr>
    <td>117</td>
    <td>33408</td>
    <td>33408</td>
    <td>211.92</td>
    <td>2085.1</td>
    <td>0.2085</td>
    <td>6.926289</td>
    <td>48.419724</td>
    <td>6° 55' 34.640"</td>
    <td>48° 25' 11.005"</td>
  </tr>
  <tr>
    <td>118</td>
    <td>35168</td>
    <td>35168</td>
    <td>579.26</td>
    <td>18669.7</td>
    <td>1.8670</td>
    <td>6.798052</td>
    <td>48.404261</td>
    <td>6° 47' 52.988"</td>
    <td>48° 24' 15.341"</td>
  </tr>
  <tr>
    <td>119</td>
    <td>35893</td>
    <td>35893</td>
    <td>428.15</td>
    <td>10147.4</td>
    <td>1.0147</td>
    <td>6.797606</td>
    <td>48.398554</td>
    <td>6° 47' 51.382"</td>
    <td>48° 23' 54.795"</td>
  </tr>
  <tr>
    <td>120</td>
    <td>36245</td>
    <td>36245</td>
    <td>224.28</td>
    <td>2418.2</td>
    <td>0.2418</td>
    <td>6.840897</td>
    <td>48.395059</td>
    <td>6° 50' 27.231"</td>
    <td>48° 23' 42.211"</td>
  </tr>
  <tr>
    <td>121</td>
    <td>36312</td>
    <td>36312</td>
    <td>300.00</td>
    <td>5400.0</td>
    <td>0.5400</td>
    <td>6.798424</td>
    <td>48.395126</td>
    <td>6° 47' 54.326"</td>
    <td>48° 23' 42.452"</td>
  </tr>
  <tr>
    <td>122</td>
    <td>36318</td>
    <td>36318</td>
    <td>243.14</td>
    <td>3019.8</td>
    <td>0.3020</td>
    <td>6.839625</td>
    <td>48.394567</td>
    <td>6° 50' 22.651"</td>
    <td>48° 23' 40.440"</td>
  </tr>
  <tr>
    <td>123</td>
    <td>36505</td>
    <td>36505</td>
    <td>212.18</td>
    <td>2112.1</td>
    <td>0.2112</td>
    <td>6.955297</td>
    <td>48.396221</td>
    <td>6° 57' 19.068"</td>
    <td>48° 23' 46.395"</td>
  </tr>
  <tr>
    <td>124</td>
    <td>36598</td>
    <td>36598</td>
    <td>675.47</td>
    <td>23152.2</td>
    <td>2.3152</td>
    <td>6.797437</td>
    <td>48.393078</td>
    <td>6° 47' 50.773"</td>
    <td>48° 23' 35.080"</td>
  </tr>
  <tr>
    <td>125</td>
    <td>36630</td>
    <td>36630</td>
    <td>322.65</td>
    <td>4785.1</td>
    <td>0.4785</td>
    <td>6.836460</td>
    <td>48.393209</td>
    <td>6° 50' 11.255"</td>
    <td>48° 23' 35.554"</td>
  </tr>
  <tr>
    <td>126</td>
    <td>36663</td>
    <td>36663</td>
    <td>300.00</td>
    <td>5400.0</td>
    <td>0.5400</td>
    <td>6.795116</td>
    <td>48.391958</td>
    <td>6° 47' 42.417"</td>
    <td>48° 23' 31.049"</td>
  </tr>
  <tr>
    <td>127</td>
    <td>36893</td>
    <td>36893</td>
    <td>316.54</td>
    <td>5611.9</td>
    <td>0.5612</td>
    <td>6.834163</td>
    <td>48.391695</td>
    <td>6° 50' 2.986"</td>
    <td>48° 23' 30.103"</td>
  </tr>
  <tr>
    <td>128</td>
    <td>37125</td>
    <td>37125</td>
    <td>439.35</td>
    <td>9668.5</td>
    <td>0.9669</td>
    <td>6.831618</td>
    <td>48.390070</td>
    <td>6° 49' 53.826"</td>
    <td>48° 23' 24.251"</td>
  </tr>
  <tr>
    <td>129</td>
    <td>37457</td>
    <td>37457</td>
    <td>120.00</td>
    <td>900.0</td>
    <td>0.0900</td>
    <td>6.807143</td>
    <td>48.385575</td>
    <td>6° 48' 25.716"</td>
    <td>48° 23' 8.070"</td>
  </tr>
  <tr>
    <td>130</td>
    <td>37496</td>
    <td>37496</td>
    <td>120.00</td>
    <td>900.0</td>
    <td>0.0900</td>
    <td>6.806750</td>
    <td>48.385298</td>
    <td>6° 48' 24.300"</td>
    <td>48° 23' 7.072"</td>
  </tr>
  <tr>
    <td>131</td>
    <td>37674</td>
    <td>37674</td>
    <td>733.46</td>
    <td>16331.9</td>
    <td>1.6332</td>
    <td>6.805255</td>
    <td>48.385085</td>
    <td>6° 48' 18.919"</td>
    <td>48° 23' 6.307"</td>
  </tr>
  <tr>
    <td>132</td>
    <td>38162</td>
    <td>38162</td>
    <td>120.00</td>
    <td>900.0</td>
    <td>0.0900</td>
    <td>6.738969</td>
    <td>48.378587</td>
    <td>6° 44' 20.287"</td>
    <td>48° 22' 42.911"</td>
  </tr>
  <tr>
    <td>133</td>
    <td>38311</td>
    <td>38311</td>
    <td>120.00</td>
    <td>900.0</td>
    <td>0.0900</td>
    <td>6.737814</td>
    <td>48.377214</td>
    <td>6° 44' 16.130"</td>
    <td>48° 22' 37.971"</td>
  </tr>
  <tr>
    <td>134</td>
    <td>38421</td>
    <td>38421</td>
    <td>120.00</td>
    <td>900.0</td>
    <td>0.0900</td>
    <td>6.720050</td>
    <td>48.375784</td>
    <td>6° 43' 12.182"</td>
    <td>48° 22' 32.821"</td>
  </tr>
</table>
<table>
  <tr>
    <td>135</td>
    <td>38734</td>
    <td>38734</td>
    <td>120.00</td>
    <td>900.0</td>
    <td>0.0900</td>
    <td>6.738386</td>
    <td>48.373447</td>
    <td>6° 44' 18.190"</td>
    <td>48° 22' 24.408"</td>
  </tr>
  <tr>
    <td>136</td>
    <td>38946</td>
    <td>38946</td>
    <td>169.77</td>
    <td>1229.0</td>
    <td>0.1229</td>
    <td>6.739053</td>
    <td>48.371588</td>
    <td>6° 44' 20.591"</td>
    <td>48° 22' 17.718"</td>
  </tr>
  <tr>
    <td>137</td>
    <td>39216</td>
    <td>39216</td>
    <td>603.36</td>
    <td>13616.3</td>
    <td>1.3616</td>
    <td>6.737833</td>
    <td>48.369714</td>
    <td>6° 44' 16.198"</td>
    <td>48° 22' 10.969"</td>
  </tr>
  <tr>
    <td>138</td>
    <td>39657</td>
    <td>39657</td>
    <td>170.22</td>
    <td>1214.9</td>
    <td>0.1215</td>
    <td>6.734896</td>
    <td>48.365836</td>
    <td>6° 44' 5.626"</td>
    <td>48° 21' 57.009"</td>
  </tr>
  <tr>
    <td>139</td>
    <td>39784</td>
    <td>39784</td>
    <td>410.27</td>
    <td>7473.3</td>
    <td>0.7473</td>
    <td>6.737050</td>
    <td>48.365014</td>
    <td>6° 44' 13.381"</td>
    <td>48° 21' 54.050"</td>
  </tr>
  <tr>
    <td>140</td>
    <td>41319</td>
    <td>41319</td>
    <td>120.00</td>
    <td>900.0</td>
    <td>0.0900</td>
    <td>7.199109</td>
    <td>48.359451</td>
    <td>7° 11' 56.791"</td>
    <td>48° 21' 34.025"</td>
  </tr>
  <tr>
    <td>141</td>
    <td>41613</td>
    <td>41613</td>
    <td>168.12</td>
    <td>1349.7</td>
    <td>0.1350</td>
    <td>7.146976</td>
    <td>48.356314</td>
    <td>7° 8' 49.112"</td>
    <td>48° 21' 22.732"</td>
  </tr>
  <tr>
    <td>142</td>
    <td>43190</td>
    <td>43190</td>
    <td>120.00</td>
    <td>900.0</td>
    <td>0.0900</td>
    <td>7.317063</td>
    <td>48.344233</td>
    <td>7° 19' 1.429"</td>
    <td>48° 20' 39.240"</td>
  </tr>
  <tr>
    <td>143</td>
    <td>46284</td>
    <td>46284</td>
    <td>212.18</td>
    <td>2112.1</td>
    <td>0.2112</td>
    <td>7.033681</td>
    <td>48.311023</td>
    <td>7° 2' 1.250"</td>
    <td>48° 18' 39.684"</td>
  </tr>
  <tr>
    <td>144</td>
    <td>46527</td>
    <td>46527</td>
    <td>598.96</td>
    <td>14064.8</td>
    <td>1.4065</td>
    <td>7.032518</td>
    <td>48.309014</td>
    <td>7° 1' 57.066"</td>
    <td>48° 18' 32.450"</td>
  </tr>
  <tr>
    <td>145</td>
    <td>46762</td>
    <td>46762</td>
    <td>1328.11</td>
    <td>54443.6</td>
    <td>3.0803</td>
    <td>7.029074</td>
    <td>48.306753</td>
    <td>7° 1' 44.667"</td>
    <td>48° 18' 24.310"</td>
  </tr>
  <tr>
    <td>146</td>
    <td>46825</td>
    <td>46825</td>
    <td>907.12</td>
    <td>24968.1</td>
    <td>0.9474</td>
    <td>7.026662</td>
    <td>48.305276</td>
    <td>7° 1' 35.982"</td>
    <td>48° 18' 18.994"</td>
  </tr>
  <tr>
    <td>147</td>
    <td>46919</td>
    <td>46919</td>
    <td>168.01</td>
    <td>1393.6</td>
    <td>0.1394</td>
    <td>7.173274</td>
    <td>48.305336</td>
    <td>7° 10' 23.786"</td>
    <td>48° 18' 19.211"</td>
  </tr>
  <tr>
    <td>148</td>
    <td>47124</td>
    <td>47124</td>
    <td>528.63</td>
    <td>15167.4</td>
    <td>1.5167</td>
    <td>7.166364</td>
    <td>48.303774</td>
    <td>7° 9' 58.910"</td>
    <td>48° 18' 13.588"</td>
  </tr>
  <tr>
    <td>149</td>
    <td>47204</td>
    <td>47204</td>
    <td>360.00</td>
    <td>7200.0</td>
    <td>0.7200</td>
    <td>6.700044</td>
    <td>48.293733</td>
    <td>6° 42' 0.158"</td>
    <td>48° 17' 37.439"</td>
  </tr>
  <tr>
    <td>150</td>
    <td>48931</td>
    <td>48931</td>
    <td>990.36</td>
    <td>27239.2</td>
    <td>1.7510</td>
    <td>6.758676</td>
    <td>48.276533</td>
    <td>6° 45' 31.232"</td>
    <td>48° 16' 35.518"</td>
  </tr>
  <tr>
    <td>151</td>
    <td>49079</td>
    <td>49079</td>
    <td>662.65</td>
    <td>16955.6</td>
    <td>1.6956</td>
    <td>6.786183</td>
    <td>48.275179</td>
    <td>6° 47' 10.259"</td>
    <td>48° 16' 30.643"</td>
  </tr>
  <tr>
    <td>152</td>
    <td>49080</td>
    <td>49080</td>
    <td>265.37</td>
    <td>3187.2</td>
    <td>0.3187</td>
    <td>6.789455</td>
    <td>48.274633</td>
    <td>6° 47' 22.038"</td>
    <td>48° 16' 28.680"</td>
  </tr>
  <tr>
    <td>153</td>
    <td>50527</td>
    <td>50527</td>
    <td>120.00</td>
    <td>900.0</td>
    <td>0.0900</td>
    <td>7.190664</td>
    <td>48.266198</td>
    <td>7° 11' 26.391"</td>
    <td>48° 15' 58.311"</td>
  </tr>
  <tr>
    <td>154</td>
    <td>50528</td>
    <td>50528</td>
    <td>233.55</td>
    <td>3314.8</td>
    <td>0.3315</td>
    <td>7.192318</td>
    <td>48.266298</td>
    <td>7° 11' 32.343"</td>
    <td>48° 15' 58.673"</td>
  </tr>
  <tr>
    <td>155</td>
    <td>50560</td>
    <td>50560</td>
    <td>505.34</td>
    <td>12019.4</td>
    <td>1.2019</td>
    <td>7.194604</td>
    <td>48.266396</td>
    <td>7° 11' 40.575"</td>
    <td>48° 15' 59.027"</td>
  </tr>
  <tr>
    <td>156</td>
    <td>50730</td>
    <td>50730</td>
    <td>265.37</td>
    <td>3187.2</td>
    <td>0.3187</td>
    <td>7.189708</td>
    <td>48.264628</td>
    <td>7° 11' 22.950"</td>
    <td>48° 15' 52.662"</td>
  </tr>
  <tr>
    <td>157</td>
    <td>50871</td>
    <td>50871</td>
    <td>348.02</td>
    <td>6707.5</td>
    <td>0.6708</td>
    <td>7.061401</td>
    <td>48.261315</td>
    <td>7° 3' 41.043"</td>
    <td>48° 15' 40.735"</td>
  </tr>
  <tr>
    <td>158</td>
    <td>50873</td>
    <td>50873</td>
    <td>355.93</td>
    <td>6448.4</td>
    <td>0.6448</td>
    <td>7.188392</td>
    <td>48.263273</td>
    <td>7° 11' 18.211"</td>
    <td>48° 15' 47.784"</td>
  </tr>
  <tr>
    <td>159</td>
    <td>53148</td>
    <td>53148</td>
    <td>227.83</td>
    <td>1949.7</td>
    <td>0.1950</td>
    <td>7.033794</td>
    <td>48.241288</td>
    <td>7° 2' 1.657"</td>
    <td>48° 14' 28.638"</td>
  </tr>
  <tr>
    <td>160</td>
    <td>53287</td>
    <td>53287</td>
    <td>698.67</td>
    <td>14754.4</td>
    <td>1.4754</td>
    <td>7.034697</td>
    <td>48.239654</td>
    <td>7° 2' 4.909"</td>
    <td>48° 14' 22.753"</td>
  </tr>
  <tr>
    <td>161</td>
    <td>53390</td>
    <td>53390</td>
    <td>251.91</td>
    <td>3318.7</td>
    <td>0.3319</td>
    <td>7.036319</td>
    <td>48.238000</td>
    <td>7° 2' 10.750"</td>
    <td>48° 14' 16.799"</td>
  </tr>
</table>
<table>
  <tr>
    <td>162</td>
    <td>53446</td>
    <td>53446</td>
    <td>168.13</td>
    <td>1349.9</td>
    <td>0.1350</td>
    <td>7.039235</td>
    <td>48.235648</td>
    <td>7° 2' 21.244"</td>
    <td>48° 14' 8.334"</td>
  </tr>
  <tr>
    <td>163</td>
    <td>53472</td>
    <td>53472</td>
    <td>371.44</td>
    <td>5772.5</td>
    <td>0.5772</td>
    <td>7.037205</td>
    <td>48.237221</td>
    <td>7° 2' 13.939"</td>
    <td>48° 14' 13.997"</td>
  </tr>
  <tr>
    <td>164</td>
    <td>53591</td>
    <td>53591</td>
    <td>168.02</td>
    <td>1393.7</td>
    <td>0.1394</td>
    <td>7.038172</td>
    <td>48.236322</td>
    <td>7° 2' 17.419"</td>
    <td>48° 14' 10.758"</td>
  </tr>
  <tr>
    <td>165</td>
    <td>56532</td>
    <td>56532</td>
    <td>1858.75</td>
    <td>73516.3</td>
    <td>0.5567</td>
    <td>6.669043</td>
    <td>48.202887</td>
    <td>6° 40' 8.556"</td>
    <td>48° 12' 10.392"</td>
  </tr>
  <tr>
    <td>166</td>
    <td>59374</td>
    <td>59374</td>
    <td>170.22</td>
    <td>1215.0</td>
    <td>0.1215</td>
    <td>7.239241</td>
    <td>48.188421</td>
    <td>7° 14' 21.269"</td>
    <td>48° 11' 18.316"</td>
  </tr>
  <tr>
    <td>167</td>
    <td>59661</td>
    <td>59661</td>
    <td>637.30</td>
    <td>24506.0</td>
    <td>2.4506</td>
    <td>7.236356</td>
    <td>48.186760</td>
    <td>7° 14' 10.881"</td>
    <td>48° 11' 12.334"</td>
  </tr>
  <tr>
    <td>168</td>
    <td>59763</td>
    <td>59763</td>
    <td>352.86</td>
    <td>5938.6</td>
    <td>0.5939</td>
    <td>7.233649</td>
    <td>48.185838</td>
    <td>7° 14' 1.138"</td>
    <td>48° 11' 9.017"</td>
  </tr>
  <tr>
    <td>169</td>
    <td>59832</td>
    <td>59832</td>
    <td>224.27</td>
    <td>2418.1</td>
    <td>0.2418</td>
    <td>7.238307</td>
    <td>48.187866</td>
    <td>7° 14' 17.904"</td>
    <td>48° 11' 16.319"</td>
  </tr>
  <tr>
    <td>170</td>
    <td>59855</td>
    <td>59855</td>
    <td>212.17</td>
    <td>2112.0</td>
    <td>0.2112</td>
    <td>7.232139</td>
    <td>48.184876</td>
    <td>7° 13' 55.701"</td>
    <td>48° 11' 5.554"</td>
  </tr>
  <tr>
    <td>171</td>
    <td>59999</td>
    <td>59999</td>
    <td>212.17</td>
    <td>2112.0</td>
    <td>0.2112</td>
    <td>7.232580</td>
    <td>48.183803</td>
    <td>7° 13' 57.286"</td>
    <td>48° 11' 1.691"</td>
  </tr>
  <tr>
    <td>172</td>
    <td>65019</td>
    <td>65019</td>
    <td>1299.96</td>
    <td>60812.4</td>
    <td>6.0812</td>
    <td>6.892683</td>
    <td>48.126341</td>
    <td>6° 53' 33.660"</td>
    <td>48° 7' 34.827"</td>
  </tr>
  <tr>
    <td>173</td>
    <td>66819</td>
    <td>66819</td>
    <td>120.00</td>
    <td>900.0</td>
    <td>0.0900</td>
    <td>7.046643</td>
    <td>48.130768</td>
    <td>7° 2' 47.914"</td>
    <td>48° 7' 50.766"</td>
  </tr>
  <tr>
    <td>174</td>
    <td>66870</td>
    <td>66870</td>
    <td>180.00</td>
    <td>1800.0</td>
    <td>0.1800</td>
    <td>7.046255</td>
    <td>48.130357</td>
    <td>7° 2' 46.519"</td>
    <td>48° 7' 49.284"</td>
  </tr>
  <tr>
    <td>175</td>
    <td>66952</td>
    <td>66952</td>
    <td>120.00</td>
    <td>900.0</td>
    <td>0.0900</td>
    <td>7.045475</td>
    <td>48.129669</td>
    <td>7° 2' 43.710"</td>
    <td>48° 7' 46.807"</td>
  </tr>
  <tr>
    <td>176</td>
    <td>67152</td>
    <td>67152</td>
    <td>240.00</td>
    <td>3600.0</td>
    <td>0.3600</td>
    <td>7.044302</td>
    <td>48.127831</td>
    <td>7° 2' 39.488"</td>
    <td>48° 7' 40.191"</td>
  </tr>
  <tr>
    <td>177</td>
    <td>67234</td>
    <td>67234</td>
    <td>390.24</td>
    <td>5818.0</td>
    <td>0.5818</td>
    <td>7.044891</td>
    <td>48.128745</td>
    <td>7° 2' 41.609"</td>
    <td>48° 7' 43.481"</td>
  </tr>
  <tr>
    <td>178</td>
    <td>67430</td>
    <td>67430</td>
    <td>335.42</td>
    <td>4966.8</td>
    <td>0.4967</td>
    <td>7.043572</td>
    <td>48.126814</td>
    <td>7° 2' 36.861"</td>
    <td>48° 7' 36.529"</td>
  </tr>
  <tr>
    <td>179</td>
    <td>67497</td>
    <td>67497</td>
    <td>240.00</td>
    <td>2700.0</td>
    <td>0.2700</td>
    <td>7.043221</td>
    <td>48.125311</td>
    <td>7° 2' 35.597"</td>
    <td>48° 7' 31.121"</td>
  </tr>
  <tr>
    <td>180</td>
    <td>67551</td>
    <td>67551</td>
    <td>180.00</td>
    <td>1800.0</td>
    <td>0.1800</td>
    <td>7.042431</td>
    <td>48.124893</td>
    <td>7° 2' 32.751"</td>
    <td>48° 7' 29.615"</td>
  </tr>
  <tr>
    <td>181</td>
    <td>67582</td>
    <td>67582</td>
    <td>224.28</td>
    <td>2418.2</td>
    <td>0.2418</td>
    <td>7.047014</td>
    <td>48.131375</td>
    <td>7° 2' 49.251"</td>
    <td>48° 7' 52.951"</td>
  </tr>
  <tr>
    <td>182</td>
    <td>67585</td>
    <td>67585</td>
    <td>224.28</td>
    <td>2418.2</td>
    <td>0.2418</td>
    <td>7.043302</td>
    <td>48.125599</td>
    <td>7° 2' 35.886"</td>
    <td>48° 7' 32.155"</td>
  </tr>
  <tr>
    <td>183</td>
    <td>67652</td>
    <td>67652</td>
    <td>267.93</td>
    <td>3090.3</td>
    <td>0.3090</td>
    <td>7.042302</td>
    <td>48.123969</td>
    <td>7° 2' 32.288"</td>
    <td>48° 7' 26.287"</td>
  </tr>
  <tr>
    <td>184</td>
    <td>67685</td>
    <td>67685</td>
    <td>163.00</td>
    <td>1255.7</td>
    <td>0.1256</td>
    <td>6.895007</td>
    <td>48.121207</td>
    <td>6° 53' 42.026"</td>
    <td>48° 7' 16.344"</td>
  </tr>
  <tr>
    <td>185</td>
    <td>67772</td>
    <td>67772</td>
    <td>120.00</td>
    <td>900.0</td>
    <td>0.0900</td>
    <td>6.893437</td>
    <td>48.120503</td>
    <td>6° 53' 36.375"</td>
    <td>48° 7' 13.811"</td>
  </tr>
  <tr>
    <td>186</td>
    <td>67898</td>
    <td>67898</td>
    <td>219.57</td>
    <td>2034.2</td>
    <td>0.2034</td>
    <td>6.887958</td>
    <td>48.119675</td>
    <td>6° 53' 16.650"</td>
    <td>48° 7' 10.830"</td>
  </tr>
  <tr>
    <td>187</td>
    <td>67931</td>
    <td>67931</td>
    <td>169.76</td>
    <td>1228.8</td>
    <td>0.1229</td>
    <td>6.884588</td>
    <td>48.119380</td>
    <td>6° 53' 4.515"</td>
    <td>48° 7' 9.769"</td>
  </tr>
  <tr>
    <td>188</td>
    <td>70038</td>
    <td>70038</td>
    <td>1846.30</td>
    <td>84529.5</td>
    <td>2.1201</td>
    <td>6.646952</td>
    <td>48.099448</td>
    <td>6° 38' 49.029"</td>
    <td>48° 5' 58.011"</td>
  </tr>
</table>
<table>
  <tr>
    <td>189</td>
    <td>70760</td>
    <td>70760</td>
    <td>644.05</td>
    <td>14818.3</td>
    <td>1.4818</td>
    <td>6.821553</td>
    <td>48.094559</td>
    <td>6° 49' 17.591"</td>
    <td>48° 5' 40.413"</td>
  </tr>
  <tr>
    <td>190</td>
    <td>70899</td>
    <td>70899</td>
    <td>347.76</td>
    <td>5710.5</td>
    <td>0.5710</td>
    <td>6.828641</td>
    <td>48.093164</td>
    <td>6° 49' 43.106"</td>
    <td>48° 5' 35.390"</td>
  </tr>
  <tr>
    <td>191</td>
    <td>70981</td>
    <td>70981</td>
    <td>217.54</td>
    <td>2316.4</td>
    <td>0.2316</td>
    <td>6.832718</td>
    <td>48.092224</td>
    <td>6° 49' 57.786"</td>
    <td>48° 5' 32.005"</td>
  </tr>
  <tr>
    <td>192</td>
    <td>71079</td>
    <td>71079</td>
    <td>312.89</td>
    <td>3487.2</td>
    <td>0.3487</td>
    <td>6.701415</td>
    <td>48.087717</td>
    <td>6° 42' 5.093"</td>
    <td>48° 5' 15.780"</td>
  </tr>
  <tr>
    <td>193</td>
    <td>71084</td>
    <td>71084</td>
    <td>168.13</td>
    <td>1349.9</td>
    <td>0.1350</td>
    <td>6.846103</td>
    <td>48.091283</td>
    <td>6° 50' 45.969"</td>
    <td>48° 5' 28.617"</td>
  </tr>
  <tr>
    <td>194</td>
    <td>71133</td>
    <td>71133</td>
    <td>318.43</td>
    <td>5590.2</td>
    <td>0.5590</td>
    <td>6.851443</td>
    <td>48.091051</td>
    <td>6° 51' 5.195"</td>
    <td>48° 5' 27.782"</td>
  </tr>
  <tr>
    <td>195</td>
    <td>71558</td>
    <td>71558</td>
    <td>162.99</td>
    <td>1255.7</td>
    <td>0.1256</td>
    <td>6.694672</td>
    <td>48.087413</td>
    <td>6° 41' 40.819"</td>
    <td>48° 5' 14.687"</td>
  </tr>
  <tr>
    <td>196</td>
    <td>71593</td>
    <td>71593</td>
    <td>213.01</td>
    <td>2132.1</td>
    <td>0.2132</td>
    <td>6.695531</td>
    <td>48.086563</td>
    <td>6° 41' 43.911"</td>
    <td>48° 5' 11.626"</td>
  </tr>
  <tr>
    <td>197</td>
    <td>72632</td>
    <td>72632</td>
    <td>375.80</td>
    <td>7542.6</td>
    <td>0.7543</td>
    <td>7.253864</td>
    <td>48.083853</td>
    <td>7° 15' 13.910"</td>
    <td>48° 5' 1.872"</td>
  </tr>
  <tr>
    <td>198</td>
    <td>72954</td>
    <td>72954</td>
    <td>120.00</td>
    <td>900.0</td>
    <td>0.0900</td>
    <td>7.254784</td>
    <td>48.080673</td>
    <td>7° 15' 17.224"</td>
    <td>48° 4' 50.424"</td>
  </tr>
  <tr>
    <td>199</td>
    <td>74046</td>
    <td>74046</td>
    <td>560.86</td>
    <td>15661.7</td>
    <td>1.5662</td>
    <td>6.480911</td>
    <td>48.058205</td>
    <td>6° 28' 51.281"</td>
    <td>48° 3' 29.539"</td>
  </tr>
  <tr>
    <td>200</td>
    <td>74214</td>
    <td>74214</td>
    <td>336.55</td>
    <td>5010.0</td>
    <td>0.5010</td>
    <td>6.474940</td>
    <td>48.055932</td>
    <td>6° 28' 29.783"</td>
    <td>48° 3' 21.357"</td>
  </tr>
  <tr>
    <td>201</td>
    <td>76429</td>
    <td>76429</td>
    <td>314.99</td>
    <td>4549.8</td>
    <td>0.4550</td>
    <td>7.092264</td>
    <td>48.050665</td>
    <td>7° 5' 32.152"</td>
    <td>48° 3' 2.395"</td>
  </tr>
  <tr>
    <td>202</td>
    <td>76445</td>
    <td>76445</td>
    <td>359.81</td>
    <td>4711.8</td>
    <td>0.4712</td>
    <td>6.724043</td>
    <td>48.043516</td>
    <td>6° 43' 26.556"</td>
    <td>48° 2' 36.658"</td>
  </tr>
  <tr>
    <td>203</td>
    <td>76839</td>
    <td>76839</td>
    <td>162.99</td>
    <td>1255.7</td>
    <td>0.1256</td>
    <td>6.626447</td>
    <td>48.037087</td>
    <td>6° 37' 35.209"</td>
    <td>48° 2' 13.513"</td>
  </tr>
  <tr>
    <td>204</td>
    <td>77043</td>
    <td>77043</td>
    <td>366.82</td>
    <td>5796.1</td>
    <td>0.5796</td>
    <td>6.713402</td>
    <td>48.037348</td>
    <td>6° 42' 48.249"</td>
    <td>48° 2' 14.453"</td>
  </tr>
  <tr>
    <td>205</td>
    <td>77091</td>
    <td>77091</td>
    <td>168.14</td>
    <td>1350.0</td>
    <td>0.1350</td>
    <td>6.715717</td>
    <td>48.037613</td>
    <td>6° 42' 56.582"</td>
    <td>48° 2' 15.406"</td>
  </tr>
  <tr>
    <td>206</td>
    <td>77113</td>
    <td>77113</td>
    <td>1375.30</td>
    <td>39462.9</td>
    <td>3.9463</td>
    <td>7.090971</td>
    <td>48.046868</td>
    <td>7° 5' 27.494"</td>
    <td>48° 2' 48.723"</td>
  </tr>
  <tr>
    <td>207</td>
    <td>77429</td>
    <td>77429</td>
    <td>224.29</td>
    <td>2418.5</td>
    <td>0.2419</td>
    <td>6.712583</td>
    <td>48.036340</td>
    <td>6° 42' 45.300"</td>
    <td>48° 2' 10.823"</td>
  </tr>
  <tr>
    <td>208</td>
    <td>77505</td>
    <td>77505</td>
    <td>211.92</td>
    <td>2085.1</td>
    <td>0.2085</td>
    <td>6.629977</td>
    <td>48.033491</td>
    <td>6° 37' 47.917"</td>
    <td>48° 2' 0.566"</td>
  </tr>
  <tr>
    <td>209</td>
    <td>77571</td>
    <td>77571</td>
    <td>632.38</td>
    <td>13184.9</td>
    <td>1.3185</td>
    <td>7.088228</td>
    <td>48.042689</td>
    <td>7° 5' 17.621"</td>
    <td>48° 2' 33.679"</td>
  </tr>
  <tr>
    <td>210</td>
    <td>77587</td>
    <td>77587</td>
    <td>172.80</td>
    <td>1295.6</td>
    <td>0.1296</td>
    <td>6.716405</td>
    <td>48.037092</td>
    <td>6° 42' 59.056"</td>
    <td>48° 2' 13.532"</td>
  </tr>
  <tr>
    <td>211</td>
    <td>77836</td>
    <td>77836</td>
    <td>360.00</td>
    <td>7200.0</td>
    <td>0.7200</td>
    <td>6.703131</td>
    <td>48.030983</td>
    <td>6° 42' 11.271"</td>
    <td>48° 1' 51.539"</td>
  </tr>
  <tr>
    <td>212</td>
    <td>77911</td>
    <td>77911</td>
    <td>311.90</td>
    <td>4424.3</td>
    <td>0.4424</td>
    <td>6.705636</td>
    <td>48.030626</td>
    <td>6° 42' 20.290"</td>
    <td>48° 1' 50.255"</td>
  </tr>
  <tr>
    <td>213</td>
    <td>77965</td>
    <td>77965</td>
    <td>858.34</td>
    <td>19206.7</td>
    <td>1.9207</td>
    <td>7.085436</td>
    <td>48.040173</td>
    <td>7° 5' 7.569"</td>
    <td>48° 2' 24.624"</td>
  </tr>
  <tr>
    <td>214</td>
    <td>77974</td>
    <td>77974</td>
    <td>224.29</td>
    <td>2418.5</td>
    <td>0.2419</td>
    <td>6.708790</td>
    <td>48.030025</td>
    <td>6° 42' 31.645"</td>
    <td>48° 1' 48.089"</td>
  </tr>
  <tr>
    <td>215</td>
    <td>78022</td>
    <td>78022</td>
    <td>180.00</td>
    <td>1800.0</td>
    <td>0.1800</td>
    <td>7.085566</td>
    <td>48.038163</td>
    <td>7° 5' 8.039"</td>
    <td>48° 2' 17.388"</td>
  </tr>
</table>
<table>
  <tr>
    <td>216</td>
    <td>78098</td>
    <td>78098</td>
    <td>219.55</td>
    <td>2033.8</td>
    <td>0.2034</td>
    <td>7.084501</td>
    <td>48.037653</td>
    <td>7° 5' 4.202"</td>
    <td>48° 2' 15.551"</td>
  </tr>
  <tr>
    <td>217</td>
    <td>78195</td>
    <td>78195</td>
    <td>376.47</td>
    <td>6022.5</td>
    <td>0.6023</td>
    <td>6.874183</td>
    <td>48.028792</td>
    <td>6° 52' 27.061"</td>
    <td>48° 1' 43.652"</td>
  </tr>
  <tr>
    <td>218</td>
    <td>78471</td>
    <td>78471</td>
    <td>361.30</td>
    <td>4965.9</td>
    <td>0.4966</td>
    <td>6.871796</td>
    <td>48.029329</td>
    <td>6° 52' 18.466"</td>
    <td>48° 1' 45.585"</td>
  </tr>
  <tr>
    <td>219</td>
    <td>78553</td>
    <td>78553</td>
    <td>1145.74</td>
    <td>48137.7</td>
    <td>4.8138</td>
    <td>6.867195</td>
    <td>48.030042</td>
    <td>6° 52' 1.902"</td>
    <td>48° 1' 48.152"</td>
  </tr>
  <tr>
    <td>220</td>
    <td>78944</td>
    <td>78944</td>
    <td>335.31</td>
    <td>5994.1</td>
    <td>0.5994</td>
    <td>6.876629</td>
    <td>48.027674</td>
    <td>6° 52' 35.865"</td>
    <td>48° 1' 39.627"</td>
  </tr>
  <tr>
    <td>221</td>
    <td>79182</td>
    <td>79182</td>
    <td>1125.30</td>
    <td>59875.2</td>
    <td>3.3453</td>
    <td>6.409031</td>
    <td>48.017559</td>
    <td>6° 24' 32.511"</td>
    <td>48° 1' 3.212"</td>
  </tr>
  <tr>
    <td>222</td>
    <td>79200</td>
    <td>79200</td>
    <td>538.76</td>
    <td>9668.9</td>
    <td>0.9669</td>
    <td>6.879022</td>
    <td>48.026555</td>
    <td>6° 52' 44.478"</td>
    <td>48° 1' 35.599"</td>
  </tr>
  <tr>
    <td>223</td>
    <td>79729</td>
    <td>79729</td>
    <td>218.93</td>
    <td>2467.5</td>
    <td>0.2468</td>
    <td>6.830954</td>
    <td>48.020969</td>
    <td>6° 49' 51.434"</td>
    <td>48° 1' 15.487"</td>
  </tr>
  <tr>
    <td>224</td>
    <td>80060</td>
    <td>80060</td>
    <td>232.80</td>
    <td>2195.7</td>
    <td>0.2196</td>
    <td>6.832337</td>
    <td>48.018405</td>
    <td>6° 49' 56.413"</td>
    <td>48° 1' 6.259"</td>
  </tr>
  <tr>
    <td>225</td>
    <td>80135</td>
    <td>80135</td>
    <td>300.00</td>
    <td>4500.0</td>
    <td>0.4500</td>
    <td>6.831154</td>
    <td>48.017957</td>
    <td>6° 49' 52.154"</td>
    <td>48° 1' 4.646"</td>
  </tr>
  <tr>
    <td>226</td>
    <td>80269</td>
    <td>80269</td>
    <td>367.54</td>
    <td>5611.9</td>
    <td>0.5612</td>
    <td>6.832784</td>
    <td>48.017249</td>
    <td>6° 49' 58.024"</td>
    <td>48° 1' 2.098"</td>
  </tr>
  <tr>
    <td>227</td>
    <td>80462</td>
    <td>80462</td>
    <td>493.93</td>
    <td>10697.1</td>
    <td>1.0697</td>
    <td>6.832227</td>
    <td>48.016283</td>
    <td>6° 49' 56.017"</td>
    <td>48° 0' 58.620"</td>
  </tr>
  <tr>
    <td>228</td>
    <td>80504</td>
    <td>80504</td>
    <td>170.21</td>
    <td>1214.7</td>
    <td>0.1215</td>
    <td>6.831521</td>
    <td>48.015411</td>
    <td>6° 49' 53.476"</td>
    <td>48° 0' 55.480"</td>
  </tr>
  <tr>
    <td>229</td>
    <td>81342</td>
    <td>81342</td>
    <td>466.54</td>
    <td>13245.5</td>
    <td>1.3245</td>
    <td>6.598190</td>
    <td>48.005240</td>
    <td>6° 35' 53.485"</td>
    <td>48° 0' 18.865"</td>
  </tr>
  <tr>
    <td>230</td>
    <td>81629</td>
    <td>81629</td>
    <td>384.24</td>
    <td>7545.9</td>
    <td>0.7546</td>
    <td>6.513222</td>
    <td>48.001015</td>
    <td>6° 30' 47.599"</td>
    <td>48° 0' 3.655"</td>
  </tr>
  <tr>
    <td>231</td>
    <td>81715</td>
    <td>81715</td>
    <td>2942.12</td>
    <td>172326.8</td>
    <td>3.6673</td>
    <td>6.679848</td>
    <td>48.004844</td>
    <td>6° 40' 47.452"</td>
    <td>48° 0' 17.439"</td>
  </tr>
  <tr>
    <td>232</td>
    <td>81762</td>
    <td>81762</td>
    <td>245.30</td>
    <td>3308.5</td>
    <td>0.3308</td>
    <td>6.601119</td>
    <td>48.001926</td>
    <td>6° 36' 4.029"</td>
    <td>48° 0' 6.932"</td>
  </tr>
  <tr>
    <td>233</td>
    <td>81871</td>
    <td>81871</td>
    <td>169.76</td>
    <td>1228.8</td>
    <td>0.1229</td>
    <td>6.456407</td>
    <td>47.997861</td>
    <td>6° 27' 23.065"</td>
    <td>47° 59' 52.300"</td>
  </tr>
  <tr>
    <td>234</td>
    <td>82586</td>
    <td>82586</td>
    <td>405.02</td>
    <td>6833.5</td>
    <td>0.6833</td>
    <td>6.596993</td>
    <td>47.996275</td>
    <td>6° 35' 49.176"</td>
    <td>47° 59' 46.590"</td>
  </tr>
  <tr>
    <td>235</td>
    <td>82935</td>
    <td>82935</td>
    <td>483.71</td>
    <td>10542.9</td>
    <td>1.0543</td>
    <td>6.597978</td>
    <td>47.994046</td>
    <td>6° 35' 52.721"</td>
    <td>47° 59' 38.566"</td>
  </tr>
  <tr>
    <td>236</td>
    <td>83349</td>
    <td>83349</td>
    <td>211.92</td>
    <td>2085.1</td>
    <td>0.2085</td>
    <td>6.685108</td>
    <td>47.991903</td>
    <td>6° 41' 6.388"</td>
    <td>47° 59' 30.850"</td>
  </tr>
  <tr>
    <td>237</td>
    <td>83422</td>
    <td>83422</td>
    <td>608.80</td>
    <td>17008.7</td>
    <td>1.7009</td>
    <td>6.657683</td>
    <td>47.991357</td>
    <td>6° 39' 27.658"</td>
    <td>47° 59' 28.885"</td>
  </tr>
  <tr>
    <td>238</td>
    <td>83523</td>
    <td>83523</td>
    <td>120.00</td>
    <td>900.0</td>
    <td>0.0900</td>
    <td>6.380711</td>
    <td>47.984162</td>
    <td>6° 22' 50.558"</td>
    <td>47° 59' 2.983"</td>
  </tr>
  <tr>
    <td>239</td>
    <td>83692</td>
    <td>83692</td>
    <td>472.93</td>
    <td>10273.4</td>
    <td>1.0273</td>
    <td>6.677225</td>
    <td>47.989840</td>
    <td>6° 40' 38.010"</td>
    <td>47° 59' 23.424"</td>
  </tr>
  <tr>
    <td>240</td>
    <td>85133</td>
    <td>85133</td>
    <td>931.19</td>
    <td>22064.8</td>
    <td>1.5282</td>
    <td>6.712040</td>
    <td>47.932766</td>
    <td>6° 42' 43.346"</td>
    <td>47° 55' 57.959"</td>
  </tr>
  <tr>
    <td>241</td>
    <td>86054</td>
    <td>86054</td>
    <td>120.00</td>
    <td>900.0</td>
    <td>0.0900</td>
    <td>6.391648</td>
    <td>47.966594</td>
    <td>6° 23' 29.934"</td>
    <td>47° 57' 59.740"</td>
  </tr>
  <tr>
    <td>242</td>
    <td>86056</td>
    <td>86056</td>
    <td>321.00</td>
    <td>4807.4</td>
    <td>0.4807</td>
    <td>6.398728</td>
    <td>47.966991</td>
    <td>6° 23' 55.420"</td>
    <td>47° 58' 1.168"</td>
  </tr>
</table>
<table><tr><td>243</td><td>86119</td><td>86119</td><td>263.03</td><td>3038.1</td><td>0.3038</td><td>6.403234</td><td>47.966478</td><td>6° 24' 11.644"</td><td>47° 57' 59.319"</td></tr><tr><td>244</td><td>86199</td><td>86199</td><td>209.30</td><td>2390.6</td><td>0.2391</td><td>6.392786</td><td>47.965736</td><td>6° 23' 34.031"</td><td>47° 57' 56.650"</td></tr><tr><td>245</td><td>86233</td><td>86233</td><td>120.00</td><td>900.0</td><td>0.0900</td><td>6.393322</td><td>47.965283</td><td>6° 23' 35.960"</td><td>47° 57' 55.018"</td></tr><tr><td>246</td><td>92585</td><td>92585</td><td>616.48</td><td>17300.4</td><td>1.7300</td><td>7.093626</td><td>47.934453</td><td>7° 5' 37.055"</td><td>47° 56' 4.030"</td></tr><tr><td>247</td><td>92626</td><td>92626</td><td>372.61</td><td>7408.7</td><td>0.7409</td><td>6.639460</td><td>47.928902</td><td>6° 38' 22.058"</td><td>47° 55' 44.047"</td></tr><tr><td>248</td><td>92813</td><td>92813</td><td>275.19</td><td>3381.6</td><td>0.3382</td><td>7.091395</td><td>47.931651</td><td>7° 5' 29.022"</td><td>47° 55' 53.944"</td></tr><tr><td>249</td><td>92847</td><td>92847</td><td>120.00</td><td>900.0</td><td>0.0900</td><td>7.090750</td><td>47.931127</td><td>7° 5' 26.699"</td><td>47° 55' 52.056"</td></tr><tr><td>250</td><td>93012</td><td>93012</td><td>300.00</td><td>4500.0</td><td>0.4500</td><td>7.086603</td><td>47.930302</td><td>7° 5' 11.770"</td><td>47° 55' 49.086"</td></tr><tr><td>251</td><td>93088</td><td>93088</td><td>221.91</td><td>2712.9</td><td>0.2713</td><td>7.089391</td><td>47.929626</td><td>7° 5' 21.806"</td><td>47° 55' 46.653"</td></tr><tr><td>252</td><td>93223</td><td>93223</td><td>1311.92</td><td>35292.2</td><td>3.5292</td><td>7.088526</td><td>47.929953</td><td>7° 5' 18.692"</td><td>47° 55' 47.830"</td></tr><tr><td>253</td><td>93291</td><td>93291</td><td>337.10</td><td>5167.5</td><td>0.5168</td><td>7.087112</td><td>47.928197</td><td>7° 5' 13.603"</td><td>47° 55' 41.511"</td></tr><tr><td>254</td><td>93874</td><td>93874</td><td>800.60</td><td>19383.3</td><td>1.9383</td><td>7.086103</td><td>47.925085</td><td>7° 5' 9.970"</td><td>47° 55' 30.304"</td></tr><tr><td>255</td><td>95141</td><td>95141</td><td>526.29</td><td>11000.6</td><td>7.5014</td><td>7.088169</td><td>47.918369</td><td>7° 5' 17.408"</td><td>47° 55' 6.127"</td></tr><tr><td>256</td><td>100805</td><td>100805</td><td>651.41</td><td>18439.0</td><td>1.8439</td><td>6.579443</td><td>47.864091</td><td>6° 34' 45.994"</td><td>47° 51' 50.728"</td></tr><tr><td>257</td><td>103531</td><td>103531</td><td>717.76</td><td>22496.0</td><td>2.2496</td><td>6.582815</td><td>47.870216</td><td>6° 34' 58.134"</td><td>47° 52' 12.777"</td></tr></table>