schain Commit - r1679:ec813df1dd48 · Jicamarca Repository

1

2

import os, sys

3

4

#from controller import *

5

from schainpy.controller import Project

6

7

desc = "EW DRIFTS Experiment"

8

filename = "EWDrifts.xml"

9

10

controllerObj = Project()

11

12

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

13

14

#Experimentos

15

16

path = "/home/soporte/minotauro/2023_06/Drifts_Faraday/main_radar/rawdata/"

17

#path = '/home/pcondor/Database/temp'

18

pathFigure = '/home/soporte/Documents/Database/ewdriftsschain2023'

19

pathFile = '/home/soporte/Documents/Database/ewdriftsschain2023'

20

21

xmin = 0

22

xmax = 24

23

#------------------------------------------------------------------------------------------------

24

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

25

path=path,

26

startDate='2023/06/23',

27

endDate='2023/06/23',

28

startTime='00:00:00',

29

endTime='23:59:59',

30

online=0,

31

walk=1)

32

33

#--------------------------------------------------------------------------------------------------

34

35

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

36

37

opObj11 = procUnitConfObj0.addOperation(name='selectChannels')

38

opObj11.addParameter(name='channelList', value='2,3,4,5', format='intlist')

39

40

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

41

opObj11.addParameter(name='profileRangeList', value='172,299', format='intlist')

42

43

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

44

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

45

46

code = [[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], [-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], [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], [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], [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], [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], [-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], [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], [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], [-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], [-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], [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], [-1, -1, 1], [1, 1, -1], [-1, -1, 1], [1, 1, -1], [-1, -1, 1], [-1, -1, 1], [-1, -1, 1], [1, 1, -1]]

47

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

48

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

49

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

50

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

51

52

opObj11 = procUnitConfObj0.addOperation(name='selectHeights')

53

opObj11.addParameter(name='minHei', value='0.0', format='float')

54

opObj11.addParameter(name='maxHei', value='960', format='float')

55

56

procUnitConfObj1 = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObj0.getId())

57

procUnitConfObj1.addParameter(name='nFFTPoints', value='128', format='int')

58

procUnitConfObj1.addParameter(name='nProfiles', value='128', format='int')

59

#procUnitConfObj1.addParameter(name='pairsList', value='(2,3),(4,5)', format='pairsList')#,(2,3)

60

procUnitConfObj1.addParameter(name='pairsList', value='(0,1),(2,3)', format='pairsList')

61

62

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

63

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

64

#opObj11.addParameter(name='timeInterval', value='300.0', format='float')

65

66

#opObj13 = procUnitConfObj1.addOperation(name='removeDC')

67

68

#opObj14 = procUnitConfObj1.addOperation(name='SpectraPlot', optype='other')

69

#opObj14.addParameter(name='id', value='65', format='int')

70

## # opObj14.addParameter(name='wintitle', value='Con interf', format='str')

71

#opObj14.addParameter(name='save', value=pathFigure, format='str')

72

##opObj14.addParameter(name='save_period', value=1, format='int')

73

#opObj14.addParameter(name='zmin', value='20', format='int')

74

#opObj14.addParameter(name='zmax', value='36', format='int')

75

#

76

77

#opObj12 = procUnitConfObj1.addOperation(name='RTIPlot', optype='other')

78

#opObj12.addParameter(name='id', value='63', format='int')

79

#opObj12.addParameter(name='wintitle', value='RTI Plot', format='str')

80

#opObj12.addParameter(name='save', value=pathFigure, format='str')

81

#opObj12.addParameter(name='save_period', value=10, format='int')

82

##opObj12.addParameter(name='figpath', value = pathFigure, format='str')

83

#opObj12.addParameter(name='xmin', value=xmin, format='float')

84

#opObj12.addParameter(name='xmax', value=xmax, format='float')

85

#opObj12.addParameter(name='zmin', value='20', format='int')

86

#opObj12.addParameter(name='zmax', value='36', format='int')

87

88

#--------------------------------------------------------------------------------------------------

89

90

procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObj1.getId())

91

opObj20 = procUnitConfObj2.addOperation(name='SpectralFitting', optype='other')

92

opObj20.addParameter(name='path', value='/home/soporte/UPDATE_SCHAIN/schain/schainpy/model/proc', format='str')

93

opObj20.addParameter(name='file', value='modelSpectralFitting', format='str')

94

opObj20.addParameter(name='groupList', value='(0,1),(2,3)',format='multiList')

95

#opObj20.addParameter(name='filec', value='weightfit', format='str')

96

97

#opObj12 = procUnitConfObj2.addOperation(name='HDF5Writer', optype='other')

98

#opObj12.addParameter(name='path', value=pathFile)

99

#opObj12.addParameter(name='blocksPerFile', value='3', format='int')

100

101

opObj21 = procUnitConfObj2.addOperation(name='EWDriftsEstimation', optype='other')

102

opObj21.addParameter(name='zenith', value='-1.64908, 3.36063', format='floatlist')

103

opObj21.addParameter(name='zenithCorrection', value='0.0', format='float')

104

opObj21.addParameter(name='fileDrifts', value=pathFile)

105

106

opObj24 = procUnitConfObj2.addOperation(name='SpectralMomentsPlot', optype='other')

107

opObj24.addParameter(name='id', value='1', format='int')

108

## # opObj14.addParameter(name='wintitle', value='Con interf', format='str')

109

opObj24.addParameter(name='save', value=pathFigure, format='str')

110

##opObj24.addParameter(name='save_period', value=1, format='int')

111

opObj24.addParameter(name='zmin', value='4', format='int')

112

opObj24.addParameter(name='zmax', value='20', format='int')

113

opObj24.addParameter(name='xaxis', value='Velocity', format='str')

114

#

115

titles=('SNR,Vertical Drifts,Zonal Drifts')

116

#titles=('Zonal Drifts,Vertical Drifts')

117

opObj23 = procUnitConfObj2.addOperation(name='GenericRTIPlot')

118

opObj23.addParameter(name='colormaps', value='jet,RdBu_r,RdBu_r')

119

opObj23.addParameter(name='attr_data', value='data_snr1,data_output')

120

#opObj23.addParameter(name='colormaps', value='RdBu,RdBu')

121

#opObj23.addParameter(name='attr_data', value='data_output')

122

opObj23.addParameter(name='wintitle', value='EW Drifts')

123

opObj23.addParameter(name='save', value=pathFigure)

124

opObj23.addParameter(name='titles', value=titles)

125

opObj23.addParameter(name='zfactors', value='1,1,1')

126

opObj23.addParameter(name='zlimits', value='(20,36),(-50,50),(-150,150)')

127

opObj23.addParameter(name='cb_labels', value='dB,m/s,m/s')

128

#opObj23.addParameter(name='titles', value=titles)

129

#opObj23.addParameter(name='zfactors', value='1,1')

130

#opObj23.addParameter(name='zlimits', value='(-150,150),(-40,40)')

131

#opObj23.addParameter(name='cb_labels', value='m/s,m/s')

132

opObj23.addParameter(name='throttle', value='1')

133

opObj23.addParameter(name='xmin', value=xmin)

134

opObj23.addParameter(name='xmax', value=xmax)

135

#--------------------------------------------------------------------------------------------------

136

137

controllerObj.start()

	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

			@@ -0,0 +1,137
		1
		2	import os, sys
		3
		4	#from controller import *
		5	from schainpy.controller import Project
		6
		7	desc = "EW DRIFTS Experiment"
		8	filename = "EWDrifts.xml"
		9
		10	controllerObj = Project()
		11
		12	controllerObj.setup(id = '191', name='test01', description=desc)
		13
		14	#Experimentos
		15
		16	path = "/home/soporte/minotauro/2023_06/Drifts_Faraday/main_radar/rawdata/"
		17	#path = '/home/pcondor/Database/temp'
		18	pathFigure = '/home/soporte/Documents/Database/ewdriftsschain2023'
		19	pathFile = '/home/soporte/Documents/Database/ewdriftsschain2023'
		20
		21	xmin = 0
		22	xmax = 24
		23	#------------------------------------------------------------------------------------------------
		24	readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader',
		25	path=path,
		26	startDate='2023/06/23',
		27	endDate='2023/06/23',
		28	startTime='00:00:00',
		29	endTime='23:59:59',
		30	online=0,
		31	walk=1)
		32
		33	#--------------------------------------------------------------------------------------------------
		34
		35	procUnitConfObj0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
		36
		37	opObj11 = procUnitConfObj0.addOperation(name='selectChannels')
		38	opObj11.addParameter(name='channelList', value='2,3,4,5', format='intlist')
		39
		40	opObj11 = procUnitConfObj0.addOperation(name='ProfileSelector', optype='other')
		41	opObj11.addParameter(name='profileRangeList', value='172,299', format='intlist')
		42
		43	opObj11 = procUnitConfObj0.addOperation(name='filterByHeights')
		44	opObj11.addParameter(name='window', value='3', format='int')
		45
		46	code = [[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], [-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], [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], [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], [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], [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], [-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], [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], [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], [-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], [-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], [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], [-1, -1, 1], [1, 1, -1], [-1, -1, 1], [1, 1, -1], [-1, -1, 1], [-1, -1, 1], [-1, -1, 1], [1, 1, -1]]
		47	opObj11 = procUnitConfObj0.addOperation(name='Decoder', optype='other')
		48	opObj11.addParameter(name='code', value=code, format='floatlist')
		49	opObj11.addParameter(name='nCode', value='128', format='int')
		50	opObj11.addParameter(name='nBaud', value='3', format='int')
		51
		52	opObj11 = procUnitConfObj0.addOperation(name='selectHeights')
		53	opObj11.addParameter(name='minHei', value='0.0', format='float')
		54	opObj11.addParameter(name='maxHei', value='960', format='float')
		55
		56	procUnitConfObj1 = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObj0.getId())
		57	procUnitConfObj1.addParameter(name='nFFTPoints', value='128', format='int')
		58	procUnitConfObj1.addParameter(name='nProfiles', value='128', format='int')
		59	#procUnitConfObj1.addParameter(name='pairsList', value='(2,3),(4,5)', format='pairsList')#,(2,3)
		60	procUnitConfObj1.addParameter(name='pairsList', value='(0,1),(2,3)', format='pairsList')
		61
		62	opObj11 = procUnitConfObj1.addOperation(name='IncohInt', optype='other')
		63	opObj11.addParameter(name='n', value='1', format='float')
		64	#opObj11.addParameter(name='timeInterval', value='300.0', format='float')
		65
		66	#opObj13 = procUnitConfObj1.addOperation(name='removeDC')
		67
		68	#opObj14 = procUnitConfObj1.addOperation(name='SpectraPlot', optype='other')
		69	#opObj14.addParameter(name='id', value='65', format='int')
		70	## # opObj14.addParameter(name='wintitle', value='Con interf', format='str')
		71	#opObj14.addParameter(name='save', value=pathFigure, format='str')
		72	##opObj14.addParameter(name='save_period', value=1, format='int')
		73	#opObj14.addParameter(name='zmin', value='20', format='int')
		74	#opObj14.addParameter(name='zmax', value='36', format='int')
		75	#
		76
		77	#opObj12 = procUnitConfObj1.addOperation(name='RTIPlot', optype='other')
		78	#opObj12.addParameter(name='id', value='63', format='int')
		79	#opObj12.addParameter(name='wintitle', value='RTI Plot', format='str')
		80	#opObj12.addParameter(name='save', value=pathFigure, format='str')
		81	#opObj12.addParameter(name='save_period', value=10, format='int')
		82	##opObj12.addParameter(name='figpath', value = pathFigure, format='str')
		83	#opObj12.addParameter(name='xmin', value=xmin, format='float')
		84	#opObj12.addParameter(name='xmax', value=xmax, format='float')
		85	#opObj12.addParameter(name='zmin', value='20', format='int')
		86	#opObj12.addParameter(name='zmax', value='36', format='int')
		87
		88	#--------------------------------------------------------------------------------------------------
		89
		90	procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObj1.getId())
		91	opObj20 = procUnitConfObj2.addOperation(name='SpectralFitting', optype='other')
		92	opObj20.addParameter(name='path', value='/home/soporte/UPDATE_SCHAIN/schain/schainpy/model/proc', format='str')
		93	opObj20.addParameter(name='file', value='modelSpectralFitting', format='str')
		94	opObj20.addParameter(name='groupList', value='(0,1),(2,3)',format='multiList')
		95	#opObj20.addParameter(name='filec', value='weightfit', format='str')
		96
		97	#opObj12 = procUnitConfObj2.addOperation(name='HDF5Writer', optype='other')
		98	#opObj12.addParameter(name='path', value=pathFile)
		99	#opObj12.addParameter(name='blocksPerFile', value='3', format='int')
		100
		101	opObj21 = procUnitConfObj2.addOperation(name='EWDriftsEstimation', optype='other')
		102	opObj21.addParameter(name='zenith', value='-1.64908, 3.36063', format='floatlist')
		103	opObj21.addParameter(name='zenithCorrection', value='0.0', format='float')
		104	opObj21.addParameter(name='fileDrifts', value=pathFile)
		105
		106	opObj24 = procUnitConfObj2.addOperation(name='SpectralMomentsPlot', optype='other')
		107	opObj24.addParameter(name='id', value='1', format='int')
		108	## # opObj14.addParameter(name='wintitle', value='Con interf', format='str')
		109	opObj24.addParameter(name='save', value=pathFigure, format='str')
		110	##opObj24.addParameter(name='save_period', value=1, format='int')
		111	opObj24.addParameter(name='zmin', value='4', format='int')
		112	opObj24.addParameter(name='zmax', value='20', format='int')
		113	opObj24.addParameter(name='xaxis', value='Velocity', format='str')
		114	#
		115	titles=('SNR,Vertical Drifts,Zonal Drifts')
		116	#titles=('Zonal Drifts,Vertical Drifts')
		117	opObj23 = procUnitConfObj2.addOperation(name='GenericRTIPlot')
		118	opObj23.addParameter(name='colormaps', value='jet,RdBu_r,RdBu_r')
		119	opObj23.addParameter(name='attr_data', value='data_snr1,data_output')
		120	#opObj23.addParameter(name='colormaps', value='RdBu,RdBu')
		121	#opObj23.addParameter(name='attr_data', value='data_output')
		122	opObj23.addParameter(name='wintitle', value='EW Drifts')
		123	opObj23.addParameter(name='save', value=pathFigure)
		124	opObj23.addParameter(name='titles', value=titles)
		125	opObj23.addParameter(name='zfactors', value='1,1,1')
		126	opObj23.addParameter(name='zlimits', value='(20,36),(-50,50),(-150,150)')
		127	opObj23.addParameter(name='cb_labels', value='dB,m/s,m/s')
		128	#opObj23.addParameter(name='titles', value=titles)
		129	#opObj23.addParameter(name='zfactors', value='1,1')
		130	#opObj23.addParameter(name='zlimits', value='(-150,150),(-40,40)')
		131	#opObj23.addParameter(name='cb_labels', value='m/s,m/s')
		132	opObj23.addParameter(name='throttle', value='1')
		133	opObj23.addParameter(name='xmin', value=xmin)
		134	opObj23.addParameter(name='xmax', value=xmax)
		135	#--------------------------------------------------------------------------------------------------
		136
		137	controllerObj.start()