一种基于极化SAR多参数的月球水冰探测方法

一种基于极化sar多参数的月球水冰探测方法
技术领域
1.本发明属于雷达技术领域，特别涉及一种基于极化sar（合成孔径雷达）的月球水冰探测方法。


背景技术：

2.极化sar是一种发射或接收多个极化状态电磁波的成像雷达。水冰作为一种低介电损耗介质，在一定条件下会产生相干后向散射反相效应(coherent backscatter opposition effect，cboe)。当雷达波在弱吸收介质（例如，水冰、风化层）中传播时，由于异质性（例如，埋藏的岩石、空隙）会发生一系列随机散射，那些沿相同路径但以相反方向散射的波相干相加，导致反向散射雷达信号具有很高的圆极化率比（circular polarization ratio, cpr）。因此极化sar被认为是探测水冰最有效的工具之一。月球两极部分区域因长时间范围内无太阳照射而形成永久阴影区（permanently shadowed region，psr）。psr内部温度极低，可以存储挥发性物质，因此月球两极可能有水冰的存在。目前基于极化sar对月球水冰探测方法主要是利用cpr。cpr是指接收的同向圆极化回波与反向圆极化回波功率的比值，当月球两极存在有水冰时，其cpr参数会因cboe而大于1。基于cpr参数，我们可以推断月球两极psr内部是否存在水冰以及存在的位置。相较于其他载荷，极化sar因其较高的分辨率可以对水冰进行较精准的定位。
3.然而，基于cpr参数的月球水冰探测方法受到了很多质疑。原因在于cpr大于1并不能准确地指向水冰。事实表明，粗糙表面相应的cpr值同样可大于1。月球陨石坑沉积物、熔岩流和岩石表面等地质目标的cpr可高达2到4。因此，仅仅依据cpr指标，我们无法将水冰和粗糙的月壤表面区分开来。当前，月球两极存在的异常陨石坑成为了研究热点。异常陨石坑是指陨石坑内部区域cpr大于1，而外部区域cpr远小于1。陨石坑内外cpr值的差异被解释为水冰的存在导致。陨石坑内部区域因存在psr而储存水冰，外部区域因太阳直射不可能有水冰存在。但也有学者将异常陨石坑解释为陨石坑内外区域的粗糙度不同，而非水冰的因素。为了准确地探测水冰，需要从cpr大于1的区域中剔除粗糙表面区域及有沉积物、石块等地质目标干扰的区域。


技术实现要素：

4.鉴于此，本发明的主要目的在于提供一种基于极化sar多参数的月球水冰探测方法，通过对异常陨石坑的极化sar数据进行分析，基于多种参数选择水冰存在的区域，剔除粗糙的月壤表面及沉积物、石块等干扰因素，提高了水冰探测精度。
5.为达到上述目的，本发明采用的技术方案如下：一种基于极化sar多参数的月球水冰探测方法，包括如下步骤：（1）选择月球两极异常陨石坑作为研究区域，并计算cpr参数，选择cpr参数大于1的区域；所述cpr参数为圆极化率比；（2）基于研究区域的极化sar数据和m-δ分解方法，计算各个散射分量的功率，选择
体散射占主导地位的区域；（3）计算研究区域的粗糙度参数r，并根据所述粗糙度参数r选择相对平坦的区域；（4）计算研究区域的后向散射系数σ，并根据所述后向散射系数σ选择表面石块较少的区域；（5）选择研究区域中同时满足步骤（1）至步骤（4）条件的区域作为最终的水冰存在位置。
6.进一步地，所述步骤（1）包括：基于stokes参数s1，s2，s3和s4，计算cpr参数，计算公式表示为：cpr=(s
1-s4)/(s1 s4)其中，s1表示接收回波的总功率，s2和s3表示接收回波中不同方向上线极化分量的功率，s4表示接收回波中左旋圆极化和右旋圆极化分量的功率和。
7.进一步地，所述步骤（2）包括：所述m-δ分解方法中包括极化度参数m，其表征雷达回波中的极化度特性，基于stokes参数，所述极化度参数m表示为：所述m-δ分解方法中包括相对相位参数δ，其为水平和垂直极化通道的相位差；基于stokes参数，相对相位参数δ表示为：δ=tan-1
(-s4/s3)利用m-δ分解方法，得到表面散射f
surface
、二次散射f
double
及体散射f
volume
的功率表达式，选择体散射占主导地位的区域，即满足：f
volume
》(f
double
f
surface
)；通过解上述表达式得到体散射占主导地位的区域其极化度参数m要满足：0≤m≤0.33。
8.进一步地，所述步骤（3）包括：根据极化sar数据的散射矩阵和平坦区域的散射矩阵进行相似度度量，得到粗糙度参数r，所述粗糙度参数r越大，表明该区域越不平坦，所述粗糙度参数r越小，表明区域越平坦；利用所述粗糙度参数r选择研究区域中平坦的部分。
9.进一步地，所述步骤（4）包括：当所述后向散射系数σ大时，目标场景的表面伴随的石块多，当所述后向散射系数σ小时，目标场景的表面伴随的石块少，利用所述后向散射系数σ选择研究区域中石块少的部分。
10.有益效果：和现有水冰探测方法相比，本发明首先将研究区域聚焦在月球两极的异常陨石坑。和其他区域相比，异常陨石坑存在水冰的可能性更大，这提高了水冰探测精度。其次，和单一地仅使用cpr参数来探测水冰的方案相比，本发明还使用了m-δ分解方法得到各个散射分量的功率，通过选择体散射占主导的区域进一步提升水冰存在的可能性。本发明首次提出了基于极化sar数据的粗糙度参数r并结合后向散射系数σ来排除粗糙表面及石块对探测结果的干扰。
附图说明
11.图1为采用本发明的基于极化sar多参数的月球水冰探测方法进行月球水冰探测的流程图；图2为月球北极异常陨石坑rozhdestvenskiy n的cpr参数图像，cpr=0为深色，cpr=1为浅色；图3为月球北极异常陨石坑rozhdestvenskiy n的极化度参数m图像，m=0为深色，m=1为浅色；图4为月球北极异常陨石坑rozhdestvenskiy n的粗糙度参数r图像，r=0为深色，r=3为浅色；图5（a）为月球北极异常陨石坑rozhdestvenskiy n的后向散射系数σ图像，陨石坑为图中圆形；图5（b）为基于本发明的基于极化sar多参数的月球水冰探测方法得到的水冰存在的区域，通过白色像素点表示。
具体实施方式
12.为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
13.本发明提出一种基于极化sar多参数的月球水冰探测方法。图1为本发明的基于极化sar多参数的月球水冰探测方法的流程图。
14.根据该流程图，所述探测方法包括如下步骤：第一步，选择极区异常陨石坑的极化sar数据，计算cpr参数。
15.所述极化sar数据来自于月球勘测轨道飞行器(lunar reconnaissance orbiter，lro)上的mini-rf载荷。mini-rf是一个混合极化sar系统，它记录了雷达回波在水平电场矢量e
lh
垂直电场矢量e
lv
的数据。这些数据可用于生成stokes矢量s，如下：其中，re()表示取实部操作，im()表示取虚部操作，*表示复共轭。stokes参数s1表示接收回波的总功率，s2和s3表示接收回波中不同方向上线极化分量的功率，s4表示接收回波中左旋圆极化和右旋圆极化分量的功率和。
16.cpr参数可以由stokes参数计算得到，如下：
因水冰在一定条件下可产生相干后向散射反相效应cboe（coherent backscatter opposition effect）从而使cpr参数大于1。因此，基于cpr参数，可以初步获得异常陨石坑中符合条件的区域。但由于cpr参数大于1不能唯一地和水冰相对应，月球熔融物，陨石坑溅射区域等粗糙表面都可使cpr参数大于1，因此，还需要对cpr参数大于1的区域进行更详细的研究。
17.第二步，基于研究区域的极化sar数据和m-δ分解方法计算研究区域各个散射分量的功率，并选择体散射占主导地位的区域。
18.m-δ分解方法中的极化度参数m是指雷达回波中极化功率与总功率的比值，并表征了后向散射的极化特性，可由stokes参数计算得到，如下：由定义可知极化度参数m取值范围为：0≤m≤1。m越小，表明目标场景越复杂，m越大，表明目标场景越单一。m-δ分解方法中的相对相位参数δ为水平和垂直极化通道的相位差，基于stokes参数，相对相位参数δ可表示：δ=tan-1
(-s4/s3)早期的研究表明：水冰存在的区域往往伴随着体散射，满足以下条件：f
volume
》(f
double
f
surface
)f
volume
,f
double
,f
surface
分别表示体散射、二次散射和表面散射分量功率。利用m-δ分解方法可知三种散射分量功率可表示为：当体散射占主导地位时，满足以下条件：求解不等式可得：0≤m≤0.33。因此，根据极化度参数m选择体散射占主导地位的区域。
19.第三步，计算粗糙度参数r，并基于该粗糙度参数r选择相对平坦的区域。
20.首先介绍极化sar数据的散射矩阵∑与stokes矢量的关系。所述散射矩阵∑和stokes矢量均为混合的极化sar数据的存储方式，二者可以等价转换，关系如下：
其中，j为虚数；由于水冰和粗糙的表面都可以使cpr参数大于1，因此，需要评估极区异常陨石坑内部粗糙程度，并选择相对平坦的区域。粗糙度参数r定义为：∑表示极化sar数据的散射矩阵，∑f表示平坦区域所对应的散射矩阵，这里取∑m=(∑ ∑f)/2表示∑和∑f的均值，表示矩阵∑m的逆。tr(
·
)表示对矩阵求迹，var(
·
)表示方差运算。
21.粗糙度参数r表示基于极化sar的散射矩阵对两个区域之间的相似度进行度量，r越大，相似度越低，反之，r越小，相似度越高。这里，选择将极区异常陨石坑的极化sar数据和平坦区域的极化sar数据进行比较。基于结果，可以对极区异常陨石坑内部的粗糙程度进行定量分析。这里采用0≤r≤1的取值范围选择极区异常陨石坑中相对平坦的区域。
22.第四步，计算后向散射系数σ，并基于后向散射系数σ选择石块较少的区域。
23.极化sar对于波长尺度大小的石块十分敏感，表面存在大量石块，往往有着很大的后向散射系数σ，石块较小的表面其回波能量很弱，后向散射系数σ往往较低。这里，对极区异常陨石坑内的后向散射系数σ进行分析，并设置石块较小区域的后向散射系数σ的取值范围为:σ≤0.03。根据该取值范围选择石块少的区域。最后，选择极区异常陨石坑中同时满足第一步至第四步条件的区域作为最终的水冰存在位置。
24.下面结合具体的极化sar的数据对本发明的技术方案作进一步详细说明。实例选择的是月球北极的一个异常陨石坑rozhdestvenskiy n（84.0
°
n, 156.5
°
w; 8.6 km 直径）。图2给出了rozhdestvenskiy n的cpr参数图像。从图2可以看出：陨石坑内部cpr参数值远远高于外部。图3给出了rozhdestvenskiy n的极化度参数m图像。从图3可知，异常陨石坑外部的m值比内部要大。图4给出了rozhdestvenskiy n的粗糙度参数r的图像。从图中可以看出，陨石坑内部的粗糙程度差别很大，左侧粗糙程度较大，右侧较小。图5（a）给出异常陨石坑rozhdestvenskiy n的后向散射系数图，图5（b）给出异常陨石坑内部白色框中最终的水冰的存在位置，用白色像素点表示。
25.以上所述，仅为本发明的某些实施例而已，并非用于限定本发明的保护范围。本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。