Minimizing Maximum Distance between 3D Point Cloud Data

4 次查看（过去 30 天）
Nicholas
Nicholas 2014-6-4
评论： Matt J 2014-6-7
I have what I'll call a "nominal point cloud" and then a "measured point cloud". My goal is that I want to minimize the maximum distance between respective points of EVERY point in the point cloud without moving points in the nominal cloud. I do not want a least-squares type fit that could sacrifice a couple points for the good of the rest of the cloud.
EDIT: Each point cloud contains the exact same amount of points and when transforming the measured point cloud, I need to move every point in relation to one another. i.e. If you move one point in a direction, every point moves in the same direction with the same magnitude.
I believe I could write a brute-force, majorly iterative, code for this process but the way I have in mind would take absolutely forever to run and that's not feasible. Looking for any advice available.
Thank you.
  7 个评论
显示 5更早的评论
José-Luis
José-Luis 2014-6-4
Also, if the distance is pairwise (the distance between one point in one cloud to one point in another cloud, not the minimum distance between a point in one cloud and all the points in other clouds) then your problem reduces to optimizing a system of 6 unknowns:
rotation (x, y and z)
translation (x, y and z)
So you need to find the six values that will minimize the distance between a pairs of points. That sounds like a job for fminsearch() and similar functions.
Nicholas
Nicholas 2014-6-5
Yes, sorry, an affine transform - like you said. Yes, the distance will be pair wise. To think of an example: I have a very thin, curved piece of metal, with holes in it. Each hole center point represents the points in my point cloud. Based off the CAD drawing, I know my center points should be at their nominal locations. However, when I fabricate the part, based on the tolerances, the holes will be in different locations when I measure the part, thus the measured point cloud. End result is I want to find a way to get the measured hole center locations to get as close to their respective nominal center locations as possible, for every hole center.
I'll look into fminsearch() like you suggested. Thank you for taking the time to consider the issue.

请先登录，再进行评论。

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

回答（1 个）

Matt J
Matt J 2014-6-4
编辑：Matt J 2014-6-4
My goal is that I want to minimize the maximum distance between respective points of EVERY point in the point cloud without moving points in the nominal cloud. I do not want a least-squares type fit that could sacrifice a couple points for the good of the rest of the cloud.
Well, the least squares method that you say you don't want is implemented here
http://www.mathworks.com/matlabcentral/fileexchange/26186-absolute-orientation-horn-s-method
I would do this to at least initialize a search using fminsearch, or better yet, fminimax if you have the Optimization Toolbox. This file may help to parameterize the rototranslation
http://www.mathworks.com/matlabcentral/fileexchange/30864-3d-rotation-about-shifted-axis
reg= absor(nominal,measured);
p0=[reg.q(2:3);0;reg.t];
options = optimoptions('fminimax','MinAbsMax',size(nominal,2));
p = fminimax(@(p) objective(p,nominal,measured),[],[],[],[],[],...
[],[],options);
function F=objective(params,nominal,measured)
x1=params(1);
x2=params(2);
x3=sqrt(1-x1^2-x2^2);
theta=params(3);
translation=params(4:6);
Tnominal=AxelRot(nominal,theta,[x1,x2,x3],translation);
F=Tnominal-measured;
  4 个评论
显示 2更早的评论
Nicholas
Nicholas 2014-6-5
Ah, I understand the flow of your answer better now. Thanks. I've read through both those functions and their uses and I believe they may help solve this problem. Again, awaiting licensure from my company to allow me to delve into the offered solution.
When the time comes, do you respond to e-mails? I may find it more beneficial to speak with you through that medium.
Regards.
Matt J
Matt J 2014-6-7
No, I don't do emails, I'm afraid. For one thing, it defeats the purpose of this forum if people move the discussion to private email exchange. The idea is to see the question resolved here so that not only you benefit as the OP, but also other people who might encounter the same problem, as well.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Point Cloud Processing 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by