Revision eac2e440
Added by Benoit Parmentier over 9 years ago
| climate/research/oregon/interpolation/mosaicing_presentation_summary.Rmd | ||
|---|---|---|
| 1 | 1 |
--- |
| 2 |
title : Mosaicing test summary |
|
| 2 |
title : Mosaicing test summary region 4,South America
|
|
| 3 | 3 |
subtitle : Global assessment interpolation |
| 4 | 4 |
author : Benoit Parmentier |
| 5 |
job : NEX job reg6
|
|
| 5 |
job : NEX job reg4
|
|
| 6 | 6 |
framework : io2012 # {io2012, html5slides, shower, dzslides, ...}
|
| 7 | 7 |
highlighter : highlight.js # {highlight.js, prettify, highlight}
|
| 8 | 8 |
hitheme : tomorrow # |
| ... | ... | |
| 12 | 12 |
--- |
| 13 | 13 |
|
| 14 | 14 |
```{r,echo=F}
|
| 15 |
#First let's set up the input director |
|
| 16 |
in_dir <- "/data/project/layers/commons/NEX_data/mosaicing_data_test/reg6/output_reg6_mosaic_run10_1500x4500_global_analyses_06212015" |
|
| 17 | 15 |
|
| 18 |
region_name <- "reg6" #PARAM 13 #reg4 South America, Africa reg5,Europe reg2, North America reg1, Asia reg3 |
|
| 16 |
#First let's set up the input parameters for the script in R: |
|
| 17 |
|
|
| 18 |
in_dir <- "/data/project/layers/commons/NEX_data/mosaicing_data_test/reg4/output_reg4_mosaic_run10_1500x4500_global_analyses_06212015" |
|
| 19 |
|
|
| 20 |
region_name <- "reg4" #PARAM 13 #reg4 South America, Africa reg5,Europe reg2, North America reg1, Asia reg3 |
|
| 19 | 21 |
mosaicing_method <- c("unweighted","use_edge_weights")
|
| 20 | 22 |
out_suffix <- paste(region_name,"_","mosaic_run10_1500x4500_global_analyses_06212015",sep="") |
| 23 |
day_to_mosaic <- c("20100831",
|
|
| 24 |
"20100901") #PARAM7 |
|
| 21 | 25 |
|
| 22 |
png_files <- list.files(path=in_dir,pattern=paste("*",out_suffix,".png",sep=""))
|
|
| 26 |
png_files <- mixedsort(list.files(path=in_dir,pattern=paste("*",out_suffix,".png",sep="")))
|
|
| 23 | 27 |
|
| 24 | 28 |
``` |
| 25 | 29 |
|
| ... | ... | |
| 35 | 39 |
|
| 36 | 40 |
--- |
| 37 | 41 |
|
| 38 |
## Region 6
|
|
| 42 |
## Region 4
|
|
| 39 | 43 |
|
| 40 |
This presentation display mosaicing results for reg6: Oceania and South East Asia
|
|
| 44 |
This presentation display mosaicing results for reg4: South America
|
|
| 41 | 45 |
|
| 42 |
Show number of tiles:
|
|
| 46 |
This regions contains the following number of tiles:
|
|
| 43 | 47 |
```{r}
|
| 44 |
length(list.files(path=in_dir,pattern="*_use_edge_prod_weights_m_20100901_reg6_mosaic_run10_1500x4500_global_analyses_06212015.tif")) |
|
| 48 |
pattern_str<- paste("*",mosaicing_method[2],"_m_",day_to_mosaic[1],"_",out_suffix,".tif",sep="")
|
|
| 49 |
length(list.files(path=in_dir,pattern=pattern_str)) |
|
| 45 | 50 |
``` |
| 46 | 51 |
|
| 47 | 52 |
--- .class #id |
| ... | ... | |
| 65 | 70 |
# Unweighted mean reg: 20100831 |
| 66 | 71 |
|
| 67 | 72 |
```{r, echo=FALSE}
|
| 68 |
image_file1 <- file.path(in_dir,"Figure2_mosaic_mean_unweighted_20100831_unweighted_20100831_20100831_reg6_mosaic_run10_1500x4500_global_analyses_06212015.png") |
|
| 73 |
#image_file1 <- file.path(in_dir,"Figure2_mosaic_mean_unweighted_20100831_unweighted_20100831_20100831_reg6_mosaic_r#un10_1500x4500_global_analyses_06212015.png") |
|
| 74 |
|
|
| 75 |
image_file1 <- file.path(in_dir,paste("Figure2_mosaic_mean_unweighted_20100831_unweighted_20100831_20100831_",out_suffix,".png",sep=""))
|
|
| 76 |
|
|
| 69 | 77 |
``` |
| 70 | 78 |
|
| 71 | 79 |
```{r image_file1, echo = F, results = 'asis'}
|
| 72 |
md_str <- paste('\n\n',sep="")
|
|
| 80 |
md_str <- paste('\n\n',sep="")
|
|
| 73 | 81 |
cat(md_str) |
| 74 | 82 |
``` |
| 75 | 83 |
|
| ... | ... | |
| 78 | 86 |
# Edge weigted mean :20100831 |
| 79 | 87 |
|
| 80 | 88 |
```{r, echo=FALSE}
|
| 81 |
image_file2 <- file.path(in_dir,"Figure2_mosaic_mean_use_edge_weights_20100831_use_edge_weights_20100831_20100831_reg6_mosaic_run10_1500x4500_global_analyses_06212015.png")
|
|
| 89 |
image_file2 <- file.path(in_dir,"Figure2_mosaic_mean_use_edge_weights_20100831_use_edge_weights_20100831_20100831_reg4_mosaic_run10_1500x4500_global_analyses_06212015.png")
|
|
| 82 | 90 |
``` |
| 83 | 91 |
|
| 84 | 92 |
```{r image_file2, echo = F, results = 'asis'}
|
| 85 |
md_str <- paste('\n\n',sep="")
|
|
| 93 |
md_str <- paste('\n\n',sep="")
|
|
| 86 | 94 |
cat(md_str) |
| 87 | 95 |
``` |
| 88 | 96 |
|
| ... | ... | |
| 91 | 99 |
# Unweighted mean : 20100901 |
| 92 | 100 |
|
| 93 | 101 |
```{r, echo=FALSE}
|
| 94 |
image_file3 <- file.path(in_dir,"Figure2_mosaic_mean_unweighted_20100901_unweighted_20100901_20100901_reg6_mosaic_run10_1500x4500_global_analyses_06212015.png")
|
|
| 102 |
image_file3 <- file.path(in_dir,"Figure2_mosaic_mean_unweighted_20100901_unweighted_20100901_20100901_reg4_mosaic_run10_1500x4500_global_analyses_06212015.png")
|
|
| 95 | 103 |
``` |
| 96 | 104 |
|
| 97 | 105 |
```{r image_file3, echo = F, results = 'asis'}
|
| 98 |
md_str <- paste('\n\n',sep="")
|
|
| 106 |
md_str <- paste('\n\n',sep="")
|
|
| 99 | 107 |
cat(md_str) |
| 100 | 108 |
``` |
| 101 | 109 |
|
| ... | ... | |
| 104 | 112 |
# Edge weigted mean : 20100901 |
| 105 | 113 |
|
| 106 | 114 |
```{r, echo=FALSE}
|
| 107 |
image_file4 <- file.path(in_dir,"Figure2_mosaic_mean_use_edge_weights_20100901_use_edge_weights_20100901_20100901_reg6_mosaic_run10_1500x4500_global_analyses_06212015.png")
|
|
| 115 |
image_file4 <- file.path(in_dir,"Figure2_mosaic_mean_use_edge_weights_20100901_use_edge_weights_20100901_20100901_reg4_mosaic_run10_1500x4500_global_analyses_06212015.png")
|
|
| 108 | 116 |
``` |
| 109 | 117 |
|
| 110 | 118 |
```{r image_file4, echo = F, results = 'asis'}
|
| 111 |
md_str <- paste('\n\n',sep="")
|
|
| 119 |
md_str <- paste('\n\n',sep="")
|
|
| 112 | 120 |
cat(md_str) |
| 113 | 121 |
``` |
| 114 | 122 |
|
| ... | ... | |
| 117 | 125 |
# Unweigted mean slope: 20100831 |
| 118 | 126 |
|
| 119 | 127 |
```{r, echo=FALSE}
|
| 120 |
image_file5 <- file.path(in_dir,"Figure2_slope_mean_unweighted_20100831_unweighted_20100831_20100831_reg6_mosaic_run10_1500x4500_global_analyses_06212015.png") |
|
| 128 |
fname <- paste("Figure2_slope_mean_unweighted_20100831_unweighted_20100831_20100831_",out_suffix,".png",sep="")
|
|
| 129 |
image_file5 <- file.path(in_dir,fname) |
|
| 121 | 130 |
``` |
| 122 | 131 |
|
| 123 | 132 |
```{r image_file5, echo = F, results = 'asis'}
|
| 124 |
md_str <- paste('\n\n',sep="")
|
|
| 133 |
#should make this a function!!! |
|
| 134 |
md_str <- paste('\n\n',sep="")
|
|
| 125 | 135 |
cat(md_str) |
| 126 | 136 |
``` |
| 127 | 137 |
|
| ... | ... | |
| 130 | 140 |
# Edge weigted mean slope: 20100831 |
| 131 | 141 |
|
| 132 | 142 |
```{r, echo=FALSE}
|
| 133 |
image_file6 <- file.path(in_dir,"Figure2_slope_mean_use_edge_weights_20100831_use_edge_weights_20100831_20100831_reg6_mosaic_run10_1500x4500_global_analyses_06212015.png") |
|
| 143 |
#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") |
|
| 144 |
fname <- paste("Figure2_slope_mean_use_edge_weights_20100831_use_edge_weights_20100831_20100831_",out_suffix,".png",sep="")
|
|
| 145 |
image_file6 <- file.path(in_dir,fname) |
|
| 146 |
|
|
| 134 | 147 |
``` |
| 135 | 148 |
|
| 136 | 149 |
```{r image_file6, echo = F, results = 'asis'}
|
| 137 |
md_str <- paste('\n\n',sep="")
|
|
| 150 |
md_str <- paste('\n\n',sep="")
|
|
| 138 | 151 |
cat(md_str) |
| 139 | 152 |
``` |
| 140 | 153 |
|
| ... | ... | |
| 143 | 156 |
# Unweigted mean slope: 20100901 |
| 144 | 157 |
|
| 145 | 158 |
```{r, echo=FALSE}
|
| 146 |
image_file7 <- file.path(in_dir,"Figure2_slope_mean_unweighted_20100901_unweighted_20100901_20100901_reg6_mosaic_run10_1500x4500_global_analyses_06212015.png") |
|
| 159 |
fname <- paste("Figure2_slope_mean_unweighted_20100901_unweighted_20100901_20100901_",out_suffix,".png",sep="")
|
|
| 160 |
|
|
| 161 |
image_file7 <- file.path(in_dir,fname) |
|
| 162 |
|
|
| 147 | 163 |
``` |
| 148 | 164 |
|
| 149 | 165 |
```{r image_file7, echo = F, results = 'asis'}
|
| 150 |
md_str <- paste('\n\n',sep="")
|
|
| 166 |
md_str <- paste('\n\n',sep="")
|
|
| 151 | 167 |
cat(md_str) |
| 152 | 168 |
``` |
| 153 | 169 |
|
| ... | ... | |
| 155 | 171 |
|
| 156 | 172 |
# Edge weigted mean slope: 20100901 |
| 157 | 173 |
|
| 174 |
|
|
| 158 | 175 |
```{r, echo=FALSE}
|
| 159 |
image_file8 <- file.path(in_dir,"Figure2_slope_mean_use_edge_weights_20100901_use_edge_weights_20100901_20100901_reg6_mosaic_run10_1500x4500_global_analyses_06212015.png") |
|
| 176 |
#fname <- paste("Figure2_slope_mean_use_edge_weights_20100901_use_edge_weights_20100901_20100901_",out_suffix,".png"#,sep="")
|
|
| 177 |
#fname <- "Figure2_slope_mean_use_edge_weights_20100901_use_edge_weights_20100901_20100901_reg4_mosaic_run10_1500x45#00_global_analyses_06212015.png" |
|
| 178 |
fname <- paste("Figure2_slope_mean_use_edge_weights_20100901_use_edge_weights_20100901_20100901_",out_suffix,".png",sep="")
|
|
| 179 |
image_file8 <- file.path(in_dir,fname) |
|
| 180 |
|
|
| 160 | 181 |
``` |
| 161 | 182 |
|
| 162 | 183 |
```{r image_file8, echo = F, results = 'asis'}
|
| 163 |
md_str <- paste('\n\n',sep="")
|
|
| 184 |
md_str <- paste('\n\n',sep="")
|
|
| 164 | 185 |
cat(md_str) |
| 165 | 186 |
``` |
| 166 | 187 |
|
| ... | ... | |
| 169 | 190 |
# Unweigted mean aspect: 20100831 |
| 170 | 191 |
|
| 171 | 192 |
```{r, echo=FALSE}
|
| 172 |
image_file9 <- file.path(in_dir,"Figure2_aspect_mean_unweighted_20100831_unweighted_20100831_20100831_reg6_mosaic_run10_1500x4500_global_analyses_06212015.png") |
|
| 193 |
|
|
| 194 |
fname <- paste("Figure2_aspect_mean_unweighted_20100831_unweighted_20100831_20100831_",out_suffix,".png",sep="")
|
|
| 195 |
image_file9 <- file.path(in_dir,fname) |
|
| 196 |
|
|
| 173 | 197 |
``` |
| 174 | 198 |
|
| 175 | 199 |
```{r image_file9, echo = F, results = 'asis'}
|
| 176 |
md_str <- paste('\n\n',sep="")
|
|
| 200 |
md_str <- paste('\n\n',sep="")
|
|
| 177 | 201 |
cat(md_str) |
| 178 | 202 |
``` |
| 179 | 203 |
|
| ... | ... | |
| 182 | 206 |
# Edge weigted mean aspect: 20100831 |
| 183 | 207 |
|
| 184 | 208 |
```{r, echo=FALSE}
|
| 185 |
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") |
|
| 209 |
#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") |
|
| 210 |
fname <- paste("Figure2_aspect_mean_use_edge_weights_20100831_use_edge_weights_20100831_20100831_",out_suffix,".png",sep="")
|
|
| 211 |
image_file10 <- file.path(in_dir,fname) |
|
| 212 |
|
|
| 186 | 213 |
``` |
| 187 | 214 |
|
| 188 | 215 |
```{r image_file10, echo = F, results = 'asis'}
|
| 189 |
md_str <- paste('\n\n',sep="")
|
|
| 216 |
md_str <- paste('\n\n',sep="")
|
|
| 190 | 217 |
cat(md_str) |
| 191 | 218 |
``` |
| 192 | 219 |
|
| ... | ... | |
| 195 | 222 |
# Unweigted mean aspect: 20100901 |
| 196 | 223 |
|
| 197 | 224 |
```{r, echo=FALSE}
|
| 198 |
image_file11 <- file.path(in_dir,"Figure2_aspect_mean_unweighted_20100901_unweighted_20100901_20100901_reg6_mosaic_run10_1500x4500_global_analyses_06212015.png") |
|
| 225 |
#image_file11 <- file.path(in_dir,"Figure2_aspect_mean_unweighted_20100901_unweighted_20100901_20100901_reg6_mosaic_#run10_1500x4500_global_analyses_06212015.png") |
|
| 226 |
fname <- paste("Figure2_aspect_mean_unweighted_20100901_unweighted_20100901_20100901_",out_suffix,".png",sep="")
|
|
| 227 |
image_file11 <- file.path(in_dir,fname) |
|
| 228 |
|
|
| 199 | 229 |
``` |
| 200 | 230 |
|
| 201 | 231 |
```{r image_file11, echo = F, results = 'asis'}
|
| 202 |
md_str <- paste('\n\n',sep="")
|
|
| 232 |
md_str <- paste('\n\n',sep="")
|
|
| 203 | 233 |
cat(md_str) |
| 204 | 234 |
``` |
| 205 | 235 |
|
| ... | ... | |
| 208 | 238 |
# Edge weigted mean aspect: 20100901 |
| 209 | 239 |
|
| 210 | 240 |
```{r, echo=FALSE}
|
| 211 |
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") |
|
| 241 |
#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") |
|
| 242 |
fname <- paste("Figure2_aspect_mean_use_edge_weights_20100901_use_edge_weights_20100901_20100901_",out_suffix,".png",sep="")
|
|
| 243 |
image_file12 <- file.path(in_dir,fname) |
|
| 212 | 244 |
``` |
| 213 | 245 |
|
| 214 | 246 |
```{r image_file12, echo = F, results = 'asis'}
|
| 215 |
md_str <- paste('\n\n',sep="")
|
|
| 247 |
md_str <- paste('\n\n',sep="")
|
|
| 216 | 248 |
cat(md_str) |
| 217 | 249 |
``` |
| 218 | 250 |
|
| ... | ... | |
| 221 | 253 |
# Unweigted - edge weighted diff: 20100831 |
| 222 | 254 |
|
| 223 | 255 |
```{r, echo=FALSE}
|
| 224 |
image_file13 <- file.path(in_dir,"Figure2_diff_raster_unweighted_use_edge_weights_20100831_reg6_mosaic_run10_1500x4500_global_analyses_06212015.png") |
|
| 256 |
#image_file13 <- file.path(in_dir,"Figure2_diff_raster_unweighted_use_edge_weights_20100831_reg6_mosaic_run10_1500x4#500_global_analyses_06212015.png") |
|
| 257 |
|
|
| 258 |
fname <- paste("Figure2_diff_raster_unweighted_use_edge_weights_20100831_",out_suffix,".png",sep="")
|
|
| 259 |
image_file13 <- file.path(in_dir,fname) |
|
| 225 | 260 |
``` |
| 226 | 261 |
|
| 227 | 262 |
```{r image_file13, echo = F, results = 'asis'}
|
| 228 |
md_str <- paste('\n\n',sep="")
|
|
| 263 |
md_str <- paste('\n\n',sep="")
|
|
| 229 | 264 |
cat(md_str) |
| 230 | 265 |
``` |
| 231 | 266 |
|
| ... | ... | |
| 234 | 269 |
# Unweigted - edge weighted diff : 20100901 |
| 235 | 270 |
|
| 236 | 271 |
```{r, echo=FALSE}
|
| 237 |
image_file14 <- file.path(in_dir,"Figure2_diff_raster_unweighted_use_edge_weights_20100901_reg6_mosaic_run10_1500x4500_global_analyses_06212015.png") |
|
| 272 |
#image_file14 <- file.path(in_dir,"Figure2_diff_raster_unweighted_use_edge_weights_20100901_reg6_mosaic_run10_1500x4#500_global_analyses_06212015.png") |
|
| 273 |
|
|
| 274 |
fname <- paste("Figure2_diff_raster_unweighted_use_edge_weights_20100901_",out_suffix,".png",sep="")
|
|
| 275 |
image_file14 <- file.path(in_dir,fname) |
|
| 276 |
|
|
| 238 | 277 |
``` |
| 239 | 278 |
|
| 240 | 279 |
```{r image_file14, echo = F, results = 'asis'}
|
| 241 |
md_str <- paste('\n\n',sep="")
|
|
| 280 |
md_str <- paste('\n\n',sep="")
|
|
| 242 | 281 |
cat(md_str) |
| 243 | 282 |
``` |
| 244 | 283 |
|
Also available in: Unified diff
automation of presentation for region 4 South America