Update url, close #9

Author: Robinlovelace

R/get_boundaries.R

@@ -19,7 +19,7 @@ get_boundaries = function(type = "la") {
     u = "https://opendata.arcgis.com/datasets/cc5941be78a8458393a03c69518b2bf9_0.geojson" # April 2020 generalised 20m (BGC)
   } else if(type == "la") {
     # u = "https://opendata.arcgis.com/datasets/3b374840ce1b4160b85b8146b610cd0c_0.geojson" # May 2020 generalised 20m (BGC)
-    u = "https://opendata.arcgis.com/datasets/0c09b7cde8b44c4ab6e2a1e47a91e400_0.geojson" # December 2011 EW generalised 20m (BGC)
+    u = "https://github.com/ITSLeeds/pct/releases/download/v0.9.7/LSOA_.Dec_2011._Boundaries_Super_Generalised_Clipped_.BSC._EW_V3.geojson" # December 2011 EW generalised 20m (BGC)
   }
   res = sf::read_sf(u)
   res

README.Rmd

@@ -74,6 +74,24 @@ jts04_gps$table_title
 jts04_gps$sheet
 ```
 
+As an example, our packages allow easy retrieval and mapping of data on the average journey time to employment centres (with 100 to 499 jobs) by public transport simply by running the following lines of code:
+
+```{r}
+jts_df = (type_code="jts05", spec="employment", sheet=2019, geo=True)
+jts.choropleth_map(
+    jts_df,
+    "500EmpPTt",
+    title="Travel time to jobs by public transport",
+    cblabel="Average travel time in minutes",
+    logscale=True,
+    cbticks=[1, 10, 100],
+    cbticklabels=["1", "10", "100"],
+    linewidth=0.001,
+)
+
+```
+
+
 Imagine you're interested in how average travel time to GP services changed between 2017 and 2019.
 You can do that as follows: