CN102252649A - 航空遥感器目标指向定位角的计算方法 - Google Patents
航空遥感器目标指向定位角的计算方法 Download PDFInfo
- Publication number
- CN102252649A CN102252649A CN 201110095512 CN201110095512A CN102252649A CN 102252649 A CN102252649 A CN 102252649A CN 201110095512 CN201110095512 CN 201110095512 CN 201110095512 A CN201110095512 A CN 201110095512A CN 102252649 A CN102252649 A CN 102252649A
- Authority
- CN
- China
- Prior art keywords
- angle
- formula
- data
- sublist
- look
- 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
Images
Landscapes
- Image Analysis (AREA)
Abstract
航空遥感器成像时目标指向定位角的计算方法,涉及一种目标指向定位角的计算方法,它解决解决因正切函数具有单调递增特性而难以用传统查表方式获得高精度的目标指向定位角的问题,该方法首先通过坐标变换的方式将待拍照目标从地理坐标系转换到载机坐标系,然后在载机坐标系采用反正切函数计算的方法求得目标指向角。为实现反正切函数计算结果的快速性和精确性,提出了一种快速查表方式,根据精度要求按角度均匀递增的顺序将各个角度对应的正切值排列在表内建立多个有序表,采用顺序查找的方式和二分法查找方法获得待拍照目标的指向定位角。本发明减少了计算时间,最终实现了航空遥感器目标指向定位角计算的快速性和精确性。
Description
技术领域
本发明涉及一种目标指向定位角的计算方法,具体涉及采用航空遥感器时待拍摄目标指向定位角的计算方法。
背景技术
航空遥感器对地面预定目标拍照时,需首先根据航空遥感器所在位置的地理坐标、待拍照目标地理坐标等数据信息计算出航空遥感器伺服指向机构的指向定位角,然后控制指向机构运动到计算出的指向角度位置,最终实现对地面预定目标成像。
为保证成像的快速性及成像位置的准确性,要求在进行目标指向定位角计算时,计算结果精度高且计算速度快。在进行指向定位角的计算过程中,需执行反正切三角函数的运算,当航空遥感器中采用的微处理器工作频率较低时,直接采用高级语言中提供的库函数进行运算的方式,其计算时间较长;利用泰勒公式展开的方式虽然可将复杂的三角函数运算转化为乘除运算,进而适当缩短计算时间,但由于其计算精度同乘除法的运算次数成反比,为减少运算次数而减少展开级数的同时,也降低了运算精度,因此,航空遥感器在进行快速反正切函数运算时,一般采用查表的工作方式。
由于正切函数具有单调递增特性,且随着角度增大,其数值递增程度越大,当角度为90°时,其对应的正切值增加到正无穷大。因此采用查表方式进行正切函数的逆运算——反正切函数的计算时,若采用传统的查表方式,按正切值成比例递增的原则设置表格结构,在表中存放角度值,则难以根据角度对应的正切值查找到均匀分布的角度,进而导致目标指向定位角计算结果的精度难以满足系统要求。
发明内容
本发明为解决因正切函数具有单调递增特性而难以用传统查表方式获得高精度的目标指向定位角的问题,提供一种航空遥感器目标指向定位角的计算方法。
航空遥感器目标指向定位角的计算方法,该方法由以下步骤实现:
步骤一、采用坐标平移的方法将待拍摄目标在地理坐标系下的坐标转换到航空遥感器所在的载机坐标系下的坐标,所述待拍摄目标在载机坐标系下的坐标计算公式用公式一表示为:
式中:γ为航向角,θ为俯仰角,为滚动角;(X0,Y0,Z0)为航空遥感器所在位置的地理坐标;(x,y,z)为待拍摄目标在地理坐标系下的坐标;(x′,y′,z′)为待拍摄目标在载机坐标系下的坐标;
步骤二、根据步骤一所述的待拍摄目标在载机坐标系下的坐标获得目标指向定位角,所述目标指向定位角用公式二和公式三表示为:
公式二:
公式三:
式中:α为垂直飞行方向的旋转角;β为沿飞行方向的旋转角;
所述目标指向定位角通过逆向查表结构和二分查找法获得,具体步骤为:
步骤a:设置总表结构,将待查数据的正切值按角度均匀递增的顺序放入总表中,将所述总表分成多个子表,确定每个子表的表首地址;
步骤b:采用顺序查找法确定待查数据所属的子表在总表中的偏移量地址;然后采用二分法进行查找,获得与待查数据相等或相近的角度的正切值在子表中的位置;
步骤c:根据步骤b所述子表在总表中的偏移量地址、与待查数据相同的数据正切值在子表中的地址及总表中设置的角度间隔计算待查数据的角度,所述待查数据的角度为目标指向定位角。
本发明的有益效果:本发明在目标指向定位角计算过程中,反正切函数的计算采用了逆向查表模式,同时在查表过程中采用了二分法,其平均查表次数为log2n(n为表长),同平均查找法的查找次数n/2相比,减少了计算时间,实现了航空遥感器目标指向定位角计算的快速性和精确性。
附图说明
图1为本发明所述的航空遥感器目标指向定位角的计算方法中确定待查数据值A所在子表的流程图;
图2为本发明所述的航空遥感器目标指向定位角的计算方法中采用二分法确定待查数据值A在子表中的相对位置的流程图。
具体实施方式
具体实施方式一、航空遥感器目标指向定位角的计算方法,该方法由以下步骤实现:
步骤一、采用坐标平移的方法将待拍摄目标在地理坐标系下的坐标转换到航空遥感器所在的载机坐标系下的坐标,所述待拍摄目标在载机坐标系下的坐标计算公式用公式一表示为:
步骤二、根据步骤一所述的待拍摄目标在载机坐标系下的坐标获得目标指向定位角,所述目标指向定位角用公式二和公式三表示为:
公式二:
公式三:
式中:α为垂直飞行方向的旋转角;β为沿飞行方向的旋转角;
所述目标指向定位角通过逆向查表结构和二分查找法获得,具体步骤为:
步骤a:设置总表结构,将待查数据的正切值按角度均匀递增的顺序放入总表中,将所述总表分成多个子表,确定每个子表的表首地址;
步骤b:采用顺序查找法确定待查数据所属的子表在总表中的偏移量地址;然后采用二分法进行查找,获得与待查数据相等或相近的角度的正切值在子表中的位置;
步骤c:根据步骤b所述子表在总表中的偏移量地址、与待查数据相同的数据正切值在子表中的地址及总表中设置的角度间隔计算待查数据的角度,所述待查数据的角度为目标指向定位角。
本实施方式所述子表中间地址的计算公式用公式四表示:
公式四:Addr中=Addr表首+L/2
式中L表示表长,计算公式用公式五表示:
公式五:L=Addr表尾-Addr表首+1。
本实施方式中步骤a所述的将待查数据的正切值按角度均匀递增的顺序放入表中,设置多个子表的具体原则为:以待查数据的正切值占用1~2个字节为原则,并将待查数据的正切值按不同的放大倍数放入不同的子表中,相同子表中的数据正切值的放大倍数相同,所述相同子表中的数据正切值的高字节相同。
本实施方式中步骤c所述的待查数据角度的计算方法为:子表在总表中的偏移量地址与与待查数据相同的数据正切值在子表中的相对位置求和后与总表中设置的角度间隔相除。
具体实施方式二、结合图1和图2说明本实施方式,本实施方式为具体实施方式一所述的航空遥感器成像时目标指向定位角的计算方法的实施例:
步骤a1、根据
计算出待拍摄目标在载机坐标系下的坐标(x′,y′,z′)后,根据公式二和公式三
公式二:
公式三:
设目标指向定位角的正切值为A,为实现计算结果的快速性和准确性,本实施例中提出的反正切函数计算方法为:
1、按目标指向定位角的计算精度要求设置查找总表的结构,将角度对应的正切值按角度递增的顺序放入总表内,根据总表的排列情况可设置多个子表,并确定各个子表的表首地址;
2、采用顺序查找法初步确定A对应的角度所在的子表地址;
3、采用二分法进行精确查找,找到同A值相等或相近的值B在子表中的对应位置;
4、根据子表在总表中的偏移量地址、B值在子表中的位置及总表设置的角度间隔计算出A值所对应的角度,即为待求的目标指向定位角。
本实施例在设置查找表时,考虑到正切函数递增较快,若将整个角度范围内的正切值按照相同的放大倍数进行放大,则大角度对应的正切值难以用1~2个字节来表示,因此根据角度递增情况,以正切值占用1~2个字节为原则,将角度对应的正切值按不同的放大倍数放入不同的子表中,相同子表中的角度正切值放大倍数相同,此外,为保证后续查表的快速性,设置同一子表中的正切值具有相同的特征,即高字节内容相同,在查表过程中,只查询低字节数据。
在确定A值所在的子表地址时,采用顺序查询的方式,首先将A值按子表1的放大倍数进行放大,查询其是否具备子表1中内容的特征,若具备则可确定A值所对应的角度值在子表1中,记录子表1在总表中的偏移量地址及子表表首和表尾地址,并继续下一工作步骤;若不具备子表1的特征则将A值按子表2的放大倍数进行放大,按上述过程直到找到A所在子表地址为止。
为进一步缩短查表时间,采取了二分法,首先将A与子表中间位置对应的数据A中值进行比较,若A=A中值,则直接执行下一工作步骤;若A>A中值,则将表首移到表中间位置的后一个地址;若A<A中值,则将表尾移到表中间位置的前一个地址,重复上述过程,直到查到与A值相等的正切值或表长为1时为止。当表长为1时,仍未找到子表中与A相等的数据,则查找表中最相近的数据位置,该数值记作B。
表中间位置的计算用公式四表示为:
Addr中=Addr表首+L/2
式中L表示表长,其计算公式用公式五表示为:
L=Addr表尾-Addr表首+1
最终确定的A值对应的角度为其所在子表在总表中的偏移量地址与对应正切值B在子表中的相对地址之和。
步骤a2、分别求步骤a1中计算出目标指向定位角的正切值A1和A2对应的角度,首先根据目标指向定位角的精度要求,按照表1所示的结构设置表格,其角度间隔为0.1°,共包含11个子表,由于89.8°、89.9°正切值较大,且90°对应的正切值为无穷大,在此进行单独处理,不列在子表中;
表1
表序号 | 角度值 | 表长 | 放大倍数 | 特征 |
子表1 | 0.0~14.0 | 8CH | 1024 | 高字节为00H |
子表2 | 14.1~26.5 | 7CH | 512 | 高字节为00H |
子表3 | 26.6~45.0 | B8H | 512 | 高字节为01H |
子表4 | 45.1~63.4 | B7H | 256 | 高字节为01H |
子表5 | 63.5~75.9 | 7CH | 128 | 高字节为01H |
子表6 | 76.0~82.8 | 28H | 64 | 高字节为01H |
子表7 | 82.9~86.4 | 24H | 32 | 高字节为01H |
子表8 | 86.5~88.2 | 12H | 16 | 高字节为01H |
子表9 | 88.3~89.1 | 09H | 8 | 高字节为01H |
子表10 | 89.2~89.5 | 04H | 4 | 高字节为01H |
子表11 | 89.6~89.7 | 02H | 2 | 高字节为01H |
步骤a3、调整待查数据A1或A2的放大倍数并查询其符合的子表特征,确定其对应的角度所在的子表地址并记录子表序号,若A1或A2不在任何子表中,则判断是否为89.8°、89.9°或90°三个特殊值,是则直接跳过步骤a5和步骤a6,输出角度结果,结合图1所示;
步骤a4、根据步骤a3中记录的子表序号确定该子表表首和表尾地址,并根据公式四和公式五计算表中间值地址和表长,将A1和A2同表中间值进行比较判断,并根据判断结果转入不同的分支:
1)若待查数据值子表中间值相等则直接执行步骤a6;
2)若待查数据值大于表中间值,则将表中间值对应的地址的下一个地址作为新表首地址,重复步骤a4,直到找到相等的数据值或表长为1时为止;
3)若待查数据值小于表中间值,则将表中间值对应的地址的前一个地址作为新表尾地址,重复步骤a4,直到找到相等的数据值或表长为1时为止;
步骤a5、当表长为1时,仍未找到相等的数据值,则将当前表首(即为表尾)所对应的数值同待查数据值进行比较判断:
1)若对应数值大于待查数据值,则求二者差值,并记录该差值,取该表首前一地址对应的数值同待查数据值进行差值,比较两个差值的大小,取差值小的地址所对应的数值地址;
2)若对应数值小于待查数据值,求二者差值,并记录该差值,然后取该表首后一地址对应的数值同待查数据值进行差值,比较两个差值的大小,取差值小的地址所对应的数值地址;
步骤a6、目标指向定位角为:(子表在总表中的偏移量地址+对应正切值在子表中的相对位置)/角度间隔0.1°。
本实施例所述的航空遥感器中负责目标指向定位角计算的微处理器选用Intel公司的MD87C52,编程语言选用汇编语言。
Claims (4)
1.航空遥感器成像时目标指向定位角的计算方法,其特征是,该方法包括以下步骤:
步骤一、采用坐标平移的方法将待拍摄目标在地理坐标系下的坐标转换到航空遥感器所在的载机坐标系下的坐标,所述待拍摄目标在载机坐标系下的坐标计算公式用公式一表示为:
步骤二、根据步骤一所述的待拍摄目标在载机坐标系下的坐标获得目标指向定位角,所述目标指向定位角用公式二和公式三表示为:
公式二:
公式三:
式中:α为垂直飞行方向的旋转角;β为沿飞行方向的旋转角;
所述目标指向定位角通过逆向查表结构和二分查找法获得,具体步骤为:
步骤a:设置总表结构,将待查数据的正切值按角度均匀递增的顺序放入总表中,将所述总表分成多个子表,确定每个子表的表首地址;
步骤b:采用顺序查找法确定待查数据所属的子表在总表中的偏移量地址;然后采用二分法进行查找,获得与待查数据相等或相近角度的正切值在子表中的地址;
步骤c:根据步骤b所述子表在总表中的偏移量地址、与待查数据相同的数据正切值在子表中的地址及总表中设置的角度间隔计算待查数据的角度,所述待查数据的角度为目标指向定位角。
2.根据权利要求1所述的航空遥感器成像时目标指向定位角的计算方法,其特征在于,所述子表中间地址的计算公式用公式四表示:
公式四:Addr中=Addr表首+L/2
式中L表示表长,计算公式用公式五表示:
公式五:L=Addr表尾-Addr表首+1。
3.根据权利要求1所述的航空遥感器成像时目标指向定位角的计算方法,其特征在于,步骤A所述的将待查数据的正切值按角度均匀递增的顺序放入表中,设置多个子表的具体原则为:以待查数据的正切值占用1~2个字节为原则,并将待查数据的正切值按不同的放大倍数放入不同的子表中,相同子表中的数据正切值的放大倍数相同,所述相同子表中的数据正切值的高字节相同。
4.根据权利要求1所述的航空遥感器成像时目标指向定位角的计算方法,其特征在于,步骤c所述的待查数据角度的计算方法为:子表在总表中的偏移量地址与与待查数据相同的数据正切值在子表中的相对位置求和后与总表中设置的角度间隔相除。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110095512A CN102252649B (zh) | 2011-04-15 | 2011-04-15 | 航空遥感器目标指向定位角的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110095512A CN102252649B (zh) | 2011-04-15 | 2011-04-15 | 航空遥感器目标指向定位角的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102252649A true CN102252649A (zh) | 2011-11-23 |
CN102252649B CN102252649B (zh) | 2012-09-19 |
Family
ID=44980066
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110095512A Expired - Fee Related CN102252649B (zh) | 2011-04-15 | 2011-04-15 | 航空遥感器目标指向定位角的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102252649B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106227291A (zh) * | 2016-07-26 | 2016-12-14 | 中国科学院自动化研究所 | 基于分段式查表法的反正切函数的实现方法及实现装置 |
CN112954136A (zh) * | 2021-01-29 | 2021-06-11 | 中国科学院长春光学精密机械与物理研究所 | 抑制航空斜视远距离成像遥感图像散粒噪声的方法及装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101790048A (zh) * | 2010-02-10 | 2010-07-28 | 深圳先进技术研究院 | 智能摄像系统及方法 |
-
2011
- 2011-04-15 CN CN201110095512A patent/CN102252649B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101790048A (zh) * | 2010-02-10 | 2010-07-28 | 深圳先进技术研究院 | 智能摄像系统及方法 |
Non-Patent Citations (2)
Title |
---|
《中国博士学位论文全文数据库》 20020615 金光 机载光电跟踪测量的目标定位误差分析和研究 , * |
《光学精密工程》 20060630 翟林培等 考虑飞机姿态角时倾斜航空相机像移速度计算 第14卷, 第3期 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106227291A (zh) * | 2016-07-26 | 2016-12-14 | 中国科学院自动化研究所 | 基于分段式查表法的反正切函数的实现方法及实现装置 |
CN112954136A (zh) * | 2021-01-29 | 2021-06-11 | 中国科学院长春光学精密机械与物理研究所 | 抑制航空斜视远距离成像遥感图像散粒噪声的方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN102252649B (zh) | 2012-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102411136B (zh) | 一种扩展基线解模糊的相位干涉仪测向方法 | |
CN105259536A (zh) | 一种基于l型天线阵列的相位干涉仪测向装置及解算方法 | |
CN103383450A (zh) | 共形阵列雷达幅相误差校正快速实现方法 | |
CN102135435A (zh) | 一种数字太阳敏感器的误差修正方法及装置 | |
CN110906924A (zh) | 一种定位初始化方法和装置、定位方法和装置及移动装置 | |
CN108663654B (zh) | 一种基于连续量子鸽群的360度全方位动态测向方法 | |
CN117347945B (zh) | 一种基于天线阵三维布局的干涉仪系统测向方法 | |
CN112051569B (zh) | 雷达目标跟踪速度修正方法及装置 | |
CN107451957B (zh) | 一种星载tdi cmos相机成像仿真方法及设备 | |
CN104535067B (zh) | 一种基于分区搜索的脉冲信号到达时间快速计算方法 | |
CN102252649B (zh) | 航空遥感器目标指向定位角的计算方法 | |
CN109781059A (zh) | 星载点波束天线对地指向精度评估系统 | |
CN111486867B (zh) | 一种视觉和惯性混合跟踪组件安装参数的标定装置及方法 | |
CN106986049B (zh) | 一种深空借力轨道精确并行优化设计方法 | |
Zhao et al. | Real-time star identification using synthetic radial pattern and its hardware implementation | |
CN107727945A (zh) | 基于无人机的大型抛物面天线面型精度测试系统 | |
CN117112980A (zh) | 一种盾构掘进过程中盾构位姿实时正解方法 | |
CN103884336A (zh) | 一种基于白昼星敏感器的红外恒星探测导航星表的建立方法 | |
CN108507569A (zh) | 用于星光/惯性复合制导的弹载恒星星库快速生成方法 | |
CN103323025B (zh) | 一种太阳光斑质心及其太阳输出角的误差修正方法 | |
CN112799058B (zh) | 一种基于四臂螺旋天线的二维测角方法 | |
CN111123974B (zh) | 用于大口径射电望远镜天文跟踪模式的轨迹规划方法 | |
CN109308078B (zh) | 定日镜控制方法、装置及计算机可读存储介质和终端设备 | |
CN109724605B (zh) | 基于边缘传感器的大口径天线面板位置姿态实时计算方法 | |
CN104240292B (zh) | 一种交会对接逼近段成像敏感器非目标成像的仿真方法 |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120919 Termination date: 20140415 |