1.数字城市单体房屋点云立面3D特征提取方法，其特征在于，基于单体房屋立面点云数据的特点，对原始点云数据构建八叉树索引结构，在近似度评估约束的条件下对点云数据进行平面分割，通过面面相交和完整性规则约束提取单体房屋立面特征线；

第一，单体房屋点云立面3D分割：一是近似度评估约束的点云分割，包括近似度评估计算、八叉树结构索引、近似度评估约束下分割算法；二是点云数据平面3D拟合，采用改进的最小二乘平面3D拟合方法；

首先在获取单体房屋点云数据之前进行预处理，通过滤波分离地面点和非地面点，对地面点分类提取单体房屋点云，采用近似度评估约束的点云分割算法，首先对原始点云数据八叉树结构化，然后进行点的近似度评估计算，基于相似度迭代分割单体房屋点云数据，获取单体房屋平面，对最小二乘平面拟合进行去噪改进，将分割后得到的平面进一步拟合，得到更精确的平面模型参数；

第二，单体房屋立面3D特征线提取：一是立面3D特征线提取方法，包括立面3D数据结构设计、特征线提取、房屋立面特征点提取方法、精简立面3D特征线；二是单体房屋3D特征线优化，包括特征线平面相交、立面特征线规则化；

通过形状语法设计单体房屋平面的立面数据结构,判定面片所属类型，建立面片之间的拓扑关系；将单个平面点云子集内的所有数据点投影至该子集对应的平面空间中，计算点云子集的最小外包络矩形，构建二维网格结构实现立面特征线的提取，通过平面相交结合立面原始特征线提取单体房屋立面特征线并进行规则化处理；

近似度评估约束下分割算法：采用迭代约束，根据点近似度对每个三维点所属的结构单元进行判定，从而不断更新每个结构单元的信息；

近似度评估约束下点云分割算法步骤为：

第一步：设基于八叉树构建的结构单元集合C＝{s

第二步：对C中每一个子集c

第三步：令c

第四步：赋予所有空间特征向量不同的权重，对于任意一点q∈Q

第五步：判定近似度的大小，如果所有c

第六步：遍历所有的3D点，当一轮遍历结束后，点q所属的结构单元发生变化，则相应信息进行更新；当遍历结束后八叉树的叶子结点为空，则从C中删除该结点对应的点云子集；

第七步：循环迭代，记录相邻点云子集之间的关系，将法向量相同且具有相邻关系的点云子集聚集为一个集合，记为一个分割平面的点云数据，直到没有3D点发生变化；

改进的最小二乘平面3D拟合方法：基于概率方法剔除噪点，获取正确的模型参数，采用随机抽样方法对点云N进行采样，假设含有噪点的概率为k,采样次数m，其中m次采样全部都抽到噪点的概率为k第1步：从分割平面N中随机的选取一个数据点Q，计算点Q的K邻域；

第2步：以得到的K邻域中的点为基础，采用最小二乘平面拟合的方法计算拟合局部平面F；

第3步：计算N中所有点到拟合平面F的距离A

第4步：重复第1步至第3步m次，选取m次中A

第5步：计算N中全部点到第4步中确定的原始平面的距离标准偏差r，计算式为：第6步：当A

第7步：把保留后的平面点集记为N

第8步：判定标准差是否满足第6步，如果不满足，循环第1步至第7步，直到得到最佳拟合平面的参数a、b、c。

2.根据权利要求1所述数字城市单体房屋点云立面3D特征提取方法，其特征在于，近似度评估计算：综合空间3D位置、法向量两个空间几何特征进行点云数据近似度评估计算，赋予不同维度的特征向量不同权重，通过权重的调配得到意向分割结果，记两点为Q其中，

其中，t点云子集的分辨率，

a

其中，α和β分别是a

基于相似度的点云数据分割算法决定性因素是不同维度特征值的权重和判定相似度的临界值大小，单体房屋平面特征突出，分配给平面位置和法向量这两个空间特征值权重不同权重，判定相似度的临界值也有重要作用，将两个空间几何特征统一到相同量纲下，用数值表示点云数据之间多维度的特征差异，将相似度临界值的范围在[0,1]之间取值，经过多次分割结果确定最终的临界值。

3.根据权利要求1所述数字城市单体房屋点云立面3D特征提取方法，其特征在于，八叉树结构索引：将单体房屋点云所在的三维空间看作是2设八叉树的深度为a,则a＝max(n,m,1)+1，八叉树索引结构中的每一个叶子结点都有已经确定的26个邻域。

4.根据权利要求1所述数字城市单体房屋点云立面3D特征提取方法，其特征在于，为拟合出更多平面，把距离临界值p

1)两个平面夹角小于临界值δ；

2)两个平面的所有数据点分别到两个平面的距离和的差值小于一定临界值△a，即将两个平面合并为一个平面；

设采用改进的最小二乘法得到的两个平面H

5.根据权利要求1所述数字城市单体房屋点云立面3D特征提取方法，其特征在于，特征线提取：通过构建二维网格结构实现立面3D特征点的提取，首先将单个立面点云子集内的所有数据点投影至该子集对应的二维平面中，计算点云子集的最小外包络矩形，将矩形划分为规则网格，把包含数据点的网格单元与不包含数据点的网格单元的公共边作为分割线，距离分割线最近的点就是单体房屋点云立面的特征点；

提取的特征点包括内部特征点和外部特征点两种，外部特征点代表的立面墙面的特征，内部特征点由立面墙面内部的空洞提取得到，对特征点连线进行直线拟合及规则化，使其能满足真实单体房屋的情况，首先计算每条折线与Y坐标轴的夹角，夹角小于临界值范围则认为折线与Y轴平行，并对这些折线进行平行矫正；如果折线与Y轴之间的夹角在垂直误差范围内，则对这些折线进行垂直矫正。

6.根据权利要求1所述数字城市单体房屋点云立面3D特征提取方法，其特征在于，房屋立面特征点提取方法：步骤一：选取点云子集集合C＝{c

步骤二：在立面点集c

步骤三：设置比例临界值e，计算原始特征与最小外包络矩形面积的比值，保留比值大于e的多边；

步骤四：设置角度临界值u,计算曲线中每段折线段与Y坐标轴的夹角，夹角小于u则认为折线与Y轴平行，并对这些折线进行平行矫正；同理,如果折线与X轴之间的夹角在[-u，+u]区间内，则对这些折线段进行平行X轴校核矫正；

步骤五：重复步骤四至步骤四，直到所有的C中所有的子集完成平行校核；

步骤六：判定完成平行校核后的立面的类型，如果立面的特征点连线较短，且位于立面内部，则这部分特征点属于内部特征点，其它的为外部特征点；

步骤七：将在二维空间中得到的立面3D特征点转化到三维空间,获取并记录对应的三维空间坐标值。

7.根据权利要求1所述数字城市单体房屋点云立面3D特征提取方法，其特征在于，精简立面3D特征线：对这些特征点连线进行直线拟合，以更加简洁的线段表达立面特征；

(1)特征点分组：对规则多边形，拐角处的方位角较大，基于拐角的方位变化将位于同一特征的点分为一组，首先在平行校核后的折线段上选取一点，判断其与相邻两点构成的两个向量的夹角，夹角超过临界值则为特征拐角点，设夹角临界值为B，点Q、Q(2)特征线拟合：采用最小二乘法将得到的不同特征线点组进行直线拟合，直线方程在笛卡尔坐标系以x或者y为自变量获取直线方程，设每个特征线点组中点的个数为m，分别计算它们的x、y均值和平方均值，基于分组的特征线点，得到直线拟合结果。

8.根据权利要求1所述数字城市单体房屋点云立面3D特征提取方法，其特征在于，单体房屋3D特征线优化：采用单体房屋不同立面之间存在相互平行、垂直拓扑关系，对原始特征线中不能满足垂直或平行要求的进行垂直和平行校核,使之满足几何形状的要求；

特征线平面相交：得到单体房屋立面直线化后的特征及立面之间拓扑关系后，首先基于立面数据结构中Nei_ID字段提取相邻立面信息，采用相交立面平面方程计算相交线并进行竖直方向调整，对单体房屋立面竖直方向特征线进行优化；

设两个平面方程Q

9.根据权利要求1所述数字城市单体房屋点云立面3D特征提取方法，其特征在于，立面特征线规则化：单体房屋的特征线是具有规则的几何形状，规则化调整立面原始特征线的长度、方位使其能准确的符合实际单体房屋形状，单体房屋的立面特征线大部分是平行或者垂直的正交关系，通过单体房屋的主方向对原始特征线进行规则化调整；

步骤1，确定主方向：将拟合边界中最长的线段作为原始主方向，然后将拟合边界集合中与原始主方向在临界值范围内的拟合边界提取出来，然后根据长度采用加权平均得出单体房屋的主方向；

设特征线段集合为H＝{H1,H2,…,H

步骤2，调整特征线：根据主方向调整边特征线段,将特征线按方位角分为三类，包括与主方向近似平行的、与主方向近似垂直的和既不平行也不垂直的；

设主方向的角度临界值为△β，获取H中方位角在(D-△β，D+△β)的特征线段集A＝{A步骤3，根据主方向调整后的特征线顺序求交，获取新特征交点，更新立面数据结构中的特征点字段，依次连接新的特征交点得到规则化后的单体房屋立面3D特征线。