CN105891535A - 一种用于返回式航天器的分离速度测量方法 - Google Patents

Abstract

本发明提出的一种用于返回式航天器的分离速度测量方法，利用单目监视相机连续成像拍摄返回式航天器与服务舱分离的图像序列，将图像序列按照从分离前至完成分离的时间先后顺序进行人工排序并进行畸变校正，通过对图像进行色彩分割、物体边缘检测、获取边缘，得到分离速度。该测量方法利用重量较轻的监视相机，实现了一种较为简便的返回式航天器在轨分离速度测量。

Description

一种用于返回式航天器的分离速度测量方法

技术领域

本发明涉及返回式航天器在分离时刻的分离速度测量方法，具体为利用经过标定的单目监视相机，实现在分离时刻的相对速度测量。

背景技术

传统航天器分离速度测量主要依赖于雷达测量或者激光测距仪。在地面使用雷达进行相对速度测量时，在航天器中安装激光测距仪进行连续距离测量时能够获得较高的精度，但激光测距仪设备复杂、占用资源多、费效比低。

发明内容

本发明的技术解决问题是利用单目监视相机，在深空背景环境中对返回式航天器(以下简称为“返回器”)分离过程进行拍摄，得到精确的相对分离速度。

本发明提出的一种用于返回式航天器(返回器)的分离速度测量方法，包括以下步骤：

1)利用单目监视相机连续成像拍摄返回式航天器与服务舱分离的图像序列；

2)图像序列按照从分离前至完成分离的时间先后顺序进行排序并通过单目监视相机畸变标定的参数校正图像畸变；

3)对图像进行色彩分割，将返回式航天器与深空背景分离；

4)通过canny算子对步骤3)获得的图像进行边缘检测，用于获得返回式航天器边缘，以便计算返回式航天器长轴在像面的长度；

5)根据所提取的返回式航天器的边缘进行霍夫变换获得返回式航天器底面椭圆形结构的长轴长度；

6)根据步骤5)获得的椭圆的长轴长度，通过小孔模型计算出每张图像对应的返回式航天器与相机的距离；

7)根据步骤2)所得到图像序列以及步骤6)获得的每张图像对应的返回式航天器与相机的距离，通过最小二乘法拟合，得到分离速度。

本发明提出的测量方法利用重量较轻的监视相机，实现了一种较为简便的返回式航天器在轨分离速度测量。

附图说明

图1是返回器分离过程图像；图1(a)-(v)为分离过程图像手动排序的序列图；

图2是返回器与深空背景颜色分割图；

图3是返回器底面椭圆形边缘提取示意图；

图4是通过霍夫变换获得返回器底面椭圆形结构图；

图5是返回器分离速度拟合结果。

具体实施方式

下面结合附图及具体实施方式对本发明进行详细说明。

对于安装在服务舱上的拍摄返回式航天器(返回器)图像所用的单目监视相机的要求如下：

(1)单目监视相机分辨率至少需要达到640*480像素。

(2)单目监视相机安装要求在分离过程中，至少有2帧图像能够完整覆盖返回器。相机的光轴指向需要指向分离场景的中心处。为了实现完整覆盖，相机视场角需要满足以下算式：

2arctan(D/2/(L₀+v_snt))＜θ＜2arctan(D/2/(L₀+v_s(n+1)t)) (0.1)

其中θ为所需要的相机视场角，D为分离场景的最大包络直径，L₀为未分离时刻分离场景中心与相机的距离，v_s为分离速度，t为连续拍摄间隔时间，n为第n帧拍摄到的分离图像。满足式(0.1)即满足了在分离过程中2帧图像能够完整覆盖返回器。

(3)单目监视相机需要经过畸变标定，通过畸变标定确定参数，以用来校正返回器图像的畸变。畸变的数学模型如式所示：

δ_x＝k₁x(x²+y²)+(p₁(3x²+y²)+2p₂xy)+s₁(x²+y²)

δ_y＝k₂x(x²+y²)+(p₂(3x²+y²)+2p₁xy)+s₂(x²+y²) (0.2)

其中δ_x为x方向的畸变，δ_y为y方向的畸变，k₁，k₂，p₁，p₂，s₁，s₂为测得的畸变参数。通过畸变校正的图像符合相机成像小孔模型，该模型在本算法中用于求解返回器与相机的距离。

返回器实体应该满足以下要求：

1.返回器底部为正圆形。

2.返回器底部在正常曝光拍摄时，获取的图像与深空背景有明确的区分，即获得的返回器底部图像与深空背景的DN值差异大于40。

拍摄时需保证曝光时间使返回器的最大外轮廓可以在图像中提取。

对返回式航天器分离速度进行测量的主要步骤如下：

第一步，利用单目监视相机连续成像拍摄返回器与服务舱分离的图像序列，如图1所示，该序列经过手动排序。

第二步，图像序列按照从分离前至完成分离的时间先后顺序进行排序并进行畸变校正，以使畸变校正的图像符合相机成像小孔模型。

第三步，对图像进行色彩分割，将返回器与深空背景分离，色彩分割采用的算法如式(0.3)所示：

$I(x,y)=\left\{\begin{array}{cc}i(x,y)& (i(x,y)>T)\\ 255& (i(x,y)\≤T)\end{array}\right.---\left(0.3\right)$

T为分割阈值，根据在轨拍摄情况选择合适的分割值，该值计算方法如所示：

T＝(A_space+A_probe)/2 (0.4)

其中A_space为深空背景的平均DN值，A_probe为返回器的平均DN值。

以某返回航天器为例，对图像进行色彩分割，阈值取128，分割结果如图2所示

其中I(x，y)是(x，y)点处色彩分割后结果图像，i(x，y)是(x，y)点处的原始图像。如果是彩色图像，则i(x，y)的计算方法如式(0.5)所示：

i(x，y)＝(R(x，y)+G(x，y)+B(x，y))/3 (0.5)

R表示红色通道DN值，G表示绿色通道DN值，B表示蓝色通道DN值。

第四步，物体边缘检测：通过canny算子(Canny算子边缘检测是一组标准方法)对上一步的结果图像进行边缘检测，用于获得返回器边缘，边缘提取如图3所示，以便计算返回器直径(长轴)在像面的长度。

第五步，根据所提取边缘返回器的边缘进行霍夫变换(霍夫变换是一组标准方法)获椭圆信息，即获取返回器底面椭圆形结构的长轴长度，如图4所示。

第六步，根据步骤五获得的椭圆的长轴长度，通过小孔模型计算出每张图像对应的返回器与相机的距离，如式所示：

r＝D_phyf/Nd (0.6)

分离监视相机焦距为f，像元尺寸为d，返回器底面圆盘直径为N个像素，d为像元尺寸，D_phy为返回器底面圆盘的物理尺寸，r即为返回器与相机的距离。

第七步，根据步骤二所得到图像序列及步骤六获得的每张图像对应的返回器与相机的距离，通过最小二乘法拟合，得到分离速度，如图5所示。

设连续成像时图像拍摄时刻为t＝(t₁，t₂…t_n)，对应的每张图像所获得返回器与相机的距离为r＝(r₁，r₂…r_n)，则分离速度v如式(0.7)所示：

$v=\frac{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}(t_i-\stackrel{\‾}{t})(r_i-\stackrel{\‾}{r})}{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{(t_i-\stackrel{\‾}{t})}^2}---\left(0.7\right)$

本发明中未说明部分属于本领域的公知技术。

Claims (5)

1.一种用于返回式航天器的分离速度测量方法，其特征在于，包括以下步骤：

1)利用单目监视相机连续成像拍摄返回式航天器与服务舱分离的图像序列；

2)图像序列按照从分离前至完成分离的时间先后顺序进行排序并通过单目监视相机畸变标定的参数校正图像畸变；

3)对图像进行色彩分割，将返回式航天器与深空背景分离；

4)通过canny算子对步骤3)获得的图像进行边缘检测，用于获得返回式航天器边缘，以便计算返回式航天器长轴在像面的长度；

5)根据所提取的返回式航天器的边缘进行霍夫变换获得返回式航天器底面椭圆形结构的长轴长度；

6)根据步骤5)获得的椭圆的长轴长度，通过小孔模型计算出每张图像对应的返回式航天器与相机的距离；

7)根据步骤2)所得到图像序列以及步骤6)获得的每张图像对应的返回式航天器与相机的距离，通过最小二乘法拟合，得到分离速度。

2.根据权利要求1所述的一种用于返回式航天器的分离速度测量方法，其特征在于，步骤3)所述的色彩分割的方法如下：

$I(x,y)=\left\{\begin{array}{cc}i(x,y)& (i(x,y)>T)\\ 255& (i(x,y)\≤T)\end{array}\right.---\left(1\right),$

T为分割阈值，根据在轨拍摄情况选择合适的分割值，该值计算方法：

T＝(A_space+A_probe)/2 (2)，

其中A_space为深空背景的平均DN值，A_probe为返回式航天器的平均DN值，

其中I(x，y)是(x，y)点处色彩分割后结果图像，i(x，y)是(x，y)点处的原始图像；如果是彩色图像，则i(x，y)的计算方法：

i(x，y)＝(R(x，y)+G(x，y)+B(x，y))/3 (3)，

其中R表示红色通道DN值，G表示绿色通道DN值，B表示蓝色通道DN值；

步骤6)中当返回式航天器底面圆盘直径为N个像素时，所述的返回式航天器距离相机光心的距离：

r＝D_phyf/Nd (4)

其中r为返回式航天器底面距离相机光心的距离，d为像元尺寸，D_phy为返回式航天器底面圆盘的物理尺寸，f分离监视相机焦距为；

步骤7)中，当连续成像时图像拍摄时刻为t＝(t₁，t₂…t_n)，对应的每张图像所获得返回式航天器与相机的距离为r＝(r₁，r₂…r_n)，则分离速度v：

$v=\frac{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}(t_i-\stackrel{\‾}{t})(r_i-\stackrel{\‾}{r})}{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{(t_i-\stackrel{\‾}{t})}^2}---\left(5\right).$

3.根据权利要求1或2任一权利要求所述的一种用于返回式航天器的分离速度测量方法，其特征在于，所述的单目监视相机满足以下要求：

单目监视相机分辨率至少需要达到640*480像素；

单目监视相机安装要求在分离过程中，至少有2帧图像能够完整覆盖返回式航天器；

单目监视相机经过畸变标定，确定用以校正图像畸变的相关参数。

4.根据权利要求3所述的一种用于返回式航天器的分离速度测量方法，其特征在于，相机视场角满足以下要求，以实现至少有2帧图像能够完整覆盖返回式航天器：

2arctan(D/2/(L₀+v_snt))＜θ＜2arctan(D/2/(L₀+v_s(n+1)t)) (6)，

其中θ为所需要的相机视场角，D为分离场景的最大包络直径，L₀为未分离时刻分离场景中心与相机的距离，v_s为分离速度，t为连续拍摄间隔时间，n为第n帧拍摄到的分离图像；

所述的单目相机畸变方法为：

δ_x＝k₁x(x²+y²)+(p₁(3x²+y²)+2p₂xy)+s₁(x²+y²)

δ_y＝k₂x(x²+y²)+(p₂(3x²+y²)+2p₁xy)+s₂(x²+y²) (7)，

其中δ_x为x方向的畸变，δ_y为y方向的畸变，k₁，k₂，p₁，p₂，s₁，s₂为测得的畸变参数。

5.根据权利要求1或2任一权利要求所述的一种用于返回式航天器的分离速度测量方法，其特征在于，获得的返回式航天器底部图像与深空背景的DN值差异大于40。