一种基于立体测绘的阻水构筑物识别方法技术

本发明专利技术提出一种基于立体测绘的阻水构筑物识别方法，属于阻水构筑物识别技术领域。根据光学卫星立体影像、干涉雷达和倾斜摄影获取数据，基于地理场景生产方式，进行多类型数据协同，建立地形与纹理一体化的三维模型，实现河湖三维场景下阻水构筑物的可视化、定量化监测；具体包括以下步骤：S1.根据光学卫星立体影像量测阻水构筑物；S2.根据干涉雷达量测阻水构筑物；S3.基于倾斜摄影数据量测阻水构筑物。解决了现有技术中存在的正射影像上难识别围堤、采砂场、垃圾堆，现场巡逻范围有限、证据收集困难的技术问题。集困难的技术问题。集困难的技术问题。

【技术实现步骤摘要】
一种基于立体测绘的阻水构筑物识别方法


[0001]本申请涉及一种识别方法，尤其涉及一种基于立体测绘的阻水构筑物识别方法，属于阻水构筑物识别


技术介绍

[0002]近年来河湖管理很重要，常态化开展河湖“清四乱”专项整治行动，对以各种名义侵占河道、围垦湖泊、非法采砂、乱堆垃圾、乱建妨碍行洪建筑物等“四乱”突出问题全面执行河湖监管。
[0003]当前“清四乱”工作主要采用遥感影像监测配合车辆巡逻的方式，面临着在正射影像上难识别围堤、采砂场、垃圾堆，现场巡逻范围有限、证据收集困难等问题，急需新的科学技术和监管手段进行河湖监测和管理。

技术实现思路

[0004]在下文中给出了关于本专利技术的简要概述，以便提供关于本专利技术的某些方面的基本理解。应当理解，这个概述并不是关于本专利技术的穷举性概述。它并不是意图确定本专利技术的关键或重要部分，也不是意图限定本专利技术的范围。其目的仅仅是以简化的形式给出某些概念，以此作为稍后论述的更详细描述的前序。
[0005]鉴于此，为解决现有技术中存在的正射影像上难识别围堤、采砂场、垃圾堆，现场巡逻范围有限、证据收集困难的技术问题，本专利技术提供一种基于立体测绘的阻水构筑物识别方法。
[0006]方案一：一种基于立体测绘的阻水构筑物识别方法，根据光学卫星立体影像、干涉雷达和倾斜摄影获取数据，基于地理场景生产方式，进行数据协同，建立地形与纹理一体化的三维模型，实现河湖三维场景下阻水构筑物的可视化、定量化监测；具体包括以下步骤：
[0007]S1.根据光学卫星立体影像量测阻水构筑物；利用河道管理范围内光学立体测绘卫星的前视、后视、正视三线阵立体像对影像，生成数字地表模型，识别阻水构筑物，实现阻水构筑物测量监测；
[0008]S2.根据干涉雷达量测阻水构筑物；利用河道管理范围内的雷达数据，生成数字地表模型，识别阻水构筑物，实现阻水构筑物测量监测；
[0009]S3.基于倾斜摄影数据量测阻水构筑物；利用倾斜航空影像、像片控制点和检查点数据,生成实景三维模型，识别阻水构筑物，实现阻水构筑物测量监测。
[0010]优选的，S1具体方法是：
[0011]S11.立体模型构建：创建河道管理范围区工程文件，导入区域网平差成果，恢复立体模型；
[0012]S12.空中三角解算：从RPC文件中获得卫星扫描坐标系与像片坐标系间的变换关系及两景影像的相对位置关系，结合地面控制点或内业加密点对应的像点坐标、卫星影像匹配结果及影像的方位元素，构成区域网进行空中三角解算，得到卫星影像的精确外方位
元素，将影像纳入大地坐标系；
[0013]S13.水平核线影像生成：基于核线几何关系，确定立体像对的同名核线，生成沿核线方向排列的数字影像，消除影像对在垂直方向上的上下视差，经重采样形成规则格网核线影像；
[0014]S14.影像匹配：采集特征点、特征线、特征面要素，以特征要素作为约束，利用协方差函数方法识别两幅影像之间的同名点，完成影像匹配，协方差函数公式如下：
[0015]R(p，q)＝∫∫
(x，y)∈D
g(x，y)g(x+p，y+q)dxdy
[0016]若满足R(p0，q0)＞R(p，q)(p≠p0，q≠q0)条件，则p0，q0为搜索区影像相对于目标影像的位移参数；
[0017]S15.基于影像匹配获取的同名点数据，利用后方交会法和移动曲面拟合处理生成数字地表模型，根据数字地表模型数据模拟地形三维场景，实现阻水构筑物测量监测。
[0018]优选的，S2具体方法是：
[0019]S21.进行干涉配准处理：对SAR影像进行图像配准；
[0020]S22.进行干涉处理：对主辅影像的复数值进行重采样，并将主影像与配准采样后的辅影像进行共轭相乘，计算得到干涉图：
[0021][0022]Δr＝Bcos(θ
‑
θ
b
)
[0023]其中，λ为波长；f为雷达传播速度，即光速；B为两天线间基线距；θ为入射角；θ
b
为天线基线与飞行水平面法线间的夹角；
[0024]S23.进行多视处理：对SLC数据方位向和/或距离计算平均值，得到多视后的强度数据；
[0025]S24.进行滤波处理：使用Lee滤波消除由于信号的衰落引起的斑点和条纹噪声；
[0026]S25.进行相位解缠处理：对去平和滤波后的位相进行相位解缠，公式如下：
[0027][0028]其中，是解缠后的相位，N是条纹数目且为实数，n是条纹级次，是N的整数部分；
[0029]S26.进行相位转高程处理：进行轨道精炼和相位偏移的计算，消除斜坡相位，对卫星轨道和相位偏移进行纠正，再对解缠后的相位转化为高程或形变值；
[0030]S27.进行地理编码处理：将SAR数据从斜距坐标系转到地理坐标系，通过地理编码生成数字地表模型。
[0031]优选的，S3具体方法是：
[0032]S31.采集倾斜摄影数据、控制点和检测点；
[0033]S32.检查S31所述倾斜摄影数据和控制点，当倾斜摄影数据满足地面分辨率不大于5厘米时，斜摄影数据合格；控制点数据布设均匀、控制点测量精度满足平面精度不大于5厘米，高程精度不大于10厘米时，控制点数据合格，执行S33，如不满足上述条件，则数据不合格，返回S31；
[0034]S33.对S31所述倾斜摄影数据、控制点和检测点进行空三计算；
[0035]S34.检查进行空三计算后的倾斜摄影数据、控制点和检测点；若三维模型平面精度不大于20厘米和高程精度不大于30厘米，则数据合格，执行S35，否则数据不合格，返回S33；
[0036]S35.生成三维模型，并检查三维模型是否合格。
[0037]优选的，S31具体方法是：
[0038]S311.采集倾斜摄影数据；无人机外业航摄，基于测区内1:1万地形图、1:5万地形，确定测区内及周边有无高大密集建筑物、10公里内无军用、民用飞机场和大型设施影响飞行安全的因素，按照无人机航摄相关标准，开展实地踏勘；
[0039]S312.进行航空拍摄，对航空拍摄数据进行质量检查，数据质量满足航摄数据地面分辨率不大于5厘米、航向重叠度不小于80％、旁向重叠度不小于75％、航向超出摄区边界线不少于10条基线、旁向超出摄区边界线不少于3条航线，若质量不合格执行S311，若质量合格执行S313；
[0040]S313.根据航空拍摄数据进行内业成果整理和编制，形成倾斜摄影数据；
[0041]S314.控制点和检测点布设，以测区四角加中心为布设基础，平均分布于测区，控制点每平方公里不少于5个，检查点不少于2个；按照不规则区域网布点，利用已有遥感影像资料进行像控点的布设，通过测区内已有的HLJCORS进行像片控点的测量；
[0042]S315.对像片控制点的测量数据进行整理和检查，控制点与检查点布设需满足S314，测量本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于立体测绘的阻水构筑物识别方法，其特征在于，根据光学卫星立体影像、干涉雷达和倾斜摄影获取数据，基于地理场景生产方式，进行数据协同，建立地形与纹理一体化的三维模型，实现河湖三维场景下阻水构筑物的可视化、定量化监测；具体包括以下步骤：S1.根据光学卫星立体影像量测阻水构筑物；利用河道管理范围内光学立体测绘卫星的前视、后视、正视三线阵立体像对影像，生成数字地表模型，识别阻水构筑物，实现阻水构筑物测量监测；S2.根据干涉雷达量测阻水构筑物；利用河道管理范围内的雷达数据，生成数字地表模型，识别阻水构筑物，实现阻水构筑物测量监测；S3.基于倾斜摄影数据量测阻水构筑物；利用倾斜航空影像、像片控制点和检查点数据,生成实景三维模型，识别阻水构筑物，实现阻水构筑物测量监测。2.根据权利要求1所述的一种基于立体测绘的阻水构筑物识别方法，其特征在于，S1具体方法是：S11.立体模型构建：创建河道管理范围区工程文件，导入区域网平差成果，恢复立体模型；S12.空中三角解算：从RPC文件中获得卫星扫描坐标系与像片坐标系间的变换关系及两景影像的相对位置关系，结合地面控制点或内业加密点对应的像点坐标、卫星影像匹配结果及影像的方位元素，构成区域网进行空中三角解算，得到卫星影像的精确外方位元素，将影像纳入大地坐标系；S13.水平核线影像生成：基于核线几何关系，确定立体像对的同名核线，生成沿核线方向排列的数字影像，消除影像对在垂直方向上的上下视差，经重采样形成规则格网核线影像；S14.影像匹配：采集特征点、特征线、特征面要素，以特征要素作为约束，利用协方差函数方法识别两幅影像之间的同名点，完成影像匹配，协方差函数公式如下：Rp,q＝∫∫
(x,y)∈D
gx,yg
′
(x+p,y+q)dxdy若满足R(p0,q0)＞R(p,q)(p≠p0,q≠q0)条件，则p0，q0为搜索区影像相对于目标影像的位移参数；S15.基于影像匹配获取的同名点数据，利用后方交会法和移动曲面拟合处理生成数字地表模型，根据数字地表模型数据模拟地形三维场景，实现阻水构筑物测量监测。3.根据权利要求2所述的一种基于立体测绘的阻水构筑物识别方法，其特征在于，S2具体方法是：S21.进行干涉配准处理：对SAR影像进行图像配准；S22.进行干涉处理：对主辅影像的复数值进行重采样，并将主影像与配准采样后的辅影像进行共轭相乘，计算得到干涉图：Δr＝Bcosθ
‑
θ
b
)其中，λ为波长；f为雷达传播速度，即光速；B为两天线间基线距；θ为入射角；θ
b
为天线基线与飞行水平面法线间的夹角；
S23.进行多视处理：对SLC数据方位向和/或距离计算平均值，得到多视后的强度数据；S24.进行滤波处理：使用Lee滤波消除由于信号的衰落引起的斑点和条纹噪声；S25.进行相位解缠处理：对去平和滤波后的位相进行相位解缠，公式如下：其中，是解缠后的相位，N是条纹数目且为实数，n是条纹级次，是N的整数部分；S26.进行相位转高程处理：进行轨道精炼和相位偏移的计算，消除斜坡相位，对卫星轨道和相位偏移进行纠正，再对解缠后的相位转化为高程或形变值；S27.进行地理编码处理：将SAR数据从斜距坐标系转到地理坐标系，通过地理编码生成数字地表模型。4.根...

【专利技术属性】
技术研发人员：吴迪，杨爱玲，王俊，郭微微，张庆全，周艳，宋昊，杨阿龙，
申请(专利权)人：黑龙江省水利科学研究院，
类型：发明
国别省市：