@@ -86,8 +86,6 def hildebrand_sekhon(data, navg): | |||||
86 |
|
86 | |||
87 | sumq += sortdata[j]**2 |
|
87 | sumq += sortdata[j]**2 | |
88 |
|
88 | |||
89 | j += 1 |
|
|||
90 |
|
||||
91 | if j > nums_min: |
|
89 | if j > nums_min: | |
92 | rtest = float(j)/(j-1) + 1.0/navg |
|
90 | rtest = float(j)/(j-1) + 1.0/navg | |
93 | if ((sumq*j) > (rtest*sump**2)): |
|
91 | if ((sumq*j) > (rtest*sump**2)): | |
@@ -96,6 +94,8 def hildebrand_sekhon(data, navg): | |||||
96 | sumq = sumq - sortdata[j]**2 |
|
94 | sumq = sumq - sortdata[j]**2 | |
97 | cont = 0 |
|
95 | cont = 0 | |
98 |
|
96 | |||
|
97 | j += 1 | |||
|
98 | ||||
99 | lnoise = sump /j |
|
99 | lnoise = sump /j | |
100 | stdv = numpy.sqrt((sumq - lnoise**2)/(j - 1)) |
|
100 | stdv = numpy.sqrt((sumq - lnoise**2)/(j - 1)) | |
101 | return lnoise |
|
101 | return lnoise |
@@ -344,16 +344,15 class VoltageReader(JRODataReader, ProcessingUnit): | |||||
344 | If nTxs > 1 then one profile is divided by nTxs and number of total |
|
344 | If nTxs > 1 then one profile is divided by nTxs and number of total | |
345 | blocks is increased by nTxs (nProfiles *= nTxs) |
|
345 | blocks is increased by nTxs (nProfiles *= nTxs) | |
346 | """ |
|
346 | """ | |
347 | if self.nTxs == 1: |
|
|||
348 |
|
|
347 | self.dataOut.flagDataAsBlock = False | |
|
348 | ||||
|
349 | if self.nTxs == 1: | |||
349 | self.dataOut.data = self.datablock[:,self.profileIndex,:] |
|
350 | self.dataOut.data = self.datablock[:,self.profileIndex,:] | |
350 | self.dataOut.profileIndex = self.profileIndex |
|
351 | self.dataOut.profileIndex = self.profileIndex | |
351 |
|
352 | |||
352 | self.profileIndex += 1 |
|
353 | self.profileIndex += 1 | |
353 |
|
354 | |||
354 | else: |
|
355 | else: | |
355 | self.dataOut.flagDataAsBlock = False |
|
|||
356 |
|
||||
357 | iniHei_ForThisTx = (self.txIndex)*int(self.processingHeaderObj.nHeights/self.nTxs) |
|
356 | iniHei_ForThisTx = (self.txIndex)*int(self.processingHeaderObj.nHeights/self.nTxs) | |
358 | endHei_ForThisTx = (self.txIndex+1)*int(self.processingHeaderObj.nHeights/self.nTxs) |
|
357 | endHei_ForThisTx = (self.txIndex+1)*int(self.processingHeaderObj.nHeights/self.nTxs) | |
359 |
|
358 | |||
@@ -374,9 +373,9 class VoltageReader(JRODataReader, ProcessingUnit): | |||||
374 | """ |
|
373 | """ | |
375 | self.dataOut.flagDataAsBlock = True |
|
374 | self.dataOut.flagDataAsBlock = True | |
376 | self.dataOut.data = self.datablock |
|
375 | self.dataOut.data = self.datablock | |
377 |
self.dataOut.profileIndex = self.processingHeaderObj.profilesPerBlock |
|
376 | self.dataOut.profileIndex = self.processingHeaderObj.profilesPerBlock | |
378 |
|
377 | |||
379 |
self.profileIndex = self.processingHeaderObj.profilesPerBlock |
|
378 | self.profileIndex = self.processingHeaderObj.profilesPerBlock | |
380 |
|
379 | |||
381 | self.dataOut.flagNoData = False |
|
380 | self.dataOut.flagNoData = False | |
382 |
|
381 |
@@ -626,7 +626,7 class Decoder(Operation): | |||||
626 |
|
626 | |||
627 | return self.datadecTime |
|
627 | return self.datadecTime | |
628 |
|
628 | |||
629 | def run(self, dataOut, code=None, nCode=None, nBaud=None, mode = 0, osamp=None): |
|
629 | def run(self, dataOut, code=None, nCode=None, nBaud=None, mode = 0, osamp=None, times=None): | |
630 |
|
630 | |||
631 | if not self.isConfig: |
|
631 | if not self.isConfig: | |
632 |
|
632 |
@@ -5,7 +5,7 import datetime | |||||
5 | path = os.path.split(os.getcwd())[0] |
|
5 | path = os.path.split(os.getcwd())[0] | |
6 | path = os.path.split(path)[0] |
|
6 | path = os.path.split(path)[0] | |
7 |
|
7 | |||
8 |
sys.path. |
|
8 | sys.path.insert(0, path) | |
9 |
|
9 | |||
10 | from schainpy.controller import Project |
|
10 | from schainpy.controller import Project | |
11 |
|
11 |
General Comments 0
You need to be logged in to leave comments.
Login now