Breitengrad Längengrad-Koordinaten, Status-Code in R

Gibt es eine schnelle Möglichkeit zum umwandeln von Längen-und Breitengrad-Koordinaten, Status-codes in R? Ich habe mit der PLZ-Paket als look-up-Tabelle, aber es ist zu langsam, wenn ich die Abfrage viel von lat/long-Werte

Wenn nicht in R gibt es eine Möglichkeit, diese mithilfe der google geocoder oder jede andere Art von schnellen Abfragen service?

Dank!

InformationsquelleAutor der Frage Michael | 2012-01-05

  1. 35

    Hier ist eine Funktion, die eine Daten.Rahmen des lat-longs in den unteren 48 Staaten, und für jeden Punkt, gibt den Status zurück, in dem es sich befindet.

    Meisten die Funktion einfach bereitet die SpatialPoints und SpatialPolygons Objekte, die erforderlich, durch die over() Funktion in der sp Paket, die nicht die eigentliche schwerarbeit der Berechnung der 'Schnittmenge' der Punkte und Polygone:

    library(sp)
library(maps)
library(maptools)

# The single argument to this function, pointsDF, is a data.frame in which:
#   - column 1 contains the longitude in degrees (negative in the US)
#   - column 2 contains the latitude in degrees

latlong2state <- function(pointsDF) {
    # Prepare SpatialPolygons object with one SpatialPolygon
    # per state (plus DC, minus HI & AK)
    states <- map('state', fill=TRUE, col="transparent", plot=FALSE)
    IDs <- sapply(strsplit(states$names, ":"), function(x) x[1])
    states_sp <- map2SpatialPolygons(states, IDs=IDs,
                     proj4string=CRS("+proj=longlat +datum=WGS84"))

    # Convert pointsDF to a SpatialPoints object 
    pointsSP <- SpatialPoints(pointsDF, 
                    proj4string=CRS("+proj=longlat +datum=WGS84"))

    # Use 'over' to get _indices_ of the Polygons object containing each point 
    indices <- over(pointsSP, states_sp)

    # Return the state names of the Polygons object containing each point
    stateNames <- sapply(states_sp@polygons, function(x) x@ID)
    stateNames[indices]
}

# Test the function using points in Wisconsin and Oregon.
testPoints <- data.frame(x = c(-90, -120), y = c(44, 44))

latlong2state(testPoints)
[1] "wisconsin" "oregon" # IT WORKS

    InformationsquelleAutor der Antwort Josh O'Brien

  2. 6

    Du kannst es in ein paar Zeilen von R.

    library(sp)
library(rgdal)
#lat and long
Lat <- 57.25
Lon <- -9.41
#make a data frame
coords <- as.data.frame(cbind(Lon,Lat))
#and into Spatial
points <- SpatialPoints(coords)
#SpatialPolygonDataFrame - I'm using a shapefile of UK counties
counties <- readOGR(".", "uk_counties")
#assume same proj as shapefile!
proj4string(points) <- proj4string(counties)
#get county polygon point is in
result <- as.character(over(points, counties)$County_Name)

    InformationsquelleAutor der Antwort tingenek

  3. 2

    Sehen ?über die im sp-Paket. Sie müssen über die staatlichen Grenzen als SpatialPolygonDataFrame.

    InformationsquelleAutor der Antwort mdsumner

Schreibe einen Kommentar