Project

General

Profile

« Previous | Next » 

Revision bd9bc508

Added by Jim Regetz over 12 years ago

  • ID bd9bc50860719d093dc4da94ae9b13e187f789a3

switched from truncating to rounding mean LST back to integer

View differences:

climate/extra/aggregate-daily-lst.py
146 146 
    gs.mapcalc('nobs_%s = %s' % (name, denominator), overwrite=overwrite)
147 147 
    numerator = '(%s)' % '+'.join(['if(isnull(%s), 0, %s)' % (m, m)
148 148 
        for m in maplist])
149 
    #gs.mapcalc('mean_%s = round(float(%s)/nobs_%s)' % (name, numerator, name),
150 
    gs.mapcalc('mean_%s = %s/nobs_%s' % (name, numerator, name),
149 
    gs.mapcalc('mean_%s = round(float(%s)/nobs_%s)' % (name, numerator, name),
151 150 
        overwrite=overwrite)
152 151 

  		




  153
  152
  
    
def load_qc_adjusted_lst(hdf):


  		




  ......




  208
  207
  
    

  		




  209
  208
  
    
# generate monthly pixelwise mean & count of high-quality daytime LST


  		




  210
  209
  
    
# values


  		




  211
  
  
    
### [atlas 17-May-2012] Wall time: 51.18 s


  		




  
  210
  
    
### [atlas 17-May-2012] Wall time: 53.79 s


  		




  212
  211
  
    
gs.os.environ['GRASS_OVERWRITE'] = '1'


  		




  213
  212
  
    
LST = [load_qc_adjusted_lst(hdf) for hdf in hdfs]


  		




  214
  213
  
    
calc_clim(LST, 'LST_%s_%d_%03d_%03d' % (tile, year, start_doy, end_doy))


  		












Also available in:  Unified diff