CN109800380A - 星载微波遥感仪器对地探测的严密成像几何模型构建方法 - Google Patents
星载微波遥感仪器对地探测的严密成像几何模型构建方法 Download PDFInfo
- Publication number
- CN109800380A CN109800380A CN201910092544.8A CN201910092544A CN109800380A CN 109800380 A CN109800380 A CN 109800380A CN 201910092544 A CN201910092544 A CN 201910092544A CN 109800380 A CN109800380 A CN 109800380A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- reflecting surface
- sight
- axis
- coordinate
- 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
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种遥感仪器建模方法技术领域的星载微波遥感仪器对地探测的严密成像几何模型构建方法,包括以下步骤:步骤一,定义仪器中各独立组件的固连坐标系;步骤二,定义各反射面和馈源的旋转和位移变形参数;步骤三,定义馈源出射视线状态向量集合;步骤四,得到出射视线状态向量集合;步骤五,构建指向探测目标的视向量簇;步骤六,将视向量簇转换到地心固连坐标系下;步骤七,计算视向量簇与地表的交点坐标;步骤九,计算交点的地理经纬度;本发明的光路模型适于含有平面、旋转抛物面和旋转双曲面反射面等组件的微波遥感仪器,可以描述星载微波遥感仪器在结构变形情况下光路从馈源到对地指向波束地理经纬度之间的传播关系。
Description
技术领域
本发明涉及遥感仪器建模方法技术领域,具体涉及一种星载微波遥感仪器对地探测的严密成像几何模型构建方法。
背景技术
图像定位配准直接反映了遥感图像信息与目标之间的空间对应关系,对卫星业务图像产品定量应用有着重要作用,如对区域情况的精确定位,准确跟踪目标以及生成云图动画等。随着遥感卫星图像的分辨率越来越高,其应用也越来越广泛和深入,对图像定位精度的要求也越来越高。遥感卫星图像定位方法很多,按原理可分为与遥感仪器模型无关的几何校正和基于遥感仪器严密成像几何模型的几何校正。
在遥感卫星对地观测的初期,由于无法获取对地遥感仪器在成像瞬间的较高精度的位置、姿态等信息,无法建立共线方程等严密的传感器成像模型,因此大都采用与遥感仪器模型无关的静态参数几何校正。该方法是在参考影像(地形图)和卫星遥感影像上选取大量的对应点(线)作为建立校正参数模型的控制点(线),选择一定的模型,如多项式模型、直接线性变换、平行光投影模型、有理函数模型(RFM)等进行影像校正,从而达到遥感卫星影像的定位功能。
随着高精度定轨和定姿技术的发展,基于遥感仪器严密成像几何模型的几何校正方法逐渐成为高分辨率卫星图像产品生产的主要方法。该方法的主要思想是:在遥感仪器的内方位元素(如光学相机的主距、主点位置、畸变参数)和外方位元素(星载仪器位置和姿态)信息已知的情况下,利用摄影光线与数字高程模型(DEM)相交获取对应的地面坐标,或采用立体模式的空间前方交会获取对应的地面坐标,从而实现卫星遥感影像对地定位。
卫星严密成像几何模型的构建是基于遥感仪器成像模型几何校正方案的关键环节。卫星平台姿态、轨道位置和热环境、力学环境存在长周期和短周期变化,使得仪器视线偏离标称方向,导致图像像素与地理位置的对应关系产生偏差。卫星平台的姿态指向变化、轨道位置变化均可精确测定后进行补偿或校正。而热环境、力学环境等因素引起的卫星平台、仪器在轨变形机理复杂,难以直接精确测量每个组件的变形,不同组件变形对图像的影响各不相同,因此各国遥感卫星制造商均针对遥感仪器的光路特性设计相应的图像定位配准和补偿方案,如美国GOES-I~M系列,GOES-NOP系列卫星,俄罗斯的Electro-L,日本的Himawari-8,欧洲第三代气象卫星MTG,SPOT-5、ALOS、IKONOS、QUICKBIRD,以及我国的资源三号、FY-4等卫星。
但是上述卫星的基于遥感仪器严密成像模型的几何校正方案均是针对光学遥感仪器,不适用于微波遥感仪器。作为我国乃至全球未来气象预报体系的重要组成部分,我国率先部署静止轨道微波遥感卫星。其中微波探测仪是主载荷。由于静止轨道微波探测频段向高频扩展,为满足空间分辨率要求,其天线反射面口径远大于常规天线口径,受限于火箭整流罩尺寸,需要将天线先折叠起来,入轨后再展开完成部署。展开动作可能造成微波各反射面的相对角度和位置偏离标称值。其次,卫星在轨运行时,展开机构的挠性振动导致反射面的角度和位置偏离标称值。与光学遥感卫星不同的是,微波遥感天线反射面不仅由平面反射面还由曲面反射面构成。平面反射面的位移不改变光路反射方向,因此光学遥感星不需要对反射面位移误差建模。曲面反射面位移直接导致光路方向变化,进而改变对地观测视向量指向,因此在微波遥感卫星光路建模中,不可忽略反射面位移误差。最后,相比于光学遥感卫星,微波遥感卫星反射面数量更多,视向量需要经过多重反射后才进入成像系统。
综上各种因素可知,当前的考虑仪器变形的卫星图像定位配准研究的研究对象多为光学遥感仪器。而微波遥感仪器成像原理不同与光学遥感仪器,成像过程比光学遥感仪器更加复杂,具有一定特殊性。为实现微波遥感仪器的高精度图像定位配准,需针对微波遥感仪器的特点,建立适用于星载微波遥感仪器对地探测的严密成像几何模型。
经对现有技术的检索,J.L.Fiorello等人编写的NOAA研究报告(0989年,编号N90-13422)介绍了图像定位配准概念和原理,提出了通过光学仪器观测恒星获取仪器变形参数的方法,但该方法不适用于微波遥感仪器。
吕旺的博士学位论文《静止气象卫星成像导航配准研究》(2017年)以风云四号辐射成像仪作为研究对象,进行了图像导航配准研究。针对由热环境、应力等因素引起的遥感仪器在轨变形问题,提出双模型建模方法。用物理模型详细描述机械变形对视向量的影响,用等效模型作为补偿算法。但该方法仅对光学遥感仪器有效。
专利号CN 104764443 A的发明专利《一种光学遥感卫星严密成像几何模型构建方法》介绍了一种利用卫星相机内外方位元素构建光学遥感卫星影像的的严密成像几何模型的方法。但该方法不适用于微波遥感仪器。
张过、祝彦敏、费文波、李德仁在《测绘通报》第五期(2009年)的文章《高分辨SAR-GEC影像严密成像几何模型及其应用研究》中,提出了构建SAR影像GEC产品严密成像几何模型的方法。但是该方法对仪器内部畸变做了简化处理,无法充分描述扫描式微波成像仪器的变形。
因此,有必要设计一种能够针对微波遥感仪器的特点,建立适用于星载微波遥感仪器对地探测的严密成像几何模型,实现微波遥感仪器的高精度图像定位配准的星载微波遥感仪器对地探测的严密成像几何模型构建方法。
发明内容
针对现有技术的缺陷,本发明的目的是提供一种星载微波遥感仪器对地探测的严密成像几何模型构建方法,本发明能够针对微波遥感仪器的特点,建立适用于星载微波遥感仪器对地探测的严密成像几何模型,实现微波遥感仪器的高精度图像定位配准。
本发明涉及一种星载微波遥感仪器对地探测的严密成像几何模型构建方法,将微波遥感仪器的独立组件分为平面反射面、旋转抛物面反射面、旋转双曲面反射面、和馈源共四类独立组件,采用视线状态向量表征微波波束视线,依据视线反射传输的路径顺序,先通过坐标转换矩阵实现视线状态向量在相邻反射面固连坐标系之间的转换,再通过反射面的状态方程对视线状态向量做变换,以描述光路的反射传输。
进一步地,采用视线状态向量表征微波波束视线,所述视线状态向量,视线状态向量X是7维向量,由三维的视线落点坐标[x,y,z]T、三维的视线方向单位矢量[u,v,w]T和1构成:
X=[x y z u v w 1]T。
进一步地,依据视线反射传输的路径顺序,通过坐标转换矩阵实现视线状态向量在相邻反射面固连坐标系之间的转换。
进一步地,所述坐标转换矩阵,以a坐标系到b坐标系的坐标转换矩阵为例,视线状态向量X在a坐标系的投影Xa,左乘a坐标系到b坐标系的坐标转换矩阵,可得到在b坐标系的投影Xb:
Xb=TbaXa
所述Ab表示a坐标系原点在b坐标系的投影,所述a坐标系到b坐标系的坐标转换矩阵,是7×7阶矩阵,形式如下:
所述Rba表示a坐标系到b坐标系的方向余弦矩阵。
进一步地,在坐标转换矩阵计算过程中,考虑光路组件的旋转变形参数和位移变形参数,以a坐标系到b坐标系的坐标转换矩阵计算为例,计算方向余弦矩阵时候需考虑a坐标系旋转变形参数θa,ψa、b坐标系旋转变形参数θb,ψb,计算a坐标系原点在b坐标系的投影时候需考虑a坐标系位移变形参数ΔA、b坐标系的位移变形参数ΔB,即具有如下形式:
所述表示标称状态下,a坐标系到b坐标系的方向余弦矩阵。
进一步地,所述状态转移方程,在当前反射面固连坐标系下,当前反射面的出射视线状态向量Xout等于当前反射面的状态转移矩阵M左乘当前反射面的入射视线状态向量Xin:
Xout=MXin。
进一步地,所述状态转移矩阵的表达式为:
其中n=[nx,ny,nz]T是视线落点处的反射面外法向单位矢量在反射面固连坐标系的投影,k是入射视线单位方向矢量与出射视线落点坐标的耦合系数,表达式如下:
当反射面为平面z=0时:
k=-zin/win
当反射面为旋转抛物面x2+y2+az=0,a<0时:
当反射面为旋转双曲面x2+y2+az2=0,z≤0,a<0时:
进一步地,根据卫星姿态信息、轨道信息和岁差、章动、格林尼治恒星时和极移信息,将仪器指向探测目标的视向量簇,转换到地球固连坐标系下,求解视向量簇与地表交点的经纬度,视线波束足迹的地理定位。
进一步地,忽略地形对地理定位的影响,以WGS-84地球参考椭球表示地球,以视向量簇与WGS-84地球参考椭球表面的交点表示波束足迹。
本发明的星载微波遥感仪器对地探测的严密成像几何模型构建方法,包括如下步骤:
步骤1,将静止轨道微波遥感仪器中的独立光路组件划分为平面反射面、旋转抛物面反射面、旋转双曲面反射面和馈源共四种类型;根据独立组件的类型,定义与独立组件固连的平面反射面固连坐标系、旋转抛物面反射面固连坐标系、旋转双曲面反射面固连坐标系、旋转轴系固连坐标系和馈源固连坐标系;
所述平面反射面固连坐标系,坐标系原点取为反射面几何中心,z轴为平面正法向,x轴在平面内指向一特征方向,y轴由右手定则确定;
所述旋转抛物面反射面固连坐标系,坐标系原点为旋转抛物面反射面的拐点,z轴与旋转抛物面反射面对称轴重合指向旋转抛物面拐点正法向方向,x轴垂直z轴指向旋转抛物面反射面一特征方向,y轴由右手定则确定;
所述旋转双曲面反射面固连坐标系,坐标系原点为旋转双曲面反射面的拐点,z轴与旋转双曲面反射面的对称轴重合指向旋转双曲面反射面拐点正法向方向,x轴垂直z轴指向旋转双曲面反射面某特征方向,y轴由右手定则确定;
所述馈源固连坐标系,馈源固连坐标系原点取为馈源入口几何中心,z轴垂直馈源入口平面指向馈源外侧,x轴垂直z轴指向馈源入口平面一特征方向,y轴由右手定则确定;
步骤2,定义平面反射面、旋转抛物面反射面、旋转双曲面反射面的旋转变形参数和位移变形参数,
所述平面反射面的旋转变形参数和位移变形参数,包括平面反射面绕自身固连坐标系X轴的旋转角、绕Y轴的旋转角,和沿Z轴的位移,
所述旋转抛物面反射面的旋转变形参数和位移变形参数,包括旋转抛物面反射面绕自身固连坐标系X轴的旋转角、绕Y轴的旋转角,和沿X轴的位移、沿Y轴的位移、沿Z轴的位移,
所述旋转双曲面反射面的旋转变形参数和位移变形参数,包括旋转双曲面反射面绕自身固连坐标系X轴的旋转角、绕Y轴的旋转角,和沿X轴的位移、沿Y轴的位移、沿Z轴的位移,
步骤3,定义从馈源出发的波束、视线矢量簇、出射视线状态向量集合,所述从馈源出发的波束、视线矢量簇、视线状态向量集合,波束从馈源固连坐标系原点出发,呈圆锥状,圆锥顶点与馈源固连坐标系原点重合,圆锥中轴与馈源系外法向重合,圆锥面为馈源入口处波束方向图-3dB包络渐近线,视线矢量簇由波束中心视线矢量和等间隔的波束边界视线矢量组成;出射视线状态向量集合是视线矢量簇对应的视线状态向量的集合;
步骤4,从馈源出发,考虑光路部件旋转变形参数和位移变形参数计算相邻光路部件之间的坐标系转换矩阵,应用坐标转换矩阵,将出射视线状态向量集合中的所有状态向量转换到光路中下一个反射面的固连坐标系下,作为该反射面的入射视线状态向量集合,根据该反射面的类型采用相应的状态转移矩阵,将入射视线状态向量集合中的所有状态向量分别代入状态转移方程,计算得到出射视线状态向量集合,循环此过程,直至获得仪器最外端指向探测目标的出射视线状态向量集合;
步骤5,根据获得仪器最外端指向探测目标的出射视线状态向量集合构建仪器最外端指向探测目标的视向量簇,并根据仪器最外端反射面的安装矩阵将视向量簇转换到卫星本体坐标系下;
步骤6,根据卫星姿态信息、轨道信息和岁差、章动、格林尼治恒星时、极移等信息,将仪器光路模型输出的本体系下视向量簇转换到地心固连坐标系下;
步骤7,在地心固连坐标系下,计算视向量簇与以WGS-84地球参考椭球表示的地表的交点坐标;
步骤8,计算交点的地理经纬度。
与现有技术相比,本发明具有如下的有益效果:
1、本发明的星载微波遥感仪器对地探测的严密成像几何模型构建方法,适用于含有旋转抛物面反射面、旋转双曲面反射面、平面反射面以及旋转轴系等多种机构的微波遥感仪器;
2、本发明的星载微波遥感仪器对地探测的严密成像几何模型构建方法,所建立的模型可以描述星载微波遥感仪器中光路从馈源到对地指向波束之间的传播关系;
3、本发明的星载微波遥感仪器对地探测的严密成像几何模型构建方法,能够针对微波遥感仪器的特点,建立适用于星载微波遥感仪器对地探测的严密成像几何模型,实现微波遥感仪器的高精度图像定位配准。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明的星载微波遥感仪器结构变形的光路建模方法流程图;
图2为本发明的馈源和出射波束示意图;
图3为本发明的平面反射面光路建模示意图;
图4为本发明的旋转双曲面反射面光路建模示意图;
图5为本发明的旋转抛物面反射面光路建模示意图;
图6为某星载微波遥感仪器内部光路示意图;
图7为本发明的相关坐标系示意图;
图8为某星载微波遥感仪器视向量簇与WGS-84地球参考椭圆交点示意图;
图9为某星载微波遥感仪器波束地面足迹图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
实施例
本实施例中,本发明的星载微波遥感仪器对地探测的严密成像几何模型构建方法,将微波遥感仪器的独立组件分为平面反射面、旋转抛物面反射面、旋转双曲面反射面、和馈源共四类独立组件,采用视线状态向量表征微波波束视线,依据视线反射传输的路径顺序,先通过坐标转换矩阵实现视线状态向量在相邻反射面固连坐标系之间的转换,再通过反射面的状态方程对视线状态向量做变换,以描述光路的反射传输。
接下来对本发明进行详细的描述。
本发明的目的是提供一种星载微波遥感仪器对地探测的严密成像几何模型构建方法,本发明能够针对微波遥感仪器的特点,建立适用于星载微波遥感仪器对地探测的严密成像几何模型,实现微波遥感仪器的高精度图像定位配准。
如图1至图7所示,本发明一较佳实施例的静止轨道微波遥感仪器光路建模方法,包括如下步骤:
步骤1,将静止轨道微波遥感仪器中的独立光路组件划分为平面反射面、旋转抛物面反射面、旋转双曲面反射面、旋转轴系和馈源共五种类型;根据独立组件的类型,定义与独立组件固连的平面反射面固连坐标系、旋转抛物面反射面固连坐标系、旋转双曲面反射面固连坐标系、旋转轴系固连坐标系和馈源固连坐标系。
定义各部件的固连坐标系如下:
馈源入口固连坐标系(下标a,坐标原点记为A):原点A取为馈源入口几何中心,Za轴取为馈源入口外法向,Xa轴垂直Za指向馈源某特征方向,Ya由右手定则确定;
快扫镜一是平面反射面,快扫镜一固连坐标系(下标b,坐标原点记为B):原点B取为快扫镜几何中心,Zb轴取为快扫镜正法向,Xb轴垂直Zb轴指向某特征方向,Yb轴由右手定则确定;
快扫镜而是平面反射面,快扫镜二固连坐标系(下标c,坐标原点记为C):原点C取为快扫镜几何中心,Zc轴取为快扫镜二正法向,Xc轴垂直Zc轴指向某特征方向,Yc轴由右手定则确定;;
第二副反是平面反射面,第二副反固连坐标系(下标d,坐标原点记为D):原点D取为第二副反几何中心,Zd轴取为第二副反正法向,Xd轴垂直Zd轴指向某特征方向,Yd轴由右手定则确定;
第一副反是旋转双曲面反射面,第一副反固连坐标系(下标e,坐标原点记为E):原点E取为第一副反拐点,Ze轴取为与第一副反对称轴重合,指向外法向方向,Xe轴垂直Zd轴指向某特征方向,Ye轴由右手定则确定;;
主反是旋转抛物面反射面主反固连坐标系(下标f,坐标原点记为F):原点F取为主反拐点,Zf轴取为与主反对称轴重合,指向主反外法向,Xf轴垂直Zf轴指向某特征方向,Yf轴由右手定则确定;
步骤2,定义平面反射面、旋转抛物面反射面、旋转双曲面反射面旋转变形参数和位移变形参数。
快扫镜一是平面反射面,标称快扫镜一固连坐标系按照XY转序旋转欧拉角θb,然后沿旋转后的固连坐标系Zb轴平移zb,得到实际快扫镜一固连坐标系。则变形前快扫镜一固连坐标系到变形后快扫镜一固连坐标系的方向余弦矩阵
快扫镜二是平面反射面,标称快扫镜二固连坐标系按照XY转序旋转欧拉角θc,然后沿旋转后的固连坐标系Zc轴平移zc,得到实际快扫镜二固连坐标系。则变形前快扫镜二固连坐标系到变形后快扫镜二固连坐标系的方向余弦矩阵
第二副反是平面反射面,标称第二副反固连坐标系按照XY转序旋转欧拉角θd,然后沿旋转后的固连坐标系Zd平移zd,得到实际第二副反固连坐标系。则变形前第二副反射面固连坐标系到变形后第二副反射面固连坐标系的方向余弦矩阵
第一副反是旋转双曲面反射面,标称第一副反固连坐标系按照XY转序旋转欧拉角θe,然后沿旋转后的固连坐标系平移xe,ye,ze,得到实际第一副反固连坐标系。则变形前第一副反射面固连坐标系到变形后第一副反射面固连坐标系的方向余弦矩阵
主反射面是旋转抛物面反射面,标称主反固连坐标系按照XY转序旋转欧拉角θf,然后沿旋转后的固连坐标系平移xf,yf,zf,得到实际第一副反固连坐标系。则变形前主反射面固连坐标系到变形后主反射面固连坐标系的方向余弦矩阵
步骤3,定义从馈源出发的波束、视线矢量簇、出射视线状态向量集合。
所述从馈源出发的波束、视线矢量簇、视线状态向量集合,其特征在于,波束从馈源固连坐标系原点出发,呈圆锥状,圆锥顶点与馈源固连坐标系原点重合,圆锥中轴与馈源系外法向重合,圆锥面为馈源入口处波束方向图-3dB包络渐近线;视线矢量簇由波束中心视线矢量和等间隔的波束边界视线矢量组成;出射视线状态向量集合{Xout_a}是视线矢量簇对应的视线状态向量Xout_a的集合。
步骤4,从馈源出发,考虑光路部件旋转变形参数和位移变形参数,计算相邻光路部件之间的坐标系转换矩阵,应用坐标转换矩阵,将出射视线状态向量集合中的所有状态向量转换到光路中下一个反射面的固连坐标系下,作为该反射面的入射视线状态向量集合。根据该反射面的类型采用相应的状态转移矩阵,将入射视线状态向量集合中的所有状态向量分别代入状态转移方程,计算得到出射视线状态向量集合。循环此过程,直至获得仪器最外端指向探测目标的出射视线状态向量集合。
首先,通过坐标变换,将馈源出射视线状态向量集合{Xout_a}转换到b系下{Xin_b}={TbaXout_a},其中:
应用快扫镜一的状态转换矩阵Mb,计算快扫镜一的出射视线状态向量集合:
{Xout_b}={MbXin_b}
然后将Xout_b转化到快扫镜二固连坐标系下,得到快扫镜二的入射视线状态向量集合{Xin_c}={TcbXin_c},其中:
应用快扫镜二的状态转移矩阵Mc,计算快扫镜二的出射视线状态向量集合:
{Xout_c}={McXin_c}
将Xout_c转换到第二副反射面固连坐标系下,得到第二副反射面的入射视线状态向量集合{Xin_d}={TdcXout_c},其中:
应用第二副反射面的状态转移矩阵Md,计算第二副反射面的出射视线状态向量集合:
{Xout_d}={MdXin_d}
将Xout_d转换到第一副反射面固连坐标系下,得到第一副反射面固连坐标系的入射视线状态向量集合{Xin_e}={TedXout_d},其中:
应用第一副反射面的状态转移矩阵Me,计算第一副反射面的出射视线状态向量集合:
{Xout_e}={MeXin_e}
将Xout_e转换到主反射面固连坐标系下,得到主反射面固连坐标系的入射视线状态向量集合{Xin_f}={TfeXout_e},其中:
应用主反射面的状态转移矩阵Mf,计算主反射面的出射视线状态向量集合:
{Xout_f}={MfXin_f}
步骤5,根据获得仪器最外端指向探测目标的出射视线状态向量集合构建仪器最外端指向探测目标的视向量簇,并根据仪器最外端反射面的安装矩阵,将视向量簇转换到卫星本体坐标系下。
根据主反射面出射视线状态向量集合,构建主反射面固连坐标系下表示的视向量簇{Vf}={Xout_f(1:3)}。然后根据反射面安装矩阵Rfm,将视向量簇转换到卫星本体坐标系下
步骤6,根据卫星姿态信息、轨道信息和岁差、章动、格林尼治恒星时、极移等信息,将本体系下视向量簇转换到地心固连坐标系下。
首先根据卫星姿态信息,将仪器光路模型输出的本体系下视向量簇转换到轨道坐标系下。已知轨道系按ZXY转序分别旋转ψ,θ,角得到卫星本体系,则视向量簇{V}在轨道系的投影是:
其次根据卫星轨道信息,将轨道坐标系下的视向量簇转换到地心惯性坐标系下。已知卫星瞬时位置矢量在地心惯性坐标系的投影ri(t)、瞬时速度矢量在地心惯性坐标系的投影vi(t)。则卫星轨道坐标系Zo轴在地心惯性坐标系的投影是Zio(t)=ri(t)/r(t),卫星轨道坐标系Yo轴在地心惯性坐标系的投影是Yio(t)=Zio(t)×v(t)/v(t),卫星轨道坐标系Xo轴在地心惯性坐标系的投影是Xio(t)=Yio(t)×Zio(t)。因此,卫星轨道坐标系到地心惯性坐标系的方向余弦矩阵是Tio=[Xio(t)Yio(t)Zio(t)]。从而,得到视向量簇{V}在地心惯性坐标系下的投影:
{Vi}=Tio{Vo}
最后根据岁差、章动、格林尼治恒星时、极移等信息,将地心惯性坐标系下的视向量簇、卫星位置矢量转换到地心固连坐标系下。
视向量簇:
{VWGS}=TWGSi{Vi}=T4(t)T3(t)T2(t)T1(t){Vi}
卫星位置矢量:
rWGS=TWGSiri=T4(t)T3(t)T2(t)T1(t)ri
其中,下标WGS表示地心固连坐标系,T1(t)表示岁差修正矩阵、T2(t)表示章动修正矩阵、T3(t)表示格林尼治恒星时旋转矩阵、T4(t)表示极移修正矩阵,TWGSi表示地心惯性坐标系到地心固连坐标系的方向余弦矩阵。
步骤7,在地心固连坐标系下,计算视向量簇与以WGS-84地球参考椭球表示的地表的交点坐标。
以WGS-84地球参考椭球表示地表,WGS-84地球参考椭球方程是:
其中a是长轴、b是短轴。视向量簇与地表的交点{pWGS}通过以下公式计算:
{pWGS}=rWGS+{sVWGS}
其中s表示卫星与交点的距离:
步骤8,计算交点的地理经纬度。
根据交点在地球固连坐标系的坐标,计算交点的地理经纬度如下:
综上所述,本发明的星载微波遥感仪器对地探测的严密成像几何模型构建方法,适用于含有旋转抛物面反射面、旋转双曲面反射面、平面反射面以及旋转轴系等多种机构的微波遥感仪器;所建立的模型可以描述星载微波遥感仪器中光路从馈源到对地指向波束之间的传播关系;能够针对微波遥感仪器的特点,建立适用于星载微波遥感仪器对地探测的严密成像几何模型,实现微波遥感仪器的高精度图像定位配准。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (10)
1.一种星载微波遥感仪器对地探测的严密成像几何模型构建方法,其特征在于,将微波遥感仪器的独立组件分为平面反射面、旋转抛物面反射面、旋转双曲面反射面、和馈源共四类独立组件,采用视线状态向量表征微波波束视线,依据视线反射传输的路径顺序,先通过坐标转换矩阵实现视线状态向量在相邻反射面固连坐标系之间的转换,再通过反射面的状态方程对视线状态向量做变换,以描述光路的反射传输。
2.根据权利要求1所述的星载微波遥感仪器对地探测的严密成像几何模型构建方法,其特征在于,采用视线状态向量表征微波波束视线,所述视线状态向量,视线状态向量X是7维向量,由三维的视线落点坐标[x,y,z]T、三维的视线方向单位矢量[u,v,w]T和1构成:
X=[x y z u v w 1]T,其中,X为视线状态向量。
3.根据权利要求1所述的星载微波遥感仪器对地探测的严密成像几何模型构建方法,其特征在于,依据视线反射传输的路径顺序,通过坐标转换矩阵实现视线状态向量在相邻反射面固连坐标系之间的转换。
4.根据权利要求3所述的星载微波遥感仪器对地探测的严密成像几何模型构建方法,其特征在于,所述坐标转换矩阵,以a坐标系到b坐标系的坐标转换矩阵为例,视线状态向量X在a坐标系的投影Xa,左乘a坐标系到b坐标系的坐标转换矩阵,可得到在b坐标系的投影Xb:
Xb=TbaXa
Ab表示a坐标系原点在b坐标系的投影,所述a坐标系到b坐标系的坐标转换矩阵,是7×7阶矩阵,形式如下:
其中,Xa为左乘视线状态向量X在a坐标系的投影,Ab为a坐标系原点在b坐标系的投影,Xb为视线状态向量X在b坐标系的投影,Rba表示a坐标系到b坐标系的坐标转换矩阵,Tba坐标系到b坐标系的坐标转换矩阵。
5.根据权利要求1所述的星载微波遥感仪器对地探测的严密成像几何模型构建方法,其特征在于,在坐标转换矩阵计算过程中,考虑光路组件的旋转变形参数和位移变形参数,以a坐标系到b坐标系的坐标转换矩阵计算为例,计算方向余弦矩阵时候需考虑a坐标系旋转变形参数θa,ψa、b坐标系旋转变形参数θb,ψb,计算a坐标系原点在b坐标系的投影时候需考虑a坐标系位移变形参数ΔA、b坐标系的位移变形参数ΔB,即具有如下形式:
所述表示标称状态下,a坐标系到b坐标系的方向余弦矩阵。
6.根据权利要求1所述的星载微波遥感仪器对地探测的严密成像几何模型构建方法,其特征在于,所述状态转移方程,在当前反射面固连坐标系下,当前反射面的出射视线状态向量Xout等于当前反射面的状态转移矩阵M左乘当前反射面的入射视线状态向量Xin:
Xout=MXin,
其中,Xout为当前为反射面的出射视线状态向量,Xout为当前反射面的状态转移矩阵,Xin为当前反射面的入射视线状态向量。
7.根据权利要求6所述的星载微波遥感仪器对地探测的严密成像几何模型构建方法,其特征在于,所述状态转移矩阵的表达式为:
其中n=[nx,ny,nz]T是视线落点处的反射面外法向单位矢量在反射面固连坐标系的投影,k是入射视线单位方向矢量与出射视线落点坐标的耦合系数,表达式如下:
当反射面为平面z=0时:
k=-zin/win
当反射面为旋转抛物面x2+y2+az=0,a<0时:
当反射面为旋转双曲面x2+y2+az2=0,z≤0,a<0时:
8.根据权利要求1所述的星载微波遥感仪器对地探测的严密成像几何模型构建方法,其特征在于,根据卫星姿态信息、轨道信息和岁差、章动、格林尼治恒星时和极移信息,将仪器指向探测目标的视向量簇,转换到地球固连坐标系下,求解视向量簇与地表交点的经纬度,视线波束足迹的地理定位。
9.根据权利要求8所述的星载微波遥感仪器对地探测的严密成像几何模型构建方法,其特征在于,忽略地形对地理定位的影响,以WGS-84地球参考椭球表示地球,以视向量簇与WGS-84地球参考椭球表面的交点表示波束足迹。
10.根据权利要求1所述的星载微波遥感仪器对地探测的严密成像几何模型构建方法,其特征在于,包括如下步骤:
步骤1,将静止轨道微波遥感仪器中的独立光路组件划分为平面反射面、旋转抛物面反射面、旋转双曲面反射面和馈源共四种类型;根据独立组件的类型,定义与独立组件固连的平面反射面固连坐标系、旋转抛物面反射面固连坐标系、旋转双曲面反射面固连坐标系、旋转轴系固连坐标系和馈源固连坐标系;
所述平面反射面固连坐标系,坐标系原点取为反射面几何中心,z轴为平面正法向,x轴在平面内指向一特征方向,y轴由右手定则确定;
所述旋转抛物面反射面固连坐标系,坐标系原点为旋转抛物面反射面的拐点,z轴与旋转抛物面反射面对称轴重合指向旋转抛物面拐点正法向方向,x轴垂直z轴指向旋转抛物面反射面一特征方向,y轴由右手定则确定;
所述旋转双曲面反射面固连坐标系,坐标系原点为旋转双曲面反射面的拐点,z轴与旋转双曲面反射面的对称轴重合指向旋转双曲面反射面拐点正法向方向,x轴垂直z轴指向旋转双曲面反射面某特征方向,y轴由右手定则确定;
所述馈源固连坐标系,馈源固连坐标系原点取为馈源入口几何中心,z轴垂直馈源入口平面指向馈源外侧,x轴垂直z轴指向馈源入口平面一特征方向,y轴由右手定则确定;
步骤2,定义平面反射面、旋转抛物面反射面、旋转双曲面反射面的旋转变形参数和位移变形参数,
所述平面反射面的旋转变形参数和位移变形参数,包括平面反射面绕自身固连坐标系X轴的旋转角、绕Y轴的旋转角,和沿Z轴的位移,
所述旋转抛物面反射面的旋转变形参数和位移变形参数,包括旋转抛物面反射面绕自身固连坐标系X轴的旋转角、绕Y轴的旋转角,和沿X轴的位移、沿Y轴的位移、沿Z轴的位移,
所述旋转双曲面反射面的旋转变形参数和位移变形参数,包括旋转双曲面反射面绕自身固连坐标系X轴的旋转角、绕Y轴的旋转角,和沿X轴的位移、沿Y轴的位移、沿Z轴的位移,
步骤3,定义从馈源出发的波束、视线矢量簇、出射视线状态向量集合,所述从馈源出发的波束、视线矢量簇、视线状态向量集合,波束从馈源固连坐标系原点出发,呈圆锥状,圆锥顶点与馈源固连坐标系原点重合,圆锥中轴与馈源系外法向重合,圆锥面为馈源入口处波束方向图-3dB包络渐近线,视线矢量簇由波束中心视线矢量和等间隔的波束边界视线矢量组成;出射视线状态向量集合是视线矢量簇对应的视线状态向量的集合;
步骤4,从馈源出发,考虑光路部件旋转变形参数和位移变形参数计算相邻光路部件之间的坐标系转换矩阵,应用坐标转换矩阵,将出射视线状态向量集合中的所有状态向量转换到光路中下一个反射面的固连坐标系下,作为该反射面的入射视线状态向量集合,根据该反射面的类型采用相应的状态转移矩阵,将入射视线状态向量集合中的所有状态向量分别代入状态转移方程,计算得到出射视线状态向量集合,循环此过程,直至获得仪器最外端指向探测目标的出射视线状态向量集合;
步骤5,根据获得仪器最外端指向探测目标的出射视线状态向量集合构建仪器最外端指向探测目标的视向量簇,并根据仪器最外端反射面的安装矩阵将视向量簇转换到卫星本体坐标系下;
步骤6,根据卫星姿态信息、轨道信息和岁差、章动、格林尼治恒星时、极移信息,将仪器光路模型输出的本体系下视向量簇转换到地心固连坐标系下;
步骤7,在地心固连坐标系下,计算视向量簇与以WGS-84地球参考椭球表示的地表的交点坐标;
步骤8,计算交点的地理经纬度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910092544.8A CN109800380B (zh) | 2019-01-30 | 2019-01-30 | 星载微波遥感仪器对地探测的严密成像几何模型构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910092544.8A CN109800380B (zh) | 2019-01-30 | 2019-01-30 | 星载微波遥感仪器对地探测的严密成像几何模型构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109800380A true CN109800380A (zh) | 2019-05-24 |
CN109800380B CN109800380B (zh) | 2023-03-17 |
Family
ID=66560692
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910092544.8A Active CN109800380B (zh) | 2019-01-30 | 2019-01-30 | 星载微波遥感仪器对地探测的严密成像几何模型构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109800380B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110278017A (zh) * | 2019-06-27 | 2019-09-24 | 广东工业大学 | 一种基于智能反射面的多天线无线能量传输系统与方法 |
CN110765944A (zh) * | 2019-10-23 | 2020-02-07 | 长光禹辰信息技术与装备(青岛)有限公司 | 基于多源遥感影像的目标识别方法、装置、设备及介质 |
CN111103001A (zh) * | 2020-01-21 | 2020-05-05 | 上海航天测控通信研究所 | 一种静止轨道微波辐射计运动扫描装置及方法 |
CN111561936A (zh) * | 2020-05-19 | 2020-08-21 | 中国科学院微小卫星创新研究院 | 旋转大幅宽光学卫星精准处理方法及系统 |
CN112068157A (zh) * | 2020-07-30 | 2020-12-11 | 国家卫星气象中心(国家空间天气监测预警中心) | 静止轨道多频太赫兹探测仪对地观测模式实现方法及装置 |
CN113063440A (zh) * | 2021-02-26 | 2021-07-02 | 上海卫星工程研究所 | 静止轨道微波探测卫星图像定位与配准全物理仿真试验方法和系统 |
CN114440885A (zh) * | 2021-12-24 | 2022-05-06 | 中国人民解放军战略支援部队信息工程大学 | 一种静止轨道遥感卫星定位方法及装置 |
CN115276760A (zh) * | 2022-06-24 | 2022-11-01 | 中国联合网络通信集团有限公司 | 波束中心的位置确定方法、装置及计算机存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104764443A (zh) * | 2015-04-24 | 2015-07-08 | 国家测绘地理信息局卫星测绘应用中心 | 一种光学遥感卫星严密成像几何模型构建方法 |
CN106526593A (zh) * | 2016-12-19 | 2017-03-22 | 国家测绘地理信息局卫星测绘应用中心 | 基于sar严密成像模型的子像素级角反射器自动定位方法 |
CN106897962A (zh) * | 2016-12-26 | 2017-06-27 | 中国科学院长春光学精密机械与物理研究所 | 大视场空间对地观测图像的圆盘投影和拼接方法 |
CN107202582A (zh) * | 2017-04-24 | 2017-09-26 | 上海卫星工程研究所 | 静止卫星在轨实时成像导航与配准方法 |
CN207992458U (zh) * | 2017-12-24 | 2018-10-19 | 航天恒星科技有限公司 | 复杂观测模式下碳卫星定位系统 |
CN108828623A (zh) * | 2018-06-15 | 2018-11-16 | 上海卫星工程研究所 | 静止气象卫星成像仪的地球固定网格映射方法 |
CN109116379A (zh) * | 2018-06-15 | 2019-01-01 | 上海卫星工程研究所 | 地表目标在静止气象卫星成像仪图像中的映射方法 |
-
2019
- 2019-01-30 CN CN201910092544.8A patent/CN109800380B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104764443A (zh) * | 2015-04-24 | 2015-07-08 | 国家测绘地理信息局卫星测绘应用中心 | 一种光学遥感卫星严密成像几何模型构建方法 |
CN106526593A (zh) * | 2016-12-19 | 2017-03-22 | 国家测绘地理信息局卫星测绘应用中心 | 基于sar严密成像模型的子像素级角反射器自动定位方法 |
CN106897962A (zh) * | 2016-12-26 | 2017-06-27 | 中国科学院长春光学精密机械与物理研究所 | 大视场空间对地观测图像的圆盘投影和拼接方法 |
CN107202582A (zh) * | 2017-04-24 | 2017-09-26 | 上海卫星工程研究所 | 静止卫星在轨实时成像导航与配准方法 |
CN207992458U (zh) * | 2017-12-24 | 2018-10-19 | 航天恒星科技有限公司 | 复杂观测模式下碳卫星定位系统 |
CN108828623A (zh) * | 2018-06-15 | 2018-11-16 | 上海卫星工程研究所 | 静止气象卫星成像仪的地球固定网格映射方法 |
CN109116379A (zh) * | 2018-06-15 | 2019-01-01 | 上海卫星工程研究所 | 地表目标在静止气象卫星成像仪图像中的映射方法 |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110278017B (zh) * | 2019-06-27 | 2023-07-25 | 广东工业大学 | 一种基于智能反射面的多天线无线能量传输系统与方法 |
CN110278017A (zh) * | 2019-06-27 | 2019-09-24 | 广东工业大学 | 一种基于智能反射面的多天线无线能量传输系统与方法 |
CN110765944A (zh) * | 2019-10-23 | 2020-02-07 | 长光禹辰信息技术与装备(青岛)有限公司 | 基于多源遥感影像的目标识别方法、装置、设备及介质 |
CN111103001A (zh) * | 2020-01-21 | 2020-05-05 | 上海航天测控通信研究所 | 一种静止轨道微波辐射计运动扫描装置及方法 |
CN111103001B (zh) * | 2020-01-21 | 2021-06-22 | 上海航天测控通信研究所 | 一种静止轨道微波辐射计运动扫描装置及方法 |
CN111561936A (zh) * | 2020-05-19 | 2020-08-21 | 中国科学院微小卫星创新研究院 | 旋转大幅宽光学卫星精准处理方法及系统 |
CN112068157A (zh) * | 2020-07-30 | 2020-12-11 | 国家卫星气象中心(国家空间天气监测预警中心) | 静止轨道多频太赫兹探测仪对地观测模式实现方法及装置 |
CN112068157B (zh) * | 2020-07-30 | 2024-04-12 | 国家卫星气象中心(国家空间天气监测预警中心) | 静止轨道多频太赫兹探测仪对地观测模式实现方法及装置 |
CN113063440A (zh) * | 2021-02-26 | 2021-07-02 | 上海卫星工程研究所 | 静止轨道微波探测卫星图像定位与配准全物理仿真试验方法和系统 |
CN114440885B (zh) * | 2021-12-24 | 2023-10-13 | 中国人民解放军战略支援部队信息工程大学 | 一种静止轨道遥感卫星定位方法及装置 |
CN114440885A (zh) * | 2021-12-24 | 2022-05-06 | 中国人民解放军战略支援部队信息工程大学 | 一种静止轨道遥感卫星定位方法及装置 |
CN115276760B (zh) * | 2022-06-24 | 2023-06-23 | 中国联合网络通信集团有限公司 | 波束中心的位置确定方法、装置及计算机存储介质 |
CN115276760A (zh) * | 2022-06-24 | 2022-11-01 | 中国联合网络通信集团有限公司 | 波束中心的位置确定方法、装置及计算机存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109800380B (zh) | 2023-03-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109800380A (zh) | 星载微波遥感仪器对地探测的严密成像几何模型构建方法 | |
US9927513B2 (en) | Method for determining the geographic coordinates of pixels in SAR images | |
CN102323571B (zh) | 综合总体参数的星载双天线sar干涉定标器布放方法 | |
CN109828248B (zh) | 星载微波遥感仪器角度变形对视向量影响的分析方法 | |
CN109975831B (zh) | 星载微波遥感仪器变形的在轨补偿方法 | |
CN101339244A (zh) | 一种机载sar图像自动目标定位方法 | |
CN108828623B (zh) | 静止气象卫星成像仪的地球固定网格映射方法 | |
CN107633108A (zh) | 静止卫星成像导航与配准的轨道运动补偿方法 | |
CN111505608B (zh) | 一种基于星载激光单片足印影像的激光指向在轨检校方法 | |
Tang et al. | Geometric accuracy analysis model of the Ziyuan-3 satellite without GCPs | |
CN109116379B (zh) | 地表目标在静止气象卫星成像仪图像中的映射方法 | |
CN102043142A (zh) | 基于数字聚束的合成孔径雷达极坐标波前弯曲补偿方法 | |
JPH0772244A (ja) | 干渉型合成開口レーダ装置および地形変動観測方式 | |
CN109085587A (zh) | 机载InSAR困难区域的DEM重建方法 | |
CN104122552A (zh) | 一种滑动式双站圆周合成孔径雷达成像方法 | |
CN109800379A (zh) | 星载微波遥感仪器光路建模方法 | |
KR100441590B1 (ko) | 간섭측정용 합성 개구 레이다의 기하학적 특성을 이용하여지형고도를 측정하기 위한 디지털 고도모형 생성방법 | |
US20230152443A1 (en) | Method for georeferencing of a digital elevation model | |
Gentile et al. | Monitoring the bidirectional relativistic jets of the radio galaxy 3C 338 | |
Luomei et al. | Demonstration of simultaneous localization and imaging with multirotor-borne MiniSAR | |
CN111521197B (zh) | 一种摆扫大幅宽光学卫星传感器校正方法 | |
Li et al. | A study of the potential attainable geometric accuracy of IKONOS satellite imagery | |
CN109975832A (zh) | 星载微波遥感仪器结构变形的描述方法 | |
CN107907881B (zh) | 一种长孔径星载sar地形估计方法 | |
CN113532372B (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 |