Project

General

Profile

« Previous | Next » 

Revision 4d0fc33b

Added by Benoit Parmentier about 8 years ago

fixing generate raster function from tile

View differences:

climate/research/oregon/interpolation/global_product_assessment_part0_functions.R
204 204
  tile_coord <- names(list_spdf)[i]
205 205
  r_ref <- list_r_ref[[i]]
206 206
    
207
  tile_id <- paste0("tile_",i)
207 208
  df_tmp <- subset(df_missing,date==date_val,select=tile_coord)
208 209
  #for each row (date)
209 210
  val <- df_tmp[[tile_coord]]
......
216 217
  tile_spdf$predicted <- val
217 218
  tile_spdf$tile_coord <- tile_coord
218 219
  tile_spdf$overlap <- 1
220
  tile_spdf$tile_id <- tile_id
219 221
    
220 222
  #r <- rasterize(tile_spdf,r_ref,"predicted")
221 223
  #r <- rasterize(tile_spdf,r_ref,col_name)
......
231 233
   	r <- init(r, fun=set1f, overwrite=TRUE)
232 234
  }
233 235
    
234
  return(r)
236
  out_dir <- "." #can set this up later
237
  out_suffix_str <- paste0(col_name,"") # can set this parameter later
238
  raster_name <- file.path(out_dir,paste("r_",tile_id,"_",tile_coord,"_",out_suffix_str,".tif",sep=""))
239
  #raster_name <- 
240
  writeRaster(r, NAflag=NA_flag_val,filename=raster_name,overwrite=TRUE)  #unweighted mean
241

  
242
  return(raster_name)
235 243
}
236 244

  
237 245
############################ END OF SCRIPT ##################################

Also available in: Unified diff