thesis/plot_boundaries.m at main · julesvanirsel/thesis · 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
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
ftn = 'Arial';
fts = 10*2;
lw = 1.4;

close all
reset(0)
set(0,'defaultFigurePaperUnits','inches')
set(0,'defaultTiledlayoutPadding','tight')
set(0,'defaultTiledlayoutTileSpacing','tight')
set(0,'defaultSurfaceEdgeColor','flat')
set(0,'defaultLineLineWidth',lw)
set(0,'defaultScatterLineWidth',lw)
set(0,'defaultQuiverLineWidth',lw*0.7)
jules.tools.setall(0,'FontName',ftn)
jules.tools.setall(0,'FontSize',fts)
jules.tools.setall(0,'Multiplier',1)

colorcet = @jules.tools.colorcet;

load('data\flux_grad.mat')
bound2 = bound;
load('data\cond_cont.mat')
bound1 = bound;

figure
set(gcf,'PaperPosition',[0,0,13.2,4.8])
tiledlayout(1,2)
ltr = 65;

nexttile
text(0.04,0.9,char(ltr),'units','normalized','FontSize',fts*0.8); ltr = ltr+1;
hold on
pcolor(X2_imag*scl.x,X3_imag*scl.x,arc.^ap*scl.arc)
plot(bound_pts*scl.x,bound1.A(bound_pts)*scl.x,'k')
plot(bound_pts*scl.x,bound1.B(bound_pts)*scl.x,'--k')
plot(bound_pts*scl.x,bound2.A(bound_pts)*scl.x,'r')
plot(bound_pts*scl.x,bound2.B(bound_pts)*scl.x,'--r')
colormap(colorcet(clm.arc))
clb = colorbar;
clb.Label.String = 'Pedersen conductance (S)';
clb.Location = 'northoutside';
xlabel(lbl.x); ylabel(lbl.y)
xlim(lim.x); ylim(lim.y)
pbaspect(ar)

nexttile
text(0.04,0.9,char(ltr),'units','normalized','FontSize',fts*0.8);
hold on
pcolor(X2_imag(2:end-1,2:end-1)*scl.x,X3_imag(2:end-1,2:end-1)*scl.x,edges)
plot(bound_pts*scl.x,bound.A(bound_pts)*scl.x,'k')
plot(bound_pts*scl.x,bound.B(bound_pts)*scl.x,'--k')
plot(bound_pts*scl.x,bound2.A(bound_pts)*scl.x,'r')
plot(bound_pts*scl.x,bound2.B(bound_pts)*scl.x,'--r')
colormap(colorcet(clm.arc))
clb = colorbar;
clb.Label.String = 'Sobel edges (a.u.)';
clb.Location = 'northoutside';
xlabel(lbl.x)
yticks([])
xlim(lim.x); ylim(lim.y)
pbaspect(ar)

filename = 'boundaries.png';
saveas(gcf,fullfile('plots','paper0',filename))
close all