schain Commit - r620:60cf5c306c72 · Jicamarca Repository

1

import os, sys

1

import os, sys

2

#import timeit

2

#import timeit

3

import datetime

3

import datetime

4

5

path = os.path.split(os.getcwd())[0]

5

path = os.path.split(os.getcwd())[0]

6

sys.path.append(path)

6

sys.path.append(path)

7

8

from controller import *

8

from controller import *

9

dt1 = datetime.datetime.now()

9

dt1 = datetime.datetime.now()

10

desc = "MST-ISR-EEJ Experiment Test"

10

desc = "MST-ISR-EEJ Experiment Test"

11

filename = "isr_blocks.xml"

11

filename = "isr_blocks.xml"

12

13

controllerObj = Project()

13

controllerObj = Project()

14

15

controllerObj.setup(id = '191', name='test01', description=desc)

15

controllerObj.setup(id = '191', name='test01', description=desc)

16

17

path = '/media/signalchain/HD-PXU2/mst_isr_eej'

17

path = '/media/signalchain/HD-PXU2/mst_isr_eej'

18

path = '/media/data/DATA/MST_ISR_EEJ'

18

19

figpath = '/home/signalchain/Pictures/mst_isr_eej/isr'

20

figpath = '/home/signalchain/Pictures/mst_isr_eej/isr'

21

figpath = '/media/DATA/mst_isr_eej/isr'

22

20

23

21

readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader',

24

readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader',

22

path=path,

25

path=path,

23

startDate='201~~4/05~~/01',

26

startDate='2015/01/01',

24

endDate='201~~4/0~~5/30',

27

endDate='2015/12/30',

25

startTime='00:00:00',

28

startTime='00:00:00',

26

endTime='23:59:59',

29

endTime='23:59:59',

27

online=1,

30

online=1,

28

delay=10,

31

delay=10,

29

walk=0,

32

walk=1,

30

getblock=1)

33

getblock=1)

31

34

32

opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')

35

opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')

33

36

34

procUnitConfObjISR = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())

37

procUnitConfObjISR = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())

35

38

36

opObj11 = procUnitConfObjISR.addOperation(name='ProfileSelector', optype='other')

39

opObj11 = procUnitConfObjISR.addOperation(name='ProfileSelector', optype='other')

37

# profileIndex = '20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99'

40

# profileIndex = '20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99'

38

# opObj11.addParameter(name='profileList', value=profileIndex, format='intlist')

41

# opObj11.addParameter(name='profileList', value=profileIndex, format='intlist')

39

opObj11.addParameter(name='profileRangeList', value='20,99', format='intlist')

42

opObj11.addParameter(name='profileRangeList', value='20,99', format='intlist')

40

opObj11.addParameter(name='byblock', value='1', format='bool')

43

opObj11.addParameter(name='byblock', value='1', format='bool')

41

44

42

# opObj11 = procUnitConfObjISR.addOperation(name='ProfileConcat', optype='other')

45

# opObj11 = procUnitConfObjISR.addOperation(name='ProfileConcat', optype='other')

43

# opObj11.addParameter(name='m', value='5', format='int')

46

# opObj11.addParameter(name='m', value='5', format='int')

44

47

45

opObj11 = procUnitConfObjISR.addOperation(name='Reshaper', optype='other') #Esta Operacion opera sobre bloques y reemplaza el ProfileConcat que opera sobre perfiles

48

opObj11 = procUnitConfObjISR.addOperation(name='Reshaper', optype='other') #Esta Operacion opera sobre bloques y reemplaza el ProfileConcat que opera sobre perfiles

46

opObj11.addParameter(name='shape', value='4,16,6750', format='intlist') # shape = (nchannels, nprofiles, nhieghts)

49

opObj11.addParameter(name='shape', value='4,16,6750', format='intlist') # shape = (nchannels, nprofiles, nhieghts)

47

50

48

opObj11 = procUnitConfObjISR.addOperation(name='filterByHeights')

51

opObj11 = procUnitConfObjISR.addOperation(name='filterByHeights')

49

opObj11.addParameter(name='window', value='20', format='int')

52

opObj11.addParameter(name='window', value='20', format='int')

50

opObj11.addParameter(name='axis', value='2', format='int')

53

#opObj11.addParameter(name='axis', value='2', format='int')

51

54

52

barker3x1 = '1,1,-1,-1,-1,1'

55

barker3x1 = '1,1,-1,-1,-1,1'

53

#barker3x5 = '1,1,1,1,1, 1,1,1,1,1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,1,1,1,1,1'

56

#barker3x5 = '1,1,1,1,1, 1,1,1,1,1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,1,1,1,1,1'

54

57

55

opObj11 = procUnitConfObjISR.addOperation(name='Decoder', optype='other')

58

opObj11 = procUnitConfObjISR.addOperation(name='Decoder', optype='other')

56

opObj11.addParameter(name='code', value=barker3x1, format='floatlist')

59

opObj11.addParameter(name='code', value=barker3x1, format='floatlist')

57

opObj11.addParameter(name='nCode', value='2', format='int')

60

opObj11.addParameter(name='nCode', value='2', format='int')

58

#opObj11.addParameter(name='nBaud', value='15', format='int')

61

#opObj11.addParameter(name='nBaud', value='15', format='int')

59

opObj11.addParameter(name='nBaud', value='3', format='int')

62

opObj11.addParameter(name='nBaud', value='3', format='int')

60

opObj11.addParameter(name='mode', value='3', format='int')

63

opObj11.addParameter(name='mode', value='3', format='int')

61

opObj11.addParameter(name='times', value='8', format='int')

64

opObj11.addParameter(name='times', value='8', format='int')

62

opObj11.addParameter(name='osamp', value='5', format='int')

65

opObj11.addParameter(name='osamp', value='5', format='int')

63

66

64

67

65

procUnitConfObjISRSpectra = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjISR.getId())

68

procUnitConfObjISRSpectra = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjISR.getId())

66

procUnitConfObjISRSpectra.addParameter(name='nFFTPoints', value='16', format='int')

69

procUnitConfObjISRSpectra.addParameter(name='nFFTPoints', value='16', format='int')

67

procUnitConfObjISRSpectra.addParameter(name='nProfiles', value='16', format='int')

70

procUnitConfObjISRSpectra.addParameter(name='nProfiles', value='16', format='int')

68

71

69

opObj11 = procUnitConfObjISRSpectra.addOperation(name='IncohInt', optype='other')

72

opObj11 = procUnitConfObjISRSpectra.addOperation(name='IncohInt', optype='other')

70

opObj11.addParameter(name='n', value='36', format='float')

73

opObj11.addParameter(name='n', value='36', format='float')

71

74

72

opObj11 = procUnitConfObjISRSpectra.addOperation(name='SpectraPlot', optype='other')

75

opObj11 = procUnitConfObjISRSpectra.addOperation(name='SpectraPlot', optype='other')

73

opObj11.addParameter(name='id', value='300', format='int')

76

opObj11.addParameter(name='id', value='300', format='int')

77

opObj11.addParameter(name='zmin', value='23', format='int')

78

opObj11.addParameter(name='zmax', value='40', format='int')

74

opObj11.addParameter(name='wintitle', value='ISR', format='str')

79

opObj11.addParameter(name='wintitle', value='ISR', format='str')

75

opObj11.addParameter(name='figpath', value=figpath, format='str')

80

opObj11.addParameter(name='figpath', value=figpath, format='str')

76

opObj11.addParameter(name='wr_period', value='5', format='int')

81

opObj11.addParameter(name='wr_period', value='5', format='int')

77

opObj11.addParameter(name='exp_code', value='20', format='int')

82

opObj11.addParameter(name='exp_code', value='20', format='int')

78

83

79

opObj11 = procUnitConfObjISRSpectra.addOperation(name='RTIPlot', optype='other')

84

opObj11 = procUnitConfObjISRSpectra.addOperation(name='RTIPlot', optype='other')

80

opObj11.addParameter(name='id', value='301', format='int')

85

opObj11.addParameter(name='id', value='301', format='int')

86

opObj11.addParameter(name='xmin', value='00', format='int')

87

opObj11.addParameter(name='xmax', value='24', format='int')

88

opObj11.addParameter(name='zmin', value='23', format='int')

89

opObj11.addParameter(name='zmax', value='40', format='int')

81

opObj11.addParameter(name='wintitle', value='ISR', format='str')

90

opObj11.addParameter(name='wintitle', value='ISR', format='str')

82

opObj11.addParameter(name='showprofile', value='0', format='int')

91

opObj11.addParameter(name='showprofile', value='0', format='int')

83

opObj11.addParameter(name='figpath', value=figpath, format='str')

92

opObj11.addParameter(name='figpath', value=figpath, format='str')

84

opObj11.addParameter(name='wr_period', value='2', format='int')

93

opObj11.addParameter(name='wr_period', value='2', format='int')

85

opObj11.addParameter(name='exp_code', value='20', format='int')

94

opObj11.addParameter(name='exp_code', value='20', format='int')

86

95

87

96

88

opObj11 = procUnitConfObjISRSpectra.addOperation(name='SendByFTP', optype='other')

97

opObj11 = procUnitConfObjISRSpectra.addOperation(name='SendByFTP', optype='other')

89

opObj11.addParameter(name='ext', value='*.png', format='str')

98

opObj11.addParameter(name='ext', value='*.png', format='str')

90

opObj11.addParameter(name='localfolder', value=figpath, format='str')

99

opObj11.addParameter(name='localfolder', value=figpath, format='str')

91

opObj11.addParameter(name='remotefolder', value='/home/wmaster/graficos', format='str')

100

opObj11.addParameter(name='remotefolder', value='/home/wmaster/graficos', format='str')

92

opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')

101

opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')

93

opObj11.addParameter(name='username', value='wmaster', format='str')

102

opObj11.addParameter(name='username', value='wmaster', format='str')

94

opObj11.addParameter(name='password', value='mst2010vhf', format='str')

103

opObj11.addParameter(name='password', value='mst2010vhf', format='str')

95

opObj11.addParameter(name='period', value='5', format='int')

104

opObj11.addParameter(name='period', value='5', format='int')

96

105

97

106

98

107

99

print "Escribiendo el archivo XML"

108

print "Escribiendo el archivo XML"

100

controllerObj.writeXml(filename)

109

controllerObj.writeXml(filename)

101

print "Leyendo el archivo XML"

110

print "Leyendo el archivo XML"

102

controllerObj.readXml(filename)

111

controllerObj.readXml(filename)

103

112

104

controllerObj.createObjects()

113

controllerObj.createObjects()

105

controllerObj.connectObjects()

114

controllerObj.connectObjects()

106

115

107

#timeit.timeit('controllerObj.run()', number=2)

116

#timeit.timeit('controllerObj.run()', number=2)

108

117

109

controllerObj.run()

118

controllerObj.run()

110

#print fib(5)

119

#print fib(5)

111

120

112

dt2 = datetime.datetime.now()

121

dt2 = datetime.datetime.now()

113

print "======================="

122

print "======================="

114

print dt2-dt1

123

print dt2-dt1

115

print "======================="

124

print "======================="

	Site-wide shortcuts
/	Use quick search box
g h	Goto home page
g g	Goto my private gists page
g G	Goto my public gists page
g 0-9	Goto bookmarked items from 0-9
n r	New repository page
n g	New gist page

	Repositories
g s	Goto summary page
g c	Goto changelog page
g f	Goto files page
g F	Goto files page with file search activated
g p	Goto pull requests page
g o	Goto repository settings
g O	Goto repository access permissions settings
t s	Toggle sidebar on some pages

             import os, sys
             #import timeit
             import datetime
             path = os.path.split(os.getcwd())[0]
             sys.path.append(path)
             from controller import *
             dt1 = datetime.datetime.now()
             desc = "MST-ISR-EEJ Experiment Test"
             filename = "isr_blocks.xml"
             controllerObj = Project()
             controllerObj.setup(id = '191', name='test01', description=desc)
             path = '/media/signalchain/HD-PXU2/mst_isr_eej'
+            path = '/media/data/DATA/MST_ISR_EEJ'
             figpath = '/home/signalchain/Pictures/mst_isr_eej/isr'
+            figpath = '/media/DATA/mst_isr_eej/isr'
             readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader',
                                                         path=path,
-                                                        startDate='2014/05/01',
+                                                        startDate='2015/01/01',
-                                                        endDate='2014/05/30',
+                                                        endDate='2015/12/30',
                                                         startTime='00:00:00',
                                                         endTime='23:59:59',
                                                         online=1,
                                                         delay=10,
-                                                        walk=0,
+                                                        walk=1,
                                                         getblock=1)
             opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
             procUnitConfObjISR = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
             opObj11 = procUnitConfObjISR.addOperation(name='ProfileSelector', optype='other')
             # profileIndex =  '20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99'
             # opObj11.addParameter(name='profileList', value=profileIndex, format='intlist')
             opObj11.addParameter(name='profileRangeList', value='20,99', format='intlist')
             opObj11.addParameter(name='byblock', value='1', format='bool')
             # opObj11 = procUnitConfObjISR.addOperation(name='ProfileConcat', optype='other')
             # opObj11.addParameter(name='m', value='5', format='int')
             opObj11 = procUnitConfObjISR.addOperation(name='Reshaper', optype='other') #Esta Operacion opera sobre bloques y reemplaza el ProfileConcat que opera sobre perfiles
             opObj11.addParameter(name='shape', value='4,16,6750', format='intlist') # shape = (nchannels, nprofiles, nhieghts)
             opObj11 = procUnitConfObjISR.addOperation(name='filterByHeights')
             opObj11.addParameter(name='window', value='20', format='int')
-            opObj11.addParameter(name='axis', value='2', format='int')
+            #opObj11.addParameter(name='axis', value='2', format='int')
             barker3x1 = '1,1,-1,-1,-1,1'
             #barker3x5 = '1,1,1,1,1, 1,1,1,1,1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,1,1,1,1,1'
             opObj11 = procUnitConfObjISR.addOperation(name='Decoder', optype='other')
             opObj11.addParameter(name='code', value=barker3x1, format='floatlist')
             opObj11.addParameter(name='nCode', value='2', format='int')
             #opObj11.addParameter(name='nBaud', value='15', format='int')
             opObj11.addParameter(name='nBaud', value='3', format='int')
             opObj11.addParameter(name='mode', value='3', format='int')
             opObj11.addParameter(name='times', value='8', format='int')
             opObj11.addParameter(name='osamp', value='5', format='int')
             procUnitConfObjISRSpectra = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjISR.getId())
             procUnitConfObjISRSpectra.addParameter(name='nFFTPoints', value='16', format='int')
             procUnitConfObjISRSpectra.addParameter(name='nProfiles', value='16', format='int')
             opObj11 = procUnitConfObjISRSpectra.addOperation(name='IncohInt', optype='other')
             opObj11.addParameter(name='n', value='36', format='float')
             opObj11 = procUnitConfObjISRSpectra.addOperation(name='SpectraPlot', optype='other')
             opObj11.addParameter(name='id', value='300', format='int')
+            opObj11.addParameter(name='zmin', value='23', format='int')
+            opObj11.addParameter(name='zmax', value='40', format='int')
             opObj11.addParameter(name='wintitle', value='ISR', format='str')
             opObj11.addParameter(name='figpath', value=figpath, format='str')
             opObj11.addParameter(name='wr_period', value='5', format='int')
             opObj11.addParameter(name='exp_code', value='20', format='int')
             opObj11 = procUnitConfObjISRSpectra.addOperation(name='RTIPlot', optype='other')
             opObj11.addParameter(name='id', value='301', format='int')
+            opObj11.addParameter(name='xmin', value='00', format='int')
+            opObj11.addParameter(name='xmax', value='24', format='int')
+            opObj11.addParameter(name='zmin', value='23', format='int')
+            opObj11.addParameter(name='zmax', value='40', format='int')
             opObj11.addParameter(name='wintitle', value='ISR', format='str')
             opObj11.addParameter(name='showprofile', value='0', format='int')
             opObj11.addParameter(name='figpath', value=figpath, format='str')
             opObj11.addParameter(name='wr_period', value='2', format='int')
             opObj11.addParameter(name='exp_code', value='20', format='int')
             opObj11 = procUnitConfObjISRSpectra.addOperation(name='SendByFTP', optype='other')
             opObj11.addParameter(name='ext', value='*.png', format='str')
             opObj11.addParameter(name='localfolder', value=figpath, format='str')
             opObj11.addParameter(name='remotefolder', value='/home/wmaster/graficos', format='str')
             opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
             opObj11.addParameter(name='username', value='wmaster', format='str')
             opObj11.addParameter(name='password', value='mst2010vhf', format='str')
             opObj11.addParameter(name='period', value='5', format='int')
             print "Escribiendo el archivo XML"
             controllerObj.writeXml(filename)
             print "Leyendo el archivo XML"
             controllerObj.readXml(filename)
             controllerObj.createObjects()
             controllerObj.connectObjects()
             #timeit.timeit('controllerObj.run()', number=2)
             controllerObj.run()
             #print fib(5)
             dt2 = datetime.datetime.now()
             print "======================="
             print dt2-dt1
             print "======================="