具体实施方式

实施例1

本发明公开了一种航空影像空中三角测量作业方法，包括以下步骤：

步骤1：根据测区范围和地形情况划分多个子测区，划分的子测区数目为N；

步骤2：通过航拍对子测区进行影像提取，依次对已经提取的子测区影像进行空中三角测量处理；

在本实施例的步骤2这一步中，影像提取和三角测量处理工作同时接力进行看，例如，第一天的航拍结束后，将该天拍摄提取的影像提交，作为第一个子测区，立即开始进行空中三角测量处理，利用当天晚上和第二天白天的时间，完成这些影像的内业处理；第二天白天继续进行航拍，航拍结束后，也将该天拍摄的影像提交给内业，作为第二个子测区，开始进行空中三角测量处理；利用当天晚上和第二天白天的时间，完成这些影像的内业处理，如此反复直至所有子测区的提取和测量作业完毕。

步骤3：将多个子测区的处理后的数据合并，得到一个完整的测区。

优选地，步骤3包括以下步骤：

步骤301：对每个子测区的三角测量结果进行整理，对每个相机、每张影像和每个三角测量连接点进行编号，且编号不重复；

步骤302：求取每张影像的邻影像数量M；每一张影像的邻影像数量的求取方法为：

利用该张影像的特征点信息，进行影像检索，对该张影像检索前k1张近邻影像；如果该张影像有GPS信息，用GPS信息计算该张影像与所在子测区内的每张影像之间的欧式距离，保留距离最近的前k2张影像做为该张影像的近邻影像，该张影像的近邻影像数量为M，M＝k1+k2；

步骤303：判断影像是否为子测区之间的边界影像；本实施例中所应用的判断影像是否为子测区之间的边界影像的方法为：对于每张影像，检查它的M张近邻影像，如果至少有一张近邻影像来源于不同的子测区，则这张影像属于边界影像。

步骤304：对所有的边界影像及其近邻影像进行影像特征匹配，形成邻接关系矩阵；

步骤305：对步骤304中得到的所有邻接关系矩阵进行分区，并对分区进行独立编号，分区数量为P，对这P个子测区进行并行三角测量处理，得到这P个子测区的三角测量结果；作为优选方案，在本实施例中，编号从N+1开始，每下一个编号依次递增1，即，编号分别为N+1、N+2、N+3......以此类推至N+P；

步骤306：对N个子测区的三角测量结果和P个边界影像及其近邻影像匹配区域进行合并，形成一个整测区的空中三角测量结果。

实施例2

基于传统方法和实施例1的方法验证效果并进行对比分析，试验方法如下：

以一个20平方公里的测区为例，计划每天航拍2平方公里，用两种作业方案处理这个测区。

传统作业方案，航拍10天，第11天开始空三内业处理，需要7天完成空三内业处理，总共需要17天；需要的计算机节点数为11台。

采用实施例1的作业方案，采用接力式空中三角测量作业方案，航拍十天，每天航拍结束后即开始处理当天的子测区影像，最后再用1天时间完成子测区的合并，总共需要12天；需要的计算机节点数为7台。

相比于以前的作业方案，可以减少4台计算机，还可以节省5天的工期。