%%?讀數據
[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