| @@ -9,7 +9,6 import numpy | |||
|
|
9 | 9 | import datetime |
|
|
10 | 10 | |
|
|
11 | 11 | from .jroheaderIO import SystemHeader, RadarControllerHeader |
|
|
12 | # from schainpy import cSchain | |
|
|
13 | 12 | |
|
|
14 | 13 | |
|
|
15 | 14 | def getNumpyDtype(dataTypeCode): |
| @@ -63,46 +62,24 def hildebrand_sekhon(data, navg): | |||
|
|
63 | 62 | navg : numbers of averages |
|
|
64 | 63 | |
|
|
65 | 64 | Return: |
|
|
66 |
|
|
|
|
67 | anoise : noise's level | |
|
|
65 | mean : noise's level | |
|
|
68 | 66 | """ |
|
|
69 | 67 | |
|
|
70 |
sort |
|
|
|
71 | lenOfData = len(sortdata) | |
|
|
72 | nums_min = lenOfData*0.2 | |
|
|
73 | ||
|
|
74 | if nums_min <= 5: | |
|
|
75 | nums_min = 5 | |
|
|
76 | ||
|
|
77 | sump = 0. | |
|
|
78 | ||
|
|
79 | sumq = 0. | |
|
|
80 | ||
|
|
81 | j = 0 | |
|
|
82 | ||
|
|
83 | cont = 1 | |
|
|
84 | ||
|
|
85 | while((cont==1)and(j<lenOfData)): | |
|
|
86 | ||
|
|
87 | sump += sortdata[j] | |
|
|
88 | ||
|
|
89 | sumq += sortdata[j]**2 | |
|
|
90 | ||
|
|
91 | if j > nums_min: | |
|
|
92 | rtest = float(j)/(j-1) + 1.0/navg | |
|
|
93 | if ((sumq*j) > (rtest*sump**2)): | |
|
|
94 | j = j - 1 | |
|
|
95 | sump = sump - sortdata[j] | |
|
|
96 | sumq = sumq - sortdata[j]**2 | |
|
|
97 | cont = 0 | |
|
|
98 | ||
|
|
99 | j += 1 | |
|
|
100 | ||
|
|
101 | lnoise = sump /j | |
|
|
102 | ||
|
|
103 | return lnoise | |
|
|
68 | sorted_spectrum = numpy.sort(data, axis=None) | |
|
|
69 | nnoise = len(sorted_spectrum) # default to all points in the spectrum as noise | |
|
|
70 | for npts in range(1, len(sorted_spectrum)+1): | |
|
|
71 | partial = sorted_spectrum[:npts] | |
|
|
72 | mean = partial.mean() | |
|
|
73 | var = partial.var() | |
|
|
74 | if var * navg < mean**2.: | |
|
|
75 | nnoise = npts | |
|
|
76 | else: | |
|
|
77 | # partial spectrum no longer has characteristics of white noise | |
|
|
78 | break | |
|
|
104 | 79 | |
|
|
105 | # return cSchain.hildebrand_sekhon(sortdata, navg) | |
|
|
80 | noise_spectrum = sorted_spectrum[:nnoise] | |
|
|
81 | mean = noise_spectrum.mean() | |
|
|
82 | return mean | |
|
|
106 | 83 | |
|
|
107 | 84 | |
|
|
108 | 85 | class Beam: |
| @@ -231,6 +208,12 class JROData(GenericData): | |||
|
|
231 | 208 | |
|
|
232 | 209 | profileIndex = None |
|
|
233 | 210 | |
|
|
211 | error = (0, '') | |
|
|
212 | ||
|
|
213 | def __str__(self): | |
|
|
214 | ||
|
|
215 | return '{} - {}'.format(self.type, self.getDatatime()) | |
|
|
216 | ||
|
|
234 | 217 | def getNoise(self): |
|
|
235 | 218 | |
|
|
236 | 219 | raise NotImplementedError |
| @@ -1,8 +1,6 | |||
|
|
1 | 1 | import sys |
|
|
2 | 2 | import numpy |
|
|
3 | 3 | from scipy import interpolate |
|
|
4 | #TODO | |
|
|
5 | #from schainpy import cSchain | |
|
|
6 | 4 | from schainpy.model.proc.jroproc_base import ProcessingUnit, MPDecorator, Operation |
|
|
7 | 5 | from schainpy.model.data.jrodata import Voltage |
|
|
8 | 6 | from schainpy.utils import log |
| @@ -45,9 +45,6 setup(name = "schainpy", | |||
|
|
45 | 45 | }, |
|
|
46 | 46 | include_package_data = False, |
|
|
47 | 47 | scripts = ['schainpy/gui/schainGUI'], |
|
|
48 | ext_modules = [ | |
|
|
49 | Extension("cSchain", ["schainpy/model/proc/extensions.c"]) | |
|
|
50 | ], | |
|
|
51 | 48 | entry_points = { |
|
|
52 | 49 | 'console_scripts': [ |
|
|
53 | 50 | 'schain = schainpy.cli.cli:main', |
| @@ -59,9 +56,6 setup(name = "schainpy", | |||
|
|
59 | 56 | "scipy >= 0.14.0", |
|
|
60 | 57 | "h5py >= 2.2.1", |
|
|
61 | 58 | "matplotlib >= 2.0.0", |
|
|
62 | "pyfits >= 3.4", | |
|
|
63 | "paramiko >= 2.1.2", | |
|
|
64 | "paho-mqtt >= 1.2", | |
|
|
65 | 59 | "zmq", |
|
|
66 | 60 | "fuzzywuzzy", |
|
|
67 | 61 | "click", |
|
|
1 | NO CONTENT: file was removed |
|
|
1 | NO CONTENT: file was removed |
General Comments 0
You need to be logged in to leave comments.
Login now