A Spatio-Temporal Local Association Query Algorithm for Multi-Source Remote Sensing Big Data

: It is extremely important to extract valuable information and achieve efﬁcient integration of remote sensing data. The multi-source and heterogeneous nature of remote sensing data leads to the increasing complexity of these relationships, and means that the processing mode based on data ontology cannot meet requirements any more. On the other hand, the multi-dimensional features of remote sensing data bring more difﬁculties in data query and analysis, especially for datasets with a lot of noise. Therefore, data quality has become the bottleneck of data value discovery, and a single batch query is not enough to support the optimal combination of global data resources. In this paper, we propose a spatio-temporal local association query algorithm for remote sensing data (STLAQ). Firstly, we design a spatio-temporal data model and a bottom-up spatio-temporal correlation network. Then, we use the method of partition-based clustering and the method of spectral clustering to measure the correlation between spatio-temporal correlation networks. Finally, we construct a spatio-temporal index to provide joint query capabilities. We carry out local association query efﬁciency experiments to verify the feasibility of STLAQ on multi-scale datasets. The results show that the STLAQ weakens the barriers between remote sensing data, and improves their application value effectively.


Research Background
With the continuous development of Earth Observation Network technology, various types of and huge amounts of remote sensing data are generated by a large number of sensors in real time. The magnitude of remote sensing data has increased from GBs to TBs and PBs, and it will continue to increase in the future [1]. These remote sensing data are collections of geographic information related to the location, based on a unified space-time reference [1][2][3]. Remote sensing big data is the fusion of big data and remote sensing data, and is active in time and space. It is big data based on a unified space-time reference, which takes the Earth as the object [1][2][3]. Remote sensing big data mainly includes space-time reference data, geodetic survey data, gravity and magnetic data, remote sensing image data, and location-related spatial media data. In addition, it has been widely used in many fields such as national defense, agriculture, water conservancy, land planning, smart cities, disaster warning, geological surveys, emergency monitoring and so on [3][4][5]. Although these massive remote sensing data come from different sources and have different structures, they are often potentially related to each other due to their own spatio-temporal 1. In view of the general problems in remote sensing big data management and retrieval, we design a general and efficient spatio-temporal local association query algorithm for remote sensing data, the STLAQ, on the basis of mining and constructing the relationships of remote sensing data from different perspectives.

2.
The STLAQ solves the association problem of high-dimensional remote sensing data, relying on a data model and correlation network, especially when global features degenerate into local features on a unified multi-scale organization model. 3.
The STLAQ weakens the barriers between remote sensing data and expands their sharing mode, and has strong versatility and practical value in the field of remote sensing data application.

Paper Organization
In Section 1, we introduce the main problems faced by remote sensing big data in applications. In Section 2, we introduce the research situation of remote sensing data index and remote sensing data correlation. Then, in Section 3, we present the main structure and concepts of the STLAQ, introducing the STLAQ from two aspects; multisource remote sensing data association and local association query based on correlation network. In Section 4, we construct the STCN on several typical remote sensing datasets and carry out association query experiments to prove the usability and universality of the STLAQ. Finally, we conclude this study and look forward to the follow-up work in Section 5.

Related Work
With the rapid growth of the volume and the increasingly rich sources and types, remote sensing data continues to be applied in various fields. In the applications of remote sensing big data, there are challenges in storage, management, analysis and visualization. It is extremely important to prove an efficient query algorithm for massive multi-source remote sensing data. Traditional remote sensing data query algorithms only provide retrieval schemes for data coming from a single-source, lacking a unified retrieval mechanism for multi-source and multi-type remote sensing data. Designing a remote sensing data index with high-performance and correlation method for multi-source remote sensing data, to solve the barriers caused by the heterogeneity of multi-source remote sensing data, is one of the focuses of the research field of remote sensing big data.
Research on spatio-temporal index is usually for the improvement and expansion of traditional spatial indexes in dimensions, such as time dimension. The earliest research can be traced back to the 1970s [9][10][11]. The Quadtree index is proposed in [12], dividing nodes into four regions in a two-dimensional space to solve the problem of combined query of data. In [13], a multi-dimensional binary tree (KD tree) in which each node is a k-dimensional point is introduced, and it could be used to store k-dimensional data; R tree is designed in [14], a hierarchical data structure based on B+ trees, used for the dynamic organization of a set of d-dimensional geometric objects. SDMR tree, a deformed R-tree index structure for multi-scale expression of spatial data was proposed in [15]. In [16], the Hilbert R tree is proposed, which uses Hilbert curve, which has better spatial proximity than other space filling curves [17], in order to perform one-dimensional sorting on k-dimensional space data to improve node storage utilization. In [18], a distributed remote sensing database with indexing techniques, GeoMesa, is presented. It uses space filling curves to map multi-dimensional data to the single lexicographic list managed by the underlying distributed database. In [19], the authors propose a design and implementation of a distributed spatial query index and query algorithm consisting of a hybrid index composed of a Quadtree, R-Tree and Hash structure. In [20], a unified index framework is proposed, which provided efficient data distribution, fault tolerance and multi-dimensional data query processing, and proposed an algorithm for range query and K-nearest neighbor query. In [21], a spatial index is designed, linked with LSM B-tree and LH R-tree, which effectively improved the operation efficiency of the entire index for update and retrieval. In [22], the authors combine grid coding, dimensional attributes of interest and HBase row keys to design a customizable multi-dimensional index structure, ST+. In [23], a remote sensing data indexing method based on Hilbert code is presented, solving the problem of fast access to remote sensing data in the virtual battlefield environment.
Data association is to mine the correlation and dependence between data, and to obtain the internal connection between multi-source geographic data [24]. Geographical data association can realize the overall holographic relationship expression of natural and human geographic elements in the region [25]. Chen Luo's team has carried out a series of studies. In [26], the authors proposea generalized spatial data correlation model for spatial data, place names and general entities, and verifies the effectiveness of the model. In [27], a multi-information correlation database is established for several types of aerospace information, such as remote sensing image information, radiation source information, geographic information, etc. In [28], the authors introduce graph theories and techniques into entity relationship research, and propose a graph-based correlation analysis technique. In [29], a data association model and a geocoding method based on space-time unity are proposed, to realize the association organization and analysis of multi-source heterogeneous data based on a unified space-time framework. However, these studies also have shortcomings. GeoStar proposes a solution of a remote sensing data platform, which mainly focuses on the connection between geographic elements and social thematic data, and establishes a unified data backplane. In these studies, only spatial, visual, and semantic data correlations are considered, with a lack of dynamic relationships in such time and events dimensions. Moreover, these correlation models are all performed at a single scale, ignoring correlation across different scales.
All of the above studies have launched the research on remote sensing data association methods from different angles. In any case, these studies only pay attention to the design and optimization of the model structure, ignoring the influence of the remote sensing data's own features on the quality of the association query results. There may be the disadvantage of incomplete association query results due to the lack of data features at a single scale. Therefore, there are still gaps to be filled in the multi-scale association query combining data features. The STDM is a digital model established to facilitate the description and analysis of remote sensing data. It consists of three basic elements: space-time reference (STR), structure description (SD) and feature description (FD). On the basis of the STDM encapsulation, we propose the spatio-temporal correlation network (STCN). It is a modelized description of spatiotemporal data and its association relationship based on the undirected hypergraph model. It includes two parts, spatio-temporal self-correlation network (STSCN) and spatio-temporal cross-associated network (STCCN). The spatio-temporal association index is an index defined according to the structure of the spatio-temporal correlation network, including two parts: a multi-dimensional feature index and a correlation index.

Spatio-Temporal Data Model
The spatio-temporal data model is a digital model proposed to facilitate data description and analyzation, as shown in Figure 2. The spatio-temporal data model can be formally expressed as ( ) ,, STM STR SD FD . Based on the spatio-temporal data model, remote sensing data can be encapsulated into two different granularities, remote sensing data unit (STDU) and remote sensing data object (STDO). Remote sensing data can be divided into one or more STDOs, according to definite classification standards. An STDU The STDM is a digital model established to facilitate the description and analysis of remote sensing data. It consists of three basic elements: space-time reference (STR), structure description (SD) and feature description (FD). On the basis of the STDM encapsulation, we propose the spatio-temporal correlation network (STCN). It is a modelized description of spatiotemporal data and its association relationship based on the undirected hypergraph model. It includes two parts, spatio-temporal self-correlation network (STSCN) and spatio-temporal cross-associated network (STCCN). The spatio-temporal association index is an index defined according to the structure of the spatio-temporal correlation network, including two parts: a multi-dimensional feature index and a correlation index.

Spatio-Temporal Data Model
The spatio-temporal data model is a digital model proposed to facilitate data description and analyzation, as shown in Figure 2. The spatio-temporal data model can be formally expressed as STM(STR, SD, FD). Based on the spatio-temporal data model, remote sensing data can be encapsulated into two different granularities, remote sensing data unit (STDU) and remote sensing data object (STDO). Remote sensing data can be divided into one or more STDOs, according to definite classification standards. An STDU is the smallest unit in remote sensing data processing, and is an element in an STDO. An STDO can be expressed as a collection of STDUs organized in a hierarchical model. Space-time Reference. Space-time reference is the basis for describing and measuring remote sensing data, including two parts; time reference and space reference. The origin and measurement scale of the datum need to be clarified in each reference. The space reference system includes several parts, such as plane coordinate system, map projection, and height datum, used in representing the spatial feature of remote sensing data and describing the spatial relationship between them. The time reference system is the basis for time information exchange and calculation. Remote sensing data can be expressed as a moment or a period in the time reference system. In [30], the Gregorian calendar, 24:00, local time or coordinated universal time are designated as the basic time reference system for information exchange. • Structure Description. The structure description of remote sensing data is a description of the underlying remote sensing data organization scheme, including multi-scale model, storage structure and index structure. The multi-scale model is a multi-levelof-detail model, established by way of breaking up different types of data according to uniform rules and dividing the data into different resolution levels. In the multi-scale model, spatio-temporal data units at different resolutions can be correlated based on features, and be stored and indexed uniformly. The storage structure defines the storage form of remote sensing data, to organize data and its characteristics in a certain structure. In order to solve the efficiency problem of massive remote sensing data in joint retrieval, it is necessary to build indexes on the basis of storage structure. The index structure describes the composition of the spatio-temporal index, including retrieval conditions, retrieval values, and keys of the remote sensing data in storage structure. • Feature Description. Feature description is the low-level features exacted from remote sensing data. The relationship between remote sensing data is essentially to calculate the similarity between the various features and then to aggregate the results. Therefore, the extraction of data features and the calculation of the correlation degree are critical parts of the construction of data correlation. This paper divides the features of remote sensing data into three categories: time feature, spatial feature, and data feature. In this case, time feature can be expressed as the production time or the life cycle of remote sensing data; spatial features include two parts, spatial position (such as latitude and longitude) and spatial form (such as points, lines, surfaces). Spatial and time features of the data can be obtained from the data itself, as well as through data analysis. Data features refer to the low-level features extracted from remote sensing data. is the smallest unit in remote sensing data processing, and is an element in an STDO. An STDO can be expressed as a collection of STDUs organized in a hierarchical model.

•
Space-time Reference. Space-time reference is the basis for describing and measuring remote sensing data, including two parts; time reference and space reference. The origin and measurement scale of the datum need to be clarified in each reference. The space reference system includes several parts, such as plane coordinate system, map projection, and height datum, used in representing the spatial feature of remote sens- An example of encapsulating digital elevation model data into a spatio-temporal data model is given in Table 1. The detailed explanations are as follows.

1.
In the part of space-time reference, the spatial reference is WGS84/EGM96, and the time reference is UTC/GMT+08:00 zone. 2.
In the part of structure description, the multi-scale model of the sample data is a tile pyramid model. The key of the storage table is composed of data location, data type, and data acquisition time. There are 2 column families in the storage table, where data is used to store tile data, meta is used to store the meta information. The index structure includes a multi-dimensional feature index and an association index. 3.
In the part of feature description, the time feature is the data acquisition time expressed in the format of "yyyy-mm-dd hh:mm:ss". The spatial feature is the geographic region of the sample data, and the data feature is the low-level features, such as the maximum, minimum, mean, standard deviation and other statistical features.

Spatio-Temporal Correlation Network
There is a certain degree of correlation between different remote sensing data, which can be described with the spatio-temporal correlation network (STCN), a model based on the undirected hypergraph model. In the STCN, vertices are used to represent spatiotemporal data objects, hyper-edges are used to represent complex relationships among them, and the weights of hyper-edges are used to define the strength of the correlation between spatio-temporal data objects. The strength of edges between spatio-temporal data objects can be expressed by the correlation degree of remote sensing data, which we will describe in Section 3.2.2.
There are two kinds of remote sensing data correlation network; spatio-temporal self-correlation network and spatio-temporal cross-correlation network. In this case, the STSCN describes the correlation among STDUs within an STDO, and can be regarded as an indivisible hyperedge in the STCN, and the STCCN describes the correlation between different STDOs, that is, the aggregation of multiple STSCN. The structure of STCN is shown in Figure 3. A remote sensing data object may belong to multiple networks, and the variables marked on the connection line represent the degree of correlation between STDOs and STCNs. For example, a lake group with a high degree of spatio-temporal correlation can be represented by constructing a STCCN. The lakes in the lake group can be represented by the STSCNs corresponding to them in the STCCN.

Preliminaries
• Space Pyramid [31,32]. The spatial pyramid technology is that where the pyramid is a multi-resolution hierarchical model. The essence of its construction process is to block and layer the image to form multiple resolution levels.

•
Hypergraph Model [33]. The graph model is a typical model that can effectively represent entities and their relationships. A graph can be expressed as  [34]. Spectral clustering turns the clustering problem into the graph partitioning problem. To solve the graph cut objective function, the properties of the Rayleigh quotient are usually utilized to map the original data points into a lower dimensional eigen-space by calculating the eigenvectors of Laplacian matrix and then conducting the clustering in the new space [34]. In spectral clustering, all data can be treated as vertices connected by edges. The edge weight value between two points farther apart is lower, and the edge weight value between two points closer together is higher. The purposes of clustering are to make the weight of the edges between the different subgraphs as low as possible, and to make the weight of the edges in the subgraphs as high as possible. The method of spectral clustering used in this paper is normalized cut. According to the implementation of spectral clustering given in [35,36], the application of spectral clustering in STCN construction can be found in Section 3.3.2.

Preliminaries
• Space Pyramid [31,32]. The spatial pyramid technology is that where the pyramid is a multi-resolution hierarchical model. The essence of its construction process is to block and layer the image to form multiple resolution levels.

•
Hypergraph Model [33]. The graph model is a typical model that can effectively represent entities and their relationships. A graph can be expressed as G = (V, E), its vertex set is V(G) = {v 1 , v 2 , · · · , v n }, and the edge set is E = {e 1 , e 2 , · · · , e m }. In the graph, each edge is a point pair (v, w), where v, w ∈ V. There are two kinds of graphs: directed graph and undirected graph [33]. If the (v, w) is ordered, G is a directed graph, otherwise G is an undirected graph. Hypergraph is a generalization of ordinary graphs and can be used to express complex relationships. For a hypergraph H = (V, E), the set of hyperedges is E = {e 1 , e 2 , · · · , e m }, and each hyperedge e is a subset of the set of vertices V, which can connect two or more vertices and meet the regulation e 1 ∪ e 2 ∪ · · · ∪ e m = V. A weighted hypergraph is to assign a weight as w(e) to each hyperedge, indicating the degree to which the vertices within belong to the hyperedge e. Adjacency matrix and incidence matrix are always used to describe structure of graph G. The adjacency matrix of graph G is a matrix of size n × n, denoted as A(G). If there is an edge between the vertex v i and v j , the value of a i,j is 1, otherwise is 0. The incidence matrix of graph G is a matrix of size n × m, denoted as R(G). If the vertex v i is an endpoint of the edge e j , the value of r i,j is 1, otherwise it is 0. • Spectral Cluster [34]. Spectral clustering turns the clustering problem into the graph partitioning problem. To solve the graph cut objective function, the properties of the Rayleigh quotient are usually utilized to map the original data points into a lower dimensional eigen-space by calculating the eigenvectors of Laplacian matrix and then conducting the clustering in the new space [34]. In spectral clustering, all data can be treated as vertices connected by edges. The edge weight value between two points farther apart is lower, and the edge weight value between two points closer together is higher. The purposes of clustering are to make the weight of the edges between the different subgraphs as low as possible, and to make the weight of the edges in the subgraphs as high as possible. The method of spectral clustering used in this paper is normalized cut. According to the implementation of spectral clustering given in [35,36], the application of spectral clustering in STCN construction can be found in Section 3.3.2.

Multi-Source Remote Sensing Data Association
As shown in Figure 4, there are two main parts in the process of multi-source remote sensing data association, spatio-temporal data model encapsulation and spatio-temporal correlation network construction, which consists of spatio-temporal self-correlation network construction and spatio-temporal cross-correlation network construction. A stepby-step feedback mechanism based on association query results and manual marking is used to adjust the values of parameters in the whole process. Then, we will introduce them in detail.

Multi-Source Remote Sensing Data Association
As shown in Figure 4, there are two main parts in the process of multi-source remote sensing data association, spatio-temporal data model encapsulation and spatio-temporal correlation network construction, which consists of spatio-temporal self-correlation network construction and spatio-temporal cross-correlation network construction. A stepby-step feedback mechanism based on association query results and manual marking is used to adjust the values of parameters in the whole process. Then, we will introduce them in detail.

Spatio-Temporal Data Model Encapsulation
To complete the encapsulation of the spatio-temporal data model, data analysis, feature extraction, and multi-scale division suitable for remote sensing data of various types and formats are required.
In order to perform remote sensing data feature extraction, firstly, we extract color features, shape features and other visual features, such as grayscale histogram, grayscale mean, directional gradient histogram, and image principal component features [37][38][39][40], to establish the data feature description. Secondly, we read the time stamp information from the metadata of the data, to establish the time feature description. Finally, we read the coordinate system information, origin coordinates and resolutions to calculate the geographic extent of the data and to establish a spatial feature description.
As for remote sensing image data, we can establish a remote sensing data grid model based on the spatial pyramid [41,42] for multi-scale division, which is the basis for later associating feature on spatio-temporal data units. Then, we split remote sensing data into data units, build storage structure and index structure on them, and perform local correlation of remote sensing data features on data units of different scales.
As shown in Figure 5, we take lake remote sensing image data as an example, and to explain, the construction of a remote sensing data grid, based on spatial pyramid, is performed.
1. Associate spatio-temporal information on the basis of the pyramidal grid, and establish the correspondence between remote sensing data and grid models according to the spatial resolution and geographic scope of the remote sensing data. 2. In the global remote sensing data grid, the size of the grid cell is fixed (such as 256 × 256), and the spatial resolution of grid cells at different levels can be obtained via

Spatio-Temporal Data Model Encapsulation
To complete the encapsulation of the spatio-temporal data model, data analysis, feature extraction, and multi-scale division suitable for remote sensing data of various types and formats are required.
In order to perform remote sensing data feature extraction, firstly, we extract color features, shape features and other visual features, such as grayscale histogram, grayscale mean, directional gradient histogram, and image principal component features [37][38][39][40], to establish the data feature description. Secondly, we read the time stamp information from the metadata of the data, to establish the time feature description. Finally, we read the coordinate system information, origin coordinates and resolutions to calculate the geographic extent of the data and to establish a spatial feature description.
As for remote sensing image data, we can establish a remote sensing data grid model based on the spatial pyramid [41,42] for multi-scale division, which is the basis for later associating feature on spatio-temporal data units. Then, we split remote sensing data into data units, build storage structure and index structure on them, and perform local correlation of remote sensing data features on data units of different scales.
As shown in Figure 5, we take lake remote sensing image data as an example, and to explain, the construction of a remote sensing data grid, based on spatial pyramid, is performed.

1.
Associate spatio-temporal information on the basis of the pyramidal grid, and establish the correspondence between remote sensing data and grid models according to the spatial resolution and geographic scope of the remote sensing data.

2.
In the global remote sensing data grid, the size of the grid cell is fixed (such as 256 × 256), and the spatial resolution of grid cells at different levels can be obtained via where lon max and lon min represent the maximum and minimum value of the longitude range of the current coordinate system, respectively, level represents the level of current grid, gridRes represents the resolution, and gridSize represents the size of current grid. 3.
Copy the remote sensing data to the relevant position in the grid based on the correspondence to obtain a grid generated from the remote sensing data.

4.
Use resampling algorithm (such as bilinear interpolation algorithm) to down-sample remote sensing data step-by-step, to obtain grid data at lower resolution levels, obtaining a complete remote sensing data pyramid.

5.
Multi-source remote sensing data is sampled and processed on a unified remote sensing data pyramid hierarchy to obtain global spatio-temporal data units with multiple levels of detail.
where max lon and min lon represent the maximum and minimum value of the longitude range of the current coordinate system, respectively, level represents the level of current grid, gridRes represents the resolution, and gridSize represents the size of current grid. 3. Copy the remote sensing data to the relevant position in the grid based on the correspondence to obtain a grid generated from the remote sensing data. 4. Use resampling algorithm (such as bilinear interpolation algorithm) to down-sample remote sensing data step-by-step, to obtain grid data at lower resolution levels, obtaining a complete remote sensing data pyramid. 5. Multi-source remote sensing data is sampled and processed on a unified remote sensing data pyramid hierarchy to obtain global spatio-temporal data units with multiple levels of detail.

Spatio-Temporal Correlation Network Construction
According to the hierarchical model of the spatio-temporal correlation network shown in Figure 3, the construction of network consists of two steps: the construction of spatio-temporal self-correlation network and of cross-correlation network. In STCN, an STDO is usually represented as a set of STDUs divided at different scales. Each unit in the collection inherits some of the features of the remote sensing data and may be related to other datasets. In order to obtain the relationships among STDOs, it is necessary to aggregate the relationships among STDUs within the objects. Moreover, on the basis of calculating the similarity of the local features of the STDU, STCN can change the conventional single method of establishing remote sensing data correlation based on the same scale, and establish cross-scale correlation.
The relationships among three STDOs, I U , J U and K U , are shown in Figure 6.
In Figure 6, spherical models of different sizes are used to represent STDUs at different scales. The correlation between spatio-temporal data units, 1 I U and 1 J U , can be calculated by Equation (6), which is an aggregation result of spatial similarity, temporal similarity, and data similarity.

Spatio-Temporal Correlation Network Construction
According to the hierarchical model of the spatio-temporal correlation network shown in Figure 3, the construction of network consists of two steps: the construction of spatiotemporal self-correlation network and of cross-correlation network. In STCN, an STDO is usually represented as a set of STDUs divided at different scales. Each unit in the collection inherits some of the features of the remote sensing data and may be related to other datasets. In order to obtain the relationships among STDOs, it is necessary to aggregate the relationships among STDUs within the objects. Moreover, on the basis of calculating the similarity of the local features of the STDU, STCN can change the conventional single method of establishing remote sensing data correlation based on the same scale, and establish cross-scale correlation.
The relationships among three STDOs, U I , U J and U K , are shown in Figure 6. In Figure 6, spherical models of different sizes are used to represent STDUs at different scales. The correlation between spatio-temporal data units, U I1 and U J1 , can be calculated by Equation (6), which is an aggregation result of spatial similarity, temporal similarity, and data similarity.
where max lon and min lon represent the maximum and minimum value of the longitude range of the current coordinate system, respectively, level represents the level of current grid, gridRes represents the resolution, and gridSize represents the size of current grid. 3. Copy the remote sensing data to the relevant position in the grid based on the correspondence to obtain a grid generated from the remote sensing data. 4. Use resampling algorithm (such as bilinear interpolation algorithm) to down-sample remote sensing data step-by-step, to obtain grid data at lower resolution levels, obtaining a complete remote sensing data pyramid. 5. Multi-source remote sensing data is sampled and processed on a unified remote sensing data pyramid hierarchy to obtain global spatio-temporal data units with multiple levels of detail.

Spatio-Temporal Correlation Network Construction
According to the hierarchical model of the spatio-temporal correlation network shown in Figure 3, the construction of network consists of two steps: the construction of spatio-temporal self-correlation network and of cross-correlation network. In STCN, an STDO is usually represented as a set of STDUs divided at different scales. Each unit in the collection inherits some of the features of the remote sensing data and may be related to other datasets. In order to obtain the relationships among STDOs, it is necessary to aggregate the relationships among STDUs within the objects. Moreover, on the basis of calculating the similarity of the local features of the STDU, STCN can change the conventional single method of establishing remote sensing data correlation based on the same scale, and establish cross-scale correlation.
The relationships among three STDOs, Figure 6.
In Figure 6, spherical models of different sizes are used to represent STDUs at different scales. The correlation between spatio-temporal data units, 1 I U and 1 J U , can be calculated by Equation (6), which is an aggregation result of spatial similarity, temporal similarity, and data similarity.  We use the adjacency tensor to represent the correlation between different STDOs in the multi-dimensional feature space, as shown in Figure 7. In Figure 7a, the adjacency tensor of the remote sensing data cross-correlation network simplified in matrix form is given. The structure of the adjacency tensor between the spatio-temporal data objects U I and U J is given in Figure 7b, identified with the red rectangles. In Figure 7c, the mapping relationship between the similarity of the remote sensing data unit in time, spatial and data feature, and the data relevance is given in β T , β S and β D , respectively, representing the contribution weight of time feature, spatial feature and data feature in the calculation of relevance, using different colors to indicate the difference in the correlation degree of spatio-temporal data units under different feature similarities.
the multi-dimensional feature space, as shown in Figure 7. In Figure 7a, the adjacency tensor of the remote sensing data cross-correlation network simplified in matrix form is given. The structure of the adjacency tensor between the spatio-temporal data objects I U and J U is given in Figure 7b, identified with the red rectangles. In Figure 7c, the mapping relationship between the similarity of the remote sensing data unit in time, spatial and data feature, and the data relevance is given in  T ,  S and  D , respectively, representing the contribution weight of time feature, spatial feature and data feature in the calculation of relevance, using different colors to indicate the difference in the correlation degree of spatio-temporal data units under different feature similarities. The Construction of the Self-Correlation Network Firstly, in the process of self-correlation network construction, we choose some features of remote sensing data as a set of vectors. Then, we used the optimized FCM algorithm [43,44] to analyze the relationship between data containing the same feature from the grid dataset, to obtain a set of remote sensing datasets containing different geographic elements. The similarity can be obtained according to the distance of the grid data in the feature space, as the correlation degree.
The similarity between two different spatio-temporal data units in characteristics can reveal their potential correlation. Therefore, the correlation degree of remote sensing data can be measured by the similarity of features. Considering that remote sensing data has multiple features of different dimensions, it is necessary to calculate the similarity in different feature spaces to construct the joint similarity. Therefore, in the self-correlation network, metrics involved in the data relationship include the similarity of data, time, and space features, the joint similarity of multi-dimensional features, and correlation degree. •

Similarity in Data Features
Considering the difference in sources and types of remote sensing data, data features can usually be extracted, and can be regarded as different high-dimensional feature spaces. In order to describe the relationship of remote sensing data in data features, a general representation method is needed to calculate the similarity degree between remote sensing data of the same type based on data features, which can usually be expressed by the distance in the data feature space of remote sensing data.
With reference to the definition of Gaussian kernel function [35] and Euclidean distance [34], the similarity in data features between remote sensing data unit u and v is defined as: The Construction of the Self-Correlation Network Firstly, in the process of self-correlation network construction, we choose some features of remote sensing data as a set of vectors. Then, we used the optimized FCM algorithm [43,44] to analyze the relationship between data containing the same feature from the grid dataset, to obtain a set of remote sensing datasets containing different geographic elements. The similarity can be obtained according to the distance of the grid data in the feature space, as the correlation degree.
The similarity between two different spatio-temporal data units in characteristics can reveal their potential correlation. Therefore, the correlation degree of remote sensing data can be measured by the similarity of features. Considering that remote sensing data has multiple features of different dimensions, it is necessary to calculate the similarity in different feature spaces to construct the joint similarity. Therefore, in the self-correlation network, metrics involved in the data relationship include the similarity of data, time, and space features, the joint similarity of multi-dimensional features, and correlation degree.

• Similarity in Data Features
Considering the difference in sources and types of remote sensing data, data features can usually be extracted, and can be regarded as different high-dimensional feature spaces. In order to describe the relationship of remote sensing data in data features, a general representation method is needed to calculate the similarity degree between remote sensing data of the same type based on data features, which can usually be expressed by the distance in the data feature space of remote sensing data.
With reference to the definition of Gaussian kernel function [35] and Euclidean distance [34], the similarity in data features between remote sensing data unit u and v is defined as: where r represents the dimension of common data features of remote sensing data unit u and v, F i (u) and F i (v) represents the value of the i-th dimension data feature of u and v respectively, F i (u) − F i (v) 2 is the Euclidean distance of u and v in the r-dimensional feature space, and α i represents the weight of the i-th dimension data feature in the similarity judgment.

• Similarity in Time Features
Remote sensing data always has time features when it is generated. Considering that the change in the time dimension is continuous, time features of a data can be quantified at different time resolutions according to the need. The similarity in time features of remote sensing data can be represented with difference of the quantified time feature, which is: where T(u) and T(v) represent the quantified time characteristic values of u and v respectively, such as the data production time expressed in milliseconds.

• Similarity in Spatial Features
Remote sensing data usually has a description of spatial features, which represent the spatial location and scope of remote sensing data. By calculating the spatial proximity of the remote sensing data, the similarity of the spatial features of the data can be obtained, and the correlation on the spatial features can be established as well.
Taking into account the continuity of remote sensing data in the spatial dimension, in order to facilitate the compression coding and distance measurement of the spatial features of the data, the space filling curve is used to fill the remote sensing data grid, converting a remote sensing data grid to a point on the space filling curve. Establish the mapping relationship between spatial features and points on the space filling curve in measuring the spatial similarity of spatio-temporal data units, then encode the points and use the size of the matched coding bits and the corresponding spatial resolution to obtain.
Define spatial similarity between spatio-temporal data units u and v as: where S(u) and S(v) represent the spatial feature coding of u and v respectively.

• Joint similarity
In order to express correlation degree of spatio-temporal data units, treat the similarity of remote sensing data in data features, temporal features, and spatial features as the attributes of the remote sensing data correlation in three-dimensional space. Therefore, the joint similarity can be expressed as: where β T , β S and β D represent the weight of the similarity of the time, space and data features of remote sensing data in the relationship of the remote sensing data, respectively. In the actual application process, β T , β S and β D can be dynamically adjusted according to the requirements of data relational retrieval, so as to construct the correlation in a single feature and in combinations of features in different dimensions.

• Correlation degree
We use the similarity of spatio-temporal data units in multiple dimensions to express the correlation degree between remote sensing data. Meanwhile, in order to reduce the cost in storage and computing, the threshold ε is specified when calculating the correlation degree. When it is less than ε, the correlation degree of two remote sensing data is considered to be 0, that is, there is no correlation between them; when it is greater than or equal to ε, the joint similarity of the spatio-temporal data model can be used to express their correlation degree. The determination of the value of ε is achieved through a feedback adjustment mechanism. Firstly, we set an initial value for ε and calculate the correlation degree and build STCN. Then we performed associative queries based on STCN, and compare the results of the associated queries with the manually marked results to obtain the direction and offset of ε in further adjustment. Next, we use the adjusted ε to re-execute the above process until the result of the associated query is similar to the result of the manual mark to the greatest extent.
In summary, we define the correlation degree of the remote sensing data unit as: where The corresponding pseudo-code can be found in Algorithm 1.

Algorithm 1.
The algorithm implementation of the multi-source remote sensing data correlation method.

STSCN Construction Algorithm
Use optimized FCM algorithm for clustering to get E ; 4. for e l ∈ E 5. e l = ∅ 6. for u i , u j ∈ e l 7.
Calculate a i,j , r i,l ; 10.

The Construction of the Cross-Correlation Network
Based on the self-correlation network, we analyze the correlation degree of the remote sensing data within different self-correlation networks to establish the correlation and obtain a wider correlation network of remote sensing data, in order to provide full-time, full-type association query capabilities for remote sensing data.
On the basis of spectral clustering, the construction process of STCN can be seen as the correlation graph cut to obtain different subgraphs as different correlation networks. In the process of graph cut, the main goal is to make the correlation degrees between spatio-temporal data objects in the same subgraph as high as possible, and the correlation degrees between spatio-temporal data objects in different subgraphs are as low as possible.
Firstly, we take spatio-temporal data objects as nodes to generate a spatio-temporal cross-correlation network based on spectral clustering. Remote sensing data is a subgraph of the remote sensing data object graph. In order to divide the remote sensing data object with the highest degree of correlation into the same correlation network, the remote sensing data object correlation graph needs to be cut to generate k subgraphs that are not connected to each other. The set of each subgraph is defined as Referring to the graph cutting process in the literature [45], the generation process of the definition is as follows.

1.
Establish similarity matrix, adjacency matrix and degree matrix for the correlation among STDOs. The adjacency tensor of STCCN can be degenerated into an adjacency matrix A = a i,j , let As for U i in STCCN, the degree d i is defined as the sum of the weights of all edges connected to it, that is The Laplacian matrix of STCCN can be calculated from the degree matrix and the adjacency matrix as L = D − W.

2.
Calculate and standardize the Laplacian matrix of the correlation graph of STDOs; 3.
Calculate the smallest eigenvalue and its corresponding eigenvector on the standardized matrix; 4.
Standardize obtained eigenvectors by rows to form a characteristic matrix of size n × k; 5.
Regard matrix M as n k-dimensional samples, perform FCM clustering, and divide the complete set of spatio-temporal data objects into c clusters, C = {c 1 , c 2 , · · · , c c }, corresponding to c STCCNs.
In the above process, a set of STDOs can be divided into c STCCNs, and each STDO belongs to only one STCCN. However, in practical applications, a piece of STDOs may belong to different STCCNs, with different correlation degrees. So, it is convenient to control the data returned meeting different association query conditions. Therefore, it is necessary to expand the STCCN obtained in the above clustering process. In summary, the expansion process of STCCN is as follows. The corresponding pseudo-code can be found in Algorithm 2.
Take a remote sensing data object o i from c j , and calculate the correlation degree of o i to the element c l in c j as where c j is the complement of c j , that is, the union of other cross-correlation networks in C except c j ; 3.
If relation(o i , c l ) obtained in step 2 is greater than the threshold ε c , it is considered that o i belongs to c l , let c l = c l ∪ {o i }; 4.
Execute step 2 to step 4 cyclically, until all the STDOs in c j have verified the correlation of each subnet; 5.
Execute step 1 to step 4 cyclically, until all the spatio-temporal cross-correlation networks in C have completed the completion operation.

Algorithm 2.
The algorithm implementation of the multi-source remote sensing data correlation method.

STCCN Construction Algorithm
10. Get the smallest k eigenvalues and eigenvectors of L; 11. Standardize the eigenvectors by rows to obtain an eigenmatrix M of size |O| × k; 12. Consider M as k-dimensional samples, and use FCM for clustering to obtain the divided cluster set C; Calculate a i,j , r i,l ; 24. A In the above two sections, multiple adjustable parameters are involved. In process of self-correlation network construction, there are adjustable parameters such as the contribution weights β T , β S , and β D in the calculation of the correlation degree, and the widths σ T , σ S , and σ D of the Gaussian kernel function, and the judgment of the correlation degree of spatio-temporal data objects threshold ε, etc. In process of cross-correlation network construction, adjustable parameters include the number of feature vectors k during spectral clustering, the number of correlation networks c, and the correlation threshold ε c , etc.
The settings of these parameters play an extremely important role in the construction of the spatio-temporal correlation network. It has a great impact on the quality of the association query results as well. In order to obtain the optimal correlation network, a feedback mechanism is introduced. We use the artificial marking results and the query results based on the correlation network as input to calculate the evaluation index to measure the quality of the association query. Then, we calculate the offset of the above parameters, and realize the step-by-step feedback from the spatio-temporal cross-correlation network to the self-correlation network. Finally, we obtain the value of each parameter meeting the optimal result of the correlation network. In Section 3, we will take threshold ε as an example, to show how different parameter values will affect the results of the associated network construction.

Spatio-Temporal Correlation Network Update
Taking into account the features of remote sensing data such as huge amount and fast update speed, it is necessary to dynamically update the correlation network to meet the demand for fast correlation retrieval. However, comparing these remote sensing data objects with all the data stored to establish correlation network after feature extraction will lead to a great consumption of computing resources and significant decline in the efficiency of data storage and access, which is contrary to the original intention of establishing the spatio-temporal correlation. Therefore, it is especially important to realize an effective dynamic construction method of remote sensing data correlation.
If there are correlations between o i and o j , and between o j and o k , then o i and o k may be related to each other. Therefore, we can take advantage of the transitivity of spatialtemporal correlation to reduce the amount of calculation in the spatio-temporal correlation network update. The spatio-temporal correlation network update process is as followed: 1.
Encapsulate the features of the remote sensing data object o into query conditions to query eligible candidate remote sensing datasets O ; 2.
Query the correlation network set C , corresponding to all spatio-temporal data objects in O ; 3.
Take an element c from C and calculate the correlation degree between o and c ; 4.
Determine whether there is relationship based on the correlation threshold. If there is a correlation, update the structure of the correlation network accordingly, that is, update the adjacency matrix and the correlation matrix stored; 5.
Repeat steps 3 and 4 until all elements are traversed.
The pseudo-code is given in Algorithm 3. Get a set of STCCNs which o i belongs to as C i ; 5.

Algorithm Implementation
The algorithm implementation of the multi-source remote sensing data association method is given in Algorithm 4. Algorithm 4. The algorithm implementation of the multi-source remote sensing data association method.

Multi-Source Remote Sensing Data Association Method Algorithm
Extract data, spatial and time features to construct FD; 6. Analyze the data structure to construct SD; 18. break 19. case UPDATE: 20. Extract data, spatial and time features to construct FD; 21. Analyze the data structure to construct SD;

Spatio-Temporal Index
Spatio-temporal correlation network is a network with hierarchical, localized, multiscale and other characteristics, defined with the hypergraph model. It can be transformed into a spatial data index structure through compression coding to optimize spatial query capabilities, and perform efficient remote sensing data retrieval. The spatio-temporal index proposed is an index established to meet multi-level, multi-scale, and multi-category joint query requirements. It solves the problem of large-scale remote sensing data from the aspects of feature correlation and spatio-temporal correlation. The spatio-temporal index consists of the multi-dimensional feature index and the association index, described as follows.

•
Multi-dimensional feature index The multi-dimensional feature index table is a secondary index table established on  the basis of data storage table, including several parts such as Rowkey, Column and version, as shown in Figure 8. Rowkey is used to search the ID of remote sensing data by their feature. There is only one Column, which is used to store the ID of the remote sensing data. version is used to identify different index versions by their corresponding creation time.
dex consists of the multi-dimensional feature index and the association index, described as follows.

•
Multi-dimensional feature index The multi-dimensional feature index table is a secondary index table established on the basis of data storage table, including several parts such as Rowkey, Column and version, as shown in Figure 8. Rowkey is used to search the ID of remote sensing data by their feature. There is only one Column, which is used to store the ID of the remote sensing data. version is used to identify different index versions by their corresponding creation time. The Rowkey in the multi-dimensional feature index table consists of the retrieval type, retrieval value, and the key of the unit in the data table, integrating the spatial, time, and data features of the remote sensing data into the index structure to ensure that similar data in different feature spaces can be stored adjacently. More detailed descriptions are as follows.
• "SPAC_geocode1_key1" means that the value of the key corresponding to the grid whose spatial feature coding value is "geocode1" is key1; • "TIME_yyyyMMdd1_key1" means that the value of the key corresponding to the grid whose time attribute value is "yyyyMMdd1" is key1; • "FEAT_name1_value1_key1" means that the value of the data feature name1 is the grid data of value1, and the value of the key stored in the data table is key1; • "SPAC", "TIME" and "FEAT" are used to distinguish spatial, time and data features.
The multi-dimensional feature index table contains column family Key and column Key, which are used to store the key of the grid data in the data table. The version of the grid data in the index table is the data production time.

Association index
The association index table is an index table established on the basis of spatio-temporal correlation network, as shown in Figure 9. The Rowkey in the multi-dimensional feature index table consists of the retrieval type, retrieval value, and the key of the unit in the data table, integrating the spatial, time, and data features of the remote sensing data into the index structure to ensure that similar data in different feature spaces can be stored adjacently. More detailed descriptions are as follows.
• "SPAC_geocode1_key1" means that the value of the key corresponding to the grid whose spatial feature coding value is "geocode1" is key1; • "TIME_yyyyMMdd1_key1" means that the value of the key corresponding to the grid whose time attribute value is "yyyyMMdd1" is key1; • "FEAT_name1_value1_key1" means that the value of the data feature name1 is the grid data of value1, and the value of the key stored in the data table is key1; • "SPAC", "TIME" and "FEAT" are used to distinguish spatial, time and data features.
The multi-dimensional feature index table contains column family Key and column Key, which are used to store the key of the grid data in the data table. The version of the grid data in the index table is the data production time.

• Association index
The association index table is an index table established on the basis of spatio-temporal correlation network, as shown in Figure 9.
Remote Sens. 2021, 13, x FOR PEER REVIEW 19 of 28 Figure 9. Association index on remote sensing data.
As shown in Figure 9, the Rowkey in the association index consists of several parts, unique identification of correlation network and unique identification of remote sensing data, etc.
• "DO_dataID1_stcnID1" indicates that the remote sensing data dataID1 belongs to the spatio-temporal correlation network stcnID1; • "CN_stcnID1_dataID1" means that the spatio-temporal correlation network stcnID1 contains the remote sensing data dataID1.
The association index table contains column family named Key, and column named Key, which is used to store the key of the grid data in the data table. The version of the grid data in the index table is the data production time.

Local Association Query
Based on the multi-dimensional feature index and association index, the spatio-tem- As shown in Figure 9, the Rowkey in the association index consists of several parts, unique identification of correlation network and unique identification of remote sensing data, etc.
• "DO_dataID1_stcnID1" indicates that the remote sensing data dataID1 belongs to the spatio-temporal correlation network stcnID1; • "CN_stcnID1_dataID1" means that the spatio-temporal correlation network stcnID1 contains the remote sensing data dataID1.
The association index table contains column family named Key, and column named Key, which is used to store the key of the grid data in the data table. The version of the grid data in the index table is the data production time.

Local Association Query
Based on the multi-dimensional feature index and association index, the spatiotemporal correlation network can be used to realize the local association query of remote sensing data, as shown in Figure 10. The detailed steps are as follows.

1.
Extract the eligible candidate remote sensing dataset O from the database according to query conditions; 2.
Query that the correlation network set C consists of all the correlation network spatio-temporal data objects in O belongs to; 3.
Encapsulate the query conditions as virtual spatio-temporal data objects o; 4.
Calculate the correlation degree between o and the element c l in set C , in turn; 5.
Put the spatio-temporal data objects in c l into the set of spatio-temporal data objects O.

6.
Remove the remote sensing data in O whose correlation degrees are less than ε, meaning that they do not meet the requirement, and sort spatio-temporal data objects by correlation degree.
As shown in Figure 9, the Rowkey in the association index consists of several parts, unique identification of correlation network and unique identification of remote sensing data, etc.
• "DO_dataID1_stcnID1" indicates that the remote sensing data dataID1 belongs to the spatio-temporal correlation network stcnID1; • "CN_stcnID1_dataID1" means that the spatio-temporal correlation network stcnID1 contains the remote sensing data dataID1.
The association index table contains column family named Key, and column named Key, which is used to store the key of the grid data in the data table. The version of the grid data in the index table is the data production time.

Local Association Query
Based on the multi-dimensional feature index and association index, the spatio-temporal correlation network can be used to realize the local association query of remote sensing data, as shown in Figure 10. The detailed steps are as follows.
1. Extract the eligible candidate remote sensing dataset O from the database according to query conditions; 2. Query that the correlation network set C  consists of all the correlation network spatio-temporal data objects in O belongs to; 3. Encapsulate the query conditions as virtual spatio-temporal data objects o ; 4. Calculate the correlation degree between o and the element l c in set C  , in turn; 5. Put the spatio-temporal data objects in l c into the set of spatio-temporal data objects O .
6. Remove the remote sensing data in O whose correlation degrees are less than  , meaning that they do not meet the requirement, and sort spatio-temporal data objects by correlation degree. The algorithm realization of local association query process on remote sensing data is given in Algorithm 5. The algorithm realization of local association query process on remote sensing data is given in Algorithm 5. c

Experimental Materials
We take remote sensing data products from the following sources as test data to effectiveness and versatility of STLAQ's applications on multi-source remote sensing data.

•
Remote sensing satellite image data The remote sensing satellite image data used in the experiments come from Microsoft's Bing map service. Data in the regions of the Hawaiian Islands and Qinghai-Tibet Plateau lakes are selected as the test dataset, in PNG format. There are 256 × 256 pixels in each image, of which the spatial resolution is 10 m.

•
Three-dimensional terrain data Three-dimensional terrain data come from the Google Earth platform. Data in the region of Qinghai-Tibet Plateau lakes are selected as the test dataset, in PNG format. In each image, there are 256 × 256 pixels to present the spatial coverage, of which the spatial resolution is 1.3 m.

•
Street view image data The street view data come from the Baidu panoramic map platform, consisting of parts of the Qinghai-Tibet Plateau, in PNG format.

•
Lake distribution data of the Qinghai-Tibet Plateau The distribution data of lakes on the Qinghai-Tibet Plateau come from the National Qinghai-Tibet Plateau Science Data Center, such as China lake dataset (1960-2020), in Shape File format.

The STCN Construction Results
We encapsulate remote sensing satellite image data into STDUs, marked with rectangles in Figure 11. Then, the STCN model proposed in this paper takes STDUs at different resolution scales as the smallest unit of remote sensing data on the basis of multi-scale division, and establishes the relationship among STDOs based on the STDU, so as to realize different remote sensing data correlations across different resolution scales.

The STCN Construction Results
We encapsulate remote sensing satellite image data into STDUs, marked with rectangles in Figure 11. Then, the STCN model proposed in this paper takes STDUs at different resolution scales as the smallest unit of remote sensing data on the basis of multi-scale division, and establishes the relationship among STDOs based on the STDU, so as to realize different remote sensing data correlations across different resolution scales. The construction results of self-associated networks and cross-associated networks are given in Figures 12 and 13, respectively. In Figure 12b,d, we use rectangles to mark data units. In particular, in Figure 12d, white rectangles and red rectangles are used to represent spatio-temporal data units of two different resolution levels, level 12 and level 13, respectively. It can be seen that in association query of remote sensing data, correlations can be established within remote sensing data of the same scale and of different scales. Therefore, the STCN model proposed in this paper is effective in cross-scale association query, which can be used to establish the relationship between remote sensing data at different resolutions. Thus, the STCN can be used to solve the problems caused by the different shapes of feature elements, and to overcome difficulties in expressing data characteristics uniformly. The settings of these parameters play an extremely important role in the construction of the spatio-temporal correlation network. We take  , the feature weight of mean, and  , the threshold of correlation, as examples, to show how they will affect the results of the correlation network construction, as given in Figure 14. Four indicators, recall rate, The construction results of self-associated networks and cross-associated networks are given in Figures 12 and 13, respectively. In Figure 12b,d, we use rectangles to mark data units. In particular, in Figure 12d, white rectangles and red rectangles are used to represent spatio-temporal data units of two different resolution levels, level 12 and level 13, respectively. It can be seen that in association query of remote sensing data, correlations can be established within remote sensing data of the same scale and of different scales. Therefore, the STCN model proposed in this paper is effective in cross-scale association query, which can be used to establish the relationship between remote sensing data at different resolutions. Thus, the STCN can be used to solve the problems caused by the different shapes of feature elements, and to overcome difficulties in expressing data characteristics uniformly.

The STCN Construction Results
We encapsulate remote sensing satellite image data into STDUs, marked with rectangles in Figure 11. Then, the STCN model proposed in this paper takes STDUs at different resolution scales as the smallest unit of remote sensing data on the basis of multi-scale division, and establishes the relationship among STDOs based on the STDU, so as to realize different remote sensing data correlations across different resolution scales. The construction results of self-associated networks and cross-associated networks are given in Figures 12 and 13, respectively. In Figure 12b,d, we use rectangles to mark data units. In particular, in Figure 12d, white rectangles and red rectangles are used to represent spatio-temporal data units of two different resolution levels, level 12 and level 13, respectively. It can be seen that in association query of remote sensing data, correlations can be established within remote sensing data of the same scale and of different scales. Therefore, the STCN model proposed in this paper is effective in cross-scale association query, which can be used to establish the relationship between remote sensing data at different resolutions. Thus, the STCN can be used to solve the problems caused by the different shapes of feature elements, and to overcome difficulties in expressing data characteristics uniformly. The settings of these parameters play an extremely important role in the construction of the spatio-temporal correlation network. We take  , the feature weight of mean, and  , the threshold of correlation, as examples, to show how they will affect the results of the correlation network construction, as given in Figure 14. Four indicators, recall rate,

The STCN Construction Results
We encapsulate remote sensing satellite image data into STDUs, marked with rectangles in Figure 11. Then, the STCN model proposed in this paper takes STDUs at different resolution scales as the smallest unit of remote sensing data on the basis of multi-scale division, and establishes the relationship among STDOs based on the STDU, so as to realize different remote sensing data correlations across different resolution scales. The construction results of self-associated networks and cross-associated networks are given in Figures 12 and 13, respectively. In Figure 12b,d, we use rectangles to mark data units. In particular, in Figure 12d, white rectangles and red rectangles are used to represent spatio-temporal data units of two different resolution levels, level 12 and level 13, respectively. It can be seen that in association query of remote sensing data, correlations can be established within remote sensing data of the same scale and of different scales. Therefore, the STCN model proposed in this paper is effective in cross-scale association query, which can be used to establish the relationship between remote sensing data at different resolutions. Thus, the STCN can be used to solve the problems caused by the different shapes of feature elements, and to overcome difficulties in expressing data characteristics uniformly. The settings of these parameters play an extremely important role in the construction of the spatio-temporal correlation network. We take  , the feature weight of mean, and  , the threshold of correlation, as examples, to show how they will affect the results of the correlation network construction, as given in Figure 14. Four indicators, recall rate, The settings of these parameters play an extremely important role in the construction of the spatio-temporal correlation network. We take α, the feature weight of mean, and ε, the threshold of correlation, as examples, to show how they will affect the results of the correlation network construction, as given in Figure 14. Four indicators, recall rate, precision rate, F1 rate and Kappa coefficient, are used to evaluate the clustering results of remote sensing data.  As shown in Figure 14a, a larger value of  always leads to lower recall rate, higher precision, higher F1, and higher Kappa coefficient. In addition, as shown in Figure 14b, if the value of  is too small, it will cause redundant data, and if the value of  is too big, it will cause greater incompleteness.

Association Query Evaluation
We take remote sensing data products introduced in Section 4.1 as test datasets, and we take "lake images in a certain geographic range" as the query target to perform local association query tests on multiple types of remote sensing datasets, as shown in Table 2. Considering the convenience of reading, the satellite image data is expressed in STCNs, and rectangular boxes are used to identify STDUs in the self-correlation network.

. Association Query Performance
We use Apache JMeter as the performance testing tool to continuously access the data query interface, to test the retrieval efficiency of three different index structures. The index structures used in the test are as follows.
• Perform remote sensing data retrieval based on the STLAQ proposed in this article; As shown in Figure 14a, a larger value of α always leads to lower recall rate, higher precision, higher F1, and higher Kappa coefficient. In addition, as shown in Figure 14b, if the value of ε is too small, it will cause redundant data, and if the value of ε is too big, it will cause greater incompleteness.

Association Query Evaluation
We take remote sensing data products introduced in Section 4.1 as test datasets, and we take "lake images in a certain geographic range" as the query target to perform local association query tests on multiple types of remote sensing datasets, as shown in Table 2. Considering the convenience of reading, the satellite image data is expressed in STCNs, and rectangular boxes are used to identify STDUs in the self-correlation network.                    We use Apache JMeter as the performance testing tool to continuously access the data query interface, to test the retrieval efficiency of three different index structures. The index structures used in the test are as follows.  Perform remote sensing data retrieval based on the STLAQ proposed in this article;  Perform remote sensing data retrieval based on Geomesa, a distributed architecture for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing data in time and event dimensions.  Perform remote sensing data retrieval based on a customizable multi-dimensional index structure named ST+, proposed in [22].  We use Apache JMeter as the performance testing tool to continuously access the data query interface, to test the retrieval efficiency of three different index structures. The index structures used in the test are as follows.  Perform remote sensing data retrieval based on the STLAQ proposed in this article;  Perform remote sensing data retrieval based on Geomesa, a distributed architecture for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing  We use Apache JMeter as the performance testing tool to continuously access the data query interface, to test the retrieval efficiency of three different index structures. The index structures used in the test are as follows.  Perform remote sensing data retrieval based on the STLAQ proposed in this article;  Perform remote sensing data retrieval based on Geomesa, a distributed architecture for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing  We use Apache JMeter as the performance testing tool to continuously access the data query interface, to test the retrieval efficiency of three different index structures. The index structures used in the test are as follows.  Perform remote sensing data retrieval based on the STLAQ proposed in this article;  Perform remote sensing data retrieval based on Geomesa, a distributed architecture for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing  We use Apache JMeter as the performance testing tool to continuously access the data query interface, to test the retrieval efficiency of three different index structures. The index structures used in the test are as follows.  Perform remote sensing data retrieval based on the STLAQ proposed in this article;  Perform remote sensing data retrieval based on Geomesa, a distributed architecture for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing   We use Apache JMeter as the performance testing tool to continuously access the data query interface, to test the retrieval efficiency of three different index structures. The index structures used in the test are as follows.  Perform remote sensing data retrieval based on the STLAQ proposed in this article;  Perform remote sensing data retrieval based on Geomesa, a distributed architecture for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing  We use Apache JMeter as the performance testing tool to continuously access the data query interface, to test the retrieval efficiency of three different index structures. The index structures used in the test are as follows.  Perform remote sensing data retrieval based on the STLAQ proposed in this article;  Perform remote sensing data retrieval based on Geomesa, a distributed architecture for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing data in time and event dimensions.  Perform remote sensing data retrieval based on a customizable multi-dimensional index structure named ST+, proposed in [22].  We use Apache JMeter as the performance testing tool to continuously access the data query interface, to test the retrieval efficiency of three different index structures. The index structures used in the test are as follows.  Perform remote sensing data retrieval based on the STLAQ proposed in this article;  Perform remote sensing data retrieval based on Geomesa, a distributed architecture for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing data in time and event dimensions.  Perform remote sensing data retrieval based on a customizable multi-dimensional index structure named ST+, proposed in [22].  We use Apache JMeter as the performance testing tool to continuously access the data query interface, to test the retrieval efficiency of three different index structures. The index structures used in the test are as follows.  Perform remote sensing data retrieval based on the STLAQ proposed in this article;  Perform remote sensing data retrieval based on Geomesa, a distributed architecture for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing data in time and event dimensions.  Perform remote sensing data retrieval based on a customizable multi-dimensional index structure named ST+, proposed in [22].  We use Apache JMeter as the performance testing tool to continuously access the data query interface, to test the retrieval efficiency of three different index structures. The index structures used in the test are as follows.  Perform remote sensing data retrieval based on the STLAQ proposed in this article;  Perform remote sensing data retrieval based on Geomesa, a distributed architecture for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing data in time and event dimensions.  Perform remote sensing data retrieval based on a customizable multi-dimensional index structure named ST+, proposed in [22].  We use Apache JMeter as the performance testing tool to continuously access the data query interface, to test the retrieval efficiency of three different index structures. The index structures used in the test are as follows.  Perform remote sensing data retrieval based on the STLAQ proposed in this article;  Perform remote sensing data retrieval based on Geomesa, a distributed architecture for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing data in time and event dimensions.  Perform remote sensing data retrieval based on a customizable multi-dimensional index structure named ST+, proposed in [22].  We use Apache JMeter as the performance testing tool to continuously access the data query interface, to test the retrieval efficiency of three different index structures. The index structures used in the test are as follows.  Perform remote sensing data retrieval based on the STLAQ proposed in this article;  Perform remote sensing data retrieval based on Geomesa, a distributed architecture for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing data in time and event dimensions.  Perform remote sensing data retrieval based on a customizable multi-dimensional index structure named ST+, proposed in [22].  We use Apache JMeter as the performance testing tool to continuously access the data query interface, to test the retrieval efficiency of three different index structures. The index structures used in the test are as follows.  Perform remote sensing data retrieval based on the STLAQ proposed in this article;  Perform remote sensing data retrieval based on Geomesa, a distributed architecture for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing data in time and event dimensions.  Perform remote sensing data retrieval based on a customizable multi-dimensional index structure named ST+, proposed in [22].  We use Apache JMeter as the performance testing tool to continuously access the data query interface, to test the retrieval efficiency of three different index structures. The index structures used in the test are as follows.  Perform remote sensing data retrieval based on the STLAQ proposed in this article;  Perform remote sensing data retrieval based on Geomesa, a distributed architecture for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing data in time and event dimensions.  Perform remote sensing data retrieval based on a customizable multi-dimensional index structure named ST+, proposed in [22].  We use Apache JMeter as the performance testing tool to continuously access the data query interface, to test the retrieval efficiency of three different index structures. The index structures used in the test are as follows.  Perform remote sensing data retrieval based on the STLAQ proposed in this article;  Perform remote sensing data retrieval based on Geomesa, a distributed architecture for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing data in time and event dimensions.  Perform remote sensing data retrieval based on a customizable multi-dimensional index structure named ST+, proposed in [22].

Association Query Performance
We use Apache JMeter as the performance testing tool to continuously access the data query interface, to test the retrieval efficiency of three different index structures. The index structures used in the test are as follows.
• Perform remote sensing data retrieval based on the STLAQ proposed in this article; • Perform remote sensing data retrieval based on Geomesa, a distributed architecture for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing data in time and event dimensions.
• Perform remote sensing data retrieval based on a customizable multi-dimensional index structure named ST+, proposed in [22].
In order to verify the retrieval efficiency of the STLAQ, we choose the remote sensing image data of layers 0-9 in the world as the test dataset, a total of 690,000. The experiment process is as follows.

1.
We perform spatial retrieval, time retrieval, and feature retrieval on tables that store different amounts of data, respectively. 2.
Then, we use query windows of different sizes to carry out retrieval in multiple dimensions such as space, time, and feature, on the table with the largest amount of stored data. 3.
Finally, we adopt two ways to query the combination of two dimensions, space-time and space-feature, to test the efficiency of STLAQ in multi-dimensional combination query.
The results of local association query on spatial, time and feature dimension are given in Figures 15-17 • Perform remote sensing data retrieval based on Geomesa, a distributed architecture for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing data in time and event dimensions. • Perform remote sensing data retrieval based on a customizable multi-dimensional index structure named ST+, proposed in [22].
In order to verify the retrieval efficiency of the STLAQ, we choose the remote sensing image data of layers 0-9 in the world as the test dataset, a total of 690,000. The experiment process is as follows.
1. We perform spatial retrieval, time retrieval, and feature retrieval on tables that store different amounts of data, respectively. 2. Then, we use query windows of different sizes to carry out retrieval in multiple dimensions such as space, time, and feature, on the table with the largest amount of stored data. 3. Finally, we adopt two ways to query the combination of two dimensions, space-time and space-feature, to test the efficiency of STLAQ in multi-dimensional combination query.
The results of local association query on spatial, time and feature dimension are given in Figures 15-17, respectively. It can be seen from Figure 15 that the STLAQ have better query performance on spatial association query, compared with Geomesa and ST+. As the amount of remote sensing data grows, the ST+ and Geomesa show a significant decline in query efficiency, while the STLAQ is relatively stable. With the expansion of the spatial query window, the STLAQ shows the slowest decline in efficiency.  for spatio-temporal fusion, proposed in [18]. We built indexes on the remote sensing data in time and event dimensions. • Perform remote sensing data retrieval based on a customizable multi-dimensional index structure named ST+, proposed in [22].
In order to verify the retrieval efficiency of the STLAQ, we choose the remote sensing image data of layers 0-9 in the world as the test dataset, a total of 690,000. The experiment process is as follows.
1. We perform spatial retrieval, time retrieval, and feature retrieval on tables that store different amounts of data, respectively. 2. Then, we use query windows of different sizes to carry out retrieval in multiple dimensions such as space, time, and feature, on the table with the largest amount of stored data. 3. Finally, we adopt two ways to query the combination of two dimensions, space-time and space-feature, to test the efficiency of STLAQ in multi-dimensional combination query.
The results of local association query on spatial, time and feature dimension are given in Figures 15-17, respectively. It can be seen from Figure 15 that the STLAQ have better query performance on spatial association query, compared with Geomesa and ST+. As the amount of remote sensing data grows, the ST+ and Geomesa show a significant decline in query efficiency, while the STLAQ is relatively stable. With the expansion of the spatial query window, the STLAQ shows the slowest decline in efficiency.  relatively stable. The STLAQ has the highest query efficiency, which slightly decreases with the expansion of the query window. As shown in Figure 17, when the feature query window is small, STLAQ has the highest query efficiency. The increase of the window size leads to a significant downward trend in query efficiency of the three indexes. Among them, the Geomesa index has the slowest decline, STLAQ's is second, and the decline of ST+ is the most severe. The results of association query on space-time and space-feature dimensions are given in Figures 18 and 19. In the experiments, the data of layer 9 is queried by default.  It can be seen from Figure 15 that the STLAQ have better query performance on spatial association query, compared with Geomesa and ST+. As the amount of remote sensing data grows, the ST+ and Geomesa show a significant decline in query efficiency, while the STLAQ is relatively stable. With the expansion of the spatial query window, the STLAQ shows the slowest decline in efficiency.
In general, the STLAQ performs better in remote sensing data association queries on time dimension, followed by Geomesa, and ST+ performs poorly, according to the result given in Figure 16. The query efficiency of the ST+ index drops sharply as the amount of stored data increases, compared to the other two indexes. As for queries in time windows of different sizes, the ST+ index fluctuates to a certain extent, while the Geomesa index is relatively stable. The STLAQ has the highest query efficiency, which slightly decreases with the expansion of the query window.
As shown in Figure 17, when the feature query window is small, STLAQ has the highest query efficiency. The increase of the window size leads to a significant downward trend in query efficiency of the three indexes. Among them, the Geomesa index has the slowest decline, STLAQ's is second, and the decline of ST+ is the most severe.
The results of association query on space-time and space-feature dimensions are given in Figure 18. In the experiments, the data of layer 9 is queried by default. In general, the STLAQ performs better in remote sensing data association queries on time dimension, followed by Geomesa, and ST+ performs poorly, according to the result given in Figure 16. The query efficiency of the ST+ index drops sharply as the amount of stored data increases, compared to the other two indexes. As for queries in time windows of different sizes, the ST+ index fluctuates to a certain extent, while the Geomesa index is relatively stable. The STLAQ has the highest query efficiency, which slightly decreases with the expansion of the query window.
As shown in Figure 17, when the feature query window is small, STLAQ has the highest query efficiency. The increase of the window size leads to a significant downward trend in query efficiency of the three indexes. Among them, the Geomesa index has the slowest decline, STLAQ's is second, and the decline of ST+ is the most severe. The results of association query on space-time and space-feature dimensions are given in Figures 18 and 19. In the experiments, the data of layer 9 is queried by default.  According to the results given in Figure 18, STLAQ is more sensitive to the size of the time query window. Under different spatial query ranges, it always shows that as the time query range increases, the query efficiency slightly decreases; while the other two indexes remain relatively stable. As the range of spatial queries increases, the gap between the STLAQ and Geomesa indexes and the ST+ index gradually increases.
It can be seen from Figure 19 that, although the spatial query ranges are different, the query performance of STLAQ and Geomesa always decreases gradually as the time query range increases. ST+ has no advantage in feature query, and the efficiency gap between it and the other two indexes gradually widens with the gradual expansion of the query window. According to the results given in Figure 18, STLAQ is more sensitive to the size of the time query window. Under different spatial query ranges, it always shows that as the time query range increases, the query efficiency slightly decreases; while the other two indexes remain relatively stable. As the range of spatial queries increases, the gap between the STLAQ and Geomesa indexes and the ST+ index gradually increases.
It can be seen from Figure 19 that, although the spatial query ranges are different, the query performance of STLAQ and Geomesa always decreases gradually as the time query range increases. ST+ has no advantage in feature query, and the efficiency gap between it and the other two indexes gradually widens with the gradual expansion of the query window.
According to the above results, the proposed STLAQ has good efficiency in association queries of massive remote sensing data in spatial, temporal, and feature dimension, especially when querying in a combination of a large spatial range and a small time and feature range. There are two main reasons: • STLAQ has established indexes on both the features and association, and maps multidimensional features to a one-dimensional space. It avoids the drawbacks of a sharp increase in the scanning range during data retrieval as the data dimension increases; • When constructing indexes on the time and feature dimensions, the sequential coding method ensures the locality of the data. This also brings about the insufficient matching accuracy when the scope of the query on these two dimensions expands, resulting in a decrease in retrieval performance.
In the future, we will optimize from two aspects.
• Firstly, under the premise of ensuring the continuity of time and attributing characteristics, we should optimize the structure of the index, improve the process of feature encoding, and achieve a more accurate match between the query range and the feature index; • Secondly, we should introduce a parallel retrieval mechanism to improve data at the system level search efficiency.
Generally speaking, the query efficiency is usually contradictory to the resource occupation. The improvement of query efficiency often leads to more consumption of resources such as CPU, memory, disk reads and writes and other resources. In order to verify the performance of STLAQ in terms of resource occupancy, we take spatial association query and time association query as examples, to compare the resource consumption of STLAQ and other indexes. The results are given in Figure 20. According to the above results, the proposed STLAQ has good efficiency in association queries of massive remote sensing data in spatial, temporal, and feature dimension, especially when querying in a combination of a large spatial range and a small time and feature range. There are two main reasons: • STLAQ has established indexes on both the features and association, and maps multidimensional features to a one-dimensional space. It avoids the drawbacks of a sharp increase in the scanning range during data retrieval as the data dimension increases; • When constructing indexes on the time and feature dimensions, the sequential coding method ensures the locality of the data. This also brings about the insufficient matching accuracy when the scope of the query on these two dimensions expands, resulting in a decrease in retrieval performance.
In the future, we will optimize from two aspects.
• Firstly, under the premise of ensuring the continuity of time and attributing characteristics, we should optimize the structure of the index, improve the process of feature encoding, and achieve a more accurate match between the query range and the feature index; • Secondly, we should introduce a parallel retrieval mechanism to improve data at the system level search efficiency.  In Figure 20, ST+ has the lowest resource occupancy rate, followed by STLAQ. Geomesa, which performs well in the previous test, consumes a lot of resources such as CPU, disk read and write. There is a slow increase in the value of resource occupancy metrics with data quantity unceasing increase. Therefore, it can be considered that although the STLAQ has a large contribution to query efficiency, it will not carry out more consumption to resources.
To sum up, the STLAQ proposed in this paper can break through the limitation of resolution scale. It establishes a cross-scale correlation between multi-source remote sensing data with spatial, time, and date features. Moreover, it makes the resource consumption of high-precision data query tend to be flat. The STLAQ case expand data sharing mode, and has high practical significance in the application of remote sensing data.

Conclusions
In view of the current issues of remote sensing big data management and retrieval, we propose a spatio-temporal local association query algorithm, STLAQ. We design a spatio-temporal data model and a bottom-up spatio-temporal correlation network. We adopt the method of the FCM clustering and the method of spectral clustering, to measure the correlation between spatio-temporal correlation networks. Then, we design the spatiotemporal index and local association query method based on spatio-temporal correlation networks, to provide multi-dimensional association query capabilities for remote sensing big data. We construct the STCN and carry out query experiments on several typical remote sensing datasets. The results show that the STLAQ can effectively establish crossscale correlations of remote sensing data, and can provide capabilities of joint retrieval in multiple dimensions, such as spatial features, time features and data features. The STLAQ provide a solution to deal with the problems of implicit correlation and fast service under large-scale remote sensing data, especially with high spatial resolution in various fields. It expands the sharing mode of remote sensing data.
Considering that multi-source remote sensing data have broad application prospects, it is necessary to improve the adaptive ability of the algorithm in the later stage. Future studies based on neural network learning will be performed, to speed up the convergence of network weights, reducing the influence of experience. With the support of parallel computing technology, a distributed framework of data index and analysis would be provided to enhance the practicality of the algorithm. On the other hand, the dynamic expansion and maintenance method of the index during spatio-temporal insertion is also one of the main contents of future research.  In Figure 20, ST+ has the lowest resource occupancy rate, followed by STLAQ. Geomesa, which performs well in the previous test, consumes a lot of resources such as CPU, disk read and write. There is a slow increase in the value of resource occupancy metrics with data quantity unceasing increase. Therefore, it can be considered that although the STLAQ has a large contribution to query efficiency, it will not carry out more consumption to resources.
To sum up, the STLAQ proposed in this paper can break through the limitation of resolution scale. It establishes a cross-scale correlation between multi-source remote sensing data with spatial, time, and date features. Moreover, it makes the resource consumption of high-precision data query tend to be flat. The STLAQ case expand data sharing mode, and has high practical significance in the application of remote sensing data.

Conclusions
In view of the current issues of remote sensing big data management and retrieval, we propose a spatio-temporal local association query algorithm, STLAQ. We design a spatio-temporal data model and a bottom-up spatio-temporal correlation network. We adopt the method of the FCM clustering and the method of spectral clustering, to measure the correlation between spatio-temporal correlation networks. Then, we design the spatiotemporal index and local association query method based on spatio-temporal correlation networks, to provide multi-dimensional association query capabilities for remote sensing big data. We construct the STCN and carry out query experiments on several typical remote sensing datasets. The results show that the STLAQ can effectively establish crossscale correlations of remote sensing data, and can provide capabilities of joint retrieval in multiple dimensions, such as spatial features, time features and data features. The STLAQ provide a solution to deal with the problems of implicit correlation and fast service under large-scale remote sensing data, especially with high spatial resolution in various fields. It expands the sharing mode of remote sensing data.
Considering that multi-source remote sensing data have broad application prospects, it is necessary to improve the adaptive ability of the algorithm in the later stage. Future studies based on neural network learning will be performed, to speed up the convergence of network weights, reducing the influence of experience. With the support of parallel computing technology, a distributed framework of data index and analysis would be provided to enhance the practicality of the algorithm. On the other hand, the dynamic expansion and maintenance method of the index during spatio-temporal insertion is also one of the main contents of future research.
Author Contributions: Y.H. proposed the concept; L.Z. and X.S. designed the STLAQ, made the first prototype of the STLAQ and wrote the original draft; X.T. and K.F. reviewed and edited the manuscript; Y.H. provided the project. All the authors contributed equally to the revisions of the manuscript. All authors have read and agreed to the published version of the manuscript.