点云滤波 一维双向扫描方法

滤波方法 一维双向扫描方法 武汉大学专业课激光遥感作业滤波方法 一维双向扫描方法 武汉大学专业课激光遥感作业滤波方法 一维双向扫描方法 武汉大学专业课激光遥感作业

%% 读数据
[fnamedirpath]=uigetfile（‘*.txt‘‘文件打开对话框‘）;
ffname=[dirpath fname];
Point = importdata（ffname）;
%% 显示原地形地貌地物
figure（1）
plot3（Point（:1）Point（:2）Point（:3））;
%显示原始地形
GPoint = find（Point（:4 ）== 1）;
%% 初始化
SizeXY=size（Point（:1:2）1）;%获取列数
%实际的地面点数‘
numRG = size（ GPoint1）
ThreString = inputdlg（‘请输入tan阈值‘‘TAN阈值‘）;
Thre = str2num（ThreString{1}）;
ThreString1 = inputdlg（‘请输入高程阈值‘‘阈值‘）;
ThreH = str2num（ThreString1{1}）;
%对组数进行一维二向 
DelPoint1 = 0;
DelPoint2 = 0;
Point（16）=1;Point（15）=1;
%% 数据处理
for i = 2:SizeXY                                              
    ThreH1 = Point（i 3） - Point（i - DelPoint1-13）;%相邻点高程差
    ThreH2 = Point（SizeXY-i+1 3） - Point（SizeXY-i+2 + DelPoint23）;
    ThreS1 = sqrt（（Point（i 2） - Point（i - DelPoint1-12））^2 ...%相邻点水平差
        + （Point（i 1） - Point（i - DelPoint1-11））^2）;
    ThreS2 = sqrt（（Point（SizeXY-i+1 2） - Point（SizeXY-i+2 + DelPoint22））^2 ...
        + （Point（SizeXY-i+1 1） - Point（SizeXY-i +2+ DelPoint21））^2）;
    ThreTan1 = ThreH1/ThreS1;%求得tan值
    ThreTan2 = ThreH2/ThreS2;
     if abs（ThreTan1）>Thre&&ThreH1        Point（i5）=0;%地面点
        DelPoint1 = DelPoint1 + 1; 
    else
        Point（i5）=1;%非地面点
     end
     if  abs（ThreTan2）>Thre&&ThreH2        Point（i6）=0;
        DelPoint2 = DelPoint2 + 1;
     else
        Point（i6）=1;
    end
end

%实际的找到的点数
Point（:7）=1;%第七列为判定值1为地面点
a = find（（Point（:5 ）== 1）&（Point（:6 ）== 1））;
Point（a7）=0;
TrueNum =  size（find（（Point（:4 ）== 1&Point（:7 ）== 1）|（Point（:4 ）== 0&Point（:7 ）== 0））1）;
rate = double（TrueNum）/double（SizeXY）;
%% 输出及拟合
msgbox（{‘最后的正确率为：‘num2str（rate）}）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-11-06 08:43  一维双向扫描滤波\
     文件        1798  2017-11-06 08:43  一维双向扫描滤波\Filtering.m
     目录           0  2017-11-06 08:32  一维双向扫描滤波\测试数据\
     文件     1178310  2011-11-27 21:38  一维双向扫描滤波\测试数据\samp11.txt
     文件     1615689  2011-11-27 21:38  一维双向扫描滤波\测试数据\samp12.txt
     文件      401760  2011-11-27 21:32  一维双向扫描滤波\测试数据\samp21.txt
     文件     1013886  2011-11-27 21:32  一维双向扫描滤波\测试数据\samp22.txt
     文件      777945  2011-11-27 21:32  一维双向扫描滤波\测试数据\samp23.txt
     文件      232252  2011-11-27 21:32  一维双向扫描滤波\测试数据\samp24.txt
     文件      894722  2011-11-27 21:32  一维双向扫描滤波\测试数据\samp31.txt
     文件      348161  2011-11-27 21:32  一维双向扫描滤波\测试数据\samp41.txt
     文件     1316570  2011-11-27 21:38  一维双向扫描滤波\测试数据\samp42.txt
     文件      553195  2011-11-27 21:32  一维双向扫描滤波\测试数据\samp51.txt
     文件      696694  2011-11-27 21:32  一维双向扫描滤波\测试数据\samp52.txt
     文件     1065718  2011-11-27 21:38  一维双向扫描滤波\测试数据\samp53.txt
     文件      266848  2011-11-27 21:32  一维双向扫描滤波\测试数据\samp54.txt
     文件     1086860  2011-11-27 21:38  一维双向扫描滤波\测试数据\samp61.txt
     文件      484995  2011-11-27 21:32  一维双向扫描滤波\测试数据\samp71.txt