CN117872366B - 一种基于相位解缠的阵列干涉sar山区点云解模糊方法 - Google Patents
一种基于相位解缠的阵列干涉sar山区点云解模糊方法 Download PDFInfo
- Publication number
- CN117872366B CN117872366B CN202311621515.9A CN202311621515A CN117872366B CN 117872366 B CN117872366 B CN 117872366B CN 202311621515 A CN202311621515 A CN 202311621515A CN 117872366 B CN117872366 B CN 117872366B
- Authority
- CN
- China
- Prior art keywords
- point cloud
- height
- image
- elevation
- array
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 59
- 238000012545 processing Methods 0.000 claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims description 29
- 238000001914 filtration Methods 0.000 claims description 24
- 238000001514 detection method Methods 0.000 claims description 6
- ODKSFYDXXFIFQN-UHFFFAOYSA-M argininate Chemical compound [O-]C(=O)C(N)CCCNC(N)=N ODKSFYDXXFIFQN-UHFFFAOYSA-M 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 230000000877 morphologic effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims 1
- 230000000875 corresponding effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000000873 masking effect Effects 0.000 description 2
- 230000006798 recombination Effects 0.000 description 2
- 238000005215 recombination Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/882—Radar or analogous systems specially adapted for specific applications for altimeters
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提出一种基于相位解缠的阵列干涉SAR山区点云解模糊方法，属于雷达技术领域，通过直接对Array‑InSAR处理生成的山区点云数据进行操作，首先将山区三维点云结果进行矩阵化处理，然后借鉴干涉处理当中相位解缠技术，利用干涉相位解缠技术实现对山区点云高度模糊的求解。本发明流程简单、易懂，可实现点云高程模糊的自动化求解。
Description
技术领域
本发明属于雷达技术领域，具体涉及一种基于相位解缠的阵列干涉SAR山区点云解模糊方法。
背景技术
阵列干涉合成孔径雷达(Array Interferometric Synthetic Aperture RadarTomography，Array-InSAR)作为干涉合成孔径雷达(Syntheic Aperture RadarInterfemetry，InSAR)技术的延伸，通过在高程向不同位置采集数据，形成目标高程方向的分辨能力，能够有效解决传统InSAR面临的叠掩难题，在山区测绘、城区三维建模等方面具有重要的应用价值和潜力。Array-InSAR在方位向通过合成孔径形成目标方位向的高分辨能力，在斜距向通过发射宽带雷达信号实现斜距方向的高分辨，在高程向通过阵列多通道技术形成获取高程向分辨能力。
阵列干涉SAR系统高程向分辨能力来源于角度分辨，根据雷达分辨理论，分辨率和天线的尺寸成正比，阵列干涉SAR系统为了保证获得足够高的高程向分辨能力，往往要求系统具有较大的阵列长度，同时考虑到系统设计成本、重量、体积等制约因素，阵列干涉SAR系统的阵元之间往往是稀疏分布的，有限的阵元数目无法满足满阵的空间采样要求，从而导致高程模糊。系统允许观测的最大不模糊高程和系统阵元间距相关。不同于城区，山区地形变化剧烈，当地形高度变化超过系统的最大不模糊高程时，重建的山区点云发生高程模糊，该现象在阵列干涉SAR山区三维重建中常常遇到。
针对上述问题，近些年来已有相关学者开展研究，按处理手段主要可分成以下几类：1、人工方式。利用人工手段对山区的点云进行手动解模糊处理；2、基于分割和重组的山区点云解模糊方法。该方法首先对模糊点云进行分类处理，然后将分类后的点云结果沿高程向重新排布从而实现点云解模糊。3、利用外源DEM数据辅助的山区点云解模糊方法。通过已知高程信息的DEM数据，利用匹配的方法进行点云高程模糊数的估计，从而完成点云解模糊。分析上述方法，其中方法1需要沿方位向逐步进行处理，处理效率低下；方法2依赖点云质量，点云质量较差时导致分类出现误差，进而影响点云解模糊效果；方法3依赖外源DEM数据的精度，当发生较大的地形变化时，该方法可能出现错误结果。
发明内容
针对机载阵列干涉SAR山区三维重建当中由于系统最大不模糊高程限制导致点云发生高程模糊问题，本发明提出一种基于相位解缠的阵列干涉SAR山区点云解模糊方法。该方法通过直接对Array-InSAR获取处理生成的山区点云进行操作，首先将山区三维点云结果进行矩阵化处理，然后借鉴干涉处理当中相位解缠技术，利用干涉相位解缠方法实现对山区点云高度模糊的求解。方法整体处理流程简单、易懂，可实现点云高程模糊的自动化求解。
本发明方法是通过下述的技术方案实现的：
一种基于相位解缠的阵列干涉SAR山区点云解模糊方法，包括以下步骤：
步骤1、提取阵列干涉SAR山区点云的点云高度矩阵；
步骤2、进行点云高度矩阵的滤波；
步骤3、进行阴影区检测和掩膜处理，获得掩膜处理后的点云高度图；
步骤4、对掩膜处理后的点云高度图利用基于枝切法的路径跟踪算法对点云高度结果进行解模糊，得到无模糊点云高度值；
步骤5、对无模糊点云高度值进行阴影区插值，获得阴影区点云高度；
步骤6、进行点云提取，得到每一个像素位置的点云的高程大小和点云强度值。
进一步地，所述步骤1中，输入阵列干涉SAR通过对雷达回波数据进行处理后生成的山区三维点云，获得需要处理的阵列干涉SAR山区模糊点云，点云坐标系为方位-斜距-高程坐标系；然后提取点云高度矩阵，机载的阵列干涉SAR生成的点云结果用一大小为M×N×K的多维数组表示，其中M、N、K分别表示方位向、斜距向和高程向的像素数目，点云高度矩阵通过对每一方位-斜距单元提取所有高程向上幅度值最大的点云所在的高程位置形成，此时形成M×N大小的图像，该图像的每一个像素元素代表点云的高程值。
进一步地，所述步骤2中，采用干涉相位滤波方法，将点云高度矩阵转化为干涉相位，此时点云高度矩阵的元素(m,n)处的点云的高程值sm,n和相位P(m,n)之间的关系如下：
其中，sm,n为点云的高程值，L为高程向的点云高度最大值；经过上式转化后，M×N大小的图像的像素元素的取值范围变为(-π,π]；采用圆周均值滤波方法对M×N大小的图像进行滤波处理，令滤波窗口为以待滤波的像素元素为中心、前后大小为(2M1+1)×(2N1+1)的矩形，M1、N1分别表示行和列方向上像素元素的数目，待滤波的像素元素为圆周均值滤波的输出值即：
上式当中，为中间参数，arg{}表示取角操作，p和q分别表示对应矩阵的行和列，j为虚数单位，exp()表示指数函数。
进一步地，所述步骤3中，利用阵列干涉SAR的首尾两通道图像，计算首尾两通道图像之间的相关系数图，利用相关系数图和阵列干涉SAR幅度图像进行图像检测并提取阴影区，并利用形态学方法使阴影区的区域边缘规整化。
进一步地，所述步骤5包括对阴影区的点云高度进行插值处理，沿图像每一方位向，采用线性插值方法对存在空洞的图像位置进行插值，补偿图像的阴影区的点云的高度值，此时得到的图像中的每一个像素元素代表该像素元素的位置处的点云的高程值大小。
有益效果：
1)本发明方法可有效解决机载阵列干涉SAR系统由于最大不模糊高程制约导致山区三维点云重建发生的高程模糊问题；
2)本发明方法和现有的基于手工方式的点云高度模糊求解手段相比，可提升点云高度模糊自动化求解的效率，减轻人力成本；和基于分割和重组的方法相比，本发明方法无需过多的参数调优过程，减少人工干预过程，同时对山区点云质量有更大容忍性；本发明方法和基于外源DEM数据辅助的点云高度模糊求解方法相比，无需担心由于地形发生大的变化导致的点云高度模糊求解出现错误的问题；
3)本发明方法本质上将阵列干涉SAR点云高度模糊求解问题归结为干涉处理当中的相位解缠过程。传统干涉相位解缠技术发展相对成熟，多种现有解缠算法可以在本发明方法中直接进行替代。
附图说明
图1是本发明的一种基于相位解缠的阵列干涉SAR山区点云解模糊方法的算法处理流程图；
图2(a)，图2(b)，图2(c)，图2(d)是本发明方法处理真实Array-InSAR点云数据的部分结果；其中，图2(a)为高度模糊点云，图2(b)为点云高度矩阵(滤波前)，图2(c)为点云高度矩阵(滤波后)，图2(d)为解模糊后点云。
具体实施方式
下面结合附图并举实例，对本发明进行详细描述。
如图1所示，本发明提供的基于相位解缠的阵列干涉SAR山区点云解模糊方法包括以下步骤：
步骤1、提取点云高度矩阵。算法输入为Array-InSAR系统通过对雷达回波数据进行处理后生成的山区三维点云，即阵列干涉SAR山区模糊点云。为便于说明，这里生成的点云坐标系为方位-斜距-高程坐标系。然后提取点云高度矩阵，机载阵列干涉SAR生成的点云结果可以用一大小为M×N×K的多维数组表示，其中M、N、K分别表示方位、斜距和高程向的像素数目，点云高度矩阵通过对每一方位-斜距单元提取所有高程向上幅度值最大的点云所在的高程位置形成，此时形成的M×N大小的图像的每一个像素元素代表了点云的高程值；
步骤2、点云高度矩阵滤波。受噪声、配准误差、重建算法等因素影响，步骤1提取的点云高度矩阵存在噪声，影响点云高度结果，故需要对点云高度矩阵进行滤波处理，得到高质量的点云高度值。这里采用干涉相位滤波方法，将点云高度矩阵转化为干涉相位，此时矩阵元素(m,n)处的点云的高度值sm,n和相位P(m,n)之间的关系如下：
其中，sm,n为点云的高程值，L为高程向的点云高度最大值。经过上述转化后，图像的像素的取值范围变为(-π,π]。采用圆周均值滤波方法对上述图像进行滤波处理，令滤波窗口为以待滤波元素为中心，前后大小为(2M1+1)×(2N1+1)的矩形，M1、N1分别表示矩阵行和列方向上像素的数目。待滤波元素为圆周均值滤波的输出值即：
上式当中，为中间参数，arg{}表示取角操作，p和q分别表示对应矩阵的行和列，j为虚数单位，exp()表示指数函数。
步骤3、阴影区检测和掩膜处理。图像中阴影区会影响解缠的结果，故需要在进行解缠之前对图像的阴影区进行掩膜处理。图像阴影区在SAR幅度图像中表现为较暗的区域，相应的点云高度呈噪声特性。利用两通道图像之间计算相关系数图，利用相关系数图和SAR幅度图像进行图像检测并提取阴影区域，并利用形态学方法使区域边缘规整化。
步骤4、利用枝切法解高度模糊。即对进行掩膜处理后的点云高度图利用基于枝切法的路径跟踪算法对点云高度进行解模糊，得到无模糊点云高度值。这里得到的点云高度值结果是以某一位置为参考的相对高度值，故要得到绝对的点云高度值需要知道该位置的高度值，可通过外源DEM辅助得到。
步骤5、阴影区高度插值。步骤4得到无模糊点云高度值，为了保留原始点云阴影区的结果，对阴影区的点云高度进行插值处理，沿图像每一方位向，采用线性插值方法对存在空洞的图像位置进行插值，补偿图像阴影区的点云的高度值，此时得到的图像中的每一个元素代表了该位置处点云的高程值大小，最终获得点云高度图。
步骤6、点云提取。步骤5得到点云高度图后，需要在原始点云中提取出正确的点云，包括点云的幅度和高程位置。首先对原始点云在高程向进行周期伸展，然后设置高程向掩膜的宽度，其次根据步骤5得到的点云高度图，按照每一个方位向，沿斜距方向逐像素进行点云提取，得到每一个像素位置点云的高程大小和点云强度值。
为了具体表现本发明的优越性，本实施例中针对本发明所提方法进行实验说明。利用机载Array-InSAR获取得到西南地区某一山区SAR图像。该区域地势由左到右高度逐渐减小，且高度落差大于系统的最大模糊高度。利用本发明所提方法对机载Array-InSAR系统获取的该区域点云数据进行处理，部分处理结果如图2(a)，图2(b)，图2(c)，图2(d)所示。由于该区域的高度差超过了系统的最大不模糊高度，通过Array-InSAR点云重建方法获取得到的点云存在高度模糊现象，如图2(a)所示，利用本发明方法对Array-InSAR山区点云进行解模糊处理，首先提取点云高度矩阵，如图2(b)，从图2(b)可发现处理得到的点云高度矩阵存在较大噪声，需要在相位解缠之前对点云高度矩阵进行滤波处理，图2(c)为进行滤波处理后的点云高度结果，此时图像结果当中噪声得到抑制，图像质量得到较大提升。图2(d)显示了完成本发明方法所有处理步骤后获取的山区点云结果，从点云对比结果可得到本发明方法能够对模糊的山区点云进行解模糊处理，获得正确的山区点云。
综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (3)
1.一种基于相位解缠的阵列干涉SAR山区点云解模糊方法，其特征在于，包括以下步骤：
步骤1、提取阵列干涉SAR山区点云的点云高度矩阵；
输入阵列干涉SAR通过对雷达回波数据进行处理后生成的山区三维点云，获得需要处理的阵列干涉SAR山区模糊点云，点云坐标系为方位-斜距-高程坐标系；然后提取点云高度矩阵，机载的阵列干涉SAR生成的点云结果用一大小为的多维数组表示，其中分别表示方位向、斜距向和高程向的像素数目，点云高度矩阵通过对每一方位-斜距单元提取所有高程向上幅度值最大的点云所在的高程位置形成，此时形成/>大小的图像，该图像的每一个像素元素代表点云的高程值；
步骤2、进行点云高度矩阵的滤波；
采用干涉相位滤波方法，将点云高度矩阵转化为干涉相位，此时点云高度矩阵的元素处的点云的高程值/>和相位/>之间的关系如下：
；
其中，为点云的高程值，/>为高程向的点云高度最大值；经过上式转化后，/>大小的图像的像素元素的取值范围变为/>；采用圆周均值滤波方法对/>大小的图像进行滤波处理，令滤波窗口为以待滤波的像素元素为中心、前后大小为/>的矩形，/>、/>分别表示行和列方向上像素元素的数目，待滤波的像素元素为圆周均值滤波的输出值/>，即：
；
上式当中，，为中间参数，arg{}表示取角操作，p和q分别表示对应矩阵的行和列，j为虚数单位，exp（）表示指数函数；
步骤3、进行阴影区检测和掩膜处理，获得掩膜处理后的点云高度图；
步骤4、对掩膜处理后的点云高度图利用基于枝切法的路径跟踪算法对点云高度结果进行解模糊，得到无模糊点云高度值；
步骤5、对无模糊点云高度值进行阴影区插值，获得阴影区点云高度；
步骤6、进行点云提取，得到每一个像素位置的点云的高程大小和点云强度值。
2.根据权利要求1所述的一种基于相位解缠的阵列干涉SAR山区点云解模糊方法，其特征在于，所述步骤3中，利用阵列干涉SAR的首尾两通道图像，计算首尾两通道图像之间的相关系数图，利用相关系数图和阵列干涉SAR幅度图像进行图像检测并提取阴影区，并利用形态学方法使阴影区的区域边缘规整化。
3.根据权利要求2所述的一种基于相位解缠的阵列干涉SAR山区点云解模糊方法，其特征在于，所述步骤5包括对阴影区的点云高度进行插值处理，沿图像每一方位向，采用线性插值方法对存在空洞的图像位置进行插值，补偿图像的阴影区的点云的高度值，此时得到的图像中的每一个像素元素代表该像素元素的位置处的点云的高程值大小。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311621515.9A CN117872366B (zh) | 2023-11-30 | 2023-11-30 | 一种基于相位解缠的阵列干涉sar山区点云解模糊方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311621515.9A CN117872366B (zh) | 2023-11-30 | 2023-11-30 | 一种基于相位解缠的阵列干涉sar山区点云解模糊方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117872366A CN117872366A (zh) | 2024-04-12 |
| CN117872366B true CN117872366B (zh) | 2024-06-07 |
Family
ID=90578017
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311621515.9A Active CN117872366B (zh) | 2023-11-30 | 2023-11-30 | 一种基于相位解缠的阵列干涉sar山区点云解模糊方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117872366B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1998002761A1 (en) * | 1996-07-11 | 1998-01-22 | Science Applications International Corporation | Terrain elevation measurement by interferometric synthetic aperture radar (ifsar) |
| CN109541593A (zh) * | 2018-10-30 | 2019-03-29 | 北京航空航天大学 | 一种改进的最小费用流InSAR相位解缠方法 |
| CN110988876A (zh) * | 2019-11-21 | 2020-04-10 | 长沙理工大学 | 闭合鲁棒双基线InSAR相位解缠方法、系统及可读存储介质 |
| CN111948654A (zh) * | 2020-08-12 | 2020-11-17 | 中国科学院空天信息创新研究院 | 机载层析sar三维点云生成方法 |
| CN115616575A (zh) * | 2022-11-28 | 2023-01-17 | 中国科学院空天信息创新研究院 | 一种星载sar立体测量辅助的干涉相位图解缠绕方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108663678B (zh) * | 2018-01-29 | 2022-01-18 | 西北农林科技大学 | 基于混合整数优化模型的多基线InSAR相位解缠算法 |
| KR102142674B1 (ko) * | 2019-08-01 | 2020-08-07 | 서울시립대학교 산학협력단 | 위성레이더 오프셋트래킹 지표변위 모델 기반 위성레이더 간섭영상의 불구속화 방법 및 그 장치 |
| CN114114181B (zh) * | 2022-01-28 | 2022-04-22 | 中国科学院空天信息创新研究院 | 基于轨道误差相位基的星载sar干涉基线矫正方法 |
| CN115421145A (zh) * | 2022-09-22 | 2022-12-02 | 中国科学院空天信息创新研究院 | 基于外源dem的阵列sar山区点云高程解模糊方法 |
-
2023
- 2023-11-30 CN CN202311621515.9A patent/CN117872366B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1998002761A1 (en) * | 1996-07-11 | 1998-01-22 | Science Applications International Corporation | Terrain elevation measurement by interferometric synthetic aperture radar (ifsar) |
| CN109541593A (zh) * | 2018-10-30 | 2019-03-29 | 北京航空航天大学 | 一种改进的最小费用流InSAR相位解缠方法 |
| CN110988876A (zh) * | 2019-11-21 | 2020-04-10 | 长沙理工大学 | 闭合鲁棒双基线InSAR相位解缠方法、系统及可读存储介质 |
| CN111948654A (zh) * | 2020-08-12 | 2020-11-17 | 中国科学院空天信息创新研究院 | 机载层析sar三维点云生成方法 |
| CN115616575A (zh) * | 2022-11-28 | 2023-01-17 | 中国科学院空天信息创新研究院 | 一种星载sar立体测量辅助的干涉相位图解缠绕方法 |
Non-Patent Citations (2)
| Title |
|---|
| 一种改进的掩膜割线相位解缠算法;钟何平;唐劲松;张森;张学波;;大地测量与地球动力学;20131015(第05期);第101页-第105页 * |
| 星载高分辨聚束模式SAR数据干涉处理;巩春和;;甘肃科技;20141030(第20期);第23页-第26页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117872366A (zh) | 2024-04-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108596103B (zh) | 基于最佳光谱指数选择的高分辨遥感影像建筑物提取方法 | |
| CN112686935B (zh) | 基于特征融合的机载测深雷达与多光谱卫星影像配准方法 | |
| CN103605119B (zh) | 一种星载合成孔径雷达条带模式方位模糊抑制方法 | |
| CN109859219B (zh) | 结合相位与光谱的高分遥感影像分割方法 | |
| CN111311644B (zh) | 一种基于视频sar的运动目标检测方法 | |
| CN110070545B (zh) | 一种城镇纹理特征密度自动提取城镇建成区的方法 | |
| CN116863357A (zh) | 一种无人机遥感堤坝影像定标与智能分割变化检测方法 | |
| CN115951350A (zh) | 永久散射体点提取方法、装置、设备及介质 | |
| CN117687014B (zh) | 一种基于二维滤波MapDrift的SAR等效雷达速度估计方法 | |
| CN115540908A (zh) | 基于小波变换的InSAR干涉条纹匹配方法 | |
| CN110095774B (zh) | 一种圆迹视频sar动目标检测方法 | |
| CN108509835B (zh) | 基于DFIC超像素的PolSAR图像地物分类方法 | |
| CN117872366B (zh) | 一种基于相位解缠的阵列干涉sar山区点云解模糊方法 | |
| CN106910178B (zh) | 一种基于色调统计特性分类的多角度sar图像融合方法 | |
| CN109615584B (zh) | 一种基于单应性约束的sar图像序列map超分辨率重建方法 | |
| CN114463417B (zh) | 一种基于空间几何比值关系的sar影像叠掩区域定位方法 | |
| CN112819909B (zh) | 基于低分辨率先验光谱图像区域分割的自适应编码方法 | |
| Zhu et al. | Global LoD-1 building model from TanDEM-X data | |
| CN115601278A (zh) | 基于子图像配准的高精度运动误差补偿的方法 | |
| CN115407282A (zh) | 一种短基线下基于干涉相位的sar有源欺骗干扰检测方法 | |
| CN111144246B (zh) | 利用多时相sar影像和光学辅助信息的道路提取方法 | |
| CN114066815A (zh) | 基于模糊c均值聚类的SAR图像变化检测方法 | |
| CN112883796A (zh) | 基于叠掩信息的sar图像多类型建筑物高度估计方法 | |
| CN116008927B (zh) | 一种基于psp-ca的低波段星载sar电离层闪烁相位误差校正方法 | |
| CN117310705B (zh) | 一种基于双极化sar影像的洪涝灾害快速检测方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant |