utrecht_district: Administrative areas of the Netherlands aroung the city of...

In goepp/graphseg: Segmentation of Graph-Based Signals

Administrative areas of the Netherlands aroung the city of Utrecht

Description

Geographical data of the Dutch administrative units "District".

Usage

utrecht_district

graph_utrecht_district

Format

utrecht_district

An sf object, with lattitude and longitude. With 7 variables and 650 zones.

graph_utrecht_district

The adjacency graph as an igraph object.

An object of class sf (inherits from data.frame) with 650 rows and 7 columns.

An object of class igraph of length 650.

Source

The data set comes from https://geodata.nationaalgeoregister.nl/.

Examples

data(utrecht_district); data(graph_utrecht_district)
coord <- sf::st_coordinates(sf::st_centroid(utrecht_district))
adj_municip <- as(as(igraph::as_adjacency_matrix(graph_utrecht_district, type = "both"),
                     "symmetricMatrix"),
                  "TsparseMatrix")
edge_list <- data.frame(adj_municip@i + 1, adj_municip@j + 1)
segment_df <- cbind(coord[edge_list[, 1], ], coord[edge_list[, 2], ])
ptmat <- as.matrix(segment_df[, 1:4])[2:nrow(segment_df), ]
linesegs <- lapply(split(ptmat, 1:nrow(ptmat)), function(x) {
    x <- matrix(x, nrow = 2, byrow = TRUE)
    x <- sf::st_linestring(x)})
final_sf <- sf::st_sf(sf::st_sfc(linesegs), 'ID' = 1:length(sf::st_sfc(linesegs)))
op <- par(mar = rep(0, 4))
plot(sf::st_geometry(utrecht_district), lwd = 0.6, border = "grey")
plot(sf::st_geometry(final_sf), lwd = 0.5, add = TRUE)
plot(sf::st_centroid(utrecht_district), add = TRUE, col = "black", pch = 20,
     cex = 0.5)
par(op)

goepp/graphseg documentation built on Nov. 23, 2023, 9:29 a.m.