schain Commit - r1264:a73eb4598037 · Jicamarca Repository

1

#!python

2

'''

3

'''

4

5

import os, sys

6

import datetime

7

import time

8

9

#path = os.path.dirname(os.getcwd())

10

#path = os.path.dirname(path)

11

#sys.path.insert(0, path)

12

13

from schainpy.controller import Project

14

15

desc = "USRP_test"

16

filename = "USRP_processing.xml"

17

controllerObj = Project()

18

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

19

20

############## USED TO PLOT IQ VOLTAGE, POWER AND SPECTRA #############

21

22

#######################################################################

23

######PATH DE LECTURA, ESCRITURA, GRAFICOS Y ENVIO WEB#################

24

#######################################################################

25

#path = '/media/data/data/vientos/57.2063km/echoes/NCO_Woodman'

26

27

28

path = '/home/soporte/data_hdf5' #### with clock 35.16 db noise

29

30

figpath = '/home/soporte/data_hdf5_imag'

31

#remotefolder = "/home/wmaster/graficos"

32

#######################################################################

33

################# RANGO DE PLOTEO######################################

34

#######################################################################

35

dBmin = '30'

36

dBmax = '60'

37

xmin = '0'

38

xmax ='24'

39

ymin = '0'

40

ymax = '600'

41

#######################################################################

42

########################FECHA##########################################

43

#######################################################################

44

str = datetime.date.today()

45

today = str.strftime("%Y/%m/%d")

46

str2 = str - datetime.timedelta(days=1)

47

yesterday = str2.strftime("%Y/%m/%d")

48

#######################################################################

49

######################## UNIDAD DE LECTURA#############################

50

#######################################################################

51

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

52

path=path,

53

startDate="2019/01/01",#today,

54

endDate="2109/12/30",#today,

55

startTime='00:00:00',

56

endTime='23:59:59',

57

delay=0,

58

#set=0,

59

online=0,

60

walk=1,

61

ippKm = 1000)

62

63

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

64

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

65

#######################################################################

66

################ OPERACIONES DOMINIO DEL TIEMPO########################

67

#######################################################################

68

69

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

70

#

71

# codigo64='1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1,1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,0,0,0,1,0,0,1,0,1,1,1,0,0,0,1,0,'+\

72

# '1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1,0,0,0,1,0,0,1,0,0,0,0,1,1,1,0,1,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1'

73

74

#opObj11 = procUnitConfObjA.addOperation(name='setRadarFrequency')

75

#opObj11.addParameter(name='frequency', value='30e6', format='float')

76

77

#opObj10 = procUnitConfObjA.addOperation(name='Scope', optype='external')

78

#opObj10.addParameter(name='id', value='10', format='int')

79

##opObj10.addParameter(name='xmin', value='0', format='int')

80

##opObj10.addParameter(name='xmax', value='50', format='int')

81

#opObj10.addParameter(name='type', value='iq')

82

#opObj10.addParameter(name='ymin', value='-5000', format='int')

83

##opObj10.addParameter(name='ymax', value='8500', format='int')

84

85

#opObj10 = procUnitConfObjA.addOperation(name='setH0')

86

#opObj10.addParameter(name='h0', value='-5000', format='float')

87

88

#opObj11 = procUnitConfObjA.addOperation(name='filterByHeights')

89

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

90

91

#codigo='1,1,-1,1,1,-1,1,-1,-1,1,-1,-1,-1,1,-1,-1,-1,1,-1,-1,-1,1,1,1,1,-1,-1,-1'

92

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

93

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

94

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

95

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

96

97

#opObj11 = procUnitConfObjA.addOperation(name='CohInt', optype='other')

98

#opObj11.addParameter(name='n', value='100', format='int')

99

100

#######################################################################

101

########## OPERACIONES DOMINIO DE LA FRECUENCIA########################

102

#######################################################################

103

procUnitConfObjSousySpectra = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjA.getId())

104

procUnitConfObjSousySpectra.addParameter(name='nFFTPoints', value='100', format='int')

105

procUnitConfObjSousySpectra.addParameter(name='nProfiles', value='100', format='int')

106

#procUnitConfObjSousySpectra.addParameter(name='pairsList', value='(0,0),(1,1),(0,1)', format='pairsList')

107

108

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

109

#opObj13.addParameter(name='mode', value='2', format='int')

110

111

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

112

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

113

#######################################################################

114

########## PLOTEO DOMINIO DE LA FRECUENCIA#############################

115

#######################################################################

116

#SpectraPlot

117

118

opObj11 = procUnitConfObjSousySpectra.addOperation(name='SpectraPlot', optype='external')

119

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

120

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

121

#opObj11.addParameter(name='xmin', value=-0.01, format='float')

122

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

123

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

124

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

125

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

126

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

127

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

128

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

129

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

130

131

132

#RTIPLOT

133

134

opObj11 = procUnitConfObjSousySpectra.addOperation(name='RTIPlot', optype='external')

135

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

136

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

137

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

138

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

139

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

140

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

141

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

142

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

143

144

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

145

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

146

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

147

148

149

# opObj11 = procUnitConfObjSousySpectra.addOperation(name='CrossSpectraPlot', optype='other')

150

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

151

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

152

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

153

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

154

# opObj11.addParameter(name='phase_cmap', value='jet', format='str')

155

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

156

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

157

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

158

# opObj11.addParameter(name='save', value=0, format='bool')

159

# opObj11.addParameter(name='pairsList', value='(0,1)', format='pairsList')

160

# #

161

# opObj11 = procUnitConfObjSousySpectra.addOperation(name='CoherenceMap', optype='other')

162

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

163

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

164

# opObj11.addParameter(name='phase_cmap', value='jet', format='str')

165

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

166

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

167

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

168

# opObj11.addParameter(name='save', value=0, format='bool')

169

# opObj11.addParameter(name='pairsList', value='(0,1)', format='pairsList')

170

#

171

#######################################################################

172

############### UNIDAD DE ESCRITURA ###################################

173

#######################################################################

174

#opObj11 = procUnitConfObjSousySpectra.addOperation(name='SpectraWriter', optype='other')

175

#opObj11.addParameter(name='path', value=wr_path)

176

#opObj11.addParameter(name='blocksPerFile', value='50', format='int')

177

print ("Escribiendo el archivo XML")

178

print ("Leyendo el archivo XML")

179

180

181

182

controllerObj.start()

183

	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,183
		1	#!python
		2	'''
		3	'''
		4
		5	import os, sys
		6	import datetime
		7	import time
		8
		9	#path = os.path.dirname(os.getcwd())
		10	#path = os.path.dirname(path)
		11	#sys.path.insert(0, path)
		12
		13	from schainpy.controller import Project
		14
		15	desc = "USRP_test"
		16	filename = "USRP_processing.xml"
		17	controllerObj = Project()
		18	controllerObj.setup(id = '191', name='Test_USRP', description=desc)
		19
		20	############## USED TO PLOT IQ VOLTAGE, POWER AND SPECTRA #############
		21
		22	#######################################################################
		23	######PATH DE LECTURA, ESCRITURA, GRAFICOS Y ENVIO WEB#################
		24	#######################################################################
		25	#path = '/media/data/data/vientos/57.2063km/echoes/NCO_Woodman'
		26
		27
		28	path = '/home/soporte/data_hdf5' #### with clock 35.16 db noise
		29
		30	figpath = '/home/soporte/data_hdf5_imag'
		31	#remotefolder = "/home/wmaster/graficos"
		32	#######################################################################
		33	################# RANGO DE PLOTEO######################################
		34	#######################################################################
		35	dBmin = '30'
		36	dBmax = '60'
		37	xmin = '0'
		38	xmax ='24'
		39	ymin = '0'
		40	ymax = '600'
		41	#######################################################################
		42	########################FECHA##########################################
		43	#######################################################################
		44	str = datetime.date.today()
		45	today = str.strftime("%Y/%m/%d")
		46	str2 = str - datetime.timedelta(days=1)
		47	yesterday = str2.strftime("%Y/%m/%d")
		48	#######################################################################
		49	######################## UNIDAD DE LECTURA#############################
		50	#######################################################################
		51	readUnitConfObj = controllerObj.addReadUnit(datatype='DigitalRFReader',
		52	path=path,
		53	startDate="2019/01/01",#today,
		54	endDate="2109/12/30",#today,
		55	startTime='00:00:00',
		56	endTime='23:59:59',
		57	delay=0,
		58	#set=0,
		59	online=0,
		60	walk=1,
		61	ippKm = 1000)
		62
		63	opObj11 = readUnitConfObj.addOperation(name='printInfo')
		64	opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
		65	#######################################################################
		66	################ OPERACIONES DOMINIO DEL TIEMPO########################
		67	#######################################################################
		68
		69	procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
		70	#
		71	# codigo64='1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1,1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,0,0,0,1,0,0,1,0,1,1,1,0,0,0,1,0,'+\
		72	# '1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1,0,0,0,1,0,0,1,0,0,0,0,1,1,1,0,1,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1'
		73
		74	#opObj11 = procUnitConfObjA.addOperation(name='setRadarFrequency')
		75	#opObj11.addParameter(name='frequency', value='30e6', format='float')
		76
		77	#opObj10 = procUnitConfObjA.addOperation(name='Scope', optype='external')
		78	#opObj10.addParameter(name='id', value='10', format='int')
		79	##opObj10.addParameter(name='xmin', value='0', format='int')
		80	##opObj10.addParameter(name='xmax', value='50', format='int')
		81	#opObj10.addParameter(name='type', value='iq')
		82	#opObj10.addParameter(name='ymin', value='-5000', format='int')
		83	##opObj10.addParameter(name='ymax', value='8500', format='int')
		84
		85	#opObj10 = procUnitConfObjA.addOperation(name='setH0')
		86	#opObj10.addParameter(name='h0', value='-5000', format='float')
		87
		88	#opObj11 = procUnitConfObjA.addOperation(name='filterByHeights')
		89	#opObj11.addParameter(name='window', value='1', format='int')
		90
		91	#codigo='1,1,-1,1,1,-1,1,-1,-1,1,-1,-1,-1,1,-1,-1,-1,1,-1,-1,-1,1,1,1,1,-1,-1,-1'
		92	#opObj11 = procUnitConfObjSousy.addOperation(name='Decoder', optype='other')
		93	#opObj11.addParameter(name='code', value=codigo, format='floatlist')
		94	#opObj11.addParameter(name='nCode', value='1', format='int')
		95	#opObj11.addParameter(name='nBaud', value='28', format='int')
		96
		97	#opObj11 = procUnitConfObjA.addOperation(name='CohInt', optype='other')
		98	#opObj11.addParameter(name='n', value='100', format='int')
		99
		100	#######################################################################
		101	########## OPERACIONES DOMINIO DE LA FRECUENCIA########################
		102	#######################################################################
		103	procUnitConfObjSousySpectra = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjA.getId())
		104	procUnitConfObjSousySpectra.addParameter(name='nFFTPoints', value='100', format='int')
		105	procUnitConfObjSousySpectra.addParameter(name='nProfiles', value='100', format='int')
		106	#procUnitConfObjSousySpectra.addParameter(name='pairsList', value='(0,0),(1,1),(0,1)', format='pairsList')
		107
		108	#opObj13 = procUnitConfObjSousySpectra.addOperation(name='removeDC')
		109	#opObj13.addParameter(name='mode', value='2', format='int')
		110
		111	#opObj11 = procUnitConfObjSousySpectra.addOperation(name='IncohInt', optype='other')
		112	#opObj11.addParameter(name='n', value='60', format='float')
		113	#######################################################################
		114	########## PLOTEO DOMINIO DE LA FRECUENCIA#############################
		115	#######################################################################
		116	#SpectraPlot
		117
		118	opObj11 = procUnitConfObjSousySpectra.addOperation(name='SpectraPlot', optype='external')
		119	opObj11.addParameter(name='id', value='1', format='int')
		120	opObj11.addParameter(name='wintitle', value='Spectra', format='str')
		121	#opObj11.addParameter(name='xmin', value=-0.01, format='float')
		122	#opObj11.addParameter(name='xmax', value=0.01, format='float')
		123	#opObj11.addParameter(name='zmin', value=dBmin, format='int')
		124	#opObj11.addParameter(name='zmax', value=dBmax, format='int')
		125	#opObj11.addParameter(name='ymin', value=ymin, format='int')
		126	#opObj11.addParameter(name='ymax', value=ymax, format='int')
		127	opObj11.addParameter(name='showprofile', value='1', format='int')
		128	opObj11.addParameter(name='save', value=figpath, format='str')
		129	opObj11.addParameter(name='save_period', value=10, format='int')
		130
		131
		132	#RTIPLOT
		133
		134	opObj11 = procUnitConfObjSousySpectra.addOperation(name='RTIPlot', optype='external')
		135	opObj11.addParameter(name='id', value='2', format='int')
		136	opObj11.addParameter(name='wintitle', value='RTIPlot', format='str')
		137	#opObj11.addParameter(name='zmin', value=dBmin, format='int')
		138	#opObj11.addParameter(name='zmax', value=dBmax, format='int')
		139	#opObj11.addParameter(name='ymin', value=ymin, format='int')
		140	#opObj11.addParameter(name='ymax', value=ymax, format='int')
		141	opObj11.addParameter(name='xmin', value=0, format='int')
		142	opObj11.addParameter(name='xmax', value=23, format='int')
		143
		144	opObj11.addParameter(name='showprofile', value='1', format='int')
		145	opObj11.addParameter(name='save', value=figpath, format='str')
		146	opObj11.addParameter(name='save_period', value=10, format='int')
		147
		148
		149	# opObj11 = procUnitConfObjSousySpectra.addOperation(name='CrossSpectraPlot', optype='other')
		150	# opObj11.addParameter(name='id', value='3', format='int')
		151	# opObj11.addParameter(name='wintitle', value='CrossSpectraPlot', format='str')
		152	# opObj11.addParameter(name='ymin', value=ymin, format='int')
		153	# opObj11.addParameter(name='ymax', value=ymax, format='int')
		154	# opObj11.addParameter(name='phase_cmap', value='jet', format='str')
		155	# opObj11.addParameter(name='zmin', value=dBmin, format='int')
		156	# opObj11.addParameter(name='zmax', value=dBmax, format='int')
		157	# opObj11.addParameter(name='figpath', value=figures_path, format='str')
		158	# opObj11.addParameter(name='save', value=0, format='bool')
		159	# opObj11.addParameter(name='pairsList', value='(0,1)', format='pairsList')
		160	# #
		161	# opObj11 = procUnitConfObjSousySpectra.addOperation(name='CoherenceMap', optype='other')
		162	# opObj11.addParameter(name='id', value='4', format='int')
		163	# opObj11.addParameter(name='wintitle', value='Coherence', format='str')
		164	# opObj11.addParameter(name='phase_cmap', value='jet', format='str')
		165	# opObj11.addParameter(name='xmin', value=xmin, format='float')
		166	# opObj11.addParameter(name='xmax', value=xmax, format='float')
		167	# opObj11.addParameter(name='figpath', value=figures_path, format='str')
		168	# opObj11.addParameter(name='save', value=0, format='bool')
		169	# opObj11.addParameter(name='pairsList', value='(0,1)', format='pairsList')
		170	#
		171	#######################################################################
		172	############### UNIDAD DE ESCRITURA ###################################
		173	#######################################################################
		174	#opObj11 = procUnitConfObjSousySpectra.addOperation(name='SpectraWriter', optype='other')
		175	#opObj11.addParameter(name='path', value=wr_path)
		176	#opObj11.addParameter(name='blocksPerFile', value='50', format='int')
		177	print ("Escribiendo el archivo XML")
		178	print ("Leyendo el archivo XML")
		179
		180
		181
		182	controllerObj.start()
		183