Help with 3D Data Contour Plot (Missing Values as Zeros) (2024)

h =

Figure (1) with properties: Number: 1 Name: '' Color: [1 1 1] Position: [671 558 577 433] Units: 'pixels' Use GET to show all properties

BB =

ColorBar with properties: Location: 'eastoutside' Limits: [1.0831 1.0848] FontSize: 18 Position: [0.8323 0.1109 0.0381 0.8152] Units: 'normalized' Use GET to show all properties

Warning: Duplicate data points have been detected and removed - corresponding values have been averaged.

h =

Figure (2) with properties: Number: 2 Name: '' Color: [1 1 1] Position: [671 558 577 433] Units: 'pixels' Use GET to show all properties

BB =

ColorBar with properties: Location: 'eastoutside' Limits: [0 1] FontSize: 18 Position: [0.8323 0.1109 0.0381 0.8152] Units: 'normalized' Use GET to show all properties

h =

Figure (3) with properties: Number: 3 Name: '' Color: [1 1 1] Position: [671 558 577 433] Units: 'pixels' Use GET to show all properties

BB =

ColorBar with properties: Location: 'eastoutside' Limits: [1.0839 1.0858] FontSize: 18 Position: [0.8323 0.1109 0.0381 0.8152] Units: 'normalized' Use GET to show all properties

y2(isnan(y2(:,4)), 4) = 0;

set(0,'DefaultaxesFontSize',20);

set(0,'DefaulttextFontsize',20);

set(0,'DefaultaxesFontName','Times-Roman');

set(0,'DefaulttextFontName','Times-Roman');

set(0,'defaulttextinterpreter','latex')

set(0, 'DefaultAxesLineWidth', 4)

set(0,'DefaultAxesFontWeight','bold')

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

T = readtable('HB2B_4298.csv','PreserveVariableNames',true);%reading the data

Vx = table2array(T(:,:));%converting to arrays

TTT = Vx(:,11:14);%Pick your columns based on the data.EG., dsapce columns number = 13

A = unique(TTT,'rows');

Acell = splitapply(@(x){x}, A,findgroups(A(:,2)));%split group based on column 2=vy

y1 = [Acell{1,:}];%first plane values

y2 = [Acell{2,:}];%second plane values

y3 = [Acell{3,:}];%third plane values

% y4 = [Acell{4,:}];%fourth plane values

y = {y1, y2, y3};

N = length(y);

rg = cell(N,1); cg = cell(N,1);

Rg = cell(N,1); Cg = cell(N,1);

Zg = cell(N,1);

h = figure;

hold on

for i= 1:N

rg{i} = linspace(min(y{i}(:,1)), max(y{i}(:,1)), 100);%x axis

cg{i} = linspace(min(y{i}(:,3)), max(y{i}(:,3)), 100);%z axis

[Rg{i}, Cg{i}] = meshgrid(rg{i}, cg{i});% x and z

F = scatteredInterpolant(y{i}(:,1), y{i}(:,3), y{i}(:,4), 'natural'); % Use natural method for interpolation

Zg{i} = F(Rg{i}, Cg{i});

contourf(Rg{i}, Cg{i}, Zg{i}, 'ZLocation',i, 'LineStyle','none');%plot contour

colormap("default");

BB = colorbar;

BB.Label.String = 'd spacing';%change

xlabel('X(mm)')%change

ylabel('Z(mm)')%change

end

view([40 40]) % adjust the view as needed

saveas(h,'FIG.png');%change the FIG name to whatever the name you need

set(0,'DefaultaxesFontSize',20);

set(0,'DefaulttextFontsize',20);

set(0,'DefaultaxesFontName','Times-Roman');

set(0,'DefaulttextFontName','Times-Roman');

set(0,'defaulttextinterpreter','latex')

set(0, 'DefaultAxesLineWidth', 4)

set(0,'DefaultAxesFontWeight','bold')

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

T = readtable('HB2B_4298.csv','PreserveVariableNames',true,'NumHeaderLines',41);%reading the data

A = T{:,11:14};%Pick your columns based on the data.EG., dsapce columns number = 13

Acell = splitapply(@(x){x}, A, findgroups(A(:,2)));%split group based on column 2=vy

y1 = [Acell{1,:}];%first plane values

y2 = [Acell{2,:}];%second plane values

y3 = [Acell{3,:}];%third plane values

% y4 = [Acell{4,:}];%fourth plane values

% Construct a new y2, the same size as y1, keeping the values of y2(:,4)

% where y2(:,1) and y2(:,3) exist in y1, making the others NaNs

[ism,idx] = ismember(y2(:,[1 3]),y1(:,[1 3]),'rows');

V = y2(ism,4);

N = size(y1,1);

y2 = [y1(:,1) y2(1,2)*ones(N,1) y1(:,3) NaN(N,1)];

y2(idx(ism),4) = V;

% make sure y1, y2, and y3 all have the same X and Z coordinates now

assert(isequal(y1(:,[1 3]),y2(:,[1 3]),y3(:,[1 3])))

% X and Z grid points

NX = numel(unique(y1(:,1)));

NZ = numel(unique(y1(:,3)));

% reshape X and Z to matrices

X = reshape(y1(:,1),NX,NZ);

Z = reshape(y1(:,3),NX,NZ);

y = {y1, y2, y3};

N = length(y);

h = figure;

hold on

for i = 1:N

% reshape values to a matrix the same size as X and Z

V = reshape(y{i}(:,4),NX,NZ);

if i == 2

% to accommodate the missing data in the second slice,

% copy the 1st and last rows so that contourf can

% render something

V([2 4],:) = V([1 5],:);

end

contourf(X,Z,V,'ZLocation',y{i}(1,2),'LineStyle','none');%plot contour

end

colormap("default");

BB = colorbar;

BB.Label.String = 'd spacing';%change

xlabel('X(mm)')%change

ylabel('Y(mm)')%change

ylabel('Z(mm)')%change

view([40 40]) % adjust the view as needed

saveas(h,'FIG.png');%change the FIG name to whatever the name you need