Question

我有一个具有以下数据格式的SQLite数据库

...
2014-02-17T11：06：22.000-05：00，Vehicle3,40.820890，-73.935900
2014-02-17T11：06：23.000-05：00，Vehicle1,40.803433，-73.945087
2014-02-17T11：06：17.000-05：00，Vehicle2,40.798135，-73.946201
2014-02-17T11：10：10.000-05：00，Vehicle3,40.820890，-73.935900
2014-02-17T11：10：07.000-05：00，Vehicle1,40.802197，-73.945343
2014-02-17T11：09：59.000-05：00，Vehicle2,40.804895，-73.941317
2014-02-17T11：13：27.000-05：00，Vehicle3,40.820890，-73.935900
2014-02-17T11：13：17.000-05：00，Vehicle1,40.794255，-73.951131
2014-02-17T11：13：09.000-05：00，Vehicle2,40.810051，-73.937497
2014-02-17T11：15：37.000-05：00，Vehicle3,40.820890，-73.935900
2014-02-17T11：15：26.000-05：00，Vehicle1,40.789557，-73.954558
2014-02-17T11：15：49.000-05：00，Vehicle2,40.813135，-73.937353
2014-02-17T11：18：49.000-05：00，Vehicle3,40.820890，-73.935900
2014-02-17T11：19：08.000-05：00，Vehicle1,40.782017，-73.960065
2014-02-17T11：19：00.000-05：00，Vehicle2,40.817062，-73.938585
2014-02-17T11：22：37.000-05：00，Vehicle3,40.820890，-73.935900
2014-02-17T11：22：20.000-05：00，Vehicle1,40.778014，-73.962983
2014-02-17T11：22：44.000-05：00，Vehicle2,40.822828，-73.937887
2014-02-17T11：25：50.000-05：00，Vehicle3,40.820890，-73.935900
2014-02-17T11：26：03.000-05：00，Vehicle1,40.774126，-73.965815
2014-02-17T11：28：33.000-05：00，Vehicle3,40.820890，-73.935900
2014-02-17T11：28：09.000-05：00，Vehicle1,40.770644，-73.968356
...

第一列是日期/时间，第二是车辆ID，第三和第四是纬度和经度。

车辆的数量不是一成不变的，而是在一整天都在变化。日期/时间取决于每辆车的实际录制时间。该数据库包含超过一百万条记录，每3分钟采样率一次。

我的基本思想是提取车辆的运行顺序（逐个车辆），对日期/时间进行排序，计算时间间隔和时间间隔之间的距离（纬度和经度）作为距离，距离和时间间隔我能够计算速度。

问题是我不知道如何将方法结构化为SQLite select语句，我感谢任何帮助。

非常感谢！

Answer 1

我已经使用awk，我认为可以在Android上使用它 - 它可以很容易地转换为Perl或C代码。

它使用Haversine公式来计算距离。

它假设您的sqlite转储位于名为locations的文件中。

#!/bin/bash
awk -F, '
   {
      # Convert date to epoch seconds for added sanity
      tstr=$1;
      cmd="gnudate --date=" tstr " +%s"
      cmd | getline epoch
      close(cmd)

      # DEBUG print epoch,$2,$3,$4

      # Pick up all fields from current record
      vehicle=$2;lat=$3;lon=$4;

      # If we have a previous record for this vehicle we are in business
      if(lats[vehicle]){
         tdiff=epoch-epochs[vehicle]
         d=haversine(lat,lon,lats[vehicle],lons[vehicle])
         speed=3600*d/tdiff
         if(speed==0)speed="0 (stationary)"
         print $1,vehicle,speed
      }

      # Update last seen lats, lons, epoch for this vehicle for next iteration
      epochs[vehicle]=epoch
      lats[vehicle]=lat
      lons[vehicle]=lon
   }

   function haversine(lat1,lon1,lat2,lon2,  a,c,dlat,dlon) {
      dlat = radians(lat2-lat1)
      dlon = radians(lon2-lon1)
      lat1 = radians(lat1)
      lat2 = radians(lat2)
      a = (sin(dlat/2))^2 + cos(lat1) * cos(lat2) * (sin(dlon/2))^2
      c = 2 * atan2(sqrt(a),sqrt(1-a))
      # 6372 = Earth radius in km, for distance in km
      return 6372 * c
   }

   function radians(degree) { # degrees to radians
      return degree * (3.1415926 / 180.)

   }' locations

<强>输出：

2014-02-17T11:10:10.000-05:00   Vehicle3 0 (stationary)
2014-02-17T11:10:07.000-05:00   Vehicle1 2.23614
2014-02-17T11:09:59.000-05:00   Vehicle2 13.8954
2014-02-17T11:13:27.000-05:00   Vehicle3 0 (stationary)
2014-02-17T11:13:17.000-05:00   Vehicle1 19.1132
2014-02-17T11:13:09.000-05:00   Vehicle2 12.4564
2014-02-17T11:15:37.000-05:00   Vehicle3 0 (stationary)
2014-02-17T11:15:26.000-05:00   Vehicle1 16.6565
2014-02-17T11:15:49.000-05:00   Vehicle2 7.72184
2014-02-17T11:18:49.000-05:00   Vehicle3 0 (stationary)
2014-02-17T11:19:08.000-05:00   Vehicle1 15.5387
2014-02-17T11:19:00.000-05:00   Vehicle2 8.46043
2014-02-17T11:22:37.000-05:00   Vehicle3 0 (stationary)
2014-02-17T11:22:20.000-05:00   Vehicle1 9.53438
2014-02-17T11:22:44.000-05:00   Vehicle2 10.349
2014-02-17T11:25:50.000-05:00   Vehicle3 0 (stationary)
2014-02-17T11:26:03.000-05:00   Vehicle1 7.97182
2014-02-17T11:28:33.000-05:00   Vehicle3 0 (stationary)
2014-02-17T11:28:09.000-05:00   Vehicle1 12.6412

备注：

单位为km / h，将代码中的6372 km地球半径更改为3959英里（单位为英里/小时）。

您的date命令可能是date，而不是第6行的gnudate。

如果您要处理离线一段时间的车辆，请将计算tdiff的行移到if语句之上并测试if tdiff<60，这样您只需计算速度如果自上次职位以来的时间不到一分钟（比如说）。

Answer 2

SQLite没有将lat / lon坐标转换为公制坐标所需的数学函数（Android不允许你添加它们）。

根据存储在SQLite数据库中的位置（纬度和经度）和时间计算速度

2 个答案: