如何拟合一个封闭的曲面

cichishia

请教大家：
# R( ~2 |+ G+ w% v) p  `已知若干个离散点的三维坐标，如何拟合成一个封闭的曲面。

程序：
%--------------------------------------------------------------------------
5 p- g7 O) T: D7 z5 V1 YA=xlsread('PointData.xls');
1 F6 {! C& Y3 G' v' g* h% 下部分网格曲面
x=A(1:85,1);   % 根据dz的大小，相应调整行的范围--及确定好哪些点属于下部分
: I" [6 I. M' P7 M8 by=A(1:85,2);
- C& e. K% h1 p2 Qz=A(1:85,3);
minx = min(x);
% {$ W9 b! x1 c6 {6 \maxx = max(x);
miny = min(y);
maxy = max(y);
xtemp=linspace(min(x),max(x),50);%这个是选择的数值样本大小，越大曲线越平滑%
ytemp=linspace(min(y),max(y),50);
[X,Y]=meshgrid(xtemp,ytemp);
: D$ f0 {1 G1 E* s# uZ=griddata(x,y,z,X,Y,'cubic');
figure('numbertitle','off','name','网格图')
suRF(X,Y,Z)
shading faceted

%--------------------------------------------------------------------------
% 上部分网格曲面
hold on
7 f$ D0 l) W0 T& v" ox=A(61:404,1);
3 s$ N: {) L# V* W2 ay=A(61:404,2);
z=A(61:404,3);
minx = min(x);
: `$ L3 m! m" J, u6 D% lmaxx = max(x);
/ `' E3 x, ]) m/ H9 Y+ a& a2 m3 M/ Iminy = min(y);
7 }! x# e7 E; a' R6 F9 Rmaxy = max(y);
7 o* H) `6 h. ]- a, M# I4 Lxtemp=linspace(min(x),max(x),50);%这个是选择的数值样本大小，越大曲线越平滑%
ytemp=linspace(min(y),max(y),50);
[X,Y]=meshgrid(xtemp,ytemp);
Z=griddata(x,y,z,X,Y,'cubic');
! Z7 G4 A- I  i1 D  c- ^% asurf(X,Y,Z)
! g% p5 q& |5 F/ p1 @shading faceted
& n8 [; }8 a; m& T) A
figure('name','散点图')
3 C, ~' _; C( H* V" c1 s; nplot3(A(:,1),A(:,2),A(:,3),'*')


我自己用上述方法分别拟合上下两部分，但是图中有缺口，如何才能得到封闭的图形呢？
2 X$ K. }7 z+ ~+ p

IBB-EUT

帮你顶一下

shuddkk

给出数据，我可以试试。

小白的白