GuestRT

Radius-Kotelniki-3

Mar 16th, 2018
88
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
SQL 3.12 KB | None | 0 0
  1. CREATE TABLE sandbox.yurbasov_radius_kotelniki3 AS
  2. -- RADIUS = 0.1 km
  3. -- ,acos(sin((90.0 - v_gps_y) * pi() / 180.0) * sin((90.0 - nextgpsy) * pi() / 180.0)
  4. --       * cos(radians(v_gps_x) - radians(nextgpsx)) +
  5. --           cos((90.0 - v_gps_y) * pi() / 180.0) * cos((90.0 - nextgpsy) * pi() / 180.0)) * 6371000 as distance
  6. -- UNITS: meters
  7. -- Волг 55.690742_37.82169
  8. -- Волг 55.690548_37.821686
  9. -- Волг 55.69926_37.7997
  10. -- Ряз 55.713745_37.809536
  11. WITH base1
  12. AS (SELECT
  13.   from_unixtime(unix_timestamp(v_time_check, 'dd.MM.yyyy HH:mm:ss'), 'yyyy-MM-dd') AS DAY,
  14.   v_regno,
  15.   CASE WHEN codd_kafka.v_gps_y > 55 THEN codd_kafka.v_gps_y ELSE codd_kafka.v_gps_x END AS camera_x,
  16.   CASE WHEN codd_kafka.v_gps_y > 55 THEN codd_kafka.v_gps_x ELSE codd_kafka.v_gps_y END AS camera_y,
  17.   lead(v_time_check, 1) OVER (partition BY v_regno ORDER BY v_time_check) AS r_previous_overall,
  18.   v_time_check
  19. FROM codd_data.codd codd_kafka
  20. WHERE (from_unixtime(unix_timestamp(v_time_check, 'dd.MM.yyyy HH:mm:ss'),'yyyy-MM-dd') <= '2018-02-28'
  21. AND from_unixtime(unix_timestamp(v_time_check, 'dd.MM.yyyy HH:mm:ss'),'yyyy-MM-dd') >= '2018-01-01')
  22. ), tgt1 AS
  23. (
  24. SELECT 55.690742 AS v_gps_x, 37.82169 AS v_gps_y
  25. UNION ALL SELECT 55.690548 AS v_gps_x, 37.821686 AS v_gps_y
  26. UNION ALL SELECT 55.69926 AS v_gps_x, 37.7997 AS v_gps_y
  27. UNION ALL SELECT 55.713745 AS v_gps_x, 37.809536 AS v_gps_y
  28. ), pre1 AS
  29. (SELECT
  30.   v_regno,
  31.   camera_x, camera_y,
  32.   MAX(abs(round((unix_timestamp(v_time_check) - unix_timestamp(r_previous_overall))/60, 0))) AS t_diff_next,
  33.   weekofyear(v_time_check) AS w_no,
  34.   COUNT(*) AS cnt
  35. FROM base1
  36. JOIN tgt1
  37. ON acos(sin((90.0 - v_gps_y) * pi() / 180.0) * sin((90.0 - camera_y) * pi() / 180.0) * cos(radians(v_gps_x) - radians(camera_x))
  38.        + cos((90.0 - v_gps_y) * pi() / 180.0) * cos((90.0 - camera_y) * pi() / 180.0)) * 6371000 <= 100
  39. GROUP BY v_regno,
  40.          camera_x, camera_y,
  41.          weekofyear(v_time_check)),
  42. f1 AS (SELECT pre1.v_regno, pre1.cnt FROM pre1 WHERE  t_diff_next >= 0),
  43. xxxx1
  44. AS (SELECT DISTINCT
  45.   v_regno
  46. FROM f1
  47. WHERE cnt >= 6),
  48. par1
  49. AS (SELECT
  50.   phoneno,
  51.   carno
  52. FROM (SELECT
  53.   phoneno,
  54.   carno,
  55.   MAX(from_unixtime(unix_timestamp(parkingstart, 'dd.MM.yyyy HH:mm:ss'), 'yyyy-MM-dd')) AS max_date,
  56.   COUNT(*) AS qty,
  57.   ROW_NUMBER() OVER (PARTITION BY carno ORDER BY COUNT(*) DESC) AS rn
  58. FROM parking_data.parking
  59. WHERE phoneno != 'NULL'
  60. AND from_unixtime(unix_timestamp(parkingstart, 'dd.MM.yyyy HH:mm:ss'), 'yyyy-MM-dd') >= '2017-01-01'
  61. GROUP BY phoneno,
  62.          carno) xx
  63. WHERE rn = 1), extra1 AS (
  64. SELECT
  65.     xxxx1.v_regno AS v_regno
  66.     , sec.regno_hash AS regno_hash
  67.     , sec.regno AS regno
  68. FROM xxxx1
  69.  JOIN secret.regno AS sec ON xxxx1.v_regno = sec.regno_hash
  70.  GROUP BY
  71.  v_regno, regno_hash, regno
  72. ), res1 AS
  73. (SELECT
  74.   extra1.v_regno,
  75.   extra1.regno,
  76.   par1.phoneno
  77. FROM extra1
  78. JOIN par1
  79.   ON extra1.v_regno = par1.carno
  80. WHERE par1.carno IS NOT NULL
  81. GROUP BY
  82. extra1.v_regno, extra1.regno,
  83. par1.phoneno), res AS
  84. (SELECT res1.phoneno, res1.regno, res1.v_regno FROM res1)
  85. SELECT s.phone, res.regno, res.v_regno FROM res JOIN secret.phones s ON s.phone_hash = res.phoneno
Add Comment
Please, Sign In to add comment