raytracing/plot_wspr.m at main · alexchartier/raytracing · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
%% load file
fname = 'data/wspr/2019-Mar-15-2200.nc';
links = read_netcdf(fname);

%% Plotting

txloc = links.txlocs;
rxloc = links.rxlocs;
home = links.home;

latI_bad = [];
lonI_bad = [];
latI_good = [];
lonI_good = [];
for l = 1:length(home)
    % Great circle calcs
    [latI,lonI] = interpm([txloc(l, 1), rxloc(l, 1)], ...
        [txloc(l, 2), rxloc(l, 2)], 1, 'gc');
    if home(l) == 0
        latI_bad = [latI_bad; latI; NaN];
        lonI_bad = [lonI_bad; lonI; NaN];
    else
        latI_good = [latI_good; latI; NaN];
        lonI_good = [lonI_good; lonI; NaN];
    end
end


earth_example
% Plot the links
CartL = sphcart([64E5 * ones(size(latI_bad)), deg2rad(latI_bad), deg2rad(lonI_bad)]);
h3 = plot3(CartL(:, 1), CartL(:, 2), CartL(:, 3), '-r', 'LineWidth', 3);


CartL = sphcart([64E5 * ones(size(latI_good)), deg2rad(latI_good), deg2rad(lonI_good)]);
h3 = plot3(CartL(:, 1), CartL(:, 2), CartL(:, 3), '-g', 'LineWidth', 3);
%CartL = sphcart([64E5 * ones(size(latI_f)), deg2rad(latI_f), deg2rad(lonI_f)]);


legend({'', 'link not predicted by SAMI', 'link predicted by SAMI'})
set(gca,'FontSize', 30)