Generic_Bugsplat/plot_nationwide_single_overlap_rev1.m at main · nicklasorte/Generic_Bugsplat · 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
function []=plot_nationwide_single_overlap_rev1(app,sim_number,cell_poly_merge,cell_overlap_poly)


tic;
f1=figure;
%AxesH = axes;
for bound_idx=1%:1:num_miti
    temp_merge=cell_poly_merge{bound_idx};
    temp_merge
    [temp_geo_poly]=convert_polyshape2geopolyshape(app,temp_merge);
    %plot(temp_merge,'FaceColor','b','EdgeColor','k')
    geoplot(temp_geo_poly,'FaceColor','b','EdgeColor','k')
    hold on;

    temp_overlap_poly=cell_overlap_poly{bound_idx};
    temp_overlap_poly
    if ~isempty(temp_overlap_poly)
        if temp_overlap_poly.NumRegions==0
        else
            [temp_geo_overlap]=convert_polyshape2geopolyshape(app,temp_overlap_poly);
            %plot(temp_overlap_poly,'FaceColor','r','EdgeColor','k','FaceAlpha',0.7)
            geoplot(temp_geo_overlap,'FaceColor','r','EdgeColor','k','FaceAlpha',0.7)
        end
    end
end
grid on;
pause(0.1)

%%%%geobasemap landcover
geobasemap streets-light%landcover
f1.Position = [100 100 1200 900];
pause(1)
print(gcf,strcat('Overlay_Single_',num2str(sim_number),'.png'),'-dpng','-r300')
pause(1)
close(f1)
toc;
end