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