Python read_DBYの例

コード例 #1

ファイル: skywave_measurements.py プロジェクト: junjunhencool/Skywaves

# -*- coding: utf-8 -*-
"""
Created on Fri Mar 18 14:20:29 2016

@author: lenz
"""
from read_DBY import read_DBY
from remove_60Hz_slope import remove_60Hz_slope
from chop_gw_ir import chop_gw_ir
import numpy as np
from lowpass import lowpass
import matplotlib.pyplot as plt
x_max = 450e-6

read_DBY_output = read_DBY(38, 1, 26.522908895, 8, x_max, 1)  #UF 15-38, RS1
#read_DBY_output=read_DBY(39,1,66.583436840,9,x_max,1) #UF 15-39, RS1
#read_DBY_output=read_DBY(40,3,20.767465080,10,x_max,1) #UF 15-40, RS3
#read_DBY_output=read_DBY(41,1,57.298446790,11,x_max,1) #UF 15-41, RS1
#read_DBY_output=read_DBY(42,4,43.058185590,12,x_max,1) #UF 15-42, RS4
#read_DBY_output=read_DBY(43,4,23.293418545,13,x_max,1) #UF 15-43, RS4
remove_60Hz_slope_output = remove_60Hz_slope(read_DBY_output, x_max)
print('UTC reference time: %s' % remove_60Hz_slope_output[2])
gw_time, gw_data, ir_time, ir_data = chop_gw_ir(remove_60Hz_slope_output)

gw_time = gw_time
gw = gw_data
ir_time = ir_time
ir = ir_data

gw_peak = gw_time[np.argmax(gw)]
print(gw_peak)

コード例 #2

ファイル: skywave_measurements.py プロジェクト: felipelenz/Skywaves

# -*- coding: utf-8 -*-
"""
Created on Fri Mar 18 14:20:29 2016

@author: lenz
"""
from read_DBY import read_DBY
from remove_60Hz_slope import remove_60Hz_slope
from chop_gw_ir import chop_gw_ir
import numpy as np
from lowpass import lowpass
import matplotlib.pyplot as plt
x_max=450e-6

read_DBY_output=read_DBY(38,1,26.522908895,8,x_max,1) #UF 15-38, RS1
#read_DBY_output=read_DBY(39,1,66.583436840,9,x_max,1) #UF 15-39, RS1
#read_DBY_output=read_DBY(40,3,20.767465080,10,x_max,1) #UF 15-40, RS3
#read_DBY_output=read_DBY(41,1,57.298446790,11,x_max,1) #UF 15-41, RS1
#read_DBY_output=read_DBY(42,4,43.058185590,12,x_max,1) #UF 15-42, RS4
#read_DBY_output=read_DBY(43,4,23.293418545,13,x_max,1) #UF 15-43, RS4
remove_60Hz_slope_output=remove_60Hz_slope(read_DBY_output,x_max)
print('UTC reference time: %s' %remove_60Hz_slope_output[2])
gw_time,gw_data,ir_time,ir_data=chop_gw_ir(remove_60Hz_slope_output)

gw_time=gw_time
gw=gw_data
ir_time=ir_time
ir=ir_data

gw_peak=gw_time[np.argmax(gw)]
print(gw_peak)

コード例 #3

ファイル: Waves_fft.py プロジェクト: junjunhencool/Skywaves

#read_DBY_output=read_DBY(41,2,57.373669615,11,x_max,10) #UF 15-41, RS2
#read_DBY_output=read_DBY(41,3,57.405116910,11,x_max,1) #UF 15-41, RS3
#read_DBY_output=read_DBY(41,4,57.555913445,11,x_max,10) #UF 15-41, RS4
#read_DBY_output=read_DBY(41,5,57.575066540,11,x_max,10) #UF 15-41, RS5
#read_DBY_output=read_DBY(42,1,42.712899355,12,x_max,10) #UF 15-42, RS1
#read_DBY_output=read_DBY(42,3,42.862766400,12,x_max,10) #UF 15-42, RS3
#read_DBY_output=read_DBY(42,4,43.058185590,12,x_max,10) #UF 15-42, RS4
#read_DBY_output=read_DBY(42,5,43.338093110,12,x_max,10) #UF 15-42, RS5
#read_DBY_output=read_DBY(42,6,43.366312590,12,x_max,10) #UF 15-42, RS6
#read_DBY_output=read_DBY(43,1,22.706011205,13,x_max,10) #UF 15-43, RS1
#read_DBY_output=read_DBY(43,2,22.934697575,13,x_max,10) #UF 15-43, RS2
#read_DBY_output=read_DBY(43,3,23.157666725,13,x_max,10) #UF 15-43, RS3
#read_DBY_output=read_DBY(43,4,23.293418545,13,x_max,10) #UF 15-43, RS4
#read_DBY_output=read_DBY(43,5,23.389957810,13,x_max,10) #UF 15-43, RS5

read_DBY_output = read_DBY(38, 1, 26.522908895, 8, x_max, 1)  #UF 15-38, RS1
remove_60Hz_slope_output = remove_60Hz_slope(read_DBY_output, x_max)
print('UTC reference time: %s' % remove_60Hz_slope_output[2])
gw_time, gw_data, ir_time, ir_data = chop_gw_ir(remove_60Hz_slope_output)
gw_time = gw_time
gw = gw_data
ir_time = ir_time
ir = ir_data
gw_peak = gw_time[np.argmax(gw)]
lpf_ir_time = ir_time
lpf_ir_data = lowpass(ir, 100e3)
#plt.plot(ir_time-gw_peak,ir/np.max(lpf_ir_data),label='raw')#,color='b')
plt.plot(lpf_ir_time - gw_peak, lpf_ir_data / np.max(lpf_ir_data),
         linewidth=2)  #,color='r')
lpf_gw_time = gw_time
lpf_gw_data = lowpass(gw, 900e3)

コード例 #4

#read_DBY_output=read_DBY(39,2,66.586552840,9,x_max,1) #UF 15-39, RS2
#read_DBY_output=read_DBY(39,3,66.633136315,9,x_max,1) #UF 15-39, RS3
#read_DBY_output=read_DBY(39,5,66.671758785,9,x_max,1) #UF 15-39, RS5
#read_DBY_output=read_DBY(40,2,20.746971600,10,x_max,1) #UF 15-40, RS2
#read_DBY_output=read_DBY(40,3,20.767465080,10,x_max,1) #UF 15-40, RS3
#read_DBY_output=read_DBY(41,1,57.298446790,11,x_max,1) #UF 15-41, RS1
#read_DBY_output=read_DBY(41,2,57.373669615,11,x_max,1) #UF 15-41, RS2
#read_DBY_output=read_DBY(41,3,57.405116910,11,x_max,1) #UF 15-41, RS3
#read_DBY_output=read_DBY(41,4,57.555913445,11,x_max,1) #UF 15-41, RS4
#read_DBY_output=read_DBY(41,5,57.575066540,11,x_max,1) #UF 15-41, RS5
#read_DBY_output=read_DBY(42,1,42.712899355,12,x_max,1) #UF 15-42, RS1
#read_DBY_output=read_DBY(42,3,42.862766400,12,x_max,1) #UF 15-42, RS3
#read_DBY_output=read_DBY(42,4,43.058185590,12,x_max,1) #UF 15-42, RS4
#read_DBY_output=read_DBY(42,5,43.338093110,12,x_max,1) #UF 15-42, RS5
#read_DBY_output=read_DBY(42,6,43.366312590,12,x_max,1) #UF 15-42, RS6
read_DBY_output=read_DBY(43,1,22.706011205,13,x_max,1) #UF 15-43, RS1
#read_DBY_output=read_DBY(43,2,22.934697575,13,x_max,1) #UF 15-43, RS2
#read_DBY_output=read_DBY(43,3,23.157666725,13,x_max,1) #UF 15-43, RS3
#read_DBY_output=read_DBY(43,4,23.293418545,13,x_max,1) #UF 15-43, RS4
#read_DBY_output=read_DBY(43,5,23.389957810,13,x_max,1) #UF 15-43, RS5
remove_60Hz_slope_output=remove_60Hz_slope(read_DBY_output,x_max)
print('UTC reference time: %s' %remove_60Hz_slope_output[2])
gw_time,gw_data,ir_time,ir_data=chop_gw_ir(remove_60Hz_slope_output)

##Uncomment for moving average filter
#n=len(remove_60Hz_slope_output[0])
#gw_time=remove_60Hz_slope_output[0][:n/2]
#gw=remove_60Hz_slope_output[1][:n/2]
#ir_time=remove_60Hz_slope_output[0][n/2:]
#ir=remove_60Hz_slope_output[1][n/2:]
#plt.plot(gw_time*1e6,gw,ir_time*1e6,ir)

コード例 #5

ファイル: Waves_fft.py プロジェクト: felipelenz/Skywaves

#read_DBY_output=read_DBY(41,4,57.555913445,11,x_max,10) #UF 15-41, RS4
#read_DBY_output=read_DBY(41,5,57.575066540,11,x_max,10) #UF 15-41, RS5
#read_DBY_output=read_DBY(42,1,42.712899355,12,x_max,10) #UF 15-42, RS1
#read_DBY_output=read_DBY(42,3,42.862766400,12,x_max,10) #UF 15-42, RS3
#read_DBY_output=read_DBY(42,4,43.058185590,12,x_max,10) #UF 15-42, RS4
#read_DBY_output=read_DBY(42,5,43.338093110,12,x_max,10) #UF 15-42, RS5
#read_DBY_output=read_DBY(42,6,43.366312590,12,x_max,10) #UF 15-42, RS6
#read_DBY_output=read_DBY(43,1,22.706011205,13,x_max,10) #UF 15-43, RS1
#read_DBY_output=read_DBY(43,2,22.934697575,13,x_max,10) #UF 15-43, RS2
#read_DBY_output=read_DBY(43,3,23.157666725,13,x_max,10) #UF 15-43, RS3
#read_DBY_output=read_DBY(43,4,23.293418545,13,x_max,10) #UF 15-43, RS4
#read_DBY_output=read_DBY(43,5,23.389957810,13,x_max,10) #UF 15-43, RS5



read_DBY_output=read_DBY(38,1,26.522908895,8,x_max,1) #UF 15-38, RS1
remove_60Hz_slope_output=remove_60Hz_slope(read_DBY_output,x_max)
print('UTC reference time: %s' %remove_60Hz_slope_output[2])
gw_time,gw_data,ir_time,ir_data=chop_gw_ir(remove_60Hz_slope_output)
gw_time=gw_time
gw=gw_data
ir_time=ir_time
ir=ir_data
gw_peak=gw_time[np.argmax(gw)]
lpf_ir_time=ir_time
lpf_ir_data=lowpass(ir,100e3)
#plt.plot(ir_time-gw_peak,ir/np.max(lpf_ir_data),label='raw')#,color='b')
plt.plot(lpf_ir_time-gw_peak,lpf_ir_data/np.max(lpf_ir_data),linewidth=2)#,color='r')
lpf_gw_time=gw_time
lpf_gw_data=lowpass(gw,900e3)
#plt.plot(gw_time-gw_peak,gw/np.max(lpf_ir_data),label='raw')#,color='b')