Как преобразовать фрейм данных в пространственные координаты

Ответ 1

Сначала вы берете столбцы lon и lat и создаете объект для coord. Затем вы вычтите их из исходного фрейма данных и создайте новый объект. Вы, наконец, используете SpatialPointsDataFrame() для создания SpatialPointsDataFrame. Когда вы создаете SpatialPointsDataFrame, вам нужно назначить proj4string. Выберите подходящий для вас.

В вашем случае у вас нет других столбцов, но lon и lat, метод не будет работать. Я намеренно оставил lon и lat @data.

DATA

mydf <- structure(list(longitude = c(128.6979, 153.0046, 104.3261, 124.9019, 
126.7328, 153.2439, 142.8673, 152.689), latitude = c(-7.4197, 
-4.7089, -6.7541, 4.7817, 2.1643, -5.65, 23.3882, -5.571)), .Names = c("longitude", 
"latitude"), class = "data.frame", row.names = c(NA, -8L))


### Get long and lat from your data.frame. Make sure that the order is in lon/lat.

xy <- mydf[,c(1,2)]

spdf <- SpatialPointsDataFrame(coords = xy, data = mydf,
                               proj4string = CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"))


#> str(spdf)
#Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots
#[email protected] data       :'data.frame':  8 obs. of  2 variables:
#.. ..$ longitude: num [1:8] 129 153 104 125 127 ...
#.. ..$ latitude : num [1:8] -7.42 -4.71 -6.75 4.78 2.16 ...
#[email protected] coords.nrs : num(0) 
#[email protected] coords     : num [1:8, 1:2] 129 153 104 125 127 ...
#.. ..- attr(*, "dimnames")=List of 2
#.. .. ..$ : NULL
#.. .. ..$ : chr [1:2] "longitude" "latitude"
#[email protected] bbox       : num [1:2, 1:2] 104.33 -7.42 153.24 23.39
#.. ..- attr(*, "dimnames")=List of 2
#.. .. ..$ : chr [1:2] "longitude" "latitude"
#.. .. ..$ : chr [1:2] "min" "max"
#[email protected] proj4string:Formal class 'CRS' [package "sp"] with 1 slot
#.. .. [email protected] projargs: chr "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"

Ответ 2

Или используя sf вместо объектов sp (подробнее о Простые функции для R или перейдите с sp на sf здесь):

require(sf)
my.df <- read.table(text="
                    longitude    latitude
                    128.6979    -7.4197
                    153.0046    -4.7089
                    104.3261    -6.7541
                    124.9019    4.7817
                    126.7328    2.1643
                    153.2439    -5.6500
                    142.8673    23.3882
                    152.6890    -5.5710",
                    header=TRUE)

my.sf.point <- st_as_sf(x = my.df, 
                        coords = c("longitude", "latitude"),
                        crs = "+proj=longlat +datum=WGS84")
# simple plot
plot(my.sf.point)
# interactive map:
require(mapview)
mapview(my.sf.point)

# convert to sp object if needed
my.sp.point <- as(my.sf.point, "Spatial")

Ответ 3

С

structure(list(longitude = c(128.6979, 153.0046, 104.3261, 124.9019, 
126.7328, 153.2439, 142.8673, 152.689), latitude = c(-7.4197, 
-4.7089, -6.7541, 4.7817, 2.1643, -5.65, 23.3882, -5.571)), .Names = c("longitude", "latitude"), class = "data.frame", row.names = c(NA, -8L))

Для преобразования в SpatialPointsDataFrame

coordinates(df) <- cbind(df$longitude , df$latitude)

Как указано @jazzurro, вам, вероятно, потребуется назначить CRS для вашего пространственного объекта.

proj4string(df) = CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0")

обратный процесс SpatialPointsDataFrame к исходному df

df <- data.frame(longitude = coordinates(df)[,1], latitude = coordinates(df)[,2])