Revision 1fe9e0ed
Added by Benoit Parmentier over 9 years ago
| climate/research/oregon/interpolation/global_run_scalingup_mosaicing_function.R | ||
|---|---|---|
| 403 | 403 |
raster_name <- file.path(out_dir_str,paste(raster_name_tmp,"_","m_",out_suffix,file_format,sep=""))#output file |
| 404 | 404 |
r_ref <- init(r_m, fun=set1f, filename=raster_name, overwrite=TRUE) |
| 405 | 405 |
#NAvalue(r_ref) <- -9999 |
| 406 |
|
|
| 407 |
cmd_str <- paste("/usr/bin/gdalwarp",inFilename,raster_name,sep=" ") #this may be a problem
|
|
| 406 |
python_cmd <- file.path(python_bin,"gdalwarp") |
|
| 407 |
|
|
| 408 |
#cmd_str <- paste("/usr/bin/gdalwarp",inFilename,raster_name,sep=" ") #this may be a problem
|
|
| 409 |
cmd_str <- paste(python_cmd,inFilename,raster_name,sep=" ") #this may be a problem |
|
| 410 |
|
|
| 411 |
|
|
| 412 |
|
|
| 413 |
|
|
| 414 |
|
|
| 408 | 415 |
#gdalwarp -t_srs '+proj=utm +zone=11 +datum=WGS84' raw_spot.tif utm11.tif |
| 409 | 416 |
|
| 410 | 417 |
system(cmd_str) |
| ... | ... | |
| 435 | 442 |
#9)out_dir: output directory, default is NULL |
| 436 | 443 |
# |
| 437 | 444 |
#OUTPUT: |
| 438 |
# Ojbec is produced with 3 components:
|
|
| 445 |
# Object is produced with 3 components:
|
|
| 439 | 446 |
# 1) mean_mosaic: list of raster files from produced mosaic , |
| 440 | 447 |
# 2) r_weights: list of raster files from weights |
| 441 | 448 |
# 3) r_weights_prod: list of raster files from product weights (weights*value) |
Also available in: Unified diff
more modifications to mosaicing function