发明名称 一种基于可见光图像的深空探测器自主着陆区选取方法
摘要 本发明公开一种基于可见光图像的深空探测器自主着陆区选取方法。重点解决在深空探测器着陆过程中利用可见光敏感器成像获得着陆区域地形数据,并判断着陆区域是否适合着陆的问题。首先根据成像分辨率、着陆器底座尺寸、相机FOV信息,通过优化自定义的全局可分和局部均匀度量,自适应地确定特征检测窗口尺寸;进而在自适应特征检测窗口尺寸下提取图像的统计、纹理和梯度特征,最后对图像纹理特性进行特征级融合,对图像统计特征、梯度特性及纹理特征信息进行决策级融合,实现对目标天体地形起伏的自适应判别,确定深空探测器自主着陆过程中的最佳自主着陆区。
申请公布号 CN106326826A 申请公布日期 2017.01.11
申请号 CN201510652643.9 申请日期 2015.10.10
申请人 北京控制与电子技术研究所 发明人 崔现军;刘辉;朱子奇;康鹏;宋小春;王云财;胡晓赛;侯建伟
分类号 G06K9/00(2006.01)I;G06K9/46(2006.01)I 主分类号 G06K9/00(2006.01)I
代理机构 核工业专利中心 11007 代理人 高尚梅
主权项 一种基于可见光图像的深空探测器自主着陆区选取方法,其特征在于具体步骤为:第一步构建基于可见光图像的深空探测器自主着陆区选取系统基于可见光图像的深空探测器自主着陆区选取系统,包括:自适应窗口大小确定模块、各纹理特征计算并融合模块和决策级融合并选取最佳着陆点模块;自适应窗口大小确定模块通过对地形描述矩阵进行分析,确定自适应窗口,以适应深空未知环境下天体表面地形判别;各纹理特征计算并融合模块利用图像四个方向上灰度共生矩阵的均值作为图像整体的纹理描述,基于此构造角二阶矩、对比度、相关性、熵和同质性五个地形判别准则从不同角度描述纹理信息;决策级融合并选取最佳着陆点模块对于各个区域分别计算图像统计特性、梯度特性得到决策指标以及各纹理特征计算并融合模块得到的特征融合得到地形判别指标向量;第二步自适应窗口大小确定模块确定自适应窗口自适应窗口大小确定模块根据图像分辨率、着陆器尺寸,通过优化全局可分度量和局部均匀度量,确定自适应窗口大小;假设图像为I<sub>M×N</sub>,窗口尺寸为l×l时计算得到的地形描述矩阵为T<sub>l</sub>={t<sup>l</sup><sub>i,j</sub>}(0≤i≤M,0≤j≤N)      (1)这样,l取不同值,可得到一系列地形描述矩阵T={T<sub>li</sub>}(i=1,2,...,num);若图像分辨率为κ,着陆器底座尺寸为M<sub>b</sub>×N<sub>b</sub>,则着陆器着陆窗口<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>d</mi><mo>=</mo><mfrac><mrow><mi>m</mi><mi>a</mi><mi>x</mi><mo>{</mo><msub><mi>M</mi><mi>b</mi></msub><mo>,</mo><msub><mi>N</mi><mi>b</mi></msub><mo>}</mo></mrow><mi>&kappa;</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000818823860000011.GIF" wi="1158" he="127" /></maths>分析满足安全着陆分析的地形数据,首先要使得在着陆窗口w中地形具有均匀性,以保证着陆器不会因地形起伏过大或局部障碍而无法着陆;其次着陆窗口之间地形具有可区分性,以保证地形矩阵可以作为着陆后有效开展科学研究的依据;定义整体可分度量whole(l)<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>w</mi><mi>h</mi><mi>o</mi><mi>l</mi><mi>e</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>M</mi><mi>N</mi></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msqrt><mrow><msup><mrow><mo>(</mo><mi>T</mi><mo>(</mo><mrow><mi>i</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>j</mi></mrow><mo>)</mo><mo>-</mo><mi>T</mi><mo>(</mo><mrow><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi></mrow><mo>)</mo><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>T</mi><mo>(</mo><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn></mrow><mo>)</mo><mo>-</mo><mi>T</mi><mo>(</mo><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000818823860000012.GIF" wi="1628" he="143" /></maths>局部均匀度量local(l,d)<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>L</mi><mi>o</mi><mi>c</mi><mi>a</mi><mi>l</mi><mrow><mo>(</mo><mi>l</mi><mo>,</mo><mi>d</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>M</mi><mi>N</mi></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><msqrt><mrow><munder><munder><mi>&Sigma;</mi><mrow><mo>|</mo><mi>m</mi><mo>-</mo><mi>i</mi><mo>|</mo><mo>&le;</mo><mi>d</mi></mrow></munder><mrow><mo>|</mo><mi>n</mi><mo>-</mo><mi>j</mi><mo>|</mo><mo>&le;</mo><mi>d</mi></mrow></munder><mi>T</mi><msup><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mi>d</mi><mn>2</mn></msup><msub><mover><mi>T</mi><mo>&OverBar;</mo></mover><mi>d</mi></msub><msup><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000818823860000021.GIF" wi="1396" he="190" /></maths>其中,l为计算地形描述矩阵采用的窗口尺寸;<img file="FDA0000818823860000022.GIF" wi="49" he="62" />为窗口内地形矩阵均值;这样,利用整体可分度量whole(l)衡量整体可分性,值越大,可分性越强;利用局部均匀度量local(l,d)衡量局部着陆窗口内的均匀性,值越小,越均匀;定义L(l,d)为地形起伏判据:L(l,d)=w<sub>1</sub>Local(l,d)‑w<sub>2</sub>Whole(l)      (5)当L(l,d)取最大值时,能够保证最好的全局可分和局部均匀;而相应的自适应窗口大小为:window_Size=l<sup>*</sup>,l<sup>*</sup>st.L(l<sup>*</sup>,d)=MA<sub>l</sub>X{L(l,d)}      (6)w<sub>1</sub>,w<sub>2</sub>分别为两种度量的贡献权重;第三步各纹理特征计算并融合模块提取纹理信息各纹理特征计算并融合模块利用图像四个方向上灰度共生矩阵的均值作为图像整体的纹理描述,基于此构造角二阶矩、对比度、相关性、熵和同质性5个地形判别准则从不同角度描述纹理信息;灰度共生矩阵描述了图像中灰度值的空间依赖性,从而刻画了邻域中灰度值的空间分布情况,反映了纹理模式内灰度的空间内在联系,从而有效地描述了纹理;根据四个方向的纹理信息,构造灰度共生矩阵<img file="FDA0000818823860000023.GIF" wi="1644" he="79" />基于灰度共生矩阵<img file="FDA0000818823860000024.GIF" wi="138" he="79" />构建特征指标对图像纹理特性进行描述,进而进行地形起伏判别;角二阶矩:<img file="FDA0000818823860000025.GIF" wi="1238" he="151" />角二阶矩是图像灰度分布均匀性的度量,由于是灰度共生矩阵元素值平方和,所以也称为能量;当灰度共生矩阵的元素分布较集中于对角线时,说明从局部区域观察图像的灰度分布是较均匀的;从图像整体来观察,纹理较粗能量较大,粗纹理含有的能量较多;反之,细纹理能量较小;当<img file="FDA0000818823860000026.GIF" wi="141" he="78" />都相等的时候,能量有最大值;能量越大,灰度分布越均匀,越利于着陆;对比度:<img file="FDA0000818823860000031.GIF" wi="1245" he="146" />纹理的对比度表示图像的清晰程度,纹理的沟纹越深则对比度值越大,越不利于着陆;相关性:<img file="FDA0000818823860000032.GIF" wi="1286" he="223" />其中<img file="FDA0000818823860000033.GIF" wi="1694" he="143" /><img file="FDA0000818823860000034.GIF" wi="550" he="143" />相关性衡量纹理在行和列方向的相似程度,值越大,区域内相关度越大,纹理越具有相似性,越利于着陆;熵:<img file="FDA0000818823860000035.GIF" wi="1270" he="140" />熵是图像信息量的度量,纹理图像中含有的信息量越大则熵值越大,越不利于着陆;同质性:<img file="FDA0000818823860000036.GIF" wi="1246" he="150" />同质性是纹理相似性的描述,同质性越大,地形起伏越小,越利于着陆;灰度共生矩阵能够很好描述小块区域基元组成的纹理;另外,位移矢量d的选择对灰度共生矩阵能否有效刻画纹理非常重要,细致纹理的分析需要较小的d值,而较大的d值更适合于粗糙的纹理描述;要满足着陆器自主着陆要求,在选择d时依赖于着陆器底座尺寸和图像分辨率信息;利用图像四个方向上灰度共生矩阵的均值作为图像整体的纹理描述,基于此构造角二阶矩、对比度、相关性、熵和同质性5个地形判别准则从不同角度描述纹理信息;对于各个区域M<sub>d</sub>(p)(p为区域中心,区域大小为d×d),图像纹理特性指标向量为F<sub>W</sub>(p);F<sub>wl</sub>(p)=[Energy,Con,Cor,Entropy,Homogeniety]      (13)I<sub>M</sub>(p)中的对比度和熵与着陆安全性能呈正相关;角二阶矩、同质性和相关性,其值与着陆安全性能呈负相关;各指标分别从纹理的粗细、清晰度、相似性和信息量的特定方面对图像纹理信息进行描述,对纹理信息的描述具有不完备和不确定性;要得到稳定有效的纹理判别,就要对各指标进行特征级融合;考虑到在着陆过程中各个纹理指标对地形判别的性能不同,设定向量Cb=[cb<sub>i</sub>],(i=0,1,…,4)对每个指标进行约束;融合后纹理判别<img file="FDA0000818823860000041.GIF" wi="68" he="78" />为<img file="FDA0000818823860000042.GIF" wi="1134" he="181" />第四步决策级融合并选取最佳着陆点模块决策级融合并选取最佳着陆点决策级融合并选取最佳着陆点模块对于各个区域M<sub>d</sub>(p)分别计算图像统计特性、梯度特性,最后通过图像统计特性、梯度特性得到的决策指标以及各纹理特征计算并融合模块得到的特征融合得到地形判别指标向量;在统计特征中,天体目标地面物质一致,因而反光属性相同;在大多情况下,地形起伏越大,灰度变化越剧烈;对于相机灰度图像,用局部标准差来描述局部地形区域内障碍分布的密集度;<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>&sigma;</mi><mo>=</mo><msqrt><mrow><munder><mo>&Sigma;</mo><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>&Element;</mo><mi>w</mi></mrow></munder><mrow><msup><mrow><mo>(</mo><mrow><msub><mi>I</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><mover><mi>I</mi><mo>&OverBar;</mo></mover></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>/</mo><mrow><mo>(</mo><mrow><mi>m</mi><mo>&times;</mo><mi>n</mi></mrow><mo>)</mo></mrow></mrow></mrow></msqrt></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msqrt><mrow><mfrac><mn>1</mn><mrow><mi>m</mi><mo>&times;</mo><mi>n</mi></mrow></mfrac><munder><mo>&Sigma;</mo><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>&Element;</mo><mi>w</mi></mrow></munder><mrow><msup><msub><mi>I</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mn>2</mn></msup></mrow><mo>-</mo><msup><mover><mi>I</mi><mo>&OverBar;</mo></mover><mn>2</mn></msup></mrow></msqrt></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000818823860000043.GIF" wi="1334" he="300" /></maths>式中:w为计算窗口,大小为m×n,;<img file="FDA0000818823860000044.GIF" wi="413" he="141" />为w中灰度均值;I<sub>i,j</sub>是相对应图像中(i,j)位置的像素灰度值;σ为整个计算子窗口的灰度标准差;σ越大,所描述的地形起伏越大,越不利于着陆;为了快速获取安全着陆候选区域,直接利用σ值加以确定;梯度特征描述图像灰度分布的特征还有梯度特征;图像中每一个像素点处x和y方向上的梯度为<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>I</mi><mi>x</mi></msub><mo>=</mo><mfrac><mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000818823860000045.GIF" wi="1245" he="125" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>I</mi><mi>y</mi></msub><mo>=</mo><mfrac><mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>17</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000818823860000046.GIF" wi="1252" he="119" /></maths>其中f(x,y)为图像上(x,y)处灰度值;定义在(x,y)处的梯度<img file="FDA0000818823860000047.GIF" wi="310" he="103" />梯度信息衡量图像灰度变化大小;地形起伏越大,梯度值越大,越不利于着陆;得到决策指标为F<sub>std</sub>(p)、F<sub>inte</sub>(p)以及经过纹理特征融合得到的<img file="FDA0000818823860000051.GIF" wi="99" he="79" />得到地形判别指标向量I<sub>M</sub>(p)<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>I</mi><mi>M</mi></msub><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mo>&lsqb;</mo><msub><mi>F</mi><mrow><mi>s</mi><mi>t</mi><mi>d</mi></mrow></msub><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>F</mi><mrow><mi>int</mi><mi>e</mi></mrow></msub><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>,</mo><msub><mover><mi>F</mi><mo>^</mo></mover><mrow><mi>w</mi><mi>l</mi></mrow></msub><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>18</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000818823860000052.GIF" wi="1180" he="78" /></maths>I<sub>M</sub>(p)中各个指标分别从图像统计信息、梯度信息和纹理信息对着陆区地形进行有效描述;但同时有各自的局限性:统计特性对图像局部噪声敏感,纹理信息在光照不均匀或存在遮挡的情况下会失真;而梯度特征值描述只能描述局部灰度变化;在着陆过程中从保证着陆器安全性能出发,对各个决策指标进行融合,最终给出统一稳定有效的联合决策;定义p<sup>*</sup><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msup><mi>p</mi><mo>*</mo></msup><mo>,</mo><mi>s</mi><mi>t</mi><mo>.</mo><mo>|</mo><mo>|</mo><msub><mi>I</mi><mi>M</mi></msub><mrow><mo>(</mo><msup><mi>p</mi><mo>*</mo></msup><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>=</mo><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>p</mi></munder><mo>|</mo><mo>|</mo><msub><mi>I</mi><mi>M</mi></msub><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>19</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000818823860000053.GIF" wi="1180" he="101" /></maths>取p<sup>*</sup>为最佳着陆点;至此,完成一种基于可见光图像的深空探测器自主着陆区选取。
地址 100038 北京市西城区木樨地北里甲51号