Show how to convert a long-lat point csv to sf

[jsta]

This issue is a follow-on from #95 where it was discussed that we never really show how to create an sf object by setting the projection with an EPSG code. My feeling is that we should do this in Ep 10 (Convert from .csv to a Shapefile in R) by showing sf object creation from a non-projected long-lat point csv:

sf_object <- sf::st_as_sf(data_frame, coords = c(longname, latname), crs = 4326)

I see two possibilities here:

• Convert data/NEON-DS-Site-Layout-Files/HARV/HARV_PlotLocations.csv to long-lat

• Bring in a different dataset. Perhaps this is an opportunity to bring-in some non-NEON data?

[ali-rawail]

This code is not working. i want to change this csv into raster

foresta<-carbonstock

head(carbonstock)

data<-data.frame(carbonstock$Longitude,carbonstock$Latitude,carbonstock$Carbon)

points from scratch

coords = cbind(carbonstock$Latitude, carbonstock$Longitude)

sp = SpatialPoints(coords)

make spatial data frame

spdf = SpatialPointsDataFrame(coords, data)

spdf = SpatialPointsDataFrame(sp, data)

promote data frame to spatial

coordinates(data) = cbind(carbonstock$Latitude, carbonstock$Longitude)

coordinates(data) = ~lon + lat

back to data

as.data.frame(data)

plot(data,)

library(raster)

dfr <- rasterFromXYZ(data) #Convert first two columns as lon-lat and third as value

plot(dfr)

dfr

library(raster)

create spatial points data frame

spg <- data

x<-carbonstock$Latitude

y<-carbonstock$Longitude

coordinates(spg) <- ~ x + y

coerce to SpatialPixelsDataFrame

gridded(spg) <- TRUE

coerce to raster

rasterDF <- raster(spg)