Introduction to pkmapr

Installation

Install pkmapr from CRAN:

install.packages("pkmapr")

Or install the development version from GitHub:

remotes::install_github("abdullahumer1101/pkmapr")

Your first map

Retrieve province boundaries and produce a map in two lines:

provinces <- get_provinces()
pk_map(provinces)

Look up names before joining

Official administrative names in the OCHA/HDX data may differ from common spellings. Use pk_dictionary() to confirm names and codes before filtering or joining:

# All provinces with their codes
pk_dictionary("provinces")

# Districts in Punjab
pk_dictionary("districts", province = "Punjab")

# Tehsils in Lahore district
pk_dictionary("tehsils", district = "Lahore")

Join your own data

pk_join() merges a data frame into an sf object by a shared code column, keeping geometries intact:

library(dplyr)

my_data <- data.frame(
  district_code = c("PK603", "PK604"),
  value         = c(42, 37)
)

districts <- get_districts() |>
  pk_join(my_data, by = "district_code")

pk_map(districts, fill = "value", title = "My Values")

Interactive maps

pk_map_interactive() produces a leaflet map with popups:

pk_map_interactive(
  districts,
  fill  = "value",
  popup = c("district_name", "value")
)

Next steps

• vignette("spatial-analysis-pkmapr") — buffers, centroids, and point-in-polygon operations
• vignette("epidemiology-pkmapr") — spatial autocorrelation, LISA clusters, and hotspot detection