Project

General

Profile

« Previous | Next » 

Revision 57d44dd9

Added by Adam Wilson about 11 years ago

Added routine to process identify regions of concern in modis cloud data

View differences:

climate/procedures/MOD35C6_Summary.r
2 2 

  		




  3
  3
  
    
setwd("~/acrobates/adamw/projects/MOD35C6")


  		




  4
  4
  
    

  		




  5
  
  
    
library(raster);beginCluster(10)


  		




  
  5
  
    
library(raster)#;beginCluster(10)


  		




  6
  6
  
    
library(rasterVis)


  		




  7
  7
  
    
library(rgdal)


  		




  8
  8
  
    
library(plotKML)


  		




  ......




  18
  18
  
    

  		




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


  		




  20
  20
  
    
  system("/usr/local/bin/pkcreatect -min 0 -max 100 -g -i data/MOD35C6_2009.tif -o data/MOD35C6_2009a.tif -ct none -co COMPRESS=LZW")


  		




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


  		




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


  		




  22
  22
  
    
}


  		




  23
  23
  
    
mod35c6=raster("data/MOD35C6_2009_v1.tif")


  		




  24
  24
  
    
names(mod35c6)="C6MOD35CF"


  		












Also available in:  Unified diff