library(terra)
## terra version 1.1.2
##
## Attaching package: 'terra'
## The following object is masked from 'package:knitr':
##
## spin
longitude <- c(-72.875691, -72.869782, -72.862467, -72.858593, -72.850833, -72.849032, -72.858058)
latitude <- c(2.212863, 2.219539, 2.22007, 2.227792, 2.230666, 2.216927, 2.2097)
lonlat <- cbind(longitude, latitude)
pts <- vect(lonlat)
class(pts)
[1] "SpatVector"
attr(,"package")
[1] "terra"
pts
class : SpatVector
geometry : points
dimensions : 7, 0 (geometries, attributes)
extent : -72.87569, -72.84903, 2.2097, 2.230666 (xmin, xmax, ymin, ymax)
coord. ref. :
geom(pts)
geom part x y hole
[1,] 1 1 -72.87569 2.212863 0
[2,] 2 1 -72.86978 2.219539 0
[3,] 3 1 -72.86247 2.220070 0
[4,] 4 1 -72.85859 2.227792 0
[5,] 5 1 -72.85083 2.230666 0
[6,] 6 1 -72.84903 2.216927 0
[7,] 7 1 -72.85806 2.209700 0
crdref <- "+proj=longlat +datum=WGS84"
pts <- vect(lonlat, crs=crdref)
pts
class : SpatVector
geometry : points
dimensions : 7, 0 (geometries, attributes)
extent : -72.87569, -72.84903, 2.2097, 2.230666 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs
crs(pts)
[1] "GEOGCRS[\"unknown\",\n DATUM[\"World Geodetic System 1984\",\n ELLIPSOID[\"WGS 84\",6378137,298.257223563,\n LENGTHUNIT[\"metre\",1]],\n ID[\"EPSG\",6326]],\n PRIMEM[\"Greenwich\",0,\n ANGLEUNIT[\"degree\",0.0174532925199433],\n ID[\"EPSG\",8901]],\n CS[ellipsoidal,2],\n AXIS[\"longitude\",east,\n ORDER[1],\n ANGLEUNIT[\"degree\",0.0174532925199433,\n ID[\"EPSG\",9122]]],\n AXIS[\"latitude\",north,\n ORDER[2],\n ANGLEUNIT[\"degree\",0.0174532925199433,\n ID[\"EPSG\",9122]]]]"
# Generate random precipitation values, same quantity as points
precipvalue <- runif(nrow(lonlat), min=100, max=1000)
df <- data.frame(ID=1:nrow(lonlat), precip=precipvalue)
#Combine Spatial vector with the Dataframe
ptv <- vect(lonlat, atts=df, crs=crdref)
ptv
class : SpatVector
geometry : points
dimensions : 7, 2 (geometries, attributes)
extent : -72.87569, -72.84903, 2.2097, 2.230666 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs
LS0tDQp0aXRsZTogIkNvb3JkcyBpbnRvIGRhdGFmcmFtZSINCm91dHB1dDogaHRtbF9ub3RlYm9vaw0KLS0tDQoNCg0KDQpgYGB7cn0NCmxpYnJhcnkodGVycmEpDQojIyB0ZXJyYSB2ZXJzaW9uIDEuMS4yDQojIw0KIyMgQXR0YWNoaW5nIHBhY2thZ2U6ICd0ZXJyYScNCiMjIFRoZSBmb2xsb3dpbmcgb2JqZWN0IGlzIG1hc2tlZCBmcm9tICdwYWNrYWdlOmtuaXRyJzoNCiMjDQojIyAgICAgc3Bpbg0KDQoNCmxvbmdpdHVkZSA8LSBjKC03Mi44NzU2OTEsIC03Mi44Njk3ODIsIC03Mi44NjI0NjcsIC03Mi44NTg1OTMsIC03Mi44NTA4MzMsIC03Mi44NDkwMzIsIC03Mi44NTgwNTgpDQpsYXRpdHVkZSA8LSBjKDIuMjEyODYzLCAyLjIxOTUzOSwgMi4yMjAwNywgMi4yMjc3OTIsIDIuMjMwNjY2LCAyLjIxNjkyNywgMi4yMDk3KQ0KbG9ubGF0IDwtIGNiaW5kKGxvbmdpdHVkZSwgbGF0aXR1ZGUpDQoNCg0KcHRzIDwtIHZlY3QobG9ubGF0KQ0KY2xhc3MocHRzKQ0KDQpwdHMNCg0KZ2VvbShwdHMpDQoNCmNyZHJlZiA8LSAiK3Byb2o9bG9uZ2xhdCArZGF0dW09V0dTODQiDQpwdHMgPC0gdmVjdChsb25sYXQsIGNycz1jcmRyZWYpDQpwdHMNCg0KY3JzKHB0cykNCg0KDQoNCg0KDQoNCg0KYGBgDQoNCg0KYGBge3J9DQoNCiMgR2VuZXJhdGUgcmFuZG9tIHByZWNpcGl0YXRpb24gdmFsdWVzLCBzYW1lIHF1YW50aXR5IGFzIHBvaW50cw0KcHJlY2lwdmFsdWUgPC0gcnVuaWYobnJvdyhsb25sYXQpLCBtaW49MTAwLCBtYXg9MTAwMCkNCmRmIDwtIGRhdGEuZnJhbWUoSUQ9MTpucm93KGxvbmxhdCksIHByZWNpcD1wcmVjaXB2YWx1ZSkNCg0KI0NvbWJpbmUgU3BhdGlhbCB2ZWN0b3Igd2l0aCB0aGUgRGF0YWZyYW1lDQpwdHYgPC0gdmVjdChsb25sYXQsIGF0dHM9ZGYsIGNycz1jcmRyZWYpDQoNCnB0dg0KDQoNCg0KYGBgDQoNCg0KYGBge3J9DQoNCg0KDQpgYGANCg0KDQpgYGB7cn0NCg0KYGBgDQoNCg0KYGBge3J9DQoNCmBgYA0KDQoNCmBgYHtyfQ0KYGBgDQoNCg==