本发明属于地形测绘技术领域，具体涉及一种基于物方反算的dem拼接方法及dsm拼接方法。
背景技术：
数字摄影测量技术被广泛地应用于从二维数字影像提取三维空间目标的空间几何信息，其中运用影像匹配技术获得密集点云，生成数字表面模型dsm或数字高程模型dem。数字高程模型dem只包括了地形高程信息，并未包括其它地表信息。数字表面模型dsm是指包括了地形高程信息和地表高程信息的地面高程模型，地表高程信息包括地表建筑物、桥梁和树木等高度的地面高程模型。数字表面模型dsm是在数字高程模型dem的基础上，进一步涵盖了除地面以外的其它地表信息的高程。数字表面模型dsm和数字高程模型dem广泛用于各行各业，在利用航空摄影获得一整片测区dem和dsm时，需要对相邻局部测区的dem和dsm进行拼接合并，由于每个立体像对之间存在系统误差等原因，使得相邻局部测区的dem和dsm之间会出现拼接痕迹，特别是在山区、丘陵等高程变化剧烈的区域会出现三角网拉伸现象，如果这种数据不剔除，拼接时无法判断边界的高程值是否为真实值，拼接后会出现明显的拼接痕迹，使得数据接边失败，影响精度。
在获取图像照片时，为了得到完整的测区图像照片，在获取图像照片时，相邻两个测区的图像照片必须保证有60％以上的重叠，相邻两个航带的图像照片必须保证30％以上的重叠，才能保证，相邻两个测区的dem和dsm至少有10％以上的重接，相邻两个航带的dem和sm至少有20％的重叠，才能避免测区遗漏。但是这样获得的小块数据拼接成整块数据，必然存在数据接边和去除冗余的问题。而在每个dsm或dem靠近边界的数据往往是精度较差的数据，很可能存在着粗差、内插的近似点等，这使得虽然dsm和dem有10％以上重叠，但是很难做到无缝对接，必须进入立体测图进行人工修改。对于大范围的测区这使得人工修改量巨大，数据生成周期变长，成本升高。这时就需要一种完全自动化不需要人工干预就可完成裁切拼接的系统来完成工作。
目前国内外的针对该问题的解决大部分是直接拼接或者按照几何距离进行裁切，如果处理的效果不好再基于人工修改。由于接边是进入立体测图进行人工修改接边，这使得数据处理效率低下不够智能化。
技术实现要素：
本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种基于物方反算的dem拼接方法及dsm拼接方法，其步骤简单，设计合理，采用dem和dsm均具有的地形表面坐标参与计算图像照片的像点坐标，计算机根据裁切边界对像点坐标进行筛选，由位于裁切边界内的像点坐标在dem中的映射构建图像照片裁切后的dem，由位于裁切边界内的像点坐标在dsm中的映射构建图像照片裁切后的dsm，对相邻的两个图像照片中具有同一水平坐标的地形表面坐标进行加权处理得到加权坐标，以加权坐标代替具有同一水平坐标的地形表面坐标，达到dem和dsm的无缝拼接，提高了拼接的精度。
为解决上述技术问题，本发明采用的技术方案是：一种基于物方反算的dem拼接方法，其特征在于，该方法包括以下步骤：
步骤一、获取图像照片：航摄相机采集n个测区的图像照片，并将得到的图像照片传输给计算机，n个所述测区中相邻两个所述测区的图像照片具有重叠区域，计算机将得到的n个所述测区的图像照片存入数据存储器中的图像信息数据库i中，其中，i＝{g1,g2,...,gh,...,gn}，gh表示第h个测区的图像照片，h＝1,2,...,n，n为不小于2的正整数；
步骤二、获取图像照片的物方坐标信息：计算机计算得到图像照片gh的像主点(xh,yh)和图像照片gh的dem，工作人员采集图像照片gh的地面影像信息(xhe-i,yhe-i,zhe-i)，其中，图像照片gh的dem包括图像照片gh的地形表面坐标(xhd-i,yhd-i,zhd-i)，i＝1,2,...,k，k为不小于2的正整数；
步骤三、根据图像照片的物方坐标信息反算图像照片的像点坐标：计算机根据图像照片gh的物方坐标信息反算图像照片gh的像点坐标(xhp-i,yhp-i)，
步骤四、根据图像照片的像点坐标进行dem裁切：
步骤401、筛选位于裁切边界内的图像照片的像点坐标：计算机根据裁切边界筛选图像照片gh的像点坐标(xhp-i,yhp-i)，裁切边界由xmin、xmax、ymax和ymin组成，其中，xmin表示横坐标最小裁切阈值，xmax表示横坐标最大裁切阈值，ymax表示纵坐标最大裁切阈值，ymin表示纵坐标最小裁切阈值，若
步骤402、构建裁切后的dem：计算机根据图像照片gh的保留像点坐标集qh中的像点坐标(xhp-i,yhp-i)在图像照片gh的地形表面坐标(xhd-i,yhd-i,zhd-i)的映射，得到图像照片gh的裁切后的dem；
步骤五、重复步骤二至步骤四，完成对图像信息数据库i中n个图像照片的dem的裁切；
步骤六、对裁切后的dem进行拼接：
步骤601、筛选相邻的两个图像照片的重复地形表面坐标：图像照片gl和图像照片gl+1表示图像信息数据库i中任意相邻的两个图像照片，l＝1,2,...,n-1，(x′ld-i,y′ld-i,z′ld-i)表示图像照片gl的裁切后的dem的地形表面坐标，(x′(l+1)d-i,y′(l+1)d-i,z′(l+1)d-i)表示图像照片gl+1的裁切后的dem的地形表面坐标，若
步骤602、计算相邻的两个图像照片的加权坐标：计算机根据公式
步骤603、计算机将图像照片gl的裁切后的dem的地形表面坐标(x′ld-i,y′ld-i,z′ld-i)和图像照片gl+1的裁切后的dem的地形表面坐标(x′(l+1)d-i,y′(l+1)d-i,z′(l+1)d-i)存储在图像照片gl和图像照片gl+1的独立地形表面坐标集
步骤604、获取拼接后的dem：计算机将步骤602中计算得到的图像照片gl和图像照片gl+1的加权坐标(x″ld-i,y″ld-i,z″ld-i)加入步骤603中图像照片gl和图像照片gl+1的独立地形表面坐标集
步骤605、重复步骤601至步骤604，完成对图像信息数据库i中任意相邻的两个图像照片的dem的拼接，得到拼接后的测区dem；
步骤七、显示拼接后的dem：计算机将步骤六中拼接后的测区dem存储在数据存储器中，并通过显示单元显示拼接后的测区dem。
上述的一种基于物方反算的dem拼接方法，其特征在于：步骤602中，wl和wl+1经公式
上述的一种基于物方反算的dem拼接方法，其特征在于：步骤一中相邻两个所述测区的图像照片具有20％～60％的重叠区域。
本发明还公开了一种基于物方反算的dsm拼接方法，其特征在于，该方法包括以下步骤：
步骤一、获取图像照片：航摄相机采集n个测区的图像照片，并将得到的图像照片传输给计算机，n个所述测区中相邻两个所述测区的图像照片具有重叠区域，计算机将得到的n个所述测区的图像照片存入数据存储器中的图像信息数据库i中，其中，i＝{g1,g2,...,gh,...,gn}，gh表示第h个测区的图像照片，h＝1,2,...,n，n为不小于2的正整数；
步骤二、获取图像照片的物方坐标信息：计算机计算得到图像照片gh的像主点(xh,yh)和图像照片gh的dsm，工作人员采集图像照片gh的地面影像信息(xhe-i,yhe-i,zhe-i)，其中，图像照片gh的dsm包括图像照片gh的地形表面坐标(xhd-i,yhd-i,zhd-i)和地表高程信息，图像照片gh的地表高程信息包括地表建筑物高度信息、地表桥梁高度信息和地表树木高度信息，i＝1,2,...,k，k为不小于2的正整数；
步骤三、根据图像照片的物方坐标信息反算图像照片的像点坐标：计算机根据图像照片gh的物方坐标信息反算图像照片gh的像点坐标(xhp-i,yhp-i)，
步骤四、根据图像照片的像点坐标进行dsm裁切：
步骤401、筛选位于裁切边界内的图像照片的像点坐标：计算机根据裁切边界筛选图像照片gh的像点坐标(xhp-i,yhp-i)，裁切边界由xmin、xmax、ymax和ymin组成，其中，xmin表示横坐标最小裁切阈值，xmax表示横坐标最大裁切阈值，ymax表示纵坐标最大裁切阈值，ymin表示纵坐标最小裁切阈值，若
步骤402、构建裁切后的dsm：计算机根据图像照片gh的保留像点坐标集qh中的像点坐标(xhp-i,yhp-i)在图像照片gh的地形表面坐标(xhd-i,yhd-i,zhd-i)和地表高程信息中的映射，得到图像照片gh的裁切后的dsm；
步骤五、重复步骤二至步骤四，完成对图像信息数据库i中n个图像照片的dsm的裁切；
步骤六、对裁切后的dsm进行拼接：
步骤601、筛选相邻的两个图像照片的重复地形表面坐标：图像照片gl和图像照片gl+1表示图像信息数据库i中任意相邻的两个图像照片，l＝1,2,...,n-1，(x′ld-i,y′ld-i,z′ld-i)表示图像照片gl的裁切后的dsm的地形表面坐标，(x′(l+1)d-i,y′(l+1)d-i,z′(l+1)d-i)表示图像照片gl+1的裁切后的dsm的地形表面坐标，若
步骤602、计算相邻的两个图像照片的加权坐标：计算机根据公式
步骤603、计算机将图像照片gl的裁切后的dsm的地形表面坐标(x′ld-i,y′ld-i,z′ld-i)和图像照片gl+1的裁切后的dsm的地形表面坐标(x′(l+1)d-i,y′(l+1)d-i,z′(l+1)d-i)存储在图像照片gl和图像照片gl+1的独立地形表面坐标集
步骤604、获取拼接后的dsm：计算机将步骤602中计算得到的图像照片gl和图像照片gl+1的加权坐标(x″ld-i,y″ld-i,z″ld-i)加入步骤603中图像照片gl和图像照片gl+1的独立地形表面坐标集
步骤605、重复步骤601至步骤604，完成对图像信息数据库i中任意相邻的两个图像照片的dsm的拼接，得到拼接后的测区dsm；
步骤七、显示拼接后的dsm：计算机将步骤六中拼接后的测区dsm存储在数据存储器中，并通过显示单元显示拼接后的测区dsm。
上述的一种基于物方反算的dsm拼接方法，其特征在于：步骤602中，wl和wl+1经公式
上述的一种基于物方反算的dsm拼接方法，其特征在于：步骤一中相邻两个所述测区的图像照片具有20％～60％的重叠区域。
本发明与现有技术相比具有以下优点：
1、本发明的步骤简单，设计合理，实现及使用操作方便。
2、本发明根据物方坐标信息反算图像照片的像点坐标，物方坐标信息包括图像照片的像主点、图像照片的地面影像信息和图像照片的地形表面坐标，由于dem和dsm均具有地形表面坐标，因此采用地形表面坐标参与计算图像照片的像点坐标，既适用于dem，又适用于dsm，实用效果好。
3、本发明根据裁切边界对像点坐标进行筛选，由位于裁切边界内的像点坐标在dem中的映射构建图像照片裁切后的dem，由位于裁切边界内的像点坐标在dsm中的映射构建图像照片裁切后的dsm，不需要人工干预，即可去掉重叠和错误数据，避免数据冗余。
4、本发明在进行dem的拼接时，筛选相邻的两个图像照片中具有同一水平坐标的地形表面坐标，对相邻的两个图像照片中具有同一水平坐标的地形表面坐标进行加权处理得到加权坐标，以加权坐标代替具有同一水平坐标的地形表面坐标，达到dem的无缝拼接，在进行dsm的拼接时，筛选相邻的两个图像照片中具有同一水平坐标的地形表面坐标，对相邻的两个图像照片中具有同一水平坐标的地形表面坐标进行加权处理得到加权坐标，以加权坐标代替具有同一水平坐标的地形表面坐标，达到dsm的无缝拼接，提高了dem和dsm的拼接精度。
综上所述，本发明步骤简单，设计合理，采用dem和dsm均具有的地形表面坐标参与计算图像照片的像点坐标，计算机根据裁切边界对像点坐标进行筛选，由位于裁切边界内的像点坐标在dem中的映射构建图像照片裁切后的dem，由位于裁切边界内的像点坐标在dsm中的映射构建图像照片裁切后的dsm，对相邻的两个图像照片中具有同一水平坐标的地形表面坐标进行加权处理得到加权坐标，以加权坐标代替具有同一水平坐标的地形表面坐标，达到dem和dsm的无缝拼接，提高了拼接的精度。
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明基于物方反算的dem拼接方法的方法流程图。
图2为本发明基于物方反算的dsm拼接方法的方法流程图。
图3为本发明的电路原理框图。
附图标记说明:
1—航摄相机；2—计算机；3—参数输入模块；
4—数据存储器；5—显示模块。
具体实施方式
实施例1
如图1和图3所示，本发明的一种基于物方反算的dem拼接方法，包括以下步骤：
步骤一、获取图像照片：航摄相机1采集n个测区的图像照片，并将得到的图像照片均传输给计算机2，n个所述测区中相邻两个所述测区的图像照片具有重叠区域，计算机2将得到的n个所述测区的图像照片存入数据存储器4中的图像信息数据库i中，其中，i＝{g1,g2,...,gh,...,gn}，gh表示第h个测区的图像照片，h＝1,2,...,n，n为不小于2的正整数。
需要说明的是，为了得到完整的测区图像照片，在获取图像照片时，相邻两个测区的图像照片具有20％～60％的重叠区域，图像照片对应的dem才能有3％～10％的重叠区域，以此避免测区遗漏。
步骤二、获取图像照片的物方坐标信息：计算机2计算得到图像照片gh的像主点(xh,yh)和图像照片gh的dem，工作人员采集图像照片gh的地面影像信息(xhe-i,yhe-i,zhe-i)，其中，图像照片gh的dem包括图像照片gh的地形表面坐标(xhd-i,yhd-i,zhd-i)，i＝1,2,...,k，k为不小于2的正整数。
图像照片gh的物方坐标信息由图像照片gh的像主点(xh,yh)、图像照片gh的地形表面坐标(xhd-i,yhd-i,zhd-i)和图像照片gh的地面影像信息(xhe-i,yhe-i,zhe-i)组成。
实际使用时，航摄相机1为数码相机，航摄相机1的最大分辨率为a×b，计算机2根据航摄相机1的最大分辨率计算图像照片gh的像主点(xh,yh)，
步骤三、根据图像照片的物方坐标信息反算图像照片的像点坐标：计算机2根据图像照片gh的物方坐标信息反算图像照片gh的像点坐标(xhp-i,yhp-i)，
需要说明的是，航向倾角
步骤四、根据图像照片的像点坐标进行dem裁切：
步骤401、筛选位于裁切边界内的图像照片的像点坐标：计算机2根据裁切边界筛选图像照片gh的像点坐标(xhp-i,yhp-i)，裁切边界由xmin、xmax、ymax和ymin组成，其中，xmin表示横坐标最小裁切阈值，xmax表示横坐标最大裁切阈值，ymax表示纵坐标最大裁切阈值，ymin表示纵坐标最小裁切阈值，若
需要说明的是，为了得到完整的测区图像照片，在获取图像照片时，相邻两个测区的图像照片具有20％～60％的重叠区域，才能保证相邻两个测区的图像照片的dem至少具有3％～10％以上的重叠数据，以此避免测区遗漏。但是由于相邻两个测区的图像照片的dem至少具有3％～10％以上的重叠数据，因此在进行dem拼接时，容易出现数据冗余的问题，所以需要对dem的边界信息进行裁切，去除重叠数据。
由横坐标最大裁切阈值xmax、横坐标最小裁切阈值xmin、纵坐标最大裁切阈值ymax和纵坐标最小裁切阈值ymin构成裁切边界，计算机2根据裁切边界对图像照片gh的像点坐标(xhp-i,yhp-i)进行筛选，若xmin≤xhp-i≤xmax且ymin≤yhp-i≤ymax，即认为图像照片gh的像点坐标(xhp-i,yhp-i)位于裁切边界内，此时，计算机2将图像照片gh的像点坐标(xhp-i,yhp-i)存储在数据存储器4中的图像照片gh的保留像点坐标集qh中。若xmin>xhp-i、xhp-i>xmax、ymin>yhp-i或yhp-i>ymax，即认为图像照片gh的像点坐标(xhp-i,yhp-i)超出裁切边界，此时计算机2对图像照片gh的像点坐标(xhp-i,yhp-i)赋予空值，即舍弃图像照片gh的像点坐标(xhp-i,yhp-i)。
步骤402、构建裁切后的dem：计算机2根据图像照片gh的保留像点坐标集qh中的像点坐标(xhp-i,yhp-i)在图像照片gh的地形表面坐标(xhd-i,yhd-i,zhd-i)的映射，得到图像照片gh的裁切后的dem；
由于图像照片gh的像点坐标(xhp-i,yhp-i)与图像照片gh的dem中的地形表面坐标(xhd-i,yhd-i,zhd-i)呈映射关系，因此对图像照片gh的像点坐标(xhp-i,yhp-i)的筛选，即相当于对图像照片gh的dem中的地形表面坐标(xhd-i,yhd-i,zhd-i)的筛选，筛选后的地形表面坐标(xhd-i,yhd-i,zhd-i)组成图像照片gh裁切后的dem。
步骤五、重复步骤二至步骤四，完成对图像信息数据库i中n个图像照片的dem的裁切。需要说明的是，由于步骤一中通过航摄获取了至少两张图像照片，因此步骤二至步骤四至少需要重复两次，完成对至少两张图像照片的dem的裁切。
步骤六、对裁切后的dem进行拼接：
步骤601、筛选相邻的两个图像照片的重复地形表面坐标：图像照片gl和图像照片gl+1表示图像信息数据库i中任意相邻的两个图像照片，l＝1,2,...,n-1，(x′ld-i,y′ld-i,z′ld-i)表示图像照片gl的裁切后的dem的地形表面坐标，(x′(l+1)d-i,y′(l+1)d-i,z′(l+1)d-i)表示图像照片gl+1的裁切后的dem的地形表面坐标，若
为了进行图像照片的dem拼接，图像照片的dem在完成裁切后，相邻图像照片的dem还存在至少3％的重叠数据。当(x′ld-i,y′ld-i,z′ld-i)和(x′(l+1)d-i,y′(l+1)d-i,z′(l+1)d-i)中的水平坐标相等时，即
步骤602、计算相邻的两个图像照片的加权坐标：计算机2根据公式
为了避免具有水平坐标的地形表面坐标带来的明显的拼接痕迹，对具有水平坐标的地形表面坐标进行加权处理，得到加权坐标(x″ld-i,y″ld-i,z″ld-i)。以加权坐标(x″ld-i,y″ld-i,z″ld-i)代替具有同一水平坐标的地形表面坐标(x′ld-i,y′ld-i,z′ld-i)和地形表面坐标(x′(l+1)d-i,y′(l+1)d-i,z′(l+1)d-i)，避免了拼接痕迹。
实际使用时，步骤602中，wl和wl+1经公式
步骤603、计算机2将图像照片gl的裁切后的dem的地形表面坐标(x′ld-i,y′ld-i,z′ld-i)和图像照片gl+1的裁切后的dem的地形表面坐标(x′(l+1)d-i,y′(l+1)d-i,z′(l+1)d-i)存储在图像照片gl和图像照片gl+1的独立地形表面坐标集
步骤604、获取拼接后的dem：计算机2将步骤602中计算得到的图像照片gl和图像照片gl+1的加权坐标(x″ld-i,y″ld-i,z″ld-i)加入步骤603中图像照片gl和图像照片gl+1的独立地形表面坐标集
由于dem只包括了地形的高程信息，并未包括其它地表信息，因此拼接地形表面坐标集
步骤605、重复步骤601至步骤604，完成对图像信息数据库i中任意相邻的两个图像照片的dem的拼接，得到拼接后的测区的dem。
步骤七、显示拼接后的dem：计算机2将步骤605中拼接后的测区dem存储在数据存储器4中，并通过显示单元5显示拼接后的测区dem。
实施例2
如图2和图3所示，本发明的一种基于物方反算的dsm拼接方法，包括以下步骤：
步骤一、获取图像照片：航摄相机1采集n个测区的图像照片，并将得到的图像照片均传输给计算机2，n个所述测区中相邻两个所述测区的图像照片具有重叠区域，计算机2将得到的n个所述测区的图像照片存入数据存储器4中的图像信息数据库i中，其中，i＝{g1,g2,...,gh,...,gn}，gh表示第h个测区的图像照片，h＝1,2,...,n，n为不小于2的正整数。
需要说明的是，为了得到完整的测区图像照片，在获取图像照片时，相邻两个测区的图像照片具有20％～60％的重叠区域，图像照片对应的dsm才能有3％～10％的重叠区域，以此避免测区遗漏。
步骤二、获取图像照片的物方坐标信息：计算机2计算得到图像照片gh的像主点(xh,yh)和图像照片gh的dsm，图像照片gh的dsm包括图像照片gh的地形表面坐标(xhd-i,yhd-i,zhd-i)和地表高程信息，地表高程信息包括地表建筑物高度信息、地表桥梁高度信息和地表树木高度信息，工作人员采集图像照片gh的地面影像信息(xhe-i,yhe-i,zhe-i)，i＝1,2,...,k，k为不小于2的正整数。
需要说明的是，dem和dsm均具有地形表面坐标(xhd-i,yhd-i,zhd-i)，因此采用地形表面坐标(xhd-i,yhd-i,zhd-i)参与计算像点坐标，既适用于dem，又适用于dsm，使用效果好。
步骤三、根据图像照片的物方坐标信息反算图像照片的像点坐标：计算机2根据图像照片gh的物方坐标信息反算图像照片gh的像点坐标(xhp-i,yhp-i)，
需要说明的是，航向倾角
步骤四、根据图像照片的像点坐标进行dsm裁切：
步骤401、筛选位于裁切边界内的图像照片的像点坐标：计算机2根据裁切边界筛选图像照片gh的像点坐标(xhp-i,yhp-i)，裁切边界由xmin、xmax、ymax和ymin组成，其中，xmin表示横坐标最小裁切阈值，xmax表示横坐标最大裁切阈值，ymax表示纵坐标最大裁切阈值，ymin表示纵坐标最小裁切阈值，若
由于相邻两个测区的图像照片的dsm至少具有3％～10％以上的重叠数据，因此在进行dsm拼接时，容易出现数据冗余的问题，所以需要对dsm的边界信息进行裁切，去除重叠数据。若xmin>xhp-i、xhp-i>xmax、ymin>yhp-i或yhp-i>ymax，即认为图像照片gh的像点坐标(xhp-i,yhp-i)超出裁切边界，此时计算机2对图像照片gh的像点坐标(xhp-i,yhp-i)赋予空值，即舍弃图像照片gh的像点坐标(xhp-i,yhp-i)。
步骤402，构建裁切后的dsm：计算机2根据图像照片gh的保留像点坐标集qh中的像点坐标(xhp-i,yhp-i)在图像照片gh的dsm的地形表面坐标(xhd-i,yhd-i,zhd-i)和地表高程信息中的映射，得到图像照片gh的裁切后的dsm。
由于图像照片gh的像点坐标(xhp-i,yhp-i)与图像照片gh的地形表面坐标(xhd-i,yhd-i,zhd-i)和地表高程信息呈映射关系，因此对图像照片gh的像点坐标(xhp-i,yhp-i)的筛选，即相当于对图像照片gh的地形表面坐标(xhd-i,yhd-i,zhd-i)和地表高程信息的筛选，筛选后的图像照片gh的地形表面坐标(xhd-i,yhd-i,zhd-i)和地表高程信息构成图像照片gh裁切后的dsm。
步骤五、重复步骤二至步骤四，完成对图像信息数据库i中n个图像照片的dem/dsm的裁切。
需要说明的是，由于步骤一中通过航摄获取了至少两张图像照片，因此步骤二至步骤四至少需要重复两次，完成对至少两张图像照片的dsm的裁切。
步骤六、对裁切后的dem/dsm进行拼接：
步骤601、筛选相邻的两个图像照片的重复地形表面坐标：图像照片gl和图像照片gl+1表示图像信息数据库i中任意相邻的两个图像照片，l＝1,2,...,n-1，(x′ld-i,y′ld-i,z′ld-i)表示图像照片gl的裁切后的dsm的地形表面坐标，(x′(l+1)d-i,y′(l+1)d-i,z′(l+1)d-i)表示图像照片gl+1的裁切后的dsm的地形表面坐标，若
为了进行图像照片的dsm拼接，图像照片的dsm在完成裁切后，相邻图像照片的dsm还存在至少3％的重叠数据。重叠数据会造成拼接痕迹，因此需要先筛选相邻的两个图像照片的重复地形表面坐标。
步骤602、计算相邻的两个图像照片的加权坐标：计算机2根据公式
为了避免具有水平坐标的地形表面坐标带来的明显的拼接痕迹，对具有水平坐标的地形表面坐标进行加权处理，得到加权坐标(x″ld-i,y″ld-i,z″ld-i)。以加权坐标(x″ld-i,y″ld-i,z″ld-i)代替具有同一水平坐标的地形表面坐标(x′ld-i,y′ld-i,z′ld-i)和地形表面坐标(x′(l+1)d-i,y′(l+1)d-i,z′(l+1)d-i)，因此在同一个水平坐标位置，只存在加权坐标(x″ld-i,y″ld-i,z″ld-i)的高程坐标，因此避免了拼接痕迹。
实际使用时，步骤602中，wl和wl+1经公式
步骤603、计算机2将图像照片gl的裁切后的dsm的地形表面坐标(x′ld-i,y′ld-i,z′ld-i)和图像照片gl+1的裁切后的dsm的地形表面坐标(x′(l+1)d-i,y′(l+1)d-i,z′(l+1)d-i)存储在图像照片gl和图像照片gl+1的独立地形表面坐标集
步骤604中，获取拼接后的dsm：计算机2根据拼接地形表面坐标集
由于dsm中的地形表面坐标(xhd-i,yhd-i,zhd-i)和地表高程信息相对应，因此根据拼接地形表面坐标集
步骤605、重复步骤601至步骤604，完成对图像信息数据库i中任意相邻的两个图像照片的dsm的拼接，得到拼接后的测区的dsm。
步骤七、拼接后的dsm显示：计算机2将步骤605中拼接后的测区dsm存储在数据存储器4中，并通过显示单元5显示拼接后的测区dsm。
以上所述，仅是本发明的实施例，并非对本发明作任何限制，凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化，均仍属于本发明技术方案的保护范围内。