# Patent

## Title
Method and system for multi-dimensional and geographic search

## Abstract

One embodiment of the present invention is directed to a multi-dimensional and geographic-search system that includes a search-region-receiving component that receives coordinates that define a region to be searched, an indexing component that generates an index for each of the received coordinates, a recursive-decomposition component that decomposes the region into sub-regions based on the indices generated by the indexing component, and a search component that selects points of interest from a points-of-interest database by comparing indices computed for the sub-regions to indices computed and stored in the database for each point of interest. Additional embodiments are directed to carrying out efficient, index-based searching for various additional types index-characterized entities.

## Background
CROSS-REFERENCE TO RELATED APPLICATION
This application claims the benefit of provisional patent application Ser. No. 61/140,553 filed Dec. 23, 2008 by the present inventor.


TECHNICAL FIELD
The present invention is directed to methods and systems for efficient search of points of interest in multi-dimensional and geographic space.
BACKGROUND OF THE INVENTION
With the advance of information technology in recent years, the problem of searching and retrieval of data records has been encountered in many different application areas. While the technology of string-based search has matured over the years, the search of multi-dimensional space remains an area of active research. In particular, web-based maps and mobile maps in cell phones have been popularized in recent years. These mapping technologies call for efficient systems to support geographic search of points of interest over geographic space. Multi-dimensional search also has applications in domains such as similarity search of database records, pattern recognition, near neighbor search, and matching problems.
SUMMARY OF THE INVENTION
One embodiment of the present invention is directed to a multi-dimensional and geographic-search system that includes a search-region-receiving component that receives coordinates that define a region to be searched, an indexing component that generates an index for each of the received coordinates, a recursive-decomposition component that decomposes the region into sub-regions based on the indices generated by the indexing component, and a search component that selects points of interest from a points-of-interest database by comparing indices computed for the sub-regions to indices computed and stored in the database for each point of interest. Additional embodiments are directed to carrying out efficient, index-based searching for various additional types index-characterized entities.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 illustrates a bit-permutation mapping F of points from a two-dimensional space to one-dimensional indices, according to one embodiment of the present invention.
FIG. 2 illustrates a function F−1, the inverse of a bit-permutation mapping for a three-dimensional space, according to one embodiment of the present invention.
FIG. 3 illustrates how the points in a 2-dimensional space are indexed with the mapping F, according to one embodiment of the present invention.
FIG. 4 illustrates how a 2-dimensional region with a pattern of (0, 1) is decomposed along a first dimension into two sub-regions, according to one embodiment of the present invention.
FIG. 5 illustrates how a 2-dimensional region with pattern (1, 0) is decomposed along a second dimension into two sub-regions, according to one embodiment of the present invention.
FIG. 6 illustrates how a 2-dimensional region with pattern (1, 1) is decomposed along both dimensions into four sub-regions, according to one embodiment of the present invention.
FIG. 7 is a diagram of a system that represents one embodiment of the present invention and that supports multi-dimensional/geographic search.



DETAILED DESCRIPTION OF THE INVENTION
The Bit Permutation Mapping
Certain embodiments of the present invention employ a bit-shuffle mapping from high dimensional space to one dimensional space and recursively decompose the region into smaller sub-regions with identifiable precisions. A search query to the given region is thus reduced to multiple searches over the linearly ordered keys derived from the mapping in the database, which is indexed.
It should be noted that method embodiments of the present invention are computational methods, encoded as computer instructions, that are necessarily carried out on an electronic computer system, or on other electronic-processor-containing devices. These methods involve far too many computations than can be feasibly undertaken by mental or hand calculation. For even simple cases, many thousands of basic arithmetic operations need to be carried out, and their results combined and processed. Hand calculation could take days or months, and would be extremely error-prone, while the problem domains in which embodiments of the present invention are applied involve real-time searching applications that are expected to produce results in seconds or fractions of seconds.
The space considered is a multi-dimensional space S, alternatively referred to as “d-dimensional,” where d is small integer greater than 1. Each point in the space has integer coordinates of the form (xd-1, . . . , x0), where 0≦xi≦(ni−1) for 0≦i≦(d−1). The shape of the space is given by the vector (nd-1, . . . , n0), the size of the space is given by the product




  
    ∏
    
      i
      =
      0
    
    
      d
      -
      1
    
  
  ⁢
  
    
  
  ⁢
  
    
      n
      i
    
    .
  




At first, a mapping from the d-dimensional space to one dimensional space is given as in Method 1 below. Note that Method 1 uses an operations Set(x, p, b), which sets the p-th least significant bit of integer x to the binary value b, and an operation Get(x, p), which returns the p-th least significant bit of integer x.











Method 1: F(xd−1,...,x0). Bit permutation mapping of a point


x=(xd−1,...,x0) with integer coordinates in d-dimensional


space to a point in a one-dimensional space.

















let



 B(xi) = (bn−1i,...,b0i) the binary form of xi, for i=0 to (d−1);



 y an integer of d × n bits



in



  for i =0 to (n−1)



   for j=0 to (d−1)



    Set (y, i × d + j , Get(xj,i))



   end



  end



  return y










Given a point x=(xd-1, . . . , x0) in a d-dimensional space, the value y=F(x) will be referred to as the index of the point.

An illustration of the method for dimensionality of two is given in FIG. 1. In FIG. 1, a point with coordinates (j, i) is mapped to an index k by permuting the binary bits of j and i to obtain the bits of k. In particular, in FIG. 1, the coordinate j is represented by an ordered set of 4 bits (j3, j2, j1, j0) 102, where j0 104 is the least significant bit. For example, k2 106, the second least significant bit of k, is obtained from the second least significant bit i1 108 of the ordered set of 4 bits (i3, i2, i1, i0) 110 that represent the i coordinate.
The bit permutation mapping F is one-one and onto, and therefore reversible. The reverse mapping of F is denoted by F−1, and can be thought as distributing the bits from the index of a point to the binary forms of the coordinates of the point in d-dimensional space from least-significant-to-most-significant order. The method for F−1 is given below:










Method 2. F−1 (y) for d-dimensional space

















let



 B(y) = (bn−1,...,b0)



 w =n/d



 xi integer of w bits, for i = 0 to (d−1)



in



 for p=0 to w



  for i = 0 to (d−1)



   Set (xi,p,Get(y,p × w + i))



  end



 end



 return (xd−1,...,x0).










The inverse mapping F−1 is illustrated in FIG. 2 for space dimensionality of three, using similar illustration conventions as used in FIG. 1. FIG. 3 illustrates the one-one correspondences, defined by the mappings of F and F−1, between points in a two-dimensional space with indices i and j that range in value from 0 to 7, and that can each therefore be represented by 3-bits, and their corresponding indices in an index k that ranges from 0 to 63, and that therefore can be represented by 6 bits. In FIG. 3, the values of index k are shown in a two-dimensional table indexed by the values of i and j.
Regions
A linear region defined by two indices k, and k2, where k1≦k2, denoted by R=[k1, k2], is the set of indices [k1, k2]={k|k1≦k≦k2}. The size of a linear region R, denoted S(R), is given by k2−k1+1. On the other hand, given two indices k1 and k2, where k1≦k2, F−1(k1)=)xd-11, . . . , x01), and F−1(k2)=(xd-12, . . . , x02), then a regular region of a d-dimensional space, where d>1, defined by the two indices k1 and k2, denoted “R=k1, k2,” is the set of indices given by:

k1,k2={F(xk-1, . . . ,x0)|min(xk-11,xk-12)≦xk-1≦max(xk-11,xk-12), . . . ,min(x01,x02)≦x0≦max(x01,x02)}

For instance, a regular region in an indexed two-dimensional space has a rectangular form. The size of a regular region R, also denoted “S(R),” is given by the product of the sizes along each of the dimensions, i.e. S(R)=(xd-12−xd-1+1)× . . . ×(x02−x01+1). For example, when d=2, [3, 15]={3, 4, 5, . . . , 15}, and (3, 15)={3, 6, 7, 9, 11, 12, 13, 14, 15}. S([3, 15])=13, S(<3, 15>)=9. It can be seen that a regular region is always a subset of the corresponding linear region when defined by the same two indices. In other words, given two indices k1 and k2, where k1≦k2, the following always holds:

k1,k2[k1,k2]

This leads to an important notion called perfect regions. A regular region in a d-dimensional space is perfect if and only the following holds:

k1,k2=[k1,k2]

In other words, a region is perfect if and only if its regular region and linear region are identical. A regular region is said to be a cube if its size along each of the dimensions are identical. A cube is said to be perfect if the region defined by the cube is perfect. A perfect cubical region in a two-dimensional space is also referred to as a perfect square. As an example, [0, 15]=0, 15={0, 1, . . . , 15}, in two-dimensional space, is a perfect square.

The precision of a regular region R=k1, k2 with respect to the corresponding linear region, denoted “P(k1, k2),” is given by

P(k1,k2)=S(k1,k2)/S([k1,k2])

In other words, the precision of a regular region with respect to its corresponding linear region is the ratio between the cardinalities of the regular region and the linear region. Clearly, the precision of a region is always a number between zero and one. It follows, from the definition, that, given two indices k1 and k2, where k1≦k2, when the regular region k1, k2 is perfect, then P(k1, k2)=1. Conversely, if P(k1, k2)=1, then k1, k2 is necessarily a perfect region.

Next, the notions of level and pattern of a region are introduced. Let k1, k2 be a d-dimensional region, where B(k1)=(bm-11, . . . , b01) and B(k2)=(Bm-12, . . . , b02) are the binary forms of k1 and k2 respectively. Let  be the bit-wise XOR operation, K=B(k1)(k2)=(bm-1, . . . , b0), bi the non-zero most significant bit of K, and Floor (x) be the operation that takes a real number x and returns the greatest integer that is smaller than or equal to x. Then Floor(i/d)+1 is referred to as the level of the region k1, k2, written Lvl(k1, k2, d). When Lvl(k1, k2, d)=L, the consecutive d bits (kLd-1, . . . , k(L-1)d) of K is referred to as the pattern of the region, written Ptn(k1, k2, d).
The following two methods compute the level of any region and pattern of a region with size greater than 1 in d-dimensional space respectively. Note that the pattern of region is used for the purpose of region decomposition only. A region of size 1 (hence with level 0) will never be decomposed, and the pattern of such a region is undefined.










Method 3: Lvl(k1,k2,d)

















if S([k1,k2]) = 1



  return 0



else



  let



   K = B(k1)  B(k2) = (bm−1,...,b0)



   bi the most significant non-zero bit in K



 in



   return Floor (i/d) + 1



end




















Method 4: Ptn(k1,k2,d)

















let



 L = Lvl(k1,k2,d)



 B(k1)  B(k2) = (bm−1,...,b0)



in



 return  (bLd−1,...,b(L−1)d)










For example, Lvl(3, 15, 2)=2, Lvl(33, 57, 2)=3, Ptn(33, 57, 2)=(0, 1), Ptn(12, 39)=(1, 0).

The following two properties of a region k1, k2 always hold once the level for the region, Lvl(k1, k2, d), is known:

    
    
        (1) S(k1, k2)=|k1, k2)| is smaller or equal to 2Ld.
        (2) The region is either identical or a sub-region of a perfect cubic region of the same level. In the case of two-dimensional space, the region is either identical or a sub-region of a perfect square of the same level.

Region Decomposition

    
    


For a given region R=k1, k2, a list of regions (R1, . . . , Rm), where Ri=k1i, k2i, is said to be an m-ary decomposition of R, denoted “(R1, . . . , Rm) ε D(R),” if and only if

    
    
        (1) Sub-Region: Ri ⊂ R for i=1 to m
        (2) Mutually Exclusive: Ri ∩ Rj=φ, if and only if i≠j, for 1≧j≧m;
        (3) Complete: R1∪ R2 . . . ∪ Rm=R.

Furthermore, a decomposition of a region R is said to be monotonic when the indices within each sub-region in the list of sub-regions that comprises the decomposition are smaller than any of the indices within sub-regions to the right of the sub-region. For example, in two dimensional space, (3, 3, 6, 7,9, 1112, 15) ε D3, 15 is 4-ary, and monotonic; (48, 51, 56, 59) ε D48, 59 is a binary and monotonic decomposition.

    
    


Note that sub-regions in a decomposition generally do not have the same precision as that of the region. For example, the region R=48, 59, in the above example, has a precision of P(48, 59)= 8/12=0.667. However the two sub-regions in the decomposition of R, as given in the above example, 48, 51 and 56, 59, both have a precision of one. Generally speaking, the sub-regions resulting from decomposition of a region have higher precision than the decomposed region.
Methods for Region Decomposition
A region can generally be decomposed in many different ways. The methods for region decomposition presented here maximize the precision of the decomposed sub-regions. To achieve this, a region is decomposed according to the pattern of the region. The number of possible patterns for a d-dimensional space is given by 2d−1. Therefore, 2d−1 sub-programs are required in the currently described region decomposition method. In the following, the method for d=2 is given in detail, which can be generalized to work with higher dimensional spaces.
The following are two auxiliary methods to be used in the decomposition:










Method 5. LB(k, L) Lower bound of an index k for a given level L.

















if L==0



 return k



else



 let B(k) = (bn,...,b2L,...,b0) in



 return B−1 (bn,...,b2L,0,...,0)



end










where B−1 is a function that takes a binary number as an argument and returns the corresponding integer.











Method 6. UB(k, L) Upper bound of an index for a given level L.

















if L==0



 Return k



else



 let B(k) = (bn,...,b2L,...,b0) in



 return B−1 (bn,...,b2L,1,...,1)



end










In the following, the decomposition of a region k1, k2 into m sub-regions (k1m-1, ksm-1, . . . , k10, k20) by a specific decomposition function D, will be denoted by D(k1, k2)=((k1m-1, k2m-1), . . . , (k10, k20)).











Method 5. D1(k1,k2) for pattern = (0,1) in 2-dimensional space.

















Let



 (x11,x01) = F−1(k1)



 (x12,x02) = F−1(k2)



 L = Lvl (k1,k2,2)



 k3 = UB(k1,L − 1)



 (x13,x03) = F−1(k3)



 k4 = F(x12,x03)



 k5 = F(x11,x03 + 1)



in



 return  ((k1,k4),(k5,k2))










For example, the 2-dimensional region 9, 15 has a pattern of (0, 1) and is decomposed by D1 along the dimension of i into sub-regions 9, 11 and 12, 15. FIG. 4, which uses the same illustration conventions as FIG. 3, shows sub-regions 9, 11 and 12, 15 within curved boundaries 402 and 404.











Method 6. D2(k1,k2) for pattern = (1,0) in 2-dimensional space.

















let



 (x11,x01) = F−1(k1)



 (x12,x02) = F−1(k2)



 L = Lvl (k1,k2,2)



 k3 = LB(k2,L − 1)



 (x13,x03) = F−1(k3)



 k4 = F(x13 − 1,x02)



 k5 = F(x13,x01)



in



 return  ((k1,k4),(k5,k2))










For example, the region 4, 15 has a pattern of (1, 0), and will be decomposed by D2 along the dimension of j into sub-regions of 4, 7 and 12, 15 FIG. 5, which uses the same illustration conventions as FIGS. 3 and 4, shows sub-regions 4, 7 and 12, 15 within curved boundaries 502 and 504.











Method 7. D3(k1,k2) for pattern = (1,1) in 2-dimensional space.

















let



 (x11,x01) = F−1(k1)



 (x12,x02) = F−1(k2)



 L = Lvl (k1,k2,2)



 m3 = UB(k1,L − 1)



 m4 = F(x11,x03 + 1)



 m5 = F(x13,x02)



 m6 = F(x13 + 1,x01)



 m7 = F(x12,x03)



 m8 = LB(k2,L − 1)



in



 return   ((k1,m3),(m4,m5),(m6,m7),(m8,k2))










For example, the region 3, 51 has a pattern of (1, 1), and is decomposed by D3 along both the dimensions of i and j into sub-regions 3, 15, 18, 27, 33, 39 and 48, 51. FIG. 6, which uses the same illustration conventions as FIGS. 3-5, shows sub-regions 3, 15, 8, 27, 33, 39 and 48, 51 within curved boundaries 602-605. Note that the three methods D1, D2, and D3 for region decomposition always produce monotonic decompositions.

With the divisions for all the three possible patterns of a region in place, in a 2-dimensional space, the following is the region decomposition method.










Method 8. D(k1,k2) decomposition of regions in 2-dimensional


space

















if precision (k1,k2) =1



 return ( (k1,k2) )



else



 let pat = Ptn(k1,k2,2) in



 case of pat



  (0,1) : D1(k1,k2);



  (1,0) : D2(k1,k2);



  (1,1) : D3(k1,k2);



 end



end










Furthermore, a region can be recursively divided, as shown in the following method:











Method 9. Div (k1,k2, Acc, Ps), where Acc is the list of sub-regions


initially set to the empty list [ ], Ps is a predicate over k1,k2 as


the termination condition of the search, which can be used to control the


precision of the search.

















if Ps (k1,k2)



 Acc += ((k1,k2))



 else



 Decom = D(k1,k2)



 for each (s,t) in Decom



  Acc += Div(s,t,Acc,Ps)



 end



  end



. return Acc










The predicate Ps can be, for example, a function that calculates the precision of the region and determines whether it is above a given threshold, a function that examines the level of the region, and stops further decomposition when the level is below a given threshold, or some combination of both. Of course, one can also take the approach of setting the threshold precision to the value of 1 as the termination condition. In that case, the region will be decomposed into a set a perfect regions that collectively cover exactly the searched region.

Let Ps-x denote a predicate that returns TRUE if and only if the precision of a given region is above or equal to the threshold value of x. Then, the following are two examples of the recursive division.

Div(3,27,[ ],Ps−0.9)=[(3,3],(6,7),(9,9),(11,11),(12,15),(18,19),(24,27)]

Div(12,51,[ ],Ps−1]=[(12,15),(24,27),(36,39),(48,51))

Regional Search

Given points of interest in a d-dimensional space without any predictable pattern of the distribution, the problem of regional search is to find all the points of interest that fall within a given region. A database table is used to hold the records related to the points of interest (POIs). For the purpose of presentation, the table will be referred to as “the TPOI table.” The system of this embodiment requires the augmentation of the table to include a column CKEY to hold the indices of the POIs, which are calculated with the bit-permutation mapping given in Method 1. Since the values in the column CKEY are integers, they can be indexed by the database system, and retrieved efficiently by the system.
A query of region search is made in the form of two corner points x1=(xd-11, . . . , x01) and x2=(xd-12, . . . , x02) of the d-dimensional space. The search is to retrieve each and every POI with coordinates (xd-1, . . . , x0) that satisfies the following condition:

min(xk-11,xk-12)≦xk-1≦max(xk-11,xk-12), . . . ,min(x01,x02)≦x0≦max(x01,x02)

The following method is used for the region search:











Method 10. Region-Search (x1,x2)

















k1 = F(x1)



k2 = F(x2)



if k1 > k2



   temp = k1



   k1 = k2



   k2 =temp



end



if precision (k1, k2) = 1



  sub-regions =[ (k1, k2) ]



else



 sub-regions = Div(k1,k2)



end



 POIS= [ ] #empty list



for each ( s, t ) in sub-regions



    POIS += Select * from TPOI where s ≦ CKEY ≦ t



end



return POIS










Geographic Search

The search of d-dimensional space has application in many different areas. Geographic search is one of them. In the following, it is assumed that the POIs are associated with latitudinal and longitudinal coordinates with a precision up to the m-th place after the decimal. A linear transformation is employed to transform a geographic point into a point in a two dimensional integer space. The following method gives the linear transformation.
Method 11. LT(lat, lng). Linear transformation of a geographic point given its latitude and longitude at the precision m.

i=(lat+90)*10m 

j=(lng+180)*10m 

return(i,j)

With this linear transformation, all points of interests can be assigned non-negative-integer coordinates in a two dimensional space. The search index of each POI can then be generated with the bit-permutation mapping F. For example:

lat=44.677198348794,lng=−122.120080823001,index=37459463583151357

lat=44.677911141839,lng=−122.405725461498,index=37458704662251054

lat=44.759668361174,lng=−122.058545902944,index=37459880683108980

The problem of geographic search is then reduced to the problem of a regional search in a two dimensional space. From the integer coordinates of points of interests, a column of indices is created by the calculation of the bit-permutation mapping. The column is indexed in the database for fast data retrieval. The following method performs geographic search given two corner points p1=(lat1, lng1) and p2(lat2, lng2), and returns all POI's in the TPOI table of the database that falls under the rectangular region with p1, p2 as the corner points.











Method 12. Geographic Search (p1,p2)

















let



 lat1 = Latitude(p1)



 lng1 = Longitude (p1)



 lat2 = Latitude(p2)



 lng2 = Longitude (p2)



 (i1,j1) = LT(lat1,lng1)



 (i2,j2) = LT(lat2,lng2)



 s1 = (i1,j1)



 s2 = (i2,j2)



in



 Region-Search (s1,s2)










As an example, let p1=(47.123456, −122.120088), p2=(47.809030, −122.450068), then the rectangular region will be decomposed into the sub-regions of

    
    
        ([37468767054828800, 37468882373443583],
        [37469020385476608, 37469138976698106],
        [37470278310237440, 37470282532782079],
        [37470286900172032, 37470291122716671],
        [37470312669975808, 37470316892520447],
        [37470321259910400, 37470394201931775],
        [37470415749190912, 37470419971735551],
        [37470424339125504, 37470428561670143],
        [37470450108929280, 37470454163035517],
        [37470462921408512, 37470522882512253],
        [37470531640885248, 37470779485766776])

which is then converted to database query based on the CKEY column values to retrieve all the POIs that falls within the queried region in the database.

Performance, Pragmatics, and Variants

    
    


FIG. 7 provides a diagram for a multi-dimensional/geographic-search system that utilizes components including mappings, linear transformations, and region decomposition and that represents one embodiment of the present invention. In FIG. 7, in step 702, the system determines whether a geographical search is to be undertaken. If so, then, in step 704, a linear transformation is undertaken to transform the real-valued longitudinal and latitudinal coordinates of two corner points that define a region to two-dimensional integer-valued coordinates. In step 706, a bit-permutation function F is applied to the integer-valued coordinates to generate one-dimensional indices. In steps 708-709, the pair of indices that define a region are placed in ascending order. In step 712, a recursive decomposition of the region, by above-described methods, is carried out to produce a list of sub-regions. In step 714, database commands are generated, based on the list of sub-regions, to retrieve POIs from the database.
The database may be a relational database managed by a relational-database-management system or other type of database that stores POI information on one or more mass-storage devices. Database management systems are generally software-implemented systems that execute on one or more computers. The multi-dimensional/geographic-search system may be a stand-alone software-implemented system, or may be incorporated, as a component, into other systems.
Given a region in a space, suppose U is the set of POIs that fall within the region, and V is the set of POIs actually returned by a searching system, then the precision and recall of the search are defined as follows:



  
    Precision
    =
    
      
        
        
          U
          ⋂
          V
        
        
      
      
        
        V
        
      
    
  




  
    Recall
    =
    
      
        
        
          U
          ⋂
          V
        
        
      
      
        
        U
        
      
    
  



The system that represents one embodiment of the present invention reduces a search of a regular region to the search of a set of sub-regions, which always collectively covers the regular region. The process can be thought as approximating a regular region with a set of linear regions. The precision of the approximation, and hence the precision of the search, is improved with each decomposition, while the recall of the search is always kept at one. The precision of the search by the system can be controlled and adjusted, by the predicate Ps of Method 9, thus allowing a trade-off between the precision and search time. When so desired, the precision can be set to the value of one.
According to embodiments of the present invention, computation of one search index of one POI is independent from computation of that of other POIs. Therefore, the contents of the CKEY column can be generated in large batches, as, for example, during initial incorporation of a large number of POIs obtained from a third-party data source, or in smaller batches or individually, as, for example, during incremental system data updates. The time taken per POI is constant, regardless of batch size. Similarly, deletions of POIs from the table are constant-time-per-POI operations.
The search method of various embodiments of the present invention does not rely on a tree-like data structure, relied on by many other searching methods. As a result, embodiments of the present invention are characterized by efficient use of storage space, simplicity, and search efficiency. Embodiments of the present invention do not traverse a tree-like data structure from top to bottom, one level at a time, as do many other searching methods. In many situations, regions can be decomposed in a number of steps logarithmically proportional to the size of the region, as opposed to the size of the search space. As an example, the region 15, 48 in the two-dimensional space can be completely decomposed in one step regardless the number of levels that actually exists in the space.
In the above presentation of Region-Search, the query to the database is not made until the list of sub-regions is completely generated. A variant scheme is to query the database as soon as a sub-region meets the standard imposed by the predicate Ps. The Region-Search can then continue while the database search is done in parallel, and possibly on a different server, and therefore further minimize the system response time. Parallelism in the search system of this embodiment can be further exploited by the observation that the queries to the database based on the sub-regions never overlap and can be done independently of each other.
There are times when it is desired to search for the POI's within a certain radius of a given center point, which is a search over a circularly-shaped region. The problem can be reduced to search of a square region inscribed by the circle. The calculation of the coordinates of the two corners of the square is a simple geometrical exercise. By the ratio of the areas of the circular-shaped region and square region, respectively, the precision of the search will be around 0.79 provided the distribution of the POI is uniform. When a precision of 1 is desired, the result can be easily filtered by a calculation of distance from the given center for each POI in the result produced by the region search based on the square. Similarly, embodiments of the present invention can be used to search regions with shapes other than rectangular or square regions by approximating the area with a set of rectangles. The irregular regions in space with dimensionality higher than two can be likewise approximated by a series of regular cubes before the region search based on regular shapes is applied.
Web-based or Mobile-based electronic map engine generally cut the map of a large space into tiles with different levels. When the map is shifted along different directions or zoomed in and out between different levels, the map needs to be updated with new tiles. In other words, when the map is moved or zoomed, an engine that powers the maps need to search for the tiles that fall within the view window, which is generally rectangular. Although map tiles are conceptually different from points of interest because, unlike a point, a tile takes up an area of space, embodiments of this invention can be used in an engine that powers the map. This is because a tile, usually shaped as a rectangle, can be represented by its four corner points, and a tile is then retrieved to update the changed map if and only if at least one of its corner points are located under the viewing region of the changed map. The embodiments of the present invention therefore can be used to power both the engine that drives an electronic map and the search engine for points of interests.

## Claims

1. A multi-dimensional and geographic-search system comprising:
a search-region-receiving component that receives a number i of sets of coordinates that define a region to be searched, the number i greater than 1 and each of the i sets of coordinates including d coordinates associated with a multidimensional space, each of the d coordinates represented by an ordered sequence of n bits;
an indexing component that generates, by a bit-permutation mapping, i integer indexes that include an integer index for each of the i sets of the received coordinates, the bit-permutation mapping grouping, for each of the n bits used to represent each of the d coordinates of a set of coordinates, the bits of the d coordinates at each bit position within the ordered-bit-sequence representations of the d coordinates together into a subsequence of bits within an ordered sequence of subsequences of bits that together compose the integer index;
a recursive-decomposition component that decomposes the region into sub-regions based on the integer indices generated by the indexing component; and
a search component that selects points of interest from a points-of-interest database by comparing integer indices computed for the regular sub-regions to indices computed and stored in the database for each point of interest.

2. The multi-dimensional and geographic-search system of claim 1 wherein the search-region-receiving component receives real-valued longitudinal and latitudinal coordinates for two points that define a geographical region.

3. The multi-dimensional and geographic-search system of claim 2 wherein the real-valued longitudinal and latitudinal coordinates for two points that define a geographical region are coordinates for two, diagonally opposed corners of a rectangular region.

4. The multi-dimensional and geographic-search system of claim 1 wherein the search-region-receiving component receives real-valued coordinates for two points that define a volume in a space of three or more dimensions.

5. The multi-dimensional and geographic-search system of claim 1 wherein the recursive-decomposition component decomposes the region into regular sub-regions based on the integer indices generated by the indexing component by:
when a termination condition is not satisfied,
decomposing the region according to a pattern computed for the region to generate two or more lower-level regular sub-regions; and
recursively decomposing the two or more regular sub-regions.

6. The multi-dimensional and geographic-search system of claim 5 wherein a region or regular sub-region is decomposed into lower-level regular sub-regions by:
computing a pattern for the region or regular sub-region;
using the computed pattern to select a decomposition method; and
applying the selected decomposition method to the region or regular sub-region.

7. The multi-dimensional and geographic-search system of claim 6 wherein a decomposition method, when applied to a region or regular sub-region:
identifies corner points for two or more regular sub-regions of the region or regular sub-region; and
computes one-dimensional integer indices, based on the identified corner points, for the two or more regular sub-regions of the region or regular sub-region.

8. The multi-dimensional and geographic-search system of claim 7 wherein the decomposition method decomposes a region or regular sub-region in such a way that the precisions of the regular sub-regions produced by the decomposition method are maximized.

9. A mapping system comprising:
a tile-specifying component that receives a number i of sets of coordinates that define region, the number i greater than 1 and each of the i sets of coordinates including d coordinates associated with a multidimensional space, each of the d coordinates represented by an ordered sequence of n bits;
an indexing component that generates, by a bit-permutation mapping, i integer indexes that include an integer index for each i sets of the coordinate that defines the region, the bit-permutation mapping grouping, for each of the n bits used to represent each of the d coordinates of a set of coordinates, the bits of the d coordinates at each bit position within the ordered-bit-sequence representations of the d coordinates together into a subsequence of bits within an ordered sequence of subsequences of bits that together compose the integer index;
a recursive-decomposition component that decomposes the region into regular sub-regions based on the integer indices generated by the indexing component; and
a tile-retrieving component that retrieves tiles corresponding to the regular sub-regions computed by the recursive-decomposition component.

10. The mapping system of claim 9 wherein the tile-specifying component receives real-valued longitudinal and latitudinal coordinates for two points that define the region.

11. The mapping system of claim 10 wherein the real-valued longitudinal and latitudinal coordinates for two points that define the region are coordinates for two, diagonally opposed corners of a rectangular region.

12. The mapping system of claim 9 wherein the recursive-decomposition component decomposes the region into regular sub-regions based on the integer indices generated by the indexing component by:
when a termination condition is not satisfied,
decomposing the region according to a pattern computed for the region to generate two or more lower-level regular sub-regions; and
recursively decomposing the two or more regular sub-regions.

13. The mapping system of claim 12 wherein a region or regular sub-region is decomposed into lower-level regular sub-regions by:
computing a pattern for the region or regular sub-region;
using the computed pattern to select a decomposition method; and
applying the selected decomposition method to the region or regular sub-region.

14. The mapping system of claim 13 wherein a decomposition method, when applied to a region or regular sub-region:
identifies corner points for two or more regular sub-regions of the region or regular sub-region; and
computes one-dimensional integer indices, based on the identified corner points, for the two or more regular sub-regions of the region or regular sub-region.

15. The mapping system of claim 14 wherein the decomposition method decomposes a region or regular sub-region in a such a way that the precisions of the regular sub-regions produced by the decomposition method are maximized.