2

私は現在、4つの別々のグラフを生成するMatlabスクリプト(以下)を実行しています。これらのグラフを2つ組み合わせて、が同時に表示されるようにします。結合したいグラフは、以下のスクリプトで図2および図4と呼ばれます。

スクリプトへの唯一の入力は、6列のtxtファイルです:x座標、y座標、および4つの変数(深さ[m]、Hsig [m]、周期[秒]、方向[度])

これに関する助けをいただければ幸いです。

% Post-process a SWAN wave model output file
%----------------------------------------------------------------------------------------------
%defaults
N_header = 7; % header lines in SWAN file
N_vars = 6;           %output variables in SWAN file
x_origin = 0;    %real world x origin
y_origin = 0;   %real world y origin
quiver_subsample = 6; %sub-sampling factor to make direction plot clearer
rot_angle = 0; %rotation angle to correct any previous rotation for SWAN
island_mask = load('island_mask.txt'); %mask for islands (set land to NaN);

%specify input file
[filename,pathname] = ...
  uigetfile('*.txt', 'Specify SWAN results file (e.g. Scilly.txt) [*.txt]');
  SWANfile = fullfile(pathname,filename);

%read (and ignore) file header lines
fid = fopen(SWANfile);
for i=1:N_header
 head = fgets(fid);
end
%and now get the data
data = fscanf(fid,'%g %g',[N_vars inf]); data = data';
fclose(fid);

%extract the datasets we want, marking any junk values (e.g. dry land)
XP = data(:,1);
YP = data(:,2);

DEPTH = data(:,3);
dudsDEPTH = (DEPTH==-99);
DEPTH(dudsDEPTH) = NaN;

HS = data(:,4);
dudsHS = (HS==-9);
HS(dudsHS) = NaN;

PER = data(:,5);
dudsPER = (PER==-9);
PER(dudsPER) = NaN;

DIR = data(:,6);
dudsDIR = (DIR==-999);
DIR(dudsDIR) = NaN;

minX = min(XP);
minY = min(YP);
maxX = max(XP);
maxY = max(YP);
cellsize = XP(2) - XP(1);

% mesh and plot data onto scaled output grids
[xp,yp] = meshgrid(minX:cellsize:maxX,minY:cellsize:maxY);

sx = size(xp); xlen = sx(2); ylen = sx(1);
depth = reshape(DEPTH,xlen,ylen);
hs = reshape(HS,xlen,ylen);
per = reshape(PER,xlen,ylen);
dir = reshape(DIR,xlen,ylen);

depth_rot = flipud(rot90(depth,1));%pcolor(depth_rot);shading flat
hs_rot = flipud(rot90(hs,1));%pcolor(hs_rot);shading flat
per_rot = flipud(rot90(per,1));%pcolor(per_rot);shading flat
dir_rot = flipud(rot90(dir,1));
%remember that actual directions also need rotating (i.e. not just matrix!)
dir_rot = dir_rot + rot_angle; %pcolor(dir_rot);shading flat

xp_rot = xp;
yp_rot = yp;

%create x and y matrices in real world co-ordinates
xp_rot = xp_rot + x_origin;
yp_rot = yp_rot + y_origin;

%and equivalent x and y vectors, in case we need these instead
grid_cells = size(xp_rot);
x_cells = grid_cells(2); %columns
y_cells = grid_cells(1); %rows
x_utm = x_origin:cellsize:x_origin + (x_cells*cellsize);
y_utm = y_origin:cellsize:y_origin + (y_cells*cellsize); 
% y_utm = fliplr(y_utm); % flip to ensure cartesian rather than image axes


%create bathymetry plot
figure(1)
if ~isempty(island_mask)
 depth_rot_plot = depth_rot;
 depth_rot_plot(island_mask) = NaN;
 imagesc(x_utm,y_utm,depth_rot_plot)
 colormap(jet(256));
 map = colormap;
 map(1,:) = 1;
 % map(2,:) = 1;
 % map(3,:) = 1;
 colormap(map); 
else
 imagesc(x_utm,y_utm,depth_rot)
end
title('Bathymetry (m)', 'fontsize', 12)
set(gca,'fontsize', 12);
axis equal
axis tight
axis xy % need this to ensure cartesian rather than image axes!
colorbar

%create direction plot
figure(2)
[U,V] = pol2cart((dir_rot) ./ (180/pi),ones(size(dir_rot)));
% and sub-sample output grid to produce clearer plot
U_subsample = nestedsubsample2(U,quiver_subsample);
V_subsample = nestedsubsample2(V,quiver_subsample);
X_subsample = nestedsubsample2(xp,quiver_subsample);
Y_subsample = nestedsubsample2(yp,quiver_subsample);
quiver(X_subsample,Y_subsample,U_subsample,V_subsample,'k');
title('Direction')
axis equal
axis tight

%visualise wave breaking by taking ratio of Hs and depth
breaking = hs_rot ./ depth_rot;
breaking(breaking > 0.7) = 0.70;

%create Hb/depth plot to show where waves are shoaling and/or breaking
figure(3)
imagesc(x_utm,y_utm,breaking)
title('Hs / depth', 'fontsize', 12)
set(gca,'fontsize', 12);
axis equal
axis tight
axis xy % need this to ensure cartesian rather than image axes!
colorbar


%create Hs plot
figure(4)
imagesc(x_utm,y_utm,hs_rot)
title('Hs (m)', 'fontsize', 12)
set(gca,'fontsize', 12);
axis equal
axis tight
axis xy % need this to ensure cartesian rather than image axes!
colorbar


function [subsampled_A] = nestedsubsample2(A,I)
%NESTEDSUBSAMPLE2(A,I)
%resample 2D matrix A by retaining every Ith element

if I < 1
 I = 1;
end

A_dim = size(A);
new_i = 1:I:A_dim(1);
new_j = 1:I:A_dim(2);

subsampled_A = A(new_i,new_j);

end


end
4

3 に答える 3

10

複数のグラフを同じ図の別々の軸に表示する場合は、 SUBPLOTコマンドを使用できます。これにより、Figure ウィンドウ内で軸を並べて表示できます。

複数のグラフを同じ軸に表示する場合は、 HOLDコマンドを使用できます。

于 2010-01-25T18:04:05.790 に答える
1

subplot(m,n,p)しばらくの間私を混乱させたので、他の誰かが
同じ精神的ブロックを持っていた場合に備えて、私は私の2セントを入れると思いました

the m,n values set up numbers of rows and columns within a figure
the p   value  determines which subplot you are accessing 
               subplots are numbered   1 ..  n in row 1
                                     n+1 .. 2n in row 2
                                    2n+1 .. 3n in row 3   etc

たとえば、3行2列の図のセットの場合、次のように各場所にアクセスします。

+-----------------------------------+
| subplot(3,2,1)  |  subplot(3,2,2) |
|-----------------+-----------------+
| subplot(3,2,3)  |  subplot(3,2,4) |
|-----------------+-----------------+
| subplot(3,2,5)  |  subplot(3,2,6) |
|-----------------------------------+

各呼び出しの後に、通常の、、、などの呼び出しsubplot(m,n,p)を行うことができ 、指定されたサブプロットの位置にプロットが描画されますplot()title()legend()

于 2012-09-30T17:56:34.733 に答える
0

これは、gnovice が推奨する subplot の使用法を簡単に実装したものです。

たとえば、図 3 と図 4 を同じ図ウィンドウに並べて表示するとします。

figure(3)

%create Hb/depth plot to show where waves are shoaling and/or breaking
subplot(1,2,1)
imagesc(x_utm,y_utm,breaking)
title('Hs / depth', 'fontsize', 12)
set(gca,'fontsize', 12);
axis equal
axis tight
axis xy % need this to ensure cartesian rather than image axes!
colorbar

%create Hs plot
subplot(1,2,2)
imagesc(x_utm,y_utm,hs_rot)
title('Hs (m)', 'fontsize', 12)
set(gca,'fontsize', 12);
axis equal
axis tight
axis xy % need this to ensure cartesian rather than image axes!
colorbar

現在、MATLAB コンピューターにアクセスできないため、確認することはできませんが、上記のコードが正常に動作すると仮定すると、これも正常に動作するはずです。

于 2010-01-25T22:03:32.673 に答える