ASCEE/lasp
ASCEE/lasp
5
1
Fork 0
Code Issues 5 Pull Requests Actions Releases 14 Activity

ignore last data point in smoothing, at Nyquist frequency

Browse Source
This commit is contained in:
Casper Jansen 2021-10-18 10:16:21 +02:00
parent cfd60f2826
commit 9aaa63512c
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
lasp/tools/tools.py
View File

@ -136,6 +136,7 @@ def smoothSpectralData(freq, M, sw: SmoothingWidth,
x0 = 1 if freq[0] == 0 else 0 # Skip first data point if zero frequency x0 = 1 if freq[0] == 0 else 0 # Skip first data point if zero frequency
Psm[0] = P[0] # Reuse old value in case first data.. Psm[0] = P[0] # Reuse old value in case first data..
# ..point is skipped. Not plotted any way. # ..point is skipped. Not plotted any way.
P[-1] = P[-2] # Last data point (Nyquist f) often invalid
# Loop through data points # Loop through data points
for x in range(x0, L): for x in range(x0, L):
@ -173,7 +174,7 @@ def smoothSpectralData(freq, M, sw: SmoothingWidth,
# Apply smoothing # Apply smoothing
Psm[x] = np.dot(g, P[xl:xu]) Psm[x] = np.dot(g, P[xl:xu])
if st == SmoothingType.levels: if st == SmoothingType.levels:
Psm = 10*np.log10(Psm) Psm = 10*np.log10(Psm)