Project

General

Profile

« Previous | Next » 

Revision 555815c9

Added by Adam Wilson about 11 years ago

Updated C5 MOD35 Processing Path code and C5MOD35_Evaluation in preparation for publication

View differences:

climate/procedures/MOD35C5_Evaluation.r
22 22 

  		




  23
  23
  
    
## mod35C6 annual


  		




  24
  24
  
    
if(!file.exists("data/MOD35C6_2009.tif")){


  		




  25
  
  
    
  system("/usr/local/gdal-1.10.0/bin/gdalbuildvrt  -a_srs '+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs' -sd 1 -b 1 data/MOD35C6.vrt `find /home/adamw/acrobates/adamw/projects/interp/data/modis/mod35/summary/ -name '*h[1]*_mean.nc'` ")


  		




  
  25
  
    
#  system("/usr/local/gdal-1.10.0/bin/gdalbuildvrt  -a_srs '+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs' -sd 1 -b 1 data/MOD35C6.vrt `find /home/adamw/acrobates/adamw/projects/interp/data/modis/mod35/summary/ -name '*h[1]*_mean.nc'` ")


  		




  
  26
  
    
  system("gdalbuildvrt data/MOD35C6.vrt `find /home/adamw/acrobates/adamw/projects/interp/data/modis/mod35/summary/ -name '*h[1]*_mean.nc'` ")


  		




  
  27
  
    

  		




  
  28
  
    
  system("/usr/local/gdal-1.10.0/bin/gdalbuildvrt -a_srs '+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs' -sd 4 -b 1 data/MOD35C6_CFday_pmiss.vrt `find /home/adamw/acrobates/adamw/projects/interp/data/modis/mod35/summary/ -name '*h[1]*.nc'` ")


  		




  
  29
  
    
  system("gdalwarp data/MOD35C6_CFday_pmiss.vrt data/MOD35C6_CFday_pmiss.tif -r bilinear")


  		




  
  30
  
    

  		




  26
  31
  
    
  system("align.sh data/MOD35C6.vrt data/MOD09_2009.tif data/MOD35C6_2009.tif")


  		




  
  32
  
    
  system("align.sh data/MOD35C6_CFday_pmiss.vrt data/MOD09_2009.tif data/MOD35C6_CFday_pmiss.tif")


  		




  27
  33
  
    
}


  		




  28
  34
  
    
mod35c6=raster("data/MOD35C6_2009_v1.tif")


  		




  29
  35
  
    
names(mod35c6)="C6MOD35CF"


  		












Also available in:  Unified diff