1 """Static Spatial Variable
2
3
4 """
5 from GP_Util import *
6
7 -def Import(study,varname,inraster):
8 """ Import a static spatial variable
9 A static spatial variable is one that is constant through the entire
10 study period. (ie elevation). Its value will be associated with
11 corresponding study points though all sample dates
12
13 Inputs
14 varname name of source data and output variable name
15 inraster raster containing source data """
16
17 vardict = {
18 'type' : 'static',
19 'inraster' : inraster,
20 }
21
22 study.AddSource(varname,vardict)
23
25 log.info('Extract Static : '+varname.ljust(30))
26 vardict = study.infodict[varname]
27 inraster = vardict['inraster']
28
29 os.chdir(study.workingdir)
30 datadir = MakeFolder(os.path.abspath(varname))
31
32 outraster = os.path.abspath(varname)
33 extraster = outraster+'_r'
34 asciifile = outraster+'.txt'
35 npyfile = outraster+'.npy'
36
37
38 gp.ExtractByMask_sa(inraster, study.studyshp, extraster)
39
40
41 gp.Resample_management(extraster, outraster, study.cellsize)
42
43
44 gp.RasterToASCII_Conversion(outraster,asciifile)
45 data = num.loadtxt(asciifile,skiprows=6).flatten()
46
47
48
49
50 colvect = data.reshape((data.shape[0],1))
51 dtable = colvect * num.ones((1,study.enddates.shape[0]))
52
53
54 num.save(npyfile,dtable)
55 vardict = {'name':varname, 'file':npyfile, 'list':[] }
56 study.AddSeries(vardict)
57
58
59 os.remove(asciifile)
60 gp.delete(extraster)
61