CN105261023B - 一种大气折射修正状态下目标重建方法 - Google Patents
一种大气折射修正状态下目标重建方法 Download PDFInfo
- Publication number
- CN105261023B CN105261023B CN201510683548.5A CN201510683548A CN105261023B CN 105261023 B CN105261023 B CN 105261023B CN 201510683548 A CN201510683548 A CN 201510683548A CN 105261023 B CN105261023 B CN 105261023B
- Authority
- CN
- China
- Prior art keywords
- camera
- target
- angle
- coordinate system
- earth
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30212—Military
Landscapes
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种大气折射修正状态下目标重建方法,属于图像处理的三维重建方法,解决现有目标重建方法未考虑大气折射对光线传播影响而导致的目标定位误差较大的问题,可用于远距离预警系统。本发明方法包括:(1)确定过目标和地心的直线方程;(2)根据相机C1成像信息确定目标位置P1';(3)根据相机C2成像信息确定目标位置P2';(4)得到准确的目标实际位置P'。本发明考虑了大气折射的影响,并针对不同地区、不同季节、不同高度时大气的折射率不同的问题,计算大气折射率时,按照不同地区、不同季节、不同高度的情况考虑,本发明重建得到的目标位置更接近目标实际位置,具有更高的目标定位精度,为高精度预警信息的获取提供了保障,可提升预警系统的整体性能。
Description
技术领域
本发明属于图像处理技术领域,更具体地,涉及一种大气折射修正状态下目标重建方法,可用于卫星导弹的预警,提高预警防御的准确性。
背景技术
高技术条件下的未来战争,突发性强,时间概念高度浓缩,对导弹的预警显得格外重要,需要能够在第一时间得到导弹弹道、方向、速度、发射点与弹着点及导弹数量等重要信息,以便于防御和反导。卫星导弹预警是国家预警系统中关键的一部分,具有极其重要的地位与作用。信息处理是卫星导弹预警系统的核心,主要包括预警卫星红外成像、弱小目标的检测与识别、轨迹重建研究等。导弹的轨迹重建是导弹卫星预警中的关键技术,根据导弹目标在像面成像位置,相机参数等信息重建出导弹的轨迹,以便导弹的拦截预警。现有的轨迹重建技术都是在不考虑大气折射的影响的情况下,按照光直线传播原理,根据目标像面位置,按照直线进行的轨迹重建。但在实际的光线传输过程中,由于大气折射的影响,光线的传输路径会发生偏折,光的传输路径并不是直线。所以现有的目标重建技术得到的目标轨迹与目标实际的轨迹之间存在着误差,不能较为准确地定位目标的位置。
发明内容
针对现有的目标重建技术的缺陷和改进需求,本发明提供一种大气折射修正状态下目标重建技术,与现有的按照光直线传播原理进行的目标重建技术相比,本发明考虑了大气折射影响,定位的目标位置更接近于目标实际位置,满足同步卫星导弹预警的要求,定位精度准确、可靠性高。
本发明所提供的一种大气折射修正状态下目标重建方法,包括确定目标和地心的直线方程,根据相机C1确定目标实际位置P1',根据相机C2确定目标实际位置P2',确定目标实际位置P'步骤,本发明包括如下步骤:
S1、根据两相机位置、由两相机发出的初始光线方向确定过目标和地心的直线方程,包括以下子步骤:
S1.1、已知相机C1的经度Lon1、纬度Lat1、地心高度H1,相机C2的经度Lon2、纬度Lat2、地心高度H2,将两相机C1、C2的经纬度、高度坐标分别转换为地心坐标系下的坐标C1=(X1,Y1,Z1)',C2=(X2,Y2,Z2)';具体转换公式如下:
Xi=Hi·cos(Lati·π/180)·cos(Loni·π/180)
Yi=Hi·cos(Lati·π/180)·sin(Loni·π/180)
Zi=Hi·sin(Lati·π/180)...i=1,2
S1.2、根据目标在相机C1成像的像面位置p1=(u1,v1)',相机光轴指向相机在地心坐标系下的坐标C1=(X1,Y1,Z1)'及相机内部参数,求得由相机C1发出并最终到达目标的光线在地心坐标系下的初始方向向量具体包括以下步骤:
S1.2.1、对相机C1,根据C1在地心坐标系下的坐标C1=(X1,Y1,Z1)',相机光轴的指向求得相机坐标系到地心坐标系的转换矩阵MC1。
S1.2.2、根据目标在相机C1成像的像面位置p1=(u1,v1)',相机C1内部参数:相机所成图像的行数row、列数col,相机的像源瞬时视场角λ。计算由相机C1发出并最终到达目标的光线在相机坐标系下的初始方向向量
S1.2.3、由相机坐标系到地心坐标系的转换矩阵MC1,由相机C1发出并最终到达目标的光线在相机坐标系下的初始方向向量求得由相机C1发出并最终到达目标的光线在地心坐标系下的初始方向向量即
S1.3、同步骤S1.2的方法,根据目标在相机C2成像的像面位置p2=(u2,v2)',相机光轴指向相机在地心坐标系下的坐标C2=(X2,Y2,Z2)'及相机内部参数,求得由相机C2发出并最终到达目标的光线初始方向在地心坐标系下的向量
S1.4、因为大气在水平方向视为均匀分布,大气折射只发生在垂直方向上,所以,由初始方向向量地心到相机C1的向量所组成的平面,即为相机C1和目标实际位置P'所在的过地心的平面。根据方向向量相机C1坐标C1=(X1,Y1,Z1)'、地心O坐标O=(0,0,0)',求得相机C1、目标实际位置P'以及地心所组成的平面OC1P'的法向量
S2、根据相机C1、过目标和地心的直线方程,确定目标实际位置P1‘;,包括以下子步骤:
S2.1、对于相机C1所在的平面OC1P',由直线OP'方程、相机C1以及地心O,求得目标实际位置P'与相机C1间的地心夹角θ1,将地心夹角θ1分为num1(100≤num1≤300)个相等的小夹角,对应的目标和相机间的大气被划分为num1层,每层地心夹角大小为dθ,dθ=θ1/num1。因为每层地心夹角很小,这使得每层大气的厚度很薄,所以可以取每层入射点位置对应的大气折射率作为该层大气的折射率。
平面坐标系到地心坐标系的转换矩阵Mplan-earth为:
S2.4、根据地心坐标系到平面坐标系的转换矩阵Mearth-plan,将地心坐标系下光线初始方向向量转换为平面坐标系中的向量
S2.5、在平面OC1P'中,以相机C1所在的位置为光线由第0层大气射向第1层大气时,在第1层的入射点位置。根据光线的初始方向向量相机坐标Cplan1,计算光线由第0层大气射向第1层大气时,第1层出射光线的方向向量,以及第2层入射光线的入射点位置。具体包括以下步骤:
S2.5.1、相机C1所在位置为光线由第0层射向第1层时的入射点位置T1;光线初始方向向量即为光线由第0层大气射向第1层大气时,第1层的入射光线方向向量根据相机C1的经纬度坐标将相机C1对应的地面位置的海拔高度、地面温度、地面压强、地面相对湿度、季节、以及C1处地心高度H1参数,代入大气折射率Hopfiled模型,计算C1处的大气折射率N1,根据大气折射指数公式n=N×10-6+1,计算第1层大气折射指数n1。考虑实际情况,第0层大气折射指数n0和第1层大气折射指数n1取相同值,n0=n1。光线由第0层大气射向第1层大气时不发生折射,第1层的出射光线方向向量等于第1层的入射光线方向向量所以求得第1层出射光线方向的单位向量
S2.5.3、在由地心O、第1层的入射点位置T1、第2层的入射点位置T2组成的三角形OT1T2中,由正弦定理,计算光线在第1层中的路径长度L1,具体公式为:
S2.6、当光由第i(i=1,2,3,...,num1-1)层大气射向第i+1层大气时,计算第i+1层出射光线方向向量以及第i+2层光线入射点位置,具体包括:
S2.6.2、根据第i+1层入射点位置的经纬度坐标,将入射点位置Ti+1对应地面处的海拔高度、地面温度、地面压强。地面相对湿度、入射点位置Ti+1的地心高度Hi+1、季节参数代入到Hopfield模型中,计算第i+1层大气折射率Ni+1,由大气折射指数公式n=N×10-6+1,计算第i+1层大气折射指数ni+1。
S2.6.3、判断第i+1层大气折射指数ni+1与第i层大气折射指数ni之间的大小关系,如果ni+1小于ni,此时可能会发生全反射,则执行步骤S2.6.4;否则直接执行步骤S2.6.5。
S2.6.4、根据全反射临界角公式计算临界角θ。判断入射角αi+1与临界角θ大小关系,如果入射角αi+1小于临界角θ,则不会发生全反射,执行步骤S2.6.5,并且每次执行步骤S2.6.5后都要重新计算临界角并判断与入射角的大小关系;如果入射角αi+1大于等于临界角θ,则光线在该层发生全反射,光线的传输方向发生变化,反射角εi+1等于入射角αi+1。此时可以将第i+1层出射角看作反射角εi+1的补角,求得第i+1层出射光线的出射角βi+1,βi+1=π-εi+1,继续执行步骤S2.6.6,此后每一层折射时不需要再判断入射角与临界角的大小关系。
大气球面分层的Snell定律为:当光线由大气折射指数为n1的大气层,射向大气折射指数为n2的大气层时,满足下面的关系:
n1sinφ1=n2sinφ2
式中φ1、φ2分别为入射角和出射角,它们分别定义为入射方向和出射方向与界面法线的夹角,在地球大气层中,界面法线为地心与入射点的连线。
S2.6.6、由第i+1层入射点位置处的法线向量出射光线的出射角βi+1,利用向量内积,计算第i+1层出射光线的单位向量具体求解方法为:
S2.7、当总层数等于num1-1,此时对应的下一层光线入射点位置即为通过相机C1求得的目标位置P1'。通过平面坐标系到地心坐标系的转换矩阵Mplan-earth将平面坐标系下的坐标转换到地心坐标系下的坐标根据地心坐标到经纬度、高度的转换公式,求得处的经纬度、高度坐标并记录,此坐标即为通过相机C1求得的目标位置P1'。
S3、根据相机C2、过目标和地心的直线方程,确定目标实际位置P2‘;
对于相机C2所在的平面OC2P'中,求目标实际位置P2'方法与步骤S2相同。最终根据相机C2求得目标实际位置为P2'。
S4、根据两相机C1、C2确定的目标实际位置P1'、P2',最终得到准确的目标实际位置P'。
为了使目标定位更准确,最终目标实际位置P'取两目标位置P1'、P2'的平均值,即P'=(P1'+P2')/2。
本发明根据目标在相机成像的像面位置,通过三维目标重建的方法获得目标的实际位置。相比于现有的目标重建方法,本发明考虑了光线传输过程中大气折射的影响,并针对不同地区、不同季节、不同高度时大气的折射率不同的问题,在计算大气折射率时,按照不同地区、不同季节、不同高度的情况考虑;并在考虑大气折射传播时,考虑了发生全反射时的光线传输路径。本发明重建的目标位置更接近于目标真实位置,提高了目标定位的精确度,可更有效地应用于导弹预警,提高预警防御的准确性。
附图说明
图1本发明大气折射修正状态下目标重建方法的流程示意图;
图2大气折射修正状态下目标重建光线传输图;
图3地心坐标系与平面坐标系之间的关系图;
图4相机、目标、地心所在的平面内,光线由第0层射向第1层传输路径图;
图5相机、目标、地心所在的平面内,光在每层大气中的传输路径图;
图6发生全反射时,光线传输路径图;
图7相机所在无人机与目标间的直线距离、两无人机间的夹角、无人机的姿态、目标距离地面的高度、相机光轴指向一定时,不同的像面误差与重建误差间的关系图;
图8相机所在无人机与目标间的直线距离、两无人机间的夹角、无人机的姿态、目标距离地面的高度、目标在像面成像位置一定时,不同的相机光轴指向误差与重建误差间的关系图;
图9相机所在无人机与目标间的直线距离、两无人机间的夹角、目标距离地面的高度、相机光轴指向、目标在像面成像位置一定时,不同的无人机姿态误差与重建误差间的关系图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
由目标发出的光线,经过大气传输,被相机接收并在像面上成像。根据光路可逆原则,由目标发出并最终到达相机的光线传输路径与由相机发出并最终到达目标的光线传输路径是相同的。所以在目标重建时,我们将光线传输路径视为由相机出发并最终到达目标位置。根据目标的像面点位置,可以得到最终入射光线方向,该方向的反方向即为从相机出发并最终到达目标的光线初始方向,大气折射修正状态下目标重建的光线传输图,如图2所示。本发明流程如图1所示,本发明提供一种折射修正状态下目标重建方法,包括以下步骤:
S1、根据两相机位置、由两相机发出的初始光线方向确定过目标和地心的直线方程;
S2、根据相机C1、过目标和地心的直线方程,确定目标实际位置P1‘;
S3、根据相机C2、过目标和地心的直线方程,确定目标实际位置P2‘;
S4、根据两相机C1、C2确定的目标实际位置P1‘、P2‘,最终得到准确的目标实际位置P'。
其中,上述步骤S1具体包括以下子步骤:
S1.1、已知相机C1的经度Lon1、纬度Lat1、地心高度H1,相机C2的经度Lon2、纬度Lat2、地心高度H2,将两相机C1、C2的经纬度、高度坐标分别转换为地心坐标系下的坐标C1=(X1,Y1,Z1)',C2=(X2,Y2,Z2)'。具体转换公式如下:
Xi=Hi·cos(Lati·π/180)·cos(Loni·π/180)
Yi=Hi·cos(Lati·π/180)·sin(Loni·π/180)
Zi=Hi·sin(Lati·π/180)...i=1,2
S1.2、根据目标在相机C1成像的像面位置p1=(u1,v1)',相机光轴指向相机在地心坐标系下的坐标C1=(X1,Y1,Z1)'及相机内部参数,求得由相机C1发出并最终到达目标的光线在地心坐标系下的初始方向向量具体包括以下步骤:
S1.2.1、对相机C1,根据C1在地心坐标系下的坐标C1=(X1,Y1,Z1)',相机光轴的指向计算相机坐标系到地心坐标系的转换矩阵MC1。
下面先介绍相机坐标系的定义,再介绍相机坐标系到地心坐标系的转换矩阵的求取方法。
(1)相机坐标系的定义
定义以相机所在位置为相机坐标系的原点OC1,相机光轴的指向为相机坐标系的ZC1轴,XC1轴垂直于相机和地心的连线OC1与相机光轴指向所组成的平面,YC1轴由右手螺旋定则确定。
(2)相机坐标系到地心坐标系的转换矩阵
S1.2.2、根据目标在相机C1成像的像面位置p1=(u1,v1)',相机C1内部参数:相机所成图像的行数row、列数col,相机的像源瞬时视场角λ;计算由相机C1发出并最终到达目标的光线在相机坐标系下的初始方向向量具体计算公式如下:
γ1=(u1-row/2)·λ
γ2=(v1-col/2)·λ
S1.2.3、由相机坐标系到地心坐标系的转换矩阵MC1,由相机C1发出并最终到达目标的光线在相机坐标系下的初始方向向量计算由相机C1发出并最终到达目标的光线在地心坐标系下的初始方向向量即
S1.3、同步骤S1.2的方法,根据目标在相机C2成像的像面位置p2=(u2,v2)',相机光轴指向相机在地心坐标系下的坐标C2=(X2,Y2,Z2)'及相机内部参数,计算由相机C2发出并最终到达目标的光线在地心坐标系下的初始方向向量
S1.4、因为大气在水平方向视为均匀分布,大气折射只发生在垂直方向上,所以,由初始方向向量地心到相机C1的向量所组成的平面,即为相机C1和目标实际位置P'所在的过地心的平面。根据方向向量相机C1坐标C1=(X1,Y1,Z1)'、地心O坐标O=(0,0,0)',计算相机C1、目标实际位置P'所在过地心的平面OC1P'的法向量求平面OC1P'的法向量具体方法如下:
在平面OC1P'中,由地心到相机的向量和光线初始方向向量的叉积可以得到平面OC1P'的法向量
S1.6、由成像原理可知,目标实际位置P'既在平面OC1P'上,又在平面OC2P'上,因此目标位于两平面的交线OP'上。根据两平面求得交线OP'的直线方程,该方程即为过目标和地心的直线方程,交线OP'的方向向量可由两平面的法向量的叉积求得,即
步骤S2具体包括以下子步骤:
S2.1、对于相机C1所在的平面OC1P',由直线OP'方程、相机C1以及地心O,计算目标实际位置P'与相机C1间的地心夹角θ1,将地心夹角θ1分为num1(100≤num1≤300)个相等的小夹角,对应的目标和相机间的大气被划分为num1层,每层地心夹角大小为dθ,dθ=θ1/num1。因为每层地心夹角很小,这使得每层大气的厚度很薄,所以可以取每层入射点位置对应的大气折射率作为该层大气的折射率。
S2.2、求地心坐标系与平面坐标系之间的转换矩阵。
下面先介绍平面坐标系的定义,以及平面坐标系与地心坐标系之间的转换矩阵:
(1)相机、目标以及地心所组成的平面坐标系定义如下:
如图3所示,设相机C在地心坐标系中坐标(xc,yc,zc)',目标P在地心坐标系中坐标(xp,yp,zp)',地心O坐标(0,0,0)'。定义平面坐标系O-XpYpZp的原点与地心直角坐标系的原点重合,设平面OCP与平面OXY的交线OXp为平面坐标系的Xp轴;Yp轴垂直平面OCP;则由右手螺旋定则确定Zp轴。
(2)地心坐标系与平面坐标系之间的转换矩阵
从而地心坐标系O-XYZ到平面坐标系O-XpYpZp的转换矩阵为:
平面坐标系O-XpYpZp到地心坐标系O-XYZ的转换矩阵为:
利用地心坐标系与平面坐标系的转换矩阵,可将任一点地心坐标系下的坐标(X,Y,Z)'与平面坐标系下的坐标(Xp,Yp,Zp)'相互转换,其转化关系如下:
对于平面OCP中的任一点,在平面坐标系中,其Yp轴坐标始终为0,即平面中一点在地心坐标系下坐标(xe,ye,ze)',经转换矩阵得到平面坐标系下坐标(xp,0,zp)'。且已知光线只在垂直方向上发生折射,在水平方向不发生偏折,所以由相机发出到达目标的光线传输路径上每一点,在平面坐标系中Yp轴坐标始终为0。
平面坐标系到地心坐标系的转换矩阵Mplan-earth为:
S2.4、根据地心坐标系到平面坐标系的转换矩阵Mearth-plan,将光线在地心坐标系下初始方向向量转换为平面坐标系中的向量
S2.5、在平面OC1P'中,以相机C1所在的位置为光线由第0层大气射向第1层大气时,在第1层大气的入射点位置。根据光线的初始方向向量相机坐标Cplan1,计算光线由第0层大气射向第1层大气时,第1层出射光线的方向向量,以及第2层入射光线的入射点位置。具体包括以下步骤:
S2.5.1、相机C1所在位置为光线由第0层大气射向第1层大气时的入射点位置T1;光线初始方向向量即为光线由第0层大气进入第1层大气时,第1层入射光线的方向向量根据相机C1的经纬度坐标将C1对应地面位置的海拔高度、地面温度、地面压强、地面相对湿度、季节参数、C1处地心高度H1,代入大气折射率Hopfiled模型,计算C1处的大气折射率N1,根据大气折射指数公式n=N×10-6+1,计算第1层大气折射指数n1。考虑实际情况,第0层大气折射指数n0和第1层大气折射指数n1取相同值,n0=n1。此时光由第0层射向第1层时不发生折射,如图4所示,第1层出射光线的方向向量等于第1层入射光线的方向向量所以求得第1层出射光线的单位向量
S2.5.3、如图4所示,在三角形OT1T2中,由正弦定理,计算光线在第1层中的路径长度L1,具体公式为:
大气折射率的Hopfield模型具体原理如下:
利用大气折射率的Hopfield模型,可以求取任意海拔高度h的大气折射率N(h),Hopfield模型考虑了大气含有水汽的情况,将大气折射率分为干、湿两项,模型的具体形式为:
S2.6、当光由第i(i=1,2,3,...,num1-1)层大气射向第i+1层大气时,计算第i+1层出射光线的方向向量以及第i+2层入射光线的入射点位置,具体包括如下子步骤:
S2.6.2、根据第i+1层入射点位置的经纬度、高度坐标,将入射点Ti+1对应地面处的海拔高度、地面温度、地面压强。地面相对湿度、入射点Ti+1处的地心高度Hi+1、季节参数代入到Hopfield模型中,计算第i+1层大气折射率Ni+1,由大气折射指数公式n=N×10-6+1,计算第i+1层大气折射指数ni+1。
S2.6.3、判断第i+1层大气折射指数ni+1与第i层大气折射指数ni之间的大小关系,如果ni+1小于ni,此时可能会发生全反射,则执行步骤S2.6.4;否则直接执行步骤S2.6.5。
S2.6.4、根据全反射临界角公式求出临界角θ。判断入射角αi+1与临界角θ大小关系,如果入射角αi+1小于临界角θ,则不会发生全反射,执行步骤S2.6.5,并且每次执行步骤S2.6.5后都要重新计算临界角并判断与入射角的大小关系;如果入射角αi+1大于等于临界角θ,则光线在该层发生全反射,光线的传输方向发生变化,反射角εi+1等于入射角αi+1,如图6所示,此时可以将第i+1层出射角看作反射角的补角,求得第i+1层出射光线的出射角βi+1,βi+1=π-εi+1,继续执行步骤S2.6.6,此后每一层折射时不需要再判断入射角与临界角的大小关系。
大气球面分层的Snell定律为:当光线由大气折射指数为n1的大气层,射向大气折射指数为n2的大气层时,满足下面的关系:
n1sinφ1=n2sinφ2
式中φ1、φ2分别为入射角和出射角,它们分别定义为入射方向和出射方向与界面法线的夹角,在地球大气层中,界面法线为地心与入射点的连线。
S2.6.6、由第i+1层入射点位置Ti+1处的法线向量出射角βi+1,利用向量内积,计算第i+1层出射光线的单位向量具体求解方法为:
S2.7、当总层数等于num1-1时,此时对应的下一层光线入射点位置即为通过相机C1求得的目标位置P1'。通过平面坐标系到地心坐标系的转换矩阵Mplan-earth将在平面坐标系下的坐标转换到地心坐标系下的坐标根据地心坐标到经纬度、高度的转换公式,求得处的经纬度、高度坐标并记录,此坐标即为通过相机C1求得的目标位置P1'。
步骤S3根据相机C2成像信息确定目标位置P2',求目标位置P2'的方法与步骤S2中根据相机C1求目标位置P1'的方法相同。最终根据相机C2求得的目标位置P2'。
步骤S4具体包括以下子步骤:
通过相机C1、C2分别确定了目标的位置P1'、P2',由于计算时会存在误差,为了使目标定位更精确,最终目标实际位置P'取两平面求得的目标位置P1'、P2'的平均值,即P'=(P1'+P2')/2。
实验结果分析:
图7~9表示分别对目标在像面位置、相机光轴指向、相机所在的无人机平台的姿态添加不同的误差时,利用本发明的方法进行目标位置重建,得到重建的目标位置与目标实际位置间的误差,称该误差为重建误差。以判断本专利方法在重建时有外界干扰情况下,对目标定位的精确度。
图7为对目标的像面成像位置分别添加3σ=0,0.5,1,1.5个像素的像面位置误差时,得到的重建误差结果;图8为对相机的光轴指向添加3σ从0.001°到0.01°间隔为0.001°的指向误差,得到的重建误差结果;图9为对相机所在的无人机平台的姿态添加3σ从0.005°到0.02°间隔0.005°的姿态误差,得到的重建误差结果。
从上面的实验结果可以看出,本发明不仅重建的目标位置精确度高,而且抗外界环境干扰能力强,在由外界误差干扰的情况下,目标重建误差在1km以内,这对提高导弹的预警防御具有重要作用。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种大气折射修正状态下的目标重建方法,其特征在于,所述方法包括步骤:
S1、根据两相机位置、由两相机发出的初始光线方向确定过目标和地心的直线方程;
S2、根据相机C1、过目标和地心的直线方程,确定目标实际位置P1';
S3、根据相机C2、过目标和地心的直线方程,确定目标实际位置P2';
S4、根据两相机C1、C2确定的目标实际位置P1'、P2',最终得到准确的目标实际位置P';
所述步骤S2包括以下子步骤:
S2.1、由直线OP'方程、相机C1以及地心O,计算目标实际位置P'与相机C1间的地心夹角θ1,并将地心夹角分为num1个相等的小夹角,100≤num1≤300;
S2.2、计算地心坐标系与平面坐标系之间的转换矩阵Mearth-plan;
S2.3、根据地心坐标系到平面坐标系的转换矩阵Mearth-plan,将相机C1在地心坐标系下的坐标C1=(X1,Y1,Z1)'转换为平面坐标系下的坐标,并记为二维平面坐标的形式,即Cplan1=(Xplan(1),Zplan(1))';
S2.4、根据地心坐标系到平面坐标系的转换矩阵Mearth-plan,将地心坐标系下光线的初始方向向量转换为平面坐标系下的初始方向向量
S2.5、在平面OC1P'中,以相机C1所在的位置为光线由第0层大气射向第1层大气时,在第1层大气的入射点位置;根据平面坐标系下光线的初始方向向量相机坐标Cplan1,计算光线由第0层大气射向第1层大气时,第1层大气的出射光线的方向向量,以及光线在第2层大气的入射点位置;
S2.6、当光由第i层大气射向第i+1层大气时,i=1,2,3,...,num1-1,计算第i+1层大气的出射光线方向向量以及光线在第i+2层大气的入射点位置,其中num1为步骤S2.1中将地心夹角等分的小夹角个数;
S2.7、当总层数等于num1-1时,此时对应的下一层的光线入射点位置即为通过相机C1求得的目标位置P1';通过平面坐标系到地心坐标系的转换矩阵Mplan-earth将平面坐标系下的坐标转换到地心坐标系下的坐标根据地心坐标到经纬度、高度的转换公式,求得处的经纬度、高度坐标并记录,此坐标即为通过相机C1求得的目标位置P1'。
2.如权利要求1所述方法,其特征在于,所述步骤S1包括以下子步骤:
S1.1、已知相机C1的经度Lon1、纬度Lat1、地心高度H1,相机C2的经度Lon2、纬度Lat2、地心高度H2,将两相机C1、C2的经纬度、高度坐标分别转换为地心坐标系下的坐标C1=(X1,Y1,Z1)',C2=(X2,Y2,Z2)';
S1.2、根据目标在相机C1成像的像面位置p1=(u1,v1)',相机光轴指向相机在地心坐标系下的坐标C1=(X1,Y1,Z1)'及相机内部参数,计算由相机C1发出并最终到达目标的光线在地心坐标系下的初始方向向量
S1.3、同步骤S1.2的方法,根据目标在相机C2成像的像面位置p2=(u2,v2)',相机光轴指向相机在地心坐标系下的坐标C2=(X2,Y2,Z2)'及相机内部参数,计算由相机C2发出并最终到达目标的光线在地心坐标系下的初始方向向量
S1.4、根据初始方向向量相机C1坐标C1=(X1,Y1,Z1)'、地心O坐标O=(0,0,0)',计算相机C1、目标实际位置P'以及地心所组成的平面OC1P'的法向量:
S1.5、同步骤S1.4的方法,根据方向向量相机C2坐标C2=(X2,Y2,Z2)'、地心O的坐标O=(0,0,0)',计算相机C2、目标实际位置P'以及地心所组成的平面OC2P'的法向量:
S1.6、根据OC1P'和OC2P'两平面法向量,计算交线OP'的直线方程,该方程即为过目标和地心的直线方程。
3.如权利要求2所述方法,其特征在于,所述步骤S1.2具体包括以下步骤:
S1.2.1、根据C1在地心坐标系下的坐标C1=(X1,Y1,Z1)',相机光轴的指向计算相机坐标系到地心坐标系的转换矩阵MC1;
S1.2.2、根据目标在相机C1成像的像面位置p1=(u1,v1)',相机C1内部参数:相机所成图像的行数row、列数col,相机的像源瞬时视场角λ;计算由相机C1发出并最终到达目标的光线在相机坐标系下的初始方向向量
S1.2.3、根据相机坐标系到地心坐标系的转换矩阵MC1,由相机C1发出并最终到达目标的光线在相机坐标系下的初始方向向量计算由相机C1发出并最终到达目标的光线在地心坐标系下的初始方向向量
4.如权利要求2所述方法,其特征在于,所述步骤S1.6中的交线OP'的方向向量由OC1P'和OC2P'两平面的法向量的叉积求得,即
5.如权利要求1所述方法,其特征在于,所述步骤S2.6中具体包括如下子步骤:
S2.6.1、当光由第i层射向第i+1层时,i=1,2,3,...,num1-1,第i+1层的入射点位置Ti+1=(Xplan(i+1),Zplan(i+1))'、第i层的大气折射指数ni以及每层地心夹角dθ已知,由第i+1层入射点位置Ti+1处的法线向量与入射光线的方向向量根据公式求得第i+1层入射光线的入射角αi+1;利用平面坐标系到地心坐标系的转换矩阵Mplan-earth,将第i+1层的入射点位置Ti+1在平面坐标系下的坐标转换到地心坐标系下的坐标Tei+1=(Xi+1,Yi+1,Zi+1)',转换公式为根据地心坐标到经纬度、高度的转换公式,求得第i+1层入射点的经度Loni+1、纬度Lati+1、高度Hi+1,并记录第i+1层的入射点位置处的经纬度、高度坐标Tli+1=(Loni+1,Lati+1,Hi+1)';
S2.6.2、根据第i+1层入射点位置的经纬度坐标,将入射点位置Ti+1对应地面处的海拔高度、地面温度、地面压强、地面相对湿度、入射点的地心高度Hi+1、季节参数代入到霍普菲尔德Hopfield大气模型中,计算第i+1层大气折射率Ni+1,由大气折射指数公式ni+1=Ni+1×10-6,计算第i+1层大气折射指数ni+1;
S2.6.3、判断第i+1层大气折射指数ni+1与第i层大气折射指数ni之间的大小关系,如果ni+1小于ni,此时发生全反射,则执行步骤S2.6.4;否则直接执行步骤S2.6.5;
S2.6.4、根据全反射临界角公式求出临界角θ;判断入射角αi+1与临界角θ大小关系,如果入射角αi+1小于临界角θ,则不会发生全反射,执行步骤S2.6.5,并且每次执行步骤S2.6.5后都要重新计算临界角并判断与入射角的大小关系;如果入射角αi+1大于等于临界角θ,则光线在该层发生全反射,光线的传输方向发生变化,反射角εi+1等于入射角αi+1,此时将第i+1层出射角看作反射角的补角,求得第i+1层出射光线的出射角βi+1,βi+1=π-εi+1,继续执行步骤S2.6.6,此后每一层折射时不需要再判断入射角与临界角的大小关系;
S2.6.5、根据大气球面分层的斯涅尔Snell定律,计算第i+1层出射光线的出射角βi+1,
S2.6.6、由第i+1层入射点Ti+1处的法线向量出射角βi+1,利用向量内积,求得第i+1层出射光线的单位向量
S2.6.7、在由地心O、第i+1层的入射点位置Ti+1、第i+2层的入射点位置Ti+2组成的三角形OTi+1Ti+2中,根据正弦定理,由公式求得光线在第i+1层中的路径长度Li+1,进而求得第i+1层出射光线的方向向量该方向向量即为第i+2层入射光线的方向向量即根据向量间的关系求得第i+2层光线入射点位置Ti+2的坐标,
6.如权利要求5所述方法,其特征在于,所述步骤S2.6.6中的具体求解方法为:
7.如权利要求1或2所述方法,其特征在于,所述步骤S4具体为:通过相机C1、C2分别确定了目标的位置,最终目标实际位置P'取两平面求得的目标位置P1'、P2'的平均值,即P'=(P1'+P2')/2。
8.如权利要求5所述方法,其特征在于,所述步骤S2.6.5中的大气球面分层的Snell定律为:当光线由折射指数为n1的大气层,射向折射指数为n2的大气层时,满足下面的关系:
n1sinφ1=n2sinφ2
式中φ1、φ2分别为入射角和出射角,它们分别定义为入射方向和出射方向与界面法线的夹角,在地球大气层中,界面法线为地心与入射点的连线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510683548.5A CN105261023B (zh) | 2015-10-20 | 2015-10-20 | 一种大气折射修正状态下目标重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510683548.5A CN105261023B (zh) | 2015-10-20 | 2015-10-20 | 一种大气折射修正状态下目标重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105261023A CN105261023A (zh) | 2016-01-20 |
CN105261023B true CN105261023B (zh) | 2018-09-21 |
Family
ID=55100695
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510683548.5A Expired - Fee Related CN105261023B (zh) | 2015-10-20 | 2015-10-20 | 一种大气折射修正状态下目标重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105261023B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108225276B (zh) * | 2018-01-11 | 2020-06-23 | 中国人民解放军战略支援部队航天工程大学 | 一种单星成像目标运动特性反演方法及系统 |
CN108535715A (zh) * | 2018-04-12 | 2018-09-14 | 西安应用光学研究所 | 一种适用于机载光电观瞄系统的大气折射下目标定位方法 |
CN109900658B (zh) * | 2019-03-06 | 2019-11-01 | 中国人民解放军国防科技大学 | 一种利用恒星观测数据校正光学成像卫星大气折射的方法 |
CN112098953A (zh) * | 2020-09-21 | 2020-12-18 | 中国人民解放军63921部队 | 一种大气折射修正量计算快速迭代方法及装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102229362A (zh) * | 2011-04-29 | 2011-11-02 | 航天东方红卫星有限公司 | 敏捷卫星成像反向推扫速度失配时的姿态补偿方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10103917A (ja) * | 1996-09-27 | 1998-04-24 | Nikon Corp | 位置計測装置 |
-
2015
- 2015-10-20 CN CN201510683548.5A patent/CN105261023B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102229362A (zh) * | 2011-04-29 | 2011-11-02 | 航天东方红卫星有限公司 | 敏捷卫星成像反向推扫速度失配时的姿态补偿方法 |
Non-Patent Citations (2)
Title |
---|
A Refraction’s Slacking in Optoelectronic Systems for Positioning of Bulky Objects Elements;Alexander V.Prokofjev等;《3rd International Symposium on Instrumentation Science and Technology》;20040801;第1-4页 * |
一种空中飞行目标空中姿态参数的算法;孙少华等;《地矿测绘》;20080325;第24卷(第1期);第6-8页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105261023A (zh) | 2016-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106403902B (zh) | 一种星地协同的光学卫星在轨实时几何定位方法及系统 | |
CN104482934B (zh) | 一种多传感器融合的超近距离自主导航装置与方法 | |
CN103323855B (zh) | 一种基线动态测量系统的精度获取方法 | |
CN105261023B (zh) | 一种大气折射修正状态下目标重建方法 | |
CN107782293B (zh) | 基于六自由度激光跟踪靶的航天器设备位姿信息测量方法 | |
CN109708649B (zh) | 一种遥感卫星的姿态确定方法及系统 | |
ES2834997T3 (es) | Un procedimiento y sistema para estimar información relacionada con el ángulo de cabeceo y/o alabeo de un vehículo | |
CN107490364A (zh) | 一种大角度倾斜成像航空相机对地目标定位方法 | |
CN111537003A (zh) | 一种基于折射面共线的星光大气折射测量校正方法 | |
CN109916406B (zh) | 一种基于无人机群的环绕目标定位方法 | |
CN104880200B (zh) | 复合制导系统初始姿态现场校准系统及方法 | |
CN110146093A (zh) | 双体小行星探测自主协同光学导航方法 | |
CN108225276B (zh) | 一种单星成像目标运动特性反演方法及系统 | |
CN107063187A (zh) | 一种全站仪与无人机影像联合的树高快速提取方法 | |
CN106489062A (zh) | 用于测量移动平台的位移的系统和方法 | |
CN105136128A (zh) | 基于两点定位的机体结构测量方法 | |
CN109959898A (zh) | 一种座底式水声被动定位基阵自校准方法 | |
Wang et al. | Atmospheric refraction calibration of geometric positioning for optical remote sensing satellite | |
CN108594255A (zh) | 一种激光测距辅助光学影像联合平差方法及系统 | |
Zhang et al. | Global Measurement Method for Large‐Scale Components Based on a Multiple Field of View Combination | |
KR20230040921A (ko) | 티베트고원(靑藏高原) 지역에서 광학적 원격탐지 위성영상의 대기굴절 위치고정 오차 수정방법 | |
CN106323271B (zh) | 基于特征奇异值的航天器相对姿态测量矢量选取方法 | |
CN105547286B (zh) | 一种复合三视场星敏感器星图模拟方法 | |
CN105424060B (zh) | 一种飞行器星敏感器与捷联惯组安装误差的测量方法 | |
CN103852079B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180921 Termination date: 20201020 |
|
CF01 | Termination of patent right due to non-payment of annual fee |