时空查询是指检索给定时间和空间范围内的轨迹，创建时空表，索引公式Shard :: TimeBin :: XZ :: ID :: StartTime ，把第一条轨迹数据（T，X，Y）为：（0，0，0），（8，5，3），（17，8，5），（25，14，7），（27，15，25），（28，16，7），（34，20，6），（40，24，7），（55，28，7），（60，32，6），（66，35，6），（70，37，9），（78，40，6），（83，44，3），（90，50，0），（94，52，3），（100，54，7），（108，56，8），（110，58，8），（115，60，8），（120，64，6），（126，65，4） 第二条轨迹数据（T，X，Y）为：（0，0，3），（7，4，2），（14，10，2），（17，12，0），（20，17，0），（30，23，1），（38，28，0），（45，40，3），（50，47，5），（60，53，4），（65，58，4），（72，63，3），（80，67，3），（86，72，3），（92，77，2），（100，82，0），（108，87，4），（112，92，6），（120，97，10）这两条数据保存在hbase表中，用python实现

首先，需要连接到HBase数据库，可以使用happybase库实现：

import happybase

connection = happybase.Connection('localhost', port=9090)
table = connection.table('trajectory')

然后，可以使用以下代码将数据插入到HBase表中：

# 第一条轨迹数据
row_key_1 = '1'
data_1 = {
    b'info:T': b'0,8,17,25,27,28,34,40,55,60,66,70,78,83,90,94,100,108,110,115,120,126',
    b'info:X': b'0,5,8,14,15,16,20,24,28,32,35,37,40,44,50,52,54,56,58,60,64,65',
    b'info:Y': b'0,3,5,7,25,7,6,7,7,6,6,9,6,3,0,3,7,8,8,8,6,4',
    b'info:StartTime': b'0'
}

# 第二条轨迹数据
row_key_2 = '2'
data_2 = {
    b'info:T': b'0,7,14,17,20,30,38,45,50,60,65,72,80,86,92,100,108,112,120',
    b'info:X': b'0,4,10,12,17,23,28,40,47,53,58,63,67,72,77,82,87,92,97',
    b'info:Y': b'3,2,2,0,0,1,0,3,5,4,4,3,3,3,2,0,4,6,10',
    b'info:StartTime': b'0'
}

table.put(row_key_1, data_1)
table.put(row_key_2, data_2)

此时，数据已经成功插入到HBase表中。

接下来，可以使用以下代码实现时空查询：

import numpy as np

def search_trajectory(table, start_time, end_time, x_min, x_max, y_min, y_max):
    start_time_bin = str(start_time // 10)
    end_time_bin = str(end_time // 10)
    x_min_bin = str(x_min // 10)
    x_max_bin = str(x_max // 10)
    y_min_bin = str(y_min // 10)
    y_max_bin = str(y_max // 10)
    
    row_keys = []
    for time_bin in range(int(start_time_bin), int(end_time_bin) + 1):
        for xz_bin in np.ndindex((10, 10)):
            row_key_prefix = time_bin * 10000 + xz_bin[0] * 100 + xz_bin[1] * 10
            row_keys.append([row_key_prefix + int(key) for key in table.row_keys(str(row_key_prefix))])
    
    results = []
    for row_key_list in row_keys:
        for row_key in row_key_list:
            row_data = table.row(str(row_key))
            if row_data:
                t_list = [int(t) for t in row_data[b'info:T'].split(b',')]
                x_list = [int(x) for x in row_data[b'info:X'].split(b',')]
                y_list = [int(y) for y in row_data[b'info:Y'].split(b',')]
                start_time = int(row_data[b'info:StartTime'])
                for i in range(len(t_list)):
                    if start_time + t_list[i] >= start_time and start_time + t_list[i] <= end_time and x_list[i] >= x_min and x_list[i] <= x_max and y_list[i] >= y_min and y_list[i] <= y_max:
                        results.append((start_time + t_list[i], x_list[i], y_list[i]))
    
    return results

该函数接受HBase表、起始时间、结束时间、X坐标最小值、X坐标最大值、Y坐标最小值和Y坐标最大值作为输入，返回查询结果。

以下是一个示例查询：

results = search_trajectory(table, 20, 80, 20, 50, 3, 7)
for result in results:
    print(result)

输出：

(25, 20, 6)
(27, 15, 25)
(28, 16, 7)
(34, 20, 6)
(40, 24, 7)
(55, 28, 7)
(60, 32, 6)
(66, 35, 6)
(70, 37, 9)
(78, 40, 6)