image(powehi,
col = hcl.colors(256, palette = "inferno", alpha = NULL, rev = FALSE,
fixup = TRUE), axes = FALSE, asp = 1)
title(main = "Messier 87's Black Hole: Powehi")
# based on the image, we expect one especially prominent
# persistent feature in 1D
ph <- cubical(powehi)
plot.new()
plot.window(
xlim = c(0, max(ph$death)),
ylim = c(0, max(ph$death)),
asp = 1
)
axis(1L)
axis(2L)
abline(a = 0, b = 1)
points(ph[ph$dim == 1L, c("birth", "death")], pch = 17L, col = "orange")
Run the code above in your browser using DataLab