What’s changed?

Andy Saunders writes:

A couple of things come to mind. The planet on the right contains, on average, around 70% fewer animals than the planet on the left.

The atmosphere of the planet on the right also contains around 100 ppm, or 30%, more CO2 than the planet on the left.

It seems to me that in many ways, humankind as a whole has not protected our Blue Marble very well in the 50-odd years between these images. Perhaps the next generation can do better.

Directional markers in R/leaflet

So you have used the excellent exiftool to extract all of the GPS-related information from a directory of photos in JPG format and write to a CSV file:

exiftool '-*GPS*' -ext jpg -csv . > outfile.csv

You’ve used R/leaflet to plot coordinates (latitude and longitude) before, but what about that tag named GPSImgDirection? It would be nice to have some kind of marker which indicates the direction in which you were facing when the photo was taken.

Continue reading

Brief thoughts on: the new Australian Bureau of Meteorology website

For those outside of Australia, the BOM website is where many Australians access weather information. It’s very popular with millions of visits per day, especially to the rain radar, and “the BOM” is even part of everyday language and culture. So much so that an attempt to prevent use of the acronym generated widespread criticism.

The recent redesign of the site, the first in over a decade, has proven controversial. My thoughts, as someone who uses many websites and even designed a few back in the day:

Continue reading

Pattern recognition in Google Maps

Headline: Second farm in the Hawkesbury region confirmed to have bird flu as biosecurity zone widened.

Quite rightly, the news organisation chooses not to name the farm, nor reveal its location. However, they do include several aerial photographs including this one.

Somewhat distinctive? Indeed, a little scrolling around the Hawkesbury region in Google Maps, with a focus on areas of likely farmland, soon reveals a candidate location.

Mind what you post. It may easier to identify than you think.

Data discovery: seasonal speed

Just writing this one quickly as it’s been hanging around my browser tabs for weeks…

I wrote Taking steps (in XML) almost 7 years ago and once in a while, I still grab Apple Health data from my phone and play around with it in R for a few minutes. Sometimes, curve fitting to a cloud of points generates a surprise.

library(tidyverse)
library(xml2)
theme_set(theme_bw())

health_data <- read_xml("~/Documents/apple_health_export/export.xml")

ws <- xml_find_all(health_data, ".//Record[@type='HKQuantityTypeIdentifierWalkingSpeed']") %>% 
    map(xml_attrs) %>% 
    map_df(as.list)

ws %>% 
    mutate(Date = ymd_hms(creationDate), 
                  value = as.numeric(value)) %>% 
    ggplot(aes(Date, value)) + 
    geom_point(size = 1, alpha = 0.2, color = "grey70", fill = "grey70") + 
    geom_smooth() + 
    labs(y = "Walking speed (km/h)", 
    title = "Walking speed data", 
    subtitle = "Apple Health 2020 - 2023")

Result:

Huh. Looks seasonal. Looks faster in the (southern) winter. Has that been reported before? Sure has.

It didn’t impress everyone but I thought it was interesting.