Project

General

Profile

« Previous | Next » 

Revision 3133718a

Added by Adam Wilson over 11 years ago

Added new script to process daily MOD35 files

View differences:

climate/procedures/MOD35_Explore.r
4 4 
library(raster)
5 5 
library(rgdal)
6 6 

  		




  7
  
  
    
f=list.files(pattern="Cloud_Mask_1_1")


  		




  
  7
  
    
f=list.files(pattern="*.hdf")


  		




  
  8
  
    

  		




  
  9
  
    
Sys.setenv(GEOL_AS_GCPS = "PARTIAL")


  		




  8
  10
  
    

  		




  9
  11
  
    
GDALinfo(f[1])


  		




  
  12
  
    
system(paste("gdalinfo",f[1]))


  		




  
  13
  
    
GDALinfo("HDF4_EOS:EOS_SWATH:\"MOD35_L2.A2000100.1445.006.2012252024758.hdf\":mod35:Cloud_Mask")


  		




  
  14
  
    
system("gdalinfo HDF4_EOS:EOS_SWATH:\"MOD35_L2.A2000100.1445.006.2012252024758.hdf\":mod35:Cloud_Mask | tail -n 200")


  		




  
  15
  
    

  		




  
  16
  
    
system("gdalwarp -overwrite -geoloc -order 2 -r near -s_srs \"EPSG:4326\" HDF4_EOS:EOS_SWATH:\"MOD35_L2.A2000100.1445.006.2012252024758.hdf\":mod35:Cloud_Mask cloudmask.tif")


  		




  
  17
  
    
system("gdalwarp -overwrite -r near -s_srs \"EPSG:4326\" HDF4_EOS:EOS_SWATH:\"MOD35_L2.A2000100.1445.006.2012252024758.hdf\":mod35:Cloud_Mask:1 cloudmask2.tif")


  		




  
  18
  
    

  		




  10
  19
  
    

  		




  11
  20
  
    
## get tile


  		




  12
  21
  
    
tile=raster("~/acrobates/projects/interp/data/modis/mod06/summary/MOD06_h09v04.nc",varname="CER")


  		












Also available in:  Unified diff