%%?讀數(shù)據(jù)
[fnamedirpath]=uigetfile（‘*.txt‘‘文件打開對話框‘）;
ffname=[dirpath?fname];
Point?=?importdata（ffname）;
%
%%?顯示原地形地貌地物
figure（1）
plot3（Point（:1）Point（:2）Point（:3））;
%顯示原始地形
GPoint?=?find（Point（:4?）==?1）;
%figure（2）
%plot3（Point（GPoint（:）1）Point（GPoint（:）2）Point（GPoint（:）3））;
%%?初始化
SizeXY=size（Point（:1:2）1）;%獲取列數(shù)
%實際的地面點數(shù)‘
numRG?=?size（?GPoint1）
ThreString?=?inputdlg（‘請輸入tan閾值‘‘TAN閾值‘）;
Thre?=?str2num（ThreString{1}）;
ThreString1?=?inputdlg（‘請輸入高程閾值‘‘閾值‘）;
ThreH?=?str2num（ThreString1{1}）;
%對組數(shù)進(jìn)行一維二向?
DelPoint1?=?0;
DelPoint2?=?0;
Point（16）=1;Point（15）=1;
%%?數(shù)據(jù)處理
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