-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathr-raster-vector-geospatial--12-time-series-raster.diff
145 lines (139 loc) · 7.81 KB
/
r-raster-vector-geospatial--12-time-series-raster.diff
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
diff --git a/home/zhian/Documents/Carpentries/Git/zkamvar--postmaul/data/rmd-repos/datacarpentry--r-raster-vector-geospatialR3/_episodes/12-time-series-raster.md b/home/zhian/Documents/Carpentries/Git/zkamvar--postmaul/data/rmd-repos/datacarpentry--r-raster-vector-geospatialR4/_episodes/12-time-series-raster.md
index e10a8a7..a1d7307 100644
--- a/home/zhian/Documents/Carpentries/Git/zkamvar--postmaul/data/rmd-repos/datacarpentry--r-raster-vector-geospatialR3/_episodes/12-time-series-raster.md
+++ b/home/zhian/Documents/Carpentries/Git/zkamvar--postmaul/data/rmd-repos/datacarpentry--r-raster-vector-geospatialR4/_episodes/12-time-series-raster.md
@@ -133,6 +133,40 @@ NDVI_HARV_stack <- stack(all_NDVI_HARV)
~~~
{: .language-r}
+
+
+~~~
+Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
+datum Unknown based on WGS84 ellipsoid in CRS definition
+Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
+datum Unknown based on WGS84 ellipsoid in CRS definition
+Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
+datum Unknown based on WGS84 ellipsoid in CRS definition
+Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
+datum Unknown based on WGS84 ellipsoid in CRS definition
+Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
+datum Unknown based on WGS84 ellipsoid in CRS definition
+Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
+datum Unknown based on WGS84 ellipsoid in CRS definition
+Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
+datum Unknown based on WGS84 ellipsoid in CRS definition
+Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
+datum Unknown based on WGS84 ellipsoid in CRS definition
+Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
+datum Unknown based on WGS84 ellipsoid in CRS definition
+Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
+datum Unknown based on WGS84 ellipsoid in CRS definition
+Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
+datum Unknown based on WGS84 ellipsoid in CRS definition
+Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
+datum Unknown based on WGS84 ellipsoid in CRS definition
+Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
+datum Unknown based on WGS84 ellipsoid in CRS definition
+Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
+datum Unknown based on WGS84 ellipsoid in CRS definition
+~~~
+{: .error}
+
We can explore the GeoTIFF tags (the embedded metadata) in a stack using the
same syntax that we used on single-band raster objects in R including: `crs()`
(coordinate reference system), `extent()` and `res()` (resolution; specifically
@@ -270,6 +304,14 @@ NDVI_HARV_stack <- NDVI_HARV_stack/10000
~~~
{: .language-r}
+
+
+~~~
+Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
+datum Unknown based on WGS84 ellipsoid in CRS definition
+~~~
+{: .error}
+
After applying our scale factor, we can recreate our plot using the same code we used above.
@@ -355,52 +397,52 @@ str(har_met_daily)
~~~
'data.frame': 5345 obs. of 46 variables:
- $ date : Factor w/ 5345 levels "2001-02-11","2001-02-12",..: 1 2 3 4 5 6 7 8 9 10 ...
+ $ date : chr "2001-02-11" "2001-02-12" "2001-02-13" "2001-02-14" ...
$ jd : int 42 43 44 45 46 47 48 49 50 51 ...
$ airt : num -10.7 -9.8 -2 -0.5 -0.4 -3 -4.5 -9.9 -4.5 3.2 ...
- $ f.airt : Factor w/ 2 levels "","E": 1 1 1 1 1 1 1 1 1 1 ...
+ $ f.airt : chr "" "" "" "" ...
$ airtmax : num -6.9 -2.4 5.7 1.9 2.4 1.3 -0.7 -3.3 0.7 8.9 ...
- $ f.airtmax: Factor w/ 2 levels "","E": 1 1 1 1 1 1 1 1 1 1 ...
+ $ f.airtmax: chr "" "" "" "" ...
$ airtmin : num -15.1 -17.4 -7.3 -5.7 -5.7 -9 -12.7 -17.1 -11.7 -1.3 ...
- $ f.airtmin: Factor w/ 2 levels "","E": 1 1 1 1 1 1 1 1 1 1 ...
+ $ f.airtmin: chr "" "" "" "" ...
$ rh : int 40 45 70 78 69 82 66 51 57 62 ...
- $ f.rh : Factor w/ 3 levels "","E","M": 1 1 1 1 1 1 1 1 1 1 ...
+ $ f.rh : chr "" "" "" "" ...
$ rhmax : int 58 85 100 100 100 100 100 71 81 78 ...
- $ f.rhmax : Factor w/ 3 levels "","E","M": 1 1 1 1 1 1 1 1 1 1 ...
+ $ f.rhmax : chr "" "" "" "" ...
$ rhmin : int 22 14 34 59 37 46 30 34 37 42 ...
- $ f.rhmin : Factor w/ 3 levels "","E","M": 1 1 1 1 1 1 1 1 1 1 ...
+ $ f.rhmin : chr "" "" "" "" ...
$ dewp : num -22.2 -20.7 -7.6 -4.1 -6 -5.9 -10.8 -18.5 -12 -3.5 ...
- $ f.dewp : Factor w/ 3 levels "","E","M": 1 1 1 1 1 1 1 1 1 1 ...
+ $ f.dewp : chr "" "" "" "" ...
$ dewpmax : num -16.8 -9.2 -4.6 1.9 2 -0.4 -0.7 -14.4 -4 0.6 ...
- $ f.dewpmax: Factor w/ 3 levels "","E","M": 1 1 1 1 1 1 1 1 1 1 ...
+ $ f.dewpmax: chr "" "" "" "" ...
$ dewpmin : num -25.7 -27.9 -10.2 -10.2 -12.1 -10.6 -25.4 -25 -16.5 -5.7 ...
- $ f.dewpmin: Factor w/ 3 levels "","E","M": 1 1 1 1 1 1 1 1 1 1 ...
+ $ f.dewpmin: chr "" "" "" "" ...
$ prec : num 0 0 0 6.9 0 2.3 0 0 0 0 ...
- $ f.prec : Factor w/ 2 levels "","E": 1 1 1 1 1 1 1 1 1 1 ...
+ $ f.prec : chr "" "" "" "" ...
$ slrt : num 14.9 14.8 14.8 2.6 10.5 6.4 10.3 15.5 15 7.7 ...
- $ f.slrt : Factor w/ 3 levels "","E","M": 1 1 1 1 1 1 1 1 1 1 ...
+ $ f.slrt : chr "" "" "" "" ...
$ part : num NA NA NA NA NA NA NA NA NA NA ...
- $ f.part : Factor w/ 3 levels "","E","M": 3 3 3 3 3 3 3 3 3 3 ...
+ $ f.part : chr "M" "M" "M" "M" ...
$ netr : num NA NA NA NA NA NA NA NA NA NA ...
- $ f.netr : Factor w/ 3 levels "","E","M": 3 3 3 3 3 3 3 3 3 3 ...
+ $ f.netr : chr "M" "M" "M" "M" ...
$ bar : int 1025 1033 1024 1016 1010 1016 1008 1022 1022 1017 ...
- $ f.bar : Factor w/ 3 levels "","E","M": 1 1 1 1 1 1 1 1 1 1 ...
+ $ f.bar : chr "" "" "" "" ...
$ wspd : num 3.3 1.7 1.7 2.5 1.6 1.1 3.3 2 2.5 2 ...
- $ f.wspd : Factor w/ 3 levels "","M","Q": 1 1 1 1 1 1 1 1 1 1 ...
+ $ f.wspd : chr "" "" "" "" ...
$ wres : num 2.9 0.9 0.9 1.9 1.2 0.5 3 1.9 2.1 1.8 ...
- $ f.wres : Factor w/ 3 levels "","M","Q": 1 1 1 1 1 1 1 1 1 1 ...
+ $ f.wres : chr "" "" "" "" ...
$ wdir : int 287 245 278 197 300 182 281 272 217 218 ...
- $ f.wdir : Factor w/ 3 levels "","M","Q": 1 1 1 1 1 1 1 1 1 1 ...
+ $ f.wdir : chr "" "" "" "" ...
$ wdev : int 27 55 53 38 40 56 24 24 31 27 ...
- $ f.wdev : Factor w/ 3 levels "","M","Q": 1 1 1 1 1 1 1 1 1 1 ...
+ $ f.wdev : chr "" "" "" "" ...
$ gspd : num 15.4 7.2 9.6 11.2 12.7 5.8 16.9 10.3 11.1 10.9 ...
- $ f.gspd : Factor w/ 4 levels "","E","M","Q": 1 1 1 1 1 1 1 1 1 1 ...
+ $ f.gspd : chr "" "" "" "" ...
$ s10t : num NA NA NA NA NA NA NA NA NA NA ...
- $ f.s10t : Factor w/ 3 levels "","E","M": 3 3 3 3 3 3 3 3 3 3 ...
+ $ f.s10t : chr "M" "M" "M" "M" ...
$ s10tmax : num NA NA NA NA NA NA NA NA NA NA ...
- $ f.s10tmax: Factor w/ 3 levels "","E","M": 3 3 3 3 3 3 3 3 3 3 ...
+ $ f.s10tmax: chr "M" "M" "M" "M" ...
$ s10tmin : num NA NA NA NA NA NA NA NA NA NA ...
- $ f.s10tmin: Factor w/ 3 levels "","E","M": 3 3 3 3 3 3 3 3 3 3 ...
+ $ f.s10tmin: chr "M" "M" "M" "M" ...
~~~
{: .output}
@@ -471,7 +513,7 @@ derive our NDVI rasters to try to understand what appear to be outlier NDVI valu
> > dimensions : 652, 696, 453792, 3 (nrow, ncol, ncell, nlayers)
> > resolution : 30, 30 (x, y)
> > extent : 230775, 251655, 4704825, 4724385 (xmin, xmax, ymin, ymax)
-> > crs : +proj=utm +zone=19 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0
+> > crs : +proj=utm +zone=19 +datum=WGS84 +units=m +no_defs
> > names : X277_HARV_landRGB.1, X277_HARV_landRGB.2, X277_HARV_landRGB.3
> > min values : 26, 29, 79
> > max values : 255, 255, 255