本发明属于遥感卫星轨道预报领域，涉及一种预测海量卫星观测需求满足度的方法，尤其涉及一种利用卫星轨道回归特性预测海量用户年度观测需求满足度的方法。
背景技术：
：我国遥感卫星大多采用太阳同步回归轨道特性，数据已广泛应用于国土、林业、水利、环境保护、农业等多个领域，发挥着重要的作用。随着用户对遥感数据应用的不断深入，用户观测需求数量不断增加，陆地观测卫星主要用户在每年年初提出一年的观测需求，主要用户年度观测需求达到上万个，需求呈现出用户多、数量多、区域大、频度高，时间跨度长等特点。经过卫星任务规划、指令生成上注、云量判读和质量检测后，从最终拍摄的数据以及有效性可以计算得到用户需求实际满足度。在用户提交年度观测需求后，数据获取前期，用户更加关注每个需求的预测满足度和有效覆盖率，以便对未来数据获取情况进行预判。目前针对卫星点目标和有限区域目标观测需求任务规划的国内外模型和算法较多，为了满足用户不同需求，有多种观测策略，针对少量观测需求的短期满足度，可以使用专用任务规划软件或者商业stk软件(satellitetoolskit)，结合卫星重访特性和回归特性，通过人机交互计算得出。但是用户的年度观测需求区域包括大区域和小区域以及国内外区域，数量大、并且时间跨度为一年，观测频度高，针对数以万计的海量大区域时间跨度长的观测需求，没有形成统一的观测需求格式和自动化快速有效的算法，靠人工参与统计，需要通过分析用户需求分布、卫星轨道特性、观测频度、天气预测计算，计算工作量大，工作效率低下，成本高。技术实现要素：本发明的目的在于克服现有技术的上述不足，提供一种预测海量卫星观测需求满足度的方法，该方法通过规定用户年度观测需求格式，设计一种自动读取用户观测需求，按照卫星轨道特性和卫星常规运行特点对全球进行区域划分，与海量区域和不同时间段需求自动匹配计算满足度，并给出指定区域有效满足度计算方法，为用户观测需求未来获取的满足度和有效性提供依据。本发明的上述目的主要是通过如下技术方案予以实现的：一种预测海量卫星观测需求满足度的方法，包括如下步骤：步骤(1)、制定统一的用户年度观测需求格式，按照所述观测需求格式填写年度观测需求，形成用户观测需求数据库；步骤(2)、计算用户观测需求数据库的每个观测需求中完成每一次观测从开始到结束的日期跨度；步骤(3)、对全球划分区域，针对卫星的每个传感器建立划分的各区域的覆盖率表，所述覆盖率表包括卫星覆盖周期内每个划分区域的覆盖率；步骤(4)、将单个需求区域与全球划分区域取交集，得到每个交集区域占所述单个需求区域面积的百分比；根据所述面积的百分比计算单个需求区域中单个传感器的满足度；步骤(5)、重复步骤(4)，得到所有需求区域中每个需求区域的单个传感器的满足度；步骤(6)、对步骤(5)中得到的所有需求区域中单个传感器的满足度取均值，得到所述单个传感器的满足度。在上述预测海量卫星观测需求满足度的方法中，所述用户年度观测需求格式包括：用户单位和具体观测信息，所述具体观测信息包括序号、观测需求名称、观测日期、观测次数、优先级、观测区域和观测数据类型。在上述预测海量卫星观测需求满足度的方法中，所述步骤(3)中不受任何约束情况下卫星覆盖周期为m天，针对卫星的每个传感器建立划分的各区域的覆盖率表中包括每个划分区域的覆盖率qi，i＝1…k,k表示划分区域的个数。在上述预测海量卫星观测需求满足度的方法中，所述步骤(4)中根据面积的百分比计算单个需求区域中单个传感器的满足度的具体公式如下；其中：b表示单个需求区域中单个传感器的满足度；c表示完成每一次观测从开始到结束的日期跨度；m表示为卫星覆盖周期；qi表示卫星覆盖周期内每个划分区域的覆盖率；pi表示每个交集区域占单个需求区域面积的百分比；n表示交集区域的个数。在上述预测海量卫星观测需求满足度的方法中，所述步骤(5)中所有需求区域中每个需求区域的单个传感器的满足度通过如下公式计算得到；其中：bj表示第j个需求区域中单个传感器的满足度；c表示完成每一次观测从开始到结束的日期跨度；m表示卫星覆盖周期；qi表示卫星覆盖周期内每个划分区域的覆盖率；pi表示每个交集区域占单个需求区域面积的百分比；n表示交集区域的个数。在上述预测海量卫星观测需求满足度的方法中，所述步骤(6)中对所有需求区域中单个传感器的满足度取均值，得到所述单个传感器的满足度的具体公式如下：其中：bs表示单个传感器的满足度；n表示需求区域的个数；bj表示第j个需求区域中单个传感器的满足度。在上述预测海量卫星观测需求满足度的方法中，单个需求区域中单个传感器的有效满足度通过如下公式计算得到：s＝(1-(1-pab)t)*100％t＝dd/m其中：pab表示a区域b季节历史数据有效比例；t表示用户需求日期跨度内可覆盖的遍数；dd表示观测从开始到结束的日期跨度；m表示卫星覆盖周期；s表示有效满足度。在上述预测海量卫星观测需求满足度的方法中，重复步骤(1)-(6)，得到用户需求的所有传感器的满足度。本发明与现有技术相比具有如下有益效果：(1)、本发明通过规定用户年度观测需求格式，设计一种自动读取用户观测需求，按照卫星轨道特性和卫星常规运行特点对全球进行区域划分，与海量区域和不同时间段需求自动匹配计算满足度，并给出指定区域有效满足度计算方法，为用户观测需求未来获取的满足度和有效性提供依据。(2)、本发明制定了用户卫星观测需求格式并建立了观测需求库，为海量用户需求自动读取、按照时间段、观测频度和卫星需求进行分解、满足度自动化计算奠定了基础，大大提高了工作效率和计算准确度，克服了人工方式效率低、耗时长的缺陷。(3)、本发明充分利用了卫星轨道回归特性，根据覆盖周期和使用约束对全球观测区域进行分割，克服使用瞬时轨道根数外推造成的精度差及计算时间效率问题，将全球分区与用户区域自动匹配，快速计算海量观测需求满足度，克服了人工计算工作量大的难题，且计算结果为预测遥感用户需求满足度，未来获取可能性提供依据。(4)、本发明方法可以同时获得需求观测满足度和需求观测有效满足度，满足了不同的用户需求。附图说明图1为本发明预测海量卫星观测需求满足度的方法流程框图。具体实施方式下面结合附图和具体实施例对本发明作进一步详细的描述：如图1所示为本发明预测海量卫星观测需求满足度的方法流程框图，本发明的具体实现方式包括如下步骤：一、确定用户观测需求申请格式。针对海量用户观测需求，定义统一的用户观测需求表格式，各用户按照统一格式自行填写年度观测需求。用户观测需求表，内容包括：用户单位和具体观测信息，其中具体观测信息包括序号、观测需求名称、观测日期、观测次数、优先级、观测区域、观测数据类型。其中观测日期表示需求的有效期，填写格式为“x月x日-x月x日”；观测次数为观测日期范围内需拍摄次数，填写数字；观测优先级由用户单位自评，分为高、中、低三档可以选择；观测区域填写区域坐标(按角点顺序填写经纬度)；观测数据类型为卫星载荷组合。下表1为用户观测需求表格式中的具体观测信息。表1用户观测需求表(具体观测信息)二、海量观测需求数据库建立。按照用户观测需求表各个输入单元，建立需求数据库，包括用户单位、需求序号、需求名称、开始日期、结束日期、观测次数、优先级、观测区域角点坐标、观测数据类型。计算用户观测需求数据库的每个观测需求中完成每一次观测从开始到结束的日期跨度；将每一次观测从开始到结束的日期跨度输入需求数据库，得到表2所示的需求数据库表结构，其中开始、结束日期格式为y-m-d，如2017年1月1日表示为2017-1-1。将收集到的各用户观测需求导入到需求数据库表中，集中管理，并用于计算输入。其中观测1次从开始到结束的日期跨度具体计算步骤如下：(1)首先计算开始结束日期日积数x，如下式：式中：m——月份数；d——日期数。(2)计算开始结束日期跨度dd，如下式：dd＝x2-x1式中：x2——结束日期日积数；x1——开始日期日积数。(3)观测1次日期跨度c，如下式：c＝[dd/p](p≠1)式中：p——观测次数。表2需求数据库表结构三、按覆盖周期划分全球区域。首先，确定不受任何约束情况下卫星覆盖周期。针对太阳同步回归轨道遥感卫星，在不受任何使用约束条件下，m天可实现星下点成像对全球(海量用户观测需求)覆盖一遍，即覆盖周期为m天，m值根据卫星传感器幅宽确定。例如，卫星回归周期41天，回归总圈次605圈，若传感器的幅宽为70公里，则覆盖周期m为41天，若传感器的幅宽为800公里，则覆盖周期m为4天，若传感器的幅宽为35公里，则覆盖周期m为82天。其次，根据使用约束将全球按照覆盖周期不同覆盖率划分区域。受卫星一天运行圈次、一圈运行时间、卫星固存容量、地面站布设、传感器幅宽等限制，以及用户需求分布特点，同一幅宽传感器在全球不同区域覆盖率不同。根据以上特点，针对每个传感器，将全球按照m天不同覆盖率划分为n类个区域：第一类区域为一个覆盖周期m天可完整覆盖区域，m天覆盖率为100％；第二类区域受各种使用约束，需两个覆盖周期交替成像可完整覆盖区域，m天覆盖率为50％；第i类区域(0≤i≤n)，m天覆盖率为qi；第n类区域为不能成像区域，m天覆盖率为0。建立卫星传感器m天全球区域划分覆盖率表，见表3。表3全球区域划分及覆盖率卫星传感器区域划分类别经纬度范围覆盖率某卫星a相机一类(经度,纬度)……1某卫星a相机二类(经度,纬度)……0.5…………………………某卫星a相机i类(经度,纬度)……qi…………………………某卫星a相机n类(经度,纬度)……0四、单个需求区域中单个传感器满足度计算。用户提出的年度观测需求，包含多个需求，每个需求中的观测数据类型为多种，计算单个需求单传感器满足度，具体计算步骤如下：(1)单个需求区域划分。将单个需求区域与步骤3中的全球划分的k个区域取交集，得到n个交集小区域(即重叠区域)，n小于或等于k，即单个需求区域与划分的k个区域均有重叠，此时n＝k，否则n＜k。使用gis软件计算出单个需求区域面积s和交集小区域的面积si(0≤i≤n)。计算每个小区域占需求区域的百分比pi，如下式：pi＝si/s式中：s——需求区域面积；si——每个交集区域面积。(2)计算单个需求区域中单个传感器满足度b，如下式：式中：其中：b表示单个需求区域中单个传感器的满足度；c表示完成每一次观测从开始到结束的日期跨度；m表示为卫星覆盖周期；qi表示卫星覆盖周期内每个划分区域的覆盖率；pi表示每个交集区域占单个需求区域面积的百分比；n表示交集区域的个数。若计算结果中b＞1，则令b＝1。五、按照步骤四的方法，计算所有需求区域中每个需求区域的单个传感器的满足度，具体公式如下；其中：bj表示第j个需求区域中单个传感器的满足度；c表示完成每一次观测从开始到结束的日期跨度；m表示为卫星覆盖周期；qi表示卫星覆盖周期内每个划分区域的覆盖率；pi表示每个交集区域占单个需求区域面积的百分比；n表示交集区域的个数。若计算结果中bj＞1，则令bj＝1。六、用户满足度计算。(1)单用户单个传感器满足度计算。将单个用户需求数据库表中观测数据类型按照传感器分类进行解析，针对筛选出卫星单个传感器的观测需求，按照步骤五计算该传感器每个需求区域的满足度bj。将该传感器的所有观测需求区域满足度取均值，即为该传感器满足度bs，如下式。其中：bs表示单个传感器的满足度；n表示需求区域的个数；bj表示第j个需求区域中单个传感器的满足度。(2)用户满足度计算。重复上述步骤，计算单个用户所有需求单颗卫星每个传感器的满足度，取均值，即为用户的满足度。七、中国区域有效满足度计算。该步骤给出用户需求在中国区域且分布在给定的同一区域同一季节计算有效满足度的方法。按照国家地理划分的7大区域：华东、华北、华中、华南、西南、西北和东北(可以按照不同区域划分)，利用卫星历史数据统计每个区域云量情况，给出每个季节各区域的有效比例，结合各区域历史数据有效比例和用户需求分布特点，单个需求单个传感器有效满足度按如下公式计算：s＝(1-(1-pab)t)*100％t＝dd/m式中：pab表示a区域b季节历史数据有效比例；t表示用户需求日期跨度内可覆盖的遍数。dd表示观测从开始到结束的日期跨度；m表示卫星覆盖周期。s表示有效满足度。以上所述，仅为本发明最佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本
技术领域：
的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。当前第1页12