過去１週間分のADIFデータです。
最後の数時間のみ。
米国局のみ。
ADIFファイルをPSK Reporterからダウンロードして比較する（４）
QSOレートではなくて、レポーティングレートの比較です。グラフは、全ての局について。
これは、JA局以外。
これは、USA局のみ。
import re
import datetime
file_name1 = 'station1.adif'
def parse_adif(fn):
raw = re.split('<eor>|<eoh>',open(fn).read() )
raw.pop(0)
raw.pop()
logbook =[]
for record in raw:
qso = {}
tags = re.findall('<(.*?):\d+.*?>([^<\t\n\r\f\v]+)',record)
for tag in tags:
qso[tag[0]] = tag[1]
logbook.append(qso)
return logbook
def mydatetime(date, time):
dt = datetime.datetime(int(date[0:4]), int(date[4:6]), int(date[6:8]), \
int(time[0:2]), int(time[2:4]), int(time[4:6]))
return dt
l1 = parse_adif(file_name1)
count = 0
dt0 = datetime.datetime(2018, 5, 23, 11, 38, 44)
dt9 = datetime.datetime(2018, 5, 23, 15, 5, 29)
for qso1 in l1:
dt1 = mydatetime(qso1['QSO_DATE'], qso1['TIME_ON'])
if dt1 >= dt0 and dt1 <=dt9:
td = dt1 - dt0
print(td.total_seconds()/3600.0, count)
count += 1
% cat stationXYZ.adif | grep FREQ:8\>7\. | grep PSKREP | grep -v Japan > station1.adif
% python3 adif_parse.py > data1.csv
% R
> data1=read.csv("data1.csv", header=F, sep="")
> data2=read.csv("data2.csv", header=F, sep="")
> plot(data1$V1, data1$V2, xlab="time [hour]", ylab="no. of stations received", pch=16, col="#00ff0080", tck=1, xlim=c(0.0, 3.5), ylim=c(0,600))
> par(new=T)
> plot(data2$V1, data2$V2, xlab="time [hour]", ylab="no. of stations received", pch=16, col="#ff000080", tck=1, xlim=c(0.0, 3.5), ylim=c(0,600))
ADIFファイルをPSK Reporterからダウンロードして比較する（３）
Station 1:
from: QSO_DATE: 20180517, TIME_ON: 033544
to: QSO_DATE: 20180523, TIME_ON: 212959
5067 lines.
Station 2:
from: QSO_DATE: 20180523, TIME_ON: 104403
to: QSO_DATE: 20180524, TIME_ON: 025314
3504 lines.
First match: QSO_DATE: 20180523, TIME_ON: 113844
Last match: QSO_DATE: 20180523, TIME_ON: 212959
909 matches are found.
２つの局での”TIME_ON”時間の差です。多くの場合、同じDX局は同時にレポートされています。
時間の差は、秒単位です。
SNRの比較です。時間差がゼロの場合（緑）と、600秒以内の場合（赤）です。
file_name1 = 'station1.adif'
file_name2 = 'station2.adif'
def parse_adif(fn):
raw = re.split('<eor>|<eoh>',open(fn).read() )
raw.pop(0)
raw.pop()
logbook =[]
for record in raw:
qso = {}
tags = re.findall('<(.*?):\d+.*?>([^<\t\n\r\f\v]+)',record)
for tag in tags:
qso[tag[0]] = tag[1]
logbook.append(qso)
return logbook
def time2sec(time):
sec = int(time[0:2])*60*60+int(time[2:4])*60+int(time[4:6])
return sec
l1 = parse_adif(file_name1)
l2 = parse_adif(file_name2)
pattern = r"7."
for qso1 in l1:
if re.match(pattern, qso1['FREQ']):
for qso2 in l2:
if (qso1['CALL'] in qso2.values()):
if (qso1['QSO_DATE'] == qso2['QSO_DATE']):
time1 = time2sec(qso1['TIME_ON'])
time2 = time2sec(qso2['TIME_ON'])
if (abs(time1 - time2) < 600):
print(qso1['APP_PSKREP_SNR'], qso2['APP_PSKREP_SNR'], time1-time2)
% python3 adif_parse.py > data.csv
% R
> data=read.csv("data.csv", header=F, sep="")
> hist(data$V3, col="#00ffff80", breaks=seq(-625,625,50))
> plot(data$V3, col="#0000ff80", xlab="index", ylab="time diff")
> plot(data$V1, data$V2, xlab="My SNR", ylab="His SNR", pch=16, col="#00ff0080", tck=1, xlim=c(-30,20), ylim=c(-30,20))
ADIFファイルをPSK Reporterからダウンロードして比較する（２）
もし、同じDX局があなたとあなたのご近所とでほぼ同時に受信されれば、あなたは両者のSN値を比較してみたいと思うかもしれません。
import re
file_name1 = 'mystation.adif'
file_name2 = 'hisstation.adif'
def parse_adif(fn):
raw = re.split('<eor>|<eoh>',open(fn).read() )
raw.pop(0)
raw.pop()
logbook =[]
for record in raw:
qso = {}
tags = re.findall('<(.*?):(\d+).*?>([^<\t\n\r\f\v]+)',record)
for tag in tags:
qso[tag[0]] = tag[2]
logbook.append(qso)
return logbook
l1 = parse_adif(file_name1)
l2 = parse_adif(file_name2)
pattern = r"7."
for qso in l1:
if re.match(pattern, qso['FREQ']):
for hisqso in l2:
if (qso['CALL'] in hisqso.values()):
if (qso['QSO_DATE'] == hisqso['QSO_DATE']):
mytime = int(qso['TIME_ON'])
histime = int(hisqso['TIME_ON'])
if (abs(mytime - histime) < 2000):
print(qso['APP_PSKREP_SNR'], ', ', hisqso['APP_PSKREP_SNR'])
% python3 adif_parse.py > data.csv
% R
> data=read.csv("data.csv")
> plot(data, tck=1, col="red", xlim=c(-30,20), ylim=c(-30,20))
ADIFファイルをPSK Reporterからダウンロードして比較する
PSK reporterのマップ表示は、あなたの受信機のパフォーマンスをあなたのローカル局のそれと比較するのに有用です。
さらに、あなたは他の局からのレポートのADIFファイルをダウンロードすることも可能です。”Show logbook -> Download ADIF files -> last 24 hours, last week”を、クリックしてみて下さい。
import re
file_name1 = 'mystation.adif'
file_name2 = 'hisstation.adif'
def parse_adif(fn):
raw = re.split('<eor>|<eoh>',open(fn).read() )
raw.pop(0)
raw.pop()
logbook =[]
for record in raw:
qso = {}
tags = re.findall('<(.*?):(\d+).*?>([^<\t\n\r\f\v]+)',record)
for tag in tags:
qso[tag[0]] = tag[2]
logbook.append(qso)
return logbook
l1 = parse_adif(file_name1)
l2 = parse_adif(file_name2)
pattern = r"7." # 7MHz band only
for qso in l1:
if re.match(pattern, qso['FREQ']):
print ('+', qso['CALL'])
for qso in l2:
if re.match(pattern, qso['FREQ']):
print ('*', qso['CALL'])
この短いpythonプログラムは、２つのADIFファイル（典型的には、あなた自身のと、あなたのローカル局の）から情報を抽出します。
% python3 adif_parse.py > ttt.txt % grep + ttt.txt | sort | uniq -c | grep -e ' [KWN]' | sort -r 22 + WD6DBM 18 + KC7UG 8 + N7TM 8 + N5DG 7 + KW6S 7 + KN7D 7 + KE8FT 6 + N7FN 6 + KF7PG 6 + KA7T 5 + KF5ZTQ 5 + KE7UIU % grep '\*' ttt.txt | sort | uniq -c | grep -e ' [KWN]' | sort -r 15 * N6QQ 14 * KC7UG 13 * KN7D 12 * N7MDW 11 * K6KY 10 * N7TR 9 * N9BD 9 * N7IY 9 * K2ANT 8 * KE7W 7 * KD7H 7 * K6TE
数字は、各DX局がレポートされた回数を示しています。
SDR用のアンテナスイッチ（４）
これは、アンテナスイッチありの状態です。
ウォータフォールで、下の方に黒いエリアがあります。これは、RF出力電力は無しで、SEND信号のみをスイッチに加えています。
あとの２つのエリアは、それぞれ送信周波数7.020MHzでのセミ・ブレークイン動作と、フル・ブレークイン動作です。
SDR用のアンテナスイッチ（３）
スモールループアンテナを用いて7.027050MHzの信号を受信しながら、7.020MHzでダミーロードに送信しています。
特に気が付くような妨害は、認めれません。
今度は同じ周波数ですが、ダイポールアンテナで送信しています。利用可能な場所が制限されているので、２つのアンテナは非常に近接して配置されていることに注意して下さい。
SDR用のアンテナスイッチ（２）
これが、私のアンテナスイッチで用いられているリレーです。
ここに書かれている動作時間と復帰時間とが正しいか調べてみましょう。
IC-7410からのSEND信号が赤で、メーク接点に印加された直流電圧が黄色で示されています。
動作時間は約5msで、復帰時間は約4msです。
SDR用のアンテナスイッチ
私は、HF+ SDRにはスモールループアンテナを、IC-7410にはダイポールを使うつもりです。しかし、２つのアンテナは非常に近接しているので、私は、私のSDRを保護するためにアンテナスイッチを用いた方が良いだろうと思いました。
IC-7410からのSEND信号が、赤いトレースで表示され、7MHzのRF信号は黄色です。
もちろん、私はメカニカルリレーによる遅延もさらにチェックする必要があります。
HF+とSpyServer
HF+は、SpyServerを用いて、ネットワークを介して使うことができます。
[Ubuntu 18.04] % ./spyserver SPY Server v2.0.1629 - http://airspy.com Reading the configuration file: spyserver.config Listening for connections on 0.0.0.0:5555 Accepted client 192.168.0.108:65062 running SDR# v1.0.0.1666 on Microsoft Windows NT 6.2.9200.0 Device was sleeping. Wake up! Acquired the device Changed streaming mode => 1 Changed display pixels for FFT stream => 1000 Changed dB offset for FFT stream => 0 Changed dB range for FFT stream => 127 Changed digital gain for IQ stream => Automatic Changed format for FFT stream => 1 Changed format for IQ stream => 2 Changed frequency for IQ stream => 100000000 Hz Changed format for IQ stream => 2 Changed streaming state => 1 Changed frequency for IQ stream => 90000000 Hz Changed frequency for IQ stream => 80000000 Hz