-
Notifications
You must be signed in to change notification settings - Fork 1
Expand file tree
/
Copy pathplot_multi_zones_reli_geoplot_rev4.m
More file actions
51 lines (45 loc) · 1.9 KB
/
plot_multi_zones_reli_geoplot_rev4.m
File metadata and controls
51 lines (45 loc) · 1.9 KB
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
42
43
44
45
46
47
48
49
50
51
function plot_multi_zones_reli_geoplot_rev4(app,temp_zones,base_protection_pts,data_label1,string_prop_model,temp_rel,grid_spacing)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Plot the multiple circles
%%%%1)Mitigation,
% %%2) Max knn dist,
% %%3)Convex Bound,
% %%4)Max Interference dB,
% %%%5)Prop Reliability
temp_zones=temp_zones(~cellfun('isempty', temp_zones(:,1)),:);
if ~isempty(temp_zones)
[num_cirlces,~]=size(temp_zones);
color_set3=plasma(num_cirlces);
f1=figure;
geoplot(base_protection_pts(:,1),base_protection_pts(:,2),'ob','LineWidth',3,'DisplayName','Federal System')
hold on;
if num_cirlces>1
for i=num_cirlces:-1:1
temp_bound=temp_zones{i,3};
if ~isempty(temp_bound)
temp_miti=temp_zones{i,1};
geoplot(temp_bound(:,1),temp_bound(:,2),'-','Color',color_set3(i,:),'LineWidth',3,'DisplayName',strcat(num2str(temp_miti),'dB'))
end
end
elseif num_cirlces==1 %%%%%%%%%%%%Probably don't have to graph if it's just one circle
temp_bound=temp_zones{3};
if ~isempty(temp_bound)
temp_miti=temp_zones{1,1};
geoplot(temp_bound(:,1),temp_bound(:,2),'-','Color',color_set3,'LineWidth',3,'DisplayName',strcat(num2str(temp_miti),'dB'))
end
end
grid on;
legend
%xlabel('Longitude')
%ylabel('Latitude')
%plot_google_map('maptype','terrain','APIKey','AIzaSyCgnWnM3NMYbWe7N4svoOXE7B2jwIv28F8') %%%Google's API key made by nick.matlab.error@gmail.com
pause(0.1)
%%%%geobasemap landcover
geobasemap streets-light%landcover
f1.Position = [100 100 1200 900];
pause(1)
filename1=strcat('Multi_Circles','_',data_label1,'_',string_prop_model,'_',num2str(temp_rel),'%_',num2str(grid_spacing),'km.png');
saveas(gcf,char(filename1))
pause(0.1);
close(f1)
end
end