Project

General

Profile

« Previous | Next » 

Revision 8f176858

Added by Benoit Parmentier over 11 years ago

LST processing function, test and modifications

View differences:

climate/research/oregon/interpolation/download_and_produce_MODIS_LST_climatology.R
111 111 
#list_tiles_modis <- c("h09v09,h10v09,h11v09,h12v09,h13v09,h14v09")
112 112 
#list_tiles_modis <- c("h09v09,h10v09,h11v09,h12v09,h13v09,h14v09,
113 113 
#                       h10v10,h11v10,h12v11,h13v12,h14v10")
114 
list_tiles_modis <-c("h30v10,h31v10,h32v10,h30v11,h31v11") #list("Queensland")
115 

  		




  114
  116
  
    
#list_tiles_modis <- c("h08v04")


  		




  115
  117
  
    
script_path<-"/home/parmentier/Data/IPLANT_project/env_layers_scripts/"


  		




  116
  118
  
    
modis_download_script <- file.path(script_path,"modis_download_05142013.py") # LST modis download python script


  		




  ......




  119
  121
  
    
#source(file.path(script_path,"download_and_produce_MODIS_LST_climatology_05162013.R"))


  		




  120
  122
  
    

  		




  121
  123
  
    
var="TMIN"


  		




  122
  
  
    
list_tiles_modis <- c("h08v04")


  		




  
  124
  
    
#list_tiles_modis <- c("h08v04")


  		




  123
  125
  
    
start_year = "2001"


  		




  124
  
  
    
#end_year = "2010"


  		




  125
  
  
    
end_year = "2002" #for testing (year included?)


  		




  
  126
  
    
end_year = "2010"


  		




  
  127
  
    
#end_year = "2002" #for testing (year included?)


  		




  126
  128
  
    
#end_month= "12"


  		




  127
  129
  
    
#start_month= "1"


  		




  128
  130
  
    
hdfdir =  '/home/layers/commons/modis/MOD11A1_tiles' #destination file where hdf files are stored locally after download.


  		




  129
  131
  
    
#hdfdir =  '/home/parmentier/Data/IPLANT_project/MOD11A1_tiles'


  		




  130
  132
  
    
download=1


  		




  131
  133
  
    
clim_calc=0


  		




  132
  
  
    
out_suffix_modis="_05132013"


  		




  
  134
  
    
out_suffix_modis="_05172013"


  		




  133
  135
  
    
#end_month= "12"


  		




  134
  136
  
    
#start_month= "1"


  		




  135
  137
  
    

  		




  136
  138
  
    
list_param_download_clim_LST_script <- list(list_tiles_modis,start_year,end_year,hdfdir,


  		




  137
  139
  
    
                                 var,grass_setting_script,modis_download_script, clim_script,


  		




  138
  
  
    
                                 night,download,clim_calc,out_suffix_modis)


  		




  
  140
  
    
                                 download,clim_calc,out_suffix_modis)


  		




  139
  141
  
    
names(list_param_download_clim_LST_script)<-c("list_tiles_modis","start_year","end_year","hdfdir",


  		




  140
  142
  
    
                                   "var","grass_setting_script","modis_download_script","clim_script",


  		




  141
  
  
    
                                   "night","download","clim_calc","out_suffix_modis")


  		




  
  143
  
    
                                   "download","clim_calc","out_suffix_modis")


  		




  142
  144
  
    
#debug(download_calculate_MODIS_LST_climatology)


  		




  143
  145
  
    

  		




  144
  
  
    
download_calculate_MODIS_LST_climatology(1,list_param_download_clim_LST_script)


  		




  145
  
  
    

  		




  
  146
  
    
#download_calculate_MODIS_LST_climatology(1,list_param_download_clim_LST_script)


  		




  146
  147
  
    

  		




  
  148
  
    
### END OF SCRIPT


  		












Also available in:  Unified diff