Project

General

Profile

« Previous | Next » 

Revision 0ad07087

Added by Benoit Parmentier over 9 years ago

slidify presentation mosaicing assessment using weighted mosaicing from edge

View differences:

climate/research/oregon/interpolation/mosaicing_presentation_summary.Rmd
1 1
---
2
title       : Mosaicing test summary region 23, Europe and Asia
2
title       : Mosaicing test summary world mosaic
3 3
subtitle    : Global assessment interpolation
4 4
author      : Benoit Parmentier
5
job         : NEX job reg23
5
job         : NEX job world
6 6
framework   : io2012        # {io2012, html5slides, shower, dzslides, ...}
7 7
highlighter : highlight.js  # {highlight.js, prettify, highlight}
8 8
hitheme     : tomorrow      # 
......
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/mosaicing_data_test/reg23/output_reg23_mosaic_run10_1500x4500_global_analyses_06292015"
18
in_dir <- "/data/project/layers/commons/NEX_data/mosaicing_data_test/mosaic_world_07012015"
19 19

  
20
region_name <- "reg23" #PARAM 13 #reg4 South America, Africa reg5,Europe reg2, North America reg1, Asia reg3
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_run10_1500x4500_global_analyses_06292015",sep="") 
22
out_suffix <- paste(region_name,"_","mosaic_07092015",sep="") 
23 23
#_mosaic_run10_1500x4500_global_analyses_06212015
24 24
day_to_mosaic <- c("20100831",
25 25
                   "20100901") #PARAM7
......
30 30

  
31 31
## Mosaicing with two methods
32 32

  
33
All regions were merged together for the two mosaicing methods:
33 34

  
34 35
1. Unweighted mean: each image a pixel has no specific weight
35 36
2. distance from edge: weights of pixels depends on distance to edges for each tiles
36 37

  
37
This was done for two dates:
38
The world mosaic was produced for two dates:
38 39
- 20100831
39 40
- 20100901
40 41

  
41 42
---
42 43

  
43
## Region 23
44

  
45
This presentation display mosaicing results for reg23: 
46

  
47
This regions contains the following number of tiles:
48
```{r}
49
pattern_str<- paste("*",mosaicing_method[2],"_m_",day_to_mosaic[1],"_",out_suffix,".tif",sep="")
50
length(list.files(path=in_dir,pattern=pattern_str))
51
```
52

  
53 44
--- .class #id 
54 45

  
55 46
<!-- Limit image width and height -->
......
71 62
# Unweighted mean reg: 20100831
72 63

  
73 64
```{r, echo=FALSE}
74
#image_file1 <- file.path(in_dir,"Figure2_mosaic_mean_unweighted_20100831_unweighted_20100831_20100831_reg6_mosaic_r#un10_1500x4500_global_analyses_06212015.png")
65
#image_file1 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_unweighted_20100831_world_mosaic_07092015")
75 66

  
76
fname <- paste("Figure2_mosaic_mean_unweighted_20100831_unweighted_20100831_20100831_",out_suffix,".png",sep="")
67
fname <- paste("Figure9_clim_mosaics_day_test_unweighted_20100831_",out_suffix,".png",sep="")
77 68

  
78 69
image_file1 <- file.path(in_dir,fname)
79 70

  
......
89 80
# Edge weigted mean :20100831
90 81

  
91 82
```{r, echo=FALSE}
92
#image_file2 <- file.path(in_dir,"Figure2_mosaic_mean_use_edge_weights_20100831_use_edge_weights_20100831_20100831_r#eg4_mosaic_run10_1500x4500_global_analyses_06212015.png")
83
#image_file2 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_edge_weighted_20100831_world_mosaic_07092015.png")
93 84

  
94 85

  
95
fname <- paste("Figure2_mosaic_mean_use_edge_weights_20100831_use_edge_weights_20100831_20100831_",out_suffix,".png",sep="")
86
fname <- paste("Figure9_clim_mosaics_day_test_edge_weighted_20100831_",out_suffix,".png",sep="")
96 87
image_file2 <- file.path(in_dir,fname)
97 88

  
98 89
```
......
107 98
# Unweighted mean : 20100901
108 99

  
109 100
```{r, echo=FALSE}
110
#image_file3 <- file.path(in_dir,"Figure2_mosaic_mean_unweighted_20100901_unweighted_20100901_20100901_reg2_mosaic_r#un10_1500x4500_global_analyses_06212015.png")
101
#image_file3 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_unweighted_20100901_world_mosaic_07092015.png")
111 102

  
112
fname <- paste("Figure2_mosaic_mean_unweighted_20100901_unweighted_20100901_20100901_",out_suffix,".png",sep="")
103
fname <- paste("Figure9_clim_mosaics_day_test_unweighted_20100901_",out_suffix,".png",sep="")
113 104

  
114 105
image_file3 <- file.path(in_dir,fname)
115 106

  
......
128 119
# Edge weigted mean : 20100901
129 120

  
130 121
```{r, echo=FALSE}
131
#image_file4 <- file.path(in_dir,"Figure2_mosaic_mean_use_edge_weights_20100901_use_edge_weights_20100901_20100901_r#eg4_mosaic_run10_1500x4500_global_analyses_06212015.png")
122
#image_file4 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_edge_weighted_20100901_world_mosaic_07092015.png")
132 123

  
133
fname <- paste("Figure2_mosaic_mean_use_edge_weights_20100901_use_edge_weights_20100901_20100901_",out_suffix,".png",sep="")
124
fname <- paste("Figure9_clim_mosaics_day_test_edge_weighted_20100901_",out_suffix,".png",sep="")
134 125
image_file4 <- file.path(in_dir,fname)
135 126

  
136 127
```
......
142 133

  
143 134
---
144 135

  
145
# Unweigted mean  slope: 20100831
146

  
147
```{r, echo=FALSE}
148
fname <- paste("Figure2_slope_mean_unweighted_20100831_unweighted_20100831_20100831_",out_suffix,".png",sep="")
149
image_file5 <- file.path(in_dir,fname)
150
```
151

  
152
```{r image_file5, echo = F, results = 'asis'}
153
#should make this a function!!!
154
md_str <- paste('\n![](',image_file5,')\n',sep="") 
155
cat(md_str)
156
```
157

  
158
---
159

  
160
# Edge weigted mean slope: 20100831
161

  
162
```{r, echo=FALSE}
163
#image_file6 <- file.path(in_dir,"Figure2_slope_mean_use_edge_weights_20100831_use_edge_weights_20100831_20100831_re#g4_mosaic_run10_1500x4500_global_analyses_06212015.png")
164
fname <- paste("Figure2_slope_mean_use_edge_weights_20100831_use_edge_weights_20100831_20100831_",out_suffix,".png",sep="")
165
image_file6 <- file.path(in_dir,fname)
166

  
167
```
168

  
169
```{r image_file6, echo = F, results = 'asis'}
170
md_str <- paste('\n![](',image_file6,')\n',sep="") 
171
cat(md_str)
172
```
173

  
174
---
175

  
176
# Unweigted mean slope: 20100901
177

  
178
```{r, echo=FALSE}
179
fname <- paste("Figure2_slope_mean_unweighted_20100901_unweighted_20100901_20100901_",out_suffix,".png",sep="")
180

  
181
image_file7 <- file.path(in_dir,fname)
182

  
183
```
184

  
185
```{r image_file7, echo = F, results = 'asis'}
186
md_str <- paste('\n![](',image_file7,')\n',sep="") 
187
cat(md_str)
188
```
189

  
190
---
191

  
192
# Edge weigted mean slope: 20100901
193

  
194

  
195
```{r, echo=FALSE}
196
#fname <- paste("Figure2_slope_mean_use_edge_weights_20100901_use_edge_weights_20100901_20100901_",out_suffix,".png"#,sep="")
197
#fname <- "Figure2_slope_mean_use_edge_weights_20100901_use_edge_weights_20100901_20100901_reg4_mosaic_run10_1500x45#00_global_analyses_06212015.png"
198
fname <- paste("Figure2_slope_mean_use_edge_weights_20100901_use_edge_weights_20100901_20100901_",out_suffix,".png",sep="")
199
image_file8 <- file.path(in_dir,fname)
200

  
201
```
202

  
203
```{r image_file8, echo = F, results = 'asis'}
204
md_str <- paste('\n![](',image_file8,')\n',sep="") 
205
cat(md_str)
206
```
207

  
208
---
209

  
210
# Unweigted mean aspect: 20100831
211

  
212
```{r, echo=FALSE}
213

  
214
fname <- paste("Figure2_aspect_mean_unweighted_20100831_unweighted_20100831_20100831_",out_suffix,".png",sep="")
215
image_file9 <- file.path(in_dir,fname)
216

  
217
```
218

  
219
```{r image_file9, echo = F, results = 'asis'}
220
md_str <- paste('\n![](',image_file9,')\n',sep="") 
221
cat(md_str)
222
```
223

  
224
---
225

  
226
# Edge weigted mean aspect: 20100831
227

  
228
```{r, echo=FALSE}
229
#image_file10 <- file.path(in_dir,"Figure2_aspect_mean_use_edge_weights_20100831_use_edge_weights_20100831_20100831_#reg6_mosaic_run10_1500x4500_global_analyses_06212015.png")
230
fname <- paste("Figure2_aspect_mean_use_edge_weights_20100831_use_edge_weights_20100831_20100831_",out_suffix,".png",sep="")
231
image_file10 <- file.path(in_dir,fname)
232

  
233
```
234

  
235
```{r image_file10, echo = F, results = 'asis'}
236
md_str <- paste('\n![](',image_file10,')\n',sep="") 
237
cat(md_str)
238
```
239

  
240
---
241

  
242
# Unweigted mean aspect: 20100901
243

  
244
```{r, echo=FALSE}
245
#image_file11 <- file.path(in_dir,"Figure2_aspect_mean_unweighted_20100901_unweighted_20100901_20100901_reg6_mosaic_#run10_1500x4500_global_analyses_06212015.png")
246
fname <- paste("Figure2_aspect_mean_unweighted_20100901_unweighted_20100901_20100901_",out_suffix,".png",sep="")
247
image_file11 <- file.path(in_dir,fname)
248

  
249
```
250

  
251
```{r image_file11, echo = F, results = 'asis'}
252
md_str <- paste('\n![](',image_file11,')\n',sep="") 
253
cat(md_str)
254
```
255

  
256
---
257

  
258
# Edge weigted mean aspect: 20100901
259

  
260
```{r, echo=FALSE}
261
#image_file12 <- file.path(in_dir,"Figure2_aspect_mean_use_edge_weights_20100901_use_edge_weights_20100901_20100901_#reg6_mosaic_run10_1500x4500_global_analyses_06212015.png")
262
fname <- paste("Figure2_aspect_mean_use_edge_weights_20100901_use_edge_weights_20100901_20100901_",out_suffix,".png",sep="")
263
image_file12 <- file.path(in_dir,fname)
264
```
265

  
266
```{r image_file12, echo = F, results = 'asis'}
267
md_str <- paste('\n![](',image_file12,')\n',sep="") 
268
cat(md_str)
269
```
270

  
271
---
272

  
273
# Unweigted - edge weighted diff: 20100831
274

  
275
```{r, echo=FALSE}
276
#image_file13 <- file.path(in_dir,"Figure2_diff_raster_unweighted_use_edge_weights_20100831_reg6_mosaic_run10_1500x4#500_global_analyses_06212015.png")
277

  
278
fname <- paste("Figure2_diff_raster_unweighted_use_edge_weights_20100831_",out_suffix,".png",sep="")
279
image_file13 <- file.path(in_dir,fname)
280
```
281

  
282
```{r image_file13, echo = F, results = 'asis'}
283
md_str <- paste('\n![](',image_file13,')\n',sep="") 
284
cat(md_str)
285
```
286

  
287
---
288

  
289
# Unweigted - edge weighted diff : 20100901
290

  
291
```{r, echo=FALSE}
292
#image_file14 <- file.path(in_dir,"Figure2_diff_raster_unweighted_use_edge_weights_20100901_reg6_mosaic_run10_1500x4#500_global_analyses_06212015.png")
293

  
294
fname <- paste("Figure2_diff_raster_unweighted_use_edge_weights_20100901_",out_suffix,".png",sep="")
295
image_file14 <- file.path(in_dir,fname)
296

  
297
```
298

  
299
```{r image_file14, echo = F, results = 'asis'}
300
md_str <- paste('\n![](',image_file14,')\n',sep="") 
301
cat(md_str)
302
```
303

  

Also available in: Unified diff