Project

General

Profile

« Previous | Next » 

Revision 9b5125b9

Added by Adam Wilson over 11 years ago

Finished function to make raster extent for any modis tile

View differences:

climate/procedures/MOD35_ExtractProcessPath.r
42 42 
stitch="sudo MRTDATADIR=\"/usr/local/heg/2.12/data\" PGSHOME=/usr/local/heg/2.12/TOOLKIT_MTD PWD=/home/adamw /usr/local/heg/2.12/bin/swtif"
43 43 
stitch="/usr/local/heg/2.12/bin/swtif"
44 44 

  		




  45
  
  
    
stitch="sudo MRTDATADIR=\"/usr/local/heg/2.11/data\" PGSHOME=/usr/local/heg/2.11/TOOLKIT_MTD PWD=/home/adamw /usr/local/heg/2.11/bin/swtif"


  		




  46
  
  
    
files=paste(getwd(),"/",list.files("swath",pattern="hdf$",full=T),sep="")


  		




  
  45
  
    
#stitch="sudo MRTDATADIR=\"/usr/local/heg/2.11/data\" PGSHOME=/usr/local/heg/2.11/TOOLKIT_MTD PWD=/home/adamw /usr/local/heg/2.11/bin/swtif"


  		




  
  46
  
    
#files=paste(getwd(),"/",list.files("swath",pattern="hdf$",full=T),sep="")


  		




  47
  47
  
    

  		




  48
  48
  
    
## vars to process


  		




  49
  49
  
    
vars=as.data.frame(matrix(c(


  		




  ......




  141
  141
  
    
file.remove(gfiles[check==0])


  		




  142
  142
  
    

  		




  143
  143
  
    
## use new gdal


  		




  144
  
  
    
system(paste("/usr/local/gdal-1.10.0/bin/gdalwarp -wm 900 -overwrite -co COMPRESS=LZW -co PREDICTOR=2 -multi -r mode ",outdir,"/*.tif MOD35_path_gdalwarp.tif",sep=""))


  		




  
  144
  
    
system(paste("nohup /usr/local/gdal-1.10.0/bin/gdalwarp -wm 900 -overwrite -co COMPRESS=LZW -co PREDICTOR=2 -multi -r mode ",outdir,"/*.tif MOD35_path_gdalwarp.tif &",sep=""))


  		




  145
  145
  
    

  		




  146
  146
  
    

  		




  147
  147
  
    
###  Merge them into a geotiff


  		












Also available in:  Unified diff