Hi,
Hereโs my latest challenge:
Iโm looking to get the 5km tiles out of this file: os_bng_grids.7z.
I extracted the content, a 200MB gpgk file of the same name as the 7z file, and can load this effortlessly using GeoIO.jl
but I cannot coax it to reveal its structure or content in a way that I can access - or at all infact.
geo = GeoIO.load("os_bng_grids.gpkg")
println(values(geo))
println((geo[geo.tile_name.=="HL", ["geometry"]]))
shows the tile_names
, but these are the names of the 100km squares. Iโm looking for the names and grid coords of the 5km squares.
(tile_name = ["HL", "HM", "HN", "HO", "HP", "HQ", "HR", "HS", "HT", "HU", "HV", "HW", "HX", "HY", "HZ", "JL", "JM", "JQ", "JR", "JV", "JW", "NA", "NB", "NC", "ND", "NE", "NF", "NG", "NH", "NJ", "NK", "NL", "NM", "NN", "NO", "NP", "NQ", "NR", "NS", "NT", "NU", "NV", "NW", "NX", "NY", "NZ", "OA", "OB", "OF", "OG", "OL", "OM", "OQ", "OR", "OV", "OW", "SA", "SB", "SC", "SD", "SE", "SF", "SG", "SH", "SJ", "SK", "SL", "SM", "SN", "SO", "SP", "SQ", "SR", "SS", "ST", "SU", "SV", "SW", "SX", "SY", "SZ", "TA", "TB", "TF", "TG", "TL", "TM", "TQ", "TR", "TV", "TW"],)
1ร1 GeoTable over 1 view(::GeometrySet, [1])
This must be simple, but itโs currently defeating me!
Can anyone point me in the right direction?
Thanks!