CN107643074A - 一种航空扫描仪摆扫成像方位预置方法 - Google Patents
一种航空扫描仪摆扫成像方位预置方法 Download PDFInfo
- Publication number
- CN107643074A CN107643074A CN201710801778.6A CN201710801778A CN107643074A CN 107643074 A CN107643074 A CN 107643074A CN 201710801778 A CN201710801778 A CN 201710801778A CN 107643074 A CN107643074 A CN 107643074A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- orientation
- mover
- 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.)
- Granted
Links
Landscapes
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Facsimile Scanning Arrangements (AREA)
Abstract
本发明属于扫描成像技术领域,具体涉及一种航空扫描仪摆扫成像方位预置方法。该方法通过补偿载机和目标相对运动引入的像移角速度,保证扫描仪成像器视轴从初始状态转移到期望状态,并且使扫描状态和预置状态切换过程中角速度超调量尽可能小。
Description
技术领域
本发明属于扫描成像技术领域,具体涉及一种航空扫描仪摆扫成像方位预置方法。
背景技术
装有线列TDI相机的航空平台扫描仪采用双向摆扫模式对地面景物进行成像,摆扫成像示意图,如图1所示。成像涉及多个过程,包括扫描仪视线(LOS)指向期望的角度位置,在适当的区域位置扫描,反向扫描时扫描仪视线指向期望的角度位置,最终将包含目标的多帧图像拼接成一幅场景图像,基本成像过程见图2。方位预置就是对扫描仪视线重定位的指向过程,它的主要任务是在接收到摆扫成像指令后,根据输入信息分析计算视线指向的角度、角速度和角加速度指令,然后发给伺服系统,驱动框架系统完成扫描仪视线的重定位。
发明内容
(一)要解决的技术问题
本发明提出一种航空扫描仪摆扫成像方位预置方法,以解决如何提高航空机载高分辨率扫描仪在扫描成像过程中视线的指向精确性的技术问题。
(二)技术方案
为了解决上述技术问题,本发明提出一种航空扫描仪摆扫成像方位预置方法,该方法包括如下步骤:
S1、采集信息数据,设定方位预置时间Tpre,并将方位预置时间Tpre分为前半段方位预置时间t1和后半段方位预置时间t2;其中,
信息数据包括方位预置目标框架角速度ωtarget、方位预置初始框架角速度ωinit、方位预置目标框架角Starget、方位预置初始框架角Sinit,以及载机和目标相对运动的像移补偿角速度ωcompensate;
前半段方位预置时间t1、后半段方位预置时间t2与与方位预置时间Tpre的关系,如公式(1):
S2、根据公式(2)和(3),分别计算前半段方位预置时间t1和后半段方位预置时间t2两个时间段的加速度微分值和
S3、根据公式(4)和(5),分别计算前半段方位预置时间t1和后半段方位预置时间t2两个时间段的角速度ω1和ω2:
S4、根据公式(6)和(7),分别计算前半段方位预置时间t1和后半段方位预置时间t2两个时间段的视轴方位预置指令角速度ωcmd1和ωcmd2:
ωcmd1=ω1+ωcompensate(0≤t≤t1) (6),
ωcmd2=ω2+ωcompensate(t1≤t≤Tpre) (7);
S5、将视轴方位预置指令角速度ωcmd1和ωcmd2作为方位预置指令发送至扫描仪的伺服控制单元。
(三)有益效果
本发明提出的航空扫描仪摆扫成像方位预置方法,通过补偿载机和目标相对运动引入的像移角速度,保证扫描仪成像器视轴从初始状态转移到期望状态,并且使扫描状态和预置状态切换过程中角速度超调量尽可能小。
附图说明
图1为现有技术中摆扫成像示意图;
图2为现有技术中基本成像过程示意图:①LOS指向初始位置,②整个条带扫描,③指向下一次反向扫描的起始位置;
图3为本发明实施例的方位预置指令图和局部细节图。
具体实施方式
为使本发明的目的、内容和优点更加清楚,下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。
本实施例提出一种航空扫描仪摆扫成像方位预置方法,该方法包括如下步骤:
S1、采集信息数据,设定方位预置时间Tpre,并将方位预置时间Tpre分为前半段方位预置时间t1和后半段方位预置时间t2;其中,
信息数据包括方位预置目标框架角速度ωtarget、方位预置初始框架角速度ωinit、方位预置目标框架角Starget、方位预置初始框架角Sinit,以及载机和目标相对运动的像移补偿角速度ωcompensate。
视轴方位预置指令对框架回转运动采用时间最优策略,预置时间二等分,即前半段方位预置时间t1、后半段方位预置时间t2与与方位预置时间Tpre的关系,如公式(1):
S2、根据公式(8)和(9),分别计算前半段方位预置时间t1和后半段方位预置时间t2两个时间段的加速度微分值和
框架回转过程是框架在初始时刻,以初始框架角速度ωinit由初始位置Sinit处开始运动,经过时间Tpre后到达Starget处,此时框架角速度为ωtarget。运动过程分为两个阶段t1和t2,两个阶段内回转运动均是匀变加速运动,其中加速度以及加速度微分分别用符号a1、a2、和表示。
这个运动过程中,由角加速度关系可以得到:
由角速度关系可以得到:
符号ω1和ω2分别为t1和t2时间段的角速度,由运动关系分解得:
由角位置关系可以得到
由关系式(1)、(2)、(3)、(4)和(7)可解得:
S3、由关系式(5)和(6),分别计算前半段方位预置时间t1和后半段方位预置时间t2两个时间段的角速度ω1和ω2,得到公式(10)和(11):
S4、根据公式(12)和(13),补偿载机和目标相对运动引起得像移后,分别计算前半段方位预置时间t1和后半段方位预置时间t2两个时间段的视轴方位预置指令角速度ωcmd1和ωcmd2:
ωcmd1=ω1+ωcompensate(0≤t≤t1) (12)
ωcmd2=ω2+ωcompensate(t1≤t≤Tpre) (13)
S5、将视轴方位预置指令角速度ωcmd1和ωcmd2作为方位预置指令发送至扫描仪的伺服控制单元。该指令及其局部细节,如图3所示。
本实施例的航空扫描仪摆扫成像方位预置方法已得到工程应用,经航空机载试验证明方法可以有效解决成像器视轴的指向问题,有效减小预置过程中角速度的超调量,并且可以对载机和目标相对运动引起的像移做出补偿。该方法中,信息数据源多是来自于扫描仪内部的传感器采样,计算量小、计算时间短,满足实时性要求;该方法适用范围广,可以普遍应用于航空平台的摆扫成像系统。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (1)
1.一种航空扫描仪摆扫成像方位预置方法,其特征在于,所述方法包括如下步骤:
S1、采集信息数据,设定方位预置时间Tpre,并将所述方位预置时间Tpre分为前半段方位预置时间t1和后半段方位预置时间t2;其中,
所述信息数据包括方位预置目标框架角速度ωtarget、方位预置初始框架角速度ωinit、方位预置目标框架角Starget、方位预置初始框架角Sinit,以及载机和目标相对运动的像移补偿角速度ωcompensate;
所述前半段方位预置时间t1、所述后半段方位预置时间t2与所述与方位预置时间Tpre的关系,如公式(1):
<mrow>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>T</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
S2、根据公式(2)和(3),分别计算所述前半段方位预置时间t1和所述后半段方位预置时间t2两个时间段的加速度微分值和
<mrow>
<msub>
<mover>
<mi>a</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>3</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>t</mi>
<mi>arg</mi>
<mi>e</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>5</mn>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mfrac>
<msub>
<mi>T</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>e</mi>
</mrow>
</msub>
<mn>2</mn>
</mfrac>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>t</mi>
<mi>arg</mi>
<mi>e</mi>
<mi>t</mi>
</mrow>
</msub>
<mfrac>
<msub>
<mi>T</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>e</mi>
</mrow>
</msub>
<mn>2</mn>
</mfrac>
</mrow>
<msup>
<mfrac>
<msub>
<mi>T</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>e</mi>
</mrow>
</msub>
<mn>4</mn>
</mfrac>
<mn>3</mn>
</msup>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mover>
<mi>a</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>11</mn>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mfrac>
<msub>
<mi>T</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>e</mi>
</mrow>
</msub>
<mn>2</mn>
</mfrac>
<mo>+</mo>
<mn>7</mn>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>t</mi>
<mi>arg</mi>
<mi>e</mi>
<mi>t</mi>
</mrow>
</msub>
<mfrac>
<msub>
<mi>T</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>e</mi>
</mrow>
</msub>
<mn>2</mn>
</mfrac>
<mo>-</mo>
<mn>9</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>t</mi>
<mi>arg</mi>
<mi>e</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msup>
<mfrac>
<msub>
<mi>T</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>e</mi>
</mrow>
</msub>
<mn>4</mn>
</mfrac>
<mn>3</mn>
</msup>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
S3、根据公式(4)和(5),分别计算所述前半段方位预置时间t1和所述后半段方位预置时间t2两个时间段的角速度ω1和ω2:
<mrow>
<msub>
<mi>&omega;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mover>
<mi>a</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>&le;</mo>
<mi>t</mi>
<mo>&le;</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>&omega;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>a</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>a</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<msup>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>a</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>a</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mi>t</mi>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mover>
<mi>a</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
<mo>,</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>&le;</mo>
<mi>t</mi>
<mo>&le;</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
S4、根据公式(6)和(7),分别计算所述前半段方位预置时间t1和所述后半段方位预置时间t2两个时间段的视轴方位预置指令角速度ωcmd1和ωcmd2:
ωcmd1=ω1+ωcompensate(0≤t≤t1) (6),
ωcmd2=ω2+ωcompensate(t1≤t≤Tpre) (7);
S5、将所述视轴方位预置指令角速度ωcmd1和ωcmd2作为方位预置指令发送至所述扫描仪的伺服控制单元。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710801778.6A CN107643074B (zh) | 2017-09-07 | 2017-09-07 | 一种航空扫描仪摆扫成像方位预置方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710801778.6A CN107643074B (zh) | 2017-09-07 | 2017-09-07 | 一种航空扫描仪摆扫成像方位预置方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107643074A true CN107643074A (zh) | 2018-01-30 |
CN107643074B CN107643074B (zh) | 2019-12-03 |
Family
ID=61110430
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710801778.6A Active CN107643074B (zh) | 2017-09-07 | 2017-09-07 | 一种航空扫描仪摆扫成像方位预置方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107643074B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111337945A (zh) * | 2020-04-14 | 2020-06-26 | 江苏域盾成鹫科技装备制造有限公司 | 一种基于螺旋线扫描的高分辨率制冷光电雷达 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103693210A (zh) * | 2013-12-20 | 2014-04-02 | 河北汉光重工有限责任公司 | 一种机载光电吊舱视轴稳定的建模方法 |
CN105444780A (zh) * | 2015-11-30 | 2016-03-30 | 上海卫星工程研究所 | 星载摆扫式光学相机图像定位验证系统及处理方法 |
CN106441591A (zh) * | 2016-09-28 | 2017-02-22 | 南京航空航天大学 | 一种星载高分辨率热红外影像几何定位方法 |
CN106444843A (zh) * | 2016-12-07 | 2017-02-22 | 北京奇虎科技有限公司 | 无人机相对方位控制方法及装置 |
CN106527457A (zh) * | 2016-11-17 | 2017-03-22 | 天津津航技术物理研究所 | 航空扫描仪扫描控制指令规划方法 |
-
2017
- 2017-09-07 CN CN201710801778.6A patent/CN107643074B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103693210A (zh) * | 2013-12-20 | 2014-04-02 | 河北汉光重工有限责任公司 | 一种机载光电吊舱视轴稳定的建模方法 |
CN105444780A (zh) * | 2015-11-30 | 2016-03-30 | 上海卫星工程研究所 | 星载摆扫式光学相机图像定位验证系统及处理方法 |
CN106441591A (zh) * | 2016-09-28 | 2017-02-22 | 南京航空航天大学 | 一种星载高分辨率热红外影像几何定位方法 |
CN106527457A (zh) * | 2016-11-17 | 2017-03-22 | 天津津航技术物理研究所 | 航空扫描仪扫描控制指令规划方法 |
CN106444843A (zh) * | 2016-12-07 | 2017-02-22 | 北京奇虎科技有限公司 | 无人机相对方位控制方法及装置 |
Non-Patent Citations (1)
Title |
---|
夏中秋: ""空间相机摆扫成像立体定位精度仿真分析"", 《中国空间科学技术》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111337945A (zh) * | 2020-04-14 | 2020-06-26 | 江苏域盾成鹫科技装备制造有限公司 | 一种基于螺旋线扫描的高分辨率制冷光电雷达 |
CN111337945B (zh) * | 2020-04-14 | 2022-05-03 | 江苏域盾成鹫科技装备制造有限公司 | 一种基于螺旋线扫描的高分辨率制冷光电雷达 |
Also Published As
Publication number | Publication date |
---|---|
CN107643074B (zh) | 2019-12-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110849362B (zh) | 一种基于车载惯性的激光雷达与视觉组合导航算法 | |
CN104574332B (zh) | 一种机载光电吊舱用图像融合方法 | |
EP3005292B1 (de) | 3d-aufnahmevorrichtung, verfahren zur erstellung eines 3d-bildes und verfahren zur einrichtung einer 3d-aufnahmevorrichtung | |
US20110013021A1 (en) | Image processing device and method, driving support system, and vehicle | |
CN102761708B (zh) | 一种线阵ccd图像扫描方法 | |
CN101710224B (zh) | 一种用于空间光学相机在轨自适应调焦的方法 | |
JP2017017696A (ja) | ワッブル型の歪みの補正を伴う無人航空機用の高解像度カメラ | |
CA3065207C (en) | Aircraft positioning on a taxiway | |
JP6381981B2 (ja) | 線路空間支障物検知システム | |
CN102450005A (zh) | 校正目标检测装置、检测校正目标的校正目标检测方法以及用于校正目标检测装置的程序 | |
US9154717B2 (en) | Aerial camera system and method for correcting distortions in an aerial photograph | |
CN106527457B (zh) | 航空扫描仪扫描控制指令规划方法 | |
CN112819711B (zh) | 一种基于单目视觉的利用道路车道线的车辆反向定位方法 | |
CN106709944A (zh) | 卫星遥感图像配准方法 | |
CN104167001B (zh) | 基于正交补偿的大视场摄像机标定方法 | |
CN109074480B (zh) | 用于检测机动车辆的环境区域的图像中的滚动快门效应的方法、计算装置、驾驶员辅助系统以及机动车辆 | |
CN106525002A (zh) | 一种tdiccd像移检测及补偿方法 | |
CN106340045A (zh) | 三维人脸重建中基于双目立体视觉的标定优化方法 | |
CN113264201A (zh) | 一种敏捷卫星主动推扫成像实现方法 | |
US20230308757A1 (en) | Imaging support device, imaging apparatus, imaging system, imaging support system, imaging support method, and program | |
CN105606123A (zh) | 一种低空航空摄影测量自动纠正数字地面高程模型的方法 | |
CN107643074A (zh) | 一种航空扫描仪摆扫成像方位预置方法 | |
CN103248838B (zh) | 基于数字域tdi算法的cmos相机图像调偏流方法 | |
CN105225211A (zh) | 一种基于振动提高视频分辨率的方法 | |
CN103196431A (zh) | 机载激光扫描点云与光学影像的整体空三方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |