Project

General

Profile

« Previous | Next » 

Revision 2705fb0e

Added by Benoit Parmentier over 9 years ago

initial commit global assessment presentation slidify, Africa 1992

View differences:

climate/research/oregon/interpolation/presentation_global_run_scalingup_assessment.Rmd
1
---
2
title       : Summary for year 1992, region 4 -Africa
3
subtitle    : Global assessment interpolation
4
author      : Benoit Parmentier and Alberto Guzman
5
job         : NEX job world
6
framework   : io2012        # {io2012, html5slides, shower, dzslides, ...}
7
highlighter : highlight.js  # {highlight.js, prettify, highlight}
8
hitheme     : tomorrow      # 
9
widgets     : []            # {mathjax, quiz, bootstrap}
10
mode        : selfcontained # {standalone, draft}
11
knit        : slidify::knit2slides
12
---
13

  
14
```{r,echo=F}
15

  
16
#First let's set up the input parameters for the script in R:
17

  
18
in_dir <- "/data/project/layers/commons/NEX_data/output_run10_1500x4500_global_analyses_pred_1992_09012015"
19

  
20
#region_name <- "world" #PARAM 13 #reg4 South America, Africa reg5,Europe reg2, North America reg1, Asia reg3
21
#mosaicing_method <- c("unweighted","use_edge_weights")
22
out_suffix <- paste(region_name,"_","pred_1992_09032015",sep="") 
23
#_mosaic_run10_1500x4500_global_analyses_06212015
24
#day_to_mosaic <- c("20100831",
25
#                   "20100901") #PARAM7
26

  
27
```
28

  
29
## Tiles of 1500x4500 km
30

  
31
There were 43 tiles for region 4 (Africa)
32

  
33
1. Tiles were run for 365 days.
34
2. Year tested was 1992.
35
3. Assessment of number of days predicted and accuracy was carried out.
36

  
37

  
38
--- .class #id 
39

  
40
<!-- Limit image width and height -->
41
<style type="text/css">
42
img {     
43
  max-height: 560px;     
44
  max-width: 964px; 
45
}
46
</style>
47
 
48

  
49

  
50
# Tiles of 1500x4500km global run 
51

  
52
```{r, echo=FALSE}
53
#image_file1 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_unweighted_20100831_world_mosaic_07092015")
54

  
55
fname <- paste("Figure1_tile_processed_region_world_run10_1500x4500_global_analyses_pred_1992_09012015.png")
56

  
57
image_file1 <- file.path(in_dir,fname)
58

  
59
```
60

  
61
```{r image_file1, echo = F, results = 'asis'}
62
#| I am text to the left  | ![Flowers](/flowers.jpeg) |
63
#md_str <- paste('\n|        |![](',image_file1,')\n',sep="") 
64
#md_str <- paste('\n->![alt text](',image_file1,'),<-\n',sep="")
65
#<center>![Alt test](http://upload.wikimedia.org/wikipedia/en/b/bc/Wiki.png)</center>
66
md_str <- paste('\n<center>![](',image_file1,')<center>\n',sep="")
67
#md_str <- paste('\n![](',image_file1,'#center)\n',sep="")
68
#fig.align = "center"
69
cat(md_str)
70
```
71

  
72
---
73

  
74
# ACCUARY PER TILES FOR MODEL 1
75

  
76
```{r, echo=FALSE}
77
#image_file2 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_edge_weighted_20100831_world_mosaic_07092015.png")
78

  
79

  
80
fname <- paste("Figure2a_boxplot_with_oultiers_by_tiles_mod1_run10_1500x4500_global_analyses_pred_1992_09012015.png")
81
image_file2 <- file.path(in_dir,fname)
82

  
83
fname <- paste("Figure2b_boxplot_scaling_by_tiles_mod1_run10_1500x4500_global_analyses_pred_1992_09012015.png")
84
image_file3 <- file.path(in_dir,fname)
85

  
86
```
87

  
88

  
89
```{r image_file2, echo = F, results = 'asis', fig.show = 'hold', out.width = '50%',out.extra='style=""'}
90

  
91
md_str <- paste('\n![](',image_file2,') ','![](',image_file3,')\n',sep="") 
92
cat(md_str)
93
#md_str <- paste('\n![](',image_file2,')\n',sep="") 
94
#cat(md_str)
95

  
96
#![alt-text-1](image1.png "title-1") ![alt-text-2](image2.png "title-2")
97

  
98

  
99
```
100

  
101
# ACCUARY OVERALL FOR MODEL 1
102

  
103
```{r, echo=FALSE}
104
#image_file2 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_edge_weighted_20100831_world_mosaic_07092015.png")
105

  
106

  
107
fname <- paste("Figure3a_boxplot_overall_region_with_oultiers_mod_kr_run10_1500x4500_global_analyses_pred_1992_09012015.png")
108
image_file2 <- file.path(in_dir,fname)
109

  
110
fname <- paste("Figure3b_boxplot_overall_region_scaling_mod_kr_run10_1500x4500_global_analyses_pred_1992_09012015.png")
111
image_file3 <- file.path(in_dir,fname)
112

  
113
```
114

  
115

  
116
```{r image_file2, echo = F, results = 'asis', fig.show = 'hold', out.width = '50%',out.extra='style=""'}
117

  
118
md_str <- paste('\n![](',image_file2,') ','![](',image_file3,')\n',sep="") 
119
cat(md_str)
120
#md_str <- paste('\n![](',image_file2,')\n',sep="") 
121
#cat(md_str)
122

  
123
#![alt-text-1](image1.png "title-1") ![alt-text-2](image2.png "title-2")
124

  
125

  
126
```
127

  
128

  
129
---
130

  
131

  
132
#########
climate/research/oregon/interpolation/presentation_summary_world.Rmd
1 1
---
2
title       : New elevation layer test summary world 
2
title       : Summary for year 1992, region 4 -Africa
3 3
subtitle    : Global assessment interpolation
4 4
author      : Benoit Parmentier and Alberto Guzman
5 5
job         : NEX job world
......
15 15

  
16 16
#First let's set up the input parameters for the script in R:
17 17

  
18
in_dir <- "/data/project/layers/commons/NEX_data/output_run10_1500x4500_global_analyses_pred_2010_testelev_07302015"
18
in_dir <- "/data/project/layers/commons/NEX_data/output_run10_1500x4500_global_analyses_pred_1992_09012015"
19 19

  
20 20
#region_name <- "world" #PARAM 13 #reg4 South America, Africa reg5,Europe reg2, North America reg1, Asia reg3
21 21
#mosaicing_method <- c("unweighted","use_edge_weights")
22
out_suffix <- paste(region_name,"_","mosaic_07092015",sep="") 
22
out_suffix <- paste(region_name,"_","pred_1992_09032015",sep="") 
23 23
#_mosaic_run10_1500x4500_global_analyses_06212015
24 24
#day_to_mosaic <- c("20100831",
25 25
#                   "20100901") #PARAM7
26 26

  
27
#png_files <- mixedsort(list.files(path=in_dir,pattern=paste("*",out_suffix,".png",sep="")))
28

  
29 27
```
30 28

  
31
## Testing elevation layers and masks
32

  
29
## Tiles of 1500x4500 km
33 30

  
34
Two tiles were predicted for 10 dates:
31
There were 43 tiles for region 4 (Africa)
35 32

  
36
1. Combined mask no buffer
37
2. Buffered 1km
38
3. Buffered 2km
33
1. Tiles were run for 365 days.
34
2. Year tested was 1992.
35
3. Assessment of number of days predicted and accuracy was carried out.
39 36

  
40 37

  
41 38
--- .class #id 
......
48 45
}
49 46
</style>
50 47
 
51
<!-- Center image on slide -->
52
<script type="text/javascript" src="http://ajax.aspnetcdn.com/ajax/jQuery/jquery-1.7.min.js"></script>
53
<script type="text/javascript">
54
$(function() {     
55
  $("p:has(img)").addClass('centered'); 
56
});
57
</script>
58 48

  
59
# Combined mask Asia: 201000101
49

  
50
# Tiles of 1500x4500km global run 
60 51

  
61 52
```{r, echo=FALSE}
62 53
#image_file1 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_unweighted_20100831_world_mosaic_07092015")
63 54

  
64
fname <- paste("Figure10_combined_combined_20100101_30_111.6_94.1_1500x4500_run10_1500x4500_global_analyses_pred_2010_testelev_07302015.png")
55
fname <- paste("Figure1_tile_processed_region_world_run10_1500x4500_global_analyses_pred_1992_09012015.png")
65 56

  
66 57
image_file1 <- file.path(in_dir,fname)
67 58

  
68 59
```
69 60

  
70 61
```{r image_file1, echo = F, results = 'asis'}
71
md_str <- paste('\n![](',image_file1,')\n',sep="") 
62
#| I am text to the left  | ![Flowers](/flowers.jpeg) |
63
#md_str <- paste('\n|        |![](',image_file1,')\n',sep="") 
64
#md_str <- paste('\n->![alt text](',image_file1,'),<-\n',sep="")
65
#<center>![Alt test](http://upload.wikimedia.org/wikipedia/en/b/bc/Wiki.png)</center>
66
md_str <- paste('\n<center>![](',image_file1,')<center>\n',sep="")
67
#md_str <- paste('\n![](',image_file1,'#center)\n',sep="")
68
#fig.align = "center"
72 69
cat(md_str)
73 70
```
74 71

  
75 72
---
76 73

  
77
# Buffered 1km mask Asia: 201000101
74
# ACCUARY PER TILES FOR MODEL 1
78 75

  
79 76
```{r, echo=FALSE}
80 77
#image_file2 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_edge_weighted_20100831_world_mosaic_07092015.png")
81 78

  
82 79

  
83
fname <- paste("Figure10_1kmBuff_1kmBuff_20100101_30_111.6_94.1_1500x4500_run10_1500x4500_global_analyses_pred_2010_testelev_07302015.png")
80
fname <- paste("Figure2a_boxplot_with_oultiers_by_tiles_mod1_run10_1500x4500_global_analyses_pred_1992_09012015.png")
84 81
image_file2 <- file.path(in_dir,fname)
85 82

  
86
```
87

  
88
```{r image_file2, echo = F, results = 'asis'}
89
md_str <- paste('\n![](',image_file2,')\n',sep="") 
90
cat(md_str)
91
```
92

  
93
---
94

  
95
# Buffered 2km mask Asia: 201000101
96

  
97
```{r, echo=FALSE}
98
#image_file3 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_unweighted_20100901_world_mosaic_07092015.png")
99

  
100
fname <- paste("Figure10_2kmBuff_2kmBuff_20100101_30_111.6_94.1_1500x4500_run10_1500x4500_global_analyses_pred_2010_testelev_07302015.png")
101

  
83
fname <- paste("Figure2b_boxplot_scaling_by_tiles_mod1_run10_1500x4500_global_analyses_pred_1992_09012015.png")
102 84
image_file3 <- file.path(in_dir,fname)
103 85

  
104

  
105
```
106

  
107
```{r image_file3, echo = F, results = 'asis'}
108
md_str <- paste('\n![](',image_file3,')\n',sep="") 
109
cat(md_str)
110
```
111

  
112
---
113

  
114
# Combined mask Asia: 201000901
115

  
116
```{r, echo=FALSE}
117

  
118
fname <- paste("Figure10_combined_combined_20100901_30_111.6_94.1_1500x4500_run10_1500x4500_global_analyses_pred_2010_testelev_07302015.png")
119
image_file4 <- file.path(in_dir,fname)
120

  
121 86
```
122 87

  
123
```{r image_file4, echo = F, results = 'asis'}
124
md_str <- paste('\n![](',image_file4,')\n',sep="") 
125
cat(md_str)
126
```
127 88

  
128
---
129

  
130
# Buffered 1km mask Asia: 201000901
131

  
132
```{r, echo=FALSE}
89
```{r image_file2, echo = F, results = 'asis', fig.show = 'hold', out.width = '50%',out.extra='style=""'}
133 90

  
134
fname <- paste("Figure10_1kmBuff_1kmBuff_20100901_30_111.6_94.1_1500x4500_run10_1500x4500_global_analyses_pred_2010_testelev_07302015.png")
135
image_file5 <- file.path(in_dir,fname)
136

  
137
```
138

  
139
```{r image_file5, echo = F, results = 'asis'}
140
md_str <- paste('\n![](',image_file5,')\n',sep="") 
91
md_str <- paste('\n![](',image_file2,') ','![](',image_file3,')\n',sep="") 
141 92
cat(md_str)
142
```
93
#md_str <- paste('\n![](',image_file2,')\n',sep="") 
94
#cat(md_str)
143 95

  
144
---
96
#![alt-text-1](image1.png "title-1") ![alt-text-2](image2.png "title-2")
145 97

  
146
# Buffered 2km mask Asia: 201000901
147 98

  
148
```{r, echo=FALSE}
149

  
150
fname <- paste("Figure10_2kmBuff_2kmBuff_20100901_30_111.6_94.1_1500x4500_run10_1500x4500_global_analyses_pred_2010_testelev_07302015.png")
151

  
152
image_file6 <- file.path(in_dir,fname)
153

  
154
```
155

  
156
```{r image_file6, echo = F, results = 'asis'}
157
md_str <- paste('\n![](',image_file6,')\n',sep="") 
158
cat(md_str)
159 99
```
160 100

  
161
---
162

  
163
# Combined mask Europe: 201000101
164

  
165
```{r, echo=FALSE}
166
#image_file1 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_unweighted_20100831_world_mosaic_07092015")
167

  
168
fname <- paste("Figure10_combined_combined_20100101_30_151.5_23.7_1500x4500_run10_1500x4500_global_analyses_pred_2010_testelev_07302015.png")
169

  
170
image_file7 <- file.path(in_dir,fname)
171

  
172
```
173

  
174
```{r image_file7, echo = F, results = 'asis'}
175
md_str <- paste('\n![](',image_file7,')\n',sep="") 
176
cat(md_str)
177
```
178

  
179
---
180

  
181
# Buffered 1km mask Europe: 201000101
101
# ACCUARY OVERALL FOR MODEL 1
182 102

  
183 103
```{r, echo=FALSE}
184 104
#image_file2 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_edge_weighted_20100831_world_mosaic_07092015.png")
185 105

  
186 106

  
187
fname <- paste("Figure10_1kmBuff_1kmBuff_20100101_30_151.5_23.7_1500x4500_run10_1500x4500_global_analyses_pred_2010_testelev_07302015.png")
188
image_file8 <- file.path(in_dir,fname)
189

  
190
```
191

  
192
```{r image_file8, echo = F, results = 'asis'}
193
md_str <- paste('\n![](',image_file8,')\n',sep="") 
194
cat(md_str)
195
```
196

  
197
---
198

  
199
# Buffered 2km mask Europe: 201000101
200

  
201
```{r, echo=FALSE}
202

  
203
fname <- paste("Figure10_2kmBuff_2kmBuff_20100101_30_151.5_23.7_1500x4500_run10_1500x4500_global_analyses_pred_2010_testelev_07302015.png")
204

  
205
image_file9 <- file.path(in_dir,fname)
206

  
207

  
208
```
209

  
210
```{r image_file9, echo = F, results = 'asis'}
211
md_str <- paste('\n![](',image_file9,')\n',sep="") 
212
cat(md_str)
213
```
214
---
215

  
216
# Combined mask Europe: 201000901
217

  
218
```{r, echo=FALSE}
219

  
220
fname <- paste("Figure10_combined_combined_20100901_30_151.5_23.7_1500x4500_run10_1500x4500_global_analyses_pred_2010_testelev_07302015.png")
221

  
222
image_file10 <- file.path(in_dir,fname)
223

  
224
```
225

  
226
```{r image_file10, echo = F, results = 'asis'}
227
md_str <- paste('\n![](',image_file10,')\n',sep="") 
228
cat(md_str)
229
```
230

  
231
---
232

  
233
# Buffered 1km mask Europe: 201000901
234

  
235
```{r, echo=FALSE}
236

  
237
fname <- paste("Figure10_1kmBuff_1kmBuff_20100901_30_151.5_23.7_1500x4500_run10_1500x4500_global_analyses_pred_2010_testelev_07302015.png")
238
image_file11 <- file.path(in_dir,fname)
239

  
240
```
241

  
242
```{r image_file11, echo = F, results = 'asis'}
243
md_str <- paste('\n![](',image_file11,')\n',sep="") 
244
cat(md_str)
245
```
246

  
247
---
248

  
249
# Buffered 2km mask Europe: 201000901
250

  
251
```{r, echo=FALSE}
252

  
253
fname <- paste("Figure10_2kmBuff_2kmBuff_20100901_30_151.5_23.7_1500x4500_run10_1500x4500_global_analyses_pred_2010_testelev_07302015.png")
254

  
255
image_file12 <- file.path(in_dir,fname)
256

  
257

  
258
```
259

  
260
```{r image_file12, echo = F, results = 'asis'}
261
md_str <- paste('\n![](',image_file12,')\n',sep="") 
262
cat(md_str)
263
```
264

  
265
---
266

  
267
# Boxplot Asia and Europe summarized
268

  
269
```{r, echo=FALSE}
270

  
271
fname <- paste("Figure8a_boxplot_overall_separated_by_region_with_oultiers_NA_run10_1500x4500_global_analyses_pred_2010_testelev_07302015.png")
272

  
273
image_file13 <- file.path(in_dir,fname)
274

  
275

  
276
```
277
```{r image_file13, echo = F, results = 'asis'}
278
md_str <- paste('\n![](',image_file13,')\n',sep="") 
279
cat(md_str)
280
```
281
---
282
---
283

  
284

  
285
# Combined Zoom in 1: Asia 
286

  
287
```{r, echo=FALSE}
288

  
289
fname <- paste("zoom_in_1_vietnam_combined_20100101.png")
290

  
291
image_file14 <- file.path(in_dir,fname)
292

  
293

  
294
```
295
```{r image_file14, echo = F, results = 'asis'}
296
md_str <- paste('\n![](',image_file14,')\n',sep="") 
297
cat(md_str)
298
```
299

  
300
---
301

  
302

  
303
# Combined Zoom in 1: Asia 
304

  
305
```{r, echo=FALSE}
306

  
307
fname <- paste("zoom_in_1_vietnam_combined_20100101.png")
308

  
309
image_file14b <- file.path(in_dir,fname)
310

  
311

  
312
```
313
```{r image_file14b, echo = F, results = 'asis'}
314
md_str <- paste('\n![](',image_file14b,')\n',sep="") 
315
cat(md_str)
316
```
317

  
318
---
319

  
320
# 1km Buffer Zoom in 1: Asia 
321

  
322
```{r, echo=FALSE}
323

  
324
fname <- paste("zoom_in_1_vietnam_1kmbuff_20100101.png")
325

  
326
image_file15 <- file.path(in_dir,fname)
107
fname <- paste("Figure3a_boxplot_overall_region_with_oultiers_mod_kr_run10_1500x4500_global_analyses_pred_1992_09012015.png")
108
image_file2 <- file.path(in_dir,fname)
327 109

  
110
fname <- paste("Figure3b_boxplot_overall_region_scaling_mod_kr_run10_1500x4500_global_analyses_pred_1992_09012015.png")
111
image_file3 <- file.path(in_dir,fname)
328 112

  
329
```
330
```{r image_file15, echo = F, results = 'asis'}
331
md_str <- paste('\n![](',image_file15,')\n',sep="") 
332
cat(md_str)
333 113
```
334 114

  
335
---
336

  
337
# Buffer 2km Zoom in 1: Asia 
338

  
339
```{r, echo=FALSE}
340

  
341
fname <- paste("zoom_in_1_vietnam_2kmbuff_20100101.png")
342

  
343
image_file16 <- file.path(in_dir,fname)
344 115

  
116
```{r image_file2, echo = F, results = 'asis', fig.show = 'hold', out.width = '50%',out.extra='style=""'}
345 117

  
346
```
347
```{r image_file16, echo = F, results = 'asis'}
348
md_str <- paste('\n![](',image_file16,')\n',sep="") 
118
md_str <- paste('\n![](',image_file2,') ','![](',image_file3,')\n',sep="") 
349 119
cat(md_str)
350
```
351

  
352
---
353

  
120
#md_str <- paste('\n![](',image_file2,')\n',sep="") 
121
#cat(md_str)
354 122

  
355
# Combined Zoom in 2: Europe
356

  
357
```{r, echo=FALSE}
358

  
359
fname <- paste("zoom_in_2_Netherlands_combined_20100101.png")
360

  
361
image_file17 <- file.path(in_dir,fname)
123
#![alt-text-1](image1.png "title-1") ![alt-text-2](image2.png "title-2")
362 124

  
363 125

  
364 126
```
365
```{r image_file17, echo = F, results = 'asis'}
366
md_str <- paste('\n![](',image_file17,')\n',sep="") 
367
cat(md_str)
368
```
369

  
370
---
371

  
372
# 1km Buffer Zoom in 2: Europe 
373

  
374
```{r, echo=FALSE}
375

  
376
fname <- paste("zoom_in_2_Netherlands_1kmbuff_20100101.png")
377

  
378
image_file18 <- file.path(in_dir,fname)
379 127

  
380 128

  
381
```
382
```{r image_file18, echo = F, results = 'asis'}
383
md_str <- paste('\n![](',image_file18,')\n',sep="") 
384
cat(md_str)
385
```
386

  
387 129
---
388 130

  
389
# Buffer 2km Zoom in 2: Europe
390

  
391
```{r, echo=FALSE}
392

  
393
fname <- paste("zoom_in_2_Netherlands_2kmbuff_20100101.png")
394

  
395
image_file19 <- file.path(in_dir,fname)
396

  
397 131

  
398
```
399
```{r image_file19, echo = F, results = 'asis'}
400
md_str <- paste('\n![](',image_file19,')\n',sep="") 
401
cat(md_str)
402
```
132
#########

Also available in: Unified diff