卫星影像定位方法和装置制造方法及图纸

本发明专利技术实施例提供一种应用于电荷藕合器件图像传感器的卫星影像定位方法和装置，采用本发明专利技术实施例提供的方法和装置，可以实现在没有地面控制点条件下，仅依据外定标获得的定轨定资数据以及卫星的相机在该定轨定资数据下拍摄的卫星影像的连接点的坐标值数据，对卫星的定轨定资数据中的空间姿态数据进行校正。具体的，利用定轨定姿数据和连接点坐标、与连接点对应的地面坐标值的对应关系，通过二次多项式构建空间姿态数据随时间变化的外方位元素方程，基于共线方程，获取连接点的地面坐标值以及利用最小二乘法迭代求解空间姿态数据。

【技术实现步骤摘要】

本专利技术涉及卫星应用
，具体的涉及一种卫星影像定位方法和装置。
技术介绍
目前，卫星遥感影像的精确定位一直依赖于地面控制点，然而获取足够数量的地面控制点通常非常困难，因此研究无地面控制点条件下的卫星遥感影像的高精度几何处理就成为一项紧迫的任务。
技术实现思路
有鉴于此，本专利技术提供了一种应用于电荷藕合器件图像传感器的卫星影像定位方法及装置，以解决现有技术中卫星遥感影像的精确定位一直依赖于地面控制点，然而获取足够数量的地面控制点通常非常困难的问题，其技术方案如下：一种应用于电荷藕合器件图像传感器的卫星影像定位方法，包括：A1、依据外定标获得卫星上的相机的每一扫描行的定轨定资数据，所述定轨定资数据包括所述相机的摄影中心的空间坐标值和所述相机的摄影中心相对于像片面的空间姿态数据；A2、获取所述相机拍摄同一被摄物体的至少两幅卫星影像，将所述至少两幅卫星影像中分别与所述被摄物体的同一点对应的像点的集合称为连接点；A3、获取所述至少两幅卫星影像的影像条带，以所述相机的扫描方向为横坐标方向，以所述相机的飞行方向作为纵坐标方向，构建以所述连接点的横坐标为自变量和以所述连接点的纵坐标为因变量的瞬时构像方程，所述瞬时构像方程包括：当前的定轨定资数据、与所述连接点对应的所述被摄物体上的点的地面坐标值与所述连接点坐标值的关系；A4、对于每一影像条带，根据所述影像条带的起始扫描时间、所述影像条带的终止扫描时间以及所述影像条带的像素个数，确定所述影像条带中每一像素的摄影时刻；A5、以所述当前的空间姿态数据为因变量，所述摄影时刻为自变量，构建描述所述当前空间姿态数据与所述摄影时刻之间的关系的外方位元素方程，所述外方位元素方程的系数为矫正后的系数，所述外方位元素方程的系数的初始值由步骤A1中的定轨定资数据确定；A6、依据所述外方位元素方程，将所述瞬时构像方程整理成以所述地面坐标值为变量的共线方程；A7、对于每一连接点，将所述连接点中各个像点的坐标值带入所述共线方程，计算出所述被摄物体中与所述连接点对应的点的地面坐标值；A8、选择所述至少两幅卫星影像中的第一幅卫星影像；A9、将所述第一幅卫星影像中各个像点分别对应的地面坐标值作为确定值，对于所述第一幅卫星影像中每一像点，将所述第一幅卫星影像中各个像点的坐标值的误差值以及外方位元素方程中各个系数的误差量带入所述共线方程，获得以所述第一幅卫星影像中各个像点的坐标值的误差值作为因变量，以所述外方位元素方程中各个系数的误差量为自变量的误差方程；A10、依据最小二乘间接平差方法，获取所述误差方程的法方程式，所述法方程式包括由所述第一幅卫星影像中各个像点的坐标值的观测权值组成的观测权值矩阵，所述观测权值矩阵的初始值为单位矩阵，所述观测权值矩阵与由前一次所述第一幅卫星影像中各个像点的坐标值的误差值组成的误差值矩阵有关；A11、通过所述法方程式获得所述外方位元素方程中各个系数的误差量；A12、根据所述外方位元素方程中各个系数的误差量以及所述误差方程，计算所述第一幅卫星影像中各个像点的坐标值的误差值；A13、计算所述第一幅卫星影像中各个像点的坐标值的误差值的中误差；A14、当所述中误差小于预设阈值时，依据所述外方位元素方程中各个系数的误差量对所述当前的定轨定资数据中的空间姿态数据进行校正，控制所述卫星中的相机以所述矫正后的定轨定资数据对所述被摄物体进行拍摄，结束；A15、当所述中误差大于或等于预设阈值时，依据所述外方位元素方程中各个系数的误差量对所述外方位元素方程中的系数进行校正，获得矫正后的空间姿态数据，以及矫正后的所述外方位元素方程的系数，返回步骤A3。一种应用于电荷藕合器件图像传感器的卫星影像定位装置，包括：第一确定模块，用于依据外定标获得卫星上的相机的每一扫描行的定轨定资数据，所述定轨定资数据包括所述相机的摄影中心的空间坐标值和所述相机的摄影中心相对于像片面的空间姿态数据；第一获取模块，用于获取所述相机拍摄同一被摄物体的至少两幅卫星影像，将所述至少两幅卫星影像中分别与所述被摄物体的同一点对应的像点的集合称为连接点；第一构建模块，用于获取所述至少两幅卫星影像的影像条带，以所述相机的扫描方向为横坐标方向，以所述相机的飞行方向作为纵坐标方向，构建以所述连接点的横坐标为自变量和以所述连接点的纵坐标为因变量的瞬时构像方程，所述瞬时构像方程包括：当前的定轨定资数据、与所述连接点对应的所述被摄物体上的点的地面坐标值与所述连接点坐标值的关系；第二确定模块，用于对于每一影像条带，根据所述影像条带的起始扫描时间、所述影像条带的终止扫描时间以及所述影像条带的像素个数，确定所述影像条带中每一像素的摄影时刻；第二构建模块，用于以所述当前的空间姿态数据为因变量，所述摄影时刻为自变量，构建描述所述当前空间姿态数据与所述摄影时刻之间的关系的外方位元素方程，所述外方位元素方程的系数为矫正后的系数，所述外方位元素方程的系数的初始值由所述第一确定模块中的定轨定资数据确定；整理模块，用于依据所述外方位元素方程，将所述瞬时构像方程整理成以所述地面坐标值为变量的共线方程；第一计算模块，用于对于每一连接点，将所述连接点中各个像点的坐标值带入所述共线方程，计算出所述被摄物体中与所述连接点对应的点的地面坐标值；选择模块，用于选择所述至少两幅卫星影像中的第一幅卫星影像；第三构建模块，用于将所述第一幅卫星影像中各个像点分别对应的地面坐标值作为确定值，对于所述第一幅卫星影像中每一像点，将所述第一幅卫星影像中各个像点的坐标值的误差值以及外方位元素方程中各个系数的误差量带入所述共线方程，获得以所述第一幅卫星影像中各个像点的坐标值的误差值作为因变量，以所述外方位元素方程中各个系数的误差量为自变量的误差方程；第二获取模块，用于依据最小二乘间接平差方法，获取所述误差方程的法方程式，所述法方程式包括由所述第一幅卫星影像中各个像点的坐标值的观测权值组成的观测权值矩阵，所述观测权值矩阵的初始值为单位矩阵，所述观测权值矩阵与由前一次所述第一幅卫星影像中各个像点的坐标值的误差值组成的误差值矩阵有关；第三获取模块，用于通过所述法方程式获得所述外方位元素方程中各个系数的误差量；第二计算模块，用于根据所述外方位元素方程中各个系数的误差量以及所述误差方程，计算所述第一幅卫星影像中各个像点的坐标值的误差值；第三计算模块，用于计算所述第一幅卫星影像中各个像点的坐标值的误差值的中误差；矫正模块，用于当所述中误差小于预设阈值时，依据所述外方位元素方程中各个系数的误差量对所述当前的定轨定资数据中的空间姿态数据进行校正，控制所述卫星中的相机以所述矫正后的定轨定资数据对所述被摄物体进行拍摄，结束；触发模块，用于当所述中误差大于或等于预设阈值时，依据所述外方位元素方程中各个系数的误差量对所述外方位元素方程中的系数进行校正，获得矫正后的空间姿态数据，以及矫正后的所述外方位元素方程的系数，触发所述第一构建模块。上述技术方案具有如下有益效果：本专利技术实施例提供的一种应用于电荷藕合器件图像传感器的卫星影像定位方法中，可以实现在没有地面控制点条件下，仅通过外定标获得的定轨定资数据以及卫星的相机在该定轨定资数据下拍摄的卫星影像数据，对卫星的定轨定资数据中的空间姿态数据进行校正。具体的，利用定轨本文档来自技高网...

【技术保护点】
一种应用于电荷藕合器件图像传感器的卫星影像定位方法，其特征在于，包括：A1、依据外定标获得卫星上的相机的每一扫描行的定轨定资数据，所述定轨定资数据包括所述相机的摄影中心的空间坐标值和所述相机的摄影中心相对于像片面的空间姿态数据；A2、获取所述相机拍摄同一被摄物体的至少两幅卫星影像，将所述至少两幅卫星影像中分别与所述被摄物体的同一点对应的像点的集合称为连接点；A3、获取所述至少两幅卫星影像的影像条带，以所述相机的扫描方向为横坐标方向，以所述相机的飞行方向作为纵坐标方向，构建以所述连接点的横坐标为自变量和以所述连接点的纵坐标为因变量的瞬时构像方程，所述瞬时构像方程包括：当前的定轨定资数据、与所述连接点对应的所述被摄物体上的点的地面坐标值与所述连接点坐标值的关系；A4、对于每一影像条带，根据所述影像条带的起始扫描时间、所述影像条带的终止扫描时间以及所述影像条带的像素个数，确定所述影像条带中每一像素的摄影时刻；A5、以所述当前的空间姿态数据为因变量，所述摄影时刻为自变量，构建描述所述当前空间姿态数据与所述摄影时刻之间的关系的外方位元素方程，所述外方位元素方程的系数为矫正后的系数，所述外方位元素方程的系数的初始值由步骤A1中的定轨定资数据确定；A6、依据所述外方位元素方程，将所述瞬时构像方程整理成以所述地面坐标值为变量的共线方程；A7、对于每一连接点，将所述连接点中各个像点的坐标值带入所述共线方程，计算出所述被摄物体中与所述连接点对应的点的地面坐标值；A8、选择所述至少两幅卫星影像中的第一幅卫星影像；A9、将所述第一幅卫星影像中各个像点分别对应的地面坐标值作为确定值，对于所述第一幅卫星影像中每一像点，将所述第一幅卫星影像中各个像点的坐标值的误差值以及外方位元素方程中各个系数的误差量带入所述共线方程，获得以所述第一幅卫星影像中各个像点的坐标值的误差值作为因变量，以所述外方位元素方程中各个系数的误差量为自变量的误差方程；A10、依据最小二乘间接平差方法，获取所述误差方程的法方程式，所述法方程式包括由所述第一幅卫星影像中各个像点的坐标值的观测权值组成的观测权值矩阵，所述观测权值矩阵的初始值为单位矩阵，所述观测权值矩阵与由前一次所述第一幅卫星影像中各个像点的坐标值的误差值组成的误差值矩阵有关；A11、通过所述法方程式获得所述外方位元素方程中各个系数的误差量；A12、根据所述外方位元素方程中各个系数的误差量以及所述误差方程，计算所述第一幅卫星影像中各个像点的坐标值的误差值；A13、计算所述第一幅卫星影像中各个像点的坐标值的误差值的中误差；A14、当所述中误差小于预设阈值时，依据所述外方位元素方程中各个系数的误差量对所述当前的定轨定资数据中的空间姿态数据进行校正，控制所述卫星中的相机以所述矫正后的定轨定资数据对所述被摄物体进行拍摄，结束；A15、当所述中误差大于或等于预设阈值时，依据所述外方位元素方程中各个系数的误差量对所述外方位元素方程中的系数进行校正，获得矫正后的空间姿态数据，以及矫正后的所述外方位元素方程的系数，返回步骤A3。...

【技术特征摘要】
1.一种应用于电荷藕合器件图像传感器的卫星影像定位方法，其特征在于，包括：A1、依据外定标获得卫星上的相机的每一扫描行的定轨定资数据，所述定轨定资数据包括所述相机的摄影中心的空间坐标值和所述相机的摄影中心相对于像片面的空间姿态数据；A2、获取所述相机拍摄同一被摄物体的至少两幅卫星影像，将所述至少两幅卫星影像中分别与所述被摄物体的同一点对应的像点的集合称为连接点；A3、获取所述至少两幅卫星影像的影像条带，以所述相机的扫描方向为横坐标方向，以所述相机的飞行方向作为纵坐标方向，构建以所述连接点的横坐标为自变量和以所述连接点的纵坐标为因变量的瞬时构像方程，所述瞬时构像方程包括：当前的定轨定资数据、与所述连接点对应的所述被摄物体上的点的地面坐标值与所述连接点坐标值的关系；A4、对于每一影像条带，根据所述影像条带的起始扫描时间、所述影像条带的终止扫描时间以及所述影像条带的像素个数，确定所述影像条带中每一像素的摄影时刻；A5、以所述当前的空间姿态数据为因变量，所述摄影时刻为自变量，构建描述所述当前空间姿态数据与所述摄影时刻之间的关系的外方位元素方程，所述外方位元素方程的系数为矫正后的系数，所述外方位元素方程的系数的初始值由步骤A1中的定轨定资数据确定；A6、依据所述外方位元素方程，将所述瞬时构像方程整理成以所述地面坐标值为变量的共线方程；A7、对于每一连接点，将所述连接点中各个像点的坐标值带入所述共线方程，计算出所述被摄物体中与所述连接点对应的点的地面坐标值；A8、选择所述至少两幅卫星影像中的第一幅卫星影像；A9、将所述第一幅卫星影像中各个像点分别对应的地面坐标值作为确定值，对于所述第一幅卫星影像中每一像点，将所述第一幅卫星影像中各个像点的坐标值的误差值以及外方位元素方程中各个系数的误差量带入所述共线方程，获得以所述第一幅卫星影像中各个像点的坐标值的误差值作为因变量，以所述外方位元素方程中各个系数的误差量为自变量的误差方程；A10、依据最小二乘间接平差方法，获取所述误差方程的法方程式，所述法方程式包括由所述第一幅卫星影像中各个像点的坐标值的观测权值组成的观测权值矩阵，所述观测权值矩阵的初始值为单位矩阵，所述观测权值矩阵与由前一次所述第一幅卫星影像中各个像点的坐标值的误差值组成的误差值矩阵有关；A11、通过所述法方程式获得所述外方位元素方程中各个系数的误差量；A12、根据所述外方位元素方程中各个系数的误差量以及所述误差方程，计算所述第一幅卫星影像中各个像点的坐标值的误差值；A13、计算所述第一幅卫星影像中各个像点的坐标值的误差值的中误差；A14、当所述中误差小于预设阈值时，依据所述外方位元素方程中各个系数的误差量对所述当前的定轨定资数据中的空间姿态数据进行校正，控制所述卫星中的相机以所述矫正后的定轨定资数据对所述被摄物体进行拍摄，结束；A15、当所述中误差大于或等于预设阈值时，依据所述外方位元素方程中各个系数的误差量对所述外方位元素方程中的系数进行校正，获得矫正后的空间姿态数据，以及矫正后的所述外方位元素方程的系数，返回步骤A3。2.根据权利要求1所述一种应用于电荷藕合器件图像传感器的卫星影像定位方法，其特征在于，所述定轨定姿数据为其中，i表示第i扫描行，所述相机的摄影中心的空间坐标值为(XSi,YSi,ZSi)，所述相机拍摄被摄物体的空间姿态数据步骤A3具体为：以所述相机的扫描方向为横坐标方向，以所述相机的飞行方向作为纵坐标方向，构建以所述连接点的横坐标为自变量和以所述连接点的纵坐标为因变量的瞬时构像方程：其中，f为所述相机的摄影中心到像片面的垂距，(x0,y0)为像主点坐标值；b1＝cosωisinκi；b2＝cosωicosκi；b3＝-sinωi；3.根据权利要求2所述一种应用于电荷藕合器件图像传感器的卫星影像定位方法，其特征在于，步骤A5具体为：以当前的空间姿态数据为因变量，所述摄影时刻为自变量，构建描述当前空间姿态数据与所述摄影时刻之间的关系的外方位元素方程：其中，q0,q1,q2,m0,m1,m2,n0,n1,n2为二次多项式系数，t为摄影时刻。4.根据权利要求3所述一种应用于电荷藕合器件图像传感器的卫星影像定位方法，其特征在于，步骤A6具体为：依据所述外方位元素方程以及所述摄影时刻，确定所述连接点对应的空间姿态数据将所述瞬时构像方程整理为：其中，l1＝fa1+(x-x0)a3,l2＝fb1+(x-x0)b3,l3＝fc1+(x-x0)c3lx＝fa1XSi+fb1YSi+fc1ZSi+(x-x0)a3XSi+(x-x0)b3YSi+(x-x0)c3ZSil4＝fa2+(y-y0)a3,l5＝fb2+(y-y0)b3,l6＝fc2+(y-y0)c3ly＝fa2XSi+fb2YSi+fc2ZSi+(y-y0)a3XSi+(y-y0)b3YSi+(y-y0)c3ZSi。5.根据权利要求4所述一种应用于电荷藕合器件图像传感器的卫星影像定位方法，其特征在于，步骤A9具体为：将所述第一幅卫星影像中各个像点对应的地面坐标值(X,Y,Z)作为确定值，把所述第一幅卫星影像中各个像点(x,y)视为变量，将所述x修正为x+vx，将所述y修正为y+vy，将q0修正为q0+Δq0，将q1修正为q1+Δq1，将q2修正为q2+Δq2，将m0修正为m0+Δm0，将m1修正为m1+Δm1，将m2修正为m2+Δm2，将n0修正为n0+Δn0，将n1修正n1+Δn1，将n2修正为n2+Δn2，将(x+vx，y+vy)、q0+Δq0、q1+Δq1、q2+Δq2、m0+Δm0、m1+Δm1、m2+Δm2、n0+Δn0、n1+Δn1、n2+Δn2带入所述共线方程，获得误差方程：V＝AT-l，其中，v＝[vxvy]T，T＝[Δq0Δq1Δq2Δm0Δm1Δm2Δn0Δn1Δn2]T，l＝[lxly]T，A矩阵为所述瞬时构像方程对所述外方位元素方程的系数求偏导获得的；步骤A10具体为：依据最小二乘间接平差方法，获取所述误差方程的法方程式：ATPAT＝ATPL，其中，P矩阵为由所述第一幅卫星影像中各个像点的坐标值的观测权值组成的观测权值矩阵，P矩阵中各个观测权值的横坐标方向的权值Px为，vx为前一次所述第一幅卫星影像中该像点的横坐标的误差值，P矩阵中各个观测权值的纵坐标方向的权值Py为，Py为前一次所述第一幅卫星影像中该像点的纵坐标的误差值，k0＝σ0；k1＝3σ0，初始权矩阵P＝E，σ0为所述第一幅卫星影像中所有像点的误差值的中误差。6.一种应用于电荷藕合器件图像传感器的卫星影像定位装置，其特征在于，包括：第一确定模块，用于依据外...

【专利技术属性】
技术研发人员：陈刚，靳笑琳，王一，赵颖芝，折晓宇，
申请(专利权)人：航天恒星科技有限公司，
类型：发明
国别省市：北京;11