CN105842259B - 星载盐度计冷空外定标方法和装置 - Google Patents
星载盐度计冷空外定标方法和装置 Download PDFInfo
- Publication number
- CN105842259B CN105842259B CN201610158460.6A CN201610158460A CN105842259B CN 105842259 B CN105842259 B CN 105842259B CN 201610158460 A CN201610158460 A CN 201610158460A CN 105842259 B CN105842259 B CN 105842259B
- Authority
- CN
- China
- Prior art keywords
- antenna
- salinometer
- spaceborne
- satellite
- external calibration
- 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
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N22/00—Investigating or analysing materials by the use of microwaves or radio waves, i.e. electromagnetic waves with a wavelength of one millimetre or more
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供了一种星载盐度计冷空外定标方法和装置，涉及卫星观测领域，所述方法包括：根据星载盐度计的天线背瓣视场内的陆地比重，确定适于冷空外定标的预选区域；根据星载盐度计的天线主瓣对应的预设观测数据，确定适于冷空外定标的时间窗口；其中，观测数据包括：太空辐射、太阳辐射和月亮辐射；根据预选区域和时间窗口对应的冷空外定标卫星实际运行轨道，计算标准天线温度；对标准天线温度和星载盐度计的实测天线温度进行统计处理，得到用于校正盐度计观测亮温偏差的冷空外定标系数；其能够定期对盐度计进行外定标以获得仪器定标系数，用于纠正盐度计获取的观测亮温的偏差，并且，上述方法系统的解决了我国星载盐度计冷空定标的问题。
Description
技术领域
本发明涉及卫星观测领域，具体而言，涉及星载盐度计冷空外定标方法和装置。
背景技术
海水盐度是描述海洋状态的重要参数之一，研究其变化和分布规律对分析海洋自身特性以及了解海洋在海-气复杂系统中的作用有重要意义。并且，海水盐度也是气象学、生态学、水文学和渔业等其他学科与应用领域重点关注的研究对象。
在对海水盐度进行研究的过程中，首先需要获取海面盐度；目前，通过卫星搭载微波辐射计直接获取目标海面的亮温数据，然后，通过观测的亮温数据反演海水盐度。
但是，上述反演海面盐度值的方法中，观测亮温的精度对反演所得海水盐度有很大影响，由于盐度计的仪器性能衰减等原因会使得其观测值有一定的偏差，为了保证盐度计观测亮温的准确定及时纠正仪器观测亮温的偏差，需要定期对盐度计进行外定标以获得仪器定标系数用于纠正偏差。
发明内容
本发明的目的在于提供一种星载盐度计冷空外定标方法和装置，能够定期对盐度计进行外定标以获得仪器定标系数，用于纠正盐度计获取的观测亮温的偏差。
第一方面，本发明实施例提供了一种星载盐度计冷空外定标方法，包括：
根据星载盐度计的天线背瓣视场内的陆地比重，确定适于冷空外定标的预选区域；
根据星载盐度计的天线主瓣对应的预设观测数据，确定适于冷空外定标的时间窗口；观测数据包括：太空辐射、太阳辐射和月亮辐射；
根据预选区域和时间窗口对应的冷空外定标卫星实际运行轨道，计算标准天线温度；
对标准天线温度和星载盐度计的实测天线温度进行统计处理，得到用于校正盐度计观测亮温偏差的冷空外定标系数。
结合第一方面，本发明实施例提供了第一方面的第一种可能的实施方式，其中，根据星载盐度计的天线背瓣视场内的陆地比重，确定适于冷空外定标的预选区域，包括：
根据卫星轨道参数和卫星姿态参数，计算星载盐度计的天线指向的地面定位位置；其中，卫星轨道参数包括：卫星坐标矢量、卫星速度矢量；卫星姿态参数包括：卫星俯仰角和卫星滚动角；
根据预设网格化地球陆地掩膜数据确定地面定位位置对应的呈现状态；其中，呈现状态包括：陆地和海洋；
根据星载盐度计的天线对应的天线方向图对地面定位位置的陆地呈现状态进行积分计算，得到星载盐度计的天线背瓣视场内的陆地比重；
根据预设阈值对天线背瓣视场内的陆地比重进行判断，确定适于冷空外定标的预选区域。
结合第一方面的第一种可能的实施方式，本发明实施例提供了第一方面的第二种可能的实施方式，其中，根据卫星轨道参数和卫星姿态参数，计算星载盐度计的天线指向的地面定位位置，包括：
利用天线的观测视角θ和天线的方位角计算天线在卫星本体坐标系的单位观测矢量W：
按照预设旋转矩阵Tecr/sc＝Tecr/locTloc/sc，将卫星本体坐标系中的单位观测矢量W转换到地心旋转坐标系中，得到单位矢量D′；其中， 其中，Tloc/sc为单位观测矢量W从卫星本体坐标系到轨道坐标系的转换对应的三维旋转矩阵，三维旋转矩阵与卫星姿态有关；分别对应Zloc轴、Yloc轴和Xloc轴在地心旋转坐标系中的单位矢量；为地心直角坐标系中的卫星位置矢量；r为位置矢量的模；为卫星的速度矢量；v为速度的模；
按照以下几何关系计算目标观测点的坐标G；其中，几何关系包括：S+D＝G；D＝d·D′；S为卫星的位置矢量；D为天线观测矢量；G为观测目标点的位置矢量；
按照以下公式计算根据观测点坐标计算地理经纬度：
结合第一方面的第二种可能的实施方式，本发明实施例提供了第一方面的第三种可能的实施方式，其中，根据星载盐度计的天线对应的天线方向图对地面定位位置的陆地呈现状态进行积分计算，得到星载盐度计的天线背瓣视场内的陆地比重，包括：
按照以下公式对地面定位位置的陆地呈现状态进行积分计算，得到星载盐度计的天线背瓣视场内的陆地比重：其中，是天线旁瓣和背瓣增益矩阵，dΩ为像素对应的陆地、海水的立体角。
结合第一方面的第三种可能的实施方式，本发明实施例提供了第一方面的第四种可能的实施方式，其中，根据星载盐度计的天线主瓣对应的预设观测数据，确定适于冷空外定标的时间窗口，包括：
按照以下公式计算天线接收到的太空辐射：其中，G为天线方向图，Ω为太空在星载盐度计视向上的立体角，为星载盐度计指向上的太空辐射；
按照以下公式计算天线接收到的太阳辐射和月亮辐射：其中，为太阳或月亮的辐射亮温，ΩS/M为太阳或月亮对应的立体角；
根据太空辐射、太阳辐射和月亮辐射与天线温度稳定性的预设关联关系，确定适于冷空外定标的时间窗口。
结合第一方面的第四种可能的实施方式，本发明实施例提供了第一方面的第五种可能的实施方式，其中，根据预选区域和时间窗口对应的冷空外定标卫星实际运行轨道，计算标准天线温度，包括：
根据卫星实际运行轨道的轨道参数和时间，计算星载盐度计的天线实际接收到的实际观测数据；
根据星载盐度计的天线对应的天线方向图和实际观测数据，计算总的天线温度。
第二方面，本发明实施例还提供了一种星载盐度计冷空外定标装置，包括：
第一确定模块，用于根据星载盐度计的天线背瓣视场内的陆地比重，确定适于冷空外定标的预选区域；
第二确定模块，用于根据星载盐度计的天线主瓣对应的预设观测数据，确定适于冷空外定标的时间窗口；观测数据包括：太空辐射、太阳辐射和月亮辐射；
计算模块，用于根据预选区域和时间窗口对应的冷空外定标卫星实际运行轨道，计算标准天线温度；
统计处理模块，用于对标准天线温度和星载盐度计的实测天线温度进行统计处理，得到用于校正盐度计观测亮温偏差的冷空外定标系数。
结合第二方面，本发明实施例提供了第二方面的第一种可能的实施方式，其中，第一确定模块，包括：
第一计算单元，用于根据卫星轨道参数和卫星姿态参数，计算星载盐度计的天线指向的地面定位位置；其中，卫星轨道参数包括：卫星坐标矢量、卫星速度矢量；卫星姿态参数包括：卫星俯仰角和卫星滚动角；
第一确定单元，用于根据预设网格化地球陆地掩膜数据确定地面定位位置对应的呈现状态；其中，呈现状态包括：陆地和海洋；
积分计算单元，用于根据星载盐度计的天线对应的天线方向图对地面定位位置的陆地呈现状态进行积分计算，得到星载盐度计的天线背瓣视场内的陆地比重；
第二确定单元，用于根据预设阈值对天线背瓣视场内的陆地比重进行判断，确定适于冷空外定标的预选区域。
结合第二方面的第一种可能的实施方式，本发明实施例提供了第二方面的第二种可能的实施方式，其中，第一计算单元，包括：
第一计算子单元，用于利用天线的观测视角θ和天线的方位角计算天线在卫星本体坐标系的单位观测矢量W：
转换子单元，用于按照预设旋转矩阵Tecr/sc＝Tecr/locTloc/sc，将卫星本体坐标系中的单位观测矢量W转换到地心旋转坐标系中，得到单位矢量D′；其中， 其中，Tloc/sc为单位观测矢量W从卫星本体坐标系到轨道坐标系的转换对应的三维旋转矩阵，三维旋转矩阵与卫星姿态有关；分别对应Zloc轴、Yloc轴和Xloc轴在地心旋转坐标系中的单位矢量；为地心直角坐标系中的卫星位置矢量；r为位置矢量的模；为卫星的速度矢量；v为速度的模；
第二计算子单元，用于按照以下几何关系计算目标观测点的坐标G；其中，几何关系包括：S+D＝G；D＝d·D′；S为卫星的位置矢量；D为天线观测矢量；G为观测目标点的位置矢量；
第三计算子单元，用于按照以下公式计算根据观测点坐标计算地理经纬度：
结合第二方面的第二种可能的实施方式，本发明实施例提供了第二方面的第三种可能的实施方式，其中，积分计算单元，包括：
积分计算子单元，用于按照以下公式对地面定位位置的陆地呈现状态进行积分计算，得到星载盐度计的天线背瓣视场内的陆地比重：其中，是天线旁瓣和背瓣增益矩阵，dΩ为像素对应的陆地、海水的立体角；
设置子单元，用于将积分计算子单元的计算结果设置为星载盐度计的天线背瓣视场内的陆地比重。
本发明实施例提供的一种星载盐度计冷空外定标方法和装置，与现有技术中的由于仪器的衰减性能使得获取的亮温数据与实际值偏差较大相比，其能够通过计算星载盐度计翻转姿态下陆地比重对定标区域进行选择，通过模拟太空背景辐射、太阳辐射和月亮辐射对定标时间进行选择，最后利用模拟的天线温度对星载盐度计实际测得的天线温度进行统计处理，得到定标系数，这系统的解决了我国星载盐度计冷空定标的技术问题，并且根据该定标系数对盐度计获取的观测亮温进行纠正，使得纠正后的观测亮温更准确。
为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1示出了本发明实施例所提供的一种星载盐度计冷空外定标方法的流程图；
图2示出了本发明实施例所提供的另一种星载盐度计冷空外定标方法的流程图；
图3示出了本发明实施例所提供的星载盐度计地理定位原理示意图；
图4示出了本发明实施例所提供的另一种星载盐度计冷空外定标方法的流程图；
图5示出了本发明实施例所提供的一种星载盐度计冷空外定标装置的结构示意图；
图6示出了本发明实施例所提供的一种星载盐度计冷空外定标装置中第一确定模块和第一计算单元的结构示意图；
图7示出了本发明实施例所提供的一种星载盐度计冷空外定标装置中积分计算单元和第二确定模块的结构示意图；
图8示出了本发明实施例所提供的一种星载盐度计冷空外定标装置中计算模块的结构示意图。
具体实施方式
下面将结合本发明实施例中附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
现有技术中直接通过观测的亮温数据反演海水盐度，由于仪器的衰减性能使得获取的亮温数据与实际值偏差较大，故使根据上述亮温数据反演得到的海水盐度与实际盐度值偏离较大。基于该问题，本发明实施例提供的一种星载盐度计冷空外定标方法和装置，与现有技术中的由于仪器的衰减性能使得获取的亮温数据与实际值偏差较大相比，定期对盐度计进行外定标以获得仪器定标系数，用于纠正盐度计获取的观测亮温的偏差。
参考图1，本发明实施例提供了一种星载盐度计冷空外定标方法，所述方法包括如下步骤：
S101、根据星载盐度计的天线背瓣视场内的陆地比重，确定适于冷空外定标的预选区域。
在上述步骤中，盐度计即传感器，该传感器上设置有天线，而星载盐度计就是在卫星上搭载盐度计；星载盐度计在正常运行时，是天线的主瓣指向地球，天线的背瓣指向太空；而本发明实施例中，是计算星载盐度计翻转状态下的天线背瓣视场内的陆地比重；星载盐度计的翻转状态下，即天线的主瓣指向太空，天线的背瓣指向地球，计算该翻转状态下天线背瓣视场内的陆地比重；并根据陆地比重筛选出适合冷空定标的预选区域。
S102、根据所述星载盐度计的天线主瓣对应的预设观测数据，确定适于冷空外定标的时间窗口；所述观测数据包括：太空辐射、太阳辐射和月亮辐射。
具体的，星载盐度计的翻转状态下，即天线的主瓣指向太空，天线的背瓣指向地球，本发明实施例通过在预测轨道上模拟太空辐射、太阳辐射和月亮辐射等预设观测数据，并根据模拟的上述预设观测数据筛选出适合冷空定标的时间窗口；其中，上述预测轨道即根据卫星的当前实际运行轨道及卫星的运行参数，预测的卫星在下一个时刻的运行轨道，在确定了上述预设轨道且计算出上述每个预测轨道上对应的预设观测数据后，根据上述预设观测数据对天线温度稳定性的影响(即上述预设观测数据与天线温度稳定性的关联关系)，确定适合冷空定标的时间窗口；其中，该时间窗口优选为一个时间段。
S103、根据所述预选区域和所述时间窗口对应的冷空外定标卫星实际运行轨道，计算标准天线温度。
步骤101中确定了适于冷空外定标的预选区域，步骤102中确定了适于冷空外定标的时间窗口，根据上述预选区域和时间窗口即确定了冷空外定标对应的卫星实际运行轨道；而在该卫星实际运行轨道中，可以根据卫星实际运行轨道的轨道参数、卫星运行时的姿态数据、太空辐射数据、太阳和月亮的同步辐射数据等综合计算出总的天线温度。
S104、对所述标准天线温度和所述星载盐度计的实测天线温度进行统计处理，得到用于校正所述盐度计观测亮温偏差的冷空外定标系数。
具体的，首先获取在步骤103确定的冷空外定标对应的卫星实际运行轨道中，获取星载盐度计直接采集的实测天线温度；将该实测天线温度与步骤103中综合计算出总的天线温度进行统计处理，即可得到冷空外定标系数，该冷空外定标系数用于校正所述盐度计的观测亮温偏差，使得校正后的观测亮温更接近于实际的观测亮温，以便后续根据该观测亮温反演的海水盐度值更接近于实际盐度值。
本发明实施例提供的一种星载盐度计冷空外定标方法和装置，与现有技术中的由于仪器的衰减性能使得获取的亮温数据与实际值偏差较大相比，其能够通过计算星载盐度计翻转姿态下陆地比重对定标区域进行选择，通过模拟太空背景辐射、太阳辐射和月亮辐射对定标时间进行选择，最后利用模拟的天线温度对星载盐度计实际测得的天线温度进行统计处理，得到定标系数，这系统的解决了我国星载盐度计冷空定标的技术问题，并且根据该定标系数对盐度计获取的观测亮温进行纠正，使得纠正后的观测亮温更准确。
进一步，参考图2，本发明实施例中进行针对上述步骤101中，确定适于冷空外定标的预选区域的具体方法包括如下步骤：
S201、根据卫星轨道参数和卫星姿态参数，计算所述星载盐度计的天线指向的地面定位位置；其中，所述卫星轨道参数包括：卫星坐标矢量、卫星速度矢量；所述卫星姿态参数包括：卫星俯仰角和卫星滚动角。
具体的，星载盐度计冷空定标区域选择方法是首先利用卫星轨道参数，包括卫星坐标矢量、卫星速度矢量；卫星姿态参数，包括卫星俯仰角、卫星滚动角，然后根据上述卫星轨道参数和卫星姿态参数计算传感器指向的地面定位位置(或者为定位点)。
S202、根据预设网格化地球陆地掩膜数据确定所述地面定位位置对应的呈现状态；其中，所述呈现状态包括：陆地和海洋。
具体的，预设网格化地球陆地掩膜数据是标准数据，根据该网格化地球陆地掩膜数据可以判断上述地面定位位置的呈现状态，即确定是陆地还是海洋。
S203、根据所述星载盐度计的天线对应的天线方向图对所述地面定位位置的陆地呈现状态进行积分计算，得到所述星载盐度计的天线背瓣视场内的陆地比重。
具体的，在确定了上述地面定位位置的呈现状态后，根据天线方向图计算传感器背瓣视场内的陆地比重。
S204、根据预设阈值对所述天线背瓣视场内的陆地比重进行判断，确定适于冷空外定标的预选区域。
具体的，设置一个陆地比重阈值，根据该陆地比重阈值对陆地比重进行判断，确定陆地比重小于上述陆地比重阈值对应的区域为冷空外定标的预选区域。
具体的，在设置陆地比重阈值时，背景场造成的天线温度误差要小于0.02K。在假设陆地和海洋最高辐射差为100K的前提下，统计得到适于进行冷空定标的阈值在0.02％-0.03％，范围对应不同的天线波束和极化。
具体的，上述步骤201中，根据卫星轨道参数和卫星姿态参数，计算所述星载盐度计的天线指向的地面定位位置，包括：
利用天线的观测视角θ和天线的方位角计算天线在卫星本体坐标系的单位观测矢量W：
按照预设旋转矩阵Tecr/sc＝Tecr/locTloc/sc，将卫星本体坐标系中的单位观测矢量W转换到地心旋转坐标系中，得到单位矢量D′；其中， 其中，Tloc/sc为单位观测矢量W从卫星本体坐标系到轨道坐标系的转换对应的三维旋转矩阵，所述三维旋转矩阵与卫星姿态有关；分别对应Zloc轴、Yloc轴和Xloc轴在地心旋转坐标系中的单位矢量；为地心直角坐标系中的卫星位置矢量；r为所述位置矢量的模；为卫星的速度矢量；v为速度的模；
按照以下几何关系计算目标观测点的坐标G；其中，所述几何关系包括：S+D＝G；D＝d·D′；S为卫星的位置矢量；D为天线观测矢量；G为观测目标点的位置矢量；
按照以下公式计算根据所述观测点坐标计算地理经纬度：
具体的，上述步骤201的过程，先计算天线在卫星本体坐标系的单位观测矢量W，公式如下：然后将星体坐标系(即上述卫星本体坐标系)中的单位观测矢量W转换到地心旋转坐标系中，得到单位矢量D′；其旋转矩阵为Tecr/sc＝Tecr/locTloc/sc。
星体坐标系到轨道坐标系的转换是由与卫星姿态有关的三维旋转矩阵实现的。卫星在运动过程中的俯仰(pitch)、滚动(roll)和偏航(yaw)三个角度是随时间变化的，因此该坐标系的转换矩阵也是也是随时间变化的。转换矩阵为：
星体坐标系到地心旋转坐标系之间的变换，首先要计算得到星体坐标系的三个坐标轴在地心旋转坐标系中的单位矢量。
Zloc轴、Yloc轴和Xloc轴在地心旋转坐标系中的单位矢量分别为：因此，局部坐标系到地心旋转坐标系的矩阵为：
其中，为地心直角坐标系中的卫星位置矢量，r为位置矢量的模，为卫星的速度矢量，v为速度的模。
至此已经求出地心旋转坐标系下的卫星位置矢量S和天线单位观测矢量D′，如果再求得观测矢量的长度d,就可以得到观测矢量D，那么根据图3所示的几何关系就可以解算出目标点的位置G。
如图3所示，根据卫星的位置矢量S、天线观测矢量D和观测目标点的位置矢量G之间的几何关系可以得出：S+D＝G，即：S+d·D′＝G；
式中存在两个未知数，必须再建立一个方程，先求出斜距d，才能得到最终的结果。展开得到：
建立椭球方程，将待求点G代入：
中，其中，a为地球的半长轴，e为地球偏心率。于是得到：
将上式展开整理可得，A·d2+B·d+C＝0；
其中，
解一元二次方程，取斜距的两个根中较小的那一个，因为这个根是观测矢量与地球近端的交点，求得d，可以计算得到目标观测点的坐标G。
考虑地球扁率的影响，通过观测点坐标计算地理经纬度如下式：
通过以上计算步骤就可以唯一确定天线视角为θ、天线方位角为处的目标观测点经纬度坐标。
具体的，上述计算得到的目标观测点经纬度坐标即确定适于冷空外定标的预选区域。
在星载盐度计冷空定标区域选择中还传感器背瓣陆地比重计算两项关键技术，本发明实施例具体通过以下方法实现，所述方法具体包括：
按照以下公式对所述地面定位位置的陆地呈现状态进行积分计算，得到所述星载盐度计的天线背瓣视场内的陆地比重：其中，是天线旁瓣和背瓣增益矩阵，dΩ为像素对应的陆地、海水的立体角。
具体的，首先，需要计算卫星平台翻转时天线背瓣视场内陆地的比重，这是为了保证视场内目标辐射相对均一且容易模拟。在实际计算中，假设地球目标辐射值相同，则陆地目标在视场中的比重为：
其中，是天线旁瓣和背瓣增益矩阵，dΩ为像素对应的陆地、海水的立体角。
由于在卫星升、降轨过程中背瓣视场不同，需要对卫星每个重复周期内的所有升降轨都进行计算。
在确定了适于冷空外定标的预选区域之后，还需要确定适于冷空外定标的时间窗口，具体的，如图3所示，确定时间窗口的具体方法具体如下：
星载盐度计冷空定标时间选择是首先利用卫星轨道参数，包括卫星坐标矢量、速度矢量，将处于卫星本体坐标系传感器指向转换到天球坐标系，然后根据天球坐标系下的传感器指向获得太空辐射数据、太阳辐射数据和月亮辐射数据，即根据天线方向图计算得到天线接收到的太空辐射、太阳辐射和月亮辐射，通过评估太空辐射对的天线温度的稳定性的影响以及太阳辐射和月亮辐射对天线温度的影响筛选出适合冷空定标的时间。
对应的，利用天线方向图计算天线接收到的太空辐射如下：
按照以下公式计算天线接收到的太空辐射：其中，G为天线方向图，Ω为太空在星载盐度计视向上的立体角，为星载盐度计指向上的太空辐射；
由于太阳和月亮可以看作点光源，所以太阳/月亮的直接辐射可以采用以下公式计算，即：
按照以下公式计算天线接收到的太阳辐射和月亮辐射：其中，为太阳或月亮的辐射亮温，ΩS/M为太阳或月亮对应的立体角；
根据所述太空辐射、所述太阳辐射和所述月亮辐射与所述天线温度稳定性的预设关联关系，确定适于冷空外定标的时间窗口。
由于盐度计卫星采用太阳同步轨道，卫星平台相对太阳的方向是不变的，但由于地球运行轨道的原因，卫星平台相对于冷空背景的角度变化大约为1°/天，变化周期为1年。针对同一定标轨道，在不同时间进行冷空定标时冷空背景辐射并不一致，需要在一年中选择合适的时间以满足冷空定标中冷空辐射场均匀的要求。为了实现精确定标，要求在定标过程中，及选定的定标时间内天线温度的变化不能超过0.1K。
具体的，根据确定的适于冷空外定标的预选区域和适于冷空外定标的时间窗口，即可确定用于冷空外定标的卫星实际运行轨道，参考图4，根据该用于冷空外定标的卫星实际运行轨道，计算标准天线温度的具体方法包括：
S301、根据所述卫星实际运行轨道的轨道参数和时间，计算星载盐度计的天线实际接收到的实际观测数据。
S302、根据星载盐度计的天线对应的天线方向图和所述实际观测数据，计算总的天线温度。
综合上述步骤301～上述步骤302，基于星载盐度计冷空定标区域选择和时间选择的确定星载盐度计进行冷空定标的轨道。在星载盐度计完成冷空定标的操作后，基于用于冷空外定标的卫星实际运行轨道的轨道参数和时间，计算传感器天线实际接收到的太空辐射、太阳辐射和月亮辐射，并根据天线方向图计算得到总的天线温度；以便于后续将技术的总的天线温度与实测的天线温度对比进行定标，最终获得冷空定标系数。
本发明实施例提供的一种星载盐度计冷空外定标方法，与现有技术中的由于仪器的衰减性能使得获取的亮温数据与实际值偏差较大相比，其能够通过计算星载盐度计翻转姿态下陆地比重对定标区域进行选择，通过模拟太空背景辐射、太阳辐射和月亮辐射对定标时间进行选择，最后利用模拟的天线温度对星载盐度计实际测得的天线温度进行统计处理，得到定标系数，这系统的解决了我国星载盐度计冷空定标的技术问题，并且根据该定标系数对盐度计获取的观测亮温进行纠正，使得纠正后的观测亮温更准确。
本发明实施例还提供了一种星载盐度计冷空外定标装置，所述装置用于执行上述星载盐度计冷空外定标方法，参考图5，所述装置包括：
第一确定模块11，用于根据星载盐度计的天线背瓣视场内的陆地比重，确定适于冷空外定标的预选区域；
第二确定模块12，用于根据星载盐度计的天线主瓣对应的预设观测数据，确定适于冷空外定标的时间窗口；观测数据包括：太空辐射、太阳辐射和月亮辐射；
计算模块13，用于根据预选区域和时间窗口对应的冷空外定标卫星实际运行轨道，计算标准天线温度；
统计处理模块14，用于对标准天线温度和星载盐度计的实测天线温度进行统计处理，得到用于校正盐度计观测亮温偏差的冷空外定标系数。
本发明实施例提供的一种星载盐度计冷空外定标装置，与现有技术中的由于仪器的衰减性能使得获取的亮温数据与实际值偏差较大相比，其能够通过计算星载盐度计翻转姿态下陆地比重对定标区域进行选择，通过模拟太空背景辐射、太阳辐射和月亮辐射对定标时间进行选择，最后利用模拟的天线温度对星载盐度计实际测得的天线温度进行统计处理，得到定标系数，这系统的解决了我国星载盐度计冷空定标的技术问题，并且根据该定标系数对盐度计获取的观测亮温进行纠正，使得纠正后的观测亮温更准确。
进一步的，该星载盐度计冷空外定标装置中，参考图6，第一确定模块11，包括：
第一计算单元111，用于根据卫星轨道参数和卫星姿态参数，计算星载盐度计的天线指向的地面定位位置；其中，卫星轨道参数包括：卫星坐标矢量、卫星速度矢量；卫星姿态参数包括：卫星俯仰角和卫星滚动角；
第一确定单元112，用于根据预设网格化地球陆地掩膜数据确定地面定位位置对应的呈现状态；其中，呈现状态包括：陆地和海洋；
积分计算单元113，用于根据星载盐度计的天线对应的天线方向图对地面定位位置的陆地呈现状态进行积分计算，得到星载盐度计的天线背瓣视场内的陆地比重；
第二确定单元114，用于根据预设阈值对天线背瓣视场内的陆地比重进行判断，确定适于冷空外定标的预选区域。
进一步的，该星载盐度计冷空外定标装置中，参考图6，第一计算单元111，包括：
第一计算子单元1111，用于利用天线的观测视角θ和天线的方位角计算天线在卫星本体坐标系的单位观测矢量W：
转换子单元1112，用于按照预设旋转矩阵Tecr/sc＝Tecr/locTloc/sc，将卫星本体坐标系中的单位观测矢量W转换到地心旋转坐标系中，得到单位矢量D′；其中， 其中，Tloc/sc为单位观测矢量W从卫星本体坐标系到轨道坐标系的转换对应的三维旋转矩阵，三维旋转矩阵与卫星姿态有关；分别对应Zloc轴、Yloc轴和Xloc轴在地心旋转坐标系中的单位矢量；为地心直角坐标系中的卫星位置矢量；r为位置矢量的模；为卫星的速度矢量；v为速度的模；
第二计算子单元1113，用于按照以下几何关系计算目标观测点的坐标G；其中，几何关系包括：S+D＝G；D＝d·D′；S为卫星的位置矢量；D为天线观测矢量；G为观测目标点的位置矢量；
第三计算子单元1114，用于按照以下公式计算根据观测点坐标计算地理经纬度：
进一步的，该星载盐度计冷空外定标装置中，参考图7，积分计算单元113，包括：
积分计算子单元1131，用于按照以下公式对地面定位位置的陆地呈现状态进行积分计算，得到星载盐度计的天线背瓣视场内的陆地比重：其中，是天线旁瓣和背瓣增益矩阵，dΩ为像素对应的陆地、海水的立体角；
设置子单元1132，用于将积分计算子单元的计算结果设置为星载盐度计的天线背瓣视场内的陆地比重。
进一步的，该星载盐度计冷空外定标装置中，参考图7，第二确定模块12包括：
第二计算单元121，用于按照以下公式计算天线接收到的太空辐射：其中，G为天线方向图，Ω为太空在星载盐度计视向上的立体角，为星载盐度计指向上的太空辐射；
第三计算单元122，用于按照以下公式计算天线接收到的太阳辐射和月亮辐射：其中，为太阳或月亮的辐射亮温，ΩS/M为太阳或月亮对应的立体角；
第三确定单元123，用于根据太空辐射、太阳辐射和月亮辐射与天线温度稳定性的预设关联关系，确定适于冷空外定标的时间窗口。
进一步的，该星载盐度计冷空外定标装置中，参考图8，计算模块13，包括：
第四计算单元131，用于根据卫星实际运行轨道的轨道参数和时间，计算星载盐度计的天线实际接收到的实际观测数据；
第五计算单元132，用于根据星载盐度计的天线对应的天线方向图和实际观测数据，计算总的天线温度。
本发明实施例提供的一种星载盐度计冷空外定标装置，与现有技术中的由于仪器的衰减性能使得获取的亮温数据与实际值偏差较大相比，其能够通过计算星载盐度计翻转姿态下陆地比重对定标区域进行选择，通过模拟太空背景辐射、太阳辐射和月亮辐射对定标时间进行选择，最后利用模拟的天线温度对星载盐度计实际测得的天线温度进行统计处理，得到定标系数，这系统的解决了我国星载盐度计冷空定标的技术问题，并且根据该定标系数对盐度计获取的观测亮温进行纠正，使得纠正后的观测亮温更准确。
本发明实施例所提供的进行星载盐度计冷空外定标的计算机程序产品，包括存储了程序代码的计算机可读存储介质，所述程序代码包括的指令可用于执行前面方法实施例中所述的方法，具体实现可参见方法实施例，在此不再赘述。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
在本申请所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。
Claims (10)
1.一种星载盐度计冷空外定标方法，其特征在于，包括：
根据星载盐度计的天线背瓣视场内的陆地比重，确定适于冷空外定标的预选区域；
根据所述星载盐度计的天线主瓣对应的预设观测数据，确定适于冷空外定标的时间窗口；所述观测数据包括：太空辐射、太阳辐射和月亮辐射；
根据所述预选区域和所述时间窗口对应的冷空外定标卫星实际运行轨道，计算标准天线温度；
对所述标准天线温度和所述星载盐度计的实测天线温度进行统计处理，得到用于校正所述盐度计观测亮温偏差的冷空外定标系数。
2.根据权利要求1所述的星载盐度计冷空外定标方法，其特征在于，所述根据星载盐度计的天线背瓣视场内的陆地比重，确定适于冷空外定标的预选区域，包括：
根据卫星轨道参数和卫星姿态参数，计算所述星载盐度计的天线指向的地面定位位置；其中，所述卫星轨道参数包括：卫星坐标矢量、卫星速度矢量；所述卫星姿态参数包括：卫星俯仰角和卫星滚动角；
根据预设网格化地球陆地掩膜数据确定所述地面定位位置对应的呈现状态；其中，所述呈现状态包括：陆地和海洋；
根据所述星载盐度计的天线对应的天线方向图对所述地面定位位置的陆地呈现状态进行积分计算，得到所述星载盐度计的天线背瓣视场内的陆地比重；
根据预设阈值对所述天线背瓣视场内的陆地比重进行判断，确定适于冷空外定标的预选区域。
3.根据权利要求2所述的星载盐度计冷空外定标方法，其特征在于，所述根据卫星轨道参数和卫星姿态参数，计算所述星载盐度计的天线指向的地面定位位置，包括：
利用天线的观测视角θ和天线的方位角计算天线在卫星本体坐标系的单位观测矢量W：
按照预设旋转矩阵Tecr/sc＝Tecr/locTloc/sc，将卫星本体坐标系中的单位观测矢量W转换到地心旋转坐标系中，得到单位矢量D′；其中，
其中，Tloc/sc为单位观测矢量W从卫星本体坐标系到轨道坐标系的转换对应的三维旋转矩阵，所述三维旋转矩阵与卫星姿态有关；分别对应Zloc轴、Yloc轴和Xloc轴在地心旋转坐标系中的单位矢量；为地心直角坐标系中的卫星位置矢量；r为所述位置矢量的模；为卫星的速度矢量；v为速度的模；
按照以下几何关系计算目标观测点的坐标G；其中，所述几何关系包括：S+D＝G；D＝d·D′；S为卫星的位置矢量；D为天线观测矢量；G为观测目标点的位置矢量；
按照以下公式计算根据所述观测点坐标计算地理经纬度：
4.根据权利要求3所述的星载盐度计冷空外定标方法，其特征在于，所述根据所述星载盐度计的天线对应的天线方向图对所述地面定位位置的陆地呈现状态进行积分计算，得到所述星载盐度计的天线背瓣视场内的陆地比重，包括：
按照以下公式对所述地面定位位置的陆地呈现状态进行积分计算，得到所述星载盐度计的天线背瓣视场内的陆地比重：其中，是天线旁瓣和背瓣增益矩阵，dΩ为像素对应的陆地、海水的立体角。
5.根据权利要求4所述的星载盐度计冷空外定标方法，其特征在于，所述根据所述星载盐度计的天线主瓣对应的预设观测数据，确定适于冷空外定标的时间窗口，包括：
按照以下公式计算天线接收到的太空辐射：其中，G为天线方向图，Ω为太空在星载盐度计视向上的立体角，为星载盐度计指向上的太空辐射；
按照以下公式计算天线接收到的太阳辐射和月亮辐射：其中，为太阳或月亮的辐射亮温，ΩS/M为太阳或月亮对应的立体角；
根据所述太空辐射、所述太阳辐射和所述月亮辐射与所述天线温度稳定性的预设关联关系，确定适于冷空外定标的时间窗口。
6.根据权利要求5所述的星载盐度计冷空外定标方法，其特征在于，所述根据所述预选区域和所述时间窗口对应的冷空外定标卫星实际运行轨道，计算标准天线温度，包括：
根据所述卫星实际运行轨道的轨道参数和时间，计算星载盐度计的天线实际接收到的实际观测数据；
根据星载盐度计的天线对应的天线方向图和所述实际观测数据，计算总的天线温度。
7.一种星载盐度计冷空外定标装置，其特征在于，包括：
第一确定模块，用于根据星载盐度计的天线背瓣视场内的陆地比重，确定适于冷空外定标的预选区域；
第二确定模块，用于根据所述星载盐度计的天线主瓣对应的预设观测数据，确定适于冷空外定标的时间窗口；所述观测数据包括：太空辐射、太阳辐射和月亮辐射；
计算模块，用于根据所述预选区域和所述时间窗口对应的冷空外定标卫星实际运行轨道，计算标准天线温度；
统计处理模块，用于对所述标准天线温度和所述星载盐度计的实测天线温度进行统计处理，得到用于校正所述盐度计观测亮温偏差的冷空外定标系数。
8.根据权利要求7所述的星载盐度计冷空外定标装置，其特征在于，所述第一确定模块，包括：
第一计算单元，用于根据卫星轨道参数和卫星姿态参数，计算所述星载盐度计的天线指向的地面定位位置；其中，所述卫星轨道参数包括：卫星坐标矢量、卫星速度矢量；所述卫星姿态参数包括：卫星俯仰角和卫星滚动角；
第一确定单元，用于根据预设网格化地球陆地掩膜数据确定所述地面定位位置对应的呈现状态；其中，所述呈现状态包括：陆地和海洋；
积分计算单元，用于根据所述星载盐度计的天线对应的天线方向图对所述地面定位位置的陆地呈现状态进行积分计算，得到所述星载盐度计的天线背瓣视场内的陆地比重；
第二确定单元，用于根据预设阈值对所述天线背瓣视场内的陆地比重进行判断，确定适于冷空外定标的预选区域。
9.根据权利要求8所述的星载盐度计冷空外定标装置，其特征在于，所述第一计算单元，包括：
第一计算子单元，用于利用天线的观测视角θ和天线的方位角计算天线在卫星本体坐标系的单位观测矢量W：
转换子单元，用于按照预设旋转矩阵Tecr/sc＝Tecr/locTloc/sc，将卫星本体坐标系中的单位观测矢量W转换到地心旋转坐标系中，得到单位矢量D′；其中，
其中，Tloc/sc为单位观测矢量W从卫星本体坐标系到轨道坐标系的转换对应的三维旋转矩阵，所述三维旋转矩阵与卫星姿态有关；分别对应Zloc轴、Yloc轴和Xloc轴在地心旋转坐标系中的单位矢量；为地心直角坐标系中的卫星位置矢量；r为所述位置矢量的模；为卫星的速度矢量；v为速度的模；
第二计算子单元，用于按照以下几何关系计算目标观测点的坐标G；其中，所述几何关系包括：S+D＝G；D＝d·D′；S为卫星的位置矢量；D为天线观测矢量；G为观测目标点的位置矢量；
第三计算子单元，用于按照以下公式计算根据所述观测点坐标计算地理经纬度：
10.根据权利要求9所述的星载盐度计冷空外定标装置，其特征在于，所述积分计算单元，包括：
积分计算子单元，用于按照以下公式对所述地面定位位置的陆地呈现状态进行积分计算，得到所述星载盐度计的天线背瓣视场内的陆地比重：其中，是天线旁瓣和背瓣增益矩阵，dΩ为像素对应的陆地、海水的立体角；
设置子单元，用于将积分计算子单元的计算结果设置为所述星载盐度计的天线背瓣视场内的陆地比重。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610158460.6A CN105842259B (zh) | 2016-03-18 | 2016-03-18 | 星载盐度计冷空外定标方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610158460.6A CN105842259B (zh) | 2016-03-18 | 2016-03-18 | 星载盐度计冷空外定标方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105842259A CN105842259A (zh) | 2016-08-10 |
| CN105842259B true CN105842259B (zh) | 2018-07-24 |
Family
ID=56588293
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610158460.6A Expired - Fee Related CN105842259B (zh) | 2016-03-18 | 2016-03-18 | 星载盐度计冷空外定标方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105842259B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110764153B (zh) * | 2019-11-15 | 2021-06-08 | 上海航天测控通信研究所 | 一种星载微波成像仪热镜背瓣在轨误差校正系统及方法 |
| CN112485738B (zh) * | 2020-11-12 | 2022-09-20 | 国家卫星气象中心(国家空间天气监测预警中心) | 检验静止轨道微波辐射计系统稳定性的方法、系统和装置 |
| CN113341238B (zh) * | 2021-06-08 | 2022-12-16 | 北方天穹信息技术（西安）有限公司 | 一种利用太阳辐射测量天线方向图的方法 |
| CN117191197B (zh) * | 2023-09-08 | 2024-03-19 | 国家卫星海洋应用中心 | 一种基于均匀面目标的观测亮温定标方法、装置及设备 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010063844A1 (fr) * | 2008-12-05 | 2010-06-10 | Thales | Procede de geo-localisation d'un objet par multitelemetrie |
| CN104155246A (zh) * | 2014-08-26 | 2014-11-19 | 中国海洋大学 | 海水盐度的检测装置与方法 |
| CN105067120A (zh) * | 2015-07-28 | 2015-11-18 | 国家卫星海洋应用中心 | 星载微波辐射计观测亮温的动态滤波重采样方法及装置 |
-
2016
- 2016-03-18 CN CN201610158460.6A patent/CN105842259B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010063844A1 (fr) * | 2008-12-05 | 2010-06-10 | Thales | Procede de geo-localisation d'un objet par multitelemetrie |
| CN104155246A (zh) * | 2014-08-26 | 2014-11-19 | 中国海洋大学 | 海水盐度的检测装置与方法 |
| CN105067120A (zh) * | 2015-07-28 | 2015-11-18 | 国家卫星海洋应用中心 | 星载微波辐射计观测亮温的动态滤波重采样方法及装置 |
Non-Patent Citations (3)
| Title |
|---|
| 基于微波遥感的海表面盐度探测机制研究进展;严明 等;《遥感信息》;20150430;第30卷(第2期);第17-25页 * |
| 海洋二号扫描微波辐射计冷空定标和地球物理参数反演研究;周武 等;《中国工程科学》;20131231;第15卷(第7期);第75-80页 * |
| 神舟四号飞船微波辐射计定标和检验(Ⅰ)—微波辐射计外定标;王振占 等;《遥感学报》;20040930;第 8卷(第 5 期);第397-403页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105842259A (zh) | 2016-08-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105842259B (zh) | 星载盐度计冷空外定标方法和装置 | |
| Reid et al. | The trigonometric parallax of Cygnus X-1 | |
| Martin et al. | Wind‐wave‐induced velocity in ATI SAR ocean surface currents: First experimental evidence from an airborne campaign | |
| US10962625B2 (en) | Celestial positioning system and method | |
| CN108225307B (zh) | 一种惯性测量信息辅助的星图匹配方法 | |
| CN105067120B (zh) | 星载微波辐射计观测亮温的动态滤波重采样方法及装置 | |
| CN107870043A (zh) | 一种海表参数同步反演优化方法 | |
| DE102011076602A1 (de) | GNSS - Atmosphärische Abschätzung mit einem Federated-Ionospärischen Filter | |
| US8837782B1 (en) | Geoposition determination using satellite ephemerides | |
| CN107238824A (zh) | 基于先验dem数据的星载sar图像几何精校正方法 | |
| CN105737806B (zh) | 基于星载微波辐射计反演海面盐度的方法和装置 | |
| CN104459728B (zh) | 一种基于gnss定位的磁偏角校准方法 | |
| CN108562305B (zh) | 一种惯性/天文深组合导航系统安装误差五位置快速粗标定方法 | |
| Tang et al. | The in-orbit calibration method based on terrain matching with pyramid-search for the spaceborne laser altimeter | |
| Perez et al. | Evolution of CIR storm on 22 July 2009 | |
| Hosokawa et al. | Estimating drift velocity of polar cap patches with all‐sky airglow imager at Resolute Bay, Canada | |
| CN106872962B (zh) | 一种用于星载激光测高仪检校的地面探测器的布设方法 | |
| Verma et al. | Expected precision of Europa Clipper gravity measurements | |
| Rivier et al. | European CO 2 fluxes from atmospheric inversions using regional and global transport models | |
| CN105444778A (zh) | 一种基于成像几何反演的星敏感器在轨定姿误差获取方法 | |
| Chisham | A new methodology for the development of high‐latitude ionospheric climatologies and empirical models | |
| Krall et al. | The effect of the thermosphere on quiet time plasmasphere morphology | |
| CN111652931A (zh) | 一种地理定位方法、装置、设备及计算机可读存储介质 | |
| Bolandi et al. | Design of a supervisory adaptive attitude control (SAAC) system for a stereo-imagery satellite based on multiple model control with switching | |
| CN107526066B (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: 20180724 Termination date: 20210318 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |