I have 3d points and I should determine a method to produce and visualize a "center line" down the central axis

4 次查看（过去 30 天）

lalik lavender 2018-8-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/413107-i-have-3d-points-and-i-should-determine-a-method-to-produce-and-visualize-a-center-line-down-the-c

评论： jordan Ashton 2022-3-28

Load and visualize the file data (attached). This is data cloud of 3d points . You will want to examine the first "field" coordinates associated with each "node"(x,y,z). The units are in mm. Load the second field "error". Determine a method to produce and visualise a "center line" down the central axis of figure. Determine a method to calculate the radius/diameter of the figure at a number of points down this center line. Provide an estimate of the overall figure "length" and "volume".

4 个评论
显示 2更早的评论隐藏 2更早的评论

jonas 2018-8-8

What is going on in this thread. Are you trying to delete the question? It's not very polite when Image Analyst gave you so much help.

Rena Berman 2018-8-29

(Answers Dev) Restored edit

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Image Analyst 2018-8-1

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/413107-i-have-3d-points-and-i-should-determine-a-method-to-produce-and-visualize-a-center-line-down-the-c#answer_331238

I'd use plot3() to plot the points. Then I'd use pca() to find the central principal component axes, then use plot3() again to plot those axes through the points.

http://www.mathworks.com/matlabcentral/answers/8626-how-do-i-get-help-on-homework-questions-on-matlab-answers

11 个评论
显示 9更早的评论隐藏 9更早的评论

Image Analyst 2018-8-3

在 MATLAB Online 中打开

dataex.xlsx

See adjusted code to plot the PC axes on the original XYZ data:

clc;    % Clear the command window.
clearvars;
close all;  % Close all figures (except those of imtool.)
workspace;  % Make sure the workspace panel is showing.
format long g;
format compact;
fontSize = 22;
% Read data.
data = xlsread('dataex.xlsx')
% Plot data.
subplot(1, 2, 1);
x = data(:, 1);
y = data(:, 2);
z = data(:, 3);
plot3(x, y, z, 'b.', 'MarkerSize', 8);
grid on;
title('Original Data', 'FontSize', 20);
xlabel('x', 'FontSize', 20);
ylabel('y', 'FontSize', 20);
zlabel('z', 'FontSize', 20);
% Enlarge figure to full screen.
set(gcf, 'Units', 'Normalized', 'OuterPosition', [0, 0.04, 1, 0.96]);
% Get the PC coefficients.
coeffs = pca(data(:, 1:3))
pcs = data(:, 1:3) * coeffs
% Get the data transformed into PC space.
pc1 = pcs(:, 1);
pc2 = pcs(:, 2);
pc3 = pcs(:, 3);
subplot(1, 2, 2);
plot3(pc1, pc2, pc3, 'r.', 'MarkerSize', 8);
grid on;
title('Transformed Data', 'FontSize', 20);
xlabel('PC #1', 'FontSize', 20);
ylabel('PC #2', 'FontSize', 20);
zlabel('PC #3', 'FontSize', 20);
% Get the max and min of each PC
pc1Min = min(pc1)
pc2Min = min(pc2)
pc3Min = min(pc3)
pc1Max = max(pc1)
pc2Max = max(pc2)
pc3Max = max(pc3)
pc1Mean = mean(pc1)
pc2Mean = mean(pc2)
pc3Mean = mean(pc3)
xMean = mean(x)
yMean = mean(y)
zMean = mean(z)
% Since pcs = xyz * coeffs, to go the other way
% xyz = pcs / coeffs
pc1Limits = [pc1Min - pc1Mean, 0, 0; pc1Max - pc1Mean, 0, 0]
pc2Limits = [0, pc2Min - pc2Mean, 0; 0, pc2Max - pc2Mean, 0]
pc3Limits = [0, 0, pc3Min - pc3Mean; 0, 0, pc3Max - pc3Mean]
xyz1Limit = pc1Limits / coeffs
xyz2Limit = pc2Limits / coeffs
xyz3Limit = pc3Limits / coeffs
% Switch back to subplot 1.
subplot(1, 2, 1);
% Plot the lines along the principal axes.
% Need to add the mean of x, y, and z back in to center the axes.
hold on;
plot3(xyz1Limit(:, 1) + xMean, xyz1Limit(:, 2) + yMean, xyz1Limit(:, 3) + zMean, 'k-', 'LineWidth', 4);
plot3(xyz2Limit(:, 1) + xMean, xyz2Limit(:, 2) + yMean, xyz2Limit(:, 3) + zMean, 'k-', 'LineWidth', 4);
plot3(xyz3Limit(:, 1) + xMean, xyz3Limit(:, 2) + yMean, xyz3Limit(:, 3) + zMean, 'k-', 'LineWidth', 4);