一种卫星遥感数据的处理方法及相关装置

本申请涉及遥感数据处理领域，尤其涉及一种卫星遥感数据的处理方法及相关装置。

背景技术：

近年来，随着各类遥感对地观测卫星的成功发射，有关地球及其各种资源环境的多分辨率的、海量的、实时的对地观测数据被不断获取。随着遥感技术的发展，其数据规模还在以惊人的速度增长。然而，由于受卫星成像方式的影响，遥感卫星在获取数据时，相邻影像间存在不同程度的航向和旁向重叠，并且其旁向重叠范围从赤道向两极逐渐增大。

目前，卫星数据的处理往往以每景影像的图幅大小作为存储与管理的最小单元。并且，采用逐景计算的方式进行卫星遥感数据的处理，其中，对于相邻两景影像的重叠区域，得到两个计算结果，需要从两个计算结果中选择一个计算结果作为最终计算结果，因此，存在如何选取重叠区域的计算结果的问题，即在计算结果的取值上存在歧义性问题。

另外，在对相邻两景影像分别进行计算的过程中，对于相邻两景影像的重叠区域，在对任一景影像的计算过程中，只有该景影像的重叠区域的数据参与运算，另一景影像的重叠区域的数据不参与计算，从而，导致重叠区域的数据利用率低的问题。

技术实现要素：

本申请提供了一种卫星遥感数据的处理方法及相关装置，目的在于解决重叠区域数据利用率低以及重叠区域的计算结果存在歧义性的问题。

为了实现上述目的，本申请提供了以下技术方案：

本申请提供了一种卫星遥感数据的处理方法，包括：

获取待研究区域各个轨道采集的各景卫星遥感图像；

对所述卫星遥感图像转投影至正弦等面积投影，得到投影后的卫星遥感图像；

按照目标度数的网格分别对投影后的卫星遥感图像进行切割，得到瓦片集合；所述目标度数小于预设度数；

依据所述瓦片集合中瓦片的经纬度，对所述瓦片进行存储。

可选的，所述依据所述瓦片集合中瓦片的经纬度，对所述瓦片进行存储，包括：

按照所述瓦片集合中瓦片的经纬度进行排列，得到第一区域；

以预设的第一尺寸为分割单元，对所述第一区域进行分割，得到第一分割单元集合；

以预设的第二尺寸为分割单元，分别对所述第一分割单元集合中的各个第一分割单元进行分割，得到第二分割单元集合；所述第二尺寸小于所述第一尺寸；

确定所述瓦片集合中每个瓦片分别在所述第一分割单元集合中所属的第一分割单元，以及在所述第二分割单元集合中所属的第二分割单元；

对所述瓦片集合中的各个瓦片以及各个瓦片的信息进行存储；其中，所述瓦片集合中任一瓦片的信息包括：该瓦片所属第一分割单元的经纬度、该瓦片所属第二分割单元经纬度和该瓦片的经纬度。

可选的，在所述获取待研究区域各个轨道采集的各景卫星遥感图像之后，还包括：

检测所述各景卫星遥感图像中的无效像元；

生成所述无效像元的掩膜，得到各景卫星遥感图像分别对应的掩膜图像；

将所述掩膜图像转投影至正弦等面积投影，得到投影后的掩膜图像；

按照所述目标度数的网格分别对投影后的掩膜图像进行切割，得到所述瓦片集合中各个瓦片分别对应的掩膜瓦片；

所述依据所述瓦片集合中瓦片的经纬度，对所述瓦片进行存储，具体为：

对所述瓦片集合中的目标瓦片进行存储；所述目标瓦片为所述瓦片集合中对应的掩膜瓦片中无效数据占比不大于预设阈值的瓦片。

可选的，在所述获取待研究区域各个轨道采集的各景卫星遥感图像之后，并且，在所述对所述卫星遥感图像转投影至正弦等面积投影，得到投影后的卫星遥感图像之前，还包括：

在所述卫星遥感图像为压缩文件的情况下，对所述压缩文件进行解压缩，得到解压缩后的卫星遥感图像；

分别对每景解压缩后的卫星遥感图像进行波段合成，得到每景波段合成后的遥感卫星图像；

分别对每景波段合成后的卫星遥感图像进行几何配准和几何校正，得到处理后的各景卫星遥感图像。

可选的，所述目标度数为0.05°。

可选的，所述瓦片集合中任一瓦片的信息还包括：采集该瓦片数据的传感器类型、该瓦片的数据类型、分级层次和成像时间。

本申请还提供了一种卫星遥感数据的处理装置，包括：

获取模块，用于获取待研究区域各个轨道采集的各景卫星遥感图像；

投影模块，用于对所述卫星遥感图像转投影至正弦等面积投影，得到投影后的卫星遥感图像；

切割模块，用于按照目标度数的网格分别对投影后的卫星遥感图像进行切割，得到瓦片集合；所述目标度数小于预设度数；

存储模块，用于依据所述瓦片集合中瓦片的经纬度，对所述瓦片进行存储。

可选的，所述存储模块，用于依据所述瓦片集合中瓦片的经纬度，对所述瓦片进行存储，包括：

所述存储模块，具体用于按照所述瓦片集合中瓦片的经纬度进行排列，得到第一区域；

以预设的第一尺寸为分割单元，对所述第一区域进行分割，得到第一分割单元集合；

以预设的第二尺寸为分割单元，分别对所述第一分割单元集合中的各个第一分割单元进行分割，得到第二分割单元集合；所述第二尺寸小于所述第一尺寸；

确定所述瓦片集合中每个瓦片分别在所述第一分割单元集合中所属的第一分割单元，以及在所述第二分割单元集合中所属的第二分割单元；

对所述瓦片集合中的各个瓦片以及各个瓦片的信息进行存储；其中，所述瓦片集合中任一瓦片的信息包括：该瓦片所属第一分割单元的经纬度、该瓦片所属第二分割单元经纬度和该瓦片的经纬度。

本申请还提供了一种存储介质，所述存储介质包括存储的程序，其中，所述程序执行上述任意一项所述的卫星遥感数据的处理方法。

本申请还提供了一种设备，所述设备包括至少一个处理器、以及与所述处理器连接的至少一个存储器、总线；其中，所述处理器、所述存储器通过所述总线完成相互间的通信；所述处理器用于调用所述存储器中的程序指令，以执行如上述任一项所述的卫星遥感数据的处理方法。

本申请所述的卫星遥感数据的处理方法及装置，获取待研究区域各个轨道采集的各景卫星遥感图像；对卫星遥感图像转投影至正弦等面积投影，由于该投影相较于经纬度、横轴墨卡托投影等在经纬度的变形是最小的，因此，得到的投影后的卫星遥感图像的质量受到最小的影响；按照目标度数的网格分别对投影后的卫星遥感图像进行切割，得到瓦片集合，依据瓦片集合中瓦片的经纬度，对瓦片进行存储。因此，本申请是以瓦片为最小存储单元进行存储。其中，由于目标度数小于预设度数，因此，切割得到的瓦片的尺寸较小。

由于本申请存储的瓦片的尺寸较小，并且，以瓦片为计算单位，一方面，使得不同经纬度的瓦片之间不会重叠，因此，对预设经纬度下的瓦片进行处理时，不存在计算结果存在歧义性的问题。另一方面，对预设经纬度下的瓦片进行处理时，是对不同景卫星遥感图像中该经纬度下的全部瓦片数据都进行了处理，所以，对于重叠区域，提高了数据利用率。

综上所述，本申请可以解决现有技术中，计算结果存在歧义性的问题，以及重叠区域数据利用率低的问题。

附图说明

为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本申请实施例公开的一种卫星遥感数据的处理方法的流程图；

图2为本申请实施例公开的又一种卫星遥感数据的处理方法的流程图；

图3为本申请实施例公开的一种三级陆地格网管理模式示意图；

图4为本申请实施例公开的现有数据存储方式与本申请实施例的三级陆地格网管理模式下的各地理位置下的时相频次分布示意图；

图5为本申请实施例公开的一种卫星遥感数据的处理装置的结构示意图；

图6为本申请实施例公开的一种设备的结构示意图。

具体实施方式

下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。

发明人在研究中发现，现有卫星遥感应用多集中于区域尺度，并且对卫星数据的处理也往往以图幅大小为处理单元而忽视相邻轨道卫星影像的重叠问题。若相邻两轨卫星数据的成像时间一致，那么逐景影像分别处理然后根据地理位置进行拼接的思路是合理的。然而，对于绝大多数卫星遥感数据而言，由于受到云、阴影等因素的影响，要保证相邻两景数据时相一致且完全清晰的前提基本无法满足。因此，若采用传统逐景处理的方式，那么重叠区域因相邻两轨数据的光谱差异而导致在最终分析取值上存在较大的歧义性。例如：zhangetal.(2018)在生产2013年中国区域30米地表覆盖产品时，以landsat逐图幅为处理单元而忽视了相邻轨道landsat的光谱差异，导致最终生产的中国区域30米产品在空间维度存在较为明显的斑块现象。

另外，发明人在研究中还发现，重叠区域因两轨数据同时覆盖，实际上具备更加丰富的信息量，若分别对每景卫星遥感图像进行数据分析和处理实际上是降低了卫星数据的使用效率，尤其是高纬度区域。以landsat卫星数据为例，kovalskyyandroy(2013)定量分析了相邻两轨卫星数据的重叠度随纬度的关系，结果表明在赤道区域其轨道重叠度约为0.196°但在南北纬80°的高纬度区域对应的轨道重叠度高达12.535°和12.826°，即高纬度区域landsat卫星数据的轨道重叠区域能够达到90％以上。

基于此，本申请实施例提供了一种卫星遥感数据的处理方法，包括以下步骤：

s101、获取待研究区域各个轨道采集的各景卫星遥感图像。

在本实施例中，各景卫星遥感图像的获取方式为现有技术，这里不再赘述。

在本步骤中，获取到的各景卫星遥感图像可能是经过压缩的卫星遥感图像。

s102、对获取的各景卫星遥感图像进行预处理。

在本实施例中，预处理过程可以包括以下步骤a1～a3：

a1、为压缩文件的情况下，对压缩文件进行解压缩，得到解压缩后的卫星遥感图像。

a2、分别对每景解压缩后的卫星遥感图像进行波段合成，得到每景波段合成后的遥感卫星图像。

a3、分别对每景波段合成后的卫星遥感图像进行几何配准和几何校正，得到处理后的各景卫星遥感图像。

上述步骤a1～步骤a3的具体实现方式都为现有技术，这里不再赘述。

需要说明的是，在实际中，本步骤是可选步骤。只是执行本步骤后，得到的预处理后的卫星遥感图像的质量，相较于不执行本步骤获取的卫星遥感图像的质量高。

在实际中，如果执行本步骤，以下步骤中的卫星遥感图像为预处理后得到的卫星遥感图像。如果不执行本步骤，以下步骤中的卫星遥感图像为解压缩后的卫星遥感图像。

s103、检测各景卫星遥感图像中的无效像元。

在本实施例中，由于无效像元可以包括：云、阴影和背景值等。

在本步骤中，可以利用云和阴影检测算法，检测各景卫星遥感图像中的无效像元，本实施例不对具体的检测算法作限定。

s104、生成无效像元的掩膜，得到各景卫星遥感图像分别对应的掩膜图像。

在本实施例中，通过将各景卫星遥感图像中，对无效像元所在区域生成掩膜，得到各景卫星遥感图像分别对应的掩膜图像。

s105、对卫星遥感图像和掩膜图像分别转投影至正弦等面积投影，得到投影后的卫星遥感图像和投影后的掩膜图像。

在本实施例中，由于正弦等面积投影能够保证全球任意区域具有相同的采样密度(该投影相较于经纬度、横轴墨卡托投影等在经纬度的变形也是最小的)，因此，在本步骤中，将卫星遥感图像和掩膜图像分别转投影至正弦等面积投影，使得得到的投影后的卫星遥感图像和投影后的掩膜图像的质量得到保证。

s106、按照目标度数的网格分别对投影后的卫星遥感图像和投影后的掩膜图像进行切割，得到瓦片集合和掩膜瓦片集合。

在本实施例中，目标度数小于预设度数，例如，目标度数可以为0.05°，当然，在实际中，目标度数也可以为其他度数。通过实验证明，在实际中，目标度数为0.05°时，对瓦片的存储效率可以达到较高。

在本步骤中，切割的具体实现方式为现有技术，这里不再赘述。

s107、确定瓦片集合中的目标瓦片。

在本步骤中，目标瓦片指：瓦片集合中对应的掩膜瓦片中无效数据占比不大于预设阈值的瓦片。

针对无效像元(云、阴影和背景值)的比例(预设阈值)选择，研究在中国区域的多时相应用中选取了30％(若研究区卫星数据较为稀疏可以适当提高该阈值)。本实施例不对预设阈值的取值作限定。

s108、依据目标瓦片的经纬度，对目标瓦片进行存储。

本实施例具有以下有益效果：

在本实施例中，通过对各景卫星遥感图像中的无效像元生成掩膜，得到各景卫星遥感图像对应的掩膜图像，对卫星遥感图像和掩膜图像分别转投影至正弦等面积投影，得到投影后的卫星遥感图像和投影后的掩膜图像，并对投影后的卫星遥感图像和投影后的掩膜图像按照目标度数进行分割，得到瓦片集合和瓦片集合中各个瓦片对应的掩膜瓦片，其中，掩膜瓦片中无效像元占比越高，表示该掩膜瓦片对应的瓦片的价值越小，因此，本实施例对瓦片集合中的目标瓦片进行存储，即不对非目标瓦片进行存储，从而，节省了存储空间。以landsat数据为例，本实施例能够节省近1/3的存储空间。

图2为本申请实施例提供的又一种卫星遥感数据的处理方法，可以包括以下步骤：

s201、获取待研究区域各个轨道采集的各景卫星遥感图像。

本步骤的含义以及具体实现方式可以参考s101，这里不再赘述。

s202、对获取的卫星遥感图像转投影至正弦等面积投影，得到投影后的卫星遥感图像。

本步骤的含义以及具体实现方式可以参考s105，这里不再赘述。

s203、按照目标度数的网格分别对投影后的卫星遥感图像进行切割，得到瓦片集合。

本步骤的含义以及具体实现方式可以参考s106，这里不再赘述。

s204、依据瓦片集合中瓦片的经纬度，对瓦片进行存储。

在本步骤中，对瓦片集合中的瓦片进行存储。

在实际中，由于分割得到的瓦片集合中瓦片总量将达到一个庞大的数量级，因此，为了提升瓦片数据的管理效率，研究提出了一种三级格网管理模式datacube。

具体的，本步骤的实现可以包括以下步骤b1～步骤b5：

b1、按照瓦片集合中瓦片的经纬度进行排列，得到第一区域。

在本步骤中，按照瓦片集合中各个瓦片的经纬度，可以将瓦片按照经纬度进行排列，从而，形成一个区域。为了描述方便，本实施例将形成的区域称为第一区域。

b2、以预设的第一尺寸为分割单元，对第一区域进行分割，得到第一分割单元集合。

在本实施例中，第一尺寸的取值可以为10°×10°，可以将全球划分为36×18个瓦片块。当然，在实际中，第一尺寸的取值可以根据实际需求确定，本实施例不对第一尺寸的具体取值作限定。

在本实施例中，为了描述方便，将对第一区域进行分割得到的各个分割区域称为第一分割单元，则分割得到的全部第一分割单元构成第一分割单元集合。

在实际中，如果第一尺寸的取值为10°×10°，目标度数为0.05°，则第一分割单包含200×200个0.05°×0.05°瓦片，因此，第一分割单元包含的瓦片数据量依旧较大。因此，本实施例对每个第一分割单元进一步分割，以更加适用于区域尺度应用。具体的，可以参考步骤b3。

b3、以预设的第二尺寸为分割单元，分别对第一分割单元集合中的各个第一分割单元进行分割，得到第二分割单元集合。

在实施例中，第二尺寸小于第一尺寸。其中，第二尺寸的取值可以为1.43°×1.43°约为150km×150km，将第一分割单元分割为7×7的单元，为了描述方便，将本步骤分割得到的每个单元，称为第二分割单元，将由分割得到的全部第二分割单元构成的集合称为第二分割单元集合。

b4、确定瓦片集合中每个瓦片分别在第一分割单元集合中所属的第一分割单元，以及在第二分割单元集合中所属的第二分割单元。

在本实施例中，由于第一区域存在经纬度，因此，对第一区域进行分割得到的第一分割单元集合中的各个第一分割单元分别存在对应的经纬度，进一步，对每个第一分割单元进行分割，得到的第二分割单元集合中的各个第二分割单元都分别对应有经纬度。

因此，在本步骤中，可以依据瓦片集合中各个瓦片的经纬度，分别确定各个瓦片在第一分割单元集合中所属的第一分割单元，以及在第二分割单元集合中所属的第二分割单元。

b5、对瓦片集合中的各个瓦片以及各个瓦片的信息进行存储。

在本步骤中，瓦片集合中任一瓦片的信息包括：该瓦片所属第一分割单元的经纬度、该瓦片所属第二分割单元经纬度和该瓦片的经纬度。

可选的，在实际中，瓦片集合中任一瓦片的信息还可以包括：采集该瓦片数据的传感器类型、该瓦片的数据类型、分级层次和成像时间。

可选的，在本实施例中，瓦片集合中每个瓦片的命名规则可以为“l8sr_t3_hh**vv**_h*v*_p**r**_yyyymmdd.gif”，其中‘l8sr’指明传感器类型和数据类型，‘t3’指代数据分级层次(例如，瓦片集合中任一瓦片的信息包括：该瓦片所属第一分割单元的经纬度、该瓦片所属第二分割单元经纬度和该瓦片的经纬度，保存了三级信息，因此，分级层次是三级)，‘hh**vv**’指代该基础瓦片在第一级格网下的位置信息(以目标度数为0.05°为例，则瓦片集合中每个瓦片的大小为0.05°×0.05°，将全球划分为36×18瓦片块，hh，vv对应的两位数编号取值范围分别为[0,36)和[0,18))。‘h*v*’指代该瓦片在第二级格网下的位置信息(其中h和v的取值范围为[0,6])；‘p**r**’指代瓦片在第三级格网中的位置信息(其中p和r的取值范围为[0,29])；最后，‘yyyymmdd’标记了瓦片数据的成像时间信息。

在本实施例中，由于对瓦片集合中的每个瓦片进行存储时，每个瓦片存储的信息包括：该瓦片所属的第一分割单元的经纬度、所属的第二分割单元的经纬度，以及该瓦片的经纬度等三级位置信息。其中，该瓦片所属的不同层级的单元的大小不同，从而，在实际中，可以更加适用于区域尺度应用，即根据待处理区域，从适用的层级上确定待处理区域包括哪些经纬度的瓦片，从而提高对瓦片的管理效率和应用效率。

为了直观展示本实施例提供的三级陆地格网管理模式示意图，给出了图3所示的示意图。图3中，标注“第三层级”的小格指示的是一个瓦片，在本实施例中，该瓦片的经纬度表示第三层级信息，标注“第二层级”的小格指示的瓦片块的经纬度表示该瓦片的第二层级信息，标注“第一层级”的小格指示的瓦片块的经纬度表示该瓦片的第一层级信息。

为了显示本申请的效果，图4中(a)为landsat现有数据存储方式下各地理位置下的时相频次分布，图4中(b)为经过本申请实施例的datacube三级陆地格网管理模式下的时相频次分布。其中，由于卫星数据是周期性观测(如：每隔16天重复观测一次，故在一年内同一轨道图幅上会有多频次的卫星数据，因此我们统计时将各景位置上的频次称作时相频次分布)。

通过图4中(a)和(b)的对比，可以看出，对于相同经纬度区域(相同地理位置的一格)，(a)比(b)的时相频次高，即本申请实施例以瓦片为单位进行处理，可以提高重叠区域的数据利用率。

图5为本申请实施例提供的一种卫星遥感数据的处理装置，可以包括：获取模块501、投影模块502、切割模块503和存储模块504，其中，

获取模块501，用于获取待研究区域各个轨道采集的各景卫星遥感图像；

投影模块502，用于对所述卫星遥感图像转投影至正弦等面积投影，得到投影后的卫星遥感图像；

切割模块503，用于按照目标度数的网格分别对投影后的卫星遥感图像进行切割，得到瓦片集合；所述目标度数小于预设度数；

存储模块504，用于依据所述瓦片集合中瓦片的经纬度，对所述瓦片进行存储。

可选的，所述存储模块504，用于依据所述瓦片集合中瓦片的经纬度，对所述瓦片进行存储，包括：

所述存储模块504，具体用于按照所述瓦片集合中瓦片的经纬度进行排列，得到第一区域；

以预设的第一尺寸为分割单元，对所述第一区域进行分割，得到第一分割单元集合；

以预设的第二尺寸为分割单元，分别对所述第一分割单元集合中的各个第一分割单元进行分割，得到第二分割单元集合；所述第二尺寸小于所述第一尺寸；

确定所述瓦片集合中每个瓦片分别在所述第一分割单元集合中所属的第一分割单元，以及在所述第二分割单元集合中所属的第二分割单元；

对所述瓦片集合中的各个瓦片以及各个瓦片的信息进行存储；其中，所述瓦片集合中任一瓦片的信息包括：该瓦片所属第一分割单元的经纬度、该瓦片所属第二分割单元经纬度和该瓦片的经纬度。

可选的，该装置还可以包括：执行模块，用于在所述获取模块501获取待研究区域各个轨道采集的各景卫星遥感图像之后，检测所述各景卫星遥感图像中的无效像元；生成所述无效像元的掩膜，得到各景卫星遥感图像分别对应的掩膜图像；将所述掩膜图像转投影至正弦等面积投影，得到投影后的掩膜图像；按照所述目标度数的网格分别对投影后的掩膜图像进行切割，得到所述瓦片集合中各个瓦片分别对应的掩膜瓦片；

所述存储模块504，用于依据所述瓦片集合中瓦片的经纬度，对所述瓦片进行存储，具体为：

所述存储模块504，具体用于对所述瓦片集合中的目标瓦片进行存储；所述目标瓦片为所述瓦片集合中对应的掩膜瓦片中无效数据占比不大于预设阈值的瓦片。

可选的，该装置还可以包括：预处理模块，用于在所述获取模块501获取待研究区域各个轨道采集的各景卫星遥感图像之后，并且，在所述投影模块502对所述卫星遥感图像转投影至正弦等面积投影，得到投影后的卫星遥感图像之前，在所述卫星遥感图像为压缩文件的情况下，对所述压缩文件进行解压缩，得到解压缩后的卫星遥感图像；分别对每景解压缩后的卫星遥感图像进行波段合成，得到每景波段合成后的遥感卫星图像；分别对每景波段合成后的卫星遥感图像进行几何配准和几何校正，得到处理后的各景卫星遥感图像。

可选的，所述目标度数为0.05°。

可选的，所述瓦片集合中任一瓦片的信息还包括：采集该瓦片数据的传感器类型、该瓦片的数据类型、分级层次和成像时间。

卫星遥感数据的处理装置包括处理器和存储器，上述获取模块501、投影模块502、切割模块503和存储模块504等均作为程序单元存储在存储器中，由处理器执行存储在存储器中的上述程序单元来实现相应的功能。

处理器中包含内核，由内核去存储器中调取相应的程序单元。内核可以设置一个或以上，通过调整内核参数来解决重叠区域数据利用率低以及重叠区域的计算结果存在歧义性的问题。

本发明实施例提供了一种存储介质，其上存储有程序，该程序被处理器执行时实现所述卫星遥感数据的处理方法。

本发明实施例提供了一种处理器，所述处理器用于运行程序，其中，所述程序运行时执行所述卫星遥感数据的处理方法。

本发明实施例提供了一种设备，如图6所示，设备包括至少一个处理器、以及与处理器连接的至少一个存储器、总线；其中，处理器、存储器通过总线完成相互间的通信；处理器用于调用存储器中的程序指令，以执行上述的卫星遥感数据的处理方法。本文中的设备可以是服务器、pc、pad、手机等。

本申请还提供了一种计算机程序产品，当在数据处理设备上执行时，适于执行初始化有如下方法步骤的程序：

获取待研究区域各个轨道采集的各景卫星遥感图像；

对所述卫星遥感图像转投影至正弦等面积投影，得到投影后的卫星遥感图像；

按照目标度数的网格分别对投影后的卫星遥感图像进行切割，得到瓦片集合；所述目标度数小于预设度数；

依据所述瓦片集合中瓦片的经纬度，对所述瓦片进行存储。

本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

在一个典型的配置中，设备包括一个或多个处理器(cpu)、存储器和总线。设备还可以包括输入/输出接口、网络接口等。

存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器(ram)和/或非易失性内存等形式，如只读存储器(rom)或闪存(flashram)，存储器包括至少一个存储芯片。存储器是计算机可读介质的示例。

计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(pram)、静态随机存取存储器(sram)、动态随机存取存储器(dram)、其他类型的随机存取存储器(ram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、快闪记忆体或其他内存技术、只读光盘只读存储器(cd-rom)、数字多功能光盘(dvd)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitorymedia)，如调制的数据信号和载波。

还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。

本领域技术人员应明白，本申请的实施例可提供为方法、系统或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。

以上仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。

本申请实施例方法所述的功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算设备可读取存储介质中。基于这样的理解，本申请实施例对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该软件产品存储在一个存储介质中，包括若干指令用以使得一台计算设备(可以是个人计算机，服务器，移动计算设备或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-onlymemory)、随机存取存储器(ram，randomaccessmemory)、磁碟或者光盘等各种可以存储程序代码的介质。

本说明书的各个实施例中记载的特征可以相互替换或者组合，每个实施例重点说明的都是与其它实施例的不同之处，各个实施例之间相同或相似部分互相参见即可。

对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本申请。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本申请的精神或范围的情况下，在其它实施例中实现。因此，本申请将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。