Skip to content

Instantly share code, notes, and snippets.

@mdsumner

mdsumner/gc_to_mtom.md

Last active March 17, 2025 13:42
Show Gist options
  • Save mdsumner/32d9222176b9d29c768178d804108441 to your computer and use it in GitHub Desktop.
Save mdsumner/32d9222176b9d29c768178d804108441 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
gc_to_mtom.md

Using function here by @h-a-graham

https://gist.github.com/h-a-graham/53a7ce9a3095245017acf93f6c674754

x <- majortom_embeds(c(xmin = 145, ymin = -43.5, xmax = 147.5, ymax = -41))
plot(wk::as_wkt(x$utm_footprint))
zlim <- range(unlist(x$embedding))
for (i in seq_len(nrow(x))) {
  ximage::ximage(matrix(x$embedding[[i]], 32, byrow = TRUE), 
                 as.numeric(wk::wk_bbox(wk::as_wkt(x$utm_footprint[[i]])))[c(1, 3, 2, 4)], 
                 add = TRUE, zlim = zlim, col = hcl.colors(128))
}

m <- reproj::reproj_xy(do.call(cbind, maps::map(plot = F)[1:2]), x$utm_crs[1], source = "EPSG:4326")
lines(m)

image

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment