𝚖𝚊𝚌𝚔𝚊 𝚜𝚣𝚎𝚌𝚑𝚗𝚘
LightNews LightNews
banner
mackaszechno.bsky.social
𝚖𝚊𝚌𝚔𝚊 𝚜𝚣𝚎𝚌𝚑𝚗𝚘
@mackaszechno.bsky.social
1.3K followers 630 following 3.4K posts
#rstats geospatial data analyst. GIS expert, of a sort. Sci-fi fan and would be writer. Black cat owner. Occasionally apocalyptical. Can't DM.
WORDS ARE QUICK AND WORDS ARE VAIN, THE SINGLE SURE AND FINAL ANSWER MUST BE PAIN. -Dan Simmons, Olympos
Posts Replies Media Videos
mackaszechno.bsky.social
Raster map of West Sussex in England from Ordnance Survey Open Data. 'tis been a long month and I didn't feel like fighting {terra} on this one. Raster map overlaid with vector data representing the border of West Sussex and the 3 relevant OS grid squares.
#30DayMapChallenge | Day 29 - Raster #QGIS
November 28, 2025 at 10:00 PM
mackaszechno.bsky.social
Wait... so it's a tree not one of these?

Puts a whole new dimension to your hunt. You know trees can't run too fast? Seriously, give them a few centuries head start if your gonna hunt them. At the very least, use a blindfold 🤣
https://en.wikipedia.org/wiki/House_sparrow
November 28, 2025 at 10:50 AM
mackaszechno.bsky.social
UK at night. Horsham and Crawley are easily identifiable.
upload.wikimedia.org/wikipedia/co...
UK at night.
November 28, 2025 at 8:56 AM
mackaszechno.bsky.social
Firefly map of West Sussex in England, showing areas of light pollution
#30DayMapChallenge | Day 28 - Dark #rstats
(code in reply)
Map of West Sussex showing streetlights on public land.
November 28, 2025 at 8:45 AM
mackaszechno.bsky.social
The Rapes Of Sussex
A rape is a traditional territorial sub-division of the county of Sussex in England, formerly used for various administrative purposes.
#30DayMapChallenge | Day 27 - Boundaries #rstats
(code in reply)
Map of the six rapes of Sussex and their Hundreds
November 27, 2025 at 3:40 PM
mackaszechno.bsky.social
Cross-border mobility between Afghanistan and its neighbours aggregated for October 2025. Much more needed, like proper borders and coloured arrows.
#30DayMapChallenge | Day 26 - Transport #rstats
(code in reply)
map of Afghanistan.. work needed
November 26, 2025 at 4:36 PM
mackaszechno.bsky.social
Proximity of addresses around target location anonymised by count inside hexbin. Rings show 0-0.5, 0.5-0.75, 0.75-1, 1-2, and 2-3 miles from target.
#30DayMapChallenge | Day 25 - Hexagons #rstats
November 25, 2025 at 3:00 PM
mackaszechno.bsky.social
Place names in West Sussex in English, Korean, and Ukrainian.
Plus a very lightly modified maplibregl js example of the same.
#30DayMapChallenge | Day 24 - Places and their names #rstats
(code in alt)
d <- read_csv("./data/ko.csv") |> st_as_sf(coords = c("lon", "lat"), crs = 4326) maplibre(style = carto_style("voyager-no-labels")) |> fit_bounds(wsx_mp) |> add_fill_layer(id = "bob", source = wsx_mp, fill_opacity = 0.1, fill_color = "navy") |> add_circle_layer(id = "locations", source = d, circle_color = "navy", circle_stroke_color = "white", circle_stroke_width = 1 ) |> add_symbol_layer(id = "names", source = d, text_field = concat( get_column("english"), "\n", get_column("korean"), "\n", get_column("ukrainian")), text_offset = c(0, 2.5), text_allow_overlap = FALSE, text_size = 10 ) |> add_control("<div style='padding:5px'> <b>Place names in West Sussex</b><br /> in English, Korean, and Ukrainian<br /> macka szechno<br />24 November 2025 </div>") <script> const map = new maplibregl.Map({ container: 'map', style: 'https://tiles.openfreemap.org/styles/bright', center: [16.05, 48], zoom: 2.9 }); document .getElementById('buttons') .addEventListener('click', (event) => { const language = event.target.id.substr('button-'.length); // Use setLayoutProperty to set the value of a layout property in a style layer. // The three arguments are the id of the layer, the name of the layout property, // and the new property value. map.setLayoutProperty('label_country_1', 'text-field', [ 'get', `name:${language}` ]); map.setLayoutProperty('label_country_2', 'text-field', [ 'get', `name:${language}` ]); map.setLayoutProperty('label_country_3', 'text-field', [ 'get', `name:${language}` ]); }); </script>
November 25, 2025 at 2:53 PM
mackaszechno.bsky.social
Another pic of Alice. She likes to wait at the door of my home office so she can follow me downstairs and demand treats.
Black cat patiently loafing #blackcat
November 22, 2025 at 10:00 PM
mackaszechno.bsky.social
My work setup. My process is constantly referring back to saved snippets of code that worked, reviewing help in RStudio and online, and deleting off all the non-working stuff until I get something I think is acceptable 😜
#30DayMapChallenge | Day 23 - Process
My setup with black cat and hopefully nothing to incriminating.
November 22, 2025 at 8:10 PM
mackaszechno.bsky.social
Ice shelf collapses in Antarctica using {mapgl}.
#30DayMapChallenge | Day 22 - Data challenge: Natural Earth #rstats
(code in reply)
Antarctica in {mapgl} with "WGS 84 / Antarctic Polar Stereographic EPSG:3031" shown on a globe as a circle. "WGS 84 / Antarctic Polar Stereographic EPSG:3031" viewed projected in {mapview} rather than as a globe.
November 22, 2025 at 7:47 PM
mackaszechno.bsky.social
Crawley and Gatwick...
Zoomed in to reveal more bus stops and Gatwick Airport train station and its 16.5m annual passengers.
November 21, 2025 at 3:09 PM
mackaszechno.bsky.social
Icons, revisiting Day 7.
#30DayMapChallenge | Day 21 - Icons #rstats
(code in reply)
November 21, 2025 at 3:07 PM
mackaszechno.bsky.social
Natural Earth to the rescue with rivers_lake_centerline, lakes, ocean, and bathy from 0-10000m.
#30DayMapChallenge | Day 20 - Water #rstats
(code in alt)
bathy_all_simp <- bathy_all |> ms_simplify(keep = 0.05) maplibre(style = carto_style("dark-matter-no-labels")) |> add_fill_layer(id = "ocean", source = ocean50, fill_color = "#008b8b") |> add_fill_layer(id = "lakes", source = lakes50, fill_color = "#e0eeee") |> add_line_layer(id = "rivers", source = rivers50, line_color = mapgl::interpolate( column = "scalerank", type = "linear", values = c(1, 12), stops = c("navy", "lightblue") ), line_width = mapgl::interpolate( column = "scalerank", type = "linear", values = c(1, 12), stops = c(2, 1) ), popup = "name") |> add_fill_layer(id = "bathy", source = bathy_all_simp, fill_color = "#f0fff0", fill_opacity = 0.4, fill_outline = mapgl::interpolate( column = "depth", values = c(0, 10000), stops = c("lightblue", "darkblue") ) ) |> add_control("<div style='padding:5px'> <b>Oceans, Rivers, and Lakes of the world</b><br /> with bathymetry from 0 to 10000m<br /> data: <a href='https://www.naturalearthdata.com/'>naturalearthdata.com</a><br /> macka szechno<br />20 November 2025 </div>") North Atlantic zoomed in
November 20, 2025 at 6:10 PM
mackaszechno.bsky.social
A day late due to man flu. West Sussex in one correct project (BNG 27700) and 3 inappropriate projections with area calculations showing how the deformation affects the total.
Note: I don't use {ggplot} much so plots are fairly basic.
#30DayMapChallenge | Day 19 - Projections #rstats
(code in alt)
library(sf) library(units) library(dplyr) library(ggplot2) library(patchwork) d <- read_sf("./OS.WESTSUSSEX_DISTRICT.shp") d <- d |> select(NAME) plot_projection <- function(d, crs, text){ dt <- d |> st_transform(crs = crs) ggplot(dt) + geom_sf(color = "white", aes(fill = NAME)) + coord_sf(crs = crs) + geom_sf_text(aes(label = NAME), size = 3) + theme(legend.position = "none") + labs(title = paste0("West Sussex - ", text, "\nArea calculation: ", st_area(dt) |> set_units("km^2") |> sum() |> round(), "km²")) } p1 <- plot_projection(d, 27700, "EPSG:27700 British National Grid") p2 <- plot_projection(d, 4326, "EPSG:4326 Latlon") p3 <- plot_projection(d, 54030, "ESRI:54030 Robinson") p4 <- plot_projection(d, 54099, "ESRI:54099 Spilhaus") p1 + p2 + p3 + p4
November 20, 2025 at 10:13 AM
mackaszechno.bsky.social
Left a message on the moon too as a bit of a running gag...
November 18, 2025 at 9:40 PM
mackaszechno.bsky.social
And a view of the moon in QGIS as a WMS layer courtesy of webmap.lroc.asu.edu. Instructions at lunaserv.im-ldi.com/getting_star....
November 18, 2025 at 9:18 PM
mackaszechno.bsky.social
Planning out interrelations between scenes, chapters, books, and associated notes.
#30DayMapChallenge | Day 18 - Out of this world #obsidian
full scale graph mid-level, I guess quite close
November 18, 2025 at 9:14 PM
mackaszechno.bsky.social
Okay, phoning this one in as got bored of all the necessary setup required to get Python spatial working, so kept to the absolute basics. Undersea cables in Python (image1) and with R (image2). Two lines of code in each.
#30DayMapChallenge | Day 17 - A new tool #python
(code in alt)
# https://py.geocompx.org/01-spatial-data import pandas as pd import shapely import geopandas as gpd gdf = gpd.read_file('./shapefile/cable.geojson') gdf.plot() library(sf) library(mapgl) cable <- read_sf("./shapefile/cable.geojson") maplibre_view(cable, column = "slug", legend = FALSE)
November 17, 2025 at 3:48 PM
mackaszechno.bsky.social
Zoomed out.
November 16, 2025 at 9:08 PM
mackaszechno.bsky.social
Reusing ONS Output Areas with pop density data from Day 8 to create cells colour-coded to pop density. Also could be considered a 10-minute map as used only maplibre_view.
#30DayMapChallenge | Day 16 - Cells #rstats
(code in alt)
library(sf) library(dplyr) library(mapgl) # read in ONS population densities by Output Areas d <- read_sf("./shapefile/Output_Areas_2021_WSx_pop.shp") # make border polygon border <- d |> summarise() # make grid of squares to cover border polygon grid <- st_make_grid(d, cellsize = 1000) |> st_as_sf(crs = 27700) grid <- grid |> mutate(id = row_number()) # get centroid for each grid square pts <- st_centroid(grid) # assign observation from ONS OAs to centroids ons <- d |> select(Obsrvtn) pts1 <- pts |> st_intersection(ons) pts1 <- pts1 |> rename(`Population density for West Sussex` = Obsrvtn) # produce map maplibre_view(pts1, column = "Population density for West Sussex", n = 10)
November 16, 2025 at 9:06 PM
mackaszechno.bsky.social
Now with megawatts as circle radius (still very basic).
November 15, 2025 at 10:11 PM
mackaszechno.bsky.social
Randomly found myself face-to-face with this culprit today following the wife's decision that we take the bus to Littlehampton during the unseasonably warm weather. Had to get a photo or two before leaving 🤣
Littlehampton Miniature Railway - end of the line (west-end) Littlehampton Miniature Railway sign
November 15, 2025 at 9:43 PM
mackaszechno.bsky.social
Okay in the end the simplest solution was to add in extra columns to the dataset, one for each decade ie S1920 for the 1920s, and then use set_filter to modify the original map.

Each section then becomes...
on_section("map", "1920", { maplibre_proxy("map") |> clear_markers() |> set_filter("all", filter = list("==", "S1920", TRUE)) |> fit_bounds(d) })
November 15, 2025 at 9:06 PM
mackaszechno.bsky.social
Maplibre in storymap mode showing changes in fuel-types for power generation in the UK.

Kudos to @kylewalker.bsky.social for the {mapgl} package and examples 😁

Data: www.gov.uk/government/c...

#30DayMapChallenge | Day 15 - Fire #rstats
November 15, 2025 at 8:57 PM