How to Convert data frame to spatial coordinates

by

First, you take the columns of lon and lat and create an object for coord. Then, you subtract them from the original data frame and create a new object. You finally use SpatialPointsDataFrame() to create a SpatialPointsDataFrame. When you create a SpatialPointsDataFrame, you need to assign proj4string. Choose an appropriate one for you.

In your case, you do not have any other columns but lon and lat, the method won’t work. I purposely left lon and 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
#..@ 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 ...
#..@ coords.nrs : num(0) 
#..@ coords     : num [1:8, 1:2] 129 153 104 125 127 ...
#.. ..- attr(*, "dimnames")=List of 2
#.. .. ..$ : NULL
#.. .. ..$ : chr [1:2] "longitude" "latitude"
#..@ 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"
#..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slot
#.. .. ..@ projargs: chr "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"

More Related Contents:

  1. Buffer (geo)spatial points in R with gbuffer
  2. How to efficiently calculate distance between pair of coordinates using data.table :=
  3. Function to calculate geospatial distance between two points (lat,long) using R [duplicate]
  4. Converting geo coordinates from degree to decimal
  5. Check if point is in spatial object which consists of multiple polygons/holes
  6. R Plot Filled Longitude-Latitude Grid Cells on Map
  7. sp::over() for point in polygon analysis
  8. How to plot a heat map on a spatial map
  9. Add raster to ggmap base map: set alpha (transparency) and fill color to inset_raster() in ggplot2
  10. Why does X[Y] join of data.tables not allow a full outer join, or a left join?
  11. Subset data frame based on multiple conditions [duplicate]
  12. UTF-8 file output in R
  13. Extract every nth element of a vector
  14. In R markdown in RStudio, how can I prevent the source code from running off a pdf page?
  15. Create columns from factors and count [duplicate]
  16. Unable to change Python path in reticulate
  17. Joining aggregated values back to the original data frame [duplicate]
  18. geom_boxplot with precomputed values
  19. Recommendations for Windows text editor for R [closed]
  20. Rscript: There is no package called …?
  21. Control the height in fluidRow in R shiny
  22. Dataframe create new column based on other columns
  23. Multiple time series in one plot
  24. Find K nearest neighbors, starting from a distance matrix
  25. Changing the order of dodged bars in ggplot2 barplot
  26. Using R to read out excel-colorinfo
  27. Existing function for seeing if a row exists in a data frame?
  28. how to plot ellipse given a general equation in R?
  29. Convert comma separated string to integer in R
  30. Interpolating timeseries

Leave a Comment