基于平面投影的多基SAR三维表面重建方法技术

本发明专利技术公开了一种基于平面投影的多基SAR三维表面重建方法，利用散射点的距离历史，结合接收机和发射机之间的定量关系，实现三维表面的重建，首先根据一发两收的多基SAR构型建立回波模型并成像，提取相关散射点的位置信息及映射关系，随后反演恢复散射点的真实三维位置信息，实现三维表面重建。本发明专利技术的方法原理较为简单，并且可以充分利用构型信息，设置特定的接收机和发射机之间的定量关系对方程进行简化，从而提高重建精度和简化运算，可以准确地求得平面投影目标的真实三维位置，实现高度恢复以及多基SAR的三维表面重建。度恢复以及多基SAR的三维表面重建。度恢复以及多基SAR的三维表面重建。

【技术实现步骤摘要】
基于平面投影的多基SAR三维表面重建方法


[0001]本专利技术属于合成孔径雷达(SAR)成像
，具体涉及一种基于平面投影的多基SAR三维表面重建方法。

技术介绍

[0002]合成孔径雷达(Synthetic Aperture Radar,SAR)是一种具有全天时、全天候、远距离等特点的高分辨微波主动成像雷达，然而传统SAR二维成像技术缺少高度维的信息，无法反映目标真实的位置和结构，也不利于后续的信息获取与处理，因此高程信息的恢复是SAR的核心问题之一。SAR的三维重建技术能够更加真实地恢复目标场景，对城市测绘、目标精确定位和灾害评估等应用具有重大价值和意义。
[0003]传统的合成孔径雷达成像是基于距离多普勒的原理，目标会聚焦于零多普勒
‑
最近斜距处，该成像原理使得SAR仅仅具备二维成像的能力，对于三维空间，SAR图像是场景向平面的投影，无法获得第三维的信息。为了解决此问题，现有提出的方法，例如立体测量与雷达成像技术结合，但该方法主要是基于灰度信息，没有利用图像的相位信息，而干涉测量，即通过对干涉相位即两个天线得到的图像之间存在的相位差进行估计，但在数据复杂的情况下相位解缠等过程较难处理，此外一些方法与压缩感知相结合，但对成像构型和场景的要求较高，计算处理也较为复杂。

技术实现思路

[0004]为解决上述技术问题，本专利技术提出了一种基于平面投影的多基SAR三维表面重建方法。
[0005]本专利技术采用的技术方案为：一种基于平面投影的多基SAR三维表面重建方法，具体步骤如下：
[0006]S1、根据一发两收的多基SAR构型建立回波模型并成像；
[0007]S2、提取相关散射点的位置信息及映射关系；
[0008]S3、反演恢复散射点的真实三维位置信息；
[0009]S4、实现三维表面重建。
[0010]进一步地，所述步骤S1中，具体如下：
[0011]S11、采用一个发射站两个接收站组成的多基SAR系统，即形成一发两收的空间几何模型；
[0012]场景中心点位于(0,0,0)，波束中心时刻发射站的位置为r
t0
＝(x
t
,y
t
,z
t
)，两个接收站的位置分别为r
r0(i)
＝(x
ri
,y
ri
,z
ri
),i＝1,2，假设发射站的速度为v
t
＝(v
xt
,v
yt
,0)，接收站的速度分别为v
ri
＝(v
xi
,v
yi
,0),i＝1,2，即发射站与接收站沿着轨迹匀速直线飞行，则在方位时刻η，发射站与第i个接收站的位置矢量分别表示为：
[0013]r
t
(η)＝r
t0
+v
t
η
[0014]r
r(i)
(η)＝r
r0(i)
+v
ri
η(i＝1,2)
[0015]S12、发射站发出chrip信号，每个发射脉冲在被场景散射点r
m
＝(x
m
,y
m
,z
m
)反射经各个接收站分别接收解调后，第i个接收站在方位时刻η处接收到的回波建模为：
[0016][0017]其中，τ表示距离快时间变量，T
p
表示信号的脉冲宽度，f
c
表示信号的载波频率，K
r
表示发射线性调频信号的斜率，A0表示目标的散射系数，η表示方位时刻，c表示光速，R
(i)
(η)表示距离历史，即在η时刻目标到发射站位置与目标到第i个接收站位置的距离之和。
[0018]S13、使用后向投影(Back Projection,BP)算法对回波处理并成像，最终可得到两幅SAR图像。
[0019]进一步地，所述步骤S2中，具体如下：
[0020]S21、对于空间目标r
p
＝(x,y,z)，由步骤S1中的构型可知其距离历史R
(i)
(η)为发射站位置到目标的距离与目标到第i个接收站位置的距离之和，则将其表示为：
[0021][0022]其中，R
t
(η)表示发射站位置到实际空间目标r
p
的距离，表示该目标到第i个接收站位置的距离。
[0023]S22、对于空间目标r
p
＝(x,y,z)所对应的成像处理后两幅SAR图像中像素点r
fi
＝(x
fi
,y
fi
,0)，将其对应的距离历史R
f(i)
(η)表示为：
[0024][0025]其中，RT
i
'(η)表示发射站位置到实际目标成像后SAR图像中目标像素点r
fi
的距离，RR
i
'(η)表示该目标像素点到第i个接收站位置的距离。
[0026]S23、接收目标r
P
＝(x,y,z)的回波并成像处理后，仅考虑地面上的目标(高度为零)，将r
p
投影至(x
j
,y
k
,0)，则有：
[0027]R(η；x,y,z)＝R(η；x
j
,y
k
,0)
[0028]其中，R(η；x,y,z)表示目标r
p
的距离历史，R(η；x
j
,y
k
,0)表示成像处理后像素点的距离历史。
[0029]S24、目标高度非零的情况，依然满足上式，对于步骤S1中的构型，则空间目标和SAR图像像素距离历史之间的关系表示为：
[0030]R
(i)
(η)＝R
f(i)
(η)
[0031]进一步地，所述步骤S3中，具体如下：
[0032]S31、由步骤S2中距离历史之间的对应关系，对其分别在方位时刻为0(η＝0)处进行泰勒展开：
[0033][0034][0035]S32、对于两幅图像中同一目标散射点，可得下列方程组：
[0036][0037]S33、由步骤S32得到一个非线性方程组，使用牛顿迭代法或最小二乘法求解，最终能够得到未知数(x,y,z)，即求得空间目标的三维位置矢量r
p
。
[0038]进一步地，所述步骤S4中，具体如下：
[0039]通过两幅双基地SAR图像中对应散射点的位置信息，求解步骤S33中得到的非线性方程组来反求目标的真实位置，将恢复得到的散射点组合形成点云，完成初步三维表面重建。
[0040]本专利技术的有益效果：本专利技术的方法利用散射点的距离历史，结合接收机和发射机之间的定量关系，实本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于平面投影的多基SAR三维表面重建方法，具体步骤如下：S1、根据一发两收的多基SAR构型建立回波模型并成像；S2、提取相关散射点的位置信息及映射关系；S3、反演恢复散射点的真实三维位置信息；S4、实现三维表面重建。2.根据权利要求1所述的一种基于平面投影的多基SAR三维表面重建方法，其特征在于，所述步骤S1中，具体如下：S11、采用一个发射站两个接收站组成的多基SAR系统，即形成一发两收的空间几何模型；场景中心点位于(0,0,0)，波束中心时刻发射站的位置为r
t0
＝(x
t
,y
t
,z
t
)，两个接收站的位置分别为r
r0(i)
＝(x
ri
,y
ri
,z
ri
),i＝1,2，假设发射站的速度为v
t
＝(v
xt
,v
yt
,0)，接收站的速度分别为v
ri
＝(v
xi
,v
yi
,0),i＝1,2，即发射站与接收站沿着轨迹匀速直线飞行，则在方位时刻η，发射站与第i个接收站的位置矢量分别表示为：r
t
(η)＝r
t0
+v
t
ηr
r(i)
(η)＝r
r0(i)
+v
ri
η(i＝1,2)S12、发射站发出chrip信号，每个发射脉冲在被场景散射点r
m
＝(x
m
,y
m
,z
m
)反射经各个接收站分别接收解调后，第i个接收站在方位时刻η处接收到的回波建模为：其中，τ表示距离快时间变量，T
p
表示信号的脉冲宽度，f
c
表示信号的载波频率，K
r
表示发射线性调频信号的斜率，A0表示目标的散射系数，η表示方位时刻，c表示光速，R
(i)
(η)表示距离历史，即在η时刻目标到发射站位置与目标到第i个接收站位置的距离之和；S13、使用后向投影算法对回波处理并成像，最终可得到两幅SAR图像。3.根据权利要求1所述的一种基于平面投影的多基SAR三维表面重建方法，其特征在于，所述步骤S2中，具体如下：S21、对于空间目标r
p
＝(x...

【专利技术属性】
技术研发人员：王雯璟，刘玲，海宇，武俊杰，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：