CN111460060B - 一种用户感兴趣区域遥感影像空间索引方法 - Google Patents

Abstract

本发明提供了一种用户感兴趣区域遥感影像空间索引方法，属于遥感影像空间索引领域。该方法通过扩充遥感影像空间范围，将遥感影像各个像素位置信息嵌入到Google S2某一层级的索引空间，实现遥感影像像素索引与Google S2空间索的引对接。利用Google S2快速的空间索引与空间覆盖算法，快速获取用户感兴趣区域遥感影像像素位置信息，并生成用户感兴趣区域遥感影像。该方法有助于摒弃传统遥感影像服务平台以“景”为单位的影像下载方法，可以按用户感兴趣区域提供所需即所得的下载服务。本发明有效减少了数据网络传输量与传输时间，降低了用户后续数据处理难度。

Description

一种用户感兴趣区域遥感影像空间索引方法

技术领域

本发明属于遥感影像空间索引领域，更具体地，涉及一种用户感兴趣区域遥感影像空间索引方法。

背景技术

随着卫星遥感技术的发展与成熟，卫星遥感数据在人类生活中起着越来越重要的作用，在地理国情普查、农业监测管理、城市规划设计、环境监测保护、经济风险评估、自然灾害监测等涉及到国计民生的诸多领域扮演者重要角色。遥感为人类社会发展带来了巨大利益，但同时自身的局限也制约着其在人类生活中的应用。其中，卫星遥感数据的获取效率是重要制约因素之一。

当前各卫星遥感数据平台通常提供“捆绑式”数据服务，即将卫星遥感数据以“景”为单位提供数据下载服务，不可避免地提供了大量冗余数据，而用户所需数据只占实际获取数据相当小的比例。极端情况如，用户只需若干像素或者较小地理空间范围的遥感数据，却不得不下载一整“景”遥感影像，而其所需数据只占实际获取数据的万分之一或者更少。大量冗余数据耗费了更多的网络资源和数据传输时间，同时也增加了后续数据处理的难度，如需要更多的数据处理时间、更大的存储资源以及更强的计算资源。由此可见，遥感数据用户有相当大一部分时间都消耗在遥感数据的获取与处理上。

因此，改变这种“捆绑式”数据获取方式，为用户提供兴趣区域所得即所需的数据服务，具有巨大的科研价值与经济效益。

发明内容

针对现有技术的以上缺陷或改进需求，本发明提供了一种用户感兴趣区域遥感影像空间索引方法，其目的是提供一种用户感兴趣区域卫星遥感数据快速索引优化方法，根据用户感兴趣区域生成卫星遥感影像数据，以解决“捆绑式”数据获取方式导致的大量冗余数据，造成网络资源和数据传输时间的不必要耗费，同时增加了后续数据处理难度的技术问题。

为实现上述目的，本发明提供了一种用户感兴趣区域遥感影像空间索引方法，包括：

S1.根据用户感兴趣区域，获取基准投影带和邻接投影带；

所述基准投影带表示属于同一投影带且与用户感兴趣区域重叠面积最大的遥感影像集所属的投影带；所述邻接投影带表示与用户感兴趣区域相交的遥感影像中，除去基准投影带遥感影像后，剩余遥感影像的投影带；

S2.利用基准投影带和邻接投影带中相邻遥感影像重叠部分的中心线，将用户感兴趣区域分割为S₁-S_n个子区域，并将基准投影带遥感影像和邻接投影带遥感影像分别嵌入Google S2索引空间，得到基准带索引空间和邻接带索引空间；其中，属于基准投影带的子区域为基准带子区域，属于邻接投影带的子区域为邻接带子区域；n表示分割得到的子区域总个数；

S3.根据单个子区域S_i边界点经纬度在基准投影带中的投影坐标，得到单个子区域S_i在基准带索引空间中覆盖的方格集合；i＝1,…,n；

S4.根据所有子区域S₁-S_n覆盖的方格集合在基准投影带中的投影坐标，构建影像生成区，并从遥感影像数据库中读取每个基准带子区域S_j覆盖的方格集合对应的像素值，将其写入影像生成区；其中，j＝1,…,J，J表示基准带子区域个数；

S5.根据所有邻接带子区域S_k覆盖的方格集合在邻接投影带中的投影坐标，构建邻接带影像生成区，并从遥感影像数据库中读取每个邻接带子区域S_k覆盖的方格集合对应的像素值，将其写入邻接带影像生成区；其中，k＝1,…,K，K表示邻接带子区域的个数；

S6.将邻接带子区域在基准投影带的投影坐标，重投影为邻接投影带的投影坐标，从邻接带影像生成区获取投影坐标对应的像素，写入影像生成区。

进一步地，步骤S1具体包括：

S1.1从遥感影像数据库中获取仅与用户感兴趣区域相交的遥感影像块；

S1.2计算用户感兴趣区域与属于同一投影带的遥感影像有效空间覆盖范围重叠的面积，将重叠面积最大的区域对应的投影带作为基准投影带，剩余投影带作为邻接投影带；其中，遥感影像有效空间覆盖范围表示属于同一投影带的遥感影像空间覆盖范围。

进一步地，步骤S2所述将基准投影带遥感影像和邻接投影带遥感影像分别嵌入Google S2索引空间，得到基准带索引空间和邻接带索引空间，具体包括：

01.获取基准投影带遥感影像有效空间覆盖范围与用户感兴趣区域并集的第一外接矩形，以及该矩形顶点在基准投影带中的投影坐标，将投影坐标在基准投影带中所构成四边形的第二外接矩形嵌入对应等级的Google S2索引空间，得到基准带索引空间；

其中，第二外接矩形嵌入过程中，基准带索引空间边长的计算公式为：

H＝max(W_rs,H_rs)

上式中，W_rs，H_rs分别为第二外接矩形的宽与高，ceil为向上取整函数，H_s2为基准带索引空间边长；

02.获取邻接投影带卫星遥感影像有效空间覆盖范围的第三外接矩形，以及该矩形顶点在邻接投影带中的投影坐标，将投影坐标在邻接投影带中所构成四边形的第四外接矩形嵌入Google S2索引空间，得到邻接带索引空间；

其中，第四外接矩形嵌入过程中，邻接带索引空间边长的计算公式为：

H'＝max(W'_rs,H'_rs)

上式中，W'_rs，H'_rs分别为第四外接矩形的宽与高，ceil为向上取整函数，H_s2adj为邻接带索引空间边长。

进一步地，步骤S3具体包括：

S3.1.将每个子区域S_i边界点经纬度转换成在基准投影带中的投影坐标；

S3.2.计算投影坐标在基准带索引空间中的行列号，将该行列号转换为Google S2对应的球面坐标，该球面坐标所构成范围即为子区域S_i对应的Google S2索引范围；

S3.3.利用Google S2空间索引算法，根据单个子区域S_i对应的球面坐标，计算单个子区域S_i在基准带索引空间中覆盖的方格集合。

进一步地，步骤S4具体包括：

S4.1.获取所有子区域S₁-S_n覆盖的方格集合的外接矩形，以及该外接矩形顶点在基准投影带的投影坐标，将投影坐标在基准投影带中所构成四边形的外接矩形区域作为影像生成区；

S4.2.根据基准带子区域S_j覆盖的方格集合所在位置，从遥感影像数据库中读取对应像素值，并写入影像生成区。

进一步地，步骤S5具体包括：

S5.1.获取所有邻接带子区域S_k覆盖的方格集合的外接矩形，以及该外接矩形在邻接投影带的投影坐标，将投影坐标在邻接投影带中所构成四边形的外接矩形区域作为邻接带影像生成区；

S5.2.根据邻接带子区域S_k覆盖的方格集合所在位置，从遥感影像数据库中读取对应像素值，并写入邻接带影像生成区。

优选地，上述遥感影像存储于矩阵对象数据库中。

总体而言，通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果。

(1)本发明最终结果是生成的用户感兴趣区域卫星遥感影像，减少了数据冗余，用户不必通过获取以“景”为单位的遥感影像来获取兴趣区域数据，减少了数据网络传输量与数据传输时间，同时减少了用户后续数据处难度，节省了数据处理时间。

(2)本发明在用户感兴趣区域卫星遥感影像生成过程中，无需将原始卫星遥感影像以“景”为单位全部读入内存，仅需要将由本算法确定的用户感兴趣区域覆盖的方格集合对应的遥感影像中的小影像块读入内存，参与用户感兴趣区域卫星遥感影像生成即可，减少了冗余数据的读入，进而减少了内存占用量。

(3)本发明将GoogleS2与对象数据库TileDB对接，根据GoogleS2对用户感兴趣区域的空间索引结果直接从TileDB像素级读取像素值，因此支持像素级索引。

(4)本发明将遥感影像像素位置信息嵌入到GoogleS2索引空间，实现了GoogleS2对像素级卫星遥感影像空间索引支持，该索引方法继承了GoogleS2空间索引的优点，能够快速的定位用户感兴趣区域，并快速生成用户感兴趣区域的空间覆盖范围。

(5)本算法并不改变卫星遥感影像原始的数据组织结构，无需对遥感影像进行重新组织即可使用该算法，即不需要对遥感影像进行分割分块以适应GoogleS2空间索引算法，实用性较强，不会增加遥感数据在数据库中的存储量。

附图说明

图1为本发明实施例提供的一种用户感兴趣区域遥感影像空间索引方法流程图；

图2为本发明实施例提供的用户感兴趣区域卫星遥感影像覆盖图；

图3本发明实施例提供的为用户感兴趣区域分割图；

图4(a)为本发明实施例提供的基准投影带遥感影像像素位置信息嵌入Google S2索引空间示意图；

图4(b)为本发明实施例提供的邻接投影带遥感影像像素位置信息嵌入Google S2索引空间示意图；

图5(a)为本发明实施例提供的用户感兴趣区域影像生成区示意图；

图5(b)为本发明实施例提供的邻接投影带影像生成区示意图；

图6为本发明实施例提供的用户感兴趣子区域影像生成示意图；

图7为本发明实施例提供的邻接投影带子区域影像生成示意图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。

本发明方法的基本思路是将遥感影像像素位置信息嵌入到GoogleS2索引空间(即GoogleS2中的某一层级的Cell集合)，利用Google S2快速的空间索引算法索引获取用户感兴趣区域的像素位置信息。其依据是GoogleS2索引空间与地理坐标系一一对应且形式可分离，用户感兴趣区域可映射到GoogleS2索引空间。利用GoogleS2空间索引覆盖算法获取的方格(GoogleS2中的Cell)位置，确定用户感兴趣区域数据在遥感影像中的位置。

如图1所示为本发明实施例提供的一种用户感兴趣区域遥感影像空间索引方法流程图，包括：

S1.根据用户感兴趣区域，获取基准投影带和邻接投影带；

具体地，本发明实施例基于对象存储数据库TileDB存储规范，存储卫星遥感影像数据，以支持像素级数据索引。图2为用户感兴趣区域卫星遥感影像覆盖图，图中虚线多边形为用户感兴趣区域，A、B、C、D、E、F为与之相交并完全覆盖该区域的Sentinel-2卫星遥感影像，其中A、B、C、D属于同一投影带，E、F属于同一投影带。多边形adgka为遥感影像A、B、C、D有效空间覆盖范围，多边形tuxyt为遥感影像E、F有效空间覆盖范围，其中，遥感影像有效空间覆盖范围表示属于同一投影带的的遥感影像空间覆盖范围。

本步骤中基准投影带表示属于同一投影带且与用户感兴趣区域重叠面积最大的遥感影像集所属的投影带；邻接投影带表示与用户感兴趣区域相交的遥感影像中，除去基准投影带遥感影像后，剩余遥感影像的投影带；结合图2，具体计算方法如下：计算用户感兴趣区域与多边形adgka、tuxyt的交集的面积，将交集边界经纬度坐标转化成A、B、C、D或E、F所在的投影带投影坐标，采用坐标法计算交集面积，选择交集面积大的部分对应的投影带为基准投影带。图2中用户感兴趣区域与多边形adgka的交集面积较大，因此选择A、B、C、D所在的投影带作为基准投影带，E、F遥感影像所在的投影带为邻接投影带。

S2.利用基准投影带和邻接投影带中相邻遥感影像重叠部分的中心线，将用户感兴趣区域分割为S₁-S_n个子区域，并将基准投影带遥感影像和邻接投影带遥感影像分别嵌入Google S2索引空间，得到基准带索引空间和邻接带索引空间；其中，属于基准投影带的子区域为基准带子区域，属于邻接投影带的子区域为邻接带子区域；n表示分割得到的子区域总个数；

具体地，图3中加粗虚线为相邻遥感影像重叠部分中间分割线，利用中间分割线与遥感影像的边缘线将用户感兴趣区域分割成6部分，即图中S1-S6，其中，S1-S4为基准带子区域，S5-S6为邻接带子区域，基准带子区域S1如图3中灰色部分所示。

进一步地，步骤S2中将基准投影带遥感影像和邻接投影带遥感影像分别嵌入Google S2索引空间，得到基准带索引空间和邻接带索引空间，具体包括：

01.获取基准投影带遥感影像有效空间覆盖范围与用户感兴趣区域并集的第一外接矩形，以及该矩形顶点在基准投影带中的投影坐标，将投影坐标在基准投影带中所构成四边形的第二外接矩形嵌入对应等级的Google S2索引空间，得到基准带索引空间；其中，第二外接矩形嵌入过程中，基准带索引空间边长的计算公式为：

H＝max(W_rs,H_rs)

上式中，W_rs，H_rs分别为第二外接矩形的宽与高，ceil为向上取整函数，H_s2为基准带索引空间边长；

图4(a)中矩形ad₁g₁ka即为上述第二外接矩形，宽为W_rs，高为H_rs，将其嵌入GoogleS2level＝log₂H_s2索引空间，索引空间正方形为ad₂g₂k₂a，该正方形边长为H_s2。

02.获取邻接投影带卫星遥感影像有效空间覆盖范围的第三外接矩形，以及该矩形顶点在邻接投影带中的投影坐标，将投影坐标在邻接投影带中所构成四边形的第四外接矩形嵌入Google S2索引空间，得到邻接带索引空间；

其中，第四外接矩形嵌入过程中，邻接带索引空间边长的计算公式为：

H'＝max(W'_rs,H'_rs)

上式中，W'_rs，H'_rs分别为第四外接矩形的宽与高，ceil为向上取整函数，H_s2adj为邻接带索引空间边长；

相似地，图4(b)中矩形tuxyt即为上述第四外接矩形，宽为W'_rs，高为H'_rs，将其嵌入GoogleS2level＝log₂H_s2adj索引空间，索引空间为正方形tu₁x₁y₁t，该正方形边长为为H_s2adj。

S3.根据单个子区域S_i边界点经纬度在基准投影带中的投影坐标，得到单个子区域S_i在基准带索引空间中覆盖的方格集合；i＝1,…,n；

具体地，步骤S3具体包括：

S3.1.将每个子区域S_i边界点经纬度转换成在基准投影带中的投影坐标；

S3.2.计算投影坐标在基准带索引空间中的行列号，将该行列号转换为Google S2对应的球面坐标，该球面坐标所构成范围即为子区域S_i对应的Google S2索引范围；

S3.3.利用Google S2空间索引算法，根据单个子区域S_i对应的球面坐标，计算单个子区域S_i在基准带索引空间中覆盖的方格集合。

S4.根据所有子区域S₁-S_n覆盖的方格集合在基准投影带中的投影坐标，构建影像生成区，并从遥感影像数据库中读取每个基准带子区域S_j覆盖的方格集合对应的像素值，将其写入影像生成区；其中，j＝1,…,J，J表示基准投子区域个数；

进一步地，步骤S4具体包括：

S4.1.获取所有子区域S₁-S_n覆盖的方格集合的外接矩形，以及该外接矩形顶点在基准投影带的投影坐标，将投影坐标在基准投影带中所构成四边形的外接矩形区域作为影像生成区；

如图5(a)中灰色小正方形为根据步骤S3利用GoogleS2空间索引算法在基准投影带索引空间得到的覆盖用户感兴趣各子区域的方格(即GoogleS2中的Cell)，这些方格覆盖用户感兴趣区域。图中外围矩形即为生成的影像生成区。

S4.2.根据基准带子区域S_j覆盖的方格集合所在位置，从遥感影像数据库中读取对应像素值，并写入影像生成区。

具体地，根据基准带子区域的覆盖的方格集合在基准带索引空间的位置，即在基准带索引空间行列号，以基准带索引空间左上角顶点为参考点，计算出其在基准带投影坐标，进而计算该像素在对应遥感影像中的位置，以及在影像生成区中的行列号，然后从TileDB数据库中读取对应像素值，写入图5(a)中的影像生成区，具体的，图6为基准带子区域S1对应像素写入影像生成区示意图。

S5.根据所有邻接带子区域S_k覆盖的方格集合在邻接投影带中的投影坐标，构建邻接带影像生成区，并从遥感影像数据库中读取每个邻接带子区域S_k覆盖的方格集合对应的像素值，将其写入邻接带影像生成区；其中，k＝1,…,K，K表示邻接带子区域的个数；

进一步地，步骤S5具体包括：

S5.1.获取所有邻接带子区域S_k覆盖的方格集合的外接矩形，以及该外接矩形在邻接投影带的投影坐标，将投影坐标在邻接投影带中所构成四边形的外接矩形区域作为邻接带影像生成区；

图5(b)中，灰色小正方形为根据步骤S3利用GoogleS2空间索引算法在邻接带索引空间得到的覆盖邻接带子区域的方格，这些方格覆盖邻接投影带子区域。图中外围矩形即为生成的影像生成区。

S5.2.根据邻接带子区域S_k覆盖的方格集合所在位置，从遥感影像数据库中读取对应像素值，并写入邻接带影像生成区。

具体地，根据邻接带子区域的覆盖的方格集合在邻接带索引空间的位置，即在邻接带索引空间行列号，以邻接带索引空间左上角顶点为参考点，计算出其在邻接带投影坐标，进而计算该像素在在对应遥感影像中的位置，以及在邻接带影像生成区中的行列号，然后从TileDB数据库中读取对应像素值，并写入S5.1中的邻接投影带影像生成区。

S6.将邻接带子区域在基准投影带的投影坐标，重投影为邻接投影带的投影坐标，从邻接带影像生成区获取投影坐标对应的像素，写入影像生成区。

具体地，根据邻接带子区域在基准带索引空间覆盖的方格集合位置，即在基准带索引空间行列号，以邻接带索引空间左上角顶点为参考点，计算其在基准投影带的坐标，并将该坐标转换到邻接投影带坐标，根据该坐标，采用最邻近采样法，从S5.2生成的影像中获取对应像素，写入S4.1中的影像生成区，具体地，图7为子区域S6对应像素写入影像生成区示意图。

本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (7)

1.一种用户感兴趣区域遥感影像空间索引方法，其特征在于，包括：

S1.根据用户感兴趣区域，获取基准投影带和邻接投影带；

所述基准投影带表示属于同一投影带且与用户感兴趣区域重叠面积最大的遥感影像集所属的投影带；所述邻接投影带表示与用户感兴趣区域相交的遥感影像中，除去基准投影带遥感影像后，剩余遥感影像的投影带；

S2.利用基准投影带和邻接投影带中相邻遥感影像重叠部分的中心线，将用户感兴趣区域分割为S₁-S_n个子区域，并将基准投影带遥感影像和邻接投影带遥感影像分别嵌入Google S2索引空间，得到基准带索引空间和邻接带索引空间；其中，属于基准投影带的子区域为基准带子区域，属于邻接投影带的子区域为邻接带子区域；n表示分割得到的子区域总个数；

S3.根据单个子区域S_i边界点经纬度在基准投影带中的投影坐标，得到单个子区域S_i在基准带索引空间中覆盖的方格集合；i＝1,…,n；

S4.根据所有子区域S₁-S_n覆盖的方格集合在基准投影带中的投影坐标，构建影像生成区，并从遥感影像数据库中读取每个基准带子区域S_j覆盖的方格集合对应的像素值，将其写入影像生成区；其中，j＝1,…,J，J表示基准带子区域个数；

S5.根据所有邻接带子区域S_k覆盖的方格集合在邻接投影带中的投影坐标，构建邻接带影像生成区，并从遥感影像数据库中读取每个邻接带子区域S_k覆盖的方格集合对应的像素值，将其写入邻接带影像生成区；其中，k＝1,…,K，K表示邻接带子区域的个数；

S6.将邻接带子区域在基准投影带的投影坐标，重投影为邻接投影带的投影坐标，从邻接带影像生成区获取投影坐标对应的像素，写入影像生成区。

2.根据权利要求1所述的一种用户感兴趣区域遥感影像空间索引方法，其特征在于，步骤S1具体包括：

S1.1从遥感影像数据库中获取仅与用户感兴趣区域相交的遥感影像块；

S1.2计算用户感兴趣区域与属于同一投影带的遥感影像有效空间覆盖范围重叠的面积，将重叠面积最大的区域对应的投影带作为基准投影带，剩余投影带作为邻接投影带；其中，遥感影像有效空间覆盖范围表示属于同一投影带的遥感影像空间覆盖范围。

3.根据权利要求1或2所述的一种用户感兴趣区域遥感影像空间索引方法，其特征在于，步骤S2所述将基准投影带遥感影像和邻接投影带遥感影像分别嵌入Google S2索引空间，得到基准带索引空间和邻接带索引空间，具体包括：

01.获取基准投影带遥感影像有效空间覆盖范围与用户感兴趣区域并集的第一外接矩形，以及该矩形顶点在基准投影带中的投影坐标，将投影坐标在基准投影带中所构成四边形的第二外接矩形嵌入对应等级的Google S2索引空间，得到基准带索引空间；

其中，第二外接矩形嵌入过程中，基准带索引空间边长的计算公式为：

H＝max(W_rs,H_rs)

上式中，W_rs，H_rs分别为第二外接矩形的宽与高，ceil为向上取整函数，H_s2为基准带索引空间边长；

02.获取邻接投影带卫星遥感影像有效空间覆盖范围的第三外接矩形，以及该矩形顶点在邻接投影带中的投影坐标，将投影坐标在邻接投影带中所构成四边形的第四外接矩形嵌入Google S2索引空间，得到邻接带索引空间；

其中，第四外接矩形嵌入过程中，邻接带索引空间边长的计算公式为：

H'＝max(W'_rs,H'_rs)

上式中，W'_rs，H'_rs分别为第四外接矩形的宽与高，ceil为向上取整函数，H_s2adj为邻接带索引空间边长。

4.根据权利要求1-3任一项所述的一种用户感兴趣区域遥感影像空间索引方法，其特征在于，步骤S3具体包括：

S3.1.将每个子区域S_i边界点经纬度转换成在基准投影带中的投影坐标；

S3.2.计算投影坐标在基准带索引空间中的行列号，将该行列号转换为Google S2对应的球面坐标，该球面坐标所构成范围即为子区域S_i对应的Google S2索引范围；

S3.3.利用Google S2空间索引算法，根据单个子区域S_i对应的球面坐标，计算单个子区域S_i在基准带索引空间中覆盖的方格集合。

5.根据权利要求1-4任一项所述的一种用户感兴趣区域遥感影像空间索引方法，其特征在于，步骤S4具体包括：

S4.1.获取所有子区域S₁-S_n覆盖的方格集合的外接矩形，以及该外接矩形顶点在基准投影带的投影坐标，将投影坐标在基准投影带中所构成四边形的外接矩形区域作为影像生成区；

S4.2.根据基准带子区域S_j覆盖的方格集合所在位置，从遥感影像数据库中读取对应像素值，并写入影像生成区。

6.根据权利要求1-5任一项所述的一种用户感兴趣区域遥感影像空间索引方法，其特征在于，步骤S5具体包括：

S5.1.获取所有邻接带子区域S_k覆盖的方格集合的外接矩形，以及该外接矩形在邻接投影带的投影坐标，将投影坐标在邻接投影带中所构成四边形的外接矩形区域作为邻接带影像生成区；

S5.2.根据邻接带子区域S_k覆盖的方格集合所在位置，从遥感影像数据库中读取对应像素值，并写入邻接带影像生成区。

7.根据权利要求1-6任一项所述的一种用户感兴趣区域遥感影像空间索引方法，其特征在于，所述遥感影像存储于对象数据库中。

Images