From cdd3f4b5ffabe5e7e9177951b72b559e7ddf49be Mon Sep 17 00:00:00 2001 From: h9b Date: Sat, 31 Aug 2024 22:45:35 +0200 Subject: [PATCH] add bboxes --- book/02_corporate_edits.ipynb | 16 +++++++++++----- 1 file changed, 11 insertions(+), 5 deletions(-) diff --git a/book/02_corporate_edits.ipynb b/book/02_corporate_edits.ipynb index d7e7e4c..12d690a 100644 --- a/book/02_corporate_edits.ipynb +++ b/book/02_corporate_edits.ipynb @@ -148,6 +148,8 @@ " 'colombia': (-78.9909352282, -4.29818694419, -66.8763258531, 12.4373031682), \n", " 'indonesia': (95.2930261576, -10.3599874813, 141.03385176, 5.47982086834),\n", " 'united_arab_emirates': (51.5795186705, 22.4969475367, 56.3968473651, 26.055464179)\n", + " #'south_america': (-93.691406,-58.263287,-22.675781,14.859850), # only do this if you have a lot of RAM!\n", + " #'africa': (-24.609375,-39.095963,56.074219,36.173357) # only do this if you have a lot of RAM!\n", "}\n", "\n", "selected_region = 'colombia'\n", @@ -180,7 +182,11 @@ " \"mapbox\": 22,\n", " \"Komoot\": 23,\n", " \"AppLogica\": 24\n", - "}" + "}\n", + "\n", + "\n", + "# Define output h3 cell resolution\n", + "h3_cell_resolution = 5" ] }, { @@ -283,8 +289,8 @@ "SELECT\n", " a.user_id,\n", " a.valid_from,\n", - " a.h3_r5,\n", - " h3_cell_to_boundary_wkt(a.h3_r5) as geometry,\n", + " h3_cell_to_parent(a.h3_r5, {h3_cell_resolution}) as h3_cell,\n", + " h3_cell_to_boundary_wkt(h3_cell) as geometry,\n", " CAST(changeset['hashtags'][1] AS VARCHAR) as hashtag_string,\n", " CASE\n", "\t\tWHEN hashtag_string ILIKE '%amap%' THEN 'amap'\n", @@ -348,12 +354,12 @@ " SELECT\n", " epoch_ms(date_trunc('month', valid_from)) as month,\n", " corporation,\n", - " h3_r5,\n", + " h3_cell,\n", " geometry,\n", " count(*) as n_edits,\n", " count(distinct user_id) as n_users\n", " FROM osm_data\n", - " GROUP by month, corporation, h3_r5, geometry;\n", + " GROUP by month, corporation, h3_cell, geometry;\n", "\"\"\"\n", "\n", "df = con.sql(map_query).df()\n",