publicImage/gprMax
1
0
派生 0
镜像自地址 https://gitee.com/sunhf/gprMax.git 已同步 2025-08-07 15:10:13 +08:00
代码 工单 软件包 项目 版本发布 百科 动态

Improved loading and selecting of correct EM field component.

浏览代码
这个提交包含在：
Craig Warren
2016-10-14 15:05:00 +01:00
父节点 000cb46c4c
当前提交 3e5515d07f
共有 1 个文件被更改，包括 34 次插入 和 35 次删除
下载 Patch 文件 下载 Diff 文件 展开所有文件 折叠所有文件

69
tools/MATLAB scripts/plot_Bscan.m
查看文件

@@ -12,43 +12,42 @@ filename = fullfile(pathname, filename);
if filename ~= 0 if filename ~= 0
iterations = double(h5readatt(filename, '/', 'Iterations')); iterations = double(h5readatt(filename, '/', 'Iterations'));
dt = h5readatt(filename, '/', 'dt'); dt = h5readatt(filename, '/', 'dt');
Ex = h5read(filename, '/rxs/rx1/Ex')';
Ey = h5read(filename, '/rxs/rx1/Ey')';
Ez = h5read(filename, '/rxs/rx1/Ez')';
end
prompt = 'Which field do you want to view? Ex, Ey, or Ez: '; prompt = 'Which field do you want to view? Ex, Ey, or Ez: ';
field = input(prompt,'s'); field = input(prompt,'s');
field = eval(field); fieldpath = strcat('/rxs/rx1/', field);
time = 0:dt:iterations * dt; field = h5read(filename, fieldpath)';
traces = 0:size(field,2); time = 0:dt:iterations * dt;
traces = 0:size(field,2);
fh1=figure('Name', filename); fh1=figure('Name', filename);
clims = [-max(max(abs(field))) max(max(abs(field)))]; clims = [-max(max(abs(field))) max(max(abs(field)))];
im = imagesc(traces, time, field, clims); im = imagesc(traces, time, field, clims);
xlabel('Trace number'); xlabel('Trace number');
ylabel('Time [s]'); ylabel('Time [s]');
c = colorbar; c = colorbar;
c.Label.String = 'Field strength [V/m]'; c.Label.String = 'Field strength [V/m]';
ax = gca; ax = gca;
ax.FontSize = 16; ax.FontSize = 16;
xlim([0 traces(end)]); xlim([0 traces(end)]);
% Options to create a nice looking figure for display and printing % Options to create a nice looking figure for display and printing
set(fh1,'Color','white','Menubar','none'); set(fh1,'Color','white','Menubar','none');
X = 60; % Paper size X = 60; % Paper size
Y = 30; % Paper size Y = 30; % Paper size
xMargin = 0; % Left/right margins from page borders xMargin = 0; % Left/right margins from page borders
yMargin = 0; % Bottom/top margins from page borders yMargin = 0; % Bottom/top margins from page borders
xSize = X - 2*xMargin; % Figure size on paper (width & height) xSize = X - 2*xMargin; % Figure size on paper (width & height)
ySize = Y - 2*yMargin; % Figure size on paper (width & height) ySize = Y - 2*yMargin; % Figure size on paper (width & height)
% Figure size displayed on screen % Figure size displayed on screen
set(fh1, 'Units','centimeters', 'Position', [0 0 xSize ySize]) set(fh1, 'Units','centimeters', 'Position', [0 0 xSize ySize])
movegui(fh1, 'center') movegui(fh1, 'center')
% Figure size printed on paper % Figure size printed on paper
set(fh1,'PaperUnits', 'centimeters') set(fh1,'PaperUnits', 'centimeters')
set(fh1,'PaperSize', [X Y]) set(fh1,'PaperSize', [X Y])
set(fh1,'PaperPosition', [xMargin yMargin xSize ySize]) set(fh1,'PaperPosition', [xMargin yMargin xSize ySize])
set(fh1,'PaperOrientation', 'portrait') set(fh1,'PaperOrientation', 'portrait')
end