Project

General

Profile

« Previous | Next » 

Revision 905db3dc

Added by Adam Wilson over 11 years ago

Update MOD35 process path extraction to use sensor angle

View differences:

climate/procedures/MOD35_ExtractProcessPath.r
124 124 

  		




  125
  125
  
    

  		




  126
  126
  
    
### run it


  		




  127
  
  
    
mclapply(files[1:200],getpath,mc.cores=10)


  		




  
  127
  
    
mclapply(files,getpath,mc.cores=10)


  		




  128
  128
  
    

  		




  129
  129
  
    
## check gdal can read all of them


  		




  130
  130
  
    
gfiles=list.files(outdir,pattern="tif$",full=T)


  		




  ......




  145
  145
  
    

  		




  146
  146
  
    

  		




  147
  147
  
    
###  Merge them into a geotiff


  		




  148
  
  
    
    system(paste("gdal_merge.py -v -init 255 -n 255 -o MOD35_ProcessPath_gdalmerge2.tif -co \"ZLEVEL=9\" -co \"COMPRESS=LZW\" -co \"PREDICTOR=2\" `ls -d -1 gridded/*.tif --sort=size `",sep=""))


  		




  
  148
  
    
    system(paste("gdal_merge.py -v -init 255 -n 255 -o ",outdir,"/../MOD35_ProcessPath_gdalmerge2.tif -co \"ZLEVEL=9\" -co \"COMPRESS=LZW\" -co \"PREDICTOR=2\" `ls -d -1 ",outdir,"/*.tif --sort=size `",sep=""))


  		




  149
  149
  
    

  		




  150
  150
  
    
#  origin(raster(gfiles[5]))


  		




  151
  151
  
    
  


  		












Also available in:  Unified diff