Dear all,
I have downloaded the NUTS 2 level data from library(“rgdal”) library(“RColorBrewer”) library(“classInt”) #library(“SmarterPoland”) library(fields) # Download Administrative Level data from EuroStat temp <- tempfile(fileext = ".zip") download.file(" http://ec.europa.eu/eurostat/cache/GISCO/geodatafiles/NUTS_2010_60M_SH.zip ", temp) unzip(temp) # Read data EU_NUTS <- readOGR(dsn = "./NUTS_2010_60M_SH/data", layer = "NUTS_RG_60M_2010") # Subset NUTS 2 level data map_nuts2 <- subset(EU_NUTS, STAT_LEVL_ == 2) I also have data for a variable by coordinates, which looks like this: structure(list(LON = c(-125.25, -124.75, -124.25, -124.25, -124.25, -124.25), LAT = c(49.75, 49.25, 42.75, 43.25, 48.75, 49.25), yr = c(2.91457704560515, 9.94774197180345, -2.71956412885765, -0.466213169185147, -36.6645659563374, 10.5168056769535)), .Names = c("LON", "LAT", "yr"), row.names = c(NA, 6L), class = "data.frame") I would like to match the coordinates to their corresponding NUTS 2 region - is this possible? Any help will be high appreciated. Thank you! Sincerely, Milu [[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo