基于数字高程模型快速更新的影像正射纠正方法及系统

摘要

本发明提供了一种基于数字高程模型快速更新的影像正射纠正方法及系统，属于地理测绘技术领域。所述方法，包括：根据当前数字表面模型以及历史数字表面模型，得到地表高差变化范围；根据与历史数字高程模型同时间获取的遥感影像以及当前遥感影像，得到影像特征变化范围；根据地表高差变化范围以及影像特征变化范围，将地面高程差大于第一高差阈值的范围确定为最终的更新区域，结合地面点云数据，得到更新区域的地面点云数据；根据更新区域的地面点云数据，得到待评估区域的数字高程模型的更新结果，根据更新结果进行未来某时刻的航空或卫星影像正射纠正；本发明提高了数字高程模型更新及影像正射纠正的效率。

权利要求

1.一种基于数字高程模型快速更新的影像正射纠正方法，其特征在于，包括以下过程：根据待评估区域的当前倾斜航空影像，得到密集匹配点云数据以及当前数字表面模型；

根据获取的密集匹配点云数据，得到待评估区域的地面点云数据；

根据当前数字表面模型以及历史数字表面模型，得到地表高差变化范围；

根据与历史数字高程模型同时间获取的遥感影像以及当前遥感影像，得到影像特征变化范围；

根据地表高差变化范围以及影像特征变化范围，将地面高程差大于第一高差阈值的范围确定为最终的更新区域，结合地面点云数据，得到更新区域的地面点云数据；

根据更新区域的地面点云数据，得到待评估区域的数字高程模型的更新结果，根据更新结果进行未来某时刻的航空或卫星影像正射纠正；

根据待评估区域的当前倾斜航空影像，得到密集匹配点云数据以及当前数字表面模型，包括：根据当前倾斜航空影像，进行倾斜影像空三加密后生成三维白模；

从多角度影像中筛选出最优角度影像，对三维白模进行纹理贴图，并对模型进行修饰，得到倾斜摄影三维模型；

根据倾斜摄影三维模型，得到密集匹配点云数据以及当前数字表面模型。

2.如权利要求1所述的基于数字高程模型快速更新的影像正射纠正方法，其特征在于，对模型进行修饰，包括：不完整建筑物的修补，悬浮、破面和结构性下沉的删除，水面异常凸起或下沉的压平，以及水域纹理痕迹的修改。

3.如权利要求1‑2任一项所述的基于数字高程模型快速更新的影像正射纠正方法，其特征在于，根据当前数字表面模型以及历史数字表面模型，得到地表高差变化范围，包括：根据当前数字表面模型以及历史数字表面模型的对比，得到高差大于第二高差阈值的所有区域，判断各个高差大于第二高差阈值的区域是否小于设定面积阈值；

如是，则判定为噪声；否则，合并相邻的大于面积阈值的区域，以所有大于设定面积阈值的区域为地表高差变化范围。

4.如权利要求1‑2任一项所述的基于数字高程模型快速更新的影像正射纠正方法，其特征在于，得到更新区域的地面点云数据，还包括如下过程：

将更新区域内旧的已有点云删除，替换为新的更新区域的地面点云数据，利用更新后点云数据的地面点，内插后得到规定格网间距的更新后的数字高程模型。

5.一种基于数字高程模型快速更新的影像正射纠正系统，其特征在于，包括：倾斜摄影三维模型生成模块，被配置为：根据待评估区域的当前倾斜航空影像，得到密集匹配点云数据以及当前数字表面模型；

地面点云提取模块，被配置为：根据获取的密集匹配点云数据，得到待评估区域的地面点云数据；

地表高差变化范围生成模块，被配置为：根据当前数字表面模型以及历史数字表面模型，得到地表高差变化范围；

影像特征变化范围生成模块，被配置为：根据与历史数字高程模型同时间获取的遥感影像以及当前遥感影像，得到影像特征变化范围；

更新区域点云更新模块，被配置为：根据地表高差变化范围以及影像特征变化范围，将地面高程差大于第一高差阈值的范围确定为最终的更新区域，结合地面点云数据，得到更新区域的地面点云数据；

影像正射纠正模块，被配置为：根据更新区域的地面点云数据，得到待评估区域的数字高程模型的更新结果，根据更新结果进行未来某时刻的航空或卫星影像正射纠正；

倾斜摄影三维模型生成模块中，包括：

根据当前倾斜航空影像，进行倾斜影像空三加密后生成三维白模；

从多角度影像中筛选出最优角度影像，对三维白模进行纹理贴图，并对模型进行修饰，得到倾斜摄影三维模型；

根据倾斜摄影三维模型，得到密集匹配点云数据以及当前数字表面模型。

6.如权利要求5所述的基于数字高程模型快速更新的影像正射纠正系统，其特征在于，对模型进行修饰，包括：不完整建筑物的修补，悬浮、破面和结构性下沉的删除，水面异常凸起或下沉的压平，以及水域纹理痕迹的修改。

7.如权利要求5‑6任一项所述的基于数字高程模型快速更新的影像正射纠正系统，其特征在于，根据当前数字表面模型以及历史数字表面模型，得到地表高差变化范围，包括：根据当前数字表面模型以及历史数字表面模型的对比，得到高差大于第二高差阈值的所有区域，判断各个高差大于第二高差阈值的区域是否小于设定面积阈值；

如是，则判定为噪声；否则，合并相邻的大于面积阈值的区域，以所有大于设定面积阈值的区域为地表高差变化范围。

8.如权利要求5‑6任一项所述的基于数字高程模型快速更新的影像正射纠正系统，其特征在于，更新区域点云更新模块中，得到更新区域的地面点云数据，还包括如下过程：将更新区域内旧的已有点云删除，替换为新的更新区域的地面点云数据，利用更新后点云数据的地面点，插值构建规定格网间距的更新后的数字高程模型。

说明书

基于数字高程模型快速更新的影像正射纠正方法及系统

技术领域

[0001] 本发明涉及地理测绘技术领域，特别涉及一种基于数字高程模型快速更新的影像正射纠正方法及系统。

背景技术

[0002] 本部分的陈述仅仅是提供了与本发明相关的背景技术，并不必然构成现有技术。

[0003] 利用已有DEM（Digital Elevation Model，数字高程模型）成果进行DOM（Digital Orthophoto Map，数字正射影像图）制作、地形测量以及相关行业应用等，可大大提高效率，节省成本，但是由于DEM时效性不够，会导致DOM、地形以及相关成果精度不能满足要求，无法实现快速的影像正射纠正处理。

[0004] 目前主要的DEM数据生成及更新方法有：传统摄影测量方法、机载激光LiDAR点云生产DEM以及外业实地补测等。传统摄影测量方法成本高、作业效率低，生产周期长；利用机载激光LiDAR点云生成DEM，存在航飞成本高、受空域、天气等影响大、点云分类工作量大等问题；而外业补测存在工作量大和经费高的问题。

发明内容

[0005] 为了解决现有技术的不足，本发明提供了一种基于数字高程模型快速更新的影像正射纠正方法及系统，解决了利用现有方法更新小面积区域数字高程模型时重新航飞成本高、生产周期长和外业实地测量工作量大的问题，实现了地理环境数字高程模型的快速更新，提高了航空或卫星影像正射纠正的效率。

[0006] 为了实现上述目的，本发明采用如下技术方案：

[0007] 本发明第一方面提供了一种基于数字高程模型快速更新的影像正射纠正方法。

[0008] 一种基于数字高程模型快速更新的影像正射纠正方法，包括以下过程：

[0009] 根据待评估区域的当前倾斜航空影像，得到密集匹配点云数据以及当前数字表面模型；

[0010] 根据获取的密集匹配点云数据，得到待评估区域的地面点云数据；

[0011] 根据当前数字表面模型以及历史数字表面模型，得到地表高差变化范围；

[0012] 根据与历史数字高程模型同时间获取的遥感影像以及当前遥感影像，得到影像特征变化范围；

[0013] 根据地表高差变化范围以及影像特征变化范围，将地面高程差大于第一高差阈值的范围确定为最终的更新区域，结合地面点云数据，得到更新区域的地面点云数据；

[0014] 根据更新区域的地面点云数据，得到待评估区域的数字高程模型的更新结果，根据更新结果进行未来某时刻的航空或卫星影像正射纠正。

[0015] 作为本发明第一方面进一步的限定，根据待评估区域的当前倾斜航空影像，得到密集匹配点云数据以及当前数字表面模型，包括：

[0016] 根据当前倾斜航空影像，进行倾斜影像空三加密后生成三维白模；

[0017] 从多角度影像中筛选出最优角度影像，对三维白模进行纹理贴图，并对模型进行修饰，得到倾斜摄影三维模型；

[0018] 根据倾斜摄影三维模型，得到密集匹配点云数据以及当前数字表面模型。

[0019] 作为本发明第一方面更进一步的限定，对模型进行修饰，包括：

[0020] 不完整建筑物的修补，悬浮、破面和结构性下沉的删除，水面异常凸起或下沉的压平，以及水域纹理痕迹的修改。

[0021] 作为本发明第一方面进一步的限定，根据当前数字表面模型以及历史数字表面模型，得到地表高差变化范围，包括：

[0022] 根据当前数字表面模型以及历史数字表面模型的对比，得到高差大于第二高差阈值的所有区域，判断各个高差大于第二高差阈值的区域是否小于设定面积阈值；

[0023] 如是，则判定为噪声；否则，合并相邻的大于面积阈值的区域，以所有大于设定面积阈值的区域为地表高差变化范围。

[0024] 作为本发明第一方面进一步的限定，得到更新区域的地面点云数据，还包括如下过程：

[0025] 将更新区域内旧的已有点云删除，替换为新的更新区域的地面点云数据，利用更新后点云数据的地面点，内插后得到规定格网间距的更新后的数字高程模型。

[0026] 本发明第二方面提供了一种基于数字高程模型快速更新的影像正射纠正系统。

[0027] 一种基于数字高程模型快速更新的影像正射纠正系统，包括：

[0028] 倾斜摄影三维模型生成模块，被配置为：根据待评估区域的当前倾斜航空影像，得到密集匹配点云数据以及当前数字表面模型；

[0029] 地面点云提取模块，被配置为：根据获取的密集匹配点云数据，得到待评估区域的地面点云数据；

[0030] 地表高差变化范围生成模块，被配置为：根据当前数字表面模型以及历史数字表面模型，得到地表高差变化范围；

[0031] 影像特征变化范围生成模块，被配置为：根据与历史数字高程模型同时间获取的遥感影像以及当前遥感影像，得到影像特征变化范围；

[0032] 更新区域点云更新模块，被配置为：根据地表高差变化范围以及影像特征变化范围，将地面高程差大于第一高差阈值的范围确定为最终的更新区域，结合地面点云数据，得到更新区域的地面点云数据；

[0033] 影像正射纠正模块，被配置为：根据更新区域的地面点云数据，得到待评估区域的数字高程模型的更新结果，根据更新结果进行未来某时刻的航空或卫星影像正射纠正。

[0034] 作为本发明第二方面进一步的限定，倾斜摄影三维模型生成模块中，包括：

[0035] 根据当前倾斜航空影像，进行倾斜影像空三加密后生成三维白模；

[0036] 从多角度影像中筛选出最优角度影像，对三维白模进行纹理贴图，并对模型进行修饰，得到倾斜摄影三维模型；

[0037] 根据倾斜摄影三维模型，得到密集匹配点云数据以及当前数字表面模型。

[0038] 作为本发明第二方面更进一步的限定，对模型进行修饰，包括：

[0039] 不完整建筑物的修补，悬浮、破面和结构性下沉的删除，水面异常凸起或下沉的压平，以及水域纹理痕迹的修改。

[0040] 作为本发明第二方面进一步的限定，根据当前数字表面模型以及历史数字表面模型，得到地表高差变化范围，包括：

[0041] 根据当前数字表面模型以及历史数字表面模型的对比，得到高差大于第二高差阈值的所有区域，判断各个高差大于第二高差阈值的区域是否小于设定面积阈值；

[0042] 如是，则判定为噪声；否则，合并相邻的大于面积阈值的区域，以所有大于设定面积阈值的区域为地表高差变化范围。

[0043] 作为本发明第二方面进一步的限定，更新区域点云更新模块中，得到更新区域的地面点云数据，还包括如下过程：

[0044] 将更新区域内旧的已有点云删除，替换为新的更新区域的地面点云数据，利用更新后点云数据的地面点，插值构建规定格网间距的更新后的数字高程模型。

[0045] 与现有技术相比，本发明的有益效果是：

[0046] 1、本发明创新性的提出了一种基于数字高程模型快速更新的影像正射纠正方法及系统，解决了利用现有方法更新小面积区域数字高程模型时重新航飞成本高、生产周期长和外业实地测量工作量大的问题，实现了地理环境数字高程模型的快速更新，提高了影像正射纠正的效率。

[0047] 2、本发明创新性的提出了一种基于数字高程模型快速更新的影像正射纠正方法及系统，利用倾斜摄影三维模型和遥感影像进行数字高程模型的更新，是最新三维成果的再利用，极大的节省了外业航飞或者人工采集原始数据的时间、人力和经费等一些列生产成本。

[0048] 3、本发明创新性的提出了一种基于数字高程模型快速更新的影像正射纠正方法及系统，结合两期数字表面模型与遥感影像进行人机交互变化检测，能快速准确地发现更新区域，极大地提高了生产效率，缩短了生产工期。

[0049] 本发明附加方面的优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。

附图说明

[0050] 构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。

[0051] 图1为本发明实施例1提供的基于数字高程模型快速更新的影像正射纠正方法的流程示意图；

[0052] 图2为本发明实施例1提供的历史遥感影像示意图；

[0053] 图3为本发明实施例1提供的当前遥感影像示意图。

具体实施方式

[0054] 下面结合附图与实施例对本发明作进一步说明。

[0055] 在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。

[0056] 实施例1：

[0057] 如图1所示，本发明实施例1提供了一种基于数字高程模型快速更新的影像正射纠正方法，包括以下过程：

[0058] S1：倾斜摄影三维模型成果制作。

[0059] 采用倾斜摄影自动化建模技术，使用三维建模软件（例如Mirauge3D、Smart3D、Context Capture Center、Photoscan等，本领域技术人员可选用获得使用授权的三维建模软件进行作业）；

[0060] 利用最新获取的倾斜航空影像，完成倾斜影像空三加密，通过三维建模软件（例如Mirauge3D、Smart3D、Context Capture Center、Photoscan等，本领域技术人员可选用获得使用授权的三维建模软件进行作业）自动生产三维白模，并从多角度影像中筛选出最优角度影像，对三维白模进行纹理贴图，并对模型进行修饰；

[0061] 本实施例所述的修饰，包括：不完整建筑物的修补，悬浮、破面、结构性下沉的删除，水面异常凸起或下沉的压平，以及水域纹理痕迹的修改等等，从而获得三维模型成果。

[0062] S2：输出密集匹配点云和数字表面模型（DSM）成果。

[0063] 基于最新高分辨率（优于0.15米）倾斜摄影三维模型成果，使用三维建模软件（例如模方、LiDARPro等，本领域技术人员可选用获得使用授权的三维建模软件进行作业），输出测区内1米间距（可根据需要设置）的密集匹配点云数据以及1米格网数字表面模型（DSM）成果。

[0064] S3：密集匹配点云处理。

[0065] 对测区内的密集匹配点云，利用点云数据处理软件（例如TerraSolid、LiDAR DP、LiDARPro等，本领域技术人员可选用获得使用授权的点云数据处理软件进行作业），根据机器性能等实际情况设置分块大小（例如1:10000比例尺标准分幅）进行自动分块，使用自动分类宏命令，试验合适的分类参数，对密集匹配点云进行自动分类，提取地面点云。

[0066] S4：更新区域确定。

[0067] S4.1：对比DSM发现地表高差变化处。

[0068] 城市的变化发展，大量建（构）筑物的建造，使得地形地表发生变化，根据一般建筑物变化的高度与面积大小，设置高差和面积阈值，对比最新三维DSM成果与旧的已有点云DSM成果，先设置地表高差阈值（即第一高差阈值，例如2.5米），认为新旧两期DSM的高度差大于地表高差阈值的地方地形发生了变化；再设置面积阈值（例如200平方米），小于面积阈值的认为是噪声，予以舍弃，然后合并相邻的大于面积阈值的变化范围，得到地表高差变化范围。

[0069] S4.2：遥感影像变化检测。

[0070] 如图2和图3所示，将旧的点云同期遥感影像（即历史点云同期遥感影像）和最新遥感影像（即当前遥感影像）导入遥感影像处理和分析软件（例如ENVI、Erdas等），使用多尺度分割算法，试验合适的分割参数，进行影像分割，再采用合适的分类方法（例如可以采用现有的监督分类或非监督分类算法）进行自动分类，然后检查分类结果，剔除噪声改正错误，根据分类结果进行变化检测，得到的变化检测结果为影像特征变化范围。

[0071] S4.3：结果对比分析检查。

[0072] 综合DSM上获取的地表高差变化范围和遥感影像上获取的特征变化范围，对所有变化范围内的新旧两期地面高程进行比较，分析检查、修正变化范围，将地面高差大于设定值（即第二高差阈值，例如0.5米，根据1：2000比例尺的DEM精度要求确定，这里的DEM为Digital Elevation Model的简写，即为数字高程模型）的范围确定为最终的更新区域。

[0073] S5：更新区域点云分类编辑。

[0074] 利用点云处理软件（例如TerraSolid、LiDAR DP、LiDARPro等，本领域技术人员可选用获得使用授权的点云数据处理软件进行作业），自动批处理裁出更新区域内自动分类后的密集匹配点云，目视检查分类效果，编辑修改错误分类点，得到更新区域内准确的地面点云。

[0075] S6：点云镶嵌。

[0076] 将更新区域内旧的已有点云成果删除，替换为分类编辑后的密集匹配点云，保证两套点云无缝拼接，更新为最新点云数据。

[0077] S7：精度质量检查。

[0078] 对更新后点云数据，采用人机交互方式进行分类精度和质量检查，具体的，通过将点云按类别、高程显示等方法，目视检查更新后点云；对有疑问处用断面图进行查询、分析。

[0079] S8：更新后DEM成果输出。

[0080] 利用更新后点云数据的地面点，内插输出规定格网间距的DEM成果，更具体的，将更新后点云中的所有地面点作为权值相同的特征点，导入点云处理软件（例如TerraSolid、LIDAR DP等获得使用授权的点云处理软件）中，采用不规则三角网内插方法，根据数据实际情况设置最大构网距离，进行DEM构建，制作规定格网间距的DEM成果。

[0081] S9：根据更新后的DEM进行新的航空或卫星影像的正射纠正。

[0082] 这里的根据DEM数据进行航空或卫星影像的正射纠正策略为现有方案的技术，例如可以采用现有的INPHO、PhotoScan、CIPS集群、GXL等获得使用授权的处理软件，直接的将影像和DEM数据导入，即可得到正射纠正后的影像。

[0083] 实施例2：

[0084] 本发明实施例2提供了一种基于数字高程模型快速更新的影像正射纠正系统，包括：

[0085] 倾斜摄影三维模型生成模块，被配置为：根据待评估区域的当前倾斜航空影像，得到密集匹配点云数据以及当前数字表面模型；

[0086] 地面点云提取模块，被配置为：根据获取的密集匹配点云数据，得到待评估区域的地面点云数据；

[0087] 地表高差变化范围生成模块，被配置为：根据当前数字表面模型以及历史数字表面模型，得到地表高差变化范围；

[0088] 影像特征变化范围生成模块，被配置为：根据与历史数字高程模型同时间获取的遥感影像以及当前遥感影像，得到影像特征变化范围；

[0089] 更新区域点云更新模块，被配置为：根据地表高差变化范围以及影像特征变化范围，将地面高程差大于第一高差阈值的范围确定为最终的更新区域，结合地面点云数据，得到更新区域的地面点云数据；

[0090] 影像正射纠正模块，被配置为：根据更新区域的地面点云数据，得到待评估区域的数字高程模型的更新结果，根据更新结果进行未来某时刻的影像正射纠正。

[0091] （1）倾斜摄影三维模型生成模块，具体的，包括：

[0092] 采用倾斜摄影自动化建模技术，使用三维建模软件；

[0093] 利用最新获取的倾斜航空影像，完成倾斜影像空三加密，通过软件自动生产三维白模，并从多角度影像中筛选出最优角度影像，对三维白模进行纹理贴图，并对模型进行修饰；

[0094] 本实施例所述的修饰，包括：不完整建筑物的修补，悬浮、破面、结构性下沉的删除，水面异常凸起或下沉的压平，以及水域纹理痕迹的修改等等，从而获得三维模型成果；

[0095] 基于最新高分辨率（优于0.15米）倾斜摄影三维模型成果，使用三维建模软件，输出测区内1米间距（可根据需要设置）的密集点云数据以及1米格网数字表面模型（DSM）成果。

[0096] （2）地面点云提取模块，具体的，包括：

[0097] 对测区内的密集匹配点云，利用点云数据处理软件，根据机器性能等实际情况设置分块大小（例如1:10000比例尺标准分幅）进行自动分块，使用自动分类宏命令，试验合适的分类参数，对密集匹配点云进行自动分类，提取地面点云。

[0098] （3）地表高差变化范围生成模块，具体的，包括：

[0099] 城市的变化发展，大量建（构）筑物的建造，使得地形地表发生变化，根据一般建筑物变化的高度与面积大小，设置高差和面积阈值，对比最新三维DSM成果与旧的已有点云DSM成果，先设置地表高差阈值（即第一高差阈值，例如2.5米），认为新旧两期DSM的高度差大于地表高差阈值的地方地形发生了变化；再设置面积阈值（例如200平方米），小于面积阈值的认为是噪声，予以舍弃，然后合并相邻的大于面积阈值的变化范围，得到地表高差变化范围。

[0100] （4）影像特征变化范围生成模块，具体的，包括：

[0101] 将旧的点云同期遥感影像（即历史点云同期遥感影像）和最新遥感影像（即当前遥感影像）导入遥感影像处理和分析软件，使用多尺度分割算法，试验合适的分割参数，进行影像分割，再采用合适的分类方法（例如可以采用现有的监督分类或非监督分类算法）进行自动分类，然后检查分类结果，剔除噪声改正错误，根据分类结果进行变化检测，得到的变化检测结果为影像特征变化范围。

[0102] 所述系统的其他工作方法与实施例1所述的基于数字高程模型快速更新的影像正射纠正方法相同，这里不再赘述。

[0103] 以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。