How to remove or exclude the intersecting part from two intersecting or overlaping spheres using Matlab?

4 次查看（过去 30 天）

M.S. Khan 2020-4-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/514469-how-to-remove-or-exclude-the-intersecting-part-from-two-intersecting-or-overlaping-spheres-using-ma

编辑： darova 2020-4-4

Hi everyone,

i have two spheres of radii R1 and R2 with centers C1 and C2 intersecting each others. i want to exclude the intersecting part. How to remove it to get new object like peanut. i used following code to get interseting spheres.

3 个评论
显示 1更早的评论隐藏 1更早的评论

darova 2020-4-2

Thank you for sharing your thoughts with us!

M.S. Khan 2020-4-2

Thanks Mr. John for your kind feedback and guidance. Warm regards!

请先登录，再进行评论。

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

采纳的回答

darova 2020-4-1

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/514469-how-to-remove-or-exclude-the-intersecting-part-from-two-intersecting-or-overlaping-spheres-using-ma#answer_423341

在 MATLAB Online 中打开

Here is what i invented

ix1 = (xc1-x2).^2+(yc1-y2).^2+(zc1-z2).^2 < R2^2;
ix2 = (xc2-x1).^2+(yc2-y1).^2+(zc2-z1).^2 < R1^2;
scatter3(xc1(ix1),yc1(ix1),zc1(ix1),20,'r','fill')
scatter3(xc2(ix2),yc2(ix2),zc2(ix2),20,'r','fill')

Of course i prepared explanations and result picture

10 个评论
显示 8更早的评论隐藏 8更早的评论

M.S. Khan 2020-4-2

在 MATLAB Online 中打开

Screen Shot 2020-04-02 at 2.03.55 PM.png

% Hi Mr. Darova, i have followed your procedure but i cannot get the same figure as  you 
% generated. Please can  you check, where i am making error. My figure is attached.
% Thanks for all your cooperation.
R1 = 10; R2 = 14; C1=[10 12 14]; C2=[16 18 20];
%************** *************************** Point cloud of Sphere-1 *******************************
% This function handle was used to create the sphere of 10,000 random
% points.
numPoints = 10000;	
r = randn(3, numPoints);
r = bsxfun(@rdivide, r, sqrt(sum(r.^2,1)));
r = R1 * r;
% Extract the Xs1, Ys1, and Zs1 coordinates from the array for sphere S1
x1=C1(1,1);
y1=C1(1,2);
z1=C1(1,3);
x2 = C2(1,1);
y2 = C2(1,2);
z2 = C2(1,3);
Xs2 = r(1,:) + x2; % Extract Xs2 from row #1.
Ys2 = r(2,:) + y2; % Extract Ys2 from row #2.
Zs2 = r(3,:) + z2; % Extrac
% Xs1, Ys1, Zs1  represent coordiantes for sphere S1
Xs1 = r(1,:) + x1 ; % Extract Xs1 from row #1.
Ys1 = r(2,:) + y1; % Extract Ys1 from row #2.
Zs1 = r(3,:) + z1; % Extract Zs1 from row #3.
% Display the shell of points
figure(1)
 ix2 = (Xs2-x1).^2+(Ys2-y1).^2+(Zs2-z1).^2 < R1^2;
 scatter3(Xs2(ix2),Ys2(ix2),Zs2(ix2),20,'r','fill', 'MarkerFaceColor',[1 0 1])
hold on
scatter3(Xs1, Ys1, Zs1,'b')
% ix1 = (Xs1-x2).^2+(Ys1-y2).^2+(Zs1-z2).^2 < R2^2;
% scatter3(Xs1(ix1),Ys1(ix1),Zs1(ix1),20,'r','fill')
axis square; 
xlabel('X', 'FontSize', 20);
ylabel('Y', 'FontSize', 20);
zlabel('Z', 'FontSize', 20);
% Enlarge figure to full screen.
set(gcf, 'Units', 'Normalized', 'OuterPosition', [0 0 1 1]);
%
%msgbox('Now use the circular arrow icon on the toolbar to rotate the sphere.');  
%% ********************** Generation of sphere-2 of same 10,000 point cloud **************************************************************
numPoints = 10000;
% Get a 3-by-numPoints list of (Xs2, Ys2, zs2) coordinates for sphere S2.
r = randn(3, numPoints);
r = bsxfun(@rdivide, r, sqrt(sum(r.^2,1)));
r = R2 * r;
%C2= input('Centre of sphere 2 =  ')
x2 = C2(1,1);
y2 = C2(1,2);
z2 = C2(1,3);
% Xs2, Ys2, Zs2  represent coordiantes for sphere S2
Xs2 = r(1,:) + x2; % Extract Xs2 from row #1.
Ys2 = r(2,:) + y2; % Extract Ys2 from row #2.
Zs2 = r(3,:) + z2; % Extract Zs2 from row #3.
% Display the shell of points
ix1 = (Xs1-x2).^2+(Ys1-y2).^2+(Zs1-z2).^2 < R2^2;
scatter3(Xs1(ix1),Ys1(ix1),Zs1(ix1),20,'r','fill', 'MarkerFaceColor',[1 0 1])
hold on
scatter3(Xs2, Ys2, Zs2,'g');
axis square; 
xlabel('X', 'FontSize', 20);
ylabel('Y', 'FontSize', 20);
zlabel('Z', 'FontSize', 20);
% Enlarge figure to full screen.
set(gcf, 'Units', 'Normalized', 'OuterPosition', [0 0 1 1]);  

darova 2020-4-3

nuhrbuh92@gmail.com

M.S. Khan 2020-4-3

Thanks sir. if you dont mind, i will keep in touch to learn from you Matlab skills. Regards and obliged for your kind guidance.

请先登录，再进行评论。

类别

MATLAB Graphics Graphics Objects Graphics Object Programming

在 Help Center 和 File Exchange 中查找有关 Graphics Object Programming 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by

How to remove or exclude the intersecting part from two intersecting or overlaping spheres using Matlab?

3 个评论
显示 1更早的评论隐藏 1更早的评论

采纳的回答

10 个评论
显示 8更早的评论隐藏 8更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

How to remove or exclude the intersecting part from two intersecting or overlaping spheres using Matlab?

3 个评论 显示 1更早的评论隐藏 1更早的评论

采纳的回答

10 个评论 显示 8更早的评论隐藏 8更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

WeChat

3 个评论
显示 1更早的评论隐藏 1更早的评论

10 个评论
显示 8更早的评论隐藏 8更早的评论