5 Spatial modeling

Let’s start to explore the data with the tools we already know. In this case, we will use the variable clay as example.

soil_data <- read.table("./Data/soil_data.txt", header = TRUE)
plot(soil_data$x, 
     soil_data$y, 
     cex=soil_data$clay*0.05)

We can produce a similar result with more specialized packages as sp

# install.packages("sp")
library(sp)

# lets make a copy to use the raw data later
soil_data_sp <- soil_data

# We can transform the data as spatial 
coordinates(soil_data_sp) = ~x + y
class(soil_data_sp)

## [1] "SpatialPointsDataFrame"
## attr(,"package")
## [1] "sp"

# Plot the spatial data frame
bubble(soil_data_sp,"clay")

Model: Limited representation of reality

Spatial modeling: Image of states or processes of locatable objects and their characteristic values

Data for spatial modeling:

• Position in space (e.g. xy-coordinate)
• Attributes/properties (e.g. sand content, carbon content, N, P, K, etc.)

5.1 Interpolation

We want to interpolate data to fill the gaps where we do not have observations, the way to do this is creating and empty raster and then filling it with the new data.

# install.packages(gstat)
library(gstat)

# Determination of the extension (xmin, xmax, ymin, ymax)
xmin <- min(soil_data$x) - 100
xmax <- max(soil_data$x) + 100
ymin <- min(soil_data$y) - 100
ymax <- max(soil_data$y) + 100

**Be careful with the variable cellsize, it determines how many calculations we will do to populate it.

# Determination of the resolution (cellsize)
cellsize = 5000

# Spanning a grid (raster format) for interpolation
grd <- expand.grid(x=seq(from = xmin, to = xmax, by = cellsize), 
                   y=seq(from = ymin, to = ymax, by = cellsize))
coordinates(grd) <- ~ x+y
gridded(grd) <- TRUE

5.1.1 IDW

Inverse distance weighted (IDW) interpolation explicitly makes the assumption that things that are close to one another are more alike than those that are farther apart. To predict a value for any unmeasured location, IDW uses the measured values surrounding the prediction location. The measured values closest to the prediction location have more influence on the predicted value than those farther away. IDW assumes that each measured point has a local influence that diminishes with distance. It gives greater weights to points closest to the prediction location, and the weights diminish as a function of distance, hence the name inverse distance weighted. Inverse distance weighted (IDW) interpolation determines cell values using a linearly weighted combination of a set of sample points.

source: How IDW works

Gstat library has for this the function

dataset.idw <- idw(soil_data$clay ~ 1, location = ~x+y, soil_data, grd)

## [inverse distance weighted interpolation]

image(dataset.idw["var1.pred"], col=topo.colors(20))
points(soil_data$x, soil_data$y)
title("IDW Interpolation of clay content")

5.1.2 Kriging

Kriging is an advanced geostatistical procedure that generates an estimated surface from a scattered set of points with z-values. Kriging is based on statistical models that include autocorrelation—that is, the statistical relationships among the measured points. Because of this, it not only have the capability of producing a prediction surface but also provide some measure of the certainty or accuracy of the predictions.

Kriging assumes that the distance or direction between sample points reflects a spatial correlation that can be used to explain variation in the surface. The Kriging tool fits a mathematical function to a specified number of points, or all points within a specified radius, to determine the output value for each location. Kriging is a multistep process:

• exploratory statistical analysis of the data,
• variogram modeling,
• creating the surface, and
• (optionally) exploring a variance surface

source: How Kriging works

# Variogram
v <- variogram(clay ~ 1, locations = ~ x + y, data = soil_data, width = cellsize)
plot(v)

# Variogram fit
v.fit <- fit.variogram(v, fit.method = TRUE, model = vgm(100, "Sph", 50000, 60))

# Output of the variogram
plot(v, model = v.fit)

5.1.2.1 Applying the model fit with the krige() function.

# ordinary kriging:
z <- krige(formula = clay ~ 1, 
           locations = ~ x + y, 
           data = soil_data, 
           newdata = grd, 
           model = v.fit, 
           nmax = 500)

## [using ordinary kriging]

image(z, col = topo.colors(20))
points(soil_data$x, soil_data$y)
title("Ordinary kriging Interpolation of clay content") 

# simple kriging:
y <- krige(clay ~ 1, 
           locations = ~ x + y, 
           data = soil_data, 
           newdata = grd, 
           model = v.fit, 
           nmax = 500, 
           beta = mean(soil_data$clay))

## [using simple kriging]

image(y, col = topo.colors(20))
points(soil_data$x, soil_data$y)
title("Simple kriging Interpolation of clay content") 

Exercise:

• Tests Simple Kriging in addition to Ordinary Kriging,
• Differences theoretically and in visual evaluation?
• Interpolate OC of the soil_data.txt using ordinary kriging
  • What model did you use?
  • What sill?
  • What range?
  • What nugget?

5.1.2.2 Model validation

To validate the data, is a good idea to separate the data in train and test data. But…

• In how many parts should I divide the data?
• Which of those divisions should I use?

An easy way to solve this is using k-fold-cross-validation, in this case, we will use 10-fold

x <- krige.cv(clay ~ 1, 
              locations = ~ x+y, 
              data = soil_data, 
              nfold = 10)

cor(x$observed, x$var1.pred) * cor(x$observed, x$var1.pred) # Rsquared

## [1] 0.04239492

sqr_error <- (x$var1.pred - x$observed )^2
sqrt(sum(sqr_error)/length(sqr_error))                  # RMSE

## [1] 10.29546

5.2 Feature extraction

We can also extract data to our points for some already existing raster

# install.packages("raster")
library(raster)

# Load a raster stack from an external files
terrain.lst <- list.files("./data/Example_Germany/SAGA", pattern="\\.sdat$", full.names = TRUE)
terrain.raster = raster::stack(terrain.lst)

# plot to see SAGA is working well
plot(terrain.raster$DEM)
plot(terrain.raster)

# set projection of the data points (transformed to spatial object)
proj4string(soil_data_sp) <- CRS("+proj=utm +zone=32 +datum=WGS84 +units=m +no_defs")

# plot the point on the raster
plot(terrain.raster$DEM, main = "DEM + Sample Points")
plot(soil_data_sp, add = TRUE, pch = 1)

# extract values
terrain_cov = raster::extract(terrain.raster,
                              soil_data_sp,
                              method = 'simple')

# combing the covariates + soil data
cov_soil = cbind(data.frame(terrain_cov), ID = soil_data_sp$OBJECTID)

# inspect the covariates + soil data 
str(cov_soil)

## 'data.frame':    346 obs. of  16 variables:
##  $ Aspect                    : num  0.429 6.025 6.283 0.501 1.102 ...
##  $ Catchment.Slope           : num  0.0403 0.0149 0.0315 0.0234 0.0269 ...
##  $ Channel.Network.Base.Level: num  399 430 589 409 535 ...
##  $ Channel.Network.Distance  : num  39.86 62.44 13.16 62.79 7.74 ...
##  $ Convergence.Index         : num  8.07 -7.47 -3.48 17.59 6.92 ...
##  $ DEM                       : num  439 493 602 471 537 ...
##  $ LS.Factor                 : num  2.032 0.792 2.588 1.472 0.484 ...
##  $ Modified.Catchment.Area   : num  7135383 14158001 4765989 5804979 13978788 ...
##  $ Plan.Curvature            : num  9.09e-06 -4.76e-06 1.50e-07 1.17e-05 5.02e-06 ...
##  $ Profile.Curvature         : num  1.53e-05 1.49e-06 -5.98e-06 1.24e-06 -4.28e-06 ...
##  $ Relative.Slope.Position   : num  0.03068 0.05149 0.0116 0.053 0.00666 ...
##  $ Slope                     : num  0.04671 0.01737 0.01944 0.02425 0.00935 ...
##  $ Topographic.Wetness.Index : num  9.21 11.06 13.54 11.19 12.46 ...
##  $ Total.Catchment.Area      : num  309775 673476 7381662 1190722 1618878 ...
##  $ Valley.Depth              : num  1260 1150 1122 1122 1155 ...
##  $ ID                        : int  1 2 3 4 5 6 7 8 9 10 ...

Aspect	Catchment.Slope	Channel.Network.Base.Level	Channel.Network.Distance	Convergence.Index	DEM	LS.Factor	Modified.Catchment.Area	Plan.Curvature	Profile.Curvature	Relative.Slope.Position	Slope	Topographic.Wetness.Index	Total.Catchment.Area	Valley.Depth	ID
0.4294031	0.0403004	399.1916	39.8607788	8.0702753	439.0524	2.0323753	7135383	9.10e-06	0.0000153	0.0306755	0.0467074	9.210298	3.097750e+05	1259.5742	1
6.0251021	0.0148807	430.4036	62.4364319	-7.4680462	492.8400	0.7922727	14158001	-4.80e-06	0.0000015	0.0514949	0.0173740	11.057698	6.734757e+05	1150.0422	2
6.2830095	0.0314850	589.4273	13.1578979	-3.4751248	602.3657	2.5878420	4765989	1.00e-07	-0.0000060	0.0115953	0.0194362	13.540216	7.381662e+06	1121.6007	3
0.5011653	0.0234480	408.5652	62.7922974	17.5941334	471.3575	1.4717088	5804979	1.17e-05	0.0000012	0.0529975	0.0242488	11.189030	1.190722e+06	1122.0229	4
1.1021508	0.0268809	534.5034	7.7376099	6.9179778	537.1065	0.4840871	13978788	5.00e-06	-0.0000043	0.0066573	0.0093498	12.459466	1.618878e+06	1154.5442	5
6.1787410	0.0085865	413.6719	52.1421204	15.2397661	465.8141	0.1153976	18593388	8.10e-06	0.0000021	0.0430102	0.0034977	13.053521	8.969104e+05	1160.1765	6
1.5702554	0.0115924	524.0745	0.0000000	-20.8493633	518.6011	0.1907610	13070727	-9.20e-06	-0.0000160	0.0000000	0.0063188	11.796611	4.198033e+05	1194.5505	7
0.0219634	0.0044136	403.9760	3.9115295	-56.5708809	407.8875	0.1198936	23013878	-6.30e-06	-0.0000059	0.0032597	0.0024526	14.657634	2.908450e+06	1196.0389	8
2.9659510	0.0194990	419.6031	67.5102539	-12.2554197	487.1133	0.6217782	21073264	-1.60e-06	0.0000034	0.0558615	0.0180725	10.284395	3.066954e+05	1141.0183	9
6.2831855	0.0352127	524.0759	0.0000000	-2.5789902	518.5703	0.0000493	14481057	-1.50e-06	-0.0000100	0.0000000	0.0000020	19.155371	1.042419e+08	1191.9161	10
4.0093184	0.0060658	421.4517	86.6167603	1.1484684	508.0685	0.1281908	6967646	-1.50e-06	0.0000047	0.0745009	0.0060658	10.976505	2.500000e+05	1076.0103	11
4.2115779	0.0248918	422.4179	80.1736145	-11.2577705	502.5915	0.8086568	5553020	2.00e-07	0.0000033	0.0687277	0.0246467	9.622850	2.527281e+05	1086.3665	12
5.6866570	0.0113736	522.6282	11.1656494	-39.8165398	532.0483	0.2205232	9518352	-7.70e-06	-0.0000008	0.0096513	0.0056647	12.623464	1.194467e+06	1145.7350	13
5.1019182	0.0301106	464.3504	100.1460876	-14.4570904	564.4965	2.2024078	7363513	-1.12e-05	0.0000024	0.0850035	0.0434108	9.722178	4.728093e+05	1077.9946	14
6.2299914	0.0195076	419.4833	75.0160217	5.0001016	494.4993	0.7321568	12600115	1.38e-05	0.0000149	0.0631425	0.0195741	10.353709	3.229835e+05	1113.0269	15
1.6470327	0.0326655	485.6871	9.7435608	5.1283569	494.5802	2.3675339	11509166	5.90e-06	0.0000099	0.0083623	0.0128723	15.069006	2.419574e+07	1155.4279	16
1.7217468	0.0205693	438.1074	26.3657532	-14.5889292	464.4732	0.7058899	10013359	-1.48e-05	-0.0000056	0.0227902	0.0207303	10.018474	2.649307e+05	1130.5249	17
1.8696905	0.0130822	468.4013	5.8935547	-2.1560748	474.2948	0.2817954	10269242	-5.40e-06	-0.0000040	0.0050130	0.0084737	11.524859	5.360345e+05	1169.7618	18
5.8745556	0.0269382	459.4295	52.4487915	4.4881229	502.9655	0.6585216	7647716	2.90e-06	0.0000114	0.0452288	0.0100270	12.931605	2.724024e+06	1107.1836	19
6.2821155	0.0231248	467.0049	0.0000000	4.5383015	463.6544	0.1422335	16389276	2.70e-06	0.0000056	0.0000000	0.0045064	12.499353	6.048746e+05	1188.0206	20
5.6739435	0.0424227	441.4896	57.1219177	-2.3343239	498.6115	9.6078730	9581768	-9.50e-06	-0.0000066	0.0491317	0.0849568	10.552262	2.268537e+06	1105.5066	21
3.8055081	0.0276956	474.5028	10.0762024	13.6967726	484.5790	2.8710864	14130641	3.00e-07	0.0000012	0.0087176	0.0353771	11.254631	1.918667e+06	1145.7643	22
1.3476894	0.0233555	475.4897	27.0938110	-0.3060751	502.5835	0.8713754	11599377	-3.20e-06	-0.0000040	0.0236048	0.0059324	15.862392	2.748267e+07	1120.7131	23
2.3694320	0.0082407	466.9540	0.0000000	-44.0294876	456.4023	0.0003811	31156224	-3.00e-07	0.0000002	0.0000000	0.0000424	14.365612	1.225612e+06	1218.0914	24
0.7568264	0.0273188	432.7445	60.5001221	-18.0090752	493.2447	1.5036186	8718220	-6.50e-06	0.0000028	0.0526677	0.0339331	9.814684	4.390969e+05	1088.2137	25
6.1815634	0.0140102	480.5701	14.8707581	4.5279641	495.4409	1.2535824	14251032	4.30e-06	0.0000005	0.0127638	0.0176535	12.137032	1.806407e+06	1150.2050	26
0.1670526	0.2397703	484.6573	9.3578491	3.1012073	494.0151	1.4383816	19979734	2.71e-05	-0.0000657	0.0077639	0.0254288	10.929847	8.179902e+05	1195.9423	27
0.4116030	0.0116828	486.7178	0.0009155	-20.2330170	483.8961	0.5170770	17571940	-6.80e-06	-0.0000062	0.0000008	0.0030708	17.356234	6.972422e+07	1156.4420	28
1.5935348	0.0975181	488.9113	0.9034424	-7.4636097	489.8148	1.2654620	38850504	-8.70e-06	-0.0000470	0.0007602	0.0138468	13.192804	3.798083e+06	1187.4652	29
1.0660473	0.0063044	450.6757	7.6308289	4.3924532	458.3065	0.1487880	43271404	4.60e-06	0.0000032	0.0066897	0.0054501	11.803857	4.953530e+05	1133.0518	30
0.9269978	0.0280522	498.1309	1.3374023	-22.2997589	499.4683	2.1418457	35491748	-7.20e-06	-0.0000025	0.0011556	0.0129917	14.779296	2.384370e+07	1155.9470	31
2.1010513	0.0211931	417.1243	13.0282593	0.5609143	430.1526	0.8065014	25194818	9.00e-07	-0.0000059	0.0112767	0.0190581	10.709024	5.837075e+05	1142.3015	32
2.4854426	0.0043096	482.8287	42.5517883	24.0456352	525.3805	0.0734831	42068284	1.54e-05	0.0000077	0.0380233	0.0038826	11.481558	2.638430e+05	1076.5466	33
5.9501805	0.0062135	511.0038	17.1245728	-2.0556815	528.1284	0.1258160	40916040	4.00e-07	0.0000009	0.0155339	0.0019346	15.786536	8.831474e+06	1085.2739	34
2.0676687	0.0240553	425.6458	0.0000000	0.4430169	425.4696	1.5522254	17004296	-1.30e-06	-0.0000241	0.0000000	0.0022404	21.444437	3.123710e+09	1150.8246	35
1.4754837	0.0148999	686.2984	75.8235474	3.5629139	762.1219	0.7067265	9036964	7.00e-06	0.0000123	0.0773636	0.0167669	10.923225	5.070329e+05	904.2698	36
1.2833834	0.0059302	388.2061	0.0007935	-6.4232211	351.1060	0.0005741	27207434	0.00e+00	0.0000000	0.0000007	0.0000051	22.297335	2.998190e+09	1155.4734	37
3.8720543	0.0312515	417.2045	62.6257019	-6.9626751	479.8302	1.6006013	6244348	-1.70e-06	-0.0000060	0.0547819	0.0357654	9.747457	4.322598e+05	1080.5574	38
0.4790111	0.0033186	445.8684	1.6858826	-24.3723431	447.5543	0.0424801	34933520	-1.50e-06	0.0000001	0.0014810	0.0019651	13.005647	5.894277e+05	1136.6678	39
1.2732722	0.0047902	468.0482	0.0000000	9.2430019	468.0482	6.8092055	44678808	1.90e-06	0.0000003	0.0000000	0.0109428	18.400211	6.696788e+08	1126.5544	40
4.7137332	0.0209352	388.2499	0.0040894	-16.9162006	374.5143	0.3417995	29454406	-1.00e-06	-0.0000019	0.0000035	0.0091248	11.692849	5.469284e+05	1162.8075	41
4.5832558	0.0162871	444.2717	38.1552124	-13.6213560	482.4269	0.2330286	18815832	-7.50e-06	0.0000158	0.0335328	0.0057110	12.726735	1.077041e+06	1099.6927	42
1.4250078	0.0273865	433.4230	55.2590332	-10.4666681	488.6821	1.2224814	9518664	-7.20e-06	0.0000097	0.0484272	0.0298781	9.838017	3.176264e+05	1085.8148	43
0.0005361	0.0414970	396.9891	21.2042847	-0.5858462	414.6100	0.2771052	8894298	-6.40e-06	-0.0000288	0.0185730	0.0037567	14.939976	5.785755e+06	1120.4670	44
0.0001652	0.0227796	388.2372	0.0000000	-12.7527924	383.9200	1.7534391	11072936	-4.90e-06	0.0000038	0.0000000	0.0086798	15.993131	3.830814e+07	1154.3915	45
2.0413733	0.0281339	413.9184	82.8390198	8.6889496	496.7574	0.9549094	6643826	0.00e+00	-0.0000034	0.0719279	0.0279688	9.501096	2.515699e+05	1068.8556	46
5.8971243	0.0079005	483.4892	5.9186401	0.3893614	489.4079	0.0816759	42919912	1.30e-06	0.0000044	0.0052979	0.0027525	13.207749	9.765109e+05	1111.2585	47
0.1971965	0.0277097	407.0197	85.8222656	1.6323705	487.8000	6.7820182	4512446	3.20e-06	-0.0000023	0.0762453	0.0420278	12.671556	7.880555e+06	1039.7849	48
0.6150197	0.0306211	389.7137	10.5391846	3.2693100	400.2528	1.6140084	16713239	4.50e-06	0.0000007	0.0090253	0.0210042	12.030229	2.455747e+06	1157.1964	49
4.6468682	0.0439678	600.0666	8.6245728	-0.5419587	608.6912	0.6536893	13840006	1.10e-06	-0.0000223	0.0082088	0.0129327	11.831675	9.457600e+05	1042.0239	50
0.8531815	0.0178096	586.4072	61.3285522	39.7424126	647.7358	0.7001847	5545159	2.57e-05	-0.0000034	0.0581078	0.0182983	10.528520	4.824184e+05	994.0991	51
3.1415927	0.0161558	388.2705	0.0010376	22.0138817	372.0182	0.0000174	23788616	1.40e-06	-0.0000021	0.0000009	0.0000020	16.603088	8.120818e+06	1157.4319	52
0.9549264	0.0032558	388.2670	0.0005188	-11.2657652	366.8014	0.0001071	26415268	0.00e+00	0.0000000	0.0000004	0.0000024	20.482401	5.477823e+08	1161.7285	53
0.7908035	0.0201335	388.2937	0.0000000	-3.3175592	382.1389	0.0000271	19779832	4.10e-06	-0.0000097	0.0000000	0.0000040	15.432438	3.562059e+06	1170.9193	54
3.7948935	0.0320218	402.9882	36.8593750	-0.7061846	439.8476	0.4344137	9834538	2.10e-06	0.0000024	0.0317436	0.0086067	12.540752	1.686203e+06	1124.2997	55
2.2145631	0.0230064	452.5968	57.8683472	44.7065887	510.4652	0.7489961	8698146	2.84e-05	0.0000036	0.0513155	0.0232286	9.683096	2.609337e+05	1069.8295	56
3.1954527	0.0341749	469.3592	31.1771240	28.7386494	500.5363	2.1179512	19923436	4.90e-06	0.0000031	0.0277658	0.0397158	10.002428	4.616741e+05	1091.6840	57
0.5890409	0.0476098	408.5839	97.4269409	3.4405394	506.0108	1.8779331	6207623	2.80e-06	0.0000033	0.0841799	0.0476098	8.931397	2.500000e+05	1059.9392	58
0.0528528	0.0193999	465.2445	0.0000000	-42.3261871	465.2445	4.3802171	10884138	-5.40e-06	-0.0000170	0.0000000	0.0096105	17.849037	2.852653e+08	1120.9910	59
1.5707334	0.0364971	414.6162	33.3876038	4.6048617	433.4400	4.6491699	3755360	-8.00e-07	-0.0000216	0.0296878	0.0223163	14.417646	2.037957e+07	1091.2374	60
3.0179920	0.0313872	396.8251	14.5644836	4.1163611	411.3896	1.0630966	13623019	9.00e-06	-0.0000123	0.0127048	0.0176634	11.722580	1.215272e+06	1131.8125	61
0.8523530	0.0214600	540.3507	27.7643433	4.9592800	568.1151	0.8079782	25420144	7.00e-07	0.0000009	0.0257892	0.0223287	10.040517	3.613986e+05	1048.8219	62
3.9086089	0.0107538	601.9206	27.5477295	6.3307309	629.4683	0.1959352	7421885	8.00e-07	-0.0000130	0.0254269	0.0084150	10.645915	2.500000e+05	1055.8628	63
0.0305921	0.0228955	476.1695	0.0000000	-36.7569733	476.1695	5.3369660	13738657	-6.00e-06	-0.0000058	0.0000000	0.0131961	16.995443	1.634383e+08	1114.3616	64
6.0818453	0.0178226	484.8479	16.0131531	-13.8505335	500.8610	1.0560383	15485285	-7.30e-06	-0.0000129	0.0144283	0.0156552	12.218871	1.870912e+06	1093.8268	65
2.6736531	0.0034597	388.3199	0.0005798	37.6333923	379.2700	0.0000986	26786464	0.00e+00	0.0000000	0.0000005	0.0000032	19.402843	1.793789e+08	1165.4951	66
0.5840763	0.0161775	390.7841	2.8949585	5.2615285	393.6791	0.3283078	22153760	6.00e-07	-0.0000001	0.0024724	0.0053613	13.852303	3.853674e+06	1168.0162	67
5.8738446	0.0127102	613.8002	45.6375122	-6.5096745	659.4377	1.6547922	7550458	-8.00e-07	0.0000030	0.0420973	0.0110028	14.840498	2.016927e+07	1038.4586	68
0.7003326	0.1181968	635.1302	0.0000000	-0.8311540	632.8980	0.0007826	23842768	6.60e-06	-0.0000584	0.0000000	0.0000168	19.179243	1.504359e+08	1155.2916	69
2.0382004	0.0902884	409.3395	80.3335571	-59.0936546	489.6731	4.6992917	4056355	-6.78e-05	-0.0000417	0.0710538	0.0934099	8.361517	2.692564e+05	1050.2677	70
6.0943661	0.0435392	589.6236	0.0000000	-3.4281051	587.6093	3.7366691	18765156	2.90e-06	0.0000014	0.0000000	0.0191134	14.529815	2.284220e+07	1098.6371	71
3.7602031	0.0832314	404.7939	38.2098389	-3.0003304	443.0038	18.0230598	3125138	-9.90e-06	-0.0000024	0.0335407	0.1738538	9.086372	1.081853e+06	1100.9983	72
4.2385325	0.0395057	398.3005	93.1581421	2.0989809	491.4586	1.5276017	4076094	3.40e-06	-0.0000130	0.0810108	0.0382859	9.341374	2.939118e+05	1056.7888	73
5.1658235	0.0145213	550.2671	47.8892212	15.3622417	598.1563	0.6454003	16497744	7.80e-06	-0.0000043	0.0448925	0.0097871	12.984212	2.849313e+06	1018.8635	74
3.4370072	0.0009616	611.8246	161.8554077	35.2689362	773.6800	0.0122794	5694243	3.80e-05	0.0000380	0.1531047	0.0009616	12.900967	2.500000e+05	895.2998	75
6.2802334	0.0239466	628.9435	72.3813477	-59.5168915	680.8800	0.0262455	13349889	-5.30e-06	-0.0000120	0.0703785	0.0004756	17.088171	1.322959e+07	956.0774	76
1.2139363	0.0359554	720.8580	25.7457275	0.1891483	746.6038	1.4523233	6637907	1.40e-06	-0.0000250	0.0262998	0.0126025	13.937304	9.155820e+06	953.1859	77
1.1512845	0.0107586	678.9570	152.3376465	17.1091232	831.2947	0.1657996	10584692	3.90e-06	-0.0000043	0.1532929	0.0060936	11.600197	4.390611e+05	841.4310	78
4.1806130	0.0721808	417.8424	98.2531433	1.5907160	516.0955	2.6669695	2848586	-1.85e-05	0.0000217	0.0896971	0.0351832	11.093617	1.584074e+06	997.1346	79
3.1480408	0.0412393	401.7710	20.3065796	12.8577614	422.0776	1.9706973	15999844	2.70e-06	0.0000011	0.0183107	0.0430846	9.476317	2.830451e+05	1088.6968	80
3.6273131	0.0440367	409.1613	77.8714600	16.9860001	487.0328	1.5980138	5413029	2.84e-05	-0.0000097	0.0703396	0.0404070	9.224892	2.771390e+05	1029.2064	81
0.2529656	0.0370695	400.8728	14.5579529	17.0072136	415.4308	2.2298334	14481657	1.12e-05	0.0000016	0.0130609	0.0415978	9.934386	5.230562e+05	1100.0636	82
4.3111601	0.0023183	401.6523	0.0175476	-15.1165419	363.4602	0.0000807	26413130	0.00e+00	-0.0000009	0.0000158	0.0000022	20.116095	3.572110e+08	1112.1362	83
2.9520261	0.0276420	478.9229	80.9860535	26.5267811	559.9090	1.0947368	4782128	2.52e-05	-0.0000026	0.0788607	0.0289302	9.699101	2.761134e+05	945.9647	84
4.2242775	0.0274940	493.9069	26.4050598	5.0188756	520.3120	1.8630586	7332767	6.00e-06	-0.0000012	0.0262747	0.0374291	9.933848	5.219535e+05	978.5579	85
4.3803592	0.0623979	374.8027	236.9539185	0.0085587	611.7567	3.1847773	3977271	-9.10e-06	0.0000099	0.2067550	0.0664151	8.837606	2.912484e+05	909.1075	86
4.3153033	0.0526704	410.9751	17.6630554	-8.8123970	428.6381	1.3459315	10883193	-2.54e-05	-0.0000099	0.0160940	0.0257030	10.718192	7.600021e+05	1079.8286	87
0.1453524	0.0417146	392.6238	114.4607849	-5.8416982	507.0846	6.4909997	2782601	1.80e-06	0.0000054	0.1017079	0.0723394	10.254725	1.167886e+06	1010.9268	88
2.9529188	0.0248701	322.1551	0.0000000	5.7107706	322.1551	0.5241104	34301732	1.54e-05	-0.0000033	0.0000000	0.0027360	17.880693	9.325815e+07	1205.3911	89
4.4562473	0.0335179	486.8655	129.2628174	3.7106555	616.1284	0.8591496	7144441	7.90e-06	0.0000000	0.1285708	0.0196506	10.737001	5.521281e+05	876.1195	90
1.7717868	0.0132856	325.8745	13.5824890	-2.6707351	339.4570	0.3794811	31873556	-1.00e-06	-0.0000010	0.0112689	0.0132332	10.374478	2.500000e+05	1191.7302	91
4.5207157	0.0368364	496.4612	15.2979736	4.0287919	511.7592	1.6599982	9477036	3.50e-06	-0.0000097	0.0153709	0.0332564	10.147640	4.978135e+05	979.9561	92
5.0104370	0.0314421	455.2575	60.1119080	-1.5332880	515.3694	3.6783929	11772932	6.90e-06	-0.0000024	0.0576308	0.0316689	12.344631	4.546881e+06	982.9407	93
0.3361260	0.0327392	405.7759	76.9106445	4.7426543	482.6866	0.1330874	3526524	-1.00e-07	0.0000089	0.0699541	0.0046920	12.161665	5.717404e+05	1022.5332	94
4.1626000	0.0093672	408.6243	118.4683838	-1.7119368	527.0927	0.2118878	2087112	-1.24e-05	-0.0000199	0.1082225	0.0087788	10.661747	2.576921e+05	976.2059	95
3.2094927	0.0301643	361.1405	156.9351196	2.0339043	518.0756	1.1272383	3609383	1.60e-05	0.0000076	0.1353440	0.0292287	9.728622	2.615544e+05	1002.5929	96
5.5194273	0.0312128	489.5268	10.1972046	0.8399060	479.3717	1.3567219	6630481	2.40e-06	-0.0000002	0.0101651	0.0187352	12.081940	2.339984e+06	992.9619	97
3.6025903	0.0619960	381.4140	65.9026794	57.5214005	447.3167	2.8705494	3477814	3.26e-05	0.0000495	0.0585944	0.0540461	9.453472	4.623801e+05	1058.8241	98
2.2358391	0.0337043	444.1265	2.0495605	-3.4112990	446.1760	1.1464005	9634529	-6.60e-06	-0.0000091	0.0019504	0.0090660	14.745733	1.613433e+07	1048.7809	99
1.6605543	0.0460184	406.7386	56.2995300	5.2219858	463.0381	2.2289720	6430317	1.40e-06	-0.0000103	0.0515108	0.0482354	9.304350	2.878563e+05	1036.6650	100
1.4652292	0.0235512	335.8281	7.3266602	2.0780013	343.1548	1.2626098	21988250	2.90e-06	-0.0000031	0.0061316	0.0208829	11.441004	1.068818e+06	1187.5719	101
4.5130167	0.0405952	377.4371	11.6819458	16.5859966	389.1190	2.1339903	9863084	5.90e-06	-0.0000010	0.0103064	0.0360968	10.427298	7.184253e+05	1121.7780	102
4.8275595	0.0743301	372.8984	20.1931152	4.8417773	393.0915	4.9963055	7598710	1.98e-05	0.0000331	0.0176629	0.0503927	11.136347	1.917842e+06	1123.0583	103
4.4047194	0.0537709	337.0520	6.5729980	-3.2103639	343.6250	1.1092265	9753377	-3.90e-06	-0.0000438	0.0055472	0.0190060	11.517440	1.199007e+06	1178.3563	104
1.2178153	0.0622412	340.2065	0.0000000	-8.6642761	340.2065	7.0708871	4730678	-8.00e-06	-0.0000446	0.0000000	0.0340615	13.668903	1.889344e+07	1188.5195	105
1.8477196	0.0380714	387.1126	12.8154297	-2.6654868	399.9280	8.6854477	4046165	9.80e-06	-0.0000001	0.0114076	0.0583026	11.899248	5.305176e+06	1110.5918	106
1.2184675	0.0253073	404.6951	35.4920654	25.3002968	440.1872	1.2144966	6750660	3.80e-06	0.0000080	0.0324278	0.0224436	11.037566	8.956504e+05	1059.0024	107
2.8426521	0.0192208	465.3434	19.1777344	-2.1798112	484.5212	0.4027589	27839356	-1.30e-06	0.0000028	0.0186549	0.0136481	10.392117	2.781452e+05	1008.8470	108
0.6488651	0.0215986	440.0932	45.7728577	14.0333395	485.8661	0.7464559	3671045	1.25e-05	-0.0000069	0.0421589	0.0227694	9.759458	2.762670e+05	1039.9489	109
3.4481411	0.0056920	468.5590	29.6457520	7.8410149	498.2047	0.1389711	26051622	7.70e-06	0.0000070	0.0289381	0.0062283	11.065987	2.500000e+05	994.8069	110
3.7625690	0.0315569	423.3007	0.0000000	-0.9709566	420.0895	2.6545441	8637834	-3.00e-06	-0.0000131	0.0000000	0.0097673	16.528202	1.026792e+08	1072.9093	111
3.6095321	0.0307710	409.9495	0.0000000	-14.0851908	398.5285	0.0000857	11944316	-1.18e-05	-0.0000075	0.0000000	0.0000032	19.052301	1.263375e+08	1090.6455	112
0.0348955	0.0247954	404.2471	0.0000000	-4.8464699	396.4838	0.5811085	9848764	-7.20e-06	-0.0000108	0.0000000	0.0108968	12.265389	1.195958e+06	1102.2064	113
1.5708560	0.0346545	404.3028	3.9148254	-6.7318835	408.2177	1.1876347	19213882	-3.20e-06	-0.0000054	0.0035622	0.0240340	10.690697	5.282218e+05	1095.0618	114
3.2958975	0.0300608	444.4190	71.1498413	-12.7837257	515.5688	4.5189872	9995168	-9.30e-06	0.0000006	0.0658859	0.0470359	11.178219	1.922886e+06	1008.7457	115
2.9281099	0.0284939	455.5407	108.3299255	54.3895569	563.8707	1.2040021	4656641	3.99e-05	0.0000194	0.1015415	0.0275648	10.142408	4.163595e+05	958.5243	116
0.0567865	0.0366037	427.1938	96.4570923	-11.2017269	523.6509	4.6785631	5864691	-5.40e-06	0.0000080	0.0901852	0.0303728	13.123512	8.023586e+06	973.0876	117
5.2879939	0.0174381	407.8834	80.9165649	-8.4786177	488.8000	0.4434487	4892179	-8.00e-07	0.0000070	0.0720263	0.0080095	12.897839	2.212739e+06	1042.5140	118
3.3765967	0.0717517	428.6212	102.3158875	-0.0010124	530.9371	5.1666064	3365104	1.03e-05	0.0000132	0.0948019	0.0783686	9.344182	5.410609e+05	976.9445	119
0.5063162	0.0376331	433.7624	103.7931519	10.9044790	537.5555	3.0222921	6451252	2.13e-05	0.0000120	0.0959827	0.0451735	10.344212	9.549064e+05	977.5803	120
5.5479264	0.0599544	434.4426	109.4284058	41.9412766	543.8710	3.7374985	4655294	2.17e-05	0.0000118	0.1005499	0.0732142	8.823666	3.518684e+05	978.8710	121
3.9702454	0.0424820	422.9643	153.7374573	43.9950485	576.7017	2.9023120	3421969	6.60e-06	0.0000076	0.1413909	0.0519884	9.645906	5.682516e+05	933.5848	122
1.4276159	0.0298516	418.5026	14.9915771	9.6874886	433.4942	1.2518033	17844076	7.50e-06	-0.0000203	0.0136556	0.0267587	10.365876	4.812863e+05	1082.8456	123
1.7424270	0.0211194	409.4274	4.5797424	5.1414452	414.0071	1.9081146	12163841	-4.40e-06	0.0000009	0.0041302	0.0110832	15.165638	2.471730e+07	1104.2616	124
2.1385789	0.0661647	265.1468	28.2730713	-4.4210529	293.4199	2.2873418	12401912	-1.26e-05	-0.0000010	0.0230631	0.0465212	9.522735	4.392934e+05	1197.6267	125
6.1085801	0.0440525	295.5614	0.0000000	-39.5098190	285.6800	32.2911530	4622434	-7.70e-06	-0.0000431	0.0000000	0.0558722	15.363013	1.522619e+08	1194.9082	126
3.7967129	0.0317779	231.5043	0.7257233	-13.5153656	232.2300	0.9396314	15833789	-1.16e-05	-0.0000056	0.0005755	0.0132297	12.642322	2.869813e+06	1260.2078	127
1.7952821	0.0326481	242.8164	0.0000000	-12.0239391	242.8164	17.2374649	15428531	-9.90e-06	-0.0000108	0.0000000	0.0308873	16.312361	2.246658e+08	1249.0038	128
2.1168747	0.0628870	290.6330	37.4869995	-1.1895900	328.1200	1.2518181	22604114	3.80e-06	-0.0000083	0.0312753	0.0196378	11.680795	1.595810e+06	1161.1270	129
2.6960533	0.0261016	254.9940	192.3505707	47.9591827	447.3445	1.0027084	4430118	4.10e-05	0.0000244	0.1545891	0.0270267	9.768824	3.150288e+05	1051.9192	130
0.0002569	0.0291693	303.2838	0.0000000	-3.6019113	293.4673	1.9821750	13054683	9.20e-06	-0.0000034	0.0000000	0.0077198	16.797785	7.618762e+07	1187.2191	131
1.6191441	0.0215267	236.9980	119.8960266	-9.5276499	356.8940	0.7818194	5216190	-1.38e-05	-0.0000118	0.0955139	0.0215483	10.109425	2.772874e+05	1135.3776	132
0.3104938	0.0228665	273.9884	106.3699341	28.7036915	380.3583	0.7550328	6108114	9.30e-06	0.0000152	0.0876253	0.0228429	9.774317	2.525280e+05	1107.5483	133
2.8585732	0.0379800	273.3537	32.1662598	-5.6642985	305.5200	1.5151381	11146017	-1.00e-05	-0.0000157	0.0264557	0.0221982	11.637231	1.558051e+06	1183.6876	134
4.0372620	0.0214835	215.7539	29.7543640	-2.6442590	245.5083	1.0824323	11657446	-1.90e-06	-0.0000002	0.0232831	0.0264659	10.049191	4.304631e+05	1248.1840	135
4.6909924	0.0124395	222.9895	19.4857178	-11.8901978	242.4752	0.3709500	47113604	6.00e-07	0.0000054	0.0152784	0.0124395	10.580486	2.500000e+05	1255.8889	136
4.6437817	0.0142128	208.6369	23.4798126	-2.1091585	232.1167	0.9771010	31040080	-1.90e-06	-0.0000061	0.0182247	0.0123111	13.045906	3.040131e+06	1264.8748	137
3.6733911	0.0380682	238.6277	55.5360260	9.3635521	294.1637	1.2658331	13819316	-1.40e-06	-0.0000057	0.0439747	0.0346727	9.292687	2.578053e+05	1207.3739	138
5.9065123	0.0056303	218.1501	32.1417084	5.5673437	250.2918	0.0883821	27481816	2.50e-06	-0.0000044	0.0250655	0.0041231	11.687637	3.186059e+05	1250.1666	139
3.0014117	0.0275107	227.5256	53.9984741	-13.5665007	281.5240	1.1426973	7715975	-1.64e-05	0.0000011	0.0426984	0.0256235	10.322105	4.401317e+05	1210.6495	140
0.5875238	0.0663242	310.3238	16.8740540	-18.8334122	327.1978	12.8942995	68954064	-1.50e-05	-0.0000122	0.0143312	0.0469105	13.810807	3.239463e+07	1160.5619	141
1.2692132	0.0204048	201.9395	66.8510590	10.3859930	268.7905	0.8433577	18501818	2.80e-06	0.0000091	0.0515954	0.0219706	10.216350	3.761459e+05	1228.8270	142
1.4898387	0.0200451	203.6484	89.6723328	13.7133560	293.3207	0.7629214	6163314	1.87e-05	0.0000243	0.0694381	0.0213627	10.085019	2.760671e+05	1201.7266	143
2.5944436	0.0168814	208.7712	38.4300232	-2.5972137	247.2012	0.4787988	11456358	-2.90e-06	-0.0000042	0.0298944	0.0155344	10.274288	3.093403e+05	1247.0961	144
0.3290379	0.0080563	205.8031	84.7248840	-26.1297226	290.5280	0.0946265	7333121	-7.20e-06	-0.0000017	0.0654843	0.0040484	11.936041	3.923168e+05	1209.0942	145
6.2238998	0.0209877	220.2083	112.3526001	31.8043404	332.5609	0.6548770	4472032	1.52e-05	0.0000101	0.0883204	0.0172693	10.607263	3.691818e+05	1159.7495	146
2.6373351	0.1320393	287.6413	337.6320801	-6.3932304	625.2734	10.4379930	1719525	-3.52e-05	0.0000073	0.2812193	0.1271524	9.047504	7.379823e+05	862.9686	147
1.5599091	0.0714339	271.5806	107.1855774	3.0367730	378.7662	4.7034159	4507146	2.60e-06	-0.0000114	0.0879904	0.0687668	9.664555	5.482442e+05	1110.9658	148
2.3398309	0.0084935	264.5526	43.8287659	-8.6430550	308.3813	0.1301304	8699155	-5.30e-06	-0.0000003	0.0353412	0.0052039	11.665416	4.285315e+05	1196.3306	149
4.5890589	0.0293110	210.0959	89.1183929	9.4945612	299.2143	1.1154171	4735673	1.19e-05	0.0000057	0.0692589	0.0281025	9.869255	3.029919e+05	1197.6238	150
5.6339970	0.0175952	219.7626	65.4285126	0.0631906	285.1911	1.3464341	10841644	2.50e-06	-0.0000002	0.0510569	0.0176551	12.315265	2.759379e+06	1216.0536	151
3.1894400	0.0463924	230.0895	56.6546936	-2.6015365	286.7442	4.0658832	18429222	-1.03e-05	-0.0000680	0.0449331	0.0519577	10.491221	9.797314e+05	1204.2125	152
2.3685338	0.0314742	207.3812	47.9874573	5.8334694	255.3687	1.9298444	3871222	7.10e-06	0.0000147	0.0372848	0.0210890	12.459912	3.844529e+06	1239.0646	153
4.6025524	0.0610740	201.2729	0.0000000	2.9994371	192.4988	6.1138220	8345019	1.12e-05	-0.0000124	0.0000000	0.0496001	11.708340	3.330302e+06	1294.3488	154
NA	NA	NA	NA	NA	NA	NA	NA	NA	NA	NA	NA	NA	NA	NA	155
1.1751279	0.0608584	251.5824	24.6713409	-12.9801855	276.2538	3.1575904	4223857	-1.19e-05	0.0000005	0.0197822	0.0623755	9.082760	3.595861e+05	1222.4780	156
4.2846923	0.0256042	204.0367	40.3245850	-20.9615650	244.3613	0.4580286	5921102	-9.40e-06	-0.0000197	0.0312003	0.0053663	14.680748	8.443778e+06	1252.1184	157
3.6912289	0.0288871	244.6763	72.8766174	-24.1510181	317.5529	1.3431267	4744866	-7.10e-06	0.0000256	0.0580900	0.0284196	10.285993	5.730424e+05	1181.6709	158
4.1501126	0.0513105	158.3926	289.5003967	76.7649689	447.8930	2.7627828	1883121	9.21e-05	0.0001222	0.2154689	0.0554406	9.249574	3.980115e+05	1054.0828	159
5.8252501	0.1153564	156.5117	161.8287811	-6.4937787	318.3405	7.5352988	3001780	6.90e-06	0.0000185	0.1210542	0.1293779	8.159244	3.044944e+05	1175.0005	160
2.5666041	0.0829851	148.6916	152.6333466	15.0731239	301.3250	5.0677505	1702084	2.43e-05	0.0000404	0.1135317	0.0717465	9.670888	7.876975e+05	1191.7778	161
1.1705346	0.0564001	118.3649	3.8913116	-5.4067960	122.2562	0.1592854	12359431	-7.90e-06	-0.0000170	0.0028238	0.0031511	14.302800	3.361553e+06	1374.1332	162
5.9487176	0.0474914	528.7094	23.8737183	6.3447204	552.5831	4.6017065	8334783	2.00e-06	-0.0000229	0.0248403	0.0483768	11.104132	2.047327e+06	937.2140	163
1.7701923	0.0320235	532.5034	51.4166260	16.3120384	583.9200	1.5261294	9094520	2.29e-05	0.0000005	0.0537685	0.0356725	9.639401	3.229007e+05	904.8434	164
3.1417832	0.0266583	496.1202	3.2051086	-11.3500967	489.7627	2.2812116	8407908	2.30e-06	-0.0000022	0.0032366	0.0105722	15.812719	3.895944e+07	987.0516	165
1.9013820	0.0324615	495.8323	68.9540405	-4.2578306	564.7864	1.6219826	8977111	-1.02e-05	-0.0000072	0.0694708	0.0300119	10.525954	7.107464e+05	923.6069	166
4.2017875	0.0792308	373.9078	0.0000000	-0.4418248	373.6087	4.1631083	19747308	7.90e-06	-0.0000076	0.0000000	0.0349951	12.229684	4.879174e+06	1116.7084	167
4.7109509	0.0310972	364.9009	17.7066345	1.6352379	382.6076	2.0791323	11774919	-5.10e-06	-0.0000198	0.0156980	0.0393177	9.998992	4.334241e+05	1110.2504	168
0.0286904	0.0573511	407.4841	199.5451355	3.2786562	607.0292	2.8214767	2058315	5.30e-06	0.0000315	0.1847196	0.0582915	9.089086	2.657762e+05	880.7147	169
1.1264369	0.0373757	355.6315	94.5248108	9.3288994	450.1563	2.2125418	4789232	5.20e-06	-0.0000241	0.0831273	0.0425080	9.822953	5.230106e+05	1042.5840	170
4.7122388	0.0349133	406.6212	18.8512268	-26.1545506	424.6638	0.6445178	14603041	-2.39e-05	-0.0000259	0.0173232	0.0093607	13.170068	2.454969e+06	1069.3579	171
4.1369362	0.0766347	304.3691	118.4113464	7.2265420	422.7805	4.0164046	6966434	3.52e-05	0.0000274	0.0999054	0.0762553	8.830722	3.614842e+05	1066.8230	172
5.0643725	0.0964543	395.6836	115.1936340	41.1953812	510.8772	3.9662921	3020132	3.18e-05	0.0000223	0.1049477	0.0816655	8.508197	2.602661e+05	982.4358	173
5.6293221	0.0503762	296.6656	7.2726135	-1.7246743	303.9382	5.2787004	24786328	-2.60e-06	-0.0000172	0.0061009	0.0385829	12.408481	6.626173e+06	1184.7838	174
0.6972474	0.0500697	389.9069	120.1149292	16.8878250	510.0218	2.6242695	3262168	1.48e-05	0.0000234	0.1081761	0.0566149	9.031928	3.339698e+05	990.2496	175
3.0376244	0.0248723	336.6650	128.4631042	-6.7371621	465.1281	0.4846936	9463982	-1.40e-06	0.0000023	0.1111375	0.0087165	12.760646	1.667076e+06	1027.4305	176
2.0106387	0.0329384	281.3593	46.0406494	16.6670723	327.4000	1.0789375	11700807	9.80e-06	0.0000082	0.0381245	0.0255465	10.191358	4.534188e+05	1161.5991	177
5.8546815	0.0641725	289.3719	235.8067017	-0.8949147	525.1786	2.8183181	2012873	-5.70e-06	0.0000074	0.1963032	0.0641725	8.677927	2.500000e+05	965.4303	178
3.4739854	0.0936796	332.0297	84.3682556	-20.5946922	416.3979	5.1705766	3285164	-2.58e-05	-0.0000066	0.0725467	0.0813064	9.189792	5.075391e+05	1078.5828	179
4.3419905	0.1004403	294.0995	57.3405151	18.0441093	351.4400	4.3105941	20844514	3.15e-05	0.0000466	0.0475755	0.0673675	9.533856	6.033519e+05	1147.9126	180
0.5569316	0.0252581	293.2041	55.6250610	6.7855310	348.8292	1.0047277	6652898	7.10e-06	0.0000025	0.0465300	0.0191583	11.236161	1.000605e+06	1139.8406	181
4.6844597	0.0899191	265.2206	129.1897278	6.1348190	394.4103	4.9080334	4823923	1.86e-05	-0.0000018	0.1053834	0.0836945	8.936565	3.277787e+05	1096.7118	182
2.0003352	0.0137973	258.7853	46.5346680	-12.1142673	305.3200	0.3783501	18250848	-1.07e-05	-0.0000209	0.0376842	0.0136382	10.238905	2.528532e+05	1188.3231	183
4.1281867	0.0206244	244.4776	3.5987549	16.4063950	248.0764	1.2283775	20625088	1.44e-05	-0.0000040	0.0028761	0.0264954	10.360671	5.800816e+05	1247.6724	184
5.8823633	0.0281540	246.5192	0.0000000	-18.2291698	246.5192	0.9358853	20677346	-2.90e-06	-0.0000041	0.0000000	0.0038040	17.929567	1.525796e+08	1247.5835	185
5.4248939	0.0136898	254.9014	41.0443115	-21.7788792	295.9457	0.4446016	9857540	-1.19e-05	0.0000182	0.0330388	0.0143488	10.426425	3.414408e+05	1201.2626	186
3.3150311	0.0402016	244.2287	103.2835846	4.0881491	347.5123	0.9570600	4787275	7.30e-06	-0.0000197	0.0827222	0.0234196	10.261115	3.877225e+05	1145.2756	187
3.4597988	0.0417778	266.8189	71.1764832	-8.6903248	337.9954	2.0309169	6317984	-3.00e-06	-0.0000144	0.0578107	0.0447895	9.386662	3.375195e+05	1160.0223	188
0.0516132	0.0173071	458.4233	41.6723328	14.8396063	500.0957	0.5835287	8757199	1.10e-05	-0.0000025	0.0370436	0.0174805	10.267204	2.641526e+05	1083.2808	189
4.9503922	0.0282415	396.3953	1.2047424	0.6727794	397.6000	1.6238108	28748238	-1.00e-06	-0.0000097	0.0010307	0.0173038	12.868989	4.055069e+06	1167.6233	190
2.6731162	0.0173940	490.3484	40.1551514	-8.0692482	530.5035	0.5341386	14837678	-5.10e-06	-0.0000047	0.0363962	0.0176008	10.016968	2.649662e+05	1063.1237	191
0.5334142	0.0231568	498.7746	32.7496948	3.6790993	531.5243	0.7983217	30778834	1.70e-06	0.0000039	0.0298475	0.0241809	9.671789	2.627302e+05	1064.4860	192
0.2791347	0.0164883	436.6657	46.9634705	-6.9347901	483.6292	0.4607596	9594235	-6.20e-06	-0.0000104	0.0412852	0.0101999	11.966116	9.924457e+05	1090.5731	193
5.6432557	0.0201462	460.0402	4.8667297	5.3702698	464.9069	0.0409994	27908340	8.50e-06	-0.0000064	0.0043385	0.0018321	13.214683	7.028931e+05	1116.8783	194
4.7224760	0.0312237	425.1406	4.5353699	-6.1203647	429.6760	3.8946359	28792468	-8.60e-06	-0.0000132	0.0039713	0.0368013	11.849185	2.602563e+06	1137.4982	195
3.4412265	0.0581236	413.1343	77.3626404	-13.1347647	490.4969	2.8928027	12410514	-5.10e-06	-0.0000144	0.0682302	0.0577270	9.192840	3.551036e+05	1056.4841	196
5.7488012	0.0019824	403.0272	4.4807434	-11.6008635	407.5079	0.0302962	26576068	-7.00e-07	-0.0000010	0.0038817	0.0019824	12.123337	2.500000e+05	1149.8374	197
5.7297049	0.0043001	450.4884	16.5894775	13.5866623	467.0779	0.0953357	24499846	4.10e-06	0.0000022	0.0146883	0.0038883	12.126189	4.940773e+05	1112.8444	198
6.2770338	0.0073883	409.0697	6.6188049	7.4710326	415.6885	0.1156888	26904310	4.20e-06	-0.0000061	0.0057482	0.0047477	11.761203	3.061542e+05	1144.8336	199
1.8748155	0.0337128	491.6411	12.9872437	-0.2408368	504.6284	1.0283564	27594166	3.50e-06	-0.0000041	0.0117792	0.0183704	11.472737	1.106135e+06	1089.5685	200
5.8307519	0.0040018	546.9897	16.5331421	11.1830244	563.5228	0.1127840	37658472	2.20e-06	0.0000005	0.0155113	0.0038921	12.542143	7.279909e+05	1049.3418	201
4.8755593	0.0125242	610.4081	16.1129761	1.5141308	626.5211	5.6955333	29301964	-7.00e-07	-0.0000104	0.0156765	0.0228373	14.827044	3.608910e+07	1011.7305	202
1.1198250	0.0516330	766.0284	73.1635132	-4.4654346	839.1920	2.2473648	4625500	-8.50e-06	-0.0000046	0.0739607	0.0535007	8.884701	2.582755e+05	916.0586	203
2.8753052	0.0292215	651.4460	30.0136719	1.6261115	681.4597	1.9573126	25579410	2.40e-06	-0.0000082	0.0299933	0.0285597	11.206527	1.291040e+06	970.6663	204
2.2981203	0.0321318	748.5365	62.8635254	-29.6638069	811.4000	0.8673154	5718247	-3.00e-07	-0.0000124	0.0639251	0.0264034	9.505325	2.503942e+05	920.5293	205
5.4933667	0.0082074	557.5358	0.6163940	-5.8482823	558.1522	0.7845772	39970788	-1.80e-06	-0.0000016	0.0005827	0.0054384	15.969561	3.314747e+07	1057.2775	206
2.2368121	0.0190516	688.5466	76.7865601	-13.6870594	765.3331	1.2255987	6659916	-5.90e-06	0.0000012	0.0774845	0.0153964	12.661965	3.410815e+06	914.2054	207
1.4955611	0.0215867	749.5330	27.7869873	-1.8632263	777.3200	0.6936781	10321102	5.80e-06	0.0000164	0.0287053	0.0159639	11.085210	5.581110e+05	940.2227	208
5.2631025	0.0341159	423.8182	25.8463745	16.2617893	449.6646	0.7347412	12123442	2.61e-05	0.0000109	0.0232471	0.0121977	12.372561	1.981743e+06	1085.9628	209
5.1061292	0.0375391	420.8898	2.7238159	0.0076863	423.6136	0.4340235	15188966	-6.20e-06	0.0000018	0.0024182	0.0100736	11.869621	9.405891e+05	1123.6627	210
1.6519448	0.0242032	409.2860	4.2119141	1.1688255	413.4979	0.1198444	25340140	-9.00e-07	-0.0000028	0.0036579	0.0020360	15.447766	5.612531e+06	1147.2560	211
2.6097248	0.0523872	500.2757	86.9176025	14.5672150	587.1933	2.1563251	5296021	3.00e-06	0.0000126	0.0789513	0.0523821	8.871098	2.556373e+05	1013.9838	212
0.9503633	0.0283960	574.0344	17.2127075	1.1010020	591.2471	3.3061759	28036872	-4.50e-06	-0.0000069	0.0163926	0.0350534	11.646445	2.795396e+06	1032.8137	213
4.3836589	0.0438987	446.9225	6.7676392	-15.1327581	453.6901	3.0271022	13639543	-1.69e-05	-0.0000063	0.0059519	0.0368866	11.209352	1.728854e+06	1130.2882	214
4.9086733	0.0131381	636.5668	34.2377319	-2.5751851	670.8045	0.3569922	28434876	2.00e-07	0.0000001	0.0337981	0.0123562	10.513165	2.673191e+05	978.7706	215
1.0992032	0.0357522	656.1976	6.8223267	-4.7597685	663.0199	4.4814758	15595127	3.40e-06	0.0000020	0.0068152	0.0375214	12.117715	4.622542e+06	994.2273	216
0.8734114	0.0467961	763.9097	57.4381104	-1.3624436	820.2256	2.3924296	6445678	5.80e-06	-0.0000198	0.0573713	0.0240031	12.447021	4.303531e+06	943.7264	217
1.2004397	0.0335407	467.9122	0.0000000	-30.5919895	467.4045	0.8767677	18301404	-1.76e-05	-0.0000074	0.0000000	0.0041820	17.363787	9.404506e+07	1126.3516	218
3.2211773	0.0156566	758.4569	89.0411987	-7.9469681	847.4981	0.5709110	10884959	-3.00e-07	0.0000124	0.0902482	0.0161547	10.547763	3.311998e+05	897.5845	219
0.4289240	0.0383854	424.3014	44.5567627	12.8497324	468.8582	2.4017200	8730902	3.30e-06	-0.0000013	0.0401377	0.0479693	9.514470	4.311861e+05	1065.5415	220
0.6257334	0.0223225	422.7468	0.5704651	-6.1506858	420.9482	0.0356192	22535390	-1.80e-06	-0.0000034	0.0005080	0.0007802	16.243025	7.910094e+06	1122.4886	221
0.8723736	0.0406082	415.9424	11.4111328	-15.9500761	427.3535	0.8257976	21666606	-1.90e-05	-0.0000241	0.0099830	0.0105783	13.270001	4.318775e+06	1131.6472	222
2.9487243	0.0150142	435.4266	68.6704712	-17.6360588	504.0971	0.4481404	8334152	-1.73e-05	0.0000031	0.0601024	0.0150142	10.253616	2.500000e+05	1073.8881	223
2.9635582	0.0369213	455.3170	37.2224121	-1.6382087	492.5394	2.5918281	10161785	-1.20e-06	0.0000009	0.0328675	0.0422071	10.248691	6.926337e+05	1095.2772	224
5.5604739	0.0463326	478.6633	10.1799622	-15.2926350	488.8432	1.6424534	11260165	-4.10e-06	-0.0000279	0.0091082	0.0344528	9.970886	5.203885e+05	1107.4857	225
0.7486501	0.0398738	520.0967	27.3995361	-29.8003597	547.4962	5.8864684	6898337	-1.06e-05	-0.0000374	0.0251865	0.0662385	10.384605	1.516752e+06	1060.4672	226
2.0561435	0.0722817	638.9504	88.9130859	-31.0996590	727.8635	4.9599171	17630590	-1.41e-05	0.0000021	0.0876907	0.0891513	8.694615	3.605274e+05	925.0272	227
0.4982191	0.0046339	420.5725	0.3473816	-21.7738514	420.9199	0.1974285	31154784	-9.00e-07	0.0000010	0.0003025	0.0062765	11.911037	6.337882e+05	1147.8453	228
0.9415824	0.0333418	685.0053	120.9772949	9.8603373	805.9826	4.9368238	12197846	2.13e-05	-0.0000004	0.1222536	0.0579639	10.511680	1.489276e+06	868.5833	229
0.4533150	0.0276808	569.9510	16.2549438	-5.2587948	586.2059	0.8741311	35080336	-3.90e-06	-0.0000057	0.0154354	0.0235839	10.004817	3.489940e+05	1036.8392	230
1.6621644	0.0465645	709.5063	113.9003906	3.0864782	823.4067	1.7535986	8739602	9.00e-07	-0.0000105	0.1153194	0.0268082	11.200732	1.066634e+06	873.7946	231
1.0793331	0.0188421	549.2617	2.9382935	-15.5331745	552.2000	0.8555444	43346888	-9.60e-06	0.0000016	0.0027510	0.0072713	14.951643	1.532797e+07	1065.1304	232
4.7141528	0.0411728	673.2536	140.8598633	-3.5922279	814.1135	5.6597176	5829434	-9.30e-06	0.0000176	0.1407636	0.0565879	10.955388	1.624750e+06	859.8240	233
0.1997989	0.0898380	737.5828	175.1514282	-25.5603523	912.7343	22.2926464	11155406	-7.44e-05	-0.0000423	0.1642551	0.1022138	11.871192	8.648636e+06	891.1866	234
4.3153896	0.0417294	606.6025	41.4783936	6.4536052	648.0809	1.7669098	22511344	6.90e-06	0.0000011	0.0401223	0.0364107	9.918600	4.840428e+05	992.3199	235
3.7587399	0.0455800	669.5507	192.0830688	-2.8455110	861.6338	2.0932271	4278624	-1.00e-06	0.0000357	0.1904532	0.0495460	9.033363	2.896092e+05	816.4747	236
0.4886341	0.0264695	446.7435	53.8656616	4.1155024	500.6091	0.9714143	23644486	5.10e-06	-0.0000048	0.0475769	0.0271499	9.670229	2.908541e+05	1078.3162	237
3.2634721	0.0462333	651.3786	180.4982910	-1.3593230	831.8769	2.1312032	4016170	-4.50e-06	-0.0000130	0.1776316	0.0472110	9.283430	2.831561e+05	835.6400	238
4.8792381	0.0350100	698.8641	46.2558594	5.4287648	745.1200	2.3713317	9369743	1.24e-05	0.0000000	0.0450192	0.0266131	11.986229	2.461753e+06	981.2136	239
4.8604302	0.0350307	536.4550	24.5302124	1.6756880	556.9904	3.4628356	7050015	8.40e-06	0.0000078	0.0227755	0.0225153	13.643413	1.077392e+07	1052.5153	240
3.0886207	0.0625205	426.4877	107.0616455	8.6604605	533.5493	3.9157984	9709757	1.29e-05	-0.0000003	0.0946007	0.0702416	9.116240	3.366974e+05	1024.6598	241
4.7125449	0.0430709	543.9439	17.4328613	1.9718819	558.5153	2.1838117	15154295	-5.20e-06	-0.0000100	0.0162658	0.0125189	14.985366	2.016914e+07	1054.3191	242
4.9776549	0.0071861	594.2603	14.2429199	-5.2092190	608.5032	0.0980270	33331160	-6.20e-06	-0.0000185	0.0136767	0.0045645	11.514355	2.804756e+05	1027.1606	243
1.3134558	0.0361383	616.8017	70.7489624	10.9230080	687.5507	0.7635341	27873736	1.42e-05	0.0000016	0.0685625	0.0162369	11.253005	7.648061e+05	961.1417	244
1.5365144	0.0200743	707.2422	28.3293457	-39.8058052	735.5716	0.6016766	8492782	-2.10e-06	-0.0000209	0.0264591	0.0179586	10.229105	2.571117e+05	1042.3542	245
1.5653869	0.0642637	693.9487	150.3799438	4.4182734	844.3287	4.1317534	4198277	1.13e-05	-0.0000052	0.1425063	0.0702679	9.248857	3.677109e+05	904.8715	246
5.2589579	0.0244886	495.4160	20.7609863	-33.9446640	516.1770	7.8633752	12413167	-4.10e-06	0.0000054	0.0188304	0.0387069	13.391194	1.740628e+07	1081.7611	247
2.5970473	0.0414126	654.7842	64.6488647	1.1996871	719.4330	1.5871735	5455541	-5.40e-06	-0.0000073	0.0635211	0.0346254	9.864050	4.573364e+05	953.1056	248
1.4007214	0.0090239	448.6846	0.0270996	-2.9303145	448.7117	0.7591063	31004854	-1.20e-06	0.0000000	0.0000240	0.0047632	16.450443	3.834682e+07	1127.2103	249
0.8245210	0.0256728	476.0648	18.1959229	0.0417716	494.2607	0.9315526	13071406	-1.00e-05	-0.0000239	0.0163459	0.0189812	11.086580	8.757262e+05	1094.9844	250
0.2387125	0.0373257	640.0788	135.1066284	9.9732132	775.1854	1.4466448	4051191	2.68e-05	0.0000057	0.1316210	0.0373257	9.313171	2.500000e+05	891.3756	251
3.1650167	0.0582976	684.5450	215.0700684	-19.3383141	899.6151	3.7874835	4600123	-2.19e-05	0.0000176	0.2033956	0.0594618	9.740731	5.176125e+05	842.3275	252
0.1539084	0.0217563	553.8563	7.6044312	-7.4368749	561.4608	0.2105589	24969802	-2.20e-06	-0.0000008	0.0071202	0.0038663	14.131175	3.025714e+06	1060.4108	253
6.2590008	0.0407595	614.0309	25.7367554	-18.7392502	639.7676	1.8140815	15211728	-2.40e-05	-0.0000230	0.0248460	0.0158342	13.523189	6.052029e+06	1010.1155	254
2.0339530	0.0081791	460.8545	4.9682007	15.8772182	465.8227	0.0549816	21630954	-5.60e-06	-0.0000325	0.0044452	0.0027241	12.262402	3.865949e+05	1112.6747	255
5.7070222	0.0487912	637.5989	100.5481567	-2.9298675	738.1470	1.0949419	4319914	-2.62e-05	-0.0000289	0.0927706	0.0285844	9.750674	3.394846e+05	983.2888	256
1.4888390	0.0328744	580.4739	213.2171021	9.9960365	793.6910	1.4115435	5959184	3.90e-06	0.0000312	0.1893965	0.0345427	9.581040	2.699700e+05	912.5541	257
2.9261842	0.0456825	504.0232	0.0000000	-18.9457607	504.0232	3.6356049	18267958	-2.14e-05	-0.0000010	0.0000000	0.0164082	15.109840	3.564255e+07	1184.2876	258
3.1759019	0.0350546	537.5117	217.7036133	-3.3054497	755.2153	1.2761049	5316965	2.79e-05	0.0000475	0.2174722	0.0282978	10.176279	3.843625e+05	783.3606	259
0.1313870	0.0552536	501.1670	57.0498352	-5.3928800	558.2169	14.3132887	5396011	-1.34e-05	0.0000003	0.0548282	0.0521164	13.624661	2.418783e+07	983.4706	260
2.7316747	0.1100575	447.5252	180.8121033	11.3664055	628.3373	12.5935574	933067	4.47e-05	0.0000407	0.1645338	0.1126272	10.031657	1.691613e+06	918.1237	261
2.1872492	0.1020515	598.5491	186.9020386	-2.6143923	785.4511	5.9461098	9659126	-1.16e-05	-0.0000257	0.2007699	0.0745255	9.909063	1.046691e+06	744.0247	262
3.3050435	0.0772666	389.1741	42.8068542	-17.9812279	431.9810	3.5073516	6945239	-2.10e-05	-0.0000119	0.0370180	0.0709730	8.796840	2.701978e+05	1113.5720	263
3.7477736	0.0490266	312.3352	49.7763672	-3.2326088	362.1115	5.9523091	7687697	-5.40e-06	-0.0000006	0.0397024	0.0610658	10.757839	1.999318e+06	1203.9609	264
5.8797941	0.0233590	419.2320	126.2909241	-5.9968710	545.5229	0.7611175	2356464	-1.04e-05	-0.0000022	0.1124241	0.0233590	9.699438	2.500000e+05	997.0525	265
0.5997809	0.0020718	312.6137	2.8774719	-3.2154458	315.4912	0.0583428	31157796	-1.50e-06	-0.0000011	0.0022743	0.0021020	13.512612	1.079055e+06	1262.3547	266
1.1262323	0.1512970	513.6517	93.0666504	-13.6733532	606.7184	11.1257133	2019027	-3.43e-05	-0.0000398	0.0913323	0.1697887	7.980660	3.340584e+05	925.9227	267
2.4886601	0.0419534	330.9816	36.6696167	7.0230460	367.6512	4.2612963	10583081	1.42e-05	-0.0000038	0.0297245	0.0528026	10.540045	1.400252e+06	1196.9796	268
2.9577181	0.0370231	329.9020	0.0000000	11.2650299	329.9020	11.0781183	44833416	5.70e-06	-0.0000261	0.0000000	0.0138967	18.601353	9.706610e+08	1248.8938	269
2.8375022	0.0328551	327.9011	5.6450806	1.7880259	333.5461	0.6686443	21435456	2.50e-06	-0.0000084	0.0044640	0.0109030	12.613762	2.054760e+06	1258.9403	270
1.0583159	0.0398316	373.8363	12.3865967	-36.1496735	378.1957	2.3610826	9008223	-6.50e-06	0.0000047	0.0105941	0.0168604	13.915173	1.268518e+07	1156.8108	271
0.9322934	0.0468024	343.7322	9.6106873	-1.6000985	353.3429	3.0464876	23240840	-5.00e-07	-0.0000244	0.0077776	0.0363120	11.292027	2.037458e+06	1226.0795	272
3.9269907	0.0350872	311.6323	0.0000916	-5.0474648	305.5102	0.0003087	35399760	-2.60e-06	-0.0000144	0.0000001	0.0000020	23.787418	1.514361e+10	1236.5005	273
5.7543297	0.0345887	345.2795	9.2392883	-25.9180870	354.5188	0.9558230	45151364	-6.50e-06	0.0000154	0.0075853	0.0236483	10.216588	4.425093e+05	1208.8158	274
0.6577826	0.0165050	352.4015	175.2011108	11.0855179	527.6026	0.4718578	4248784	1.91e-05	0.0000165	0.1421374	0.0165050	9.980199	2.500000e+05	1057.4170	275
2.5984280	0.1435135	381.3238	113.9295959	-1.6743664	495.2534	7.7457023	4698559	-4.19e-05	-0.0000782	0.0987924	0.1041623	9.148223	6.744113e+05	1039.2928	276
NA	0.0090858	311.6484	-0.0010071	20.7552757	305.8179	0.0000000	35409160	8.50e-06	-0.0000017	-0.0000008	0.0000000	24.439259	2.839498e+10	1233.2993	277
5.1151824	0.0597004	357.0955	86.8029480	2.3721509	443.8985	4.2398844	5890493	1.48e-05	-0.0000108	0.0708415	0.0731015	9.145505	4.502294e+05	1138.5099	278
1.5893495	0.0026678	311.6640	-0.0012817	19.2195740	309.2482	0.1238878	33475094	2.00e-06	-0.0000006	-0.0000010	0.0001842	24.049967	1.417907e+10	1230.1896	279
4.3595552	0.0272923	325.9868	9.8588257	-7.0028462	335.7126	0.3650776	47015672	-7.20e-06	-0.0000125	0.0080422	0.0054195	14.071792	4.495572e+06	1216.0282	280
5.9763045	0.0210300	360.6233	77.1538391	7.2080526	437.7771	1.5876853	7066355	5.40e-06	0.0000095	0.0637890	0.0243997	11.352305	1.304532e+06	1132.3621	281
0.8573329	0.0043097	314.5449	8.7010803	-0.3690335	323.2460	0.0832206	37591464	-6.00e-07	-0.0000014	0.0071113	0.0043310	11.328163	2.539254e+05	1214.8505	282
1.6665773	0.0323149	408.5267	81.0752258	3.4861221	489.6020	1.2909887	5994325	-1.20e-06	-0.0000035	0.0686844	0.0329227	9.561972	2.553754e+05	1099.3270	283
1.0653715	0.0022306	321.2417	9.0560303	9.4712629	330.2977	0.0809368	41410204	7.00e-07	0.0000001	0.0074140	0.0030230	12.786618	7.342363e+05	1212.4133	284
1.5498948	0.0132680	315.1377	0.0443115	9.1459770	315.1820	0.5591200	38846460	4.20e-06	-0.0000095	0.0000363	0.0131656	11.365131	5.796166e+05	1219.0415	285
0.9260806	0.0294345	385.8515	87.9745789	16.0797234	473.8261	1.0015531	12336873	9.70e-06	0.0000080	0.0746546	0.0294345	9.403262	2.500000e+05	1090.4470	286
0.0370776	0.0258187	408.8651	0.5224304	-54.9695892	409.3875	2.6188071	26198200	-2.25e-05	-0.0000047	0.0004439	0.0093561	16.677107	8.479550e+07	1176.3275	287
4.7861953	0.1046766	367.9671	62.6986389	-0.6526507	430.6658	7.0564413	6014227	7.00e-06	0.0000027	0.0541598	0.0712607	10.527358	1.426677e+06	1094.9611	288
0.4292768	0.0285943	365.2354	46.4255981	-9.4552326	411.6610	3.3958178	36735084	-2.05e-05	0.0000000	0.0389758	0.0374228	11.435380	2.296172e+06	1144.7129	289
0.7673782	0.0320601	405.1554	48.9249268	3.0057271	454.0804	1.3730307	10813427	4.10e-06	-0.0000081	0.0416314	0.0315554	9.896261	4.431221e+05	1126.2682	290
2.8939526	0.0106069	333.0764	8.2748413	-4.7707438	341.3512	0.0388033	36267092	-1.20e-06	-0.0000016	0.0068130	0.0019300	12.855759	4.489074e+05	1206.2843	291
2.0894802	0.0388121	374.6116	53.2908630	3.4914296	427.9024	1.4638753	37895364	5.00e-06	-0.0000156	0.0450224	0.0370263	9.376992	2.984946e+05	1130.3613	292
3.8672376	0.0273431	329.6507	26.0717773	-19.5631504	355.7225	1.7790630	23850142	-6.80e-06	0.0000168	0.0214615	0.0294501	10.837351	1.058091e+06	1188.7429	293
2.0235631	0.0262687	404.1501	83.6476135	8.7460480	487.7977	0.9892258	5407213	3.40e-06	-0.0000024	0.0714268	0.0264329	9.829383	3.281677e+05	1087.4485	294
5.1068153	0.0041811	318.7838	3.7287903	-3.1722131	322.5126	0.0687538	41424516	-2.00e-07	-0.0000010	0.0030766	0.0036706	11.553903	2.500000e+05	1208.2504	295
3.3221855	0.0024824	318.1884	0.8097839	-9.2063131	318.9982	0.0152457	41766672	0.00e+00	0.0000000	0.0006671	0.0011075	12.880982	2.530154e+05	1213.0762	296
1.1817948	0.0265178	350.3470	2.5467834	-1.2200981	352.8937	0.0775257	28072192	-1.70e-06	-0.0000192	0.0021180	0.0024742	13.530375	1.213493e+06	1199.8973	297
6.1442108	0.0225864	338.3526	80.2185364	-6.1125450	418.5712	0.7639796	19692434	-4.90e-06	-0.0000016	0.0666632	0.0222994	9.906095	2.524343e+05	1123.1223	298
1.4508719	0.0395412	377.0728	58.3180237	35.8902550	435.3909	2.9711459	6140913	2.28e-05	0.0000017	0.0493005	0.0455515	10.266127	7.287490e+05	1124.5908	299
2.8245804	0.0294298	434.9175	16.7626343	-15.3329592	451.6801	5.0184517	8011175	-1.53e-05	-0.0000038	0.0144723	0.0403689	12.089800	4.537448e+06	1141.4974	300
5.7868400	0.0391004	373.6391	49.6718445	9.6533585	423.3110	1.9342201	26251616	9.30e-06	0.0000133	0.0423954	0.0422328	9.514456	3.882107e+05	1121.9606	301
1.4800488	0.0343777	380.6569	59.4610596	3.7438395	440.1180	2.0046475	13955069	-4.10e-06	-0.0000086	0.0509214	0.0359285	10.290842	5.752881e+05	1108.2422	302
2.0585692	0.0264874	382.3417	5.2873840	-8.6284599	387.6291	0.5196146	21703370	-4.20e-06	-0.0000079	0.0044894	0.0059489	14.558105	8.450985e+06	1172.4481	303
4.2950392	0.0319944	379.7975	71.7590637	6.0672054	451.5566	1.1560996	24227254	1.11e-05	-0.0000026	0.0616437	0.0315041	9.473256	2.704420e+05	1092.3357	304
1.7185161	0.0302136	424.5866	69.9255676	8.3657217	494.5121	1.3356496	5623618	6.10e-06	0.0000116	0.0606496	0.0321455	9.748511	3.129356e+05	1083.0179	305
3.9584420	0.0297578	403.4018	19.2138062	-10.2802114	422.6156	0.2550426	9798272	-2.10e-06	0.0000089	0.0167796	0.0030633	15.599713	1.289244e+07	1125.8589	306
3.5767550	0.0538813	401.0066	88.8211670	-15.5327320	489.8278	2.2438462	5797071	-1.05e-05	0.0000082	0.0767548	0.0538813	8.850667	2.500000e+05	1068.3859	307
3.3773143	0.0518283	343.7562	78.4984741	0.3278200	422.2547	3.9533837	4166312	7.20e-06	0.0000229	0.0660514	0.0701885	9.143332	3.964212e+05	1109.9462	308
1.9529816	0.0233103	411.9419	93.6799011	23.0334778	505.6218	0.8039615	5271857	1.70e-06	0.0000190	0.0813599	0.0238567	9.746747	2.653577e+05	1057.7465	309
1.0198419	0.0320270	373.8167	18.1132507	12.8461714	391.9299	0.8382998	9417522	4.10e-06	-0.0000084	0.0155490	0.0211133	10.370479	4.633404e+05	1146.7994	310
3.3498497	0.0347173	392.4326	24.7784729	8.6248407	417.2111	2.2132897	9632746	4.40e-06	-0.0000005	0.0214842	0.0391554	10.172898	6.078402e+05	1128.5552	311
4.7124381	0.0268604	414.0009	41.2359314	-45.8984375	449.3200	10.8053122	4580521	-1.05e-05	-0.0000192	0.0360917	0.0409518	13.946174	2.334821e+07	1101.2974	312
1.5580069	0.0774138	373.2288	58.3240356	-2.6912863	431.5528	3.7250218	10407621	-3.20e-06	0.0000003	0.0500315	0.0710138	8.944929	2.762238e+05	1107.4227	313
5.4409776	0.0445295	404.8543	165.0777283	-7.6068397	569.9320	1.8766335	3497986	-2.80e-06	0.0000125	0.1496257	0.0460664	9.069694	2.827442e+05	938.1934	314
3.3008642	0.1118476	350.8338	114.7787476	6.4261456	465.6125	6.3101029	5069278	2.29e-05	0.0000305	0.0992126	0.1034139	8.666373	3.451610e+05	1042.1176	315
1.4400885	0.0068443	401.9545	13.1899719	-0.4702837	415.1445	0.1754715	14974216	-9.00e-07	-0.0000041	0.0117836	0.0070146	11.143725	2.719940e+05	1106.1604	316
2.6907213	0.0162341	442.0473	130.3380737	50.9975243	572.3854	0.4694458	5937160	3.92e-05	0.0000021	0.1200282	0.0136793	10.765446	4.326945e+05	955.5574	317
3.6168063	0.0315609	310.8580	57.1759949	-13.3556309	368.0340	1.0686938	7260449	-5.70e-06	0.0000293	0.0478910	0.0145417	12.562228	2.796399e+06	1136.7015	318
1.7721144	0.0052579	411.1426	63.3243103	-3.8661056	474.4669	0.0188706	5864252	1.00e-06	0.0000164	0.0568570	0.0012435	12.921738	3.000964e+05	1050.4230	319
4.4074984	0.0193537	408.5789	48.4671021	-6.1330667	457.0460	0.6366290	13471446	-1.24e-05	-0.0000047	0.0433716	0.0126815	11.848929	1.112778e+06	1069.0170	320
4.2075830	0.0176322	443.1326	86.1841431	12.8542957	529.3168	0.1132470	4350826	1.61e-05	0.0000271	0.0789177	0.0040810	12.350964	6.410441e+05	1005.8926	321
0.3940128	0.0129430	417.1843	34.8763123	24.1077709	452.0606	0.2793030	11980740	1.58e-05	0.0000076	0.0312920	0.0101134	10.750885	3.085303e+05	1079.6692	322
4.0617018	0.0101842	295.4418	33.5981750	-0.0986546	329.0400	0.2275746	10337817	1.50e-06	-0.0000105	0.0277441	0.0077921	11.347039	4.625406e+05	1177.4059	323
1.2054230	0.0221615	304.1122	61.3278503	33.7436752	365.4400	0.8547414	6270944	2.48e-05	0.0000111	0.0508863	0.0191413	10.835736	6.279414e+05	1143.8651	324
2.9422362	0.0667855	429.3797	138.9470520	67.5464096	568.3267	3.6781981	2867567	3.79e-05	0.0000443	0.1256415	0.0705697	8.939951	3.177699e+05	966.9543	325
1.4431311	0.0275921	395.0110	10.3090515	6.7359524	405.3200	1.4179286	7208180	7.50e-06	-0.0000202	0.0091614	0.0194761	12.027444	1.823387e+06	1114.9631	326
1.8318641	0.0241384	414.5231	30.7442932	16.0769787	445.2674	1.7177018	14602073	1.92e-05	0.0000032	0.0276528	0.0268306	11.145473	1.137585e+06	1081.0516	327
3.6635244	0.0304742	268.5699	40.5500793	2.7787406	309.1200	0.6967074	10768462	1.75e-05	0.0000160	0.0329330	0.0178096	10.631125	5.034690e+05	1190.7389	328
0.3342191	0.0225697	423.8070	11.5087891	8.1416149	435.3158	0.7485024	20814114	6.90e-06	-0.0000010	0.0104184	0.0207088	10.169434	3.439068e+05	1093.1564	329
1.5593089	0.0065369	370.9563	74.8837280	-42.0582619	445.8400	0.0926008	5995710	-1.17e-05	-0.0000008	0.0654060	0.0034802	12.524646	4.840523e+05	1070.0225	330
2.8873963	0.0121421	327.2518	48.5842896	-9.3623762	375.8361	0.4122616	7595339	-4.90e-06	0.0000029	0.0410955	0.0081734	12.629429	1.521933e+06	1133.6447	331
1.4493849	0.0317811	345.8135	19.2011108	-9.5260048	365.0146	7.3057508	5286377	-3.20e-06	-0.0000066	0.0164681	0.0250804	15.051330	4.807216e+07	1146.7543	332
2.3094857	0.0039922	398.1082	46.3889465	-11.7856216	444.4972	0.0237033	19577646	-4.90e-06	-0.0000026	0.0413690	0.0014771	12.760317	3.632004e+05	1074.9572	333
4.9947395	0.0339115	424.6793	30.3173523	38.8196487	454.9966	3.8403807	11229707	1.48e-05	0.0000105	0.0274660	0.0432207	11.130878	1.828470e+06	1073.4968	334
1.3712677	0.0299273	304.1875	86.3724670	10.1909161	390.5600	1.4645642	6853189	3.40e-06	-0.0000058	0.0718553	0.0330846	9.856502	3.720996e+05	1115.6613	335
2.8796389	0.0486747	293.0000	0.0000000	-5.4007034	293.0000	15.3985491	5898558	7.40e-06	0.0000035	0.0000000	0.0315047	15.946224	1.625288e+08	1210.6652	336
6.2598667	0.0281482	293.5083	19.7316589	-10.5386744	313.2400	1.4714558	12470167	-5.90e-06	-0.0000033	0.0163044	0.0308698	10.162689	4.093885e+05	1190.4697	337
0.8019695	0.0430367	390.0986	39.8400879	-28.7221737	429.9387	1.2018919	8904691	-2.80e-06	-0.0000278	0.0355077	0.0250125	10.550935	6.759058e+05	1082.1726	338
4.1139579	0.0278110	422.6501	1.1092529	-16.6053028	422.1921	0.2489745	25430700	-1.90e-05	-0.0000154	0.0010127	0.0019414	17.477892	5.254351e+07	1094.1910	339
2.7905071	0.0091433	316.3464	1.3756409	-18.5778313	317.1600	1.2064096	36850112	-5.60e-06	-0.0000102	0.0011540	0.0073243	15.779985	3.350427e+07	1190.7104	340
0.8121135	0.0143237	389.5445	49.4335022	-31.2736912	434.7037	0.5196902	15207513	-3.00e-06	-0.0000020	0.0440043	0.0076469	13.491356	3.908805e+06	1073.9458	341
1.5712389	0.0303363	442.1251	15.0446167	22.0060616	452.9889	0.9384912	9798053	8.60e-06	-0.0000143	0.0139250	0.0063439	15.762863	2.224587e+07	1065.3601	342
0.7165681	0.0196909	306.8919	78.5992737	23.1158810	385.4912	1.0493386	14079120	1.63e-05	0.0000064	0.0654889	0.0236146	10.455995	5.790263e+05	1121.5930	343
2.9996293	0.0129644	311.0177	57.1249695	-1.1237273	368.1427	1.6552866	8762238	1.60e-06	-0.0000011	0.0477261	0.0105905	15.003557	1.965581e+07	1139.8094	344
3.2424917	0.0150626	304.2055	48.8230591	2.8778739	353.0286	0.5871220	10180995	4.00e-06	0.0000015	0.0405821	0.0088408	13.179790	2.564778e+06	1154.2452	345
0.6110631	0.0200919	354.8913	104.6942139	29.9164944	459.5855	0.6697589	5919100	6.20e-06	0.0000018	0.0903933	0.0200749	10.023645	3.153393e+05	1053.5138	346

5.3 Work with vector data

Two type of data exist in the G.I.S. world raster and vector. if we already saw the first one quickly with covariates of Germany here is an example of importing either ESRI Shapefile or .GPKG files into R and handling them the package sf and sp are useful to manage vector data.

To import a shapefile st_read("./file_location/yourfile.shp") command can be used and add the , layer = "layer_name" parameter when importing a gpkg file.

#Import the area shapefile

# install.packages(sf, mapview)
library(sf)
library(mapview)

Area = st_read("./Data/Example_Irak/Area.shp")
plot(Area)
mapview(Area)

5.3.1 Coordinate systems

You can also check the coordinate system with crs() (Raster) or st_crs() (Vector) set one if their is none proj4string(vector) <- CRS() or transformed it with st_transform(vector, new_crs)) from sf package for vector or using project(raster, new_crs) for raster layers.

library(raster)
st_crs(Area)

## Coordinate Reference System:
##   User input: WGS 84 / UTM zone 38N 
##   wkt:
## PROJCRS["WGS 84 / UTM zone 38N",
##     BASEGEOGCRS["WGS 84",
##         DATUM["World Geodetic System 1984",
##             ELLIPSOID["WGS 84",6378137,298.257223563,
##                 LENGTHUNIT["metre",1]]],
##         PRIMEM["Greenwich",0,
##             ANGLEUNIT["degree",0.0174532925199433]],
##         ID["EPSG",4326]],
##     CONVERSION["UTM zone 38N",
##         METHOD["Transverse Mercator",
##             ID["EPSG",9807]],
##         PARAMETER["Latitude of natural origin",0,
##             ANGLEUNIT["Degree",0.0174532925199433],
##             ID["EPSG",8801]],
##         PARAMETER["Longitude of natural origin",45,
##             ANGLEUNIT["Degree",0.0174532925199433],
##             ID["EPSG",8802]],
##         PARAMETER["Scale factor at natural origin",0.9996,
##             SCALEUNIT["unity",1],
##             ID["EPSG",8805]],
##         PARAMETER["False easting",500000,
##             LENGTHUNIT["metre",1],
##             ID["EPSG",8806]],
##         PARAMETER["False northing",0,
##             LENGTHUNIT["metre",1],
##             ID["EPSG",8807]]],
##     CS[Cartesian,2],
##         AXIS["(E)",east,
##             ORDER[1],
##             LENGTHUNIT["metre",1]],
##         AXIS["(N)",north,
##             ORDER[2],
##             LENGTHUNIT["metre",1]],
##     ID["EPSG",32638]]

# Transformed it into WG84

Area_WGS84 <- st_transform(Area, CRS("+init=epsg:4326"))

## Warning in CPL_crs_from_input(x): GDAL Message 1: +init=epsg:XXXX syntax is
## deprecated. It might return a CRS with a non-EPSG compliant axis order.

st_crs(Area_WGS84)

## Coordinate Reference System:
##   User input: GEOGCRS["WGS 84",
##     ENSEMBLE["World Geodetic System 1984 ensemble",
##         MEMBER["World Geodetic System 1984 (Transit)",
##             ID["EPSG",1166]],
##         MEMBER["World Geodetic System 1984 (G730)",
##             ID["EPSG",1152]],
##         MEMBER["World Geodetic System 1984 (G873)",
##             ID["EPSG",1153]],
##         MEMBER["World Geodetic System 1984 (G1150)",
##             ID["EPSG",1154]],
##         MEMBER["World Geodetic System 1984 (G1674)",
##             ID["EPSG",1155]],
##         MEMBER["World Geodetic System 1984 (G1762)",
##             ID["EPSG",1156]],
##         MEMBER["World Geodetic System 1984 (G2139)",
##             ID["EPSG",1309]],
##         ELLIPSOID["WGS 84",6378137,298.257223563,
##             LENGTHUNIT["metre",1],
##             ID["EPSG",7030]],
##         ENSEMBLEACCURACY[2.0],
##         ID["EPSG",6326]],
##     PRIMEM["Greenwich",0,
##         ANGLEUNIT["degree",0.0174532925199433],
##         ID["EPSG",8901]],
##     CS[ellipsoidal,2],
##         AXIS["longitude",east,
##             ORDER[1],
##             ANGLEUNIT["degree",0.0174532925199433,
##                 ID["EPSG",9122]]],
##         AXIS["latitude",north,
##             ORDER[2],
##             ANGLEUNIT["degree",0.0174532925199433,
##                 ID["EPSG",9122]]],
##     USAGE[
##         SCOPE["unknown"],
##         AREA["World."],
##         BBOX[-90,-180,90,180]]] 
##   wkt:
## GEOGCRS["WGS 84",
##     ENSEMBLE["World Geodetic System 1984 ensemble",
##         MEMBER["World Geodetic System 1984 (Transit)",
##             ID["EPSG",1166]],
##         MEMBER["World Geodetic System 1984 (G730)",
##             ID["EPSG",1152]],
##         MEMBER["World Geodetic System 1984 (G873)",
##             ID["EPSG",1153]],
##         MEMBER["World Geodetic System 1984 (G1150)",
##             ID["EPSG",1154]],
##         MEMBER["World Geodetic System 1984 (G1674)",
##             ID["EPSG",1155]],
##         MEMBER["World Geodetic System 1984 (G1762)",
##             ID["EPSG",1156]],
##         MEMBER["World Geodetic System 1984 (G2139)",
##             ID["EPSG",1309]],
##         ELLIPSOID["WGS 84",6378137,298.257223563,
##             LENGTHUNIT["metre",1],
##             ID["EPSG",7030]],
##         ENSEMBLEACCURACY[2.0],
##         ID["EPSG",6326]],
##     PRIMEM["Greenwich",0,
##         ANGLEUNIT["degree",0.0174532925199433],
##         ID["EPSG",8901]],
##     CS[ellipsoidal,2],
##         AXIS["longitude",east,
##             ORDER[1],
##             ANGLEUNIT["degree",0.0174532925199433,
##                 ID["EPSG",9122]]],
##         AXIS["latitude",north,
##             ORDER[2],
##             ANGLEUNIT["degree",0.0174532925199433,
##                 ID["EPSG",9122]]],
##     USAGE[
##         SCOPE["unknown"],
##         AREA["World."],
##         BBOX[-90,-180,90,180]]]

5.3.2 Resample, crop and mask

One major complication of raster data is when using stack raster the overlaying layers must have the same extent and pixel size. One line of command resample(raster to resample , raster with final parameters , method) simplify this process if your data have already a similar CRS. The method parameter can be “bilinear” when dealing with continuous data or “ngb” (nearest neighbor) mostly for discrete data.
When your want to crop your covariates to a study area extent you can use the crop(raster, extent) command and for masking the raster or stack raster mask(raster, extent)`

library(raster)

# Import and check the CRS of the DEM
DEM <- stack("./Data/Example_Irak/Download/DEM_GLO30.tif")
crs(DEM)

## Coordinate Reference System:
## Deprecated Proj.4 representation:
##  +proj=utm +zone=38 +datum=WGS84 +units=m +no_defs 
## WKT2 2019 representation:
## PROJCRS["unknown",
##     BASEGEOGCRS["unknown",
##         DATUM["World Geodetic System 1984",
##             ELLIPSOID["WGS 84",6378137,298.257223563,
##                 LENGTHUNIT["metre",1]],
##             ID["EPSG",6326]],
##         PRIMEM["Greenwich",0,
##             ANGLEUNIT["degree",0.0174532925199433],
##             ID["EPSG",8901]]],
##     CONVERSION["UTM zone 38N",
##         METHOD["Transverse Mercator",
##             ID["EPSG",9807]],
##         PARAMETER["Latitude of natural origin",0,
##             ANGLEUNIT["degree",0.0174532925199433],
##             ID["EPSG",8801]],
##         PARAMETER["Longitude of natural origin",45,
##             ANGLEUNIT["degree",0.0174532925199433],
##             ID["EPSG",8802]],
##         PARAMETER["Scale factor at natural origin",0.9996,
##             SCALEUNIT["unity",1],
##             ID["EPSG",8805]],
##         PARAMETER["False easting",500000,
##             LENGTHUNIT["metre",1],
##             ID["EPSG",8806]],
##         PARAMETER["False northing",0,
##             LENGTHUNIT["metre",1],
##             ID["EPSG",8807]],
##         ID["EPSG",16038]],
##     CS[Cartesian,2],
##         AXIS["(E)",east,
##             ORDER[1],
##             LENGTHUNIT["metre",1,
##                 ID["EPSG",9001]]],
##         AXIS["(N)",north,
##             ORDER[2],
##             LENGTHUNIT["metre",1,
##                 ID["EPSG",9001]]]]

proj4string(DEM) <- CRS("+init=epsg:32638")
crs(DEM)

## Coordinate Reference System:
## Deprecated Proj.4 representation:
##  +proj=utm +zone=38 +datum=WGS84 +units=m +no_defs 
## WKT2 2019 representation:
## PROJCRS["WGS 84 / UTM zone 38N",
##     BASEGEOGCRS["WGS 84",
##         ENSEMBLE["World Geodetic System 1984 ensemble",
##             MEMBER["World Geodetic System 1984 (Transit)"],
##             MEMBER["World Geodetic System 1984 (G730)"],
##             MEMBER["World Geodetic System 1984 (G873)"],
##             MEMBER["World Geodetic System 1984 (G1150)"],
##             MEMBER["World Geodetic System 1984 (G1674)"],
##             MEMBER["World Geodetic System 1984 (G1762)"],
##             MEMBER["World Geodetic System 1984 (G2139)"],
##             ELLIPSOID["WGS 84",6378137,298.257223563,
##                 LENGTHUNIT["metre",1]],
##             ENSEMBLEACCURACY[2.0]],
##         PRIMEM["Greenwich",0,
##             ANGLEUNIT["degree",0.0174532925199433]],
##         ID["EPSG",4326]],
##     CONVERSION["UTM zone 38N",
##         METHOD["Transverse Mercator",
##             ID["EPSG",9807]],
##         PARAMETER["Latitude of natural origin",0,
##             ANGLEUNIT["degree",0.0174532925199433],
##             ID["EPSG",8801]],
##         PARAMETER["Longitude of natural origin",45,
##             ANGLEUNIT["degree",0.0174532925199433],
##             ID["EPSG",8802]],
##         PARAMETER["Scale factor at natural origin",0.9996,
##             SCALEUNIT["unity",1],
##             ID["EPSG",8805]],
##         PARAMETER["False easting",500000,
##             LENGTHUNIT["metre",1],
##             ID["EPSG",8806]],
##         PARAMETER["False northing",0,
##             LENGTHUNIT["metre",1],
##             ID["EPSG",8807]],
##         ID["EPSG",16038]],
##     CS[Cartesian,2],
##         AXIS["(E)",east,
##             ORDER[1],
##             LENGTHUNIT["metre",1,
##                 ID["EPSG",9001]]],
##         AXIS["(N)",north,
##             ORDER[2],
##             LENGTHUNIT["metre",1,
##                 ID["EPSG",9001]]],
##     USAGE[
##         SCOPE["unknown"],
##         AREA["Between 42°E and 48°E, northern hemisphere between equator and 84°N, onshore and offshore. Armenia. Azerbaijan. Djibouti. Eritrea. Ethiopia. Georgia. Islamic Republic of Iran. Iraq. kazakhstan. Kuwait. Russian Federation. Saudi Arabia. Somalia. Türkiye (Turkey). Yemen."],
##         BBOX[0,42,84,48]]]

st_crs(Area)

## Coordinate Reference System:
##   User input: WGS 84 / UTM zone 38N 
##   wkt:
## PROJCRS["WGS 84 / UTM zone 38N",
##     BASEGEOGCRS["WGS 84",
##         DATUM["World Geodetic System 1984",
##             ELLIPSOID["WGS 84",6378137,298.257223563,
##                 LENGTHUNIT["metre",1]]],
##         PRIMEM["Greenwich",0,
##             ANGLEUNIT["degree",0.0174532925199433]],
##         ID["EPSG",4326]],
##     CONVERSION["UTM zone 38N",
##         METHOD["Transverse Mercator",
##             ID["EPSG",9807]],
##         PARAMETER["Latitude of natural origin",0,
##             ANGLEUNIT["Degree",0.0174532925199433],
##             ID["EPSG",8801]],
##         PARAMETER["Longitude of natural origin",45,
##             ANGLEUNIT["Degree",0.0174532925199433],
##             ID["EPSG",8802]],
##         PARAMETER["Scale factor at natural origin",0.9996,
##             SCALEUNIT["unity",1],
##             ID["EPSG",8805]],
##         PARAMETER["False easting",500000,
##             LENGTHUNIT["metre",1],
##             ID["EPSG",8806]],
##         PARAMETER["False northing",0,
##             LENGTHUNIT["metre",1],
##             ID["EPSG",8807]]],
##     CS[Cartesian,2],
##         AXIS["(E)",east,
##             ORDER[1],
##             LENGTHUNIT["metre",1]],
##         AXIS["(N)",north,
##             ORDER[2],
##             LENGTHUNIT["metre",1]],
##     ID["EPSG",32638]]

# Crop the DEM
DEM_crop <- crop(DEM,Area)

# Mask the covariates
DEM_mask <- mask(DEM_crop, Area)

Plot the three different rasters: original, cropped and masked on the same map.

mapview(DEM$DEM_GLO30) + mapview(DEM_crop$DEM_GLO30) + mapview(DEM_mask$DEM_GLO30)