Package firetools :: Module FocalSeries
[hide private]
[frames] | no frames]

Source Code for Module firetools.FocalSeries

 1  """Derived Data Series based on neighborhood statistics 
 2   
 3   
 4  """ 
 5  from GP_Util import *  
 6   
7 -def FocalSeries(study, varname, baseseries, shapeopt, statopt, dataopt='DATA'):
8 """Create a Focal series based on an existing series""" 9 os.chdir(study.workingdir) 10 workingdir = os.path.abspath(varname) 11 MakeFolder(workingdir) 12 gp.Workspace = workingdir 13 14 rasterlist = [os.path.splitext(f)[0] for f in study.seriesdict[baseseries]['list'] ] 15 16 focallist = [os.path.abspath(varname+str(i)) for (i,f) in enumerate(rasterlist)] 17 18 prog = ProgDots(len(rasterlist),title='FocalStat') 19 for (inraster,focraster) in zip(rasterlist,focallist): 20 gp.FocalStatistics_sa(inraster,focraster,shapeopt,statopt,dataopt) 21 prog.step() 22 23 focNpy = RasterToNpySeries(focallist) 24 25 joinfile = JoinDateRasters(focNpy,'SAGE') 26 27 vardict = {'name':varname,'file':joinfile,'list':focNpy} 28 29 study.AddSeries(vardict)
30