CN109992636A

CN109992636A - Space-time code method, temporal index and querying method and device

Info

Publication number: CN109992636A
Application number: CN201910221266.1A
Authority: CN
Inventors: 童晓冲; 雷毅; 程承旗; 吴翔宇; 赖广陵; 李贺; 郭从洲; 张勇; 陈波
Original assignee: Information Engineering University of PLA Strategic Support Force
Current assignee: Information Engineering University of PLA Strategic Support Force
Priority date: 2019-03-22
Filing date: 2019-03-22
Publication date: 2019-07-09
Anticipated expiration: 2039-03-22
Also published as: CN109992636B

Abstract

This specification one or more embodiment discloses a kind of space-time code method, temporal index and querying method and device, to realize the high-efficiency tissue to massive multisource remote sensing data, management and inquiry.The described method includes: carrying out subdivision according to longitude and latitude space of the specified partition patterns to the earth, multiple grids are obtained；Wherein, longitudinal is identical with corresponding grid number on latitude direction；The specified partition patterns include specified subdivision level；Integer coding is carried out to the corresponding coordinate information of each grid, obtains corresponding first encoded radio of each grid；Establish the incidence relation between the grid and corresponding first encoded radio of each grid and the remote sensing image data.The technical solution can be associated with to realize with the position of multi- source Remote Sensing Data data using the multiple dimensioned grid integer coding in the whole world, and realize the temporal index and inquiry of multi- source Remote Sensing Data data.

Description

Space-time code method, temporal index and querying method and device

Technical field

This specification is related to technical field of data processing more particularly to a kind of space-time code method, temporal index and inquiry Method and device.

Background technique

Global earth observation technology, system and ability are increasingly complete, form more means, multi-platform stereopsis body The novel observation such as system and developing PNRTC (positioning, time service, remote sensing, communicates integrated space-based information service system at navigation) System.In recent years microsatellite is in the fast development stage, and if dove groups of stars' seat is up to more than 200, earth observation data resource is very big It is abundant.Remotely-sensed data shows polyphyly, the features such as high spatial, time, spectral resolution and Global coverage, data volume at present Amount forms the development trend of PB grades, EB grades.High-efficiency tissue, management and the inquiry for how realizing massive multisource remote sensing data, have become For one of related service mechanism and application department's urgent problem to be solved.

Existing remotely-sensed data storage management system or center are broadly divided into following two categories according to Method of Data Organization:

(1) it is based on multi-resolution pyramid image tiled manner

Had using the canonical system of which: Word Wind, TerraServer, Bing Maps, Google Maps/ Earth, day map etc..As shown in Figure 1, such mode be mainly utilize series of rules, it is seamless, have multiple dimensioned level The grid tile of structure completely continuously covers earth surface space, and the Spatial Cable of image is established by the station location marker of tile Draw.Which is mainly used in seamless organization and the visualization (photomap) of remotely-sensed data, allow users to rapidly obtain, The characteristics remotely-sensed datas such as multi-source, the multidate inquiring the satellite image data in some region and showing, but lack the same area Unified management.Further, since image subdivision and pyramid construction, the remotely-sensed data organizational form based on global subdivision grid will A large amount of tile small documents are generated, about 1/3 or so data increment is caused, becomes larger so as to cause data space.If It is stored using distributed file system, it will generate a large amount of mapped file and journal file, certain nodes is caused single-point occur Failure is unfavorable for the storage and management of distributed file system.

(2) based on when null record system satellite orbit band or scape mode

Had using the canonical system of which: NASAEOS, data center, European Space Agency, remote sensing fairground, Chinese Resources satellite are answered With center, State Satellite Meterological Center, National Satellite Ocean Application Service etc..As shown in Fig. 2, such mode presses image data According to band or divide scape imaged tissue, metadata is established frequently with business database management using the spatial information of image Its spatial index is mainly used for managing received raw video data and data product.Which is easy to operate, and remote sensing is defended Star negligible amounts, the lesser data organization and management of the data scale of construction can satisfy the management and application demand for achieving data.So And as remotely-sensed data source category increases, since the track strip band between each remotely-sensed data production unit lacks unification with scape Segmentation standard and position it is corresponding, product data mark learns connotation, the multi-source of areal, multiple dimensioned, multidate number with lacking Lack space scale between to be associated with position, causes the different data product relevance of the same area poor；The same area it is more Source remotely-sensed data is often recorded in different track strip bands, to multi-source, the multidate of one specific region of trans-departmental integration Data, it is very time-consuming, to bring the inconvenience of remotely-sensed data retrieval, management and integration.

As shown in the above, it is based on multi-resolution pyramid image tiled manner, main purpose is solved based on image Real world truly expressed and presentation, lack the system of the characteristics remotely-sensed datas such as the multi-source, multiple dimensioned, multidate of the same area One management.Based on when null record system satellite orbit band or scape mode, the position relevance between multi- source Remote Sensing Data data is poor, And different model satellite often also uses respective Method of Data Organization and recording mode, so that data information island phenomenon is tight Weight leads to the retrieval of the same area remotely-sensed data, integration and shared difficulty, reduces data effective utilization.As it can be seen that both modes It is each has something to recommend him, but the demands such as remotely-sensed data Organization And Management that can not all fully meet magnanimity multi-source instantly.

Summary of the invention

The purpose of this specification one or more embodiment is to provide a kind of space-time code method, temporal index and issuer Method and device, to realize the high-efficiency tissue to massive multisource remote sensing data, management and inquiry.

In order to solve the above technical problems, this specification one or more embodiment is achieved in that

On the one hand, this specification one or more embodiment provides a kind of space-time code method, comprising:

Subdivision is carried out according to longitude and latitude space of the specified partition patterns to the earth, obtains multiple grids；Wherein, longitudinal It is identical with grid number corresponding on latitude direction；The specified partition patterns include specified subdivision level；

Integer coding is carried out to the corresponding coordinate information of each grid, each grid corresponding first is obtained and compiles Code value；

Establish the grid and corresponding first encoded radio of each grid and the remote sensing image data it Between incidence relation.

In one embodiment, described to carry out subdivision according to longitude and latitude space of the specified partition patterns to the earth, it obtains more A grid, comprising:

Latitude space is extended, so that the spatial dimension one of the spatial dimension in the latitude space and longitude spaces It causes, the longitude and latitude space after being expanded；

According to quaternary tree mode, the longitude and latitude space after the extension is divided into the specified subdivision level grid, Obtain the multiple grid；And determine the coordinate origin of the multiple grid.

In one embodiment, it is described establish the grid and corresponding first encoded radio of each grid with Incidence relation between the remote sensing image data, comprising:

Calculate the boundary rectangle of the corresponding imagery zone of the remote sensing image data；

Gridding is carried out to the boundary rectangle using the grid of lowest hierarchical level, is obtained corresponding with the boundary rectangle The first grid；

Judge whether first grid meets preset condition；Wherein, the preset condition includes at least one of the following: institute The sum for stating the first kind grid and the second class grid in the first grid is greater than or equal to preset group/cording quantity upper limit value, described The number of grid of second class grid is zero；

If so, determining that first grid is grid corresponding with the imagery zone；

If it is not, execution following steps are then recycled, until the grid obtained meets the preset condition: for each described Second grid of highest priority in second class grid, the sub-grid of the adjacent level of level where determining second grid； The sub-grid is divided into the first kind grid and/or the second class grid, and updates the excellent of the second class grid First grade.

In one embodiment, the first kind grid includes: the grid for including by the boundary rectangle, or, clathrum Grade is greater than or equal to the grid of preset level upper limit value；

The second class grid include: intersect with the boundary rectangle and trellis stage be less than the level upper limit value Grid, or, the grid comprising the influence area.

In one embodiment, the priority orders of the second class grid are true according at least one in following rule It is fixed:

Clathrum series is smaller, and the priority is higher；

Specified type number of grid in the sub-grid that the grid is included is more, and the priority is higher；The finger Determining type cell includes the grid that region corresponding with the remote sensing image data mutually separates；

First kind number of grid in the sub-grid that the grid is included is fewer, and the priority is higher.

On the other hand, this specification one or more embodiment provides a kind of temporal index and querying method, comprising:

Create the data directory for managing the source data of remote sensing image data；It include following in the data directory At least one field: the source data of the remote sensing image data, Data Identification, data obtaining time coding, data correlation Grid coding set；

According to each field in the data directory, the spatial index for establishing the remote sensing image data is created Grid coding concordance list；It include earth subdivision grid coding and Data Identification collection in the grid coding concordance list；Wherein, described Earth subdivision grid coding is to carry out encoded radio corresponding to each grid obtained after subdivision to the longitude and latitude space of the earth；

According to the data directory and the grid coding concordance list, the remote sensing image data to be checked is carried out Inquiry.

In one embodiment, the data directory created for managing the source data of remote sensing image data, comprising:

Determine the field name and data type of each field in the data directory；

Read the source data corresponding with the field name and the data type, and the source that will be read Data import in the data directory；

Determine each field in the image density look-up table of the whole world；The whole world image density look-up table includes following at least one Item field: the degree of covering of the remote sensing image data in subdivision grid coding, each grid, number of grid upper limit value multiplying power；

It is encoded according to the data obtaining time, establishes the time encoding index of the remote sensing image data.

In one embodiment, the data obtaining time coding in the data directory is true in the following way It is fixed:

The integer coding that specified level is carried out to the temporal information of the remote sensing image data, obtains the remote sensing image number According to the data obtaining time coding.

In one embodiment, the integer that the temporal information to the remote sensing image data carries out specified level is compiled Code obtains the data obtaining time coding of the remote sensing image data, comprising:

The temporal information is decomposed into the integer of multiple specified time scales；

Integer on each specified time scale is separately encoded the binary number to specify digit, obtains the time Binary code value of the information on each specified time scale；

The temporal information is subjected to the connection on bit field in the binary code value on each specified time scale, is obtained Corresponding second encoded radio of the temporal information；

Second encoded radio is moved to left one, the data obtaining time for obtaining the remote sensing image data is compiled Code.

In one embodiment, the specified time scale include year, month, day, hour, minute, second, millisecond, in microsecond At least one of；

Correspondingly, the integer by each specified time scale is separately encoded the binary number to specify digit, It includes at least one of the following:

The binary number for being 17bit by the integer coding in this time scale of the year；

The binary number for being 4bit by the integer coding in this time scale of the moon；

The binary number for being 5bit by the integer coding in this time scale of the day；

The binary number for being 5bit by the integer coding in this time scale of the hour；

The binary number for being 6bit by the integer coding in described point of this time scale；

The binary number for being 6bit by the integer coding in this time scale of the second；

The binary number for being 10bit by the integer coding in this time scale of the millisecond；

The binary number for being 10bit by the integer coding in this time scale of the microsecond.

In one embodiment, the method also includes:

When receiving the specified operational order to the first data in the data directory, first data are held Row corresponding operating；And the time encoding updated in the data directory indexes；

Wherein, the specified operational order includes deleting instruction and/or inserting instruction.

In one embodiment, each field according in the data directory is created for establishing the remote sensing The grid coding concordance list of the spatial index of image data, comprising:

According to the data directory and the global image density look-up table, each number in the data directory is determined According to corresponding data correlation grid coding set；

The data correlation grid coding set is inserted into the grid coding concordance list；

The data correlation grid coding set being inserted into the grid coding concordance list is ranked up, to establish State the one-dimensional space index of grid coding concordance list.

In one embodiment, the method also includes:

When the data in the data directory change, judge whether the data volume of the variation reaches present count According to amount；Wherein, described change includes insertion data and/or deletion data；

If so, re-creating the grid coding concordance list.

In one embodiment, described according to the data directory and the grid coding concordance list, to be checked The remote sensing image data is inquired, comprising:

Determine the query region in the earth subdivision grid；The query region is made of multiple latitude and longitude coordinates；

According to the latitude and longitude coordinates and the earth subdivision grid coding, the latitude and longitude coordinates corresponding first are determined Gridding coding；

The first sub-grid that corresponding each grid is encoded for first gridding, in the grid coding concordance list The data correlation grid coding set for belonging to first sub-grid is searched, as the first query result；

It determines that first gridding encodes father's grid of corresponding each grid, and determines the corresponding earth of father's grid Subdivision grid coding is that the second gridding of the query region encodes；

The second sub-grid that corresponding each grid is encoded for second gridding, in the grid coding concordance list The data correlation grid coding set for belonging to second sub-grid is searched, and screens and meets the described of the first preset condition The corresponding data correlation grid coding set of second sub-grid is as the second query result；The first preset condition packet It includes: there is the gridding coding for belonging to second query result in the first gridding coding；

It is extracted from the grid coding concordance list corresponding with first query result and second query result The first Data Identification collection, obtain the corresponding space querying result of the remote sensing image data to be checked.

Determine the query argument of the remote sensing image data to be checked, the query argument include the query region and Query time range；The query time range includes initial time and termination time；

Time encoding is carried out to the initial time and the termination time respectively, obtains initial time encoded radio and termination Time encoding value；

For the first Data Identification collection, filtered out from the data directory and the first Data Identification collection pair Field that is answering and meeting the second preset condition, as the corresponding spatial-temporal query of the remote sensing image data to be checked As a result；Wherein, second preset condition includes: that the corresponding data obtaining time coding of the first Data Identification collection is located at institute It states initial time encoded radio and terminates between time encoding value.

In another aspect, this specification one or more embodiment provides a kind of space-time code device, comprising:

Subdivision module obtains multiple grids for carrying out subdivision according to longitude and latitude space of the specified partition patterns to the earth； Wherein, longitudinal is identical with corresponding grid number on latitude direction；The specified partition patterns include specified subdivision layer Grade；

Coding module obtains each grid point for carrying out integer coding to the corresponding coordinate information of each grid Not corresponding first encoded radio；

Relating module, for establish the grid and corresponding first encoded radio of each grid with it is described distant Feel the incidence relation between image data.

In another aspect, this specification one or more embodiment provides a kind of temporal index and inquiry unit, comprising:

First creation module, for creating the data directory of the source data for managing remote sensing image data；The number According to including at least one of following field: the source data of the remote sensing image data, Data Identification, data acquisition in concordance list Time encoding, data correlation grid coding set；

Second creation module, for creating for establishing the remote sensing shadow according to each field in the data directory As the grid coding concordance list of the spatial index of data；It include earth subdivision grid coding and number in the grid coding concordance list According to identification sets；Wherein, the earth subdivision grid coding is each grid for obtain after subdivision to the longitude and latitude space of the earth Corresponding encoded radio；

Enquiry module, for according to the data directory and the grid coding concordance list, to be checked described distant Sense image data is inquired.

Using the technical solution of this specification one or more embodiment, cutd open by the longitude and latitude space to the earth Point, multiple grids are obtained, and integer coding is carried out to the corresponding coordinate information of each grid.Obtain the corresponding coding of each grid Value, and then establish the incidence relation between grid and the corresponding encoded radio of each grid and remote sensing image data.As it can be seen that the skill Art scheme can be associated with to realize with the position of multi- source Remote Sensing Data data using the multiple dimensioned grid integer coding in the whole world, and then realized and used One-dimensional integer coding carrys out the spatial positional information of Unify legislation image.

Further, which can be formatted the space-time code method of tissue based on World Wide Web, establish multi-source remote sensing The spatial index of image data, and unify using multiple dimensioned period integer coding method the temporal information of remote sensing image data, With this settling time index, it is achieved that the high-efficiency tissue of massive multisource remote sensing data, management and inquiry.

Detailed description of the invention

In order to illustrate more clearly of this specification one or more embodiment or technical solution in the prior art, below will A brief introduction will be made to the drawings that need to be used in the embodiment or the description of the prior art, it should be apparent that, it is described below Attached drawing is only some embodiments recorded in this specification one or more embodiment, and those of ordinary skill in the art are come It says, without any creative labor, is also possible to obtain other drawings based on these drawings.

Fig. 1 is according to a kind of signal of the tile division mode based on multi-resolution pyramid image in the prior art Figure；

Fig. 2 be according to it is in the prior art it is a kind of based on when null record system satellite orbit strip data tissue model Schematically；

Fig. 3 is the schematic flow chart according to a kind of space-time code method of one embodiment of this specification；

Fig. 4 is according to a kind of the schematic of single scale lattice encoding method of one embodiment of this specification；

Fig. 5 is according to a kind of the schematic of multiple dimensioned lattice encoding method of one embodiment of this specification；

Fig. 6 is the schematic diagram of earth subdivision grid in a kind of space-time code method according to one embodiment of this specification；

Fig. 7 is the corresponding relationship in a kind of space-time code method according to one embodiment of this specification between data and grid Schematic diagram；

Fig. 8 is in a kind of space-time code method according to one embodiment of this specification between grid and image overlay area The schematic diagram of spatial relationship；

Fig. 9 is the schematic flow chart of a kind of temporal index and querying method according to one embodiment of this specification；

Figure 10 is the time encoding conversion side in a kind of temporal index and querying method according to one embodiment of this specification The schematic diagram of method；

Figure 11 is creation grid coding index in a kind of temporal index and querying method according to one embodiment of this specification The schematic diagram of table GICT；

Figure 12 is the schematic block diagram according to a kind of space-time code device of one embodiment of this specification；

Figure 13 is the schematic block diagram of a kind of temporal index and inquiry unit according to one embodiment of this specification.

Specific embodiment

This specification one or more embodiment provides a kind of space-time code method, temporal index and querying method and dress It sets, to realize the high-efficiency tissue to massive multisource remote sensing data, management and inquiry.

In order to make those skilled in the art more fully understand the technical solution in this specification one or more embodiment, Below in conjunction with the attached drawing in this specification one or more embodiment, to the technology in this specification one or more embodiment Scheme is clearly and completely described, it is clear that and described embodiment is only this specification a part of the embodiment, rather than Whole embodiments.Based on this specification one or more embodiment, those of ordinary skill in the art are not making creativeness The model of this specification one or more embodiment protection all should belong in every other embodiment obtained under the premise of labour It encloses.

Fig. 3 be according to a kind of schematic flow chart of space-time code method of one embodiment of this specification, as shown in figure 3, This method comprises:

S302 carries out subdivision according to longitude and latitude space of the specified partition patterns to the earth, obtains multiple grids；Wherein, it passes through It is identical with corresponding grid number on latitude direction to spend direction；Specified partition patterns include specified subdivision level.

S304 carries out integer coding to the corresponding coordinate information of each grid, obtains corresponding first coding of each grid Value.

S306 establishes being associated between grid and corresponding first encoded radio of each grid and remote sensing image data System.

Each step in above-described embodiment is described in detail below.

Firstly, carrying out subdivision according to longitude and latitude space of the specified partition patterns to the earth.

The multiple dimensioned lattice encoding method in the whole world in the present embodiment, replaces conventional two-dimensional coordinate (longitude and latitude using integer coding Degree coordinate) come multiple dimensioned grid is identified, organization and management, in order to realize the coding and height of multiple dimensioned grid The coding of effect calculates and spatial retrieval.Illustrate the integer coding method of two-dimensional space single scale grid, and then first below with this Based on, illustrate the process for establishing the integer coding of the multiple dimensioned grid of two-dimensional space.

When carrying out the integer coding of two-dimensional space single scale grid, entire two-dimensional space is subjected to four forks first and is divided, Two-dimensional space is divided into four identical subspaces, then every sub-spaces are continued to be divided into the sub empty of four higher levels Between, recurrence according to the method, until the subspace of highest level as defined in obtaining.For example, sharing 2 in m grades of indexes^m×2^m A grid, the specific coding rule of every level-one is as shown in Figure 4.

As seen from Figure 4, which is made of one-dimensional Z-shaped curve a series of integers, so that two-dimensional mesh coordinate can It is mapped in the one-dimensional space.Integer coding is carried out to mesh coordinate, needs to establish the corresponding relationship of mesh coordinate and encoded radio.For Computer memory space is sufficiently used, meets and encodes simple and highly efficient requirement, devise mesh coordinate knot as shown in Table 1 Structure.Computer is in the environment of x64, and when not considering the case where memory splices, the maximum integer that can store is 64, gives X, Y Two coordinates respectively distribute 31 (concrete reason will be described hereinafter), then the value range of each coordinate be 0~ 2147483647。

Table 1

Meaning	X	Y
			Digit	31	31
Value range	0~2147483647	0~2147483647
			Storage value	X-direction mesh coordinate	Y-direction mesh coordinate

The transverse and longitudinal coordinate of grid is converted to integer coding using the method for intersecting fetch bit by the characteristics of according to Z-shaped coding.Example Such as, mesh coordinate is that steps are as follows for (20,17) corresponding single scale integer coding calculating: (1) mesh coordinate (20,17) is corresponding Binary coding position (10100,10001)；(2) intersect fetch bit and obtain 1100010010；(3) the single scale integer coding of grid is obtained 786。

Above-mentioned single scale integer coding cannot express multi-scale information simultaneously.Since the foundation of multiple dimensioned integer coding is Based on single scale integer coding, therefore by single scale integer coding value to moving to left one, it can be obtained multiple dimensioned integer and compile The encoded radio of maximum level (the 31st level) in code.If the X, Y coordinates of grid respectively take 32, then reconvert is whole at single scale Number encoder value obtains multiple dimensioned integer coding value since the value can not move to left one, so in single scale integer coding method, The X, Y coordinates of grid respectively take 31, in order to reserve 2 to record dimensional information.

Obviously, the integer coding value in the 31st level is even number, and the encoded radio of other levels is all using this level as base Plinth and generate.It will be averaged in 31st level per 4 adjacent integer coding values, the integer that the 30th level can be obtained is compiled Code value, and be odd number, and so on, as shown in figure 5, totally 31 grades of integer coding value can be obtained, form the quaternary tree of handstand.

According to the characteristics of above-mentioned multiple dimensioned grid integer coding and the characteristic of terrestrial space, can cut open as follows Divide earth grid: firstly, being extended to latitude space, so that the spatial dimension of the spatial dimension in latitude space and longitude spaces Unanimously, the longitude and latitude space after being expanded；Secondly, the longitude and latitude space after extension is divided into finger according to quaternary tree mode Determine subdivision level grid, obtains multiple grids；And determine the coordinate origin of multiple grids.

In the present embodiment, be using longitude and latitude mode to the reason of earth progress subdivision: longitude and latitude mode has conversion The minimum feature of cost, and the coordinate basis of earth subdivision is in longitude and latitude space.The latitude and longitude coordinates system can be extensively Adopted WGS84 coordinate system.

Further, since multiple dimensioned coding is based on the rectangular mesh of quaternary tree, therefore according to two sides of quaternary tree To equal principle, global longitude range (- 180 °, 180 °) when the whole subdivision of (- 90 °, 90 °) progress of latitude scope, need Dimensional space (i.e. on latitude direction) is extended, so that it is consistent with longitudinal span, such as the white grid regions in Fig. 4 Domain.The longitude and latitude regional code of extension exists, but geographically and is not present, and the part is continuous.

In one embodiment, the coordinate origin of multiple grids can be determined according to latitude and longitude coordinates.Latitude and longitude coordinates Choice of Origin under the line with the point of intersection of the first meridian, but due to multiple dimensioned integer coding support be Morton code friendship Fork-shaped formula does not support the processing of complement code directly, and what Morton code was supported in other words is signless integer, it is therefore desirable to net The Choice of Origin of lattice carries out planning processing, so that all grids are all in positive direction.The present embodiment selects (- 180 °, -180 °) Origin coordinates origin as earth subdivision grid.The characteristics of according to multiple dimensioned coding, coordinate translation are equivalent to integer translation, because This needs ± 0x80000000 on the basis of integer coding.

According to quaternary tree mode to the longitude and latitude space after extension is divided after, the grid of obtained each level meets Below:

The corresponding whole world of 0 grade of grid, range are 360 ° × 180 °；

Corresponding 1/4 earth of 1 grade of grid, range are 180 ° × 90 °；

Corresponding 1/8 earth of 2 grades of grids, range are 90 ° × 90 °；

Corresponding 1/32 earth of 3 grades of grids, range are 45 ° × 45 °；

Corresponding 1/128 earth of 4 grades of grids, range are 22.5 ° × 22.5 °；

…

31 grades of grids are correspondingA earth, range are

By weight of longitude and latitude, the earth subdivision grid of the 31st level under the line on scale be about: 1.86cm × 1.86cm (terrestrial equator radius is calculated according to 6378140m).The scale supports the smallest dimension letter of all kinds of remotely-sensed datas at present Breath, and support from the entire earth to the grid of totally 31 kinds of different scales 1.86cm.

To sum up, earth subdivision grid as shown in FIG. 6 can be obtained.

Secondly, carrying out integer coding to the corresponding coordinate information of each grid, corresponding first coding of each grid is obtained Value.

The detailed process of multiple dimensioned integer coding is detailed in above-described embodiment, details are not described herein again.Each grid pair It is as shown in Figure 5 that the coordinate information answered carries out the encoded radio obtained after integer coding.

Finally, establishing grid and corresponding first encoded radio of each grid (encoded radio i.e. shown in fig. 5) and remote sensing shadow As the incidence relation between data.

In the step, using a small amount of and closed closest to the grid of image capturing range and its encoded radio to remote sensing image data Connection, so that the subsequent query of remote sensing image data and statistical work become inquiry and statistics to space lattice.Since multi-source is distant The scape (different type, different brackets) for feeling image has different scales, position and coverage area, and earth subdivision grid is rigid Property multiple dimensioned (fixed multiple dimensioned) grid framework, therefore there is the not corresponding situation of can completely between the two, such as A~F in Fig. 7 The remote sensing image data of representative, respectively situation corresponding between the earth subdivision grid of each level.

The spatial relationship of multiple dimensioned earth subdivision grid and remote sensing image overlay area as shown in figure 8, as seen from Figure 8, Remote sensing image overlay area 70 includes that remote sensing image overlay area 70 includes grid a, remote sensing image by the grid 71 of lowest hierarchical level Overlay area 70 is intersected with grid b, remote sensing image overlay area 70 and grid c phase from, therefore, grid and the remote sensing image area of coverage Spatial relationship between domain may include following four: grid includes by image comprising image, grid, grid intersects with image, net Lattice and image mutually from, be briefly referred to as include, by comprising, intersection, mutually from.

According to above-mentioned analysis, the technical program is closed using the top-down space for judging grid and image overlay area step by step System, the gridding association of Lai Shixian image, and its coding result need to meet following three qualifications: level is encoded after gridding Lower limit N_min, level upper limit N is encoded after gridding_max, group/cording quantity upper limit S after gridding_max.By above three qualifications, User can adjust gridding association results according to image data concrete condition, to obtain optimal result.

Therefore, following steps realization can be used in the gridding association of image:

Step A1, the boundary rectangle of the corresponding imagery zone of remote sensing image data is calculated.

Step A2, gridding is carried out to boundary rectangle using the grid of lowest hierarchical level, obtains corresponding with boundary rectangle the One grid.

That is, being N using level_minGrid gridding, the concrete mode of gridding and above-mentioned S306 are carried out to boundary rectangle The gridding method of middle image data is identical, and details are not described herein.

Step A3, judge whether the first grid meets preset condition；Wherein, preset condition includes at least one of the following: The sum of first kind grid and the second class grid in one grid is greater than or equal to preset group/cording quantity upper limit value (i.e. gridding Group/cording quantity upper limit S afterwards_max), the number of grid of the second class grid be zero.

Wherein, first kind grid includes: the grid for including by boundary rectangle, or, trellis stage is more than or equal to preset Level upper limit value (encodes level upper limit N i.e. after gridding_max) grid.Second class grid include: intersect with boundary rectangle and Trellis stage is less than level upper limit value and (encodes level upper limit N i.e. after gridding_max) grid, or, the net comprising influence area Lattice.

If step A4, the first grid meets preset condition, it is determined that the first grid is grid corresponding with imagery zone.

If the first grid is unsatisfactory for preset condition, execution following steps are recycled, until the grid finally obtained meets State preset condition: for the second grid of highest priority in each second class grid, level is adjacent where determining the second grid The sub-grid of level；Sub-grid is divided into first kind grid and/or the second class grid, and updates the preferential of the second class grid Grade.

Above-mentioned steps are specific as follows: the adjacent level of level where first calculating the grid of highest priority in the second class grid 4 sub-grids, and by 4 sub- grid dividings be first kind grid and/or the second class grid, while update the second class grid In priority orders；Then judge whether the first kind grid obtained after 4 sub- grid dividings and/or the second class grid are full The above-mentioned preset condition of foot, if satisfied, grid corresponding with imagery zone is then exported, if not satisfied, then redefining out second The grid of highest priority in class grid, and repeat the above steps for the grid determined, until the grid obtained meets State preset condition.

Wherein, the priority orders of the second class grid can be determined according at least one in following rule:

(1) clathrum series is smaller, the higher of priority

(2) the specified type number of grid in the sub-grid that grid is included is more, and priority is higher；The specified type net Lattice include the grid that region corresponding with remote sensing image data mutually separates.

(3) first kind number of grid in the sub-grid that grid is included is fewer, and priority is higher.

Fig. 9 is the schematic flow chart of a kind of temporal index and querying method according to one embodiment of this specification, such as Fig. 9 It is shown, this method comprises:

S902 creates the data directory for managing the source data of remote sensing image data；Include in data directory with At least one of lower field: the source data of remote sensing image data, Data Identification, data obtaining time coding, data correlation grid are compiled Code collection is closed.

S904 creates the net for establishing the spatial index of remote sensing image data according to each field in data directory Trellis coding concordance list；It include earth subdivision grid coding and Data Identification collection in grid coding concordance list.

Wherein, earth subdivision grid coding is carried out corresponding to each grid obtained after subdivision to the longitude and latitude space of the earth Encoded radio.

S906 inquires remote sensing image data to be checked according to data directory and grid coding concordance list.

Using the technical solution of this specification one or more embodiment, by creating for managing remote sensing image data The data directory of source data, and according to each field in data directory, create the space for establishing remote sensing image data The grid coding concordance list of index, and then according to the data directory of creation and grid coding concordance list, to remote sensing to be checked Image data is inquired.The space-time code method of tissue as it can be seen that the technical solution can be formatted based on World Wide Web, establishes multi-source The spatial index of remote sensing image data, and believed using the time that multiple dimensioned period integer coding method unifies remote sensing image data Breath, with this settling time index, it is achieved that the high-efficiency tissue of massive multisource remote sensing data, management and inquiry.

Each step in above-described embodiment is described in detail below.

Firstly, creation for manage the source data of remote sensing image data data directory (Data Index Table, DIT)。

The present embodiment realizes the Organization And Management to remotely-sensed data by the way of file and database.Specifically, remote sensing Image data is stored in data-storage system in the form of a file, and the source data of remote sensing image is managed with database, and with This completes the time of remotely-sensed data and the creation and update of spatial index.User can select according to self-demand and appointed condition It selects data-storage system type (such as distributed file storage system, cloud storage system) and type of database is (such as pure relational Database Oracle, non-relational database HBase etc.).

Below by taking oracle database as an example, creation data directory DIT is illustrated how.

Creation data directory DIT is the concordance list to construct the attribute information of remote sensing image data, to multi-source remote sensing number According to attribute data once standardized, by the data in source data table set (Source DataTable Set, SDTS) Carry out unified management.Data directory DIT is in addition to further including with lower word comprising the data field in source data table set SDTS Section: Data Identification ImageID, the i.e. unique identification of remote sensing image data, data type is integer, as data directory The major key of DIT table, and be set as increasing certainly, in order to unifying identifier remote sensing image data improves remote sensing image data mark Search efficiency；Data obtaining time encodes TimeCode, and data type is 64 signless integers, for recording remote sensing shadow As the acquisition time of data, it can be used for settling time index, in order to unified time format and the inquiry effect for improving the time Rate；Data correlation grid coding set GridCodeSet, data type is 64 signless integer arrays, distant for recording Feel the corresponding association code set of image data.

Based on the field contents in above-mentioned data directory DIT, in one embodiment, number can be created in accordance with the following steps According to concordance list DIT:

Step B1, the field name and data type of each field in data directory DIT are determined.

Table 2 schematically lists the field name and data type of each field in data directory DIT in an embodiment, It should be noted that in practical applications, in data directory DIT in addition to the field that table 2 is enumerated, can also increase as needed Other relevant fields.

Step B2, reading source data corresponding with field name and data type, and the source data read is imported In data directory DIT.

Wherein, source data can be stored in source data table set (Source DataTable Set, SDTS), source data table Set SDTS is to indicate source data table all in data management system, it is generally the case that the image data attribute of different satellites Information is stored in different source data tables.Therefore, it can be read from source data table set SDTS and field name and data class The corresponding source data of type imports in data directory DIT.

In table 2, field " image unique identification " i.e. Data Identification；Field " acquisition time coding " i.e. data obtaining time is compiled Code；Field " association code set " i.e. data correlation grid coding set.Wherein, the acquisition time coding of data can pass through the time Code conversion method is converted to the temporal information of remote sensing image data, such as carries out the integer coding of specified level, when Between code conversion method will be described in detail later, in the present embodiment, time encoding level (i.e. specified level) is 63 grades, the time Precision is 1 microsecond.Be associated with code set can through the foregoing embodiment in be described in detail grid relationship method to remote sensing image data Location information be converted to, details are not described herein again.

Table 2

In one embodiment, parallel work-flow can be used to import multiple source datas in data directory DIT, to improve number According to importing efficiency.

In one embodiment, it is safeguarded to prevent the big data quantity in data directory DIT to be inserted into data with inquiry, table It is affected with update and database performance, following strategy can be used: when the disk occupied space of data directory DIT is more than pre- If when threshold value (such as 1GB), carrying out table division operation, several points of tables are formed, make every table is divided to be no more than preset threshold (such as 1GB).It include: firstly, data directory DIT is in logic to the advantages of data directory DIT carry out table subregion in the present embodiment On be still a table, realize the unified management of data；Secondly, can be mapped to several disks defeated to realize for several points of tables Enter/the equilibrium of output port I/O, improves the performance of system；Again, the reliability of data directory DIT is enhanced, for example, working as When some table subregion failure, other table subregions are still available, need to only repair the table subregion of failure.

Step B3, each field in global image density look-up table is determined.

Wherein, global image density look-up table includes at least one of following field: distant in subdivision grid coding, each grid Feel degree of covering, the number of grid upper limit value multiplying power of image data.

Global image density look-up table (Global Image Density Look-up Table, GIDLUT) is for storing The distribution density etc. of global image, purpose include: the distribution situation that (1) makes the global image of user's grasp；(2) according to global shadow As density look-up table GIDLUT, user can optimize restrictive condition S when carrying out gridding association to remote sensing image data_max(i.e. net The group/cording quantity upper limit after formatting).Earth subdivision trellis stage used by global image density look-up table GIDLUT is N_GIDLUT (if without particular/special requirement, which is 5, then earth surface is divided into 512 grids), which is not more than remote sensing image data grid Change restrictive condition N when association_min(level lower limit is encoded i.e. after gridding).

Field in global image density look-up table GIDLUT is as shown in table 3, in which: Code is subdivision grid coding, table Show N_GIDLUTGrade earth subdivision grid coding, the value initialize before data import data directory DIT, base area net lattice Subdivision coding method (carries out the coding of integer coding to the corresponding coordinate information of grid described in S304 i.e. in above-described embodiment Method) it is calculated；Count is image degree of covering, indicates that remote sensing image covers in the corresponding grid of subdivision grid coding Code The number of lid, the calculating of the value is synchronous when data import data directory DIT to be carried out, according in data directory DIT GridCodeset field calculates；Alpha is number of grid upper limit value multiplying power, indicates that subdivision grid coding Code is corresponded in grid Remote sensing image data gridding be associated with when restrictive condition S_max(i.e. the group/cording quantity upper limit after gridding) and image length-width ratio Ratio between example k, i.e. S_max=Alphak.The value calculates after data import data directory DIT, and calculation formula is as follows:

Wherein, Alpha (i) indicates the corresponding image number of grid upper limit value multiplying power of i-th of grid coding Code (i)； Count_minWith Count_maxIt respectively indicates in global image density look-up table GIDLUT in Count field (i.e. image degree of covering) Minimum value and maximum value；The value user of m can select as needed, and value range is [0,2¹⁰]；α₀Value user It can be selected as needed, and value range is [4,2¹⁰]。

Table 3

Serial number	Field name	Field meanings	Data type	Remarks
					1	Code	Subdivision grid coding	NUMERIC(20,0)	Major key
2	Count	Image degree of covering	NUMERIC(38,0)
					3	Alpha	Number of grid upper limit value multiplying power	NUMERIC(10,0)

Step B4, it is encoded according to data obtaining time, establishes the time encoding index of remote sensing image data.

The following detailed description of how being converted by temporal information of the time encoding conversion method to remote sensing image data, To obtain the acquisition time coding of remote sensing image data.In time encoding conversion method described below, time encoding level is 63 grades, the temporal information of the enough accurate expression remote sensing image datas of the level,

Firstly, the temporal information of remote sensing image data to be decomposed into the integer of multiple specified time scales.

Wherein, specified time scale may include year (A), the moon (B), day (C), hour (D), point (E), the second (F), millisecond (G), Microsecond (H) etc..

Secondly, the integer on each specified time scale is separately encoded the binary number to specify digit, remote sensing shadow is obtained As binary code value of the temporal information on each specified time scale of data.

Wherein, each specified time scale respectively corresponds respective specified digit, and finger corresponding to each specified time scale Positioning number can be identical, can also be different.Also, it, can be by a high position " 0 " polishing for each binary code obtained after conversion.

Specifically, the corresponding specified digit of specified time scale " year " is 17bit, the corresponding finger of specified time scale " moon " Positioning number is 4bit, and the corresponding specified digit of specified time scale " day " is 5bit, the corresponding finger of specified time scale " hour " Positioning number is 5bit, and the corresponding specified digit of specified time scale " dividing " is 6bit, and specified time scale " second " is corresponding specified Digit is 6bit, and the corresponding specified digit of specified time scale " millisecond " is 10bit, the corresponding finger of specified time scale " microsecond " Positioning number is 10bit.Based on this, the integer on each specified time scale can be separately encoded as follows to specify digit Binary number:

The binary number for being 17bit by the integer coding in this time scale of year；

The binary number for being 5bit by the integer coding in this time scale of day；

The binary number for being 5bit by the integer coding in this time scale of hour；

The binary number for being 6bit by the integer coding in point this time scale；

The binary number for being 6bit by the integer coding in this time scale of second；

The binary number for being 10bit by the integer coding in this time scale of millisecond；

The binary number for being 10bit by the integer coding in this time scale of microsecond.

Again, the temporal information of remote sensing image data is carried out on bit field in the binary code value on each specified time scale Connection, obtain corresponding second encoded radio of temporal information of remote sensing image data.

The second encoded radio that this step obtains is single scale time encoding corresponding to the temporal information of remote sensing image data Value.

Finally, the second encoded radio is moved to left one, the data obtaining time for obtaining remote sensing image data is encoded, i.e., the 63rd The multiple dimensioned time encoding value of level.

Figure 10 is schematically shown to time encoding conversion method performed by remote sensing image data.

As shown in Figure 10, the temporal information of remote sensing image data are as follows: year (A): the 2018, moon (B): 10, day (C): 1, hour (D): 13, dividing (E): 30, second (F): 29, millisecond (G): 300, microsecond (H): 0.

It is obtained after carrying out Binary Conversion to the temporal information: year (A): the 00000011111100010, moon (B): 1010, Day (C): 0001, hour (D): 01101, divide (E): 011110, second (F): 011101, millisecond (G): 0100101100, microsecond (H):0000000000。

After above-mentioned binary code value is spliced, when obtaining single scale corresponding to the temporal information of remote sensing image data Between encoded radio: 71024248425328640.

Obtained single scale time encoding value is moved to left one again, the multiple dimensioned time encoding of the 63rd level can be obtained Value: 142048496850657280.

The above-mentioned creation method that data directory DIT is described in detail, data directory DIT may be used also after the completion of creation It is updated.Specifically, when receiving the specified operational order to the first data in data directory DIT, to this First data execute corresponding operating, and update the index of the time encoding in data directory DIT.Wherein, operational order packet is specified Include deletion instruction and/or inserting instruction.

For example, when needing to delete the data in data directory DIT, the data deleted needed for deleting, and it is right Field " acquisition time coding " in data directory DIT rebuilds a secondary index, while updating global image density and searching Table GIDLUT, with the update of complete paired data concordance list DIT.

After having created data directory DIT, according to each field in data directory DIT, grid coding rope is created Draw table.Grid coding concordance list (Grid Code Index Table, GCIT) indicates the one-dimensional coding rope formed by grid coding Draw table, for storing the spatial index of remote sensing image data.

Establishing grid coding concordance list GICT is for two-dimensional spatial index is converted to one-dimensional coding index (with B-tree etc. Based on one-dimensional index), inherently improve building and the search efficiency of index.As shown in table 4, in the present embodiment, grid is compiled Include at least following field in code concordance list GICT: CodeIndex, i.e. earth subdivision grid coding are grid coding concordance lists The major key of GICT, the sequence of the value can establish one-dimensional coding index in conjunction with B-tree；ImageIDSet, i.e. data directory DIT The set of middle field ImageID composition, indicates that field CodeIndex corresponds to the Data Identification set of the covered image of grid, should Field stores in binary form, so that the splicing of ImageID and fractionation have high efficiency, to guarantee the query and search effect of data Rate.

Table 4

In one embodiment, grid coding concordance list GICT can be created according to step as shown in figure 11.

Step C1, according to data directory DIT and global image density look-up table GIDLUT, data directory DIT is determined In each data corresponding to data correlation grid coding set (i.e. the coding of image association gridding shown in Figure 11).

It step C2, will be in data correlation grid coding set interpenetration network code index table GICT.

Step C3, the data correlation grid coding set in interpenetration network code index table GICT is ranked up, to build The one-dimensional space index of vertical grid coding concordance list.

In the present embodiment, after being ranked up to data correlation grid coding set, one-dimensional space index can be established, such as B-tree indexed.

After tested, under identical environment (Oracle11g, Intel Xeon X5650@2.67GHz), it is assumed that data directory It is stored in 10,000,000 datas in table DIT, then it is time-consuming respectively to construct grid coding concordance list GICT, Oracle Spatial index 5.6 minutes, 36.9 minutes, it can be seen that building grid coding concordance list GICT has very high efficiency.

The above-mentioned creation method that grid coding concordance list GICT is described in detail, grid coding concordance list GICT are being created Cheng Hou can also be updated it.Depend on whether sending out in data directory DIT whether the update of grid coding concordance list GICT Raw data movement.Specifically, judging whether changed data volume reaches when the data in data directory DIT change To preset data amount；If so, re-creating grid coding concordance list GICT；If it is not, then without updating grid coding concordance list GICT.Wherein, the variation includes insertion data and/or deletion data.

For example, if insert data in data directory DIT, and it (has been more than preset data that the data volume be inserted into is larger Amount), then rebuild grid coding concordance list GICT, otherwise only need to according to building grid coding concordance list GICT method will be new Data interpenetration network code index table GICT in.If the data volume for deleting data in data directory DIT, and deleting Larger (being more than preset data amount), then rebuild grid coding concordance list GICT, otherwise do not operate.

The present embodiment only rebuilds grid coding concordance list GICT in the larger Shi Caihui of the data volume of variation, and reason is When being inserted into or delete fairly large data, efficiency is lower than the efficiency for rebuilding grid coding concordance list GICT.

It, can be right according to data directory DIT and grid coding concordance list GICT after creating grid coding concordance list GICT Remote sensing image data to be checked is inquired.

As can be seen from the above embodiments, remote sensing image and query region are primarily present four kinds of spatial relationships: the two was both non-intersecting Also mutually do not include, both intersection, query region include image, image include query region.Therefore, remote sensing image and interrogation zone There are following three kinds of relationships for the grid formed after the gridding of domain: the two grid mutually from, area grid be image grid father's unit, Area grid is the subelement of image grid.Based on this, remote sensing image data can be inquired according to following any mode.

Mode one, inquires remote sensing image data by the way of space querying, comprising the following steps:

Step D1, determine that the query region in earth subdivision grid, the query region are made of multiple latitude and longitude coordinates.

In the step, user can input inquiry region boundary vector coordinate data, query region can be any polygon Shape region.

Step D2, according to latitude and longitude coordinates and earth subdivision grid coding, corresponding first grid of latitude and longitude coordinates is determined Change coding.

Step D3, the first sub-grid that corresponding each grid is encoded for the first gridding, in grid coding concordance list The data correlation grid coding set for belonging to the first sub-grid is searched in GICT, as the first query result.

Wherein, the first sub-grid of grid is obtained sub-grid after the division to grid progress more high-level.Net Pair between the corresponding latitude and longitude coordinates of each grid and data correlation grid coding set is stored in trellis coding concordance list GICT It should be related to, therefore, when lookup belongs to the data correlation grid coding set of the first sub-grid, can first determine the first sub-grid pair The latitude and longitude coordinates answered, and then find out data correlation grid coding collection corresponding to the corresponding latitude and longitude coordinates of the first sub-grid It closes.

It should be noted that the corresponding latitude and longitude coordinates of the first sub-grid may be a coordinate range, at this point, being searched Data correlation grid coding set corresponding to the corresponding latitude and longitude coordinates of the first sub-grid arrived is corresponding with coordinate range Code set.

Step D4, it determines that the first gridding encodes father's grid of corresponding each grid, and determines father's grid correspondingly Ball subdivision grid coding is that the second gridding of query region encodes.

Specifically, the encoded radio of all father's grids for encoding corresponding grid in the first gridding coding is determined first, Then duplicate encoded radio is deleted, the second gridding coding of query region can be obtained.Wherein, level where father's grid isI.e. image gridding when minimum level.

Step D5, the second sub-grid that corresponding each grid is encoded for the second gridding, in grid coding concordance list The data correlation grid coding set for belonging to the second sub-grid is searched in GICT, and screens the second son for meeting the first preset condition The corresponding data correlation grid coding set of grid is as the second query result.

Wherein, the first preset condition includes: the gridding volume for existing in the first gridding coding and belonging to the second query result Code.

Step D6, it is extracted from grid coding concordance list GICT corresponding with the first query result and the second query result First Data Identification collection obtains the corresponding space querying result of remote sensing image data to be checked.

Specifically, can be extracted from grid coding concordance list GICT corresponding with the first query result and the second query result Binary code value (i.e. grid coding value), then by the binary code value extracted disassemble at ImageID (i.e. Data Identification) Set, and delete repetition values and the first Data Identification collection corresponding with the first query result and the second query result can be obtained.Its In, when binary code value dismantling is gathered at ImageID, multithreading progress can be used, to improve dismantling efficiency.

After tested, under identical environment (Oracle11g, Intel Xeon X5650@2.67GHz), it is assumed that data directory In table DIT be stored in 10,000,000 datas, then when inquiring each province's region overlay image, using above-mentioned spacing query method with Traditional Oracle Spatial is compared, and efficiency can averagely promote 10 times or more.

Mode two, inquires remote sensing image data by the way of spatial-temporal query, comprising the following steps:

Step E1, determine that the query argument of remote sensing image data to be checked, query argument include query region and inquiry Time range；Query time range includes initial time and termination time.

Step E2, time encoding is carried out to initial time and termination time respectively, obtains initial time encoded radio and termination Time encoding value.

In the step, when carrying out time encoding respectively to initial time and termination time, it can be described in detail according to above-described embodiment Time encoding conversion method carry out, details are not described herein again.

Step E3, it is directed to the first Data Identification collection, is filtered out from data directory DIT corresponding with the first Data Identification collection And meet the second preset condition field, as the corresponding spatial-temporal query result of remote sensing image data to be checked.

Wherein, when the second preset condition includes: that the corresponding data obtaining time coding of the first Data Identification collection is located at starting Between encoded radio and terminate time encoding value between.

First Data Identification collection described in the step is mentioned from grid coding concordance list GICT in as above-mentioned steps D6 The the first Data Identification collection corresponding with the first query result and the second query result taken out.As it can be seen that spatial-temporal query method need to be built It stands on the basis of spacing query method.

By the various embodiments described above it is found that the technical effect of the technical solution is including at least as follows: utilizing the more rulers of Global Grid It spends integer coding method, realizes that the multi-source of areal, multiple dimensioned, space scale and position are associated between multi-temporal data, It is advantageously implemented the integration of the data such as trans-regional, trans-departmental and shares；In addition, the spatial index based on grid coding constructs efficiency Height, Dynamic Maintenance are simple, and are suitble to parallel building, therefore the index that can be suitable for magnanimity, high dynamic remotely-sensed data be built It is vertical；Furthermore the image space search efficiency based on grid coding and index is high, is suitble to large-scale data query；Finally, the skill Art scheme is applicable not only to the Organization And Management of point scape image and product, is also applied for multi-resolution pyramid image.

To sum up, the specific embodiment of this theme is described.Other embodiments are in the appended claims In range.In some cases, the movement recorded in detail in the claims can execute and still in a different order Desired result may be implemented.In addition, process depicted in the drawing not necessarily requires the particular order shown or continuous suitable Sequence, to realize desired result.In some embodiments, multitasking and parallel processing can be advantageous.

The above are this specification one or more embodiment provide space-time code method, temporal index and querying method, Based on same thinking, this specification one or more embodiment also provide a kind of space-time code device, a kind of temporal index and Inquiry unit.

Figure 12 is the schematic block diagram according to a kind of space-time code device of one embodiment of this specification.As shown in figure 12, Space-time code device 1200 includes:

Subdivision module 1210 obtains multiple for carrying out subdivision according to longitude and latitude space of the specified partition patterns to the earth Grid；Wherein, longitudinal is identical with corresponding grid number on latitude direction；Specified partition patterns include specified subdivision layer Grade；

It is right respectively to obtain each grid for carrying out integer coding to the corresponding coordinate information of each grid for coding module 1220 The first encoded radio answered；

Relating module 1230, for establish grid and corresponding first encoded radio of each grid and remote sensing image data it Between incidence relation.

In one embodiment, subdivision module 1210 includes:

Expanding element, for being extended to latitude space, so that the sky of the spatial dimension in latitude space and longitude spaces Between range it is consistent, the longitude and latitude space after being expanded；

Division unit, for according to quaternary tree mode, the longitude and latitude space after extension to be divided into specified subdivision level Grid obtains multiple grids；And determine the coordinate origin of multiple grids.

In one embodiment, relating module 1230 includes:

Computing unit, for calculating the boundary rectangle of the corresponding imagery zone of remote sensing image data；

Gridding unit carries out gridding to boundary rectangle for the grid using lowest hierarchical level, obtains and boundary rectangle Corresponding first grid；

Judging unit, for judging whether the first grid meets preset condition；Wherein, preset condition includes following at least one : the sum of first kind grid and the second class grid in the first grid is greater than or equal to preset group/cording quantity upper limit value, the The number of grid of two class grids is zero；

First determination unit, if meeting preset condition for the first grid, it is determined that the first grid be and imagery zone phase Corresponding grid；

Execution unit recycles execution following steps if being unsatisfactory for preset condition for the first grid, until the net obtained Lattice meet preset condition: for the second grid of highest priority in each second class grid, level where determining the second grid The sub-grid of adjacent level；Sub-grid is divided into first kind grid and/or the second class grid, and updates the excellent of the second class grid First grade.

In one embodiment, first kind grid includes: the grid for including by boundary rectangle, or, trellis stage be greater than or Equal to the grid of preset level upper limit value；

Second class grid include: intersect with boundary rectangle and trellis stage be less than level upper limit value grid, or, comprising The grid of influence area.

In one embodiment, the priority orders of the second class grid are determined according at least one in following rule:

Clathrum series is smaller, and priority is higher；

Specified type number of grid in the sub-grid that grid is included is more, and priority is higher；Specified type grid packet Include the grid that region corresponding with remote sensing image data mutually separates；

First kind number of grid in the sub-grid that grid is included is fewer, and priority is higher.

Using the device of this specification one or more embodiment, subdivision is carried out by the longitude and latitude space to the earth, is obtained Integer coding is carried out to multiple grids, and to the corresponding coordinate information of each grid.The corresponding encoded radio of each grid is obtained, into And establish the incidence relation between grid and the corresponding encoded radio of each grid and remote sensing image data.As it can be seen that the technical side Case can be associated with to realize with the position of multi- source Remote Sensing Data data using the multiple dimensioned grid integer coding in the whole world, and then be realized with one-dimensional Integer coding carry out the spatial positional information of Unify legislation image.

Figure 13 is the schematic block diagram of a kind of temporal index and inquiry unit according to one embodiment of this specification.Such as Figure 13 Shown, temporal index and inquiry unit 1300 include:

First creation module 1310, for creating the data directory of the source data for managing remote sensing image data；Institute Stating in data directory includes at least one of following field: the source data, Data Identification, the data of the remote sensing image data Acquisition time coding, data correlation grid coding set；

Second creation module 1320, for creating described distant for establishing according to each field in the data directory Feel the grid coding concordance list of the spatial index of image data；It include earth subdivision grid coding in the grid coding concordance list And Data Identification collection；Wherein, the earth subdivision grid coding be obtained after carrying out subdivision to the longitude and latitude space of the earth it is each Encoded radio corresponding to grid；

Enquiry module 1330 is used for according to the data directory and the grid coding concordance list, to institute to be checked Remote sensing image data is stated to be inquired.

In one embodiment, the first creation module 1310 includes:

Second determination unit, for determining the field name and data type of each field in the data directory；

Reading unit, for reading the source data corresponding with the field name and the data type, and will The source data read imports in the data directory；

Third determination unit, for determining each field in global image density look-up table；The whole world image density is looked into Looking for table includes at least one of following field: degree of covering, the net of the remote sensing image data in subdivision grid coding, each grid Lattice the upper limit of the number value multiplying power；

Unit is established, for encoding according to the data obtaining time, establishes the time encoding of the remote sensing image data Index.

In one embodiment, the first creation module 1310 determines described in the data directory in the following way Data obtaining time coding:

In one embodiment, the first creation module 1310 includes:

Decomposition unit, for the temporal information to be decomposed into the integer of multiple specified time scales；

Coding unit, for the integer on each specified time scale to be separately encoded to the binary system to specify digit Number, obtains binary code value of the temporal information on each specified time scale；

Connection unit, for binary code value of the temporal information on each specified time scale to be carried out bit field On connection, obtain corresponding second encoded radio of the temporal information；

Shift cells left obtains the number of the remote sensing image data for second encoded radio to be moved to left one It is encoded according to acquisition time.

Correspondingly, the coding unit be also used to execute it is at least one of following:

In one embodiment, device 1300 further include:

Execution and update module, for working as the specified operational order received to the first data in the data directory When, corresponding operating is executed to first data；And the time encoding updated in the data directory indexes；

In one embodiment, the second creation module 1320 includes:

4th determination unit, described in determining according to the data directory and the global image density look-up table Data correlation grid coding set corresponding to each data in data directory；

It is inserted into unit, for the data correlation grid coding set to be inserted into the grid coding concordance list；

Sequencing unit, for being carried out to the data correlation grid coding set being inserted into the grid coding concordance list Sequence, to establish the one-dimensional space index of the grid coding concordance list.

In one embodiment, device 1300 further include:

Judgment module, for judging the data volume of the variation when the data in the data directory change Whether preset data amount is reached；Wherein, described change includes insertion data and/or deletion data；

Third creation module re-creates the grid coding if the data volume for variation reaches preset data amount Concordance list.

In one embodiment, the enquiry module 1330 includes:

5th determination unit, for determining the query region in the earth subdivision grid；The query region is by more A latitude and longitude coordinates are constituted；

6th determination unit, for determining the warp according to the latitude and longitude coordinates and the earth subdivision grid coding The corresponding first gridding coding of latitude coordinate；

First searching unit, for encoding the first sub-grid of corresponding each grid for first gridding, in institute The data correlation grid coding set searched in grid coding concordance list and belong to first sub-grid is stated, is looked into as first Ask result；

7th determination unit for determining that first gridding encodes father's grid of corresponding each grid, and determines institute The second gridding that the corresponding earth subdivision grid coding of father's grid is the query region is stated to encode；

Second searching unit, for encoding the second sub-grid of corresponding each grid for second gridding, in institute The data correlation grid coding set searched in grid coding concordance list and belong to second sub-grid is stated, and screens satisfaction The corresponding data correlation grid coding set of second sub-grid of first preset condition is as the second query result；Institute Stating the first preset condition includes: the gridding coding for existing in the first gridding coding and belonging to second query result；

Extraction unit, for being extracted from the grid coding concordance list and first query result and described second The corresponding first Data Identification collection of query result obtains the corresponding space querying knot of the remote sensing image data to be checked Fruit.

In one embodiment, the enquiry module 1330 includes:

8th determination unit, for determining the query argument of the remote sensing image data to be checked, the query argument Including the query region and query time range；The query time range includes initial time and termination time；

Time encoding unit is risen for carrying out time encoding respectively to the initial time and the termination time Beginning time encoding value and termination time encoding value；

Screening unit filters out and described for being directed to the first Data Identification collection from the data directory One Data Identification collection is corresponding and meets the field of the second preset condition, as the remote sensing image data to be checked Corresponding spatial-temporal query result；Wherein, second preset condition includes: the corresponding data acquisition of the first Data Identification collection Time encoding is located at the initial time encoded radio and terminates between time encoding value.

Using the device of this specification one or more embodiment, the space-time code side for the tissue that can be formatted based on World Wide Web Method establishes the spatial index of Methods on Multi-Sensors RS Image, and unifies remote sensing image using multiple dimensioned period integer coding method The temporal information of data, with this settling time index, it is achieved that the high-efficiency tissue of massive multisource remote sensing data, managing and looking into It askes.

It should be understood that above-mentioned space-time code device can be used to realize that previously described space-time is compiled Code method, and above-mentioned temporal index and inquiry unit can be used to realize previously described temporal index and querying method, wherein Datail description should be described with method part above it is similar, it is cumbersome to avoid, do not repeat separately herein.

System, device, module or the unit that above-described embodiment illustrates can specifically realize by computer chip or entity, Or it is realized by the product with certain function.It is a kind of typically to realize that equipment is computer.Specifically, computer for example may be used Think personal computer, laptop computer, cellular phone, camera phone, smart phone, personal digital assistant, media play It is any in device, navigation equipment, electronic mail equipment, game console, tablet computer, wearable device or these equipment The combination of equipment.

For convenience of description, it is divided into various units when description apparatus above with function to describe respectively.Certainly, implementing this The function of each unit can be realized in the same or multiple software and or hardware when specification one or more embodiment.

It should be understood by those skilled in the art that, this specification one or more embodiment can provide for method, system or Computer program product.Therefore, complete hardware embodiment can be used in this specification one or more embodiment, complete software is implemented The form of example or embodiment combining software and hardware aspects.Moreover, this specification one or more embodiment can be used one It is a or it is multiple wherein include computer usable program code computer-usable storage medium (including but not limited to disk storage Device, CD-ROM, optical memory etc.) on the form of computer program product implemented.

This specification one or more embodiment is referring to according to the method for the embodiment of the present application, equipment (system) and meter The flowchart and/or the block diagram of calculation machine program product describes.It should be understood that can be realized by computer program instructions flow chart and/ Or the combination of the process and/or box in each flow and/or block and flowchart and/or the block diagram in block diagram.It can These computer program instructions are provided at general purpose computer, special purpose computer, Embedded Processor or other programmable datas The processor of equipment is managed to generate a machine, so that holding by the processor of computer or other programmable data processing devices Capable instruction generates for realizing in one or more flows of the flowchart and/or one or more blocks of the block diagram The device of specified function.

These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates, Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or The function of being specified in multiple boxes.

These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one The step of function of being specified in a box or multiple boxes.

In a typical configuration, calculating equipment includes one or more processors (CPU), input/output interface, net Network interface and memory.

Memory may include the non-volatile memory in computer-readable medium, random access memory (RAM) and/or The forms such as Nonvolatile memory, such as read-only memory (ROM) or flash memory (flash RAM).Memory is computer-readable medium Example.

Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by any method Or technology come realize information store.Information can be computer readable instructions, data structure, the module of program or other data. The example of the storage medium of computer includes, but are not limited to phase change memory (PRAM), static random access memory (SRAM), moves State random access memory (DRAM), other kinds of random access memory (RAM), read-only memory (ROM), electric erasable Programmable read only memory (EEPROM), flash memory or other memory techniques, read-only disc read only memory (CD-ROM) (CD-ROM), Digital versatile disc (DVD) or other optical storage, magnetic cassettes, tape magnetic disk storage or other magnetic storage devices Or any other non-transmission medium, can be used for storage can be accessed by a computing device information.As defined in this article, it calculates Machine readable medium does not include temporary computer readable media (transitory media), such as the data-signal and carrier wave of modulation.

It should also be noted that, the terms "include", "comprise" or its any other variant are intended to nonexcludability It include so that the process, method, commodity or the equipment that include a series of elements not only include those elements, but also to wrap Include other elements that are not explicitly listed, or further include for this process, method, commodity or equipment intrinsic want Element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that including described want There is also other identical elements in the process, method of element, commodity or equipment.

This specification one or more embodiment can computer executable instructions it is general on It hereinafter describes, such as program module.Generally, program module includes executing particular task or realization particular abstract data type Routine, programs, objects, component, data structure etc..The application can also be practiced in a distributed computing environment, at these In distributed computing environment, by executing task by the connected remote processing devices of communication network.In distributed computing In environment, program module can be located in the local and remote computer storage media including storage equipment.

All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method Part explanation.

The foregoing is merely this specification one or more embodiments, are not limited to this specification.For this For the technical staff of field, this specification one or more embodiment can have various modifications and variations.It is all in this specification one Any modification, equivalent replacement, improvement and so within the spirit and principle of a or multiple embodiments, should be included in this explanation Within the scope of the claims of book one or more embodiment.

Claims

1. a kind of space-time code method, comprising:

Subdivision is carried out according to longitude and latitude space of the specified partition patterns to the earth, obtains multiple grids；Wherein, longitudinal and latitude It is identical to spend corresponding grid number on direction；The specified partition patterns include specified subdivision level；

Integer coding is carried out to the corresponding coordinate information of each grid, obtains corresponding first coding of each grid Value；

It establishes between the grid and corresponding first encoded radio of each grid and the remote sensing image data Incidence relation.

2. according to the method described in claim 1, described carry out subdivision to the longitude and latitude space of the earth according to specified partition patterns, Obtain multiple grids, comprising:

Latitude space is extended, so that the spatial dimension in the latitude space is consistent with the spatial dimension of longitude spaces, is obtained Longitude and latitude space after to extension；

According to quaternary tree mode, the longitude and latitude space after the extension is divided into the specified subdivision level grid, is obtained The multiple grid；And determine the coordinate origin of the multiple grid.

3. according to the method described in claim 1, described establish the grid and each grid corresponding described first Incidence relation between encoded radio and the remote sensing image data, comprising:

Gridding is carried out to the boundary rectangle using the grid of lowest hierarchical level, obtains corresponding with the boundary rectangle the One grid；

Judge whether first grid meets preset condition；Wherein, the preset condition includes at least one of the following: described The sum of first kind grid and the second class grid in one grid is greater than or equal to preset group/cording quantity upper limit value, described second The number of grid of class grid is zero；

If it is not, execution following steps are then recycled, until the grid obtained meets the preset condition: for each described second Second grid of highest priority in class grid, the sub-grid of the adjacent level of level where determining second grid；By institute It states sub-grid and is divided into the first kind grid and/or the second class grid, and update the priority of the second class grid.

4. according to the method described in claim 3, the first kind grid includes: the grid for including by the boundary rectangle, or, Trellis stage is greater than or equal to the grid of preset level upper limit value；

5. according to the method described in claim 3, the priority orders of the second class grid according to it is following rule at least One determination:

Clathrum series is smaller, and the priority is higher；

Specified type number of grid in the sub-grid that the grid is included is more, and the priority is higher；The specified class Type grid includes the grid that region corresponding with the remote sensing image data mutually separates；

6. a kind of temporal index and querying method, comprising:

Create the data directory for managing the source data of remote sensing image data；In the data directory include it is following at least One field: the source data of the remote sensing image data, Data Identification, data obtaining time coding, data correlation grid Code set；

According to each field in the data directory, the grid for establishing the spatial index of the remote sensing image data is created Code index table；It include earth subdivision grid coding and Data Identification collection in the grid coding concordance list；Wherein, the earth Subdivision grid coding is to carry out encoded radio corresponding to each grid obtained after subdivision to the longitude and latitude space of the earth；

According to the data directory and the grid coding concordance list, the remote sensing image data to be checked is looked into It askes.

7. according to the method described in claim 6, the data directory created for managing the source data of remote sensing image data Table, comprising:

Determine the field name and data type of each field in the data directory；

Read the source data corresponding with the field name and the data type, and the source data that will be read It imports in the data directory；

Determine each field in the image density look-up table of the whole world；The whole world image density look-up table includes at least one of following word Section: the degree of covering of the remote sensing image data in subdivision grid coding, each grid, number of grid upper limit value multiplying power；

8. according to the method described in claim 6, the data obtaining time in the data directory is encoded according to following Mode determines:

The integer coding that specified level is carried out to the temporal information of the remote sensing image data, obtains the remote sensing image data The data obtaining time coding.

9. according to the method described in claim 8, the temporal information to the remote sensing image data carries out specified level Integer coding obtains the data obtaining time coding of the remote sensing image data, comprising:

Integer on each specified time scale is separately encoded the binary number to specify digit, obtains the temporal information Binary code value on each specified time scale；

The temporal information is subjected to the connection on bit field in the binary code value on each specified time scale, is obtained described Corresponding second encoded radio of temporal information；

Second encoded radio is moved to left one, obtains the data obtaining time coding of the remote sensing image data.

10. according to the method described in claim 9, the specified time scale include year, month, day, hour, minute, second, millisecond, At least one of in microsecond；

Correspondingly, the integer by each specified time scale is separately encoded the binary number to specify digit, including At least one of below:

11. according to the method described in claim 7, further include:

When receiving the specified operational order to the first data in the data directory, phase is executed to first data It should operate；And the time encoding updated in the data directory indexes；

12. according to the method described in claim 6, each field according in the data directory, creates for establishing The grid coding concordance list of the spatial index of the remote sensing image data, comprising:

According to the data directory and the global image density look-up table, each data institute in the data directory is determined Corresponding data correlation grid coding set；

The data correlation grid coding set being inserted into the grid coding concordance list is ranked up, to establish the net The one-dimensional space of trellis coding concordance list indexes.

13. according to the method described in claim 6, further include:

When the data in the data directory change, judge whether the data volume of the variation reaches preset data Amount；Wherein, described change includes insertion data and/or deletion data；

If so, re-creating the grid coding concordance list.

14. right according to the method described in claim 6, described according to the data directory and the grid coding concordance list The remote sensing image data to be checked is inquired, comprising:

According to the latitude and longitude coordinates and the earth subdivision grid coding, corresponding first grid of the latitude and longitude coordinates is determined Change coding；

The first sub-grid that corresponding each grid is encoded for first gridding is searched in the grid coding concordance list The data correlation grid coding set for belonging to first sub-grid, as the first query result；

It determines that first gridding encodes father's grid of corresponding each grid, and determines the corresponding earth subdivision of father's grid Grid coding is that the second gridding of the query region encodes；

The second sub-grid that corresponding each grid is encoded for second gridding is searched in the grid coding concordance list Belong to the data correlation grid coding set of second sub-grid, and screens and meet described the second of the first preset condition The corresponding data correlation grid coding set of sub-grid is as the second query result；First preset condition includes: institute It states and there is the gridding coding for belonging to second query result in the first gridding coding；

Corresponding with first query result and second query result is extracted from the grid coding concordance list One Data Identification collection obtains the corresponding space querying result of the remote sensing image data to be checked.

15. it is according to the method for claim 14, described according to the data directory and the grid coding concordance list, it is right The remote sensing image data to be checked is inquired, comprising:

Determine that the query argument of the remote sensing image data to be checked, the query argument include the query region and inquiry Time range；The query time range includes initial time and termination time；

Time encoding is carried out to the initial time and the termination time respectively, obtain initial time encoded radio and terminates the time Encoded radio；

For the first Data Identification collection, filtered out from the data directory corresponding with the first Data Identification collection And meet the second preset condition field, as the corresponding spatial-temporal query knot of the remote sensing image data to be checked Fruit；Wherein, second preset condition includes: described in the corresponding data obtaining time coding of the first Data Identification collection is located at Between initial time encoded radio and termination time encoding value.

16. a kind of space-time code device, comprising:

Subdivision module obtains multiple grids for carrying out subdivision according to longitude and latitude space of the specified partition patterns to the earth；Its In, longitudinal is identical with corresponding grid number on latitude direction；The specified partition patterns include specified subdivision level；

It is right respectively to obtain each grid for carrying out integer coding to the corresponding coordinate information of each grid for coding module The first encoded radio answered；

Relating module, for establishing the grid and corresponding first encoded radio of each grid and the remote sensing shadow As the incidence relation between data.

17. a kind of temporal index and inquiry unit, comprising:

First creation module, for creating the data directory of the source data for managing remote sensing image data；The data rope Drawing in table includes at least one of following field: the source data of the remote sensing image data, Data Identification, data obtaining time Coding, data correlation grid coding set；

Second creation module, for creating for establishing the remote sensing image number according to each field in the data directory According to spatial index grid coding concordance list；It include earth subdivision grid coding and data mark in the grid coding concordance list Know collection；Wherein, the earth subdivision grid coding is that each grid obtained after subdivision institute is carried out to the longitude and latitude space of the earth is right The encoded radio answered；

Enquiry module is used for according to the data directory and the grid coding concordance list, to the remote sensing shadow to be checked As data are inquired.