# Load required libraries library(sp) library(raster) # Read raster image raster_img <- raster("path_to_raster_image.tif") # Read point data point_data <- read.csv("path_to_point_data.csv") # Create SpatialPoints object coordinates(point_data) <- c("Longitude", "Latitude") # Convert to SpatialPointsDataFrame object sp_df <- SpatialPointsDataFrame(point_data, data = point_data) # Extract values corresponding to raster image values <- extract(raster_img, sp_df) # Add extracted values to point data point_data$Raster_Values <- values # Output to table write.csv(point_data, "output_file.csv", row.names = FALSE)