CN104502918B - 一种低轨卫星星载图谱关联探测方法与载荷 - Google Patents

一种低轨卫星星载图谱关联探测方法与载荷 Download PDF

Info

Publication number
CN104502918B
CN104502918B CN201410851144.8A CN201410851144A CN104502918B CN 104502918 B CN104502918 B CN 104502918B CN 201410851144 A CN201410851144 A CN 201410851144A CN 104502918 B CN104502918 B CN 104502918B
Authority
CN
China
Prior art keywords
target
moving
image
dynamic phenomenon
satellite
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
Application number
CN201410851144.8A
Other languages
English (en)
Other versions
CN104502918A (zh
Inventor
张天序
杨智慧
刘珈羽
周雨田
姚守悝
张俊青
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN201410851144.8A priority Critical patent/CN104502918B/zh
Priority to US15/106,696 priority patent/US9759605B2/en
Priority to PCT/CN2015/072669 priority patent/WO2016106954A1/zh
Publication of CN104502918A publication Critical patent/CN104502918A/zh
Application granted granted Critical
Publication of CN104502918B publication Critical patent/CN104502918B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28—Investigating the spectrum
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/02—Details
    • G01J3/0205—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows
    • G01J3/0208—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows using focussing or collimating elements, e.g. lenses or mirrors; performing aberration correction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/02—Details
    • G01J3/0289—Field-of-view determination; Aiming or pointing of a spectrometer; Adjusting alignment; Encoding angular position; Size of measurement area; Position tracking
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/02—Details
    • G01J3/0291—Housings; Spectrometer accessories; Spatial arrangement of elements, e.g. folded path arrangements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28—Investigating the spectrum
    • G01J3/2823—Imaging spectrometer
    • H—ELECTRICITY
    • H01—ELECTRIC ELEMENTS
    • H01Q—ANTENNAS, i.e. RADIO AERIALS
    • H01Q15/00—Devices for reflection, refraction, diffraction or polarisation of waves radiated from an antenna, e.g. quasi-optical devices
    • H01Q15/14—Reflecting surfaces; Equivalent structures

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Image Analysis (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)

Abstract

本发明公开了一种低轨卫星星载图谱关联探测方法与载荷。该方法包括：(1)具有基于像素偏移补偿方法的检测并跟踪动目标和动态现象；(2)对动目标和动态现象的红外光谱进行多维度特性分析，以识别动目标和动态现象。该载荷包括二维伺服转台、红外反射镜、多谱段红外光学系统、红外成像单元、宽波段红外测谱单元、数据处理单元和控制单元。本发明实现了红外成像光路和短/中/长波红外测谱光路共光轴，能同时探测动目标和动态现象的红外图像信息及其红外光谱，实现对场景中的多个动目标和动态现象进行自动检测、跟踪、测谱与识别，且识别效率高，跟踪定位精确度高。

Description

一种低轨卫星星载图谱关联探测方法与载荷
技术领域
本发明属于卫星星载遥感探测技术领域，更具体地，涉及一种低轨卫星星载图谱关联探测方法与载荷。
背景技术
低轨卫星是指运行在地球上空500km～10000km的椭圆轨道卫星，其运行轨道示意图如图1所示，不同于同步轨道卫星，低轨卫星轨道高度低，和地球间相对运动速度较快，由stk软件导出的美国MSX卫星相对地球的运动速度就在7km/s～7.59km/s之间。对临边和地球背景中的强干扰物体(如卷云、结冰湖面、雷电等)在成像过程会产生较大的帧间像素偏移，形成虚假目标，因此需要有针对于卫星运动、地球运动和物体运动补偿的算法。
一切物体都会辐射能量且其光谱曲线唯一，可利用光谱唯一性对检测出的目标进行识别。现有国内外低轨遥感载荷均没有既能探测动目标和时变对象的红外图像信息，又能探测其红外光谱的能力，如美国的MSX卫星，不能对不同背景区域下动态变化的动目标和动态现象进行实时红外光谱数据采集和遥感探测。
发明内容
针对现有技术的以上缺陷或改进需求，本发明提供了一种低轨卫星星载图谱关联探测方法与载荷，能同时探测动目标和动态现象的红外图像信息及其红外光谱，实现对场景中的多个动目标和动态现象进行自动检测、跟踪、测谱与识别，且识别效率高，跟踪定位精确度高。
为实现上述目的，本发明提供了一种低轨卫星星载图谱关联探测方法，其特征在于，包括如下步骤：
(1)检测并跟踪动目标和动态现象；进一步包括如下子步骤：
(1-1)对红外图像进行背景杂波抑制；
(1-2)将抑制后的图像进行门限分割，提取疑似动目标和动态现象；
(1-3)以疑似动目标和动态现象为处理对象，通过执行步骤(A1)～(A4)，检测出真实的动目标和动态现象，并将视场中心置于动目标或动态现象的预测位置对其进行跟踪；
(1-4)测量视场中心周围区域的红外宽光谱，得到与动目标或动态现象对应的区域光谱；
(1-5)以步骤(1-3)跟踪的动目标或动态现象为处理对象，重复执行步骤(A1)～(A4)，实现对动目标或动态现象的跟踪；
其中，所述步骤(A1)～(A4)分别为：
(A1)对多帧图像的处理对象进行坐标变换的像素偏移补偿；
(A2)对像素偏移补偿后的多帧图像进行关联，检测出真实的处理对象；
(A3)在处理对象为动目标时，计算其像面速度特征，得到其预测位置；在处理对象为动态现象时，计算其灰度特征，得到其预测位置；
(A4)定位下一帧处理对象的位置到其预测位置，将视场中心置于处理对象的预测位置；
(2)对动目标和动态现象的红外光谱进行多维度特性分析，以识别动目标和动态现象；进一步包括如下子步骤：
(2-1)定位动目标和动态现象并测量其红外光谱；
(2-2)对动目标和动态现象的红外光谱进行补偿校正；
(2-3)从校正后的动目标和动态现象的红外光谱中剔除与其对应的区域光谱，进而提取动目标和动态现象的红外光谱指纹；
(2-4)将动目标和动态现象的红外光谱指纹与典型的红外光谱库进行匹配，以识别动目标和动态现象。
优选地，所述步骤(A1)中，对相邻两帧图像的处理对象进行坐标变换的像素偏移补偿通过如下方法实现：
(B1)根据探测卫星和太阳的位置，确定地心固联坐标系到卫星坐标系的转换关系矩阵M1，确定卫星坐标系到图像坐标系的转换关系矩阵M2；
其中，地心固联坐标系中，地心为原点O，X轴指向本初子午线，Z轴指向北极点，对X轴和Z轴按照右手螺旋法得到的方向为Y轴方向；卫星坐标系中，探测卫星指向地心方向l1为Z′轴方向，探测卫星指向太阳的方向为l2，由l1和l2叉乘得到的向量方向为Y′轴方向，对Y′轴和Z′轴按照右手螺旋法则得到的方向为X′轴方向；图像中的单个像素点在图像坐标系中的坐标是该像素点所在的行数和列数，图像u坐标轴与Y′轴同方向，图像v坐标轴与Z′OX′平面平行，与X′轴和Z′轴均成45度夹角；
(B2)根据当前帧图像和前一帧图像成像时刻卫星和太阳的位置，分别得到当前帧图像对应的地心固联坐标系到卫星坐标系的转换关系矩阵M1′和卫星坐标系到图像坐标系的转换关系矩阵M2′，以及前一帧图像对应的地心固联坐标系到卫星坐标系的转换关系矩阵M1″和卫星坐标系到图像坐标系的转换关系矩阵M2″，求出同一处理对象在当前帧和前一帧像面上出现位置之间的偏差
(B3)对之前获得的多帧图像中的同一处理对象进行基于卫星和地球运动的像素偏移补偿，得到处理对象在补偿后的多帧图像中的位置，分别计算同一处理对象在相邻两帧图像间的像面位移；根据处理对象的帧间位移和帧间成像时间间隔得到处理对象在像面上的移动速度；对多帧图像得到的处理对象在像面上的移动速度进行平均，得到处理对象的平均像面速度计算得到下一帧图像中处理对象由其自身运动引起的像素偏移其中，t为相邻两帧的时间间隔。
(B4)对卫星和地球运动以及处理对象运动引起的像素偏移进行补偿，预测下一帧处理对象在图像坐标系中的位置为 ( u 1 , v 1 ) = ( u 0 + Δ u 0 → + Δ u 1 → , v 0 + Δv 0 → + Δ v 1 → ) .
优选地，所述步骤(2-4)中，对任一动目标或动态现象，若典型的红外光谱库中不存在其红外光谱指纹，则将其红外光谱指纹添加进典型的红外光谱库中；若典型的红外光谱库中存在其红外光谱指纹，则识别出动目标或动态现象，并根据识别结果做出相应的决策与处理。
按照本发明的另一方面，提供了一种低轨卫星星载图谱关联探测载荷，其特征在于，包括二维伺服转台、红外反射镜、多谱段红外光学系统、红外成像单元、宽波段红外测谱单元、数据处理单元和控制单元；所述红外反射镜置于所述二维伺服转台上，用于将动目标和动态现象辐射的红外光线反射到所述多谱段红外光学系统中；所述多谱段红外光学系统用于将红外入射光分为两路红外出射光，一路到达所述红外成像单元，进行红外成像，另一路到达所述宽波段红外测谱单元，用于对动目标和动态现象进行测谱，提取其红外光谱指纹；所述数据处理单元用于对所述红外成像单元得到的红外图像和所述宽波段红外测谱单元提取的红外光谱指纹进行实时处理，得到动目标和动态现象的预测位置以及其识别结果；所述控制单元用于根据动目标和动态现象的预测位置，控制所述二维伺服转台转动，使动目标和动态现象的影像始终落在所述载荷的视场中心，实现视场扫描和目标跟踪。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，具有以下有益效果：
(1)本发明实现了红外成像光路和短/中/长波红外测谱光路共光轴，不仅使得本发明具有同时探测动目标和动态现象的红外图像信息以及其红外光谱的能力，还减小了装备空间的体积，便于卫星的携带，具有性价比高的特点。
(2)通过基于卫星和地球运动补偿以及动目标和动态现象运动补偿的像素偏移补偿，能实现对场景中的多个动目标和动态现象的精确定位和跟踪，具有响应时间短，识别效率高的特点。
附图说明
图1是低轨卫星运行的轨道示意图；
图2是本发明实施例的低轨卫星星载图谱关联探测方法流程图；
图3是坐标系示意图；
图4是本发明实施例的低轨卫星星载图谱关联探测载荷示意图；
图5是多谱段红外光学系统的布局图；
图6是多谱段红外光学系统采用的卡塞格林系统示意图；
图7是低轨卫星星载图谱关联探测载荷的工作示意图；
图8是低轨卫星星载图谱关联探测载荷进行多目标捕获、跟踪与测谱的实例图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图2所示，本发明实施例的低轨卫星星载图谱关联探测方法包括如下步骤：
(1)检测并跟踪动目标和动态现象。进一步包括如下子步骤：
(1-1)对红外图像进行背景杂波抑制；
(1-2)将抑制后的图像进行门限分割，提取疑似动目标和动态现象；
(1-3)以疑似动目标和动态现象为处理对象，通过下述步骤(A1)～(A4)，检测出真实的动目标和动态现象，并跟踪至动目标或动态现象的预测位置。
(A1)对多帧图像的处理对象进行坐标变换的像素偏移补偿；
首先通过坐标变换求取处理对象在同一坐标系下不同时刻的位置，然后求出处理对象在红外图像上的偏移量，进行像素偏移补偿。
涉及如下三个坐标系：地心固联坐标系、卫星坐标系和图像坐标系。
地心固联坐标系：以地心为原点O，X轴指向本初子午线，Z轴指向北极点，对X轴和Z轴按照右手螺旋法得到的方向为Y轴方向。
卫星坐标系：探测卫星指向北极点的方向l1为Z′轴方向，探测卫星指向太阳的方向为l2，由l1和l2叉乘得到的向量方向为Y′轴方向，对Y′轴和Z′轴按照右手螺旋法则得到的方向为X′轴方向。
图像坐标系：卫星相机的图像以像素为单位，图像中的单个像素点在图像坐标系中的坐标是该像素点所在的行数和列数。
(A2)对像素偏移补偿后的多帧图像进行关联，检测出真实的处理对象；
(A3)在处理对象为动目标时，计算其像面速度特征，得到其预测位置；在处理对象为动态现象时，计算其灰度特征，得到其预测位置；
(A4)定位下一帧处理对象的位置到其预测位置，将视场中心置于处理对象的预测位置。
(1-4)由于此时定位精度较低，可先对动目标或动态现象对应的背景区域进行测谱，具体地，测量视场中心周围区域(如5*5像素区域)的红外宽光谱，得到与动目标或动态现象对应的区域光谱；
(1-5)以步骤(1-3)跟踪的动目标或动态现象为处理对象，重复执行上述步骤(A1)～(A4)，实现对动目标或动态现象的精确跟踪。
上述步骤(A1)中，对相邻两帧图像的处理对象进行坐标变换的像素偏移补偿通过如下方法实现：
(B1)根据探测卫星和太阳的位置，由STK软件确定地心固联坐标系到卫星坐标系的转换关系矩阵M1，确定卫星坐标系到图像坐标系的转换关系矩阵M2；
其中，地心固联坐标系中，地心为原点O，X轴指向本初子午线，Z轴指向北极点，对X轴和Z轴按照右手螺旋法得到的方向为Y轴方向；卫星坐标系中，探测卫星指向地心方向l1为Z′轴方向，探测卫星指向太阳的方向为l2，由l1和l2叉乘得到的向量方向为Y′轴方向，对Y′轴和Z′轴按照右手螺旋法则得到的方向为X′轴方向；图像中的单个像素点在图像坐标系中的坐标是该像素点所在的行数和列数，图像u坐标轴与Y′轴同方向，图像v坐标轴与Z′OX′平面平行，与X′轴和Z′轴均成45度夹角；
(B2)根据当前帧图像和前一帧图像成像时刻卫星和太阳的位置，分别得到当前帧图像对应的地心固联坐标系到卫星坐标系的转换关系矩阵M1′和卫星坐标系到图像坐标系的转换关系矩阵M2′，以及前一帧图像对应的地心固联坐标系到卫星坐标系的转换关系矩阵M1″和卫星坐标系到图像坐标系的转换关系矩阵M2″，求出同一处理对象在当前帧和前一帧像面上出现位置之间的偏差
(B3)对之前获得的多帧图像中的同一处理对象进行基于卫星和地球运动的像素偏移补偿，得到处理对象在补偿后的多帧图像中的位置，分别计算同一处理对象在相邻两帧图像间的像面位移；根据处理对象的帧间位移和帧间成像时间间隔得到处理对象在像面上的移动速度；对多帧图像得到的处理对象在像面上的移动速度进行平均，得到处理对象的平均像面速度计算得到下一帧图像中处理对象由其自身运动引起的像素偏移其中，t为相邻两帧的时间间隔；
(B4)对卫星和地球运动以及处理对象运动引起的像素偏移进行补偿，预测下一帧处理对象在图像坐标系中的位置为 ( u 1 , v 1 ) = ( u 0 + Δ u 0 → + Δ u 1 → , v 0 + Δv 0 → + Δ v 1 → ) .
(2)对动目标和动态现象的红外光谱进行多维度特性分析，以识别动目标和动态现象。进一步包括如下子步骤：
(2-1)通过3～5帧图像精确定位动目标和动态现象并测量其红外光谱；
(2-2)对动目标和动态现象的红外光谱进行补偿校正；
(2-3)从校正后的动目标和动态现象的红外光谱中剔除与其对应的区域光谱，进而提取动目标和动态现象的红外光谱指纹；
(2-4)将动目标和动态现象的红外光谱指纹与典型的红外光谱库进行匹配，以识别动目标和动态现象。
具体地，对任一动目标或动态现象，若典型的红外光谱库中不存在其红外光谱指纹，则将其红外光谱指纹添加进典型的红外光谱库中；若典型的红外光谱库中存在其红外光谱指纹，则识别出动目标或动态现象，并根据识别结果做出相应的决策与处理。
如图4所示，实现上述探测方法的低轨卫星星载图谱关联探测载荷包括二维伺服转台、红外反射镜、多谱段红外光学系统、红外成像单元、宽波段红外测谱单元、数据处理单元和控制单元。
红外反射镜置于二维伺服转台上，二维伺服转台在伺服电机和陀螺(包括方位混转测量陀螺和光纤陀螺)的控制下转动，带动红外反射镜转动，将动目标和动态现象辐射的红外光线反射到多谱段红外光学系统中，由多谱段红外光学系统得到两路红外出射光，一路到达红外成像单元，进行红外成像，另一路到达宽波段红外测谱单元，对动目标和动态现象进行测谱。
多谱段红外光学系统的布局如图5所示，它采用如图6所示的卡塞格林系统，其由一个抛物面主镜和一个双曲面次镜组成，并由若干透镜组校正像差。此系统实现了红外成像光路和短/中/长波红外测谱光路共轴，抛物面的焦点和双曲面的虚焦点重合，再经双曲面成理想像于另一个焦点。为了降低镜头本身的辐射对探测造成的影响，进行无热化设计。分光镜镀双层增透膜，使其对短、中波红外光反射率高，对长波红外光则具有半透半反的功能。
红外成像单元为焦平面阵列(FPA)，其上排列着感光元件阵列，可将发光源发射的红外线经过上述多谱段红外光学系统成像在焦平面上的这些感光元件上，然后将接收到的光信号转换为电信号并进行积分放大、采样保持，通过输出缓冲和多路传输系统，最终形成图像。
宽波段红外测谱单元用于对动目标和动态现象进行测谱，提取其红外光谱指纹。
数据处理单元可对FPA所成的红外图像和宽谱段红外测谱单元所提取的红外光谱指纹进行实时处理，具体为一方面根据红外图像检测动目标和动态现象，对动目标和动态现象的位置进行像素偏移补偿，并且通过多帧关联得到动目标的坐标和像面速度信息或者动态现象的坐标和灰度信息，进而得到动目标和动态现象的预测位置并将其发送至控制单元；另一方面根据提取出的动目标和动态现象的红外光谱指纹与典型的红外光谱库进行匹配，根据光谱匹配结果与红外图像对动目标和动态现象进行识别，并将识别结果传递给控制单元。
控制单元接收数据处理单元传输的数据，得到各动目标和动态现象的预测位置，进而控制二维伺服电机以及方位混转测量陀螺和光纤陀螺使二维伺服转台转动，带动红外反射镜实现偏转和倾斜两个维度的运动，使动目标和动态现象的影像始终落在低轨卫星星载图谱关联探测载荷的视场中心，实现视场扫描和目标跟踪。
为使本领域技术人员更好地理解本发明，下面结合具体实施例，对采用本发明的探测载荷实现上述低轨卫星星载图谱关联探测方法的过程进行详细说明。
如图7所示，按照卫星移动轨迹在固定视场内连续三帧检测提取疑似目标。为了减少图谱关联探测载荷的摆动次数，固定视场，按照卫星移动轨迹对连续三帧图像进行检测；对获取的每一帧图像进行背景杂波抑制，然后进行门限分割提取疑似动目标和动态现象，将高于门限值的像素点灰度值置为最高，将低于门限值的像素点灰度值置为0，进而提取疑似动目标和动态现象。
通过像素偏移补偿和多帧关联从所有疑似动目标和动态现象中确定真实的动目标和动态现象，按照连通域规则标记出真实的动目标和动态现象，然后计算每个动目标或动态现象的质心位置，通过坐标关联初步确定各真实动目标或动态现象的位置和像面速度信息或灰度信息，从而降低虚警。
计算所有动目标和动态现象的平均像面速度和灰度变化，判断其类型。若动目标或动态现象的位置具有不变性，且灰度特征具有时变性则可确定其为动态现象，否则就确定其为动目标；若确定为动态现象类型，则进入专注模式，即控制扫描镜中心光轴对准其位置进行长时间测谱。
因为低轨卫星运行速度较快，同时为了避免二维伺服转台转动太大的角度，那么我们可从实际情况出发，从已发现的真实动目标和动态现象中选择离视场中心最近的5个动目标或动态现象进行跟踪测谱，若不足5个，则对所有的动目标或动态现象都进行跟踪测谱，对筛选出的动目标和动态现象按照从下到上，从左到右的原则测量它们的红外光谱并进行数据处理和目标识别。如图8所示，利用第①～③帧图像进行动目标跟踪，白色框表示整个观测视场，黑色框表示动目标所处的区域。
包括如下步骤：
由红外图像检测得到动目标和动态现象的坐标和运动特征，对动目标和动态现象的位置进行卫星和地球运动补偿以及动目标和动态现象运动补偿预测定位；
根据当前动目标或动态现象的位置和像面速度信息判断转镜是否能跟踪上此动目标或动态现象。若可以跟踪上，则进入下面的步骤；若跟踪不上，则根据优先级对下一动目标或动态现象进行跟踪测谱；
控制二维伺服转台，带动红外反射镜指向此动目标或动态现象的预测坐标位置；
对当前视场中心附近区域进行3～5帧精检测、跟踪，同时测量动目标或动态现象的背景光谱，计算像素偏移补偿，控制多谱段红外光学系统的中心精对准此动目标或动态现象；
锁定当前的动目标或动态现象，并对其测谱；
对测到的谱线数据进行在线补偿矫正处理，若光谱指纹库中无此光谱数据，则添加此动目标或动态现象的红外光谱特征数据；否则，即可对光谱进行匹配，进而实现动目标或动态现象的识别与确认；
对下一动目标或动态现象重复上述步骤的操作，直至对所有筛选出来的动目标和动态现象都进行了红外光谱采集和匹配。
转动二维伺服转台，恢复到初始方向，重复上述所有操作。
实施示例为，将载荷放置在某一低轨卫星(如MSX)上，则卫星运行最大速度为7.6km/s，当前超音速飞机速度在1km/s以内，则我们的装置采用1024*1024的图像，测量光谱区域大小5*5，测谱积分时间为0.1s，则需要的分辨率大小为400m*400m，图像产生的帧率为50fps。因此，本载荷具有高性价比，高捕获率的特点。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (3)

1.一种低轨卫星星载图谱关联探测方法，其特征在于，包括如下步骤：
(1)检测并跟踪动目标和动态现象；进一步包括如下子步骤：
(1-1)对红外图像进行背景杂波抑制；
(1-2)将抑制后的图像进行门限分割，提取疑似动目标和动态现象；
(1-3)以疑似动目标和动态现象为处理对象，通过执行步骤(A1)～(A4)，检测出真实的动目标和动态现象，并将视场中心置于动目标或动态现象的预测位置对其进行跟踪；
(1-4)测量视场中心周围区域的红外宽光谱，得到与动目标或动态现象对应的区域光谱；
(1-5)以步骤(1-3)跟踪的动目标或动态现象为处理对象，重复执行步骤(A1)～(A4)，实现对动目标或动态现象的跟踪；
其中，所述步骤(A1)～(A4)分别为：
(A1)对多帧图像的处理对象进行坐标变换的像素偏移补偿；
(A2)对像素偏移补偿后的多帧图像进行关联，检测出真实的处理对象；
(A3)在处理对象为动目标时，计算其像面速度特征，得到其预测位置；在处理对象为动态现象时，计算其灰度特征，得到其预测位置；
(A4)定位下一帧处理对象的位置到其预测位置，将视场中心置于处理对象的预测位置；
(2)对动目标和动态现象的红外光谱进行多维度特性分析，以识别动目标和动态现象；进一步包括如下子步骤：
(2-1)定位动目标和动态现象并测量其红外光谱；
(2-2)对动目标和动态现象的红外光谱进行补偿校正；
(2-3)从校正后的动目标和动态现象的红外光谱中剔除与其对应的区域光谱，进而提取动目标和动态现象的红外光谱指纹；
(2-4)将动目标和动态现象的红外光谱指纹与典型的红外光谱库进行匹配，以识别动目标和动态现象。
2.如权利要求1所述的低轨卫星星载图谱关联探测方法，其特征在于，所述步骤(A1)中，对相邻两帧图像的处理对象进行坐标变换的像素偏移补偿通过如下方法实现：
(B1)根据探测卫星和太阳的位置，确定地心固联坐标系到卫星坐标系的转换关系矩阵M1，确定卫星坐标系到图像坐标系的转换关系矩阵M2；
其中，地心固联坐标系中，地心为原点O，X轴指向本初子午线，Z轴指向北极点，对X轴和Z轴按照右手螺旋法得到的方向为Y轴方向；卫星坐标系中，探测卫星指向地心方向l1为Z′轴方向，探测卫星指向太阳的方向为l2，由l1和l2叉乘得到的向量方向为Y′轴方向，对Y′轴和Z′轴按照右手螺旋法则得到的方向为X′轴方向；图像中的单个像素点在图像坐标系中的坐标是该像素点所在的行数和列数，图像u坐标轴与Y′轴同方向，图像v坐标轴与Z′OX′平面平行，与X′轴和Z′轴均成45度夹角；
(B2)根据当前帧图像和前一帧图像成像时刻卫星和太阳的位置，分别得到当前帧图像对应的地心固联坐标系到卫星坐标系的转换关系矩阵M1′和卫星坐标系到图像坐标系的转换关系矩阵M2′，以及前一帧图像对应的地心固联坐标系到卫星坐标系的转换关系矩阵M1″和卫星坐标系到图像坐标系的转换关系矩阵M2″，求出同一处理对象在当前帧和前一帧像面上出现位置之间的偏差
(B3)对之前获得的多帧图像中的同一处理对象进行基于卫星和地球运动的像素偏移补偿，得到处理对象在补偿后的多帧图像中的位置，分别计算同一处理对象在相邻两帧图像间的像面位移；根据处理对象的帧间位移和帧间成像时间间隔得到处理对象在像面上的移动速度；对多帧图像得到的处理对象在像面上的移动速度进行平均，得到处理对象的平均像面速度计算得到下一帧图像中处理对象由其自身运动引起的像素偏移其中，t为相邻两帧的时间间隔；
(B4)对卫星和地球运动以及处理对象运动引起的像素偏移进行补偿，预测下一帧处理对象在图像坐标系中的位置为
3.如权利要求1或2所述的低轨卫星星载图谱关联探测方法，其特征在于，所述步骤(2-4)中，对任一动目标或动态现象，若典型的红外光谱库中不存在其红外光谱指纹，则将其红外光谱指纹添加进典型的红外光谱库中；若典型的红外光谱库中存在其红外光谱指纹，则识别出动目标或动态现象，并根据识别结果做出相应的决策与处理。
CN201410851144.8A 2014-12-30 2014-12-30 一种低轨卫星星载图谱关联探测方法与载荷 Active CN104502918B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201410851144.8A CN104502918B (zh) 2014-12-30 2014-12-30 一种低轨卫星星载图谱关联探测方法与载荷
US15/106,696 US9759605B2 (en) 2014-12-30 2015-02-10 Low-orbit satellite-borne image-spectrum associated detection method and payload
PCT/CN2015/072669 WO2016106954A1 (zh) 2014-12-30 2015-02-10 一种低轨卫星星载图谱关联探测方法与载荷

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410851144.8A CN104502918B (zh) 2014-12-30 2014-12-30 一种低轨卫星星载图谱关联探测方法与载荷

Publications (2)

Publication Number Publication Date
CN104502918A CN104502918A (zh) 2015-04-08
CN104502918B true CN104502918B (zh) 2016-09-07

Family

ID=52944329

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410851144.8A Active CN104502918B (zh) 2014-12-30 2014-12-30 一种低轨卫星星载图谱关联探测方法与载荷

Country Status (3)

Country Link
US (1) US9759605B2 (zh)
CN (1) CN104502918B (zh)
WO (1) WO2016106954A1 (zh)

Families Citing this family (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105628201A (zh) * 2015-12-30 2016-06-01 华中科技大学 一种时敏目标图谱结合的测谱方法
CN107272746A (zh) * 2017-06-30 2017-10-20 中国科学院长春光学精密机械与物理研究所 一种中低轨卫星的跟踪方法
CN107666262B (zh) * 2017-08-30 2024-06-18 中国科学院上海技术物理研究所 一种星载二维指向机构永磁同步电机极点检测装置及方法
CN107886498B (zh) * 2017-10-13 2021-04-13 中国科学院上海技术物理研究所 一种基于星载图像序列的空间目标检测跟踪方法
WO2019083533A1 (en) * 2017-10-26 2019-05-02 Bae Systems Information And Electronic Systems Integration Inc. SYSTEM AND METHOD FOR DETECTING SUBSTANCES FROM SPECTRAL SIGNATURES
CN108305290B (zh) * 2017-12-29 2020-12-08 华中科技大学 一种动目标的精确测谱方法
CN108804621A (zh) * 2018-05-31 2018-11-13 四川华控图形科技有限公司 装备标准的关联图谱构建方法
CN109087378A (zh) * 2018-09-11 2018-12-25 首都师范大学 图像处理方法及系统
CN109255353B (zh) * 2018-09-12 2022-06-28 首都师范大学 一种动目标检测方法、装置、电子设备及存储介质
CN112953632A (zh) * 2018-12-29 2021-06-11 长沙天仪空间科技研究院有限公司 一种激光捕获系统及方法
CN109993768B (zh) * 2018-12-29 2020-10-30 华中科技大学 一种提高伺服跟踪实时性和精准性的空中目标测谱方法
CN109828234B (zh) * 2019-01-18 2020-12-15 中国久远高新技术装备公司 一种异源光学平台对空间目标融合定位方法
CN109800379B (zh) * 2019-01-30 2023-03-24 上海卫星工程研究所 星载微波遥感仪器光路建模方法
CN109886993B (zh) * 2019-02-27 2022-12-30 中国科学院上海技术物理研究所 一种地球静止轨道高帧频大动态范围动目标探测方法
CN110780327B (zh) * 2019-10-29 2022-04-08 中国人民解放军军事科学院国防科技创新研究院 基于星载ais与红外相机的海上目标协同定位方法
CN111161304B (zh) * 2019-12-16 2023-06-30 北京空间机电研究所 一种快速背景估计的遥感视频目标轨迹跟踪方法
WO2021163040A1 (en) * 2020-02-11 2021-08-19 Ast & Science, Llc Geolocation of radio frequency devices using spaceborne phased arrays
CN111708381A (zh) * 2020-06-29 2020-09-25 西安方元明科技股份有限公司 一种双轴跟瞄转台跟踪精度测试方法
CN111913190B (zh) * 2020-07-24 2022-07-29 北京控制工程研究所 一种基于彩色红外光谱共孔径成像的临近空间暗弱目标定向装置
CN112330739B (zh) * 2020-10-27 2022-05-27 北京航空航天大学 一种基于光学概率图模型的卫星探测方法
CN112960145B (zh) * 2021-02-04 2022-12-13 上海卫星工程研究所 一种遥感卫星对地姿态机动扫描的轨迹规划方法及系统
CN113063437B (zh) * 2021-02-26 2023-03-31 上海卫星工程研究所 适用于全物理仿真的卫星视线指向模拟方法及系统
CN113038498B (zh) * 2021-03-02 2022-09-23 航天科工空间工程发展有限公司 一种通过低轨卫星波位设计实现波束对地覆盖的方法
CN114004862A (zh) * 2021-10-12 2022-02-01 北京航天长征飞行器研究所 一种长波红外弱小目标检测识别方法
CN114166357B (zh) * 2021-11-02 2024-07-02 中航洛阳光电技术有限公司 一种红外探测器横滚角测量方法
CN114679541B (zh) * 2022-03-11 2024-06-18 中国科学院长春光学精密机械与物理研究所 一种星上运动目标跟踪方法
CN114694064B (zh) * 2022-03-23 2023-05-02 成都智元汇信息技术股份有限公司 一种基于目标识别的切图方法与系统
CN114792392B (zh) * 2022-05-18 2023-04-18 西南交通大学 一种基于马尔科夫链的雷电先导发展路径的模拟方法
CN116500779B (zh) * 2023-04-27 2024-04-26 中国科学院长春光学精密机械与物理研究所 基于天基平台与周转旋镜联动的高频宽幅成像方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7106511B2 (en) * 2004-03-02 2006-09-12 Agilent Technologies, Inc. Imaging system with large depth of field
EP2051051B1 (en) 2007-10-16 2020-06-03 Cambridge Research & Instrumentation, Inc. Spectral imaging system with dynamic optical correction
JP2009139352A (ja) 2007-12-11 2009-06-25 Shimadzu Corp フーリエ変換型赤外分光光度計
CN101303291B (zh) 2008-06-05 2010-06-30 苏州大学 基于数字微镜器件的多通道多目标超光谱成像方法及系统
CN102564589B (zh) * 2011-12-20 2013-07-24 华中科技大学 一种多波段动目标光谱特征探测识别方法和装置
US20140125811A1 (en) * 2012-11-08 2014-05-08 Lockheed Martin Corporation System, method and computer program product for identifying chemical composition of optically lifted latent prints and contaminants
CN103196862B (zh) * 2013-02-25 2015-01-21 北京师范大学 基于ASAR和Hyperion数据反演植被覆盖下土壤水分的方法及系统
CN103776540B (zh) * 2013-12-30 2015-07-08 华中科技大学 一种多波段共光路图谱联合遥感测量系统及方法
CN103777348B (zh) 2013-12-30 2016-01-13 华中科技大学 一种多波段灵巧红外光学系统

Also Published As

Publication number Publication date
US20160356650A1 (en) 2016-12-08
US9759605B2 (en) 2017-09-12
WO2016106954A1 (zh) 2016-07-07
CN104502918A (zh) 2015-04-08

Similar Documents

Publication Publication Date Title
CN104502918B (zh) 一种低轨卫星星载图谱关联探测方法与载荷
US20230148273A1 (en) Coded localization systems, methods and apparatus
US11778289B2 (en) Multi-camera imaging systems
Nouri et al. Cloud height and tracking accuracy of three all sky imager systems for individual clouds
CN105676305B (zh) 一种共口径多视场图谱协同探测系统与方法
CN104501956A (zh) 一种超宽波段图谱关联探测装置与探测方法
CN204439211U (zh) 一种超宽波段图谱关联探测装置
CN102927982B (zh) 双光谱的自主导航敏感器及其设计方法
CN105182678B (zh) 一种基于多通道相机观测空间目标的系统及方法
Beekmans et al. Cloud photogrammetry with dense stereo for fisheye cameras
CN108830811A (zh) 一种飞行参数与相机内参相结合的航空影像实时校正方法
CN103129752A (zh) 一种基于地面导航的光学遥感卫星姿态角误差动态补偿方法
Ofek et al. The large array survey telescope—system overview and performances
CN105181137A (zh) 用于地基对月观测的宽波段高光谱分辨率成像系统
CN106569206A (zh) 一种基于微波光学复合的目标探测方法
CN112924028A (zh) 一种面向海面溢油的光场偏振成像探测系统
CN103234552A (zh) 光学导航目标星模拟仿真图像生成方法
US20200096388A1 (en) Cubesat infrared atmospheric sounder (ciras)
CN116343057A (zh) 一种联合sar与光学图像的舰船目标检测与识别方法
CN110111294A (zh) 一种sar图像质量检测方法
CN109241554A (zh) 一种临边全波段超光谱图像仿真方法
US11994365B2 (en) Projectile tracking and 3D traceback method
Serres et al. Passive Polarized Vision for Autonomous Vehicles: A Review
CN115933170A (zh) 适用于高光谱视频成像的大通量编码孔径组件
Cavanaugh Electro-Optical and Infrared Design for Uncrewed Aerial System Collision Avoidance

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant