forked from r-spatial/stars
-
Notifications
You must be signed in to change notification settings - Fork 0
/
sp.Rout.save
185 lines (161 loc) · 7.23 KB
/
sp.Rout.save
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
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
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
R version 4.0.2 (2020-06-22) -- "Taking Off Again"
Copyright (C) 2020 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> options(rgdal_show_exportToProj4_warnings = "none")
> library(sp)
> demo(meuse)
demo(meuse)
---- ~~~~~
> require(sp)
> crs = CRS("+init=epsg:28992")
> data("meuse")
> coordinates(meuse) <- ~x+y
> proj4string(meuse) <- crs
> data("meuse.grid")
> coordinates(meuse.grid) <- ~x+y
> gridded(meuse.grid) <- TRUE
> proj4string(meuse.grid) <- crs
> data("meuse.riv")
> meuse.riv <- SpatialPolygons(list(Polygons(list(Polygon(meuse.riv)),"meuse.riv")))
> proj4string(meuse.riv) <- crs
> data("meuse.area")
> meuse.area = SpatialPolygons(list(Polygons(list(Polygon(meuse.area)), "area")))
> proj4string(meuse.area) <- crs
> # remove +init=xxx
> meuse@proj4string = CRS(
+ "+proj=sterea +lat_0=52.1561605555556 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +no_defs"
+ )
> meuse.grid@proj4string = CRS(
+ "+proj=sterea +lat_0=52.1561605555556 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +no_defs"
+ )
>
> suppressPackageStartupMessages(library(stars))
> (m = st_as_stars(meuse))
stars object with 1 dimensions and 12 attributes
attribute(s):
cadmium copper lead zinc
Min. : 0.200 Min. : 14.00 Min. : 37.0 Min. : 113.0
1st Qu.: 0.800 1st Qu.: 23.00 1st Qu.: 72.5 1st Qu.: 198.0
Median : 2.100 Median : 31.00 Median :123.0 Median : 326.0
Mean : 3.246 Mean : 40.32 Mean :153.4 Mean : 469.7
3rd Qu.: 3.850 3rd Qu.: 49.50 3rd Qu.:207.0 3rd Qu.: 674.5
Max. :18.100 Max. :128.00 Max. :654.0 Max. :1839.0
elev dist om ffreq soil lime
Min. : 5.180 Min. :0.00000 Min. : 1.000 1:84 1:97 0:111
1st Qu.: 7.546 1st Qu.:0.07569 1st Qu.: 5.300 2:48 2:46 1: 44
Median : 8.180 Median :0.21184 Median : 6.900 3:23 3:12
Mean : 8.165 Mean :0.24002 Mean : 7.478
3rd Qu.: 8.955 3rd Qu.:0.36407 3rd Qu.: 9.000
Max. :10.520 Max. :0.88039 Max. :17.000
NA's :2
landuse dist.m
W :50 Min. : 10.0
Ah :39 1st Qu.: 80.0
Am :22 Median : 270.0
Fw :10 Mean : 290.3
Ab : 8 3rd Qu.: 450.0
(Other):25 Max. :1000.0
NA's : 1
dimension(s):
from to offset delta refsys point
geometry 1 155 NA NA +proj=sterea +lat_0=52.15... TRUE
values
geometry POINT (181072 333611),...,POINT (180627 330190)
> (gr = st_as_stars(meuse.grid))
stars object with 2 dimensions and 5 attributes
attribute(s):
part.a part.b dist soil ffreq
Min. :0.000 Min. :0.000 Min. :0.000 1 :1665 1 : 779
1st Qu.:0.000 1st Qu.:0.000 1st Qu.:0.119 2 :1084 2 :1335
Median :0.000 Median :1.000 Median :0.272 3 : 354 3 : 989
Mean :0.399 Mean :0.601 Mean :0.297 NA's:5009 NA's:5009
3rd Qu.:1.000 3rd Qu.:1.000 3rd Qu.:0.440
Max. :1.000 Max. :1.000 Max. :0.993
NA's :5009 NA's :5009 NA's :5009
dimension(s):
from to offset delta refsys point values x/y
x 1 78 178440 40 +proj=sterea +lat_0=52.15... NA NULL [x]
y 1 104 333760 -40 +proj=sterea +lat_0=52.15... NA NULL [y]
>
> m2 = as(m, "Spatial")
> summary(m2)
Object of class SpatialPointsDataFrame
Coordinates:
min max
coords.x1 178605 181390
coords.x2 329714 333611
Is projected: TRUE
proj4string :
[+proj=sterea +lat_0=52.1561605555556 +lon_0=5.38763888888889
+k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +no_defs]
Number of points: 155
Data attributes:
cadmium copper lead zinc
Min. : 0.200 Min. : 14.00 Min. : 37.0 Min. : 113.0
1st Qu.: 0.800 1st Qu.: 23.00 1st Qu.: 72.5 1st Qu.: 198.0
Median : 2.100 Median : 31.00 Median :123.0 Median : 326.0
Mean : 3.246 Mean : 40.32 Mean :153.4 Mean : 469.7
3rd Qu.: 3.850 3rd Qu.: 49.50 3rd Qu.:207.0 3rd Qu.: 674.5
Max. :18.100 Max. :128.00 Max. :654.0 Max. :1839.0
elev dist om ffreq soil lime
Min. : 5.180 Min. :0.00000 Min. : 1.000 1:84 1:97 0:111
1st Qu.: 7.546 1st Qu.:0.07569 1st Qu.: 5.300 2:48 2:46 1: 44
Median : 8.180 Median :0.21184 Median : 6.900 3:23 3:12
Mean : 8.165 Mean :0.24002 Mean : 7.478
3rd Qu.: 8.955 3rd Qu.:0.36407 3rd Qu.: 9.000
Max. :10.520 Max. :0.88039 Max. :17.000
NA's :2
landuse dist.m
W :50 Min. : 10.0
Ah :39 1st Qu.: 80.0
Am :22 Median : 270.0
Fw :10 Mean : 290.3
Ab : 8 3rd Qu.: 450.0
(Other):25 Max. :1000.0
NA's : 1
> gr2 = as(gr, "Spatial")
> summary(gr2)
Object of class SpatialGridDataFrame
Coordinates:
min max
x 178440 181560
y 329600 333760
Is projected: TRUE
proj4string :
[+proj=sterea +lat_0=52.1561605555556 +lon_0=5.38763888888889
+k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +no_defs]
Grid attributes:
cellcentre.offset cellsize cells.dim
x 178460 40 78
y 329620 40 104
Data attributes:
part.a part.b dist soil ffreq
Min. :0.000 Min. :0.000 Min. :0.000 1 :1665 1 : 779
1st Qu.:0.000 1st Qu.:0.000 1st Qu.:0.119 2 :1084 2 :1335
Median :0.000 Median :1.000 Median :0.272 3 : 354 3 : 989
Mean :0.399 Mean :0.601 Mean :0.297 NA's:5009 NA's:5009
3rd Qu.:1.000 3rd Qu.:1.000 3rd Qu.:0.440
Max. :1.000 Max. :1.000 Max. :0.993
NA's :5009 NA's :5009 NA's :5009
>
> all.equal(meuse, m2, check.attributes = FALSE)
[1] TRUE
> fullgrid(meuse.grid) = TRUE
> all.equal(meuse.grid, gr2, check.attributes = FALSE)
[1] TRUE
>
> #all.equal(meuse, m2)
> #all.equal(meuse.grid, gr2)
>
> proc.time()
user system elapsed
0.970 0.081 1.043