【技术实现步骤摘要】
本专利技术属于雷达，具体涉及一种基于相位解缠的阵列干涉sar山区点云解模糊方法。
技术介绍
1、阵列干涉合成孔径雷达(array interferometric synthetic aperture radartomography，array-insar)作为干涉合成孔径雷达(syntheic aperture radarinterfemetry，insar)技术的延伸，通过在高程向不同位置采集数据，形成目标高程方向的分辨能力，能够有效解决传统insar面临的叠掩难题，在山区测绘、城区三维建模等方面具有重要的应用价值和潜力。array-insar在方位向通过合成孔径形成目标方位向的高分辨能力，在斜距向通过发射宽带雷达信号实现斜距方向的高分辨，在高程向通过阵列多通道技术形成获取高程向分辨能力。
2、阵列干涉sar系统高程向分辨能力来源于角度分辨，根据雷达分辨理论，分辨率和天线的尺寸成正比，阵列干涉sar系统为了保证获得足够高的高程向分辨能力，往往要求系统具有较大的阵列长度，同时考虑到系统设计成本、重量、体积等制约因素，阵列干涉sar系统的阵元之间往往是稀疏分布的，有限的阵元数目无法满足满阵的空间采样要求，从而导致高程模糊。系统允许观测的最大不模糊高程和系统阵元间距相关。不同于城区，山区地形变化剧烈，当地形高度变化超过系统的最大不模糊高程时，重建的山区点云发生高程模糊，该现象在阵列干涉sar山区三维重建中常常遇到。
3、针对上述问题，近些年来已有相关学者开展研究，按处理手段主要可分成以下几类：1、人工方式。利
技术实现思路
1、针对机载阵列干涉sar山区三维重建当中由于系统最大不模糊高程限制导致点云发生高程模糊问题，本专利技术提出一种基于相位解缠的阵列干涉sar山区点云解模糊方法。该方法通过直接对array-insar获取处理生成的山区点云进行操作，首先将山区三维点云结果进行矩阵化处理，然后借鉴干涉处理当中相位解缠技术，利用干涉相位解缠方法实现对山区点云高度模糊的求解。方法整体处理流程简单、易懂，可实现点云高程模糊的自动化求解。
2、本专利技术方法是通过下述的技术方案实现的：
3、一种基于相位解缠的阵列干涉sar山区点云解模糊方法，包括以下步骤：
4、步骤1、提取阵列干涉sar山区点云的点云高度矩阵；
5、步骤2、进行点云高度矩阵的滤波；
6、步骤3、进行阴影区检测和掩膜处理，获得掩膜处理后的点云高度图；
7、步骤4、对掩膜处理后的点云高度图利用基于枝切法的路径跟踪算法对点云高度结果进行解模糊，得到无模糊点云高度值；
8、步骤5、对无模糊点云高度值进行阴影区插值，获得阴影区点云高度；
9、步骤6、进行点云提取，得到每一个像素位置的点云的高程大小和点云强度值。
10、进一步地，所述步骤1中，输入阵列干涉sar通过对雷达回波数据进行处理后生成的山区三维点云，获得需要处理的阵列干涉sar山区模糊点云，点云坐标系为方位-斜距-高程坐标系；然后提取点云高度矩阵，机载的阵列干涉sar生成的点云结果用一大小为m×n×k的多维数组表示，其中m、n、k分别表示方位向、斜距向和高程向的像素数目，点云高度矩阵通过对每一方位-斜距单元提取所有高程向上幅度值最大的点云所在的高程位置形成，此时形成m×n大小的图像，该图像的每一个像素元素代表点云的高程值。
11、进一步地，所述步骤2中，采用干涉相位滤波方法，将点云高度矩阵转化为干涉相位，此时点云高度矩阵的元素(m,n)处的点云的高程值sm,n和相位p(m,n)之间的关系如下：
12、
13、其中，sm,n为点云的高程值，l为高程向的点云高度最大值；经过上式转化后，m×n大小的图像的像素元素的取值范围变为(-π,π]；采用圆周均值滤波方法对m×n大小的图像进行滤波处理，令滤波窗口为以待滤波的像素元素为中心、前后大小为(2m1+1)×(2n1+1)的矩形，m1、n1分别表示行和列方向上像素元素的数目，待滤波的像素元素为圆周均值滤波的输出值即：
14、
15、上式当中，为中间参数，arg{}表示取角操作，p和q分别表示对应矩阵的行和列，j为虚数单位，exp()表示指数函数。
16、进一步地，所述步骤3中，利用阵列干涉sar的首尾两通道图像，计算首尾两通道图像之间的相关系数图，利用相关系数图和阵列干涉sar幅度图像进行图像检测并提取阴影区，并利用形态学方法使阴影区的区域边缘规整化。
17、进一步地，所述步骤5包括对阴影区的点云高度进行插值处理，沿图像每一方位向，采用线性插值方法对存在空洞的图像位置进行插值，补偿图像的阴影区的点云的高度值，此时得到的图像中的每一个像素元素代表该像素元素的位置处的点云的高程值大小。
18、有益效果：
19、1)本专利技术方法可有效解决机载阵列干涉sar系统由于最大不模糊高程制约导致山区三维点云重建发生的高程模糊问题；
20、2)本专利技术方法和现有的基于手工方式的点云高度模糊求解手段相比，可提升点云高度模糊自动化求解的效率，减轻人力成本；和基于分割和重组的方法相比，本专利技术方法无需过多的参数调优过程，减少人工干预过程，同时对山区点云质量有更大容忍性；本专利技术方法和基于外源dem数据辅助的点云高度模糊求解方法相比，无需担心由于地形发生大的变化导致的点云高度模糊求解出现错误的问题；
21、3)本专利技术方法本质上将阵列干涉sar点云高度模糊求解问题归结为干涉处理当中的相位解缠过程。传统干涉相位解缠技术发展相对成熟，多种现有解缠算法可以在本专利技术方法中直接进行替代。
本文档来自技高网...【技术保护点】
1.一种基于相位解缠的阵列干涉SAR山区点云解模糊方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的一种基于相位解缠的阵列干涉SAR山区点云解模糊方法，其特征在于，所述步骤1中，输入阵列干涉SAR通过对雷达回波数据进行处理后生成的山区三维点云，获得需要处理的阵列干涉SAR山区模糊点云，点云坐标系为方位-斜距-高程坐标系；然后提取点云高度矩阵，机载的阵列干涉SAR生成的点云结果用一大小为M×N×K的多维数组表示，其中M、N、K分别表示方位向、斜距向和高程向的像素数目，点云高度矩阵通过对每一方位-斜距单元提取所有高程向上幅度值最大的点云所在的高程位置形成，此时形成M×N大小的图像，该图像的每一个像素元素代表点云的高程值。
3.根据权利要求2所述的一种基于相位解缠的阵列干涉SAR山区点云解模糊方法，其特征在于，所述步骤2中，采用干涉相位滤波方法，将点云高度矩阵转化为干涉相位，此时点云高度矩阵的元素(m,n)处的点云的高程值sm,n和相位P(m,n)之间的关系如下：
4.根据权利要求3所述的一种基于相位解缠的阵列干涉SAR山区点云解模糊方法，其
5.根据权利要求4所述的一种基于相位解缠的阵列干涉SAR山区点云解模糊方法，其特征在于，所述步骤5包括对阴影区的点云高度进行插值处理，沿图像每一方位向，采用线性插值方法对存在空洞的图像位置进行插值，补偿图像的阴影区的点云的高度值，此时得到的图像中的每一个像素元素代表该像素元素的位置处的点云的高程值大小。
...【技术特征摘要】
1.一种基于相位解缠的阵列干涉sar山区点云解模糊方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的一种基于相位解缠的阵列干涉sar山区点云解模糊方法，其特征在于，所述步骤1中，输入阵列干涉sar通过对雷达回波数据进行处理后生成的山区三维点云，获得需要处理的阵列干涉sar山区模糊点云，点云坐标系为方位-斜距-高程坐标系；然后提取点云高度矩阵，机载的阵列干涉sar生成的点云结果用一大小为m×n×k的多维数组表示，其中m、n、k分别表示方位向、斜距向和高程向的像素数目，点云高度矩阵通过对每一方位-斜距单元提取所有高程向上幅度值最大的点云所在的高程位置形成，此时形成m×n大小的图像，该图像的每一个像素元素代表点云的高程值。
3.根据权利要求2所述的一种基于相位解缠的阵列干涉sar山区点云解模糊方法，其特征在于，所述步骤2中...
【专利技术属性】
技术研发人员：郭其昌，万阳良，张新，
申请(专利权)人：中国科学院空天信息创新研究院，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。