##// END OF EJS Templates
schain-jc
RSS
Noise.py:...
Victor Sarmiento -
r93:1e5ff7f4bb15
parent child
Show More
Show file before | Show file after | Hide comments
@@ -1,7 +1,5
1 import numpy
1 import numpy
2
2 from Model.Spectra import Spectra
3 from Model.JROHeader import *
4 from Model.Voltage import Voltage
5
3
6 def hildebrand_sekhon(Data, navg=1 ):
4 def hildebrand_sekhon(Data, navg=1):
7 """
5 """
@@ -73,7 +71,7 class Noise():
73 m_DataObj = None
71 m_DataObj = None
74
72
75
73
76 def __init__(self, m_Voltage=None):
74 def __init__(self, m_Spectra=None):
77 """
75 """
78 Inicializador de la clase Noise para la la determinacion del nivel de ruido en un Spectro Doppler.
76 Inicializador de la clase Noise para la la determinacion del nivel de ruido en un Spectro Doppler.
79
77
@@ -83,13 +81,13 class Noise():
83 Return:
81 Return:
84 None
82 None
85 """
83 """
86 # if m_Voltage == None:
84 if m_Spectra == None:
87 # m_Voltage = Voltage()
85 m_Spectra = Spectra()
88 #
86
89 # if not(isinstance(m_Voltage, Voltage)):
87 if not(isinstance(m_Spectra, Spectra)):
90 # raise ValueError, "in Noise class, m_Voltage must be an Voltage class object"
88 raise ValueError, "in Noise class, m_Spectra must be an Spectra class object"
91
89
92 self.m_DataObj = m_Voltage
90 self.m_DataObj = m_Spectra
93
91
94
92
95 def getNoiseLevelByHildebrandSekhon(self):
93 def getNoiseLevelByHildebrandSekhon(self):
@@ -99,10 +97,16 class Noise():
99 Return:
97 Return:
100 noise level
98 noise level
101 """
99 """
102 data = self.m_DataObj.data # heights x perfiles
100 data = self.m_DataObj.data_spc
103 #heights = numpy.transpose( Data, (2,0,1) ) # channel x profile x height
101 daux = None
104 #data = Data[0,0,:]
102
103 for channel in range(self.m_DataObj.nChannels):
104 daux = data[channel,:,:]
105 noiselevel = hildebrand_sekhon(daux)
106 print noiselevel
105
107
106 noiselevel = hildebrand_sekhon(data)
107
108
109 for pair in range(self.m_DataObj.nPairs):
110 daux = data[pair,:,:]
111 noiselevel = hildebrand_sekhon(daux)
108 print noiselevel
112 print noiselevel
General Comments 0
You need to be logged in to leave comments. Login now