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