Matlab 版 (精华区)

发信人: zjliu (秋天的萝卜), 信区: Matlab
标  题: 想找出一个(平面)阵列（n*2）的边界 v2.0
发信站: 哈工大紫丁香 (Fri Apr 22 15:42:00 2005), 转信

续前贴：
load data.mat;
X=data_org(:,1);
Y=data_org(:,2);
m=round((X-min(X))/0.3910+1);
n=round((Y-min(Y))/0.3910+1);
A=zeros(max(n),max(m));
AA=zeros(size(A)+2);
A(sub2ind([size(A)],n,m))=1;
AA(2:end-1,2:end-1)=A;
C=(AA(2:end-1,1:end-2)+AA(2:end-1,3:end)+...
    AA(1:end-2,2:end-1)+AA(3:end,2:end-1)).*AA(2:end-1,2:end-1);
p=find(C>0&C<4);
[y,x]=ind2sub([size(C)],p);
Px=(x-1)*0.3910+min(X);
Py=(y-1)*0.3910+min(Y);
plot(X,Y,'o')
hold on
plot(Px,Py,'*')

--
  我的友情测试更新了，欢迎测试!
    有事找我可以到P_zjliu版给我留言，那里是我的地盘
╔═══════════════════╗
║★★★★★友谊第一  比赛第二★★★★★║
╚═══════════════════╝


※ 来源:·哈工大紫丁香 bbs.hit.edu.cn·[FROM: 202.118.229.162]