激光雷达点云滤波的topscan算法实现

激光雷达点云滤波的topscan算法实现，能计算精度，并绘制出地面点效果图。

A=load（‘sample （1）-多层较大建筑物.txt‘）;
SA=size（A）;
for numm=60:-4:10%将窗口大小在60*60到10*10之间循环变化
 a=topscan（nummA）;%调用topscan函数，得到特定窗口大小和阈值求得的地面点集合
 m=0;
 if （numm==60）
     c=a;%第一次循环的结果赋给c
 end
 if （numm~=60）
 c=setdiff（intersect（ac）0）;%每次循环的结果求交集，得到迭代到当下时的地面点集合
 end
 b=setdiff（[1:SA（1）]c）;%得到非地面点集合
 Sc=size（c）;
 Sb=size（b）;
 %地面点正确划分的个数和非地面点正确划分的个数之和m%
for i=1:Sc（2）
    if（A（c（i）4）==0）
        m=m+1;
    end
end
for i=1:Sb（2）
   if（A（b（i）4）==1）
      m=m+1;
   end
end
%求解精度，并判断精度是否满足条件，满足则停止迭代%
 PRE=m/SA（1）;
  if（PRE>=0.98）
    break;
 end
end
j=1;
%绘制计算出来的地面点%
for i=1:Sc（2）
x（j）=A（c（i）1）;
y（j）=A（c（i）2）;
z（j）=A（c（i）3）;
j=j+1;
end
plot3（xyz‘r‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         850  2014-04-10 20:43  LIARA.m
     文件     1316570  2007-01-27 22:09  sample （1）-多层较大建筑物.txt
     文件         797  2014-04-10 20:37  topscan.m